Controlled quantum treatment of the outer region, interior-exterior matching, and interior accuracy of the exchange-corrected Thomas-Fermi (Dirac) atom

Dirac's 1930 derivation of the exchange correction to the Thomas-Fermi (TF) atom reduced the antisymmetrized many-electron problem by a semiclassical phase-space mapping that treats coordinates and momenta as commuting classical variables and drops the spin variables. The resulting Thomas-Fermi-Dirac (TFD) equation - the origin of the local-density approximation (LDA) for exchange still used throughout density-functional theory - gives nonphysical negative, oscillating solutions in the atom's exterior, where the occupied phase-space volumes become comparable to the elementary quantum cell $(2\pi\hbar)^3$ and the discrete level structure matters. The open problem is threefold and unified by one objective: (1) a quantum-mechanically controlled description of the outer region where the semiclassical picture fails; (2) a controlled (error-bounded or rigorously asymptotic) matching between the interior semiclassical exchange-corrected solution and the exterior quantum regime, replacing the phenomenological cutoffs in use since the 1950s; (3) a quantification of the interior accuracy of the TFD equation that Dirac expected but did not quantify. Later literature resolved the leading-order questions (TF exactness; the Dirac term as the exact leading exchange contribution; the energy expansion through order $Z^{5/3}$) and produced purely practical exterior repairs, but the controlled matching and the beyond-LDA exchange asymptotics (leading $Z\ln Z$ correction known only numerically, 2022-2024) remain open.

Background

In the Thomas-Fermi (TF) statistical model of a heavy atom, the electrons are treated as a zero-temperature Fermi gas in the self-consistent electrostatic potential of the nucleus (nuclear charge $Z$) and of the electron cloud: the lowest-energy region of one-particle phase space is assumed saturated, with two electrons of opposite spin in each elementary cell of volume $(2\pi\hbar)^3$ (the Pauli exclusion principle), and the rest empty. Because an exact solution of the many-electron Schrödinger equation is out of reach for large $Z$, this semiclassical model and its refinements form the standard asymptotic limit of atomic structure.

Dirac (1930) supplied the exchange correction to this model. Exchange is a purely quantum effect of antisymmetry: electrons are fermions, so the many-electron wavefunction changes sign under interchange of any two of them; the Hartree-Fock approximation implements this exactly through a determinant of one-electron orbitals, and the resulting exchange term lowers the electrostatic self-energy of same-spin electrons. Dirac evaluated the exchange kernel in a semiclassical limit - treating position and momentum as commuting classical variables and neglecting the spin variables in the phase-space mapping (the spin degeneracy factor 2 survives only in the saturation count) - and found that exchange contributes a local energy density proportional to $-\rho^{4/3}$, where $\rho$ is the electron number density. This is the original local-density approximation (LDA) for exchange, the term that still enters essentially every density-functional exchange functional. Inserted into the TF self-consistency condition it yields the Thomas-Fermi-Dirac (TFD) equation, which in Dirac's words differs from the corresponding Thomas equation 'only on account of its having a term linear in $P$', where $P(r)$ is the screening function through which the self-consistent potential is expressed.

Dirac himself stated the limits of the construction. For sufficiently large $r$ the extra term drives $P$ negative, after which it oscillates with decreasing amplitude - 'a negative value for P has, of course, no physical meaning' - and he attributed the failure to the breakdown, in the outer region of the atom, of the approximation 'of regarding p as a function of commuting q's and p's': far from the nucleus the density is so low that the occupied phase-space volumes are comparable to the elementary cell $(2\pi\hbar)^3$, the quantum levels are discrete and sparsely spaced, and the continuous classical phase-space picture is untenable. He nevertheless expected the corrected equation to be more accurate than Thomas' equation in the interior of the atom, without quantifying this.

The subsequent literature establishes what is and is not settled. Rigorously, TF theory is the exact leading description of large atoms for the ground-state energy (Lieb and Simon 1977); Scott (1952) identified the leading quantum (near-nuclear) correction, and Schwinger (1980, 1981) computed the leading quantum and exchange corrections heuristically; Fefferman and Seco (1994) proved that Dirac's term is the exact leading contribution to the exchange energy of a large atom and carried the rigorous energy expansion through order $Z^{5/3}$, including the Dirac and Schwinger corrections. The spin question at leading order is likewise settled: the spin-polarized local exchange is standard (von Barth and Hedin 1972), and the $-\int\rho^{4/3}$ form is the shape of a rigorous lower bound on the indirect Coulomb energy (Lieb and Oxford 1981). In the exterior, by contrast, all existing treatments are repairs rather than controlled theories: the Latter truncation (1955) and the Gombás-type cutoff (Csavinszky 1969) patch the potential or density; the Englert-Schwinger program (1982, 1984, 1985) constructs a matched model (TF interior, quantum treatment of strongly bound electrons, edge treatment with a shell-averaged cut-off) whose cut-off energy is chosen semi-empirically; Patil (1984) builds an asymptotically correct density tail $\rho(r)\sim r^{\alpha}e^{-br}$ with $b=\sqrt{2I}$, $I$ the experimental ionization energy. The unphysical oscillations of TFD-type densities were analyzed by Goodisman (1973), and in the TFDW model (TFD plus a von Weizsäcker inhomogeneity, i.e. gradient, correction to the kinetic energy) the Dirac term destabilizes the exterior so strongly that minimizers fail to exist in some regimes (Lu and Otto 2014; Frank, Nam and Van Den Bosch 2018). What no later work supplies is a controlled answer to the limitations Dirac acknowledged: a quantum treatment of the outer region where his semiclassical approximation breaks down, a controlled matching to the interior, and a quantification of the interior accuracy he expected but did not establish - a question that modern work approaches through the beyond-LDA exchange asymptotics, whose leading $Z\ln Z$ correction is known only numerically (Argaman et al 2022; Redd et al 2024).

