Context. Estimating properties of star clusters from unresolved broadband photometry is a challenging problem that is classically tackled using spectral energy distribution (SED) fitting methods that are based on simple stellar population models. However, grid-based methods suffer from computational limitations. Because of their exponential scaling, they can become intractable when the number of inference parameters grows. In addition, nuisance parameters in the model can make the computation of the likelihood function intractable. These limitations can be overcome by modern generative deep learning methods that offer flexible and powerful tools for modeling high-dimensional posterior distributions and fast inference from learned data. Aims. We present a normalizing flow approach for the inference of cluster age, mass, and reddening parameters from Hubble Space Telescope broadband photometry. In particular, we explore our network's behavior when dealing with an inference problem that has been analyzed in previous works. Methods. We used the SED modeling code CIGALE to create a dataset of synthetic photometric observations for 5 & times; 106 mock star clusters. Subsequently, this dataset was used to train a coupling-based flow in the form of a conditional invertible neural network to predict posterior probability distributions for cluster age, mass, and reddening from photometric observations. Results. We predicted cluster parameters for the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) Data Release 3 catalog. To evaluate the capabilities of the network, we compared our results to the publicly available PHANGS estimates and found that the estimates agree reasonably well. Conclusions. We demonstrate that normalizing flow methods can be a viable tool for the inference of cluster parameters, and argue that this approach is especially useful when nuisance parameters make the computation of the likelihood intractable and in scenarios that require efficient density estimation.

Light reprocessed by dust grains emitting in the infrared enables the study of the physics at play in dusty embedded regions, where ultraviolet and optical wavelengths are attenuated. Infrared telescopes such as JWST have made it possible to study the earliest feedback phases, when stars are shielded by cocoons of gas and dust. Comprehending this phase is crucial for unravelling the effects of feedback from young stars that leads to their emergence and the dispersal of their host molecular clouds. Here we show that the transition from the embedded to the exposed phase of star formation is short ($< 4$ Myr) and sometimes almost absent ($< 1$ Myr) across a sample of 37 nearby star-forming galaxies covering a wide range of morphologies, from massive barred spirals to irregular dwarfs. The short duration of the dust-clearing timescales suggests a predominant role of pre-supernova feedback mechanisms in revealing newborn stars, confirming previous results on smaller samples and allowing, for the first time, a statistical analysis of their dependencies. We find that the timescales associated with mid-infrared emission at 21 μm, tracing a dust-embedded feedback phase, are controlled by a complex interplay between giant molecular cloud properties (masses and velocity dispersions) and galaxy morphology. We report relatively longer durations of the embedded phase of star formation in barred spiral galaxies, while this phase is significantly reduced in low-mass irregular dwarf galaxies. We discuss tentative trends with gas-phase metallicity, which may favor faster cloud dispersal at low metallicities.

Red supergiant (RSG) stars represent the final evolutionary phase of the majority of massive stars and hold a unique role in testing the physics of stellar models. Eighty-eight RSGs in the Small Magellanic Cloud (SMC) were recently found to have an ultraviolet (UV) excess, which was attributed to a B-type companion. We present follow-up Hubble Space Telescope (HST) UV (1700-3000 & Aring;) spectroscopy for 16 of these stars to investigate the nature of the UV excess and confirm the presence of a hot companion. In all cases, we are able to confirm that the companion is a main-sequence B-type star based on the near-UV continuum. We determine effective temperatures, radii, and luminosities from fitting the UV continuum with TLUSTY models and find stellar parameters in the expected range for SMC B-type stars. We display these results on a Hertzsprung-Russell diagram and assess the previously determined stellar parameters using UV photometry alone. From this comparison, we conclude that UV photometric surveys are vital to identify such companions and that UV spectroscopy is similarly vital to characterise the hot companions. From a comparison with IUE spectra of 32 Cyg, a well-known RSG binary system in the Galaxy where the companion is embedded in the wind of the RSG, four targets display evidence of also being embedded in such a wind. The ages of seven targets, determined via the stellar parameters of the hot companions, are found to be in tension with the ages determined for the RSG. An explanation for the discrepancy could be unaccounted for binary mass-transfer or red straggler stars.

