Researcher profile

John P. Perdew

John P. Perdew contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

Incorporation of density scaling constraint in density functional design via contrastive representation learning

In a data-driven paradigm, machine learning (ML) is the central component for developing accurate and universal exchange-correlation (XC) functionals in density functional theory (DFT). It is well known that XC functionals must satisfy several exact conditions and physical constraints, such as density scaling, spin scaling, and derivative discontinuity. In this work, we demonstrate that contrastive learning is a computationally efficient and flexible method to incorporate a physical constraint in ML-based density functional design. We propose a schematic approach to incorporate the uniform density scaling property of electron density for exchange energies by adopting contrastive representation learning during the pretraining task. The pretrained hidden representation is transferred to the downstream task to predict the exchange energies calculated by DFT. The electron density encoder transferred from the pretraining task based on contrastive learning predicts exchange energies that satisfy the scaling property, while the model trained without using contrastive learning gives poor predictions for the scaling-transformed electron density systems. Furthermore, the model with pretrained encoder gives a satisfactory performance with only small fractions of the whole augmented dataset labeled, comparable to the model trained from scratch using the whole dataset. The results demonstrate that incorporating exact constraints through contrastive learning can enhance the understanding of density-energy mapping using neural network (NN) models with less data labeling, which will be beneficial to generalizing the application of NN-based XC functionals in a wide range of scenarios that are not always available experimentally but theoretically justified. This work represents a viable pathway toward the machine learning design of a universal density functional via representation learning.

preprint2022arXiv

Laplacian-level meta-generalized gradient approximation for solid and liquid metals

We derive and motivate a Laplacian-level, orbital-free meta-generalized-gradient approximation (LL-MGGA) for the exchange-correlation energy, targeting accurate ground-state properties of $sp$ and $sd$ metallic condensed matter, in which the density functional for the exchange-correlation energy is only weakly nonlocal due to perfect long-range screening. Our model for the orbital-free kinetic energy density restores the fourth-order gradient expansion for exchange to the r$^2$SCAN meta-GGA [Furness et al., J. Phys. Chem. Lett. 11, 8208 (2020)], yielding a LL-MGGA we call OFR2. OFR2 matches the accuracy of SCAN for prediction of common lattice constants and improves the equilibrium properties of alkali metals, transition metals, and intermetallics that were degraded relative to the PBE GGA values by both SCAN and r$^2$SCAN. We compare OFR2 to the r$^2$SCAN-L LL-MGGA [D. Mejia-Rodriguez and S.B. Trickey, Phys. Rev. B 102, 121109 (2020)] and show that OFR2 tends to outperform r$^2$SCAN-L for the equilibrium properties of solids, but r$^2$SCAN-L much better describes the atomization energies of molecules than OFR2 does. For best accuracy in molecules and non-metallic condensed matter, we continue to recommend SCAN and r$^2$SCAN. Numerical performance is discussed in detail, and our work provides an outlook to machine learning.

preprint2022arXiv

Predictive Power of the Exact Constraints and Appropriate Norms in Density Functional Theory

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 twenty 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: (1) the satisfaction of more exact constraints and appropriate norms makes a functional more predictive over the immense space of many-electron systems; (2) fitting to bonded systems yields an interpolative DFA that may not extrapolate well to systems unlike those in the fitting set. We discuss how the class of well-described systems has grown along with constraint satisfaction, and the possibilities for future functional development.

preprint2022arXiv

Testing the r$^2$SCAN density functional for the thermodynamic stability of solids with and without a van der Waals correction

A central aim of materials discovery is an accurate and numerically reliable description of thermodynamic properties, such as the enthalpies of formation and decomposition. The r$^2$SCAN revision of the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (meta-GGA) balances numerical stability with high general accuracy. To assess the r$^2$SCAN description of solid-state thermodynamics, we evaluate the formation and decomposition enthalpies, equilibrium volumes, and fundamental bandgaps of more than 1,000 solids using r$^2$SCAN, SCAN, and PBE, as well as two dispersion-corrected variants, SCAN+rVV10 and r$^2$SCAN+rVV10. We show that r$^2$SCAN achieves accuracy comparable to SCAN and often improves upon SCAN's already excellent accuracy. Whereas SCAN+rVV10 is often observed to worsen the formation enthalpies of SCAN, and makes no substantial correction to SCAN's cell volume predictions, r$^2$SCAN+rVV10 predicts marginally less-accurate formation enthalpies than r$^2$SCAN, and slightly more-accurate cell volumes than r$^2$SCAN. The average absolute errors in predicted formation enthalpies are found to decrease by a factor of 1.5 to 2.5 from the GGA level to the meta-GGA level. Smaller decreases in error are observed for decomposition enthalpies. For formation enthalpies r$^2$SCAN improves over SCAN for intermetallic systems. For a few classes of systems -- transition metals, intermetallics, weakly-bound solids, and enthalpies of decomposition into compounds -- GGAs are comparable to meta-GGAs. In total, r$^2$SCAN and r$^2$SCAN+rVV10 can be recommended as stable, general-purpose meta-GGAs for materials discovery.

