Search

Models of planetary core growth by either planetesimal or pebble accretion are traditionally disconnected from the models of dust evolution and formation of the first gravitationally-bound planetesimals. The state-of-the-art models typically start with massive planetary cores already present. We aim to study the formation and growth of planetary cores in a pressure bump, motivated by the annular structures observed in protoplanetary disks, starting with sub-micron-sized dust grains. We connect the models of dust coagulation and drift, planetesimal formation in the streaming instability, gravitational interactions between planetesimals, pebble accretion, and planet migration, into one uniform framework. We find that planetesimals forming early at the massive end of the size distribution grow quickly dominantly by pebble accretion. These few massive bodies grow on the timescales of ~100 000 years and stir the planetesimals formed later preventing the emergence of further planetary cores. Additionally, a migration trap occurs allowing for retention of the growing cores. Pressure bumps are favourable locations for the emergence and rapid growth of planetary cores by pebble accretion as the dust density and grain size are increased and the pebble accretion onset mass is reduced compared to a smooth-disk model., 15 pages, 11 figures. Replaced with published version

Many processes during the evolution of protoplanetary disks and during planet formation are highly sensitive to the sizes of dust particles that are present in the disk: The efficiency of dust accretion in the disk and volatile transport on dust particles, gravoturbulent instabilities leading to the formation of planetesimals, or the accretion of pebbles onto large planetary embryos to form giant planets are typical examples of processes that depend on the sizes of the dust particles involved. Furthermore, radiative properties like absorption or scattering opacities depend on the particle sizes. To interpret observations of dust in protoplanetary disks, a proper estimate of the dust particle sizes is needed. We present DustPy - A Python package to simulate dust evolution in protoplanetary disks. DustPy solves gas and dust transport including viscous advection and diffusion as well as collisional growth of dust particles. DustPy is written with a modular concept, such that every aspect of the model can be easily modified or extended to allow for a multitude of research opportunities., Accepted for publication in ApJ. Installation via "pip install dustpy"

Recent sub-arcsecond resolution surveys of the dust continuum emission from nearby protoplanetary disks showed a strong correlation between the sizes and luminosities of the disks. We aim to explain the origin of the (sub-)millimeter size-luminosity relation (SLR) between the $68\%$ effective radius ($r_{eff}$) of disks with their continuum luminosity ($L_{mm}$), with models of gas and dust evolution in a simple viscous accretion disk and radiative transfer calculations. We use a large grid of models ($10^{5}$ simulations) with and without planetary gaps, varying the initial conditions of the key parameters. We calculate the disk continuum emission and the effective radius for all models as a function of time. By selecting those simulations that continuously follow the SLR, we can derive constraints on the input parameters of the models. We confirm previous results that models of smooth disks in the radial drift regime are compatible with the observed SLR ($L_{mm}\propto r_{eff}^{2}$) but only smooth disks cannot be the reality. We show that the SLR is more widely populated if planets are present. However they tend to follow a different relation than smooth disks, potentially implying that a mixture of smooth and sub-structured disks are present in the observed sample. We derive a SLR ($L_{mm}\propto r_{eff}^{5/4}$) for disks with strong sub-structure. To be compatible with the SLR, models need to have an initially high disk mass ($\geq 2.5 \cdot 10^{-2}M_{\star}$) and low turbulence-parameter $\alpha$ values ($\leq 10^{-3}$). Furthermore, we find that the grain composition and porosity drastically affects the evolution of disks on the size-luminosity diagram where relatively compact grains that include amorphous carbon are favoured. Moreover, a uniformly optically thick disk with high albedo ($0.9$) that follows the SLR cannot be formed from an evolutionary procedure., Comment: Replaced with the published version in A&A. 19+11 pages, 11+10 figures

The mass of protoplanetary disks is arguably one of their most important quantities shaping their evolution toward planetary systems, but it remains a challenge to determine this quantity. Using the high spatial resolution now available on telescopes such as the Atacama Large Millimeter/submillimeter Array (ALMA), recent studies derived a relation between the disk surface density and the location of the "dust lines". This is a new concept in the field, linking the disk size at different continuum wavelengths with the radial distribution of grain populations of different sizes. We aim to use a dust evolution model to test the dependence of the dust line location on disk gas mass. In particular, we are interested in the reliability of the method for disks showing radial substructures, as recent high-resolution observations revealed. Our models show that the determination of the dust line location is a promising approach to the mass estimate of protoplanetay disks, but the exact relation between the dust line location and disk mass depends on the structure of the particular disk. We calibrated the relation for disks without evidence of radial structures, while for more complex structures we ran a simple dust evolution model. However, this method fails when there is evidence of strong dust traps. It is possible to reveal when dust evolution is dominated by traps, providing the necessary information for when the method should be applied with caution., 14 pages, 18 figures, accepted by Astronomy & Astrophysics

