About this Journal Submit a Manuscript Table of Contents
The Scientific World Journal
Volume 2013 (2013), Article ID 890215, 8 pages
http://dx.doi.org/10.1155/2013/890215
Research Article

Effect of Electronic Acceptor Segments on Photophysical Properties of Low-Band-Gap Ambipolar Polymers

1College of Science, Northeast Forestry University, Harbin 150040, China
2School of Physics and Optoelectronic Technology, Dalian University of Technology, Dalian 116024, China
3College of Science, Jiamusi University, Jiamusi 154001, China
4Department of Physics, Liaoning University, Shenyang 110036, China

Received 6 November 2012; Accepted 17 December 2012

Academic Editors: J. Karwowski, A. Kochur, and B. Zygelman

Copyright © 2013 Yuanzuo Li et al. 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.

Abstract

Stimulated by a recent experimental report, charge transfer and photophysical properties of donor-acceptor ambipolar polymer were studied with the quantum chemistry calculation and the developed 3D charge difference density method. The effects of electronic acceptor strength on the structure, energy levels, electron density distribution, ionization potentials, and electron affinities were also obtained to estimate the transporting ability of hole and electron. With the developed 3D charge difference density, one visualizes the charge transfer process, distinguishes the role of molecular units, and finds the relationship between the role of DPP and excitation energy for the three polymers during photo-excitation.

1. Introduction

Ambipolar devices are very attractive organic semiconductors due to their high chemical stability and uncommon versatility; particularly, ambipolar copolymer films have been applied as active materials in the field of organic optoelectronic materials [16]. For example, ambipolar organic field-effect transistors, which are capable of both p- and n-channel operations by changing the polarity of the gate voltage, are gaining attention as an alternative approach to mimicking complementary metal-oxide semiconductor (CMOS) digital integrated circuits for achieving high-performance and cost-effective circuits in organic electronics. On one hand, the semiconductors suitable as single component channel materials should facilitate the formation of an exciton through the cation and anion radicals, and they show both stable hole and electron transporting characteristics; on the other hand, the semiconductors are required to have a well match between ambipolar structures and metal electrodes to balance the charge injection barriers between the relative positions of HOMO and LUMO energy levels and the work function of electrodes.

The 1, 4-diketo-3, 6-diarylpyrrolo pyrroles diketopyrrolopyrroles (DPPs) are chromophoric systems that combine in a rigid planar structural frame. Generally, the DPP unit possesses the strong electron-withdrawing capability. When combining with appropriate building blocks (such as thiophene, phenylene, and benzofuran building blocks), the certain DPP-containing polymers show ambipolar characteristic. Meanwhile, this structural adjustment should result in tunable band gap, energy levels, and molecular packing. Therefore, DPP-containing polymers show wide application prospect in the field of organic photovoltaics [7, 8], electroluminescence device [9, 10], organic field-effect transistors [11, 12], and logic circuits [13, 14]. Recently, electronic, physical, and transistor properties of a family of donor-acceptor polymers (which consist of DPP coupled with neutral benzene (B), the weakly accepting benzothiadiazole (BT) and the strongly accepting benzobisthiadiazole (BBT)) have been experimented [15]. The report found that the rational design can realize the conversion from unipolar to ambipolar donor-acceptor polymer. Although experimental study made tremendous progress, it is important to understand the relationship between chemical structures and the optical and electronic properties of the organic ambipolar layer from the molecular level. Theoretical calculations not only provide some parameters (such as electron transition, ionization potential (IP), and electron affinity (EA)) affecting the optical response and charge transport ability, but also reveal the microscopic mechanism behind the experiment and give theoretical basis to rational design of new functional materials. In this paper, we attempt to study the electronic structure and photo-physical nature from the viewpoint of theory, to deeply insight into photo-induced charge transfer process with the developed 3D cube representations [1619]. This method can visually distinguish the role of molecular units during the photo-excitation. The article is organized as follows: firstly, the ground and excited state geometries are optimized with density functional theory. Secondly, molecular orbital shape and energy, IP and EA are estimated. Lastly, charge transfer and excited states properties are studied, and role of molecular units upon excitation is distinguished by the visualized analysis.

2. Methods

