Time-Dependent Long-Range-Corrected Double-Hybrid Density Functionals with Spin-Component and Spin-Opposite Scaling: A Comprehensive Analysis of Singlet-Singlet and Singlet-Triplet Excitation Energies

How do spin-scaled double hybrids designed for excitation energies perform for noncovalent excited-state interactions? An investigation on aromatic excimer models

Time-dependent double hybrids with spin-component or spin-opposite scaling to their second-order perturbative correlation correction have demonstrated competitive robustness in the computation of electronic excitation energies. Some of the most robust are those recently published by our group (M. Casanova-Páez, L. Goerigk, J. Chem. Theory Comput. 2021, 20, 5165). So far, the implementation of these functionals has not allowed correctly calculating their ground-state total energies. Herein, we define their correct spin-scaled ground-state energy expressions which enables us to test our methods on the noncovalent excited-state interaction energies of four aromatic excimers. A range of 22 double hybrids with and without spin scaling are compared to the reasonably accurate wavefunction reference from our previous work (A. C. Hancock, L. Goerigk, RSC Adv. 2023, 13, 35964). The impact of spin scaling is highly dependent on the underlying functional expression, however, the smallest overall errors belong to spin-scaled functionals with range separation: SCS- and SOS- ω PBEPP86, and SCS-RSX-QIDH. We additionally determine parameters for DFT-D3(BJ)/D4 ground-state dispersion corrections of these functionals, which reduce errors in most cases. We highlight the necessity of dispersion corrections for even the most robust TD-DFT methods but also point out that ground-state based corrections are insufficient to completely capture dispersion effects for excited-state interaction energies.

Cobalt-Mediated Photochemical C-H Arylation of Pyrroles

Precious metal complexes remain ubiquitous in photoredox catalysis (PRC) despite concerted efforts to find more earth-abundant catalysts and replacements based on 3d metals in particular. Most otherwise plausible 3d metal complexes are assumed to be unsuitable due to short-lived excited states, which has led researchers to prioritize the pursuit of longer excited-state lifetimes through careful molecular design. However, we report herein that the C-H arylation of pyrroles and related substrates (which are benchmark reactions for assessing the efficacy of photoredox catalysts) can be achieved using a simple and readily accessible octahedral bis(diiminopyridine) cobalt complex, [1-Co](PF6)2. Notably, [1-Co]2+ efficiently functionalizes both chloro- and bromoarene substrates despite the short excited-state lifetime of the key photoexcited intermediate *[1-Co]2+ (8 ps). We present herein the scope of this C-H arylation protocol and provide mechanistic insights derived from detailed spectroscopic and computational studies. These indicate that, despite its transient existence, reduction of *[1-Co]2+ is facilitated via pre-assembly with the NEt3 reductant, highlighting an alternative strategy for the future development of 3d metal-catalyzed PRC.

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.

Determining Excited-State Absorption Properties of a Quinoid Flavin by Polarization-Resolved Transient Spectroscopy

As important naturally occurring chromophores, photophysical/chemical properties of quinoid flavins have been extensively studied both experimentally and theoretically. However, little is known about the transition dipole moment (TDM) orientation of excited-state absorption transitions of these important compounds. This aspect is of high interest in the fields of photocatalysis and quantum control studies. In this work, we employ polarization-associated spectra (PAS) to study the excited-state absorption transitions and the underlying TDM directions of a standard quinoid flavin compound. As compared to transient absorption anisotropy (TAA), an analysis based on PAS not only avoids diverging signals but also retrieves the relative angle for ESA transitions with respect to known TDM directions. Quantum chemical calculations of excited-state properties lead to good agreement with TA signals measured in magic angle configuration. Only when comparing experiment and theory for TAA spectra and PAS, do we find deviations when and only when the S0 → S1 of flavin is used as a reference. We attribute this to the vibronic coupling of this transition to a dark state. This effect is only observed in the employed polarization-controlled spectroscopy and would have gone unnoticed in conventional TA.

Tris-Azo Triangular Paraphenylenes: Synthesis and Reversible Interconversion into Radial π-Conjugated Macrocycles

We report the synthesis of cycloparaphenylene derivatives featuring tris-azo groups. The smaller derivative, [3]cycloazobenzene, adopts a triangular all-cis form and exhibits thermally and photochemically stable characteristics due to significant ring strain as well as symmetric Kagome-patterned crystal packing. In contrast, the as-synthesized [3]cycloazobenzene with three biphenylene bridges adopts a triangular all-cis form, which undergoes photoinduced isomerization, leading to a photostationary state. Interestingly, the addition of an excess of acid selectively leads to the formation of an all-trans form. DFT calculations reveal that the interconversion from a triangular to a circular shape correlates with an increase in HOMO and a decrease in LUMO, characteristics intrinsic to radial π-conjugated systems.

Benzazole-Based ESIPT Fluorophores: Proton Transfer from the Chalcogen Perspective. A Combined Theoretical and Experimental Study

In this study, we present a comprehensive photophysical investigation of ESIPT-reactive benzazole derivatives in both solution and the solid state. These derivatives incorporate different chalcogen atoms (O, S, and Se) into their structures, and we explore how these variations impact their electronic properties in both ground and excited states. Changes in the UV-Vis absorption and fluorescence emission spectra were analyzed and correlated with the chalcogen atom and solvent polarity. In general, the spectral band of the benzazole derivative containing selenium was redshifted in both the ground and excited states compared to that of its oxygen and sulfur counterparts. Furthermore, we observed that the solvent played a distinctive role in influencing the ESIPT process within these compounds, underscoring once again the significant influence of the chalcogen atom on their photophysical behavior. Theoretical calculations provided a deeper understanding of the molecular dynamics, electronic structures, and photophysical properties of these compounds. These calculations highlighted the effect of chalcogen atoms on the molecular geometry, absorption and emission characteristics, and intramolecular hydrogen bonding, revealing intricate details of the ESIPT mechanism. The integration of experimental and computational data offers a detailed view of the structural and electronic factors governing the photophysical behavior of benzazole derivatives.

