Higher order and infinite Trotter-number extrapolations in path integral Monte Carlo

Abstract

Improvements beyond the primitive approximation in the path integral Monte Carlo method are explored both in a model problem and in real systems. Two different strategies are studied: the Richardson extrapolation on top of the path integral Monte Carlo data and the Takahashi-Imada action. The Richardson extrapolation, mainly combined with the primitive action, always reduces the number-of-beads dependence, helps in determining the approach to the dominant power law behavior, and all without additional computational cost. The Takahashi-Imada action has been tested in two hard-core interacting quantum liquids at low temperature. The results obtained show that the fourth-order behavior near the asymptote is conserved, and that the use of this improved action reduces the computing time with respect to the primitive approximation.

19 pages, RevTex, to appear in J. Chem. Phys

References in corpus (4)

• Gradient Symplectic Algorithms for Solving the Schroedinger Equation with Time-Dependent Potentials
• An ab initio path integral Monte Carlo simulation method for molecules and clusters: application to Li_4 and Li_5^+
• Optimal Energy Estimation in Path-Integral Monte Carlo Simulations
• Comparison of two non-primitive methods for path integral simulations: Higher-order corrections vs. an effective propagator approach

Cited by in corpus (26)

• The Fermion Sign Problem in Path Integral Monte Carlo Simulations: Quantum Dots, Ultracold Atoms, and Warm Dense Matter
• Electronic Density Response of Warm Dense Matter
• High order Chin actions in path integral Monte Carlo
• Nonlinear Electronic Density Response in Warm Dense Matter
• Path Integral Ground State with a Fourth-Order Propagator: Application to Condensed Helium
• Mitigating algorithmic errors in Hamiltonian simulation
• Path integral virial estimator based on the scaling of fluctuation coordinates: Application to quantum clusters with fourth-order propagators
• Dynamic properties of the warm dense electron gas: an ab initio path integral Monte Carlo approach
• Nonlinear density response from imaginary-time correlation functions: Ab initio path integral Monte Carlo simulations of the warm dense electron gas
• Attenuating the fermion sign problem in path integral Monte Carlo simulations using the Bogoliubov inequality and thermodynamic integration
• Finite-size effects in the reconstruction of dynamic properties from ab initio path integral Monte-Carlo simulations
• Electronic Density Response of Warm Dense Hydrogen: Ab initio Path Integral Monte Carlo Simulations
• Any order imaginary time propagation method for solving the Schrodinger equation
• The Relevance of Electronic Perturbations in the Warm Dense Electron Gas
• Benchmarking Exchange-Correlation Functionals in the Spin-Polarized Inhomogeneous Electron Gas under Warm Dense Conditions
• Effective electronic forces and potentials from ab initio path integral Monte Carlo simulations
• Ab initio path integral Monte Carlo simulation of the Uniform Electron Gas in the High Energy Density Regime
• Equation of state of an interacting Bose gas at finite temperature: a Path Integral Monte Carlo study
• Lattice stability and high pressure melting mechanism of dense hydrogen up to 1.5 TPa
• Monte Carlo evaluation of the equilibrium isotope effects using the Takahashi-Imada factorization of the Feynman path integral
• Monte Carlo studies of quantum and classical annealing on a double-well
• Fluctuation-induced interactions between dielectrics in general geometries
• Accelerating ab initio path integral molecular dynamics with multilevel sampling of potential surface
• On the dynamics of Simulated Quantum Annealing in random Ising chains
• Abnormal superfluid fraction and structural properties of electrons in 2D and 3D quantum dots: an ab initio path-integral Monte Carlo study
• Energy response and spatial alignment of the perturbed electron gas