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.

Despite the existence of co-orbital bodies in the solar system, and the prediction of the formation of co-orbital planets by planetary system formation models, no co-orbital exoplanets (also called trojans) have been detected thus far. Here we study the signature of co-orbital exoplanets in transit surveys when two planet candidates in the system orbit the star with similar periods. Such pair of candidates could be discarded as false positives because they are not Hill-stable. However, horseshoe or long libration period tadpole co-orbital configurations can explain such period similarity. This degeneracy can be solved by considering the Transit Timing Variations (TTVs) of each planet. We then focus on the three planet candidates system TOI-178: the two outer candidates of that system have similar orbital period and had an angular separation near $\pi/3$ during the TESS observation of sector 2. Based on the announced orbits, the long-term stability of the system requires the two close-period planets to be co-orbitals. Our independent detrending and transit search recover and slightly favour the three orbits close to a 3:2:2 resonant chain found by the TESS pipeline, although we cannot exclude an alias that would put the system close to a 4:3:2 configuration. We then analyse in more detail the co-orbital scenario. We show that despite the influence of an inner planet just outside the 2:3 mean-motion resonance, this potential co-orbital system can be stable on the Giga-year time-scale for a variety of planetary masses, either on a trojan or a horseshoe orbit. We predict that large TTVs should arise in such configuration with a period of several hundred days. We then show how the mass of each planet can be retrieved from these TTVs.

In the framework of the planar and circular restricted three-body problem, we consider an asteroid that orbits the Sun in quasi-satellite motion with a planet. A quasi-satellite trajectory is a heliocentric orbit in co-orbital resonance with the planet, characterized by a non zero eccentricity and a resonant angle that librates around zero. Likewise, in the rotating frame with the planet it describes the same trajectory as the one of a retrograde satellite even though the planet acts as a perturbator. In the last few years, the discoveries of asteroids in this type of motion made the term "quasi-satellite" more and more present in the literature. However, some authors rather use the term "retrograde satellite" when referring to this kind of motion in the studies of the restricted problem in the rotating frame. In this paper we intend to clarify the terminology to use, in order to bridge the gap between the perturbative co-orbital point of view and the more general approach in the rotating frame. Through a numerical exploration of the co-orbital phase space, we describe the quasi-satellite domain and highlight that it is not reachable by low eccentricities by averaging process. We will show that the quasi-satellite domain is effectively included in the domain of the retrograde satellites and neatly defined in terms of frequencies. Eventually, we highlight a remarkable high eccentric quasi-satellite orbit corresponding to a frozen ellipse in the heliocentric frame. We extend this result to the eccentric case (planet on an eccentric motion) and show that two families of frozen ellipses originate from this remarkable orbit.

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.

Co-orbital bodies are the byproduct of planet formation and evolution, as we know from the Solar System. Although planet-size co-orbitals do not exists in our planetary system, dynamical studies show that they can remain stable for long periods of time in the gravitational well of massive planets. Should they exist, their detection is feasible with the current instrumentation. In this paper, we present new ground-based observations searching for these bodies co-orbiting with nine close-in (P<5 days) planets, using different observing techniques. The combination of all of them allows us to restrict the parameter space of any possible trojan in the system. We use multi-technique observations (radial velocity, precision photometry and transit timing variations), both newly acquired in the context of the TROY project and publicly available, to constrain the presence of planet-size trojans in the Lagrangian points of nine known exoplanets. We find no clear evidence of trojans in these nine systems through any of the techniques used down to the precision of the observations. However, this allows us to constrain the presence of any potential trojan in the system, specially in the trojan mass/radius versus libration amplitude plane. In particular, we can set upper mass limits in the super-Earth mass regime for six of the studied systems.

