Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics

High-accuracy Trotter-formula method for path integrals

Kevin E. Schmidt, Michael A. Lee

Published Jun 1, 1995• 38 citations
Publisher

Abstract

Path integrals are a powerful method for calculating real time, finite temperature, and ground state properties of quantum systems. By exploiting some remarkable properties of the symmetric Trotter formula and the discrete Fourier transform, we arrive at a high-accuracy method for removing ``time slice'' errors in Trotter-approximated propagators. We provide an explicit demonstration of the method applied to the two-body density matrix of $^{4}\mathrm{He}$. Our method is simultaneously fast, high precision, and computationally simple and can be applied to a wide variety of quantum propagators.

DOI
10.1103/physreve.51.5495

Related papers

Improved Feynman’s path integral method with a large time step: Formalism and applications

Alexander N. Drozdov

We describe an efficient path integral scheme for calculating the propagator of an arbitrary quantum system, as well as that of a stochastic system in special cases where the Fokker–Planck equation obeys strict detailed balance. The basic idea is to split the respective Hamiltonian into two exactly solvable parts and then to employ a symmetric decomposition of the time evolution operator, which is exact up to a high order in the time step. The resulting single step propagator allows rather large time steps in a path integral and leads to convergence with fewer time slices. Because it involves no system-specific reference system, the algorithm is amenable to all known numerical schemes available for evaluating quantum path integrals. In this way one obtains a highly accurate method, which is simultaneously fast, stable, and computationally simple. Numerical applications to the real time quantum dynamics in a double well and to the stochastic dynamics of a bistable system coupled to a harmonic mode show our method to be superior over the approach developed by the Makri group in their quasiadiabatic propagator representation, to say nothing about the propagation scheme based on the standard Trotter splitting.

The Journal of Chemical Physics•Apr 22, 1998•14 citations

Discretized path integral method and properties of a quantum system

Ashok Sethia, Sudip Sanyal, Yashwant Singh

Using the discretized path integral formalism we develop a numerically convenient and accurate method for evaluating the finite time propagator (and density matrix) of a given system. The eigenfunctions and energies of large number of states of the system are also found. We demonstrate the accuracy of the method by computing eigenvalues, eigenfunctions, position and flux autocorrelation functions and complex time propagators for a large number of systems and comparing these results with those available from other sources. In all the cases we find very good agreement. The method thus provides a very simple and accurate procedure for the study of static and dynamic properties of a quantum system.

The Journal of Chemical Physics•Nov 15, 1990•40 citations

The Feynman path integral

Mark E. Tuckerman

Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations

•Aug 1, 2023•0 citations

Algebraic Generalization

Tarō Kashiwa, Yoshio Ohnuki, Masuo Suzuki

Abstract The path integral method discussed up to the preceding section is based on Trotter’s formula (1.10) and its accuracy is of the order of Δt for the time slice Δt. This is sufficient for discussing quantum systems qualitatively. Higher accuracy of the Quantum Monte method is required for studying quantum systems numerically. In the present chapter, we explain the algebraic generalization of the path integral method. By applying the integral representations given in chapters 1−4 to such generalized product formulae, some new kinds of path integral of higher accuracy can be obtained. This is crucial when the essence of a phenomenon can be found only by increasing the accuracy of calculations. It should be emphasized that such higher order decomposition formulae to be explained in the present chapter can also be applied not only to the path integral method but to many other problems such as non-linear dynamics in Hamiltonian systems.

•Jul 3, 1997•0 citations

Path integral density matrix dynamics: A method for calculating time-dependent properties in thermal adiabatic and non-adiabatic systems

Scott Habershon

We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency.

The Journal of Chemical Physics•Sep 12, 2013•8 citations

On the approximation of Feynman-Kac path integrals for quantum statistical mechanics

Stephen Bond, Brian B. Laird, Benedict Leimkuhler

Discretizations of the Feynman-Kac path integral representation of the quantum mechanical density matrix are investigated. Each infinite-dimensional path integral is approximated by a Riemann integral over a finite-dimensional function space, by restricting the integration to a subspace of all admissible paths. Using this process, a wide class of methods can be derived, with each method corresponding to a different choice for the approximating subspace. The traditional ``short-time'' approximation and ``Fourier discretization'' can be recovered from this approach, using linear and spectral basis functions respectively. As an illustration, a novel method is formulated using cubic elements and is shown to have improved convergence properties when applied to a simple model problem.

arXiv:cond-mat/0007112•Jul 6, 2000•0 citations
Path integrals are a powerful method for calculating real time, finite temperature, and ground state properties of quantum systems. By exploiting some remarkable properties of the symmetric Trotter formula and the discrete Fourier transform, we arrive at a high-accuracy method for removing ``time slice'' errors in Trotter-approximated propagators. We provide an explicit demonstration of the method applied to the two-body density matrix of 4He^{4}\mathrm{He}4He. Our method is simultaneously fast, high precision, and computationally simple and can be applied to a wide variety of quantum propagators.
We describe an efficient path integral scheme for calculating the propagator of an arbitrary quantum system, as well as that of a stochastic system in special cases where the Fokker–Planck equation obeys strict detailed balance. The basic idea is to split the respective Hamiltonian into two exactly solvable parts and then to employ a symmetric decomposition of the time evolution operator, which is exact up to a high order in the time step. The resulting single step propagator allows rather large time steps in a path integral and leads to convergence with fewer time slices. Because it involves no system-specific reference system, the algorithm is amenable to all known numerical schemes available for evaluating quantum path integrals. In this way one obtains a highly accurate method, which is simultaneously fast, stable, and computationally simple. Numerical applications to the real time quantum dynamics in a double well and to the stochastic dynamics of a bistable system coupled to a harmonic mode show our method to be superior over the approach developed by the Makri group in their quasiadiabatic propagator representation, to say nothing about the propagation scheme based on the standard Trotter splitting.
Using the discretized path integral formalism we develop a numerically convenient and accurate method for evaluating the finite time propagator (and density matrix) of a given system. The eigenfunctions and energies of large number of states of the system are also found. We demonstrate the accuracy of the method by computing eigenvalues, eigenfunctions, position and flux autocorrelation functions and complex time propagators for a large number of systems and comparing these results with those available from other sources. In all the cases we find very good agreement. The method thus provides a very simple and accurate procedure for the study of static and dynamic properties of a quantum system.
Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations
Abstract The path integral method discussed up to the preceding section is based on Trotter’s formula (1.10) and its accuracy is of the order of Δt for the time slice Δt. This is sufficient for discussing quantum systems qualitatively. Higher accuracy of the Quantum Monte method is required for studying quantum systems numerically. In the present chapter, we explain the algebraic generalization of the path integral method. By applying the integral representations given in chapters 1−4 to such generalized product formulae, some new kinds of path integral of higher accuracy can be obtained. This is crucial when the essence of a phenomenon can be found only by increasing the accuracy of calculations. It should be emphasized that such higher order decomposition formulae to be explained in the present chapter can also be applied not only to the path integral method but to many other problems such as non-linear dynamics in Hamiltonian systems.
We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency.
Discretizations of the Feynman-Kac path integral representation of the quantum mechanical density matrix are investigated. Each infinite-dimensional path integral is approximated by a Riemann integral over a finite-dimensional function space, by restricting the integration to a subspace of all admissible paths. Using this process, a wide class of methods can be derived, with each method corresponding to a different choice for the approximating subspace. The traditional ``short-time'' approximation and ``Fourier discretization'' can be recovered from this approach, using linear and spectral basis functions respectively. As an illustration, a novel method is formulated using cubic elements and is shown to have improved convergence properties when applied to a simple model problem.