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.

Planetary debris is observed in the atmospheres of over 1,000 white dwarfs, and two white dwarfs are now observed to contain orbiting minor planets. Exoasteroids and planetary core fragments achieve orbits close to the white dwarf through scattering with major planets. However, the architectures that allow for this scattering to take place are time-consuming to explore with $N$-body simulations lasting $\sim 10^{10}$ yr; these long-running simulations restrict the amount of phase space that can be investigated. Here we use planar and three-dimensional (spatial) elliptic periodic orbits, as well as chaotic indicators through dynamical stability maps, as quick scale-free analytic alternatives to $N$-body simulations in order to locate and predict instability in white dwarf planetary systems that consist of one major and one minor planet on very long timescales. We then classify the instability according to ejection versus collisional events. We generalized our previous work by allowing eccentricity and inclination of the periodic orbits to increase, thereby adding more realism but also significantly more degrees of freedom to our architectures. We also carried out a suite of computationally expensive 10 Gyr $N$-body simulations to provide comparisons with chaotic indicators in a limited region of phase space. We compute dynamical stability maps that are specific to white dwarf planetary systems and that can be used as tools in future studies to quickly estimate pollution prospects and timescales for one-planet architectures. We find that these maps also agree well with the outcomes of our $N$-body simulations. As observations of metal-polluted white dwarfs mount exponentially, particularly in the era of Gaia, tools such as periodic orbits can help infer dynamical histories for ensembles of systems.

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.

Applying the method of analytical continuation of periodic orbits, we study quasi-satellite motion in the framework of the three-body problem. In the simplest, yet not trivial model, namely the planar circular restricted problem, it is known that quasi-satellite motion is associated with a family of periodic solutions, called family $f$, which consists of 1:1 resonant retrograde orbits. In our study, we determine the critical orbits of family $f$ that are continued both in the elliptic and in the spatial model and compute the corresponding families that are generated and consist the backbone of the quasi-satellite regime in the restricted model. Then, we show the continuation of these families in the general three-body problem, we verify and explain previous computations and show the existence of a new family of spatial orbits. The linear stability of periodic orbits is also studied. Stable periodic orbits unravel regimes of regular motion in phase space where 1:1 resonant angles librate. Such regimes, which exist even for high eccentricities and inclinations, may consist dynamical regions where long-lived asteroids or co-orbital exoplanets can be found.

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.

Librational motion in celestial mechanics is generally associated with the existence of stable resonant configurations and signified by the existence of stable periodic solutions and oscillation of critical (resonant) angles. When such an oscillation takes place around a value different than 0 or $\pi$, the libration is called asymmetric. In the context of the planar circular restricted three-body problem (CRTBP), asymmetric librations have been identified for the exterior mean-motion resonances (MMRs) 1:2, 1:3 etc. as well as for co-orbital motion (1:1). In exterior MMRs the massless body is the outer one. In this paper, we study asymmetric librations in the 3-dimensional space. We employ the computational approach of Markellos (1978) and compute families of asymmetric periodic orbits and their stability. Stable, asymmetric periodic orbits are surrounded in phase space by domains of initial conditions which correspond to stable evolution and librating resonant angles. Our computations were focused on the spatial circular restricted three-body model of the Sun-Neptune-TNO system (TNO= trans-Neptunian object). We compare our results with numerical integrations of observed TNOs, which reveal that some of them perform 1:2-resonant, inclined asymmetric librations. For the stable 1:2 TNOs librators, we find that their libration seems to be related with the vertically stable planar asymmetric orbits of our model, rather than the 3-dimensional ones found in the present study.

