Quantum Statistical Calculations and Symplectic Corrector Algorithms

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 Greens 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 give rises 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.

9 pages, no figure, corrected typos, mostly missing right brackets

References in corpus (3)

• Gradient Symplectic Algorithms for Solving the Schroedinger Equation with Time-Dependent Potentials
• On the construction of high-order force gradient algorithms for integration of motion in classical and quantum systems
• Short time evolved wave functions for solving quantum many-body problems

Cited by in corpus (18)

• High order Chin actions in path integral Monte Carlo
• 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
• High-order Time Expansion Path Integral Ground State
• Forward Symplectic Integrators and the Long Time Phase Error in Periodic Motions
• Berezinskii-Kosterlitz-Thouless Transition in Two-Dimensional Dipolar Stripes
• The physics of symplectic integrators: perihelion advances and symplectic corrector algorithms
• Monte Carlo evaluation of the equilibrium isotope effects using the Takahashi-Imada factorization of the Feynman path integral
• Phase diagram of Rydberg atoms with repulsive van der Waals interaction
• Accelerating quantum instanton calculations of the kinetic isotope effects
• A Fundamental Theorem on the Structure of Symplectic Integrators
• The Complete Characterization of Fourth-Order Symplectic Integrators with Extended-Linear Coefficients
• A quantum Monte Carlo based density functional for Dysprosium dipolar system
• Anatomy of Path Integral Monte Carlo: algebraic derivation of the harmonic oscillator's universal discrete imaginary-time propagator and its sequential optimization
• Reply to the Comment on "Berezinskii-Kosterlitz-Thouless Transition in Two-Dimensional Dipolar Stripes"
• Negative time splitting is stable
• Fokker-Planck equation for coupled anharmonic oscillators: symplectic integration and application to mean-field model of cooperative behavior
• Forward and non-forward symplectic integrators in solving classical dynamics problems