NANOSYSTEMS: PHYSICS, CHEMISTRY, MATHEMATICS, 2013, 4 (5), P. 657–669
ACCURATE ENERGY CONSERVATION IN MOLECULAR DYNAMICS SIMULATION
O. A. Zolotov – Siberian Federal University, 79 Svobodny Prospect, Krasnoyarsk 660041, Russia; ozolot_@mail.ru
V. E. Zalizniak – Siberian Federal University, 79 Svobodny Prospect, Krasnoyarsk 660041, Russia; vzalizniak@sfu-kras.ru
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