Showing 1–50 of 76 results for author: Koch, E

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)

PHANGS-MeerKAT and MHONGOOSE HI observations of nearby spiral galaxies: physical drivers of the molecular gas fraction, $R_{\mathrm{mol}}$

Authors: Cosima Eibensteiner, Jiayi Sun, Frank Bigiel, Adam K. Leroy, Eva Schinnerer, Erik Rosolowsky, Sushma Kurapati, D. J. Pisano, W. J. G de Blok, Ashley T. Barnes, Mallory Thorp, Dario Colombo, Eric W. Koch, I-Da Chiang, Eve C. Ostriker, Eric J. Murphy, Nikki Zabel, Sebstian Laudage, Filippo M. Maccagni, Julia Healy, Srikrishna Sekhar, Dyas Utomo, Jakob den Brok, Yixian Cao, Mélanie Chevance , et al. (14 additional authors not shown)

Abstract: The molecular-to-atomic gas ratio is crucial to the evolution of the interstellar medium in galaxies. We investigate the balance between the atomic ($Σ_{\rm HI}$) and molecular gas ($Σ_{\rm H2}$) surface densities in eight nearby star-forming galaxies using new high-quality observations from MeerKAT and ALMA (for HI and CO, respectively). We define the molecular gas ratio as… ▽ More The molecular-to-atomic gas ratio is crucial to the evolution of the interstellar medium in galaxies. We investigate the balance between the atomic ($Σ_{\rm HI}$) and molecular gas ($Σ_{\rm H2}$) surface densities in eight nearby star-forming galaxies using new high-quality observations from MeerKAT and ALMA (for HI and CO, respectively). We define the molecular gas ratio as $R_{\rm mol} = Σ_{\rm H2} / Σ_{\rm HI}$ and measure how it depends on local conditions in the galaxy disks using multi-wavelength observations. We find that, depending on the galaxy, HI is detected at $>3σ$ out to 20-120 kpc in galactocentric radius ($r_{\rm gal}$). The typical radius at which $Σ_{\rm HI}$ reaches 1~$\rm M_\odot~pc^{-2}$ is $r_{\rm HI}\approx22$~kpc, which corresponds to 1-3 times the optical radius ($r_{25}$). $R_{\rm mol}$ correlates best with the dynamical equilibrium pressure, P$_{\rm DE}$, among potential drivers studied, with a median correlation coefficient of $<ρ>=0.89$. Correlations between $R_{\rm mol}$ and star formation rate, total gas and stellar surface density, metallicity, and $Σ_{\rm SFR}$/P$_{\rm DE}$ are present but somewhat weaker. Our results also show a direct correlation between P$_{\rm DE}$ and $Σ_{\rm SFR}$, supporting self-regulation models. Quantitatively, we measure similar scalings as previous works and attribute the modest differences that we find to the effect of varying resolution and sensitivity. At $r_{\rm gal} {\gtrsim}0.4~r_{25}$, atomic gas dominates over molecular gas, and at the balance of these two gas phases, we find that the baryon mass is dominated by stars, with $Σ_{*} > 5~Σ_{\rm gas}$. Our study constitutes an important step in the statistical investigation of how local galaxy properties impact the conversion from atomic to molecular gas in nearby galaxies. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: accepted for publication in A&A; 20 pages, 12 Figures (+4 appendix pages)

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER). VI. The High-Mass Stellar Initial Mass Function of M33

Authors: Tobin M. Wainer, Benjamin F. Williams, L. Clifton Johnson, Daniel R. Weisz, Julianne J. Dalcanton, Anil C. Seth, Andrew Dolphin, Meredith J. Durbin, Eric F. Bell, Zhuo Chen, Puragra Guhathakurta, Eric W. Koch, Christina W. Lindberg, Erik Rosolowsky, Karin M. Sandstrom, Evan D. Skillman, Adam Smercina, Estephani E. TorresVillanueva

Abstract: We measure the high-mass stellar initial mass function (IMF) from resolved stars in M33 young stellar clusters. Leveraging \textit{Hubble Space Telescope's} high resolving power, we fully model the IMF probabilistically. We first model the optical CMD of each cluster to constrain its power-law slope $Γ$, marginalized over other cluster parameters in the fit (e.g., cluster age, mass, and radius). W… ▽ More We measure the high-mass stellar initial mass function (IMF) from resolved stars in M33 young stellar clusters. Leveraging \textit{Hubble Space Telescope's} high resolving power, we fully model the IMF probabilistically. We first model the optical CMD of each cluster to constrain its power-law slope $Γ$, marginalized over other cluster parameters in the fit (e.g., cluster age, mass, and radius). We then probabilistically model the distribution of MF slopes for a highly strict cluster sample of 9 clusters more massive than log(Mass/M$_{\odot}$)=3.6; above this mass, all clusters have well-populated main sequences of massive stars and should have accurate recovery of their MF slopes, based on extensive tests with artificial clusters. We find the ensemble IMF is best described by a mean high-mass slope of $\overlineΓ = 1.49\pm0.18$, with an intrinsic scatter of $σ^{2}_Γ = 0.02^{+0.16}_{0.00}$, consistent with a universal IMF. We find no dependence of the IMF on environmental impacts such as the local star formation rate or galactocentric radius within M33, which serves as a proxy for metallicity. This $\overlineΓ$ measurement is consistent with similar measurements in M31, despite M33 having a much higher star formation rate intensity. While this measurement is formally consistent with the canonical Kroupa ($Γ= 1.30$) IMF, as well as the Salpeter ($Γ= 1.35)$) value, it is the second Local Group cluster sample to show evidence for a somewhat steeper high-mass IMF slope. We explore the impacts a steeper IMF slope has on a number of astronomical sub-fields. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: Accepted for publication in ApJ. 9 Figures, 1 Table

The Fraction of Dust Mass in the Form of PAHs on 10-50 pc Scales in Nearby Galaxies

