The thermal structure of a protoplanetary disc is regulated by the opacity that dust grains provide. However, previous works have often considered simplified prescriptions for the dust opacity in hydrodynamical disc simulations, e.g. by considering only a single particle size. In the present work we perform 2D hydrodynamical simulations of protoplanetary discs where the opacity is self-consistently calculated for the dust population, taking into account the particle size, composition and abundance. We first compare simulations using single grain sizes to two different multi-grain size distributions at different levels of turbulence strengths, parameterized through the $\alpha$-viscosity, and different gas surface densities. Assuming a single dust size leads to inaccurate calculations of the thermal structure of discs, because the grain size dominating the opacity increases with orbital radius. Overall the two grain size distributions, one limited by fragmentation only and the other determined from a more complete fragmentation-coagulation equilibrium, give similar results for the thermal structure. We find that both grain size distributions give less steep opacity gradients that result in less steep aspect ratio gradients, in comparison to discs with only micrometer sized dust. Moreover, in the discs with a grain size distribution, the innermost outward migration region is removed and planets embedded is such discs experience lower migration rates. We also investigate the dependency of the water iceline position on the alpha-viscosity, the initial gas surface density at 1 AU and the dust-to-gas ratio and find $r_{ice} \propto \alpha^{0.61} \Sigma_{g,0}^{0.8} f_{DG}^{0.37}$ independently of the distribution used. The inclusion of the feedback loop between grain growth, opacities and disc thermodynamics allows for more self-consistent simulations of accretion discs and planet formation.

We conduct a pebble-driven planet population synthesis study to investigate the formation of planets around very low-mass stars and brown dwarfs, in the (sub)stellar mass range between $0.01 \ M_{\odot}$ and $0.1 \ M_{\odot}$. Based on the extrapolation of numerical simulations of planetesimal formation by the streaming instability, we obtain the characteristic mass of the planetesimals and the initial masses of the protoplanets (largest bodies from the planetesimal size distributions), in either the early self-gravitating phase or the later non-self-gravitating phase of the protoplanetary disk evolution. We find that the initial protoplanets form with masses that increase with host mass, orbital distance and decrease with disk age. Around late M-dwarfs of $0.1 \ M_{\odot}$, these protoplanets can grow up to Earth-mass planets by pebble accretion. However, around brown dwarfs of $0.01 \ M_{\odot}$, planets do not grow larger than Mars mass when the initial protoplanets are born early in self-gravitating disks, and their growth stalls at around $0.01$ Earth-mass when they are born late in non-self-gravitating disks. Around these low mass stars and brown dwarfs, we find no channel for gas giant planet formation because the solid cores remain too small. When the initial protoplanets form only at the water-ice line, the final planets typically have ${\gtrsim} 15\%$ water mass fraction. Alternatively, when the initial protoplanets form log-uniformly distributed over the entire protoplanetary disk, the final planets are either very water-rich (water mass fraction ${\gtrsim}15\%$) or entirely rocky (water mass fraction ${\lesssim}5\%$).

Recent surveys have revealed that protoplanetary disks typically have dust masses that appear to be insufficient to account for the high occurrence rate of exoplanet systems. We demonstrate that this observed dust depletion is consistent with the radial drift of pebbles. Using a Monte Carlo method we simulate the evolution of a cluster of protoplanetary disks, using a 1D numerical method to viscously evolve each gas disk together with the radial drift of dust particles that have grown to 100 $\mu$m in size. For a 2 Myr old cluster of stars, we find a slightly sub-linear scaling between the gas disk mass and the gas accretion rate ($M_\mathrm{g}\propto\dot{M}^{0.9}$). However, for the dust mass we find that evolved dust disks have a much weaker scaling with the gas accretion rate, with the precise scaling depending on the age at which the cluster is sampled and the intrinsic age spread of the disks in the cluster. Ultimately, we find that the dust mass present in protoplanetary disk is on the order of 10-100 Earth masses in 1-3 Myr old star-forming regions, a factor of 10 to 100 depleted from the original dust budget. As the dust drains from the outer disk, pebbles pile up in the inner disk and locally increase the dust-to-gas ratio by a factor of up to 4 above the initial value. In these high dust-to-gas ratio regions we find conditions that are favourable for planetesimal formation via the streaming instability and subsequent growth by pebble accretion. We also find the following scaling relations with stellar mass within a 1-2 Myr old cluster: a slightly super-linear scaling between the gas accretion rate and stellar mass ($\dot{M}\propto M_\star^{1.4}$), a slightly super-linear scaling between the gas disk mass and the stellar mass ($M_\mathrm{g}\propto M_\star^{1.4}$) and a super-linear relation between the dust disk mass and stellar mass ($M_\mathrm{d}\propto M_\star^{1.4-4.1}$).