The ground state geometries of PBDPP, PBTDPP, and PBBTDPP (see Figure 1) were optimized with density functional theory (DFT) [20], B3LYP functional [2123], and 6-31G (d) basis set. The alkyl groups on five members were replaced with methoxyl groups due to the fact that the alkyl groups do not significantly affect the equilibrium and optical property of fluorene-based polymer [24]. The excited-state geometries were optimized by time-dependent DFT (TD-DFT) [25], by using the 6-31G (d) basis set. The electronic transition of them for absorption and fluorescence was calculated at both ground-state and excited-state optimized geometries by means of TD-DFT, Cam-B3LYP functional [26] with the same basis set. The photoinduced intramolecular charge transfer process was visualized with the 3D cube representations. 3D charge difference density indicates that the electronic redistribution involving the whole structure takes place upon excitation [16, 17, 27, 28]. The charge difference density is defined as where is the uth eigenvector of the single configuration interaction (CI) Hamiltonian in the basis of the occupied Hartree-Fock molecular orbital and the unoccupied orbital [16, 27], and in this equation the first and the second terms stand for hole and electron, respectively. All the quantum chemical calculations were done with Gaussian 09 software [29]. The optical absorption spectra of them were done with GaussSum 2.1 [30].

890215.fig.001
Figure 1: Molecular structures of PBDPP, PBTDPP, and PBBTDPP.

3. Results and Discussion

3.1. Ground Optimization

The ground state structures of three molecular structures (PBDPP, PBTDPP, and PBBTDPP) have been optimized, and selected optimized bond length and dihedral angles of them were listed in Table 1. Comparing the dihedral angles of three polymers, it found that PBDPP, PBTDPP, and PBBTDPP displayed an obvious twisted configuration along the main molecular skeleton. For example, thiophene-diketopyrrolopyrrole (T-DPP) of PBDPP, PBTDPP, and PBBTDPP are 13.76026, 11.78848, and 10.44882 degree, respectively. The dihedral angles of thiophene-benzene (T-B) in PBDPP, thiophene-benzothiadiazole (T-BT) in PBTDPP, and thiophene-benzobisthiadiazole (T-BBT) are 21.51214, −1.82834, and −0.27066 degree, respectively. Although the PBBTDPP also has a certain twisted structure, it exhibits the well conplanar structure than the other two. The order of coplanarity is increased from PBDPP, PBTDPP to PBBTDPP. The date of intermolecular bond length demonstrated the conjugated degree of PBDPP is lower than that of PBTDPP and PBBTDPP. Therefore, introduction of BT and BBT improves the molecular coplanarity, which will lead to enhancement of electronic delocalization and improvement of their optical response.

tab1
Table 1: Selected optimized bond length and dihedral angles of three polymers in the ground and excited state.

It will be important to analyse the energy levels and the electron density distribution of molecular orbitals (MOs), since the change of electron density of MOs upon photo-excitation provides a reasonable qualitative indication of the excitation properties [31]. The energy levels of three polymers are indicated in Figure 2, and contour plot surfaces of HOMO and LUMO of them are inserted in Figure 2. As shown, for PBDPP electron densities of HOMO are mainly located in the two thiophens of first monomer, the DPP and two thiophens of second monomer. As compared with HOMO, the LUMO of PBDPP also resides the right two thiophens and the DPP and two thiophens, and a few of electron density is in the backbone chain of first monomer. For PBTDPP, electron density of the HOMO mainly resides in the whole dimer and electron density of the LUMO in the first monomer. Upon photo-excitation, there is a decreasing trend of electron density in the second monomer, and an increasing trend of electron density occurs in the BT unit of first monomer. The introduction of BT not only induces the change of electron density distribution, but also destabilizes the energy of LUMO and stabilizes the energy of HOMO, so PBTDPP has a lower energy gap than PBDPP. Among the three molecules, PBBTDPP has narrowest energy gap about 0.89 eV closed to the experiment value of 0.65 eV [15]. In addition, from the electron density distribution, it found that electron density of HOMO distributes over the whole backbone, while electron density of LUMO is mainly located in the BBT unit of the first monomer. So the character of MOs for PBBTDPP displays an obvious electron moving from the second monomer to the BBT unit.

890215.fig.002
Figure 2: The energy levels and the electron density distribution of molecular orbitals.

