Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1012.0347 · Physical Review B

Gradient corrections to the kinetic energy density functional of a two-dimensional Fermi gas at finite temperature

Brandon P. van Zyl, K. Berkane, K. Bencheikh, Aaron Farrell

Published May 31, 2011• 6 citations
arXivHTMLPDF

Abstract

We examine the leading-order semiclassical gradient corrections to the noninteracting kinetic-energy density functional of a two-dimensional Fermi gas by applying the extended Thomas-Fermi theory at finite temperature. We find a nonzero von Weizs\"acker-like gradient correction, which in the high-temperature limit goes over to the functional form $(\ensuremath{\hbar}{}^{2}/24m)(\ensuremath{\nabla}\ensuremath{\rho}){}^{2}/\ensuremath{\rho}$. Our work provides a theoretical justification for the inclusion of gradient corrections in applications of density-functional theory to inhomogeneous two-dimensional Fermi systems at any finite temperature.

DOI
10.1103/physrevb.83.195136

Related papers

Nonlocal kinetic energy functional for an inhomogeneous two-dimensional Fermi gas

Brandon P. van Zyl, Aaron Farrell, E. Zaremba +3 more

The average-density approximation is used to construct a nonlocal kinetic energy functional for an inhomogeneous two-dimensional Fermi gas. This functional is then used to formulate a Thomas-Fermi-von Weizs\"acker-like theory for the description of the ground-state properties of the system. The quality of the kinetic energy functional is tested by performing a fully self-consistent calculation for an ideal, harmonically confined, two-dimensional system. Good agreement with exact results are found, with the number and kinetic energy densities exhibiting oscillatory structure associated with the nonlocality of the energy functional. Most importantly, this functional shows a marked improvement over the two-dimensional Thomas-Fermi von Weizs\"acker theory, particularly in the vicinity of the classically forbidden region.

arXiv:1311.5608•Physical Review A•Feb 3, 2014•12 citations

Kirzhnits gradient expansion for aD-dimensional Fermi gas

Luca Salasnich

For an ideal D-dimensional Fermi gas under generic external confinement we derive the correcting coefficient $(D-2)/3D$ of the von Weizsacker term in the kinetic energy density. To obtain this coefficient we use the Kirzhnits semiclassical expansion of the number operator up to the second order in the Planck constant $\hbar$. Within this simple and direct approach we determine the differential equation of the density profile and the density functional of the Fermi gas. In the case D=2 we find that the Kirzhnits gradient corrections vanish to all order in $\hbar$.

arXiv:0707.0833•Journal of Physics A Mathematical and Theoretical•Aug 1, 2007•46 citations

Extended Thomas-Fermi Model at Finite Temperature

M. Brack

We derive for the first time second-order gradient corrections to the Thomas-Fermi kinetic-energy and entropy density functionals for a fermion system at finite temperature. Density variational calculations for a heated nucleus with a Skyrme interaction lead to excellent agreement with results of Hartree-Fock calculations.

Physical Review Letters•Jul 9, 1984•47 citations

Testing the nonlocal kinetic energy functional of an inhomogeneous, two-dimensional degenerate Fermi gas within the average density approximation

Joseph Towers, Brandon P. van Zyl, W. Kirkby

In a recent paper [B. P. van Zyl et al., Phys. Rev. A 89, 022503 (2014)], the average density approximation (ADA) was implemented to develop a parameter-free, nonlocal kinetic energy functional to be used in the orbital-free density functional theory of an inhomogeneous, two-dimensional (2D) Fermi gas. In this work, we provide a detailed comparison of self-consistent calculations within the ADA with the exact results of the Kohn-Sham density functional theory and the elementary Thomas-Fermi (TF) approximation. We demonstrate that the ADA for the 2D kinetic energy functional works very well under a wide variety of confinement potentials, even for relatively small particle numbers. Remarkably, the TF approximation for the kinetic energy functional, without any gradient corrections, also yields good agreement with the exact kinetic energy for all confining potentials considered, although at the expense of the spatial and kinetic energy densities exhibiting poor pointwise agreement, particularly near the TF radius. Our findings illustrate that the ADA kinetic energy functional yields accurate results for both the local and global equilibrium properties of an inhomogeneous 2D Fermi gas, without the need for any fitting parameters.

arXiv:1507.03118•Physical Review B•Aug 20, 2015•1 citation

