Investigation of the Exchange-Correlation Potentials of Functionals Based on the Adiabatic Connection Interpolation

Adiabatic connection interaction strength interpolation method made accurate for the uniform electron gas

The adiabatic connection interaction strength interpolation (ISI)-like method provides a high-level expression for the correlation energy, being, in principle, exact not only in the weak-interaction limit, where it recovers the second-order Görling-Levy perturbation term, but also in the strong-interaction limit that is described by the strictly correlated electron approach. In this work, we construct a genISI functional made accurate for the uniform electron gas, a solid-state physics paradigm that is a very difficult test for ISI-like correlation functionals. We assess the genISI functional for various jellium spheres with the number of electrons Z ≤ 912 and for the non-relativistic noble atoms with Z ≤ 290. For the jellium clusters, the genISI is remarkably accurate, while for the noble atoms, it shows a good performance, similar to other ISI-like methods. Then, the genISI functional can open the path using the ISI-like method in solid-state calculations.

Semilocal Meta-GGA Exchange-Correlation Approximation from Adiabatic Connection Formalism: Extent and Limitations

The incorporation of a strong-interaction regime within the approximate semilocal exchange-correlation functionals still remains a very challenging task for density functional theory. One of the promising attempts in this direction is the recently proposed adiabatic connection semilocal correlation (ACSC) approach [Constantin, L. A.; Phys. Rev. B 2019, 99, 085117] allowing one to construct the correlation energy functionals by interpolation of the high and low-density limits for the given semilocal approximation. The current study extends the ACSC method to the meta-generalized gradient approximations (meta-GGA) level of theory, providing some new insights in this context. As an example, we construct the correlation energy functional on the basis of the high- and low-density limits of the Tao-Perdew-Staroverov-Scuseria (TPSS) functional. Arose in this way, the TPSS-ACSC functional is one-electron self-interaction free and accurate for the strictly correlated and quasi-two-dimensional regimes. Based on simple examples, we show the advantages and disadvantages of ACSC semilocal functionals and provide some new guidelines for future developments in this context.

Møller-Plesset and Density-Fixed Adiabatic Connections for a Model Diatomic System at Different Correlation Regimes

In recent years, adiabatic connection (AC) interpolations developed within density functional theory (DFT) have been found to provide good performances in the calculation of interaction energies when used with Hartree-Fock (HF) ingredients. The physical and mathematical reasons for such unanticipated performance have been clarified, to some extent, by studying the strong-interaction limit of the Møller-Plesset (MP) AC. In this work, we calculate both the MP and the DFT AC integrand for the asymmetric Hubbard dimer, which allows for a systematic investigation of different correlation regimes by varying two simple parameters in the Hamiltonian: the external potential, Δv, and the interaction strength, U. Notably, we find that, while the DFT AC integrand appears to be convex in the full parameter space, the MP integrand may change curvature twice. Furthermore, we discuss different aspects of the second-order expansion of the correlation energy in each AC, and we demonstrate why the derivative of the λ-dependent density in the MP AC at λ = 0 (i.e., at the HF density) is zero in the model. Concerning the strong-interaction limit of both ACs in the Hubbard dimer setting, we show that the asymptotic value of the MP AC, W∞HF, is lower than (or equal to) its DFT analogue, W∞KS, if the two are compared at a given density, just like in real space. However, we also show that this is not always the case if the two quantities are compared at a given external potential.

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.

Nonadiabatic Exchange-Correlation Potential for Strongly Correlated Materials in the Weak and Strong Interaction Limits

