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.
We investigate the dynamic structure factor of a system of Bose particles at zero temperature using quantum Monte Carlo methods. Interactions are modeled using a hard-sphere potential of size
a and simulations are performed for values of the gas parameter
na3 ranging from the dilute regime up to densities
n where the thermodynamically stable phase is a solid. With increasing density, we observe a crossover of the dispersion of elementary excitations from a Bogoliubov-type spectrum to a phonon-maxon-roton curve and the emergence of a broad multiphonon contribution accompanying the single-quasiparticle peak. In particular, for
na3=0.2138, which corresponds to superfluid
4He at equilibrium density, the extracted spectrum turns out to be in good agreement with the experimental energy-momentum dispersion relation in the roton region and for higher momenta. The behavior of the spectral function at the same density in the stable solid and metastable gas phase above the freezing point is also discussed.
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.
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.
With the pair Hamiltonian model as the starting point, perturbation theory calculations are performed for the radial distribution function and the structure factor through second order in the two-body interaction potential for a many-boson system below the Bose-Einstein transition temperature. For the special case of a dilute hard-sphere Bose gas at T=0° K, the structure factor is calculated in the low-momentum limit keeping terms through second order in the hard-sphere gas parameter (n a3)1/2, where n is the gas density and a is the hard-sphere diameter. The resulting expression for the structure factor is shown to satisfy the Feynman-Bijl relation explicitly.
The zero-temperature dynamical structure factor $S(q,\ensuremath{\omega})$ of one-dimensional hard rods is computed using state-of-the-art quantum Monte Carlo and analytic continuation techniques, complemented by a Bethe ansatz analysis. As the density increases, $S(q,\ensuremath{\omega})$ reveals a crossover from the Tonks-Girardeau gas to a quasisolid regime, along which the low-energy properties are found in agreement with the nonlinear Luttinger liquid theory. Our quantitative estimate of $S(q,\ensuremath{\omega})$ extends beyond the low-energy limit and confirms a theoretical prediction regarding the behavior of $S(q,\ensuremath{\omega})$ at specific wave vectors ${\mathcal{Q}}_{n}=n2\ensuremath{\pi}/a$, where
a is the core radius, resulting from the interplay of the particle-hole boundaries of suitably rescaled ideal Fermi gases. We observe significant similarities between hard rods and one-dimensional
4He at high density, suggesting that the hard-rods model may provide an accurate description of dense one-dimensional liquids of quantum particles interacting through a strongly repulsive, finite-range potential.