If we want to understand planetesimal formation, the only data set we have is our own Solar System. It is particularly interesting as it is so far the only planetary system we know of that developed life. Understanding the conditions under which the Solar Nebula evolved is crucial in order to understand the different processes in the disk and the subsequent dynamical interaction between (proto-)planets, once the gas disk is gone. Protoplanetary disks provide a plethora of different parameters to explore. The question is whether this parameter space can be constrained, allowing simulations to reproduce the Solar System. Models and observations of planet formation provide constraints on the initial planetesimal mass in certain regions of the Solar Nebula. By making use of pebble flux-regulated planetesimal formation, we perform a parameter study with nine different disk parameters like the initial disk mass, initial disk size, initial dust-to-gas ratio, turbulence level, and more. We find that the distribution of mass in planetesimals in the disk depends on the planetesimal formation timescale and the pebbles' drift timescale. Multiple disk parameters can influence pebble properties and thus planetesimal formation. However, it is still possible to draw some conclusions on potential parameter ranges. Pebble flux-regulated planetesimal formation seems to be very robust, allowing simulations with a wide range of parameters to meet the initial planetesimal constraints for the Solar Nebula. I.e., it does not require a lot of fine tuning.

In this work, we study how the dust coagulation/fragmentation will influence the evolution and observational appearances of vortices induced by a massive planet embedded in a low viscosity disk by performing global 2D high-resolution hydrodynamical simulations. Within the vortex, due to its higher gas surface density and steeper pressure gradients, dust coagulation, fragmentation and drift (to the vortex center) are all quite efficient, producing dust particles ranging from micron to $\sim 1.0\ {\rm cm}$, as well as overall high dust-to-gas ratio (above unity). In addition, the dust size distribution is quite non-uniform inside the vortex, with the mass weighted average dust size at the vortex center ($\sim 4.0$ mm) being a factor of $\sim10$ larger than other vortex regions. Both large ($\sim$ mm) and small (tens of micron) particles contribute strongly to affect the gas motion within the vortex. As such, we find that the inclusion of dust coagulation has a significant impact on the vortex lifetime and the typical vortex lifetime is about 1000 orbits. After the initial gaseous vortex is destroyed, the dust spreads into a ring with a few remaining smaller gaseous vortices with a high dust concentration and a large maximum size ($\sim$ mm). At late time, the synthetic dust continuum images for the coagulation case show as a ring inlaid with several hot spots at 1.33 mm band, while only distinct hot spots remain at 7.0 mm.

Context. The water snowline divides dry and icy solid material in protoplanetary disks, and has been thought to significantly affect planet formation at all stages. If dry particles break up more easily than icy ones, then the snowline causes a traffic jam, because small grains drift inward at lower speeds than larger pebbles. Aims. We aim to evaluate the effect of high dust concentrations around the snowline onto the gas dynamics. Methods. Using numerical simulations, we model the global radial evolution of an axisymmetric protoplanetary disk. Our model includes particle growth, evaporation and recondensation of water, and the back-reaction of dust onto the gas, taking into account the vertical distribution of dust particles. Results. We find that the dust back-reaction can stop and even reverse the net flux of gas outside the snowline, decreasing the gas accretion rate onto the star to under $50\%$ of its initial value. At the same time the dust accumulates at the snowline, reaching dust-to-gas ratios of $\epsilon \gtrsim 0.8$, and delivers large amounts of water vapor towards the inner disk, as the icy particles cross the snowline. However, the accumulation of dust at the snowline and the decrease in the gas accretion rate only take place if the global dust-to-gas ratio is high ($\varepsilon_0 \gtrsim 0.03$), if the viscous turbulence is low ($\alpha_\nu \lesssim 10^{-3} $), if the disk is large enough ($r_c \gtrsim 100\, \textrm{au}$), and only during the early phases of the disk evolution ($t \lesssim 1\, \textrm{Myr}$). Otherwise the dust back-reaction fails to perturb the gas motion.

