Self-consistent Kohn-Sham method based on the adiabatic-connection fluctuation-dissipation theorem and the exact-exchange kernel

Basis Set Requirements of σ-Functionals for Gaussian- and Slater-Type Basis Functions and Comparison with Range-Separated Hybrid and Double Hybrid Functionals

σ-Functionals belong to the class of Kohn-Sham (KS) correlation functionals based on the adiabatic-connection fluctuation-dissipation theorem and are technically closely related to the random phase approximation (RPA). They have the same computational demand as the latter, with the computational effort of an energy evaluation for both methods being lower than that of a preceding hybrid DFT calculation for typical systems but yield much higher accuracy, reaching chemical accuracy of 1 kcal/mol for quantities such as reactions and transition energies in main group chemistry. In previous work on σ-functionals, rather large Gaussian basis sets have been used. Here, we investigate the actual basis set requirements of σ-functionals and present three setups that employ smaller Gaussian basis sets ranging from quadruple-ζ (QZ) to triple-ζ (TZ) quality and represent a good compromise between accuracy and computational efficiency. Furthermore, we introduce an implementation of σ-functionals based on Slater-type basis sets and present two setups of QZ and TZ quality for this implementation. We test the accuracy of these setups on a large database of various physical properties and types of reactions, as well as equilibrium geometries and vibrational frequencies. As expected, the accuracy of σ-functional calculations becomes somewhat lower with a decreasing basis set size. However, for all setups considered here, calculations with σ-functionals are clearly more accurate than those within the RPA and even more so than those of the conventional KS methods. For the smallest setup using Gaussian-type basis functions and Slater-type basis functions, we introduce a reparametrization that reduces the loss in accuracy due to the basis set error to some extent. A comparison with the range-separated hybrid ωB97X-V and the double hybrid DSD-BLYP-D3 shows that σ functionals outperform in accuracy both of these accurate and, for their class, representative functionals.

Post-Kohn-Sham Random-Phase Approximation and Correction Terms in the Expectation-Value Coupled-Cluster Formulation

Using expectation-value coupled-cluster theory and many-body perturbation theory (MBPT), we formulate a series of corrections to the post-Kohn-Sham (post-KS) random-phase approximation (RPA) energy. The beyond-RPA terms are of two types: those accounting for the non-Hartree-Fock reference and those introducing the coupled-cluster doubles non-ring contractions. The contributions of the former type, introduced via the semicanonical orbital basis, drastically reduce the binding strength in noncovalent systems. The good accuracy is recovered by the attractive third-order doubles correction referred to as Ec2g. The existing RPA approaches based on KS orbitals neglect most of the proposed corrections but can perform well thanks to error cancellation. The proposed method accounts for every contribution in the state-of-the-art renormalized second-order perturbation theory (rPT2) approach but adds additional terms which initially contribute in the third order of MBPT. The cost of energy evaluation scales as noniterative O ( N 4 ) in the implementation with low-rank tensor decomposition. The numerical tests of the proposed approach demonstrate accurate results for noncovalent dimers of polar molecules and for the challenging many-body noncovalent cluster of CH4···(H2O)20.

Strengths and limitations of the adiabatic exact-exchange kernel for total energy calculations

We investigate the adiabatic approximation to the exact-exchange kernel for calculating correlation energies within the adiabatic-connection fluctuation-dissipation framework of time-dependent density functional theory. A numerical study is performed on a set of systems having bonds of different character (H2 and N2 molecules, H-chain, H2-dimer, solid-Ar, and the H2O-dimer). We find that the adiabatic kernel can be sufficient in strongly bound covalent systems, yielding similar bond lengths and binding energies. However, for non-covalent systems, the adiabatic kernel introduces significant errors around equilibrium geometry, systematically overestimating the interaction energy. The origin of this behavior is investigated by studying a model dimer composed of one-dimensional, closed-shell atoms, interacting via soft-Coulomb potentials. The kernel is shown to exhibit a strong frequency dependence at small to intermediate atomic separation that affects both the low-energy spectrum and the exchange-correlation hole obtained from the corresponding diagonal of the two-particle density matrix.

Assessment of the Second-Order Statically Screened Exchange Correction to the Random Phase Approximation for Correlation Energies

