Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 5:00 PM ET until 8:00 PM ET on Friday, September 11 due to maintenance. We apologize for the inconvenience.


This content will become publicly available on March 12, 2027

Title: Where Do Stars Explode in the ISM?—The Distribution of Dense Gas around Evolved Massive Stars in M33
Abstract The effect of supernovae (SNe) on star formation in the interstellar medium (ISM) depends sensitively on where SNe explode with respect to ISM clouds. Observationally, SN ISM environments characterized by spatially resolved gas maps can empirically guide the placement of SNe in subgrid models, but unfortunately such measurements remain scarce, as SNe are rare and often distant. Here we demonstrate a new approach—mapping the ISM around evolved massive stars that are soon to explode. These provide a substantially larger sample of “explosion sites” (than just historical SNe) in nearby galaxies that have high-resolution atomic and molecular ISM maps from the Jansky Very Large Array and Atacama Large Millimeter/submillimeter Array. We demonstrate this technique in the well-resolved Local Group spiral M33 by analyzing the 50 pc scale projected ISM densities around red supergiants (RSGs; 8–30Mstars) Wolf–Rayet stars (W-Rs; >30Mstars), and supernova remnants. We find amass-dependentcorrelation between stars and gas clouds, with at least 45% of W-Rs and up to 77% of RSGs having no detectable H2at their pixel locations. In the sample with H2detections, we find that more-massive younger progenitors are coincident with denser gas. We show that the density distributions for stars >15Mare statistically distinct from random alignment of stars and gas in M33. Our work provides the firstobservationally derivedestimate of the fraction of the SN-producing stellar population correlated with ISM density peaks. We demonstrate how this can be compared with galaxy simulations, and advocate similar comparisons to the community for constraining subgrid models.  more » « less
Award ID(s):
2205631 2205628
PAR ID:
10703064
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; « less
Publisher / Repository:
IOP
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
1000
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
70
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Star formation in galaxies is regulated by turbulence, outflows, gas heating and cloud dispersal -- processes which depend sensitively on the properties of the interstellar medium (ISM) into which supernovae (SNe) explode. Unfortunately, direct measurements of ISM environments around SNe remain scarce, as SNe are rare and often distant. Here we demonstrate a new approach: mapping the ISM around the massive stars that are soon to explode. This provides a much larger census of explosion sites than possible with only SNe, and allows comparison with sensitive, high-resolution maps of the atomic and molecular gas from the Jansky VLA and ALMA. In the well-resolved Local Group spiral M33, we specifically observe the environments of red supergiants (RSGs, progenitors of Type II SNe), Wolf-Rayet stars (WRs, tracing stars >30 M⊙, and possibly future stripped-envelope SNe), and supernova remnants (SNRs, locations where SNe have exploded). We find that massive stars evolve not only in dense, molecular-dominated gas (with younger stars in denser gas), but also a substantial fraction (∼45\% of WRs; higher for RSGs) evolve in lower-density, atomic-gas-dominated, inter-cloud media. We show that these measurements are consistent with expectations from different stellar-age tracer maps, and can be useful for validating SN feedback models in numerical simulations of galaxies. Along with the discovery of a 20-pc diameter molecular gas cavity around a WR, these findings re-emphasize the importance of pre-SN/correlated-SN feedback evacuating the dense gas around massive stars before explosion, and the need for high-resolution (down to pc-scale) surveys of the multi-phase ISM in nearby galaxies. 
    more » « less
  2. Abstract We measure the molecular gas environment near recent (<100 yr old) supernovae (SNe) using ∼1″ or ≤150 pc resolution CO (2–1) maps from the PHANGS–Atacama Large Millimeter/submillimeter Array (ALMA) survey of nearby star-forming galaxies. This is arguably the first such study to approach the scales of individual massive molecular clouds (Mmol≳ 105.3M). Using the Open Supernova Catalog, we identify 63 SNe within the PHANGS–ALMA footprint. We detect CO (2–1) emission near ∼60% of the sample at 150 pc resolution, compared to ∼35% of map pixels with CO (2–1) emission, and up to ∼95% of the SNe at 1 kpc resolution, compared to ∼80% of map pixels with CO (2–1) emission. We expect the ∼60% of SNe within the same 150 pc beam, as a giant molecular cloud will likely interact with these clouds in the future, consistent with the observation of widespread SN–molecular gas interaction in the Milky Way, while the other ∼40% of SNe without strong CO (2–1) detections will deposit their energy in the diffuse interstellar medium, perhaps helping drive large-scale turbulence or galactic outflows. Broken down by type, we detect CO (2–1) emission at the sites of ∼85% of our 9 stripped-envelope SNe (SESNe), ∼40% of our 34 Type II SNe, and ∼35% of our 13 Type Ia SNe, indicating that SESNe are most closely associated with the brightest CO (2–1) emitting regions in our sample. Our results confirm that SN explosions are not restricted to only the densest gas, and instead exert feedback across a wide range of molecular gas densities. 
    more » « less
  3. Context.Feedback from stars in the form of radiation, stellar winds, and supernovae is crucial to regulating the star formation activity of galaxies. Dwarf galaxies are especially susceptible to these processes, making them an ideal test bed for studying the effects of stellar feedback in detail. Recent numerical models have aimed to resolve the interstellar medium (ISM) in dwarf galaxies with a very high resolution of several solar masses. However, when it comes to modeling the radiative feedback from stars, many models opt for simplified approaches instead of explicitly solving radiative transfer (RT) because of the computational complexity involved. Aims.We introduce the Realistic ISM modeling in Galaxy Evolution and Lifecycles (RIGEL) model, a novel framework to self-consistently model the effects of stellar feedback in the multiphase ISM of dwarf galaxies with explicit RT on a star-by-star basis. Methods.The RIGEL model integrates detailed implementations of feedback from individual massive stars into the state-of-the-art radiation-hydrodynamics code,AREPO-RT. It forms individual massive stars from the resolved multiphase ISM by sampling the initial mass function and tracks their evolution individually. The lifetimes, photon production rates, mass-loss rates, and wind velocities of these stars are determined by their initial masses and metallicities based on a library that incorporates a variety of stellar models. The RT equations are solved explicitly in seven spectral bins accounting for the infrared to He IIionizing bands, using a moment-base scheme with the M1 closure relation. The thermochemistry model tracks the nonequilibrium H, He chemistry as well as the equilibrium abundance of C I, C II, O I, O II, and CO in the irradiated ISM to capture the thermodynamics of all ISM phases, from cold molecular gas to hot ionized gas. Results.We evaluated the performance of the RIGEL model using 1 Mresolution simulations of isolated dwarf galaxies. We found that the star formation rate (SFR) and interstellar radiation field (ISRF) show strong positive correlations with the metallicity of the galaxy. Photoionization and photoheating can reduce the SFR by an order of magnitude by removing the available cold, dense gas fuel for star formation. The presence of ISRF also significantly changes the thermal structure of the ISM. Radiative feedback occurs immediately after the birth of massive stars and rapidly disperses the molecular clouds within 1 Myr. As a consequence, radiative feedback reduces the age spread of star clusters to less than 2 Myr, prohibits the formation of massive star clusters, and shapes the cluster initial mass function to a steep power-law form with a slope of ∼ − 2. The mass-loading factor (measured atz = 1 kpc) of the fiducial galaxy has a median ofηM ∼ 50, while turning off radiative feedback reduces this factor by an order of magnitude. Conclusions.We demonstrate that RIGEL effectively captures the nonlinear coupling of early radiative feedback and supernova feedback in the multiphase ISM of dwarf galaxies. This novel framework enables the utilization of a comprehensive stellar feedback and ISM model in cosmological simulations of dwarf galaxies and various galactic environments spanning a wide dynamic range in both space and time. 
    more » « less
  4. Abstract Unveiling the chemical fingerprints of the first (Population III, hereafter Pop III) stars is crucial for indirectly studying their properties and probing their massive nature. In particular, very massive Pop III stars explode as energetic pair-instability supernovae (PISNe), allowing their chemical products to escape in the diffuse medium around galaxies, opening the possibility to observe their fingerprints in distant gas clouds. Recently, threez> 6.3 absorbers with abundances consistent with an enrichment from PISNe have been observed with JWST. In this Letter, we present novel chemical diagnostics to uncover environments mainly imprinted by PISNe. Furthermore, we revise the JWST low-resolution measurements by analyzing the publicly available high-resolution X-Shooter spectra for two of these systems. Our results reconcile the chemical abundances of these absorbers with those from literature, which are found to be consistent with an enrichment dominated (>50% metals) by normal Pop II SNe. We show the power of our novel diagnostics in isolating environments uniquely enriched by PISNe from those mainly polluted by other Pop III and Pop II SNe. When the subsequent enrichment from Pop II SNe is included, however, we find that the abundances of PISN-dominated environments partially overlap with those predominantly enriched by other Pop III and Pop II SNe. We dub these areas confusion regions. Yet, the odd–even abundance ratios [Mg,Si/Al] are extremely effective in pinpointing PISN-dominated environments and allowed us to uncover, for the first time, an absorber consistent with a combined enrichment by a PISN and another Pop III SN for all the six measured elements. 
    more » « less
  5. Abstract We introduce new high-resolution galaxy simulations accelerated by a surrogate model that reduces the computation cost by approximately 75%. Massive stars with a zero-age main-sequence mass of more than about 10Mexplode as core-collapse supernovae (CCSNe), which play a critical role in galaxy formation. The energy released by CCSNe is essential for regulating star formation and driving feedback processes in the interstellar medium (ISM). However, the short integration time steps required for SN feedback have presented significant bottlenecks in astrophysical simulations across various scales. Overcoming this challenge is crucial for enabling star-by-star galaxy simulations, which aim to capture the dynamics of individual stars and the inhomogeneous shell’s expansion within the turbulent ISM. To address this, our new framework combines direct numerical simulations and surrogate modeling, including machine learning and Gibbs sampling. The star formation history and the time evolution of outflow rates in the galaxy match those obtained from resolved direct numerical simulations. Our new approach achieves high-resolution fidelity while reducing computational costs, effectively bridging the physical scale gap and enabling multiscale simulations. 
    more » « less