Finite-size-scaling analysis of a simulation of the 4He superfluid transition

Density functional theory of superfluid helium at finite temperatures

A density functional theory-based method is developed to describe the static and dynamic response of superfluid helium at finite temperatures. The model relies on the well-established 0 K Orsay-Trento functional, which has been extensively used to study the response of bulk superfluid helium as well as superfluid helium droplets. By including a phenomenological stochastic term in this model, it is possible to obtain thermodynamic equilibrium corresponding to a given temperature by propagating the system in imaginary time. The temperature dependence of thermodynamic quantities, such as the internal energy and entropy, is computed and is compared well with experimental reference data for the bulk liquid up to about 2 K, but begins to gradually deviate above that temperature. Inspection of pseudovorticity during real-time evolution of the system near 2 K reveals the presence of roton quasiparticles, which are suggested to be precursors for quantized vortex rings (Onsager's ghosts), as well as weaker analogs of extended vortex loops.

Numerical Simulation of Finite-Temperature Field Theory for Interacting Bosons

We report a stable and efficient complex Langevin sampling scheme for performing approximation-free numerical simulations directly on the path-integral coherent-states field theory for an assembly of interacting bosons. We apply the method to generate the λ line of critical phase transitions associated with Bose-Einstein condensation in a model ϕ^{4} scalar field theory. The new approach enjoys near-linear scaling in the resolved (d+1) spatial and imaginary-time dimensions and should be particularly efficient for the study of dense systems at low temperature.

Momentum Distribution in the Unitary Bose Gas from First Principles

We consider a realistic bosonic N-particle model with unitary interactions relevant for Efimov physics. Using quantum Monte Carlo methods, we find that the critical temperature for Bose-Einstein condensation is decreased with respect to the ideal Bose gas. We also determine the full momentum distribution of the gas, including its universal asymptotic behavior, and compare this crucial observable to recent experimental data. Similar to the experiments with different atomic species, differentiated solely by a three-body length scale, our model only depends on a single parameter. We establish a weak influence of this parameter on physical observables. In current experiments, the thermodynamic instability of our model from the atomic gas towards an Efimov liquid could be masked by the dynamical instability due to three-body losses.

Critical temperature of interacting Bose gases in periodic potentials

The superfluid transition of a repulsive Bose gas in the presence of a sinusoidal potential which represents a simple-cubic optical lattice is investigated using quantum Monte Carlo simulations. At the average filling of one particle per well the critical temperature has a nonmonotonic dependence on the interaction strength, with an initial sharp increase and a rapid suppression at strong interactions in the vicinity of the Mott transition. In an optical lattice the positive shift of the transition is strongly enhanced compared to the homogenous gas. By varying the lattice filling we find a crossover from a regime where the optical lattice has the dominant effect to a regime where interactions dominate and the presence of the lattice potential becomes almost irrelevant.

Onset temperature of Bose-Einstein condensation in incommensurate solid (4)He

The temperature dependence of the one-body density matrix in (4)He crystals presenting vacancies is computed with path integral Monte Carlo simulations. The main purpose of this study is to estimate the onset temperature T(0) of Bose-Einstein condensation in these systems. We see that T(0) depends on the vacancy concentration X(v) of the simulated system, but not following the law T(0) ~ X(v)(2/3) obtained assuming noninteracting vacancies. For the lowest X(v) we study, that is X(v)= 1/256, we get T(0) = 0.15 ± 0.05 K, close to the temperatures at which a finite fraction of nonclassical rotational inertia is experimentally observed. Below T(0), vacancies do not act as classical point defects becoming completely delocalized entities.

Critical temperature of interacting Bose gases in two and three dimensions

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.

Worm algorithm and diagrammatic Monte Carlo: a new approach to continuous-space path integral Monte Carlo simulations

A detailed description is provided of a new worm algorithm, enabling the accurate computation of thermodynamic properties of quantum many-body systems in continuous space, at finite temperature. The algorithm is formulated within the general path integral Monte Carlo (PIMC) scheme, but also allows one to perform quantum simulations in the grand canonical ensemble, as well as to compute off-diagonal imaginary-time correlation functions, such as the Matsubara Green function, simultaneously with diagonal observables. Another important innovation consists of the expansion of the attractive part of the pairwise potential energy into elementary (diagrammatic) contributions, which are then statistically sampled. This affords a complete microscopic account of the long-range part of the potential energy, while keeping the computational complexity of all updates independent of the size of the simulated system. The computational scheme allows for efficient calculations of the superfluid fraction and off-diagonal correlations in space-time, for system sizes which are orders of magnitude larger than those accessible to conventional PIMC. We present illustrative results for the superfluid transition in bulk liquid 4He in two and three dimensions, as well as the calculation of the chemical potential of hcp 4He.

Electron bubbles in helium clusters. II. Probing superfluidity

In this paper we present calculations of electron tunneling times from the ground electronic state of excess electron bubbles in ((4)He)(N) clusters (N=6500-10(7), cluster radius R=41.5-478 A), where the equilibrium bubble radius varies in the range R(b)=13.5-17.0 A. For the bubble center located at a radial distance d from the cluster surface, the tunneling transition probability was expressed as A(0)phi(d,R)exp(-betad), where beta approximately 1 A(-1) is the exponential parameter, A(0) is the preexponential factor for the bubble located at the cluster center, and phi(d,R) is a correction factor which accounts for cluster curvature effects. Electron tunneling dynamics is grossly affected by the distinct mode of motion of the electron bubble in the image potential within the cluster, which is dissipative (i.e., tau(D)<tau(0)) in normal fluid ((4)He)(N) and ((3)He)(N) clusters, while it is undamped (i.e., tau(D)>>tau(0)) in superfluid ((4)He)(N) clusters, where tau(D) is the bubble motional damping time (tau(D) approximately 4 x 10(-12) s for normal fluid clusters and tau(D) approximately 10 s for superfluid clusters), while tau(0) approximately 10(-9)-10(-10) s is the bubble oscillatory time. Exceedingly long tunneling lifetimes, which cannot be experimentally observed, are manifested from bubbles damped to the center of the normal fluid cluster, while for superfluid clusters electron tunneling occurs from bubbles located in the vicinity of the initial distance d near the cluster boundary. Model calculations of the cluster size dependence of the electron tunneling time (for a fixed value of d=38-39 A), with lifetimes increasing in the range of 10(-3)-0.3 s for N=10(4)-10(7), account well for the experimental data [M. Farnik and J. P. Toennies, J. Chem. Phys. 118, 4176 (2003)], manifesting cluster curvature effects on electron tunneling dynamics. The minimal cluster size for the dynamic stability of the bubble was estimated to be N=3800, which represents the threshold cluster size for which the excess electron bubble in ((4)He)(N) (-) clusters is amenable to experimental observation.

Reentrant phenomenon in the quantum phase transitions of a gas of bosons trapped in an optical lattice

We calculate the location of the quantum phase transitions of a Bose gas trapped in an optical lattice as a function of effective scattering length a(eff) and temperature T. Knowledge of recent high-loop results on the shift of the critical temperature at weak couplings is used to locate a nose in the phase diagram above the free Bose-Einstein critical temperature T((0))(c), thus predicting the existence of a reentrant transition above T((0))(c), where a condensate should form when increasing a(eff). At zero temperature, the transition to the normal phase produces the experimentally observed Mott insulator.

Effect of disorder on the critical temperature of a dilute hard-sphere Gas

We have performed path integral Monte Carlo calculations to determine the effect of quenched disorder on the superfluid density of a dilute 3D hard-sphere gas. The disorder was introduced by locating hard cylinders randomly inside the simulation cell. Our results indicate that the disorder does not strongly affect the superfluid critical temperature. There is a reduction of rho(s)/rho with increasing disorder and with excluded volume for similar disorders and a possible change of universality class (as evidenced by the correlation length exponent) at high disorder. Comparison to experiments of helium in Vycor is made.