Upon the existence of short-time approximations of any polynomial order for the computation of density matrices by path integral methods

Abstract

In this article, I provide significant mathematical evidence in support of the existence of short-time approximations of any polynomial order for the computation of density matrices of physical systems described by arbitrarily smooth and bounded from below potentials. While for Theorem 2, which is ``experimental'', I only provide a ``physicist's'' proof, I believe the present development is mathematically sound. As a verification, I explicitly construct two short-time approximations to the density matrix having convergence orders 3 and 4, respectively. Furthermore, in the Appendix, I derive the convergence constant for the trapezoidal Trotter path integral technique. The convergence orders and constants are then verified by numerical simulations. While the two short-time approximations constructed are of sure interest to physicists and chemists involved in Monte Carlo path integral simulations, the present article is also aimed at the mathematical community, who might find the results interesting and worth exploring. I conclude the paper by discussing the implications of the present findings with respect to the solvability of the dynamical sign problem appearing in real-time Feynman path integral simulations.

19 pages, 4 figures; the discrete short-time approximations are now treated as independent from their continuous version; new examples of discrete short-time approximations of order three and four are given; a new appendix containing a short review on Brownian motion has been added; also, some additional explanations are provided here and there; this is the last version; to appear in Phys. Rev. E

References in corpus (6)

• Heat capacity estimators for random series path-integral methods by finite-difference schemes
• Optimal series representations for numerical path integral simulations
• Numerical implementation of some reweighted path integral methods
• Energy estimators for random series path-integral methods
• Random Series and Discrete Path Integral methods: The Levy-Ciesielski implementation
• The Partial Averaging method

Cited by in corpus (12)

• Thermodynamics and equilibrium structure of Ne_38 cluster: Quantum Mechanics versus Classical
• Extrapolated High-Order Propagators for Path Integral Monte Carlo Simulations
• Path integral virial estimator based on the scaling of fluctuation coordinates: Application to quantum clusters with fourth-order propagators
• Heat capacity estimators for random series path-integral methods by finite-difference schemes
• Sampling diffusive transition paths
• Recursive Schr\" odinger Equation Approach to Faster Converging Path Integrals
• Fast Convergence of Path Integrals for Many-body Systems
• Onset temperature of Bose-Einstein condensation in incommensurate solid 4He
• Reconstruction of thermally-symmetrized quantum autocorrelation functions from imaginary-time data
• The fast sampling algorithm for Lie-Trotter products
• On the efficient Monte Carlo implementation of path integrals
• Moments of spectral functions: Monte Carlo evaluation and verification