We propose a semianalytical method for the calculation of widths, libration centers and small amplitude libration periods of the mean motion resonances k_p:k in the framework of the circular restricted three body problem valid for arbitrary eccentricities and inclinations. Applying the model to the trans Neptunian region (TNR) we obtain several atlas of resonances between 30 and 100 au showing their domain in the plane (a,e) for different orbital inclinations. The resonance width may change substantially when varying the argument of the perihelion of the resonant object and in order to take into account these variations we introduce the concept of resonance fragility. Resonances 1:k and 2:k are the widest, strongest, most isolated ones and with lower fragility for all interval of inclinations and eccentricities. We discuss about the existence of high k_p:k resonances. We analyze the distribution of the resonant populations inside resonances 1:1, 2:3, 3:5, 4:7, 1:2 and 2:5. We found that the populations are in general located near the regions of the space (e,i) where the resonances are wider and less fragile with the notable exception of the population inside the resonance 4:7 and in a lesser extent the population inside 3:5 which are shifted to lower eccentricities.

Applying the technique of dynamical maps we study the orbital stability of test particles in the Solar System in the space (a,e,i) defined by 0.1<a<38 au, 0<e<0.9 and 0<i<180 identifying the unstable and stable regions. We find stable niches where small bodies can survive even for very high eccentricities. Mean motion resonances play a fundamental role providing stability against the planetary perturbations specially for high inclination orbits. A stability stripe around i=150 is present all along the Solar System. We found that the population of objects with semimajor axes between 10 and 30 au is evolving inside a highly unstable region according to our maps. For the inner Solar System we found that the region between the Hildas and Jupiter is more stable for high eccentricity orbits than for low eccentricity ones.

Jupiter Trojan asteroids are located around L4 and L5 Lagrangian points on relatively stable orbits, in 1:1 MMR with Jupiter. However, not all of them lie in orbits that remain stable over the age of the Solar System. Unstable zones allow some Trojans to escape in time scales shorter than the Solar System age. This may contribute to populate other small body populations. In this paper, we study this process by performing long-term numerical simulations of the observed Trojans, focusing on the trajectories of those that leave the resonance. The orbits of current Trojans are taken as initial conditions and their evolution is followed under the gravitational action of the Sun and the planets. We find the rate of escape of Trojans from L5, ~1.1 times greater than from L4. The majority of escaped Trojans have encounters with Jupiter although they have encounters with the other planets too. Almost all escaped Trojans reach the comet zone, ~90% cross the Centaur zone and only L4 Trojans reach the transneptunian zone. Considering the real asymmetry between L4 and L5, we show that 18 L4 Trojans and 14 L5 Trojans with diameter D > 1 km are ejected from the resonance every Myr. The contribution of the escaped Trojans to other minor body populations would be negligible, being the contribution from L4 and L5 to JFCs and no-JFCs almost the same, and the L4 contribution to Centaurs and TNOs, orders of magnitude greater than that of L5. Considering the collisional removal, besides the dynamical one, and assuming that Trojans that escape due to collisions follow the same dynamical behavior that the ones removed by dynamics, we would have a minor contribution of Trojans to comets and Centaurs. However, there would be some specific regions were escaped Trojans could be important such as ACOs, Encke-type comets, S-L 9-type impacts on Jupiter and NEOs.

In the framework of the circular restricted three body problem we show that the numerically computed strength SR(e,i,w) is a good indicator of the strength and width of the mean-motion resonances in the full space (e,i,w). We present a survey of strengths in the space (e,i) for typical interior and exterior resonances. The resonance strength is highly dependent on (e,i,w) except for exterior resonances of the type 1:k for which the dependence with (i,w) is softer. Such resonances are thus strong even for retrograde orbits. All other resonances are weaker at very-high eccentricities for w ~ 90 or 270 and 60 < i < 120. We explore the resonance structure in the space (a,i) by means of dynamical maps and we find structures similar to those of space (a,e).