Contrasting the excited state properties of different conformers of trans- and cis-2,2'-bipyridine oligomers in the gas phase

In this article, we present conformation-dependent photophysical and excited state properties of trans- and cis- BPY oligomers. Oligomers up to tetramers for three conformers, namely, o-, m-, and p-, are constructed and optimized at the B3LYP-D3/def2-SVPD level. The photophysical and excited state properties are interpreted in terms of UV and CD spectra at the RI-ADC(2)/def2-TZVPD level. The UV spectra of oligomers of the m-conformer show high-intensity and red-shifted UV bands compared to o- and p-oligomers. The CD spectra of p-oligomers show intense CD bands compared to o- and p-oligomers in the case of trans-structures. In contrast, oligomers of each conformer of cis-structures show high-intensity CD bands. The excited states of (BPY)2 and (BPY)4 are also characterized by analysis of one-electron transition density matrix considering three descriptors: ωCT, dexc, and PRNTO. The ωCT values of dimers are in the range of 0.06-0.32, which indicates the excited states are mainly LE states, whereas, for (BPY)4, the ωCT values range from 0.17 to 0.53, indicating the possibility of partial CT in the excited states. These observations are also explained using the NTOs and e-h correlation plots.

Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models

Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for time-dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.

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.

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.

Benchmarking Density Functional Theory Methods for Metalloenzyme Reactions: The Introduction of the MME55 Set

We present a new benchmark set of metalloenzyme model reaction energies and barrier heights that we call MME55. The set contains 10 different enzymes, representing eight transition metals, both open and closed shell systems, and system sizes of up to 116 atoms. We use four DLPNO-CCSD(T)-based approaches to calculate reference values against which we then benchmark the performance of a range of density functional approximations with and without dispersion corrections. Dispersion corrections improve the results across the board, and triple-ζ basis sets provide the best balance of efficiency and accuracy. Jacob's ladder is reproduced for the whole set based on averaged mean absolute (percent) deviations, with the double hybrids SOS0-PBE0-2-D3(BJ) and revDOD-PBEP86-D4 standing out as the most accurate methods for the MME55 set. The range-separated hybrids ωB97M-V and ωB97X-V also perform well here and can be recommended as a reliable compromise between accuracy and efficiency; they have already been shown to be robust across many other types of chemical problems, as well. Despite the popularity of B3LYP in computational enzymology, it is not a strong performer on our benchmark set, and we discourage its use for enzyme energetics.

Reference Vertical Excitation Energies for Transition Metal Compounds

To enrich and enhance the diversity of the quest database of highly accurate excitation energies [Véril, M.; et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], we report vertical transition energies in transition metal compounds. Eleven diatomic molecules with a singlet or doublet ground state containing a fourth-row transition metal (CuCl, CuF, CuH, ScF, ScH, ScO, ScS, TiN, ZnH, ZnO, and ZnS) are considered, and the corresponding excitation energies are computed using high-level coupled-cluster (CC) methods, namely, CC3, CCSDT, CC4, and CCSDTQ, as well as multiconfigurational methods such as CASPT2 and NEVPT2. In many cases, to provide more comprehensive benchmark data, we also provide full configuration interaction estimates computed with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. Based on these calculations, theoretical best estimates of the transition energies are established in both the aug-cc-pVDZ and aug-cc-pVTZ basis sets. This allows us to accurately assess the performance of the CC and multiconfigurational methods for this specific set of challenging transitions. Furthermore, comparisons with experimental data and previous theoretical results are also reported.

Theoretical design and evaluation of efficient small donor molecules for organic solar cells

CONTEXT

The development of high-efficiency photovoltaic devices is the need of time with increasing demand for energy. Herein, we designed seven small molecule donors (SMDs) with A-π-D-π-A backbones containing various acceptor groups for high-efficiency organic solar cells (OSCs). Molecular engineering was performed by substituting the acceptor group in the synthesized compound (BPR) with another highly efficient acceptor group to improve the photoelectric performance of the molecule.

METHOD

The photovoltaic, optoelectronic, and photophysical properties of the proposed compounds (BP1-BP7) were investigated in comparison to BPR using DFT and TD-DFT at MPW1PW91/6-311G(d,p) level of theory. All molecules we designed have red-shifted absorption spectra. The modification of the acceptor fragment of the BPR resulted in a reduced HOMO-LUMO energy gap; thus, the designed compounds (BP1-BP7) had improved optoelectronic responses as compared with the BPR molecule. Various key factors that are crucial for efficient SMDs such as exciton binding energy, frontier molecular orbitals (FMOs), absorption maximum (λmax), open circuit voltage (VOC), dipole moment (μ), excitation charge mobilities, and the transition density matrix of (BPR, BP1-BP7) have also been studied. Low reorganizational energy (holes and electrons) values provide high charge mobility, and all the designed compounds are efficient in this regard. Here, BP6 exhibits low excitation energy (1.66 eV), highest open circuit voltage (2.00 V), normalized VOC (77.23), and fill factor (0.931). Consequently, the superiority of the designed molecules advises experimenters to envision future developments in extremely effective OSC devices.