ALMA surveys have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass in exoplanets, under the assumption that the mm disk continuum emission is optically thin. This optically thin assumption seems to be supported by recent DSHARP observations where the measured optical depths of spatially resolved disks are mostly less than one. However, we point out that dust scattering can considerably reduce the emission from an optically thick region. If that scattering is ignored, the optical depth will be considerably underestimated. An optically thick disk with scattering can be misidentified as an optically thin disk. Dust scattering in more inclined disks can reduce the intensity even further, making the disk look even fainter. The measured optical depth of $\sim$0.6 in several DSHARP disks can be naturally explained by optically thick dust with an albedo of $\sim$0.9 at 1.25 mm. Using the DSHARP opacity, this albedo corresponds to a dust population with the maximum grain size ($s_{max}$) of 0.1-1 mm. For optically thick scattering disks, the measured spectral index $\alpha$ can be either larger or smaller than 2 depending on if the dust albedo increases or decreases with wavelength. Using the DSHARP opacity, $\alpha, Comment: 15 pages, 11 figures, Accepted for publication in ApJL

The Disk Substructures at High Angular Resolution Project (DSHARP) provides a large sample of protoplanetary disks having substructures which could be induced by young forming planets. To explore the properties of planets that may be responsible for these substructures, we systematically carry out a grid of 2-D hydrodynamical simulations including both gas and dust components. We present the resulting gas structures, including the relationship between the planet mass and 1) the gaseous gap depth/width, and 2) the sub/super-Keplerian motion across the gap. We then compute dust continuum intensity maps at the frequency of the DSHARP observations. We provide the relationship between the planet mass and 1) the depth/width of the gaps at millimeter intensity maps, 2) the gap edge ellipticity and asymmetry, and 3) the position of secondary gaps induced by the planet. With these relationships, we lay out the procedure to constrain the planet mass using gap properties, and study the potential planets in the DSHARP disks. We highlight the excellent agreement between observations and simulations for AS 209 and the detectability of the young Solar System analog. Finally, under the assumption that the detected gaps are induced by young planets, we characterize the young planet population in the planet mass-semimajor axis diagram. We find that the occurrence rate for $>$ 5 $M_J$ planets beyond 5-10 au is consistent with direct imaging constraints. Disk substructures allow us probe a wide-orbit planet population (Neptune to Jupiter mass planets beyond 10 au) that is not accessible to other planet searching techniques., 36 pages, 22 figures, ApJ Letters, in press

We present a detailed analysis of new ALMA observations of the disk around the T-Tauri star HD 143006, which at 46 mas (7.6 au) resolution reveal new substructures in the 1.25 mm continuum emission. The disk resolves into a series of concentric rings and gaps together with a bright arc exterior to the rings that resembles hydrodynamics simulations of a vortex, and a bridge-like feature connecting the two innermost rings. Although our $^{12}$CO observations at similar spatial resolution do not show obvious substructure, they reveal an inner disk depleted of CO emission. From the continuum emission and the CO velocity field we find that the innermost ring has a higher inclination than the outermost rings and the arc. This is evidence for either a small ($\sim8^{\circ}$) or moderate ($\sim41^{\circ}$) misalignment between the inner and outer disk, depending on the specific orientation of the near/far sides of the inner/outer disk. We compare the observed substructures in the ALMA observations with recent scattered light data from VLT/SPHERE of this object. In particular, the location of narrow shadow lanes in the SPHERE image combined with pressure scale height estimates, favor a large misalignment of about $41^{\circ}$. We discuss our findings in the context of a dust-trapping vortex, planet-carved gaps, and a misaligned inner disk due to the presence of an inclined companion to HD 143006., 15 pages, 8 figures, accepted for publication at ApJ Letters

