Aims: ARCiS, a novel code for the analysis of exoplanet transmission and emission spectra is presented. The aim of the modelling framework is to provide a tool able to link observations to physical models of exoplanet atmospheres. Methods: The modelling philosophy chosen in this paper is to use physical and chemical models to constrain certain parameters while keeping free the parts where our physical understanding is still more limited. This approach, in between full physical modelling and full parameterisation, allows us to use the processes we understand well and parameterise those less understood. A Bayesian retrieval framework is implemented and applied to the transit spectra of a set of 10 hot Jupiters. The code contains chemistry and cloud formation and has the option for self consistent temperature structure computations. Results: The code presented is fast and flexible enough to be used for retrieval and for target list simulations for e.g. JWST or the ESA Ariel missions. We present results for the retrieval of elemental abundance ratios using the physical retrieval framework and compare this to results obtained using a parameterised retrieval setup. Conclusions: We conclude that for most of the targets considered the current dataset is not constraining enough to reliably pin down the elemental abundance ratios. We find no significant correlations between different physical parameters. We confirm that planets in our sample with a strong slope in the optical transmission spectrum are the planets where we find cloud formation to be most active. Finally, we conclude that with ARCiS we have a computationally efficient tool to analyse exoplanet observations in the context of physical and chemical models.

Proto-planets embedded in their natal disks acquire hot envelopes as they grow and accrete solids. This ensures that the material they accrete - pebbles, as well as (small) planetesimals - will vaporize to enrich their atmospheres. Enrichment modifies an envelope's structure and significantly alters its further evolution. Our aim is to describe the formation of planets with polluted envelopes from the moment that impactors begin to sublimate to beyond the disk's eventual dissipation. We constructed an analytical interior structure model, characterized by a hot and uniformly mixed high-Z vapor layer surrounding the core, located below the usual unpolluted radiative-convective regions. The evolution of planets with uniformly mixed polluted envelopes follows four potential phases. Initially, the central core grows directly through impacts and rainout until the envelope becomes hot enough to vaporize and absorb all incoming solids. We find that a planet reaches runaway accretion when the sum of its core and vapor mass exceeds a value that we refer to as the critical metal mass - a criterion that supersedes the traditional critical core mass. It scales positively with both the pollutant's evaporation temperature and with the planet's core mass. Hence, planets at shorter orbital separations require the accretion of more solids to reach runaway as they accrete less volatile materials. If the solids accretion rate dries up, we identify the decline of the mean molecular weight - dilution - as a mechanism to limit gas accretion during a polluted planet's embedded cooling phase. When the disk ultimately dissipates, the envelope's inner temperature declines and its vapor eventually rains out, augmenting the mass of the core. The energy release that accompanies this does not result in significant mass-loss, as it only occurs after the planet has substantially contracted.

The detection of Earth-size exoplanets around low-mass stars -- in stars such as Proxima Centauri and TRAPPIST-1 -- provide an exceptional chance to improve our understanding of the formation of planets around M stars and brown dwarfs. We explore the formation of such planets with a population synthesis code based on a planetesimal-driven model previously used to study the formation of the Jovian satellites. Because the discs have low mass and the stars are cool, the formation is an inefficient process that happens at short periods, generating compact planetary systems. Planets can be trapped in resonances and we follow the evolution of the planets after the gas has dissipated and they undergo orbit crossings and possible mergers. We find that formation of planets above Mars mass and in the planetesimal accretion scenario, is only possible around stars with masses $M_{\star} \ge 0.07 M_{sun}$ and discs of $M_{disc} \ge 10^{-2}~M_{sun} $. We find that planets above Earth-mass form around stars with masses larger than $0.15 M_{\oplus}$ while planets larger than $5 M_{\oplus}$do not form in our model, even not under the most optimal conditions (massive disc), showing that planets such as GJ 3512b form with another, more efficient mechanism. Our results show that the majority of planets form with a significant water fraction; that most of our synthetic planetary systems have 1, 2, or 3 planets, but those with 4, 5, 6, and 7 planets are also common, confirming that compact planetary systems with many planets should be a relatively common outcome of planet formation around small stars.

It is generally accepted that the four major (Galilean) satellites formed out of the gas disk that accompanied Jupiter's formation. However, understanding the specifics of the formation process is challenging as both small particles (pebbles) as well as the satellites are subject to fast migration processes. Here, we hypothesize a new scenario for the origin of the Galilean system, based on the capture of several planetesimal seeds and subsequent slow accretion of pebbles. To halt migration, we invoke an inner disk truncation radius, and other parameters are tuned for the model to match physical, dynamical, compositional, and structural constraints. In our scenario it is natural that Ganymede's mass is determined by pebble isolation. Our slow-pebble-accretion scenario then reproduces the following characteristics: (1) the mass of all the Galilean satellites; (2) the orbits of Io, Europa, and Ganymede captured in mutual 2:1 mean motion resonances; (3) the ice mass fractions of all the Galilean satellites; (4) the unique ice-rock partially differentiated Callisto and the complete differentiation of the other satellites. Our scenario is unique to simultaneously reproduce these disparate properties.

Context: Planet formation by pebble accretion is an alternative to planetesimal-driven core accretion. In this scenario, planets grow by accreting cm-to-m-sized pebbles instead of km-sized planetesimals. One of the main differences with planetesimal-driven core accretion is the increased thermal ablation experienced by pebbles. This provides early enrichment to the planet's envelope, which changes the process of core growth. Aims: We aim to predict core masses and envelope compositions of planets that form by pebble accretion and compare mass deposition of pebbles to planetesimals. Methods: We model the early growth of a proto-planet by calculating the structure of its envelope, taking into account the fate of impacting pebbles or planetesimals. The region where high-Z material can exist in vapor form is determined by the vapor pressure. We include enrichment effects by locally modifying the mean molecular weight. Results: In the pebble case, three phases of core growth can be identified. In the first phase, pebbles impact the core without significant ablation. During the second phase, ablation becomes increasingly severe. A layer of high-Z vapor starts to form around the core that absorbs a small fraction of the ablated mass. The rest either rains out to the core or mixes outwards instead, slowing core growth. In the third phase, the vapor inner region expands outwards, absorbing an increasing fraction of the ablated material as vapor. Rainout ends before the core mass reaches 0.6 M_Earth, terminating direct core growth. In the case of icy H2O pebbles, this happens before 0.1 M_Earth. Conclusions: Our results indicate that pebble accretion can directly form rocky cores up to only 0.6 M_Earth, and is unable to form similarly sized icy cores. Subsequent core growth can proceed indirectly when the planet cools, provided it is able to retain its high-Z material.

Most major planetary bodies in the solar system rotate in the same direction as their orbital motion: their spin is prograde. Theoretical studies to explain the direction as well as the magnitude of the spin vector have had mixed success. When the accreting building blocks are $\sim$ km-size planetesimals -- as predicted by the classical model -- the accretion process is so symmetric that it cancels out prograde with retrograde spin contributions, rendering the net spin minute. For this reason, the currently-favored model for the origin of planetary rotation is the giant impact model, in which a single collision suffices to deliver a spin, which magnitude is close to the breakup rotation rate. However, the giant impact model does not naturally explain the preference for prograde spin. Similarly, an increasing number of spin-vector measurement of asteroids also shows that the spin vector of large (primordial) asteroids is not isotropic. Here, we re-assess the viability of smaller particles to bestow planetary bodies with a net spin, focusing on the pebble accretion model in which gas drag and gravity join forces to accrete small particles at a large cross section. Similar to the classical calculation for planetesimals, we integrate the pebble equation of motion and measure the angular momentum transfer at impact. We consider a variety of disk conditions and pebble properties and conduct our calculations in the limits of 2D (planar) and 3D (homogeneous) pebble distributions. We find that in certain regions of the parameter space the angular momentum transfer is significant, much larger than with planetesimals and on par with or exceeding the current spin of planetary bodies.

