In the quest for a stable triplet state in small polyaromatic hydrocarbons: an in silico tool for rational design and prediction.

Combining the roles of spin frustration and geometry of odd and even numbered rings in polyaromatic hydrocarbons (PAHs), we design small molecules that show exceedingly small singlet-triplet gaps and stable triplet ground states. Furthermore, a computationally efficient protocol with a model spin Hamiltonian is shown to be capable of qualitative agreement with respect to high level multireference calculations and therefore, can be used for fast molecular discovery and screening. This journal is © The Royal Society of Chemistry 2019.

Introduction

Oligoacenes,1–7 have formed the basis of intense research for a few decades due to their tunable optical gaps. Oligoacenes, up to undecacene,8–10 have been synthesized. Due to their fascinating properties, the theoretical aspects related to their π-bonding topology have been studied extensively,2,4,5,11–13 and their possible usage in singlet fission phenomena14–18 has been investigated.

From a theoretical perspective, the electronic structure of oligoacenes has been the subject of intense debate. Stable triplet states in longer polyacenes have been predicted.2 Multi-reference calculations,11,19,20 however, did not observe triplet ground states, which have been corroborated by experiment.8 However, the exact nature of decay of the singlet triplet (ST) gap was extremely relevant to the synthesis of materials with functionalities important for organic electronics.21–27 Theoretical studies until recently showed that the optical gap decreases monotonically for long polyacenes.28,29 However, a recent study reports the non-monotonicity of the decay in the optical gap with the size of polyacenes.30 Therefore, our understanding of the electronic structure of polyacenes is incomplete to this day.

Azulene, consisting of fused five and seven membered rings, has been much less discussed.31 It is known that azulene has a significantly lower ST gap as compared to its acene isomer, naphthalene. Since the five and seven membered rings do not satisfy the (4n + 2) rule, the singlet state of azulene is significantly less stable, thereby reducing the ST gap. In fact, the dipole moment of the singlet state of azulene points towards charge (or electron) transfer from the seven to five membered ring, to attain a stable electronic state. Repeated fused azulene moieties show low ST gaps and a crossover of the triplet state as the ground state (as opposed to the singlet state).32 This observation is fundamentally different from that of the acenes and was also predicted from the UB3LYP level of theory.33 A previous study by Malrieu and co-workers investigated the development of spin bistability in donor–acceptor complexes and they have proposed a few succinct rules to achieve the goal.34 Much work has also been done in the understanding of the inherent aromaticity of such systems.35 Valence bond theory has been used to understand the root of bistability and radicaloid properties of hydrocarbons.36

In our work, we will show that the UB3LYP method gives rise to certain artifacts and multireference calculations can change the quantitative findings quite significantly. We have designed a class of fused acene–azulenes (Fig. 1) that can reduce the ST gap much more effectively than only repeated azulene units. What is even more significant is that this class of fused acene–azulenes is as stable as polyazulenes, and yet shows a much lower ST gap. We also show that the determining factor behind this large ST gap reduction is multiple spin frustrated pairs or bonds. Thus, the number of spin frustrated bonds can form a simple design principle to effectively tune the ST gap. We have developed a simple protocol for parameterization of a model Hamiltonian that is capable of fast molecular screening and discovery.

Fig. 1

(a) Polyacenes, (b) fused acene–azulenes and (c) polyazulenes.

Computational details

The geometries of polyacenes are taken from earlier work by Hachmann et al.11 and Bendikov et al.2 Geometries of the fused acene–azulenes and polyazulenes are optimized at the (2e,2o)CASCI/RB3LYP/6-31G(d), (10e,10o) CASSCF/6-31G(d) and UB3LYP/6-31G(d) levels of theory. All the optimized geometries for both singlet and triplet states are given in the ESI.† It was noticed that UB3LYP optimizations gave rise to qualitatively different geometries (different bond length alternation (BLA) patterns) as compared to full valence CASSCF optimizations for small PAHs (2 and 3 fused rings) (shown in Fig. 2(a)–(d)). On the other hand, the scheme of obtaining RB3LYP orbitals followed by a minimal CASCI (2e,2o) optimization (referred to as (2e,2o) CASCI/RB3LYP/6-31G(d)) gives the correct BLA pattern (shown in Fig. 2(e) and (f)).

Fig. 2

Comparison of geometries (BLA) of small PAHs at different levels of theory.

Therefore, this computationally affordable scheme of optimization has been used for larger PAHs. In the case of polyacenes, Yang and co-workers37 have also noticed that different levels of theory can give rise to different patterns of BLA and subsequent symmetries. With the minimal active space DFT-CASCI optimized geometries, DMRG (full valence 2p active space),38–42 CASSCF (14e,14o)43 and CASPT2 (14e,14o)44 calculations have been performed using the cc-pVDZ basis set.45 CASSCF and CASPT2 calculations have been performed using the quantum chemistry software package Molpro 2018.2.46 We have used Block-1.5.3 47 for DMRG calculations and Q-Chem 5.1 48 for DFT optimizations.

