Application of the dielectric-dependent screened exchange potential approach to organic photocell materials

This paper discusses the fundamental features of the dielectric-dependent screened exchange potential approach for organic molecules and photocell materials. The energy diﬀerence (gap) between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) is a key factor when designing organic photocell materials. The magnitude of this energy gap strongly depends on the ratio of the Hartree–Fock (HF) exchange term combined with hybrid density functional theory (DFT) functionals. In ordinary hybrid DFT functionals, the HF exchange term ratio is often determined empirically and a constant value is employed for all materials. Conversely, the dielectric-dependent exchange potential approach employs a system-dependent value for this parameter, which is proportional to the inverse of the dielectric constant. The organic materials examined in this paper take different dielectric constant values according to their material types and unit lengths; therefore, their HF exchange ratios are also different. This paper discusses the relationship between the energy gaps, the dielectric constants, and the HF exchange ratios in the dielectric-dependent screened exchange potential approach. For a series of acene-type compounds, it was found that the HOMO–LUMO gap decreased when their conjugation systems are extended. The dielectric-dependent screened exchange potential approach demonstrated that the values of the dielectric constants of the compounds increased; consequently, the HF exchange ratio decreased. Similar results were obtained for a series of thiophene-based donor molecules. The calculated values were compared with those obtained via experimental measurements. We found that although the theoretical calculations of the HOMO–LUMO gap usually overestimate experimental optical gaps, a better agreement between the experimental and calculated values can be obtained if we correct for the exciton binding energy.


I. Introduction
In recent years, organic electronic devices such as electroluminescent displays, transistors, and photocells have been extensively studied because of their useful properties such as light weight, thin film structures, flexibility, low manufacturing costs, and design.In order to more widely utilize organic materials in electronic devices, a deeper understanding of their electronic structure properties is required.For example, the highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) gap of the donor material is an important indicator when fabricating organic photocell devices because the energy of the photon absorbed by the device is closely related to this gap. 1,2The energy of the light absorbed is roughly equivalent to the value obtained by subtracting the binding energy of an electron-hole pair (exciton) from the HOMO-LUMO energy gap.The dielectric constants of organic materials are usually smaller than those of inorganic semiconductors; therefore, the electron and the hole generated by photon absorption are strongly bound in organic materials.3][4][5][6][7][8][9] The energy offset at the donoracceptor junction can give excess energy to the exciton and assist the dissociation process required to generate a free electron and a hole. 8,9Several experiments have demonstrated that the energy conversion efficiency of an organic photocell device is suppressed in small offset cases.Thus, the ability to determine energy difference properties, such as the HOMO-LUMO gap, is essential in the development of organic photocell devices.
Unfortunately, it is not a simple matter to precisely predict the HOMO-LUMO gap of materials using density functional theory (DFT).For example, it is well known that local density approximation used in DFT calculations underestimates the HOMO-LUMO (band) gap as a result of the lack of nonlocal exchange interaction. 10,11Hybrid DFT methods, which possess partial nonlocal Hartree-Fock (HF) exchange terms, can reduce the bandgap underestimation, but ordinary hybrid DFT functionals sometimes fail to accurately reproduce the energy difference. 12This may be because the ratio of the HF exchange term is empirically determined in most hybrid DFT functionals.For example, in the B3LYP functional, the HF exchange ratio is empirically determined to reproduce the experimental results of molecules and atoms. 13To overcome this deficiency, we have proposed the dielectric-dependent screened exchange potential approach, in which the HF exchange term ratio can change according to the target system. 14,15In the dielectric-dependent exchange potential, the HF exchange ratio is proportional to the inverse of the dielectric constant.In contrast, an ordinary DFT functional usually employs a constant value for all materials.The dielectric-dependent approach has also been applied by other research groups and has been shown to give good descriptions of the energy gap for several systems. 12,16-28Skone et al. recently reported that their dielectric-dependent rangeseparated hybrid DFT functional can describe the electronic structure of molecular systems well. 29Our dielectric-dependent exchange potentials have been mainly applied to inorganic materials, although we previously demonstrated simple polymer calculations together with the Gaussian-based range-separation method. 30,31This paper discusses the dielectric-dependent exchange potential approach to organic molecules.We focus specifically on several donor materials suitable for use in organic photocell devices, because their energy gap properties are especially important for solar cells.In the next section, we briefly explain the calculation methodology for the screened exchange potential approach; calculation results are presented and discussed in Section III.In the final section, we give a summary of this paper.