We consider planetary systems evolving under the effect of a Stokes-type dissipative force mimicking the outcome of a type II migration process. As inward migration proceeds and the planets follow the circular family (they start on circular orbits) and even though they are initially almost coplanar, resonance capture can be realized. Then, at the \textit{vertical critical orbits} (VCOs), that the circular family possesses, the inclination excitation can abruptly take place. The planets are now guided by the spatial elliptic families, which bifurcate from those critical orbits. We herein, perform a direct link of mutually inclined stable planetary systems on circular orbits trapped in \textit{mean-motion resonance} (MMR) with the existence of VCOs of high values of multiplicity. It is shown that the more the multiplicity of the periodic orbits of the circular family increases, the more VCOs (corresponding to more MMRs) appear. In this way, we can provide a justification for the existence of resonant planets on circular orbits, which could, even further to that, evolve stably if they were mutually inclined.

Nowadays, many extrasolar planetary systems possessing at least one planet on a highly eccentric orbit have been discovered. In this work, we study the possible long-term stability of such systems. We consider the general three body problem as our model. Highly eccentric orbits are out of the Hill stability regions. However, mean motion resonances can provide phase protection and orbits with long-term stability exist. We construct maps of dynamical stability based on the computation of chaotic indicators and we figure out regions in phase space, where the long-term stability is guaranteed. We focus on regions where at least one planet is highly eccentric and attempt to associate them with the existence of stable periodic orbits. The values of the orbital elements, which are derived from observational data, are often given with very large deviations. Generally, phase space regions of high eccentricities are narrow and thus, our dynamical analysis may restrict considerably the valid domain of the system's location.

Mounting discoveries of debris discs orbiting newly-formed stars and white dwarfs (WDs) showcase the importance of modeling the long-term evolution of small bodies in exosystems. WD debris discs are in particular thought to form from very long-term (0.1-5.0 Gyr) instability between planets and asteroids. However, the time-consuming nature of $N$-body integrators which accurately simulate motion over Gyrs necessitates a judicious choice of initial conditions. The analytical tools known as \textit{periodic orbits} can circumvent the guesswork. Here, we begin a comprehensive analysis directly linking periodic orbits with $N$-body integration outcomes with an extensive exploration of the planar circular restricted three-body problem (CRTBP) with an outer planet and inner asteroid near or inside of the $2$:$1$ mean motion resonance. We run nearly 1000 focused simulations for the entire age of the Universe (14 Gyr) with initial conditions mapped to the phase space locations surrounding the unstable and stable periodic orbits for that commensurability. In none of our simulations did the planar CRTBP architecture yield a long-timescale ($\gtrsim 0.25$% of the age of the Universe) asteroid-star collision. The pericentre distance of asteroids which survived beyond this timescale ($\approx 35$ Myr) varied by at most about 60%. These results help affirm that collisions occur too quickly to explain WD pollution in the planar CRTBP $2$:$1$ regime, and highlight the need for further periodic orbit studies with the eccentric and inclined TBP architectures and other significant orbital period commensurabilities.

Many of exoplanetary systems consist of more than one planet and the study of planetary orbits with respect to their long-term stability is very interesting. Furthermore, many exoplanets seem to be locked in a mean-motion resonance (MMR), which offers a phase protection mechanism, so that, even highly eccentric planets can avoid close encounters. However, the present estimation of their initial conditions, which may change significantly after obtaining additional observational data in the future, locate most of the systems in chaotic regions and consequently, they are destabilized. Hence, dynamical analysis is imperative for the derivation of proper planetary orbital elements. We utilize the model of spatial general three body problem, in order to simulate such resonant systems through the computation of families periodic orbits. In this way, we can figure out regions in phase space, where the planets in resonances should be ideally hosted in favour of long-term stability and therefore, survival. In this review, we summarize our methodology and showcase the fact that stable resonant planetary systems evolve being exactly centered at stable periodic orbits. We apply this process to co-orbital motion and systems HD 82943, HD 73526, HD 128311, HD 60532, HD 45364 and HD 108874.

