Next Article in Journal
The Use of Flavylium Salts as Dynamic Inhibitor Moieties for Human Cb5R
Next Article in Special Issue
Manipulation of Morphology, Particle Size of Barium Sulfate and the Interacting Mechanism of Methyl Glycine Diacetic Acid
Previous Article in Journal
Enantioseparation of 3-Hydroxycarboxylic Acids via Diastereomeric Salt Formation by 2-Amino-1,2-diphenylethanol (ADPE) and Cinchonidine
Previous Article in Special Issue
Some Theoretical and Experimental Evidence for Particularities of the Siloxane Bond
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intermolecular Interactions and Charge Resonance Contributions to Triplet and Singlet Exciton States of Oligoacene Aggregates

1
Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, 40126 Bologna, Italy
2
Donostia International Physics Center (DIPC), 20018 Donostia-San Sebastian, Euskadi, Spain
3
Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU), Manuel Lardizabal 3, 20018 Donostia-San Sebastian, Euskadi, Spain
4
IKERBASQUE, Basque Foundation for Science, 48009 Bilbao, Euskadi, Spain
5
INSTM UdR Bologna, 40126 Bologna, Italy
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(1), 119; https://doi.org/10.3390/molecules28010119
Submission received: 12 November 2022 / Revised: 19 December 2022 / Accepted: 20 December 2022 / Published: 23 December 2022
(This article belongs to the Special Issue Covalent and Noncovalent Interactions in Crystal Chemistry)

Abstract

:
Intermolecular interactions modulate the electro-optical properties of molecular materials and the nature of low-lying exciton states. Molecular materials composed by oligoacenes are extensively investigated for their semiconducting and optoelectronic properties. Here, we analyze the exciton states derived from time-dependent density functional theory (TDDFT) calculations for two oligoacene model aggregates: naphthalene and anthracene dimers. To unravel the role of inter-molecular interactions, a set of diabatic states is selected, chosen to coincide with local (LE) and charge-transfer (CT) excitations within a restricted orbital space including two occupied and two unoccupied orbitals for each molecular monomer. We study energy profiles and disentangle inter-state couplings to disclose the (CT) character of singlet and triplet exciton states and assess the influence of inter-molecular orientation by displacing one molecule with respect to the other along the longitudinal translation coordinate. The analysis shows that (CT) contributions are relevant, although comparably less effective for triplet excitons, and induce a non-negligible mixed character to the low-lying exciton states for eclipsed monomers and for small translational displacements. Such (CT) contributions govern the La/Lb state inversion occurring for the low-lying singlet exciton states of naphthalene dimer and contribute to the switch from H- to J-aggregate type of the strongly allowed Bb transition of both oligoacene aggregates.

1. Introduction

Electronic and electro-optical materials based on organic molecules have been the subject of numerous investigations ranging from basic materials science to possible technological applications [1,2,3]. Oligoacenes have attracted great interest for energy and charge transport and have emerged as model organic systems for low cost, flexible, large-scale optoelectronic devices [4,5]. The properties of such devices strongly depend on the photophysical behavior of the organic molecular material, ultimately governed by the nature of their low-lying exciton states [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22]. Thus, to exploit the electro-optical properties of organic semiconductors and to tune device performances, it is desirable to fully understand the character of their electronic excitations.
Although the assessment of the nature of singlet exciton states in crystals and aggregates of organic chromophores has been the subject of several studies [23,24,25,26,27,28,29,30,31,32,33,34], the low-lying triplet exciton states have been addressed in comparably fewer investigations [35,36,37,38,39,40,41,42,43], even though triplet excitons play a fundamental role in several photoinduced processes. For instance, they are often responsible for detrimental non-radiative charge recombination in organic photovoltaic devices [44] or can be relevant for application in photocatalysis and photodynamic therapy [45,46,47].
The exciton states of molecular materials are superpositions of local (intra-molecular) excitations (LEs) and charge transfer (CT) (i.e., inter-molecular) excitations. Such intra- and inter-molecular excitations can be considered as a suitable basis of diabatic states describing electron promotions between occupied molecular orbitals to unoccupied molecular orbitals of the same (neutral) or neighboring (ionic) molecules (or sites), respectively [26,34,36,48,49,50,51,52,53,54,55,56,57,58,59].
Numerous investigations have underscored the relevance of CT states in several photoinduced processes: for example, they favor intersystem crossing through spin orbit coupling [60] or mediate singlet fission [61,62,63,64,65,66,67]. CT states have also a crucial role in the formation of excimers [37,39,42] and in exciton-dissociation and charge-separation in hetero- and homo-junctions [68,69,70].
To analyze the character of excitonic states predicted by quantum-chemical (QC) calculations, a diabatization procedure can be used to determine the superposition of LE and CT diabatic states in each adiabatic exciton state [26,29,36,55,57]. Recently, a simple diabatization approach has been developed and applied to disentangle the nature of exciton states of perylene di-imide (PDI) aggregates computed with time dependent density functional theory (TDDFT) [30,38].
In the present work, we seek to provide a better understanding of the photophysical properties of two oligoacene aggregates, naphthalene and anthracene, and employ the same diabatization procedure to characterize singlet and triplet exciton states in terms of LE and CT contributions. The approach is however extended to include a larger orbital space, required to correctly describe the orbital nature of the low-lying oligoacene excited states. Oligoacene dimers in their triplet and singlet exciton states are perfect model systems for the investigation and understanding of triplet-triplet and singlet-singlet interactions between π-conjugated molecules, and how these interactions are modulated by different intermolecular organizations. Specifically, here we focus on the longitudinal translation coordinate (Figure 1).
The paper is organized as follows: first we discuss the selection of the orbital space, mandatory to correctly define the diabatic basis for the lowest energy exciton states of oligoacene aggregates. Within the chosen orbital space, we perform a diabatization of singlet and triplet exciton states and discuss their nature, with specific attention to the magnitude of CT contributions and, for the L E components, to the parentage with isolated molecule low-lying excited states. The results of our protocol are compared with those obtained through the decomposition of the one-electron transition density matrix performed with the TheoDORE package [71]. Finally, the modulation of the absorption spectrum of oligoacene aggregates along the longitudinal translation coordinate is discussed.

2. Results and Discussion

2.1. Orbital Space and Diabatic States for the Analysis of Exciton States

