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.
Introduction Symplectic integrators are designed for the numerical solution of Hamilton's equations, which read
p ˙ = − ∂ H ∂ q and q ˙ = ∂ H ∂ p , {\displaystyle {\dot {p}}=-{\frac {\partial H}{\partial q}}\quad {\mbox{and}}\quad {\dot {q}}={\frac {\partial H}{\partial p}},}
where q {\displaystyle q} denotes the position coordinates, p {\displaystyle p} the momentum coordinates, and H {\displaystyle H} is the Hamiltonian. The set of position and momentum coordinates ( q , p ) {\displaystyle (q,p)} are called canonical coordinates. (See Hamiltonian mechanics for more background.) The time evolution of Hamilton's equations is a symplectomorphism, meaning that it conserves the symplectic 2-form d p ∧ d q {\displaystyle dp\wedge dq} . A numerical scheme is a symplectic integrator if it also conserves this 2-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, such as the primitive Euler scheme and the classical Runge–Kutta scheme, are not symplectic integrators.
Methods for constructing symplectic algorithms
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 z = ( q , p ) {\displaystyle z=(q,p)} to denote the canonical coordinates including both the position and momentum coordinates. Then, the set of the Hamilton's equations given in the introduction can be expressed in a single expression as
where { ⋅ , ⋅ } {\displaystyle \{\cdot ,\cdot \}} is a Poisson bracket. Furthermore, by introducing an operator D H ⋅ = { ⋅ , H } {\displaystyle D_{H}\cdot =\{\cdot ,H\}} , which returns a Poisson bracket of the operand with the Hamiltonian, the expression of the Hamilton's equation can be further simplified to
z ˙ = D H z . {\displaystyle {\dot {z}}=D_{H}z.}
The formal solution of this set of equations at time t {\displaystyle t} is given as a matrix exponential:
Note the positivity of t D H {\displaystyle tD_{H}} in the matrix exponential. When the Hamiltonian has the form of equation (1), the solution (3) is equivalent to
The SI scheme approximates the time-evolution operator exp [ t ( D T + D V ) ] {\displaystyle \exp[t(D_{T}+D_{V})]} in the formal solution (4) by a product of operators as
… excerpt ends here. Continue reading the full article.