The long-term stability of the evolution of two-planet systems is considered by using the general three body problem (GTBP). Our study is focused on the stability of systems with adjacent orbits when at least one of them is highly eccentric. In these cases, in order for close encounters, which destabilize the planetary systems, to be avoided, phase protection mechanisms should be considered. Additionally, since the GTBP is a non-integrable system, chaos may also cause the destabilization of the system after a long time interval. By computing dynamical maps, based on Fast Lyapunov Indicator, we reveal regions in phase space with stable orbits even for very high eccentricities (e>0.5). Such regions are present in mean motion resonances (MMR). We can determine the position of the exact MMR through the computation of families of periodic orbits in a rotating frame. Elliptic periodic orbits are associated with the presence of apsidal corotation resonances (ACR). When such solutions are stable, they are associated with neighbouring domains of initial conditions that provide long-term stability. We apply our methodology so that the evolution of planetary systems of highly eccentric orbits is assigned to the existence of such stable domains. Particularly, we study the orbital evolution of the extrasolar systems HD 82943, HD 3651, HD 7449, HD 89744 and HD 102272 and discuss the consistency between the orbital elements provided by the observations and the dynamical stability.

We herein utilize the general three-body problem (GTBP) as a model, in order to simulate resonant systems consisting of a star and two planets, where at least one of them is highly eccentric. We study them in terms of their long-term stability, via the construction of maps of dynamical stability and the computation of the corresponding families of periodic orbits. We identify the way their survival is connected with the regions of regular motion in phase space, which, in turn, were created by stable resonant periodic orbits in their vicinity. Consequently, a phase protection mechanism is provided and the planets avoid close encounters and collisions even on long timescales. We apply our methodology to the extrasolar system HD 82943.

We study the dynamics of a two-planet system, which evolves being in a $1/1$ mean motion resonance (co-orbital motion) with non-zero mutual inclination. In particular, we examine the existence of bifurcations of periodic orbits from the planar to the spatial case. We find that such bifurcations exist only for planetary mass ratios $\rho=\frac{m_2}{m_1}<0.0205$. For $\rho$ in the interval $0<\rho<0.0205$, we compute the generated families of spatial periodic orbits and their linear stability. These spatial families form bridges, which start and end at the same planar family. Along them the mutual planetary inclination varies. We construct maps of dynamical stability and show the existence of regions of regular orbits in phase space.

In dynamical systems of few degrees of freedom, periodic solutions consist the backbone of the phase space and the determination and computation of their stability is crucial for understanding the global dynamics. In this paper we study the classical three body problem in three dimensions and use its dynamics to assess the long-term evolution of extrasolar systems. We compute periodic orbits, which correspond to exact resonant motion, and determine their linear stability. By computing maps of dynamical stability we show that stable periodic orbits are surrounded in phase space with regular motion even in systems with more than two degrees of freedom, while chaos is apparent close to unstable ones. Therefore, families of stable periodic orbits, indeed, consist backbones of the stability domains in phase space.

We consider a two-planet system, which migrates under the influence of dissipative forces that mimic the effects of gas-driven (Type II) migration. It has been shown that, in the planar case, migration leads to resonant capture after an evolution that forces the system to follow families of periodic orbits. Starting with planets that differ slightly from a coplanar configuration, capture can, also, occur and, additionally, excitation of planetary inclinations has been observed in some cases. We show that excitation of inclinations occurs, when the planar families of periodic orbits, which are followed during the initial stages of planetary migration, become vertically unstable. At these points, {\em vertical critical orbits} may give rise to generating stable families of $3D$ periodic orbits, which drive the evolution of the migrating planets to non-coplanar motion. We have computed and present here the vertical critical orbits of the $2/1$ and $3/1$ resonances, for various values of the planetary mass ratio. Moreover, we determine the limiting values of eccentricity for which the "inclination resonance" occurs.

