Showing 1–50 of 148 results for author: Meidt, S

Neutral atomic and molecular gas dynamics in the nearby spiral galaxies NGC 1512, NGC 4535, and NGC 7496

Authors: Sebastian Laudage, Cosima Eibensteiner, Frank Bigiel, Adam K. Leroy, Sharon Meidt, Eva Schinnerer, W. J. G. de Blok, Miguele Querejeta, Sophia Stuber, Dario Colombo, Erik Rosolowsky, D. J. Pisano, Dyas Utomo, Rebecca C. Levy, Ralf Klessen, Yixian Cao, Eric W. Koch, Sushma Kurapati, Patricia Sanchez-Blazquez, Justus Neumann, Lukas Neumann, Hsi-An Pan, Thomas G. Williams

Abstract: Neutral atomic gas (HI) effectively traces galactic dynamics across mid to large galactocentric radii. However, its limitations in observing small-scale changes within the central few kiloparsecs, coupled with the often observed HI deficit in galactic centers, necessitates using molecular gas emission as a preferred tracer in these regions. Understanding the dynamics of both neutral atomic and mol… ▽ More Neutral atomic gas (HI) effectively traces galactic dynamics across mid to large galactocentric radii. However, its limitations in observing small-scale changes within the central few kiloparsecs, coupled with the often observed HI deficit in galactic centers, necessitates using molecular gas emission as a preferred tracer in these regions. Understanding the dynamics of both neutral atomic and molecular gas is crucial for a more complete understanding of how galaxies evolve, funnel gas from the outer disk into their central parts, and eventually form stars. In this work we aim to quantify the dynamics of both, the neutral atomic and molecular gas, in the nearby spiral galaxies NGC 1512, NGC 4535, and NGC 7496 using new MeerKAT-HI observations together with ALMA CO (2-1) observations from the PHANGS collaboration. We use the analysis tool 3D-Barolo to fit tilted ring models to the HI and CO observations. A combined approach of using the HI to constrain the true disk orientation parameters before applying these to the CO datasets is tested. This paper sets expectations for the results of the upcoming high-resolution HI coverage of many galaxies in the PHANGS-ALMA sample using MeerKAT or VLA, to establish a robust methodology for characterizing galaxy orientations and deriving dynamics from combining new HI with existing CO data. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: accepted for publication in A&A; 13 pages, 9 Figures (+2 appendix pages)

A 260 pc resolution ALMA map of HCN(1-0) in the galaxy NGC 4321

Authors: Lukas Neumann, Frank Bigiel, Ashley T. Barnes, Molly J. Gallagher, Adam Leroy, Antonio Usero, Erik Rosolowsky, Ivana Bešlić, Médéric Boquien, Yixian Cao, Mélanie Chevance, Dario Colombo, Daniel A. Dale, Cosima Eibensteiner, Kathryn Grasha, Jonathan D. Henshaw, María J. Jiménez-Donaire, Sharon Meidt, Shyam H. Menon, Eric J. Murphy, Hsi-An Pan, Miguel Querejeta, Toshiki Saito, Eva Schinnerer, Sophia K. Stuber , et al. (2 additional authors not shown)

Abstract: The star formation rate (SFR) is tightly connected to the amount of dense gas in molecular clouds. However, it is not fully understood how the relationship between dense molecular gas and star formation varies within galaxies and in different morphological environments. In this work, we study dense gas and star formation in the nearby spiral galaxy NGC 4321 to test how the amount of dense gas and… ▽ More The star formation rate (SFR) is tightly connected to the amount of dense gas in molecular clouds. However, it is not fully understood how the relationship between dense molecular gas and star formation varies within galaxies and in different morphological environments. In this work, we study dense gas and star formation in the nearby spiral galaxy NGC 4321 to test how the amount of dense gas and its ability to form stars varies with environmental properties at 260 pc scales. We present new ALMA observations of HCN(1-0) line emission. Combined with existing CO(2-1) observations from ALMA, and H-alpha from MUSE, as well as F2100W from JWST to trace the SFR, we measure the HCN/CO line ratio, a proxy for the dense gas fraction and SFR/HCN, a proxy for the star formation efficiency of the dense gas. Towards the centre of the galaxy, HCN/CO systematically increases while SFR/HCN decreases, but these ratios stay roughly constant throughout the disc. Spiral arms, interarm regions, and bar ends show similar HCN/CO and SFR/HCN. On the bar, there is a significantly lower SFR/HCN at a similar HCN/CO. We conclude that the centres of galaxies show the strongest environmental influence on dense gas and star formation, suggesting either that clouds couple strongly to the surrounding pressure or that HCN is tracing more of the bulk molecular gas that is less efficiently converted into stars. On the contrary, across the disc of NGC 4321, where the ISM pressure is typically low, SFR/HCN does not show large variations (< 0.3 dex) in agreement with Galactic observations of molecular clouds. Despite the large variations across environments and physical conditions, HCN/CO is a good predictor of the mean molecular gas surface density at 260 pc scales. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 18 pages, 9 figures, accepted for pub in A&A, Jun 13, 2024

