Adrienn Ruzsinszky

Restoring the density-gradient expansion for exchange in solids and surfaces.

Popular modern generalized gradient approximations are biased toward the description of free-atom energies. Restoration of the first-principles gradient expansion for exchange over a wide range of density gradients eliminates this bias. We introduce a revised Perdew-Burke-Ernzerhof generalized gradient approximation that improves equilibrium properties of densely packed solids and their surfaces.Successful modern generalized gradient approximations (GGAs) are biased toward atomic energies. Restoration of the first-principles gradient expansion for the exchange energy over a wide range of density gradients eliminates this bias. With many collaborators, I introduce PBEsol, a revised Perdew-Burke-Ernzerhof GGA that improves equilibrium properties of densely-packed solids and their surfaces.

Assessing the performance of recent density functionals for bulk solids

We assess the performance of recent density functionals for the exchange-correlation energy of a nonmolecular solid, by applying accurate calculations with the GAUSSIAN, BAND, and VASP codes to a test set of 24 solid metals and nonmetals. The functionals tested are the modified Perdew-Burke-Ernzerhof generalized gradient approximation PBEsol GGA, the second-order GGA SOGGA, and the Armiento-Mattsson 2005 AM05 GGA. For completeness, we also test more standard functionals: the local density approximation, the original PBE GGA, and the Tao-Perdew-Staroverov-Scuseria meta-GGA. We find that the recent density functionals for solids reach a high accuracy for bulk properties lattice constant and bulk modulus. For the cohesive energy, PBE is better than PBEsol overall, as expected, but PBEsol is actually better for the alkali metals and alkali halides. For fair comparison of calculated and experimental results, we consider the zeropoint phonon and finite-temperature effects ignored by many workers. We show how GAUSSIAN basis sets and inaccurate experimental reference data may affect the rating of the quality of the functionals. The results show that PBEsol and AM05 perform somewhat differently from each other for alkali metal, alkaline-earth metal, and alkali halide crystals where the maximum value of the reduced density gradient is about 2, but perform very similarly for most of the other solids where it is often about 1. Our explanation for this is consistent with the importance of exchange-correlation nonlocality in regions of core-valence overlap.

Accurate first-principles structures and energies of diversely bonded systems from an efficient density functional

One atom or molecule binds to another through various types of bond, the strengths of which range from several meV to several eV. Although some computational methods can provide accurate descriptions of all bond types, those methods are not efficient enough for many studies (for example, large systems, ab initio molecular dynamics and high-throughput searches for functional materials). Here, we show that the recently developed non-empirical strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (meta-GGA) within the density functional theory framework predicts accurate geometries and energies of diversely bonded molecules and materials (including covalent, metallic, ionic, hydrogen and van der Waals bonds). This represents a significant improvement at comparable efficiency over its predecessors, the GGAs that currently dominate materials computation. Often, SCAN matches or improves on the accuracy of a computationally expensive hybrid functional, at almost-GGA cost. SCAN is therefore expected to have a broad impact on chemistry and materials science.

Exchange and correlation in open systems of fluctuating electron number

While the exact total energy of a separated open system varies linearly as a function of average electron number between adjacent integers, the energy predicted by semilocal density-functional approximations is concave up and the exact-exchange-only or Hartree-Fock energy is concave down. As a result, semilocal density functionals fail for separated open systems of fluctuating electron number, as in stretched molecular ions A{sub 2}{sup +} and in solid transition-metal oxides. We develop an exact-exchange theory and an exchange-hole sum rule that explain these failures and we propose a way to correct them via a local hybrid functional.

Scaling down the Perdew-Zunger self-interaction correction in many-electron regions

Semilocal density functional approximations (DFAs) for the exchange-correlation energy suffer from self-interaction error, which is believed to be the cause of many of the failures of common DFAs, such as poor description of charge transfer and transition states of chemical reactions. The standard self-interaction correction (SIC) of Perdew and Zunger mends some of these failures but spoils such essential properties as thermochemistry and equilibrium bond lengths. The Perdew-Zunger SIC seems to overcorrect many-electron systems. In this paper, we propose a modified SIC, which is scaled down in many-electron regions. The new SIC has an improved performance for many molecular properties, including total energies, atomization energies, barrier heights of chemical reactions, ionization potentials, electron affinities, and bond lengths. The local spin-density approximation (LSDA) benefits from SIC more than higher-level functionals do. The scaled-down SIC has only one adjustable parameter. Rationalization of the optimal value of this parameter enables us to construct an almost-nonempirical version of the scaled-down SIC-LSDA, which is significantly better than uncorrected LSDA and even better than the uncorrected generalized gradient approximation. We present an analysis of the formal properties of the scaled-down SIC and define possible directions for further improvements. In particular, we find that exactness for all one-electron densities does not guarantee correct asymptotics for the exchange-correlation potential of a many-electron system.

