Research Article  Open Access
Ram S. Bhatta, Mesfin Tsige, "Structural Dependence of Electronic Properties in AADAAType Organic Solar Cell Material", International Journal of Photoenergy, vol. 2015, Article ID 708048, 7 pages, 2015. https://doi.org/10.1155/2015/708048
Structural Dependence of Electronic Properties in AADAAType Organic Solar Cell Material
Abstract
Small conjugated molecules (SCMs) are promising candidates for organic photovoltaic (OPV) devices because of their structural simplicity, well control over synthetic reproducibility, and low purification cost. However, industrial development of SCMbased OPV devices requires improving their performance, which in turn relies on the fundamental understanding of structural dependence of electronic properties of SCMs. Herein, we report the structural and electronic properties of the BCNDTS molecule as a model system for acceptoracceptordonoracceptoracceptor (AADAA) type SCMs, using density functional theory (DFT) and timedependent DFT methods. Systematic calculations of twodimensional potential energy surfaces, molecular electrostatic potential surfaces, ground state frontier molecular orbital energies, and the vertical excitation energies are performed. We found that the lowest energy conformation of the BCNDTS molecule is planar. The planar conformation favors the lowest ground state and the excited state energies as well as the strongest oscillator strength. The present results suggest that SCMs containing central dithienosilole cores connected with 2,1,3benzothiadiazole groups have potential to be an efficient electron donor for OPV devices.
1. Introduction
Organic photovoltaic (OPV) devices have gained wide interest in recent years because of the world’s growing demand for clean, renewable, and sustainable energy [1–3]. The promise of OPV devices derives not only from their light weight, flexibility, and environmental friendliness, but also from an expected drastic reduction of the cost associated with the electricity generation [4]. OPV devices contain two active materials: electron acceptors and electron donors. Fullerene derivatives are commonly used as electron acceptors. Electron donors are either conjugated polymers (CPs) or small conjugated molecules (SCMs). Although the current PCE of singlejunction CPbased OPV devices has reached ~9% [5], CPs are suffering from various problems such as lack of synthetic reproducibility, high purification cost, and residual terminal groups of polymers [6]. A potential option to overcome these problems is the replacement of CPs by soluble SCMs [7]. Therefore, research on SCMbased OPV devices has grown very rapidly reaching the current PCE of ~7%, which is comparable to the PCE of the best CPbased OPV devices [8]. However, the PCE of SCMbased OPV devices needs to be further increased for the industrial development of OPV devices.
Various experimental studies [8–14] have been concentrated on the synthesis of SCMs and their OPV performances in recent years. Among SCMs, molecules containing central dithienosilole cores connected with 2,1,3benzothiadiazole groups have been reported as the promising electron donors for OPV devices [8, 9]. In SCMs, typical charge carriers are holes and electrons in porbitals [15]. Electron transfer from SCMs to fullerene derivatives depends on the energy difference between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). Furthermore, HOMO and LUMO energies are affected by the accessible molecular conformations. Although the synthesis and the performance of dithienosilolebased molecule (Figure 1) are studied experimentally, the fundamental factor such as the conformational dependence of electronic properties that affect the performance of these molecules remained unexplored. Thus, a necessary step of understanding the performance of SCMs in OPV application is to investigate the relative variation of conformational energies caused by the torsional defects, HOMOLUMO gaps, and the vertical excitation energies.
In the present study, we perform firstprinciples calculations to investigate the structural and electronic properties of BCNDTS molecule. This work is organized as follows. Section 2 contains computational methods employed in this study. The results and discussion are presented in Section 3. Twodimensional potential energy surface for the BCNDTS molecule is discussed in Section 3.1. Then, in Sections 3.2 and 3.3, we present surfaces of molecular electrostatic potential and HOMOLUMO gap surfaces, respectively. Section 3.4 presents the vertical excitation energies, followed by the main conclusions of the present work in Section 4.
2. Computational Methods
The molecular structure of BCNDTS molecule studied is shown in Figure 1. The two torsional angles ( and ) are defined in terms of the relevant interring C–C–C–C dihedral angles in such a way that = = 0° represents the molecular conformation as shown in Figure 1. Two types of geometry minimizations were performed. The first is the fully relaxed geometry optimization of BCNDTS molecule without any constraint. The second is sets of geometry optimizations at constrained values of and in the range of 0° ≤ ≥ 180° and 0° ≤ ≥ 180°. All geometry optimizations were performed using density functional theory (DFT) [16] combined with Becke’s threeparameter LeeYangParr exchangecorrelations functional (B3LYP) [17] as implemented in the Gaussian 09 program package [18]. The dispersion correction was implemented using Grimme approach [19]. The 631G(d) basis set was used with diffuse functions on all heavy atoms. The dispersioncorrected B3LYP functional was chosen because our previous study [20] on CPs has shown that it yields accurate geometry and torsional potentials. The robustness of B3LYP functional to produce accurate geometries and energies of SCMs and CPs has also been validated in similar studies [21–25].
The HOMOLUMO gap in the ground state was computed as the energy difference between HOMO and LUMO. The vertical excitation energy from the ground state to the first singletexcited state was computed using timedependent density functional theory (TDDFT) [26] method starting from the ground state optimized geometry.
In the context of the present work, some computations of electronic properties were also performed using the coulomb attenuated functional (CAMB3LYP) [27] and longrage corrected functionals (LCBLYP [28] and wB97XD [29]) for the sake of comparison. Table 1 summarizes the computed band gaps at different functionals for fully relaxed BCNDTS molecule. Available experimental values are also included for comparison. Different functionals give different HOMO, LUMO, and vertical excitation energies. The dispersioncorrected B3LYP provides overall good agreement with the experimental values, whereas coulomb attenuated and longrange separated functionals overestimate the values. Therefore, all the results discussed hereafter are computed at the dispersioncorrected B3LYP level unless otherwise indicated.

