# Compact and high excitation molecular clumps in the extended ultraviolet disk of M83

Jin Koda, Francoise Combes, Mónica Rubio, Morten Andersen, Frank Bigiel, Armando Gil de Paz, Junais Snc (sans Nom Connu/ No Family Name), Amanda Lee, Jennifer Donovan Meyer, Kana Morokuma-Matsui, et al.

► To cite this version:

Jin Koda, Francoise Combes, Mónica Rubio, Morten Andersen, Frank Bigiel, et al. Compact and high excitation molecular clumps in the extended ultraviolet disk of M83. Astronomy & Astrophysics
- A&A, 2024, 691, pp.A197. 10.1051/0004-6361/202451506. hal-04996957

HAL Id: hal-04996957

https://college-de-france.hal.science/hal-04996957v1

Submitted on 20 Mar 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Compact and high excitation molecular clumps in the extended ultraviolet disk of M83

Jin Koda<fnref n="1" />,<fnref n="2" />,*, Françoise Combes<fnref n="3" />,*, Monica Rubio<fnref n="4" />,*, Morten Andersen<fnref n="5" />,*, Frank Bigiel<fnref n="6" />,*, Armando Gil de Paz<fnref n="7" />,<fnref n="8" />,*, Junais<fnref n="9" />,*, Amanda M Lee<fnref n="10" />,*, Jennifer Donovan Meyer<fnref n="11" />,*, Kana Morokuma-Matsui<fnref n="12" />,*, Masafumi Yagi<fnref n="13" />,*, and Annie Zavagno<fnref n="14" />

Received 15 July 2024 / Accepted 23 September 2024

## ABSTRACT

*Context.* The extended ultraviolet (XUV) disks of nearby galaxies show ongoing massive-star formation, but their parental molecular clouds remain mostly undetected despite searches in CO(1–0) and CO(2–1). The recent detection of 23 clouds in the higher excitation transition CO(3–2) within the XUV disk of M83 thus requires an explanation.

*Aims.* We test the hypothesis introduced to explain the non-detections and recent detection simultaneously: The clouds in XUV disks have a clump-envelope structure similar to those in Galactic star-forming clouds, having dense star-forming clumps (or concentrations of multiple clumps) at their centers, which predominantly contribute to the CO(3–2) emission and are surrounded by less dense envelopes, where CO molecules are photo-dissociated due to the low-metallicity environment there.

*Methods.* We utilize new high-resolution ALMA CO(3–2) observations of a subset (11) of the 23 clouds in the XUV disk of M83.

*Results.* We confirm the compactness of the CO(3–2)-emitting dense clumps (or their concentrations), finding clump diameters below the spatial resolution of 6–9 pc. This is similar to the size of the dense gas region in the Orion A molecular cloud, a local star-forming cloud with massive-star formation.

*Conclusions.* The dense star-forming clumps are common between normal and XUV disks. This may also indicate that once the cloud structure is set, the process of star formation is governed by the cloud’s internal physics rather than by external triggers. This simple model explains the current observations of clouds with ongoing massive-star formation, although it may require some adjustment, for example including the effect of cloud evolution, to describe star formation in molecular clouds more generally.

**Key words.** ISM: clouds – ISM: molecules – galaxies: ISM – galaxies: star formation

## 1. Introduction

The Galaxy Evolution Explorer (GALEX) satellite has enabled the discovery of a surprisingly large number of massive-star formation (SF) sites beyond the optical radii ($R_{25}$) of galactic disks (often ~2–4 $R_{25}$; Gil de Paz et al. 2005; Thilker et al. 2005; Bigiel et al. 2010). Numerous bright UV sources (OB stars) are present beyond the edges of the optical disks and extend into the low-density HI gas disks (see the example of M83 in Bigiel et al. 2010; Koda et al. 2012, 2022; Eibensteiner et al. 2023; Rautio et al. 2024). Such extended ultraviolet (XUV)

disks are found in about one-third of local disk galaxies (Thilker et al. 2007). Thus, this mode of SF is fairly common despite the extreme conditions. It is important to study the parental sites of this SF activity: the molecular gas clouds.

Efforts have been made to detect CO emission (molecular gas) in XUV disks but have rarely succeeded (see Watson & Koda 2017 for a review, including unpublished non-detections). Such studies utilized CO($J=1-0$) or CO($J=2-1$) emission, widely used tracers of the bulk molecular gas in normal molecular clouds, as these lines are excited easily at the average density ($n_{H_2} \sim 300 \text{ cm}^{-3}$) and temperature ($T_k \sim 10 \text{ K}$) of the bulk gas. Until recently, CO was detected in the outer disk environment of only four galaxies (Braine & Herpin 2004;

<footnote n="1">¹ Observatoire de Paris, LERMA, PSL Univ., CNRS, Sorbonne Univ., Paris, France</footnote>

<footnote n="2">² Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794-3800, USA</footnote>

<footnote n="3">³ Observatoire de Paris, LERMA, Collège de France, PSL Univ., CNRS, Sorbonne Univ., Paris, France</footnote>

<footnote n="4">⁴ Departamento de Astronomia, Universidad de Chile, Casilla 36-D, Santiago, Chile</footnote>

<footnote n="5">⁵ European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany</footnote>

<footnote n="6">⁶ Argelander-Institut für Astronomie, Universität Bonn, Auf dem Hügel 71, D-53121 Bonn, Germany</footnote>

<footnote n="7">⁷ Departamento de Física de la Tierra y Astrofísica, Facultad de CC. Físicas, Universidad Complutense de Madrid, 28040 Madrid, Spain</footnote>

<footnote n="8">⁸ Instituto de Física de Partículas y del Cosmos IPARCOS, Facultad de CC. Físicas, Universidad Complutense de Madrid, 28040 Madrid, Spain</footnote>

<footnote n="9">⁹ National Centre for Nuclear Research, Pasteura 7, 02-093 Warsaw, Poland</footnote>

<footnote n="10">¹⁰ Department of Astronomy, University of Massachusetts Amherst, 710 North Pleasant Street, Amherst, MA 01003, USA</footnote>

<footnote n="11">¹¹ National Radio Astronomy Observatory (NRAO), 520 Edgemont Road, Charlottesville, VA 22903, USA</footnote>

<footnote n="12">¹² Institute of Astronomy, Graduate School of Science, The University of Tokyo, 2-21-1 Osawa, Mitaka, Tokyo 181-0015, Japan</footnote>

<footnote n="13">¹³ National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan</footnote>

<footnote n="14">¹⁴ Aix Marseille Univ., CNRS, CNES, Laboratoire d'Astrophysique de Marseille, Marseille, France</footnote>

<footnote>* Corresponding author; jin.koda@stonybrook.edu</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Clump-envelope structure of a molecular cloud in an environment of (a) high metallicity and (b) low metallicity, suggested by Koda et al. (2022) based on previous studies (e.g., Maloney & Black 1988; Wolfire et al. 2010). In the low-metallicity environment (b), CO molecules are selectively photo-dissociated in the envelope (orange) by the external UV radiation field, producing a thick CO-deficient H₂ layer (gray), which substantially reduces the CO(2–1) and CO(1–0) emission. The dense star-forming clumps (brown and their immediate surroundings) reside at the hearts of the clouds; the CO molecules are protected there, and the clumps can remain bright in CO(3–2).</figcaption></figure>

Braine et al. 2007, 2010; Dessauges-Zavadsky et al. 2014), and only one of the four is identified as a single molecular cloud (Braine et al. 2012). Even a sensitive search in CO(2–1) with the Atacama Large Millimeter/submillimeter Array (ALMA) over a $2 \times 4$ kpc² region in the XUV disk of M83 resulted in a non-detection (Bicalho et al. 2019). The rarity of detections in CO(1–0) and CO(2–1) has been at odds with the abundance of UV sources (indicating massive SF) across these XUV disks.

A breakthrough has come recently. ALMA revisited the XUV disk of M83, a 1 kpc² region (Koda et al. 2022) within the larger region of the non-detection in CO(2–1) (Bicalho et al. 2019). It detected 23 molecular clouds in CO(3–2) emission. While these clouds are spatially unresolved at 1" resolution (20 pc), the detection led to a hypothesis, presented in Koda et al. (2022), regarding the parental clouds of the outer disk massive SF (see Sect. 2). This work aims to test this hypothesis with higher-resolution CO(3–2) observations, in conjunction with a consideration of the previous CO(2–1) non-detection.

This hypothesis assumes that the clouds have a warm, dense gas concentration surrounded by an envelope of more diffuse molecular gas (Sect. 2). We refer to this dense gas concentration as a “clump”, though it could be the clump often discussed in resolved Galactic studies with a typical size of 0.3–3 pc (Bergin & Tafalla 2007) or a concentration of multiple such clumps of a larger size.

M83 is one of the closest XUV disks at $d = 4.5$ mpc (Thim et al. 2003). This prototype XUV disk has been a target of optical spectroscopy, having the low metallicity of ~1/3 Z⊙ throughout its XUV disk (Gil de Paz et al. 2007; Bresolin et al. 2009).

## 2. Hypothesis and predictions

Figure 1 presents the hypothesis (from Koda et al. 2022, with some additional explanations). It assumes that clouds in XUV disks have an internal mass distribution similar to those of Galactic star-forming clouds (such as the Orion A molecular cloud; e.g., Ikeda et al. 1999; Nakamura et al. 2019), where dense star-forming clumps, or their concentrations, are embedded in a large envelope of bulk molecular gas (Figs. 1a,b). In the cloud envelopes, CO molecules can be photo-dissociated, especially in low-metallicity environments such as those of XUV disks (Fig. 1b; see Maloney & Black 1988; van Dishoeck & Black

1988; Wolfire et al. 2010). The CO(3–2) excitation has to overcome the high $J=3$ level energy of $E_J/k \sim 33$ K and high critical density of $\ge 10^3$ cm$^{-3}$ (even after photon trapping reduces the critical density). Thus, the majority of the CO(3–2) emission must be from the dense clumps at the hearts of the clouds, which can remain CO-rich both in high and low metallicities (Fig. 1a,b). On the contrary, CO(2–1) can be excited in relatively low $n_{H_2}$ and $T_k$ in the cloud envelopes (Fig. 1a). However, the envelopes are CO-deficient (Fig. 1b) and are dark in CO(2–1) in the low metallicity of XUV disks. This can explain the CO(3–2) detection and CO(2–1) non-detection simultaneously.

This hypothesis is testable by two means. First, the sizes of the dense clumps would be smaller than the typical size of molecular clouds. For example, the Orion A cloud has a diameter of 20 pc (Nakamura et al. 2019). This size could not be resolved at the 1" (20 pc) resolution of the previous CO(3–2) observations in the M83 XUV disk (Koda et al. 2022), but it can be resolved at a higher-resolution in CO(3–2) with the new observations presented in this paper.

Second, the dense clumps should emit CO(2–1) as well, but it is expected to be faint. If the gas is thermalized at the high $n_{H_2}$ and $T_k$ required for the substantial CO(3–2) excitation, the CO 3–2/2–1 ratio in surface brightness temperature would be close to one, and the CO 3–2/2–1 intensity ratio would be $I_{CO(3-2)}/I_{CO(2-1)} = 2.25$ (see more detailed calculations in Sect. 6). This can be tested with deep CO(2–1) observations in conjunction with the measured CO(3–2) fluxes. Such CO(2–1) observations were approved in ALMA Cy9, but were not executed due largely to the loss of observing times in small array configurations as a result of the cyberattack in 2022–2023. In this work, we reevaluate the significance of the upper limit in CO(2–1) flux from the previous non-detection.

## 3. Observations and data reduction

