The superfluid to normal fluid transition of dipolar bosons in two dimensions is studied in a broad density range by using path integral Monte Carlo simulations and summarized in the phase diagram. While at low densities we find good agreement with the universal results depending only on the scattering length a{s}, at moderate and high densities the transition temperature is strongly affected by interactions and the excitation spectrum of quasiparticles. The results are expected to be of relevance to dipolar atomic and molecular systems and indirect excitons in quantum wells.
The superfluid to normal fluid transition of dipolar bosons in two dimensions is studied in a broad density range by using path integral Monte Carlo simulations and summarized in the phase diagram. While at low densities we find good agreement with the universal results depending only on the scattering length a{s}, at moderate and high densities the transition temperature is strongly affected by interactions and the excitation spectrum of quasiparticles. The results are expected to be of relevance to dipolar atomic and molecular systems and indirect excitons in quantum wells.
A two-dimensional quantum system of dipoles, with a polarization angle not perpendicular to the plane, shows a transition from a gas to a stripe phase. We have studied the thermal properties of these two phases using the path-integral Monte Carlo (PIMC) method. By simulating the thermal density matrix, PIMC provides exact results for magnitudes of interest such as the superfluid fraction and the one-body density matrix. As it is well known, in two dimensions the superfluid-to-normal phase transition follows the Berezinskii-Kosterlitz-Thouless (BKT) scenario. Our results show that both the anisotropic gas and the stripe phases follow the BKT scaling laws. At fixed density and increasing the tilting angle, the transition temperature decreases in going from the gas to the stripe phase. Superfluidity in the perpendicular direction to the stripes is rather small close to the critical temperature but it becomes larger at lower temperatures, mainly close to the transition to the gas. Our results are in qualitative agreement with the supersolidity observed recently in a quasi-one-dimensional array of dipolar droplets.
The Berezinskii-Kosterlitz-Thouless transition in two-dimensional dipolar systems has been studied recently by path integral Monte Carlo simulations [A. Filinov et al., Phys. Rev. Lett. 105, 070401 (2010)]. Here, we complement this analysis and study temperature-coupling strength dependence of the density (particle-hole) and single-particle (SP) excitation spectra both in superfluid and normal phases. The dynamic structure factor, $S(q,\ensuremath{\omega})$, of the longitudinal excitations is rigorously reconstructed with full information on damping. The SP spectral function, $A(q,\ensuremath{\omega})$, is worked out from the one-particle Matsubara Green's function. A stochastic optimization method is applied for reconstruction from imaginary times. In the superfluid regime sharp energy resonances are observed in both the density and SP excitations. The involved hybridization of both spectra is discussed. In contrast, in the normal phase, when there is no coupling, the density modes, beyond acoustic phonons, are significantly damped. Our results generalize previous zero-temperature analyses based on variational many-body wave functions [F. Mazzanti et al., Phys. Rev. Lett. 102, 110405 (2009); D. Hufnagl et al., Phys. Rev. Lett. 107, 065303 (2011)], where the underlying physics of the excitation spectrum and the role of the condensate has not been addressed.
The supersolid phase of soft-core bosons in two dimensions is investigated using the self-consistent Hartree-Fock and quantum Monte Carlo methods. An approximate phase diagram at finite temperatures is initially constructed using the mean-field approach, which is subsequently validated through precise path-integral simulations, enabling a microscopic characterization of the various phases. Superfluid and melting/freezing transitions are analyzed through the superfluid density and the long-range behavior of correlation functions associated with positional and orientational order, in accordance with the general picture of Berezinskii-Kosterlitz-Thouless transitions. A broad region at low temperatures is identified where the supersolid phase exists, separating the uniform superfluid phase from the normal quasicrystal phase. Additionally, a potential intermediate hexatic phase with quasi-long-range orientational order is identified in a narrow region between the normal solid and fluid phases. These findings establish self-consistent Hartree-Fock theory beyond the local density approximation as an effective tool, complementary to computationally intensive quantum Monte Carlo simulations, for investigating the melting of the supersolid phase and the possible emergence of the hexatic superfluid phase in bosonic systems with various interaction potentials.
Long-range and anisotropic dipolar interactions induce complex order in quantum systems, particularly interesting in 2D where superfluidity emerges via the Berezinskii-Kosterlitz-Thouless (BKT) mechanism. We observe a superfluid phase with algebraically decaying correlation in a quasi-2D dipolar Bose gas of erbium atoms and further identify the superfluid phase transition by monitoring extended coherence. We also measure equations of state with tunable dipolar interactions. Our findings show a transition point shift due to dipole orientation, aligning with BKT transition with effective short-range interaction. However, we observe that in-plane tilted dipoles show nonlocal effects in the superfluid regime, preventing the effective short-range interaction from revealing universal behavior near the BKT critical point. We also measure anisotropic atom number fluctuations in the superfluid regime, highlighting its dipolar nature. Our results provide a foundation for understanding dipolar bosons in 2D and pave the way for investigating complex orders in dipolar superfluids.
We study the universal critical behavior of two-dimensional (2D) lattice bosonic gases at the Berezinskii-Kosterlitz-Thouless (BKT) transition, which separates the low-temperature superfluid phase from the high-temperature normal phase. For this purpose, we perform quantum Monte Carlo simulations of the hard-core Bose-Hubbard (BH) model at zero chemical potential. We determine the critical temperature by using a matching method that relates finite-size data for the BH model with corresponding data computed in the classical XY model. In this approach, the neglected scaling corrections decay as inverse powers of the lattice size L, and not as powers of 1/lnL, as in more standard approaches, making the estimate of the critical temperature much more reliable. Then, we consider the BH model in the presence of a trapping harmonic potential, and we verify the universality of the trap-size dependence at the BKT critical point. This issue is relevant for experiments with quasi-2D trapped cold atoms.
We calculate the superfluid transition temperature of homogeneous interacting Bose gases in three and two spatial dimensions using large-scale path integral Monte Carlo simulations (with up to N=10;{5} particles). In 3D we investigate the limits of the universal critical behavior in terms of the scattering length alone by using different models for the interatomic potential. We find that this type of universality sets in at small values of the gas parameter na3 < or approximately 10(-4). This value is different from the estimate na3 < or approximately 10(-6) for the validity of the asymptotic expansion in the limit of vanishing na3. In 2D we study the Berezinskii-Kosterlitz-Thouless transition of a gas with hard-core interactions. For this system we find good agreement with the classical lattice |psi|4 model up to very large densities. We also explain the origin of the existing discrepancy between previous studies of the same problem.