By using exact quantum Monte Carlo methods we calculate the ground-state properties of the liquid phase in one-dimensional Bose mixtures with contact interactions. We find that the liquid state can be formed if the ratio of coupling strengths between interspecies attractive and intraspecies repulsive interactions exceeds a critical value. As a function of this ratio we determine the density where the energy per particle has a minimum and the one where the compressibility diverges, thereby identifying the equilibrium density and the spinodal point in the phase diagram of the homogeneous liquid. Furthermore, in the stable liquid state, we calculate the chemical potential, the speed of sound, as well as structural and coherence properties, such as the pair correlation function, the static structure factor, and the one-body density matrix, thus providing a detailed description of the bulk region in self-bound droplets.
The method of self-consistent-field approximation is applied for the study of the ground-state properties of a one-dimensional system of bosons, of density $\ensuremath{\rho}$, first interacting through finite rigid rods, of length
l, and again via a repulsive $\ensuremath{\delta}$-function potential. We obtain the energy spectrum and the sound velocity. The results are close to the exact one for $\ensuremath{\rho}l<0.5$. The energy spectrum obtained shows a rotonlike region which is qualitatively similar to the two- and three-dimensional versions. For a Bose gas interacting through a repulsive $\ensuremath{\epsilon}$ potential of strength
2c, we derive the energy spectrum, the sound velocity, and the ground-state energy. The results are compared with Girardeau and Lieb's exact expressions. In the small-coupling limit our calculations show the same results as obtained by Bogoliubov's perturbation theory. The ground-state energy obtained for small $\ensuremath{\gamma}=\frac{c}{\ensuremath{\rho}}$ is found to agree closely with several different works found in the literature. Unfortunately our theory diverges for the ground-state energy in the limit $\ensuremath{\gamma}\ensuremath{\rightarrow}\ensuremath{\infty}$. The sound velocity in this limit is finite but twice as small as from Girardeau's exact results.
We continue the analysis of the one-dimensional gas of Bose particles interacting via a repulsive delta function potential by considering the excitation spectrum. Among other things we show that: (i) the elementary excitations are most naturally thought of as a double spectrum, not a single one; (ii) the velocity of sound derived from the macroscopic compressibility is shown to agree with the velocity of sound derived from microscopic considerations, i.e., from the phonon spectrum. We also introduce a distinction between elementary excitations and quasiparticles, on the basis of which we give some heuristic reasons for expecting the double spectrum to be a general feature, even in three dimensions, and not an exception.
We study the correlation effects in a one-dimensional Bose gas with repulsive delta-function interaction. The correlation effects are described by a local-field correction which takes into account the short-range correlations. We find that the ground state energy is in good agreement with the exact result up to intermediate coupling strengths, showing an improvement over the Bogoliubov and other related approximations. The velocity of sound, the static structure factor, and pair-correlation function are also calculated within the present approximation.
The finite-temperature properties of Bose quantum fluids are studied by applying a variational density-matrix approach. Our ansatz for the density matrix is similar to that introduced by Campbell, K\"urten, Ristig, and Senger [Phys. Rev. 30, 3728 (1984)] but has been generalized to include three-body correlations. An expression for the finite-temperature elementary-excitation energy, \ensuremath{\epsilon}(k), is derived from an ansatz for the variational entropy. At zero temperature \ensuremath{\epsilon}(k) reduces to the same expression obtained by Brillouin-Wigner perturbation theory. The hypernetted-chain (HNC) equations provide the relations between the n-particle distribution functions and the correlation functions. A substantial set of elementary diagrams have been included to minimize the loss of precision associated with the HNC approximation. We report values for the static structure function, radial distribution function, excitation energies, particle-hole interaction, Helmholtz free energy, entropy, and isothermal sound velocity. We present a detailed study of the spectra, three-body correlations, and elementary diagram dependence of the thermodynamic spinodal curve. For relevant densities, we find that three-body correlations tend to increase the liquid-gas critical temperature
Tc while adding elementary diagrams tends to reduce
Tc. Including multiphonon terms in the spectrum produces several features in the variational entropy that were not observed at the Feynman level of excitations. Entropy isotherms show a rather abrupt transition from a low-density regime, where multiphonon contributions are insignificant, to a high-density regime where the isotherms are strongly influenced by multiphonon scattering processes.
We report diffusion Monte Carlo (DMC) and path integral Monte Carlo (PIMC) calculations of the properties of a one-dimensional (1D) Bose quantum fluid. The equation of state, the superfluid fraction ${\ensuremath{\rho}}_{S}/{\ensuremath{\rho}}_{0}$, the one-body density matrix
n(x), the pair distribution function
g(x), and the static structure factor
S(q) are evaluated. The aim is to test Luttinger liquid (LL) predictions for 1D fluids over a wide range of fluid density and LL parameter
K. The 1D Bose fluid examined is a single chain of
4He atoms confined to a line in the center of a narrow nanopore. The atoms cannot exchange positions in the nanopore, the criterion for 1D. The fluid density is varied from the spinodal density where the 1D liquid is unstable to droplet formation to the density of bulk liquid
4He. In this range,
K varies from
K>2 at low density, where a robust superfluid is predicted, to
K<0.5, where fragile 1D superflow and solidlike peaks in
S(q) are predicted. For uniform pore walls, the ${\ensuremath{\rho}}_{S}/{\ensuremath{\rho}}_{0}$ scales as predicted by LL theory. The
n(x) and
g(x) show long range oscillations and decay with
x as predicted by LL theory. The amplitude of the oscillations is large at high density (small
K) and small at low density (large
K). The
K values obtained from different properties agree well verifying the internal structure of LL theory. In the presence of disorder, the ${\ensuremath{\rho}}_{S}/{\ensuremath{\rho}}_{0}$ does not scale as predicted by LL theory. A single
vJ parameter in the LL theory that recovers LL scaling was not found. The one body density matrix (OBDM) in disorder is well predicted by LL theory. The ``dynamical'' superfluid fraction, ${\ensuremath{\rho}}_{S}^{D}/{\ensuremath{\rho}}_{0}$, is determined. The physics of the deviation from LL theory in disorder and the ``dynamical'' ${\ensuremath{\rho}}_{S}^{D}/{\ensuremath{\rho}}_{0}$ are discussed.