Many extrasolar systems possessing planets in mean-motion resonance or resonant chain have been discovered to date. The transit method coupled with transit timing variation analysis provides an insight into the physical and orbital parameters of the systems, but suffers from observational limitations. When a (near-)resonant planetary system resides in the dynamical neighbourhood of a stable periodic orbit, its long-term stability, and thus survival, can be guaranteed. We use the intrinsic property of the periodic orbits, namely their linear horizontal and vertical stability, to validate or further constrain the orbital elements of detected two-planet systems. We computed the families of periodic orbits in the general three-body problem for several two-planet Kepler and K2 systems. The dynamical neighbourhood of the systems is unveiled with maps of dynamical stability. Additional validations or constraints on the orbital elements of K2-21, K2-24, Kepler-9, and (non-coplanar) Kepler-108 near-resonant systems were achieved. While a mean-motion resonance locking protects the long-term evolution of the systems K2-21 and K2-24, such a resonant evolution is not possible for the Kepler-9 system, whose stability is maintained through an apsidal anti-alignment. For the Kepler-108 system, we find that the stability of its mutually inclined planets could be justified either solely by a mean-motion resonance, or in tandem with an inclination-type resonance. Going forward, dynamical analyses based on periodic orbits could yield better constrained orbital elements of near-resonant extrasolar systems when performed in parallel to the fitting of the observational data.

Extrasolar systems with planets on eccentric orbits close to or in mean-motion resonances are common. The classical low-order resonant Hamiltonian expansion is unfit to describe the long-term evolution of these systems. We extend the Laplace-Lagrange secular approximation for coplanar systems with two planets by including (near-)resonant harmonics, and realize an expansion at high order in the eccentricities of the resonant Hamiltonian both at orders one and two in the masses. We show that the expansion at first order in the masses gives a qualitative good approximation of the dynamics of resonant extrasolar systems with moderate eccentricities, while the second order is needed to reproduce more accurately their orbital evolutions. The resonant approach is also required to correct the secular frequencies of the motion given by the Laplace-Lagrange secular theory in the vicinity of a mean-motion resonance. The dynamical evolutions of four (near-)resonant extrasolar systems are discussed, namely GJ 876 (2:1 resonance), HD 60532 (3:1), HD 108874 and GJ 3293 (close to 4:1).

Aims. To date, more than 600 multi-planetary systems have been discovered. Due to the limitations of the detection methods, our knowledge of the systems is usually far from complete. In particular, for planetary systems discovered with the radial velocity (RV) technique, the inclinations of the orbital planes, and thus the mutual inclinations and planetary masses, are unknown. Our work aims to constrain the spatial configuration of several RV-detected extrasolar systems that are not in a mean-motion resonance. Methods. Through an analytical study based on a first-order secular Hamiltonian expansion and numerical explorations performed with a chaos detector, we identified ranges of values for the orbital inclinations and the mutual inclinations, which ensure the long-term stability of the system. Our results were validated by comparison with n-body simulations, showing the accuracy of our analytical approach up to high mutual inclinations (approx. 70{\deg}-80{\deg}). Results. We find that, given the current estimations for the parameters of the selected systems, long-term regular evolution of the spatial configurations is observed, for all the systems, i) at low mutual inclinations (typically less than 35{\deg}) and ii) at higher mutual inclinations, preferentially if the system is in a Lidov-Kozai resonance. Indeed, a rapid destabilisation of highly mutually inclined orbits is commonly observed, due to the significant chaos that develops around the stability islands of the Lidov-Kozai resonance. The extent of the Lidov-Kozai resonant region is discussed for ten planetary systems (HD 11506, HD 12661, HD 134987, HD 142, HD 154857, HD 164922, HD 169830, HD 207832, HD 4732, and HD 74156).

The quest of exo-Earths has become a prominent field. In this work, we study the stability of non-coplanar planetary configurations consisting of an inclined inner terrestrial planet in mean-motion resonance with an outer giant planet. We examine the families of circular and elliptic symmetric periodic orbits with respect to the vertical stability, and identify the vertical critical orbits from which the spatial families emanate. We showcase that stable spatial periodic orbits can exist for both prograde and retrograde motion in 3/2, 2/1, 5/2, 3/1, 4/1 and 5/1 resonances for broad ranges of inclinations, when the giant evolves on a circular orbit. When the orbit of the giant is elliptic, only the 2/1 resonance has stable periodic orbits up to high inclinations, while the 3/1, 4/1 and 5/1 resonances possess segments of stability for low inclinations. Furthermore, we show that regular motion can also take place in the vicinity of both horizontally and vertically stable planar periodic orbits, even for very high inclinations. Finally, the results are discussed in the context of asteroid dynamics.

