Search

The accretion/ejection processes in T-Tauri stars are fundamental to their physical evolution, while also impacting the properties and evolution of the circumstellar material at a time when planet formation takes place. To this date, characterization of ongoing accretion processes in stellar pairs at 5-50\,au scales has been challenging, high angular resolution spectrographs are required to extract the spectral features of each component. We present the analysis of spectroscopic observations of the tight (160mas, 25au) T-Tauri system HT Lup A/B, obtained with MUSE at VLT in March and July of 2021. We focus on constraining the accretion/ejection processes and variability of the secondary component HT Lup B, by searching for accretion tracers applying High-Resolution Spectral Differential Imaging techniques. We retrieve strong (SNR $>$ 5) $H\alpha, H\beta$ and [OI]$\lambda6300$ emission in both epochs. The $H\alpha$ and $H\beta$ line fluxes showcase high variability, with variations up to 400-500\% between epochs. The fluxes are consistent with accretion rates of $8\times10^{-9} M_\odot \, yr^{-1}$ and $2\times10^{-9} M_\odot \, yr^{-1}$ for the first and second epoch, respectively. We attribute the increased accretion activity during the first night to a "burst" like event, followed by a relaxation period more representative of the common accretion activity of the system. The [OI]$\lambda6300$ line profiles remain relatively similar between epochs and suggest ejection rates on the order of $10^{-9}-10^{-10} M_\odot \, yr^{-1}$, compatible with moderate disk winds emission. Our results also indicate that the accretion processes of HT Lup B are compatible with Classical T Tauri Stars, unlike previous classifications, Comment: 28 pages, 13 fgures, Accepted by ApJ

The physical origin of the large cavities observed in transition disks is to date still unclear. Different physical mechanisms (e.g., a companion, dead zones, enhanced grain growth) produce disk cavities of different depth, and the expected spatial distribution of gas and solids in each mechanism is not the same. In this work, we analyze the multi-wavelength interferometric visibilities of dust continuum observations obtained with ALMA and VLA for six transition disks: CQTau, UXTau A, LkCa15, RXJ1615, SR24S, and DMTau, and calculate brightness radial profiles, where diverse emission morphology is revealed at different wavelengths. The multi-wavelength data is used to model the spectral energy distribution and compute constraints on the radial profile of the dust surface density, maximum grain size, and dust temperature in each disk. They are compared with the observational signatures expected from various physical mechanisms responsible for disk cavities. The observational signatures suggest that the cavities observed in the disks around UXTau A, LkCa15, and RXJ1615 could potentially originate from a dust trap created by a companion. Conversely, in the disks around CQTau, SR24S, DMTau, the origin of the cavity remains unclear, although it is compatible with a pressure bump and grain growth within the cavity.

Detecting planet signatures in protoplanetary disks is fundamental to understanding how and where planets form. In this work, we report dust and gas observational hints of planet formation in the disk around 2MASS-J16120668-301027, as part of the ALMA Large Program "AGE-PRO: ALMA survey of Gas Evolution in Protoplanetary disks". The disk was imaged with the Atacama Large Millimeter/submillimeter Array (ALMA) at Band 6 (1.3 mm) in dust continuum emission and four molecular lines: $^{12}$CO(J=2-1), $^{13}$CO(J=2-1), C$^{18}$O(J=2-1), and H$_2$CO(J=3$_{(3,0)}$-2$_{(2,0)}$). Resolved observations of the dust continuum emission (angular resolution of $\sim 150$ mas, 20 au) show a ring-like structure with a peak at $0.57 ^{\prime \prime}$ (75 au), a deep gap with a minimum at 0.24$^{\prime \prime}$ (31 au), an inner disk, a bridge connecting the inner disk and the outer ring, along with a spiral arm structure, and a tentative detection (to $3\sigma$) of a compact emission at the center of the disk gap, with an estimated dust mass of $\sim 2.7-12.9$ Lunar masses. We also detected a kinematic kink (not coincident with any dust substructure) through several $^{12}$CO channel maps (angular resolution $\sim$ 200 mas, 30 au), located at a radius of $\sim 0.875^{\prime \prime}$ (115.6 au). After modeling the $^{12}$CO velocity rotation around the protostar, we identified a tentative rotating-like structure at the kink location with a geometry similar to that of the disk. We discuss potential explanations for the dust and gas substructures observed in the disk, and their potential connection to signatures of planet formation., Comment: 24 pages, 15 figures

Rings and gaps are among the most widely observed forms of substructure in protoplanetary disks. A gap-ring pair may be formed when a planet carves a gap in the disk, which produces a local pressure maximum following the gap that traps inwardly drifting dust grains and appears as a bright ring due to the enhanced dust density. A dust-trapping ring would provide a promising environment for solid growth and possibly planetesimal production via the streaming instability. We present evidence of dust trapping in the bright ring of the planet-hosting disk Elias 2-24, from the analysis of 1.3 mm and 3 mm ALMA observations at high spatial resolution (0.029 arcsec, 4.0 au). We leverage the high spatial resolution to demonstrate that larger grains are more efficiently trapped and place constraints on the local turbulence ($8 \times 10^{-4} < \alpha_\mathrm{turb} < 0.03$) and the gas-to-dust ratio ($\Sigma_g / \Sigma_d < 30$) in the ring. Using a scattering-included marginal probability analysis we measure a total dust disk mass of $M_\mathrm{dust} = 13.8^{+0.7}_{-0.5} \times 10^{-4} \ M_\odot$. We also show that at the orbital radius of the proposed perturber, the gap is cleared of material down to a flux contrast of 10$^{-3}$ of the peak flux in the disk., Comment: 18 pages, 12 figures, accepted to ApJ

We demonstrate the existence of semivortex (SV) solitons, with vorticities $0$ and $1$ in the two components, in a two-dimensional (2D) fermionic spinor system under the action of the Rashba-type spin-orbit coupling in the combination with the Zeeman splitting (ZS). In the ``heavy-atom" approximation, which was previously elaborated for the bosonic system, the usual kinetic energy is neglected, which gives rise to a linear spectrum with a bandgap. The model includes the effective Pauli self-repulsion with power $7/3$, as produced by the density-functional theory of Fermi superfluids. In the general case, the inter-component contact repulsion is included too. We construct a family of gap solitons of the SV type populating the spectral bandgap. A stability region is identified for the SV solitons, by means of systematic simulations, in the parameter plane of the cross-repulsion strength and chemical potential. The stability region agrees with the prediction of the anti-Vakhitov-Kolokolov criterion, which is a relevant necessary stability condition for systems with self-repulsive nonlinearities. We also test the stability of the SV solitons against a sudden change of the ZS strength, which initiates robust oscillations in the spin state of the soliton due to transfer of particles between the system's components., Comment: Chaos, Solitons & Fractals, in press