We study the dynamics and growth of dust particles in circumstellar disks of different masses that are prone to gravitational instability during the critical first Myr of their evolution. The dust component is made up of two different components: micron-sized dust and grown dust of evolving size. For the dust component, we considered the dust coagulation, fragmentation, momentum exchange with the gas, and dust self-gravity. We found that the micron-sized dust particles grow rapidly in the circumstellar disk, reaching a few cm in size in the inner 100 au of the disk during less than 100 kyr after the disk formation, provided that fragmentation velocity is $30\rm~ms^{-1}$. Due to the accretion of micron dust particles from the surrounding envelope, which serves as a micron dust reservoir, the approximately cm-sized dust particles continue to be present in the disk for more than 900 kyr after the disk formation and maintain a dust-to-gas ratio close to 0.01. We show that a strong correlation exists between the gas and pebble fluxes in the disk. We find that radial surface density distribution of pebbles in the disk shows power-law distribution with an index similar to that of the Minimum-mass solar nebula (MMSN) regardless the disk mass. We also show that the gas surface density in our models agrees well with measurements of dust in protoplanetary disks of AS 209, HD 163296, and DoAr 25 systems. Pebbles are formed during the very early stages of protoplanetary disk evolution. They play a crucial role in the planet formation process. Our disc simulations reveal the early onset ($<10^5$ yr) of an inwards-drifting flux of pebble-sized particles that makes up approximately between one hundredth and one tenth of the gas mass flux, which appears consistent with mm-observations of discs. (Abridged)

Circumplanetary discs (CPDs) form around young gas giants and are thought to be the sites of moon formation as well as an intermediate reservoir of gas that feeds the growth of the gas giant. How the physical properties of such CPDs are affected by the planetary mass and the overall opacity is relatively poorly understood. In order to clarify this, we use the global radiation hydrodynamics code FARGOCA, with a grid structure that allows resolving the planetary gravitational potential sufficiently well for a CPD to form. We then study the gas flows and density/temperature structures that emerge as a function of planet mass, opacity and potential depth. Our results indicate interesting structure formation for Jupiter-mass planets at low opacities, which we subsequently analyse in detail. Using an opacity level that is 100 times lower than that of ISM dust, our Jupiter-mass protoplanet features an envelope that is sufficiently cold for a CPD to form, and a free-fall region separating the CPD and the circumstellar disc emerges. Interestingly, this free-fall region appears to be a result of supersonic erosion of outer envelope material, as opposed to the static structure formation that one would expect at low opacities. Our analysis reveals that the planetary spiral arms seem to pose a significant pressure barrier that needs to be overcome through radiative cooling in order for gas to free-fall onto the CPD. The circulation inside the CPD is near-keplerian and modified by the presence of CPD spiral arms. For high opacities we recover results from the literature, finding an essentially featureless hot envelope. With this work, we demonstrate the first simulation and analysis of a complete detachment process of a protoplanet from its parent disc in a 3D radiation hydrodynamics setting.

We develop a pebble-driven model to study the formation and evolution of planets around stars in the mass range of 0.08 and 1 solar mass. The growth and migration of a large number of individual protoplanetary embryos are simulated in a population synthesis manner. We test two hypotheses for the birth locations of embryos: at the water ice line or log-uniformly distributed over entire protoplanetary disks. Two types of disks with different turbulent viscous parameters alpha of 1e-3 and 1e-4 are investigated, to shed light on the role of outward migration of protoplanets. The forming planets are compared with the observed exoplanets in terms of masses, semimajor axes, metallicities, and water contents. We find that gas giant planets are likely to form when the characteristic disk sizes are larger, the disk accretion rates are higher, the disks are more metal-rich and/or their stellar hosts are more massive. Our model shows that 1) the characteristic mass of super-Earth is set by the pebble isolation mass. Super-Earth masses increase linearly with the mass of its stellar host, corresponding to one Earth mass around a late M-dwarf star and 20 Earth masses around a solar-mass star. 2) The low-mass planets up to 20 Earth masses can form around stars with a wide range of metallicities, while massive gas giant planets are preferred to grow around metal-rich stars. 3) Super-Earth planets that are mainly composed of silicates, with relatively low water fractions can form from protoplanetary embryos at the water ice line in weakly turbulent disks where outward migration is suppressed. However, if the embryos are formed over a wide range of radial distances, the super-Earths would end up having a distinctive, bimodal composition in water mass. Our model succeeds in quantitatively reproducing several important observed properties of exoplanets and correlations with their stellar hosts.

