# The Disk Substructures at High Angular Resolution Project (DSHARP). I. Motivation, Sample, Calibration, and Overview

Sean M. Andrews, Jane Huang, Laura M. Pérez, Andrea Isella, Cornelis P.

Dullemond, Nicolás T. Kurtovic, Viviana V. Guzmán, John M. Carpenter, David J. Wilner, Shangjia Zhang, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Sean M. Andrews, Jane Huang, Laura M. Pérez, Andrea Isella, Cornelis P. Dullemond, et al. The Disk Substructures at High Angular Resolution Project (DSHARP). I. Motivation, Sample, Calibration, and Overview. The Astrophysical Journal Letters, 2018, 869, 10.3847/2041-8213/aaf741. insu-03635024

HAL Id: insu-03635024

https://insu.hal.science/insu-03635024v1

Submitted on 1 Sep 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# The Disk Substructures at High Angular Resolution Project (DSHARP).
## I. Motivation, Sample, Calibration, and Overview

Sean M. Andrews¹, Jane Huang¹, Laura M. Pérez², Andrea Isella³, Cornelis P. Dullemond⁴, Nicolás T. Kurtovic²,
Viviana V. Guzmán⁵,⁶, John M. Carpenter⁵, David J. Wilner¹, Shangjia Zhang⁷, Zhaohuan Zhu⁷, Tilman Birnstiel⁸,¹
Xue-Ning Bai⁹, Myriam Benisty¹⁰,¹¹, A. Meredith Hughes¹², Karin I. Öberg¹, and Luca Ricci¹³

¹Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA

²Departamento de Astronomía, Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile

³Department of Physics and Astronomy, Rice University, 6100 Main Street, Houston, TX 77005, USA

⁴Zentrum für Astronomie, Heidelberg University, Albert Ueberle Str. 2, D-69120 Heidelberg, Germany

⁵Joint ALMA Observatory, Avenida Alonso de Córdova 3107, Vitacura, Santiago, Chile

⁶Instituto de Astrofísica, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile

⁷Department of Physics and Astronomy, University of Nevada, Las Vegas, 4505 S. Maryland Parkway, Las Vegas, NV 89154, USA

⁸University Observatory, Faculty of Physics, Ludwig-Maximilians-Universität München, Scheinerstr. 1, D-81679 Munich, Germany

⁹Institute for Advanced Study and Tsinghua Center for Astrophysics, Tsinghua University, Beijing 100084, People's Republic of China

¹⁰Unidad Mixta Internacional Franco-Chilena de Astronomía, CNRS/INSU UMI 3386, Departamento de Astronomía,
Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile

¹¹Univ. Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France

¹²Department of Astronomy, Van Vleck Observatory, Wesleyan University, 96 Foss Hill Drive, Middletown, CT 06459, USA

¹³Department of Physics and Astronomy, California State University Northridge, 18111 Nordhoff Street, Northridge, CA 91130, USA

Received 2018 September 26; revised 2018 November 13; accepted 2018 November 21; published 2018 December 26

## Abstract

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.

**Key words:** circumstellar matter – planets and satellites: formation – protoplanetary disks

**Supporting material:** figure set, machine-readable tables

## 1. Introduction and Motivation

There is a long-standing desire to link the properties of circumstellar disks with the initial conditions of planetary systems. The theoretical aspiration in the field is to develop a deterministic framework that takes a set of measured disk properties (e.g., the spatial distribution of densities and temperatures; Andrews et al. 2009, 2010; Isella et al. 2009, 2010) and predicts the key characteristics of the exoplanet population (e.g., masses, orbital architectures, and atmosphere compositions; Ida & Lin 2004, 2008; Alibert et al. 2005; Mordasini et al. 2009). A quality reproduction in this population synthesis context requires the tuning of increasingly sophisticated models for the formation of planetary systems, their interactions with disk material, and their subsequent long-term dynamical evolution (see Benz et al. 2014).

The most crucial obstacle in the planet formation process is the assembly of planetesimals (see Johansen et al. 2014). The formation of terrestrial planets and giant planet cores hinges on the rapid agglomeration of small particles into these much

larger (≥kilometer-sized) bodies. Astronomers have worked on this topic and its pitfalls for more than 50 yr, although without much observational guidance. Fortunately, that is changing. Resolved observations of the continuum emission from millimeter/centimeter-sized particles in disks measure how the solids are distributed. Resolved variations in the continuum spectrum shape have been interpreted as radial gradients in the particle size distributions (larger solids closer to the star; Isella et al. 2010; Guilloteau et al. 2011; Pérez et al. 2012, 2015; Menu et al. 2014; Tazzari et al. 2016; Tripathi et al. 2018). Pronounced discrepancies between the spatial distributions of continuum and spectral line emission have led to suggestions that the mass ratio of solids relative to gas also varies with radius (higher closer to the star; Panić et al. 2009; Andrews et al. 2012; de Gregorio-Monsalvo et al. 2013; Rosenfeld et al. 2013; Zhang et al. 2014; Facchini et al. 2017; Ansdell et al. 2018). Those results provide strong qualitative support for evolutionary models of solids early in the
planetesimal assembly process (e.g., Birnstiel & Andrews 2014; Testi et al. 2014; Birnstiel et al. 2016).

Despite that progress, there is still considerable tension regarding planetesimal formation timescales for the default assumption of a smooth gas disk (with pressure, *P*, decreasing monotonically with radius, *r*). This tension is associated with radial drift, the inward migration of solids toward the global *P* maximum that occurs when they decouple from the sub-Keplerian gas flow (Adachi et al. 1976; Weidenschilling 1977; Nakagawa et al. 1986). The predicted drift rates for millimeter/centimeter solids located tens of astronomical units from the host star are fast enough to severely limit planetesimal growth (Takeuchi & Lin 2002, 2005; Brauer et al. 2007, 2008) and are in conflict with routine observations of emission from those particles at *r* ≈ 10–100 au (e.g., Barenfeld et al. 2017; Tazzari et al. 2017; Tripathi et al. 2017; Andrews et al. 2018).

This contradiction indicates that the *P(r)* profiles in disks are likely not smooth. Localized *P* modulations can slow or trap drifting solids (Whipple 1972; Pinilla et al. 2012a), perhaps concentrating them enough to trigger gravitational and/or streaming instabilities that rapidly convert pebbles to planetesimals (e.g., Youdin & Shu 2002; Youdin & Goodman 2005; Johansen et al. 2009). Such particle traps or other migration bottlenecks could be produced by the dynamics associated with how gas, dust, and magnetic fields are coupled (e.g., Dzyurkevich et al. 2013; Bai & Stone 2014; Dipierro et al. 2015; Flock et al. 2015; Lyra et al. 2015; Béthune et al. 2017; Dullemond & Penzlin 2018; Suriano et al. 2018) or by strong gradients in material properties (e.g., Okuzumi et al. 2012; Armitage et al. 2016; Estrada et al. 2016; Pinilla et al. 2017; Stammler et al. 2017; but see Long et al. 2018; van Terwisga et al. 2018). These small-scale material concentrations—substructures—are largely absent in contemporary models of planet formation, but they would likely play fundamental roles in nearly all aspects of the formation process.

If such substructures were prominent in disks at early evolution stages, it is possible that planetesimals and even entire planetary systems were created much more efficiently than is expected in the traditional models (e.g., Greaves & Rice 2010; Najita & Kenyon 2014; Nixon et al. 2018). In that scenario, the typical approximately million-year-old disk may harbor a “second generation” of substructures created by the dynamical interactions between young planets and their nascent disk material (see Lin & Papaloizou 1993; Kley & Nelson 2012), which in turn can affect the orbital architectures of those burgeoning planetary systems (e.g., Coleman & Nelson 2016).

In any case, observations of disk substructures are essential. Direct constraints on small-scale gas pressure variations in disks based on high resolution measurements of molecular line emission are a formidable challenge. However, the particle trapping capabilities of even modest pressure maxima should substantially amplify the associated local millimeter/centimeter-sized particle density (e.g., Paardekooper & Mellemela 2006; Rice et al. 2006; Pinilla et al. 2012b; Zhu et al. 2012), generating a bright signature in the broadband (sub) millimeter continuum that is much easier to measure on the smallest scales.

The initial foray into such work came from the “transition” disks (Strom et al. 1989; Skrutskie et al. 1990; Calvet et al. 2002), which show dense particle rings at *r* ≈ tens of astronomical units, outside depleted central cavities (e.g., Andrews et al. 2011;

Pinilla et al. 2018; van der Marel et al. 2018). Observations with sufficient resolution reveal that these particle traps exhibit complex substructures, including azimuthal asymmetries (Casassus et al. 2013; Isella et al. 2013; van der Marel et al. 2013; Pérez et al. 2014), additional rings (Fedele et al. 2017; van der Plas et al. 2017), warped geometries (and/or radial inflows; Rosenfeld et al. 2012, 2014; Marino et al. 2015; Casassus et al. 2018), and spiral arms (Christiaens et al. 2014; Boehler et al. 2018; Dong et al. 2018). Similar features have been identified from the IR starlight scattered off the disk atmospheres (e.g., Muto et al. 2012; Grady et al. 2013; Quanz et al. 2013; Avenhaus et al. 2014; Benisty et al. 2015; Rapson et al. 2015; Akiyama et al. 2016; de Boer et al. 2016; Ginski et al. 2016).

Some serendipitous discoveries at modest (~10–20 au) resolution hint that the more general disk population frequently exhibits substructures in the forms of rings/gaps (Cieza et al. 2016, 2017; Isella et al. 2016; Zhang et al. 2016; Cox et al. 2017; Huang et al. 2017; Loomis et al. 2017; Dipierro et al. 2018; Fedele et al. 2018; van Terwisga et al. 2018) and spirals (Pérez et al. 2016). The richness of these substructures becomes clear for the few individual cases that have had their continuum emission probed at resolutions of only a few astronomical units (HL Tau, ALMA Partnership et al. 2015; TW Hya, Andrews et al. 2016; MWC 758, Dong et al. 2018). Again, similar conclusions are being drawn from complementary measurements of scattered light from small dust grains (e.g., van Boekel et al. 2017; Avenhaus et al. 2018).

All of these observations suggest that substructures are common, and therefore are likely significant factors in many disk evolution and planet formation processes. Moreover, they demonstrate a tremendous opportunity: high resolution millimeter continuum measurements can quantify the forms, prevalence, and diversity (e.g., in scales, locations, and amplitudes) of disk substructures, and thereby help develop a more robust theoretical framework for characterizing the early evolution of planetary systems. The next step along that path is to move from a serendipitous discovery-space to a principled survey specifically designed to study these features.

In this article, we introduce a new survey that moves in this direction. The Disk Substructures at High Angular Resolution Project (DSHARP) was conducted as one of the first ALMA Large Programs. DSHARP measures the 240 GHz continuum emission at ~35 mas (5 au) resolution for 20 disks, to help better understand the evolution of solid particles during the planet formation process. Having motivated the project, this article also describes the DSHARP survey design and sample (Section 2), the ALMA observations (Section 3) and their calibration (Section 4), along with some basic observational results and the DSHARP data release (Section 5). We conclude with an overview of the highlights from a series of accompanying articles (Section 6).

## 2. Survey Design and Sample

The DSHARP survey was designed to optimize the spatial resolution and contrast sensitivity to continuum emission substructures. Secondarily, measurements of CO line emission were also of interest as a preliminary opportunity to identify corresponding gas structures and infer other relevant bulk disk properties (e.g., geometry). We defined two criteria to guide the survey design, based on previous observations and theoretical expectations for the origins of disk substructures.
The first criterion was access to a wide range of spatial scales down to an FWHM resolution of ~5 au. Such high resolution was essential for identifying the disk substructures in the sharpest ALMA continuum images available to date (ALMA Partnership et al. 2015; Andrews et al. 2016). Moreover, it is comparable to the (disk-averaged) pressure scale height, $h_p$ (where $h_p/r \approx 0.1$; Kenyon & Hartmann 1987), a benchmark size that is directly related to the *P* deviations generated by turbulent zonal flows (e.g., Johansen et al. 2009), vortices (e.g., Barge & Sommeria 1995), or planetary gaps (e.g., Bryden et al. 1999). At 5 au resolution, $h_p$-sized features in radius or azimuth are resolved in the outer disk, and detectable down to $r \approx 10$ au (for sufficient contrast).

