Computational Prediction of Electronic and Photovoltaic Properties of Anthracene-Based Organic Dyes for Dye-Sensitized Solar Cells

Three kinds of anthracene-based organic dyes for dye-sensitized solar cells (DSSCs) were studied, and their structures are based on a push–pull framework with anthracenyl diphenylamine as the donor connected to a carboxyphenyl or carboxyphenylbromothiazole (BTZ) as the acceptor via an acetylene bridge. The photoelectric properties of the three dyes were investigated using density functional theory (DFT). The simulations indicate that the improvement of anthracene-based dyes (the addition of BTZ and the change of alkyl groups to alkoxy chains) can reduce the energy gap and produce a red shift. This structural modification also improves the light capturing and the electron injection capability, making it excellent in photoelectric conversion efficiency (PCE). In addition, twelve molecules have been designed to regulate photovoltaic performance.


Introduction
With the depletion of traditional fossil fuels and environmental pollution, green energy has aroused widespread concern in academia [1]. Therefore, nonpolluting solar energy has become the most promising alternative energy source [2]. Compared with traditional inorganic solar cells based on silicon crystal, dye-sensitized solar cells (DSSCs) have the advantages of easy synthesis, low cost, and high conversion efficiency [3,4]. Since the first report in 1991, DSSCs have a high PCE [5]. In general, a typical DSSC device consists of a titania semiconductor film, a dye sensitizer, a redox electrolyte, a counterelectrode, and a transparent conductive substrate [6][7][8]. The dye is mainly divided into metalcontaining ruthenium dyes [9], porphyrin dyes [10], and metal-free organic dyes [11]. As an important part of DSSCs, sensitizers play an important role in capturing sunlight and the electron transfer. Among them, ruthenium (II) polypyridyl complexes are considered to be efficient and stable sensitizers with a power conversion efficiency (PCE) above 11% at AM1.5G [12,13]. However, the scarcity, high cost, and toxicity of ruthenium metal limit the widespread use of such sensitizers in DSSCs. In addition, the zinc porphyrin dye is more than 12% efficient in Co II /Co III electrolytes under standard conditions, which is considered to be a very promising sensitizer [14,15]. In recent years, perovskite solar cells have become another potential photovoltaic approach with efficiency of over 20% under AM1.5G light sources and dim light irradiation [16,17]. However, solution of the instability of devices and pollution for environment caused from raw materials are still a challenge [18][19][20]. Metal-free organic dyes are characterized by low cost, ease of purification, and flexible molecular design [21], and metal-free sensitizers are designed with donor-π-acceptor (D-π-A), D-π-A-A, or D-A-π-A, which can lead to light-induced charge separation, improvement in stability, and optimization of the energy level of the dye from structure modifications [22][23][24][25]. To date, PCE of such kind of organic dyes has reportedly reached 14% [26]. Due to a large number of functional groups available for molecular design, there still is much work for DSSC performance improvement.
Molecular materials with the anthracene structure have good stability and special luminescent properties, showing a bright blue electroluminescence [27][28][29]. However, there are relatively few studies on the photoelectric properties of anthracene-based molecules. There are few metal-free sensitizers featuring a 9,10-disubsituted anthracene entity as a conjugated spacer between electron donor and acceptor moieties, and the best PCE is 7.03% [30][31][32][33]. Recently, a dye containing 2,6-conjugated anthracene showed a photoconversion efficiency of 9.11% [34] at 1 sun condition, which is the highest PCE reported by anthracene dyes. Mai and coworkers [35] have synthesized a simple D-π-A sensitizer (MS3), which is based on a 9,10-disubstituted anthracene entity with an optical efficiency of 5.84% at AM1.5G. Based on the MS3 sensitizer, TY3 (D-π-A) and TY6 (D-A-π-A) were also synthesized. TY6 has the best PCE (up to 8.80%) [36]. To study the relationship between structure and properties, we used the theory of density functional theory (DFT) and time-dependent density functional theory (TD-DFT) to calculate the three-molecular-geometry, electron injection, dye regeneration, and optical properties, and results confirmed that its excellent performance was due to its excellent J SC and V OC characteristics. In addition, a series of design molecules based on TY6 were investigated to improve optical response and electron injection.

