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

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, et al.

► To cite this version:

Karin M. Sandstrom, Eric W. Koch, Adam K. Leroy, Erik Rosolowsky, Eric Emsellem, et al. PHANGS-JWST First Results: Tracing the Diffuse Interstellar Medium with JWST Imaging of Polycyclic Aromatic Hydrocarbon Emission in Nearby Galaxies. The Astrophysical journal letters, 2023, 944, 10.3847/2041-8213/aca972 . insu-04473228

HAL Id: insu-04473228

https://insu.hal.science/insu-04473228

Submitted on 22 Feb 2024

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# PHANGS-JWST First Results: Tracing the Diffuse Interstellar Medium with JWST Imaging of Polycyclic Aromatic Hydrocarbon Emission in Nearby Galaxies

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, et al.

► To cite this version:

Karin M. Sandstrom, Eric W. Koch, Adam K. Leroy, Erik Rosolowsky, Eric Emsellem, et al. PHANGS-JWST First Results: Tracing the Diffuse Interstellar Medium with JWST Imaging of Polycyclic Aromatic Hydrocarbon Emission in Nearby Galaxies. The Astrophysical Journal Letters, 2023, 944, 10.3847/2041-8213/aca972. insu-04473228

HAL Id: insu-04473228

https://insu.hal.science/insu-04473228v1

Submitted on 22 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# PHANGS–JWST First Results: Tracing the Diffuse Interstellar Medium with JWST Imaging of Polycyclic Aromatic Hydrocarbon Emission in Nearby Galaxies

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²⁸, Kathryn Grasha²⁹,³⁰,

Brent Groves³¹, Hamid Hassani⁴, Jonathan D. Henshaw¹⁰,³², Annie Hughes³³, Jaeyeon Kim³⁴, Ralf S. Klessen²⁸,³⁵,

Kathryn Kreckel⁸, J. M. Diederik Kruijssen²⁴, Laura A. Lopez³⁶,³⁷,³⁸, Daizhong Liu²¹, Sharon E. Meidt²²,

Eric J. Murphy³⁹, Hsi-An Pan⁴⁰, Miguel Querejeta⁴¹, Toshiki Saito⁴², Amy Sardone³,³⁷, Mattia C. Sormani⁴³,

Jessica Sutter¹, Antonio Usero⁴¹, and Elizabeth J. Watkins⁸

¹Center for Astrophysics & Space Sciences, Department of Physics, University of California, San Diego, 9500 Gilman Drive, San Diego, CA 92093, USA
kmsandstrom@ucsd.edu

²Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA

³Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA

⁴Department of Physics, University of Alberta, Edmonton, AB, T6G 2E1, Canada

⁵European Southern Observatory, Karl-Schwarzschild-Straße 2, D-85748 Garching, Germany

⁶Univ Lyon, Univ Lyon1, ENS de Lyon, CNRS, Centre de Recherche Astrophysique de Lyon UMR5574, F-69230 Saint-Genis-Laval, France

⁷Jodrell Bank Centre for Astrophysics, Department of Physics and Astronomy, University of Manchester, Oxford Road, Manchester, M13 9PL, UK

⁸Astronomisches Rechen-Institut, Zentrum für Astronomie der Universität Heidelberg, Mönchhofstraße 12-14, D-69120 Heidelberg, Germany

⁹Sub-department of Astrophysics, Department of Physics, University of Oxford, Keble Road, Oxford, OX1 3RH, UK

¹⁰Max-Planck-Institut für Astronomie, Königstuhl 17, D-69117, Heidelberg, Germany

¹¹AURA for the European Space Agency (ESA), Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA

¹²Gemini Observatory/NSF's NOIRLab, 950 N. Cherry Avenue, Tucson, AZ 85721, USA

¹³Steward Observatory, University of Arizona, 933 N. Cherry Avenue, Tucson, AZ 85721, USA

¹⁴Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, MD 21218, USA

¹⁵Argelander-Institut für Astronomie, Universität Bonn, Auf dem Hügel 71, D-53121, Bonn, Germany

¹⁶INAF—Arcetri Astrophysical Observatory, Largo E. Fermi 5, I-50125, Florence, Italy

¹⁷The Observatories of the Carnegie Institution for Science, 131 Santa Barbara St., Pasadena 91101, CA, USA

¹⁸Departamento de Astronomía, Universidad de Chile, Camino del Observatorio 1515, Las Condes, Santiago, Chile

¹⁹Department of Astronomy and Joint Space-Science Institute, University of Maryland, College Park, MD 20742, USA

²⁰Centro de Astronomía (CITEVA), Universidad de Antofagasta, Avenida Angamos 601, Antofagasta, Chile

²¹Max-Planck-Institut für Extraterrestrische Physik (MPE), Giessenbachstr. 1, D-85748 Garching, Germany

²²Sterrenkundig Observatory, Universiteit Gent, Krijgslaan 281 S9, B-9000 Gent, Belgium

²³Institut für Theoretische Astrophysik, Zentrum für Astronomie der Universität Heidelberg, Albert-Ueberle-Strasse 2, D-69120 Heidelberg, Germany

²⁴Cosmic Origins of Life (COOL) Research DAO, coolresearch.io

²⁵Institute of Astronomy and Astrophysics, Academia Sinica, No. 1, Section 4, Roosevelt Road, Taipei 10617, Taiwan

²⁶Department of Physics and Astronomy, University of Wyoming, Laramie, WY 82071, USA

²⁷Department of Physics, University of Connecticut, 196A Auditorium Road, Unit 3046, Storrs, CT 06269, USA

²⁸Universität Heidelberg, Zentrum für Astronomie, Institut für Theoretische Astrophysik, Albert-Ueberle-Straße 2, D-69120 Heidelberg, Germany

²⁹Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611, Australia

³⁰ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Australia

³¹International Centre for Radio Astronomy Research, University of Western Australia, 7 Fairway, Crawley, 6009 WA, Australia

³²Astrophysics Research Institute, Liverpool John Moores University, 146 Brownlow Hill, Liverpool, L3 5RF, UK

³³IRAP, Université de Toulouse, CNRS, CNES, UPS, (Toulouse), France

³⁴Zentrum für Astronomie der Universität Heidelberg, Institut für Theoretische Astrophysik, Albert-Ueberle-Str. 2, D-69120 Heidelberg, Germany

³⁵Universität Heidelberg, Interdisziplinäres Zentrum für Wissenschaftliches Rechnen, Im Neuenheimer Feld 205, D-69120 Heidelberg, Germany

³⁶Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA

³⁷Center for Cosmology and Astroparticle Physics, 191 West Woodruff Avenue, Columbus, OH 43210, USA

³⁸Flatiron Institute, Center for Computational Astrophysics, NY 10010, USA

³⁹National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903, USA

⁴⁰Department of Physics, Tamkang University, No.151, Yingzhuan Road, Tamsui District, New Taipei City 251301, Taiwan

⁴¹Observatorio Astronómico Nacional (IGN), C/Alfonso XII, 3, E-28014 Madrid, Spain

⁴²National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo, 181-8588, Japan

⁴³Universität Heidelberg, Zentrum für Astronomie, Institut für Theoretische Astrophysik, Albert-Ueberle-Str 2, D-69120 Heidelberg, Germany

Received 2022 October 22; revised 2022 December 1; accepted 2022 December 7; published 2023 February 16

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

<figure><img src="image_2.png" /></figure>