Problem Statement

The problem posed by the acknowledged limitations of Dirac's 1930 semiclassical reduction of exchange in the Thomas atom is to replace the present patched structure by a controlled quantum-mechanical theory of the exchange-corrected statistical atom. Its components, which share this single objective, are:

(1) Outer region. Give a quantum-mechanical treatment of the region sufficiently far from the nucleus that the number of occupied single-particle quantum states per phase-space cell of volume $(2\pi\hbar)^3$ is no longer large - equivalently, where the phase-space volume occupied per electron is comparable to the elementary cell, so that treating coordinates and momenta as commuting classical variables and ignoring the discreteness of the quantum levels ceases to be valid - that reproduces the correct physical behavior there: a nonnegative density with an exponential falloff governed by the ionization energy, rather than the negative, oscillating solutions of the TFD screening equation. Such a treatment must be controlled (explicit error bounds or a rigorous large-$Z$ asymptotic statement) for all sufficiently large $Z$, and reproducing the known ionization-energy tail behaviour is a necessary check on it, not by itself a resolution of this component.

(2) Matching. Construct, with controlled errors (explicit error bounds or a rigorous large-$Z$ asymptotic statement), the crossover between the interior semiclassical exchange-corrected solution and the exterior quantum regime, in place of the phenomenological cutoffs and matchings in current use (Latter truncation, Gombás-type density cutoff, Englert-Schwinger-type shell-averaged boundary matching, ionization-energy tails).

(3) Interior accuracy. Quantify the accuracy of the exchange-corrected Thomas equation in the interior of heavy atoms, for which Dirac could state only the unquantified expectation that it improves on Thomas' equation. The leading-order part of this question is settled (TF exactness; the Dirac term as the exact leading exchange contribution); what survives is the accuracy beyond the leading term, which the modern literature approaches through the large-$Z$ expansion of the exchange energy beyond the leading Dirac local-density term and whose first correction is known at present only numerically or conjecturally. A resolution of this component requires a controlled determination (explicit error bounds or a rigorous large-$Z$ asymptotic statement) of that correction; restating its known numerical or conjectural value does not by itself resolve the component.

References

