Abstract We use a combinatorial approach to obtain exact expressions for the many-body density of states of fermionic and bosonic gases with equally spaced single-particle spectra. We identify a mapping that reveals a remarkable property, namely, fermionic and bosonic gases have the same many-body density of states, up to a shift corresponding to ground state energy. Additionally, we show that there is a regime, comprising the validity range of the Bethe approximation, where the many-body density of states becomes independent of the number of particles.
As a preliminary to discussing the properties of Bose–Einstein and Fermi–Dirac gases, this chapter presents the basic quantum equations that underlie both. The process of building N-body states from single-particle states is described, along with notation necessary to simplify the problem. The average number of particles in a given energy state and chemical potential is derived for fermions and bosons. Although all atoms are either bosons or fermions, and therefore indistinguishable, there are nevertheless real systems that are composed of distinguishable particles. Although all atoms are either bosons or fermions, and therefore indistinguishable, there are nevertheless real systems that are composed of distinguishable particles.
The ground state energy and elementary excitations of a charged gas of bosons at high densities are examined by use of the method developed by Bogoliubov for boson gases. It is conjectured, but not herein established, that this method yields exact results in the high-density limit analogous to those obtained by Gell-Mann and Brueckner, and by Sawada, in the corresponding case of a charged fermion gas. The ground state energy is essentially correlation energy, and is therefore negative, and its magnitude varies as the one-fourth power of the density at high densities. The elementary excitations have for low momenta the energy appropriate to plasma waves, and for high momenta the energy appropriate to single-particle excitation. There is therefore an energy gap, suggesting that the gas is both a superfluid and a superconductor at low temperatures. At low densities the behavior of a charged gas is independent of statistics; hence, such a gap must disappear as the system is expanded.
Bose–Einstein and Fermi–Dirac statistics versus Boltzmann statistics It was remarked in Chapter 2, §2.1, that there are circumstances in which the “Boltzmann statistics” for an ideal gas – the approximation in Eq. (2.6), in which the total partition function is decomposed into a product of single-molecule partition functions, with division by N 1 ! N 2 ! … to correct for particle identity – is inaccurate even if the gas is still ideal (non-interacting particles). That happens when the number of single-particle states within the energy kT of the ground state is no longer very much greater than the number of particles. (Recall Fig. 2.1 and the accompanying discussion.) We noted that when the Boltzmann statistics is no longer an adequate approximation one must distinguish the particles as fermions, no two of which can be in the same single-particle state (counting states of different spin as different), or as bosons, for which there is no restriction on multiple occupancy of the singleparticle states. Fermions are particles of half-odd-integer spin quantum number; they could be composite particles consisting of an odd number of spin-½ particles. Protons, neutrons, and electrons, all of which are of spin ½, are fermions. Bosons are particles of integer spin; they could be spinless, or of spin 1, etc., or could be composites consisting of an even number of fermions. The nuclei of helium atoms of mass number 4 consist of two protons and two neutrons, and the neutral atom has two electrons outside the nucleus; 4 He atoms are bosons.
One-dimensional spinless Bose and Fermi gases with contact interactions have the close interrelation via Girardeau's Bose-Fermi mapping, leading to the correspondences in their energy spectra and thermodynamics. However, correlation functions are in general not identical between these systems. We derive in both systems the exact universal relations for correlation functions, which hold for any energy eigenstate and any statistical ensemble of the eigenstates with or without a trapping potential. These relations include the large-momentum behaviors of static structure factors and of momentum distributions as well as energy relations, which connect the sums of kinetic and interaction energies to the momentum distributions. The relations involve two- and three-body contacts, which are the integrals of local pair and triad correlations, respectively. We clarify how the relations for bosons and fermions differ and are connected with each other. In particular, we find that the three-body contact makes no contribution to the bosonic energy relation, but it plays a crucial role in the fermionic one. In addition, we compute the exact momentum distribution for any number of fermions in the unitary limit.
Many useful properties of dilute Bose gases at ultralow temperature are predicted precisely by the (mean-field) product-state Ansatz, in which all particles are in the same quantum state. Yet, in situations where particle-particle correlations become important, the product Ansatz fails. To include correlations nonperturbatively, we consider a new set of states: the particle-correlated state of $N=l\ifmmode\times\else\texttimes\fi{}n$ bosons is derived by symmetrizing the
n-fold product of an
l-particle quantum state. Quantum correlations of the
l-particle state ``spread out'' to any subset of the
N bosons by symmetrization. The particle-correlated states can be simulated efficiently for large
N, because their parameter spaces, which depend on
l, do not grow with
n. Here we formulate and develop in great detail the pure-state case for
l=2, where the many-body state is constructed from a two-particle pure state. These paired wave functions, which we call pair-correlated states (PCS), were introduced by A. J. Leggett [Rev. Mod. Phys. 73, 307 (2001)] as a particle-number-conserving version of the Bogoliubov approximation. We present an iterative algorithm that solves for the reduced (marginal) density matrices (RDMs), i.e., the correlation functions, associated with PCS in time
O(N). The RDMs can also be derived from the normalization factor of PCS, which is derived analytically in the large-
N limit. To test the efficacy of PCS, we analyze the ground state of the two-site Bose-Hubbard model by minimizing the energy of the PCS state, both in its exact form and in its large-
N approximate form, and comparing with the exact ground state. For
N=1000, the relative errors of the ground-state energy for both cases are within ${10}^{\ensuremath{-}5}$ over the entire parameter region from a single condensate to a Mott insulator. We present numerical results that suggest that PCS might be useful for describing the dynamics in the strongly interacting regime.