Original content from this work may be used under the terms of the [Creative Commons Attribution 4.0 licence](https://creativecommons.org/licenses/by/4.0/). Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
# PHANGS–JWST First Results: Tracing the Diffuse Interstellar Medium with JWST Imaging of Polycyclic Aromatic Hydrocarbon Emission in Nearby Galaxies

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²⁸, Kathryn Grasha²⁹,³⁰,

Brent Groves³¹, Hamid Hassani⁴, Jonathan D. Henshaw¹⁰,³², Annie Hughes³³, Jaeyeon Kim³⁴, Ralf S. Klessen²⁸,³⁵,

Kathryn Kreckel⁸, J. M. Diederik Kruijssen²⁴, Laura A. Lopez³⁶,³⁷,³⁸, Daizhong Liu²¹, Sharon E. Meidt²²,

Eric J. Murphy³⁹, Hsi-An Pan⁴⁰, Miguel Querejeta⁴¹, Toshiki Saito⁴², Amy Sardone³,³⁷, Mattia C. Sormani⁴³,

Jessica Sutter¹, Antonio Usero⁴¹, and Elizabeth J. Watkins⁸

¹Center for Astrophysics & Space Sciences, Department of Physics, University of California, San Diego, 9500 Gilman Drive, San Diego, CA 92093, USA
kmsandstrom@ucsd.edu

²Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA

³Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA

⁴Department of Physics, University of Alberta, Edmonton, AB, T6G 2E1, Canada

⁵European Southern Observatory, Karl-Schwarzschild-Straße 2, D-85748 Garching, Germany

⁶Univ Lyon, Univ Lyon1, ENS de Lyon, CNRS, Centre de Recherche Astrophysique de Lyon UMR5574, F-69230 Saint-Genis-Laval, France

⁷Jodrell Bank Centre for Astrophysics, Department of Physics and Astronomy, University of Manchester, Oxford Road, Manchester, M13 9PL, UK

⁸Astronomisches Rechen-Institut, Zentrum für Astronomie der Universität Heidelberg, Mönchhofstraße 12-14, D-69120 Heidelberg, Germany

⁹Sub-department of Astrophysics, Department of Physics, University of Oxford, Keble Road, Oxford, OX1 3RH, UK

¹⁰Max-Planck-Institut für Astronomie, Königstuhl 17, D-69117, Heidelberg, Germany

¹¹AURA for the European Space Agency (ESA), Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA

¹²Gemini Observatory/NSF's NOIRLab, 950 N. Cherry Avenue, Tucson, AZ 85721, USA

¹³Steward Observatory, University of Arizona, 933 N. Cherry Avenue, Tucson, AZ 85721, USA

¹⁴Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, MD 21218, USA

¹⁵Argelander-Institut für Astronomie, Universität Bonn, Auf dem Hügel 71, D-53121, Bonn, Germany

¹⁶INAF—Arcetri Astrophysical Observatory, Largo E. Fermi 5, I-50125, Florence, Italy

¹⁷The Observatories of the Carnegie Institution for Science, 131 Santa Barbara St., Pasadena 91101, CA, USA

¹⁸Departamento de Astronomía, Universidad de Chile, Camino del Observatorio 1515, Las Condes, Santiago, Chile

¹⁹Department of Astronomy and Joint Space-Science Institute, University of Maryland, College Park, MD 20742, USA

²⁰Centro de Astronomía (CITEVA), Universidad de Antofagasta, Avenida Angamos 601, Antofagasta, Chile

²¹Max-Planck-Institut für Extraterrestrische Physik (MPE), Giessenbachstr. 1, D-85748 Garching, Germany

²²Sterrenkundig Observatory, Universiteit Gent, Krijgslaan 281 S9, B-9000 Gent, Belgium

²³Institut für Theoretische Astrophysik, Zentrum für Astronomie der Universität Heidelberg, Albert-Ueberle-Strasse 2, D-69120 Heidelberg, Germany
²⁴Cosmic Origins of Life (COOL) Research DAO, coolresearch.io

²⁵Institute of Astronomy and Astrophysics, Academia Sinica, No. 1, Section 4, Roosevelt Road, Taipei 10617, Taiwan

²⁶Department of Physics and Astronomy, University of Wyoming, Laramie, WY 82071, USA

²⁷Department of Physics, University of Connecticut, 196A Auditorium Road, Unit 3046, Storrs, CT 06269, USA

²⁸Universität Heidelberg, Zentrum für Astronomie, Institut für Theoretische Astrophysik, Albert-Ueberle-Straße 2, D-69120 Heidelberg, Germany

²⁹Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611, Australia

³⁰ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Australia

³¹International Centre for Radio Astronomy Research, University of Western Australia, 7 Fairway, Crawley, 6009 WA, Australia

³²Astrophysics Research Institute, Liverpool John Moores University, 146 Brownlow Hill, Liverpool, L3 5RF, UK

³³IRAP, Université de Toulouse, CNRS, CNES, UPS, (Toulouse), France

³⁴Zentrum für Astronomie der Universität Heidelberg, Institut für Theoretische Astrophysik, Albert-Ueberle-Str. 2, D-69120 Heidelberg, Germany

³⁵Universität Heidelberg, Interdisziplinäres Zentrum für Wissenschaftliches Rechnen, Im Neuenheimer Feld 205, D-69120 Heidelberg, Germany

³⁶Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA

³⁷Center for Cosmology and Astroparticle Physics, 191 West Woodruff Avenue, Columbus, OH 43210, USA

³⁸Flatiron Institute, Center for Computational Astrophysics, NY 10010, USA

³⁹National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903, USA

⁴⁰Department of Physics, Tamkang University, No.151, Yingzhuan Road, Tamsui District, New Taipei City 251301, Taiwan

⁴¹Observatorio Astronómico Nacional (IGN), C/Alfonso XII, 3, E-28014 Madrid, Spain

⁴²National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo, 181-8588, Japan

⁴³Universität Heidelberg, Zentrum für Astronomie, Institut für Theoretische Astrophysik, Albert-Ueberle-Str. 2, D-69120 Heidelberg, Germany

Received 2022 October 22; revised 2022 December 1; accepted 2022 December 7; published 2023 February 16

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

<figure><img src="image_2.png" /></figure>

Original content from this work may be used under the terms of the [Creative Commons Attribution 4.0 licence](https://creativecommons.org/licenses/by/4.0/). Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
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 Mid-Infrared Instrument observations of IC 5332, NGC 628, NGC 1365, and NGC 7496 from the Physics at High Angular resolution in Nearby GalaxieS-JWST Treasury. Using scaling relationships calibrated in Leroy et al., scaled F1130W provides 10–40 pc resolution and 3σ sensitivity of $\Sigma_{\text{gas}} \sim 2 M_{\odot} \text{pc}^{-2}$. We characterize the surface densities of structures seen at $<7 M_{\odot} \text{pc}^{-2}$ in our targets, where we expect the gas to be H I-dominated. We highlight the existence of filaments, interarm emission, and holes in the diffuse ISM at these low surface densities. Below $\sim 10 M_{\odot} \text{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 the 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.

*Unified Astronomy Thesaurus concepts: Polycyclic aromatic hydrocarbons (1280); Interstellar atomic gas (833); Interstellar medium (847); Astronomical simulations (1857)*

# 1. Introduction

With currently available radio telescopes, observations of galaxies outside the Local Group do not have sufficient resolution and sensitivity to characterize the structure and column densities of the atomic-gas-dominated diffuse interstellar medium (ISM) on scales of less than hundreds of parsecs. Observations in the Milky Way show that at small spatial scales this phase of the ISM has structure and characteristics that reflect the interaction of turbulence, feedback, galactic dynamics, magnetic fields, gas-phase transitions, and more (Ehlerová & Palouš 2005; Kalberla & Kerp 2009; Suad et al. 2014; Li et al. 2021; Soler et al. 2022). The insights provided by studying the diffuse ISM in the Milky Way and Local Group at high resolution inform our understanding of key processes in ISM physics, but are limited by distance ambiguities in the Milky Way and the mainly nonrepresentative sample of galaxies available in the Local Group (i.e., a nearly quiescent massive spiral, two dwarf spirals, and several low-mass dwarf irregulars). It is therefore of great interest to characterize the properties of the diffuse ISM at $\lesssim 100$ pc resolution in other massive star-forming galaxies, but until the Square Kilometer Array or the next-generation Very Large Array (ngVLA) are operational this is not possible directly with H I 121 cm observations.

One approach to studying the diffuse ISM is to use dust as a tracer of gas. With a known dust-to-gas ratio (D/G), assumed to apply across ISM phases, a measured dust surface density can be converted to a gas surface density (e.g., Bot et al. 2007; Draine et al. 2007; Peretto & Fuller 2009; Leroy et al. 2011; Eales et al. 2012; Roman-Duval et al. 2014; Scoville et al. 2014; Groves et al. 2015). Such a technique does not generally overcome the resolution limitations, however, since most tracers of the dust surface density are observed at coarse resolution. For instance, far-IR images of galaxies from the Herschel Space Observatory ($\sim 11''$ at 160 µm) still have resolution lower than what can be achieved with Very Large Array (VLA) H I observations.

Observations of polycyclic aromatic hydrocarbon (PAH) emission from the JWST open up new pathways to probe the diffuse ISM at $\lesssim 100$ pc scales in nearby galaxies because of their high angular resolution and sensitivity. Tracing gas with PAH emission should be possible if (i) dust and PAHs are mixed with gas with approximately constant ratios of D/G and PAH-to-dust (also known as “PAH fraction”; Draine & Li 2007), and (ii) the PAHs are illuminated by an approximately constant or known radiation field intensity and

spectrum. Because PAHs experience single-photon heating over a wide range of radiation field intensities (Li & Draine 2001), their resulting emission intensity varies linearly with the incident radiation field. If that radiation field is known (or can be calibrated locally), the PAH emission then linearly traces the PAH surface density, which could then be related to gas mass with a PAH-to-gas ratio, as described above. Variations in the PAH fraction and PAH-to-gas ratio are active areas of research (including Leroy et al. 2023; Chastenet et al. 2023; Egorov et al. 2023, in this Issue). Outside the vicinity of H II regions, PAH abundance appears to be relatively constant in the cold atomic and molecular ISM, though dependent on metallicity (Chastenet et al. 2023). Future work using JWST observations will be critical to characterizing the PAH fraction across a variety of galactic environments.

For the radiation field, the diffuse ISM in particular is expected to be dominated by the average diffuse interstellar radiation field (ISRF), which we will refer to as the diffuse ISRF. In the Milky Way, the diffuse ISRF is found to be powered predominantly by B stars (Mathis et al. 1983) and does not vary strongly on small scales (though it does change radially in the Galaxy). The diffuse ISRF dominates dust heating outside the sphere of influence of star-forming regions ($\sim 100$ pc in the Magellanic Clouds; Lawton et al. 2010) and outside the shielding interiors of dense clouds. By isolating PAH emission in regions where dust heating is dominated by the diffuse ISRF we can hope to construct a tracer for gas column based on the PAH emission.

Previous generations of IR telescopes did not have the resolution to separate out the diffuse-ISRF-powered PAH emission spatially outside the Local Group. At the typical $\sim 100–1000$ pc spatial resolution of Spitzer in nearby galaxies ($D \sim 10$ Mpc), the sphere of influence of individual star-forming regions is not resolved, leaving the measurement of the fraction of PAH emission being powered by the diffuse ISRF reliant on modeling. However, IR spectral energy distribution (SED) modeling of the mid- to far-IR emission from galaxies suggests that a substantial fraction of the dust and PAH emission in galaxies arises from diffuse ISRF heating (e.g., Draine et al. 2007; Aniano et al. 2012, 2020; Galliano et al. 2018).

The idea of tracing diffuse ISM gas with PAH emission has a long heritage from studies of the high-latitude “cirrus” emission in the Milky Way. Observations with IRAS found that the PAH-dominated broadband emission at 12 µm was well correlated with H I at high Galactic latitudes (Boulanger et al. 1985; Puget et al. 1985; Boulanger & Perault 1988).
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 Mid-Infrared Instrument observations of IC 5332, NGC 628, NGC 1365, and NGC 7496 from the Physics at High Angular resolution in Nearby GalaxieS-JWST Treasury. Using scaling relationships calibrated in Leroy et al., scaled F1130W provides 10–40 pc resolution and 3σ sensitivity of $\Sigma_{\text{gas}} \sim 2 M_{\odot} \text{pc}^{-2}$. We characterize the surface densities of structures seen at $<7 M_{\odot} \text{pc}^{-2}$ in our targets, where we expect the gas to be H I-dominated. We highlight the existence of filaments, interarm emission, and holes in the diffuse ISM at these low surface densities. Below $\sim 10 M_{\odot} \text{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 the 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.

*Unified Astronomy Thesaurus concepts: Polycyclic aromatic hydrocarbons (1280); Interstellar atomic gas (833); Interstellar medium (847); Astronomical simulations (1857)*

# 1. Introduction

With currently available radio telescopes, observations of galaxies outside the Local Group do not have sufficient resolution and sensitivity to characterize the structure and column densities of the atomic-gas-dominated diffuse interstellar medium (ISM) on scales of less than hundreds of parsecs. Observations in the Milky Way show that at small spatial scales this phase of the ISM has structure and characteristics that reflect the interaction of turbulence, feedback, galactic dynamics, magnetic fields, gas-phase transitions, and more (Ehlerová & Palouš 2005; Kalberla & Kerp 2009; Suad et al. 2014; Li et al. 2021; Soler et al. 2022). The insights provided by studying the diffuse ISM in the Milky Way and Local Group at high resolution inform our understanding of key processes in ISM physics, but are limited by distance ambiguities in the Milky Way and the mainly nonrepresentative sample of galaxies available in the Local Group (i.e., a nearly quiescent massive spiral, two dwarf spirals, and several low-mass dwarf irregulars). It is therefore of great interest to characterize the properties of the diffuse ISM at $\lesssim 100$ pc resolution in other massive star-forming galaxies, but until the Square Kilometer Array or the next-generation Very Large Array (ngVLA) are operational this is not possible directly with H I 121 cm observations.

One approach to studying the diffuse ISM is to use dust as a tracer of gas. With a known dust-to-gas ratio (D/G), assumed to apply across ISM phases, a measured dust surface density can be converted to a gas surface density (e.g., Bot et al. 2007; Draine et al. 2007; Peretto & Fuller 2009; Leroy et al. 2011; Eales et al. 2012; Roman-Duval et al. 2014; Scoville et al. 2014; Groves et al. 2015). Such a technique does not generally overcome the resolution limitations, however, since most tracers of the dust surface density are observed at coarse resolution. For instance, far-IR images of galaxies from the Herschel Space Observatory ($\sim 11''$ at 160 µm) still have resolution lower than what can be achieved with Very Large Array (VLA) H I observations.

Observations of polycyclic aromatic hydrocarbon (PAH) emission from the JWST open up new pathways to probe the diffuse ISM at $\lesssim 100$ pc scales in nearby galaxies because of their high angular resolution and sensitivity. Tracing gas with PAH emission should be possible if (i) dust and PAHs are mixed with gas with approximately constant ratios of D/G and PAH-to-dust (also known as “PAH fraction”; Draine & Li 2007), and (ii) the PAHs are illuminated by an approximately constant or known radiation field intensity and

spectrum. Because PAHs experience single-photon heating over a wide range of radiation field intensities (Li & Draine 2001), their resulting emission intensity varies linearly with the incident radiation field. If that radiation field is known (or can be calibrated locally), the PAH emission then linearly traces the PAH surface density, which could then be related to gas mass with a PAH-to-gas ratio, as described above. Variations in the PAH fraction and PAH-to-gas ratio are active areas of research (including Leroy et al. 2023; Chastenet et al. 2023; Egorov et al. 2023, in this Issue). Outside the vicinity of H II regions, PAH abundance appears to be relatively constant in the cold atomic and molecular ISM, though dependent on metallicity (Chastenet et al. 2023). Future work using JWST observations will be critical to characterizing the PAH fraction across a variety of galactic environments.

For the radiation field, the diffuse ISM in particular is expected to be dominated by the average diffuse interstellar radiation field (ISRF), which we will refer to as the diffuse ISRF. In the Milky Way, the diffuse ISRF is found to be powered predominantly by B stars (Mathis et al. 1983) and does not vary strongly on small scales (though it does change radially in the Galaxy). The diffuse ISRF dominates dust heating outside the sphere of influence of star-forming regions ($\sim 100$ pc in the Magellanic Clouds; Lawton et al. 2010) and outside the shielding interiors of dense clouds. By isolating PAH emission in regions where dust heating is dominated by the diffuse ISRF we can hope to construct a tracer for gas column based on the PAH emission.

Previous generations of IR telescopes did not have the resolution to separate out the diffuse-ISRF-powered PAH emission spatially outside the Local Group. At the typical $\sim 100–1000$ pc spatial resolution of Spitzer in nearby galaxies ($D \sim 10$ Mpc), the sphere of influence of individual star-forming regions is not resolved, leaving the measurement of the fraction of PAH emission being powered by the diffuse ISRF reliant on modeling. However, IR spectral energy distribution (SED) modeling of the mid- to far-IR emission from galaxies suggests that a substantial fraction of the dust and PAH emission in galaxies arises from diffuse ISRF heating (e.g., Draine et al. 2007; Aniano et al. 2012, 2020; Galliano et al. 2018).

The idea of tracing diffuse ISM gas with PAH emission has a long heritage from studies of the high-latitude “cirrus” emission in the Milky Way. Observations with IRAS found that the PAH-dominated broadband emission at 12 µm was well correlated with H I at high Galactic latitudes (Boulanger et al. 1985; Puget et al. 1985; Boulanger & Perault 1988).
<table><caption>**Table 1**
Galaxy Properties</caption><thead><tr><th>Target</th><th>Distance<br>(Mpc)</th><th>log M<sub>*</sub><br>(M<sub>☉</sub>)</th><th>SFR<sub>tot</sub><br>(M<sub>☉</sub> yr<sup>-1</sup>)</th><th>12 +log(O/H)<br>(dex)</th><th>i<br>(°)</th><th>R<sub>e</sub><br>(kpc)</th><th>F<sub>770W</sub><br>l<sub>CO,2-1</sub><br>(MJy sr<sup>-1</sup> (K km s<sup>-1</sup>)<sup>-1</sup>)</th><th>F<sub>1130W</sub><br>l<sub>CO,2-1</sub><br>(MJy sr<sup>-1</sup> (K km s<sup>-1</sup>)<sup>-1</sup>)</th></tr></thead><tbody><tr><td>IC 5332</td><td>9.0</td><td>9.7</td><td>0.4</td><td>8.30</td><td>27</td><td>3.6</td><td>0.79, 1.07</td><td></td></tr><tr><td>NGC 628</td><td>9.8</td><td>10.3</td><td>1.7</td><td>8.53</td><td>9</td><td>3.9</td><td>0.93, 1.29</td><td></td></tr><tr><td>NGC 1365</td><td>19.6</td><td>11.0</td><td>17</td><td>8.67</td><td>55</td><td>2.8</td><td>0.72, 1.07</td><td></td></tr><tr><td>NGC 7496</td><td>18.7</td><td>10.0</td><td>2.2</td><td>8.51</td><td>36</td><td>3.8</td><td>1.29, 1.78</td><td></td></tr></tbody></table>

**Note.** Properties adopted following Lee et al. (2023), which draws distances from Anand et al. (2021); stellar masses (log M*), star formation rates (SFRtot), and effective radius (R*), values from Leroy et al. (2021); inclinations (i) from Lang et al. (2020); and metallicities (12+log(O/H)) from Groves et al. (2023). Gas-phase metallicities are presented in the S-cal system (Pilyugin & Grebel 2016) estimated at the effective radius, R*. The F770W/lCO,2-1 and F1130W/lCO,2-1 values are taken from Leroy et al. (2023).

Observations with the AROME balloon experiment mapped the 3.3 µm PAH feature specifically and found it to be well correlated with the 12 µm emission from IRAS and H I (Giard et al. 1988, 1994). COBE/DIRBE observations from 3.5 to 240 µm allowed measurement of the mean spectrum of the H I-correlated high-latitude “cirrus” dust emission, which includes prominent PAH features (Dwek et al. 1997). Using Spitzer and Herschel observations, Compiègne et al. (2010) found that IRAC 8 µm emission linearly scaled with a combination of PAH abundance, ISM column density, and the intensity of the ISRF, as would be expected for grains undergoing single-photon heating. At this point, the association of PAH emission with diffuse H I at high latitudes in the Milky Way is well enough understood that it is a key constraint on essentially all dust models (see recent summary in Hensley et al. 2022).

In this Letter, we present a case study for using PAH emission observed with JWST to trace low-surface-density gas at 10–40 pc resolution in nearby galaxies. We show maps of the diffuse ISM at $\Sigma_{\text{gas}} < 7 M_{\odot} \text{ pc}^{-2}$, where we expect the gas to be primarily atomic (Krumholz et al. 2009; Sternberg et al. 2014; Bialy & Sternberg 2016). We identify low-surface-density gas structures in interarm regions, bubbles, and the outskirts of galaxies that are detected at high signal-to-noise ratio (S/N) and resolution in the JWST maps. While for this study we extrapolate a simple scaling determined in the H$_2$-dominated regions of the galaxy, in future work more detailed modeling will be needed to address the uncertainties inherent in the PAH emission to gas surface density conversion. It is clear, however, from these early JWST observations that the gas surface density sensitivity of Mid-Infrared Instrument (MIRI) observations in nearby galaxies far exceeds that of existing millimeter and radio observations.

We proceed as follows. In Section 2 we briefly describe the JWST-MIRI observations. We describe using PAH emission as a gas tracer in Section 3. In Section 4 we present maps of the gas surface in our target galaxies. We characterize the topology of the diffuse ISM using the *genus* statistic (Coles 1988; Koch et al. 2019) in Section 5. In Section 6 we discuss the implications of these findings for PAH emission as a gas tracer and highlight the utility of the extremely high-resolution view of the atomic ISM available from JWST that will not be equaled before the ngVLA or Square Kilometer Array come online.

## 2. Observations

We observed NGC 628, NGC 1365, NGC 7496, and IC 5332 with the MIRI imager between 2022 July 6 to August 13 using the F770W, F1000W, F1130W, and F2100W filters. Table 1 lists the relevant properties of our targets. F770W and F1130W are both

centered specifically on strong PAH features at 7.7 and 11.3 µm. Because of the relatively narrow widths of these filters ($\Delta\lambda \sim 2.2$ and 0.7 µm for F770W and F1130W, respectively), which are comparable or smaller than the widths of the 7.7 and 11.3 µm PAH features themselves, the fluxes in our targets are generally dominated by PAHs for typical radiation field intensities, with only a minimal contribution from the underlying hot dust continuum (see models from Draine & Li 2007, for example). Stellar continuum is only an important contribution to the F770W and F1130W bands toward specific mid-IR bright stars or low-ISM, high-stellar-surface-density regions like the nucleus of NGC 628 (Hoyer et al. 2023), which do not make up a significant fraction of the sight lines in our maps. Given the dominance of PAH emission at the F770W and F1130W filters in our targets, we do not attempt any removal of the hot dust continuum using F1000W (which may also include a strong or dominant contribution from PAH emission) or F2100W.

The observations and data reduction for the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) MIRI images are described in Lee et al. (2023). Validation of the background levels of the images, which is of particular importance for studying faint emission, follows the methods described in Appendices A and B of Leroy et al. (2023). Noise levels in the native resolution F770W and F1130W maps are $\sim 0.1-0.2$ MJy sr$^{-1}$.

## 3. Tracing the Diffuse Interstellar Medium with Polycyclic Aromatic Hydrocarbon Emission

### 3.1. Our Approach and Caveats

As discussed above, tracing gas surface density with PAH emission involves several assumptions: (i) a PAH-to-gas ratio (which can be separated into D/G and PAH fraction components), and (ii) a radiation field intensity and spectrum (in the regime where PAHs experience single-photon heating so the scaling with radiation field intensity is linear). In this Letter, we make two main assumptions that allow us use PAH emission to estimate the gas column density. These assumptions are likely too simplistic and should be greatly improved in future work, but they let us perform a pilot study for this approach to tracing gas.

First, we assume that the PAH emission to gas surface density ratio does not vary with cold gas phase. We can then pick regions where we know the gas surface density to empirically estimate this ratio. For this purpose, we select regions that are likely to be H$_2$-dominated where we also have an estimate of the H$_2$ surface density based on Atacama Large Millimeter/submillimeter Array (ALMA) $^{12}$CO J = (2–1)
<table><caption>**Table 1**
Galaxy Properties</caption><thead><tr><th>Target</th><th>Distance<br>(Mpc)</th><th>log M<sub>*</sub><br>(M<sub>☉</sub>)</th><th>SFR<sub>tot</sub><br>(M<sub>☉</sub> yr<sup>-1</sup>)</th><th>12 +log(O/H)<br>(dex)</th><th>i<br>(°)</th><th>R<sub>e</sub><br>(kpc)</th><th>F<sub>770W</sub><br>l<sub>CO,2-1</sub><br>(MJy sr<sup>-1</sup> (K km s<sup>-1</sup>)<sup>-1</sup>)</th><th>F<sub>1130W</sub><br>l<sub>CO,2-1</sub><br>(MJy sr<sup>-1</sup> (K km s<sup>-1</sup>)<sup>-1</sup>)</th></tr></thead><tbody><tr><td>IC 5332</td><td>9.0</td><td>9.7</td><td>0.4</td><td>8.30</td><td>27</td><td>3.6</td><td>0.79, 1.07</td><td></td></tr><tr><td>NGC 628</td><td>9.8</td><td>10.3</td><td>1.7</td><td>8.53</td><td>9</td><td>3.9</td><td>0.93, 1.29</td><td></td></tr><tr><td>NGC 1365</td><td>19.6</td><td>11.0</td><td>17</td><td>8.67</td><td>55</td><td>2.8</td><td>0.72, 1.07</td><td></td></tr><tr><td>NGC 7496</td><td>18.7</td><td>10.0</td><td>2.2</td><td>8.51</td><td>36</td><td>3.8</td><td>1.29, 1.78</td><td></td></tr></tbody></table>

Note. Properties adopted following Lee et al. (2023), which draws distances from Anand et al. (2021); stellar masses (log M*), star formation rates (SFRtot), and effective radius (R_e) values from Leroy et al. (2021); inclinations (i) from Lang et al. (2020); and metallicities (12+log(O/H)) from Groves et al. (2023). Gas-phase metallicities are presented in the S-cal system (Pilyugin & Grebel 2016) estimated at the effective radius, R_e. The F770W/lCO,2-1 and F1130W/lCO,2-1 values are taken from Leroy et al. (2023).

Observations with the AROME balloon experiment mapped the 3.3 µm PAH feature specifically and found it to be well correlated with the 12 µm emission from IRAS and H I (Giard et al. 1988, 1994). COBE/DIRBE observations from 3.5 to 240 µm allowed measurement of the mean spectrum of the H I-correlated high-latitude “cirrus” dust emission, which includes prominent PAH features (Dwek et al. 1997). Using Spitzer and Herschel observations, Compiègne et al. (2010) found that IRAC 8 µm emission linearly scaled with a combination of PAH abundance, ISM column density, and the intensity of the ISRF, as would be expected for grains undergoing single-photon heating. At this point, the association of PAH emission with diffuse H I at high latitudes in the Milky Way is well enough understood that it is a key constraint on essentially all dust models (see recent summary in Hensley et al. 2022).

In this Letter, we present a case study for using PAH emission observed with JWST to trace low-surface-density gas at 10–40 pc resolution in nearby galaxies. We show maps of the diffuse ISM at $\Sigma_{\text{gas}} < 7 M_{\odot} \text{ pc}^{-2}$, where we expect the gas to be primarily atomic (Krumholz et al. 2009; Sternberg et al. 2014; Bialy & Sternberg 2016). We identify low-surface-density gas structures in interarm regions, bubbles, and the outskirts of galaxies that are detected at high signal-to-noise ratio (S/N) and resolution in the JWST maps. While for this study we extrapolate a simple scaling determined in the H₂-dominated regions of the galaxy, in future work more detailed modeling will be needed to address the uncertainties inherent in the PAH emission to gas surface density conversion. It is clear, however, from these early JWST observations that the gas surface density sensitivity of Mid-Infrared Instrument (MIRI) observations in nearby galaxies far exceeds that of existing millimeter and radio observations.

We proceed as follows. In Section 2 we briefly describe the JWST-MIRI observations. We describe using PAH emission as a gas tracer in Section 3. In Section 4 we present maps of the gas surface in our target galaxies. We characterize the topology of the diffuse ISM using the *genus* statistic (Coles 1988; Koch et al. 2019) in Section 5. In Section 6 we discuss the implications of these findings for PAH emission as a gas tracer and highlight the utility of the extremely high-resolution view of the atomic ISM available from JWST that will not be equaled before the ngVLA or Square Kilometer Array come online.

## 2. Observations

We observed NGC 628, NGC 1365, NGC 7496, and IC 5332 with the MIRI imager between 2022 July 6 to August 13 using the F770W, F1000W, F1130W, and F2100W filters. Table 1 lists the relevant properties of our targets. F770W and F1130W are both

centered specifically on strong PAH features at 7.7 and 11.3 µm. Because of the relatively narrow widths of these filters ($\Delta\lambda \sim 2.2$ and 0.7 µm for F770W and F1130W, respectively), which are comparable or smaller than the widths of the 7.7 and 11.3 µm PAH features themselves, the fluxes in our targets are generally dominated by PAHs for typical radiation field intensities, with only a minimal contribution from the underlying hot dust continuum (see models from Draine & Li 2007, for example). Stellar continuum is only an important contribution to the F770W and F1130W bands toward specific mid-IR bright stars or low-ISM, high-stellar-surface-density regions like the nucleus of NGC 628 (Hoyer et al. 2023), which do not make up a significant fraction of the sight lines in our maps. Given the dominance of PAH emission at the F770W and F1130W filters in our targets, we do not attempt any removal of the hot dust continuum using F1000W (which may also include a strong or dominant contribution from PAH emission) or F2100W.

The observations and data reduction for the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) MIRI images are described in Lee et al. (2023). Validation of the background levels of the images, which is of particular importance for studying faint emission, follows the methods described in Appendices A and B of Leroy et al. (2023). Noise levels in the native resolution F770W and F1130W maps are $\sim 0.1-0.2$ MJy sr$^{-1}$.

## 3. Tracing the Diffuse Interstellar Medium with Polycyclic Aromatic Hydrocarbon Emission

### 3.1. Our Approach and Caveats

As discussed above, tracing gas surface density with PAH emission involves several assumptions: (i) a PAH-to-gas ratio (which can be separated into D/G and PAH fraction components), and (ii) a radiation field intensity and spectrum (in the regime where PAHs experience single-photon heating so the scaling with radiation field intensity is linear). In this Letter, we make two main assumptions that allow us use PAH emission to estimate the gas column density. These assumptions are likely too simplistic and should be greatly improved in future work, but they let us perform a pilot study for this approach to tracing gas.

First, we assume that the PAH emission to gas surface density ratio does not vary with cold gas phase. We can then pick regions where we know the gas surface density to empirically estimate this ratio. For this purpose, we select regions that are likely to be H₂-dominated where we also have an estimate of the H₂ surface density based on Atacama Large Millimeter/submillimeter Array (ALMA) $^{12}$CO $J=(2-1)$
<figure><img src="image_3.png" /><figcaption>**Figure 1.** IC 5332 gas surface density from scaled F1130W emission. A single black contour highlights the brightness level where the F1130W emission implies $\Sigma_{\text{gas}} = 7 M_{\odot} \text{ pc}^{-2}$. This figure and the subsequent figures for NGC 628, 1365, and 7496 (Figures 2, 3, and 4) are all shown with the same color bar. We show a 500 pc scale bar for reference.</figcaption></figure>

maps. To do this, we convert the CO emission to gas surface density with an assumed CO-to-H$_2$ conversion factor and find the ratio with PAH emission.

The second assumption we make is that the diffuse ISRF responsible for heating the dust is approximately constant outside the vicinity of star-forming regions. This allows us to apply the PAH emission to gas surface density ratio measured for the molecular gas to the rest of the ISM.

In practice, we adopt the median ratio of PAH emission to CO (2–1) emission measured by Leroy et al. (2023) for moderate mid-IR intensities ~0.5–30 MJy sr$^{-1}$ and 1"/7 ≈ 70–160 pc resolution in each target. Over this range of mid-IR intensities, the PAH tracing F770W and F1130W emission correlates strongly with CO, and an approximately linear relationship relates the PAH emission to CO emission.<fnref n="44" /> These PAH emission to CO ratios resemble those expected for typical diffuse ISRF values in inner-galaxy regions ($U \sim$ few $U_{\text{MMP}}^{45}$), PAH fractions ($q_{\text{PAH}} \sim$ 4%–5%), and dust-to-gas ratios (D/G

~ 0.01; see their Equations (4) and (6)). But we emphasize that the value we adopt is fundamentally empirical and anchored to the ALMA CO (2–1) data, not derived from models.

At the high resolution of our data we expect this median PAH to CO emission ratio to capture mostly diffuse ISM emission. Massive star-forming regions with high radiation fields are certainly present and associated with bright PAH and continuum emission (e.g., Chastenet et al. 2023; Egorov et al. 2023; Hassani et al. 2023), but as Figures 1 through 4 show, such bright regions do not make up most of the area in our maps and typically only dominate the dust heating on scales of ~100 pc (e.g., Lawton et al. 2010).

Adopting these two assumptions and this empirical ratio of PAH emission to CO emission, we calibrate a relationship between PAH emission and gas surface density. Then we extrapolate this scaling to the diffuse ISM, where we expect the gas to be predominantly atomic. We note that we do not mask star-forming regions in our resulting maps, but our focus on low-surface-density gas (or, equivalently, low-surface-brightness PAH emission) means such regions are generally not included in the diffuse gas we consider.

For the purpose of this pilot study, these simple assumptions and this fundamentally empirically approach allow a straightforward scaling of PAH emission to gas surface density, but a number of potential effects will need to be addressed in future

<footnote n="44">⁴⁴ The relationship between CO emission and PAH emission is not exactly linear. We adopt the median ratio for simplicity. If we had chosen to evaluate the power-law fits of CO to PAH emission in Leroy et al. (2023) at ~1 MJy sr$^{-1}$ instead, the ratio of PAH emission to CO emission would be ~20%–50% higher. However, at this low surface brightness, HI could already make important contributions to the total gas surface density.</footnote>

<footnote n="45">⁴⁵ $U_{\text{MMP}}$ is the ISRF intensity normalized to the Mathis et al. (1983) solar neighborhood ISRF.</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** IC 5332 gas surface density from scaled F1130W emission. A single black contour highlights the brightness level where the F1130W emission implies $\Sigma_{\text{gas}} = 7 M_{\odot} \text{ pc}^{-2}$. This figure and the subsequent figures for NGC 628, 1365, and 7496 (Figures 2, 3, and 4) are all shown with the same color bar. We show a 500 pc scale bar for reference.</figcaption></figure>

maps. To do this, we convert the CO emission to gas surface density with an assumed CO-to-H$_2$ conversion factor and find the ratio with PAH emission.

The second assumption we make is that the diffuse ISRF responsible for heating the dust is approximately constant outside the vicinity of star-forming regions. This allows us to apply the PAH emission to gas surface density ratio measured for the molecular gas to the rest of the ISM.

In practice, we adopt the median ratio of PAH emission to CO (2–1) emission measured by Leroy et al. (2023) for moderate mid-IR intensities ~0.5–30 MJy sr$^{-1}$ and 1"/7 ≈ 70–160 pc resolution in each target. Over this range of mid-IR intensities, the PAH tracing F770W and F1130W emission correlates strongly with CO, and an approximately linear relationship relates the PAH emission to CO emission.<fnref n="44" /> These PAH emission to CO ratios resemble those expected for typical diffuse ISRF values in inner-galaxy regions ($U \sim$ few $U_{\text{MMP}}^{45}$), PAH fractions ($q_{\text{PAH}} \sim$ 4%–5%), and dust-to-gas ratios (D/G

~ 0.01; see their Equations (4) and (6)). But we emphasize that the value we adopt is fundamentally empirical and anchored to the ALMA CO (2–1) data, not derived from models.