The most relevant frontier molecular orbitals (MOs) of naphthalene and anthracene are collected in Figure 2, while the lowest excited states of triplet and singlet spin multiplicity (computed with TDDFT calculations in the Tamm-Dancoff approximation (TDA) [72]) for the isolated monomers, are collected in Tables S1 and S2. As it is well known according to the perimeter model [73], the lowest lying excited states of oligoacenes [73,74,75] are determined by excitations within an orbital space that includes at least two occupied (HOMO and HOMO-1) and two unoccupied (LUMO and LUMO+1) MOs leading to four excited states known as the L a ,   L b ,   B a , B b in Platt’s notation [73]. Tables S1 and S2 show that the chosen functional (ωB97X-D, see Section 3) correctly predicts the S 1 state of naphthalene to be the L b state ( B 2 u symmetry, dominated by a combination of the HOMO-1→LUMO and HOMO→LUMO+1 excitations), while in the case of anthracene the lowest state is L a ( B 1 u symmetry, dominated by the HOMO→LUMO excitation). Computed excitation energies are generally overestimated, as is typical of the long-range corrected functional employed and of the TDA [76,77].
Notably, in between the four L a ,   L b ,   B a ,   B b excited states, there are states originated by excitations also involving the HOMO-2 and the LUMO+2 orbitals. However, these states are relatively high in energy and their omission in the following discussion of exciton states will not affect the analysis of intermolecular interactions leading to the lowest exciton states of the two oligoacenes. Accordingly, the orbital space selected in the following analysis of exciton states includes two unoccupied and two occupied MOs for each monomer and, consequently, four occupied and four unoccupied orbitals for their dimers. Due to inter-molecular interactions, the energy levels associated to these dimer orbitals show distinctive oscillations along the longitudinal translation coordinate (Figures S1 and S2). The number of spin-adapted singly excited configurations (and states) originated by an orbital space of n doubly occupied and n unoccupied orbitals is n 2 , such that sixteen exciton states have been considered for each spin multiplicity and their wavefunctions expressed in terms of diabatic states, with the diabatization procedure outlined in Section 3.
The protocol employed to analyze the intermolecular interactions leading ultimately to the character of exciton states, transforms the amplitudes of TDDFT calculations carried out on the aggregate, from the basis of single excitations between aggregate’s orbitals (the delocalized excitation (DE) basis) to the basis of single excitations between molecular (monomers A and B ) site orbitals. Localized orbitals on monomers A and B allow to define diabatic states of four different types as follows:
L E A ( 1 ) = | 1 A = ( H 1 ) A ( L + 1 ) A L E B ( 1 ) = | 1 B = ( H 1 ) B ( L + 1 ) B C T A B ( 1 ) = | 1 A B = ( H 1 ) A ( L + 1 ) B C T B A ( 1 ) = | 1 B A = ( H 1 ) B ( L + 1 ) A L E A ( 2 ) = | 2 A = ( H 1 ) A ( L ) A L E B ( 2 ) = | 2 B = ( H 1 ) B ( L ) B C T A B ( 2 ) = | 2 A B = ( H 1 ) A ( L ) B C T B A ( 2 ) = | 2 B A = ( H 1 ) B ( L ) A L E A ( 3 ) = | 3 A = ( H ) A ( L ) A L E B ( 3 ) = | 3 B = ( H ) B ( L ) B C T A B ( 3 ) = | 3 A B = ( H ) A ( L ) B C T B A ( 3 ) = | 3 B A = ( H ) B ( L ) A L E A ( 4 ) = | 4 A = ( H ) A ( L + 1 ) A L E B ( 4 ) = | 4 B = ( H ) B ( L + 1 ) B C T A B ( 4 ) = | 4 A B = ( H ) A ( L + 1 ) B C T B A ( 4 ) = | 4 B A = ( H ) B ( L + 1 ) A
where HOMO and LUMO are abbreviated as H and L , subscripts A and B indicate different monomers, and each type of excitation is differentiated by a superscript.
The diabatization procedure discussed in Section 3 provides the Hamiltonian matrix in the diabatic basis representation, H d i a . In this work, compared to previous investigations in which only the HOMO and LUMO orbitals were included for each monomer, [29,30,36,38] the H d i a matrix (dimension 16 × 16) includes not only interactions between diabatic states of the same type ( L E A ( n ) , L E B ( n ) , C T A B ( n ) , C T B A ( n ) ), such as excitonic interactions V e ( n ) and super-exchange interactions [78] D e / h ( n ) (Table S3), but also interactions between diabatic states of different type ( L E / C T ( n ) , L E / C T ( p ) ), hereafter labeled with two superscript numbers, e.g., V e ( n , p ) (Table S4). In addition, intramolecular interactions between diabatic states localized on the same monomer ( L E A ( n ) , L E A ( p ) ) are also uncovered by the diabatization procedure ( H ( n , p ) in Table S4).
For aggregates characterized by a symmetric arrangement of chromophores, as is the case for those investigated here, when the molecules approach each other, intermolecular interactions mix the L E states to form a symmetry adapted (SA) superposition of (neutral) L E states, that is, Frenkel excitons ( F E ). Similarly, C T states form delocalized charge resonance ( C R ) states [48,79,80] of appropriate symmetry. The symmetry point group of an oligoacene aggregate, when intermolecular displacements along the longitudinal translation coordinate are considered, (Figure 1) is C 2 h . As a result, the most relevant π π * exciton states along with F E and C R diabatic states, all belong to A g , A u ,   B g and B u symmetry representations. In the following Figures, each symmetry will be distinguished by a specific color code used throughout this work: yellow for A g , orange for A u , blue for B g and green for B u .
The SA diabatic states are then obtained as linear combinations of L E and C T states as shown in Figure 3 and the corresponding diabatic matrices H d i a S A are collected in Tables S5 and S6. The energy profiles of the diabatic states and the most relevant interactions extracted from the diabatization procedure are collected in Figures S3–S9.
In the following sections we analyze the TDA computed singlet and triplet exciton states of oligoacene dimers to determine their character ( C T / L E ) and disentangle the role of interactions between diabatic states along the longitudinal translation coordinate.

2.2. Singlet Exciton States of Naphthalene and Anthracene Dimers

The TDA computed excitation energy profiles of the singlet exciton states of naphthalene and anthracene dimers are collected in Figure 4 and Figure S10. We note that B 1 u ( L a ) excited states of the isolated molecules originate two L E states in the dimer, belonging to A u and B g symmetries of the C 2 h point group, while B 2 u ( L b ) monomer excited states originate two L E states in the dimer belonging to A g and B u symmetries. Because the lowest excited states of naphthalene and anthracene are, respectively, L b and L a , one might expect the lowest exciton state with different nature and symmetry for the two aggregates. In contrast, both naphthalene and anthracene dimers show that, for small displacements from the eclipsed geometry, the lowest energy exciton state belongs to B g symmetry, which therefore originates from the L a ( B 1 u ) monomer state for both dimers. This result is expected for anthracene while it implies an inversion of the L b / L a states for naphthalene, when moving from the monomer to the aggregate. Such an inversion has been documented in previous studies at the equilibrium intermolecular distance of the singlet state excimer [81,82], which is considerably shorter than the intermolecular distance of 3.4 Å considered here. Interestingly, the stabilization of the lowest B g state leading to the formation of the naphthalene excimer is predicted by our TDA-ωB97X-D calculations even at such large intermolecular distances.
Although the 1 B g exciton state of the oligoacene dimers derives from the L a state and therefore owns a considerable L E contribution, its character is indeed mixed, showing more that 40% C T contribution (Figure 4c,d) at the eclipsed geometry. This is true not only for the 1 B g state, but also for other low-lying exciton states, with the C T contribution slightly increasing from naphthalene to anthracene dimers. Such a mixed character can be appreciated by graphically representing the wavefunction of the two lowest exciton states (1 B g and 1 A g ) in terms of their SA diabatic states. The contributions of C R states (red and brown lines in Figure 5) emerge clearly not only for the eclipsed dimer configuration, but up to about 4 Å longitudinal displacements, even though the weight of F E states (green and blue lines) is dominant.
As shown in previous investigations [29,30,36,38], the modulation of adiabatic energy profiles along the longitudinal translation coordinate can be rationalized in terms of inter-state interactions between SA diabatic states. The different effects of inter-state and intermolecular interactions can be appreciated by comparing SA diabatic and resulting adiabatic energy profiles of A g (Figure 6) and B g (Figure 7) symmetry. In both cases, the adiabatic exciton energy profiles result from the combination of i) the interaction between F E and C R SA diabatic states (grey lines in panels (c,d)), which is maximum at the eclipsed geometry and oscillates along the translation, with ii) the interaction between the two F E SA diabatic states (dark-turquoise lines in panels (c,d)), less dramatically changing along the translational coordinate. Such an interaction is larger for A g symmetry states and more effective, given the quasi degeneracy of F E ( 2 ) and F E ( 4 ) energy profiles, while it is far less effective for B g states due to the large energy difference between F E ( 1 ) and F E ( 3 ) states. This explains the remarkable energy lowering of the 1 A g adiabatic exciton state (yellow squares in Figure 6a,b) compared to the SA diabatic states (green) for all translational displacements, which is not observed for B g states except for specific displacement ranges corresponding to large F E / C R interactions. Similar considerations justify the excitation energy profiles of adiabatic exciton states belonging to the remaining two symmetries (Figures S11 and S12).

2.3. Triplet Exciton States of Naphthalene and Anthracene Dimers

