{'fraction_non_alphanumeric': 0.05013275315288738, 'fraction_numerical': 0.04104284976768198, 'mean_word_length': 4.072196764238287, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 7, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
First Detection of the Molecular Cloud Population in the Extended Ultraviolet (XUV) Disk of M83

We report a CO(3-2) detection of 23 molecular clouds in the extended ultraviolet (XUV) disk of the spiral galaxy M83 with ALMA. The observed 1kpc^2 region is at about 1.24 times the optical radius (R25) of the disk, where CO(2-1) was previously not detected. The detection and non-detection, as well as the level of star formation (SF) activity in the region, can be explained consistently if the clouds have the mass distribution common among Galactic clouds, such as Orion A -- with star-forming dense clumps embedded in thick layers of bulk molecular gas, but in a low-metallicity regime where their outer layers are CO-deficient and CO-dark. The cloud and clump masses, estimated from CO(3-2), range from 8.2x10^2 to 2.3x10^4 Msun and from 2.7x10^2 to 7.5x10^3 Msun, respectively. The most massive clouds appear similar to Orion A in star formation activity as well as in mass, as expected if the cloud mass structure is universal. The overall low SF activity in the XUV disk could be due to the relative shortage of gas in the molecular phase. The clouds are distributed like chains up to 600 pc (or longer) in length, suggesting that the trigger of cloud formation is on large scales. The universal cloud mass structure also justifies the use of high-J CO transitions to trace the total gas mass of clouds, or galaxies, even in the high-z universe. This study is the first demonstration that CO(3-2) is an efficient tracer of molecular clouds even in low-metallicity environments.


INTRODUCTION
The Galaxy Evolution Explorer (GALEX) satellite found massive star formation (SF) in the far outskirts of galactic disks (Gil de Paz et al. 2005;Thilker et al. 2005). Bright ultraviolet (UV) sources are distributed beyond the optical radius R 25 . They reveal abundant and recent (< 100 Myr) SF (Gil de Paz et al. 2005;Thilker et al. 2005;Lemonias et al. 2011, see Figure 1 for the example of M83). These extended UV disks, dubbed XUV disks, are fairly common among local disk galaxies (Thilker et al. 2007). They offer an opportunity to study SF in extreme conditions, in particular, at a low average gas density and molecular fraction. XUV disks often exhibit lower Hα-to-far UV (FUV) flux ratios compared to the optical disk, which likely demonstrates the importance of stochastic sampling of the initial mass function (IMF) and/or bursty SF histories in these regions (e.g., Alberts et al. 2011;Koda et al. 2012a;Watson et al. 2016).
Molecular clouds host virtually all SF within the optical disks of local galaxies. It is crucial to study whether the same is true in XUV disks. The Milky Way's outer disk hosts small molecular clouds (10 2−4 M ; Sun et al. 2015Sun et al. , 2017, but unlike XUV disks there is not much SF above late B-type stars (Izumi et al. 2017). Numerous efforts have been made to detect CO emission in XUV disks, however, they have rarely succeeded (Watson & Koda 2017, for review including unpublished efforts that resulted in non-detections). Only four galaxies with XUV disks permitted CO detection at a few positions in their outskirts: NGC 4414 (Braine & Herpin 2004), NGC 6946 (Braine et al. 2007), M33 (Braine et al. 2010), and M63 (Dessauges-Zavadsky et al. 2014). Most detections could not reveal individual clouds in the vicinity of UV emission due to low sensitivity and resolution ( 10 5−6 M and 300-500 pc). Only M33, the closest to the Milky Way, allowed detection down to 4 × 10 4 M in two regions beyond the R 25 radius, and only one of them is identified as a single cloud (Braine et al. 2012). The rarity of the CO detection is at odds with the abundance of the UV sources across the XUV disks.
M83 is at a distance of D = 4.5 Mpc (Thim et al. 2003) and is one of the nearest XUV disks. Previous ALMA observations of CO(2-1) in the XUV disk of M83 resulted in non-detection (Bicalho et al. 2019) even at a high sensitivity of 2.2 × 10 4 M (3σ; calculated with the low CO-to-H 2 conversion factor in the Large Magellanic Cloud (LMC) from Fukui et al. 2008). The observed field is a relatively large area of ∼ 1.5 × 3 (2 × 4 kpc 2 ) at a galactic radius of r gal ∼ 1.24R 25 with R 25 = 6.44 (= D 25 /2) from the optical center at (RA, DEC) J2000 = (13:37:00.4, -29:52:04.1) (de Vaucouleurs et al. 1991). We revisit a smaller area of ∼ 0.75 × 0.85 within the CO(2-1) field at a higher mass sensitivity by a factor of about 10 in the CO(3-2) line emission, as well as in the Band 7 dust continuum emission.
M83 presents a prototype XUV disk (Thilker et al. 2005). Deep Hα imaging with the Subaru telescope revealed HII regions over the XUV disk, including some in the region of this study (Koda et al. 2012a). The metallicity is about constant across the XUV disk and is as low as 12 + log(O/H) ∼ 8.2 ("preferred abundances" by Bresolin et al. 2009, from their analyses with multiple metallicity indicators) 1 . This is ∼ 0.4Z by adopting the solar abundance of 12 + log(O/H) =8.66 (Asplund et al. 2005), which is similar to that of the LMC (∼ 0.3-0.5Z , Westerlund 1997). With Spitzer infrared images, Dong et al. (2008) concluded that the SF has been ongoing for at least 1 Gyr. The SF efficiency derived with atomic hydrogen (HI) data indicates a gas consumption timescale much longer than the Hubble time (Bigiel et al. 2010). This extremely low SF efficiency could be due to physical processes either on large scales (e.g., low HI-to-H 2 phase transition), or on small scales (e.g., low SF rate within molecular clouds). Observations and detection of molecular gas would improve our understanding of SF in this extreme environment.

ALMA Observations
We observed a region of ∼ 0.75 × 0.85 (∼ 0.98 × 1.11 kpc 2 ; approximately ∼ 1 kpc 2 ) in the outskirts of M83 with ALMA's 12m array in CO(3-2) and Band 7 dust continuum emission. Figure 1 shows the target region: the left panel displays the global context of the galaxy, and the right panels show the observed region with magenta boxes (smaller boxes). The outer (cyan) contour in the left panel shows the extent of the HI disk, which also encloses the XUV disk and numerous UV blobs (the color image). The red boxes in all panels show the region of the previous CO(2-1) and Band 6 continuum emission data (Bicalho et al. 2019). Our region is selected around the maximum HI 21 cm intensity position within the CO(2-1) observation. It includes bright UV peaks located at the galactocentric distance of ∼ 1.24R 25 (∼ 8.0 , ∼ 10.4 kpc), and is outside the traditional optical disk (green contour).
The CO(3-2) line (ν CO(3−2) = 345.79599 GHz) and continuum emission in Band 7 were observed with ALMA (project # 2017.1.00065.S). It mosaicked the target region with 17 pointing positions around the map center at (RA, DEC) J2000 = (13:37:05.8, -29:59:57.4). One spectral window (SPW) was configured for the line emission with a bandwidth of 1.875 GHz (1626 km s −1 ) with 1920 channels of 976.6 kHz width (0.8466 km s −1 ).  The inner (green) and outer (cyan) contours are the edge of the optical disk and the extent of HI gas at an HI surface density of 1.5 × 10 20 cm −2 (see Koda et al. 2012a). The red box is the 1.5 × 3 (2.0 × 3.9 kpc 2 ) field from the ALMA CO(J=2-1) observations (Bicalho et al. 2019). The three panels on the right are zoom-in to the CO(2-1) field: (b) the GALEX color composite, as in Panel a, (c) the continuum-subtracted Hα image from the Subaru telescope, and (d) the HI 21cm emission image from Walter et al. (2008). The area of our CO(3-2) observations is indicated by the inner magenta box of ∼ 0.75 × 0.85 (∼ 0.98 × 1.11 kpc 2 ). The bottom two panels are zoom-in to the CO(3-2) field: (e) The Hα image with visually-identified HII regions (circles The other three SPWs were configured to cover different sky frequencies for the continuum emission with a bandwidth of 2 GHz with 128 channels per each SPW. The central frequency of the continuum emission is 349.498 GHz (hereafter, 349 GHz continuum emission). The full mosaic observations were repeated six times, resulting in six execution blocks. The final uv-coverage at the CO(3-2) frequency extends over the angular scales of 0.57-13.7 (∼12-299 pc). The CO(2-1) line (ν CO(2−1) = 230.53800 GHz) and continuum emission in Band 6 were taken in project # 2013.1.00861.S and were re-reduced for consistency within this study. These data have been studied by Bicalho et al. (2019), and the details of the observations are found there. The rectangle 1.5 × 3 region includes the area of the CO(3-2) and Band 7 continuum data. The bandwidth and channel width are 1.875 GHz (2, 438 km s −1 ) and 1.953 MHz (2.540 km s −1 ) for the CO(2-1) emission. The central frequency for the continuum emission is 224.516 GHz (hereafter, 225 GHz continuum emission). The covered angular scales are 0.49-10.7 (∼11-233 pc) at the CO(2-1) frequency.
These data are calibrated using the data reduction scripts provided by the ALMA observatory using the Common Astronomy Software Application (CASA: Mc-Mullin et al. 2007). The amplitude and phase of bandpass and gain calibrators are confirmed to be flat over time and frequency after the calibrations.
For the CO(3-2) and CO(2-1) lines we generate data cubes with the tclean task in a standard way. We use a cell size of 0.15 and channel width of 2.54 km s −1 (i.e., the channel width of the CO(2-1) data) over a velocity range of ∼400 to 700 km s −1 . While we use the 2.54 km s −1 cube throughout this study unless otherwise specified, we also make a separate CO(3-2) data cube with a channel width of 1.00 km s −1 . This cube is used solely for a measurement of velocity dispersions of molecular clouds. The 349 and 225 GHz continuum data are imaged with the multi frequency synthesis (MFS) mode of tclean. The CO(2-1) and 225 GHz continuum data are regridded to the image format of the CO(3-2) and 349 GHz continuum data for comparisons. The parameters of the data are listed in Table 1, including the total band width for continuum, the channel width for lines, the cell/pixel size, beam size (its major and minor axis diameters, b maj , b min , and position angle PA), and root-mean-square (RMS) noise.

Supplementary Data
We obtained the HI 21cm data with the Green Bank Telescope (GBT; project 20A-432) and combined it with the archival Karl G. Jansky Very Large Array (VLA) data (projects 13B-194 and 14B-192). We reduced the data in the standard way with GBTIDL and CASA, respectively. The GBT and VLA data are combined after imaging with CASA's feather task. The beam size of the GBT+VLA data is 15 . The details of the observations and data reduction will be presented in a separate publication.
The Hα narrow-band image was obtained from narrow (NA656) and broad (Rc) band photometry with the Subaru Prime Focus Camera (Suprime-Cam) on the Subaru telescope. The seeing was about 1 . The observations and data reduction were similar to those presented in Koda et al. (2012a), who used the privatelyowned old Hα filter (NA659). M83 was re-observed with the observatory-owned filter in the context of a separate, larger Hα survey of galaxies for consistency. The photometric zero points were calibrated using the Pan-STARRS1 Surveys (PS1) Data Release 1 catalog (Chambers et al. 2016). We followed the procedure described in Yagi et al. (2013) and converted from the PS1 to the Suprime-Cam photometric system using stellar photometry. The flux errors are about 10% or less, and the old and new measurements are consistent within the errors. Note that discussions in this paper require a flux accuracy of only a factor of ∼ 2. The details of the Hα data will be presented in a separate paper.

HII Regions
A classification of HII regions is not within the scope of this paper, but for the purpose of rough evaluation of SF activity in the observed region, we visually identified significant Hα peaks in the Subaru H α image. Table 2 lists the ID, coordinates, aperture radius r ap , and Hα luminosity L Hα . Figures 1e and 1f show their locations. The L Hα ranges from 1 × 10 35 to 3 × 10 37 erg s −1 . 2 We do not apply an extinction correction since it is expected to be small (Gil de Paz et al. 2007a) and is negligible for the discussion here which requires only order of magnitude estimations. Our HII regions (B, F, J) correspond to (#09, #08, #07) in Gil de Paz et al. (2007a) and (#39, #41, #40) in Bresolin et al. (2009). The total Hα luminosity in our 1 kpc 2 region is ∼ 6.8 × 10 37 erg s −1 .  Note-LH α is a Hα luminosity within an aperture of the radius rap which is set by visual judgement (see Figure 1e).
The identified HII regions contain about 80% of the total luminosity. As a reference, a model calculation by Sternberg et al. (2003) suggests that a single B0-type star can produce an HII region of L Hα = 1.4 × 10 36 erg s −1 assuming the Case B recombination and electron temperature of 10 4 K. The HII regions brighter than this likely have Otype stars, but they have to be relatively late O-type, as an HII region around a single O6-type star can already be as bright as 2.7 × 10 37 erg s −1 , similar to the brightest in our XUV region. As we see below, the observed XUV region has SF activity similar to, or slightly less than, that of the Orion Nebula, a relatively minor OB star forming region in the Milky Way.

The Orion A Molecular Cloud in the MW
The measured L Hα are similar to those of the Galactic HII regions, M42 and M43. Both of them are in the Orion Nebula, which is located in the OMC-1 dense gas clump in the Orion A molecular cloud 3 . They show L Hα = 7.1 × 10 36 and 3.5 × 10 35 erg s −1 , respectively (Scoville et al. 2001, from radio continuum observations by Schraml & Mezger 1969). M42 is ionized predominantly by a couple of late O-type stars (O7 & O9) with small contributions from a few B-type stars (Hillenbrand 1997), and M43 is excited by a single B0.5-type star (Simón-Díaz et al. 2011). Hence, individual HII regions in the XUV disk have similar SF activity to that of the Orion Nebula. There are also numerous low-mass stars forming across the Orion A cloud (Großschedl et al. 2019), while their counterparts are not easily detected at the distance of M83.
The Orion A molecular cloud is one of the nearest clouds with massive SF and arguably the most studied high-mass star forming region. Its internal structure, i.e, a dense gas clump surrounded by a large volume of bulk molecular gas, is common among Galactic molecular clouds; not only among high-mass star forming clouds (e.g., W49 and W51; Nakamura et al. 1984;Watanabe et al. 2017;Barnes et al. 2020), but also among low-mass star forming clouds (e.g., Perseus, Serpens, and Ophiuchus; Enoch et al. 2008). Therefore, we use Orion A as our reference cloud. If the Orion A cloud (and Orion Nebula) were in the XUV disk of M83, it would appear as one of brightest UV and H α blobs. As a gauge of the activity, compared to other Galactic OB associations, Orion A is "not particularly impressive as a massive star-forming region" (Hillenbrand 1997).  Table 3. Each circle has a diameter of 2 . The beam size of 0.96 × 0.82 (PA=−84.7 • ) is shown on the bottom-left corner, corresponding to 21 × 18 pc 2 at 4.5 Mpc.
The most recent estimation of bulk molecular gas mass in Orion A is M cloud = 4.0×10 4 M , and its dense clump, OMC-1, has a clump mass of M clump = 7.4 × 10 3 M (Nakamura et al. 2019). Thus, the OMC-1 clump contains about 20% of the cloud mass. We adopt these masses for our discusson, but note that an often-quoted older value of M cloud is about two times larger and is M cloud = 1.0 × 10 5 M after correction for the distance of 414 pc (Dame et al. 1987).
The star formation rate (SFR) within Orion A is estimated with several tracer emissions by Pabst et al. (2021) (their Table 5). For comparisons to the XUV disk later, the most relevant is SFR∼ 7.7 × 10 −5 M yr −1 , based on Hα and 24µm emission. The gas consumption timescale by the SF is M cloud /SFR ∼ 520 Myr. A caveat is that the SFR estimation in a single cloud is quite uncertain; with the tracers closely related to highmass SF (Hα, 24µm, far and total infrared emissions), Pabst et al. (2021) found a range of SFR=1.6 × 10 −5 to 1.0 × 10 −4 M yr −1 . Keeping this uncertainty in mind, we will compare SFR in the XUV disk from Hα on an assumption of negligible extinction, and that in Orion A from Hα and 24µm to take extinction into account. Both of these trace the amount of UV photons from recently-formed massive stars.

Cloud Identification and Parameters
We identify 23 clouds in the CO(3-2) data in the following way and list them in Table 3. In the data cube, we first identify pixels with > 5σ peaks and find their envelopes by expanding their volumes to adjacent pixels down to 3σ significance. In some cases the edges of the 3σ envelopes of two peaks (two clouds) are connected. We use the CLUMPFIND algorithm (Williams et al. 1994) to split them -the partition of the two clouds is determined simply by assigning each pixel to the nearest of the two peaks.
We visually inspect the results and manually merge them back to a single cloud when the intensity-weighted positions of the split clouds are closer than the beam size (0.96 ) and one velocity pixel (2.54 km s −1 ). Clouds 4, 8, and 20 are merged back. Cloud 4 might consist of two clouds based on visual inspection, but we strictly apply the same merge criteria for all cases. Figure 2d shows the locations of the identified clouds with 2"-diameter circles (about double the beam size). Figure 4 shows the integrated intensity and peak intensity maps of individual clouds. Figure 5 shows the spectra of the clouds. The previous CO(2-1) study did not detect these clouds (Bicalho et al. 2019). In fact, none of them show > 3σ emission in the CO(2-1) cube (see Section 5.3.1). Table 3 lists the cloud properties. The first three columns are their intensity-weighted coordinates in RA, DEC, and recession velocity V . The FWHMx and FWHMy are the full width at half maximums in the RA (=x) and DEC (=y) directions and are calculated from intensity-weighted dispersions σ X (where X =x, y). The σ X s are calculated without any model profile, but we convert them to FWHM by assuming a Gaussian distribution (i.e., FWHM = √ 8 ln 2σ X ), so that they can be directly compared to the beam size in FWHM. We note the major axis of the beam is roughly in the RA (x) direction, and the minor axis is in the DEC (y) direction. N xy is the number of spatial pixels with > 3σ emission in the sky projection. The beam area amounts to N xy = 40 in FWHM. σ v is the intensity-weighted velocity dispersion calculated with the 1 km s −1 cube. Only the pixels with > 3σ emission in the 1 km s −1 cube (red in Figure  5) are included for this calculation. When the number of > 3σ pixels is less than four, we place only an upper limit in σ v , assuming the upper limit of the FWHM width to be 4 km s −1 (i.e., σ v = 1.7 km s −1 ). We do not apply deconvolution for the 1 km s −1 channel width since it is negligible (i.e., the second moment of the 1 km s −1 boxcar comes to only σ = 1/ √ 12 ∼ 0.29 km s −1 , and its quadratic subtraction from σ v makes only a negligible difference). The peak intensity, S peak CO(3−2) , and total flux, S CO(3−2) dv, are also listed in the table. The error in total flux is calculated by shifting the cloud volume within the 3D data cube in random directions 300 times and by calculating the standard deviation among the shifted volumes.

Sizes
Some clouds apparently show internal structures in Figure 4 and are on the verge of being spatially resolved. On the other hand, most have FWHM sizes similar to, or smaller than, the beam size. Thus, all clouds, except Cloud 4 and 8, are spatially unresolved or only very marginally resolved at the resolution of 0.96 (21 pc). Clouds 3, 5, and 22 show FWHMx and/or FH-WMy larger than the beam size (Table 3), but they suffer from low level emission (Figure 4). For example, Cloud 3 has a large peculiar extension at around 3-4σ significance, which increases FWHMx and FHWMy, but appears spurious. The extended components of Cloud 5 and 22 are also at low level (3-4σ). We also note that FWHMy is measured in about the direction of the beam minor axis. It is often smaller than the minor axis size of 0.82 , again because the measurements suffer from the low level emission.
Cloud 4 and 8 show two peaks in Figure 4, which are merged into a single cloud in our cloud identification process since they are closer than the resolution. We consider them as slightly resolved single clouds with a caveat that they could be a blend of two unresolved clouds.
Even in the marginally-resolved cases, measuring cloud sizes is difficult. The cloud radius R may be estimated with the equation R = 1.92σ r (Solomon et al. 1987) using the beam-deconvolved average standard deviation in radius Since the measured size is sometimes smaller than the beam size (Table 3), we consider σ r and R measurable only when the measured size is at least 20% larger than the beam size (assuming that FWHMx and FWHMy are measured with 20% accuracy). This sets the minimum measurable R to be R = 0.54 b maj b min = 0.48 (10.5 pc). Otherwise, we have only an upper limit of 10.5 pc. In this way, only one cloud with the spurious extension (Cloud 3; Figure 4) gives a measurable R, while, as discussed, it suffers from the error. Therefore, the CO(3-2)-emitting parts in the detected clouds (their diameters) are smaller than the resolution of 21 pc. As a reference, the Orion A cloud has a whole extent of ∼ 20 pc when it is observed in the CO(1-0) or CO(2-1) emission which trace the bulk molecular gas (Sakamoto et al. 1994;Nakamura et al. 2019;Kong et al. 2018). Its star-forming dense clump OMC-1 occupies only a small area (∼ 2 pc) at the heart of the cloud and is bright in CO(3-2) (Ikeda et al. 1999). It is possible that our detected CO(3-2) peaks are mainly from dense clumps within molecular clouds, and hence, are unresolved.

Velocity Widths
The velocity widths of most clouds are resolved in the 1 km s −1 cube (see σ v in Table 3). Cloud 3 has a relatively large σ v due to the spurious extension. Cloud 4 has a wide width of σ v = 6.1 km s −1 with a caveat that it could be a blend of two clouds.
Again using Orion A as a reference, the velocity width of the bulk gas, measured in CO(1-0) and averaged over the whole cloud, is about 5 km s −1 at 50% of the peak brightness (from Figure 5 of Nakamura et al. 2019). The full CO(3-2) velocity width over the OMC-1 clump (size of about 2 pc) is not found in the literature; toward Orion KL (the brightest region), it is about 7 km s −1 at 50% of the peak brightness in a 134 (∼ 0.3 pc) beam (from Figure 15 of Masui et al. 2021). Hence, the dense part has a broader linewidth than the entire cloud. These are similar to the velocity widths of Cloud 4 measured in CO(3-2). Figure 6 plots σ v against the radius R (mostly the upper limits; see Section 4.2). It also shows clouds in some other galaxies (Solomon et al. 1987;Bolatto et al. 2008;Wong et al. 2011;Rubio et al. 2015), with notes that they are measured in CO(1-0) or CO(2-1) while our XUV clouds are in CO(3-2) and that the different CO transitions may trace different parts of clouds (Section 5). For their measured σ v , our XUV clouds generally have smaller (or similar) sizes compared to the clouds in the other galaxies. This makes sense if the CO(3-2) emission is primarily tracing dense clumps within the clouds, as seen in Orion A.

Distribution
The 23 clouds are distributed over a large area across the 1 kpc 2 field of view (Figures 2d and 3). Their overall distribution shows some coherent (not random) structures, which appear like chains of clouds. These chains extend over 600 pc in length, or potentially longer beyond the field of view. The cloud distribution also shows a large-scale velocity gradient in the channel maps (Figure 3). From the low to high velocity channels, it spatially shifts from the northeast to the southwest side. Figure 7 shows the cloud locations on the Hα, UV, and HI images. Cloud 4, 7, 22, 23, and possibly 5, 16, 19, are directly associated with HII regions B, E, J, K, C, F, and H, respectively. Among these, B, J, K, and C are four of the most prominent HII regions in the field. Cloud 4 is the brightest in CO(3-2), while, interestingly, the other three clouds are in the bottom 2/3 of the clouds in total CO(3-2) flux. All the other clouds are in the proximity of, but outside, the HII regions. The variations could be due to an evolutionary sequence of molecular clouds, or to the stochasticity of massive SF in this low SF environment (i.e., when the total stellar mass produced in a region is low, clouds may not always produce stars over the full mass spectrum. The stellar initial mass function may be filled stochastically, but not completely; Koda et al. 2012a). With respect to the atomic gas distribution, most clouds and HII regions are around, but not on, the peak of the HI emission.

Molecular Clouds in Low Metallicity
We are using the Orion A molecular cloud as a fiducial reference cloud given the similarity of its SF activity to the (brightest) clouds in the XUV disk. Orion A has an internal mass distribution common among other large and small Galactic molecular clouds, having small starforming dense clump(s) embedded deeply in thick layers of bulk molecular gas (Nakamura et al. 1984;Enoch et al. 2008;Watanabe et al. 2017;Barnes et al. 2020). Figure 8a schematically illustrates this mass distribution. In terms of CO emissions, CO(1-0) and CO(2-1) can be excited easily at an average density and temperature in molecular clouds (∼ 10 2 cm −3 , ∼ 10 K; Scoville & Sanders 1987), and hence, the bulk gas contributes significantly to their total CO(1-0) and CO(2-1) luminosities (Kong et al. 2018;Nakamura et al. 2019). By having the J = 3 level temperature of E J=3 /k ∼ 33.2 K and critical density of > 10 3 cm −3 , the CO(3-2) emission requires higher density and temperature than the averages for excitation and is radiated predominantly from the dense clumps (Ikeda et al. 1999). Once it is excited in thermalized gas, CO(3-2) can become brighter than CO(1-0) and CO(2-1) in total flux.
We hypothesize that the clouds in the XUV disk share the common clump-envelope mass distribution of the Galactic molecular clouds. The common mass structure is expected for gravitationally-bound clouds as it is determined predominantly by the internal physics rather than the environment.
Even when the mass distribution is the same, the chemical structure and appearance of the clouds in the CO emissions depend on the metal abundance (e.g., Maloney & Black 1988;van Dishoeck & Black 1988;Israel 1997;Pak et al. 1998;Wolfire et al. 2010). In high metallicity (∼ 1Z ; Figure 8a), the bulk gas in clouds can be detected in the low excitation transitions (J =1-0 or 2-1) of CO -the brightest and second-most-abundant molecule after H 2 . In low metallicity (< 1Z ; ; Figure  8b), the dust extinction is low, and the ambient stellar radiation can penetrate deeper into the clouds. This radiation dissociates more CO in their outer layers through CO line absorptions, while the abundant H 2 molecules can be self-shielded. This differential photo-dissociation makes the cloud outer layers CO-deficient and CO-dark while still H 2 -abundant (Figure 8b).
The ambient radiation can come from the UVemitting sources across the 1 kpc 2 region ( Figure 1f) and be potentially as strong as that in the solar neighborhood. There are several Orion-class SF regions within the region (Table 2, Figure 7). In the solar neighborhood, the Orion cloud is the closest massive SF region at a distance of 414 pc.
In addition to this photo-dissociation model for the cloud envelopes, we hypothesize that the dense clumps reside in the self-shielded, CO-abundant part of the clouds. The clouds in the XUV disk may suffer from this differential photo-dissociation due to the low metallicity. With this hypothesis, we can consistently explain the detection of CO(3-2), which is mainly from the dense clumps; non-detection of CO(2-1); as well as Hα luminosities, tracing SF activity.

Dense Clumps at Hearts of Clouds
In our hypothesis, the dense clumps can be shielded from photo-dissociation and radiate CO emission even in the low-metallicity environment (Figure 8b). The fractions of the CO-abundant part and dense clumps over entire clouds can be crudely estimated from a combination of previous measurements of the mass-to-light ratio M cloud /L CO -between the cloud's total mass M cloud and CO luminosity L CO . The M cloud /L CO is often denoted as the CO-to-H 2 conversion factor X CO and X CO ∝ M cloud /L CO . 4 The following estimation supports the hypothesis that the dense clumps are indeed protected even in low metallicity, and thus, can radiate CO(3-2) emission as much as their counterparts in high metallicity.
The clouds of the same mass M cloud in high and low metallicity have different amounts and distributions of CO molecules (Figure 8a,b). The ones in low metallicity show a smaller L CO , and hence larger X CO . This difference is detected using the virial mass measurements of molecular clouds with CO(1-0): X CO(1−0) = 3 × 10 20 in the Milky Way (∼ 1Z ; Scoville et al. 2016) and 7×10 20 in the LMC (∼ 0.4Z ; Fukui et al. 2008). Therefore, for the same mass, clouds in the LMC emit only about 40% of L CO(1−0) compared to their Galactic counterparts of the same M cloud . In other words, only about 40% of molecular gas, presumably at cloud centers, would have CO molecules at the LMC metallicity.
The 40% of molecular gas that has CO likely hosts dense clumps at the cloud centers (Figure 8b). The CO(3-2) emission requires higher density and/or kinetic temperature for excitation than CO(1-0) and CO(2-1),  Note that the beam-deconvolved radius R is not listed, since it is only an upper limit (R < 0.48 ) for all clouds except Cloud 3 (R = 1.05 ) whose measurement suffers from the spurious extension (see Section 4.2). and preferentially traces the dense clumps rather than the bulk gas (e.g., Komugi et al. 2007;Wilson et al. 2009). Wilson et al. (2009) compared the CO(3-2) luminosity L CO(3−2) against the bulk mass M cloud in nearby galaxies at solar metallicity (∼ 1Z ). Their result suggests that the detected L CO(3−2) are only 34% of those expected from L CO(1−0) on an assumption of thermalized gas (i.e., the CO 3-2/1-0 ratio in the brightness temperature is about unity). Thus, only about 34% of molecular gas in clouds radiates CO(3-2) emission at the high metallicity (more accurately, the molecular gas that emits about 34% of L CO(1−0) radiates CO(3-2) emission). Therefore, at LMC metallicity, we expect the dense clumps emitting at CO(3-2) (34% of the molecu-lar gas) are located in the same region as the gas that is protected from photo-dissociation and has CO molecules (40% of the molecular gas), as in Figure 8b. Leroy et al. (2022) recently found a CO 3-2/1-0 ratio of 0.31 (0.20-0.42 at 16th to 84th percentiles) as galaxyintegrated mean values among 20 massive solar metallicity galaxies. This is consistent with our adopted ratio of 0.34 by Wilson et al. (2009). Our discussions in this paper are for typical, Orion-class SF regions (Figure 3.1), but the ratio may increase mildly in more intense SF regions (Onodera et al. 2012;Miura et al. 2012;Vlahakis et al. 2013;Miura et al. 2014;Morokuma-Matsui & Muraoka 2017). spectra of the clouds produced with the 1 km s −1 cube. The spatial pixels and volume that each cloud occupies in the (RA, DEC, Vel) space are determined with the 2.54 km s −1 cube. The gray spectra show the sum over all the spatial pixels including noises in each 1 km s −1 channel. The red spectra include only the pixels with > 3σ emission in the cloud volume. Figure 6. Cloud radius-velocity dispersion plot. M83 XUV clouds are blue (and the blue arrows are the upper limits in radius, while they are resolved in velocity). Reference points in the background are for the MW's inner disk (Solomon et al. 1987), outer disk (r gal ∼12-14 kpc; Sun et al. 2017), and far outer disk (∼14-22 kpc; Sun et al. 2015). Large and Small Magellanic Clouds (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).

Masses of Molecular Clouds and Clumps
Even when the cloud outer layers become CO-dark in low metallicity, there are still abundant H 2 molecules there, and M cloud does not depend on the metallicity. In our hypothesis, the dense clumps are protected from photo-dissociation and would radiate a similar amount of L CO(3−2) in high and low metallicities. Hence, the calibration of M cloud /L CO(3−2) in high metallicity should be applicable in low metallicity, and the conversion equation from L CO(3−2) to M cloud derived in high metallicity can be used here. We adopt the calibration by Wilson et al. (2009), where the CO(3-2) flux S CO(3−2) dv is related to the luminosity L CO(3−2) = 4πD 2 S CO(3−2) dv. R 3−2/1−0 is the CO 3-2/1-0 brightness temperature ratio averaged over galaxies with high metallicity. The X CO(1−0) term is from the calibration made in the high metallicity, and hence the X CO(1−0) in high metallicity should be used here even when it is applied in low metallicity. M cloud from CO(3-2) is listed in Table 3 and ranges from 8.2 × 10 2 M to 2.3 × 10 4 M . The derived masses are much lower than the previous detections in XUV disks (Braine & Herpin 2004;Braine et al. 2007Braine et al. , 2010Braine et al. , 2012Dessauges-Zavadsky et al. 2014). From their high abundance in our region, the low-mass molecular clouds are likely a more common cloud population in XUV disks than previously detected. The most massive cloud (Cloud 4) has a mass similar to that of Orion A. Its SF activity measured in L Hα is also similar to that of the Orion Nebula in Orion A (Section 3.2). Note again that among massive star-forming regions in the MW, Orion A is one of the least impressive. All the other clouds in our study have HII regions in their proximity, but are less massive and less active in SF. The least massive cloud of ∼ 800M , only 4% of Orion A's mass, is unlikely to form OB stars. The SF in the XUV disk appears to be limited by the cloud masses.
The CO(3-2) emission is mainly from dense clumps in the clouds since it requires high density and temperature for excitation. We attempt a crude estimation of the clump masses M clump under a set of assumptions: (1) the clumps do not suffer from photo-dissociation and have CO molecules, (2) X CO can be applied to a part (clump) of a cloud, and (3) the gas is thermalized (the CO 3-2/1-0 brightness temperature ratio within the clump, R 3−2/1−0 , is unity). In this case, we should adopt X CO = 3 × 10 20 from the Galactic measurement even for the low metallicity, because no part of the clump is photo-dissociated. It gives The 1σ sensitivity in the 2.54 km s −1 cube of CO(3-2) is 88M . M clump ranges from M clump = 2.7 × 10 2 M to 7.5 × 10 3 M (Table 3). Given the assumptions, we expect the total uncertainties to be not smaller than a factor of 2-3. Again, the clump mass in Cloud 4 is about the same as that of the OMC-1 clump in the Orion A cloud.  Table 2. Blue or yellow circles are molecular clouds and are the same as the ones in Figure 2d with the labels of '1', '2', '3', ...  Figure 8. Schematic illustrations of the physical and chemical structures of a molecular cloud (a) in high-metallicity environment, and (b) in low-metallicity environment. The physical structures (i.e., mass distribution and SF activity) are common between (a) and (b), but their chemical structures (i.e., CO abundance) are different. The ambient radiation penetrates into the cloud, dissociates CO molecules in the cloud outer layers, and makes them CO-deficient and CO-dark, while more abundant H2 molecules can self-shield themselves and remain (gray). These CO-dark layers thicken at low metallicity, since a lower dust abundance results in less efficient extinction and allows the radiation to penetrate deeper into the cloud (e.g., Maloney & Black 1988;van Dishoeck & Black 1988). Still, the CO molecules can remain near the cloud center (orange). We hypothesize that the dense clumps (brown), often associated with SF (yellow), are located at the heart of the cloud and likely remain intact with abundant CO molecules. The CO(3-2) emission originates from the dense clump region.
Note that the two assumptions we adopted above -CO being confined only in the clump and being thermalized -minimize the predicted CO(2-1) flux for the detected CO(3-2) flux. If the CO(2-1) emission is also from a cloud envelope, it should add to the CO(2-1) flux. If the gas is not thermalized, the CO 3-2/2-1 flux ratio would be lower, and the predicted CO(2-1) flux would be higher. Either way, the CO(2-1) emission from Cloud 4 should have been detected. In other words, it is difficult to deviate significantly from the conditions we assumed.
The CO(2-1) non-detection also excludes the possibility of a CO-bright envelope as large as Orion A's (∼ 20 pc). The CO(2-1) data can resolve such an envelope at the beam size of ∼ 17 × 12 pc 2 at a high sensitivity of 58 mK (1σ) ( Table 1). The non-detection suggests that its average surface brightness is < 0.17 K (3σ). Hence, the envelope is CO(2-1)-dark.
An independent constraint on M cloud may be derived from the dust continuum emission. Using the calibration by Scoville et al. (2016) in a range of environments and adopting the power-law dependence of the dust emissivity on frequency ν as κ ν ∝ ν β , we obtain where and ν 850µm ≈ 353 GHz. The X CO term is due to the calibration in high metallicity. We adopt a dust emissivity index of β = 1.8 (Planck Collaboration et al. 2011) and dust temperature of T d = 25 K (Scoville et al. 2016). The detection limits from the 225 and 349 GHz dust continuum measurements correspond to 1.5 × 10 6 M and 4.5 × 10 4 M (3σ), respectively. The M cloud from CO(3-2) is consistent with these limits (Table 3).

Star Formation On 1-kpc Scale
On a 1 kpc-scale average, the SF efficiency in total gas content (HI+H 2 ) is much lower in XUV disks than in the optical disks of nearby galaxies. The average gas consumption timescale by SF, i.e., the inverse of SF efficiency, is about 2 Gyr in the optical disks (Bigiel et al. 2008;Leroy et al. 2008; see also Kennicutt 1998;Kennicutt & Evans 2012), while it is as long as 100 Gyr over the XUV disk of M83 (Bigiel et al. 2010). As another reference, the gas consumption timescale in the Orion A cloud is about 520 Myr (see Section 3.2 including the caveat on this estimation).
The region of this study is 1 kpc 2 (more accurately, 0.99 kpc 2 by taking into account the round corners of the observed region). The total HI mass in the region is M HI ∼ 7.4 × 10 6 M from the VLA+GBT data. The total H 2 mass is derived by summing up all the clouds (Table 3) and is M H2 ∼ 8.4 × 10 4 M , which includes the CO-dark H 2 gas. Only ∼ 1% of the gas is in the molecular phase. The stellar initial mass function (IMF) may not be fully populated in this low SF region (Koda et al. 2012a). Nevertheless, if we simplistically adopt the standard calibration that assumes a fully populated IMF (Kennicutt & Evans 2012), L Hα ∼ 6.8 × 10 37 erg s −1 (Section 3.1) translates to a SFR of ∼ 3.7 × 10 −4 M yr −1 .
The atomic and molecular gas in total, including Helium (applying a factor of 1.36 to the mass), would be consumed in about 27 Gyr at this SFR. This is a few times shorter than the average across M83's XUV disk, but is an order of magnitude longer than the average of 2 Gyr in the optical disks and is longer than the Hubble time. On the other hand, the consumption timescale of the molecular gas alone would be as short as 310 Myr, which is similar to that in Orion A (Section 3.2).
Therefore, the molecular clouds in the XUV disk on average appear to be forming stars as much as a typical cloud in the MW, but the XUV disk as a whole is not producing stars as much. This supports the hypothesis that the clouds in the XUV disk and Galactic disk share a similar mass structure, and hence, SF activity. The key to the low SF activity in XUV disks should be in the conversion efficiency from the atomic gas to molecular clouds. This is a possibility suggested by Bigiel et al. (2010), and the resolved analysis of the clouds in this study supports it.

Comparisons to Previous Studies
CO emission was previously detected in the outskirts of four galaxies with XUV disks. Three of them, NGC 4414, NGC 6946, and M63, were observed with a singledish telescope alone at a relatively low sensitivity and resolution ( 10 5−6 M and 300-500 pc, Braine & Herpin 2004;Braine et al. 2007;Dessauges-Zavadsky et al. 2014). If our 1 kpc 2 region in M83 were to be observed at these resolutions, multiple clouds would be in a single telescope beam (Figure 2). A detection in M33 (Braine et al. 2010) was resolved into a single cloud of ∼ 4×10 4 M with an interferometer (Braine et al. 2012). Our study, together with Braine et al. (2012), demonstrates that high resolution and sensitivity are necessary to identify individual clouds and to investigate the environment of SF in the XUV disks.
There are large-scale CO(1-0) surveys in the outer disk of the MW (e.g., Wouterloot & Brand 1989;Digel et al. 1994;Brand & Wouterloot 1994;Snell et al. 2002;Yang et al. 2002;Brunt et al. 2003;Nakagawa et al. 2005;Sun et al. 2015). The optical edge of the MW disk is around R 25 ∼ 13.4 kpc (Goodwin et al. 1997(Goodwin et al. , 1998 with considerable uncertainty due to its edge-on projection. The galactocentric radius of 1.24R 25 , corresponding to the 1 kpc 2 region in M83, is r gal ∼ 16.7 kpc for the MW. The largest, most recent CO survey by Sun et al. (2015) detected clouds in the far outer MW disk (r gal 14 kpc), which include 51 clouds of 10 2−4 M in the extreme outer disk (r gal 16 up to 21 kpc) with 11 of them in the most massive end of 10 4 M . They are spread over a large area of ∼ 80 kpc 2 in the Galactic longitude of l ∼100-150 • . Therefore, the cloud density is low. The majority are associated with a possible extension of the Scutum-Centaurus spiral arm traced in HI (Sun et al. 2015). SF exists even in the extreme outer disk, but they appear to be mostly low-mass SF (later than B0type) according to an analysis of the Wide-field Infrared Survey Explorer (WISE) data (Izumi et al. 2017, see their Figure 15), while much obscured higher-mass SF has yet to be found.
Our 1 kpc 2 region around an HI peak in the XUV disk of M83 has 23 clouds of 10 3−4 M , including three 10 4 M clouds. Compared to the extreme outer MW disk, the cloud density is high. Obviously, a much larger sample is required to make a general conclusion, but the new detection appears to indicate that molecular clouds, especially their abundance, are one of the key factors for the XUV activity. Figure 6 shows clouds in the inner, outer, and far outer MW disk traced in CO(1-0) (Solomon et al. 1987;Sun et al. 2015Sun et al. , 2017. The XUV clouds in M83, traced in CO(3-2), generally show larger velocity dispersions than the MW clouds. This can be explained if the CO(3-2) emission primarily traces the dense clumps embedded in the clouds (Section 4.3).

Clouds and Star Formation in XUV Disks
The ALMA CO(3-2) observations constrain the properties of 23 molecular clouds in the XUV disk, having relatively small bulk gas masses of M cloud = 8.2 × 10 2 to 2.3 × 10 4 M , and dense clump masses of M clump = 2.7 × 10 2 to 7.5 × 10 3 M . The ubiquity suggests that these relatively low-mass clouds are a common population of molecular clouds and are the main drivers of the SF activity in the XUV disk.
Numerous previous CO(1-0) and CO(2-1) searches in XUV disks resulted in non detection (Watson & Koda 2017). A small number of previous detections could not isolate individual clouds due to a combination of low sensitivity and resolution ( 10 5−6 M and 300-500 pc), and one successful case, in the very close galaxy M33, identified only one single cloud of 4.3 × 10 4 M (Braine et al. 2012). The common cloud population turns out to be less massive than the previous detections. This study is the first demonstration of the utility of CO(3-2) to detect molecular clouds.
A comparison to Galactic clouds provides insights into the clouds and SF in the XUV disk. For example, a cloud like Orion A -having M cloud = 4.0 × 10 4 M and M clump = 7.4 × 10 3 M -can be one of the most massive clouds in the observed region. The SF activity in the region, traced by Hα, is also similar to that of the Orion Nebula in Orion A. Orion A is one of the least impressive massive star-forming regions in the MW, but would be the most impressive in the XUV disk. Clouds smaller than Orion A are also found in the XUV disk.
It is likely, and we hypothesize, that the internal mass distributions of molecular clouds are similar between the MW disk and XUV disks. The Galactic clouds, including Orion A, typically have star-forming dense clumps surrounded by thick layers of bulk molecular gas. In the case of low metallicity of XUV disks, the outer layers should have CO molecules photo-dissociated and become CO-deficient and CO-dark (Maloney & Black 1988;van Dishoeck & Black 1988;Wolfire et al. 2010). We further hypothesize that the dense clumps at the hearts of the clouds remain intact and maintain a high CO abundance even in the low metallicity (Figure 8). This structure consistently explains the detection of the CO(3-2) emission mainly from the dense clumps, and previous non-detection of the CO(2-1) emission. The common clump-envelope structure between the Galactic and XUV clouds is expected if they are governed primarily by their internal physics (e.g., by self-gravity) independent of their environment. This hypothesis can be tested with future observations, e.g., higher-resolution CO(3-2) observations for resolving the compactness of the clumps, deeper CO(2-1) or CO(1-0) observations for confirming the thermalization, and deeper dust continuum observations for independent mass estimation.
The common mass structure of clouds suggests that average clouds of the same mass should form stars at a similar rate independent of their environment. In other words, the SF efficiency on cloud scales, on average, would be universal. This is supported by the similar gas consumption timescales between the XUV clouds and Orion A.
This study finds only relatively low-mass molecular clouds in the XUV disk of M83. The absence of very bright HII regions in XUV disks may be due to the lack of more massive clouds. Koda et al. (2012a) took the standard stellar initial mass function (Salpeter 1955) as the probability distribution function of stellar mass, and calculated that a star cluster should have a total stellar mass of at least 10 3 M to populate one 40M star (O6 or O7-type; Sternberg et al. 2003). The most massive cloud in our sample has M cloud = 2.3×10 4 M , requiring a cloud-to-star mass conversion efficiency of ∼ 4% to produce a 10 3 M cluster. This is about the efficiency in Galactic molecular clouds (Enoch et al. 2007). The other clouds have even smaller M cloud and are unlikely to form very massive stars and prominent HII regions.
The low SF efficiency on a 1 kpc scale has to originate on scales larger than the clouds; for example, in the conversion efficiency of the extended HI gas into molecular clouds (Bigiel et al. 2010). The mechanism for atomic to molecular gas conversion should also explain the cloud mass function and the lack of massive molecular clouds. Further studies are required to elucidate these aspects, but some clue may already be in the distribution of the CO(3-2) clouds/clumps. The long, chain-like distributions over 600 pc, and potentially the velocity gradient across our field, indicate that their formation is triggered on large scales, such as by large-scale galactic flows.
The hypothesis of the common cloud mass structure (i.e., the clump-envelope structure) in conjunction with the selective photo-dissociation explains the clouds and SF activity over the optical and XUV disks. This simplistic view could be modified in future studies. For example, this hypothesis exhibits merely an average picture over cloud populations and evolutionary sequences (e.g., see Vázquez-Semadeni et al. 2018). In other extreme environments, such as in the dense Galactic center region, the cloud properties show environmental influences (e.g., Oka et al. 2001) although such regions are relatively small in volume and mass.

Application to High-z and Other Environments
The common cloud structure ( Figure 8) suggests a constant M cloud /L CO(3−2) between high-and low-metallicity environments, even when CO molecules in cloud outer layers are photo-dissociated. This has an important implication for CO observations of distant galaxies. The red-shifted high-J transitions of CO (higher than J=1 or 2) are often used to trace their total gas mass (Carilli & Walter 2013;Combes 2018). Those high-J CO emissions likely capture only starforming dense clumps, but not the bulk cloud gas. In fact, the CO(1-0) and CO(2-1) luminosities are often beyond those expected from higher-J transition analysis (Ivison et al. 2011;Carilli et al. 2011;Riechers et al. 2011;Bothwell et al. 2013;Daddi et al. 2015) Even then, the high-J CO emissions can trace the total gas mass if the self-gravitating clouds have the same internal structure between the local and distant universe, and if the average ratio between the dense clumps and their parent clouds is about constant independent of environment and redshift as expected for self-gravitating clouds.

SUMMARY
The new ALMA CO(3-2) observations detected a large number of molecular clouds across a 1 kpc 2 area of the XUV disk of M83. Their high abundance suggests that the 23 clouds represent a common population of molecular clouds and are likely the main drivers of the SF activity in XUV disks.
We hypothesized that these clouds share, on average, the same common structure (mass distribution) as Galactic clouds, such as Orion A; having starforming dense clumps embedded in thick layers of bulk molecular gas. In the low-metallicity regime of XUV disks, the outer layers of the clouds are likely COdeficient and CO-dark, while still abundant in H 2 , due to selective photo-dissociation. This hypothesis consistently explains the CO(3-2) detection and CO(2-1) nondetection. The CO(1-0) and CO(2-1) emission should be reduced significantly due to the lack of CO molecules in the photo-dissociated outer layers where a large fraction of the CO(1-0) and CO(2-1) emission normally originate. The dense clumps have sufficiently high density and temperature for CO(3-2) excitation, reside at the hearts of the clouds, and thus, are protected from the dissociation. The common cloud mass structure independent of environment is conceivable if it is governed by their internal physics, e.g., by self-gravity, while their chemical structure could still depend on the environment (e.g., metallicity and radiation field).
The common cloud structure constrains the geometry and physical conditions within molecular clouds, which permits us to estimate their masses. The total masses of the clouds, including the CO-dark H 2 component, ranges from M cloud = 8.2 × 10 2 to 2.3 × 10 4 M . Their star-forming dense clumps have masses of M clump = 2.7 × 10 2 to 7.5 × 10 3 M . The most massive clouds and their dense clumps are similar to Orion A and its dense clump OMC-1. The dense clumps should also radiate CO(2-1) emission, but the flux is lower than the detection limit of the previous CO(2-1) study. This study demonstrates, for the first time, the utility of CO(3-2) emission to find molecular clouds.
The SF activity and efficiency are also similar between the clouds in the XUV disk and Orion A. The XUV disk as a whole shows a much lower SF activity, which is not due to peculiar cloud properties, but due to the lack of massive clouds. The abundance and mass function of molecular clouds should be regulated by the (unknown) trigger of cloud formation from the abundant atomic gas. Further studies are necessary to elucidate this mechanism, but a clue is already found in the distribution of the clouds. The long chain-like distribution ( 600 pc in length) indicates that the trigger is on large scales.
The common mass structure and SF efficiency among molecular clouds should justify the use of high-J CO transitions as tracers of bulk molecular gas mass in distant galaxies, assuming that the bulk molecular gas mass is proportional to the dense clump mass traced by the high-J transitions.
We thank the referees for useful comments, in particular, the first referee for suggesting the discussion in Section 5.5 and the second referee for the advice on the title. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2013.1.00861.S and 2017.1.00065.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. The Green Bank Observatory and the National Radio Astronomy Observatory are facilities of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.. This research is based in part on data collected at the Subaru Telescope, which is operated by the National Astronomical Observatory of Japan (NAOJ). We are honored and grateful for the opportunity of observing the Universe from Maunakea, which has the cultural, historical, and natural significance in Hawaii. Data analysis is in part carried out on the Multi-wavelength Data Analysis System operated by the Astronomy Data Center (ADC) in NAOJ. JK acknowledges support from NSF through grants AST-1812847 and AST-2006600 and from NASA grant NNX14AF74G. MR wishes to acknowledge support from ANID(CHILE) through FONDECYT grant No1190684.