This paper is the second in a series presenting the catalogs and properties of the largest sample to date of ∼100,000 star clusters and compact associations, in 38 spiral galaxies observed by the PHANGS-HST Treasury survey. Here, we present spectral energy distribution (SED) fitting techniques used to compute the age, mass, and reddening for each object. Our decision-tree-based strategy incorporates categorical priors on model age, reddening, and metallicity determined from additional observed parameters: localized H α emission, source morphology, and demographic-specific locations in the UBVI color–color diagram. This approach is implemented to mitigate model degeneracies, particularly between young dusty clusters and old clusters with minimal dust, which can have identical optical colors. Results based on H α narrowband imaging from the ground and from Hubble Space Telescope are intercompared, and contrasted with previous SED-fitting efforts. The fraction of the population that is subject to such priors is ∼14%, and of this subset, ∼63% of old globular clusters (GCs) have ages that change by a factor of 10 or more relative to unconstrained fits with single metallicity ( Z _⊙ ) simple stellar population models. The demographics of the population are examined through age–mass and age–reddening diagrams (for individual galaxies as well as aggregated over the sample), and the GC mass function. We demonstrate relationships between cluster age–mass diagrams and properties of parent galaxies (galaxy morphology and location relative to the galaxy main sequence). We outline continuing efforts to improve the inference of physical properties, including the incorporation of JWST infrared photometry and updated synthesis models.

We present a multiwavelength catalog of mid-infrared-selected compact sources in 19 nearby galaxies, combining JWST NIRCam/MIRI, HST UV-optical broadband, Hα narrow-band, and ALMA CO observations. We detect 24,945 compact sources at 21 μm and 55,581 at 10 μm. Artificial star tests show 50

The large number of star clusters in nearby galaxies permits us to statistically test the predictions of stellar, dust, and gas models. Using Hubble Space Telescope (HST) broadband plus H α imaging combined with JWST near-infrared imaging, we use a total of 10 filters spanning near-ultraviolet through near-infrared wavelengths to model key physical parameters, including age, mass, and reddening, of 6130 star clusters in 16 nearby spiral galaxies from the Physics at High Angular resolution in Nearby GalaxieS sample, focusing on their ages, masses, and reddenings. We find that HST/H α and JWST/NIRCam 2–3.6 μ m photometry significantly improves our ability to disentangle the age–reddening degeneracy between young, gas- and dust-rich clusters and older, dustless clusters. The near-infrared data provide strong constraints on hot continuum dust and small polycyclic aromatic hydrocarbon emission for populations where gas and dust are present. These hot dust constraints demonstrate that Bruzual & Charlot stellar population models do not align with the observed near-ultraviolet-optical-near-infrared spectral energy distributions of star clusters in the first 10 Myr. We note that for old and low-metallicity globular clusters, the inclusion of narrowband H α and/or broadband near-infrared data does not improve the determination of age and reddening parameters, due to the lack of stars capable of heating dust in the near-infrared regime.

The Multiphase Astrophysics to Unveil the Virgo Environment (MAUVE) project is a multi-facility programme exploring how dense environments transform galaxies. Combining a VLT/MUSE P110 Large Programme and ALMA observations of 40 late-type Virgo Cluster galaxies, MAUVE resolves star formation, kinematics, and chemical enrichment within their molecular gas discs. A key goal is to track the evolution of cold gas that survives in the inner regions of satellites after entering the cluster, and how it evolves across different infall stages. With its high spatial resolution – probing down to the physical scales of giant molecular cloud complexes – and multiphase synergy, MAUVE aims to offer a time-resolved view of environmental quenching and set a new benchmark for cluster galaxy studies.

We present a comparison of observed polycyclic aromatic hydrocarbon (PAH) feature ratios in 19 nearby galaxies with a grid of theoretical expectations for near- and mid-infrared dust emission. The PAH feature ratios are drawn from Cycle 1 JWST observations and are measured for 7224 stellar clusters and 29,176 stellar associations for which we have robust ages and mass estimates from Hubble Space Telescope five-band photometry. Though there are galaxy-to-galaxy variations, the observed PAH feature ratios largely agree with the theoretical models, particularly those that are skewed toward more ionized and larger PAH size distributions. For each galaxy we also extract PAH feature ratios for 200 pc wide circular regions in the diffuse interstellar medium, which serve as a noncluster/association control sample. Compared to what we find for stellar clusters and associations, the 3.3 μ m/7.7 μ m and 3.3 μ m/11.3 μ m ratios from the diffuse interstellar medium are ∼0.10–0.15 dex smaller. When the observed PAH feature ratios are compared to the radiation field hardness as probed by the [O iii ]/H β ratio, we find anticorrelations for nearly all galaxies in the sample. These results together suggest that the PAH feature ratios are driven by the shape and intensity of the radiation field and that the smallest PAHs—observed via JWST F335M imaging—are increasingly “processed” or destroyed in regions with the most intense and hard radiation fields.

