skip to main content

Title: Momentum injection by clustered supernovae: testing subgrid feedback prescriptions

Using a 1D Lagrangian code specifically designed to assess the impact of multiple, time-resolved supernovae (SNe) from a single-star cluster on the surrounding medium, we test three commonly used feedback recipes: delayed cooling (e.g. used in the gasoline-2 code), momentum-energy injection (a resolution-dependent transition between momentum-dominated feedback and energy-dominated feedback used, e.g. in the fire-2 code), and simultaneous energy injection (e.g. used in the EAGLE simulations). Our work provides an intermediary test for these recipes: we analyse a setting that is more complex than the simplified scenarios for which many were designed, but one more controlled than a full galactic simulation. In particular, we test how well these models reproduce the enhanced momentum efficiency seen for an 11 SN cluster simulated at high resolution (0.6 pc; a factor of 12 enhancement relative to the isolated SN case) when these subgrid recipes are implemented in low resolution (20 pc) runs. We find that: (1) the delayed cooling model performs well – resulting in 9 times the momentum efficiency of the fiducial isolated SN value – when SNe are clustered and 1051 erg are injected per SN, while clearly overpredicting the momentum efficiency in the single SN test case; (2) the momentum-energy model always more » achieves good results, with a factor of 5 boost in momentum efficiency; and (3) injecting the energy from all SNe simultaneously does little to prevent overcooling and greatly underproduces the momentum deposited by clustered SNe, resulting in a factor of 3 decrease in momentum efficiency on the average.