The three kingdoms-Photoinduced electron transfer cascades controlled by electronic couplings

Excited states are the key species in photocatalysis, while the critical parameters that govern their applications are (i) excitation energy, (ii) accessibility, and (iii) lifetime. However, in molecular transition metal-based photosensitizers, there is a design tension between the creation of long-lived excited (triplet), e.g., metal-to-ligand charge transfer (3MLCT) states and the population of such states. Long-lived triplet states have low spin-orbit coupling (SOC) and hence their population is low. Thus, a long-lived triplet state can be populated but inefficiently. If the SOC is increased, the triplet state population efficiency is improved-coming at the cost of decreasing the lifetime. A promising strategy to isolate the triplet excited state away from the metal after intersystem crossing (ISC) involves the combination of transition metal complex and an organic donor/acceptor group. Here, we elucidate the excited state branching processes in a series of Ru(II)-terpyridyl push-pull triads by quantum chemical simulations. Scalar-relativistic time-dependent density theory simulations reveal that efficient ISC takes place along 1/3MLCT gateway states. Subsequently, competitive electron transfer (ET) pathways involving the organic chromophore, i.e., 10-methylphenothiazinyl and the terpyridyl ligands are available. The kinetics of the underlying ET processes were investigated within the semiclassical Marcus picture and along efficient internal reaction coordinates that connect the respective photoredox intermediates. The key parameter that governs the population transfer away from the metal toward the organic chromophore either by means of ligand-to-ligand (3LLCT; weakly coupled) or intra-ligand charge transfer (3ILCT; strongly coupled) states was determined to be the magnitude of the involved electronic coupling.

Vertical Ionization Potentials and Electron Affinities at the Double-Hybrid Density Functional Level

The double-hybrid (DH) time-dependent density functional theory is extended to vertical ionization potentials (VIPs) and electron affinities (VEAs). Utilizing the density fitting approximation, efficient implementations are presented for the genuine DH ansatz relying on the perturbative second-order correction, while an iterative analogue is also elaborated using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH approach. The favorable computational requirements of the present schemes are discussed in detail. The performance of the recently proposed spin-component-scaled and spin-opposite-scaled (SOS) range-separated (RS) and long-range corrected (LC) DH functionals is comprehensively assessed, while popular hybrid and global DH approaches are also discussed. For the benchmark calculations, up-to-date test sets are selected with high-level coupled-cluster references. Our results show that the ADC(2)-based SOS-RS-PBE-P86 approach is the most accurate and robust functional. This method consistently outperforms the excellent SOS-ADC(2) approach for VIPs, although the results are somewhat less satisfactory for VEAs. Among the genuine DH functionals, the SOS-ωPBEPP86 approach is also recommended for describing ionization processes, but its performance is even less reliable for electron-attached states. In addition, surprisingly good results are attained by the LC hybrid ωB97X-D functional, where the corresponding occupied (unoccupied) orbital energies are retrieved as VIPs (VEAs) within the present formalism.

Enantiopure J-Aggregate of Quaterrylene Bisimides for Strong Chiroptical NIR-Response

Chiral polycyclic aromatic hydrocarbons can be tailored for next-generation photonic materials by carefully designing their molecular as well as supramolecular architectures. Hence, excitonic coupling can boost the chiroptical response in extended aggregates but is still challenging to achieve by pure self-assembly. Whereas most reports on these potential materials cover the UV and visible spectral range, systems in the near infrared (NIR) are underdeveloped. We report a new quaterrylene bisimide derivative with a conformationally stable twisted π-backbone enabled by the sterical congestion of a fourfold bay-arylation. Rendering the π-subplanes accessible by small imide substituents allows for a slip-stacked chiral arrangement by kinetic self-assembly in low polarity solvents. The well dispersed solid-state aggregate reveals a sharp optical signature of strong J-type excitonic coupling in both absorption (897 nm) and emission (912 nm) far in the NIR region and reaches absorption dissymmetry factors up to 1.1 × 10^-2. The structural elucidation was achieved by atomic force microscopy and single-crystal X-ray analysis which we combined to derive a structural model of a fourfold stranded enantiopure superhelix. We could deduce that the role of phenyl substituents is not only granting stable axial chirality but also guiding the chromophore into a chiral supramolecular arrangement needed for strong excitonic chirality.

Performance of Common Density Functionals for Excited States of Tetraphenyldibenzoperiflanthene

Time-dependent density functional theory is the method of choice to efficiently calculate excitation spectra with the functional and basis set choice allowing one to compromise between accuracy and computational cost. In this work, the performance of different functionals as well as the second-order approximate coupled cluster singles and doubles model CC2 is evaluated by comparing the results to experimental results of the example molecule tetraphenyldibenzoperiflanthene (DBP). The choice of the functional has a significant impact on the calculated spectrum of DBP. The performance of a number of different functionals was evaluated, quantified, and, where possible, discussed. The best functional, tuned-CAM-B3LYP, is used to investigate DBP on a surface of hexagonal boron nitride (h-BN). The resulting spectrum shows excellent agreement with experimental results for a monolayer of DBP on h-BN.

Thermal site energy fluctuations in photosystem I: new insights from MD/QM/MM calculations

