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.

Aims. We developed a six-part collisional evolution model of the main asteroid belt (MB) and used it to study the contribution of the different regions of the MB to the near-earth asteroids (NEAs). Methods. We built a statistical code called ACDC that simulates the collisional evolution of the MB split into six regions (namely Inner, Middle, Pristine, Outer, Cybele and High-Inclination belts) according to the positions of the major resonances present there ($\nu_{6}$, 3:1J, 5:2J, 7:3J and 2:1J). We consider the Yarkovsky effect and the mentioned resonances as the main mechanism that removes asteroids from the different regions of the MB and delivers them to the NEA region. We calculated the evolution of the NEAs coming from the different source regions by considering the bodies delivered by the resonances and mean dynamical timescales in the NEA population. Results. Our model is in agreement with the major observational constraints associated with the MB, such as the size distributions of the different regions of the MB and the number of large asteroid families. It is also able to reproduce the observed NEAs with H < 16 and agrees with recent estimations for H < 20, but deviates for smaller sizes. We find that most sources make a significant contribution to the NEAs; however the Inner and Middle belts stand out as the most important source of NEAs followed by the Outer belt. The contributions of the Pristine and Cybele regions are minor. The High-Inclination belt is the source of only a fraction of the actual observed NEAs with high inclination, as there are dynamical processes in that region that enable asteroids to increase and decrease their inclinations.

Recent observational results show that very low mass stars and brown dwarfs are able to host close-in rocky planets. Low-mass stars are the most abundant stars in the Galaxy and the formation efficiency of their planetary systems is relevant in the computation of a global probability of finding Earth-like planets inside habitable zones. Tidal forces and relativistic effects are relevant in the latest dynamical evolution of planets around low-mass stars and their effect on the planetary formation efficiency still needs to be addressed. Our goal is to evaluate the impact of tidal forces and relativistic effects on the formation of rocky planets around a star close to the substellar mass limit, in terms of the resulting planetary architectures and its distribution according to the corresponding evolving habitable zone. Thus, we performed a set of $N$-body simulations spanning the first 100~Myr of the evolution of two systems composed respectively by 224 embryos with a total mass 0.25M$_\oplus$ and 74 embryos with a total mass 3 M$_\oplus$ around a central object of 0.08~M$_\odot$. For these two scenarios, we compared the planetary architectures that result from simulations that are purely gravitational with those from simulations that include the early contraction and spin-up of the central object, the distortions and dissipation tidal terms and general relativistic effects. We found that including these effects allows the formation and survival of a close-in population located in the habitable zone of the system. This means that both effects are relevant during the formation of rocky planets and their early evolution around stars close to the substellar mass limit, in particular when low-mass planetary embryos are involved.

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.

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

In the present research, we study the effects of a single giant planet in the dynamical evolution of water-rich embryos and planetesimals, located beyond the snow line of systems around Sun-like stars, in order to determine what kind of terrestrial-like planets could be formed in the habitable zone (hereafter HZ) of these systems. To do this, we carry out N-body simulations of planetary accretion, considering that the gas has been already dissipated from the disk and a single giant planet has been formed beyond the snow line of the system, at 3 au. We find that a giant planet with a value of mass between Saturn-mass and Jupiter-mass, represents a limit from which the amount of water-rich embryos that moves inward from beyond the snow line starts to decrease. From this, our research suggests that giant planets more massive than one Jupiter-mass become efficient dynamical barriers to inward-migrating water-rich embryos. Moreover, we infer that the number of these embryos that survive in the HZ significantly decreases for systems that host a giant planet more massive than one Jupiter-mass. This result has important consequences concerning the formation of terrestrial-like planets in the HZ with very high water contents and could provide a selection criteria in the search of potentially habitable exoplanets in systems that host a gaseous giant around solar-type stars.

This work is the second part of a project that attempts to analyze the formation of Solar system analogues (SSAs) from the gaseous to the post-gas phase, in a self-consistently way. In the first paper (PI) we presented our model of planet formation during the gaseous phase which provided us with embryo distributions, planetesimal surface density, eccentricity and inclination profiles of SSAs, considering different planetesimal sizes and type I migration rates at the time the gas dissipates. In this second work we focus on the late accretion stage of SSAs using the results obtained in PI as initial conditions to carry out N-body simulations. One of our interests is to analyze the formation of rocky planets and their final water contents within the habitable zone. Our results show that the formation of potentially habitable planets (PHPs) seems to be a common process in this kind of scenarios. However, the efficiency in forming PHPs is directly related to the size of the planetesimals. The smaller the planetesimals, the greater the efficiency in forming PHPs. We also analyze the sensitivity of our results to scenarios with type I migration rates and gap-opening giants, finding that both phenomena act in a similar way. These effects seem to favor the formation of PHPs for small planetesimal scenarios and to be detrimental for scenarios formed from big planetesimals. Finally, another interesting result is that the formation of water-rich PHPs seems to be more common than the formation of dry PHPs.

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.

