This paper aims to give some explanation and examples to show that numerical schemes applied to nonlinear problems can have convergent modified (or interpolating) Hamiltonians (MHs). It is pointed out, by means of specific and rather simple examples, that there are systems where the MH has a closed form. These systems are related to integrable mappings. For a special illustrative example given here the MHs, of the numerical approximations arising from the symplectic Euler and leap-frog methods, are presented in closed form. This MH is also given (in orders of the time step) by Yoshida's scheme (which uses in an essential way the Baker–Campbell–Hausdorff series). It is shown that all symplectic integrations of this example, which belong to the class presented by Yoshida, preserve the original energy of the system and admit closed form MHs. A second Hamiltonian example is also introduced. Upon application of the symplectic Euler method to this system, we obtain an integrable mapping which exactly conserves a biquadratic invariant. Hence it is also possible to derive a closed form for the MH in this case. In the light of these examples, the potential importance of integrable mappings to the field of symplectic integration is discussed.
Read full abstract