Figure 8 and Figure S10 collect the excitation energy profiles of triplet excitons of oligoacene aggregates. In contrast with singlet excitons, we note that the lowest two triplet exciton states belong to B g and A u symmetries, while the lowest A g exciton state is found at much higher energy. Similar to singlet excitons, a few low-lying states (dark-blue/dark-red in Figure 8) were not included in the diabatization procedure, since their influence is negligible due to their different orbital nature and energy separation from the lowest exciton states. The C T character analysis (Figure 8c,d) demonstrates that both the 1 B g and 1 A g states have a mixed C T / L E character with the C T contribution slightly increasing for the longer acene dimer.
The significant C T contribution to 1 B g and 1 A g states, contrasts with the almost negligible contribution to 1 B u and 1 A u states. Such differences can be rationalized by comparing, for instance, the SA diabatic energy profiles ( B g and A u symmetry) with F E and C R character (green and red curves, respectively, in Figure 9 and Figure 10) and their inter-state interactions, depicted in the bottom part of the above figures. Similar considerations hold for the other symmetry species (Figures S13 and S14). The F E / C R energy differences are very similar for the two sets of SA diabatic states, but their couplings (grey lines) are much larger for the B g symmetry states. Specifically, the largest interaction amounts to more than 0.8 eV for B g and it does not exceed 0.3 eV for A u states. As a result, 1 B g adiabatic states display a non-negligible C T character for small longitudinal displacements, while 1 A u states do not.
The analysis of interactions between SA diabatic states uncovers also specific differences between singlet and triplet excitons. One well-known distinctive element, previously documented for PDI aggregates [36,38], is the larger energy separation between C R and F E triplet diabatic states compared to singlet states, accounting for the reduced effect of super-exchange interactions, ultimately leading to a less marked C T character of lowest-lying triplet exciton states. An additional distinction can be appreciated in the case of oligoacenes and concerns the couplings (Tables S5 and S6) between SA diabatic states of F E type (dark-turquoise lines in Figure 6 and Figure 7 for singlet and Figure 9 and Figure 10 for triplet excitons). These are larger than 0.8 eV for singlet SA states of A g / B u symmetry while they do not exceed 0.35 eV for the triplet SA states of any symmetry. As a result, the adiabatic triplet exciton states of all symmetry species almost overlap with the SA diabatic states for large portions of the energy profiles, in contrast with those of singlet spin multiplicity discussed in the previous section. These interactions result from the combination of an intra-molecular contribution (top part of Figures S8 and S9) with the inter-molecular exciton coupling (bottom part of Figures S8 and S9) both of which are much larger for singlet states.

2.4. CT Character of Singlet and Triplet Excitons of Oligoacenes

The C T character analysis obtained from the above discussed diabatization procedure (Figure 4c,d and Figure 8c,d) can be compared with the results of other analysis tools. To visualize the nature of the lowest singlet and triplet exciton states, for the eclipsed dimer configurations of the two oligoacene dimers, we carried out a fragment-based analysis via electron–hole correlation plots, using TheoDORE [71] (Figure 11). The two selected fragments correspond to the two molecules forming the dimer. Exciton states are identified by the non-vanishing elements of the 2 × 2 matrix (the Ω-matrix [71]) represented by different levels of grey. Locally excited contributions appear in Figure 11 off-diagonally (going from lower left to upper right), while C T contributions appear on the main diagonal. In agreement with the character analysis shown in Figure 4c,d and Figure 8c,d, also Figure 11 shows that for the eclipsed configuration, the character of the lowest singlet and triplet exciton states is mixed L E / C T as indicated by the light-grey squares in the electron-hole correlation plot. As discussed in previous sections, the C T contribution is larger for singlet exciton states than for triplet exciton states (main diagonal squares for the singlet excitons are dark grey while those of triplet exciton states are light grey) as a result of the larger energy separation between C T and L E diabatic triplet states (compare Figures S3 and S5 or Figures S4 and S6). Overall, both approaches provide the same results for the lowest triplet and singlet exciton states as regard the C T contributions, although our diabatization analysis also provides detailed information on relevant intermolecular and interstate interactions.

2.5. Absorption Spectrum and H-/J- Character Switch along the Longitudinal Translation Coordinate

In previous sections we focused the attention on low-lying exciton states of the oligoacene dimers. Because the absorption spectrum of oligoacenes is dominated by the transition to the higher energy B 2 u state (the B b state in Platt’s notation), it is interesting to discuss the evolution (along the displacement coordinate) of the dipole allowed exciton state bearing the largest parentage with the B b state and how its excitation energy modulation influences the appearance of the absorption spectrum. The computed absorption spectrum for increasingly large translational displacements (Figure 12) shows that the intense absorption occurs at higher energies than the isolated molecule, suggesting an H- type aggregate behavior for small displacements from the eclipsed configurations. For larger displacements, the most intense absorption band moves abruptly to lower excitation energies (J- type) after a narrow intermediate region in which two bands of similar intensity appear. The evolution of the excitation energy of such strongly allowed exciton state is shown in black in Figure S15. This state corresponds to the 4 B u state for both naphthalene and anthracene at the eclipsed geometry and for displacements up to 4 or 5 Å, then it switches to the 2 B u state. Such a switch is associated with a sudden change in the exciton state character, acquiring a dominant C T contribution for larger displacements (Figure S16), suggesting that the H- to J- type transition along the longitudinal displacement is assisted not only by the sign change of the exciton interaction (Figure S7) but also by the interaction between F E and C R states.
Finally, for large displacements (8 Å) the strongest absorption peak has moved back almost toward the isolated molecule value, although not completely owing to some still non-negligible exciton interactions. Indeed the V e ( 2 ) and V e ( 4 ) couplings are still not vanishing at such distance (Figure S7), in contrast with the V e ( 1 ) and V e ( 3 ) interactions that have already become negligible.

3. Computational Models

The ground state monomer structure of naphthalene and anthracene was optimized at the ωB97X-D/6-31G* level of theory. The distance between the planes of different monomers was set to 3.4 Å, as used in previous investigations on dimers of other polycyclic aromatic hydrocarbons [29,30]. Exciton states were computed for the eclipsed aggregates and for displacements of 0.5 Å up to 8.0 Å, along the longitudinal translation coordinate (x) (Figure 1). Excitation energies were determined with TDDFT calculations with the TDA [72], using the ωB97X-D functional [83], previously shown to provide a reliable description of C T character in singlet excitons of PDI dimers [30,34,38], and the 6-31G* basis set. All QC calculations were carried out with the Gaussian16 suite of programs [84].
To analyze the exciton character, we followed the approach described in previous works [30,38] and we expressed each relevant exciton state in terms of L E s . To this end, we selected the orbital subspace corresponding to relevant π π * exciton states. As discussed in Section 2, for oligoacene aggregates this must include at least the HOMO/HOMO-1 and LUMO/LUMO+1 of each monomer [27,75] and represents the minimal orbital space (MIOS) sufficient to reliably describe low-lying excited states of each monomer. Each aggregate’s orbital obtained from QC calculations is then expressed as linear combination of monomer orbitals. These linear combination coefficients C i , j A G G R _ M O B form the C A G G R _ M O B matrix describing each aggregate’s orbital in the monomer orbital basis (MOB) and are obtained as [25,30,85]:
C A G G R _ M O B = C M O N _ A O B t · S M O N _ A O B · C A G G R _ A O B    
where the C M O N _ A O B matrix is a block diagonal matrix containing the MOs coefficients in the atomic orbital basis (AOB) of each monomer, with off-diagonal blocks set to zero and S M O N _ A O B is the overlap matrix of the monomers in the AOB.
Since monomer orbitals belonging to two different molecules are non-orthogonal to each other, aggregate’s orbitals C A G G R _ M O B L expressed in terms of orthogonalized monomer orbitals are obtained as:
C A G G R _ M O B L = S A G G R _ _ M O B 1 2 · C A G G R _ M O B    
where superscript L indicates Löwdin’s orthogonalization [86], and the overlap matrix S A G G R _ M O B = C M O N _ A O B t · S A G G R _ A O B · C M O N _ A O B is obtained from the coefficients of monomer’s orbitals C M O N _ A O B and the overlap of the atomic orbitals in the aggregate configuration S A G G R _ A O B .
From the results of TDDFT calculations on the aggregate, the subset of n 2 exciton states originated from the MIOS of the aggregate (including n occupied and n unoccupied orbitals) are then selected out of the full set of computed eigenstates. TDDFT amplitudes are expressed on the basis of delocalized excitations (DEs), namely excitations between aggregate’s orbitals, and form the columns of the B D E a d i a matrix. Thus, each DE must be expanded in terms of excitations between monomer orbitals (diabatic L E and C T states). With aggregate’s orbitals expressed in terms of monomer orbitals via the C A G G R _ M O B L matrix, each DE( i j ) from an occupied i to an empty j aggregate’s orbital can be expressed as a linear combination of diabatic ( L E and C T ) excitations ( k l ) from an occupied k to an empty l monomer orbital, with expansion coefficients given by
U k l , i j D E d i a = C k , i A G G R _ M O B ,   L   · C l , j A G G R _ M O B , L  
these coefficients form the columns of the unitary matrix U D E   d i a .
Exciton states are then readily expressed in the diabatic basis as
B d i a a d i a = U D E d i a · B D E a d i a  
and the character of each exciton state is obtained by summing up the contributions from C T and L E states.
The corresponding n 2 eigenvalues (excitation energies of the selected adiabatic excitons) form the diagonal H a d i a matrix, from which the Hamiltonian in the diabatic L E / C T basis, H d i a , can be obtained as [26,54,57,87]
H d i a = B d i a a d i a · H a d i a · B d i a a d i a t
Finally, the H d i a is rotated in the SA diabatic basis formed by F E and C R states, to obtain a block diagonal matrix with four sub-matrices H d i a S A   (for B u , A g ,   B g and A u states, (Tables S5 and S6)) whose off-diagonal elements are the interactions between C R and F E states that ultimately govern the modulation of adiabatic exciton state energies along the longitudinal translation coordinate.

