Showing 1–50 of 294 results for author: Kruijssen, J M D

Death of the Immortal Molecular Cloud: Resolution Dependence of the Gas-Star Formation Relation Rules out Decoupling by Stellar Drift

Authors: J. M. Diederik Kruijssen, Mélanie Chevance, Steven N. Longmore, Adam Ginsburg, Lise Ramambason, Andrea Romanelli

Abstract: Recent observations have demonstrated that giant molecular clouds (GMCs) are short-lived entities, surviving for the order of a dynamical time before turning a few percent of their mass into stars and dispersing, leaving behind an isolated young stellar population. The key question has been whether this GMC dispersal actually marks a point of GMC destruction by stellar feedback from the new-born s… ▽ More Recent observations have demonstrated that giant molecular clouds (GMCs) are short-lived entities, surviving for the order of a dynamical time before turning a few percent of their mass into stars and dispersing, leaving behind an isolated young stellar population. The key question has been whether this GMC dispersal actually marks a point of GMC destruction by stellar feedback from the new-born stars, or if GMCs might be `immortal' and only dynamically decouple from their nascent stars due to stellar drift. We address this question in six nearby galaxies, by quantifying how the gas-star formation relation depends on the spatial scale for scales between the GMC diameter and the GMC separation length, i.e. the scales where an excess of GMCs would be expected to be found in the stellar drift scenario. Our analysis reveals a consistent dearth of GMCs near young stellar populations regardless of the spatial scale, discounting the notion of `immortal' GMCs that decouple from their nascent stars through stellar drift. Instead, our findings demonstrate that stellar feedback destroys most GMCs at the end of their lifecycle. Employing a variety of statistical techniques to test both hypotheses, we find that the probability that stellar feedback concludes the GMC lifecycle is about 2,000 times higher than the probability that stellar drift separates GMCs and young stellar regions. This observation strengthens the emerging picture that galaxies consist of dynamic building blocks undergoing vigorous, feedback-driven lifecycles that collectively regulate star formation and drive the baryon cycle within galaxies. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 10 pages, 4 figures, 3 tables; submitted to The Open Journal of Astrophysics (April 24, 2024); Figure 1 explains the experiment, Figure 2 shows the measurements, Table 3 shows the main result

Origin of the correlation between stellar kinematics and globular cluster system richness in ultra-diffuse galaxies

Authors: Joel Pfeffer, Steven R. Janssens, Maria Luisa Buzzo, Jonah S. Gannon, Nate Bastian, Kenji Bekki, Jean P. Brodie, Warrick J. Couch, Robert A. Crain, Duncan A. Forbes, J. M. Diederik Kruijssen, Aaron J. Romanowsky

Abstract: Observational surveys have found that the dynamical masses of ultra-diffuse galaxies (UDGs) correlate with the richness of their globular cluster (GC) system. This could be explained if GC-rich galaxies formed in more massive dark matter haloes. We use simulations of galaxies and their GC systems from the E-MOSAICS project to test whether the simulations reproduce such a trend. We find that GC-ric… ▽ More Observational surveys have found that the dynamical masses of ultra-diffuse galaxies (UDGs) correlate with the richness of their globular cluster (GC) system. This could be explained if GC-rich galaxies formed in more massive dark matter haloes. We use simulations of galaxies and their GC systems from the E-MOSAICS project to test whether the simulations reproduce such a trend. We find that GC-rich simulated galaxies in galaxy groups have enclosed masses that are consistent with the dynamical masses of observed GC-rich UDGs. However, simulated GC-poor galaxies in galaxy groups have higher enclosed masses than those observed. We argue that GC-poor UDGs with low stellar velocity dispersions are discs observed nearly face on, such that their true mass is underestimated by observations. Using the simulations, we show that galactic star-formation conditions resulting in dispersion-supported stellar systems also leads to efficient GC formation. Conversely, conditions leading to rotationally-supported discs leads to inefficient GC formation. This result may explain why early-type galaxies typically have richer GC systems than late-type galaxies. This is also supported by comparisons of stellar axis ratios and GC specific frequencies in observed dwarf galaxy samples, which show GC-rich systems are consistent with being spheroidal, while GC-poor systems are consistent with being discs. Therefore, particularly for GC-poor galaxies, rotation should be included in dynamical mass measurements from stellar dynamics. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: 16 pages, 8 figures. Accepted for publication in MNRAS

The properties and kinematics of HCN emission across the closest starburst galaxy NGC 253 observed with ALMA

Authors: Ivana Beslic, Ashley T. Barnes, Frank Bigiel, Maria Jesus Jimenez-Donaire, Antonio Usero, Jonathan D. Henshaw, Christopher Faesi, Adam K. Leroy, Erik Rosolowsky, Jakob S. den Brok, Melanie Chevance, Cosima Eibensteiner, Kathryn Grasha, Ralf S. Klessen, J. M. Diedrerik Kruijssen, Daizhong Liu, Sharon Meidt, Justus Neumann, Lukas Neumann, Hsi-An Pan, Johannes Puschnig, Miguel Querejeta, Eva Schinnerer, Thomas G. Williams

Abstract: Studying molecular gas in nearby galaxies using hydrogen cyanide (HCN) as a tracer for higher densities than CO emission still poses a significant challenge. Even though several galaxies have HCN maps on a few kpc scales, higher-resolution maps are still required. Our goal is to examine the contrast in intensity between two tracers that probe different density regimes - HCN(1-0)/CO(2-1) ratio - an… ▽ More Studying molecular gas in nearby galaxies using hydrogen cyanide (HCN) as a tracer for higher densities than CO emission still poses a significant challenge. Even though several galaxies have HCN maps on a few kpc scales, higher-resolution maps are still required. Our goal is to examine the contrast in intensity between two tracers that probe different density regimes - HCN(1-0)/CO(2-1) ratio - and their kinematics across NGC 253. By utilizing the advanced capabilities of the Atacama Large Millimeter/submillimeter Array (ALMA), we can map these features at high resolution across a large field of view and uncover the nature of such dense gas in extragalactic systems. We present new ALMA Atacama Compact Array and Total Power (ACA+TP) observations of the HCN emission across NGC 253, covering the inner 8.6' of the galaxy disk at 300 pc scales. We analyze the integrated intensity and mean velocity of HCN and CO along each line of sight and use SCOUSE software to perform spectral decomposition, which considers each velocity component separately. Molecular gas traced by HCN piles up in a ring-like structure at a radius of 2 kpc. The HCN emission is enhanced by 2 orders of magnitude in the central 2 kpc regions, beyond which its intensity decreases with increasing galactocentric distance. The number of components in the HCN spectra shows a robust environmental dependence, with multiple velocity features across the center and bar. We have identified an increase in the HCN/CO ratio in these regions, corresponding to a velocity component likely associated with a molecular outflow. We have also discovered that the ratio between the total infrared luminosity and dense gas mass, which indicates the star formation efficiency of dense gas, is anti-correlated with the molecular gas surface density up to approximately 200 Msul/pc^2. In contrast, beyond this point, the ratio starts to increase. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: Accepted for publication to Astronomy and Astrophysics

Magnetic field morphology and evolution in the Central Molecular Zone and its effect on gas dynamics

Authors: R. G. Tress, M. C. Sormani, P. Girichidis, S. C. O. Glover, R. S. Klessen, R. J. Smith, E. Sobacchi, L. Armillotta, A. T. Barnes, C. Battersby, K. R. J. Bogue, N. Brucy, L. Colzi, C. Federrath, P. García, A. Ginsburg, J. Göller, H P. Hatchfield, C. Henkel, P. Hennebelle, J. D. Henshaw, M. Hirschmann, Y. Hu, J. Kauffmann, J. M. D. Kruijssen , et al. (12 additional authors not shown)

Abstract: The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We f… ▽ More The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We find that: (1) The magnetic field is conveniently decomposed into a regular time-averaged component and an irregular turbulent component. The regular component aligns well with the velocity vectors of the gas everywhere, including within the bar lanes. (2) The field geometry transitions from parallel to the Galactic plane near $z=0$ to poloidal away from the plane. (3) The magneto-rotational instability (MRI) causes an in-plane inflow of matter from the CMZ gas ring towards the central few parsecs of $0.01-0.1$ M$_\odot$ yr$^{-1}$ that is absent in the unmagnetised simulations. However, the magnetic fields have no significant effect on the larger-scale bar-driven inflow that brings the gas from the Galactic disc into the CMZ. (4) A combination of bar inflow and MRI-driven turbulence can sustain a turbulent vertical velocity dispersion of $σ_z \simeq 5$ km s$^{-1}$ on scales of $20$ pc in the CMZ ring. The MRI alone sustains a velocity dispersion of $σ_z \simeq 3$ km s$^{-1}$. Both these numbers are lower than the observed velocity dispersion of gas in the CMZ, suggesting that other processes such as stellar feedback are necessary to explain the observations. (5) Dynamo action driven by differential rotation and the MRI amplifies the magnetic fields in the CMZ ring until they saturate at a value that scales with the average local density as $B \simeq 102 (n/10^3 cm^{-3})^{0.33}$ $μ$G. Finally, we discuss the implications of our results within the observational context in the CMZ. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: submitted to A&A, comments welcome

CMZoom IV. Incipient High-Mass Star Formation Throughout the Central Molecular Zone

Authors: H Perry Hatchfield, Cara Battersby, Ashley T. Barnes, Natalie Butterfield, Adam Ginsburg, Jonathan D. Henshaw, Steven N. Longmore, Xing Lu, Brian Svoboda, Daniel Walker, Daniel Callanan, Elisabeth A. C. Mills, Luis C. Ho, Jens Kauffmann, J. M. Diederik Kruijssen, Jürgen Ott, Thushara Pillai, Qizhou Zhang

Abstract: In this work, we constrain the star-forming properties of all possible sites of incipient high-mass star formation in the Milky Way's Galactic Center. We identify dense structures using the CMZoom 1.3mm dust continuum catalog of objects with typical radii of $\sim$0.1pc, and measure their association with tracers of high-mass star formation. We incorporate compact emission at 8, 21, 24, 25, and 70… ▽ More In this work, we constrain the star-forming properties of all possible sites of incipient high-mass star formation in the Milky Way's Galactic Center. We identify dense structures using the CMZoom 1.3mm dust continuum catalog of objects with typical radii of $\sim$0.1pc, and measure their association with tracers of high-mass star formation. We incorporate compact emission at 8, 21, 24, 25, and 70um from MSX, Spitzer, Herschel, and SOFIA, catalogued young stellar objects, and water and methanol masers to characterize each source. We find an incipient star formation rate (SFR) for the CMZ of ~0.08 Msun yr^{-1} over the next few 10^5 yr. We calculate upper and lower limits on the CMZ's incipient SFR of ~0.45 Msun yr^{-1} and ~0.05 Msun yr^{-1} respectively, spanning between roughly equal to and several times greater than other estimates of CMZ's recent SFR. Despite substantial uncertainties, our results suggest the incipient SFR in the CMZ may be higher than previously estimated. We find that the prevalence of star formation tracers does not correlate with source volume density, but instead ~75% of high-mass star formation is found in regions above a column density ratio (N_{SMA}/N_{Herschel}) of ~1.5. Finally, we highlight the detection of ``atoll sources'', a reoccurring morphology of cold dust encircling evolved infrared sources, possibly representing HII regions in the process of destroying their envelopes. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: Accepted for publication in ApJ

The PHANGS-AstroSat Atlas of Nearby Star Forming Galaxies

Authors: Hamid Hassani, Erik Rosolowsky, Eric W. Koch, Joseph Postma, Joseph Nofech, Harrisen Corbould, David Thilker, Adam K. Leroy, Eva Schinnerer, Francesco Belfiore, Frank Bigiel, Mederic Boquien, Melanie Chevance, Daniel A. Dale, Oleg V. Egorov, Eric Emsellem, Simon C. O. Glover, Kathryn Grasha, Brent Groves, Kiana Henny, Jaeyeon Kim, Ralf S. Klessen, Kathryn Kreckel, J. M. Diederik Kruijssen, Janice C. Lee , et al. (7 additional authors not shown)

Abstract: We present the Physics at High Angular resolution in Nearby GalaxieS (PHANGS)-AstroSat atlas, which contains ultraviolet imaging of 31 nearby star-forming galaxies captured by the Ultraviolet Imaging Telescope (UVIT) on the AstroSat satellite. The atlas provides a homogeneous data set of far- and near-ultraviolet maps of galaxies within a distance of 22 Mpc and a median angular resolution of 1.4 a… ▽ More We present the Physics at High Angular resolution in Nearby GalaxieS (PHANGS)-AstroSat atlas, which contains ultraviolet imaging of 31 nearby star-forming galaxies captured by the Ultraviolet Imaging Telescope (UVIT) on the AstroSat satellite. The atlas provides a homogeneous data set of far- and near-ultraviolet maps of galaxies within a distance of 22 Mpc and a median angular resolution of 1.4 arcseconds (corresponding to a physical scale between 25 and 160 pc). After subtracting a uniform ultraviolet background and accounting for Milky Way extinction, we compare our estimated flux densities to GALEX observations, finding good agreement. We find candidate extended UV disks around the galaxies NGC 6744 and IC 5332. We present the first statistical measurements of the clumping of the UV emission and compare it to the clumping of molecular gas traced with ALMA. We find that bars and spiral arms exhibit the highest degree of clumping, and the molecular gas is even more clumped than the FUV emission in galaxies. We investigate the variation of the ratio of observed FUV to H$α$ in different galactic environments and kpc-sized apertures. We report that $\sim 65 \%$ varation of the $\log_{10}$(FUV/H$α$) can be described through a combination of dust attenuation with star formation history parameters. The PHANGS-AstroSat atlas enhances the multi-wavelength coverage of our sample, offering a detailed perspective on star formation. When integrated with PHANGS data sets from ALMA, VLT-MUSE, HST and JWST, it develops our comprehensive understanding of attenuation curves and dust attenuation in star-forming galaxies. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: 35 pages, 16 figures. The survey webpage is available at https://sites.google.com/view/phangs/home/data/astrosat and the data archive can be accessed at https://www.canfar.net/storage/vault/list/phangs/RELEASES/PHANGS-AstroSat/v1p0

