×

Preservation and destruction of periodic orbits by symplectic integrators. (English) Zbl 1196.65191

The paper is interested in what happens to periodic orbits and lower dimensional invariant tori of a Hamiltonian system when they are discretized by a symplectic integrator. The main question addressed is whether the numerical solution preserve the periodic orbits of the original vectorfield or not.
The results are directly related with two cases studied in sections 2 and 3 of the paper. The first case considers non resonant step sizes, this is, when none of the integer multiples of the step size of the numerical integrator coincides with the period of the periodic orbit. In this case the periodic orbits generally persist but slightly perturbed. Similar to what happens with invariant tori in KAM theory, section 2 starts commenting a previous result of Shang where for full dimensional invariant tori the frequencies of the numerical algorithm are the same ones of the original system but multiplied by the step size. Then it continues with its goal of developing a similar result for the preservation of periodic orbits, and what the authors prove, is that under suitable hypothesis (essentially of analyticity, reducibility and isotropy) the result is the same as in the full dimensional case but with the set of frequencies shifted from the original ones. To prove this result, the main tool that the authors use is an embedding of the symplectic map produced by the numerical integrator into an analytic symplectic flow and a KAM result of Á. Jorba and J. Villanueva [J. Nonlinear Sci., 7, 427–473 (1997; Zbl 0898.58044)].
The resonant case is studied in section 3. Here it is assumed that the period of a periodic orbit divided by the step size of the numerical integrator is an integer (\(T/h=n\)). As it should be expected from the results of the previous section, in this case the periodic orbit is destroyed by the numerical integrator leaving \(n\) elliptic and \(n\) hyperbolic periodic points under the numerical map. This fact is clearly illustrated numerically considering some periodic orbits in the two degrees of freedom Hénon-Heiles system, but a full description of the numerical flow in the neighborhood of a periodic orbit under resonant symplectic discretization, and a proof of the behaviors conjectured in the paper remains to be undertaken.

MSC:

65P10 Numerical methods for Hamiltonian systems including symplectic integrators

Citations:

Zbl 0898.58044
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Arnold, V.I.: Proof of a theorem of A.N. Kolmogorov on the preservation of conditionally periodic motions under a small perturbation of the Hamiltonian. Ups. Mat. Nauk 18(5), 91–192 (1963)
[2] Banyaga, A.: The structure of classical diffeomorphism groups. In: Mathematics and Its Applications, vol. 400. Kluwer Academic, Boston (1997) · Zbl 0874.58005
[3] Benettin, G., Giorgilli, A.: On the Hamiltonian interpolation of near-to-the-identity symplectic mappings with application to symplectic integration algorithms. J. Stat. Phys. 74(5/6), 1117–1143 (1994) · Zbl 0842.58020 · doi:10.1007/BF02188219
[4] Chierchia, L.: Absolutely continuous spectra of quasi-periodic Schrödinger operators. J. Math. Phys. 28, 2891–2898 (1987) · Zbl 0642.34018 · doi:10.1063/1.527690
[5] Eliasson, L.H.: Floquet solutions for the 1-dimensional quasi-periodic Schrödinger equation. Commun. Math. Phys. 146(3), 447–482 (1992) · Zbl 0753.34055 · doi:10.1007/BF02097013
[6] Hairer, E., Lubich, C.: The life-span of backward error analysis for numerical integrators. Numer. Math. 76, 441–462 (1997) · Zbl 0874.65061 · doi:10.1007/s002110050271
[7] Hairer, E., Lubich, C., Wanner, G.: Geometric Numerical Integration: Structure-preserving Algorithms for Ordinary Differential Equations, 2nd edn., pp. 179–236. Springer, Berlin (2006) · Zbl 1094.65125
[8] Johnson, R.A., Sell, G.R.: Smoothness of spectral subbundles and reducibility of quasi-periodic linear differential equations. J. Differ. Equ. 41, 262–288 (1981) · Zbl 0459.34021 · doi:10.1016/0022-0396(81)90062-0
[9] Jorba, À., Simó, C.: On the reducibility of linear differential equations with quasi-periodic coefficients. J. Differ. Equ. 98(1), 111–124 (1992) · Zbl 0761.34026 · doi:10.1016/0022-0396(92)90107-X
[10] Jorba, À., Simó, C.: On quasi-periodic perturbations of elliptic equilibrium points. SIAM J. Math. Anal. 27(6), 1704–1737 (1996) · Zbl 0863.34043 · doi:10.1137/S0036141094276913
[11] Jorba, À., Villanueva, J.: On the persistence of lower dimensional invariant tori under quasi-periodic perturbations. J. Nonlinear Sci. 7, 427–473 (1997) · Zbl 0898.58044 · doi:10.1007/s003329900036
[12] Jorba, À.: A methodology for the numerical computation of normal forms, centre manifolds, and first integrals of Hamiltonian systems. Exp. Math. 8(2), 155–195 (1999) · Zbl 0983.37107
[13] Kolmogorov, A.N.: On conservation of conditionally periodic motions under small perturbations of the Hamiltonian. Dokl. Akad. Nauk SSSR 98, 527–530 (1954) · Zbl 0056.31502
[14] de la Llave, R.: A tutorial on KAM theory. In: Proceedings of the AMS Summer Research Institute on Smooth Ergodic Theory and its Applications, Seattle, 1999, Proc. Sympos. Pure Math., vol. 69, Amer. Math. Soc., vol. 175292 (2001). Expanded version available at ftp.ma.utexas.edu/pub/papers/llave/tutorial.pdf
[15] MacKay, R.S., Meiss, J.D.: Hamiltonian Dynamical Systems: A Reprint Selection. Adam Hilger, Bristol (1987) · Zbl 0628.01045
[16] Moan, P.C.: On the KAM and Nekhoroschev theorems for symplectic integrators and implications for error growth. Nonlinearity 17, 67–83 (2004) · Zbl 1068.37044 · doi:10.1088/0951-7715/17/1/005
[17] Moan, P.-C.: On rigorous modified equations for discretizations of ODEs. (2005). Available at http://www.focm.net/gi/gips/2005/3.html · Zbl 1094.65068
[18] Moan, P.-C.: On modified equations for discretizations of ODEs. J. Phys. A Math. Gen. 39, 5545–5561 (2006) · Zbl 1094.65068 · doi:10.1088/0305-4470/39/19/S13
[19] Moser, J.: On invariant curves of area-preserving mappings of an annulus. Nachr. Akad. Wiss. Gott. Math.-Phys. Kl. 2, 1–20 (1962) · Zbl 0107.29301
[20] Moser, J., Pöschel, J.: An extension of a result by Dinaburg and Sinai on quasiperiodic potentials. Comment. Math. Helv. 59(1), 39–85 (1984) · Zbl 0533.34023 · doi:10.1007/BF02566337
[21] Palis, J.: Vector fields generate few diffeomorphisms. Bull. Am. Math. Soc. 80(3), 503–505 (1974) · Zbl 0296.57008 · doi:10.1090/S0002-9904-1974-13470-1
[22] Poincaré, H.: Les Méthodes Nouvelles de la Mécanique Céleste, pp. 366–401. NASA English translation of the Dover publication (1957)
[23] Pronin, A.V., Treschev, D.V.: On the inclusion of analytic symplectic maps into analytic flows. Regul. Chaotic Dyn. 2(2), 14–24 (1997) · Zbl 0921.58024
[24] Rychlik, M.: Renormalization of cocycles and linear ODE with almost periodic coefficients. Invent. Math. 110(1), 173–206 (1992) · Zbl 0771.58013 · doi:10.1007/BF01231330
[25] Shang, Z.-J.: KAM theorm of symplectic algorithms for Hamiltonian systems. Numer. Math. 83, 477–496 (1999) · Zbl 0938.65149 · doi:10.1007/s002110050460
[26] Shang, Z.-J.: Resonant and diophantine step sizes in computing invariant tori of Hamiltonian systems. Nonlinearity 13, 299–308 (2000) · Zbl 1008.65092 · doi:10.1088/0951-7715/13/1/314
[27] Shang, Z.-J.: A note on the KAM theorem for symplectic mappings. J. Dyn. Differ. Equ. 12(2), 357–383 (2000) · Zbl 0970.37051 · doi:10.1023/A:1009068425415
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.