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 observations of the (J = 2-1) and (1-0) lines of 13 CO and C 18 O and the (J = 1-0) line of C 17 O for five protoplanetary disks. We retrieve CO gas density distributions, using three independent methods: (1) a thermochemical modeling framework based on the CO data, the broadband spectral energy distribution, and the millimeter 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 10-100 times lower than the interstellar medium ratio. The CO gas distributions between 150 and 400 au match well with models of viscous disks, supporting the longstanding 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.

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 0.2. C 2 H and HCN have lower opacities and signal-to-noise ratios, making surface fitting more difficult, and could only be reliably constrained in AS209, HD163296, and MWC480, with  z r 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 thermochemical 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.

The European FP7 project DIANA has performed a coherent analysis of a large set of observational data of protoplanetary disks by means of thermo-chemical disk models. The collected data include extinction-corrected stellar UV and X-ray input spectra (as seen by the disk), photometric fluxes, low and high resolution spectra, interferometric data, emission line fluxes, line velocity profiles and line maps, which probe the dust, polycyclic aromatic hydrocarbons (PAHs) and the gas in these objects. We define and apply a standardized modeling procedure to fit these data by state-of-the-art modeling codes (ProDiMo, MCFOST, MCMax), solving continuum and line radiative transfer (RT), disk chemistry, and the heating and cooling balance for both the gas and the dust. 3D diagnostic RT tools (e.g., FLiTs) are eventually used to predict all available observations from the same disk model, the DIANA-standard model. Our aim is to determine the physical parameters of the disks, such as total gas and dust masses, the dust properties, the disk shape, and the chemical structure in these disks. We allow for up to two radial disk zones to obtain our best-fitting models that have about 20 free parameters. This approach is novel and unique in its completeness and level of consistency. It allows us to break some of the degeneracies arising from pure Spectral Energy Distribution (SED) modeling. In this paper, we present the results from pure SED fitting for 27 objects and from the all inclusive DIANA-standard models for 14 objects. Our analysis shows a number of Herbig Ae and T Tauri stars with very cold and massive outer disks which are situated at least partly in the shadow of a tall and gas-rich inner disk. The disk masses derived are often in excess to previously published values, since these disks are partially optically thick even at millimeter wavelength and so cold that they emit less than in the Rayleigh–Jeans limit. We fit most infrared to millimeter emission line fluxes within a factor better than 3, simultaneously with SED, PAH features and radial brightness profiles extracted from images at various wavelengths. However, some line fluxes may deviate by a larger factor, and sometimes we find puzzling data which the models cannot reproduce. Some of these issues are probably caused by foreground cloud absorption or object variability. Our data collection, the fitted physical disk parameters as well as the full model output are available to the community through an online database (http://www.univie.ac.at/diana).

The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a detailed, high-resolution (∼10–20 au) view of molecular line emission in five protoplanetary disks at spatial scales relevant for planet formation. Here we present a systematic analysis of chemical substructures in 18 molecular lines toward the MAPS sources: IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. We identify more than 200 chemical substructures, which are found at nearly all radii where line emission is detected. A wide diversity of radial morphologies—including rings, gaps, and plateaus—is observed both within each disk and across the MAPS sample. This diversity in line emission profiles is also present in the innermost 50 au. Overall, this suggests that planets form in varied chemical environments both across disks and at different radii within the same disk. Interior to 150 au, the majority of chemical substructures across the MAPS disks are spatially coincident with substructures in the millimeter continuum, indicative of physical and chemical links between the disk midplane and warm, elevated molecular emission layers. Some chemical substructures in the inner disk and most chemical substructures exterior to 150 au cannot be directly linked to dust substructure, however, which indicates that there are also other causes of chemical substructures, such as snowlines, gradients in UV photon fluxes, ionization, and radially varying elemental ratios. This implies that chemical substructures could be developed into powerful probes of different disk characteristics, in addition to influencing the environments within which planets assemble. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.

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 ∼0 15 and with a channel spacing of 200 m s −1 , we find perturbations from Keplerian rotation in the projected velocity fields of both disks (5% of the local Keplerian velocity), suggestive of large-scale (tens of astronomical units 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 f , v r , and v z . Using the optically thick 12 CO emission as a probe of the gas temperature, local variations of ≈3 K (≈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, ∼1 M Jup , at a radius of ≈245 au. In the disk of HD 163296, the kinematic substructures were consistent with previous studies of Pinte et al. and Teague et al. advocating for multiple ∼1 M 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.