Open Access
Issue
A&A
Volume 695, March 2025
Article Number A202
Number of page(s) 10
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202553891
Published online 24 March 2025

© The Authors 2025

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

M 82 is one of the nearest starburst galaxies (∼3.6 Mpc; Freedman et al. 1994), experiencing strong galactic winds and serving as a superb analog of high-redshift star-forming galaxies. Detections of extraplanar polycyclic aromatic hydrocarbons (PAHs) in M 82 provide a new window to uniquely observe the cool pockets of interstellar medium (ISM) material, which can be ejected from the disk through galactic winds and is capable of surviving the harsh circumgalactic medium (CGM) environment (e.g., Bolatto et al. 2024; Chastenet et al. 2024). Consequently, starburst outflows can potentially explain quenching mechanisms shutting down the production of new stars in galaxies by either gas depletion or removal (e.g., Oppenheimer et al. 2010; Page et al. 2012).

In extreme environments such as actively star-forming galaxies at low (e.g., Zschaechner et al. 2018) and high redshift (e.g., Swinbank et al. 2011), unbound molecular gas not immediately available for star formation and not traced by CO emission could dominate both the mass and volume of the ISM (e.g., Papadopoulos et al. 2018). Moreover, ionization cones of active galactic nuclei (AGN) can be faint in CO even in the presence of significant emission of warm rotational (and also ro-vibrational) H2 lines (e.g., Davies et al. 2024). This could produce severe biases on the derived H2 content in outflows from starburst galaxies, hiding the real impact of galactic winds on quenching the star formation activity due to gas removal.

In these scenarios, mid-infrared (mid-IR) data due to PAH emission are useful to account for the H2 not traced by CO rotational transitions. PAH mid-IR emission arises from small grains heated mostly by ultraviolet (UV) photons as a consequence of young stars producing strong radiation fields near sites of star formation. In particular, PAH emission is well represented by features at λ = 7.7 μm and λ = 11.3 , μm (e.g., Tielens 2008; Draine 2011) with dust capturing UV photons and reprocessing them at IR wavelengths. Consequently, the dust IR emission has been shown to correlate tightly with Hα and UV (e.g., Kennicutt & Evans 2012; Calzetti 2013). In addition, the emission from dust reflects the total surface density (i.e., dust+gas+stars) of the region where dust is located. Since star formation and cold gas tend to track one another well at large scales, the PAH emission has also been widely used as a tracer of the gas mass (atomic and molecular) or the gas distribution in galaxies (e.g., Gao et al. 2019; Chown et al. 2021; Leroy et al. 2021a). Several studies have demonstrated the close relationship between different CO rotational transitions and PAH emission features along the mid-IR in low (e.g., Regan et al. 2006; Leroy et al. 2015; Alonso-Herrero et al. 2020; Leroy et al. 2023) and high-z galaxies (e.g., Pope et al. 2013; Cortzen et al. 2019). Moreover, spatially resolved studies have shown that PAH emission arises primarily from diffuse ISM and H II regions; however, only a small fraction of PAH emission could be originating from the ionized gas (indicating rapid processing and/or destruction of PAHs; e.g., Sutter et al. 2024). PAH emission also disappears in very dense regions where they may congregate in larger grains (e.g., Chastenet et al. 2023), thus do not produce the characteristic broad, mid-IR features since the UV that illuminates the PAH features is shielded away. Since PAHs are bright and detectable at high redshift (e.g., Spilker et al. 2023; Shivaei et al. 2024; Shivaei & Boogaard 2024), understanding how the properties mentioned above and their abundance changes in local galaxy populations yields important gains for the study of the distant universe.

M 82 provides an ideal laboratory to verify how CO rotational transitions couple with mid-IR emission from PAHs in an extreme environment. In this work, we test the relations between CO(J = 1–0) emission line and the λ = 7.7, 11.3 μm PAH features using archival NOrthern Extended Millimeter Array (NOEMA) and new James Webb Space Telescope (JWST) data for M 82. To do so, we focus on the analysis of CO clouds since these structures are usually related to self-shielded, dense molecular clouds capable of surviving in outflow streams. Although is not yet clear if H2 clouds are actually able to move with or only exposed to outflows, if these structures are capable to escape with the material may carry out a significant fraction of molecular gas available for star formation in galaxies (i.e., Recchi & Hensler 2013; Bezanson et al. 2022), leading eventually to star-formation quenching in the post-starburst stage (e.g., Suess et al. 2017). The comparison between the scaling relations we obtain for M 82 and those for low-z main-sequence field galaxies at sub-kpc scales could provide helpful information about how the heating mechanisms of PAHs change in extreme environments and valuable constraints for future numerical simulations.

The present work is structured as follows: Section 2 presents the observations and the processing of NOEMA and JWST data. In Section 3 we explain the methods applied to analyze the data and the equations used to derive the key physical quantities. In Section 4 we present our results and discussion, and in Section 5 we summarize the main conclusions.

2. Observations and data reduction

2.1. JWST data

The F770W and F1130W images (centered on the λ = 7.7 μm and λ = 11.3 μm PAH features, respectively; see Fig.1) for M 82 used here were acquired with JWST as part of the Cycle 1 GO project 1701 (P.I. Alberto Bolatto; see Bolatto et al. 2024 for more details) using the MIRI iMager (MIRIM). The details of the data synthesis will be presented in Cronin et al. (in prep.); in this section, we only summarize the main features of the data.

thumbnail Fig. 1.

From left to right: M 82 NOEMA CO(J = 1 − 0), MIRI-F770W, and MIRI-F1130W images in cutouts of 5′×4.2′. The left panel also contains the molecular clouds identified in Section 3.2 (white contours), the centroid of the CO(J = 1 − 0) and the different regions identified in the starburst as described in Krieger et al. (2021). The inset in the left panel correspond to the NOEMA beamsize at ν = 115 GHz. The green crosses in the middle and right panels are indicate the optical center of M 82. The black contours in the middle and right panels correspond to the regions covered by the MIRI-MRS and the yellow boxes are the regions with Spitzer IRS data from Beirão et al. (2008); both are used in this work to apply a first-order continuum subtraction (see Section 3.1).

The data presented here are a combination of images taken on 2023-12-31 between 12:59:59 and 15:28:19 UT in two modes. We used a 4 × 1 mosaic of the FULL MIRIM detector array along the minor axis of the galaxy. These images used a 4-point CYCLING dither and the FASTR1 readout pattern for a total exposure time of 333 s per filter. In addition, we required a single SUB128 pointing to avoid saturation over the center of the galaxy. These subarray images used a 4-Point-Sets dither optimized for extended sources and the FASTR1 readout for a total exposure time of 2.4 s per filter. The final subarray images yield F770W and F1130W maps (see middle and right panels of Fig. 1) with angular resolutions of 0 . 4 $ {\sim}{0}{{\overset{\prime\prime}{.}}}4 $ and 0 . 5 $ {\sim}{0}{{\overset{\prime\prime}{.}}}5 $ (given by the FWHM of the PSFs in each case).

We note that the corresponding MIRI background pointings failed to execute due to a guide star acquisition problem, and have not yet been repeated. We, therefore, cannot perform a global background subtraction for these data. Using the JWST background model (JWST_BACKGROUND1) the background levels are predicted to be 7 MJy sr−1 for F770W and 21 MJy sr−1 for F1130W at the time of these observations. These corrections are applied to the data used in this work.