Cyanobacterial photosystem I (PSI) is one of the most efficient photosynthetic machineries found in nature. Due to the large scale and complexity of the system, the energy transfer mechanism from the antenna complex to the reaction center is still not fully understood. A central element is the accurate evaluation of the individual chlorophyll excitation energies (site energies). Such an evaluation must include a detailed treatment of site specific environmental influences on structural and electrostatic properties, but also their evolution in the temporal domain, because of the dynamic nature of the energy transfer process. In this work, we calculate the site energies of all 96 chlorophylls in a membrane-embedded model of PSI. The employed hybrid QM/MM approach using the multireference DFT/MRCI method in the QM region allows to obtain accurate site energies under explicit consideration of the natural environment. We identify energy traps and barriers in the antenna complex and discuss their implications for energy transfer to the reaction center. Going beyond previous studies, our model also accounts for the molecular dynamics of the full trimeric PSI complex. Via statistical analysis we show that the thermal fluctuations of single chlorophylls prevent the formation of a single prominent energy funnel within the antenna complex. These findings are also supported by a dipole exciton model. We conclude that energy transfer pathways may form only transiently at physiological temperatures, as thermal fluctuations overcome energy barriers. The set of site energies provided in this work sets the stage for theoretical and experimental studies on the highly efficient energy transfer mechanisms in PSI.

Inverse molecular design and parameter optimization with Hückel theory using automatic differentiation

Semiempirical quantum chemistry has recently seen a renaissance with applications in high-throughput virtual screening and machine learning. The simplest semiempirical model still in widespread use in chemistry is Hückel's π-electron molecular orbital theory. In this work, we implemented a Hückel program using differentiable programming with the JAX framework based on limited modifications of a pre-existing NumPy version. The auto-differentiable Hückel code enabled efficient gradient-based optimization of model parameters tuned for excitation energies and molecular polarizabilities, respectively, based on as few as 100 data points from density functional theory simulations. In particular, the facile computation of the polarizability, a second-order derivative, via auto-differentiation shows the potential of differentiable programming to bypass the need for numeric differentiation or derivation of analytical expressions. Finally, we employ gradient-based optimization of atom identity for inverse design of organic electronic materials with targeted orbital energy gaps and polarizabilities. Optimized structures are obtained after as little as 15 iterations using standard gradient-based optimization algorithms.

Double-Hybrid Density Functional Theory for Core Excitations: Theory and Benchmark Calculations

The double-hybrid (DH) time-dependent density functional theory is extended to core excitations. Two different DH formalisms are presented utilizing the core-valence separation (CVS) approximation. First, a CVS-DH variant is introduced relying on the genuine perturbative second-order correction, while an iterative analogue is also presented using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH ansatz. The performance of the new approaches is tested for the most popular DH functionals using the recently proposed XABOOM [J. Chem. Theory Comput.2021, 17, 1618] benchmark set. In order to make a careful comparison, the accuracy and precision of the methods are also inspected. Our results show that the genuine approaches are highly competitive with the more advanced CVS-ADC(2)-based methods if only excitation energies are required. In contrast, as expected, significant differences are observed in oscillator strengths; however, the precision is acceptable for the genuine functionals as well. Concerning the performance of the CVS-DH approaches, the PBE0-2/CVS-ADC(2) functional is superior, while its spin-opposite-scaled variant is also recommended as a cost-effective alternative. For these approaches, significant improvements are realized in the error measures compared with the popular CVS-ADC(2) method.

One-electron self-interaction error and its relationship to geometry and higher orbital occupation

Density Functional Theory (DFT) sees prominent use in computational chemistry and physics; however, problems due to the self-interaction error (SIE) pose additional challenges to obtaining qualitatively correct results. As an unphysical energy an electron exerts on itself, the SIE impacts most practical DFT calculations. We conduct an in-depth analysis of the one-electron SIE in which we replicate delocalization effects for simple geometries. We present a simple visualization of such effects, which may help in future qualitative analysis of the one-electron SIE. By increasing the number of nuclei in a linear arrangement, the SIE increases dramatically. We also show how molecular shape impacts the SIE. Two- and three-dimensional shapes show an even greater SIE stemming mainly from the exchange functional with some error compensation from the one-electron error, which we previously defined [D. R. Lonsdale and L. Goerigk, Phys. Chem. Chem. Phys. 22, 15805 (2020)]. Most tested geometries are affected by the functional error, while some suffer from the density error. For the latter, we establish a potential connection with electrons being unequally delocalized by the DFT methods. We also show how the SIE increases if electrons occupy higher-lying atomic orbitals; seemingly one-electron SIE free methods in a ground are no longer SIE free in excited states, which is an important insight for some popular, non-empirical density functional approximations (DFAs). We conclude that the erratic behavior of the SIE in even the simplest geometries shows that robust DFAs are needed. Our test systems can be used as a future benchmark or contribute toward DFT development.

Excitation energies of polycylic aromatic hydrocarbons by double-hybrid functionals: Assessing the PBE0-DH and PBE-QIDH models and their range-separated versions

A family of non-empirical double-hybrid (DH) density functionals, such as Perdew-Burke-Ernzerhof (PBE)0-DH, PBE-QIDH, and their range-separated exchange (RSX) versions RSX-0DH and RSX-QIDH, all using Perdew-Burke-Ernzerhof(PBE) exchange and correlationfunctionals, is applied here to calculate the excitation energies for increasingly longer linear and cyclic acenes as part of their intense benchmarking for excited states of all types. The energies for the two lowest-lying singlet ¹L_a and ¹L_b states of linear oligoacenes as well as the triplet ³L_a and ³L_b states, are calculated and compared with experimental results. These functionals clearly outperform the results obtained from hybrid functionals and favorably compare with other double-hybrid expressions also tested here, such as B2-PLYP, B2GP-PLYP, ωB2-PLYP, and ωB2GP-PLYP. The study is complemented by the computation of adiabatic S₀-T₁ singlet-triplet energy difference for linear acenes as well as the extension of the study to strained cyclic oligomers, showing how the family of non-empirical expressions robustly leads to competitive results.

