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

Journal of Physics A Mathematical and General

The determination of the Dirac density matrix of the d-dimensional harmonic oscillator for an arbitrary number of closed shells

I. A. Howard, N. H. March, Luis Miguel Nieto

Published Jun 8, 2002• 8 citations
Publisher

Abstract

In 1959, March and Young ( Nucl. Phys. 12 237) rewrote the equation of motion for the Dirac density matrix γ( x , x 0 ) in terms of sum and difference variables. Here, γ( , 0 ) for the d -dimensional isotropic harmonic oscillator for an arbitrary number of closed shells is shown to satisfy, using the variables ∣ + 0 ∣/2 and ∣ − 0 ∣/2, a generalized partial differential equation embracing the March–Young equation for d = 1. As applications, we take in turn the cases d = 1, 2, 3 and 4, and obtain both the density matrix γ( , 0 ) and the diagonal density ϱ( r ) = γ( , 0 )∣ 0 = , this diagonal element already being known to satisfy a third-order linear homogeneous differential equation for d = 1 through 3. Some comments are finally made on the d -dimensional kinetic energy density, which is important for first-principles density functional theory in allowing one to bypass one-particle Schrödinger equations (the so-called Slater–Kohn–Sham equations).

DOI
10.1088/0305-4470/35/24/302

Related papers

Density and near-diagonal Dirac density matrix for closed shells of isotropic harmonically confined fermions in three dimensions

I. A. Howard, Norman Henry March

In early work, Lawes and March obtained a differential equation for an arbitrary number of independent harmonically confined fermions in one dimension, and very recently this result has been generalized to apply to three-dimensional (3D) isotropic harmonic confinement. Here, an exact solution of this 3D equation for the fermion particle density ρ( r ) is constructed, and the near-diagonal form of the Dirac density matrix is also obtained.

Journal of Physics A Mathematical and General•Sep 3, 2001•11 citations

Single-particle kinetic energy in density-functional theory: Harmonic confinement in two and three dimensions

N. H. March

Using the known perturbation series for the idempotent Dirac density matrix in powers of a given one-body potential $V(\mathbf{r})$, Stoddart and March (SM) generated a corresponding series for the kinetic energy density. While the general term of the SM series is known, a summation has not been achieved to date. A contribution to solve this problem is made here by exhibiting an explicit form of the above kinetic energy density for Fermions filling an arbitrary number of closed shells, when the confinement is harmonic. This example is of considerable current interest because of ongoing experiments on ultracold atomic gases of Fermions.

Physical Review A•Apr 20, 2006•2 citations

Differential equation for the Dirac single-particle first-order density matrix in terms of the ground-state electron density

Norman Henry March, Tamás Gál, I. A. Howard

