Spatially resolved maps of gas-phase metallicity provide key constraints on the chemical enrichment and mixing processes that drive galaxy evolution, but measurements based only on strong lines remain highly uncertain and dependent on emission-line coverage. In this work, we present a joint analysis of SITELLE observations, covering the [O II]λλ3726,3729 doublet, with PHANGS-MUSE spectroscopy covering 4800-9300 Angstroms, including Hβ, [O III]λ4959,5007, [N II]λ6584, Hα, [S II]λλ6716,6731, and [S III]λ9069, within five nearby spiral galaxies. By combining these data, we construct a homogeneous catalog of emission-line fluxes for 604 ionized nebulae, 556 of which are classified as H II regions. This enables a comparison of eight widely used strong-line metallicity calibrations, five new strong-line calibrations, and an investigation of ionization-parameter diagnostics. We recover known systematic offsets among calibrations, but also find that many exhibit very low scatter, less than 0.03-0.04 dex, in radial metallicity gradients. We find that [S III]/[S II] exhibits minimal secondary dependence on metallicity or extinction, suggesting that it may be a more robust tracer of ionization parameter than [O III]/[O II]. No significant outliers are identified in O/H or N/O within the sampled regions, indicating internally consistent abundance trends across the inner disks probed by our data. We provide a publicly available catalog of all measured emission-line fluxes, designed to support future investigations, including temperature modeling and strong-line abundance calibrations.

We present PHANGS-H alpha, a narrowband imaging survey that maps H alpha emission over a sample of 65 nearby massive star-forming galaxies. The data were obtained using the MPG-ESO 2.2-meter telescope at La Silla and the du Pont 2.5-meter telescope at Las Campanas Observatory, within the framework of the multi-wavelength cloud-scale (50-100 pc) resolution mapping of molecular gas and star formation conducted by the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) collaboration. PHANGS-H alpha complements the published PHANGS-ALMA, PHANGS-MUSE, PHANGS-HST, and PHANGS-JWST surveys, providing an anchor point for the photometric and astrometric calibration of these datasets, as well as samples of H II regions and star formation rate maps for the bulk of the PHANGS sample. We describe the observations, data processing, and calibration of the PHANGS-H alpha dataset, as well as the procedures used to derive emission-line fluxes from narrowband imaging. A subset of galaxies with available spectroscopic Ha mapping from the PHANGS-MUSE survey serves as the basis of a detailed comparison with the narrowband photometry presented here. This comparison informs a set of best practices for the processing of narrowband H alpha imaging, which we subsequently apply to the full dataset.

Newly formed stars profoundly affect their environment by depositing energy and momentum into the surrounding gas. However, only a fraction of the stellar feedback is retained in the cloud, and observational constraints are needed to improve our understanding of this process. In a sample of 19 nearby galaxies, we matched H II regions from PHANGS–MUSE to their ionizing stellar source from PHANGS–HST and measured the percentage of ionizing radiation that leaks into the surrounding diffuse ionised gas (DIG). Based on a catalogue in which each H II region is powered by a single young and massive stellar association, we measure a photon escape fraction of fesc = 82−24+12 %. We obtain comparable results when using different procedures to match the ionised gas to its source. All samples in our study contain a substantial fraction of objects (up to 20%), in which the stellar source is insufficient to produce the H α flux observed from the nebula. Many of these cases are probably related to uncertain age estimates, but we also find numerous regions for which a significant fraction of the ionising photon budget is contributed by stars that reside outside the boundaries of the H II region. This finding motivates the use of an alternative galaxy-wide approach in which we include all H II regions and stellar sources, not just those that show a clear overlap. When we sum the ionisation budget over entire galaxies, we measure slightly lower, but consistent values.

We present PHANGS-H α , a narrowband imaging survey that maps H α emission over a sample of 65 nearby massive star-forming galaxies. The data were obtained using the MPG-ESO 2.2-meter telescope at La Silla and the du Pont 2.5-meter telescope at Las Campanas Observatory, within the framework of the multi-wavelength cloud-scale (50–100 pc) resolution mapping of molecular gas and star formation conducted by the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) collaboration. PHANGS-H α complements the published PHANGS-ALMA, PHANGS-MUSE, PHANGS-HST, and PHANGS-JWST surveys, providing an anchor point for the photometric and astrometric calibration of these datasets, as well as samples of H II regions and star formation rate maps for the bulk of the PHANGS sample. We describe the observations, data processing, and calibration of the PHANGS-H α dataset, as well as the procedures used to derive emission-line fluxes from narrowband imaging. A subset of galaxies with available spectroscopic Ha mapping from the PHANGS-MUSE survey serves as the basis of a detailed comparison with the narrowband photometry presented here. This comparison informs a set of best practices for the processing of narrowband H α imaging, which we subsequently apply to the full dataset.

Newly formed stars profoundly affect their environment by depositing energy and momentum into the surrounding gas. However, only a fraction of the stellar feedback is retained in the cloud, and observational constraints are needed to improve our understanding of this process. In a sample of 19 nearby galaxies, we matched H II regions from PHANGS-MUSE to their ionizing stellar source from PHANGS-HST and measured the percentage of ionizing radiation that leaks into the surrounding diffuse ionised gas (DIG). Based on a catalogue in which each H II region is powered by a single young and massive stellar association, we measure a photon escape fraction of f(esc) = 82(-24)(+12) %. We obtain comparable results when using different procedures to match the ionised gas to its source. All samples in our study contain a substantial fraction of objects (up to 20%), in which the stellar source is insufficient to produce the H alpha flux observed from the nebula. Many of these cases are probably related to uncertain age estimates, but we also find numerous regions for which a significant fraction of the ionising photon budget is contributed by stars that reside outside the boundaries of the H II region. This finding motivates the use of an alternative galaxy-wide approach in which we include all H II regions and stellar sources, not just those that show a clear overlap. When we sum the ionisation budget over entire galaxies, we measure slightly lower, but consistent values.