Computational Details
The ground-state geometries of three molecules were optimized by DFT//B3LYP/6-31G(d) level [37][38][39][40]. In order to simulate the more realistic performance of dye-sensitized solar cells, the related calculations were performed in the solvent condition (THF) by using the C-PCM [41] method. Frequency calculations showed the minima on the potential energy surface for optimization. The bond lengths, dihedral angles, energy gaps, frontier molecular orbitals, electron injection, and recombination of the optimized molecules were calculated. The absorption spectra, transition energies, and oscillator strengths of molecules were obtained with TD-DFT [42] by using the CAM-B3LYP [43] functional at the same basis set as the ground state. Three excited states are calculated, including the first excited state (S1), the second excited state (S2), and the third excited state (S3). The natural bond orbital (NBO) analysis [44] for the charge difference between the ground state and the excited state was carried out at the B3LYP/6-31G(d) level using the NBO 3.1 program. By introduction of an electron-withdrawing group in the acceptor, it hoped that the molecular modification can attract electrons and promote an intramolecular charge transfer from donor to acceptor, further leading to better electron injection into the conduction band of TiO 2 . Therefore, twelve new dye molecules were designed by introducing CN, F, and CF 3 into the acceptor of the TY6 molecule, and the correlation calculations were made using the same method as the original dye molecule. All calculations are made through the Gaussian 09 package [45]. Figure 1 shows the optimized ground-state molecular structure. As shown in Figure 1(a), MS3 is an original molecule, and TY3 is obtained by converting the C-6 alkyl chain of MS3 to the carbon alloy group. Based on the TY3 molecule, benzotriazole (BTZ) was introduced between the acetylene bridge and benzoic acid. In order to reduce the aggregation of sensitizer and improve the performance, a long alkyl chain was introduced into the N position of BTZ to obtain TY6. Figure 1(b) shows the ground-state structures of the three optimized molecules in the THF solvent. In order to facilitate the calculation, the long carbon chains on the donor were pruned appropriately. Both MS3 and TY3 are typical D-π-A structures; the donor and the π-bridge are the amino donor and the acetylene bridge, and the acceptor is benzoic acid. TY6 is the D-A-π-A structure, in which an additional acceptor is added between the acetylene bridge and benzoic acid. Table 1 shows the bond length and dihedral angle in gas and solvent (THF), respectively. For example, in gas, the dihedral angles of MS3 ∠C1-C2-N3-C4 and ∠C2-N3-C4-C5 were 33.16°and 70.94°, and the average value of the two warped dihedral angles is 52.05°. In the same way, the calculated mean values of the donor dihedral angles for TY3 and TY6 are 52.60°and 52.40°. It shows that the donor has a distorted structure. Dyes TY3 and TY6 are more largely distorted than the original molecule MS3, which can reduce the aggregation of dye molecules. At the same time, the three molecules have similar results in solvent: TY3 (52.66°) and TY6 (52.33°) are greater than the donor dihedral angle of MS3 (51.88°). In gas, the bond lengths of the three molecules N3-C4 and C10-C11 are smaller than those of the single bond (i.e., C-C: 1.530 Å [46], N-C: 1.471 Å [47]). In general, the stability of a molecule can be judged by the length of the bond [48]. The shorter the bond length is, the more stable the molecule becomes. By comparison, the bond lengths of TY3 and TY6 are less than that of MS3, and the same results are also observed in solvent condition; therefore, TY3 and TY6 molecules are more stable. Table 2 shows the HOMO, HOMO-1, LUMO, LUMO + 1, and energy gaps (Δ H−L ) for MS3, TY3, and TY6. In the gas phase, the HOMO order of the molecule is as follows: TY6 (−4.59 eV) > TY3 (−4.67 eV) > MS3 (−4.91 eV). LUMO energy is arranged in the following order: TY3 (−2.26 eV) > TY6 (−2.30 eV) > MS3 (−2.33 eV). The energy gap is arranged in the following order: MS3 (2.58 eV) > TY3 (2.41 eV) > TY6 (2.29 eV). As shown in Figure 2, we can find that LUMO of the three molecules has little change, and thus, the gradual increase in the HOMO energy level leads to the decrease in energy gap. The reduction in the energy gap also favors the red-shifted absorption peak in the UV absorption spectrum (data in Table 3 and spectra nature in Figure 3). Therefore, the absorption peak of TY6 will have a significantly red shift.