The Disk Substructures at High Angular Resolution Project (DSHARP) is the largest homogeneous high-resolution ($\sim 0.035$ arcsec, or $\sim$ 5 au) disk continuum imaging survey with ALMA so far. In the coming years, many more disks will be mapped with ALMA at similar resolution. Interpreting the results in terms of the properties and quantities of the emitting dusty material is, however, a very non-trivial task. This is in part due to the uncertainty in the dust opacities, an uncertainty which is not likely to be resolved any time soon. It is also partly due to the fact that, as the DSHARP survey has shown, these disk often contain regions of intermediate to high optical depth, even at millimeter wavelengths and at relatively large radius in the disk. This makes the interpretation challenging, in particular if the grains are large and have a large albedo. On the other hand, the highly structured features seen in the DSHARP survey, of which strong indications were already seen in earlier observations, provide a unique opportunity to study the dust growth and dynamics. To provide continuity within the DSHARP project, its follow-up projects, and projects by other teams interested in these data, we present here the methods and opacity choices used within the DSHARP collaboration to link the measured intensity $I_\nu$ to dust surface density $\Sigma_d$., Comment: ApJ Letters, in press

We present a sub-arcsecond resolution survey of the 340 GHz dust continuum emission from 50 nearby protoplanetary disks, based on new and archival observations with the Submillimeter Array. The observed visibility data were modeled with a simple prescription for the radial surface brightness profile. The results were used to extract intuitive, empirical estimates of the emission "size" for each disk, $R_{\rm eff}$, defined as the radius that encircles a fixed fraction of the total continuum luminosity, $L_{\rm mm}$. We find a significant correlation between the sizes and luminosities, such that $R_{\rm eff} \propto L_{\rm mm}^{0.5}$, providing a confirmation and quantitative characterization of a putative trend that was noted previously. This correlation suggests that these disks have roughly the same average surface brightness interior to their given effective radius, ~0.2 Jy arcsec$^{-2}$ (or 8 K in brightness temperature). The same trend remains, but the 0.2dex of dispersion perpendicular to this relation essentially disappears, when we account for the irradiation environment of each disk with a crude approximation of the dust temperatures based on the stellar host luminosities. We consider two (not mutually exclusive) explanations for the origin of this size-luminosity relationship. Simple models of the growth and migration of disk solids can account for the observed trend for a reasonable range of initial conditions, but only on timescales that are much shorter than the nominal ages present in the sample. An alternative scenario invokes optically thick emission concentrated on unresolved scales, with filling factors of a few tens of percent, that are perhaps manifestations of localized particle traps., Accepted for publication in ApJ

We present the next step in a series of papers devoted to connecting the composition of the atmospheres of forming planets with the chemistry of their natal evolving protoplanetary disks. The model presented here computes the coupled chemical and dust evolution of the disk and the formation of three planets per disk model. Our three canonical planet traps produce a Jupiter near 1 AU, a Hot Jupiter and a Super-Earth. We study the dependency of the final orbital radius, mass, and atmospheric chemistry of planets forming in disk models with initial disk masses that vary by 0.02 $M_\odot$ above and below our fiducial model ($M_{disk,0} = 0.1 ~M_\odot$). We compute C/O and C/N for the atmospheres formed in our 3 models and find that C/O$_{\rm planet}\sim$ C/O$_{\rm disk}$, which does not vary strongly between different planets formed in our model. The nitrogen content of atmospheres can vary in planets that grow in different disk models. These differences are related to the formation history of the planet, the time and location that the planet accretes its atmosphere, and are encoded in the bulk abundance of NH$_3$. These results suggest that future observations of atmospheric NH$_3$ and an estimation of the planetary C/O and C/N can inform the formation history of particular planetary systems., Accepted for publication in MNRAS

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., Comment: Accepted to A&A. 16 pages, 11 figures

We present long-baseline Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 870 micron continuum emission from the nearest gas-rich protoplanetary disk, around TW Hya, that trace millimeter-sized particles down to spatial scales as small as 1 AU (20 mas). These data reveal a series of concentric ring-shaped substructures in the form of bright zones and narrow dark annuli (1-6 AU) with modest contrasts (5-30%). We associate these features with concentrations of solids that have had their inward radial drift slowed or stopped, presumably at local gas pressure maxima. No significant non-axisymmetric structures are detected. Some of the observed features occur near temperatures that may be associated with the condensation fronts of major volatile species, but the relatively small brightness contrasts may also be a consequence of magnetized disk evolution (the so-called zonal flows). Other features, particularly a narrow dark annulus located only 1 AU from the star, could indicate interactions between the disk and young planets. These data signal that ordered substructures on ~AU scales can be common, fundamental factors in disk evolution, and that high resolution microwave imaging can help characterize them during the epoch of planet formation., published in Astrophysical Journal Letters; image files are available at http://www.cfa.harvard.edu/~sandrews/data/twhya/