At the high resolution of our data we expect this median PAH to CO emission ratio to capture mostly diffuse ISM emission. Massive star-forming regions with high radiation fields are certainly present and associated with bright PAH and continuum emission (e.g., Chastenet et al. 2023; Egorov et al. 2023; Hassani et al. 2023), but as Figures 1 through 4 show, such bright regions do not make up most of the area in our maps and typically only dominate the dust heating on scales of ~100 pc (e.g., Lawton et al. 2010).

Adopting these two assumptions and this empirical ratio of PAH emission to CO emission, we calibrate a relationship between PAH emission and gas surface density. Then we extrapolate this scaling to the diffuse ISM, where we expect the gas to be predominantly atomic. We note that we do not mask star-forming regions in our resulting maps, but our focus on low-surface-density gas (or, equivalently, low-surface-brightness PAH emission) means such regions are generally not included in the diffuse gas we consider.

For the purpose of this pilot study, these simple assumptions and this fundamentally empirically approach allow a straightforward scaling of PAH emission to gas surface density, but a number of potential effects will need to be addressed in future

<footnote n="44">⁴⁴ The relationship between CO emission and PAH emission is not exactly linear. We adopt the median ratio for simplicity. If we had chosen to evaluate the power-law fits of CO to PAH emission in Leroy et al. (2023) at ~1 MJy sr$^{-1}$ instead, the ratio of PAH emission to CO emission would be ~20%–50% higher. However, at this low surface brightness, HI could already make important contributions to the total gas surface density.</footnote>

<footnote n="45">⁴⁵ $U_{\text{MMP}}$ is the ISRF intensity normalized to the Mathis et al. (1983) solar neighborhood ISRF.</footnote>
work to rigorously calibrate PAH emission as a tracer. There is evidence from measurements of depletion that D/Gs vary with ISM density (Jenkins 2009; Roman-Duval et al. 2021) by a factor of up to 2, though this depends on whether the observed depleted oxygen is incorporated into dust (Whittet 2010; Choban et al. 2022). Changes in the D/G between ISM phases at fixed metallicity have also been observed using far-IR, H I, and CO observations (Roman-Duval et al. 2014; Chiang et al. 2018; Vilchez et al. 2019). From high-spatial-resolution studies of Local Group galaxies (Chastenet et al. 2019) and the early PHANGS-JWST targets (Chastenet et al. 2023; Egorov et al. 2023), the PAH fraction appears to be relatively constant in the atomic and molecular phases, but does drop sharply in H II regions, and has a metallicity dependence. Milky Way observations have suggested that PAH abundance can vary between the warm neutral medium (WNM) and cold neutral medium (CNM) (Hensley et al. 2022). The diffuse ISRF in galaxies, judged from the $U_{min}$ component in far-IR SED modeling (e.g., Draine et al. 2007; Aniano et al. 2012), appears to vary smoothly and primarily radially in highly resolved Local Group galaxies (Draine et al. 2014; Chastenet et al. 2019). In the Milky Way, the diffuse ISRF varies primarily radially and is dominated by B stars (Mathis et al. 1983). To summarize, we expect there are a number of effects that can alter the PAH emission to gas surface density conversion, not all of which will be mitigated by our local calibration of the scaling in each galaxy.

## 3.2. Polycyclic Aromatic Hydrocarbon Emission to Gas Surface Density Scaling

We adopt PAH emission to CO ratios for each galaxy from Leroy et al. (2023). We note that these values are specific to the four galaxies we target and could change if D/G, $q_{PAH}$, or the diffuse ISRF intensities vary. The average value across all galaxies for the ratio of PAH emission to CO (2–1) integrated intensity are as follows:

$$ \frac{\text{F770W}}{\text{MJy sr}^{-1}} = 0.91 \frac{I_{\text{CO},2-1}}{\text{K km s}^{-1}}, \qquad (1) $$

$$ \frac{\text{F1130W}}{\text{MJy sr}^{-1}} = 1.26 \frac{I_{\text{CO},2-1}}{\text{K km s}^{-1}}. \qquad (2) $$

These ratios vary slightly between galaxies. The ratio for F770W is 0.79, 0.93, 0.72, and 1.29 in IC 5332, NGC 628, 1365, and 7496, respectively. For F1130W, it ranges between 1.07, 1.29, 1.07, and 1.78 for the same galaxies (see Leroy et al. 2023, Table 3, and Table 1 above). Assuming the Milky Way CO-to-H$_2$ conversion factor, including helium and metals, $\alpha_{\text{CO 1-0}}^{\text{MW}} = 4.35 M_{\odot} \text{ pc}^{-2}$ (K km s$^{-1}$; Bolatto et al. 2013), and a CO J = (2–1) to (1–0) line ratio $R_{21} = 0.65$ (den Brok et al. 2021; Leroy et al. 2022), we can convert these ratios into PAH emission per gas surface density, assuming all gas is molecular. These values, for the average ratios shown in Equations (1) and (2), are as follows:

$$ \frac{\Sigma_{\text{gas}}}{M_{\odot} \text{ pc}^2} = 7.4 \times \left( \frac{\text{F770W}}{\text{MJy sr}^{-1}} \right) \left( \frac{\alpha_{\text{CO 1-0}}}{\alpha_{\text{CO 1-0}}^{\text{MW}}} \right) \left( \frac{0.65}{R_{21}} \right), \qquad (3) $$

$$ \frac{\Sigma_{\text{gas}}}{M_{\odot} \text{ pc}^2} = 5.3 \times \left( \frac{\text{F1130W}}{\text{MJy sr}^{-1}} \right) \left( \frac{\alpha_{\text{CO 1-0}}}{\alpha_{\text{CO 1-0}}^{\text{MW}}} \right) \left( \frac{0.65}{R_{21}} \right). \qquad (4) $$

Given that the noise levels are comparable in the F770W and F1130W maps at native resolution (between 0.1 and 0.2 MJy sr$^{-1}$), it is clear that our most sensitive tracer is F1130W, so we proceed in the following to use that as our primary observable. F1130W’s point-spread function (PSF) has a FWHM of 0.36, yielding 15.7, 17.1, 34.2, and 32.6 pc resolution at the distances of IC 5332, NGC 628, NGC 1365, and NGC 7496, respectively. We note that at this resolution, which resolves giant molecular cloud complexes and is significantly smaller than the expected scale height of the gas disk, correcting for the inclination of the galaxy introduces an additional assumption about the extent of the emitting structures along the line of sight. The appropriate inclination correction may differ between clouds or a diffuse medium. Therefore, we proceed in the following without an inclination correction, but note that, if one were applied, the largest effect would be for NGC 1365, which has an inclination of $i = 55^{\circ}$.

For this pilot study, we use a single scaling factor to convert PAH emission to gas surface density. Several of the components in the conversion of PAH emission to gas surface density are likely to change radially in galaxies. The CO-to-H$_2$ conversion factor likely increases with radius, tracking decreasing metallicity (Wolfire et al. 2010; Leroy et al. 2011; Glover & Clark 2012; Bolatto et al. 2013; Madden et al. 2020) and also can decrease in the centers of galaxies (Sandstrom et al. 2013), likely due to excitation and optical depth effects (Teng et al. 2022). Dust-to-gas ratios and PAH fractions may decrease with metallicity, and therefore radius (we note that the observed PAH fraction in our targets does vary with metallicity, but at a relatively minor level given the $\sim Z_{\odot}$ metallicities; Chastenet et al. 2023). The diffuse ISRF also is likely to vary radially, tracing a combination of the stellar mass surface density and star formation rate surface density of the galaxies (Mathis et al. 1983; Wolfire et al. 2003). In future work, we aim to calibrate a radially dependent PAH-to-gas surface density tracer by finding diffuse-ISRF-dominated, CO-bright regions to measure the scaling factor as a function of radius. We note that some of the factors will likely vary with opposing radial trends, potentially leading to a less variable PAH-to-gas conversion than the individual components.

## 4. Results