Ionization potentials (IPs) and electron affinities (EAs) as the criterion of evaluating the injection ability of the hole and electron have been obtained by calculating adiabatic potential on the basis of neutral and charged optimization. For the three polymers, the value of ionization potentials for PBDPP, PBTDPP, and PBBTDPP is 5.259, 5.214, and 5.032 eV, respectively, that is, it decreases in this order ; this means that for PBBTDPP with the lowest IP the hole injection is more easily than the others. Comparing the electron affinities of PBDPP, PBTDPP, and PBBTDPP, one found that they have an increasing trend, that is, EA value is 2.023 (PBDPP), 2.367 (PBTDPP), and 3.103 eV (PBBTDPP), respectively; so for the PBBTDPP the electron injection from the cathode to the electron transporting layer is more easily than the others.

3.2. Absorption Spectra and Charge Transfer Process

Simulated absorption spectra of PBDPP, PBTDPP, and PBBTDPP are shown in Figure 3, and transition energies and the oscillator strengths are listed in Table 2. It is clearly seen from Figures 3(a) and 3(b) that the absorption peaks of PBDPP and PBTDPP are mainly on the range of visible and ultraviolet light, and the absorption peak is red shifted ~60 nm when the B is replaced by BT. By comparing Figures 3(a) and 3(c), one can see that shape of absorption spectra of PBBTDPP does not change, but the lowest strong absorption peak takes place further significantly red-shifted to the infrared region (to ~1084 nm). For the three polymers, the strongest absorption all corresponding to the S1 excited state that is composed of electron transition HOMO→LUMO (H→L). As shown in Table 2, for PBDPP the second absorption peak is at 375.63 nm, which is close to the second absorption peak of PBTDPP; but for PBBTDPP, the second absorption peak make red-shifted to be 522.83 nm.

tab2
Table 2: Calculated transition energies (eV, nm), CI coefficients, and oscillator strengths () for PBDPP, PBTDPP, and PBBTDPP.
fig3
Figure 3: Absorption spectra of PBDPP, PBTDPP, and PBBTDPP.

