The ejected mass distribution of Type Ia supernovae: a significant rate of non-Chandrasekhar-mass progenitors

The ejected mass distribution of Type Ia supernovae (SNe Ia) directly probes progenitor evolutionary history and explosion mechanisms, with implications for their use as cosmological probes. Although the Chandrasekhar mass is a natural mass scale for the explosion of white dwarfs as SNe Ia, models allowing SNe Ia to explode at other masses have attracted much recent attention. Using an empirical relation between the ejected mass and the light-curve width, we derive ejected masses Mej and 56Ni masses MNi for a sample of 337 SNe Ia with redshifts z <0.7 used in recent cosmological analyses. We use hierarchical Bayesian inference to reconstruct the joint Mej-MNi distribution, accounting for measurement errors. The inferred marginal distribution of Mej has a long tail towards sub-Chandrasekhar masses, but cuts off sharply above 1.4 M⊙. Our results imply that 25-50 per cent of normal SNe Ia are inconsistent with Chandrasekhar-mass explosions, with almost all of these being sub-Chandrasekhar mass; super-Chandrasekhar-mass explosions make up no more than 1 per cent of all spectroscopically normal SNe Ia. We interpret the SN Ia width-luminosity relation as an underlying relation between Mej and MNi, and show that the inferred relation is not naturally explained by the predictions of any single known explosion mechanism.

The white dwarf's carbon fraction as a secondary parameter of Type Ia supernovae

Context. Binary stellar evolution calculations predict thatChandrasekhar-mass carbon/oxygen white dwarfs (WDs) show a radiallyvarying profile for the composition with a carbon depleted core. Manyrecent multi-dimensional simulations of Type Ia supernovae (SNe Ia),however, assume the progenitor WD has a homogeneous chemicalcomposition.
Aims: In this work, we explore the impact ofdifferent initial carbon profiles of the progenitor WD on the explosionphase and on synthetic observables in the Chandrasekhar-mass delayeddetonation model. Spectra and light curves are compared to observationsto judge the validity of the model.
Methods: The explosion phaseis simulated using the finite volume supernova code Leafs, which isextended to treat different compositions of the progenitor WD. Thesynthetic observables are computed with the Monte Carlo radiativetransfer code Artis. Results: Differences in binding energies ofcarbon and oxygen lead to a lower nuclear energy release for carbondepleted material; thus, the burning fronts that develop are weaker andthe total nuclear energy release is smaller. For otherwise identicalconditions, carbon depleted models produce less 56Ni.Comparing different models with similar 56Ni yields showslower kinetic energies in the ejecta for carbon depleted models, butonly small differences in velocity distributions and line velocities inspectra. The light curve width-luminosity relation (WLR) obtained formodels with differing carbon depletion is roughly perpendicular to theobserved WLR, hence the carbon mass fraction is probably only asecondary parameter in the family of SNe Ia.
Tables 3 and 4 are available in electronic form at http://www.aanda.org

Locations of peculiar supernovae as a diagnostic of their origins

We put constraints on the properties of the progenitors of peculiar calcium-rich transients using the distribution of locations within their host galaxies. We confirm that this class of transients do not follow the galaxy stellar mass profile and are more likely to be found in remote locations of their apparent hosts. We test the hypothesis that these transients are from low-metallicity progenitors by comparing their spatial distributions with the predictions of self-consistent cosmological simulations that include star formation and chemical enrichment. We find that while metal-poor stars and our transient sample show a consistent preference for large offsets, metallicity alone cannot explain the extreme cases. Invoking a lower age limit on the progenitor helps to improve the match, indicating these events may result from a very old metal-poor population. We also investigate the radial distribution of globular cluster systems, and show that they too are consistent with the class of calcium-rich transients. Because photometric upper limits exist for globular clusters for some members of the class, a production mechanism related to the dense environment of globular clusters is not favoured for the calcium-rich events. However, the methods developed in this paper may be used in the future to constrain the effects of low metallicity on radially distant core-collapse events or help establish a correlation with globular clusters for other classes of peculiar explosions.

Gamma-ray diagnostics of Type Ia supernovae. Predictions of observables from three-dimensional modeling