We report kinematic and thermal signatures associated with the directly imaged protoplanet candidate in the Elias 2-24 disc. Using the DSHARP ALMA observations of the $^{12}$CO J=2-1 line, we show that the disc kinematics are perturbed, with a detached CO emission spot at the location of the planet candidate and traces of spiral wakes, and also that the observed CO emission intensities require local heating. While the foreground extinction hides the velocity channels associated with the planet, preventing a planet mass estimate, the level of gas heating implied by the CO emission indicates the presence of a warm, embedded giant planet. Comparison with models show this could either be a $\gtrsim 5$M$_\mathrm{Jup}$, or a lower mass ( $\gtrsim 2$M$_\mathrm{Jup}$) but accreting proto-planet., Comment: 6 pages, 3 figures, accepted for publication in MNRAS Letters

Transition disks are protoplanetary disks with inner cavities possibly cleared by massive companions, which makes them prime targets to observe at high resolution to map their velocity structure. We present ALMA Band 6 dust and gas observations of the transition disk around RXJ1604.3-2130 A, known to feature nearly symmetric shadows in scattered light. We studied the $^{12}$CO line channel maps and moment maps of the line of sight velocity and peak intensity. We fitted a Keplerian model of the channel-by-channel emission to study line profile differences and produced deprojected radial profiles for all velocity components. The $^{12}$CO emission shows a cavity inwards of $\sim$56 au and within the dust continuum ring at 81 au. The azimuthal brightness variations in the $^{12}$CO line and dust continuum are broadly aligned with the shadows detected in scattered-light observations. We find a strong localized non-Keplerian feature toward the west within the continuum ring (at $R=41\pm 10$ au and $PA=280\pm 2 ^\circ$). A tightly wound spiral is also detected which extends over 300$^\circ$ in azimuth, possibly connected to the localized non-Keplerian feature. Finally, a bending of the iso-velocity contours within the gas cavity indicates a highly perturbed inner region, possibly related to the presence of a misaligned inner disk. While broadly aligned with the scattered-light shadows, the localized non-Keplerian feature cannot be solely due to changes in temperature. Instead, we interpret the kinematical feature as tracing a massive companion located at the edge of the dust continuum ring. We speculate that the spiral is caused by buoyancy resonances driven by planet-disk interactions. However, this potential planet at $\sim$41 au cannot explain the gas-depleted cavity, the low accretion rate, and the misaligned inner disk, which suggests the presence of another companion closer in., Comment: 13 pages, 14 figures. To be published in Astronomy & Astrophysics Letters

High spatial resolution observations of CO isotopologue line emission in protoplanetary disks at mid-inclinations (${\approx}$30-75{\deg}) allow us to characterize the gas structure in detail, including radial and vertical substructures, emission surface heights and their dependencies on source characteristics, and disk temperature profiles. By combining observations of a suite of CO isotopologues, we can map the 2D (r, z) disk structure from the disk upper atmosphere, as traced by CO, to near the midplane, as probed by less abundant isotopologues. Here, we present high angular resolution (${\lesssim}$0."1 to ${\approx}$0."2; ${\approx}$15-30 au) observations of CO, $^{13}$CO, and C$^{18}$O in either or both J=2-1 and J=3-2 lines in the transition disks around DM Tau, Sz 91, LkCa 15, and HD 34282. We derived line emission surfaces in CO for all disks and in $^{13}$CO for the DM Tau and LkCa 15 disks. With these observations, we do not resolve the vertical structure of C$^{18}$O in any disk, which is instead consistent with C$^{18}$O emission originating from the midplane. Both the J=2-1 and J=3-2 lines show similar heights. Using the derived emission surfaces, we computed radial and vertical gas temperature distributions for each disk, including empirical temperature models for the DM Tau and LkCa 15 disks. After combining our sample with literature sources, we find that $^{13}$CO line emitting heights are also tentatively linked with source characteristics, e.g., stellar host mass, gas temperature, disk size, and show steeper trends than seen in CO emission surfaces., Comment: 31 pages, 15 figures, accepted for publication in ApJ. Image cubes available at https://zenodo.org/record/7430257