In order to understand the nature of absorption spectra, we study the change of charge density during excitation for the calculated ten excited states, and photoinduced intramolecular charge transfer in PBDPP, PBTDPP, and PBBTDPP on electronic transitions can be seen from Figure 4. For the S1 of PBDPP, electron and hole mainly in turn cover the two thiophens and the DPP and two thiophens of second monomer, and hence this state is characterized as * electron transition. Simultaneously, we further checked the other states (see supporting material available online at http://dx.doi.org/10.1155/2013/890215). For the states (S2, S3, and S4) electron and hole are in turn distributed in the one monomer or the two monomers with the character of * electron transition.

890215.fig.004
Figure 4: Charge different density of PBDPP, PBTDPP, and PBBTDPP for absorption (where red and green stand for electron and hole, resp.).

While for S5–S10 states closed to the different excitation energies, they have different charge transfer characters. For example, for S6 state electron move to DPP units and hole in the thiophens and benzenes, so in this state DPP serves as electron acceptor; for the S7–S9 states, DPP serves as electron donor owing to the more electron moving from it to the adjacent thiophens and benzenes. It is worth noting that for the S10 state, it is an obvious intramolecular charge transfer excited state, where electron moves from the DPP of right monomer to the DPP of left monomer, and DPP simultaneously serves as electron donor and acceptor during excitation. So it can be found that for PBDPP (containing neutral benzene (B)) the providing and accepting electronic rules of DPP have the characteristics of wavelength dependence.

For PBTDPP, charge different density (see Figure 4) shows that photo-induced red electron is transferred to the BT subunit (where BT serves as electron acceptor), so this state is an intramolecular charge transfer state (ICT). While for the second absorption peak (S7), excitation takes place on the whole segments of PBTDPP. It is found in Figure 4 that red electron moves to two BT units and green hole in DPP units, so BT unit and DPP unit act as the electronic acceptor and donor group, respectively.

For the PBBTDPP, the S1 state is an ICT state, where electron is transferred from DPP to BBT; therefore, for this strong absorption peak DPP serves as electronic donor and BBT as acceptor. Moreover, there are three excited states with closing oscillator strength, that is, S5 state (564 nm, ), S6 state (522.83 nm, ), and S7 state (497 nm, ), but they have different charge transfer characters. As shown in Figure 4, for S5 state red electron moves to the BBT unit and neighboring thiophen units of left monomer, and green hole resides in the two DPP units and neighboring thiophen units of two monomers, so this state is also ICT state. Exaction of S6 state also occurs on the leading backbone, where more electrons have arisen in the two BBT units. Different with the above-excited states, S7 state comes from the excitation of the single monomer, and electron moves from right thiophen unit of left monomer to the left thiophen unit. Then one can conclude that the second absorption band of PBBTDPP, including the S5, S6, and S7 states, strongly depends on the different ICT modes. Furthermore, for the above three states, BBT units all serve as electron acceptor displaying the strong ability of accepting electron.

We therefore studied in detail the charge transfer of three polymers and visualized the function of every unit in photo-induction. When comparing the ICT states of three polymers, it is found that for the different excited states corresponding to different excitations, DPP, of PBDPP, PBTDPP, and PBBTDPP exhibits different abilities of electronic donor and acceptor gain or loss. Along with the increasing strength of donor groups (from the neutral benzene (B) to the strongly accepting benzobisthiadiazole (BBT)), the rules of electronic donor or acceptor of DPP for three polymers have a change from the initial wavelength dependence to the final accepting groups dependence. Especially for PBBTDPP, DPP nearly acts as the donor segment providing electron and their spectral region cover from 420 nm to 1084 nm.

3.3. Excited State Optimization and Fluorescence Properties

The geometries parameters of excited state optimization by using TDDFT method were listed in Table 1. It can be seen from Table 1 that the bond lengths of adjacent units for the excited state geometries have been mildly shorten about 0.02 angstrom (A) compared with the corresponding bond lengths of ground state geometries. However, the large changes calculated with TDDFT method have taken place in the dihedral angles of adjacent units. For example, for ground state the dihedral angle of T-DPP is 18.96641 (for PBDPP), 19.37309 (for PBTDPP), and 18.69428 (for PBBTDPP); and for excited state it is 18.09190 (for PBDPP), 18.40059 (for PBTDPP), and 17.93294 (for PBBTDPP), respectively. For T-B (for PBDPP), T-BT (for PBTDPP), and T-BBT (for PBBTDPP), the dihedral angles are changed from 21.51214, −1.82834, and −0.27066 to 1.39836, −0.47006, and −0.08738. Therefore, excited state geometries of all oligomers are more coplanar in comparison with ground-state geometries.

The fluorescence energies of all oligomers were calculated on the basis of excited state optimization, as shown in Table 3, and fluorescence lifetime was calculated with the Einstein transition probabilities [32] as follows: where, , and stand for the velocity of light, fluorescence energy, and oscillator strength, respectively. As shown in Table 3, fluorescence peaks of PBDPP and PBTDPP are 660.93 nm and 724.45 nm. Compared with the two formers, fluorescence peaks of PBBTDPP make red-shifted into the infrared region, and the oscillator strength of PBBTDPP is smaller than the others, which result in the longer radiative lifetime according to (2). Therefore, introduction of PBBTDPP not only widens the absorption spectra range but also improves the efficiency of emission.

tab3
Table 3: Fluorescence energies (eV, nm), main configurations, and radiative lifetime of PBDPP, PBTDPP, and PBBTDPP.

4. Conclusion

The stable ground structures of PBDPP, PBTDPP, and PBBTDPP have been optimized with density functional theory, and the energy levels, the electron density distribution, ionization potentials, and electron affinities have also been estimated in current calculations. The results show that the degree of coplanarity is increased from PBDPP, PBTDPP to PBBTDPP; for PBBTDPP it has the lowest IP and highest EA values, indicating that the well electron and hole transporting abilities among three oligomers. The absorption spectra of PBDPP and PBTDPP are mainly on the range of visible and ultraviolet light, while absorption spectra of PBBTDPP along with the introduction of the BBT unit make red-shifted to the infrared region. 3D charge difference density was used to analyze the calculated states, which indicated that the S1 state of PBBTDPP is an intramolecular charge transfer state (ICT), and change of electron density has obvious difference with that of PBDPP and PBTDPP. For the calculated excited states of three polymers, DPP of PBDPP and PBTDPP takes the role of electron donor or acceptor, and the certain role has wavelength dependence; but for PBBTDPP with the BBT unit, the DPP takes the role of electronic donor for most excited states, and PBBTDPP does not display the dependence of excitation energy.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 11274149), Heilongjiang Provincial Youth Science Foundation (Grant no. QC2011C054), the Natural Science Foundation of Liaoning Province (Grant no. 20111035), Liaoning Provincial Department of Education Project (Grant no. L20111003), and the Fundamental Research Funds for the Central Universities (Grant no. DL12BB19).