preprint2022arXiv

The Lieb-Oxford Lower Bounds on the Coulomb Energy, Their Importance to Electron Density Functional Theory, and a Conjectured Tight Bound on Exchange

Lieb and Oxford (1981) derived rigorous lower bounds, in the form of local functionals of the electron density, on the indirect part of the Coulomb repulsion energy. The greatest lower bound for a given electron number N depends monotonically upon N, and the N-> infinity limit is a bound for all N. These bounds have been shown to apply to the exact density functionals for the exchange- and exchange-correlation energies that must be approximated for an accurate and computationally efficient description of atoms, molecules, and solids. A tight bound on the exact exchange energy has been derived therefrom for two-electron ground states, and is conjectured to apply to all spin-unpolarized electronic ground states. Some of these and other exact constraints have been used to construct two generations of non-empirical density functionals beyond the local density approximation: the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), and the strongly constrained and appropriately normed (SCAN) meta-GGA.

preprint2022arXiv

Understanding density driven errors for reaction barrier heights

Delocalization errors, such as charge-transfer and some self-interaction errors, plague computationally-efficient and otherwise-accurate density functional approximations (DFAs). Evaluating a semi-local DFA non-self-consistently on the Hartree-Fock (HF) density is often recommended as a computationally cheap remedy for delocalization errors. For sophisticated meta-GGAs like SCAN, this approach can achieve remarkable accuracy. When this HF-DFT (or DFA@HF) significantly improves over the DFA, it is often presumed that the HF density is more accurate than the self-consistent DFA density. By applying the metrics of density-corrected density functional theory (DFT), we show that HF-DFT works for barrier heights by making a localizing charge transfer error or density over-correction, thereby producing a somewhat-reliable cancellation of density- and functional-driven errors for the energy. A quantitative analysis of the charge transfer errors in a few transition states confirms this trend. We do not have the exact functional and exact densities that are needed to evaluate the exact density- and functional-driven errors for the large BH76 database of barrier heights. Instead, we have identified and used three non-local proxy functionals (the SCAN 50% global hybrid, the range-separated hybrid LC-$ω$PBE, and SCAN-FLOSIC) and their self-consistent densities. These functionals yield reasonably accurate self-consistent barrier heights, and their self-consistent total energies are nearly piecewise linear in fractional electron number - two important points of similarity to the exact functional. We argue that density-driven errors of the energy in a self-consistent density functional calculation are second-order in the density error, and that large density-driven errors arise primarily from incorrect electron transfers over length scales larger than the diameter of an atom.

preprint2022arXiv

Workhorse minimally-empirical dispersion-corrected density functional, with tests for weakly-bound systems: r$^{2}$SCAN+rVV10

SCAN+rVV10 has been demonstrated to be a versatile van der Waals (vdW) density functional that delivers good predictions of both energetic and structural properties for many types of bonding. Recently, the r$^{2}$SCAN functional has been devised as a revised form of SCAN with improved numerical stability. In this work, we refit the rVV10 functional to optimize the r$^{2}$SCAN+rVV10 vdW density functional, and test its performance for molecular interactions and layered materials. Our molecular tests demonstrate that r$^{2}$SCAN+rVV10 outperforms its predecessor SCAN+rVV10 in both efficiency (numerical stability) and accuracy. This good performance is also found in lattice constant predictions. In comparison with benchmark results from higher-level theories or experiments, r$^{2}$SCAN+rVV10 yields excellent interlayer binding energies and phonon dispersions for layered materials.

preprint2021arXiv

Construction of meta-GGA functionals through restoration of exact constraint adherence to regularized SCAN functionals

