Molecules with ALMA at Planet-forming Scales (MAPS) VI: Distribution of the small organics HCN, C2H, and H2CO

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.


INTRODUCTION
Young stars have flattened disks of gas and dust rotating around them, a natural outcome of the star formation process. The material available in these protoplanetary disks will eventually be accreted onto the star, blown out by winds, or be incorporated into planets (see reviews by Williams & Cieza 2011;Andrews 2020;Öberg & Bergin 2021). The chemical composition of nascent planets is therefore directly linked to the chemical reservoir of disks. Small organics are some of the main carriers of C, N, and O in protoplanetary disks. However, the ability of planets to access these organic reservoirs depends on the distribution of this material, particularly in the inner few tens of au where rocky planets are expected to form.
Although the main carrier of N in disks is expected to be N 2 , this molecule has no observable lines in cold gas. HCN is a secondary N carrier, but its rotational lines in the millimeter are readily detected in disks (Öberg et al. 2011;Chapillon et al. 2012;Bergner et al. 2019). In addition, observations of mid-infrared HCN lines have shown that it is abundant in the inner ∼ 10 au of disks (Salyk et al. 2011;Najita et al. 2013Najita et al. , 2018. HCN is therefore one of our best tracers of the N budget in disks and allows us to characterize the spatial distribution of N. Moreover, HCN is thought to be a crucial molecule to the emergence of life, as it is the starting point in the chemistry leading to proteins and RNA (Powner et al. 2009;Patel et al. 2015;Becker et al. 2019).
In addition to revealing the distribution of the organic and elemental reservoir, molecular lines provide us key information on the physical and chemical conditions prevailing in disks. The distribution of HCN and C 2 H has been shown to depend on the UV radiation field and on the abundances of C and O (Du et al. 2015;Bergin et al. 2016;Cleeves et al. 2018). The distribution of these small organics are also potentially sensitive to density and temperature variations. For example, HCN has been proposed to be a useful tool to study the formation of giant planets in disks, as they will heat up the * NASA Hubble Fellowship Program Sagan Fellow † NASA Hubble Fellow surrounding gas and result in enhanced HCN emission around the forming planet . H 2 CO is a third commonly observed small organic in disks that is connected to the oxygen cycle. It is an important precursor of O-bearing Complex Organic Molecules (COMs), which as methanol, acetaldehyde, ethylene glycol, and formamide (Watanabe & Kouchi 2002;Chuang et al. 2017). It is hypothesised that Obearing COMs are preferentially locked in ices and removed from the gas phase as the dust grows and settles to the midplane (e.g, Bergin et al. 2010;Hogerheijde et al. 2011;Du et al. 2017). This would explain why O-bearing COMs have been more difficult to detect in disks compared to nitriles or hydrocarbons (e.g., Chapillon et al. 2012;Walsh et al. 2016;Bergner et al. 2018;Favre et al. 2018;Loomis et al. 2018;Carney et al. 2019). H 2 CO lines, however, are much brighter in disks than lines from more complex Obearing species because H 2 CO has a simpler molecular structure and therefore the internal energy is distributed in fewer molecular lines (Aikawa et al. 2003;Pegues et al. 2020). Observations of H 2 CO therefore provide a view of the spatial distribution of the O-bearing molecules in disks.
While rich substructure has been observed in dust continuum emission of numerous disks (e.g., Andrews et al. 2018;Long et al. 2018;Pérez et al. 2020), the chemical structure of these disks has been less explored at high angular resolution. The combination of sensitivity and spatial resolution needed to spatially resolve molecular line emission on planet-forming scales has been the main limitation so far in characterizing the distribution of most key organic molecules. However, sensitive observations at moderate angular resolution have demonstrated the presence of a rich substructure in the emission of various molecular lines (e.g., Henning & Semenov 2008;Mathews et al. 2013;Bergin et al. 2016;Huang et al. 2017;Öberg et al. 2017;Salinas et al. 2017;Cazzoletti et al. 2018;Kastner et al. 2018;Miotello et al. 2019;Pegues et al. 2020;Teague & Loomis 2020).
In this paper, we present observations from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. Here, we focus on molecular line emission from the small organic molecules HCN, C 2 H, and H 2 CO in five nearby and otherwise wellstudied protoplanetary disks. HCN, C 2 H, and H 2 CO have been previously observed in these disks, but at lower angular resolution (∼0. 5). In particular, Bergner et al. (2019) presented observations of HCN and C 2 H toward four of these disks, and Pegues et al. (2020) presented observations of H 2 CO toward the five disks. The MAPS data presented in this paper provide a combination of high angular resolution with excellent brightness sensitivity, which enables us to identify substructure in the line emission of these small organics. The observations are briefly described in Section 2. In Section 3, we present the spatial distribution of the emission. In Section 4, we retrieve the column density and excitation temperature profiles. We discuss the results in Section 5 and summarize our conclusions in Section 6.

OBSERVATIONS
The observations are part of the MAPS ALMA Large Program (2018.1.01055.L), where a large number of molecular lines were observed at high spatial resolution to search for links between chemistry and planet formation in protoplanetary disks. SeeÖberg & MAPS team (2021) for a detailed description of the scope and aims of MAPS. The MAPS sample includes disks around both T Tauri (IM Lup, GM Aur, and AS 209) and Herbig Ae (HD 163296 and MWC 480) stars, with a wide variety of dust substructures (e.g., rings, gaps, and spirals). A detailed description of the sources can be found inÖberg & MAPS team (2021), and a brief summary is given in Table 1.
Four spectral settings, two in Band 6 and two in Band 3, were observed between 2019 and 2020. The Band 6 setting covered the HCN J = 3 − 2, C 2 H N = 3 − 2, and the H 2 CO 3 03 − 2 02 lines. The Band 3 setting covered the HCN J = 1 − 0 lines and the C 2 H N = 1 − 0 ladder. See Table 2 for a full list of QNs or transitions observed.
We subsequently just refer to these lines as 3 − 2 and 1 − 0 for simplicity. The data include observations using two different array configurations, to provide sufficient dynamic range in baseline length to recover both largescale and small-scale emission.
The data were first calibrated by the ALMA staff using standard routines, and then self-calibrated by the MAPS team to improve the signal-to-noise ratio of the data. A detailed description of the observations and calibration process can be found inÖberg & MAPS team (2021).
The imaging of the observed visibilities was done with the tclean task of the CASA 6 software. A Keplerian mask was created during the CLEANing process, by selecting regions in each channel where line emission is expected given the Keplerian rotation of each disk. Novel techniques were applied to improve the image quality, in particular the scaling of the CLEAN residual map, which can affect substantially the data quality when combining data from different configuration arrays. We encourage the reader to read a more detailed description of the imaging strategy in Czekala & MAPS team (2021).
We make use of the fiducial MAPS images, which have a circular beam of 0. 15 (corresponding to linear spatial scales of 15 − 24 au) for the 3 − 2 lines, and a circular beam of 0. 3 (corresponding to linear spatial scales of 30 − 48 au) for the 1 − 0 lines. In addition to these fiducial images, we used a second set of images for the 3 − 2 lines, that were produced by applying a taper to match the resolution of the 1 − 0 lines (see Section 6.2 in Czekala & MAPS team 2021 for a detailed description of how these images were produced). This set of maps are available on the MAPS project homepage (https://www.alma-maps.info). The H 2 CO line is weaker per beam than that of the other small organics, resulting in images with low S/N at high angular (0. 15) resolution. GM Aur is an exception here, as it presents bright H 2 CO emission that is well-detected in the 0. 15 resolution image. For consistency, here we use the tapered version of the H 2 CO images (with 0. 3 resolution) for all sources, which are more sensitive to low surface brightness and extended features. Disk integrated line fluxes for each line can be found in Table A.1.

Zeroth-moment maps
The velocity-integrated emission, or "zeroth moment", maps of the HCN, C 2 H, and H 2 CO 3 − 2 and 1 − 0 lines are shown in Fig. 1, the dust continuum emission at 260 GHz (with ∼0. 1 angular resolution) is shown in the upper row for comparison. A description of the The spectroscopic constants are taken from the CDMS database (Müller et al. 2001(Müller et al. , 2005 for HCN, and from the JPL database (Pickett et al. 1998) for C2H and H2CO. Measurements are provided by Ahrens et al. (2002) for HCN, by Sastry et al. (1981) and Padovani et al. (2009) for C2H, and by Bocquet et al. (1996) for H2CO. a Velocity shift from the brightest hyperfine component in the group.
b Relative intensity with respect to the brightest hyperfine component in the group.
continuum emission images can be found in Sierra & MAPS team (2021). For HCN J = 3 − 2 and J = 1 − 0, all hyperfine lines were included to generate the zerothmoment maps, while for C 2 H only the two brightest (and blended) lines are included in the N = 3 − 2 line, and the first and brightest component in the N = 1 − 0 line. For more details on the zeroth-moment map generation process, see Law & MAPS team (2021a). The 3−2 lines have higher fluxes than the lower-energy 1 − 0 lines for all disks. HCN J = 3 − 2 is strongly detected in all disks, and presents a variety of substructures, including central holes, single, and multiple rings. HCN J = 1 − 0 presents similar spatial distributions to the J = 3−2 line in AS 209, HD 163296, and MWC 480. Several rings are seen toward HD 163296, including a faint ring in the outer disk (centered at ∼ 400 au), far beyond the dust continuum edge. Toward IM Lup, the HCN J = 1−0 line emission is fainter and more extended than in the other disks, resulting in a lower-S/N zerothmoment map. Toward GM Aur, the HCN J = 1 − 0 line is mostly detected in the central regions of the disk, although an outer diffuse ring is seen in the radial profile (see also Law & MAPS team 2021a).
The C 2 H N = 3 − 2 and N = 1 − 0 lines have spatial distributions and radial extents similar to HCN in AS 209, HD 163296, and MWC 480. Toward IM Lup and GM Aur, the N = 3 − 2 line is fainter than in the other disks, resulting in lower-S/N zeroth-moment maps, and the N = 1 − 0 line is mostly detected in the central regions of the disks, similar to HCN.
Fewer substructures are observed for H 2 CO compared to HCN and C 2 H, although a few ring-like structures are seen, in particular in HD 163296. The H 2 CO line emission is more extended than HCN and C 2 H, and presents a central hole of varying depths in all the disks except GM Aur and MWC 480, where only a tentative dip is observed. GM Aur presents exceptionally bright emission compared to the other disks.
Normally, for molecules other than 12 CO and 13 CO, it is not possible to derive the height of the emission directly from the spatially resolved image cubes (e.g., Rosenfeld et al. 2013;Pinte et al. 2018a). This is be- cause other lines are typically fainter and arise in less elevated disk layers. However, thanks to the higher in- clination and relative proximity of the HD 163296 disk, and the high angular resolution and sensitivity of the MAPS data, we can distinguish emission arising in a flat disk from emission arising in the elevated layers above the disk midplane. Indeed, the HCN and C 2 H 3 − 2 lines observed toward HD 163296 show a clear asymmetry along the disk minor axis that cannot be explained with emission from a flat disk. Part of the asymmetry, also seen as projected ellipses that look shifted to the south with respect to the disk center, is due to the fact that we are observing the emission from the near cone that is facing us. This effect is clearly seen in the zerothmoment map, and corresponds to the heights estimated in Law & MAPS team (2021b), i.e., z/r 0.1. See Law & MAPS team (2021b) for a detailed description of the emission surfaces for HCN and C 2 H in the MAPS disk sample.

Radial profiles
To better characterize the spatial distribution of the emission we used the radial profile function of the Python package GoFish (Teague 2019) to produce deprojected radial intensity profiles from the zerothmoment maps, considering the disk inclination and po-sition angle (see Table 1). Although it is possible to take into account the vertical height of the emission to derive the radial intensity profiles (Law & MAPS team 2021a), here we assume a flat disk but produce the deprojected profiles by averaging the emission within a 30 • wedge along the disk major axis. The C 2 H lines were not detected at sufficiently high S/N toward IM Lup and GM Aur, so we averaged the emission azimuthally to produce the deprojected radial profiles. For additional details on the deprojected radial profiles, see Law & MAPS team (2021a). The resulting profiles for the 3 − 2 (solid) and 1 − 0 (dashed) lines are shown in Fig. 2. For comparison, the profiles of the 260 GHz and 90 GHz dust continuum emission are shown in the first row.
In general, the HCN and C 2 H emission is more compact than that of H 2 CO -their radial extents are smaller by a factor ∼ 2. Moreover, their radial extents are similar to that of the dust continuum emission, although faint HCN and C 2 H emission is observed up to ∼ 400 au in HD 163296. IM Lup is an exception, as it presents similar radial extents in all molecular lines. One possibility to explain the lack of HCN and C 2 H emission in the outermost disk regions where H 2 CO is observed is that the H 2 gas density is not high enough to excite the lines. Indeed, for gas temperatures between 10 and 50 K, the HCN 3 − 2 line has a critical density of (6 − 10) × 10 6 cm −3 , while the H 2 CO 3 − 2 line has a critical density that is 10 times lower (Shirley 2015). However, models of these disks find that gas densities are > 10 7 cm −3 in the outer disk midplanes, enough to excite these lines (e.g., Zhang & MAPS team 2021). This suggests that the abundance of HCN and C 2 H decreases in the outermost disk regions where H 2 CO is still abundant. We note that no additional emission arises in the outermost disk regions in the lower-resolution (0. 3) images that have better S/N.
The radial profiles are consistent with profiles previously observed at lower angular resolution (0. 5) for HCN and C 2 H in IM Lup, AS 209, HD 163296, and MWC 480 (Bergner et al. 2019), and for H 2 CO in all sources (Pegues et al. 2020). However, our higher angular resolution observations reveal additional radial substructures, in particular for HCN and C 2 H. We focus on the 3 − 2 lines which show more structure than the 1 − 0 lines, due to the higher signal-to-noise ratios and angular resolution (see Law & MAPS team 2021a for additional substructure analysis).

HCN
HCN shows a central hole in four of five disks, the exception being MWC 480 where the emission is centrally peaked. A main ring of emission is seen for IM Lup, GM Aur, and AS 209, although the rings are not Gaussian but present one or more emission shoulders that could result from blended ring components not resolved in the MAPS data. IM Lup shows an unusually broad, almost plateau-like ring. Two bright rings are seen for HD 163296, centered at ∼ 30 and 110 au, and a fainter plateau that extends up to ∼ 400 au. In addition, a shallow gap is seen at ∼ 327 au along with a faint ring centered at ∼ 366 au (Law & MAPS team 2021a). The ring is better seen in the zeroth-moment map (see Fig.1). GM Aur and AS 209 also present faint HCN emission in the outer disk. In addition to the central emission component, MWC 480 shows an emission shoulder at ∼ 30 au, and a partially blended ring at ∼ 90 au (see also Law & MAPS team 2021a). HCN radial profiles observed at similar linear spatial resolution toward other disks around T Tauri stars generally show centrally peaked emission profiles with emission bumps in the outer disk (e.g., Hily-Blant et al. 2019;Kastner et al. 2018).

C2H
C 2 H shows a central hole for all disks, which is generally larger than the hole seen for HCN. Several unresolved ring components are seen toward IM Lup and GM Aur, similar to what is observed in HCN. AS 209 presents a single ring of C 2 H emission, very similar in size to the ring seen in HCN. Two rings are seen toward HD 163296, centered roughly at the same locations of the HCN rings. The two faint HCN rings seen in the outer disk are also observed for C 2 H in the zerothmoment map (see Fig. 1), but are harder to see in the radial profiles. The presence of these faint C 2 H outer rings is confirmed with the tapered (0. 3 resolution) image, which has higher S/N (Law & MAPS team 2021a). Toward MWC 480, a single C 2 H ring is observed centered at the same location of the outermost HCN ring. Interestingly, the innermost C 2 H ring seen in HD 163296 and the rings seen in AS 209 and MWC 480 are all located near a gap in the dust continuum emission. Similar ring-like structures with large central cavities have been observed for C 2 H toward other disks, like TW Hya ) and V4046 Sgr (Kastner et al. 2018).

H2CO
H 2 CO radial profiles present more complicated substructure, such as broad ring-like emission components with several shoulders and/or local maxima. In the inner disk, IM Lup, AS 209, and HD 163296 present a central depression in H 2 CO emission, while GM Aur and MWC 480 show centrally peaked emission profiles. In the outer disks (> 100 au), the slope at which the H 2 CO intensity decreases is much more shallow compared to HCN and C 2 H, in particular for the disks around HD 163296 and MWC 480.
The radial profiles show spatial links between H 2 CO substructures, namely line emission peaks, and the outer edge of the millimeter continuum disk. H 2 CO peaks are aligned with the outer edge of the continuum in IM Lup and HD 163296. Additionally, the higherresolution (0. 15) images of AS 209, as noted in Law & MAPS team (2021a), also show the presence of an additional H 2 CO ring that is radially coincident with the dust edge. Toward HD 163296, a spatial link between a H 2 CO line emission peak and outer dust edge was seen in lower-resolution data from Carney et al. (2017). Spatial links between millimeter dust edges and H 2 CO line emission have also been observed toward other disks, like TW Hya (Öberg et al. 2017), and V4046 Sgr (Kastner et al. 2018), as well as in four out of a sample of 15 disks in the survey presented by Pegues et al. (2020). All MAPS disks were included in those surveyed by Pegues et al. (2020), but a close association between an outer H 2 CO ring and millimeter dust edge in the MAPS disks, as reported here, was not clearly seen in Pegues et al. (2020), due to these outer H 2 CO rings not being clearly spatially resolved. Thus, when considering these previous survey results with the new MAPS observations, this suggests that at least ∼50% of disks may show some spatial links between H 2 CO substructures and the outer edge of the continuum disk.
In summary, we find some similarities between the HCN and C 2 H radial profiles. Furthermore, several of the H 2 CO peaks spatially coincide with the HCN and C 2 H peaks, for example the first H 2 CO peak seen toward AS 209. However, most of the H 2 CO features do not have direct counterparts in HCN or C 2 H, and the H 2 CO radial extents are larger.

COLUMN DENSITIES
To derive column density profiles, we first produce deprojected spectra at different radii using the radial spectra function of GoFish. The radial bins are set to 1/4 of the beam major axis, and the disk is assumed to be vertically flat. Before averaging the spectra in each radial bin, the velocity is shifted according to the Keplerian rotation of the disk, so that all the spectra are aligned. The stellar masses needed to compute the Keplerian velocities are taken fromÖberg & MAPS team (2021) and are estimated from gas dynamics in Teague & MAPS team (2021). The spectra are azimuthally averaged in those disks with low S/N (C 2 H for IM Lup and GM Aur, and H 2 CO for all sources), and averaged within a 30 • wedge along the disk major axis for the others. The main motivation for using a 30 • wedge when possible is to avoid the uncertainty of the emission surface height, which strongly influences the spatial location of emission near the minor axis.

Column density retrieval
To derive the column densities, we assume that the emission is in local thermal equilibrium (LTE) and that the emission fills the beam. The LTE assumption is reasonable, considering the high densities present in protoplanetary disks. Cataldi & MAPS team (2021) performed a non-LTE fit to the HCN lines and found similar column densities, suggesting our LTE assumption is a good approximation, at least for HCN. The line intensity is computed as where B ν is the Planck function, T ex is the excitation temperature, T bg is the background temperature, and τ ν is the line optical depth. The background temperature T bg is fixed to the maximum of the continuum brightness temperature or the cosmic microwave background (CMB) temperature (2.73 K).
At the line center, the optical depth is given by where A ul is the Einstein coefficient, ν is the frequency of the line, σ is the Gaussian line width and N u is the upper-level column density, Here, N tot corresponds to the total column density of the molecule, g u is the upper-level degeneracy, E u is the upper-level energy, and Q(T ex ) is the partition function. The spectroscopic constants and partition functions were taken from the CDMS and JPL catalogs and are listed in Table 2. The free parameters in the fit are the total column density N tot , the FWHM ∆V = 2 √ 2 ln 2σ, and the excitation temperature in the case of HCN and C 2 H, as their lines present hyperfine structure. For H 2 CO, the excitation temperature is fixed because only a single line is available. We also allow for a small offset from the source velocity v lsr , constrained between −0.25 and 0.25 km s −1 . We explore the parameter space using the emcee package (Foreman-Mackey et al. 2013). The best fit is taken as the 50th percentile of the samples in the marginalized distributions, and the associated errors are computed from the 16th and 84th percentiles.

HCN and C 2 H fitting method
We follow a method similar to the one implemented in Bergner et al. (2019). The HCN and C 2 H lines present hyperfine structure that is spectrally resolved. We fit the observed spectra assuming that all the hyperfine lines have the same excitation temperature. Each hyperfine line is assumed to be a Gaussian, such that the total line opacity is given by where τ i,0 is the line center opacity of the ith component, v i is the relative velocity of each component with respect to the brightest line, and v lsr is the source velocity. The relative intensities r i as well as the relative velocities v i are listed in Table 2. The width of spectral lines is dominated by thermal broadening, which depends on the gas temperature: where m is the molecule's mass. Under the assumption of LTE, T gas = T ex . However, we find that the observed line widths for HCN and C 2 H are broader than their thermal widths would imply. This is most likely produced by beam smearing, where line emission from different radii in the disk that have different velocities are combined in a single beam. In addition to beam smearing, it is possible that emission from both the back and front sides of the flared disk, which have slightly different projected velocities, are being combined in the observed spectra. The inner disk is particularly affected by these issues, which makes it difficult to fit the lines without artificially decreasing the line opacity or increasing the excitation temperature. A way to alleviate this problem is to first compute the line intensity assuming a purely thermal line width, then compute the total flux, and finally redistribute this flux in the line but now for the observed and broader line width. With this method, the fit to the lines is better, and the inferred line opacity-and hence the column density-is more realistic. The same solution was applied to fit the CN and DCN lines in the MAPS disks (Bergner & MAPS team 2021;Cataldi & MAPS team 2021). The fit to the observed spectra was done in two steps. First, we fit the 1 − 0 and 3 − 2 lines simultaneously to obtain a better constraint on the excitation temperature. The tapered (0. 3) versions of the 3 − 2 line cubes were used for this step, to match the resolution of the 1 − 0 line cubes. We included priors on the excitation temperature informed by the gas temperature traced by the optically thick 12 CO 2-1 line. Law & MAPS team (2021b) found that the HCN and C 2 H emission arises from relatively flat surfaces z/r 0.1 in AS 209, HD 163296, and MWC 480, while the CO isotopologues arise in more elevated disk layers. We thus restrict the rotational temperature between 10 K and the gas temperature inferred from the peak brightness temperature of the optically thick 12 CO 2-1 line, T gas ( 12 CO) (Law & MAPS team 2021b), for both HCN and C 2 H. The results for the simultaneous fit to the 1 − 0 and 3 − 2 lines are shown in Appendix B. The resulting excitation temperature and column density profiles are shown in Fig. B.1, and the gray curves show the T gas ( 12 CO) profiles for each disk. HCN presents rotational temperatures between 15 and 40 K, while C 2 H presents higher values of 15 − 60 K. The column densities are on the order of 10 13 − 10 15 cm −2 for both molecules. Despite the lower angular resolution of 0. 3 of the data used in these analyses, two column density rings are resolved in both HCN and C 2 H for HD 163296, and a single ring or central component for the other sources.
In the second step, and to obtain a better-resolved column density profile, we fit the high angular resolution 3 − 2 lines only. We used the excitation temperature found in the previous step (T ex,0 ) as an initial parameter, and constrain T ex,0 − 10 K < T ex < T ex,0 + 10 K for both HCN and C 2 H. Figure 3 shows a selection of observed spectra and fit for each species and disk. Note that the radial extent is different for each disk. The lines are fit well in most of the disk, but poorly fit in the inner disk, where the lines are substantially broader, resulting in the individual hyperfine components being completely blended. For this reason, we used the 12 CO or 13 CO temperature as the excitation temperature in regions where the fit does not work properly. For C 2 H, the T gas ( 12 CO) was used for all sources. For HCN, the T gas ( 12 CO) was used for IM Lup, GM Aur, and AS 209. For the two disks around Herbig Ae stars, however, the T gas ( 13 CO) was used instead, since the T gas ( 12 CO) produced a large discontinuity in the column density profile.
An independent retrieval of the HCN column density and excitation temperature is presented in Bergner & MAPS team (2021) and Cataldi & MAPS team (2021), using a method similar to the one presented here but with different constraints on the excitation temperature. A comparison between the HCN column densities presented in these three papers is shown in Appendix C. The resulting column density profiles and excitation temperatures are consistent with the ones presented here (see Fig. C.1).

H 2 CO fitting method
For H 2 CO, we fit a single Gaussian to the J = 3 − 2 line in each radial bin. In order to constrain the H 2 CO excitation temperature, a combination of lines with different upper level energies is needed. Since we only observed one line, the excitation temperature is fixed and the column density was computed for a range of excitation temperatures between 20 and 50 K so that we are not a priori assuming the vertical location of the H 2 CO emitting layer. The range of temperatures is based on previous estimates. For example, Pegues et al. (2020) recently showed that the bulk H 2 CO emission in four Class II disks arises in 11−40 K gas, using multi-line H 2 CO observations with similar energy ranges (E u = 20 − 70 K). In particular, for MWC 480 they found a low excitation temperature of ∼ 21 K. Similarly, toward HD 163296 the disk-averaged H 2 CO excitation temperature has been estimated to be ∼ 24 K (Carney et al. 2017;Guzmán et al. 2018 We found that the column density is not very sensitive to T ex ; the difference in column density is only a factor Figure 4. Column density profiles for HCN (top), C2H (middle), and H2CO (bottom). The horizontal bar at the upper right marks the spatial resolution of the profiles. Gray shaded areas mark the regions where the excitation temperature was assumed to be equal to the 12 CO or 13 CO temperature. Two column density profiles are shown for H2CO, assuming a fixed excitation temperature of 20 and 50 K.
∼ 2 for this range of temperatures. The observed spectra with the best-fit overlaid are shown in the bottom panels in Fig. 3. As was found for HCN and C 2 H, the lines are fit well in most of the disk, but poorly fit in the inner ∼ 30 au, especially in HD 163296 and MWC 480. This is probably due to contamination from the back side of the disk.

HCN, C 2 H, and H 2 CO column density profiles
The resulting column density profiles for HCN, C 2 H, and H 2 CO are shown in Fig. 4. Two profiles are shown for H 2 CO, corresponding to T ex of 20 and 50 K. The gray shaded area marks the inner regions where T gas ( 12 CO) or T gas ( 13 CO) was used as the excitation temperature in the fit. The horizontal bars show the spatial resolution of the profiles for each molecule and disk.
The HCN and C 2 H column densities range from 10 12 to 10 15 cm −2 , although the HCN column density may reach 10 16 cm −2 in the inner 30 au of GM Aur, HD 163296, and MWC 480. Overall, HCN and C 2 H present similar column density profiles, except for that HCN has a central component in GM Aur, HD 163296, and MWC 480, while the inner disk seems to be depleted of C 2 H in all sources but toward GM Aur. GM Aur presents a central dust cavity (see Fig. 2), which could result in enhanced C 2 H abundance in the inner disk. However, there is evidence that the total H 2 gas is also depleted inside the dust cavity (Dutrey et al. 2008;Huang et al. 2020;Bosman & MAPS team 2021b). Another difference between HCN and C 2 H is in the outer > 150 au disk. The HCN column density seems to reach a plateau of a few times 10 12 cm −2 , which is not present for C 2 H.
The H 2 CO column density profiles are consistent with the ones derived by Pegues et al. (2020) for all sources using different H 2 CO transitions and at lower angular resolution. H 2 CO is present at higher column density at larger disk radii compared to HCN and C 2 H. Although the spatial resolution for the H 2 CO column density profiles is lower than for the smaller organics, a few distinct rings can be seen in the profiles. GM Aur and AS 209 present less structure, while at least two rings or local maxima are seen for HD 163296 and MWC 480. The H 2 CO column density varies between the sources, ranging between 10 12 and 10 14 cm −2 in the inner 200 au. The largest column densities are seen toward GM Aur and IM Lup, which together with AS 209 are the coldest disks in the sample if we consider their stellar mass and  Figure 5 shows the excitation temperature profiles inferred for HCN and C 2 H from fitting the higher angular resolution 3 − 2 lines only, informed by the T ex obtained from the simultaneous fit to the 1 − 0 and 3 − 2 lines. The T gas ( 12 CO) and T gas ( 13 CO) profiles are shown in gray. At least for 3/5 sources, the C 2 H excitation temperatures are higher than those derived for HCN. HCN temperatures range from 20 to 30 K, while C 2 H temperatures range from 20 to 60 K. The C 2 H excitation temperatures in the outer disk regions (> 100 au) of AS 209, HD 163296, and MWC 480 are consistent with the temperatures derived by Bergner et al. (2019) at lower angular resolution. Excitation temperatures for HCN do not vary substantially within the disk sample. However, the inferred T ex values for C 2 H are higher for HD 163296 and MWC 480 than for the other three disks. This is consistent with disks around Herbig Ae stars being warmer than disks around lower-mass T Tauri stars, and suggest that C 2 H arises in warmer elevated disk layers compared to HCN.

HCN and C 2 H excitation temperatures
The inferred excitation temperatures are low (20 − 30 K) for both HCN and C 2 H in IM Lup and GM Aur, and do not vary substantially across the disks. Toward AS 209, the C 2 H excitation temperature profile shows a ring-like distribution similar to the column density ring seen in Fig. 4, reaching a peak T ex ∼ 40 K in one of the dust gaps. HCN, on the other hand, shows a decreasing T ex profile with radius. Toward HD 163296, the C 2 H and HCN excitation temperatures show two rings, similar to the distribution of the column densities, although the rings are more pronounced for C 2 H than for HCN. Toward MWC 480, the T ex of HCN and C 2 H decrease with radius in a smoother profile compared to HD 163296. Finally, in AS 209 and HD 163296, the T ex peaks of C 2 H spatially coincide with the column density peaks. We do not find the same correlation for the other sources or for HCN.

Location of the HCN and C 2 H emission surfaces
Chemical models predict HCN and C 2 H emission arises in elevated disk layers where UV photons drive the gas to temperatures of 40 − 50 K and dominate the chemistry Cleeves et al. 2018). However, a recent analysis of observations toward a sample of 16 disks at lower (>0. 5) angular resolution resulted in low excitation temperatures of 5 − 10 K and 10−40 K for HCN and C 2 H, respectively (Bergner et al. 2019). Bergner et al. (2019) suggested that the origin of the lower-than-expected excitation temperatures could be due to subthermal emission or beam dilution due to unresolved chemical substructure. The higher angular resolution MAPS data reveals these substructures across a sample of disks for the first time. In addition, the 1 − 0 line observations now allow us to better constrain the molecular excitation temperature. While we do find HCN excitation temperatures that are higher than inferred by Bergner et al. (2019) toward these disks, they remain relatively low (< 30 K). For C 2 H, on the other hand, we find higher excitation temperatures that are more consistent with the gas temperature at the expected height of the emission, in particular for the warmer disks HD 163296 and MWC 480, which are also disks that present the most substructure in their emission.
The inferred excitation temperatures for HD 163296 are consistent with the brightness temperature distributions of the HCN and C 2 H 3 − 2 lines found by Law & MAPS team (2021b). These brightness temperatures provide a lower limit to the gas temperature because these lines are partly optically thin. These results, com- bined with the emission heights estimated by Law & MAPS team (2021b) of z/r 0.1, suggest that the HCN emission arises in vertically deeper and colder layers than C 2 H.
For IM Lup and GM Aur, the T ex of C 2 H and HCN are very similar, which suggests that their emission arises from similar disk layers. However, given our uncertainties, we cannot rule out a vertical stratification. Indeed, detailed chemical models of IM Lup suggest that HCN arises in deeper layers compared to C 2 H. . Finally, it is worth mentioning that the C 2 H excitation temperature is lower for the colder disks around T Tauri stars, which show fewer substructures compared to the disks around Herbig Ae stars. This could be due to intrinsic lower gas temperatures or to small-scale substructure not resolved in our observations.

Comparison to CO column densities
CO is widely used to trace the gas in disks because H 2 does not have strong transitions in the low temperature range that characterizes disks. Thus, we compare the column density profiles of the small organics with that of CO. Figure 6 shows the normalized column density profiles of HCN, C 2 H, and H 2 CO, as well as the CO column densities derived from the C 18 O 2-1 line by Zhang & MAPS team (2021). In general, the HCN and CO column densities tend to be more centrally peaked than that of C 2 H, suggestive of an active warm cyanide chemistry in the inner disk. In the outer disk, HCN and C 2 H show similar radial distributions up to ∼ 150 au, where the C 2 H column density starts to drop faster than HCN and CO. There is notably good spatial correlation between the profiles toward HD 163296 between 50 < r < 150 au.
In general, the HCN and CO column densities correlate well, except in AS 209, where significant differences are observed. In particular, an anticorrelation is observed between CO and both HCN and C 2 H. Alarcón & MAPS team (2021) used chemical models to investigate the gas structure in AS 209, and showed that this anticorrelation can be explained by CO chemical processing, which produces a decrease in the CO abundance and an elevated C/O ratio (> 2) in the gas that boosts the formation of C 2 H. Bosman & MAPS team (2021a) also found that elevated C/O ratios (∼ 2) are needed to explain the observed C 2 H column densities for AS 209, MWC 480, and HD 163296. C 2 H is therefore very sensitive to the gas-phase C/O abundance ratio in the disk atmosphere Cleeves et al. 2018;Miotello et al. 2019).

Photochemistry of HCN and C 2 H
If the formation of HCN and C 2 H is the result of a UV-dominated chemistry, we might expect to see enhanced column densities inside the dust gaps. The depletion of millimeter dust at these locations should allow the deeper penetration of UV photons, and thus efficient formation of C 2 H and HCN. However, the penetration of UV photons will also depend on the geometry of the source (for example, the flaring of the disk). It is nevertheless worth comparing the structures observed in mm dust and molecular emission. The position of the dust continuum gaps, as measured by Law & MAPS team (2021a ; but see references therein), are marked by the gray dashed vertical lines in Fig. 6. We only include the gaps that are resolved at ∼ 0. 1 resolution. However, several other gaps are seen in higher angular resolution observations of the dust continuum emission (Huang et al. , 2020. The peak of the C 2 H column density seen toward IM Lup, AS 209, MWC 480, and the first peak seen toward HD 163296, coincides with the location of the innermost dust gaps, although the C 2 H structures are generally broad compared to the dust gaps. The second column density peak seen at ∼ 110 au toward HD 163296 does not spatially coincide with any of the dust gaps, revealing that there is not a universal connection between dust depletion and C 2 H or HCN formation. The weak correlation between dust structures and HCN and C 2 H chemical substructures could be due to the HCN and C 2 H emission arising in disk elevated layers, which are not always impacted by the dust millimeter substructure that trace the midplane. Moreover, recent studies have shown that the H 2 density seems to be much less altered than the dust inside the dust gaps, and chemical effects play a more important role in setting molecular abundances (e.g., Teague et al. 2018b;Alarcón & MAPS team 2021;Zhang & MAPS team 2021).
If photochemistry drives the formation of HCN and C 2 H, we could also expect to see enhanced emission in the outer disk where the gas density starts to drop allowing UV photons to penetrate the disk. The flaring of the outer disk could also help to capture more UV photons from the star. Although the HCN and C 2 H emission is compact compared to H 2 CO, some faint emission is indeed detected at large radii, far beyond the dust millimeter edge in IM Lup, GM Aur, and HD 163296. As an example, Fig. 7 shows radial structures seen in the tapered images for HD 163296, the disk that presents the most substructure in molecular line emission. A power-law stretch was used in the y-axis, to enhance the low-S/N structure seen in the outer disk. A faint ring that is seen in both HCN and C 2 H near 400 au is also seen in CN (Bergner & MAPS team 2021). The ring is not clearly seen in H 2 CO in the radial profile, but a ring structure is visible at large radii in the zerothmoment map, in particular along the disk minor axis (see Fig. 1). No millimeter dust grains are present at these larger radii, but small dust grains are present, as suggested by scattered light images (e.g., Monnier et al. 2017;Avenhaus et al. 2018;Rich et al. 2020). A small dust grain population is also present in the disk elevated layers where HCN and C 2 H emission seems to arise. This small grain population is likely regulating the HCN and C 2 H chemistry. An increased amount of UV photons will also contribute to enhanced HCN destruction in the outer disk. Indeed, Bergner & MAPS team (2021) find an increasing CN/HCN profile with radius, indicative of increased HCN photodissociation relative to CN in the lower-density outer disk (see also Guzmán et al. 2015). All this suggests that an active photochemistry is at play in the outer disk beyond the millimeter dust emission around HD 163296. Even so, as mentioned in the previous section, the C/O ratio in the gas may still be more important than UV exposure for the formation of HCN and C 2 H in the outer disk.

Links to planet formation
Several efforts have been made in recent years to find planets forming in protoplanetary disks. The presence of several planets has been proposed to explain the gaps seen in the dust continuum emission and the deviation from Keplerian rotation in CO observed in a few disks, including HD 163296 (Teague et al. 2018a;Pinte et al. 2018bPinte et al. , 2020Teague & MAPS team 2021) and MWC 480 (Teague & MAPS team 2021). While dust depletion in the midplane can be inferred from the millimeter dust continuum observations, it has been harder to infer the gaps are also depleted in gas, since all millimeter lines, including CO isotopologues, may reflect chemistry variations rather than true gas density variations. Indeed, chemistry has a nonlinear response to gas density. See Rab et al. (2020) for detailed thermochemical modeling of HD 163296 in particular. Here, we compare the locations of the planets proposed for HD 163296 and MWC 480 with the structures observed in HCN, C 2 H, and H 2 CO.
In Fig. 7, we compare the emission of HCN, C 2 H, and H 2 CO with the locations of the putative planets in HD 163296 (dashed vertical lines). Besides the two bright rings seen in the inner 150 au, two fainter rings are resolved in the outer disk for both HCN and C 2 H (see also Law & MAPS team 2021a). The planet predicted at 83 au by Teague et al. (2018a) and Isella et al. (2016) is located inside a dust millimeter gap, which coincides with the gap seen in HCN and C 2 H. This gap is also recovered in the column density profiles (Fig. 4). The correlation between the planet location and a decrease in the C 2 H column density was previously hinted by Bergner et al. (2019). Although the gaps seen for HCN and C 2 H are much shallower than what is observed for the millimeter dust, these observations suggest the gap is also depleted of gas.
Toward MWC 480, Teague & MAPS team (2021) proposed that a planet near 245 au could explain the observed CO deviations from Keplerian rotation. However, HCN and C 2 H are not detected at these large radii. At smaller radii (< 200 au), no perturbations could be confirmed in the CO gas kinematics. However, a gap is seen in HCN at around 70 au, which roughly corresponds to the location of a millimeter dust gap. The gap is not seen in C 2 H. Hydrodynamical simulations suggest that a ∼ 2.3 M J planet located at an orbital radius of ∼ 78 au could produce the observed dust gap (Liu et al. 2019). If there is a corresponding gas gap, then this suggests that HCN may be more sensitive to gas density perturbations than C 2 H, and could potentially be used as a diagnostic to identify massive planets in protoplanetary disks.

Origin of the H 2 CO emission
H 2 CO could be a better tracer of the colder gas in the outer disk than HCN and C 2 H, which seem to trace warmer elevated layers in the disks. The use of H 2 CO as a tracer of cold gas depends on its dominant formation pathway and on the efficiency of desorption mechanisms. H 2 CO can form from both gas phase and grain surface chemistry. The dominant reaction in the gas-phase at low ( 100 K) temperatures is the neutral-neutral reaction The grain surface chemistry pathway involves the successive hydrogenation of CO ices (Hidaka et al. 2004;Watanabe et al. 2004;Fuchs et al. 2009) This pathway can only be efficient in regions where CO has started to freeze out onto dust grains, i.e., beyond the CO snowline and below and around the CO snow surface, depending on whether the CO hydrogenation is faster than the CO desorption. The light blue vertical lines in the bottom panels of Fig. 6 mark the estimated location of the CO snowline in each disk as measured by Zhang & MAPS team (2021). CO snowlines are located at smaller radii in the colder disks around T Tauri stars, compared to the warmer HD 163296 and MWC 480 disks, where the CO snowline is well-resolved by our observations. The first H 2 CO ring seen toward HD 163296 and MWC 480 spatially coincides with the estimated location of the CO snowline, suggesting that CO hydrogenation is an important (if not dominant) formation pathway of H 2 CO in disks.
In the warm inner disk, H 2 CO is expected to form mainly in the gas phase, since CO cannot freeze out onto dust grains. The inner regions of the MAPS disks all show central H 2 CO depressions, although they are less pronounced for GM Aur and AS 209. This does not mean that H 2 CO is depleted in the inner disk; it could purely be a result of the observed line being more sensitive to cold gas. Indeed, higher-energy lines have shown centrally peaked emission profiles suggesting a warmer H 2 CO component is present in the inner disk Öberg et al. 2017).
In the outer disk, both gas-phase and CO hydrogenation could contribute to the overall H 2 CO emission. In Section 3.2.3 we showed that at least ∼50% of disks may show some spatial links between H 2 CO substructures and the outer edge of the continuum disk. We also observe spatial links between the H 2 CO column density profiles and the outer edge of the continuum disk, which are marked by the gray vertical lines in Fig. 6. In particular, the second ring or bump component seen in the H 2 CO column density profiles of HD 163296 and MWC 480 spatially coincides with the millimeter dust edge. A similar coincidence is seen toward IM Lup and GM Aur, where a shoulder-type structure is seen near the continuum edge. This suggests that the release of ices at large distances becomes more efficient, either through thermal desorption due to thermal inversion (Cleeves 2016), or by nonthermal desorption due to the penetration of UV photons (e.g., Podio et al. 2020;van't Hoff et al. 2020). Another possibility is that other ice fragments are released into the gas phase in the outer disk, which then boost the gas-phase chemistry enhancing the abundance of H 2 CO.
The H 2 CO traced by our observations could also form predominantly in the gas phase through reaction (6), with little contribution from the CO hydrogenation  (7). This could happen if nonthermal desorption of H 2 CO ices is not efficient and H 2 CO remains frozen onto dust grains. In this scenario, the emission would arise in more elevated and warmer disk layers. Although chemical models have shown that pure gasphase chemistry cannot account for the observed H 2 CO abundances in the outer disk around DM Tau , observations of edge-on Class I disks have shown that H 2 CO emission arises in the warmer surface layers, as revealed by X-shaped patterns in their zeroth-moment maps (Podio et al. 2020;van't Hoff et al. 2020). These observations, however, targeted the H 2 CO 3 12 − 2 11 line which has a relatively large upper-level energy of ∼ 33 K, and is therefore more sensitive to the warmer elevated disk layers. Lower-energy lines could trace colder gas present at intermediate disk layers. In order to confirm the origin of the H 2 CO emission, multiline observations are needed (e.g., Pegues et al. 2020). The results toward Class II disks found by Pegues et al. (2020) suggest a possible trend between the location of the H 2 CO emitting layer and the age of the disk, or perhaps a difference between disks around T Tauri and Herbig stars. Future multi-line observations toward a larger sample of disks with a consistent set of lines may help to elucidate whether this trend is real.

Small organics reservoir in the inner 50 au
We can estimate the amount of HCN, C 2 H, and H 2 CO present in the inner disk regions. The gas masses within 50 au are listed in Table 3. Because of the uncertainties of the fit in the inner disk, we extrapolated the column density profiles inward, assuming a constant column density starting from the radius where we obtain a good fit (see Fig. 4). This means that the HCN masses estimated for GM Aur, HD 163296, and MWC 480, which tend to have centrally peaked N profiles, could be larger than what we have estimated here. The total amount of HCN, C 2 H, and H 2 CO is, however, much larger because a large fraction of the organic reservoir is expected to reside in icy mantles. Indeed, as discussed before, the emission of these small organics arises in disk layers with temperatures of 20 − 60 K, which are lower than the freeze-out temperature of these molecules. Assuming a conservative ice-to-gas ratio of 1000 (these ratios could be larger, as discussed inÖberg & MAPS team 2021), we find a large amount (10 23 − 10 26 g) of gas phase plus ice mantle HCN, C 2 H, and H 2 CO in the inner disks. Interestingly, the relative amounts of these organics varies between the disks. For example, H 2 CO is ∼ 5 times more abundant than HCN and C 2 H in IM Lup, while H 2 CO is less abundant than HCN and C 2 H in the rest of the disks. HCN is more abundant (by a factor ∼ 10 − 50) than both C 2 H and H 2 CO in GM Aur, HD 163296, and MWC 480, and AS 209 has a slightly larger C 2 H mass than HCN.
We can also compare the amount of organics with respect to water ice in these disks. The total water mass inside 50 au was estimated from the models presented in Zhang & MAPS team (2021), and are listed in Table 1 ofÖberg & MAPS team (2021). The abundance of total gas and ice HCN, C 2 H, and H 2 CO relative to water are listed in the right columns of Table 3. The fractions range from 0.5% to 1.1% and 0.02% to 1.8% for HCN and C 2 H, respectively, not taking into account IM Lup, which has much lower organic fractions of < 0.001%. The fractions for H 2 CO are lower, ranging from 0.001% to 0.08%. Typical abundances with respect to water in cometary ices range from 0.1% to 0.6% for HCN and from 0.1 − 1% for H 2 CO (Mumma & Charnley 2011), which are consistent with our organic fraction estimates.
The MAPS disks are therefore rich in organic species, suggesting that future comets formed in these disks could efficiently deliver water and other key organics to rocky planets forming in the inner disk.

CONCLUSIONS
We have presented spatially resolved observations of the small organics HCN, C 2 H, and H 2 CO toward a sam-ple of five protoplanetary disks from the MAPS program. The main conclusions are the following: 1. The line emission reveals substantial substructure, in particular toward the warmer disks around HD 163296 and MWC 480, where several rings are observed in HCN and C 2 H.
2. We derive column density profiles and find the distributions vary across disks, suggesting that the specific physical conditions in each disk play a major role in setting the organic distribution in disks.
3. We find similarities between the HCN and C 2 H distribution in the outer disk regions, suggesting a similar chemistry drives their formation, most likely photochemistry in the temperate disk layers.
4. The HCN and C 2 H excitation temperatures, estimated by the 1 − 0 and 3 − 2 lines as well as the hyperfine structure of the lines, range between 20 and 60 K. The temperatures are higher for C 2 H, however, suggesting that the HCN emission arises in a colder and deeper layer than C 2 H, and is thus more sensitive to density variations.
5. HD 163296 and MWC 480 show gaps in their HCN column density profiles, which spatially coincide with the location of planets suggested in previous studies. This opens the possibility of using HCN to trace gas density variations in disks, and finding the massive planets that could be carving these gaps.
6. Contrary to HCN and C 2 H, H 2 CO emits from a much more extended region of the disks, far beyond the dust millimeter edge. H 2 CO is slightly more abundant in the colder disks around T Tauri stars than in disks around Herbig stars, but is particularly abundant in the GM Aur disk.
7. We find enhanced H 2 CO column densities at the CO snowline location and another column density enhancement around the millimeter dust edge.
8. We estimated the total gas plus ice mass inside 50 au for the MAPS disks, and find large amounts (10 23 − 10 26 g) of organic material, with fractions with respect to water ice that are consistent with fractions observed in comets of our Solar System.
We have shown that high angular resolution (down to ∼ 10 au) and sensitivity observations from these small organics provide a window into the C, O, and N budget of disks. The inner 50 au regions of the MAPS disks are rich in organic species, so comets forming in these disks could potentially deliver water and other key organics to nascent rocky planets. Future observations toward other disks with different stellar masses, ages, and disk properties will provide a more complete picture of the organic reservoir in disks, and the impact these properties may have on the possible outcomes of planet formation.

ACKNOWLEDGMENTS
We thank the anonymous referee for valuable comments that improved this manuscript.
This paper makes use of the following ALMA data: ADS/JAO.ALMA#2018.1.01055.L. 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. The National Radio Astronomy Observatory is a facility of the National   Table A.1. For HCN and C 2 H, the integrated line fluxes include several hyperfine components that are blended in velocity. The fluxes were computed using elliptical masks that were constructed using the disk inclinations and position angles listed in Table 1, and an outer radius based on the observed extent of the emission of each line. The uncertainties were estimated by computing integrated fluxes using the same elliptical mask but centered at 500 random positions that are free from line emission in the map. The listed uncertainties correspond to the standard deviation of these line-free integrated fluxes.
B. SIMULTANEOUS FIT TO THE 3 − 2 AND 1 − 0 LINES Figure B.1 shows the HCN and C 2 H column density and excitation temperature profiles found when fitting the 1 − 0 and 3 − 2 lines simultaneously. The spatial resolution of the profiles is 0. 3 which corresponds to the resolution of the tapered version of the 1 − 0 line images. An example with the resulting best-fit spectra found for AS 209 are shown in Fig. B.2. The weakest C 2 H 3 − 2 hyperfine line listed in Table 2 and shown in the bottom panel of the figure is detected in all disks except IM Lup.

C. COMPARISON OF HCN COLUMN DENSITY PROFILES
The MAPS data were used to derive HCN column density profiles in two other papers. Bergner & MAPS team (2021) derived CN and HCN column density profiles using a column density retrieval method similar to the one presented here. Cataldi & MAPS team (2021) also derived HCN column density profiles, but using the lower angular resolution images because the goal was to compare them to the profiles of DCN. The Band 6 DCN line is faint and was not detected at high S/N in the high angular resolution images. The resulting HCN column density profiles from these different analyses are shown in Fig. C.1. The profiles are consistent for all disks, and differences are only observed in the inner disk where the fit is known to behave poorly.  Best-fit rotational temperature and column density profiles to the 1 − 0 and 3 − 2 lines. The two upper and two bottom rows correspond to HCN and C2H, respectively. The gray curves correspond to the gas temperatures derived from the optically thick 12 CO 2-1 line (Law & MAPS team 2021b) and was used as the maximum rotational temperature allowed in the fit. Gray shaded areas mark the regions where the lines are poorly fit. The resolution of the profiles is 0. 3. Figure B.2. Deprojected spectra averaged over radial annuli in AS 209. HCN and C2H lines are shown in the top and bottom panels, respectively. The 3 − 2 lines correspond to the tapered images with 0. 3 resolution. The observed spectra are shown in black, and the best-fit are shown in color. Several hyperfine lines are shown for C2H: (a) and (b) correspond to the F = 2 − 1 and F = 1 − 0 components of the 1 − 0 line; (c) corresponds to the F = 4 − 3 and F = 3 − 2 components associated with the J = 7/2 − 5/2 group of the 3 − 2 line; (d) corresponds to the F = 3 − 2 and F = 2 − 1 components associated with the J = 5/2 − 3/2 group of the 3 − 2 line; and (e) corresponds to the F = 2 − 2 component associated with the J = 5/2 − 3/2 group of the 3 − 2 line.