Symplectic algorithm for constant-pressure molecular dynamics using a Nosé–Poincaré thermostat

We present a new algorithm for isothermal–isobaric molecular-dynamics simulation. The method uses an extended Hamiltonian with an Andersen piston combined with the Nose–Poincare thermostat, recently developed by Bond, Leimkuhler, and Laird [J. Comp. Phys. 151, 114 (1999)]. This Nose–Poincare–Andersen (NPA) formulation has advantages over the Nose-Hoover-Andersen approach in that the NPA is Hamiltonian and can take advantage of symplectic integration schemes, which lead to enhanced stability for long-time simulations. The equations of motion are integrated using a generalized leapfrog algorithm (GLA) and the method is easy to implement, symplectic, explicit, and time reversible. To demonstrate the superior stability of the method we show results for test simulations using a model for aluminum and compare it to a recently developed time-reversible algorithm for Nose–Hoover–Anderson. In addition, an extension of the NPA to multiple time steps is outlined and a symplectic and time-reversible integration algorithm, based on the GLA, is given.

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

[2]  J. M. Sanz-Serna,et al.  Numerical Hamiltonian Problems , 1994 .

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

[4]  J. M. Sanz-Serna,et al.  Runge-kutta schemes for Hamiltonian systems , 1988 .

[5]  Denis J. Evans,et al.  Non-Newtonian molecular dynamics , 1984 .

[6]  J. M. Sanz-Serna,et al.  Symplectic integrators for Hamiltonian problems: an overview , 1992, Acta Numerica.

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

[8]  Mark E. Tuckerman,et al.  Explicit reversible integrators for extended systems dynamics , 1996 .

[9]  Hoover,et al.  Constant-pressure equations of motion. , 1986, Physical review. A, General physics.

[10]  Shuichi Nosé,et al.  Constant Temperature Molecular Dynamics Methods , 1991 .

[11]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[12]  Davenport,et al.  Free-energy calculations and the melting point of Al. , 1992, Physical review. B, Condensed matter.

[13]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[14]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

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

[16]  Wilma K. Olson,et al.  The dependence of DNA tertiary structure on end conditions: Theory and implications for topological transitions , 1994 .

[17]  Y. Suris,et al.  The canonicity of mappings generated by Runge-Kutta type methods when integrating the systems x¨ = - 6 U/ 6 x , 1989 .

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

[19]  W. C. Swope,et al.  A computer simulation method for the calculation of equilibrium constants for the formation of physi , 1981 .

[20]  Ray,et al.  Isothermal molecular-dynamics ensembles. , 1988, Physical review. A, General physics.

[21]  F. Lasagni.  Canonical Runge-Kutta methods , 1988 .

[22]  J. H. R. Clarke,et al.  A comparison of constant energy, constant temperature and constant pressure ensembles in molecular dynamics simulations of atomic liquids , 1984 .