5. Integration method for Hamiltonian system.
Consider the integration of Hamiltonian system in which the energy is conserved.
• In many of formulas (e.g. the classical RK4), the errors in conserved
quantities (energy, angular momentum) accumulate in time.
• A numerical integration scheme that conserves the energy may besuitable for solving particular problems.
ex) A long-term stability of the planetary system.
• Symplectic integration method and Symmetric integration method are
known to have this property.

• Comparison of 2nd order symplectic method and Runge-Kutta method• Hamiltonian system has pure imaginary eigenvalue,• When a stable numerical method is applied, an oscillatory solution often
tends to be dumped or diverges.
• If the region of stability of the method includes the imaginary axis, thenumerical solution oscillates correctly.
Cf.) Some methods suitable for 2nd order ODEs.
• Stormer-Cowell method.

Extremely high precision formulas O(h14), O(h16) are used in thecelestial mechanics. Amount of floating operations, memory, areabout a half of Adams methods.
Ex) Verlet Method :• Runge-Kutta-Nystroem method.
Better accuracy for the same level RK formula, less memory.
Exc 5-1) Show that the Verlet method and leap flog method are the same.

Hamiltonian system.

Hamilton’sequation.

• Rewriting Hamilton’s equation,a formal solution is written
The operator exp[tDH] is well defined in the following sense;

Symplectic formula.

Exponential law for the non-commutative operators.

• Symplectic integration formula approximate this infinite products byfinite products as

The number of product, k, is taken large enough to have a desirable ordertruncation error O(hn+1).
(For this formula, the local truncation error is τ = O(hn).)
Exc 5-2) Show that the 2nd order symplectic formula is the same as
the Verlet method.
Exc 5-3) Check if the symplectic formulas generated by the symplectic maps
are the same.

Why the symplectic formula conserves energy?

This means that the symplectic formula is associated with the Hamiltonianof the different dynamical system.
Exc 5-4) Integrate harmonic oscillator using some symplectic formulas, aswell as the other non-symplectic formula of the same order.
Then compare the conservation of the energy in time.

Some systematic decomposition method for the exponential products.

(1)

(2)

(3)

Formulas (2) and (3) are symmetric, Sn(h) Sn(-h) = 1

Symplectic formula with symmetric decomposition.

Formulas (2) and (3) becomes(2)(3)

Symmetric integration formulas
A few drawback of the symplectic formula– Difficult to change the time step size h.
Once change the step size, the formula is no more symplectic.
– floating operation is way more than the multistep method ofcomparable order.

Ex) Trapezoidal formula

Symmetric formula– the error in the energy is bounded similar to symplectic formula.
– No difficulty in changing the time step size h.

Symmetric integration formulas (continue)– An implicit formula. Symmetric.
– Direct iteration can be used.
– relatively simple coding and decent accuracy.
Exc 5-5) Derive the above 2nd order Hermite formula– Method for changing step size in symmetric formula.
– Stability of symmetric formula – linear stability and P-stability.