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.

We present the first 2D hydrodynamical finite volume simulations in which dust is fully coupled with the gas, including its back-reaction onto it, and at the same time the dust size is evolving according to coagulation and fragmentation based on a sub-grid model. The aim of this analysis is to present the differences occurring when dust evolution is included relative to simulations with fixed dust size, with and without an embedded Jupiter-mass planet that triggers gap formation. We use the two-fluid polar Godunov-type code RoSSBi developed by Surville et al. 2016 combined with a new local sub-grid method for dust evolution based on the model by Birnstiel et al. 2012. We find striking differences between simulations with variable and fixed dust sizes. The timescales for dust depletion differ significantly and yield a completely different evolution of the dust surface density. In general sharp features such as pile-ups of dust in the inner disk and near gap edges, when a massive planet is present, become much weaker. This has important implications on the interpretation of observed substructure in disks, suggesting that the presence of a massive planet does not necessarily cause sharp gaps and rings in the dust component. Also, particles with different dust sizes show a different distribution, pointing to the importance of multi-wavelength synthetic observations in order to compare with observations by ALMA and other instruments. We also find that simulations adopting fixed intermediate particle sizes, in the range $10^{-2} - 10^{-1}$ cm, best approximate the surface density evolution seen in simulations with dust evolution.

It is believed that satellites of giant planets form in circumplanetary disks. Many of the previous contributions assumed that their formation process proceeds similarly to rocky planet formation, via accretion of the satellite seeds, called satellitesimals. However, the satellitesimal formation itself poses a nontrivial problem as the dust evolution in the circumplanetary disk is heavily impacted by fast radial drift and thus dust growth to satellitesimals is hindered. To address this problem, we connected state-of-the-art hydrodynamical simulations of a circumplanetary disk around a Jupiter-mass planet with dust growth and drift model in a post-processing step. We found that there is an efficient pathway to satellitesimal formation if there is a dust trap forming within the disk. Thanks to the natural existence of an outward gas flow region in the hydrodynamical simulation, a significant dust trap arises at the radial distance of 85~R$_{\rm J}$ from the planet, where the dust-to-gas ratio becomes high enough to trigger streaming instability. The streaming instability leads to the efficient formation of the satellite seeds. Because of the constant infall of material from the circumstellar disk and the very short timescale of dust evolution, the circumplanetary disk acts as a satellitesimal factory, constantly processing the infalling dust to pebbles that gather in the dust trap and undergo the streaming instability.

Models of dust coagulation and subsequent planetesimal formation are usually computed on the backdrop of an already fully formed protoplanetary disk model. At the same time, observational studies suggest that planetesimal formation should start early, possibly even before the protoplanetary disk is fully formed. In this paper, we investigate under which conditions planetesimals already form during the disk buildup stage, in which gas and dust fall onto the disk from its parent molecular cloud. We couple our earlier planetesimal formation model at the water snow line to a simple model of disk formation and evolution. We find that under most conditions planetesimals only form after the buildup stage when the disk becomes less massive and less hot. However, there are parameters for which planetesimals already form during the disk buildup. This occurs when the viscosity driving the disk evolution is intermediate ($\alpha_v \sim 10^{-3}-10^{-2}$) while the turbulent mixing of the dust is reduced compared to that ($\alpha_t \lesssim 10^{-4}$), and with the assumption that water vapor is vertically well-mixed with the gas. Such $\alpha_t \ll \alpha_v$ scenario could be expected for layered accretion, where the gas flow is mostly driven by the active surface layers, while the midplane layers, where most of the dust resides, are quiescent.

Planetesimal formation stage represents a major gap in our understanding of the planet formation process. The late-stage planet accretion models typically make arbitrary assumptions about planetesimals and pebbles distribution while the dust evolution models predict that planetesimal formation is only possible at some orbital distances. We want to test the importance of water snow line for triggering formation of the first planetesimals during the gas-rich phase of protoplanetary disk, when cores of giant planets have to form. We connect prescriptions for gas disk evolution, dust growth and fragmentation, water ice evaporation and recondensation, as well as transport of both solids and water vapor, and planetesimal formation via streaming instability into a single, one-dimensional model for protoplanetary disk evolution. We find that processes taking place around the snow line facilitate planetesimal formation in two ways. First, due to the change of sticking properties between wet and dry aggregates, there is a "traffic jam" inside of the snow line that slows down the fall of solids onto the star. Second, ice evaporation and outward diffusion of water followed by its recondensation increases the abundance of icy pebbles that trigger planetesimal formation via streaming instability just outside of the snow line. Planetesimal formation is hindered by growth barriers and radial drift and thus requires particular conditions to take place. Snow line is a favorable location where planetesimal formation is possible for a wide range of conditions, but still not in every protoplanetary disk model. This process is particularly promoted in large, cool disks with low intrinsic turbulence and increased initial dust-to-gas ratio.