Authors: Jessica Sutter, Karin Sandstrom, Jérémy Chastenet, Adam K. Leroy, Eric W. Koch, Thomas G. Williams, Ryan Chown, Francesco Belfiore, Frank Bigiel, Médéric Boquien, Yixian Cao, Mélanie Chevance, Daniel A. Dale, Oleg V. Egorov, Simon C. O. Glover, Brent Groves, Ralf S. Klessen, Kathryn Kreckel, Kirsten L. Larson, Elias K. Oakes, Debosmita Pathak, Lise Ramambason, Erik Rosolowsky, Elizabeth J. Watkins

Abstract: Polycyclic aromatic hydrocarbons (PAHs) are a ubiquitous component of the interstellar medium (ISM) in z~0 massive, star-forming galaxies and play key roles in ISM energy balance, chemistry, and shielding. Wide field of view, high resolution mid-infrared (MIR) images from JWST provides the ability to map the fraction of dust in the form of PAHs and the properties of these key dust grains at 10-50… ▽ More Polycyclic aromatic hydrocarbons (PAHs) are a ubiquitous component of the interstellar medium (ISM) in z~0 massive, star-forming galaxies and play key roles in ISM energy balance, chemistry, and shielding. Wide field of view, high resolution mid-infrared (MIR) images from JWST provides the ability to map the fraction of dust in the form of PAHs and the properties of these key dust grains at 10-50 pc resolution in galaxies outside the Local Group. We use MIR JWST photometric observations of a sample of 19 nearby galaxies from the "Physics at High Angular Resolution in Nearby GalaxieS" (PHANGS) survey to investigate the variations of the PAH fraction. By comparison to lower resolution far-IR mapping, we show that a combination of the MIRI filters (R$_{\rm{PAH}}$ = [F770W+F1130W]/F2100W) traces the fraction of dust by mass in the form of PAHs (i.e., the PAH fraction, or q$_{\rm{PAH}}$). Mapping R$_{\rm{PAH}}$ across the 19 PHANGS galaxies, we find that the PAH fraction steeply decreases in HII regions, revealing the destruction of these small grains in regions of ionized gas. Outside HII regions, we find R$_{\rm{PAH}}$ is constant across the PHANGS sample with an average value of 3.43$\pm$0.98, which, for an illuminating radiation field of intensity 2-5 times that of the radiation field in the solar neighborhood, corresponds to q$_{\rm{PAH}}$ values of 3-6%. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: Accepted at ApJ, 39 pages, 25 figures

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

PHANGS-ML: dissecting multiphase gas and dust in nearby galaxies using machine learning

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

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

Submitted 6 February, 2024; originally announced February 2024.

Comments: Main results in figures 6 and 12. Submitted to ApJ, and comments are welcome!

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

JWST Reveals Star Formation Across a Spiral Arm in M33

Authors: Joshua Peltonen, Erik Rosolowsky, Thomas G. Williams, Eric W. Koch, Andrew Dolphin, Jeremy Chastenet, Julianne J. Dalcanton, Adam Ginsburg, L. Clifton Johnson, Adam K. Leroy, Theo Richardson, Karin M. Sandstrom, Sumit K. Sarbadhicary, Adam Smercina, Tobin Wainer, Benjamin F. Williams

Abstract: Young stellar objects (YSOs) are the gold standard for tracing star formation in galaxies but have been unobservable beyond the Milky Way and Magellanic Clouds. But that all changed when the James Webb Space Telescope was launched, which we use to identify YSOs in the Local Group galaxy M33, marking the first time that individual YSOs have been identified at these large distances. We present MIRI… ▽ More Young stellar objects (YSOs) are the gold standard for tracing star formation in galaxies but have been unobservable beyond the Milky Way and Magellanic Clouds. But that all changed when the James Webb Space Telescope was launched, which we use to identify YSOs in the Local Group galaxy M33, marking the first time that individual YSOs have been identified at these large distances. We present MIRI imaging mosaics at 5.6 and 21 microns that cover a significant portion of one of M33's spiral arms that has existing panchromatic imaging from the Hubble Space Telescope and deep ALMA CO measurements. Using these MIRI and Hubble Space Telescope images, we identify point sources using the new DOLPHOT MIRI module. We identify 793 candidate YSOs from cuts based on colour, proximity to giant molecular clouds (GMCs), and visual inspection. Similar to Milky Way GMCs, we find that higher mass GMCs contain more YSOs and YSO emission, which further shows YSOs identify star formation better than most tracers that cannot capture this relationship at cloud scales. We find evidence of enhanced star formation efficiency in the southern spiral arm by comparing the YSOs to the molecular gas mass. △ Less

Submitted 18 December, 2023; v1 submitted 14 December, 2023; originally announced December 2023.

Comments: 12 pages, 10 figures, 1 tables, accepted for publication at MNRAS

The PHANGS-AstroSat Atlas of Nearby Star Forming Galaxies

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

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

Submitted 10 December, 2023; originally announced December 2023.

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

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

Resolved Measurements of the CO-to-H$_2$ Conversion Factor in 37 Nearby Galaxies

Authors: I-Da Chiang, Karin M. Sandstrom, Jeremy Chastenet, Alberto D. Bolatto, Eric W. Koch, Adam K. Leroy, Jiayi Sun, Yu-Hsuan Teng, Thomas G. Williams

