Recently the Molecules with ALMA at Planetforming Scales MAPS ALMA Large Program reported a high number of lineemission substructures coincident with dust rings and gaps in the continuum emission suggesting a causal link between these axisymmetric lineemission and dustcontinuum substructures. To test the robustness of the claimed correlation we compare the observed spatial overlap fraction in substructures with that from the null hypothesis in which the overlap is assumed to arise from the random placement of lineemission substructures. Our results reveal that there is no statistically significant evidence for a universal correlation between lineemission and continuum substructures questioning the frequently made link between continuum rings and pressure bumps. The analysis also clearly identifies outliers. The chemical rings and the dust gaps in MWC 480 appear to be strongly correlated gt4 and the gaps in the CO isotopologues tend to moderately 3 correlate with dust rings.
['1300', '75', '492', 'Astrophysics - Earth and Planetary Astrophysics']
No Significant Correlation between Lineemission and Continuum Substructures in the Molecules with ALMA at Planetforming Scales Program
{'No Header': 'No Significant Correlation between Line-emission and Continuum Substructures in the Molecules with ALMA at Planet-forming Scales Program  \nHaochang Jiang ( GLYPH<139> GLYPH<10>GLYPH<12> ) , 1 Wei Zhu ( ] GLYPH<31> ) , 1 and Chris W. Ormel 1  \n1 Department of Astronomy, Tsinghua University, Haidian DS 100084, Beijing, China', 'ABSTRACT': 'Recently, the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program reported a high number of line-emission substructures coincident with dust rings and gaps in the continuum emission, suggesting a causal link between these axisymmetric line-emission and dust-continuum substructures. To test the robustness of the claimed correlation, we compare the observed spatial overlap fraction in substructures with that from the null hypothesis, in which the overlap is assumed to arise from the random placement of line-emission substructures. Our results reveal that there is no statistically significant evidence for a universal correlation between line-emission and continuum substructures, questioning the frequently made link between continuum rings and pressure bumps. The analysis also clearly identifies outliers. The chemical rings and the dust gaps in MWC 480 appear to be strongly correlated ( > 4 σ ), and the gaps in the CO isotopologues tend to moderately ( ∼ 3 σ ) correlate with dust rings.  \nKeywords: Protoplanetary disks - Astrochemistry - ISM: molecules - ISM: dust - Methods: statistical - Exoplanet formation', '1. INTRODUCTION': 'In the past decade, the Atacama Large Millimeter/submillimeter Array (ALMA) has provided us with unprecedented imagery on protoplanetary disks starting with HL Tau (ALMA Partnership et al. 2015). Ubiquitous substructures-including inner cavities, rings, gaps, plateaus, spiral arms, arcs, etc.-are shown in the dust continuum emission in protoplanetary disks surveys (e.g. Andrews et al. 2018; Long et al. 2018; Cieza et al. 2019; Francis & van der Marel 2020). Among these substructures, annular rings and gaps are particularly abundant (Huang et al. 2018; van der Marel et al. 2019), whereas nonaxisymmetric substructures are relatively rare (Andrews 2020).  \nSeveral formation theories of the dust ring and gap pairs predict that gas and dust substructures should correlate. One such mechanism is the pressure bump model, where the continuum ring is formed by dust trapping inside a local pressure maxima (e.g. Pinilla et al. 2012; Dullemond et al. 2018). Such pressure bumps  \njhc19@mails.tsinghua.edu.cn, weizhu@tsinghua.edu.cn, chrisormel@tsinghua.edu.cn  \ncan be formed by gap-opening planets (e.g. Lin & Papaloizou 1986), dead zone boundary (e.g. Flock et al. 2015), zonal flow (e.g. Bai & Stone 2014). Consequently, when the molecular line emission is optically thin the gas surface density enhancement would appear as a chemical ring, leading to an alignment of line-emission and continuum substructures. Such correlations between molecular emission lines and continuum rings are observed at the outer edges of inner cavities in many transition disks, which is thought to indicate the presence of giant planets inside the cavities (e.g. van der Marel et al. 2016, 2018b; Dong et al. 2017; Boehler et al. 2017; Fedele et al. 2017). In particular, Facchini et al. (2021) report this association in the famous two-giant-planethosting PDS 70 disk. Another mechanism that also expects a correlation between gas and dust substructures is the iceline scenario. Around icelines, evaporation (e.g. Zhang et al. 2015) or sintering of dust aggregates (e.g. Okuzumi et al. 2016) can contribute to the appearance of a dust ring, while chemical rings and gaps appear when the molecules freeze out beyond the icelines. However, whether the iceline coincides with the dust ring is still a matter of debate (Long et al. 2018; Huang et al. 2018; van der Marel et al. 2019; Zhang et al. 2021a). On the other hand, there are also mechanisms of dust  \nFigure 1. top: Fraction of line-emission substructures that spatially overlap with continuum substructures in individual disks and all disks combined.  \n<!-- image -->  \nbottom: The same as the top but only counting line-emission substructures with radial locations < 150 au. Circles are observational results. Colored violin plots show the PDF of null hypothesis (line-emission substructures follow random distribution independent of continuum). The error bars show the [15 . 9 , 50 , 84 . 1] percentiles of the probability distribution.  \nring formation, e.g., the Secular Gravitational Instability (Tominaga et al. 2020) and the Clumpy Ring Model (Jiang & Ormel 2021), which do not necessarily require that gas and dust substructures correlate spatially.  \nRecently, the Molecules with ALMA at Planetforming Scales (MAPS) ALMA Large Program allows studies on the gas component of protoplanetary disks at unprecedented resolution and sensitivity ( Oberg et al. 2021). The MAPS program surveyed approximately 50 molecular lines from 20 different chemical species around five sources-IM Lup, GM Aur, AS 209, HD 163296, and MWC 480-where dust substructures have been detected in their continuum counterparts (Long et al. 2018; Andrews et al. 2018; Huang et al. 2020). In total, more than 200 chemical substructures corresponding to certain molecular lines are identified at high spatial resolutions (7-30 au) (Law et al. 2021). Almost all of these line-emission substructures show signs of an axisymmetric structure (Law et al. 2021). In other words, rings and gaps are also the most abundant morphology in line-emission substructures. These rings and gaps are found at nearly all radii in the observed disks. Many of these line-emission substructures coincide with the dust rings and gaps in the continuum emission, suggesting possible spatial associations between chemical lines and dust continuum (Law et al. 2021). However, statistical tests are needed in order to robustly establish the physical nature of such associations. Because the sub-st  \nuctures are so abundant, a random placement of them may also lead to some level of spatial association.  \nIn this work, we quantify the apparent correlation between chemistry and continuum features by means of a null-hypothesis test. We test whether the observed apparent correlation arises solely from randomly distributing the chemical lines. We search for correlations within individual disks (Section 2.2), and also conduct test on chemical groups (Section 2.3).', '2. METHOD AND RESULT': 'A null hypothesis assumes that there is no universal spatial association between dust and line-emission substructures. If the spatial correlation is physical, then the observed overlap fractions should show a significant difference from the overlap fractions from the null hypothesis.  \nThe real line-emission distribution varies depending on the disks and the chemical species. To quantify the spatial correlation between the dust and line-emission substructures, we build two groups of statistic tests. We summarize the statistical setup as follows.', '2.1. Quantifying the Spatial Correlation': "We generally follow the approach of Law et al. (2021) in quantifying the correlation between line-emission and continuum substructures. We define an overlap as when  \nTable 1. Model Parameters  \nNote -Model parameters of null hypothesis distributions. (1) sample group name; (2) number of rings; (3) inner boundary of the ring distribution; (4) outer boundary of the ring distribution; (5) number of gaps; (6) inner boundary of the gap distribution; (7) outer boundary of the gap distribution.  \nthe radial location of the line-emission substructure, r 0 , chem 1 , falls within the width of the annular continuum substructures  \nr 0 , dust -w dust 2 < r 0 , chem < r 0 , dust + w dust 2 . (1)  \nHere r 0 , dust and w dust are the radial location and the width of the continuum substructure respectively. 2 The distance between any two dust rings(gaps) is always larger than the sum of their half-width and thus one chemical feature won't fall in two different dust rings(gaps). However, because the dust ring-gap pair could spatially overlap, one chemical feature may fall into one dust ring and one dust gap at the same time.  \nThe overlap fraction is given by  \nF (A | B) = #of feature B overlap with feature A #of feature B (2)  \nwhere features A and B are chemical rings and gaps or dust continuum rings and gaps. Four combinations are available. The observed values of F (dust | chem . )  \n(hereafter 'chem.' and 'dust' inside formula referring to all four combinations) for individual disks are shown in Fig. 1 as open circles, which reproduces the results in the Fig. 21 of Law et al. (2021). In addition, we calculate the combined overlap fractions-the sum of overlapped features among 5 disks over the total number of chemical features, which is colored in purple. The ratios range from 0 to ∼ 70%, depending on the disk and the combination of substructures of chosen. Following Law et al. (2021), we also calculate the overlap fraction among rings/gaps within 150 au, within which the majority of dust substructures are located. The results are shown as the open circles in the lower panel of Fig. 1.", '2.2. Group by disks': 'We firstly test the correlation separately for each disk. We construct the substructure distribution of the null hypothesis as follows. For each disk:  \n- 1. The chemical rings or gaps appear between a certain range [ r in , r out ]. Here r in and r out are chosen such that the boundary locations of a finite number of randomly generated substructures better match those detected in observations. As a result, these boundary values are slightly smaller (larger) than the location of the observed inner-  \nTable 2. Significance of the Observed Overlap Fractions  \nNote -The percentiles of the observed overlap fractions in randomly generated distributions. The corresponding significances are inside brackets by translation via a standard normal distribution table. (1) sample group name; (2) overlap fraction between chemical gap and continuum gap; (3) overlap fraction between chemical gap and continuum ring; (4) overlap fraction between chemical ring and continuum ring; (5) overlap fraction between chemical ring and continuum gap; (6)-(9) are the same as (2)-(5), but only counting for substructures located within 150 au.  \nmost (outermost) feature. 3 As the rings and gaps appear at different disk radii, gaps are on average closer to the host star than rings. Therefore, we take different [ r in , r out ] for rings and gaps. The values we take for each source are listed in Table 1. We tested different values and find that the results are largely unchanged.  \nlog r in = log r obs in -1 2 L N -1 (3)  \nand  \n- 2. The number density of rings or gaps follows the distribution n ( r )d r ∝ r -p d r . We take p = 1 for all disks. This power-law approximately matches the distribution of the radial locations of substructure shown in Fig. 17 of Law et al. (2021). We have also varied the slope between 0.5 and 1.5, which did not affect our conclusions.\n- 3. For each synthetic system, the total numbers of randomly generated chemical rings or gaps are equal to the numbers of rings and gaps observed in the same system by MAPS, respectively. The numbers are listed in Table 1.  \nWe randomly generate substructures following this setup and calculate the overlap fractions. We repeat this process for 10 5 times for each disk to get the probability density functions (PDFs) of the overlap fraction. The results are shown in Fig. 1.  \nWe evaluate the statistical significance of the observed overlap fractions based on the corresponding PDFs from  \nlog r out = log r obs out + 1 2 L N -1 (4)  \nfor the inner and outer boundary of the random distribution. As the sample size N increases, these numbers will thus get closer to the observed locations. We have verified that our conclusions are unaffected for any choice of the numerical prefactor of the L/ ( N -1) term in the range between 0 and 1.  \nFigure 2. top: Fraction of line-emission substructures that spatially overlap with continuum substructures in each chemical species group.  \n<!-- image -->  \nbottom: The same as the top but only counting lineemission substructures with radial locations < 150 au.  \nCircles are observational results. Colored violin plots show the PDF of null hypothesis (line-emission substructures follow random distribution independent of continuum). The error bars show the [15 . 9 , 50 , 84 . 1] percentiles of the probability distribution.  \nthe simulation. We define the statistical significance as if the probability distribution is of Gaussian shape. Specifically, a detection has a significance of < 1 σ , < 2 σ , and < 3 σ if the observed value falls within (15 . 9 , 84 . 1), (2 . 3 , 97 . 7), and (0 . 1 , 99 . 9) percentiles of the probability distribution derived from the simulation, respectively. The percentile values and the corresponding statistical significances of the observed overlap fraction relative to the simulated PDFs are provided in Table 2.  \nWe find that statistically there is no significant evidence for a universal correlation between line-emission and continuum substructures. That is, the presence of a line-emission substructure does not enhance the probability of finding a dust ring (or gap) in its vicinity. More than half of the observed overlap fractions show statistical significances of < 1 σ . For individual disks, the observed overlap fractions all fall within the 2 σ ranges of the corresponding PDFs from simulations, except for the overlap fractions, F (dust gap | chemical ring), of MWC 480 ( > 4 σ ) and GM Aur (2 . 5 σ ). When all disks are combined, there seems to be a statistically significant ( ∼ 4 σ ) correlation between chemical rings and dust gaps. However, such a detection is primarily driven by the single  \nFigure 3. top: Fraction of continuum substructures that spatially overlap with line-emission substructures in each chemical species group.  \n<!-- image -->  \nbottom: The same as the top but only counting lineemission substructures with radial locations < 150 au.  \nCircles are observational results. Colored violin plots show the PDF of null hypothesis (line-emission substructures follow random distribution independent of continuum). The error bars show the [15 . 9 , 50 , 84 . 1] percentiles of the probability distribution.  \nsystem of MWC 480. Specifically, there are 13 chemical rings located inside the D76 dust gap in MWC 480, making up one third of the (chem. ring, dust gap) associations of all five disks. If the MWC 480 system is removed from the analysis, the statistical significance of F (dust | chem . ) is reduced to below 3 σ . Therefore, we conclude that there is no universal correlation between chemical rings and dust gaps.', '2.3. Group by chemical species': 'Even though line-emission substructures as a whole are distributed randomly for most disks, it is still possible that specific molecular lines may correlate with dust rings or gaps more strongly than others. For example, if a certain molecule X triggers the formation of a continuum ring or gap, it may not easily show up in the previous test where substructures from species X are blended with substructures from all other species. We therefore further investigate the F (dust | chem . ) by grouping the data with different chemical species. Following Law et al. (2021) (see their Section 3.5), we select and group the chemical species into three groups depending on their chemical similarity-CO isotopologues, nitriles, and hydrocarbons. We collect the chemical fea-  \ntures depending on their species group from every disk. Each chemical group is assumed to follow its own random distribution. Therefore, the generation of synthetic systems needs to be modified:  \n- 1. For each group of chemical species, the spatial range of the randomly placed substructure spans from the inner- and outer-most locations where substructures from this group of species are detected regardless of disks.\n- 2. For each disk, the number of randomly generated rings/gaps is the same as the number of observed rings/gaps associated with the same group of species. The substructures from different disks are then combined to compute the overlap fraction of the given group of species.  \nThe model parameters are summarized in Table 1. We first calculate F (dust | chem . ) and show the results in Fig. 2. Similar to Fig. 1, the open circles indicate the fractions from observations, the violin plots are used to illustrate the PDFs of the fractions derived from the simulation, and the error bars indicate medians and 1 σ ranges. The percentiles of observed overlap are also reported in Table 2.  \nAs shown in Fig. 2 and listed in Table 2, the majority of the observed overlap fractions of the three groups of species are statistically insignificant ( < 3 σ ), regardless of the combinations of substructures. The only exception is the chemical gap and dust ring alignment in CO isotopologues, which has a statistical significance of 3 . 3 σ (See Sect. 3). As has been discussed previously, the majority of the line-emission substructures in MWC 480 coincide with the D76 dust gap. This single outlier leads to the marginal detections ( ∼ 2 σ ) of correlations between dust gaps and chemical rings associated with all three groups of species.  \nBy combining the substructures of the same group of species from different disks, we now have a large enough sample to assess the overlap fraction F (chem . | dust). This directly informs the fraction of dust substructures associated with certain chemical species. The results are shown in Fig. 3 and listed in Table 2.  \nThe observed overlap fractions are statistically insignificant ( glyph[lessorsimilar] 2 σ ), regardless of the combination of substructures or the group of chemical species. In other words, the presence of a continuum substructure does not significantly enhance, or decrease, the probability of observing a line-emission substructure in its vicinity. The results are largely unaffected when only substructures within 150 au are included. This analysis strengthens our conclusion that line-emission and dust substructures are generally independent.', '3. CONCLUSION AND DISCUSSION': "In this work, we revisit the spatial correlation between line-emission and continuum substructures in the MAPS survey. We have calculated the observed overlap fractions of different combinations of substructures for individual disks and individual group of chemical species. We assess the statistical significance by comparing the observed fraction with that from synthetic systems that assume no correlation. Statistically, the current observational results suggest that the line-emission substructures distribution and dust substructures distribution are overall independent of each other. No significant line-emission and continuum substructures correlation is found among the MAPS disks, regardless of accounting for the entire disk or limiting our analysis to the inner 150 au. We report our key findings and briefly discuss them as below:  \n- 1. No statistically significant correlation is found between the locations of molecular line emission and dust substructures in at least four out of the five MAPS disks, regardless of which combination of substructure types is considered. In particular, a positive correlation between dust rings and chemical rings is absent. Adopting the simplifying assumption that pressure bumps reflect gas surface density maxima, we would naively expect concomitant peaks in the line emission. Although in reality the actual physical-chemical response to a pressure bump may be more complex such that pressure maxima do not necessarily lead to local chemical emission lines, the lack of a clear continuum-line correlation may alternatively suggest that pressure maxima are perhaps not always the driving mechanism behind continuum rings. If this interpretation is right, it would imply that the gas density profile is smooth around continuum and line substructures. Indeed Alarc'on et al. (2021) found that for the AS 209 disk the H 2 gas likely follows a smooth surface-density profile even though there are substructures in certain chemical species.\n- 2. Among the five MAPS systems, MWC 480 is a clear outlier with a significant correlation ( > 4 σ ) between chemical rings and dust gaps. This strong correlation is primarily due to the large accumulation of line-emission substructures ( > 13 rings) at its D76 dust gap. Active chemical processing, e.g., gas flow onto the accreting planet(s) (e.g. Liu et al. 2019; Dong et al. 2019; Teague et al. 2019), could be ongoing inside the D76 dust gap in MWC 480.  \n- 3. In addition, a weak ( ∼ 2 σ ) trend of association between line-emission rings and dust gaps appears when focusing on CO isotopologues. This meets the scenario demonstrated by some physicalchemical modeling, where the CO emission is increased inside the gaps due to the combination of increased temperature (e.g. van der Marel et al. 2018a; Kim et al. 2020) and contributions of the backside of the disk (e.g. Rab et al. 2020). Specifically, Rab et al. (2020) reproduce this association for both CO 2-1 and C 18 O 2-1 emission around the D49 continuum gap in HD 163296. Similar processes may operate in GM Aur to produce the observed trend( ∼ 2 . 5 σ ) between dust gaps and CO isotopologues rings.\n- 4. We report a moderate ( ∼ 3 σ ) correlation between line-emission gaps of CO isotopologues and continuum rings. One possible general explanation has been highlighted by Law et al. (2021)-the overdense regions of dust particles may absorb the optically thick lines and lead to gaps in molecular line observations (see also Weaver et al. 2018). A naive expectation of this mechanism would be that the same correlation should get stronger with closerin substructures, because the regions closer to the star should preferentially be optically thicker (Bosman et al. 2021). This is not observed in our results. In fact, the significance of the same correlation decreases when we only count substructures within 150 au. Another explanation is that excess cooling of grains inside the dust ring can cause a temperature dip and may lead to a lineemission gap (see Sec. 5.2 of Zhang et al. 2021b). Furthermore, optically thick dust rings can absorb  \nthe line emission from the back side of the disk (Rab et al. 2020). However, both of these two mechanisms should apply to other species rather than only CO. Finally, with the concentration of pebbles at the midplane, more efficient physical sequestration of CO onto pebbles and local chemical processing may happen inside continuum rings, which could contribute to such a correlation (See Krijt et al. 2020 for a more comprehensive study on CO depletion). More data is needed to confirm this correlation and further analysis and modeling are required to understand its cause.  \nGiven the limitations of current resolution and sample size, one cannot rule out the possibility that some spatial correlations between chemical and dust features may hold for specific disks or specific chemical species. Better conclusions require higher resolution and a larger sample. There is increasing evidence that the current chemical inventory in disks may be strongly affected by dust transport, which will lead to significant correlation between chemical components and continuum disk globally (e.g. Bosman & Banzatti 2019; Banzatti et al. 2020; van der Marel et al. 2021). Yet, at present, a universal local continuum-emission-line correlation does not stand.  \nWe thank the referee, Nienke van der Marel, for her thoughtful and constructive comments, which significantly improve the quality of this manuscript. H.J. appreciates helpful discussion and comments from Gregory J. Herczeg, Sebastiaan Krijt, Charles J. Law, and Feng Long.  \nSoftware: Matplotlib (Hunter 2007), Scipy (Virtanen et al. 2020), Numpy (Harris et al. 2020)", 'REFERENCES': "Alarc'on, F., Bosman, A. D., Bergin, E. A., et al. 2021, ApJS, 257, 8, doi: 10.3847/1538-4365/ac22ae ALMA Partnership, Brogan, C. L., P'erez, L. M., et al. 2015, ApJL, 808, L3, doi: 10.1088/2041-8205/808/1/L3 Andrews, S. M. 2020, ARA&A, 58, 483, doi: 10.1146/annurev-astro-031220-010302 Andrews, S. M., Huang, J., P'erez, L. M., et al. 2018, ApJL, 869, L41, doi: 10.3847/2041-8213/aaf741 Bai, X.-N., & Stone, J. M. 2014, ApJ, 796, 31, doi: 10.1088/0004-637X/796/1/31 Banzatti, A., Pascucci, I., Bosman, A. D., et al. 2020, ApJ,  \n903, 124, doi: 10.3847/1538-4357/abbc1a  \n```\nBoehler, Y., Weaver, E., Isella, A., et al. 2017, ApJ, 840, 60, doi: 10.3847/1538-4357/aa696c Bosman, A. D., & Banzatti, A. 2019, A&A, 632, L10, doi: 10.1051/0004-6361/201936638 Bosman, A. D., Bergin, E. A., Loomis, R. A., et al. 2021, ApJS, 257, 15, doi: 10.3847/1538-4365/ac1433 Cieza, L. A., Ru'ız-Rodr'ıguez, D., Hales, A., et al. 2019, MNRAS, 482, 698, doi: 10.1093/mnras/sty2653 Dong, R., Liu, S.-Y., & Fung, J. 2019, ApJ, 870, 72, doi: 10.3847/1538-4357/aaf38e Dong, R., van der Marel, N., Hashimoto, J., et al. 2017,\n```  \nApJ, 836, 201, doi: 10.3847/1538-4357/aa5abf"}