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.

The goal of this work is to investigate Jupiter's growth focusing on the amount of heavy elements accreted by the planet, and its comparison with recent structure models. Our model assumes an initial core growth dominated by pebble accretion, and a second growth phase that is characterized by a moderate accretion of both planetesimals and gas. The third phase is dominated by runaway gas accretion during which the planet becomes detached from the disk. The second and third phases are computed in detail, considering two different prescriptions for the planetesimal accretion and fits from hydrodynamical studies to compute the gas accretion in the detached phase. In order for Jupiter to consist of $\sim$20-40 $M_{\oplus}$ of heavy elements as suggested by structure models, we find that Jupiter's formation location is preferably at an orbital distance of $1\lesssim a \lesssim 10$ au once the accretion of planetesimals dominates. We find that Jupiter could accrete between $\sim$1 and $\sim$15 $M_{\oplus}$ of heavy elements during runaway gas accretion, depending on the assumed initial surface density of planetesimals and the prescription used to estimate the heavy-element accretion during the final stage of the planetary formation. This would yield an envelope metallicity of $\sim$0.5 to $\sim$3 times solar. By computing the solid accretion during the detached phase, we infer a planetary mass-metallicity ($M_P$-$M_Z$) relation of $M_Z \sim M_P^{2/5}$ when a gap in the planetesimal disk is created, and of $M_Z \sim M_P^{1/6}$ without a planetesimal gap. Our hybrid pebble-planetesimal model can account for Jupiter's bulk and atmospheric enrichment. The high bulk metallicity inferred for many giant exoplanets is difficult to explain from standard formation models. This might suggest a migration history for such highly enriched giant exoplanets and/or giant impacts after the disk's dispersal.

Computing the mass of planetary envelopes and the critical mass beyond which planets accrete gas in a runaway fashion is important when studying planet formation, in particular for planets up to the Neptune mass range. This computation requires in principle solving a set of differential equations, the internal structure equations, for some boundary conditions (pressure, temperature in the protoplanetary disk where a planet forms, core mass and accretion rate of solids by the planet). Solving these equations in turn proves being time consuming and sometimes numerically unstable. We developed a method to approximate the result of integrating the internal structure equations for a variety of boundary conditions. We compute a set of planet internal structures for a very large number (millions) of boundary conditions, considering two opacities,(ISM and reduced). This database is then used to train Deep Neural Networks in order to predict the critical core mass as well as the mass of planetary envelopes as a function of the boundary conditions. We show that our neural networks provide a very good approximation (at the level of percents) of the result obtained by solving interior structure equations, but with a much smaller required computer time. The difference with the real solution is much smaller than the one obtained using some analytical formulas available in the literature which at best only provide the correct order of magnitude. We compare the results of the DNN with other popular machine learning methods (Random Forest, Gradient Boost, Support Vector Regression) and show that the DNN outperforms these methods by a factor of at least two. We show that some analytical formulas that can be found in various papers can severely overestimate the mass of planets, therefore predicting the formation of planets in the Jupiter-mass regime instead of the Neptune-mass regime.

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.

Constraining the planetary composition is essential for exoplanetary characterization. In this paper, we use a statistical analysis to determine the characteristic maximum (threshold) radii for various compositions for exoplanets with masses up to 25 Earth masses (M$_\oplus$). We confirm that most planets with radii larger than 1.6 Earth radius (R$_\oplus$) are not rocky, and must consist of lighter elements, as found by previous studies. We find that planets with radii above 2.6 R$_\oplus$ cannot be pure-water worlds, and must contain significant amounts of hydrogen and helium (H-He). We find that planets with radii larger than about 3 R$_\oplus$, 3.6 R$_\oplus$, and 4.3 R$_\oplus$ are expected to consist of 2%, 5% and 10% of H-He, respectively. We investigate the sensitivity of the results to the assumed internal structure, the planetary temperature and albedo, and the accuracy of the mass and radius determination. We show that the envelope's metallicity, the percentage of H-He and the distribution of the elements play a significant role in the determination of the threshold radius. Finally, we conclude that despite the degenerate nature of the problem, it is possible to put limits on the possible range of compositions for planets with well-measured mass and radius.

