Context. The strong tidal dissipation in the couple Jupiter-Io is spread to all the moons involved in the Laplace resonance (Io, Europa, and Ganymede), leading to a migration of their orbits. Aims. We aim to characterize the future behavior of the Galilean satellites over the Solar System lifetime and to quantify the stability of the Laplace resonance. Since tidal dissipation makes possible the exit from the current resonances or capture into new ones, we investigate the capture of Callisto into resonance. Methods. We perform hundreds of propagations using an improved version of a recent semi-analytical model. As Ganymede moves outwards, it approaches the 2:1 resonance with Callisto, inducing a temporary chaotic motion in the system. For this reason, we draw a statistical picture of the outcome of the resonant encounter. Results. The system can settle into two distinct outcomes: A) a chain of three 2:1 two-body resonances (Io-Europa, Europa-Ganymede and Ganymede-Callisto), or B) a resonant chain involving the 2:1 two-body resonance Io-Europa plus at least one pure 4:2:1 three-body resonance, most frequently between Europa, Ganymede and Callisto. In case A (56\% of the simulations), the Laplace resonance is always preserved and the eccentricities remain confined to small values below 0.01. In case B (44\% of the simulations), the Laplace resonance is generally disrupted and the eccentricities of Ganymede and Callisto can increase up to about 0.1, making this configuration unstable and driving the system into new resonances. Conclusion. From our results, the capture of Callisto into resonance appears to be extremely likely (100\% of our simulations). Assuming the most recent estimate of the dissipation between Io and Jupiter, the resonant encounter happens at about 1.5 Gyrs from now. Therefore, the stability of the Laplace resonance is guaranteed at least up to about 1.5 Gyrs.

This article reviews the different mechanisms affecting the orbits of trans-Neptunian objects, ranging from internal perturbations (planetary scattering, mean-motion resonances, secular effects) to external perturbations (galactic tides, passing stars). We outline the theoretical tools that can be used to model and study them, focussing on analytical approaches. We eventually compare these mechanisms to the observed distinct populations of trans-Neptunian objects and conclude on how they participate to the sculpting of the whole distribution.

Context: Distant trans-Neptunian objects are subject to planetary perturbations and galactic tides. The former decrease with the distance, while the latter increase. In the intermediate regime where they have the same order of magnitude (the 'inert Oort cloud'), both are weak, resulting in very long evolution timescales. To date, three observed objects can be considered to belong to this category. Aims: We aim to provide a clear understanding of where this transition occurs, and to characterise the long-term dynamics of small bodies in the intermediate regime: relevant resonances, chaotic zones (if any), and timescales at play. Results: There exists a tilted equilibrium plane (Laplace plane) about which orbits precess. The dynamics is integrable in the low and high semi-major axis regimes, but mostly chaotic in between. From 800 to 1100 au, the chaos covers almost all the eccentricity range. The diffusion timescales are large, but not to the point of being indiscernible in a 4.5 Gyrs duration: the perihelion distance can actually vary from tens to hundreds of au. Orbital variations are favoured in specific ranges of inclination corresponding to well-defined resonances. Starting from uniform distributions, the orbital angles cluster after 4.5 Gyrs for semi-major axes larger than 500 au, because of a very slow differential precession. Conclusions: Even if it is characterised by very long timescales, the inert Oort cloud is much less inert than it appears. Orbits can be considered inert over 4.5 Gyrs only in small portions of the space of orbital elements, which include (90377) Sedna and 2012VP113. Effects of the galactic tides are discernible down to semi-major axes of about 500 au. We advocate including the galactic tides in simulations of distant trans-Neptunian objects, especially when studying the formation of detached bodies or the clustering of orbital elements.

Context: Seasonal variations and climate stability of a planet are very sensitive to the planet obliquity and its evolution. This is of particular interest for the emergence and sustainability of land-based life, but orbital and rotational parameters of exoplanets are still poorly constrained. Numerical explorations usually realised in this situation are thus in heavy contrast with the uncertain nature of the available data. Aims: We aim to provide an analytical formulation of the long-term spin-axis dynamics of exoplanets, linking it directly to physical and dynamical parameters, but still giving precise quantitative results if the parameters are well known. Together with bounds for the poorly constrained parameters of exoplanets, this analysis is designed to allow a quick and straightforward exploration of the spin-axis dynamics. Methods: The long-term orbital solution is decomposed in quasi-periodic series and the spin-axis Hamiltonian is expanded in powers of eccentricity and inclination. Chaotic zones are measured by the resonance overlap criterion. Bounds for the poorly known parameters of exoplanets are obtained from physical grounds (rotational breakup) and dynamical considerations (equipartition of AMD). Results: This method gives accurate results when the orbital evolution is well known. The chaotic zones for planets of the Solar System can be retrieved in details from simple analytical formulas. For less constrained planetary systems, the maximal extent of the chaotic regions can be computed, requiring only the mass, the semi-major axis and the eccentricity of the planets present in the system. Additionally, some estimated bounds of the precession constant allow to classify which observed exoplanets are necessarily out of major spin-orbit secular resonances (unless the precession rate is affected by the presence of massive satellites).

