Worm Algorithm and Diagrammatic Monte Carlo: A New Approach to Continuous-Space Path Integral Monte Carlo Simulations
arXiv:physics/0605225 · doi:10.1103/PhysRevE.74.036701
Abstract
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 helium-four in two and three dimensions, as well as the calculation of the chemical potential of solid helium-four.
17 pages, 13 figures
References in corpus (3)
Cited by in corpus (31)
- A superfluid-droplet crystal and a free-space supersolid in a dipole-blockaded gas
- Luttinger Liquid in the Core of Screw Dislocation in Helium-4
- Structure, superfluidity, and quantum melting of hydrogen clusters
- Helium-4 Luttinger liquids in nanopores
- Phase diagram of 4He adsorbed on graphite
- A molecular superfluid: non-classical rotations in doped para-hydrogen clusters
- Local Superfluidity of Parahydrogen Clusters
- Interacting bosons in one dimension and Luttinger liquid theory
- Exact ground state Monte Carlo method for Bosons without importance sampling
- Bose-Einstein-condensed gases in arbitrarily strong random potentials
- Momentum distribution of 1D Bose gases at the quasi-condensation crossover: theoretical and experimental investigation
- Supersolid phases of cold atom assemblies
- Collective and single-particle excitations in 2D dipolar Bose gases
- Role of Bose Statistics in Crystallization and Quantum Jamming
- A path integral Monte Carlo method for Rényi entanglement entropies
- Bose-Einstein-condensed gases with arbitrary strong interactions
- Quantum statistics and the momentum distribution of liquid para-hydrogen
- A semi-classical field method for the equilibrium Bose gas and application to thermal vortices in two dimensions
- Binding of a 3He impurity to a screw dislocation in solid 4He
- Superfluidity of strongly correlated bosons in two- and three-dimensional traps
- Superfluidity of isotopically doped parahydrogen clusters
- Bose-Einstein condensation in self-consistent mean-field theory
- Path Integral Monte Carlo study confirms a highly ordered snowball in 4He nanodroplets doped with an Ar+ ion
- Disorder and the elusive superfluid phase of para-hydrogen
- Observation of crystallization slowdown in supercooled para-hydrogen and ortho-deuterium quantum liquid mixtures
- Molecular hydrogen isotopes adsorbed on krypton-preplated graphite: Quantum Monte Carlo simulations
- Enhanced superfluid response of parahydrogen in nanoscale confinement
- Quantum Gibbs ensemble Monte Carlo
- Gas-liquid coexistence for the bosons square-well fluid and the $\mbox{}^4$He binodal anomaly
- Evaluation of the Permutational Structure of the Canonical Ensemble Partition Function
- Superfluid response of two-dimensional parahydrogen clusters in confinement