skip to main content


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
NSF-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. We present optical photometry and spectroscopy of the Type II supernova ASASSN-14jb, together with Very Large Telescope (VLT) Multi Unit Spectroscopic Explorer (MUSE) integral field observations of its host galaxy and a nebular-phase spectrum. This supernova, in the nearby galaxy ESO 467-G051 ( z  = 0.006), was discovered and followed-up by the all-sky automated survey for supernovae (ASAS-SN). We obtained well-sampled las cumbres network (LCOGTN) B V g r i and Swift w 2 m 1 w 1 u b v optical, near-UV/optical light curves, and several optical spectra in the early photospheric phases. The transient ASASSN-14jb exploded ∼2 kpc above the star-forming disk of ESO 467-G051, an edge-on disk galaxy. The large projected distance from the disk of the supernova position and the non-detection of any H II region in a 1.4 kpc radius in projection are in conflict with the standard environment of core-collapse supernova progenitors and suggests the possible scenario that the progenitor received a kick in a binary interaction. We present analysis of the optical light curves and spectra, from which we derived a distance of 25 ± 2 Mpc using state-of-the-art empirical methods for Type II SNe, physical properties of the SN explosion ( 56 Ni mass, explosion energy, and ejected mass), and properties of the progenitor; namely the progenitor radius, mass, and metallicity. Our analysis yields a 56 Ni mass of 0.0210  ±  0.0025  M ⊙ , an explosion energy of ≈0.25 × 10 51 ergs, and an ejected mass of ≈6  M ⊙ . We also constrained the progenitor radius to be R *  = 580  ±  28  R ⊙ which seems to be consistent with the sub-Solar metallicity of 0.3  ±  0.1  Z ⊙ derived from the supernova Fe II λ 5018 line. The nebular spectrum constrains strongly the progenitor mass to be in the range 10–12 M ⊙ . From the Spitzer data archive we detect ASASSN-14jb ≈330 days past explosion and we derived a total dust mass of 10 −4   M ⊙ from the 3.6 μ m and 4.5 μ m photometry. Using the F U V , N U V , B V g r i , K s , 3.6 μ m, and 4.5 μ m total magnitudes for the host galaxy, we fit stellar population synthesis models, which give an estimate of M *  ≈ 1 × 10 9   M ⊙ , an age of 3.2 Gyr, and a SFR ≈0.07  M ⊙ yr −1 . We also discuss the low oxygen abundance of the host galaxy derived from the MUSE data, having an average of 12 + log(O/H) = 8.27 +0.16 −0.20 using the O 3 N 2 diagnostic with strong line methods. We compared it with the supernova spectra, which is also consistent with a sub-Solar metallicity progenitor. Following recent observations of extraplanar H II regions in nearby edge-on galaxies, we derived the metallicity offset from the disk, being positive, but consistent with zero at 2 σ , suggesting enrichment from disk outflows. We finally discuss the possible scenarios for the unusual environment for ASASSN-14jb and conclude that either the in-situ star formation or runaway scenario would imply a low-mass progenitor, agreeing with our estimate from the supernova nebular spectrum. Regardless of the true origin of ASASSN-14jb, we show that the detailed study of the environment roughly agree with the stronger constraints from the observation of the transient. 
    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 derive empirical constraints on the nucleosynthetic yields of nitrogen by incorporating N enrichment into our previously developed and empirically tuned multizone galactic chemical evolution model. We adopt a metallicity-independent (‘primary’) N yield from massive stars and a metallicity-dependent (‘secondary’) N yield from AGB stars. In our model, galactic radial zones do not evolve along the observed [N/O]–[O/H] relation, but first increase in [O/H] at roughly constant [N/O], then move upward in [N/O] via secondary N production. By t ≈ 5 Gyr, the model approaches an equilibrium [N/O]–[O/H] relation, which traces the radial oxygen gradient. Reproducing the [N/O]–[O/H] trend observed in extragalactic systems constrains the ratio of IMF-averaged N yields to the IMF-averaged O yield of core-collapse supernovae. We find good agreement if we adopt $y_\text{N}^\text{CC}/y_\text{O}^\text{CC}=0.024$ and $y_\text{N}^\text{AGB}/y_\text{O}^\text{CC} = 0.062(Z/Z_\odot)$. For the theoretical AGB yields we consider, simple stellar populations release half their N after only ∼250 Myr. Our model reproduces the [N/O]–[O/H] relation found for Milky Way stars in the APOGEE survey, and it reproduces (though imperfectly) the trends of stellar [N/O] with age and [O/Fe]. The metallicity-dependent yield plays the dominant role in shaping the gas-phase [N/O]–[O/H] relation, but the AGB time-delay is required to match the stellar age and [O/Fe] trends. If we add ∼40 per cent oscillations to the star formation rate, the model reproduces the scatter in the gas phase [N/O]–[O/H] relation observed in external galaxies by MaNGA. We discuss implications of our results for theoretical models of N production by massive stars and AGB stars.

     
    more » « less
  5. 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