The outer Solar System object (486958) Arrokoth (provisional designation 2014 MU$_{69}$) has been largely undisturbed since its formation. We study its surface composition using data collected by the New Horizons spacecraft. Methanol ice is present along with organic material, which may have formed through radiation of simple molecules. Water ice was not detected. This composition indicates hydrogenation of carbon monoxide-rich ice and/ or energetic processing of methane condensed on water ice grains in the cold, outer edge of the early Solar System. There are only small regional variations in color and spectra across the surface, suggesting Arrokoth formed from a homogeneous or well-mixed reservoir of solids. Microwave thermal emission from the winter night side is consistent with a mean brightness temperature of 29$\pm$5 K.

Protoplanetary disks are dynamic objects, within which dust grains and gas are expected to be redistributed over large distances. Evidence for this redistribution is seen both in other protoplanetary disks and in our own Solar System, with high-temperature materials thought to originate close to the central star found in the cold, outer regions of the disks. While models have shown this redistribution is possible through a variety of mechanisms, these models have generally ignored the possible growth of solids via grain-grain collisions that would occur during transit. Here we investigate the interplay of coagulation and radial and vertical transport of solids in protoplanetary disks, considering cases where growth is limited by bouncing or by fragmentation. We find that in all cases, growth effectively limits the ability for materials to be carried outward or preserved at large distances from the star. This is due to solids being incorporated into large aggregates which drift inwards rapidly under the effects of gas drag. We discuss the implications for mixing in protoplanetary disks, and how the preservation of high temperature materials in outer disks may require structures or outward flow patterns to avoid them being lost via radial drift.

CO is the most widely used gas tracer of protoplanetary disks. Its abundance is usually assumed to be an interstellar ratio throughout the warm molecular layer of the disk. But recent observations of low CO gas abundance in many protoplanetary disks challenge our understanding of physical and chemical evolutions in disks. Here we investigate the CO abundance structures in four well-studied disks and compare their structures with predictions of chemical processing of CO and transport of CO ice-coated dust grains in disks. We use spatially resolved CO isotopologue line observations and detailed thermo-chemical models to derive CO abundance structures. We find that the CO abundance varies with radius by an order of magnitude in these disks. We show that although chemical processes can efficiently reduce the total column of CO gas within 1 Myr under an ISM level of cosmic-ray ionization rate, the depletion mostly occurs at the deep region of a disk. Without sufficient vertical mixing, the surface layer is not depleted enough to reproduce weak CO emissions observed. The radial profiles of CO depletion in three disks are qualitatively consistent with predictions of pebble formation, settling, and drifting in disks. But the dust evolution alone cannot fully explain the high depletion observed in some disks. These results suggest that dust evolution may play a significant role in transporting volatile materials and a coupled chemical-dynamical study is necessary to understand what raw materials are available for planet formation at different distances from the central star.

There is evidence that protoplanetary disks--including the protosolar one--contain crystalline dust grains on spatial scales where the dust temperature is lower than the threshold value for their formation through thermal annealing of amorphous interstellar silicates. We interpret these observations in terms of an extended, magnetocentrifugally driven disk wind that transports grains from the inner disk--where they are thermally processed by the stellar radiation after being uplifted from the disk surfaces--to the outer disk regions. For any disk radius $r$ there is a maximum grain size $a_\mathrm{max}(r)$ that can be uplifted from that location: grains of size $a$$\ll$$a_\mathrm{max}$ are carried away by the wind, whereas those with $a$$\lesssim$$a_\mathrm{max}$ reenter the disk at larger radii. A significant portion of the reentering grains converge to--and subsequently accumulate in--a narrow region just beyond $r_\mathrm{max}(a)$, the maximum radius from which grains of size $a$ can be uplifted. We show that this model can account for the inferred crystallinity fractions in classical T Tauri and Herbig Ae disks and for their indicated near constancy after being established early in the disk evolution. It is also consistent with the reported radial gradients in the mean grain size, crystallinity, and crystal composition. In addition, this model yields the properties of the grains that remain embedded in the outflows from protoplanetary disks and naturally explains the inferred persistence of small grains in the surface layers of these disks.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

