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.

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.

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 develop a semi-analytical model for computing planetary system formation with the aim of explaining the observed diversity of planetary systems architectures and relate this primordial diversity with the initial properties of the disc where they were born. We adopt different initial conditions based on recent results in protoplanetary discs observations, to generate a variety of planetary systems and analyze them statistically. We explore the relevance of the mass and size of the disc, its metallicity, the mass of the central star and the time-scale of gaseous disc dissipation, in defining the architecture of the planetary system. We also test different values of some key parameters of our model, to find out which factors best reproduce the diverse sample of observed planetary systems. According to this, we predict which systems are the most common in the solar neighbourhood. Our results show that planetary systems with only terrestrial planets are the most common, being the only planetary systems formed when considering low metallicity discs and which also represent the best environment for the developing of rocky, potentially habitable planets. We also found that planetary systems like our own are not rare in the solar neighbourhood, being its formation favoured in massive discs where there is not a large accumulation of solids in the inner region of the disc. Regarding the planetary systems that harbor hot and warm Jupiter planets, we found that this systems are born in very massive, metal-rich discs. Also a fast migration rate is required in order to form these systems. According to our results, most of the hot and warm Jupiter systems are composed by only one giant planet, which is also a tendency of the current observational data.

Phoebe is one of the irregular satellites of Saturn; the images taken by Cassini-Huygens spacecraft allowed us to analyze its surface and the craters on it. We study the craters on Phoebe produced by Centaur objects from the Scattered Disk (SD) and plutinos escaped from the 3:2 mean motion resonance with Neptune and compare our results with the observations by Cassini. We use previous simulations on trans-Neptunian Objects and a method that allows us to obtain the number of craters and the cratering rate on Phoebe. We obtain the number of craters and the greatest crater on Phoebe produced by Centaurs in the present configuration of the Solar System. Moreover, we obtain a present normalized rate of encounters of Centaurs with Saturn of $\dot F = 7.1 \times 10^{-11}$ per year, from which we can infer the current cratering rate on Phoebe for each crater diameter. Our study and the comparison with the observations suggest that the main crater features on Phoebe are unlikely to have been produced in the present configuration of the Solar System and that they must have been acquired when the SD were depleted in the early Solar System. If this is what happened and the craters were produced when Phoebe was a satellite of Saturn, then it had to be captured, very early in the evolution of the Solar System.

In the last few years, the so-called "Nice model" has got a significant importance in the study of the formation and evolution of the solar system. According to this model, the initial orbital configuration of the giant planets was much more compact than the one we observe today. We study the formation of the giant planets in connection with some parameters that describe the protoplanetary disk. The aim of this study is to establish the conditions that favor their simultaneous formation in line with the initial configuration proposed by the Nice model. We focus in the conditions that lead to the simultaneous formation of two massive cores, corresponding to Jupiter and Saturn, able to achieve the cross-over mass (where the mass of the envelope of the giant planet equals the mass of the core, and gaseous runway starts) while Uranus and Neptune have to be able to grow to their current masses. We compute the in situ planetary formation, employing the numerical code introduced in our previous work, for different density profiles of the protoplanetary disk. Planetesimal migration is taken into account and planetesimals are considered to follow a size distribution between $r_p^{min}$ (free parameter) and $r_p^{max}= 100$ km. The core's growth is computed according to the oligarchic growth regime. The simultaneous formation of the giant planets was successfully completed for several initial conditions of the disk. We find that for protoplanetary disks characterized by a power law ($\Sigma \propto r^{-p}$), smooth surface density profiles ($p \leq 1.5$) favor the simultaneous formation. However, for steep slopes ($p\sim 2$, as previously proposed by other authors) the simultaneous formation of the solar system giant planets is unlikely ...

In order to explain the main characteristics of the observed population of extrasolar planets and the giant planets in the Solar System, we need to get a clear understanding of which are the initial conditions that allowed their formation. To this end we develop a semi-analytical model for computing planetary systems formation based on the core instability model for the gas accretion of the embryos and the oligarchic growth regime for the accretion of the solid cores. With this model we explore not only different initial discs profiles motivated by similarity solutions for viscous accretion discs, but we also consider different initial conditions to generate a variety of planetary systems assuming a large range of discs masses and sizes according to the last results in protoplanetary discs observations. We form a large population of planetary systems in order to explore the effects in the formation of assuming different discs and also the effects of type I and II regimes of planetary migration, which were found to play fundamental role in reproducing the distribution of observed exoplanets. Our results show that the observed population of exoplanets and the giant planets in the Solar System are well represented when considering a surface density profile with a power law in the inner part characterized by an exponent of -1, which represents a softer profile when compared with the case most similar to the MMSN model case.

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.

