Do any types of double-hybrid models render the correct order of excited state energies in inverted singlet–triplet emitters?

Benchmark computations of nearly degenerate singlet and triplet states of N-heterocyclic chromophores. II. Density-based methods

In this paper, we demonstrate the performance of several density-based methods in predicting the inversion of S1 and T1 states of a few N-heterocyclic triangulene based fused ring molecules (popularly known as INVEST molecules) with an eye to identify a well performing but cost-effective preliminary screening method. Both conventional linear-response time-dependent density functional theory (LR-TDDFT) and ΔSCF methods (namely maximum overlap method, square-gradient minimization method, and restricted open-shell Kohn-Sham) are considered for excited state computations using exchange-correlation (XC) functionals from different rungs of Jacob's ladder. A well-justified systematism is observed in the performance of the functionals when compared against fully internally contracted multireference configuration interaction singles and doubles and/or equation of motion coupled-cluster singles and doubles (EOM-CCSD), with the most important feature being the capture of spin-polarization in the presence of correlation. A set of functionals with the least mean absolute error is proposed for both the approaches, LR-TDDFT and ΔSCF, which can be more cost-effective alternatives for computations on synthesizable larger derivatives of the templates studied here. We have based our findings on extensive studies of three cyclazine-based molecular templates, with additional studies on a set of six related templates. Previous benchmark studies for subsets of the functionals were conducted against the domain-based local pair natural orbital-similarity transformed EOM-CCSD (STEOM-CCSD), which resulted in an inadequate evaluation due to deficiencies in the benchmark theory. The role of exact-exchange, spin-contamination, and spin-polarization in the context of DFT comes to the forefront in our studies and supports the numerical evaluation of XC functionals for these applications. Suitable connections are drawn to two and three state exciton models, which identify the minimal physics governing the interactions in these molecules.

Machine Learning-Corrected Simplified Time-Dependent DFT for Systems With Inverted Single-t-o-Triplet Gaps of Interest for Photocatalytic Water Splitting

Hydrogen gas (H2) can be produced via entirely solar-driven photocatalytic water splitting (PWS). A promising set of organic materials for facilitating PWS are the so-called inverted singlet-triplet, INVEST, materials. Inversion of the singlet (S1) and triplet (T1) energies reduces the population of triplet states, which are otherwise destructive under photocatalytic conditions. Moreover, when INVEST materials possess dark S1 states, the excited state lifetimes are maximized, facilitating energy transfer to split water. In the context of solar-driven processes, it is also desirable that these INVEST materials absorb near the solar maximum. Many aza-triangulenes possess the desired INVEST property, making it beneficial to describe an approach for systematically and efficiently predicting the INVEST property as well as properties that make for efficient photocatalytic water splitting, while exploring the large chemical space of the aza-triangulenes. Here, we utilize machine learning to generate post hoc corrections to simplified Tamm-Dancoff approximation density functional theory (sTDA-DFT) for singlet and triplet excitation energies that are within 28-50 meV of second-order algebraic diagrammatic construction, ADC(2), as well as the singlet-to-triplet, ΔES1T1, gaps of PWS systems. Our Δ-ML model is able to recall 85% of the systems identified by ADC(2) as candidates for PWS. Further, with a modest database of ADC(2) excitation energies of 4025 aza-triangulenes, we identified 78 molecules suitable for entirely solar-driven PWS.

Benchmark computations of nearly degenerate singlet and triplet states of N-heterocyclic chromophores. I. Wavefunction-based methods

In this paper, we investigate the role of electron correlation in predicting the S1-S0 and T1-S0 excitation energies and, hence, the singlet-triplet gap (ΔEST) in a set of cyclazines, which act as templates for potential candidates for fifth generation organic light emitting diode materials. This issue has recently garnered much interest with the focus being on the inversion of the ΔEST, although experiments have indicated near degenerate levels with both positive and negative being within the experimental error bar [J. Am. Chem. Soc. 102, 6068 (1980), J. Am. Chem. Soc. 108, 17(1986)]. We have carried out a systematic and exhaustive study of various excited state electronic structure methodologies and identified the strengths and shortcomings of the various approaches and approximations in view of this challenging case. We have found that near degeneracy can be achieved either with a proper balance of static and dynamic correlation in multireference theories or with state-specific orbital corrections, including its coupling with correlation. The role of spin contamination is also discussed. Eventually, this paper seeks to produce benchmark numbers for establishing cost-effective theories, which can then be used for screening derivatives of these templates with desirable optical and structural properties. Additionally, we would like to point out that the use of domain-based local pair natural orbital-similarity transformed EOM-coupled cluster singles and doubles as the benchmark for ΔEST [as used in J. Phys. Chem. A 126(8), 1378 (2022), Chem. Phys. Lett. 779, 138827 (2021)] is not a suitable benchmark for these classes of molecules.

Design of Anti-Hund Organic Emitters Based on Heptazine

