Construction of first-principles density functional approximations and their applications to materials

Abstract

Kohn-Sham density functional theory is a rigorous formulation of many-electron quantum mechanics which, for practical purposes, requires approximation of one term in its total energy expression: the exchange-correlation energy. This work elucidates systematic methods for constructing approximations to the exchange-correlation energy solely from first-principles physics. We review the constraints that can be built into approximate density functionals, and use thermochemical data to argue that satisfaction of these constraints permits a more general description of electronic matter. Contact with semiclassical physics is made by studying the turning surfaces of Kohn-Sham potentials in solids. Perfect metals and covalently-bound, narrow-gap insulators do not have turning surfaces at equilibrium, but do under expansive strain. Wide-gap insulators, ionic crystals, and layered solids tend to have turning surfaces at equilibrium. Chemical bonds in solids are classified using the turning surface radii of its constituent atoms. Depletion of the charge density, such as near a monovacancy in platinum, is shown to produce a turning surface. Further, this work demonstrates why generalized gradient approximations (GGAs) are often able to describe some properties of sp-bonded narrow-gap insulators well. A Laplacian-level pure-density functional is developed with the goal of describing metallic condensed matter. This functional is derived from the r2SCAN orbital-dependent meta-GGA, and reduces its tendency to over-magnetize ferromagnets; improves its description of the equation of state properties of alkali metals; and improves its description of intermetallic thermodynamics. It is constructed to enforce the fourth-order exchange gradient expansion constraint (not satisfied by r2SCAN), and a few free parameters are fitted to paradigmatic metallic systems: jellium surfaces and closed-shell jellium clusters. Last, we modify an exchange-correlation kernel that describes the density-density response of jellium to better satisfy known frequency sum rules. We also constrain the kernel to reproduce the correlation energies of jellium, and compare it to a wide variety of common kernels in use for linear response, time-dependent density functional theory calculations.