The consistency of planet formation models suffers from the disconnection between the regime of small and large bodies. This is primarily caused by so-called growth barriers: the direct growth of larger bodies is halted at centimetre-sized objects and particular conditions are required for the formation of larger, gravitationally bound planetesimals. We aim to connect models of dust evolution and planetesimal formation to identify regions of protoplanetary discs that are favourable for the formation of kilometre-sized bodies and the first planetary embryos. We combine semi-analytical models of viscous protoplanetary disc evolution, dust growth and drift including backreaction of the dust particles on the gas, and planetesimal formation via the streaming instability into one numerical code. We investigate how planetesimal formation is affected by the mass of the protoplanetary disc, its initial dust content, and the stickiness of dust aggregates. We find that the dust growth and drift leads to a global redistribution of solids. The pile-up of pebbles in the inner disc provides local conditions where the streaming instability is effective. Planetesimals form in an annulus with its inner edge lying between 0.3 AU and 1 AU and its width ranging from 0.3 AU to 3 AU. The resulting surface density of planetesimals follows a radial profile that is much steeper than the initial disc profile. These results support formation of terrestrial planets in the solar system from a narrow annulus of planetesimals, which reproduces their peculiar mass ratios.

Streaming instability can be a very efficient way of overcoming growth and drift barriers to planetesimal formation. However, it was shown that strong clumping, which leads to planetesimal formation, requires a considerable number of large grains. State-of-the-art streaming instability models do not take into account realistic size distributions resulting from the collisional evolution of dust. We investigate whether a sufficient quantity of large aggregates can be produced by sticking and what the interplay of dust coagulation and planetesimal formation is. We develop a semi-analytical prescription of planetesimal formation by streaming instability and implement it in our dust coagulation code based on the Monte Carlo algorithm with the representative particles approach. We find that planetesimal formation by streaming instability may preferentially work outside the snow line, where sticky icy aggregates are present. The efficiency of the process depends strongly on local dust abundance and radial pressure gradient, and requires a super-solar metallicity. If planetesimal formation is possible, the dust coagulation and settling typically need ~100 orbits to produce sufficiently large and settled grains and planetesimal formation lasts another ~1000 orbits. We present a simple analytical model that computes the amount of dust that can be turned into planetesimals given the parameters of the disk model.

Simple toy models are often not sufficient to cover the complexity of the dust coagulation process, and a number of numerical approaches are therefore used, among which integration of the Smoluchowski equation and various versions of Monte Carlo algorithm are the most popular. In this paper, we directly compare the Smoluchowski and Monte Carlo approaches and we find a general agreement for most of the coagulation problems. However, for the sweep-up growth driven by the "lucky" breakthrough mechanism, the methods exhibit very different resolution dependencies. With too few mass bins, the Smoluchowski algorithm tends to overestimate the growth rate and the probability of breakthrough. The Monte Carlo method is less resolution dependent in the growth timescale aspect but it tends to underestimate the breakthrough chance due to its limited dynamic mass range. We discuss the features and drawbacks of both the approaches, which may limit their astrophysical applications.

The early stages of planet formation are still not well understood. Coagulation models have revealed numerous obstacles to the dust growth, such as the bouncing, fragmentation and radial drift barriers. We study the interplay between dust coagulation and drift in order to determine the conditions in protoplanetary disk that support the formation of planetesimals. We focus on planetesimal formation via sweep-up and investigate whether it can take place in a realistic protoplanetary disk. We have developed a new numerical model that resolves spatial distribution of dust in the radial and vertical dimension. The model uses representative particles approach to follow the dust evolution in protoplanetary disk. The coagulation and fragmentation of solids is taken into account using Monte Carlo method. A collision model adopting the mass transfer effect, that can occur for different-sized dust aggregate collisions, is implemented. We focus on a protoplanetary disk including a pressure bump caused by a steep decline of turbulent viscosity around the snow line. Our results show that sufficient resolution of the vertical disk structure in dust coagulation codes is necessary to obtain adequately short growth timescales, especially in the case of a low turbulence region. We find that a sharp radial variation of the turbulence strength at the inner edge of dead zone promotes planetesimal formation in several ways. It provides a pressure bump that efficiently prevents the dust from drifting inwards. It also causes a radial variation in the size of aggregates at which growth barriers occur, favoring the growth of large aggregates via sweeping up of small particles. In our model, by employing an ad hoc alpha viscosity change near the snow line, it is possible to grow planetesimals by incremental growth on timescales of approximately 10^5 years.