Inversion Theory Leveling as a New Methodological Approach to Antioxidant Thermodynamics: A Case Study on Phenol

Antioxidants are various types of compounds that represent a link between biology and chemistry. With the development of theoretical and computational methods, antioxidants are now being studied theoretically. Here, a novel method is presented that aims to reduce the estimated wall times for DFT calculations that result in the same or higher degree of accuracy in the second derivatives over energy than is the case with the regular computational route (i.e., optimizing the reaction system at a lower model and then recalculating the energies at a higher level of theory) by applying the inversion of theory level to the universal chemical scavenger model, i.e., phenol. The resulting accuracy and wall time obtained with such a methodological setup strongly suggest that this methodology could be generally applied to antioxidant thermodynamics for some costly DFT methods with relative absolute deviation.

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.

Double Hybrid Density Functionals for the Electronic Excitation Energies of Linear Cyanines

The lowest bright electronic excitations of seven model linear cyanines (CN3-CN15) using 28 double-hybrid (DH) density functionals are benchmarked against accurate and recent CC3 results. Some of these DH functionals are recently designed specifically for excited electronic state calculations. In addition, CIS, CIS(D), SCS-CIS(D), and SOS-CIS(D) were also tested. Four different basis sets were used for the vertical electronic excitation calculations: cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis. Augmented basis sets (e.g. aug-cc-pVDZ and aug-cc-pVTZ) are found to be required for accurate and consistent results using DH functionals. The DH functionals tested in this work are classified into four main groups: global double-hybrids (GDH), range-separated double-hybrids (RSDH), spin-component and spin-opposite scaling global double-hybrids (SCS/SOS-GDH), and spin-component and spin-opposite scaling range-separated double-hybrids (SCS/SOS-RSDH). Within these groups, the SCS/SOS-RSDH group of functionals is found to provide the lowest mean absolute error (MAE) values (in the range 0.020-0.148 eV) in comparison with the GDH group (0.195-0.441 eV), the RSDH group (0.186-0.511 eV), and the SCS/SOS-GDH group (0.079-0.461 eV). Of all the DH functionals and ab initio methods investigated in the present contribution, the following functionals are found to be the most accurate and consistent: SCS-ωB2GPPLYP (MAE = 0.036 eV), SOS-ωB2GPPLYP (MAE = 0.020 eV), SOS-ωB88PP86 (MAE = 0.035 eV), and SOS-ωPBEPP86 (MAE = 0.037 eV). In general, the ab initio methods tested here show mediocre performance as compared to many DH functionals.

Interexcited State Photophysics I: Benchmarking Density Functionals for Computing Nonadiabatic Couplings and Internal Conversion Rate Constants

We present the first benchmarking study of nonadiabatic matrix coupling elements (NACMEs) calculated using different density functionals. Using the S₁ → S₀ transition in perylene solvated in toluene as a case study, we calculate the photophysical properties and corresponding rate constants for a variety of density functionals from each rung of Jacob's ladder. The singlet photoluminescence quantum yield (sPLQY) is taken as a measure of accuracy, measured experimentally here as 0.955. Important quantum chemical parameters such as geometries, absorption, emission, and adiabatic energies, NACMEs, Hessians, and transition dipole moments were calculated for each density functional basis set combination (data set) using density functional theory based multireference configuration interaction (DFT/MRCI) and compared to experiment where possible. We were able to derive simple relations between the TDDFT and DFT/MRCI photophysical properties; with semiempirical damping factors of ∼0.843 ± 0.017 and ∼0.954 ± 0.064 for TDDFT transition dipole moments and energies to DFT/MRCI level approximations, respectively. NACMEs were dominated by out-of-plane derivative components belonging to the center-most ring atoms with weaker contributions from perturbations along the transverse and longitudinal axes. Calculated theoretical spectra compared well to both experiment and literature, with fluorescence lifetimes between 7.1 and 12.5 ns, agreeing within a factor of 2 with experiment. Internal conversion (IC) rates were then calculated and were found to vary wildly between 10⁶-10¹⁶ s^-1 compared with an experimental rate of the order 10⁷ s^-1. Following further testing by mixing data sets, we found a strong dependence on the method used to obtain the Hessian. The 5 characterized data sets ranked in order of most promising are PBE0/def2-TZVP, ωB97XD/def2-TZVP, HCTH407/TZVP, PBE/TZVP, and PBE/def2-TZVP.

Benefits of Range-Separated Hybrid and Double-Hybrid Functionals for a Large and Diverse Data Set of Reaction Energies and Barrier Heights

To better understand the thermochemical kinetics and mechanism of a specific chemical reaction, an accurate estimation of barrier heights (forward and reverse) and reaction energies is vital. Because of the large size of reactants and transition state structures involved in real-life mechanistic studies (e.g., enzymatically catalyzed reactions), density functional theory remains the workhorse for such calculations. In this paper, we have assessed the performance of 91 density functionals for modeling the reaction energies and barrier heights on a large and chemically diverse data set (BH9) composed of 449 organic chemistry reactions. We have shown that range-separated hybrid functionals perform better than the global hybrids for BH9 barrier heights and reaction energies. Except for the PBE-based range-separated nonempirical double hybrids, range separation of the exchange term helps improve the performance for barrier heights and reaction energies. The 16-parameter Berkeley double hybrid, ωB97M(2), performs remarkably well for both properties. However, our minimally empirical range-separated double hybrid functionals offer marginally better accuracy than ωB97M(2) for BH9 barrier heights and reaction energies.