II. Calculation method
In this section, we discuss several dielectric-dependent screened exchange potentials with Coulomb hole (COH) and screened exchange (SEX) interactions. 32The dielectric-dependent screened exchange potential can be derived from the dielectric function e(q) = 1 + [(e s À 1) À1 + b(q 2 /q TF 2 )] À1 . 14Here, q is the wave vector, q TF is the Thomas-Fermi wave number, and e s is the electronic part of the static dielectric constant.Bechstedt et al. reported that the value of b is independent of the target semiconductor, and in this paper we employ their suggested value of b = 1.563. 33The physical concept of the dielectricdependent exchange approach is summarized in the following equations: 34 (1.1) Here, r(r) is the electron density.The Fermi momentum k F (r) is written as k F ðrÞ ¼ 3p 2 rðrÞ À Á 1 3 .r 12 is the distance between two electrons.This paper employs the simple potential description.The precise one is found in the literature. 34,35Eqn (1.3) is used to correct the behavior of the dielectric function in the large q region. 33,34,36The first term gives the local COH interaction, and the second and third terms express the nonlocal SEX contributions.The third term may be the most important, because it represents the partial fragment of the HF exchange interaction in the potential.In the dielectric-dependent screened exchange potential, the HF exchange ratio is proportional to the inverse of the dielectric constant.In semiconductors and insulators, the screening effect is not complete, and therefore the bare interaction of 1/r appears in the third term of eqn (1.1), but the interaction is modified by the inverse of the dielectric constant. 14In the limit of e s = 1, the potential of eqn (1.1) reduces to the ordinary HF exchange interaction.In a metallic system in which e s -N, the SEX part of the dielectric-dependent exchange potential reduces to a Thomas-Fermi-like screening interaction. 36On the other hand, the nonlocal behavior of the second term, which is described by the Yukawa-type potential, is not strong because of the large q ˜TF values for molecules.Thus, we can approximate the weak nonlocal term by the local potential based on the dielectric-dependent Slater term as follows: 34,37 (2.1) (2.2) 3) The first and second terms of eqn (2.1) are local, but the third term is nonlocal.In this paper, the Yukawa-type potential, exp(ÀZr)/r, is replaced by the complementary error function erfc(2Zr/3)/r because of their similar behaviors. 14,38This paper also examines the following dielectric-dependent DFT potential to compare with eqn (1.1) and (2.1).
Here, V VWN is the Vosko-Wilk-Nusair correlation DFT potential. 39he dielectric-dependent exchange potential of eqn (3) may be more easily used in ordinary quantum chemistry, because the terms V Slater and V VWN are already prepared in many quantum chemistry program packages.Finally, we describe the self-consistent procedure used to automatically determine the dielectric constant. 21,40In our approach, the dielectric constant is determined by the following equation. 30,34,36,41 Here, the indices i and a represent the occupied and virtual (unoccupied) orbitals, respectively; E i and E a are the orbital energies for occupied and virtual orbitals of f i and f a , respectively; V is the volume of the target molecule; and an electron density of 0.01 a.u. is used to determine the boundary of the volume.The integral value of hf i (r)|exp(iqÁr)|f a (r)i is calculated as described in the previous study. 35We estimate the dielectric constant from the isotropically averaged value at the limit q -0.
In the self-consistent dielectric constant approach, the e s value used in eqn (1.1), (2.1), and ( 3) is renewed for each self-consistent field (SCF) step.Thus, information about the virtual orbital (conduction band) is taken into account through eqn ( 4) to obtain the electronic structure of the materials.We will show the calculation results based on these self-consistent dielectricdependent screened exchange potentials in the next section.