The second criterion was the ability to detect a ~10% contrast out to solar system size scales ($r \approx 40$ au). This is roughly the contrast measured for the weaker substructures in the HL Tau and TW Hya disks (e.g., Akiyama et al. 2016; Huang et al. 2018c). It is also sufficient to detect the continuum emission that (indirectly) traces the ~20% pressure variations produced by $\gtrsim 0.1 M_{\text{Jup}}$ planets (Fung et al. 2014), zonal flows (Simon & Armitage 2014), or weak vortices (e.g., Goodman et al. 1987), even if (contrary to expectations) there is no accompanying amplification in the concentration of the solids (presuming the emission is optically thin).

The combination of these criteria and ALMA technical restrictions meant that the optimal observing frequency was in the vicinity of 240 GHz (Band 6). Higher frequency observations at comparable (or better) resolution were not permitted for Cycle 4 Large Programs, and the resolution and sensitivity options at lower frequencies were both insufficient for our goals.

The resolution criterion drove planning for the survey sample. The Cycle 4 configuration schedule was set to provide the requisite resolution (with baseline lengths out to 6.8–12.6 km) during 2017 June and July. We targeted disks that are nearby enough to give the required spatial resolution for those configurations, and that transit at high elevations at night during this period. This limited the sample pool to the Oph (Wilking et al. 2008), Lup (Comerón 2008), and Upper Sco (Preibisch & Mamajek 2008) regions, plus a few isolated targets. The field was narrowed to focus on Class II sources to avoid confusion with envelope emission. We excluded “transition” disks, because they are already known to exhibit substructures (by definition).

Those criteria leave ~200 viable targets. A more severe cut was then made to meet the contrast criterion. A general framing of that criterion is somewhat arbitrary, but we chose some fiducial numbers as a guide. For a target at 140 pc and with a synthesized beam FWHM of 35 mas, we aimed to measure a 10% deviation from an otherwise smooth brightness profile (at SNR $\ge 2$ per beam) out at $r = 40$ au (~0.3). This metric requires previous continuum observations at modest (0.3) resolution for selection (Andrews et al. 2009, 2010; Ansdell et al. 2016; Barenfeld et al. 2016). For reasonable assumptions about the shape of the brightness profile,<fnref n="14" /> this criterion can be met with a cut on the 0.3 peak brightness. Experimentation with simulated data suggested a peak brightness cut at 20 mJy per 0.3 beam (4.8 K) is appropriate, implying an objective noise level of 17 µJy per 35 mas beam (0.3 K).<fnref n="15" /> The caveat

is that much of the available data at 0.3 resolution were taken at 340 GHz; in the applicable cases, we assumed that $I_\nu \propto \nu^{2.5}$ (see Andrews & Williams 2005).

While that brightness cut substantially reduces the pool, the sample size was ultimately set by ALMA restrictions. Only ~30 hr in the LST ranges of interest were set aside for Large Programs in each of the two relevant array configurations. The desired noise could be reached in ~1 hr of integration per target, but the factor of three overhead costs meant that the sample size was limited to 10 targets per configuration. We selected 10 targets (mostly) in Oph for the more compact of the two configurations (C40-8, $\le 6.8$ km baselines; 50 mas resolution), based on their nominally closer distances (125 pc; de Geus et al. 1989; Loinard et al. 2008).<fnref n="16" /> Ten more targets (primarily in Lup) were chosen for C40-9 ($\le 12.6$ km baselines; 35 mas resolution).

The resulting sample and its stellar host properties are compiled in Table 1. Target distances ($d$) were derived from Gaia DR2 parallax measurements (Gaia Collaboration et al. 2018), following Astraatmadja & Bailer-Jones (2016) for a flat $d$ prior. Literature estimates of the effective temperatures ($T_{\text{eff}}$) and luminosities ($L_*$; rescaled for the appropriate $d$) were adopted to derive masses ($M_*$) and ages ($t_*$) based on the MIST models (Choi et al. 2016), following the methodology described by Andrews et al. (2018). Accretion rates ($\dot{M}_*$) were calculated from those host parameters and literature measurements of accretion luminosities (scaled for $d$; see Table 1). The sample hosts exhibit a range of young star properties, with $M_* \approx 0.2–2 M_\odot$ and nearly two decades spanned in both $L_*$ and $\dot{M}_*$. The mean age is 1 Myr, although with considerable individual uncertainties (and various untreated systematics; see Soderblom et al. 2014).

The broadband spectral energy distributions (SEDs) for the sample are shown together in Figure 1. Relative to the median SED (normalized at 1.5 µm) of larger samples of Class II targets (e.g., Ribas et al. 2017), RU Lup, AS 205, and AS 209 are in the top quartile (i.e., are over-luminous); the SEDs for HD 143006, SR 4, and DoAr 25 are relatively low in the near-infrared and high in the far-infrared (similar to, though not nearly as pronounced as, the typical transition disk SED); and the SEDs for DoAr 33 and WaOph 6 are in the bottom quartile. This diversity in the SEDs is one potential basis for future explorations of how the resolved emission distributions vary with relevant “bulk” parameters (e.g., the amount of dust settling toward the disk midplane).

While these sample targets do cover a range of properties, it is worth emphasizing that this range is not representative of the general population. The sample hosts tend to have earlier spectral types, and are accordingly more massive, luminous, and accreting more vigorously than stars at the peak of the initial mass function. This host bias enters implicitly with the sensitivity criterion, because we required a previous resolved measurement. The studies that provided those data are biased toward brighter continuum sources, which permeates to the host properties because the continuum luminosity scales steeply with $M_*$ (Andrews et al. 2013; Mohanty et al. 2013) and $\dot{M}_*$ (Manara et al. 2016; Mulders et al. 2017). The same is true for multiple star systems: these were not explicitly

<footnote n="14">¹⁴ We conservatively assumed a face-on orientation with $I_\nu \propto r^{-0.5}-r^{-1}$ (see Tripathi et al. 2017; Andrews et al. 2018).</footnote>

<footnote n="15">¹⁵ Unless specified otherwise, DSHARP brightness temperatures are calculated assuming the standard Rayleigh–Jeans relation.</footnote>

<footnote n="16">¹⁶ Note that the 125 pc distance used to motivate the slightly coarser resolution for Oph targets was inappropriate. However, this aspect of the survey design was ignored anyway, due to unforeseen alterations in the configuration schedule.</footnote>
<table><caption>**Table 1**
DSHARP Sample: Host Star Properties</caption><thead><tr><th>Name<br>(1)</th><th>Region<br>(2)</th><th>2MASS<br>Designation<br>(3)</th><th>d<br>(pc)<br>(4)</th><th>SpT<br>(5)</th><th>log T<sub>eff</sub><br>(K)<br>(6)</th><th>log L<sub>*</sub><br>(L<sub>☉</sub>)<br>(7)</th><th>log M<sub>*</sub><br>(M<sub>☉</sub>)<br>(8)</th><th>log t<sub>*</sub><br>(yr)<br>(9)</th><th>log M&#x0307;<sub>*</sub><br>M<sub>☉</sub> yr<sup>-1</sup><br>(10)</th><th>References<br>(11)</th></tr></thead><tbody><tr><td>HT Lup<sup>a</sup></td><td>Lup I</td><td>J15451286-3417305</td><td>154 ± 2</td><td>K2</td><td>3.69 ± 0.02</td><td>0.74 ± 0.20</td><td>0.23<sup>+0.06</sup><sub>-0.13</sub></td><td>5.9 ± 0.3</td><td>&lt;-8.4</td><td>1, 1, 1</td></tr><tr><td>GW Lup</td><td>Lup I</td><td>J15464473-3430354</td><td>155 ± 3</td><td>M1.5</td><td>3.56 ± 0.02</td><td>-0.48 ± 0.20</td><td>-0.34<sup>+0.10</sup><sub>-0.17</sub></td><td>6.3 ± 0.4</td><td>-9.0 ± 0.4</td><td>1, 1, 1</td></tr><tr><td>IM Lup</td><td>Lup II</td><td>J15560921-3756057</td><td>158 ± 3</td><td>K5</td><td>3.63 ± 0.03</td><td>0.41 ± 0.20</td><td>-0.05<sup>+0.09</sup><sub>-0.13</sub></td><td>5.7 ± 0.4</td><td>-7.9 ± 0.4</td><td>1, 1, 1</td></tr><tr><td>RU Lup</td><td>Lup II</td><td>J15564230-3749154</td><td>159 ± 3</td><td>K7</td><td>3.61 ± 0.02</td><td>0.16 ± 0.20</td><td>-0.20<sup>+0.12</sup><sub>-0.11</sub></td><td>5.7 ± 0.4</td><td>-7.1 ± 0.3</td><td>1, 1, 1</td></tr><tr><td>Sz 114</td><td>Lup III</td><td>J16090185-3905124</td><td>162 ± 3</td><td>M5</td><td>3.50 ± 0.01</td><td>-0.69 ± 0.20</td><td>-0.76<sup>+0.08</sup><sub>-0.07</sub></td><td>6.0<sup>+0.1</sup><sub>-0.8</sub></td><td>-9.1 ± 0.3</td><td>1, 1, 1</td></tr><tr><td>Sz 129</td><td>Lup IV</td><td>J15591647-4157102</td><td>161 ± 3</td><td>K7</td><td>3.61 ± 0.02</td><td>-0.36 ± 0.20</td><td>-0.08<sup>+0.03</sup><sub>-0.15</sub></td><td>6.6 ± 0.4</td><td>-8.3 ± 0.3</td><td>1, 1, 1</td></tr><tr><td>MY Lup<sup>b</sup></td><td>Lup IV</td><td>J16004452-4155310</td><td>156 ± 3</td><td>K0</td><td>3.71 ± 0.02</td><td>-0.06 ± 0.20</td><td>0.09<sup>+0.03</sup><sub>-0.13</sub></td><td>7.0<sup>+0.6</sup><sub>-0.3</sub></td><td>&lt;-9.6</td><td>1, 1, 1</td></tr><tr><td>HD 142666</td><td>Upper Sco</td><td>J15564002-2201400</td><td>148 ± 2</td><td>A8</td><td>3.88 ± 0.02</td><td>0.96 ± 0.21</td><td>0.20<sup>+0.04</sup><sub>-0.01</sub></td><td>7.1 ± 0.3</td><td>&lt;-8.4</td><td>2, 2, 2</td></tr><tr><td>HD 143006</td><td>Upper Sco</td><td>J15583692-2257153</td><td>165 ± 5</td><td>G7</td><td>3.75 ± 0.02</td><td>0.58 ± 0.15</td><td>0.25<sup>+0.05</sup><sub>-0.08</sub></td><td>6.6 ± 0.3</td><td>-8.1 ± 0.4</td><td>3, 4, 5</td></tr><tr><td>AS 205<sup>a</sup></td><td>Upper Sco</td><td>J16113134-1838259</td><td>128 ± 2</td><td>K5</td><td>3.63 ± 0.03</td><td>0.33 ± 0.15</td><td>-0.06<sup>+0.07</sup><sub>-0.05</sub></td><td>5.8 ± 0.3</td><td>-7.4 ± 0.4</td><td>3, 4, 6</td></tr><tr><td>SR 4</td><td>Oph L1688</td><td>J16255615-2420481</td><td>134 ± 2</td><td>K7</td><td>3.61 ± 0.02</td><td>0.07 ± 0.20</td><td>-0.17<sup>+0.11</sup><sub>-0.14</sub></td><td>5.9 ± 0.4</td><td>-6.9 ± 0.5</td><td>7, 8, 9</td></tr><tr><td>Elias 20</td><td>Oph L1688</td><td>J16261886-2428196</td><td>138 ± 5</td><td>M0</td><td>3.59 ± 0.03</td><td>0.35 ± 0.20</td><td>-0.32<sup>+0.12</sup><sub>-0.07</sub></td><td>&lt;5.9</td><td>-6.9 ± 0.5</td><td>9, 10, 9</td></tr><tr><td>DoAr 25</td><td>Oph L1688</td><td>J16262367-2443138</td><td>138 ± 3</td><td>K5</td><td>3.63 ± 0.03</td><td>-0.02 ± 0.20</td><td>-0.02<sup>+0.04</sup><sub>-0.19</sub></td><td>6.3 ± 0.4</td><td>-8.3 ± 0.5</td><td>11, 10, 12</td></tr><tr><td>Elias 24</td><td>Oph L1688</td><td>J16262407-2416134</td><td>136 ± 3</td><td>K5</td><td>3.63 ± 0.03</td><td>0.78 ± 0.20</td><td>-0.11<sup>+0.16</sup><sub>-0.08</sub></td><td>5.3 ± 0.4</td><td>-6.4 ± 0.5</td><td>11, 8, 9</td></tr><tr><td>Elias 27</td><td>Oph L1688</td><td>J16264502-2423077</td><td>116<sup>+19</sup><sub>-10</sub></td><td>M0</td><td>3.59 ± 0.03</td><td>-0.04 ± 0.23</td><td>-0.31<sup>+0.15</sup><sub>-0.11</sub></td><td>5.9 ± 0.5</td><td>-7.2 ± 0.5</td><td>7, 10, 9</td></tr><tr><td>DoAr 33</td><td>Oph L1688</td><td>J16273901-2358187</td><td>139 ± 2</td><td>K4</td><td>3.65 ± 0.03</td><td>0.18 ± 0.20</td><td>0.04<sup>+0.05</sup><sub>-0.17</sub></td><td>6.2 ± 0.4</td><td>...</td><td>13, 8</td></tr><tr><td>WSB 52</td><td>Oph L1688</td><td>J16273942-2439155</td><td>136 ± 3</td><td>M1</td><td>3.57 ± 0.03</td><td>-0.15 ± 0.20</td><td>-0.32<sup>+0.13</sup><sub>-0.17</sub></td><td>5.8 ± 0.5</td><td>-7.6 ± 0.5</td><td>7, 8, 9</td></tr><tr><td>WaOph 6</td><td>Oph N 3a</td><td>J16484562-1416359</td><td>123 ± 2</td><td>K6</td><td>3.62 ± 0.03</td><td>0.46 ± 0.20</td><td>-0.17<sup>+0.17</sup><sub>-0.09</sub></td><td>5.5 ± 0.5</td><td>-6.6 ± 0.5</td><td>14, 10, 14</td></tr><tr><td>AS 209</td><td>Oph N 3a</td><td>J16491530-1422087</td><td>121 ± 2</td><td>K5</td><td>3.63 ± 0.03</td><td>0.15 ± 0.20</td><td>-0.08<sup>+0.11</sup><sub>-0.14</sub></td><td>6.0 ± 0.4</td><td>-7.3 ± 0.5</td><td>15, 10, 6</td></tr><tr><td>HD 163296</td><td>isolated?</td><td>J17562128-2157218</td><td>101 ± 2</td><td>A1</td><td>3.97 ± 0.03</td><td>1.23 ± 0.30</td><td>0.31<sup>+0.05</sup><sub>-0.03</sub></td><td>7.1 ± 0.6</td><td>-7.4 ± 0.3</td><td>2, 2, 2</td></tr></tbody></table>

