Accurate adiabatic connection curve beyond the physical interaction strength

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.

Approximate bounds and temperature dependence of adiabatic connection integrands for the uniform electron gas

Thermal density functional theory is commonly used in simulations of warm dense matter, a highly energetic phase characterized by substantial thermal effects and by correlated electrons demanding quantum mechanical treatment. Methods that account for temperature dependence, such as Mermin–Kohn–Sham finite-temperature density functional theory and free energy density functional theory, are now employed with more regularity and available in many standard code packages. However, approximations from zero-temperature density functional theory are still often used in temperature-dependent simulations using thermally weighted electronic densities as an input to exchange–correlation functional approximations, a practice known to miss temperature-dependent effects in the exchange–correlation free energy of these systems. In this work, the temperature-dependent adiabatic connection is demonstrated and analyzed using a well-known parameterization of the uniform electron gas free energy. Useful tools based on this formalism for analyzing and constraining approximations of the exchange–correlation at zero temperature are leveraged for the finite-temperature case. Inspired by the Lieb–Oxford inequality, which provides a lower bound for the ground-state exchange–correlation energy, bounds for the exchange–correlation at finite temperatures are approximated for various degrees of electronic correlation.

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.

Describing static correlation in bond dissociation by Kohn–Sham density functional theory

We show that density functional theory within the RPA (random phase approximation for the exchange-correlation energy) provides a correct description of bond dissociation in H(2) in a spin-restricted Kohn-Sham formalism, i.e., without artificial symmetry breaking. We present accurate adiabatic connection curves both at equilibrium and beyond the Coulson-Fisher point. The strong curvature at large bond length implies important static (left-right) correlation, justifying modern hybrid functional constructions but also demonstrating their limitations. Although exact at infinite separation and accurate near the equilibrium bond length, the RPA dissociation curve displays unphysical repulsion at larger but finite bond lengths. Going beyond the RPA by including the exact exchange kernel (RPA+X), we find a similar repulsion. We argue that this deficiency is due to the absence of double excitations in adiabatic linear response theory. Further analyzing the H(2) dissociation limit we show that the RPA+X is not size consistent, in contrast to the RPA.