Gas-giant planets, like Jupiter and Saturn, acquire massive gaseous envelopes during the approximately 3 Myr-long lifetimes of protoplanetary discs. In the core accretion scenario, the formation of a solid core of around 10 Earth masses triggers a phase of rapid gas accretion. Previous 3D grid-based hydrodynamical simulations found runaway gas accretion rates corresponding to approximately 10 to 100 Jupiter masses per Myr. Such high accretion rates would result in all planets with larger-than-10-Earth-mass cores forming Jupiter-like planets, in clear contrast to the ice giants in the Solar System and the observed exoplanet population. In this work, we use 3D hydrodynamical simulations, that include radiative transfer, to model the growth of the envelope on planets with different masses. We find that gas flows rapidly through the outer part of the envelope, but this flow does not drive accretion. Instead, gas accretion is the result of quasi-static contraction of the inner envelope, which can be orders of magnitude smaller than the mass flow through the outer atmosphere. For planets smaller than Saturn, we measure moderate gas accretion rates that are below 1 Jupiter mass per Myr. Higher mass planets, however, accrete up to 10 times faster and do not reveal a self-driven mechanism that can halt gas accretion. Therefore, the reason for the final masses of Saturn and Jupiter remains difficult to understand, unless their completion coincided with the dissipation of the Solar Nebula.

Context. The dynamics of a low-mass protoplanet accreting solids is influenced by the heating torque, which was found to suppress inward migration in protoplanetary disks with constant opacities. Aims. We investigate the differences of the heating torque between disks with constant and temperature-dependent opacities. Methods. Interactions of a super-Earth-sized protoplanet with the gas disk are explored using 3D radiation hydrodynamic simulations. Results. Accretion heating of the protoplanet creates a hot underdense region in the surrounding gas, leading to misalignment of the local density and pressure gradients. As a result, the 3D gas flow is perturbed and some of the streamlines form a retrograde spiral rising above the protoplanet. In the constant-opacity disk, the perturbed flow reaches a steady state and the underdense gas responsible for the heating torque remains distributed in accordance with previous studies. If the opacity is non-uniform, however, the differences in the disk structure can lead to more vigorous streamline distortion and eventually to a flow instability. The underdense gas develops a one-sided asymmetry which circulates around the protoplanet in a retrograde fashion. The heating torque thus strongly oscillates in time and does not on average counteract inward migration. Conclusions. The torque variations make the radial drift of the protoplanet oscillatory, consisting of short intervals of alternating rapid inward and outward migration. We speculate that transitions between the positive and oscillatory heating torque may occur in specific disk regions susceptible to vertical convection, resulting in the convergent migration of multiple planetary embryos.

Giant planets migrate though the protoplanetary disc as they grow. We investigate how the formation of planetary systems depends on the radial flux of pebbles through the protoplanetary disc and on the planet migration rate. Our N-body simulations confirm previous findings that Jupiter-like planets in orbits outside the water ice line originate from embryos starting out at 20-40 AU when using nominal type-I and type-II migration rates and a pebble flux of 100-200 Earth masses per million years, enough to grow Jupiter within the lifetime of the solar nebula. The planetary embryos placed up to 30AU migrate into the inner system (r<1AU) and form super-Earths or hot and warm gas giants, producing systems that are inconsistent with the configuration of the solar system, but consistent with some exoplanetary systems. We also explore slower migration rates which allow the formation of gas giants from embryos originating from the 5-10AU region, which are stranded exterior to 1 AU at the end of the gas-disc phase. We identify a pebble flux threshold below which migration dominates and moves the planetary core to the inner disc, where the pebble isolation mass is too low for the planet to accrete gas efficiently. Giant planet growth requires a sufficiently-high pebble flux to enable growth to out-compete migration. Even higher pebble fluxes produce systems with multiple gas giants. We show that planetary embryos starting interior to 5AU do not grow into gas giants, even if migration is slow and the pebble flux is large. Instead they grow to the mass regime of super-Earths. This stunted growth is caused by the low pebble isolation mass in the inner disc and is independent of the pebble flux. Additionally we show that the long term evolution of our formed planetary systems can produce systems with hot super-Earths and outer gas giants as well as systems of giants on eccentric orbits (abridged).