Many triple-star systems have an inner pair with an orbital period of a few days only. A common mechanism to explain the short-period pile-up present in the observations is the migration through Lidov-Kozai cycles combined with tidal friction. Here, we revisit this mechanism and aim to determine the initial orbital configurations leading to this process. We show that the mutual inclination of the triple-star system is not the only critical parameter, since the eccentricity as well as the argument of the pericenter of the inner orbit also play an important role in the establishment of the Lidov-Kozai migration. Our framework is the secular hierarchical three-body problem (octupole order approximation) with general relativity corrections, including the effects of tides, stellar oblateness and magnetic spin-down braking. Both the orbital and the spin evolutions are considered. Extensive numerical simulations with uniform and non-uniform distributions of the initial orbital parameters are carried out, and unbiased initial conditions leading to Lidov-Kozai migration are revealed. Finally, we highlight the importance of the initial "Kozai constant" $h=\sqrt{(1-e^2)}\cos{i}$ in the dynamical evolution of triple-star systems, by showing that phase portraits at given $h$-values unveil different evolution paths.

We consider a planetary system consisting of two primaries, namely a star and a giant planet, and a massless secondary, say a terrestrial planet or an asteroid, which moves under their gravitational attraction. We study the dynamics of this system in the framework of the circular and elliptic restricted TBP, when the motion of the giant planet describes circular and elliptic orbits, respectively. Originating from the circular family, families of symmetric periodic orbits in the 3/2, 5/2, 3/1, 4/1 and 5/1 mean-motion resonances are continued in the circular and the elliptic problems. New bifurcation points from the circular to the elliptic problem are found for each of the above resonances and thus, new families, continued from these points are herein presented. Stable segments of periodic orbits were found at high eccentricity values of the already known families considered as whole unstable previously. Moreover, new isolated (not continued from bifurcation points) families are computed in the elliptic restricted problem. The majority of the new families mainly consist of stable periodic orbits at high eccentricities. The families of the 5/1 resonance are investigated for the first time in the restricted three-body problems. We highlight the effect of stable periodic orbits on the formation of stable regions in their vicinity and unveil the boundaries of such domains in phase space by computing maps of dynamical stability. The long-term stable evolution of the terrestrial planets or asteroids is dependent on the existence of regular domains in their dynamical neighbourhood in phase space, which could host them for long time spans. This study, besides other celestial architectures that can be efficiently modelled by the circular and elliptic restricted problems, is particularly appropriate for the discovery of terrestrial companions among the single-giant planet systems discovered so far.

Hundreds of giant planets have been discovered so far and the quest of exo-Earths in giant planet systems has become intriguing. In this work, we aim to address the question of the possible long-term coexistence of a terrestrial companion on an orbit interior to a giant planet, and explore the extent of the stability regions for both non-resonant and resonant configurations. Our study focuses on the restricted three-body problem, where an inner terrestrial planet (massless body) moves under the gravitational attraction of a star and an outer massive planet on a circular or elliptic orbit. Using the Detrended Fast Lyapunov Indicator as a chaotic indicator, we constructed maps of dynamical stability by varying both the eccentricity of the outer giant planet and the semi-major axis of the inner terrestrial planet, and identify the boundaries of the stability domains. Guided by the computation of families of periodic orbits, the phase space is unravelled by meticulously chosen stable periodic orbits, which buttress the stability domains. We provide all possible stability domains for coplanar symmetric configurations and show that a terrestrial planet, either in mean-motion resonance or not, can coexist with a giant planet, when the latter moves on either a circular or an (even highly) eccentric orbit. New families of symmetric and asymmetric periodic orbits are presented for the 2/1 resonance. It is shown that an inner terrestrial planet can survive long time spans with a giant eccentric outer planet on resonant symmetric orbits, even when both orbits are highly eccentric. For 22 detected single-planet systems consisting of a giant planet with high eccentricity, we discuss the possible existence of a terrestrial planet. This study is particularly suitable for the research of companions among the detected systems with giant planets, and could assist with refining observational data.