2.2. CO NOEMA data

We use the 12CO(J = 1 − 0) observations presented in Krieger et al. (2021). The data comprise a 25 arcmin2 NOEMA mosaic, covering a physical area ∼7.9 × 2.9 kpc2. The observations were designed to include the disk and outflow regions similar to those covered by single dish observations that were analyzed by Leroy et al. (2015). The main features of the CO NOEMA data and their processing are described in Krieger et al. (2021); here we summarize their most important details for the purpose of the present work.

We use the moment 0 maps generated by Krieger et al. (2021) for the NOEMA CO(1–0) emission line. The latter was derived from a cube with a uniform beamsize of 2 . 08 × 1 . 65 $ {2}{{\overset{\prime\prime}{.}}}08 \times {1}{{\overset{\prime\prime}{.}}}65 $ (with a position angle PA = 51°), a spectral coverage of ∼800 km s−1 centered at the CO(1–0) emission line, and a channel width of 5 km s−1. The median rms noise per pixel in a 5.0 km s−1 wide channel of the resulting cube was 138 mK (5.15 mJy beam−1), achieving a noise level of 100–150 mK across the entire map (see Krieger et al. 2021 for more details).

The resulting CO(1–0) moment 0 map (left panel of Fig. 1) was produced from the cube discussed above after blanking channels below a signal-to-noise ratio SNR = 5 threshold, and collapsing the cube in the range between 130 and 375 km s−1 (i.e., within a spectral window centered at the line peak around ∼210 km s−1). We make use of this map to extract the CO integrated intensities for the clouds identified in Section 3.2.

3. Methods and products

3.1. Mid-IR continuum subtraction and band corrections

We implement a first-order continuum subtraction correction to the MIRI-F770W and F1130W filter intensities of the CO clouds identified in M 82 (see Section 3.2) using the MIRI medium resolution spectroscopy (MRS) data available for a fraction of the disk region (see black rectangles in the middle and right panels of Fig. 1). The data are also part of the Cycle 1 GO project 1701 (P.I. Alberto Bolatto). The details of this analysis are described in Appendix A. We find that the contribution of the continuum in the disk ranges from ∼30 to 53% of the total intensity at 7.7 μm (with a median of ∼48 ± 5%), and at 11.3 μm ranges from ∼7 to 35% of the total intensity (with a median of ∼35). In order to estimate the continuum contribution in the regions outside the disk, we apply a similar procedure using the spectra available for regions A and B in Figure 1 included in Beirão et al. (2008), which correspond to data taken with the Infrared Spectrograph spectrometer (IRS; Houck et al. 2004) on board of the Spitzer Space Telescope. In these regions, we find that the continuum contribution accounts for ∼28% and 39% of the total emission in the F770W and F1130W filter maps, respectively. To make sure that this approach is accurate enough, we compare our first-order continuum estimations with those derived by using the Continuum And Feature Extraction tool, CAFE (Díaz-Santos et al., in prep.; Marshall et al. 2007) for six different location within M 82, including the region A+B and five different spaxels in the region covered by the MIRI-MRS. These estimations are in agreement with our first-order approach, yielding continuum contributions accounting for ∼35% and 47% of the total emission in the F770W and F1130W filter maps, respectively, in region A+B, and average continuum contributions of 48% and 41% in the F770W and F1130W filter maps, respectively, for the MIRI-MRS spaxels. We therefore apply our first-order corrections to the MIRI intensities of clouds located in the streamers east and west, and outflows north and south (see Fig. A.2 for more details).

In addition to the continuum emission, the PAH 7.7 and 11.3 μm features are conformed by several Drude profiles with centers at 7.42, 7.60, 7.85 μm and at 11.23 and 11.33 μm, respectively (Smith et al. 2007). PAH features other than those at 7.7 and 11.3 μm can thus contribute significantly to the total emission of the MIRI-F770W and F1130W bands. To mitigate these effects, we apply the corrections for the MIRI-F770W and F1130W bands included in Donnelly et al. (2025), which are derived using CAFE. Donnelly et al. (2025) estimate the missing contribution of the broad wings of PAH profiles in each band (Cwing, 7.7 μm and Cwing, 11.3 μm), and the bandwidth and shape of the PAH-dominated filter (W7.7 μm and W11.3 μm; see also Gordon et al. 2022). To obtain the final intensities, we use the following equation:

I MIRI , j / [ MJy sr 1 ] = W MIRI j C wing , MIRI j I cont _ subt , MIRI j / [ MJy sr 1 ] , $$ \begin{aligned} I_{\mathrm{MIRI}, j}/[\mathrm{MJy\,sr^{-1}}] = W_{\mathrm{MIRI}j} C_{\mathrm{wing,MIRI}j} I_{\mathrm{cont\_subt, MIRI}j}/[\mathrm{MJy\,sr^{-1}}], \end{aligned} $$(1)

where Icont_subt, j are the intensities of the PAH features after applying the continuum subtraction, and j = 770, 1130. As prescribed by Donnelly et al. (2025), we adopt the values Cwing, 7.7 μm = 1.15, W7.7 μm = 0.217 and Cwing, 11.3 μm = 1.16, W11.3 μm = 0.066 for the MIRI-F770W and F1130W bands, respectively. The final intensities derived from Equation (1) are those used in this work.

3.2. Cloud identification

We use QUICKCLUMP2 (Sidorin 2017) to identify molecular cloud structures in the NOEMA CO(J = 1 − 0) moment 0 map in the inner ∼2 kpc from the central starburst, which is based primarily on DENDROFIND (Wünsch et al. 2012; widely used by studies on molecular clouds; e.g. Rosolowsky et al. 2008; Indebetouw et al. 2020). The details of this analysis are described in the appendices. Clouds can be identified in the three-dimensional distribution along the CO NOEMA datacube (i.e., identifying their RA, declination, and velocity; Krieger et al. 2021). However, we perform their identification using the NOEMA moment 0 map; although this could introduce biases in the analysis since some of the clouds may overlap along a common line-of-sight, we implement this methodology in order to match the 2D-projected properties of both CO clouds and PAH structures in the MIRI F770W and F1130W maps. Since the PAH maps lack in velocity information, we look for cloud footprints that can be applied to both NOEMA and MIRI data. We validate 306 CO clouds (white contours in the left panel of Fig. 1), with sizes ranging between ∼31 and 270 pc (see Fig. 2). When we compare the main features of our CO clouds to those reported in Krieger et al. (2021), who characterize them using the NOEMA CO(1–0) datacube and identify ∼2000 CO clouds, we note a good agreement in the size distributions. They find CO clouds with radii in the range of 40 − 60 pc, noting that larger clouds are found in the disk (median ∼60 pc, with values up to ∼150 pc) compared to the outflows and streamers (medians 40 − 45 pc).

thumbnail Fig. 2.

Distribution of sizes of the 306 CO clouds we identify in the NOEMA map ( θ mean = 1 . 8 31 $ \theta_{\mathrm{mean}}={1}{{\overset{\prime\prime}{.}}}8 \approx 31 $ pc) using QUICKCLUMP as described in Section 3.2. The red-dashed line marks the median size for the cloud sizes. The inset includes the distributions of clouds respect to the five regions we use to characterize M 82 (see Section 3.2).