The JWST Galactic Center Survey -- A White Paper

Authors: Rainer Schoedel, Steve Longmore, Jonny Henshaw, Adam Ginsburg, John Bally, Anja Feldmeier, Matt Hosek, Francisco Nogueras Lara, Anna Ciurlo, Mélanie Chevance, J. M. Diederik Kruijssen, Ralf Klessen, Gabriele Ponti, Pau Amaro-Seoane, Konstantina Anastasopoulou, Jay Anderson, Maria Arias, Ashley T. Barnes, Cara Battersby, Giuseppe Bono, Lucía Bravo Ferres, Aaron Bryant, Miguel Cano Gonzáalez, Santi Cassisi, Leonardo Chaves-Velasquez , et al. (85 additional authors not shown)

Abstract: The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of… ▽ More The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of the most well-studied regions in astrophysics. Due to its proximity, we can study the center of our Galaxy on scales down to a few hundred AU, a hundred times better than in similar Local Group galaxies and thousands of times better than in the nearest active galaxies. The Galactic Center (GC) is therefore of outstanding astrophysical interest. However, in spite of intense observational work over the past decades, there are still fundamental things unknown about the GC. JWST has the unique capability to provide us with the necessary, game-changing data. In this White Paper, we advocate for a JWST NIRCam survey that aims at solving central questions, that we have identified as a community: i) the 3D structure and kinematics of gas and stars; ii) ancient star formation and its relation with the overall history of the Milky Way, as well as recent star formation and its implications for the overall energetics of our galaxy's nucleus; and iii) the (non-)universality of star formation and the stellar initial mass function. We advocate for a large-area, multi-epoch, multi-wavelength NIRCam survey of the inner 100\,pc of the Galaxy in the form of a Treasury GO JWST Large Program that is open to the community. We describe how this survey will derive the physical and kinematic properties of ~10,000,000 stars, how this will solve the key unknowns and provide a valuable resource for the community with long-lasting legacy value. △ Less

Submitted 14 March, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: This White Paper will be updated when required (e.g. new authors joining, editing of content). Most recent update: 24 Oct 2023

The WISDOM of power spectra: how the galactic gravitational potential impacts a galaxy's central gas reservoir in simulations and observations

Authors: Jindra Gensior, Timothy A. Davis, Martin Bureau, J. M. Diederik Kruijssen, Michele Cappellari, Ilaria Ruffa, Thomas G. Williams

Abstract: Observations indicate that the central gas discs are smoother in early-type galaxies than their late-type counterparts, while recent simulations predict that the dynamical suppression of star formation in spheroid-dominated galaxies is preceded by the suppression of fragmentation of their interstellar media. The mass surface density power spectrum is a powerful tool to constrain the degree of stru… ▽ More Observations indicate that the central gas discs are smoother in early-type galaxies than their late-type counterparts, while recent simulations predict that the dynamical suppression of star formation in spheroid-dominated galaxies is preceded by the suppression of fragmentation of their interstellar media. The mass surface density power spectrum is a powerful tool to constrain the degree of structure within a gas reservoir. Specifically here, we focus on the power spectrum slope and aim to constrain whether the shear induced by a dominant spheroidal potential can induce sufficient turbulence to suppress fragmentation, resulting in the smooth central gas discs observed. We compute surface density power spectra for the nuclear gas reservoirs of fourteen simulated isolated galaxies and twelve galaxies observed as part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project. Both simulated and observed galaxies range from disc-dominated galaxies to spheroids, with central stellar mass surface densities, a measure of bulge dominance, varying by more than an order of magnitude. For the simulations, the power spectra steepen with increasing central stellar mass surface density, thereby clearly linking the suppression of fragmentation to the shear-driven turbulence induced by the spheroid. The WISDOM observations show a different (but potentially consistent) picture: while there is no correlation between the power spectrum slopes and the central stellar mass surface densities, the slopes scatter around a value of 2.6. This is similar to the behaviour of the slopes of the simulated galaxies with high central stellar mass surface densities, and could indicate that high shear eventually drives incompressible turbulence. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 14 pages, 7 figures + appendices (7 pages, 9 figures); accepted by MNRAS (October 09th, 2023)

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)

HI discs of L$_{\ast}$ galaxies as probes of the baryonic physics of galaxy evolution

Authors: Jindra Gensior, Robert Feldmann, Marta Reina-Campos, Sebastian Trujillo-Gomez, Lucio Mayer, Benjamin W. Keller, Andrew Wetzel, J. M. Diederik Kruijssen, Philip F. Hopkins, Jorge Moreno

Abstract: Understanding what shapes the cold gas component of galaxies, which both provides the fuel for star formation and is strongly affected by the subsequent stellar feedback, is a crucial step towards a better understanding of galaxy evolution. Here, we analyse the HI properties of a sample of 46 Milky Way halo-mass galaxies, drawn from cosmological simulations (EMP-Pathfinder and FIREbox). This set o… ▽ More Understanding what shapes the cold gas component of galaxies, which both provides the fuel for star formation and is strongly affected by the subsequent stellar feedback, is a crucial step towards a better understanding of galaxy evolution. Here, we analyse the HI properties of a sample of 46 Milky Way halo-mass galaxies, drawn from cosmological simulations (EMP-Pathfinder and FIREbox). This set of simulations comprises galaxies evolved self-consistently across cosmic time with different baryonic sub-grid physics: three different star formation models [constant star formation efficiency (SFE) with different star formation eligibility criteria, and an environmentally-dependent, turbulence-based SFE] and two different feedback prescriptions, where only one sub-sample includes early stellar feedback. We use these simulations to assess the impact of different baryonic physics on the HI content of galaxies. We find that the galaxy-wide HI properties agree with each other and with observations. However, differences appear for small-scale properties. The thin HI discs observed in the local Universe are only reproduced with a turbulence-dependent SFE and/or early stellar feedback. Furthermore, we find that the morphology of HI discs is particularly sensitive to the different physics models: galaxies simulated with a turbulence-based SFE have discs that are smoother and more rotationally symmetric, compared to those simulated with a constant SFE; galaxies simulated with early stellar feedback have more regular discs than supernova-feedback-only galaxies. We find that the rotational asymmetry of the HI discs depends most strongly on the underlying physics model, making this a promising observable for understanding the physics responsible for shaping the interstellar medium of galaxies. △ Less

Submitted 7 May, 2024; v1 submitted 2 October, 2023; originally announced October 2023.

Comments: 13 pages, 10 figures + appendices (7 pages, 10 figures); updated to version accepted by MNRAS

Star Cluster Classification using Deep Transfer Learning with PHANGS-HST

Authors: Stephen Hannon, Bradley C. Whitmore, Janice C. Lee, David A. Thilker, Sinan Deger, E. A. Huerta, Wei Wei, Bahram Mobasher, Ralf Klessen, Mederic Boquien, Daniel A. Dale, Melanie Chevance, Kathryn Grasha, Patricia Sanchez-Blazquez, Thomas Williams, Fabian Scheuermann, Brent Groves, Hwihyun Kim, J. M. Diederick Kruijssen, the PHANGS-HST Team

Abstract: Currently available star cluster catalogues from HST imaging of nearby galaxies heavily rely on visual inspection and classification of candidate clusters. The time-consuming nature of this process has limited the production of reliable catalogues and thus also post-observation analysis. To address this problem, deep transfer learning has recently been used to create neural network models which ac… ▽ More Currently available star cluster catalogues from HST imaging of nearby galaxies heavily rely on visual inspection and classification of candidate clusters. The time-consuming nature of this process has limited the production of reliable catalogues and thus also post-observation analysis. To address this problem, deep transfer learning has recently been used to create neural network models which accurately classify star cluster morphologies at production scale for nearby spiral galaxies (D < 20 Mpc). Here, we use HST UV-optical imaging of over 20,000 sources in 23 galaxies from the Physics at High Angular Resolution in Nearby GalaxieS (PHANGS) survey to train and evaluate two new sets of models: i) distance-dependent models, based on cluster candidates binned by galaxy distance (9-12 Mpc, 14-18 Mpc, 18-24 Mpc), and ii) distance-independent models, based on the combined sample of candidates from all galaxies. We find that the overall accuracy of both sets of models is comparable to previous automated star cluster classification studies (~60-80 per cent) and show improvement by a factor of two in classifying asymmetric and multi-peaked clusters from PHANGS-HST. Somewhat surprisingly, while we observe a weak negative correlation between model accuracy and galactic distance, we find that training separate models for the three distance bins does not significantly improve classification accuracy. We also evaluate model accuracy as a function of cluster properties such as brightness, colour, and SED-fit age. Based on the success of these experiments, our models will provide classifications for the full set of PHANGS-HST candidate clusters (N ~ 200,000) for public release. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 16 pages, 10 figures

On the origin of globular clusters in a hierarchical Universe

Authors: Gabriella De Lucia, J. M. Diederik Kruijssen, Sebastian Trujillo-Gomez, Michaela Hirschmann, Lizhi Xie

Abstract: We present an end-to-end description of the formation process of globular clusters (GCs) which combines a treatment for their formation and dynamical evolution within galaxy haloes with a state-of-the-art semi-analytic simulation of galaxy formation. Our approach allows us to obtain exquisite statistics to study the effect of the environment and assembly history of galaxies, while still allowing a… ▽ More We present an end-to-end description of the formation process of globular clusters (GCs) which combines a treatment for their formation and dynamical evolution within galaxy haloes with a state-of-the-art semi-analytic simulation of galaxy formation. Our approach allows us to obtain exquisite statistics to study the effect of the environment and assembly history of galaxies, while still allowing a very efficient exploration of the parameter space of star cluster physics. Our reference model, including both efficient cluster disruption during galaxy mergers and a model for the dynamical friction of GCs within the galactic potential, accurately reproduces the observed correlation between the total mass in GCs and the parent halo mass. A deviation from linearity is predicted at low halo masses, which is driven by a strong dependence on morphological type: bulge-dominated galaxies tend to host larger masses of GCs than their later-type counterparts of similar stellar mass. While the significance of the difference might be affected by resolution at the lowest halo masses considered, this is a robust prediction of our model and represents a natural consequence of the assumption that cluster migration from the disk to the halo is triggered by galaxy mergers. Our model requires an environmental dependence of GC radii to reproduce the observed mass distribution of GCs in our Galaxy at the low-mass end. At GC masses $>10^6\,{\rm M}_\odot$, our model predicts fewer GCs than observed due to an overly aggressive treatment of dynamical friction. The metallicity distribution measured for Galactic GCs is well reproduced by our model, even though it predicts systematically younger GCs than observed. We argue that this adds further evidence for an anomalously early formation of the stars in our Galaxy. △ Less

Submitted 11 April, 2024; v1 submitted 5 July, 2023; originally announced July 2023.

Comments: 18 pages, 16 figures, accepted for publication in MNRAS

The Gas Morphology of Nearby Star-Forming Galaxies

Authors: S. K. Stuber, E. Schinnerer, T. G. Williams, M. Querejeta, S. Meidt, E. Emsellem, A. Barnes, R. S. Klessen, A. K. Leroy, J. Neumann, M. C. Sormani, F. Bigiel, M. Chevance, D. Dale, C. Faesi, S. C. O. Glover, K. Grasha, J. M. D. Kruijssen, D. Liu, H. Pan, J. Pety, F. Pinna, T. Saito, A. Usero, E. J. Watkins

Abstract: The morphology of a galaxy stems from secular and environmental processes during its evolutionary history. Thus galaxy morphologies have been a long used tool to gain insights on galaxy evolution. We visually classify morphologies on cloud-scales based on the molecular gas distribution of a large sample of 79 nearby main-sequence galaxies, using 1'' resolution CO(2-1) ALMA observations taken as pa… ▽ More The morphology of a galaxy stems from secular and environmental processes during its evolutionary history. Thus galaxy morphologies have been a long used tool to gain insights on galaxy evolution. We visually classify morphologies on cloud-scales based on the molecular gas distribution of a large sample of 79 nearby main-sequence galaxies, using 1'' resolution CO(2-1) ALMA observations taken as part of the PHANGS survey. To do so, we devise a morphology classification scheme for different types of bars, spiral arms (grand-design, flocculent, multi-arm and smooth), rings (central and non-central rings) similar to the well-established optical ones, and further introduce bar lane classes. In general, our cold gas based morphologies agree well with the ones based on stellar light. Both our bars as well as grand-design spiral arms are preferentially found at the higher mass end of our sample. Our gas-based classification indicates a potential for misidentification of unbarred galaxies in the optical when massive star formation is present. Central or nuclear rings are present in a third of the sample with a strong preferences for barred galaxies (59%). As stellar bars are present in 45$\pm$5% of our sample galaxies, we explore the utility of molecular gas as tracer of bar lane properties. We find that more curved bar lanes have a shorter radial extent in molecular gas and reside in galaxies with lower molecular to stellar mass ratios than those with straighter geometries. Galaxies display a wide range of CO morphology, and this work provides a catalogue of morphological features in a representative sample of nearby galaxies. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: 17 pages, 14 figures (+ Appendix 9 pages, 4 figures). Accepted for publication in A&A

Journal ref: A&A 676, A113 (2023)

Fuelling the nuclear ring of NGC 1097