A Mountaineering Strategy to Excited States: Revising Reference Values with EOM-CC4

In the framework of the computational determination of highly accurate vertical excitation energies in small organic compounds, we explore the possibilities offered by the equation-of-motion formalism relying on the approximate fourth-order coupled-cluster (CC) method, CC4. We demonstrate, using an extended set of more than 200 reference values based on CC including up to quadruples excitations (CCSDTQ), that CC4 is an excellent approximation to CCSDTQ for excited states with a dominant contribution from single excitations with an average deviation as small as 0.003 eV. We next assess the accuracy of several additive basis set correction schemes, in which vertical excitation energies obtained with a compact basis set and a high-order CC method are corrected with lower-order CC calculations performed in a larger basis set. Such strategies are found to be overall very beneficial, though their accuracy depends significantly on the actual scheme. Finally, CC4 is employed to improve several theoretical best estimates of the QUEST database for molecules containing between four and six (nonhydrogen) atoms, for which previous estimates were computed at the CCSDT level.

Can domain-based local pair natural orbitals approaches accurately predict phosphorescence energies?

Since the discovery of the peculiar conducting and optical properties of aromatics, many efforts have been made to characterize and predict their phosphorescence. This physical process is exploited in modern Organic Emitting Light Diodes (OLEDs), and it is also one of the processes decreasing the efficiency of Dye-sensitized solar cells (DSSCs). Herein, we propose a computational strategy for the accurate calculation of singlet-triplet gaps of aromatic compounds, which provides results that are in excellent agreement with available experimental data. Our approach relies on the domain-based local pair natural orbital (DLPNO) variant of the "gold standard" CCSD(T) method. The convergence of our results with respect to the key technical parameters of the calculation, such as the basis set used, the approximations employed in the perturbative triples correction, and the dimension of the PNOs space, was thoroughly discussed.

Assessing Recent Time-Dependent Double-Hybrid Density Functionals on Doublet-Doublet Excitations

This work is the first thorough investigation of time-dependent double-hybrid density functionals (DHDFs) for the calculation of doublet-doublet excitation energies. It sheds light on the current state-of-the-art techniques in the field and clarifies if there is still room for future improvements. Overall, 29 hybrid functionals and DHDFs are investigated. We separately analyze the individual impacts of the Tamm-Dancoff approximation (TDA), range separation, and spin-component/opposite scaling (SCS/SOS) on 45 doublet-doublet excitations in 23 radicals before concluding with an overarching analysis that includes and excludes challenging excitations with double-excitation or multireference character. Our results show again that so-called "nonempirical" DHDFs are outperformed by semiempirical ones. While the best assessed functionals are DHDFs, some of the worst are also DHDFs and outperformed by all assessed hybrids. SCS/SOS is particularly beneficial for range-separated DHDFs. Spin-scaled, range-separated DHDFs paired with the TDA belong to the best tested methods here, and we particularly highlight SCS-ωB2GP-PLYP, SOS-ωB2PLYP, SOS-ωB2GP-PLYP, SOS-ωB88PP86, SOS-RSX-QIDH, and SOS-ωPBEPP86. When comparing our functional rankings with previous studies on singlet-singlet and singlet-triplet excitations, we recommend TDA-SOS-ωB88PP86 and TDA-SOS-ωPBEPP86 as robust methods for excitation energies in general until further improvements have been achieved that surpass the chemical accuracy threshold for challenging open-shell excitations without increasing the computational effort.

Activity-Based Photosensitizers with Optimized Triplet State Characteristics Toward Cancer Cell Selective and Image Guided Photodynamic Therapy

Activity-based theranostic photosensitizers are highly attractive in photodynamic therapy as they offer enhanced therapeutic outcome on cancer cells with an imaging opportunity at the same time. However, photosensitizers (PS) cores that can be easily converted to activity-based photosensitizers (aPSs) are still quite limited in the literature. In this study, we modified the dicyanomethylene-4H-chromene (DCM) core with a heavy iodine atom to get two different PSs (DCM_O-I, I-DCM_O-Cl) that can be further converted to aPS after simple modifications. The effect of iodine positioning on singlet oxygen generation capacity was also evaluated through computational studies. DCM_O-I showed better performance in solution experiments and further proved to be a promising phototheranostic scaffold via cell culture studies. Later, a cysteine (Cys) activatable PS based on the DCM_O-I core (DCM_O-I-Cys) was developed, which induced selective photocytotoxicity along with a fluorescence turn-on response in Cys rich cancer cells.

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

In this paper, the performance of more than 40 popular or recently developed density functionals is assessed for the calculation of 463 vertical excitation energies against the large and accurate QuestDB benchmark set. For this purpose, the Tamm-Dancoff approximation offers a good balance between computational efficiency and accuracy. The functionals ωB97X-D and BMK are found to offer the best performance overall with a root-mean square error (RMSE) of around 0.27 eV, better than the computationally more demanding CIS(D) wave function method with a RMSE of 0.36 eV. The results also suggest that Jacob's ladder still holds for time-dependent density functional theory excitation energies, though hybrid meta generalized-gradient approximations (meta-GGAs) are not generally better than hybrid GGAs. Effects of basis set convergence, gauge invariance correction to meta-GGAs, and nonlocal correlation (VV10) are also studied, and practical basis set recommendations are provided.

Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models

Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for Time-Dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.

Do Double-Hybrid Functionals Benefit from Regularization in the PT2 Term? Observations from an Extensive Benchmark