The formation history of giant planets determines their primordial structure and consequent evolution. We simulate various formation paths of Jupiter to determine its primordial entropy, and find that a common outcome is for proto-Jupiter to have non-convective regions in its interior. We use planet formation models to calculate how the entropy and post-formation luminosity depend on model properties such as the solid accretion rate and opacity, and show that the gas accretion rate and its time evolution play a key role in determining the entropy profile. The predicted luminosity of Jupiter shortly after formation varies by a factor of $2$--$3$ for different choices of model parameters. We find that entropy gradients inside Jupiter persist for $\sim 10\ {\rm Myr}$ after formation. We suggest that these gradients should be considered together with heavy-element composition gradients when modeling Jupiter's evolution and internal structure.

Mini-Neptunes seem to be common planets. In this work we investigate the possible formation histories and predicted occurrence rates of mini-Neptunes assuming the planets form beyond the iceline. We consider pebble and planetesimal accretion accounting for envelope enrichment and two different opacity conditions. We find that the formation of mini-Neptunes is a relatively frequent output when envelope enrichment by volatiles is included, and that there is a "sweet spot" for mini-Neptune formation with a relatively low solid accretion rate of ~10^{-6} Earth masses per year. This rate is typical for low/intermediate-mass protoplanetary disks and/or disks with low metallicities. With pebble accretion, envelope enrichment and high opacity favor the formation of mini-Neptunes, with more efficient formation at large semi-major axes (~30 AU) and low disk viscosity. For planetesimal accretion, such planets can form also without enrichment, with the opacity being a key aspect in the growth history and favorable formation location. Finally, we show that the formation of Neptune-like planets remains a challenge for planet formation theories.

Using a generalized Bayesian inference method, we aim to explore the possible interior structures of six selected exoplanets for which planetary mass and radius measurements are available in addition to stellar host abundances: HD~219134b, Kepler-10b, Kepler-93b, CoRoT-7b, 55~Cnc~e, and HD~97658b. We aim to investigate the importance of stellar abundance proxies for the planetary bulk composition (namely Fe/Si and Mg/Si) on prediction of planetary interiors. We performed a full probabilistic Bayesian inference analysis to formally account for observational and model uncertainties while obtaining confidence regions of structural and compositional parameters of core, mantle, ice layer, ocean, and atmosphere. We determined how sensitive our parameter predictions depend on (1) different estimates of bulk abundance constraints and (2) different correlations of bulk abundances between planet and host star. [...] Although the possible ranges of interior structures are large, structural parameters and their correlations are constrained by the sparse data. The probability for the tested exoplanets to be Earth-like is generally very low. Furthermore, we conclude that different estimates of planet bulk abundance constraints mainly affect mantle composition and core size.

We aim to present a generalized Bayesian inference method for constraining interiors of super Earths and sub-Neptunes. Our methodology succeeds in quantifying the degeneracy and correlation of structural parameters for high dimensional parameter spaces. Specifically, we identify what constraints can be placed on composition and thickness of core, mantle, ice, ocean, and atmospheric layers given observations of mass, radius, and bulk refractory abundance constraints (Fe, Mg, Si) from observations of the host star's photospheric composition. We employed a full probabilistic Bayesian inference analysis that formally accounts for observational and model uncertainties. Using a Markov chain Monte Carlo technique, we computed joint and marginal posterior probability distributions for all structural parameters of interest. We included state-of-the-art structural models based on self-consistent thermodynamics of core, mantle, high-pressure ice, and liquid water. Furthermore, we tested and compared two different atmospheric models that are tailored for modeling thick and thin atmospheres, respectively. First, we validate our method against Neptune. Second, we apply it to synthetic exoplanets of fixed mass and determine the effect on interior structure and composition when (1) radius, (2) atmospheric model, (3) data uncertainties, (4) semi-major axes, (5) atmospheric composition (i.e., a priori assumption of enriched envelopes versus pure H/He envelopes), and (6) prior distributions are varied. Our main conclusions are: [...]

