PHANGS-ML: Dissecting Multiphase Gas and Dust in Nearby Galaxies Using Machine Learning

The Astrophysical Journal American Astronomical Society 968:1 (2024) 24

Authors:

Dalya Baron, Karin M Sandstrom, Erik Rosolowsky, Oleg V Egorov, Ralf S Klessen, Adam K Leroy, Médéric Boquien, Eva Schinnerer, Francesco Belfiore, Brent Groves, Jérémy Chastenet, Daniel A Dale, Guillermo A Blanc, José E Méndez-Delgado, Eric W Koch, Kathryn Grasha, Mélanie Chevance, David A Thilker, Dario Colombo, Thomas G Williams, Debosmita Pathak, Jessica Sutter, Toby Brown, John F Wu

Abstract:

The PHANGS survey uses Atacama Large Millimeter/submillimeter Array, Hubble Space Telescope, Very Large Telescope, and JWST to obtain an unprecedented high-resolution view of nearby galaxies, covering millions of spatially independent regions. The high dimensionality of such a diverse multiwavelength data set makes it challenging to identify new trends, particularly when they connect observables from different wavelengths. Here, we use unsupervised machine-learning algorithms to mine this information-rich data set to identify novel patterns. We focus on three of the PHANGS-JWST galaxies, for which we extract properties pertaining to their stellar populations; warm ionized and cold molecular gas; and polycyclic aromatic hydrocarbons (PAHs), as measured over 150 pc scale regions. We show that we can divide the regions into groups with distinct multiphase gas and PAH properties. In the process, we identify previously unknown galaxy-wide correlations between PAH band and optical line ratios and use our identified groups to interpret them. The correlations we measure can be naturally explained in a scenario where the PAHs and the ionized gas are exposed to different parts of the same radiation field that varies spatially across the galaxies. This scenario has several implications for nearby galaxies: (i) The uniform PAH ionized fraction on 150 pc scales suggests significant self-regulation in the interstellar medium, (ii) the PAH 11.3/7.7 μm band ratio may be used to constrain the shape of the non-ionizing far-ultraviolet to optical part of the radiation field, and (iii) the varying radiation field affects line ratios that are commonly used as PAH size diagnostics. Neglecting this effect leads to incorrect or biased PAH sizes.

SYREN-HALOFIT: A fast, interpretable, high-precision formula for the ΛCDM nonlinear matter power spectrum

Astronomy & Astrophysics EDP Sciences 686 (2024) ARTN A150

Authors:

Deaglan J Bartlett, Benjamin D Wandelt, Matteo Zennaro, Pedro G Ferreira, Harry Desmond

Abstract:

<jats:p><jats:italic>Context.</jats:italic>Rapid and accurate evaluation of the nonlinear matter power spectrum,<jats:italic>P</jats:italic>(<jats:italic>k</jats:italic>), as a function of cosmological parameters and redshift is of fundamental importance in cosmology. Analytic approximations provide an interpretable solution, yet current approximations are neither fast nor accurate relative to numerical emulators.</jats:p><jats:p><jats:italic>Aims.</jats:italic>We aim to accelerate symbolic approximations to<jats:italic>P</jats:italic>(<jats:italic>k</jats:italic>) by removing the requirement to perform integrals, instead using short symbolic expressions to compute all variables of interest. We also wish to make such expressions more accurate by re-optimising the parameters of these models (using a larger number of cosmologies and focussing on cosmological parameters of more interest for present-day studies) and providing correction terms.</jats:p><jats:p><jats:italic>Methods.</jats:italic>We use symbolic regression to obtain simple analytic approximations to the nonlinear scale,<jats:italic>k</jats:italic><jats:sub><jats:italic>σ</jats:italic></jats:sub>, the effective spectral index,<jats:italic>n</jats:italic><jats:sub>eff</jats:sub>, and the curvature,<jats:italic>C</jats:italic>, which are required for the<jats:sc>HALOFIT</jats:sc>model. We then re-optimise the coefficients of<jats:sc>HALOFIT</jats:sc>to fit a wide range of cosmologies and redshifts. We then again exploit symbolic regression to explore the space of analytic expressions to fit the residuals between<jats:italic>P</jats:italic>(<jats:italic>k</jats:italic>) and the optimised predictions of<jats:sc>HALOFIT</jats:sc>. Our results are designed to match the predictions of<jats:sc>EUCLIDEMULATOR</jats:sc>2, but we validate our methods against<jats:italic>N</jats:italic>-body simulations.</jats:p><jats:p><jats:italic>Results.</jats:italic>We find symbolic expressions for<jats:italic>k</jats:italic><jats:sub><jats:italic>σ</jats:italic></jats:sub>,<jats:italic>n</jats:italic><jats:sub>eff</jats:sub>and<jats:italic>C</jats:italic>which have root mean squared fractional errors of 0.8%, 0.2% and 0.3%, respectively, for redshifts below 3 and a wide range of cosmologies. We provide re-optimised<jats:sc>HALOFIT</jats:sc>parameters, which reduce the root mean squared fractional error (compared to<jats:sc>EUCLIDEMULATOR</jats:sc>2) from 3% to below 2% for wavenumbers<jats:italic>k</jats:italic> = 9 × 10<jats:sup>−3</jats:sup> − 9 <jats:italic>h</jats:italic> Mpc<jats:sup>−1</jats:sup>. We introduce<jats:sc>SYREN-HALOFIT</jats:sc>(symbolic-regression-enhanced<jats:sc>HALOFIT</jats:sc>), an extension to<jats:sc>HALOFIT</jats:sc>containing a short symbolic correction which improves this error to 1%. Our method is 2350 and 3170 times faster than current<jats:sc>HALOFIT</jats:sc>and<jats:sc>HMCODE</jats:sc>implementations, respectively, and 2680 and 64 times faster than<jats:sc>EUCLIDEMULATOR</jats:sc>2 (which requires running<jats:sc>CLASS</jats:sc>) and the<jats:sc>BACCO</jats:sc>emulator. We obtain comparable accuracy to<jats:sc>EUCLIDEMULATOR</jats:sc>2 and the<jats:sc>BACCO</jats:sc>emulator when tested on<jats:italic>N</jats:italic>-body simulations.</jats:p><jats:p><jats:italic>Conclusions.</jats:italic>Our work greatly increases the speed and accuracy of symbolic approximations to<jats:italic>P</jats:italic>(<jats:italic>k</jats:italic>), making them significantly faster than their numerical counterparts without loss of accuracy.</jats:p>

WISDOM Project – XXI. Giant molecular clouds in the central region of the barred spiral galaxy NGC 613: a steep size – linewidth relation

Monthly Notices of the Royal Astronomical Society Oxford University Press (OUP) (2024) stae1394-stae1394

Authors:

Woorak Choi, Martin Bureau, Lijie Liu, Michele Cappellari, Timothy A Davis, Jindra Gensior, Fu-Heng Liang, Anan Lu, Sanghyuk Moon, Ilaria Ruffa, Thomas G Williams, Aeree Chung

Abstract:

<jats:title>Abstract</jats:title> <jats:p>NGC 613 is a nearby barred spiral galaxy with a nuclear ring. Exploiting high spatial resolution (≈20 pc) Atacama Large Millimeter/sub-millimeter Array 12CO(1-0) observations, we study the giant molecular clouds (GMCs) in the nuclear ring and its vicinity, identifying 158 spatially- and spectrally-resolved GMCs. The GMC sizes (Rc) are comparable to those of the clouds in the Milky Way (MW) disc, but their gas masses, observed linewidths (σobs, los) and gas mass surface densities are larger. The GMC size – linewidth relation ($\sigma _{\mathrm{obs,los}}\propto R_{\mathrm{c}}^{0.77}$) is steeper than that of the clouds of the MW disc and centre, and the GMCs are on average only marginally gravitationally bound (with a mean virial parameter 〈αobs, vir〉 ≈ 1.7). We discuss the possible origins of the steep size – linewidth relation and enhanced observed linewidths of the clouds and suggest that a combination of mechanisms such as stellar feedback, gas accretion and cloud-cloud collisions, as well as the gas inflows driven by the large-scale bar, may play a role.</jats:p>

WISDOM Project – XXI. Giant molecular clouds in the central region of the barred spiral galaxy NGC 613: a steep size – linewidth relation

Monthly Notices of the Royal Astronomical Society Oxford University Press 531:4 (2024) 4045-4059

Authors:

Woorak Choi, Martin Bureau, Lijie Liu, Michele Cappellari, Timothy A Davis, Jindra Gensior, Fu-Heng Liang, Anan Lu, Sanghyuk Moon, Ilaria Ruffa, Thomas G Williams, Aeree Chung

Abstract:

NGC 613 is a nearby barred spiral galaxy with a nuclear ring. Exploiting high spatial resolution (≈20 pc) Atacama Large Millimeter/sub-millimeter Array 12CO(1-0) observations, we study the giant molecular clouds (GMCs) in the nuclear ring and its vicinity, identifying 158 spatially- and spectrally-resolved GMCs. The GMC sizes (Rc) are comparable to those of the clouds in the Milky Way (MW) disc, but their gas masses, observed linewidths (σobs, los) and gas mass surface densities are larger. The GMC size – linewidth relation ($\sigma _{\mathrm{obs,los}}\propto R_{\mathrm{c}}^{0.77}$) is steeper than that of the clouds of the MW disc and centre, and the GMCs are on average only marginally gravitationally bound (with a mean virial parameter 〈αobs, vir〉 ≈ 1.7). We discuss the possible origins of the steep size – linewidth relation and enhanced observed linewidths of the clouds and suggest that a combination of mechanisms such as stellar feedback, gas accretion and cloud-cloud collisions, as well as the gas inflows driven by the large-scale bar, may play a role.

WISDOM project XX. – Strong shear tearing molecular clouds apart in NGC 524

Monthly Notices of the Royal Astronomical Society Oxford University Press (OUP) 531:4 (2024) stae1395-stae1395

Authors:

Anan Lu, Daryl Haggard, Martin Bureau, Jindra Gensior, Sarah Jeffreson, Carmelle Robert, Thomas G Williams, Fu-Heng Liang, Woorak Choi, Timothy A Davis, Sara Babic, Hope Boyce, Benjamin Cheung, Laurent Drissen, Jacob S Elford, Lijie Liu, Thomas Martin, Carter Rhea, Laurie Rousseau-Nepton, Ilaria Ruffa

Abstract:

<jats:title>ABSTRACT</jats:title> <jats:p>Early-type galaxies (ETGs) are known to harbour dense spheroids of stars but scarce star formation (SF). Approximately a quarter of these galaxies have rich molecular gas reservoirs yet do not form stars efficiently. We study here the ETG NGC 524, with strong shear suspected to result in a smooth molecular gas disc and low star-formation efficiency (SFE). We present new spatially resolved observations of the 12CO(2-1)-emitting cold molecular gas from the Atacama Large Millimeter/sub-millimeter Array (ALMA) and of the warm ionized-gas emission lines from SITELLE at the Canada–France–Hawaii Telescope. Although constrained by the resolution of the ALMA observations (≈37 pc), we identify only 52 GMCs with radii ranging from 30 to 140 pc, a low mean molecular gas mass surface density 〈Σgas〉 ≈ 125 M⊙ pc−2 and a high mean virial parameter 〈αobs, vir〉 ≈ 5.3. We measure spatially resolved molecular gas depletion times (τdep ≡ 1/SFE) with a spatial resolution of ≈100 pc within a galactocentric distance of 1.5 kpc. The global depletion time is ≈2.0 Gyr but τdep increases towards the galaxy centre, with a maximum τdep, max ≈ 5.2 Gyr. However, no pure H ii region is identified in NGC 524 using ionized-gas emission-line ratio diagnostics, so the τdep inferred are in fact lower limits. Measuring the GMC properties and dynamical states, we conclude that shear is the dominant mechanism shaping the molecular gas properties and regulating SF in NGC 524. This is supported by analogous analyses of the GMCs in a simulated ETG similar to NGC 524.</jats:p>