Stochastic Dynamics with Correct Sampling for Constrained Systems.

In this paper we discuss thermostatting using stochastic methods for molecular simulations where constraints are present. For so-called impulsive thermostats, like the Andersen thermostat, the equilibrium temperature will differ significantly from the imposed temperature when a limited number of particles are picked and constraints are applied. We analyze this problem and give two rigorous solutions for it. A correct general treatment of impulsive stochastic thermostatting, including pairwise dissipative particle dynamics and stochastic forcing in the presence of constraints, is given and it is shown that the constrained canonical distribution is sampled rigorously. We discuss implementation issues such as second order Trotter expansions. The method is shown to rigorously maintain the correct temperature for the case of extended simple point charge (SPC/E) water simulations.

[1]  H. Berendsen,et al.  Efficient Algorithms for Langevin and DPD Dynamics. , 2012, Journal of chemical theory and computation.

[2]  Gregory H. Miller,et al.  A Second-Order Strong Method for the Langevin Equations with Holonomic Constraints , 2011, SIAM J. Sci. Comput..

[3]  M. Parrinello,et al.  Accurate sampling using Langevin dynamics. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[4]  Eric Vanden-Eijnden,et al.  Second-order integrators for Langevin equations with holonomic constraints , 2006 .

[5]  P. Coveney,et al.  A stochastic Trotter integration scheme for dissipative particle dynamics , 2006, Math. Comput. Simul..

[6]  Eajf Frank Peters,et al.  Elimination of time step effects in DPD , 2004 .

[7]  Tony Shardlow,et al.  Splitting for Dissipative Particle Dynamics , 2002, SIAM J. Sci. Comput..

[8]  C. Lowe.  An alternative approach to dissipative particle dynamics , 1999 .

[9]  Roger B. Sidje,et al.  Expokit: a software package for computing matrix exponentials , 1998, TOMS.

[10]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

[11]  B. Leimkuhler,et al.  Symplectic Numerical Integrators in Constrained Hamiltonian Systems , 1994 .

[12]  J. Koelman,et al.  Dynamic simulations of hard-sphere suspensions under steady shear , 1993 .

[13]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[14]  H. Berendsen,et al.  A LEAP-FROG ALGORITHM FOR STOCHASTIC DYNAMICS , 1988 .

[15]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[16]  Giovanni Ciccotti,et al.  Andersen's canonical-ensemble molecular dynamics for molecules with constraints , 1986 .

[17]  Herman J. C. Berendsen,et al.  ALGORITHMS FOR BROWNIAN DYNAMICS , 1982 .

[18]  H. Berendsen,et al.  STOCHASTIC DYNAMICS FOR MOLECULES WITH CONSTRAINTS BROWNIAN DYNAMICS OF NORMAL-ALKANES , 1981 .

[19]  D. Ermak,et al.  Numerical integration of the Langevin equation: Monte Carlo simulation , 1980 .

[20]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[21]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[22]  Berk Hess,et al.  P-LINCS:  A Parallel Linear Constraint Solver for Molecular Simulation. , 2008, Journal of chemical theory and computation.