We put to the test a recent suggestion [Shee, J., et al. J. Phys. Chem. Lett. 2021, 12 (50), 12084-12097] that MP2 regularization might improve the performance of double-hybrid density functionals. Using the very large and chemically diverse GMTKN55 benchmark, we find that κ-regularization is indeed beneficial at lower percentages of Hartree-Fock exchange, especially if spin-component scaling is not applied [such as in B2GP-PLYP or ωB97M(2)]. This benefit dwindles for DSD and DOD functionals and vanishes entirely in the ∼70% HF exchange region optimal for them.

EPR Spectroscopy of Cu(II) Complexes: Prediction of g-Tensors Using Double-Hybrid Density Functional Theory

Computational electron paramagnetic resonance (EPR) spectroscopy is an important field of applied quantum chemistry that contributes greatly to connecting spectroscopic observations with the fundamental description of electronic structure for open-shell molecules. However, not all EPR parameters can be predicted accurately and reliably for all chemical systems. Among transition metal ions, Cu(II) centers in inorganic chemistry and biology, and their associated EPR properties such as hyperfine coupling and g-tensors, pose exceptional difficulties for all levels of quantum chemistry. In the present work, we approach the problem of Cu(II) g-tensor calculations using double-hybrid density functional theory (DHDFT). Using a reference set of 18 structurally and spectroscopically characterized Cu(II) complexes, we evaluate a wide range of modern double-hybrid density functionals (DHDFs) that have not been applied previously to this problem. Our results suggest that the current generation of DHDFs consistently and systematically outperform other computational approaches. The B2GP-PLYP and PBE0-DH functionals are singled out as the best DHDFs on average for the prediction of Cu(II) g-tensors. The performance of the different functionals is discussed and suggestions are made for practical applications and future methodological developments.

Origin of the Failure of Density Functional Theories in Predicting Inverted Singlet-Triplet Gaps

Recent experimental and theoretical studies have shown several new organic molecules that violate Hund’s rule and have the first singlet excited state lower in energy than the first triplet excited state. While many correlated single reference wave function methods have successfully predicted excited-state energetics of these low-lying states, conventional linear-response time-dependent density functional theory (TDDFT) fails to predict the correct excited-state energy ordering. In this article, we have explored the performance of combined DFT and wave function methods like doubles-corrected TDDFT and multiconfiguration pair-density functional theory for the calculation of inverted singlet–triplet gaps. We have also tested the performance of the excited-state DFT (eDFT) method for this problem. Our results have shown that it is possible to obtain inverted singlet–triplet gaps both by using doubles-corrected TDDFT with a proper choice of double-hybrid functionals or by using eDFT.

Charge-Transfer Excitations within Density Functional Theory: How Accurate Are the Most Recommended Approaches?

The performance of the most recent density functionals is assessed for charge-transfer (CT) excitations using comprehensive intra- and intermolecular CT benchmark sets with high-quality reference values. For this comparison, the state-of-the-art range-separated (RS) and long-range-corrected (LC) double hybrid (DH) approaches are selected, and global DH and LC hybrid functionals are also inspected. The correct long-range behavior of the exchange–correlation (XC) energy is extensively studied, and various CT descriptors are compared as well. Our results show that the most robust performance is attained by RS-PBE-P86/SOS-ADC(2), as it is suitable to describe both types of CT excitations with outstanding accuracy. Furthermore, concerning the intramolecular transitions, unexpectedly excellent results are obtained for most of the global DHs, but their limitations are also demonstrated for bimolecular complexes. Despite the outstanding performance of the LC-DH methods for common intramolecular excitations, serious deficiencies are pointed out for intermolecular CT transitions, and the wrong long-range behavior of the XC energy is revealed. The application of LC hybrids to such transitions is not recommended in any respect.

Theoretical Study of Light-Induced Crosslinking Reaction Between Pyrimidine DNA Bases and Aromatic Amino Acids

Low-lying electronic excited states and their relaxation pathways as well as energetics of the crosslinking reaction between uracil as a model system for pyrimidine-type building blocks of DNA and RNA and benzene as a model system for aromatic groups of tyrosine (Tyr) and phenylalanine (Phe) amino acids have been studied in the framework of density functional theory. The equilibrium geometries of the ground and electronic excited states as well as the crossing points between the potential energy surfaces of the uracil–benzene complex were computed. Based on these results, different relaxation pathways of the electronic excited states that lead to either back to the initial geometry configuration or the dimerization between the six-membered rings of the uracil–benzene complex have been identified, and the energetic conditions for their occurrence are discussed. It can be concluded that the DNA–protein crosslinking reaction can be induced by the external electromagnetic field via the dimerization reaction between the six-membered rings of the uracil–benzene pair at the electronic excited-state level of the complex. In the case of the uracil–phenol complex, the configuration of the cyclic adduct (dimerized) conformation is less likely to be formed.

Accurate Spectral Properties within Double-Hybrid Density Functional Theory: A Spin-Scaled Range-Separated Second-Order Algebraic-Diagrammatic Construction-Based Approach

