Numerical integration of variational equations for Hamiltonian systems with long range interactions

Helen Christodoulidi, Tassos Bountis, Lambros Drossos

Research output: Contribution to journalArticlepeer-review

3 Citations (Scopus)


We study numerically classical 1-dimensional Hamiltonian lattices involving inter-particle long range interactions that decay with distance like 1/rα, for α≥0. We demonstrate that although such systems are generally characterized by strong chaos, they exhibit an unexpectedly organized behavior when the exponent α<1. This is shown by computing dynamical quantities such as the maximal Lyapunov exponent, which decreases as the number of degrees of freedom increases. We also discuss our numerical methods of symplectic integration implemented for the solution of the equations of motion together with their associated variational equations. The validity of our numerical simulations is estimated by showing that the total energy of the system is conserved within an accuracy of 4 digits (with integration step τ=0.02), even for as many as N=8000 particles and integration times as long as 106 units.

Original languageEnglish
Pages (from-to)158-165
Number of pages8
JournalApplied Numerical Mathematics
Publication statusPublished - Jun 1 2016


  • Hamiltonian systems
  • Long range interactions
  • Symplectic integration
  • Variational equations

ASJC Scopus subject areas

  • Numerical Analysis
  • Computational Mathematics
  • Applied Mathematics

Fingerprint Dive into the research topics of 'Numerical integration of variational equations for Hamiltonian systems with long range interactions'. Together they form a unique fingerprint.

Cite this