A sensitive, high-resolution, wide-field IRAM NOEMA CO(1-0) survey of the very nearby spiral galaxy IC 342

We present a new wide-field 10.75 x 10.75 arcmin^2 (~11x11 kpc^2), high-resolution (theta = 3.6"~ 60 pc) NOEMA CO(1-0) survey of the very nearby (d=3.45 Mpc) spiral galaxy IC 342. The survey spans out to about 1.5 effective radii and covers most of the region where molecular gas dominates the cold interstellar medium. We resolved the CO emission into>600 individual giant molecular clouds and associations. We assessed their properties and found that overall the clouds show approximate virial balance, with typical virial parameters of alpha_vir=1-2. The typical surface density and line width of molecular gas increase from the inter-arm region to the arm and bar region, and they reach their highest values in the inner kiloparsec of the galaxy (median Sigma_mol~80, 140, 160, and 1100 M_sun/pc^2, sigma_CO~6.6, 7.6, 9.7, and 18.4 km/s for inter-arm, arm, bar, and center clouds, respectively). Clouds in the central part of the galaxy show an enhanced line width relative to their surface densities and evidence of additional sources of dynamical broadening. All of these results agree well with studies of clouds in more distant galaxies at a similar physical resolution. Leveraging our measurements to estimate the density and gravitational free-fall time at 90 pc resolution, averaged on 1.5 kpc hexagonal apertures, we estimate a typical star formation efficiency per free-fall time of 0.45% with a 16-84% variation of 0.33-0.71% among such 1.5 kpc regions. We speculate that bar-driven gas inflow could explain the large gas concentration in the central kiloparsec and the buildup of the massive nuclear star cluster. This wide-area CO map of the closest face-on massive spiral galaxy demonstrates the current mapping power of NOEMA and has many potential applications. The data and products are publicly available.


Introduction
Giant molecular clouds (GMCs) represent the immediate sites of star formation and set the boundary conditions for star formation (e.g., Kennicutt & Evans 2012).Far from universal, the properties of GMCs strongly depend on the local galactic environment (e.g., Gratier et al. 2010;Hughes et al. 2013;Colombo et al. 2014;Rosolowsky et al. 2021).In theory, the properties of molecular clouds play a key role in setting the rate at which molecular gas can form stars (e.g., Padoan et al. 2012;Federrath & Klessen 2013;Krumholz et al. 2019).These theoretical predictions have spurred a number of recent observational attempts to quantitatively link the properties of individual molecular clouds to the rate of star formation in galaxies or parts of galaxies (e.g., Barnes et al. 2017;Leroy et al. 2017;Ochsendorf et al. 2017;Hirota et al. 2018;Kreckel et al. 2018;Utomo et al. 2018;Schruba et al. 2019;Pessa et al. 2021).
To make such measurements, one must survey the properties of molecular gas with high resolution and high completeness across large parts of galaxies.Most often, this is done by observing low-J CO line emission, the standard tracer of molecular gas in galaxies (e.g., Bolatto et al. 2013).Though CO is among the brightest millimeter-wave lines, CO surveys still represent a major technical challenge, requiring a focused effort using the best millimeter and submillimeter telescopes in the world.
At a distance of only d = 3.45 Mpc (Anand et al. 2021), IC 342 represents one of the closest, massive, star-forming spiral galaxies beyond the Local Group (see Table 1).Among massive, vigorously star-forming disk galaxies, its proximity is rivalled only by Andromeda, NGC 253, and NGC 4945 and all of these systems are far more inclined than IC 342.This proximity and face-on orientation (i ∼ 30 • ) translate into a sharper physical resolution and less confusion due to galaxy projection across the entire electromagnetic spectrum.Even among nearby northern galaxies, IC 342 stands out as it is a prime example of a bulgeless spiral galaxy hosting a massive, compact nuclear star cluster (Carson et al. 2015;Neumayer et al. 2020) which is among the brightest ones known in nearby late-type galaxies (Georgiev & Böker 2014).Indeed, were it not for its position behind the Galactic plane (b = 10.6 • ), IC 342 would certainly be among the best-known galaxies in the sky.Even given this situation, the gas-rich nuclear region of IC 342 has been a key target of millimeter-wave studies for decades (e.g., Lo et al. 1984;Eckart et al. 1990;Downes et al. 1992;Sakamoto et al. 1999;Meier & Turner 2005;Pan et al. 2014).This rich region hosts one of the nearest nuclear star clusters that appears to heavily impact its surrounding gas reservoir (Schinnerer et al. 2003) and the nuclear region of IC 342 is one of the closest starbursts (we estimated a star formation rate inside the central 1.5 kpc equal to ∼0.2 M ⊙ yr −1 ; see section 2.2).Thus, there are some similarities with the less massive, very nearby galaxies M33 and NGC 5253 Fig. 1.Molecular, atomic, and ionized gas in IC 342.This is a false-color composite of IC 342 showing molecular gas traced by our new NOEMA CO (1−0) survey (green), along with stellar light from DSS2 (white), atomic gas traced by the 21 cm line using the VLA (blue), and Hα emission from ionized gas (red).The dashed rectangle delimits the field of view observed with NOEMA and solid ellipses show galactocentric radii at 5, 10, 15, and 20 kpc.
regarding the presence of nuclear stellar clusters and nuclear starbursts.In any case, the proximity, orientation, and active star formation that allowed these earlier studies make IC 342 an ideal target for wider area CO surveys using more modern facilities.
Here, we report on a new NOrthern Extended Millimeter Array (NOEMA) wide-area survey of CO emission of IC 342.We used NOEMA, supported by the IRAM 30-m telescope, to observe CO (1−0) emission from ∼1000 individual pointings across the inner region of IC 342.In total, we map 10.75 × 10.75 arcmin 2 , that is, about 11 × 11 kpc 2 , at an angular resolution of ∼3.6 ′′ (or ∼60 pc).Here we describe this new survey, present a catalog of 637 molecular clouds from the galaxy, and compare their properties across different environments.We also assess how well environmental variations in the mean surface density of molecular clouds predict apparent variations in the molecular gas depletion time, testing some of the theoretical predictions described above.
IRAM NOEMA is the most powerful millimeter-wave interferometer in the northern hemisphere and one of the few telescopes in the world with the ability to conduct such CO surveys.With twelve 15-m dishes in the French Alps, recently upgraded receivers, and efficient observing, NOEMA can rapidly survey molecular line emission over hundreds or thousands of individual pointings.While the capabilities of ALMA are widely appreciated, NOEMA can survey CO emission at high completeness and high resolution from crucial northern targets such as the prototypical starburst M82, the iconic spirals M51 and M101, the Local Group galaxies M31 and M33, and the "fire- The images show spiral arms, abundant inter-arm emission, and a velocity field that mostly reflects a regularly rotating gas disk.The bottom panels show the molecular gas traced by CO in the context of atomic gas and stars.In these panels, the CO and H I data are shown at a common ∼350 pc resolution.The bottom left panel shows a map of the total neutral gas surface density adding molecular gas surface density to atomic gas surface density from Chiang et al. (2021).The white ellipses show 1, 2, 3 × R e and the gray dashed ellipse indicates 1 R 25 .The lower right panel shows the azimuthally averaged mass surface density profiles for atomic gas, molecular gas, and stellar mass estimated from the near-infrared.The stellar mass distribution using HST imaging is consistent with a nuclear star cluster and an exponential disk (Carson et al. 2015), while the presented profile suffers from the resolution of the WISE data.Our NOEMA survey covers the inner molecule-dominated region where the H I emission is depressed, including the CO-bright center.
works" galaxy NGC 6946 (all northern with d < 10 Mpc).The Plateau de Bure Interferometer, NOEMA's predecessor, demonstrated these capabilities with the PAWS survey of M51 (Schinnerer et al. 2013;Pety et al. 2013) and a 2-3 mm molecular line survey of NGC 6946 (Eibensteiner et al. 2022).Krieger et al. (2021) recently used NOEMA to produce a wide-field, high resolution map of CO emission in M82, and Bešlić et al. (2021) mapped with NOEMA the J = 1−0 lines of two CO isotopologues ( 13 CO and C 18 O) and of high critical density molecular Notes.Morphological Hubble type, distance, and kinematic orientation parameters of IC 342 according to the reference provided, the NASA/IPAC Extragalactic Database (NED) a , or LEDA b .R 25 is the 25 mag/arcsec 2 B-band isophotal radius, while the effective radius contains half of the light.E(B− V) represents optical extinction, 12 + log(O/H) R 0 is the gas-phase metallicity, and ∆log(O/H)/R gal is the metallicity radial gradient (where R gal is galactocentric radius).M star , M atom , M mol , and M dust represent the total stellar, atomic, molecular, and dust mass in this galaxy.SFR is the total star formation rate.All masses were scaled to our adopted distance. (a) Oxygen abundance gradient constrained by only five HII regions (McCall et al. 1985). (b) We estimated the molecular gas mass using our short-spacing-corrected NOEMA data (total CO luminosity multiplied by the Galactic α CO ).Since the NOEMA field of view is limited, we extrapolated to the whole galaxy according to the proportion of the total WISE 22 µm luminosity that falls within the NOEMA field of view; as WISE 22 µm is a SFR tracer, it should roughly scale according to the molecular mass at each radius.
a https://ned.ipac.caltech.edub https://leda.univ-lyon1.frtracers (HCN, HNC, and HCO + ) towards the nearby strongly barred galaxy NGC 3627.The data processing (calibration, merging, imaging, and deconvolution) followed the standard procedures of both observatories.Details are given in Appendix A. Briefly, data were calibrated using GILDAS/CLIC.We used reference quasars as phase and amplitude calibrators, and bright Galactic sources as flux calibrators.Around 20% of the data were flagged due to unstable weather.Imaging was carried out using GILDAS/MAPPING.First, we constructed pseudo-visibilities out of the IRAM 30m data and they were merged with the NOEMA data in the uv plane using the task uvshort.The imaging process involved cleaning with the Högbom algorithm using a broad mask based on a smoothed and clipped version of the IRAM 30m cube.The resulting data cube has an rms of 0.11 K over 5 km s −1 channels.We constructed data cubes with 5, 10, and 20 km s −1 spectral resolution.These three cubes were deconvolved separately.Because the deconvolution becomes more accurate when the signal-to-noise ratio of the data increases, we were able to deconvolve more flux in our 20 km s −1 resolution cube.The three cubes share the same angular resolution of 4 × 3.25 arcsec 2 with a position angle of 91 • .Table A.1 lists the salient features of the observations and reports the properties of our 5 km s −1 resolution data cube.

Data
We compared the flux between the NRO 45m data at 20 ′′resolution published in the Nobeyama CO Atlas of Nearby Spiral Galaxies (Kuno et al. 2007) with the IRAM 30m data at 22.5 ′′ resolution described here.We used the same channel spacing of 5 km s −1 .We computed the flux in a field of view of 600 ′′ × 600 ′′ centered on the galaxy center and within the velocity range [−80, +140] km s −1 .For a distance of 3.45 Mpc, we found a total luminosity of 7.04 and 7.00 K km s −1 pc 2 , for NRO and IRAM, respectively.The difference amounts to only 0.6%.We then compared the flux between the IRAM 30m data cube, the NOEMA+30m combined cube, and the NOEMA-only cube.We here used the same channel spacing of 5 km s −1 and the same velocity range as above, but we restricted the field of view to 500 ′′ × 500 ′′ to avoid edge effects.We first regrid the 30m data onto the NOEMA spatial grid, and we smoothed the NOEMA+30m and NOEMA-only cubes to the spatial resolution of the 30m data, that is, 22.5 ′′ .We found a luminosity of 6.20, 6.14, and 0.95 K km s −1 pc 2 for the IRAM 30m, NOEMA+30m, and NOEMA-only cubes.The flux scales for the combined and single-dish data cubes thus agree within 1%.Only 15% of the total flux is recovered with the interferometer.
From the deconvolved data cubes, we produced a set of derived products including cubes with round beams at a series of fixed angular (5, 15, 21 arcsec) and physical (90, 150, 500, 1000, 1500 pc) resolutions, noise estimates, masks, and accompanying moment maps and uncertainties.The moment maps were produced using a dilated mask approach that follows the strategy presented in Leroy et al. (2021b) which was applied to the products delivered for PHANGS-ALMA (Leroy et al. 2021a).Line widths were measured using the effective width approach (following Heyer et al. 2001), which is often more robust than second-order moment maps: σ = I CO /( √ 2πT peak ).Details are provided in Appendix B. Data cubes and associated data products are publicly available 1 .
The conversion from CO(1-0) integrated intensity to molecular gas surface density follows where α CO represents the conversion factor from CO(1-0) to H 2 surface density (including a factor 1.36 to account for helium and heavier elements).The factor cos i corrects the surface densities for inclination.Throughout this paper, we apply a fixed α CO = 4.35 M ⊙ pc −2 (K km s −1 ) −1 to IC 342 whenever converting from CO intensity to molecular gas mass or surface density.This value is appropriate for the Milky Way and typical of the disks of massive spiral galaxies (Bolatto et al. 2013).We do expect α CO to vary across IC 342, approaching a lower value in the center (e.g., Downes & Solomon 1998;Sandstrom et al. 2013;Chiang et al. 2021) and to anticorrelate with the known metallicity gradient of IC 342 (Pilyugin et al. 2014, and see Table 1).
Recent dust-based analysis points to an α CO value in the central ∼kpc which is ∼4 times below the Galactic value (I.-D.Chiang, priv. comm.;Chiang et al. 2023).According to Meier & Turner (2001) and Pan et al. (2014), in the central few hundred parsec of IC 342, α CO is ∼2-3 times lower than the Galactic value; Pan et al. (2014) suggest that the disk of IC 342 reaches α CO values two times above the Galactic α CO (albeit with significant fluctuations).This means that we can expect up to a factor of ∼5 bulk difference in α CO between disk and center.Based on this, we note the cases below where deviation from the adopted constant conversion factor likely influences our results.

Supporting data and measurements
We used VLA H I 21 cm line observations from Chiang et al. (2021) to trace neutral atomic gas.Specifically, we applied the following equation to transform H I observations to atomic gas surface density, which assumes optically thin 21 cm emission: where Σ atom includes the 1.36 factor to account for helium and heavier elements.The cos i term again corrects for galaxy inclination.
We used a WISE 3.4 µm image to trace stellar mass and a WISE 22 µm image to estimate the local recent star formation rate (SFR) surface density (Wright et al. 2010).The WISE images were processed into an estimate of stellar mass following the procedures in Leroy et al. (2019Leroy et al. ( , 2021c)).The WISE 3.4 µm image has significant contamination by foreground stars because IC 342 lies behind the Galactic plane.But since the only application of these data in this paper is large regional averages and radial profiles, we do not expect this to represent a large concern.The 22 µm image is not heavily affected by the Galactic foreground, and we proceed with the WISE 22 µm alone as our SFR tracer.We used the following prescription to convert the WISE 22 µm map to SFR: This prescription follows Leroy et al. (2021a) and assumes a Chabrier initial mass function (Chabrier 2003), which agrees with extinction-corrected Hα measurements of SFR from the PHANGS-MUSE survey (Belfiore et al. 2022) within ∼20−30%.In this case, we also use a Mayall MOSAIC Hα image kindly provided by K. Herrmann (private communication) for a visual comparison.
We also made heavy use of an environmental "mask" based on the near-IR morphology of the galaxy following Querejeta et al. (2021).Appendix B.3 describes the construction of this mask in more detail and explains how we identify distinct "center," "bar," "arm," and "inter-arm" regions.
We aggregated information from the cloud catalogs, SFR estimates, and molecular gas mass estimates into a set of regional property estimates.This aggregation follows the methods described in Sun et al. (2022) and resembles those deployed in Leroy et al. (2017) and Utomo et al. (2018).Briefly, we calculated the mass-weighted average of cloud properties as well as the aperture average of SFR and molecular gas surface densities in a series of 1.5 kpc hexagonal apertures, which together tile the entire footprint of the NOEMA+IRAM observations.Details about this data aggregation scheme can be found in Sun et al. (2022).We used these derived products to compare populationaveraged estimates of the cloud gravitational free-fall time to the regional molecular gas depletion time.These measurements can be found in Table A.3.

Results
Figure 1 shows the molecular gas distribution as mapped by NOEMA+IRAM 30-m, together with the distributions of Hα line emission tracing gas ionized by young, massive stars (K.Herrmann private communication), and neutral atomic gas traced by the 21 cm line (Chiang et al. 2021).Together, these three tracers reveal an intricate network of multiple spiral arms that emanate from the bright center.The inner arms visible in CO emission and dominated by molecular gas appear prominent inside a galactocentric radius of 5 kpc.These molecular gas arms connect seamlessly to the outer arms traced by 21 cm emission and dominated by atomic gas, which can be traced all the way out to radii of ∼20 kpc.The ionized gas emission follows these spiral arms, showing a patchy distribution that reflects the presence of many individual H II regions and indicating that massive stars are actively forming along the gas arms.
Figure 2 shows our CO (1−0) survey in more detail.The top left panel illustrates the overall structure of CO emission in the galaxy for the observed field of view.Though the galaxy has been mapped before by lower-resolution, lower-sensitivity, or smaller field-of-view observations (including Crosthwaite et al. 2001;Helfer et al. 2003;Kuno et al. 2007;Hirota et al. 2011), the NOEMA map offers the sharpest, most complete view to date of the CO distribution in IC 342.The peak intensity map in Fig. 2 shows multiple spiral arms, fine inter-arm filaments, and a bright, elongated inner structure, which is reminiscent of gas lanes along a stellar bar and encompasses the nuclear starburst (e.g., Downes et al. 1992;Schinnerer et al. 2003;Meier & Turner 2005).The spiral arm east of the center has a significantly higher contrast similar to the (eastern) H I spiral arms further out in the disk (see Fig. 1).We speculate that this might be related to a mild interaction with a neighboring galaxy, which has also been evoked to explain the warping of the outer H I disk (e.g., Rots 1979;Crosthwaite et al. 2000;Stil et al. 2005).The CO velocity field shown in the top right panel follows a gentle S-shape which is consistent with streaming motions associated with a bar potential on top of a regularly rotating gas disk.
The bottom row of Fig. 2 places our new CO survey in context.We combine the distribution of the atomic gas from Chiang et al. (2021) with our new molecular gas map at a matched resolution of ∼350 pc to show the total neutral gas surface density distribution (left panel) and to construct radial profiles of atomic gas, molecular gas, and stellar surface density (right panel).The effect of changing angular resolution is particularly apparent for the central bright molecular gas distribution where the position angle changes from a small offset to the east at ∼4 ′′ resolution to a small western offset at 20 ′′ resolution.This change is caused by lower surface brightness material outside the central 60 ′′ which becomes apparent at the lower angular resolution.This might be a consequence of the bar in IC 342 (see Sect. 4).
The radial profiles show the basic distribution of the cold ISM in IC 342.Within 6 kpc, the molecular gas and stellar mass profiles have a nearly constant offset in logarithmic scale, implying a molecular gas to stellar fraction Σ mol /(Σ ⋆ + Σ mol ) ≈ 0.20−0.25 for a constant Galactic α CO ; if we assume a 2-3 times lower α CO in the center and twice higher α CO in the disk, as suggested by Pan et al. (2014), the molecular to stellar fraction would drop to ∼10% for the center and reach ∼40% throughout the disk.The molecular gas surface density exceeds the atomic gas surface density out to at least r gal = 6.7 kpc.Outside approximately this radius, the atomic gas dominates the cold ISM.As shown in Fig. 1 and 2, the structure of the molecular gas appears to connect smoothly with the atomic gas, consistent with a density-or pressure-dependent phase change.This largescale picture agrees well with previous low-resolution work by Crosthwaite et al. (2001) and Pan et al. (2014).However, in contrast to previous work, our map captures the sharp concentration of gas at the galaxy center (r gal < 1 kpc) and resolves the CO emission from the galaxy into discrete molecular clouds.We explore the implications of these higher resolution observations in the next sections.

Giant molecular clouds in IC 342
We used the python implementation of the CPROPS algorithm (Rosolowsky & Leroy 2006) to generate a catalog of molecular clouds in IC 342.We apply the algorithm to the 90 pc and 5 km s −1 resolution cube and use the same methods developed by Rosolowsky et al. (2021) for the PHANGS-ALMA survey; see Appendix B.2 for details.These methods yield radius (R), line width (σ), and mass (M) for each cloud.We relied on the extrapolated and deconvolved GMC properties that account for the spatial and spectral response of NOEMA.When computing cloud radii, we followed R = η √ σ maj,d σ min,d with η = 1.18 as in Rosolowsky et al. (2021).We adopted a fixed, Galactic α CO as stated above.
The CPROPS catalog contains 637 clouds with masses between 10 5 and 4 × 10 7 M ⊙ (out of which 602 have radius measurements).The catalog is publicly available2 and an excerpt showing the first ten rows can be seen in Table A.2.Given the θ = 90 pc resolution and this mass range, these objects likely span the range between giant molecular clouds and giant molecular associations.As we see below, they appear approximately gravitationally bound.Fig. 3 visualizes their sizes and locations in the galaxy, with the clouds colored according to their dynamical environment (see Appendix B.3).
The top right and bottom left panels of Fig. 3 illustrate the dynamical state of these molecular clouds using plots of σ 2 /R versus Σ mol ≡ M/(2πR 2 ) space.Following Heyer et al. (2009) and Field et al. (2011), variants of this plot have become a standard way to assess the properties of a molecular cloud population (e.g., see Leroy et al. 2015;Sun et al. 2018Sun et al. , 2020b;;Rosolowsky et al. 2021).In the log-log plot of σ as a function of Σ mol (bottom right panel), clouds with fixed virial parameter α vir = 2KE/UE, that is, a fixed ratio of kinetic energy (KE) to self-gravitating potential energy (UE), follow a straight line with slope of 1. Conversely, clouds with fixed turbulent pressure follow straight lines with slope of −1 (Sun et al. 2018(Sun et al. , 2020b)). Figure 3 shows the IC 342 clouds, separated by environment (as illustrated in the top left panel).As a reference we show a large set of GMC properties derived at identical spatial resolution and similar sensitivity from PHANGS-ALMA (Rosolowsky et al. 2021, A. Hughes et al., in preparation).We note that PHANGS-ALMA relies on the CO(2-1) transition, while the IC 342 observations presented here mapped CO(1-0); the ratio of both lines, R 21 , shows some mild variations within galaxies (e.g., den Brok et al. 2021;Leroy et al. 2021b).The upper right panel shows individual clouds, while the lower left panel shows running medians in σ 2 /R as a function of Σ mol .In the lower left panel, we plot the mass-weighted median and 16−84% range in surface density for each dynamical environment in IC 342.
In the bottom right panel of Fig. 3, we show a complementary view of the molecular ISM that captures similar information to the top right and bottom left panels.Here, we use a lineof-sight based approach following Sun et al. (2018Sun et al. ( , 2020b)).We again take the θ = 90 pc data and now measure the intensity and line width along each line of sight.We sample the map using a hexagonal grid with spacing equal to the beam FWHM and consider only lines of sight where CO is detected over at least three 5 km s −1 channels.Then we show σ at this fixed spatial scale, which will be analogous to σ 2 /R for the clouds, as a function of Σ mol , also at the fixed scale.Again we show results for PHANGS-ALMA, now from Sun et al. (2020b) for reference.
Here we also plot the mass-weighted median and 16−84% range for each different environment.We plot typical error bars as the median fractional error associated with the NOEMA intensity map (14%) and linewidth map (18%) for Σ mol and σ, respectively; we use a bootstrap approach (with N = 1000 iterations) for the uncertainty in σ 2 /R, taking the difference between the CPROPS radius with and without deconvolution as a representative error bar for R.
The two approaches and all three panels tell a consistent story.The clouds seen in IC 342 have properties broadly consistent with other clouds seen throughout the local galaxy population.Some of the flattening in line width at low surface density reflects the 2× coarser velocity resolution of the IC 342 data compared to PHANGS-ALMA (in gray).If we consider the dynamical state in only the high surface density, Σ mol > 100 M ⊙ pc −2 clouds, where the CO line is typically well-resolved by NOEMA, then we find a median virial parameter α vir = 2KE/UE ≈ 1.2 and 1.5 in the arm and inter-arm regions, respectively.The median virial parameter is higher (α vir ≈ 2.0) within the bar region, but then drops again for clouds in the very center (α vir ≈ 1.4).In any case, we find a range of virial parameters across environments (the global 16th-84th percentile range goes from 0.8 to 4.0, with a median α vir ≈ 1.6).As the figure shows, these arm and inter-arm results are well in-line with the PHANGS-ALMA clouds.The bar, on the other hand, shows enhanced σ 2 /R compared to arm or inter-arm values at matched molecular gas surface density.For the center, we note two subtleties.First, the highest surface density regions associated with the most massive clouds or pixels in the lower panels of Fig. 3 actually have α vir lower than 2.0, indicated by lower σ at fixed Σ mol .In particular, the galaxy center itself appears, at our 90 pc resolution, roughly virialized, though we caution that 90 pc is quite coarse for such a study (e.g., see Schinnerer et al. 2006).Second, we caution that, as mentioned above, the center of IC 342 likely has a lower α CO than the Galactic value that we adopt (Israel 2020;Chiang et al. 2021).This will have the effect of moving the center points in all panels to the left, increasing their virial parameter by decreasing the surface density and thus the gravitational potential.
The surface density and line width increase moving from the inter-arm to the arm or bar, and then to the center regions.Weighting by mass, the median Σ mol is ∼80, 140, 160, and 1100 M ⊙ pc −2 for inter-arm, arm, bar, and center clouds, respectively, while the corresponding median line widths are 6.6, 7.6, 9.7 and 18.4 km s −1 .The beam-by-beam statistics show similar trends, with mass-weighted median Σ mol of 50, 80, 90, Fig. 3. Giant molecular clouds or associations in IC 342.The top left panel shows the location of the molecular clouds that we catalog at 90 pc resolution.The size of each circle corresponds to the deconvolved radius and the color indicates which dynamical environment the cloud is assigned to (red shows interarm clouds, green arm cloud, and blue center clouds).The top right and bottom left panels show the cloud properties in σ 2 /R vs. Σ mol space, in which clouds with a fixed virial parameter follow a diagonal line with a slope of unity in log-log space.Σ mol is the average surface density within the FWHM size of each cloud, Σ mol = M CO /(2πR 2 ).We show the IC 342 clouds for each environment and clouds at the same resolution and sensitivity, but better velocity resolution, from PHANGS-ALMA (Rosolowsky et al. 2021, and A. Hughes et al. in preparation).The bottom left panel replaces individual clouds with running medians.The bottom right panel adopts a beam-wise approach in which each independent line of sight at a fixed 90 pc resolution supplies a measurement of line width and surface density (Sun et al. 2018(Sun et al. , 2020b)).Again we compare the IC 342 points to those from a large sample of PHANGS-ALMA galaxies and here we mark the channel width and approximate sensitivity limit for the IC 342 data with dashed lines (for a Gaussian CO line, FWHM = 2.35σ, so we always have more than one channel across the FWHM of the emission line).The bottom row also indicates the mass-weighted median (position of the color circles) and 16−84% range (span of the color horizontal or vertical line) for surface densities and line widths.These refer to the corresponding horizontal or vertical axis, with an arbitrary positioning along the perpendicular direction.The mass-weighted averages tend to lie at higher values than the bulk of the individual measurements because much of the molecular gas mass resides in a few high-mass clouds or associations.All panels tell a consistent story: arm clouds show mildly enhanced surface density and line width compared to inter-arm clouds, and the center shows significant enhancements in surface density and line width.The elevated line widths in the center indicate high virial parameters suggesting clouds with additional contributions to their line widths.The deviations from self-gravity virialization would be even more extreme if we adopted a CO-to-H 2 conversion factor below the Galactic value used to construct these plots.The black crosses in the right panels show representative error bars for our measurements in IC 342, as explained in the main text.and 820 M ⊙ pc −1 and mass-weighted median line widths of 5.2, 6.2, 7.7, and 17.2 km s −1 in the inter-arm, arm, bar, and center regions.According to Pan et al. (2014), α CO could be 2-3 times lower in what we define as the center region (innermost ∼kpc), which would lead to more moderate center surface densities (∼500 and ∼300 M ⊙ pc −2 for GMCs and beam-by-beam, Fig. 4. Star formation efficiency per gravitational free-fall time in 1.5 kpc regions.Each hexagonal point shows the integrated molecular gas depletion time for a 1.5 kpc diameter region as a function of the mass-weighted average τ ff (see Eq. 4) for 90 pc resolution clouds in that region.A straight, diagonal line in this space, such as the dark red one, corresponds to a fixed ϵ ff , which could be expected if density variations represent the primary drivers of depletion time variations.We color the regions by galactocentric radius and calculate the implied ϵ ff for each region.We estimate a median ϵ ff of 0.45% with a 16−84% range of 0.33−0.71%,and we illustrate these with the solid red line and shaded pink region.For the central 1.5 kpc region, we illustrate the effect of switching from our adopted Galactic α CO to a starburst conversion factor.The black cross shows a representative error bar as explained in the main text.The plotted values can be found in Table A.3. respectively).A two-sided Kolmogorov-Smirnov test confirms that the arm-interarm differences are statistically significant (pvalue < 0.1%).At a similar level, the center is clearly different from the rest.The bar surface densities are not statistically distinct from the spiral arm values (p-value 36.8% for GMCs; 2.2% pixel-by-pixel), but bars do deviate from spiral arms in terms of σ 2 /R and σ (p-value < 0.1%).
Overall, these trends show an increase in surface density and line width moving from inter-arm to arm regions but with an approximately similar dynamical state.The center shows an increased surface density and line width.Many clouds in the central region, though not in the immediate center itself, show high line width relative to their surface density and evidence of broadening by forces other than self-gravity.A depressed α CO value would even emphasize this effect for the center.These results agree well with those found for larger populations of more distant galaxies in PHANGS-ALMA (Sun et al. 2020b;Rosolowsky et al. 2021).

Star formation efficiency per free-fall time
The rate at which stars form out of gas is often claimed to be primarily driven by gas density, which determines the local freefall time, τ ff (e.g., see reviews in Lada & Lada 2003;Kennicutt & Evans 2012;Federrath & Klessen 2013;Padoan et al. 2014).Therefore, models of star formation often treat τ ff as a characteristic timescale for star formation, and express the fraction of gas converted to stars per gravitational free-fall time via ϵ ff = SFR/(M mol /τ ff ), a.k.a. the "efficiency per free-fall time" (e.g., McKee & Ostriker 2007;Krumholz et al. 2019).Some work has suggested ϵ ff to be approximately constant across a wide range of scales and systems, so that density variations represent a main physical driver of variations in the molecular gas depletion time (e.g., see Krumholz & Thompson 2007;Krumholz et al. 2019).In other work, ϵ ff represents a key quantity but varies depending on the local physical state of the gas (e.g., Utomo et al. 2018;Sun et al. 2023).
We estimated the molecular cloud free-fall time τ ff in IC 342 from our high-resolution CO data.Our CPROPS analysis provides measurements of mass and size for each identified molecular cloud.From these measurements, we estimated the mean volume density via ρ = 3M/4πR 3 and the associated gravitational free-fall time via τ ff = 3π/(32Gρ).We then averaged over all clouds in each 1.5 kpc-diameter hexagonal region of IC 342, calculating the mass-weighted mean reciprocal of τ ff via where i is the index of clouds in that kpc-sized aperture.We follow Utomo et al. (2018) and Sun et al. (2022Sun et al. ( , 2023) ) in calculating the mass-weighted harmonic mean of τ ff , which ensures that the result is representative of the free-fall time for the entire cloud population in each hexagonal aperture.This ensures that ⟨τ ff ⟩ −1 is the correct quantity to compare to τ mol dep to test the hypothesis of fixed ϵ ff .We estimated typical error bars on our measurements following a bootstrap approach with N = 1000 iterations.Within each hexagonal aperture, we perturbed the surface density and radius of each cloud following a Gaussian with standard deviation equal to 14% of each cloud's surface density and 16% of its radius, which are the typical error bars that we have derived across our sample.In each iteration, we recomputed ρ and τ ff for each cloud with these perturbed values, and calculated the mass-weighted average of τ −1 ff within each hexagonal aperture.Based on these 1000 iterations, we obtained the fractional error on τ ff for each aperture, and computed the median across all apertures, which is 6%.For τ mol dep , we assume the same 14% typical uncertainty on Σ mol and take the median difference between the SFR calculated from WISE 4 and WISE 3 (17%) as a representative error bar.These combine to yield a typical error bar of 22% for τ mol dep .The median ⟨τ ff ⟩ in our sample of regions is 11 Myr with a 16−84% range of 8.5−14.8Myr.This is slightly longer than the median value of τ ff = 6.4 +3.5 −2.5 found for other galaxies at this scale by Rosolowsky et al. (2021), but still in agreement within 2σ given the measured galaxy-to-galaxy variation.For comparison, Kim et al. (2021), also using these data, made a statistical estimate of the cloud lifetime in IC 342 and found ∼20 Myr.This implies that clouds live ∼2τ ff in this galaxy, in good agreement with the results of Kruijssen et al. (2019), Chevance et al. (2020), andKim et al. (2022) for a broader sample of galaxies.
In addition to ⟨τ ff ⟩, we also calculated the integrated molecular gas depletion time, τ mol dep , in each region via Here, M mol and SFR are the integrated molecular gas mass and star formation rate in that aperture (measured from WISE 22 µm using Eq.3).A high τ mol dep corresponds to a low rate of star formation per unit molecular gas and implies that star formation will take a long time to consume the available gas reservoir.
We combined ⟨τ ff ⟩ and τ mol dep to estimate the mean efficiency per free-fall time in that aperture via ϵ ff = ⟨τ ff ⟩/τ mol dep .By using 1.5 kpc-sized apertures, we average over many clouds at different evolutionary stages, which makes our results robust to stochasticity and the evolutionary state of individual clouds.This approach follows similar work by Leroy et al. (2017), Utomo et al. (2018) and Schruba et al. (2019), and the hexagonal sampling scheme and measurements follow Sun et al. (2020a).
In Fig. 4, we compare ⟨τ ff ⟩ and τ mol dep for each 1.5 kpc region.Each hexagon symbol shows results for one aperture, with the points color-coded according to the deprojected galactocentric radius at the aperture center.The red line shows the expectation for τ dep = ⟨τ ff ⟩/ϵ ff at the median measured ϵ ff and the pink region shows the 16−84% range of estimated ϵ ff .If ϵ ff is indeed constant and if density variations drive variations in τ mol dep then we would expect the data to follow that red line.
For our 90 pc resolution cloud catalog and using CO (1−0) with a fixed conversion factor and our WISE 22 µm-based SFR estimates, we find a median of ϵ ff = 0.45% with a 16−84% range of 0.33−0.71%.This value is near the frequently invoked value 0.5% (e.g., McKee & Ostriker 2007) and broadly consistent with previous estimates in the literature (e.g., Vutisalchavakul et al. 2016;Barnes et al. 2017;Leroy et al. 2017;Utomo et al. 2018;Schruba et al. 2019).In detail, the value is towards the low end of the distribution for PHANGS-ALMA galaxies found by Utomo et al. (2018) and intermediate between their 0.7% mean value and the 0.3% found at 40 pc resolution for PAWS by Leroy et al. (2017).All of these values tend to be somewhat lower than those found for Galactic objects.
We can wonder if our IC 342 measurements support the idea of a fixed ϵ ff .On one hand, within the disk of the galaxy, that is, excluding the central aperture, there does not appear to be any significant relation between ⟨τ ff ⟩ and τ mol dep .The Spearman's rank correlation relating the two is about −0.14, which goes opposite the predicted direction for a constant ϵ ff , and the relationship has an insignificant p-value of 0.26.The disk of IC 342 appears to exhibit a similar value of ϵ ff to other cases, but any internal changes in τ mol dep do not show clear evidence of being driven by density variations.
On the other hand, the contrast between the central 1.5 kpc averaged measurement, indicated by a bright yellow hexagon, and the cloud of regions in the disk does offer some support for the ideas that ϵ ff is fixed and that density variations drive τ mol dep variations.Unfortunately, for the center, α CO complicates the picture.We have adopted a single fixed α CO for Fig. 4 but expect that the center, at least, shows a depressed α CO .For instance, Meier et al. (2011), Pan et al. (2014), Israel (2020), and Chiang et al. (2021), all find low α CO for the inner part of IC 342.This can have a dramatic effect on ϵ ff , which we illustrate by also plotting results for the central 1.5 kpc of IC 342 using a "starburst" α CO = 0.8 M ⊙ pc −2 (K km s −1 ) −1 conversion factor (Downes & Solomon 1998;Bolatto et al. 2013).We do not expect α CO to be necessarily this low in the center of IC 342, but this implies a difference of a factor of ∼5 between center and disk, which is probably realistic.For example, Pan et al. (2014) found a central α CO 2-3 below the Galactic value, while for the disk α CO was found to be typically two times above the Galactic value, amounting to a similar difference of a factor ∼5. Applying this five times lower α CO in the center implies a significantly higher ϵ ff for the central region, ∼4.2% instead of 0.3%.This suggests that ϵ ff varies between the disk and the central region and highlights the importance of accurate α CO to these estimates.

Discussion: Central molecular gas reservoir and nuclear star cluster
Figure 5 shows three zooms of the IC 342 center after deprojection of its inclination on the plane of sky.Several properties of IC 342's central region are similar to those in our Milky Way: the high molecular gas concentration, its geometry, and the massive nuclear star cluster (NSC, e.g., Schinnerer et al. 2003;Neumayer et al. 2020).A striking feature of the new CO map is the sharp increase in molecular gas surface density in the central ∼1 kpc (see Fig. 2).The average molecular mass surface density in the central 1.5 kpc measured at 90 pc resolution is about 600 M ⊙ pc −2 (for a Galactic α CO ), that is, a factor of 10× higher than in the immediate surroundings.Such sharp contrasts in CO surface brightness between center and disk are only seen in nearby galaxies that host a prominent stellar bar (see, e.g., CO atlas from PHANGS-ALMA, Leroy et al. 2021c).
The non-axisymmetric gas morphology in the center appear related to the underlying dynamical structure in IC 342.However, the existence of a stellar bar in IC 342 has been debated in the literature.Buta & McCall (1999) did not discuss a potential bar, but they carried out ellipse fitting of V and I band images that shows changes in position angle and ellipticity which we find suggestive of a large-scale bar component.Based on ellipse fitting of deep near infrared images, Fathi et al. (2009) reported a large-scale stellar bar of 307 ′′ (or ∼5.1 kpc) semi-major axis oriented roughly north-south, while Hernandez et al. (2005) found the bar to be not well-defined.Crosthwaite et al. (2000) reported evidence for a "fat" or "boxy" bar potential based on the analysis of their H I kinematic data.
We examined high-quality NIR imaging from Spitzer IRAC 3.6 µm and found clear evidence for the presence of a stellar bar with semi-major length of 120 ′′ (or ∼2 kpc) and an orientation of P.A. = 160 • (Appendix B.3).In addition to this photometric evidence, the presence of a bar is also suggested by CO morphology in the central 60 ′′ (see for a summary and sketch Meier & Turner 2005).The distribution of the CO line emission shows two slightly curved gas lanes ending towards the center in a ringlike shape which coincides with a ∼100 pc diameter star-forming pseudoring (or nuclear spiral), and shock tracers appear preferentially on the leading sides of the lanes close to where they connect to the ring (a morphology that is clearly revealed by our data, see Fig. 5.).In the whole PHANGS-ALMA CO morphology catalog, such a setup is only observed in barred galaxies (Stuber et al. 2023).Moreover, CO morphology suggests that it is a relatively weak bar, because the bar lanes are curved and do not extend all the way out to the bar ends (e.g., Athanassoula 1992;Comerón et al. 2009).The CO kinematics are also consistent with the gas responding to an underlying non-axisymmetric (bar) potential.The presence of bar lanes is also supported by observations of the magnetic field (Beck 2015).Finally, the NSC in IC 342 is at the same time extremely massive and compactconsistent with a picture where bar-driven gas inflow triggers in situ star formation that drives the stellar densities in NSCs to extremely high values (Neumayer et al. 2020;Pechetti et al. 2020).
In this paper, we have presented measurements showing elevated gas velocity dispersion within the center environment (Sect.3.1).The current central gas concentration and bar lanes, which are connected to a pseudoring (or nuclear spiral), can be interpreted as a consequence of that high velocity dispersion.Indeed, theoretical studies suggest that local gas properties such as the effective sound speed play a key role in determining the gas morphology within bars, beyond the existence or not of an inner Lindblad resonance (ILR; e.g., Englmaier & Gerhard 1997;Fig. 5. Zooms of the molecular surface density, peak temperature, and line width towards the galaxy center.The images have been deprojected from the galaxy inclination on the plane of sky, and then converted to kpc using a distance of 3.45 Mpc.The yellow lines show our spiral mask, while the red and white ellipses delimit the extent of the bar and center environment, respectively (see Appendix B.3 for the definition of these environment).Englmaier & Shlosman 2000;Sormani et al. 2018).The expectation is that warmer gas with high velocity dispersion (higher sound speed) will result in a more open nuclear spiral, which also critically depends on the central mass concentration.
Simulations show that the morphology of these central gas structures varies quickly over time (e.g., Emsellem et al. 2015;Renaud et al. 2015;Sormani et al. 2020).Therefore, the present configuration is probably not long-lived and was different in the past.We speculate that there might have been a stable nuclear ring associated with the bar in the past, which prevented gas inflow; the gas in this ring could have experienced dynamical and thermal heating due to feedback, giving way to the current open spiral morphology.This could signify that gas has only recently been able to reach the center, where it could contribute to the growth of the massive NSC.This process could have occurred intermittently, tied to rapid changes in the central gas morphology (i.e., from ring to open spiral and back again) that lead to periodic gas inflow to the center.In this case, bar-driven inflow episodes like this over the course of IC342's recent history would be responsible for building the NSC.Indeed, the nuclei of barred galaxies typically host much more elevated gas surface densities than unbarred galaxies (e.g., Sakamoto et al. 1999;Sun et al. 2020a).It should be kept in mind, however, that it is difficult to determine precisely the past evolutionary history from the present gas distribution (Renaud et al. 2015;Sormani et al. 2020).

Conclusions
We have presented a new, large, sensitive, ∼1000-pointing CO (1−0) line mosaic of the very nearby galaxy IC 342, carried out using the IRAM NOEMA and 30-m telescopes.The survey covers the inner ∼11 × 11 kpc 2 region, which corresponds to roughly 1.5 effective radii and extends out to near the radius where the galaxy becomes atomic gas dominated.With a linear resolution of ∼60 pc, our beam is about the size of a massive giant molecular cloud.We resolve the galaxy into a series of bright spiral arms, inter-arm filaments, and the well-known bright central molecular structure which is consistent with gas flowing along a stellar bar.The well-known bright starburst center of IC 342 appears prominent in all of our data products.
We applied the CPROPS algorithm to identify >600 massive molecular objects that are either molecular clouds or molecular associations.We did this at 90 pc resolution, so that IC 342 can be placed in the context of a large set of molecular cloud properties recently measured by PHANGS-ALMA.We also conducted a parallel "beam-wise" analysis in which we estimated the line width and molecular surface density for each independent line of sight where CO is detected over at least three consecutive channels.Both analyses highlight that the surface density and line width increase from the inter-arm region to the arm and bar region, and gas in the center of the galaxy has the highest surface density and line width.For clouds, we find a mass-weighted median Σ mol ≈ 80, 140, 160, and 1100 M ⊙ pc −2 , and σ CO ≈ 6.6, 7.6, 9.7 and 18.4 km s −1 for inter-arm, arm, bar, and center, respectively.Our beam-by-beam statistics yield similar results, Σ mol ≈ 50, 80, 90, and 820 M ⊙ pc −1 , and σ CO ≈ 5.2, 6.2, 7.7, and 17.2 km s −1 , respectively.A 2-3 times lower α CO in the center yields median surface densities of ∼500 and ∼300 M ⊙ pc −2 for GMCs and beam-by-beam, respectively.Moreover, for a Galactic α CO clouds appear to be in rough virial equilibrium in both the cloud-based and "beam-wise" analysis.The clouds in the bar and central region show elevated line widths, consistent with dynamical broadening by forces other than self-gravity and in good agreement with previous observations of galaxy centers using both techniques.If we had adopted a lower conversion factor, as is likely appropriate for the inner region, the dynamical state of the gas in the center would depart even more from simple virial equilibrium where turbulence balances self-gravity.
We used the cloud properties measured at 90 pc resolution to estimate the density, gravitational free-fall times, and star formation efficiency per free-fall time.We find median ⟨τ ff ⟩ ≈ 11 Myr and ϵ ff = 0.45% on average with a 16−84% range of 0.33−0.71%.This is in agreement with literature estimates of this quantity, though towards the low end of such estimates.The center shows shorter molecular gas depletion time and gravitational free-fall time.Whether it has the same ϵ ff compared to the disk depends sensitively on the adopted α CO .Combining these values with recent statistical estimates of the GMC lifetime in IC 342 implies that clouds in the galaxy live for typically ∼2 times the free-fall time.
Given the proximity of IC 342, the CO observations and GMC catalog presented in this paper could have a wide range of applications for studies of the nearby universe.The data, along with our data products, are publicly available.acknowledge funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 694343).JS acknowledges support by the Natural Sciences and Engineering Research Council of Canada (NSERC) through a Canadian Institute for Theoretical Astrophysics (CITA) National Fellowship.MC and JMDK gratefully acknowledge funding from the DFG through an Emmy Noether Grant (grant number KR4801/1-1), as well as from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme via the ERC Starting Grant MUSTANG (grant agreement number 714907).MC gratefully acknowledges funding from the DFG through an Emmy Noether Grant (grant number CH2137/1-1).COOL Research DAO is a Decentralised Autonomous Organisation supporting research in astrophysics aimed at uncovering our cosmic origins.MC, JK, and JMDK gratefully acknowledge funding from the German Research Foundation (DFG) through the DFG Sachbeihilfe (grant number KR4801/2-1).CE gratefully acknowledges funding from the Deutsche Forschungsgemeinschaft (DFG) Sachbeihilfe, grant number BI1546/3-1.KK gratefully acknowledges funding from the German Research Foundation (DFG) in the form of an Emmy Noether Research Group (grant number KR4598/2-1).RSK acknowledges financial support from the European Research Council (ERC) via the ERC Synergy Grant "ECOGAL" (grant 855130) as well as from the Heidelberg cluster of excellence EXC 2181 (Project-ID 390900948) "STRUCTURES" funded by the German Excellence Strategy.RSK also thanks DFG for support via the collaborative research center (SFB 881, Project-ID 138713538) "The Milky Way System" (subprojects A1, B1, B2, and B8).Notes. (a) The spectral resolutions and associated RMS noise are given for the best spectral resolution of our products, close to the WIDEX channel spacing.The channel spacing achived at the IRAM-30m was about 10 times smaller. (b) NOEMA had 8 antennas in 2016.With 12 antennas and the On-The-Fly observing mode, the same data quality would have been achieved with 26 hours of NOEMA observatory time.Notes.The columns correspond to the cloud coordinates, velocity centroid, deconvolved cloud radius, position angle of the cloud major axis, cloud eccentricity, line width, and cloud luminous and virial mass.The cloud eccentricity is calculated as e = 1 − σ 2 min,d /σ 2 maj,d .M CO is obtained as the cloud CO luminosity (flux in K km s −1 pc 2 ) multiplied by our assumed Galactic CO-to-H 2 conversion factor.The virial mass is estimated as Notes.The columns correspond to the aperture galactocentric radius, ensemble-averaged free-fall time, molecular gas and SFR surface densities, and mean depletion time within hexagonal apertures of 1.5 kpc diameter.
at each resolution, and a broad complete mask composed of the union of these masks, a set of moments (peak temperature, line integrated area, centroid velocity, line width, equivalent width) and their associated uncertainties.For the beam-by-beam analysis, we employ the CO intensity and line width derived consistently using the strict mask.The creation of these products followed the PHANGS standard procedure described in Leroy et al. (2021c), except that we worked at 5, 10, 20 km s −1 channel resolutions instead of 2.5 km s −1 used for the PHANGS-ALMA data.