The search for extrasolar planets in the past decades has shown that planets abound in the Solar neighborhood. While we are still missing an Earth twin, the forthcoming space missions and ground-based instrumentation are already driven to achieve this goal. But, in order to fully understand the conditions for life appearing in the Solar System, we still miss some pieces of the planetary system jigsaw puzzle, namely a deeper understanding of the minor bodies. Trojans, moons, and comets are tracers of the formation and evolution processes of planetary systems. These missing pieces are also critical to understand the emergence and evolution of life over millions of years. With the large crop of planetary systems discovered so far and yet to be detected with the forthcoming missions, the hunt for minor bodies in extrasolar systems is a natural continuation of our search for real Solar System- and, in particular, Earth- analogs. This white paper is focused on detection of these minor components and their relevance in the emergence, evolution and survival of life.

We study the phase space of eccentric coplanar co-orbitals in the non-restricted case. Departing from the quasi-circular case, we describe the evolution of the phase space as the eccentricities increase. We find that over a given value of the eccentricity, around $0.5$ for equal mass co-orbitals, important topological changes occur in the phase space. These changes lead to the emergence of new co-orbital configurations and open a continuous path between the previously distinct trojan domains near the $L_4$ and $L_5$ eccentric Lagrangian equilibria. These topological changes are shown to be linked with the reconnection of families of quasi-periodic orbits of non-maximal dimension.

We present an analytical formalism to study the secular dynamics of a system consisting of N-2 planets orbiting a binary star in outer orbits. We introduce a canonical coordinate system and expand the disturbing function in terms of canonical elliptic elements, combining both Legendre polynomials and Laplace coefficients, to obtain a general formalism for the secular description of this type of configuration. With a quadractic approximation of the development, we present a simplified analytical solution for the planetary orbits for both the single planet and the two-planet cases. From the two-planet model, we show that the inner planet accelerates the precession rate of the binary pericenter, which, in turn, may enter in resonance with the secular frequency of the outer planet, characterizing a secular resonance. We calculate an analytical expression for the approximate location of this resonance and apply it to known circumbinary systems, where we show that it can occur at relatively close orbits, for example at 2.4 au for the Kepler-38 system. With a more refined model, we analyse the dynamics of this secular resonance and we show that a bifurcation of the corresponding fixed points can affect the long term evolution and stability of planetary systems. By comparing our results with complete integrations of the exact equations of motion, we verified the accuracy of our analytical model.

The detection of Earth-like planets, exocomets or Kuiper belts show that the different components found in the solar system should also be present in other planetary systems. Trojans are one of these components and can be considered fossils of the first stages in the life of planetary systems. Their detection in extrasolar systems would open a new scientific window to investigate formation and migration processes. In this context, the main goal of the TROY project is to detect exotrojans for the first time and to measure their occurrence rate (eta-Trojan). In this first paper, we describe the goals and methodology of the project. Additionally, we used archival radial velocity data of 46 planetary systems to place upper limits on the mass of possible trojans and investigate the presence of co-orbital planets down to several tens of Earth masses. We used archival radial velocity data of 46 close-in (P<5 days) transiting planets (without detected companions) with information from high-precision radial velocity instruments. We took advantage of the time of mid-transit and secondary eclipses (when available) to constrain the possible presence of additional objects co-orbiting the star along with the planet. This, together with a good phase coverage, breaks the degeneracy between a trojan planet signature and signals coming from additional planets or underestimated eccentricity. We identify nine systems for which the archival data provide 1-sigma evidence for a mass imbalance between L4 and L5. Two of these systems provide 2-sigma detection, but no significant detection is found among our sample. We also report upper limits to the masses at L4/L5 in all studied systems and discuss the results in the context of previous findings.

We study the spin evolution of close-in planets in compact multi-planetary systems. The rotation period of these planets is often assumed to be synchronous with the orbital period due to tidal dissipation. Here we show that planet-planet perturbations can drive the spin of these planets into non-synchronous or even chaotic states. In particular, we show that the transit timing variation (TTV) is a very good probe to study the spin dynamics, since both are dominated by the perturbations of the mean longitude of the planet. We apply our model to KOI-227b and Kepler-88b, which are both observed undergoing strong TTVs. We also perform numerical simulations of the spin evolution of these two planets. We show that for KOI-227b non-synchronous rotation is possible, while for Kepler-88b the rotation can be chaotic.

