Permutation sampling in Path Integral Monte Carlo

A simple algorithm is described to sample permutations of identical particles in Path Integral Monte Carlo (PIMC) simulations of continuum many-body systems. The sampling strategy illustrated here is fairly general, and can be easily incorporated in any PIMC implementation based on the staging algorithm. Although it is similar in spirit to an existing prescription, it differs from it in some key aspects. It allows one to sample permutations efficiently, even if long paths (e.g., hundreds, or thousands of slices) are needed. We illustrate its effectiveness by presenting results of a PIMC calculation of thermodynamic properties of superfluid He-4, in which a very simple approximation for the high-temperature density matrix was utilized.