Continuum limit and preconditioned Langevin sampling of the path integral molecular dynamics
arXiv:1811.10995
Abstract
We investigate the continuum limit that the number of beads goes to infinity in the ring polymer representation of thermal averages. Studying the continuum limit of the trajectory sampling equation sheds light on possible preconditioning techniques for sampling ring polymer configurations with large number of beads. We propose two preconditioned Langevin sampling dynamics, which are shown to have improved stability and sampling accuracy. We present a careful mode analysis of the preconditioned dynamics and show their connections to the normal mode, the staging coordinate and the Matsubara mode representation for ring polymers. In the case where the potential is quadratic, we show that the continuum limit of the preconditioned mass modified Langevin dynamics converges to its equilibrium exponentially fast, which suggests that the finite-dimensional counterpart has a dimension-independent convergence rate. In addition, the preconditioning techniques can be naturally applied to the multi-level quantum systems in the nonadiabatic regime, which are compatible with various numerical approaches.
References in corpus (8)
- Efficient stochastic thermostatting of path integral molecular dynamics
- Nuclear quantum effects in solids using a colored-noise thermostat
- Boltzmann-conserving classical dynamics in quantum time-correlation functions: Matsubara dynamics
- Path-Integral Isomorphic Hamiltonian for Including Nuclear Quantum Effects in Non-adiabatic Dynamics
- The motion of a random string
- Path integral molecular dynamics with surface hopping for thermal equilibrium sampling of nonadiabatic systems
- Path Integral Molecular Dynamics for Exact Quantum Statistics of Multi-Electronic-State Systems
- Adaptive Dimension Reduction to Accelerate Infinite-Dimensional Geometric Markov Chain Monte Carlo