Co-orbital planets have not yet been discovered, although they constitute a frequent by-product of planetary formation and evolution models. This lack may be due to observational biases, since the main detection methods are unable to spot co-orbital companions when they are small or near the Lagrangian equilibrium points. However, for a system with one known transiting planet (with mass $m_1$), we can detect a co-orbital companion (with mass $m_2$) by combining the time of mid-transit with the radial-velocity data of the star. Here, we propose a simple method that allows the detection of co-orbital companions, valid for eccentric orbits, that relies on a single parameter $\alpha$, which is proportional to the mass ratio $m_2/m_1$. Therefore, when $\alpha$ is statistically different from zero, we have a strong candidate to harbour a co-orbital companion. We also discuss the relevance of false positives generated by different planetary configurations.

Several celestial bodies in co-orbital configurations exist in the solar system. However, co-orbital exoplanets have not yet been discovered. This lack may result from a degeneracy between the signal induced by co-orbital planets and other orbital configurations. Here we determine a criterion for the detectability of quasi-circular co-orbital planets and develop a demodulation method to bring out their signature from the observational data. We show that the precision required to identify a pair of co-orbital planets depends only on the libration amplitude and on the planet's mass ratio. We apply our method to synthetic radial velocity data, and show that for tadpole orbits we are able to determine the inclination of the system to the line of sight. Our method is also valid for planets detected through the transit and astrometry techniques.

We investigate the resonant rotation of co-orbital bodies in eccentric and planar orbits. We develop a simple analytical model to study the impact of the eccentricity and orbital perturbations on the spin dynamics. This model is relevant in the entire domain of horseshoe and tadpole orbit, for moderate eccentricities. We show that there are three different families of spin-orbit resonances, one depending on the eccentricity, one depending on the orbital libration frequency, and another depending on the pericenter's dynamics. We can estimate the width and the location of the different resonant islands in the phase space, predicting which are the more likely to capture the spin of the rotating body. In some regions of the phase space the resonant islands may overlap, giving rise to chaotic rotation.

We analyse the secular dynamics of planets on S-type coplanar orbits in tight binary systems, based on first- and second-order analytical models, and compare their predictions with full N-body simulations. The perturbation parameter adopted for the development of these models depends on the masses of the stars and on the semimajor axis ratio between the planet and the binary. We show that each model has both advantages and limitations. While the first-order analytical model is algebraically simple and easy to implement, it is only applicable in regions of the parameter space where the perturbations are sufficiently small. The second-order model, although more complex, has a larger range of validity and must be taken into account for dynamical studies of some real exoplanetary systems such as $\gamma$-Cephei and HD 41004A. However, in some extreme cases, neither of these analytical models yields quantitatively correct results, requiring either higher-order theories or direct numerical simulations. Finally, we determine the limits of applicability of each analytical model in the parameter space of the system, giving an important visual aid to decode which secular theory should be adopted for any given planetary system in a close binary.

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.

We investigate the resonant rotation of circumbinary bodies in planar quasi-circular orbits. Denoting $n_b$ and $n$ the orbital mean motion of the inner binary and of the circumbinary body, respectively, we show that spin-orbit resonances exist at the frequencies $n\pm k\nu/2$, where $\nu = n_b - n$, and $k$ is an integer. Moreover, when the libration at natural frequency has the same magnitude as $\nu$, the resonances overlap and the rotation becomes chaotic. We apply these results to the small satellites in the Pluto-Charon system, and conclude that their rotations are likely chaotic. However, the rotation can also be stable and not synchronous for small axial asymmetries.