Recently, seven Earth-sized planets were discovered around the M-dwarf star TRAPPIST-1. Thanks to transit-timing variations, the masses and therefore the bulk densities of the planets have been constrained, suggesting that all TRAPPIST-1 planets are consistent with water mass fractions on the order of 10%. These water fractions, as well as the similar planet masses within the system, constitute strong constraints on the origins of the TRAPPIST-1 system. In a previous work, we outlined a pebble-driven formation scenario. In this paper we investigate this formation scenario in more detail. We used a Lagrangian smooth-particle method to model the growth and drift of pebbles and the conversion of pebbles to planetesimals through the streaming instability. We used the N-body code \texttt{MERCURY} to follow the composition of planetesimals as they grow into protoplanets by merging and accreting pebbles. This code is adapted to account for pebble accretion, type-I migration, and gas drag. In this way, we modelled the entire planet formation process (pertaining to planet masses and compositions, not dynamical configuration). We find that planetesimals form in a single, early phase of streaming instability. The initially narrow annulus of planetesimals outside the snowline quickly broadens due to scattering. Our simulation results confirm that this formation pathway indeed leads to similarly-sized planets and is highly efficient in turning pebbles into planets ($\sim$50% solids-to-planets conversion efficiency). [...] The water content of planets resulting from our simulations is on the order of 10%, and our results predict a `V-shaped' trend in the planet water fraction with orbital distance: from relatively high (innermost planets) to relatively low (intermediate planets) to relatively high (outermost planets).

Streaming instability is a key mechanism in planet formation, clustering pebbles into planetesimals. It is triggered at a particular disk location where the local volume density of solids exceeds that of the gas. After their formation, planetesimals can grow by feeding from other planetesimals in the birth ring as well as by accreting inwardly drifting pebbles from the outer disk. To investigate the growth of planetesimals at a single location by the streaming instability, we test the conditions under which super-Earths are able to form within the lifetime of the gaseous disk. We modify the \texttt{Mercury} N-body code to trace the growth and dynamical evolution of a swarm of planetesimals at the ice line for a solar-mass star. Three distributions of planetesimal sizes are investigated: (i) a mono-dispersed population of 400 km radius planetesimals, (ii) a poly-dispersed populations of planetesimals from 200 km up to 1000 km, (iii) a bimodal distribution with a single runaway body and a swarm of smaller, 100 km size planetesimals. The mono-disperse population of 400 km size planetesimals cannot form $\gtrsim$ Earth mass protoplanets. Their velocity dispersions are quickly excited, which suppresses both planetesimal and pebble accretion. Planets can form from the poly-dispersed and bimodal distributions. In these circumstances, the two-component nature damps the random velocity of the large embryo by small planetesimals' dynamical friction, allowing the embryo to accrete pebbles efficiently when it approaches $10^{-2}$ Earth mass. We find that super-Earth planets are preferred to form when the pebble mass flux is higher, the disk turbulence is lower, or the Stokes number of the pebbles is higher.

Understanding of clouds is instrumental in interpreting current and future spectroscopic observations of exoplanets. Modelling clouds consistently is complex, since it involves many facets of chemistry, nucleation theory, condensation physics, coagulation, and particle transport. We develop a simple physical model for cloud formation and transport, efficient and versatile enough that it can be used in modular fashion for parameter optimization searches of exoplanet atmosphere spectra. The transport equations are formulated in 1D, accounting for sedimentation and diffusion. The grain size is obtained through a moment method. For simplicity, only one cloud species is considered and the nucleation rate is parametrized. From the resulting physical profiles we simulate transmission spectra covering the visual to mid-IR wavelength range. We apply our models towards KCl clouds in the atmosphere of GJ1214 b and towards MgSiO3 clouds of a canonical hot-Jupiter. We find that larger cloud diffusivity $K_{zz}$ increases the thickness of the cloud, pushing the $\tau=1$ surface to a lower pressure layer higher in the atmosphere. A larger nucleation rate also increases the cloud thickness while it suppresses the grain size. Coagulation is most important at high nuclei injection rates ($\dot\Sigma_n$) and low $K_{zz}$. We find that the investigated combinations of $K_{zz}$ and $\dot\Sigma_n$ greatly affect the transmission spectra in terms of the slope at near-IR wavelength (a proxy for grain size), the molecular features seen at ~1\micr (which disappear for thick clouds, high in the atmosphere), and the 10\micr silicate feature, which becomes prominent for small grains high in the atmosphere. The result of our hybrid approach -- aimed to provide a good balance between physical consistency and computational efficiency -- is ideal towards interpreting (future) spectroscopic observations of exoplanets.

Sub-Neptune planets are a very common type of planets. They are inferred to harbour a primordial (H/He) envelope, on top of a (rocky) core, which dominates the mass. Here, we investigate the long-term consequences of the core properties on the planet mass-radius relation. We consider the role of various core energy sources resulting from core formation, its differentiation, its solidification (latent heat), core contraction and radioactive decay. We divide the evolution of the rocky core into three phases: the formation phase, which sets the initial conditions, the magma ocean phase, characterized by rapid heat transport, and the solid state phase, where cooling is inefficient. We find that for typical sub-Neptune planets of ~2-10 Earth masses and envelope mass fractions of 0.5-10% the magma ocean phase lasts several Gyrs, much longer than for terrestrial planets. The magma ocean phase effectively erases any signs of the initial core thermodynamic state. After solidification, the reduced heat flux from the rocky core causes a significant drop in the rocky core surface temperature, but its effect on the planet radius is limited. In the long run, radioactive heating is the most significant core energy source in our model. Overall, the long term radius uncertainty by core thermal effects is up to 15%.

We introduce a new Lagrangian smooth-particle method to model the growth and drift of pebbles in protoplanetary disks. The Lagrangian nature of the model makes it especially suited to follow characteristics of individual (groups of) particles, such as their composition. In this work we focus on the water content of solid particles. Planetesimal formation via streaming instability is taken into account, partly based on previous results on streaming instability outside the water snowline that were presented in Schoonenberg & Ormel (2017). We validate our model by reproducing earlier results from the literature and apply our model to steady-state viscous gas disks (with constant gas accretion rate) around stars with different masses. We also present various other models where we explore the effects of pebble accretion, the fragmentation velocity threshold, the global metallicity of the disk, and a time-dependent gas accretion rate. We find that planetesimals preferentially form in a local annulus outside the water snowline, at early times in the lifetime of the disk ($\lesssim$$10^{5} \: \rm{yr}$), when the pebble mass fluxes are high enough to trigger the streaming instability. During this first phase in the planet formation process, the snowline location hardly changes due to slow viscous evolution, and we conclude that assuming a constant gas accretion rate is justified in this first stage. The efficiency of converting the solids reservoir of the disk to planetesimals depends on the location of the water snowline. Cooler disks with a closer-in water snowline are more efficient at producing planetesimals than hotter disks where the water snowline is located further away from the star. Therefore, low-mass stars tend to form planetesimals more efficiently, but any correlation may be overshadowed by the spread in disk properties.

The Earth is strongly depleted in carbon compared to the dust in the ISM, implying efficient removal of refractory carbon before parent body formation. It has been argued that grains get rid of their carbon through oxidation and photolysis in the exposed upper disk layers. We assess the efficacy of these C-removal mechanisms accounting for the vertical and radial transport of grains. We obtain the carbon and carbon free mass budget of solids by solving two 1D advection-diffusion equations, accounting for the dust grain size distribution and radial transport. The carbon removal acts on the fraction of the grains that are in the exposed layer and requires efficient vertical transport. In models without radial transport, oxidation and photolysis can destroy most of the refractory carbon in terrestrial planet formation region. But it only reaches the observed depletion levels for extreme parameter combinations and requires that parent body formation was delayed by 1 Myr. Adding radial transport of solids prevents the depletion entirely, leaving refractory carbon equally distributed throughout the disk. It is unlikely that the observed carbon depletion can ultimately be attributed to mechanisms operating on small grains in the disk surface layers. Other mechanisms need to be studied, for example flash heating events or FU Ori outbursts in order to remove carbon quickly and deeply. However, a sustained drift barrier or strongly reduced radial grain mobility are necessary to prevent replenishment of carbon from the outer disk.

We present nested-grid, high-resolution hydrodynamic simulations of gas and particle dynamics in the vicinity of Mars- to Earth-mass planetary embryos. The simulations extend from the surface of the embryos to a few vertical disk scale heights, with \rev{a spatial} dynamic range \rev{of} $\sim\! 1.4\times 10^5$. Our results confirm that "pebble"-sized particles are readily accreted, with accretion rates continuing to increase up to metre-size "boulders" for a 10\% MMSN surface density model. The gas mass flux in and out of the Hill sphere is consistent with the Hill rate, $\Sigma\Omega R_\mathrm{H}^2 = 4\, 10^{-3}$ M$_\oplus$ yr$^{-1}$. While smaller size particles mainly track the gas, a net accretion rate of $\approx 2\,10^{-5}$ M$_\oplus$ yr$^{-1}$ is reached for 0.3--1 cm particles, even though a significant fraction leaves the Hill sphere again. Effectively all pebble-sized particles that cross the Bondi sphere are accreted. The resolution of these simulations is sufficient to resolve accretion-driven convection. Convection driven by a nominal accretion rate of $10^{-6}$ M$_\oplus$ yr$^{-1}$ does not significantly alter the pebble accretion rate. We find that, due to cancellation effects, accretion rates of pebble-sized particles are nearly independent of disk surface density. As a result, we can estimate accurate growth times for specified particle sizes. For 0.3--1 cm size particles, the growth time from a small seed is $\sim$0.15 million years for an Earth mass planet at 1 AU and $\sim$0.1 million years for a Mars mass planet at 1.5 AU.

Coagulation theory predicts that micron-sized dust grains grow into pebbles which drift inward towards the star, when they reach sizes of mm-cm. When they cross the orbit of a planet, a fraction of these drifting pebbles will be accreted. In the pebble accretion mechanism, the combined effects of the planet's gravitational attraction and gas drag greatly increase the accretion rate. We calculate the pebble accretion efficiency $\varepsilon_{2D}$ -- the probability a pebble is accreted by the planet -- in the 2D limit (pebbles reside in the midplane). In particular, we investigate the dependence of $\varepsilon_{2D}$ on the planet eccentricity and its implications for planet formation models. We conduct N-body simulations to calculate the pebble accretion efficiency in both the local frame and the global frame. With the global method we investigate the pebble accretion efficiency when the planet is on an eccentric orbit. We find that the local and the global methods generally give consistent results. The efficiency increases with the planet's eccentricity once the relative velocity between the pebble and the planet is determined by the planet's eccentric velocity. At large eccentricities, however, the relative velocity becomes too large for pebble accretion. The efficiency then drops significantly and the accretion enters the ballistic regime. We present general expressions for $\varepsilon_{2D}$. Applying the obtained formula to the formation of a secondary planet, in resonance with an already-formed giant planet, we find that the embryo grows quickly due to its larger eccentricity. The maximum $\varepsilon_{2D}$ for a planet on an eccentric orbit is several times higher than for a planet on a circular orbit, but this increase gives the planet an important headstart and boosts its following mass growth.

Turbulence plays a key role in the transport of pebble-sized particles. It also affects the ability of pebbles to be accreted by protoplanets, because it stirs pebbles out of the disk midplane. In addition, turbulence can suppress pebble accretion once the relative velocities become too large for the settling mechanism to be viable. Following Paper I, we aim to quantify these effects by calculating the pebble accretion efficiency $\varepsilon$ using three-body simulations. To model the effect of turbulence on the pebbles, we derive a stochastic equation of motion (SEOM) applicable to stratified disk configurations. In the strong coupling limit (ignoring particle inertia) the limiting form of this equation agrees with previous works. We conduct a parameter study and calculate $\varepsilon$ in 3D, varying pebble and gas (turbulence) properties and accounting for the planet inclination. We find that strong turbulence suppresses pebble accretion through turbulent diffusion, agreeing within factors of order unity with previous works. Another reduction of $\varepsilon$ occurs when the turbulent rms motions are large and the settling mechanism fails. Efficiency-wise, the outer disk regions are more affected by turbulence than the inner regions. At the location of the H$_2$O iceline, planets around low-mass stars achieve much higher efficiencies. Including the results from Paper I, we present a framework to obtain $\varepsilon$ under general circumstances.

Sub-Neptune planets are very common in our galaxy and show a large diversity in their mass-radius relation. In sub-Neptunes most of the planet mass is in the rocky part (hereafter core) which is surrounded by a modest hydrogen-helium envelope. As a result, the total initial heat content of such a planet is dominated by that of the core. Nonetheless, most studies contend that the core cooling will only have a minor effect on the radius evolution of the gaseous envelope, because the core's cooling is in sync with the envelope, i.e., most of the initial heat is released early on timescales of about 10-100 Myr. In this Letter we examine the importance of the core cooling rate for the thermal evolution of the envelope. Thus, we relax the early core cooling assumption and present a model where the core is characterized by two parameters: the initial temperature and the cooling time. We find that core cooling can significantly enhance the radius of the planet when it operates on a timescale similar to the observed age, i.e. several Gyr. Consequently, the interpretation of sub-Neptunes' mass-radius observations depends on the assumed core thermal properties and the uncertainty therein. The degeneracy of composition and core thermal properties can be reduced by obtaining better estimates of the planet ages (in addition to their radii and masses) as envisioned by future observations.

Recently, an Atacama Large Millimeter/submillimeter Array (ALMA) observation of the water snow line in the protoplanetary disk around the FU Orionis star V883 Ori was reported. The radial variation of the spectral index at mm-wavelengths around the snow line was interpreted as being due to a pileup of particles interior to the snow line. However, radial transport of solids in the outer disk operates on timescales much longer than the typical timescale of an FU Ori outburst ($10^{1}$--$10^{2}$ yr). Consequently, a steady-state pileup is unlikely. We argue that it is only necessary to consider water evaporation and re-coagulation of silicates to explain the recent ALMA observation of V883 Ori because these processes are short enough to have had their impact since the outburst. Our model requires the inner disk to have already been optically thick before the outburst, and our results suggest that the carbon content of pebbles is low.

Context.The Kepler mission has provided a large sample to statistically analyze the orbital properties of the super-Earth planets. We hypothesize that these planets formed early and consider the problem of matching planet formation theory to the current observations. Two scenarios, disk migration and in-situ formation, have been proposed to explain their origin. In the migration scenario planets migrate inward due to planet-disk interaction, whereas in the in-situ scenario planets assemble locally. Therefore, planets formed by migration are expected to end up in resonances, whereas those formed in-situ are expected to stay in short period ratios and in non-resonant orbits. Both predictions are at odds with observations. Aims. We investigate whether a preferred formation scenario can be identified through a comparison between the magnetospheric rebound model and the Kepler data. Methods. We conduct N-body simulations of two-planet systems during the disk dispersal phase, and make a statistical comparison between the simulations and the Kepler observations. Results. Comparing the two scenarios, we find that magnetospheric rebound tends to erase the difference in the orbital configuration that was initially presented. After disk dispersal, not all planets are in resonance in the migration scenario, whereas planets do not remain in compact configurations in the in-situ scenario. In both scenarios, the orbits of planets increase with the cavity expansion, and their period ratios have a wider distribution. Conclusions. From a statistical perspective, the magnetospheric rebound model reproduces several observed properties of Kepler planets, such as the significant number of planets are not in resonances and planet pairs can end up at large period ratios. The disparity in orbital configuration between the two formation scenarios is substantially reduced after disk dispersal.

The population of close-in super-Earths, with gas mass fractions of up to 10% represents a challenge for planet formation theory: how did they avoid runaway gas accretion and collapsing to hot Jupiters despite their core masses being in the critical range of $M_\mathrm{c} \simeq 10 M_\mathrm{\oplus}$? Previous three-dimensional (3D) hydrodynamical simulations indicate that atmospheres of low-mass planets cannot be considered isolated from the protoplanetary disc, contrary to what is assumed in 1D-evolutionary calculations. This finding is referred to as the recycling hypothesis. In this Paper we investigate the recycling hypothesis for super-Earth planets, accounting for realistic 3D radiation hydrodynamics. Also, we conduct a direct comparison in terms of the evolution of the entropy between 1D and 3D geometries. We clearly see that 3D atmospheres maintain higher entropy: although gas in the atmosphere loses entropy through radiative cooling, the advection of high entropy gas from the disc into the Bondi/Hill sphere slows down Kelvin-Helmholtz contraction, potentially arresting envelope growth at a sub-critical gas mass fraction. Recycling, therefore, operates vigorously, in line with results by previous studies. However, we also identify an "inner core" -- in size $\approx$ 25% of the Bondi radius -- where streamlines are more circular and entropies are much lower than in the outer atmosphere. Future studies at higher resolutions are needed to assess whether this region can become hydrodynamically-isolated on long time-scales.

TRAPPIST-1 is a nearby 0.08 M M-star, which was recently found to harbor a planetary system of at least seven Earth-mass planets, all within 0.1 au. The configuration confounds theorists as the planets are not easily explained by either in situ or migration models. In this Paper we present a scenario for the formation and orbital architecture of the TRAPPIST-1 system. In our model, planet formation starts at the H2O iceline, where pebble-size particles -- whose origin is the outer disk -- concentrate to trigger streaming instabilities. After their formation, planetary embryos quickly mature by pebble accretion. Planet growth stalls at Earth masses, where the planet's gravitational feedback on the disk keeps pebbles at bay. Planets are transported by Type I migration to the inner disk, where they stall at the magnetospheric cavity and end up in mean motion resonances. During disk dispersal, the cavity radius expands and the inner-most planets escape resonance. We argue that the model outlined here can also be applied to other compact systems and that the many close-in super-Earth systems are a scaled-up version of TRAPPIST-1. We also hypothesize that few close-in compact systems harbor giant planets at large distances, since they would have stopped the pebble flux from the outer disk.

The formation of planetesimals in protoplanetary disks is not well-understood. Streaming instability is a promising mechanism to directly form planetesimals from pebble-sized particles, provided a high enough solids-to-gas ratio. However, local enhancements of the solids-to-gas ratio are difficult to realize in a smooth disk, which motivates the consideration of special disk locations such as the snowline. In this article we investigate the viability of planetesimal formation by streaming instability near the snowline due to water diffusion and condensation. We adopt a viscous disk model, and numerically solve the transport equations for vapor and solids on a cylindrical, 1D grid. We take into account radial drift of solids, gas accretion on to the central star, and turbulent diffusion. We study the importance of the back-reaction of solids on the gas and of the radial variation of the mean molecular weight of the gas. We find that water diffusion and condensation can locally enhance the ice surface density by a factor 3--5 outside the snowline. Assuming that icy pebbles contain many micron-sized silicate grains that are released during evaporation, the enhancement is increased by another factor $\sim$$2$. In this `many-seeds' model, the solids-to-gas ratio interior to the snowline is enhanced as well, but not as much as just outside the snowline. In the context of a viscous disk, the diffusion-condensation mechanism is most effective for high values of the turbulence parameter $\alpha$ ($10^{-3}$--$10^{-2}$). Therefore, assuming young disks are more vigorously turbulent than older disks, planetesimals near the snowline can form in an early stage of the disk. In highly turbulent disks, tens of Earth masses can be stored in an annulus outside the snowline, which can be identified with recent ALMA observations.