At least 30% of main sequence stars host planets with sizes between 1 and 4 Earth radii and orbital periods of less than 100 days. We use N-body simulations including a model for gas-assisted pebble accretion and disk-planet tidal interaction to study the formation of super-Earth systems. We show that the integrated pebble mass reservoir creates a bifurcation between hot super-Earths or hot-Neptunes ($\lesssim15M_{\oplus}$) and super-massive planetary cores potentially able to become gas giant planets ($\gtrsim15M_{\oplus}$). Simulations with moderate pebble fluxes grow multiple super-Earth-mass planets that migrate inwards and pile up at the disk's inner edge forming long resonant chains. We follow the long-term dynamical evolution of these systems and use the period ratio distribution of observed planet-pairs to constrain our model. Up to $\sim$95% of resonant chains become dynamically unstable after the gas disk dispersal, leading to a phase of late collisions that breaks the resonant configuration. Our simulations match observations if we combine a dominant fraction ($\gtrsim95\%$) of unstable systems with a sprinkling ($\lesssim5\%$) of stable resonant chains (the Trappist-1 system represents one such example). Our results demonstrate that super-Earth systems are inherently multiple (${\rm N\geq2}$) and that the observed excess of single-planet transits is a consequence of the mutual inclinations excited by the planet-planet instability. In simulations in which planetary seeds are initially distributed in the inner and outer disk, close-in super-Earths are systematically ice-rich. This contrasts with the interpretation that most super-Earths are rocky based on bulk density measurements of super-Earths and photo-evaporation modeling of their bimodal radius distribution. We investigate the conditions needed to form rocky super-Earths. The formation of rocky super-Earths (abridged)

Super-Earths are found in tighter orbits than the Earth's around more than one third of main sequence stars. It has been proposed that super-Earths are scaled-up terrestrial planets that formed similarly, through mutual accretion of planetary embryos, but in discs much denser than the solar protoplanetary disc. We argue instead that terrestrial planets and super-Earths have two distinct formation pathways that are regulated by the disc's pebble reservoir. Through numerical integrations, which combine pebble accretion and N-body gravity between embryos, we show that a difference of a factor of two in the pebble mass-flux is enough to change the evolution from the terrestrial to the super-Earth growth mode. If the pebble mass-flux is small, then the initial embryos within the ice line grow slowly and do not migrate substantially, resulting in a widely spaced population of Mars-mass embryos when the gas disc dissipates. Without gas being present, the embryos become unstable and a small number of terrestrial planets are formed by mutual collisions. The final terrestrial planets are at most 5 Earth masses. Instead, if the pebble mass-flux is high, then the initial embryos within the ice line rapidly become sufficiently massive to migrate through the gas disc. Embryos concentrate at the inner edge of the disc and growth accelerates through mutual merging. This leads to the formation of a system of closely spaced super-Earths in the 5 to 20 Earth-mass range, bounded by the pebble isolation mass. Generally, instabilities of these super-Earth systems after the disappearance of the gas disc trigger additional merging events and dislodge the system from resonant chains. The pebble flux - which controls the transition between the two growth modes - may be regulated by the initial reservoir of solids in the disc or the presence of more distant giant planets that can halt the radial flow of pebbles.

