skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Exploring metallicity-dependent rates of Type Ia supernovae and their impact on galaxy formation
Abstract Type Ia supernovae are critical for feedback and elemental enrichment in galaxies. Recent surveys like the All-Sky Automated Survey for Supernova (ASAS-SN) and the Dark Energy Survey (DES) find that the specific supernova Ia rate at z ∼ 0 may be ≲ 20 − 50 × higher in lower-mass galaxies than at Milky Way-mass. Independently, observations show that the close-binary fraction of solar-type Milky Way stars is higher at lower metallicity. Motivated by these observations, we use the FIRE-2 cosmological zoom-in simulations to explore the impact of metallicity-dependent rate models on galaxies of $$M_* \sim 10^7\, \rm {M}_{\odot }-10^{11}\, \rm {M}_{\odot }$$. First, we benchmark our simulated star-formation histories (SFHs) against observations, and show that assumed stellar mass functions play a major role in determining the degree of tension between observations and metallicity-independent rate models, potentially causing ASAS-SN and DES observations to agree more than might appear. Models in which the supernova Ia rate increases with decreasing metallicity ($$\propto Z^{-0.5 \; \rm {to} \; -1}$$) provide significantly better agreement with observations. Encouragingly, these rate increases (≳ 10 × in low-mass galaxies) do not significantly impact galaxy masses and morphologies, which remain largely unaffected except for our most extreme models. We explore implications for both [Fe/H] and [$$\alpha /\rm {Fe}$$] enrichment; metallicity-dependent rate models can improve agreement with the observed stellar mass-metallicity relations in low-mass galaxies. Our results demonstrate that a range of metallicity-dependent rate models are viable for galaxy formation and motivate future work.  more » « less
Award ID(s):
2045928 2108318 2108230 1715216 1715847 1713353
PAR ID:
10351389
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
ISSN:
0035-8711
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The scaling of the specific Type Ia supernova (SN Ia) rate with host galaxy stellar mass $$\dot{\text{N}}_\text{Ia} / \text{M}_\star \sim \text{M}_\star ^{-0.3}$$ as measured in ASAS-SN and DES strongly suggests that the number of SNe Ia produced by a stellar population depends inversely on its metallicity. We estimate the strength of the required metallicity dependence by combining the average star formation histories (SFHs) of galaxies as a function of their stellar mass with the mass–metallicity relation (MZR) for galaxies and common parametrizations for the SN Ia delay-time distribution. The differences in SFHs can account for only ∼30 per cent of the increase in the specific SN Ia rate between stellar masses of M⋆ = 1010 and 107.2 M⊙. We find that an additional metallicity dependence of approximately ∼Z−0.5 is required to explain the observed scaling. This scaling matches the metallicity dependence of the close binary fraction observed in APOGEE, suggesting that the enhanced SN Ia rate in low-mass galaxies can be explained by a combination of their more extended SFHs and a higher binary fraction due to their lower metallicities. Due to the shape of the MZR, only galaxies below M⋆ ≈ 3 × 109 M⊙ are significantly affected by the metallicity-dependent SN Ia rates. The $$\dot{\text{N}}_\text{Ia} / \text{M}_\star \sim \text{M}_\star ^{-0.3}$$ scaling becomes shallower with increasing redshift, dropping by factor of ∼2 at 107.2 M⊙ between z = 0 and 1 with our ∼Z−0.5 scaling. With metallicity-independent rates, this decrease is a factor of ∼3. We discuss the implications of metallicity-dependent SN Ia rates for one-zone models of galactic chemical evolution. 
    more » « less
  2. ABSTRACT We explore the effect of including progenitor mass- and metallicity-dependent yields, supernova rates and energetics on variations in elemental abundance ratios (particularly [α/Fe]) in dwarf galaxies. To understand how the scatter and overall trends in [α/Fe] are affected by including variable metal yields from a discretely sampled initial mass function, we run FIRE simulations of a dwarf galaxy (M⋆(z = 0$$) \sim 10^6\rm \, M_{\odot })$$ using nucleosynthetic yields from the NuGrid data base that depend on the stellar progenitor mass and metallicity. While NuGrid exhibits lower aggregate α-element production than default FIRE yields, we find that its explicit mass dependence, even when including turbulent metal diffusion, substantially widens the intrinsic scatter in the simulated [Fe/H]-[α/Fe] – a phenomenon visible in some observations of dwarf galaxies. 
    more » « less
  3. ABSTRACT We develop a hybrid model of galactic chemical evolution that combines a multiring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradient in outflow mass loading tuned to reproduce the observed metallicity gradient. With this approach, the model reproduces many of the striking qualitative features of the Milky Way disc’s abundance structure: (i) the dependence of the [O/Fe]–[Fe/H] distribution on radius Rgal and mid-plane distance |z|; (ii) the changing shapes of the [O/H] and [Fe/H] distributions with Rgal and |z|; (iii) a broad distribution of [O/Fe] at sub-solar metallicity and changes in the [O/Fe] distribution with Rgal, |z|, and [Fe/H]; (iv) a tight correlation between [O/Fe] and stellar age for [O/Fe] > 0.1; (v) a population of young and intermediate-age α-enhanced stars caused by migration-induced variability in the Type Ia supernova rate; (vi) non-monotonic age–[O/H] and age–[Fe/H] relations, with large scatter and a median age of ∼4 Gyr near solar metallicity. Observationally motivated models with an enhanced star formation rate ∼2 Gyr ago improve agreement with the observed age–[Fe/H] and age–[O/H] relations, but worsen agreement with the observed age–[O/Fe] relation. None of our models predict an [O/Fe] distribution with the distinct bimodality seen in the observations, suggesting that more dramatic evolutionary pathways are required. All code and tables used for our models are publicly available through the Versatile Integrator for Chemical Evolution (VICE; https://pypi.org/project/vice). 
    more » « less
  4. ABSTRACT We present improved photometric measurements for the host galaxies of 206 spectroscopically confirmed type Ia supernovae discovered by the Dark Energy Survey Supernova Program (DES-SN) and used in the first DES-SN cosmological analysis. For the DES-SN sample, when considering a 5D (z, x1, c, α, β) bias correction, we find evidence of a Hubble residual ‘mass step’, where SNe Ia in high-mass galaxies (>1010M⊙) are intrinsically more luminous (after correction) than their low-mass counterparts by $$\gamma =0.040\pm 0.019$$ mag. This value is larger by 0.031 mag than the value found in the first DES-SN cosmological analysis. This difference is due to a combination of updated photometric measurements and improved star formation histories and is not from host-galaxy misidentification. When using a 1D (redshift-only) bias correction the inferred mass step is larger, with $$\gamma =0.066\pm 0.020$$ mag. The 1D−5D γ difference for DES-SN is $$0.026\pm 0.009$$ mag. We show that this difference is due to a strong correlation between host galaxy stellar mass and the x1 component of the 5D distance-bias correction. Including an intrinsic correlation between the observed properties of SNe Ia, stretch and colour, and stellar mass in simulated SN Ia samples, we show that a 5D fit recovers γ with −9 mmag bias compared to a +2 mmag bias for a 1D fit. This difference can explain part of the discrepancy seen in the data. Improvements in modelling correlations between galaxy properties and SN is necessary to ensure unbiased precision estimates of the dark energy equation of state as we enter the era of LSST. 
    more » « less
  5. ABSTRACT Motivated by the early excess of bright galaxies seen by JWST, we run zoom-in cosmological simulations of a massive galaxy at Cosmic Dawn, in a halo of $$10^{11} {\rm M}_\odot$$ at $z = 9$, using the hydro-gravitational code ramses at an effective resolution $$\sim 10~{\rm pc}$$. We investigate physical mechanisms that enhance the star formation efficiencies (SFEs) at the high gas densities of the star-forming regions in this galaxy ($$\sim 3\times 10^3~{\rm cm^{-3}}$$, $$\sim 10^4~{\rm M}_\odot \,{\rm pc^{-2}}$$). Our fiducial star formation recipe uses a physically motivated, turbulence-based, multi-freefall model, avoiding ad hoc extrapolation from lower redshifts. By $z = 9$, our simulated galaxy is a clumpy, thick, rotating disc with a high stellar mass $$\sim 3\times 10^9~{\rm M}_\odot$$ and high star formation rate $$\sim 50~{\rm M}_\odot \,{\rm yr^{-1}}$$. The high gas density makes supernova (SN) feedback less efficient, producing a high local SFE $$\gtrsim 10~{{\ \rm per\ cent}}$$. The global SFE is set by feedback-driven outflows and only weakly correlated with the local SFE. Photoionization heating makes SN feedback more efficient, but the integrated SFE always remains high. Intense accretion at Cosmic Dawn seeds turbulence that reduces local SFE, but this only weakly affects the global SFE. The star formation histories of our simulated galaxies are similar to observed massive galaxies at Cosmic Dawn, despite our limited resolution. We set the stage for future simulations which treat radiation self-consistently and use a higher effective resolution $$\sim 1~{\rm pc}$$ that captures the physics of star-forming clouds. 
    more » « less