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.

Discoveries of exoplanets orbiting evolved stars motivate critical examinations of the dynamics of $N$-body systems with mass loss. Multi-planet evolved systems are particularly complex because of the mutual interactions between the planets. Here, we study the underlying dynamical mechanisms which can incite planetary escape in two-planet post-main sequence systems. Stellar mass loss alone is unlikely to be rapid and high enough to eject planets at typically-observed separations. However, the combination of mass loss and planet-planet interactions can prompt a shift from stable to chaotic regions of phase space. Consequently, when mass loss ceases, the unstable configuration may cause escape. By assuming a constant stellar mass loss rate, we utilize maps of dynamical stability to illustrate the distribution of regular and chaotic trajectories in phase space. We show that chaos can drive the planets to undergo close encounters, leading to the ejection of one planet. Stellar mass loss can trigger the transition of a planetary system from a stable to chaotic configuration, subsequently causing escape. We find that mass loss non-adiabatically affects planet-planet interaction for the most massive progenitor stars which avoid the supernova stage. For these cases, we present specific examples of planetary escape.

The recent discovery of free-floating planets and their theoretical interpretation as celestial bodies, either condensed independently or ejected from parent stars in tight clusters, introduced an intriguing possibility. Namely, that some exoplanets are not condensed from the protoplanetary disk of their parent star. In this novel scenario a free-floating planet interacts with an already existing planetary system, created in a tight cluster, and is captured as a new planet. In the present work we study this interaction process by integrating trajectories of planet-sized bodies, which encounter a binary system consisting of a Jupiter-sized planet revolving around a Sun-like star. To simplify the problem we assume coplanar orbits for the bound and the free-floating planet and an initially parabolic orbit for the free-floating planet. By calculating the uncertainty exponent, a quantity that measures the dependence of the final state of the system on small changes of the initial conditions, we show that the interaction process is a fractal classical scattering. The uncertainty exponent is in the range (0.2-0.3) and is a decreasing function of time. In this way we see that the statistical approach we follow to tackle the problem is justified. The possible final outcomes of this interaction are only four, namely flyby, planet exchange, capture or disruption. We give the probability of each outcome as a function of the incoming planet's mass. We find that the probability of exchange or capture (in prograde as well as retrograde orbits and for very long times) is non-negligible, a fact that might explain the possible future observations of planetary systems with orbits that are either retrograde or tight and highly eccentric.

The motion of a satellite around a planet can be studied by the Hill model, which is a modification of the restricted three body problem pertaining to motion of a satellite around a planet. Although the dynamics of the circular Hill model have been extensively studied in the literature, only few results about the dynamics of the elliptic model were known up to now, namely the equations of motion and few unstable families of periodic orbits. In the present study we extend these results by computing a large set of families of periodic orbits and their linear stability and classify them according to their resonance condition. Although most of them are unstable, we were able to find a considerable number of stable ones. By computing appropriate maps of dynamical stability, we study the effect of the planetary eccentricity on the stability of satellite orbits. We see that, even for large values of the planetary eccentricity, regular orbits can be found in the vicinity of stable periodic orbits. The majority of irregular orbits are escape orbits.

A model, plane symmetric, 3-D potential, which preserves some features of galactic problems,is used in order to examine the phase space structure through the study of the properties of orbits crossing perpendicularly the plane of symmetry. It is found that the lines formed by periodic orbits, belonging to Farey sequences, are not smooth neither continuous. Instead they are deformed and broken in regions characterised by high Lyapunov Characteristic Numbers (LCN's). It is suggested that these lines are an incomplete form of a fractal repeller, as discussed by Gaspard and Baras (1995), and are thus closely associated to the ``quasi-barriers'' discussed by Varvoglis et al. (1997). There are numerical indications that the contour lines of constant LCN's possess fractal properties. Finally it is shown numerically that some of the periodic orbits -members of the lines- belong to true irregular families. It is argued that the fractal properties of the phase space should affect the transport of trajectories in phase or action space and,therefore, play a certain role in the chaotic motion of stars in more realistic galactic potentials.