In the circular case of the coplanar Restricted Three-body Problem, we studied how the family of quasi-satellite (QS) periodic orbits allows to define an associated libration center. Using the averaged problem, we highlighted a validity limit of this one: for QS orbits with low eccentricities, the averaged problem does not correspond to the real problem. We do the same procedure to L 3 , L 4 and L 5 emerging periodic orbits families and remarked that for very high eccentricities F L4 and F L5 merge with F L3 which bifurcates to a stable family.

The presence of a co-orbital companion induces the splitting of the well known Keplerian spin-orbit resonances. It leads to chaotic rotation when those resonances overlap.

The rotation of asymmetric bodies in eccentric Keplerian orbits can be chaotic when there is some overlap of spin-orbit resonances. Here we show that the rotation of two coorbital bodies (two planets orbiting a star or two satellites of a planet) can also be chaotic even for quasi-circular orbits around the central body. When dissipation is present, the rotation period of a body on a nearly circular orbit is believed to always end synchronous with the orbital period. Here we demonstrate that for coorbital bodies in quasi-circular orbits, stable non-synchronous rotation is possible for a wide range of mass ratios and body shapes. We further show that the rotation becomes chaotic when the natural rotational libration frequency, due to the axial asymmetry, is of the same order of magnitude as the orbital libration frequency.

Coorbital bodies are observed around the Sun sharing their orbits with the planets, but also in some pairs of satellites around Saturn. The existence of coorbital planets around other stars has also been proposed. For close-in planets and satellites, the rotation slowly evolves due to dissipative tidal effects until some kind of equilibrium is reached. When the orbits are nearly circular, the rotation period is believed to always end synchronous with the orbital period. Here we demonstrate that for coorbital bodies in quasi-circular orbits, stable non-synchronous rotation is possible for a wide range of mass ratios and body shapes. We show the existence of an entirely new family of spin-orbit resonances at the frequencies $n\pm k\nu/2$, where $n$ is the orbital mean motion, $\nu$ the orbital libration frequency, and $k$ an integer. In addition, when the natural rotational libration frequency due to the axial asymmetry, $\sigma$, has the same magnitude as $\nu$, the rotation becomes chaotic. Saturn coorbital satellites are synchronous since $\nu\ll\sigma$, but coorbital exoplanets may present non-synchronous or chaotic rotation. Our results prove that the spin dynamics of a body cannot be dissociated from its orbital environment. We further anticipate that a similar mechanism may affect the rotation of bodies in any mean-motion resonance.

We develop an analytical Hamiltonian formalism adapted to the study of the motion of two planets in co-orbital resonance. The Hamiltonian, averaged over one of the planetary mean longitude, is expanded in power series of eccentricities and inclinations. The model, which is valid in the entire co-orbital region, possesses an integrable approximation modeling the planar and quasi-circular motions. First, focusing on the fixed points of this approximation, we highlight relations linking the eigenvectors of the associated linearized differential system and the existence of certain remarkable orbits like the elliptic Eulerian Lagrangian configurations, the Anti-Lagrange (Giuppone et al., 2010) orbits and some second sort orbits discovered by Poincar\'e. Then, the variational equation is studied in the vicinity of any quasi-circular periodic solution. The fundamental frequencies of the trajectory are deduced and possible occurrence of low order resonances are discussed. Finally, with the help of the construction of a Birkhoff normal form, we prove that the elliptic Lagrangian equilateral configurations and the Anti-Lagrange orbits bifurcate from the same fixed point L4.