Hund's rule, which is powerful in governing the first excited states of closed-shell organic materials, can hardly be violated to get inverted singlet-triplet gap (INVEST) molecules with negative singlet-triplet energy gaps (ΔEST), although INVEST materials have shown extraordinary photophysical properties and promising device performance especially in light-emitting diodes. Here, we propose a facile strategy to construct emissive INVEST molecules by introducing different types of substituents to heptazine in various modes, which can effectively tune the ΔEST to be negative with the enlarged oscillator strength (f) for the high fluorescence rate of the heptazine derivatives. Systematic computational studies show that the double substitution of electron-donating units with another nonconjugated substituent in hybrid substitution mode is the most favorable way in achieving slightly negative ΔEST and large f values; the conjugated substituent will compete with heptazine to make the molecule deviate from the INVEST feature. Especially, a series of high-performance heptazine-based INVEST emitters were constructed, exhibiting ΔEST low to -0.362 eV, f up to 0.0436, as well as a wide range emission color from 339 to 716 nm. Also, the designed molecules were predicted to have fluorescence radiative rates up to 106 s-1, along with efficient reverse intersystem crossing rates reaching 108 s-1. Importantly, the figure of merit (FM) was first proposed as a parameter to wholly evaluate the performance of INVEST emitters, and the highest FM of 0.198 was found in the triazine and double nonconjugated amine-substituted heptazine. These results highlight the great potential of the heptazine chromophore in constructing INVEST emitters, revealing fundamental structure-property understandings for the material design of efficient anti-Hund organic molecules with improved emission properties.

Reverse Intersystem Crossing Dynamics in Vibronically Modulated Inverted Singlet-Triplet Gap System: A Wigner Phase Space Study

We inspect the origin of the inverted singlet-triplet gap (INVEST) and slow change in the reverse intersystem crossing (rISC) rate with temperature, as recently observed. A Wigner phase space study reveals that, though INVEST is found at equilibrium geometry, variation in the exchange interaction and the doubles-excitation for other geometries in the harmonic region leads to non-INVEST behavior. This highlights the importance of nuclear degrees of freedom for the INVEST phenomenon, and in this case, geometric puckering of the studied molecule determines INVEST and the associated rISC dynamics.

Spectroscopic Quantification of the Inverted Singlet-Triplet Gap in Pentaazaphenalene

We report the direct and accurate spectroscopic quantification of the inverted singlet-triplet gap in 1,3,4,6,9b-pentaazaphenalene. This measurement is achieved by directly probing the lowest singlet and triplet states via high-resolution cryogenic anion photoelectron spectroscopy. The assignment of the first excited singlet state is confirmed by visible absorption spectroscopy in an argon matrix at 20 K. Our measurements yield an inverted singlet-triplet gap with ΔEST= -0.047(7) eV. The accurate quantification of the singlet-triplet gap presented here allows for direct evaluation of various computational electronic structure methods and highlights the critical importance of the proper description of the double excitation character of these electronic states. Overall, this study validates the idea that despite Hund's multiplicity rule, useful organic chromophores can have inherently inverted singlet-triplet gaps.

Resilience of Hund's rule in the chemical space of small organic molecules

We embark on a quest to identify small molecules in the chemical space that can potentially violate Hund's rule. Utilizing twelve TDDFT approximations and the ADC(2) many-body method, we report the energies of S1 and T1 excited states of 12 880 closed-shell organic molecules within the bigQM7ω dataset with up to 7 CONF atoms. In this comprehensive dataset, none of the molecules, in their minimum energy geometry, exhibit a negative S1-T1 energy gap at the ADC(2) level while several molecules display values <0.1 eV. The spin-component-scaled double-hybrid method, SCS-PBE-QIDH, demonstrates the best agreement with ADC(2). Yet, at this level, a few molecules with a strained sp3-N center turn out as false-positives with the S1 state lower in energy than T1. We investigate a prototypical cage molecule with an energy gap <-0.2 eV, which a closer examination revealed as another false positive. We conclude that in the chemical space of small closed-shell organic molecules, it is possible to identify geometric and electronic structural features giving rise to S1-T1 degeneracy; still, there is no evidence of a negative gap. We share the dataset generated for this study as a module, to facilitate seamless molecular discovery through data mining.

Heptazine, Cyclazine, and Related Compounds: Chemically-Accurate Estimates of the Inverted Singlet-Triplet Gap

Molecules that violate Hund's rule and exhibit an inverted gap between the lowest singlet S1 and triplet T1 excited states have attracted considerable attention due to their potential applications in optoelectronics. Among these molecules, the triangular-shaped heptazine, and its derivatives, have been in the limelight. However, conflicting reports have arisen regarding the relative energies of S1 and T1. Here, we employ highly accurate levels of theory, such as CC3, to not only resolve the debate concerning the sign but also quantify the magnitude of the S1-T1 gap. We also determined the 0-0 energies to evaluate the significance of the vertical approximation. In addition, we compute reference S1-T1 gaps for a series of 10 related molecules. This enables us to benchmark lower-order methods for future applications in larger systems within the same family of compounds. This contribution can serve as a foundation for the design of triangular-shaped molecules with enhanced photophysical properties.

Computational Design of Multiple Resonance-Type BN Molecules for Inverted Singlet and Triplet Excited States

