Showing 1–50 of 120 results for author: Schruba, A

A sensitive, high-resolution, wide-field IRAM NOEMA CO(1-0) survey of the very nearby spiral galaxy IC 342

Authors: M. Querejeta, J. Pety, A. Schruba, A. K. Leroy, C. N. Herrera, I-D. Chiang, S. E. Meidt, E. Rosolowsky, E. Schinnerer, K. Schuster, J. Sun, K. A. Herrmann, A. T. Barnes, I. Beslic, F. Bigiel, Y. Cao, M. Chevance, C. Eibensteiner, E. Emsellem, C. M. Faesi, A. Hughes, J. Kim, R. S. Klessen, K. Kreckel, J. M. D. Kruijssen , et al. (9 additional authors not shown)

Abstract: We present a new wide-field 10.75 x 10.75 arcmin^2 (~11x11 kpc^2), high-resolution (theta = 3.6" ~ 60 pc) NOEMA CO(1-0) survey of the very nearby (d=3.45 Mpc) spiral galaxy IC 342. The survey spans out to about 1.5 effective radii and covers most of the region where molecular gas dominates the cold interstellar medium. We resolved the CO emission into >600 individual giant molecular clouds and ass… ▽ More We present a new wide-field 10.75 x 10.75 arcmin^2 (~11x11 kpc^2), high-resolution (theta = 3.6" ~ 60 pc) NOEMA CO(1-0) survey of the very nearby (d=3.45 Mpc) spiral galaxy IC 342. The survey spans out to about 1.5 effective radii and covers most of the region where molecular gas dominates the cold interstellar medium. We resolved the CO emission into >600 individual giant molecular clouds and associations. We assessed their properties and found that overall the clouds show approximate virial balance, with typical virial parameters of alpha_vir=1-2. The typical surface density and line width of molecular gas increase from the inter-arm region to the arm and bar region, and they reach their highest values in the inner kiloparsec of the galaxy (median Sigma_mol~80, 140, 160, and 1100 M_sun/pc^2, sigma_CO~6.6, 7.6, 9.7, and 18.4 km/s for inter-arm, arm, bar, and center clouds, respectively). Clouds in the central part of the galaxy show an enhanced line width relative to their surface densities and evidence of additional sources of dynamical broadening. All of these results agree well with studies of clouds in more distant galaxies at a similar physical resolution. Leveraging our measurements to estimate the density and gravitational free-fall time at 90 pc resolution, averaged on 1.5 kpc hexagonal apertures, we estimate a typical star formation efficiency per free-fall time of 0.45% with a 16-84% variation of 0.33-0.71% among such 1.5 kpc regions. We speculate that bar-driven gas inflow could explain the large gas concentration in the central kiloparsec and the buildup of the massive nuclear star cluster. This wide-area CO map of the closest face-on massive spiral galaxy demonstrates the current mapping power of NOEMA and has many potential applications. The data and products are publicly available. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 16 pages, 6 figures, accepted for publication in A&A

Journal ref: A&A 680, A4 (2023)

The ALMOND Survey: Molecular cloud properties and gas density tracers across 25 nearby spiral galaxies with ALMA

Authors: Lukas Neumann, Molly J. Gallagher, Frank Bigiel, Adam K. Leroy, Ashley T. Barnes, Antonio Usero, Jakob S. den Brok, Francesco Belfiore, Ivana Bešlić, Yixian Cao, Mélanie Chevance, Daniel A. Dale, Cosima Eibensteiner, Simon C. O. Glover, Kathryn Grasha, Jonathan D. Henshaw, María J. Jiménez-Donaire, Ralf S. Klessen, J. M. Diederik Kruijssen, Daizhong Liu, Sharon Meidt, Jérôme Pety, Johannes Puschnig, Miguel Querejeta, Erik Rosolowsky , et al. (6 additional authors not shown)

Abstract: We use new HCN(1-0) data from the ALMOND (ACA Large-sample Mapping Of Nearby galaxies in Dense gas) survey to trace the kpc-scale molecular gas density structure and CO(2-1) data from 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 SFR/HCN ratio to trace the star formation efficiency i… ▽ More We use new HCN(1-0) data from the ALMOND (ACA Large-sample Mapping Of Nearby galaxies in Dense gas) survey to trace the kpc-scale molecular gas density structure and CO(2-1) data from 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 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. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: 48 pages, 40 figures

Multi-Scale Stellar Associations across the Star Formation Hierarchy in PHANGS-HST Nearby Galaxies: Methodology and Properties

Authors: Kirsten L. Larson, Janice C. Lee, David A. Thilker, Bradley C. Whitmore, Sinan Deger, James Lilly, Rupali Chandar, Daniel A. Dale, Frank Bigiel, Kathryn Grasha, Brent Groves, Ralf S. Klessen, Kathryn Kreckel, J. M. Diederik Kruijssen, Adam K. Leroy, Hsi-An Pan, Erik Rosolowsky, Eva Schinnerer, Andreas Schruba, Elizabeth J. Watkins, Thomas G. Williams

Abstract: We develop a method to identify and determine the physical properties of stellar associations using Hubble Space Telescope (HST) NUV-U-B-V-I imaging of nearby galaxies from the PHANGS-HST survey. We apply a watershed algorithm to density maps constructed from point source catalogues Gaussian smoothed to multiple physical scales from 8 to 64 pc. We develop our method on two galaxies that span the d… ▽ More We develop a method to identify and determine the physical properties of stellar associations using Hubble Space Telescope (HST) NUV-U-B-V-I imaging of nearby galaxies from the PHANGS-HST survey. We apply a watershed algorithm to density maps constructed from point source catalogues Gaussian smoothed to multiple physical scales from 8 to 64 pc. We develop our method on two galaxies that span the distance range in the PHANGS-HST sample: NGC 3351 (10 Mpc), NGC 1566 (18 Mpc). We test our algorithm with different parameters such as the choice of detection band for the point source catalogue (NUV or V), source density image filtering methods, and absolute magnitude limits. We characterise the properties of the resulting multi-scale associations, including sizes, number of tracer stars, number of associations, photometry, as well as ages, masses, and reddening from Spectral Energy Distribution fitting. Our method successfully identifies structures that occupy loci in the UBVI colour-colour diagram consistent with previously published catalogues of clusters and associations. The median ages of the associations increases from log(age/yr) = 6.6 to log(age/yr) = 6.9 as the spatial scale increases from 8 pc to 64 pc for both galaxies. We find that the youngest stellar associations, with ages < 3 Myr, indeed closely trace H ii regions in H$α$ imaging, and that older associations are increasingly anti-correlated with the H$α$ emission. Owing to our new method, the PHANGS-HST multi-scale associations provide a far more complete census of recent star formation activity than found with previous cluster and compact association catalogues. The method presented here will be applied to the full sample of 38 PHANGS-HST galaxies. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Submitted to MNRAS. Referee report received with minor comments, and "request to clarify if the smaller associations are always included in the larger ones and how this may affect the photometric fitting of the larger association if the groups have different ages." Revision in progress

CI and CO in Nearby Spiral Galaxies -- I. Line Ratio and Abundance Variations at ~ 200 pc Scales

Authors: Daizhong Liu, Eva Schinnerer, Toshiki Saito, Erik Rosolowsky, Adam Leroy, Antonio Usero, Karin Sandstrom, Ralf S. Klessen, Simon C. O. Glover, Yiping Ao, Ivana Bešlić, Frank Bigiel, Yixian Cao, Jérémy Chastenet, Mélanie Chevance, Daniel A. Dale, Yu Gao, Annie Hughes, Kathryn Kreckel, J. M. Diederik Kruijssen, Hsi-An Pan, Jérôme Pety, Dragan Salak, Francesco Santoro, Andreas Schruba , et al. (3 additional authors not shown)

Abstract: We present new neutral atomic carbon [CI](3P1-3P0) mapping observations within the inner ~7 kpc and ~4 kpc of the disks of NGC3627 and NGC4321 at a spatial resolution of 190 pc and 270 pc, respectively, using the ALMA Atacama Compact Array (ACA). We combine these with the CO(2-1) data from PHANGS-ALMA, and literature [CI] and CO data for two other starburst and/or active galactic nucleus (AGN) gal… ▽ More We present new neutral atomic carbon [CI](3P1-3P0) mapping observations within the inner ~7 kpc and ~4 kpc of the disks of NGC3627 and NGC4321 at a spatial resolution of 190 pc and 270 pc, respectively, using the ALMA Atacama Compact Array (ACA). We combine these with the CO(2-1) data from PHANGS-ALMA, and literature [CI] and CO data for two other starburst and/or active galactic nucleus (AGN) galaxies (NGC1808, NGC7469), to study: a) the spatial distributions of CI and CO emission; b) the observed line ratio RCICO = I_[CI](1-0)/I_CO(2-1) as a function of various galactic properties; and c) the abundance ratio of [CI/CO]. We find excellent spatial correspondence between CI and CO emission and nearly uniform RCICO ~0.1 across the majority of the star-forming disks of NGC3627 and NGC4321. However, RCICO strongly varies from ~0.05 at the centre of NGC4321 to >0.2-0.5 in NGC1808's starburst centre and NGC7469's centre with an X-ray AGN. Meanwhile, RCICO does not obviously vary with $U$, similar to the prediction of PDR models. We also find a mildly decreasing RCICO with an increasing metallicity over 0.7-0.85 solar metallicity, consistent with the literature. Assuming various typical ISM conditions representing GMCs, active star-forming regions and strong starbursting environments, we calculate the LTE radiative transfer and estimate the [CI/CO] abundance ratio to be ~0.1 across the disks of NGC3627 and NGC4321, similar to previous large-scale findings in Galactic studies. However, this abundance ratio likely has a substantial increase to ~1 and >1-5 in NGC1808's starburst and NGC7469's strong AGN environments, respectively, in line with the expectations for cosmic-ray dominated region (CRDR) and X-ray dominated region (XDR) chemistry. Finally, we do not find a robust evidence for a generally CO-dark, CI-bright gas in the disk areas we probed. (abbreviated) △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: 23 pages, 13 figures and one table in total (17 pages and 9 figures in main text). Accepted for publication in A&A. For associated data cubes and moment maps, see https://www.canfar.net/storage/vault/list/phangs/RELEASES/DZLIU_etal_2022

Journal ref: A&A 672, A36 (2023)

PDFchem: A new fast method to determine ISM properties and infer environmental parameters using probability distributions

Authors: Thomas G. Bisbas, Ewine F. van Dishoeck, Chia-Yu Hu, Andreas Schruba

Abstract: Determining the atomic and molecular content of the interstellar medium (ISM) as a function of environmental parameters is of fundamental importance to understand the star-formation process across the epochs. Although there exist various three-dimensional hydro-chemical codes modelling the ISM at different scales and redshifts, they are computationally expensive and inefficient for studies over a… ▽ More Determining the atomic and molecular content of the interstellar medium (ISM) as a function of environmental parameters is of fundamental importance to understand the star-formation process across the epochs. Although there exist various three-dimensional hydro-chemical codes modelling the ISM at different scales and redshifts, they are computationally expensive and inefficient for studies over a large parameter space. Building on our earlier approach, we present PDFchem, a novel algorithm that models the cold ISM at moderate and large scales using functions connecting the quantities of the local ($A_{\rm V,eff}$) and the observed ($A_{\rm V,obs}$) visual extinctions, and the local number density, $n_{\rm H}$, with probability density functions (PDF) of $A_{\rm V,obs}$ on cloud scales typically tens-to-hundreds of pc as an input. For any given $A_{\rm V,obs}$-PDF, sampled with thousands of clouds, the algorithm instantly computes the average abundances of the most important species (HI, H$_2$, CII, CI, CO, OH, OH$^+$, H$_2$O$^+$, CH, HCO$^+$) and performs radiative transfer calculations to estimate the average emission of the most commonly observed lines ([CII]~$158μ$m, both [CI] fine-structure lines and the first five rotational transitions of $^{12}$CO). We examine two $A_{\rm V,obs}$-PDFs corresponding to a non star-forming and a star-forming ISM region, under a variety of environmental parameters combinations. These cover FUV intensities in the range of $χ/χ_0=10^{-1}-10^3$, cosmic-ray ionization rates in the range of $ζ_{\rm CR}=10^{-17}-10^{-13}\,{\rm s}^{-1}$ and metallicities in the range of $Z=0.1-2\,{\rm Z}_{\odot}$. PDFchem is fast, easy to use, reproduces the PDR quantities of the time-consuming hydrodynamical models and can be used directly with observed data to understand the evolution of the cold ISM chemistry. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 33 pages, 27 figures. Accepted in MNRAS. The algorithm can be found in: https://github.com/tbisbas/PDFchem. Comments welcome!

PHANGS: Constraining Star Formation Timescales Using the Spatial Correlations of Star Clusters and Giant Molecular Clouds

Authors: Jordan A. Turner, Daniel A. Dale, James Lilly, Mederic Boquien, Sinan Deger, Janice C. Lee, Bradley C. Whitmore, Gagandeep S. Anand, Samantha M. Benincasa, Frank Bigiel, Guillermo A. Blanc, Melanie Chevance, Eric Emsellem, Christopher M. Faesi, Simon C. O. Glover, Kathryn Grasha, Annie Hughes, Ralf S. Klessen, Kathryn Kreckel, J. M. Diederik Kruijssen, Adam K. Leroy, Hsi-An Pan, Erik Rosolowsky, Andreas Schruba, Thomas G. Williams

Abstract: In the hierarchical view of star formation, giant molecular gas 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.… ▽ More In the hierarchical view of star formation, giant molecular gas 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 ($\leq$ 10 Myr) star clusters are more strongly autocorrelated on kpc and smaller spatial scales than the >10 Myr stellar populations, indicating that the hierarchical structure dissolves over time. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: 15 pages, 11 figures, 4 tables. Accepted to MNRAS Sept 6 2022

Environmental dependence of the molecular cloud lifecycle in 54 main sequence galaxies

Authors: Jaeyeon Kim, Mélanie Chevance, J. M. Diederik Kruijssen, Adam K. Leroy, Andreas Schruba, Ashley T. Barnes, Frank Bigiel, Guillermo A. Blanc, Yixian Cao, Enrico Congiu, Daniel A. Dale, Christopher M. Faesi, Simon C. O. Glover, Kathryn Grasha, Brent Groves, Annie Hughes, Ralf S. Klessen, Kathryn Kreckel, Rebecca McElroy, Hsi-An Pan, Jérôme Pety, Miguel Querejeta, Alessandro Razza, Erik Rosolowsky, Toshiki Saito , et al. (5 additional authors not shown)

Abstract: The processes of star formation and feedback, regulating the cycle of matter between gas and stars on the scales of giant molecular clouds (GMCs; $\sim$100pc), play a major role in governing galaxy evolution. Measuring the time-scales of GMC evolution is important to identify and characterise the specific physical mechanisms that drive this transition. By applying a robust statistical method to hi… ▽ More The processes of star formation and feedback, regulating the cycle of matter between gas and stars on the scales of giant molecular clouds (GMCs; $\sim$100pc), play a major role in governing galaxy evolution. Measuring the time-scales of GMC evolution is important to identify and characterise the specific physical mechanisms that drive this transition. By applying a robust statistical method to high-resolution CO and narrow-band H$α$ imaging from the PHANGS survey, we systematically measure the evolutionary timeline from molecular clouds to exposed young stellar regions on GMC scales, across the discs of an unprecedented sample of 54 star-forming main-sequence galaxies (excluding their unresolved centres). We find that clouds live for about $1{-}3$ GMC turbulence crossing times ($5{-}30$Myr) and are efficiently dispersed by stellar feedback within $1{-}5$Myr once the star-forming region becomes partially exposed, resulting in integrated star formation efficiencies of $1{-}8$%. These ranges reflect physical galaxy-to-galaxy variation. In order to evaluate whether galactic environment influences GMC evolution, we correlate our measurements with average properties of the GMCs and their local galactic environment. We find several strong correlations that can be physically understood, revealing a quantitative link between galactic-scale environmental properties and the small-scale GMC evolution. Notably, the measured CO-visible cloud lifetimes become shorter with decreasing galaxy mass, mostly due to the increasing presence of CO-dark molecular gas in such environment. Our results represent a first step towards a comprehensive picture of cloud assembly and dispersal, which requires further extension and refinement with tracers of the atomic gas, dust, and deeply-embedded stars. △ Less

Submitted 22 August, 2022; v1 submitted 20 June, 2022; originally announced June 2022.

Comments: 24 pages, 8 figures, 5 tables; accepted for publication in MNRAS (August 16, 2022)

Molecular Cloud Populations in the Context of Their Host Galaxy Environments: A Multiwavelength Perspective

Authors: Jiayi Sun, Adam K. Leroy, Erik Rosolowsky, Annie Hughes, Eva Schinnerer, Andreas Schruba, Eric W. Koch, Guillermo A. Blanc, I-Da Chiang, Brent Groves, Daizhong Liu, Sharon Meidt, Hsi-An Pan, Jerome Pety, Miguel Querejeta, Toshiki Saito, Karin Sandstrom, Amy Sardone, Antonio Usero, Dyas Utomo, Thomas G. Williams, Ashley T. Barnes, Samantha M. Benincasa, Frank Bigiel, Alberto D. Bolatto , et al. (13 additional authors not shown)

Abstract: We present a rich, multiwavelength, multiscale database built around the PHANGS-ALMA CO$\,$(2-1) survey and ancillary data. We use this database to present the distributions of molecular cloud populations and sub-galactic environments in 80 PHANGS galaxies, to characterize the relationship between population-averaged cloud properties and host galaxy properties, and to assess key timescales relevan… ▽ More We present a rich, multiwavelength, multiscale database built around the PHANGS-ALMA CO$\,$(2-1) survey and ancillary data. We use this database to present the distributions of molecular cloud populations and sub-galactic environments in 80 PHANGS galaxies, to characterize the relationship between population-averaged cloud properties and host galaxy properties, and to assess key timescales relevant to molecular cloud evolution and star formation. We show that PHANGS probes a wide range of kpc-scale gas, stellar, and star formation rate (SFR) surface densities, as well as orbital velocities and shear. The population-averaged cloud properties in each aperture correlate strongly with both local environmental properties and host galaxy global properties. Leveraging a variable selection analysis, we find that the kpc-scale surface densities of molecular gas and SFR tend to possess the most predictive power for the population-averaged cloud properties. Once their variations are controlled for, galaxy global properties contain little additional information, which implies that the apparent galaxy-to-galaxy variations in cloud populations are likely mediated by kpc-scale environmental conditions. We further estimate a suite of important timescales from our multiwavelength measurements. The cloud-scale free-fall time and turbulence crossing time are ${\sim}5{-}20$ Myr, comparable to previous cloud lifetime estimates. The timescales for orbital motion, shearing, and cloud-cloud collisions are longer, ${\sim}100$ Myr. The molecular gas depletion time is $1{-}3$ Gyr and shows weak to no correlations with the other timescales in our data. We publish our measurements online and expect them to have broad utility to future studies of molecular clouds and star formation. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: 32 pages + 6 appendices. AJ in press. Data products available at https://www.canfar.net/storage/list/phangs/RELEASES/Sun_etal_2022 . Associated software package available at https://doi.org/10.5281/zenodo.6584842 . Slides summarizing the key results can be found at https://www.dropbox.com/s/7llnrob9wi0isfk/Sun_et_PHANGS_2022.pptx?dl=0

Journal ref: AJ, 164, 43 (2022)

Variations in the $Σ_{\rm SFR} {-} Σ_{\rm mol} {-} Σ_{\rm \star}$ plane across galactic environments in PHANGS galaxies

Authors: I. Pessa, E. Schinnerer, A. Leroy, E. Koch, E. Rosolowsky, T. Williams, H. -A. Pan, A. Schruba, A. Usero, F. Belfiore, F. Bigiel, G. Blanc, M. Chevance, D. Dale, E. Emsellem, J. Gensior, S. Glover, K. Grasha, B. Groves, R. Klessen, K. Kreckel, J. M. D. Kruijssen, D. Liu, S. E. Meidt, J. Pety , et al. (4 additional authors not shown)

Abstract: There exists some consensus that stellar mass surface density ($Σ_{*}$) and molecular gas mass surface density ($Σ_{\rm mol}$) are the main quantities responsible for locally setting the star formation rate. This regulation is inferred from locally resolved scaling relations between these two quantities and the star formation rate surface density ($Σ_{\rm SFR}$). However, the universality of these… ▽ More There exists some consensus that stellar mass surface density ($Σ_{*}$) and molecular gas mass surface density ($Σ_{\rm mol}$) are the main quantities responsible for locally setting the star formation rate. This regulation is inferred from locally resolved scaling relations between these two quantities and the star formation rate surface density ($Σ_{\rm SFR}$). However, the universality of these relations is debated. Here, we probe the interplay between these three quantities across different galactic environments at a spatial resolution of 150 pc. We perform a hierarchical Bayesian linear regression to find the best set of parameters $C_{*}$, $C_{\rm mol}$, and $C_{\rm norm}$ that describe the star-forming plane conformed by these quantities, such that $\log Σ_{\rm SFR} = C_{*} \log Σ_{*} + C_{\rm mol} \log Σ_{\rm mol} + C_{\rm norm}$, and explore variations in the determined parameters across galactic environments, focusing our analysis on the $C_{*}$ and $C_{\rm mol}$ slopes. We find signs of variations in the posterior distributions of $C_{*}$ and $C_{\rm mol}$ across different galactic environments. Bars show the most negative value of $C_{*}$, a sign of longer depletion times, while spiral arms show the highest $C_{*}$ among all environments. We conclude that systematic variations in the interplay of $Σ_{*}$, $Σ_{\rm mol}$ and $Σ_{\rm SFR}$ across galactic environments exist at a spatial resolution of 150 pc, and we interpret these variations as produced by an additional mechanism regulating the formation of stars that is not captured by either $Σ_{*}$ or $Σ_{\rm mol}$. We find that these variations correlate with changes in the star formation efficiency across environments, which could be linked to the dynamical state of the gas that prevents it from collapsing and forming stars, or to changes in the molecular gas fraction. △ Less

Submitted 16 June, 2022; v1 submitted 22 March, 2022; originally announced March 2022.

Comments: 21 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 663, A61 (2022)

WISDOM Project -- X. The morphology of the molecular ISM in galaxy centres and its dependence on galaxy structure

Authors: Timothy A. Davis, Jindra Gensior, Martin Bureau, Michele Cappellari, Woorak Choi, Jacob S. Elford, J. M. Diederik Kruijssen, Federico Lelli, Fu-Heng Liang, Lijie Liu, Ilaria Ruffa, Toshiki Saito, Marc Sarzi, Andreas Schruba, Thomas G. Williams

Abstract: We use high-resolution maps of the molecular interstellar medium (ISM) in the centres of eighty-six 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 galaxi… ▽ More We use high-resolution maps of the molecular interstellar medium (ISM) in the centres of eighty-six 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 idealised 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. △ Less

Submitted 7 March, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: 20 pages, 9 figures. Accepted for publication in MNRAS

Stellar feedback in M83 as observed with MUSE -- I. Overview, an unprecedented view of the stellar and gas kinematics and evidence of outflowing gas

Authors: Lorenza Della Bruna, Angela Adamo, Philippe Amram, Erik Rosolowsky, Christopher Usher, Mattia Sirressi, Andreas Schruba, Eric Emsellem, Adam Leroy, Arjan Bik, William P. Blair, Anna F. McLeod, Göran Östlin, Florent Renaud, Carmelle Robert, Laurie Rousseau-Nepton, Linda J. Smith