**Notes.** Column 1: target name. Column 2: associated star-forming region. The Lup subcloud regions are as designated by Cambrésy (1999). Upper Sco memberships were made following Luhman et al. (2018). AS 209 and WaOph 6 are located well northeast of the main Oph region in the Oph N 3a complex. They are most closely associated with the L163 and L162 dark clouds, respectively. Column 3: the 2MASS designations, to aid in catalog cross-referencing. Column 4: distance (computed from the Gaia DR2 parallaxes). Column 5: spectral type from the literature (first reference entries in column 11). Column 6: Effective temperatures from the literature (second reference entries in column 11). Column 7: stellar luminosities from the literature, scaled according to the appropriate *d* in column 4 (second reference entries in column 11). Columns 8 and 9: stellar masses and ages. Column 9: accretion rates, inferred from (properly scaled) accretion luminosities (third reference entries in column 11). All quoted measurements correspond to the peak of the marginalized posterior distributions. Uncertainties reflect the 68.3% confidence interval; limits are taken at the 95.5% confidence level.

**References.** In column 11, the references for the quoted SpT, {T<sub>eff</sub>, L<sub>*</sub>}, and accretion luminosity measurements, respectively: (1) Alcalá et al. (2017); (2) Fairlamb et al. (2015); (3) Luhman & Mamajek (2012); (4) Barenfeld et al. (2016); (5) Rigliaco et al. (2015); (6) Salyk et al. (2013); (7) Luhman & Rieke (1999); (8) Andrews et al. (2010); (9) Natta et al. (2006); (10) Andrews et al. (2009); (11) Wilking et al. (2005); (12) Muzerolle et al. (1998); (13) Bouvier & Appenzeller (1992); (14) Eisner et al. (2005); (15) Herbig & Bell (1988).

<footnote><sup>a</sup> HT Lup (Sz 68) and AS 205 (V866 Sco) are triple systems. See Kurtovic et al. (2018) for details.</footnote>

<footnote><sup>b</sup> The MY Lup disk is inclined and flared enough that it likely extincts the host: the *L*<sub>*</sub> and *t*<sub>*</sub> estimates may be too faint and old, respectively.</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Broadband SEDs for the DSHARP targets. The ordinate is $L_{\nu} = 4\pi d^2\nu F_{\nu}$ in $L_{\odot}$ units. These SEDs have been dereddened using the extinction values quoted by the references in column 11 of Table 1 (second entries) and the prescription described by Andrews et al. (2013). Blue curves show the NEXTGEN/ BT-sett1 photosphere models (Allard et al. 2003, 2011) corresponding to the stellar parameters listed in Table 1. Red curves show the *Spitzer* IRS spectra. Note that the SEDs for HT Lup and AS 205 include contributions from multiple components. Optical photometry was collected from a range of sources (Vrba et al. 1993; Herbst et al. 1994; Hughes et al. 1994; Oudmaijer et al. 2001; Eisner et al. 2005; Gras-Velázquez & Ray 2005; Wilking et al. 2005; Padgett et al. 2006; Grankin et al. 2007; Merín et al. 2008; Mendigutía et al. 2012); infrared data were culled from 2MASS (Skrutskie et al. 2006), WISE (Wright et al. 2010), Spitzer imaging surveys (Carpenter et al. 2008; Evans et al. 2009), AKARI (Ishihara et al. 2010), and Herschel (IRSA); (sub)millimeter data come from various sources (Andre & Montmerle 1994; Mannings & Emerson 1994; Sylvester et al. 1996; Mannings & Sargent 1997; Nuernberger et al. 1997; Dent et al. 1998; Henning et al. 1998; Natta et al. 2004; Stanke et al. 2006; Andrews & Williams 2007; Isella et al. 2007, 2009; Lommen et al. 2007, 2009; Pinte et al. 2008; Andrews et al. 2009, 2018; Roccatagliata et al. 2009; Ricci et al. 2010; Öberg et al. 2011; Sandell et al. 2011; Pérez et al. 2012, 2015; Qi et al. 2015; Ansdell et al. 2016; Barenfeld et al. 2016; Cleeves et al. 2016; Cox et al. 2017; Huang et al. 2017; Tripathi et al. 2017; Ubach et al. 2017). These SEDs are available in the DSHARP data release.</figcaption></figure>

excluded, but the selection criteria bias against them because close companions tend to reduce the system continuum emission (e.g., Jensen et al. 1994; Harris et al. 2012).

The bias in favor of targets with brighter continuum emission also translates into a preferential selection of larger disks, given the observed size–luminosity correlation (Tazzari et al. 2017; Tripathi et al. 2017; Andrews et al. 2018). This corresponding size bias is decidedly beneficial for achieving the DSHARP goals discussed in Section 1. As we noted above, the general theoretical predictions for substructure sizes are comparable to the gas pressure scale height ($h_P$), which increases roughly linearly with disk radius. For a fixed resolution, it should be easier to identify and characterize the larger substructures expected at larger disk radii.

To roughly quantify these biases, we can make a comparison between targets that are more representative of the general disk population and the average member of the DSHARP sample. A “typical” target has a host star mass near the peak of the mass function ($M_* \approx 0.3 M_\odot$, or spectral type M3–M4) and continuum emission from its disk that is both relatively faint ($F_\nu \approx 10–15$ mJy; Ansdell et al. 2016; Cieza et al. 2019) and compact ($R_{eff} \approx 10–20$ au, with the effective radius defined by Tripathi et al. 2017; see also Andrews et al. 2018). The DSHARP averages are $M_* \approx 0.8 M_\odot$ (spectral type K7),

$F_\nu \approx 150$ mJy, and $R_{eff} \approx 50$ au; only the sample extremes stretch down toward “typical” values.

These biases are difficult to mitigate for studies focused on finding and characterizing disk substructures, presuming their size scales are usually comparable to $h_P$. The “typical” disk is compact enough that $h_P$ for the radii where there is still continuum emission is smaller than the best resolutions available with ALMA. If this is the case, then we could be left probing only the extreme large end of substructures in “typical” disks, making any assessments of prevalence difficult (i.e., failed searches for substructures would still permit plenty of $h_P$-size features to be present on subresolution scales). One option is to push to higher frequencies and thereby better resolution, but then high optical depths would limit the discovery-space to substructures in the form of dramatic depletions (e.g., very deep gaps) only.

### 3. Observations

The DSHARP ALMA observations were conducted in 2017 May–November as part of program 2016.1.00484.L. All measurements used the Band 6 receivers and correlated data from four spectral windows (SPWs) in the dual polarization mode. The continuum was sampled in three SPWs, centered at 232.6, 245.0, and 246.9 GHz, each with 128 channels spanning
<table><caption>**Table 2**
DSHARP Observing Log (ALMA Program 2016.1.00484.L)</caption><thead><tr><th>Name<br>(1)</th><th>UTC Date<br>(2)</th><th>Config.<br>(3)</th><th>Baselines<br>(4)</th><th>N<sub>ant</sub><br>(5)</th><th>E/°<br>(6)</th><th>PWV/mm<br>(7)</th><th>Calibrators<br>(8)</th></tr></thead><tbody><tr><td rowspan="4">HT Lup</td><td>2017 May 14–04:11</td><td>C40-5</td><td>15 m–1.1 km</td><td>43</td><td>76–77</td><td>1.00–1.15</td><td>J1517–2422, J1427–4206, J1610–3958, J1540–3906</td></tr><tr><td>2017 May 17–02:12</td><td>C40-5</td><td>15 m–1.1 km</td><td>49</td><td>58–67</td><td>0.90–1.05</td><td>J1517–2422, J1517–2422, J1610–3958, J1540–3906</td></tr><tr><td>2017 Sep 24–17:39</td><td>C40-8/9</td><td>41 m–12.1 km</td><td>39</td><td>59–70</td><td>0.60–1.15</td><td>J1517–2422, J1517–2422, J1534–3526, J1536–3151</td></tr><tr><td>2017 Sep 24–19:12</td><td>C40-8/9</td><td>41 m–12.1 km</td><td>39</td><td>75–78</td><td>0.65–1.05</td><td>J1517–2422, J1427–4206, J1534–3526, J1536–3151</td></tr><tr><td>GW Lup</td><td>2017 May 14–04:11</td><td>C40-5</td><td>15 m–1.1 km</td><td>43</td><td>69–72</td><td>1.00–1.15</td><td>J1517–2422, J1427–4206, J1610–3958, J1540–3906</td></tr></tbody></table>