1. P. A. M. Dirac, 'Note on Exchange Phenomena in the Thomas Atom', Mathematical Proceedings of the Cambridge Philosophical Society 26(3), 376-385 (1930). DOI 10.1017/S0305004100016108, https://doi.org/10.1017/S0305004100016108
2. J. M. C. Scott, 'The binding energy of the Thomas-Fermi Atom', Philosophical Magazine 43, 859-867 (1952). DOI 10.1080/14786440808520234, https://doi.org/10.1080/14786440808520234
3. Richard Latter, 'Atomic Energy Levels for the Thomas-Fermi and Thomas-Fermi-Dirac Potential', Physical Review 99, 510-519 (1955). DOI 10.1103/PhysRev.99.510, https://doi.org/10.1103/PhysRev.99.510
4. P. Csavinszky, 'Approximate Analytical Solutions of the Thomas–Fermi–Dirac and Thomas–Fermi–Dirac–Gombás Equations', The Journal of Chemical Physics 50, 2476-2485 (1969). DOI 10.1063/1.1671405, https://doi.org/10.1063/1.1671405
5. U. von Barth and L. Hedin, 'A local exchange-correlation potential for the spin polarized case. I', Journal of Physics C: Solid State Physics 5, 1629-1642 (1972). DOI 10.1088/0022-3719/5/13/012, https://doi.org/10.1088/0022-3719/5/13/012
6. Jerry Goodisman, 'Oscillations in the Thomas-Fermi-Dirac electron density', Theoretica Chimica Acta 31, 101-102 (1973). DOI 10.1007/BF00527444, https://doi.org/10.1007/BF00527444
7. Elliott H. Lieb and Barry Simon, 'The Thomas-Fermi theory of atoms, molecules and solids', Advances in Mathematics 23, 22-116 (1977). DOI 10.1016/0001-8708(77)90108-6, https://doi.org/10.1016/0001-8708(77)90108-6
8. Julian Schwinger, 'Thomas-Fermi model: The leading correction', Physical Review A 22, 1827-1832 (1980). DOI 10.1103/PhysRevA.22.1827, https://doi.org/10.1103/PhysRevA.22.1827
9. Julian Schwinger, 'Thomas-Fermi model: The second correction', Physical Review A 24, 2353-2361 (1981). DOI 10.1103/PhysRevA.24.2353, https://doi.org/10.1103/PhysRevA.24.2353
10. Elliott H. Lieb and Stephen Oxford, 'Improved lower bound on the indirect Coulomb energy', International Journal of Quantum Chemistry 19, 427-439 (1981). DOI 10.1002/qua.560190306, https://doi.org/10.1002/qua.560190306
11. Berthold-Georg Englert and Julian Schwinger, 'Thomas-Fermi revisited: The outer regions of the atom', Physical Review A 26, 2322-2329 (1982). DOI 10.1103/PhysRevA.26.2322, https://doi.org/10.1103/PhysRevA.26.2322
12. Joseph G. Conlon, 'Semi-classical limit theorems for Hartree-Fock theory', Communications in Mathematical Physics 88, 133-150 (1983). DOI 10.1007/BF01206884, https://doi.org/10.1007/BF01206884
13. Berthold-Georg Englert and Julian Schwinger, 'Statistical atom: Handling the strongly bound electrons', Physical Review A 29, 2331-2338 (1984). DOI 10.1103/PhysRevA.29.2331, https://doi.org/10.1103/PhysRevA.29.2331
14. Berthold-Georg Englert and Julian Schwinger, 'Statistical atom: Some quantum improvements', Physical Review A 29, 2339-2352 (1984). DOI 10.1103/PhysRevA.29.2339, https://doi.org/10.1103/PhysRevA.29.2339
15. S. H. Patil, 'Asymptotic region of a Thomas-Fermi atom', The Journal of Chemical Physics 80, 5073-5075 (1984). DOI 10.1063/1.446576, https://doi.org/10.1063/1.446576
16. Berthold-Georg Englert and Julian Schwinger, 'Semiclassical atom', Physical Review A 32, 26-35 (1985). DOI 10.1103/PhysRevA.32.26, https://doi.org/10.1103/PhysRevA.32.26
17. Berthold-Georg Englert, 'Semiclassical Theory of Atoms', Lecture Notes in Physics vol. 300, Springer, Berlin (1988). DOI 10.1007/3-540-19204-2, https://doi.org/10.1007/3-540-19204-2
18. C. Fefferman and L. A. Seco, 'On the Dirac and Schwinger Corrections to the Ground-State Energy of an Atom', Advances in Mathematics 107, 1-185 (1994). DOI 10.1006/aima.1994.1060, https://doi.org/10.1006/aima.1994.1060
19. Jianfeng Lu and Felix Otto, 'Nonexistence of a Minimizer for Thomas–Fermi–Dirac–von Weizsäcker Model', Communications on Pure and Applied Mathematics 67, 1605-1617 (2014). DOI 10.1002/cpa.21477, https://doi.org/10.1002/cpa.21477
20. Jouko Lehtomäki and Olga Lopez-Acevedo, 'Self-consistent assessment of Englert-Schwinger model on atomic properties', The Journal of Chemical Physics 147, 234102 (2017). DOI 10.1063/1.5000908, https://doi.org/10.1063/1.5000908; preprint arXiv:1711.05639, https://arxiv.org/abs/1711.05639
21. Victor Ivrii, 'Asymptotics of the ground state energy in the relativistic settings', arXiv:1707.07014 (2017). https://arxiv.org/abs/1707.07014
22. Rupert L. Frank, Phan Thành Nam and Hanne Van Den Bosch, 'The Ionization Conjecture in Thomas–Fermi–Dirac–von Weizsäcker Theory', Communications on Pure and Applied Mathematics 71, 577-614 (2018). DOI 10.1002/cpa.21717, https://doi.org/10.1002/cpa.21717; preprint arXiv:1606.07355, https://arxiv.org/abs/1606.07355
23. Nathan Argaman, Jeremy Redd, Antonio C. Cancio and Kieron Burke, 'Leading Correction to the Local Density Approximation for Exchange in Large-Z Atoms', Physical Review Letters 129, 153001 (2022). DOI 10.1103/PhysRevLett.129.153001, https://doi.org/10.1103/PhysRevLett.129.153001; preprint arXiv:2204.01030, https://arxiv.org/abs/2204.01030
24. Jeremy J. Redd, Antonio C. Cancio, Nathan Argaman and Kieron Burke, 'Investigations of the exchange energy of neutral atoms in the large-Z limit', The Journal of Chemical Physics 160, 044101 (2024). DOI 10.1063/5.0179278, https://doi.org/10.1063/5.0179278; preprint arXiv:2310.03982, https://arxiv.org/abs/2310.03982