An official website of the United States government
In this work, we introduce the concept of a tunable noninteracting free-energy density functional and present two examples realized: (i) via a simple one-parameter convex combination of two existing functionals and (ii) via the construction of a generalized gradient approximation (GGA) enhancement factor that contains one free parameter and is designed to satisfy a set of incorporated constraints. Functional (i), constructed as a combination of the local Thomas–Fermi and a pseudopotential-adapted GGA for the noninteracting free-energy, has already demonstrated its practical usability for establishing the high temperature end of the equation of state of deuterium [Phys. Rev. B 104, 144104 (2021)] and CHON resin [Phys. Rev. E 106, 045207 (2022)] for inertial confinement fusion applications. Hugoniot calculations for liquid deuterium are given as another example of how the application of computationally efficient orbital-free density functional theory (OF-DFT) can be utilized with the employment of the developed functionals. Once the functionals have been tuned such that the OF-DFT Hugoniot calculation matches the Kohn–Sham solution at some low-temperature point, agreement with the reference Kohn–Sham results for the rest of the high temperature Hugoniot path is very good with relative errors for compression and pressure on the order of 2% or less.
more »
« less
Introductory lecture: when the density of the noninteracting reference system is not the density of the physical system in density functional theory
A major challenge in density functional theory (DFT) is the development of density functional approximations (DFAs) to overcome errors in existing DFAs, leading to more complex functionals. For such functionals, we consider roles of the noninteracting reference systems. The electron density of the Kohn–Sham (KS) reference with a local potential has been traditionally defined as being equal to the electron density of the physical system. This key idea has been applied in two ways: the inverse calculation of such a local KS potential for the reference from a given density and the direct calculation of density and energy based on given DFAs. By construction, the inverse calculation can yield a KS reference with the density equal to the input density of the physical system. In application of DFT, however, it is the direct calculation of density and energy from a DFA that plays a central role. For direct calculations, we find that the self-consistent density of the KS reference defined by the optimized effective potential (OEP), is not the density of the physical system, when the DFA is dependent on the external potential. This inequality holds also for the density of generalized KS (GKS) or generalized OEP reference, which allows a nonlocal potential, when the DFA is dependent on the external potential. Instead, the density of the physical system, consistent with a given DFA, is given by the linear response of the total energy with respect to the variation of the external potential. This is a paradigm shift in DFT on the use of noninteracting references: the noninteracting KS or GKS references represent the explicit computational variables for energy minimization, but not the density of the physical system for external potential-dependent DFAs. We develop the expressions for the electron density so defined through the linear response for general DFAs, demonstrate the results for orbital functionals and for many-body perturbation theory within the second-order and the random-phase approximation, and explore the connections to developments in DFT.
more »
« less
- Award ID(s):
- 1900338
- PAR ID:
- 10258454
- Date Published:
- Journal Name:
- Faraday Discussions
- Volume:
- 224
- ISSN:
- 1359-6640
- Page Range / eLocation ID:
- 9 to 26
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The natural determinant reference (NDR) or principal natural determinant is the Slater determinant comprised of the N most strongly occupied natural orbitals of an N-electron state of interest. Unlike the Kohn–Sham (KS) determinant, which yields the exact ground-state density, the NDR only yields the best idempotent approximation to the interacting one-particle reduced density matrix, but it is well-defined in common atom-centered basis sets and is representation-invariant. We show that the under-determination problem of prior attempts to define a ground-state energy functional of the NDR is overcome in a grand-canonical ensemble framework at the zero-temperature limit. The resulting grand potential functional of the NDR ensemble affords the variational determination of the ground state energy, its NDR (ensemble), and select ionization potentials and electron affinities. The NDR functional theory can be viewed as an “exactification” of orbital optimization and empirical generalized KS methods. NDR functionals depending on the noninteracting Hamiltonian do not require troublesome KS-inversion or optimized effective potentials.more » « less
-
Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. We further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.more » « less
-
Ground-state Kohn-Sham density functional theory provides, in principle, the exact ground-state energy and electronic spin densities of real interacting electrons in a static external potential. In practice, the exact density functional for the exchange-correlation (xc) energy must be approximated in a computationally efficient way. About 20 mathematical properties of the exact xc functional are known. In this work, we review and discuss these known constraints on the xc energy and hole. By analyzing a sequence of increasingly sophisticated density functional approximations (DFAs), we argue that ( a) the satisfaction of more exact constraints and appropriate norms makes a functional more predictive over the immense space of many-electron systems and ( b) fitting to bonded systems yields an interpolative DFA that may not extrapolate well to systems unlike those in the fitting set. We discuss both how the class of well-described systems has grown along with constraint satisfaction and the possibilities for future functional development.more » « less
-
We reformulate the thermally assisted-occupation density functional theory (TAO-DFT) into the Kohn–Sham single-determinant framework and construct two new post-self-consistent field (post-SCF) static correlation correction schemes, named rTAO and rTAO-1. In contrast to the original TAO-DFT with the density in an ensemble form, in which each orbital density is weighted with a fractional occupation number, the ground-state density is given by a single-determinant wavefunction, a regular Kohn–Sham (KS) density, and total ground state energy is expressed in the normal KS form with a static correlation energy formulated in terms of the KS orbitals. In post-SCF calculations with rTAO functionals, an efficient energy scanning to quantitatively determine θ is also proposed. The rTAOs provide a promising method to simulate systems with strong static correlation as original TAO, but simpler and more efficient. We show that both rTAO and rTAO-1 is capable of reproducing most results from TAO-DFT without the additional functional Eθ used in TAO-DFT. Furthermore, our numerical results support that, without the functional Eθ, both rTAO and rTAO-1 can capture correct static correlation profiles in various systems.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/D0FD00102C
