A path integral Monte Carlo method based on the worm algorithm has been developed to compute the chemical potential of interacting bosonic quantum fluids. By applying it to finite-sized systems of helium-4 atoms, we have confirmed that the chemical p
otential scales inversely with the number of particles to lowest order. The introduction of a simple scaling form allows for the extrapolation of the chemical potential to the thermodynamic limit, where we observe excellent agreement with known experimental results for helium-4 at saturated vapor pressure. We speculate on future applications of the proposed technique, including its use in studies of confined quantum fluids.
We have performed quantum Monte Carlo simulations measuring the finite size and temperature superfluid response of helium-4 to the linear and rotational motion of the walls of a nanopore. Within the two-fluid model, the portion of the normal liquid d
ragged along with the boundaries is dependent on the type of motion and the resulting anisotropic superfluid density saturates far below unity at T=0.5 K. The origin of the saturation is uncovered by computing the spatial distribution of superfluidity, with only the core of the nanopore exhibiting any evidence of phase coherence. The superfluid core displays scaling behavior consistent with Luttinger liquid theory, thereby providing an experimental test for the emergence of a one dimensional quantum liquid.