Context. Although the question of progenitor systems and detailed explosion mechanisms still remains a matter of discussion, it is commonly believed that Type Ia supernovae (SNe Ia) are production sites of large amounts of radioactive nuclei. Even though the gamma-ray emission due to radioactive decays is responsible for powering the light curves of SNe Ia, gamma rays themselves are of particular interest as a diagnostic tool because they directly lead to deeper insight into the nucleosynthesis and the kinematics of these explosion events. Aims: We study the evolution of gamma-ray line and continuum emission of SNe Ia with the objective of analyzing the relevance of observations in this energy range. We seek to investigate the chances for the success of future MeV missions regarding their capabilities for constraining the intrinsic properties and the physical processes of SNe Ia. Methods: Focusing on two of the most broadly discussed SN Ia progenitor scenarios - a delayed detonation in a Chandrasekhar-mass white dwarf (WD) and a violent merger of two WDs - we used three-dimensional explosion models and performed radiative transfer simulations to obtain synthetic gamma-ray spectra. Both chosen models produce the same mass of 56Ni and have similar optical properties that are in reasonable agreement with the recently observed supernova SN 2011fe. We examine the gamma-ray spectra with respect to their distinct features and draw connections to certain characteristics of the explosion models. Applying diagnostics, such as line and hardness ratios, the detection prospects for future gamma-ray missions with higher sensitivities in the MeV energy range are discussed. Results: In contrast to the optical regime, the gamma-ray emission of our two chosen models proves to be quite different. The almost direct connection of the emission of gamma rays to fundamental physical processes occurring in SNe Ia permits additional constraints concerning several explosion model properties that are not easily accessible within other wavelength ranges. Proposed future MeV missions such as GRIPS will resolve all spectral details only for nearby SNe Ia, but hardness ratio and light curve measurements still allow for a distinction of the two different models at 10 Mpc and 16 Mpc for an exposure time of 106 s. The possibility of detecting the strongest line features up to the Virgo distance will offer the opportunity to build up a first sample of SN Ia detections in the gamma-ray energy range and underlines the importance of future space observatories for MeV gamma rays.

PTF11iqb: cool supergiant mass-loss that bridges the gap between Type IIn and normal supernovae

PTF11iqb was initially classified as a TypeIIn event caught very early after explosion. It showed narrow Wolf-Rayet (WR) spectral features on day 2, but the narrow emission weakened quickly and the spectrum morphed to resemble those of Types II-L and II-P. At late times, Halpha emission exhibited a complex, multipeaked profile reminiscent of SN1998S. In terms of spectroscopic evolution, we find that PTF11iqb was a near twin of SN~1998S, although with weaker interaction with circumstellar material (CSM) at early times, and stronger CSM interaction at late times. We interpret the spectral changes as caused by early interaction with asymmetric CSM that is quickly (by day 20) enveloped by the expanding SN ejecta photosphere, but then revealed again after the end of the plateau when the photosphere recedes. The light curve can be matched with a simple model for weak CSM interaction added to the light curve of a normal SN~II-P. This plateau requires that the progenitor had an extended H envelope like a red supergiant, consistent with the slow progenitor wind speed indicated by narrow emission. The cool supergiant progenitor is significant because PTF11iqb showed WR features in its early spectrum --- meaning that the presence of such WR features in an early SN spectrum does not necessarily indicate a WR-like progenitor. [abridged] Overall, PTF11iqb bridges SNe~IIn with weaker pre-SN mass loss seen in SNe II-L and II-P, implying a continuum between these types.

Type Ia supernovae from violent mergers of carbon-oxygen white dwarfs: polarization signatures

The violent merger of two carbon-oxygen white dwarfs has been proposed as a viable progenitor for some Type Ia supernovae. However, it has been argued that the strong ejecta asymmetries produced by this model might be inconsistent with the low degree of polarization typically observed in Type Ia supernova explosions. Here, we test this claim by carrying out a spectropolarimetric analysis for the model proposed by Pakmor et al. for an explosion triggered during the merger of a 1.1 and 0.9 M⊙ carbon-oxygen white dwarf binary system. Owing to the asymmetries of the ejecta, the polarization signal varies significantly with viewing angle. We find that polarization levels for observers in the equatorial plane are modest (≲1 per cent) and show clear evidence for a dominant axis, as a consequence of the ejecta symmetry about the orbital plane. In contrast, orientations out of the plane are associated with higher degrees of polarization and departures from a dominant axis. While the particular model studied here gives a good match to highly polarized events such as SN 2004dt, it has difficulties in reproducing the low polarization levels commonly observed in normal Type Ia supernovae. Specifically, we find that significant asymmetries in the element distribution result in a wealth of strong polarization features that are not observed in the majority of currently available spectropolarimetric data of Type Ia supernovae. Future studies will map out the parameter space of the merger scenario to investigate if alternative models can provide better agreement with observations.

