Transit Timing Variations, or TTVs, can be a very efficient way of constraining masses and eccentricities of multi-planet systems. Recent measurements of the TTVs of TRAPPIST-1 led to an estimate of the masses of the planets, enabling an estimate of their densities. A recent TTV analysis using data obtained in the past two years yields a 34% and 13% increase in mass for TRAPPIST-1b and c, respectively. In most studies to date, a Newtonian N-body model is used to fit the masses of the planets, while sometimes general relativity is accounted for. Using the Posidonius N-body code, we show that in the case of the TRAPPIST-1 system, non-Newtonian effects might be also relevant to correctly model the dynamics of the system and the resulting TTVs. In particular, using standard values of the tidal Love number $k_2$ (accounting for the tidal deformation) and the fluid Love number $k_{2f}$ (accounting for the rotational flattening) leads to differences in the TTVs of TRAPPIST-1b and c similar to the differences caused by general relativity. We also show that relaxing the values of tidal Love number $k_2$ and the fluid Love number $k_{2f}$ can lead to TTVs which differ by as much as a few 10~s on a $3-4$-year timescale, which is a potentially observable level. The high values of the Love numbers needed to reach observable levels for the TTVs could be achieved for planets with a liquid ocean, which, if detected, might then be interpreted as a sign that TRAPPIST-1b and TRAPPIST-1c could have a liquid magma ocean. For TRAPPIST-1 and similar systems, the models to fit the TTVs should potentially account for general relativity, for the tidal deformation of the planets, for the rotational deformation of the planets and, to a lesser extent, for the rotational deformation of the star, which would add up to 7x2+1 = 15 additional free parameters in the case of TRAPPIST-1.

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. In this paper we investigate how a pair of co-orbital exoplanets would fare during their migration in a protoplanetary disc. To this end, we computed a stability criterion of the Lagrangian equilibria L4 and L5 under generic dissipation and slow mass evolution. Depending on the strength and shape of these perturbations, the system can either evolve towards the Lagrangian equilibrium, or tend to increase its amplitude of libration, possibly all the way to horseshoe orbits or even exiting the resonance. We estimated the various terms of our criterion using a set of hydrodynamical simulations, and show that the dynamical coupling between the disc perturbations and both planets have a significant impact on the stability: the structures induced by each planet in the disc perturb the dissipative forces applied on the other planets over each libration cycle. Amongst our results on the stability of co-orbitals, several are of interest to constrain the observability of such configurations: long-distance inward migration and smaller leading planets tend to increase the libration amplitude around the Lagrangian equilibria, while leading massive planets and belonging to a resonant chain tend to stabilise it. We also show that, depending on the strength of the dissipative forces, both the inclination and the eccentricity of the smaller of the two co-orbitals can be significantly increased during the inward migration of the co-orbital pair, which can have a significant impact on the detectability by transit of such configurations.

We present a study on the formation of planetary systems around low mass stars similar to Trappist-1, through the accretion of either planetesimals or pebbles. The aim is to determine if the currently observed systems around low mass stars favour one scenario over the other. We ran numerous N-body simulations, coupled to a thermally evolving viscous disc model, including prescriptions for planet migration and photoevaporation. We examine the differences between the pebble and planetesimal accretion scenarios, but also look at the influences of disc mass, planetesimal size, and the percentage of solids locked up within pebbles. When comparing the resulting planetary systems to Trappist-1, we find that a wide range of initial conditions for both accretion scenarios can form planetary systems similar to Trappist-1, in terms of planet mass, periods, and resonant configurations. Typically these planets formed exterior to the water iceline and migrated in resonant convoys to close to the central star. When comparing the planetary systems formed from pebbles to those formed from planetesimals, we find a large number of similarities, including average planet masses, eccentricities, inclinations and period ratios. One major difference was that of the water content of the planets. When including the effects of ablation and full recycling of the planets envelope with the disc, planets formed from pebbles were extremely dry, whilst those formed from planetesimals were extremely wet. If the water content is not fully recycled and instead falls to the planets core, or if ablation of the water is neglected, then the planets formed from pebbles are extremely wet, similar to those formed from planetesimals. Should the water content of the Trappist-1 planets be determined accurately, this could point to a preferred formation pathway for planetary systems, or to specific physics that may be at play.