Abstract: We present a large VLT/MUSE mosaic (3.8 x 3.8 kpc) of the nearby spiral galaxy M83, with a spatial resolution ~20 pc. We obtained the kinematics of the stars and ionised gas, and compared them with molecular gas kinematics from ALMA CO(2-1). We separated the ionised gas into HII regions and diffuse ionised gas (DIG) and determined the fraction of Ha luminosity originating from the DIG (f_DIG). We… ▽ More We present a large VLT/MUSE mosaic (3.8 x 3.8 kpc) of the nearby spiral galaxy M83, with a spatial resolution ~20 pc. We obtained the kinematics of the stars and ionised gas, and compared them with molecular gas kinematics from ALMA CO(2-1). We separated the ionised gas into HII regions and diffuse ionised gas (DIG) and determined the fraction of Ha luminosity originating from the DIG (f_DIG). We observe that both stars and gas trace the galactic disk rotation, as well as a fast-rotating nuclear component, likely connected to secular processes driven by the galactic bar. In the gas kinematics, we observe a stream east of the nucleus, redshifted with respect to the disk. The stream is surrounded by an extended ionised gas region with enhanced velocity dispersion and a high ionisation state, which is largely consistent with being ionised by slow shocks. We interpret this feature as either the superposition of the disk and an extraplanar layer of DIG, or as a bar-driven inflow of shocked gas. A double Gaussian component fit to the Ha line also reveals the presence of a nuclear biconic structure whose axis of symmetry is perpendicular to the bar. The two cones appear blue- and redshifted along the line of sight and stand out for having an Ha emission separated by up to 200 km s-1 from that of the disk, and a high velocity dispersion ~80-200 km s-1. At the far end of the cones, we observe that the gas is consistent with being ionised by shocks. These features had never been observed before in M83; we postulate that they are tracing a starburst-driven outflow shocking into the surrounding ISM. Finally, we obtain f_DIG ~ 13% in our field of view. We inspect the emission of the HII regions and DIG in `BPT' diagrams, finding that in HII regions photoionisation accounts for 99.8% of the Ha flux, whereas the DIG has a mixed contribution from photoionisation (94.9%) and shocks (5.1%). [abridged] △ Less

Submitted 3 February, 2022; originally announced February 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 660, A77 (2022)

A CO isotopologue Line Atlas within the Whirlpool galaxy Survey (CLAWS)

Authors: Jakob S. den Brok, Frank Bigiel, Kazimierz Sliwa, Toshiki Saito, Antonio Usero, Eva Schinnerer, Adam K. Leroy, María J. Jiménez-Donaire, Erik Rosolowsky, Ashley T. Barnes, Johannes Puschnig, Jérôme Pety, Andreas Schruba, Ivana Bešlić, Yixian Cao, Cosima Eibensteiner, Simon C. O. Glover, Ralf S. Klessen, J. M. Diederik Kruijssen, Sharon E. Meidt, Lukas Neumann, Neven Tomičić, Hsi-An Pan, Miguel Querejeta, Elizabeth Watkins , et al. (2 additional authors not shown)

Abstract: We present the CO isotopologue Line Atlas within the Whirpool galaxy Survey (CLAWS) based on an IRAM 30-m large programme which provides a benchmark study of numerous, faint CO isotopologues in the mm-wavelength regime across the full disc of M51 (NGC 5194). The survey's core goal is to use the low-J CO isotopologue lines to constrain CO excitation and chemistry, and therefrom the local physical c… ▽ More We present the CO isotopologue Line Atlas within the Whirpool galaxy Survey (CLAWS) based on an IRAM 30-m large programme which provides a benchmark study of numerous, faint CO isotopologues in the mm-wavelength regime across the full disc of M51 (NGC 5194). The survey's core goal is to use the low-J CO isotopologue lines to constrain CO excitation and chemistry, and therefrom the local physical conditions of the gas. In this survey paper, we describe the CLAWS observing and data reduction strategies. We map the J=1-0 and 2-1 transitions of the CO isotopologues $^{12}$CO,$^{13}$CO, C$^{18}$O and C$^{17}$O, as well as several supplementary lines within the 1 mm and 3 mm window (CN(1-0), CS(2-1), CH$_3$OH(2-1), N$_2$H$^+$(1-0), HC$_3$N(10-9)) at ~1 kpc resolution. A total observation time of 149 h offers unprecedented sensitivity. We use these data to explore several CO isotopologue line ratios in detail, study their radial (and azimuthal) trends and investigate whether changes in line ratios stem from changes in ISM properties such as gas temperatures, densities or chemical abundances. For example, we find negative radial trends for the $^{13}$CO}/$^{12}$CO, C$^{18}$O/$^{12}$CO and C$^{18}$O/$^{13}$CO line ratios in their J=1-0 transitions. We also find variations with local environment, such as higher $^{12}$CO(2-1)/(1-0) or $^{13}$CO/$^{12}$CO(1-0) line ratios in interarm regions compared to spiral arm regions. We propose that these aforementioned variations of CO line ratios are most likely due to a variation of the optical depth, while abundance variations due to selective nucleosynthesis on a galaxy-wide scale could also play a role. We also study the CO spectral line energy distribution (SLED) using archival JCMT $^{12}$CO(3-2) data and find a variation of the SLED shape with local environmental parameters further underlying changes in optical depth, gas temperatures or densities. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: 29 pages, 22 figures, accepted for publication in A&A

Journal ref: A&A 662, A89 (2022)

Planetary Nebula Luminosity Function distances for 19 galaxies observed by PHANGS-MUSE

Authors: Fabian Scheuermann, Kathryn Kreckel, Gagandeep S. Anand, Guillermo A. Blanc, Enrico Congiu, Francesco Santoro, Schuyler D. Van Dyk, Ashley T. Barnes, Frank Bigiel, Simon C. O. Glover, Brent Groves, Ralf S. Klessen, J. M. Diederik Kruijssen, Erik Rosolowsky, Eva Schinnerer, Andreas Schruba, Elizabeth J. Watkins, Thomas G. Williams

Abstract: We provide new planetary nebula luminosity function (PNLF) distances to 19 nearby spiral galaxies that were observed with VLT/MUSE by the PHANGS collaboration. Emission line ratios are used to separate planetary nebulae (PNe) from other bright [OIII] emitting sources like compact supernovae remnants (SNRs) or HII regions. While many studies have used narrowband imaging for this purpose, the detail… ▽ More We provide new planetary nebula luminosity function (PNLF) distances to 19 nearby spiral galaxies that were observed with VLT/MUSE by the PHANGS collaboration. Emission line ratios are used to separate planetary nebulae (PNe) from other bright [OIII] emitting sources like compact supernovae remnants (SNRs) or HII regions. While many studies have used narrowband imaging for this purpose, the detailed spectral line information provided by integral field unit (IFU) spectroscopy grants a more robust way of categorising different [OIII] emitters. We investigate the effects of SNR contamination on the PNLF and find that we would fail to classify all objects correctly, when limited to the same data narrowband imaging provides. However, the few misclassified objects usually do not fall on the bright end of the luminosity function, and only in three cases does the distance change by more than $1σ$. We find generally good agreement with literature values from other methods. Using metallicity constraints that have also been derived from the same IFU data, we revisit the PNLF zero point calibration. Over a range of $8.34 < 12 + \log(\mathrm{O}/\mathrm{H}) < 8.59$, our sample is consistent with a constant zero point and yields $M^*=-4.542^{+0.103}_{-0.059}\, \mathrm{mag}$, within $1σ$ of other literature values. MUSE pushes the limits of PNLF studies and makes galaxies beyond $20\, \mathrm{Mpc}$ accessible for this kind of analysis. This approach to the PNLF shows great promise for leveraging existing archival IFU data on nearby galaxies. △ Less

Submitted 12 January, 2022; originally announced January 2022.

Comments: 23 pages, 31 figures. Accepted for publication in MNRAS

Dependence of $X_{\rm CO}$ on metallicity, intensity, and spatial scale in a self-regulated interstellar medium

Authors: Chia-Yu Hu, Andreas Schruba, Amiel Sternberg, Ewine F. van Dishoeck

Abstract: We study the CO(1-0)-to-H$_2$ conversion factor ($X_{\rm CO}$) and the line ratio of CO(2-1)-to-CO(1-0) ($R_{21}$) across a wide range of metallicity ($0.1 \leq Z/Z_\odot \leq 3$) in high-resolution (~0.2 pc) hydrodynamical simulations of a self-regulated multiphase interstellar medium. We construct synthetic CO emission maps via radiative transfer and systematically vary the "observational" beam… ▽ More We study the CO(1-0)-to-H$_2$ conversion factor ($X_{\rm CO}$) and the line ratio of CO(2-1)-to-CO(1-0) ($R_{21}$) across a wide range of metallicity ($0.1 \leq Z/Z_\odot \leq 3$) in high-resolution (~0.2 pc) hydrodynamical simulations of a self-regulated multiphase interstellar medium. We construct synthetic CO emission maps via radiative transfer and systematically vary the "observational" beam size to quantify the scale dependence. We find that the kpc-scale $X_{\rm CO}$ can be over-estimated at low $Z$ if assuming steady-state chemistry or assuming that the star-forming gas is H$_2$-dominated. On parsec scales, $X_{\rm CO}$ varies by orders of magnitude from place to place, primarily driven by the transition from atomic carbon to CO. The pc-scale $X_{\rm CO}$ drops to the Milky Way value of $2\times 10^{20}\ {\rm cm^{-2}~(K~km~s^{-1})^{-1}}$ once dust shielding becomes effective, independent of $Z$. The CO lines become increasingly optically thin at lower $Z$, leading to a higher $R_{21}$. Most cloud area is filled by diffuse gas with high $X_{\rm CO}$ and low $R_{21}$, while most CO emission originates from dense gas with low $X_{\rm CO}$ and high $R_{21}$. Adopting a constant $X_{\rm CO}$ strongly over- (under-)estimates H$_2$ in dense (diffuse) gas. The line intensity negatively (positively) correlates with $X_{\rm CO}$ ($R_{21}$) as it is a proxy of column density (volume density). On large scales, $X_{\rm CO}$ and $R_{21}$ are dictated by beam averaging, and they are naturally biased towards values in dense gas. Our predicted $X_{\rm CO}$ is a multivariate function of $Z$, line intensity, and beam size, which can be used to more accurately infer the H$_2$ mass. △ Less

Submitted 23 May, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: Published in ApJ (typos in the abstract fixed). Code for interpolating Lagrangian (particle) data onto an adaptive mesh (and auxiliary scripts for RADMC-3D) available at https://github.com/huchiayu/ParticleGridMapper.jl

A 2-3 mm high-resolution molecular line survey towards the centre of the nearby spiral galaxy NGC 6946

Authors: Cosima Eibensteiner, Ashley T. Barnes, Frank Bigiel, Eva Schinnerer, Daizhong Liu, David S. Meier, Antonio Usero, Adam K. Leroy, Erik Rosolowsky, Johannes Puschnig, Ilin Lazar, Jérôme Pety, Laura A. Lopez, Eric Emsellem, Ivana Bešlić, Miguel Querejeta, Eric J. Murphy, Jakob den Brok, Andreas Schruba, Mélanie Chevance, Simon C. O. Glover, Yu Gao, Kathryn Grasha, Hamid Hassani, Jonathan D. Henshaw , et al. (8 additional authors not shown)

Abstract: The complex physical, kinematic, and chemical properties of galaxy centres make them interesting environments to examine with molecular line emission. We present new $2-4$" (${\sim}75{-}150$ pc at $7.7$ Mpc) observations at 2 and 3 mm covering the central $50$" (${\sim}1.9$ kpc) of the nearby double-barred spiral galaxy NGC 6946 obtained with the IRAM Plateau de Bure Interferometer. We detect spec… ▽ More The complex physical, kinematic, and chemical properties of galaxy centres make them interesting environments to examine with molecular line emission. We present new $2-4$" (${\sim}75{-}150$ pc at $7.7$ Mpc) observations at 2 and 3 mm covering the central $50$" (${\sim}1.9$ kpc) of the nearby double-barred spiral galaxy NGC 6946 obtained with the IRAM Plateau de Bure Interferometer. We detect spectral lines from ten molecules: CO, HCN, HCO$^+$, HNC, CS, HC$_3$N, N$_2$H$^+$, C$_2$H, CH$_3$OH, and H$_2$CO. We complemented these with published 1mm CO observations and 33 GHz continuum observations to explore the star formation rate surface density ${Σ_{\mathrm{SFR}}}$ on 150 pc scales. In this paper, we analyse regions associated with the inner bar of NGC 6946 $-$ the nuclear region (NUC), the northern (NBE), and southern inner bar end (SBE) and we focus on short-spacing corrected bulk (CO) and dense gas tracers (HCN, HCO$^+$, and HNC). We find that HCO$^+$ correlates best with ${Σ_{\mathrm{SFR}}}$, but the dense gas fraction ($f_{\mathrm{dense}}$) and star formation efficiency of the dense gas (${\mathrm{SFE_{dense}}}$) fits show different behaviours than expected from large-scale disc observations.The SBE has a higher ${Σ_{\mathrm{SFR}}}$, $f_{\mathrm{dense}}$, and shocked gas fraction than the NBE. We examine line ratio diagnostics and find a higher CO(2-1)/CO(1-0) ratio towards NBE than for the NUC. Moreover, comparison with existing extragalactic datasets suggests that using the HCN/HNC ratio to probe kinetic temperatures is not suitable on kiloparsec and sub-kiloparsec scales in extragalactic regions. Lastly, our study shows that the HCO$^+$/HCN ratio might not be a unique indicator to diagnose AGN activity in galaxies. △ Less

Submitted 10 January, 2022; v1 submitted 6 January, 2022; originally announced January 2022.

Comments: 38 pages, 23 figures, accepted for publication in A&A

Journal ref: A&A 659, A173 (2022)

The Gas-Star Formation Cycle in Nearby Star-forming Galaxies II. Resolved Distributions of CO and H$α$ Emission for 49 PHANGS Galaxies

Authors: Hsi-An Pan, Eva Schinnerer, Annie Hughes, Adam Leroy, Brent Groves, Ashley Thomas Barnes, Francesco Belfiore, Frank Bigiel, Guillermo A. Blanc, Yixian Cao, Melanie Chevance, Enrico Congiu, Daniel A. Dale, Cosima Eibensteiner, Eric Emsellem, Christopher M. Faesi, Simon C. O. Glover, Kathryn Grasha, Cinthya N. Herrera, I-Ting Ho, Ralf S. Klessen, J. M. Diederik Kruijssen, Philipp Lang, Daizhong Liu, Rebecca McElroy , et al. (14 additional authors not shown)

Abstract: The relative distribution of molecular gas and star formation in galaxies gives insight into the physical processes and timescales of the cycle between gas and stars. In this work, we track the relative spatial configuration of CO and H$α$ emission at high resolution in each of our galaxy targets, and use these measurements to quantify the distributions of regions in different evolutionary stages… ▽ More The relative distribution of molecular gas and star formation in galaxies gives insight into the physical processes and timescales of the cycle between gas and stars. In this work, we track the relative spatial configuration of CO and H$α$ emission at high resolution in each of our galaxy targets, and use these measurements to quantify the distributions of regions in different evolutionary stages of star formation: from molecular gas without star formation traced by H$α$ to star-forming gas, and to HII regions. The large sample, drawn from the Physics at High Angular resolution in Nearby GalaxieS ALMA and narrowband H$α$ (PHANGS-ALMA and PHANGS-H$α$) surveys, spans a wide range of stellar mass and morphological types, allowing us to investigate the dependencies of the gas-star formation cycle on global galaxy properties. At a resolution of 150 pc, the incidence of regions in different stages shows a dependence on stellar mass and Hubble type of galaxies over the radial range probed. Massive and/or earlier-type galaxies exhibit a significant reservoir of molecular gas without star formation traced by H$α$, while lower-mass galaxies harbor substantial HII regions that may have dispersed their birth clouds or formed from low-mass, more isolated clouds. Galactic structures add a further layer of complexity to relative distribution of CO and H$α$ emission. Trends between galaxy properties and distributions of gas traced by CO and H$α$ are visible only when the observed spatial scale is $\ll$ 500 pc, reflecting the critical resolution requirement to distinguish stages of star formation process. △ Less

Submitted 4 January, 2022; originally announced January 2022.

Comments: 60 pages, 22 figures, 9 tables, accepted for publication in The Astrophysical Journal

A tale of two DIGs: The relative role of HII regions and low-mass hot evolved stars in powering the diffuse ionised gas (DIG) in PHANGS-MUSE galaxies

Authors: Francesco Belfiore, Francesco Santoro, Brent Groves, Eva Schinnerer, Kathryn Kreckel, Simon C. O. Glover, Ralf S. Klessen, Eric Emsellem, Guillermo A. Blanc, Enrico Congiu, Ashley T. Barnes, Médéric Boquien, Mélanie Chevance, Daniel A. Dale, J. M. Diederik Kruijssen, Adam K. Leroy, Hsi-An Pan, Ismael Pessa, Andreas Schruba, Thomas G. Williams

Abstract: We use integral field spectroscopy from the PHANGS-MUSE survey, which resolves the ionised interstellar medium at ${\sim}50$ pc resolution in 19 nearby spiral galaxies, to study the origin of the diffuse ionised gas (DIG). We examine the physical conditions of the diffuse gas by first removing morphologically-defined HII regions and then binning the low-surface-brightness areas to achieve signific… ▽ More We use integral field spectroscopy from the PHANGS-MUSE survey, which resolves the ionised interstellar medium at ${\sim}50$ pc resolution in 19 nearby spiral galaxies, to study the origin of the diffuse ionised gas (DIG). We examine the physical conditions of the diffuse gas by first removing morphologically-defined HII regions and then binning the low-surface-brightness areas to achieve significant detections of the key nebular lines. A simple model for the leakage and propagation of ionising radiation from HII regions is able to reproduce the observed distribution of H$α$ in the DIG. Leaking radiation from HII regions also explains the observed decrease in line ratios of low-ionisation species ([SII]/H$α$, [NII]/H$α$ and [OI]/H$α$) with increasing H$α$ surface brightness ($Σ_{Hα}$). Emission from hot low-mass evolved stars, however, is required to explain: (1) the enhanced low-ionisation line ratios observed in the central regions of some galaxies; (2) the observed trends of a flat or decreasing [OIII]/H$β$ with $Σ_{Hα}$; and (3) the offset of some DIG regions from the locus of HII regions in the Baldwin-Phillips-Terlevich (BPT) diagram, extending into the area of low-ionisation (nuclear) emission-line regions (LI[N]ERs). Hot low-mass evolved stars make a small contribution to the energy budget of the DIG (2% of the galaxy-integrated H$α$ emission), but their harder spectra make them fundamental contributors to [OIII] emission. The DIG might result from a superposition of two components, an energetically dominant contribution from young stars and a more diffuse background of harder ionising photons from old stars. This unified framework bridges observations of the Milky Way DIG with LI(N)ER-like emission observed in nearby galaxy bulges. △ Less

Submitted 13 February, 2022; v1 submitted 29 November, 2021; originally announced November 2021.

Comments: accepted by A&A, posting version after proof corrections

Journal ref: A&A 659, A26 (2022)

PHANGS-MUSE: the HII region luminosity function of local star-forming galaxies

Authors: Francesco Santoro, Kathryn Kreckel, Francesco Belfiore, Brent Groves, Enrico Congiu, David A. Thilker, Guillermo A. Blanc, Eva Schinnerer, I-Ting Ho, J. M. Diederik Kruijssen, Sharon Meidt, Ralf S. Klessen, Andreas Schruba, Miguel Querejeta, Ismael Pessa, Mélanie Chevance, Jaeyeon Kim, Eric Emsellem, Rebecca McElroy, Ashley T. Barnes, Frank Bigiel, Médéric Boquien, Daniel A. Dale, Simon C. O. Glover, Kathryn Grasha , et al. (8 additional authors not shown)

Abstract: We use an unprecedented sample of about 23,000 HII regions detected at an average physical resolution of 67pc in the PHANGS-MUSE sample to study the extragalactic HII region Ha luminosity function (LF). Our observations probe the star-forming disk of 19 nearby spiral galaxies with low inclination and located close to the star formation main sequence at z=0. The mean LF slope $α$ in our sample is =… ▽ More We use an unprecedented sample of about 23,000 HII regions detected at an average physical resolution of 67pc in the PHANGS-MUSE sample to study the extragalactic HII region Ha luminosity function (LF). Our observations probe the star-forming disk of 19 nearby spiral galaxies with low inclination and located close to the star formation main sequence at z=0. The mean LF slope $α$ in our sample is =1.73 with a $σ$ of 0.15. We find that $α$ decreases with the galaxy's star formation rate surface density and argue that this is driven by an enhanced clustering of young stars at high gas surface densities. Looking at the HII regions within single galaxies we find that no significant variations occur between the LF of the inner and outer part of the star-forming disk, whereas the LF in the spiral arm areas is shallower than in the inter-arm areas for six out of the 13 galaxies with clearly visible spiral arms. We attribute these variations to the spiral arms increasing the molecular clouds arm--inter-arm mass contrast and find suggestive evidence that they are more evident for galaxies with stronger spiral arms. Furthermore, we find systematic variations in $α$ between samples of HII regions with high and low ionization parameter q and argue that they are driven by the aging of HII regions. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: Accepted in A&A. Main text: 19 pages, 11 Figures, 4 Tables. Catalog including 23301 HII regions will be available in electronic form at the CDS upon official publication

Journal ref: A&A 658, A188 (2022)

Molecular Gas Properties and CO-to-H2 Conversion Factors in the Central Kiloparsec of NGC 3351

Authors: Yu-Hsuan Teng, Karin M. Sandstrom, Jiayi Sun, Adam K. Leroy, L. Clifton Johnson, Alberto D. Bolatto, J. M. Diederik Kruijssen, Andreas Schruba, Antonio Usero, Ashley T. Barnes, Frank Bigiel, Guillermo A. Blanc, Brent Groves, Frank P. Israel, Daizhong Liu, Erik Rosolowsky, Eva Schinnerer, J. D. Smith, Fabian Walter

Abstract: The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is critical to studying molecular gas and star formation in galaxies. The value of $α_\rm{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 $\sim$kpc scales revealed low values of $α_\rm{CO}$ in the centers of some barred sp… ▽ More The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is critical to studying molecular gas and star formation in galaxies. The value of $α_\rm{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 $\sim$kpc scales revealed low values of $α_\rm{CO}$ in the centers of some barred spiral galaxies, including NGC 3351. We present new ALMA Band 3, 6, and 7 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales in the inner $\sim$2 kpc of NGC 3351. Using multi-line 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\times10^3$ $\rm{cm^{-3}}$ in the central ${\sim}$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\times10^{-4}$, our analysis yields $α_\rm{CO}{\sim}0.5{-}2.0$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$ with a decreasing trend with galactocentric radius in the central $\sim$1 kpc. The inflows show a substantially lower $α_\rm{CO} < 0.1$ $\rm{M_\odot\,(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 $α_\rm{CO}$ of ${\sim}1.5$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, which is similar to previous dust modeling based results at kpc scales. This suggests that low $α_\rm{CO}$ on kpc scales in the centers of some barred galaxies may be due to the contribution of low optical depth CO emission in bar-driven inflows. △ Less

Submitted 28 December, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: 31 pages, 17 figures; accepted for publication in ApJ; font size of text and figures adjusted

Investigating variations in the dust emissivity index in the Andromeda galaxy

Authors: G. Athikkat-Eknath, S. A. Eales, M. W. L. Smith, A. Schruba, K. A. Marsh, A. P. Whitworth

Abstract: Over the past decade, studies of dust in the Andromeda galaxy (M31) have shown radial variations in the dust emissivity index ($β$). Understanding the astrophysical reasons behind these radial variations may give clues about the chemical composition of dust grains, their physical structure, and the evolution of dust. We use $^{12}$CO(J=1-0) observations taken by the Combined Array for Research in… ▽ More Over the past decade, studies of dust in the Andromeda galaxy (M31) have shown radial variations in the dust emissivity index ($β$). Understanding the astrophysical reasons behind these radial variations may give clues about the chemical composition of dust grains, their physical structure, and the evolution of dust. We use $^{12}$CO(J=1-0) observations taken by the Combined Array for Research in Millimeter Astronomy (CARMA) and dust maps derived from \textit{Herschel} images, both with an angular resolution of 8" and spatial resolution of 30 pc, to study variations in $β$ across an area of $\approx$ 18.6 kpc$^2$ in M31. We extract sources, which we identify as molecular clouds, by applying the astrodendro algorithm to the $^{12}$CO and dust maps, which as a byproduct allows us to compare continuum emission from dust and CO emission as alternative ways of finding molecular clouds. We then use these catalogues to investigate whether there is evidence that $β$ is different inside and outside molecular clouds. Our results confirm the radial variations of $β$ seen in previous studies. However, we find little difference between the average $β$ inside molecular clouds compared to outside molecular clouds, in disagreement with models which predict an increase of $β$ in dense environments. Finally, we find some clouds traced by dust with very little CO which may be either clouds dominated by atomic gas or clouds of molecular gas that contain little CO. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 14 pages, 11 figures, accepted for publication in MNRAS

Bright, Relatively Isolated Star Clusters in PHANGS-HST Galaxies: Aperture Corrections, Quantitative Morphologies, and Comparison with Synthetic Stellar Population Models

Authors: Sinan Deger, Janice C. Lee, Bradley C. Whitmore, David A. Thilker, Médéric Boquien, Rupali Chandar, Daniel A. Dale, Leonardo Ubeda, Rick White, Kathryn Grasha, Simon C. O. Glover, Andreas Schruba, Ashley T. Barnes, Ralf Klessen, J. M. Diederik Kruijssen, Erik Rosolowsky, Thomas G. Williams

Abstract: Using PHANGS-HST NUV-U-B-V-I imaging of 17 nearby spiral galaxies, we study samples of star clusters and stellar associations, visually selected to be bright and relatively isolated, for three purposes: to compute aperture corrections for star cluster photometry, to explore the utility of quantitative morphologies in the analysis of clusters and associations, and to compare to synthetic stellar po… ▽ More Using PHANGS-HST NUV-U-B-V-I imaging of 17 nearby spiral galaxies, we study samples of star clusters and stellar associations, visually selected to be bright and relatively isolated, for three purposes: to compute aperture corrections for star cluster photometry, to explore the utility of quantitative morphologies in the analysis of clusters and associations, and to compare to synthetic stellar population models. We provide a technical summary of our procedures to determine aperture corrections, a standard step in the production of star cluster candidate catalogs, and compare to prior work. We also use this specialized sample to launch an analysis into the measurement of star cluster light profiles. We focus on one measure, $M_{20}$ (normalized second order moment of the brightest 20% of pixels), applied previously to study the morphologies of galaxies. We find that $M_{20}$ in combination with UB-VI colors, yields a parameter space where distinct loci are formed by single-peaked symmetric clusters, single-peaked asymmetric clusters, and multi-peaked associations. We discuss the potential applications for using $M_{20}$ to gain insight into the formation and evolution of clusters and associations. Finally, we compare the color distributions of this sample with various synthetic stellar population models. One finding is that the standard procedure of using a single-metallicity SSP track to fit the entire population of clusters in a given galaxy should be revisited, as the oldest globular clusters will be more metal-poor compared to clusters formed recently. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: 20 figures, 4 tables, 23 pages. Accepted for publication in MNRAS

The Two-Dimensional Metallicity Distribution and Mixing Scales of Nearby Galaxies

Authors: Thomas G. Williams, Kathryn Kreckel, Francesco Belfiore, Brent Groves, Karin Sandstrom, Francesco Santoro, Guillermo A. Blanc, Frank Bigiel, Médéric Boquien, Mélanie Chevance, Enrico Congiu, Eric Emsellem, Simon C. O. Glover, Kathryn Grasha, Ralf S. Klessen, Eric Koch, J. M. Diederik Kruijssen, Adam K. Leroy, Daizhong Liu, Sharon Meidt, Hsi-An Pan, Miguel Querejeta, Erik Rosolowsky, Toshiki Saito, Patricia Sánchez-Blázquez , et al. (3 additional authors not shown)

Abstract: Understanding the spatial distribution of metals within galaxies allows us to study the processes of chemical enrichment and mixing in the interstellar medium (ISM). In this work, we map the two-dimensional 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 part of the… ▽ More Understanding the spatial distribution of metals within galaxies allows us to study the processes of chemical enrichment and mixing in the interstellar medium (ISM). In this work, we map the two-dimensional 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 part of the PHANGS-MUSE survey. We find that 12 of our 19 galaxies show significant two-dimensional metallicity variation. Those without significant variations typically have fewer metallicity measurements, indicating this is due to the dearth of HII 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. △ Less

Submitted 21 December, 2021; v1 submitted 20 October, 2021; originally announced October 2021.

Comments: 20 pages, 14 Figures, 2 Appendices. Accepted for publication in MNRAS. Updated to include link to publicly available metallicity maps, have a good holiday and a restful new year!

Comparing the pre-SNe feedback and environmental pressures for 6000 HII regions across 19 nearby spiral galaxies

Authors: A. T. Barnes, S. C. O. Glover, K. Kreckel, E. C. Ostriker, F. Bigiel, F. Belfiore, I. Bešlić, G. A. Blanc, M. Chevance, D. A. Dale, O. Egorov, C. Eibensteiner, E. Emsellem, K. Grasha, B. A. Groves, R. S. Klessen, J. M. D. Kruijssen, A. K. Leroy, S. N. Longmore, L. Lopez, R. McElroy, S. E. Meidt, E. J. Murphy, E. Rosolowsky, T. Saito , et al. (6 additional authors not shown)

Abstract: The feedback from young stars (i.e. pre-supernova) is thought to play a crucial role in molecular cloud destruction. In this paper, we assess the feedback mechanisms acting within a sample of 5810 HII regions identified from the PHANGS-MUSE survey of 19 nearby ($<$ 20 Mpc) star-forming, main sequence spiral galaxies (log($M_\star$/M$_\odot$)= 9.4 $-$ 11). These optical spectroscopic maps are essen… ▽ More The feedback from young stars (i.e. pre-supernova) is thought to play a crucial role in molecular cloud destruction. In this paper, we assess the feedback mechanisms acting within a sample of 5810 HII regions identified from the PHANGS-MUSE survey of 19 nearby ($<$ 20 Mpc) star-forming, main sequence spiral galaxies (log($M_\star$/M$_\odot$)= 9.4 $-$ 11). These optical spectroscopic maps are essential to constrain the physical properties of the HII regions, which we use to investigate their internal pressure terms. We estimate the photoionised gas ($P_\mathrm{therm}$), direct radiation ($P_\mathrm{rad}$), and mechanical wind pressure ($P_\mathrm{wind}$), which we compare to the confining pressure of their host environment ($P_\mathrm{de}$). The HII regions remain unresolved within our ${\sim}50{-}100$ pc resolution observations, so we place upper ($P_\mathrm{max}$) and lower ($P_\mathrm{min}$) limits on each of the pressures by using a minimum (i.e. clumpy structure) and maximum (i.e. smooth structure) size, respectively. We find that the $P_\mathrm{max}$ measurements are broadly similar, and for $P_\mathrm{min}$ the $P_\mathrm{therm}$ is mildly dominant. We find that the majority of HII regions are over-pressured, $P_\mathrm{tot}/P_\mathrm{de} = (P_\mathrm{therm}+P_\mathrm{wind}+P_\mathrm{rad})/P_\mathrm{de} > 1$, and expanding, yet there is a small sample of compact HII regions with $P_\mathrm{tot,max}/P_\mathrm{de} < 1$ ($\sim$1% of the sample). These mostly reside in galaxy centres ($R_\mathrm{gal}<1$kpc), or, specifically, environments of high gas surface density; log($Σ_\mathrm{gas}/\mathrm{M_\odot} \mathrm{pc}^{-2}$)$\sim$2.5 (measured on kpc-scales). Lastly, we compare to a sample of literature measurements for $P_\mathrm{therm}$ and $P_\mathrm{rad}$ to investigate how dominant pressure term transitions over around 5dex in spatial dynamic range and 10 dex in pressure. △ Less

Submitted 25 October, 2021; v1 submitted 11 October, 2021; originally announced October 2021.

Comments: 28 pages total, 16+3 figures, and 3 tables. Accepted for publication in MNRAS

The PHANGS-MUSE survey -- Probing the chemo-dynamical evolution of disc galaxies

Authors: Eric Emsellem, Eva Schinnerer, Francesco Santoro, Francesco Belfiore, Ismael Pessa, Rebecca McElroy, Guillermo A. Blanc, Enrico Congiu, Brent Groves, I-Ting Ho, Kathryn Kreckel, Alessandro Razza, Patricia Sanchez-Blazquez, Oleg Egorov, Chris Faesi, Ralf S. Klessen, Adam K. Leroy, Sharon Meidt, Miguel Querejeta, Erik Rosolowsky, Fabian Scheuermann, Gagandeep S. Anand, Ashley T. Barnes, Ivana Bešlić, Frank Bigiel , et al. (23 additional authors not shown)

Abstract: We present the PHANGS-MUSE survey, a programme using the MUSE IFS at the ESO VLT to map 19 massive $(9.4 < \log(M_{*}/M_\odot) < 11.0)$ nearby (D < 20 Mpc) star-forming disc galaxies. The survey consists of 168 MUSE pointings (1'x1' each), a total of nearly 15 Million spectra, covering ~1.5 Million independent spectra. PHANGS-MUSE provides the first IFS view of star formation across different loca… ▽ More We present the PHANGS-MUSE survey, a programme using the MUSE IFS at the ESO VLT to map 19 massive $(9.4 < \log(M_{*}/M_\odot) < 11.0)$ nearby (D < 20 Mpc) star-forming disc galaxies. The survey consists of 168 MUSE pointings (1'x1' each), a total of nearly 15 Million spectra, covering ~1.5 Million independent spectra. PHANGS-MUSE provides the first IFS view of star formation across different local environments (including galaxy centres, bars, spiral arms) in external galaxies at a median resolution of 50~pc, better than the mean inter-cloud distance in the ionised interstellar medium. This `cloud-scale' resolution allows detailed demographics and characterisations of HII regions and other ionised nebulae. PHANGS-MUSE further delivers a unique view on the associated gas and stellar kinematics, and provides constraints on the star formation history. The PHANGS-MUSE survey is complemented by dedicated ALMA CO(2-1) and multi-band HST observations, therefore allowing us to probe the key stages of the star formation process from molecular clouds to HII regions and star clusters. This paper describes the scientific motivation, sample selection, observational strategy, data reduction and analysis process of the PHANGS-MUSE survey. We present our bespoke automated data-reduction framework, which is built on the reduction recipes provided by ESO, but additionally allows for mosaicking and homogenisation of the point spread function. We further present a detailed quality assessment and a brief illustration of the potential scientific applications of the large set of PHANGS-MUSE data products generated by our data analysis framework. The data cubes and analysis data products described in this paper represent the basis for the first PHANGS-MUSE public data release and are available in the ESO archive and via the Canadian Astronomy Data Centre. △ Less

