BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Kieron Burke

Kieron Burke appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.chem-phcond-mat.mtrl-scicond-mat.otherquant-phphysics.comp-phcond-mat.str-elMachine Learningphysics.atom-phcond-mat.mes-hallmath-phmath.MPphysics.atm-cluscond-mat.stat-mech

Catalog footprint

What is connected

70works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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.

Published work

70 published item(s)

preprint2022arXiv

Conditional probability density functional theory

We present conditional probability (CP) density functional theory (DFT) as a formally exact theory. In essence, CP-DFT determines the ground-state energy of a system by finding the CP density from a series of independent Kohn-Sham (KS) DFT calculations. By directly calculating CP densities, we bypass the need for an approximate XC energy functional. In this work we discuss and derive several key properties of the CP density and corresponding CP-KS potential. Illustrative examples are used throughout to help guide the reader through the various concepts and theory presented. We explore a suitable CP-DFT approximation and discuss exact conditions, limitations, and results for selected examples.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Correlation energy of the uniform gas determined by ground state conditional probability density functional theory

Conditional-probability density functional theory (CP-DFT) is a formally exact method for finding correlation energies from Kohn-Sham DFT without evaluating an explicit energy functional. We present details on how to generate accurate exchange-correlation energies for the ground-state uniform gas. We also use the exchange hole in a CP antiparallel spin calculation to extract the high-density limit. We give a highly accurate analytic solution to the Thomas-Fermi model for this problem, showing its performance relative to Kohn-Sham and may be useful at high temperatures. We explore several approximations to the CP potential. Results are compared to accurate parameterizations for both exchange-correlation energies and holes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

How Well Does Kohn-Sham Regularizer Work for Weakly Correlated Systems?

Kohn-Sham regularizer (KSR) is a differentiable machine learning approach to finding the exchange-correlation functional in Kohn-Sham density functional theory (DFT) that works for strongly correlated systems. Here we test KSR for weak correlation. We propose spin-adapted KSR (sKSR) with trainable local, semilocal, and nonlocal approximations found by minimizing density and total energy loss. We assess the atoms-to-molecules generalizability by training on one-dimensional (1D) H, He, Li, Be, Be$^{++}$ and testing on 1D hydrogen chains, LiH, BeH$_2$, and helium hydride complexes. The generalization error from our semilocal approximation is comparable to other differentiable approaches, but our nonlocal functional outperforms any existing machine learning functionals, predicting ground-state energies of test systems with a mean absolute error of 2.7 milli-Hartrees.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Improving results by improving densities: Density-corrected density functional theory

DFT calculations have become widespread in both chemistry and materials, because they usually provide useful accuracy at much lower computational cost than wavefunction-based methods. All practical DFT calculations require an approximation to the unknown exchange-correlation energy, which is then used self-consistently in the Kohn-Sham scheme to produce an approximate energy from an approximate density. Density-corrected DFT is simply the study of the relative contributions to the total energy error. In the vast majority of DFT calculations, the error due to the approximate density is negligible. But with certain classes of functionals applied to certain classes of problems, the density error is sufficiently large as to contribute to the energy noticeably, and its removal leads to much better results. These problems include reaction barriers, torsional barriers involving $π$-conjugation, halogen bonds, radicals and anions, most stretched bonds, etc. In all such cases, use of a more accurate density significantly improves performance, and often the simple expedient of using the Hartree-Fock density is enough. This article explains what DC-DFT is, where it is likely to improve results, and how DC-DFT can produce more accurate functionals. We also outline challenges and prospects for the field.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Lieb's most useful contribution to density functional theory?

The importance of the Lieb-Simon proof of the relative exactness of Thomas-Fermi theory in the large-Z limit to modern density functional theory (DFT) is explored. The principle, that there is a specific semiclassical limit in which functionals become local, implies that there exist well-defined leading functional corrections to local approximations that become relatively exact for the error in local approximations in this limit. It is argued that this principle might be used to greatly improve the accuracy of the thousand or so DFT calculations that are now published each week. A key question is how to find the leading corrections to any local density approximation as this limit is approached. These corrections have been explicitly derived in ridiculously simple model systems to ridiculously high order, yielding ridiculously accurate energies. Much analytic work is needed to use this principle to improve realistic calculations of molecules and solids.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Explaining and Fixing DFT Failures for Torsional Barriers

