Explicit symplectic integrators of molecular dynamics algorithms for rigid-body molecules in the canonical, isobaric-isothermal, and related ensembles.

The authors propose explicit symplectic integrators of molecular dynamics (MD) algorithms for rigid-body molecules in the canonical and isobaric-isothermal ensembles. They also present a symplectic algorithm in the constant normal pressure and lateral surface area ensemble and that combined with the Parrinello-Rahman algorithm. Employing the symplectic integrators for MD algorithms, there is a conserved quantity which is close to Hamiltonian. Therefore, they can perform a MD simulation more stably than by conventional nonsymplectic algorithms. They applied this algorithm to a TIP3P pure water system at 300 K and compared the time evolution of the Hamiltonian with those by the nonsymplectic algorithms. They found that the Hamiltonian was conserved well by the symplectic algorithm even for a time step of 4 fs. This time step is longer than typical values of 0.5-2 fs which are used by the conventional nonsymplectic algorithms.

[1]  N. Matubayasi,et al.  Reversible molecular dynamics for rigid bodies and hybrid Monte Carlo , 1999 .

[2]  J. W. Humberston.  Classical mechanics , 1980, Nature.

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

[4]  Yuko Okamoto,et al.  Molecular dynamics simulations in the multibaric–multithermal ensemble , 2004 .

[5]  C. Poole,et al.  Classical Mechanics, 3rd ed. , 2002 .

[6]  S. Nosé.  A molecular dynamics method for simulations in the canonical ensemble , 1984 .

[7]  Yuko Okamoto,et al.  Multibaric–multithermal ensemble molecular dynamics simulations , 2006, J. Comput. Chem..

[8]  H. Yoshida.  Construction of higher order symplectic integrators , 1990 .

[9]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[10]  C. W. Gear,et al.  Numerical initial value problem~ in ordinary differential eqttations , 1971 .

[11]  G. W. Swift,et al.  Effect of the Quasiparticle Mean Free Path on Poiseuille Flow in Normal Liquid 3He , 1980 .

[12]  S. Nosé,et al.  Constant pressure molecular dynamics for molecular systems , 1983 .

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

[14]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[15]  Mitsunori Ikeguchi,et al.  Partial rigid‐body dynamics in NPT, NPAT and NPγT ensembles for proteins and membranes , 2004, J. Comput. Chem..

[16]  Thomas F. Miller,et al.  Symplectic quaternion scheme for biophysical molecular dynamics , 2002 .

[17]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[18]  Stephen D. Bond,et al.  The Nosé-Poincaré Method for Constant Temperature Molecular Dynamics , 1999 .

[19]  S. Nosé.  An Improved Symplectic Integrator for Nosé-Poincaré Thermostat , 2001 .