**Note.** Column 1: target name. Column 2: UTC date and time at the start of the observations. Column 3: ALMA configuration. Column 4: minimum and maximum baseline lengths. Column 5: number of antennas available. Column 6: target elevation range. Column 7: range of precipitable water vapor levels. Column 8: from left to right, the quasars observed for calibrating the bandpass, amplitude scale, phase variations, and checking the phase transfer. Additional archival observations used in our analysis are compiled in Table 3.

(This table is available in its entirety in machine-readable form.)

<table><caption>**Table 3**
Archival ALMA Data Sets Used by DSHARP</caption><thead><tr><th>Name<br>(1)</th><th>UTC Date<br>(2)</th><th>Config.<br>(3)</th><th>Baselines<br>(4)</th><th>N<sub>ant</sub><br>(5)</th><th>Calibrators<br>(6)</th><th>Program<br>(7)</th><th>References<br>(8)</th></tr></thead><tbody><tr><td rowspan="5">IM Lup</td><td>2014 Jul 06–22:18</td><td>C34-4</td><td>20–650 m</td><td>31</td><td>J1427–4206, Titan, J1534–3526, J1626–2951</td><td>2013.1.00226.S</td><td>1</td></tr><tr><td>2014 Jul 17–01:38</td><td>C34-4</td><td>20–650 m</td><td>32</td><td>J1427–4206, Titan, J1534–3526, ...</td><td>2013.1.00226.S</td><td>1</td></tr><tr><td>2015 Jan 29–09:48</td><td>C34-2/1</td><td>15–349 m</td><td>40</td><td>J1517–2422, Titan, J1610–3958, ...</td><td>2013.1.00694.S</td><td>2</td></tr><tr><td>2015 May 13–08:30</td><td>C34-3/4</td><td>21–558 m</td><td>36</td><td>J1517–2422, Titan, J1610–3958, ...</td><td>2013.1.00694.S</td><td>2</td></tr><tr><td>2015 Jun 09–23:42</td><td>C34-5</td><td>21–784 m</td><td>37</td><td>J1517–2422, Titan, J1610–3958, J1614–3543</td><td>2013.1.00798.S</td><td>3</td></tr></tbody></table>

**Note.** Column 1: target name. Column 2: UTC date and time at the start of the observations. Column 3: ALMA configuration. Column 4: range of baseline lengths. Column 5: number of antennas available. Column 6: from left to right, the quasars observed for calibrating the bandpass, amplitude scale, phase variations, and checking the phase transfer. An entry of “...” indicates no calibrator was observed for checking the phase transfer. Column 7: ALMA program ID. Column 8: original references for these data sets.

**References.** (1) Öberg et al. (2015); (2) Cleeves et al. (2017); (3) Pinte et al. (2018); (4) Salyk et al. (2014); (5) Dipierro et al. (2018); (6) Pérez et al. (2016); (7) Huang et al. (2016); (8) Fedele et al. (2018); (9) Flaherty et al. (2015); (10) Isella et al. (2016).

(This table is available in its entirety in machine-readable form.)

1.875 GHz (31.25 MHz per channel). The remaining SPW was centered at the $^{12}$CO $J=2-1$ rest frequency (230.538 GHz) and covered a bandwidth of 938 MHz in 3840 channels (488 kHz channel spacing, 0.64 km s$^{-1}$ velocity resolution). The plan was to observe each target briefly in the C40-5 (hereafter “compact”) configuration, and also for ~1 hr in the C40-8 or C40-9 (hereafter “extended”) configurations. The compact observations are necessary to recover emission on the larger angular scales that are not sampled in the extended configurations. The actual observing log is provided in Table 2.

The compact observations used an array with baseline lengths from 15 m to 1.1 km (a resolution of ~0.25′′). The FWHM continuum and CO (per channel) emission scales are ≲2′′, so spatial filtering should be negligible (see Wilner & Welch 1994). These observations cycled between nearby targets and totaled ~12 minutes of integration time per source. A nearby phase calibrator was observed every 6 minutes; an additional “check” calibrator (to assess the quality of phase transfer) was observed every 30 minutes. A bandpass and amplitude calibrator (sometimes the same quasar) were observed during each observing block. The log in Table 2 includes information about the observing conditions and calibrators.

We relied on archival ALMA observations of five targets (IM Lup, HD 142666, Elias 24, Elias 27, and HD 163296) instead of obtaining new compact data, and folded in archival data for three other targets (HD 143006, AS 205, and AS 209). Information about these data sets are compiled in Table 3.

The setups, observing strategies, and weather conditions are described in the listed references.

Due to a long stretch of inclement weather, the extended configuration observations were delayed until 2017 September, and continued through November. Despite the nonoptimal scheduling for the DSHARP sample, nearly all of the targets were observed for two executions (often in different configurations) in good conditions, each with ~35 minutes of on-source integration time. Sz 114, AS 205, and DoAr 25 each had only a single successful execution. The spectral setup was the same as that for the compact data sets. Observations cycled between a single target and a nearby phase calibrator on 1 minute intervals, with a “check” calibrator visited every 30 minutes. Bandpass and amplitude calibrators were observed in each execution block (see Table 2).

## 4. Calibration and Imaging

Since the DSHARP survey is among the first to collect a large volume of ALMA data on such long baselines, a substantial effort was made to explore various calibration strategies to enhance the data quality. The standard methodology we adopted is described here. The specific details on the calibration of data sets for individual targets (i.e., calibration scripts) are available in the DSHARP data release (see Section 5). All calibration tasks are performed with the CASA package (McMullin et al. 2007) and a small supplement of `python` tasks.
## 4.1. Pipeline Calibration

The first step was a standard ALMA pipeline calibration. This procedure was performed by ALMA staff separately for the compact and extended data, using CASA v4.7.2 or v5.1.1 for data sets that were processed before or after 2017 November, respectively. The pipeline imports the raw data and flags problematic scans, channels, or antennas. It then derives a table of system temperatures ($T_{\text{sys}}$). Most of the DSHARP data have $T_{\text{sys}} \approx 60–80$ K; in the poorest conditions it reached 130 K, and in the best cases it was 50 K. Next, the pipeline adjusts the visibility phases according to water vapor radiometer (WVR) measurements. For the extended data, the WVR corrections improved the median rms phase variations by a factor of ~1.7, although individual data sets saw improvements between 1.2 and 3. The corrected rms phase variations (far from the reference antenna) were typically 30° (with a range ~15°–50°). The compact observations saw similar improvement factors (1.5–3.0) and rms phase variations (~10°).

The pipeline then performs a bandpass calibration, using the first quasar in the calibrator list in Table 2. It continues by setting the amplitude scale, using measurements of the second quasar in the Table 2 list. The flux density in each SPW for that quasar is determined from a power-law spectral model based on bi-monthly monitoring in ALMA Bands 3 and 7 (~100 and 340 GHz) that is tied to primary calibrators (planets or moons). Finally, the gain variations with time are corrected, with reference to repeated measurements of the nearby quasar listed third in the Table 2 calibrator list.

## 4.2. Self-calibration

We next performed some substantial post-processing, with particular emphasis on combining data sets (from different array configurations and observations) and self-calibrating the visibilities. We generally followed the homogenized strategy described below, using CASA v5.1.1 and a set of custom `python` routines.

The procedure started with the compact data. A pseudo-continuum data set was created by flagging data within $\pm 25 \text{ km s}^{-1}$ from the CO $J = 2-1$ line center and averaging into 125 MHz channels. The visibilities corresponding to each individual observation were imaged (Section 4.3) and checked to ensure consistent astrometric registration and flux calibration (if necessary, they are corrected; see Section 4.4). The individual data sets were then recombined. Next, we performed a series of phase-only self-calibration iterations, stepping down the solution interval (60, 30, 18, and 6 s). Reference antennas were selected based on data quality and proximity to the array center. When possible, we avoided combining SPWs (or scans) to correct for SPW-dependent gain variations. After each iteration, the data were imaged. A noise estimate was made in an annular region within a 4"-25-radius circle centered on the target but excluding the image mask. This self-calibration sequence is stopped after reaching a solution interval on the record length (6 s) or if the peak SNR does not increase by >5% from the previous iteration. Finally, we performed one iteration of amplitude self-calibration (for each SPW independently) on a scan interval (~6 minutes). The (phase + amplitude) self-calibration provided a dramatic improvement in quality. The typical peak SNR increased by a factor of 3; the

resulting noise was $30 \ \mu\text{Jy beam}^{-1}$ (10 mK) for a ~0"-25 beam. The same procedure was applied to archival data sets.

Next, we prepared the extended data as was described above, with an additional time-averaging to 6 s integrations (from the original 2 s records). The data for each individual extended observation were imaged and checked for misalignments and flux discrepancies. Once those are corrected (if necessary; see Section 4.4), the compact (already self-calibrated) and extended data sets were combined. The phases for this combined data set were iteratively self-calibrated on solution intervals of {900, 360, 180, 60, 30 s} (usually only the latter three are necessary). The SPWs were combined in this case to enhance the SNR on longer baselines. For antenna pairings with SNR $\le 1.5$ on these intervals, the self-calibration solutions were not applied but the corresponding data were not flagged (`applymode = 'calonly'` in the `applycal` task). The sequence was stopped when the peak SNR does not increase by >5% and the map quality does not visually improve. One iteration of amplitude self-calibration was attempted on the starting interval of the phase self-calibration sequence.

This self-calibration of the combined data sets resulted in a typical improvement of 40% in the peak SNR, although there is a large range in benefits across the sample. The improvements are generally smaller here because the compact data were already self-calibrated and the extended data were taken in excellent conditions. The typical noise measured in the combined, self-calibrated data sets is $10–20 \ \mu\text{Jy beam}^{-1}$ (0.1–0.5 K).

Once the continuum self-calibration was satisfactory, the same gain tables are applied to the nonspectrally averaged visibilities (after any required astrometric and flux calibration adjustments) to obtain a corresponding calibrated measurement set for the region of the spectrum around the CO $J = 2-1$ emission line.

## 4.3. Imaging during Self-calibration

Self-calibration uses continuum emission models assembled from the “clean” components derived from interferometric imaging. We adopted a set of imaging standards to homogenize that process. These were informed by considerable experimentation with the associated parameter choices. We explored alternative sets of deconvolution scales, clean thresholds, masks, and pixel sizes and found that reasonable other options had negligible influence on the end products of self-calibration.

All imaging was performed with the `tclean` task. For the compact data, we imaged out to the primary beam FWHM ($26''$) with 30 mas pixels (~10 per synthesized beam FWHM, $\theta_b$) to check for problematic background sources. Finding nothing of concern, we used 9"-wide images with 3 mas pixels (again, ~10 pixels per $\theta_b$) for the combined data sets. We used the multiscale, multifrequency synthesis (assuming a flat spectrum) deconvolution mode (Cornwell 2008) with a Briggs `robust = 0.5` weighting scheme. Elliptical masks were designed to reflect the target geometry (aspect ratio and position angle) and pad the outer reaches of the emission distribution. The adopted (Gaussian) deconvolution scales are target-dependent, but always include a point-like contribution and scales comparable to $\theta_b$ and $2-3 \times \theta_b$; additional scales (increasing by factors of 2–3) could be selected up to the mask radius. The algorithm was halted on thresholds; $3 \times$ the noise early in the self-calibration sequence, and $2 \times$ the noise for the last phase-only step and the amplitude self-calibration.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Illustration of the effects of spatial frequency coverage and visibility weighting on PSF structure and image morphology, for the RU Lup disk. The top panels show the observed *u*, *v* coverage. The middle panels show the corresponding PSF structures, annotated with the corresponding *robust* weighting parameter. The bottom panels show the corresponding images, on the same $T_b$ scale, created from subsets of the end product of the self-calibration. FWHM beam dimensions are marked in the lower left corners of each image.</figcaption></figure>