In Figures 1 to 4 we show the F1130W emission rescaled to gas surface density using the scalings quoted above for each galaxy. In regions without bright CO emission, the diffuse ISM traced by our scaled PAH maps may include a combination of (i) CO-bright H$_2$ missed due to the current CO map sensitivity, (ii) a CO-dark H$_2$ component, or (iii) a mixture of atomic gas in the WNM and CNM phases.

Shielding models for H$_2$ formation (Krumholz et al. 2009; Sternberg et al. 2014; Bialy & Sternberg 2016) predict the transition between H$_2$- or H I-dominated ISM at $\Sigma_{\text{HI}} \approx 10(Z/Z_{\odot})^{-1} M_{\odot} \text{ pc}^{-2}$. In the models, the transition point varies with the assumed geometry and role of self-shielding. In simulations which take into account the filamentary, turbulent nature of the gas, the column density of the H I/H$_2$ transition can vary significantly (e.g., Glover & Mac Low 2007; Glover & Smith 2016). Transitions at $\sim 10 M_{\odot} \text{ pc}^{-2}$ have been measured in the Milky Way (Stanimirović et al. 2014; Lee et al. 2015). This transition is found at higher surface
work to rigorously calibrate PAH emission as a tracer. There is evidence from measurements of depletion that D/Gs vary with ISM density (Jenkins 2009; Roman-Duval et al. 2021) by a factor of up to 2, though this depends on whether the observed depleted oxygen is incorporated into dust (Whittet 2010; Choban et al. 2022). Changes in the D/G between ISM phases at fixed metallicity have also been observed using far-IR, H I, and CO observations (Roman-Duval et al. 2014; Chiang et al. 2018; Vilchez et al. 2019). From high-spatial-resolution studies of Local Group galaxies (Chastenet et al. 2019) and the early PHANGS-JWST targets (Chastenet et al. 2023; Egorov et al. 2023), the PAH fraction appears to be relatively constant in the atomic and molecular phases, but does drop sharply in H II regions, and has a metallicity dependence. Milky Way observations have suggested that PAH abundance can vary between the warm neutral medium (WNM) and cold neutral medium (CNM) (Hensley et al. 2022). The diffuse ISRF in galaxies, judged from the $U_{min}$ component in far-IR SED modeling (e.g., Draine et al. 2007; Aniano et al. 2012), appears to vary smoothly and primarily radially in highly resolved Local Group galaxies (Draine et al. 2014; Chastenet et al. 2019). In the Milky Way, the diffuse ISRF varies primarily radially and is dominated by B stars (Mathis et al. 1983). To summarize, we expect there are a number of effects that can alter the PAH emission to gas surface density conversion, not all of which will be mitigated by our local calibration of the scaling in each galaxy.

## 3.2. Polycyclic Aromatic Hydrocarbon Emission to Gas Surface Density Scaling

We adopt PAH emission to CO ratios for each galaxy from Leroy et al. (2023). We note that these values are specific to the four galaxies we target and could change if D/G, $q_{PAH}$, or the diffuse ISRF intensities vary. The average value across all galaxies for the ratio of PAH emission to CO (2–1) integrated intensity are as follows:

$$ \frac{\text{F770W}}{\text{MJy sr}^{-1}} = 0.91 \frac{I_{\text{CO},2-1}}{\text{K km s}^{-1}}, \qquad (1) $$

$$ \frac{\text{F1130W}}{\text{MJy sr}^{-1}} = 1.26 \frac{I_{\text{CO},2-1}}{\text{K km s}^{-1}}. \qquad (2) $$

These ratios vary slightly between galaxies. The ratio for F770W is 0.79, 0.93, 0.72, and 1.29 in IC 5332, NGC 628, 1365, and 7496, respectively. For F1130W, it ranges between 1.07, 1.29, 1.07, and 1.78 for the same galaxies (see Leroy et al. 2023, Table 3, and Table 1 above). Assuming the Milky Way CO-to-H$_2$ conversion factor, including helium and metals, $\alpha_{\text{CO 1-0}}^{\text{MW}} = 4.35 M_{\odot} \text{ pc}^{-2}$ (K km s$^{-1}$; Bolatto et al. 2013), and a CO J = (2–1) to (1–0) line ratio $R_{21} = 0.65$ (den Brok et al. 2021; Leroy et al. 2022), we can convert these ratios into PAH emission per gas surface density, assuming all gas is molecular. These values, for the average ratios shown in Equations (1) and (2), are as follows:

$$ \frac{\Sigma_{\text{gas}}}{M_{\odot} \text{ pc}^2} = 7.4 \times \left( \frac{\text{F770W}}{\text{MJy sr}^{-1}} \right) \left( \frac{\alpha_{\text{CO 1-0}}}{\alpha_{\text{CO 1-0}}^{\text{MW}}} \right) \left( \frac{0.65}{R_{21}} \right), \qquad (3) $$

$$ \frac{\Sigma_{\text{gas}}}{M_{\odot} \text{ pc}^2} = 5.3 \times \left( \frac{\text{F1130W}}{\text{MJy sr}^{-1}} \right) \left( \frac{\alpha_{\text{CO 1-0}}}{\alpha_{\text{CO 1-0}}^{\text{MW}}} \right) \left( \frac{0.65}{R_{21}} \right). \qquad (4) $$

Given that the noise levels are comparable in the F770W and F1130W maps at native resolution (between 0.1 and 0.2 MJy sr$^{-1}$), it is clear that our most sensitive tracer is F1130W, so we proceed in the following to use that as our primary observable. F1130W’s point-spread function (PSF) has a FWHM of 0.36, yielding 15.7, 17.1, 34.2, and 32.6 pc resolution at the distances of IC 5332, NGC 628, NGC 1365, and NGC 7496, respectively. We note that at this resolution, which resolves giant molecular cloud complexes and is significantly smaller than the expected scale height of the gas disk, correcting for the inclination of the galaxy introduces an additional assumption about the extent of the emitting structures along the line of sight. The appropriate inclination correction may differ between clouds or a diffuse medium. Therefore, we proceed in the following without an inclination correction, but note that, if one were applied, the largest effect would be for NGC 1365, which has an inclination of $i = 55^{\circ}$.

For this pilot study, we use a single scaling factor to convert PAH emission to gas surface density. Several of the components in the conversion of PAH emission to gas surface density are likely to change radially in galaxies. The CO-to-H$_2$ conversion factor likely increases with radius, tracking decreasing metallicity (Wolfire et al. 2010; Leroy et al. 2011; Glover & Clark 2012; Bolatto et al. 2013; Madden et al. 2020) and also can decrease in the centers of galaxies (Sandstrom et al. 2013), likely due to excitation and optical depth effects (Teng et al. 2022). Dust-to-gas ratios and PAH fractions may decrease with metallicity, and therefore radius (we note that the observed PAH fraction in our targets does vary with metallicity, but at a relatively minor level given the $\sim Z_{\odot}$ metallicities; Chastenet et al. 2023). The diffuse ISRF also is likely to vary radially, tracing a combination of the stellar mass surface density and star formation rate surface density of the galaxies (Mathis et al. 1983; Wolfire et al. 2003). In future work, we aim to calibrate a radially dependent PAH-to-gas surface density tracer by finding diffuse-ISRF-dominated, CO-bright regions to measure the scaling factor as a function of radius. We note that some of the factors will likely vary with opposing radial trends, potentially leading to a less variable PAH-to-gas conversion than the individual components.

## 4. Results

In Figures 1 to 4 we show the F1130W emission rescaled to gas surface density using the scalings quoted above for each galaxy. In regions without bright CO emission, the diffuse ISM traced by our scaled PAH maps may include a combination of (i) CO-bright H$_2$ missed due to the current CO map sensitivity, (ii) a CO-dark H$_2$ component, or (iii) a mixture of atomic gas in the WNM and CNM phases.

Shielding models for H$_2$ formation (Krumholz et al. 2009; Sternberg et al. 2014; Bialy & Sternberg 2016) predict the transition between H$_2$- or H I-dominated ISM at $\Sigma_{\text{HI}} \approx 10(Z/Z_{\odot})^{-1} M_{\odot} \text{ pc}^{-2}$. In the models, the transition point varies with the assumed geometry and role of self-shielding. In simulations which take into account the filamentary, turbulent nature of the gas, the column density of the H I/H$_2$ transition can vary significantly (e.g., Glover & Mac Low 2007; Glover & Smith 2016). Transitions at $\sim 10 M_{\odot} \text{ pc}^{-2}$ have been measured in the Milky Way (Stanimirović et al. 2014; Lee et al. 2015). This transition is found at higher surface
<figure><img src="image_4.png" /><figcaption>**Figure 2.** NGC 628 gas surface density from scaled F1130W emission. A single black contour highlights the brightness level where the F1130W emission implies $\Sigma_{\text{gas}} = 7 M_{\odot} \text{ pc}^{-2}$.</figcaption></figure>

density in the Magellanic Clouds (Jameson et al. 2016), following the expected trend with metallicity. Similarly, Schruba et al. (2018) find a $\propto(Z/Z_{\odot})^{-1}$ dependence sampling a range of metallicities in 70 nearby galaxies, primarily on $>100$ pc scales.

The galaxies we target are all at approximately solar metallicity, aside from IC 5332, which has a gradient extending down to $Z \sim 0.5 Z_{\odot}$. We expect therefore that $H\text{I}/H\text{\textsc{2}}$ transitions should occur above the noise floor of our maps ($3\sigma \sim 2 M_{\odot} \text{ pc}^{-2}$), suggesting that much of the faint emission we detect should be arising from atomic gas. We select a threshold $\Sigma_{\text{gas}} = 7 M_{\odot} \text{ pc}^{-2}$ as a conservative choice for where the gas should become atomic-dominated given the predictions from shielding models at the relevant metallicities. The black contour on each image corresponds to this representative value of $\Sigma_{\text{gas}} = 7 M_{\odot} \text{ pc}^{-2}$.

In the four galaxies we see a variety of behavior at low gas surface densities. The maps of NGC 628 and NGC 1365 cover only parts of the galaxy disks, and the majority of pixels in these two maps have surface densities above $7 M_{\odot} \text{ pc}^{-2}$. In NGC 7496 the map extends into the outskirts of the disk, so we can observe the faint, low-surface-density structures in the outskirts. IC 5332 has low surface density throughout its disk, and a large fraction of the emission in the map is below $7 M_{\odot} \text{ pc}^{-2}$. This is in good agreement with the faint CO emission observed in this galaxy (Leroy et al. 2021).

Many of the faint structures in all of the targets appear filamentary. Filaments in the PAH emission are explored in much greater detail in Thilker et al. (2023) and Meidt et al. (2023) in this Issue. In some cases, filamentary structures are observed to have widths comparable to a few resolution elements, and lengths that are substantially longer. For

instance, in a few cases in NGC 628, we measure lengths of up to 400 pc for filamentary structures with $\Sigma_{\text{gas}} < 7 M_{\odot} \text{ pc}^{-2}$, yielding a 10:1 aspect ratio. Many of these filamentary structures appear as spurs or low-surface-density continuations of spiral-arm structures. While such morphology may be expected as a consequence of the spiral structure of the galaxies, previous H I observations have not had the resolution to characterize these outside the Local Group. Throughout the diffuse ISM, structure is evident down to the noise level of the maps ($3\sigma \sim 2 M_{\odot} \text{ pc}^{-2}$).

At present, comparing our predicted surface densities from scaled F1130W to measured CO and H I is limited by the available gas observations. The ALMA CO observations from the PHANGS survey have noise limits near this surface density (Leroy et al. 2021; Leroy et al. 2023). More limiting is the availability and resolution of H I data. NGC 628 was observed as part of the THINGS survey (Walter et al. 2008) and the “robust” weighted map has $\sim 6''$ resolution. NGC 7496 has recent MeerKAT observations at $12''$ resolution. Future observations with higher resolution, higher S/N H I and CO maps for the PHANGS-JWST targets will be necessary to measure gas surface densities for comparison to the F1130W emission at these faint levels.

**5. Measuring Diffuse Interstellar Medium Topology with the Genus Statistic**

One of the distinctive features of the faint regions of the F770W and F1130W maps is the clear presence of holes in the ISM distribution (see Figures 1 to 4). These holes may be due to dynamical effects or feedback, topics explored in more detail in Watkins et al. (2023), Barnes et al. (2023), and Meidt et al. (2023). The pervasive emission in the PAH-tracing bands is one of the reasons these holes appear so clearly in the JWST maps. To quantify the topology of the faint emission, we use the *genus* statistic. The genus statistic is the difference between the number of connected regions above and below a brightness threshold. A positive genus corresponds to a topology reflecting compact regions or sources that are bright compared to the threshold, surrounded by connected networks of lower brightness regions. Positive genus topology is referred to as source-dominated or “meatball”-like. A negative genus is characterized by connected networks of emission above the threshold level and isolated holes below the threshold, which are referred to as a “hole-dominated or a “Swiss cheese” topology (Coles 1988). By computing the genus through a range of intensity values, the resulting genus curve characterizes variations in the topology. When the genus curve crosses zero, the topology changes from “Swiss cheese”-like to “meatball”-like. The genus statistic is a well-known measure of spatial topology that was developed to study the topology of the cosmic microwave background (Coles 1988) and has since been used to quantify the topology of the turbulent ISM, both in simulations and observations (e.g., Kowal & Lazarian 2007; Chepurnov et al. 2008; Burkhart et al. 2010; Koch et al. 2017; Burkhart 2021).

Here we use the implementation of the genus statistic in the TurbuStat<fnref n="46" /> package (Koch et al. 2019). We compute the genus curve for each scaled F1130W image at native resolution using a set of intensities chosen by percentiles in the images. We discard all regions smaller than 20 pixels, 4× the F1130W PSF area, in the images to minimize the effect of noise on the

<footnote n="46">⁴⁶ v1.3; turbustat.readthedocs.io.</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** NGC 628 gas surface density from scaled F1130W emission. A single black contour highlights the brightness level where the F1130W emission implies $\Sigma_{\text{gas}} = 7 M_{\odot} \text{ pc}^{-2}$.</figcaption></figure>

density in the Magellanic Clouds (Jameson et al. 2016), following the expected trend with metallicity. Similarly, Schruba et al. (2018) find a $\propto(Z/Z_{\odot})^{-1}$ dependence sampling a range of metallicities in 70 nearby galaxies, primarily on $>100$ pc scales.

The galaxies we target are all at approximately solar metallicity, aside from IC 5332, which has a gradient extending down to $Z \sim 0.5 Z_{\odot}$. We expect therefore that $H\text{I}/H\text{\textsc{2}}$ transitions should occur above the noise floor of our maps ($3\sigma \sim 2 M_{\odot} \text{ pc}^{-2}$), suggesting that much of the faint emission we detect should be arising from atomic gas. We select a threshold $\Sigma_{\text{gas}} = 7 M_{\odot} \text{ pc}^{-2}$ as a conservative choice for where the gas should become atomic-dominated given the predictions from shielding models at the relevant metallicities. The black contour on each image corresponds to this representative value of $\Sigma_{\text{gas}} = 7 M_{\odot} \text{ pc}^{-2}$.

In the four galaxies we see a variety of behavior at low gas surface densities. The maps of NGC 628 and NGC 1365 cover only parts of the galaxy disks, and the majority of pixels in these two maps have surface densities above $7 M_{\odot} \text{ pc}^{-2}$. In NGC 7496 the map extends into the outskirts of the disk, so we can observe the faint, low-surface-density structures in the outskirts. IC 5332 has low surface density throughout its disk, and a large fraction of the emission in the map is below $7 M_{\odot} \text{ pc}^{-2}$. This is in good agreement with the faint CO emission observed in this galaxy (Leroy et al. 2021).

Many of the faint structures in all of the targets appear filamentary. Filaments in the PAH emission are explored in much greater detail in Thilker et al. (2023) and Meidt et al. (2023) in this Issue. In some cases, filamentary structures are observed to have widths comparable to a few resolution elements, and lengths that are substantially longer. For

instance, in a few cases in NGC 628, we measure lengths of up to 400 pc for filamentary structures with $\Sigma_{\text{gas}} < 7 M_{\odot} \text{ pc}^{-2}$, yielding a 10:1 aspect ratio. Many of these filamentary structures appear as spurs or low-surface-density continuations of spiral-arm structures. While such morphology may be expected as a consequence of the spiral structure of the galaxies, previous H I observations have not had the resolution to characterize these outside the Local Group. Throughout the diffuse ISM, structure is evident down to the noise level of the maps ($3\sigma \sim 2 M_{\odot} \text{ pc}^{-2}$).

At present, comparing our predicted surface densities from scaled F1130W to measured CO and H I is limited by the available gas observations. The ALMA CO observations from the PHANGS survey have noise limits near this surface density (Leroy et al. 2021; Leroy et al. 2023). More limiting is the availability and resolution of H I data. NGC 628 was observed as part of the THINGS survey (Walter et al. 2008) and the “robust” weighted map has $\sim 6''$ resolution. NGC 7496 has recent MeerKAT observations at $12''$ resolution. Future observations with higher resolution, higher S/N H I and CO maps for the PHANGS-JWST targets will be necessary to measure gas surface densities for comparison to the F1130W emission at these faint levels.

## 5. Measuring Diffuse Interstellar Medium Topology with the Genus Statistic