4. Conclusions

Intermolecular interactions determine the nature of exciton states which ultimately govern the optoelectronic properties and the outcome of photoinduced processes in molecular materials. In this work we have analyzed the character and modulation of singlet and triplet exciton states of two oligoacene homo-dimers formed by naphthalene and anthracene, along the intermolecular longitudinal translation coordinate.
The character of exciton states computed with TDDFT and the relevant inter-molecular and inter-state interactions were determined with a diabatization procedure successfully employed in previous investigations of PDI and here extended to include a larger orbital space, mandatory to correctly describe the low-lying excited states of oligoacenes.
The analysis in terms of localized excitations shows that C T contributions are relevant, although comparably less effective for triplet excitons, and induce a non-negligible mixed character to the low-lying exciton states for eclipsed monomers and for small translational displacements compatible with the formation of excimers not only for singlet but also for triplet excitons.
Concerning the singlet spin manifold, the study shows that such C T contributions drive the L a / L b state inversion of the lowest-lying exciton state of naphthalene dimer and assist the switch from H- to J-aggregate type of the strongly allowed B b transition of both oligoacene aggregates.
We believe that this study provides useful insights on the magnitude of inter-molecular interactions occurring in molecular materials and determining the nature of exciton states of both singlet and triplet spin multiplicity, paving the way to future investigations on more complex aggregates and inter-molecular organizations.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28010119/s1, Figure S1. Energy profiles of the frontier molecular orbitals of naphthalene. Figure S2. Energy profiles of the frontier molecular orbitals of anthracene. Figures S3–S6. Energy profiles of the singlet and triplet diabatic exciton states of naphthalene and anthracene. Figure S7. Energy profiles of the excitonic and super-exchange interactions. Figures S8 and S9. Energy profiles of interactions between LE diabatic singlet and triplet states. Figure S10. Adiabatic energy profiles of singlet and triplet adiabatic exciton states of naphthalene and anthracene. Figures S11–S14. Analysis of the excitation energy profiles. Figure S15. The evolution of the strongly allowed exciton state of naphthalene and anthracene. Figure S16. CT character of the singlet adiabatic exciton states of Bu symmetry of naphthalene and anthracene. Figure S17. Angle between the transition dipole moments of the monomers forming naphthalene and anthracene dimers. Table S1. Triplet excited states of naphthalene and anthracene monomers. Table S2. Singlet excited states of naphthalene and anthracene monomer. Table S3. Matrix elements of H d i a between diabatic states of the same type. Table S4. Matrix elements of H d i a between diabatic states of different types. Tables S5 and S6. Matrix elements of H d i a S A .

Author Contributions

Conceptualization, F.N.; investigation and formal analysis, F.N., D.C., Y.D., A.C. and M.Z.-U.; writing—original draft preparation, F.N.; writing—review and editing, F.N., D.C. and Y.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the Supplementary Material.

Acknowledgments

The authors acknowledge financial support from “Valutazione della Ricerca di Ateneo” (VRA)—University of Bologna. Y.D. acknowledges Ministero dell’Università e della Ricerca for her PhD fellowship.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not applicable.