Aims. We study the formation and water delivery of planets in the habitable zone (HZ) around solar-type stars. In particular, we study different dynamical environments that are defined by the most massive body in the system. Methods. First of all, a semi-analytical model was used to define the mass of the protoplanetary disks that produce each of the five dynamical scenarios of our research. Then, we made use of the same semi-analytical model to describe the evolution of embryos and planetesimals during the gaseous phase. Finally, we carried out N-body simulations of planetary accretion in order to analyze the formation and water delivery of planets in the HZ in the different dynamical environments. Results. Water worlds are efficiently formed in the HZ in different dynamical scenarios. In systems with a giant planet analog to Jupiter or Saturn around the snow line, super-Earths tend to migrate into the HZ from outside the snow line as a result of interactions with other embryos and accrete water only during the gaseous phase. In systems without giant planets, Earths and super-Earths with high water by mass contents can either be formed in situ in the HZ or migrate into it from outer regions, and water can be accreted during the gaseous phase and in collisions with water-rich embryos and planetesimals. Conclusions. The formation of planets in the HZ with very high water by mass contents seems to be a common process around Sun- like stars. Our research suggests that such planets are still very efficiently produced in different dynamical environments. Moreover, our study indicates that the formation of planets in the HZ with masses and water contents similar to those of Earth seems to be a rare process around solar-type stars in the systems under consideration.

Several observational works have shown the existence of Jupiter-mass planets covering a wide range of semi-major axes around Sun-like stars. We aim to analyse the planetary formation processes around Sun-like stars that host a Jupiter-mass planet at intermediate distances ranging from $\sim$1 au to 2 au. Our study focusses on the formation and evolution of terrestrial-like planets and water delivery in the habitable zone (HZ) of the system. Our goal is also to analyse the long-term dynamical stability of the resulting systems. A semi-analytic model was used to define the properties of a protoplanetary disk that produces a Jupiter-mass planet around the snow line, which is located at $\sim$2.7 au for a solar-mass star. Then, it was used to describe the evolution of embryos and planetesimals during the gaseous phase up to the formation of the Jupiter-mass planet, and we used the results as the initial conditions to carry out N-body simulations of planetary accretion. Our simulations produce three different classes of planets in the HZ: 'water worlds', with masses between 2.75 $M_{\oplus}$ and 3.57 $M_{\oplus}$ and water contents of 58% and 75% by mass, terrestrial-like planets, with masses ranging from 0.58 $M_{\oplus}$ to 3.8 $M_{\oplus}$ and water contents less than 1.2% by mass, and 'dry worlds', simulations of which show no water. A relevant result suggests the efficient coexistence in the HZ of a Jupiter-mass planet and a terrestrial-like planet with a percentage of water by mass comparable to the Earth. Moreover, our study indicates that these planetary systems are dynamically stable for at least 1 Gyr. Systems with a Jupiter-mass planet located at 1.5 au - 2 au around solar-type stars are of astrobiological interest. These systems are likely to harbour terrestrial-like planets in the HZ with a wide diversity of water contents.

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.

Population synthesis models of planetary systems developed during the last $\sim$15 years could reproduce several of the observables of the exoplanet population, and also allowed to constrain planetary formation models. We present our planet formation model, which calculates the evolution of a planetary system during the gaseous phase. The code incorporates relevant physical phenomena for the formation of a planetary system, like photoevaporation, planet migration, gas accretion, water delivery in embryos and planetesimals, a detailed study of the orbital evolution of the planetesimal population, and the treatment of the fusion between embryos, considering their atmospheres. The main goal of this work, unlike other works of planetary population synthesis, is to find suitable scenarios and physical parameters of the disc to form solar system analogs. We are specially interested in the final planet distributions, and in the final surface density, eccentricity and inclination profiles for the planetesimal population. These final distributions will be used as initial conditions for N-body simulations, to study the post-oligarchic formation in a second work. We then consider different formation scenarios, with different planetesimal sizes and different type I migration rates. We find that solar system analogs are favored in massive discs, with low type I migration rates, and small planetesimal sizes. Besides, those rocky planets within their habitables zones are dry when discs dissipate. At last, the final configurations of solar system analogs include information about the mass and semimajor-axis of the planets, water contents, and the properties of the planetesimal remnants.

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.

