RUS  ENG
Full version
JOURNALS // Nanosystems: Physics, Chemistry, Mathematics // Archive

Nanosystems: Physics, Chemistry, Mathematics, 2013 Volume 4, Issue 5, Pages 657–669 (Mi nano805)

Accurate energy conservation in molecular dynamics simulation

O. A. Zolotov, V. E. Zalizniak

Siberian Federal University, 79 Svobodny Prospect, Krasnoyarsk 660041, Russia

Abstract: In molecular dynamics, Hamiltonian systems of differential equations are numerically integrated using some symplectic method. Symplectic integrators are simple algorithms that appear to be wellsuited for large scale simulations. One feature of these simulations is that there is an unphysical drift in the energy of the system over long integration periods. A drift in the energy is more obvious when a relatively long time step is used. In this article, a special approach, based on symplectic discretization and momenta corrections, is presented. The proposed method conserves the total energy of the system over the interval of simulation for any acceptable time step. A new approach to perform a constant-temperature molecular dynamics simulation is also presented. Numerical experiments illustrating these approaches are described.

Keywords: Hamiltonian systems, symplectic numerical methods, energy conservation, molecular dynamics.

PACS: 02.70.Ns, 45.20.dh, 45.50.Jf

Language: English



Bibliographic databases:


© Steklov Math. Inst. of RAS, 2024