The ring system and small satellites of Neptune were discovered during Voyager 2 flyby in 1989 (Smith et al.1989). In this work we analyse the diffusion maps which can give an overview of the system. As a result we found the width of unstable and stable regions close to each satellite. The innermost Galle ring, which is further from the satellites, is located in a stable region, while Lassel ring (W = 4000 km) has its inner border in a stable region depending on its eccentricity. The same happens to the Le Verrier and Adams rings , they are stable for small values of the eccentricity. They can survive to the close satellites perturbation only for values of e < 0.012. When the solar radiation force is taken into account the rings composed by 1{$\mu$}m sized particles have a lifetime of about $10^4$ years while larger particles (10{$\mu$}m in radius) can survive up to 10 5 years. The satellites Naiad, Thalassa and Despina can help replenish the lost particles of the Le Verrier, Arago and Lassel rings, while the ejecta produced by Galatea, Larissa and Proteus do not have enough velocity to escape from the satellite gravity.

We explore the long-term evolution of a bias-free orbital representation of the cometary nuclei (with diameters above 2 km) of the Kuiper belt, using the so-called L7 synthetic model from CFEPS, which consists of three dynamical sub-populations: the Classical, the Resonant, and the Scattering. The dynamical evolution of belt particles is studied under the gravitational influence of the Sun and the four giant planets, as well as of the 34 largest known trans-Neptunian objects (TNOs with $H_V < 4$). Here we indistinctly call Dwarf Planets (DPs) to the full sample of 34 large TNOs. Over a 1 Gyr timescale, we analyze the secular influence of the DPs over Kuiper belt disk particles and their contribution to the injection rate of new visible Jupiter Family Comets (JFCs). We found that DPs globally increase the number of JFCs by 12.6% when compared with the comets produced by the giant planets alone. When considering each population separately, we find that the increment produced by DPs is 17%, 12%, and 3% for the Classical, Resonant, and Scattering populations, respectively. Given the rate of escapes from the Kuiper belt, we find upper limits to the number of objects in each population required to maintain the JFCs in steady-state; the results are $55.9\times10^6$, $78.5\times10^6$, and $274.3\times10^6$ for the Scattering, Resonant, and Classical populations, respectively. Finally, we find that the Plutinos are the most important source of comets which were originally in a resonant configuration, where the presence of Pluto alone enhances by 10% the number of JFCs.

In this work we explore a new dynamical path for the delivery of low-inclination comets. In a configuration formed by an interior giant planet and an exterior massive debris disk, where the mass is accounted for by the 50 largest objects in the disk, the strongest mean motion resonances of the giant, located along the belt, are replenished with new material (test particles) due to the influence of the 50 massive objects. Once in resonance, slow chaotic diffusion stirs the orbital elements of the cometary nuclei enough to encounter the giant and to be scattered by it. When the disk is massive enough, both resonant and non-resonant particles are stirred quickly to encounter the giant and form an scattered disk component, greatly increasing the rate for the delivery of cometary material to the inner part of the system. This mechanism is applicable both to the solar system and extrasolar systems in general. Preliminary results, using a disk as massive as the classical Kuiper belt, indicate that the mechanism here proposed can account for about a tenth of the required injection rate to maintain the population of ecliptic comets in steady state. In a more massive belt of 0.25 M$_\oplus$, an estimated rate of around 0.6 new comets per year is found. Such a high rate would pose a serious risk for the habitability of rocky interior planets, yet would resemble the late heavy bombardment that was present in the early solar system.

Aegaeon, Methone, Anthe, and Pallene are four Saturnian small moons, discovered by the Cassini spacecraft. Although their orbital characterization has been carried on by a number of authors, their long-term evolution has not been studied in detail so far. In this work, we numerically explore the long-term evolution, up to $10^5$ yr, of the small moons in a system formed by an oblate Saturn and the five largest moons close to the region: Janus, Epimetheus, Mimas, Enceladus, and Tethys. By using frequency analysis we determined the stability of the small moons and characterize, through diffusion maps, the dynamical behavior of a wide region of geometric phase space, $a$ vs $e$, surrounding them. Those maps could shed light on the possible initial number of small bodies close to Mimas, and help to better understand the dynamical origin of the small satellites. We found that the four small moons are long-term stable and no mark of chaos is found for them. Aegaeon, Methone, and Anthe could remain unaltered for at least $\sim0.5$Myr, given the current configuration of the system. They remain well-trapped in the corotation eccentricity resonances with Mimas in which they currently librate. However, perturbations from nearby resonances, such as Lindblad eccentricity resonances with Mimas, seem responsible for largest variations observed for Methone and Anthe. Pallene remains in a non-resonant orbit and it is the more stable, at least for 64 Myr. Nonetheless, it is affected by a quasi-resonance with Mimas, which induces long-term orbital oscillations of its eccentricity and inclination.

We have explored the evolution of a cold debris disk under the gravitational influence of dwarf planet sized objects (DPs), both in the presence and absence of an interior giant planet. Through detailed long-term numerical simulations, we demonstrate that, when the giant planet is not present, DPs can stir the eccentricities and inclinations of disk particles, in linear proportion to the total mass of the DPs; on the other hand, when the giant planet is included in the simulations, the stirring is approximately proportional to the mass squared. This creates two regimes: below a disk mass threshold (defined by the total mass of DPs), the giant planet acts as a stabilizing agent of the orbits of cometary nucleii, diminishing the effect of the scatterers; above the threshold, the giant contributes to the dispersion of the particles.

With the use of long-term numerical simulations, we study the evolution and orbital behavior of cometary nuclei in cold Kuiper belt-like debris disks under the gravitational influence of dwarf planets (DPs); we carry out these simulations with and without the presence of a Neptune-like giant planet. This exploratory study shows that in the absence of a giant planet, 10 DPs are enough to induce strong radial and vertical heating on the orbits of belt particles. On the other hand, the presence of a giant planet close to the debris disk, acts as a stability agent reducing the radial and vertical heating. With enough DPs, even in the presence of a Neptune-like giant planet some radial heating remains; this heating grows steadily, re-filling resonances otherwise empty of cometary nuclei. Specifically for the solar system, this secular process seems to be able to provide material that, through resonant chaotic diffusion, increase the rate of new comets spiraling into the inner planetary system, but only if more than the $\sim10$ known DP sized objects exist in the trans-Neptunian region.

The orbital elements of comet Halley are known to a very high precision, suggesting that the calculation of its future dynamical evolution is straightforward. In this paper we seek to characterize the chaotic nature of the present day orbit of comet Halley and to quantify the timescale over which its motion can be predicted confidently. In addition, we attempt to determine the timescale over which its present day orbit will remain stable. Numerical simulations of the dynamics of test particles in orbits similar to that of comet Halley are carried out with the Mercury 6.2 code. On the basis of these we construct survival time maps to assess the absolute stability of Halley's orbit, frequency analysis maps, to study the variability of the orbit and we calculate the Lyapunov exponent for the orbit for variations in initial conditions at the level of the present day uncertainties in our knowledge of its orbital parameters. On the basis of our calculations of the Lyapunov exponent for comet Halley, the chaotic nature of its motion is demonstrated. The e-folding timescale for the divergence of initially very similar orbits is approximately 70 years. The sensitivity of the dynamics on initial conditions is also evident in the self-similarity character of the survival time and frequency analysis maps in the vicinity of Halley's orbit, which indicates that, on average, it is unstable on a timescale of hundreds of thousands of years. The chaotic nature of Halley's present day orbit implies that a precise determination of its motion, at the level of the present day observational uncertainty, is difficult to predict on a timescale of approximately 100 years. Furthermore, we also find that the ejection of Halley from the solar system or its collision with another body could occur on a timescale as short as 10,000 years.