Aims. Comparing theoretical models of exoplanet formation with observations allows one to make key step forward towards an understanding of planetary systems. It however requires a model able to (i) predict all the necessary observable quantities (not only masses and orbits, but also radii, luminosities, magnitudes, or evaporation rates) and (ii) address the large range in relevant planetary masses (from Mars mass to super-Jupiters) and distances (from stellar-grazing to wide orbits). Methods. We have developed a combined global planetary formation and evolution model that solves as directly as possible the underlying differential equations for the structure and evolution of the gas disc, the dynamical state of the planetesimals, the internal structure of the planets yielding their gas accretion rate and internal structure, the accretion rate of planetesimals, disc-driven orbital migration, and the gravitational interaction of concurrently forming planets via a full N-body calculation. Importantly, the model now also follows the long-term evolution of the planets on Gigayear timescales after formation including the effects of cooling and contraction, atmospheric escape, bloating, and stellar tides. Results. To test the model, we compared it with classical scenarios of Solar System formation. For the terrestrial planets, we find that we obtain a giant impact phase provided enough embryos (~100) are initially emplaced in the disc. For the giant planets, we find that Jupiter-mass planets must accrete their core shortly before the dispersal of the gas disc to prevent strong inward migration that would bring them to the inner edge of the disc. Conclusions. he model can form planetary systems with a wide range of properties. We find that giant-planet bearing systems are more diverse than the ones with only terrestrial planets.

Recent observational findings have suggested a positive correlation between the occurrence rates of inner super-Earths and outer giant planets. These results raise the question if the trend can be reproduced and explained by planet formation theory. Here, we investigate the properties of inner super-Earths and outer giant planets that form according to the core accretion picture. The aim of this research was to study mutual relations between these planet species in synthetic planet populations and to compare them to the observed exoplanet population. We invoke the Generation 3 Bern model of planet formation and evolution to simulate 1000 multi-planet systems. We then confront these synthetic systems with the observed sample, taking into account the detection bias that distorts the observed demographics. The formation of warm super-Earths and cold Jupiters in the same system is enhanced compared to the individual appearances, albeit weaker than proposed by observations. We attribute the discrepancy to dynamically active giant planets, which frequently disrupt inner super-Earth systems in the synthetic population. This result could be interpreted as an overestimation of the migration rates predicted by giant planet migration theory. We further find differences in the volatile content of planets in different system architectures and predict that high-density super-Earths are more likely to host an outer giant. This correlation can be tested observationally. Finally, super-Earth occurrence shows a negative correlation with metallicity in giant-hosting systems. The reason for this is that giant planets around high-metallicity stars are frequently dynamically active, which poses a threat to inner planetary systems.

We want to understand the global observable consequences of different physical processes and the initial system properties on the demographics of the planetary population. We select distributions of initial conditions that are representative of known protoplanetary discs and we use the Bern model to synthesise five planetary populations around 1 Sun-mass stars. Each population has a different initial number of Moon-mass embryos per disc to test for convergence: 1, 10, 20, 50, and 100. The last one is our nominal planetary population. The properties of giant planets do not change much starting with 10 embryos per system. For terrestrial planets inside 1 au, only the 100-embryos population is able to attain the giant-impact stage. In that population, each system contains 8 planets larger than 1 $M_\oplus$ on average. The fraction of systems with giants planets at all orbital distances is 18%, but only 1.6% are at > 10 au. Systems with giants contain on average 1.6 such planets. Inside of 1 au, the planet type with the highest occurrence rate and multiplicity are super-Earths (2.4 and 3.7), followed by terrestrial planets (1.6 and 2.8). Their frequency peaks at a stellar [Fe/H] of 0.0 and -0.2 respectively. The frequency of more massive planet types (Neptunian, giants) increases monotonically with [Fe/H]. Low-mass planets in the habitable zone are found around 44% of the stars with a multiplicity of 1.3, with a mean [Fe/H] of -0.11. This is one of the most comprehensive simulations of planetary system formation and evolution to date. For observations, the syntheses provide a large data set to search for comparison synthetic systems that show how they formed. For theory, they provide the framework to observationally test the consequences of specific physical processes. This is a important ingredient towards the development of a standard model of planetary formation and evolution.