Special effort was made to verify that sidelobes in the point-
spread function (PSF, or “dirty” beam) do not corrupt the self-
calibration. The extended ALMA configurations place antennas
along three distinct arms (set by the site topography). The
corresponding spatial frequency coverage generates compli-
cated PSF features, with sidelobes up to ~30%. Figure 2
illustrates the impact, showing the connections between the
sampling function (*u*, *v* coverage), PSF, and image for different
configurations and weighting schemes. We vetted the effects of
those PSF features on self-calibration by repeating the process
for different combinations of weighting schemes and tapers.
Coupling lower *robust* values with tapers can mitigate PSF
artifacts while maintaining resolution, but at a substantial SNR
cost. Direct comparisons (of both visibilities and images)
between these variants and the standard methodology outlined
above demonstrated that the PSF features had negligible impact
on the self-calibration.<fnref n="17" />

While the effects on self-calibration are minimal, the
resulting images can still exhibit PSF-related artifacts. One of
the more interesting is the imprint of a hexagonal structure on

emission rings (e.g., second image from left, bottom row of
Figure 2), produced by convolution with a “spoked” PSF (a
consequence of the extended ALMA configuration arms). As
demonstrated in the bottom right panel, this can usually be
minimized with an appropriate visibility weighting and/or
tapering.

4.4. Astrometric and Flux Scale Alignment

Half the sample targets show clear spatial offsets between
their emission centers in different observations. For the larger
of these shifts (~100 mas), the cause is proper motion
(especially when using archival data); in other cases, smaller
(10–30 mas) mismatches might instead be attributed to
instrumental or atmospheric artifacts. Combining these data
sets without correcting these shifts creates blurred (or even
double) images, which is problematic when they are used as
initial self-calibration models. The solution is to simply adjust
the visibility phases to shift into alignment. We measure
emission centroid positions with Gaussian fits in the image
plane for each individual observation and calculate the offsets
relative to the highest quality extended data set. The *fixvis*
task then implements the appropriate phase adjustments. In
cases where the observations have different pointing centers,
we manually reconcile them with the *fixplanets* task.

<footnote n="17">¹⁷ The HD 163296 disk is the one exception (albeit a quite modest one): we find ~10% SNR improvements (relative to the standard) when self-calibration is conducted for images with robust = -0.5, due to the combination of the target emission distribution and the unusual spatial frequency coverage from the archival data.</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Gallery of 240 GHz (1.25 mm) continuum emission images for the disks in the DSHARP sample. Beam sizes and 10 au scalebars are shown in the lower left and right corners of each panel, respectively. All images are shown with an asinh stretch to reduce the dynamic range (accentuate fainter details without oversaturating the bright emission peaks). For more quantitative details regarding the image dimensions and intensity scales, see Huang et al. (2018a) and Kurtovic et al. (2018).</figcaption></figure>
<table><caption>**Table 4**
DSHARP Fiducial Continuum Image Properties</caption><thead><tr><th>Name<br>(1)</th><th>ν<br>(GHz)<br>(2)</th><th>θ<sub>b</sub>, PA<sub>b</sub><br>(mas, °)<br>(3)</th><th>Rms Noise<br>(μJy beam<sup>-1</sup>, K)<br>(4)</th><th>Peak I<sub>ν</sub>, T<sub>b</sub><br>(mJy beam<sup>-1</sup>, K)<br>(5)</th><th>F<sub>ν</sub><br>(mJy)<br>(6)</th><th>robust<br>(7)</th><th>θ<sub>tap</sub>, PA<sub>tap</sub><br>(mas, °)<br>(8)</th><th>References<br>(9)</th></tr></thead><tbody><tr><td>HT Lup</td><td>239.0</td><td>38 × 33, 61</td><td>14, 0.24</td><td>8.25, 140</td><td>77</td><td>0.5</td><td>...</td><td>IV</td></tr><tr><td>GW Lup</td><td>239.0</td><td>45 × 43, 1</td><td>15, 0.17</td><td>3.35, 37</td><td>89</td><td>0.5</td><td>35 × 15, 0</td><td>II</td></tr><tr><td>IM Lup</td><td>239.0</td><td>44 × 43, 115</td><td>14, 0.16</td><td>7.11, 80</td><td>253</td><td>0.5</td><td>33 × 26, 138</td><td>II, III</td></tr><tr><td>RU Lup</td><td>239.0</td><td>25 × 24, 129</td><td>21, 0.73</td><td>3.45, 123</td><td>203</td><td>-0.5</td><td>22 × 10, 174</td><td>II</td></tr><tr><td>Sz 114</td><td>239.0</td><td>67 × 28, 92</td><td>19, 0.22</td><td>3.36, 38</td><td>49</td><td>0.5</td><td>...</td><td>II</td></tr><tr><td>Sz 129</td><td>239.0</td><td>44 × 31, 94</td><td>15, 0.24</td><td>0.96, 15</td><td>86</td><td>0.0</td><td>...</td><td>II</td></tr><tr><td>MY Lup</td><td>239.0</td><td>44 × 43, 122</td><td>16, 0.18</td><td>1.78, 20</td><td>79</td><td>0.0</td><td>39 × 15, 163</td><td>II</td></tr><tr><td>HD 142666</td><td>231.9</td><td>32 × 22, 62</td><td>13, 0.35</td><td>1.28, 41</td><td>130</td><td>0.5</td><td>...</td><td>II</td></tr><tr><td>HD 143006</td><td>239.0</td><td>46 × 45, 51</td><td>15, 0.15</td><td>0.67, 7</td><td>59</td><td>0.0</td><td>42 × 20, 172</td><td>II, X</td></tr><tr><td>AS 205</td><td>233.7</td><td>38 × 25, 95</td><td>16, 0.38</td><td>6.15, 145</td><td>358</td><td>0.5</td><td>...</td><td>IV</td></tr><tr><td>SR 4</td><td>239.0</td><td>34 × 34, 10</td><td>25, 0.46</td><td>3.40, 63</td><td>69</td><td>-0.5</td><td>35 × 10, 0</td><td>II</td></tr><tr><td>Elias 20</td><td>239.0</td><td>32 × 23, 76</td><td>15, 0.44</td><td>2.59, 75</td><td>104</td><td>0.0</td><td>...</td><td>II</td></tr><tr><td>DoAr 25</td><td>239.0</td><td>41 × 22, 70</td><td>13, 0.31</td><td>1.35, 32</td><td>246</td><td>0.5</td><td>...</td><td>II</td></tr><tr><td>Elias 24</td><td>231.9</td><td>37 × 34, 82</td><td>19, 0.49</td><td>4.63, 119</td><td>352</td><td>0.0</td><td>35 × 10, 166</td><td>II</td></tr><tr><td>Elias 27</td><td>231.9</td><td>49 × 47, 47</td><td>14, 0.14</td><td>4.83, 48</td><td>330</td><td>0.5</td><td>40 × 20, 173</td><td>II, III</td></tr><tr><td>DoAr 33</td><td>239.0</td><td>37 × 24, 75</td><td>17, 0.41</td><td>1.89, 46</td><td>35</td><td>0.0</td><td>20 × 10, 167</td><td>II</td></tr><tr><td>WSB 52</td><td>239.0</td><td>33 × 27, 74</td><td>16, 0.38</td><td>2.60, 62</td><td>67</td><td>0.0</td><td>...</td><td>II, III</td></tr><tr><td>WaOph 6</td><td>239.0</td><td>58 × 54, 84</td><td>17, 0.12</td><td>8.67, 59</td><td>161</td><td>0.0</td><td>55 × 10, 10</td><td>II, III</td></tr><tr><td>AS 209</td><td>239.0</td><td>38 × 36, 68</td><td>19, 0.30</td><td>1.83, 29</td><td>288</td><td>-0.5</td><td>37 × 10, 162</td><td>II, VIII</td></tr><tr><td>HD 163296</td><td>239.0</td><td>48 × 38, 82</td><td>23, 0.27</td><td>4.26, 50</td><td>715</td><td>-0.5</td><td>...</td><td>II, IX</td></tr></tbody></table>

**Note.** Column 1: target name. Column 2: mean frequency. Column 3: synthesized beam FWHM and position angle. Column 4: rms noise in the map, as described in Section 4.3. Column 5: peak intensity in the map. Note that noise and peak brightness temperatures are calculated assuming the Rayleigh–Jeans limit. Column 6: integrated flux density inside the image mask. Column 7: Briggs robust value. Column 8: FWHM and position angle of the taper (if applicable).

**References.** (II) Huang et al. (2018a); (III) Huang et al. (2018b); (IV) Kurtovic et al. (2018); (VIII) Guzmán et al. (2018); (IX) Isella et al. (2018); (X) Pérez et al. (2018).

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Deprojected, azimuthally averaged radial brightness temperature profile for the 240 GHz continuum emission from the AS 209 disk (see Guzmán et al. 2018; Huang et al. 2018a, for more details). The corresponding image is shown in the bottom row of Figure 3, second from the right. The PSF profile (resolution) is marked in black in the upper right corner, along with a gray Gaussian profile that has FWHM = 10 au, to illustrate that the disk substructures typically have compact dimensions.</figcaption></figure>

We also routinely found mismatches in the amplitude scales among different observations of a target. Some experimentation showed that noticeably improved self-calibration results were obtained if the relative flux scales between observations were consistent within 5%. To quantify any mismatches, we inspected the deprojected (according to the Gaussian fit geometries noted above), azimuthally averaged visibilities from different data sets on 200–500 kλ baseline lengths (at lower spatial frequencies, the

extended configuration data are too sparse, and at higher frequencies the averages are more strongly affected by low SNR and phase noise).

These mismatches are caused by inaccurate flux calibration. The claimed calibration accuracy is ~10%, although the adopted methodology for estimating calibrator fluxes (interpolation in time and frequency) can lead to some added uncertainty. About a third of the sample had 5%–10% mismatches, but the majority exhibited 15%–25% discrepancies for at least one data set. In some cases, these were tracked down to a bookkeeping issue: the data were pipeline-processed before a relevant calibrator catalog update. Some 2017 November data sets that used J1427–4206 as the calibrator were problematic. There is no obvious error in the calibrator catalog, so the issue must be with the interpolation: perhaps this quasar flared or changed its spectrum between catalog entries. Regardless of the cause, these misalignments were rectified. We selected a reference data set and used the gaincal task to rescale the outlier data sets.

### 4.5. Fiducial Images

After the calibration was complete, we synthesized a set of fiducial images for further analysis. The continuum imaging followed the methodology outlined in Section 4.3, but was tailored to individual sources with the aim of minimizing PSF artifacts. In many cases, this involved adopting a visibility weighting scheme that traded SNR for resolution, as well as a visibility taper to improve the PSF symmetry. Table 4 lists the basic parameters and resulting properties of these fiducial images. A gallery of the continuum images are shown in Figure 3. Small-scale substructures are notable in all of the
<table><caption>**Table 5**