Tuning the J parameter

For the calculations of the model Hamiltonian, we have considered each C atom as a spin-1/2 site which is antiferromagnetically coupled (favorable interaction between opposite spins) to its nearest neighboring sites with equal coupling parameter J. This essentially corresponds to bonding interactions between the nearest neighbor bonded C atoms. J parameters in the Heisenberg spin Hamiltonian of the polyacenes have been tuned to reproduce the adiabatic DMRG/cc-pVDZ ST gap of the molecular Hamiltonian. These tuned J values are used for isomeric polyazulenes and fused acene–azulenes.

Non-empirical forms of the Heisenberg Hamiltonian have erstwhile49,50 been used to investigate the effect of ground state geometry on the excited state manifold of polyenes. In this work, we have included a simplified and empirical form of the Hamiltonian, in which the absolute energies of the singlet and triplet states have not been used. On the other hand the adiabatic ST gap (obtained from ab initio calculations) has been used to fit the J parameter for each acene. The same J parameter is used for all bonds for all isomers of acenes, i.e., polyazulene and the fused acene–azulene systems. Therefore, all the bonds in the model Hamiltonian are treated as equivalent and the only direct interactions that are considered are those across the bonds in the moiety.

Results and discussion

Adiabatic ST gap

Three types of structures (shown in Fig. 1) are considered and compared – (a) polyacenes, (b) fused acene–azulenes and (c) polyazulenes. The adiabatic ST gaps, calculated at the DMRG/cc-pVDZ level of theory, for the fused rings of size 2–7 are shown in Fig. 3(a). The complete 2p valence space is chosen as the active space. As expected, the azulene ST gap is significantly lower than that of its isomer, naphthalene. However, we notice that the rate of decay of the adiabatic ST gap for polyazulene is much lower than that of the polyacenes. Therefore, contrary to what has been observed at the model Hamiltonian32 and UDFT level of theory,33 accurate multireference calculations show that the polyazulenes do not form stable triplet states within 6 fused rings.

Fig. 3

Variation of the ST gap for the acene, azulene and fused acene–azulene systems with (a) the DMRG/cc-pVDZ level of theory and (b) spin-Hamiltonian model calculations. (c) Correlation between the ST gaps with model and molecular Hamiltonians. The black circles refer to the acenes (used for tuning) and red circles are for other PAHs. The fitness measure or correlation coefficient is 0.96 for all the PAHs.

Surprisingly a combined acene–azulene fused system (Fig. 1(b)) shows the lowest adiabatic ST gaps. Within 5 fused rings the ST gap of the fused acene–azulene system is close to the computational accuracy of the method (∼0.05 eV). The ST gaps with various levels of multireference theory (CASSCF and CASPT2) show the same qualitative trend (given in the ESI†). It is important to note that when the acenes and azulenes are fused in different orders,51i.e., structures with fused 7–5–6 or 7–5–6–6 rings, they show much higher (1.58 and 1.40 eV respectively) adiabatic ST gaps than 7–6–5 or 7–6–6–5 motifs (0.97 and 0.46 eV respectively) calculated at the DMRG/cc-pVDZ level of theory. Furthermore, a comparison of the stabilization energies per C atom of the isomeric systems of PAHs shows that the fused acene–azulenes systems (8–9 kcal mol–1) are comparable in stability with the polyazulenes (8–9 kcal mol–1). But both of these systems are somewhat less stable than the acene–azulene systems that have already been synthesized (10–11 kcal mol–1), such as 7–5–6 and 7–5–6–6 moieties.51 These stabilization energies were calculated with the help of isodesmic reactions of the form,and the extra stabilization energy with respect to an equal number of C–C single and double bonds and C–H bonds is computed. The details of the stabilization energies of all the systems are included in the ESI.†

The single reference UDFT calculations when compared to the multireference calculations show that for most of these systems the UDFT level of theory gives rise to quantitatively wrong results (preferential stabilization of triplet states). On the other hand, the multireference calculations are computationally expensive and cannot be routinely used for molecular design and screening. Therefore, we try to develop a model Hamiltonian based protocol to predict the ST gaps in a computationally affordable way. It is important to note that such model Hamiltonians have been suggested in the context of polyacenes in previous studies by a few groups.50,52,53 The ST gaps are calculated with a model Heisenberg Hamiltonian, given by,where S denotes the spin at the ith site (C atom) and J is the coupling constant between neighboring spins (on C atoms), which signifies the effective bonding interaction. J is parameterized with respect to ST gaps of acenes.54 Since here the effective Hamiltonian thus formed is dependent only on the nearest neighbor interactions between bonded p orbitals, the computational cost is significantly reduced. The ST gaps calculated with this model Hamiltonian (Fig. 3(b)) show the same qualitative trends as the ab initio calculations (Fig. 3(a)). The correlation between these calculations (Fig. 3(c)) shows that the ST gaps estimated with the model Hamiltonian are in excellent agreement with those of the molecular Hamiltonian and can therefore be used as an effective screening and prediction tool.