The Kepler mission has discovered that multiple close-in super-Earth planets are common around solar-type stars, but their period ratios do not show strong pile-ups near mean motion resonances (MMRs). One scenario is that super-Earths form in a gas-rich disk, and they interact gravitationally with the surrounding gas, inducing their orbital migration. Disk migration theory predicts, however, that planets would end up at resonant orbits due to their differential migration speed. Motivated by the discrepancy between observation and theory, we seek for a mechanism that moves planets out of resonances. We examine the orbital evolution of planet pairs near the magnetospheric cavity during the gas disk dispersal phase. Our study determines the conditions under which planets can escape resonances. We perform two-planet N-body simulations, varying the planet masses, stellar magnetic field strengths, disk accretion rates and gas disk depletion timescales. As planets migrate outward with the expanding magnetospheric cavity, their dynamical configurations can be rearranged. Migration of planets is substantial (minor) in a massive (light) disk. When the outer planet is more massive than the inner planet, the period ratio of two planets increases through outward migration. On the other hand, when the inner planet is more massive, the final period ratio tends to remain similar to the initial one. Larger stellar magnetic field strengths result in planets stopping their migration at longer periods. We highlight \textit{magnetospheric rebound} as an important ingredient able to reconcile disk migration theory with observations. Even when planets are trapped into MMR during the early gas-rich stage, subsequent cavity expansion would induce substantial changes to their orbits, moving them out of resonance.

