Path-integral Monte Carlo worm algorithm for Bose systems with periodic boundary conditions

Abstract

We provide a detailed description of the path-integral Monte Carlo worm algorithm used to exactly calculate the thermodynamics of Bose systems in the canonical ensemble. The algorithm is fully consistent with periodic boundary conditions, that are applied to simulate homogeneous phases of bulk systems, and it does not require any limitation in the length of the Monte Carlo moves realizing the sampling of the probability distribution function in the space of path configurations. The result is achieved adopting a representation of the path coordinates where only the initial point of each path is inside the simulation box, the remaining ones being free to span the entire space. Detailed balance can thereby be ensured for any update of the path configurations without the ambiguity of the selection of the periodic image of the particles involved. We benchmark the algorithm using the non-interacting Bose gas model for which exact results for the partition function at finite number of particles can be derived. Convergence issues and the approach to the thermodynamic limit are also addressed for interacting systems of hard spheres in the regime of high density.

arXiv commentv1: 18 pages, 6 figures. v2: Fixed typo in eq.(30) and (31), minor changes, matches published version

References in corpus (21)

• Worm Algorithm and Diagrammatic Monte Carlo: A New Approach to Continuous-Space Path Integral Monte Carlo Simulations
• Worm Algorithm for Continuous-space Path Integral Monte Carlo Simulations
• Exact, Complete, and Universal Continuous-Time Worldline Monte Carlo Approach to the Statistics of Discrete Quantum Systems
• Worm algorithms for classical statistical models
• A superfluid-droplet crystal and a free-space supersolid in a dipole-blockaded gas
• Superglass Phase of Helium-four
• The fate of vacancy-induced supersolidity in 4He
• Defect-induced supersolidity with soft-core bosons
• Phase diagram of 4He adsorbed on graphite
• High order Chin actions in path integral Monte Carlo
• Phase diagram of soft-core bosons in two dimensions
• Permutation sampling in Path Integral Monte Carlo
• Superfluid filaments of dipolar bosons in free space
• Dynamic structure factor of liquid 4He across the normal-superfluid transition
• Universal Superfluid Transition and Transport Properties of Two-Dimensional Dirty Bosons
• Quantum statistics and the momentum distribution of liquid para-hydrogen
• Phases of dipolar bosons in a bilayer geometry
• Superfluidity of metastable bulk glass para-hydrogen at low temperature
• Berezinskii-Kosterlitz-Thouless Transition in Two-Dimensional Dipolar Stripes
• Path Integral Monte Carlo calculation of momentum distribution in solid \^4\He
• Thermodynamics of a dilute Bose gas: A path-integral Monte Carlo study

Cited by in corpus (11)

• Ab initio Calculation of Fluid Properties for Precision Metrology
• Normal liquid ${}^3$He studied by Path Integral Monte Carlo with a parametrized partition function
• Attractive Solution of Binary Bose Mixtures: Liquid-Vapor Coexistence and Critical Point
• Quadratic scaling path integral molecular dynamics for fictitious identical particles and its application to fermion systems
• Phase separation in binary Bose mixtures at finite temperature
• Quantum droplets in two-dimensional Bose mixtures at finite temperature
• Revisiting the properties of superfluid and normal liquid ${}^4$He using ab initio potentials
• Periodic Boundary Conditions for Bosonic Path Integral Molecular Dynamics
• Supersolid phase in two-dimensional soft-core bosons at finite temperature
• A Pseudo-Fermion Propagator Approach to the Fermion Sign Problem
• Combining Harmonic Sampling with the Worm Algorithm to Improve the Efficiency of Path Integral Monte Carlo