Tidal interactions between the embedded planets and their surrounding protoplanetary disks are often postulated to produce the observed complex dust substructures, including rings, gaps, and asymmetries. In this Letter, we explore the consequences of dust coagulation on the dust dynamics and ring morphology. Coagulation of dust grains leads to dust size growth which, under typical disk conditions, produces faster radial drifts, potentially threatening the dust ring formation. Utilizing 2D hydrodynamical simulations of protoplanetary disks which include a full treatment of dust coagulation, we find that if the planet does not open a gap quickly enough, the formation of an inner ring is impeded due to dust coagulation and subsequent radial drift. Furthermore, we find that a "buildup" of sub-mm sized grains often appears in the dust emission at the outer edge of the dust disk.

Dust growth is often neglected when building models of protoplanetary disks due to its complexity and computational expense. However, it does play a major role in shaping the evolution of protoplanetary dust and planet formation. In this paper, we present a numerical model coupling 2-D hydrodynamic evolution of a protoplanetary disk, including a Jupiter-mass planet, and dust coagulation. This is obtained by including multiple dust fluids in a single grid-based hydrodynamic simulation and solving the Smoluchowski equation for dust coagulation on top of solving for the hydrodynamic evolution. We find that fragmentation of dust aggregates trapped in a pressure bump outside of the planetary gap leads to an enhancement in density of small grains. We compare the results obtained from the full coagulation treatment to the commonly used, fixed dust size approach and to previously applied, less computationally intensive methods for including dust coagulation. We find that the full coagulation results cannot be reproduced using the fixed-size treatment, but some can be mimicked using a relatively simple method for estimating the characteristic dust size in every grid cell.

Recent high-resolution interferometric observations of protoplanetary disks at (sub-)millimeter wavelengths reveal omnipresent substructures, such as rings, spirals, and asymmetries. A detailed investigation of eight rings detected in five disks by the DSHARP survey came to the conclusion that all rings are just marginally optically thick with optical depths between 0.2 and 0.5 at a wavelength of 1.25 mm. This surprising result could either be coincidental or indicate that the optical depth in all of the rings is regulated by the same process. We investigated if ongoing planetesimal formation could explain the "fine-tuned" optical depths in the DSHARP rings by removing dust and transforming it into "invisible" planetesimals. We performed a one-dimensional simulation of dust evolution in the second dust ring of the protoplanetary disk around HD 163296, including radial transport of gas and dust, dust growth and fragmentation, and planetesimal formation via gravitational collapse of sufficiently dense pebble concentrations. We show that planetesimal formation can naturally explain the observed optical depths if streaming instability regulates the midplane dust-to-gas ratio to unity. Furthermore, our simple monodisperse analytical model supports the hypothesis that planetesimal formation in dust rings should universally limit their optical depth to the observed range.

The growth of solids from sub-micron to millimeter and centimeter sizes is the early step toward the formation of planets inside protoplanetary disks (PPDs). However, such processes and their potential impact on the later stages of solid growth are still poorly understood. In this work, we test the hypothesis that most disks contain at least one ringed structure with a relative small radius. We have carried out a large family of 1D two-fluid (gas+dust) hydrodynamical simulations by evolving the gas and dust motion self-consistently while allowing dust size to evolve via coagulation and fragmentation. We investigate the joint effects of ringed structures and dust size growth on the overall sub-millimeter and millimeter (mm) flux and spectral index of PPDs. Ringed structures slow down the dust radial drift and speed up the dust growth. In particular, we find that those unresolved disks with a high fragmentation velocity ($\sim10\ {\rm m~s^{-1}}$) and a high dust surface density ($\sim10\ {\rm g\ cm^{-2}}$ in the ring) can have mm spectral indices as low as $\sim2.0$, consistent with mm observations of faint disks in nearby star forming regions. Furthermore, disks with more than one ringed structure can potentially reproduce brighter disks with spectral indices lower than $\sim2.5$. Future multi-wavelength high-resolution observations of these low spectral index sources can be used to test the existence of the ringed structures in the unresolved disks and differentiate the effects of dust size growth from optical depth.