With increasing interelectronic distance, the screening of the electron–electron interaction by the presence of other electrons becomes the dominant source of electron correlation. This effect is described by the random phase approximation (RPA) which is therefore a promising method for the calculation of weak interactions. The success of the RPA relies on the cancellation of errors, which can be traced back to the violation of the crossing symmetry of the 4-point vertex, leading to strongly overestimated total correlation energies. By the addition of second-order screened exchange (SOSEX) to the correlation energy, this issue is substantially reduced. In the adiabatic connection (AC) SOSEX formalism, one of the two electron–electron interaction lines in the second-order exchange term is dynamically screened (SOSEX(W, v_c)). A related SOSEX expression in which both electron–electron interaction lines are statically screened (SOSEX(W(0), W(0))) is obtained from the G3W2 contribution to the electronic self-energy. In contrast to SOSEX(W, v_c), the evaluation of this correlation energy expression does not require an expensive numerical frequency integration and is therefore advantageous from a computational perspective. We compare the accuracy of the statically screened variant to RPA and RPA+SOSEX(W, v_c) for a wide range of chemical reactions. While both methods fail for barrier heights, SOSEX(W(0), W(0)) agrees very well with SOSEX(W, v_c) for charged excitations and noncovalent interactions where they lead to major improvements over RPA.

Scaled σ-functionals for the Kohn-Sham correlation energy with scaling functions from the homogeneous electron gas

The recently introduced σ-functionals constitute a new type of functionals for the Kohn-Sham (KS) correlation energy. σ-Functionals are based on the adiabatic-connection fluctuation-dissipation theorem, are computationally closely related to the well-known direct random phase approximation (dRPA), and are formally rooted in many-body perturbation theory along the adiabatic connection. In σ-functionals, the function of the eigenvalues σ of the Kohn-Sham response matrix that enters the coupling constant and frequency integration in the dRPA is replaced by another function optimized with the help of reference sets of atomization, reaction, transition state, and non-covalent interaction energies. σ-Functionals are highly accurate and yield chemical accuracy of 1 kcal/mol in reaction or transition state energies, in main group chemistry. A shortcoming of σ-functionals is their inability to accurately describe processes involving a change of the electron number, such as ionizations or electron attachments. This problem is attributed to unphysical self-interactions caused by the neglect of the exchange kernel in the dRPA and σ-functionals. Here, we tackle this problem by introducing a frequency- and σ-dependent scaling of the eigenvalues σ of the KS response function that models the effect of the exchange kernel. The scaling factors are determined with the help of the homogeneous electron gas. The resulting scaled σ-functionals retain the accuracy of their unscaled parent functionals but in addition yield very accurate ionization potentials and electron affinities. Moreover, atomization and total energies are found to be exceptionally accurate. Scaled σ-functionals are computationally highly efficient like their unscaled counterparts.

Self-Consistent Implementation of Kohn-Sham Adiabatic Connection Models with Improved Treatment of the Strong-Interaction Limit

Adiabatic connection models (ACMs), which interpolate between the limits of weak and strong interaction, are powerful tools to build accurate exchange–correlation functionals. If the exact weak-interaction expansion from the second-order perturbation theory is included, a self-consistent implementation of these functionals is challenging and still absent in the literature. In this work, we fill this gap by presenting a fully self-consistent-field (SCF) implementation of some popular ACM functionals. While using second-order perturbation theory at weak interactions, we have also introduced new generalized gradient approximations (GGAs), beyond the usual point-charge-plus-continuum model, for the first two leading terms at strong interactions, which are crucial to ensure robustness and reliability. We then assess the SCF–ACM functionals for molecular systems and for prototypical strong-correlation problems. We find that they perform well for both the total energy and the electronic density and that the impact of SCF orbitals is directly connected to the accuracy of the ACM functional form. For the H₂ dissociation, the SCF–ACM functionals yield significant improvements with respect to standard functionals also thanks to the use of the new GGAs for the strong-coupling functionals.

Chemical accuracy with σ-functionals for the Kohn-Sham correlation energy optimized for different input orbitals and eigenvalues