Submitted 5 January, 2022; v1 submitted 7 October, 2021; originally announced October 2021.

Comments: 48 pages, 28 figures, 4 tables, accepted for pub in A&A, Dec 8, 2021. Data products released publicly at this address https://archive.eso.org/scienceportal/home?data_collection=PHANGS

Journal ref: A&A 659, A191 (2022)

The HASHTAG project: The First Submillimeter Images of the Andromeda Galaxy from the Ground

Authors: Matthew W. L. Smith, Stephen A. Eales, Thomas G. Williams, Bumhyun Lee, Zongnan Li, Pauline Barmby, Martin Bureau, Scott Chapman, Brian S. Cho, Aeree Chung, Eun Jung Chung, Hui-Hsuan Chung, Christopher J. R. Clark, David L. Clements, Timothy A. Davis, Ilse De Looze, David J. Eden, Gayathri Athikkat-Eknath, George P. Ford, Yu Gao, Walter Gear, Haley L. Gomez, Richard de Grijs, Jinhua He, Luis C. Ho , et al. (24 additional authors not shown)

Abstract: Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing… ▽ More Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing the galaxy's large-scale structure. We have developed a technique for producing high-resolution submillimeter images of galaxies of large angular size by using the telescope on the ground to determine the small-scale structure (the large Fourier components) and a space telescope (Herschel or Planck) to determine the large-scale structure (the small Fourier components). Using this technique, we are carrying out the HARP and SCUBA-2 High Resolution Terahertz Andromeda Galaxy Survey (HASHTAG), an international Large Program on the James Clerk Maxwell Telescope, with one aim being to produce the first high-fidelity high-resolution submillimeter images of Andromeda. In this paper, we describe the survey, the method we have developed for combining the space-based and ground-based data, and present the first HASHTAG images of Andromeda at 450 and 850um. We also have created a method to predict the CO(J=3-2) line flux across M31, which contaminates the 850um band. We find that while normally the contamination is below our sensitivity limit, the contamination can be significant (up to 28%) in a few of the brightest regions of the 10 kpc ring. We therefore also provide images with the predicted line emission removed. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Comments: 26 pages, 19 figures. Submitted to ApJS June 2021, Accepted September 2021

