We present closed analytical expressions for the particle and kinetic-energy spatial densities at finite temperatures for a system of noninteracting fermions (bosons) trapped in a d-dimensional harmonic-oscillator potential. For $d=2$ and $3,$ exact expressions for the N-particle densities are used to calculate perturbatively the temperature dependence of the splittings of the energy levels in a given shell due to a very weak interparticle interaction in a dilute Fermi gas. In two dimensions, we obtain analytically the surprising result that the l degeneracy in a harmonic-oscillator shell is not lifted in the lowest order even when the exact, rather than the Thomas-Fermi expression for the particle density is used. We also demonstrate rigorously (in two dimensions) the reduction of the exact zero-temperature fermionic expressions to the Thomas-Fermi form in the large-N limit.We present closed analytical expressions for the particle and kinetic-energy spatial densities at finite temperatures for a system of noninteracting fermions (bosons) trapped in a d-dimensional harmonic-oscillator potential. For $d=2$ and $3,$ exact expressions for the N-particle densities are used to calculate perturbatively the temperature dependence of the splittings of the energy levels in a given shell due to a very weak interparticle interaction in a dilute Fermi gas. In two dimensions, we obtain analytically the surprising result that the l degeneracy in a harmonic-oscillator shell is not lifted in the lowest order even when the exact, rather than the Thomas-Fermi expression for the particle density is used. We also demonstrate rigorously (in two dimensions) the reduction of the exact zero-temperature fermionic expressions to the Thomas-Fermi form in the large-N limit.

arXiv:cond-mat/0209460Physical Review AFeb 24, 200336 citations

Asymptotically exact solutions for Fermi systems

I.V. Simenog

Theoretical and Mathematical PhysicsAug 1, 19743 citations

Variational Thomas-Fermi theory of a nonuniform Bose condensate at zero temperature

Eddy Timmermans, Paolo Tommasini, Kerson Huang

We derive a description of the spatially inhomogeneous Bose-Einstein condensate, treating the system locally as a homogeneous Bose gas. This approach, similar to the Thomas-Fermi model for the inhomogeneous many-particle fermion system, is well suited to describe the atomic Bose-Einstein condensates that have recently been obtained experimentally through atomic trapping and cooling. In this paper, we confine our attention to the zero-temperature case, although the treatment can be generalized to finite temperatures, as we shall discuss elsewhere. Several features of this approach, which we shall call the Thomas-Fermi-Bogoliubov description, are very attractive. (i) It is simpler than the Hartree-Fock-Bogoliubov technique. We can obtain analytical results in the case of weakly interacting bosons for quantities such as the chemical potential, the local depletion, pairing, pressure, and density of states. (ii) The method provides an estimate for the error due to the inhomogeneity of the Bose-condensed system. This error is a local quantity so that the validity of the description for a given trap and a given number of trapped atoms can be tested as a function of position. We see, for example, that at the edge of the condensate the Thomas-Fermi-Bogoliubov theory always breaks down. (iii) The Thomas-Fermi-Bogoliubov description can be generalized to treat the statistical mechanics of the Bose gas at finite temperatures.We derive a description of the spatially inhomogeneous Bose-Einstein condensate, treating the system locally as a homogeneous Bose gas. This approach, similar to the Thomas-Fermi model for the inhomogeneous many-particle fermion system, is well suited to describe the atomic Bose-Einstein condensates that have recently been obtained experimentally through atomic trapping and cooling. In this paper, we confine our attention to the zero-temperature case, although the treatment can be generalized to finite temperatures, as we shall discuss elsewhere. Several features of this approach, which we shall call the Thomas-Fermi-Bogoliubov description, are very attractive. (i) It is simpler than the Hartree-Fock-Bogoliubov technique. We can obtain analytical results in the case of weakly interacting bosons for quantities such as the chemical potential, the local depletion, pairing, pressure, and density of states. (ii) The method provides an estimate for the error due to the inhomogeneity of the Bose-condensed system. This error is a local quantity so that the validity of the description for a given trap and a given number of trapped atoms can be tested as a function of position. We see, for example, that at the edge of the condensate the Thomas-Fermi-Bogoliubov theory always breaks down. (iii) The Thomas-Fermi-Bogoliubov description can be generalized to treat the statistical mechanics of the Bose gas at finite temperatures.

