The small ($\le$ 135 km mean radius) saturnian satellites are closely related to the rings and together they constitute a complex dynamical system where formation and destruction mechanisms compete against each other. We model the cratering process by Centaur objects on the small saturnian satellites, and compare our results with the observational crater counts obtained from the Voyager and Cassini missions. Using a theoretical model previously developed we calculate the crater production on these satellites considering two slopes of the size-frequency distribution (SFD) for the smaller objects of the Centaur population and compare our results with the available observations. In addition, we consider the case of catastrophic collisions between these satellites and Centaur objects and calculate the age of formation of those satellites that suffer one or more disruptions. In general we find that the observed crater distributions are best modeled by the crater size distribution corresponding to the $s_2 = 3.5$ index of the SFD of impactors with diameters smaller than 60 km. However, for crater diameters $D \lesssim 3-8$ km (which correspond to impactor diameters $d \sim 0.04 - 0.15$ km), the observed distributions become flatter and deviate from our results, which may evidence processes of erosion and/or crater saturation at small crater sizes or a possible break in the SFD of impactors at $d \sim 0.04 - 0.15$ km to a much shallower differential slope of $\sim -1.5$. Our results suggest that Pan, Daphnis, Atlas, Aegaeon, Methone, Anthe, Pallene, Calypso and Polydeuces suffered one or more catastrophic collisions over the age of the Solar System, the younger being associated to arcs with ages of $\sim10^8$ years. We have also calculated the satellites' surface ages, which indicate ongoing resurfacing processes.

We analyze the behavior of the argument of pericenter $\omega_2$ of an outer particle in the elliptical restricted three-body problem, focusing on the inverse Lidov-Kozai resonance. First, we calculate the contribution of the terms of quadrupole, octupole, and hexadecapolar order of the secular approximation of the potential to the outer particle's $\omega_2$ precession rate $(d\omega_2/d\tau)$. Then, we derive analytical criteria that determine the vanishing of the $\omega_2$ quadrupole precession rate $(d\omega_2/d\tau)_{\text{quad}}$ for different values of the inner perturber's eccentricity $e_1$. Finally, we use such analytical considerations and describe the behavior of $\omega_2$ of outer particles extracted from N-body simulations. Our analytical study indicates that the values of the inclination $i_2$ and the ascending node longitude $\Omega_2$ associated with the outer particle that vanish $(d\omega_2/d\tau)_{\text{quad}}$ strongly depend on the eccentricity $e_1$ of the inner perturber. In fact, if $e_1 <$ 0.25 ($>$ 0.40825), $(d\omega_2/d\tau)_{\text{quad}}$ is only vanished for particles whose $\Omega_2$ circulates (librates). For $e_1$ between 0.25 and 0.40825, $(d\omega_2/d\tau)_{\text{quad}}$ can be vanished for any particle for a suitable selection of pairs ($\Omega_2$, $i_2$). Our analysis of the N-body simulations shows that the inverse Lidov-Kozai resonance is possible for small, moderate and high values of $e_1$. Moreover, such a resonance produces distinctive features in the evolution of a particle in the ($\Omega_2$, $i_2$) plane. In fact, if $\omega_2$ librates and $\Omega_2$ circulates, the extremes of $i_2$ at $\Omega_2 =$ 90$^{\circ}$ and 270$^{\circ}$ do not reach the same value, while if $\omega_2$ and $\Omega_2$ librate, the evolutionary trajectory of the particle in the ($\Omega_2$, $i_2$) plane evidences an asymmetry respect to $i_2 =$ 90$^{\circ}$.

Recent studies have analyzed the dynamical evolution of outer small body populations under the effects of an eccentric inner massive perturber. Such outer reservoirs are composed of particles on prograde and retrograde orbits, as well as particles whose orbit flips from prograde to retrograde and back again showing a coupling between the inclination i and the ascending node longitude \Omega (Type-F particles). We analyze the role of the General Relativity (GR) on the dynamics of outer particles under the influence of an inner eccentric Jupiter-mass planet produced by a planetary scattering event. In particular, we study how the GR affects the dynamical evolution of the outer Type-F particles, which experience an eccentric Lidov-Kozai mechanism. We carry out N-body simulations with and without GR effects. When the GR is included, the extreme values of \Omega are obtained for retrograde inclinations, while the minimum and maximum inclinations allowed for Type-F particles increase in comparison with that derived without GR effects. According to this, if the GR is included in the simulations, the range of prograde (retrograde) inclinations of the libration region is reduced (increased) respect to that obtained in absence of GR. We find two new class of particles when the GR effects are included in the simulations. On the one hand, particles whose orbital plane flips from prograde to retrograde and back again without experiencing a coupling between i and \Omega. On the other hand, retrograde particles that show a strong coupling between i and \Omega. We infer that the GR may significantly modify the dynamical properties of the outer reservoirs that evolve under the effects of an eccentric inner perturber.

We analyze the process of planetary scattering around M0-type stars. To do this, we carry out N-body simulations with three Jupiter-mass planets close to their instability limit together with an outer planetesimal disk. This paper focuses on the analysis of systems in which a single Jupiter-mass planet survives after the dynamical instability event. The small body reservoirs show different dynamical behaviors. In fact, our simulations produce particles on prograde and retrograde orbits, as well as particles whose orbital plane flips from prograde to retrograde and back again along their evolution. Such particles are called "Type-F particles". We find strong correlations between the inclination $i$ and the ascending node longitude $\Omega$ of such particles. First, $\Omega$ librates around 90$^{\circ}$ or/and 270$^{\circ}$. This property is very important since it represents a necessary and sufficient condition for the flipping of an orbit. Moreover, the libration periods of $i$ and $\Omega$ are equal and they are out to phase by a quarter period. We also remark that the larger the libration amplitude of $i$, the larger the libration amplitude of $\Omega$. Finally, we analyze the initial conditions of Type-F particles of all our simulations immediately after the dynamical instability event, when a single Jupiter-mass planet survives in the system. We carry out this study with the goal to determine the parameter space that lead to the flipping of an orbit. Our results suggest that the orbit of a test particle can flip for any value of its initial eccentricity, although we found only two Type-F particles with initial inclinations $i <$ 17$^{\circ}$. Moreover, our study indicates that the minimum value of the inclination of the Type-F particles in a given system decreases with an increase in the eccentricity of the giant planet.

The secular approximation of the hierarchical three body systems has been proven to be very useful in addressing many astrophysical systems, from planets, stars to black holes. In such a system two objects are on a tight orbit, and the tertiary is on a much wider orbit. Here we study the dynamics of a system by taking the tertiary mass to zero and solve the hierarchical three body system up to the octupole level of approximation. We find a rich dynamics that the outer orbit undergoes due to gravitational perturbations from the inner binary. The nominal result of the precession of the nodes is mostly limited for the lowest order of approximation, however, when the octupole-level of approximation is introduced the system becomes chaotic, as expected, and the tertiary oscillates below and above 90deg, similarly to the non-test particle flip behavior (e.g., Naoz 2016). We provide the Hamiltonian of the system and investigate the dynamics of the system from the quadrupole to the octupole level of approximations. We also analyze the chaotic and quasi-periodic orbital evolution by studying the surfaces of sections. Furthermore, including general relativity, we show case the long term evolution of individual debris disk particles under the influence of a far away interior eccentric planet. We show that this dynamics can naturally result in retrograde objects and a puffy disk after a long timescale evolution (few Gyr) for initially aligned configuration.

The observations of the surfaces of the mid sized Saturnian satellites made by Cassini Huygens mission have shown a variety of features that allows study of the processes that took place and are taking place on those worlds. Research of the Saturnian satellite surfaces has clear implications for Saturn history and surroundings. In a recent paper, the production of craters on the mid sized Saturnian satellites by Centaur objects was calculated considering the current Solar System. We have compared our results with crater counts from Cassini images and we have noted that the number of observed small craters is less than our calculated number. In this paper we estimate the age of the surface for each observed terrain on each mid sized satellite of Saturn. We have noticed that since there are less observed small craters than calculated (except on Iapetus), this results in younger ages. This could be the result of efficient endogenous or exogenous process(es) for erasing small craters and or crater saturation at those sizes. The size limit from which the observed number of smaller craters is less than the calculated is different for each satellite, possibly indicating processes that are unique to each, but other potential common explanations would be crater saturation and or deposition of E ring particles. These processes are also suggested by the findings that the smaller craters are being preferentially removed, and the erasure process is gradual. On Enceladus, only mid and high latitude plains have remnants of old terrains; the other regions could be young; the regions near the South Polar Terrain could be as young as 50 Myr old. On the contrary for Iapetus, all the surface is old and it notably registers a primordial source of craters. As the crater size is decreased, it would be perceived to approach saturation until D less than 2 km craters, where saturation is complete.