A. Acene molecules
First, we employed several acene molecules such as naphthalene, anthracene, tetracene, and pentacene to explain the fundamental behaviors of the dielectric-dependent screened exchange potential when applied to conjugated molecular systems.We consider n = 1-8, where n is the number of benzene rings, as shown in the inset of Fig. 1a.We investigated the electronic structure of these acenes based on the dielectric-dependent potentials given in eqn (1.1), (2.1), and (3).We also adopted the SVWN, B3LYP, and half-and-half (HandH) functionals for comparison purposes. 13,42hese functionals are described as follows: Here, V B88 is the Becke88 (B88)-type gradient correction. 43V LYP is the Lee-Yang-Parr correlation potential. 44The SVWN functional does not have any nonlocal exchange interaction, but the B3LYP and HandH methods contain 20% and 50% of the HF exchange term, respectively.We employ a dynamic-language framework for the calculations. 45The B3LYP/6-31G* method was used to obtain the optimized structures of these molecules.
The HOMO-LUMO gaps and the dielectric constants calculated using the 6-31G* and 6-311G* basis sets are summarized in Table 1.We can confirm that the HOMO-LUMO gap is narrowed when conjugated systems are extended.For example, in the case of naphthalene, B3LYP gave a HOMO-LUMO gap of 4.79 eV for the 6-311G* basis set.For the tetracene molecule, the HOMO-LUMO gap was reduced to 2.75 eV.When the dielectric-dependent screened exchange potential of eqn (1.1) was employed in conjunction with the 6-311G* basis set, the HOMO-LUMO gap for naphthalene was 6.16 eV and that for tetracene was 3.45 eV.From Table 1, it can be seen that the dielectric-dependent potentials of eqn (2.1) and (3) provide very similar results to eqn (1.1).Fig. 1a shows the HOMO-LUMO gaps calculated using the 6-311G* basis set as a function of the length of the conjugated system.Compared with the other methods, the SVWN method largely underestimated the HOMO-LUMO gaps because of the lack of nonlocal HF interaction.The HandH method yielded larger gaps, and the B3LYP method gave values intermediate between the SVWN and HandH methods.The dielectric-dependent approach gave similar HOMO-LUMO gaps to the HandH method for smaller unit lengths.However, at longer unit lengths, the HOMO-LUMO gaps calculated using eqn (1.1) gradually approached those obtained in the B3LYP calculations.Fig. 1b shows the HOMO and LUMO energies calculated using B3LYP, HandH, and the dielectric-dependent exchange potential approach of eqn (1.1).When the conjugated systems were extended, the HOMO (LUMO) energies gradually increased (decreased), and the HOMO-LUMO gaps narrowed.In the dielectric-dependent exchange potential approach, the HOMO and LUMO energies changed from values close to those obtained using the HandH method to values more closely corresponding to those generated by the B3LYP method.We will discuss these behaviors in the dielectric-dependent potential later from the viewpoint of the HF exchange ratio.
The energy differences between occupied and virtual (unoccupied) molecular orbitals appear in the denominator of the equation, as seen in eqn (4), so we can expect an inverse relationship between the HOMO-LUMO gap and the dielectric constant.Fig. 1c shows the calculated dielectric constants as a function of the unit lengths of the acene molecules.We note that the dielectric constants are automatically determined using the self-consistent procedure in the dielectric-dependent screened exchange potentials.In contrast, the dielectric constants for the SVWN, B3LYP, and HandH methods are estimated after SCF convergence.Fig. 1c shows that small systems such as benzene and naphthalene give smaller dielectric constants.When the conjugation lengths of the acene molecules were extended, an increase in the dielectric constants was observed.For example, in the case of naphthalene, the dielectric-dependent exchange potential of eqn (1.1) yielded a value of 2.51 for the dielectric constant at the 6-31G* level, as seen in Table 1.For tetracene, this value was 3.12.Thus, we can confirm that molecular systems with larger HOMO-LUMO band gaps give smaller dielectric constants.On the other hand, the SVWN and B3LYP methods calculated larger dielectric constants than those obtained using the dielectricdependent exchange potentials of eqn (1.1), (2.1), and (3).This means that the SVWN and B3LYP methods provide smaller HOMO-LUMO gaps than the dielectric-dependent exchange potential approach.In the dielectric-dependent exchange potential approach, the ratio of the exact HF exchange term is determined by the inverse of the dielectric constant value.We summarize the HF exchange ratio in Fig. 1d.In the case of the dielectric-dependent exchange potentials, the HF exchange ratio decreased when the conjugation lengths were extended.For example, eqn (1.1) gives a dielectric constant of 2.20, and therefore about 45% of the HF exchange interaction is employed for benzene at the 6-31G* level.In the case of tetracene, the dielectric constant and the HF exchange ratio were 3.12 and 32%, respectively.Conversely, the B3LYP method always employs 20% of the HF exchange term, which was empirically determined by Becke to reproduce the experimental data of several atoms and molecules. 13The HF exchange ratio of the HandH functional is 50%. 42Fig.1d shows that the dielectricdependent exchange potential approach employs a higher ratio for the partial HF exchange term at smaller unit lengths, which is closer to that of the HandH method; however, when the conjugation length is extended, the dielectric-dependent approach can adopt an HF exchange ratio similar to the B3LYP method.In the n = 8 case shown in Fig. 1d, the HF exchange ratio of eqn (1.1) was reduced to about 0.23, where the dielectric-dependent exchange approach gave similar calculation results to the B3LYP method.
In this paper, we employed acenes to simply explain basic behaviors of the dielectric-dependent approach according to the conjugation size.Here, we discuss several studies.Rangel et al. reported the G 0 W 0 and evGW calculations for acene crystals. 46The G 0 W 0 method depends on the DFT starting point, but the evGW method can partly compensate for the starting-point effect by using the self-consistent scheme.They reported the change in the gap properties according to the conjugation size.For example, the G 0 W 0 and evGW method yields the fundamental gap properties of 4.0 and 4.5 eV for anthracene crystals, and 2.9 and 3.5 eV for tetracene crystals, respectively.The dielectric-dependent method gives the fundamental gap of 3.91 and 2.80 eV for anthracene and tetracene crystals, respectively. 29The dielectric-dependent approach is closely related to the GW method, 12,14,34 and can give comparable results to the more expensive GW method. 24On the other hand, Stein et al. calculated the fundamental gap of acenes and Si nanocrystals based on the range-separated hybrid functional within the generalized Kohn-Sham approach. 47The rangeseparated hybrid DFT function possesses a system-dependent parameter.They showed that the system-dependent parameter is essential to consider the quantum size effect such as the conjugation length of acene oligomers.When the electron delocalization increases with the system size, the necessary weight of exact exchange term becomes smaller.In this paper, we also presented the change of the HF exchange term ratio in dielectric-dependent potentials according to the conjugation length.The behavior of the long-range HF exchange term discussed in this paper is similar to that of Stein's approach, although the systemdependent parameters used in two methods are different.

