Derivative discontinuities in density functional theory

Machine learning the derivative discontinuity of density-functional theory

Machine learning is a powerful tool to design accurate, highly non-local, exchange-correlation functionals for density functional theory. So far, most of those machine learned functionals are trained for systems with an integer number of particles. As such, they are unable to reproduce some crucial and fundamental aspects, such as the explicit dependency of the functionals on the particle number or the infamous derivative discontinuity at integer particle numbers. Here we propose a solution to these problems by training a neural network as the universal functional of density-functional theory that (a) depends explicitly on the number of particles with a piece-wise linearity between the integer numbers and (b) reproduces the derivative discontinuity of the exchange-correlation energy. This is achieved by using an ensemble formalism, a training set containing fractional densities, and an explicitly discontinuous formulation.

From Kohn-Sham to Many-Electron Energies via Step Structures in the Exchange-Correlation Potential

Accurately describing excited states within Kohn–Sham (KS) density functional theory (DFT), particularly those which induce ionization and charge transfer, remains a great challenge. Common exchange-correlation (xc) approximations are unreliable for excited states owing, in part, to the absence of a derivative discontinuity in the xc energy (Δ), which relates a many-electron energy difference to the corresponding KS energy difference. We demonstrate, analytically and numerically, how the relationship between KS and many-electron energies leads to the step structures observed in the exact xc potential in four scenarios: electron addition, molecular dissociation, excitation of a finite system, and charge transfer. We further show that steps in the potential can be obtained also with common xc approximations, as simple as the LDA, when addressed from the ensemble perspective. The article therefore highlights how capturing the relationship between KS and many-electron energies with advanced xc approximations is crucial for accurately calculating excitations, as well as the ground-state density and energy of systems which consist of distinct subsystems.

Generalized Gradient Approximation Exchange Energy Functional with Near-Best Semilocal Performance

We develop and validate a nonempirical generalized gradient approximation (GGA) exchange (X) density functional that performs as well as the SCAN (strongly constrained and appropriately normed) meta-GGA on standard thermochemistry tests. Additionally, the new functional (NCAP, nearly correct asymptotic potential) yields Kohn-Sham eigenvalues that are useful approximations of the density functional theory (DFT) ionization potential theorem values by inclusion of a systematic derivative discontinuity shift of the X potential. NCAP also enables time-dependent DFT (TD-DFT) calculations of good-quality polarizabilities, hyper-polarizabilities, and one-Fermion excited states without modification (calculated or ad hoc) of the long-range behavior of the exchange potential or other patches. NCAP is constructed by reconsidering the imposition of the asymptotic correctness of the X potential (-1/ r) as a constraint. Inclusion of derivative discontinuity and approximate integer self-interaction correction treatments along with first-principles determination of the effective second-order gradient expansion coefficient yields a major advance over our earlier correct asymptotic potential functional [CAP; J. Chem. Phys. 2015 , 142 , 054105 ]. The new functional reduces a spurious bump in the CAP atomic exchange potential and moves it to distances irrelevantly far from the nucleus (outside the tail of essentially all practical basis functions). It therefore has nearly correct atomic exchange-potential behavior out to rather large finite distances r from the nucleus but eventually goes as - c/ r with an estimated value for the constant c of around 0.3, so as to achieve other important properties of exact DFT exchange within the restrictions of the GGA form. We illustrate the results with the Ne atom optimized effective potentials and with standard molecular benchmark test data sets for thermochemical, structural, and response properties.

How Interatomic Steps in the Exact Kohn-Sham Potential Relate to Derivative Discontinuities of the Energy

Accurate density functional calculations hinge on reliable approximations to the unknown exchange-correlation (xc) potential. The most popular approximations usually lack features of the exact xc potential that are important for an accurate prediction of the fundamental gap and the distribution of charge in complex systems. Two principal features in this regard are the spatially uniform shift in the potential, as the number of electrons infinitesimally surpasses an integer, and the spatial steps that form, for example, between the atoms of stretched molecules. Although both aforementioned concepts are well known, the exact relationship between them remained unclear. Here we establish this relationship via an analytical derivation. We support our result by numerically solving the many-electron Schrödinger equation to extract the exact Kohn-Sham potential and directly observe its features. Spatial steps in the exact xc potential of a full configuration-interaction (FCI) calculation of a molecule are presented in three dimensions.

Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems

Real-space grids are a powerful alternative for the simulation of electronic systems. One of the main advantages of the approach is the flexibility and simplicity of working directly in real space where the different fields are discretized on a grid, combined with competitive numerical performance and great potential for parallelization. These properties constitute a great advantage at the time of implementing and testing new physical models. Based on our experience with the Octopus code, in this article we discuss how the real-space approach has allowed for the recent development of new ideas for the simulation of electronic systems. Among these applications are approaches to calculate response properties, modeling of photoemission, optimal control of quantum systems, simulation of plasmonic systems, and the exact solution of the Schrödinger equation for low-dimensionality systems.