Spin frustration

The model Hamiltonian wavefunction gives us a simple paradigm to understand the physics behind ST gaps in these molecules. It is important to note that the spins in adjacent sites (bonded C atoms) are antiferromagnetically coupled, so as to model the bonding interactions between them. Thus, opposite spins on adjacent sites refer to bonding and therefore, stabilization, while same spins on adjacent sites destabilize the configurations. The even-membered rings can contain configurations with all anti-parallel adjacent spins (i.e. stable configurations), while this is not possible in the case of odd-membered rings. This leads to spin frustration,55 in the odd-membered rings. Examples of such configurations are shown in Fig. 4(a) for polyazulene and fused acene–azulene. There is effectively one frustrated pair for each azulene unit (2 fused rings). In case of fused acene–azulene, this situation may be exacerbated. It should be noted that in the fused acene–azulene system, configurations with at least 2 spin frustrations are possible. These observations can also be explained within generalizations of Ovchinnikov's rule.

Fig. 4

(a) Example configurations of spins on the polyazulene and fused acene–azulene systems. (b) Probabilities of same spins across bonds.

The probability of finding same spins in adjacent sites can, therefore, form a measure of destabilization and is shown in Fig. 4(b). Since the acenes are made of even-membered rings only, there is a low probability of encountering the same spins on adjacent sites. However, in the case of both polyazulene and fused acene–azulene there are regions of high probabilities, especially across certain common bonds (bonds shared by two fused rings). To measure the instability in each system, the probabilities of the same spins across any bonds are computed. It is 0.22 for tetracene, 0.24 for 4-azulene and 0.24 for the 4-fused acene–azulene system. We notice that in the case of polyazulenes, each azulene (7–5 membered fused ring) forms a unit and the common bond between the 7 and 5 membered rings show a higher probability of the same spin configuration. In the case of fused acene–azulene, this trend is spread across all common bonds and sometimes even across the rings. These observations point towards the presence of spin frustration in polyazulenes and fused acene–azulenes, as opposed to acenes.

From the comparison of the probabilities it is expected that the polyazulenes and fused acene–azulenes would be of comparable stability in their singlet state and both these compounds would be less stable than acenes. This is indeed noticed from the ab initio calculations. Furthermore, one would expect that the regions of high same-spin probability would try to reduce the instability by introducing lattice degrees of freedom, i.e., change in molecular geometry. In line with our expectations, the geometries of the fused acene–azulene and polyazulene systems (shown in Fig. 5) show elongated bond lengths across the common bonds with high same-spin probability (Fig. 4(b)). This is akin to Peierls distortion in polyenes and polyacenes.

Fig. 5

Optimized geometries of singlet states of (a) fused acene–azulene, and (b) polyazulene.

In the polyacenes we notice the expected BLA, where the fused bonds also participate in the BLA pattern.56 In the case of polyazulene (Fig. 5(a)), we notice that each azulene moiety forms a unit and the common bond in this unit is elongated. This common bond is almost a single bond in nature (1.48 Å), i.e., only minimally participates in conjugation. Thus, each azulene forms a pseudo ten membered conjugated ring pattern. However, the common bond between the two azulene moieties is shorter than a single bond (1.39 Å) and takes part in conjugation. The situation is quite different when we notice the bond lengths of the fused acene–azulene system (Fig. 5(b)). Here, all the common bonds are elongated to almost single bond lengths (1.47–1.48 Å) and all of these bonds are minimally present in conjugation. These observations are in accordance with that of the spin probabilities calculated in the model system.

Stability of the triplet state

The discussion till now considered only the relative stability of the singlet states and showed that both polyazulenes and fused acene–azulenes are expected to have low ST gaps. However, in the case of polyazulenes, the triplet state is destabilized with respect to that of fused acene–azulene (Fig. 6). This can be explained from the spin probabilities being in the same line of argument as those for singlets. The probabilities of finding the same spins in adjacent sites for the triplet states are 0.255 in polyazulene and 0.24 in fused acene–azulene, in isomeric C18H12. Thus, the ST gap of polyazulenes is not as low as expected due to its unstable singlet state. Polyazulenes do not show a ST crossover within 6 fused rings.