The SCAN meta-GGA exchange-correlation functional [Phys. Rev. Lett. 115, 036402 (2015)] is constructed as a chemical environment-determined interpolation between two separate energy densities: one describes single orbital electron densities accurately, and another describes slowly-varying densities accurately. To conserve constraints known for the exact exchange-correlation functional, the derivatives of this interpolation vanish in the slowly-varying limit. While theoretically convenient, this choice introduces numerical challenges that degrade the functional's efficiency. We have recently reported a modification to the SCAN functional, termed r$^2$SCAN [J. Phys. Chem. Lett. 11, 8208 (2020)] that introduces two regularizations into SCAN which improve its numerical performance at the expense of not recovering the fourth order term of the slowly-varying density gradient expansion for exchange. Here we show the derivation of a progression of functionals (rSCAN, r++SCAN, r$^2$SCAN, and r$^4$SCAN) with increasing adherence to exact conditions while maintaining a smooth interpolation. The greater smoothness of r$^2$SCAN seems to lead to better general accuracy than the additional exact constraint of SCAN or r$^4$SCAN does.

preprint2021arXiv

First-principles wavevector- and frequency-dependent exchange-correlation kernel for jellium at all densities

We propose a spatially and temporally nonlocal exchange-correlation (xc) kernel for the spin-unpolarized fluid phase of ground-state jellium, for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact xc kernel, to accurately describe the correlation energies of bulk jellium, and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel [A. Ruzsinszky, et al., Phys. Rev. B 101, 245135 (2020)] are maintained in the new revised MCP07 (rMCP07) kernel, while others are added. The revision $f_\mathrm{xc}^\mathrm{rMCP07}(q,ω)$ differs from MCP07 only for non-zero frequencies $ω$. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wavevector $q$ tends to zero, the kernel has a $-4πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.}

preprint2020arXiv

A Step in the Direction of Resolving the Paradox of Perdew-Zunger Self-interaction Correction. II. Gauge Consistency of the Energy Density at Three Levels of Approximation

The Perdew-Zunger(PZ) self-interaction correction (SIC) was designed to correct the one-electron limit of any approximate density functional for the exchange-correlation (xc) energy, while yielding no correction to the exact functional. Unfortunately, it spoils the slowly-varying-in-space limits of the uncorrected approximate functionals, where those functionals are right by construction. The right limits can be restored by locally scaling down the energy density of the PZ SIC in many-electron regions, but then a spurious correction to the exact functional would be found unless the self-Hartree and exact self-xc terms of the PZ SIC energy density were expressed in the same gauge. Only the local density approximation satisfies the same-gauge condition for the energy density, which explains why the recent local-scaling SIC (LSIC) is found here to work excellently for atoms and molecules only with this basic approximation, and not with the more advanced generalized gradient approximations (GGAs) and meta-GGAs, which lose the Hartree gauge via simplifying integrations by parts. The transformation of energy density that achieves the Hartree gauge for the exact xc functional can also be applied to approximate functionals. Doing so leads to a simple scaled-down self-interaction (sdSIC) correction that is typically much more accurate than PZ SIC in tests for many molecular properties (including equilibrium bond lengths). The present work shows unambiguously that the largest errors of PZ SIC applied to standard functionals at three levels of approximation can be removed by restoring their correct slowly-varying-density limits. It also confirms the relevance of these limits to atoms and molecules.

preprint2020arXiv

Accurate and numerically efficient r$^2$SCAN meta-generalized gradient approximation

The recently proposed rSCAN functional [J. Chem. Phys. 150, 161101 (2019)] is a regularized form of the SCAN functional [Phys. Rev. Lett. 115, 036402 (2015)] that improves SCAN's numerical performance at the expense of breaking constraints known from the exact exchange-correlation functional. We construct a new meta-generalized gradient approximation by restoring exact constraint adherence to rSCAN. The resulting functional maintains rSCAN's numerical performance while restoring the transferable accuracy of SCAN.

preprint2020arXiv

Classical turning surfaces in solids: When do they occur, and what do they mean?

Classical turning surfaces of Kohn-Sham potentials, separating classically-allowed regions (CARs) from classically-forbidden regions (CFRs), provide a useful and rigorous approach to understanding many chemical properties of molecules. Here we calculate such surfaces for several paradigmatic solids. Our study of perfect crystals at equilibrium geometries suggests that CFRs are absent in metals, rare in covalent semiconductors, but common in ionic and molecular crystals. A CFR can appear at a monovacancy in a metal. In all materials, CFRs appear or grow as the internuclear distances are uniformly expanded. Calculations with several approximate density functionals and codes confirm these behaviors. A classical picture of conduction suggests that CARs should be connected in metals, and disconnected in wide-gap insulators. This classical picture is confirmed in the limits of extreme uniform compression of the internuclear distances, where all materials become metals without CFRs, and extreme expansion, where all materials become insulators with disconnected and widely-separated CARs around the atoms.