Protoplanets of Super-Earth sizes may get trapped in convergence zones for planetary migration and form gas giants there. These growing planets undergo accretion heating, which triggers a hot-trail effect that can reverse migration directions, increase eccentricities and prevent resonant captures (Chrenko et al. 2017). We study populations of embryos accreting pebbles using Fargo-Thorin 2D hydrocode. We find that embryos in a disk with high surface density ($\Sigma_0 = 990\,{\rm g}\,{\rm cm}^{-2}$) undergo `unsuccessful' two-body encounters which do not lead to a merger. Only when a 3rd protoplanet arrives to the convergence zone, three-body encounters lead to mergers. For a low-viscosity disk ($\nu = 5\times10^{13}\,{\rm cm}^2\,{\rm s}^{-1}$) a massive coorbital is a possible outcome, for which a pebble isolation develops and the coorbital is stabilised. For more massive protoplanets ($5\,M_\oplus$), the convergence radius is located further out, in the ice-giant zone. After a series of encounters, there is an evolution driven by a dynamical torque of a tadpole region, which is systematically repeated several times, until the coorbital configuration is disrupted and planets merge. This may be a pathway how to solve the problem that coorbitals often form in simulations but they are not observed in nature. In contrast, the joint evolution of 120 low-mass protoplanets ($0.1\,M_\oplus$) reveals completely different dynamics. The evolution is no longer smooth, but rather a random walk. This is because the spiral arms, developed in the gas disk due to Lindblad resonances, overlap with each other and affect not only a single protoplanet but several in the surroundings. Our hydrodynamical simulations may have important implications for N-body simulations of planetary migration that use simplified torque prescriptions and are thus unable to capture protoplanet dynamics in its full glory.

Forming macroscopic solid bodies in circumstellar discs requires local dust concentration levels significantly higher than the mean. Interactions of the dust particles with the gas must serve to augment local particle densities, and facilitate growth past barriers in the metre size range. Amongst a number of mechanisms that can amplify the local density of solids, aerodynamic streaming instability (SI) is one of the most promising. This work tests the physical assumptions of models that lead to SI in protoplanetary disks (PPDs). We conduct laboratory experiments in which we track the three-dimensional motion of spherical solid particles fluidized in a low-pressure, laminar, incompressible, gas stream. The particle sizes span the Stokes-Epstein drag regime transition and the overall dust-to-gas mass density ratio is close to unity. Lambrechts et al. (2016) established the similarity of the laboratory flow to a simplified PPD model flow. We obtain experimental results suggesting an instability due to particle-gas interaction: i) there exist variations in particle concentration in the direction of the mean drag forces; ii) the particles have a tendency to 'catch up' to one another when they are in proximity; iii) particle clumping occurs on very small scales, which implies local enhancements above the background dust-to-gas mass density ratio by factors of several tens; v) the presence of these density enhancements occurs for a mean dust-to-gas mass density ratio approaching or greater than 1; v) we find evidence for collective particle drag reduction when the local particle number density becomes high and when the background gas pressure is high so that the drag is in the continuum regime. The experiments presented here are precedent-setting for observing SI under controlled conditions and may lead to a deeper understanding of how it operates in nature.

The Wide Field Infrared Survey Telescope (WFIRST) was the top ranked large space mission in the 2010 New Worlds, New Horizons decadal survey, and it was formed by merging the science programs of 3 different mission concepts, including the Microlensing Planet Finder (MPF) concept (Bennett \etal\ 2010). The WFIRST science program (Spergel \etal\ 2015) consists of a general observer program, a wavefront controlled technology program, and two targeted science programs: a program to study dark energy, and a statistical census of exoplanets with a microlensing survey, which uses nearly one quarter of WFIRST's observing time in the current design reference mission. The New Worlds, New Horizons (decadal survey) midterm assessment summarizes the science case for the WFIRST exoplanet microlensing survey with this statement: "WFIRST's microlensing census of planets beyond 1 AU will perfectly complement Kepler's census of compact systems, and WFIRST will also be able to detect free-floating planets unbound from their parent stars\rlap."

The growth of a planetary core by pebble accretion stops at the so called pebble isolation mass, when the core generates a pressure bump that traps drifting pebbles outside its orbit. If the isolation mass is very small, then gas accretion is protracted and the planet remains at a few Earth masses with a mainly solid composition. For larger values of the pebble isolation mass, the planet might be able to accrete gas from the protoplanetary disc and grow into a gas giant. Previous works have determined a scaling of the pebble isolation mass with cube of the disc aspect ratio. Here we expand on previous measurements and explore the dependency of the pebble isolation mass on all relevant parameters of the protoplanetary disc. We use 3D hydrodynamical simulations to measure the pebble isolation mass and derive a simple scaling law that captures the dependence on the local disc structure and the turbulent viscosity parameter $\alpha$. We find that small pebbles, coupled to the gas, with Stokes number $\tau_{\rm f}<0.005$ can drift through the partial gap at pebble isolation mass. However, as the planetary mass increases, particles must be decreasingly smaller to penetrate through the pressure bump. Turbulent diffusion of particles, however, can lead to an increase of the pebble isolation mass by a factor of two, depending on the strength of the background viscosity and on the pebble size. We finally explore the implications of the new scaling law of the pebble isolation mass on the formation of planetary systems by numerically integrating the growth and migration pathways of planets in evolving protoplanetary discs. Compared to models neglecting the dependence of the pebble isolation mass on the $\alpha$-viscosity, our models including this effect result in larger core masses for giant planets. These larger core masses are more akin to the core masses of the giant planets in the Solar System.

A large fraction of giant planets have gaseous envelopes that are limited to about 10 % of their total mass budget. Such planets are present in the Solar System (Uranus, Neptune) and are frequently observed in short periods around other stars (the so-called Super-Earths). In contrast to these observations, theoretical calculations based on the evolution of hydrostatic envelopes argue that such low mass envelopes cannot be maintained around cores exceeding five Earth masses. Instead, under nominal disc conditions, these planets would acquire massive envelopes through runaway gas accretion within the lifetime of the protoplanetary disc. In this work, we show that planetary envelopes are not in hydrostatic balance, which slows down envelope growth. A series of 3-dimensional, global, radiative hydrodynamical simulations reveal a steady state gas flow, which enters through the poles and exits in the disc midplane. Gas is pushed through the outer envelope in about 10 orbital timescales. In regions of the disc that are not significantly dust-depleted, envelope accretion onto cores of about five Earth masses can get stalled as the gas flow enters the deep interior. Accreted solids sublimate deep in the convective interior, but small opacity-providing grains are trapped in the flow and do not settle, which further prevents rapid envelope accretion. The transition to runaway gas accretion can however be reached when cores grow larger than typical Super-Earths, beyond 15 Earth masses, and preferably when disc opacities are below kappa=1 cm^2/g. These findings offer an explanation for the typical low-mass envelopes around the cores of Super-Earths.

Context: Planetary embryos can continue to grow by pebble accretion until they become giant planet cores. Simultaneously, these embryos mutually interact and also migrate due to torques arising from the protoplanetary disk. Aims: Our aim is to investigate how pebble accretion alters the orbital evolution of embryos undergoing the Type-I migration. In particular, we study whether they establish resonant chains, whether these chains are prone to instabilities and if giant planet cores form through embryo merging, thus occurring more rapidly than by pebble accretion alone. Methods: For the first time, we perform self-consistent global-scale radiative hydrodynamic simulations of a two-fluid protoplanetary disk consisting of gas and pebbles, the latter being accreted by embedded embryos. Accretion heating, along with other radiative processes, is accounted for to correctly model the Type-I migration. Results: We track the evolution of four super-Earth-like embryos, initially located in a region where the disk structure allows for a convergent migration. Generally, embryo merging is facilitated by rapidly increasing embryo masses and breaks the otherwise oligarchic growth. Moreover, we find that the orbital eccentricity of each embryo is considerably excited ($\simeq0.03$) due to the presence of an asymmetric underdense lobe of gas, a so-called `hot trail', produced by accretion heating of the embryo's vicinity. Eccentric orbits lead the embryos to frequent close encounters and make resonant locking more difficult. Conclusions: Embryo merging typically produces one massive core ($\gtrsim 10\,\mathrm{M_{E}}$) in our simulations, orbiting near $10\,\mathrm{AU}$. Pebble accretion is naturally accompanied by occurrence of eccentric orbits which should be considered in future efforts to explain the structure of exoplanetary systems.