Amorphous ice has long been invoked as a means for trapping extreme volatiles into solids, explaining the abundances of these species in comets and planetary atmospheres. Experiments have shown that such trapping is possible and have been used to estimate the abundances of each species in primitive ices after they formed. However, these experiments have been carried out at deposition rates which exceed those expected in a molecular cloud or solar nebula by many orders of magnitude. Here we develop a numerical model which reproduces the experimental results and apply it to those conditions expected in molecular clouds and protoplanetary disks. We find that two regimes of ice trapping exist: `burial trapping' where the ratio of trapped species to water in the ice reflects that same ratio in the gas and `equilibrium trapping' where the ratio in the ice depends only on the partial pressure of the trapped species in the gas. The boundary between these two regimes is set by both the temperature and rate of ice deposition. Such effects must be accounted for when determining the source of trapped volatiles during planet formation.

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.

Current models of (exo)planet formation often rely on a large influx of so-called `pebbles' from the outer disk into the planet formation region. In this paper, we investigate how the formation of pebbles in the cold outer regions of protoplanetary disks and their subsequent migration to the inner disk can alter the gas-phase CO distribution both interior and exterior to the midplane CO snowline. By simulating the resulting CO abundances in the midplane as well as the warm surface layer, we identify observable signatures of large-scale pebble formation and migration that can be used as `smoking guns' for these important processes. Specifically, we find that after $1\mathrm{~Myr}$, the formation and settling of icy pebbles results in the removal of up to $80\%$ of the CO vapor in the warm ($T>22\mathrm{~K}$) disk layers outside the CO snowline, while the radial migration of pebbles results in the generation of a plume of CO vapor interior the snowline, increasing the CO abundance by a factor ${\sim}2{-}6$ depending on the strength of the turbulence and the sizes of the individual pebbles. The absence of this plume of CO vapor in young nearby disks could indicate efficient conversion of CO into a more refractory species, or a reduction in the radial mass flux of pebbles by, for example, disk inhomogeneities or early planetesimal formation.

In a protoplanetary disk, dust aggregates in the $\mu$m to mm size range possess mean collision velocities of 10 to 60 ms$^{-1}$ with respect to dm- to m-size bodies. We performed laboratory collision experiments to explore this parameter regime and found a size- and velocity-dependent threshold between erosion and growth. By using a local Monte Carlo coagulation calculation and complementary a simple semi-analytical timescale approach, we show that erosion considerably limits particle growth in protoplanetary disks and leads to a steady-state dust-size distribution from $\mu$m to dm sized particles.

With several short-period, Earth-mass planets in the habitable zone, the TRAPPIST-1 system potentially allows litho-panspermia to take place on very short timescales. We investigate the efficiency and speed of inter-planetary material transfer resulting from impacts onto the habitable zone planets. By simulating trajectories of impact ejecta from their moment of ejection until (re-)accretion, we find that transport between the habitable zone planets is fastest for ejection velocities around and just above planetary escape velocity. At these ejection velocities, ${\sim}10\%$ of the ejected material reaches another habitable zone planet within $10^2\mathrm{~yr}$, indicating litho-panspermia can be 4 to 5 orders of magnitude faster in TRAPPIST-1 than in the Solar System.

The processes that govern the evolution of dust and water (in the form of vapor or ice) in protoplanetary disks are intimately connected. We have developed a model that simulates dust coagulation, dust dynamics (settling, turbulent mixing), vapor diffusion, and condensation/sublimation of volatiles onto grains in a vertical column of a protoplanetary disk. We employ the model to study how dust growth and dynamics influence the vertical distribution of water vapor and water ice in the region just outside the radial snowline. Our main finding is that coagulation (boosted by the enhanced stickiness of icy grains) and the ensuing vertical settling of solids results in water vapor being depleted, but not totally removed, from the region above the snowline on a timescale commensurate with the vertical turbulent mixing timescale. Depending on the strength of the turbulence and the temperature, the depletion can reach factors of up to ${\sim}50$ in the disk atmosphere. In our isothermal column, this vapor depletion results in the vertical snowline moving closer to the midplane (by up to 2 gas scale heights) and the gas-phase C/O ratio above the vertical snowline increasing. Our findings illustrate the importance of dynamical effects and the need for understanding coevolutionary dynamics of gas and solids in planet-forming environments.