Authors: Mattia C. Sormani, Ashley T. Barnes, Jiayi Sun, Sophia K. Stuber, Eva Schinnerer, Eric Emsellem, Adam K. Leroy, Simon C. O. Glover, Jonathan D. Henshaw, Sharon E. Meidt, Justus Neumann, Miguel Querejeta, Thomas G. Williams, Frank Bigiel, Cosima Eibensteiner, Francesca Fragkoudi, Rebecca C. Levy, Kathryn Grasha, Ralf S. Klessen, J. M. Diederik Kruijssen, Nadine Neumayer, Francesca Pinna, Erik W. Rosolowsky, Rowan J. Smith, Yu-Hsuan Teng , et al. (2 additional authors not shown)

Abstract: Galactic bars can drive cold gas inflows towards the centres of galaxies. The gas transport happens primarily through the so-called bar ``dust lanes'', which connect the galactic disc at kpc scales to the nuclear rings at hundreds of pc scales much like two gigantic galactic rivers. Once in the ring, the gas can fuel star formation activity, galactic outflows, and central supermassive black holes.… ▽ More Galactic bars can drive cold gas inflows towards the centres of galaxies. The gas transport happens primarily through the so-called bar ``dust lanes'', which connect the galactic disc at kpc scales to the nuclear rings at hundreds of pc scales much like two gigantic galactic rivers. Once in the ring, the gas can fuel star formation activity, galactic outflows, and central supermassive black holes. Measuring the mass inflow rates is therefore important to understanding the mass/energy budget and evolution of galactic nuclei. In this work, we use CO datacubes from the PHANGS-ALMA survey and a simple geometrical method to measure the bar-driven mass inflow rate onto the nuclear ring of the barred galaxy NGC~1097. The method assumes that the gas velocity in the bar lanes is parallel to the lanes in the frame co-rotating with the bar, and allows one to derive the inflow rates from sufficiently sensitive and resolved position-position-velocity diagrams if the bar pattern speed and galaxy orientations are known. We find an inflow rate of $\dot{M}=(3.0 \pm 2.1)\, \rm M_\odot\, yr^{-1}$ averaged over a time span of 40 Myr, which varies by a factor of a few over timescales of $\sim$10 Myr. Most of the inflow appears to be consumed by star formation in the ring which is currently occurring at a rate of ${\rm SFR}\simeq~1.8$-$2 \rm M_\odot\, yr^{-1}$, suggesting that the inflow is causally controlling the star formation rate in the ring as a function of time. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: Accepted in MNRAS

Kinematics of Galactic Centre clouds shaped by shear-seeded solenoidal turbulence

Authors: Maya A. Petkova, J. M. Diederik Kruijssen, Jonathan D. Henshaw, Steven N. Longmore, Simon C. O. Glover, Mattia C. Sormani, Lucia Armillotta, Ashley T. Barnes, Ralf S. Klessen, Francisco Nogueras-Lara, Robin G. Tress, Jairo Armijos-Abendaño, Laura Colzi, Christoph Federrath, Pablo García, Adam Ginsburg, Christian Henkel, Sergio Martín, Denise Riquelme, Víctor M. Rivilla