arXiv:cond-mat/9609234Physical Review AMay 1, 199774 citations

A sharp gradient-corrected Lieb--Thirring inequality in one dimension with applications to density functional theory

Thiago Carvalho Corso, Louis Garrigue

We first prove a sharp Lieb-Thirring inequality with semi-classical constant and gradient correction in one dimension. In the periodic setting, the corresponding inequality also includes a parity dependent shell-correction term, which appears to be new. Moreover, we also characterize the optimizers of these inequalities. In particular, the optimizing density profiles form a two-dimensional manifold, which we call here the Poincaré density manifold. We then apply these results to study the Thomas--Fermi--von Weizsäcker (TFW) approximation of the kinetic energy in Density Functional Theory. More precisely, we first obtain an explicit and non-asymptotic bound on the global and local error made by the TFW approximation. These bounds are controlled by the distance of the normalized density profile from the Poincaré density manifold times the number of particles. In particular, these bounds show that the TFW approximation is exact on the Poincaré manifold. We then show that, in the periodic setting, the distance between the Kohn-Sham potential and the inverse TFW potential in the natural dual Sobolev space can also be controlled by explicit density dependent constants times the distance to the Poincaré manifold. As a corollary, we prove that for a fixed density profile, the Kohn-Sham potential converges towards the inverse Thomas--Fermi potential in the high-density limit $N \rightarrow \infty$ with asymptotic rate of order $N^{-2}$ . This rate is shown to be optimal.We first prove a sharp Lieb-Thirring inequality with semi-classical constant and gradient correction in one dimension. In the periodic setting, the corresponding inequality also includes a parity dependent shell-correction term, which appears to be new. Moreover, we also characterize the optimizers of these inequalities. In particular, the optimizing density profiles form a two-dimensional manifold, which we call here the Poincaré density manifold. We then apply these results to study the Thomas--Fermi--von Weizsäcker (TFW) approximation of the kinetic energy in Density Functional Theory. More precisely, we first obtain an explicit and non-asymptotic bound on the global and local error made by the TFW approximation. These bounds are controlled by the distance of the normalized density profile from the Poincaré density manifold times the number of particles. In particular, these bounds show that the TFW approximation is exact on the Poincaré manifold. We then show that, in the periodic setting, the distance between the Kohn-Sham potential and the inverse TFW potential in the natural dual Sobolev space can also be controlled by explicit density dependent constants times the distance to the Poincaré manifold. As a corollary, we prove that for a fixed density profile, the Kohn-Sham potential converges towards the inverse Thomas--Fermi potential in the high-density limit $N \rightarrow \infty$ with asymptotic rate of order $N^{-2}$. This rate is shown to be optimal.

math-phCross-listedarXiv:2610.30887Oct 23, 2026

Surface corrections to the Fermi-Thomas statistical theory

M.A. Naqvi

Nuclear PhysicsFeb 1, 195914 citations

d=2

and

3,

exact expressions for the N-particle densities are used to calculate perturbatively the temperature dependence of the splittings of the energy levels in a given shell due to a very weak interparticle interaction in a dilute Fermi gas. In two dimensions, we obtain analytically the surprising result that the l degeneracy in a harmonic-oscillator shell is not lifted in the lowest order even when the exact, rather than the Thomas-Fermi expression for the particle density is used. We also demonstrate rigorously (in two dimensions) the reduction of the exact zero-temperature fermionic expressions to the Thomas-Fermi form in the large-N limit.

N \rightarrow \infty

with asymptotic rate of order

N^{-2}

. This rate is shown to be optimal.