Quantum statistical calculations and symplectic corrector algorithms
- 30 April 2004
- journal article
- research article
- Published by American Physical Society (APS) in Physical Review E
- Vol. 69 (4), 046118
- https://doi.org/10.1103/physreve.69.046118
Abstract
The quantum partition function at finite temperature requires computing the trace of the imaginary time propagator. For numerical and Monte Carlo calculations, the propagator is usually split into its kinetic and potential parts. A higher-order splitting will result in a higher-order convergent algorithm. At imaginary time, the kinetic energy propagator is usually the diffusion Green’s function. Since diffusion cannot be simulated backward in time, the splitting must maintain the positivity of all intermediate time steps. However, since the trace is invariant under similarity transformations of the propagator, one can use this freedom to “correct” the split propagator to higher order. This use of similarity transforms classically gives rise to symplectic corrector algorithms. The split propagator is the symplectic kernel and the similarity transformation is the corrector. This work proves a generalization of the Sheng-Suzuki theorem: no positive time step propagators with only kinetic and potential operators can be corrected beyond second order. Second-order forward propagators can have fourth-order traces only with the inclusion of an additional commutator. We give detailed derivations of four forward correctable second-order propagators and their minimal correctors.Keywords
All Related Versions
This publication has 23 references indexed in Scilit:
- Short-time-evolved wave functions for solving quantum many-body problemsPhysical Review B, 2003
- Symplectic analytically integrable decomposition algorithms: classification, derivation, and application to molecular dynamics, quantum and celestial mechanics simulationsComputer Physics Communications, 2003
- Splitting methodsActa Numerica, 2002
- A fourth-order real-space algorithm for solving local Schrödinger equationsThe Journal of Chemical Physics, 2001
- High order numerical integrators for differential equations using composition and processing of low order methodsApplied Numerical Mathematics, 2001
- Fourth-order diffusion Monte Carlo algorithms for solving quantum many-body problemsPhysical Review B, 2001
- Fourth-order algorithms for solving the multivariable Langevin equation and the Kramers equationPhysical Review E, 2000
- Hybrid exponential product formulas for unbounded operators with possible applications to Monte Carlo simulationsPhysics Letters A, 1995
- Quadratic diffusion Monte Carlo algorithms for solving atomic many-body problemsPhysical Review A, 1990
- Helium at zero temperature with hard-sphere and other forcesPhysical Review A, 1974