Recently, a new type of orbital-dependent functional for the Kohn-Sham (KS) correlation energy, σ-functionals, was introduced. Technically, σ-functionals are closely related to the well-known direct random phase approximation (dRPA). Within the dRPA, a function of the eigenvalues σ of the frequency-dependent KS response function is integrated over purely imaginary frequencies. In σ-functionals, this function is replaced by one that is optimized with respect to reference sets of atomization, reaction, transition state, and non-covalent interaction energies. The previously introduced σ-functional uses input orbitals and eigenvalues from KS calculations with the generalized gradient approximation (GGA) exchange-correlation functional of Perdew, Burke, and Ernzerhof (PBE). Here, σ-functionals using input orbitals and eigenvalues from the meta-GGA TPSS and the hybrid-functionals PBE0 and B3LYP are presented and tested. The number of reference sets taken into account in the optimization of the σ-functionals is larger than in the first PBE based σ-functional and includes sets with 3d-transition metal compounds. Therefore, also a reparameterized PBE based σ-functional is introduced. The σ-functionals based on PBE0 and B3LYP orbitals and eigenvalues reach chemical accuracy for main group chemistry. For the 10 966 reactions from the highly accurate W4-11RE reference set, the B3LYP based σ-functional exhibits a mean average deviation of 1.03 kcal/mol compared to 1.08 kcal/mol for the coupled cluster singles doubles perturbative triples method if the same valence quadruple zeta basis set is used. For 3d-transition metal chemistry, accuracies of about 2 kcal/mol are reached. The computational effort for the post-self-consistent evaluation of the σ-functional is lower than that of a preceding PBE0 or B3LYP calculation for typical systems.

Numerically stable optimized effective potential method with standard Gaussian basis sets

We present a numerically stable optimized effective potential (OEP) method based on Gaussian basis sets. The key point of the approach is a sequence of preprocessing steps of the auxiliary basis set used to represent exchange or correlation potentials, the Kohn-Sham (KS) response function, and the right-hand side of the OEP equation in conjunction with a representation of exchange or correlation potentials via exchange or correlation charge densities whose electrostatic potentials generate the potentials. Due to the preprocessing, standard Gaussian basis sets from basis set libraries can be used in OEP calculations. As examples, we present numerical stable computational setups based on aux-cc-pwCVXZ basis sets with X = T, Q, 5 for the orbitals and aux-cc-pVDZ/mp2fit and aux-cc-pVTZ/mp2fit auxiliary basis sets and use them to calculate KS exchange potentials with the exact exchange-only KS method for various atoms and molecules. The resulting exchange potentials not only are numerically stable and physically reasonable but also show convergence with increasing quality of the orbital basis sets. The effect of incorporating exact conditions that the KS exchange potential has to obey is discussed. Moreover, it is briefly demonstrated that the presented approach not only works for KS exchange potentials but equally well for correlation potentials within the direct random phase approximation. Besides for OEP methods, the introduced preprocessing of auxiliary basis sets should also be beneficial in procedures to calculate back effective KS potentials from given electron densities.

Optimized effective potentials from the random-phase approximation: Accuracy of the quasiparticle approximation

The optimized effective potential (OEP) method presents an unambiguous way to construct the Kohn-Sham potential corresponding to a given diagrammatic approximation for the exchange-correlation functional. The OEP from the random-phase approximation (RPA) has played an important role ever since the conception of the OEP formalism. However, the solution of the OEP equation is computationally fairly expensive and has to be done in a self-consistent way. So far, large scale solid state applications have, therefore, been performed only using the quasiparticle approximation (QPA), neglecting certain dynamical screening effects. We obtain the exact RPA-OEP for 15 semiconductors and insulators by direct solution of the linearized Sham-Schlüter equation. We investigate the accuracy of the QPA on Kohn-Sham bandgaps and dielectric constants, and comment on the issue of self-consistency.

Toward chemical accuracy at low computational cost: Density-functional theory with σ-functionals for the correlation energy

We introduce new functionals for the Kohn-Sham correlation energy that are based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem and are named σ-functionals. Like in the well-established direct random phase approximation (dRPA), σ-functionals require as input exclusively eigenvalues σ of the frequency-dependent KS response function. In the new functionals, functions of σ replace the σ-dependent dRPA expression in the coupling-constant and frequency integrations contained in the ACFD theorem. We optimize σ-functionals with the help of reference sets for atomization, reaction, transition state, and non-covalent interaction energies. The optimized functionals are to be used in a post-self-consistent way using orbitals and eigenvalues from conventional Kohn-Sham calculations employing the exchange-correlation functional of Perdew, Burke, and Ernzerhof. The accuracy of the presented approach is much higher than that of dRPA methods and is comparable to that of high-level wave function methods. Reaction and transition state energies from σ-functionals exhibit accuracies close to 1 kcal/mol and thus approach chemical accuracy. For the 10 966 reactions of the W4-11RE reference set, the mean absolute deviation is 1.25 kcal/mol compared to 3.21 kcal/mol in the dRPA case. Non-covalent binding energies are accurate to a few tenths of a kcal/mol. The presented approach is highly efficient, and the post-self-consistent calculation of the total energy requires less computational time than a density-functional calculation with a hybrid functional and thus can be easily carried out routinely. σ-Functionals can be implemented in any existing dRPA code with negligible programming effort.

