Yu
Harabuchi
ab,
Tetsuya
Taketsugu
ab and
Satoshi
Maeda
*ab
aDepartment of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan. E-mail: smaeda@mail.sci.hokudai.ac.jp
bCREST, Japan Science and Technology Agency, Tokyo 102-8666, Japan
First published on 30th July 2015
Minimum energy conical intersection (MECI) geometries were searched for polycyclic aromatic hydrocarbons (PAHs) consisting of up to 26 atoms. The energy barriers to the obtained MECIs showed a correlation with their fluorescence quantum yields. This provides a theoretical rationale for the size dependence of the fluorescence quantum yields seen in these PAHs.
The CI region also plays a role in non-ultrafast photoreactions. For example, fluorescent molecules stay on the electronic excited-state potential energy surface (PES) on sub-nano- to micro-second timescales. In this case, radiative and non-radiative decays compete with each other, and their competition determines the fluorescence quantum yield. When the energy barriers to reach CI regions are sufficiently high, and are thus difficult to access from the FC region, the fluorescence process is dominant. However, AIMD simulations of the electronic excited-state PES on sub-nano- to micro-second timescales remain highly demanding, even for small molecules. Thus, in this study, we tackled the problem by directly finding minimum energy conical intersections (MECIs).
The MECI is the lowest energy point within a CI hyperspace, and MECI geometries have been optimized as a critical point, or a representative geometry, of the nonadiabatic transition.11–18 However, preparing initial geometries for MECI optimizations is not trivial, because MECI geometries are generally distant from the stable geometries in the ground state. Motivated by this, the seam model function (SMF) approach was proposed to find the all-important MECIs.19 In the SMF, automated stable structure search methods are initially applied to a model function on which approximate MECI structures can be obtained as local minima. These approximate structures are then optimized using a standard MECI optimization technique. This approach has been used extensively in the studies of the photodissociation mechanisms of small molecules.20 Recently, its applicability has been improved dramatically by the use of the spin–flip time dependent density functional theory (SF-TDDFT)21–23 in the electronic structure calculations and the use of the single-component artificial force induced reaction (SC-AFIR)24 method in the automated stable structure search.25,26
In the present study, MECI geometries for the lowest two singlet electronic states (denoted S1/S0-MECI) were systematically searched using the SMF/SC-AFIR/SF-TDDFT approach. These searches were applied to the five polycyclic aromatic hydrocarbons (PAHs) shown in Fig. 1 and initiated from their FC points. Transition state (TS) geometries between the FC region and the obtained MECIs were also determined. Based on the energies of the obtained MECI and TS structures, the size dependence of the fluorescence quantum yields seen in these PAHs is discussed.
Fig. 1 Molecular geometries of five PAHs: benzene, naphthalene, phenanthrene, anthracene, and pyrene. The experimental fluorescence quantum yields (measured in cyclohexane) are also shown for each molecule.34 |
The S1/S0-CI regions and the S1/S0-MECI structures of benzene have been studied thoroughly by Robb and co-workers.35–37 However, the S1/S0-MECIs for the larger PAHs have not been studied in detail. In this study, we therefore systematically explored the MECIs of the larger PAHs using the SMF/SC-AFIR/SF-TDDFT approach, and we discuss here the efficiency of the non-radiative decay.
The S1MIN, the three most stable MECIs, and the related TSs for each PAH are summarized in Fig. 2. A low-energy MECI with a CH bond dissociation for pyrene was excluded, because a high-energy TS providing a route to the CI region was found. The most stable MECI geometry for benzene (shown in Fig. 2) was consistent with that determined in previous studies.35–37 The FC point and the three most stable MECIs for benzene were further optimized using the CASPT2/cc-pVDZ method. The S1 energies determined using the CASPT2 method for the optimized FC, MECI1, MECI2 and MECI3 geometries of benzene were 4.84 eV, 4.97 eV, 5.07 eV, and 5.17 eV, respectively, relative to the ground state benzene. The geometries and the order of the relative energies were consistent with the SF-BHHLYP results shown in Fig. 2 (see the ESI† for the geometries).
Fig. 2 S1-minimum geometries near the FC point (S1MIN), the three most stable MECI geometries, and the TS geometries to reach the lowest three MECIs from S1MIN for the five PAHs benzene, naphthalene, phenanthrene, anthracene, and pyrene. All the geometries were optimized using SF-BHHLYP/6-31G*. For each system, the potential energy of each geometry is given relative to the S0 energy in the ground state equilibrium structure. All the obtained geometries are summarized in the ESI.† The S1MIN of phenanthrene could not be optimized at SF-BHHLYP, and the TS3 of anthracene did not exist, so these entities are shown as “N/A”. |
As shown in Fig. 2, the S1MIN geometries were planar for all of the PAHs. In Fig. 2, the label “N/A” for the S1MIN of phenanthrene indicates that the corresponding local minimum was not found. That is, the optimizers implemented in both GRRM and GAMESS showed oscillating behavior, and ultimately failed to reach a converged geometry. This oscillation was caused by the instability of the SF-TDDFT calculations for geometries in which the reference triplet state was nearly degenerate with the other triplet state. Although most of the MECIs were accessed from the S1MIN in one step via a single TS, there were two cases showing two barriers, and one case in which there was no barrier. The MECI2 of benzene and the MECI3 of naphthalene fell into the two-barrier classification; for these two, only the higher TS is shown in Fig. 2. In the two-step path to the MECI2 of benzene, MECI3 was the intermediate, and MECI2 and MECI3 were connected by a small barrier. Among the two TSs along this path, the first TS between S1MIN and MECI3 was higher than the second TS between MECI3 and MECI2. Thus, in Fig. 2, the first TS, which was identical to TS2, is presented as TS3. Similarly, for the MECI3 of naphthalene, S1MIN and MECI3 were connected by two steps via the fourth stable MECI (see MECI4 in the ESI†). The TS connecting S1MIN and MECI4 was higher than that between MECI4 and MECI3, and the former TS is shown as TS3 in Fig. 2. The MECI3 of anthracene corresponded to the case with no barrier; the corresponding TS3 is indicated by “N/A” in Fig. 2. This was confirmed by the meta-IRC (the mass-weighted steepest descent path) determined from the MECI3 of anthracene; the meta-IRC led directly to S1MIN.
From the extensive list of MECI structures shown here, it was found that the stable MECI geometries for the PAHs could be categorized into three types characterized by the three stable MECI structures determined for benzene. In other words, the structural deformations in the stable MECIs were localized on a single ring in these small PAHs. Consequently, the relative energies for the most stable MECIs fell within a narrow range of 4.94–5.57 eV, as shown in Fig. 2.
Here, we discuss the magnitude of the fluorescence quantum yields based on the shape of potential energy surfaces. Fig. 3A and B show schematic illustrations of two different potential profiles: one illustrates the case in which the non-radiative decay is dominant (Fig. 3A), and the other illustrates the case in which the radiative decay with fluorescence is dominant (Fig. 3B). Fig. 3A shows an MECI in a lower energy region, relative to the FC point; in this case a non-radiative, ultrafast decay via the related CI region can take place. In the case of Fig. 3B, the lowest energy MECI is much higher in energy than the FC point. Hence, the access frequency of a molecule to the CI region is low, and the non-radiative decay should thus be slow. In other words, the relative importance of the radiative decay is high. In this manner, one can discuss the relative importance of the radiative decay by comparing the energy difference between the most stable MECI and the FC point, i.e., ELowest–MECI − EFC. For a more quantitative prediction, the S1MIN and the TS connecting the S1MIN and the CI region are required, as illustrated in Fig. 3B. Based on the transition state theory,38 the Gibbs free energy gap between the S1MIN and the TS, i.e. GS1TS – GS1MIN, can be compared directly with the access frequency of a molecule to the CI region. In this study, the ELowest–MECI – EFC and GS1TS – GS1MIN gaps were compared with the experimental fluorescence quantum yields, where GS1TS – GS1MIN corresponds to the Gibbs free-energy gap at 298.15 K and 1 atm.
With this information in mind, the experimental fluorescence quantum yields were plotted against the two calculated energy differences, i.e., the ELowest–MECI – EFC gap and the GS1TS – GS1MIN gap, for the five PAHs, as shown in Fig. 4. The GS1TS – GS1MIN gap for phenanthrene was excluded from Fig. 4, because the S1MIN geometry was not found at the SF-BHHLYP level, as noted above. As expected, the GS1TS – GS1MIN gap showed a clear correlation with the experimental fluorescence quantum yields. For these PAHs, the ELowest–MECI – EFC gap also correlated well with the fluorescence quantum yields. This was because the TSs were close to the MECIs, and the S1MIN were also close to the FC point. Thus, in the following, the origin of the size dependence of fluorescence quantum yields is discussed based on the nature of the FC point and the MECIs.
Fig. 4 Diagram illustrating the correlation between the fluorescence quantum yields measured in cyclohexane34 and the two energy gaps, ELowest–MECI – EFC (cross marks) and GS1TS – GS1MIN (circles). Each PAH is indicated with the name of the PAH. |
For these PAHs, EFC or the experimental excitation energy decreases with increase in the molecular size because of the strong stabilization of the π–π* states that results from the delocalization of the π and π* orbitals. However, the observation that ELowest–MECI did not differ very much for the different PAHs was new; this phenomenon occurred because the structural deformations in the stable MECIs were localized on a single ring, as shown in Fig. 2. Based on this result, we proposed that the delocalized and localized characters of the S1 at the FC point, and the S1 at the lowest MECI geometry, allowed the larger PAHs to emit stronger fluorescence for the PAHs studied here.
We note that in this study we did not consider a number of factors, including the rate of radiative decay, the contributions of the triplet states, and the weak vibronic coupling around the S1 local minimum. A more quantitative analysis considering all of these factors together will be performed in a future study working toward the quantitative prediction of fluorescent quantum yields. Finally, we emphasize that the present approach is promising for illustrating the impact of the nonadiabatic transition near the CI in fluorescence quenching.
Based on the extensive MECI structural database determined here, an interesting correlation between the fluorescence quantum yields and the energy barrier to the CI region was discovered for these five small PAHs. This suggested that two processes—specifically, radiative decay with the emission of fluorescence and nonradiative decay near the conical intersection—were in competition. This correlation provided a theoretical explanation for the observed experimental trend that the larger PAHs showed stronger fluorescence. Further applications of the present approach to other fluorescent systems are in progress to investigate the generality of this correlation.
Footnote |
† Electronic supplementary information (ESI) available: All of the MECI geometries and the relative energies are shown. See DOI: 10.1039/c5cp02103k |
This journal is © the Owner Societies 2015 |