
Figure Chart 1. Molecular structures of the studied molecules
Effect of Intermolecular Excited-state Interaction on Vibrationally Resolved Optical Spectra in Organic Molecular Aggregates
Wenqiang Li , Qian Peng , Yujun Xie , Tian Zhang , Zhigang Shuai
Organic materials have been a research hotspot owing to their potential applications in promising optoelectronic devices, such as organic light emitting diodes, [1] lasing, [2] field-effect transistors, [3] sensors, [4] solar cells, [5] and so on. The solid-state photophysics is different from solution owing to a variety of intermolecular interactions. The measurable optical spectrum is an effective means to reveal the molecularly oriented packing feature and the luminescent mechanism of organic molecules in aggregates. For instance, comparing with solution-phase optical spectra, the rod-shaped molecules always show a blue shift in absorption but red shift and quenching in emission for typical H-aggregates while their absorption and emission both are red-shifted and the emission intensity is noticeably enhanced in J-aggregates. Differently, for the typical aggregation induced emission (AIE) compounds, when going from solution to aggregate, the absorption peak is almost unchanged but emission is blue-shifted.[6] The fine structure of the vibrationally resolved optical spectrum can provide much information of molecular geometrical modification between the ground and excited states, electronic transition property, intermolecular interaction and resonance energy transfer (intermolecular excited-state interaction or exciton coupling) by comparing the solid-phase spectrum with the gas-phase or solution-phase ones.[7] Therefore, it is imperial to describe the vibrationally resolved spectrum by considering both the intramolecular and intermolecular interaction as well as the exciton coupling from first principles, and figure out the origin of the spectral peaks, which helps to deeply understand the mechanism of luminescence and to identify important molecular parameters governing the luminescent property of organic molecules in aggregates.
Kasha[8] has sketched out the spectral feature for H/J-aggregates by taking resonance energy transfer process into consideration on the basis of a simple model from a purely electronic state view. For H-aggregates, the absorption peak goes blue shift while the fluorescence is quenched when compared with those in solution. For J-aggregates, both absorption and emission are red-shifted and swelled in intensity. Recently, Spanoet al.[9] revealed the relationship between the intensity of 0-0 transition and H or J-aggregate conformation. Gierschner et al.[10] found the intermolecular Herzberg-Teller vibronic transitions by calculating the excited state property of the cluster consisting of several molecules at time-dependent density functional theory (TDDFT) level. Liang et al.[11] attributed the large red-shifted emission of PTCDI to the intermolecular charge transfer character. For flexible AIE systems, Wu et al.[12]claimed that the spectral character with unchanged absorption and blue-shifted emission upon crystallization are mainly caused by the intramolecular vibronic coupling and intermolecular electrostatic interaction at the level of quantum mechanics and molecular mechanics (QM/MM) theory.
In this contribution, we unveil the effect of exciton coupling on the vibrationally resolved solid-state optical spectral for organic molecules aggregates in crystalline state for the first time considering the intramolecular vibronic coupling, intermolecular excitonic coupling, and the intermolecular electrostatic interaction. Based on the Frenkel exciton model, we analyze the competitive relationship between the intramolecular vibronic coupling and intermolecular exciton coupling. And then we quantitatively calculate the optical spectra of a series of crystalline systems ranging from conventional non-AIEgens to novel AIEgens and reveal the mechanism of different spectrum nature in the non-AIEgens and AIEgens' aggregates, which can provide a clear relationship between structure and optical property.
To demonstrate the general behavior, we further investigated the optical spectra of a series of fluorophores, including four conventional non-AIEgens (DSB, [15] 6T, [16] Anthracene, [17] Rubrene[18]) and five AIEgens (DCDPP, [19] CB, [15a, 20] HPS, [15a, 21, 22a, 24] BFTPS, [22] BTPES[22a, 23]), see Chart 1. As known from the analysis above, the reorganization energy (λ) and exciton coupling (J) are two important competing factors in governing the spectral feature in aggregate. We first calculated the two parameters for all the studied molecular crystalline aggregates and present the results in Figure 4 (detailed data can be found in Table S1 of SI). Herein, the maximum exciton coupling value among different dimers in crystal is selected to characterize representatively the intermolecular interaction (detailed exciton coupling values of different dimers for all the studied systems can be found in Figure S2~S10 of SI).
It is interesting to observe from Figure 4 that the ratio J/λ for all the AIEgens is always less than 0.17 but for non-AIEgens J/λ is always greater than 0.17. Especially for DSB and 6T with typical aggregation caused quenching (ACQ) property, the J/λ ratio is well above the critical value.
Figure 5, Figure 6 and Figure S11 plot the calculated crystalline optical spectra with and without exciton coupling for the non-AIEgens and AIEgens at room temperature, respectively. In the spectrum calculations, the used parameters were given in Table S2~S5 of SI. As expected, the spectral fine structures of the non-AIEgens change apparently after the exciton coupling is taken into account, while the ones of AIEgens are almost unchanged. This fact confirms our previous observation that the influence of the exciton coupling becomes weak with the decrease of J/λ and tends to disappear as J/λ < 0.17. Among the non-AIEgens, the spectra of DSB and 6T experience the largest modification owing to very large J/λ, and the ones of anthracene and rubrene show slight excitonic effect because of small J/λ. The absorption spectra of DSB and 6T are sharply blue-shifted and their intensities are significantly enhanced, and the emissions are red-shifted with reduced intensity, which are in agreement with the typical H-type aggregation of DSB and 6T's crystals. For the AIEgens, the crystal packing is usually loose and the molecular geometry is usually flexible, leading to small J and large λ concomitantly, in another word J/λ < 0.17. Thus, their spectra are independent of the exciton coupling. This justifies our previous work on the AIEgens in which the radiative and nonradiative decay rates were calculated by neglecting the exciton couplings.[12,15a, 22a, 24]
In order to get deep insight into the effect of exciton interaction on the optical spectrum, we took DSB and HPS as examples to dissect the origin of changes in the emission spectrum by comparing the optical spectra in gas phase with those in cluster, taking intermolecular electrostatic interaction and excitonic coupling into accounts in Figure 7 and Figure 8 where the available experimental results are also depicted. Here, the intermolecular electrostatic interaction is considered during the geometry optimizations and frequency calculations by using QM/MM approach. That is, the intramolecular vibronic coupling in cluster is obtained based on the electronic structure information, which includes the influence of electrostatic interaction.
Figure 7. The calculated emission spectra of DSB in the gas phase, in a cluster with and without considering J at 298 K, as well as the experimental one in ref. 15b[15b]. The 921 0 means the 0-1 transition of mode 92, ω92=1647 cm-1 in the gas phase, ω92=1658 cm-1 in the crystal
Figure 8. The calculated emission spectra of HPS in the gas phase, in a cluster with and without considering J at 298 K, as well as the experimental one in ref. 24b[24b]. The 151 0 means the 0-1 transition of mode 15, ω2=21 cm-1, ω6=54 cm-1and ω11=70 cm-1 in the gas phase, ω15=137 cm-1and ω19=188 cm-1 in the crystal
It is obvious in Figure 7 that for DSB the emission spectra in three cases consist of two significant peaks. The two peaks are 0-0 peak and 0-1 peak arising from the normal mode of ca. 1650 cm-1 for the gas-phase DSB and cluster without considering J, respectively. The emission spectrum is red shifted from gas-phase to cluster without considering J owing to the decrease of 0-0 transition energy (△E0-0). When taking the exciton coupling into consideration, the 0-0 peak noticeably slumps owing to the offset of the two oppositely oriented transition dipole moment. And the transitions of lower FE state mixing electronic and vibrational states become dominant, which leads to more pronounced red shift, reproducing the experimental results.
Differently from DSB, the emission spectra of HPS (in Figure 8) are one broad band without any fine structure. The band is assigned as mixed 2-8, 2-9, 2-10, 0-3, 0-2, 1-2, and 1-1 transitions of low-frequency modes ( < 100 cm-1) for gas-phase HPS, and mixed 0-1, 0-2 and 1-0 transitions of low-frequency ( < 200 cm-1) for cluster without considering J. More strikingly, the emission spectrum is sharply enhanced in intensity and blue shifted from 495 nm to 484 nm when going from gas phase to cluster without J. This is easily understood that the intermolecular electrostatic interaction restricts molecular structural relaxation, namely out-of-plane twists, and decreases the vibronic coupling.[12a] In addition, the exciton coupling has no effect on the emission spectrum of HPS, as its J/λ < 0.17. This indicates that it is enough to reproduce the experimental result by considering the intermolecular electrostatic interaction. It is learned that the exciton coupling is necessary to be considered when predicting the optical spectra of the similar planar π-conjugated molecules with large J/λ. And the intermolecular electrostatic interaction always needs to be taken into consideration for the flexible AIE systems.
The vibronic coupling characterizes the intramolecular geometrical organization upon absorbing or emitting a phonon, which is always directly measured by the reorganization energy (λ,
Taking a dimer as an example, the basis function of the Frenkel exciton without considering vibration states can be written as
Here J is the exciton coupling; ε is the molecular excitation energy. By diagonalizing the matrix, we can get the eigenvalue and eigenvector,
The constructed eigenstates are called as the Frenkel exciton (FE) states. The dependence of optically spectral feature on J is numerously examined by setting ε=3.0 eV at 298 K in Figure 1. For the H-aggregates, the transition dipole moment is arranged parallels and the resultant exciton coupling J is positive, so the lowest FE state is a dipole-forbidden transition and the fluorescence is quenched according to Kasha rule.[13] The absorption processes appear in higher dipole-allowed FE states and exhibit remarkable blue shift. While, for the J-aggregates, the J is negative, resulting in the dipole-allowed transition of the lower FE state. Consequently, and the absorption and emission are both red-shifted and the emission intensity is sharply enhanced. It is also seen that the shifts of the spectra become very serious when the J>20 meV.
As known, the vibronic coupling is essential in determining the spectral shape.[14] For simplification, a dimer system with one normal mode with two quantum number (ν=0 and 1) for each molecule is considered under the vibronic coupled exciton model. The model exciton Hamiltonian is constructed as
Here, the
In order to figure out the influence of the vibronic coupling, we calculated the optical spectra of the model dimer as well as that of the monomer for comparison by setting ε=3.0 eV, J=±50 meV, ω=1400 cm-1, S=1, T=298 K in Figure 2. It is clearly seen that the absorption and emission spectra of the monomer both exhibit two main peaks which correspond to the 0-0 and 0-1 transition, respectively. Whereas the absorption and emission spectra of dimer show distinct feature for the different molecular packing. The H-aggregate produces the blue-shifted absorption with maximum peak assigned as 4th FE state, while the J-aggregate generates the red-shifted one with maximum peak dressed as 1st FE state. The emission spectra of the H-aggregate is largely red shifted (
In this work, we systematically investigated the role of intramolecular vibronic coupling and intermolecular interaction in determining the optical spectrum features. Starting from vibronic-coupled exciton models, we reveal the role of the vibronic coupling in the fine structure of the optical spectrum. Then we examined the competitive relationship of the intramolecular vibronic coupling (λ) and the intermolecular exciton coupling (J) to dominate the crystalline state spectral characters. It is found that the impact of the exciton coupling on the spectra decreases with the ratio of J/λ, and goes to null when J/λ < 0.17 at room temperature.
Furthermore, we compare the crystalline optical spectra of a series of non-AIEgens with AIEgens at room temperature. It is found that the spectra of the non-AIEgens are strongly dependent on the intermolecular exciton coupling but for AIEgens the ratio J/λ is always less than 0.17, leading to the independence of optical spectra on the exciton coupling. The spectra of typical non-AIEgens DSB and 6T undergo remarkable modification from the gas phase to aggregate owing to their very large J/λ, and for anthracene and rubrene, the modification is small for the small J/λ value close to 0.17. But for AIEgens, the optical spectra are remarkably insensitive to the exciton coupling and the intramolecular vibronic coupling with intermolecular electrostatic interaction dominate the photophysical processes. These results help to understand the differences between AIEgens and non-AIEgens.
In the organic molecule aggregates, the Frenkel exciton (FE) model is usually adopted to describe the exciton coupling[8, 9]. The basis function of the aggregate (containing p monomers) at the excited state can be expressed as:
with
Assuming the harmonic oscillator and one-particle approximation (in which the vibration states of the excited molecule are active while the ones of other molecules are frozen[9]), the elements in Eq. (4) are recast as,
Here, m/n is the index of the excited monomer;
Here
By diagonalizing the matrix Eq. (4), we can obtain the eigenvalue (eigenenergy) and eigenvector (eigenstate) of the excited states. Then, we can calculate the absorption and emission spectra based on Eq. (8) and Eq. (9)[14]:
Here,
The quantum calculations of all the studied compounds are performed out by using QM/MM approach, including geometrical optimization, electronic structure, vibration frequency, transition property and so on. The QM/MM models are set up by cutting out a cluster from the X-ray crystal structure and prescribing the central molecule as QM part and its surrounding ones as MM part. Taking distyrylbenzene (DSB) as an example, the QM/MM cluster (Chart 2) size is 5×5×3 and the thirdly nearest molecules are taken into account in the compact herringbone directions (total 105 molecules). (TD) B3LYP[26] and 6-31G (d) were used for QM calculation. General Amber force field (GAFF) was adopted to describe the MM molecules.[27] The QM/MM calculation was accomplished by using ChemShell 3.5 package, [28] which interfaces the Turbomole 6.5 program[29] for QM and DL_POLY program[30] for MM. It should be noted that the QM molecule was active while the MM molecules were kept frozen and the electrostatic embedding scheme with QM polarization was adopted. In addition, the exciton couplings (J) were calculated at the level of CAM-B3LYP[31]/6-31G (d) by using our home-built MOMAP program[32] coupling with the NWchem 6.3 program.[33]
A Frenkel exciton computational model consists of a central monomer and its most adjacent molecules in the crystal (8~9 molecules, shown in Figure S2~S10 of SI). The optical spectra were evaluated by employing the vibronic coupled Frenkel exciton theory in MOMAP program. Herein, the Lorentzian broadening was applied with a constant γ=1100 cm-1 for non-AIEgens, [34] and 2000 cm-1 for AIEgens to reproduce the spectrum with fully vibrational states obtained by the vibration correlation function approach at room temperature.[24a] To balance the computational cost and accuracy, the important normal modes were selected (see details in Table S2~S5), and its vibrational states were cut off when the FC integral fvk-0 less than 1% of the maximum.
Diagonalization of model exciton Hamiltonian in Eq. (3), the reorganization energy and exciton coupling of the selected fluorophores; The temperature dependencies of optical spectra in the H-and J-aggregates; the Frenkel exciton computational model cutting off from crystal and the relative exciton coupling values; optical spectra of the AIEgens (BFTPS and BTPES) in cluster with and without exciton coupling at 298 K; the selected normal modes and their quantum number of the Non-AIEgens and AIEgens.
Tang, C. W.; VanSlyke, S. A. Appl. Phys. Lett. 1987, 51, 913; (b) Burroughes, J. H.; Bradley, D. D. C.; Brown, A. R.; Marks, R. N.; Mackay, K.; Friend, R. H.; Burns, P. L.; Holmes, A. B. Nature 1990, 347, 539; (c) Malissa, H.; Kavand, M.; Waters, D. P.; van Schooten, K. J.; Burn, P. L.; Vardeny, Z. V.; Saam, B.; Lupton, J. M.; Boehme, C. Science 2014, 345, 1487; (d) Lee, J.; Chen, H.-F.; Batagoda, T.; Coburn, C.; Djurovich, P. I.; Thompson, M. E.; Forrest, S. R. Nat. Mater. 2016, 15, 92.
Schäfer, F. P.; Schmidt, W.; Volze, J. Appl. Phys. Lett. 1966, 9, 306; (b) Morales-Vidal, M.; Boj, P. G.; Villalvilla, J. M.; Quintana, J. A.; Yan, Q.; Lin, N.-T.; Zhu, X.; Ruangsupapichat, N.; Casado, J.; Tsuji, H.; Nakamura, E.; Diaz-Garcia, M. A. Nat. Commun. 2015, 6, 8458.
Horowitz, G. Adv. Mater. 1998, 10, 365; (b) Liu, J.; Zhang, H.; Dong, H.; Meng, L.; Jiang, L.; Jiang, L.; Wang, Y.; Yu, J.; Sun, Y.; Hu, W.; Heeger, A. J. Nat. Commun. 2015, 6, 10032.
Gaylord, B. S.; Heeger, A. J.; Bazan, G. C. PNAS 2002, 99, 10954; (b) Rana, S.; Elci, S. G.; Mout, R.; Singla, A. K.; Yazdani, M.; Bender, M.; Bajaj, A.; Saha, K.; Bunz, U. H. F.; Jirik, F. R.; Rotello, V. M. J. Am. Chem. Soc. 2016, 138, 4522.
Yu, G.; Gao, J.; Hummelen, J. C.; Wudl, F.; Heeger, A. J. Science 1995, 270, 1789; (b) Page, Z. A.; Liu, Y.; Duzhko, V. V.; Russell, T. P.; Emrick, T. Science 2014, 346, 441; (c) Holliday, S.; Ashraf, R. S.; Wadsworth, A.; Baran, D.; Yousaf, S. A.; Nielsen, C. B.; Tan, C. H.; Dimitrov, S. D.; Shang, Z. R.; Gasparini, N.; Alamoudi, M.; Laquai, F.; Brabec, C. J.; Salleo, A.; Durrant, J. R.; McCulloch, I. Nat. Commun. 2016, 7, 11; (d) Rong, Y.; Mei, A.; Liu, L.; Li, X.; Han, H. Acta Chim. Sinica 2015, 73, 237. (荣耀光, 梅安意, 刘林峰, 李雄, 韩宏伟, 化学学报, 2015, 73, 237.) (e) Fu, Y.-T.; Yi, Y.; Coropceanu, V.; Risko, C.; Aziz, S. G.; Brédas, J.-L. Sci. China Chem. 2014, 57, 1330.
Mei, J.; Leung, N. L. C.; Kwok, R. T. K.; Lam, J. W. Y.; Tang, B. Z. Chem. Rev. 2015, 115, 11718. doi: 10.1021/acs.chemrev.5b00263
Valeur, B.; Berberan-Santos, M. N. Molecular Fluorescence, Wiley-VCH, Weinheim, 2012, pp. 141~179.
Kasha, M. Radiat. Res. 1963, 20, 55; (b) Kasha, M.; Rawls, H.; El-Bayoumi, M. A. Pure Appl. Chem. 1965, 11, 371.
Spano, F. C. Acc. Chem. Res. 2010, 43, 429; (b) Spano, F. C. J. Chem. Phys. 2003, 118, 981; (c) Spano, F. C. Phys. Rev. B 2005, 71, 094110; (d) Spano, F. C. Annu. Rev. Phys. Chem. 2006, 57, 217.
Wykes, M.; Parambil, R.; Beljonne, D.; Gierschner, J. J. Chem. Phys. 2015, 143, 114116; (b) Gierschner, J.; Ehni, M.; Egelhaaf, H.-J.; Milián Medina, B.; Beljonne, D.; Benmansour, H.; Bazan, G. C. J. Chem. Phys. 2005, 123, 144914.
Gao, F.; Liang, W.; Zhao, Y. Sci. China Chem. 2010, 53, 297. doi: 10.1007/s11426-010-0075-2
Wu, Q.; Zhang, T.; Peng, Q.; Wang, D.; Shuai, Z. Phys. Chem. Chem. Phys. 2014, 16, 5545. (b) Wu, Q.; Peng, Q.; Zhang, T.; Shuai, Z. Sci. China Chem. 2013, 43, 1078.
Kasha, M. Discuss. Faraday Soc. 1950, 9, 14. doi: 10.1039/df9500900014
Niu, Y.; Peng, Q.; Deng, C.; Gao, X.; Shuai, Z. J. Phys. Chem. A 2010, 114, 7817;
Zhang, T.; Peng, Q.; Quan, C.; Nie, H.; Niu, Y.; Xie, Y.; Zhao, Z.; Tang, B. Z.; Shuai, Z. Chem. Sci. 2016, 7, 5573; (b) Wu, C. C.; Korovyanko, O. J.; Delong, M. C.; Vardeny, Z. V.; Ferraris, J. P. Synth. Met. 2003, 139, 735.
Yassar, A.; Horowitz, G.; Valat, P.; Wintgens, V.; Hmyene, M.; Deloffre, F.; Srivastava, P.; Lang, P.; Garnier, F. J. Phys. Chem. 1995, 99, 9155; (b) Stradomska, A.; Petelenz, P. J. Chem. Phys. 2009, 130, 094705.
Mason, R. Acta Crystallogr. 1964, 17, 547; (b) Pope, M.; Kallmann, H. P.; Magnante, P. J. Chem. Phys. 1963, 38, 2042; (c) Li, H.; Duan, L.; Zhang, D.; Dong, G.; Wang, L.; Qiu, Y. Sci. China Ser. B:Chem. 2009, 52, 181.
Gao, F.; Liang, W. Z.; Zhao, Y. J. Phys. Chem. A 2009, 113, 12847; (b) Mitrofanov, O.; Kloc, C.; Siegrist, T.; Lang, D. V.; So, W.-Y.; Ramirez, A. P. Appl. Phys. Lett. 2007, 91, 212106.
Wu, Q.; Deng, C.; Peng, Q.; Niu, Y.; Shuai, Z. J. Comput. Chem. 2012, 33, 1862; (b) Qin, A.; Lam, J. W. Y.; Mahtab, F.; Jim, C. K. W.; Tang, L.; Sun, J.; Sung, H. H. Y.; Williams, I. D.; Tang, B. Z. Appl. Phys. Lett. 2009, 94, 253308.
Dong, Y.; Lam, J. W. Y.; Qin, A.; Sun, J.; Liu, J.; Li, Z.; Sun, J.; Sung, H. H. Y.; Williams, I. D.; Kwok, H. S.; Tang, B. Z. Chem. Commun. 2007, 31, 3255.
Chen, J.; Law, C. C. W.; Lam, J. W. Y.; Dong, Y.; Lo, S. M. F.; Williams, I. D.; Zhu, D.; Tang, B. Z. Chem. Mater. 2003, 15, 1535. doi: 10.1021/cm021715z
Xie, Y.; Zhang, T.; Li, Z.; Peng, Q.; Yi, Y.; Shuai, Z. Chem. Asian J. 2015, 10, 2154; (b) Zhan, X.; Haldi, A.; Risko, C.; Chan, C. K.; Zhao, W.; Timofeeva, T. V.; Korlyukov, A.; Antipin, M. Y.; Montgomery, S.; Thompson, E.; An, Z.; Domercq, B.; Barlow, S.; Kahn, A.; Kippelen, B.; Bredas, J.-L.; Marder, S. R. J. Mater. Chem. 2008, 18, 3157.
Zhao, Z.; Liu, D.; Mahtab, F.; Xin, L.; Shen, Z.; Yu, Y.; Chan, C. Y. K.; Lu, P.; Lam, J. W. Y.; Sung, H. H. Y.; Williams, I. D.; Yang, B.; Ma, Y.; Tang, B. Z. Chem.-Eur. J. 2011, 17, 5998. doi: 10.1002/chem.v17.21
Zhang, T.; Jiang, Y.; Niu, Y.; Wang, D.; Peng, Q.; Shuai, Z. J. Phys. Chem. A 2014, 118, 9094. (b) Zhan, X. W.; Risko, C.; Amy, F.; Chan, C.; Zhao, W.; Barlow, S.; Kahn, A.; Bredas, J.-L.; Marder, S. R. J. Am. Chem. Soc. 2005, 127, 9021.
Hsu, C.-P.; You, Z.-Q.; Chen, H.-C. J. Phys. Chem. C 2008, 112, 1204. doi: 10.1021/jp076512i
Becke, A. D. J. Chem. Phys. 1993, 98, 5648; (b) Lee, C.; Yang, W.; Parr, R. G. Phys. Rev. B 1988, 37, 785.
Wang, J.; Wolf, R. M.; Caldwell, J. W.; Kollman, P. A.; Case, D. A. J. Comput. Chem. 2004, 25, 1157. doi: 10.1002/(ISSN)1096-987X
Sherwood, P.; de Vries, A. H.; Guest, M. F.; Schreckenbach, G.; Catlow, C. R. A.; French, S. A.; Sokol, A. A.; Bromley, S. T.; Thiel, W.; Turner, A. J.; Billeter, S.; Terstegen, F.; Thiel, S.; Kendrick, J.; Rogers, S. C.; Casci, J.; Watson, M.; King, F.; Karlsen, E.; Sjøvoll, M.; Fahmi, A.; Schäfer, A.; Lennartz, C. J. Mol. Struct. THEOCHEM 2003, 632, 1. doi: 10.1016/S0166-1280(03)00285-9
TURBOMOLE V6.52013, University of Karlsruhe and of the Forschungszentrum Karlsruhe GmbH, 1989-2007; TURBOLE GmbH, since 2007(accessed May 23, 2013); (b) Ahlrichs, R.; Bär, M.; Häser, M.; Horn, H.; Kölmel, C. Chem. Phys. Lett. 1989, 162, 165.
Smith, W.; Forester, T. R. J. Mol. Graphics 1996, 14, 136. doi: 10.1016/S0263-7855(96)00043-4
Yanai, T.; Tew, D. P.; Handy, N. C. Chem. Phys. Lett. 2004, 393, 51. doi: 10.1016/j.cplett.2004.06.011
Shuai, Z.; Peng, Q.; Niu, Y.; Geng, H.; MOMAP, Revision 0.3.001 ed.; MOMAP:a free and open-source molecular materials property prediction package; avaliable online:http://www.shuaigroup.net, Beijing, China, 2016.
Valiev, M.; Bylaska, E. J.; Govind, N.; Kowalski, K.; Straatsma, T. P.; Van Dam, H. J. J.; Wang, D.; Nieplocha, J.; Apra, E.; Windus, T. L.; de Jong, W. A. Comput. Phys. Commun. 2010, 181, 1477. doi: 10.1016/j.cpc.2010.04.018
Gierschner, J.; Mack, H. G.; Egelhaaf, H. J.; Schweizer, S.; Doser, B.; Oelkrug, D. Synth. Met. 2003, 138, 311. doi: 10.1016/S0379-6779(03)00030-4
Figure 7 The calculated emission spectra of DSB in the gas phase, in a cluster with and without considering J at 298 K, as well as the experimental one in ref. 15b[15b]. The 921 0 means the 0-1 transition of mode 92, ω92=1647 cm-1 in the gas phase, ω92=1658 cm-1 in the crystal
Figure 8 The calculated emission spectra of HPS in the gas phase, in a cluster with and without considering J at 298 K, as well as the experimental one in ref. 24b[24b]. The 151 0 means the 0-1 transition of mode 15, ω2=21 cm-1, ω6=54 cm-1and ω11=70 cm-1 in the gas phase, ω15=137 cm-1and ω19=188 cm-1 in the crystal