A Monte Carlo method for quantum Boltzmann statistical mechanics using Fourier representations of path integrals

Abstract
By expanding Feynman path integrals in a Fourier series a practical Monte Carlo method is developed to calculate the thermodynamic properties of interacting systems obeying quantum Boltzmann statistical mechanics. Working expressions are developed to calculate internal energies, heat capacities, and quantum corrections to free energies. The method is applied to the harmonic oscillator, a double-well potential, and clusters of Lennard-Jones atoms parametrized to mimic the behavior of argon. The expansion of the path integrals in a Fourier series is found to be rapidly convergent and the computational effort for quantum calculations is found to be within an order of magnitude of the corresponding classical calculations. Unlike other related methods no special techniques are required to handle systems with strong short-range repulsive forces.