In this work, nonadiabatic exchange-correlation (XC) potentials for time-dependent density-functional theory (TDDFT) for strongly correlated materials are derived in the limits of strong and weak correlations. After summarizing some essentials of the available dynamical mean-field theory (DMFT) XC potentials valid for these systems, we present details of the Sham–Schluter equation approach that we use to obtain, in principle, an exact XC potential from a many-body theory solution for the nonequilibrium electron self-energy. We derive the XC potentials for the one-band Hubbard model in the limits of weak and strong on-site Coulomb repulsion. To test the accuracy of the obtained potentials, we compare the TDDFT results obtained with these potentials with the corresponding nonequilibrium DMFT solution for the one-band Hubbard model and find that the agreement between the solutions is rather good. We also discuss possible directions to obtain a universal XC potential that would be appropriate for the case of intermediate interaction strengths, i.e., a nonadiabatic potential that can be used to perform TDDFT analysis of nonequilibrium phenomena, such as transport and other ultrafast properties of materials with any strength of electron correlation at any value in the applied perturbing field.

Numerically Precise Benchmark of Many-Body Self-Energies on Spherical Atoms

We investigate the performance of beyond-GW approaches in many-body perturbation theory by addressing atoms described within the spherical approximation via a dedicated numerical treatment based on B-splines and spherical harmonics. We consider the GW, second Born (2B), and GW + second order screened exchange (GW+SOSEX) self-energies and use them to obtain ionization potentials from the quasi-particle equation (QPE) solved perturbatively on top of independent-particle calculations. We also solve the linearized Sham–Schlüter equation (LSSE) and compare the resulting xc potentials against exact data. We find that the LSSE provides consistent starting points for the QPE but does not present any practical advantage in the present context. Still, the features of the xc potentials obtained with it shed light on possible strategies for the inclusion of beyond-GW diagrams in the many-body self-energy. Our findings show that solving the QPE with the GW+SOSEX self-energy on top of a PBE or PBE0 solution is a viable scheme to go beyond GW in finite systems, even in the atomic limit. However, GW shows a comparable performance if one agrees to use a hybrid starting point. We also obtain promising results with the 2B self-energy on top of Hartree–Fock, suggesting that the full time-dependent Hartree–Fock vertex may be another viable beyond-GW scheme for finite systems.

Improving Results by Improving Densities: Density-Corrected Density Functional Theory

Density functional theory (DFT) calculations have become widespread in both chemistry and materials, because they usually provide useful accuracy at much lower computational cost than wavefunction-based methods. All practical DFT calculations require an approximation to the unknown exchange-correlation energy, which is then used self-consistently in the Kohn-Sham scheme to produce an approximate energy from an approximate density. Density-corrected DFT is simply the study of the relative contributions to the total energy error. In the vast majority of DFT calculations, the error due to the approximate density is negligible. But with certain classes of functionals applied to certain classes of problems, the density error is sufficiently large as to contribute to the energy noticeably, and its removal leads to much better results. These problems include reaction barriers, torsional barriers involving π-conjugation, halogen bonds, radicals and anions, most stretched bonds, etc. In all such cases, use of a more accurate density significantly improves performance, and often the simple expedient of using the Hartree-Fock density is enough. This Perspective explains what DC-DFT is, where it is likely to improve results, and how DC-DFT can produce more accurate functionals. We also outline challenges and prospects for the field.

Gradient Expansions for the Large-Coupling Strength Limit of the Møller-Plesset Adiabatic Connection

The adiabatic connection that has, as weak-interaction expansion, the Møller-Plesset perturbation series has been recently shown to have a large coupling-strength expansion, in terms of functionals of the Hartree-Fock density with a clear physical meaning. In this work, we accurately evaluate these density functionals and we extract second-order gradient coefficients from the data for neutral atoms, following ideas similar to the ones used in the literature for exchange, with some modifications. These new gradient expansions will be the key ingredient for performing interpolations that have already been shown to reduce dramatically MP2 errors for large noncovalent complexes. As a byproduct, our investigation of neutral atoms with large number of electrons N indicates that the second-order gradient expansion for exchange grows as N log(N) rather than as N, as often reported in the literature.

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.

Kohn-Sham equations with functionals from the strictly-correlated regime: investigation with a spectral renormalization method