Once we identify the clouds, we compute both their integrated CO and the JWST-MIRI F770W and F1130W filter intensities. First, we convolved the F770W and F1130W maps to match the angular resolution of the NOEMA CO(1–0) moment 0 map using the PYTHON package astropy.convolution. Then, we regrid the JWST-MIRI F770W and F1130W maps to also match the pixel size of the NOEMA CO moment 0 map ( 0 . 4 $ {0}{{\overset{\prime\prime}{.}}}4 $). Next, we sum up the intensity contribution in each map of all the pixels contained in a given CO cloud. We use the conversion factor given by Krieger et al. (2021) for the CO(J = 1 − 0) line (i.e., Jpk ≈ 26.8 Jy K−1). In addition, we also compute the projected distance of the clouds to central starburst of M 82 using its optical center taken from the NASA-IPAC Extragalactic Database3 (see green crosses in Fig. 1). Finally, we adopt the region characterization for M 82 performed by Krieger et al. (2021) to classify the clouds in five different zones (delimited by yellow lines in the left panel of Fig. 1): disk (DISK; 22 clouds), streamer east (SE; 46 clouds), streamer west (SW; 51 clouds), outflow north (ON; 83 clouds), and outflow south (OS; 104 clouds). As mentioned above, even if we refer to these structures as “clouds”, we warn that the choice of deconvolving the space and geometry of M 82 (with inclination i ∼ 80°; McKeith et al. 1993) in a 2D projection introduces a bias in the total number of structures that can be identified. Although we prefer to call them with the more generalized term “structures”, we use the term “clouds” hereafter for simplicity.

3.3. Control sample

Chown et al. (2024) conduct a detailed comparison between CO(2–1) data from the Atacama Large Millimeter-submillimeter Array (ALMA) and 7.7, 11.3 μm data from JWST-MIRI for 66 galaxies selected from the Physics at High Angular resolution in Nearby Galaxies survey (PHANGS; Leroy et al. 2021b). The sample consists primarily of star-forming main-sequence spirals (log(SFR/[M yr−1]) > 2.8), with relatively high total stellar masses (log(M/M) > 9.73). Hereafter, and for simplicity, the latter sample and the results of its analysis included in Chown et al. (2024) are referred to as C24.

To properly compare with results from C24, we use R21 = LCO(2 − 1)/LCO(1 − 0) = 1.0 by Weiß et al. (2005) (which corresponds to the average value found across the streamer and outflows of M 82; Walter et al. 2002) to scale the relationships included in Chown et al. (2024) from LCO(2 − 1) to LCO(1 − 0). We use these resulting relations to compare with those we obtain in this work using the NOEMA CO(1–0) moment 0 (beamsize of 2 . 08 × 1 . $ {2}{{\overset{\prime\prime}{.}}}08 \times {1}{{\overset{\prime\prime}{.}}} $65) and the MIRI F770W and F1130W maps.

4. Results and discussion

4.1. CO versus PAH relations

We compute the relations between CO and mid-IR PAH emission for the 306 molecular clouds identified in Section 3.2. For the two relations, we perform an ordinary least square (OLS) linear fitting using the model y = mx + b (in log-log scale). The top left panel of Figure 3 shows the CO–F770W relation, colorcoded by the region in which clouds are located. The best-linear fit relation for the data binned in IMIRI770 values, is given by:

thumbnail Fig. 3.

ICO(1 − 0) versus IMIRI770 (top), and IMIRI1130 (bottom). Points in left panels are colored in the five regions we adopt to classify the clouds in M 82 (see Section 3.2). We also include the best-fit linear relations we obtain for the whole sample after binning the points in the x-axis and performing an ordinary least square (OLS) linear fitting using the model y = mx + b (black-dashed lines), and the Pearson’s r-value for the different regions in the right-bottom corners. Error bars correspond to the scatter of the points in the y-axis for the binned data. The left panels show the same relations that those in right ones, but colorcoded by the size of clouds, where symbols are split in clouds larger (circles) and shorter (stars) than 100 pc.

log ( I CO ( 1 0 ) / [ K km / s ] ) = ( 0.74 ± 0.19 ) × log ( I MIRI 770 / [ MJy sr 1 ] ) ( 1.08 ± 0.39 ) . $$ \begin{aligned} \log (I_{\rm CO(1{-}0)}/[\mathrm{K \, km/s}])&= (0.74 \pm 0.19)\times \log (I_{\rm MIRI770} /[\mathrm{MJy \, sr}^{-1}])\nonumber \\&- (1.08 \pm 0.39). \end{aligned} $$(2)

We note that the slope for the entire cloud sample (m = 0.74 ± 0.19) is lower than that in C24 (m = 0.93 ± 0.05; dashed-doted gray line in the top panel of Fig. 3). When looking at the Pearson rp values we obtain for the cloud sample divided by regions, we note moderate correlations between CO and F770W emission in clouds belonging to the disk, streamer west, and north and south outflows (with Pearson rp covering the range 0.49 − 0.81, p-value ≪0.01; see Table 1). However, the CO–F770W relation for clouds from the streamer east show a noticeable poor correlation (Pearson rp = 0.14; see Table 1). We warn that the moderate correlations in some region could be due to the geometry of M 82; the latter gives us access only to the almost edge-on portions of M 82.

Table 1.

Correlation coefficients for the relations between CO and F770W, F1130W, and F770W/F1130W.

Similarly, the bottom left panel of Figure 3 shows the ICO(1 − 0)–F1130W relation for the molecular clouds detected. We find the best-linear relation for the data binned in IMIRI1130 values as it follows:

log ( I CO ( 1 0 ) / [ K km / s ] ) = ( 0.81 ± 0.25 ) × log ( I MIRI 1130 / [ MJy sr 1 ] ) ( 0.87 ± 0.42 ) . $$ \begin{aligned} \log (I_{\rm CO(1-0)}/[\mathrm{K \, km/s}])&= (0.81 \pm 0.25)\times \log (I_{\rm MIRI1130} /[\mathrm{MJy \, sr^{-1}}])\nonumber \\&- (0.87 \pm 0.42). \end{aligned} $$(3)

When comparing this result with that found for the CO–F770W relation, we note that the CO–F1130W best fit has a similar slope (m = 0.81 ± 0.25; although also below unity). Nevertheless, we remark again that such a slope is still significantly below than that found for star-forming spirals in C24 (m = 1.00 ± 0.08). Likewise, when splitting the clouds into the regions of M 82, we note that the CO–F1130W relation show a moderate or strong correlation for all the groups (Pearson rp = 0.47 − 0.81; p-value ≪0.01).

