No Arabic abstract
Metallicity is expected to influence not only the lives of massive stars but also the outcome of their deaths as supernovae (SNe) and as gamma-ray bursts (GRBs). However, there are surprisingly few direct measurements of the local metallicities of different flavors of core-collapse SNe. Here we present the largest existing set of host-galaxy spectra with H II region emission lines at the sites of 35 stripped-envelope core-collapse SNe. We derive local oxygen abundances in a robust manner in order to constrain the SN Ib/c progenitor population. We obtain spectra at the SN sites, include SNe from targeted and untargeted surveys, and perform the abundance determinatinos using three different oxygen-abundance calibrations. The sites of SNe Ic (the demise of the most heavily stripped stars having lost both the H and He layers) are systematically more metal rich than those of SNe Ib (arising from stars that retained their He layer) in all calibrations. A Kolmogorov-Smirnov test yields the very low probability of 1% that SN Ib and SN Ic environment abundances, which are different on average by ~0.2 dex (in the Pettini & Pagel scale), are drawn from the same parent population. Broad-lined SNe Ic (without GRBs) occur at metallicities between those of SNe Ib and SNe Ic. Lastly, we find that the host-galaxy central oxygen abundance is not a good indicator of the local SN metallicity; hence, large-scale SN surveys need to obtain local abundance measurements in order to quantify the impact of metallicity on stellar death.
We present modelling of line polarization to study multi-dimensional geometry of stripped-envelope core-collapse supernovae (SNe). We demonstrate that a purely axisymmetric, two-dimensional geometry cannot reproduce a loop in the Stokes Q-U diagram, i.e., a variation of the polarization angles along the velocities associated with the absorption lines. On the contrary, three-dimensional (3D) clumpy structures naturally reproduce the loop. The fact that the loop is commonly observed in stripped-envelope SNe suggests that SN ejecta generally have a 3D structure. We study the degree of line polarization as a function of the absorption depth for various 3D clumpy models with different clump sizes and covering factors. Comparison between the calculated and observed degree of line polarization indicates that a typical size of the clump is relatively large, >~ 25 % of the photospheric radius. Such large-scale clumps are similar to those observed in the SN remnant Cassiopeia A. Given the small size of the observed sample, the covering factor of the clumps is only weakly constrained (~ 5-80 %). The presence of large-scale clumpy structure suggests that the large-scale convection or standing accretion shock instability takes place at the onset of the explosion.
We investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 successful explosions during several seconds. We present a broad study based on three progenitors (11.2 $M_odot$, 15 $M_odot$, and 27 $M_odot$), different neutrino-heating efficiencies, and various rotation rates. We show that the first seconds after shock revival determine the final explosion energy, remnant mass, and properties of ejected matter. Our results suggest that a continued mass accretion increases the explosion energy even at late times. We link the late-time mass accretion to initial conditions such as rotation strength and shock deformation at explosion time. Only some of our simulations develop a neutrino-driven wind that survives for several seconds. This indicates that neutrino-driven winds are not a standard feature expected after every successful explosion. Even if our neutrino treatment is simple, we estimate the nucleosynthesis of the exploding models for the 15 $M_odot$ progenitor after correcting the neutrino energies and luminosities to get a more realistic electron fraction.
We infer the progenitor mass distribution for 22 historic core-collapse supernovae (CCSNe) using a Bayesian hierarchical model. For this inference, we use the local star formation histories to estimate the age for each supernova (SN). These star formation histories often show multiple bursts of star formation; our model assumes that one burst is associated with the SN progenitor and the others are random bursts of star formation. The primary inference is the progenitor age distribution. Due to the limited number of historic SNe and highly uncertain star formation at young ages, we restrict our inference to the slope of the age distribution and the maximum age for CCSNe. Using single-star evolutionary models, we transform the progenitor age distribution into a progenitor mass distribution. Under these assumptions, the minimum mass for CCSNe is ${M_textrm{min}}~=~8.60^{+0.37}_{-0.41} M_odot$ and the slope of the progenitor mass distribution is $alpha = -2.61^{+1.05}_{-1.18}$. The power-law slope for the progenitor mass distribution is consistent with the standard Salpeter initial mass function ($alpha = -2.35$). These values are consistent with previous estimates using precursor imaging and the age-dating technique, further confirming that using stellar populations around SN and supernova remnants is a reliable way to infer the progenitor masses.
We age-date the stellar populations associated with 12 historic nearby core-collapse supernovae (CCSNe) and 2 supernova impostors, and from these ages, we infer their initial masses and associated uncertainties. To do this, we have obtained new HST imaging covering these CCSNe. Using these images, we measure resolved stellar photometry for the stars surrounding the locations of the SNe. We then fit the color-magnitude distributions of this photometry with stellar evolution models to determine the ages of any young existing populations present. From these age distributions, we infer the most likely progenitor mass for all of the SNe in our sample. We find ages between 4 and 50 Myr, corresponding to masses from 7.5 to 59 solar masses. There were no SNe that lacked a young population within 50~pc. Our sample contains 4 type Ib/c SNe; their masses have a wide range of values, suggesting that the progenitors of stripped-envelope SNe are binary systems. Both impostors have masses constrained to be $lesssim$7.5 solar masses. In cases with precursor imaging measurements, we find that age-dating and precursor imaging give consistent progenitor masses. This consistency implies that, although the uncertainties for each technique are significantly different, the results of both are reliable to the measured uncertainties. We combine these new measurements with those from our previous work and find that the distribution of 25 core-collapse SNe progenitor masses is consistent with a standard Salpeter power-law mass function, no upper mass cutoff, and an assumed minimum mass for core-collapse of 7.5~M$_{odot}$.
Recent multi-dimensional simulations of core-collapse supernovae are producing successful explosions and explosion-energy predictions. In general, the explosion-energy evolution is monotonic and relatively smooth, suggesting a possible analytic solution. We derive analytic solutions for the expansion of the gain region under the following assumptions: spherical symmetry, one-zone shell, and powered by neutrinos and $alpha$ particle recombination. We consider two hypotheses: I) explosion energy is powered by neutrinos and $alpha$ recombination, II) explosion energy is powered by neutrinos alone. Under these assumptions, we derive the fundamental dimensionless parameters and analytic scalings. For the neutrino-only hypothesis (II), the asymptotic explosion energy scales as $E_{infty} approx 1.5 M_g v_0^2 eta^{2/3}$, where $M_g$ is the gain mass, $v_0$ is the free-fall velocity at the shock, and $eta$ is a ratio of the heating and dynamical time scales. Including both neutrinos and recombination (hypothesis I), the asymptotic explosion energy is $E_{infty} approx M_g v_0^2 (1.5eta^{2/3} + beta f(rho_0))$, where $beta$ is the dimensionless recombination parameter. We use Bayesian inference to fit these analytic models to simulations. Both hypotheses fit the simulations of the lowest progenitor masses that tend to explode spherically. The fits do not prefer hypothesis I or II; however, prior investigations suggest that $alpha$ recombination is important. As expected, neither hypothesis fits the higher-mass simulations that exhibit aspherical explosions. In summary, this explosion-energy theory is consistent with the spherical explosions of low progenitor masses; the inconsistency with higher progenitor-mass simulations suggests that a theory for them must include aspherical dynamics.