Symplectic Runge-kutta Schemes I: Order Conditions

Much recent work has indicated that considerable beneet arises from the use of symplectic algorithms when numerically integrating Hamiltonian systems of diierential equations. Runge-Kutta schemes are symplectic subject to a simple algebraic condition. Starting with Butcher's formalism it is shown that there exists a more natural basis for the set of necessary and suucient order conditions for these methods; involving only s(s + 1)=2 free parameters for a symplectic s stage scheme. A graph theoretical process for determining the new order conditions is outlined. Furthermore, it is shown that any rooted tree arising from the same free tree enforces the same algebraic constraint on the parametrised coeecients. When coupled with the standard simplifying assumptions for implicit schemes the number of order conditions may be further reduced. In the new framework a simple symmetry of the parameter matrix yields (not necessarily symplectic) self-adjoint methods. In this case the order conditions associated with even trees become redundant.

[1]  Arieh Iserles,et al.  A unified approach to spurious solutions introduced by time discretization. Part I: basic theory , 1991 .

[2]  Kevin Burrage,et al.  A special family of Runge-Kutta methods for solving stiff differential equations , 1978 .

[3]  T. Chow,et al.  Classical Mechanics, Second Edition , 1995 .

[4]  W. Oevel,et al.  Symplectic Runge-kutta Schemes Iii: Canonical Elementary Differentials , 1993 .

[5]  E. Hairer Backward analysis of numerical integrators and symplectic methods , 1994 .

[6]  J. M. Sanz-Serna,et al.  Symplectic integrators for Hamiltonian problems: an overview , 1992, Acta Numerica.

[7]  M. Sofroniou Symbolic and numerical methods for Hamiltonian systems , 1994 .

[8]  V. Arnold Mathematical Methods of Classical Mechanics , 1974 .

[9]  J. Verwer,et al.  Stability of Runge-Kutta Methods for Stiff Nonlinear Differential Equations , 1984 .

[10]  R. Russell,et al.  Unitary integrators and applications to continuous orthonormalization techniques , 1994 .

[11]  R. Scherer REDUCIBILITY AND CHARACTERIZATION OF SYMPLECTIC RUNGE{KUTTA METHODS , 1994 .

[12]  P. Albrecht The Runge-Kutta Theory in a Nutshell , 1996 .

[13]  Mark Sofroniou,et al.  Symbolic Derivation of Runge-Kutta Methods , 1994, J. Symb. Comput..

[14]  John Guckenheimer,et al.  Numerical Analysis of Dynamical Systems , 1999 .

[15]  F. Lasagni Canonical Runge-Kutta methods , 1988 .

[16]  Geetha Ramaswami,et al.  Perturbed collocation and symplectic RKN methods , 1995, Adv. Comput. Math..

[17]  M. E. Hosea A new recurrence for computing Runge-Kutta truncation error coefficients , 1995 .

[18]  J. Butcher Implicit Runge-Kutta processes , 1964 .

[19]  H. Munthe-Kaas Lie-Butcher theory for Runge-Kutta methods , 1995 .

[20]  J. M. Sanz-Serna,et al.  Order conditions for canonical Runge-Kutta schemes , 1991 .

[21]  Butcher's simplifying assumption for symplectic integrators , 1992 .