DSHARP Fiducial CO Datacube Properties</caption><thead><tr><th>Name<br>(1)</th><th>θ<sub>b</sub>, PA<sub>b</sub><br>(mas, °)<br>(2)</th><th>Rms Noise<br>(mJy beam<sup>-1</sup>, K)<br>(3)</th><th>Peak I<sub>t</sub>, T<sub>b</sub><br>(mJy beam<sup>-1</sup>, K)<br>(4)</th><th>robust<br>(5)</th><th>θ<sub>tap</sub>, PA<sub>tap</sub><br>(mas, °)<br>(6)</th><th>Comments<br>(7)</th><th>References<br>(8)</th></tr></thead><tbody><tr><td>HT Lup</td><td>53 × 50, 66</td><td>1.2, 10.1</td><td>12.8, 111</td><td>0.5</td><td>...</td><td>cloud (severe); &lt;150 m filter</td><td>IV</td></tr><tr><td>GW Lup</td><td>109 × 81, 97</td><td>1.4, 3.5</td><td>20.9, 54</td><td>1.0</td><td>40 × 40, ...</td><td>...</td><td>...</td></tr><tr><td>IM Lup</td><td>122 × 115, 47</td><td>1.9, 3.2</td><td>34.4, 56</td><td>0.0</td><td>100 × 100, ...</td><td>...</td><td>III</td></tr><tr><td>RU Lup</td><td>95 × 83, 72</td><td>1.2, 3.4</td><td>49.7, 146</td><td>1.0</td><td>40 × 40, ...</td><td>cloud (mild), complex outflow</td><td>...</td></tr><tr><td>Sz 114</td><td>130 × 90, 105</td><td>2.0, 4.0</td><td>26.8, 53</td><td>1.0</td><td>40 × 40, ...</td><td>cloud (moderate)</td><td>...</td></tr><tr><td>Sz 129</td><td>110 × 83, 76</td><td>1.0, 2.6</td><td>15.3, 38</td><td>1.0</td><td>40 × 40, ...</td><td>...</td><td>...</td></tr><tr><td>MY Lup</td><td>100 × 82, 79</td><td>1.1, 3.1</td><td>15.4, 43</td><td>1.0</td><td>40 × 40, ...</td><td>cloud (mild)</td><td>...</td></tr><tr><td>HD 142666</td><td>77 × 61, 81</td><td>1.3, 6.3</td><td>12.4, 61</td><td>1.0</td><td>40 × 40, ...</td><td>...</td><td>...</td></tr><tr><td>HD 143006</td><td>66 × 49, 84</td><td>1.0, 7.1</td><td>11.2, 80</td><td>0.8</td><td>20 × 20, ...</td><td>...</td><td>X</td></tr><tr><td>AS 205</td><td>115 × 92, 93</td><td>1.4, 3.2</td><td>75.9, 166</td><td>1.0</td><td>40 × 40, ...</td><td>...</td><td>IV</td></tr><tr><td>SR4</td><td>111 × 87, 90</td><td>1.5, 3.5</td><td>29.9, 71</td><td>1.0</td><td>40 × 40, ...</td><td>cloud (moderate)</td><td>...</td></tr><tr><td>Elias 20</td><td>102 × 72, 88</td><td>1.8, 5.6</td><td>27.0, 85</td><td>1.0</td><td>40 × 40, ...</td><td>cloud (severe), outflow</td><td>...</td></tr><tr><td>DoAr 25</td><td>101 × 78, 87</td><td>1.3, 3.9</td><td>18.9, 55</td><td>1.0</td><td>40 × 40, ...</td><td>cloud (moderate)</td><td>...</td></tr><tr><td>Elias 24</td><td>94 × 60, 90</td><td>1.5, 6.2</td><td>23.6, 97</td><td>1.0</td><td>40 × 40, ...</td><td>cloud (severe)</td><td>...</td></tr><tr><td>Elias 27</td><td>132 × 111, 123</td><td>1.6, 2.5</td><td>44.9, 71</td><td>1.0</td><td>100 × 70, 145</td><td>cloud (moderate), envelope?</td><td>III</td></tr><tr><td>DoAr 33</td><td>103 × 79, 88</td><td>1.3, 3.6</td><td>15.8, 45</td><td>1.0</td><td>40 × 40, ...</td><td>cloud (mild)</td><td>...</td></tr><tr><td>WSB 52</td><td>114 × 80, 90</td><td>1.1, 2.8</td><td>31.4, 79</td><td>1.0</td><td>40 × 40, ...</td><td>cloud (severe), complex outflow</td><td>...</td></tr><tr><td>WaOph 6</td><td>126 × 115, 100</td><td>1.3, 2.1</td><td>41.0, 65</td><td>0.5</td><td>100 × 30, 17</td><td>cloud (mild)</td><td>III</td></tr><tr><td>AS 209</td><td>95 × 72, 96</td><td>0.9, 2.9</td><td>21.4, 72</td><td>1.0</td><td>25 × 10, 10</td><td>cloud (mild)</td><td>VIII</td></tr><tr><td>HD 163296</td><td>104 × 95, 100</td><td>0.8, 1.9</td><td>40.7, 95</td><td>0.5</td><td>...</td><td>...</td><td>IX</td></tr></tbody></table>

**Note.** Column 1: target name. Column 2: synthesized beam FWHM and position angle. Column 3: rms noise per channel, measured as described in Section 4.3. HD 143006 and HD 163296 are imaged with 0.32 km s⁻¹ channels; for all other targets, we used 0.35 km s⁻¹ channels. Column 4: peak intensity. Note that noise and peak brightness temperatures are calculated assuming the Rayleigh-Jeans limit. Column 5: Briggs robust value. Column 6: FWHM and position angle of the adopted taper (if applicable). Column 7: comments on issues with the channel maps, including degree of contamination from the ambient molecular cloud and the presence of nondisk features.

**References.** (III) Huang et al. (2018b); (IV) Kurtovic et al. (2018); (VIII) Guzmán et al. (2018); (IX) Isella et al. (2018); (X) Pérez et al. (2018).

DSHARP targets, often with compact (FWHM ≲ 10 au) dimensions. Figure 4 emphasizes the utility of pushing the ALMA resolution for recovering such features in one particularly illustrative example.

We also synthesized channel maps of the CO $J = 2-1$ emission following the basic steps outlined above. The self-calibrated CO visibilities were continuum-subtracted and imaged in LSRK velocity channels at roughly the native channel spacing (0.35 km s⁻¹; the actual velocity resolution is about two channels, due to Hanning smoothing in the ALMA correlator). The DSHARP data are generally not sensitive enough to reconstruct useful channel maps of the emission line at the best available resolution. We compromised by increasing the relative weight of shorter baselines and employing a modest taper. Table 5 lists the imaging parameters, and Figure 5 shows the channel maps. For many of the targets, the CO channel maps exhibit partially recovered large-scale emission structures from the ambient cloud material. These are noted in Table 5 to prevent confusion in the interpretation of extended emission features in some cases (e.g., Elias 24 and WSB 52 are particularly problematic cases).

## 5. Data Release

One key inspiration for conducting the DSHARP survey was to provide a set of resources to the community that can seed and develop a range of related work. To that end, we have released a suite of data products that go beyond the standard contents in the ALMA archive. This release is available online at https://almascience.org/alma-data/lp/DSHARP. It includes (1) CASA scripts and associated python modules used to calibrate and

image the data; (2) fully calibrated continuum and CO measurement sets (visibility datafiles); (3) continuum images and CO channel maps; and (4) some secondary products (radial intensity profiles, and SED data). With this data release and the standard ALMA archive products, the community has the access needed to both reproduce and expand on the efforts detailed in the initial series of DSHARP articles.

## 6. Overview: Initial DSHARP Results

This article has detailed the scientific motivations behind DSHARP, introduced the survey strategy and sample, described the observations and calibration process, and presented the resulting products as part of our data release. It is also the first in a series of articles that explore and analyze the data in more detail. The principal DSHARP conclusions can be summarized as follows:

1. Continuum substructures are ubiquitous in this sample, as can be deduced from Figure 3. Small-scale emission features are found at effectively any disk radius, from 5 au out to more than 150 au.

2. The most common form of these substructures are concentric bright rings and dark gaps. There are no obvious patterns in their distributions or connections to the stellar host properties. There are hints of ring/gap substructures that are obfuscated due to their smaller size scales (relative to the DSHARP resolution) and/or their modest amplitudes with respect to an optically thick background in the inner disk. Measurements of the rings and gaps, as well as a more detailed exploration of their
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Channel maps of the $^{12}$CO J = 2–1 line emission from the HT Lup disk. Beam dimensions are shown in the lower left corner of each panel. The LSRK velocity is marked in the lower right corner of each panel. All the figures are available in the figure set.</figcaption></figure>

(The complete figure set (20 images) is available.)
potential origins and associated issues, are presented by Huang et al. (2018a).

3. While less common, the spiral morphologies identified for a subset of disks in the DSHARP sample are striking. For the cases with apparently single host stars (IM Lup, Elias 27, and WaOph 6), the spiral patterns are complex and appear to be superposed with rings and gaps. Their emission distributions and potential origins are characterized by Huang et al. (2018b).

4. For the two known multiple star systems in the DSHARP sample, HT Lup and AS 205, the disks around the primary stars show clear two-armed spirals and complicated CO distributions that are indicative of strong dynamical interactions. The circumstellar material in these systems is studied by Kurtovic et al. (2018).

5. Azimuthal asymmetries are rare in this sample. Substantial deviations from axisymmetry (or point symmetry for the spirals) are only identified in two cases. The disks around HD 143006 and HD 163296 show small, arc-shaped features in otherwise emission-depleted regions (i.e., beyond the continuum disk edge and in a gap, respectively). The properties and potential origins of these special cases are scrutinized by Pérez et al. (2018) and Isella et al. (2018), respectively.

6. In some cases, the continuum emission can be decomposed into *only* small-scale substructures. The AS 209 disk is a particularly compelling example. Guzmán et al. (2018) quantify its substructures and highlight an important point: there are analogous features lurking in the gas (even as traced by optically thick $^{12}$CO), at radii well beyond the extent of the continuum emission.

7. The ring substructure sizes and amplitudes suggest that these features can be understood as dust trapped in axisymmetric gas pressure bumps. Dullemond et al. (2018) demonstrate this conclusion and derive a lower limit on the strength of the turbulence in the disk. These and other related analyses are guided by a fiducial dust model developed by Birnstiel et al. (2018).

8. A new suite of hydrodynamics simulations by Zhang et al. (2018) suggest that dynamical interactions between low-mass (sub-Jupiter) planets and their local disk material are plausible explanations of the observed ring/gap substructures. Assuming this is the case, those simulations are used to reconstruct the associated planet population in the mass—semimajor axis plane.

There is, of course, much more to learn from the DSHARP data set. Our hope is that this preliminary foray not only provides useful results and motivation for many other studies, but also lays some technical groundwork for designing and calibrating future ALMA surveys of disks at very high angular resolution.

We thank Erica Keller and Tony Remijan for their pipeline calibration efforts, the ALMA staff for pushing through the observations even after the reconfiguration delays, Todd Hunter and Crystal Brogan for some useful technical consultations, Ian Czekala for his advice on measuring stellar parameters, and an anonymous reviewer for thoughtful suggestions. S.A. and J.H. acknowledge support from the National Aeronautics and Space Administration under grant No. 17-XRP17_2-0012 issued through the Exoplanets Research Program. J.H. acknowledges support from the

National Science Foundation Graduate Research Fellowship under grant No. DGE-1144152. L.P. acknowledges support from CONICYT project Basal AFB-170002 and from FCFM/U. de Chile Fondo de Instalación Académica. A.I. acknowledges support from the National Aeronautics and Space Administration under grant No. NNX15AB06G issued through the Origins of Solar Systems program, and from the National Science Foundation under grant No. AST-1715719. C.P.D. acknowledges support by the German Science Foundation (DFG) Research Unit FOR 2634, grants DU 414/22-1 and DU 414/23-1. V.V.G. and J.C. acknowledge support from the National Aeronautics and Space Administration under grant No. 15XRP15_20140 issued through the Exoplanets Research Program. T.B. acknowledges funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme under grant agreement No. 714769. M.B. acknowledges funding from ANR of France under contract number ANR-16-CE31-0013 (Planet Forming disks). Z.Z. and S.Z. acknowledge support from the National Aeronautics and Space Administration through the Astrophysics Theory Program with grant No. NNX17AK40G and the Sloan Research Fellowship. L.R. acknowledges support from the ngVLA Community Studies program, coordinated by the National Radio Astronomy Observatory, which is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This paper makes use of the following ALMA data: ADS/JAO.ALMA #2016.1.00484.L, ADS/JAO.ALMA #2011.0.00531.S, ADS/JAO.ALMA #2012.1.00694.S, ADS/JAO.ALMA #2013.1.00226.S, ADS/JAO.ALMA #2013.1.00366.S, ADS/JAO.ALMA #2013.1.00498.S, ADS/JAO.ALMA #2013.1.00631.S, ADS/JAO.ALMA #2013.1.00798.S, ADS/JAO.ALMA #2015.1.00486.S, ADS/JAO.ALMA #2015.1.00964.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ. This work presents results from the European Space Agency (ESA) space mission Gaia. Gaia data are being processed by the Gaia Data Processing and Analysis Consortium (DPAC). Funding for the DPAC is provided by national institutions, in particular, the institutions participating in the Gaia MultiLateral Agreement (MLA). The Gaia mission website is https://www.cosmos.esa.int/gaia. The Gaia archive website is https://archives.esac.esa.int/gaia.

