Hybrid Monte Carlo implementation of the Fourier path integral algorithm

Abstract

This paper formulates a hybrid Monte Carlo implementation of the Fourier path integral (FPI-HMC) approach with partial averaging. Such a hybrid Monte Carlo approach allows one to generate collective moves through configuration space using molecular dynamics while retaining the computational advantages associated with the Fourier path integral Monte Carlo method. In comparison with the earlier Metropolis Monte Carlo implementations of the FPI algorithm, the present HMC method is shown to be significantly more efficient for quantum Lennard-Jones solids and suggests that such algorithms may prove useful for efficient simulations of a range of atomic and molecular systems.