« less
 ;  ;  
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 1243-1256
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: Gizmo , Arepo , Gadget , and Ramses . All codes adopt the same physical model, which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to subgrid models, eliminating one of the major uncertainties in cosmological simulations. We find reasonable agreement on the time-averaged star formation rates as well as the joint density–temperature distributions between all codes. However, the Lagrangian codes show significantly burstier star formation, larger SN-driven bubbles, and stronger galactic outflows compared to the Eulerian code. This is caused by the behavior in the dense, collapsing gas clouds when the Jeans length becomes unresolved: Gas in Lagrangian codes collapses to much higher densities than that in Eulerian codes, as the latter is stabilized by the minimal cell size. Therefore, more of the gas cloud is converted to stars and SNe are much more clustered in the Lagrangian models, amplifying their dynamical impact. The differences between Lagrangian and Eulerian codes can be reduced by adopting a higher star formation efficiency in Eulerian codes, which significantly enhances SN clustering in themore »latter. Adopting a zero SN delay time reduces burstiness in all codes, resulting in vanishing outflows as SN clustering is suppressed.« less
  2. Abstract Momentum feedback from isolated supernova remnants (SNRs) have been increasingly recognized by modern cosmological simulations as a resolution-independent means to implement the effects of feedback in galaxies, such as turbulence and winds. However, the integrated momentum yield from SNRs is uncertain due to the effects of SN clustering and interstellar medium (ISM) inhomogeneities. In this paper, we use spatially resolved observations of the prominent 10 kpc star-forming ring of M31 to test models of mass-weighted ISM turbulence driven by momentum feedback from isolated, nonoverlapping SNRs. We use a detailed stellar age distribution (SAD) map from the Panchromatic Hubble Andromeda Treasury survey, observationally constrained SN delay-time distributions, and maps of the atomic and molecular hydrogen to estimate the mass-weighted velocity dispersion using the Martizzi et al. ISM turbulence model. Our estimates are within a factor of two of the observed mass-weighted velocity dispersion in most of the ring, but exceed observations at densities ≲0.2 cm −3 and SN rates >2.1 × 10 −4 SN yr −1 kpc −2 , even after accounting for plausible variations in SAD models and ISM scale height assumptions. We conclude that at high SN rates the momentum deposited is most likely suppressed by the nonlinearmore »effects of SN clustering, while at low densities, SNRs reach pressure equilibrium before the cooling phase. These corrections should be introduced in models of momentum-driven feedback and ISM turbulence.« less
  3. Context. Recent developments in time domain astronomy, such as Zwicky Transient Facility (ZTF), have made it possible to conduct daily scans of the entire visible sky, leading to the discovery of hundreds of new transients every night. Among these detections, 10 to 15 of these objects are supernovae (SNe), which have to be classified prior to cosmological use. The spectral energy distribution machine (SEDM) is a low-resolution ( ℛ ~ 100) integral field spectrograph designed, built, and operated with the aim of spectroscopically observing and classifying targets detected by the ZTF main camera. Aims. As the current pysedm pipeline can only handle isolated point sources, it is limited by contamination when the transient is too close to its host galaxy core. This can lead to an incorrect typing and ultimately bias the cosmological analyses, affecting the homogeneity of the SN sample in terms of local environment properties. We present a new scene modeler to extract the transient spectrum from its structured background, with the aim of improving the typing efficiency of the SEDM. Methods. H yper G al is a fully chromatic scene modeler that uses archival pre-transient photometric images of the SN environment to generate a hyperspectral model ofmore »the host galaxy. It is based on the cigale SED fitter used as a physically-motivated spectral interpolator. The galaxy model, complemented by a point source for the transient and a diffuse background component, is projected onto the SEDM spectro-spatial observation space and adjusted to observations, and the SN spectrum is ultimately extracted from this multi-component model. The full procedure, from scene modeling to transient spectrum extraction and typing, is validated on 5000 simulated cubes built from actual SEDM observations of isolated host galaxies, covering a broad range of observing conditions and scene parameters. Results. We introduce the contrast, c , as the transient-to-total flux ratio at the SN location, integrated over the ZTF r -band. From estimated contrast distribution of real SEDm observations, we show that H yper G al correctly classifies ~95% of SNe Ia, and up to 99% for contrast c ≳ 0.2, representing more than 90% of the observations. Compared to the standard point-source extraction method (without the hyperspectral galaxy modeling step), H yper G al correctly classifies 20% more SNe Ia between 0.1 < c < 0.6 (50% of the observation conditions), with less than 5% of SN Ia misidentifications. The false-positive rate is less than 2% for c > 0.1 (> 99% of the observations), which represents half as much as the standard extraction method. Assuming a similar contrast distribution for core-collapse SNe, H yper G al classifies 14% additional SNe II and 11% additional SNe Ibc. Conclusions. H yper G al has proven to be extremely effective in extracting and classifying SNe in the presence of strong contamination by the host galaxy, providing a significant improvement with respect to the single point-source extraction.« less
  4. Abstract The quenching “maintenance” and “cooling flow” problems are important from the Milky Way through massive cluster elliptical galaxies. Previous work has shown that some source of energy beyond that from stars and pure magnetohydrodynamic processes is required, perhaps from AGN, but even the qualitative form of this energetic input remains uncertain. Different scenarios include thermal “heating,” direct wind or momentum injection, cosmic ray heating or pressure support, or turbulent “stirring” of the intra-cluster medium (ICM). We investigate these in 1012 − 1014 M⊙ halos using high-resolution non-cosmological simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, including simplified toy energy-injection models, where we arbitrarily vary the strength, injection scale, and physical form of the energy. We explore which scenarios can quench without violating observational constraints on energetics or ICM gas. We show that turbulent stirring in the central ∼100 kpc, or cosmic-ray injection, can both maintain a stable low-SFR halo for >Gyr timescales with modest energy input, by providing a non-thermal pressure which stably lowers the core density and cooling rates. In both cases, associated thermal-heating processes are negligible. Turbulent stirring preserves cool-core features while mixing condensed core gas into the hotter halo and is by far the mostmore »energy efficient model. Pure thermal heating or nuclear isotropic momentum injection require vastly larger energy, are less efficient in lower-mass halos, easily over-heat cores, and require fine-tuning to avoid driving unphysical temperature gradients or gas expulsion from the halo center.« less
  5. Context. The inclusion of molecular physics is an important piece that tends to be missing from the puzzle when modeling the spectra of supernovae (SNe). Molecules have both a direct impact on the spectra, particularly in the infrared, and an indirect one as a result of their influence on certain physical conditions, such as temperature. Aims. In this paper, we aim to investigate molecular formation and non-local thermodynamic equilibrium (NLTE) cooling, with a particular focus on CO, the most commonly detected molecule in supernovae. We also aim to determine the dependency of supernova chemistry on physical parameters and the relative sensitivity to rate uncertainties. Methods. We implemented a chemical kinetic description of the destruction and formation of molecules into the SN spectral synthesis code SUMO . In addition, selected molecules were coupled into the full NLTE level population framework and, thus, we incorporated molecular NLTE cooling into the temperature equation. We produced a test model of the CO formation in SN 1987A between 150 and 600 days and investigated the sensitivity of the resulting molecular masses to the input parameters. Results. We find that there is a close inter-dependency between the thermal evolution and the amount of CO formed, mainlymore »through an important temperature-sensitive CO destruction process with O + . After a few hundred days, CO completely dominates the cooling of the oxygen-carbon zone of the supernova which, therefore, contributes little optical emission. The uncertainty of the calculated CO mass scales approximately linearly with the typical uncertainty factor for individual rates. We demonstrate how molecular masses can potentially be used to constrain various physical parameters of the supernova.« less