Dr Thomas Williams

Abstract:

We present 50–100 pc-resolution JWST/MIRI and NIRCam measurements of mid-infrared (mid-IR) color variations in the diffuse interstellar medium (ISM) of 71 nearby star-forming galaxies from the PHANGS-JWST survey. Mid-IR emission traces the dust column density, intensity (U) and hardness of the interstellar radiation field, and the physical state (charge, size) and abundance of polycyclic aromatic hydrocarbons (PAHs). Mid-IR colors that trace PAH band ratios remain fairly constant in the diffuse ISM of star-forming disks. However, they show stark variations in extreme environments: highly star-forming central molecular zones (CMZs) and star formation deserts/quiescent bulges. In CMZs, PAH-to-continuum (3.3/21, 7.7/21, and 11.3/21 μm) and the 10/21 μm continuum colors are 0.2–0.4 dex lower than in normal disks. We attribute this to higher U based on the far-infrared dust colors and the high 21 μm/ΣMol, which we suggest to be a good tracer of U outside star-forming regions. Meanwhile, star formation deserts show low 7.7 μm PAH emission, resulting in low 7.7/21 μm and 7.7/11.3 μm, while all other mid-IR colors remain typical. This suggests the presence of more neutral PAHs in star formation deserts, where low 7.7 μm likely reflects ISM conditions similar to early-type and elliptical galaxies. All environments form part of a continuous trend in 7.7/11.3 μm versus specific star formation rate.

Abstract:

We use JWST NIRCam and MIRI imaging at 2, 4, 10, and 21 μm to study young, dusty compact sources in four nearby galaxies at distances of ∼1–5 Mpc (M33, NGC 300, NGC 7793, and NGC 5068). This work bridges well-characterized massive young stellar objects (MYSOs) in the Magellanic Clouds from the Spitzer Space Telescope SAGE survey to new studies of embedded clusters in more distant galaxies with JWST. Guided by the SAGE-LMC catalog, we define JWST color–magnitude selection criteria (F1000W versus F1000W − F2100W) and test them using resolution-degradation experiments. We identify 216, 32, 80, and 139 dusty young objects in the four galaxies, respectively. The selected population spans sources from systems dominated by a single MYSO to compact marginally resolved sources hosting multiple MYSOs. The color selection remains stable across 1–5 Mpc, and the 10 μm luminosity function retains a slope of α ∼ −2. However, blending and surface-brightness dilution remove fainter sources, leading to incompleteness of up to ∼50% at 5.2 Mpc and biasing the sample toward brighter objects (F1000W < 19 mag). The sample spans approximate stellar masses of ∼10–2 × 10⁵ M⊙. Spatial resolution affects the interpretation of mid-IR emission: clustering increases the fraction of emission attributed to compact sources in active regions, while blending into diffuse emission dominates in quiescent environments. Comparisons with polycyclic aromatic hydrocarbon (PAH)-selected young clusters in the PHANGS galaxy NGC 5068 show that our selection recovers ∼80% of the PAH-selected sources. We show that the practical limit for studying individual MYSOs with JWST is ∼3 Mpc. The resulting catalog provides a foundation for future resolved studies of star formation rates and early cluster evolution.

Abstract:

Context. High-resolution JWST-MIRI images allow us to resolve the multi-scale nature of the emission in nearby star-forming galaxies in great detail, from compact star-forming regions to large-scale diffuse emission, providing new insights into dust emission, its composition, and the surrounding interstellar medium (ISM). Aims. We aim to understand at which scale the different processes driving dust emission in mid-infrared (7.7−21 μm) wavelengths take place and whether we can disentangle emission from dense regions from the one linked to a more diffuse component. Methods. We use and enhance the Constrained Diffusion Decomposition algorithm (CDD), an alternative to wavelet transform decomposition, to disentangle the emission from compact regions from the emission originating from diffuse sources. This allows us to quantify the mid-IR spectral properties of the ISM at intervals within a continuum of physical scales. Results. We find a transition scale for Polycyclic Aromatic Hydrocarbons (PAH) emission around 300 pc, with a weaker PAH fraction at smaller scales, highlighting the destruction of PAHs in H II regions. We also show variations in the PAH fraction in different morphological environments, with a lower fraction in bright and star-forming environments. By studying and comparing the probability distribution functions (PDFs) of H II regions and of diffuse ISM with the PDFs at different scales, we find a similar separation scale around 200 pc at which we observe a transition from a power-law PDF for dense structures to a log-normal distribution for the diffuse ISM. Conclusions. We infer the destruction of PAHs in nebular environments and the existence of a transition scale between stellar feedback in nebular regions and heating in the diffuse ISM. We also show that the CDD enables future applications to study the physical scales of emission.

Abstract:

Context. The spatial distribution of giant molecular clouds (GMCs) at sub-kiloparsec scales encodes information about cloud formation and evolution. However, we still lack a general quantitative characterisation of molecular gas structure at this scale. Aims. We aim to provide a quantitative description of molecular gas structure at 150–1000 pc for a typical star-forming main sequence galaxy. We analyse how GMCs cluster together and how CO emission is spatially correlated with bright GMCs using a sample of 8984 GMCs from 40 galaxies observed by PHANGS-ALMA. Methods. We homogenized our data to a common spatial resolution of 150 pc and a mass sensitivity of 2.5 M ⊙ pc −2 to remove observational bias. We then calculated the nearest neighbour distances, neighbour number density, and two-point correlation functions (2PCFs) for the catalogued GMCs in each galaxy. When analysing the 2PCFs, we generated several control samples that reflect different null hypotheses on large spatial scales. We stacked integrated intensity CO emission profiles around the position of catalogued GMCs to probe the gas distribution on scales between the observational resolution and the typical GMC-GMC spacing. Results. Our measurements of cloud spacing and the number of neighbours show that GMC clustering follows the large-scale gas distribution. Once we accounted for this contribution, the peak excess clustering relative to the null hypothesis in the 2PCF dropped from 1 + ω ∼ 2.3 to 1.3, with the power-law slope flattened from −0.25 to 0. Stacks of CO intensity around local maxima show a strong clustering signal on scales smaller than the typical GMC-GMC separation. We show that this is largely the same signal captured by the ‘GMC size’ measured by CPROPS , with an additional ∼20% of the flux in an extended component beyond 500 pc. We find that our stacked profiles can be fit with a double Gaussian function plus a constant offset. The broad Gaussian component accounts for 70% of the over-density power above the constant background and is stronger around massive and gravitationally bound GMCs. Conclusions. Our measurements yield a general statistical description of the structure of CO emission from ≈150 pc to galactic scales that can serve as a benchmark for simulations of molecular cloud formation and destruction in galaxy disks. Our results indicate that galactic structure exerts a strong influence on the GMC distribution in galaxy disks and that the formation of massive gravitationally bound GMCs is related to strong local gas clustering.