We present PHANGS-Hα, a narrowband imaging survey that maps Hα emission over a sample of 65 nearby massive star-forming galaxies. The data were obtained using the MPG-ESO 2.2-meter telescope at La Silla and the du Pont 2.5-meter telescope at Las Campanas Observatory, within the framework of the multi-wavelength cloud-scale (50–100 pc) resolution mapping of molecular gas and star formation conducted by the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) collaboration. PHANGS-Hα complements the published PHANGS-ALMA, PHANGS-MUSE, PHANGS-HST, and PHANGS-JWST surveys, providing an anchor point for the photometric and astrometric calibration of these datasets, as well as samples of H II regions and star formation rate maps for the bulk of the PHANGS sample. We describe the observations, data processing, and calibration of the PHANGS-Hα dataset, as well as the procedures used to derive emission-line fluxes from narrowband imaging. A subset of galaxies with available spectroscopic Ha mapping from the PHANGS-MUSE survey serves as the basis of a detailed comparison with the narrowband photometry presented here. This comparison informs a set of best practices for the processing of narrowband Hα imaging, which we subsequently apply to the full dataset.

JWST observations, when combined with Hubble Space Telescope (HST) data, promise to improve age estimates of star clusters in nearby spiral galaxies. However, feedback from young cluster stars pushes out the natal gas and dust, making cluster formation and evolution a challenge to model. Here, we use JWST+ HST observations of the nearby spiral galaxy NGC 628 to produce spectral energy distribution (SED) templates of compact star clusters spanning 275 nm through 21 μ m. These preliminary SEDs capture the cluster stars and associated gas and dust within radii of ≈0 . ″ 12–0 . ″ 67 (corresponding to ≈6–33 pc at the distance of NGC 628). One important finding is that the SEDs of 1, 2, 3, and 4 Myr clusters can be differentiated in the infrared. Another is that, in 80%–90% of the cases we study, the polycyclic aromatic hydrocarbon (PAH) and H α emission track one another, with the dust responsible for the 3.3 μ m PAH emission largely removed by 4 Myr, consistent with pre-supernova stellar feedback acting quickly on the surrounding gas and dust. Nearly embedded cluster candidates have infrared SEDs that are quite similar to optically visible 1–3 Myr clusters. In nearly all cases, we find there is a young star cluster within a few tenths of an arcsec (10–30 pc) of the nearly embedded cluster, suggesting the formation of the cluster was triggered by its presence. The resulting age estimates from the empirical templates are compatible both with dynamical estimates based on CO superbubble expansion velocities, as well as the TODDLERS models, which track spherical evolution of homogeneous gas clouds around young stellar clusters.

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

In this work, we present a new catalogue of >40 000 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 main 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. In fact, 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), which is 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 in 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 diffuse-ionised-gas-corrected line fluxes prompts the use of a higher theoretical H alpha/H beta ratio (3.03) than what is commonly used when recovering the E(B - V) via the Balmer decrement technique in massive star-forming galaxies.

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.

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{\alpha} 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.

Connecting the gas in H II regions to the underlying source of the ionizing radiation can help us constrain the physical processes of stellar feedback and how H II regions evolve over time. With PHANGS-MUSE, we detect nearly 24 000 H II regions across 19 galaxies and measure the physical properties of the ionized gas (e.g. metallicity, ionization parameter, and density). We use catalogues of multiscale stellar associations from PHANGS-HST to obtain constraints on the age of the ionizing sources. We construct a matched catalogue of 4177 H II regions that are clearly linked to a single ionizing association. A weak anticorrelation is observed between the association ages and the H a equi v alent width EW (H a), the H a/ 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 a) and log q show the most consistent trends and appear to be most reliable tracers for the age of an H II region.

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 PHANGSx2013MUSE 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 PHANGSx2013HST 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.

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.

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.

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

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.

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

Aims: Temperature uncertainties plague our understanding of abundance variations within the ISM. Using the PHANGS-MUSE large program, we develop and apply a new technique to model the strong emission lines arising from HII regions in 19 nearby spiral galaxies at 50 pc resolution and infer electron temperatures for the nebulae. Methods: Due to the charge-exchange coupling of the ionization fraction of the atomic oxygen to that of hydrogen, the emissivity of the observed [OI]6300/Ha line ratio can be modeled as a function of gas phase oxygen abundance (O/H), ionization fraction (f_ion) and electron temperature (T_e). We measure (O/H) using a strong line metallicity calibration, and identify a correlation between f_ion and [SIII]9069/[SII]6716,6730, tracing ionization parameter variations. Results: We solve for T_e, and test the method by reproducing direct measurements of T_e([NII]5755) based on auroral line detections to within 600 K. We apply this charge-exchange method of calculating T_e to 4,129 HII regions across 19 PHANGS-MUSE galaxies. We uncover radial temperature gradients, increased homogeneity on small scales, and azimuthal temperature variations in the disks that correspond to established abundance patterns. This new technique for measuring electron temperatures leverages the growing availability of optical integral field unit spectroscopic maps across galaxy samples, increasing the statistics available compared to direct auroral line detections.