Abstract: We measure the CO-to-H$_2$ conversion factor ($α_\mathrm{CO}$) in 37 galaxies at 2~kpc resolution, using dust surface density inferred from far-infrared emission as a tracer of the gas surface density and assuming a constant dust-to-metals ratio. In total, we have $\sim790$ and $\sim610$ independent measurements of $α_\mathrm{CO}$ for CO (2-1) and (1-0), respectively. The mean values for… ▽ More We measure the CO-to-H$_2$ conversion factor ($α_\mathrm{CO}$) in 37 galaxies at 2~kpc resolution, using dust surface density inferred from far-infrared emission as a tracer of the gas surface density and assuming a constant dust-to-metals ratio. In total, we have $\sim790$ and $\sim610$ independent measurements of $α_\mathrm{CO}$ for CO (2-1) and (1-0), respectively. The mean values for $α_\mathrm{CO~(2-1)}$ and $α_\mathrm{CO~(1-0)}$ are $9.3^{+4.6}_{-5.4}$ and $4.2^{+1.9}_{-2.0}~M_\odot~pc^{-2}~(K~km~s^{-1})^{-1}$, respectively. The CO-intensity-weighted mean for $α_\mathrm{CO~(2-1)}$ is 5.69, and 3.33 for $α_\mathrm{CO~(1-0)}$. We examine how $α_\mathrm{CO}$ scales with several physical quantities, e.g.\ star-formation rate (SFR), stellar mass, and dust-mass-weighted average interstellar radiation field strength ($\overline{U}$). Among them, $\overline{U}$, $Σ_{\rm SFR}$, and integrated CO intensity ($W_\mathrm{CO}$) have the strongest anti-correlation with spatially resolved $α_\mathrm{CO}$. We provide linear regression results to \aco for all quantities tested. At galaxy integrated scales, we observe significant correlations between $α_\mathrm{CO}$ and $W_\mathrm{CO}$, metallicity, $\overline{U}$, and $Σ_{\rm SFR}$. We also find that the normalized $α_\mathrm{CO}$ decreases with stellar mass surface density ($Σ_\star$) in the high surface density regions ($Σ_\star\geq100~{\rm M_\odot~pc^{-2}}$), following the power-law relations $α_\mathrm{CO~(2-1)}\proptoΣ_\star^{-0.5}$ and $α_\mathrm{CO~(1-0)}\proptoΣ_\star^{-0.2}$. The power-law index is insensitive to the assumed dust-to-metals ratio. (abridged) △ Less

Submitted 26 December, 2023; v1 submitted 1 November, 2023; originally announced November 2023.

Comments: 27 pages, 13 figures, 4 tables, submitted to ApJ

Resolved low-J $^{12}$CO excitation at 190 parsec resolution across NGC 2903 and NGC 3627

Authors: J. S. den Brok, A. K. Leroy, A. Usero, E. Schinnerer, E. Rosolowsky, E. W. Koch, M. Querejeta, D. Liu, F. Bigiel, A. T. Barnes, M. Chevance, D. Colombo, D. A. Dale, S. C. O. Glover, M. J. Jimenez-Donaire, Y. -H. Teng, T. G. Williams

Abstract: The low-$J$ rotational transitions of $^{12}$CO are commonly used to trace the distribution of molecular gas in galaxies. Their ratios are sensitive to excitation and physical conditions in the molecular gas. Spatially resolved studies of CO ratios are still sparse and affected by flux calibration uncertainties, especially since most do not have high angular resolution or do not have short-spacing… ▽ More The low-$J$ rotational transitions of $^{12}$CO are commonly used to trace the distribution of molecular gas in galaxies. Their ratios are sensitive to excitation and physical conditions in the molecular gas. Spatially resolved studies of CO ratios are still sparse and affected by flux calibration uncertainties, especially since most do not have high angular resolution or do not have short-spacing information and hence miss any diffuse emission. We compare the low-$J$ CO ratios across the disk of two massive, star-forming spiral galaxies NGC2903 and NGC3627 to investigate whether and how local environments drive excitation variations at GMC scales. We use Atacama Large Millimeter Array (ALMA) observations of the three lowest-$J$ CO transitions at a common angular resolution of 4$''$ (190pc). We measure median line ratios of $R_{21}=0.67^{+0.13}_{-0.11}$, $R_{32}=0.33^{+0.09}_{-0.08}$, and $R_{31}=0.24^{+0.10}_{-0.09}$ across the full disk of NGC3627. We see clear CO line ratio variation across the galaxy consistent with changes in temperature and density of the molecular gas. In particular, toward the center, $R_{21}$, $R_{32}$, and $R_{31}$ increase by 35\%, 50\%, and 66\%, respectively compared to their average disk values. The overall line ratio trends suggest that CO(3-2) is more sensitive to changes in the excitation conditions than the two lower-$J$ transitions. Furthermore, we find a similar radial $R_{32}$ trend in NGC2903, albite a larger disk-wide average of $\langle R_{32}\rangle=0.47^{+0.14}_{-0.08}$. We conclude that the CO low-$J$ line ratios vary across environments in such a way that they can trace changes in the molecular gas conditions, with the main driver being changes in temperature. △ Less

Submitted 27 October, 2023; originally announced October 2023.

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

Where do stars explode in the ISM? -- The distribution of dense gas around massive stars and supernova remnants in M33

Authors: Sumit K. Sarbadhicary, Jordan Wagner, Eric W. Koch, Ness Mayker Chen, Adam K. Leroy, Natalia Lahén, Erik Rosolowsky, Kathryn F. Neugent, Chang-Goo Kim, Laura Chomiuk, Julianne J. Dalcanton, Laura A. Lopez, Nickolas M. Pingel, Remy Indebetouw, Thomas G. Williams, Elizabeth Tarantino, Jennifer Donovan Meyer, Evan D. Skillman, Adam Smercina, Amanda A. Kepley, Eric J. Murphy, Jay Strader, Tony Wong, Snežana Stanimirović, Vicente Villanueva , et al. (5 additional authors not shown)

Abstract: Star formation in galaxies is regulated by turbulence, outflows, gas heating and cloud dispersal -- processes which depend sensitively on the properties of the interstellar medium (ISM) into which supernovae (SNe) explode. Unfortunately, direct measurements of ISM environments around SNe remain scarce, as SNe are rare and often distant. Here we demonstrate a new approach: mapping the ISM around th… ▽ More Star formation in galaxies is regulated by turbulence, outflows, gas heating and cloud dispersal -- processes which depend sensitively on the properties of the interstellar medium (ISM) into which supernovae (SNe) explode. Unfortunately, direct measurements of ISM environments around SNe remain scarce, as SNe are rare and often distant. Here we demonstrate a new approach: mapping the ISM around the massive stars that are soon to explode. This provides a much larger census of explosion sites than possible with only SNe, and allows comparison with sensitive, high-resolution maps of the atomic and molecular gas from the Jansky VLA and ALMA. In the well-resolved Local Group spiral M33, we specifically observe the environments of red supergiants (RSGs, progenitors of Type II SNe), Wolf-Rayet stars (WRs, tracing stars $>$30 M$_{\odot}$, and possibly future stripped-envelope SNe), and supernova remnants (SNRs, locations where SNe have exploded). We find that massive stars evolve not only in dense, molecular-dominated gas (with younger stars in denser gas), but also a substantial fraction ($\sim$45\% of WRs; higher for RSGs) evolve in lower-density, atomic-gas-dominated, inter-cloud media. We show that these measurements are consistent with expectations from different stellar-age tracer maps, and can be useful for validating SN feedback models in numerical simulations of galaxies. Along with the discovery of a 20-pc diameter molecular gas cavity around a WR, these findings re-emphasize the importance of pre-SN/correlated-SN feedback evacuating the dense gas around massive stars before explosion, and the need for high-resolution (down to pc-scale) surveys of the multi-phase ISM in nearby galaxies. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: 34 pages, 14 figures. Submitted to ApJ. Comments welcome! The density distributions will be made publicly available after journal acceptance of manuscript. Please feel free to contact us in the meantime if you would like to use them

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER). V. The Structure of M33 in Resolved Stellar Populations