Energy Levels.
In solvent condition, for MS3 the energy of HOMO (−5.01 eV) is reduced by 0.10 eV compared with that in the gas condition, and the energy of LUMO (−2.44 eV) is 2 International Journal of Photoenergy  International Journal of Photoenergy reduced by 0.11 eV compared with that in the gas condition. Therefore, the energy gap (2.58 eV) in solvent is lower than the energy gap in gas (2.57 eV). Similarly, dyes TY3 and TY6 also show the same change, and the results in solvent are better than those in the gas phase, and the energy gap is arranged in the following order: MS3 (2.57 eV) > TY3 (2.39 eV) > TY6 (2.28 eV). TY6 still has a smaller energy gap. The frontier molecular orbitals of three molecules are shown in Figure 4, which indicates the HOMO and LUMO of the three molecules and the distribution area of the electron density. As a whole, the electron density of HOMO is mainly located near the donor and the π-bridge, and the electron density of LUMO is mainly located near the π-bridge and the acceptor. HOMO-1 is lower than HOMO, and the electron density is dispersed in Figure 2. LUMO+1 is higher than LUMO, and its electron density is mainly concentrated on the π-bridge and the acceptor. Therefore, it can be concluded that the intramolecular charge transfer (ICT) exists in the dye molecule when the electron is transferred from the donor to the acceptor. Figure 2 shows the orbital energy levels of three molecules. In the gas phase, the HOMO energy (−4.91 eV) of MS3 is lower than the energy of I − /I − 3 (−4.85 eV), while the HOMO energy of TY3 (−4.67 eV) and TY6 (−4.59 eV) is higher than that of I − /I − 3 . This shows that TY3 and TY6 can more easily recover electrons from electrolytes. The energies of LUMO of three dye molecules are lower than the conduction band energy of TiO 2 (−4.00 eV), which indicates that electrons can be successfully injected into TiO 2 from the excited state of dye molecules. One of the most important characteristics in the excellent ICT is charge separation. The charge distribution of HOMO and LUMO can promote the transfer of electrons (see Figure 4). In order to investigate the characteristics of ICT, we used the charge density difference (CDD) to show the change of the electron density between the ground state and the excited state. It can characterize CT in organic molecular systems [49,50], which clearly reflects the direction of the electrons. Taking the first excited state of MS3 as an example, the donor is covered with green holes, and the acceptor is covered with red electron (see Figure 3). Hole and electron are alternately distributed in the π-bridge moiety, and TY3 also shows similar results. For TY6, the S1 of CDD showed electron transfer is from donor to π-bridge, and for S2 week electron move to auxiliary acceptor BTZ and acceptor benzoic acid; wherefore, the S3 of CDD showed that the green holes make a shift from the donor to the π-bridge, and the red electrons transfer from the π-bridge to auxiliary acceptor BTZ and acceptor benzoic acid. The results show that there is a significant charge separation between the donor and the acceptor, which is considered as an ICT feature. Table 3 lists the calculated absorption peaks, transition energies, and oscillator strengths (only discuss the state of f > 0 1). In the UV-Vis spectral region, the main absorption band was found to be the first excited state (S1 state). The first excited state (S1) of MS3 corresponds to the electron transition from HOMO to LUMO. It can be seen from Figure 4 that electrons transfer from the amino donor to the benzoic acid acceptor. In gas, the maximum absorption peak is 452 nm (463 nm in solvent), and the oscillator strength is 0.5880 (0.7021 in solvent). For the higher excited state (S2), the absorption intensity is lower than in S1, and the maximum absorption peak is 375 nm (f = 0 1800), which shows the electron transition process from HOMO-1 to LUMO. For TY3, the maximum absorption peak in gas S1 is 470 nm (481 nm in solvent), and the oscillator strength is 0.5112 (0.6182 in solvent), which shows the electron transition process from HOMO to LUMO. For the second excited state (S2), the maximum absorption peak is 384 nm (f = 0 2814), which shows the electron transition process from HOMO-1 to LUMO. For TY6, the maximum absorption peak of S1 in gas is 481 nm (494 nm in solvent). The oscillator strength is 0.8886 (0.9799 in solvent), which indicates the electron transition process from HOMO to LUMO. For the higher excited state (S2 and S3), the maximum absorption peak of S2 is 395 nm (f = 0 3714), showing an electron transition from HOMO-1 to LUMO. The maximum absorption peak of S3 is 344 nm (f = 0 1763), which shows the electron transition process from HOMO to LUMO+1. According to the data in Table 3, the maximum absorption peak value of the three molecules in gas is as follows: TY6 (481 nm) > TY3 (470 nm) > MS3 (452 nm). The UV-Vis spectra are given in Figure 5. Compared with   (29 nm). This result indicates that the reduction in the energy gap favors the red-shifted absorption peak. The peak ranges are between 500 and 700 nm, which is helpful for effectively absorbing sunlight. It can be seen from Figure 5 that TY6 has the most pronounced red-shift absorption with the highest molar absorption coefficient. From the gas to solvent condition, absorption peaks of the red shift is approximately 11 nm, and the molar extinction coefficient is also increased ( Figure 5). Performance in solvent is better than in gas. In summary, TY6 absorbs sunlight more efficiently than other molecules do, which may lead to higher PCE.