Do spiral arms enhance star formation efficiency?

Authors: Miguel Querejeta, Adam K. Leroy, Sharon E. Meidt, Eva Schinnerer, Francesco Belfiore, Eric Emsellem, Ralf S. Klessen, Jiayi Sun, Mattia Sormani, Ivana Bešlic, Yixian Cao, Mélanie Chevance, Dario Colombo, Daniel A. Dale, Santiago García-Burillo, Simon C. O. Glover, Kathryn Grasha, Brent Groves, Eric. W. Koch, Lukas Neumann, Hsi-An Pan, Ismael Pessa, Jérôme Pety, Francesca Pinna, Lise Ramambason , et al. (10 additional authors not shown)

Abstract: Spiral arms are some of the most spectacular features in disc galaxies, and also present in our own Milky Way. It has been argued that star formation should proceed more efficiently in spiral arms as a result of gas compression. Yet, observational studies have so far yielded contradictory results. Here we examine arm/interarm surface density contrasts at ~100 pc resolution in 28 spiral galaxies fr… ▽ More Spiral arms are some of the most spectacular features in disc galaxies, and also present in our own Milky Way. It has been argued that star formation should proceed more efficiently in spiral arms as a result of gas compression. Yet, observational studies have so far yielded contradictory results. Here we examine arm/interarm surface density contrasts at ~100 pc resolution in 28 spiral galaxies from the PHANGS survey. We find that the arm/interarm contrast in stellar mass surface density (Sigma_*) is very modest, typically a few tens of percent. This is much smaller than the contrasts measured for molecular gas (Sigma_mol) or star formation rate (Sigma_SFR) surface density, which typically reach a factor of ~2-3. Yet, Sigma_mol and Sigma_SFR contrasts show a significant correlation with the enhancement in Sigma_*, suggesting that the small stellar contrast largely dictates the stronger accumulation of gas and star formation. All these contrasts increase for grand-design spirals compared to multi-armed and flocculent systems (and for galaxies with high stellar mass). The median star formation efficiency (SFE) of the molecular gas is 16% higher in spiral arms than in interarm regions, with a large scatter, and the contrast increases significantly (median SFE contrast 2.34) for regions of particularly enhanced stellar contrast (Sigma_* contrast >1.97). The molecular-to-atomic gas ratio (Sigma_mol/Sigma_atom) is higher in spiral arms, pointing to a transformation of atomic to molecular gas. In conclusion, the boost in the star formation efficiency of molecular gas in spiral arms is generally modest or absent, except for locations with exceptionally large stellar contrasts. (abridged) △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 26 pages, 16 figures. Accepted for publication in A&A

Simulating nearby disc galaxies on the main star formation sequence I. Bar formation and the building of the central gas reservoir

Authors: Pierrick Verwilghen, Eric Emsellem, Florent Renaud, Milena Valentini, Jiayi Sun, Sarah Jeffreson, Ralf S. Klessen, Mattia C. Sormani, Ashley. T. Barnes, Klaus Dolag, Kathryn Grasha, Fu-Heng Liang, Sharon Meidt, Justus Neumann, Miguel Querejeta, Eva Schinnerer, Thomas G. Williams

Abstract: Past studies have long emphasised the key role played by galactic stellar bars in the context of disc secular evolution, via the redistribution of gas and stars, the triggering of star formation, and the formation of prominent structures such as rings and central mass concentrations. However, the exact physical processes acting on those structures, as well as the timescales associated with the bui… ▽ More Past studies have long emphasised the key role played by galactic stellar bars in the context of disc secular evolution, via the redistribution of gas and stars, the triggering of star formation, and the formation of prominent structures such as rings and central mass concentrations. However, the exact physical processes acting on those structures, as well as the timescales associated with the building and consumption of central gas reservoirs are still not well understood. We are building a suite of hydro-dynamical RAMSES simulations of isolated, low-redshift galaxies that mimic the properties of the PHANGS sample. The initial conditions of the models reproduce the observed stellar mass, disc scale length, or gas fraction, and this paper presents a first subset of these models. Most of our simulated galaxies develop a prominent bar structure, which itself triggers central gas fuelling and the building of an over-density with a typical scale of 100-1000 pc. We confirm that if the host galaxy features an ellipsoidal component, the formation of the bar and gas fuelling are delayed. We show that most of our simulations follow a common time evolution, when accounting for mass scaling and the bar formation time. In our simulations, the stellar mass of $10^{10}$~M$_{\odot}$ seems to mark a change in the phases describing the time evolution of the bar and its impact on the interstellar medium. In massive discs (M$_{\star} \geq 10^{10}$~M$_{\odot}$), we observe the formation of a central gas reservoir with star formation mostly occurring within a restricted starburst region, leading to a gas depletion phase. Lower-mass systems (M$_{\star} < 10^{10}$~M$_{\odot}$) do not exhibit such a depletion phase, and show a more homogeneous spread of star-forming regions along the bar structure, and do not appear to host inner bar-driven discs or rings. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 22 pages, 17 figures