Authors: Adam Smercina, Julianne J. Dalcanton, Benjamin F. Williams, Meredith J. Durbin, Margaret Lazzarini, Eric F. Bell, Yumi Choi, Andrew Dolphin, Karoline Gilbert, Puragra Guhathakurta, Eric W. Koch, Hans-Walter Rix, Erik Rosolowsky, Anil Seth, Evan D. Skillman, Daniel R. Weisz

Abstract: We present a detailed analysis of the the structure of the Local Group flocculent spiral galaxy M33, as measured using the Panchromatic Hubble Andromeda Treasury Triangulum Extended Region (PHATTER) survey. Leveraging the multiwavelength coverage of PHATTER, we find that the oldest populations are dominated by a smooth exponential disk with two distinct spiral arms and a classical central bar $-$… ▽ More We present a detailed analysis of the the structure of the Local Group flocculent spiral galaxy M33, as measured using the Panchromatic Hubble Andromeda Treasury Triangulum Extended Region (PHATTER) survey. Leveraging the multiwavelength coverage of PHATTER, we find that the oldest populations are dominated by a smooth exponential disk with two distinct spiral arms and a classical central bar $-$ completely distinct from what is seen in broadband optical imaging, and the first-ever confirmation of a bar in M33. We estimate a bar extent of $\sim$1 kpc. The two spiral arms are asymmetric in orientation and strength, and likely represent the innermost impact of the recent tidal interaction responsible for M33's warp at larger scales. The flocculent multi-armed morphology for which M33 is known is only visible in the young upper main sequence population, which closely tracks the morphology of the ISM. We investigate the stability of M33's disk, finding $Q{\sim}1$ over the majority of the disk. We fit multiple components to the old stellar density distribution and find that, when considering recent stellar kinematics, M33's bulk structure favors the inclusion of an accreted halo component, modeled as a broken power-law. The best-fit halo model has an outer power-law index of $-$3 and accurately describes observational evidence of M33's stellar halo from both resolved stellar spectroscopy in the disk and its stellar populations at large radius. Integrating this profile yields a total halo stellar mass of ${\sim}5{\times}10^8\ M_{\odot}$, giving a total stellar halo mass fraction of 16%, most of which resides in the innermost 2.5 kpc. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 31 pages, 17 figures, 5 tables, Accepted for publication in ApJ

RomAndromeda: The Roman Survey of the Andromeda Halo

Authors: Arjun Dey, Joan Najita, Carrie Filion, Jiwon Jesse Han, Sarah Pearson, Rosemary Wyse, Adrien C. R. Thob, Borja Anguiano, Miranda Apfel, Magda Arnaboldi, Eric F. Bell, Leandro Beraldo e Silva, Gurtina Besla, Aparajito Bhattacharya, Souradeep Bhattacharya, Vedant Chandra, Yumi Choi, Michelle L. M. Collins, Emily C. Cunningham, Julianne J. Dalcanton, Ivanna Escala, Hayden R. Foote, Annette M. N. Ferguson, Benjamin J. Gibson, Oleg Y. Gnedin , et al. (28 additional authors not shown)

Abstract: As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the ha… ▽ More As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the halo (10$σ$ detection in F146, F062 of 26.5, 26.1AB mag respectively) and yield proper motions to $\sim$25 microarcsec/year (i.e., $\sim$90 km/s) for all stars brighter than F146 $\approx 23.6$ AB mag (i.e., reaching the red clump stars in the Andromeda halo). This survey will yield (through averaging) high-fidelity proper motions for all satellites and compact substructures in the Andromeda halo and will enable statistical searches for clusters in chemo-dynamical space. Adding a third epoch during the extended mission will improve these proper motions by $\sim t^{-1.5}$, to $\approx 11$ km/s, but this requires obtaining the first epoch in Year 1 of Roman operations. In combination with ongoing and imminent spectroscopic campaigns with ground-based telescopes, this Roman survey has the potential to yield full 3-d space motions of $>$100,000 stars in the Andromeda halo, including (by combining individual measurements) robust space motions of its entire globular cluster and most of its dwarf galaxy satellite populations. It will also identify high-velocity stars in Andromeda, providing unique information on the processes that create this population. These data offer a unique opportunity to study the immigration history, halo formation, and underlying dark matter scaffolding of a galaxy other than our own. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: Submitted in response to the call for Roman Space Telescope Core Community Survey white papers

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)

Clusters, Clouds, and Correlations: Relating Young Clusters to Giant Molecular Clouds in M33 and M31

Authors: Joshua Peltonen, Erik Rosolowsky, L. Clifton Johnson, Anil C. Seth, Julianne Dalcanton, Eric F. Bell, Jonathan Braine, Eric W. Koch, Margaret Lazzarini, Adam K. Leroy, Evan D. Skillman, Adam Smercina, Tobin Wainer, Benjamin F. Williams