Chemical Reactivity Parameters.
Ionization potential (IP) and electron affinity (EA) are important data for measuring the injection ability of holes and electrons [51][52][53]. The calculated IPs and EAs are listed in Table 4. In gas, the IPs of the three molecules are arranged in the order TY6 (5.49 eV) < TY3 (5.65 eV) < MS3 (5.89 eV). In addition, the reduction in IP means the increase in contribution ability. In gas, the EAs of the three molecules are arranged in the order TY6 (1.40 eV) > MS3 (1.34 eV) > TY3 (1.22 eV). The greater the value of EA becomes, the stronger the ability to receive the electronic have [54]. Therefore, TY6 has a higher ability to accept electrons. In solvent, the IP ordering of the three molecules is consistent with that in gas, and each  molecule decreases by about 0.9 eV. At the same time, EA increased by 1.2 eV compared with the gas phase. In general, the effect in solvent is better than that in gas. It is obvious that the electron transport performance of TY6 is better than that of TY3 and MS3. Table 5 lists the chemical hardness (h), electrophilicity (ω), and electroaccepting power (ω + ). h represents the strength of resistance to ICT [55,56]; a small h is helpful in reducing resistance to ICT. ω + represents the higher ability to accept electronics [57]. Many studies have shown that lower h and higher ω + can lead to higher short-circuit currents. In gas, h is arranged in the order MS3 (3.62 eV) > TY6 (3.44 eV) > TY3 (3.43 eV); ω + is arranged in the following order: TY6 (1.43 eV) > MS3 (1.35 eV) > TY3 (1.22 eV). There was no significant difference between TY3 (3.43 eV) and TY6 (3.44 eV) in h. However, TY6 (1.43 eV) was significantly higher than MS3 (1.35 eV) in terms of ω + . In solvent, the orders of h and ω + are the same as those in the gas phase. However, ω + has a significant improvement, and MS3 increases by 2.75 eV, TY3 increases by 2.69 eV, and TY6 increases by 3.30 eV. The promotion of TY6 is the most obvious, and it represents a higher receiving capacity. ω represents the stability of the dye molecule system [58]. In gas, ω is arranged in the following order: TY6 (2.89 eV) > MS3 (2.87 eV) > TY3 (2.66 eV). In solvent conditions, MS3 increased by 2.94 eV, TY3 increased by 2.90 eV, and TY6 increased by 3.53 eV. It is clear that TY6 has not only the highest ω but also the highest increase in solvent. Therefore, TY6 has the highest energetic stability. MS3 (S1) TY3 (S1) TY6 (S1) TY6 (S2) TY6 (S3) Figure 3: The charge difference density (CDD) for MS3, TY3, and TY6 in solvent (green and red stand for the hole and electron, respectively). 6 International Journal of Photoenergy incident solar energy (P in ). Calculated efficiency can be written as follows [59]: From the formula, it can be seen that high V OC and J SC are the basis for producing photoelectric conversion efficiency. The J SC in DSSCs can be calculated by the following equation: where LHE λ is a light harvesting efficiency at maximum wavelength, Φ inject is defined as the dye molecule exciting electron injection efficiency, and η collect is the charge collection efficiency. LHE can be expressed as [60] where f is the oscillator strength of the dye molecules; a large oscillator strength can contribute to the improvement of LHE. The electron injection-free energy (ΔG inject ) can be represented as [61] E CB is the reduction potential of the TiO 2 semiconductor and is equal to 4.0 eV [62] (versus vacuum) in this work. E dye * OX is the oxidation reduction potential of the dye in the excited state. E dye * OX can be expressed as [63] E dye * where E dye OX is the oxidation potential energy of the dye in the ground state, while E 00 is the electronic vertical transition energy corresponding to λ max .
The calculated values of LHE and ΔG inject are listed in Table 6. ΔG inject is an important factor affecting the electron    Higher LUMO can lead to a higher ΔG inject (absolute value). As shown in Table 2, LUMO of TY3 (−2.26 eV) and TY6 (−2.30 eV) is higher than that of MS3 (−2.33 eV). The corresponding ΔG inject is higher than that of MS3. Therefore, the improvement in DSSC's PCE is due to the high ΔG inject , further resulting in high electron injection efficiency. The higher ΔG inject is good for increasing the J SC of the DSSC. In solvent conditions, ΔG inject has a small change. However, TY6 still has the highest ΔG inject (absolute value). To sum up, TY6 will have higher J SC . Three molecules of light harvesting efficiency can be obtained from Table 6. In gas, the data is sorted in the order TY6 (0.871) > MS3 (0.742) > TY3 (0.692). The introduction of the side chains on BTZ will increase LHE, and a larger oscillator strength will lead to a higher LHE. From Table 3, it can be found that the oscillator strength of the three molecules in solvent was higher than that in gas. So the LHE in solvent is greatly improved, and TY6 (0.895) still has the highest LHE. Therefore, under the excitation of light, TY6 has higher solar light utilization, producing more photocurrent, and this result indicates that TY6 will have higher J SC Dye regeneration free energy (ΔG regen dye ) is a significant factor affecting photoelectric conversion efficiency which can be written as follows [64]: is the redox potential of I − /I − 3 (−4.85 eV) [65]. Table 6 shows ΔG regen dye of three molecules. In gas conditions, the data is sorted in the order TY6 (0.26 eV) > TY3 (0.18 eV) > MS3 (−0.06 eV). In solvent conditions, the decrease in HOMO results in the increase in ΔG regen dye (see Table 2). So, it has a great influence on TY6, which has the largest ΔG regen dye (0.08 eV). The larger ΔG regen dye can promote dye regeneration and increase J SC . This means that TY6 will have a better performance.
Another important factor affecting the efficiency of electron transfer is the lifetime (t) of the first excited state. If the molecule has a longer life span, it will contribute to the charge transfer of the molecule [66]. The lifetime (t) can be obtained by the latter formula: t = 1 499/ f E 2 , where f is the oscillator strength and E (cm −1 ) is the excitation energy of the different electronic states [67]. Table 6 shows the corresponding data, and the result in the gas phase is in the order TY3 (6.48 ns) > MS3 (5.20 ns) > TY6 (3.91 ns). In solvent, the order does not change, but the gap of them decreases. TY6 has a lower t, which is owing to the fact that TY6 has a very high f ( Table 3). As known, there are many complicated factors that affect J SC . According to (5), although TY6 has a lower t, it has the highest LHE and Δ G inject , and it does not affect the highest J SC . As for V OC in DSSCs, it is described by [68] Here, E CB is the conduction band edge of TiO 2 , q is the unit charge, KT is the thermal energy,n c is the number of electrons in the CB, N CB is the accessible density of CB states, and E redox is the redox potential of the electrolyte. Equation (7) shows that the energy E CB and the number of electrons in the CB are important factors affecting V OC . ΔE CB is the displacement of CB when the dye absorbs on the surface of TiO 2 , which can be expressed as [69] ΔE CB = qμ normal γ εε 0 8 γ is the concentration of dyes in the surface, μ normal is the dipole moment component perpendicular to the direction of the TiO 2 surface (where μ normal is the x-axis direction), ε is the dielectric constant of the organic monolayer, and ε 0 is the dielectric constant of the vacuum. Figure 6 shows the values of the vertical dipole moment (μ normal ). TY3 (5.7641 D) 3.6. Natural Bond Orbital Analysis. In order to understand the mechanism of photoexcitation, we simulated the natural bond orbital (NBO) of the optimized structure of the ground state (S0) and the first excited state (S1) of MS3, TY3, and TY6. The acceptor of TY6 includes benzoic acid and the auxiliary acceptor (BTZ group), and the detailed data is listed in Table 7. The amount of charge difference (Δ q ) from S0 to S1 in the donor group shows that TY3 (−0.3151) can contribute more electrons compared to MS3 (−0.2808) and TY3 (−0.2749). But Δ q in the anthracene group and acceptor groups shows that more of the contributing electrons in TY3 are concentrated in the anthracene group. So, only a small amount of electrons reach the acceptor part in TY3. Eventually, the acceptor group of TY6 has more electrons. It showed that the auxiliary acceptor (BTZ) in TY6 increased the total amount of electron in the acceptor. At the same time, Δ q of MS3 (0.0104), TY3 (0.0190), and TY6 (0.0103) in the π region indicates that the π-linker is only a channel for the charge transfer. To sum up, TY6 can provide more efficient excitation electrons in the photoexcitation mechanism.

