Metric-tensor flexible-cell algorithm for isothermal-isobaric molecular dynamics simulations

Abstract

An extended Hamiltonian approach to conduct isothermal-isobaric molecular dynamics simulations with full cell flexibility is presented. The components of the metric tensor are used as the fictitious degrees of freedom for the cell, thus avoiding the problem of spurious cell rotations and artificial symmetry breaking effects present in the original Parrinello-Rahman scheme. This is complemented by the Nose-Poincare approach for isothermal sampling. The combination of these two approaches leads to equations of motion that are Hamiltonian in structure, and which can therefore be solved numerically using recently developed powerful symplectic integrators. One such integrator, the generalised leap-frog, is employed to provide a numerical algorithm for integrating the isothermal-isobaric equations of motion obtained.

REVTeX 4.0, 10 pages, 7 figures, to appear in Journal of Chemical Physics

References in corpus (1)

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

Cited by in corpus (6)

• Temperature of the inner-core boundary of the Earth: Melting of iron at high pressure from first-principles coexistence simulations
• Proton Transfer in Phase IV of Solid Hydrogen and Deuterium
• Room-Temperature Structures of Solid Hydrogen at High Pressures
• Phase diagram of silicon from atomistic simulations
• High-temperature behaviour of supported graphene: electron-phonon coupling and substrate-induced doping
• Ab initio Simulations of Superionic H2O, H2O2, and H9O4 Compounds