The b parameter of the best-linear relations derived above, in combination with the m values, can provide a more panoramic view of the extreme conditions in M 82. This parameter is a more stable result than the slope (in terms of the fitting routines used) and helpful to identify PAH emission mechanisms in different contexts. When analyzing the m values as an indicator of the CO per PAH emission for the two relations, we note that CO in PHANGS galaxies is ∼25% brighter (with respect to PAH emission) than in M 82. This is also confirmed when we compare the x0 = −b/m values (i.e., the x value where the best linear relation intercepts the x-axis) for the two relations; we find that M 82 has x0 values ∼1–1.5 orders of magnitude larger than those for PHANGS galaxies. Although these results suggest an extreme contrast in such two cases, they appear to mainly reflect differences in the mechanisms of the PAH emission. Comparing the integrated intensities of the 3.3 μm PAH feature and the Paschen-α line (Pa-α), Fisher et al. (2024) find that the Pa-α per PAH emission in M 82 is ∼5–6 times brighter than in PHANGS galaxies. Such results are consistent with our findings, and they could be indicating a brighter PAH emission per mass (across several features) in starburst environments (and maybe outflows) than in nearby star-forming galaxies. We warn, however, that the best-linear fit parameters we find for the CO–PAH relations are strongly limited by the sensitivity of the NOEMA CO moment 0 map (∼0.2 K km/s). The MIRI-F770W and F1130W data are more sensitive and have a larger dynamical range than the CO maps, so their combination has the potential to bias the resulting slopes. Moreover, although we have access to (comparatively) lower PAH intensities, we do not have access to CO intensities lower than the noise threshold used to define the CO clouds.

We find that, on average, the slopes for the relations we analyze in this work are lower than what is available in the literature for normal, star-forming galaxies. However, the scaling relations between CO and mid-IR have been shown to be dependent on the CO rotational transition used. For instance, Leroy et al. (2023) note that the best-fit power laws using CO(1–0) data have slightly shallower slopes compared to those for CO(2–1) when comparing CO–PAH relations for four PHANGS galaxies, with typical offsets of ∼0.2 − 0.3. In this sense, we obtain slopes lower than those derived in C24; this is not only the case for the CO–F770W and CO–F1130W relations, but also in the slopes offset for such relations (∼0.2). Similarly to C24, Gao et al. (2019) find a slope slightly below unity (m = 0.98) when computing the CO(1–0)–L12 μm relation for 31 galaxies selected from the Mapping Nearby Galaxies at Apache Point Observatory survey (MaNGA; Bundy et al. 2015). Several studies have also found that this offset depends on the local SFR surface densities (i.e., IR luminosity densities since ΣSFR ∝ ΣLIR), which reflect changes in the CO(2–1)-to-CO(1–0) luminosity ratio described by R 21 Σ SFR 0.15 $ R_{21}\propto \Sigma^{0.15}_{\mathrm{SFR}} $ (e.g., den Brok et al. 2021; Yajima et al. 2021; Egusa et al. 2022; Leroy et al. 2022; Keenan et al. 2024). M 82 hosts a large number of compact, massive stellar clusters representative of a significant fraction of the total star formation activity (e.g., McCrady et al. 2003; Smith et al. 2006; Mayya et al. 2008; Levy et al. 2024). Consequently, the difference in the CO(1–0) and CO(2–1) power laws may be reflecting the dissimilar physical conditions of the gas in the different regions of M 82, which could vary strongly with the sizes and locations of the CO clouds. Indeed, Lopez et al. (in prep.) find that the point where the best-power laws intercept the x-axis for the CO–F770W relation in M 82 depends strongly on vertical distance and whether the emission comes from the north or south outflows. In these opposite regions, they show that the differences in the intercept can vary even up to 1.0 dex (m ∼ 0.6 − 1.6).

Interestingly, the right panels of Figure 3 show that when splitting our sample in two size groups, we find a stronger correlation in the CO–F770W (top) and CO–F1130W (bottom) relations in clouds larger than 100 pc (Pearson rp ≳ 0.7) compared to that of those with shorter diameters (Pearson rp ≲ 0.4). This result reflects that larger clouds have more significant CO detections compared to those for smaller structures, therefore covering a broader CO intensity range not limited by the sensitivity of the NOEMA CO moment 0 map. A more careful analysis of the CO-F770W relation will be performed in Lopez et al. (in prep.), who based on a pixel-by-pixel basis will mitigate some of these effects.

The best-linear parameters of the different power laws found for the CO–F770W and CO–F1130W relations, in particular the normalization parameter, suggest that mid-IR is brighter than the CO in M 82. On the one hand, this could be due to effects directly related to the physical conditions of dust grains. PAH emission depends on the radiation field strength (e.g. Draine & Li 2007; Leroy et al. 2023; Chown et al. 2024), which has been shown to be significantly higher in M 82 than the star-forming spiral disks in C24. In addition, using [CII] 158 μm line data taken with the Stratospheric Observatory for Infrared Astronomy (SOFIA), Levy et al. (2023) find that the UV field varies up to ∼2–3 orders of magnitude from the outflow to the disk of M 82. On the other hand, the results could be indicating that the PAH features arise from (warm) gas phases not completely traced by CO(1–0), perhaps due to photoionization produced by the outflow that may be suppressing the CO emission (e.g., Shimizu et al. 2019) or turbulent gas producing changes in the strength of CO emission (e.g., Teng et al. 2022). Low values of CO-to-H2 conversion factor (αCO) are also expected for M 82 and differ to those for disks of PHANGS galaxies (Leroy et al. 2022). In this sense, several gas phases have already been detected in outflows aside of molecular gas (e.g. Martini et al. 2018; Levy et al. 2021), including PAH mixed with atomic gas (e.g. Yun et al. 1994). Finally, another interpretation is that the offset in the CO-PAH relationships shown in Figure 2 could be just a consequence of the Kennicutt–Schmidt relation (Schmidt 1959; Kennicutt 1998), where the molecular gas content is traced by the ICO (y-axes) and star-formation rate on the (x-axes), since PAHs can trace the SFR (e.g., Shipley et al. 2016). In this context, the offset may be simply reflecting M 82’s starburst nature, where the star-formation efficiency is higher rather than a suppression of CO.

4.2. CO versus F770W/F1130W

PAH-band ratios are commonly used to estimate the strength of nearby radiation fields since they have been shown to be very sensitive to the ionization level and size of PAH grains. Particularly, the 7.7 μm/11.3 μm PAH ratio has been shown to be sensitive to the PAH state, with a lower value reflecting more neutral PAHs (e.g., Tielens 2008; Draine 2011; Draine et al. 2021; Chastenet et al. 2023; Egorov et al. 2023). Using the 3.3 μm PAH feature, Bolatto et al. (2024) detect the presence of small PAH grains in the innermost regions of M 82 very close to the ionizing starburst, which extend toward larger and more neutral grains at larger distances. If the conditions of the outflow energetics are enough to ionize and eventually destroy a significant fraction of PAH grains linked to H2 clouds: can the CO–IMIRI770/IMIRI1130 relation be used as a predictor of the physical changes of the cold gas (e.g., the traceability of H2 through CO)?