Hyperpolarizabilities and Reorganization Energies.
The total static first hyperpolarisability can be expressed as [70] β tot = β 2 The static component is calculated by the following equation: where β ijk i, j, k = x, y, z are tenser components of hyperpolarizability. Finally, the equation is written as For DSSCs, the ICT process facilitates the aggregation of electrons in the acceptor moiety, and the enhanced electron density in the acceptor moiety can enhance the electronic coupling effect between acceptor and semiconductor. The first β is directly proportional to the transition dipole moment ((oscillator strength) (μ eg )) and the difference in the dipole moment between the ground and excited orbitals Δμ eg , and it is inversely proportional to the transition energy (E eg ). β can be expressed as follows [71]: where Δμ eg and μ eg are difference in the dipole moment for ground state and excited state, and the transition dipole, E eg , is the transition energy. The first hyperpolarizabilities are listed in Table 8. β xxx (along the coordinate axis of the molecule) is a negative value, and the negative charge is far away from the nuclear charge of the molecule. β tot of TY6 is much larger than MS3 and TY3 values. Table 3 and Figure 6 support the results of hyperpolarizabilities. Figure 6 shows that TY3 and TY6 have higher dipole moments, and TY6 has higher oscillator strengths (in Table 3); so, TY6 has a higher μ eg . At the same time, TY6 has less excitation energy in Table 2. So, TY6 has larger hyperpolarizabilities with an obvious charge transfer. The reorganization energy can affect CT, which is more beneficial to the improvement of CT [72]. On the basis of the Marcus theory, CT can be calculated by [73] where λ is the reorganization energy, A is the electronic coupling, K B is the Boltzmann constant, and T is the temperature. Hole or electron reorganization energy is determined by the following equation [74]: where E 0 − (E + 0 ) represents the energy of the neutral molecule calculated at the anionic (cationic) state. E − 0 (E + 0 ) represents the anion (cation) energy calculated with the optimized structure of the neutral molecule. E − (E + ) represents the anion (cation) energy calculated with the optimized anion (cation) structure. E 0 represents the energy at the neutral molecule at the ground state. λ h , λ e , and λ total values are listed in Table 9. In the gas phase, the λ h /λ e values of MS3, TY3,  The results show that the λ h values of three molecules are lower than those of λ e , which indicates that the electron transfer rate is lower than the hole transfer rate. Obviously, MS3 has a higher hole and electron transfer rate, followed by TY6. From the gas phase to solvent, λ h and λ e have decreased. For λ total , MS3, TY3, and TY6 are reduced to be 0.05 eV, 0.10 eV, and 0.11 eV (see Table 9). The above data shows that values of λ h and λ e have been reduced in the THF solvent. TY6 has the greatest reduction in solvent and is very close to MS3 (0.40 eV). As a result, MS3 and TY6 will lead to better hole/ electron transport and efficient luminescent materials.
3.8. Analysis of Electrostatic Potential Distribution on the Molecular Surface. In order to determine the position of electrolyte ions, we can determine the reactive sites of dye molecules by the molecular surface electrostatic potential (ESP) [75]. The ESP of three molecules is listed in Figure 7 In order to show the molecular surface area of different electrostatic potential intervals, the quantitative distribution chart of the electrostatic potential on the surface of molecules is plotted, which is listed in Figure 7(b). As shown, the distribution of the surface electrostatic potential at the maximum and minimum points of the three molecules is very small in this region. The electrostatic potential area of the larger region is as follows: MS3 (−17.83, 13.83 kcal/mol), TY3 (−18.07, 13.27 kcal/mol), and TY6 (−24.67, 18.27 kcal/mol). It seems that the electrostatic potential distribution of TY6 is wider and more homogeneous than that of MS3 and TY3. Therefore, the surface reaction region of TY6 is the largest, and the overall reaction activity is stronger.
3.9. Molecular Design. After the above discussion, the TY6 molecule has the best photoelectric conversion efficiency. On the basis of this molecule, a series of molecules were designed by inserting different electron-withdrawing groups (EWGs (−CF 3 (A), -CN (B), and -F (C))) into number 1, 2, 3, and 4 positions of the original molecule acceptor, respectively. The designed molecular structures are shown in Figure 8, and those dyes are named as TY6-X (X = 1A, 2A, 3A, 4A, 1B, …,4C). The following calculation results are obtained in the gas phase. The calculated bond lengths and dihedral angle are listed in Table 10. The dihedral angle (∠C1-C2-N3-C4 and ∠C2-N3-C4-C5) on the donor is not significantly changed compared with the original molecule, and bond length has almost no change. However, the dihedral angle (∠C6-C7-C8-C9) of the acceptor is changed obviously, and the change of the molecule with the insertion of the -CN group has the largest change. It is worth noting that the dihedral angles of TY6-1X and TY6-4X (∠C6-C7-C8-C9) have a significant increase, while TY6-2X and TY6-3X have no significant changes. It is due to the fact that the close distance between the introduction of EWGs and the long alkyl chain introduced on BTZ will result in mutual exclusion, and TY6-1X (X = A, B, C, D) is larger than TY6-4X (X = A, B, C, D). For example, the dihedral angles of TY6-1A and TY6-4A are 65.82°and 49.99°, respectively. The twisted structure of the donor blocks the movement of electrons on the donor; therefore, TY6-2X and TY6-3X may exhibit a better performance.
The MOs of the designed molecules are shown in Table 11. Compared with the original molecule TY6   The -CN group can reduce the molecular LUMO and lead to the decrease in the energy gap. Through the analysis, it can be concluded that the -CN group is the most conducive to the reduction of the molecular energy gap, followed by  11 International Journal of Photoenergy the -CF 3 group, and least influence corresponding to the -F group. The electronegativity of the 3 groups increased gradually, which was the same as that of the LUMO. Therefore, the greater the electronegativity of the inserted EWG is, the higher the energy of the LUMO is. Table 12 As shown, 2/3 of the designed molecules are higher than the original molecules. The maximum value appears in TY6-2C (0.9726), which is the insertion of the -F group in acceptor position 2. The minimum value appears in TY6-1A (0.7379), which is the insertion of the -CF 3 group in acceptor position 4. It is also found that the designed molecules of the insertion group in acceptor position 1 are less than the original molecules, and the performance is poor. However, the designed molecules inserted at positions 2 and 3 of the acceptor are higher than the original molecules, which shows a better performance. This phenomenon is consistent with the previous conclusions that positions 2 and 3 are better than positions 1 and 4. Besides, the magnitude of oscillator strength is directly related to the light harvesting efficiency (LHE). It can be seen from Table 13 that the order of the LHE of the designed molecules is consistent with the oscillator strength (TY6, 0.871). What is important is that a higher LHE leads to a better PCE, and it comes to the following conclusion: Molecules with groups inserted at positions 2 and 3 will have a better performance than those at positions 1 and 4. Table 13 shows some of the chemical reaction parameters of the design molecules. Compared with the original molecule TY6, most of the designed molecules have significantly increased the values of chemical reaction parameters, and the effect is obvious. For EA, the designed molecules are arranged in the following order: TY6-3B = TY6-2B = TY6-3ATY6- The maximum value of EA appears at TY6-3B/TY6-2B/TY6-3A (1.64 eV), which adds a -CN/− CF 3 group on position 3/2 of the acceptor; the minimum value of EA occurs at TY6-1C (1.35 eV), which adds an -F group on position 1 of the acceptor. For ω + , the designed molecules are arranged in the following order: TY6-3A > TY6- The maximum value of ω + appears at TY6-3A (1.77 eV), which adds a -CF 3 group on position 3 of the acceptor; the minimum value of ω + occurs at TY6-1C (1.37 eV), which adds an -F group on position 1 of the acceptor. For ω, the designed molecules are arranged in the following order: The maximum value of ω appears at TY6-3B/TY6-2B (3.31 eV), which adds a -CN group on position 3/2 of the acceptor; the minimum value of ω occurs at TY6-1C (2.83 eV), which adds an -F group on position 1 of the acceptor. As mentioned above, we can get some general conclusions as follows. First, for the three kinds of chemical reactivity parameters, the order of the designed molecules is almost consistent, and the order of EA and ω is exactly the same. Second, the maximum value of the chemical reactivity parameters is found in TY6-3X (position 3 of the acceptor); the corresponding minimum value appears at TY6-1C, and its value is smaller than the original molecular TY6. As a whole, the molecular properties designed at positions 2 and 3 of the acceptor are better than those designed at position 1 and 4 of the acceptor. Finally, all the molecules inserted by the -CN group showed better performance than the original molecule TY6 did. Furthermore, the introduction of the -CN group not only can effectively improve the electronic transmission performance and the ability to receive electrons but also can improve the stability of energy.
The total static first hyperpolarizability of the designed molecules is listed in Table 13. Compared with TY6 (747115), the designed molecules have obvious changes and are arranged in the following order: TY6-3B > TY6- The maximum value of β tot appears at TY6-3B (1068163), which adds a -CN group on position 3 of the acceptor; the minimum value of β tot occurs at TY6-4A (529611), which adds a -CF 3 group on position 4 of the acceptor. Through the above analysis, it can be concluded that, for different positions of the acceptor, the molecules at acceptor positions 2 and 3 are designed to be better than the acceptors designed at positions 1 and 4; for different groups, the designed molecules inserted by the -CN group have the highest β tot , while the dye inserted by the -CF 3 group has the worst performance. Therefore, by inserting the -CN group, the molecules can have more obvious ICT. The lifetime (t) of the designed molecules is listed in Table 13. Compared to the original molecule TY6 (3.91 ns), the majority of the molecules behaved better and are arranged in the following order: TY6- The maximum value of t appears at TY6-1A (4.57 ns), which adds a -CF 3 group on position 3 of the acceptor; the minimum value of t occurs at TY6-2C (3.48 ns), which adds a -F group on position 2 of the acceptor. Through the above analysis, we can get some conclusions: Molecules designed at position 1 of the acceptor have a longer lifetime. The lifetime of the designed molecules inserted at position 2 of the acceptor is lower than that of the original molecule (TY6), which is contrary to the conclusion that positions 2 and 3 are better than positions 1 and 4. However, there is no significant difference between the different groups.

