We present a sensitive, multiwavelength submillimeter continuum survey of 153 young stellar objects in the Taurus-Auriga star formation region. The submillimeter detection rate is 61% to a completeness limit of ∼10 mJy (3-σ) at 850 µm. The inferred circumstellar disk masses are log-normally distributed with a mean mass of ∼ 5 × 10 −3 M ⊙ and a large dispersion (0.5 dex). Roughly one third of the submillimeter sources have disk masses larger than the minimal nebula from which the solar system formed. The median disk to star mass ratio is 0.5%. The empirical behavior of the submillimeter continuum is best described as F ν ∝ ν 2.0±0.5 between 350 µm and 1.3 mm, which we argue is due to the combined effects of the fraction of optically thick emission and a flatter frequency behavior of the opacity compared to the interstellar medium. This latter effect could be due to a substantial population of large dust grains, which presumably would have grown through collisional agglomeration. In this sample, the only stellar property that is correlated with the outer disk is the presence of a companion. We find evidence for significant decreases in submillimeter flux densities, disk masses, and submillimeter continuum slopes along the canonical infrared spectral energy distribution evolution sequence for young stellar objects. The fraction of objects detected in the submillimeter is essentially identical to the fraction with excess nearinfrared emission, suggesting that dust in the inner and outer disk are removed nearly simultaneously.

We present a substantial extension of the millimeter (mm) wave continuum photometry catalog for circumstellar dust disks in the Taurus star-forming region, based on a new "snapshot" λ = 1.3 mm survey with the Submillimeter Array. Combining these new data with measurements in the literature, we construct a mm-wave luminosity distribution, f (L mm ), for Class II disks that is statistically complete for stellar hosts with spectral types earlier than M8.5 and has a 3σ depth of roughly 3 mJy. The resulting census eliminates a longstanding selection bias against disks with late-type hosts, and thereby demonstrates that there is a strong correlation between L mm and the host spectral type. By translating the locations of individual stars in the Hertzsprung-Russell diagram into masses and ages, and adopting a simple conversion between L mm and the disk mass, M d , we confirm that this correlation corresponds to a statistically robust relationship between the masses of dust disks and the stars that host them. A Bayesian regression technique is used to characterize these relationships in the presence of measurement errors, data censoring, and significant intrinsic scatter: the best-fit results indicate a typical 1.3 mm flux density of ∼25 mJy for 1 M hosts and a power-law scaling L mm ∝ M 1.5−2.0 * . We suggest that a reasonable treatment of dust temperature in the conversion from L mm to M d favors an inherently linear M d ∝ M * scaling, with a typical disk-to-star mass ratio of ∼0.2%-0.6%. The measured rms dispersion around this regression curve is ±0.7 dex, suggesting that the combined effects of diverse evolutionary states, dust opacities, and temperatures in these disks imprint a full width at half-maximum range of a factor of ∼40 on the inferred M d (or L mm ) at any given host mass. We argue that this relationship between M d and M * likely represents the origin of the inferred correlation between giant planet frequency and host star mass in the exoplanet population, and provides some basic support for the core accretion model for planet formation. Moreover, we caution that the effects of incompleteness and selection bias must be considered in comparative studies of disk evolution, and illustrate that fact with statistical comparisons of f (L mm ) between the Taurus catalog presented here and incomplete subsamples in the Ophiuchus, IC 348, and Upper Sco young clusters.

Circumstellar disks are thought to experience a rapid "transition" phase in their evolution that can have a considerable impact on the formation and early development of planetary systems. We present new and archival high angular resolution (0. ′′ 3 ≈ 40-75 AU) Submillimeter Array (SMA) observations of the 880 µm (340 GHz) dust continuum emission from 12 such transition disks in nearby star-forming regions. In each case, we directly resolve a dust-depleted disk cavity around the central star. Using two-dimensional Monte Carlo radiative transfer calculations, we interpret these dust disk structures in a homogeneous, parametric model framework by reproducing their SMA continuum visibilities and spectral energy distributions. The cavities in these disks are large (R cav = 15-73 AU) and substantially depleted of small (∼µm-sized) dust grains, although their mass contents are still uncertain. The structures of the remnant material at larger radii are comparable to normal disks. We demonstrate that these large cavities are relatively common among the millimeter-bright disk population, comprising at least 1 in 5 (20%) of the disks in the bright half (and ≥26% of the upper quartile) of the millimeter luminosity (disk mass) distribution. Utilizing these results, we assess some of the physical mechanisms proposed to account for transition disk structures. As has been shown before, photoevaporation models do not produce the large cavity sizes, accretion rates, and disk masses representative of this sample. It would be difficult to achieve a sufficient decrease of the dust optical depths in these cavities by particle growth alone: substantial growth (to meter sizes or beyond) must occur in large (tens of AU) regions of low turbulence without also producing an abundance of small particles. Given those challenges, we suggest instead that the observations are most commensurate with dynamical clearing due to tidal interactions with low-mass companions -young brown dwarfs or giant planets on long-period orbits. Subject headings: circumstellar matter -protoplanetary disks -planet-disk interactions -planets and satellites: formation -submillimeter: planetary systems

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 au 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.

We present the results of a high angular resolution (0. 3 ≈ 40 AU) Submillimeter Array survey of the 345 GHz (870 μm) thermal continuum emission from nine of the brightest, and therefore most massive, circumstellar disks in the ∼1 Myr-old Ophiuchus star-forming region. Using two-dimensional radiative transfer calculations, we simultaneously fit the observed continuum visibilities and broadband spectral energy distribution for each disk with a parametric structure model. Compared to previous millimeter studies, this survey includes significant upgrades in modeling, data quality, and angular resolution that provide improved constraints on key structure parameters, particularly those that characterize the spatial distribution of mass in the disks. In the context of a surface density profile motivated by similarity solutions for viscous accretion disks,, the bestfit models for the sample disks have characteristic radii R c ≈ 20-200 AU, high disk masses M d ≈ 0.005-0.14 M (a sample selection bias), and a narrow range of radial Σ gradients (γ ≈ 0.4-1.0) around a median γ = 0.9. These density structures are used in conjunction with accretion rate estimates from the literature to help characterize the viscous evolution of the disk material. Using the standard prescription for disk viscosities, those combined constraints indicate that α ≈ 0.0005-0.08. Three of the sample disks show large (R ≈ 20-40 AU) central cavities in their continuum emission morphologies, marking extensive zones where dust has been physically removed and/or has significantly diminished opacities. Based on the current requirements of planet formation models, these emission cavities and the structure constraints for the sample as a whole suggest that these young disks may eventually produce planetary systems, and have perhaps already started.