A unified thermostat scheme for efficient configurational sampling for classical/quantum canonical ensembles via molecular dynamics

Abstract

We show a unified second-order scheme for constructing simple, robust and accurate algorithms for typical thermostats for configurational sampling for the canonical ensemble. When Langevin dynamics is used, the scheme leads to the BAOAB algorithm that has been recently investigated. We show that the scheme is also useful for other types of thermostat, such as the Andersen thermostat and Nosé-Hoover chain. Two 1-dimensional models and three typical realistic molecular systems that range from the gas phase, clusters, to the condensed phase are used in numerical examples for demonstration. Accuracy may be increased by an order of magnitude for estimating coordinate-dependent properties in molecular dynamics (when the same time interval is used), irrespective of which type of thermostat is applied. The scheme is especially useful for path integral molecular dynamics, because it consistently improves the efficiency for evaluating all thermodynamic properties for any type of thermostat.

arXiv commentManuscript + Supplementary Material. A unified second-order scheme for constructing simple, robust and accurate algorithms for typical thermostats for configurational sampling for the canonical ensemble. (Accepted by J. Chem. Phys.)

References in corpus (12)

• Canonical sampling through velocity-rescaling
• Accurate sampling using Langevin dynamics
• Efficient stochastic thermostatting of path integral molecular dynamics
• A simple and effective Verlet-type algorithm for simulating Langevin dynamics
• Robust and efficient configurational molecular sampling via Langevin 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
• Stochastic resonance-free multiple time-step algorithm for molecular dynamics with very large time steps
• Optimized Forest-Ruth- and Suzuki-like algorithms for integration of motion in many-body systems
• New Langevin and Gradient Thermostats for Rigid Body Dynamics
• High order three part split symplectic integrators: Efficient techniques for the long time simulation of the disordered discrete nonlinear Schroedinger equation

Cited by in corpus (15)

• Cayley modification for strongly stable path-integral and ring-polymer molecular dynamics
• Dimension-free path-integral molecular dynamics without preconditioning
• Critical role of quantum dynamical effects in the Raman spectroscopy of liquid water
• A generalized class of strongly stable and dimension-free T-RPMD integrators
• Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics
• Path Integral Molecular Dynamics for Exact Quantum Statistics of Multi-Electronic-State Systems
• Understanding Molecular Dynamics with Stochastic Processes via Real or Virtual Dynamics
• Stability of velocity-Verlet- and Liouville-operator-derived algorithms to integrate non-Hamiltonian systems
• Path Integral Brownian Chain Molecular Dynamics: A Simple Approximation of Quantum Vibrational Dynamics
• Electronic Friction Near Metal Surface: Incorporating Nuclear Quantum Effect with Ring Polymer Molecular Dynamics
• The Bayesian Inversion Problem for Thermal Average Sampling of Quantum Systems
• On the relation between the velocity- and position-Verlet integrators
• On the Numerical Stationary Distribution of Overdamped Langevin Equation in Harmonic System
• Geometric Phase Effect in Thermodynamic Properties and in the Imaginary-Time Multi-Electronic-State Path Integral Formulation
• Path-Integral Approach to Quantum Acoustics