How to remove the spurious resonances from ring polymer molecular dynamics

Abstract

Two of the most successful methods that are presently available for simulating the quantum dynamics of condensed phase systems are centroid molecular dynamics (CMD) and ring polymer molecular dynamics (RPMD). Despite their conceptual differences, practical implementations of these methods differ in just two respects: the choice of the Parrinello-Rahman mass matrix and whether or not a thermostat is applied to the internal modes of the ring polymer during the dynamics. Here we explore a method which is halfway between the two approximations: we keep the path integral bead masses equal to the physical particle masses but attach a Langevin thermostat to the internal modes of the ring polymer during the dynamics. We justify this by showing analytically that the inclusion of an internal mode thermostat does not affect any of the desirable features of RPMD: thermostatted RPMD (TRPMD) is equally valid with respect to everything that has actually been proven about the method as RPMD itself. In particular, because of the choice of bead masses, the resulting method is still optimum in the short-time limit, and the transition state approximation to its reaction rate theory remains closely related to the semiclassical instanton approximation in the deep quantum tunneling regime. In effect, there is a continuous family of methods with these properties, parameterised by the strength of the Langevin friction. Here we explore numerically how the approximation to quantum dynamics depends on this friction, with a particular emphasis on vibrational spectroscopy. We find that a broad range of frictions approaching optimal damping give similar results, and that these results are immune to both the resonance problem of RPMD and the curvature problem of CMD.

References in corpus (13)

• Competing quantum effects in the dynamics of a flexible water model
• Efficient stochastic thermostatting of path integral molecular dynamics
• i-PI: A Python interface for ab initio path integral molecular dynamics simulations
• Colored-noise thermostats à la carte
• Full dimensional (15D) quantum-dynamical simulation of the protonated water-dimer II: infrared spectrum and vibrational dynamics
• Direct simulation of electron transfer using ring polymer molecular dynamics: Comparison with semiclassical instanton theory and exact quantum methods
• Derivation of a true (t->0+) quantum transition-state theory. Part I: Uniqueness and equivalence to ring-polymer molecular dynamics transition-state-theory
• Quantum fluctuations can promote or inhibit glass formation
• Kinetically Constrained Ring-Polymer Molecular Dynamics for Non-adiabatic Chemical Reactions
• Derivation of a true (t->0+) quantum transition-state theory. II. Recovery of the exact quantum rate in the absence of recrossing
• Isomorphic classical molecular dynamics model for an excess electron in a supercritical fluid
• Effective potential analytic continuation approach for real time quantum correlation functions involving nonlinear operators
• On the uniqueness of t->0+ quantum transition-state theory

Cited by in corpus (89)

