In this paper, we study the long-term dynamical evolution of highly-elliptical orbits (HEOs) in the medium-Earth orbit (MEO) region around the Earth. The real population consists primarily of Geosynchronous Transfer Orbits (GTOs), launched at specific inclinations, Molniya-type satellites and related debris. We performed a suite of long-term numerical integrations (up to 200 years) within a realistic dynamical model, aimed primarily at recording the dynamical lifetime of such orbits (defined as the time needed for atmospheric reentry) and understanding its dependence on initial conditions and other parameters, such as the area-to-mass ratio (A/m). Our results are presented in the form of 2-D lifetime maps, for different values of inclination, A/m, and drag coefficient. We find that the majority of small debris (> 70%, depending on the inclination) can naturally reenter within 25-90 years, but these numbers are significantly less optimistic for large debris (e.g., upper stages), with the notable exception of those launched from high latitude (Baikonur). We estimate the reentry probability and mean dynamical lifetime for different classes of GTOs and we find that both quantities depend primarily and strongly on initial perigee altitude. Atmospheric drag and higher A/m values extend the reentry zones, especially at low inclinations. For high inclinations, this dependence is weakened, as the primary mechanisms leading to reentry are overlapping lunisolar resonances. This study forms part of the EC-funded (H2020) "ReDSHIFT" project.

We have carried out a numerical investigation of the coupled gravitational and non-gravitational perturbations acting on Earth satellite orbits in an extensive grid, covering the whole circumterrestrial space, using an appropriately modified version of the SWIFT symplectic integrator, which is suitable for long-term (120 years) integrations of the non-averaged equations of motion. Hence, we characterize the long-term dynamics and the phase-space structure of the Earth-orbiter environment, starting from low altitudes (400 km) and going up to the GEO region and beyond. This investigation was done in the framework of the EC-funded "ReDSHIFT" project, with the purpose of enabling the definition of passive debris removal strategies, based on the use of physical mechanisms inherent in the complex dynamics of the problem (i.e., resonances). Accordingly, the complicated interactions among resonances, generated by different perturbing forces (i.e., lunisolar gravity, solar radiation pressure, tesseral harmonics in the geopotential) are accurately depicted in our results, where we can identify the regions of phase space where the motion is regular and long-term stable and regions for which eccentricity growth and even instability due to chaotic behavior can emerge. The results are presented in an "atlas" of dynamical stability maps for different orbital zones, with a particular focus on the (drag-free) range of semimajor axes, where the perturbing effects of the Earth's oblateness and lunisolar gravity are of comparable order. In some regions, the overlapping of the predominant lunisolar secular and semi-secular resonances furnish a number of interesting disposal hatches at moderate to low eccentricity orbits. All computations were repeated for an increased area-to-mass ratio, simulating the case of a satellite equipped with an on-board, area-augmenting device.

The Medium Earth Orbit (MEO) region hosts satellites for navigation, communication, and geodetic/space environmental science, among which are the Global Navigation Satellites Systems (GNSS). Safe and efficient removal of debris from MEO is problematic due to the high cost for maneuvers needed to directly reach the Earth (reentry orbits) and the relatively crowded GNSS neighborhood (graveyard orbits). Recent studies have highlighted the complicated secular dynamics in the MEO region, but also the possibility of exploiting these dynamics, for designing removal strategies. In this paper, we present our numerical exploration of the long-term dynamics in MEO, performed with the purpose of unveiling the set of reentry and graveyard solutions that could be reached with maneuvers of reasonable DV cost. We simulated the dynamics over 120-200 years for an extended grid of millions of fictitious MEO satellites that covered all inclinations from 0 to 90deg, using non-averaged equations of motion and a suitable dynamical model that accounted for the principal geopotential terms, 3rd-body perturbations and solar radiation pressure (SRP). We found a sizeable set of usable solutions with reentry times that exceed ~40years, mainly around three specific inclination values: 46deg, 56deg, and 68deg; a result compatible with our understanding of MEO secular dynamics. For DV <= 300 m/s (i.e., achieved if you start from a typical GNSS orbit and target a disposal orbit with e<0.3), reentry times from GNSS altitudes exceed ~70 years, while low-cost (DV ~= 5-35 m/s) graveyard orbits, stable for at lest 200 years, are found for eccentricities up to e~0.018. This investigation was carried out in the framework of the EC-funded "ReDSHIFT" project.

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.

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.