RW Aur A has experienced various dimming events in the last years, decreasing its brightness by $\sim 2\ \textrm{mag}$ for periods of months to years. Multiple observations indicate that a high concentration of dust grains, from the protoplanetary disk's inner regions, is blocking the starlight during these events. We propose a new mechanism that can send large amounts of dust close to the star on short timescales, through the reactivation of a dead zone in the protoplanetary disk. Using numerical simulations we model the accretion of gas and dust, along with the growth and fragmentation of particles in this scenario. We find that after the reactivation of the dead zone, the accumulated dust is rapidly accreted towards the star in around 15 years, at rates of $\dot{M}_\textrm{d} = 6 \times 10^{-6}\, \textrm{M}_\odot/\textrm{yr}$ and reaching dust-to-gas ratios of $\epsilon \approx 5$, preceding an increase in the gas accretion by a few years. This sudden rise of dust accretion can provide the material required for the dimmings, although the question of how to put the dust into the line of sight remains open to speculation.

Context. Despite the recent discovery of spiral-shaped features in protoplanetary discs in the near-infrared and millimetric wavelengths, there is still an active discussion to understand how they formed. In fact, the spiral waves observed in discs around young stars can be due to different physical mechanisms: planet/companion torques, gravitational perturbations or illumination effects. Aims. We study the spirals formed in the gaseous phase due to two diametrically opposed shadows cast at fixed disc locations. The shadows are created by an inclined non-precessing disc inside the cavity, which is assumed to be optically thick. In particular, we analyse the effect of these spirals on the dynamics of the dust particles and discuss their detectability in transition discs. Methods. We perform gaseous hydrodynamical simulations with shadows, then we compute the dust evolution on top of the gaseous distribution, and finally we produce synthetic ALMA observations of the dust emission based on radiative transfer calculations. Results. Our main finding is that mm- to cm-sized dust particles are efficiently trapped inside the shadow-triggered spirals. We also observe that particles of various sizes starting at different stellocentric distances are well mixed inside these pressure maxima. This dynamical effect would favour grain growth and affect the resulting composition of planetesimals in the disc. In addition, our radiative transfer calculations show spiral patterns in the disc at 1.6 {\mu}m and 1.3 mm. Due to their faint thermal emission (compared to the bright inner regions of the disc) the spirals cannot be detected with ALMA. Our synthetic observations prove however that shadows are observable as dips in the thermal emission.

Recent high angular resolution observations of protoplanetary disks at different wavelengths have revealed several kinds of structures, including multiple bright and dark rings. Embedded planets are the most used explanation for such structures, but there are alternative models capable of shaping the dust in rings as it has been observed. We assume a disk around a Herbig star and investigate the effect that ice lines have on the dust evolution, following the growth, fragmentation, and dynamics of multiple dust size particles, covering from 1 $\mu$m to 2 m sized objects. We use simplified prescriptions of the fragmentation velocity threshold, which is assumed to change radially at the location of one, two, or three ice lines. We assume changes at the radial location of main volatiles, specifically H$_2$O, CO$_2$, and NH$_3$. Radiative transfer calculations are done using the resulting dust density distributions in order to compare with current multiwavelength observations. We find that the structures in the dust density profiles and radial intensities at different wavelengths strongly depend on the disk viscosity. A clear gap of emission can be formed between ice lines and be surrounded by ring-like structures, in particular between the H$_2$O and CO$_2$ (or CO). The gaps are expected to be shallower and narrower at millimeter emission than at near-infrared, opposite to model predictions of particle trapping. In our models, the total gas surface density is not expected to show strong variations, in contrast to other gap-forming scenarios such as embedded giant planets or radial variations of the disk viscosity.