References

  1. Ahmad, S. Organic semiconductors for device applications: Current trends and future prospects. J. Polym. Eng. 2014, 34, 279–338. [Google Scholar] [CrossRef]
  2. Shirota, Y.; Kageyama, H. Charge carrier transporting molecular materials and their applications in devices. Chem. Rev. 2007, 107, 953–1010. [Google Scholar] [CrossRef] [PubMed]
  3. Friend, R. Organic Materials for Large Area Electronics. Mater. Sci. Forum 2008, 608, 159–179. [Google Scholar] [CrossRef]
  4. Lee, N.E.; Zhou, J.J.; Agapito, L.A.; Bernardi, M. Charge transport in organic molecular semiconductors from first principles: The bandlike hole mobility in a naphthalene crystal. Phys. Rev. B 2018, 97, 115203. [Google Scholar] [CrossRef] [Green Version]
  5. Mohajeri, A.; Farmani, M. Toward More Efficient Organic Semiconductors: The Relationship between Morphology, Charge Transport, and Photophysical Properties. ACS Appl. Electron. Mater. 2022, 4, 246–258. [Google Scholar] [CrossRef]
  6. Ostroverkhova, O. Organic Optoelectronic Materials: Mechanisms and Applications. Chem. Rev. 2016, 116, 13279–13412. [Google Scholar] [CrossRef]
  7. Xu, Y.; Xu, P.; Hu, D.; Ma, Y. Recent progress in hot exciton materials for organic light-emitting diodes. Chem. Soc. Rev. 2021, 50, 1030–1069. [Google Scholar] [CrossRef]
  8. Yu, P.; Zhen, Y.; Dong, H.; Hu, W. Crystal Engineering of Organic Optoelectronic Materials. Chem 2019, 5, 2814–2853. [Google Scholar] [CrossRef]
  9. Würthner, F.; Saha-Möller, C.R.; Fimmel, B.; Ogi, S.; Leowanawat, P.; Schmidt, D. Perylene Bisimide Dye Assemblies as Archetype Functional Supramolecular Materials. Chem. Rev. 2016, 116, 962–1052. [Google Scholar] [CrossRef]
  10. Hecht, M.; Würthner, F. Supramolecularly Engineered J-Aggregates Based on Perylene Bisimide Dyes. Acc. Chem. Res. 2021, 54, 642–653. [Google Scholar] [CrossRef]
  11. Hoeben, F.J.M.; Jonkheijm, P.; Meijer, E.W.; Schenning, A.P.H.J. About supramolecular assemblies of pi-conjugated systems. Chem. Rev. 2005, 105, 1491–1546. [Google Scholar] [CrossRef] [PubMed]
  12. Su, F.; Chen, G.; Korevaar, P.A.; Pan, F.; Liu, H.; Guo, Z.; Schenning, A.P.H.J.; Zhang, H.; Lin, J.; Jiang, Y. Discrete π-Stacks from Self-Assembled Perylenediimide Analogues. Angew. Chem. Int. Ed. Engl. 2019, 58, 15273–15277. [Google Scholar] [CrossRef] [PubMed]
  13. Yang, L.; Langer, P.; Davies, E.S.; Baldoni, M.; Wickham, K.; Besley, N.A.; Besley, E.; Champness, N.R. Synthesis and characterisation of rylene diimide dimers using molecular handcuffs. Chem. Sci. 2019, 10, 3723–3732. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Kaufmann, C.; Bialas, D.; Stolte, M.; Würthner, F. Discrete π-Stacks of Perylene Bisimide Dyes within Folda-Dimers: Insight into Long- and Short-Range Exciton Coupling. J. Am. Chem. Soc. 2018, 140, 9986–9995. [Google Scholar] [CrossRef]
  15. Kaufmann, C.; Kim, W.; Nowak-Król, A.; Hong, Y.; Kim, D.; Würthner, F. Ultrafast Exciton Delocalization, Localization, and Excimer Formation Dynamics in a Highly Defined Perylene Bisimide Quadruple π-Stack. J. Am. Chem. Soc. 2018, 140, 4253–4258. [Google Scholar] [CrossRef] [PubMed]
  16. Samanta, S.; Chaudhuri, D. Suppressing Excimers in H-Aggregates of Perylene Bisimide Folda-Dimer: Role of Dimer Conformation and Competing Assembly Pathways. J. Phys. Chem. Lett. 2017, 8, 3427–3432. [Google Scholar] [CrossRef] [PubMed]
  17. Hong, Y.; Kim, J.; Kim, W.; Kaufmann, C.; Kim, H.; Würthner, F.; Kim, D. Efficient Multiexciton State Generation in Charge-Transfer-Coupled Perylene Bisimide Dimers via Structural Control. J. Am. Chem. Soc. 2020, 142, 7845–7857. [Google Scholar] [CrossRef]
  18. Haldar, R.; Mazel, A.; Krstić, M.; Zhang, Q.; Jakoby, M.; Howard, I.A.; Richards, B.S.; Jung, N.; Jacquemin, D.; Diring, S.; et al. A de novo strategy for predictive crystal engineering to tune excitonic coupling. Nat. Commun. 2019, 10, 2048. [Google Scholar] [CrossRef] [Green Version]
  19. Haldar, R.; Mazel, A.; Joseph, R.; Adams, M.; Howard, I.A.; Richards, B.S.; Tsotsalas, M.; Redel, E.; Diring, S.; Odobel, F.; et al. Excitonically Coupled States in Crystalline Coordination Networks. Chem. A Eur. J. 2017, 23, 14316–14322. [Google Scholar] [CrossRef]
  20. Spano, F.C. The Spectral Signatures of Frenkel Polarons in H- and J-Aggregates. Acc. Chem. Res. 2010, 43, 429–439. [Google Scholar] [CrossRef]
  21. Hestand, N.J.; Spano, F.C. Molecular Aggregate Photophysics beyond the Kasha Model: Novel Design Principles for Organic Materials. Acc. Chem. Res. 2017, 50, 341–350. [Google Scholar] [CrossRef] [PubMed]
  22. Hestand, N.J.; Spano, F.C. Expanded Theory of H- and J-Molecular Aggregates: The Effects of Vibronic Coupling and Intermolecular Charge Transfer. Chem. Rev. 2018, 118, 7069–7163. [Google Scholar] [CrossRef] [PubMed]
  23. Yamagata, H.; Norton, J.; Hontz, E.; Olivier, Y.; Beljonne, D.; Brédas, J.L.; Silbey, R.J.; Spano, F.C. The nature of singlet excitons in oligoacene molecular crystals. J. Chem. Phys. 2011, 134, 204703. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Hestand, N.J.; Yamagata, H.; Xu, B.; Sun, D.; Zhong, Y.; Harutyunyan, A.R.; Chen, G.; Dai, H.-L.; Rao, Y.; Spano, F.C. Polarized Absorption in Crystalline Pentacene: Theory vs Experiment. J. Phys. Chem. C 2015, 119, 22137–22147. [Google Scholar] [CrossRef]
  25. Kim, D. A Theoretical Analysis of the Excited State of Oligoacene Aggregates: Local Excitation vs. Charge-Transfer Transition. Bull. Korean Chem. Soc. 2015, 36, 2284–2289. [Google Scholar] [CrossRef]
  26. Shirai, S.; Iwata, S.; Tani, T.; Inagaki, S. Ab Initio Studies of Aromatic Excimers Using Multiconfiguration Quasi-Degenerate Perturbation Theory. J. Phys. Chem. A 2011, 115, 7687–7699. [Google Scholar] [CrossRef]
  27. Fliegl, H.; You, Z.; Hsu, C.; Sundholm, D. The Excitation Spectra of Naphthalene Dimers: Frenkel and Charge-Transfer Excitons. J. Chinese Chem. Soc. 2016, 63, 20–32. [Google Scholar] [CrossRef]
  28. Puschnig, P.; Ambrosch-Draxl, C. Excitons in organic semiconductors. Comptes Rendus Phys. 2009, 10, 504–513. [Google Scholar] [CrossRef]
  29. Liu, W.; Canola, S.; Köhn, A.; Engels, B.; Negri, F.; Fink, R.F. A model hamiltonian tuned toward high level ab initio calculations to describe the character of excitonic states in perylenebisimide aggregates. J. Comput. Chem. 2018, 39, 1979–1989. [Google Scholar] [CrossRef]
  30. Canola, S.; Bagnara, G.; Dai, Y.; Ricci, G.; Calzolari, A.; Negri, F. Addressing the Frenkel and charge transfer character of exciton states with a model Hamiltonian based on dimer calculations: Application to large aggregates of perylene bisimide. J. Chem. Phys. 2021, 154, 124101. [Google Scholar] [CrossRef]
  31. Liu, W.; Settels, V.; Harbach, P.H.P.; Dreuw, A.; Fink, R.F.; Engels, B. Assessment of TD-DFT- and TD-HF-based approaches for the prediction of exciton coupling parameters, potential energy curves, and electronic characters of electronically excited aggregates. J. Comput. Chem. 2011, 32, 1971–1981. [Google Scholar] [CrossRef] [PubMed]
  32. Stehr, V.; Engels, B.; Deibel, C.; Fink, R.F. Anisotropy of singlet exciton diffusion in organic semiconductor crystals from ab initio approaches. J. Chem. Phys. 2014, 140, 024503. [Google Scholar] [CrossRef] [PubMed]
  33. Engels, B.; Engel, V. The dimer-approach to characterize opto-electronic properties of and exciton trapping and diffusion in organic semiconductor aggregates and crystals. Phys. Chem. Chem. Phys. 2017, 19, 12604–12619. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Walter, C.; Krämer, V.; Engels, B. On the applicability of time-dependent density functional theory (TDDFT) and semiempirical methods to the computation of excited-state potential energy surfaces of perylene-based dye-aggregates. Int. J. Quantum Chem. 2017, 117, e25337. [Google Scholar] [CrossRef]
  35. Jortner, J.; Rice, S.A.; Katz, J.L.; Choi, S. Triplet Excitons in Crystals of Aromatic Molecules. J. Chem. Phys. 1965, 42, 309–323. [Google Scholar] [CrossRef]
  36. Zubiria-Ulacia, M.; Matxain, J.M.; Casanova, D. The role of CT excitations in PDI aggregates. Phys. Chem. Chem. Phys. 2020, 22, 15908–15918. [Google Scholar] [CrossRef]
  37. Diaz-Andres, A.; Casanova, D. Benzene Excimer and Excited Multimers: Electronic Character, Interaction Nature, and Aromaticity. J. Phys. Chem. Lett. 2021, 12, 7400–7408. [Google Scholar] [CrossRef]
  38. Dai, Y.; Zubiria-Ulacia, M.; Casanova, D.; Negri, F. Impact of Charge-Resonance Excitations on CT-Mediated J-Type Aggregation in Singlet and Triplet Exciton States of Perylene Di-Imide Aggregates: A TDDFT Investigation. Computation 2022, 10, 18. [Google Scholar] [CrossRef]
  39. Krishnan, A.; Diaz-Andres, A.; Sudhakaran, K.P.; John, A.T.; Hariharan, M.; Casanova, D. Deciphering the role of (anti)aromaticity in cofacial excimers of linear acenes. J. Phys. Org. Chem. 2022, 36, e4438. [Google Scholar] [CrossRef]
  40. Hartzler, D.A.; Slipchenko, L.V.; Savikhin, S. Triplet-Triplet Coupling in Chromophore Dimers: Theory and Experiment. J. Phys. Chem. A 2018, 122, 6713–6723. [Google Scholar] [CrossRef]
  41. You, Z.-Q.; Hsu, C.-P.; Fleming, G.R. Triplet-triplet energy-transfer coupling: Theory and calculation. J. Chem. Phys. 2006, 124, 044506. [Google Scholar] [CrossRef]
  42. East, A.L.L.; Lim, E.C. Naphthalene dimer: Electronic states, excimers, and triplet decay. J. Chem. Phys. 2000, 113, 8981–8994. [Google Scholar] [CrossRef] [Green Version]
  43. Pabst, M.; Lunkenheimer, B.; Köhn, A. The Triplet Excimer of Naphthalene: A Model System for Triplet−Triplet Interactions and Its Spectral Properties. J. Phys. Chem. C 2011, 115, 8335–8344. [Google Scholar] [CrossRef]
  44. Gillett, A.J.; Privitera, A.; Dilmurat, R.; Karki, A.; Qian, D.; Pershin, A.; Londi, G.; Myers, W.K.; Lee, J.; Yuan, J.; et al. The role of charge recombination to triplet excitons in organic solar cells. Nature 2021, 597, 666–671. [Google Scholar] [CrossRef]
  45. Cakmak, Y.; Kolemen, S.; Duman, S.; Dede, Y.; Dolen, Y.; Kilic, B.; Kostereli, Z.; Yildirim, L.T.; Dogan, A.L.; Guc, D.; et al. Designing excited states: Theory-guided access to efficient photosensitizers for photodynamic action. Angew. Chem. Int. Ed. 2011, 50, 11937–11941. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Bassan, E.; Gualandi, A.; Cozzi, P.G.; Ceroni, P. Design of BODIPY dyes as triplet photosensitizers: Electronic properties tailored for solar energy conversion, photoredox catalysis and photodynamic therapy. Chem. Sci. 2021, 12, 6607–6628. [Google Scholar] [CrossRef] [PubMed]
  47. Gibbons, D.J.; Farawar, A.; Mazzella, P.; Leroy-Lhez, S.; Williams, R.M. Making triplets from photo-generated charges: Observations, mechanisms and theory. Photochem. Photobiol. Sci. 2020, 19, 136–158. [Google Scholar] [CrossRef] [Green Version]
  48. Bardeen, C.J. The Structure and Dynamics of Molecular Excitons. Annu. Rev. Phys. Chem. 2014, 65, 127–148. [Google Scholar] [CrossRef] [PubMed]
  49. Mewes, S.A.; Dreuw, A. Density-based descriptors and exciton analyses for visualizing and understanding the electronic structure of excited states. Phys. Chem. Chem. Phys. 2019, 21, 2843–2856. [Google Scholar] [CrossRef]
  50. Mewes, S.A.; Plasser, F.; Krylov, A.; Dreuw, A. Benchmarking Excited-State Calculations Using Exciton Properties. J. Chem. Theory Comput. 2018, 14, 710–725. [Google Scholar] [CrossRef] [PubMed]
  51. Jurinovich, S.; Cupellini, L.; Guido, C.A.; Mennucci, B. EXAT: EXcitonic analysis tool. J. Comput. Chem. 2018, 39, 279–286. [Google Scholar] [CrossRef] [PubMed]
  52. Accomasso, D.; Persico, M.; Granucci, G. Diabatization by Localization in the Framework of Configuration Interaction Based on Floating Occupation Molecular Orbitals (FOMO−CI). ChemPhotoChem 2019, 3, 933–944. [Google Scholar] [CrossRef]
  53. Tamura, H. Diabatization for Time-Dependent Density Functional Theory: Exciton Transfers and Related Conical Intersections. J. Phys. Chem. A 2016, 120, 9341–9347. [Google Scholar] [CrossRef] [PubMed]
  54. Liu, W.; Lunkenheimer, B.; Settels, V.; Engels, B.; Fink, R.F.; Köhn, A. A general ansatz for constructing quasi-diabatic states in electronically excited aggregated systems. J. Chem. Phys. 2015, 143, 084106. [Google Scholar] [CrossRef]
  55. Carreras, A.; Uranga-Barandiaran, O.; Castet, F.; Casanova, D. Photophysics of Molecular Aggregates from Excited State Diabatization. J. Chem. Theory Comput. 2019, 15, 2320–2330. [Google Scholar] [CrossRef]
  56. Darghouth, A.A.M.H.M.; Correa, G.C.; Juillard, S.; Casida, M.E.; Humeniuk, A.; Mitrić, R. Davydov-type excitonic effects on the absorption spectra of parallel-stacked and herringbone aggregates of pentacene: Time-dependent density-functional theory and time-dependent density-functional tight binding. J. Chem. Phys. 2018, 149, 134111. [Google Scholar] [CrossRef] [Green Version]
  57. Casanova, D. Theoretical investigations of the perylene electronic structure: Monomer, dimers, and excimers. Int. J. Quantum Chem. 2015, 115, 442–452. [Google Scholar] [CrossRef]
  58. Mao, Y.; Montoya-Castillo, A.; Markland, T.E. Accurate and efficient DFT-based diabatization for hole and electron transfer using absolutely localized molecular orbitals. J. Chem. Phys. 2019, 151, 164114. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Mao, Y.; Montoya-Castillo, A.; Markland, T.E. Excited state diabatization on the cheap using DFT: Photoinduced electron and hole transfer. J. Chem. Phys. 2020, 153, 244111. [Google Scholar] [CrossRef] [PubMed]
  60. Lv, M.; Lu, X.; Jiang, Y.; Sandoval-Salinas, M.E.; Casanova, D.; Sun, H.; Sun, Z.; Xu, J.; Yang, Y.; Chen, J. Near-Unity Triplet Generation Promoted via Spiro-Conjugation. Angew. Chem. Int. Ed. 2021, 61, e202113190. [Google Scholar] [CrossRef]
  61. Margulies, E.A.; Logsdon, J.L.; Miller, C.E.; Ma, L.; Simonoff, E.; Young, R.M.; Schatz, G.C.; Wasielewski, M.R. Direct Observation of a Charge-Transfer State Preceding High-Yield Singlet Fission in Terrylenediimide Thin Films. J. Am. Chem. Soc. 2017, 139, 663–671. [Google Scholar] [CrossRef] [PubMed]
  62. Margulies, E.A.; Miller, C.E.; Wu, Y.; Ma, L.; Schatz, G.C.; Young, R.M.; Wasielewski, M.R. Enabling singlet fission by controlling intramolecular charge transfer in π-stacked covalent terrylenediimide dimers. Nat. Chem. 2016, 8, 1120–1125. [Google Scholar] [CrossRef]
  63. Young, R.M.; Wasielewski, M.R. Mixed Electronic States in Molecular Dimers: Connecting Singlet Fission, Excimer Formation, and Symmetry-Breaking Charge Transfer. Acc. Chem. Res. 2020, 53, 1957–1968. [Google Scholar] [CrossRef] [PubMed]
  64. Beljonne, D.; Yamagata, H.; Brédas, J.L.; Spano, F.C.; Olivier, Y. Charge-Transfer Excitations Steer the Davydov Splitting and Mediate Singlet Exciton Fission in Pentacene. Phys. Rev. Lett. 2013, 110, 226402. [Google Scholar] [CrossRef] [Green Version]
  65. Berkelbach, T.C.; Hybertsen, M.S.; Reichman, D.R. Microscopic theory of singlet exciton fission. II. Application to pentacene dimers and the role of superexchange. J. Chem. Phys. 2013, 138, 114103. [Google Scholar] [CrossRef] [Green Version]
  66. Zeng, T.; Hoffmann, R.; Ananth, N. The Low-Lying Electronic States of Pentacene and Their Roles in Singlet Fission. J. Am. Chem. Soc. 2014, 136, 5755–5764. [Google Scholar] [CrossRef] [PubMed]
  67. Casanova, D. Theoretical Modeling of Singlet Fission. Chem. Rev. 2018, 118, 7164–7207. [Google Scholar] [CrossRef] [PubMed]
  68. Han, G.; Yi, Y. Local Excitation/Charge-Transfer Hybridization Simultaneously Promotes Charge Generation and Reduces Nonradiative Voltage Loss in Nonfullerene Organic Solar Cells. J. Phys. Chem. Lett. 2019, 10, 2911–2918. [Google Scholar] [CrossRef]
  69. Coropceanu, V.; Chen, X.-K.; Wang, T.; Zheng, Z.; Brédas, J.-L. Charge-transfer electronic states in organic solar cells. Nat. Rev. Mater. 2019, 4, 689–707. [Google Scholar] [CrossRef]
  70. Dong, Y.; Nikolis, V.C.; Talnack, F.; Chin, Y.-C.; Benduhn, J.; Londi, G.; Kublitski, J.; Zheng, X.; Mannsfeld, S.C.B.; Spoltore, D.; et al. Orientation dependent molecular electrostatics drives efficient charge generation in homojunction organic solar cells. Nat. Commun. 2020, 11, 4617. [Google Scholar] [CrossRef]
  71. Plasser, F. TheoDORE: A toolbox for a detailed and automated analysis of electronic excited state computations. J. Chem. Phys. 2020, 152, 084108. [Google Scholar] [CrossRef] [Green Version]
  72. Hirata, S.; Head-Gordon, M. Time-dependent density functional theory within the Tamm–Dancoff approximation. Chem. Phys. Lett. 1999, 314, 291–299. [Google Scholar] [CrossRef]
  73. Platt, J.R. Classification of Spectra of Cata-Condensed Hydrocarbons. J. Chem. Phys. 1949, 17, 484–495. [Google Scholar] [CrossRef]
  74. Negri, F.; Zgierski, M.Z. Vibronic structure of the emission spectra from single vibronic levels of the S1 manifold in naphthalene: Theoretical simulation. J. Chem. Phys. 1996, 104, 3486. [Google Scholar] [CrossRef]
  75. Knippenberg, S.; Starcke, J.H.; Wormit, M.; Dreuw, A. The low-lying excited states of neutral polyacenes and their radical cations: A quantum chemical study employing the algebraic diagrammatic construction scheme of second order. Mol. Phys. 2010, 108, 2801–2813. [Google Scholar] [CrossRef]
  76. Liang, J.; Feng, X.; Hait, D.; Head-Gordon, M. Revisiting the Performance of Time-Dependent Density Functional Theory for Electronic Excitations: Assessment of 43 Popular and Recently Developed Functionals from Rungs One to Four. J. Chem. Theory Comput. 2022, 18, 3460–3473. [Google Scholar] [CrossRef]
  77. Jacquemin, D.; Mennucci, B.; Adamo, C. Excited-state calculations with TD-DFT: From benchmarks to simulations in complex environments. Phys. Chem. Chem. Phys. 2011, 13, 16987. [Google Scholar] [CrossRef]
  78. Yamagata, H.; Pochas, C.M.; Spano, F.C. Designing J- and H-Aggregates through Wave Function Overlap Engineering: Applications to Poly(3-hexylthiophene). J. Phys. Chem. B 2012, 116, 14494–14503. [Google Scholar] [CrossRef]
  79. Plasser, F.; Lischka, H. Analysis of excitonic and charge transfer interactions from quantum chemical calculations. J. Chem. Theory Comput. 2012, 8, 2777–2789. [Google Scholar] [CrossRef]
  80. Casanova, D.; Krylov, A.I. Quantifying local exciton, charge resonance, and multiexciton character in correlated wave functions of multichromophoric systems. J. Chem. Phys. 2016, 144, 014102. [Google Scholar] [CrossRef]
  81. Bao, P.; Hettich, C.P.; Shi, Q.; Gao, J. Block-Localized Excitation for Excimer Complex and Diabatic Coupling. J. Chem. Theory Comput. 2021, 17, 240–254. [Google Scholar] [CrossRef]
  82. Shirai, S.; Kurashige, Y.; Yanai, T. Computational Evidence of Inversion of 1La and 1Lb-Derived Excited States in Naphthalene Excimer Formation from ab Initio Multireference Theory with Large Active Space: DMRG-CASPT2 Study. J. Chem. Theory Comput. 2016, 12, 2366–2372. [Google Scholar] [CrossRef]
  83. Chai, J.-D.; Head-Gordon, M. Long-range corrected hybrid density functionals with damped atom–atom dispersion corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615. [Google Scholar] [CrossRef] [Green Version]
  84. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  85. Norton, J.E.; Brédas, J.-L. Theoretical characterization of titanyl phthalocyanine as a p-type organic semiconductor: Short intermolecular π-π interactions yield large electronic couplings and hole transport bandwidths. J. Chem. Phys. 2008, 128, 034701. [Google Scholar] [CrossRef] [Green Version]
  86. Löwdin, P. On the Non-Orthogonality Problem Connected with the Use of Atomic Wave Functions in the Theory of Molecules and Crystals. J. Chem. Phys. 1950, 18, 365–375. [Google Scholar] [CrossRef]
  87. Nottoli, M.; Jurinovich, S.; Cupellini, L.; Gardiner, A.T.; Cogdell, R.; Mennucci, B. The role of charge-transfer states in the spectral tuning of antenna complexes of purple bacteria. Photosynth. Res. 2018, 137, 215–226. [Google Scholar] [CrossRef]