Planetesimals in gaseous protoplanetary disks may grow by collecting dust particles. Hydrodynamical studies show that small particles generally avoid collisions with the planetesimals because they are entrained by the flow around them. This occurs when $St$, the Stokes number, defined as the ratio of the dust stopping time to the planetesimal crossing time, becomes much smaller than unity. However, these studies have been limited to the laminar case, whereas these disks are believed to be turbulent. We want to estimate the influence of gas turbulence on the dust-planetesimal collision rate and on the impact speeds. We used three-dimensional direct numerical simulations of a fixed sphere (planetesimal) facing a laminar and turbulent flow seeded with small inertial particles (dust) subject to a Stokes drag. A no-slip boundary condition on the planetesimal surface is modeled via a penalty method. We find that turbulence can significantly increase the collision rate of dust particles with planetesimals. For a high turbulence case (when the amplitude of turbulent fluctuations is similar to the headwind velocity), we find that the collision probability remains equal to the geometrical rate or even higher for $St\geq 0.1$, i.e., for dust sizes an order of magnitude smaller than in the laminar case. We derive expressions to calculate impact probabilities as a function of dust and planetesimal size and turbulent intensity.

The journey from dust particle to planetesimal involves physical processes acting on scales ranging from micrometers (the sticking and restructuring of aggregates) to hundreds of astronomical units (the size of the turbulent protoplanetary nebula). Considering these processes simultaneously is essential when studying planetesimal formation. We develop a novel, global, semi-analytical model for the evolution of the mass-dominating dust particles in a turbulent protoplanetary disk that takes into account the evolution of the dust surface density while preserving the essential characteristics of the porous coagulation process. This panoptic model is used to study the growth from submicron to planetesimal sizes in disks around Sun-like stars. For highly porous ices, unaffected by collisional fragmentation and erosion, rapid growth to planetesimal sizes is possible in a zone stretching out to ${\sim}10\mathrm{~AU}$ for massive disks. When porous coagulation is limited by erosive collisions, the formation of planetesimals through direct coagulation is not possible, but the creation of a large population of aggregates with Stokes numbers close to unity might trigger the streaming instability (SI). However, we find that reaching conditions necessary for SI is difficult and limited to dust-rich disks, (very) cold disks, or disks with weak turbulence. Behind the snow-line, porosity-driven aggregation of icy grains results in rapid (${\sim}10^{4}\mathrm{~yr}$) formation of planetesimals. If erosive collisions prevent this, SI might be triggered for specific disk conditions. The numerical approach introduced in this work is ideally suited for studying planetesimal formation and pebble delivery simultaneously and will help build a coherent picture of the start of the planet formation process.

Pebble accretion is a new mechanism to quickly grow the cores of planets. In pebble accretion, gravity and gas drag conspire to yield large collisional cross sections for small particles in protoplanetary disks. However, before pebble accretion commences, aerodynamical deflection may act to prevent planetesimals from becoming large, because particles tend to follow gas streamlines. We derive the planetesimal radius where pebble accretion is initiated and determine the growth timescales of planetesimals by sweepup of small particles. We obtain the collision efficiency factor as the ratio of the numerically-obtained collisional cross section to the planetesimal surface area, from which we obtain the growth timescales. Integrations are conducted in the potential flow limit (steady, inviscid) and in the Stokes flow regime (steady, viscid). Only particles of stopping time $t_s \ll t_X$ where $t_X\approx10^3$ s experience aerodynamic deflection. Even in that case, the planetesimal's gravity always ensures positive collision factors. The maximum growth timescale occurs typically at around $R\approx100 \ \mathrm{km}$, but is less for colder disks, corresponding to interactions shifting to the Safronov focusing regime. For particles $t_s \gg t_X$ pebble accretion commences only after this phase and is characterized by a steep drop in growth timescales. Consequently, at distances beyond ~10 AU sweepup growth timescales are always longer than $10$ Myr, while in the inner disk (~<3 AU) the viability of the sweepup scenario is determined by the outcome of pebble-planetesimal collisions in the geometric regime. We present analytical fits for the collision efficiency factors and the minimum planetesimal size needed for pebble accretion. (Abridged)

The coagulation of microscopic dust into planetesimals is the first step towards planet formation. The size and shape of the growing aggregates determine the efficiency of this early growth. It has been proposed that fluffy ice aggregates can grow very efficiently, suffering less from the bouncing and radial drift barriers. While the collision velocity between icy aggregates of similar size is thought to stay below the fragmentation threshold, they may nonetheless lose mass from collisions with much smaller projectiles. We investigate the effect of these erosive collisions on the ability of porous ice aggregates to cross the radial drift barrier. We develop a Monte Carlo code that calculates the evolution of the growing aggregates, while resolving the entire mass distribution at all times. The aggregate's porosity is treated independently of its mass, and is determined by collisions, gas compaction, and eventually self-gravity compaction. For erosion threshold velocities of 20-40 m/s, high-velocity collisions with small projectiles prevent the largest aggregates from growing when they start to drift. In these cases, our local simulations result in a steady-state distribution, with the majority of the dust mass in particles with Stokes numbers close to unity. Only for the highest erosion threshold considered (60 m/s), do porous aggregates manage to cross the radial drift barrier in the inner 10 AU of MMSN-like disks. Erosive collisions are more effective in limiting the growth than fragmentary collisions between similar-size particles. Conceivably, erosion limits the growth before the radial drift barrier, although the robustness of this statement depends on (uncertain) material properties of icy aggregates. If erosion inhibits planetesimal formation through direct sticking, the sea of ${\sim}10^9$ g, highly porous particles appears well-suited for triggering streaming instability.

Accumulation of dust and ice particles into planetesimals is an important step in the planet formation process. Planetesimals are the seeds of both terrestrial planets and the solid cores of gas and ice giants forming by core accretion. Left-over planetesimals in the form of asteroids, trans-Neptunian objects and comets provide a unique record of the physical conditions in the solar nebula. Debris from planetesimal collisions around other stars signposts that the planetesimal formation process, and hence planet formation, is ubiquitous in the Galaxy. The planetesimal formation stage extends from micrometer-sized dust and ice to bodies which can undergo run-away accretion. The latter ranges in size from 1 km to 1000 km, dependent on the planetesimal eccentricity excited by turbulent gas density fluctuations. Particles face many barriers during this growth, arising mainly from inefficient sticking, fragmentation and radial drift. Two promising growth pathways are mass transfer, where small aggregates transfer up to 50% of their mass in high-speed collisions with much larger targets, and fluffy growth, where aggregate cross sections and sticking probabilities are enhanced by a low internal density. A wide range of particle sizes, from mm to 10 m, concentrate in the turbulent gas flow. Overdense filaments fragment gravitationally into bound particle clumps, with most mass entering planetesimals of contracted radii from 100 to 500 km, depending on local disc properties. We propose a hybrid model for planetesimal formation where particle growth starts unaided by self-gravity but later proceeds inside gravitationally collapsing pebble clumps to form planetesimals with a wide range of sizes.

Close-in super-Earths having radii 1--4 $R_\oplus$ may possess hydrogen atmospheres comprising a few percent by mass of their rocky cores. We determine the conditions under which such atmospheres can be accreted by cores from their parent circumstellar disks. Accretion from the nebula is problematic because it is too efficient: we find that 10-$M_\oplus$ cores embedded in solar metallicity disks tend to undergo runaway gas accretion and explode into Jupiters, irrespective of orbital location. The threat of runaway is especially dire at $\sim$0.1 AU, where solids may coagulate on timescales orders of magnitude shorter than gas clearing times; thus nascent atmospheres on close-in orbits are unlikely to be supported against collapse by planetesimal accretion. The time to runaway accretion is well approximated by the cooling time of the atmosphere's innermost convective zone, whose extent is controlled by where H$_2$ dissociates. Insofar as the temperatures characterizing H$_2$ dissociation are universal, timescales for core instability tend not to vary with orbital distance --- and to be alarmingly short for 10-$M_\oplus$ cores. Nevertheless, in the thicket of parameter space, we identify two scenarios, not mutually exclusive, that can reproduce the preponderance of percent-by-mass atmospheres for super-Earths at $\sim$0.1 AU, while still ensuring the formation of Jupiters at $\gtrsim 1$ AU. Scenario (a): planets form in disks with dust-to-gas ratios that range from $\sim$20$\times$ solar at 0.1 AU to $\sim$2$\times$ solar at 5 AU. Scenario (b): the final assembly of super-Earth cores from mergers of proto-cores --- a process that completes quickly at $\sim$0.1 AU once begun --- is delayed by gas dynamical friction until just before disk gas dissipates completely.