In this work we pose the possibility that, at an early stage, the migration of a proto--giant planet caused by the presence of a gaseous circumstellar disk could explain the continuous feeding of small bodies into its orbit. Particularly, we study the probability of capture and permanence in co--orbital resonance of these small bodies, as planets of diverse masses migrate by interaction with the gaseous disk, and the drag induced by this disk dissipates energy from these small objects, making capture more likely. Also, we study the relevance of the circumplanetary disk, a structure formed closely around the planet where the gas density is enhanced, in the process of capture. It is of great interest for us to study the capture of small bodies in 1:1 resonance because it could account for the origin of the Trojan population, which has been proposed \citep{2011Icar..215..669K} as a mechanism of quasi-satellites and irregular satellites capture.

We numerically study the dynamical evolution of observed samples of active and inactive Centaurs and clones that reach the Jupiter-Saturn region. Our aim is to compare the evolution between active and inactive Centaurs, their end states and their transfer to Jupiter family comets and Halley-type comets. We find that the median lifetime of inactive Centaurs is about twice longer than that for active Centaurs, suggesting that activity is related to the residence time in the region. This view is strengthened by the observation that high-inclination and retrograde Centaurs (Tisserand parameters with respect to Jupiter $T_J < 2$) which have the longest median dynamical lifetime ($=1.37 \times 10^6$ yr) are all inactive. We also find that the perihelion distances of some active, comet-like Centaurs have experienced drastic drops of a few au in the recent past ($\sim 10^2-10^3$ yr), while such drops are not found among inactive Centaurs. Inactive Centaurs with $T_J \lsim 2.5$ usually evolve to Halley-type comets, whereas inactive Centaurs with $T_J \gsim 2.5$ and active Centaurs (that also have $T_J \gsim 2.5$) evolve almost always to Jupiter family comets and very seldom to Halley type comets. Inactive Centaurs are also more prone to end up as sungrazers, and both inactive and active Centaurs transit through different mean motion resonances (generally with Jupiter) during their evolution.

Mean motion resonances play a fundamental role in the dynamics of the small bodies of the Solar System. The last decades of the 20th century gave us a detailed description of the dynamics as well as the process of capture of small bodies in coplanar or small inclination resonant orbits. More recently, semianalytical or numerical methods allowed us to explore the behavior of resonant motions for arbitrary inclination orbits. The emerging dynamics is very rich, including large orbital changes due to secular effects inside mean motion resonances. The process of capture in highly inclined or retrograde resonant orbits was addressed showing that the capture in retrograde resonances is more efficient than in direct ones. A new terminology appeared in order to characterize the properties of the resonances. Numerical explorations in the transneptunian region showed the relevance and the particular dynamics of the exterior resonances with Neptune which can account for some of the known high perihelion orbits in the scattered disk. Moreover, several asteroids evolving in resonance with planets other than Jupiter or Neptune were found and a large number of asteroids in three-body resonances were identified.

We analyze a sample of 73 old long-period comets (LPCs) (orbital periods $200 < P < 1000$ yr) with perihelion distances $q < 2.5$ au, discovered in the period 1850-2014. We cloned the observed comets and also added fictitious LPCs with perihelia in the Jupiter's zone. We consider both a purely dynamical evolution and a physico-dynamical one with different physical lifetimes. We can fit the computed energy distribution of comets with $q < 1.3$ au to the observed one only within the energy range $0.01 < x < 0.04$ au$^{-1}$ (or periods $125 < P < 1000$ yr), where the "energy" is taken as the inverse of the semimajor axis $a$, namely $x \equiv 1/a$. The best results are obtained for physical lifetimes of about 200-300 revolutions (for a comet with a standard $q = 1$ au). We find that neither a purely dynamical evolution, nor a physico-dynamical one can reproduce the long tail of larger binding energies ($x \gsim 0.04$ au$^{-1}$) that correspond to most Halley-type comets (HTCs) and Jupiter-family comets. We conclude that most HTCs are not the end states of the evolution of LPCs, but come from a different source, a flattened one that we identify with the Centaurs that are scattered to the inner planetary region from the trans-Neptunian belt. These results also show that the boundary between LPCs and HTCs should be located at an energy $x \sim 0.04$ au$^{-1}$ ($P \sim 125$ yr), rather than the conventional classical boundary at $P = 200$ yr.