The late-stage formation of giant planetary systems is rich in interesting dynamical mechanisms. Previous simulations of three giant planets initially on quasi-circular and quasi-coplanar orbits in the gas disc have shown that highly mutually inclined configurations can be formed, despite the strong eccentricity and inclination damping exerted by the disc. Much attention has been directed to inclination-type resonance, asking for large eccentricities to be acquired during the migration of the planets. Here we show that inclination excitation is also present at small to moderate eccentricities in two-planet systems that have previously experienced an ejection or a merging and are close to resonant commensurabilities at the end of the gas phase. We perform a dynamical analysis of these planetary systems, guided by the computation of planar families of periodic orbits and the bifurcation of families of spatial periodic orbits. We show that inclination excitation at small to moderate eccentricities can be produced by (temporary) capture in inclination-type resonance and the possible proximity of the non-coplanar systems to spatial periodic orbits contributes to maintaining their mutual inclination over long periods of time.

Evidence of mutually inclined planetary orbits has been reported for giant planets these last years. Here we aim to study the impact of eccentric and inclined massive giant planets on the terrestrial planet formation process, and investigate whether it can possibly lead to the existence of inclined terrestrial planets. We have performed 126 simulations of the late-stage planetary accretion in eccentric and inclined giant planet systems. The physical and orbital parameters of the giant planet systems result from n-body simulations of three giant planets in the late stage of the gas disc, under the combined action of Type II migration and planet-planet scattering. Fourteen two- and three-planet configurations have been selected, with diversified masses, semi-major axes (resonant configurations or not), eccentricities and inclinations (including coplanar systems) at the dispersal of the gas disc. We have then followed the gravitational interactions of these systems with an inner disc of planetesimals and embryos (9 runs per system), studying in detail the final configurations of the formed terrestrial planets. While coplanar giant systems harbour several terrestrial planets, generally as massive as the Earth and mainly on low eccentric and low inclined orbits, terrestrial planets formed in systems with mutually inclined giant planets are usually fewer, less massive (<0.5 M_{\Earth}) and with larger eccentricities and inclinations. This work shows that terrestrial planets can form on stable inclined orbits through the classical accretion theory, even in coplanar giant planet systems emerging from the disc phase.

We have found many Earth-sized worlds but we have no way of determining if their surfaces are Earth-like. This makes it impossible to quantitatively compare habitability, and pretending we can risks damaging the field.

We aim to investigate the influence of the eccentricity and inclination damping due to planet-disc interactions on the final configurations of the systems, generalizing previous studies on the combined action of the gas disc and planet-planet scattering during the disc phase. Instead of the simplistic $K$-prescription, our n-body simulations adopt the damping formulae for eccentricity and inclination provided by the hydrodynamical simulations of our companion paper. We follow the evolution of $11000$ numerical experiments of three giant planets in the late stage of the gas disc, exploring different initial configurations, planetary mass ratios and disc masses. The dynamical evolutions of the planetary systems are studied along the simulations, with emphasis on the resonance captures and inclination-growth mechanisms. Most of the systems are found with small inclinations ($\le10^{\circ}$) at the dispersal of the disc. Even though many systems enter an inclination-type resonance during the migration, the disc usually damps the inclinations on a short timescale. Although the majority of the multiple systems in our results are quasi-coplanar, $\sim5\%$ of them end up with high mutual inclinations ($\ge10^{\circ}$). Half of these highly mutually inclined systems result from two- or three-body MMR captures, the other half being produced by orbital instability and/or planet-planet scattering. When considering the long-term evolution over $100$ Myr, destabilization of the resonant systems is common, and the percentage of highly mutually inclined systems still evolving in resonance drops to $30\%$. Finally, the parameters of the final system configurations are in very good agreement with the semi-major axis and eccentricity distributions in the observations, showing that planet-planet interactions during the disc phase could have played an important role in sculpting planetary systems.

Contrary to Earth, which has a small orbital eccentricity, some exoplanets discovered in the insolation habitable zone (HZ) have high orbital eccentricities (e.g., up to an eccentricity of $\sim0.97$ for HD~20782~b). This raises the question of the capacity of these planets to host surface liquid water. In order to assess the habitability of an eccentric planet, the mean flux approximation is often used. It states that a planet on an eccentric orbit is called habitable if it receives on average a flux compatible with the presence of surface liquid water. However, as the planets do experience important insolation variations over one orbit and even spend some time outside the HZ for high eccentricities, the question of their habitability might not be as straightforward. We performed a set of simulations using the Global Climate Model LMDz, exploring the limits of the mean flux approximation when varying the luminosity of the host star and the eccentricity of the planet. We computed the climate of tidally locked ocean covered planets with orbital eccentricity from 0 to 0.9 receiving a mean flux equal to Earth's, around stars of luminosity ranging from $L_\odot$ to $10^{-4}~L_\odot$. Using here a definition of habitability based on the presence of surface liquid water, we find that most of the planets considered can sustain surface liquid water on the dayside with an ice cap on the nightside. However, for high eccentricity and high luminosity, planets cannot sustain surface liquid water during the whole orbital period. They completely freeze at apoastron and when approaching periastron an ocean appears around the substellar point. We conclude that the higher the eccentricity and the higher the luminosity of the star, the less reliable the mean flux approximation.