**Facility:** ALMA.

**Software:** CASA (McMullin et al. 2007), Numpy (Van Der Walt et al. 2011), Matplotlib (Hunter 2007), Astropy (Astropy Collaboration et al. 2018), ScottiePippen (Czekala et al. 2016, https://github.com/iancze/ScottiePippen).

## ORCID iDs

Sean M. Andrews † https://orcid.org/0000-0003-2253-2270

Jane Huang † https://orcid.org/0000-0001-6947-6072

Laura M. Pérez † https://orcid.org/0000-0002-1199-9564

Andrea Isella † https://orcid.org/0000-0001-8061-2207

Cornelis P. Dullemond † https://orcid.org/0000-0002-7078-5910

Viviana V. Guzmán † https://orcid.org/0000-0003-4784-3040

John M. Carpenter † https://orcid.org/0000-0003-2251-0602

David J. Wilner † https://orcid.org/0000-0003-1526-7587
Shangjia Zhang https://orcid.org/0000-0002-8537-9114

Zhaohuan Zhu https://orcid.org/0000-0003-3616-6822

Tilman Birnstiel https://orcid.org/0000-0002-1899-8783

Xue-Ning Bai https://orcid.org/0000-0001-6906-9549

Myriam Benisty https://orcid.org/0000-0002-7695-7605

Karin I. Öberg https://orcid.org/0000-0001-8798-1347

## References

Adachi, I., Hayashi, C., & Nakazawa, K. 1976, PThPh, 56, 1756

Akiyama, E., Hasegawa, Y., Hayashi, M., & Iguchi, S. 2016, ApJ, 818, 158

Alcalá, J. M., Manara, C. F., Natta, A., et al. 2017, A&A, 600, A20

Alibert, Y., Mordasini, C., Benz, W., & Winisdoerffer, C. 2005, A&A, 434, 343

Allard, F., Guillot, T., Ludwig, H.-G., et al. 2003, in IAU Symp. 211, Brown Dwarfs, ed. E. Martín (San Francisco, CA: ASP), 325

Allard, F., Homeier, D., & Freytag, B. 2011, in ASP Conf. Ser. 448, 16th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, ed. C. Johns-Krull, M. K. Browning, & A. A. West (San Francisco, CA: ASP), 91

ALMA Partnership, Brogan, C. L., Pérez, L. M., et al. 2015, ApJL, 808, L3

Andre, P., & Montmerle, T. 1994, ApJ, 420, 837

Andrews, S. M., Rosenfeld, K. A., Kraus, A. L., & Wilner, D. J. 2013, ApJ, 771, 129

Andrews, S. M., & Williams, J. P. 2005, ApJ, 631, 1134

Andrews, S. M., & Williams, J. P. 2007, ApJ, 671, 1800

Andrews, S. M., Wilner, D. J., Espaillat, C., et al. 2011, ApJ, 732, 42

Andrews, S. M., Wilner, D. J., Hughes, A. M., et al. 2012, ApJ, 744, 162

Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C., & Dullemond, C. P. 2009, ApJ, 700, 1502

Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C., & Dullemond, C. P. 2010, ApJ, 723, 1241

Andrews, S. M., Wilner, D. J., Zhu, Z., et al. 2016, ApJL, 820, L40

Andrews, S. M., Wilner, D. J., Zhu, Z., et al. 2018, ApJL, 820, L40

Ansdell, M., Williams, J. P., Trapman, L., et al. 2018, ApJ, 859, 21

Ansdell, M., Williams, J. P., van der Marel, N., et al. 2016, ApJ, 828, 46

Armitage, P. J., Eisner, J. A., & Simon, J. B. 2016, ApJL, 828, L2

Astraatmadja, T. L., & Bailer-Jones, C. A. L. 2016, ApJ, 832, 137

Astropy Collaboration, Price-Whelan, A. M., Sipőcz, B. M., et al. 2018, AJ, 156, 123

Avenhaus, H., Quanz, S. P., Garufi, A., et al. 2018, ApJ, 863, 44

Avenhaus, H., Quanz, S. P., Schmid, H. M., et al. 2014, ApJ, 781, 87

Bai, X.-N., & Stone, J. M. 2014, ApJ, 796, 31

Barenfeld, S. A., Carpenter, J. M., Ricci, L., & Isella, A. 2016, ApJ, 827, 142

Barenfeld, S. A., Carpenter, J. M., Ricci, L., & Isella, A. 2017, ApJ, 851, 85

Barge, P., & Sommeria, J. 1995, A&A, 295, L1

Benisty, M., Juhasz, A., Boccaletti, A., et al. 2015, A&A, 578, L6

Benz, W., Ida, S., Alibert, Y., Lin, D., & Mordasini, C. 2014, in Protostars and Planets VI, ed. H. Beuther et al. (Tucson, AZ: Univ. Arizona Press), 691

Béthune, W., Lesur, G., & Ferreira, J. 2017, A&A, 600, A75

Birnstiel, T., & Andrews, S. M. 2014, ApJ, 780, 153

Birnstiel, T., Dullemond, C. P., Zhu, Z., et al. 2018, ApJL, 869, L45

Birnstiel, T., Fang, M., & Johansen, A. 2016, SSRv, 205, 41

Boehler, Y., Ricci, L., Weaver, E., et al. 2018, ApJ, 853, 162

Bouvier, J., & Appenzeller, I. 1992, A&AS, 92, 481

Brauer, F., Dullemond, C. P., & Henning, T. 2008, A&A, 480, 859

Brauer, F., Dullemond, C. P., Johansen, A., et al. 2007, A&A, 469, 1169

Bryden, G., Chen, X., Lin, D. N. C., Nelson, R. P., & Papaloizou, J. C. B. 1999, ApJ, 514, 344

Calvet, N., D'Alessio, P., Hartmann, L., et al. 2002, ApJ, 568, 1008

Cambrésy, L. 1999, A&A, 345, 965

Carpenter, J. M., Bouwman, J., Silverstone, M. D., et al. 2008, ApJS, 179, 423

Casassus, S., Avenhaus, H., Pérez, S., et al. 2018, MNRAS, 477, 5104

Casassus, S., van der Plas, G., Sebastain Preze, M., et al. 2013, Natur, 493, 191

Choi, J., Dotter, A., Conroy, C., et al. 2016, ApJ, 823, 102

Christiaens, V., Casassus, S., Perez, S., van der Plas, G., & Ménard, F. 2014, ApJL, 785, L12

Cieza, L. A., Casassus, S., Pérez, S., et al. 2017, ApJL, 851, L23

Cieza, L. A., Casassus, S., Tobin, J., et al. 2016, Natur, 535, 258

Cieza, L. A., Ruíz-Rodríguez, D., Hales, A., et al. 2019, MNRAS, 482, 698

Cleeveres, L. I., Bergin, E. A., Öberg, K. I., et al. 2017, ApJL, 843, L3

Cleeveres, L. I., Öberg, K. I., Wilner, D. J., et al. 2016, ApJ, 832, 110

Coleman, G. A. L., & Nelson, R. P. 2016, MNRAS, 460, 2779

Comerón, F. 2008, in Handbook of Star Forming Regions, Volume II: The Southern Sky, ed. B. Reipurth (San Francisco, CA: ASP), 295

Cornwell, T. J. 2008, ISTSP, 2, 793

Cox, E. G., Harris, R. J., Looney, L. W., et al. 2017, ApJ, 851, 83

Czekala, I., Andrews, S. M., Torres, G., et al. 2016, ApJ, 818, 156

de Boer, J., Salter, G., Benisty, M., et al. 2016, A&A, 595, A114

de Geus, E. J., de Zeeuw, P. T., & Lub, J. 1989, A&A, 216, 44

de Gregorio-Monsalvo, I., Ménard, F., Dent, W., et al. 2013, A&A, 557, A133

Dent, W. R. F., Matthews, H. E., & Ward-Thompson, D. 1998, MNRAS, 301, 1049

Dipierro, G., Pinilla, P., Lodato, G., & Testi, L. 2015, MNRAS, 451, 974

Dipierro, G., Ricci, L., Pérez, L., et al. 2018, MNRAS, 475, 5296

Dong, R., Liu, S.-y., Eisner, J., et al. 2018, ApJ, 860, 124

Dullemond, C. P., Birnstiel, T., Huang, J., et al. 2018, ApJL, 869, L46

Dullemond, C. P., & Penzlin, A. B. T. 2018, A&A, 609, A50

Dzyurkevich, N., Turner, N. J., Henning, T., & Kley, W. 2013, ApJ, 765, 114

Eisner, J. A., Hillenbrand, L. A., White, R. J., Akeson, R. L., & Sargent, A. I. 2005, ApJ, 623, 952

Estrada, P. R., Cuzzi, J. N., & Morgan, D. A. 2016, ApJ, 818, 200

Evans, N. J., II, Dunham, M. M., Jørgensen, J. K., et al. 2009, ApJS, 181, 321

Facchini, S., Birnstiel, T., Bruderer, S., & van Dishoeck, E. F. 2017, A&A, 605, A16

Fairlamb, J. R., Oudmaijer, R. D., Mendigutía, I., Ilee, J. D., & van den Ancker, M. E. 2015, MNRAS, 453, 976

Fedele, D., Carney, M., Hogerheijde, M. R., et al. 2017, A&A, 600, A72

Fedele, D., Tazzari, M., Booth, R., et al. 2018, A&A, 610, A24

Flaherty, K. M., Hughes, A. M., Rosenfeld, K. A., et al. 2015, ApJ, 813, 99

Flock, M., Ruge, J. P., Dzyurkevich, N., et al. 2015, A&A, 574, A68

Fung, J., Shi, J.-M., & Chiang, E. 2014, ApJ, 782, 88

Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, A&A, 616, A1

Ginski, C., Stolker, T., Pinilla, P., et al. 2016, A&A, 595, A112

Goodman, J., Narayan, R., & Goldreich, P. 1987, MNRAS, 225, 695

Grady, C. A., Muto, T., Hashimoto, J., et al. 2013, ApJ, 762, 48

Grankin, K. N., Melnikov, S. Y., Bouvier, J., Herbst, W., & Shevchenko, V. S. 2007, A&A, 461, 183

Gras-Velázquez, À., & Ray, T. P. 2005, A&A, 443, 541

Greaves, J. S., & Rice, W. K. M. 2010, MNRAS, 407, 1981

Guilloteau, S., Dutrey, A., Piétu, V., & Boehler, Y. 2011, A&A, 529, A105

Guzmán, V. V., Huang, J., Andrews, S. M., et al. 2018, ApJ, 869, L48

Harris, R. J., Andrews, S. M., Wilner, D. J., & Kraus, A. L. 2012, ApJ, 751, 115

Henning, T., Burkert, A., Launhardt, R., Leinert, C., & Stecklum, B. 1998, A&A, 336, 565

Herbig, G. H., & Bell, K. R. 1988, Third Catalog of Emission-Line Stars of the Orion Population, Vol. 3 (Santa Cruz, CA: Lick Observatory)

Herbst, W., Herbst, D. K., Grossman, E. J., & Weinstein, D. 1994, AJ, 108, 1906

Huang, J., Andrews, S. M., Dullemond, C. P., et al. 2018a, ApJL, 869, L42

Huang, J., Andrews, S. M., Pérez, L. M., et al. 2018b, ApJL, 869, L43

Huang, J., Andrews, S. M., Cleeveres, L. I., et al. 2018c, ApJ, 852, 122

Huang, J., Öberg, K. I., & Andrews, S. M. 2016, ApJ, 823, L18

Huang, J., Öberg, K. I., Qi, C., et al. 2017, ApJ, 835, 231

Hughes, J., Hartigan, P., Krautter, J., & Kelemen, J. 1994, AJ, 108, 1071

Hunter, J. D. 2007, CSE, 9, 90

Ida, S., & Lin, D. N. C. 2004, ApJ, 604, 388

Ida, S., & Lin, D. N. C. 2008, ApJ, 673, 487

Isella, A., Carpenter, J. M., & Sargent, A. I. 2009, ApJ, 701, 260

Isella, A., Carpenter, J. M., & Sargent, A. I. 2010, ApJ, 714, 1746

Isella, A., Guidi, G., Testi, L., et al. 2016, PhRvL, 117, 251101

Isella, A., Huang, J., Andrews, S. M., et al. 2018, ApJL, 869, L49

Isella, A., Pérez, L. M., Carpenter, J. M., et al. 2013, ApJ, 775, 30

Isella, A., Testi, L., Natta, A., et al. 2007, A&A, 469, 213

Ishihara, D., Onaka, T., Kataza, H., et al. 2010, A&A, 514, A1

Jensen, E. L. N., Mathieu, R. D., & Fuller, G. A. 1994, ApJL, 429, L29

Johansen, A., Blum, J., Tanaka, H., et al. 2014, in Protostars & Planets VI, ed. H. Beuther et al. (Tucson: Univ. Arizona Press) in press, (arXiv:1402.1344)

Johansen, A., Youdin, A., & Klahr, H. 2009, ApJ, 697, 1269

Kenyon, S. J., & Hartmann, L. 1987, ApJ, 323, 714

Kley, W., & Nelson, R. P. 2012, ARA&A, 50, 211

Kurtovic, N., Pérez, L. M., Benisty, M., et al. 2018, ApJL, 869, L44

Lin, D. N. C., & Papaloizou, J. C. B. 1993, in Protostars and Planets III, ed. E. H. Levy, J. I. Lunine (Tucson, AZ: Univ. Arizona Press), 749

Loinard, L., Torres, R. M., Mioduszewski, A. J., & Rodríguez, L. F. 2008, ApJL, 675, L29

Lommen, D., Maddison, S. T., Wright, C. M., et al. 2009, A&A, 495, 869

Lommen, D., Wright, C. M., Maddison, S. T., et al. 2007, A&A, 462, 211
Long, F., Pinilla, P., Herczeg, G. J., et al. 2018, ApJ, 869, 17

Loomis, R. A., Öberg, K. I., Andrews, S. M., & MacGregor, M. A. 2017, ApJ, 840, 23

Luhman, K. L., Herrmann, K. A., Mamajek, E. E., Esplin, T. L., & Pecaut, M. J. 2018, AJ, 156, 76

Luhman, K. L., & Mamajek, E. E. 2012, ApJ, 758, 31

Luhman, K. L., & Rieke, G. H. 1999, ApJ, 525, 440

Lyra, W., Turner, N. J., & McNally, C. P. 2015, A&A, 574, A10

Manara, C. F., Rosotti, G., Testi, L., et al. 2016, A&A, 591, L3

Mannings, V., & Emerson, J. P. 1994, MNRAS, 267, 361

Mannings, V., & Sargent, A. I. 1997, ApJ, 490, 792

Marino, S., Perez, S., & Casassus, S. 2015, ApJL, 798, L44

McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, in ASP Conf. Ser. 376, Astronomical Data Analysis Software and Systems XVI, ed. R. A. Shaw, F. Hill, & D. J. Bell (San Francisco, CA: ASP), 127

Mendigutía, I., Mora, A., Montesinos, B., et al. 2012, A&A, 543, A59

Menu, J., van Boekel, R., Henning, T., et al. 2014, A&A, 564, A93

Merín, B., Jørgensen, J., Spezzi, L., et al. 2008, ApJS, 177, 551

Mohanty, S., Greaves, J., Mortlock, D., et al. 2013, ApJ, 773, 168

Mordasini, C., Alibert, Y., & Benz, W. 2009, A&A, 501, 1139

Mulders, G. D., Pascucci, I., Manara, C. F., et al. 2017, ApJ, 847, 31

Muto, T., Grady, C. A., Hashimoto, J., et al. 2012, ApJL, 748, L22

Muzerolle, J., Hartmann, L., & Calvet, N. 1998, AJ, 116, 2965

Najita, J. R., & Kenyon, S. J. 2014, MNRAS, 445, 3315

Nakagawa, Y., Sekiya, M., & Hayashi, C. 1986, Icar, 67, 375

Natta, A., Testi, L., Neri, R., Shepherd, D. S., & Wilner, D. J. 2004, A&A, 416, 179

Natta, A., Testi, L., & Randich, S. 2006, A&A, 452, 245

Nixon, C. J., King, A. R., & Pringle, J. E. 2018, MNRAS, 477, 3273

Nuernberger, D., Chini, R., & Zinnecker, H. 1997, A&A, 324, 1036

Öberg, K. I., Furuya, K., Loomis, R., et al. 2015, ApJ, 810, 112

Öberg, K. I., Qi, C., Fogel, J. K. J., et al. 2011, ApJ, 734, 98

Okuzumi, S., Tanaka, H., Kobayashi, H., & Wada, K. 2012, ApJ, 752, 106

Oudmaijer, R. D., Palacios, J., Eiroa, C., et al. 2001, A&A, 379, 564

Paardekooper, S.-J., & Mellema, G. 2006, A&A, 453, 1129

Padgett, D. L., Cieza, L., Stapelfeldt, K. R., et al. 2006, ApJ, 645, 1283

Panić, O., Hogerheijde, M. R., Wilner, D., & Qi, C. 2009, A&A, 501, 269

Pérez, L. M., Benisty, M., Andrews, S. M., et al. 2018, ApJL, 869, L50

Pérez, L. M., Carpenter, J. M., Andrews, S. M., et al. 2016, Sci, 353, 1519

Pérez, L. M., Carpenter, J. M., Chandler, C. J., et al. 2012, ApJL, 760, L17

Pérez, L. M., Isella, A., Carpenter, J. M., & Chandler, C. J. 2014, ApJL, 783, L13

Pérez, L. M., Isella, A., Carpenter, J. M., & Chandler, C. J. 2015, ApJL, 783, L13

Pinilla, P., Benisty, M., & Birnstiel, T. 2012a, A&A, 545, A81

Pinilla, P., Birnstiel, T., Ricci, L., et al. 2012b, A&A, 538, A114

Pinilla, P., Pohl, A., Stammler, S. M., & Birnstiel, T. 2017, ApJ, 845, 68

Pinilla, P., Tazzari, M., Pascucci, I., et al. 2018, ApJ, 859, 32

Pinte, C., Ménard, F., Duchêne, G., et al. 2018, A&A, 609, A47

Pinte, C., Padgett, D. L., Ménard, F., et al. 2008, A&A, 489, 633

Preibisch, T., & Mamajek, E. 2008, in Handbook of Star Forming Regions, Volume II: The Southern Sky, ed. B. Reipurth (San Francisco, CA: ASP), 235

Qi, C., Öberg, K. I., Andrews, S. M., et al. 2015, ApJ, 813, 128

Quanz, S. P., Avenhaus, H., Buenzli, E., et al. 2013, ApJL, 766, L2

Rapson, V. A., Kastner, J. H., Andrews, S. M., et al. 2015, ApJL, 803, L10

Ribas, Á, Espaillat, C. C., Macías, E., et al. 2017, ApJ, 849, 63

Ricci, L., Testi, L., Natta, A., & Brooks, K. J. 2010, A&A, 521, A66

Rice, W. K. M., Armitage, P. J., Wood, K., & Lodato, G. 2006, MNRAS, 373, 1619

Rigliaco, E., Pascucci, I., Duchene, G., et al. 2015, ApJ, 801, 31

Roccatagliata, V., Henning, T., Wolf, S., et al. 2009, A&A, 497, 409

Rosenfeld, K. A., Andrews, S. M., Wilner, D. J., Kastner, J. H., & McClure, M. K. 2013, ApJ, 775, 136

Rosenfeld, K. A., Chiang, E., & Andrews, S. M. 2014, ApJ, 782, 62

Rosenfeld, K. A., Qi, C., Andrews, S. M., et al. 2012, ApJ, 757, 129

Salyk, C., Herczeg, G. J., Brown, J. M., et al. 2013, ApJ, 769, 21

Salyk, C., Pontoppidan, K., Corder, S., et al. 2014, ApJ, 792, 68

Sandell, G., Weintraub, D. A., & Hamidouche, M. 2011, ApJ, 727, 26

Simon, J. B., & Armitage, P. J. 2014, ApJ, 784, 15

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163

Skrutskie, M. F., Dutkevitch, D., Strom, S. E., et al. 1990, AJ, 99, 1187

Soderblom, D. R., Hillenbrand, L. A., Jeffries, R. D., Mamajek, E. E., & Naylor, T. 2014, in Protostars and Planets VI, ed. H. Beuther et al. (Tucson, AZ: Univ. Arizona Press), 219

Stammler, S. M., Birnstiel, T., Panić, O., Dullemond, C. P., & Dominik, C. 2017, A&A, 600, A140

Stanke, T., Smith, M. D., Gredel, R., & Khanzadyan, T. 2006, A&A, 447, 609

Strom, K. M., Strom, S. E., Edwards, S., Cabrit, S., & Skrutskie, M. F. 1989, AJ, 97, 1451

Suriano, S. S., Li, Z.-Y., Krasnopolsky, R., & Shang, H. 2018, MNRAS, 477, 1239

Sylvester, R. J., Skinner, C. J., Barlow, M. J., & Mannings, V. 1996, MNRAS, 279, 915

Takeuchi, T., & Lin, D. N. C. 2002, ApJ, 581, 1344

Takeuchi, T., & Lin, D. N. C. 2005, ApJ, 623, 482

Tazzari, M., Testi, L., Ercolano, B., et al. 2016, A&A, 588, A53

Tazzari, M., Testi, L., Natta, A., et al. 2017, A&A, 606, 88

Testi, L., Birnstiel, T., Ricci, L., et al. 2014, in Protostars and Planets VI, ed. H. Beuther et al. (Tucson, AZ: Univ. Arizona Press), 339

Tripathi, A., Andrews, S. M., Birnstiel, T., et al. 2018, ApJ, 861, 64

Tripathi, A., Andrews, S. M., Birnstiel, T., & Wilner, D. J. 2017, ApJ, 845, 44

Ubach, C., Maddison, S. T., Wright, C. M., et al. 2017, MNRAS, 466, 4083

van Boekel, R., Henning, T., Menu, J., et al. 2017, ApJ, 837, 132

van der Marel, N., van Dishoeck, E. F., Bruderer, S., et al. 2013, Sci, 340, 1199

van der Marel, N., Williams, J. P., Ansdell, M., et al. 2018, ApJ, 854, 177

van der Plas, G., Wright, C. M., Ménard, F., et al. 2017, A&A, 597, A32

Van Der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, arXiv:1102.1523

van Terwisga, S. E., van Dishoeck, E. F., Ansdell, M., et al. 2018, A&A, 616, A88

Vrba, F. J., Chugainov, P. F., Weaver, W. B., & Stauffer, J. S. 1993, AJ, 106, 1608

Weidenschilling, S. J. 1977, MNRAS, 180, 57

Whipple, F. L. 1972, in From Plasma to Planet, ed. A. Elvius (New York: Wiley Interscience), 211

Wilking, B. A., Gagné, M., & Allen, L. E. 2008, in Handbook of Star Forming Regions, Volume II: The Southern Sky, ed. B. Reipurth (San Francisco, CA: ASP), 351

Wilking, B. A., Meyer, M. R., Robinson, J. G., & Greene, T. P. 2005, AJ, 130, 1733

Wilner, D. J., & Welch, W. J. 1994, ApJ, 427, 898

Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868

Youdin, A. N., & Goodman, J. 2005, ApJ, 620, 459

Youdin, A. N., & Shu, F. H. 2002, ApJ, 580, 494

Zhang, K., Bergin, E. A., Blake, G. A., et al. 2016, ApJL, 818, L16

Zhang, K., Isella, A., Carpenter, J. M., & Blake, G. A. 2014, ApJ, 791, 42

Zhang, S., Zhu, Z., Huang, J., et al. 2018, ApJL, 869, L47

Zhu, Z., Nelson, R. P., Dong, R., Espaillat, C., & Hartmann, L. 2012, ApJ, 755, 6