In the core accretion paradigm of planet formation, gas giants only form a massive atmosphere after their progenitors exceeded a threshold mass: the critical core mass. Most (exo)planets, being smaller and rock/ice-dominated, never crossed this line. Nevertheless, they were massive enough to attract substantial amounts of gas from the disc, while their atmospheres remained in pressure-equilibrium with the disc. Our goal is to characterise the hydrodynamical properties of the atmospheres of such embedded planets and their implication for their (long-term) evolution. In this paper -- the first in series -- we start to investigate the properties of an isothermal and inviscid flow past a small, embedded planet by conducting local, 2D hydrodynamical simulations. Using the PLUTO code we confirm that the flow is steady and bound. This steady outcome is most apparent for the log-polar grid (with the grid spacing proportional to the distance from the planet). For low-mass planets, Cartesian grids are somewhat less efficient as they have difficulty to follow the circular, large speeds in the deep atmosphere. Relating the amount of rotation to the gas fraction of the atmosphere, we find that more massive atmospheres rotate faster -- a finding consistent with Kelvin's circulation theorem. Rotation therefore limits the amount of gas that planets can acquire from the nebula. Dependent on the Toomre-Q parameter of the circumstellar disc, the planet's atmosphere will reach Keplerian rotation before self-gravity starts to become important.

Following Paper I we investigate the properties of atmospheres that form around small protoplanets embedded in a protoplanetary disc by conducting hydrodynamical simulations. These are now extended to three dimensions, employing a spherical grid centred on the planet. Compression of gas is shown to reduce rotational motions. Contrasting the 2D case, no clear boundary demarcates bound atmospheric gas from disc material; instead, we find an open system where gas enters the Bondi sphere at high latitudes and leaves through the midplane regions, or, vice versa, when the disc gas rotates sub-Keplerian. The simulations do not converge to a time-independent solution; instead, the atmosphere is characterized by a time-varying velocity field. Of particular interest is the timescale to replenish the atmosphere by nebular gas, $t_\mathrm{replenish}$. It is shown that the replenishment rate, $M_\mathrm{atm}/t_\mathrm{replenish}$, can be understood in terms of a modified Bondi accretion rate, $\sim$$R_\mathrm{Bondi}^2\rho_\mathrm{gas}v_\mathrm{Bondi}$, where $v_\mathrm{Bondi}$ is set by the Keplerian shear or the magnitude of the sub-Keplerian motion of the gas, whichever is larger. In the inner disk, the atmosphere of embedded protoplanets replenishes on a timescale that is shorter than the Kelvin-Helmholtz contraction (or cooling) timescale. As a result, atmospheric gas can no longer contract and the growth of these atmospheres terminates. Future work must confirm whether these findings continue to apply when the (thermodynamical) idealizations employed in this study are relaxed. But if shown to be broadly applicable, replenishment of atmospheric gas provides a natural explanation for the preponderance of gas-rich but rock-dominant planets like super-Earths and mini-Neptunes.