We combined hydrodynamical simulations of planet-disk interactions with dust evolution models that include coagulation and fragmentation of dust grains over a large range of radii and derived observational properties using radiative transfer calculations. We studied the role of the snow line in the survival of the inner disk of transition disks. Inside the snow line, the lack of ice mantles in dust particles decreases the sticking efficiency between grains. As a consequence, particles fragment at lower collision velocities than in regions beyond the snow line. This effect allows small particles to be maintained for up to a few Myrs within the first astronomical unit. These particles are closely coupled to the gas and do not drift significantly with respect to the gas. For lower mass planets (1$M_{\rm{Jup}}$), the pre-transition appearance can be maintained even longer because dust still trickles through the gap created by the planet, moves invisibly and quickly in the form of relatively large grains through the gap, and becomes visible again as it fragments and gets slowed down inside of the snow line. The global study of dust evolution of a disk with an embedded planet, including the changes of the dust aerodynamics near the snow line, can explain the concentration of millimetre-sized particles in the outer disk and the survival of the dust in the inner disk if a large dust trap is present in the outer disk. This behaviour solves the conundrum of the combination of both near-infrared excess and ring-like millimetre emission observed in several transition disks., Accepted for publication in A&A (including acknowledgments)

Molecular snow lines in protoplanetary disks have been studied theoretically for decades because of their importance in shaping planetary architectures and compositions. The water snow line lies in the planet formation region at < 10 AU, and so far its location has been estimated only indirectly from spatially-unresolved spectroscopy. This work presents a proof-of-concept method to directly image the water snow line in protoplanetary disks through its physical and chemical imprint in the local dust properties. We adopt a physical disk model that includes dust coagulation, fragmentation, drift, and a change in fragmentation velocities of a factor 10 between dry silicates and icy grains as found by laboratory work. We find that the presence of a water snow line leads to a sharp discontinuity in the radial profile of the dust emission spectral index {\alpha}_mm, due to replenishment of small grains through fragmentation. We use the ALMA simulator to demonstrate that this effect can be observed in protoplanetary disks using spatially-resolved ALMA images in two continuum bands. We explore the model dependence on the disk viscosity and find that the spectral index reveals the water snow line for a wide range of conditions, with opposite trends when the emission is optically thin rather than thick. If the disk viscosity is low ({\alpha}_visc < 10^-3) the snow line produces a ring-like structure with a minimum at {\alpha}_mm ~ 2 in the optically thick regime, possibly similar to what has been measured with ALMA in the innermost region of the HL Tau disk., Comment: 6 pages, 4 figures, accepted for publication in ApJ Letters

Recent imaging of protoplanetary disks with high resolution and contrast have revealed a striking variety of substructure. Of particular interest are cases where near-infrared scattered light images show evidence for low-intensity annular "gaps." The origins of such structures are still uncertain, but the interaction of the gas disk with planets is a common interpretation. We study the impact that the evolution of the solid material can have on the observable properties of disks in a simple scenario without any gravitational or hydrodynamical disturbances to the gas disk structure. Even with a smooth and continuous gas density profile, we find that the scattered light emission produced by small dust grains can exhibit ring-like depressions similar to those presented in recent observations. The physical mechanisms responsible for these features rely on the inefficient fragmentation of dust particles. The occurrence and position of the proposed "gap" features depend most strongly on the dust-to-gas ratio, the fragmentation threshold velocity, the strength of the turbulence, and the age of the disk, and should be generic (at some radius) for typically adopted disk parameters. The same physical processes can affect the thermal emission at optically thin wavelengths ($\sim$1 mm), although the behavior can be more complex; unlike for disk-planet interactions, a "gap" should not be present at these longer wavelengths., 6 pages, 4 figures, accepted for publication in ApJL, fixed typos