The disk mass is among the most important input parameter for every planet formation model to determine the number and masses of the planets that can form. We present an ALMA 887micron survey of the disk population around objects from 2 to 0.03Msun in the nearby 2Myr-old Chamaeleon I star-forming region. We detect thermal dust emission from 66 out of 93 disks, spatially resolve 34 of them, and identify two disks with large dust cavities of about 45AU in radius. Assuming isothermal and optically thin emission, we convert the 887micron flux densities into dust disk masses, hereafter Mdust. We find that the Mdust-Mstar relation is steeper than linear with power law indices 1.3-1.9, where the range reflects two extremes of the possible relation between the average dust temperature and stellar luminosity. By re-analyzing all millimeter data available for nearby regions in a self-consistent way, we show that the 1-3 Myr-old regions of Taurus, Lupus, and Chamaeleon I share the same Mdust-Mstar relation, while the 10Myr-old Upper Sco association has a steeper relation. Theoretical models of grain growth, drift, and fragmentation reproduce this trend and suggest that disks are in the fragmentation-limited regime. In this regime millimeter grains will be located closer in around lower-mass stars, a prediction that can be tested with deeper and higher spatial resolution ALMA observations.

In protoplanetary disks, the distribution and abundance of small (sub)micron grains are important for a range of physical and chemical processes. For example, they dominate the optical depth at short wavelengths and their surfaces are the sites of many important chemical reactions such as the formation of water. Based on their aerodynamical properties (i.e., their strong dynamical coupling with the surrounding gas) it is often assumed that these small grains are well-mixed with the gas. Our goal is to study the vertical (re)distribution of grains taking into account settling, turbulent diffusion, as well as collisions with other dust grains. Assuming a fragmentation-limited background dust population, we developed a Monte Carlo approach that follows single monomers as they move through a vertical column of gas and become incorporated in different aggregates as they undergo sticking and fragmenting collisions. We find that (sub)micron grains are not necessarily well-mixed vertically, but can become trapped in a thin layer with a scale-height that is significantly smaller than that of the gas. This collisional trapping occurs when the timescale for diffusion is comparable to or longer than the collision timescale in the midplane and its effect is strongest when the most massive particles in the size-distribution show significant settling. Based on simulations and analytical considerations we conclude that for typical dust-to-gas ratios and turbulence levels, the collisional trapping of small grains should be a relatively common phenomenon. The absence of trapping could then indicate a low dust-to-gas ratio, possibly because a large portion of the dust mass has been removed through radial drift or is locked up in planetesimals.

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.

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.

Observations with the Spitzer Space Telescope and the WISE satellite have revealed a prominent arc-like structure at 50" ($\simeq0.1$ pc) from the O9.5V/B0.5V system $\sigma$ Ori AB. We attribute this dust structure to the interaction of radiation pressure from the star with dust carried along by the IC 434 photo-evaporative flow of ionized gas from the dark cloud L1630. We have developed a quantitative model for the interaction of a dusty ionized flow with nearby (massive) stars where radiation pressure stalls dust, piling it up at an appreciable distance (> 0.1 pc), and force it to flow around the star. The model demonstrates that for the conditions in IC 434, the gas will decouple from the dust and will keep its original flow lines. We argue that this dust structure is the first example of a dust wave created by a massive star moving through the interstellar medium. Dust waves (and bow waves) stratify dust grains according to their radiation pressure opacity, which reflects the size distribution and composition of the grain material. Comparison of our model with observations implies that dust-gas coupling through Coulomb interaction is less important than previously thought, challenging our understanding of grain dynamics in hot, ionized regions of space. We describe the difference between dust (and bow) waves and classical bow shocks. We conclude that dust waves and bow waves should be common around stars showing the weak-wind phenomenon and that these structures are best observed at mid-IR to FIR wavelengths. In particular, dust waves and bow waves are most efficiently formed around weak-wind stars moving through a high density medium. Moreover, they provide a unique opportunity to study the direct interaction between a (massive) star and its immediate surroundings.