One of the distinctive features of the faint regions of the F770W and F1130W maps is the clear presence of holes in the ISM distribution (see Figures 1 to 4). These holes may be due to dynamical effects or feedback, topics explored in more detail in Watkins et al. (2023), Barnes et al. (2023), and Meidt et al. (2023). The pervasive emission in the PAH-tracing bands is one of the reasons these holes appear so clearly in the JWST maps. To quantify the topology of the faint emission, we use the *genus* statistic. The genus statistic is the difference between the number of connected regions above and below a brightness threshold. A positive genus corresponds to a topology reflecting compact regions or sources that are bright compared to the threshold, surrounded by connected networks of lower brightness regions. Positive genus topology is referred to as source-dominated or “meatball”-like. A negative genus is characterized by connected networks of emission above the threshold level and isolated holes below the threshold, which are referred to as a “hole-dominated or a “Swiss cheese” topology (Coles 1988). By computing the genus through a range of intensity values, the resulting genus curve characterizes variations in the topology. When the genus curve crosses zero, the topology changes from “Swiss cheese”-like to “meatball”-like. The genus statistic is a well-known measure of spatial topology that was developed to study the topology of the cosmic microwave background (Coles 1988) and has since been used to quantify the topology of the turbulent ISM, both in simulations and observations (e.g., Kowal & Lazarian 2007; Chepurnov et al. 2008; Burkhart et al. 2010; Koch et al. 2017; Burkhart 2021).

Here we use the implementation of the genus statistic in the TurbuStat<fnref n="46" /> package (Koch et al. 2019). We compute the genus curve for each scaled F1130W image at native resolution using a set of intensities chosen by percentiles in the images. We discard all regions smaller than 20 pixels, 4× the F1130W PSF area, in the images to minimize the effect of noise on the

<footnote n="46">⁴⁶ v1.3; turbustat.readthedocs.io.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** NGC 1365 gas surface density from scaled F1130W emission. A single black contour highlights the brightness level where the F1130W emission implies $\Sigma_{\text{gas}} = 7 M_{\odot} \text{ pc}^{-2}$. In NGC 1365, saturation in the center leads to an artifact where the columns of pixels with saturation show low values. This is evident in the faint regions of our map.</figcaption></figure>

genus curves, where a connected region includes pixels adjacent along both the edges and corners. We note that the zero-crossing point is not strongly affected by the choice of minimum area (for 10–40 pixels). There is no restriction on the morphology of a region (e.g., round versus elongated).

We additionally impose three masking constraints on the data when computing the genus:

1. We mask the PSF artifacts from the active galactic nuclei in NGC 1365 and NGC 7496. We use PSFs generated with WebbPSF and impose a masking cutoff at 1/1000th the maximum of the normalized PSF.

2. Two negative strips remain in the NGC 1365 (vertical through center; see Figure 3) and NGC 7496 (horizontal through center; see Figure 4) mosaics after our calibration. We impose by-hand masks to remove the visible negative regions.

3. Finally, we apply a radial mask at $R_{\text{gal}} > 0.4 R_{25}$ to provide a consistent basis for comparison between the four targets. This corresponds to radial cutoffs of 3.2, 5.7, 13.7, and 3.6 kpc in IC 5332, NGC 628, NGC 1365, and NGC 7496, respectively. This cutoff was chosen based on the range of $R_{\text{gal}}$ that is well sampled across all four galaxies, and where diffuse PAH emission is visually obvious in Figures 1 to 4.

Masks applied to the data will reduce the value of the genus statistic, roughly proportional to the total area, but will not change the relative shape within a single closed contour.

Because of the wide dynamic range, we empirically choose 600 thresholds linearly sampled from the 0.01th to 98th percentiles, and 50 logarithmically sampled thresholds from the 98th to 99.99th percentiles. This ensures we sample the genus curve well near zero crossings, tracing transitions from “Swiss cheese” to “meatball” topologies in the faint emission tail, while still characterizing the high-end intensity tail.

Figure 5 shows the genus curves in the F1130W band<fnref n="47" /> computed with a cut $R_{\text{gal}} < 0.4 R_{25}$. There is a clear transition in the topology traced by the 11.3 µm emission that varies by galaxy. Table 2 provides the zero crossing points from the genus curves. The genus crossing corresponds to a gas surface density ~8–10 $M_{\odot} \text{ pc}^{-2}$ (scaled from Equation (2)) for NGC 628, NGC 1365, and NGC 7496. This closely matches the expected H$_2$/H I transition point of ~10 $M_{\odot} \text{ pc}^{-2}$ based on shielding models at solar metallicity.

The genus zero crossing point in IC 5332 is substantially lower, at 1.6 $M_{\odot} \text{ pc}^{-2}$ (0.3 MJy sr$^{-1}$), ~2σ significance above the background level of the map, and thus affected by noise fluctuations altering the genus curve. A field of Gaussian noise will produce a similar transition from negative to positive genus, with a zero crossing at the mean field value (Coles 1988). The significant detected structure alters the shape of the genus, as we see in Figure 5, but for the

<footnote n="47">⁴⁷ We note the F770W band has a similar topology and the genus analysis gives similar results, though with lower surface density sensitivity (see Section 2).</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** NGC 1365 gas surface density from scaled F1130W emission. A single black contour highlights the brightness level where the F1130W emission implies $\Sigma_{\text{gas}} = 7 M_{\odot} \text{ pc}^{-2}$. In NGC 1365, saturation in the center leads to an artifact where the columns of pixels with saturation show low values. This is evident in the faint regions of our map.</figcaption></figure>

genus curves, where a connected region includes pixels adjacent along both the edges and corners. We note that the zero-crossing point is not strongly affected by the choice of minimum area (for 10–40 pixels). There is no restriction on the morphology of a region (e.g., round versus elongated).

We additionally impose three masking constraints on the data when computing the genus:

1. We mask the PSF artifacts from the active galactic nuclei in NGC 1365 and NGC 7496. We use PSFs generated with WebbPSF and impose a masking cutoff at 1/1000th the maximum of the normalized PSF.

2. Two negative strips remain in the NGC 1365 (vertical through center; see Figure 3) and NGC 7496 (horizontal through center; see Figure 4) mosaics after our calibration. We impose by-hand masks to remove the visible negative regions.

3. Finally, we apply a radial mask at $R_{\text{gal}} > 0.4 R_{25}$ to provide a consistent basis for comparison between the four targets. This corresponds to radial cutoffs of 3.2, 5.7, 13.7, and 3.6 kpc in IC 5332, NGC 628, NGC 1365, and NGC 7496, respectively. This cutoff was chosen based on the range of $R_{\text{gal}}$ that is well sampled across all four galaxies, and where diffuse PAH emission is visually obvious in Figures 1 to 4.

Masks applied to the data will reduce the value of the genus statistic, roughly proportional to the total area, but will not change the relative shape within a single closed contour.

Because of the wide dynamic range, we empirically choose 600 thresholds linearly sampled from the 0.01th to 98th percentiles, and 50 logarithmically sampled thresholds from the 98th to 99.99th percentiles. This ensures we sample the genus curve well near zero crossings, tracing transitions from “Swiss cheese” to “meatball” topologies in the faint emission tail, while still characterizing the high-end intensity tail.

Figure 5 shows the genus curves in the F1130W band<fnref n="47" /> computed with a cut $R_{\text{gal}} < 0.4 R_{25}$. There is a clear transition in the topology traced by the 11.3 µm emission that varies by galaxy. Table 2 provides the zero crossing points from the genus curves. The genus crossing corresponds to a gas surface density ~8–10 $M_{\odot} \text{ pc}^{-2}$ (scaled from Equation (2)) for NGC 628, NGC 1365, and NGC 7496. This closely matches the expected H$_2$/H I transition point of ~10 $M_{\odot} \text{ pc}^{-2}$ based on shielding models at solar metallicity.

The genus zero crossing point in IC 5332 is substantially lower, at 1.6 $M_{\odot} \text{ pc}^{-2}$ (0.3 MJy sr$^{-1}$), ~2σ significance above the background level of the map, and thus affected by noise fluctuations altering the genus curve. A field of Gaussian noise will produce a similar transition from negative to positive genus, with a zero crossing at the mean field value (Coles 1988). The significant detected structure alters the shape of the genus, as we see in Figure 5, but for the

<footnote n="47">⁴⁷ We note the F770W band has a similar topology and the genus analysis gives similar results, though with lower surface density sensitivity (see Section 2).</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 4.** NGC 7496 gas surface density from scaled F1130W emission. A single black contour highlights the brightness level where the F1130W emission implies $\Sigma_{\text{gas}} = 7 M_{\odot} \text{ pc}^{-2}$. In NGC 7496, saturation in the center leads to an artifact where the columns of pixels with saturation show low values. This is evident in the faint regions of our map.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Genus curves for the F1130W band, scaled to gas surface density (Equation (2)). The gray region shows where the noise dominates the statistic up to $3\sigma$ significance above the background, which corresponds to $\Sigma_{\text{gas}} \approx 2 M_{\odot} \text{ pc}^{-2}$. For a consistent comparison, we impose a radial cut at 0.4 $R_{25}$ in each galaxy, corresponding to 3.2, 5.7, 13.7, and 3.6 kpc in IC 5332, NGC 628, NGC 1365, and NGC 7496, respectively. The 11.3 µm emission traces a clear conversion from a “Swiss cheese” (genus < 0) to “meatball” topology (genus > 0) when going from low to high gas surface density.</figcaption></figure>

<table><caption>**Table 2**
Genus Zero Crossings</caption><thead><tr><th>Target</th><th>Intensity<br>(MJy sr<sup>-1</sup>)</th><th>Σ<sub>gas</sub><br>(M<sub>☉</sub> pc<sup>-2</sup>)</th></tr></thead><tbody><tr><td>IC 5332</td><td>0.3</td><td>1.6</td></tr><tr><td>NGC 628</td><td>2.5</td><td>13.0</td></tr><tr><td>NGC 1365</td><td>1.5</td><td>8.0</td></tr><tr><td>NGC 7496</td><td>1.9</td><td>10.1</td></tr></tbody></table>

**Note.** Zero crossings computed from the genus curves in Figure 5. We scale to a gas surface density using Equation (4).

sensitivity-limited IC 5332, the genus zero crossing is qualitatively consistent with a noise-like field. The diffuse PAH emission in IC 5332 being sensitivity-limited is a substantially lower filling fraction in the PAH emission compared to the other three galaxies (Figure 1). Given that IC 5332 is the only dwarf spiral in this sample, with a lower stellar mass, star formation rate, and metallicity relative to the other galaxies, this lower filling fraction at a common sensitivity limit is expected (e.g., Jameson et al. 2018 and Chevance et al. 2020, in the Magellenic Clouds). The genus zero crossing point for IC 5332 then indicates that the genus is limited to measuring topology changes for tracers with a high filling fraction (see Section 6.2).
<figure><img src="image_6.png" /><figcaption>**Figure 4.** NGC 7496 gas surface density from scaled F1130W emission. A single black contour highlights the brightness level where the F1130W emission implies $\Sigma_{\text{gas}} = 7 M_{\odot} \text{ pc}^{-2}$. In NGC 7496, saturation in the center leads to an artifact where the columns of pixels with saturation show low values. This is evident in the faint regions of our map.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Genus curves for the F1130W band, scaled to gas surface density (Equation (2)). The gray region shows where the noise dominates the statistic up to $3\sigma$ significance above the background, which corresponds to $\Sigma_{\text{gas}} \approx 2 M_{\odot} \text{ pc}^{-2}$. For a consistent comparison, we impose a radial cut at $0.4 R_{25}$ in each galaxy, corresponding to 3.2, 5.7, 13.7, and 3.6 kpc in IC 5332, NGC 628, NGC 1365, and NGC 7496, respectively. The 11.3 µm emission traces a clear conversion from a “Swiss cheese” (genus < 0) to “meatball” topology (genus > 0) when going from low to high gas surface density.</figcaption></figure>

<table><caption>**Table 2**
Genus Zero Crossings</caption><thead><tr><th>Target</th><th>Intensity<br>(MJy sr<sup>-1</sup>)</th><th>Σ<sub>gas</sub><br>(M<sub>☉</sub> pc<sup>-2</sup>)</th></tr></thead><tbody><tr><td>IC 5332</td><td>0.3</td><td>1.6</td></tr><tr><td>NGC 628</td><td>2.5</td><td>13.0</td></tr><tr><td>NGC 1365</td><td>1.5</td><td>8.0</td></tr><tr><td>NGC 7496</td><td>1.9</td><td>10.1</td></tr></tbody></table>

**Note.** Zero crossings computed from the genus curves in Figure 5. We scale to a gas surface density using Equation (4).

sensitivity-limited IC 5332, the genus zero crossing is qualitatively consistent with a noise-like field. The diffuse PAH emission in IC 5332 being sensitivity-limited is a substantially lower filling fraction in the PAH emission compared to the other three galaxies (Figure 1). Given that IC 5332 is the only dwarf spiral in this sample, with a lower stellar mass, star formation rate, and metallicity relative to the other galaxies, this lower filling fraction at a common sensitivity limit is expected (e.g., Jameson et al. 2018 and Chevance et al. 2020, in the Magellenic Clouds). The genus zero crossing point for IC 5332 then indicates that the genus is limited to measuring topology changes for tracers with a high filling fraction (see Section 6.2).
# 6. Discussion

## 6.1. Morphology of the Faint Polycyclic Aromatic Hydrocarbon Emission

In the regions of our F1130W maps corresponding to gas surface densities <7 $M_\odot$ pc$^{-2}$, we see a variety of structures down to our detection limit of ~2 $M_\odot$ pc$^{-2}$, including filaments, clumps, bubbles, interarm spurs, and faint structures in the outskirts of galaxies. The 10–40 pc scales of PAH emission and gas column we explore here are comparable only to previous observations of the Milky Way and Local Group. In the Milky Way, the 8 µm Spitzer/IRAC map from GLIMPSE similarly showed filamentary structures, pervasive PAH emission from diffuse gas, and a porous, bubble-filled structure for the ISM. PAH emission in the Milky Way is observed to be associated with both molecular and atomic gas phases (Churchwell et al. 2009).

The external perspective provided by nearby galaxies at JWST’s resolution may also enable identification of filamentary structures in both molecular- and atomic-dominated gas. For example, at 10–40 pc resolution, our maps reveal filamentary structures with aspect ratios of 10:1 or more and lengths of several 100 pc in both interarm regions (e.g., NGC 628) and in the outskirts of galaxies (NGC 7496, IC 5332; see also Thilker et al. 2023; Meidt et al. 2023). While these filaments are wider (>10–40 pc) and likely associated with diffuse atomic gas, the filament aspect ratios are similar to giant molecular filaments identified in the Milky Way (Ragan et al. 2014; Zucker et al. 2018; Hacar et al. 2022). These may be similar to giant H I filaments (Syed et al. 2022), which have proven difficult to identify even in the Milky Way due to line-of-sight confusion (Soler et al. 2020).

One exciting possibility with gas tracers at 10–40 pc scales extending to H I-dominated gas is structurally separating the CNM and WNM. The CNM is dense atomic gas where $H_2$ formation proceeds, and therefore the balance between the WNM and CNM is central to the regulation of star formation. With existing H I observations, separation of the WNM and CNM outside the Local Group relies on Gaussian decomposition of H I emission spectra, which is uncertain at low spatial resolution (Koch et al. 2021). On ~10 pc scales, the WNM and CNM are expected to have distinct spatial structures: Galactic studies separating the WNM and CNM phases in the Milky Way primarily find a diffuse, widespread WNM component and small-scale CNM filamentary networks (e.g., Clark et al. 2014; Smith et al. 2014; Kalberla & Haud 2018; Marchal et al. 2019; Kalberla et al. 2021). This morphological separation is limited by the H I resolution and sensitivity restricting these studies to the Milky Way, but early GASKAP work points to recovering a pervasive filamentary structure on 10 pc scales in the Magellanic Clouds (Pingel et al. 2022). Given that different phases of the neutral medium have different morphologies, the PAH emission map could allow morphological separation of CNM/WNM components. The emission morphology sets its topology, so the genus statistic and related analyses may offer a new way to characterize ISM phase structure.

## 6.2. The Transition in Interstellar Medium Topology at ~10 $M_\odot$ pc$^{-2}$

A striking visual impression of the JWST F1130W (and other PAH tracing maps) is the clear presence of bubbles throughout the ISM (see Watkins et al. 2023; Barnes et al. 2023,

this Issue). These bubbles are clearly visible because they occur in a pervasive background of PAH emission throughout the galaxies. The genus statistic allows a quantification of this visually evident topology. At faint levels in the maps, the genus statistic is negative, indicating more closed structures below the threshold than above, or a “Swiss cheese”-like topology. From our analysis, we found that the PAH emission transitions from a “Swiss cheese” to a “meatball”-like topology (i.e., sources rather than holes) at ~10 $M_\odot$ pc$^{-2}$ in all galaxies where the topology transition is well measured (NGC 628, 1365, and 7496). Interestingly, this threshold is close to what is expected for the H I/H$_2$ transition.

