Showing 1–50 of 265 results for author: Dale, D A

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

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

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

Submitted 1 July, 2024; originally announced July 2024.

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

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

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

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

Submitted 17 June, 2024; originally announced June 2024.

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

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

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

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

Submitted 23 May, 2024; originally announced May 2024.

Comments: Accepted at ApJ, 39 pages, 25 figures

The Metallicity Dependence of PAH Emission in Galaxies I: Insights from Deep Radial Spitzer Spectroscopy

Authors: Cory M. Whitcomb, J. -D. T. Smith, Karin Sandstrom, Carl A. Starkey, Grant P. Donnelly, Bruce T. Draine, Evan D. Skillman, Daniel A. Dale, Lee Armus, Brandon S. Hensley, Thomas S. -Y. Lai, Robert C. Kennicutt

Abstract: We use deep Spitzer mid-infrared spectroscopic maps of radial strips across three nearby galaxies with well-studied metallicity gradients (M101, NGC 628, and NGC 2403) to explore the physical origins of the observed deficit of polycyclic aromatic hydrocarbons (PAHs) at sub-solar metallicity (i.e. the PAH-metallicity relation or PZR). These maps allow us to trace the evolution of all PAH features f… ▽ More We use deep Spitzer mid-infrared spectroscopic maps of radial strips across three nearby galaxies with well-studied metallicity gradients (M101, NGC 628, and NGC 2403) to explore the physical origins of the observed deficit of polycyclic aromatic hydrocarbons (PAHs) at sub-solar metallicity (i.e. the PAH-metallicity relation or PZR). These maps allow us to trace the evolution of all PAH features from 5-18 $μ$m as metallicity decreases continuously from solar ($Z_\odot$) to 0.2 $Z_\odot$. The total PAH to dust luminosity ratio remains relatively constant until reaching a threshold of $\sim$$\frac{2}{3}$$Z_\odot$, below which it declines smoothly but rapidly. The PZR has been attributed to preferential destruction of the smallest grains in the hard radiation environments found at low metallicity. In this scenario, a decrease in emission from the shortest wavelength PAH features is expected. In contrast, we find a steep decline in long wavelength power below $Z_\odot$, especially in the 17 $μ$m feature, with the shorter wavelength PAH bands carrying an increasingly large fraction of power at low metallicity. We use newly developed grain models to reproduce the observed PZR trends, including these variations in fractional PAH feature strengths. The model that best reproduces the data employs an evolving grain size distribution that shifts to smaller sizes as metallicity declines. We interpret this as a result of inhibited grain growth at low metallicity, suggesting continuous replenishment in the interstellar medium is the dominant process shaping the PAH grain population in galaxies. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 22 pages, 13 figures, submitted for publication in ApJ

Do spiral arms enhance star formation efficiency?

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

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

Submitted 8 May, 2024; originally announced May 2024.

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

JWST Observations of Starbursts: Cold Clouds and Plumes Launching in the M82 Outflow

Authors: Deanne B. Fisher, Alberto D. Bolatto, John Chisholm, Drummond Fielding, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Serena A. Cronin, Laura A. Lopez, J. D. Smith, Danielle A. Berg, Sebastian Lopez, Sylvain Veilleux, Paul P. van der Werf, Torsten Böker, Leindert A. Boogaard, Laura Lenkić, Simon C. O. Glover, Vicente Villanueva, Divakara Mayya, Thomas S. -Y. Lai, Daniel A. Dale, Kimberly L. Emig, Fabian Walter, Monica Relaño , et al. (6 additional authors not shown)