Appendix B.2: Cloud catalog
To identify molecular clouds, we ran the CPROPS algorithm (Rosolowsky & Leroy 2006;Rosolowsky et al. 2021) on the 5 km s −1 data cube after convolving it to θ = 90 pc (FWHM) resolution.We used the implementation by Rosolowsky et al. (2021) and choose θ = 90 pc to exactly match the physical resolution of their analysis.This allows a fair comparison to the measurements of GMC properties for the subset of PHANGS-ALMA galaxies that have resolution ≤ 90 pc by Rosolowsky et al. (2021) and A. Hughes et al. (in preparation), though we do note that the PHANGS-ALMA cubes have velocity resolu-

Fig. 2 .
Fig. 2. NOEMA survey of CO (1−0) emission from IC 342.The distribution of the CO (1−0) peak intensity at the native angular resolution (top left) and the intensity-weighted mean CO velocity at 90 pc resolution (top right) as revealed by our new CO survey.The images show spiral arms, abundant inter-arm emission, and a velocity field that mostly reflects a regularly rotating gas disk.The bottom panels show the molecular gas traced by CO in the context of atomic gas and stars.In these panels, the CO and H I data are shown at a common ∼350 pc resolution.The bottom left panel shows a map of the total neutral gas surface density adding molecular gas surface density to atomic gas surface density from Chiang et al. (2021).The white ellipses show 1, 2, 3 × R e and the gray dashed ellipse indicates 1 R 25 .The lower right panel shows the azimuthally averaged mass surface density profiles for atomic gas, molecular gas, and stellar mass estimated from the near-infrared.The stellar mass distribution using HST imaging is consistent with a nuclear star cluster and an exponential disk(Carson et al. 2015), while the presented profile suffers from the resolution of the WISE data.Our NOEMA survey covers the inner molecule-dominated region where the H I emission is depressed, including the CO-bright center.

Table 1 .
Global properties of IC 342.

Table A .
2. Example catalog of GMC properties of IC 342.The full table is available online.

Table A .
3.Example table of free-fall time and depletion time in IC 342 (Fig.4).The full table is available online.