Most torsional barriers are predicted to high accuracy (about 1kJ/mol) by standard semilocal functionals, but a small subset has been found to have much larger errors. We create a database of almost 300 carbon-carbon torsional barriers, including 12 poorly behaved barriers, all stemming from Y=C-X group, where X is O or S, and Y is a halide. Functionals with enhanced exchange mixing (about 50%) work well for all barriers. We find that poor actors have delocalization errors caused by hyperconjugation. These problematic calculations are density sensitive (i.e., DFT predictions change noticeably with the density), and using HF densities (HF-DFT) fixes these issues. For example, conventional B3LYP performs as accurately as exchange-enhanced functionals if the HF density is used. For long-chain conjugated molecules, HF-DFT can be much better than exchange-enhanced functionals. We suggest that HF-PBE0 has the best overall performance.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bypassing the energy functional in density functional theory: Direct calculation of electronic energies from conditional probability densities

Density functional calculations can fail for want of an accurate exchange-correlation approximation. The energy can instead be extracted from a sequence of density functional calculations of conditional probabilities (CP-DFT). Simple CP approximations yield usefully accurate results for two-electron ions, the hydrogen dimer, and the uniform gas at all temperatures. CP-DFT has no self-interaction error for one electron, and correctly dissociates H2, both major challenges. For warm dense matter, classical CP-DFT calculations can overcome the convergence problems of Kohn-Sham DFT.

0 citations0 reviewsNo code linkedOpen 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Confirmation of the PPLB derivative discontinuity: Exact chemical potential at finite temperatures of a model system

The landmark 1982 paper of Perdew, Parr, Levy, and Balduz (often called PPLB) laid the foundation for our modern understanding of the role of the derivative discontinuity in density functional theory, which drives much development to account for its effects. A simple model for the chemical potential at vanishing temperature played a crucial role in their argument. We investigate the validity of this model in the simplest non-trivial system to which it can be applied and which can be easily solved exactly, the Hubbard dimer. We find exact agreement in the crucial zero-temperature limit, and show the model remains accurate for a significant range of temperatures. We identify how this range depends on the strength of correlations. We extend the model to approximate free energies accounting for the derivative discontinuity, a feature missing in standard semi-local approximations. We provide a correction to this approximation to yield even more accurate free energies. We discuss the relevance of these results for warm dense matter.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Deriving approximate functionals with asymptotics

Modern density functional approximations achieve moderate accuracy at low computational cost for many electronic structure calculations. Some background is given relating the gradient expansion of density functional theory to the WKB expansion in one dimension, and modern approaches to asymptotic expansions. A mathematical framework for analyzing asymptotic behavior for the sums of energies unites both corrections to the gradient expansion of DFT and hyperasymptotics of sums. Simple examples are given for the model problem of orbital-free DFT in one dimension. In some cases, errors can be made as small as 10$^{-32}$ Hartree suggesting that, if these new ingredients can be applied, they might produce approximate functionals that are much more accurate than those in current use. A variation of the Euler-Maclaurin formula generalizes previous results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Kohn-Sham equations as regularizer: building prior knowledge into machine-learned physics

Including prior knowledge is important for effective machine learning models in physics, and is usually achieved by explicitly adding loss terms or constraints on model architectures. Prior knowledge embedded in the physics computation itself rarely draws attention. We show that solving the Kohn-Sham equations when training neural networks for the exchange-correlation functional provides an implicit regularization that greatly improves generalization. Two separations suffice for learning the entire one-dimensional H$_2$ dissociation curve within chemical accuracy, including the strongly correlated region. Our models also generalize to unseen types of molecules and overcome self-interaction error.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

MAP: an MP2 accuracy predictor for weak interactions from adiabatic connection theory

Second order Møller-Plesset perturbation theory (MP2) approximates the exact Hartree-Fock (HF) adiabatic connection (AC) curve by a straight line. Thus by using the deviation of the exact curve from the linear behaviour, we construct an indicator for the accuracy of MP2. We then use an interpolation along the HF AC to transform the exact form of our indicator into a highly practical MP2 accuracy predictor (MAP) that comes at negligible additional computational cost. We show that this indicator is already applicable to systems that dissociate into fragments with a non-degenerate ground state, and we illustrate its usefulness by applying it to the S22 and S66 datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring density-driven errors using Kohn-Sham inversion

Kohn-Sham inversion, that is, the finding of the exact Kohn-Sham potential for a given density, is difficult in localized basis sets. We study the precision and reliability of several inversion schemes, finding estimates of density-driven errors at a useful level of accuracy. In typical cases of substantial density-driven errors, HF-DFT is almost as accurate as DFT evaluated on CCSD(T) densities. A simple approximation in practical HF-DFT also makes errors much smaller than the density-driven errors being calculated. Two paradigm examples, stretched NaCl and the HO$\cdot$Cl$^-$ radical, illustrate just how accurate HF-DFT is.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantifying and understanding errors in molecular geometries