Abstract: We use young clusters and giant molecular clouds (GMCs) in the galaxies M33 and M31 to constrain temporal and spatial scales in the star formation process. In M33, we compare the PHATTER catalogue of 1214 clusters with ages measured via colour-magnitude diagram (CMD) fitting to 444 GMCs identified from a new 35 pc resolution ALMA $^{12}$CO(2-1) survey. In M31, we compare the PHAT catalogue of 1249… ▽ More We use young clusters and giant molecular clouds (GMCs) in the galaxies M33 and M31 to constrain temporal and spatial scales in the star formation process. In M33, we compare the PHATTER catalogue of 1214 clusters with ages measured via colour-magnitude diagram (CMD) fitting to 444 GMCs identified from a new 35 pc resolution ALMA $^{12}$CO(2-1) survey. In M31, we compare the PHAT catalogue of 1249 clusters to 251 GMCs measured from a CARMA $^{12}$CO(1-0) survey with 20 pc resolution. Through two-point correlation analysis, we find that young clusters have a high probability of being near other young clusters, but correlation between GMCs is suppressed by the cloud identification algorithm. By comparing the positions, we find that younger clusters are closer to GMCs than older clusters. Through cross-correlation analysis of the M33 cluster data, we find that clusters are statistically associated when they are $\leq$10 Myr old. Utilizing the high precision ages of the clusters, we find that clusters older than $\approx 18$ Myr are uncorrelated with the molecular ISM. Using the spatial coincidence of the youngest clusters and GMCs in M33, we estimate that clusters spend $\approx$4-6 Myr inside their parent GMC. Through similar analysis, we find that the GMCs in M33 have a total lifetime of $\approx 11$-15 Myr. We also develop a drift model and show that the above correlations can be explained if the clusters in M33 have a 5-10 km s$^{-1}$ velocity dispersion relative to the molecular ISM. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: 14 pages, 13 figures, 1 tables, accepted for publication at MNRAS

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)

Wide-field CO isotopologue emission and the CO-to-H$_2$ factor across the nearby spiral galaxy M101

Authors: Jakob S. den Brok, Frank Bigiel, Jérémy Chastenet, Karin Sandstrom, Adam Leroy, Antonio Usero, Eva Schinnerer, Erik W. Rosolowsky, Eric W. Koch, I-Da Chiang, Ashley T. Barnes, Johannes Puschnig, Toshiki Saito, Ivana Bešlić, Melanie Chevance, Daniel A. Dale, Cosima Eibensteiner, Simon Glover, María J. Jiménez-Donaire, Yu-Hsuan Teng, Thomas G. Williams

Abstract: Carbon monoxide (CO) emission is the most widely used tracer of the bulk molecular gas in the interstellar medium (ISM) in extragalactic studies. The CO-to-H$_2$ conversion factor, $α_{\rm CO}$, links the observed CO emission to the total molecular gas mass. However, no single prescription perfectly describes the variation of $α_{\rm CO}$ across all environments across galaxies as a function of me… ▽ More Carbon monoxide (CO) emission is the most widely used tracer of the bulk molecular gas in the interstellar medium (ISM) in extragalactic studies. The CO-to-H$_2$ conversion factor, $α_{\rm CO}$, links the observed CO emission to the total molecular gas mass. However, no single prescription perfectly describes the variation of $α_{\rm CO}$ across all environments across galaxies as a function of metallicity, molecular gas opacity, line excitation, and other factors. Using resolved spectral line observations of CO and its isotopologues, we can constrain the molecular gas conditions and link them to a variation in the conversion factor. We present new IRAM 30-m 1mm and 3mm line observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O} across the nearby galaxy M101. Based on the CO isotopologue line ratios, we find that selective nucleosynthesis and opacity changes are the main drivers of the variation in the line emission across the galaxy. Furthermore, we estimated $α_{\rm CO(1-0)}$ using different approaches, including (i) the dust mass surface density derived from far-IR emission as an independent tracer of the total gas surface density and (ii) LTE-based measurements using the optically thin $^{13}$CO(1-0) intensity. We find an average value of $α_{\rm CO}=4.4{\pm}0.9\rm\,M_\odot\,pc^{-2}(K\,km\,s^{-1})^{-1}$ across the galaxy, with a decrease by a factor of 10 toward the 2 kpc central region. In contrast, we find LTE-based values are lower by a factor of 2-3 across the disk relative to the dust-based result. Accounting for $α_{\rm CO}$ variations, we found significantly reduced molecular gas depletion time by a factor 10 in the galaxy's center. In conclusion, our result suggests implications for commonly derived scaling relations, such as an underestimation of the slope of the Kennicutt Schmidt law, if $α_{\rm CO}$ variations are not accounted for. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: Accepted for publication in A&A, 25 pages, 15 figures

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

Kpc-scale properties of dust temperature in terms of dust mass and star formation activity

Authors: I-Da Chiang, Hiroyuki Hirashita, Jeremy Chastenet, Eric W. Koch, Adam K. Leroy, Erik Rosolowsky, Karin M. Sandstrom, Amy Sardone, Jiayi Sun, Thomas G. Williams

Abstract: We investigate how the dust temperature is affected by local environmental quantities, especially dust surface density ($Σ_\mathrm{dust}$), dust-to-gas ratio (D/G) and interstellar radiation field. We compile multi-wavelength observations in 46 nearby galaxies, uniformly processed with a common physical resolution of $2~$kpc. A physical dust model is used to fit the infrared dust emission spectral… ▽ More We investigate how the dust temperature is affected by local environmental quantities, especially dust surface density ($Σ_\mathrm{dust}$), dust-to-gas ratio (D/G) and interstellar radiation field. We compile multi-wavelength observations in 46 nearby galaxies, uniformly processed with a common physical resolution of $2~$kpc. A physical dust model is used to fit the infrared dust emission spectral energy distribution (SED) observed with WISE and Herschel. The star formation rate (SFR) is traced with GALEX ultraviolet data corrected by WISE infrared. We find that the dust temperature correlates well with the SFR surface density ($Σ_{\rm SFR}$), which traces the radiation from young stars. The dust temperature decreases with increasing D/G at fixed $Σ_{\rm SFR}$ as expected from stronger dust shielding at high D/G, when $Σ_\mathrm{SFR}$ is higher than $\sim 2\times 10^{-3}~\rm M_\odot~yr^{-1}~kpc^{-2}$. These measurements are in good agreement with the dust temperature predicted by our proposed analytical model. Below this range of $Σ_\mathrm{SFR}$, the observed dust temperature is higher than the model prediction and is only weakly dependent on D/G, which is possibly due to the dust heating from old stellar population or the variation of SFR within the past $10^{10}~$yr. Overall, the dust temperature as a function of $Σ_\mathrm{SFR}$ and $Σ_\mathrm{dust}$ predicted by our analytical model is consistent with observations. We also notice that at fixed gas surface density, $Σ_{\rm SFR}$ tends to increase with D/G, i.e. we can empirically modify the Kennicutt-Schmidt law with a dependence on D/G to better match observations. △ Less