The top panel of Figure 4 presents the relation between the CO and the IMIRI770/IMIRI1130 ratio (hereafter the CO–F770W/F1130W relation) for the 306 clouds identified in this work. Except for clouds in the streamer east and the south outflow, the rest of the groups have moderate correlations for the CO–F770W/F1130W relation (Pearson rp = 0.32 − 0.62; see column 4 in Table 1), showing a constant behavior of the IMIRI770/IMIRI1130 ratio with increasing ICO(1 − 0). We also note that there is a point at log(IMIRI770/IMIRI1130)∼0.56 where the relation (with a significant fraction of clouds on the disk) shows an increasing in the scatter. This is even more clear when we take data bins in the x-axis (unfilled-black pentagons in Fig. 4). Similarly to the top panel, the middle panel of Figure 4 shows the CO-F770W/F1130W relationship colorcoded by the size of the CO clouds. We note signs of two groups of clouds split by their sizes, where clouds with sizes larger than 100 pc (i.e., typical sizes of giant molecular clouds, GMC) tend to have log(ICO(1 − 0))≳1. Beyond interpretations of the physical processes behind the apparent increase in the scatter of the relationship, we note that this effect may be induced by the low number of clouds identified with ratios greater than log(IMIRI770/IMIRI1130)≳0.56.

thumbnail Fig. 4.

ICO(1 − 0) versus IMIRI770/IMIRI1130 ratio. Top: symbols follow the same convention as in Fig. 3. Middle: symbols are colorcoded by the size of the CO cloud derived in Section 3.2; while circles are CO clouds with size larger than 100 pc, stars correspond to those smaller than this limit. The legend on the bottom panel contains the values of the Pearson rp for the two groups. The vertical dashed-dotted red lines mark the point log(IMIRI770/IMIRI1130)≈0.56 where the scatter of the CO–F770W/F1130W relationship increases significantly (see Section 4.2). Bottom: same as middle panel but colorcoded by the projected distance of the CO cloud respect to the optical center.

The three panels in Figure 4 show some physical differences between distinct locations of the CO clouds, which seem to depend not only on their location within the galaxy, but also on their sizes and the projected distance to the center of M 82. Similarly to the right panels of Figure 3, the middle panel of Figure 4 shows that when splitting our sample into two size groups, we find a stronger correlation in the CO–F770W/F1130W relation in clouds larger than 100 pc (Pearson rp = 0.51) compared to that for those with shorter diameters (Pearson rp ∼ 0.01). Analogously to Section 4.1, these results seem to reflect that larger CO clouds are both brighter and cover a broader range of ICO values (not limited by sensitivity) than smaller clouds. Figure 4 also includes the projected distances of the CO clouds respect to the optical center (bottom panel). Interestingly, we note that small clouds (sizes ≲50 pc) from the streamer east and outflow south (i.e., those with low rp in top panel of Fig. 4), and at distances larger than 1 kpc from the optical center contribute significantly to the scatter of the CO–F770W/F1130W relationship at log(IMIRI770/IMIRI1130) ∼ 0.56. These results seem to confirm the noticeable variability in the heating mechanism of PAHs in the different regions of M 82, particularly in the smallest CO clouds exposed to hard radiation fields. Such effects could ultimately produce a mild decoupling between CO and H2 emission (and possibly between H2 and PAHs; e.g., Berné et al. 2009; Hill & Zakamska 2014).

Summarizing, our results suggest that the relations linking the CO to both PAH emission features and their ratios vary significantly in such extreme environments. Changes in the heating mechanism of PAHs and sizes of the CO clouds, or an increase in the H2 content not traced by CO are perhaps plausible scenarios in this context. Upcoming studies will thus focus on a careful analysis to address some of these topics more in deep, either verifying the dependency of the CO-7.7 μm PAH feature relation on vertical distance (e.g., Lopez et al. in prep.), including the 3.3 μm PAH feature and strength of the surrounding radiation fields (e.g., Arens et al., in prep.), and a complete inspection of the F1130M-to-F770M ratio in the entire area covered by the MIRI data for M 82 (Cronin et al., in prep.).

5. Summary and conclusions

We present a study of the relation between new, very high-resolution JWST-MIRI data and archival CO(1–0) NOEMA observations of the local starburst galaxy M 82. We focus on the dependence of the integrated CO(1–0) line intensity on the MIRI F770W and F1130W band emission. To do so, we identify 306 CO clouds in the inner 2 kpc from the central starburst and we classify them into regions as described in Krieger et al. (2021) and by their sizes. Our main conclusions are:

  1. We compute the relations between the integrated CO(1–0) line intensities and the F770W and F1130W band emission. On average, we find that except for the smaller CO clouds and those from the streamer east, most of the cloud groups show a moderate or strong correlation. We compare the best-linear parameters found for nearby spirals (e.g., m ∼ 0.9–1.0; PHANGS) with those we obtain for M 82. On average, we find both lower global slopes and amplitudes for the CO–F770W and CO–F1130W relations (after binning the data by the F770W and F1130W values in the x-axes) than what is expected for local main-sequence galaxies.

  2. We derive the relation between CO(1–0) and the F770W-to-F1130W ratio. In general, we see a moderate correlation of the CO–F770W/F1130W relation for most of the cloud groups analyzed in this work (Pearson rp ∼ 0.3 − 0.6). Although we note a constant behavior of the CO intensities as a function of F770W/F1130W, the scatter of this relation increases significantly above log(IMIRI770/IMIRI1130)∼0.56. Since most of the CO emission is concentrated in large clouds from the disk region, this point seems to be linked with significant changes in the physical conditions of the CO. Based on the widespread use of I7.7 μm/I11.3 μm as a prescription of the ionization and/or destruction level of PAHs, our results support a scenario with more UV photons available to either break down or suppress CO emission in small clouds located outside the central region of M 82, due to strong radiation fields.

The best-linear parameters of the different power laws found for the CO–PAH relations presented in this work suggest that mid-IR is brighter than the CO in M 82 compared to local star-forming galaxies and that the heating mechanisms of PAH change dramatically depending on the region analyzed. However, when compared to the native angular resolution of the MIRI data, the larger beamsize of the CO NOEMA observations could be introducing biases in our interpretations. This effect could impact the integrated intensities that we estimate for the CO as a result of potential dilution of the CO emission in a larger area than that constrained by PAH’s regions. Moreover, the relations included in this work could be severely limited by the CO data sensitivity; consequently, we could be missing fainter clouds in our CO clouds identification routine. This could drive biases in the derivations of the power-law parameters, e.g. producing artificially shallower slopes for the global relations compared to those for local main-sequence galaxies.

Future studies should therefore focus on both high-resolution and sensitivity CO rotational transitions or even direct H2 ro-vibrational line data. In addition, analyses of excitation conditions derived from multi-J CO transitions could also help to address some of these questions. These will be crucial to verify in detail the actual connection between CO, H2, and PAH features in the extreme environment of M 82 at physical scales comparable to those of even smaller molecular clouds (≲20 pc).


Acknowledgments