Context. Circumstellar disks are known to contain a significant mass in dust ranging from micron to centimeter size. Meteorites are evidence that individual grains of those sizes were collected and assembled into planetesimals in the young solar system. Aims. We assess the efficiency of dust collection of a swarm of non-drifting planetesimals {\rev with radii ranging from 1 to $10^3$\,km and beyond. Methods. We calculate the collision probability of dust drifting in the disk due to gas drag by planetesimal accounting for several regimes depending on the size of the planetesimal, dust, and orbital distance: the geometric, Safronov, settling, and three-body regimes. We also include a hydrodynamical regime to account for the fact that small grains tend to be carried by the gas flow around planetesimals. Results. We provide expressions for the collision probability of dust by planetesimals and for the filtering efficiency by a swarm of planetesimals. For standard turbulence conditions (i.e., a turbulence parameter $\alpha=10^{-2}$), filtering is found to be inefficient, meaning that when crossing a minimum-mass solar nebula (MMSN) belt of planetesimals extending between 0.1 AU and 35 AU most dust particles are eventually accreted by the central star rather than colliding with planetesimals. However, if the disk is weakly turbulent ($\alpha=10^{-4}$) filtering becomes efficient in two regimes: (i) when planetesimals are all smaller than about 10 km in size, in which case collisions mostly take place in the geometric regime; and (ii) when planetary embryos larger than about 1000 km in size dominate the distribution, have a scale height smaller than one tenth of the gas scale height, and dust is of millimeter size or larger in which case most collisions take place in the settling regime. These two regimes have very different properties: we find that the local filtering efficiency $x_{filter,MMSN}$ scales with $r^{-7/4}$ (where $r$ is the orbital distance) in the geometric regime, but with $r^{-1/4}$ to $r^{1/4}$ in the settling regime. This implies that the filtering of dust by small planetesimals should occur close to the central star and with a short spread in orbital distances. On the other hand, the filtering by embryos in the settling regime is expected to be more gradual and determined by the extent of the disk of embryos. Dust particles much smaller than millimeter size tend only to be captured by the smallest planetesimals because they otherwise move on gas streamlines and their collisions take place in the hydrodynamical regime. Conclusions. Our results hint at an inside-out formation of planetesimals in the infant solar system because small planetesimals in the geometrical limit can filter dust much more efficiently close to the central star. However, even a fully-formed belt of planetesimals such as the MMSN only marginally captures inward-drifting dust and this seems to imply that dust in the protosolar disk has been filtered by planetesimals even smaller than 1 km (not included in this study) or that it has been assembled into planetesimals by other mechanisms (e.g., orderly growth, capture into vortexes). Further refinement of our work concerns, among other things: a quantitative description of the transition region between the hydro and settling regimes; an assessment of the role of disk turbulence for collisions, in particular in the hydro regime; and the coupling of our model to a planetesimal formation model.

We present a method to include the evolution of the grain size and grain opacity $\kappa_\mathrm{gr}$ in the equations describing the structure of protoplanetary atmospheres. The key assumption of this method is that a single grain size dominates the grain size distribution at any height $r$. In addition to following grain growth, the method accounts for mass deposition by planetesimals and grain porosity. We illustrate this method by computation of a simplified atmosphere structure model. In agreement with previous works, grain coagulation is seen to be very efficient. The opacity drops to values much below the often-used `ISM-opacities' ($\sim$$1\ \mathrm{cm^2\ g}^{-1}$) and the atmosphere structure profiles for temperature and density resemble that of the grain-free case. Deposition of planetesimals in the radiative part of the atmosphere hardly influences this outcome as the added surface is quickly coagulated away. We observe a modest dependence on the internal structure (porosity), but show that filling factors cannot become too large because of compression by gas drag.

High levels of exozodiacal dust have been observed in the inner regions of a large fraction of main sequence stars. Given the short lifetime of the observed small dust grains, these 'exozodis' are difficult to explain, especially for old (>100 Myr) stars. The exozodiacal dust may be observed as excess emission in the mid-infrared, or using interferometry. We hypothesise that exozodi are sustained by planetesimals scattered by planets inwards from an outer planetesimal belt, where collision timescales are long. In this work, we use N-body simulations to show that the outwards migration of a planet into a belt, driven by the scattering of planetesimals, can increase, or sustain, the rate at which planetesimals are scattered from the outer belt to the exozodi region. We hypothesise that this increase is sufficient to sustain the observed exozodi on Gyr timescales. No correlation between observations of an outer belt and an exozodi is required for this scenario to work, as the outer belt may be too faint to detect. If planetesimal driven migration does explain the observed exozodi, this work suggests that the presence of an exozodi indicates the presence of outer planets and a planetesimal belt.

The recently discovered coreshine effect can aid in exploring the core properties and in probing the large grain population of the ISM. We discuss the implications of the coreshine detected from the molecular cloud core L1506C in the Taurus filament for the history of the core and the existence of a primitive ISM component of large grains becoming visible in cores. The coreshine surface brightness of L1506C is determined from IRAC Spitzer images at 3.6 micron. We perform grain growth calculations to estimate the grain size distribution in model cores similar in gas density, radius, and turbulent velocity to L1506C. Scattered light intensities at 3.6 micron are calculated for a variety of MRN and grain growth distributions to compare with the observed coreshine. For a core with the overall physical properties of L1506C, no detectable coreshine is predicted for an MRN size distribution. Extending the distribution to grain radii of about 0.65 $\mu$m allows to reproduce the observed surface brightness level in scattered light. Assuming the properties of L1506C to be preserved, models for the growth of grains in cores do not yield sufficient scattered light to account for the coreshine within the lifetime of the Taurus complex. Only increasing the core density and the turbulence amplifies the scattered light intensity to a level consistent with the observed coreshine brightness. The grains could be part of primitive omni-present large grain population becoming visible in the densest part of the ISM, could grow under the turbulent dense conditions of former cores, or in L1506C itself. In the later case, L1506C must have passed through a period of larger density and stronger turbulence. This would be consistent with the surprisingly strong depletion usually attributed to high column densities, and with the large-scale outward motion of the core envelope observed today.

Turbulence in protoplanetary disks affects planet formation in many ways. While small dust particles are mainly affected by the aerodynamical coupling with turbulent gas velocity fields, planetesimals and larger bodies are more affected by gravitational interaction with gas density fluctuations. For the latter process, a number of numerical simulations have been performed in recent years, but a fully parameter-independent understanding has not been yet established. In this study, we present simple scaling relations for the planetesimal stirring rate in turbulence driven by magnetorotational instability (MRI), taking into account the stabilization of MRI due to Ohmic resistivity. We begin with order-of-magnitude estimates of the turbulence-induced gravitational force acting on solid bodies and associated diffusion coefficients for their orbital elements. We then test the predicted scaling relations using the results of recent Ohmic-resistive MHD simulations by Gressel et al. We find that these relations successfully explain the simulation results if we properly fix order-of-unity uncertainties within the estimates. We also update the saturation predictor for the density fluctuation amplitude in MRI-driven turbulence originally proposed by Okuzumi & Hirose. Combination of the scaling relations and saturation predictor allows to know how the turbulent stirring rate of planetesimals depends on disk parameters such as the gas column density, distance from the central star, vertical resistivity distribution, and net vertical magnetic flux. In Paper II, we apply our recipe to planetesimal accretion to discuss its viability in turbulent disks.

A critical phase in the standard model for planet formation is the runaway growth phase. During runaway growth bodies in the 0.1--100 km size range (planetesimals) quickly produce a number of much larger seeds. The runaway growth phase is essential for planet formation as the emergent planetary embryos can accrete the leftover planetesimals at large gravitational focusing factors. However, torques resulting from turbulence-induced density fluctuations may violate the criterion for the onset of runaway growth, which is that the magnitude of the planetesimals' random (eccentric) motions are less than their escape velocity. This condition represents a more stringent constraint than the condition that planetesimals survive their mutual collisions. To investigate the effects of MRI turbulence on the viability of the runaway growth scenario, we apply our semi-analytical recipes of Paper I, which we augment by a coagulation/fragmentation model for the dust component. We find that the surface area-equivalent abundance of 0.1 micron particles is reduced by factors 10^2--10^3, which tends to render the dust irrelevant to the turbulence. We express the turbulent activity in the midplane regions in terms of a size s_run above which planetesimals will experience runaway growth. We find that s_run is mainly determined by the strength of the vertical net field that threads the disks and the disk radius. At disk radii beyond 5 AU, s_run becomes larger than ~100 km and the collision times among these bodies longer than the duration of the nebula phase. Our findings imply that the classical, planetesimal-dominated, model for planet formation is not viable in the outer regions of a turbulent disk.

Gravitating bodies significantly alter the flow pattern (density and velocity) of the gas that attempts to stream past. Still, small protoplanets in the Mars--super-Earth range can only bind limited amounts of nebular gas; until the so-called critical core mass has been reached (~1-10 Earth masses) this gas is in near hydrostatic equilibrium with the nebula. Here we aim for a general description of the flow pattern surrounding these low-mass, embedded planets. Using various simplifying assumptions (subsonic, 2D, inviscid flow, etc), we reduce the problem to a partial differential equation that we solve numerically as well as approximate analytically. It is found that the boundary between the atmosphere and the nebula gas strongly depends on the value of the disc headwind (deviation from Keplerian rotation). With increasing headwind the atmosphere decreases in size and also becomes more asymmetrical. Using the derived flow pattern for the gas, trajectories of small solid particles, which experience both gas drag and gravitational forces, are integrated numerically. Accretion rates for small particles (dust) are found to be low, as they closely follow the streamlines, which curl away from the planet. However, pebble-size particles achieve large accretion rates, in agreement with previous numerical and analytical works.

It is unknown how far dust growth can proceed by coagulation. Obstacles to collisional growth are the fragmentation and bouncing barriers. However, in all previous simulations of the dust-size evolution in protoplanetary disks, only the mean collision velocity has been considered, neglecting that a small but possibly important fraction of the collisions will occur at both much lower and higher velocities. We study the effect of the probability distribution of impact velocities on the collisional dust growth barriers. Assuming a Maxwellian velocity distribution for colliding particles to determine the fraction of sticking, bouncing, and fragmentation, we implement this in a dust-size evolution code. We also calculate the probability of growing through the barriers and the growth timescale in these regimes. We find that the collisional growth barriers are not as sharp as previously thought. With the existence of low-velocity collisions, a small fraction of the particles manage to grow to masses orders of magnitude above the main population. A particle velocity distribution softens the fragmentation barrier and removes the bouncing barrier. It broadens the size distribution in a natural way, allowing the largest particles to become the first seeds that initiate sweep-up growth towards planetesimal sizes.

Planets migrate due to the recoil they experience from scattering solid (planetesimal) bodies. To first order, the torques exerted by the interior and exterior disks cancel, analogous to the cancellation of the torques from the gravitational interaction with the gas (type I migration). Assuming the dispersion-dominated regime and power-laws characterized by indices {\alpha} and {\beta} for the surface density and eccentricity profiles, we calculate the net torque on the planet. We consider both distant encounters and close (orbit-crossing) encounters. We find that the close and distant encounter torques have opposite signs with respect to their {\alpha} and {\beta} dependences; and that the torque is especially sensitive to the eccentricity gradient ({\beta}). Compared to type-I migration due to excitation of density waves, the planetesimal-driven migration rate is generally lower due to the lower surface density of solids in gas-rich disk, although this may be partially or fully offset when their eccentricity and inclination are small. Allowing for the feedback of the planet on the planetesimal disk through viscous stirring, we find that under certain conditions a self-regulated migration scenario emerges, in which the planet migrates at a steady pace that approaches the rate corresponding to the one-sided torque. If the local planetesimal disk mass to planet mass ratio is low, however, migration stalls. We quantify the boundaries separating the three migration regimes.

We have investigated Saturn's core formation at a radial pressure maximum in a protoplanetary disk, which is created by gap opening by Jupiter. A core formed via planetesimal accretion induces the fragmentation of surrounding planetesimals, which generally inhibits further growth of the core by removal of the resulting fragments due to radial drift caused by gas drag. However, the emergence of the pressure maximum halts the drift of the fragments, while their orbital eccentricities and inclinations are efficiently damped by gas drag. As a result, the core of Saturn rapidly grows via accretion of the fragments near the pressure maximum. We have found that in the minimum-mass solar nebula, kilometer sized planetesimals can produce a core exceeding 10 Earth masses within two million years. Since Jupiter may not have undergone significant type II inward migration, it is likely that Jupiter's formation was completed when the local disk mass has already decayed to a value comparable to or less than Jovian mass. The expected rapid growth of Saturn's core on a timescale comparable to or shorter than observationally inferred disk lifetime enables Saturn to acquire the current amount of envelope gas before the disk gas is completely depleted. The high heat energy release rate onto the core surface due to the rapid accretion of the fragments delays onset of runaway gas accretion until the core mass becomes somewhat larger than that of Jupiter, which is consistent with the estimate based on interior modeling. Therefore, the rapid formation of Saturn induced by gap opening of Jupiter can account for the formation of multiple gas giants (Jupiter and Saturn) without significant inward migration and larger core mass of Saturn than that of Jupiter.

The formation of giant planets requires accumulation of ~10 Earth mass in solids; but how do protoplanets acquire their mass? There are many, often competing processes that regulate the accretion rate of protoplanets. To assess their effects we present a new, publicly-available toy model. The rationale behind the toy model is that it encompasses as many physically-relevant processes as possible, but at the same time does not compromise its simplicity, speed, and physical insight. The toy model follows a modular structure, where key features -- e.g. planetesimal fragmentation, radial orbital decay, nebula turbulence -- can be switched on or off. Our model assumes three discrete components (fragments, planetesimals, and embryos) and is zero dimensional in space. We have tested the outcomes of the toy model against literature results and generally find satisfactory agreement. We include, for the first time, model features that capture the three-way interactions among small particles, gas, and protoplanets. Collisions among planetesimals will result in fragmentation, transferring a substantial amount of the solid mass to small particles, which couple strongly to the gas. Our results indicate that the efficiency of the accretion process then becomes very sensitive to the gas properties -- especially to the turbulent state and the magnitude of the disk headwind (the decrease of the orbital velocity of the gas with respect to Keplerian) -- as well as to the characteristic fragment size.

We investigate dust growth due to settling in a 1D vertical column of a protoplanetary disk. It is known from the observed 10 micron feature in disk SEDs, that small micron-sized grains are present at the disk atmosphere throughout the lifetime of the disk. We hope to explain such questions as what process can keep the disk atmospheres dusty for the lifetime of the disk and how does the particle properties change as a function of height above the midplane. We use a Monte Carlo code to follow the mass and porosity evolution of the particles in time. The used collision model is based on laboratory experiments performed on dust aggregates. As the experiments cannot cover all possible collision scenarios, the largest uncertainty of our model is the necessary extrapolations we had to perform. We simultaneously solve for the particle growth and motion. Particles can move vertically due to settling and turbulent mixing. We assume that the vertical profile of the gas density is fixed in time and only the solid component evolves. We find that the used collision model strongly influences the masses and sizes of the particles. The laboratory experiment based collision model greatly reduces the particle sizes compared to models that assume sticking at all collision velocities. We find that a turbulence parameter of alpha = 10^-2 is needed to keep the dust atmospheres dusty, but such strong turbulence can produce only small particles at the midplane which is not favorable for planetesimal formation models. We also see that the particles are larger at the midplane and smaller at the upper layers of the disk. At 3-4 pressure scale heights micron-sized particles are produced. These particle sizes are needed to explain the 10 micron feature of disk SEDs. Turbulence may therefore help to keep small dust particles in the disk atmosphere.

The dust size distribution in molecular clouds can be strongly affected by ice-mantle formation and (subsequent) grain coagulation. Following previous work where the dust size distribution has been calculated from a state-of-the art collision model for dust aggregates that involves both coagulation and fragmentation (Paper I), the corresponding opacities are presented in this study. The opacities are calculated by applying the effective medium theory assuming that the dust aggregates are a mix of 0.1{\mu}m silicate and graphite grains and vacuum. In particular, we explore how the coagulation affects the near-IR opacities and the opacity in the 9.7{\mu}m silicate feature. We find that as dust aggregates grow to {\mu}m-sizes both the near-IR color excess and the opacity in the 9.7 {\mu}m feature increases. Despite their coagulation, porous aggregates help to prolong the presence of the 9.7{\mu}m feature. We find that the ratio between the opacity in the silicate feature and the near-IR color excess becomes lower with respect to the ISM, in accordance with many observations of dark clouds. However, this trend is primarily a result of ice mantle formation and the mixed material composition of the aggregates, rather than being driven by coagulation. With stronger growth, when most of the dust mass resides in particles of size 10{\mu}m or larger, both the near-IR color excess and the 9.7{\mu}m silicate feature significantly diminish. Observations at additional wavelengths, in particular in the sub-mm range, are essential to provide quantitative constraints on the dust size distribution within dense cores. Our results indicate that the sub-mm index {\beta} will increase appreciably, if aggregates grow to ~100{\mu}m in size.

We report the first spectroscopic mapping of an atomic carbon line in an infrared dark cloud (IRDC). By observing the spatial distribution of the [CI] emission in an IRDC, comparing it with the 13CO emission and the known distribution of internal heating sources, we can quantify the role of internal and external UV irradiation in the production of atomic carbon. We used the 2x4 pixel SMART receiver of the KOSMA observatory on Gornergrat to map the [CI] 3P1 - 3P0 line in the IRDC G48.65-0.29 and compared the resulting spectra with data from the BU-FCRAO 13CO 1-0 Galactic Ring Survey. The [CI]/13CO effective beam temperature ratio falls at about 0.3 with local deviations by less than a factor two. All velocity components seen in 13CO are also detected in[CI]. We find, however, significant differences in the morphology of the brightest regions seen in the two tracers. While 13CO basically follows the column density distribution derived from the near-infrared (NIR) extinction and the submm continuum, the [CI] emission peaks at the locations of the two known NIR point sources. We find CI/CO abundance ratios between 0.07 and 0.13, matching the lower end of the range previously measured in star-forming regions. Evaluating the relative importance of the irradiation by embedded sources and by the Galactic interstellar UV field, we find that in G48.65-0.29 most [CI] emission can be attributed to externally illuminated surfaces. Embedded sources have a significant impact on the overall abundance distribution of atomic carbon as soon as they are found in an evolved state with noticeable NIR flux.

Context. Grains in circumstellar disks are believed to grow by mutual collisions and subsequent sticking due to surface forces. Results of many fields of research involving circumstellar disks, such as radiative transfer calculations, disk chemistry, magneto-hydrodynamic simulations largely depend on the unknown grain size distribution. Aims. As detailed calculations of grain growth and fragmentation are both numerically challenging and computationally expensive, we aim to find simple recipes and analytical solutions for the grain size distribution in circumstellar disks for a scenario in which grain growth is limited by fragmentation and radial drift can be neglected. Methods. We generalize previous analytical work on self-similar steady-state grain distributions. Numerical simulations are carried out to identify under which conditions the grain size distributions can be understood in terms of a combination of power-law distributions. A physically motivated fitting formula for grain size distributions is derived using our analytical predictions and numerical simulations. Results. We find good agreement between analytical results and numerical solutions of the Smoluchowski equation for simple shapes of the kernel function. The results for more complicated and realistic cases can be fitted with a physically motivated "black box" recipe presented in this paper. Our results show that the shape of the dust distribution is mostly dominated by the gas surface density (not the dust-to-gas ratio), the turbulence strength and the temperature and does not obey an MRN type distribution.

Planetary bodies form by accretion of smaller bodies. It has been suggested that a very efficient way to grow protoplanets is by accreting particles of size <<km (e.g., chondrules, boulders, or fragments of larger bodies) as they can be kept dynamically cold. We investigate the effects of gas drag on the impact radii and the accretion rates of these particles. As simplifying assumptions we restrict our analysis to 2D settings, a gas drag law linear in velocity, and a laminar disk characterized by a smooth (global) pressure gradient that causes particles to drift in radially. These approximations, however, enable us to cover an arbitrary large parameter space. The framework of the circularly restricted three body problem is used to numerically integrate particle trajectories and to derive their impact parameters. Three accretion modes can be distinguished: hyperbolic encounters, where the 2-body gravitational focusing enhances the impact parameter; three-body encounters, where gas drag enhances the capture probability; and settling encounters, where particles settle towards the protoplanet. An analysis of the observed behavior is presented; and we provide a recipe to analytically calculate the impact radius, which confirms the numerical findings. We apply our results to the sweepup of fragments by a protoplanet at a distance of 5 AU. Accretion of debris on small protoplanets (<50 km) is found to be slow, because the fragments are distributed over a rather thick layer. However, the newly found settling mechanism, which is characterized by much larger impact radii, becomes relevant for protoplanets of ~10^3 km in size and provides a much faster channel for growth.

(abridged) When preplanetary bodies reach proportions of ~1 km or larger in size, their accretion rate is enhanced due to gravitational focusing (GF). We have developed a new numerical model to calculate the collisional evolution of the gravitationally-enhanced growth stage. We validate our approach against existing N-body and statistical codes. Using the numerical model, we explore the characteristics of the runaway growth and the oligarchic growth accretion phases starting from an initial population of single planetesimal radius R_0. In models where the initial random velocity dispersion (as derived from their eccentricity) starts out below the escape speed of the planetesimal bodies, the system experiences runaway growth. We find that during the runaway growth phase the size distribution remains continuous but evolves into a power-law at the high mass end, consistent with previous studies. Furthermore, we find that the largest body accretes from all mass bins; a simple two component approximation is inapplicable during this stage. However, with growth the runaway body stirs up the random motions of the planetesimal population from which it is accreting. Ultimately, this feedback stops the fast growth and the system passes into oligarchy, where competitor bodies from neighboring zones catch up in terms of mass. Compared to previous estimates, we find that the system leaves the runaway growth phase at a somewhat larger radius. Furthermore, we assess the relevance of small, single-size fragments on the growth process. In classical models, where the initial velocity dispersion of bodies is small, these do not play a critical role during the runaway growth; however, in models that are characterized by large initial relative velocities due to external stirring of their random motions, a situation can emerge where fragments dominate the accretion.

Context. Observations at sub-millimeter and mm wavelengths will in the near future be able to resolve the radial dependence of the mm spectral slope in circumstellar disks with a resolution of around a few AU at the distance of the closest star-forming regions. Aims. We aim to constrain physical models of grain growth and fragmentation by a large sample of (sub-)mm observations of disks around pre-main sequence stars in the Taurus-Auriga and Ophiuchus star-forming regions. Methods. State-of-the-art coagulation/fragmentation and disk-structure codes are coupled to produce steady-state grain size distributions and to predict the spectral slopes at (sub-)mm wavelengths. Results. This work presents the first calculations predicting the mm spectral slope based on a physical model of grain growth. Our models can quite naturally reproduce the observed mm-slopes, but a simultaneous match to the observed range of flux levels can only be reached by a reduction of the dust mass by a factor of a few up to about 30 while keeping the gas mass of the disk the same. This dust reduction can either be due to radial drift at a reduced rate or during an earlier evolutionary time (otherwise the predicted fluxes would become too low) or due to efficient conversion of dust into larger, unseen bodies.

Accretion among macroscopic bodies of ~km size or larger is enhanced significantly due to gravitational focusing. Two regimes can be distinguished. Initially, the system experiences runaway growth, in which the gravitational focusing factors increase, and bodies at the high-mass tail of the distribution grow fastest. However, at some point the runaway body dynamically heats its environment, gravitational focusing factors decrease, and runaway growth passes into oligarchic growth. Based on the results of recent simulations, we reconsider the runaway growth-oligarchy transition. In contrast to oligarchy, we find that runaway growth cannot be approximated with a two component model (of small and large bodies) and that the criterion of Ida & Makino (1993), which is frequently adopted as the start of oligarchy, is not a sufficient condition to signify the transition. Instead, we propose a new criterion based on timescale arguments. We then find a larger value for the runaway growth-oligarchy transition: from several hundreds of km in the inner disk regions up to ~10^3 km for the outer disk. These findings are consistent with the view that runaway growth has been responsible for the size distribution of the present day Kuiper belt objects. Our finding furthermore outlines the proper initial conditions at the start of the oligarchy stage.

The sticking of micron sized dust particles due to surface forces in circumstellar disks is the first stage in the production of asteroids and planets. The key ingredients that drive this process are the relative velocity between the dust particles in this environment and the complex physics of dust aggregate collisions. Here we present the results of a collision model, which is based on laboratory experiments of these aggregates. We investigate the maximum aggregate size and mass that can be reached by coagulation in protoplanetary disks. We model the growth of dust aggregates at 1 AU at the midplane at three different gas densities. We find that the evolution of the dust does not follow the previously assumed growth-fragmentation cycles. Catastrophic fragmentation hardly occurs in the three disk models. Furthermore we see long lived, quasi-steady states in the distribution function of the aggregates due to bouncing. We explore how the mass and the porosity change upon varying the turbulence parameter and by varying the critical mass ratio of dust particles. Particles reach Stokes numbers of roughly 10^-4 during the simulations. The particle growth is stopped by bouncing rather than fragmentation in these models. The final Stokes number of the aggregates is rather insensitive to the variations of the gas density and the strength of turbulence. The maximum mass of the particles is limited to approximately 1 gram (chondrule-sized particles). Planetesimal formation can proceed via the turbulent concentration of these aerodynamically size-sorted chondrule-sized particles.

The growth processes from protoplanetary dust to planetesimals are not fully understood. Laboratory experiments and theoretical models have shown that collisions among the dust aggregates can lead to sticking, bouncing, and fragmentation. However, no systematic study on the collisional outcome of protoplanetary dust has been performed so far so that a physical model of the dust evolution in protoplanetary disks is still missing. We intend to map the parameter space for the collisional interaction of arbitrarily porous dust aggregates. This parameter space encompasses the dust-aggregate masses, their porosities and the collision velocity. With such a complete mapping of the collisional outcomes of protoplanetary dust aggregates, it will be possible to follow the collisional evolution of dust in a protoplanetary disk environment. We use literature data, perform own laboratory experiments, and apply simple physical models to get a complete picture of the collisional interaction of protoplanetary dust aggregates. In our study, we found four different types of sticking, two types of bouncing, and three types of fragmentation as possible outcomes in collisions among protoplanetary dust aggregates. We distinguish between eight combinations of porosity and mass ratio. For each of these cases, we present a complete collision model for dust-aggregate masses between 10^-12 and 10^2 g and collision velocities in the range 10^-4 to 10^4 cm/s for arbitrary porosities. This model comprises the collisional outcome, the mass(es) of the resulting aggregate(s) and their porosities. We present the first complete collision model for protoplanetary dust. This collision model can be used for the determination of the dust-growth rate in protoplanetary disks.

In dense molecular clouds collisions between dust grains alter the ISM-dust size distribution. We study this process by inserting the results from detailed numerical simulations of two colliding dust aggregates into a coagulation model that computes the dust size distribution with time. All collisional outcomes -- sticking, fragmentation (shattering, breakage, and erosion) -- are included and the effects on the internal structure of the aggregates are also tabulated. The dust aggregate evolution model is applied to an homogeneous and static cloud of temperature 10 K and gas densities between 10^3 and 10^7 cm^-3. The coagulation is followed locally on timescales of ~10^7 yr. We find that the growth can be divided into two stages: a growth dominated phase and a fragmentation dominated phase. Initially, the mass distribution is relatively narrow and shifts to larger sizes with time. At a certain point, dependent on the material properties of the grains as well as on the gas density, collision velocities will become sufficiently energetic to fragment particles, halting the growth and replenishing particles of lower mass. Eventually, a steady state is reached, where the mass distribution is characterized by a mass spectrum of approximately equal amount of mass per logarithmic size bin. The amount of growth that is achieved depends on the cloud's lifetime. If clouds exist on free-fall timescales the effects of coagulation on the dust size distribution are very minor. On the other hand, if clouds have long-term support mechanisms, the impact of coagulation is important, resulting in a significant decrease of the opacity on timescales longer than the initial collision timescale between big grains.

We present a method which extends Monte Carlo studies to situations that require a large dynamic range in particle number. The underlying idea is that, in order to calculate the collisional evolution of a system, some particle interactions are more important than others and require more resolution, while the behavior of the less important, usually of smaller mass, particles can be considered collectively. In this approximation groups of identical particles, sharing the same mass and structural parameters, operate as one unit. The amount of grouping is determined by the zoom factor -- a free parameter that determines on which particles the computational effort is focused. Two methods for choosing the zoom factors are discussed: the `equal mass method,' in which the groups trace the mass density of the distribution, and the `distribution method,' which additionally follows fluctuations in the distribution. Both methods achieve excellent correspondence with analytic solutions to the Smoluchowski coagulation equation. The grouping method is furthermore applied to simulations involving runaway kernels, where the particle interaction rate is a strong function of particle mass, and to situations that include catastrophic fragmentation. For the runaway simulations previous predictions for the decrease of the runaway timescale with the initial number of particles ${\cal N}$ are reconfirmed, extending ${\cal N}$ to $10^{160}$. Astrophysical applications include modeling of dust coagulation, planetesimal accretion, and the dynamical evolution of stars in large globular clusters. The proposed method is a powerful tool to compute the evolution of any system where the particles interact through discrete events, with the particle properties characterized by structural parameters.

We present a mechanism for chondrules to stick together by means of compaction of a porous dust rim they sweep up as they move through the dusty nebula gas. It is shown that dust aggregates formed out of micron-sized grains stick to chondrules, forming a porous dust rim. When chondrules collide, this dust can be compacted by means of rolling motions within the porous dust layer. This mechanism dissipates the collisional energy, compacting the rim and allowing chondrules to stick. The structure of the obtained chondrule-dust agglomerates (referred to as compounds) then consists of three phases: chondrules, porous dust, and dust that has been compacted by collisions. Subsequently, these compounds accrete their own dust and collide with other compounds. The evolution of the compound size distribution and the relative importance of the phases is calculated by a Monte Carlo code. Growth ends, and a simulation is terminated when all the dust in the compounds has been compacted. Numerous runs are performed, reflecting the uncertainty in the physical conditions at the chondrule formation time. It is found that compounds can grow by 1-2 orders of magnitudes in radius, upto dm-sizes when turbulence levels are low. However, relative velocities associated with radial drift form a barrier for further growth. Earlier findings that the dust sweep-up by chondrules is proportional to their sizes are confirmed. We contrast two scenarios regarding how this dust evolved further towards the densely packed rims seen in chondrites.

In this note we present complete, closed-form expressions for random relative velocities between colliding particles of arbitrary size in nebula turbulence. These results are exact for very small particles (those with stopping times much shorter than the large eddy overturn time) and are also surprisingly accurate in complete generality (that is, also apply for particles with stopping times comparable to, or much longer than, the large eddy overturn time). We note that some previous studies may have adopted previous simple expressions, which we find to be in error regarding the size dependence in the large particle regime.

Context: Sticking of colliding dust particles through van der Waals forces is the first stage in the grain growth process in protoplanetary disks, eventually leading to the formation of comets, asteroids and planets. A key aspect of the collisional evolution is the coupling between dust and gas motions, which depends on the internal structure (porosity) of aggregates. Aims: To quantify the importance of the internal structure on the collisional evolution of particles, and to create a new coagulation model to investigate the difference between porous and compact coagulation in the context of a turbulent protoplanetary disk. Methods: We have developed simple prescriptions for the collisional evolution of porosity of grain-aggregates in grain-grain collisions. Three regimes can then be distinguished: `hit-and-stick' at low velocities, with an increase in porosity; compaction at intermediate velocities, with a decrease of porosity; and fragmentation at high velocities. (..) Results: (..) We can discern three different stages in the particle growth process (..) We find that when compared to standard, compact models of coagulation, porous growth delays the onset of settling, because the surface area-to-mass ratio is higher, a consequence of the build-up of porosity during the initial stages. As a result, particles grow orders of magnitudes larger in mass before they rain-out to the mid-plane. Depending on the turbulent viscosity and on the position in the nebula, aggregates can grow to (porous) sizes of ~ 10 cm in a few thousand years. We also find that collisional energies are higher than in the limited PCA/CCA fractal models, thereby allowing aggregates to restructure. It is concluded that the microphysics of collisions plays a key role in the growth process.

This paper presents results from modelling 450 micron and 850 micron continuum and HCO+ line observations of three distinct cores of an infrared dark cloud (IRDC) directed toward the W51 GMC. In the sub-mm continuum these cores appear as bright, isolated emission features. One of them coincides with the peak of 8.3 micron extinction as measured by the Midcourse Space Experiment satellite. Detailed radiative transfer codes are applied to constrain the cores' physical conditions to address the key question: Do these IRDC-cores harbour luminous sources? The results of the continuum model, expressed in the $\chi^2$ quality-of-fit parameter, are also constrained by the absence of 100 micron emission from IRAS. For the sub-mm emission peaks this shows that sources of 300 solar luminosities are embedded within the cores. For the extinction peak, the combination of continuum and HCO+ line modelling indicates that a heating source is present as well. Furthermore, the line model provides constraints on the clumpiness of the medium. All three cores have similar masses of about 70-150 solar masses and similar density structures. The extinction peak differs from the other two cores by hosting a much weaker heating source, and the sub-mm emission core at the edge of the IRDC deviates from the other cores by a higher internal clumpiness.