Researcher profile

Kieron Burke

Kieron Burke contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
10topics
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

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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

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.

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.