Manifestly Hermitian semiclassical expansion for the one-particle density matrix of a two-dimensional Fermi gas

K. Bencheikh, Brandon P. van Zyl, K. Berkane

The semiclassical $\ensuremath{\hbar}$ expansion of the one-particle density matrix for a two-dimensional Fermi gas is calculated within the Wigner transform method of B. Grammaticos and A. Voros [Ann. Phys. (N.Y.) 123, 359 (1979)], originally developed in the context of nuclear physics. The method of Grammaticos and Voros has the virtue of preserving both the Hermiticity and idempotency of the density matrix to all orders in the $\ensuremath{\hbar}$ expansion. As a topical application, we use our semiclassical expansion to go beyond the local-density approximation for the construction of the total dipole-dipole interaction energy functional of a two-dimensional, spin-polarized dipolar Fermi gas. We find a finite, second-order gradient correction to the Hartree-Fock energy, which takes the form $\ensuremath{\varepsilon}{(\mathbf{\ensuremath{\nabla}}\ensuremath{\rho})}^{2}/\sqrt{\ensuremath{\rho}}$, with $\ensuremath{\varepsilon}$ being small ($|\ensuremath{\varepsilon}|\ensuremath{\ll}1$) and negative. We test the quality of the corrected energy by comparing it with the exact results available for harmonic confinement. Even for small particle numbers, the gradient correction to the dipole-dipole energy provides a significant improvement over the local-density approximation.

arXiv:1605.08014•Physical review. B./Physical review. B•Aug 15, 2016•3 citations

Some tests of an approximate density functional for the ground-state kinetic energy of a fermion system

Conyers Herring, Manoj Chopra

Any proposed approximation to the ground-state kinetic energy of a system of noninteracting fermions in terms of the particle density \ensuremath{\rho}(r) can be used to estimate the potential field v(r) that will give rise to a given \ensuremath{\rho}(r), or the \ensuremath{\rho}(r) that will result from placing a given number of particles in a given v(r). Comparison with exact quantum-mechanical calculations thus gives two possible types of tests for any proposed kinetic-energy functional. This paper reports such tests for a recently proposed nonlocal functional applicable to one-dimensional systems, comparing its predictions for several simple problems not only with the correct answers but also with the predictions of Thomas-Fermi and Thomas-Fermi-Weizs\"acker functionals, to which the new functional proves considerably superior. The comparisons yield useful insights on the virtues and defects of the new functional, and on the directions in which improvements should be sought. Particular attention is devoted to Friedel oscillations and shell structure. It is shown that functionals of the Thomas-Fermi-Weizs\"acker type can never predict multiple maxima in \ensuremath{\rho} if v has only a single minimum and no maxima; the new functional does not have this defect, though it yields density oscillations rather weaker than the exact quantum-mechanical calculations. A natural general inference from the present tests is that any satisfactory kinetic-energy functional must (as ours does) replace the entire Thomas-Fermi term by a nonlocal expression.

