Abstract
The discrete dynamics obtained by a simple second-order map is shown to lead to an exact integration of harmonic modes with a conserved quantity: a ‘‘shadow Hamiltonian.’’ The exact relations for ‘‘velocity’’ and energy are used to correct the energy obtained by the second-order algorithm commonly used in molecular dynamics simulations for complex many-body systems and lead to an improved energy conservation with a factor of 4 to 5.