Journal ref: A&A 687, A53 (2024)

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

PHANGS-JWST: Data Processing Pipeline and First Full Public Data Release

Authors: Thomas G. Williams, Janice C. Lee, Kirsten L. Larson, Adam K. Leroy, Karin Sandstrom, Eva Schinnerer, David A. Thilker, Francesco Belfiore, Oleg V. Egorov, Erik Rosolowsky, Jessica Sutter, Joseph DePasquale, Alyssa Pagan, Travis A. Berger, Gagandeep S. Anand, Ashley T. Barnes, Frank Bigiel, Médéric Boquien, Yixian Cao, Jérémy Chastenet, Mélanie Chevance, Ryan Chown, Daniel A. Dale, Sinan Deger, Cosima Eibensteiner , et al. (33 additional authors not shown)

Abstract: The exquisite angular resolution and sensitivity of JWST is opening a new window for our understanding of the Universe. In nearby galaxies, JWST observations are revolutionizing our understanding of the first phases of star formation and the dusty interstellar medium. Nineteen local galaxies spanning a range of properties and morphologies across the star-forming main sequence have been observed as… ▽ More The exquisite angular resolution and sensitivity of JWST is opening a new window for our understanding of the Universe. In nearby galaxies, JWST observations are revolutionizing our understanding of the first phases of star formation and the dusty interstellar medium. Nineteen local galaxies spanning a range of properties and morphologies across the star-forming main sequence have been observed as part of the PHANGS-JWST Cycle 1 Treasury program at spatial scales of $\sim$5-50pc. Here, we describe pjpipe, an image processing pipeline developed for the PHANGS-JWST program that wraps around and extends the official JWST pipeline. We release this pipeline to the community as it contains a number of tools generally useful for JWST NIRCam and MIRI observations. Particularly for extended sources, pjpipe products provide significant improvements over mosaics from the MAST archive in terms of removing instrumental noise in NIRCam data, background flux matching, and calibration of relative and absolute astrometry. We show that slightly smoothing F2100W MIRI data to 0.9" (degrading the resolution by about 30 percent) reduces the noise by a factor of $\approx$3. We also present the first public release (DR1.1.0) of the pjpipe processed eight-band 2-21 $μ$m imaging for all nineteen galaxies in the PHANGS-JWST Cycle 1 Treasury program. An additional 55 galaxies will soon follow from a new PHANGS-JWST Cycle 2 Treasury program. △ Less

Submitted 9 May, 2024; v1 submitted 26 January, 2024; originally announced January 2024.

Comments: 49 pages (27 in Appendices), 54 Figures (39 in Appendices), 3 Tables. Accepted for publication in ApJS. Updated to match accepted version. Data available at https://archive.stsci.edu/hlsp/phangs/phangs-jwst

Hidden Gems on a Ring: Infant Massive Clusters and Their Formation Timeline Unveiled by ALMA, HST, and JWST in NGC 3351

Authors: Jiayi Sun, Hao He, Kyle Batschkun, Rebecca C. Levy, Kimberly Emig, M. Jimena Rodriguez, Hamid Hassani, Adam K. Leroy, Eva Schinnerer, Eve C. Ostriker, Christine D. Wilson, Alberto D. Bolatto, Elisabeth A. C. Mills, Erik Rosolowsky, Janice C. Lee, Daniel A. Dale, Kirsten L. Larson, David A. Thilker, Leonardo Ubeda, Bradley C. Whitmore, Thomas G. Williams, Ashley. T. Barnes, Frank Bigiel, Melanie Chevance, Simon C. O. Glover , et al. (16 additional authors not shown)

Abstract: We study young massive clusters (YMCs) in their embedded "infant" phase with $\sim0.\!^{\prime\prime}1$ ALMA, HST, and JWST observations targeting the central starburst ring in NGC 3351, a nearby Milky Way analog galaxy. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on… ▽ More We study young massive clusters (YMCs) in their embedded "infant" phase with $\sim0.\!^{\prime\prime}1$ ALMA, HST, and JWST observations targeting the central starburst ring in NGC 3351, a nearby Milky Way analog galaxy. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on the ALMA continuum and molecular line data, as well as ancillary measurements for the HST and JWST counterparts, we identify 14 sources as infant star clusters with high stellar and/or gas masses (${\sim}10^5\;\mathrm{M_\odot}$), small radii (${\lesssim}\,5\;\mathrm{pc}$), large escape velocities ($6{-}10\;\mathrm{km/s}$), and short free-fall times ($0.5{-}1\;\mathrm{Myr}$). Their multiwavelength properties motivate us to divide them into four categories, likely corresponding to four evolutionary stages from starless clumps to exposed HII region-cluster complexes. Leveraging age estimates for HST-identified clusters in the same region, we infer an evolutionary timeline going from $\sim$1-2 Myr before cluster formation as starless clumps, to $\sim$4-6 Myr after as exposed HII region-cluster complexes. Finally, we show that the YMCs make up a substantial fraction of recent star formation across the ring, exhibit an non-uniform azimuthal distribution without a very coherent evolutionary trend along the ring, and are capable of driving large-scale gas outflows. △ Less