Our observation of the topology transition point matches well with the expectation that the pervasive, relatively smooth gas distribution in which holes and bubbles are evident is atomic-dominated. Previous observations of nearby galaxies have shown that atomic and molecular gas show distinct spatial distributions. Molecular gas observed in CO emission is found to be “clumpy,” whereas atomic gas shows smoother distributions as traced by 21 cm H I emission on the >100 pc scales that are representative of the WNM (Leroy et al. 2013). However, part of this contrast in structure stems from observational limitations of the 21 cm H I resolution and the CO sensitivity, and so represents the largest possible difference in topologies between the phases. Physically one might also expect this topology transition at ~10 $M_\odot$ pc$^{-2}$ if it is also related to gravitational or thermal instability above this column causing gas to collapse into clumps (Dobbs et al. 2014), which the genus statistic should measure through a transition in topology.

## 6.3. Comparison to Simulations

High-resolution ISM column density maps from PAHs open new opportunities for comparing ISM topology between observations and simulations (Chepurnov et al. 2008; Burkhart et al. 2010; Koch et al. 2017). We briefly explore this application with the genus statistic by comparing between two total neutral gas surface density maps of simulated galaxies, one from an M51-like galaxy (Tress et al. 2020; Treß et al. 2021) and an NGC 628-like galaxy (E. Emsellem et al. 2022, in preparation). Briefly, the NGC 628-like simulation was initialized using known properties of NGC 628 in terms of stellar mass, gas content and distribution and viewing angles, and run with the adaptive-mesh refinement RAMSES code (Teyssier 2002) at a maximum sampling of about 3.7 pc; it includes atomic and molecular cooling, a UV radiation background, star formation with a fixed efficiency per freefall time above a given density threshold, and feedback from Type Ia and Type II SNe, stellar winds and radiation pressure (more details in E. Emsellem et al. 2022, in preparation). The M51-like simulation was initialized to have properties similar to M51. Tress et al. (2020) and Treß et al. (2021) ran simulations both with and without an interacting companion, but here we just consider the simulation without a companion, in which the galaxy develops into a flocculent spiral. The simulation was run with the AREPO moving mesh code (Springel 2010) and achieves a maximum spatial resolution of around 0.5 pc. It includes a similar range of physical processes to the RAMSES calculation, but does not include stellar winds or radiation pressure. Further details can be found in Tress et al. (2020). Figure 6 shows the genus curves and gas surface density maps of the two snapshots, with the contour at the genus zero-crossing point.
# 6. Discussion

## 6.1. Morphology of the Faint Polycyclic Aromatic Hydrocarbon Emission

In the regions of our F1130W maps corresponding to gas surface densities <7 $M_\odot$ pc$^{-2}$, we see a variety of structures down to our detection limit of ~2 $M_\odot$ pc$^{-2}$, including filaments, clumps, bubbles, interarm spurs, and faint structures in the outskirts of galaxies. The 10–40 pc scales of PAH emission and gas column we explore here are comparable only to previous observations of the Milky Way and Local Group. In the Milky Way, the 8 µm Spitzer/IRAC map from GLIMPSE similarly showed filamentary structures, pervasive PAH emission from diffuse gas, and a porous, bubble-filled structure for the ISM. PAH emission in the Milky Way is observed to be associated with both molecular and atomic gas phases (Churchwell et al. 2009).

The external perspective provided by nearby galaxies at JWST’s resolution may also enable identification of filamentary structures in both molecular- and atomic-dominated gas. For example, at 10–40 pc resolution, our maps reveal filamentary structures with aspect ratios of 10:1 or more and lengths of several 100 pc in both interarm regions (e.g., NGC 628) and in the outskirts of galaxies (NGC 7496, IC 5332; see also Thilker et al. 2023; Meidt et al. 2023). While these filaments are wider (>10–40 pc) and likely associated with diffuse atomic gas, the filament aspect ratios are similar to giant molecular filaments identified in the Milky Way (Ragan et al. 2014; Zucker et al. 2018; Hacar et al. 2022). These may be similar to giant H I filaments (Syed et al. 2022), which have proven difficult to identify even in the Milky Way due to line-of-sight confusion (Soler et al. 2020).

One exciting possibility with gas tracers at 10–40 pc scales extending to H I-dominated gas is structurally separating the CNM and WNM. The CNM is dense atomic gas where $H_2$ formation proceeds, and therefore the balance between the WNM and CNM is central to the regulation of star formation. With existing H I observations, separation of the WNM and CNM outside the Local Group relies on Gaussian decomposition of H I emission spectra, which is uncertain at low spatial resolution (Koch et al. 2021). On ~10 pc scales, the WNM and CNM are expected to have distinct spatial structures: Galactic studies separating the WNM and CNM phases in the Milky Way primarily find a diffuse, widespread WNM component and small-scale CNM filamentary networks (e.g., Clark et al. 2014; Smith et al. 2014; Kalberla & Haud 2018; Marchal et al. 2019; Kalberla et al. 2021). This morphological separation is limited by the H I resolution and sensitivity restricting these studies to the Milky Way, but early GASKAP work points to recovering a pervasive filamentary structure on 10 pc scales in the Magellanic Clouds (Pingel et al. 2022). Given that different phases of the neutral medium have different morphologies, the PAH emission map could allow morphological separation of CNM/WNM components. The emission morphology sets its topology, so the genus statistic and related analyses may offer a new way to characterize ISM phase structure.

## 6.2. The Transition in Interstellar Medium Topology at ~10 $M_\odot$ pc$^{-2}$

A striking visual impression of the JWST F1130W (and other PAH tracing maps) is the clear presence of bubbles throughout the ISM (see Watkins et al. 2023; Barnes et al. 2023,

this Issue). These bubbles are clearly visible because they occur in a pervasive background of PAH emission throughout the galaxies. The genus statistic allows a quantification of this visually evident topology. At faint levels in the maps, the genus statistic is negative, indicating more closed structures below the threshold than above, or a “Swiss cheese”-like topology. From our analysis, we found that the PAH emission transitions from a “Swiss cheese” to a “meatball”-like topology (i.e., sources rather than holes) at ~10 $M_\odot$ pc$^{-2}$ in all galaxies where the topology transition is well measured (NGC 628, 1365, and 7496). Interestingly, this threshold is close to what is expected for the H I/H$_2$ transition.

Our observation of the topology transition point matches well with the expectation that the pervasive, relatively smooth gas distribution in which holes and bubbles are evident is atomic-dominated. Previous observations of nearby galaxies have shown that atomic and molecular gas show distinct spatial distributions. Molecular gas observed in CO emission is found to be “clumpy,” whereas atomic gas shows smoother distributions as traced by 21 cm H I emission on the >100 pc scales that are representative of the WNM (Leroy et al. 2013). However, part of this contrast in structure stems from observational limitations of the 21 cm H I resolution and the CO sensitivity, and so represents the largest possible difference in topologies between the phases. Physically one might also expect this topology transition at ~10 $M_\odot$ pc$^{-2}$ if it is also related to gravitational or thermal instability above this column causing gas to collapse into clumps (Dobbs et al. 2014), which the genus statistic should measure through a transition in topology.

## 6.3. Comparison to Simulations

High-resolution ISM column density maps from PAHs open new opportunities for comparing ISM topology between observations and simulations (Chepurnov et al. 2008; Burkhart et al. 2010; Koch et al. 2017). We briefly explore this application with the genus statistic by comparing between two total neutral gas surface density maps of simulated galaxies, one from an M51-like galaxy (Tress et al. 2020; Treß et al. 2021) and an NGC 628-like galaxy (E. Emsellem et al. 2022, in preparation). Briefly, the NGC 628-like simulation was initialized using known properties of NGC 628 in terms of stellar mass, gas content and distribution and viewing angles, and run with the adaptive-mesh refinement RAMSES code (Teyssier 2002) at a maximum sampling of about 3.7 pc; it includes atomic and molecular cooling, a UV radiation background, star formation with a fixed efficiency per freefall time above a given density threshold, and feedback from Type Ia and Type II SNe, stellar winds and radiation pressure (more details in E. Emsellem et al. 2022, in preparation). The M51-like simulation was initialized to have properties similar to M51. Tress et al. (2020) and Treß et al. (2021) ran simulations both with and without an interacting companion, but here we just consider the simulation without a companion, in which the galaxy develops into a flocculent spiral. The simulation was run with the AREPO moving mesh code (Springel 2010) and achieves a maximum spatial resolution of around 0.5 pc. It includes a similar range of physical processes to the RAMSES calculation, but does not include stellar winds or radiation pressure. Further details can be found in Tress et al. (2020). Figure 6 shows the genus curves and gas surface density maps of the two snapshots, with the contour at the genus zero-crossing point.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Genus curves (top panels) from the gas surface density maps for two simulation snapshots (bottom panels). We show a gray vertical line at $\Sigma_{\text{gas}} = 10 M_{\odot} pc^{-2}$ as visual reference point to compare the two genus curves. The red contour indicates the surface density at the topology change measured by the genus zero crossing point. Left: M51-like simulation (Tress et al. 2020; Treß et al. 2021). Right: NGC628-like simulation (E. Emsellem et al. 2022, in preparation).</figcaption></figure>

We recover the two behaviors we observed using the scaled PAH emission in Figure 5: (i) the M51-like simulation has a lower filling fraction of diffuse gas that leads to a lower genus crossing point around $\Sigma_{\text{gas}} \sim 2 M_{\odot} pc^{-2}$, (ii) the NGC628-like simulation has a broad diffuse component and a genus crossing point at $\Sigma_{\text{gas}} \sim 10 M_{\odot} pc^{-2}$, remarkably similar to the transition we find for NGC628, 1365, and 7496. These crossing points remain at similar surface densities after varying the resolution of the two simulated maps to match the resolution of the JWST-MIRI data. However, more work remains to map from variations in the ISM morphology to the underlying physics contained in the simulations. The small transition value for the M51-like simulation differs clearly from our observations and the other simulations. However, this simulation reaches higher resolution than the observations (subparsec in dense regions) and the contour structures appear notably different in Figure 6. This difference suggests that this approach would require testing the effects of changing resolution. We will be able to test this effect and others further by making similar comparisons using a suite of simulations and the complete PHANGS-JWST sample.

Methods like the genus provide a promising route to improve our interpretation of the neutral ISM morphology at high resolution in both observations and simulations: the high-resolution JWST PAH maps provide new constraints for simulations to reproduce the low-surface-density morphology, while close correspondence with topology measures in simulations can help us better interpret how the topology varies across the molecular to atomic ISM phases in observations.

## 7. Conclusions

In this Letter, we present a test case for using PAH emission observed with JWST-MIRI as a high-resolution, high-sensitivity, cold ISM phase-agnostic tracer of gas. We use scalings determined in H₂-dominated (as traced by CO) but diffuse-ISRF-heated gas, and apply those to the faint PAH emission in the F1130W map. With these simplistic assumptions, we trace gas surface densities down to $\Sigma_{\text{gas}} \sim 2 M_{\odot} pc^{-2}$ at $3\sigma$ at 10–40 pc resolution. This approach is enabled by the high angular resolution of JWST, which allows ISM gas outside the sphere of influence of star-forming regions, heated by the diffuse ISRF, to be isolated to calibrate our gas tracer.

Using these gas maps, we isolate regions at $<7 M_{\odot} pc^{-2}$ which we expect to be atomic gas dominated. The existing H I and CO maps are not sensitive or high enough resolution to determine definitively if the faint PAH emission arises from atomic gas, CO-bright H₂ below the detection threshold, or CO-dark H₂. Given the predictions from H₂ formation models and the ~solar metallicity of our targets, we expect that gas at $<7 M_{\odot} pc^{-2}$ is primarily atomic.

At these low surface densities, we see a wealth of structures, from filaments, clumps, interarm emission, bubbles, and faint galaxy outskirts. These structures are inaccessible at this resolution in our targets with current H I and CO observations. We highlight the possibility of using PAH emission to separate the CNM and WNM gas by its spatial structure, which would make the WNM/CNM balance an accessible observable in galaxies outside the Local Group. We also note that the aspect ratios and lengths of filamentary structures evident at low gas surface densities in our maps suggest we may be seeing the equivalent of giant H I filaments in other galaxies.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Genus curves (top panels) from the gas surface density maps for two simulation snapshots (bottom panels). We show a gray vertical line at $\Sigma_{\text{gas}} = 10 M_{\odot} pc^{-2}$ as visual reference point to compare the two genus curves. The red contour indicates the surface density at the topology change measured by the genus zero crossing point. Left: M51-like simulation (Tress et al. 2020; Treß et al. 2021). Right: NGC628-like simulation (E. Emsellem et al. 2022, in preparation).</figcaption></figure>

We recover the two behaviors we observed using the scaled PAH emission in Figure 5: (i) the M51-like simulation has a lower filling fraction of diffuse gas that leads to a lower genus crossing point around $\Sigma_{\text{gas}} \sim 2 M_{\odot} pc^{-2}$, (ii) the NGC628-like simulation has a broad diffuse component and a genus crossing point at $\Sigma_{\text{gas}} \sim 10 M_{\odot} pc^{-2}$, remarkably similar to the transition we find for NGC628, 1365, and 7496. These crossing points remain at similar surface densities after varying the resolution of the two simulated maps to match the resolution of the JWST-MIRI data. However, more work remains to map from variations in the ISM morphology to the underlying physics contained in the simulations. The small transition value for the M51-like simulation differs clearly from our observations and the other simulations. However, this simulation reaches higher resolution than the observations (subparsec in dense regions) and the contour structures appear notably different in Figure 6. This difference suggests that this approach would require testing the effects of changing resolution. We will be able to test this effect and others further by making similar comparisons using a suite of simulations and the complete PHANGS-JWST sample.

Methods like the genus provide a promising route to improve our interpretation of the neutral ISM morphology at high resolution in both observations and simulations: the high-resolution JWST PAH maps provide new constraints for simulations to reproduce the low-surface-density morphology, while close correspondence with topology measures in simulations can help us better interpret how the topology varies across the molecular to atomic ISM phases in observations.

## 7. Conclusions

In this Letter, we present a test case for using PAH emission observed with JWST-MIRI as a high-resolution, high-sensitivity, cold ISM phase-agnostic tracer of gas. We use scalings determined in H₂-dominated (as traced by CO) but diffuse-ISRF-heated gas, and apply those to the faint PAH emission in the F1130W map. With these simplistic assumptions, we trace gas surface densities down to $\Sigma_{\text{gas}} \sim 2 M_{\odot} pc^{-2}$ at $3\sigma$ at 10–40 pc resolution. This approach is enabled by the high angular resolution of JWST, which allows ISM gas outside the sphere of influence of star-forming regions, heated by the diffuse ISRF, to be isolated to calibrate our gas tracer.

Using these gas maps, we isolate regions at $<7 M_{\odot} pc^{-2}$ which we expect to be atomic gas dominated. The existing H I and CO maps are not sensitive or high enough resolution to determine definitively if the faint PAH emission arises from atomic gas, CO-bright H₂ below the detection threshold, or CO-dark H₂. Given the predictions from H₂ formation models and the ~solar metallicity of our targets, we expect that gas at $<7 M_{\odot} pc^{-2}$ is primarily atomic.

At these low surface densities, we see a wealth of structures, from filaments, clumps, interarm emission, bubbles, and faint galaxy outskirts. These structures are inaccessible at this resolution in our targets with current H I and CO observations. We highlight the possibility of using PAH emission to separate the CNM and WNM gas by its spatial structure, which would make the WNM/CNM balance an accessible observable in galaxies outside the Local Group. We also note that the aspect ratios and lengths of filamentary structures evident at low gas surface densities in our maps suggest we may be seeing the equivalent of giant H I filaments in other galaxies.
We also measure the topology of the PAH emission using the genus statistic. We find that the PAH emission at faint levels shows a “Swiss cheese”-like morphology, suggesting relatively smooth, diffuse emission with holes and bubbles. At around $10 M_\odot\text{pc}^{-2}$ we find a transition from “Swiss cheese”-like to “meatball”-like or source-dominated topology. The fact that this transition occurs at a similar surface density to the $H_1/H_2$ transition agrees with previous observations that $H_2$ tends to have a clumpy distribution while $H_1$ is smoother.

This initial investigation shows the potential for using JWST observations of PAH emission to reveal key aspects of the structure and phase of ISM gas in nearby galaxies. Future work that calibrates these relationships with higher resolution and higher sensitivity $H_1$ and CO will be necessary. In addition, deeper understanding of the variation of D/Gs, PAH fractions, and ISRF intensity and spectra will be crucial to fully making use of this new tool.

We thank the anonymous referee for helpful feedback that improved the paper. This work is based on observations made with the NASA/ESA/CSA James Webb Space Telescope. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with program 2107. The specific observations analyzed can be accessed via 10.17909/9bdf-jn24

K.S. acknowledges funding support from grant support by JWST-GO-02107.006-A. E.W.K. acknowledges support from the Smithsonian Institution as a Submillimeter Array (SMA) Fellow and the Natural Sciences and Engineering Research Council of Canada. H.A.P. acknowledges support by the National Science and Technology Council of Taiwan under grant No. 110-2112-M-032-020-MY3. J.M.D.K. gratefully acknowledges funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program via the ERC Starting Grant MUSTANG (grant agreement No. 714907). COOL Research DAO is a Decentralized Autonomous Organization supporting research in astrophysics aimed at uncovering our cosmic origins. M.B. acknowledges support from FONDECYT regular grant No. 1211000 and by the ANID BASAL project FB210003. I.C. thanks the National Science and Technology Council for support through grant Nos. 108-2112-M-001-007-MY3 and 111-2112-M-001-038-MY3, and the Academia Sinica for Investigator Award AS-IA-109-M02. E.J.W. acknowledges the funding provided by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—Project-ID 138713538—SFB 881 (“The Milky Way System,” subproject P1). K.K. gratefully acknowledges funding from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) in the form of an Emmy Noether Research Group (grant number KR4598/2-1; PI: Kreckel). M.C. gratefully acknowledges funding from the DFG through an Emmy Noether Research Group (grant No. CH2137/1-1). T.G.W. and E.S. acknowledge funding from the European Research Council (ERC) under the European Unions Horizon 2020 research and innovation program (grant agreement No. 694343). F.B. would like to acknowledge funding from the European Research Council (ERC) under the European Unions Horizon 2020 research and innovation program (grant agreement No. 726384/Empire).

R.S.K. acknowledges financial support from the European Research Council via the ERC Synergy Grant “ECOGAL” (project ID 855130), from the Deutsche Forschungs-gemeinschaft (DFG) via the Collaborative Research Center “The Milky Way System” (SFB 881—funding ID 138713538—subprojects A1, B1, B2 and B8) and from the Heidelberg Cluster of Excellence (EXC 2181-390900948) “STRUC-TURES,” funded by the German Excellence Strategy. R.S.K. also thanks the German Ministry for Economic Affairs and Climate Action for funding in project “MAINN” (funding ID 50002206). M.Q. acknowledges support from the Spanish grant PID2019-106027GA-C44, funded by MCIN/AEI/10.13039/501100011033. E.R. and H.H. acknowledge the support of the Natural Sciences and Engineering Research Council of Canada (NSERC), funding reference No. RGPIN-2022-03499. This research was supported by the Excellence Cluster ORIGINS which is funded by the Deutsche Forschungs-gemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy—EXC-2094-390783311.

Some of the simulations in this paper have been carried out on the computing facilities of the Computational Center for Particle and Astrophysics (C2PAP). E.E. would like to thank Alexey Krukau and Margarita Petkova for their support through C2PAP. K.G. is supported by the Australian Research Council through the Discovery Early Career Researcher Award (DECRA) Fellowship DE220100766 funded by the Australian Govern-ment. K.G. is supported by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. J.K. gratefully acknowledges funding from the Deutsche Forschungs-gemeinschaft (DFG, German Research Foundation) through the DFG Sachbeihilfe (grant No. KR4801/2-1). A.S. is supported by an NSF Astronomy and Astrophysics Postdoctoral Fellow-ship under award AST-1903834. A.K.L. gratefully acknowl-edges support by grant Nos. 1653300 and 2205628 from the National Science Foundation, by award JWST-GO-02107.009-A, and by a Humboldt Research Award from the Alexander von Humboldt Foundation. G.A.B. acknowledges the support from ANID Basal project FB210003.

**Facilities:** JWST, ALMA.

**Software:** astropy (Astropy Collaboration et al. 2013,2018, 2022), TurbuStat (Koch et al. 2019), WebbPSF (Oschmann et al. 2014).

## ORCID iDs

Karin M. Sandstrom https://orcid.org/0000-0002-4378-8534

Eric W. Koch https://orcid.org/0000-0001-9605-780X

Adam K. Leroy https://orcid.org/0000-0002-2545-1700

Erik Rosolowsky https://orcid.org/0000-0002-5204-2259

Eric Emsellem https://orcid.org/0000-0002-6155-7166

Rowan J. Smith https://orcid.org/0000-0002-0820-1814

Oleg V. Egorov https://orcid.org/0000-0002-4755-118X

Thomas G. Williams https://orcid.org/0000-0002-0012-2142

Kirsten L. Larson https://orcid.org/0000-0003-3917-6460

Janice C. Lee https://orcid.org/0000-0002-2278-9407

Eva Schinnerer https://orcid.org/0000-0002-3933-7677

David A. Thilker https://orcid.org/0000-0002-8528-7340

Ashley T. Barnes https://orcid.org/0000-0003-0410-4504

Francesco Belfiore https://orcid.org/0000-0002-2545-5752

F. Bigiel https://orcid.org/0000-0003-0166-9745
We also measure the topology of the PAH emission using the genus statistic. We find that the PAH emission at faint levels shows a “Swiss cheese”-like morphology, suggesting relatively smooth, diffuse emission with holes and bubbles. At around $10 M_\odot\text{pc}^{-2}$ we find a transition from “Swiss cheese”-like to “meatball”-like or source-dominated topology. The fact that this transition occurs at a similar surface density to the $H_1/H_2$ transition agrees with previous observations that $H_2$ tends to have a clumpy distribution while $H_1$ is smoother.

This initial investigation shows the potential for using JWST observations of PAH emission to reveal key aspects of the structure and phase of ISM gas in nearby galaxies. Future work that calibrates these relationships with higher resolution and higher sensitivity $H_1$ and CO will be necessary. In addition, deeper understanding of the variation of D/Gs, PAH fractions, and ISRF intensity and spectra will be crucial to fully making use of this new tool.

We thank the anonymous referee for helpful feedback that improved the paper. This work is based on observations made with the NASA/ESA/CSA James Webb Space Telescope. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with program 2107. The specific observations analyzed can be accessed via 10.17909/9bdf-jn24

K.S. acknowledges funding support from grant support by JWST-GO-02107.006-A. E.W.K. acknowledges support from the Smithsonian Institution as a Submillimeter Array (SMA) Fellow and the Natural Sciences and Engineering Research Council of Canada. H.A.P. acknowledges support by the National Science and Technology Council of Taiwan under grant No. 110-2112-M-032-020-MY3. J.M.D.K. gratefully acknowledges funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program via the ERC Starting Grant MUSTANG (grant agreement No. 714907). COOL Research DAO is a Decentralized Autonomous Organization supporting research in astrophysics aimed at uncovering our cosmic origins. M.B. acknowledges support from FONDECYT regular grant No. 1211000 and by the ANID BASAL project FB210003. I.C. thanks the National Science and Technology Council for support through grant Nos. 108-2112-M-001-007-MY3 and 111-2112-M-001-038-MY3, and the Academia Sinica for Investigator Award AS-IA-109-M02. E.J.W. acknowledges the funding provided by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—Project-ID 138713538—SFB 881 (“The Milky Way System,” subproject P1). K.K. gratefully acknowledges funding from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) in the form of an Emmy Noether Research Group (grant number KR4598/2-1; PI: Kreckel). M.C. gratefully acknowledges funding from the DFG through an Emmy Noether Research Group (grant No. CH2137/1-1). T.G.W. and E.S. acknowledge funding from the European Research Council (ERC) under the European Unions Horizon 2020 research and innovation program (grant agreement No. 694343). F.B. would like to acknowledge funding from the European Research Council (ERC) under the European Unions Horizon 2020 research and innovation program (grant agreement No. 726384/Empire).

R.S.K. acknowledges financial support from the European Research Council via the ERC Synergy Grant “ECOGAL” (project ID 855130), from the Deutsche Forschungs-gemeinschaft (DFG) via the Collaborative Research Center “The Milky Way System” (SFB 881—funding ID 138713538—subprojects A1, B1, B2 and B8) and from the Heidelberg Cluster of Excellence (EXC 2181-390900948) “STRUC-TURES,” funded by the German Excellence Strategy. R.S.K. also thanks the German Ministry for Economic Affairs and Climate Action for funding in project “MAINN” (funding ID 50002206). M.Q. acknowledges support from the Spanish grant PID2019-106027GA-C44, funded by MCIN/AEI/10.13039/501100011033. E.R. and H.H. acknowledge the support of the Natural Sciences and Engineering Research Council of Canada (NSERC), funding reference No. RGPIN-2022-03499. This research was supported by the Excellence Cluster ORIGINS which is funded by the Deutsche Forschungs-gemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy—EXC-2094-390783311.

Some of the simulations in this paper have been carried out on the computing facilities of the Computational Center for Particle and Astrophysics (C2PAP). E.E. would like to thank Alexey Krukau and Margarita Petkova for their support through C2PAP. K.G. is supported by the Australian Research Council through the Discovery Early Career Researcher Award (DECRA) Fellowship DE220100766 funded by the Australian Government. K.G. is supported by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. J.K. gratefully acknowledges funding from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through the DFG Sachbeihilfe (grant No. KR4801/2-1). A.S. is supported by an NSF Astronomy and Astrophysics Postdoctoral Fellowship under award AST-1903834. A.K.L. gratefully acknowledges support by grant Nos. 1653300 and 2205628 from the National Science Foundation, by award JWST-GO-02107.009-A, and by a Humboldt Research Award from the Alexander von Humboldt Foundation. G.A.B. acknowledges the support from ANID Basal project FB210003.

**Facilities:** JWST, ALMA.

**Software:** astropy (Astropy Collaboration et al. 2013,2018, 2022), TurbuStat (Koch et al. 2019), WebbPSF (Oschmann et al. 2014).

## ORCID iDs

Karin M. Sandstrom https://orcid.org/0000-0002-4378-8534

Eric W. Koch https://orcid.org/0000-0001-9605-780X

Adam K. Leroy https://orcid.org/0000-0002-2545-1700

Erik Rosolowsky https://orcid.org/0000-0002-5204-2259

Eric Emsellem https://orcid.org/0000-0002-6155-7166

Rowan J. Smith https://orcid.org/0000-0002-0820-1814

Oleg V. Egorov https://orcid.org/0000-0002-4755-118X

Thomas G. Williams https://orcid.org/0000-0002-0012-2142

Kirsten L. Larson https://orcid.org/0000-0003-3917-6460

Janice C. Lee https://orcid.org/0000-0002-2278-9407

Eva Schinnerer https://orcid.org/0000-0002-3933-7677

David A. Thilker https://orcid.org/0000-0002-8528-7340

Ashley T. Barnes https://orcid.org/0000-0003-0410-4504

Francesco Belfiore https://orcid.org/0000-0002-2545-5752

F. Bigiel https://orcid.org/0000-0003-0166-9745
Guillermo A. Blanc https://orcid.org/0000-0003-4218-3944

Alberto D. Bolatto https://orcid.org/0000-0002-5480-5686

Médéric Boquien https://orcid.org/0000-0003-0946-6176

Yixian Cao https://orcid.org/0000-0001-5301-1326

Jérémy Chastenet https://orcid.org/0000-0002-5235-5589

Mélanie Chevance https://orcid.org/0000-0002-5635-5180

I-Da Chiang (江宜達) https://orcid.org/0000-0003-2551-7148

Daniel A. Dale https://orcid.org/0000-0002-5782-9093

Christopher M. Faesi https://orcid.org/0000-0001-5310-467X

Simon C. O. Glover https://orcid.org/0000-0001-6708-1317

Kathryn Grasha https://orcid.org/0000-0002-3247-5321

Brent Groves https://orcid.org/0000-0002-9768-0246

Hamid Hassani https://orcid.org/0000-0002-8806-6308

Jonathan D. Henshaw https://orcid.org/0000-0001-9656-7682

Annie Hughes https://orcid.org/0000-0002-9181-1161

Jaeyeon Kim https://orcid.org/0000-0002-0432-6847

Ralf S. Klessen https://orcid.org/0000-0002-0560-3172

Kathryn Kreckel https://orcid.org/0000-0001-6551-3091

J. M. Diederik Kruijssen https://orcid.org/0000-0002-8804-0212

Laura A. Lopez https://orcid.org/0000-0002-1790-3148

Daizhong Liu https://orcid.org/0000-0001-9773-7479

Sharon E. Meidt https://orcid.org/0000-0002-6118-4048

Eric J. Murphy https://orcid.org/0000-0001-7089-7325

Hsi-An Pan https://orcid.org/0000-0002-1370-6964

Miguel Querejeta https://orcid.org/0000-0002-0472-1011

Toshiki Saito https://orcid.org/0000-0001-9016-2641

Amy Sardone https://orcid.org/0000-0002-5783-145X

Mattia C. Sormani https://orcid.org/0000-0001-6113-6241

Jessica Sutter https://orcid.org/0000-0002-9183-8102

Antonio Usero https://orcid.org/0000-0003-1242-505X

Elizabeth J. Watkins https://orcid.org/0000-0002-7365-5791

## References

Anand, G. S., Lee, J. C., Van Dyk, S. D., et al. 2021, MNRAS, 501, 3621

Aniano, G., Draine, B. T., Calzetti, D., et al. 2012, ApJ, 756, 138

Aniano, G., Draine, B. T., Hunt, L. K., et al. 2020, ApJ, 889, 150

Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167

Astropy Collaboration, Price-Whelan, A. M., Sipocz, B. M., et al. 2018, AJ, 156, 123

Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33

Barnes, A., Watkins, E., Meidt, S. E., et al. 2023, ApJL, 944, L22

Bialy, S., & Sternberg, A. 2016, ApJ, 822, 83

Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207

Bot, C., Boulanger, F., Rubio, M., & Rantakyro, F. 2007, A&A, 471, 103

Boulanger, F., Baud, B., & van Albada, G. D. 1985, A&A, 144, L9

Boulanger, F., & Perault, M. 1988, ApJ, 330, 964

Burkhart, B. 2021, PASP, 133, 102001

Burkhart, B., Stanimirovic, S., Lazarian, A., & Kowal, G. 2010, ApJ, 708, 1204

Chastenet, J., Sandstrom, K., Chiang, I.-D., et al. 2019, ApJ, 876, 62

Chastenet, J., Sutter, J., Sandstrom, K., et al. 2023, ApJL, 944, L11

Chepurnov, A., Gordon, J., Lazarian, A., & Stanimirovic, S. 2008, ApJ, 688, 1021

Chevance, M., Madden, S. C., Fischer, C., et al. 2020, MNRAS, 494, 5279

Chiang, I.-D., Sandstrom, K. M., Chastenet, J., et al. 2018, ApJ, 865, 117

Choban, C. R., Kereš, D., Hopkins, P. F., et al. 2022, MNRAS, 514, 4506

Churchwell, E., Babler, B. L., Meade, M. R., et al. 2009, PASP, 121, 213

Clark, S. E., Peek, J. E. G., & Putman, M. E. 2014, ApJ, 789, 82

Coles, P. 1988, MNRAS, 234, 509

Compiègne, M., Flagey, N., Noriega-Crespo, A., et al. 2010, ApJL, 724, L44

den Brok, J. S., Chatzigiannakis, D., Bigiel, F., et al. 2021, MNRAS, 504, 3221

Dobbs, C. L., Krumholz, M. R., Ballesteros-Paredes, J., et al. 2014, in Protostars and Planets VI, ed. H. Beuther et al. (Tucson, AZ: Univ. of Arizona Press), 3

Draine, B. T., Aniano, G., Krause, O., et al. 2014, ApJ, 780, 172

Draine, B. T., Dale, D. A., Bendo, G., et al. 2007, ApJ, 663, 866

Draine, B. T., & Li, A. 2007, ApJ, 657, 810

Dwek, E., Arendt, R. G., Fixsen, D. J., et al. 1997, ApJ, 475, 565

Eales, S., Smith, M. W. L., Auld, R., et al. 2012, ApJ, 761, 168

Ehlerová, S., & Palouš, J. 2005, A&A, 437, 101

Egorov, Kreckel, K., Sandstrom, K. M., et al. 2023, ApJL, 944, L16

Galliano, F., Galametz, M., & Jones, A. P. 2018, ARA&A, 56, 673

Giard, M., Lamarre, J. M., Pajot, F., & Serra, G. 1994, A&A, 286, 203

Giard, M., Pajot, F., Lamarre, J. M., et al. 1988, A&A, 201, L1

Glover, S. C. O., & Clark, P. C. 2012, MNRAS, 426, 377

Glover, S. C. O., & Mac Low, M.-M. 2007, ApJ, 659, 1317

Glover, S. C. O., & Smith, R. J. 2016, MNRAS, 462, 3011

Groves, B. A., Schinnerer, E., Leroy, A., et al. 2015, ApJ, 799, 96

Groves, B., Kreckel, K., Santoro, F., et al. 2023, arXiv:2301.03811

Hacar, A., Clark, S., Heitsch, F., et al. 2022, arXiv:2203.09562

Hassani, H., Rosolowsky, E. W., Leroy, A. K., et al. 2023, ApJL, 944, L21

Hensley, B. S., Murray, C. E., & Dodici, M. 2022, ApJ, 929, 23

Hoyer, N., Pinna, F., Kamlah, A. W. H., et al. 2023, ApJL, 944, L25

Jameson, K. E., Bolatto, A. D., Leroy, A. K., et al. 2016, ApJ, 825, 12

Jameson, K. E., Bolatto, A. D., Wolfire, M., et al. 2018, ApJ, 853, 111

Jenkins, E. B. 2009, ApJ, 700, 1299

Kalberla, P. M. W., & Haud, U. 2018, A&A, 619, A58

Kalberla, P. M. W., & Kerp, J. 2009, ARA&A, 47, 27

Kalberla, P. M. W., Kerp, J., & Haud, U. 2021, A&A, 654, A91

Koch, E. W., Rosolowsky, E. W., Boyden, R. D., et al. 2019, AJ, 158, 1

Koch, E. W., Rosolowsky, E. W., Leroy, A. K., et al. 2021, MNRAS, 504, 1801

Koch, E. W., Ward, C. G., Offner, S., Loeppky, J. L., & Rosolowsky, E. W. 2017, MNRAS, 471, 1506

