Search for: All records

We use high-resolution maps of the molecular interstellar medium (ISM) in the centres of 86 nearby galaxies from the millimetre-Wave Interferometric Survey of Dark Object Masses (WISDOM) and Physics at High Angular Resolution in Nearby GalaxieS (PHANGS) surveys to investigate the physical mechanisms setting the morphology of the ISM at molecular cloud scales. We show that early-type galaxies tend to have smooth, regular molecular gas morphologies, while the ISM in spiral galaxy bulges is much more asymmetric and clumpy when observed at the same spatial scales. We quantify these differences using non-parametric morphology measures (Asymmetry, Smoothness, and Gini), and compare these measurements with those extracted from idealized galaxy simulations. We show that the morphology of the molecular ISM changes systematically as a function of various large-scale galaxy parameters, including galaxy morphological type, stellar mass, stellar velocity dispersion, effective stellar mass surface density, molecular gas surface density, star formation efficiency, and the presence of a bar. We perform a statistical analysis to determine which of these correlated parameters best predicts the morphology of the ISM. We find the effective stellar mass surface (or volume) density to be the strongest predictor of the morphology of the molecular gas, while star formation and bars maybe be important secondary drivers. We find that gas self-gravity is not the dominant process shaping the morphology of the molecular gas in galaxy centres. Instead effects caused by the depth of the potential well, such as shear, suppression of stellar spiral density waves, and/or inflow, affect the ability of the gas to fragment.

 

ABSTRACT In the centres of the Milky Way and M83, the global environmental properties thought to control star formation are very similar. However, M83’s nuclear star formation rate (SFR), as estimated by synchrotron and H α emission, is an order of magnitude higher than the Milky Way’s. To understand the origin of this difference we use ALMA observations of HCN (1 − 0) and HCO+ (1 − 0) to trace the dense gas at the size scale of individual molecular clouds (0.54 arcsec, 12 pc) in the inner ∼500 pc of M83, and compare this to gas clouds at similar resolution and galactocentric radius in the Milky Way. We find that both the overall gas distribution and the properties of individual clouds are very similar in the two galaxies, and that a common mechanism may be responsible for instigating star formation in both circumnuclear rings. Given the considerable similarity in gas properties, the most likely explanation for the order of magnitude difference in SFR is time variability, with the Central Molecular Zone (CMZ) currently being at a more quiescent phase of its star formation cycle. We show M83’s SFR must have been an order of magnitude higher 5–7 Myr ago. M83’s ‘starburst’ phase was highly localized, both spatially and temporally, greatly increasing the feedback efficiency and ability to drive galactic-scale outflows. This highly dynamic nature of star formation and feedback cycles in galaxy centres means (i) modelling and interpreting observations must avoid averaging over large spatial areas or time-scales, and (ii) understanding the multiscale processes controlling these cycles requires comparing snapshots of a statistical sample of galaxies in different evolutionary stages. 

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. 

In the hierarchical view of star formation, giant molecular clouds (GMCs) undergo fragmentation to form small-scale structures made up of stars and star clusters. Here we study the connection between young star clusters and cold gas across a range of extragalactic environments by combining the high resolution (1″) PHANGS–ALMA catalogue of GMCs with the star cluster catalogues from PHANGS–HST. The star clusters are spatially matched with the GMCs across a sample of 11 nearby star-forming galaxies with a range of galactic environments (centres, bars, spiral arms, etc.). We find that after 4 − 6 Myr the star clusters are no longer associated with any gas clouds. Additionally, we measure the autocorrelation of the star clusters and GMCs as well as their cross-correlation to quantify the fractal nature of hierarchical star formation. Young (≤10 Myr) star clusters are more strongly autocorrelated on kpc and smaller spatial scales than the $\gt \, 10$ Myr stellar populations, indicating that the hierarchical structure dissolves over time.

 