Submitted 10 April, 2024; v1 submitted 25 January, 2024; originally announced January 2024.

Comments: 27 pages, 12 figures; ApJ accepted

Bottom's Dream and the amplification of filamentary gas structures and stellar spiral arms

Authors: Sharon E. Meidt, Arjen van der Wel

Abstract: Theories of spiral structure traditionally separate into tight-winding Lin-Shu spiral density waves and the swing-amplified material patterns of Goldreich & Lynden-Bell and Julian & Toomre. In this paper we consolidate these two types of spirals into a unified description, treating density waves beyond the tight-winding limit, in the regime of shearing and non-steady open spirals. This 'shearing w… ▽ More Theories of spiral structure traditionally separate into tight-winding Lin-Shu spiral density waves and the swing-amplified material patterns of Goldreich & Lynden-Bell and Julian & Toomre. In this paper we consolidate these two types of spirals into a unified description, treating density waves beyond the tight-winding limit, in the regime of shearing and non-steady open spirals. This 'shearing wave' scenario novelly captures swing amplification that enables structure formation above conventional Q thresholds. However, it also highlights the fundamental role of spiral forcing on the amplification process in general, whether the wave is shearing or not. Thus it captures resonant and non-resonant mode growth through the donkey effect described by Lynden-Bell & Kalnajs and, critically, the cessation of growth when donkey behavior is no longer permitted. Our calculations predict growth exclusive to trailing spirals above the Jeans length, the prominence of spirals across a range of orientations that increases with decreasing arm multiplicity, and a critical orientation where growth is fastest that is the same for both modes and material patterns. Predicted structures are consistent with highly regular, high-multiplicity gaseous spur features and long filaments spaced close to the Jeans scale in spirals and bars. Applied to stellar disks, conditions favor low multiplicity (m<5) open trailing spirals with pitch angles in the observed range $10 deg$<$i_p$<$50 deg$. The results of this work serve as a basis for describing spirals as a unified class of transient waves, abundantly stimulated but narrowly selected for growth depending on local conditions. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: Accepted for publication in ApJ, 30 pages, 4 figures

A Two-Component Probability Distribution Function Describes the mid-IR Emission from the Disks of Star-Forming Galaxies

Authors: Debosmita Pathak, Adam K. Leroy, Todd A. Thompson, Laura A. Lopez, Francesco Belfiore, Mederic Boquien, Daniel A. Dale, Simon C. O. Glover, Ralf S. Klessen, Eric W. Koch, Erik Rosolowsky, Karin M. Sandstrom, Eva Schinnerer, Rowan Smith, Jiayi Sun, Jessica Sutter, Thomas G. Williams, Frank Bigiel, Yixian Cao, Jeremy Chastenet, Melanie Chevance, Ryan Chown, Eric Emsellem, Christopher M. Faesi, Kirsten L. Larson , et al. (6 additional authors not shown)

Abstract: High-resolution JWST-MIRI images of nearby spiral galaxies reveal emission with complex substructures that trace dust heated both by massive young stars and the diffuse interstellar radiation field. We present high angular (0."85) and physical resolution (20-80 pc) measurements of the probability distribution function (PDF) of mid-infrared (mid-IR) emission (7.7-21 $μ$m) from 19 nearby star-formin… ▽ More High-resolution JWST-MIRI images of nearby spiral galaxies reveal emission with complex substructures that trace dust heated both by massive young stars and the diffuse interstellar radiation field. We present high angular (0."85) and physical resolution (20-80 pc) measurements of the probability distribution function (PDF) of mid-infrared (mid-IR) emission (7.7-21 $μ$m) from 19 nearby star-forming galaxies from the PHANGS-JWST Cycle-1 Treasury. The PDFs of mid-IR emission from the disks of all 19 galaxies consistently show two distinct components: an approximately log-normal distribution at lower intensities and a high-intensity power-law component. These two components only emerge once individual star-forming regions are resolved. Comparing with locations of HII regions identified from VLT/MUSE H$α$-mapping, we infer that the power-law component arises from star-forming regions and thus primarily traces dust heated by young stars. In the continuum-dominated 21 $μ$m band, the power-law is more prominent and contains roughly half of the total flux. At 7.7-11.3 $μ$m, the power-law is suppressed by the destruction of small grains (including PAHs) close to HII regions while the log-normal component tracing the dust column in diffuse regions appears more prominent. The width and shape of the log-normal diffuse emission PDFs in galactic disks remain consistent across our sample, implying a log-normal gas column density $N$(H)$\approx10^{21}$cm$^{-2}$ shaped by supersonic turbulence with typical (isothermal) turbulent Mach numbers $\approx5-15$. Finally, we describe how the PDFs of galactic disks are assembled from dusty HII regions and diffuse gas, and discuss how the measured PDF parameters correlate with global properties such as star-formation rate and gas surface density. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: 30 pages without appendix, 17 figures, (with appendix images of full sample: 56 pages, 39 figures), accepted in AJ

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)