To study the terrestrial-type planet formation during the post oligarchic growth, the initial distributions of planetary embryos and planetesimals used in N-body simulations play an important role. Most of these studies typically use ad hoc initial distributions based on theoretical and numerical studies. We analyze the formation of planetary systems without gas giants around solar-type stars focusing on the sensitivity of the results to the particular initial distributions of planetesimals and embryos. The formation of terrestrial planets in the habitable zone (HZ) and their final water contents are topics of interest. We developed two different sets of N-body simulations from the same protoplanetary disk. The first set assumes ad hoc initial distributions for embryos and planetesimals and the second set obtains these distributions from the results of a semi-analytical model which simulates the evolution of the gaseous phase of the disk. Both sets form planets in the HZ. Ad hoc initial conditions form planets in the HZ with masses from $0.66M_{\oplus}$ to $2.27M_{\oplus}$. More realistic initial conditions obtained from a semi-analytical model, form planets with masses between $1.18M_{\oplus}$ and $2.21M_{\oplus}$. Both sets form planets in the HZ with water contents between 4.5% and 39.48% by mass. Those planets with the highest water contents respect to those with the lowest, present differences regarding the sources of water supply. We suggest that the number of planets in the HZ is not sensitive to the particular initial distribution of embryos and planetesimals and thus, the results are globally similar between both sets. However, the main differences are associated to the accretion history of the planets in the HZ. These discrepancies have a direct impact in the accretion of water-rich material and in the physical characteristics of the resulting planets.

Models of planet formation are mainly focused on the accretion and dynamical processes of the planets, neglecting their chemical composition. In this work, we calculate the condensation sequence of the different chemical elements for a low-mass protoplanetary disk around a solar-type star. We incorporate this sequence of chemical elements (refractory and volatile elements) in our semi-analytical model of planet formation which calculates the formation of a planetary system during its gaseous phase. The results of the semi-analytical model (final distributions of embryos and planetesimals) are used as initial conditions to develope N-body simulations that compute the post-oligarchic formation of terrestrial-type planets. The results of our simulations show that the chemical composition of the planets that remain in the habitable zone has similar characteristics to the chemical composition of the Earth. However, exist differences that can be associated to the dynamical environment in which they were formed.

In general, most of the studies of terrestrial-type planet formation typically use ad hoc initial conditions. In this work we improved the initial conditions described in Ronco & de El\'ia (2014) starting with a semi-analytical model wich simulates the evolution of the protoplanetary disk during the gas phase. The results of the semi-analytical model are then used as initial conditions for the N-body simulations. We show that the planetary systems considered are not sensitive to the particular initial distribution of embryos and planetesimals and thus, the results are globally similar to those found in the previous work.

In the standard model of core accretion, the cores of the giant planets form by the accretion of planetesimals. In this scenario, the evolution of the planetesimal population plays an important role in the formation of massive cores. Recently, we studied the role of planetesimal fragmentation in the in situ formation of a giant planet. However, the exchange of angular momentum between the planet and the gaseous disk causes the migration of the planet in the disk. In this new work, we incorporate the migration of the planet and globally study the role of planet migration in the formation of a massive core when the population of planetesimals evolves by planet accretion, migration due to the nebular drag, and fragmentation due to planetesimal collisions.

Several studies, observational and theoretical, suggest that planetary systems with only rocky planets should be the most common in the Universe. We study the diversity of planetary systems that might form around Sun-like stars in low-mass disks without giant planets. We focus on the formation process of terrestrial planets in the habitable zone (HZ) and analyze their water contents with the goal to determine systems of astrobiological interest. Besides, we study the formation of planets on wide orbits because they can be detected with the microlensing technique. N-body simulations of high resolution (embryos + planetesimals) are developed for a wide range of surface density profiles. The surface density profile combines a power law to the inside of the disk of the form r^{-gamma}, with an exponential decay to the outside. We adopt a disk of 0.03M_sun and values of gamma = 0.5, 1 and 1.5. All our simulations form planets in the HZ with different masses and final water contents depending on the 3 profiles. For gamma = 0.5, we produce 3 planets in the HZ with masses between 0.03 M_e to 0.1 M_e and water contents between 0.2 and 16 Earth oceans. For gamma = 1, 3 planets form in the HZ with masses between 0.18 M_e and 0.52 M_e and water contents from 34 to 167 Earth oceans. For gamma = 1.5, we find 4 planets in the HZ with masses from 0.66 M_e to 2.21 M_e and water contents between 192 and 2326 Earth oceans. This profile shows distinctive results because it is the only one of those studied here that leads to the formation of water worlds. Since planetary systems with gamma = 1 and 1.5 present planets in the HZ with suitable masses to retain a long-live atmosphere and to maintain plate tectonics, they seem to be the most outstanding candidates to be potentially habitable. Particularly, these systems form Earths and Super-Earths near the snow line which can be discovered by microlensing.