A low-mass companion to the two-solar mass star HIP65426 has recently been detected by SPHERE at around 100 au from its host. Explaining the presence of super-Jovian planets at large separations, as revealed by direct imaging, is currently an open question. We want to derive statistical constraints on the mass and initial entropy of HIP65426b and to explore possible formation pathways of directly imaged objects within the core-accretion paradigm, focusing on HIP65426b. Constraints on the planet's mass and post-formation entropy are derived from its age and luminosity combined with cooling models. For the first time, the results of population synthesis are also used to inform the results. Then, a formation model that includes N-body dynamics with several embryos per disc is used to study possible formation histories and the properties of possible additional companions. Finally, the outcomes of two- and three-planet scattering in the post-disc phase are analysed, taking tides into account. The mass of HIP65426b is found to be Mp = 9.9 +1.1 -1.8 MJ using the hot population and Mp = 10.9 +1.4 -2.0 MJ with the cold-nominal population. Core formation at small separations from the star followed by outward scattering and runaway accretion at a few hundred AU succeeds in reproducing the mass and separation of HIP65426b. Alternatively, systems having two or more giant planets close enough to be on an unstable orbit at disc dispersal are likely to end up with one planet on a wide HIP65426b-like orbit with a relatively high eccentricity (>~ 0.5). If this scattering scenario explains its formation, HIP65426b is predicted to have a high eccentricity and to be accompanied by one or several roughly Jovian-mass planets at smaller semi-major axes, which also could have a high eccentricity. This could be tested by further direct-imaging as well as radial-velocity observations.

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.

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 Cassini spacecraft revealed the spectacular, highly irregular shapes of the small inner moons of Saturn, ranging from the unique "ravioli-like" forms of Pan and Atlas to the highly elongated structure of Prometheus. Closest to Saturn, these bodies provide important clues regarding the formation process of small moons in close orbits around their host planet, but their range of irregular shapes has not been explained yet. Here we show that the spectrum of shapes among Saturn's small moons is a natural outcome of merging collisions among similar-sized moonlets possessing physical properties and orbits that are consistent with those of the current moons. A significant fraction of such merging collisions take place either at the first encounter or after 1-2 hit-and-run events, with impact velocities in the range of 1-5 times the mutual escape velocity. Close to head-on mergers result in flattened objects with large equatorial ridges, as observed on Atlas and Pan. With slightly more oblique impact angles, collisions lead to elongated, Prometheus-like shapes. These results suggest that the current forms of the small moons provide direct evidence of the processes at the final stages of their formation, involving pairwise encounters of moonlets of comparable size. Finally, we show that this mechanism may also explain the formation of Iapetus' equatorial ridge, as well as its oblate shape.

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.

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.

This work focuses on the dynamics and the detection methods of co-orbital exoplanets. We call "co-orbital" any configuration in which two planets orbit with the same mean mean-motion around the same star. First, we revisit the results of the circular coplanar case. We also recall that the manifold associated to the coplanar case and the manifold corresponding to the circular case are both invariant by the flow of the averaged Hamiltonian. We hence study these two particular cases. We focus mainly on the coplanar case (eccentric), where we study the evolution of families of non-maximal quasi-periodic orbits parametrized by the eccentricity of the planets. We show that the geometry of these families is highly dependent on the eccentricity, which causes significant topology changes across the space of phases as the latter increases. A chapter is dedicated to the detection of co-orbital exoplanets. We recall the different detection methods adapted to the co-orbital case. We focus on the radial velocity technique, and the combination of radial velocity and transit measurements. Finally, we describe a method to study the effect of orbital perturbations on the spin-orbit resonances for a rigid body. We apply this method in two cases: the eccentric co-orbital case and the circumbinary case.

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 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.

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.