Wide-field CO isotopologue emission and the CO-to-H2_2 factor across the nearby spiral galaxy M101

Carbon monoxide (CO) emission is the most widely used tracer of the bulk molecular gas in the interstellar medium (ISM) in extragalactic studies. The CO-to-H$_2$ conversion factor, $\alpha_{\rm CO}$, links the observed CO emission to the total molecular gas mass. However, no single prescription perfectly describes the variation of $\alpha_{\rm CO}$ across all environments across galaxies as a function of metallicity, molecular gas opacity, line excitation, and other factors. Using resolved spectral line observations of CO and its isotopologues, we can constrain the molecular gas conditions and link them to a variation in the conversion factor. We present new IRAM 30-m 1mm and 3mm line observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O} across the nearby galaxy M101. Based on the CO isotopologue line ratios, we find that selective nucleosynthesis and opacity changes are the main drivers of the variation in the line emission across the galaxy. Furthermore, we estimated $\alpha_{\rm CO(1-0)}$ using different approaches, including (i) the dust mass surface density derived from far-IR emission as an independent tracer of the total gas surface density and (ii) LTE-based measurements using the optically thin $^{13}$CO(1-0) intensity. We find an average value of $\alpha_{\rm CO}=4.4{\pm}0.9\rm\,M_\odot\,pc^{-2}(K\,km\,s^{-1})^{-1}$ across the galaxy, with a decrease by a factor of 10 toward the 2 kpc central region. In contrast, we find LTE-based values are lower by a factor of 2-3 across the disk relative to the dust-based result. Accounting for $\alpha_{\rm CO}$ variations, we found significantly reduced molecular gas depletion time by a factor 10 in the galaxy's center. In conclusion, our result suggests implications for commonly derived scaling relations, such as an underestimation of the slope of the Kennicutt Schmidt law, if $\alpha_{\rm CO}$ variations are not accounted for.Comment: Accepted for publication in A&A, 25 pages, 15 figure

Calibrating mid-infrared emission as a tracer of obscured star formation on HII-region scales in the era of JWST

Measurements of the star formation activity on cloud scales are fundamental to uncovering the physics of the molecular cloud, star formation, and stellar feedback cycle in galaxies. Infrared (IR) emission from small dust grains and polycyclic aromatic hydrocarbons (PAHs) are widely used to trace the obscured component of star formation. However, the relation between these emission features and dust attenuation is complicated by the combined effects of dust heating from old stellar populations and an uncertain dust geometry with respect to heating sources. We use images obtained with NIRCam and MIRI as part of the PHANGS--JWST survey to calibrate dust emission at 21$\rm \mu m$, and the emission in the PAH-tracing bands at 3.3, 7.7, 10, and 11.3$\rm \mu m$ as tracers of obscured star formation. We analyse $\sim$ 20000 optically selected HII regions across 19 nearby star-forming galaxies, and benchmark their IR emission against dust attenuation measured from the Balmer decrement. We model the extinction-corrected H$\alpha$ flux as the sum of the observed H$\alpha$ emission and a term proportional to the IR emission, with $a_{IR}$ as the proportionality coefficient. A constant $a_{IR}$ leads to extinction-corrected H$\alpha$ estimates which agree with those obtained with the Balmer decrement with a scatter of $\sim$ 0.1 dex for all bands considered. Among these bands, 21$\rm \mu m$ emission is demonstrated to be the best tracer of dust attenuation. The PAH-tracing bands underestimate the correction for bright HII regions, since in these environments the ratio of PAH-tracing bands to 21$\rm \mu m$ decreases, signalling destruction of the PAH molecules. For fainter HII regions all bands suffer from an increasing contamination from the diffuse infrared background.Comment: accepted for publication in A&

The Physical Drivers and Observational Tracers of CO-to-H2 Conversion Factor Variations in Nearby Barred Galaxy Centers