A prescription for the fragment size distribution resulting from dust grain collisions is essential when modelling a range of astrophysical systems, such as debris disks and planetary rings. While the slope of the fragment size distribution and the size of the largest fragment are well known, the behaviour of the distribution at the small size end is theoretically and experimentally poorly understood. This leads debris disk codes to generally assume a limit equal to, or below, the radiation blow-out size. We use energy conservation to analytically derive a lower boundary of the fragment size distribution for a range of collider mass ratios. Focussing on collisions between equal-sized bodies, we apply the method to debris disks. For a given collider mass, the size of the smallest fragments is found to depend on collision velocity, material parameters, and the size of the largest fragment. We provide a physically motivated recipe for the calculation of the smallest fragment, which can be easily implemented in codes for modelling collisional systems. For plausible parameters, our results are consistent with the observed predominance of grains much larger than the blow-out size in Fomalhaut's main belt and in the Herschel cold debris disks.

Aims: The aim of this work is to gain a deeper insight into how much different aggregate types are affected by erosion. Especially, it is important to study the influence of the velocity of the impacting projectiles. We also want to provide models for dust growth in protoplanetary disks with simple recipes to account for erosion effects. Methods: To study the erosion of dust aggregates we employed a molecular dynamics approach that features a detailed micro-physical model of the interaction of spherical grains. For the first time, the model has been extended by introducing a new visco-elastic damping force which requires a proper calibration. Afterwards, different sample generation methods were used to cover a wide range of aggregate types. Results: The visco-elastic damping force introduced in this work turns out to be crucial to reproduce results obtained from laboratory experiments. After proper calibration, we find that erosion occurs for impact velocities of 5 m/s and above. Though fractal aggregates as formed during the first growth phase are most susceptible to erosion, we observe erosion of aggregates with rather compact surfaces as well. Conclusions: We find that bombarding a larger target aggregate with small projectiles results in erosion for impact velocities as low as a few m/s. More compact aggregates suffer less from erosion. With increasing projectile size the transition from accretion to erosion is shifted to higher velocities. This allows larger bodies to grow through high velocity collisions with smaller aggregates.

A small number of young stellar objects show signs of a halo-like structure of optically thin dust. This halo or torus is located within a few AU of the star, but its origin has not yet been understood. A dynamically excited cloud of planetesimals colliding to eventually form dust could produce such a structure. The cause of the dynamical excitation could be one or more planets. This work investigates an inwardly migrating planet that is dynamically scattering planetesimals as a possible cause for the observed structures. If this mechanism is responsible, the observed halo-like structure could be used to infer the existence of planets in these systems. We present analytical estimates on the maximum inclination reached owing to dynamical interactions between planetesimals and a migrating planet. A symplectic integrator is used to simulate the effect of a migrating planet on a population of planetesimals. It is found that an inwardly migrating planet is only able to scatter the material it encounters to highly-inclined orbits if the material is on an eccentric orbit. Such eccentric orbits can be the result of resonance trapping and eccentricity pumping. Simulations show that for a certain range of migration rates and planet masses, resonance capture combined with planetary migration indeed causes the planetesimals to reach eccentric orbits and subsequently get scattered to highly-inclined orbits. The size distribution of the resulting dust is calculated determined to find the total mass and optical depth, which are found to compare reasonably well with the observed structures. Dynamical scattering of planetesimals caused by a planet migrating in, followed by the grinding down of these planetesimals to dust grains, appears to be a promising explanation for the inferred circumstellar dust clouds. Further study is needed to see if the haloes can be used to infer the presence of planets.