Janus and Epimetheus are two moons of Saturn with very peculiar motions. As they orbit around Saturn on quasi-coplanar and quasi-circular trajectories whose radii are only 50 km apart (less than their respective diameters), every four (terrestrial) years the bodies approach each other and their mutual gravitational influence lead to a swapping of the orbits: the outer moon becomes the inner one and vice-versa. This behavior generates horseshoe-shaped trajectories depicted in an appropriate rotating frame. In spite of analytical theories and numerical investigations developed to describe their long-term dynamics, so far very few rigorous long-time stability results on the "horseshoe motion" have been obtained even in the restricted three-body problem. Adapting the idea of Arnol'd (1963) to a resonant case (the co-orbital motion is associated with trajectories in 1:1 mean motion resonance), we provide a rigorous proof of existence of 2-dimensional elliptic invariant tori on which the trajectories are similar to those followed by Janus and Epimetheus. For this purpose, we apply KAM theory to the planar three-body problem.

Popular science article associated with the work `On the co-orbital motion in the three-body problem: existence of quasi-periodic horseshoe-shaped orbits" (arXiv:1806.07262) from the same authors. Janus and Epimetheus are two moons of Saturn which exhibit a really peculiar dynamics. As they orbit on circular trajectories whose radii are only 50 km apart (less than their respective diameters), every four (terrestrial) years the bodies are getting closer and their mutual gravitational influence leads to a swapping of the orbits: the outer moon becoming the inner one and vice-versa. In this article, we describe how, from this specific astronomical problem to the KAM theory, we came to prove the existence of perpetually stable trajectories associated with the Janus and Epimetheus orbits.

We develop a rigorous analytical Hamiltonian formalism adapted to the study of the motion of two planets in co-orbital resonance. By constructing a complex domain of holomorphy for the planetary Hamilto-nian, we estimate the size of the transformation that maps this Hamil-tonian to its first order averaged over one of the fast angles. After having derived an integrable approximation of the averaged problem, we bound the distance between this integrable approximation and the averaged Hamiltonian. This finally allows to prove rigorous theorems on the behavior of co-orbital motions over a finite but large timescale.