• Nuclear quantum effects enter the mainstream
• i-PI 2.0: A Universal Force Engine for Advanced Molecular Simulations
• Committee neural network potentials control generalization errors and enable active learning
• Chemical Reaction Rates from Ring Polymer Molecular Dynamics: Theory and Practical Applications
• Boltzmann-conserving classical dynamics in quantum time-correlation functions: Matsubara dynamics
• A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat
• Relation of centroid molecular dynamics and ring-polymer molecular dynamics to exact quantum dynamics
• Elucidating the NuclearQuantum Dynamics of Intramolecular Double Hydrogen Transfer in Porphycene
• Ab initio molecular dynamics with nuclear quantum effects at classical cost: ring polymer contraction for density functional theory
• SchNetPack 2.0: A neural network toolbox for atomistic machine learning
• On the Consistency of Approximate Quantum Dynamics Simulation Methods for Vibrational Spectra in the Condensed Phase
• Fine Tuning Classical and Quantum Molecular Dynamics using a Generalized Langevin Equation
• Path-Integral Isomorphic Hamiltonian for Including Nuclear Quantum Effects in Non-adiabatic Dynamics
• Accurate molecular dynamics and nuclear quantum effects at low cost by multiple steps in real and imaginary time: using density functional theory to accelerate wavefunction methods
• Competing quantum effects in the free energy profiles and diffusion rates of hydrogen and deuterium molecules through clathrate hydrates
• Universal Machine Learning for the Response of Atomistic Systems to External Fields
• Simultaneous Deep Tunneling and Classical Hopping for Hydrogen Diffusion on Metals
• High Order Path Integrals Made Easy
• Stability of Complex Biomolecular Structures: Vander Waals, Hydrogen Bond Cooperativity, and Nuclear Quantum Effects
• i-PI 3.0: a flexible and efficient framework for advanced atomistic simulations
• Cayley modification for strongly stable path-integral and ring-polymer molecular dynamics
• Path-integral dynamics of water using curvilinear centroids
• An analysis of nonadiabatic ring-polymer molecular dynamics and its application to vibronic spectra
• Inexpensive modelling of quantum dynamics using path integral generalized Langevin equation thermostats
• Dimension-free path-integral molecular dynamics without preconditioning
• All-Electron, Real-Space Perturbation Theory for Homogeneous Electric Fields: Theory, Implementation, and Application within DFT
• Nuclear quantum effects in water exchange around lithium and fluoride ions
• "Divide and Conquer" Semiclassical Molecular Dynamics: A practical method for Spectroscopic calculations of High Dimensional Molecular Systems
• Quantum dynamics using path integral coarse-graining
• Hydrogen Diffusion and Trapping in α-Iron: The Role of Quantum and Anharmonic Fluctuations
• Critical role of quantum dynamical effects in the Raman spectroscopy of liquid water
• Decoding the spectroscopic features and timescales of aqueous proton defects
• Fully First-Principles Surface Spectroscopy with Machine Learning
• Incorporating Nuclear Quantum Effects in Molecular Dynamics with a Constrained Minimized Energy Surface
• Simplified Approach to the Mixed Time-averaging Semiclassical Initial Value Representation for the Calculation of Dense Vibrational Spectra
• Herman-Kluk propagator is free from zero-point energy leakage
• A generalized class of strongly stable and dimension-free T-RPMD integrators
• Vibrational strong coupling in liquid water from cavity molecular dynamics
• Progress and Challenges in Ab Initio Simulations of Quantum Nuclei in Weakly Bonded Systems
• Two-dimensional infrared-Raman spectroscopy as a probe of water's tetrahedrality
• Fast quasi-centroid molecular dynamics
• Nuclear Quantum Effects in liquid water at near classical computational cost using the adaptive Quantum Thermal Bath
• Reduced Rovibrational Coupling Cartesian Dynamics for Semiclassical Calculations: Application to the Spectrum of the Zundel Cation
• Highly efficient path-integral molecular dynamics simulations with GPUMD using neuroevolution potentials: Case studies on thermal properties of materials
• Routine Molecular Dynamics Simulations Including Nuclear Quantum Effects: from Force Fields to Machine Learning Potentials
• PYSED: A tool for extracting kinetic-energy-weighted phonon dispersion and lifetime from molecular dynamics simulations
• State Dependent Ring Polymer Molecular Dynamics for Investigating Excited Nonadiabatic Dynamics
• Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches
• Mean-field Matsubara dynamics: Analysis of path-integral curvature effects in rovibrational spectra
• NEP-MB-pol: A unified machine-learned framework for fast and accurate prediction of water's thermodynamic and transport properties
• Should Thermostatted Ring Polymer Molecular Dynamics be used to calculate thermal reaction rates?
• How to obtain thermostatted ring polymer molecular dynamics from exact quantum dynamics and when to use it
• Thermal quantum time-correlation functions from classical-like dynamics
• Testing the quasicentroid molecular dynamics method on gas-phase ammonia
• Developing machine-learned potentials to simultaneously capture the dynamics of excess protons and hydroxide ions in classical and path integral simulations
• Improving Condensed Phase Water Dynamics with Explicit Nuclear Quantum Effects: the Polarizable Q-AMOEBA Force Field
• Mode-coupling theory of lattice dynamics for classical and quantum crystals
• The dynamic nature of high pressure ice VII
• Multi-time Formulation of Matsubara Dynamics
• From Classical to Quantum and Back: Hamiltonian Adaptive Resolution Path Integral, Ring Polymer, and Centroid Molecular Dynamics
• Dissipative Tunneling Rates through the Incorporation of First-Principles Electronic Friction in Instanton Rate Theory I: Theory
• Effective Mass Path Integral Simulations of Quasiparticles in Condensed Phases
• Frustrated hops in Ring Polymer Surface Hopping: Real-time dynamics and detailed balance
• Simulating vibronic spectra via Matsubara dynamics: coping with the sign problem
• Fast quasi-centroid molecular dynamics for water and ice
• Improved torque estimator for condensed-phase quasicentroid molecular dynamics
• Path Integral Brownian Chain Molecular Dynamics: A Simple Approximation of Quantum Vibrational Dynamics
• Proton-transfer spectroscopy beyond the normal-mode scenario
• Hydrogen Bonding and Nuclear Quantum Effects in Clays
• Which Algorithm Best Propagates the Meyer-Miller-Stock-Thoss Mapping Hamiltonian for Non-Adiabatic Dynamics?
• Performant Automatic Differentiation of Local Coupled Cluster Theories: Response Properties and Ab Initio Molecular Dynamics
• A Fully Quantum-Mechanical Treatment for Kaolinite
• Comparison of Matsubara dynamics with exact quantum dynamics for an oscillator coupled to a dissipative bath
• The quantum nature of ubiquitous vibrational features revealed for ethylene glycol
• Consistency between the Green-Kubo formula and Lorentz model for predicting the infrared dielectric function of polar materials
• Predicting neutron experiments from first principles: A workflow powered by machine learning
• 2D spectroscopies from condensed phase dynamics: Accessing third-order response properties from equilibrium multi-time correlation functions
• Computing excited eigenstates using inexact Lanczos methods and tree tensor network states
• DL_POLY Quantum 2.1 software: A suite of real-time path integral methods for the simulation of dynamical properties and vibrational spectra
• Rényi entropy of quantum anharmonic chain at non-zero temperature
• Neural Network Potential Energy Surface for the low temperature Ring Polymer Molecular Dynamics of the H2CO + OH reaction
• Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab-initio water with near-zero added cost
• A tensor network representation of path integrals: Implementation and analysis
• Quantitative Molecular Simulations
• Accurate, full-dimensional computations of thousands of complex vibrational eigenstates with tree tensor network states
• Nuclear Quantum Effects in Multi-Step Condensed Matter Chemistry: A Path Integral Molecular Dynamics Study of Thermal Decomposition
• On the electronic path integral normal modes of the Meyer-Miller-Stock-Thoss representation of nonadiabatic dynamics
• A trigonometric integrator for the constrained ring polymer Hamiltonian dynamics
• Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy