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. 

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 report CO(5 → 4) and CO(6 → 5) line observations in the dusty starbursting galaxy CRLE (z= 5.667) and the main-sequence (MS) galaxy HZ10 (z= 5.654) with the Northern Extended Millimeter Array. CRLE is the most luminousz> 5 starburst in the COSMOS field and HZ10 is the most gas-rich “normal” galaxy currently known atz> 5. We find line luminosities for CO(5 → 4) and CO(6 → 5) of (4.9 ± 0.5) and (3.8 ± 0.4) × 10¹⁰K km s⁻¹pc²for CRLE and upper limits of < 0.76 and < 0.60 × 10¹⁰K km s⁻¹pc²for HZ10, respectively. The CO excitation of CRLE appears comparable to otherz> 5 dusty star-forming galaxies. For HZ10, these line luminosity limits provide the first significant constraints of this kind for an MS galaxy atz> 5. We find the upper limit of$L_{5\to 4}^{\prime }/L_{2\to 1}^{\prime }$in HZ10 could be similar to the average value for MS galaxies aroundz≈ 1.5, suggesting that MS galaxies with comparable gas excitation may already have existed one billion years after the Big Bang. For CRLE we determine the most likely values for the H₂density, kinetic temperature, and dust temperature based on excitation modeling of the CO line ladder. We also derive a total gas mass of (7.1 ± 1.3) × 10¹⁰M_⊙. Our findings provide some of the currently most detailed constraints on the gas excitation that sets the conditions for star formation in a galaxy protocluster environment atz> 5.

 

We measure the molecular gas environment near recent (<100 yr old) supernovae (SNe) using ∼1″ or ≤150 pc resolution CO (2–1) maps from the PHANGS–Atacama Large Millimeter/submillimeter Array (ALMA) survey of nearby star-forming galaxies. This is arguably the first such study to approach the scales of individual massive molecular clouds (M_mol≳ 10^5.3M_⊙). Using the Open Supernova Catalog, we identify 63 SNe within the PHANGS–ALMA footprint. We detect CO (2–1) emission near ∼60% of the sample at 150 pc resolution, compared to ∼35% of map pixels with CO (2–1) emission, and up to ∼95% of the SNe at 1 kpc resolution, compared to ∼80% of map pixels with CO (2–1) emission. We expect the ∼60% of SNe within the same 150 pc beam, as a giant molecular cloud will likely interact with these clouds in the future, consistent with the observation of widespread SN–molecular gas interaction in the Milky Way, while the other ∼40% of SNe without strong CO (2–1) detections will deposit their energy in the diffuse interstellar medium, perhaps helping drive large-scale turbulence or galactic outflows. Broken down by type, we detect CO (2–1) emission at the sites of ∼85% of our 9 stripped-envelope SNe (SESNe), ∼40% of our 34 Type II SNe, and ∼35% of our 13 Type Ia SNe, indicating that SESNe are most closely associated with the brightest CO (2–1) emitting regions in our sample. Our results confirm that SN explosions are not restricted to only the densest gas, and instead exert feedback across a wide range of molecular gas densities.

 

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 The CO-to-H 2 conversion factor ( α CO ) is critical to studying molecular gas and star formation in galaxies. The value of α CO has been found to vary within and between galaxies, but the specific environmental conditions that cause these variations are not fully understood. Previous observations on ~kiloparsec scales revealed low values of α CO in the centers of some barred spiral galaxies, including NGC 3351. We present new Atacama Large Millimeter/submillimeter Array Band 3, 6, and 7 observations of 12 CO, 13 CO, and C 18 O lines on 100 pc scales in the inner ∼2 kpc of NGC 3351. Using multiline radiative transfer modeling and a Bayesian likelihood analysis, we infer the H 2 density, kinetic temperature, CO column density per line width, and CO isotopologue abundances on a pixel-by-pixel basis. Our modeling implies the existence of a dominant gas component with a density of 2–3 × 10 3 cm −3 in the central ∼1 kpc and a high temperature of 30–60 K near the nucleus and near the contact points that connect to the bar-driven inflows. Assuming a CO/H 2 abundance of 3 × 10 −4 , our analysis yields α CO ∼ 0.5–2.0 M ⊙ (K km s −1 pc 2 ) −1 with a decreasing trend with galactocentric radius in the central ∼1 kpc. The inflows show a substantially lower α CO ≲ 0.1 M ⊙ (K km s −1 pc 2 ) −1 , likely due to lower optical depths caused by turbulence or shear in the inflows. Over the whole region, this gives an intensity-weighted α CO of ∼1.5 M ⊙ (K km s −1 pc 2 ) −1 , which is similar to previous dust-modeling-based results at kiloparsec scales. This suggests that low α CO on kiloparsec scales in the centers of some barred galaxies may be due to the contribution of low-optical-depth CO emission in bar-driven inflows. 

We use PHANGS–James Webb Space Telescope (JWST) data to identify and classify 1271 compact 21μm sources in four nearby galaxies using MIRI F2100W data. We identify sources using a dendrogram-based algorithm, and we measure the background-subtracted flux densities for JWST bands from 2 to 21μm. Using the spectral energy distribution (SED) in JWST and HST bands plus ALMA and MUSE/VLT observations, we classify the sources by eye. Then we use this classification to define regions in color–color space and so establish a quantitative framework for classifying sources. We identify 1085 sources as belonging to the ISM of the target galaxies with the remainder being dusty stars or background galaxies. These 21μm sources are strongly spatially associated with Hiiregions (>92% of sources), while 74% of the sources are coincident with a stellar association defined in the HST data. Using SED fitting, we find that the stellar masses of the 21μm sources span a range of 10²–10⁴M_⊙with mass-weighted ages down to 2 Myr. There is a tight correlation between attenuation-corrected Hαand 21μm luminosity forL_ν,F2100W> 10¹⁹W Hz⁻¹. Young embedded source candidates selected at 21μm are found below this threshold and haveM_⋆< 10³M_⊙.

 

Active galactic nucleus (AGN) feedback is postulated as a key mechanism for regulating star formation within galaxies. Studying the physical properties of the outflowing gas from AGNs is thus crucial for understanding the coevolution of galaxies and supermassive black holes. Here we report 55 pc resolution ALMA neutral atomic carbon [Ci]³P₁−³P₀observations toward the central 1 kpc of the nearby Type 2 Seyfert galaxy NGC 1068, supplemented by 55 pc resolution CO(J= 1−0) observations. We find that [Ci] emission within the central kiloparsec is strongly enhanced by a factor of >5 compared to the typical [Ci]/CO intensity ratio of ∼0.2 for nearby starburst galaxies (in units of brightness temperature). The most [Ci]-enhanced gas (ratio > 1) exhibits a kiloparsec-scale elongated structure centered at the AGN that matches the known biconical ionized gas outflow entraining molecular gas in the disk. A truncated, decelerating bicone model explains well the kinematics of the elongated structure, indicating that the [Ci] enhancement is predominantly driven by the interaction between the ISM in the disk and the highly inclined ionized gas outflow (which is likely driven by the radio jet). Our results strongly favor the “CO dissociation scenario” rather than the “in situ C formation” one, which prefers a perfect bicone geometry. We suggest that the high-[Ci]/CO intensity ratio gas in NGC 1068 directly traces ISM in the disk that is currently dissociated and entrained by the jet and the outflow, i.e., the “negative” effect of the AGN feedback.