Dr Thomas Williams

Abstract:

We present the first spatially resolved infrared images of supernova remnants (SNRs) in M33 with the unprecedented sensitivity and resolution of JWST. We analyze 40 SNRs in four JWST fields: two covering central and southern M33 with separate NIRCam (F335M, F444W) and MIRI (F560W, F2100W) observations, one ∼5 kpc-long radial strip observed with MIRI F770W, and one covering the giant H II region NGC 604 with multiple NIRCam and MIRI broad/narrowband filters. Of the 21 SNRs in the MIRI (F560W+F2100W) field, we found three clear detections (i.e., identical infrared and Hα morphologies), and six partial-detections, implying a detection fraction of 43% in these bands. One of the SNRs in this field, L10-080, is a potential candidate for having freshly formed ejecta dust, based on its size and centrally concentrated 21 μm emission. In contrast, only one SNR (out of 16) is detectable in the NIRCam F335M+F444W field. Two SNRs near NGC 604 have strong evidence of molecular (H2) emission at 4.7 μm, making them the farthest known SNRs with visible molecular shocks. Five SNRs have F770W observations, with the smaller younger objects showing tentative signs of emission, while the older, larger ones have voids. Multiwavelength data indicate that the clearly detected SNRs are also among the smallest, brightest at other wavelengths (Hα, radio, and X-ray), have the broadest line widths (Hα FWHM ∼ 250–350 km s−1), and the densest environments. No correlation between the JWST-detectability and local star formation history of the SNRs is apparent.

Abstract:

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.

Abstract:

New extragalactic measurements of the cloud population-averaged star formation efficiency per free-fall time, ϵ ff , from PHANGS show little sign of a theoretically predicted dependence on the gas virial level and weak variation with cloud-scale gas velocity dispersion. We explore ways to bring theory into consistency with the observations, particularly by highlighting systematic variations in internal density structure that must accompany an increase in virial parameter typically found toward denser galaxy centers. To introduce these variations into conventional turbulence-regulated star formation models, we adopted three adjustments, all motivated by the expectation that the background host galaxy has an influence on the cloud scale: (1) We incorporate self-gravity and an internal density distribution that contains a broad power-law (PL) component and resembles the structure observed in local resolved clouds; (2) We allow the internal gas kinematics to include motion in the background potential and let this regulate the onset of self-gravitation; (3) We assume that the distribution of gas densities is in a steady state for only a fraction of a cloud free-fall time. In practice, these changes significantly reduce the efficiencies predicted in multi-free-fall (MFF) scenarios compared to purely lognormal probability density functions (PDFs) and tie efficiency variations to variations in the slope of the PL α . We fit the model to PHANGS measurements of ϵ ff to identify the PL slopes that yield an optimal match. These slopes vary systematically with galactic environment in the sense that gas that sits furthest from virial balance contains fractionally more gas at high density. We relate this to the equilibrium response of gas in the presence of the galactic gravitational potential, which forces more gas to high density than characteristic of fully self-gravitating clouds. Viewing the efficiency variations as originating with time evolution in the PL slope, our findings would alternatively imply coordination of the cloud evolutionary stage within environment. With this “galaxy regulation” behavior included, our preferred “self-gravitating” multi-freefall sgMFF models function similarly to the original, roughly “virialized cloud” single-free-fall models. However, outside the environment of disks with their characteristic regulation, the flexible MFF models may be better suited.

Abstract:

Recent hydrodynamical simulations of isolated barred disc galaxies have suggested a structural change in the distribution of the interstellar medium (ISM) around a stellar mass M * of 10 10 M ⊙ . In the higher-mass regime ( M ∗ ≥ 10 10 M ⊙ ), we observe the formation of a central gas and stellar disc with a typical size of a few hundred parsecs connected through lanes to the ends of the stellar bar. In the lower-mass regime ( M ∗ < 10 10 M ⊙ ), such an inner disc is absent and the gas component exhibits a more chaotic distribution. Observations of nearby star-forming galaxies support the existence of such a change. These inner gas discs may represent an important intermediate scale connecting the large kiloparsec-scale structures with the nuclear (sub-parsec) region, transporting gas inwards to fuel the central supermassive black hole (SMBH). For this work we used an extended set of high-resolution hydrodynamical simulations of isolated disc galaxies with initial properties (i.e. stellar mass, gas fraction, stellar disc scale length, and the bulge mass fraction) with properties covering the range of galaxies in the PHANGS sample to investigate this change of regime. We studied the physical properties of the star-forming ISM in both stellar mass regimes and extracted a few physical tracers: the inner Lindblad resonance (ILR), the probability distribution function (PDF), the virial parameter, and the Mach number. In line with observations, we confirm a structure transition in the simulations that occurs between a stellar mass of 10 9.5 and 10 10 M ⊙ . We show that the physical origin of this change of regime is driven by stellar feedback and its contribution relative to the underlying gravitational potential. With their shallower potential and typically higher gas mass fraction, lower-mass disc PHANGS galaxies combine two ingredients that significantly delay or even prevent the formation of a central gas (and stellar) disc. These results shed some light on the observed properties of star-forming barred galaxies and have implications for the star formation regimes, the growth of central structures, and the overall secular evolution of disc galaxies.

Abstract:

The Sculptor galaxy, NGC 253, is the southern massive star-forming disk galaxy that is closest to the Milky Way. We present a new 103-pointing MUSE mosaic of this galaxy that covers most of its star-forming disk up to 0.75 × R 25 . With an area of ∼20 × 5 arcmin 2 (∼20 × 5 kpc 2 , projected) and a physical resolution of ∼15 pc, this mosaic constitutes one of the largest integral field spectroscopy surveys with the highest physical resolution of any star-forming galaxy to date. We exploited the mosaic to identify a sample of ∼500 planetary nebulae (the sample is ∼20 times larger than in previous studies) to build the planetary nebula luminosity function (PNLF) and obtain a new estimate of the distance to NGC 253. The value we obtained is 17% higher than the estimates returned by other reliable measurements, which were mainly obtained via the top of the red giant branch method. The PNLF also varies between the centre ( r < 4 kpc) and the disk of the galaxy. The distance derived from the PNLF of the outer disk is comparable to that of the full sample, while the PNLF of the centre returns a distance that is larger by ∼0.9 Mpc. Our analysis suggests that extinction related to the dust-rich interstellar medium and edge-on view of the galaxy (the average E ( B − V ) across the disk is ∼0.35 mag) plays a major role in explaining both the larger distance recovered from the full PNLF and the difference between the PNLFs in the centre and the disk.