We report the discovery of a circumplanetary disk (CPD) candidate embedded in the circumstellar disk of the T Tauri star AS 209 at a radial distance of about 200 au (on-sky separation of 1."4 from the star at a position angle of $161^\circ$), isolated via $^{13}$CO $J=2-1$ emission. This is the first instance of CPD detection via gaseous emission capable of tracing the overall CPD mass. The CPD is spatially unresolved with a $117\times82$ mas beam and manifests as a point source in $^{13}$CO, indicating that its diameter is $\lesssim14$ au. The CPD is embedded within an annular gap in the circumstellar disk previously identified using $^{12}$CO and near-infrared scattered light observations, and is associated with localized velocity perturbations in $^{12}$CO. The coincidence of these features suggests that they have a common origin: an embedded giant planet. We use the $^{13}$CO intensity to constrain the CPD gas temperature and mass. We find that the CPD temperature is $\gtrsim35$ K, higher than the circumstellar disk temperature at the radial location of the CPD, 22 K, suggesting that heating sources localized to the CPD must be present. The CPD gas mass is $\gtrsim 0.095 M_{\rm Jup} \simeq 30 M_{\rm Earth}$ adopting a standard $^{13}$CO abundance. From the non-detection of millimeter continuum emission at the location of the CPD ($3\sigma$ flux density $\lesssim26.4~\mu$Jy), we infer that the CPD dust mass is $\lesssim 0.027 M_{\rm Earth} \simeq 2.2$ lunar masses, indicating a low dust-to-gas mass ratio of $\lesssim9\times10^{-4}$. We discuss the formation mechanism of the CPD-hosting giant planet on a wide orbit in the framework of gravitational instability and pebble accretion., Comment: Accepted for publication in the ApJ Letters (July 7, 2022), 19 pages, 13 figures, interactive figures (Figure 7, 8, 9) are available at http://jaehanbae.com/as209/

We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of transitional, pre-transitional, and full disks. Our approach explicitly did not favor well-known, "famous" disks or those observed by ALMA, resulting in a less-biased sample suitable to probe the major stages of disk evolution during planet formation. Our optimized data reduction allowed polarized flux as low as 0.002% of the stellar light to be detected, and we report polarized scattered light around 80% of our targets. We detected point-like companions for 47% of the targets, including 3 brown dwarfs (2 confirmed, 1 new), and a new super-Jupiter mass candidate around V1295 Aql. We searched for correlations between the polarized flux and system parameters, finding a few clear trends: presence of a companion drastically reduces the polarized flux levels, far-IR excess correlates with polarized flux for non-binary systems, and systems hosting disks with ring structures have stellar masses $<$ 3 Msun. Our sample also included four hot, dusty "FS CMa" systems and we detected large-scale ($>100$ au) scattered light around each, signs of extreme youth for these enigmatic systems. Science-ready images are publicly available through multiple distribution channels using a new FITS file standard jointly developed with members of the VLT/SPHERE team., Comment: 51 pages, 31 figures, 7 tables, accepted to AJ

High spatial resolution CO observations of mid-inclination (30-75{\deg}) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure of the local gas temperature. Our analysis of ALMA archival data yields CO emission surfaces, dynamically-constrained stellar host masses, and disk atmosphere gas temperatures for the disks around: HD 142666, MY Lup, V4046 Sgr, HD 100546, GW Lup, WaOph 6, DoAr 25, Sz 91, CI Tau, and DM Tau. These sources span a wide range in stellar masses (0.50-2.10 M$_{\odot}$), ages (${\sim}$0.3-23 Myr), and CO gas radial emission extents (${\approx}$200-1000 au). This sample nearly triples the number of disks with mapped emission surfaces and confirms the wide diversity in line emitting heights ($z/r\approx0.1$ to ${\gtrsim}0.5$) hinted at in previous studies. We compute radial and vertical CO gas temperature distributions for each disk. A few disks show local temperature dips or enhancements, some of which correspond to dust substructures or the proposed locations of embedded planets. Several emission surfaces also show vertical substructures, which all align with rings and gaps in the millimeter dust. Combining our sample with literature sources, we find that CO line emitting heights weakly decline with stellar mass and gas temperature, which, despite large scatter, is consistent with simple scaling relations. We also observe a correlation between CO emission height and disk size, which is due to the flared structure of disks. Overall, CO emission surfaces trace ${\approx}2$-$5\times$ gas pressure scale heights (H$_{\rm{g}}$) and could potentially be calibrated as empirical tracers of H$_{\rm{g}}$., Comment: 31 pages, 17 figures, accepted for publication in ApJ. Image cubes available at https://zenodo.org/record/6410045

Multi-wavelengths high-resolution imaging of protoplanetary disks has revealed the presence of multiple, varied substructures in their dust and gas components which might be signposts of young, forming planetary systems. AB Aurigae bears an emblematic (pre)transitional disk showing spiral structures observed in the inner cavity of the disk in both the sub-millimeter (ALMA; 1.3mm, $^{12}$CO) and near-infrared (SPHERE; 1.5-2.5$\mu$m) wavelengths which have been claimed to arise from dynamical interactions with a massive companion. In this work, we present new deep $K_s$ (2.16$\mu$m) and L' (3.7$\mu$m) band images of AB Aurigae obtained with LMIRCam on the Large Binocular Telescope, aimed for the detection of both planetary companions and extended disk structures. No point source is recovered, in particular at the outer regions of the disk, where a putative candidate ($\rho = 0.681", PA = 7.6^{\circ}$) had been previously claimed. The nature of a second innermost planet candidate ($\rho = 0.16'', PA = 203.9^{\circ}$) can not be investigated by the new data. We are able to derive 5$\sigma$ detection limits in both magnitude and mass for the system, going from 14 \Mjup at 0.3'' (49 au) down to 3-4 \Mjup at 0.6'' (98 au) and beyond, based on the ATMO 2020 evolutionary models. We detect the inner spiral structures (< 0.5'') resolved in both CO and polarimetric H-band observations. We also recover the ring structure of the system at larger separation (0.5-0.7") showing a clear south-east/north-west asymmetry. This structure, observed for the first time at L'-band, remains interior to the dust cavity seen at ALMA, suggesting an efficient dust trapping mechanism at play in the disk., Comment: 13 pages, 5 figures, accepted on ApJ

Constraining dust properties of planet-forming disks via high angular resolution observations is fundamental to understanding how solids are trapped in substructures and how dust growth may be favored or accelerated therein. We use ALMA dust continuum observations of the Molecules with ALMA at Planet-forming Scales (MAPS) disks and explore a large parameter space to constrain the radial distribution of solid mass and maximum grain size in each disk, including or excluding dust scattering. In the nonscattering model, the dust surface density and maximum grain size profiles decrease from the inner disks to the outer disks, with local maxima at the bright ring locations, as expected from dust trapping models. The inferred maximum grain sizes from the inner to outer disks decrease from ~1 cm to 1 mm. For IM Lup, HD 163296, and MWC 480 in the scattering model, two solutions are compatible with their observed inner disk emission: one solution corresponding to a maximum grain size of a few millimeters (similar to the nonscattering model), and the other corresponding to a few hundred micrometer sizes. Based on the estimated Toomre parameter, only IM Lup -- which shows a prominent spiral morphology in millimeter dust -- is found to be gravitationally unstable. The estimated maximum Stokes number in all the disks lies between 0.01 and 0.3, and the estimated turbulence parameters in the rings of AS 209 and HD 163296 are close to the threshold where dust growth is limited by turbulent fragmentation. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.

Small organic molecules, such as C2H, HCN, and H2CO, are tracers of the C, N, and O budget in protoplanetary disks. We present high angular resolution (10-50 au) observations of C2H, HCN, and H2CO lines in five protoplanetary disks from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We derive column density and excitation temperature profiles for HCN and C2H, and find that the HCN emission arises in a temperate (20-30 K) layer in the disk, while C2H is present in relatively warmer (20-60 K) layers. In the case of HD 163296, we find a decrease in column density for HCN and C2H inside one of the dust gaps near 83 au, where a planet has been proposed to be located. We derive H2CO column density profiles assuming temperatures between 20 and 50 K, and find slightly higher column densities in the colder disks around T Tauri stars than around Herbig Ae stars. The H2CO column densities rise near the location of the CO snowline and/or millimeter dust edge, suggesting an efficient release of H2CO ices in the outer disk. Finally, we find that the inner 50 au of these disks are rich in organic species, with abundances relative to water that are similar to cometary values. Comets could therefore deliver water and key organics to future planets in these disks, similar to what might have happened here on Earth. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: 23 pages, 10 figures, 4 tables. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement

Planets form and obtain their compositions in dust and gas-rich disks around young stars, and the outcome of this process is intimately linked to the disk chemical properties. The distributions of molecules across disks regulate the elemental compositions of planets, including C/N/O/S ratios and metallicity (O/H and C/H), as well as access to water and prebiotically relevant organics. Emission from molecules also encodes information on disk ionization levels, temperature structures, kinematics, and gas surface densities, which are all key ingredients of disk evolution and planet formation models. The Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program was designed to expand our understanding of the chemistry of planet formation by exploring disk chemical structures down to 10 au scales. The MAPS program focuses on five disks - around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480 - in which dust substructures are detected and planet formation appears to be ongoing. We observed these disks in 4 spectral setups, which together cover ~50 lines from over 20 different species. This paper introduces the ApJS MAPS Special Issue by presenting an overview of the program motivation, disk sample, observational details, and calibration strategy. We also highlight key results, including discoveries of links between dust, gas, and chemical sub-structures, large reservoirs of nitriles and other organics in the inner disk regions, and elevated C/O ratios across most disks. We discuss how this collection of results is reshaping our view of the chemistry of planet formation., Comment: Accepted for publication in the ApJS MAPS Special Issue. v2 has updated MAPS references and a correction to Fig. 3

Here we present high resolution (15-24 au) observations of CO isotopologue lines from the Molecules with ALMA on Planet-forming Scales (MAPS) ALMA Large Program. Our analysis employs $^{13}$CO and C$^{18}$O ($J$=2-1), (1-0), and C$^{17}$O (1-0) line observations of five protoplanetary disks. We retrieve CO gas density distributions, using three independent methods: (1) a thermo-chemical modeling framework based on the CO data, the broadband spectral energy distribution, and the mm-continuum emission; (2) an empirical temperature distribution based on optically thick CO lines; and (3) a direct fit to the C$^{17}$O hyperfine lines. Results from these methods generally show excellent agreement. The CO gas column density profiles of the five disks show significant variations in the absolute value and the radial shape. Assuming a gas-to-dust mass ratio of 100, all five disks have a global CO-to-H$_2$ abundance of 10-100 times lower than the ISM ratio. The CO gas distributions between 150-400 au match well with models of viscous disks, supporting the long-standing theory. CO gas gaps appear to be correlated with continuum gap locations, but some deep continuum gaps do not have corresponding CO gaps. The relative depths of CO and dust gaps are generally consistent with predictions of planet-disk interactions, but some CO gaps are 5-10 times shallower than predictions based on dust gaps. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. 36 pages, 21 figures, accepted for publication in ApJS

The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a unique opportunity to study the vertical distribution of gas, chemistry, and temperature in the protoplanetary disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. By using the asymmetry of molecular line emission relative to the disk major axis, we infer the emission height ($z$) above the midplane as a function of radius ($r$). Using this method, we measure emitting surfaces for a suite of CO isotopologues, HCN, and C$_2$H. We find that $^{12}$CO emission traces the most elevated regions with $z/r > 0.3$, while emission from the less abundant $^{13}$CO and C$^{18}$O probes deeper into the disk at altitudes of $z/r \lesssim 0.2$. C$_2$H and HCN have lower opacities and SNRs, making surface fitting more difficult, and could only be reliably constrained in AS 209, HD 163296, and MWC 480, with $z/r \lesssim 0.1$, i.e., relatively close to the planet-forming midplanes. We determine peak brightness temperatures of the optically thick CO isotopologues and use these to trace 2D disk temperature structures. Several CO temperature profiles and emission surfaces show dips in temperature or vertical height, some of which are associated with gaps and rings in line and/or continuum emission. These substructures may be due to local changes in CO column density, gas surface density, or gas temperatures, and detailed thermo-chemical models are necessary to better constrain their origins and relate the chemical compositions of elevated disk layers with those of planet-forming material in disk midplanes. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: 31 pages, 20 figures, accepted for publication in ApJS, MAPS cross-references updated

Constraining the distribution of gas and dust in the inner 20 au of protoplanetary disks is difficult. At the same time, this region is thought to be responsible for most planet formation, especially around the water ice line at 3-10 au. Under the assumption that the gas is in a Keplerian disk, we use the exquisite sensitivity of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA large program to construct radial surface brightness profiles with a ~3 au effective resolution for the CO isotopologue J=2-1 lines using the line velocity profile. IM Lup reveals a central depression in 13CO and C18O that is ascribed to a pileup of ~500 $M_\oplus$ of dust in the inner 20 au, leading to a gas-to-dust ratio of around <10. This pileup is consistent with efficient drift of grains ($\gtrsim$ 100 $M_\oplus$ Myr$^{-1}$) and a local gas-to-dust ratio that suggests that the streaming instability could be active. The CO isotopologue emission in the GM Aur disk is consistent with a small (~15 au), strongly depleted gas cavity within the ~40 au dust cavity. The radial surface brightness profiles for both the AS 209 and HD 163296 disks show a local minimum and maximum in the C18O emission at the location of a known dust ring (~14 au) and gap (~10 au), respectively. This indicates that the dust ring has a low gas-to-dust ratio ($>$ 10) and that the dust gap is gas-rich enough to have optically thick C18O., Comment: 22 pages, 14 figures, accepted by ApJS. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. Update with correct references to other MAPS papers

We explore the dynamical structure of the protoplanetary disks surrounding HD 163296 and MWC 480 as part of the Molecules with ALMA at Planet Forming Scales (MAPS) large program. Using the $J = 2-1$ transitions of $^{12}$CO, $^{13}$CO and C$^{18}$O imaged at spatial resolutions of $\sim 0.^{\prime \prime}15$ and with a channel spacing of $200$ ${\rm m\,s^{-1}}$, we find perturbations from Keplerian rotation in the projected velocity fields of both disks ($\lesssim\!5\%$ of the local Keplerian velocity), suggestive of large-scale (10s of au in size), coherent flows. By accounting for the azimuthal dependence on the projection of the velocity field, the velocity fields were decomposed into azimuthally averaged orthogonal components, $v_{\phi}$, $v_r$ and $v_z$. Using the optically thick $^{12}$CO emission as a probe of the gas temperature, local variations of $\approx\! 3$ K ($\approx\! 5 \%$ relative changes) were observed and found to be associated with the kinematic substructures. The MWC 480 disk hosts a suite of tightly wound spiral arms. The spirals arms, in conjunction with the highly localized perturbations in the gas velocity structure (kinematic planetary signatures), indicate a giant planet, $\sim\! 1$ $M_{\rm Jup}$, at a radius of $\approx 245$ au. In the disk of HD 163296, the kinematic substructures were consistent with previous studies of Pinte et al. (2018a) and Teague et al. (2018a) advocating for multiple $\sim\! 1$ $M_{\rm Jup}$ planets embedded in the disk. These results demonstrate that molecular line observations that characterize the dynamical structure of disks can be used to search for the signatures of embedded planets. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: 32 pages, 21 figures, accepted for publication in ApJS, MAPS cross-references updated

The Molecules with ALMA at Planet-forming Scales large program (MAPS LP) surveyed the chemical structures of five protoplanetary disks across more than 40 different spectral lines at high angular resolution (0.15" and 0.30" beams for Bands 6 and 3, respectively) and sensitivity (spanning 0.3 - 1.3 mJy/beam and 0.4 - 1.9 mJy/beam for Bands 6 and 3, respectively). In this article, we describe our multi-stage workflow -- built around the CASA tclean image deconvolution procedure -- that we used to generate the core data product of the MAPS LP: the position-position-velocity image cubes for each spectral line. Owing to the expansive nature of the survey, we encountered a range of imaging challenges; some are familiar to the sub-mm protoplanetary disk community, like the benefits of using an accurate CLEAN mask, and others less well-known, like the incorrect default flux scaling of the CLEAN residual map first described in Jorsater & van Moorsel 1995 (the "JvM effect"). We distill lessons learned into recommended workflows for synthesizing image cubes of molecular emission. In particular, we describe how to produce image cubes with accurate fluxes via the "JvM correction," a procedure that is generally applicable to any image synthesized via CLEAN deconvolution but is especially critical for low S/N emission. We further explain how we used visibility tapering to promote a common, fiducial beam size and contextualize the interpretation of signal to noise ratio when detecting molecular emission from protoplanetary disks. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: ApJS accepted. Part of the MAPS ALMA large program series: http://www.alma-maps.info/. Updated bibliography with MAPS LP arXiv references

PDS70 is a unique system in which two protoplanets, PDS70b and c, have been discovered within the dust-depleted cavity of their disk, at $\sim$22 and 34au respectively, by direct imaging at infrared wavelengths. Subsequent detection of the planets in the H$\alpha$ line indicates that they are still accreting material through circumplanetary disks. In this Letter, we present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the dust continuum emission at 855$\mu$m at high angular resolution ($\sim$20mas, 2.3au) that aim to resolve the circumplanetary disks and constrain their dust masses. Our observations confirm the presence of a compact source of emission co-located with PDS70c, spatially separated from the circumstellar disk and less extended than $\sim$1.2au in radius, a value close to the expected truncation radius of the cicumplanetary disk at a third of the Hill radius. The emission around PDS70c has a peak intensity of $\sim$86$\pm$16 $\mu \mathrm{Jy}~\mathrm{beam}^{-1}$ which corresponds to a dust mass of $\sim$0.031M$_{\oplus}$ or $\sim$0.007M$_{\oplus}$, assuming that it is only constituted of 1 $\mu$m or 1 mm sized grains, respectively. We also detect extended, low surface brightness continuum emission within the cavity near PDS70b. We observe an optically thin inner disk within 18au of the star with an emission that could result from small micron-sized grains transported from the outer disk through the orbits of b and c. In addition, we find that the outer disk resolves into a narrow and bright ring with a faint inner shoulder., Comment: ApJ Letters, in press; 19 pages, 9 figures

Small-sized systems exhibit a finite number of routes to chaos. However, in extended systems, not all routes to complex spatiotemporal behavior have been fully explored. Starting from the sine-Gordon model of parametrically driven chain of damped nonlinear oscillators, we investigate a route to spatiotemporal chaos emerging from standing waves. The route from the stationary to the chaotic state proceeds through quasiperiodic dynamics. The standing wave undergoes the onset of oscillatory instability, which subsequently exhibits a different critical frequency, from which the complexity originates. A suitable amplitude equation, valid close to the parametric resonance, makes it possible to produce universe results. The respective phase-space structure and bifurcation diagrams are produced in a numerical form. We characterize the relevant dynamical regimes by means of the largest Lyapunov exponent, the power spectrum, and the evolution of the total intensity of the wave field., Comment: to be published in Chaos, Solitons & Fractals

We present a detailed analysis for a subset of the high resolution (~35 mas, or 5 au) ALMA observations from the Disk Substructures at High Angular Resolution Project (DSHARP) to search for faint 1.3 mm continuum emission associated with dusty circumplanetary material located within the narrow annuli of depleted emission (gaps) in circumstellar disks. This search used the Jennings et al. (2020) $\tt{frank}$ modeling methodology to mitigate contamination from the local disk emission, and then deployed a suite of injection-recovery experiments to statistically characterize point-like circumplanetary disks in residual images. While there are a few putative candidates in this sample, they have only marginal local signal-to-noise ratios and would require deeper measurements to confirm. Associating a 50% recovery fraction with an upper limit, we find these data are sensitive to circumplanetary disks with flux densities $\gtrsim 50-70$ $\mu$Jy in most cases. There are a few examples where those limits are inflated ($\gtrsim 110$ $\mu$Jy) due to lingering non-axisymmetric structures in their host circumstellar disks, most notably for a newly identified faint spiral in the HD 143006 disk. For standard assumptions, this analysis suggests that these data should be sensitive to circumplanetary disks with dust masses $\gtrsim 0.001-0.2$ M$_\oplus$. While those bounds are comparable to some theoretical expectations for young giant planets, we discuss how plausible system properties (e.g., relatively low host planet masses or the efficient radial drift of solids) could require much deeper observations to achieve robust detections., Comment: ApJ, in press; 30 pages, 15 figures (one is a 4-part figure set)

We present high resolution millimeter continuum and CO line observations for the circumbinary disk around V892 Tau to constrain the stellar and disk properties. The total mass of the two near-equal-mass A stars is estimated to be $6.0\pm0.2\,M_{\odot}$ based on our models of the Keplerian-dominated gas disk rotation. The detection of strong ionized gas emission associated with the two stars at 8 mm, when combined with previous astrometric measurements in the near-infrared, provides an updated view of the binary orbit with $a=7.1\pm0.1$ au, $e=0.27\pm0.1$, and $P=7.7\pm0.2$ yr, which is about half of a previously reported orbital period. The binary orbital plane is proposed to be near coplanar to the circumbinary disk plane (with a mutual inclination of only $\Delta=8\pm4.2$ deg; another solution with $\Delta=113$ deg is less likely given the short re-alignment timescale). An asymmetric dust disk ring peaking at a radius of 0.''2 is detected at 1.3 mm and its fainter counterparts are also detected at the longer 8 and 9.8 mm. The CO gas disk, though dominated by Keplerian rotation, presents a mild inner and outer disk misalignment, such that the inner disk to the SW and outer disk to the NE appear brighter than their counterparts at the opposite disk sides. The radial extension of the disk, its asymmetric dust ring, and the presence of a disk warp could all be explained by the interaction between the eccentric binary and the circumbinary disk, which we assume were formed with non-zero mutual inclination. Some tentatively detected gas spirals in the outer disk are likely produced by interactions with the low mass tertiary component located 4'' to the northeast. Our analyses demonstrate the promising usage of V892 Tau as an excellent benchmark system to study the details of binary--disk interactions., Comment: Accepted for publication in ApJ

Recent multi-wavelength ALMA observations of the protoplanetary disk orbiting around Elias 2-27 revealed a two armed spiral structure. The observed morphology together with the young age of the star and the disk-to-star mass ratio estimated from dust continuum emission make this system a perfect laboratory to investigate the role of self-gravity in the early phases of star formation. This is particularly interesting if we consider that gravitational instabilities could be a fundamental first step for the formation of planetesimals and planets. In this Letter, we model the rotation curve obtained by CO data of Elias 2-27 with a theoretical rotation curve including both the disk self-gravity and the star contribution to the gravitational potential. We compare this model with a purely Keplerian one and with a simple power-law function. We find that (especially for the $^{13}$CO isotopologue) the rotation curve is better described by considering not only the star, but also the disk self-gravity. We are thus able to obtain for the first time a dynamical estimate of the disk mass of $0.08\pm0.04\,M_{\odot}$ and the star mass of $0.46\pm0.03\,M_{\odot}$ (in the more general case), the latter being comparable with previous estimates. From these values, we derive that the disk is 17$\%$ of the star mass, meaning that it could be prone to gravitational instabilities. This result would strongly support the hypothesis that the two spiral arms are generated by gravitational instabilities., Comment: 14 pages, 9 figures, 2 tables, accepted for publication on ApJL

The "Disk Substructures at High Angular Resolution Project" (DSHARP) has revealed an abundance and ubiquity of rings and gaps over a large sample of young planet-forming disks, which are hypothesised to be induced by the presence of forming planets. In this context, we present the first attempt to directly image these young companions for 10 of the DSHARP disks, by using NaCo/VLT high contrast observations in L'-band instrument and angular differential imaging techniques. We report the detection of a point-like source candidate at 1.1" (174.9 au) for RU Lup, and at 0.42" (55 AU) for Elias 24. In the case of RU Lup, the proper motion of the candidate is consistent with a stationary background contaminant, based on the astrometry derived from our observations and available archival data. For Elias 24 the point-like source candidate is located in one of the disk gaps at 55 AU. Assuming it is a planetary companion, our analysis suggest a mass ranging from $0.5 M_J$ up to $5 M_J$, depending on the presence of a circumplanetary disk and its contribution to the luminosity of the system. However, no clear confirmation is obtained at this stage, and follow-up observations are mandatory to verify if the proposed source is physical, comoving with the stellar host, and associated with a young massive planet sculpting the gap observed at 55\,AU. For all the remaining systems, the lack of detections suggests the presence of planetary companions with masses lower than $5M_J$, based on our derived mass detection limits. This is consistent with predictions of both hydrodynamical simulations and kinematical signatures on the disk, and allows us to set upper limits on the presence of massive planets in these young disks., Comment: 19 pages, 10 figures; Published by AJ

While protoplanetary disks often appear to be compact and well-organized in millimeter continuum emission, CO spectral line observations are increasingly revealing complex behavior at large distances from the host star. We present deep ALMA maps of the $J=2-1$ transition of $^{12}$CO, $^{13}$CO, and C$^{18}$O, as well as the $J=3-2$ transition of DCO$^+$, toward the T Tauri star RU Lup at a resolution of $\sim0.3''$ ($\sim50$ au). The CO isotopologue emission traces four major components of the RU Lup system: a compact Keplerian disk with a radius of $\sim120$ au, a non-Keplerian ``envelope-like'' structure surrounding the disk and extending to $\sim260$ au from the star, at least five blueshifted spiral arms stretching up to 1000 au, and clumps outside the spiral arms located up to 1500 au in projection from RU Lup. We comment on potential explanations for RU Lup's peculiar gas morphology, including gravitational instability, accretion of material onto the disk, or perturbation by another star. RU Lup's extended non-Keplerian CO emission, elevated stellar accretion rate, and unusual photometric variability suggest that it could be a scaled-down Class II analog of the outbursting FU Ori systems., Comment: 32 pages, 13 figures, accepted by ApJ. Typos corrected. Data are available at https://zenodo.org/record/3928862

The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ~50 mas (~8 au), as well as HCO$^+$ $J=3-2$ observations at a resolution of ~100 mas. The dust continuum shows at least three rings atop faint, extended emission. Unresolved emission is detected at the center of the disk cavity at both wavelengths, likely due to a combination of dust and free-free emission. Compared to the 1.1 mm image, the 2.1 mm image shows a more pronounced "shoulder" near R~40 au, highlighting the utility of longer-wavelength observations for characterizing disk substructures. The spectral index $\alpha$ features strong radial variations, with minima near the emission peaks and maxima near the gaps. While low spectral indices have often been ascribed to grain growth and dust trapping, the optical depth of GM Aur's inner two emission rings renders their dust properties ambiguous. The gaps and outer disk ($R>100$ au) are optically thin at both wavelengths. Meanwhile, the HCO$^+$ emission indicates that the gas cavity is more compact than the dust cavity traced by the millimeter continuum, similar to other disks traditionally classified as "transitional.", Comment: 31 pages, 15 figures, accepted by ApJ. Self-calibration and imaging scripts, images, and self-calibrated visibilities available at https://zenodo.org/record/3628656#.XjEXyRNKiYV

While numerical simulations have been playing a key role in the studies of planet-disk interaction, testing numerical results against observations has been limited so far. With the two directly imaged protoplanets embedded in its circumstellar disk, PDS 70 offers an ideal testbed for planet-disk interaction studies. Using two-dimensional hydrodynamic simulations we show that the observed features can be well explained with the two planets in formation, providing strong evidence that previously proposed theories of planet-disk interaction are in action, including resonant migration, particle trapping, size segregation, and filtration. Our simulations suggest that the two planets are likely in 2:1 mean motion resonance and can remain dynamically stable over million-year timescales. The growth of the planets at $10^{-8}-10^{-7}~M_{\rm Jup}~{\rm yr}^{-1}$, rates comparable to the estimates from H$\alpha$ observations, does not destabilize the resonant configuration. Large grains are filtered at the gap edge and only small, (sub-)$\mu$m grains can flow to the circumplanetary disks and the inner circumstellar disk. With the sub-millimeter continuum ring observed outward of the two directly imaged planets, PDS 70 provides the first observational evidence of particle filtration by gap-opening planets. The observed sub-millimeter continuum emission at the vicinity of the planets can be reproduced when (sub-)$\mu$m grains survive over multiple circumplanetary disk gas viscous timescales and accumulate therein. One such possibility is if (sub-)$\mu$m grains grow in size and remain trapped in pressure bumps, similar to what we find happening in circumstellar disks. We discuss potential implications to planet formation in the solar system and mature extrasolar planetary systems., Comment: Accepted for publication in ApJL, 15 pages, 8 figures, An animated version of M_c = 2.5 M_Jup model in Figure 1 is temporarily available at http://jaehanbae.com/movies/pds70/

We present the discovery of a spatially unresolved source of sub-millimeter continuum emission ($\lambda=855$ $\mu$m) associated with a young planet, PDS 70 c, recently detected in H$\alpha$ emission around the 5 Myr old T Tauri star PDS 70. We interpret the emission as originating from a dusty circumplanetary disk with a dust mass between $2\times10^{-3}$ and $4.2 \times 10^{-3}$ Earth masses. Assuming a standard gas-to-dust ratio of 100, the ratio between the total mass of the circumplanetary disk and the mass of the central planet would be between $10^{-4}-10^{-5}$. Furthermore, we report the discovery of another compact continuum source located $0.074''\pm0.013''$ South-West of a second known planet in this system, PDS 70 b, that was previously detected in near-infrared images. We speculate that the latter source might trace dust orbiting in proximity of the planet, but more sensitive observations are required to unveil its nature., Comment: Accepted for publication on ApJL

We present new determinations of disk surface density, independent of an assumed dust opacity, for a sample of 7 bright, diverse protoplanetary disks using measurements of disk dust lines. We develop a robust method for determining the location of dust lines by modeling disk interferometric visibilities at multiple wavelengths. The disks in our sample have newly derived masses that are 9-27% of their host stellar mass, substantially larger than the minimum mass solar nebula. All are stable to gravitational collapse except for one which approaches the limit of Toomre-Q stability. Our mass estimates are 2-15 times larger than estimates from integrated optically thin dust emission. We derive depleted dust-to-gas ratios with typical values of ~$10^{-3}$ in the outer disk. Using coagulation models we derive dust surface density profiles that are consistent with millimeter dust observations. In these models, the disks formed with an initial dust mass that is a factor of ~10 greater than is presently observed. Of the three disks in our sample with resolved CO line emission, the masses of HD 163296, AS 209, and TW Hya are roughly 3, 115, and 40 times more massive than estimates from CO respectively. This range indicates that CO depletion is not uniform across different disks and that dust is a more robust tracer of total disk mass. Our method of determining surface density using dust lines is robust even if particles form as aggregates and is useful even in the presence of dust substructure caused by pressure traps. The low Toomre-Q values observed in this sample indicate that at least some disks do not accrete efficiently., Comment: 29 pages, 12 figures, submitted to ApJ, revised following comments from referee

We report new Atacama Large Millimeter/sub-millimeter Array (ALMA) Band 3 observations at 2.75 mm of the TD around SR 24S with an angular resolution of $\sim$0.11''$\times$ 0.09'' and a peak signal-to-noise ratio of $\sim24$. We detect an inner disk and a mostly symmetric ring-like structure that peaks at $\sim$0.32'', that is $\sim$37 au at a distance of $\sim$114.4 pc. The full width at half maximum of this ring is $\sim$28 au. We analyze the observed structures by fitting the dust continuum visibilities using different models for the intensity profile, and compare with previous ALMA observations of the same disk at 0.45 mm and 1.30 mm. We qualitatively compare the results of these fits with theoretical predictions of different scenarios for the formation of a cavity or large gap. The comparison of the dust continuum structure between different ALMA bands indicates that photoevaporation and dead zone can be excluded as leading mechanisms for the cavity formation in SR 24S disk, leaving the planet scenario (single or multiple planets) as the most plausible mechanism. We compared the 2.75 mm emission with published (sub-)centimeter data and find that the inner disk is likely tracing dust thermal emission. This implies that any companion in the system should allow dust to move inwards throughout the gap and replenish the inner disk. In the case of one single planet, this puts strong constraints on the mass of the potential planet inside the cavity and the disk viscosity of about $\lesssim$5 $M_{\rm{Jup}}$ and $\alpha\sim10^{-4}-10^{-3}$, respectively., Comment: Accepted to ApJ

ALMA surveys have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass in exoplanets, under the assumption that the mm disk continuum emission is optically thin. This optically thin assumption seems to be supported by recent DSHARP observations where the measured optical depths of spatially resolved disks are mostly less than one. However, we point out that dust scattering can considerably reduce the emission from an optically thick region. If that scattering is ignored, the optical depth will be considerably underestimated. An optically thick disk with scattering can be misidentified as an optically thin disk. Dust scattering in more inclined disks can reduce the intensity even further, making the disk look even fainter. The measured optical depth of $\sim$0.6 in several DSHARP disks can be naturally explained by optically thick dust with an albedo of $\sim$0.9 at 1.25 mm. Using the DSHARP opacity, this albedo corresponds to a dust population with the maximum grain size ($s_{max}$) of 0.1-1 mm. For optically thick scattering disks, the measured spectral index $\alpha$ can be either larger or smaller than 2 depending on if the dust albedo increases or decreases with wavelength. Using the DSHARP opacity, $\alpha<2$ corresponds to $s_{max}$ of 0.03-0.3 mm. We describe how this optically thick scattering scenario could explain the observed scaling between submm continuum sizes and luminosities, and might help ease the tension between the dust size constraints from polarization and dust continuum measurements. We suggest that a significant amount of disk mass can be hidden from ALMA observations at short millimeter wavelengths. For compact disks smaller than 30 au, we can easily underestimate the dust mass by more than a factor of 10. Longer wavelength observations (e.g. VLA or SKA) are desired to probe the dust mass in disks., Comment: 15 pages, 11 figures, Accepted for publication in ApJL

We present an analysis of ALMA 1.25 millimeter continuum observations of spiral structures in three protoplanetary disks from the Disk Substructures at High Angular Resolution Project. The disks around Elias 27, IM Lup, and WaOph 6 were observed at a resolution of $\sim40-60$ mas ($\sim6-7$ au). All three disks feature $m=2$ spiral patterns in conjunction with annular substructures. Gas kinematics established by $^{12}$CO $J=2-1$ observations indicate that the continuum spiral arms are trailing. The arm-interarm intensity contrasts are modest, typically less than 3. The Elias 27 spiral pattern extends throughout much of the disk, and the arms intersect the gap at $R\sim69$ au. The spiral pattern in the IM Lup disk is particularly complex-it extends about halfway radially through the disk, exhibiting pitch angle variations with radius and interarm features that may be part of ring substructures or spiral arm branches. Spiral arms also extend most of the way through the WaOph 6 disk, but the source overall is much more compact than the other two disks. We discuss possible origins for the spiral structures, including gravitational instability and density waves induced by a stellar or planetary companion. Unlike the millimeter continuum counterparts of many of the disks with spiral arms detected in scattered light, these three sources do not feature high-contrast crescent-like asymmetries or large ($R>20$ au) emission cavities. This difference may point to multiple spiral formation mechanisms operating in disks., Comment: 19 pages, 9 figures, minor typos corrected

ALMA observations of protoplanetary disks acquired by the Disk Substructure at High Angular Resolution Project (DSHARP) resolve the dust and gas emission on angular scales as small as 3 astronomical units, offering an unprecedented detailed view of the environment where planets form. In this article, we present and discuss observations of the HD 163296 protoplanetary disk that imaged the 1.25 mm dust continuum and $^{12}$CO J=2-1 rotational line emission at a spatial resolution of 4 and 10 au, respectively. The continuum observations resolve and allow us to characterize the previously discovered dust rings at radii of 67 and 100 au. They also reveal new small scale structures, such as a dark gap at 10 au, a bright ring at 15 au, a dust crescent at a radius of 55 au, and several fainter azimuthal asymmetries. The observations of the CO and dust emission inform about the vertical structure of the disk and allow us to directly constrain the dust extinction optical depth at the dust rings. Furthermore, the observed asymmetries in the dust continuum emission corroborate to the hypothesis that the complex structure of the HD 163296 disk is the result of the gravitational interaction with yet unseen planets., Comment: ApJ Letters, in press. 19 pages, 11 figures

A large fraction of the protoplanetary disks observed with ALMA display multiple well-defined and nearly perfectly circular rings in the continuum, in many cases with substantial peak-to-valley contrast. The DSHARP campaign shows that several of these rings are very narrow in radial extent. In this paper we test the hypothesis that these dust rings are caused by dust trapping in radial pressure bumps, and if confirmed, put constraints on the physics of the dust trapping mechanism. We model this process analytically in 1D, assuming axisymmetry. By comparing this model to the data, we find that all rings are consistent with dust trapping. Based on a plausible model of the dust temperature we find that several rings are narrower than the pressure scale height, providing strong evidence for dust trapping. The rings have peak absorption optical depth in the range between 0.2 and 0.5. The dust masses stored in each of these rings is of the order of tens of Earth masses, though much ambiguity remains due to the uncertainty of the dust opacities. The dust rings are dense enough to potentially trigger the streaming instability, but our analysis cannot give proof of this mechanism actually operating. Our results show, however, that the combination of very low alpha_turb << 5e-4 and very large grains a_grain >> 0.1 cm can be excluded by the data for all the rings studied in this paper., Comment: Accepted for publication in ApJ Letters

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

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

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

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

We present a detailed analysis of the high-angular resolution (0.''037, corresponding to 5 au) observations of the 1.25 mm continuum and $^{12}$CO $2-1$ emission from the disk around the T Tauri star AS 209. AS 209 hosts one of the most unusual disks from the DSHARP sample, the first high angular resolution ALMA survey of disks (Andrews et al. 2018), as nearly all of the emission can be explained with concentric Gaussian rings. In particular, the dust emission consists of a series of narrow and closely spaced rings in the inner $\sim$60 au, two well-separated bright rings in the outer disk, centered at 74 and 120 au, and at least two fainter emission features at 90 and 130 au. We model the visibilities with a parametric representation of the radial surface brightness profile, consisting of a central core and 7 concentric Gaussian rings. Recent hydro-dynamical simulations of low viscosity disks show that super-Earth planets can produce the multiple gaps seen in AS 209 millimeter continuum emission. The $^{12}$CO line emission is centrally peaked and extends out to $\sim$300 au, much farther than the millimeter dust emission. We find axisymmetric, localized deficits of CO emission around four distinct radii, near 45, 75, 120 and 210 au. The outermost gap is located well beyond the edge of the millimeter dust emission, and therefore cannot be due to dust opacity and must be caused by a genuine CO surface density reduction, due either to chemical effects or depletion of the overall gas content., Comment: ApJ Letters, in press

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

Emission from protostars at centimeter radio wavelengths has been shown to trace the free-free emission arising from ionizing shocks as a result of jets and outflows driven by protostars. Therefore, measuring properties of protostars at radio frequencies can provide valuable insights into the nature of their outflows and jets. We present a C-band (4.1 cm and 6.4 cm) survey of all known protostars (Class 0 and Class I) in Perseus as part of the VLA Nascent Disk and Multiplicity (VANDAM) Survey. We examine the known correlations between radio flux density and protostellar parameters such as bolometric luminosity and outflow force, for our sample. We also investigate the relationship between radio flux density and far-infrared line luminosities from Herschel. We show that free-free emission originates most likely from J-type shocks; however, the large scatter indicates that those two types of emission probe different time and spatial scales. Using C-band fluxes, we removed an estimation of free-free contamination from the corresponding Ka-band (9 mm) flux densities that primarily probe dust emission from embedded disks. We find that the compact ($<$~1") dust emission is lower for Class I sources (median dust mass 96 M$_{\oplus}$) relative to Class 0 (248 M$_{\oplus}$), but several times higher than in Class II (5-15 M$_{\oplus}$). If this compact dust emission is tracing primarily the embedded disk, as is likely for many sources, this result provides evidence for decreasing disk masses with protostellar evolution, with sufficient mass for forming giant planet cores primarily at early times., Comment: 90 pages, 21 figures, 10 tables; accepted for publication in ApJS