B. Organic photocell materials
Organic photocell devices have heterojunctions that are composed of donor and acceptor materials.For the acceptor region, fullerene derivatives such as phenyl-C 61 -butyric acid methyl ester are frequently employed.In the donor region, various types of molecules and polymers have been explored to improve the energy conversion efficiency of the organic photocell.In this paper, we examine the three materials depicted in Fig. 2 48 We investigate the electronic structure of these materials using the dielectric-dependent screened exchange potentials of eqn (1.1), (2.1), and (3) and the 6-31G* and 6-311G* basis sets.We also adopt the SVWN, B3LYP, and HandH methods for comparison.The optimized structures used in our previous study were employed here. 49We summarize the calculation results for TPD, BTz, and NTz in Tables 2, 3, and 4, respectively.
In the organic photocell materials, the HOMO-LUMO gaps decreased as their oligomer lengths were extended.Conversely, at longer conjugation (unit) lengths, those molecules gave smaller dielectric constants because of the inverse relationship between the magnitude of the HOMO-LUMO gap and the dielectric constant.For example, the HOMO-LUMO gap and the dielectric constant of the TPD monomer were 3.45 eV and 3.16, respectively, based on eqn (1.1) with the 6-31G* basis set, as seen in Table 2. Trimeric TPD exhibited a HOMO-LUMO gap and a dielectric constant of 2.52 eV and 3.75, respectively.Thus, the conjugation (unit) lengths of these materials strongly influenced both the HOMO-LUMO gaps and the dielectric constants.Conversely, at large unit lengths, the HOMO-LUMO gaps and the dielectric constants become saturated.For example, the HOMO-LUMO gaps for trimeric and tetrameric NTz were 2.09 Fig. 2 Thiophene-based polymers used in this study.R is replaced with a methyl group.and 2.04 eV, respectively, when eqn (1.1) was employed together with the 6-311G* basis set, as seen in Table 4.The B3LYP method and the dielectric-dependent exchange potential approach gave similar calculation results.It should be noted that, compared with both the B3LYP method and the dielectric-dependent approach, the SVWN method underestimated the HOMO-LUMO gap and overestimated the dielectric constant.This occurs as a result of the lack of nonlocal HF exchange interaction in the SVWN method.Conversely, the HandH method overestimated the gap and underestimated the dielectric constant.
We summarize the HOMO-LUMO gaps, HOMO and LUMO energies, and dielectric constants of the organic photocell materials in Fig. 3a, b, and c, respectively.In Fig. 3d, we show the inverse values of dielectric constants, which represent the HF exchange ratio for each target system.To generate these figures, we employed eqn (1.1) with the 6-311G* basis set.Fig. 3a shows that the HOMO-LUMO gaps decreased when the unit lengths were extended.The HOMO and LUMO energies changed according to the unit length, as seen in Fig. 3b, but those influences were saturated in the tetramer cases.In addition, the HOMO-LUMO gaps narrowed when the conjugation lengths were extended.Fig. 3c shows that the values of the dielectric constants increased in parallel with the increase in the conjugation (unit) lengths.However, the calculated dielectric constant values were different and depended on the type of material and the unit length.For example, the TPD monomer exhibited a dielectric constant of 3.21 in eqn (1.1) with the 6-311G* basis set.Conversely, for the NTz trimer, this value was estimated to be 4.00.Thus, organic photocell materials can take different dielectric constant values.These different dielectric constants give different HF exchange ratios in the dielectric-dependent exchange potentials, as seen in Fig. 1d.Next, we discuss comparisons between theoretical calculations and experiments.For organic materials, reliable experimental measurements of the HOMO-LUMO gap may be more difficult than analogous experiments on inorganic bulk materials.In the case of polymeric materials, changes in the polymer type and/or experimental conditions can lead to the formation of several different crystal structures, which may include amorphous and semi-crystalline states.This disturbance and disorder in a crystal structure causes the localization of electrons in the organic material, and may make it more difficult to accurately estimate the HOMO-LUMO gap.In addition, the dielectric constants of organic materials are usually lower than those of inorganic materials.Therefore, when a photon is absorbed by an organic photocell material, a strongly bound exciton is generated instead of a free electron and hole.Thus, the light absorption energy (optical gap) is usually lower than the HOMO-LUMO gap energy; the difference corresponds to the exciton binding energy.The HOMO-LUMO gap expresses the energy difference between a free electron and a hole.Thus, we need some correction to estimate the HOMO-LUMO gap from photon absorption experiments of organic materials.Vanossi et al. reported the correction of DE EC = 0.53 + 0.99DE opt from electrochemical and optical experimental measurements on several organic photocell materials, 50 where DE opt is the optical bandgap and DE EC is the electrochemically determined bandgap property; the latter parameter is more suitable to compare with the theoretical HOMO-LUMO gap.In Table 5, we show the optical gaps measured by Takimiya et al. 48and the corrected experimental values calculated using the Vanossi method. 50Also presented are theoretical calculations using the 6-311G* basis set for the tetrameric materials.It can be easily seen that there are large differences between the optical gap and the theoretical HOMO-LUMO gap resulting from the strong exciton binding energy.Conversely, the Vanossi correction gives values closer to the theoretical calculations.For example, in the NTz case, the optical band gap was found to be 1.56 eV.The corrected experimental gap was 2.07 eV.The dielectric-dependent exchange potential of eqn (1.1) gave a HOMO-LUMO gap of 2.04 eV.The B3LYP method produced a similar calculation result of 1.84 eV, although it slightly underestimated the experimental value.Conversely, the larger HF exchange ratio employed in the HandH method meant that the HOMO-LUMO gap was overestimated to be 3.33 eV.The SVWN method yielded a HOMO-LUMO gap of 0.94 eV; this method largely underestimated the experimental values because of the lack of the HF exchange term.Thus, the DFT methods gave different values according to the HF exchange ratio combined with the DFT functional.The B3LYP method, in which a suitable HF exchange ratio was used for the TPD, BTz, and NTz materials, resulted in better predictions for the HOMO-LUMO gap.The self-consistent dielectric-dependent approach can automatically determine a suitable HF exchange ratio according to the target system.Thus, we have shown that the dielectric-dependent exchange approach can produce good agreement between theoretical and experimental results.On the other hand, it has also been demonstrated that the physical concept of the dielectric-dependent approach, where the HF exchange ratio becomes the inverse value of the dielectric constant, may be helpful in selecting an appropriate DFT functional to investigate organic photocell materials.
The dielectric-dependent approach can reproduce experimental gaps well, for the fullerene molecular crystal, 29 and for other molecular crystals and semiconductors. 21,30Thus, the validity of the dielectric-dependent approach for various systems has been investigated by several research groups. 19,21,24,34onversely, it is not a simple matter to consider polymer materials, because of their ''dirty'' material quality.The experimental absorption energies summarized in Table 5 are for the polymer materials.This means that the data are modified by various effects, such as electrostatic and polarized interactions among polymers, deformations and amorphous structures of polymers, impurities, and so on. 51In simulation studies, it may be difficult to consider all of those.So, this paper focused on the conjugation length, because it is one of the most influential effects for the gap property.Polymer materials usually are composed of a plurality of conjugated parts.Therefore, the relationship between the conjugation size and the band gap property is an important factor to investigate actual polymers.In the NTz case, the monomer unit gave the HOMO-LUMO gap   5.For other organic photocell materials, we can observe similar behaviors.These results suggest that the trimer and tetramer units can be used to discuss more complex effects on the gap property.On the other hand, Rafaely-Abramson et al. discussed the gap renormalization of molecular crystals based on the DFT approach, 20 and Neaton et al. studied the renormalization at metal-molecular interfaces based on the GW method. 52These theoretical considerations may be useful for organic photocells.
Here, we discuss the HOMO-LUMO gap and the fundamental gap. 47,504][55] Thus, the HOMO-LUMO gap tends to underestimate the fundamental gap properties, especially in small molecular systems.The DFT-based many body perturbation theory (MBPT) may be used to correct such errors in the molecular orbital energies. 56Conversely, in the infinite limit under the periodic boundary condition, the energy difference between quasiparticles is frequently employed to estimate the fundamental gap. 36This means that the evaluation based on quasiparticle energies can give fairly acceptable results for the gap property, when molecular systems become larger and orbitals more delocalized. 57Thus, the HOMO-LUMO gap is suitable to evaluate the fundamental gap of organic photocell materials with large conjugated systems.
The fundamental gap can be estimated from the total energy difference of neutral and ionized molecules in the DSCF scheme.In other words, the DSCF approach is based on the total energy expression.Conversely, this paper focused on the quasiparticle description, because it is essential to understand the energy conversion mechanism of the organic photocell device.It should be noted that the HOMO-LUMO gap gives a good estimation for the fundamental gap in large conjugated systems.Therefore, the quasiparticle description is often employed to analyze organic photocell materials. 48,58The dielectric-dependent screened exchange potentials were originally developed to improve the quasiparticle description, and therefore this paper paid more attention to it.In addition, the description is required to develop more realistic theoretical models for organic photocells. 8,9Accurate quasiparticles can give a better starting point for advanced methods such as the MBPT calculation.The dielectric-dependent quasiparticle description will be useful to investigate organic photocell materials further.

