Exact first-order effect of interactions on the ground-state energy of harmonically-confined fermions
Pierre Le Doussal,Naftali R. Smith,Nathan Argaman
DOI: https://doi.org/10.21468/scipostphys.17.2.038
2024-08-07
SciPost Physics
Abstract:We consider a system of N N spinless fermions, interacting with each other via a power-law interaction \epsilon/r^n ε / r n , and trapped in an external harmonic potential V(r) = r^2/2 V ( r ) = r 2 / 2 , in d=1,2,3 d = 1 , 2 , 3 dimensions. For any 0 < n < d+2 0 < n < d + 2 , we obtain the ground-state energy E_N E N of the system perturbatively in \epsilon ε , E_{N}=E_{N}^{≤ft(0)}+\epsilon E_{N}^{≤ft(1)}+O≤ft(\epsilon^{2}) E N = E N ≤ f t ( 0 ) + ε E N ≤ f t ( 1 ) + O ≤ f t ( ε 2 ) . We calculate E_{N}^{≤ft(1)} E N ≤ f t ( 1 ) exactly, assuming that N N is such that the “outer shell” is filled. For the case of n=1 n = 1 (corresponding to a Coulomb interaction for d=3 d = 3 ), we extract the N \gg 1 N ≫ 1 behavior of E_{N}^{≤ft(1)} E N ≤ f t ( 1 ) , focusing on the corrections to the exchange term with respect to the leading-order term that is predicted from the local density approximation applied to the Thomas-Fermi approximate density distribution. The leading correction contains a logarithmic divergence, and is of particular importance in the context of density functional theory. We also study the effect of the interactions on the fermions’ spatial density. Finally, we find that our result for E_{N}^{≤ft(1)} E N ≤ f t ( 1 ) significantly simplifies in the case where n n is even.
physics, multidisciplinary