How to obtain thermostatted ring polymer molecular dynamics from exact quantum dynamics and when to use it

Abstract

We obtain thermostatted ring polymer molecular dynamics (TRPMD) from exact quantum dynamics via Matsubara dynamics, a recently-derived form of linearization which conserves the quantum Boltzmann distribution. Performing a contour integral in the complex quantum Boltzmann distribution of Matsubara dynamics, replacement of the imaginary Liouvillian which results with a Fokker-Planck term gives TRPMD. We thereby provide error terms between TRPMD and quantum dynamics and predict the systems in which they are likely to be small. Using a harmonic analysis we show that careful addition of friction causes the correct oscillation frequency of the higher ring-polymer normal modes in a harmonic well, which we illustrate with calculation of the position-squared autocorrelation function. However, no physical friction parameter will produce the correct fluctuation dynamics for a parabolic barrier. The results in this paper are consistent with previous numerical studies and advise the use of TRPMD for the computation of spectra.

12 pages, 2 figures. Submitted to Molecular Physics

References in corpus (10)

• Accurate sampling using Langevin dynamics
• Efficient stochastic thermostatting of path integral molecular dynamics
• How to remove the spurious resonances from ring polymer molecular dynamics
• Colored-noise thermostats à la carte
• Boltzmann-conserving classical dynamics in quantum time-correlation functions: Matsubara dynamics
• Relation of centroid molecular dynamics and ring-polymer molecular dynamics to exact quantum dynamics
• On the Consistency of Approximate Quantum Dynamics Simulation Methods for Vibrational Spectra in the Condensed Phase
• Exact quantum statistics for electronically nonadiabatic systems using continuous path variables
• Non-oscillatory flux correlation functions for efficient nonadiabatic rate theory
• Should Thermostatted Ring Polymer Molecular Dynamics be used to calculate thermal reaction rates?

Cited by in corpus (11)

• Fine Tuning Classical and Quantum Molecular Dynamics using a Generalized Langevin Equation
• Competing quantum effects in the free energy profiles and diffusion rates of hydrogen and deuterium molecules through clathrate hydrates
• Deriving the exact nonadiabatic quantum propagator in the mapping variable representation
• Progress and Challenges in Ab Initio Simulations of Quantum Nuclei in Weakly Bonded Systems
• State Dependent Ring Polymer Molecular Dynamics for Investigating Excited Nonadiabatic Dynamics
• Thermal quantum time-correlation functions from classical-like dynamics
• Probing anharmonic phonons by quantum correlators: A path integral approach
• Simulating vibronic spectra via Matsubara dynamics: coping with the sign problem
• Which Algorithm Best Propagates the Meyer-Miller-Stock-Thoss Mapping Hamiltonian for Non-Adiabatic Dynamics?
• 2D spectroscopies from condensed phase dynamics: Accessing third-order response properties from equilibrium multi-time correlation functions
• On the electronic path integral normal modes of the Meyer-Miller-Stock-Thoss representation of nonadiabatic dynamics