The core accretion mechanism is presently the most widely accepted cause of the formation of giant planets. For simplicity, most models presently assume that the growth of planetary embryos occurs in isolation. We explore how the simultaneous growth of two embryos at the present locations of Jupiter and Saturn affects the outcome of planetary formation. We model planet formation on the basis of the core accretion scenario and include several key physical ingredients. We consider a protoplanetary gas disk that exponentially decays with time. For planetesimals, we allow for a distribution of sizes from 100~m to 100~km with most of the mass in the smaller objects. We include planetesimal migration as well as different profiles for the surface density $\Sigma$ of the disk. The core growth is computed in the framework of the oligarchic growth regime and includes the viscous enhancement of the planetesimal capture cross-section. Planet migration is ignored. By comparing calculations assuming formation of embryos in isolation to calculations with simultaneous embryo growth, we find that the growth of one embryo generally significantly affects the other. This occurs in spite of the feeding zones of each planet never overlapping. The results may be classified as a function of the gas surface density profile $\Sigma$: if $\Sigma \propto r^{-3/2}$ and the protoplanetary disk is rather massive, Jupiter's formation inhibits the growth of Saturn. If $\Sigma \propto r^{-1}$ isolated and simultaneous formation lead to very similar outcomes; in the the case of $\Sigma \propto r^{-1/2}$ Saturn grows faster and induces a density wave that later acclerates the formation of Jupiter. Our results indicate that the simultaneous growth of several embryos impacts the final outcome and should be taken into account by planet formation models.

It was shown in previous works the existence of weakly chaotic orbits in the plutino population that diffuse very slowly. These orbits correspond to long-term plutino escapers and then represent the plutinos that are escaping from the resonance at present. In this paper we perform numerical simulations in order to explore the dynamical evolution of plutinos recently escaped from the resonance. The numerical simulations were divided in two parts. In the first one we evolved 20,000 test particles in the resonance in order to detect and select the long-term escapers. In the second one, we numerically integrate the selected escaped plutinos in order to study their dynamical post escaped behavior. Our main results include the characterization of the routes of escape of plutinos and their evolution in the Centaur zone. We obtained a present rate of escape of plutinos between 1 and 10 every 10 years. The escaped plutinos have a mean lifetime in the Centaur zone of 108 Myr and their contribution to the Centaur population would be a fraction of less than 6 % of the total Centaur population. In this way, escaped plutinos would be a secondary source of Centaurs.

We develop a simple model of planetary formation, focusing our attention on those planets with masses less than 10 Earth masses and studying particularly the primordial spin parameters of planets resulting from the accretion of planetesimals and produced by the collisions between the embryos. As initial conditions, we adopt the oligarchic growth regime of protoplanets in a disc where several embryos are allowed to form. We take different initial planetary system parameters and for each initial condition, we consider an evolution of 20 millon of years of the system. We perform simulations for 1000 different discs, and from their results we derive the statistical properties of the assembled planets. We have taken special attention to the planetary obliquities and rotation periods, such as the information obtained from the mass and semi major axis diagram, which reflects the process of planetary formation. The distribution of obliquities was found to be isotropic, which means that planets can rotate in direct or indirect sense, regardless of their mass. Our results regarding the primordial rotation periods show that they are dependent on the region where the embryo was formed and evolved. According to our results, most of the planets have rotation periods between 10 and 10000 hours and there are also a large population of planets similar to terrestrial planets in the Solar System.

Giant planet formation process is still not completely understood. The current most accepted paradigm, the core instability model, explains several observed properties of the solar system's giant planets but, to date, has faced difficulties to account for a formation time shorter than the observational estimates of protoplanetary disks' lifetimes, especially for the cases of Uranus and Neptune. In the context of this model, and considering a recently proposed primordial solar system orbital structure, we performed numerical calculations of giant planet formation. Our results show that if accreted planetesimals follow a size distribution in which most of the mass lies in 30-100 meter sized bodies, Jupiter, Saturn, Uranus and Neptune may have formed according to the nucleated instability scenario. The formation of each planet occurs within the time constraints and they end up with core masses in good agreement with present estimations.

The equation of state calculated by Saumon and collaborators has been adopted in most core-accretion simulations of giant-planet formation performed to date. Since some minor errors have been found in their original paper, we present revised simulations of giant-planet formation that considers a corrected equation of state. We employ the same code as Fortier and collaborators in repeating our previous simulations of the formation of Jupiter. Although the general conclusions of Fortier and collaborators remain valid, we obtain significantly lower core masses and shorter formation times in all cases considered. The minor errors in the previously published equation of state have been shown to affect directly the adiabatic gradient and the specific heat, causing an overestimation of both the core masses and formation times.