Kowal, G., & Lazarian, A. 2007, ApJL, 666, L69

Krumholz, M. R., McKee, C. F., & Tumlinson, J. 2009, ApJ, 693, 216

Lang, P., Meidt, S. E., Rosolowsky, E., et al. 2020, ApJ, 897, 122

Lawton, B., Gordon, K. D., Babler, B., et al. 2010, ApJ, 716, 453

Lee, J. C., Sandstrom, K. M., Leroy, A. K., et al. 2023, ApJL, 944, L17

Lee, M.-Y., Stanimirović, S., Murray, C. E., Heiles, C., & Miller, J. 2015, ApJ, 809, 56

Leroy, A. K., Bolatto, A., Gordon, K., et al. 2011, ApJ, 737, 12

Leroy, A. K., Lee, C., Schruba, A., et al. 2013, ApJL, 769, L12

Leroy, A. K., Rosolowsky, E., Usero, A., et al. 2022, ApJ, 927, 149

Leroy, A. K., Sandstrom, K., Rosolowsky, E., et al. 2023, ApJL,

Leroy, A. K., Schinnerer, E., Hughes, A., et al. 2021, ApJS, 257, 43

Li, A., & Draine, B. T. 2001, ApJ, 554, 778

Li, C., Qiu, K., Hu, B., & Cao, Y. 2021, ApJL, 918, L2

Madden, S. C., Cormier, D., Hony, S., et al. 2020, A&A, 643, A141

Marchal, A., Miville-Deschênes, M.-A., Orieux, F., et al. 2019, A&A, 626, A101

Mathis, J. S., Mezger, P. G., & Panagia, N. 1983, A&A, 128, 212

Meidt, S. E., Rosolowsky, E., Sun, J., et al. 2023, ApJL, 944, L18

Peretto, N., & Fuller, G. A. 2009, A&A, 505, 405

Oschmann, J. M., Clampin, M., Fazio, G. G., et al. 2014, Proc. SPIE, 9143, 91433X

Pilyugin, L. S., & Grebel, E. K. 2016, MNRAS, 457, 3678

Pingel, N. M., Dempsey, J., McClure-Griffiths, N. M., et al. 2022, PASA, 39, e005

Puget, J. L., Leger, A., & Boulanger, F. 1985, A&A, 142, L19

Ragan, S. E., Henning, T., Tackenberg, J., et al. 2014, A&A, 568, A73

Roman-Duval, J., Gordon, K. D., Meixner, M., et al. 2014, ApJ, 797, 86

Roman-Duval, J., Jenkins, E. B., Tchernyshyov, K., et al. 2021, ApJ, 910, 95

Sandstrom, K. M., Leroy, A. K., Walter, F., et al. 2013, ApJ, 777, 5

Schruba, A., Bialy, S., & Sternberg, A. 2018, ApJ, 862, 110

Scoville, N., Aussel, H., Sheth, K., et al. 2014, ApJ, 783, 84

Smith, R. J., Glover, S. C. O., & Klessen, R. 2014, MNRAS, 445, 2900

Soler, J. D., Beuther, H., Syed, J., et al. 2020, A&A, 642, A163

Soler, J. D., Miville-Deschênes, M. A., Molinari, S., et al. 2022, A&A, 662, A96

Springel, V. 2010, MNRAS, 401, 791

Stanimirović, S., Murray, C. E., Lee, M.-Y., Heiles, C., & Miller, J. 2014, ApJ, 793, 132

Sternberg, A., Le Petit, F., Roueff, E., & Le Bourlot, J. 2014, ApJ, 790, 10

Suad, L. A., Caiafa, C. F., Arnal, E. M., & Cichowolski, S. 2014, A&A, 564, A116
Guillermo A. Blanc https://orcid.org/0000-0003-4218-3944

Alberto D. Bolatto https://orcid.org/0000-0002-5480-5686

Médéric Boquien https://orcid.org/0000-0003-0946-6176

Yixian Cao https://orcid.org/0000-0001-5301-1326

Jérémy Chastenet https://orcid.org/0000-0002-5235-5589

Mélanie Chevance https://orcid.org/0000-0002-5635-5180

I-Da Chiang (江宜達) https://orcid.org/0000-0003-2551-7148

Daniel A. Dale https://orcid.org/0000-0002-5782-9093

Christopher M. Faesi https://orcid.org/0000-0001-5310-467X

Simon C. O. Glover https://orcid.org/0000-0001-6708-1317

Kathryn Grasha https://orcid.org/0000-0002-3247-5321

Brent Groves https://orcid.org/0000-0002-9768-0246

Hamid Hassani https://orcid.org/0000-0002-8806-6308

Jonathan D. Henshaw https://orcid.org/0000-0001-9656-7682

Annie Hughes https://orcid.org/0000-0002-9181-1161

Jaeyeon Kim https://orcid.org/0000-0002-0432-6847

Ralf S. Klessen https://orcid.org/0000-0002-0560-3172

Kathryn Kreckel https://orcid.org/0000-0001-6551-3091

J. M. Diederik Kruijssen https://orcid.org/0000-0002-8804-0212

Laura A. Lopez https://orcid.org/0000-0002-1790-3148

Daizhong Liu https://orcid.org/0000-0001-9773-7479

Sharon E. Meidt https://orcid.org/0000-0002-6118-4048

Eric J. Murphy https://orcid.org/0000-0001-7089-7325

Hsi-An Pan https://orcid.org/0000-0002-1370-6964

Miguel Querejeta https://orcid.org/0000-0002-0472-1011

Toshiki Saito https://orcid.org/0000-0001-9016-2641

Amy Sardone https://orcid.org/0000-0002-5783-145X

Mattia C. Sormani https://orcid.org/0000-0001-6113-6241

Jessica Sutter https://orcid.org/0000-0002-9183-8102

Antonio Usero https://orcid.org/0000-0003-1242-505X

Elizabeth J. Watkins https://orcid.org/0000-0002-7365-5791

## References

Anand, G. S., Lee, J. C., Van Dyk, S. D., et al. 2021, MNRAS, 501, 3621

Aniano, G., Draine, B. T., Calzetti, D., et al. 2012, ApJ, 756, 138

Aniano, G., Draine, B. T., Hunt, L. K., et al. 2020, ApJ, 889, 150

Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167

Astropy Collaboration, Price-Whelan, A. M., Sipocz, B. M., et al. 2018, AJ, 156, 123

Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33

Barnes, A., Watkins, E., Meidt, S. E., et al. 2023, ApJL, 944, L22

Bialy, S., & Sternberg, A. 2016, ApJ, 822, 83

Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207

Bot, C., Boulanger, F., Rubio, M., & Rantakyro, F. 2007, A&A, 471, 103

Boulanger, F., Baud, B., & van Albada, G. D. 1985, A&A, 144, L9

Boulanger, F., & Perault, M. 1988, ApJ, 330, 964

Burkhart, B. 2021, PASP, 133, 102001

Burkhart, B., Stanimirovic, S., Lazarian, A., & Kowal, G. 2010, ApJ, 708, 1204

Chastenet, J., Sandstrom, K., Chiang, I.-D., et al. 2019, ApJ, 876, 62

Chastenet, J., Sutter, J., Sandstrom, K., et al. 2023, ApJL, 944, L11

Chepurnov, A., Gordon, J., Lazarian, A., & Stanimirovic, S. 2008, ApJ, 688, 1021

Chevance, M., Madden, S. C., Fischer, C., et al. 2020, MNRAS, 494, 5279

Chiang, I.-D., Sandstrom, K. M., Chastenet, J., et al. 2018, ApJ, 865, 117

Choban, C. R., Kereš, D., Hopkins, P. F., et al. 2022, MNRAS, 514, 4506

Churchwell, E., Babler, B. L., Meade, M. R., et al. 2009, PASP, 121, 213

Clark, S. E., Peek, J. E. G., & Putman, M. E. 2014, ApJ, 789, 82

Coles, P. 1988, MNRAS, 234, 509

Compiègne, M., Flagey, N., Noriega-Crespo, A., et al. 2010, ApJL, 724, L44

den Brok, J. S., Chatzigiannakis, D., Bigiel, F., et al. 2021, MNRAS, 504, 3221

Dobbs, C. L., Krumholz, M. R., Ballesteros-Paredes, J., et al. 2014, in Protostars and Planets VI, ed. H. Beuther et al. (Tucson, AZ: Univ. of Arizona Press), 3

Draine, B. T., Aniano, G., Krause, O., et al. 2014, ApJ, 780, 172

Draine, B. T., Dale, D. A., Bendo, G., et al. 2007, ApJ, 663, 866

Draine, B. T., & Li, A. 2007, ApJ, 657, 810

Dwek, E., Arendt, R. G., Fixsen, D. J., et al. 1997, ApJ, 475, 565

Eales, S., Smith, M. W. L., Auld, R., et al. 2012, ApJ, 761, 168

Ehlerová, S., & Palouš, J. 2005, A&A, 437, 101

Egorov, Kreckel, K., Sandstrom, K. M., et al. 2023, ApJL, 944, L16

Galliano, F., Galametz, M., & Jones, A. P. 2018, ARA&A, 56, 673

Giard, M., Lamarre, J. M., Pajot, F., & Serra, G. 1994, A&A, 286, 203

Giard, M., Pajot, F., Lamarre, J. M., et al. 1988, A&A, 201, L1

Glover, S. C. O., & Clark, P. C. 2012, MNRAS, 426, 377

Glover, S. C. O., & Mac Low, M.-M. 2007, ApJ, 659, 1317

Glover, S. C. O., & Smith, R. J. 2016, MNRAS, 462, 3011

Groves, B. A., Schinnerer, E., Leroy, A., et al. 2015, ApJ, 799, 96

Groves, B., Kreckel, K., Santoro, F., et al. 2023, arXiv:2301.03811

Hacar, A., Clark, S., Heitsch, F., et al. 2022, arXiv:2203.09562

Hassani, H., Rosolowsky, E. W., Leroy, A. K., et al. 2023, ApJL, 944, L21

Hensley, B. S., Murray, C. E., & Dodici, M. 2022, ApJ, 929, 23

Hoyer, N., Pinna, F., Kamlah, A. W. H., et al. 2023, ApJL, 944, L25

Jameson, K. E., Bolatto, A. D., Leroy, A. K., et al. 2016, ApJ, 825, 12

Jameson, K. E., Bolatto, A. D., Wolfire, M., et al. 2018, ApJ, 853, 111

Jenkins, E. B. 2009, ApJ, 700, 1299

Kalberla, P. M. W., & Haud, U. 2018, A&A, 619, A58

Kalberla, P. M. W., & Kerp, J. 2009, ARA&A, 47, 27

Kalberla, P. M. W., Kerp, J., & Haud, U. 2021, A&A, 654, A91

Koch, E. W., Rosolowsky, E. W., Boyden, R. D., et al. 2019, AJ, 158, 1

Koch, E. W., Rosolowsky, E. W., Leroy, A. K., et al. 2021, MNRAS, 504, 1801

Koch, E. W., Ward, C. G., Offner, S., Loeppky, J. L., & Rosolowsky, E. W. 2017, MNRAS, 471, 1506

Kowal, G., & Lazarian, A. 2007, ApJL, 666, L69

Krumholz, M. R., McKee, C. F., & Tumlinson, J. 2009, ApJ, 693, 216

Lang, P., Meidt, S. E., Rosolowsky, E., et al. 2020, ApJ, 897, 122

Lawton, B., Gordon, K. D., Babler, B., et al. 2010, ApJ, 716, 453

Lee, J. C., Sandstrom, K. M., Leroy, A. K., et al. 2023, ApJL, 944, L17

Lee, M.-Y., Stanimirović, S., Murray, C. E., Heiles, C., & Miller, J. 2015, ApJ, 809, 56

Leroy, A. K., Bolatto, A., Gordon, K., et al. 2011, ApJ, 737, 12

Leroy, A. K., Lee, C., Schruba, A., et al. 2013, ApJL, 769, L12

Leroy, A. K., Rosolowsky, E., Usero, A., et al. 2022, ApJ, 927, 149

Leroy, A. K., Sandstrom, K., Rosolowsky, E., et al. 2023, ApJL,

Leroy, A. K., Schinnerer, E., Hughes, A., et al. 2021, ApJS, 257, 43

Li, A., & Draine, B. T. 2001, ApJ, 554, 778

Li, C., Qiu, K., Hu, B., & Cao, Y. 2021, ApJL, 918, L2

Madden, S. C., Cormier, D., Hony, S., et al. 2020, A&A, 643, A141

Marchal, A., Miville-Deschênes, M.-A., Orieux, F., et al. 2019, A&A, 626, A101

Mathis, J. S., Mezger, P. G., & Panagia, N. 1983, A&A, 128, 212

Meidt, S. E., Rosolowsky, E., Sun, J., et al. 2023, ApJL, 944, L18

Peretto, N., & Fuller, G. A. 2009, A&A, 505, 405

Oschmann, J. M., Clampin, M., Fazio, G. G., et al. 2014, Proc. SPIE, 9143, 91433X

Pilyugin, L. S., & Grebel, E. K. 2016, MNRAS, 457, 3678

Pingel, N. M., Dempsey, J., McClure-Griffiths, N. M., et al. 2022, PASA, 39, e005

Puget, J. L., Leger, A., & Boulanger, F. 1985, A&A, 142, L19

Ragan, S. E., Henning, T., Tackenberg, J., et al. 2014, A&A, 568, A73

Roman-Duval, J., Gordon, K. D., Meixner, M., et al. 2014, ApJ, 797, 86

Roman-Duval, J., Jenkins, E. B., Tchernyshyov, K., et al. 2021, ApJ, 910, 95

Sandstrom, K. M., Leroy, A. K., Walter, F., et al. 2013, ApJ, 777, 5

Schruba, A., Bialy, S., & Sternberg, A. 2018, ApJ, 862, 110

Scoville, N., Aussel, H., Sheth, K., et al. 2014, ApJ, 783, 84

Smith, R. J., Glover, S. C. O., & Klessen, R. 2014, MNRAS, 445, 2900

Soler, J. D., Beuther, H., Syed, J., et al. 2020, A&A, 642, A163

Soler, J. D., Miville-Deschênes, M. A., Molinari, S., et al. 2022, A&A, 662, A96

Springel, V. 2010, MNRAS, 401, 791

Stanimirović, S., Murray, C. E., Lee, M.-Y., Heiles, C., & Miller, J. 2014, ApJ, 793, 132

Sternberg, A., Le Petit, F., Roueff, E., & Le Bourlot, J. 2014, ApJ, 790, 10

Suad, L. A., Caiafa, C. F., Arnal, E. M., & Cichowolski, S. 2014, A&A, 564, A116
Syed, J., Soler, J. D., Beuther, H., et al. 2022, A&A, 657, A1

Teng, Y.-H., Sandstrom, K. M., Sun, J., et al. 2022, ApJ, 925, 72

Teyssier, R. 2002, A&A, 385, 337

Thilkier, D. A., Lee, J. C., Deger, S., et al. 2023, ApJL, 944, L13

Treß, R. G., Sormani, M. C., Smith, R. J., et al. 2021, MNRAS, 505, 5438

Tress, R. G., Smith, R. J., Sormani, M. C., et al. 2020, MNRAS, 492, 2973

Vilchez, J. M., Relaño, M., Kennicutt, R., et al. 2019, MNRAS, 483, 4968

Walter, F., Brinks, E., de Blok, W. J. G., et al. 2008, AJ, 136, 2563

Watkins, E. J., Barnes, A., Henny, K. F., et al. 2023, ApJL, 944, L24

Whittet, D. C. B. 2010, ApJ, 710, 1009

Wolfire, M. G., Hollenbach, D., & McKee, C. F. 2010, ApJ, 716, 1191

Wolfire, M. G., McKee, C. F., Hollenbach, D., & Tielens, A. G. G. M. 2003, ApJ, 587, 278

Zucker, C., Battersby, C., & Goodman, A. 2018, ApJ, 864, 153
Syed, J., Soler, J. D., Beuther, H., et al. 2022, A&A, 657, A1

Teng, Y.-H., Sandstrom, K. M., Sun, J., et al. 2022, ApJ, 925, 72

Teyssier, R. 2002, A&A, 385, 337

Thilker, D. A., Lee, J. C., Deger, S., et al. 2023, ApJL, 944, L13

Treß, R. G., Sormani, M. C., Smith, R. J., et al. 2021, MNRAS, 505, 5438

Tress, R. G., Smith, R. J., Sormani, M. C., et al. 2020, MNRAS, 492, 2973

Vilchez, J. M., Relaño, M., Kennicutt, R., et al. 2019, MNRAS, 483, 4968

Walter, F., Brinks, E., de Blok, W. J. G., et al. 2008, AJ, 136, 2563

Watkins, E. J., Barnes, A., Henny, K. F., et al. 2023, ApJL, 944, L24

Whittet, D. C. B. 2010, ApJ, 710, 1009

Wolfire, M. G., Hollenbach, D., & McKee, C. F. 2010, ApJ, 716, 1191

Wolfire, M. G., McKee, C. F., Hollenbach, D., & Tielens, A. G. G. M. 2003, ApJ, 587, 278

Zucker, C., Battersby, C., & Goodman, A. 2018, ApJ, 864, 153