In the core accretion scenario, a massive core forms first and then accretes an envelope. When discussing how this core forms some divergences appear. First scenarios of planet formation predict the accretion of km-sized bodies, called planetesimals, while more recent works suggest growth by accretion of pebbles, which are cm-sized objects. These two accretion models are often discussed separately and we aim here at comparing the outcomes of the two models with identical initial conditions. We use two distinct codes: one computing planetesimal accretion, the other pebble accretion. Using a population synthesis approach, we compare planet simulations and study the impact of the two solid accretion models, focussing on the formation of single planets. We find that the planetesimal model predicts the formation of more giant planets, while the pebble accretion model forms more super-Earth mass planets. This is due to the pebble isolation mass concept, which prevents planets formed by pebble accretion to accrete gas efficiently before reaching Miso. This translates into a population of planets that are not heavy enough to accrete a consequent envelope but that are in a mass range where type I migration is very efficient. We also find higher gas mass fractions for a given core mass for the pebble model compared to the planetesimal one caused by luminosity differences. This also implies planets with lower densities which could be confirmed observationally. Focusing on giant planets, we conclude that the sensitivity of their formation differs: for the pebble accretion model, the time at which the embryos are formed, as well as the period over which solids are accreted strongly impact the results, while for the planetesimal model it depends on the planetesimal size and on the splitting in the amount of solids available to form planetesimals.

Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought.

Context. The composition of a protoplanetary disk at a given location does not only depend on temperature and pressure but also on the time dependent transport of matter, such as radial drift of solid bodies, which could release water and other volatile species before disintegration or accretion onto a larger body with potentially considerable implications for the composition of planets. Aims. We perform a parameter study focused on the water depletion of different sized bodies able to cross the water snowline by gas induced radial drift. Methods. Either the analytical Hertz-Knudsen-Langmuir sublimation formula assuming equilibrium temperature within the body or a more involved, numerical model for the internal thermal evolution is coupled with an alpha-disk model. Different properties of the disk and the embedded body are explored. Results. Bodies with radii up to 100 meter drift faster towards the central star than the water snowline, hence, cross it. The region that can be reached before complete disintegration - and is therefore polluted with H2O ice - extends to 10 percent closer to the star than the snowline location. The extent of this polluted region could be multiple times larger in the presence of a dust mantle, which is, however, unlikely to form due to frequent collisions with smaller-than centimeter sized objects. Conclusions. Given a significant abundance of meter sized boulders in protoplanetary disks, the transport of water by radial drift of these bodies towards regions closer to the star than the snowline is not negligible and this flux of volatiles can be estimated for a given distribution of solid body sizes and compositions. A simple expression for surface sublimation is applicable for a homogeneous body consisting of only dust and water ice without a dust mantle.

In contrast to the water-poor inner solar system planets, stochasticity during planetary formation and order of magnitude deviations in exoplanet volatile contents suggest that rocky worlds engulfed in thick volatile ice layers are the dominant family of terrestrial analogues among the extrasolar planet population. However, the distribution of compositionally Earth-like planets remains insufficiently constrained, and it is not clear whether the solar system is a statistical outlier or can be explained by more general planetary formation processes. Here we employ numerical models of planet formation, evolution, and interior structure, to show that a planet's bulk water fraction and radius are anti-correlated with initial $^{26}$Al levels in the planetesimal-based accretion framework. The heat generated by this short-lived radionuclide rapidly dehydrates planetesimals prior to accretion onto larger protoplanets and yields a system-wide correlation of planet bulk abundances, which, for instance, can explain the lack of a clear orbital trend in the water budgets of the TRAPPIST-1 planets. Qualitatively, our models suggest two main scenarios of planetary systems' formation: high-$^{26}$Al systems, like our solar system, form small, water-depleted planets, whereas those devoid of $^{26}$Al predominantly form ocean worlds, where the mean planet radii between both scenarios deviate by up to about 10%.

The standard model for giant planet formation is based on the accretion of solids by a growing planetary embryo, followed by rapid gas accretion once the planet exceeds a so-called critical mass. The dominant size of the accreted solids (cm-size particles named pebbles or km to hundred km-size bodies named planetesimals) is, however, unknown. Recently, high-precision measurements of isotopes in meteorites provided evidence for the existence of two reservoirs in the early Solar System. These reservoirs remained separated from ~1 until ~ 3 Myr after the beginning of the Solar System's formation. This separation is interpreted as resulting from Jupiter growing and becoming a barrier for material transport. In this framework, Jupiter reached ~20 Earth masses within ~1 Myr and slowly grew to ~50 Earth masses in the subsequent 2 Myr before reaching its present-day mass. The evidence that Jupiter slowed down its growth after reaching 20 Earth masses for at least 2 Myr is puzzling because a planet of this mass is expected to trigger fast runaway gas accretion. Here, we use theoretical models to describe the conditions allowing for such a slow accretion and show that Jupiter grew in three distinct phases. First, rapid pebble accretion brought the major part of Jupiter's core mass. Second, slow planetesimal accretion provided the energy required to hinder runaway gas accretion during 2 Myr. Third, runaway gas accretion proceeded. Both pebbles and planetesimals therefore have an important role in Jupiter's formation.