Natural satellite systems present a large variety of orbital configurations in the solar system. While some are clearly the result of known processes, others still have largely unexplained eccentricity and inclination values. Iapetus has a still unexplained 3% orbital eccentricity and its orbital plane is tilted with respect to its local Laplace plane. On the other hand, astrometric measurements of saturnian moons have revealed high tidal migration rates, corresponding to a quality factor Q of Saturn of around 1600 for the mid-sized icy moons. We show how a past crossing of the 5:1 mean motion resonance between Titan and Iapetus may be a plausible scenario to explain Iapetus' orbit. We have carried out numerical simulations of the resonance crossing using an N-Body code as well as using averaged equations of motion. A large span of migration rates were explored for Titan and Iapetus was started on its local Laplace plane with a circular orbit. The resonance crossing can trigger a chaotic evolution of the eccentricity and the inclination of Iapetus. The outcome of the resonance is highly dependent on the migration rate. For a quality factor Q of over around 2000, the chaotic evolution of Iapetus in the resonance leads in most cases to its ejection, while simulations with a quality factor between 100 and 2000 show a departure from the resonance with post-resonant eccentricities spanning from 0 up to 15%, and free inclinations capable of reaching 11 degrees. Usually high inclinations come with high eccentricities but some simulations show elements compatible with Iapetus' current orbit. A quality factor between 100 and 2000 at the frequency of Titan would bring Titan and Iapetus into a 5:1 resonance, which would perturb Iapetus' eccentricity and inclination to values observed today. Such rapid tidal migration would have avoided Iapetus' ejection around 40 to 800 million years ago.

Aims. Observations of solar protons near comet 67P/Churyumov-Gerasimenko (67P) by the Rosetta spacecraft can be modelled by the planar motion in an effective magnetic field proportional to 1/r^2. We aim to provide a thorough study of such dynamics, with a clear description of the behaviour of an incoming flux of particles. We will be able, then, to calibrate the free parameters of the model to Rosetta observations. Methods. Basic tools of dynamical analysis are used. They lead to a definition of the relevant parameters for the system and a classification of the possible types of trajectories. Using the so-obtained formalism, the structures formed by a flux of particles coming from infinity can be studied. Results. All the trajectories are parametrised by two characteristic radii, rE and rC, derived from first integrals. There are three different types of motion possible divided by a separatrix corresponding to rE = rC. An analytical expression of the trajectories, defined by an integral, is developed. Using this formalism, the application to a flux of particles coming from infinity (modelling the incident solar wind) gives one free parameter only, the radius rE, which scales the problem. A circular cavity of radius 0.28 rE is created, as well as an overdensity curve (analogous to a caustic in optics). At each observation time, rE can be calibrated to Rosetta plasma measurements, giving a qualitative understanding of the solar particle dynamics (incoming direction, cavity and density map). We also deduce that, in order to properly capture the essence of the dynamics, numerical simulations of the solar wind around a comet must use simulation boxes much larger than rE and grids much finer than rE.

We aim at analytically modelling the solar wind proton trajectories during their interaction with a partially ionised cometary atmosphere, not in terms of bulk properties of the flow but in terms of single particle dynamics. We first derive a generalised gyromotion, in which the electric field is reduced to its motional component. Steady-state is assumed, and simplified models of the cometary density and of the electron fluid are used to express the force experienced by individual solar wind protons during the interaction. A three-dimensional (3D) analytical expression of the gyration of two interacting plasma beams is obtained. Applying it to a comet case, the force on protons is always perpendicular to their velocity and has an amplitude proportional to 1/r^2. The solar wind deflection is obtained at any point in space. The resulting picture presents a caustic of intersecting trajectories, and a circular region is found that is completely free of particles. The particles do not lose any kinetic energy and this absence of deceleration, together with the solar wind deflection pattern and the presence of a solar wind ion cavity, is in good agreement with the general results of the Rosetta mission. The qualitative match between the model and the in situ data highlights how dominant the motional electric field is throughout most of the interaction region for the solar wind proton dynamics. The model provides a simple general kinetic description of how momentum is transferred between these two collisionless plasmas. It also shows the potential of this semi-analytical model for a systematic quantitative comparison to the data.

The first 1000 km of the ion tail of comet 67P/Churyumov-Gerasimenko were explored by the European Rosetta spacecraft, 2.7 au away from the Sun. We characterised the dynamics of both the solar wind and the cometary ions on the night-side of the comet's atmosphere. We analysed in situ ion and magnetic field measurements and compared the data to a semi-analytical model. The cometary ions are observed flowing close to radially away from the nucleus during the entire excursion. The solar wind is deflected by its interaction with the new-born cometary ions. Two concentric regions appear, an inner region dominated by the expanding cometary ions and an outer region dominated by the solar wind particles. The single night-side excursion operated by Rosetta revealed that the near radial flow of the cometary ions can be explained by the combined action of three different electric field components, resulting from the ion motion, the electron pressure gradients, and the magnetic field draping. The observed solar wind deflection is governed mostly by the motional electric field.