We compute, for the first time, self-consistent models of planet growth including the effect of envelope enrichment. The change of envelope metallicity is assumed to be the result of planetesimal disruption or icy pebble sublimation. We solve internal structure equations taking into account global energy conservation for the envelope to compute in-situ planetary growth. We consider different opacities and equations of state suited for a wide range of metallicities. We find that envelope enrichment speeds up the formation of gas giants. It also explains naturally the formation of low and intermediate mass objects with large fractions of H-He (~ 20 - 30 % in mass). High opacity models explain well the metallicity of the giant planets of the solar system, whereas low opacity models are suited for forming small mass objects with thick H-He envelopes and gas giants with sub-solar envelope metallicities. We find good agreement between our models and the estimated water abundance for WASP-43b. For HD 189733b, HD 209458b and WASP-12b we predict fractions of water larger than what is estimated from observations, by at least a factor ~ 2. Envelope enrichment by icy planetesimals is the natural scenario to explain the formation of a large variety of objects, ranging from mini-Neptunes, to gas giants. We predict that the total and envelope metallicity decrease with planetary mass.

Context. Within the core accretion scenario of planetary formation, most simulations performed so far always assume the accreting envelope to have a solar composition. From the study of meteorite showers on Earth and numerical simulations, we know that planetesimals must undergo thermal ablation and disruption when crossing a protoplanetary envelope. Once the protoplanet has acquired an atmosphere, the primordial envelope gets enriched in volatiles and silicates from the planetesimals. This change of envelope composition during the formation can have a significant effect in the final atmospheric composition and on the formation timescale of giant planets. Aims. To investigate the physical implications of considering the envelope enrichment of protoplanets due to the disruption of icy planetesimals during their way to the core. Particular focus is placed on the effect on the critical core mass for envelopes where condensation of water can occur. Methods. Internal structure models are numerically solved with the implementation of updated opacities for all ranges of metallicities and the software CEA to compute the equation of state. CEA computes the chemical equilibrium for an arbitrary mixture of gases and allows the condensation of some species, including water. This means that the latent heat of phase transitions is consistently incorporated in the total energy budget. Results. The critical core mass is found to decrease significantly when an enriched envelope composition is considered in the internal structure equations. A particular strong reduction of the critical core mass is obtained for planets whose envelope metallicity is larger than Z=0.45 when the outer boundary conditions are suitable for condensation of water to occur in the top layers of the atmosphere. We show that this effect is qualitatively preserved when the atmosphere is out of chemical equilibrium.

At present, approximately 1500 asteroids are known to evolve inside or sticked to the exterior 1:2 resonance with Mars at a = 2.418 AU, being (142) Polana the largest member of this group. The effect of the forced secular modes superposed to the resonance gives rise to a complex dynamical evolution. Chaotic diffusion, collisions, close encounters with massive asteroids and mainly orbital migration due to the Yarkovsky effect generate continuous captures to and losses from the resonance, with a fraction of asteroids remaining captured over long time scales and generating a concentration in the semimajor axis distribution that exceeds by 20% the population of background asteroids. The Yarkovsky effect induces different dynamics according to the asteroid size, producing an excess of small asteroids inside the resonance. The evolution in the resonance generates a signature on the orbits, mainly in eccentricity, that depends on the time the asteroid remains captured inside the resonance and on the magnitude of the Yarkovsky effect. The greater the asteroids, the larger the time they remain captured in the resonance, allowing greater diffusion in eccentricity and inclination. The resonance generates a discontinuity and mixing in the space of proper elements producing misidentification of dynamical family members, mainly for Vesta and Nysa-Polana families. The half-life of resonant asteroids large enough for not being affected by the Yarkovsky effect is about 1 Gyr. From the point of view of taxonomic classes, the resonant population does not differ from the background population and the excess of small asteroids is confirmed.

Using Gauss' averaged equations, we compute the secular relativistic effects generated by the Sun on the argument of the perihelion and the mean anomaly of an orbit. Then we test different alternative simpler models that have been proposed to reproduce the secular relativistic effects in the orbital elements. Generally, models introduce artificial perturbations that are velocity-independent but that depend on the heliocentric distance. If these perturbations are set as an impulse in a constant timestep integrator, when the particle approaches perihelion the generated impulse could be very strong and badly sampled, originating a spurious orbital evolution. In order to overcome this setback, we propose two new models based on a constant, distance-independent, perturbation. With these models we obtain the correct secular drift in the argument of perihelion and the expected secular orbital evolution is reproduced. We also discuss with some detail the secular effect generated on the mean anomaly by different models.