Submitted 1 February, 2023; originally announced February 2023.

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

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

Cloud-scale Radio Surveys of Star Formation and Feedback in Triangulum Galaxy M 33: VLA Observations

Authors: F. S. Tabatabaei, W. Cotton, E. Schinnerer, R. Beck, A. Brunthaler, K. M. Menten, J. Braine, E. Corbelli, C. Kramer, J. E. Beckman, J. H. Knapen, R. Paladino, E. Koch, A. Camps Farina

Abstract: Studying the interplay between massive star formation and the interstellar medium (ISM) is paramount to understand the evolution of galaxies. Radio continuum (RC) emission serves as an extinction-free tracer of both massive star formation and the energetic components of the interstellar medium. We present a multi-band radio continuum survey of the local group galaxy M 33 down to ~30 pc linear reso… ▽ More Studying the interplay between massive star formation and the interstellar medium (ISM) is paramount to understand the evolution of galaxies. Radio continuum (RC) emission serves as an extinction-free tracer of both massive star formation and the energetic components of the interstellar medium. We present a multi-band radio continuum survey of the local group galaxy M 33 down to ~30 pc linear resolution observed with the Karl G. Jansky Very Large Array (VLA). We calibrate the star-formation rate surface density and investigate the impact of diffuse emission on this calibration using a structural decomposition. Separating the thermal and nonthermal emission components, the correlation between different phases of the interstellar medium and the impact of massive star formation are also investigated. Radio sources with sizes <~ 200 pc constitute about 36% (46%) of the total RC emission at 1.5 GHz (6.3 GHz) in the inner 18' x 18' (or 4kpc x 4kpc) disk of M 33. The nonthermal spectral index becomes flatter with increasing star-formation rate surface density, indicating the escape of cosmic ray electrons {from their birth places}. The magnetic field strength also increases with star-formation rate following a bi-modal relation, indicating that the small-scale turbulent dynamo acts more efficiently at higher luminosities and star-formation rates. Although the correlations are tighter in star-forming regions, the nonthermal emission is correlated also with the more quiescent molecular gas in the ISM. An almost linear molecular star-formation law exists in M 33 when excluding diffuse structures. Massive star formation amplifies the magnetic field and increases the number of high-energy cosmic ray electrons, which can help the onset of winds and outflows. △ Less

Submitted 3 September, 2022; originally announced September 2022.

Comments: 20 pages, 12 figures, accepted for publication in MNRAS journal

Testing the Momentum-driven Supernova Feedback Paradigm in M31

Authors: Sumit K. Sarbadhicary, Davide Martizzi, Enrico Ramirez-Ruiz, Eric Koch, Katie Auchettl, Carles Badenes, Laura Chomiuk

Abstract: Momentum feedback from isolated supernova remnants (SNRs) have been increasingly recognized by modern cosmological simulations as a resolution-independent means to implement the effects of feedback in galaxies, such as turbulence and winds. However, the integrated momentum yield from SNRs is uncertain due to the effects of SN clustering and interstellar medium (ISM) inhomogeneities. In this paper,… ▽ More Momentum feedback from isolated supernova remnants (SNRs) have been increasingly recognized by modern cosmological simulations as a resolution-independent means to implement the effects of feedback in galaxies, such as turbulence and winds. However, the integrated momentum yield from SNRs is uncertain due to the effects of SN clustering and interstellar medium (ISM) inhomogeneities. In this paper, we use spatially-resolved observations of the prominent 10-kpc star-forming ring of M31 to test models of mass-weighted ISM turbulence driven by momentum feedback from isolated, non-overlapping SNRs. We use a detailed stellar-age distribution (SAD) map from the Panchromatic Hubble Andromeda Treasury (PHAT) survey, observationally-constrained SN delay-time distributions, and maps of the atomic and molecular hydrogen to estimate the mass-weighted velocity dispersion using the Martizzi et al. ISM turbulence model. Our estimates are within a factor of 2 of the observed mass-weighted velocity dispersion in most of the ring, but exceed observations at densities $\lesssim 0.2$ cm$^{-3}$ and SN rates $>2.1\times 10^{-4}$ SN yr$^{-1}$ kpc$^{-2}$, even after accounting for plausible variations in stellar-age distribution models and ISM scale height assumptions. We conclude that at high SN rates the momentum deposited is most likely suppressed by the non-linear effects of SN clustering, while at low densities, SNRs reach pressure equilibrium before the cooling phase. These corrections should be introduced in models of momentum-driven feedback and ISM turbulence. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 14 pages, 7 figures, published in ApJ

Journal ref: March 2022, The Astrophysical Journal, Volume 928, Issue 1, id.54, 10 pp

The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package

Authors: The Astropy Collaboration, Adrian M. Price-Whelan, Pey Lian Lim, Nicholas Earl, Nathaniel Starkman, Larry Bradley, David L. Shupe, Aarya A. Patil, Lia Corrales, C. E. Brasseur, Maximilian Nöthe, Axel Donath, Erik Tollerud, Brett M. Morris, Adam Ginsburg, Eero Vaher, Benjamin A. Weaver, James Tocknell, William Jamieson, Marten H. van Kerkwijk, Thomas P. Robitaille, Bruce Merry, Matteo Bachetti, H. Moritz Günther, Thomas L. Aldcroft , et al. (111 additional authors not shown)

Abstract: The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as… ▽ More The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates for the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: 43 pages, 5 figures. To appear in ApJ. The author list has two parts: the authors that made significant contributions to the writing and/or coordination of the paper, followed by maintainers of and contributors to the Astropy Project. The position in the author list does not correspond to contributions to the Astropy Project as a whole

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)

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!

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

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

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

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

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

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

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

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

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