Astronomy and Astrophysics 573 (2015): A19 reproduced with permission from Astronomy & Astrophysics, Context. The Cep OB2 region, with its two intermediate-aged clusters Tr 37 and NGC7160, is a paradigm of sequential star formation and an ideal site for studies of protoplanetary disk evolution. Aims. We use Herschel data to study the protoplanetary disks and the star formation history of the region. Methods. Herschel/PACS observations at 70 and 160 μm probe the disk properties (mass, dust sizes, structure) and the evolutionary state of a large number of young stars. Far-IR data also trace the remnant cloud material and small-scale cloud structure. Results. We detect 95 protoplanetary disks at 70 μm, 41 at 160 μm, and obtain upper limits for more than 130 objects. The detection fraction at 70 μm depends on the spectral type (88% for K4 or earlier stars, 17% for M3 or later stars) and on the disk type (∼50% for full and pre-transitional disks, ∼35% for transitional disks, no low-excess/depleted disks detected). Non-accreting disks are not detected, suggesting significantly lower masses. Accreting transition and pre-transition disks have systematically higher 70 μm excesses than full disks, suggestive of more massive, flared and/or thicker disks. Herschel data also reveal several mini-clusters in Tr 37, which are small, compact structures containing a few young stars surrounded by nebulosity. Conclusions. Far-IR data are an excellent probe of the evolution of disks that are too faint for sub-millimetre observations. We find a strong link between far-IR emission and accretion, and between the inner and outer disk structure. Herschel confirms the dichotomy between accreting and non-accreting transition disks. Accretion is a powerful measure of global disk evolution: substantial mass depletion and global evolution need to occur to shut down accretion in a protoplanetary disk, even if the disk has inner holes. Disks likely follow different evolutionary paths: low disk masses do not imply opening inner holes, and having inner holes does not require low disk masses. The mini-clusters reveal multi-episodic star formation in Tr 37. The long survival of mini-clusters suggest that they formed from the fragmentation of the same core. Their various morphologies favour different formation/triggering mechanisms acting within the same cluster. The beads-on-a-string structure in one mini-cluster is consistent with gravitational fragmentation or gravitational focusing, acting on very small scales (solar-mass stars in ∼0.5 pc filaments). Multi-episodic star formation could also produce evolutionary variations between disks in the same region. Finally, Herschel also unveils what could be the first heavy mass loss episode of the O6.5 star HD206267 in Tr 37, A.S.A. acknowledges support by the Spanish MICINN/MINECO “Ramón y Cajal” program, grant number RYC-2010-06164. A.S.A. and M.F. acknowledge support by the action “Proyectos de Investigación fundamental no orientada”, grant number AYA2012-35008. C.E. is partly supported by Spanish MICINN/MINECO grant AYA2011-26202. V.R. is supported by the DLR grant number 50 OR 1109 and by the Bayerischen Gleichstellungsförderung (BGF). T.B. acknowledges support from NASA Origins of Solar Systems grant NNX12AJ04G. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France

We present sub-arcsecond resolution observations of continuum emission associated with the GG Tau quadruple star system at wavelengths of 1.3, 2.8, 7.3, and 50 mm. These data confirm that the GG Tau A binary is encircled by a circumbinary ring at a radius of 235 AU with a FWHM width of ~60 AU. We find no clear evidence for a radial gradient in the spectral shape of the ring, suggesting that the particle size distribution is spatially homogeneous on angular scales of ~0.1". A central point source, likely associated with the primary component (GG Tau Aa), exhibits a composite spectrum from dust and free-free emission. Faint emission at 7.3 mm is observed toward the low-mass star GG Tau Ba, although its origin remains uncertain. Using these measurements of the resolved, multifrequency emission structure of the GG Tau A system, models of the far-infrared to radio spectrum are developed to place constraints on the grain size distribution and dust mass in the circumbinary ring. The non-negligible curvature present in the ring spectrum implies a maximum particle size of 1-10 mm, although we are unable to place strong constraints on the distribution shape. The corresponding dust mass is 30-300 M_earth, at a temperature of 20-30 K. We discuss how this significant concentration of relatively large particles in a narrow ring at a large radius might be produced in a local region of higher gas pressures (i.e., a particle "trap") located near the inner edge of the circumbinary disk., ApJ, in press