TODDLERS: A new UV-mm emission library for star-forming regions. 1. Integration with SKIRT and public release

Authors: Anand Utsav Kapoor, Maarten Baes, Arjen van der Wel, Andrea Gebek, Peter Camps, Angelos Nersesian, Sharon E. Meidt, Aaron Smith, Sebastien Vicens, Francesco D'Eugenio, Marco Martorano, Daniela Barrientos, Nina Sanches Sartorio

Abstract: We present and publicly release a new star-forming regions emission library TODDLERS (Time evolution of Observables including Dust Diagnostics and Line Emission from Regions containing young Stars) for the publicly available radiative transfer code SKIRT. The library generation involves the spherical evolution of a homogeneous gas cloud around a young stellar cluster that accounts for stellar feed… ▽ More We present and publicly release a new star-forming regions emission library TODDLERS (Time evolution of Observables including Dust Diagnostics and Line Emission from Regions containing young Stars) for the publicly available radiative transfer code SKIRT. The library generation involves the spherical evolution of a homogeneous gas cloud around a young stellar cluster that accounts for stellar feedback processes including stellar winds, supernovae, and radiation pressure, as well as the gravitational forces on the gas. The semi-analytical evolution model is coupled with the photoionization code Cloudy to calculate time-dependent UV-mm spectral energy distributions (SEDs) from star-forming regions of varying metallicity, star-formation efficiency, birth-cloud density, and mass. The calculated SEDs include the stellar, nebular, and dust continuum emission along with a wide range of emission lines originating from H ii, photo-dissociation, and molecular gas regimes tabulated at high resolution. The SEDs incorporated in SKIRT are generated by calculating a stellar-mass normalized luminosity, which assumes that each emission source is composed of a power-law population of star-forming clouds. When compared to the previous treatment of star-forming regions in SKIRT, TODDLERS shows a better agreement with low-redshift observational data in the IR wavelength range while offering a more comprehensive line-emission support. This paves the way for a variety of applications using simulated galaxies at low and high redshift. △ Less

Submitted 20 October, 2023; v1 submitted 30 September, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS on 27/09/23

Quantifying the energy balance between the turbulent ionised gas and young stars

Authors: Oleg V. Egorov, Kathryn Kreckel, Simon C. O. Glover, Brent Groves, Francesco Belfiore, Eric Emsellem, Ralf S. Klessen, Adam K. Leroy, Sharon E. Meidt, Sumit K. Sarbadhicary, Eva Schinnerer, Elizabeth J. Watkins, Brad C. Whitmore, Ashley T. Barnes, Enrico Congiu, Daniel A. Dale, Kathryn Grasha, Kirsten L. Larson, Janice C. Lee, J. Eduardo Méndez-Delgado, David A. Thilker, Thomas G. Williams

Abstract: We investigate the ionised gas morphology, excitation properties, and kinematics in 19 nearby star-forming galaxies from the PHANGS-MUSE survey. We directly compare the kinetic energy of expanding superbubbles and the turbulent motions in the interstellar medium with the mechanical energy deposited by massive stars in the form of winds and supernovae, with the aim to answer whether the stellar fee… ▽ More We investigate the ionised gas morphology, excitation properties, and kinematics in 19 nearby star-forming galaxies from the PHANGS-MUSE survey. We directly compare the kinetic energy of expanding superbubbles and the turbulent motions in the interstellar medium with the mechanical energy deposited by massive stars in the form of winds and supernovae, with the aim to answer whether the stellar feedback is responsible for the observed turbulent motions and to quantify the fraction of mechanical energy retained in the superbubbles. Based on the distribution of the flux and velocity dispersion in the H$α$ line, we select 1484 regions of locally elevated velocity dispersion ($σ$(H$α$)>45 km/s), including at least 171 expanding superbubbles. We analyse these regions and relate their properties to those of the young stellar associations and star clusters identified in PHANGS-HST data. We find a good correlation between the kinetic energy of the ionised gas and the total mechanical energy input from supernovae and stellar winds from the stellar associations, with a typical coupling efficiency of 10-20%. The contribution of mechanical energy by the supernovae alone is not sufficient to explain the measured kinetic energy of the ionised gas, which implies that pre-supernova feedback in the form of radiation/thermal pressure and winds is necessary. We find that the gas kinetic energy decreases with metallicity for our sample covering Z=0.5-1.0 Zsun, reflecting the lower impact of stellar feedback. For the sample of superbubbles, we find that about 40% of the young stellar associations are preferentially located in their rims. We also find a slightly higher (by ~15%) fraction of the youngest (<3 Myr) stellar associations in the rims of the superbubbles than in the centres, and the opposite for older associations, which implies possible propagation or triggering of star formation. △ Less