Our second-order algebraic-diagrammatic construction [ADC(2)]-based double-hybrid (DH) ansatz (J. Chem. Theory Comput. 2019, 15, 4440. DOI: 10.1021/acs.jctc.9b00391) is combined with range-separation techniques. In the present scheme, both the exchange and the correlation contributions are range-separated, while spin-scaling approaches are also applied. The new methods are thoroughly tested for the most popular benchmark sets including 250 singlet and 156 triplet excitations, as well as 80 oscillator strengths. It is demonstrated that the range separation for the correlation contributions is highly recommended for both the genuine and the ADC(2)-based DH approaches. Our results show that the latter scheme slightly but consistently outperforms the former one for single excitation dominated transitions. Furthermore, states with larger fractions of double excitations are assessed as well, and challenging charge-transfer excitations are also discussed, where the recently proposed spin-scaled long-range corrected DHs fail. The suggested iterative fourth-power scaling RS-PBE-P86/SOS-ADC(2) method, using only three adjustable parameters, provides the most robust and accurate excitation energies within the DH theory. In addition, the relative error of the oscillator strengths is reduced by 65% compared to the best genuine DH functionals.

Accurate predictions of the electronic excited states of BODIPY based dye sensitizers using spin-component-scaled double-hybrid functionals: a TD-DFT benchmark study

The vertical excitation energies of 13 BODIPY based dye sensitizers are benchmarked by means of TD-DFT, using 36 functionals from different DFT rungs. Most TD-DFT results were found to overestimate the excitation energies, and show mean absolute error (MAE) values in the range 0.2–0.5 eV. The dispersion-corrected, spin-component-scaled, double-hybrid (DSD) functionals DSD-BLYP and DSD-PBEP86 were found to have the smallest MAE values of 0.083 eV and 0.106 eV, respectively, which is close to the range of average errors found in the more expensive coupled-cluster methods. Moreover, DSD-BLYP and DSD-PBEP86 functionals show excellent consistency and quality of results (standard deviation = 0.048 eV and 0.069 eV respectively). However, the range separated hybrid (RSH) and the range separated double hybrid (RSDH) functionals were found to provide the best predictability (linear determination coefficient R² > 0.97 eV).

The excitation energies of 13 BODIPY dye sensitizers are benchmarked by means of TD-DFT, using 36 functionals. Spin-component-scaled double-hybrid (DSD) functionals are found to show the best performance.

A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Bicyclic Systems

Pursuing our efforts to define highly accurate estimates of the relative energies of excited states in organic molecules, we investigate, with coupled-cluster methods including iterative triples (CC3 and CCSDT), the vertical excitation energies of 10 bicyclic molecules (azulene, benzoxadiazole, benzothiadiazole, diketopyrrolopyrrole, furofuran, phthalazine, pyrrolopyrrole, quinoxaline, tetrathiafulvalene, and thienothiophene). In total, we provide aug-cc-pVTZ reference vertical excitation energies for 91 excited states of these relatively large systems. We use these reference values to benchmark various wave function methods, i.e., CIS(D), EOM-MP2, CC2, CCSD, STEOM-CCSD, CCSD(T)(a)*, CCSDR(3), CCSDT-3, ADC(2), ADC(2.5), and ADC(3), as well as some spin-scaled variants of both CC2 and ADC(2). These results are compared to those obtained previously on smaller molecules. It turns out that while the accuracy of some methods is almost unaffected by system size, e.g., CIS(D) and CC3, the performance of others can significantly deteriorate as the systems grow, e.g., EOM-MP2 and CCSD, whereas others, e.g., ADC(2) and CC2, become more accurate for larger derivatives.

Surprisingly Good Performance of XYG3 Family Functionals Using a Scaled KS-MP3 Correlation

By adding a GLPT3 (third-order Görling-Levy perturbation theory, or KS-MP3) term E₃ to the XYG7 form for a double hybrid, we are able to bring down WTMAD2 (weighted total mean absolute deviation) for the very large and chemically diverse GMTKN55 benchmark to an unprecedented 1.17 kcal/mol, competitive with much costlier composite wave function ab initio approaches. Intriguingly, (a) the introduction of E₃ makes an empirical dispersion correction redundant; (b) generalized gradient approximation (GGA) or meta-GGA semilocal correlation functionals offer no advantage over the local density approximation (LDA) in this framework; (c) if a dispersion correction is retained, then simple Slater exchange leads to no significant loss in accuracy. It is possible to create a six-parameter functional with WTMAD2 = 1.42 that has no post-LDA density functional theory components and no dispersion correction in the final energy.

Assessment of hybrid functionals for singlet and triplet excitations: Why do some local hybrid functionals perform so well for triplet excitation energies?

The performance of various hybrid density functionals is assessed for 105 singlet and 105 corresponding triplet vertical excitation energies from the QUEST database. The overall lowest mean absolute error is obtained with the local hybrid (LH) functional LH12ct-SsirPW92 with individual errors of 0.11 eV (0.11 eV) for singlet (triplet) n → π* excitations and 0.29 eV (0.17 eV) for π → π* excitations. This is slightly better than with the overall best performing global hybrid M06-2X [n → π*: 0.13 eV (0.17 eV), π → π*: 0.30 eV (0.20 eV)], while most other global and range-separated hybrids and some LHs suffer from the "triplet problem" of time-dependent density functional theory. This is exemplified by correlating the errors for singlet and triplet excitations on a state-by-state basis. The excellent performance of LHs based on a common local mixing function, i.e., an LMF constructed from the spin-summed rather than the spin-resolved semilocal quantities, is systematically investigated by the introduction of a spin-channel interpolation scheme that allows us to continuously modulate the fraction of opposite-spin terms used in the LMF. The correlation of triplet and singlet errors is systematically improved for the n → π* excitations when larger fractions of the opposite-spin-channel are used in the LMF, whereas this effect is limited for the π → π* excitations. This strongly supports a previously made hypothesis that attributes the excellent performance of LHs based on a common LMF to cross-spin-channel nondynamical correlation terms.