The authors acknowledge the comments from the anonymous referee, which helped to improve this manuscript. V. V. acknowledges support from the ALMA-ANID Postdoctoral Fellowship under the award ASTRO21-0062 and from the ANID BASAL project FB210003. I.D.L. acknowledges funding from the Belgian Science Policy Office (BELSPO) through the PRODEX project “JWST/MIRI Science exploitation” (C4000142239) and funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program DustOrigin (ERC-2019-StG-851622). R.H.-C. thanks the Max Planck Society for support under the Partner Group project “The Baryon Cycle in Galaxies” between the Max Planck for Extraterrestrial Physics and the Universidad de Concepción. R.H-C. also gratefully acknowledge financial support from ANID – MILENIO – NCN2024_112 and ANID BASAL FB210003. M.R. acknowledges support from project PID2023-150178NB-I00 (and PID2020-114414GB-I00*), financed by MICIU/AEI/10.13039/501100011033, and by FEDER, UE. L.A.L. and S.L. were supported by NASA’s Astrophysics Data Analysis Pro- gram under grant No. 80NSSC22K0496SL, and L.A.L. also acknowledges support through the Heising-Simons Foundation grant 2022-3533. A.D.B. and S.A.C. acknowledge support from grant JWST-GO-01701.001-A by STScI. This work is based on observations made with the NASA/ESA/CSA James Webb Space Telescope. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with program JWST-GO- 01701. Support for program JWST-GO-01701 is provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127. This work is also based on observations carried out under project Nos. w18by and 107-19 with the IRAM NOEMA Interferometer and the IRAM 30 m telescope, respectively. IRAM is supported by INSU/CNRS (France), MPG (Germany), and IGN (Spain). This research has made use of NASA’s Astrophysics Data System Bibliographic Services. Software: Astropy (Astropy Collaboration 2018), MatPlotLib (Hunter 2007), NumPy (Harris et al. 2020), SciPy (Virtanen et al. 2020), seaborn (Waskom 2021), Scikit-learn (Pedregosa et al. 2011).

References

  1. Alonso-Herrero, A., Pereira-Santaella, M., Rigopoulou, D., et al. 2020, A&A, 639, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Astropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123 [Google Scholar]
  3. Beirão, P., Brandl, B. R., Appleton, P. N., et al. 2008, ApJ, 676, 304 [CrossRef] [Google Scholar]
  4. Berné, O., Fuente, A., Goicoechea, J. R., et al. 2009, ApJ, 706, L160 [CrossRef] [Google Scholar]
  5. Bezanson, R., Spilker, J. S., Suess, K. A., et al. 2022, ApJ, 925, 153 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bolatto, A. D., Levy, R. C., Tarantino, E., et al. 2024, ApJ, 967, 63 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bundy, K., Bershady, M. A., Law, D. R., et al. 2015, ApJ, 798, 7 [Google Scholar]
  8. Calzetti, D. 2013, in Secular Evolution of Galaxies, eds. J. Falcón-Barroso, & J. H. Knapen, 419 [Google Scholar]
  9. Chastenet, J., Sutter, J., Sandstrom, K., et al. 2023, ApJ, 944, L11 [NASA ADS] [CrossRef] [Google Scholar]
  10. Chastenet, J., De Looze, I., Relaño, M., et al. 2024, A&A, 690, A348 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Chown, R., Li, C., Parker, L., et al. 2021, MNRAS, 500, 1261 [Google Scholar]
  12. Chown, R., Leroy, A. K., Sandstrom, K., et al. 2024, ApJ, submitted [arXiv:2410.05397] [Google Scholar]
  13. Cortzen, I., Garrett, J., Magdis, G., et al. 2019, MNRAS, 482, 1618 [NASA ADS] [CrossRef] [Google Scholar]
  14. Davies, R., Shimizu, T., Pereira-Santaella, M., et al. 2024, A&A, 689, A263 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. den Brok, J. S., Chatzigiannakis, D., Bigiel, F., et al. 2021, MNRAS, 504, 3221 [NASA ADS] [CrossRef] [Google Scholar]
  16. Donnelly, G. P., Lai, T. S. Y., Armus, L., et al. 2025, ApJ, submitted [arXiv:2501.19397] [Google Scholar]
  17. Draine, B. T. 2011, ApJ, 732, 100 [NASA ADS] [CrossRef] [Google Scholar]
  18. Draine, B. T., & Li, A. 2007, ApJ, 657, 810 [CrossRef] [Google Scholar]
  19. Draine, B. T., Li, A., Hensley, B. S., et al. 2021, ApJ, 917, 3 [CrossRef] [Google Scholar]
  20. Egorov, O. V., Kreckel, K., Sandstrom, K. M., et al. 2023, ApJ, 944, L16 [NASA ADS] [CrossRef] [Google Scholar]
  21. Egusa, F., Gao, Y., Morokuma-Matsui, K., Liu, G., & Maeda, F. 2022, ApJ, 935, 64 [NASA ADS] [CrossRef] [Google Scholar]
  22. Fisher, D. B., Bolatto, A. D., Chisholm, J., et al. 2024, MNRAS, submitted [arXiv:2405.03686] [Google Scholar]
  23. Förster Schreiber, N. M., Sauvage, M., Charmandaris, V., et al. 2003, A&A, 399, 833 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Freedman, W. L., Hughes, S. M., Madore, B. F., et al. 1994, ApJ, 427, 628 [NASA ADS] [CrossRef] [Google Scholar]
  25. Gao, Y., Xiao, T., Li, C., et al. 2019, ApJ, 887, 172 [Google Scholar]
  26. Gordon, K. D., Bohlin, R., Sloan, G. C., et al. 2022, AJ, 163, 267 [NASA ADS] [CrossRef] [Google Scholar]
  27. Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357 [NASA ADS] [CrossRef] [Google Scholar]
  28. Hill, M. J., & Zakamska, N. L. 2014, MNRAS, 439, 2701 [Google Scholar]
  29. Houck, J. R., Roellig, T. L., van Cleve, J., et al. 2004, ApJS, 154, 18 [NASA ADS] [CrossRef] [Google Scholar]
  30. Hunter, J. D. 2007, Computing in Science& Engineering, 9, 90 [NASA ADS] [CrossRef] [Google Scholar]
  31. Indebetouw, R., Wong, T., Chen, C. H. R., et al. 2020, ApJ, 888, 56 [NASA ADS] [CrossRef] [Google Scholar]
  32. Keenan, R. P., Marrone, D. P., & Keating, G. K. 2024, ApJ, 979, 16 [Google Scholar]
  33. Kennicutt, R. C., Jr 1998, ApJ, 498, 541 [Google Scholar]
  34. Kennicutt, R. C., & Evans, N. J. 2012, ARA&A, 50, 531 [NASA ADS] [CrossRef] [Google Scholar]
  35. Krieger, N., Walter, F., Bolatto, A. D., et al. 2021, ApJ, 915, L3 [NASA ADS] [CrossRef] [Google Scholar]
  36. Leroy, A. K., Walter, F., Martini, P., et al. 2015, ApJ, 814, 83 [Google Scholar]
  37. Leroy, A. K., Hughes, A., Liu, D., et al. 2021a, ApJS, 255, 19 [NASA ADS] [CrossRef] [Google Scholar]
  38. Leroy, A. K., Schinnerer, E., Hughes, A., et al. 2021b, ApJS, 257, 43 [NASA ADS] [CrossRef] [Google Scholar]
  39. Leroy, A. K., Rosolowsky, E., Usero, A., et al. 2022, ApJ, 927, 149 [NASA ADS] [CrossRef] [Google Scholar]
  40. Leroy, A. K., Bolatto, A. D., Sandstrom, K., et al. 2023, ApJ, 944, L10 [NASA ADS] [CrossRef] [Google Scholar]
  41. Levy, R. C., Bolatto, A. D., Leroy, A. K., et al. 2021, ApJ, 912, 4 [NASA ADS] [CrossRef] [Google Scholar]
  42. Levy, R. C., Bolatto, A. D., Tarantino, E., et al. 2023, ApJ, 958, 109 [NASA ADS] [CrossRef] [Google Scholar]
  43. Levy, R. C., Bolatto, A. D., Mayya, D., et al. 2024, ApJ, 973, L55 [NASA ADS] [CrossRef] [Google Scholar]
  44. Marshall, J. A., Herter, T. L., Armus, L., et al. 2007, ApJ, 670, 129 [Google Scholar]
  45. Martini, P., Leroy, A. K., Mangum, J. G., et al. 2018, ApJ, 856, 61 [Google Scholar]
  46. Mayya, Y. D., Romano, R., Rodríguez-Merino, L. H., et al. 2008, ApJ, 679, 404 [NASA ADS] [CrossRef] [Google Scholar]
  47. McCrady, N., Gilbert, A. M., & Graham, J. R. 2003, ApJ, 596, 240 [NASA ADS] [CrossRef] [Google Scholar]
  48. McKeith, C. D., Castles, J., Greve, A., & Downes, D. 1993, A&A, 272, 98 [NASA ADS] [Google Scholar]
  49. Oppenheimer, B. D., Davé, R., Kereš, D., et al. 2010, MNRAS, 406, 2325 [NASA ADS] [CrossRef] [Google Scholar]
  50. Page, M. J., Symeonidis, M., Vieira, J. D., et al. 2012, Nature, 485, 213 [NASA ADS] [CrossRef] [Google Scholar]
  51. Papadopoulos, P. P., Bisbas, T. G., & Zhang, Z.-Y. 2018, MNRAS, 478, 1716 [NASA ADS] [CrossRef] [Google Scholar]
  52. Pedregosa, F., Varoquaux, G., Gramfort, A., et al. 2011, Journal of Machine Learning Research, 12, 2825 [Google Scholar]
  53. Pope, A., Wagg, J., Frayer, D., et al. 2013, ApJ, 772, 92 [NASA ADS] [CrossRef] [Google Scholar]
  54. Recchi, S., & Hensler, G. 2013, A&A, 551, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Regan, M. W., Thornley, M. D., Vogel, S. N., et al. 2006, ApJ, 652, 1112 [NASA ADS] [CrossRef] [Google Scholar]
  56. Rosolowsky, E. W., Pineda, J. E., Kauffmann, J., & Goodman, A. A. 2008, ApJ, 679, 1338 [Google Scholar]
  57. Schmidt, M. 1959, ApJ, 129, 243 [NASA ADS] [CrossRef] [Google Scholar]
  58. Shimizu, T. T., Davies, R. I., Lutz, D., et al. 2019, MNRAS, 490, 5860 [Google Scholar]
  59. Shipley, H. V., Papovich, C., Rieke, G. H., Brown, M. J. I., & Moustakas, J. 2016, ApJ, 818, 60 [NASA ADS] [CrossRef] [Google Scholar]
  60. Shivaei, I., & Boogaard, L. A. 2024, A&A, 691, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Shivaei, I., Alberts, S., Florian, M., et al. 2024, A&A, 690, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Sidorin, V. 2017, Astrophysics Source Code Library [record ascl:1704.006] [Google Scholar]
  63. Smith, L. J., Westmoquette, M. S., Gallagher, J. S., et al. 2006, MNRAS, 370, 513 [NASA ADS] [CrossRef] [Google Scholar]
  64. Smith, J. D. T., Draine, B. T., Dale, D. A., et al. 2007, ApJ, 656, 770 [Google Scholar]
  65. Spilker, J. S., Phadke, K. A., Aravena, M., et al. 2023, Nature, 618, 708 [NASA ADS] [CrossRef] [Google Scholar]
  66. Suess, K. A., Bezanson, R., Spilker, J. S., et al. 2017, ApJ, 846, L14 [NASA ADS] [CrossRef] [Google Scholar]
  67. Sutter, J., Sandstrom, K., Chastenet, J., et al. 2024, ApJ, 971, 178 [NASA ADS] [CrossRef] [Google Scholar]
  68. Swinbank, A. M., Papadopoulos, P. P., Cox, P., et al. 2011, ApJ, 742, 11 [NASA ADS] [CrossRef] [Google Scholar]
  69. Teng, Y.-H., Sandstrom, K. M., Sun, J., et al. 2022, ApJ, 925, 72 [NASA ADS] [CrossRef] [Google Scholar]
  70. Tielens, A. G. G. M. 2008, ARA&A, 46, 289 [NASA ADS] [CrossRef] [Google Scholar]
  71. Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261 [CrossRef] [Google Scholar]
  72. Walter, F., Weiss, A., & Scoville, N. 2002, ApJ, 580, L21 [Google Scholar]
  73. Waskom, M. L. 2021, Journal of Open Source Software, 6, 3021 [CrossRef] [Google Scholar]
  74. Weiß, A., Walter, F., & Scoville, N. Z. 2005, A&A, 438, 533 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Williams, J. P., de Geus, E. J., & Blitz, L. 1994, ApJ, 428, 693 [Google Scholar]
  76. Wünsch, R., Jáchym, P., Sidorin, V., et al. 2012, A&A, 539, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Yajima, Y., Sorai, K., Miyamoto, Y., et al. 2021, PASJ, 73, 257 [NASA ADS] [CrossRef] [Google Scholar]
  78. Yun, M. S., Ho, P. T. P., & Lo, K. Y. 1994, Nature, 372, 530 [NASA ADS] [CrossRef] [Google Scholar]
  79. Zschaechner, L. K., Bolatto, A. D., Walter, F., et al. 2018, ApJ, 867, 111 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Continuum subtraction