The planetary dynamics of $4/3$, $3/2$, $5/2$, $3/1$ and $4/1$ mean motion resonances is studied by using the model of the general three body problem in a rotating frame and by determining families of periodic orbits for each resonance. Both planar and spatial cases are examined. In the spatial problem, families of periodic orbits are obtained after analytical continuation of vertical critical orbits. The linear stability of orbits is also examined. Concerning initial conditions nearby stable periodic orbits, we obtain long-term planetary stability, while unstable orbits are associated with chaotic evolution that destabilizes the planetary system. Stable periodic orbits are of particular importance in planetary dynamics, since they can host real planetary systems. We found stable orbits up to $60^\circ$ of mutual planetary inclination, but in most families, the stability does not exceed $20^\circ$-$30^\circ$, depending on the planetary mass ratio. Most of these orbits are very eccentric. Stable inclined circular orbits or orbits of low eccentricity were found in the $4/3$ and $5/2$ resonance, respectively.

Over the last decades, there has been a tremendous increase in research on extrasolar planets. Many exosolar systems, which consist of a Star and two inclined Planets, seem to be locked in 4/3, 3/2, 2/1, 5/2, 3/1 and 4/1 mean motion resonance (MMR). We herewith present the model used to simulate three dimensional planetary systems and provide planar families of periodic orbits (PO), which belong to all possible configurations that each MMR has, along with their linear horizontal and vertical stability. We focus on depicting stable spatial families (most of them up to mutual inclination of $60^\circ$) generated by PO of planar circular families, because the trapping in MMR could be a consequence of planetary migration process. We attempt to connect the linear stability of PO with long-term stability of a planetary system close to them. This can stimulate the search of real planetary systems in the vicinity of stable spatial PO-counterbalanced by the planets' orbital elements, masses and MMR; all of which could constitute a suitable environment convenient to host them.

Many exo-solar systems discovered in the last decade consist of planets orbiting in resonant configurations and consequently, their evolution should show long-term stability. However, due to the mutual planetary interactions a multi-planet system shows complicated dynamics with mostly chaotic trajectories. We can determine possible stable configurations by computing resonant periodic trajectories of the general planar three body problem, which can be used for modeling a two-planet system. In this work, we review our model for both the planar and the spatial case. We present families of symmetric periodic trajectories in various resonances and study their linear horizontal and vertical stability. We show that around stable periodic orbits there exist regimes in phase space where regular evolution takes place. Unstable periodic orbits are associated with the existence of chaos and planetary destabilization.

We consider the general spatial three body problem and study the dynamics of planetary systems consisting of a star and two planets which evolve into 2/1 mean motion resonance and into inclined orbits. Our study is focused on the periodic orbits of the system given in a suitable rotating frame. The stability of periodic orbits characterize the evolution of any planetary system with initial conditions in their vicinity. Stable periodic orbits are associated with long term regular evolution, while unstable periodic orbits are surrounded by regions of chaotic motion. We compute many families of symmetric periodic orbits by applying two schemes of analytical continuation. In the first scheme, we start from the 2/1 (or 1/2) resonant periodic orbits of the restricted problem and in the second scheme, we start from vertical critical periodic orbits of the general planar problem. Most of the periodic orbits are unstable, but many stable periodic orbits have been, also, found with mutual inclination up to $50^\circ$ - $60^\circ$, which may be related with the existence of real planetary systems.

We consider the planar three body problem of planetary type and we study the generation and continuation of periodic orbits and mainly of asymmetric periodic orbits. Asymmetric orbits exist in the restricted circular three body problem only in particular resonances called "asymmetric resonances". However, numerical studies showed that in the general three body problem asymmetric orbits may exist not only for asymmetric resonances, but for other kinds, too. In this work, we show the existence of asymmetric periodic orbits in the elliptic restricted problem. These orbits are continued and clarify the origin of many asymmetric periodic orbits in the general problem. Also, we illustrate how the families of periodic orbits of the restricted circular problem and those of the elliptic one join smoothly and form families in the general problem, verifying in this way the scenario described firstly by Bozis and Hadjidemetriou (1976).