preprint2020arXiv

Constraint-based Wavevector- and Frequency-dependent Exchange-Correlation Kernel of the Uniform Electron Gas

According to time-dependent density functional theory, the exact exchange-correlation kernel f$_{xc}$(n, q, $ω$) determines not only the ground-state energy but also the excited-state energies/lifetimes and time-dependent linear density response of an electron gas of uniform density n $=$ 3/(4$π$r$^3_s$). Here we propose a parametrization of this function based upon the satisfaction of exact constraints. For the static ($ω$ = 0) limit, we modify the model of Constantin and Pitarke at small wavevector q to recover the known second-order gradient expansion, plus other changes. For all frequencies $ω$ at q $=$ 0, we use the model of Gross, Kohn, and Iwamoto. A Cauchy integral extends this model to complex $ω$ and implies the standard Kramers-Kronig relations. A scaling relation permits closed forms for not only the imaginary but also the real part of f$_{xc}$ for real $ω$. We then combine these ingredients by damping out the $ω$ dependence at large q in the same way that the q dependence is damped. Away from q $=$ 0 and $ω$ $=$ 0, the correlation contribution to the kernel becomes dominant over exchange, even at r$_s$ $=$ 4, the valence electron density of metallic sodium. The resulting correlation energy from integration over imaginary $ω$ is essentially exact. The plasmon pole of the density response function is found by analytic continuation of f$_{xc}$ to $ω$ just below the real axis, and the resulting plasmon lifetime first decreases from infinity and then increases as q grows from 0 toward the electron-hole continuum. A static charge-density wave is found for r$_s$ $>$ 69, and shown to be associated with softening of the plasmon mode. The exchange-only version of our static kernel confirms Overhauser's 1968 prediction that correlation enhances the charge-density wave.

preprint2020arXiv

r2SCAN-D4: Dispersion corrected meta-generalized gradient approximation for general chemical applications

We combine a regularized variant of the strongly constrained and appropriately normed semilocal density functional [J. Sun, A. Ruzsinszky, and J. P. Perdew, Phys. Rev. Lett. 115, 036402 (2015)] with the latest generation semi-classical London dispersion correction. The resulting density functional approximation r2SCAN-D4 has the speed of generalized gradient approximations while approaching the accuracy of hybrid functionals for general chemical applications. We demonstrate its numerical robustness in real-life settings and benchmark molecular geometries, general main group and organo-metallic thermochemistry, as well as non-covalent interactions in supramolecular complexes and molecular crystals. Main group and transition metal bond lengths have errors of just 0.8%, which is competitive with hybrid functionals for main group molecules and outperforms them for transition metal complexes. The weighted mean absolute deviation (WTMAD2) on the large GMTKN55 database of chemical properties is exceptionally small at 7.5 kcal/mol. This also holds for metal organic reactions with an MAD of 3.3 kcal/mol. The versatile applicability to organic and metal-organic systems transfers to condensed systems, where lattice energies of molecular crystals are within chemical accuracy (errors <1 kcal/mol).

preprint2020arXiv

Simple hydrogenic estimates for the exchange and correlation energies of atoms and atomic ions, with implications for density functional theory