The CO-to-H$_2$ conversion factor (\alpha_\rm{CO}) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower \alpha_\rm{CO} in the centers of some barred galaxies on kpc scales. To unveil the physical drivers of such variations, we obtained ALMA Band 3, 6, and 7 observations toward the inner 2 kpc of NGC 3627 and NGC 4321 tracing $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales. Our multi-line modeling and Bayesian likelihood analysis of these datasets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of \alpha_\rm{CO}. The central 300 pc nuclei in both galaxies show strong enhancement of temperature T_\rm{k}>100 K and density n_\rm{H_2}>10^3 cm$^{-3}$. Assuming a CO-to-H$_2$ abundance of $3\times10^{-4}$, we derive 4-15 times lower \alpha_\rm{CO} than the Galactic value across our maps, which agrees well with previous kpc-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of \alpha_\rm{CO} with low-J $^{12}$CO optical depths (\tau_\rm{CO}), as well as an anti-correlation with T_\rm{k}. The \tau_\rm{CO} correlation explains most of the \alpha_\rm{CO} variation in the three galaxy centers, whereas changes in T_\rm{k} influence \alpha_\rm{CO} to second order. Overall, the observed line width and $^{12}$CO/$^{13}$CO 2-1 line ratio correlate with \tau_\rm{CO} variation in these centers, and thus they are useful observational indicators for \alpha_\rm{CO} variation. We also test current simulation-based \alpha_\rm{CO} prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations.Comment: Accepted for publication in ApJ; 30 pages of main text + 3 appendice

Fuelling the nuclear ring of NGC 1097

Galactic bars can drive cold gas inflows towards the centres of galaxies. The gas transport happens primarily through the so-called bar ``dust lanes'', which connect the galactic disc at kpc scales to the nuclear rings at hundreds of pc scales much like two gigantic galactic rivers. Once in the ring, the gas can fuel star formation activity, galactic outflows, and central supermassive black holes. Measuring the mass inflow rates is therefore important to understanding the mass/energy budget and evolution of galactic nuclei. In this work, we use CO datacubes from the PHANGS-ALMA survey and a simple geometrical method to measure the bar-driven mass inflow rate onto the nuclear ring of the barred galaxy NGC~1097. The method assumes that the gas velocity in the bar lanes is parallel to the lanes in the frame co-rotating with the bar, and allows one to derive the inflow rates from sufficiently sensitive and resolved position-position-velocity diagrams if the bar pattern speed and galaxy orientations are known. We find an inflow rate of $\dot{M}=(3.0 \pm 2.1)\, \rm M_\odot\, yr^{-1}$ averaged over a time span of 40 Myr, which varies by a factor of a few over timescales of $\sim$10 Myr. Most of the inflow appears to be consumed by star formation in the ring which is currently occurring at a rate of ${\rm SFR}\simeq~1.8$-$2 \rm M_\odot\, yr^{-1}$, suggesting that the inflow is causally controlling the star formation rate in the ring as a function of time.Comment: Accepted in MNRA

Kinematic analysis of the super-extended HI disk of the nearby spiral galaxy M 83

Funding: CE, FB, AB, IB, JdB and JP acknowledge funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No.726384/Empire). TGW acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 694343). JMDK gratefully acknowledges funding 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). SCOG acknowledges funding from the European Research Council via the ERC Synergy Grant “ECOGAL – Understanding our Galactic ecosystem: From the disk of the Milky Way to the formation sites of stars and planets” (project ID 855130). WJGdB received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 882793 ‘MeerGas’).We present new HI observations of the nearby massive spiral galaxy M83, taken with the VLA at 21″ angular resolution (≈500 pc) of an extended (1.5 deg2) 10-point mosaic combined with GBT single dish data. We study the super-extended HI disk of M83 (∼50 kpc in radius), in particular disc kinematics, rotation and the turbulent nature of the atomic interstellar medium. We define distinct regions in the outer disk (rgal > central optical disk), including ring, southern area, and southern and northern arm. We examine HI gas surface density, velocity dispersion and non-circular motions in the outskirts, which we compare to the inner optical disk. We find an increase of velocity dispersion (σv) towards the pronounced HI ring, indicative of more turbulent HI gas. Additionally, we report over a large galactocentric radius range (until rgal ∼ 50 kpc) that σv is slightly larger than thermal (i.e. > 8 km s-1). We find that a higher star formation rate (as traced by FUV emission) is not always necessarily associated with a higher HI velocity dispersion, suggesting that radial transport could be a dominant driver for the enhanced velocity dispersion. We further find a possible branch that connects the extended HI disk to the dwarf irregular galaxy UGCA365, that deviates from the general direction of the northern arm. Lastly, we compare mass flow rate profiles (based on 2D and 3D tilted ring models) and find evidence for outflowing gas at rgal ∼ 2 kpc, inflowing gas at rgal ~ 5.5 kpc and outflowing gas at rgal ~ 14 kpc. We caution that mass flow rates are highly sensitive to the assumed kinematic disk parameters, in particular, to the inclination.Publisher PDFPeer reviewe

