The PHANGS collaboration has been building a reference data set for the multiscale, multiphase study of star formation and the interstellar medium (ISM) 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 (∼5–50 pc). In Cycle 1, PHANGS is conducting an eight-band imaging survey from 2 to 21μ m of 19 nearby spiral galaxies. Optical integral field spectroscopy, CO(2–1) mapping, 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 ∼10,000 star clusters, MUSE spectroscopic mapping of ∼20,000 Hiiregions, and ∼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.

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 and 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–HSTB -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 Hiiregions, young star clusters, and associations within the DFN. We examine the dependence of these quantities on 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 cospatial 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 the star formation in Hiiregions and 60% of the mass in star clusters younger than 5 Myr are contained within the DFN. However, the DFN at this scale is anticorrelated with looser associations of stars younger than 5 Myr identified using PHANGS–HST near-UV imaging. We discuss the impact of these findings on studies of star formation and the interstellar medium, and the broad range of new investigations enabled by multiscale maps of the DFN.

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 Hiiregions. 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 Hiiregions 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 Hiiregions. This relation becomes steeper for more luminous Hiiregions. The metallicity of Hiiregions 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 Hiiregions—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.

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 Mid-Infrared Instrument observations of IC 5332, NGC 628, NGC 1365, and NGC 7496 from the Physics at High Angular resolution in Nearby GalaxieS-JWST Treasury. Using scaling relationships calibrated in Leroy et al., scaled F1130W provides 10–40 pc resolution and 3 σsensitivity of Σ_gas ∼ 2M _⊙ pc⁻². We characterize the surface densities of structures seen at <7M _⊙ pc⁻² in our targets, where we expect the gas to be Hi -dominated. We highlight the existence of filaments, interarm emission, and holes in the diffuse ISM at these low surface densities. Below ∼10M _⊙ pc⁻² for NGC 628, NGC 1365, and NGC 7496 the gas distribution shows a "Swiss cheese"-like topology due to holes and bubbles pervading the relatively smooth distribution of the diffuse ISM. Comparing to recent galaxy simulations, we observe similar topology for the low-surface-density gas, though with notable variations between simulations with different setups and resolution. Such a comparison of high-resolution, low-surface-density gas with simulations is not possible with existing atomic and molecular gas maps, highlighting the unique power of JWST maps of PAH emission.

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 indexp= −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.

We present new 0.3–21μ m photometry of SN 2021aefx in the spiral galaxy NGC 1566 at +357 days afterB -band maximum, including the first detection of any Type Ia supernova (SN Ia) at >15μ m. These observations follow earlier JWST observations of SN 2021aefx at +255 days after the time of maximum brightness, allowing us to probe the temporal evolution of the emission properties. We measure the fraction of flux emerging at different wavelengths and its temporal evolution. Additionally, the integrated 0.3–14μ m decay rate of Δ m _0.3–14 = 1.35 ± 0.05 mag/100 days is higher than the decline rate from the radioactive decay of ⁵⁶Co of ∼1.2 mag/100 days. The most plausible explanation for this discrepancy is that flux is shifting to >14μ m, and future JWST observations of SNe Ia will be able to directly test this hypothesis. However, models predicting nonradiative energy loss cannot be excluded with the present data.