Rajat
Walia
,
Zexiang
Deng
and
Jun
Yang
*
Department of Chemistry, The University of Hong Kong, Pokfulam Road, Hong Kong, P. R. China. E-mail: juny@hku.hk
First published on 30th August 2021
Singlet fission duplicates triplet excitons for improving light harvesting efficiency. The presence of the interaction between electronic and nuclear degrees of freedom complicates the interpretation of correlated triplet pairs. We report a quantum chemistry study on the significance and subtleties of multistate and multimode pathways in forming triplet pair states of the pentacene dimer through a six-state vibronic-coupling Hamiltonian derived from many-electron adiabatic wavefunctions of an ab initio density matrix renormalization group. The resulting spin values of the singlet manifolds on each pentacene center are computed, and the varying spin nature can be distinguished clearly with respect to dimer stacking and vibronic progression. Our monomer spin assignments reveal the coexistence of both lower-lying weak and higher-lying strong charge transfer states which interact vibronically with the triplet pair state, providing important implications for its generation and separation occurring in vibronic regions. This work conveys the importance of the many-electron process requiring close low-lying singlet manifolds to determine the subtle fission details, and represents an important step for understanding vibronically resolved spin states and conversions underlying efficient singlet fission.
A set of general theoretical guidelines was proposed for the low-lying excited states of SF chromophores, i.e., E(S) − 2E(T) ≥ 0 to make the whole process isoergic (or slightly exoergic), and E(T2) − 2E(T) > 0 to compete with an inverse phenomenon, triplet–triplet annihilation.5 A few organic compounds have been identified to produce triplet yields more than 100%, such as polyacenes (anthracene, tetracene, pentacene, and hexacene), 1,3-diphenylisobenzofuran, perylene derivatives, and biradicaloids.6–9
As pointed out by Casanova in a recent seminal review article,10 “in singlet fission, it is commonly assumed that if the optical state is not the lowest excited singlet (S1), the system rapidly decays to S1 (Kasha's rule) before it undergoes fission, …”. The most interesting state in these electronically excited molecules for SF applications is the low-lying state which has a predominant double excitation character. This state is vaguely understood as a composition of two triplets (on different monomers) coupled in a singlet and usually called a correlated triplet pair state as 1TT, that is generated at an ultrafast time scale (∼100 fs) and has a tendency to split into two distinct triplets.11 Although the energetic criteria for pentacene and other SF chromophores have already been clearly established,12,13 a quantitative understanding of the photophysical nature of TT states arising in SF materials has been difficult from experiments alone, which causes confusion in unambiguous assignments and proper analysis of TT spectral signatures, and complicates the interpretation of the SF mechanism.14
Pentacene is the smallest polyacene that undergoes exoergic SF and has thus been extensively investigated in both theoretical12,15–20 and experimental SF studies.11,21–27 These studies have demonstrated diverse mechanistic pathways in the SF generation process which require an active interaction between single- and multi-excitonic states in one way or another, such as the inter-conversion of a local excitation (LE) singlet S1 state with a nearby dark state through the excimer formation,16 a superexchange mechanism via a virtual charge-transfer (CT) state,28 and an indirect quantum coherent mechanism involving low-lying CT and multiexcitonic diabatic states.19,29,30 Although the CT-mediated mechanism appears to be the largely accepted explanation for SF in tetracene and pentacene crystals,10 more detailed studies have also revealed various vibronic origins for making a coherent TT population, ranging from a single vibrational mode31 through a collection of vibrations in both low- and high-frequency regions,32,33 to the vibrational coherency.34–36 Interestingly, the vibronic coupling (VC) between two CT states has been shown to enhance SF dynamics by reducing the two CTs' destructive interference.37 The crystalline polarizable environment was deduced to stabilize vibronically tuned CT states, causing a significant acceleration in SF generation rates with a slight diminution in the TT yield.38
Conceptually, the presence of strongly correlated triplet pairs in SF chromophores requires a high-level ab initio multiconfigurational wavefunction39,40 for which the number of significant determinants grows exponentially with correlated electrons. Dealing with such strongly correlated systems becomes intractable, which limits correlation to only a small number of electrons in these molecules, leaving many profound questions unresolved despite already intriguing findings. For instance, the part of the molecule which can be statically correlated with traditional complete active space self-consistent field (CASSCF) methods remains very small in large molecules particularly in acenes larger than naphthalene and these methods may not provide an accurate picture of the SF pathway involving two or more such coupled chromophores. Although the ideal CASSCF-based calculations should treat complete π-valence electrons, due to the aforementioned challenges, a practical choice of active space should be at least sufficiently large to account for most valence π → π* excitation spectra of SF chromophores. The use of ab initio density matrix renormalization group (ab-DMRG) method provides a way of correlating a larger number of excited electrons in only polynomial costs, particularly effective in obtaining numerically exact wavefunctions for strongly correlated one- and quasi-one-dimensional systems,41,42e.g., molecules such as polyenes and metallic hydrogen chains by using up to (100o, 100e) active spaces.43
In this work, we applied the ab-DMRG algorithm to correlate a large part of the π-valence space and investigate the low-lying excited states of pentacene. These electronic states in the pentacene dimer are characterized using state-average density matrix renormalization group self-consistent field (DMRGSCF)44 calculations. Furthermore, we implemented an extensive wavefunction analysis tool based on DMRG one- and two-particle state and transition density matrices to reveal the role of these low-lying adiabatic states in explaining the essence of the SF mechanism. Specifically, we used DMRG wavefunctions to calculate the projected local spin distribution for monomer fragments in the dimer system, the multipole couplings to quantify the role of intermediate CT states electrostatically, and a VC analysis to understand molecular vibrations facilitating SF state conversion. This study serves to not only reconcile the need of multiple excitonic states and vibrational modes which are shown to be critical to revealing the SF mechanistic details in pentacene, but also to extend a general high-level computational approach and modeling tool for systematically analyzing potential candidates that harvest efficient SF.
In what follows, to avoid ambiguity, the mechanistic discussions on the formation and separation of the correlated triplet–triplet (TT) pair state are differentiated, while the abbreviation “SF” is generally applied to the whole process, as recommended by Smith and Michl.5 Moreover, we denote the lowest adiabatic excited singlet manifold as the dimeric S1 without linking to its character, the next singlet state as the dimeric S2 and so on, as recommended by Casanova.10
Fig. 1 The reference dimer geometry for SF applications: two monomers are used to form the dimer model similar to the crystal structure in pentacene polymorphs. |
An analysis of the DMRGSCF/NEVPT2 energetics of the dimer model over a range of orbitals reveals that the (22π,22e) active space is suitable to effectively replicate excited state energies of low-lying singlets as compared to those from the GW/BSE method. The DMRGSCF(22π,22e)/NEVPT2 results account for the CAS-based static correlation, corrected with the NEVPT2 dynamic correlations on top of the DMRGSCF(12π,12e) wavefunction due to the complexity of NEVPT2 calculations for the cc-pVTZ basis set at the DMRG(22π,22e) reference wavefunction. These energies for the first five excited singlet states of the pentacene dimer and solid are listed in Table 1. The detailed analyses of excited state energies and the choice of active space in both the monomer and dimer are given in Section S1 of the ESI.† The first four excited states of the pentacene dimer are close-lying within an energy window of ∼0.3 eV, indicating that the wavefunction delocalizes over the dimer when two monomers approach each other. Similar excited state energies were obtained by Zimmerman et al. using the restricted active space double spin-flip (RAS-2SF) method for a pentacene dimer, where RC–C = 5.65 Å.18 The delocalization of the wavefunction is more enhanced in a crystal environment and this energy window further reduces to ∼0.2 eV.
State | Dimer (22π,22e) | Pentacene crystal | Character | 〈Ŝ2〉 |
---|---|---|---|---|
S1 | 1.85 | 1.86 | TT | 1.82 |
S2 | 1.93 | 1.88 | CT1 | 0.29 |
S3 | 1.99 | 1.97 | LE | 0.00 |
S4 | 2.16 | 2.05 | CT2 | 0.62 |
S5 | 2.84 | 2.06 | LE | 0.00 |
These SF states in the dimer possess distinct spin and spatial characteristics and interact with one another differently. Luzanov et al. determined the spin nature for various dimers of SF relevance (e.g., 1,3-diphenylisobenzofuran (DPBF) and tetracene) based on the weights of LE, charge-resonance (CR) and TT spin-adapted configurations employing a fragment spin-correlator.59 In this work, we reformulate the well-established60–62 local spin decomposition of total spin square operator 〈Ŝ2〉 in terms of the local one- and two-particle density matrices for exact adiabatic wavefunctions (see eqn (S3) in the ESI†). Fig. 2(a) shows the 〈Ŝ2〉 distribution for a fragment in the pentacene dimer as a function of inter-dimer separations. 〈Ŝ2〉 is quite sensitive to the overlap between monomers; at smaller separations, the spin is equally distributed as the wavefunctions on two fragments overlap very strongly. However, as the monomers drift apart from each other, the spin changes rapidly, the S1 state shows considerable TT character with 〈Ŝ2〉 = ∼2.0 at separations >5.90 Å and behaves as a pure TT configuration at ∼6.10–6.50 Å. We will show later that this TT state is further stabilized and becomes bound with respect to free triplets upon vibronic stretching, which provides computational evidence supporting previous experimental findings for the pentacene dimer.63 The S2 state has an 〈Ŝ2〉 value of >0.75 on each fragment and is expected to retain some TT character at small separations. However, the spin drops rapidly at larger distances and this state can be primarily termed as a weak CT state (CT1). The S4 excited state remains a strong CT (CT2) state irrespective of the distance and behaves as a pure CT state at >5.90 Å (where the TT formation begins for S) with a significant doublet character (〈Ŝ2〉 = 0.75). The distinction in the extent of charge transfer is clearly evident for CT1 and CT2 states in their corresponding density difference plots in Fig. 2(b) and (c) (see Fig. S5† for natural transition orbital analysis). The moderate energetic separation of 0.23 eV and the drastic difference in character between the two CT states arise due to the asymmetric alignment in the dimer (Fig. 1).19,64 On the other hand, it has also been argued that the energy gap of “CT” states (diabatic) can be further substantially reduced from dimeric 0.8 eV to crystalline 0.04 eV due to the balance between the charge–quadrupole interaction and the crystalline dielectric environment.65,66
S3 and S5 primarily remain locally excited states termed as LE (see Section S2 in the ESI† for mathematical details). Although the overall local spin density for each of these states is zero because of their singlet characters, the spin densities on fragments can be projected. One thing of particular interest is that the spin of many electronic states except S4 drops rapidly as soon as the S1 state acquires a pure TT phenomenon, implying that the spin wavefunctions in this system interact in a way to maximize the population in the S1 state. It also suggests that the interaction between S1 and S4, i.e., TT–CT2 interaction, should be investigated more carefully as the spin manifold for only the S4 state survives during triplet pair formation. The local molecular orbitals employed in this spin analysis are visualized in Fig. S4.†
To quantify the overall intermolecular electron–electron interactions, we extract the important coulombic interactions among these SF states from multipolar moments such as dipole, quadrupole, hexapole and octupole moments using DMRG one-particle transition density matrices.67–69 The detailed mathematical expressions are included in Section S3 of the ESI.† These interactions drop rapidly as the inter-dimer separation approaches the TT formation region as predicted by the spin analysis. Most of the states weakly couple with each other, except for S1 → S4 (TT–CT2) interactions with an effective magnitude of 100–400 meV in the region of interest (Fig. 3(b)). The other CT1 state couples weakly with the TT state and the S1 → S2 couplings drop rapidly with inter-dimer separation. Other couplings of the ground and LE states with the TT state are much smaller in magnitude because of the two-electron nature of these transitions, validating the previous findings that these interactions are smaller in comparison to one-electron transitions (Fig. 3(a)).1,19 It should also be noted that the multipolar moments are one-electron operators and cannot reflect the true multielectronic nature of the SF states. Our analysis reveals strong couplings, particularly for TT–CT2 interactions, as compared to previously reported values of <200 meV at the Hartree–Fock level of theory.28 Such interactions are restored by employing DMRG ansätze which enables the coulombic coupling of the multireference electronic density among a larger number of π-electrons by the inclusion of higher-order multipolar contributions. For instance, at RC–C = 6.30 Å, more than 22% of the total coulombic coupling of the S1 → S4 transition originates from hexapole based interactions between the monomers.
Although direct electron–electron coulombic couplings exist between adiabatic S1 and S2/S4 states, these coulombic couplings do not seem to be very relevant to the formation of the TT character in S1. This is clearly seen by comparing Fig. 2 and 3: for RC–C = 6.1–6.5 Å, S1 states exhibit increasingly perfect TT characters with respect to increasing distance, while S1–S4 coulombic couplings are rapidly attenuated to zero at RC–C = 6.5–6.6 Å; at shorter RC–C ≤ 5.9 Å, on the other hand, the S1–S4 coulombic couplings are significantly intensified and surprisingly destroy the TT nature by creating a strong CT character.
In the spectrum, hereafter, the in-phase and out-of-phase vibrations are distinguished with superscripts i and o, respectively. First, we analyze the Holstein vibronic signatures of the S1 (TT) state from the active tuning modes by considering the difference density of its parent S1 state. The calculated HCs are shown in Fig. 4(a), indicating large HCs depending upon the vibrational modulation in each monomer and a few significant peaks are observed at in-phase vibrations of 101i cm−1, 476i cm−1, 778o cm−1, 1204o cm−1 and 1615i cm−1 with amplitudes of 175, 491, 403, 233 and 215 meV, respectively. This significant geometric relaxation is a result of in-plane ring bending, in-plane ring breathing, out-of-plane ring breathing, out-of-plane CC stretching and in-plane CC stretching motions. In principle, the VCs for out-of-phase modes should vanish in a perfect TT state as previously reported in a slip-stacked tetracene dimer by Ito et al.73 However, the pentacene units in our dimer are not parallel and tilted to an angle resembling more the crystal structure of pentacene in solid films. A recent study by Schnedermann et al.74 observed strong vibrational coherence in the frequency region >1000 cm−1 for the intramolecular pentacene SF system, with a perfect intensity match for the 1207 cm−1 mode. Interestingly, our calculations reveal a near-perfect energy matching for both in- and out-of-phase modes at 1204 cm−1 with HCs of 123 and 233 meV, respectively, which can effectively compensate for the energy difference [E(TT) – E(LE)] of 100–200 meV in pentacene or tetracene.
In a locally excited state S3, the majority of similar vibronic patterns originates from in-phase vibronic interactions, validating the non-CT nature of these electronic states (Fig. 4(b)). A few significant HCs for the S4 state are observed at 476 cm−1, 1204 cm−1 and 1615 cm−1, similar frequencies to the S1 state. The large HCs in S4 are attributed to its strong CT nature modulated by evident interactions through out-of-phase vibronic motion. The strong couplings of 356 meV, 1170 meV and 788 meV originate from very large out-of-plane interactions at 101o cm−1, 778o cm−1 and 1204o cm−1, respectively, and appears to activate the very crucial CT-channel (Fig. 4(c)).
The PC spectra for the significant transitions in the dimer are shown in Fig. 4(d)–(f), which relates to the vibronic modulation from the coupling modes between electronic states, i.e., S1 → S3 (TT–LE), S3 → S4 (LE–CT2) and S1 → S4 (TT–CT2). Our calculations establish the importance of PCs (∼100–300 meV) for the pentacene dimer, which is linked to the small energy differences of these low-lying states derived from DMRG/NEVPT2 computation, as well as the herringbone geometry. This is in sharp contrast to negligible PCs (<50 meV) for a model tetracene as compared to HCs in previous vibronic exciton modeling.73,75 We realize that the importance of the Peierls-type interaction is highly controversial and sensitive to the dimer packing geometry and crystalline bulk environment.35 The importance of PCs associated with an intermolecular mode was discovered early for the pentacene dimer by Zimmerman et al.,18 and an intramolecular mode of the covalent dimer was also proposed to facilitate SF via a direct coupling mechanism.76 In some other simulations, PCs were entirely neglected for pentacene,77,78 and also found to be much weaker than HCs by at least an order of magnitude in a recent DFT/MM study in which more surrounding pentacene units were included.21 Although the role of PCs appears debatable due to model approximations,31,37 there exists strong experimental evidence pointing to the significance of symmetry-breaking intermolecular modes leading to considerable couplings between singlet and TT manifolds.32,79,80
The PCs for the S1 → S3 transition have large peaks for many out-of-plane modes at 101i cm−1, 476i cm−1, 778i/o cm−1, 1204i cm−1 and 1615i/o cm−1 with amplitudes of 0.05–0.2 eV. The PCs for S3 → S4 and S1 → S4 transitions have strong PCs at similar vibrational modes; however, the spectra in these cases are dominated by several out-of-phase couplings as shown in Fig. 4(e) and (f). In all the vibronic spectra involving the strong CT2 state, we observe a characteristic peak at 1332 cm−1 for a CC ring stretching vibrational mode; indeed, it has been shown that crystalline pentacene has a similar fundamental vibrational frequency at 1380 cm−1 contributing towards the VC,81 in the condensed phase. A VC density (VCD) analysis is performed to understand the microscopic origin or VCs in Section S5† and the important normal modes are sketched in Fig. S11 of the ESI.† The other significant HCs and PCs are plotted in Fig. S8 and S9 of the ESI,† respectively.
According to the theoretical analysis by Feng et al.,56 a multiexcitonic state can be populated via either nonadiabatic transition from a source state or by coherently creating the population in the target state, both of which are governed by nonadiabatic couplings. While nonadiabatic transition is commonly promoted at the geometry where two adiabatic states intersect, coherent population can take place in any region where nonadiabatic couplings are strong. Hence the vibrational modes for promoting nonadiabatic transition may be different from the modes strengthening nonadiabatic couplings. In our vibronic coupling analysis, these modes associated with intensive PCs are therefore expected to play a major role in coherently admixing the TT configuration with CT1 and CT2 states, and these modes with intensive HCs are much more meaningful in helping tune the energies of crossing states than well separated states, depending on the symmetries of vibrational modes and the states involved (see Table S4† for the herringbone pentacene dimer).
The population transfer is evident in the composition of final SF states shown in Fig. 6, as a result of interstate VCs. At Q = 0.10, the population of the LE basis state in S′1 drops (Fig. 6(a)), while the contribution from the CT2 basis state rises rapidly with a concomitant rise of TT in S′4 (Fig. 6(c)). Consequently, a significant TT population is observed not only in S′1 but in S′4 as well. In contrast to this, the CT1 contribution to S′2 plays a crucial role when the dimer is compressed along with the vibronic modulation. The TT basis state destabilizes as a result of the crossing with the CT1 state at about Q = −0.11 and LE state at Q = −0.05 (Fig. 5) and a rather strong mixing among electronic states becomes evident at these displaced geometries (Fig. 6(a)). At Q = −0.11, the overall TT character is mainly distributed among S′1 (46%), S′2 (46%) and S′3 (5%), in which TT, CT1 and LE basis states are admixed (Fig. 6). This paves two indirect pathways for SF in the regions of vibrational stretching and compression, with the former being energetically favorable. The S′5 does not appear to play any significant role in SF in any region.
Upon the vibronic stretching, the TT of the S1 state stabilizes as E(S′1) = ∼1.71 eV at Q = +0.10 and E(S′1) = ∼1.61 eV at Q = +0.20, which becomes energetically bound from a free triplet pair (2 E(T1) = 2 × 0.83 eV of the pentacene monomer82). The TT stabilization also accounts for the evident exothermicity of ∼110 meV in the formation of the correlated triplet pair in pentacene.83 As shown in Fig. 5 and 8, this stabilization of the TT state is a result of vibronic mixing of LE (in S′3) and CT2 (in S′4) configurations, which has also been observed for pentacene and other polyenes.63,84,85 At Q = +0.10, the CT2 and LE states cross each other (Fig. 5 and 6(c)), indicating a rapid population transfer between LE and CT2 manifolds by nonadiabatic transition.
Miyata et al. also observed a similar phenomenon of TT stabilization (destabilization) in the corresponding regions of vibrational stretching (compression) for the tetracene dimer and related it to an interpolated reaction coordinate (QSF), which is dominated by the presence of high-frequency CC stretching modes.80 Interestingly, the S′1 state does not present any possibility of state crossings (or avoided crossings) with other states in the stretching region, indicating that the coherent channel involving S′1 alone dominates the SF pathway of TT radiationless conversion. As a result, at Q = +0.10, the two states, S′1 and S′4, retain a significant TT population of ∼60% and ∼40%, respectively. This indicates that the strong coherent PC interaction admixing the high-lying CT2 and TT characters leads to a multimode vibronic relaxation in S′1 (Fig. 5) by the stretching modes, towards forming a bound triplet pair state. Our computational evidence implies that the TT dissociation channel can be tuned and facilitated independently from intermolecular interactions.86
Although our dimer essentially incorporates the most significant electronic and structural features for SF, yet the mixing between electronic states in crystalline pentacene can be different owing to other solid state effects and different dielectric environments. Keeping this in mind, we adapt a similar approach to construct an effective e–ph interaction Hamiltonian for the pentacene crystal using the GW/BSE method (see Section S7 of the ESI† for computational details). The e–ph interactions in the pentacene crystal further stabilize the S1 state to ∼1.6 eV (Table S5†), so that this state becomes near-degenerate with the 2 × E(T1) for pentacene, similar to the phenomenon observed in the vibronic stretching regions of the pentacene dimer to reduce the binding energy of a TT pair along with the vibrational modulation. Although there is an evident difference in how different SF basis states mix with each other between the dimer and the crystal (see Fig. S12 of the ESI† for SF mixing in the pentacene crystal), it is straightforward to identify S′1 and high-lying S* states with a considerable TT population, provided the fact that the basis S1 state bears a TT pair in both the dimer and crystal environments. Such an S* state has been interpreted as a geometrically twisted conformation of the TT state in carotenoids.87 Kobori et al. highlighted the conformational changes in the covalent pentacene dimer as a result of the vibronic interactions induced by low-frequency modes.88
Fig. 7 The exciton sizes for SF states, representing the delocalization of the S′1 state caused by out-of-phase vibronic modulations. |
The prevailing intermediate mechanism mediating TT generation is not justified in the pentacene dimer. The TT character merely results from the superposition of adiabatic electronic configurations associated with double electron excitations, without assuming an intermediate state model. The direct electron–electron coulombic interaction exists between adiabatic S1 (TT) and S2 (weak CT)/S4 (strong CT) states. However, coulombic interaction is not found to be strongly relevant to the formation of TT pair state S1, since the S1 state of perfect TT character also appears in the pentacene dimer at a larger range of intermolecular distances (6.3–6.6 Å) where interstate coulombic interaction is rapidly attenuated to zero, while at short intermolecular distances where coulombic interaction is intensified, the singlet TT nature is annihilated with a strong charge-transfer character created.
To reveal the vibronic details of the SF process, we further include vibronic interactions in our multistate–multimode Hamiltonian, revealing that there are significant indirect couplings between the electronic states, as compared to the limited direct coupling involving states of local excitation. Through various wavefunction analyses, it is found that the initial bright singlet state exhibits the character of low-lying weak CT1 which however does not favor strong SF through the CT-mediated coherent pathway alone. Both CT1 and CT2 states must be effectively engaged in SF through a vibronic state admixture with TT pair states, causing significant TT population transfer to CT1 and CT2 states in the vibrational region of compression and stretching, respectively. The vibronic compression makes local TT characters photophysically accessible by mixing with the lower-lying weak CT1, and the vibronic stretching facilitates substantial delocalization and stabilization of TT pair states in an effective mixture with the higher-lying strong CT2. Moreover, the presence of strong VCs with some low-frequency vibrational modes provides qualitative evidence for a thermally activated SF process. The vibrational motion can be tuned by structural modification to achieve enhanced admixture between electronic states for desirable TT character and exciton energies. A comparison between the pentacene dimer and crystal environments highlights the structural features and provides a link between the crystal and vibronically tuned dimer structure. These findings brick a multistate–multimode landscape for understanding the fundamental and complex SF process in pentacene.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1sc01703a |
This journal is © The Royal Society of Chemistry 2021 |