Lieb-Oxford bound and pair correlation functions for density-functional methods based on the adiabatic-connection fluctuation-dissipation theorem

Compliance with the Lieb-Oxford bound for the indirect Coulomb energy and for the exchange-correlation energy is investigated for a number of density-functional methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem to treat correlation. Furthermore, the correlation contribution to the pair density resulting from these methods is compared with highly accurate reference values for the helium atom and for the hydrogen molecule at several bond distances. For molecules, the Lieb-Oxford bound is obeyed by all considered methods. For the homogeneous electron gas, it is violated by all methods for low electron densities. The simplest considered ACFD method, the direct random phase approximation (dRPA), violates the Lieb-Oxford bound much earlier than more advanced ACFD methods that, in addition to the simple Hartree kernel, take into account the exchange kernel and an approximate correlation kernel in the calculation of the correlation energy. While the dRPA yields quite poor correlation contributions to the pair density, those from more advanced ACFD methods are physically reasonable but still leave room for improvements, particularly in the case of the stretched hydrogen molecule.

Efficient and accurate description of adsorption in zeolites

Accurate theoretical methods are needed to correctly describe adsorption on solid surfaces or in porous materials. The random phase approximation (RPA) with singles corrections scheme and the second order Møller-Plesset perturbation theory (MP2) are two schemes, which offer high accuracy at affordable computational cost. However, there is little knowledge about their applicability and reliability for different adsorbates and surfaces. Here, we calculate adsorption energies of seven different molecules in zeolite chabazite to show that RPA with singles corrections is superior to MP2, not only in terms of accuracy but also in terms of computer time. Therefore, RPA with singles is a suitable scheme for obtaining highly accurate adsorption energies in porous materials and similar systems.

Directly patching high-level exchange-correlation potential based on fully determined optimized effective potentials

The key element in Kohn-Sham (KS) density functional theory is the exchange-correlation (XC) potential. We recently proposed the exchange-correlation potential patching (XCPP) method with the aim of directly constructing high-level XC potential in a large system by patching the locally computed, high-level XC potentials throughout the system. In this work, we investigate the patching of the exact exchange (EXX) and the random phase approximation (RPA) correlation potentials. A major challenge of XCPP is that a cluster's XC potential, obtained by solving the optimized effective potential equation, is only determined up to an unknown constant. Without fully determining the clusters' XC potentials, the patched system's XC potential is "uneven" in the real space and may cause non-physical results. Here, we developed a simple method to determine this unknown constant. The performance of XCPP-RPA is investigated on three one-dimensional systems: H₂₀, H₁₀Li₈, and the stretching of the H₁₉-H bond. We investigated two definitions of EXX: (i) the definition based on the adiabatic connection and fluctuation dissipation theorem (ACFDT) and (ii) the Hartree-Fock (HF) definition. With ACFDT-type EXX, effective error cancellations were observed between the patched EXX and the patched RPA correlation potentials. Such error cancellations were absent for the HF-type EXX, which was attributed to the fact that for systems with fractional occupation numbers, the integral of the HF-type EXX hole is not -1. The KS spectra and band gaps from XCPP agree reasonably well with the benchmarks as we make the clusters large.

Power Series Approximation for the Correlation Kernel Leading to Kohn-Sham Methods Combining Accuracy, Computational Efficiency, and General Applicability

A power series approximation for the correlation kernel of time-dependent density-functional theory is presented. Using this approximation in the adiabatic-connection fluctuation-dissipation (ACFD) theorem leads to a new family of Kohn-Sham methods. The new methods yield reaction energies and barriers of unprecedented accuracy and enable a treatment of static (strong) correlation with an accuracy of high-level multireference configuration interaction methods but are single-reference methods allowing for a black-box-like handling of static correlation. The new methods exhibit a better scaling of the computational effort with the system size than rivaling wave-function-based electronic structure methods. Moreover, the new methods do not suffer from the problem of singularities in response functions plaguing previous ACFD methods and therefore are applicable to any type of electronic system.