Electronic structure calculations are ubiquitous in most branches of chemistry, but all have errors in both energies and equilibrium geometries. Quantifying errors in possibly dozens of bond angles and bond lengths is a Herculean task. A single natural measure of geometric error is introduced, the geometry energy offset (GEO). GEO links many disparate aspects of geometry errors: a new ranking of different methods, quantitative insight into errors in specific geometric parameters, and insight into trends with different methods. GEO can also reduce the cost of high-level geometry optimizations and shows when geometric errors distort the overall error of a method. Results, including some surprises, are given for both covalent and weak interactions.

0 citations0 reviewsNo code linkedOpen 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-$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$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Explicit corrections to the gradient expansion for the kinetic energy in one dimension

A mathematical framework is constructed for the sum of the lowest N eigenvalues of a potential. Exactness is illustrated on several model systems (harmonic oscillator, particle in a box, and Poschl-Teller well). Its order-by-order semiclassical expansion reduces to the gradient expansion for slowly-varying densities, but also yields a correction when the system is finite and the spectrum discrete. Some singularities can be avoided when evaluating the correction to the leading term. Explicit corrections to the gradient expansion to the kinetic energy in one dimension are found which, in simple cases, greatly improve accuracy. We discuss the relevance to practical density functional calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Can ensemble density functional theory yield accurate double excitations?

The recent use of a new ensemble in density functional theory (DFT) to yield direct corrections to the Kohn-Sham transitions yields the elusive double excitations that are missed by time-dependent DFT with the standard adiabatic approximation. But accuracies are lower than for single excitations, and formal arguments suggest that direct corrections at the exchange level should not be sufficient. We show that in principle, EDFT with direct corrections can yield accurate doubles and explain the error in naive formal arguments about TDDFT. Exact calculations and analytic results on a simple model, the Hubbard dimer, illustrate the results, showing that the answer to the title question is typically yes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Connection formulas for thermal density functional theory

The adiabatic connection formula of ground-state density functional theory relates the correlation energy to a coupling-constant integral over a purely potential contribution, and is widely used to understand and improve approximations. The corresponding formula for thermal density functional theory is cast as an integral over temperatures instead, ranging upward from the system's physical temperature. We also show how to relate different different correlation components to each other, either in terms of temperature- or coupling-constant integrations. We illustrate our results on the uniform electron gas.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Exact conditions on the temperature dependence of density functionals

Universal exact conditions guided the construction of most ground-state density functional approximations in use today. We derive the relation between the entropy and Mermin free energy density functionals for thermal density functional theory. Both the entropy and sum of kinetic and electron-electron repulsion functionals are shown to be monotonically increasing with temperature, while the Mermin functional is concave downwards. Analogous relations are found for both exchange and correlation. The importance of these conditions is illustrated in two extremes: the Hubbard dimer and the uniform gas.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Exact thermal density functional theory for a model system: Correlation components and accuracy of the zero-temperature exchange-correlation approximation

Thermal density functional theory calculations often use the Mermin-Kohn-Sham scheme, but employ ground-state approximations to the exchange-correlation (XC) free energy. In the simplest solvable nontrivial model, an asymmetric Hubbard dimer, we calculate the exact many-body energies and the exact Mermin-Kohn-Sham functionals for this system and extract the exact XC free energy. For moderate temperatures and weak correlation, we find this approximation to be excellent. We extract various exact free-energy correlation components and the exact adiabatic connection formula.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Locality of correlation in density functional theory

The Hohenberg-Kohn density functional was long ago shown to reduce to the Thomas-Fermi approximation in the non-relativistic semiclassical (or large-$Z$) limit for all matter, i.e, the kinetic energy becomes local. Exchange also becomes local in this limit. Numerical data on the correlation energy of atoms supports the conjecture that this is also true for correlation, but much less relevant to atoms. We illustrate how expansions around large particle number are equivalent to local density approximations and their strong relevance to density functional approximations. Analyzing highly accurate atomic correlation energies, we show that the correlation energy tends to $-A_c Z ln Z + B_c Z$ as $Z$ tends to infinity, where $Z$ is the atomic number, $A_c$ is known, and we estimate $B_c$ to be about 37 millihartrees. The local density approximation yields $A_c$ exactly, but a very incorrect value for $B_c$, showing that the local approximation is less relevant for correlation alone. This limit is a benchmark for the non-empirical construction of density functional approximations. We conjecture that, beyond atoms, the leading correction to the local density approximation in the large-$Z$ limit generally takes this form, but with $B_c$ a functional of the TF density for the system. The implications for construction of approximate density functionals are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

One Dimensional Mimicking of Electronic Structure: The Case for Exponentials