We observed a region with an approximate diameter of 30" with the ALMA 12-m array in CO(3–2) and dust continuum emission (project # 2022.1.00359.S in Cycle 9, hereafter Cy9). It is located at a galactocentric distance of ~1.24R$_{25}$ (~8.0', ~10.4 kpc). The observed region is covered with 3 pointings in CO(3–2). Figure 2 shows the region with the sensitivity contours at 70, 50, and 30% of the peak sensitivity at the center, which is (RA, Dec)$_{J2000}$ = (13:37:05.18, -29:59:53.6). This region is
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Field of view (FoV) of Cy9 (white contours) with respect to that of Cy7 (green), displayed on the CO(3–2) peak intensity map of the Cy7 data imaged with natural weighting (CY7NA). The FoV center of Cy9 is at (RA, Dec)12000 = (13:37:05.18, −29:59:53.6). The sensitivity decreases outward from the centers for both FoVs, and the contours show the locations of 70, 50, and 30% of the peak sensitivity for the Cy7 (green) and 9 data (white). The 11 clouds within the 50% white contour are the focus of this paper. The yellow circles have 2'' (44 pc) diameters and mark the locations of the 23 molecular clouds in Koda et al. (2022). The colorbar is in units of Jy/beam.</figcaption></figure>

selected based on previous ALMA observations: starting from the CO(2–1) observations of a $2 \times 4 \text{ kpc}^2$ region (Bicalho et al. 2019, # 2013.1.00861.S in Cycle 3, Cy3), and the $1 \text{ kpc}^2$ sub-region studied in CO(3–2) (the gray scale image in Fig. 2; Koda et al. 2022, # 2017.1.00065.S in Cycle 7, Cy7). Within the $1 \text{ kpc}^2$ region, a smaller region is selected for this work. This region has bright UV and $H\alpha$ peaks and multiple molecular clouds (Fig. 2).

We analyzed the Cy3, Cy7, and Cy9 data. The details of the Cy3 and Cy7 observations are in Koda et al. (2022). In the new Cy9 observations, we configured one spectral window for the CO(3–2) line emission with band and channel widths of 937.5 MHz (812.8 km s⁻¹) and 564.5 kHz (0.4894 km s⁻¹). The other three spectral windows were configured to cover different sky frequencies for the continuum emission with a bandwidth of 1.875 GHz with 128 channels per spectral window. The central frequency of the continuum emission in the combination of the Cy7 and Cy9 data is 351.3 GHz (hereafter, 351 GHz continuum emission). The final uv-coverage at the CO(3–2) frequency ranges over angular scales of 0.23–12.2'' (~5.0–266 pc).

The visibility data are calibrated using the data reduction scripts provided by the ALMA Observatory using the Common Astronomy Software Application (CASA; McMullin et al. 2007; CASA Team 2022). The amplitude and phase of the bandpass and gain calibrators are confirmed to be flat over time and frequency after the calibrations.

For imaging, we used the Multichannel Image Reconstruction, Image Analysis, and Display package (MIRIAD; Sault et al. 1995, 1996). The CO(3–2) data from Cy7 and Cy9 were imaged separately and together. We used the natural (NA) and uniform (UN) weightings. The parameters for the imaging and final data cubes are listed in Table 1. Depending on the dataset and weighting, we refer to each data cube as CY7NA, CY7+9NA, CY7+9UN, and so on (see Table 1). The RMS noise per channel is calculated using the dirty cubes before deconvolution in the areas where the primary beam attenuation is practically negligible (PB > 0.95). Since the sensitivity and beam shape vary in an asymmetric way across the field of the Cy7 and Cy9 combined data, CY7+9NA and CY7+9UN, they are used only for reference, but not for measurements.

Table 1 includes the coefficient to convert intensity/flux units from Jy/beam to K. In this paper we mainly use Jy/beam since the sources are not resolved (only marginally resolved at best; see Sect. 5).

## 4. Cloud and clump identifications and parameters

We identified clouds and clumps hierarchically, first finding objects in the CY7NA cube and then searching for their sub-structures in the CY9NA cube by limiting the search within the volumes of the CY7NA objects, doing the same with the CY9UN cube within the volumes of the CY9NA objects. We adopted this approach because the CY7NA data have the highest sensitivity, albeit the lowest spatial resolution (see Table 1). Once the volumes of the objects were restricted by the CY7NA analysis, we relaxed search thresholds in the CY9NA and CY9UN data that have higher resolutions but lower sensitivities (see the details below). To remove the effects of sensitivity variations across the observed fields, we generated and used signal-to-noise ratio (S/N) cubes for the identifications. Unless stated otherwise, we adopted the cubes with a 2.54 km s⁻¹ channel width.

We refer to the objects identified in the CY7NA cube as clouds and those in the CY9NA and CY9UN cubes as clumps. These names reflect only the spatial resolutions of the data used for their identifications. The CY7NA cube has a resolution of 0.89'' (19 pc), approximately the diameter of a molecular cloud, and the CY9NA and CY9UN cubes have 0.40'' (8.6 pc) and 0.28'' (6.1 pc), respectively, tracing sub-cloud scales.

### 4.1. Clouds

Twenty-three molecular clouds were previously identified by Koda et al. (2022) in the same Cy7 data with an earlier data reduction. For the sake of consistency, we repeated the same identification procedure on the new CY7NA data (2.54 km s⁻¹ channel width). In the S/N cube, we first find pixels with $>5\sigma$ significance and extend their volumes down to pixels with $3\sigma$ per channel. In some cases, the envelopes of two clouds overlap, which are split with the watershed algorithm (CLUMPFIND; Williams et al. 1994). We find the same set of 23 molecular clouds as in Koda et al. (2022).

This cloud identification defines a volume of each cloud in the three-dimensional data cubes, which is used as a mask for the clump identification in the next subsection.

Hereafter, we discuss only the clouds within the Cy9 field of view (to the 50% level of the primary beam), which includes eleven clouds (Clouds 4, 7, 8, 9, 10, 11, 12, 13, 14, 15, and 17; see Fig. 2). Figures A.1 and A.2 show their cutout images in integrated intensity and peak intensity, respectively. The clouds are arranged from left to right, with their ID#s in the top-left cor-
<table><caption>**Table 1.** Parameters of the reduced data cubes.</caption><thead><tr><th rowspan="2">(1)<br>Data</th><th rowspan="2">(2)<br>Name</th><th rowspan="2">(3)<br>Δv</th><th rowspan="2">(4)<br>Cell</th><th rowspan="2">(5)<br>Wt.</th><th colspan="2">(6)<br>Beam size</th><th colspan="2">(7)<br>RMS (1σ)</th><th>(8)<br>C<sub>K</sub></th></tr><tr><th>b<sub>maj</sub>, b<sub>min</sub>, PA<br>"×", deg</th><th>√b<sub>maj</sub>b<sub>min</sub><br>pc</th><th>mJy/bm</th><th>K/[Jy/bm]</th></tr></thead><tbody><tr><td>Cy3 CO(2–1)</td><td>CY3NA</td><td>2.54</td><td>0.15</td><td>NA</td><td>0.75, 0.54, -87.47</td><td>0.64</td><td>14</td><td>10.7</td><td>56.8</td></tr><tr><td rowspan="2">Cy7 CO(3–2)</td><td rowspan="2">CY7NA</td><td>2.54</td><td>0.05</td><td>NA</td><td>0.96, 0.82, -84.05</td><td>0.89</td><td>19</td><td>1.05</td><td>13.0</td></tr><tr><td>0.85</td><td>0.05</td><td>NA</td><td>0.96, 0.82, -84.05</td><td>0.89</td><td>19</td><td>1.36</td><td>13.0</td></tr><tr><td>Cy7&9 CO(3–2)</td><td>CY7+9NA</td><td>2.54</td><td>0.05</td><td>NA</td><td>0.73, 0.67, -76.31</td><td>0.70</td><td>15</td><td>0.75</td><td>20.9</td></tr><tr><td>Cy7&9 CO(3–2)</td><td>CY7+9UN</td><td>2.54</td><td>0.05</td><td>UN</td><td>0.62, 0.57, -73.89</td><td>0.59</td><td>13</td><td>1.19</td><td>28.9</td></tr><tr><td rowspan="2">Cy9 CO(3–2)</td><td rowspan="2">CY9NA</td><td>2.54</td><td>0.05</td><td>NA</td><td>0.41, 0.38, -49.83</td><td>0.40</td><td>8.6</td><td>0.98</td><td>65.6</td></tr><tr><td>0.85</td><td>0.05</td><td>NA</td><td>0.41, 0.38, -49.83</td><td>0.40</td><td>8.6</td><td>1.44</td><td>65.6</td></tr><tr><td rowspan="2">Cy9 CO(3–2)</td><td rowspan="2">CY9UN</td><td>2.54</td><td>0.05</td><td>UN</td><td>0.30, 0.26, -44.73</td><td>0.28</td><td>6.1</td><td>1.29</td><td>131.</td></tr><tr><td>0.85</td><td>0.05</td><td>UN</td><td>0.30, 0.26, -44.73</td><td>0.28</td><td>6.1</td><td>1.90</td><td>131.</td></tr><tr><td>Cy7&9 Cont.</td><td>CY7+9CONT</td><td></td><td>0.05</td><td>NA</td><td>0.77, 0.70, -77.56</td><td>0.73</td><td>16</td><td>2.35E-2</td><td>18.4</td></tr></tbody></table>

**Notes.** (1) Data used. (2) Data name. (3) Velocity channel width. (4) Spatial cell/pixel size. (5) Adopted weighting scheme (NA=Natural, UN=Uniform weighting). (6) Convolution beam size. (7) RMS noise. The RMS is calculated in region where the primary beam attenuation, with respect to its peak across the field, is >0.95. (8) Conversion coefficient from Jy/beam to K.

ners. They are also shown from several different imaging runs in CY7NA, CY7+9NA, CY7+9UN, CY9NA, and CY9UN. These 11 clouds are listed in Table A.1 (the full list of the 23 clouds is in Table A.2). Cloud ID is labeled in the format “XX”.

None of the eleven clouds is above the 3σ detection limit in the CO(2–1) line and 351 GHz continuum emission. Figure 3 shows the S/N images calculated from the peak intensity maps: (a)(b) CO(3–2), (c) CO(2–1), and (d) 351 GHz continuum. The CO(2–1) emission is not convincingly detected even in a stack analysis of the 11 cutout images (see Appendix B).

## 4.2. Clumps (Level A and B)

We first searched the CY9NA cube, finding clumps within the cloud volumes identified in CY7NA (in the CY7NA mask), and then searched the CY9UN cube to find clumps within the volumes of the CY9NA clumps. To differentiate them, the first group of CY9NA clumps is called “Level A clumps”, and the second group, from CY9UN, is “Level B clumps”. Level B clumps reside within the volumes of Level A clumps and are their substructures. We adopted the same identification procedure as for the clouds (Sect. 4.1) but limited the search volumes.

For Level A clumps, we started from the CY9NA cube and applied the CY7NA mask to limit the search range within the volume of the CY7NA clouds. Since this is the volume where a cloud exists at a higher significance, it is less likely to pick up noise even with a lenient detection threshold. Hence, we started from >4σ peaks and extended them to the 3σ envelopes. In the CY9NA cube, we detect one or two clumps within all the clouds except Clouds 11 and 17. These two clouds are near the edge of the Cy9 field of view (Fig. 2), where the sensitivities are compromised due to the primary beam attenuation. In Table A.1, the eleven CY9NA clumps are labeled in the format “XX-YY” with “XX” as the parental cloud ID and “YY” as the clump ID within the cloud.

For Level B clumps, we did the same as for the CY9UN cube, using the same detection thresholds, but within the volumes of the CY9NA clumps as masks. The CY9UN cube has the lowest sensitivity among the CO(3–2) data cubes (Table 1) and can identify only the brightest clumps, and yet at the highest spatial resolution. We find nine CY9UN clumps. They are listed in Table A.1 with the IDs labeled in the format “XX-YY-ZZ”.

Clouds 4 and 10 are resolved to two clumps 4-1 and 4-2, and 10-1 and 10-2, respectively. Clump 4-2 is resolved further into three clumps, 4-2-1, 4-2-2, and 4-2-3.

## 4.3. Parameters

Tables A.1 and A.2 lists the measured parameters of the clouds and clumps, including their ID; positions in RA, Dec, and velocity; full width half maximum (FWHM) sizes in RA (=x), Dec (=y), and recession velocity (=v) directions; peak brightness temperature T<sub>p</sub>, peak intensity I<sub>p</sub> and its uncertainty ΔI<sub>p</sub>; the integrated flux F and its uncertainty ΔF; the velocity dispersion σ<sub>v</sub> and its uncertainty Δσ<sub>v</sub>; and the primary beam attenuation (the fraction with respect to the peak sensitivity). We used the CY7NA, CY9NA, and CY9UN cubes and masks for the measurements of the CY7NA clouds and the CY9NA and CY9UN clumps, respectively. We adopted the cubes of the 2.54 km s<sup>-1</sup> channel width for these measurements except for FWHM<sub>v</sub> (and σ<sub>v</sub>), for which we used 0.85 km s<sup>-1</sup> cubes.

As explained, cloud and clump ID are labeled in a hierarchical manner. The substructures inherit the labels of their parental structures: “XX” for the CY7NA clouds, “XX-YY” for the CY9NA (Level A) clumps, and “XX-YY-ZZ” for the CY9UN (Level B) clumps.

The coordinates are the intensity-weighted centroid positions.

The N<sub>pix</sub> is the number of the pixels identified in a cloud or clump in the RA-Dec-Vel volume (i.e., x-y-v volume). N<sub>xy</sub> is the number of the pixels in the RA-Dec projection (x-y plane). Some clumps show N<sub>pix</sub> = N<sub>xy</sub>, meaning that they are detected in one velocity channel. For reference, the areas of the beams (Table 1) correspond to N<sub>xy</sub> = 356.8 (CY7NA), 70.6 (CY9NA), and 35.4 (CY9UN), respectively.

The FWHMs are calculated from the dispersions σ<sub>i</sub><sup>obs</sup> (where i = x, y, v) within the cloud volume in the cube, assuming Gaussian profiles in space and velocity, as

$$ \text{FWHM}_i = \sqrt{8 \ln 2} \sigma_i^{\text{obs}}. \qquad (1) $$

Table A.1 lists FWHM<sub>i</sub>s instead of σ<sub>i</sub><sup>obs</sup>, because they are directly comparable to the beam sizes. We used the 0.85 km s<sup>-1</sup> cubes for FWHM<sub>v</sub>; we re-gridded the masks from the
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** S/N images in peak intensity: (a) CY9NA – CO(3–2), (b) CY7+9NA – CO(3–2), (c) CY3NA – CO(2–1), and (d) CY7+9CONT – 351 GHz continuum. The white contours show the locations of 70, 50, and 30% the peak sensitivity for the Cy9 data. The magenta circles have diameters of 2'' (44 pc) and mark the locations of the clouds; the cloud IDs are given above each circle.</figcaption></figure>

2.54 km s⁻¹ cubes and applied them to the 0.85 km s⁻¹ cubes before taking the measurements.

We note that the beam sizes are different for the CY7NA clouds, CY9NA clumps, and CY9UN clumps (Table 1). FWHMs in Tables A.1 and A.2 are the measured values before dilution corrections for the beam size and velocity resolution.

The *I*<sub>p</sub> and *F* are the values corrected for primary beam attenuation; *I*<sub>p</sub> is also written in the form of the brightness temperature, *T*<sub>p</sub>, in kelvins, converted with the Rayleigh-Jeans equation. Δ*I*<sub>p</sub> is from the RMS noise in Table 1. Δ*F* was calculated with the standard noise propagation starting from Δ*I*<sub>p</sub> by taking into account the number of beams in the volume of objects.

The S/N (*I*<sub>p</sub>/Δ*I*<sub>p</sub>) can be calculated from Table A.1. Level A (CY9NA) clumps are detected at peak S/N ratios of 4–15, mostly at S/N > 5 with the exceptions of Clumps 10-2, 15-1, and 7-1 (S/N = 4–5). Level B (CY9UN) clumps are detected less significantly, mostly in the range S/N = 4–6 except 4-1-1 (S/N = 9.9). Figures A.1 and A.2 show their integrated intensity and peak intensity maps.

From these measured parameters, we calculated the radius, *R*, and velocity dispersion, σ<sub>v</sub>, of the clouds and clumps. The observed σ<sub>v</sub><sup>obs</sup> and FWHM<sub>i</sub> contain the beam size or velocity channel width. To remove their effects, we calculated the second
moments in radius and velocity as

$$ \sigma_r = \sqrt{\sigma_x^{\text{obs}} \sigma_y^{\text{obs}} - \frac{b_{\text{maj}} b_{\text{min}}}{8 \ln 2}}, \quad (2) $$

$$ \sigma_v = \sqrt{(\sigma_v^{\text{obs}})^2 - \sigma_{\text{ch}}^2}, \quad (3) $$

and

$$ R = 1.92 \sigma_r. \quad (4) $$

The $\sigma_{\text{ch}}^2$ is the second moment of the channel window function (channel shape), and $\sigma_{\text{ch}} = w/\sqrt{12}$ for a top-hat function with a full width of $w$. We measured the $\sigma_v$ and FWHM$_v$ with the cubes of $w = 0.85 \text{ km s}^{-1}$. For $R$, we used Eq. (4), the definition by Solomon et al. (1987), which is supposed to represent the whole radial extent of the object.

Equation (2) requires the measured FWHM$_x$ and FWHM$_y$ to be greater than $b_{\text{maj}}$ and $b_{\text{min}}$ (beam size) and the $\sigma_r$ and $R$ to be positive. This condition is met only for Cloud 4 and Clump 13-1, and even their $R$ (=0.35" and 0.20" respectively) are less than half the beam sizes. Therefore, we consider that none of the clouds nor the clumps are spatially resolved in any of the measurements. We followed Koda et al. (2022) to set the upper limit of their size and considered $R$ to be measurable only when the measured size is at least 20% larger than the beam size. The minimum measurable radius is then $R_{\text{min}} = 0.54 \sqrt{b_{\text{maj}} b_{\text{min}}}$. This is 0.50" (10.9 pc), 0.22" (4.8 pc), and 0.16" (3.5 pc) for the clouds (CY7NA), Level A clumps (CY9NA), and Level B clumps (CY9UN), respectively. Obviously, these minimum radii correspond to the diameters of 1.00" (21.8 pc), 0.44" (9.6 pc), and 0.32" (7.0 pc), respectively.

The $R$ (upper limit) and $\sigma_v$ are listed in Table A.1.

# 5. Small sizes of CO(3–2) clumps

Our hypothesis is that the CO(3–2) emitting regions are dense clumps embedded in the CO-deficient cloud envelopes (Fig. 1 and Sect. 2). If this is correct, the CO(3–2) emitting regions should be compact compared to typical cloud sizes. In what follows, we show that the CO(3–2) emitting regions are indeed compact and are not spatially resolved (or at most, are only marginally resolved) in the new data, based on a consideration of the measured sizes, peak intensities, and total fluxes of the clouds and clumps.

## 5.1. Size measurements

Figures A.1 and A.2 show the clouds and clumps, from the top to bottom panels, at the lowest (0.89", 19 pc for CY7NA) to highest resolutions (0.40", 8.7 pc for CY9NA and 0.28", 6.1 pc for CY9UN). The yellow circles have a diameter of 2" (44 pc), which is close to a typical diameter of giant molecular clouds, 40 pc, in the Milky Way (Scoville & Sanders 1987). We recall that the clouds are identified in the CY7NA data, and Level A and B clumps are in the CY9NA and CY9UN data, respectively.

Clouds 4 and 10 are formally resolved into smaller clumps by the identification criteria in Sect. 4.2. None of these clumps are resolved at the higher resolutions (Table A.1). By eye, clouds 13 and 17 could also be split into smaller clumps, but are not formally resolved with the identification criteria. Their apparent clumps appear to be as small as those in cloud 4 (Figs. A.1 and A.2). The other clouds have only a single, unresolved clump (Table A.1).

Quantitatively, in Table A.1, the FWHM$_x$ and FWHM$_y$ of the clouds are similar to the beam size of the CY7NA data (0.89", 19 pc). This linear spatial scale is close to the sizes of small molecular clouds in the Milky Way disk. FWHM$_x$ and FWHM$_y$ of Level A and B clumps are also similar to the beam sizes of CY9NA and CY9UN (0.40", 9 pc and 0.28", 6 pc, respectively). Hence, they are not spatially resolved even at the higher resolutions.

From this, we conclude that the CO(3–2) emitting clumps in each cloud are not spatially resolved at the beam sizes of CY9NA and CY9UN. We note that our size measurements are not very accurate because our clump detection thresholds are not high: S/N > 4 and >3 at the peak and envelope. Within the uncertainties, the CO(3–2) sources appear to be as compact as the beam sizes of ~6–9 pc or possibly smaller.

## 5.2. Implications of the peak intensities and fluxes

Table 2 compares the peak intensities (top part) and total fluxes (bottom) between the parental clouds and their clumps (both Level A and B clumps). Each row gives measurements within each cloud. The peak intensity of the clumps, $\text{max}(I_p)$, is that of the brightest clump within the cloud. The total flux of the clumps, $\sum F$, is the sum of the fluxes of all the identified clumps in the cloud. Table 2 also lists their ratios (i.e., the fractions of emission in the clumps over that in the cloud).

Figure 4a shows the peak intensity ratio between the brightest clump in a cloud and the parental cloud (columns 4 and 6 in the top part of Table 2). The red squares are for Level A clumps, and the blue circles are for Level B clumps. In all cases, this ratio is as high as ~50–100% with the averages of 77% for Level A and 68% for Level B. The peak intensity $I_p$ represents a flux within a beam area at the peak position within an object; and the beam sizes are different between the cloud (0.89") and clump measurements (0.40" or 0.28"). The high ratios indicate that approximately the same amount of flux is enclosed within the large and small beams, and hence, the CO(3–2) sources are not much spatially resolved even with the small beams. This supports the hypothesis that the CO(3–2) emission comes from a small region within a cloud. The resolved clouds (especially #10 and 13) have relatively low fractions because the cloud flux is split into the fluxes of the multiple unresolved clumps, and this measurements take only one of them (the brightest) into account.

Figure 4b shows the fraction of the total flux of all clumps over that in their parental cloud (columns 4 and 6 in the bottom part of Table 2). An overall trend of this fraction appears different between Level A and B clumps. For Level A (red squares), the fraction is high and is mostly ~50–100%, with exceptions of clouds 8 and 15 where it is 30 and 21%, respectively. The areal ratio of the CY7NA and CY9NA beams is $(0.40"/0.89")^2 \sim 0.2$, and hence, the fraction would be 20% if the emission is uniformly distributed within the CY7NA’s 0.89" beam. The measured fractions of $\gg$20% suggest that the CO(3–2)-emitting regions are concentrated mostly in the CY9NA’s 0.40" beam. For Level B (blue circles), the fraction decreases to ~10–50% except for cloud 12. The decreasing trend from the low to high resolutions may indicate that the clumps are starting to be resolved between the resolutions of CY9NA (0.40", 8.7 pc) and CY9UN (0.28", 6.1 pc). In fact, the fraction of Level B over Level A in total flux is ~30–70% (column 7 in Table 2 bottom). This is similar to the areal ratio of the CY9NA and CY9UN beams, $(0.28"/0.40")^2 \sim 0.5$. Therefore, the CO(3–2)-emitting clumps may be being resolved at the resolutions of ~0.28–0.40" (6–9 pc).
We should note a couple of caveats. The above discussion is to capture the average characteristics of the observed clouds and clumps, and there are clearly some variations among them. For example, the small total flux fractions of Level A clumps over clouds 8 and 15 may indicate that there are some CO(3–2) emission outside the unresolved clumps. We also caution that interferometer imaging without short spacing data (or with the uniform weighting) often reduces the apparent flux. This is because the emission contained in those data is missed, and also because they carry the flux information and their weight reduction causes systematic uncertainties in the imaging process. However, the overall trends discussed above suggest that, on average, the CO(3–2)-emitting objects are smaller than the typical sizes of Galactic molecular clouds.

We draw the same conclusion from the peak surface brightness temperature $T_p$. Table A.1 shows that $T_p$ generally increases from the low-resolution (CY7NA) to high-resolution data (CY9UN), indicating that the objects are unresolved, but are approaching to be resolved. The CO emission easily becomes optically thick on theoretical and empirical bases (Goldreich & Kwan 1974; Solomon et al. 1979), and thus, we assumed that it is optically thick. Using the beam filling factor $f$ (i.e., the areal fraction of an object over beam), $T_p$ is related to the excitation temperature $T_{32}$ as $T_p = [f\Gamma(\nu_{32}, T_{32})]T_{32}$ (see Eq. (6) for $\Gamma$). $\Gamma$ is an increasing function of $T_{32}$ and can be evaluated numerically. This gives $T_p = (3.9f - 22.5f)$ K for $T_{32} = 10-30$ K with $f$ as an unknown factor. Compared to $T_p$ in Table A.1 (the brightest is Clump 4-1-1 with $T_p = 2.49$ K), the filling factor is $f = 0.64$ or smaller, again indicating unresolved or marginally resolved objects. The CO emission rarely becomes optically thin, except under some very special conditions (e.g., line wings and extremely low metallicity); thus, we did not consider it further but briefly note that an optically thin case would reduce the apparent $T_p$ and, hence, make the above $f$ a lower limit.

## 6. Non-detection in CO(2–1)

The combination of the detection in CO(3–2) and non-detection in CO(2–1) also constrains the cloud structure and physical condition. The new data reduction of the CO(2–1) also supports this.

Koda et al. (2022) assumed, as in Fig. 1b, that [1] the cloud envelope is CO-deficient due to the low metallicity and does not emit much CO emission, while the dense clump(s) in the cloud is still CO-rich and emit CO(3–2) as well as CO(1–0) and CO(2–1), and [2] the CO-rich part of the cloud (i.e., the dense clumps) is close to being thermalized, and hence, the CO 3–2/2–1 ratio in surface brightness temperature is relatively high and is approximately unity. Under these assumptions, the CO(2–1) flux can be predicted from the observed CO(3–2) flux. Koda et al. (2022) concluded that the predicted CO(2–1) flux of the brightest cloud is just below the detection limit of the previous CO(2–1) observations (Bicalho et al. 2019). We note that in this section, we discuss “flux” instead of “surface brightness temperature” because the emission sources are not resolved.

The assumptions that Koda et al. (2022) adopted predict the minimum of the CO(2–1) flux. In other words, if any of the assumptions deviate significantly from the reality, the predicted CO(2–1) flux would be higher and should have been detected in the previous CO(2–1) study.

The CO(2–1) flux of a CO(3–2)-emitting dense clump, without the cloud envelope, can be estimated from its CO(3–2) flux. We assumed that the beam filling factors of the two transitions are equal. Using the flux $F_i$, excitation temperature $T_i$, and fre-

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>Fig. 4. Clumps over cloud ratios for Level A (red square) and B clumps (blue circle) measured in CO(3–2). (a) Peak intensity ratio. The peak intensity of the brightest clump, $\max(I_p)$ in mJy/beam(1), is compared with the peak intensity of its parental cloud, $I_p$, in mJy/beam(2). Note that the two beams (1 and 2) have different sizes. (b) Total flux ratio. On the right edge of each panel, the triangles pointing left show the averages of the Level A (red) and B (blue) clumps.</figcaption></figure>

quency $v_i$ of the CO(2–1) and CO(3–2) emission ($i = 21$ or 32, respectively), the CO 3–2/2–1 flux ratio is

$$ \frac{F_{21}}{F_{32}} = \left(\frac{v_{21}}{v_{32}}\right)^2 \left(\frac{T_{21}}{T_{32}}\right) \left(\frac{\Gamma(\nu_{21}, T_{21})}{\Gamma(\nu_{32}, T_{32})}\right) = \frac{4}{9} \left(\frac{T_{21}}{T_{32}}\right) \alpha(T_{21}, T_{32}), \quad (5) $$

where

$$ \Gamma(\nu, T) = \frac{h\nu/kT}{\exp(h\nu/kT) - 1}, \quad (6) $$

and we defined

$$ \alpha(T_{21}, T_{32}) = \frac{\Gamma(\nu_{21}, T_{21})}{\Gamma(\nu_{32}, T_{32})}. \quad (7) $$

The “$T_{21}/T_{32}$” term in Eq. (5) takes the minimum under the thermalized condition ($T_{21} = T_{32}$) with respect to non-thermalized conditions ($T_{21} > T_{32}$). The “$\alpha$” term can be evaluated numerically and is always $\alpha(T_{21} > T_{32}) > \alpha(T_{21} = T_{32})$, taking the minimum in the thermalized condition. Therefore, using $\alpha_{th}(T) \equiv \alpha(T \equiv T_{21} = T_{32})$, the thermalized condition always
<table><caption>**Table 2.** Comparison of clouds and Level A and B clumps in terms of peak intensity and total flux.</caption><thead><tr><th>(1)<br>Cloud ID</th><th>(2)<br>Cloud</th><th>(3)<br>Level A</th><th>(4)<br>(3)/(2)</th><th>(5)<br>Level B</th><th>(6)<br>(5)/(2)</th><th>(7)<br>(5)/(3)</th></tr><tr><td></td><td><i>I</i><sub>p</sub><br>mJy/bm</td><td><i>max</i><sub>A</sub>(<i>I</i><sub>p</sub>)<br>(mJy/bm)</td><td><i>max</i><sub>A</sub>(<i>I</i><sub>p</sub>)/<i>I</i><sub>p</sub></td><td><i>max</i><sub>B</sub>(<i>I</i><sub>p</sub>)<br>mJy/bm</td><td><i>max</i><sub>B</sub>(<i>I</i><sub>p</sub>)/<i>I</i><sub>p</sub></td><td><i>max</i><sub>B</sub>(<i>I</i><sub>p</sub>)/<i>max</i><sub>A</sub>(<i>I</i><sub>p</sub>)</td></tr></thead><tbody><tr><td>4</td><td>26.47 ± 1.25</td><td>22.21 ± 1.47</td><td>0.84 ± 0.07</td><td>19.04 ± 1.92</td><td>0.72 ± 0.08</td><td>0.86 ± 0.10</td></tr><tr><td>7</td><td>5.53 ± 1.08</td><td>4.92 ± 1.06</td><td>0.89 ± 0.26</td><td></td><td></td><td></td></tr><tr><td>8</td><td>7.26 ± 1.17</td><td>5.21 ± 1.04</td><td>0.72 ± 0.18</td><td></td><td></td><td></td></tr><tr><td>9</td><td>7.80 ± 1.08</td><td>6.79 ± 1.02</td><td>0.87 ± 0.18</td><td>6.16 ± 1.33</td><td>0.79 ± 0.20</td><td>0.91 ± 0.23</td></tr><tr><td>10</td><td>11.64 ± 1.17</td><td>7.56 ± 1.00</td><td>0.65 ± 0.11</td><td>6.49 ± 1.31</td><td>0.56 ± 0.13</td><td>0.86 ± 0.19</td></tr><tr><td>12</td><td>9.74 ± 1.11</td><td>7.59 ± 1.00</td><td>0.78 ± 0.14</td><td>6.70 ± 1.30</td><td>0.69 ± 0.15</td><td>0.88 ± 0.20</td></tr><tr><td>13</td><td>17.57 ± 1.16</td><td>11.15 ± 1.01</td><td>0.63 ± 0.07</td><td>7.63 ± 1.32</td><td>0.43 ± 0.08</td><td>0.68 ± 0.13</td></tr><tr><td>14</td><td>6.81 ± 1.22</td><td>6.64 ± 1.09</td><td>0.98 ± 0.24</td><td>5.99 ± 1.42</td><td>0.88 ± 0.26</td><td>0.90 ± 0.27</td></tr><tr><td>15</td><td>12.09 ± 1.34</td><td>6.38 ± 1.42</td><td>0.53 ± 0.13</td><td></td><td></td><td></td></tr><tr><td></td><td><i>F</i><br>mJy km s<sup>-1</sup></td><td><i>&Sigma;</i><sub>A</sub><i>F</i><br>mJy km s<sup>-1</sup></td><td></td><td><i>&Sigma;</i><sub>B</sub><i>F</i><br>mJy km s<sup>-1</sup></td><td></td><td></td></tr><tr><td>4</td><td>258.20 ± 11.20</td><td>196.00 ± 12.60</td><td>0.76 ± 0.06</td><td>104.50 ± 11.10</td><td>0.40 ± 0.05</td><td>0.53 ± 0.06</td></tr><tr><td>7</td><td>5.80 ± 2.00</td><td>4.40 ± 1.80</td><td>0.76 ± 0.41</td><td></td><td></td><td></td></tr><tr><td>8</td><td>30.50 ± 4.70</td><td>9.30 ± 2.50</td><td>0.30 ± 0.09</td><td></td><td></td><td></td></tr><tr><td>9</td><td>30.20 ± 4.40</td><td>17.00 ± 3.20</td><td>0.56 ± 0.13</td><td>8.20 ± 2.80</td><td>0.27 ± 0.10</td><td>0.48 ± 0.18</td></tr><tr><td>10</td><td>49.60 ± 5.70</td><td>23.60 ± 3.80</td><td>0.48 ± 0.09</td><td>6.30 ± 2.30</td><td>0.13 ± 0.05</td><td>0.27 ± 0.10</td></tr><tr><td>12</td><td>27.40 ± 4.00</td><td>28.00 ± 3.90</td><td>1.02 ± 0.21</td><td>19.00 ± 4.00</td><td>0.69 ± 0.18</td><td>0.68 ± 0.17</td></tr><tr><td>13</td><td>78.60 ± 6.40</td><td>74.50 ± 6.10</td><td>0.95 ± 0.11</td><td>36.90 ± 5.70</td><td>0.47 ± 0.08</td><td>0.50 ± 0.09</td></tr><tr><td>14</td><td>13.30 ± 3.30</td><td>9.70 ± 2.50</td><td>0.73 ± 0.26</td><td>5.20 ± 2.30</td><td>0.39 ± 0.20</td><td>0.54 ± 0.27</td></tr><tr><td>15</td><td>35.80 ± 5.20</td><td>7.40 ± 2.70</td><td>0.21 ± 0.08</td><td></td><td></td><td></td></tr></tbody></table>

**Notes.** Peak and integrated intensities for clouds and clumps. (1) Cloud ID. (2) Peak intensity $I_p$ or total flux $F$, the same as those in Table A.1. (3)(5) Maximum $I_p$ among clumps or total flux $F$ of all clumps for Level A and B clumps. (4)(6)(7) Ratios of Level A/Cloud, Level B/Cloud, and Level B/Level A. The beam sizes for Cloud (~0.89"), Level A clumps (~0.40"), and Level B clumps (~0.28") are different.

predicts the minimum CO(2–1) flux of $F_{21} = (4/9\alpha_{th}(T))F_{32}$, where $\alpha_{th} \approx 1.10, 1.16$, and 1.40 for $T = 30, 20$, and 10 K, respectively. Here, we neglected the contribution of the cosmic microwave background radiation since its temperature is low for the CO(3–2)-emitting gas.

The brightest cloud (Cloud 4 in Table A.1) has the CO(3–2) flux of $F_{32} = 258.2$ mJy km s⁻¹. In the thermalized condition, this translates to the minimum expected CO(2–1) flux of $F_{21} = 126.2, 133.1$, and $160.7$ mJy km s⁻¹ at 30, 20, and 10 K (Eq. (5)). This cloud has FWHM$_v = 6.21$ km s⁻¹ in CO(3–2), and if this is also the width in CO(2–1), the minimum predicted CO(2–1) intensity is $I_{21} = 20.3, 21.4$, and $25.9$ mJy. These are at 1.9, 2.0, and $2.4\sigma$ significance, respectively. Thus, it is slightly below the detection limit of the existing CY3NA data.

If the clump is not thermalized, the predicted CO(2–1) flux becomes higher and should have been detected. For example, if the CO(3–2)-emitting region has $T_{32}/T_{21} \sim 0.5$ (i.e., the average in the main disks of nearby galaxies; Leroy et al. 2022), $\alpha(T_{21}, T_{32}) = \alpha(T, T/2) = 1.51, 1.92$, and $4.39$ at $T = 30, 20$, and 10 K. Cloud 4 would have $F_{21} = 346.6, 440.7$, and $1007.6$ mJy km s⁻¹ and $I_{21} = 55.8, 71.0$, and $162.2$ mJy. It should have been detected in CO(2–1) at 5.2, 6.6, and $15.2\sigma$ in the CY3NA data. Therefore, the excitation condition of the clump cannot deviate significantly from the thermalized condition, given the non-detection in CO(2–1).

In addition, if the envelope of the cloud is filled with CO, the CO(2–1) flux would be higher since the CO(2–1) emission is easily excited in the average density and temperature of typical cloud envelopes. This should also have made the clump detectable in the previous study.

In summary, the non-detection in CO(2–1) supports the cloud structure presented in Fig. 1b. This is obviously a sim-

plistic picture, and the details will need some modification in the future. However, as a first order approximation, Fig. 1b likely represents the cloud structure in the XUV disk.

## 7. The size–velocity dispersion relation

Figure 5 plots $\sigma_v$ against the upper limits of $R$ for the clouds (arrows at $\log R = 1.03$), Level A (0.68) and B clumps (0.54). For reference, it also shows clouds in the Milky Way, Large and Small Magellanic Clouds (LMC and SMC), and the dwarf galaxy WLM (Solomon et al. 1987; Bolatto et al. 2008; Wong et al. 2011; Rubio et al. 2015). The clouds in the XUV disk of M83 appear consistent with those in the dwarfs, while their radii are only upper limits.

There is also a caveat in this comparison, since these reference clouds are observed in the lower-excitation transitions, CO(1–0) or CO(2–1). However, we note that the three CO transitions ($J = 1–0, 2–1$, and $3–2$) should be emitted from a similar central region (Fig. 1b) if the clouds and clumps share a similar physical and chemical structure in the low-metallicity environments (i.e., the Milky Way outskirts, LMC, SMC, and other dwarfs, as well as the XUV disk). In fact, molecular clouds in LMC and SMC appear compact in recent CO(1–0) and CO(2–1) measurements (Harada et al. 2019; Tokuda et al. 2021; Ohno et al. 2023). Therefore, the consistency in Fig. 5 supports our hypothesis (Fig. 1), of course, within the limitation that we have only the upper limits of $R$.

## 8. Further discussion

Koda et al. (2022) presented a hypothesis on the mass structure and CO chemistry in molecular clouds in the XUV disk
<figure><img src="image_7.png" /><figcaption>Fig. 5. $\sigma_v-R$ plot. M83 XUV disk clouds and clumps are in blue (and the blue arrows are the upper limits in radius). The clouds have the largest radii, followed by the Level A and then the Level B clumps. The reference points in the background are for the Milky Way's inner disk (Solomon et al. 1987), outer disk ($r_{gal} \sim 12-14$ kpc; Sun et al. 2017), and far outer disk ($\sim 14-22$ kpc; Sun et al. 2015), the LMC and SMC (Wong et al. 2011; Bolatto et al. 2008), and WLM (Rubio et al. 2015). The solid line is a fit to the inner disk clouds (Solomon et al. 1987).</figcaption></figure>

(Fig. 1). The Orion A cloud (i.e., a nearby Galactic molecular cloud) exhibits SF activity similar to those around the XUV disk clouds, and it consists of a relatively low-density envelope of ~20–30 pc with a dense CO(3–2)-bright region of ~10 pc in size (Ikeda et al. 1999; Nakamura et al. 2019). The hypothesis suggests that the XUV clouds also have a similar mass structure, but in their envelopes, CO molecules are photo-dissociated and deficient due to the low metallicity condition. The dissociating UV radiation is in the XUV disk (see Koda et al. 2022; there are multiple Orion-like sites of massive SF-traced by UV and Hα sources – in the 1 kpc² region, i.e., the Cy7 field of view). This study confirms the compactness of the CO(3–2)-emitting regions (≤6–9 pc in diameter) in the XUV disk clouds. In the Orion A cloud, the dense CO(3–2)-bright region consists of multiple smaller clumps (see Ikeda et al. 1999), while our resolution is not high enough to resolve such substructures. This study also shows that the CO(3–2)-emitting regions in the XUV disk clouds must be close to the thermalized condition, since otherwise the expected CO(2–1) flux would be brighter and should have been detected in the past ALMA observations.

## 8.1. Mass and radius of the CO(3–2)-emitting regions

The masses of the clumps (i.e., CO(3–2)-emitting regions) can be crudely estimated (Koda et al. 2022), assuming that (1) the clumps have CO molecules (as in Fig. 1b), (2) the Galactic CO-to-H₂ conversion factor $X_{CO(1-0)} = 3 \times 10^{20} \text{ cm}^{-2} [\text{K km s}^{-1}]^{-1}$ can be applied to the clumps (i.e., parts of the clouds)<fnref n="1" />, and (3) the gas is thermalized and hence the CO 3–2/1–0 brightness temperature ratio is unity ($R_{3-2/1-0} = 1$). With these assumptions,

the clump mass, including helium and other heavier elements (i.e., multiplying the total hydrogen mass by 1.37), is

$$ M_{\text{clump}} = 3.6 \times 10^3 \text{ M}_{\odot} \left( \frac{F_{\text{CO}(3-2)} dv}{100 \text{ mJy km s}^{-1}} \right) \left( \frac{D}{4.5 \text{ Mpc}} \right)^2 \times \left( \frac{R_{3-2/1-0}}{1.0} \right)^{-1} \left( \frac{X_{\text{CO}(1-0)}}{3.0 \times 10^{20}} \right). \quad (8) $$

The measured CO(3–2) fluxes (Table A.1) give $M_{\text{clump}} = 2.1 \times 10^2-9.3 \times 10^3 \text{ M}_{\odot}$.

These $M_{\text{clump}}$ can be translated, also crudely, to the radius of clumps (or clump concentrations). The CO(3–2) excitation requires a density of $n_{H_2} \sim 10^3 \text{ cm}^{-3}$, to within an order of magnitude. By assuming this density as an average density in a spherical CO(3–2)-emitting region, the clump mass within its radius, including helium and other heavier elements, is

$$ M_{\text{clump}} = 3.8 \times 10^3 \text{ M}_{\odot} \left( \frac{R}{3 \text{ pc}} \right)^3 \left( \frac{n_{\text{H}_2}}{10^3 \text{ cm}^{-3}} \right). \quad (9) $$

The range of the clump masses estimated above then corresponds to $R = 1–4$ pc. This is consistent with the new measurements in this study. The clumps (or their concentrations) are compact and not spatially resolved at the 6–9 pc resolutions, while there is an indication that they may be starting to be resolved at the 6 pc resolution.

These mass estimations easily have a factor of 2–3 uncertainty. It is also worth noting that, within such uncertainties, they are roughly consistent with the virial mass $M_{vir} (= 9/2(R\sigma_v^2/G)$ for a radial density profile of $\propto 1/r$; see Lee et al. (2024), which is

$$ M_{\text{vir}} = 3 \times 10^3 \text{ M}_{\odot} \left( \frac{R}{3 \text{ pc}} \right) \left( \frac{\sigma_v}{1 \text{ km s}^{-1}} \right)^2, \quad (10) $$

using rounded values of $R = 3$ pc and $\sigma_v = 1 \text{ km s}^{-1}$ (see Table A.1).

## 8.2. Cloud evolution and distribution

The cloud structure in Fig. 1 does not take into account the evolution of molecular clouds. It is possible that the clouds evolve and develop dense clumps only at a later stage of their lifetime. In this case, the CO(3–2) observations would selectively trace the clouds at the evolved stage when they are ready for SF. In fact, in the inner disk of M83, the populations of star-forming molecular clouds are mostly localized around the galactic center, bar, and spiral arms. Molecular gas and clouds are abundant even in the interarm regions, but appear to be dormant in SF (Koda et al. 2023). If this is also the case in XUV disks, there could also be (undetected) clouds at an earlier stage of their evolution before the dense clumps develop in their hearts. In order to understand the SF in XUV disks and in general, it is important to confirm the presence or absence of such dormant cloud populations and to characterize their structures and distribution as a function of their environments (e.g., the amount of the surrounding HI gas, spiral arms, UV radiation field, stellar feedback, etc.).

The large-scale distribution of the clouds may also provide a clue on their formation and evolution (e.g., Dobbs & Baba 2014). Two additional clouds emerge at the highest sensitivity part of the CY7+9NA data (see Fig. 3b), while we focused on the clouds and clumps discovered in the CY7NA data. One is detected at $6.3\sigma$ between Clouds 9 and 12 at a peak position

<footnote n="1">¹ We assumed this $X_{CO(1-0)}$ for clumps. If their parental clouds have the structure in Fig. 1b, their envelopes should have a fair amount of mass while no/little CO emission. Thus, $X_{CO(1-0)}$ for clouds should be larger than the assumed, which is the often-discussed metallicity dependence of $X_{CO(1-0)}$. Recently, Lee et al. (2024) also pointed out a possible dependence of $X_{CO(1-0)}$ on cloud population.</footnote>
of (RA, Dec) = (13:37:05.29, −29:59:54.1), and the other is at
5.3σ between Clouds 7 and 8 at (13:37:05.18, −29:59:50.7).
The CY7+9NA data, however, have a nonuniform, asymmet-
ric sensitivity distribution (Fig. 2), and thus, we did not include
them in the sample for our analysis. Nevertheless, they exem-
plify the fact that more CO(3−2) sources could emerge with
longer integration. Koda et al. (2022) noted that the CO(3−2)
sources are distributed along filamentary structures as chains
(like seeds in edamame green soybeans). The two CY7+9NA
clouds appear to support this view. The distribution may indicate
that their formation requires compression due to large-scale gas
flows.

9. Summary

Massive SF sites have been detected in XUV disks, but their parental gas clouds remained mostly undetected despite searches in CO(1–0) and CO(2–1). Therefore, the recent discovery of 23 molecular clouds in CO(3–2) in the XUV disk of M83 is startling. To explain the previous non-detections and the new detections, Koda et al. (2022) introduced the hypothesis that clouds in XUV disks, like Galactic star-forming clouds, have dense star-forming clumps embedded in a large envelope of bulk molecular gas, where CO molecules are photo-dissociated due to the lack of CO self-shielding and dust extinction in the low-metallicity environment of XUV disks (Figs. 1a,b). CO(3–2) can be excited preferentially in the dense clumps but not as much in the less dense envelopes. This can explain why CO(3–2) is detected when CO(1–0) and CO(2–1) levels are significantly reduced.

This paper confirms the compactness of the dense, CO(3–2)-emitting clumps in the clouds as predicted by this hypothesis. In the new ALMA data, they are not resolved at a 0.40" (9 pc) resolution (natural weighting), nor at a 0.28" (6 pc) resolution (uniform weighting), but there is some indication that they may be starting to be resolved at the 6 pc resolution. These are similar to the sizes of the dense central clumps (or their concentration) in the Orion A molecular cloud in the solar neighborhood, where a similar amount of massive SF is ongoing.

The CO(3–2) to CO(2–1) and CO(3–2) to CO(1–0) line ratios in the dense clumps must be higher than the average ratio observed in normal galactic disks (e.g., T32/T10 ~ 1/3 and T32/T21 ~ 1/2; Wilson et al. 2009; Vlahakis et al. 2013; Leroy et al. 2022). Otherwise, the previous study (Bicalho et al. 2019) would have detected the CO(2–1) emission. Likewise, the molecular envelopes around the dense clumps within the clouds cannot have much CO, as the CO(1–0) and CO(2–1) emission should be excited easily in the envelopes and, thus, would have been detected. The CO(3–2)-emitting regions in the XUV clouds must be close to thermalized since, otherwise, the CO(2–1) flux would be brighter and thus would have been detected in past observations.

The cloud structure hypothesized here is, of course, simplis-
tic and will need some adjustments in the future. For example,
it does not include the possible evolution of clouds from a non-
star-forming phase to a star-forming phase. The clouds detected
in CO(3-2) may trace only the clouds in the star-forming phase.
For the clouds in the XUV disk discussed in this paper, this
model appears to capture the reality within the limitation of
the current observations.

**Acknowledgements.** JK thanks his colleagues at Paris Observatory for their hospitality, both professionally and personally, during his stay on sabbatical, and colleagues at the Institut d'Astrophysique de Paris for stimulating discussions. We also thank the anonymous referee. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2013.1.00861.S, 2017.1.00065.S, and 2022.1.00359.S. 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. Data analysis is in part carried out on the Multi-wavelength Data Analysis System operated by the Astronomy Data Center (ADC) in NAOJ. J.K. acknowledges support from NSF through grants AST-1812847 and AST-2006600. M.R. wishes to acknowledge support from ANID(CHILE) through FONDE-CYT grant No1190684.

References

Bergin, E. A., & Tafalla, M. 2007, ARA&A, 45, 339

Bicalho, I. C., Combes, F., Rubio, M., Verdugo, C., & Salome, P. 2019, A&A, 623, A66

Bigiel, F., Leroy, A., Seibert, M., et al. 2010, ApJ, 720, L31

Bolatto, A. D., Leroy, A. K., Rosolowsky, E., Walter, F., & Blitz, L. 2008, ApJ, 686, 948

Braine, J., & Herpin, F. 2004, Nature, 432, 369

Braine, J., Ferguson, A. M. N., Bertoldi, F., & Wilson, C. D. 2007, ApJ, 669, L73

Braine, J., Gratier, P., Kramer, C., et al. 2010, A&A, 518, L69

Braine, J., Gratier, P., Contreras, Y., Schuster, K. F., & Brouillet, N. 2012, A&A, 548, A52

Bresolin, F., Ryan-Weber, E., Kennicutt, R. C., & Goddard, Q. 2009, ApJ, 695, 580

CASA Team (Bean, B., et al.) 2022, PASP, 134, 114501

Dessauges-Zavadsky, M., Verdugo, C., Combes, F., & Pfenniger, D. 2014, A&A, 566, A147

Dobbs, C., & Baba, J. 2014, PASA, 31, e035

Eibensteiner, C., Bigiel, F., Leroy, A. K., et al. 2023, A&A, 675, A37

Gil de Paz, A., Madore, B. F., Boissier, S., et al. 2005, ApJ, 627, L29

Gil de Paz, A., Boissier, S., Madore, B. F., et al. 2007, ApJS, 173, 185

Goldreich, P., & Kwan, J. 1974, ApJ, 189, 441

Harada, R., Onishi, T., Tokuda, K., et al. 2019, PASJ, 71, 44

Ikeda, M., Maezawa, H., Ito, T., et al. 1999, ApJ, 527, L59

Koda, J., Yagi, M., Boissier, S., et al. 2012, ApJ, 749, 20

Koda, J., Watson, L., Combes, F., et al. 2022, ApJ, 941, 3

Koda, J., Hirota, A., Egusa, F., et al. 2023, ApJ, 949, 108

Lee, A. M., Koda, J., Hirota, A., Egusa, F., & Heyer, M. 2024, ApJ, 968, 97

Leroy, A. K., Rosolowsky, E., Usero, A., et al. 2022, ApJ, 927, 149

Maloney, P., & Black, J. H. 1988, ApJ, 325, 389

McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, ASP Conf. Ser., 376, 127

Nakamura, F., Ishii, S., Dobashi, K., et al. 2019, PASJ, 71, S3

Ohno, T., Tokuda, K., Konishi, A., et al. 2023, ApJ, 949, 63

Rautio, R. P. V., Watkins, A. E., Salo, H., et al. 2024, A&A, 681, A76

Rubio, M., Elmegreen, B. G., Hunter, D. A., et al. 2015, Nature, 525, 218

Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, ASP Conf. Ser., 77, 433

Sault, R. J., Staveley-Smith, L., & Brouw, W. N. 1996, A&AS, 120, 375

Scoville, N. Z., & Sanders, D. B. 1987, Astrophys. Space Sci. Lib., 134, 21

Solomon, P. M., Sanders, D. B., & Scoville, N. Z. 1979, ApJ, 232, L89

Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 319, 730

Sun, Y., Xu, Y., Yang, J., et al. 2015, ApJ, 798, L27

Sun, Y., Su, Y., Zhang, S.-B., et al. 2017, ApJS, 230, 17

Thilker, D. A., Bianchi, L., Boissier, S., et al. 2005, ApJ, 619, L79

Thilker, D. A., Bianchi, L., Meurer, G., et al. 2007, ApJS, 173, 538

Thim, F., Tammann, G. A., Saha, A., et al. 2003, ApJ, 590, 256

Tokuda, K., Kondo, H., Ohno, T., et al. 2021, ApJ, 922, 171

van Dishoeck, E. F., & Black, J. H. 1988, ApJ, 334, 771

Vlahakis, C., van der Werf, P., Israel, F. P., & Tilanus, R. P. J. 2013, MNRAS, 433, 1837

Watson, L. C., & Koda, J. 2017, in Outskirts of Galaxies, eds. J. H. Knapen, J. C. Lee, & A. Gil de Paz, 434, 175

Williams, J. P., de Geus, E. J., & Blitz, L. 1994, ApJ, 428, 693

Wilson, C. D., Warren, B. E., Israel, F. P., et al. 2009, ApJ, 693, 1736

Wolfire, M. G., Hollenbach, D., & McKee, C. F. 2010, ApJ, 716, 1191

Wong, T., Hughes, A., Ott, J., et al. 2011, ApJS, 197, 16
**Appendix A: Parameters and cutout images of
clouds and clumps**

Table A.1 lists the measured parameters of the clouds and clumps (see Sect. 4.3 for explanations). Figures A.1 and A.2 show their cutout images in integrated intensity and peak intensity, respectively.

Table A.2 is the same as the "Clouds" section of Table A.1,
but lists all 23 clouds identified in the CY7NA cube. This table
does not include the peak brightness temperature Tp. The Ip and
Tp are related as Tp = CKIp (see Table 1 for CK).
<table><caption>Table A.1. Measured parameters of clouds and Level A and B clumps in CO(3-2)</caption><thead><tr><th>(1)</th><th>(2)</th><th>(3)</th><th>(4)</th><th>(5)</th><th>(6)</th><th>(7)</th><th>(8)</th><th>(9)</th><th>(10)</th><th>(11)</th><th>(12)</th><th>(13)</th><th>(14)</th><th>(15)</th><th>(16)</th><th>(17)</th></tr><tr><th>ID</th><th>RA<br>J2000</th><th>Dec<br>J2000</th><th>V<br>km s<sup>-1</sup></th><th>N<sub>pix</sub></th><th>N<sub>xy</sub></th><th>x<br>"</th><th>y<br>"</th><th>v<br>km s<sup>-1</sup></th><th>T<sub>p</sub><br>K</th><th>I<sub>p</sub><br>km y/bm</th><th>ΔI<sub>p</sub><br>km y/bm</th><th>F<br>km y km s<sup>-1</sup></th><th>ΔF<br>km y km s<sup>-1</sup></th><th>R<br>"</th><th>σ<sub>v</sub><br>km s<sup>-1</sup></th><th>Δσ<sub>v</sub><br>km s<sup>-1</sup></th><th>PB</th></tr><tr><th colspan="18">Clouds (CY7NA)</th></tr></thead><tbody><tr><td>4</td><td>13:37:04.984</td><td>-29:59:46.15</td><td>564.4</td><td>3973</td><td>1438</td><td>1.02</td><td>0.96</td><td>6.21</td><td>0.34</td><td>26.47</td><td>1.25</td><td>258.2</td><td>11.2</td><td>&lt;0.50</td><td>2.62</td><td>0.06</td><td>0.79</td></tr><tr><td>7</td><td>13:37:05.337</td><td>-29:59:50.60</td><td>559.8</td><td>192</td><td>192</td><td>0.50</td><td>0.40</td><td>1.38</td><td>0.07</td><td>5.53</td><td>1.08</td><td>5.8</td><td>2.0</td><td>&lt;0.50</td><td>0.52</td><td>0.15</td><td>0.97</td></tr><tr><td>8</td><td>13:37:05.054</td><td>-29:59:51.21</td><td>565.9</td><td>856</td><td>684</td><td>1.07</td><td>0.58</td><td>3.23</td><td>0.09</td><td>7.26</td><td>1.17</td><td>30.5</td><td>4.7</td><td>&lt;0.50</td><td>1.35</td><td>0.08</td><td>0.88</td></tr><tr><td>9</td><td>13:37:05.352</td><td>-29:59:52.77</td><td>561.2</td><td>904</td><td>698</td><td>0.81</td><td>0.84</td><td>2.82</td><td>0.10</td><td>7.80</td><td>1.08</td><td>30.2</td><td>4.4</td><td>&lt;0.50</td><td>1.17</td><td>0.08</td><td>0.97</td></tr><tr><td>10</td><td>13:37:05.040</td><td>-29:59:53.04</td><td>567.5</td><td>1266</td><td>935</td><td>0.80</td><td>1.04</td><td>3.36</td><td>0.15</td><td>11.64</td><td>1.17</td><td>49.6</td><td>5.7</td><td>&lt;0.50</td><td>1.41</td><td>0.10</td><td>0.87</td></tr><tr><td>11</td><td>13:37:05.917</td><td>-29:59:55.36</td><td>564.4</td><td>752</td><td>572</td><td>0.71</td><td>0.71</td><td>2.21</td><td>0.12</td><td>9.13</td><td>1.06</td><td>27.1</td><td>3.9</td><td>&lt;0.50</td><td>0.91</td><td>0.13</td><td>0.99</td></tr><tr><td>12</td><td>13:37:05.186</td><td>-29:59:55.41</td><td>565.5</td><td>700</td><td>494</td><td>0.74</td><td>0.58</td><td>2.60</td><td>0.13</td><td>9.74</td><td>1.11</td><td>27.4</td><td>4.0</td><td>&lt;0.50</td><td>1.07</td><td>0.11</td><td>0.93</td></tr><tr><td>13</td><td>13:37:05.084</td><td>-29:59:55.69</td><td>567.6</td><td>1607</td><td>748</td><td>0.82</td><td>0.69</td><td>3.64</td><td>0.23</td><td>17.57</td><td>1.16</td><td>78.6</td><td>6.4</td><td>&lt;0.50</td><td>1.52</td><td>0.08</td><td>0.88</td></tr><tr><td>14</td><td>13:37:04.978</td><td>-29:59:56.98</td><td>567.8</td><td>388</td><td>313</td><td>0.59</td><td>0.49</td><td>2.34</td><td>0.09</td><td>6.81</td><td>1.22</td><td>13.3</td><td>3.3</td><td>&lt;0.50</td><td>0.96</td><td>0.17</td><td>0.85</td></tr><tr><td>15</td><td>13:37:04.820</td><td>-29:59:59.49</td><td>570.4</td><td>763</td><td>607</td><td>0.75</td><td>0.78</td><td>2.32</td><td>0.16</td><td>12.09</td><td>1.34</td><td>35.8</td><td>5.2</td><td>&lt;0.50</td><td>0.94</td><td>0.11</td><td>0.75</td></tr><tr><td>17</td><td>13:37:04.738</td><td>-30:00:01.06</td><td>569.6</td><td>754</td><td>626</td><td>0.74</td><td>0.83</td><td>2.06</td><td>0.13</td><td>9.76</td><td>1.42</td><td>35.6</td><td>5.5</td><td>&lt;0.50</td><td>0.85</td><td>0.10</td><td>0.70</td></tr><tr><th colspan="18">Level A Clumps (CY7NA)</th></tr><tr><td>4-1</td><td>13:37:05.010</td><td>-29:59:46.20</td><td>563.0</td><td>230</td><td>154</td><td>0.33</td><td>0.39</td><td>3.21</td><td>1.46</td><td>22.21</td><td>1.47</td><td>79.3</td><td>7.1</td><td>&lt;0.22</td><td>1.34</td><td>0.09</td><td>0.63</td></tr><tr><td>4-2</td><td>13:37:04.971</td><td>-29:59:46.22</td><td>565.8</td><td>452</td><td>250</td><td>0.32</td><td>0.66</td><td>5.23</td><td>0.88</td><td>13.40</td><td>1.49</td><td>116.7</td><td>10.3</td><td>&lt;0.22</td><td>2.21</td><td>0.07</td><td>0.61</td></tr><tr><td>7-1</td><td>13:37:05.325</td><td>-29:59:50.43</td><td>559.8</td><td>31</td><td>31</td><td>0.20</td><td>0.17</td><td>1.46</td><td>0.32</td><td>4.92</td><td>1.06</td><td>4.4</td><td>1.8</td><td>&lt;0.22</td><td>0.57</td><td>0.16</td><td>0.92</td></tr><tr><td>8-1</td><td>13:37:05.049</td><td>-29:59:51.22</td><td>565.0</td><td>64</td><td>61</td><td>0.23</td><td>0.28</td><td>1.90</td><td>0.34</td><td>5.21</td><td>1.04</td><td>9.3</td><td>2.5</td><td>&lt;0.22</td><td>0.77</td><td>0.21</td><td>0.94</td></tr><tr><td>9-1</td><td>13:37:05.363</td><td>-29:59:52.46</td><td>559.8</td><td>103</td><td>103</td><td>0.40</td><td>0.27</td><td>1.43</td><td>0.45</td><td>6.79</td><td>1.02</td><td>17.0</td><td>3.2</td><td>&lt;0.22</td><td>0.55</td><td>0.08</td><td>0.95</td></tr><tr><td>10-1</td><td>13:37:05.055</td><td>-29:59:52.89</td><td>567.2</td><td>94</td><td>91</td><td>0.29</td><td>0.32</td><td>2.37</td><td>0.50</td><td>7.56</td><td>1.00</td><td>15.9</td><td>2.9</td><td>&lt;0.22</td><td>0.97</td><td>0.17</td><td>0.97</td></tr><tr><td>10-2</td><td>13:37:05.032</td><td>-29:59:53.00</td><td>565.7</td><td>61</td><td>49</td><td>0.22</td><td>0.27</td><td>3.12</td><td>0.27</td><td>4.17</td><td>1.01</td><td>7.7</td><td>2.4</td><td>&lt;0.22</td><td>1.30</td><td>0.20</td><td>0.97</td></tr><tr><td>12-1</td><td>13:37:05.195</td><td>-29:59:55.34</td><td>565.7</td><td>169</td><td>122</td><td>0.33</td><td>0.32</td><td>2.69</td><td>0.50</td><td>7.59</td><td>1.00</td><td>28.0</td><td>3.9</td><td>&lt;0.22</td><td>1.11</td><td>0.10</td><td>0.98</td></tr><tr><td>13-1</td><td>13:37:05.086</td><td>-29:59:55.76</td><td>567.6</td><td>391</td><td>246</td><td>0.45</td><td>0.49</td><td>3.17</td><td>0.73</td><td>11.15</td><td>1.01</td><td>74.5</td><td>6.1</td><td>&lt;0.22</td><td>1.33</td><td>0.10</td><td>0.95</td></tr><tr><td>14-1</td><td>13:37:04.982</td><td>-29:59:56.87</td><td>567.4</td><td>57</td><td>57</td><td>0.21</td><td>0.27</td><td>1.30</td><td>0.44</td><td>6.64</td><td>1.09</td><td>9.7</td><td>2.5</td><td>&lt;0.22</td><td>0.49</td><td>0.11</td><td>0.90</td></tr><tr><td>15-1</td><td>13:37:04.806</td><td>-29:59:59.42</td><td>570.0</td><td>40</td><td>40</td><td>0.16</td><td>0.26</td><td>1.08</td><td>0.42</td><td>6.38</td><td>1.42</td><td>7.4</td><td>2.7</td><td>&lt;0.22</td><td>0.38</td><td>0.23</td><td>0.68</td></tr><tr><th colspan="18">Level B Clumps (CY9UN)</th></tr><tr><td>4-1-1</td><td>13:37:05.010</td><td>-29:59:46.22</td><td>562.7</td><td>77</td><td>61</td><td>0.26</td><td>0.22</td><td>2.32</td><td>2.49</td><td>19.04</td><td>1.92</td><td>55.0</td><td>7.4</td><td>&lt;0.16</td><td>0.96</td><td>0.09</td><td>0.65</td></tr><tr><td>4-2-1</td><td>13:37:04.970</td><td>-29:59:46.30</td><td>564.9</td><td>39</td><td>39</td><td>0.18</td><td>0.22</td><td>1.62</td><td>1.29</td><td>9.83</td><td>1.94</td><td>21.3</td><td>5.3</td><td>&lt;0.16</td><td>0.64</td><td>0.04</td><td>0.65</td></tr><tr><td>4-2-2</td><td>13:37:04.968</td><td>-29:59:46.69</td><td>564.9</td><td>22</td><td>22</td><td>0.14</td><td>0.16</td><td>1.55</td><td>1.12</td><td>8.54</td><td>1.87</td><td>10.9</td><td>3.8</td><td>&lt;0.16</td><td>0.61</td><td>0.06</td><td>0.68</td></tr><tr><td>4-2-3</td><td>13:37:04.969</td><td>-29:59:45.99</td><td>567.4</td><td>34</td><td>25</td><td>0.16</td><td>0.14</td><td>5.42</td><td>1.15</td><td>8.76</td><td>2.02</td><td>17.3</td><td>5.1</td><td>&lt;0.16</td><td>2.29</td><td>0.14</td><td>0.63</td></tr><tr><td>9-1-1</td><td>13:37:05.357</td><td>-29:59:52.41</td><td>559.8</td><td>23</td><td>23</td><td>0.17</td><td>0.14</td><td>1.49</td><td>0.81</td><td>6.16</td><td>1.33</td><td>8.2</td><td>2.8</td><td>&lt;0.16</td><td>0.59</td><td>0.11</td><td>0.96</td></tr><tr><td>10-1-1</td><td>13:37:05.059</td><td>-29:59:52.83</td><td>567.4</td><td>17</td><td>17</td><td>0.13</td><td>0.13</td><td>1.34</td><td>0.85</td><td>6.49</td><td>1.31</td><td>6.3</td><td>2.3</td><td>&lt;0.16</td><td>0.52</td><td>0.14</td><td>0.98</td></tr><tr><td>12-1-1</td><td>13:37:05.194</td><td>-29:59:55.37</td><td>565.9</td><td>51</td><td>31</td><td>0.18</td><td>0.17</td><td>2.80</td><td>0.88</td><td>6.70</td><td>1.30</td><td>19.0</td><td>4.0</td><td>&lt;0.16</td><td>1.16</td><td>0.15</td><td>0.98</td></tr><tr><td>13-1-1</td><td>13:37:05.089</td><td>-29:59:55.76</td><td>567.7</td><td>97</td><td>80</td><td>0.27</td><td>0.41</td><td>2.63</td><td>1.00</td><td>7.63</td><td>1.32</td><td>36.9</td><td>5.7</td><td>&lt;0.16</td><td>1.10</td><td>0.16</td><td>0.96</td></tr><tr><td>14-1-1</td><td>13:37:04.983</td><td>-29:59:56.90</td><td>567.4</td><td>14</td><td>14</td><td>0.09</td><td>0.15</td><td>1.42</td><td>0.78</td><td>5.99</td><td>1.42</td><td>5.2</td><td>2.3</td><td>&lt;0.16</td><td>0.53</td><td>0.14</td><td>0.90</td></tr></tbody></table>

Notes. (1) Cloud ID. (2)(3)(4) Coordinate and recession velocity V. (5)(6) Numbers of pixels identified within cloud in the RA-Dec-Vel space $N_{pix}$ and in the RA-Dec projection $N_{xy}$. (7)(8)(9) FWHM widths in RA (=x), Dec (=y), and velocity (=v) directions, calculated from intensity-weighted dispersions. The FWHMv is calculated with the 0.85 km s⁻¹-channel cube, while the rest is with the 2.54 km s⁻¹-channel cube. (10)(11) Peak intensity $I_p$ and uncertainty $\Delta I_p$ after primary beam correction. $I_p$ is also given in the form of brightness temperature $T_p$. (12)(13) Total flux $F$ and uncertainty $\Delta F$ after primary beam correction, calculated by summing up the fluxes of all the pixels identified in cloud or clump. (14) Radius after beam dilution correction. Unlike $\sigma_r$, $R$ is supposed to trace the whole radial extent of object. (15)(16) Velocity dispersion and uncertainty after channel dilution correction. (17) Primary beam attenuation.
<figure><img src="image_8.png" /><figcaption>**Fig. A.1.** Cutout images of the clouds and clumps in CO(3–2) integrated intensity. The velocity range of each cloud is integrated. The primary beam correction is applied. From left to right: Different clouds, with their ID# in the top-left corners. From top to bottom: CY7NA, CY7+9NA, CY7+9UN, CY9NA, and CY9UN (see Table 1). The yellow circles have a diameter of 2" (~ 44 pc).</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. A.2.** Same as Fig. A.1, but in peak intensity.</figcaption></figure>
<table><caption>Table A.2. Cloud parameter measurements in CO(3-2) with the CY7NA cube.</caption><thead><tr><th rowspan="2">ID</th><th rowspan="2">RA<br>J2000</th><th rowspan="2">Dec<br>J2000</th><th rowspan="2">V<br>km s<sup>-1</sup></th><th rowspan="2">N<sub>pix</sub></th><th rowspan="2">N<sub>xy</sub></th><th rowspan="2">x<br>"</th><th rowspan="2">y<br>"</th><th colspan="3">FWHM</th><th rowspan="2">I<sub>p</sub><br>km s<sup>-1</sup></th><th rowspan="2">ΔI<sub>p</sub><br>kmJy/bm</th><th rowspan="2">F<br>kmJy km s<sup>-1</sup></th><th rowspan="2">ΔF<br>kmJy km s<sup>-1</sup></th><th rowspan="2">R<br>"</th><th rowspan="2">σ<sub>v</sub><br>km s<sup>-1</sup></th><th rowspan="2">Δσ<sub>v</sub><br>km s<sup>-1</sup></th><th rowspan="2">PB</th></tr><tr><th>v<br>km s<sup>-1</sup></th><th>mJy/bm</th><th>mJy</th></tr></thead><tbody><tr><td>1</td><td>13:37:05.044</td><td>-29:59:33.96</td><td>558.6</td><td>679</td><td>456</td><td>0.85</td><td>0.48</td><td>3.24</td><td>18.37</td><td>3.50</td><td>67.5</td><td>13.7</td><td>&lt;0.50</td><td>1.36</td><td>0.03</td><td>0.27</td></tr><tr><td>2</td><td>13:37:05.035</td><td>-29:59:36.15</td><td>555.3</td><td>364</td><td>281</td><td>0.55</td><td>0.49</td><td>3.50</td><td>13.10</td><td>2.51</td><td>24.0</td><td>6.8</td><td>&lt;0.50</td><td>1.46</td><td>0.10</td><td>0.40</td></tr><tr><td>3</td><td>13:37:04.814</td><td>-29:59:40.69</td><td>561.6</td><td>1471</td><td>945</td><td>0.94</td><td>0.86</td><td>3.55</td><td>17.71</td><td>1.90</td><td>100.2</td><td>10.8</td><td>&lt;0.50</td><td>1.49</td><td>0.06</td><td>0.50</td></tr><tr><td>4</td><td>13:37:04.984</td><td>-29:59:46.15</td><td>564.4</td><td>3973</td><td>1438</td><td>1.02</td><td>0.96</td><td>6.21</td><td>26.47</td><td>1.25</td><td>258.2</td><td>11.2</td><td>&lt;0.50</td><td>2.62</td><td>0.06</td><td>0.79</td></tr><tr><td>5</td><td>13:37:06.467</td><td>-29:59:49.13</td><td>562.2</td><td>618</td><td>574</td><td>0.82</td><td>0.83</td><td>1.96</td><td>8.22</td><td>1.12</td><td>22.7</td><td>3.8</td><td>&lt;0.50</td><td>0.80</td><td>0.12</td><td>0.92</td></tr><tr><td>6</td><td>13:37:06.179</td><td>-29:59:50.41</td><td>552.2</td><td>508</td><td>508</td><td>1.00</td><td>0.54</td><td>1.52</td><td>6.37</td><td>1.07</td><td>16.5</td><td>3.2</td><td>&lt;0.50</td><td>0.60</td><td>0.04</td><td>0.98</td></tr><tr><td>7</td><td>13:37:05.337</td><td>-29:59:50.60</td><td>559.8</td><td>192</td><td>192</td><td>0.50</td><td>0.40</td><td>1.38</td><td>5.53</td><td>1.08</td><td>5.8</td><td>2.0</td><td>&lt;0.50</td><td>0.52</td><td>0.15</td><td>0.97</td></tr><tr><td>8</td><td>13:37:05.054</td><td>-29:59:51.21</td><td>565.9</td><td>856</td><td>684</td><td>1.07</td><td>0.58</td><td>3.23</td><td>7.26</td><td>1.17</td><td>30.5</td><td>4.7</td><td>&lt;0.50</td><td>1.35</td><td>0.08</td><td>0.88</td></tr><tr><td>9</td><td>13:37:05.352</td><td>-29:59:52.77</td><td>561.2</td><td>904</td><td>698</td><td>0.81</td><td>0.84</td><td>2.82</td><td>7.80</td><td>1.08</td><td>30.2</td><td>4.4</td><td>&lt;0.50</td><td>1.17</td><td>0.08</td><td>0.97</td></tr><tr><td>10</td><td>13:37:05.040</td><td>-29:59:53.04</td><td>567.5</td><td>1266</td><td>935</td><td>0.80</td><td>1.04</td><td>3.36</td><td>11.64</td><td>1.17</td><td>49.6</td><td>5.7</td><td>&lt;0.50</td><td>1.41</td><td>0.10</td><td>0.87</td></tr><tr><td>11</td><td>13:37:05.917</td><td>-29:59:55.36</td><td>564.4</td><td>752</td><td>572</td><td>0.71</td><td>0.71</td><td>2.21</td><td>9.13</td><td>1.06</td><td>27.1</td><td>3.9</td><td>&lt;0.50</td><td>0.91</td><td>0.13</td><td>0.99</td></tr><tr><td>12</td><td>13:37:05.186</td><td>-29:59:55.41</td><td>565.5</td><td>700</td><td>494</td><td>0.74</td><td>0.58</td><td>2.60</td><td>9.74</td><td>1.11</td><td>27.4</td><td>4.0</td><td>&lt;0.50</td><td>1.07</td><td>0.11</td><td>0.93</td></tr><tr><td>13</td><td>13:37:05.084</td><td>-29:59:55.69</td><td>567.6</td><td>1607</td><td>748</td><td>0.82</td><td>0.69</td><td>3.64</td><td>17.57</td><td>1.16</td><td>78.6</td><td>6.4</td><td>&lt;0.50</td><td>1.52</td><td>0.08</td><td>0.88</td></tr><tr><td>14</td><td>13:37:04.978</td><td>-29:59:56.98</td><td>567.8</td><td>388</td><td>313</td><td>0.59</td><td>0.49</td><td>2.34</td><td>6.81</td><td>1.22</td><td>13.3</td><td>3.3</td><td>&lt;0.50</td><td>0.96</td><td>0.17</td><td>0.85</td></tr><tr><td>15</td><td>13:37:04.820</td><td>-29:59:59.49</td><td>570.4</td><td>763</td><td>607</td><td>0.75</td><td>0.78</td><td>2.32</td><td>12.09</td><td>1.34</td><td>35.8</td><td>5.2</td><td>&lt;0.50</td><td>0.94</td><td>0.11</td><td>0.75</td></tr><tr><td>16</td><td>13:37:06.684</td><td>-29:59:59.73</td><td>564.9</td><td>513</td><td>387</td><td>0.54</td><td>0.69</td><td>3.79</td><td>6.18</td><td>1.22</td><td>16.6</td><td>3.8</td><td>&lt;0.50</td><td>1.59</td><td>0.17</td><td>0.85</td></tr><tr><td>17</td><td>13:37:04.738</td><td>-30:00:01.06</td><td>569.6</td><td>754</td><td>626</td><td>0.74</td><td>0.83</td><td>2.06</td><td>9.76</td><td>1.42</td><td>35.6</td><td>5.5</td><td>&lt;0.50</td><td>0.85</td><td>0.10</td><td>0.70</td></tr><tr><td>18</td><td>13:37:06.407</td><td>-30:00:01.43</td><td>564.0</td><td>465</td><td>329</td><td>0.54</td><td>0.57</td><td>2.71</td><td>7.40</td><td>1.09</td><td>15.0</td><td>3.2</td><td>&lt;0.50</td><td>1.12</td><td>0.14</td><td>0.96</td></tr><tr><td>19</td><td>13:37:06.006</td><td>-30:00:02.05</td><td>562.3</td><td>410</td><td>410</td><td>0.62</td><td>0.74</td><td>1.58</td><td>6.00</td><td>1.06</td><td>13.1</td><td>2.9</td><td>&lt;0.50</td><td>0.62</td><td>0.04</td><td>1.00</td></tr><tr><td>20</td><td>13:37:04.848</td><td>-30:00:03.57</td><td>569.0</td><td>681</td><td>544</td><td>0.88</td><td>0.62</td><td>2.68</td><td>7.43</td><td>1.30</td><td>24.3</td><td>4.6</td><td>&lt;0.50</td><td>1.12</td><td>0.10</td><td>0.79</td></tr><tr><td>21</td><td>13:37:05.567</td><td>-30:00:06.58</td><td>567.4</td><td>230</td><td>230</td><td>0.64</td><td>0.36</td><td>1.54</td><td>6.07</td><td>1.07</td><td>7.3</td><td>2.2</td><td>&lt;0.50</td><td>0.60</td><td>0.06</td><td>0.99</td></tr><tr><td>22</td><td>13:37:05.770</td><td>-30:00:11.22</td><td>570.5</td><td>547</td><td>533</td><td>1.03</td><td>0.58</td><td>2.65</td><td>8.34</td><td>1.16</td><td>19.7</td><td>3.7</td><td>&lt;0.50</td><td>1.10</td><td>0.12</td><td>0.89</td></tr><tr><td>23</td><td>13:37:05.209</td><td>-30:00:15.15</td><td>567.1</td><td>230</td><td>198</td><td>0.49</td><td>0.41</td><td>2.20</td><td>8.05</td><td>1.56</td><td>9.8</td><td>3.3</td><td>&lt;0.50</td><td>0.90</td><td>0.18</td><td>0.65</td></tr></tbody></table>

Notes. (1) Cloud ID. (2)(3)(4) Coordinate and recession velocity V. (5)(6) Numbers of pixels identified within cloud in the RA-Dec-Vel space $N_{pix}$ and in the RA-Dec projection $N_{xy}$. (7)(8)(9) FWHM widths in RA (=x), Dec (=y), and velocity (=v) directions, calculated from intensity-weighted dispersions. The FWHMv is calculated with the 0.85 km s⁻¹-channel cube, while the rest is with the 2.54 km s⁻¹-channel cube. (10)(11) Peak intensity $I_p$ and uncertainty $\Delta I_p$ after primary beam correction. (12)(13) Total flux F and uncertainty $\Delta F$ after primary beam correction, calculated by summing up the fluxes of all the pixels identified in cloud. (14) Radius after beam dilution correction. Unlike $\sigma_r$, R is supposed to trace the whole radial extent of object. (15)(16) Velocity dispersion and uncertainty after channel dilution correction. (17) Primary beam attenuation.
**Appendix B: The stacked cutout images**

A stacking analysis does not provide a significant CO(2–1) detection in this case. The clouds have a large range of flux (Table A.1), and only the brightest ones contribute to the stacked flux, while all the cloud images contribute to the stacked noise.

Figure B.1 shows the stacked CO(3–2) and CO(2–1) images, that is, the sum of the integrated intensity cutout images of the 11 clouds (see Figs. 2 and A.1 and Table A.1). The same velocity ranges are integrated for CO(3–2) and CO(2–1). The beam sizes are different between the two panels (Table 1). If we assume that the (stacked) object is not resolved with either beam, the total flux of the object is the values displayed in these images.

Panel (a) shows the peak CO(3–2) flux of $F_{32} = 550$ mJy km s$^{-1}$. The noise of this image is 16 mJy km s$^{-1}$. The corresponding CO(2–1) flux, predicted under the thermalized condition, is $F_{21} = (4/9\alpha_{th})F_{32} = 269$ mJy km s$^{-1}$ ($\alpha_{th} \approx 1.1$ at T=30 K; see Sect. 6 and Eq. 5).

Panel (b) shows a tempting emission peak, but there is also a similarly deep negative just north of it. Indeed, the peak flux is $F_{21} = 263$ mJy km s$^{-1}$ in CO(2–1), which is apparently consistent with the prediction under the thermalized condition. However, the noise in this image is 140 mJy km s$^{-1}$. Thus, it is only at a $1.9\sigma$ significance. The consistency is supportive of the model (Fig. 1) but is not conclusive with this CO(2–1) data.

<figure><img src="image_10.png" /><figcaption>**Fig. B.1.** Stacked (summed) images of the 11 clouds in the Cy9 FoV (see Fig. 2 and Table A.1). (a) CO(3–2) image with CY7NA (with an RMS noise of 16 mJy/beam km s$^{-1}$). (b) CO(2–1) image with CY3NA (140 mJy/beam km s$^{-1}$). The white circles have a diameter of 2". Note that the beam sizes are different in the two panels (see Table 1). Assuming that the emission is not resolved, the fluxes in these images indicate the total fluxes of the (stacked) object.</figcaption></figure>