Low-J CO Line Ratios From Single Dish CO Mapping Surveys and PHANGS-ALMA

Authors: Adam K. Leroy, Erik Rosolowsky, Antonio Usero, Karin Sandstrom, Eva Schinnerer, Andreas Schruba, Alberto D. Bolatto, Jiayi Sun, Ashley. T. Barnes, Francesco Belfiore, Frank Bigiel, Jakob S. den Brok, Yixian Cao, I-Da Chiang, Mélanie Chevance, Daniel A. Dale, Cosima Eibensteiner, Christopher M. Faesi, Simon C. O. Glover, Annie Hughes, Maria J. Jiménez Donaire, Ralf S. Klessen, Eric W. Koch, J. M. Diederik Kruijssen, Daizhong Liu , et al. (10 additional authors not shown)

Abstract: We measure the low-J CO line ratio R21=CO(2-1)/CO(1-0), R32=CO(3-2)/CO(2-1), and R31 = 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 JCMT NGLS and APEX LASMA mapping. Altogether this yields 76, 47, and 29 maps of R21,… ▽ More We measure the low-J CO line ratio R21=CO(2-1)/CO(1-0), R32=CO(3-2)/CO(2-1), and R31 = 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 JCMT NGLS and APEX LASMA mapping. Altogether this yields 76, 47, and 29 maps of R21, R32, and R31 at 20" \sim 1.3 kpc resolution, covering 43, 34, and 20 galaxies. Disk galaxies with high stellar mass, log10 M_* [Msun]=10.25-11 and star formation rate, SFR=1-5 Msun/yr, dominate the sample. We find galaxy-integrated mean values and 16%-84% range of R21 = 0.65 (0.50-0.83), R32=0.50 (0.23-0.59), and R31=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 \sim 0.18^{+0.09}_{-0.14} dex for R21, 0.27^{+0.13}_{-0.15} dex for R31, and 0.08^{+0.11}_{-0.09} dex for R32. All three line ratios anti-correlate with galactocentric radius and positively correlate with the local star formation rate surface density and specific star formation rate, 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 ISM. Because the line ratios are expected to anti-correlate with the CO(1-0)-to-H_2 conversion factor, alphaCO^(1-0), these results have general implications for the interpretation of CO emission from galaxies. △ Less

Submitted 23 September, 2021; originally announced September 2021.

Comments: 34 pages; 12 figures; accepted for publication in the Astrophysical Journal; data tables available http://www.astronomy.ohio-state.edu/~leroy.42/int_corat_tab.ecsv and http://www.astronomy.ohio-state.edu/~leroy.42/co_model.ecsv before publication

The impact of pre-supernova feedback and its dependence on environment

Authors: Anna F. Mcleod, Ahmad A. Ali, Mélanie Chevance, Lorenza Della Bruna, Andreas Schruba, Heloise F. Stevance, Angela Adamo, J. M. Diederik Kruijssen, Steven N. Longmore, Daniel R. Weisz, Peter Zeidler

Abstract: Integral field units enable resolved studies of a large number of star-forming regions across entire nearby galaxies, providing insight on the conversion of gas into stars and the feedback from the emerging stellar populations over unprecedented dynamic ranges in terms of spatial scale, star-forming region properties, and environments. We use the VLT/MUSE legacy data set covering the central $35$… ▽ More Integral field units enable resolved studies of a large number of star-forming regions across entire nearby galaxies, providing insight on the conversion of gas into stars and the feedback from the emerging stellar populations over unprecedented dynamic ranges in terms of spatial scale, star-forming region properties, and environments. We use the VLT/MUSE legacy data set covering the central $35$ arcmin$^{2}$ (${\sim}12$ kpc$^{2}$) of the nearby galaxy NGC 300 to quantify the effect of stellar feedback as a function of the local galactic environment. We extract spectra from emission line regions identified within dendrograms, combine emission line ratios and line widths to distinguish between HII regions, planetary nebulae, and supernova remnants, and compute their ionised gas properties, gas-phase oxygen abundances, and feedback-related pressure terms. For the HII regions, we find that the direct radiation pressure ($P_\mathrm{dir}$) and the pressure of the ionised gas ($P_{HII}$) weakly increase towards larger galactocentric radii, i.e. along the galaxy's (negative) abundance and (positive) extinction gradients. While the increase of $P_{HII}$ with galactocentric radius is likely due to higher photon fluxes from lower-metallicity stellar populations, we find that the increase of $P_\mathrm{dir}$ is likely driven by the combination of higher photon fluxes and enhanced dust content at larger galactocentric radii. In light of the above, we investigate the effect of increased pre-supernova feedback at larger galactocentric distances (lower metallicities and increased dust mass surface density) on the ISM, finding that supernovae at lower metallicities expand into lower-density environments, thereby enhancing the impact of supernova feedback. △ Less

Submitted 17 September, 2021; originally announced September 2021.

Comments: Accepted for publication in MNRAS

Stellar structures, molecular gas, and star formation across the PHANGS sample of nearby galaxies

Authors: M. Querejeta, E. Schinnerer, S. Meidt, J. Sun, A. K. Leroy, E. Emsellem, R. S. Klessen, J. C. Munoz-Mateos, H. Salo, E. Laurikainen, I. Beslic, G. A. Blanc, M. Chevance, D. A. Dale, C. Eibensteiner, C. Faesi, A. Garcia-Rodriguez, S. C. O. Glover, K. Grasha, J. Henshaw, C. Herrera, A. Hughes, K. Kreckel, J. M. D. Kruijssen, D. Liu , et al. (10 additional authors not shown)

Abstract: We identify stellar structures in the PHANGS sample of 74 nearby galaxies and construct morphological masks of sub-galactic environments based on Spitzer 3.6 micron images. At the simplest level, we distinguish centres, bars, spiral arms, interarm and discs without strong spirals. Slightly more sophisticated masks include rings and lenses, publicly released but not explicitly used in this paper. W… ▽ More We identify stellar structures in the PHANGS sample of 74 nearby galaxies and construct morphological masks of sub-galactic environments based on Spitzer 3.6 micron images. At the simplest level, we distinguish centres, bars, spiral arms, interarm and discs without strong spirals. Slightly more sophisticated masks include rings and lenses, publicly released but not explicitly used in this paper. We examine trends using PHANGS-ALMA CO(2-1) intensity maps and tracers of star formation. The interarm regions and discs without strong spirals dominate in area, whereas molecular gas and star formation are quite evenly distributed among the five basic environments. We reproduce the molecular Kennicutt-Schmidt relation with a slope compatible with unity within the uncertainties, without significant slope differences among environments. In contrast to early studies, we find that bars are not always deserts devoid of gas and star formation, but instead they show large diversity. Similarly, spiral arms do not account for most of the gas and star formation in disc galaxies, and they do not have shorter depletion times than the interarm regions. Spiral arms accumulate gas and star formation, without systematically boosting the star formation efficiency. Centres harbour remarkably high surface densities and on average shorter depletion times than other environments. Centres of barred galaxies show higher surface densities and wider distributions compared to the outer disc; yet, depletion times are similar to unbarred galaxies, suggesting highly intermittent periods of star formation when bars episodically drive gas inflow, without enhancing the central star formation efficiency permanently. In conclusion, we provide quantitative evidence that stellar structures in galaxies strongly affect the organisation of molecular gas and star formation, but their impact on star formation efficiency is more subtle. △ Less

Submitted 7 October, 2021; v1 submitted 9 September, 2021; originally announced September 2021.

Comments: 28 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 656, A133 (2021)

Rotation Curves in z~1-2 Star-Forming Disks: Comparison of Dark Matter Fractions and Disk Properties for Different Fitting Methods

Authors: S. H. Price, T. T. Shimizu, R. Genzel, H. Übler, N. M. Förster Schreiber, L. J. Tacconi, R. I. Davies, R. T. Coogan, D. Lutz, S. Wuyts, E. Wisnioski, A. Nestor, A. Sternberg, A. Burkert, R. Bender, A. Contursi, R. L. Davies, R. Herrera-Camus, M. -J. Lee, T. Naab, R. Neri, A. Renzini, R. Saglia, A. Schruba, K. Schuster

Abstract: We present a follow-up analysis examining the dynamics and structures of 41 massive, large star-forming galaxies at z~0.67-2.45 using both ionized and molecular gas kinematics. We fit the galaxy dynamics with models consisting of a bulge, a thick, turbulent disk, and a NFW dark matter halo, using code that fully forward models the kinematics, including all observational and instrumental effects. W… ▽ More We present a follow-up analysis examining the dynamics and structures of 41 massive, large star-forming galaxies at z~0.67-2.45 using both ionized and molecular gas kinematics. We fit the galaxy dynamics with models consisting of a bulge, a thick, turbulent disk, and a NFW dark matter halo, using code that fully forward models the kinematics, including all observational and instrumental effects. We explore the parameter space using Markov Chain Monte Carlo (MCMC) sampling, including priors based on stellar and gas masses and disk sizes. We fit the full sample using extracted 1D kinematic profiles. For a subset of 14 well-resolved galaxies, we also fit the 2D kinematics. The MCMC approach robustly confirms the results from least-squares fitting presented in Paper I (Genzel et al. 2020): the sample galaxies tend to be baryon-rich on galactic scales (within one effective radius). The 1D and 2D MCMC results are also in good agreement for the subset, demonstrating that much of the galaxy dynamical information is captured along the major axis. The 2D kinematics are more affected by the presence of non-circular motions, which we illustrate by constructing a toy model with constant inflow for one galaxy that exhibits residual signatures consistent with radial motions. This analysis, together with results from Paper I and other studies, strengthens the finding that massive, star-forming galaxies at z~1-2 are baryon-dominated on galactic scales, with lower dark matter fractions towards higher baryonic surface densities. Finally, we present details of the kinematic fitting code used in this analysis. △ Less

Submitted 6 September, 2021; originally announced September 2021.

Comments: Accepted for publication in ApJ (23 pages, 8 figures, 5 tables)

Star Cluster Classification in the PHANGS-HST Survey: Comparison between Human and Machine Learning Approaches

