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.

The birth and death of planets may be affected by mass outflows from their parent stars during the T-Tauri or post-main-sequence phases of stellar evolution. These outflows are often modelled to be isotropic, but this assumption is not realistic for fast rotators, bipolar jets and supernovae. Here we derive the general equations of motion for the time evolution of a single planet, brown dwarf, comet or asteroid perturbed by anisotropic mass loss in terms of a complete set of planetary orbital elements, the ejecta velocity, and the parent star's co-latitude and longitude. We restrict our application of these equations to 1) rapidly rotating giant stars, and 2) arbitrarily-directed jet outflows. We conclude that the isotropic mass-loss assumption can safely be used to model planetary motion during giant branch phases of stellar evolution within distances of hundreds of au. In fact, latitudinal mass loss variations anisotropically affect planetary motion only if the mass loss is asymmetric about the stellar equator. Also, we demonstrate how constant-velocity, asymmetric bipolar outflows in young systems incite orbital inclination changes. Consequently, this phenomenon readily tilts exoplanetary orbits external to a nascent disc on the order of degrees.

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.

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.

We present families of symmetric and asymmetric periodic orbits at the 1/1 resonance, for a planetary system consisting of a star and two small bodies, in comparison to the star, moving in the same plane under their mutual gravitational attraction. The stable 1/1 resonant periodic orbits belong to a family which has a planetary branch, with the two planets moving in nearly Keplerian orbits with non zero eccentricities and a satellite branch, where the gravitational interaction between the two planets dominates the attraction from the star and the two planets form a close binary which revolves around the star. The stability regions around periodic orbits along the family are studied. Next, we study the dynamical evolution in time of a planetary system with two planets which is initially trapped in a stable 1/1 resonant periodic motion, when a drag force is included in the system. We prove that if we start with a 1/1 resonant planetary system with large eccentricities, the system migrates, due to the drag force, {\it along the family of periodic orbits} and is finally trapped in a satellite orbit. This, in principle, provides a mechanism for the generation of a satellite system: we start with a planetary system and the final stage is a system where the two small bodies form a close binary whose center of mass revolves around the star.

The continuation of resonant periodic orbits from the restricted to the general three body problem is studied in a systematic way. Starting from the Keplerian unperturbed system we obtain the resonant families of the circular restricted problem. Then we find all the families of the resonant elliptic restricted three body problem which bifurcate from the circular model. All these families are continued to the general three body problem, and in this way we can obtain a global picture of all the families of periodic orbits of a two-planet resonant system. We consider planar motion only. We show that the continuation follows a scheme proposed by Bozis and Hadjidemetriou (1976) for symmetric orbits. Our study includes also asymmetric periodic orbits, which exist in cases of external resonances. The families formed by passing from the restricted to the general problem are continued within the framework of the general problem by varying the planetary mass ratio $\rho$. We obtain bifurcations which are caused either due to collisions of the families in the space of initial conditions or due to the vanishing of bifurcation points. Our study refers to the whole range of planetary mass ratio values ($\rho \in (0,\infty)$) and, therefore we include the passage from external to internal resonances. In the present work, our numerical study includes the case of the 2/1 and 1/2 resonance. The same method can be used to study all other planetary resonances.

We study the dynamics of planetary systems with two planets moving in the same plane, when frictional forces act on the two planets, in addition to the gravitational forces. The model of the general three-body problem is used. Different laws of friction are considered. The topology of the phase space is essential in understanding the evolution of the system. The topology is determined by the families of stable and unstable periodic orbits, both symmetric and non symmetric. It is along the stable families, or close to them, that the planets migrate when dissipative forces act. At the critical points where the stability along the family changes, there is a bifurcation of a new family of stable periodic orbits and the migration process changes route and follows the new stable family up to large eccentricities or to a chaotic region. We consider both resonant and non resonant planetary systems. The 2/1, 3/1 and 3/2 resonances are studied. The migration to larger or smaller eccentricities depends on the particular law of friction. Also, in some cases the semimajor axes increase and in other cases they are stabilized. For particular laws of friction and for special values of the parameters of the frictional forces, it is possible to have partially stationary solutions, where the eccentricities and the semimajor axes are fixed.