Submitted 17 August, 2023; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: Accepted by A&A. The abstract is abridged. 31 pages (including 8 pages in appendix), 20 figures, 2 tables

Journal ref: A&A 678, A153 (2023)

On the Tremaine-Weinberg method: how much can we trust gas tracers to measure pattern speeds?

Authors: Olga Borodina, Thomas G. Williams, Mattia C. Sormani, Sharon Meidt, Eva Schinnerer

Abstract: Pattern speeds are a fundamental parameter of the dynamical features (e.g. bars, spiral arms) of a galaxy, setting resonance locations. Pattern speeds are not directly observable, so the Tremaine-Weinberg (TW) method has become the most common method used to measure them in galaxies. However, it has not been tested properly whether this method can straightforwardly be applied to gas tracers, despi… ▽ More Pattern speeds are a fundamental parameter of the dynamical features (e.g. bars, spiral arms) of a galaxy, setting resonance locations. Pattern speeds are not directly observable, so the Tremaine-Weinberg (TW) method has become the most common method used to measure them in galaxies. However, it has not been tested properly whether this method can straightforwardly be applied to gas tracers, despite this being widely done in the literature. When applied to observations, the TW method may return invalid results, which are difficult to diagnose due to a lack of ground truth for comparison. Although some works applying the TW method to simulated galaxies exist, only stellar populations have been tested. Therefore, here we explore the applicability of the TW method for gas gracers, by applying it to hydrodynamical simulations of galaxies, where we know the true value of the bar pattern speed. We perform some simple tests to see if the TW method has a physically reasonable output. We add different kinds of uncertainties (e.g. in position angle or flux) to the data to mock observational errors based on the magnitude of uncertainty present in the observations. Second, we test the method on 3D simulations with chemical networks. We show that in general, applying TW to observations of gas will not recover the true pattern speed. These results have implications for many "pattern speeds" reported in the literature, and based on these tests we also give some best practices for measuring pattern speeds using gas tracers going forwards. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: 9 pages, 8 figures, submitted to 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

The impact of HII regions on Giant Molecular Cloud properties in nearby galaxies sampled by PHANGS ALMA and MUSE

Authors: Antoine Zakardjian, Jérôme Pety, Cinthya N. Herrera, Annie Hughes, Elias Oakes, Kathryn Kreckel, Chris Faesi, Simon C. O. Glover, Brent Groves, Ralf S. Klessen, Sharon Meidt, Ashley Barnes, Francesco Belfiore, Ivana Bešlić, Frank Bigiel, Guillermo A. Blanc, Mélanie Chevance, Daniel A. Dale, Jakob den Brok, Cosima Eibensteiner, Eric Emsellem, Axel García-Rodríguez, Kathryn Grasha, Eric W. Koch, Adam K. Leroy , et al. (14 additional authors not shown)

Abstract: We identify giant molecular clouds (GMCs) associated with HII regions for a sample of 19 nearby galaxies using catalogs of GMCs and H regions released by the PHANGS-ALMA and PHANGS-MUSE surveys, using the overlap of the CO and Hα emission as the key criterion for physical association. We compare the distributions of GMC and HII region properties for paired and non-paired objects. We investigate co… ▽ More We identify giant molecular clouds (GMCs) associated with HII regions for a sample of 19 nearby galaxies using catalogs of GMCs and H regions released by the PHANGS-ALMA and PHANGS-MUSE surveys, using the overlap of the CO and Hα emission as the key criterion for physical association. We compare the distributions of GMC and HII region properties for paired and non-paired objects. We investigate correlations between GMC and HII region properties among galaxies and across different galactic environments to determine whether GMCs that are associated with HII regions have significantly distinct physical properties to the parent GMC population. We identify trends between the Hα luminosity of an HII region and the CO peak brightness and the molecular mass of GMCs that we tentatively attribute to a direct physical connection between the matched objects, and which arise independently of underlying environmental variations of GMC and HII region properties within galaxies. The study of the full sample nevertheless hides a large variability galaxy by galaxy. Our results suggests that at the ~100 pc scales accessed by the PHANGS-ALMA and PHANGS-MUSE data, pre-supernova feedback mechanisms in HII regions have a subtle but measurable impact on the properties of the surrounding molecular gas, as inferred from CO observations. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Journal ref: A&A 678, A171 (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)

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

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

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

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: 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: 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

Molecular Clouds as Gravitational Instabilities in Rotating Disks: A Modified Stability Criterion