Depending on the region of M 82, mid-IR data may be susceptible to non-negligible contribution of continuum emission from dust and star (especially for the F770W map). Using the ISOCAM instrument on board the Infrared Space Observatory, Förster Schreiber et al. (2003) showed that the continuum at λ = 7.7 μm has been shown to represent up to ∼20% of the total emission in the core region in M 82. We address this by implementing a first-order continuum subtraction correction to the MIRI-F770W and F1130W filter intensities. We perform the continuum subtraction using the spectra from regions A and B (yellow boxes), and the disk (black rectangles) included in the middle and right panels of Figure 1. For the disk, we first fit the continuum emission given by the MIRI medium resolution spectroscopy (MRS) data available for a portion of the disk region of M 82, and we cover the wavelength range from λ ≈ [5.0, 18] μm. Next, we extract the continuum emission at the PAH-free parts of the spectra from all the spaxels available at 5.0, 15.0, and 16.0 μm, thus avoiding the 10 μm region (since its susceptibility to silicate absorption up to 2.5 mag; see Beirão et al. 2008). Then, we fit them using the model ImidIR − Cont = A × λ + B (see the top panel of Fig.A.1). Finally, we integrate both the best-fit continuum function and the spectra intensities along the wavelength range covered by the F770W and F1130W filters4 ([6.581,8.687] μm and [10.953,11.667] μm, respectively), and compute the ratio between the resulting continuum and total intensity. We follow this procedure to also perform the continuum subtraction for the other regions using the combined spectrum from areas A and B (see Fig. 1), originally included in Beirão et al. (2008). The data consists of combined short-low and long-low spectra taken with the IRS of the Spitzer Space Telescope, covering the spectral range between 5 and 38 μm (see bottom panel of Fig. A.1).

