BZPEER
More

Discover

About BZPEERSearchGraph

Research tools

HomeCollectionsResearchCollaborationPublish

Network

ExpertsCommunities

Opportunities

Listings

Account

Log inCreate account
Log inCreate account

Source author record

Aaron D. Kaplan

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.chem-phphysics.comp-phcond-mat.othercond-mat.quant-gascond-mat.softphysics.atm-clus

Catalog footprint

Research overview

11works
7topics
4close collaborators

Actions

Is this your work?

Log in to claim
Open graphBrowse papers

Published work

11 publications

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 spsp and sdsd 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 r2^2SCAN 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 r2^2SCAN. We compare OFR2 to the r2^2SCAN-L LL-MGGA [D. Mejia-Rodriguez and S.B. Trickey, Phys. Rev. B 102, 121109 (2020)] and show that OFR2 tends to outperform r2^2SCAN-L for the equilibrium properties of solids, but r2^2SCAN-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 r2^2SCAN. Numerical performance is discussed in detail, and our work provides an outlook to machine learning.

Open access
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.

Open access
preprint2022arXiv

Testing the r2^2SCAN 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 r2^2SCAN 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 r2^2SCAN description of solid-state thermodynamics, we evaluate the formation and decomposition enthalpies, equilibrium volumes, and fundamental bandgaps of more than 1,000 solids using r2^2SCAN, SCAN, and PBE, as well as two dispersion-corrected variants, SCAN+rVV10 and r2^2SCAN+rVV10. We show that r2^2SCAN 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, r2^2SCAN+rVV10 predicts marginally less-accurate formation enthalpies than r2^2SCAN, and slightly more-accurate cell volumes than r2^2SCAN. 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 r2^2SCAN 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, r2^2SCAN and r2^2SCAN+rVV10 can be recommended as stable, general-purpose meta-GGAs for materials discovery.

Open access
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.

Open access
preprint2022arXiv

Workhorse minimally-empirical dispersion-corrected density functional, with tests for weakly-bound systems: r2^{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 r2^{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 r2^{2}SCAN+rVV10 vdW density functional, and test its performance for molecular interactions and layered materials. Our molecular tests demonstrate that r2^{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, r2^{2}SCAN+rVV10 yields excellent interlayer binding energies and phonon dispersions for layered materials.

Open access
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 r2^2SCAN [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, r2^2SCAN, and r4^4SCAN) with increasing adherence to exact conditions while maintaining a smooth interpolation. The greater smoothness of r2^2SCAN seems to lead to better general accuracy than the additional exact constraint of SCAN or r4^4SCAN does.

Open access
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 fxcrMCP07(q,ω)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 qq tends to zero, the kernel has a −4πα(ω)/q2-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.}

Open access
preprint2020arXiv

Accurate and numerically efficient r2^2SCAN 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.

Open access
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.

Open access
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).

Open access
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-ZZ expansions for the corresponding energies of neutral atoms with atomic number ZZ and electron number N=ZN=Z, which are correct to leading order (−0.221Z5/3-0.221 Z^{5/3} and −0.021Zln⁡Z-0.021 Z \ln Z respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to Ex(N,Z)≈−0.354N2/3ZE_x(N,Z) \approx -0.354 N^{2/3} Z (as known before only for Z≫N≫1Z \gg N \gg 1) and Ec≈−0.02Nln⁡NE_c \approx -0.02 N \ln N. These asymptotic estimates are most correct for atomic ions with large NN and Z≫NZ \gg N, but we find that they are qualitatively and semi-quantitatively correct even for small NN and for N≈ZN \approx Z. The large-NN asymptotic behavior of the energy is pre-figured in small-NN 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 ZZ in the Z→∞Z\to \infty limit for any fixed NN.

Open access

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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