The Cassini spacecraft collects high resolution images of the saturnian satellites and reveals the surface of these new worlds. The shape and rotation of the satellites can be determined from the Cassini Imaging Science Subsystem data, employing limb coordinates and stereogrammetric control points. This is the case for Epimetheus (Tiscareno et al. 2009) that opens elaboration of new rotational models (Tiscareno et al. 2009; Noyelles 2010; Robutel et al. 2011). Especially, Epimetheus is characterized by its horseshoe shape orbit and the presence of the swap is essential to introduce explicitly into rotational models. During its journey in the saturnian system, Cassini spacecraft accumulates the observational data of the other satellites and it will be possible to determine the rotational parameters of several of them. To prepare these future observations, we built rotational models of the coorbital (also called Trojan) satellites Telesto, Calypso, Helene, and Polydeuces, in addition to Janus and Epimetheus. Indeed, Telesto and Calypso orbit around the L_4 and L_5 Lagrange points of Saturn-Tethys while Helene and Polydeuces are coorbital of Dione. The goal of this study is to understand how the departure from the Keplerian motion induced by the perturbations of the coorbital body, influences the rotation of these satellites. To this aim, we introduce explicitly the perturbation in the rotational equations by using the formalism developed by Erdi (1977) to represent the coorbital motions, and so we describe the rotational motion of the coorbitals, Janus and Epimetheus included, in compact form.

Janus and Epimetheus are famously known for their distinctive horseshoe-shaped orbits resulting from a 1:1 orbital resonance. Every four years these two satellites swap their orbits by a few tens of kilometers as a result of their close encounter. Recently Tiscareno et al. (2009) have proposed a model of rotation based on images from the Cassini orbiter. These authors inferred the amplitude of rotational librational motion in longitude at the orbital period by fitting a shape model to the recent Cassini ISS images. By a quasiperiodic approximation of the orbital motion, we describe how the orbital swap impacts the rotation of the satellites. To that purpose, we have developed a formalism based on quasi-periodic series with long and short-period librations. In this framework, the amplitude of the libration at the orbital period is found proportional to a term accounting for the orbital swap. We checked the analytical quasi-periodic development by performing a numerical simulation and find both results in good agreement. To complete this study, the results regarding the short-period librations are studied with the help of an adiabatic-like approach.

We hereby study the stability of a massless probe orbiting around an oblate central body (planet or planetary satellite) perturbed by a third body, assumed to lie in the equatorial plane (Sun or Jupiter for example) using an Hamiltonian formalism. We are able to determine, in the parameters space, the location of the frozen orbits, namely orbits whose orbital elements remain constant on average, to characterize their stability/unstability and to compute the periods of the equilibria. The proposed theory is general enough, to be applied to a wide range of probes around planet or natural planetary satellites. The BepiColombo mission is used to motivate our analysis and to provide specific numerical data to check our analytical results. Finally, we also bring to the light that the coefficient J_2 is able to protect against the increasing of the eccentricity due to the Kozai-Lidov effect.

In a previous paper, we have found that the resonance structure of the present Jupiter Trojan swarms could be split up into four different families of resonances. Here, in a first step, we generalize these families in order to describe the resonances occurring in Trojan swarms embedded in a generic planetary system. The location of these families changes under a modification of the fundamental frequencies of the planets and we show how the resonant structure would evolve during a planetary migration. We present a general method, based on the knowledge of the fundamental frequencies of the planets and on those that can be reached by the Trojans, which makes it possible to predict and localize the main events arising in the swarms during migration. In particular, we show how the size and stability of the Trojan swarms are affected by the modification of the frequencies of the planets. Finally, we use this method to study the global dynamics of the Jovian Trojan swarms when Saturn migrates outwards. Besides the two resonances found by Morbidelli et al (2005) which could have led to the capture of the current population just after the crossing of the 2:1 orbital resonance, we also point out several sequences of chaotic events that can influence the Trojan population.

We present a class of symplectic integrators adapted for the integration of perturbed Hamiltonian systems of the form $H=A+\epsilon B$. We give a constructive proof that for all integer $p$, there exists an integrator with positive steps with a remainder of order $O(\tau^p\epsilon +\tau^2\epsilon^2)$, where $\tau$ is the stepsize of the integrator. The analytical expressions of the leading terms of the remainders are given at all orders. In many cases, a corrector step can be performed such that the remainder becomes $O(\tau^p\epsilon +\tau^4\epsilon^2)$. The performances of these integrators are compared for the simple pendulum and the planetary 3-Body problem of Sun-Jupiter-Saturn.