We study the secular evolution of several exoplanetary systems by extending the Laplace-Lagrange theory to order two in the masses. Using an expansion of the Hamiltonian in the Poincar\'e canonical variables, we determine the fundamental frequencies of the motion and compute analytically the long-term evolution of the keplerian elements. Our study clearly shows that, for systems close to a mean-motion resonance, the second order approximation describes their secular evolution more accurately than the usually adopted first order one. Moreover, this approach takes into account the influence of the mean anomalies on the secular dynamics. Finally, we set up a simple criterion that is useful to discriminate between three different categories of planetary systems: (i) {\it secular} systems (HD 11964, HD 74156, HD 134987, HD 163607, HD 12661 and HD 147018); (ii) systems {\it near a mean-motion resonance} (HD 11506, HD 177830, HD 9446, HD 169830 and $\upsilon$ Andromedae); (iii) systems {\it really close to} or {\it in a mean-motion resonance} (HD 108874, HD 128311 and HD 183263).

In the Solar System, planets have a small inclination with respect to the equatorial plane of the Sun, but there is evidence that in extrasolar systems the inclination can be very high. This spin-orbit misalignment is unexpected, as planets form in a protoplanetary disc supposedly aligned with the stellar spin. Planet-planet interactions are supposed to lead to a mutual inclination, but the effects of the protoplanetary disc are still unknown. We investigate therefore planet-disc interactions for planets above 1M_Jup. We check the influence of the inclination i, eccentricity e, and mass M_p of the planet. We perform 3D numerical simulations of protoplanetary discs with embedded high-mass planets. We provide damping formulae for i and e as a function of i, e, and M_p that fit the numerical data. For highly inclined massive planets, the gap opening is reduced, and the damping of i occurs on time-scales of the order of 10^-4 deg/yr M_disc/(0.01 M_star) with the damping of e on a smaller time-scale. While the inclination of low planetary masses (<5M_Jup) is always damped, large planetary masses with large i can undergo a Kozai-cycle with the disc. These Kozai-cycles are damped in time. Eccentricity is generally damped, except for very massive planets (M_p = 5M_Jup) where eccentricity can increase for low inclinations. The dynamics tends to a final state: planets end up in midplane and can then, over time, increase their eccentricity as a result of interactions with the disc. The interactions with the disc lead to damping of i and e after a scattering event of high-mass planets. If i is sufficiently reduced, the eccentricity can be pumped up because of interactions with the disc. If the planet is scattered to high inclination, it can undergo a Kozai-cycle with the disc that makes it hard to predict the exact movement of the planet and its orbital parameters at the dispersal of the disc.

Transit timing variations (TTVs) are useful to constrain the existence of perturbing planets, especially in resonant systems where the variations are strongly enhanced. Here we focus on Laplace-resonant three-planet systems, and assume the inner planet transits the star. A dynamical study is performed for different masses of the three bodies, with a special attention to terrestrial planets. We consider a maximal time-span of ~ 100 years and discuss the shape of the inner planet TTVs curve. Using frequency analysis, we highlight the three periods related to the evolution of the system: two periods associated with the Laplace-resonant angle and the third one with the precession of the pericenters. These three periods are clearly detected in the TTVs of an inner giant planet perturbed by two terrestrial companions. Only two periods are detected for a Jupiter-Jupiter-Earth configuration (the ones associated with the giant interactions) or for three terrestrial planets (the Laplace periods). However, the latter system can be constrained from the inner planet TTVs. We finally remark that the TTVs of resonant three or two Jupiter systems mix up, when the period of the Laplace resonant angle matches the pericenter precession of the two-body configuration. This study highlights the importance of TTVs long-term observational programs for the detection of multiple-planet resonant systems.