Fig. 6

Molecular and model Hamiltonian energies for the singlet and triplet states are compared for 4-fused ring systems, i.e., tetracene, polyazulene and fused acene–azulene.

As the size of the molecule grows, it is the triplet stabilization that takes over. For example, in isomers of C26H16, i.e., containing 6 rings, the singlet state of fused acene–azulene shows lower spin frustration than polyazulene. Since, fused acene–azulenes can form configurations with at least 2 spin frustrations as opposed to polyazulenes with 3 spin frustrations, this is to be expected. This is depicted in the same spin probabilities of the singlet states as 0.246 for polyazulene versus 0.240 for fused acene–azulene, in C26H16. Therefore, one would expect a lower ST gap for polyazulene from the consideration of the singlet state only. However, at these system sizes, it is the relative triplet state stabilities that result in the ST gaps of fused acene–azulenes being much lower than those of the polyazulenes. The same spin probabilities for the triplet states are 0.251 for polyazulene versus 0.240 for fused acene–azulene, in C26H16. Thus, we can notice that the probabilities of finding the same adjacent spins remain almost constant in the case of fused acene–azulene systems across system sizes, while in polyazulenes they increase as expected due to more number of spin frustrations.

Symmetry of the molecule also plays a minor role in the ST gaps. The difference in ST gaps of straight and kinked acenes has been reported before.57 The same is noticed here. The bent PAHs, i.e., polyazulenes, have higher ST gaps than the straight isomers, fused acene–azulene. It was further noticed that the fused acene–azulenes systems which were arranged in a bent topology also show relatively higher ST gaps. Therefore, if one were to design a low ST gap molecule, there are three considerations – (i) maximize spin frustrations in singlet states, (ii) minimize spin frustrations in triplet states and (iii) more symmetric molecules. With these considerations in mind, a large number of differently fused acene–azulene systems have been designed and screened. We notice that the lowest ST gap moiety consists of a fused ring motif of 7–6···6–5 membered rings.

Bond length alternation

Bond length alternation (BLA) patterns of the PAHs have been found to be quite different using different optimization schemes (details given in the ESI†). Furthermore, it is known in the context of polyenes that inclusion of only static correlation can give rise to more pronounced BLAs. In this study, we have used CASSCF and CASCI based approaches to ascertain the geometry and therefore, the BLA might be over-estimated. Therefore, it is important to understand the effect of lower BLAs on the ST gaps. Fig. 7(a) shows the effect of reduction in the BLAs on the vertical ST gaps of the hydrocarbons. The BLA is reduced to about half the amount obtained from DFT-CASCI optimization to understand the effect of reduction of BLA on the ST gaps. Furthermore, Fig. 7(b) shows the effect of change of BLA to different values on the tetracene (6–6–6–6), 7–5–7–5 and 7–6–6–5 systems, as a representative example. It can be seen from both the figures that as expected, the reduction in BLA stabilizes the triplet state and therefore, reduces the vertical ST gap. The comparison of the triplet and singlet state geometries further shows that low BLA is favored by triplet states and high BLA is favored by singlet states. Furthermore, it is interesting to note that the effect of BLA reduction on the ST gap of fused acene–azulenes is markedly more extensive than in the case of polyacenes and polyazulenes.

Fig. 7

(a) Variation of the vertical ST gap with change in BLA in PAHs, (b) reduction of the vertical ST gap with the decrease of the BLA parameter for 7–5–7–5, 6–6–6–6, and 7–6–6–5 systems.

Conclusions

In conclusion, we have developed a novel model Hamiltonian based approach for the design of PAHs with tunable ST gaps. This analogy between the molecular system and model spin lattice can be expressed as a generalization of Ovchinnikov's rule, which states the existence of ferromagnetic phases in infinite hydrocarbon systems of appropriate alternating or non-alternating patterns. While this rule discusses the existence of high spin states in infinite or large systems and can be easily correlated with odd membered ring systems such as poly-azulenes, we find it fascinating that fused acene–azulenes which can indeed be thought of as an acene system attached to two odd membered rings in its extremity can exacerbate the effect and give rise to high spin stable states with much smaller sizes.

We have noticed that the ST gaps predicted by the model Heisenberg Hamiltonian approach are well correlated with the ab initio calculations. Therefore, this approach has been used to predict the smallest PAH with a negative ST gap. The stability of these species was tested with the help of isodesmic reactions and they were shown to be reasonably stable with respect to their acene analogues. We have further been able to propose a few succinct design principles which can be used to engineer low ST gap species. In this study, we have tested only the systems that are isomeric with polyacenes; however, we believe that this method can be used for other PAH systems as well.

Conflicts of interest

There are no conflicts to declare.

Click here for additional data file.