The March-Suhai (MS) partial differential equation for the Dirac density matrix ${\ensuremath{\gamma}}_{s}(\stackrel{\ensuremath{\rightharpoonup}}{r},{\stackrel{\ensuremath{\rightharpoonup}}{r}}^{\ensuremath{'}})$, proved for one- and two-level occupancies, involves both the ground-state density $n(\stackrel{\ensuremath{\rightharpoonup}}{r})$, with its low-order derivatives, and the positive definite kinetic energy density ${t}_{s}(\stackrel{\ensuremath{\rightharpoonup}}{r})$. Here, we examine the relation between the equation of motion for ${\ensuremath{\gamma}}_{s}(\stackrel{\ensuremath{\rightharpoonup}}{r},{\stackrel{\ensuremath{\rightharpoonup}}{r}}^{\ensuremath{'}})$, with input now being the one-body potential of density-functional theory, and the MS equation. The important link is the differential virial theorem, which can be used to remove ${t}_{s}(\stackrel{\ensuremath{\rightharpoonup}}{r})$ from the MS differential equation. For multiple occupancy, the Pauli potential enters in an important manner. In one dimension, however, the appearance of the Pauli potential can be avoided, obtaining a necessary condition for ${\ensuremath{\gamma}}_{s}(x,{x}^{\ensuremath{'}})$ to satisfy for arbitrary level occupancy, in the form of a MS-type differential equation.

Physical Review A•Jun 22, 2010•1 citation

Complete functional theory for the fermion density of independent particles subject to harmonic confinement in d dimensions for an arbitrary number of closed shells

I. A. Howard, N. H. March, Luis Miguel Nieto

In earlier work, expressions have been constructed for the single-particle kinetic-energy functional ${T}_{s}[\ensuremath{\rho}]$ for independent fermions subject to harmonic confinement in low dimensions, with $\ensuremath{\rho}$ the particle density. Here, the differential equation for $\ensuremath{\rho}$ is first obtained in d dimensions for an arbitrary number of closed shells. Then, by using the known Euler-Lagrange equation, the functional derivative $\ensuremath{\delta}{T}_{s}/\ensuremath{\delta}\ensuremath{\rho}(\mathbf{r})$ is constructed. ${T}_{s}[\ensuremath{\rho}]$ itself is proved to take the form of a linear combination of three pieces: (i) a von Weizs\"acker inhomogeneity kinetic energy, but with the original coefficient reduced by a dimensionality factor $1/d,$ (ii) a Thomas-Fermi kinetic energy, and (iii) a truly nonlocal contribution that, however, is shown to involve only the density $\ensuremath{\rho}$ itself and its first derivative. Thus, for this model, which is currently highly relevant to the interpretation of experiments on the evaporative cooling of dilute, and hence almost noninteracting, fermions, a complete density-functional theory now exists.

Physical Review A•Nov 22, 2002•25 citations

Dirac density matrix and the Legendre transform of the kinetic energy generated by one-dimensional model potentials

I. A. Howard, Norman Henry March

In density functional theory, the single-particle kinetic energy is still not known in an orbital-free form. The natural tool to calculate such kinetic energy for a given one-body potential V is the Dirac density matrix y. Here, by first taking solvable one-dimensional examples, and in particular the sech 2 (x) potential, forms for the Dirac density matrix are exhibited in terms of the potential.This in turn generates the Legendre transform of the single-particle kinetic energy. Finally, one-dimensional perturbation theory for the kinetic energy density, summed to all orders in V, is presented in the appendix.

Journal of Physics A Mathematical and General•Dec 21, 2005•0 citations

Harmonic oscillator thermal density matrix: First-order differential equations for the position representation

Luis F. Barragán-Gil, R. Walser

A first-order partial differential equation is derived whose solution enables us to find straightforwardly the off-diagonal matrix elements in the position representation of the harmonic oscillator density operator. This approach constitutes an alternative to techniques that require advanced knowledge of mathematical and quantum mechanical results.

American Journal of Physics•Dec 19, 2017•12 citations
In 1959, March and Young ( Nucl. Phys. 12 237) rewrote the equation of motion for the Dirac density matrix γ( x , x 0 ) in terms of sum and difference variables. Here, γ( , 0 ) for the d -dimensional isotropic harmonic oscillator for an arbitrary number of closed shells is shown to satisfy, using the variables ∣ + 0 ∣/2 and ∣ − 0 ∣/2, a generalized partial differential equation embracing the March–Young equation for d = 1. As applications, we take in turn the cases d = 1, 2, 3 and 4, and obtain both the density matrix γ( , 0 ) and the diagonal density ϱ( r ) = γ( , 0 )∣ 0 = , this diagonal element already being known to satisfy a third-order linear homogeneous differential equation for d = 1 through 3. Some comments are finally made on the d -dimensional kinetic energy density, which is important for first-principles density functional theory in allowing one to bypass one-particle Schrödinger equations (the so-called Slater–Kohn–Sham equations).
Using the known perturbation series for the idempotent Dirac density matrix in powers of a given one-body potential V(r)V(\mathbf{r})V(r), Stoddart and March (SM) generated a corresponding series for the kinetic energy density. While the general term of the SM series is known, a summation has not been achieved to date. A contribution to solve this problem is made here by exhibiting an explicit form of the above kinetic energy density for Fermions filling an arbitrary number of closed shells, when the confinement is harmonic. This example is of considerable current interest because of ongoing experiments on ultracold atomic gases of Fermions.
The March-Suhai (MS) partial differential equation for the Dirac density matrix ${\ensuremath{\gamma}}_{s}(\stackrel{\ensuremath{\rightharpoonup}}{r},{\stackrel{\ensuremath{\rightharpoonup}}{r}}^{\ensuremath{'}})$, proved for one- and two-level occupancies, involves both the ground-state density $n(\stackrel{\ensuremath{\rightharpoonup}}{r})$, with its low-order derivatives, and the positive definite kinetic energy density ${t}_{s}(\stackrel{\ensuremath{\rightharpoonup}}{r})$. Here, we examine the relation between the equation of motion for ${\ensuremath{\gamma}}_{s}(\stackrel{\ensuremath{\rightharpoonup}}{r},{\stackrel{\ensuremath{\rightharpoonup}}{r}}^{\ensuremath{'}})$, with input now being the one-body potential of density-functional theory, and the MS equation. The important link is the differential virial theorem, which can be used to remove ${t}_{s}(\stackrel{\ensuremath{\rightharpoonup}}{r})$ from the MS differential equation. For multiple occupancy, the Pauli potential enters in an important manner. In one dimension, however, the appearance of the Pauli potential can be avoided, obtaining a necessary condition for ${\ensuremath{\gamma}}_{s}(x,{x}^{\ensuremath{'}})$ to satisfy for arbitrary level occupancy, in the form of a MS-type differential equation.
In earlier work, expressions have been constructed for the single-particle kinetic-energy functional ${T}_{s}[\ensuremath{\rho}]$ for independent fermions subject to harmonic confinement in low dimensions, with $\ensuremath{\rho}$ the particle density. Here, the differential equation for $\ensuremath{\rho}$ is first obtained in d dimensions for an arbitrary number of closed shells. Then, by using the known Euler-Lagrange equation, the functional derivative $\ensuremath{\delta}{T}_{s}/\ensuremath{\delta}\ensuremath{\rho}(\mathbf{r})$ is constructed. ${T}_{s}[\ensuremath{\rho}]$ itself is proved to take the form of a linear combination of three pieces: (i) a von Weizs\"acker inhomogeneity kinetic energy, but with the original coefficient reduced by a dimensionality factor 1/d,1/d,1/d, (ii) a Thomas-Fermi kinetic energy, and (iii) a truly nonlocal contribution that, however, is shown to involve only the density $\ensuremath{\rho}$ itself and its first derivative. Thus, for this model, which is currently highly relevant to the interpretation of experiments on the evaporative cooling of dilute, and hence almost noninteracting, fermions, a complete density-functional theory now exists.
In density functional theory, the single-particle kinetic energy is still not known in an orbital-free form. The natural tool to calculate such kinetic energy for a given one-body potential V is the Dirac density matrix y. Here, by first taking solvable one-dimensional examples, and in particular the sech 2 (x) potential, forms for the Dirac density matrix are exhibited in terms of the potential.This in turn generates the Legendre transform of the single-particle kinetic energy. Finally, one-dimensional perturbation theory for the kinetic energy density, summed to all orders in V, is presented in the appendix.