Efficient stochastic thermostatting of path integral molecular dynamics

Abstract

The path integral molecular dynamics (PIMD) method provides a convenient way to compute the quantum mechanical structural and thermodynamic properties of condensed phase systems at the expense of introducing an additional set of high-frequency normal modes on top of the physical vibrations of the system. Efficiently sampling such a wide range of frequencies provides a considerable thermostatting challenge. Here we introduce a simple stochastic path integral Langevin equation (PILE) thermostat which exploits an analytic knowledge of the free path integral normal mode frequencies. We also apply a recently-developed colored-noise thermostat based on a generalized Langevin equation (GLE), which automatically achieves a similar, frequency-optimized sampling. The sampling efficiencies of these thermostats are compared with that of the more conventional Nosé-Hoover chain (NHC) thermostat for a number of physically relevant properties of the liquid water and hydrogen-in-palladium systems. In nearly every case, the new PILE thermostat is found to perform just as well as the NHC thermostat while allowing for a computationally more efficient implementation. The GLE thermostat also proves to be very robust delivering a near-optimum sampling efficiency in all of the cases considered. We suspect that these simple stochastic thermostats will therefore find useful application in many future PIMD simulations.

Accepted for publication on JCP

References in corpus (7)

• Canonical sampling through velocity-rescaling
• Competing quantum effects in the dynamics of a flexible water model
• Accurate sampling using Langevin dynamics
• Nuclear quantum effects in solids using a colored-noise thermostat
• Colored-noise thermostats à la carte
• Langevin equation with colored noise for constant-temperature molecular dynamics simulations
• Stochastic thermostats: comparison of local and global schemes

Cited by in corpus (112)