The expectation that aerodynamic drag will force the solids in a gas-rich protoplanetary disk to spiral in toward the host star on short timescales is one of the fundamental problems in planet formation theory. The nominal efficiency of this radial drift process is in conflict with observations, suggesting that an empirical calibration of solid transport mechanisms in a disk is highly desirable. However, the fact that both radial drift and grain growth produce a similar particle size segregation in a disk (such that larger particles are preferentially concentrated closer to the star) makes it difficult to disentangle a clear signature of drift alone. We highlight a new approach, by showing that radial drift leaves a distinctive "fingerprint" in the dust surface density profile that is directly accessible to current observational facilities. Using an analytical framework for dust evolution, we demonstrate that the combined effects of drift and (viscous) gas drag naturally produce a sharp outer edge in the dust distribution (or, equivalently, a sharp decrease in the dust-to-gas mass ratio). This edge feature forms during the earliest phase in the evolution of disk solids, before grain growth in the outer disk has made much progress, and is preserved over longer timescales when both growth and transport effects are more substantial. The key features of these analytical models are reproduced in detailed numerical simulations, and are qualitatively consistent with recent millimeter-wave observations that find gas/dust size discrepancies and steep declines in dust continuum emission in the outer regions of protoplanetary disks., 12 pages, 5 figures, accepted to ApJ

Context. Particle trapping in local or global pressure maxima in protoplanetary disks is one of the new paradigms in the theory of the first stages of planet formation. However, finding observational evidence for this effect is not easy. Recent work suggests that the large ring-shaped outer disks observed in transition disk sources may in fact be lopsided and constitute large banana-shaped vortices. Aims. We wish to investigate how effective dust can accumulate along the azimuthal direction. We also want to find out if the size- sorting resulting from this can produce a detectable signatures at millimeter wavelengths. Methods. To keep the numerical cost under control we develop a 1+1D method in which the azimuthal variations are treated sepa- rately from the radial ones. The azimuthal structure is calculated analytically for a steady-state between mixing and azimuthal drift. We derive equilibration time scales and compare the analytical solutions to time-dependent numerical simulations. Results. We find that weak, but long-lived azimuthal density gradients in the gas can induce very strong azimuthal accumulations of dust. The strength of the accumulations depends on the P\'eclet number, which is the relative importance of advection and diffusion. We apply our model to transition disks and our simulated observations show that this effect would be easily observable with ALMA and in principle allows to put constraints on the strength of turbulence and the local gas density., Comment: 4 pages, 4 figures, accepted for publication in A&A Letters

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., 4 pages, 3 figures. Accepted for publication as a Letter in Astronomy and Astrophysics

In order to explain grain growth to mm sized particles and their retention in outer regions of protoplanetary disks, as it is observed at sub-mm and mm wavelengths, we investigate if strong inhomogeneities in the gas density profiles can slow down excessive radial drift and can help dust particles to grow. We use coagulation/fragmentation and disk-structure models, to simulate the evolution of dust in a bumpy surface density profile which we mimic with a sinusoidal disturbance. For different values of the amplitude and length scale of the bumps, we investigate the ability of this model to produce and retain large particles on million years time scales. In addition, we introduced a comparison between the pressure inhomogeneities considered in this work and the pressure profiles that come from magnetorotational instability. Using the Common Astronomy Software Applications ALMA simulator, we study if there are observational signatures of these pressure inhomogeneities that can be seen with ALMA. We present the favorable conditions to trap dust particles and the corresponding calculations predicting the spectral slope in the mm-wavelength range, to compare with current observations. Finally we present simulated images using different antenna configurations of ALMA at different frequencies, to show that the ring structures will be detectable at the distances of the Taurus Auriga or Ophiucus star forming regions., Pages 15, Figures 14. Accepted for publication in Astronomy and Astrophysics

