Comparison of two non-primitive methods for path integral simulations: Higher-order corrections vs. an effective propagator approach

Abstract

Two methods are compared that are used in path integral simulations. Both methods aim to achieve faster convergence to the quantum limit than the so-called primitive algorithm (PA). One method, originally proposed by Takahashi and Imada, is based on a higher-order approximation (HOA) of the quantum mechanical density operator. The other method is based upon an effective propagator (EPr). This propagator is constructed such that it produces correctly one and two-particle imaginary time correlation functions in the limit of small densities even for finite Trotter numbers P. We discuss the conceptual differences between both methods and compare the convergence rate of both approaches. While the HOA method converges faster than the EPr approach, EPr gives surprisingly good estimates of thermal quantities already for P = 1. Despite a significant improvement with respect to PA, neither HOA nor EPr overcomes the need to increase P linearly with inverse temperature. We also derive the proper estimator for radial distribution functions for HOA based path integral simulations.

17 pages, latex, 6 postscript figures

Cited by in corpus (5)

• Higher order and infinite Trotter-number extrapolations in path integral Monte Carlo
• Path integral virial estimator based on the scaling of fluctuation coordinates: Application to quantum clusters with fourth-order propagators
• Systematically Accelerated Convergence of Path Integrals
• Systematic Speedup of Path Integrals of a Generic $N$-fold Discretized Theory
• Quantum creep and quantum creep transitions in 1D sine-Gordan chains