We study the establishment of three-planet resonances -similar to the Laplace resonance in the Galilean satellites- and their effects on the mutual inclinations of the orbital planes of the planets, assuming that the latter undergo migration in a gaseous disc. In particular, we examine the resonance relations that occur, by varying the physical and initial orbital parameters of the planets (mass, initial semi-major axis and eccentricity) as well as the parameters of the migration forces (migration rate and eccentricity damping rate), which are modeled here through a simplified analytic prescription. We find that, in general, for planetary masses below 1.5 M_J, multiple-planet resonances of the form n3:n2:n1=1:2:4 and 1:3:6 are established, as the inner planets, m1 and m2, get trapped in a 1:2 resonance and the outer planet m3 subsequently is captured in a 1:2 or 1:3 resonance with m2. For mild eccentricity damping, the resonance pumps the eccentricities of all planets on a relatively short time-scale, to the point where they enter an inclination-type resonance (as in Libert & Tsiganis 2011); then mutual inclinations can grow to ~35{\deg}, thus forming a "3-D system". On the other hand, we find that trapping of m2 in a 2:3 resonance with m1 occurs very rarely, for the range of masses used here, so only two cases of capture in a respective three-planet resonance were found. Our results suggest that trapping in a three-planet resonance can be common in exoplanetary systems, provided that the planets are not very massive. Inclination pumping could then occur relatively fast, provided that eccentricity damping is not very efficient so that at least one of the inner planets acquires an orbital eccentricity higher than e=0.3.

Assuming that giant planets are formed in thin protoplanetary discs, a '3D' system can form, provided that the mutual inclination is excited by some dynamical mechanism. Resonant interactions and close planetary encounters are thought to be the primary inclination-excitation mechanisms, resulting in a resonant and non-resonant system, respectively. Here we propose an alternative formation scenario, starting from a system composed of three giant planets in a nearly coplanar configuration. As was recently shown for the case of the Solar system, planetary migration in the gas disc (Type II migration) can force the planets to become trapped in a multiply resonant state. We simulate this process, assuming different values for the planetary masses and mass ratios. We show that such a triple resonance generally becomes unstable as the resonance excites the eccentricities of all planets and planet-planet scattering sets in. One of the three planets is typically ejected from the system, leaving behind a dynamically 'hot' (but stable) two-planet configuration. The resulting two-planet systems typically have large values of semimajor axial ratios (a1/a2 < 0.3), while the mutual inclination can be as high as 70{\deg}, with a median of \sim30{\deg}. A small fraction of our two-planet systems (\sim5 per cent) ends up in the stability zone of the Kozai resonance. In a few cases, the triple resonance can remain stable for long times and a '3D' system can form by resonant excitation of the orbital inclinations; such a three-planet system could be stable if enough eccentricity damping is exerted on the planets. Finally, in the single-planet resulting systems, which are formed when two planets are ejected from the system, the inclination of the planet's orbital plane with respect to the initial invariant plane -presumably the plane perpendicular to the star's spin axis- can be as large as \sim40{\deg}.

We study the dynamics of the 3-D three-body problem of a small body moving under the attractions of a star and a giant planet which orbits the star on a much wider and elliptic orbit. In particular, we focus on the influence of an eccentric orbit of the outer perturber on the dynamics of a small highly inclined inner body. Our analytical study of the secular perturbations relies on the classical octupole hamiltonian expansion (third-order theory in the ratio of the semi-major axes), as third-order terms are needed to consider the secular variations of the outer perturber and potential secular resonances between the arguments of the pericenter and/or longitudes of the node of both bodies. Short-period averaging and node reduction (Laplace plane) reduce the problem to two degrees of freedom. The four-dimensional dynamics is analyzed through representative planes which identify the main equilibria of the problem. As in the circular problem (i.e. perturber on a circular orbit), the "Kozai-bifurcated" equilibria play a major role in the dynamics of an inner body on quasi-circular orbit: its eccentricity variations are very limited for mutual inclination between the orbital planes smaller than ~40^{\deg}, while they become large and chaotic for higher mutual inclination. Particular attention is also given to a region around 35^{\deg} of mutual inclination, detected numerically by Funk et al. (2011) and consisting of long-time stable and particularly low eccentric orbits of the small body. Using a 12th-order Hamiltonian expansion in eccentricities and inclinations, in particular its action-angle formulation obtained by Lie transforms in Libert & Henrard (2008), we show that this region presents an equality of two fundamental frequencies and can be regarded as a secular resonance. Our results also apply to binary star systems where a planet is revolving around one of the two stars.

In this work we propose a new numerical approach to distinguish between regular and chaotic orbits in Hamiltonian systems, based on the simultaneous integration of both the orbit and the deviation vectors using a symplectic scheme, hereby called {\em global symplectic integrator}. In particular, the proposed method allows us to recover the correct orbits character with very large integration time steps, small energy losses and short CPU times. To illustrate the numerical performances of the global symplectic integrator we will apply it to two well-known and widely studied problems: the H\'enon-Heiles model and the restricted three-body problem.