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 ${}^4\mathrm{He}$ in two and three dimensions, as well as the calculation of the chemical potential of hcp ${}^4\mathrm{He}$.

The ongoing experimental and theoretical effort aimed at understanding nonclassical rotational inertia in solid helium has sparked renewed interest in the supersolid phase of matter, its microscopic origin and character, and its experimental detection. The purpose of this Colloquium is to provide a general theoretical framework for the phenomenon of supersolidity, review some of the experimental evidence for solid He4, and discuss its possible interpretation in terms of physical effects underlain by extended defects (such as dislocations). Quantitative support to our theoretical scenarios by means of first-principle numerical simulations is provided. Alternate avenues for the observation of the supersolid phase, not involving helium but rather assemblies of ultracold atoms, are also discussed.

We present a new approach to path integral Monte Carlo (PIMC) simulations based on the worm algorithm, originally developed for lattice models and extended here to continuous-space many-body systems. The scheme allows for efficient computation of thermodynamic properties, including winding numbers and off-diagonal correlations, for systems of much greater size than that accessible to conventional PIMC simulations. As an illustrative application of the method, we simulate the superfluid transition of ${}^4\mathrm{He}$ in two dimensions.

We discuss techniques to generate long-range interactions in a gas of ground state alkali atoms, by weakly admixing excited Rydberg states with laser light. This provides a tool to engineer strongly correlated phases with reduced decoherence from inelastic collisions and spontaneous emission. As an illustration, we discuss the quantum phases of dressed atoms with dipole-dipole interactions confined in a harmonic potential, as relevant to experiments. We show that residual spontaneous emission from the Rydberg state acts as a heating mechanism, leading to a quantum-classical crossover.

A novel supersolid phase is predicted for an ensemble of Rydberg atoms in the dipole-blockade regime, interacting via a repulsive dipolar potential softened at short distances. Using exact numerical techniques, we study the low-temperature phase diagram of this system, and observe an intriguing phase consisting of a crystal of mesoscopic superfluid droplets. At low temperature, phase coherence throughout the whole system, and the ensuing bulk superfluidity, are established through tunnelling of identical particles between neighboring droplets.

We study properties of the supersolid phase observed for hard-core bosons on the triangular lattice near half-integer filling factor, and the phase diagram of the system at finite temperature. We find that the solid order is always of the $(2m,-m^{\prime },-m^{\prime })$ with $m$ changing discontinuously from positive to negative values at half filling, in contrast with phases observed for Ising spins in a transverse magnetic field. At finite temperature we find two intersecting second-order transition lines: one in the 3-state Potts universality class and the other of the Kosterlitz-Thouless type.