References

  1. A. Dodabalapur, H. E. Katz, L. Torsi, and R. C. Haddon, “Organic heterostructure field-effect transistors,” Science, vol. 269, no. 5230, pp. 1560–1562, 1995. View at Scopus
  2. H. Yan, Z. H. Chen, Y. Zheng et al., “A high-mobility electron-transporting polymer for printed transistors,” Nature, vol. 457, no. 7230, pp. 679–686, 2009. View at Publisher · View at Google Scholar · View at Scopus
  3. R. Martel, V. Derycke, C. Lavoie et al., “Ambipolar electrical transport in semiconducting single-wall carbon nanotubes,” Physical Review Letters, vol. 87, no. 25, Article ID 256805, 4 pages, 2001.
  4. K. S. Novoselov, A. K. Geim, S. V. Morozov et al., “Electric field in atomically thin carbon films,” Science, vol. 306, no. 5696, pp. 666–669, 2004. View at Publisher · View at Google Scholar · View at Scopus
  5. J. C. Bijleveld, R. A. M. Verstrijden, M. M. Wienk, and R. A. J. Janssen, “Copolymers of diketopyrrolopyrrole and thienothiophene for photovoltaic cells,” Journal of Materials Chemistry, vol. 21, no. 25, pp. 9224–9231, 2011. View at Publisher · View at Google Scholar · View at Scopus
  6. F. Algi and A. Cihaner, “An electroactive polymeric material and its voltammetric response towards alkali metal cations in neat water,” Tetrahedron Letters, vol. 49, no. 21, pp. 3530–3533, 2008. View at Publisher · View at Google Scholar · View at Scopus
  7. D. Sahu, H. Padhy, D. Patra et al., “Synthesis and applications of novel acceptor-donor-acceptor organic dyes with dithienopyrrole- and fluorene-cores for dye-sensitized solar cells,” Tetrahedron, vol. 67, no. 2, pp. 303–311, 2011. View at Publisher · View at Google Scholar · View at Scopus
  8. M. C. Chen, D. J. Liaw, Y. C. Huang, H. Y. Wu, and Y. Tai, “Improving the efficiency of organic solar cell with a novel ambipolar polymer to form ternary cascade structure,” Solar Energy Materials and Solar Cells, vol. 95, no. 9, pp. 2621–2627, 2011. View at Publisher · View at Google Scholar · View at Scopus
  9. F. Algi and A. Cihaner, “An ambipolar low band gap material based on BODIPY and EDOT,” Organic Electronics, vol. 10, no. 3, pp. 453–458, 2009. View at Publisher · View at Google Scholar · View at Scopus
  10. M. Karakus, A. Balan, D. Baran, L. Toppare, and A. Cirpan, “Electrochemical and optical properties of solution processable benzotriazole and benzothiadiazole containing copolymers,” Synthetic Metals, vol. 162, no. 1-2, pp. 79–84, 2012.
  11. T. D. Anthopoulos, C. Tanase, S. Setayesh et al., “Ambipolar organic field-effect transistors based on a solution-processed methanofullerene,” Advanced Materials, vol. 16, no. 23-24, pp. 2174–2179, 2004. View at Publisher · View at Google Scholar · View at Scopus
  12. Y. Kunugi, K. Takimiya, N. Negishi, T. Otsubo, and Y. Aso, “An ambipolar organic field-effect transistor using oligothiophene incorporated with two [60]fullerenes,” Journal of Materials Chemistry, vol. 14, no. 19, pp. 2840–2841, 2004. View at Publisher · View at Google Scholar · View at Scopus
  13. J. C. Bijleveld, A. P. Zoombelt, S. G. J. Mathijssen et al., “Poly(diketopyrrolopyrroleĀterthiophene) for ambipolar logic and photovoltaics,” Journal of the American Chemical Society, vol. 131, no. 46, pp. 16616–16617, 2009. View at Publisher · View at Google Scholar · View at Scopus
  14. F. S. Kim, E. Ahmed, S. Subramaniyan, and S. A. Jenekhe, “Air-stable ambipolar field-effect transistors and complementary logic circuits from solution-processed n/p polymer heterojunctions,” ACS Applied Materials and Interfaces, vol. 2, no. 11, pp. 2974–2977, 2010. View at Publisher · View at Google Scholar · View at Scopus
  15. J. D. Yuen, J. Fan, J. Seifter et al., “High performance weak donor-acceptor polymers in thin film transistors: effect of the acceptor on electronic properties, ambipolar conductivity, mobility, and thermal stability,” Journal of the American Chemical Society, vol. 133, no. 51, pp. 20799–20807, 2011.
  16. M. T. Sun, “Control of structure and photophysical properties by protonation and subsequent intramolecular hydrogen bonding,” Journal of Chemical Physics, vol. 124, no. 5, Article ID 054903, 6 pages, 2006. View at Publisher · View at Google Scholar · View at Scopus
  17. Y. Z. Li, T. Pullerits, M. Y. Zhao, and M. T. Sun, “Theoretical characterization of the PC60BM: PDDTT mode for an organic solar cell,” The Journal of Physical Chemistry C, vol. 115, no. 44, pp. 21865–21873, 2011.
  18. M. Sun, T. Pullerits, P. Kjellberg, W. J. D. Beenken, and K. L. Han, “Control of emission by intermolecular fluorescence resonance energy transfer and intermolecular charge transfer,” Journal of Physical Chemistry A, vol. 110, no. 19, pp. 6324–6328, 2006. View at Publisher · View at Google Scholar · View at Scopus
  19. M. T. Sun and H. X. Xu, “A novel application of plasmonics: plasmon-driven surface-catalyzed reactions,” Small, vol. 8, no. 18, pp. 2777–2786, 2012. View at Publisher · View at Google Scholar
  20. J. M. R. Dreizler and E. K. U. Gross, Density Functional Theory, Springer, Heidelberg, Germany, 1990.
  21. A. D. Becke and R. G. Parr, “Density-functional exchange-energy approximation with correct asymptotic behavior,” Physical Review A, vol. 38, no. 6, pp. 3098–3100, 1988. View at Publisher · View at Google Scholar
  22. A. D. Becke, “Density-functional thermochemistry. III. The role of exact exchange,” The Journal of Chemical Physics, vol. 98, no. 7, pp. 5648–5652, 1993.
  23. C. Lee and W. Yang, “Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density,” vol. 37, no. 2, pp. 785–789, 1988. View at Publisher · View at Google Scholar
  24. Y. Z. Li, F. C. Ma, B. Dong, J. Li, and M. D. Chen, “Theoretical study of charge transfer mechanism in fullerene-phenylphenothiazine compound: a real-space analysis,” Dyes Pigments, vol. 92, no. 3, pp. 1344–1350, 2012.
  25. R. E. Stratmann, G. E. Scuseria, and M. J. Frisch, “An efficient implementation of time-dependent density-functional theory for the calculation of excitation energies of large molecules,” Journal of Chemical Physics, vol. 109, no. 19, pp. 8218–8224, 1998. View at Publisher · View at Google Scholar · View at Scopus
  26. T. Yanai, D. P. Tew, and N. C. Handy, “A new hybrid exchange-correlation functional using the Coulomb-attenuating method (CAM-B3LYP),” Chemical Physics Letters, vol. 393, no. 1–3, pp. 51–57, 2004. View at Publisher · View at Google Scholar · View at Scopus
  27. W. J. D. Beenken and T. Pullerits, “Spectroscopic units in conjugated polymers: a quantum chemically founded concept?” Journal of Physical Chemistry B, vol. 108, no. 20, pp. 6164–6169, 2004. View at Scopus
  28. M. T. Sun, P. Kjellberg, W. J. D. Beenken, and T. Pullerits, “Comparison of the electronic structure of PPV and its derivative DIOXA-PPV,” Chemical Physics, vol. 327, no. 2-3, pp. 474–484, 2006. View at Publisher · View at Google Scholar · View at Scopus
  29. M. J. Frisch, G. W. Trucks, H. B. Schlegel, et al., Gaussian 09, Revision A.02, Gaussian, Wallingford, Conn, USA, 2009.
  30. N. M. O’Boyle and J. G. Vos, “GaussSum 1. 0,” Dublin City University, 2005, http://gausssum.sourceforge.net/.
  31. R. L. Martin, “Natural transition orbitals,” Journal of Chemical Physics, vol. 118, no. 11, pp. 4775–4777, 2003. View at Publisher · View at Google Scholar · View at Scopus
  32. D. Sundholm, S. Taubert, and F. Pichierri, “Calculation of absorption and emission spectra of [n]cycloparaphenylenes: the reason for the large Stokes shift,” Physical Chemistry Chemical Physics, vol. 12, no. 11, pp. 2751–2757, 2010. View at Publisher · View at Google Scholar · View at Scopus