We propose a semianalytical method to compute the strengths on each of the three massive bodies participating in a three body mean motion resonance (3BR). Applying this method we explore the dependence of the strength on the masses, the orbital parameters and the order of the resonance and we compare with previous studies. We confirm that for low eccentricity low inclination orbits zero order resonances are the strongest ones; but for excited orbits higher order 3BRs become also dynamically relevant. By means of numerical integrations and the construction of dynamical maps we check some of the predictions of the method. We numerically explore the possibility of a planetary system to be trapped in a 3BR due to a migrating scenario. Our results suggest that capture in a chain of two body resonances is more probable than a capture in a pure 3BR. When a system is locked in a 3BR and one of the planets is forced to migrate the other two can react migrating in different directions. We exemplify studying the case of the Galilean satellites where we show the relevance of the different resonances acting on the three innermost satellites.

We analyze a sample of 139 near-Earth asteroids (NEAs), defined as those that reach perihelion distances $q < 1.3$ au, and that also fulfill the conditions of approaching or crossing Jupiter's orbit (aphelion distances $Q > 4.8$ au), having Tisserand parameters $2 < T < 3$ and orbital periods $P < 20$ yr. In order to compare the dynamics, we also analyze a sample of 42 Jupiter family comets (JFCs) in near-Earth orbits, i.e. with $q < 1.3$ au. We integrated the orbits of these two samples for $10^4$ yr in the past and in the future. We find that the great majority of the NEAs move on stable orbits during the considered period, and that a large proportion of them are in one of the main mean motion resonances with Jupiter, in particular the 2:1. We find a strong coupling between the perihelion distance and the inclination in the motion of most NEAs, due to Kozai mechanism, that generates many sungrazers. On the other hand, most JFCs are found to move on very unstable orbits, showing large variations in their perihelion distances in the last few $10^2 - 10^3$ yr, which suggests a rather recent capture in their current near-Earth orbits. Even though most NEAs of our sample move in typical 'asteroidal' orbits, we detect a small group of NEAs whose orbits are highly unstable, resembling those of the JFCs. These are: 1997 SE5, 2000 DN1, 2001 XQ, 2002 GJ8, 2002 RN38, 2003 CC11, 2003 WY25, 2009 CR2, and 2011 OL51. These objects might be inactive comets, and indeed 2003 WY25 has been associated with comet Blanpain, and it is now designed as comet 289P/Blanpain. Under the assumption that these objects are inactive comets, we can set an upper limit of $\sim 0.17$ to the fraction of NEAs with $Q > 4.8$ au of cometary origin, but it could be even lower if the NEAs in unstable orbits listed before turn out to be {\it bona fide} asteroids from the main belt.

We present a numerical method to estimate the strengths of arbitrary three body mean motion resonances between two planets in circular coplanar orbits and a massless particle in an arbitrary orbit. This method allows us to obtain an atlas of the three body resonances in the Solar System showing where are located and how strong are thousands of resonances involving all the planets from 0 to 1000 au. This atlas confirms the dynamical relevance of the three body resonances involving Jupiter and Saturn in the asteroid belt but also shows the existence of a family of relatively strong three body resonances involving Uranus and Neptune in the far Trans-Neptunian region and relatively strong resonances involving terrestrial and jovian planets in the inner planetary system. We calculate the density of relevant resonances along the Solar System resulting that the main asteroid belt is located in a region of the planetary system with the lowest density of three body resonances. The method also allows the location of the equilibrium points showing the existence of asymmetric librations (sigma different from 0 or 180 degrees). We obtain the functional dependence of the resonance's strength with the order of the resonance and the eccentricity and inclination of the particle's orbit. We identify some objects evolving in or very close to three body resonances with Earth-Jupiter, Saturn-Neptune and Uranus-Neptune apart from Jupiter-Saturn, in particular the NEA 2009 SJ18 is evolving in the resonance 1-1E-1J and the centaur 10199 Chariklo is evolving under the influence of the resonance 5-2S-2N.