Abstract: In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority (… ▽ More In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority ($\sim$70\%) of PAH emission is associated with the plumes. We show that those structures contain smaller scale "clouds" with widths that are $\sim$5-15 pc, and they are morphologically similar to the results of "cloud-crushing" simulations. We estimate the cloud-crushing time-scales of $\sim$0.5-3 Myr, depending on assumptions. We show this time scale is consistent with a picture in which these observed PAH clouds survived break-out from the disk rather than being destroyed by the hot wind. The PAH emission in both the midplane and the outflow is shown to tightly correlate with that of Pa$α$ emission (from HST/NICMOS data), at the scale of both plumes and clouds, though the ratio of PAH-to-Pa$α$ increases at further distances from the midplane. Finally, we show that the outflow PAH emission is suppressed in regions of the M82 wind that are bright in X-ray emission. Overall, our results are broadly consistent with a picture in which cold gas in galactic outflows is launched via hierarchically structured plumes, and those small scale clouds are more likely to survive the wind environment when collected into the larger plume structure. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: Submitted to MNRAS

H-alpha emission and HII regions at the locations of recent supernovae in nearby galaxies

Authors: Ness Mayker Chen, Adam K. Leroy, Sumit K. Sarbadhicary, Laura A. Lopez, Todd A. Thompson, Ashley T. Barnes, Eric Emsellem, Brent Groves, Rupali Chandar, Mélanie Chevance, Ryan Chown, Daniel A. Dale, Oleg V. Egorov, Simon C. O. Glover, Kathryn Grasha, Ralf S. Klessen, Kathryn Kreckel, Jing Li, J. Eduardo Méndez-Delgado, Eric J. Murphy, Debosmita Pathak, Eva Schinnerer, David A. Thilker, Leonardo Úbeda, Thomas G. Williams

Abstract: We present a statistical analysis of the local, approximately 50-100 pc scale, H-alpha emission at the locations of recent (less than 125 years) supernovae (SNe) in nearby star-forming galaxies. Our sample consists of 32 SNe in 10 galaxies that are targets of the PHANGS-MUSE survey. We find that 41% (13/32) of these SNe occur coincident with a previously identified HII region. For comparison, HII… ▽ More We present a statistical analysis of the local, approximately 50-100 pc scale, H-alpha emission at the locations of recent (less than 125 years) supernovae (SNe) in nearby star-forming galaxies. Our sample consists of 32 SNe in 10 galaxies that are targets of the PHANGS-MUSE survey. We find that 41% (13/32) of these SNe occur coincident with a previously identified HII region. For comparison, HII regions cover 32% of the area within 1 kpc of any recent SN. Contrasting this local covering fraction with the fraction of SNe coincident with HII regions, we find a statistical excess of 7.6% +/- 8.7% of all SNe to be associated with HII regions. This increases to an excess of 19.2% +/- 10.4% when considering only core-collapse SNe. These estimates appear to be in good agreement with qualitative results from new, higher resolution HST H-alpha imaging, which also suggest many CCSNe detonate near but not in HII regions. Our results appear consistent with the expectation that only a modest fraction of stars explode during the first 5 Myr of the life of a stellar population, when H-alpha emission is expected to be bright. Of the HII region associated SNe, 8% (11/13) also have associated detected CO(2-1) emission, indicating the presence of molecular gas. The HII region associated SNe have typical Av extinctions approximately equal to 1 mag, consistent with a significant amount of pre-clearing of gas from the region before the SNe explode. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ; 33 pages, 13 figures, 3 tables in two-column AASTEX63 format

Monitoring AGNs with H$β$ Asymmetry. IV. First Reverberation Mapping Results of 14 AGNs

Authors: T. E. Zastrocky, Michael S. Brotherton, Pu Du, Jacob N. McLane, Kianna A. Olson, D. A. Dale, H. A. Kobulnicky, Jaya Maithil, My L. Nguyen, William T. Chick, David H. Kasper, Derek Hand, C. Adelman, Z. Carter, G. Murphree, M. Oeur, T. Roth, S. Schonsberg, M. J. Caradonna, J. Favro, A. J. Ferguson, I. M. Gonzalez, L. M. Hadding, H. D. Hagler, C. J. Rogers , et al. (19 additional authors not shown)

Abstract: We report first-time reverberation mapping results for 14 AGNs from the ongoing Monitoring AGNs with H$β$ Asymmetry campaign (MAHA). These results utilize optical spectra obtained with the Long Slit Spectrograph on the Wyoming Infrared 2.3m Telescope between 2017 November-2023 May. MAHA combines long-duration monitoring with high cadence. We report results from multiple observing seasons for 9 of… ▽ More We report first-time reverberation mapping results for 14 AGNs from the ongoing Monitoring AGNs with H$β$ Asymmetry campaign (MAHA). These results utilize optical spectra obtained with the Long Slit Spectrograph on the Wyoming Infrared 2.3m Telescope between 2017 November-2023 May. MAHA combines long-duration monitoring with high cadence. We report results from multiple observing seasons for 9 of the 14 objects. These results include H$β$ time lags, supermassive black hole masses, and velocity-resolved time lags. The velocity-resolved lags allow us to investigate the kinematics of the broad-line region. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 35 pages, 19 figures, accepted for publication in ApJ Supplement

PHANGS-HST: Globular Cluster Systems in 17 Nearby Spiral Galaxies

Authors: Matthew Floyd, Rupali Chandar, Bradley C. Whitmore, David A. Thilker, Janice C. Lee, Rachel E. Pauline, Zion L. Thomas, William J. Berschback, Kiana F. Henny, Daniel A. Dale, Ralf S. Klessen, Eva Schinnerer, Kathryn Grasha, Mederic Boquien, Kirsten L. Larson, Sinan Deger, Ashley T. Barnes, Adam K. Leroy, Erik Rosolowsky, Thomas G. Williams, Leonardo Ubeda

Abstract: We present new catalogs of likely globular clusters (GCs) in 17 nearby spiral galaxies studied as part of the PHANGS-HST Treasury Survey. The galaxies were imaged in five broad-band filters from the near-ultraviolet through the $I$ band. PHANGS-HST has produced catalogs of stellar clusters of all ages by selecting extended sources (from multiple concentration index measurements) followed by morpho… ▽ More We present new catalogs of likely globular clusters (GCs) in 17 nearby spiral galaxies studied as part of the PHANGS-HST Treasury Survey. The galaxies were imaged in five broad-band filters from the near-ultraviolet through the $I$ band. PHANGS-HST has produced catalogs of stellar clusters of all ages by selecting extended sources (from multiple concentration index measurements) followed by morphological classification (centrally concentrated and symmetric or asymmetric, multiple peaks, contaminant) by visually examining the V-band image and separately by a machine-learning algorithm which classified larger samples to reach fainter limits. From both cluster catalogs, we select an initial list of candidate GCs to have $B-V \geq 0.5$ and $V-I \geq 0.73$~mag, then remove likely contaminants (including reddened young clusters, background galaxies misclassified by the neural network, and chance superpositions/blends of stars) after a careful visual inspection. We find that $\approx86$ % of the color-selected candidates classified as spherically symmetric, and $\approx68$ of those classified as centrally concentrated but asymmetric are likely to be GCs. The luminosity functions of the GC candidates in 2 of our 17 galaxies, NGC 628 and NGC 3627, are atypical, and continue to rise at least 1~mag fainter than the expected turnover near $M_V \sim -7.4$. These faint candidate GCs have more extended spatial distributions than their bright counterparts, and may reside in the disk rather than the bulge/halo, similar to faint GCs previously discovered in M101. These faint clusters may be somewhat younger since the age-metallicity degeneracy makes it difficult to determine precise cluster ages from integrated colors once they reach $\approx1$~Gyr. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 21 pages, 12 figures

Journal ref: The Astronomical Journal, Volume 167, Issue 3, id.95, 16 pp. March 2024

PHANGS-HST catalogs for $\sim$100,000 star clusters and compact associations in 38 galaxies: I. Observed properties

Authors: Daniel Maschmann, Janice C. Lee, David A. Thilker, Bradley C. Whitmore, Sinan Deger, Mederic Boquien, Rupali Chandar, Daniel A. Dale, Aida Wofford, Stephen Hannon, Kirsten L. Larson, Adam K. Leroy, Eva Schinnerer, Erik W. Rosolowsky, Leonardo Ubeda, Ashley Barnes, Eric Emsellem, Kathryn Grasha, Brent Groves, Hwihyun Kim, Ralf S. Klessen, Kathryn Kreckel, Rebecca C. Levy, Francesca Pinna, Jimena Rodriguez , et al. (2 additional authors not shown)

Abstract: We present the largest catalog to-date of star clusters and compact associations in nearby galaxies. We have performed a V-band-selected census of clusters across the 38 spiral galaxies of the PHANGS-HST Treasury Survey, and measured integrated, aperture-corrected NUV-U-B-V-I photometry. This work has resulted in uniform catalogs that contain $\sim$20,000 clusters and compact associations which ha… ▽ More We present the largest catalog to-date of star clusters and compact associations in nearby galaxies. We have performed a V-band-selected census of clusters across the 38 spiral galaxies of the PHANGS-HST Treasury Survey, and measured integrated, aperture-corrected NUV-U-B-V-I photometry. This work has resulted in uniform catalogs that contain $\sim$20,000 clusters and compact associations which have passed human inspection and morphological classification, and a larger sample of $\sim$100,000 classified by neural network models. Here, we report on the observed properties of these samples, and demonstrate that tremendous insight can be gained from just the observed properties of clusters, even in the absence of their transformation into physical quantities. In particular, we show the utility of the UBVI color-color diagram, and the three principal features revealed by the PHANGS-HST cluster sample: the young cluster locus, the middle-age plume, and the old globular cluster clump. We present an atlas of maps of the 2D spatial distribution of clusters and compact associations in the context of the molecular clouds from PHANGS-ALMA. We explore new ways of understanding this large dataset in a multi-scale context by bringing together once-separate techniques for the characterization of clusters (color-color diagrams and spatial distributions) and their parent galaxies (galaxy morphology and location relative to the galaxy main sequence). A companion paper presents the physical properties: ages, masses, and dust reddenings derived using improved spectral energy distribution (SED) fitting techniques. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 48 pages, 26 figures, article in review at APJ

The Impact of an AGN on PAH Emission in Galaxies: the Case of Ring Galaxy NGC 4138

Authors: G. P. Donnelly, J. D. T. Smith, B. T. Draine, A. Togi, T. S. -Y. Lai, L. Armus, D. A. Dale, V. Charmandaris

Abstract: We present a focused study of radially-resolved varying PAH emission in the low-luminosity AGN-host NGC 4138 using deep Spitzer/IRS spectral maps. Using new model PAH spectra, we investigate whether these variations could be associated with changes to the PAH grain size distribution due to photodestruction by the AGN. Separately, we model the effects of the varying radiation field within NGC 4138,… ▽ More We present a focused study of radially-resolved varying PAH emission in the low-luminosity AGN-host NGC 4138 using deep Spitzer/IRS spectral maps. Using new model PAH spectra, we investigate whether these variations could be associated with changes to the PAH grain size distribution due to photodestruction by the AGN. Separately, we model the effects of the varying radiation field within NGC 4138, and we use this model to predict the corresponding changes in the PAH emission spectrum. We find that PAH band ratios are strongly variable with radius in this galaxy with short-to-long wavelength band ratios peaking in the starburst ring. The changing mix of starlight appears to have a considerable effect on the trends in these band ratios, and our radiation model predicts the shapes of these trends. However, the amplitude of observed variation is ~2.5 times larger than predicted for some ratios. A cutoff of small grains in the PAH size distribution, as has been suggested for AGN, together with changes in PAH ionization fraction could explain the behavior of the shorter bands, but this model fails to reproduce longer band behaviors. Additionally, we find that short-to-long wavelength PAH band ratios increase slightly within ~270pc of the center, suggesting that the AGN may directly influence PAH emission there. △ Less

Submitted 23 April, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: 19 pages, 13 figures, published in ApJ. Updated values in radiation model to the same as the published version

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

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

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

Submitted 6 February, 2024; originally announced February 2024.

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

JWST Observations of Starbursts: Polycyclic Aromatic Hydrocarbon Emission at the Base of the M 82 Galactic Wind

Authors: Alberto D. Bolatto, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Deanne B. Fisher, Adam K. Leroy, Serena A. Cronin, Ralf S. Klessen, J. D. Smith, Dannielle A. Berg, Torsten Boeker, Leindert A. Boogaard, Eve C. Ostriker, Todd A. Thompson, Juergen Ott, Laura Lenkic, Laura A. Lopez, Daniel A. Dale, Sylvain Veilleux, Paul P. van der Werf, Simon C. O. Glover, Karin M. Sandstrom, Evan D. Skillman, John Chisholm, Vicente Villanueva , et al. (15 additional authors not shown)

Abstract: We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find promine… ▽ More We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find prominent plumes of PAH emission extending outward from the central starburst region, together with a network of complex filamentary substructure and edge-brightened bubble-like features. The structure of the PAH emission closely resembles that of the ionized gas, as revealed in Paschen alpha and free-free radio emission. We discuss the origin of the structure, and suggest the PAHs are embedded in a combination of neutral, molecular, and photoionized gas. △ Less

Submitted 21 April, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted to The Astrophysical Journal

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

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

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

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

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

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

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

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

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

Comments: 27 pages, 12 figures; ApJ accepted

ALMA-LEGUS II: The Influence of Sub-Galactic Environment on Molecular Cloud Properties

Authors: Molly K. Finn, Kelsey E. Johnson, Remy Indebetouw, Allison H. Costa, Angela Adamo, Alessandra Aloisi, Lauren Bittle, Daniela Calzetti, Daniel A. Dale, Clare L. Dobbs, Jennifer Donovan Meyer, Bruce G. Elmegreen, Debra M. Elmegreen, Michele Fumagalli, J. S. Gallagher, Kathryn Grasha, Eva K. Grebel, Robert C. Kennicutt, Mark R. Krumholz, Janice C. Lee, Matteo Messa, Preethi Nair, Elena Sabbi, Linda J. Smith, David A. Thilker , et al. (2 additional authors not shown)

Abstract: We compare the molecular cloud properties in sub-galactic regions of two galaxies, barred spiral NGC 1313, which is forming many massive clusters, and flocculent spiral NGC 7793, which is forming significantly fewer massive clusters despite having a similar star formation rate to NGC 1313. We find that there are larger variations in cloud properties between different regions within each galaxy tha… ▽ More We compare the molecular cloud properties in sub-galactic regions of two galaxies, barred spiral NGC 1313, which is forming many massive clusters, and flocculent spiral NGC 7793, which is forming significantly fewer massive clusters despite having a similar star formation rate to NGC 1313. We find that there are larger variations in cloud properties between different regions within each galaxy than there are between the galaxies on a global scale, especially for NGC 1313. There are higher masses, linewidths, pressures, and virial parameters in the arms of NGC 1313 and center of NGC 7793 than in the interarm and outer regions of the galaxies. The massive cluster formation of NGC 1313 may be driven by its greater variation in environments, allowing more clouds with the necessary conditions to arise, although no one parameter seems primarily responsible for the difference in star formation. Meanwhile NGC 7793 has clouds that are as massive and have as much kinetic energy as clouds in the arms of NGC 1313, but have densities and pressures more similar to the interarm regions and so are less inclined to collapse and form stars. The cloud properties in NGC 1313 and NGC 7793 suggest that spiral arms, bars, interarm regions, and flocculent spirals each represent distinct environments with regard to molecular cloud populations. We see surprisingly little difference in surface densities between the regions, suggesting that the differences in surface densities frequently seen between arm and interarm regions of lower-resolution studies are indicative of the sparsity of molecular clouds, rather than differences in their true surface density. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 24 pages, 12 figures, accepted for publication in ApJ

ALMA-LEGUS I: The Influence of Galaxy Morphology on Molecular Cloud Properties

Authors: Molly K. Finn, Kelsey E. Johnson, Remy Indebetouw, Allison H. Costa, Angela Adamo, Alessandra Aloisi, Lauren Bittle, Daniela Calzetti, Daniel A. Dale, Clare L. Dobbs, Jennifer Donovan Meyer, Bruce G. Elmegreen, Debra M. Elmegreen, Michele Fumagalli, J. S. Gallagher, Kathryn Grasha, Eva K. Grebel, Robert C. Kennicutt, Mark R. Krumholz, Janice C. Lee, Matteo Messa, Preethi Nair, Elena Sabbi, Linda J. Smith, David A. Thilker , et al. (2 additional authors not shown)

Abstract: We present a comparative study of the molecular gas in two galaxies from the LEGUS sample: barred spiral NGC 1313 and flocculent spiral NGC 7793. These two galaxies have similar masses, metallicities, and star formation rates, but NGC 1313 is forming significantly more massive star clusters than NGC 7793, especially young massive clusters (<10 Myr, >10^4 Msol). Using ALMA CO(2-1) observations of t… ▽ More We present a comparative study of the molecular gas in two galaxies from the LEGUS sample: barred spiral NGC 1313 and flocculent spiral NGC 7793. These two galaxies have similar masses, metallicities, and star formation rates, but NGC 1313 is forming significantly more massive star clusters than NGC 7793, especially young massive clusters (<10 Myr, >10^4 Msol). Using ALMA CO(2-1) observations of the two galaxies with the same sensitivities and resolutions of 13 pc, we directly compare the molecular gas in these two similar galaxies to determine the physical conditions responsible for their large disparity in cluster formation. By fitting size-linewidth relations for the clouds in each galaxy, we find that NGC 1313 has a higher intercept than NGC 7793, implying that its clouds have higher kinetic energies at a given size scale. NGC 1313 also has more clouds near virial equilibrium than NGC 7793, which may be connected to its higher rate of massive cluster formation. However, these virially bound clouds do not show a stronger correlation with young clusters than that of the general cloud population. We find surprisingly small differences between the distributions of molecular cloud populations in the two galaxies, though the largest of those differences are that NGC 1313 has higher surface densities and lower free-fall times. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 25 pages, 10 figures, accepted for publication in ApJ

Surveying the Whirlpool at Arcseconds with NOEMA (SWAN)- I. Mapping the HCN and N$_2$H$^+$ 3mm lines

Authors: Sophia K. Stuber, Jerome Pety, Eva Schinnerer, Frank Bigiel, Antonio Usero, Ivana Beslić, Miguel Querejeta, María J. Jiménez-Donaire, Adam Leroy, Jakob den Brok, Lukas Neumann, Cosima Eibensteiner, Yu-Hsuan Teng, Ashley Barnes, Mélanie Chevance, Dario Colombo, Daniel A. Dale, Simon C. O. Glover, Daizhong Liu, Hsi-An Pan

Abstract: We present the first results from "Surveying the Whirlpool at Arcseconds with NOEMA" (SWAN), an IRAM Northern Extended Millimetre Array (NOEMA)+30m large program that maps emission from several molecular lines at 90 and 110 GHz in the iconic nearby grand-design spiral galaxy M~51 at cloud-scale resolution ($\sim$3\arcsec=125\,pc). As part of this work, we have obtained the first sensitive cloud-sc… ▽ More We present the first results from "Surveying the Whirlpool at Arcseconds with NOEMA" (SWAN), an IRAM Northern Extended Millimetre Array (NOEMA)+30m large program that maps emission from several molecular lines at 90 and 110 GHz in the iconic nearby grand-design spiral galaxy M~51 at cloud-scale resolution ($\sim$3\arcsec=125\,pc). As part of this work, we have obtained the first sensitive cloud-scale map of N$_2$H$^+$(1-0) of the inner $\sim5\,\times 7\,$kpc of a normal star-forming galaxy, which we compare to HCN(1-0) and CO(1-0) emission to test their ability in tracing dense, star-forming gas. The average N$_2$H$^+$-to-HCN line ratio of our total FoV is $0.20\pm0.09$, with strong regional variations of a factor of $\gtrsim 2$ throughout the disk, including the south-western spiral arm and the center. The central $\sim1\,$kpc exhibits elevated HCN emission compared to N$_2$H$^+$, probably caused by AGN-driven excitation effects. We find that HCN and N$_2$H$^+$ are strongly super-linearily correlated in intensity ($ρ_\mathrm{Sp}\sim 0.8$), with an average scatter of $\sim0.14\,$dex over a span of $\gtrsim 1.5\,$dex in intensity. When excluding the central region, the data is best described by a power-law of exponent $1.2$, indicating that there is more N$_2$H$^+$ per unit HCN in brighter regions. Our observations demonstrate that the HCN-to-CO line ratio is a sensitive tracer of gas density in agreement with findings of recent Galactic studies which utilize N$_2$H$^+$. The peculiar line ratios present near the AGN and the scatter of the power-law fit in the disk suggest that in addition to a first-order correlation with gas density, second-order physics (such as optical depth, gas temperature) or chemistry (abundance variations) are encoded in the N$_2$H$^+$/CO, HCN/CO and N$_2$H$^+$/HCN ratios. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: Accepted for publication in A&A. 6 pages, 3 figures (+ Appendix 4 pages, 2 figures)

The PHANGS-AstroSat Atlas of Nearby Star Forming Galaxies

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

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

Submitted 10 December, 2023; originally announced December 2023.

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

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

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

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

Submitted 29 November, 2023; originally announced November 2023.

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

The spatially resolved star formation history of the dwarf spiral galaxy NGC 5474

Authors: G. Bortolini, M. Cignoni, E. Sacchi, M. Tosi, F. Annibali, R. Pascale, M. Bellazzini, D. Calzetti, A. Adamo, Daniel. A. Dale, M. Fumagalli, John. S. Gallagher, K. Grasha, Kelsey E. Johnson, Sean. T. Linden, M. Messa, G. Östlin, E. Sabbi, A. Wofford

Abstract: We study the resolved stellar populations and derive the star formation history of NGC 5474, a peculiar star-forming dwarf galaxy at a distance of $\sim 7$ Mpc, using Hubble Space Telescope Advanced Camera for Surveys data from the Legacy Extragalactic UV Survey (LEGUS) program. We apply an improved colour-magnitude diagram fitting technique based on the code SFERA and use the latest PARSEC-COLIBR… ▽ More We study the resolved stellar populations and derive the star formation history of NGC 5474, a peculiar star-forming dwarf galaxy at a distance of $\sim 7$ Mpc, using Hubble Space Telescope Advanced Camera for Surveys data from the Legacy Extragalactic UV Survey (LEGUS) program. We apply an improved colour-magnitude diagram fitting technique based on the code SFERA and use the latest PARSEC-COLIBRI stellar models. Our results are the following. The off-centre bulge-like structure, suggested to constitute the bulge of the galaxy, is dominated by star formation (SF) activity initiated $14$ Gyr ago and lasted at least up to $1$ Gyr ago. Nevertheless, this component shows clear evidence of prolonged SF activity (lasting until $\sim 10$ Myr ago). We estimate the total stellar mass of the bulge-like structure to be $(5.0 \pm 0.3) \times 10^{8}$ \MSUN. Such a mass is consistent with published suggestions that this structure is in fact an independent system orbiting around and not within NGC 5474's disc. The stellar over-density located to the South-West of the bulge-like structure shows a significant SF event older than $1$ Gyr, while it is characterised by two recent peaks of SF, around $\sim10$ and $\sim100$ Myr ago. In the last Gyr, the behavior of the stellar disc is consistent with what is known in the literature as `gasping'. The synchronised burst at $10-35$ Myr in all components might hint to the recent gravitational interaction between the stellar bulge-like structure and the disc of NGC 5474. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: 18 pages, 12 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Extragalactic Magnetism with SOFIA (SALSA Legacy Program). VII. A Tomographic View of Far-infrared and Radio Polarimetric Observations through MHD Simulations of Galaxies

Authors: Sergio Martin-Alvarez, Enrique Lopez-Rodriguez, Tara Dacunha, Susan E. Clark, Alejandro S. Borlaff, Rainer Beck, Francisco Rodríguez Montero, S. Lyla Jung, Julien Devriendt, Adrianne Slyz, Julia Roman-Duval, Evangelia Ntormousi, Mehrnoosh Tahani, Kandaswamy Subramanian, Daniel A. Dale, Pamela M. Marcum, Konstantinos Tassis, Ignacio del Moral-Castro, Le Ngoc Tram, Matt J. Jarvis

Abstract: The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy t… ▽ More The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy to identify and quantify the regions, scales, and interstellar medium (ISM) phases probed at FIR and radio wavelengths. Our studied suite of magnetohydrodynamical cosmological zoom-in simulations features high-resolutions (10 pc full-cell size) and multiple magnetization models. Our synthetic observations have a striking resemblance to those of observed galaxies. We find that the total and polarized radio emission extends to approximately double the altitude above the galactic disk (half-intensity disk thickness of $h_\text{I radio} \sim h_\text{PI radio} = 0.23 \pm 0.03$ kpc) relative to the total FIR and polarized emission that are concentrated in the disk midplane ($h_\text{I FIR} \sim h_\text{PI FIR} = 0.11 \pm 0.01$ kpc). Radio emission traces magnetic fields at scales of $\gtrsim 300$ pc, whereas FIR emission probes magnetic fields at the smallest scales of our simulations. These scales are comparable to our spatial resolution and well below the spatial resolution ($<300$ pc) of existing FIR polarimetric measurements. Finally, we confirm that synchrotron emission traces a combination of the warm neutral and cold neutral gas phases, whereas FIR emission follows the densest gas in the cold neutral phase in the simulation. These results are independent of the ISM magnetic field strength. The complementarity we measure between radio and FIR wavelengths motivates future multiwavelength polarimetric observations to advance our knowledge of extragalactic magnetism. △ Less

Submitted 24 March, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: Accepted for publication in the ApJ. Replaced to match final version. 35 pages, 17 figures, 4 tables

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

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

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

Submitted 27 October, 2023; originally announced October 2023.

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

Halfway to the peak: Spatially resolved star formation and kinematics in a z=0.54 dusty galaxy with JWST/MIRI

Authors: Jason Young, Alexandra Pope, Anna Sajina, Lin Yan, Thiago S Goncalves, Miriam Eleazer, Stacey Alberts, Lee Armus, Matteo Bonato, Daniel A. Dale, Duncan Farrah, Carl Ferkinhoff, Christopher C. Hayward, Jed McKinney, Eric J. Murphy, Nicole Nesvadba, Patrick Ogle, Leonid Sajkov, Sylvain Veilleux

Abstract: We present JWST/MIRI/MRS observations of an infrared luminous disk galaxy, FLS1, at z=0.54. With a lookback time of 5 Gyr, FLS1 is chronologically at the midpoint between the peak epoch of star formation and the present day. The MRS data provide maps of the atomic fine structure lines [Ar II]6.99 micron, [Ar III]8.99 micron, [Ne II]12.81 micron, and [Ne III]15.55 micron, polycyclic aromatic hydroc… ▽ More We present JWST/MIRI/MRS observations of an infrared luminous disk galaxy, FLS1, at z=0.54. With a lookback time of 5 Gyr, FLS1 is chronologically at the midpoint between the peak epoch of star formation and the present day. The MRS data provide maps of the atomic fine structure lines [Ar II]6.99 micron, [Ar III]8.99 micron, [Ne II]12.81 micron, and [Ne III]15.55 micron, polycyclic aromatic hydrocarbon (PAH) features at 3.3 micron, 6.2 micron, and 11.3 micron, and the warm molecular gas indicators H2S(5) and H2S(3); all these emission features are spatially resolved. We find that the PAH emission is more extended along the Northern side of the galaxy when compared to the well-studied star-formation tracer [Ne II]. The H2 rotational lines, which are shock indicators, are strongest and most extended on the Southern side of the galaxy. [Ar II] is the second brightest fine structure line detected in FLS1 and we show that it is a useful kinematic probe which can be detected with JWST out to z=3. Velocity maps of [Ar II] show a rotating disk with signs of turbulence. Our results provide an example of how spatially resolved mid-infrared spectroscopy can allow us to better understand the star formation and ISM conditions in a galaxy halfway back to the peak epoch of galaxy evolution. △ Less

Submitted 12 October, 2023; v1 submitted 10 October, 2023; originally announced October 2023.

Comments: 10 pages, 5 figures, 1 table

Investigating the Drivers of Electron Temperature Variations in HII Regions with Keck-KCWI and VLT-MUSE

Authors: Ryan J. Rickards Vaught, Karin M. Sandstrom, Francesco Belfiore, Kathryn Kreckel, J. Eduardo Méndez-Delgado, Eric Emsellem, Brent Groves, Guillermo A. Blanc, Daniel A. Dale, Oleg V. Egorov, Simon C. O. Glover, Kathryn Grasha, Ralf S. Klessen, Justus Neumann, Thomas G. Williams

Abstract: HII region electron temperatures are a critical ingredient in metallicity determinations and recent observations reveal systematic variations in the temperatures measured using different ions. We present electron temperatures ($T_e$) measured using the optical auroral lines ([NII]$\lambda5756$, [OII]$λ\lambda7320,7330$, [SII]$λ\lambda4069,4076$, [OIII]$\lambda4363$, and [SIII]$\lambda6312$) for a… ▽ More HII region electron temperatures are a critical ingredient in metallicity determinations and recent observations reveal systematic variations in the temperatures measured using different ions. We present electron temperatures ($T_e$) measured using the optical auroral lines ([NII]$\lambda5756$, [OII]$λ\lambda7320,7330$, [SII]$λ\lambda4069,4076$, [OIII]$\lambda4363$, and [SIII]$\lambda6312$) for a sample of HII regions in seven nearby galaxies. We use observations from the Physics at High Angular resolution in Nearby Galaxies survey (PHANGS) obtained with integral field spectrographs on Keck (Keck Cosmic Web Imager; KCWI) and the Very Large Telescope (Multi-Unit Spectroscopic Explorer; MUSE). We compare the different $T_e$ measurements with HII region and interstellar medium environmental properties such as electron density, ionization parameter, molecular gas velocity dispersion, and stellar association/cluster mass and age obtained from PHANGS. We find that the temperatures from [OII] and [SII] are likely over-estimated due to the presence of electron density inhomogeneities in HII regions. We observe that differences between [NII] and [SIII] temperatures are weakly correlated with stellar association mass and molecular gas velocity dispersion. We measure high [OIII] temperatures in a subset of regions with high molecular gas velocity dispersion and low ionization parameter, which may be explained by the presence of low-velocity shocks. In agreement with previous studies, the $T_{\rm{e}}$--$T_{\rm{e}}$ between [NII] and [SIII] temperatures have the lowest observed scatter and generally follow predictions from photoionization modeling, which suggests that these tracers reflect HII region temperatures across the various ionization zones better than [OII], [SII], and [OIII]. △ Less

Submitted 7 March, 2024; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: Accepted to ApJ, 66 pages, 36 figures

A constant N$_2$H$^+$(1-0)-to-HCN(1-0) ratio on kiloparsec scales

Authors: M. J. Jiménez-Donaire, A. Usero, I. Bešlić, M. Tafalla, A. Chacón-Tanarro, Q. Salomé, C. Eibensteiner, A. García-Rodríguez, A. Hacar, A. T. Barnes, F. Bigiel, M. Chevance, D. Colombo, D. A. Dale, T. A. Davis, S. C. O. Glover, J. Kauffmann, R. S. Klessen, A. K. Leroy, L. Neumann, H. Pan, J. Pety, M. Querejeta, T. Saito, E. Schinnerer , et al. (2 additional authors not shown)

Abstract: Nitrogen hydrides such as NH$_3$ and N$_2$H$^+$ are widely used by Galactic observers to trace the cold dense regions of the interstellar medium. In external galaxies, because of limited sensitivity, HCN has become the most common tracer of dense gas over large parts of galaxies. We provide the first systematic measurements of N$_2$H$^+$(1-0) across different environments of an external spiral gal… ▽ More Nitrogen hydrides such as NH$_3$ and N$_2$H$^+$ are widely used by Galactic observers to trace the cold dense regions of the interstellar medium. In external galaxies, because of limited sensitivity, HCN has become the most common tracer of dense gas over large parts of galaxies. We provide the first systematic measurements of N$_2$H$^+$(1-0) across different environments of an external spiral galaxy, NGC6946. We find a strong correlation ($r>0.98,p<0.01$) between the HCN(1-0) and N$_2$H$^+$(1-0) intensities across the inner $\sim8\mathrm{kpc}$ of the galaxy, at kiloparsec scales. This correlation is equally strong between the ratios N$_2$H$^+$(1-0)/CO(1-0) and HCN(1-0)/CO(1-0), tracers of dense gas fractions ($f_\mathrm{dense}$). We measure an average intensity ratio of N$_2$H$^+$(1-0)/HCN(1-0)$=0.15\pm0.02$ over our set of five IRAM-30m pointings. These trends are further supported by existing measurements for Galactic and extragalactic sources. This narrow distribution in the average ratio suggests that the observed systematic trends found in kiloparsec-scale extragalactic studies of $f_\mathrm{dense}$ and the efficiency of dense gas (SFE$_\mathrm{dense}$) would not change if we employed N$_2$H$^+$(1-0) as a more direct tracer of dense gas. At kiloparsec scales our results indicate that the HCN(1-0) emission can be used to predict the expected N$_2$H$^+$(1-0) over those regions. Our results suggest that, even if HCN(1-0) and N$_2$H$^+$(1-0) trace different density regimes within molecular clouds, subcloud differences average out at kiloparsec scales, yielding the two tracers proportional to each other. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

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

Star Cluster Classification using Deep Transfer Learning with PHANGS-HST

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

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

Submitted 27 July, 2023; originally announced July 2023.

Comments: 16 pages, 10 figures

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

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

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

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

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

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

Calibrating mid-infrared emission as a tracer of obscured star formation on HII-region scales in the era of JWST

Authors: Francesco Belfiore, Adam K. Leroy, Thomas G. Williams, Ashley T. Barnes, Frank Bigiel, Médéric Boquien, Yixian Cao, Jérémy Chastenet, Enrico Congiu, Daniel A. Dale, Oleg V. Egorov, Cosima Eibensteiner, Eric Emsellem, Simon C. O. Glover, Brent Groves, Hamid Hassani, Ralf S. Klessen, Kathryn Kreckel, Lukas Neumann, Justus Neumann, Miguel Querejeta, Erik Rosolowsky, Patricia Sanchez-Blazquez, Karin Sandstrom, Eva Schinnerer , et al. (3 additional authors not shown)

Abstract: Measurements of the star formation activity on cloud scales are fundamental to uncovering the physics of the molecular cloud, star formation, and stellar feedback cycle in galaxies. Infrared (IR) emission from small dust grains and polycyclic aromatic hydrocarbons (PAHs) are widely used to trace the obscured component of star formation. However, the relation between these emission features and dus… ▽ More Measurements of the star formation activity on cloud scales are fundamental to uncovering the physics of the molecular cloud, star formation, and stellar feedback cycle in galaxies. Infrared (IR) emission from small dust grains and polycyclic aromatic hydrocarbons (PAHs) are widely used to trace the obscured component of star formation. However, the relation between these emission features and dust attenuation is complicated by the combined effects of dust heating from old stellar populations and an uncertain dust geometry with respect to heating sources. We use images obtained with NIRCam and MIRI as part of the PHANGS--JWST survey to calibrate dust emission at 21$\rm μm$, and the emission in the PAH-tracing bands at 3.3, 7.7, 10, and 11.3$\rm μm$ as tracers of obscured star formation. We analyse $\sim$ 20000 optically selected HII regions across 19 nearby star-forming galaxies, and benchmark their IR emission against dust attenuation measured from the Balmer decrement. We model the extinction-corrected H$α$ flux as the sum of the observed H$α$ emission and a term proportional to the IR emission, with $a_{IR}$ as the proportionality coefficient. A constant $a_{IR}$ leads to extinction-corrected H$α$ estimates which agree with those obtained with the Balmer decrement with a scatter of $\sim$ 0.1 dex for all bands considered. Among these bands, 21$\rm μm$ emission is demonstrated to be the best tracer of dust attenuation. The PAH-tracing bands underestimate the correction for bright HII regions, since in these environments the ratio of PAH-tracing bands to 21$\rm μm$ decreases, signalling destruction of the PAH molecules. For fainter HII regions all bands suffer from an increasing contamination from the diffuse infrared background. △ Less

Submitted 1 September, 2023; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: accepted for publication in A&A

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

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

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

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

Submitted 5 May, 2023; originally announced May 2023.

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

Spectral Energy Distributions for 258 Local Volume Galaxies

Authors: Daniel A. Dale, Mederic Boquien, Jordan A. Turner, Daniela Calzetti, Robert C. Kennicutt Jr., Janice C. Lee

Abstract: We present model spectral energy distribution (SED) fits to ultraviolet/optical/infrared observations for the 258 nearby galaxies in the Local Volume Legacy survey, a sample dominated by lower-luminosity dwarf irregular systems. The data for each galaxy include up to 26 spatially-integrated broadband and narrowband fluxes from the Galaxy Evolution Explorer, Spitzer Space Telescope, and Infrared As… ▽ More We present model spectral energy distribution (SED) fits to ultraviolet/optical/infrared observations for the 258 nearby galaxies in the Local Volume Legacy survey, a sample dominated by lower-luminosity dwarf irregular systems. The data for each galaxy include up to 26 spatially-integrated broadband and narrowband fluxes from the Galaxy Evolution Explorer, Spitzer Space Telescope, and Infrared Astronomical Satellite space-based platforms and from the Sloan Digital Sky Survey, Two Micron All Sky Survey, and other ground-based efforts. The CIGALE SED fitting package is employed using a delayed star formation history with an optional late burst or quenching episode to constrain 11 different free parameters that characterize the properties of each galaxy's stellar and dust emission, with the overriding constraint that the ultraviolet/optical emission absorbed by interstellar dust grains is emitted in equal energy portions at infrared wavelengths. The main results are: i) 94% of the SED fits yield reduced chi^2 values less than 3; ii) the modeled stellar masses agree with those derived from 3.6um-based measures with a scatter of 0.07 dex; iii) for a typical galaxy in the sample the SED-derived star formation rate averaged over the past 100 Myr is about 88% of the value derived from standard hybrid indicators on similar timescales; and iv) there is a statistically significant inverse relation between the stellar mass fraction appearing in the late burst and the total stellar mass. These results build upon prior SED modeling efforts in the local volume and lay the groundwork for future studies of more distant low-metallicity galaxies with JWST. △ Less

Submitted 24 April, 2023; originally announced April 2023.

Comments: Accepted for publication in AJ

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

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

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

Submitted 4 April, 2023; originally announced April 2023.

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

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

Resolved stellar population properties of PHANGS-MUSE galaxies

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

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

Submitted 23 March, 2023; originally announced March 2023.

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

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

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- V: First results on the magnetic field orientation of galaxies

Authors: Alejandro S. Borlaff, Enrique Lopez-Rodriguez, Rainer Beck, Susan E. Clark, Evangelia Ntormousi, Konstantinos Tassis, Sergio Martin-Alvarez, Mehrnoosh Tahani, Daniel A. Dale, Ignacio del Moral Castro, Julia Roman-Duval, Pamela M. Marcum, John E. Beckman, Kandaswamy Subramanian, Sarah Eftekharzadeh, Leslie Proudfit

Abstract: We present the analysis of the magnetic field ($B$-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby ($<20$ Mpc) galaxies with resolved (5 arcsec-18 arcsec; $90$ pc--$1$ kpc) imaging polarimetric observations using HAWC+/SOFIA from… ▽ More We present the analysis of the magnetic field ($B$-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby ($<20$ Mpc) galaxies with resolved (5 arcsec-18 arcsec; $90$ pc--$1$ kpc) imaging polarimetric observations using HAWC+/SOFIA from $53$ to $214$ \um. We compute the magnetic pitch angle ($Ψ_{B}$) profiles as a function of the galactrocentric radius. We introduce a new magnetic alignment parameter ($ζ$) to estimate the disordered-to-ordered $B$-field ratio in spiral $B$-fields. We find FIR and radio wavelengths to not generally trace the same $B$-field morphology in galaxies. The $Ψ_{B}$ profiles tend to be more ordered with galactocentric radius in radio ($ζ_{\rm{6cm}} = 0.93\pm0.03$) than in FIR ($ζ_{\rm{154μm}} = 0.84\pm0.14$). For spiral galaxies, FIR $B$-fields are $2-75$\% more turbulent than the radio $B$-fields. For starburst galaxies, we find that FIR polarization is a better tracer of the $B$-fields along the galactic outflows than radio polarization. Our results suggest that the $B$-fields associated with dense, dusty, turbulent star-forming regions, those traced at FIR, are less ordered than warmer, less-dense regions, those traced at radio, of the interstellar medium. The FIR $B$-fields seem to be more sensitive to the activity of the star-forming regions and the morphology of the molecular clouds within a vertical height of few hundred pc in the disk of spiral galaxies than the radio $B$-fields. △ Less

Submitted 13 June, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

Comments: 28 pages, 14 figures. Accepted for publication in The Astrophysical Journal

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

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

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

Submitted 21 March, 2023; originally announced March 2023.

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

Star Formation Laws and Efficiencies across 80 Nearby Galaxies

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

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

Submitted 23 February, 2023; originally announced February 2023.

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

Journal ref: ApJL, 945, L19 (2023)

Quantifying the energetics of molecular superbubbles in PHANGS galaxies

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

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

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

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

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

PHANGS-MUSE: Detection and Bayesian classification of ~40000 ionised nebulae in nearby spiral galaxies

Authors: Enrico Congiu, Guillermo A. Blanc, Francesco Belfiore, Francesco Santoro, Fabian Scheuermann, Kathryn Kreckel, Eric Emsellem, Brent Groves, Hsi-An Pan, Frank Bigiel, Daniel A. Dale, Simon C. O. Glover, Kathryn Grasha, Oleg V. Egorov, Adam Leroy, Eva Schinnerer, Elizabeth J. Watkins, Thomas G. Williams

Abstract: In this work, we present a new catalogue of >40000 ionised nebulae distributed across the 19 galaxies observed by the PHANGS-MUSE survey. The nebulae have been classified using a new model-comparison-based algorithm that exploits the odds ratio principle to assign a probabilistic classification to each nebula in the sample. The resulting catalogue is the largest catalogue containing complete spect… ▽ More In this work, we present a new catalogue of >40000 ionised nebulae distributed across the 19 galaxies observed by the PHANGS-MUSE survey. The nebulae have been classified using a new model-comparison-based algorithm that exploits the odds ratio principle to assign a probabilistic classification to each nebula in the sample. The resulting catalogue is the largest catalogue containing complete spectral and spatial information for a variety of ionised nebulae available so far in the literature. We developed this new algorithm to address some of the limitations of the traditional classification criteria, such as their binarity, the sharpness of the involved limits, and the limited amount of data they rely on for the classification. The analysis of the catalogue shows that the algorithm performs well when selecting H II regions. We can recover their luminosity function, and its properties are in line with what is available in the literature. We also identify a rather significant population of shock-ionised regions (mostly composed of supernova remnants), an order of magnitude larger than any other homogeneous catalogue of supernova remnants currently available in the literature. The number of supernova remnants we identify per galaxy is in line with results in our Galaxy and other very nearby sources. However, limitations in the source detection algorithm result in an incomplete sample of planetary nebulae, even though their classification seems robust. Finally, we demonstrate how applying a correction for the contribution of the diffuse ionised gas to the nebulae's spectra is essential to obtain a robust classification of the objects and how a correct measurement of the extinction using DIG-corrected line fluxes prompts the use of a higher theoretical Ha/Hb ratio (3.03) than what is commonly used when recovering the E(B-V) via the Balmer decrement technique in massive star-forming galaxies. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: 58 pages, 46 figures. Paper accepted for pubblications in A&A. The catalogue will be available via the CDS or at the following link: http://dx.doi.org/10.11570/23.0006

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

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

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

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

Submitted 6 February, 2023; originally announced February 2023.

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

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

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

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

Abstract: We use new HCN(1-0) data from the ALMOND (ACA Large-sample Mapping Of Nearby galaxies in Dense gas) survey to trace the kpc-scale molecular gas density structure and CO(2-1) data from PHANGS-ALMA to trace the bulk molecular gas across 25 nearby, star-forming galaxies. At 2.1 kpc scale, we measure the density-sensitive HCN/CO line ratio and the SFR/HCN ratio to trace the star formation efficiency i… ▽ More We use new HCN(1-0) data from the ALMOND (ACA Large-sample Mapping Of Nearby galaxies in Dense gas) survey to trace the kpc-scale molecular gas density structure and CO(2-1) data from PHANGS-ALMA to trace the bulk molecular gas across 25 nearby, star-forming galaxies. At 2.1 kpc scale, we measure the density-sensitive HCN/CO line ratio and the SFR/HCN ratio to trace the star formation efficiency in the denser molecular medium. At 150 pc scale, we measure structural and dynamical properties of the molecular gas via CO(2-1) line emission, which is linked to the lower resolution data using an intensity-weighted averaging method. We find positive correlations (negative) of HCN/CO (SFR/HCN) with the surface density, the velocity dispersion and the internal turbulent pressure of the molecular gas. These observed correlations agree with expected trends from turbulent models of star formation, which consider a single free-fall time gravitational collapse. Our results show that the kpc-scale HCN/CO line ratio is a powerful tool to trace the 150 pc scale average density distribution of the molecular clouds. Lastly, we find systematic variations of the SFR/HCN ratio with cloud-scale molecular gas properties, which are incompatible with a universal star formation efficiency. Overall, these findings show that mean molecular gas density, molecular cloud properties and star formation are closely linked in a coherent way, and observations of density-sensitive molecular gas tracers are a useful tool to analyse these variations, linking molecular gas physics to stellar output across galaxy discs. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: 48 pages, 40 figures

Sub-kiloparsec empirical relations and excitation conditions of HCN and HCO+ J=3-2 in nearby star-forming galaxies

Authors: Axel Garcia-Rodriguez, Antonio Usero, Adam K. Leroy, Frank Bigiel, Maria Jesus Jimenez-Donaire, Daizhong Liu, Miguel Querejeta, Toshiki Saito, Eva Schinnerer, Ashley Barnes, Francesco Belfiore, Ivana Beslic, Yixian Cao, Melanie Chevance, Daniel A. Dale, Jakob S. den Brok, Cosima Eibensteiner, Santiago Garcia-Burillo, Simon C. O. Glover, Ralf S. Klessen, Jerome Pety, Johannes Puschnig, Erik Rosolowsky, Karin Sandstrom, Mattia C. Sormani , et al. (2 additional authors not shown)

Abstract: We present new HCN and HCO$^+$ ($J$=3-2) images of the nearby star-forming galaxies (SFGs) NGC 3351, NGC 3627, and NGC 4321. The observations, obtained with the Morita ALMA Compact Array, have a spatial resolution of $\sim$290-440 pc and resolve the inner $R_\textrm{gal} \lesssim$ 0.6-1 kpc of the targets, as well as the southern bar end of NGC 3627. We complement this data set with publicly avail… ▽ More We present new HCN and HCO$^+$ ($J$=3-2) images of the nearby star-forming galaxies (SFGs) NGC 3351, NGC 3627, and NGC 4321. The observations, obtained with the Morita ALMA Compact Array, have a spatial resolution of $\sim$290-440 pc and resolve the inner $R_\textrm{gal} \lesssim$ 0.6-1 kpc of the targets, as well as the southern bar end of NGC 3627. We complement this data set with publicly available images of lower excitation lines of HCN, HCO$^+$, and CO and analyse the behaviour of a representative set of line ratios: HCN(3-2)/HCN(1-0), HCN(3-2)/HCO$^+$(3-2), HCN(1-0)/CO(2-1), and HCN(3-2)/CO(2-1). Most of these ratios peak at the galaxy centres and decrease outwards. We compare the HCN and HCO$^+$ observations with a grid of one-phase, non-local thermodynamic equilibrium (non-LTE) radiative transfer models and find them compatible with models that predict subthermally excited and optically thick lines. We study the systematic variations of the line ratios across the targets as a function of the stellar surface density ($Σ_\textrm{star}$), the intensity-weighted CO(2-1) ($\langle I_\text{CO}\rangle$), and the star formation rate surface density ($Σ_\text{SFR}$). We find no apparent correlation with $Σ_\text{SFR}$, but positive correlations with the other two parameters, which are stronger in the case of $\langle I_\text{CO}\rangle$. The HCN/CO-$\langle I_\text{CO}\rangle$ relations show $\lesssim$0.3 dex galaxy-to-galaxy offsets, with HCN(3-2)/CO(2-1)-$\langle I_\text{CO}\rangle$ being $\sim$2 times steeper than HCN(1-0)/CO(2-1). In contrast, the HCN(3-2)/HCN(1-0)-$\langle I_\text{CO}\rangle$ relation exhibits a tighter alignment between galaxies. We conclude that the overall behaviour of the line ratios cannot be ascribed to variations in a single excitation parameter (e.g. density or temperature). △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: Accepted for publication in A&A. 14 pages, 8 figures

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

Broad-line region in NGC 4151 monitored by two decades of reverberation mapping campaigns. I. Evolution of structure and kinematics

Authors: Yong-Jie Chen, Dong-Wei Bao, Shuo Zhai, Feng-Na Fang, Chen Hu, Pu Du, Sen Yang, Zhu-Heng Yao, Yan-Rong Li, Michael S. Brotherton, Jacob N. McLane, T. E. Zastrocky, Kianna A. Olson, Edi Bon, Hua-Rui Bai, Yi-Xin Fu, Jun-Rong Liu, Yi-Lin Wang, Jaya Maithil, H. A. Kobulnicky, D. A. Dale, C. Adelman, M. J. Caradonna, Z. Carter, J. Favro , et al. (11 additional authors not shown)

Abstract: We report the results of long-term reverberation mapping (RM) campaigns of the nearby active galactic nuclei (AGN) NGC 4151, spanning from 1994 to 2022, based on archived observations of the FAST Spectrograph Publicly Archived Programs and our new observations with the 2.3m telescope at the Wyoming Infrared Observatory. We reduce and calibrate all the spectra in a consistent way, and derive light… ▽ More We report the results of long-term reverberation mapping (RM) campaigns of the nearby active galactic nuclei (AGN) NGC 4151, spanning from 1994 to 2022, based on archived observations of the FAST Spectrograph Publicly Archived Programs and our new observations with the 2.3m telescope at the Wyoming Infrared Observatory. We reduce and calibrate all the spectra in a consistent way, and derive light curves of the broad H$β$ line and 5100\,Å continuum. Continuum light curves are also constructed using public archival photometric data to increase sampling cadences. We subtract the host galaxy contamination using {\it HST} imaging to correct fluxes of the calibrated light curves. Utilizing the long-term archival photometric data, we complete the absolute flux-calibration of the AGN continuum. We find that the H$β$ time delays are correlated with the 5100\,Å luminosities as $τ_{\rm Hβ}\propto L_{5100}^{0.46\pm0.16}$. This is remarkably consistent with Bentz et al. (2013)'s global size-luminosity relationship of AGNs. Moreover, the data sets for five of the seasons allow us to obtain the velocity-resolved delays of the H$β$ line, showing diverse structures (outflows, inflows and disks). Combining our results with previous independent measurements, we find the measured dynamics of the H$β$ broad-line region (BLR) are possibly related to the long-term trend of the luminosity. There is also a possible additional $\sim$1.86 years time lag between the variation in BLR radius and luminosity. These results suggest that dynamical changes in the BLR may be driven by the effects of radiation pressure. △ Less

Submitted 15 January, 2023; originally announced January 2023.

Comments: Accepted for publication in MNRAS; comments welcome!

The PHANGS-MUSE Nebula Catalogue

Authors: B. Groves, K. Kreckel, F. Santoro, F. Belfiore, E. Zavodnik, E. Congiu, O. V. Egorov, E. Emsellem, K. Grasha, A. Leroy, F. Scheuermann, E. Schinnerer, E. J. Watkins, A. T. Barnes, F. Bigiel, D. A. Dale, S. C. O. Glover, I. Pessa, P. Sanchez-Blazquez, T. G. Williams

Abstract: Ionized nebulae provide critical insights into the conditions of the interstellar medium (ISM). Their bright emission lines enable the measurement of physical properties, such as the gas-phase metallicity, across galaxy disks and in distant galaxies. The PHANGS--MUSE survey has produced optical spectroscopic coverage of the central star-forming discs of 19 nearby main-sequence galaxies. Here, we u… ▽ More Ionized nebulae provide critical insights into the conditions of the interstellar medium (ISM). Their bright emission lines enable the measurement of physical properties, such as the gas-phase metallicity, across galaxy disks and in distant galaxies. The PHANGS--MUSE survey has produced optical spectroscopic coverage of the central star-forming discs of 19 nearby main-sequence galaxies. Here, we use the Hα morphology from this data to identify 30,790 distinct nebulae, finding thousands of nebulae per galaxy. For each nebula, we extract emission line fluxes and, using diagnostic line ratios, identify the dominant excitation mechanism. A total of 23,244 nebulae (75%) are classified as HII regions. The dust attenuation of every nebulae is characterised via the Balmer decrement and we use existing environmental masks to identify their large scale galactic environment (centre, bar, arm, interarm and disc). Using strong-line prescriptions, we measure the gas-phase oxygen abundances (metallicity) and ionization parameter for all HII regions. With this new catalogue, we measure the radial metallicity gradients and explore second order metallicity variations within each galaxy. By quantifying the global scatter in metallicity per galaxy, we find a weak negative correlation with global star formation rate and stronger negative correlation with global gas velocity dispersion (in both ionized and molecular gas). With this paper we release the full catalogue of strong line fluxes and derived properties, providing a rich database for a broad variety of ISM studies. △ Less

Submitted 10 January, 2023; originally announced January 2023.

Comments: 27 pages, accepted by MNRAS

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

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

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

Submitted 9 January, 2023; originally announced January 2023.

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

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

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

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

Submitted 2 January, 2023; originally announced January 2023.

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

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

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

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

Submitted 2 January, 2023; originally announced January 2023.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Submitted 21 December, 2022; originally announced December 2022.

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

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

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

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

Submitted 21 December, 2022; originally announced December 2022.

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