Constant pressure and temperature discrete-time Langevin molecular dynamics

Abstract

We present a new and improved method for simultaneous control of temperature and pressure in molecular dynamics simulations with periodic boundary conditions. The thermostat-barostat equations are build on our previously developed stochastic thermostat, which has been shown to provide correct statistical configurational sampling for any time step that yields stable trajectories. Here, we extend the method and develop a set of discrete-time equations of motion for both particle dynamics and system volume in order to seek pressure control that is insensitive to the choice of the numerical time step. The resulting method is simple, practical, and efficient. The method is demonstrated through direct numerical simulations of two characteristic model systems - a one dimensional particle chain for which exact statistical results can be obtained and used as benchmarks, and a three dimensional system of Lennard-Jones interacting particles simulated in both solid and liquid phases. The results, which are compared against the method of Kolb & Dunweg, show that the new method behaves according to the objective, namely that acquired statistical averages and fluctuations of configurational measures are accurate and robust against the chosen time step applied to the simulation.

arXiv comment14 pages, 12 figures

References in corpus (3)

• A simple and effective Verlet-type algorithm for simulating Langevin dynamics
• Application of the G-JF Discrete-Time Thermostat for Fast and Accurate Molecular Simulations
• Optimized Constant Pressure Stochastic Dynamics

Cited by in corpus (12)

• Pressure control using stochastic cell rescaling
• Complete Set of Stochastic Verlet-Type Thermostats for Correct Langevin Simulations
• New Langevin and Gradient Thermostats for Rigid Body Dynamics
• A random batch Ewald method for charged particles in the isothermal-isobaric ensemble
• Sampling the isothermal-isobaric ensemble by Langevin dynamics
• Elongational viscosity of weakly entangled polymer melt via coarse-grained molecular dynamics simulation
• Accurate Configurational and Kinetic Statistics in Discrete-Time Langevin Systems
• On the connection between dissipative particle dynamics and the Itô-Stratonovich dilemma
• Comparing four hard-sphere approximations for the low-temperature WCA melting line
• The G-JF Thermostat for Accurate Configurational Sampling in Soft-Matter Simulations
• Stochastic sampling of the isothermal-isobaric ensemble: phase diagram of crystalline solids from molecular dynamics simulation
• Comparing zero-parameter theories for the WCA and harmonic-repulsive melting lines