3. Results and Discussion
3.1. TwoDimensional Potential Energy Surface
In the pristine BCNDTS bulk, chains conformations are controlled by intermolecular packing interactions that favor planar backbone structures. However, in the practical OPV devices, BCNDTS is heterogeneously mixed with fullerene derivatives. It is interesting to explore if there is a propensity of backbone conformational disorder in such environment. Different conformations are connected by the interring torsional motion of the backbone. Hence, computation of potential energy surface as functions of the backbone torsional angles ( and ) allow us to probe the propensities of nonplanar conformations that have relatively close in energy, kT at room temperature.
To obtain the potential energy surface, fortynine grid points were computed at constrained values of and in the ranges 0 ≤ ≤ 180° and 0 ≤ ≤ 180°. Computed energies are summarized in Table 2.

Then, the computed energies were fitted using (1) with onedimensional and coupling terms: In (1), is the constant and the rootmeansquare residual of the final fit was 0.6 meV. The contribution of coupling terms in the fit is negligible indicating that the nearestneighbor torsional coupling in BCNDTS molecule is negligible. Figure 2 shows the twodimensional potential energy surface of the BCNDTS molecule as a function of the two angles and . The contour plot of the same data is provided in the supporting information (see Figure S1 in Supplementary Material available online at http://dx.doi.org/10.1155/2015/708048) for the sake of better visualization.
The twodimensional potential energy surface (Figure 2) has a minimum at . The surface has three local minima at and , and , and and , respectively. At , the surface has a maximum. The maximum is connected with other four local maxima at and , and , and , and and , respectively. These local maxima are connected by other local minima. For example, the local minimum at and is connected with another local minimum at and via the maximum at and . The energies of local minima at and and and are ~21 meV and are easily accessible at room temperature. However, the energy of local minima at and is ~2 kT at room temperature. These results show that the BCNDTS molecule favors more planar conformation compared to CPs as found in our previous study [20].
The form of the potential energy surface (Figure 2) can be understood in terms of the competitions among various interactions such as piconjugation, nonbonded repulsive interaction, and hydrogen bonding. The molecular conformation at and is stabilized by maximum effective piconjugation across the backbone chain and by hydrogen bonding between hydrogen and nitrogen atoms. The conformational energy at and is ~2 kT higher than the one at and because of the absence of hydrogen bonding and the nonbonded repulsion between nitrogen and sulfur atoms. The maximum of the potential energy surface at and is accounted for the breakage of piconjugation across the backbone chain.
3.2. MEP Surfaces
BCNDTS molecule is composed of three units. Dithienosilole is the central unit, which is connected to cyano groups through 2,1,3benzothiadiazole moieties. Cyano cores are strong electronwithdrawing groups. They pull electron density towards them by exhibiting inductive effect when bonded to 2,1,3benzothiadiazole units. Consequently, the electronaccepting 2,1,3benzothiadiazole units also pull electron density from the central dithienosilole unit. As a result, electron densities of cyano units are higher compared to dithienosilole and 2,1,3benzothiadiazole units. These electron densities can be visualized by constructing threedimensional molecular electrostatic potential (MEP) surfaces.
Figure 3(a) shows the MEP surfaces of fully relaxed BCNDTS molecule. The region of the MEP surface with blue color represents a lower electron density and has a positive charge. Red color on the MEP surface corresponds to a high electron density and has a negative charge. In the absence of torsional defects, the electron density is uniformly distributed across the BCNDTS molecule with sharp blue color on dithienosilole region and sharp red on cyano region. However, the electron density distribution is affected in the presence of the torsional defect.
(a)
(b)
(c)
(d)
At and (Figure 3(b)), electron density on one 2,1,3benzothiadiazole region is distinct with more extended red color. Similarly, at and (Figure 3(c)), the MEP surface of BCNDTS molecule represents the extension of red color on both 2,1,3benzothiadiazole units compared to the MEP surface without the torsional defect. At and (Figure 3(d)), MEP surface resembles closely the surface at and . Hence, the extension of red color in the 2,1,3benzothiadiazole units increases with the breakage of conjugation across the backbone chain.
3.3. HOMOLUMO Gap
Torsional defect affects not only the energy (Figure 2) and electron density (Figure 3), but also the HOMOLUMO gap of BCNDTS molecule. As discussed in our previous studies [20, 22] of poly(3alkylthiophene) and thieno[3,4b]thiophenealtbenzodithiophene copolymer, a usual approach for estimating torsional dependence of the HOMOLUMO gap is to compute the band gap variation when and are varied in the range of 0° to 180°. The HOMOLUMO gap is minimum for planar conformation. It increases with the increase in the torsional angle and reaches the maximum when the two planarhalves are orthogonal to each other.
Figure 4 shows the contour plot of HOMOLUMO gap for BCNDTS molecule as a function of the two angles, and . The twodimensional HOMOLUMO gap surface is provided in supporting information (Figure S2). The plot was generated from the computed fortynine grid points using the fit function as described in Section 3.1. The rootmeansquare residual of the final fit was 0.02 eV. Like the potential energy surface, only onedimensional terms are significant in the fit of HOMOLUMO gap.
Figure 4 has a minimum at and a maximum at . The maximum is connected with four local maxima at and , and , and , and and , respectively. The maximum is also connected with four minima at , and , and , and and , respectively. The HOMOLUMO gap at and and at and are 0.02 eV and 0.04 eV higher, respectively, compared to the HOMOLUMO gap at and .
The torsional dependence of band gap can be understood in terms of HOMO and LUMO energies. At or , the planar geometry favors the effective overlap of porbitals resulting in the extended conjugation across the HOMO and LUMO (Figures 5(a) and 5(d)). The overlapping of porbitals is decreased as the values of and are increased. When and/or approach , the extended conjugation is broken as shown in Figures 5(b) and 5(c). Consequently, the energies of HOMO and LUMO get more separated and the band gap of BCNDTS molecule is increased as depicted in Figure 4.
(a)
(b)
(c)
(d)
3.4. Excitation Energies
The TDDFT calculations were performed on three representative conformations with and , and , and and to investigate the effect of torsional defect on excitation energies and the oscillator strength. Figure 6 shows the computed UV/VIS spectra for the representative conformations with their oscillator strengths. The planar conformation () has the vertical excitation energy of ~2 eV from the ground (S_{0}) state to the first excited singlet (S_{1}) state with the highest oscillator strength of ~1.2. At and , the spectrum is blue shifted by ~ 0.3 eV. At and , the spectrum is further blue shifted with the reduced oscillator strength of ~0.8.
The transition from S_{0} to S_{1} represents the promotion of an electron from HOMO to LUMO. Stronger peak in Figure 6 for the planar conformation represents the electronic transition from HOMO to LUMO and a weaker peak represents the electronic transition from HOMO to LUMO+2 (with the oscillator strength of ~0.3). At and , BCNDTS molecule shows two peaks with the oscillator strength of ~0.5. These peaks represent the electronic transition from HOMO to LUMO+1 and HOMO1 to LUMO+1, respectively. At the maximum torsional defect (), BCNDTS molecule shows a peak with the oscillator strength of ~0.7, which represents the transition from HOMO2 to LUMO or HOMO to LUMO+2. There are two important effects of torsional defects: (i) the oscillator strength decreases with the increase in torsional defects and (ii) the spectrum is blue shifted with the increase in torsional defects. These results indicate that the torsional defects present in BCNDTS molecule decrease the OPV performance, thereby reducing the strength of ππ^{*} transition.
4. Conclusions
DFT and TDDFT calculations were performed to investigate the structural and electronic properties of BCNDTS molecule. Twodimensional potential energy surface as function of interring torsional angles and energies of frontier molecular orbitals were computed using dispersion corrected B3LYP functional. The ground state electronic properties calculated at the dispersion corrected B3LYP were found to be in good agreement with the available experimental values. The vertical excitation energies were calculated starting from the ground state optimized geometries at the dispersion corrected B3LYP level.
In the present investigation, the lowest energy structure of BCNDTS molecule was found to have planar conformation. The planar conformation was found to favor the lowest HOMOLUMO gap and excitation energies but the highest oscillator strength. The oscillator strength was found to be decreased and the absorption spectrum was found to be blue shifted for higher energy conformations caused by the torsional defects. The present study suggests that small conjugated molecules containing central dithienosilole cores connected with 2,1,3benzothiadiazole groups, which favor the lowest energy planar conformation, have potential to be an efficient electron donor for OPV devices.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors are grateful to the National Science Foundation (Grant no. DMR0847580) for financial support. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation Grant no. ACI1053575. The authors thank additional support from the University of Akron.
Supplementary Materials
Figure S1 shows the contour plot of potential energy as functions of torsional angles a and β. The contour plot has a minimum at a = β = 0° and a maximum at a = β = 90°. The maximum is connected with other three local minima via local maxima. Refer to the text for detailed discussion.
Figure S2 shows the twodimensional HOMOLUMO surface as functions of torsional angles a and β. The HOMOLUMO surface has a minimum at a = β = 0° and a maximum at a = β = 90°. There are three local minima at a = 0° and β = 180°, a = 180° and β = 0°, and a = 180° and β = 180°, respectively. These local maxima are connected by three local minima. Refer to the text for detailed discussion.
References
 C. J. Brabec, N. S. Sariciftci, and J. C. Hummelen, “Plastic solar cells,” Advanced Functional Materials, vol. 11, no. 1, pp. 15–26, 2001. View at: Publisher Site  Google Scholar
 H. Hoppe and N. S. Sariciftci, “Organic solar cells: an overview,” Journal of Materials Research, vol. 19, no. 7, pp. 1924–1945, 2004. View at: Publisher Site  Google Scholar
 N. Serdar Sariciftci, “Plastic photovoltaic devices,” Materials Today, vol. 7, no. 8, pp. 36–40, 2004. View at: Publisher Site  Google Scholar
 S. B. Darling and F. You, “The case for organic photovoltaics,” RSC Advances, vol. 3, no. 39, pp. 17633–17648, 2013. View at: Publisher Site  Google Scholar
 Z. He, C. Zhong, S. Su, M. Xu, H. Wu, and Y. Cao, “Enhanced powerconversion efficiency in polymer solar cells using an inverted device structure,” Nature Photonics, vol. 6, no. 9, pp. 591–595, 2012. View at: Publisher Site  Google Scholar
 J. Roncali, P. Leriche, and P. Blanchard, “Molecular materials for organic photovoltaics: small is beautiful,” Advanced Materials, vol. 26, no. 23, pp. 3821–3838, 2014. View at: Publisher Site  Google Scholar
 J. Roncali, P. Frère, P. Blanchard et al., “Molecular and supramolecular engineering of $\pi $conjugated systems for photovoltaic conversion,” Thin Solid Films, vol. 511512, pp. 567–575, 2006. View at: Publisher Site  Google Scholar
 Y. Sun, G. C. Welch, W. L. Leong, C. J. Takacs, G. C. Bazan, and A. J. Heeger, “Solutionprocessed smallmolecule solar cells with 6.7% efficiency,” Nature Materials, vol. 11, no. 1, pp. 44–48, 2012. View at: Publisher Site  Google Scholar
 L.Y. Lin, C.W. Lu, W.C. Huang, Y.H. Chen, H.W. Lin, and K.T. Wong, “New AADAAtype electron donors for small molecule organic solar cells,” Organic Letters, vol. 13, no. 18, pp. 4962–4965, 2011. View at: Publisher Site  Google Scholar
 S. Roquet, A. Cravino, P. Leriche, O. Alévêque, P. Frère, and J. Roncali, “Triphenylaminethienylenevinylene hybrid systems with internal charge transfer as donor materials for heterojunction solar cells,” Journal of the American Chemical Society, vol. 128, no. 10, pp. 3459–3466, 2006. View at: Publisher Site  Google Scholar
 J. Zhou, X. Wan, Y. Liu et al., “Small molecules based on benzo[1,2b:4,5b′]dithiophene unit for highperformance solutionprocessed organic solar cells,” Journal of the American Chemical Society, vol. 134, no. 39, pp. 16345–16351, 2012. View at: Publisher Site  Google Scholar
 H. Lu, C. Lu, Y. Lee et al., “New molecular donors with dithienopyrrole as the electrondonating group for efficient smallmolecule organic solar cells,” Chemistry of Materials, vol. 26, no. 15, pp. 4361–4367, 2014. View at: Publisher Site  Google Scholar
 H.C. Ting, Y.H. Chen, L.Y. Lin et al., “Benzochalcogenodiazolebased donoracceptoracceptor molecular donors for organic solar cells,” ChemSusChem, vol. 7, no. 2, pp. 457–465, 2014. View at: Publisher Site  Google Scholar
 D. Demeter, S. Mohamed, A. Diac, I. Grosu, and J. Roncali, “Small molecular donors for organic solar cells obtained by simple and clean synthesis,” ChemSusChem, vol. 7, no. 4, pp. 1046–1050, 2014. View at: Publisher Site  Google Scholar
 A. Mishra and P. Bäuerle, “Small molecule organic semiconductors on the move: promises for future solar energy technology,” Angewandte Chemie  International Edition, vol. 51, no. 9, pp. 2020–2067, 2012. View at: Publisher Site  Google Scholar
 P. Hohenberg and W. Kohn, “Inhomogeneous electron gas,” Physical Review, vol. 136, pp. B864–B871, 1964. View at: Google Scholar  MathSciNet
 A. D. Becke, “Densityfunctional thermochemistry. III. The role of exact exchange,” The Journal of Chemical Physics, vol. 98, no. 7, pp. 5648–5652, 1993. View at: Publisher Site  Google Scholar
 M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian 09, Revision D.01, Gaussian Inc., Wallingford, Conn, USA, 2009.
 S. Grimme, “Semiempirical GGAtype density functional constructed with a longrange dispersion correction,” Journal of Computational Chemistry, vol. 27, no. 15, pp. 1787–1799, 2006. View at: Publisher Site  Google Scholar
 R. S. Bhatta, D. S. Perry, and M. Tsige, “Nanostructures and electronic properties of a highefficiency electrondonating polymer,” The Journal of Physical Chemistry A, vol. 117, no. 47, pp. 12628–12634, 2013. View at: Publisher Site  Google Scholar
 R. S. Bhatta and D. S. Perry, “Correlated backbone torsional potentials in poly(3methylthiophene),” Computational and Theoretical Chemistry, vol. 1008, pp. 90–95, 2013. View at: Publisher Site  Google Scholar
 R. S. Bhatta, M. Tsige, and D. S. Perry, “Torsionallyinduced blue shift of the band gap in poly(3hexylthiophene),” Journal of Computational and Theoretical Nanoscience, vol. 11, no. 10, pp. 2157–2164, 2014. View at: Publisher Site  Google Scholar
 M. D'Auria, L. Emanuele, and R. Racioppi, “Regio and stereoselectivity in the PaternòBüchi reaction on furan derivatives,” International Journal of Photoenergy, vol. 2006, Article ID 84645, 11 pages, 2006. View at: Publisher Site  Google Scholar
 T. Chokbunpiam, P. Thamyongkit, O. Saengsawang, and S. Hannongbua, “Molecular structure and electronic properties of porphyrinthiopheneperylene using quantum chemical calculation,” International Journal of Photoenergy, vol. 2010, Article ID 492313, 8 pages, 2010. View at: Publisher Site  Google Scholar
 T. Yakhanthip, N. Kungwan, J. Jitonnom, P. Anuragudom, S. Jungsuttiwong, and S. Hannongbua, “Theoretical investigation on the electronic and optical properties of poly(fluorenevinylene) derivatives as lightemitting materials,” International Journal of Photoenergy, vol. 2011, Article ID 570103, 9 pages, 2011. View at: Publisher Site  Google Scholar
 E. Runge and E. K. U. Gross, “Densityfunctional theory for timedependent systems,” Physical Review Letters, vol. 52, no. 12, pp. 997–1000, 1984. View at: Publisher Site  Google Scholar
 T. Yanai, D. P. Tew, and N. C. Handy, “A new hybrid exchangecorrelation functional using the Coulombattenuating method (CAMB3LYP),” Chemical Physics Letters, vol. 393, no. 1–3, pp. 51–57, 2004. View at: Publisher Site  Google Scholar
 H. Iikura, T. Tsuneda, T. Yanai, and K. Hirao, “A longrange correction scheme for generalizedgradientapproximation exchange functionals,” Journal of Chemical Physics, vol. 115, no. 8, pp. 3540–3544, 2001. View at: Publisher Site  Google Scholar
 J.D. Chai and M. HeadGordon, “Longrange corrected hybrid density functionals with damped atomatom dispersion corrections,” Physical Chemistry Chemical Physics, vol. 10, no. 44, pp. 6615–6620, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Ram S. Bhatta and Mesfin Tsige. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.