An exponential interaction is constructed so that one-dimensional atoms and chains of atoms mimic the general behavior of their three-dimensional counterparts. Relative to the more commonly used soft-Coulomb interaction, the exponential greatly diminishes the computational time needed for calculating highly accurate quantities with the density matrix renormalization group. This is due to the use of a small matrix product operator and to exponentially vanishing tails. Furthermore, its more rapid decay closely mimics the screened Coulomb interaction in three dimensions. Choosing parameters to best match earlier calculations, we report results for the one dimensional hydrogen atom, uniform gas, and small atoms and molecules both exactly and in the local density approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pure density functional for strong correlations and the thermodynamic limit from machine learning

We use density-matrix renormalization group, applied to a one-dimensional model of continuum Hamiltonians, to accurately solve chains of hydrogen atoms of various separations and numbers of atoms. We train and test a machine-learned approximation to $F[n]$, the universal part of the electronic density functional, to within quantum chemical accuracy. Our calculation (a) bypasses the standard Kohn-Sham approach, avoiding the need to find orbitals, (b) includes the strong correlation of highly-stretched bonds without any specific difficulty (unlike all standard DFT approximations) and (c) is so accurate that it can be used to find the energy in the thermodynamic limit to quantum chemical accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Importance of being consistent

We review the role of self-consistency in density functional theory. We apply a recent analysis to both Kohn-Sham and orbital-free DFT, as well as to Partition-DFT, which generalizes all aspects of standard DFT. In each case, the analysis distinguishes between errors in approximate functionals versus errors in the self-consistent density. This yields insights into the origins of many errors in DFT calculations, especially those often attributed to self-interaction or delocalization error. In many classes of problems, errors can be substantially reduced by using `better' densities. We review the history of these approaches, many of their applications, and give simple pedagogical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Thermal density functional theory: Time-dependent linear response and approximate functionals from the fluctuation-dissipation theorem

The van Leeuwen proof of linear-response time-dependent density functional theory (TDDFT) is generalized to thermal ensembles. This allows generalization to finite temperatures of the Gross-Kohn relation, the exchange-correlation kernel of TDDFT, and fluctuation dissipation theorem for DFT. This produces a natural method for generating new thermal exchange-correlation (XC) approximations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accurate atomic quantum defects from particle-particle random phase approximation

The accuracy of calculations of atomic Rydberg excitations cannot be judged by the usual measures, such as mean unsigned errors of many transitions. We show how to use quantum defect theory to (a) separate errors due to approximate ionization potentials, (b) extract smooth quantum defects to compare with experiment, and (c) quantify those defects with a few characteristic parameters. The particle-particle random phase approximation (pp-RPA) produces excellent Rydberg transitions that are an order of magnitude more accurate than those of time-dependent density functional theory with standard approximations. We even extract reasonably accurate defects from the lithium Rydberg series, despite the reference being open-shell. Our methodology can be applied to any Rydberg series of excitations with 4 transitions or more to extract the underlying threshold energy and characteristic quantum defect parameters. Our pp-RPA results set a demanding challenge for other excitation methods to match.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Corrections to Thomas-Fermi densities at turning points and beyond

Uniform semiclassical approximations for the number and kinetic-energy densities are derived for many non-interacting fermions in one-dimensional potentials with two turning points. The resulting simple, closed-form expressions contain the leading corrections to Thomas-Fermi theory, involve neither sums nor derivatives, are spatially uniform approximations, and are exceedingly accurate.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Density functional description of Coulomb blockade: Adiabatic or dynamic exchange-correlation?

Above the Kondo temperature, the Kohn-Sham zero-bias conductance of an Anderson junction has been shown to completely miss the Coulomb blockade. Within a standard model for the spectral function, we deduce a parameterization for both the onsite exchange-correlation potential and the bias drop as a function of the site occupation that applies for all correlation strengths. We use our results to sow doubt on the common interpretation of such corrections as arising from dynamical exchange-correlation contributions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Hubbard Dimer: A density functional case study of a many-body problem

This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of Density Functional Theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy-Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green's function, demonstrating the `failure' of the Kohn-Sham method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree-Fock and Bethe Ansatz Local Density Approximation) to show their successes and limitations. We also discuss and illustrate the concept of the derivative discontinuity. Useful appendices include analytic expressions for Density Functional energy components, several limits of the exact functional (weak- and strong-coupling, symmetric and asymmetric), the Kohn-Sham hopping energy functional for 3 sites, various adiabatic connection results, proofs of exact conditions for this model, and the origin of the Hubbard model from a minimal basis model for stretched H$_2$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Understanding Kernel Ridge Regression: Common behaviors from simple functions to density functionals

