Symplectic algorithm for constant-pressure molecular dynamics using a Nose-Poincare thermostat

Abstract

We present a new algorithm for isothermal-isobaric molecular-dynamics simulation. The method uses an extended Hamiltonian with an Andersen piston combined with the Nos'e-Poincar'e thermostat, recently developed by Bond, Leimkuhler and Laird [J. Comp. Phys., 151, (1999)]. This Nos'e-Poincar'e-Andersen (NPA) formulation has advantages over the Nos'e-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 and the method is easy to implement, symplectic, explicit and time reversible. To demonstrate the stability of the method we show results for test simulations using a model for aluminum.

7 pages

Cited by in corpus (17)

• Pressure control using stochastic cell rescaling
• Irreversible reorganization in a supercooled liquid originates from localised soft modes
• Langevin Dynamics in Constant Pressure Extended Systems
• Adjusting the melting point of a model system via Gibbs-Duhem integration: application to a model of Aluminum
• Explicit symplectic approximation of nonseparable Hamiltonians: algorithm and long time performance
• Localised soft modes and the supercooled liquid's irreversible passage through its configuration space
• Metric-tensor flexible-cell algorithm for isothermal-isobaric molecular dynamics simulations
• Free Volume cannot Explain the Spatial Heterogeneity of Debye-Waller factors in a Glass-Forming Binary Alloy
• A random batch Ewald method for charged particles in the isothermal-isobaric ensemble
• Generating Generalized Distributions from Dynamical Simulation
• Ab initio mass tensor molecular dynamics
• The central role of thermal collective strain in the relaxation of structure in a supercooled liquid
• Semiexplicit Symplectic Integrators for Non-separable Hamiltonian Systems
• Quantitative Theory of a Relaxation Function in a Glass-Forming System
• Preservation of Quadratic Invariants by Semiexplicit Symplectic Integrators for Non-separable Hamiltonian Systems
• Structural Phases in Non-Additive Soft-Disk Mixtures: Glasses, Substitutional Order, and Random Tilings
• Explicit Symplectic Integrators for Massive Point Vortex Dynamics in Binary Mixture of Bose--Einstein Condensates