The streaming instability is often invoked as solution to the fragmentation and drift barriers in planetesimal formation, catalyzing the aggregation of dust on kyr timescales to grow km-sized cores. However there remains a lack of consensus on the physical mechanism(s) responsible for initiating it. One potential avenue is disc photoevaporation, wherein the preferential removal of relatively dust-free gas increases the disc metallicity. Late in the disc lifetime, photoevaporation dominates viscous accretion, creating a gradient in the depleted gas surface density near the location of the gap. This induces a local pressure maximum that collects drifting dust particles, which may then become susceptible to the streaming instability. Using a one-dimensional viscous evolution model of a disc subject to internal X-ray photoevaporation, we explore the efficacy of this process to build planetestimals. Over a range of parameters we find that the amount of dust mass converted into planetesimals is often < 1 M_Earth and at most a few M_Earth spread across tens of AU. We conclude that photoevaporation may at best be relevant for the formation of debris discs, rather than a common mechanism for the formation of planetary cores. Our results are in contrast to a recent, similar investigation that considered an FUV-driven photoevaporation model and reported the formation of tens of M_Earth at large (> 100 AU) disc radii. The discrepancies are primarily a consequence of the different photoevaporation profiles assumed. Until observations more tightly constrain photoevaporation models, the relevance of this process to the formation of planets remains uncertain., 10 pages, 5 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

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., Comment: Accepted for publication in ApJL, 10 pages, 5 figures

We introduce the Disk Substructures at High Angular Resolution Project (DSHARP), one of the initial Large Programs conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal of DSHARP is to find and characterize substructures in the spatial distributions of solid particles for a sample of 20 nearby protoplanetary disks, using very high resolution (0.035 arcsec, or 5 au FWHM) observations of their 240 GHz (1.25 mm) continuum emission. These data provide a first homogeneous look at the small-scale features in disks that are directly relevant to the planet formation process, quantifying their prevalence, morphologies, spatial scales, spacings, symmetry, and amplitudes, for targets with a variety of disk and stellar host properties. We find that these substructures are ubiquitous in this sample of large, bright disks. They are most frequently manifested as concentric, narrow emission rings and depleted gaps, although large-scale spiral patterns and small arc-shaped azimuthal asymmetries are also present in some cases. These substructures are found at a wide range of disk radii (from a few au to more than 100 au), are usually compact ($, Comment: ApJ Letters, in press; some figure resolutions reduced, full-res version at https://www.cfa.harvard.edu/~sandrews/DSHARP_I.pdf ; Data Release at https://almascience.org/alma-data/lp/DSHARP

The Disk Substructures at High Angular Resolution Project used ALMA to map the 1.25 millimeter continuum of protoplanetary disks at a spatial resolution of ~5 au. We present a systematic analysis of annular substructures in the 18 single-disk systems targeted in this survey. No dominant architecture emerges from this sample; instead, remarkably diverse morphologies are observed. Annular substructures can occur at virtually any radius where millimeter continuum emission is detected and range in widths from a few au to tens of au. Intensity ratios between gaps and adjacent rings range from near-unity to just a few percent. In a minority of cases, annular substructures co-exist with other types of substructures, including spiral arms (3/18) and crescent-like azimuthal asymmetries (2/18). No clear trend is observed between the positions of the substructures and stellar host properties. In particular, the absence of an obvious association with stellar host luminosity (and hence the disk thermal structure) suggests that substructures do not occur preferentially near major molecular snowlines. Annular substructures like those observed in DSHARP have long been hypothesized to be due to planet-disk interactions. A few disks exhibit characteristics particularly suggestive of this scenario, including substructures in possible mean-motion resonance and "double gap" features reminiscent of hydrodynamical simulations of multiple gaps opened by a planet in a low-viscosity disk., 33 pages, 12 figures, captions corrected

We present the first kinematical detection of embedded protoplanets within a protoplanetary disk. Using archival ALMA observations of HD 163296, we demonstrate a new technique to measure the rotation curves of CO isotopologue emission to sub-percent precision relative to the Keplerian rotation. These rotation curves betray substantial deviations caused by local perturbations in the radial pressure gradient, likely driven by gaps carved in the gas surface density by Jupiter-mass planets. Comparison with hydrodynamic simulations shows excellent agreement with the gas rotation profile when the disk surface density is perturbed by two Jupiter mass planets at 83 au and 137 au. As the rotation of the gas is dependent on the pressure of the total gas component, this method provides a unique probe of the gas surface density profile without incurring significant uncertainties due to gas-to-dust ratios or local chemical abundances which plague other methods. Future analyses combining both methods promise to provide the most accurate and robust measures of embedded planetary mass. Furthermore, this method provides a unique opportunity to explore wide-separation planets beyond the mm continuum edge and to trace the gas pressure profile essential in modelling grain evolution in disks., Comment: Accepted by ApJL