IV. Summary
This paper focused on the fundamental behaviors of the dielectric-dependent screened exchange potential when applied to conjugated molecular systems.In particular, the HOMO-LUMO gap of the donor material of an organic photocell device was investigated, because this gap is an essential factor in the development of organic photosolar devices.When conjugation systems are extended, the HOMO-LUMO gaps decrease, but the dielectric constants increase because of the inverse relationship between these properties.In the dielectric-dependent exchange potential approach, the ratio of the exact HF exchange term combined with the potentials is determined by the inverse of the dielectric constant.Therefore, when the conjugation (unit) lengths are extended, the HF exchange ratio decreases in the dielectric-dependent approach.The molecules and materials examined in this paper take different dielectric constant values in accordance with the material types and unit lengths; thus, the HF exchange ratios are also different among the dielectricdependent exchange potentials.We also compared our calculated results with experimental measurements.The theoretical HOMO-LUMO gap usually overestimates experimental optical gaps as a result of the strong binding energy between an electron and a hole (exciton) in organic materials; however, if we correct for the exciton binding energy, we can observe better agreement between experimental and calculated values.To develop organic photocell materials, information regarding the gap property is essential.
The dielectric-dependent exchange potential approach will be useful to investigate organic photocell materials.Additional studies in this area are ongoing in our laboratory.