Context. Ice lines are suggested to play a significant role in grain growth and planetesimal formation in protoplanetary disks. Evaporation fronts directly influence the gas and ice abundances of volatile species in the disk and therefore the coagulation physics and efficiency and the chemical composition of the resulting planetesimals. Aims. In this work we investigate the influence of the existence of the CO ice line on the particle growth and on the distribution of CO in the disk. Methods. We include the possibility of tracking the CO content and/or other volatiles in particles and in the gas in our existing dust coagulation and disk evolution model and developed a method for evaporation and condensation of CO using the Hertz-Knudsen equation. Our model does not include fragmentation, yet, which will be part of further investigations. Results. We find no enhanced grain growth just outside the ice line where the particle size is limited by radial drift. Instead we find a depletion of solid material inside the ice line which is solely due to evaporation of the CO. Such a depression inside the ice line may be observable and may help to quantify the processes described in this work. Furthermore, we find that the viscosity and diffusivity of the gas heavily influence the re-distribution of vaporized CO at the ice line and can lead to an increase in the CO abundance by up to a factors of a few in the region just inside the ice line. Depending on the strength of the gaseous transport mechanisms the position of the ice line in our model can change by up to 10 AU and consequently, the temperature at that location can range from 21 K to 23 K.

Solving the mystery of the origin of chondrules is one of the most elusive goals in the field of meteoritics. Recently the idea of planet(esimal) collisions releasing splashes of lava droplets, long considered out of favor, has been reconsidered as a possible origin of chondrules by several papers. One of the main problems with this idea is the lack of quantitative and simple models that can be used to test this scenario by directly comparing to the many known observables of chondrules. In Paper I of this series we presented a simple thermal evolution model of a spherically symmetric expanding cloud of molten lava droplets that is assumed to emerge from a collision between two planetesimals. In the present paper, number II of this series, we use this model to calculate whether or not volatile elements such as Na and K will remain abundant in these droplets or whether they will get depleted due to evaporation. The high density of the droplet cloud (e.g. small distance between adjacent droplets) causes the vapor to quickly reach saturation pressure and thus shutting down further evaporation. We show to which extent, and under which conditions, this keeps the abundances of these elements high, as is seen in chondrules. We find that for most parameters of our model (cloud mass, expansion velocity, initial temperature) the volatile elements Mg, Si and Fe remain entirely in the chondrules. The Na and K abundances inside the droplets will initially stay mostly at their initial values due to the saturation of the vapor pressure, but at some point start to drop due to the cloud expansion. However, as soon as the temperature starts to decrease, most or all of the vapor recondenses again. At the end the Na and K elements retain most of their initial abundances, albeit occasionally somewhat reduced, depending on the parameters of the expanding cloud model.

The formation of chondrules is one of the oldest unsolved mysteries in meteoritics and planet formation. Recently an old idea has been revived: the idea that chondrules form as a result of collisions between planetesimals in which the ejected molten material forms small droplets which solidify to become chondrules. Pre-melting of the planetesimals by radioactive decay of 26Al would help producing sprays of melt even at relatively low impact velocity. In this paper we study the radiative cooling of a ballistically expanding spherical cloud of chondrule droplets ejected from the impact site. We present results from a numerical radiative transfer models as well as analytic approximate solutions. We find that the temperature after the start of the expansion of the cloud remains constant for a time t_cool and then drops with time t approximately as T ~ T_0[(3/5)t/t_cool+ 2/5]^(-5/3) for t>t_cool. The time at which this temperature drop starts t_cool depends via an analytical formula on the mass of the cloud, the expansion velocity and the size of the chondrule. During the early isothermal expansion phase the density is still so high that we expect the vapor of volatile elements to saturate so that no large volatile losses are expected.

In recent years many models of chondrule formation have been proposed. One of those models is the processing of dust in shock waves in protoplanetary disks. In this model, the dust and the chondrule precursors are overrun by shock waves, which heat them up by frictional heating and thermal exchange with the gas. In this paper we reanalyze the nebular shock model of chondrule formation and focus on the downstream boundary condition. We show that for large-scale plane-parallel chondrule-melting shocks the postshock equilibrium temperature is too high to avoid volatile loss. Even if we include radiative cooling in lateral directions out of the disk plane into our model (thereby breaking strict plane-parallel geometry) we find that for a realistic vertical extent of the solar nebula disk the temperature decline is not fast enough. On the other hand, if we assume that the shock is entirely optically thin so that particles can radiate freely, the cooling rates are too high to produce the observed chondrules textures. Global nebular shocks are therefore problematic as the primary sources of chondrules.