The Disk Substructures at High Angular Resolution Project (DSHARP) is the largest homogeneous high-resolution (∼0 035, or ∼5 au) disk continuum imaging survey with the Atacama Large ...Millimeter/submillimeter Array (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 that 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 to dust surface density d.
The Disk Substructures at High Angular Resolution Project (DSHARP) provides a large sample of protoplanetary disks with substructures that 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 2D hydrodynamical simulations, including both gas and dust components. We present the resulting gas structures, including the relationship between the planet mass, as well as (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, as well as (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 MJ planets beyond 5-10 au is consistent with direct imaging constraints. Disk substructures allow us to probe a wide-orbit planet population (Neptune to Jupiter mass planets beyond 10 au) that is not accessible to other planet searching techniques.
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, 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 astronomical units to more than 100 au), are usually compact ( 10 au), and show a wide range of amplitudes (brightness contrasts). Here we discuss the motivation for the project, describe the survey design and the sample properties, detail the observations and data calibration, highlight some basic results, and provide a general overview of the key conclusions that are presented in more detail in a series of accompanying articles. The DSHARP data-including visibilities, images, calibration scripts, and more-are released for community use at https://almascience.org/alma-data/lp/DSHARP.
Abstract Turbulence is crucial for protoplanetary disk dynamics, and vertical shear instability (VSI) is a promising mechanism in outer disk regions to generate turbulence. We use the Athena++ ...radiation module to study VSI in full and transition disks, accounting for radiation transport and stellar irradiation. We find that the thermal structure and cooling timescale significantly influence VSI behavior. The inner rim location and radial optical depth affect disk kinematics. Compared with previous vertically isothermal simulations, our full disk and transition disks with small cavities have a superheated atmosphere and cool midplane with long cooling timescales, which suppresses the corrugation mode and the associated meridional circulation. This temperature structure also produces a strong vertical shear at τ * = 1, producing an outgoing flow layer at τ * < 1 on top of an ingoing flow layer at τ * ∼ 1. The midplane becomes less turbulent, while the surface becomes more turbulent with effective α reaching ∼10 −2 at τ * ≲ 1. This large surface stress drives significant surface accretion, producing substructures. Using temperature and cooling time measured/estimated from radiation-hydro simulations, we demonstrate that less computationally intensive simulations incorporating simple orbital cooling can almost reproduce radiation-hydro results. By generating synthetic images, we find that substructures are more pronounced in disks with larger cavities. The higher velocity dispersion at the gap edge could also slow particle settling. Both properties are consistent with recent near-IR and Atacama Large Millimeter/submillimeter Array (ALMA) observations. Our simulations predict that regions with significant temperature changes are accompanied by significant velocity changes, which can be tested by ALMA kinematics/chemistry observations.
A large fraction of the protoplanetary disks observed with ALMA display multiple well-defined and nearly perfectly circular rings in the continuum, in many cases with substantial peak-to-valley ...contrast. The DSHARP campaign shows that several of these rings are very narrow in radial extent. In this Letter we test the hypothesis that these dust rings are caused by dust trapping in radial pressure bumps, and if confirmed, put constraints on the physics of the dust trapping mechanism. We model this process analytically in 1D, assuming axisymmetry. By comparing this model to the data, we find that all rings are consistent with dust trapping. Based on a plausible model of the dust temperature we find that several rings are narrower than the pressure scale height, providing strong evidence for dust trapping. The rings have peak absorption optical depth in the range between 0.2 and 0.5. The dust masses stored in each of these rings is of the order of tens of Earth masses, though much ambiguity remains due to the uncertainty of the dust opacities. The dust rings are dense enough to potentially trigger the streaming instability, but our analysis cannot give proof of this mechanism actually operating. Our results show, however, that the combination of very low and very large grains can be excluded by the data for all the rings studied in this Letter.
The Disk Substructures at High Angular Resolution Project (DSHARP) used ALMA to map the 1.25 mm 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 astronomical units to tens of astronomical units. Intensity ratios between gaps and adjacent rings range from near-unity to just a few percent. In a minority of cases, annular substructures coexist 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.
Abstract
The polarimetric observations of the protoplanetary disk around HL Tau have shown the scattering-induced polarization at ALMA Band 7, which indicates that the maximum dust size is ∼100
μ
m, ...while the spectral energy distribution (SED) has suggested that the maximum dust size is approximately a millimeter. To solve the contradiction, we investigate the impact of differential settling of dust grains on the SED and polarization. If the disk is optically thick, a longer observing wavelength traces more interior layers, which would be dominated by larger grains. We find that the SED of the center part of the HL Tau disk can be explained with millimeter-sized grains for a broad range of turbulence strength, while 160
μ
m–sized grains cannot be explained unless the turbulence strength parameter
α
t
is lower than 10
−5
. We also find that the observed polarization fraction can be potentially explained with a maximum dust size of 1 mm if
α
t
≲ 10
−5
, although models with 160
μ
m–sized grains are also acceptable. However, if the maximum dust size is ∼3 mm, the simulated polarization fraction is too low to explain the observations even if the turbulence strength is extremely small, indicating a maximum dust size of ≲1 mm. The degeneracy between 100
μ
m– and millimeter-sized grains can be solved by improving the ALMA calibration accuracy or polarimetric observations at (sub)centimeter wavelengths.
We present an analysis of Atacama Large Millimeter/submillimeter Array 1.25 mm continuum observations of spiral structures in three protoplanetary disks from the Disk Substructures at High Angular ...Resolution Project. The disks around Elias 27, IM Lup, and WaOph 6 were observed at a resolution of ∼40-60 mas (∼6-7 au). All three disks feature m = 2 spiral patterns in conjunction with annular substructures. Gas kinematics established by 12CO J = 2−1 observations indicate that the continuum spiral arms are trailing. The arm-interarm intensity contrasts are modest, typically less than 3. The Elias 27 spiral pattern extends throughout much of the disk, and the arms intersect the gap at R ∼ 69 au. The spiral pattern in the IM Lup disk is particularly complex-it extends about halfway radially through the disk, exhibiting pitch angle variations with radius and interarm features that may be part of ring substructures or spiral arm branches. Spiral arms also extend most of the way through the WaOph 6 disk, but the source overall is much more compact than the other two disks. We discuss possible origins for the spiral structures, including gravitational instability and density waves induced by a stellar or planetary companion. Unlike the millimeter continuum counterparts of many of the disks with spiral arms detected in scattered light, these three sources do not feature high-contrast crescent-like asymmetries or large (R > 20 au) emission cavities. This difference may point to multiple spiral formation mechanisms operating in disks.
The Atacama Large Millimeter/submillimeter Array observations of protoplanetary disks acquired by the Disk Substructure at High Angular Resolution Project resolve the dust and gas emission on angular ...scales as small as 3 astronomical units, offering an unprecedented detailed view of the environment where planets form. In this Letter, we present and discuss observations of the HD 163296 protoplanetary disk that imaged the 1.25 mm dust continuum and 12CO J = 2-1 rotational line emission at a spatial resolution of 4 and 10 au, respectively. The continuum observations resolve and allow us to characterize the previously discovered dust rings at radii of 68 and 100. They also reveal new small-scale structures, such as a dark gap at 10 au, a bright ring at 15 au, a dust crescent at a radius of 55 au, and several fainter azimuthal asymmetries. The observations of the CO and dust emission provide information about the vertical structure of the disk and allow us to directly constrain the dust extinction optical depth at the dust rings. Furthermore, the observed asymmetries in the dust continuum emission corroborate the hypothesis that the complex structure of the HD 163296 disk is the result of the gravitational interaction with yet-unseen planets.
The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence ...on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ∼50 mas (∼8 au), as well as HCO+ J = 3 − 2 observations at a resolution of ∼100 mas. The dust continuum shows at least three rings atop faint, extended emission. Unresolved emission is detected at the center of the disk cavity at both wavelengths, likely due to a combination of dust and free-free emission. Compared to the 1.1 mm image, the 2.1 mm image shows a more pronounced "shoulder" near R ∼ 40 au, highlighting the utility of longer-wavelength observations for characterizing disk substructures. The spectral index features strong radial variations, with minima near the emission peaks and maxima near the gaps. While low spectral indices have often been ascribed to grain growth and dust trapping, the optical depth of GM Aur's inner two emission rings renders their dust properties ambiguous. The gaps and outer disk (R > 100 au) are optically thin at both wavelengths. Meanwhile, the HCO+ emission indicates that the gas cavity is more compact than the dust cavity traced by the millimeter continuum, similar to other disks traditionally classified as "transitional."