Submitted 30 September, 2021; originally announced October 2021.

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

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

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

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

Submitted 23 September, 2021; originally announced September 2021.

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

On the Three-Dimensional Structure of Local Molecular Clouds

Authors: Catherine Zucker, Alyssa Goodman, João Alves, Shmuel Bialy, Eric W. Koch, Joshua S. Speagle, Michael M. Foley, Douglas Finkbeiner, Reimar Leike, Torsten Enßlin, Joshua E. G. Peek, Gordian Edenhofer

Abstract: We leverage the 1 pc spatial resolution of the Leike et al. 2020 3D dust map to characterize the three-dimensional structure of nearby molecular clouds ($d \lesssim 400$ pc). We start by "skeletonizing" the clouds in 3D volume density space to determine their "spines," which we project on the sky to constrain cloud distances with $\approx 1\%$ uncertainty. For each cloud, we determine an average r… ▽ More We leverage the 1 pc spatial resolution of the Leike et al. 2020 3D dust map to characterize the three-dimensional structure of nearby molecular clouds ($d \lesssim 400$ pc). We start by "skeletonizing" the clouds in 3D volume density space to determine their "spines," which we project on the sky to constrain cloud distances with $\approx 1\%$ uncertainty. For each cloud, we determine an average radial volume density profile around its 3D spine and fit the profiles using Gaussian and Plummer functions. The radial volume density profiles are well-described by a two-component Gaussian function, consistent with clouds having broad, lower-density outer envelopes and narrow, higher-density inner layers. The ratio of the outer to inner envelope widths is $\approx 3:1$. We hypothesize that these two components may be tracing a transition between atomic and diffuse molecular gas or between the unstable and cold neutral medium. Plummer-like models can also provide a good fit, with molecular clouds exhibiting shallow power-law wings with density, $n$, falling off like $n^{-2}$ at large radii. Using Bayesian model selection, we find that parameterizing the clouds' profiles using a single Gaussian is disfavored. We compare our results with 2D dust extinction maps, finding that the 3D dust recovers the total cloud mass from integrated approaches with fidelity, deviating only at higher levels of extinction ($A_V \gtrsim 2 - 3$ mag). The 3D cloud structure described here will enable comparisons with synthetic clouds generated in simulations, offering unprecedented insight into the origins and fates of molecular clouds in the interstellar medium. △ Less

Submitted 27 September, 2021; v1 submitted 20 September, 2021; originally announced September 2021.

Comments: Published in ApJ. For data behind the figures see https://dataverse.harvard.edu/dataverse/3D_cloud_structure . For interactive gallery of local clouds see https://faun.rc.fas.harvard.edu/czucker/Paper_Figures/3D_Cloud_Topologies/gallery.html . Companion paper to Bialy et al. 2021 (arXiv:2109.09763). For press see https://www.cfa.harvard.edu/news/gigantic-cavity-space-sheds-new-light-how-stars-form

Journal ref: ApJ 919 35 (2021)

Local Environments of Low-Redshift Supernovae

Authors: Serena A. Cronin, Dyas Utomo, Adam K. Leroy, Erica A. Behrens, Jeremy Chastenet, Tyler Holland-Ashford, Eric W. Koch, Laura A. Lopez, Karin M. Sandstrom, Thomas G. Williams

Abstract: We characterize the local (2-kpc sized) environments of Type Ia, II, and Ib/c supernovae (SNe) that have recently occurred in nearby ($d\lesssim50$ Mpc) galaxies. Using ultraviolet (UV, from GALEX) and infrared (IR, from WISE) maps of 359 galaxies and a sample of 472 SNe, we measure the star formation rate surface density ($Σ_{\rm SFR}$) and stellar mass surface density ($Σ_\star$) in a 2-kpc beam… ▽ More We characterize the local (2-kpc sized) environments of Type Ia, II, and Ib/c supernovae (SNe) that have recently occurred in nearby ($d\lesssim50$ Mpc) galaxies. Using ultraviolet (UV, from GALEX) and infrared (IR, from WISE) maps of 359 galaxies and a sample of 472 SNe, we measure the star formation rate surface density ($Σ_{\rm SFR}$) and stellar mass surface density ($Σ_\star$) in a 2-kpc beam centered on each SN site. We show that core-collapse SNe are preferentially located along the resolved galactic star-forming main sequence, whereas Type Ia SNe are extended to lower values of $Σ_{\rm SFR}$ at fixed $Σ_\star$, indicative of locations inside quiescent galaxies or quiescent regions of galaxies. We also test how well the radial distribution of each SN type matches the radial distributions of UV and IR light in each host galaxy. We find that, to first order, the distributions of all types of SNe mirror that of both near-IR light (3.4 and 4.5 microns, tracing the stellar mass distribution) and mid-IR light (12 and 22 microns, tracing emission from hot, small grains), and also resemble our best-estimate $Σ_{\rm SFR}$. All types of SNe appear more radially concentrated than the UV emission of their host galaxies. In more detail, the distributions of Type II SNe show small statistical differences from that of near-IR light. We attribute this overall structural uniformity to the fact that within any individual galaxy, $Σ_{\rm SFR}$ and $Σ_\star$ track one another well, with variations in $Σ_{\rm SFR}/Σ_\star$ most visible when comparing between galaxies. △ Less

Submitted 15 September, 2021; originally announced September 2021.

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

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

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

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

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

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

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

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

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

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

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

PHANGS-ALMA Data Processing and Pipeline

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

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

Submitted 14 April, 2021; originally announced April 2021.

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

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

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

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

Submitted 6 April, 2021; originally announced April 2021.

Comments: Accepted in MNRAS. 25 pg, 13 figures

Benchmarking Dust Emission Models in M101

Authors: Jeremy Chastenet, Karin Sandstrom, I-Da Chiang, Brandon S. Hensley, Bruce T. Draine, Karl D. Gordon, Eric W. Koch, Adam K. Leroy, Dyas Utomo, Thomas G. Williams