The formation of planetary cores must proceed rapidly in order for the giant planets to accrete their gaseous envelopes before the dissipation of the protoplanetary gas disc (<3 Myr). In orbits beyond 10 AU, direct accumulation of planetesimals by the cores is too slow. Fragments of planetesimals could be accreted faster, but planetesimals are likely too large for fragmentation to be efficient, and resonant trapping poses a further hurdle. Here we instead investigate the accretion of small pebbles (mm-cm sizes) that are the natural outcome of an equilibrium between the growth and radial drift of particles. We construct a simplified analytical model of dust coagulation and pebble drift in the outer disc, between 5 AU and 100 AU, which gives the temporal evolution of the solid surface density and the dominant particle size. These two key quantities determine how core growth proceeds at various orbital distances. We find that pebble surface densities are sufficiently high to achieve the inside-out formation of planetary cores within the disc lifetime. The overall efficiency by which dust gets converted to planets can be high, close to 50 % for planetary architectures similar to the Solar System. Growth by pebble accretion in the outer disc is sufficiently fast to overcome catastrophic Type I migration of the cores. These results require protoplanetary discs with large radial extent (~100 AU) and assume a low number of initial seed embryos. Our findings imply that protoplanetary discs with low disc masses, as expected around low-mass stars (<1 M_sun), or with sub-solar dust-to-gas ratios, do not easily form gas-giant planets (M > 100 M_E), but preferentially form Neptune-mass planets or smaller (< 10 M_E). This is consistent with exoplanet surveys which show that gas giants are relatively uncommon around stars of low mass or low metallicity.

