Langevin Thermostat for Rigid Body Dynamics

Abstract

We present a new method for isothermal rigid body simulations using the quaternion representation and Langevin dynamics. It can be combined with the traditional Langevin or gradient (Brownian) dynamics for the translational degrees of freedom to correctly sample the NVT distribution in a simulation of rigid molecules. We propose simple, quasi-symplectic second-order numerical integrators and test their performance on the TIP4P model of water. We also investigate the optimal choice of thermostat parameters.

15 pages, 13 figures, 1 table

References in corpus (1)

• Accurate sampling using Langevin dynamics

Cited by in corpus (12)

• DNA hybridization kinetics: zippering, internal displacement and sequence dependence
• Coarse-graining DNA for simulations of DNA nanotechnology
• Stable and fast semi-implicit integration of the stochastic Landau-Lifshitz equation
• Brownian Dynamics of Confined Rigid Bodies
• New Langevin and Gradient Thermostats for Rigid Body Dynamics
• DNA nanotechnology: understanding and optimisation through simulation
• Discretization errors in molecular dynamics simulations with deterministic and stochastic thermostats
• Designing the Optimal Bit: Balancing Energetic Cost, Speed and Reliability
• Probing the microscopic structure and flexibility of oxidized DNA by molecular simulations
• Geometric integrator for Langevin systems with quaternion-based rotational degrees of freedom and hydrodynamic interactions
• Refinement of Thermostated Molecular Dynamics Using Backward Error Analysis
• Effects of Structural Inhomogeneity on Equilibration Processes in Langevin Dynamics