Fig. 1
Fig. 1 (a) HOMO-LUMO gaps, (b) HOMO and LUMO energies, (c) dielectric constants, and (d) HF exchange ratios combined with DFT functionals with respect to unit lengths for acene molecules.The white triangles show the Half and Half (HandH) method, and black ones indicate the SVWN method.The black circles express the B3LYP method.The white circles represent the dielectric-dependent screened exchange potential of eqn (1.1).

Fig. 3
Fig. 3 (a) HOMO-LUMO gaps, (b) HOMO and LUMO energies, (c) dielectric constants, and (d) HF exchange ratios combined with the potential of eqn (1.1) in accordance with unit lengths for organic photocell materials.The black and white circles show calculation results for TPD and BTz materials, respectively.The black triangles represent the NTz material.

Table 1
Dielectric constants and HOMO-LUMO gaps [eV] for acene molecules This journal is © the Owner Societies 2016

Table 2
Dielectric constants and HOMO-LUMO gaps [eV] for TPD units

Table 3
Dielectric constants and HOMO-LUMO gaps [eV] for BTz units

Table 4
Dielectric constants and HOMO-LUMO gaps [eV] for NTz units

Table 5
Comparisons between experimental and theoretical gap energies [eV] for organic photocell materials ) with the 6-311G* basis set, and the trimer and tetramer units yielded 2.10 and 2.05 eV, respectively.The experimental value corrected by the Vanossi method was 2.07 eV, as shown in Table