Authors: Sharon E. Meidt

Abstract: Molecular gas disks are generally Toomre stable ($Q_T>$1) and yet clearly gravitationally unstable to structure formation as evidenced by the existence of molecular clouds and ongoing star formation. This paper adopts a 3D perspective to obtain a general picture of instabilities in flattened rotating disks, using the 3D dispersion relation to describe how disks evolve when perturbed over their ver… ▽ More Molecular gas disks are generally Toomre stable ($Q_T>$1) and yet clearly gravitationally unstable to structure formation as evidenced by the existence of molecular clouds and ongoing star formation. This paper adopts a 3D perspective to obtain a general picture of instabilities in flattened rotating disks, using the 3D dispersion relation to describe how disks evolve when perturbed over their vertical extents. By explicitly adding a vertical perturbation to an unperturbed equilibrium disk, stability is shown to vary with height above the mid-plane. Near to $z$=0 where the equilibrium density is roughly constant, instability takes on a Jeans-like quality, occurring on scales larger than the Jeans length and subject to a threshold $Q_M=κ^2/(4πGρ)=1$ or roughly $Q_T\approx 2$. Far from the mid-plane, on the other hand, stability is pervasive, and the threshold for the total disk (out to $z=\pm\infty$) to be stabilized is lowered to $Q_T=1$ as a consequence. In this new framework, gas disks are able to fragment through partial 3D instability even where total 2D instability is suppressed. The growth rates of the fragments formed via 3D instability are comparable to, or faster than, Toomre instabilities. The rich structure in molecular disks on the scale of 10s of pc can thus be viewed as a natural consequence of their 3D nature and their exposure to a variety of vertical perturbations acting on roughly a disk scale height, i.e. due to their situation within the more extended galaxy potential, participation in the disk-halo flow, and exposure to star formation feedback. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: Accepted for publication in ApJ, 23 pages, 3 figures

Investigating the link between inner gravitational potential and star-formation quenching in CALIFA galaxies

Authors: V. Kalinova, D. Colombo, S. F. Sánchez, E. Rosolowsky, K. Kodaira, R. García-Benito, S. E. Meidt, T. A. Davis, A. B. Romeo, S. -Y. Yu, R. González Delgado, E. A. D. Lacerda

Abstract: It has been suggested that the gravitational potential can have a significant role in suppressing the star formation in the nearby galaxies. To establish observational constrains on this scenario, we investigate the connection between the dynamics, through the circular velocity curves (CVCs) as a proxy of the inner gravitational potential, and star formation quenching in 215 non-active galaxies ac… ▽ More It has been suggested that the gravitational potential can have a significant role in suppressing the star formation in the nearby galaxies. To establish observational constrains on this scenario, we investigate the connection between the dynamics, through the circular velocity curves (CVCs) as a proxy of the inner gravitational potential, and star formation quenching in 215 non-active galaxies across Hubble sequence from the Calar Alto Legacy Integral Field Area (CALIFA) survey. Our results show that galaxies with similar CVCs tend to have a certain star-formation quenching pattern. To explore these findings in more details, we construct kpc-resolved relations of the equivalent width of the H$α$ ($W_{{\rm H}α}$) versus the amplitude ($V_c$) and shape ($β= d\ln V_c/ d\ln R$) of the circular velocity at given radius. We find that the $W_{{\rm H}α}-V_c$ is a declining relationship, where the retired regions of the galaxies (the ones with $W_{{\rm H}α}$ values below 3 Å) tend to have higher $V_c$. Differently, $W_{{\rm H}α}-β$ is a bi-modal relationship, characterised by two peaks: concentration of the star forming regions at a positive $β$ (rising CVC) and another one of the retired regions with a negative $β$ (declining CVC). Our results show that both the amplitude of the CVC, driven by the mass of the galaxies, and the shape of the CVC, reflecting the internal structure of the galaxies, play an important role in galaxy's quenching history. △ Less

Submitted 8 July, 2022; originally announced July 2022.

Comments: 16 pages, 12 figures, Accepted for publication in A&A

Journal ref: A&A 665, A90 (2022)

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

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

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

Submitted 14 June, 2022; originally announced June 2022.

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

Journal ref: AJ, 164, 43 (2022)

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

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

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

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

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

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

LEGA-C: analysis of dynamical masses from ionized gas and stellar kinematics at z~0.8

Authors: Caroline M. S. Straatman, Arjen van der Wel, Josha van Houdt, Rachel Bezanson, Eric F. Bell, Pieter van Dokkum, Francesco D'Eugenio, Marijn Franx, Anna Gallazzi, Anna de Graaff, Michael Maseda, Sharon E. Meidt, Adam Muzzin, David Sobral, Po-Feng Wu