Authors: Bradley C. Whitmore, Janice C. Lee, Rupali Chandar, David A. Thilker, Stephen Hannon, Wei Wei, E. A. Huerta, Frank Bigiel, Médéric Boquien, Mélanie Chevance, Daniel A. Dale, Sinan Deger, Kathryn Grasha, Ralf S. Klessen, J. M. Diederik Kruijssen, Kirsten L. Larson, Angus Mok, Erik Rosolowsky, Eva Schinnerer, Andreas Schruba, Leonardo Ubeda, Schuyler D. Van Dyk, Elizabeth Watkins, Thomas Williams

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 (N… ▽ More 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$\%$ for Class 1 clusters (symmetric, centrally concentrated) and 40 to 70$\%$ for Class 2 clusters (asymmetric, centrally concentrated). If Class 1 and 2 are considered together the agreement is 82 $\pm$ 3$\%$. 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. △ Less

Submitted 27 July, 2021; originally announced July 2021.

Comments: 28 pages, 25 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society. Version with full resolution figures found at https://sites.google.com/view/phangs/publications

Frequency and nature of central molecular outflows in nearby star-forming disk galaxies

Authors: Sophia K. Stuber, Toshiki Saito, Eva Schinnerer, Eric Emsellem, Miguel Querejeta, Thomas G. Williams, Ashley T. Barnes, Frank Bigiel, Guillermo Blanc, Daniel A. Dale, Kathryn Grasha, Ralf Klessen, J. M. Diederik Kruijssen, Adam K. Leroy, Sharon Meidt, Hsi-An Pan, Erik Rosolowsky, Andreas Schruba, Jiayi Sun, Antonio Usero

Abstract: Central molecular outflows in spiral galaxies are assumed to modulate their host galaxy's star formation rate by removing gas from the inner region of the galaxy. Outflows consisting of different gas phases appear to be a common feature in local galaxies, yet, little is known about the frequency of molecular outflows in main sequence galaxies in the nearby universe. We develop a rigorous set of se… ▽ More Central molecular outflows in spiral galaxies are assumed to modulate their host galaxy's star formation rate by removing gas from the inner region of the galaxy. Outflows consisting of different gas phases appear to be a common feature in local galaxies, yet, little is known about the frequency of molecular outflows in main sequence galaxies in the nearby universe. We develop a rigorous set of selection criteria, which allow the reliable identification of outflows in large samples of galaxies. Our criteria make use of central spectra, position-velocity diagrams and velocity-integrated intensity maps (line-wing maps). We use this method on high-angular resolution CO(2-1) observations from the PHANGS-ALMA survey, which provides observations of the molecular gas for a homogeneous sample of 90 nearby main sequence galaxies at a resolution of ${\sim}100\,$pc. We find correlations between the assigned outflow confidence and stellar mass or global star formation rate (SFR). We determine the frequency of central molecular outflows to be $25\pm2$% considering all outflow candidates, or $20\pm2$% for secure outflows only. Our resulting outflow candidate sample of $16{-}20$ galaxies shows an overall enhanced fraction of active galactic nuclei (AGN) (50%) and bars (89%) compared to the full sample (galaxies with AGN: 24%, with bar: 61%). We extend the trend between mass outflow rates and SFR known for high outflow rates down to lower values ($\log_{10}{\dot{\rm M}_{\rm out}}\,[\mathrm{M}_\odot~\mathrm{yr}^{-1}]<0$). Mass loading factors are of order unity, indicating that these outflows are not efficient in quenching the SFR in main sequence galaxies. △ Less

Submitted 3 September, 2021; v1 submitted 23 July, 2021; originally announced July 2021.

Comments: 18 pages, 11 pages appendix, 13 figures, 6 tables, accepted for publication in A&A

Journal ref: A&A 653, A172 (2021)

PHANGS-HST: New Methods for Star Cluster Identification in Nearby Galaxies

Authors: David A. Thilker, Bradley C. Whitmore, Janice C. Lee, Sinan Deger, Rupali Chandar, Kirsten L. Larson, Stephen Hannon, Leonardo Ubeda, Daniel A. Dale, Simon C. O. Glover, Kathryn Grasha, Ralf S. Klessen, J. M. Diederik Kruijssen, Erik Rosolowsky, Andreas Schruba, Richard L. White, Thomas G. Williams

Abstract: We present an innovative and widely applicable approach for the detection and classification of stellar clusters, developed for the PHANGS-HST Treasury Program, an $NUV$-to-$I$ band imaging campaign of 38 spiral galaxies. Our pipeline first generates a unified master source list for stars and candidate clusters, to enable a self-consistent inventory of all star formation products. To distinguish c… ▽ More We present an innovative and widely applicable approach for the detection and classification of stellar clusters, developed for the PHANGS-HST Treasury Program, an $NUV$-to-$I$ band imaging campaign of 38 spiral galaxies. Our pipeline first generates a unified master source list for stars and candidate clusters, to enable a self-consistent inventory of all star formation products. To distinguish cluster candidates from stars, we introduce the Multiple Concentration Index (MCI) parameter, and measure inner and outer MCIs to probe morphology in more detail than with a single, standard concentration index (CI). We improve upon cluster candidate selection, jointly basing our criteria on expectations for MCI derived from synthetic cluster populations and existing cluster catalogues, yielding model and semi-empirical selection regions (respectively). Selection purity (confirmed clusters versus candidates, assessed via human-based classification) is high (up to 70$\%$) for moderately luminous sources in the semi-empirical selection region, and somewhat lower overall (outside the region or fainter). The number of candidates rises steeply with decreasing luminosity, but pipeline-integrated Machine Learning (ML) classification prevents this from being problematic. We quantify the performance of our PHANGS-HST methods in comparison to LEGUS for a sample of four galaxies in common to both surveys, finding overall agreement with $50-75\%$ of human verified star clusters appearing in both catalogues, but also subtle differences attributable to specific choices adopted by each project. The PHANGS-HST ML-classified Class 1 or 2 catalogues reach $\sim1$ magnitude fainter, $\sim2\times$ lower stellar mass, and are $2{-}5\times$ larger in number, than attained in the human classified samples. △ Less

Submitted 29 October, 2021; v1 submitted 24 June, 2021; originally announced June 2021.

Comments: 36 pages, 18 figures, updated to v2, accepted to MNRAS

Dense molecular gas properties on 100 pc scales across the disc of NGC 3627

Authors: I. Bešlić, A. T. Barnes, F. Bigiel, J. Puschnig, J. Pety, C. Herrera Contreras, A. K. Leroy, A. Usero, E. Schinnerer, S. E. Meidt, E. Emsellem, A. Hughes, C. Faesi, K. Kreckel, F. M. C. Belfiore, M. Chevance, J. S. den Brok, C. Eibensteiner, S. C. O. Glover, K. Grasha, M. J. Jimenez-Donaire, R. S. Klessen, J. M. D. Kruijssen, D. Liu, I. Pessa , et al. (7 additional authors not shown)

Abstract: It is still poorly constrained how the densest phase of the interstellar medium varies across galactic environment. A large observing time is required to recover significant emission from dense molecular gas at high spatial resolution, and to cover a large dynamic range of extragalactic disc environments. We present new NOrthern Extended Millimeter Array (NOEMA) observations of a range of high cri… ▽ More It is still poorly constrained how the densest phase of the interstellar medium varies across galactic environment. A large observing time is required to recover significant emission from dense molecular gas at high spatial resolution, and to cover a large dynamic range of extragalactic disc environments. We present new NOrthern Extended Millimeter Array (NOEMA) observations of a range of high critical density molecular tracers (HCN, HNC, HCO+) and CO isotopologues (13CO, C18O) towards the nearby (11.3 Mpc), strongly barred galaxy NGC 3627. These observations represent the current highest angular resolution (1.85"; 100 pc) map of dense gas tracers across a disc of a nearby spiral galaxy, which we use here to assess the properties of the dense molecular gas, and their variation as a function of galactocentric radius, molecular gas, and star formation. We find that the HCN(1-0)/CO(2-1) integrated intensity ratio does not correlate with the amount of recent star formation. Instead, the HCN(1-0)/CO(2-1) ratio depends on the galactic environment, with differences between the galaxy centre, bar, and bar end regions. The dense gas in the central 600 pc appears to produce stars less efficiently despite containing a higher fraction of dense molecular gas than the bar ends where the star formation is enhanced. In assessing the dynamics of the dense gas, we find the HCN(1-0) and HCO+(1-0) emission lines showing multiple components towards regions in the bar ends that correspond to previously identified features in CO emission. These features are co-spatial with peaks of Halpha emission, which highlights that the complex dynamics of this bar end region could be linked to local enhancements in the star formation. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Comments: 25 pages, 12 figures; accepted for publication in MNRAS

The centres of M83 and the Milky Way: opposite extremes of a common star formation cycle

Authors: Daniel Callanan, Steven N. Longmore, J. M. Diederik Kruijssen, Andreas Schruba, Adam Ginsburg, Mark R. Krumholz, Nate Bastian, Joao Alves, Jonathan D. Henshaw, Johan H. Knapen, Melanie Chevance

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-alpha 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 g… ▽ More 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-alpha 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", 12pc) 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 localised, 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) modeling and interpreting observations must avoid averaging over large spatial areas or timescales, and (ii) understanding the multi-scale processes controlling these cycles requires comparing snapshots of a statistical sample of galaxies in different evolutionary stages. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: 31 pages, 22 figures, accepted by MNRAS

Star formation scaling relations at ~100 pc from PHANGS: Impact of completeness and spatial scale

Authors: I. Pessa, E. Schinnerer, F. Belfiore, E. Emsellem, A. K. Leroy, A. Schruba, J. M. D. Kruijssen, H. -A. Pan, G. A. Blanc, P. Sanchez-Blazquez, F. Bigiel, M. Chevance, E. Congiu, D. Dale, C. M. Faesi, S. C. O. Glover, K. Grasha, B. Groves, I. Ho, M. Jiménez-Donaire, R. Klessen, K. Kreckel, E. W. Koch, D. Liu, S. Meidt , et al. (9 additional authors not shown)

Abstract: Aims: The complexity of star formation at the physical scale of molecular clouds is not yet fully understood. We investigate the mechanisms regulating the formation of stars in different environments within nearby star-forming galaxies from the PHANGS sample. Methods: Integral field spectroscopic data and radio-interferometric observations of 18 galaxies were combined to explore the existence of t… ▽ More Aims: The complexity of star formation at the physical scale of molecular clouds is not yet fully understood. We investigate the mechanisms regulating the formation of stars in different environments within nearby star-forming galaxies from the PHANGS sample. Methods: Integral field spectroscopic data and radio-interferometric observations of 18 galaxies were combined to explore the existence of the resolved star formation main sequence (rSFMS), resolved Kennicutt-Schmidt relation (rKS), and resolved molecular gas main sequence (rMGMS), and we derived their slope and scatter at spatial resolutions from 100 pc to 1 kpc (under various assumptions). Results: All three relations were recovered at the highest spatial resolution (100 pc). Furthermore, significant variations in these scaling relations were observed across different galactic environments. The exclusion of non-detections has a systematic impact on the inferred slope as a function of the spatial scale. Finally, the scatter of the $Σ_\mathrm{mol. gas + stellar}$ versus $Σ_\mathrm{SFR}$ correlation is smaller than that of the rSFMS, but higher than that found for the rKS. Conclusions: The rMGMS has the tightest relation at a spatial scale of 100 pc (scatter of 0.34 dex), followed by the rKS (0.41 dex) and then the rSFMS (0.51 dex). This is consistent with expectations from the timescales involved in the evolutionary cycle of molecular clouds. Surprisingly, the rKS shows the least variation across galaxies and environments, suggesting a tight link between molecular gas and subsequent star formation. The scatter of the three relations decreases at lower spatial resolutions, with the rKS being the tightest (0.27 dex) at a spatial scale of 1 kpc. Variation in the slope of the rSFMS among galaxies is partially due to different detection fractions of $Σ_\mathrm{SFR}$ with respect to $Σ_\mathrm{stellar}$. △ Less

Submitted 23 June, 2021; v1 submitted 19 April, 2021; originally announced April 2021.

Comments: 27 pages, 29 figures, Accepted in A&A

Journal ref: A&A 650, A134 (2021)

PHANGS-ALMA: Arcsecond CO(2-1) Imaging of Nearby Star-Forming Galaxies

Authors: Adam K. Leroy, Eva Schinnerer, Annie Hughes, Erik Rosolowsky, Jérôme Pety, Andreas Schruba, Antonio Usero, Guillermo A. Blanc, Mélanie Chevance, Eric Emsellem, Christopher M. Faesi, Cinthya N. Herrera, Daizhong Liu, Sharon E. Meidt, Miguel Querejeta, Toshiki Saito, Karin M. Sandstrom, Jiayi Sun, Thomas G. Williams, Gagandeep S. Anand, Ashley T. Barnes, Erica A. Behrens, Francesco Belfiore, Samantha M. Benincasa, Ivana Bešlić , et al. (47 additional authors not shown)

Abstract: We present PHANGS-ALMA, the first survey to map CO J=2-1 line emission at ~1" ~ 100pc spatial resolution from a representative sample of 90 nearby (d<~20 Mpc) galaxies that lie on or near the z=0 "main sequence" of star-forming galaxies. CO line emission traces the bulk distribution of molecular gas, which is the cold, star-forming phase of the interstellar medium. At the resolution achieved by PH… ▽ More We present PHANGS-ALMA, the first survey to map CO J=2-1 line emission at ~1" ~ 100pc spatial resolution from a representative sample of 90 nearby (d<~20 Mpc) galaxies that lie on or near the z=0 "main sequence" of star-forming galaxies. CO line emission traces the bulk distribution of molecular gas, which is the cold, star-forming phase of the interstellar medium. At the resolution achieved by PHANGS-ALMA, each beam reaches the size of a typical individual giant molecular cloud (GMC), so that these data can be used to measure the demographics, life-cycle, and physical state of molecular clouds across the population of galaxies where the majority of stars form at z=0. This paper describes the scientific motivation and background for the survey, sample selection, global properties of the targets, ALMA observations, and characteristics of the delivered ALMA data and derived data products. As the ALMA sample serves as the parent sample for parallel surveys with VLT/MUSE, HST, AstroSat, VLA, and other facilities, we include a detailed discussion of the sample selection. We detail the estimation of galaxy mass, size, star formation rate, CO luminosity, and other properties, compare estimates using different systems and provide best-estimate integrated measurements for each target. We also report the design and execution of the ALMA observations, which combine a Cycle~5 Large Program, a series of smaller programs, and archival observations. Finally, we present the first 1" resolution atlas of CO emission from nearby galaxies and describe the properties and contents of the first PHANGS-ALMA public data release. △ Less

Submitted 28 April, 2021; v1 submitted 15 April, 2021; originally announced April 2021.

Comments: 76 pages, 33 figures. Accepted for publication in the Astrophysical Journal Supplement series. Full resolution version and the image atlas to appear as a figure set in the published version can be found https://sites.google.com/view/phangs/publications . Data release coming soon to the ALMA archive and CADC temporarily available at http://phangs.org/data

PHANGS-ALMA Data Processing and Pipeline

Authors: Adam K. Leroy, Annie Hughes, Daizhong Liu, Jerome Pety, Erik Rosolowsky, Toshiki Saito, Eva Schinnerer, Andreas Schruba, Antonio Usero, Christopher M. Faesi, Cinthya N. Herrera, Melanie Chevance, Alexander P. S. Hygate, Amanda A. Kepley, Eric W. Koch, Miguel Querejeta, Kazimierz Sliwa, David Will, Christine D. Wilson, Gagandeep S. Anand, Ashley Barnes, Francesco Belfiore, Ivana Beslic, Frank Bigiel, Guillermo A. Blanc , et al. (43 additional authors not shown)

Abstract: We describe the processing of the PHANGS-ALMA survey and present the PHANGS-ALMA pipeline, a public software package that processes calibrated interferometric and total power data into science-ready data products. PHANGS-ALMA is a large, high-resolution survey of CO J=2-1 emission from nearby galaxies. The observations combine ALMA's main 12-m array, the 7-m array, and total power observations and… ▽ More We describe the processing of the PHANGS-ALMA survey and present the PHANGS-ALMA pipeline, a public software package that processes calibrated interferometric and total power data into science-ready data products. PHANGS-ALMA is a large, high-resolution survey of CO J=2-1 emission from nearby galaxies. The observations combine ALMA's main 12-m array, the 7-m array, and total power observations and use mosaics of dozens to hundreds of individual pointings. We describe the processing of the u-v data, imaging and deconvolution, linear mosaicking, combining interferometer and total power data, noise estimation, masking, data product creation, and quality assurance. Our pipeline has a general design and can also be applied to VLA and ALMA observations of other spectral lines and continuum emission. We highlight our recipe for deconvolution of complex spectral line observations, which combines multiscale clean, single scale clean, and automatic mask generation in a way that appears robust and effective. We also emphasize our two-track approach to masking and data product creation. We construct one set of "broadly masked" data products, which have high completeness but significant contamination by noise, and another set of "strictly masked" data products, which have high confidence but exclude faint, low signal-to-noise emission. Our quality assurance tests, supported by simulations, demonstrate that 12-m+7-m deconvolved data recover a total flux that is significantly closer to the total power flux than the 7-m deconvolved data alone. In the appendices, we measure the stability of the ALMA total power calibration in PHANGS--ALMA and test the performance of popular short-spacing correction algorithms. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: Accepted for publication in the Astrophysical Journal Supplement series. 65 pages, 33 figures. Software available at https://github.com/akleroy/phangs_imaging_scripts . For a full resolution version see https://sites.google.com/view/phangs/publications

A lack of constraints on the cold opaque HI mass: HI spectra in M31 and M33 prefer multi-component models over a single cold opaque component

Authors: Eric W. Koch, Erik W. Rosolowsky, Adam K. Leroy, Jeremy Chastenet, I-Da Chiang, Julianne Dalcanton, Amanda A. Kepley, Karin M. Sandstrom, Andreas Schruba, Snezana Stanimirovic, Dyas Utomo, Thomas G. Williams

Abstract: Previous work has argued that atomic gas mass estimates of galaxies from 21 cm HI emission are systematically low due to a cold opaque atomic gas component. If true, this opaque component necessitates a ~35% correction factor relative to the mass from assuming optically-thin HI emission. These mass corrections are based on fitting HI spectra with a single opaque component model that produces a dis… ▽ More Previous work has argued that atomic gas mass estimates of galaxies from 21 cm HI emission are systematically low due to a cold opaque atomic gas component. If true, this opaque component necessitates a ~35% correction factor relative to the mass from assuming optically-thin HI emission. These mass corrections are based on fitting HI 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 HI VLA observations of M31 and M33 to test if these top-hat profiles are instead superpositions of multiple HI components along the line-of-sight. We fit both models and find that >80% of the spectra strongly prefer a multi-component Gaussian model while <2% 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 HI. Our findings are enabled by the improved spectral resolution (0.42 km/s), 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 HI surface density with the dust surface density and extinction has significantly more scatter when the inferred HI 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 HI mass is best constrained by HI absorption studies. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: Accepted in MNRAS. 25 pg, 13 figures

The organization of cloud-scale gas density structure: high resolution CO vs. 3.6 $μ$m brightness contrasts in nearby galaxies

Authors: Sharon E. Meidt, Adam K. Leroy, Miguel Querejeta, Eva Schinnerer, Jiayi Sun, Arjen van der Wel, Eric Emsellem, Jonathan Henshaw, Annie Hughes, J. M. Diederik Kruijssen, Erik Rosolowsky, Andreas Schruba, Ashley Barnes, Frank Bigiel, Guillermo A. Blanc, Melanie Chevance, Yixian Cao, Daniel A. Dale, Christopher Faesi, Simon C. O. Glover, Kathryn Grasha, Brent Groves, Cinthya Herrera, Ralf S. Klessen, Kathryn Kreckel , et al. (7 additional authors not shown)

Abstract: In this paper we examine the factors that shape the distribution of molecular gas surface densities on the 150 pc scale across 67 morphologically diverse star-forming galaxies in the PHANGS-ALMA CO (2-1) survey. Dividing each galaxy into radial bins, we measure molecular gas surface density contrasts, defined here as the ratio between a fixed high percentile of the CO distribution and a fixed refe… ▽ More In this paper we examine the factors that shape the distribution of molecular gas surface densities on the 150 pc scale across 67 morphologically diverse star-forming galaxies in the PHANGS-ALMA CO (2-1) survey. Dividing each galaxy into radial bins, we measure molecular gas surface density contrasts, defined here as the ratio between a fixed high percentile of the CO distribution and a fixed reference level in each bin. This reference level captures the level of the faint CO floor that extends between bright filamentary features, while the intensity level of the higher percentile probes the structures visually associated with bright, dense ISM features like spiral arms, bars, and filaments. We compare these contrasts to matched percentile-based measurements of the 3.6 $μ$m emission measured using Spitzer/IRAC imaging, which trace the underlying stellar mass density. We find that the logarithms of CO contrasts on 150 pc scales are 3-4 times larger than, and positively correlated with, the logarithms of 3.6 $μ$m contrasts probing smooth non-axisymmetric stellar bar and spiral structures. The correlation appears steeper than linear, consistent with the compression of gas as it flows supersonically in response to large-scale stellar structures, even in the presence of weak or flocculent spiral arms. Stellar dynamical features appear to play an important role in setting the cloud-scale gas density in our galaxies, with gas self-gravity perhaps playing a weaker role in setting the 150 pc-scale distribution of gas densities. △ Less

Submitted 24 March, 2021; originally announced March 2021.

Comments: Accepted for publication in ApJ, 16 pages, 7 figures

New Constraints on the $^{12}$CO(2-1)/(1-0) Line Ratio Across Nearby Disc Galaxies

Authors: J. S. den Brok, D. Chatzigiannakis, F. Bigiel, J. Puschnig, A. T. Barnes, A. K. Leroy, M. J. Jiménez-Donaire, A. Usero, E. Schinnerer, E. Rosolowsky, C. M. Faesi, K. Grasha, A. Hughes, J. M. D. Kruijssen, D. Liu, L. Neumann, J. Pety, M. Querejeta, T. Saito, A. Schruba, S. Stuber

Abstract: Both the CO(2-1) and CO(1-0) lines are used to trace the mass of molecular gas in galaxies. Translating the molecular gas mass estimates between studies using different lines requires a good understanding of the behaviour of the CO(2-1)-to-CO(1-0) ratio, $R_{21}$. We compare new, high quality CO(1-0) data from the IRAM 30-m EMPIRE survey to the latest available CO(2-1) maps from HERACLES, PHANGS-A… ▽ More Both the CO(2-1) and CO(1-0) lines are used to trace the mass of molecular gas in galaxies. Translating the molecular gas mass estimates between studies using different lines requires a good understanding of the behaviour of the CO(2-1)-to-CO(1-0) ratio, $R_{21}$. We compare new, high quality CO(1-0) data from the IRAM 30-m EMPIRE survey to the latest available CO(2-1) maps from HERACLES, PHANGS-ALMA, and a new IRAM 30-m M51 Large Program. This allows us to measure $R_{21}$ across the full star-forming disc of nine nearby, massive, star-forming spiral galaxies at 27" (${\sim} 1{-}2$ kpc) resolution. We find an average $R_{21} = 0.64\pm0.09$ when we take the luminosity-weighted mean of all individual galaxies. This result is consistent with the mean ratio for disc galaxies that we derive from single-pointing measurements in the literature, $R_{\rm 21, lit}~=~0.59^{+0.18}_{-0.09}$. The ratio shows weak radial variations compared to the point-to-point scatter in the data. In six out of nine targets the central enhancement in $R_{21}$ with respect to the galaxy-wide mean is of order $\sim 10{-}20\%$. We estimate an azimuthal scatter of $\sim$20% in $R_{21}$ at fixed galactocentric radius but this measurement is limited by our comparatively coarse resolution of 1.5 kpc. We find mild correlations between $R_{21}$ and CO brightness temperature, IR intensity, 70-to-160$ μ$m ratio, and IR-to-CO ratio. All correlations indicate that $R_{21}$ increases with gas surface density, star formation rate surface density, and the interstellar radiation field. △ Less

Submitted 18 March, 2021; originally announced March 2021.

Comments: 27 pages, 18 figures, 5 tables; accepted for publication in MNRAS

Applying the Tremaine-Weinberg Method to Nearby Galaxies: Stellar Mass-Based Pattern Speeds, and Comparisons with ISM Kinematics

Authors: Thomas G. Williams, Eva Schinnerer, Eric Emsellem, Sharon Meidt, Miguel Querejeta, Francesco Belfiore, Ivana Bešlić, Frank Bigiel, Mélanie Chevance, Daniel A. Dale, Simon C. O. Glover, Kathryn Grasha, Ralf S. Klessen, J. M. Diederik Kruijssen, Adam K. Leroy, Hsi-An Pan, Jérôme Pety, Ismael Pessa, Erik Rosolowsky, Toshiki Saito, Francesco Santoro, Andreas Schruba, Mattia C. Sormani, Jiayi Sun, Elizabeth J. Watkins

Abstract: We apply the Tremaine-Weinberg method to 19 nearby galaxies using stellar mass surface densities and velocities derived from the PHANGS-MUSE survey, to calculate (primarily bar) pattern speeds ($Ω_{\rm P}$). After quality checks, we find that around half (10) of these stellar mass-based measurements are reliable. For those galaxies, we find good agreement between our results and previously publish… ▽ More We apply the Tremaine-Weinberg method to 19 nearby galaxies using stellar mass surface densities and velocities derived from the PHANGS-MUSE survey, to calculate (primarily bar) pattern speeds ($Ω_{\rm P}$). After quality checks, we find that around half (10) of these stellar mass-based measurements are reliable. For those galaxies, we find good agreement between our results and previously published pattern speeds, and use rotation curves to calculate major resonance locations (co-rotation radii and Lindblad resonances). We also compare these stellar-mass derived pattern speeds with H$α$ (from MUSE) and CO($J=2{-}1$) emission from the PHANGS-ALMA survey. We find that in the case of these clumpy ISM tracers, this method erroneously gives a signal that is simply the angular frequency at a representative radius set by the distribution of these clumps ($Ω_{\rm clump}$), and that this $Ω_{\rm clump}$ is significantly different to $Ω_{\rm P}$ ($\sim$20% in the case of H$α$, and $\sim$50% in the case of CO). Thus, we conclude that it is inadvisable to use "pattern speeds" derived from ISM kinematics. Finally, we compare our derived pattern speeds and co-rotation radii, along with bar properties, to the global parameters of these galaxies. Consistent with previous studies, we find that galaxies with a later Hubble type have a larger ratio of co-rotation radius to bar length, more molecular-gas rich galaxies have higher $Ω_{\rm P}$, and more bulge-dominated galaxies have lower $Ω_{\rm P}$. Unlike earlier works, however, there are no clear trends between the bar strength and $Ω_{\rm P}$, nor between the total stellar mass surface density and the pattern speed. △ Less

Submitted 15 March, 2021; v1 submitted 1 February, 2021; originally announced February 2021.

Comments: 22 pages, 13 figures, plus an additional 24 pages and 19 figures in appendices. Accepted for publication in AJ. Include DOI, as well as link to GitHub page containing .fits machine-readable version of Table 2

Giant Molecular Cloud Catalogues for PHANGS-ALMA: Methods and Initial Results

Authors: Erik Rosolowsky, Annie Hughes, Adam K. Leroy, Jiayi Sun, Miguel Querejeta, Andreas Schruba, Antonio Usero, Cinthya N. Herrera, Daizhong Liu, Jérôme Pety, Toshiki Saito, Ivana Bešlić, Frank Bigiel, Guillermo Blanc, Mélanie Chevance, Daniel A. Dale, Sinan Deger, Christopher M. Faesi, Simon C. O. Glover, Jonathan D. Henshaw, Ralf S. Klessen, J. M. Diederik Kruijssen, Kirsten Larson, Janice Lee, Sharon Meidt , et al. (4 additional authors not shown)

Abstract: We present improved methods for segmenting CO emission from galaxies into individual molecular clouds, providing an update to the CPROPS algorithms presented by Rosolowsky & Leroy (2006; arXiv:astro-ph/0601706 ). The new code enables both homogenization of the noise and spatial resolution among data, which allows for rigorous comparative analysis. The code also models the completeness of the data… ▽ More We present improved methods for segmenting CO emission from galaxies into individual molecular clouds, providing an update to the CPROPS algorithms presented by Rosolowsky & Leroy (2006; arXiv:astro-ph/0601706 ). The new code enables both homogenization of the noise and spatial resolution among data, which allows for rigorous comparative analysis. The code also models the completeness of the data via false source injection and includes an updated segmentation approach to better deal with blended emission. These improved algorithms are implemented in a publicly available python package, PYCPROPS. We apply these methods to ten of the nearest galaxies in the PHANGS-ALMA survey, cataloguing CO emission at a common 90 pc resolution and a matched noise level. We measure the properties of 4986 individual clouds identified in these targets. We investigate the scaling relations among cloud properties and the cloud mass distributions in each galaxy. The physical properties of clouds vary among galaxies, both as a function of galactocentric radius and as a function of dynamical environment. Overall, the clouds in our target galaxies are well-described by approximate energy equipartition, although clouds in stellar bars and galaxy centres show elevated line widths and virial parameters. The mass distribution of clouds in spiral arms has a typical mass scale that is 2.5x larger than interarm clouds and spiral arms clouds show slightly lower median virial parameters compared to interarm clouds (1.2 versus 1.4). △ Less

Submitted 12 January, 2021; originally announced January 2021.

Comments: Accepted to MNRAS. Catalogues, code, and atlas available at http://phangs.org/data

The PHANGS-HST Survey: Physics at High Angular resolution in Nearby GalaxieS with the Hubble Space Telescope

Authors: Janice C. Lee, Bradley C. Whitmore, David A. Thilker, Sinan Deger, Kirsten L. Larson, Leonardo Ubeda, Gagandeep S. Anand, Mederic Boquien, Rupali Chandar, Daniel A. Dale, Eric Emsellem, Adam K. Leroy, Erik Rosolowsky, Eva Schinnerer, Judy Schmidt, Jordan Turner, Schuyler Van Dyk, Richard L. White, Ashley T. Barnes, Francesco Belfiore, Frank Bigiel, Guillermo A. Blanc, Yixian Cao, Melanie Chevance, Enrico Congiu , et al. (29 additional authors not shown)

Abstract: The PHANGS program is building the first dataset to enable the multi-phase, multi-scale study of star formation across the nearby spiral galaxy population. This effort is enabled by large survey programs with ALMA, VLT/MUSE, and HST, with which we have obtained CO(2-1) imaging, optical spectroscopic mapping, and high resolution UV-optical imaging, respectively. Here, we present PHANGS-HST, which h… ▽ More The PHANGS program is building the first dataset to enable the multi-phase, multi-scale study of star formation across the nearby spiral galaxy population. This effort is enabled by large survey programs with ALMA, VLT/MUSE, and HST, with which we have obtained CO(2-1) imaging, optical spectroscopic mapping, and high resolution UV-optical imaging, respectively. Here, we present PHANGS-HST, which has obtained five band NUV-U-B-V-I imaging of the disks of 38 spiral galaxies at distances of 4-23 Mpc, and parallel V and I band imaging of their halos, to provide a census of tens of thousands of compact star clusters, and multi-scale stellar associations. The combination of HST, ALMA, and VLT/MUSE observations will yield an unprecedented joint catalog of the observed and physical properties of ~100,000 star clusters, associations, HII regions, and molecular clouds. With these basic units of star formation, PHANGS will systematically chart the evolutionary cycling between gas and stars, across a diversity of galactic environments found in nearby galaxies. We discuss the design of the PHANGS-HST survey, and provide an overview of the HST data processing pipeline and first results, highlighting new methods for selecting star cluster candidates, morphological classification of candidates with convolutional neural networks, and identification of stellar associations over a range of physical scales with a watershed algorithm. We describe the cross-observatory imaging, catalogs, and software products to be released. These high-level science products will seed a broad range of investigations, in particular, the study of embedded stellar populations and dust with JWST, for which a PHANGS Cycle 1 Treasury program to obtain eight band 2-21 $μ$m imaging has been approved. △ Less

Submitted 19 October, 2021; v1 submitted 8 January, 2021; originally announced January 2021.

Comments: ApJS in press. Accepted draft includes moderate revision and updates to the originally posted submitted version. 27 pages, 11 figures, 1 table. Data products released at https://archive.stsci.edu/hlsp/phangs-hst/

PHANGS-HST: Star Cluster Spectral Energy Distribution Fitting with CIGALE

Authors: Jordan A. Turner, Daniel A. Dale, Janice C. Lee, Mederic Boquien, Rupali Chandar, Sinan Deger, Kirsten L. Larson, Angus Mok, David A. Thilker, Leonardo Ubeda, Bradley C. Whitmore, Francesco Belfiore, Frank Bigiel, Guillermo A. Blanc, Eric Emsellem, Kathryn Grasha, Brent Groves, Ralf S. Klessen, Kathryn Kreckel, J. M. Diederik Kruijssen, Adam K. Leroy, Erik Rosolowsky, Patricia Sanchez-Blazquez, Eva Schinnerer, Andreas Schruba , et al. (2 additional authors not shown)

Abstract: The sensitivity and angular resolution of photometric surveys executed by the Hubble Space Telescope (HST) enable studies of individual star clusters in galaxies out to a few tens of megaparsecs. The fitting of spectral energy distributions (SEDs) of star clusters is essential for measuring their physical properties and studying their evolution. We report on the use of the publicly available Code… ▽ More The sensitivity and angular resolution of photometric surveys executed by the Hubble Space Telescope (HST) enable studies of individual star clusters in galaxies out to a few tens of megaparsecs. The fitting of spectral energy distributions (SEDs) of star clusters is essential for measuring their physical properties and studying their evolution. We report on the use of the publicly available Code Investigating GALaxy Emission (CIGALE) SED fitting package to derive ages, stellar masses, and reddenings for star clusters identified in the Physics at High Angular resolution in Nearby GalaxieS-HST (PHANGS-HST) survey. Using samples of star clusters in the galaxy NGC 3351, we present results of benchmark analyses performed to validate the code and a comparison to SED fitting results from the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). We consider procedures for the PHANGS-HST SED fitting pipeline, e.g., the choice of single stellar population models, the treatment of nebular emission and dust, and the use of fluxes versus magnitudes for the SED fitting. We report on the properties of clusters in NGC 3351 and find, on average, the clusters residing in the inner star-forming ring of NGC 3351 are young ($< 10$ Myr) and massive ($10^{5} M_{\odot}$) while clusters in the stellar bulge are significantly older. Cluster mass function fits yield $β$ values around -2, consistent with prior results with a tendency to be shallower at the youngest ages. Finally, we explore a Bayesian analysis with additional physically-motivated priors for the distribution of ages and masses and analyze the resulting cluster distributions. △ Less

Submitted 6 January, 2021; originally announced January 2021.

Comments: Accepted for publication in MNRAS (2021-01-06), 22 pages, 18 figures, 2 tables

The KMOS$^{\rm 3D}$ Survey: Investigating the Origin of the Elevated Electron Densities in Star-Forming Galaxies at $1\lesssim{z}\lesssim{3}$

Authors: Rebecca L. Davies, N. M. Förster Schreiber, R. Genzel, T. T. Shimizu, R. I. Davies, A. Schruba, L. J. Tacconi, H. Übler, E. Wisnioski, S. Wuyts, M. Fossati, R. Herrera-Camus, D. Lutz, J. T. Mendel, T. Naab, S. H. Price, A. Renzini, D. Wilman, A. Beifiori, S. Belli, A. Burkert, J. Chan, A. Contursi, M. Fabricius, M. M. Lee , et al. (2 additional authors not shown)

Abstract: We investigate what drives the redshift evolution of the typical electron density ($n_e$) in star-forming galaxies, using a sample of 140 galaxies drawn primarily from KMOS$^{\rm 3D}$ ($0.6\lesssim{z}\lesssim{2.6}$) and 471 galaxies from SAMI ($z<0.113$). We select galaxies that do not show evidence of AGN activity or outflows, to constrain the average conditions within H II regions. Measurements… ▽ More We investigate what drives the redshift evolution of the typical electron density ($n_e$) in star-forming galaxies, using a sample of 140 galaxies drawn primarily from KMOS$^{\rm 3D}$ ($0.6\lesssim{z}\lesssim{2.6}$) and 471 galaxies from SAMI ($z<0.113$). We select galaxies that do not show evidence of AGN activity or outflows, to constrain the average conditions within H II regions. Measurements of the [SII]$λ$6716/[SII]$λ$6731 ratio in four redshift bins indicate that the local $n_e$ in the line-emitting material decreases from 187$^{+140}_{-132}$ cm$^{-3}$ at $z\sim$ 2.2 to 32$^{+4}_{-9}$ cm$^{-3}$ at $z\sim$ 0; consistent with previous results. We use the H$α$ luminosity to estimate the root-mean-square (rms) $n_e$ averaged over the volumes of star-forming disks at each redshift. The local and volume-averaged $n_e$ evolve at similar rates, hinting that the volume filling factor of the line-emitting gas may be approximately constant across $0\lesssim{z}\lesssim{2.6}$. The KMOS$^{\rm 3D}$ and SAMI galaxies follow a roughly monotonic trend between $n_e$ and star formation rate, but the KMOS$^{\rm 3D}$ galaxies have systematically higher $n_e$ than the SAMI galaxies at fixed offset from the star-forming main sequence, suggesting a link between the $n_e$ evolution and the evolving main sequence normalization. We quantitatively test potential drivers of the density evolution and find that $n_e$(rms) $\simeq{n_{H_2}}$, suggesting that the elevated $n_e$ in high-$z$ H II regions could plausibly be the direct result of higher densities in the parent molecular clouds. There is also tentative evidence that $n_e$ could be influenced by the balance between stellar feedback, which drives the expansion of H II regions, and the ambient pressure, which resists their expansion. △ Less

Submitted 18 December, 2020; originally announced December 2020.

Comments: Main text 24 pages, 8 figures. Accepted for publication in ApJ

Distances to PHANGS Galaxies: New Tip of the Red Giant Branch Measurements and Adopted Distances

Authors: Gagandeep S. Anand, Janice C. Lee, Schuyler D. Van Dyk, Adam K. Leroy, Erik Rosolowsky, Eva Schinnerer, Kirsten Larson, Ehsan Kourkchi, Kathryn Kreckel, Fabian Scheuermann, Luca Rizzi, David Thilker, R. Brent Tully, Frank Bigiel, Guillermo A. Blanc, Médéric Boquien, Rupali Chandar, Daniel Dale, Eric Emsellem, Sinan Deger, Simon C. O. Glover, Kathryn Grasha, Brent Groves, Ralf S. Klessen, J. M. Diederik Kruijssen , et al. (7 additional authors not shown)

Abstract: PHANGS-HST is an ultraviolet-optical imaging survey of 38 spiral galaxies within ~20 Mpc. Combined with the PHANGS-ALMA, PHANGS-MUSE surveys and other multiwavelength data, the dataset will provide an unprecedented look into the connections between young stars, HII regions, and cold molecular gas in these nearby star-forming galaxies. Accurate distances are needed to transform measured observables… ▽ More PHANGS-HST is an ultraviolet-optical imaging survey of 38 spiral galaxies within ~20 Mpc. Combined with the PHANGS-ALMA, PHANGS-MUSE surveys and other multiwavelength data, the dataset will provide an unprecedented look into the connections between young stars, HII regions, and cold molecular gas in these nearby star-forming galaxies. Accurate distances are needed to transform measured observables into physical parameters (e.g., brightness to luminosity, angular to physical sizes of molecular clouds, star clusters and associations). PHANGS-HST has obtained parallel ACS imaging of the galaxy halos in the F606W and F814W bands. Where possible, we use these parallel fields to derive tip of the red giant branch (TRGB) distances to these galaxies. In this paper, we present TRGB distances for 11 galaxies from ~4 to ~15 Mpc, based on the first year of PHANGS-HST observations. Five of these represent the first published TRGB distance measurements (IC 5332, NGC 2835, NGC 4298, NGC 4321, and NGC 4328), and eight of which are the best available distances to these targets. We also provide a compilation of distances for the 118 galaxies in the full PHANGS sample, which have been adopted for the first PHANGS-ALMA public data release. △ Less

Submitted 1 December, 2020; originally announced December 2020.

Comments: 20 pages, 8 figures, accepted to MNRAS

On the duration of the embedded phase of star formation

Authors: Jaeyeon Kim, Mélanie Chevance, J. M. Diederik Kruijssen, Andreas Schruba, Karin Sandstrom, Ashley T. Barnes, Frank Bigiel, Guillermo A. Blanc, Yixian Cao, Daniel A. Dale, Christopher M. Faesi, Simon C. O. Glover, Kathryn Grasha, Brent Groves, Cinthya Herrera, Ralf S. Klessen, Kathryn Kreckel, Janice C. Lee, Adam K. Leroy, Jérôme Pety, Miguel Querejeta, Eva Schinnerer, Jiayi Sun, Antonio Usero, Jacob L. Ward , et al. (1 additional authors not shown)

Abstract: Feedback from massive stars plays a key role in molecular cloud evolution. After the onset of star formation, the young stellar population is exposed by photoionization, winds, supernovae, and radiation pressure from massive stars. Recent observations of nearby galaxies have provided the evolutionary timeline between molecular clouds and exposed young stars, but the duration of the embedded phase… ▽ More Feedback from massive stars plays a key role in molecular cloud evolution. After the onset of star formation, the young stellar population is exposed by photoionization, winds, supernovae, and radiation pressure from massive stars. Recent observations of nearby galaxies have provided the evolutionary timeline between molecular clouds and exposed young stars, but the duration of the embedded phase of massive star formation is still ill-constrained. We measure how long massive stellar populations remain embedded within their natal cloud, by applying a statistical method to six nearby galaxies at 20-100 pc resolution, using CO, Spitzer 24$\rm\,μm$, and H$α$ emission as tracers of molecular clouds, embedded star formation, and exposed star formation, respectively. We find that the embedded phase (with CO and 24$\rm\,μm$ emission) lasts for $2{-}7$ Myr and constitutes $17{-}47\%$ of the cloud lifetime. During approximately the first half of this phase, the region is invisible in H$α$, making it heavily obscured. For the second half of this phase, the region also emits in H$α$ and is partially exposed. Once the cloud has been dispersed by feedback, 24$\rm\,μm$ emission no longer traces ongoing star formation, but remains detectable for another $2{-}9$ Myr through the emission from ambient CO-dark gas, tracing star formation that recently ended. The short duration of massive star formation suggests that pre-supernova feedback (photoionization and winds) is important in disrupting molecular clouds. The measured timescales do not show significant correlations with environmental properties (e.g. metallicity). Future JWST observations will enable these measurements routinely across the nearby galaxy population. △ Less

Submitted 26 March, 2021; v1 submitted 30 November, 2020; originally announced December 2020.

Comments: 23 pages, 9 figures, 4 tables; accepted for publication in MNRAS (March 19, 2021)

Resolving the Dust-to-Metals Ratio and CO-to-H$_2$ Conversion Factor in the Nearby Universe

Authors: I-Da Chiang, Karin M. Sandstrom, Jérémy Chastenet, Cinthya N. Herrera, Eric W. Koch, Kathryn Kreckel, Adam K. Leroy, Jérôme Pety, Andreas Schruba, Dyas Utomo, Thomas Williams

Abstract: We investigate the relationship between the dust-to-metals ratio (D/M) and the local interstellar medium environment at ~2 kpc resolution in five nearby galaxies: IC342, M31, M33, M101, and NGC628. A modified blackbody model with a broken power-law emissivity is used to model the dust emission from 100 to 500 um observed by Herschel. We utilize the metallicity gradient derived from auroral line me… ▽ More We investigate the relationship between the dust-to-metals ratio (D/M) and the local interstellar medium environment at ~2 kpc resolution in five nearby galaxies: IC342, M31, M33, M101, and NGC628. A modified blackbody model with a broken power-law emissivity is used to model the dust emission from 100 to 500 um observed by Herschel. We utilize the metallicity gradient derived from auroral line measurements in HII regions whenever possible. Both archival and new CO rotational line and HI 21 cm maps are adopted to calculate gas surface density, including new wide field CO and HI maps for IC342 from IRAM and the VLA, respectively. We experiment with several prescriptions of CO-to-H$_2$ conversion factor, and compare the resulting D/M-metallicity and D/M-density correlations, both of which are expected to be non-negative from depletion studies. The D/M is sensitive to the choice of the conversion factor. The conversion factor prescriptions based on metallicity only yield too much molecular gas in the center of IC342 to obtain the expected correlations. Among the prescriptions tested, the one that yields the expected correlations depends on both metallicity and surface density. The 1-$σ$ range of the derived D/M spans 0.40-0.58. Compared to chemical evolution models, our measurements suggest that the dust growth time scale is much shorter than the dust destruction time scale. The measured D/M is consistent with D/M in galaxy-integrated studies derived from infrared dust emission. Meanwhile, the measured D/M is systematically higher than the D/M derived from absorption, which likely indicates a systematic offset between the two methods. △ Less

Submitted 22 January, 2021; v1 submitted 20 November, 2020; originally announced November 2020.

Comments: 22 pages, 8 figures, 4 tables; accepted by ApJ

Pre-supernova feedback mechanisms drive the destruction of molecular clouds in nearby star-forming disc galaxies

Authors: Mélanie Chevance, J. M. Diederik Kruijssen, Mark R. Krumholz, Brent Groves, Benjamin W. Keller, Annie Hughes, Simon C. O. Glover, Jonathan D. Henshaw, Cinthya N. Herrera, Jenny J. Kim, Adam K. Leroy, Jérôme Pety, Alessandro Razza, Erik Rosolowsky, Eva Schinnerer, Andreas Schruba, Ashley T. Barnes, Frank Bigiel, Guillermo A. Blanc, Eric Emsellem, Christopher M. Faesi, Kathryn Grasha, Ralf S. Klessen, Kathryn Kreckel, Daizhong Liu , et al. (6 additional authors not shown)

Abstract: It is a major open question which physical processes stop the accretion of gas onto giant molecular clouds (GMCs) and limit the efficiency at which gas is converted into stars within these GMCs. While feedback from supernova explosions has been the popular feedback mechanism included in simulations of galaxy formation and evolution, `early' feedback mechanisms such as stellar winds, photoionisatio… ▽ More It is a major open question which physical processes stop the accretion of gas onto giant molecular clouds (GMCs) and limit the efficiency at which gas is converted into stars within these GMCs. While feedback from supernova explosions has been the popular feedback mechanism included in simulations of galaxy formation and evolution, `early' feedback mechanisms such as stellar winds, photoionisation and radiation pressure are expected to play an important role in dispersing the gas after the onset of star formation. These feedback processes typically take place on small scales ($\sim 10-100$ pc) and their effects have therefore been difficult to constrain in environments other than the Milky Way. We apply a novel statistical method to $\sim 1$" resolution maps of CO and Ha emission across a sample of nine nearby disc galaxies, in order to measure the time over which GMCs are dispersed by feedback from young, high-mass stars, as a function of the galactic environment. We find that GMCs are typically dispersed within $\sim$ 3 Myr after the emergence of unembedded high-mass stars, showing no significant trend with galactocentric radius. Comparison with analytical predictions demonstrates that, independently of the environment, early feedback mechanisms (particularly photoionisation and stellar winds) play a crucial role in dispersing GMCs and limiting their star formation efficiency in nearby galaxies. Finally, we show that the efficiency at which the energy injected by these early feedback mechanisms couples with the parent GMC is relatively low (a few tens of per cent), such that the vast majority of momentum and energy emitted by the young stellar populations escapes the parent GMC. △ Less

Submitted 26 October, 2020; originally announced October 2020.

Comments: 17 pages, 8 figures, 3 tables; submitted to MNRAS (October 23, 2020)

Measuring the mixing scale of the ISM within nearby spiral galaxies

Authors: Kathryn Kreckel, I-Ting Ho, Guillermo A. Blanc, Simon C. O. Glover, Brent Groves, Erik Rosolowsky, Frank Bigiel, Mederic Boquien, Melanie Chevance, Daniel A. Dale, Sinan Deger, Eric Emsellem, Kathryn Grasha, Jenny J. Kim, Ralf S. Klessen, J. M. Diederik Kruijssen, Janice C. Lee, Adam K. Leroy, Daizhong Liu, Rebecca McElroy, Sharon E. Meidt, Ismael Pessa, Patricia Sanchez-Blazquez, Karin Sandstrom, Francesco Santoro , et al. (6 additional authors not shown)

Abstract: The spatial distribution of metals reflects, and can be used to constrain, the processes of chemical enrichment and mixing. Using PHANGS-MUSE optical integral field spectroscopy, we measure the gas phase oxygen abundances (metallicities) across 7,138 HII regions in a sample of eight nearby disc galaxies. In Paper I (Kreckel et al. 2019) we measure and report linear radial gradients in the metallic… ▽ More The spatial distribution of metals reflects, and can be used to constrain, the processes of chemical enrichment and mixing. Using PHANGS-MUSE optical integral field spectroscopy, we measure the gas phase oxygen abundances (metallicities) across 7,138 HII regions in a sample of eight nearby disc galaxies. In Paper I (Kreckel et al. 2019) we measure and report linear radial gradients in the metallicities of each galaxy, and qualitatively searched for azimuthal abundance variations. Here, we examine the two-dimensional variation in abundances once the radial gradient is subtracted, Delta(O/H), in order to quantify the homogeneity of the metal distribution and to measure the mixing scale over which HII region metallicities are correlated. We observe low (0.03--0.05 dex) scatter in Delta(O/H) globally in all galaxies, with significantly lower (0.02--0.03 dex) scatter on small (<600 pc) spatial scales. This is consistent with the measurement uncertainties, and implies the two-dimensional metallicity distribution is highly correlated on scales of <600 pc. We compute the two point correlation function for metals in the disc in order to quantify the scale lengths associated with the observed homogeneity. This mixing scale is observed to correlate better with the local gas velocity dispersion (of both cold and ionized gas) than with the star formation rate. Selecting only HII regions with enhanced abundances relative to a linear radial gradient, we do not observe increased homogeneity on small scales. This suggests that the observed homogeneity is driven by the mixing introducing material from large scales rather than by pollution from recent and on-going star formation. △ Less

Submitted 4 September, 2020; originally announced September 2020.

Comments: 17 pages, 14 figures. Accepted for publication in MNRAS