We re-adapt a spectral renormalization method, introduced in nonlinear optics, to solve the Kohn-Sham (KS) equations of density functional theory, with a focus on functionals based on the strictly-correlated electrons (SCE) regime, which are particularly challenging to converge. Important aspects of the method are: (i) the eigenvalues and the density are computed simultaneously; (ii) it converges using randomized initial guesses; (iii) easy to implement. Using this method we could converge for the first time the Kohn-Sham equations with functionals that include the next leading term in the strong-interaction limit of density functional theory, the so called zero-point energy (ZPE) functional as well as with an interaction-strength-interpolation functional that includes both the exact SCE and ZPE terms. This work is the first building block for future studies on quantum systems confined in low dimensions with different statistics and long-range repulsions, such as localization properties of fermions and bosons with strong long-range repulsive interactions in the presence of a random external potential.

Modified Interaction-Strength Interpolation Method as an Important Step toward Self-Consistent Calculations

The modified point charge plus continuum (mPC) model [ConstantinL. A.; Phys. Rev. B2019, 99, 085117] solves the important failures of the original counterpart, namely, the divergences when the reduced gradient of the density is large, such as in the tail of the density and in quasi-dimensional density regimes. The mPC allows us to define a modified interaction-strength interpolation (mISI) method inheriting these good features, which are important steps toward the full self-consistent treatment. Here, we provide an assessment of mISI for molecular systems (i.e., considering thermochemistry properties, correlation energies, vertical ionization potentials, and several noncovalent interactions), harmonium atoms, and functional derivatives in the strong-interaction limit. For all our tests, mISI provides a systematic improvement over the original ISI method. Semilocal approximations of the second-order Görling–Levy (GL2) perturbation theory are also considered in the mISI method, showing considerable worsening of the results. Possible further development of mISI is briefly discussed.

Unveiling the Physics Behind Hybrid Functionals

We show that accurate exchange-correlation hybrid functionals give very physically optimized effective-correlation potentials, capable of correctly describing the quantum oscillations of atoms and molecules. Based on this analysis and on understanding the error cancellation between semilocal exchange and correlation functionals, we propose a very simple, semilocal correlation potential model compatible with the exact exchange of density functional theory, which performs remarkably well for charge densities and orbital energies.

Kinetic Correlation Functionals from the Entropic Regularization of the Strictly Correlated Electrons Problem

In this work, we study the entropic regularization of the strictly correlated electrons formalism, discussing the implications for density functional theory and establishing a link with earlier works on quantum kinetic energy and classical entropy. We carry out a very preliminary investigation (using simplified models) on the use of the solution of the entropic regularized problem to build approximations for the kinetic correlation functional at large coupling strengths. We also analyze lower and upper bounds to the Hohenberg-Kohn functional using the entropic regularized strictly correlated electrons problem.

Adventures in DFT by a wavefunction theorist

The attraction density functional theory (DFT) has for electronic structure theory is that it is easier to do computationally than ab initio, correlated wavefunction methods, due to its effective one-particle structure. On the contrary, ab initio theorists insist on the ability to converge to the right answer in appropriate limits, but this requires a treatment of the reduced two-particle density matrix. DFT avoids that by appealing to an "existence" theorem (not a constructive one) that all its effects are subsummed into a DFT functional of the one-particle density. However, the existence of thousands of DFT functionals emphasizes that there is no satisfactory way to systematically improve the Kohn-Sham (KS) version as most changes in parameterization or formulation seldom lead to a new functional that is genuinely better than others. Some researchers in the DFT community try to address this issue by imposing conditions rigorously derived from exact DFT considerations, but to date, no one has shown how this route will ever lead to converged results even for the ground state, much less for all the other electronic states obtained from time-dependent DFT that are critically important for chemistry. On the contrary, coupled-cluster (CC) theory and its equation-of-motion extensions provide rigorous results for both that KS-DFT methods are attempting to emulate. How to use them and their exact formal properties to tie CC theory to an effective one-particle form is the target of this perspective. This route addresses the devil's triangle of KS-DFT problems: the one-particle spectrum, self-interaction, and the integer discontinuity.