Density functionals that recognize covalent, metallic, and weak bonds.

Computationally efficient semilocal approximations of density functional theory at the level of the local spin density approximation (LSDA) or generalized gradient approximation (GGA) poorly describe weak interactions. We show improved descriptions for weak bonds (without loss of accuracy for strong ones) from a newly developed semilocal meta-GGA (MGGA), by applying it to molecules, surfaces, and solids. We argue that this improvement comes from using the right MGGA dimensionless ingredient to recognize all types of orbital overlap.

Communication: Self-interaction correction with unitary invariance in density functional theory

Standard spin-density functionals for the exchange-correlation energy of a many-electron ground state make serious self-interaction errors which can be corrected by the Perdew-Zunger self-interaction correction (SIC). We propose a size-extensive construction of SIC orbitals which, unlike earlier constructions, makes SIC computationally efficient, and a true spin-density functional. The SIC orbitals are constructed from a unitary transformation that is explicitly dependent on the non-interacting one-particle density matrix. When this SIC is applied to the local spin-density approximation, improvements are found for the atomization energies of molecules.

Communication: Effect of the orbital-overlap dependence in the meta generalized gradient approximation.

We study for the first time the effect of the dependence of meta generalized gradient approximation (MGGA) for the exchange-correlation energy on its input, the kinetic energy density, through the dimensionless inhomogeneity parameter, α, that characterizes the extent of orbital overlap. This leads to a simple MGGA exchange functional, which interpolates between the single-orbital regime, where α = 0, and the slowly varying density regime, where α ≈ 1, and then extrapolates to α → ∞. When combined with a variant of the Perdew-Burke-Ernzerhof GGA correlation, the resulting MGGA performs equally well for atoms, molecules, surfaces, and solids.

Semilocal density functional obeying a strongly tightened bound for exchange

Significance Efficient calculation of the properties of atoms, molecules, and solids on the computer requires a semilocal approximation to the density functional for the exchange-correlation energy, which becomes thereby a single integral over three-dimensional space. A recent, strongly tightened lower bound on the exchange energy has been built into the approximation “meta-generalized gradient approximations made very simple,” or MGGA-MVS, with accurate results for heats of formation, energy barriers, and weak interactions of molecules, and for lattice constants of solids. This would not have been possible without the use of a third ingredient (the local kinetic energy density) in addition to the standard two (the local electron density and its gradient). This third ingredient permits accurate energies even with the drastically tightened bound. Because of its useful accuracy and efficiency, density functional theory (DFT) is one of the most widely used electronic structure theories in physics, materials science, and chemistry. Only the exchange-correlation energy is unknown, and needs to be approximated in practice. Exact constraints provide useful information about this functional. The local spin-density approximation (LSDA) was the first constraint-based density functional. The Lieb–Oxford lower bound on the exchange-correlation energy for any density is another constraint that plays an important role in the development of generalized gradient approximations (GGAs) and meta-GGAs. Recently, a strongly and optimally tightened lower bound on the exchange energy was proved for one- and two-electron densities, and conjectured for all densities. In this article, we present a realistic “meta-GGA made very simple” (MGGA-MVS) for exchange that respects this optimal bound, which no previous beyond-LSDA approximation satisfies. This constraint might have been expected to worsen predicted thermochemical properties, but in fact they are improved over those of the Perdew–Burke–Ernzerhof GGA, which has nearly the same correlation part. MVS exchange is however radically different from that of other GGAs and meta-GGAs. Its exchange enhancement factor has a very strong dependence upon the orbital kinetic energy density, which permits accurate energies even with the drastically tightened bound. When this nonempirical MVS meta-GGA is hybridized with 25% of exact exchange, the resulting global hybrid gives excellent predictions for atomization energies, reaction barriers, and weak interactions of molecules.

Accurate van der Waals coefficients from density functional theory

The van der Waals interaction is a weak, long-range correlation, arising from quantum electronic charge fluctuations. This interaction affects many properties of materials. A simple and yet accurate estimate of this effect will facilitate computer simulation of complex molecular materials and drug design. Here we develop a fast approach for accurate evaluation of dynamic multipole polarizabilities and van der Waals (vdW) coefficients of all orders from the electron density and static multipole polarizabilities of each atom or other spherical object, without empirical fitting. Our dynamic polarizabilities (dipole, quadrupole, octupole, etc.) are exact in the zero- and high-frequency limits, and exact at all frequencies for a metallic sphere of uniform density. Our theory predicts dynamic multipole polarizabilities in excellent agreement with more expensive many-body methods, and yields therefrom vdW coefficients C6, C8, C10 for atom pairs with a mean absolute relative error of only 3%.