ABSTRACT Previous work has argued that atomic gas mass estimates of galaxies from 21-cm H i emission are systematically low due to a cold opaque atomic gas component. If true, this opaque component necessitates a $\sim 35{{\ \rm per\ cent}}$ correction factor relative to the mass from assuming optically thin H i emission. These mass corrections are based on fitting H i spectra with a single opaque component model that produces a distinct ‘top-hat’ shaped line profile. Here, we investigate this issue using deep, high spectral resolution H i VLA observations of M31 and M33 to test if these top-hat profiles are instead superpositions of multiple H i components along the line of sight. We fit both models and find that ${\gt}80{{\ \rm per\ cent}}$ of the spectra strongly prefer a multicomponent Gaussian model while ${\lt}2{{\ \rm per\ cent}}$ prefer the single opacity-corrected component model. This strong preference for multiple components argues against previous findings of lines of sight dominated by only cold H i. Our findings are enabled by the improved spectral resolution (0.42 ${\rm km\, s^{-1}}$), whereas coarser spectral resolution blends multiple components together. We also show that the inferred opaque atomic ISM mass strongly depends on the goodness-of-fit definition and is highly uncertain when the inferred spin temperature has a large uncertainty. Finally, we find that the relation of the H i surface density with the dust surface density and extinction has significantly more scatter when the inferred H i opacity correction is applied. These variations are difficult to explain without additionally requiring large variations in the dust properties. Based on these findings, we suggest that the opaque H i mass is best constrained by H i absorption studies. 

ABSTRACT Understanding the spatial distribution of metals within galaxies allows us to study the processes of chemical enrichment and mixing in the interstellar medium. In this work, we map the 2D distribution of metals using a Gaussian Process Regression (GPR) for 19 star-forming galaxies observed with the Very Large Telescope/Multi Unit Spectroscopic Explorer (VLT–MUSE) as a part of the PHANGS–MUSE survey. We find that 12 of our 19 galaxies show significant 2D metallicity variation. Those without significant variations typically have fewer metallicity measurements, indicating this is due to the dearth of ${\rm H\, {\small II}}$ regions in these galaxies, rather than a lack of higher-order variation. After subtracting a linear radial gradient, we see no enrichment in the spiral arms versus the disc. We measure the 50 per cent correlation scale from the two-point correlation function of these radially subtracted maps, finding it to typically be an order of magnitude smaller than the fitted GPR kernel scale length. We study the dependence of the two-point correlation scale length with a number of global galaxy properties. We find no relationship between the 50 per cent correlation scale and the overall gas turbulence, in tension with existing theoretical models. We also find more actively star-forming galaxies, and earlier type galaxies have a larger 50 per cent correlation scale. The size and stellar mass surface density do not appear to correlate with the 50 per cent correlation scale, indicating that perhaps the evolutionary state of the galaxy and its current star formation activity is the strongest indicator of the homogeneity of the metal distribution. 

ABSTRACT When completed, the PHANGS–HST project will provide a census of roughly 50 000 compact star clusters and associations, as well as human morphological classifications for roughly 20 000 of those objects. These large numbers motivated the development of a more objective and repeatable method to help perform source classifications. In this paper, we consider the results for five PHANGS–HST galaxies (NGC 628, NGC 1433, NGC 1566, NGC 3351, NGC 3627) using classifications from two convolutional neural network architectures (RESNET and VGG) trained using deep transfer learning techniques. The results are compared to classifications performed by humans. The primary result is that the neural network classifications are comparable in quality to the human classifications with typical agreement around 70 to 80 per cent for Class 1 clusters (symmetric, centrally concentrated) and 40 to 70 per cent for Class 2 clusters (asymmetric, centrally concentrated). If Class 1 and 2 are considered together the agreement is 82 ± 3 per cent. Dependencies on magnitudes, crowding, and background surface brightness are examined. A detailed description of the criteria and methodology used for the human classifications is included along with an examination of systematic differences between PHANGS–HST and LEGUS. The distribution of data points in a colour–colour diagram is used as a ‘figure of merit’ to further test the relative performances of the different methods. The effects on science results (e.g. determinations of mass and age functions) of using different cluster classification methods are examined and found to be minimal.