Accurate approximations to density functionals have recently been obtained via machine learning (ML). By applying ML to a simple function of one variable without any random sampling, we extract the qualitative dependence of errors on hyperparameters. We find universal features of the behavior in extreme limits, including both very small and very large length scales, and the noise-free limit. We show how such features arise in ML models of density functionals.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Almost Exact Exchange At Almost No Cost

A recently developed semiclassical approximation to exchange in one dimension is shown to be almost exact, with essentially no computational cost. The variational stability of this approximation is tested, and its far greater accuracy relative to local density functional calculations demonstrated. Even a fully orbital-free potential-functional calculation (no orbitals of any kind) yields little error relative to exact exchange, for more than one orbital.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Atomic correlation energies and the generalized gradient approximation

Careful extrapolation of atomic correlation energies suggests that $E_c$ tends to $-AZ\log{Z} + BZ$ as $Z$ tends to infinity, where $Z$ is the atomic number, $A$ is known, and $B$ is about 38 milliHartrees. The coefficients roughly agree with those of the high-density limit of the real-space construction of the generalized gradient approximation. An asymptotic coefficient, missed by previous derivations, is included in a revised approximation. The exchange is also corrected, reducing atomic errors considerably.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

DFT: A Theory Full of Holes?

This article is a rough, quirky overview of both the history and present state of the art of density functional theory. The field is so huge that no attempt to be comprehensive is made. We focus on the underlying exact theory, the origin of approximations, and the tension between empirical and non-empirical approaches. Many ideas are illustrated on the exchange energy and hole. Features unique to this article include how approximations can be systematically derived in a non-empirical fashion and a survey of warm dense matter.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exact and approximate Kohn-Sham potentials in ensemble density-functional theory

We construct exact Kohn-Sham potentials for the ensemble density-functional theory (EDFT) from the ground and excited states of helium. The exchange-correlation (XC) potential is compared with the quasi-local-density approximation and both single determinant and symmetry eigenstate ghost-corrected exact exchange approximations. Symmetry eigenstate Hartree-exchange recovers distinctive features of the exact XC potential and is used to calculate the correlation potential. Unlike the exact case, excitation energies calculated from these approximations depend on ensemble weight, and it is shown that only the symmetry eigenstate method produces an ensemble derivative discontinuity. Differences in asymptotic and near-ground-state behavior of exact and approximate XC potentials are discussed in the context of producing accurate optical gaps.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Excitations and benchmark ensemble density functional theory for two electrons

