Solving Linearized Equations of the $N$-body Problem Using the Lie-integration Method

Abstract

Several integration schemes exits to solve the equations of motion of the $N$-body problem. The Lie-integration method is based on the idea to solve ordinary differential equations with Lie-series. In the 1980s this method was applied for the $N$-body problem by giving the recurrence formula for the calculation of the Lie-terms. The aim of this works is to present the recurrence formulae for the linearized equations of motion of $N$-body systems. We prove a lemma which greatly simplifies the derivation of the recurrence formulae for the linearized equations if the recurrence formulae for the equations of motions are known. The Lie-integrator is compared with other well-known methods. The optimal step size and order of the Lie-integrator are calculated. It is shown that a fine-tuned Lie-integrator can be 30%-40% faster than other integration methods.

accepted for publication in MNRAS (13 pages, 4 figures); see http://cm.elte.hu/lie (cm.elte.hu/lie) for software

Cited by in corpus (7)

• Refined stellar, orbital and planetary parameters of the eccentric HAT-P-2 planetary system
• Numerical integration of dynamical systems with Lie series: Relativistic acceleration and non-gravitational forces
• Main-Belt Asteroids in the K2 Engineering Field of View
• Analysis of radial velocity variations in multiple planetary systems
• Dynamical stability of the Gliese 581 exoplanetary system
• Lie-series for orbital elements -- I. The planar case
• Lie-series for orbital elements -- II. The spatial case