The behaviour of sedimenting particles depends on the dust-to-gas ratio of the fluid. Linear stability analysis shows that solids settling in the Epstein drag regime would remain homogeneously distributed in non-rotating incompressible fluids, even when dust-to-gas ratios reach unity. However, the non-linear evolution has not been probed before. Here, we present numerical calculations indicating that in a particle-dense mixture solids spontaneously mix out of the fluid and form swarms overdense in particles by at least a factor 10. The instability is caused by mass-loaded regions locally breaking the equilibrium background stratification. The driving mechanism depends on non-linear perturbations of the background flow and shares some similarity to the streaming instability in accretion discs. The resulting particle-rich swarms may stimulate particle growth by coagulation. In the context of protoplanetary discs, the instability could be relevant for aiding small particles to settle to the midplane in the outer disc. Inside the gas envelopes of protoplanets, enhanced settling may lead to a reduced dust opacity, which facilitates the contraction of the envelope. We show that the relevant physical set up can be recreated in a laboratory setting. This will allow our numerical calculations to be investigated experimentally in the future.

The terrestrial planets and the asteroids dominant in the inner asteroid belt are water poor. However, in the protoplanetary disk the temperature should have decreased below water condensation level well before the disk was photoevaporated. Thus, the global water depletion of the inner Solar System is puzling. We show that, even if the inner disk becomes cold, there cannot be direct condensation of water. This is because the snowline moves towards the Sun more slowly than the gas itself. The appearance of ice in a range of heliocentric distances swept by the snowline can only be due to the radial drift of icy particles from the outer disk. However, if a sufficiently massive planet is present, the radial drift of particles is interrupted, because the disk acquires a superKeplerian rotation just outside of the planetary orbit. From this result, we propose that the precursor of Jupiter achieved about 20 Earth masses when the snowline was still around 3 AU. This effectively fossilized the snowline at that location. Although cooling, the disk inside of the Jovian orbit remained ice-depleted because the flow of icy particles from the outer system was intercepted by the planet. This scenario predicts that planetary systems without giant planets should be much more rich in water in their inner regions than our system. We also show that the inner edge of the planetesimal disk at 0.7AU, required in terrestrial planet formation models to explain the small mass of Mercury and the absence of planets inside of its orbit, could be due to the silicate condensation line, fossilized at the end of the phase of streaming instability that generated the planetesimal seeds. Thus, when the disk cooled, silicate particles started to drift inwards of 0.7AU without being sublimated, but they could not be accreted by any pre-existing planetesimals.

The formation of planets depends on the underlying protoplanetary disc structure, which influences both the accretion and migration rates of embedded planets. The disc itself evolves on time-scales of several Myr during which both temperature and density profiles change as matter accretes onto the central star. Here we use a detailed model of an evolving disc to determine the growth of planets by pebble accretion and their migration through the disc. Cores that reach their pebble isolation mass accrete gas to finally form giant planets with extensive gas envelopes, while planets that do not reach pebble isolation mass are stranded as ice giants and ice planets containing only minor amounts of gas in their envelopes. Unlike earlier population synthesis models, our model works without any artificial reductions in migration speed and for protoplanetary discs with gas and dust column densities similar to those inferred from observations. We find that in our nominal disc model the emergence of planetary embryos preferably occurs after approximately 2 Myr in order to not exclusively form gas giants, but also ice giants and smaller planets. The high pebble accretion rates ensure that critical core masses for gas accretion can be reached at all orbital distances. Gas giant planets nevertheless experience significant reduction in semi-major axes by migration. Considering instead planetesimal accretion for planetary growth, we show that formation time-scales are too long to compete with the migration time-scales and the dissipation time of the protoplanetary disc. Altogether, we find that pebble accretion overcomes many of the challenges in the formation of ice and gas giants in evolving protoplanetary discs.

The basic structure of the solar system is set by the presence of low-mass terrestrial planets in its inner part and giant planets in its outer part. This is the result of the formation of a system of multiple embryos with approximately the mass of Mars in the inner disk and of a few multi-Earth-mass cores in the outer disk, within the lifetime of the gaseous component of the protoplanetary disk. What was the origin of this dichotomy in the mass distribution of embryos/cores? We show in this paper that the classic processes of runaway and oligarchic growth from a disk of planetesimals cannot explain this dichotomy, even if the original surface density of solids increased at the snowline. Instead, the accretion of drifting pebbles by embryos and cores can explain the dichotomy, provided that some assumptions hold true. We propose that the mass-flow of pebbles is two-times lower and the characteristic size of the pebbles is approximately ten times smaller within the snowline than beyond the snowline (respectively at heliocentric distance $r<r_{ice}$ and $r>r_{ice}$, where $r_{ice}$ is the snowline heliocentric distance), due to ice sublimation and the splitting of icy pebbles into a collection of chondrule-size silicate grains. In this case, objects of original sub-lunar mass would grow at drastically different rates in the two regions of the disk. Within the snowline these bodies would reach approximately the mass of Mars while beyond the snowline they would grow to $\sim 20$ Earth masses. The results may change quantitatively with changes to the assumed parameters, but the establishment of a clear dichotomy in the mass distribution of protoplanets appears robust, provided that there is enough turbulence in the disk to prevent the sedimentation of the silicate grains into a very thin layer.

