In mathematics, a symplectic integrator (SI) is a numerical integration scheme for Hamiltonian systems. Symplectic integrators form the subclass of geometric integrators which, by definition, are canonical transformations. They are widely used in nonlinear dynamics, molecular dynamics, discrete element methods, accelerator physics, plasma physics, quantum physics, and celestial mechanics.
Contents
Introduction
Symplectic integrators are designed for the numerical solution of Hamilton's equations, which read
where
The time evolution of Hamilton's equations is a symplectomorphism, meaning that it conserves the symplectic two-form
Symplectic integrators possess, as a conserved quantity, a Hamiltonian which is slightly perturbed from the original one. By virtue of these advantages, the SI scheme has been widely applied to the calculations of long-term evolution of chaotic Hamiltonian systems ranging from the Kepler problem to the classical and semi-classical simulations in molecular dynamics.
Most of the usual numerical methods, like the primitive Euler scheme and the classical Runge-Kutta scheme, are not symplectic integrators.
Splitting methods for separable Hamiltonians
A widely used class of symplectic integrators is formed by the splitting methods.
Assume that the Hamiltonian is separable, meaning that it can be written in the form
This happens frequently in Hamiltonian mechanics, with T being the kinetic energy and V the potential energy.
For the notational simplicity, let us introduce the symbol
where
The formal solution of this set of equations is given as a matrix exponential:
Note the positivity of
When the Hamiltonian has the form of eq. (1), the solution (3) is equivalent to
The SI scheme approximates the time-evolution operator
where
Since
By using a Taylor series,
where
In concrete terms,
and
Note that both of these maps are practically computable.
Splitting methods for general nonseparable Hamiltonians
General nonseparable Hamiltonians can also be explicitly and symplectically integrated.
To do so, Tao introduced a restraint that binds two copies of phase space together to enable an explicit splitting of such systems. The idea is, instead of
The new Hamiltonian is advantageous for explicit symplectic integration, because it can be split into the sum of three sub-Hamiltonians,
Examples
Several simplectic integrators are given below. An illustrative way to use them is to consider a particle with position
To apply a timestep with values
A first-order example
The symplectic Euler method is the first-order integrator with
Note that the algorithm above does not work if time-reversibility is needed. The algorithm has to be implemented in two parts, one for positive time steps, one for negative time steps.
A second-order example
The Verlet method is the second-order integrator with
Since
A third-order example
A third-order symplectic integrator (with
A fourth-order example
A fourth-order integrator (with
To determine these coefficients, the Baker–Campbell–Hausdorff formula can be used. Yoshida, in particular, gives an elegant derivation of coefficients for higher-order integrators. Later on, Blanes and Moan further developed partitioned Runge–Kutta methods for the integration of systems with separable Hamiltonians with very small error constants.