On the accuracy of Stoermer/Verlet method as the numerical integrator of the n-body problem – Application to solar system


The Newton-Störmer/Verlet-leapfrog method (S/V) is a symplectic and symmetric one of order two, which, when applied to separable Hamiltonian dynamical systems, becomes explicit and conserves quadratic first integrals, e.g., the angular momentum in the N-body problem.

As for high accuracy long-time computations required by dynamical systems coming from astronomy, the order two is too low, we consider composed S/V methods up to order 10. Beyond all these remarkable qualities, our numerical experiments on the outer as well as inner Solar system (N-body problem for N larger than 5) show that, in practice, the conservation of the Hamiltonian of the flow, up to a prescribed order, still remains an open task. Specifically, the order at which the Hamiltonian is conserved seems to be quasi independent of the order of the composed method in spite of the fact that the compensated summation, as a specific technique to reduce the round-off errors, was used. In order to clear up this aspect we perform a round-off error analysis of S/V method as a one-step method applied to an arbitrary vector field.

Despite of the fact that we succeeded in improving the classical round-off error estimates, numerical experiments lead to the opinion that the accumulation of round-off errors bears the main responsibility for the lack of accuracy in long time energy conservation.


Calin-Ioan Gheorghiu,
Adrian-Cristian Muresan


N-body problem; composed Störmer/Verlet method.


See the expanding block below.

Cite this paper as

C.I. Gheorghiu, A.C. Mureşan, On the cccuracy of Stoermer/Verlet method as the numerical integrator of the n-body problem – Application to solar system, Rom. Astron. J., 16 (2006) 93-102.


About this paper


Romanian Astronomical Journal

Publisher Name

Romanian Academy

Print ISSN


Online ISSN


Google Scholar Profile

[1] Applegate, J. H. Douglas, M. R. Guersel, Y. Sussman, G. J. Wisdom, J: 1986,  Astron, J.,  92, 176
[2] Benetin, G., Giorgilli, A. : 1994, J. Statist. Phys.,74, 1117
[3] Butcher, J.: 2003,  Numerical Methods for Ordinary Differential Equations, J. Wiley, Chichester
[4] Dekker, T.J.: 1971, Numer. Math. 18, 224.
[5] Fukushima, T.: 2001, Astron. J., 121, 1768
[6] Gladman, B., Duncan, M., Candy, J: 1991, Celest. Mech., 52, 221
[7] Heirer, E., Lubicj, C. Wanner, G.: 2002, Geometric Numerical Integration – Structure Preserving Algorithms for Ordinary Differential Equations, Springer Verlag, Berlin, Heidelberg.
[8] Hairer, E., Mairer, M.: 2002, http://www.unige.ch/math/folks
[9] Hairer, E., Lubich, C., Wanner, G.: Acta Numerica, 1
[10] Henrici, P: 1962, Discrete Variable Methods for Ordinary Differential Equations, J. Wiley, New York, London
[11] Higham, N.J.: 1993, SIAM J. Sci. Comput. 14, 787
[12] Huang, W., Leimkuhler, B.: 1997, SIAM J. Sci. Comput., 18, 239
[13] Quinn, T., Tremaine, S: 1990, Astron. J., 99, 1016
[14] Quinn, T., Tremaine, S., Duncan, M: 1991, Asstron. J., 101, 2287
[15] Sanz-Serna, J.M., Calvo, P.M.: 1994, Numerical Hamiltonian Problems, Chapman and Hall, London.
[16] Stuart, A.M., Humphries, A.R.: 1996, Dynamical systems and Numerical Analysis, Cambridge University Press. Cambridge, UK.
[17] Sussman, G.J., Wisdom, J: 1992, Science, 257, 56
[18] Wisdom, J., Holman, M.: Astron. J., 102, 1528
[19] Yosida, H.: 1990, Phys. Lett., A, 150, 262.
[20] Yosida, H: 1993, Celest Mech. Dyn. Astron., 56, 27
[21] Zhang, M., Skeel, R. D.: 1997, Appl. Num. Meth., 25, 297


Related Posts