Designing a paradigm for parameter-free double-hybrid density functionals through the adiabatic connection path

A general justification for hybrid functionals in DFT by means of linear response theory

In the present work, resorting to linear response theory, we examine the plausibility of postulating Kohn-Sham (KS)-type equations which contain, by definition, an effective hybrid potential made up by some arbitrary mixture of local and non-local terms. In this way a general justification for the construction of hybrid functionals is provided without resorting to arguments based on the adiabatic connection, the generalized KS theory or the Levy's constrained search (or its variations). In particular, we examine the cases of single-hybrid functionals, derived from non-local exchange and of double-hybrid functionals, emerging from non-local second-order expressions obtained from the KS perturbation theory. A further generalization for higher-order hybrid functionals is also included.

Stability of the polyynic form of C18, C22, C26, and C30 nanorings: a challenge tackled by range-separated double-hybrid density functionals

We calculate the relative energy between the cumulene and polyyne structures of a set of C_4k+2 (k = 4-7) rings (C₁₈, C₂₂, C₂₆, and C₃₀ prompted by the recent synthesis of the cyclo[18]carbon (or simply C₁₈) compounds. Reference results were obtained by a costly Quantum Monte-Carlo (QMC) approach, providing thus very accurate values allowing to systematically compare the performance of a variety of wavefunction methods [(i.e., MP2, SCS-MP2, SOS-MP2, DLPNO-CCSD, and DLPNO-CCSD(T)] as well as DFT approaches, applying for the latter a diversity of density functionals covering global and range-separated hybrid and double-hybrid models. The influence of the use of a range-separation scheme for density functionals, for both hybrid and double-hybrid expressions, is discussed according to its key role. Overall, range-separated double-hybrid functionals (e.g., RSX-QIDH) behave very accurately and provide competitive results compared with DLPNO-CCSD(T), at a more reasonable computational cost.

Spin-Opposite-Scaled Range-Separated Exchange Double-Hybrid Models (SOS-RSX-DHs): Marriage Between DH and RSX/SOS-RSX Is Not Always a Happy Match

The range-separated version of double-hybrid density functional theory (DH-DFT) with a remarkable efficiency for both ground-state and excited-state characteristics has recently come into spotlight. In this work, based on theoretical arguments, several variants of spin-opposite-scaled range-separated exchange double-hybrid models (SOS-RSX-DHs) have been proposed and validated. More specifically, we first extend the RSX-DHs to design some other related models. Next, the SOS version of the resulting approximations is constructed and thoroughly evaluated using standard benchmark compilations of various properties. It is shown that although there are properties for which the RSX-DH and SOS-RSX-DH frameworks are rival, there are still some problems particularly prone to the self-interaction error issues where our proposed models seem to be beneficial. Furthermore, some of the presented models devoid of any additional corrections can also surpass the recently proposed approximations from different rungs of "Jacob's Ladder". Nonetheless, perusing the results of different methods and detailed comparisons with the predecessors discloses that all things may not necessarily be well with the RSX and SOS-RSX schemes, where the parent DHs as well as their SOS counterparts can still come into play.

Generalizing Double-Hybrid Density Functionals: Impact of Higher-Order Perturbation Terms

Connections between the Görling-Levy (GL) perturbation theory and the parameters of double-hybrid (DH) density functional are established via adiabatic connection formalism. Moreover, we present a more general DH density functional theory, where the higher-order perturbation terms beyond the second-order GL2 one, such as GL3 and GL4, also contribute. It is shown that a class of DH functionals including previously proposed ones can be formed using the present construction. Based on the proposed formalism, we assess the performance of higher-order DH and long-range corrected DH formed on the Perdew-Burke-Ernzerhof (PBE) semilocal functional and second-order GL2 correlation energy. The underlying construction of DH functionals based on the generalized many-body perturbation approaches is physically appealing in terms of the development of the non-local forms using more accurate and sophisticated semilocal functionals.

Appraising spin-state energetics in transition metal complexes using double-hybrid models: accountability of SOS0-PBESCAN0-2(a) as a promising paradigm

Through a comprehensive survey, reliable double-hybrid models have been validated and proposed for spin-state energetics in transition metal complexes.

Semi-empirical or non-empirical double-hybrid density functionals: which are more robust?

The development of non-empirical double-hybrid density functionals (DHDFs) is a very active research area with the number of approaches in this field having increased rapidly. At the same time, there is a lack of published work that provides a fair assessment and comparison between non-empirical and semi-empirical DHDFs on an equal footing. Herein, we close this gap and present a thorough analysis of both classes of DHDFs on the large GMTKN55 benchmark database for general main-group thermochemistry, kinetics, and noncovalent interactions [Goerigk et al., Phys. Chem. Chem. Phys., 2017, 19, 32184-32215]. In total, 115 variations of dispersion-corrected and -uncorrected DHDFs are tested, which will be condensed to an in-depth assessment of 31 methods: 19 non-empirical and 12 semi-empirical DHDFs. As such, our study represents the largest DHDF study ever conducted and can serve as an important benchmark informing method developers and users alike. Our results show that semi-empirical DHDFs are the most robust density functional approximations and more reliable and accurate than non-empirical ones. In fact, some non-empirical approaches are even outperformed by hybrid approaches or even dispersion-corrected and -uncorrected MP2 and SCS-MP2. SOS0-PBE0-2-D3(BJ) is the only exception and the only non-empirical DHDF that we can safely recommend for general applicability. However, it is still outperformed by six semi-empirical DHDFs, of which we would like to particularly recommend the following five: ωB97X-2-D3(BJ), DSD-BLYP-D3(BJ), DSD-PBEP86-D3(BJ), B2NC-PLYP-D3(BJ), and B2GPPLYP-D3(BJ). Our findings seriously question current trends in the field and they highlight that novel strategies have to be found in order to outperform the currently best density functional theory methods on the market. We hope that our study can function as an important cornerstone inspiring such a change of direction in the field.

Dissecting the accountability of parameterized and parameter-free single-hybrid and double-hybrid functionals for photophysical properties of TADF-based OLEDs

Organic light emitting diodes (OLEDs) based on thermally activated delayed fluorescence (TADF) emitters are an attractive category of materials that have witnessed a booming development in recent years. In the present contribution, we scrutinize the accountability of parameterized and parameter-free single-hybrid (SH) and double-hybrid (DH) functionals through the two formalisms, full time-dependent density functional theory (TD-DFT) and Tamm-Dancoff approximation (TDA), for the estimation of photophysical properties like absorption energy, emission energy, zero-zero transition energy, and singlet-triplet energy splitting of TADF molecules. According to our detailed analyses on the performance of SHs based on TD-DFT and TDA, the TDA-based parameter-free SH functionals, PBE0 and TPSS0, with one-third of exact-like exchange turned out to be the best performers in comparison to other functionals from various rungs to reproduce the experimental data of the benchmarked set. Such affordable SH approximations can thus be employed to predict and design the TADF molecules with low singlet-triplet energy gaps for OLED applications. From another perspective, considering this point that both the nonlocal exchange and correlation are essential for a more reliable description of large charge-transfer excited states, applicability of the functionals incorporating these terms, namely, parameterized and parameter-free DHs, has also been evaluated. Perusing the role of exact-like exchange, perturbative-like correlation, solvent effects, and other related factors, we find that the parameterized functionals B2π-PLYP and B2GP-PLYP and the parameter-free models PBE-CIDH and PBE-QIDH have respectable performance with respect to others. Lastly, besides the recommendation of reliable computational protocols for the purpose, hopefully this study can pave the way toward further developments of other SHs and DHs for theoretical explorations in the field of OLEDs technology.

Role of exact exchange in thermally-assisted-occupation density functional theory: A proposal of new hybrid schemes

We propose hybrid schemes incorporating exact exchange into thermally assisted-occupation-density functional theory (TAO-DFT) [J.-D. Chai, J. Chem. Phys. 136, 154104 (2012)] for an improved description of nonlocal exchange effects. With a few simple modifications, global and range-separated hybrid functionals in Kohn-Sham density functional theory (KS-DFT) can be combined seamlessly with TAO-DFT. In comparison with global hybrid functionals in KS-DFT, the resulting global hybrid functionals in TAO-DFT yield promising performance for systems with strong static correlation effects (e.g., the dissociation of H₂ and N₂, twisted ethylene, and electronic properties of linear acenes), while maintaining similar performance for systems without strong static correlation effects. Besides, a reasonably accurate description of noncovalent interactions can be efficiently achieved through the inclusion of dispersion corrections in hybrid TAO-DFT. Relative to semilocal density functionals in TAO-DFT, global hybrid functionals in TAO-DFT are generally superior in performance for a wide range of applications, such as thermochemistry, kinetics, reaction energies, and optimized geometries.

SCAN-based hybrid and double-hybrid density functionals from models without fitted parameters

By incorporating the nonempirical strongly constrained and appropriately normed (SCAN) semilocal density functional [J. Sun, A. Ruzsinszky, and J. P. Perdew, Phys. Rev. Lett. 115, 036402 (2015)] in the underlying expression of four existing hybrid and double-hybrid models, we propose one hybrid (SCAN0) and three double-hybrid (SCAN0-DH, SCAN-QIDH, and SCAN0-2) density functionals, which are free from any fitted parameters. The SCAN-based double-hybrid functionals consistently outperform their parent SCAN semilocal functional for self-interaction problems and noncovalent interactions. In particular, SCAN0-2, which includes about 79% of Hartree-Fock exchange and 50% of second-order Møller-Plesset correlation, is shown to be reliably accurate for a very diverse range of applications, such as thermochemistry, kinetics, noncovalent interactions, and self-interaction problems.

Seeking for Spin-Opposite-Scaled Double-Hybrid Models Free of Fitted Parameters

On the basis of theoretical arguments, a new exchange-correlation energy expression free of any fitted parameter has been proposed for spin-opposite-scaled double-hybrid density functionals (SOS0-DHs). Employing the recently presented DHs, the working expressions for SOS0-DH functionals are obtained and benchmarked numerically against several standard databases. Our test calculations show that for some cases such as interaction energies and barrier heights the SOS0-DHs without dispersion corrections perform better than their non-SOS counterparts. On the other hand, for other properties like atomization energies, the conventional DHs provide reliable results. We hope that the findings of this work can excite further developments of DH functionals in the framework of SOS scheme for a wide variety of applications resolving the failures at a reasonable computational cost. It seems that a bright future lies ahead in this arena.

Beyond energies: geometry predictions with the XYG3 type of doubly hybrid density functionals

The scaled mean absolute deviations (s-MADs) of the optimized geometric parameters for covalent bondings (the CCse set), nonbonded interactions (the S22G30 set) and the transition state structures (the TSG36 set), with Tot referring to the averaged s-MAD for general performances.