Search for: All records

Context.Mapping molecular line emission beyond the bright low-JCO transitions is still challenging in extragalactic studies, even with the latest generation of (sub-)millimetre interferometers, such as ALMA and NOEMA.

Aims.We summarise and test a spectral stacking method that has been used in the literature to recover low-intensity molecular line emission, such as HCN(1−0), HCO⁺(1−0), and even fainter lines in external galaxies. The goal is to study the capabilities and limitations of the stacking technique when applied to imaged interferometric observations.

Methods.The core idea of spectral stacking is to align spectra of the low S/N spectral lines to a known velocity field calculated from a higher S/N line expected to share the kinematics of the fainter line (e.g. CO(1−0) or 21 cm emission). Then these aligned spectra can be coherently averaged to produce potentially high S/N spectral stacks. Here we used imaged simulated interferometric and total power observations at different S/N levels, based on real CO observations.

Results.For the combined interferometric and total power data, we find that the spectral stacking technique is capable of recovering the integrated intensities even at low S/N levels across most of the region where the high S/N prior is detected. However, when stacking interferometer-only data for low S/N emission, the stacks can miss up to 50% of the emission from the fainter line.

Conclusions.A key result of this analysis is that the spectral stacking method is able to recover the true mean line intensities in low S/N cubes and to accurately measure the statistical significance of the recovered lines. To facilitate the application of this technique we provide a public Python package, called PYSTACKER.

 

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-H 2 conversion factor, α 12 CO(1−0) , links the observed CO emission to the total molecular gas mass. However, no single prescription perfectly describes the variation of α 12 CO(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 α 12 CO(1−0) . Here, we present new, wide-field (10 × 10 arcmin 2 ) IRAM 30-m telescope 1 mm and 3 mm line observations of 12 CO, 13 CO, and C 18 O 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 α 12 CO(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 13 CO(1–0) intensity. We find an average value of ⟨ α 12 CO(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 α 12 CO(1−0) values are lower by a factor of 2–3 across the disk relative to the dust-based result. Accounting for α 12 CO(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 α 12 CO(1−0) variations are not accounted for. 

Determining how the galactic environment, especially the high gas densities and complex dynamics in bar-fed galaxy centers, alters the star formation efficiency (SFE) of molecular gas is critical to understanding galaxy evolution. However, these same physical or dynamical effects also alter the emissivity properties of CO, leading to variations in the CO-to-H₂conversion factor (α_CO) that impact the assessment of the gas column densities and thus of the SFE. To address such issues, we investigate the dependence ofα_COon the local CO velocity dispersion at 150 pc scales using a new set of dust-basedα_COmeasurements and propose a newα_COprescription that accounts for CO emissivity variations across galaxies. Based on this prescription, we estimate the SFE in a sample of 65 galaxies from the PHANGS–Atacama Large Millimeter/submillimeter Array survey. We find increasing SFE toward high-surface-density regions like galaxy centers, while using a constant or metallicity-basedα_COresults in a more homogeneous SFE throughout the centers and disks. Our prescription further reveals a mean molecular gas depletion time of 700 Myr in the centers of barred galaxies, which is overall three to four times shorter than in nonbarred galaxy centers or the disks. Across the galaxy disks, the depletion time is consistently around 2–3 Gyr, regardless of the choice ofα_COprescription. All together, our results suggest that the high level of star formation activity in barred centers is not simply due to an increased amount of molecular gas, but also to an enhanced SFE compared to nonbarred centers or disk regions.

 

Abstract The CO-to-H 2 conversion factor ( α 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 α CO in the centers of some barred galaxies on kiloparsec scales. To unveil the physical drivers of such variations, we obtained Atacama Large Millimeter/submillimeter Array bands (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 multiline modeling and Bayesian likelihood analysis of these data sets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of α CO . The central 300 pc nuclei in both galaxies show strong enhancement of temperature T k ≳ 100 K and density n H 2 > 10 3 cm −3 . Assuming a CO-to-H 2 abundance of 3 × 10 −4 , we derive 4–15 times lower α CO than the Galactic value across our maps, which agrees well with previous kiloparsec-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of α CO with low- J 12 CO optical depths ( τ CO ), as well as an anticorrelation with T k . The τ CO correlation explains most of the α CO variation in the three galaxy centers, whereas changes in T k influence α CO to second order. Overall, the observed line width and 12 CO/ 13 CO 2–1 line ratio correlate with τ CO variation in these centers, and thus they are useful observational indicators for α CO variation. We also test current simulation-based α CO prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations. 

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 on to 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 time-scales of ∼10 Myr. Most of the inflow appears to be consumed by star formation in the ring, which is currently occurring at a star formation rate (SFR) of $\simeq\!1.8\!-\!2 \, \rm M_\odot \, yr^{-1}$, suggesting that the inflow is causally controlling the SFR in the ring as a function of time.

 

We use new HCN(1–0) data from the ACA Large-sample Mapping Of Nearby galaxies in Dense gas (ALMOND) survey to trace the kpc-scale molecular gas density structure and CO(2–1) data from the Physics at High Angular resolution in Nearby GalaxieS–Atacama Large Millimeter/submillimeter Array (PHANGS–ALMA) to trace the bulk molecular gas across 25 nearby star-forming galaxies. At 2.1 kpc scale, we measure the density-sensitive HCN/CO line ratio and the star formation rate (SFR)/HCN ratio to trace the star formation efficiency in the denser molecular medium. At 150 pc scale, we measure structural and dynamical properties of the molecular gas via CO(2–1) line emission, which is linked to the lower resolution data using an intensity-weighted averaging method. We find positive correlations (negative) of HCN/CO (SFR/HCN) with the surface density, the velocity dispersion, and the internal turbulent pressure of the molecular gas. These observed correlations agree with expected trends from turbulent models of star formation, which consider a single free-fall time gravitational collapse. Our results show that the kpc-scale HCN/CO line ratio is a powerful tool to trace the 150 pc scale average density distribution of the molecular clouds. Lastly, we find systematic variations of the SFR/HCN ratio with cloud-scale molecular gas properties, which are incompatible with a universal star formation efficiency. Overall, these findings show that mean molecular gas density, molecular cloud properties, and star formation are closely linked in a coherent way, and observations of density-sensitive molecular gas tracers are a useful tool to analyse these variations, linking molecular gas physics to stellar output across galaxy discs.

 

Abstract We compare mid-infrared (mid-IR), extinction-corrected H α , and CO (2–1) emission at 70–160 pc resolution in the first four PHANGS–JWST targets. We report correlation strengths, intensity ratios, and power-law fits relating emission in JWST’s F770W, F1000W, F1130W, and F2100W bands to CO and H α . At these scales, CO and H α each correlate strongly with mid-IR emission, and these correlations are each stronger than the one relating CO to H α emission. This reflects that mid-IR emission simultaneously acts as a dust column density tracer, leading to a good match with the molecular-gas-tracing CO, and as a heating tracer, leading to a good match with the H α . By combining mid-IR, CO, and H α at scales where the overall correlation between cold gas and star formation begins to break down, we are able to separate these two effects. We model the mid-IR above I ν = 0.5 MJy sr −1 at F770W, a cut designed to select regions where the molecular gas dominates the interstellar medium (ISM) mass. This bright emission can be described to first order by a model that combines a CO-tracing component and an H α -tracing component. The best-fitting models imply that ∼50% of the mid-IR flux arises from molecular gas heated by the diffuse interstellar radiation field, with the remaining ∼50% associated with bright, dusty star-forming regions. We discuss differences between the F770W, F1000W, and F1130W bands and the continuum-dominated F2100W band and suggest next steps for using the mid-IR as an ISM tracer. 

Abstract We measure the low- J CO line ratios R 21 ≡ CO (2–1)/CO (1–0), R 32 ≡ CO (3–2)/CO (2–1), and R 31 ≡CO (3–2)/CO (1–0) using whole-disk CO maps of nearby galaxies. We draw CO (2–1) from PHANGS-ALMA, HERACLES, and follow-up IRAM surveys; CO (1–0) from COMING and the Nobeyama CO Atlas of Nearby Spiral Galaxies; and CO (3–2) from the James Clerk Maxwell Telescope Nearby Galaxy Legacy Survey and Atacama Pathfinder Experiment Large APEX Sub-Millimetre Array mapping. All together, this yields 76, 47, and 29 maps of R 21 , R 32 , and R 31 at 20″ ∼ 1.3 kpc resolution, covering 43, 34, and 20 galaxies. Disk galaxies with high stellar mass, log ( M ⋆ / M ⊙ ) = 10.25 – 11 , and star formation rate (SFR) = 1–5 M ⊙ yr −1 , dominate the sample. We find galaxy-integrated mean values and a 16%–84% range of R 21 = 0.65 (0.50–0.83), R 32 = 0.50 (0.23–0.59), and R 31 = 0.31 (0.20–0.42). We identify weak trends relating galaxy-integrated line ratios to properties expected to correlate with excitation, including SFR/ M ⋆ and SFR/ L CO . Within galaxies, we measure central enhancements with respect to the galaxy-averaged value of ∼ 0.18 − 0.14 + 0.09 dex for R 21 , 0.27 − 0.15 + 0.13 dex for R 31 , and 0.08 − 0.09 + 0.11 dex for R 32 . All three line ratios anticorrelate with galactocentric radius and positively correlate with the local SFR surface density and specific SFR, and we provide approximate fits to these relations. The observed ratios can be reasonably reproduced by models with low temperature, moderate opacity, and moderate densities, in good agreement with expectations for the cold interstellar medium. Because the line ratios are expected to anticorrelate with the CO (1–0)-to-H 2 conversion factor, α CO 1 − 0 , these results have general implications for the interpretation of CO emission from galaxies. 

Large-scale bars can fuel galaxy centers with molecular gas, often leading to the development of dense ringlike structures where intense star formation occurs, forming a very different environment compared to galactic disks. We pair ∼0.″3 (30 pc) resolution new JWST/MIRI imaging with archival ALMA CO(2–1) mapping of the central ∼5 kpc of the nearby barred spiral galaxy NGC 1365 to investigate the physical mechanisms responsible for this extreme star formation. The molecular gas morphology is resolved into two well-known bright bar lanes that surround a smooth dynamically cold gas disk (R_gal∼ 475 pc) reminiscent of non-star-forming disks in early-type galaxies and likely fed by gas inflow triggered by stellar feedback in the lanes. The lanes host a large number of JWST-identified massive young star clusters. We find some evidence for temporal star formation evolution along the ring. The complex kinematics in the gas lanes reveal strong streaming motions and may be consistent with convergence of gas streamlines expected there. Indeed, the extreme line widths are found to be the result of inter-“cloud” motion between gas peaks;ScousePydecomposition reveals multiple components with line widths of 〈σ_CO,scouse〉 ≈ 19 km s⁻¹and surface densities of$〈{\mathrm{\Sigma }}_{{\mathrm{H}}_2,\mathrm{scouse}}〉\approx 800M_{\odot }{\mathrm{pc}}^{-2}$, similar to the properties observed throughout the rest of the central molecular gas structure. Tailored hydrodynamical simulations exhibit many of the observed properties and imply that the observed structures are transient and highly time-variable. From our study of NGC 1365, we conclude that it is predominantly the high gas inflow triggered by the bar that is setting the star formation in its CMZ.

 

We present a high-resolution view of bubbles within the Phantom Galaxy (NGC 628), a nearby (∼10 Mpc), star-forming (∼2M_⊙yr⁻¹), face-on (i∼ 9°) grand-design spiral galaxy. With new data obtained as part of the Physics at High Angular resolution in Nearby GalaxieS (PHANGS)-JWST treasury program, we perform a detailed case study of two regions of interest, one of which contains the largest and most prominent bubble in the galaxy (the Phantom Void, over 1 kpc in diameter), and the other being a smaller region that may be the precursor to such a large bubble (the Precursor Phantom Void). When comparing to matched-resolution Hαobservations from the Hubble Space Telescope, we see that the ionized gas is brightest in the shells of both bubbles, and is coincident with the youngest (∼1 Myr) and most massive (∼10⁵M_⊙) stellar associations. We also find an older generation (∼20 Myr) of stellar associations is present within the bubble of the Phantom Void. From our kinematic analysis of the HI, H₂(CO), and Hiigas across the Phantom Void, we infer a high expansion speed of around 15 to 50 km s⁻¹. The large size and high expansion speed of the Phantom Void suggest that the driving mechanism is sustained stellar feedback due to multiple mechanisms, where early feedback first cleared a bubble (as we observe now in the Precursor Phantom Void), and since then supernovae have been exploding within the cavity and have accelerated the shell. Finally, comparison to simulations shows a striking resemblance to our JWST observations, and suggests that such large-scale, stellar-feedback-driven bubbles should be common within other galaxies.