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.

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.

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.

This paper describes several applications in astronomy and cosmology that are addressed using probabilistic modelling and statistical inference.

This paper tackles important aspects of comets dynamics from a statistical point of view. Existing methodology uses numerical integration for computing planetary perturbations for simulating such dynamics. This operation is highly computational. It is reasonable to wonder whenever statistical simulation of the perturbations can be much more easy to handle. The first step for answering such a question is to provide a statistical study of these perturbations in order to catch their main features. The statistical tools used are order statistics and heavy tail distributions. The study carried out indicated a general pattern exhibited by the perturbations around the orbits of the important planet. These characteristics were validated through statistical testing and a theoretical study based on Opik theory.

We present Monte Carlo simulations of the dynamical evolution of the Oort cloud over the age of the Solar System, using an initial sample of one million test comets without any cloning. Our model includes perturbations due to the Galactic tide (radial and vertical) and passing stars. We present the first detailed analysis of the injection mechanism into observable orbits by comparing the complete model with separate models for tidal and stellar perturbations alone. We find that a fundamental role for injecting comets from the region outside the loss cone (perihelion distance q > 15 AU) into observable orbits (q < 5 AU) is played by stellar perturbations. These act in synergy with the tide such that the total injection rate is significantly larger than the sum of the two separate rates. This synergy is as important during comet showers as during quiescent periods and concerns comets with both small and large semi-major axes. We propose different dynamical mechanisms to explain the synergies in the inner and outer parts of the Oort Cloud. We find that the filling of the observable part of the loss cone under normal conditions in the present-day Solar System rises from <1 % for a < 20 000 AU to about 100 % for a >= 100 000 AU.