• SchNet - a deep learning architecture for molecules and materials
• Nuclear quantum effects enter the mainstream
• i-PI 2.0: A Universal Force Engine for Advanced Molecular Simulations
• How to remove the spurious resonances from ring polymer molecular dynamics
• i-PI: A Python interface for ab initio path integral molecular dynamics simulations
• Accelerating the convergence of path integral dynamics with a generalized Langevin equation
• Committee neural network potentials control generalization errors and enable active learning
• Simulation and understanding of quantum crystals
• Automated Fitting of Neural Network Potentials at Coupled Cluster Accuracy: Protonated Water Clusters as Testing Ground
• A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat
• Efficient methods and practical guidelines for simulating isotope effects
• Pressure and stress tensor of complex anharmonic crystals within the stochastic self-consistent harmonic approximation
• Unraveling quantum mechanical effects in water using isotopic fractionation
• Quantum fluctuations can promote or inhibit glass formation
• Efficient multiple time scale molecular dynamics: using colored noise thermostats to stabilize resonances
• 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
• A unified thermostat scheme for efficient configurational sampling for classical/quantum canonical ensembles via molecular dynamics
• On the Consistency of Approximate Quantum Dynamics Simulation Methods for Vibrational Spectra in the Condensed Phase
• 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
• High Order Path Integrals Made Easy
• Equation of state of fluid methane from first principles with machine learning potentials
• TurboRVB: a many-body toolkit for {\it ab initio} electronic simulations by quantum Monte Carlo
• 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
• Frontiers, challenges, and solutions in modeling of swift heavy ion effects in materials
• Efficient and versatile model for vibrational STEM-EELS
• Path Integral Molecular Dynamics within the Grand Canonical-like Adaptive Resolution Technique: Simulation of Liquid Water
• From Classical to Quantum Glasses with Ultracold Polar Molecules
• Dimension-free path-integral molecular dynamics without preconditioning
• Nuclear quantum effects in water exchange around lithium and fluoride ions
• The phase diagram of water from quantum simulations
• Simulation of quantum zero-point effects in water using a frequency-dependent thermostat
• Solvent Fluctuations and Nuclear Quantum Effects Modulate the Molecular Hyperpolarizability of Water
• The inefficiency of re-weighted sampling and the curse of system size in high order path integration
• Theory and simulations of quantum glass forming liquids
• Analytic continuation of Wolynes theory into the Marcus inverted regime
• Unravelling the influence of quantum proton delocalization on electronic charge transfer through the hydrogen bond
• Hydrogen Diffusion and Trapping in α-Iron: The Role of Quantum and Anharmonic Fluctuations
• Direct path integral estimators for isotope fractionation ratios
• Decoding the spectroscopic features and timescales of aqueous proton defects
• On the calculation of quantum mechanical electron transfer rates
• Converged Colored Noise Path Integral Molecular Dynamics Study of the Zundel Cation down to Ultra-low Temperatures at Coupled Cluster Accuracy
• A generalized class of strongly stable and dimension-free T-RPMD integrators
• The MD17 Datasets from the Perspective of Datasets for Gas-Phase "Small" Molecule Potentials
• Two-dimensional infrared-Raman spectroscopy as a probe of water's tetrahedrality
• Progress and Challenges in Ab Initio Simulations of Quantum Nuclei in Weakly Bonded Systems
• Modeling the structural and thermal properties of loaded metal-organic frameworks. An interplay of quantum and anharmonic fluctuations
• Highly efficient path-integral molecular dynamics simulations with GPUMD using neuroevolution potentials: Case studies on thermal properties of materials
• Phonon-induced localization of excitons in molecular crystals from first principles
• Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics
• Accelerated ab-initio Molecular Dynamics: probing the weak dispersive forces in dense liquid hydrogen
• State Dependent Ring Polymer Molecular Dynamics for Investigating Excited Nonadiabatic Dynamics
• Path integral molecular dynamics with surface hopping for thermal equilibrium sampling of nonadiabatic systems
• Mean-field Matsubara dynamics: Analysis of path-integral curvature effects in rovibrational spectra
• Path Integral Molecular Dynamics for Exact Quantum Statistics of Multi-Electronic-State Systems
• UV-Visible Absorption Spectra of Solvated Molecules by Quantum Chemical Machine Learning
• Thermal quantum time-correlation functions from classical-like dynamics
• 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
• 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
• Ab Initio Generalized Langevin Equation
• Quasi-classical approaches to vibronic spectra revisited
• Reducing the cost of neural network potential generation for reactive molecular systems
• Nuclear Quantum Effects on Autoionization of Water Isotopologues Studied by Ab Initio Path Integral Molecular Dynamics
• Accelerating ab initio path integral molecular dynamics with multilevel sampling of potential surface
• Isotope quantum effects in the metallization transition in liquid hydrogen
• Probing anharmonic phonons by quantum correlators: A path integral approach
• From Classical to Quantum and Back: Hamiltonian Adaptive Resolution Path Integral, Ring Polymer, and Centroid Molecular Dynamics
• The Role of Quantum Fluctuations in the Hexatic Phase of Cold Polar Molecules
• Effects of chain length on Rouse modes and non-Gaussianity in linear and ring polymer melts
• Fast-Forward Langevin Dynamics with Momentum Flips
• Nuclear quantum effects in thermal conductivity from centroid molecular dynamics
• Frustrated hops in Ring Polymer Surface Hopping: Real-time dynamics and detailed balance
• Simulating nuclear and electronic quantum effects in enzymes
• Accurate and efficient machine learning interatomic potentials for finite temperature modeling of molecular crystals
• Metatensor and metatomic: foundational libraries for interoperable atomistic machine learning
• Uncertainty quantification for generalized Langevin dynamics
• Deciphering High-order Structural Correlations within Fluxional Molecules from Classical and Quantum Configurational Entropy
• Hydrogen Bonding and Nuclear Quantum Effects in Clays
• Stability and distortion of fcc-LaH${}_{10}$ with path-integral molecular dynamics
• Accelerated sampling by infinite swapping of path integral molecular dynamics with surface hopping
• Coarse-graining conformational dynamics with multi-dimensional generalized Langevin equation: how, when, and why
• Electronic Friction Near Metal Surface: Incorporating Nuclear Quantum Effect with Ring Polymer Molecular Dynamics
• Coexistence of Superconductivity and Superionicity in Li${}_2$MgH${}_{16}$
• A Fully Quantum-Mechanical Treatment for Kaolinite
• Comparison of Matsubara dynamics with exact quantum dynamics for an oscillator coupled to a dissipative bath
• Efficient Sampling of Thermal Averages of Interacting Quantum Particle Systems with Random Batches
• Revisiting the properties of superfluid and normal liquid ${}^4$He using ab initio potentials
• Continuum limit and preconditioned Langevin sampling of the path integral molecular dynamics
• Stochastic path-integral approach for predicting the superconducting temperatures of anharmonic solids
• Evaluating functions of positive-definite matrices using colored noise thermostats
• Harmonic Oscillator Staging Coordinates for Efficient Path Integral Simulations of Quantum Oscillators and Crystals
• When Theory Meets Experiment: What Does it Take to Accurately Predict ${}^1$H NMR Dipolar Relaxation Rates in Neat Liquid Water from Theory?
• Non-adiabatic Ring Polymer Molecular Dynamics in the Phase Space of the SU(N) Lie Group
• h-CMD: An efficient hybrid fast centroid and quasi-centroid molecular dynamics method for the simulation of vibrational spectra
• DL_POLY Quantum 2.1 software: A suite of real-time path integral methods for the simulation of dynamical properties and vibrational spectra
• Non-Markovian Effects in Quantum Rate Calculations of Hydrogen Diffusion with Electronic Friction
• Ergostatting and Thermostatting at a Fixed Point
• Self-learning path integral hybrid Monte Carlo with mixed ab initio and machine learning potentials for modeling nuclear quantum effects in water
• Anharmonicity and Coulomb pseudopotential effects on superconductivity in YH${}_6$ and YH${}_9$
• Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab-initio water with near-zero added cost
• 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
• Hydrogen phase-IV characterization by full account of quantum anharmonicity
• Multi-level Monte Carlo path integral molecular dynamics for thermal average calculation in the nonadiabatic regime
• NQCDynamics.jl: A Julia Package for Nonadiabatic Quantum Classical Molecular Dynamics in the Condensed Phase
• Open quantum systems beyond equilibrium: Lindblad equation and path integral molecular dynamics
• A trigonometric integrator for the constrained ring polymer Hamiltonian dynamics