We use a secular model to describe the non-resonant dynamics of trans-Neptunian objects in the presence of an external ten-earth-mass perturber. The secular dynamics is analogous to an "eccentric Kozai mechanism" but with both an inner component (the four giant planets) and an outer one (the eccentric distant perturber). By the means of Poincar\'e sections, the cases of a non-inclined or inclined outer planet are successively studied, making the connection with previous works. In the inclined case, the problem is reduced to two degrees of freedom by assuming a non-precessing argument of perihelion for the perturbing body. The size of the perturbation is typically ruled by the semi-major axis of the small body: we show that the classic integrable picture is still valid below about 70 AU, but it is progressively destroyed when we get closer to the external perturber. In particular, for a>150 AU, large-amplitude orbital flips become possible, and for a>200 AU, the Kozai libration islands are totally submerged by the chaotic sea. Numerous resonance relations are highlighted. The most large and persistent ones are associated to apsidal alignments or anti-alignments with the orbit of the distant perturber.

Aims. Numerous trans-Neptunian objects are known to be in mean-motion resonance with Neptune. We aim to describe their long-term orbital evolution (both past and future) by means of a one-degree-of-freedom secular model. In this paper, we focus only on objects with a semi-major axis larger than 50 astronomical units (au). Methods. For each resonant object considered, a 500 000-year numerical integration is performed. The output is digitally filtered to get the parameters of the resonant secular model. Their long-term (Giga-year) orbital evolution is then represented by the level curves of the secular Hamiltonian. Results. For the majority of objects considered, the mean-motion resonance has little impact on the long-term trajectories (the secular dynamics is similar to a non-resonant one). However, a subset of objects is strongly affected by the resonance, producing moderately-high-amplitude oscillations of the perihelion distance and/or libration of the argument of perihelion around a fixed centre. Moreover, the high perihelion distance of the object 2015 FJ345 is plainly explained by long-term resonant dynamics, allowing us to also deduce its orbital elements at the time of capture in resonance (at least 15 million years ago). The same type of past evolution is expected for 2014 FZ71.

Two semi-analytical one-degree-of-freedom secular models are presented for the motion of small bodies beyond Neptune. A special attention is given to trajectories entirely exterior to the planetary orbits. The first one is the well-known non-resonant model of Kozai (1962) adapted to the transneptunian region. Contrary to previous papers, the dynamics is fully characterized with respect to the fixed parameters. A maximum perihelion excursion possible of 16.4 AU is determined. The second model handles the occurrence of a mean-motion resonance with one of the planets. In that case, the one-degree-of-freedom integrable approximation is obtained by postulating the adiabatic invariance, and is much more general and accurate than previous secular models found in the literature. It brings out in a plain way the possibility of perihelion oscillations with a very high amplitude. Such a model could thus be used in future studies to deeper explore that kind of motion. For complex resonant orbits (especially of type 1:k), a segmented secular description is necessary since the trajectories are only "integrable by parts". The two models are applied to the Solar System but the notations are kept general so that it could be used for any quasi-circular and coplanar planetary system.

We use a secular representation to describe the long-term dynamics of transneptunian objects in mean-motion resonance with Neptune. The model applied is thoroughly described in Saillenfest et al. (2016). The parameter space is systematically explored, showing that the secular trajectories depend little on the resonance order. High-amplitude oscillations of the perihelion distance are reported and localised in the space of the orbital parameters. In particular, we show that a large perihelion distance is not a sufficient criterion to declare that an object is detached from the planets. Such a mechanism, though, is found unable to explain the orbits of Sedna or 2012VP113, which are insufficiently inclined (considering their high perihelion distance) to be possibly driven by such a resonant dynamics. The secular representation highlights the existence of a high-perihelion accumulation zone due to resonances of type 1:k with Neptune. That region is found to be located roughly at semi-major axes in [100;300] AU, perihelion distances in [50;70] AU and inclinations [30;50]{\deg}. In addition to the flux of objects directly coming from the Scattered Disc, numerical simulations show that the Oort Cloud is also a substantial source for such objects. Naturally, as that mechanism relies on fragile captures in high-order resonances, our conclusions break down in the case of a significant external perturber. The detection of such a reservoir could thus be an observational constraint to probe the external Solar System.

We revisit the dynamics of Atlas. Using Cassini ISS astrometric observations spanning February 2004 to August 2013, Cooper et al. (2015) found evidence that Atlas is currently perturbed by both a 54:53 corotation eccentricity resonance (CER) and a 54:53 Lindblad eccentricity resonance (LER) with Prometheus. They demonstrated that the orbit of Atlas is chaotic, with a Lyapunov time of order 10 years, as a direct consequence of the coupled resonant interaction (CER/LER) with Prometheus. Here we investigate the interactions between the two resonances using the CoraLin analytical model (El Moutamid et al. 2014), showing that the chaotic zone fills almost all the corotation sites occupied by the satellite's orbit. Four 70:67 apse-type mean motion resonances with Pandora are also overlapping, but these resonances have a much weaker effect. Frequency analysis allows us to highlight the coupling between the 54:53 resonances, and confirms that a simplified system including the perturbations due to Prometheus and Saturn's oblateness only captures the essential features of the dynamics.