Figure 1. The oligoacene dimers considered in this work. Singlet and triplet exciton states have been determined at the eclipsed configuration and along the interchromophore longitudinal (x axis) translation coordinate.
Figure 1. The oligoacene dimers considered in this work. Singlet and triplet exciton states have been determined at the eclipsed configuration and along the interchromophore longitudinal (x axis) translation coordinate.
Molecules 28 00119 g001
Figure 2. Frontier molecular orbital levels and shapes of (left) naphthalene and (right) anthracene, from ωB97X-D/6-31G* calculations (HOMO and LUMO abbreviated as H and L).
Figure 2. Frontier molecular orbital levels and shapes of (left) naphthalene and (right) anthracene, from ωB97X-D/6-31G* calculations (HOMO and LUMO abbreviated as H and L).
Molecules 28 00119 g002
Figure 3. Definition of SA diabatic states ( F E , C R , each symmetry with a specific color code used throughout this work: yellow for A g , orange for A u , blue for B g and green for B u ) defined as linear combinations of diabatic states (black, L E , C T ) and employed to analyze the nature and the effect of interactions on the adiabatic exciton states of oligoacene dimers.
Figure 3. Definition of SA diabatic states ( F E , C R , each symmetry with a specific color code used throughout this work: yellow for A g , orange for A u , blue for B g and green for B u ) defined as linear combinations of diabatic states (black, L E , C T ) and employed to analyze the nature and the effect of interactions on the adiabatic exciton states of oligoacene dimers.
Molecules 28 00119 g003
Figure 4. Lowest energy singlet exciton states of (left) naphthalene and (right) anthracene: (a,b) adiabatic excitation energy profiles of the low-lying exciton states depicted with different color codes for different symmetries: yellow for A g , orange for A u , blue for B g and green for B u . The lowest energy exciton states not included in the diabatization procedure are also shown (dark blue, B g symmetry, dark red, A u symmetry), (c,d) C T character of the lowest four exciton states.
Figure 4. Lowest energy singlet exciton states of (left) naphthalene and (right) anthracene: (a,b) adiabatic excitation energy profiles of the low-lying exciton states depicted with different color codes for different symmetries: yellow for A g , orange for A u , blue for B g and green for B u . The lowest energy exciton states not included in the diabatization procedure are also shown (dark blue, B g symmetry, dark red, A u symmetry), (c,d) C T character of the lowest four exciton states.
Molecules 28 00119 g004
Figure 5. Wavefunction composition of the adiabatic lowest exciton states of (left) naphthalene and (right) anthracene in terms of the SA diabatic states defined in Figure 3. (a,b) 1 A g state and (c,d) 1 B g state. Red and brown lines represent the contributions to the total wavefunction of C T states, while green and blue lines represent the weight of L E states.
Figure 5. Wavefunction composition of the adiabatic lowest exciton states of (left) naphthalene and (right) anthracene in terms of the SA diabatic states defined in Figure 3. (a,b) 1 A g state and (c,d) 1 B g state. Red and brown lines represent the contributions to the total wavefunction of C T states, while green and blue lines represent the weight of L E states.
Molecules 28 00119 g005
Figure 6. Analysis of the excitation energy profiles of singlet exciton states ( A g symmetry) for (left) naphthalene and (right) anthracene dimer (TDA-ωB97X-D/6-31G*) in terms of SA diabatic states (green for F E states, red for C R states) and their interactions. (a,b) Computed adiabatic and SA diabatic excitation energy profiles. (c,d) Magnitude and modulation along the longitudinal translation coordinate of the (grey) D e ( 2 ) D h ( 2 ) , D e ( 4 ) D h ( 4 ) interactions, coupling F E and C R states, and of the (dark-turquoise) H ( 2 , 4 ) V e ( 2 , 4 ) interactions mixing F E ( 2 ) and F E ( 4 ) states.
Figure 6. Analysis of the excitation energy profiles of singlet exciton states ( A g symmetry) for (left) naphthalene and (right) anthracene dimer (TDA-ωB97X-D/6-31G*) in terms of SA diabatic states (green for F E states, red for C R states) and their interactions. (a,b) Computed adiabatic and SA diabatic excitation energy profiles. (c,d) Magnitude and modulation along the longitudinal translation coordinate of the (grey) D e ( 2 ) D h ( 2 ) , D e ( 4 ) D h ( 4 ) interactions, coupling F E and C R states, and of the (dark-turquoise) H ( 2 , 4 ) V e ( 2 , 4 ) interactions mixing F E ( 2 ) and F E ( 4 ) states.
Molecules 28 00119 g006
Figure 7. Analysis of the excitation energy profiles of singlet exciton states ( B g symmetry) for (left) naphthalene and (right) anthracene dimer (TDA-ωB97X-D/6-31G*) in terms of SA diabatic states (green for F E states, red for C R states) and their interactions. (a,b) Computed adiabatic and SA diabatic excitation energy profiles. (c,d) Magnitude and modulation along the longitudinal translation coordinate of the (grey) D e ( 1 ) D h ( 1 ) , D e ( 3 ) D h ( 3 ) interactions, coupling F E and C R states, and of the (dark-turquoise) H ( 1 , 3 ) V e ( 1 , 3 ) interactions mixing F E ( 1 ) and F E ( 3 ) states.
Figure 7. Analysis of the excitation energy profiles of singlet exciton states ( B g symmetry) for (left) naphthalene and (right) anthracene dimer (TDA-ωB97X-D/6-31G*) in terms of SA diabatic states (green for F E states, red for C R states) and their interactions. (a,b) Computed adiabatic and SA diabatic excitation energy profiles. (c,d) Magnitude and modulation along the longitudinal translation coordinate of the (grey) D e ( 1 ) D h ( 1 ) , D e ( 3 ) D h ( 3 ) interactions, coupling F E and C R states, and of the (dark-turquoise) H ( 1 , 3 ) V e ( 1 , 3 ) interactions mixing F E ( 1 ) and F E ( 3 ) states.
Molecules 28 00119 g007
Figure 8. Lowest energy triplet exciton states of (left) naphthalene and (right) anthracene: (a,b) adiabatic excitation energy profiles of the low-lying exciton states depicted with different color codes for different symmetries: yellow for A g , orange for A u , blue for B g and green for B u . The lowest energy exciton states not included in the diabatization procedure are also shown (dark blue, B g symmetry, dark red, A u symmetry); (c,d) C T character of the lowest four exciton states.
Figure 8. Lowest energy triplet exciton states of (left) naphthalene and (right) anthracene: (a,b) adiabatic excitation energy profiles of the low-lying exciton states depicted with different color codes for different symmetries: yellow for A g , orange for A u , blue for B g and green for B u . The lowest energy exciton states not included in the diabatization procedure are also shown (dark blue, B g symmetry, dark red, A u symmetry); (c,d) C T character of the lowest four exciton states.
Molecules 28 00119 g008
Figure 9. Analysis of the excitation energy profiles of triplet exciton states ( B g symmetry) for (left) naphthalene and (right) anthracene dimer (TDA-ωB97X-D/6-31G*) in terms of SA diabatic states (green for F E states, red for C R states) and their interactions. (a,b) Computed adiabatic and SA diabatic excitation energy profiles. (c,d) Magnitude and modulation along the longitudinal translation coordinate of the (grey) D e ( 1 ) D h ( 1 ) , D e ( 3 ) D h ( 3 ) interactions, coupling F E and C R states, and of the (dark-turquoise) H ( 1 , 3 ) V e ( 1 , 3 ) interactions mixing F E ( 1 ) and F E ( 3 ) states.
Figure 9. Analysis of the excitation energy profiles of triplet exciton states ( B g symmetry) for (left) naphthalene and (right) anthracene dimer (TDA-ωB97X-D/6-31G*) in terms of SA diabatic states (green for F E states, red for C R states) and their interactions. (a,b) Computed adiabatic and SA diabatic excitation energy profiles. (c,d) Magnitude and modulation along the longitudinal translation coordinate of the (grey) D e ( 1 ) D h ( 1 ) , D e ( 3 ) D h ( 3 ) interactions, coupling F E and C R states, and of the (dark-turquoise) H ( 1 , 3 ) V e ( 1 , 3 ) interactions mixing F E ( 1 ) and F E ( 3 ) states.
Molecules 28 00119 g009
Figure 10. Analysis of the excitation energy profiles of triplet exciton states ( A u symmetry) for (left) naphthalene and (right) anthracene dimer (TDA-ωB97X-D/6-31G*) in terms of SA diabatic states (green for F E states, red for C R states) and their interactions. (a,b) Computed adiabatic and SA diabatic excitation energy profiles. (c,d) Magnitude and modulation along the longitudinal translation coordinate of the (grey) D e ( 1 ) + D h ( 1 ) , D e ( 3 ) + D h ( 3 ) interactions, coupling F E and C R states, and of the (dark-turquoise) H ( 1 , 3 ) + V e ( 1 , 3 ) interactions mixing F E ( 1 ) and F E ( 3 ) states.
Figure 10. Analysis of the excitation energy profiles of triplet exciton states ( A u symmetry) for (left) naphthalene and (right) anthracene dimer (TDA-ωB97X-D/6-31G*) in terms of SA diabatic states (green for F E states, red for C R states) and their interactions. (a,b) Computed adiabatic and SA diabatic excitation energy profiles. (c,d) Magnitude and modulation along the longitudinal translation coordinate of the (grey) D e ( 1 ) + D h ( 1 ) , D e ( 3 ) + D h ( 3 ) interactions, coupling F E and C R states, and of the (dark-turquoise) H ( 1 , 3 ) + V e ( 1 , 3 ) interactions mixing F E ( 1 ) and F E ( 3 ) states.
Molecules 28 00119 g010
Figure 11. Lowest energy singlet and triplet exciton-state analysis via electron-hole correlation plots [71] for naphthalene and anthracene dimers in their eclipsed configuration. The grey scale used is shown on the right panel. From TDA-ωB97X-D/6-31G* calculations. The magnitude of the C T contribution is shown on the bottom part of each panel.
Figure 11. Lowest energy singlet and triplet exciton-state analysis via electron-hole correlation plots [71] for naphthalene and anthracene dimers in their eclipsed configuration. The grey scale used is shown on the right panel. From TDA-ωB97X-D/6-31G* calculations. The magnitude of the C T contribution is shown on the bottom part of each panel.
Molecules 28 00119 g011
Figure 12. Comparison between vertical absorption spectra predicted for (a) naphthalene and (b) anthracene dimer, from TDA-ωB97X-D/6-31G* calculations. Modulation along the longitudinal shift from 0 Å to 8 Å. Different colors are used to plot spectra computed for different displacements. Black bar: absorption peak of the B b state of the monomer, calculated at TDA-ωB97X-D/6-31G* level.
Figure 12. Comparison between vertical absorption spectra predicted for (a) naphthalene and (b) anthracene dimer, from TDA-ωB97X-D/6-31G* calculations. Modulation along the longitudinal shift from 0 Å to 8 Å. Different colors are used to plot spectra computed for different displacements. Black bar: absorption peak of the B b state of the monomer, calculated at TDA-ωB97X-D/6-31G* level.
Molecules 28 00119 g012
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dai, Y.; Calzolari, A.; Zubiria-Ulacia, M.; Casanova, D.; Negri, F. Intermolecular Interactions and Charge Resonance Contributions to Triplet and Singlet Exciton States of Oligoacene Aggregates. Molecules 2023, 28, 119. https://doi.org/10.3390/molecules28010119

AMA Style

Dai Y, Calzolari A, Zubiria-Ulacia M, Casanova D, Negri F. Intermolecular Interactions and Charge Resonance Contributions to Triplet and Singlet Exciton States of Oligoacene Aggregates. Molecules. 2023; 28(1):119. https://doi.org/10.3390/molecules28010119

Chicago/Turabian Style

Dai, Yasi, Alessandro Calzolari, Maria Zubiria-Ulacia, David Casanova, and Fabrizia Negri. 2023. "Intermolecular Interactions and Charge Resonance Contributions to Triplet and Singlet Exciton States of Oligoacene Aggregates" Molecules 28, no. 1: 119. https://doi.org/10.3390/molecules28010119

APA Style

Dai, Y., Calzolari, A., Zubiria-Ulacia, M., Casanova, D., & Negri, F. (2023). Intermolecular Interactions and Charge Resonance Contributions to Triplet and Singlet Exciton States of Oligoacene Aggregates. Molecules, 28(1), 119. https://doi.org/10.3390/molecules28010119

Article Metrics

Back to TopTop