The rectilinear elliptic restricted Three Body Problem (TBP) is the limiting case of the elliptic restricted TBP when the motion of the primaries is described by a Keplerian ellipse with eccentricity $e'=1$, but the collision of the primaries is assumed to be a non-singular point. The rectilinear model has been proposed as a starting model for studying the dynamics of motion around highly eccentric binary systems. Broucke (1969) explored the rectilinear problem and obtained isolated periodic orbits for mass parameter $\mu=0.5$ (equal masses of the primaries). We found that all orbits obtained by Broucke are linearly unstable. We extend Broucke's computations by using a finer search for symmetric periodic orbits and computing their linear stability. We found a large number of periodic orbits, but only eight of them were found to be linearly stable and are associated with particular mean motion resonances. These stable orbits are used as generating orbits for continuation with respect to $\mu$ and $e'<1$. Also, continuation of periodic solutions with respect to the mass of the small body can be applied by using the general TBP. FLI maps of dynamical stability show that stable periodic orbits are surrounded in phase space with regions of regular orbits indicating that systems of very highly eccentric orbits can be found in stable resonant configurations. As an application we present a stability study for the planetary system HD7449.

Periodic solutions of the three body problem are very important for understanding its dynamics either in a theoretical framework or in various applications in celestial mechanics. In this paper we discuss the computation and continuation of periodic orbits for planetary systems. The study is restricted to coplanar motion. Staring from known results of two-planet systems around single stars, we perform continuation of solutions with respect to the mass and approach periodic orbits of single planets in two-star systems. Also, families of periodic solutions can be computed for fixed masses of the primaries. When they are linearly stable, we can conclude about the existence of phase space domains of long-term orbital stability.

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.

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.

Migration of planetary systems caused by the action of dissipative forces may lead the planets to be trapped in a resonance. In this work we study the conditions and the dynamics of such resonant trapping. Particularly, we are interested in finding out whether resonant capture ends up in a long-term stable planetary configuration. For two planet systems we associate the evolution of migration with the existence of families of periodic orbits in the phase space of the three-body problem. The family of circular periodic orbits exhibits a gap at the 2:1 resonance and an instability and bifurcation at the 3:1 resonance. These properties explain the high probability of 2:1 and 3:1 resonant capture at low eccentricities. Furthermore, we study the resonant capture of three-planet systems. We show that such a resonant capture is possible and can occur under particular conditions. Then, from the migration path of the system, stable three-planet configurations, either symmetric or asymmetric, can be determined.

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.

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

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.

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

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.

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.

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.

The acceleration of charged particles in the presence of a magnetic field and gravitational waves is under consideration. It is shown that the weak gravitational waves can cause the acceleration of low energy particles under appropriate conditions. Such conditions may be satisfied close to the source of the gravitational waves if the magnetized plasma is in a turbulent state.

In the framework of the planar restricted three body problem we study a considerable number of resonances associated to the Kuiper Belt dynamics and located between 30 and 48 a.u. Our study is based on the computation of resonant periodic orbits and their stability. Stable periodic orbits are surrounded by regular librations in phase space and in such domains the capture of trans-Neptunian object is possible. All the periodic orbits found are symmetric and there is evidence for the existence of asymmetric ones only in few cases. In the present work first, second and third order resonances are under consideration. In the planar circular case we found that most of the periodic orbits are stable. The families of periodic orbits are temporarily interrupted by collisions but they continue up to relatively large values of the Jacobi constant and highly eccentric regular motion exists for all cases. In the elliptic problem and for a particular eccentricity value of the primary bodies the periodic orbits are isolated. The corresponding families, where they belong to, bifurcate from specific periodic orbits of the circular problem and seem to continue up to the rectilinear problem. Both stable and unstable orbits are obtained for each case. In the elliptic problem the unstable orbits found are associated with narrow chaotic domains in phase space. The evolution of the orbits, which are located in such chaotic domains, seems to be practically regular and bounded for long time intervals.

We investigate the non-linear interaction of a strong Gravitational Wave with the plasma during the collapse of a massive magnetized star to form a black hole, or during the merging of neutron star binaries (central engine). We found that under certain conditions this coupling may result in an efficient energy space diffusion of particles. We suggest that the atmosphere created around the central engine is filled with 3-D magnetic neutral sheets (magnetic nulls). We demonstrate that the passage of strong pulses of Gravitational Waves through the magnetic neutral sheets accelerates electrons to very high energies. Superposition of many such short lived accelerators, embedded inside a turbulent plasma, may be the source for the observed impulsive short lived bursts. We conclude that in several astrophysical events, gravitational pulses may accelerate the tail of the ambient plasma to very high energies and become the driver for many types of astrophysical bursts.