We obtain ab-initio estimations of the dynamic structure factor, S ( q ,ω), of Bose gases at zero temperature. More precisely, we use the Genetic Inversion via Falsification of Theories (GIFT) algorithm to perform analytic continuations of imaginary time correlation functions computed via an exact Path Integral projector method. Using the hard-sphere potential to model the two-body interactions between the atoms, we compute S ( q ,ω) changing the gas parameter from the dilute regime ( na 3 = 10 −4 ) up to the density corresponding to superfluid 4 He at equilibrium ( na 3 = 0.2138). With increasing density, we observe the emergence of a broad multiphonon contribution accompanying the quasiparticle peak and a crossover of the dispersion of elementary excitations from a Bogoliubov-like spectrum to a phonon-maxon- roton curve. Apart from the low wave vector region, for na 3 = 0.2138 the energy-momentum dispersion relation and the static density response function, χ( q ), turns out to be in good agreement with the superfluid 4 He experimental data at equilibrium density.
The properties of Bose-Einstein-condensed systems of hard spheres at high density are studied by taking into account higher-order corrections beyond the Hartree-Fock-Bogoliubov approximation. Based on a diagrammatic technique of the perturbation expansion starting from the Bogoliubov approximation, we show that the effective interactions separate into two distinct parts: the screened and condensate-mediated interactions. The latter is described by the exchange of virtual quasiparticles between noncondensate bosons via the Bose-Einstein condensate and is attractive for small energy transfer. We calculate the condensation temperature, excitation spectrum, and condensate density in the whole region of the gas parameter. It is shown that the condensate temperature decreases compared to that of an ideal Bose gas due to hard-sphere interactions. We find that the condensate-mediated attraction between noncondensate bosons leads to an enhancement of the anomalous self-energy and, as a consequence, to the emergence of a roton minimum in the excitation spectrum and to the strong suppression of the condensate fraction in the region of high density.
The thermodynamics of a free Bose gas with effective temperature scale
T~ and hard-sphere Bose gas with the
T~ scale are studied.
T~ arises as the temperature experienced by a single particle in a quantum gas with 2-body harmonic oscillator interaction
Vosc, which at low temperatures is expected to simulate, almost correctly, the attractive part of the interatomic potential
VHe between
4He atoms. The repulsive part of
VHe is simulated by a hard-sphere (HS) potential. The thermodynamics of this system of HS bosons, with the
T~ temperature scale (HSET), is investigated, first, by the Bogoliubov-Huang method and next by a modified version of this method, which takes approximate account of those terms of the 2-body repulsion which are linear in the zero-momentum Bose operators
a0,a0∗ (originally rejected by Bogoliubov). Theoretical heat capacity
CV(T) exhibits good agreement, below 2.1 K, with the experimental heat capacity graph observed in
4He at saturated vapour pressure. The phase transition to the low-temperature phase, with a Bose-Einstein condensate, occurs in the HSET at $T_λ=$2.17 K, and is accompanied, in the modified HSET version, by a singularity of
CV(T). Other thermal properties of HSET, such as the momentum distribution function, the fraction of atoms in the momentum condensate and normal fluid density, agree qualitatively with those of
4He, but improve those of the free Bose gas.
Using a method due to the author which correctly reproduces the scattering of two hard spheres for all relative momenta and all partial waves, we study the properties of a low-density hard-sphere boson gas at
T=0\ifmmode^\circ\else\textdegree\fi{}K. The system is described by a Hamiltonian
H defined for all relative distances of the
N particles. The terms of
H corresponding to the forward motion of a particle through the system lead to the notion of an effective mass
m∗ for each particle in excess of its actual mass
m, with
mm∗ given by a monotonically increasing function of $\ensuremath{\rho}{a}^{3}$, where $\ensuremath{\rho}$ and
a are the number of particles per unit volume and the hard-sphere diameter, respectively. The excitation spectrum is linear at low momenta, displays a roton-like behavior for momenta $p\ensuremath{\gtrsim}\frac{\ensuremath{\hbar}}{a}$, possesses further extrema at higher momenta, and asymptotically approaches
2m∗p2. The pair distribution function vanishes in the low-density limit for particle separations $r\ensuremath{\le}a$.
The energy and structure of dilute hard- and soft-sphere Bose gases are systematically studied in the framework of several many-body approaches, such as the variational correlated theory, the Bogoliubov model, and the uniform limit approximation, valid in the weak-interaction regime. When possible, the results are compared with the exact diffusion Monte Carlo ones. Jastrow-type correlation provides a good description of the systems, both hard- and soft-spheres, if the hypernetted chain energy functional is freely minimized and the resulting Euler equation is solved. The study of the soft-sphere potentials confirms the appearance of a dependence of the energy on the shape of the potential at gas paremeter values of $x\ensuremath{\sim}0.001.$ For quantities other than the energy, such as the radial distribution functions and the momentum distributions, the dependence appears at any value of x. The occurrence of a maximum in the radial distribution function, in the momentum distribution, and in the excitation spectrum is a natural effect of the correlations when x increases. The asymptotic behaviors of the functions characterizing the structure of the systems are also investigated. The uniform limit approach is very easy to implement and provides a good description of the soft-sphere gas. Its reliability improves when the interaction weakens.
We study a bosonic gas of hard spheres by using the exact zero-temperature path-integral ground-state (PIGS) Monte Carlo method and the equations of superfluid hydrodynamics. The PIGS method is implemented to calculate for the bulk system the energy per particle and the condensate fraction through a large range of the gas parameter
na3 (with
n the number density and
a the
s-wave scattering length), going from the dilute gas into the solid phase. The Maxwell construction is then adopted to determine the freezing at $n{a}^{3}=0.264\ifmmode\pm\else\textpm\fi{}0.003$ and the melting at $n{a}^{3}=0.290\ifmmode\pm\else\textpm\fi{}0.003$. In the liquid phase, where the condensate fraction is finite, the equations of superfluid hydrodynamics, based on the PIGS equation of state, are used to find other relevant quantities as a function of the gas parameter: the chemical potential, the pressure, and the sound velocity. In addition, within Feynman's approximation, from the PIGS static structure factor we determine the full excitation spectrum, which displays a maxon-roton behavior when the gas parameter is close to the freezing value. Finally, the equations of superfluid hydrodynamics with the PIGS equation of state are solved for the bosonic system under axially symmetric harmonic confinement obtaining its collective breathing modes.