In the standard scenario of planet formation, terrestrial planets and the cores of the giant planets are formed by accretion of planetesimals. As planetary embryos grow the planetesimal velocity dispersion increases due to gravitational excitations produced by embryos. The increase of planetesimal relative velocities causes the fragmentation of them due to mutual collisions. We study the role of planetesimal fragmentation on giant planet formation. We analyze how planetesimal fragmentation modifies the growth of giant planet's cores for a wide range of planetesimal sizes and disk masses. We incorporate a model of planetesimal fragmentation into our model of in situ giant planet formation. We calculate the evolution of the solid surface density (planetesimals plus fragments) due to the accretion by the planet, migration and fragmentation. The incorporation of planetesimal fragmentation significantly modifies the process of planetary formation. If most of the mass loss in planetesimal collisions is distributed in the smaller fragments, planetesimal fragmentation inhibits the growth of the embryo for initial planetesimals of radii lower than 10 km. Only for initial planetesimals of 100 km of radius, and disks greater than 0.06 Msun, embryos achieve masses greater than the mass of the Earth. However, even for such big planetesimals and massive disks, planetesimal fragmentation induces the quickly formation of massive cores only if most of the mass loss in planetesimal collisions is distributed in the bigger fragments. Planetesimal fragmentation seems to play an important role in giant planet formation. The way in which the mass loss in planetesimal collisions is distributed leads to different results, inhibiting or favoring the formation of massive cores.

We study the impactor flux and cratering on Ceres and Vesta caused by the collisional and dynamical evolution of the asteroid Main Belt. We develop a statistical code based on a well-tested model for the simultaneous evolution of the Main Belt and NEA size distributions. This code includes catastrophic collisions and noncollisional removal processes such as the Yarkovsky effect and the orbital resonances. The model assumes that the dynamical depletion of the early Main Belt was very strong, and owing to that, most Main Belt comminution occurred when its dynamical structure was similar to the present one. Our results indicate that the number of D > 1 km Main Belt asteroids striking Ceres and Vesta over the Solar System history are approximately 4 600 and 1 100 respectively. The largest Main Belt asteroids expected to have impacted Ceres and Vesta had diameters of 71.7 km and 21.1 km. The number of D > 0.1 km craters on Ceres is \sim 3.4 \times 10^8 and 6.2 \times 10^7 on Vesta. The number of craters with D > 100 km are 47 on Ceres and 8 on Vesta. Our study indicates that the D = 460 km crater observed on Vesta had to be formed by the impact of a D \sim 66.2 km projectile, which has a probability of occurr \sim 30% over the Solar System history. If significant discrepancies between our results about the cratering on Ceres and Vesta and data obtained from the Dawn Mission were found, they should be linked to a higher degree of collisional evolution during the early Main Belt and/or the existence of the late heavy bombardment. An increase in the collisional activity in the early phase may be provided for an initial configuration of the giant planets consistent with, for example, the Nice model. From this, the Dawn Mission would be able to give us clues about the initial configuration of the early Solar System and its subsequent dynamical evolution.

We study the impactor flux and cratering on Pluto and Charon due to the collisional evolution of Plutinos. Plutinos are those trans-Neptunian objects located at 39.5 AU, in the 3:2 mean motion resonance with Neptune. To do this, we develop a statistical code that includes catastrophic collisions and cratering events, and takes into account the stability and instability zones of the 3:2 mean motion resonance with Neptune. We proposes different initial populations that account for the uncertainty in the size distribution of Plutinos at small sizes. Depending on the initial population, our results indicate the following. The number of D > 1 km Plutinos streaking Pluto over 3.5 Gyr is between 1271 and 5552. For Charon, the number of D > 1 km Plutino impactors is between 354 and 1545. The number of D > 1 km craters on Pluto produced by Plutinos during the last 3.5 Gyr is between 43076 and 113879. For Charon, the number of D > 1 km craters is between 20351 and 50688.The largest Plutino impactor onto Pluto has a diameter between 17 and 23 km, which produces a craterwith a diameter of 31 - 39 km. The largest Plutino impactor onto Charon has a diameter between 10 and 15 km, which produces a crater with a diameter of 24 - 33 km. We test if 2 Pluto-sized objects are assumed in the 3:2 Neptune resonance, then the total number of Plutino impactors both onto Pluto as Charon with diameters D > 1 km is a factor of 1.6 - 1.8 larger if considering 1 Pluto-sized object. Given the dynamical structure of the trans-Neptunian region, it is necessary to study in detail the contribution of all the potential sources of impactors on the Pluto-Charon system, to obtain the main contributor and the whole production of craters. Then, we will be able to contrast those studies with observations which will help us to understand the geological processes and history of the surface of those worlds.