Abstract: The Central Molecular Zone (CMZ; the central ~ 500 pc of the Galaxy) is a kinematically unusual environment relative to the Galactic disc, with high velocity dispersions and a steep size-linewidth relation of the molecular clouds. In addition, the CMZ region has a significantly lower star formation rate (SFR) than expected by its large amount of dense gas. An important factor in explaining the low… ▽ More The Central Molecular Zone (CMZ; the central ~ 500 pc of the Galaxy) is a kinematically unusual environment relative to the Galactic disc, with high velocity dispersions and a steep size-linewidth relation of the molecular clouds. In addition, the CMZ region has a significantly lower star formation rate (SFR) than expected by its large amount of dense gas. An important factor in explaining the low SFR is the turbulent state of the star-forming gas, which seems to be dominated by rotational modes. However, the turbulence driving mechanism remains unclear. In this work, we investigate how the Galactic gravitational potential affects the turbulence in CMZ clouds. We focus on the CMZ cloud G0.253+0.016 (`the Brick'), which is very quiescent and unlikely to be kinematically dominated by stellar feedback. We demonstrate that several kinematic properties of the Brick arise naturally in a cloud-scale hydrodynamics simulation that takes into account the Galactic gravitational potential. These properties include the line-of-sight velocity distribution, the steepened size-linewidth relation, and the predominantly solenoidal nature of the turbulence. Within the simulation, these properties result from the Galactic shear in combination with the cloud's gravitational collapse. This is a strong indication that the Galactic gravitational potential plays a crucial role in shaping the CMZ gas kinematics, and is a major contributor to suppressing the SFR by inducing predominantly solenoidal turbulent modes. △ Less

Submitted 1 August, 2023; v1 submitted 21 April, 2023; originally announced April 2023.

Comments: 7 pages, 8 figures; accepted to MNRAS (July 24th 2023)

Kinematic analysis of the super-extended HI disk of the nearby spiral galaxy M83

Authors: Cosima Eibensteiner, Frank Bigiel, Adam K. Leroy, Eric W. Koch, Erik Rosolowsky, Eva Schinnerer, Amy Sardone, Sharon Meidt, W. J. G de Blok, David Thilker, D. J. Pisano, Jürgen Ott, Ashley Barnes, Miguel Querejeta, Eric Emsellem, Johannes Puschnig, Dyas Utomo, Ivana Bešlic, Jakob den Brok, Shahram Faridani, Simon C. O. Glover, Kathryn Grasha, Hamid Hassani, Jonathan D. Henshaw, Maria J. Jiménez-Donaire , et al. (11 additional authors not shown)

Abstract: We present new HI observations of the nearby massive spiral galaxy M83, taken with the VLA at $21^{\prime\prime}$ angular resolution ($\approx500$ pc) of an extended ($\sim$1.5 deg$^2$) 10-point mosaic combined with GBT single dish data. We study the super-extended HI disk of M83 (${\sim}$50 kpc in radius), in particular disc kinematics, rotation and the turbulent nature of the atomic interstellar… ▽ More We present new HI observations of the nearby massive spiral galaxy M83, taken with the VLA at $21^{\prime\prime}$ angular resolution ($\approx500$ pc) of an extended ($\sim$1.5 deg$^2$) 10-point mosaic combined with GBT single dish data. We study the super-extended HI disk of M83 (${\sim}$50 kpc in radius), in particular disc kinematics, rotation and the turbulent nature of the atomic interstellar medium. We define distinct regions in the outer disk ($r_{\rm gal}>$central optical disk), including ring, southern area, and southern and northern arm. We examine HI gas surface density, velocity dispersion and non-circular motions in the outskirts, which we compare to the inner optical disk. We find an increase of velocity dispersion ($σ_v$) towards the pronounced HI ring, indicative of more turbulent HI gas. Additionally, we report over a large galactocentric radius range (until $r_{\rm gal}{\sim}$50 kpc) that $σ_v$ is slightly larger than thermal (i.e. $>8$km s$^{-1}$ ). We find that a higher star formation rate (as traced by FUV emission) is not always necessarily associated with a higher HI velocity dispersion, suggesting that radial transport could be a dominant driver for the enhanced velocity dispersion. We further find a possible branch that connects the extended HI disk to the dwarf irregular galaxy UGCA365, that deviates from the general direction of the northern arm. Lastly, we compare mass flow rate profiles (based on 2D and 3D tilted ring models) and find evidence for outflowing gas at r$_{\rm gal}$ $\sim$2 kpc, inflowing gas at r$_{\rm gal}$ $\sim$5.5 kpc and outflowing gas at r$_{\rm gal}$ $\sim$14 kpc. We caution that mass flow rates are highly sensitive to the assumed kinematic disk parameters, in particular, to the inclination. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: accepted for publication in A&A; 16 pages, 12 figures (+8 pages appendix)

Journal ref: A&A 675, A37 (2023)

Resolved stellar population properties of PHANGS-MUSE galaxies

Authors: I. Pessa, E. Schinnerer, P. Sanchez-Blazquez, F. Belfiore, B. Groves, E. Emsellem, J. Neumann, A. K. Leroy, F. Bigiel, M. Chevance, D. A. Dale, S. C. O. Glover, K. Grasha, R. S. Klessen, K. Kreckel, J. M. D. Kruijssen, F. Pinna, M. Querejeta, E. Rosolowsky, T. G. Williams

Abstract: Analyzing resolved stellar populations across the disk of a galaxy can provide unique insights into how that galaxy assembled its stellar mass over its lifetime. Previous work at ~1 kpc resolution has already revealed common features in the mass buildup (e.g., inside-out growth of galaxies). However, even at approximate kpc scales, the stellar populations are blurred between the different galactic… ▽ More Analyzing resolved stellar populations across the disk of a galaxy can provide unique insights into how that galaxy assembled its stellar mass over its lifetime. Previous work at ~1 kpc resolution has already revealed common features in the mass buildup (e.g., inside-out growth of galaxies). However, even at approximate kpc scales, the stellar populations are blurred between the different galactic morphological structures such as spiral arms, bars and bulges. Here we present a detailed analysis of the spatially resolved star formation histories (SFHs) of 19 PHANGS-MUSE galaxies, at a spatial resolution of ~100 pc. We show that our sample of local galaxies exhibits predominantly negative radial gradients of stellar age and [Z/H], consistent with previous findings, and a radial structure that is primarily consistent with local star formation, and indicative of inside-out formation. In barred galaxies, we find flatter [Z/H] gradients along the semi-major axis of the bar than along the semi-minor axis, as is expected from the radial mixing of material along the bar. In general, the derived assembly histories of the galaxies in our sample tell a consistent story of inside-out growth, where low-mass galaxies assembled the majority of their stellar mass later in cosmic history than high-mass galaxies. We also show how stellar populations of different ages exhibit different kinematics, with younger stellar populations having lower velocity dispersions than older stellar populations at similar galactocentric distances, which we interpret as an imprint of the progressive dynamical heating of stellar populations as they age. Finally, we explore how the time-averaged star formation rate evolves with time, and how it varies across galactic disks. This analysis reveals a wide variation of the SFHs of galaxy centers and additionally shows that structural features become less pronounced with age. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: 52 pages, 48 figures, accepted for publication in A&A

Journal ref: A&A 673, A147 (2023)

Stellar associations powering HII regions $\unicode{x2013}$ I. Defining an evolutionary sequence

Authors: Fabian Scheuermann, Kathryn Kreckel, Ashley T. Barnes, Francesco Belfiore, Brent Groves, Stephen Hannon, Janice C. Lee, Rebecca Minsley, Erik Rosolowsky, Frank Bigiel, Guillermo A. Blanc, Médéric Boquien, Daniel A. Dale, Sinan Deger, Oleg V. Egorov, Eric Emsellem, Simon C. O. Glover, Kathryn Grasha, Hamid Hassani, Sarah Jeffreson, Ralf S. Klessen, J. M. Diederik Kruijssen, Kirsten L. Larson, Adam K. Leroy, Laura Lopez , et al. (8 additional authors not shown)

Abstract: Connecting the gas in HII regions to the underlying source of the ionizing radiation can help us constrain the physical processes of stellar feedback and how HII regions evolve over time. With PHANGS$\unicode{x2013}$MUSE we detect nearly 24,000 HII regions across 19 galaxies and measure the physical properties of the ionized gas (e.g. metallicity, ionization parameter, density). We use catalogues… ▽ More Connecting the gas in HII regions to the underlying source of the ionizing radiation can help us constrain the physical processes of stellar feedback and how HII regions evolve over time. With PHANGS$\unicode{x2013}$MUSE we detect nearly 24,000 HII regions across 19 galaxies and measure the physical properties of the ionized gas (e.g. metallicity, ionization parameter, density). We use catalogues of multi-scale stellar associations from PHANGS$\unicode{x2013}$HST to obtain constraints on the age of the ionizing sources. We construct a matched catalogue of 4,177 HII regions that are clearly linked to a single ionizing association. A weak anti-correlation is observed between the association ages and the H$α$ equivalent width EW(H$α$), the H$α$/FUV flux ratio and the ionization parameter, log q. As all three are expected to decrease as the stellar population ages, this could indicate that we observe an evolutionary sequence. This interpretation is further supported by correlations between all three properties. Interpreting these as evolutionary tracers, we find younger nebulae to be more attenuated by dust and closer to giant molecular clouds, in line with recent models of feedback-regulated star formation. We also observe strong correlations with the local metallicity variations and all three proposed age tracers, suggestive of star formation preferentially occurring in locations of locally enhanced metallicity. Overall, EW(H$α$) and log q show the most consistent trends and appear to be most reliable tracers for the age of an HII region. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: 15 pages, 12 figures. Accepted for publication in MNRAS

Star Formation Laws and Efficiencies across 80 Nearby Galaxies

Authors: Jiayi Sun, Adam K. Leroy, Eve C. Ostriker, Sharon Meidt, Erik Rosolowsky, Eva Schinnerer, Christine D. Wilson, Dyas Utomo, Francesco Belfiore, Guillermo A. Blanc, Eric Emsellem, Christopher Faesi, Brent Groves, Annie Hughes, Eric W. Koch, Kathryn Kreckel, Daizhong Liu, Hsi-An Pan, Jerome Pety, Miguel Querejeta, Alessandro Razza, Toshiki Saito, Amy Sardone, Antonio Usero, Thomas G. Williams , et al. (15 additional authors not shown)

Abstract: We measure empirical relationships between the local star formation rate (SFR) and properties of the star-forming molecular gas on 1.5 kpc scales across 80 nearby galaxies. These relationships, commonly referred to as "star formation laws," aim at predicting the local SFR surface density from various combinations of molecular gas surface density, galactic orbital time, molecular cloud free-fall ti… ▽ More We measure empirical relationships between the local star formation rate (SFR) and properties of the star-forming molecular gas on 1.5 kpc scales across 80 nearby galaxies. These relationships, commonly referred to as "star formation laws," aim at predicting the local SFR surface density from various combinations of molecular gas surface density, galactic orbital time, molecular cloud free-fall time, and the interstellar medium dynamical equilibrium pressure. Leveraging a multiwavelength database built for the PHANGS survey, we measure these quantities consistently across all galaxies and quantify systematic uncertainties stemming from choices of SFR calibrations and the CO-to-H$_2$ conversion factors. The star formation laws we examine show 0.3-0.4 dex of intrinsic scatter, among which the molecular Kennicutt-Schmidt relation shows a $\sim$10% larger scatter than the other three. The slope of this relation ranges $β\approx0.9{-}1.2$, implying that the molecular gas depletion time remains roughly constant across the environments probed in our sample. The other relations have shallower slopes ($β\approx0.6{-}1.0$), suggesting that the star formation efficiency (SFE) per orbital time, the SFE per free-fall time, and the pressure-to-SFR surface density ratio (i.e., the feedback yield) may vary systematically with local molecular gas and SFR surface densities. Last but not least, the shapes of the star formation laws depend sensitively on methodological choices. Different choices of SFR calibrations can introduce systematic uncertainties of at least 10-15% in the star formation law slopes and 0.15-0.25 dex in their normalization, while the CO-to-H$_2$ conversion factors can additionally produce uncertainties of 20-25% for the slope and 0.10-0.20 dex for the normalization. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 10 pages main text + 2 appendices. ApJL in press. Data products available at https://www.canfar.net/storage/list/phangs/RELEASES/Sun_etal_2023 . Slides summarizing key results can be found at https://www.dropbox.com/s/5gsegexeo9n0t05/Sun_et_PHANGS_2023.pptx?dl=0

Journal ref: ApJL, 945, L19 (2023)

Quantifying the energetics of molecular superbubbles in PHANGS galaxies

Authors: E. J. Watkins, K. Kreckel, B. Groves, S. C. O. Glover, B. C. Whitmore, A. K. Leroy, E. Schinnerer, S. E. Meidt, O. V. Egorov, A. T. Barnes, J. C. Lee, F. Bigiel, M. Boquien, R. Chandar, M. Chevance, D. A. Dale, K. Grasha, R. S. Klessen, J. M. D. Kruijssen, K. L. Larson, J. Li, J. E. Méndez-Delgado, I. Pessa, T. Saito, P. Sanchez-Blazquez , et al. (4 additional authors not shown)

Abstract: Star formation and stellar feedback are interlinked processes that redistribute energy and matter throughout galaxies. When young, massive stars form in spatially clustered environments, they create pockets of expanding gas termed superbubbles. As these processes play a critical role in shaping galaxy discs and regulating the baryon cycle, measuring the properties of superbubbles provides importan… ▽ More Star formation and stellar feedback are interlinked processes that redistribute energy and matter throughout galaxies. When young, massive stars form in spatially clustered environments, they create pockets of expanding gas termed superbubbles. As these processes play a critical role in shaping galaxy discs and regulating the baryon cycle, measuring the properties of superbubbles provides important input for galaxy evolution models. With wide coverage and high angular resolution (50-150 pc) of the PHANGS-ALMA $^{12}$CO (2-1) survey, we can now resolve and identify a statistically representative number of superbubbles with molecular gas in nearby galaxies. We identify superbubbles by requiring spatial correspondence between shells in CO with stellar populations identified in PHANGS-HST, and combine the properties of the stellar populations with CO to constrain feedback models and quantify their energetics. We visually identify 325 cavities across 18 PHANGS-ALMA galaxies, 88 of which have clear superbubble signatures (unbroken shells, central clusters, kinematic signatures of expansion). We measure their radii and expansion velocities using CO to dynamically derive their ages and the mechanical power driving the bubbles, which we use to compute the expected properties of the parent stellar populations driving the bubbles. We find consistency between the predicted and derived stellar ages and masses of the stellar populations if we use a supernova blast wave model that injects energy with a coupling efficiency of 10%, whereas continuous models fail to explain stellar ages we measure. Not only does this confirm molecular gas accurately traces superbubble properties, but it also provides key observational constraints for superbubble models. We also find evidence that the bubbles sweep up gas as they expand and speculate that these sites have the potential to host new generations of stars. △ Less

Submitted 14 June, 2023; v1 submitted 7 February, 2023; originally announced February 2023.

Comments: 21 pages, 15 figures, 3 tables. Accepted to A&A. Abstract abridged for arXiv

Journal ref: A&A 676, A67 (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

In-situ or accreted? Using deep learning to infer the origin of extragalactic globular clusters from observables

Authors: Sebastian Trujillo-Gomez, J. M. Diederik Kruijssen, Joel Pfeffer, Marta Reina-Campos, Robert A. Crain, Nate Bastian, Ivan Cabrera-Ziri

Abstract: Globular clusters (GCs) are powerful tracers of the galaxy assembly process, and have already been used to obtain a detailed picture of the progenitors of the Milky Way. Using the E-MOSAICS cosmological simulation of a (34.4 Mpc)$^3$ volume that follows the formation and co-evolution of galaxies and their star cluster populations, we develop a method to link the origin of GCs to their observable p… ▽ More Globular clusters (GCs) are powerful tracers of the galaxy assembly process, and have already been used to obtain a detailed picture of the progenitors of the Milky Way. Using the E-MOSAICS cosmological simulation of a (34.4 Mpc)$^3$ volume that follows the formation and co-evolution of galaxies and their star cluster populations, we develop a method to link the origin of GCs to their observable properties. We capture this complex link using a supervised deep learning algorithm trained on the simulations, and predict the origin of individual GCs (whether they formed in the main progenitor or were accreted from satellites) based solely on extragalactic observables. An artificial neural network classifier trained on $\sim50,000$ GCs hosted by $\sim 700$ simulated galaxies successfully predicts the origin of GCs in the test set with a mean accuracy of $89$ per cent for the objects with [Fe/H]<-0.5 that have unambiguous classifications. The network relies mostly on the alpha-element abundances, metallicities, projected positions, and projected angular momenta of the clusters to predict their origin. A real-world test using the known progenitor associations of the Milky Way GCs achieves up to $90$ per cent accuracy, and successfully identifies as accreted most of the GCs in the inner Galaxy associated to the Kraken progenitor, as well as all the Gaia-Enceladus GCs. We demonstrate that the model is robust to observational uncertainties, and develop a method to predict the classification accuracy across observed galaxies. The classifier can be optimized for available observables (e.g. to improve the accuracy by including GC ages), making it a valuable tool to reconstruct the assembly histories of galaxies in upcoming wide-field surveys. △ Less

Submitted 13 January, 2023; originally announced January 2023.

Comments: 21 pages, 22 figures. Submitted to MNRAS. Comments welcome

Improving Star Cluster Age Estimates in PHANGS-HST Galaxies and the Impact on Cluster Demographics in NGC 628

Authors: Bradley C. Whitmore, Rupali Chandar, Janice C. Lee, Matthew Floyd, Sinan Deger, James Lilly, Rebecca Minsley, David A. Thilker, Médéric Boquien, Daniel A. Dale, Kiana Henny, Fabian Scheuermann, Ashley T. Barnes, Frank Bigiel, Eric Emsellem, Simon Glover, Kathryn Grasha, Brent Groves, Stephen Hannon, Ralf S. Klessen, Kathryn Kreckel, J. M. Diederik Kruijssen, Kirsten L. Larson, Adam Leroy, Angus Mok , et al. (7 additional authors not shown)

Abstract: A long-standing problem when deriving the physical properties of stellar populations is the degeneracy between age, reddening, and metallicity. When a single metallicity is used for all star clusters in a galaxy, this degeneracy can result in $`$catastrophic$'$ errors for old globular clusters. Typically, approximately 10 - 20 % of all clusters detected in spiral galaxies can have ages that are in… ▽ More A long-standing problem when deriving the physical properties of stellar populations is the degeneracy between age, reddening, and metallicity. When a single metallicity is used for all star clusters in a galaxy, this degeneracy can result in $`$catastrophic$'$ errors for old globular clusters. Typically, approximately 10 - 20 % of all clusters detected in spiral galaxies can have ages that are incorrect by a factor of ten or more. In this paper we present a pilot study for four galaxies (NGC 628, NGC 1433, NGC 1365, and NGC 3351) from the PHANGS-HST survey. We describe methods to correct the age-dating for old globular clusters, by first identifying candidates using their colors, and then reassigning ages and reddening based on a lower metallicity solution. We find that young $`$interlopers$'$ can be identified from their Halpha flux. CO (2-1) intensity or the presence of dust can also be used, but our tests show that they do not work as well. Improvements in the success fraction are possible at the $\sim$ 15 % level (reducing the fraction of catastrophic age-estimates from between 13 - 21 % to 3 - 8 %). A large fraction of the incorrectly age-dated globular clusters are systematically given ages around 100 Myr, polluting the younger populations as well. Incorrectly age-dated globular clusters significantly impact the observed cluster age distribution in NGC 628, which affects the physical interpretation of cluster disruption in this galaxy. For NGC 1365, we also demonstrate how to fix a second major age-dating problem, where very dusty young clusters with E(B-V) $>$ 1.5 mag are assigned old, globular-cluster like ages. Finally, we note the discovery of a dense population of $\sim$ 300 Myr clusters around the central region of NGC 1365. and discuss how this results naturally from the dynamics in a barred galaxy. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Comments: 27 pages, 17 figures (NOTE: 1, 5, 9, 12, 14 are lower resolution than in the journal)

PHANGS-JWST First Results: The Dust Filament Network of NGC 628 and its Relation to Star Formation Activity

Authors: David A. Thilker, Janice C. Lee, Sinan Deger, Ashley T. Barnes, Frank Bigiel, Médéric Boquien, Yixian Cao, Mélanie Chevance, Daniel A. Dale, Oleg V. Egorov, Simon C. O. Glover, Kathryn Grasha, Jonathan D. Henshaw, Ralf S. Klessen, Eric Koch, J. M. Diederik Kruijssen, Adam K. Leroy, Ryan A. Lessing, Sharon E. Meidt, Francesca Pinna, Miguel Querejeta, Erik Rosolowsky, Karin M. Sandstrom, Eva Schinnerer, Rowan J. Smith , et al. (14 additional authors not shown)

Abstract: PHANGS-JWST mid-infrared (MIR) imaging of nearby spiral galaxies has revealed ubiquitous filaments of dust emission in intricate detail. We present a pilot study to systematically map the dust filament network (DFN) at multiple scales between 25-400 pc in NGC 628. MIRI images at 7.7, 10, 11.3 and 21$μ$m of NGC 628 are used to generate maps of the filaments in emission, while PHANGS-HST B-band imag… ▽ More PHANGS-JWST mid-infrared (MIR) imaging of nearby spiral galaxies has revealed ubiquitous filaments of dust emission in intricate detail. We present a pilot study to systematically map the dust filament network (DFN) at multiple scales between 25-400 pc in NGC 628. MIRI images at 7.7, 10, 11.3 and 21$μ$m of NGC 628 are used to generate maps of the filaments in emission, while PHANGS-HST B-band imaging yields maps of dust attenuation features. We quantify the correspondence between filaments traced by MIR thermal continuum / polycyclic aromatic hydrocarbon (PAH) emission and filaments detected via extinction / scattering of visible light; the fraction of MIR flux contained in the DFN; and the fraction of HII regions, young star clusters and associations within the DFN. We examine the dependence of these quantities with the physical scale at which the DFN is extracted. With our highest resolution DFN maps (25 pc filament width), we find that filaments in emission and attenuation are co-spatial in 40% of sight lines, often exhibiting detailed morphological agreement; that ~30% of the MIR flux is associated with the DFN; and that 75-80% of HII regions and 60% of star clusters younger than 5 Myr are contained within the DFN. However, the DFN at this scale is anti-correlated with looser associations of stars younger than 5 Myr identified using PHANGS-HST near-UV imaging. We discuss the impact of these findings for studies of star formation and the ISM, and the broad range of new investigations enabled with multi-scale maps of the DFN. △ Less

Submitted 2 January, 2023; originally announced January 2023.

Comments: 21 pages, 14 figures, accepted for publication as part of PHANGS-JWST ApJL Focus Issue

PHANGS-JWST First Results: Mapping the 3.3 micron Polycyclic Aromatic Hydrocarbon Vibrational Band in Nearby Galaxies with NIRCam Medium Bands

Authors: Karin Sandstrom, Jérémy Chastenet, Jessica Sutter, Adam K. Leroy, Oleg V. Egorov, Thomas G. Williams, Alberto D. Bolatto, Médéric Boquien, Yixian Cao, Daniel A. Dale, Janice C. Lee, Erik Rosolowsky, Eva Schinnerer, Ashley. T. Barnes, Francesco Belfiore, F. Bigiel, Mélanie Chevance, Kathryn Grasha, Brent Groves, Hamid Hassani, Annie Hughes, Ralf S. Klessen, J. M. Diederik Kruijssen, Kirsten L. Larson, Daizhong Liu , et al. (6 additional authors not shown)

Abstract: We present maps of the 3.3 micron polycyclic aromatic hydrocarbon (PAH) emission feature in NGC 628, NGC 1365, and NGC 7496 as observed with the Near-Infrared Camera (NIRCam) imager on JWST from the PHANGS-JWST Cycle 1 Treasury project. We create maps that isolate the 3.3 micron PAH feature in the F335M filter (F335M$_{\rm PAH}$) using combinations of the F300M and F360M filters for removal of sta… ▽ More We present maps of the 3.3 micron polycyclic aromatic hydrocarbon (PAH) emission feature in NGC 628, NGC 1365, and NGC 7496 as observed with the Near-Infrared Camera (NIRCam) imager on JWST from the PHANGS-JWST Cycle 1 Treasury project. We create maps that isolate the 3.3 micron PAH feature in the F335M filter (F335M$_{\rm PAH}$) using combinations of the F300M and F360M filters for removal of starlight continuum. This continuum removal is complicated by contamination of the F360M by PAH emission and variations in the stellar spectral energy distribution slopes between 3.0 and 3.6 micron. We modify the empirical prescription from Lai et al. (2020) to remove the starlight continuum in our highly resolved galaxies, which have a range of starlight- and PAH-dominated lines-of-sight. Analyzing radially binned profiles of the F335M$_{\rm PAH}$ emission, we find that between 5-65% of the F335M intensity comes from the 3.3 micron feature within the inner 0.5 $r_{25}$ of our targets. This percentage systematically varies from galaxy to galaxy, and shows radial trends within the galaxies related to each galaxy's distribution of stellar mass, interstellar medium, and star formation. The 3.3 micron emission is well correlated with the 11.3 micron PAH feature traced with the MIRI F1130W filter, as is expected, since both features arise from C-H vibrational modes. The average F335M$_{\rm PAH}$/F1130W ratio agrees with the predictions of recent models by Draine et al. (2021) for PAHs with size and charge distributions shifted towards larger grains with normal or higher ionization. △ Less

Submitted 2 January, 2023; originally announced January 2023.

Comments: submitted to AAS journals and revised according to referee comments, part of a PHANGS-JWST Focus Issue to appear in ApJ

PHANGS-JWST First Results: Variations in PAH Fraction as a Function of ISM Phase and Metallicity

Authors: Jérémy Chastenet, Jessica Sutter, Karin Sandstrom, Francesco Belfiore, Oleg V. Egorov, Kirsten L. Larson, Adam K. Leroy, Daizhong Liu, Erik Rosolowsky, David A. Thilker, Elizabeth J. Watkins, Thomas G. Williams, Ashley T. Barnes, Frank Bigiel, Médéric Boquien, Mélanie Chevance, I-Da Chiang, Daniel A. Dale, J. M. Diederik Kruijssen, Eric Emsellem, Kathryn Grasha, Brent Groves, Hamid Hassani, Annie Hughes, Kathryn Kreckel , et al. (4 additional authors not shown)

Abstract: We present maps tracing the fraction of dust in the form of polycyclic aromatic hydrocarbons (PAHs) in IC 5332, NGC 628, NGC 1365, and NGC 7496 from JWST/MIRI observations. We trace the PAH fraction by combining the F770W ($7.7~μ$m) and F1130W ($11.3~μ$m) filters to track ionized and neutral PAH emission, respectively, and comparing the PAH emission to F2100W which traces small, hot dust grains. W… ▽ More We present maps tracing the fraction of dust in the form of polycyclic aromatic hydrocarbons (PAHs) in IC 5332, NGC 628, NGC 1365, and NGC 7496 from JWST/MIRI observations. We trace the PAH fraction by combining the F770W ($7.7~μ$m) and F1130W ($11.3~μ$m) filters to track ionized and neutral PAH emission, respectively, and comparing the PAH emission to F2100W which traces small, hot dust grains. We find average $R{\rm_{PAH} = (F770W+F1130W)/F2100W}$ values of 3.3, 4.7, 5.1, and 3.6 in IC 5332, NGC 628, NGC 1365, and NGC 7496, respectively. We find that H II regions traced by MUSE H$α$ show a systematically low PAH fraction. The PAH fraction remains relatively constant across other galactic environments, with slight variations. We use CO + H I + H$α$ to trace the interstellar gas phase and find that the PAH fraction decreases above a value of I$_{Hα}/Σ_{H~I+H_2}$ $\sim~10^{37.5}$ erg s$^{-1}$ kpc$^{-2}$ (M$_\odot$ pc$^{-2}$)$^{-1}$, in all four galaxies. Radial profiles also show a decreasing PAH fraction with increasing radius, correlated with lower metallicity, in line with previous results showing a strong metallicity dependence to the PAH fraction. Our results suggest that the process of PAH destruction in ionized gas operates similarly across the four targets. △ Less

Submitted 3 January, 2023; v1 submitted 2 January, 2023; originally announced January 2023.

Comments: 13 pages, 4 figures. Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

PHANGS-JWST First Results: Massive Young Star Clusters and New Insights from JWST Observations of NGC 1365

Authors: Bradley C. Whitmore, Rupali Chandar, M. Jimena Rodríguez, Janice C. Lee, Eric Emsellem, Matthew Floyd, Hwihyun Kim, J. M. Diederik Kruijssen, Angus Mok, Mattia C. Sormani, Médéric Bodquien, Daniel A. Dale, Christopher M. Faesi, Kiana F. Henny, Stephen Hannon, David A. Thilker, Richad L. White, Ashley T. Barnes, F. Bigiel, Mélanie Chevance, Jonathan D. Henshaw, Ralf S. Klessen, Adam K. Leroy, Daizhong Liu, Daniel Maschmann , et al. (6 additional authors not shown)

Abstract: A primary new capability of JWST is the ability to penetrate the dust in star forming galaxies to identify and study the properties of young star clusters that remain embedded in dust and gas. In this paper we combine new infrared images taken with JWST with our optical HST images of the star-bursting barred (Seyfert2) spiral galaxy NGC 1365. We find that this galaxy has the richest population of… ▽ More A primary new capability of JWST is the ability to penetrate the dust in star forming galaxies to identify and study the properties of young star clusters that remain embedded in dust and gas. In this paper we combine new infrared images taken with JWST with our optical HST images of the star-bursting barred (Seyfert2) spiral galaxy NGC 1365. We find that this galaxy has the richest population of massive young clusters of any known galaxy within 30 Mpc, with $\sim$ 30 star clusters that are more massive than 10$^6$ Msolar and younger than 10 Myr. Sixteen of these clusters are newly discovered from our JWST observations. An examination of the optical images reveals that 4 of 30 ($\sim$13$\%$) are so deeply embedded that they cannot be seen in the I band (AV $\gt$ 10 mag), and that 11 of 30 ($\sim$37$\%$) are missing in the HST B band, so age and mass estimates from optical measurements alone are challenging. These numbers suggest that massive clusters in NGC 1365 remain obscured in the visible for $\sim$ 1.3 $\pm$ 0.7 Myr, and are either completely or partially obscured for $\sim$ 3.7 $\pm$ 1.1 Myr. We also use the JWST observations to gain new insights into the triggering of star cluster formation by the collision of gas and dust streamers with gas and dust in the bar. The JWST images reveal previously unknown structures (e.g., bridges and overshoot regions from stars that form in the bar) that help us better understand the orbital dynamics of barred galaxies and associated star-forming rings. Finally, we note that the excellent spatial resolution of the NIRCAM F200W filter provides a better way to separate barely resolved compact clusters from individual stars based on their sizes. △ Less

Submitted 29 December, 2022; v1 submitted 22 December, 2022; originally announced December 2022.

Comments: 17 pages, 9 figures, accepted for publications as part of PHANGS-JWST ApJL Focus Issue

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

PHANGS-JWST First Results: Tracing the Diffuse ISM with JWST Imaging of Polycyclic Aromatic Hydrocarbon Emission in Nearby Galaxies

Authors: Karin M. Sandstrom, Eric W. Koch, Adam K. Leroy, Erik Rosolowsky, Eric Emsellem, Rowan J. Smith, Oleg V. Egorov, Thomas G. Williams, Kirsten L. Larson, Janice C. Lee, Eva Schinnerer, David A. Thilker, Ashley. T. Barnes, Francesco Belfiore, F. Bigiel, Guillermo A. Blanc, Alberto D. Bolatto, Médéric Boquien, Yixian Cao, Jérémy Chastenet, Mélanie Chevance, I-Da Chiang, Daniel A. Dale, Christopher M. Faesi, Simon C. O. Glover , et al. (21 additional authors not shown)

Abstract: JWST observations of polycyclic aromatic hydrocarbon (PAH) emission provide some of the deepest and highest resolution views of the cold interstellar medium (ISM) in nearby galaxies. If PAHs are well mixed with the atomic and molecular gas and illuminated by the average diffuse interstellar radiation field, PAH emission may provide an approximately linear, high resolution, high sensitivity tracer… ▽ More JWST observations of polycyclic aromatic hydrocarbon (PAH) emission provide some of the deepest and highest resolution views of the cold interstellar medium (ISM) in nearby galaxies. If PAHs are well mixed with the atomic and molecular gas and illuminated by the average diffuse interstellar radiation field, PAH emission may provide an approximately linear, high resolution, high sensitivity tracer of diffuse gas surface density. We present a pilot study that explores using PAH emission in this way based on MIRI observations of IC 5332, NGC 628, NGC 1365, and NGC 7496 from the PHANGS-JWST Treasury. Using scaling relationships calibrated in Leroy et al. (2022), scaled F1130W provides 10--40 pc resolution and 3$σ$ sensitivity of $Σ_{\rm gas} \sim 2$ M$_\odot$ pc$^{-2}$. We characterize the surface densities of structures seen at $< 7$ M$_\odot$ pc$^{-2}$ in our targets, where we expect the gas to be HI-dominated. We highlight the existence of filaments, inter-arm emission, and holes in the diffuse ISM at these low surface densities. Below $\sim 10$ M$_\odot$ pc$^{-2}$ for NGC 628, NGC 1365, and NGC 7496 the gas distribution shows a ``Swiss cheese''-like topology due to holes and bubbles pervading the relatively smooth distribution of diffuse ISM. Comparing to recent galaxy simulations, we observe similar topology for the low surface density gas, though with notable variations between simulations with different setups and resolution. Such a comparison of high resolution, low surface density gas with simulations is not possible with existing atomic and molecular gas maps, highlighting the unique power of JWST maps of PAH emission. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: 18 pages, 6 figures, accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

PHANGS-JWST First Results: Mid-infrared emission traces both gas column density and heating at 100 pc scales

Authors: Adam K. Leroy, Karin Sandstrom, Erik Rosolowsky, Francesco Belfiore, Alberto D. Bolatto, Yixian Cao, Eric W. Koch, Eva Schinnerer, Ashley. T. Barnes, Ivana Bešlić, F. Bigiel, Guillermo A. Blanc, Jérémy Chastenet, Ness Mayker Chen, Mélanie Chevance, Ryan Chown, Enrico Congiu, Daniel A. Dale, Oleg V. Egorov, Eric Emsellem, Cosima Eibensteiner, Christopher M. Faesi, Simon C. O. Glover, Kathryn Grasha, Brent Groves , et al. (26 additional authors not shown)

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

Submitted 6 January, 2023; v1 submitted 20 December, 2022; originally announced December 2022.

Comments: 49 pages, 17 figures, Section 8 provides a detailed summary, accepted for publication in ApJ Letters, part of a PHANGS-JWST Focus Issue to appear in ApJ

PHANGS-JWST First Results: Measuring PAH Properties across the multiphase ISM

Authors: Jérémy Chastenet, Jessica Sutter, Karin Sandstrom, Francesco Belfiore, Oleg V. Egorov, Kirsten L. Larson, Adam K. Leroy, Daizhong Liu, Erik Rosolowsky, David A. Thilker, Elizabeth J. Watkins, Thomas G. Williams, Ashley T. Barnes, Frank Bigiel, Médéric Boquien, Mélanie Chevance, Daniel A. Dale, J. M. Diederik Kruijssen, Eric Emsellem, Kathryn Grasha, Brent Groves, Hamid Hassani, Annie Hughes, Kathryn Kreckel, Sharon E. Meidt , et al. (4 additional authors not shown)

Abstract: Ratios of polycyclic aromatic hydrocarbon (PAH) vibrational bands are a promising tool for measuring the properties of the PAH population and their effect on star formation. The photometric bands of the MIRI and NIRCam instruments on JWST provide the opportunity to measure PAH emission features across entire galaxy disks at unprecedented resolution and sensitivity. Here we present the first result… ▽ More Ratios of polycyclic aromatic hydrocarbon (PAH) vibrational bands are a promising tool for measuring the properties of the PAH population and their effect on star formation. The photometric bands of the MIRI and NIRCam instruments on JWST provide the opportunity to measure PAH emission features across entire galaxy disks at unprecedented resolution and sensitivity. Here we present the first results of this analysis in a sample of three nearby galaxies: NGC 628, NGC 1365, and NGC 7496. Based on the variations observed in the 3.3, 7.7, and 11.3 $μ$m features, we infer changes to the average PAH size and ionization state across the different galaxy environments. High values of F335M$_{\rm PAH}$/F1130W and low values of F1130W/F770W are measured in H II regions in all three galaxies. This suggests that these regions are populated by hotter PAHs, and/or that the PAH ionization fraction is larger. We see additional evidence of heating and/or changes in PAH size in regions with higher molecular gas content as well as increased ionization in regions with higher H$α$ intensity. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: 12 pages, 6 figures. Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

PHANGS-JWST First Results: A Global and Moderately Resolved View of Mid-Infrared and CO Line Emission from Galaxies at the Start of the JWST Era

Authors: Adam K. Leroy, Alberto D. Bolatto, Karin Sandstrom, Erik Rosolowsky, Ashley. T. Barnes, F. Bigiel, Médéric Boquien, Jakob S. den Brok, Yixian Cao, Jérémy Chastenet, Mélanie Chevance, I-Da Chiang, Ryan Chown, Dario Colombo, Sara L. Ellison, Eric Emsellem, Kathryn Grasha, Jonathan D. Henshaw, Annie Hughes, Ralf S. Klessen, Eric W. Koch, Jaeyeon Kim, Kathryn Kreckel, J. M. Diederik Kruijssen, Kirsten L. Larson , et al. (19 additional authors not shown)

Abstract: We explore the relationship between mid-infrared (mid-IR) and CO rotational line emission from massive star-forming galaxies, which is one of the tightest scalings in the local universe. We assemble a large set of unresolved and moderately ($\sim 1$ kpc) spatially resolved measurements of CO (1-0) and CO (2-1) intensity, $I_{\rm CO}$, and mid-IR intensity, $I_{\rm MIR}$, at 8, 12, 22, and 24$μ$m.… ▽ More We explore the relationship between mid-infrared (mid-IR) and CO rotational line emission from massive star-forming galaxies, which is one of the tightest scalings in the local universe. We assemble a large set of unresolved and moderately ($\sim 1$ kpc) spatially resolved measurements of CO (1-0) and CO (2-1) intensity, $I_{\rm CO}$, and mid-IR intensity, $I_{\rm MIR}$, at 8, 12, 22, and 24$μ$m. The $I_{\rm CO}$ vs. $I_{\rm MIR}$ relationship is reasonably described by a power law with slopes $0.7{-}1.2$ and normalization $I_{\rm CO} \sim 1$ K km s$^{-1}$ at $I_{\rm MIR} \sim 1$ MJy sr$^{-1}$. Both the slopes and intercepts vary systematically with choice of line and band. The comparison between the relations measured for CO~(1-0) and CO (2-1) allow us to infer that $R_{21} \propto I_{\rm MIR}^{0.2}$, in good agreement with other work. The $8μ$m and $12μ$m bands, with strong PAH features, show steeper CO vs. mid-IR slopes than the $22μ$m and $24μ$m, consistent with PAH emission arising not just from CO-bright gas but also from atomic or CO-dark gas. The CO-to-mid-IR ratio correlates with global galaxy stellar mass ($M_\star$) and anti-correlates with SFR/$M_\star$. At $\sim 1$ kpc resolution, the first four PHANGS-JWST targets show CO to mid-IR relationships that are quantitatively similar to our larger literature sample, including showing the steep CO-to-mid-IR slopes for the JWST PAH-tracing bands, although we caution that these initial data have a small sample size and span a limited range of intensities. △ Less

Submitted 27 December, 2022; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 29 pages, 13 figures, key quantitative results in Table 3, Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

Comparing the Locations of Supernovae to CO (2-1) Emission in their Host Galaxies

Authors: Ness Mayker Chen, Adam K. Leroy, Laura A. Lopez, Samantha Benincasa, Mélanie Chevance, Simon C. O. Glover, Annie Hughes, Kathryn Kreckel, Sumit Sarbadhicary, Jiayi Sun, Todd A. Thompson, Dyas Utomo, Frank Bigiel, Guillermo A. Blanc, Daniel A. Dale, Kathryn Grasha, J. M. Diederik Kruijssen, Hsi-An Pan, Miguel Querejeta, Eva Schinnerer, Elizabeth J. Watkins, Thomas G. Williams

Abstract: We measure the molecular gas environment near recent ($< 100$ yr old) supernovae (SNe) using $\sim1''$ or $\leq 150$pc resolution CO (2-1) maps from the PHANGS-ALMA survey of nearby star-forming galaxies. This is arguably the first such study to approach the scales of individual massive molecular clouds ($M_{\rm mol} \gtrsim 10^{5.3}$ M$_{\odot}$). Using the Open Supernova Catalog (OSC), we identi… ▽ More We measure the molecular gas environment near recent ($< 100$ yr old) supernovae (SNe) using $\sim1''$ or $\leq 150$pc resolution CO (2-1) maps from the PHANGS-ALMA survey of nearby star-forming galaxies. This is arguably the first such study to approach the scales of individual massive molecular clouds ($M_{\rm mol} \gtrsim 10^{5.3}$ M$_{\odot}$). Using the Open Supernova Catalog (OSC), we identify 63 SNe within the PHANGS-ALMA footprint. We detect CO (2-1) emission near $\sim60\%$ of the sample at 150pc resolution, compared to $\sim35\%$ of map pixels with CO (2-1) emission, and up to $\sim95\%$ of the SNe at 1kpc resolution compared to $\sim80\%$ of map pixels with CO (2-1) emission. We expect the $\sim60\%$ of SNe within the same 150pc beam as a GMC will likely interact with these clouds in the future, consistent with the observation of widespread SN-molecular gas interaction in the Milky Way, while the other $\sim40\%$ of SNe without strong CO (2-1) detections will deposit their energy in the diffuse interstellar medium (ISM), perhaps helping drive large-scale turbulence or galactic outflows. Broken down by type, we detect CO (2-1) emission at the sites of $\sim85\%$ of our 9 stripped-envelope SNe (SESNe), $\sim40\%$ of our 34 Type II SNe, and $\sim35\%$ of our 13 Type Ia SNe, indicating that SESNe are most closely associated with the brightest CO (2-1) emitting regions in our sample. Our results confirm that SN explosions are not restricted to only the densest gas, and instead exert feedback across a wide range of molecular gas densities. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: Accepted to ApJ; 38 pages, 12 figures, 7 tables in two-column AASTEX63 format

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)

PHANGS-JWST First Results: Stellar Feedback-Driven Excitation and Dissociation of Molecular Gas in the Starburst Ring of NGC 1365?

Authors: Daizhong Liu, Eva Schinnerer, Yixian Cao, Adam Leroy, Antonio Usero, Erik Rosolowsky, Eric Emsellem, J. M. Diederik Kruijssen, Mélanie Chevance, Simon C. O. Glover, Mattia C. Sormani, Alberto D. Bolatto, Jiayi Sun, Sophia K. Stuber, Yu-Hsuan Teng, Frank Bigiel, Ivana Bešlić, Kathryn Grasha, Jonathan D. Henshaw, Ashley. T. Barnes, Jakob S. den Brok, Toshiki Saito, Daniel A. Dale, Elizabeth J. Watkins, Hsi-An Pan , et al. (14 additional authors not shown)

Abstract: We compare embedded young massive star clusters (YMCs) to (sub-)millimeter line observations tracing the excitation and dissociation of molecular gas in the starburst ring of NGC 1365. This galaxy hosts one of the strongest nuclear starbursts and richest populations of YMCs within 20 Mpc. Here we combine near-/mid-IR PHANGS-JWST imaging with new ALMA multi-J CO (1-0, 2-1 and 4-3) and [CI](1-0) map… ▽ More We compare embedded young massive star clusters (YMCs) to (sub-)millimeter line observations tracing the excitation and dissociation of molecular gas in the starburst ring of NGC 1365. This galaxy hosts one of the strongest nuclear starbursts and richest populations of YMCs within 20 Mpc. Here we combine near-/mid-IR PHANGS-JWST imaging with new ALMA multi-J CO (1-0, 2-1 and 4-3) and [CI](1-0) mapping, which we use to trace CO excitation via R42 = I_CO(4-3)/I_CO(2-1) and R21 = I_CO(2-1)/I_CO(1-0) and dissociation via RCICO = I_[CI](1-0)/I_CO(2-1) at 330 pc resolution. We find that the gas flowing into the starburst ring from northeast to southwest appears strongly affected by stellar feedback, showing decreased excitation (lower R42) and increased signatures of dissociation (higher RCICO) in the downstream regions. There, radiative transfer modeling suggests that the molecular gas density decreases and temperature and [CI/CO] abundance ratio increase. We compare R42 and RCICO with local conditions across the regions and find that both correlate with near-IR 2 um emission tracing the YMCs and with both PAH (11.3 um) and dust continuum (21 um) emission. In general, RCICO exhibits ~ 0.1 dex tighter correlations than R42, suggesting CI to be a more sensitive tracer of changing physical conditions in the NGC 1365 starburst than CO (4-3). Our results are consistent with a scenario where gas flows into the two arm regions along the bar, becomes condensed/shocked, forms YMCs, and then these YMCs heat and dissociate the gas. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: 19 pages, 7 figures and 2 tables in total (12 pages and 6 figures in main text). Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJL

PHANGS-JWST First Results: Rapid Evolution of Star Formation in the Central Molecular Gas Ring of NGC1365

Authors: Eva Schinnerer, Eric Emsellem, Jonathan D. Henshaw, Daizhong Liu, Sharon E. Meidt, Miguel Querejeta, Florent Renaud, Mattia C. Sormani, Jiayi Sun, Oleg V. Egorov, Kirsten L. Larson, Adam K. Leroy, Erik Rosolowsky, Karin M. Sandstrom, T. G. Williams, Ashley T. Barnes, F. Bigiel, Melanie Chevance, Yixian Cao, Rupali Chandar, Daniel A. Dale, Cosima Eibensteiner, Simon C. O. Glover, Kathryn Grasha, Stephen Hannon , et al. (14 additional authors not shown)

Abstract: Large-scale bars can fuel galaxy centers with molecular gas, often leading to the development of dense ring-like structures where intense star formation occurs, forming a very different environment compared to galactic disks. We pair ~0.3" (30pc) resolution new JWST/MIRI imaging with archival ALMA CO(2-1) mapping of the central ~5kpc of the nearby barred spiral galaxy NGC1365, to investigate the p… ▽ More Large-scale bars can fuel galaxy centers with molecular gas, often leading to the development of dense ring-like structures where intense star formation occurs, forming a very different environment compared to galactic disks. We pair ~0.3" (30pc) resolution new JWST/MIRI imaging with archival ALMA CO(2-1) mapping of the central ~5kpc of the nearby barred spiral galaxy NGC1365, to investigate the physical mechanisms responsible for this extreme star formation. The molecular gas morphology is resolved into two well-known bright bar lanes that surround a smooth dynamically cold gas disk (R_gal ~ 475pc) reminiscent of non-star-forming disks in early type galaxies and likely fed by gas inflow triggered by stellar feedback in the lanes. The lanes host a large number of JWST-identified massive young star clusters. We find some evidence for temporal star formation evolution along the ring. The complex kinematics in the gas lanes reveal strong streaming motions and may be consistent with convergence of gas streamlines expected there. Indeed, the extreme line-widths are found to be the result of inter-`cloud' motion between gas peaks; ScousePy decomposition reveals multiple components with line widths of <sigma_CO,scouse> ~ 19km/s and surface densities of <Sigma_H2,scouse> ~ 800M_sun/pc^2, similar to the properties observed throughout the rest of the central molecular gas structure. Tailored hydro-dynamical simulations exhibit many of the observed properties and imply that the observed structures are transient and highly time-variable. From our study of NGC1365, we conclude that it is predominantly the high gas inflow triggered by the bar that is setting the star formation in its CMZ. △ Less

Submitted 18 December, 2022; originally announced December 2022.

Comments: 24 pages, 8 figures, accepted for publications as part of PHANGS-JWST ApJL Focus Issue

PHANGS-JWST First Results: Destruction of the PAH molecules in HII regions probed by JWST and MUSE

Authors: Oleg V. Egorov, Kathryn Kreckel, Karin M. Sandstrom, Adam K. Leroy, Simon C. O. Glover, Brent Groves, J. M. Diederik Kruijssen, Ashley. T. Barnes, Francesco Belfiore, F. Bigiel, Guillermo A. Blanc, Médéric Boquien, Yixian Cao, Jérémy Chastenet, Mélanie Chevance, Enrico Congiu, Daniel A. Dale, Eric Emsellem, Kathryn Grasha, Ralf S. Klessen, Kirsten L. Larson, Daizhong Liu, Eric J. Murphy, Hsi-An Pan, Ismael Pessa , et al. (8 additional authors not shown)

Abstract: Polycyclic aromatic hydrocarbons (PAHs) play a critical role in the reprocessing of stellar radiation and in balancing the heating and cooling processes in the interstellar medium (ISM), but appear to be destroyed in HII regions. However, the mechanisms driving their destruction are still not completely understood. Using PHANGS-JWST and PHANGS-MUSE observations, we investigate how the PAH fraction… ▽ More Polycyclic aromatic hydrocarbons (PAHs) play a critical role in the reprocessing of stellar radiation and in balancing the heating and cooling processes in the interstellar medium (ISM), but appear to be destroyed in HII regions. However, the mechanisms driving their destruction are still not completely understood. Using PHANGS-JWST and PHANGS-MUSE observations, we investigate how the PAH fraction changes in about 1500 HII regions across four nearby star-forming galaxies (NGC 628, NGC 1365, NGC 7496, IC 5332). We find a strong anti-correlation between the PAH fraction and the ionization parameter (the ratio between the ionizing photon flux and the hydrogen density) of HII regions. This relation becomes steeper for more luminous HII regions. The metallicity of HII regions has only a minor impact on these results in our galaxy sample. We find that the PAH fraction decreases with the H$α$ equivalent width - a proxy for the age of the HII regions - although this trend is much weaker than the one identified using the ionization parameter. Our results are consistent with a scenario where hydrogen-ionizing UV radiation is the dominant source of PAH destruction in star-forming regions. △ Less

Submitted 18 December, 2022; originally announced December 2022.

Comments: 16 pages, 7 figures, 1 table. Accepted for publication in ApJL as part of a PHANGS-JWST First Results Focus issue

PHANGS--JWST First Results: ISM structure on the turbulent Jeans scale in four disk galaxies observed by JWST and ALMA

Authors: Sharon E. Meidt, Erik Rosolowsky, Jiayi Sun, Eric W. Koch, Ralf S. Klessen, Adam K. Leroy, Eva Schinnerer, Ashley. T. Barnes, Simon C. O. Glover, Janice C. Lee, Arjen van der Wel, Elizabeth J. Watkins, Thomas G. Williams, Frank Bigiel, Médéric Boquien, Guillermo A. Blanc, Yixian Cao, Mélanie Chevance, Daniel A. Dale, Oleg V. Egorov, Eric Emsellem, Kathryn Grasha, Jonathan D. Henshaw, J. M. Diederik Kruijssen, Kirsten L. Larson , et al. (9 additional authors not shown)

Abstract: JWST/MIRI imaging of the nearby galaxies IC 5332, NGC 628, NGC 1365 and NGC 7496 from PHANGS reveals a richness of gas structures that in each case form a quasi-regular network of interconnected filaments, shells and voids. We examine whether this multi-scale network of structure is consistent with the fragmentation of the gas disk through gravitational instability. We use FilFinder to detect the… ▽ More JWST/MIRI imaging of the nearby galaxies IC 5332, NGC 628, NGC 1365 and NGC 7496 from PHANGS reveals a richness of gas structures that in each case form a quasi-regular network of interconnected filaments, shells and voids. We examine whether this multi-scale network of structure is consistent with the fragmentation of the gas disk through gravitational instability. We use FilFinder to detect the web of filamentary features in each galaxy and determine their characteristic radial and azimuthal spacings. These spacings are then compared to estimates of the most Toomre-unstable length (a few kpc), the turbulent Jeans length (a few hundred pc) and the disk scale height (tens of pc) reconstructed using PHANGS-ALMA observations of the molecular gas as a dynamical tracer. Our analysis of the four galaxies targeted in this work indicates that Jeans-scale structure is pervasive. Future work will be essential for determining how the structure observed in gas disks impacts not only the rate and location of star formation but also how stellar feedback interacts positively or negatively with the surrounding multi-phase gas reservoir. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: 15 pages, 4 figures. Accepted as part of a PHANGS-JWST ApJL Focus Issue

The PHANGS-JWST Treasury Survey: Star Formation, Feedback, and Dust Physics at High Angular resolution in Nearby GalaxieS

Authors: Janice C. Lee, Karin M. Sandstrom, Adam K. Leroy, David A. Thilker, Eva Schinnerer, Erik Rosolowsky, Kirsten L. Larson, Oleg V. Egorov, Thomas G. Williams, Judy Schmidt, Eric Emsellem, Gagandeep S. Anand, Ashley T. Barnes, Francesco Belfiore, Ivana Beslic, Frank Bigiel, Guillermo A. Blanc, Alberto D. Bolatto, Mederic Boquien, Jakob den Brok, Yixian Cao, Rupali Chandar, Jeremy Chastenet, Melanie Chevance, I-Da Chiang , et al. (52 additional authors not shown)

Abstract: The PHANGS collaboration has been building a reference dataset for the multi-scale, multi-phase study of star formation and the interstellar medium in nearby galaxies. With the successful launch and commissioning of JWST, we can now obtain high-resolution infrared imaging to probe the youngest stellar populations and dust emission on the scales of star clusters and molecular clouds ($\sim$5-50 pc)… ▽ More The PHANGS collaboration has been building a reference dataset for the multi-scale, multi-phase study of star formation and the interstellar medium in nearby galaxies. With the successful launch and commissioning of JWST, we can now obtain high-resolution infrared imaging to probe the youngest stellar populations and dust emission on the scales of star clusters and molecular clouds ($\sim$5-50 pc). In Cycle 1, PHANGS is conducting an 8-band imaging survey from 2-21$μ$m of 19 nearby spiral galaxies. CO(2-1) mapping, optical integral field spectroscopy, and UV-optical imaging for all 19 galaxies have been obtained through large programs with ALMA, VLT/MUSE, and Hubble. PHANGS-JWST enables a full inventory of star formation, accurate measurement of the mass and age of star clusters, identification of the youngest embedded stellar populations, and characterization of the physical state of small dust grains. When combined with Hubble catalogs of $\sim$10,000 star clusters, MUSE spectroscopic mapping of $\sim$20,000 HII regions, and $\sim$12,000 ALMA-identified molecular clouds, it becomes possible to measure the timescales and efficiencies of the earliest phases of star formation and feedback, build an empirical model of the dependence of small dust grain properties on local ISM conditions, and test our understanding of how dust-reprocessed starlight traces star formation activity, all across a diversity of galactic environments. Here we describe the PHANGS-JWST Treasury survey, present the remarkable imaging obtained in the first few months of science operations, and provide context for the initial results presented in the first series of PHANGS-JWST publications. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: Re-submitted after addressing minor comments from referee. To be published as part of PHANGS-JWST ApJL Focus Issue

PHANGS-JWST First Results: The 21 $μ$m Compact Source Population

Authors: Hamid Hassani, Erik Rosolowsky, Adam K. Leroy, Mederic Boquien, Janice C. Lee, Ashley. T. Barnes, Francesco Belfiore, Frank Bigiel, Yixian Cao, Melanie Chevance, Daniel A. Dale, Oleg V. Egorov, Eric Emsellem, Christopher M. Faesi, Kathryn Grasha, Jaeyeon Kim, Ralf S. Klessen, Kathryn Kreckel, J. M. Diederik Kruijssen, Kirsten L. Larson, Sharon E. Meidt, Karin M. Sandstrom, Eva Schinnerer, David A. Thilker, Elizabeth J. Watkins , et al. (2 additional authors not shown)

Abstract: We use PHANGS-JWST data to identify and classify 1271 compact 21 $μ$m sources in four nearby galaxies using MIRI F2100W data. We identify sources using a dendrogram-based algorithm, and we measure the background-subtracted flux densities for JWST bands from 2 $μ$m to 21 $μ$m. Using the SED in JWST as well as HST bands, plus ALMA and MUSE/VLT observations, we classify the sources by eye. Then we us… ▽ More We use PHANGS-JWST data to identify and classify 1271 compact 21 $μ$m sources in four nearby galaxies using MIRI F2100W data. We identify sources using a dendrogram-based algorithm, and we measure the background-subtracted flux densities for JWST bands from 2 $μ$m to 21 $μ$m. Using the SED in JWST as well as HST bands, plus ALMA and MUSE/VLT observations, we classify the sources by eye. Then we use this classification to define regions in color-color space, and so establish a quantitative framework for classifying sources. We identify 1085 sources as belonging to the ISM of the target galaxies with the remainder being dusty stars or background galaxies. These 21 $μ$m sources are strongly spatially associated with HII regions ($>92\%$ of sources), while 74$\%$ of sources are coincident with a stellar association defined in the HST data. Using SED fitting, we find that the stellar masses of the 21 $μ$m sources span a range of 10$^{2}$ to 10$^{4}~M_\odot$ with mass-weighted ages down to 2 Myr. There is a tight correlation between attenuation-corrected H$α$ and 21 $μ$m luminosity for $L_{ν,\mathrm{F2100W}}>10^{19}~\mathrm{W~Hz}^{-1}$. Young embedded source candidates selected at 21 $μ$m are found below this threshold and have $M_\star < 10^{3}~M_\odot$. △ Less

Submitted 2 December, 2022; originally announced December 2022.

Comments: 15 pages, 5 figures, Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

PHANGS-JWST First Results: Multi-wavelength view of feedback-driven bubbles (The Phantom Voids) across NGC 628

Authors: Ashley T. Barnes, Elizabeth J. Watkins, Sharon E. Meidt, Kathryn Kreckel, Mattia C. Sormani, Robin G. Tress, Simon C. O. Glover, Frank Bigiel, Rupali Chandar, Eric Emsellem, Janice C. Lee, Adam K. Leroy, Karin M. Sandstrom, Eva Schinnerer, Erik W. Rosolowsky, Francesco Belfiore, Guillermo Blanc, Mederic Boquien, Jakob S. den Brok, Yixian Cao, Mélanie Chevance, Daniel A. Dale, Oleg Egorov, Cosima Eibensteiner, Kathryn Grasha , et al. (29 additional authors not shown)

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

Submitted 1 December, 2022; originally announced December 2022.

Comments: 12 pages total, 8 figures, and 1 table. Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

PHANGS-JWST First Results: A statistical view on bubble evolution in NGC628

Authors: Elizabeth J. Watkins, Ashley Barnes, Kiana F. Henny, Hwihyun Kim, Kathryn Kreckel, Sharon E. Meidt, Ralf S. Klessen, Simon C. O. Glover, Thomas G. Williams, B. W. Keller, Adam K. Leroy, Erik W. Rosolowsky, Mederic Boquien, Gagandeep S. Anand, Francesco Belfiore, Frank Bigiel, Guillermo Blanc, Yixian Cao, Rupali Chandar, Ness Mayker Chen, Mélanie Chevance, Enrico Congiu, Daniel A. Dale, Sinan Deger, Oleg Egorov , et al. (27 additional authors not shown)

Abstract: The first JWST observations of nearby galaxies have unveiled a rich population of bubbles that trace the stellar feedback mechanisms responsible for their creation. Studying these bubbles therefore allows us to chart the interaction between stellar feedback and the interstellar medium, and the larger galactic flows needed to regulate star formation processes globally. We present the first catalog… ▽ More The first JWST observations of nearby galaxies have unveiled a rich population of bubbles that trace the stellar feedback mechanisms responsible for their creation. Studying these bubbles therefore allows us to chart the interaction between stellar feedback and the interstellar medium, and the larger galactic flows needed to regulate star formation processes globally. We present the first catalog of bubbles in NGC628, visually identified using MIRI F770W PHANGS-JWST observations, and use them to statistically evaluate bubble characteristics. We classify 1694 structures as bubbles with radii between 6-552 pc. Of these, 31% contain at least one smaller bubble at their edge, indicating that previous generations of star formation have a local impact on where new stars form. On large scales, most bubbles lie near a spiral arm, and their radii increase downstream compared to upstream. Furthermore, bubbles are elongated in a similar direction to the spiral arm ridge-line. These azimuthal trends demonstrate that star formation is intimately connected to the spiral arm passage. Finally, the bubble size distribution follows a power-law of index $p=-2.2\pm0.1$, which is slightly shallower than the theoretical value by 1-3.5$σ$ that did not include bubble mergers. The fraction of bubbles identified within the shells of larger bubbles suggests that bubble merging is a common process. Our analysis therefore allows us to quantify the number of star-forming regions that are influenced by an earlier generation, and the role feedback processes have in setting the global star formation rate. With the full PHANGS-JWST sample, we can do this for more galaxies. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Comments: 25 pages total, 13 Figures and 1 Table. Accepted for publication in ApJL as part of a PHANGS-JWST First Results Focus issue

PHANGS-JWST First Results: The Influence of Stellar Clusters on PAHs in Nearby Galaxies

Authors: Daniel A. Dale, Médéric Boquien, Ashley T. Barnes, Francesco Belfiore, F. Bigiel, Yixian Cao, Rupali Chandar, Jérémy Chastenet, Mélanie Chevance, Sinan Deger, Oleg V. Egorov, Kathryn Grasha, Brent Groves, Hamid Hassani, Kiana F. Henny, Ralf S. Klessen, Kathryn Kreckel, J. M. Diederik Kruijssen, Kirsten L. Larson, Janice C. Lee, Adam K. Leroy, Daizhong Liu, Eric J. Murphy, Erik Rosolowsky, Karin Sandstrom , et al. (6 additional authors not shown)

Abstract: We present a comparison of theoretical predictions of dust continuum and polycyclic aromatic hydrocarbon (PAH) emission with new JWST observations in three nearby galaxies: NGC 628, NGC 1365, and NGC 7496. Our analysis focuses on a total of 1063 compact stellar clusters and 2654 stellar associations previously characterized by HST in the three galaxies. We find that the distributions and trends in… ▽ More We present a comparison of theoretical predictions of dust continuum and polycyclic aromatic hydrocarbon (PAH) emission with new JWST observations in three nearby galaxies: NGC 628, NGC 1365, and NGC 7496. Our analysis focuses on a total of 1063 compact stellar clusters and 2654 stellar associations previously characterized by HST in the three galaxies. We find that the distributions and trends in the observed PAH-focused infrared colors generally agree with theoretical expectations, and that the bulk of the observations is more aligned with models of larger, ionized PAHs. These JWST data usher in a new era of probing interstellar dust and studying how the intense radiation fields near stellar clusters and associations play a role in shaping the physical properties of PAHs. △ Less

Submitted 30 November, 2022; originally announced December 2022.

Comments: 8 pages, accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

PHANGS-JWST First Results: Spurring on Star Formation: JWST Reveals Localised Star Formation in a Spiral Arm Spur of NGC 628

Authors: Thomas G. Williams, Jiayi Sun, Ashley T. Barnes, Eva Schinnerer, Jonathan D. Henshaw, Sharon E. Meidt, Miguel Querejeta, Elizabeth J. Watkins, Frank Bigiel, Guillermo A. Blanc, Médéric Boquien, Yixian Cao, Mélanie Chevance, Oleg V. Egorov, Eric Emsellem, Simon C. O. Glover, Kathryn Grasha, Hamid Hassani, Sarah Jeffreson, María J. Jiménez-Donaire, Jaeyeon Kim, Ralf S. Klessen, Kathryn Kreckel, J. M. Diederik Kruijssen, Kirsten L. Larson , et al. (12 additional authors not shown)

Abstract: We combine JWST observations with ALMA CO and VLT-MUSE H$α$ data to examine off-spiral arm star formation in the face-on, grand-design spiral galaxy NGC 628. We focus on the northern spiral arm, around a galactocentric radius of 3-4 kpc, and study two spurs. These form an interesting contrast, as one is CO-rich and one CO-poor, and they have a maximum azimuthal offset in MIRI 21$μ$m and MUSE H$α$… ▽ More We combine JWST observations with ALMA CO and VLT-MUSE H$α$ data to examine off-spiral arm star formation in the face-on, grand-design spiral galaxy NGC 628. We focus on the northern spiral arm, around a galactocentric radius of 3-4 kpc, and study two spurs. These form an interesting contrast, as one is CO-rich and one CO-poor, and they have a maximum azimuthal offset in MIRI 21$μ$m and MUSE H$α$ of around 40$^\circ$ (CO-rich) and 55$^\circ$ (CO-poor) from the spiral arm. The star formation rate is higher in the regions of the spurs near to spiral arms, but the star formation efficiency appears relatively constant. Given the spiral pattern speed and rotation curve of this galaxy and assuming material exiting the arms undergoes purely circular motion, these offsets would be reached in 100-150 Myr, significantly longer than the 21$μ$m and H$α$ star formation timescales (both <10 Myr). The invariance of the star formation efficiency in the spurs versus the spiral arms indicates massive star formation is not only triggered in spiral arms, and cannot simply occur in the arms and then drift away from the wave pattern. These early JWST results show that in-situ star formation likely occurs in the spurs, and that the observed young stars are not simply the `leftovers' of stellar birth in the spiral arms. The excellent physical resolution and sensitivity that JWST can attain in nearby galaxies will well resolve individual star-forming regions and help us to better understand the earliest phases of star formation. △ Less

Submitted 2 March, 2023; v1 submitted 30 November, 2022; originally announced December 2022.

Comments: 11 pages, 5 Figures, accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ. Updated to include missing author and published paper references

PHANGS-JWST First Results: Duration of the early phase of massive star formation in NGC628

Authors: Jaeyeon Kim, Mélanie Chevance, J. M. Diederik Kruijssen, Ashley. T. Barnes, Frank Bigiel, Guillermo A. Blanc, Médéric Boquien, Yixian Cao, Enrico Congiu, Daniel A. Dale, Oleg V. Egorov, Christopher M. Faesi, Simon C. O. Glover, Kathryn Grasha, Brent Groves, Hamid Hassani, Annie Hughes, Ralf S. Klessen, Kathryn Kreckel, Kirsten L. Larson, Janice C. Lee, Adam K. Leroy, Daizhong Liu, Steven N. Longmore, Sharon E. Meidt , et al. (11 additional authors not shown)

Abstract: The earliest stages of star formation, when young stars are still deeply embedded in their natal clouds, represent a critical phase in the matter cycle between gas clouds and young stellar regions. Until now, the high-resolution infrared observations required for characterizing this heavily obscured phase (during which massive stars have formed, but optical emission is not detected) could only be… ▽ More The earliest stages of star formation, when young stars are still deeply embedded in their natal clouds, represent a critical phase in the matter cycle between gas clouds and young stellar regions. Until now, the high-resolution infrared observations required for characterizing this heavily obscured phase (during which massive stars have formed, but optical emission is not detected) could only be obtained for a handful of the most nearby galaxies. One of the main hurdles has been the limited angular resolution of the Spitzer Space Telescope. With the revolutionary capabilities of the JWST, it is now possible to investigate the matter cycle during the earliest phases of star formation as a function of the galactic environment. In this Letter, we demonstrate this by measuring the duration of the embedded phase of star formation and the implied time over which molecular clouds remain inert in the galaxy NGC628 at a distance of 9.8Mpc, demonstrating that the cosmic volume where this measurement can be made has increased by a factor of $>100$ compared to Spitzer. We show that young massive stars remain embedded for $5.1_{-1.4}^{+2.7}$Myr ($2.3_{-1.4}^{+2.7}$Myr of which being heavily obscured), representing $\sim20\%$ of the total cloud lifetime. These values are in broad agreement with previous measurements in five nearby ($D < 3.5$Mpc) galaxies and constitute a proof of concept for the systematic characterization of the early phase of star formation across the nearby galaxy population with the PHANGS-JWST survey. △ Less

Submitted 5 December, 2022; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: 15 pages, 4 figures, 2 tables; accepted for publication in ApJL as part of a PHANGS-JWST First Results Focus issue

PHANGS-JWST First Results: A combined HST and JWST analysis of the nuclear star cluster in NGC 628

Authors: Nils Hoyer, Francesca Pinna, Albrecht W. H. Kamlah, Francisco Nogueras-Lara, Anja Feldmeier-Krause, Nadine Neumayer, Mattia C. Sormani, Médéric Boquien, Eric Emsellem, Anil C. Seth, Ralf S. Klessen, Thomas G. Williams, Eva Schinnerer, Ashley T. Barnes, Adam K. Leroy, Silvia Bonoli, J. M. Diederik Kruijssen, Justus Neumann, Patricia Sánchez-Blázquez, Daniel A. Dale, Elizabeth J. Watkins, David A. Thilker, Erik Rosolowsky, Frank Bigiel, Kathryn Grasha , et al. (6 additional authors not shown)

Abstract: We combine archival HST and new JWST imaging data, covering the ultraviolet to mid-infrared regime, to morphologically analyze the nuclear star cluster (NSC) of NGC 628, a grand-design spiral galaxy. The cluster is located in a 200 pc x 400 pc cavity, lacking both dust and gas. We find roughly constant values for the effective radius (r_eff ~ 5 pc) and ellipticity (ε ~ 0.05), while the Sérsic inde… ▽ More We combine archival HST and new JWST imaging data, covering the ultraviolet to mid-infrared regime, to morphologically analyze the nuclear star cluster (NSC) of NGC 628, a grand-design spiral galaxy. The cluster is located in a 200 pc x 400 pc cavity, lacking both dust and gas. We find roughly constant values for the effective radius (r_eff ~ 5 pc) and ellipticity (ε ~ 0.05), while the Sérsic index (n) and position angle (PA) drop from n ~ 3 to ~ 2 and PA ~ 130° to 90°, respectively. In the mid-infrared, r_eff ~ 12pc, ε ~ 0.4, and n ~ 1-1.5, with the same PA ~ 90°. The NSC has a stellar mass of log10 (M_nsc / M_Sun) = 7.06 +- 0.31, as derived through B-V, confirmed when using multi-wavelength data, and in agreement with the literature value. Fitting the spectral energy distribution, excluding the mid-infrared data, yields a main stellar population's age of (8 +- 3) Gyr with a metallicity of Z = 0.012 +- 0.006. There is no indication of any significant star formation over the last few Gyr. Whether gas and dust were dynamically kept out or evacuated from the central cavity remains unclear. The best-fit suggests an excess of flux in the mid-infrared bands, with further indications that the center of the mid-infrared structure is displaced with respect to the optical center of the NSC. We discuss five potential scenarios, none of them fully explaining both the observed photometry and structure. △ Less

Submitted 10 January, 2023; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: 26 pages, 10 figures, 6 tables. Accepted for publication by ApJL

Globular cluster metallicity distributions in the E-MOSAICS simulations

Authors: Joel Pfeffer, J. M. Diederik Kruijssen, Nate Bastian, Robert A. Crain, Sebastian Trujillo-Gomez

Abstract: The metallicity distributions of globular cluster (GC) systems in galaxies are a critical test of any GC formation scenario. In this work, we investigate the predicted GC metallicity distributions of galaxies in the MOdelling Star cluster population Assembly In Cosmological Simulations within EAGLE (E-MOSAICS) simulation of a representative cosmological volume ($L = 34.4$ comoving Mpc). We find th… ▽ More The metallicity distributions of globular cluster (GC) systems in galaxies are a critical test of any GC formation scenario. In this work, we investigate the predicted GC metallicity distributions of galaxies in the MOdelling Star cluster population Assembly In Cosmological Simulations within EAGLE (E-MOSAICS) simulation of a representative cosmological volume ($L = 34.4$ comoving Mpc). We find that the predicted GC metallicity distributions and median metallicities from the fiducial E-MOSAICS GC formation model agree well the observed distributions, except for galaxies with masses $M_\ast \sim 2 \times 10^{10}$ M$_\odot$, which contain an overabundance of metal-rich GCs. The predicted fraction of galaxies with bimodal GC metallicity distributions ($37 \pm 2$ per cent in total; $45 \pm 7$ per cent for $M_\ast > 10^{10.5}$ M$_\odot$) is in good agreement with observed fractions ($44^{+10}_{-9}$ per cent), as are the mean metallicities of the metal-poor and metal-rich peaks. We show that, for massive galaxies ($M_\ast > 10^{10}$ M$_\odot$), bimodal GC distributions primarily occur as a result of cluster disruption from initially-unimodal distributions, rather than as a result of cluster formation processes. Based on the distribution of field stars with GC-like abundances in the Milky Way, we suggest that the bimodal GC metallicity distribution of Milky Way GCs also occurred as a result of cluster disruption, rather than formation processes. We conclude that separate formation processes are not required to explain metal-poor and metal-rich GCs, and that GCs can be considered as the surviving analogues of young massive star clusters that are readily observed to form in the local Universe today. △ Less

Submitted 24 January, 2023; v1 submitted 1 October, 2022; originally announced October 2022.

Comments: 19 pages, 16 figures. Published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 519, Issue 4, March 2023, Pages 5384-5401

Towards a multi-tracer timeline of star formation in the LMC -- II. The formation and destruction of molecular clouds

Authors: Jacob L. Ward, J. M. Diederik Kruijssen, Mélanie Chevance, Jaeyeon Kim, Steven N. Longmore

Abstract: The time-scales associated with various stages of the star formation process represent major unknowns in our understanding of galactic evolution, as well as of star and planet formation. This is the second paper in a series aiming to establish a multi-tracer time-line of star formation in the Large Magellanic Cloud (LMC), focusing on the lifecycle of molecular clouds. We use a statistical method t… ▽ More The time-scales associated with various stages of the star formation process represent major unknowns in our understanding of galactic evolution, as well as of star and planet formation. This is the second paper in a series aiming to establish a multi-tracer time-line of star formation in the Large Magellanic Cloud (LMC), focusing on the lifecycle of molecular clouds. We use a statistical method to determine a molecular cloud lifetime in the LMC of $t_{\text{CO}}=11.8^{+2.7}_{-2.2}$ Myr. This short time-scale is similar to the cloud dynamical time, and suggests that molecular clouds in the LMC are largely decoupled from the effects of galactic dynamics and have lifetimes set by internal processes. This provides a clear contrast to atomic clouds in the LMC, of which the lifetimes are correlated with galactic dynamical time-scales. We additionally derive the time-scale for which molecular clouds and HII regions co-exist as $t_{\text{fb}}=1.2^{+0.3}_{-0.2}$ Myr, implying an average feedback front expansion velocity of 12 km s$^{-1}$, consistent with expansion velocities of HII regions in the LMC observed directly using optical spectroscopy. Taken together, these results imply that the molecular cloud lifecycle in the LMC proceeds rapidly and is regulated by internal dynamics and stellar feedback. We conclude by discussing our measurements in the context of previous work in the literature, which reported considerably longer lifetimes for molecular clouds in the LMC, and find that these previous findings resulted from a subjective choice in timeline calibration that is avoided by our statistical methodology. △ Less

Submitted 12 September, 2022; originally announced September 2022.

Comments: 18 pages (including appendices), 10 figures, 2 tables; MNRAS in press (accepted August 15, 2022). Figures 3, 4, 5, and 7 show the main results of the paper

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

Monochromatic globular clusters as a critical test of formation models for the dark matter deficient galaxies NGC1052-DF2 and NGC1052-DF4

Authors: Pieter van Dokkum, Zili Shen, Aaron J. Romanowsky, Roberto Abraham, Charlie Conroy, Shany Danieli, Dhruba Dutta Chowdhury, Michael A. Keim, J. M. Diederik Kruijssen, Joel Leja, Sebastian Trujillo-Gomez

Abstract: It was recently proposed that the dark matter-deficient ultra-diffuse galaxies DF2 and DF4 in the NGC1052 group could be the products of a "bullet dwarf" collision between two gas-rich progenitor galaxies. In this model DF2 and DF4 formed at the same time in the immediate aftermath of the collision, and a strong prediction is that their globular clusters should have nearly identical stellar popula… ▽ More It was recently proposed that the dark matter-deficient ultra-diffuse galaxies DF2 and DF4 in the NGC1052 group could be the products of a "bullet dwarf" collision between two gas-rich progenitor galaxies. In this model DF2 and DF4 formed at the same time in the immediate aftermath of the collision, and a strong prediction is that their globular clusters should have nearly identical stellar populations. Here we test this prediction by measuring accurate F606W-F814W colors from deep HST/ACS imaging. We find that the clusters are extremely homogeneous. The mean color difference between the globular clusters in DF2 and DF4 is $-0.003\pm 0.005$ mag and the observed scatter for the combined sample of 18 clusters with $M_V<-8.6$ in both galaxies is $0.015 \pm 0.002$ mag. After accounting for observational uncertainties and stochastic cluster-to-cluster variation in the number of red giants, the remaining scatter is $0.008^{+0.005}_{-0.006}$ mag. Both the color difference and the scatter are an order of magnitude smaller than in other dwarf galaxies, and we infer that the bullet scenario passes an important test that could have falsified it. No other formation models have predicted this extreme uniformity of the globular clusters in the two galaxies. We find that the galaxies themselves are slightly redder than the clusters, consistent with a previously-measured metallicity difference. Numerical simulations have shown that such differences are expected in the bullet scenario, as the galaxies continued to self-enrich after the formation of the globular clusters. △ Less

Submitted 21 December, 2022; v1 submitted 14 July, 2022; originally announced July 2022.

Comments: Published in ApJ Letters. It is difficult to visualize the uniformity of these strange clusters. Our best attempt is Fig. 2, where we compare them to globular clusters in Virgo galaxies

Journal ref: ApJ, 940, L9 (2022)