Abstract: We compare dynamical mass estimates based on spatially extended stellar and ionized gas kinematics ($\mathrm{M_{dyn,*}}$ and $\mathrm{M_{dyn,eml}}$, respectively) of 157 star forming galaxies at $0.6\leq z<1$. Compared to $z\sim0$, these galaxies have enhanced star formation rates, with stellar feedback likely affecting the dynamics of the gas. We use LEGA-C DR3, the highest redshift dataset provi… ▽ More We compare dynamical mass estimates based on spatially extended stellar and ionized gas kinematics ($\mathrm{M_{dyn,*}}$ and $\mathrm{M_{dyn,eml}}$, respectively) of 157 star forming galaxies at $0.6\leq z<1$. Compared to $z\sim0$, these galaxies have enhanced star formation rates, with stellar feedback likely affecting the dynamics of the gas. We use LEGA-C DR3, the highest redshift dataset providing sufficiently deep measurements of a $K_s-$band limited sample. For $\mathrm{M_{dyn,*}}$ we use Jeans Anisotropic Multi-Gaussian Expansion models. For $\mathrm{M_{dyn,eml}}$ we first fit a custom model of a rotating exponential disk with uniform dispersion, whose light is projected through a slit and corrected for beam smearing. We then apply an asymmetric drift correction based on assumptions common in the literature to the fitted kinematic components to obtain the circular velocity, assuming hydrostatic equilibrium. Within the half-light radius, $\mathrm{M_{dyn,eml}}$ is on average lower than $\mathrm{M_{dyn,*}}$, with a mean offset of $-0.15\pm0.016$ dex and galaxy-to-galaxy scatter of $0.19$ dex, reflecting the combined random uncertainty. While data of higher spatial resolution are needed to understand this small offset, it supports the assumption that the galaxy-wide ionized gas kinematics do not predominantly originate from disruptive events such as star formation driven outflows. However, a similar agreement can be obtained without modeling from the integrated emission line dispersions for axis ratios $q<0.8$. This suggests that our current understanding of gas kinematics is not sufficient to efficiently apply asymmetric drift corrections to improve dynamical mass estimates compared to observations lacking the $S/N$ required for spatially extended dynamics. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: 19 pages, 7 figures, accepted for publication by ApJ, in press

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

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

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

Submitted 13 January, 2022; originally announced January 2022.

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

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

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

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

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

Submitted 4 January, 2022; originally announced January 2022.

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

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

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

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

Submitted 17 November, 2021; originally announced November 2021.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Abstract: We measure the low-J CO line ratio R21=CO(2-1)/CO(1-0), R32=CO(3-2)/CO(2-1), and R31 = CO(3-2)/CO(1-0) using whole-disk CO maps of nearby galaxies. We draw CO(2-1) from PHANGS--ALMA, HERACLES, and follow-up IRAM surveys; CO(1-0) from COMING and the Nobeyama CO Atlas of Nearby Spiral Galaxies; and CO(3-2) from the JCMT NGLS and APEX LASMA mapping. Altogether this yields 76, 47, and 29 maps of R21,… ▽ More We measure the low-J CO line ratio R21=CO(2-1)/CO(1-0), R32=CO(3-2)/CO(2-1), and R31 = CO(3-2)/CO(1-0) using whole-disk CO maps of nearby galaxies. We draw CO(2-1) from PHANGS--ALMA, HERACLES, and follow-up IRAM surveys; CO(1-0) from COMING and the Nobeyama CO Atlas of Nearby Spiral Galaxies; and CO(3-2) from the JCMT NGLS and APEX LASMA mapping. Altogether this yields 76, 47, and 29 maps of R21, R32, and R31 at 20" \sim 1.3 kpc resolution, covering 43, 34, and 20 galaxies. Disk galaxies with high stellar mass, log10 M_* [Msun]=10.25-11 and star formation rate, SFR=1-5 Msun/yr, dominate the sample. We find galaxy-integrated mean values and 16%-84% range of R21 = 0.65 (0.50-0.83), R32=0.50 (0.23-0.59), and R31=0.31 (0.20-0.42). We identify weak trends relating galaxy-integrated line ratios to properties expected to correlate with excitation, including SFR/M_* and SFR/L_CO. Within galaxies, we measure central enhancements with respect to the galaxy-averaged value of \sim 0.18^{+0.09}_{-0.14} dex for R21, 0.27^{+0.13}_{-0.15} dex for R31, and 0.08^{+0.11}_{-0.09} dex for R32. All three line ratios anti-correlate with galactocentric radius and positively correlate with the local star formation rate surface density and specific star formation rate, and we provide approximate fits to these relations. The observed ratios can be reasonably reproduced by models with low temperature, moderate opacity, and moderate densities, in good agreement with expectations for the cold ISM. Because the line ratios are expected to anti-correlate with the CO(1-0)-to-H_2 conversion factor, alphaCO^(1-0), these results have general implications for the interpretation of CO emission from galaxies. △ Less

Submitted 23 September, 2021; originally announced September 2021.

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

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

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

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

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

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

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

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

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

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

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

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

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