Exact density functionals for the exchange and correlation energies are approximated in practical calculations for the ground-state electronic structure of a many-electron system. An important exact constraint for the construction of approximations is to recover the correct non-relativistic large-$Z$ expansions for the corresponding energies of neutral atoms with atomic number $Z$ and electron number $N=Z$, which are correct to leading order ($-0.221 Z^{5/3}$ and $-0.021 Z \ln Z$ respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to $E_x(N,Z) \approx -0.354 N^{2/3} Z$ (as known before only for $Z \gg N \gg 1$) and $E_c \approx -0.02 N \ln N$. These asymptotic estimates are most correct for atomic ions with large $N$ and $Z \gg N$, but we find that they are qualitatively and semi-quantitatively correct even for small $N$ and for $N \approx Z$. The large-$N$ asymptotic behavior of the energy is pre-figured in small-$N$ atoms and atomic ions, supporting the argument that widely-predictive approximate density functionals should be designed to recover the correct asymptotics. It is shown that the exact Kohn-Sham correlation energy, when calculated from the pure ground-state wavefunction, should have no contribution proportional to $Z$ in the $Z\to \infty$ limit for any fixed $N$.

preprint2019arXiv

A simple self-interaction correction to RPA-like correlation energies

The random phase approximation (RPA) is exact for the exchange energy of a many-electron ground state, but RPA makes the correlation energy too negative by about 0.5 eV/electron. That large short-range error, which tends to cancel out of iso-electronic energy differences, is largely corrected by an exchange-correlation kernel, or (as in RPA+) by an additive local or semilocal correction. RPA+ is by construction exact for the homogeneous electron gas, and it is also accurate for the jellium surface. RPA+ often gives realistic total energies for atoms or solids in which spin-polarization corrections are absent or small. RPA and RPA+ also yield realistic singlet binding energy curves for H2 and N2, and thus RPA+ yields correct total energies even for spin-unpolarized atoms with fractional spins and strong correlation, as in stretched H2 or N2. However, RPA and RPA+ can be very wrong for spin-polarized one-electron systems (especially for stretched H2+), and also for the spin-polarization energies of atoms. The spin-polarization energy is often a small part of the total energy of an atom, but important for ionization energies, electron affinities, and the atomization energies of molecules. Here we propose a computationally efficient generalized RPA+ (gRPA+) that changes RPA+ only for spin-polarized systems by making gRPA+ exact for all one-electron densities, in the same simple semilocal way that the correlation energy densities of many meta-generalized gradent approximations are made self-correlation free. By construction, gRPA+ does not degrade the exact RPA+ description of jellium. gRPA+ is found to greatly improve upon RPA and RPA+ for the ionization energies and electron affinities of light atoms. Many versions of RPA with an approximate exchange-correlation kernel fail to be exact for all one-electron densities, and they can also be self-interaction corrected in this way.

preprint2019arXiv

A step in the direction of resolving the paradox of Perdew-Zunger self-interaction correction

Self-interaction (SI) error, which results when exchange-correlation contributions to the total energy are approximated, limits the reliability of many density functional approximations. The Perdew-Zunger SI correction (PZSIC), when applied in conjunction with the local spin density approximation (LSDA), improves the description of many properties, but overall, this improvement is limited. Here we propose a modification to PZSIC that uses an iso-orbital indicator to identify regions where local SI corrections should be applied. Using this local-scaling SIC (LSIC) approach with LSDA, we analyze predictions for a wide range of properties including, for atoms, total energies, ionization potentials, and electron affinities, and for molecules, atomization energies, dissociation energy curves, reaction energies, and reaction barrier heights. LSIC preserves the results of PZSIC-LSDA for properties where it is successful and provides dramatic improvements for many of the other properties studied. Atomization energies calculated using LSIC are better than those of the Perdew, Burke, and Ernzerhof (PBE) generalized gradient approximation (GGA) and close to those obtained with the Strongly Constrained and Appropriately Normed (SCAN) meta-GGA. LSIC also restores the uniform gas limit for the exchange energy that is lost in PZSIC-LSDA. Further performance improvements may be obtained by an appropriate combination or modification of the local scaling factor and the particular density functional approximation.

preprint2017arXiv

Properties of real metallic surfaces: Effects of density functional semilocality and van der Waals nonlocality

We have computed the surface energies, work functions, and interlayer surface relaxations of clean (111), (110), and (100) surfaces of Al, Cu, Ru, Rh, Pd, Ag, Pt, and Au. Many of these metallic surfaces have technological or catalytic applications. We compare experimental reference values to those of the local density approximation (LDA), the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), the PBEsol (PBE for solids) GGA, the SCAN meta-GGA, and SCAN+rVV10 (SCAN with a long-range van der Waals or vdW correction). The closest agreement with uncertain experimental values is achieved by the simplest density functional (LDA) and by the most sophisticated general-purpose one (SCAN+rVV10). The long-range vdW interaction increases the surface energies by about 10%, and the work functions by about 1%. LDA works for metal surfaces through a stronger-than-usual error cancellation. PBE yields the most-underestimated and presumably least accurate surface energies and work functions. Surface energies within the random phase approximation (RPA) are also reported. Interlayer relaxations from different functionals are in reasonable agreement with one another, and usually with experiment.