Abstract: We present a comparative study of four physical dust models and two single-temperature modified blackbody models by fitting them to the resolved WISE, Spitzer, and Herschel photometry of M101 (NGC 5457). Using identical data and a grid-based fitting technique, we compare the resulting dust and radiation field properties derived from the models. We find that the dust mass yielded by the different m… ▽ More We present a comparative study of four physical dust models and two single-temperature modified blackbody models by fitting them to the resolved WISE, Spitzer, and Herschel photometry of M101 (NGC 5457). Using identical data and a grid-based fitting technique, we compare the resulting dust and radiation field properties derived from the models. We find that the dust mass yielded by the different models can vary by up to factor of 3 (factor of 1.4 between physical models only), although the fits have similar quality. Despite differences in their definition of the carriers of the mid-IR aromatic features, all physical models show the same spatial variations for the abundance of that grain population. Using the well determined metallicity gradient in M101 and resolved gas maps, we calculate an approximate upper limit on the dust mass as a function of radius. All physical dust models are found to exceed this maximum estimate over some range of galactocentric radii. We show that renormalizing the models to match the same Milky Way high latitude cirrus spectrum and abundance constraints can reduce the dust mass differences between models and bring the total dust mass below the maximum estimate at all radii. △ Less

Submitted 23 February, 2021; v1 submitted 22 January, 2021; originally announced January 2021.

Comments: Accepted for publication in ApJ -- 28 pp (12 figures) + appendices

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

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

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

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

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

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) I. Ultraviolet to Infrared Photometry of 22 Million Stars in M33

Authors: Benjamin F. Williams, Meredith J. Durbin, Julianne J. Dalcanton, Dustin Lang, Leo Girardi, Adam Smercina, Andrew Dolphin, Daniel R. Weisz, Yumi Choi, Eric F. Bell, Erik Rosolowsky, Evan Skillman, Eric W. Koch, Christine W. Lindberg, Lea Hagen, Karl D. Gordon, Anil Seth, Karoline Gilbert, Puragra Guhathakurta, Tod Lauer, Luciana Bianchi

Abstract: We present panchromatic resolved stellar photometry for 22 million stars in the Local Group dwarf spiral Triangulum (M33), derived from Hubble Space Telescope (HST) observations with the Advanced Camera for Surveys (ACS) in the optical (F475W, F814W), and the Wide Field Camera 3 (WFC3) in the near ultraviolet (F275W, F336W) and near-infrared (F110W, F160W) bands. The large, contiguous survey area… ▽ More We present panchromatic resolved stellar photometry for 22 million stars in the Local Group dwarf spiral Triangulum (M33), derived from Hubble Space Telescope (HST) observations with the Advanced Camera for Surveys (ACS) in the optical (F475W, F814W), and the Wide Field Camera 3 (WFC3) in the near ultraviolet (F275W, F336W) and near-infrared (F110W, F160W) bands. The large, contiguous survey area covers $\sim$14 square kpc and extends to 3.5 kpc (14 arcmin, or 1.5-2 scale lengths) from the center of M33. The PHATTER observing strategy and photometry technique closely mimic those of the Panchromatic Hubble Andromeda Treasury (PHAT), but with updated photometry techniques that take full advantage of all overlapping pointings (aligned to within $<$5-10 milliarcseconds) and improved treatment of spatially-varying point spread functions. The photometry reaches a completeness-limited depth of F475W$\sim$28.5 in the lowest surface density regions observed in M33 and F475W$\sim$26.5 in the most crowded regions found near the center of M33. We find the young populations trace several relatively tight arms, while the old populations show a clear, looser two-armed structure. We present extensive analysis of the data quality including artificial star tests to quantify completeness, photometric uncertainties, and flux biases. This stellar catalog is the largest ever produced for M33, and is publicly available for download by the community. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Comments: 38 pages, 6 tables, 25 figures, accepted for publication in ApJS

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

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

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

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

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

CHIMPS2: Survey description and $^{12}$CO emission in the Galactic Centre

Authors: D. J. Eden, T. J. T. Moore, M. J. Currie, A. J. Rigby, E. Rosolowsky, Y. Su, Kee-Tae Kim, H. Parsons, O. Morata, H. -R. Chen, T. Minamidani, Geumsook Park, S. E. Ragan, J. S. Urquhart, R. Rani, K. Tahani, S. J. Billington, S. Deb, C. Figura, T. Fujiyoshi, G. Joncas, L. W. Liao, T. Liu, H. Ma, P. Tuan-Anh , et al. (81 additional authors not shown)

Abstract: The latest generation of Galactic-plane surveys is enhancing our ability to study the effects of galactic environment upon the process of star formation. We present the first data from CO Heterodyne Inner Milky Way Plane Survey 2 (CHIMPS2). CHIMPS2 is a survey that will observe the Inner Galaxy, the Central Molecular Zone (CMZ), and a section of the Outer Galaxy in $^{12}$CO, $^{13}$CO, and C… ▽ More The latest generation of Galactic-plane surveys is enhancing our ability to study the effects of galactic environment upon the process of star formation. We present the first data from CO Heterodyne Inner Milky Way Plane Survey 2 (CHIMPS2). CHIMPS2 is a survey that will observe the Inner Galaxy, the Central Molecular Zone (CMZ), and a section of the Outer Galaxy in $^{12}$CO, $^{13}$CO, and C$^{18}$O $(J = 3\rightarrow2)$ emission with the Heterodyne Array Receiver Program on the James Clerk Maxwell Telescope (JCMT). The first CHIMPS2 data presented here are a first look towards the CMZ in $^{12}$CO J = 3$\rightarrow$2 and cover $-3^{\circ}\leq\,\ell\,\leq\,5^{\circ}$ and $\mid$b$\mid \leq 0.5^{\circ}$ with angular resolution of 15 arcsec, velocity resolution of 1 km s$^{-1}$, and rms $ΔT_A ^\ast =$ 0.58 K at these resolutions. Such high-resolution observations of the CMZ will be a valuable data set for future studies, whilst complementing the existing Galactic Plane surveys, such as SEDIGISM, the Herschel infrared Galactic Plane Survey, and ATLASGAL. In this paper, we discuss the survey plan, the current observations and data, as well as presenting position-position maps of the region. The position-velocity maps detect foreground spiral arms in both absorption and emission. △ Less

Submitted 10 September, 2020; originally announced September 2020.

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