A new method for extracting ensemble Kohn-Sham potentials from accurate excited state densities is applied to a variety of two electron systems, exploring the behavior of exact ensemble density functional theory. The issue of separating the Hartree energy and the choice of degenerate eigenstates is explored. A new approximation, spin eigenstate Hartree-exchange (SEHX), is derived. Exact conditions that are proven include the signs of the correlation energy components, the virial theorem for both exchange and correlation, and the asymptotic behavior of the potential for small weights of the excited states. Many energy components are given as a function of the weights for two electrons in a one-dimensional flat box, in a box with a large barrier to create charge transfer excitations, in a three-dimensional harmonic well (Hooke's atom), and for the He atom singlet-triplet ensemble, singlet-triplet-singlet ensemble, and triplet bi-ensemble.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gedanken Densities and Exact Constraints in Density Functional Theory

Approximations to the exact density functional for the exchange-correlation energy of a many-electron ground state can be constructed by satisfying constraints that are universal, i.e., valid for all electron densities. Gedanken densities are designed for the purpose of this construction, but need not be realistic. The uniform electron gas is an old gedanken density. Here, we propose a spherical two-electron gedanken density in which the dimensionless density gradient can be an arbitrary positive constant wherever the density is non-zero. The Lieb-Oxford lower bound on the exchange energy can be satisfied within a generalized gradient approximation (GGA) by bounding its enhancement factor or simplest GGA exchange-energy density. This enhancement-factor bound is well known to be sufficient, but our gedanken density shows that it is also necessary. The conventional exact exchange-energy density satisfies no such local bound, but energy densities are not unique, and the simplest GGA exchange-energy density is not an approximation to it. We further derive a strongly and optimally tightened bound on the exchange enhancement factor of a two-electron density, which is satisfied by the local density approximation but is violated by all published GGA's or meta-GGA's. Finally, some consequences of the non-uniform density-scaling behavior for the asymptotics of the exchange enhancement factor of a GGA or meta-GGA are given.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ions in solution: Density Corrected Density Functional Theory (DC-DFT)

Standard density functional approximations often give questionable results for odd-electron radical complexes, with the error typically attributed to self-interaction. In density corrected density functional theory (DC-DFT), certain classes of density functional theory calculations are significantly improved by using densities more accurate than the self-consistent densities. We discuss how to identify such cases, and how DC-DFT applies more generally. To illustrate, we calculate potential energy surfaces of HO$\cdot$Cl$^-$ and HO$\cdot$H$_2$O complexes using various common approximate functionals, with and without this density correction. Commonly used approximations yield wrongly shaped surfaces and/or incorrect minima when calculated self consistently, while yielding almost identical shapes and minima when density corrected. This improvement is retained even in the presence of implicit solvent.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Kohn-Sham calculations with the exact functional

As a proof of principle, self-consistent Kohn--Sham calculations are performed with the exact exchange-correlation functional. Finding the exact functional for even one trial density requires solving the interacting Schrödinger equation many times. The density matrix renormalization group method makes this possible for one-dimensional, real-space systems of more than two interacting electrons. We illustrate and explore the convergence properties of the exact KS scheme for both weakly and strongly correlated systems. We also explore the spin-dependent generalization and densities for which the functional is ill defined.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum critical benchmark for density functional theory

Two electrons at the threshold of ionization represent a severe test case for electronic structure theory. A pseudospectral method yields a very accurate density of the two-electron ion with nuclear charge close to the critical value. Highly accurate energy components and potentials of Kohn-Sham density functional theory are given, as well as a useful parametrization of the critical density. The challenges for density functional approximations and the strength of correlation are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Testing and Using the Lewin-Lieb Bound in Density Functional Theory

Lewin and Lieb have recently proven several new bounds on the exchange-correlation energy that complement the Lieb-Oxford bound. We test these bounds for atoms, for slowly-varying gases, and for Hooke's atom, finding them usually less strict than the Lieb-Oxford bound. However, we also show that, if a generalized gradient approximation (GGA) is to guarantee satisfaction of the new bounds for all densities, new restrictions on the the exchange-correlation enhancement factor are implied.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Understanding Machine-learned Density Functionals

Kernel ridge regression is used to approximate the kinetic energy of non-interacting fermions in a one-dimensional box as a functional of their density. The properties of different kernels and methods of cross-validation are explored, and highly accurate energies are achieved. Accurate {\em constrained optimal densities} are found via a modified Euler-Lagrange constrained minimization of the total energy. A projected gradient descent algorithm is derived using local principal component analysis. Additionally, a sparse grid representation of the density can be used without degrading the performance of the methods. The implications for machine-learned density functional approximations are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Guaranteed convergence of the Kohn-Sham equations

A sufficiently damped iteration of the Kohn-Sham equations with the exact functional is proven to always converge to the true ground-state density, regardless of the initial density or the strength of electron correlation, for finite Coulomb systems. We numerically implement the exact functional for one-dimensional continuum systems and demonstrate convergence of the damped KS algorithm. More strongly correlated systems converge more slowly.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Non-existence of Taylor expansion in time due to cusps

In the usual treatment of electronic structure, all matter has cusps in the electronic density at nuclei. Cusps can produce non-analytic behavior in time, even in response to perturbations that are time-analytic. We analyze these non-analyticities in a simple case from many perspectives. We describe a method, the s-expansion, that can be used in several such cases, and illustrate it with a variety of examples. These include both the sudden appearance of electric fields and disappearance of nuclei, in both one and three dimensions. When successful, the s-expansion yields the dominant short-time behavior, no matter how strong the external electric field, but agrees with linear response theory in the weak limit. We discuss the relevance of these results to time-dependent density functional theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Potential functionals versus density functionals

Potential functional approximations are an intriguing alternative to density functional approximations. The potential functional that is dual to the Lieb density functional is defined and properties given. The relationship between Thomas-Fermi theory as a density functional and as a potential functional is derived. The properties of several recent semiclassical potential functionals are explored, especially in their approach to the large particle number and classical continuum limits. The lack of ambiguity in the energy density of potential functional approximations is demonstrated. The density-density response function of the semiclassical approximation is calculated and shown to violate a key symmetry condition.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermal Density Functional Theory in Context

This chapter introduces thermal density functional theory, starting from the ground-state theory and assuming a background in quantum mechanics and statistical mechanics. We review the foundations of density functional theory (DFT) by illustrating some of its key reformulations. The basics of DFT for thermal ensembles are explained in this context, as are tools useful for analysis and development of approximations. We close by discussing some key ideas relating thermal DFT and the ground state. This review emphasizes thermal DFT's strengths as a consistent and general framework.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Accuracy of density functionals for molecular electronics: the Anderson junction

The exact ground-state exchange-correlation functional of Kohn-Sham density functional theory yields the exact transmission through an Anderson junction at zero bias and temperature. The exact impurity charge susceptibility is used to construct the exact exchange-correlation potential. We analyze the successes and limitations of various types of approximations, including smooth and discontinuous functionals of the occupation, as well as symmetry-broken approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

One-dimensional Continuum Electronic Structure with the Density Matrix Renormalization Group and Its Implications For Density Functional Theory

We extend the density matrix renormalization group to compute exact ground states of continuum many-electron systems in one dimension with long-range interactions. We find the exact ground state of a chain of 100 strongly correlated artificial hydrogen atoms. The method can be used to simulate 1d cold atom systems and to study density functional theory in an exact setting. To illustrate, we find an interacting, extended system which is an insulator but whose Kohn-Sham system is metallic.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Perspective on density functional theory

Density functional theory (DFT) is an incredible success story. The low computational cost, combined with useful (but not yet chemical) accuracy, has made DFT a standard technique in most branches of chemistry and materials science. Electronic structure problems in a dazzling variety of fields are currently being tackled. However, DFT has many limitations in its present form: Too many approximations, failures for strongly correlated systems, too slow for liquids, etc. This perspective reviews some recent progress and ongoing challenges.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Reference electronic structure calculations in one dimension

Large strongly correlated systems provide a challenge to modern electronic structure methods, because standard density functionals usually fail and traditional quantum chemical approaches are too demanding. The density-matrix renormalization group method, an extremely powerful tool for solving such systems, has recently been extended to handle long-range interactions on real-space grids, but is most efficient in one dimension where it can provide essentially arbitrary accuracy. Such 1d systems therefore provide a theoretical laboratory for studying strong correlation and developing density functional approximations to handle strong correlation, {\em if} they mimic three-dimensional reality sufficiently closely. We demonstrate that this is the case, and provide reference data for exact and standard approximate methods, for future use in this area.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Avoiding Unbound Anions in Density Functional Calculations

Converged approximate density functional calculations usually do not bind anions, due to large self-interaction error. But Hartree-Fock calculations have no such prob- lem, producing negative HOMO energies. A recently proposed scheme for calculating DFT energies on HF densities is shown to work very well for molecules, better than the common practice of restricting the basis set, except for cases like CN, where the HF density is too inaccurate due to spin contamination.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Electronic structure via potential functional approximations

The universal functional of Hohenberg-Kohn is given as a coupling-constant integral over the density as a functional of the potential. Conditions are derived under which potential-functional approximations are variational. Construction via this method and imposition of these conditions are shown to greatly improve the accuracy of the non-interacting kinetic energy needed for orbital-free Kohn-Sham calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Finding Density Functionals with Machine Learning

Machine learning is used to approximate density functionals. For the model problem of the kinetic energy of non-interacting fermions in 1d, mean absolute errors below 1 kcal/mol on test densities similar to the training set are reached with fewer than 100 training densities. A predictor identifies if a test density is within the interpolation region. Via principal component analysis, a projected functional derivative finds highly accurate self-consistent densities. Challenges for application of our method to real electronic structure problems are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kondo effect given exactly by density functional theory

Transport through an Anderson junction (two macroscopic electrodes coupled to an Anderson impurity) is dominated by a Kondo peak in the spectral function at zero temperature. The exact single-particle Kohn-Sham potential of density functional theory reproduces the linear transport exactly, despite the lack of a Kondo peak in its spectral function. Using Bethe ansatz techniques, we calculate this potential exactly for all coupling strengths, including the cross-over from mean-field behavior to charge quantization caused by the derivative discontinuity. A simple and accurate interpolation formula is also given.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Role of Exact Conditions in TDDFT

This chapter is devoted to exact conditions in time-dependent density functional theory. Many conditions have been derived for the exact ground-state density functional, and several have played crucial roles in the construction of popular approximations. We believe that the reliability of the most fundamental approximation of any density functional theory, the local density approximation (LDA), is due to the exact conditions that it satisfies. Improved approximations should satisfy at least those conditions that LDA satisfies, plus others. (Which others is part of the art of functional approximation). In the time-dependent case, as we shall see, the adiabatic LDA (ALDA) plays the same role as LDA in the ground-state case, as it satisfies many exact conditions. But we do not have a generally applicable improvement beyond ALDA that includes nonlocality in time. For TDDFT, we have a surfeit of exact conditions, but that only makes finding those that are useful to impose an even more demanding task.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ionization potentials in the limit of large atomic number

By extrapolating the energies of non-relativistic atoms and their ions with up to 3000 electrons within Kohn-Sham density functional theory, we find that the ionization potential remains finite and increases across a row, even as $Z\rightarrow\infty$. The local density approximation becomes chemically accurate (and possibly exact) in some cases. Extended Thomas-Fermi theory matches the shell-average of both the ionization potential and density change. Exact results are given in the limit of weak electron-electron repulsion.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Leading corrections to local approximations

For the kinetic energy of 1d model finite systems the leading corrections to local approximations as a functional of the potential are derived using semiclassical methods. The corrections are simple, non-local functionals of the potential. Turning points produce quantum oscillations leading to energy corrections, which are completely different from the gradient corrections that occur in bulk systems with slowly-varying densities. Approximations that include quantum corrections are typically much more accurate than their local analogs. The consequences for density functional theory are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the density-potential mapping in time-dependent density functional theory

The key questions of uniqueness and existence in time-dependent density functional theory are usually formulated only for potentials and densities that are analytic in time. Simple examples, standard in quantum mechanics, lead however to non-analyticities. We reformulate these questions in terms of a non-linear Schrödinger equation with a potential that depends non-locally on the wavefunction.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Adiabatic Connection in the Low-Density Limit

In density functional theory (DFT), the exchange-correlation functional can be exactly expressed by the adiabatic connection integral. It has been noticed that as lambda goes to infinity, the lambda^(-1) term in the expansion of W(lambda) vanishes. We provide a simple but rigorous derivation to this exact condition in this work. We propose a simple parametric form for the integrand, satisfying this condition, and show that it is highly accurate for weakly-correlated two-electron systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Must Kohn-Sham Oscillator Strengths be Accurate at Threshold?

The exact ground-state Kohn-Sham(KS) potential for the helium atom is known from accurate wavefunction calculations of the ground-state density. The threshold for photoabsorption from this potential matches the physical system exactly. By carefully studying its absorption spectrum, we show the answer to the title question is no. To address this problem in detail, we generate a highly accurate simple fit of a two-electron spectrum near the threshold, and apply the method to both the experimental spectrum and that of the exact ground-state Kohn-Sham potential.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Non-empirical 'derivation' of B88 exchange functional

The B88 exchange energy density functional (created by Becke in 1988) is a crucial part of the most popular density functional in use today, B3LYP. B88 contains one empirical parameter which was fitted to Hartree-Fock exchange energies for the noble gas atoms. We show how local approximations to exchange become relatively exact under a very specific approach to the limit of large numbers, but that the usual gradient expansion does not. The leading corrections can be captured by generalized gradient approximations, producing a non-empirical derivation of the parameter in B88.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Charge Transfer in Partition Theory

The recently proposed Partition Theory (PT) [J.Phys.Chem.A 111, 2229 (2007)] is illustrated on a simple one-dimensional model of a heteronuclear diatomic molecule. It is shown that a sharp definition for the charge of molecular fragments emerges from PT, and that the ensuing population analysis can be used to study how charge redistributes during dissociation and the implications of that redistribution for the dipole moment. Interpreting small differences between the isolated parts' ionization potentials as due to environmental inhomogeneities, we gain insight into how electron localization takes place in H2+ as the molecule dissociates. Furthermore, by studying the preservation of the shapes of the parts as different parameters of the model are varied, we address the issue of transferability of the parts. We find good transferability within the chemically meaningful parameter regime, raising hopes that PT will prove useful in chemical applications.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Partition theory: A very simple illustration

We illustrate the main features of a recently proposed method based on ensemble density functional theory to divide rigorously a complex molecular system into its parts [M.H. Cohen and A. Wasserman, J. Phys. Chem. A 111, 2229 (2007)]. The illustrative system is an analog of the hydrogen molecule for which analytic expressions for the densities of the parts (hydrogen "atoms") are found along with the "reactivity potential" that enters the theory. While previous formulations of Chemical Reactivity Theory lead to zero, or undefined, values for the chemical hardness of the isolated parts, we demonstrate they can acquire a finite and positive hardness within the present formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A new challenge for time-dependent density-functional theory

Time-dependent density functional theory is thought to work well for the test cases of He and Be atoms. We perform a quantum defect analysis of the s, p, and d Rydberg states of Be with accurate ground state Kohn-Sham potentials. The s- and p-quantum defects are well described by the ALDA, but fails badly for the d-quantum defect. The same failure is observed in case of He. This provides a new challenge for functional development in time-dependent density functional theory.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Current-Density Functional Theory of the Response of Solids

The response of an extended periodic system to a homogeneous field (of wave-vector $q=0$) cannot be obtained from a $q=0$ time-dependent density functional theory (TDDFT) calculation, because the Runge-Gross theorem does not apply. Time-dependent {\em current}-density functional theory is needed and demonstrates that one key ingredient missing from TDDFT is the macroscopic current. In the low-frequency limit, in certain cases, density polarization functional theory is recovered and a formally exact expression for the polarization functional is given.

0 citations0 reviewsNo code linkedOpen access