Conclusion
In this study, the DFT and TD-DFT methods were used to calculate the properties of the ground and excited states for MS3, TY3, and TY6. The calculated results show that TY6 has the highest HOMO energy compared to MS3 and TY3, resulting in the smallest energy gap. Smaller energy gaps favor the red-shifted absorption; so, TY6 has the most pronounced red-shifted absorption and the highest molar extinction coefficient, which results in more effective absorption of sunlight. TY6 not only has the largest EA and the lowest IP but also has lower hole/electron reorganization International Journal of Photoenergy energies for higher hole and electron transfer. Through orbital and NBO analysis, the D-A-π-A (TY6) structure enhances the donoring electron ability and acceptor ability of receiving electrons. At the same time, TY6 also has the highest β and the most obvious ICT. Through ESP analysis, we can find that TY6 has better reactivity. TY6 has the lowest h and the highest LHE, ΔG inject , ΔG regen dye , and ω + , which leads to higher J SC . A higher μ normal results in a higher V OC . This is consistent with the experimental results that TY6 has the highest J SC and V OC .
The calculated results for the design molecules show that the overall effect is improved to some extent. For different groups, by introducing the -CN group, the LUMO energy and energy gap can be reduced, leading to an obvious redshifted absorption. Chemical reactivity parameters and β were significantly improved. For the different positions of the molecular acceptors, the positions away from the acceptors (positions 2 and 3) performed better than the other positions did. Especially, the oscillator strength led to a significant increase in LHE. Therefore, it will be of some reference to molecular design.

Data Availability
The three molecules (MS3, TY3, and TY6) used in this article and the experimental data mentioned were obtain from [36,37], and other data were calculated using software such as the Gaussian 09 package. In addition, other data generated or analyzed during this study are available from the first and corresponding authors on reasonable request.

Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this article.