We develop a simple model for computing planetary formation based on the core instability model for the gas accretion and the oligarchic growth regime for the accretion of the solid core. In this model several planets can form simultaneously in the disc, a fact that has important implications specially for the changes in the dynamic of the planetesimals and the growth of the cores since we consider the collision between them as a source of potential growth. The type I and II migration of the embryos and the migration of the planetesimals due to the interaction with the disc of gas are also taken into account. With this model we consider different initial conditions to generate a variety of planetary systems and analyse them statistically. We explore the effects of using different type I migration rates on the final number of planets formed per planetary system such as on the distribution of masses and semimajor axis of extrasolar planets, where we also analyse the implications of considering different gas accretion rates. A particularly interesting result is the generation of a larger population of habitable planets when the gas accretion rate and type I migration are slower.

It is widely accepted that the large obliquity of Uranus is the result of a great tangential collision (GC) with an Earth size proto-planet at the end of the accretion. We attempt to constraint the GC scenario as the cause of Uranus' obliquity as well as on the mechanisms able to give origin to the Uranian irregulars. Different capture mechanisms for irregulars operate at different stages on the giant planets formation process. The mechanisms able to capture the uranian irregulars before and after the GC are analysed. Assuming that they were captured before the GC, we calculate the orbital transfer of the nine irregulars by the impulse imparted by the GC. If their orbital transfer results dynamically implausible, they should have originated after the GC. We investigate and discuss the dissipative mechanisms able to operate later. In particular Prospero could not exist at the time of the GC. Different capture mechanisms for Prospero after the GC are investigated. Gas drag by Uranus'envelope and pull-down capture are not plausible mechanisms. Capture of Prospero through a collisionless interaction seems to be difficult. The GC itself provides a mechanism of permanent capture. However, the capture of Prospero by the GC is a low probable event. Catastrophic collisions could be a possible mechanism for the birth of Prospero and the other irregulars after the GC. Orbital and physical clusterings should then be expected. Either Prospero had to originate after the GC or the GC did not occur. In the former case, the mechanism for the origin of Prospero after the GC remains an open question. In the latter case, another theory to account for Uranus' obliquity and the formation of the Uranian regular satellites on the equatorial plane of the planet would be needed.

Aims. In the context of the core instability model, we present calculations of in situ giant planet formation. The oligarchic growth regime of solid protoplanets is the model adopted for the growth of the core. Methods. The full differential equations of giant planet formation were numerically solved with an adaptation of a Henyey-type code. The planetesimals accretion rate was coupled in a self-consistent way to the envelope's evolution. Results. We performed several simulations for the formation of a Jupiter-like object by assuming various surface densities for the protoplanetary disc and two different sizes for the accreted planetesimals. We find that the atmospheric gas drag gives rise to a major enhancement on the effective capture radius of the protoplanet, thus leading to an average timescale reduction of 30% -- 55% and ultimately to an increase by a factor of 2 of the final mass of solids accreted as compared to the situation in which drag effects are neglected. With regard to the size of accreted planetesimals, we find that for a swarm of planetesimals having a radius of 10 km, the formation time is a factor 2 to 3 shorter than that of planetesimals of 100 km, the factor depending on the surface density of the nebula. Moreover, planetesimal size does not seem to have a significant impact on the final mass of the core.

We have withdrawn this paper because we have presented the basic results in Brunini & Cionco (2005), Icarus 177, 264

The sample of known exoplanets is strongly biased to masses larger than the ones of the giant gaseous planets of the solar system. Recently, the discovery of two extrasolar planets of considerably lower masses around the nearby stars GJ 436 and $\rho$ Cancri was reported. They are like our outermost icy giants, Uranus and Neptune, but in contrast, these new planets are orbiting at only some hundredth of the Earth-Sun distance from their host stars, raising several new questions about their origin and constitution. Here we report numerical simulations of planetary accretion that show, for the first time through N-body integrations that the formation of compact systems of Neptune-like planets close to the hosts stars could be a common by-product of planetary formation. We found a regime of planetary accretion, in which orbital migration accumulates protoplanets in a narrow region around the inner edge of the nebula, where they collide each other giving rise to Neptune-like planets. Our results suggest that, if a protoplanetary solar environment is common in the galaxy, the discovery of a vast population of this sort of 'hot cores' should be expected in the near future.