Dimension-free path-integral molecular dynamics without preconditioning

Abstract

Convergence with respect to imaginary-time discretization is an essential part of any path-integral-based calculation. However, an unfortunate property of existing non-preconditioned numerical integration schemes for path-integral molecular dynamics (PIMD) - including ring-polymer molecular dynamics (RPMD) and thermostatted RPMD (T-RPMD) - is that for a given MD timestep, the overlap between the exact ring-polymer Boltzmann-Gibbs distribution and that sampled using MD becomes zero in the infinite-bead limit. This has clear implications for hybrid Metropolis Monte-Carlo/MD sampling schemes. We show that these problems can be avoided through the introduction of "dimension-free" numerical integration schemes for which the sampled ring-polymer position distribution has non-zero overlap with the exact distribution in the infinite-bead limit for the case of a harmonic potential. We show that dimension freedom can be achieved via mollification of the forces from the physical potential and with the BCOCB integration scheme. The dimension-free numerical integration schemes yield finite error bounds for a given MD timestep as the number of beads is taken to infinity; these conclusions are proven for harmonic potential and borne out numerically for anharmonic systems, including water. The numerical results for BCOCB are particularly striking, allowing for three-fold increases in the stable timestep for liquid water with respect to the Bussi-Parrinello (OBABO) and Leimkuhler (BAOAB) integrators while introducing negligible errors in the statistical properties and absorption spectrum. Importantly, the dimension-free, non-preconditioned integration schemes introduced here preserve ergodicity and global second-order accuracy, and they remain simple, black-box methods that avoid additional computational costs, tunable parameters, or system-specific implementations.

arXiv comment13 pages, 6 figures. New results (for the 2 fs time-step) are added in Figure 1(d), y-axes of figure 6 are rescaled and a few typos are corrected

References in corpus (15)

• Canonical sampling through velocity-rescaling
• Competing quantum effects in the dynamics of a flexible water model
• Efficient stochastic thermostatting of path integral molecular dynamics
• How to remove the spurious resonances from ring polymer molecular dynamics
• Accelerating the convergence of path integral dynamics with a generalized Langevin equation
• Rational Construction of Stochastic Numerical Methods for Molecular Sampling
• A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat
• Long-Run Accuracy of Variational Integrators in the Stochastic Context
• A unified thermostat scheme for efficient configurational sampling for classical/quantum canonical ensembles via molecular dynamics
• Geometric integrators and the Hamiltonian Monte Carlo method
• Fine Tuning Classical and Quantum Molecular Dynamics using a Generalized Langevin Equation
• Cayley modification for strongly stable path-integral and ring-polymer molecular dynamics
• Pathwise Accuracy and Ergodicity of Metropolized Integrators for SDEs
• Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics
• Modified trigonometric integrators

Cited by in corpus (14)

• Use the force! Reduced variance estimators for densities, radial distribution functions and local mobilities in molecular simulations
• A generalized class of strongly stable and dimension-free T-RPMD integrators
• Routine Molecular Dynamics Simulations Including Nuclear Quantum Effects: from Force Fields to Machine Learning Potentials
• Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches
• Convergence of unadjusted Hamiltonian Monte Carlo for mean-field models
• Confirming the role of nuclear tunnelling in aqueous ferrous-ferric electron transfer
• Effective Mass Path Integral Simulations of Quasiparticles in Condensed Phases
• Two-scale coupling for preconditioned Hamiltonian Monte Carlo in infinite dimensions
• Nuclear quantum effects in thermal conductivity from centroid molecular dynamics
• New Stable and Fast Ring-Polymer Molecular Dynamics for Calculating Bimolecular Rate Coefficients with Example of OH + CH_4
• Efficient Sampling of Thermal Averages of Interacting Quantum Particle Systems with Random Batches
• Harmonic Oscillator Staging Coordinates for Efficient Path Integral Simulations of Quantum Oscillators and Crystals
• NQCDynamics.jl: A Julia Package for Nonadiabatic Quantum Classical Molecular Dynamics in the Condensed Phase
• Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy