Planets orbiting low-mass stars such as M dwarfs are now considered a cornerstone in the search for life-harbouring planets. GJ273 is a planetary system orbiting an M dwarf only 3.75 pc away, composed of two confirmed planets, GJ273b and GJ273c, and two promising candidates, GJ273d and GJ273e. Planet GJ273b resides in the habitable zone. Currently, due to a lack of observed planetary transits, only the minimum masses of the planets are known.Despite being an interesting system, the GJ273 planetary system is still poorly studied. We aim at precisely determine the physical parameters of the individual planets, in particular to break the mass--inclination degeneracy to accurately determine the mass of the planets. Moreover, we present thorough characterisation of planet GJ273b in terms of its potential habitability. We explored the planetary formation and hydration phases of GJ273 during the first 100 Myr. Then, we analysed the stability of the system. We also searched for regions which may harbour minor bodies such as an asteroid belt and Kuiper belt analogues. We found that the four-planet configuration of the system allows us to break the mass--inclination degeneracy with the following masses: $2.89\leq M_{\mathrm{b}}\leq3.03\,\mathrm{M}_\oplus$, $1.18\leq M_{\mathrm{c}}\leq1.24\,\mathrm{M}_\oplus$, $10.80\leq M_{\mathrm{d}}\leq11.35\,\mathrm{M}_\oplus$, and $9.30\leq M_{\mathrm{e}}\leq9.70\,\mathrm{M}_\oplus$. That is an Earth-mass planet, a super-Earth and two mini-Neptunes. Moreover, GJ273b is found to be an efficient water captor and GJ273c likely a dry planet. Several stable regions are predicted where minor bodies might reside. We comprehensively discuss the habitability of GJ273b.

The goal of this research is to study how the fragmentation of planetary embryos can affect the physical and dynamical properties of terrestrial planets around solar-type stars. Our work focuses on the formation and evolution of planets and water delivery in the Habitable Zone. We distinguish class A and class B HZ planets, which have an accretion seed initially located inside and beyond the snowline, respectively. We develop an Nbody integrator that incorporates fragmentation and hit-and-run collisions, which is called D3 N-body code. From this, we perform 46 numerical simulations of planetary accretion in systems that host two gaseous giants like Jupiter and Saturn. We compare two sets of 23 N-body simulations, one of which includes a realistic collisional treatment and the other one models all impacts as perfect mergers. The final masses of the HZ planets formed in runs with fragmentation are about 15\%-20\% smaller than those obtained without fragmentation. As for the class A HZ planets, those formed in simulations without fragmentation experience very significant increases in mass respect to their initial values, while the growth of those produced in runs with fragmentation is less relevant. We remark that the fragments play a secondary role in the masses of the class A HZ planets, providing less than 30\% of their final values. In runs without fragmentation, the final fraction of water of the class A HZ planets keeps the initial value since they do not accrete water-rich embryos. In runs with fragmentation, the final fraction of water of such planets strongly depends on the model used to distribute the water after each collision. The class B HZ planets do not show significant differences concerning their final water contents in runs with and without fragmentation. From this, we find that the collisional fragmentation is not a barrier to the survival of water worlds in the HZ.

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