In about 6 Giga years our Sun will evolve into a red giant and finally end its life as a white dwarf. This stellar metamorphosis will occur to virtually all known host stars of exo-planetary systems and is therefore crucial for their final fate. It is clear that the innermost planets will be engulfed and evaporated during the giant phase and that planets located farther out will survive. However, the destiny of planets in-between, at ~ 1 - 10 au, has not yet been investigated with a multi-planet tidal treatment. We here combine for the first time multi-planet interactions, stellar evolution, and tidal effects in an N-body code to study the evolution of a Neptune-Jupiter planetary system. We report that the fate of the Neptune-mass planet, located closer to the star than the Jupiter-mass planet, can be very different from the fate of a single Neptune. The simultaneous effects of gravitational interactions, mass loss and tides can drive the planetary system towards mean motion resonances. Crossing these resonances affects particularly the eccentricity of the Neptune and thereby also its fate, which can be engulfment, collision with the Jupiter-mass planet, ejection from the system, or survival at a larger separation.

The diversity in mass and composition of planetary atmospheres stems from the different building blocks present in protoplanetary discs and from the different physical and chemical processes that these experience during the planetary assembly and evolution. This review aims to summarise, in a nutshell, the key concepts and processes operating during planet formation, with a focus on the delivery of volatiles to the inner regions of the planetary system.

Recent observations of protoplanetary disks have revealed ring-like structures that can be associated to pressure maxima. Pressure maxima are known to be dust collectors and planet migration traps. Most of planet formation works are based either on the pebble accretion model or on the planetesimal accretion model. However, recent studies proposed the possible formation of Jupiter by the hybrid accretion of pebbles and planetesimals. We aim to study the full process of planet formation consisting of dust evolution, planetesimal formation and planet growth at a pressure maximum in a protoplanetary disk. We compute, through numerical simulations, the gas and dust evolution, including dust growth, fragmentation, radial drift and particle accumulation at a pressure bump. We also consider the formation of planetesimals by streaming instability and the formation of a moon-size embryo that grows into a giant planet by the hybrid accretion of pebbles and planetesimals. We find that pressure maxima in protoplanetary disks are efficient collectors of dust drifting inwards. The condition of planetesimal formation by streaming instability is fulfilled due to the large amount of dust accumulated at the pressure bump. Then, a massive core is quickly formed (in $\sim 10^4$ yr) by the accretion of pebbles. After the pebble isolation mass is reached, the growth of the core slowly continues by the accretion of planetesimals. The energy released by planetesimal accretion delays the onset of runaway gas accretion, allowing a gas giant to form after $\sim$1 Myr of disk evolution. The pressure maximum also acts as a migration trap. Pressure maxima in protoplanetary disks are preferential locations for dust traps, planetesimal formation by streaming instability and planet migration traps. All these conditions allow the fast formation of a giant planet by the hybrid accretion of pebbles and planetesimals.

As planets grow the exchange of angular momentum with the gaseous component of the protoplanetary disc produces a net torque resulting in a variation of the semi-major axis of the planet. For low-mass planets not able to open a gap in the gaseous disc this regime is known as type I migration. Pioneer works studied this mechanism in isothermal discs finding fast inward type I migration rates that were unable to reproduce the observed properties of extrasolar planets. In the last years, several improvements have been made in order to extend the study of type I migration rates to non-isothermal discs. Moreover, it was recently shown that if the planet's luminosity due to solid accretion is taken into account, inward migration could be slowed down and even reversed. In this work, we study the planet formation process incorporating, and comparing, updated type I migration rates for non-isothermal discs and the role of planet's luminosity over such rates. We find that the latter can have important effects on planetary evolution, producing a significant outward migration for the growing planets.

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.

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.

Some recent studies exposed rather strong statistical evidence of in-vacuo-dispersion-like spectral lags for gamma-ray bursts (GRBs), a linear correlation between time of observation and energy of GRB particles. Those results focused on testing in-vacuo dispersion for the most energetic GRB particles, and in particular only included photons with energy at emission greater than 40 GeV. We here extend the window of the statistical analysis down to 5 GeV and find results that are consistent with what had been previously noticed at higher energies.

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 discovery of giant planets in wide orbits represents a major challenge for planet formation theory. In the standard core accretion paradigm planets are expected to form at radial distances $\lesssim 20$ au in order to form massive cores (with masses $\gtrsim 10~\textrm{M}_{\oplus}$) able to trigger the gaseous runaway growth before the dissipation of the disc. This has encouraged authors to find modifications of the standard scenario as well as alternative theories like the formation of planets by gravitational instabilities in the disc to explain the existence of giant planets in wide orbits. However, there is not yet consensus on how these systems are formed. In this letter, we present a new natural mechanism for the formation of giant planets in wide orbits within the core accretion paradigm. If photoevaporation is considered, after a few Myr of viscous evolution a gap in the gaseous disc is opened. We found that, under particular circumstances planet migration becomes synchronised with the evolution of the gap, which results in an efficient outward planet migration. This mechanism is found to allow the formation of giant planets with masses $M_p\lesssim 1 M_{\rm Jup}$ in wide stable orbits as large as $\sim$130 au from the central star.

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.

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.