The formation of planets with gaseous envelopes takes place in protoplanetary accretion discs on time-scales of several millions of years. Small dust particles stick to each other to form pebbles, pebbles concentrate in the turbulent flow to form planetesimals and planetary embryos and grow to planets, which undergo substantial radial migration. All these processes are influenced by the underlying structure of the protoplanetary disc, specifically the profiles of temperature, gas scale height and density. The commonly used disc structure of the Minimum Mass Solar Nebular (MMSN) is a simple power law in all these quantities. However, protoplanetary disc models with both viscous and stellar heating show several bumps and dips in temperature, scale height and density caused by transitions in opacity, which are missing in the MMSN model. These play an important role in the formation of planets, as they can act as sweet spots for the formation of planetesimals via the streaming instability and affect the direction and magnitude of type-I-migration. We present 2D simulations of accretion discs that feature radiative cooling, viscous and stellar heating, and are linked to the observed evolutionary stages of protoplanetary discs and their host stars. These models allow us to identify preferred planetesimal and planet formation regions in the protoplanetary disc as a function of the disc's metallicity, accretion rate and lifetime. We derive simple fitting formulae that feature all structural characteristics of protoplanetary discs during the evolution of several Myr. These fits are straightforward to apply for modelling any growth stage of planets where detailed knowledge of the underlying disc structure is required.

In the Solar System giant planets come in two flavours: 'gas giants' (Jupiter and Saturn) with massive gas envelopes and 'ice giants' (Uranus and Neptune) with much thinner envelopes around their cores. It is poorly understood how these two classes of planets formed. High solid accretion rates, necessary to form the cores of giant planets within the life-time of protoplanetary discs, heat the envelope and prevent rapid gas contraction onto the core, unless accretion is halted. We find that, in fact, accretion of pebbles (~ cm-sized particles) is self-limiting: when a core becomes massive enough it carves a gap in the pebble disc. This halt in pebble accretion subsequently triggers the rapid collapse of the super-critical gas envelope. As opposed to gas giants, ice giants do not reach this threshold mass and can only bind low-mass envelopes that are highly enriched by water vapour from sublimated icy pebbles. This offers an explanation for the compositional difference between gas giants and ice giants in the Solar System. Furthermore, as opposed to planetesimal-driven accretion scenarios, our model allows core formation and envelope attraction within disc life-times, provided that solids in protoplanetary discs are predominantly in pebbles. Our results imply that the outer regions of planetary systems, where the mass required to halt pebble accretion is large, are dominated by ice giants and that gas-giant exoplanets in wide orbits are enriched by more than 50 Earth masses of solids.

The observed lifetimes of gaseous protoplanetary discs place strong constraints on gas and ice giant formation in the core accretion scenario. The approximately 10-Earth-mass solid core responsible for the attraction of the gaseous envelope has to form before gas dissipation in the protoplanetary disc is completed within 1-10 million years. Building up the core by collisions between km-sized planetesimals fails to meet this time-scale constraint, especially at wide stellar separations. Nonetheless, gas-giant planets are detected by direct imaging at wide orbital distances. In this paper, we numerically study the growth of cores by the accretion of cm-sized pebbles loosely coupled to the gas. We measure the accretion rate onto seed masses ranging from a large planetesimal to a fully grown 10-Earth-mass core and test different particle sizes. The numerical results are in good agreement with our analytic expressions, indicating the existence of two accretion regimes, one set by the azimuthal and radial particle drift for the lower seed masses and the other, for higher masses, by the velocity at the edge of the Hill sphere. In the former, the optimally accreted particle size increases with core mass, while in the latter the optimal size is centimeters, independent of core mass. We discuss the implications for rapid core growth of gas-giant and ice-giant cores. We conclude that pebble accretion can resolve the long-standing core accretion time-scale conflict. This requires a near-unity dust-to-gas ratio in the midplane, particle growth to mm and cm and the formation of massive planetesimals or low radial pressure support. The core growth time-scale is shortened by a factor 30-1,000 at 5 AU and by a factor 100-10,000 at 50 AU, compared to the gravitationally focused accretion of, respectively, low-scale-height planetesimal fragments or standard km-sized planetesimals.