Searching for swept-up hydrogen and helium in the late-time spectra of 11 nearby Type Ia supernovae

The direct detection of a stellar system that explodes as a Type Ia supernova (SN Ia) has not yet been successful. Various indirect methods have been used to investigate SN Ia progenitor systems but none have produced conclusive results. A prediction of single-degenerate models is that H- (or He-) rich material from the envelope of the companion star should be swept up by the SN ejecta in the explosion. Seven SNe Ia have been analysed to date looking for signs of H-rich material in their late-time spectra and none were detected. We present results from new late-time spectra of 11 SNe Ia obtained at the Very Large Telescope using XShooter and FORS2. We present the tentative detection of Hα emission for SN 2013ct, corresponding to ∼0.007 M_⊙ of stripped/ablated companion star material (under the assumptions of the spectral modelling). This mass is significantly lower than expected for single-degenerate scenarios, suggesting that >0.1 M_⊙ of H-rich is present but not observed. We do not detect Hα emission in the other 10 SNe Ia. This brings the total sample of normal SNe Ia with non-detections (<0.001–0.058 M_⊙) of H-rich material to 17 events. The simplest explanation for these non-detections is that these objects did not result from the explosion of a CO white dwarf accreting matter from a H-rich companion star via Roche lobe overflow or symbiotic channels. However, further spectral modelling is needed to confirm this. We also find no evidence of He-emission features, but models with He-rich companion stars are not available to place mass limits.

Progenitor-Dependent Explosion Dynamics in Self-Consistent, Axisymmetric Simulations of Neutrino-Driven Core-Collapse Supernovae

We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11–28 M ⊙, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si–O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si–O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection timescales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.

Non-radial instabilities and progenitor asphericities in core-collapse supernovae

Since core-collapse supernova simulations still struggle to produce robust neutrino-driven explosions in 3D, it has been proposed that asphericities caused by convection in the progenitor might facilitate shock revival by boosting the activity of non-radial hydrodynamic instabilities in the post-shock region. We investigate this scenario in depth using 42 relativistic 2D simulations with multigroup neutrino transport to examine the effects of velocity and density perturbations in the progenitor for different perturbation geometries that obey fundamental physical constraints (like the anelastic condition). As a framework for analysing our results, we introduce semi-empirical scaling laws relating neutrino heating, average turbulent velocities in the gain region, and the shock deformation in the saturation limit of non-radial instabilities. The squared turbulent Mach number, 〈Ma2〉, reflects the violence of aspherical motions in the gain layer, and explosive runaway occurs for 〈Ma2〉 ≳ 0.3, corresponding to a reduction of the critical neutrino luminosity by ∼25∼25 per cent compared to 1D. In the light of this theory, progenitor asphericities aid shock revival mainly by creating anisotropic mass flux on to the shock: differential infall efficiently converts velocity perturbations in the progenitor into density perturbations δρ/ρ at the shock of the order of the initial convective Mach number Maprog. The anisotropic mass flux and ram pressure deform the shock and thereby amplify post-shock turbulence. Large-scale (ℓ = 2, ℓ = 1) modes prove most conducive to shock revival, whereas small-scale perturbations require unrealistically high convective Mach numbers. Initial density perturbations in the progenitor are only of the order of Ma2progMaprog2 and therefore play a subdominant role.

Applying the expanding photosphere and standardized candle methods to Type II-Plateau supernovae at cosmologically significant redshifts . The distance to SN 2013eq

Based on optical imaging and spectroscopy of the Type II-Plateau SN 2013eq, we present a comparative study of commonly used distance determination methods based on Type II supernovae. The occurrence of SN 2013eq in the Hubble flow (z = 0.041 ± 0.001) prompted us to investigate the implications of the difference between "angular" and "luminosity" distances within the framework of the expanding photosphere method (EPM) that relies upon a relation between flux and angular size to yield a distance. Following a re-derivation of the basic equations of the EPM for SNe at non-negligible redshifts, we conclude that the EPM results in an angular distance. The observed flux should be converted into the SN rest frame and the angular size, θ, has to be corrected by a factor of (1 + z)2. Alternatively, the EPM angular distance can be converted to a luminosity distance by implementing a modification of the angular size. For SN 2013eq, we find EPM luminosity distances of DL = 151 ± 18 Mpc and DL = 164 ± 20 Mpc by making use of different sets of dilution factors taken from the literature. Application of the standardized candle method for Type II-P SNe results in an independent luminosity distance estimate (DL = 168 ± 16 Mpc) that is consistent with the EPM estimate. Spectra of SN 2013eq are available in the Weizmann Interactive Supernova data REPository (WISeREP): http://wiserep.weizmann.ac.il