We study the Kozai dynamics affecting the orbital evolution of transneptunian objects being captured or not in MMR with Neptune. We provide energy level maps of the type ({\omega},q) describing the possible orbital paths from Neptune up to semimajor axis of hundreds of AU. The dynamics for non resonant TNOs with perihelion distances, q, outside Neptune's orbit, a_N, is quite different from the dynamics of TNOs with q < a_N, already studied in previous works. While for the last case there are stable equilibrium points at {\omega} = 0\circ, 90\circ, 180\circ and 270\circ in a wide range of orbital inclinations, for the former case it appears a family of stable equilibrium points only at a specific value of the orbital inclination, i \sim 62\circ, that we call critical inclination. We show this family of equilibrium points is generated by a mechanism analogue to which drives the dynamics of an artificial satellite perturbed by an oblate planet. The planetary system also generates an oscillation in the longitude of the perihelion of the TNOs with i \sim 46\circ, being Eris a paradigmatic case. We discuss how the resonant condition with Neptune modify the energy level curves and the location of equilibrium points. The asymmetric librations of resonances of the type 1:N generate a distortion in the energy level curves and in the resulting location of the equilibrium points in the phase space ({\omega},q). We study the effect on the Kozai dynamics due to the diffusion process in "a" that occurs in the Scattered Disk. We show that a minimum orbital inclination is required to allow substantial variations in perihelion distances once the object is captured in MMR and that minimum inclination is greater for greater semimajor axis.

At present, approximately 1500 asteroids are known to evolve inside or sticked to the exterior 1:2 resonance with Mars at a = 2.418 AU, being (142) Polana the largest member of this group. The effect of the forced secular modes superposed to the resonance gives rise to a complex dynamical evolution. Chaotic diffusion, collisions, close encounters with massive asteroids and mainly orbital migration due to the Yarkovsky effect generate continuous captures to and losses from the resonance, with a fraction of asteroids remaining captured over long time scales and generating a concentration in the semimajor axis distribution that exceeds by 20% the population of background asteroids. The Yarkovsky effect induces different dynamics according to the asteroid size, producing an excess of small asteroids inside the resonance. The evolution in the resonance generates a signature on the orbits, mainly in eccentricity, that depends on the time the asteroid remains captured inside the resonance and on the magnitude of the Yarkovsky effect. The greater the asteroids, the larger the time they remain captured in the resonance, allowing greater diffusion in eccentricity and inclination. The resonance generates a discontinuity and mixing in the space of proper elements producing misidentification of dynamical family members, mainly for Vesta and Nysa-Polana families. The half-life of resonant asteroids large enough for not being affected by the Yarkovsky effect is about 1 Gyr. From the point of view of taxonomic classes, the resonant population does not differ from the background population and the excess of small asteroids is confirmed.

Using Gauss' averaged equations, we compute the secular relativistic effects generated by the Sun on the argument of the perihelion and the mean anomaly of an orbit. Then we test different alternative simpler models that have been proposed to reproduce the secular relativistic effects in the orbital elements. Generally, models introduce artificial perturbations that are velocity-independent but that depend on the heliocentric distance. If these perturbations are set as an impulse in a constant timestep integrator, when the particle approaches perihelion the generated impulse could be very strong and badly sampled, originating a spurious orbital evolution. In order to overcome this setback, we propose two new models based on a constant, distance-independent, perturbation. With these models we obtain the correct secular drift in the argument of perihelion and the expected secular orbital evolution is reproduced. We also discuss with some detail the secular effect generated on the mean anomaly by different models.

There is a growing population of relativistically relevant minor bodies in the Solar System and a growing population of massive extrasolar planets with orbits very close to the central star where relativistic effects should have some signature. Our purpose is to review how general relativity affects the orbital dynamics of the planetary systems and to define a suitable relativistic correction for Solar System orbital studies when only point masses are considered. Using relativistic formulae for the N body problem suited for a planetary system given in the literature we present a series of numerical orbital integrations designed to test the relevance of the effects due to the general theory of relativity in the case of our Solar System. Comparison between different algorithms for accounting for the relativistic corrections are performed. Relativistic effects generated by the Sun or by the central star are the most relevant ones and produce evident modifications in the secular dynamics of the inner Solar System. The Kozai mechanism, for example, is modified due to the relativistic effects on the argument of the perihelion. Relativistic effects generated by planets instead are of very low relevance but detectable in numerical simulations.

The characteristics of the resonant disturbing function for an asteroid perturbed by a planet in circular orbit are discussed. The location of the libration centers and their dependence with the orbital elements of the resonant orbit are analyzed. A proposed numerical method (Gallardo 2006a) for computing the strengths of the resonances is revised and applied to the region of the main belt of asteroids showing the relevance of several mean motion resonances (MMR) with several planets.