thumbnail Fig. A.1.

Top: MIRI-MRS for a spaxel within the black region in middle and right panels of Fig. 1. While green crosses represent the regions where continuum emission is extracted, the red-dashed and gray-dashed lines correspond to the best-fit function for such points and the continuum subtraction applied to the spectrum, respectively. Bottom:Spitzer IRS spectrum taken from regions A and B in Fig. 1 and included originally in Beirão et al. (2008). Conventions are as in top panel.

Figure A.2 shows the continuum-to-total emission fraction at 7.7 μm (top) and at 11.3 μm (bottom) as a function of IMIRI770/IMIRI1130 (before continuum subtraction) in the spaxels within the region covered by the MIRI-MRS. For both the F770W and F1130W filter intensities, we apply a first-order correction approximating the dependency of the continuum-contribution fraction on the F770M/F1130W ratio (before subtraction) by the function f(IMIRI770/IMIRI1130) = α − β/(1 + eγ(IMIRI770/IMIRI1130 − ρ)). While for the [7.7CONTμm]/[7.7PAHμm] we find the best-fit parameters α = 0.16, β = 0.06, γ = 90.1, and ρ = 0.59, for the [11.3CONTμm]/[11.3PAHμm] we obtain α = 0.19, β = 0.1, γ = 147, and ρ = 0.55. Finally, we use these best-fit functions to obtain the continuum contribution (i.e., the continuum-to-total emission fraction) in the disk area based on the IMIRI770/IMIRI1130 values.

thumbnail Fig. A.2.

Continuum-to-total emission fraction at 7.7 μm (top) and at 11.3 μm (bottom) versus IMIRI770/IMIRI1130. While empty pentagons correspond to the data binned in the x-axis, the black-dashed lines are the best-fit functions to such binned data. Although with some scatter, at first order the distributions of the points are well described by the model α − β/(1 + eγ(IMIRI770/IMIRI1130 − ρ)). The top and right histograms in both panels correspond to the distribution of the values for the respective axes, and the red-dashed lines mark the median of such distributions.

Appendix B: Cloud identification

We use QUICKCLUMP (Sidorin 2017) to identify CO clouds in the CO(1-0) NOEMA moment 0 map. In a nutshell, QUICKCLUMP was developed to identify clouds standing on the hierarchical structure of the molecular gas in galaxies. QUICKCLUMP creates dendrograms as a natural byproduct of the clump-finding process. The algorithm is similar to the well-known CLUMPFIND (Williams et al. 1994); its results, however, are less dependent on technical parameters (such as the temperature difference between contours). The algorithm is primarily controlled by three parameters: NPXMIN (minimum number of pixels of a clump), DTLEAF (minimum difference between the peak temperature of the cloud and the temperature at which it connects to other clouds), and TCUTOFF (minimum ICO(1 − 0) considered for assignment to clouds).

With this in mind, we run QUICKCLUMP to find clouds at least 3σ CO detected and with a minimal diameter of 6 pixels of the CO(J = 1 − 0) map (i.e., equivalent to the mean between the minor and major axes of the beamsize θmean = 1.8″≈31 pc∼5.2 pixels). We discard clouds identified in the CO map outside the coverage of the F770W and F1130W maps, validating 306 CO clouds.

All Tables

Table 1.

Correlation coefficients for the relations between CO and F770W, F1130W, and F770W/F1130W.

All Figures

thumbnail Fig. 1.

From left to right: M 82 NOEMA CO(J = 1 − 0), MIRI-F770W, and MIRI-F1130W images in cutouts of 5′×4.2′. The left panel also contains the molecular clouds identified in Section 3.2 (white contours), the centroid of the CO(J = 1 − 0) and the different regions identified in the starburst as described in Krieger et al. (2021). The inset in the left panel correspond to the NOEMA beamsize at ν = 115 GHz. The green crosses in the middle and right panels are indicate the optical center of M 82. The black contours in the middle and right panels correspond to the regions covered by the MIRI-MRS and the yellow boxes are the regions with Spitzer IRS data from Beirão et al. (2008); both are used in this work to apply a first-order continuum subtraction (see Section 3.1).

In the text
thumbnail Fig. 2.

Distribution of sizes of the 306 CO clouds we identify in the NOEMA map ( θ mean = 1 . 8 31 $ \theta_{\mathrm{mean}}={1}{{\overset{\prime\prime}{.}}}8 \approx 31 $ pc) using QUICKCLUMP as described in Section 3.2. The red-dashed line marks the median size for the cloud sizes. The inset includes the distributions of clouds respect to the five regions we use to characterize M 82 (see Section 3.2).

In the text
thumbnail Fig. 3.

ICO(1 − 0) versus IMIRI770 (top), and IMIRI1130 (bottom). Points in left panels are colored in the five regions we adopt to classify the clouds in M 82 (see Section 3.2). We also include the best-fit linear relations we obtain for the whole sample after binning the points in the x-axis and performing an ordinary least square (OLS) linear fitting using the model y = mx + b (black-dashed lines), and the Pearson’s r-value for the different regions in the right-bottom corners. Error bars correspond to the scatter of the points in the y-axis for the binned data. The left panels show the same relations that those in right ones, but colorcoded by the size of clouds, where symbols are split in clouds larger (circles) and shorter (stars) than 100 pc.

In the text
thumbnail Fig. 4.

ICO(1 − 0) versus IMIRI770/IMIRI1130 ratio. Top: symbols follow the same convention as in Fig. 3. Middle: symbols are colorcoded by the size of the CO cloud derived in Section 3.2; while circles are CO clouds with size larger than 100 pc, stars correspond to those smaller than this limit. The legend on the bottom panel contains the values of the Pearson rp for the two groups. The vertical dashed-dotted red lines mark the point log(IMIRI770/IMIRI1130)≈0.56 where the scatter of the CO–F770W/F1130W relationship increases significantly (see Section 4.2). Bottom: same as middle panel but colorcoded by the projected distance of the CO cloud respect to the optical center.

In the text
thumbnail Fig. A.1.

Top: MIRI-MRS for a spaxel within the black region in middle and right panels of Fig. 1. While green crosses represent the regions where continuum emission is extracted, the red-dashed and gray-dashed lines correspond to the best-fit function for such points and the continuum subtraction applied to the spectrum, respectively. Bottom:Spitzer IRS spectrum taken from regions A and B in Fig. 1 and included originally in Beirão et al. (2008). Conventions are as in top panel.

In the text
thumbnail Fig. A.2.

Continuum-to-total emission fraction at 7.7 μm (top) and at 11.3 μm (bottom) versus IMIRI770/IMIRI1130. While empty pentagons correspond to the data binned in the x-axis, the black-dashed lines are the best-fit functions to such binned data. Although with some scatter, at first order the distributions of the points are well described by the model α − β/(1 + eγ(IMIRI770/IMIRI1130 − ρ)). The top and right histograms in both panels correspond to the distribution of the values for the respective axes, and the red-dashed lines mark the median of such distributions.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.