Physical Review A•Jan 1, 1988•54 citations
We examine the leading-order semiclassical gradient corrections to the noninteracting kinetic-energy density functional of a two-dimensional Fermi gas by applying the extended Thomas-Fermi theory at finite temperature. We find a nonzero von Weizs"acker-like gradient correction, which in the high-temperature limit goes over to the functional form (\ensuremathℏ2/24m)(\ensuremath∇\ensuremathρ)2/\ensuremathρ(\ensuremath{\hbar}{}^{2}/24m)(\ensuremath{\nabla}\ensuremath{\rho}){}^{2}/\ensuremath{\rho}(\ensuremathℏ2/24m)(\ensuremath∇\ensuremathρ)2/\ensuremathρ. Our work provides a theoretical justification for the inclusion of gradient corrections in applications of density-functional theory to inhomogeneous two-dimensional Fermi systems at any finite temperature.
The average-density approximation is used to construct a nonlocal kinetic energy functional for an inhomogeneous two-dimensional Fermi gas. This functional is then used to formulate a Thomas-Fermi-von Weizs\"acker-like theory for the description of the ground-state properties of the system. The quality of the kinetic energy functional is tested by performing a fully self-consistent calculation for an ideal, harmonically confined, two-dimensional system. Good agreement with exact results are found, with the number and kinetic energy densities exhibiting oscillatory structure associated with the nonlocality of the energy functional. Most importantly, this functional shows a marked improvement over the two-dimensional Thomas-Fermi von Weizs\"acker theory, particularly in the vicinity of the classically forbidden region.
For an ideal D-dimensional Fermi gas under generic external confinement we derive the correcting coefficient (D−2)/3D(D-2)/3D(D−2)/3D of the von Weizsacker term in the kinetic energy density. To obtain this coefficient we use the Kirzhnits semiclassical expansion of the number operator up to the second order in the Planck constant ℏ\hbarℏ. Within this simple and direct approach we determine the differential equation of the density profile and the density functional of the Fermi gas. In the case D=2 we find that the Kirzhnits gradient corrections vanish to all order in ℏ\hbarℏ.
In a recent paper [B. P. van Zyl et al., Phys. Rev. A 89, 022503 (2014)], the average density approximation (ADA) was implemented to develop a parameter-free, nonlocal kinetic energy functional to be used in the orbital-free density functional theory of an inhomogeneous, two-dimensional (2D) Fermi gas. In this work, we provide a detailed comparison of self-consistent calculations within the ADA with the exact results of the Kohn-Sham density functional theory and the elementary Thomas-Fermi (TF) approximation. We demonstrate that the ADA for the 2D kinetic energy functional works very well under a wide variety of confinement potentials, even for relatively small particle numbers. Remarkably, the TF approximation for the kinetic energy functional, without any gradient corrections, also yields good agreement with the exact kinetic energy for all confining potentials considered, although at the expense of the spatial and kinetic energy densities exhibiting poor pointwise agreement, particularly near the TF radius. Our findings illustrate that the ADA kinetic energy functional yields accurate results for both the local and global equilibrium properties of an inhomogeneous 2D Fermi gas, without the need for any fitting parameters.
The semiclassical $\ensuremath{\hbar}$ expansion of the one-particle density matrix for a two-dimensional Fermi gas is calculated within the Wigner transform method of B. Grammaticos and A. Voros [Ann. Phys. (N.Y.) 123, 359 (1979)], originally developed in the context of nuclear physics. The method of Grammaticos and Voros has the virtue of preserving both the Hermiticity and idempotency of the density matrix to all orders in the $\ensuremath{\hbar}$ expansion. As a topical application, we use our semiclassical expansion to go beyond the local-density approximation for the construction of the total dipole-dipole interaction energy functional of a two-dimensional, spin-polarized dipolar Fermi gas. We find a finite, second-order gradient correction to the Hartree-Fock energy, which takes the form $\ensuremath{\varepsilon}{(\mathbf{\ensuremath{\nabla}}\ensuremath{\rho})}^{2}/\sqrt{\ensuremath{\rho}}$, with $\ensuremath{\varepsilon}$ being small ($|\ensuremath{\varepsilon}|\ensuremath{\ll}1$) and negative. We test the quality of the corrected energy by comparing it with the exact results available for harmonic confinement. Even for small particle numbers, the gradient correction to the dipole-dipole energy provides a significant improvement over the local-density approximation.
Any proposed approximation to the ground-state kinetic energy of a system of noninteracting fermions in terms of the particle density \ensuremath{\rho}(r) can be used to estimate the potential field v(r) that will give rise to a given \ensuremath{\rho}(r), or the \ensuremath{\rho}(r) that will result from placing a given number of particles in a given v(r). Comparison with exact quantum-mechanical calculations thus gives two possible types of tests for any proposed kinetic-energy functional. This paper reports such tests for a recently proposed nonlocal functional applicable to one-dimensional systems, comparing its predictions for several simple problems not only with the correct answers but also with the predictions of Thomas-Fermi and Thomas-Fermi-Weizs\"acker functionals, to which the new functional proves considerably superior. The comparisons yield useful insights on the virtues and defects of the new functional, and on the directions in which improvements should be sought. Particular attention is devoted to Friedel oscillations and shell structure. It is shown that functionals of the Thomas-Fermi-Weizs\"acker type can never predict multiple maxima in \ensuremath{\rho} if v has only a single minimum and no maxima; the new functional does not have this defect, though it yields density oscillations rather weaker than the exact quantum-mechanical calculations. A natural general inference from the present tests is that any satisfactory kinetic-energy functional must (as ours does) replace the entire Thomas-Fermi term by a nonlocal expression.