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.
more »
« less
Exact and approximate energy sums in potential wells
Sums of the
- Award ID(s):
- 1856165
- PAR ID:
- 10303286
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Journal of Physics A: Mathematical and Theoretical
- Volume:
- 53
- Issue:
- 9
- ISSN:
- 1751-8113
- Page Range / eLocation ID:
- Article No. 095203
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We argue that the success of DFT can be understood in terms of a semiclassical expansion around a very specific limit. This limit was identified long ago by Lieb and Simon for the total electronic energy of a system. This is a universal limit of all electronic structure: atoms, molecules, and solids. For the total energy, Thomas-Fermi theory becomes relatively exact in the limit. The limit can also be studied for much simpler model systems, including non-interacting fermions in a one-dimensional well, where the WKB approximation applies for individual eigenvalues and eigenfunctions. Summation techniques lead to energies and densities that are functionals of the potential. We consider several examples in one dimension (fermions in a box, in a harmonic well, in a linear half-well, and in the Pöschl-Teller well. The effects of higher dimension are also illustrated with the three-dimensional harmonic well and the Bohr atom, non-interacting fermions in a Coulomb well. Modern density functional calculations use the Kohn-Sham scheme almost exclusively. The same semiclassical limit can be studied for the Kohn-Sham kinetic energy, for the exchange energy, and for the correlation energy. For all three, the local density approximation appears to become relatively exact in this limit. Recent work, both analytic and numerical, explores how this limit is approached, in an effort to deduce the leading corrections to the local approximation. A simple scheme, using the Euler-Maclaurin summation formula, is the result of many different attempts at this problem. In very simple cases, the correction formulas are much more accurate than standard density functionals. Several functionals are already in widespread use in both chemistry and materials that incorporate these limits, and prospects for the future are discussed.more » « less
-
Abstract Recent work has shown a deep connection between semilocal approximations in density functional theory and the asymptotics of the sum of the Wentzel–Kramers–Brillouin (WKB) semiclassical expansion for the eigenvalues. However, all examples studied to date have potentials with only real classical turning points. But systems with complex turning points generate subdominant (SD) terms beyond those in the WKB series. The simplest case is a pure quartic oscillator. We show how to generalize the asymptotics of eigenvalue sums to include SD contributions to the sums, if they are known for the eigenvalues. These corrections to WKB greatly improve accuracy for eigenvalue sums, especially for many levels. We obtain further improvements to the sums through hyperasymptotics. For the lowest level, our summation method has error below 2 × 10 −4 . For the sum of the lowest ten levels, our error is less than 10 −22 . We report all results to many digits and include copious details of the asymptotic expansions and their derivation.more » « less
-
more » « less
A<sc>bstract</sc> Spacetime wormholes can provide non-perturbative contributions to the gravitational path integral that make the actual number of states
e S in a gravitational system much smaller than the number of states$$ {e}^{S_{\textrm{p}}} $$ N =O (e S ) of the possible$$ {e}^{S_{\textrm{p}}} $$ N ≪e S (in which case no quantum extremal surfaces can arise). In particular, even for smallN , the insertion of generic operators into the path integral can cause the non-perturbative time evolution to differ dramatically from perturbative expectations. On the other hand, this discrepancy is small for the special case where the inserted operators are non-trivial only in a subspace of dimensionD ≪e S . We thus study this latter case in detail. We also discuss potential implications for more realistic gravitational systems. -
more » « less
Abstract Consider a quantum cat map
M associated with a matrix$$A\in {{\,\textrm{Sp}\,}}(2n,{\mathbb {Z}})$$ M on any nonempty open set in the position–frequency space satisfies a lower bound which is uniform in the semiclassical limit, under two assumptions: (1) there is a unique simple eigenvalue ofA of largest absolute value and (2) the characteristic polynomial ofA is irreducible over the rationals. This is similar to previous work (Dyatlov and Jin in Acta Math 220(2):297–339, 2018; Dyatlov et al. in J Am Math Soc 35(2):361–465, 2022) on negatively curved surfaces and (Schwartz in The full delocalization of eigenstates for the quantized cat map, 2021) on quantum cat maps with$$n=1$$