The PHANGS project is assembling a comprehensive, multiwavelength data set of nearby (∼5–20 Mpc), massive star-forming galaxies to enable multiphase, multiscale investigations into the processes that drive star formation and galaxy evolution. To date, large survey programs have provided molecular gas (CO) cubes with the Atacama Large Millimeter/submillimeter Array, optical integral field unit (IFU) spectroscopy with the Very Large Telescope/Multi-Unit Spectroscopic Explorer (MUSE), high-resolution near-ultraviolet–optical imaging in five broadband filters with Hubble Space Telescope (HST), and infrared imaging in NIRCAM+MIRI filters with JWST. Here we present PHANGS-HST-H α , which has obtained high-resolution (∼2–10 pc), narrowband imaging in the F658N or F657N filters with the HST/WFC3 camera of the warm ionized gas in the first 19 nearby galaxies observed in common by all four of the PHANGS large programs. We summarize our data reduction process, with a detailed discussion of the production of flux-calibrated, Milky Way extinction-corrected, continuum-subtracted H α maps. PHANGS-MUSE IFU spectroscopy data are used to background-subtract the HST-H α maps and to determine the [N ii ] correction factors for each galaxy. We describe our public data products (the data released as part of this work include the reduced drizzled narrowband images and the flux-calibrated, continuum-subtracted H α maps for each galaxy; these images are available for download via MAST at https://archive.stsci.edu/hlsp/phangs.html , as well as at the Canadian Astronomy Data Centre as part of the PHANGS archive at https://www.canfar.net/storage/vault/list/phangs/RELEASES ) and highlight a few key science cases enabled by the PHANGS-HST-H α observations.

Red supergiants (RSGs) are important for our understanding of supernova progenitors, stellar populations, stellar evolution, mass loss and dust production. Extragalactic surveys of RSGs have a long history in the Local Group, but few studies exist beyond that due to the limited resolution and sensitivity of ground-based and previous space-based infrared observatories. Here we demonstrate the combined power of HST and JWST to push systematic searches of RSGs out to ∼20 Mpc. We introduce a catalog of 97057 RSGs – the largest single-survey release of RSGs – with masses ≳10 M_⊙ in 19 galaxies from the PHANGS HST+JWST Treasury program. We use HST F814W and JWST F200W photometry to select stars as RSGs based on predicted colors and magnitudes from PARSEC isochrones. The spatial distribution of our recovered RSGs follow the familiar pattern of mostly being concentrated in active star-forming regions such as spiral arms and central starburst rings. The RSG number density on kpc-scales is strongly correlated (r_s∼0.82) with local star-formation rate density (Σ_SFR) traced by extinction-corrected far-ultraviolet (FUV) from GALEX+WISE, and weakly correlated (r_s∼0.57) with the total stellar mass density (Σ_*), traced by near-infrared emission from WISE+Spitzer. The number of RSGs per mass of stellar populations with ages 6-30 Myr (the likely age range of RSGs >10 M_⊙) is ∼1 per 10^3.77±0.27 M_⊙, assuming constant star-formation rates from FUV+W4. Our sample will be a useful resource for tracking progenitors and feedback sites of future supernovae in PHANGS, age-dating stellar populations, and more.

Polycyclic aromatic hydrocarbons (PAHs) are widespread in the interstellar medium (ISM) of near solar metallicity galaxies, where they play a critical role in ISM heating, cooling, and reprocessing stellar radiation. The PAH fraction, the abundance of PAHs relative to total dust mass, is a key parameter in ISM physics. Using JWST and MUSE observations of 42 galaxies from the PHANGS survey, we analyzed the PAH fraction in over 17 000 H II regions spanning a gas-phase oxygen abundance of 12 + log(O/H) = 8.0-8.8 (Z similar to 0.2-1.3 Z(circle dot)), and similar to 400 isolated supernova remnants (SNRs). We find a significantly lower PAH fraction toward H II regions compared to a reference sample of diffuse ISM areas at matched metallicity. At 12 + log(O/H) > 8.2, the PAH fraction toward H II regions is strongly anti-correlated with the local ionization parameter, suggesting that PAH destruction is correlated with ionized gas and/or hydrogen-ionizing UV radiation. At lower metallicities, the PAH fraction declines steeply in H II regions and in the diffuse ISM, likely reflecting less efficient PAH formation in metal-poor environments. Carefully isolating dust emission from the vicinity of optically identified supernova remnants, we see evidence of selective PAH destruction from measurements of lower PAH fractions, which is, however, indistinguishable at similar to 50 pc scales. Overall, our results point to ionizing radiation as the dominant agent of PAH destruction within H II regions; metallicity plays a key role in their global abundance in galaxies.

We present the variations in far-ultraviolet (FUV) and H α star formation rates (SFR), SFR _UV and SFR _H _α , respectively, at subkiloparsec scales in 11 galaxies as part of the Deciphering the Interplay between the Interstellar Medium, Stars, and the Circumgalactic medium survey. Using archival GALEX FUV imagery and H α +[N ii ] narrowband images obtained with the Vatican Advanced Technology Telescope, we detect a total of 1335 (FUV-selected) and 1474 (H α -selected) regions of recent high-mass star formation, respectively. We find the H α -to-FUV SFR ratios tend to be lower primarily for FUV-selected regions, where SFR _H _α generally underestimates the SFR by an average factor of 2–3, for SFR < 10 ^−4 M _⊙ yr ^−1 . In contrast, the SFRs are generally observed to be consistent for H α -selected regions. This discrepancy arises from morphological differences between the two indicators. Extended FUV morphologies and larger areas covered by FUV-only regions, along with decreasing overlap between FUV clumps and compact H ii regions with R / R _25 suggest that stochastic sampling of the initial mass function may be more pronounced in the outer regions of galaxies. Our observed H α -to-FUV SFR ratios are also consistent with stochastic star formation model predictions. However, using larger apertures that include diffuse FUV emission results in an offset of 1 dex between SFR _H _α and SFR _UV , suggesting that the observed low H α -to-FUV SFR ratios in galaxies are likely caused by diffuse FUV emission, which can contribute ∼60%–90% to the total FUV flux.

We characterize two candidate cool galactic outflows in two relatively low-mass, highly inclined Virgo cluster galaxies: NGC 4424 and NGC 4694. Previous analyses of observations using the Atacama Large Millimeter/submillimeter Array carbon monoxide (CO) line emission maps did not classify these sources as cool outflow hosts. Using new high-sensitivity, high-spatial-resolution, JWST mid-infrared photometry in the polycyclic aromatic hydrocarbon (PAH)-tracing F770W band, we identify extended structures present off of the stellar disk. The identified structures are bright in the MIRI F770W and F2100W bands, suggesting they include PAHs as well as other dust grains. As PAHs have been shown to be destroyed in hot, ionized gas, these structures are likely to be outflows of cool (T <= 104 K) gas. This work represents an exciting possibility for using mid-infrared observations to identify and measure outflows in lower-mass, lower star formation galaxies.

The earliest stages of star and cluster formation are hidden within dense cocoons of gas and dust, limiting their detection at optical wavelengths. With the unprecedented infrared capabilities of JWST, we can now observe dust-enshrouded star formation with ∼10 pc resolution out to ∼20 Mpc. Early findings from PHANGS-JWST suggest that 3.3 μ m polycyclic aromatic hydrocarbon (PAH) emission can identify star clusters in their dust-embedded phases. Here, we extend this analysis to 19 galaxies from the PHANGS-JWST Cycle 1 Treasury survey, providing the first characterization of compact sources exhibiting 3.3 μ m PAH emission across a diverse sample of nearby star-forming galaxies. We establish a selection criteria based on a median color threshold of F300M − F335M = 0.67 at F335M = 20 and identify 1816 sources. These sources are predominantly located in dust lanes, spiral arms, rings, and galaxy centers, with ∼87% showing concentration indices (CIs) similar to optically detected star clusters. Comparison with the PHANGS-HST catalogs suggests that PAH emission fades within ∼3 Myr. The H α equivalent width of PAH emitters is 1–2.8 times higher than that of young PHANGS-HST clusters, providing evidence that PAH emitters are on average younger. Analysis of the bright portions of luminosity functions (which should not suffer from incompleteness) shows that young dusty clusters may increase the number of optically visible ≤3 Myr old clusters in PHANGS-HST by a factor between ∼1.8× and 8.5×.

Combining Atacama Large Millimeter/submillimeter Array CO(2–1) mapping and JWST near- and mid-infrared imaging, we characterize the relationship between CO(2–1) and polycyclic aromatic hydrocarbon (PAH) emission at ≈100 pc resolution in 70 nearby star-forming galaxies. Leveraging a new Cycle 2 JWST Treasury program targeting nearby galaxies, we expand the sample size by more than an order of magnitude compared to previous ≈100 pc resolution CO–PAH comparisons. In regions of galaxies where most of the gas is likely to be molecular, we find strong correlations between CO(2–1) and 3.3 μ m, 7.7 μ m, and 11.3 μ m PAH emission, estimated from JWST’s F335M, F770W, and F1130W filters. We derive power-law relations between CO(2–1) and PAH emission, with indices in the range 0.8–1.3, implying relatively weak variations in the observed CO-to-PAH ratios across our sample. We find that CO-to-PAH ratios and scaling relationships near H ii regions are similar to those in diffuse sight lines. The main difference between the two types of regions is that sight lines near H ii regions show higher intensities in all tracers. Galaxy centers show higher overall intensities and enhanced CO-to-PAH ratios compared to galaxy disks. Individual galaxies show 0.19 dex scatter in the normalization of CO at fixed I _PAH , and this normalization anticorrelates with specific star formation rate and correlates with stellar mass. We provide a prescription that accounts for galaxy-to-galaxy variations, representing our best current empirical predictor to estimate CO(2–1) intensity from PAH emission, allowing one to take advantage of JWST’s excellent sensitivity and resolution to trace cold gas.

Recent JWST mid-infrared (mid-IR) images, tracing polycyclic aromatic hydrocarbons (PAHs) and dust continuum emission, provide detailed views of the interstellar medium (ISM) in nearby galaxies. Leveraging PHANGS-JWST Cycle 1 and PHANGS-MUSE data, we measure the PAH and dust continuum emission lifetimes of gas clouds across 17 nearby star-forming galaxies by analyzing the relative spatial distributions of mid-IR (7.7–11.3 μ m) and H α emission at various scales. We find that the mid-IR emitting timescale of gas clouds in galactic disks (excluding centers) ranges from 10–30 Myr. After star formation is detected in H α , mid-IR emission persists for 3–7 Myr during the stellar feedback phase, covering 70%–80% of the H α emission. This significant overlap is due to intense radiation from star-forming regions, illuminating the surrounding PAHs and dust grains. In most galaxies, the mid-IR time-scale closely matches the molecular cloud lifetime measured with CO. Although mid-IR emission is complex, as influenced by ISM distribution, radiation, and abundances of dust and PAHs, the similarity between the two timescales suggests that once gas clouds form with compact mid-IR emission, they quickly provide sufficient shielding for stable CO formation. This is likely due to our focus on molecular gas-rich regions of galaxies with near-solar metallicity. Finally, we find that the mid-IR emitting timescale is longer in galaxies with well-defined H ii regions and less structured backgrounds, allowing photons to more efficiently heat the ambient ISM surrounding the H ii regions, rather than contributing to diffuse emission. This suggests that the shape of the ISM also influences mid-IR emission.

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 exquisite angular resolution and sensitivity of JWST are 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 ∼5–50 pc. 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%) reduces the noise by a factor of ≈3. We also present the first public release (DR1.1.0) of the pjpipe processed eight-band 2–21 μ m imaging for all 19 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.

We compare the molecular cloud properties in subgalactic 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, line widths, pressures, and virial parameters in the arms of NGC 1313 and the 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 environment, allowing more clouds with the necessary conditions to emerge, 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 the clouds in the arms of NGC 1313, but have densities and pressures more similar to those in 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 density between the regions, suggesting that the differences in surface densities frequently seen between arm and interarm regions in lower-resolution studies are indicative of the sparsity of molecular clouds, rather than differences in their true surface density.