A Cloud-Scale View of the Molecular Gas Disk in the Whirlpool Galaxy and Beyond

The nearby galaxy M51 (also known as the Whirlpool galaxy) hosts an iconic grand-design spiral pattern and both IRAM facilities conducted the first cloud-scale (∼50 pc resolution) survey of the molecular gas reservoir across a the disk of a massive star-forming galaxy (PAWS, PdBI+30m Arcsecond Whirlpool galaxy Survey) using the CO(1-0) line emission. PAWS showed that the various properties of the giant molecular cloud (GMC) population vary with galactic environment (center/bar, spiral arms, inter-arm). Recent observations of a ∼1000 pointing mosaic of the nearby late-type spiral galaxy IC342 using NOEMA resolved its GMC population at ∼70 pc resolution and find consistent trends. Investigation of the dense molecular gas phase at cloud-scales using tracers such as HCN(1-0) confirms the trends seen in kpc-scale surveys, namely that the dense gas star formation efficiency inn general apparently anticorrelates with the inferred dense gas fraction. Multi-line studies of the molecular gas in the galactic disks of nearby galaxies such as the ongoing large NOEMA+30m program to map the dense molecular gas phase in GMCs in the central part of M51 will allow for gaining new insights of the properties of this important molecular gas phase

Calibration of hybrid resolved star formation rate recipes based on PHANGS–MUSE H

Mapping star-formation rates (SFR) within galaxies is key to unveiling their assembly and evolution. Calibrations exist for computing the SFR from a combination of ultraviolet and infrared bands for galaxies as integrated systems, but their applicability to sub-galactic (kiloparsec) scales remains largely untested. We used integral field spectroscopy of 19 nearby (D   − 9.9). Moreover, they lead to negligible bias and < 0.16 dex scatter when compared to our reference attenuation-corrected SFR from Hα. These calibrations prove reliable as a function of physical scale. In particular, they agree within 10% with the attenuation corrections computed from the Balmer decrement on 100 pc scales. Despite small quantitative differences, our calibrations are also applicable to integrated galaxy scales probed by the MaNGA survey, but with a larger scatter (up to 0.22 dex). Observations with JWST open up the possibility to calibrate these relations in nearby galaxies with cloud-scale (∼100 pc) resolution mid-IR imaging

Wide-field CO isotopologue emission and the CO-to-H

Carbon monoxide (CO) emission constitutes the most widely used tracer of the bulk molecular gas in the interstellar medium (ISM) in extragalactic studies. The CO-to-H2 conversion factor, α12CO(1−0), links the observed CO emission to the total molecular gas mass. However, no single prescription perfectly describes the variation of α12CO(1−0) across all environments within and across galaxies as a function of metallicity, molecular gas opacity, line excitation, and other factors. Using spectral line observations of CO and its isotopologues mapped across a nearby galaxy, we can constrain the molecular gas conditions and link them to a variation in α12CO(1−0). Here, we present new, wide-field (10 × 10 arcmin2) IRAM 30-m telescope 1 mm and 3 mm line observations of 12CO, 13CO, and C18O across the nearby, grand-design, spiral galaxy M101. From the CO isotopologue line ratio analysis alone, we find that selective nucleosynthesis and changes in the opacity are the main drivers of the variation in the line emission across the galaxy. In a further analysis step, we estimated α12CO(1−0) using different approaches, including (i) via the dust mass surface density derived from far-IR emission as an independent tracer of the total gas surface density and (ii) local thermal equilibrium (LTE) based measurements using the optically thin 13CO(1–0) intensity. We find an average value of ⟨α12CO(1 − 0)⟩ = 4.4  ±  0.9 M⊙ pc−2 (K km s−1)−1 across the disk of the galaxy, with a decrease by a factor of 10 toward the 2 kpc central region. In contrast, we find LTE-based α12CO(1−0) values are lower by a factor of 2–3 across the disk relative to the dust-based result. Accounting for α12CO(1−0) variations, we found significantly reduced molecular gas depletion time by a factor 10 in the galaxy’s center. In conclusion, our result suggests implications for commonly derived scaling relations, such as an underestimation of the slope of the Kennicutt Schmidt law, if α12CO(1−0) variations are not accounted for