Predicting polarization signatures for double-detonation and delayed-detonation models of Type Ia supernovae

Calculations of synthetic spectropolarimetry are one means to test multidimensional explosion models for Type Ia supernovae. In a recent paper, we demonstrated that the violent merger of a 1.1 and 0.9 M⊙ white dwarf binary system is too asymmetric to explain the low polarization levels commonly observed in normal Type Ia supernovae. Here, we present polarization simulations for two alternative scenarios: the sub-Chandrasekhar mass double-detonation and the Chandrasekhar mass delayed-detonation model. Specifically, we study a 2D double-detonation model and a 3D delayed-detonation model, and calculate polarization spectra for multiple observer orientations in both cases. We find modest polarization levels (<1 per cent) for both explosion models. Polarization in the continuum peaks at ∼0.1–0.3 per cent and decreases after maximum light, in excellent agreement with spectropolarimetric data of normal Type Ia supernovae. Higher degrees of polarization are found across individual spectral lines. In particular, the synthetic Si II λ6355 profiles are polarized at levels that match remarkably well the values observed in normal Type Ia supernovae, while the low degrees of polarization predicted across the O I λ7774 region are consistent with the non-detection of this feature in current data. We conclude that our models can reproduce many of the characteristics of both flux and polarization spectra for well-studied Type Ia supernovae, such as SN 2001el and SN 2012fr. However, the two models considered here cannot account for the unusually high level of polarization observed in extreme cases such as SN 2004dt.

A study of the Type II-P supernova 2003gd in M74

We present photometric and spectroscopic data of the Type II-P supernova (SN II-P) 2003gd, which was discovered in M74 close to the end of its plateau phase. SN 2003gd is the first Type II supernova ( SN) to have a directly confirmed red supergiant ( RSG) progenitor. We compare SN 2003gd to SN 1999 em, a similar SN II-P, and estimate an explosion date of 2003 March 18. We determine a reddening towards the SN of E(B-V) = 0.14 +/- 0.06, using three different methods. We also calculate three new distances to M74 of 9.6 +/- 2.8, 7.7 +/- 1.7 and 9.6 +/- 2.2 Mpc. The former was estimated using the standard candle method (SCM), for Type II supernovae (SNe II), and the latter two using the brightest supergiants method (BSM). When combined with existing kinematic and BSM distance estimates, we derive a mean value of 9.3 +/- 1.8 Mpc. SN 2003gd was found to have a lower tail luminosity compared with other normal Type II-P supernovae ( SNe II-P) bringing into question the nature of this SN. We present a discussion concluding that this is a normal SN II-P, which is consistent with the observed progenitor mass of 8(-2)(+4) M-circle dot.

Detection of a Red Supergiant Progenitor Star of a Type II-Plateau Supernova

We present the discovery of a red supergiant star that exploded as supernova 2003gd in the nearby spiral galaxy M74. The Hubble Space Telescope (HST) and the Gemini Telescope imaged this galaxy 6 to 9 months before the supernova explosion, and subsequent HST images confirm the positional coincidence of the supernova with a single resolved star that is a red supergiant of 8+4-2 solar masses. This confirms both stellar evolution models and supernova theories predicting that cool red supergiants are the immediate progenitor stars of type II-plateau supernovae.

The massive binary companion star to the progenitor of supernova 1993J

The massive star that underwent a collapse of its core to produce supernova (SN)1993J was subsequently identified as a non-variable red supergiant star in images of the galaxy M81 taken before explosion(1, 2). It showed an excess in ultraviolet and B-band colours, suggesting either the presence of a hot, massive companion star or that it was embedded in an unresolved young stellar association1. The spectra of SN1993J underwent a remarkable transformation from the signature of a hydrogen-rich type II supernova to one of a helium-rich (hydrogen-deficient) type Ib(3, 4). The spectral and photometric peculiarities were best explained by models in which the 13�20 solar mass supergiant had lost almost its entire hydrogen envelope to a close binary companion(5, 6, 7), producing a 'type IIb' supernova, but the hypothetical massive companion stars for this class of supernovae have so far eluded discovery. Here we report photometric and spectroscopic observations of SN1993J ten years after the explosion. At the position of the fading supernova we detect the unambiguous signature of a massive star: the binary companion to the progenitor.