Coherent State Path Integral Monte Carlo

Abstract

We propose a new quantum simulation method for a many body quantum liquid of identical particles at finite (non-zero) temperature. The new scheme expands the high temperature density matrix on the overcomplete set of single particles coherent states of John Rider Klauder instead of the usual plane waves as in conventional path integral methods. One is free to tune the elastic constant and/or the mass of the harmonic oscillator subtending the coherent states so as to maximize the computational efficiency of the algorithm. We prove that in the limit of an extremely stiff harmonic oscillator the results for the internal energy tends towards the correct expected values. Moreover we suggest that a stiff harmonic oscillator could allow the use of larger (imaginary) timesteps. This additional degree of freedom is the characteristic feature of our new algorithm and is not available in more conventional path integral methods.

19 pages, 2 tables, 3 figures

References in corpus (12)

• Worm Algorithm for Continuous-space Path Integral Monte Carlo Simulations
• Path Integral Monte Carlo Simulation of the Warm-Dense Homogeneous Electron Gas
• The Uniform Electron Gas at Warm Dense Matter Conditions
• Jellium at finite temperature using the restricted worm algorithm
• Radial distribution function in a diffusion Monte Carlo simulation of a Fermion fluid between the ideal gas and the Jellium model
• Quantum Gibbs ensemble Monte Carlo
• One-component fermion plasma on a sphere at finite temperature
• Two component boson-fermion plasma at finite temperature
• Gas-liquid coexistence for the bosons square-well fluid and the $\mbox{}^4$He binodal anomaly
• Two phase coexistence for the hydrogen-helium mixture
• One-component fermion plasma on a sphere at finite temperature. The anisotropy in the paths conformations
• Supercooled superfluids in Monte Carlo simulations