Abstract
Based on Feng’s theory of formal vector fields and formal flows, we study the convergence problem of the formal energies of symplectic methods for Hamiltonian systems and give the clear growth of the coefficients in the formal energies. With the help of B-series and Bernoulli functions, we prove that in the formal energy of the mid-point rule, the coefficient sequence of the merging products of an arbitrarily given rooted tree and the bushy trees of height 1 (whose subtrees are vertices), approaches 0 as the number of branches goes to ∞; in the opposite direction, the coefficient sequence of the bushy trees of height m (m ⩾ 2), whose subtrees are all tall trees, approaches ∞ at large speed as the number of branches goes to +∞. The conclusion extends successfully to the modified differential equations of other Runge-Kutta methods. This disproves a conjecture given by Tang et al. (2002), and implies: (1) in the inequality of estimate given by Benettin and Giorgilli (1994) for the terms of the modified formal vector fields, the high order of the upper bound is reached in numerous cases; (2) the formal energies/formal vector fields are nonconvergent in general case.
Similar content being viewed by others
References
Benettin G, Giorgilli A. On the Hamiltonian interpolation of near to the identity symplectic mappings with application to symplectic integration algorithms. J Statist Phys, 1994, 74: 1117–1143
Calvo M P, Murua A, Sanz-Serna J M. Modified equations for ODEs. Contemp Math, 1994, 172: 63–74
Channell P J, Scovel J C. Symplectic integration of Hamiltonian systems. Nonlinearity, 1990, 3: 231–259
Feng K. Formal power series and numerical algorithms for dynamical systems. In: Proceedings of International Conference on Scientific Computation. Singapore: World Scientific, 1992, 28–35
Feng K. The calculus of generating functions and the formal energy for Hamiltonian algorithms. Collected Works of Feng Kang (II). Beijing: National Defence Industry Press, 1995, 284–302
Field C M, Nijhoff F W. A note on modified Hamiltonians for numerical integrations admitting an exact invariant. Nonlinearity, 2003, 16: 1673–1683
Hairer E, Wanner G. On the Butcher group and general multi-value methods. Computing, 1974, 13: 1–15
Hairer E. Backward analysis of numerical integrators and symplectic methods. Ann Numer Math, 1994, 1: 107–132
Hairer E, Lubich C, Wanner G. Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations. New York: Springer, 2002
Kuang J C. Applied Inequalities. Shandong: Shandong Science and Technology Press, 2010
Moser J K. Lectures on Hamiltonian systems. Mem Amer Math, 1968, 81: 1–60
Murua A. Formal series and numerical integrators, part i: Systems of ODEs and symplectic integrators. Appl Numer Math, 1999, 29: 221–251
Neishtadt A I. The separation of motions in systems with rapidly rotating phase. J Appl Math Mech, 1984, 48: 133–139
Sanz-Serna J M, Calvo M P. Numerical Hamiltonian Problems. London: Chapman and Hall, 1994
Tang Y F. Formal energy of a symplectic scheme for Hamiltonian systems and its application (I). Comput Math Appl, 1994, 27: 31–39
Tang Y F, Vázquez L, Zhang F, et al. Symplectic methods for the nonlinear schrödinger Equation. Comput Math Appl, 1996, 32: 73–83
Tang Y F, Xiao A G, Chen J B. Is the formal energy of the mid-point rule convergent? Comput Math Appl, 2002, 43: 1171–1181
Wang Z X, Guo D R. Introduction to special function (in Chinese). Beijing: Peking University Press, 2000
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Zhang, R., Tang, Y., Zhu, B. et al. Convergence analysis of the formal energies of symplectic methods for Hamiltonian systems. Sci. China Math. 59, 379–396 (2016). https://doi.org/10.1007/s11425-015-5003-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11425-015-5003-7