A computational design of linearly extended multiple resonance (MR)-type BN molecules based on DABNA-1 is proposed herein in the quest to find potential candidates that exhibit a negative singlet-triplet gap (ΔEST) and a large oscillator strength value. The impact of a proper account of the electron correlation in the lowest singlet and triplet excited states is systematically investigated by using double-hybrid functionals within the TD-DFT framework, as well as wavefunction-based methods (EOM-CCSD and SCS-CC2), since this contribution plays an essential role in driving the magnitude of the ΔEST in MR-TADF and inverted singlet-triplet gap compounds. Our results point out a gradual reduction of the ΔEST gap with respect to the increasing sum of the number of B and N atoms, reaching negative ΔEST values for some molecules as a function of their size. The double-hybrid functionals reproduce the gap with only slight deviation compared to available experimental data for DABNA-1, ν-DABNA, and mDBCz and nicely agree with high-level quantum mechanical methods (e.g., EOM-CCSD and SCS-CC2). Larger oscillator strengths are found compared to the azaphenalene-type molecules, also exhibiting the inversion of their singlet and triplet excited states. We hope this study can serve as a motivation for further design of the molecules showing negative ΔEST based on boron- and nitrogen-doped polyaromatic hydrocarbons.

Inverted Lowest Singlet and Triplet Excitation Energy Ordering of Graphitic Carbon Nitride Flakes

In organic light-emitting diodes (OLEDs), only 25% of electrically generated excitons are in a singlet state, S1, and the remaining 75% are in a triplet state, T1. In thermally activated delayed fluorescence (TADF) chromophores the transition from the nonradiative T1 state to the radiative S1 state can be thermally activated, which improves the efficiency of OLEDs. Chromophores with inverted energy ordering of S1 and T1 states, S1 < T1, are superior to TADF chromophores, thanks to the absence of an energy barrier for the transition from T1 to S1. We benchmark the performance of time-dependent density functional theory using different exchange-correlation functionals and find that scaled long-range corrected double-hybrid functionals correctly predict the inverted singlet-triplet gaps of N-substituted phenalene derivatives. We then show that the inverted energy ordering of S1 and T1 is an intrinsic property of graphitic carbon nitride flakes. A design strategy of new chromophores with inverted singlet-triplet gaps is proposed. The color of emitted light can be fine-tuned through flake size and amine substitution on flake vertices.

Building on the strengths of a double-hybrid density functional for excitation energies and inverted singlet-triplet energy gaps

It is demonstrated that a double hybrid density functional approximation, ωB88PTPSS, that incorporates equipartition of density functional theory and the non-local correlation, however with a meta-generalized gradient approximation correlation functional, as well as with the range-separated exchange of ωB2PLYP, provides accurate excitation energies for conventional systems, as well as correct prescription of negative singlet-triplet gaps for non-conventional systems with inverted gaps, without any necessity for parametric scaling of the same-spin and opposite-spin non-local correlation energies. Examined over "safe" excitations of the QUESTDB set, ωB88PTPSS performs quite well for open-shell systems, correctly and fairly accurately [relative to equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) reference] predicts negative gaps for 50 systems with inverted singlet-triplet gaps, and is one of the leading performers for intramolecular charge-transfer excitations and achieves near-second-order approximate coupled cluster (CC2) and second-order algebraic diagrammatic construction quality for the Q1 and Q2 subsets. Subsequently, we tested ωB88PTPSS on two sets of real-life examples from recent computational chemistry literature-the low energy bands of chlorophyll a (Chl a) and a set of thermally activated delayed fluorescence (TADF) systems. For Chl a, ωB88PTPSS qualitatively and quantitatively achieves DLPNO-STEOM-CCSD-level performance and provides excellent agreement with experiment. For TADF systems, ωB88PTPSS agrees quite well with spin-component-scaled CC2 (SCS-CC2) excitation energies, as well as experimental values, for the gaps between the S₁ and T₁ excited states.

Organic molecules with inverted singlet-triplet gaps

According to Hund's multiplicity rule, the energy of the lowest excited triplet state (T1) is always lower than that of the lowest excited singlet state (S1) in organic molecules, resulting in a positive singlet-triplet energy gap (ΔE ST). Therefore, the up-converted reverse intersystem crossing (RISC) from T1 to S1 is an endothermic process, which may lead to the quenching of long-lived triplet excitons in electroluminescence, and subsequently the reduction of device efficiency. Interestingly, organic molecules with inverted singlet-triplet (INVEST) gaps in violation of Hund's multiplicity rule have recently come into the limelight. The unique feature has attracted extensive attention in the fields of organic optoelectronics and photocatalysis over the past few years. For an INVEST molecule possessing a higher T1 with respect to S1, namely a negative ΔE ST, the down-converted RISC from T1 to S1 does not require thermal activation, which is possibly conducive to solving the problems of fast efficiency roll-off and short lifetime of organic light-emitting devices. By virtue of this property, INVEST molecules are recently regarded as a new generation of organic light-emitting materials. In this review, we briefly summarized the significant progress of INVEST molecules in both theoretical calculations and experimental studies, and put forward suggestions and expectations for future research.