TY - JOUR
T1 - Upper bound to the gradient-based kinetic energy density of noninteracting electrons in an external potential
AU - Witt, William C.
AU - Jiang, Kaili
AU - Carter, Emily Ann
N1 - Publisher Copyright:
© 2019 Author(s).
PY - 2019/8/14
Y1 - 2019/8/14
N2 - We examine a simple upper bound to the gradient-based kinetic energy density (KED) of noninteracting electrons in an external potential: t(r) ≤ [μ-v(r)]n(r)+14λ2n(r), where t(r) is the gradient-based (non-negative) KED, μ is the Fermi energy, v(r) is the external potential, and n(r) is the electron density. The bound emerges naturally from a well-known expression for t(r), leading to an intuitive physical interpretation. For example, t(r) approaches the upper bound in regions where the electron density consists mainly of contributions from states with energies close to the Fermi energy. This upper bound complements the orbital-free lower bound provided by the gradient form of the von Weizsäcker (vW) KED, which is also non-negative. Both bounds yield t(r) exactly for single-orbital systems, and accordingly, they merge in single-orbital regions of more general systems. We demonstrate the universality of the two bounds over a wide range of test systems, including model potentials, atoms, diatomic molecules, and pseudopotential approximations of crystals. We also show that the exact t(r) frequently exceeds the sum of the vW and Thomas-Fermi KEDs, rendering that sum unsuitable as a strict upper bound to the gradient-based KED.
AB - We examine a simple upper bound to the gradient-based kinetic energy density (KED) of noninteracting electrons in an external potential: t(r) ≤ [μ-v(r)]n(r)+14λ2n(r), where t(r) is the gradient-based (non-negative) KED, μ is the Fermi energy, v(r) is the external potential, and n(r) is the electron density. The bound emerges naturally from a well-known expression for t(r), leading to an intuitive physical interpretation. For example, t(r) approaches the upper bound in regions where the electron density consists mainly of contributions from states with energies close to the Fermi energy. This upper bound complements the orbital-free lower bound provided by the gradient form of the von Weizsäcker (vW) KED, which is also non-negative. Both bounds yield t(r) exactly for single-orbital systems, and accordingly, they merge in single-orbital regions of more general systems. We demonstrate the universality of the two bounds over a wide range of test systems, including model potentials, atoms, diatomic molecules, and pseudopotential approximations of crystals. We also show that the exact t(r) frequently exceeds the sum of the vW and Thomas-Fermi KEDs, rendering that sum unsuitable as a strict upper bound to the gradient-based KED.
UR - http://www.scopus.com/inward/record.url?scp=85070663591&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85070663591&partnerID=8YFLogxK
U2 - 10.1063/1.5108896
DO - 10.1063/1.5108896
M3 - Article
AN - SCOPUS:85070663591
SN - 0021-9606
VL - 151
JO - Journal of Chemical Physics
JF - Journal of Chemical Physics
IS - 6
M1 - 064113
ER -