Synthetic light curves and spectra for three-dimensional delayed-detonation models of type ia supernovae

In a companion paper, Seitenzahl et al. have presented a set of three-dimensional delayed detonation models for thermonuclear explosions of near-Chandrasekhar-mass white dwarfs (WDs). Here,we present multidimensional radiative transfer simulations that provide synthetic light curves and spectra for those models. The model sequence explores both changes in the strength of the deflagration phase (which is controlled by the ignition configuration in our models) and the WD central density. In agreement with previous studies, we find that the strength of the deflagration significantly affects the explosion and the observables. Variations in the central density also have an influence on both brightness and colour, but overall it is a secondary parameter in our set of models. In many respects, the models yield a good match to the observed properties of normal Type Ia supernovae (SNe Ia): peak brightness, rise/decline time-scales and synthetic spectra are all in reasonable agreement. There are, however, several differences. In particular, the models are systematically too red around maximum light, manifest spectral line velocities that are a little too high and yield I-band light curves that do not match observations. Although some of these discrepancies may simply relate to approximations made in the modelling, some pose real challenges to the models. If viewed as a complete sequence, our models do not reproduce the observed light-curve width- luminosity relation (WLR) of SNe Ia: all our models show rather similar B-band decline rates, irrespective of peak brightness. This suggests that simple variations in the strength of the deflagration phase in Chandrasekhar-mass deflagration-to-detonation models do not readily explain the observed diversity of normal SNe Ia. This may imply that some other parameter within the Chandrasekhar-mass paradigm is key to the WLR, or that a substantial fraction of normal SNe Ia arise from an alternative explosion scenario.

Multi-Dimensional Simulations of the Expanding Supernova Remnant of SN 1987A

The expanding remnant from SN 1987A is an excellent laboratory for investigating the physics of supernovae explosions. There is still a large number of outstanding questions, such as the reason for the asymmetric radio morphology, the structure of the pre-supernova environment, and the efficiency of particle acceleration at the supernova shock. We explore these questions using three-dimensional simulations of the expanding remnant between days 820 and 10,000 after the supernova. We combine a hydrodynamical simulation with semi-analytic treatments of diffusive shock acceleration and magnetic field amplification to derive radio emission as part of an inverse problem. Simulations show that an asymmetric explosion, combined with magnetic field amplification at the expanding shock, is able to replicate the persistent one-sided radio morphology of the remnant. We use an asymmetric Truelove & McKee progenitor with an envelope mass of 10 M-circle dot and an energy of 1.5 x 10(44) J. A termination shock in the progenitor's stellar wind at a distance of 0 ''.43-0 ''.51 provides a good fit to the turn on of radio emission around day 1200. For the H II region, a minimum distance of 0 ''.63 +/- 0 ''.01 and maximum particle number density of (7.11 +/- 1.78) x 10(7) m(-3) produces a good fit to the evolving average radius and velocity of the expanding shocks from day 2000 to day 7000 after explosion. The model predicts a noticeable reduction, and possibly a temporary reversal, in the asymmetric radio morphology of the remnant after day 7000, when the forward shock left the eastern lobe of the equatorial ring.

Toward Characterization Of The Type IIP Supernova Progenitor Population: A Statistical Sample Of Light Curves From Pan-STARRS1

In recent years, wide-field sky surveys providing deep multi-band imaging have presented a new path for indirectly characterizing the progenitor populations of core-collapse supernovae (SN): systematic light curve studies. We assemble a set of 76 grizy-band Type IIP SN light curves from Pan-STARRS1, obtained over a constant survey program of 4 years and classified using both spectroscopy and machine learning-based photometric techniques. We develop and apply a new Bayesian model for the full multi-band evolution of each light curve in the sample. We find no evidence of a sub-population of fast-declining explosions (historically referred to as "Type IIL" SNe). However, we identify a highly significant relation between the plateau phase decay rate and peak luminosity among our SNe IIP. These results argue in favor of a single parameter, likely determined by initial stellar mass, predominantly controlling the explosions of red supergiants. This relation could also be applied for supernova cosmology, offering a standardizable candle good to an intrinsic scatter of 0.2 mag. We compare each light curve to physical models from hydrodynamic simulations to estimate progenitor initial masses and other properties of the Pan-STARRS1 Type IIP SN sample. We show that correction of systematic discrepancies between modeled and observed SN IIP light curve properties and an expanded grid of progenitor properties, are needed to enable robust progenitor inferences from multi-band light curve samples of this kind. This work will serve as a pathfinder for photometric studies of core-collapse SNe to be conducted through future wide field transient searches.

Supersolar Ni/Fe production in the Type IIP SN 2012ec

SN 2012ec is a Type IIP supernova (SN) with a progenitor detection and comprehensive photospheric phase observational coverage. Here, we present Very Large Telescope and Public ESO Spectroscopic Survey of Transient Objects observations of this SN in the nebular phase. We model the nebular [O I] lambda lambda 6300, 6364 lines and find their strength to suggest a progenitor main-sequence mass of 13-15 M-circle dot. SN2012ec is unique among hydrogen-rich SNe in showing a distinct line of stable nickel [Ni II] lambda 7378. This line is produced by Ni-58, a nuclear burning ash whose abundance is a sensitive tracer of explosive burning conditions. Using spectral synthesis modelling, we use the relative strengths of [Ni II] lambda 7378 and [Fe II] lambda 7155 (the progenitor of which is Ni-56) to derive a Ni/Fe production ratio of 0.20 +/- 0.07 (by mass), which is a factor 3.4 +/- 1.2 times the solar value. High production of stable nickel is confirmed by a strong [Ni II] 1.939 mu m line. This is the third reported case of a core-collapse SN producing a Ni/Fe ratio far above the solar value, which has implications for core-collapse explosion theory and galactic chemical evolution models.

Observational Constraints on the Progenitors of Core-Collapse Supernovae: The Case for Missing High-Mass Stars

Over the last 15 years, the supernova community has endeavoured to directly identify progenitor stars for core-collapse supernovae discovered in nearby galaxies. These precursors are often visible as resolved stars in high-resolution images from space-and ground-based telescopes. The discovery rate of progenitor stars is limited by the local supernova rate and the availability and depth of archive images of galaxies, with 18 detections of precursor objects and 27 upper limits. This review compiles these results (from 1999 to 2013) in a distance-limited sample and discusses the implications of the findings. The vast majority of the detections of progenitor stars are of type II-P, II-L, or IIb with one type Ib progenitor system detected and many more upper limits for progenitors of Ibc supernovae (14 in all). The data for these 45 supernovae progenitors illustrate a remarkable deficit of high-luminosity stars above an apparent limit of log L/L-circle dot similar or equal to 5.1 dex. For a typical Salpeter initial mass function, one would expect to have found 13 high-luminosity and high-mass progenitors by now. There is, possibly, only one object in this time-and volume-limited sample that is unambiguously high-mass (the progenitor of SN2009ip) although the nature of that supernovae is still debated. The possible biases due to the influence of circumstellar dust, the luminosity analysis, and sample selection methods are reviewed. It does not appear likely that these can explain the missing high-mass progenitor stars. This review concludes that the community's work to date shows that the observed populations of supernovae in the local Universe are not, on the whole, produced by high-mass (M greater than or similar to 18 M-circle dot) stars. Theoretical explosions of model stars also predict that black hole formation and failed supernovae tend to occur above an initial mass of M similar or equal to 18 M-circle dot. The models also suggest there is no simple single mass division for neutron star or black-hole formation and that there are islands of explodability for stars in the 8-120 M-circle dot range. The observational constraints are quite consistent with the bulk of stars above M similar or equal to 18 M-circle dot collapsing to form black holes with no visible supernovae.

Massive stars exploding in a He-rich circumstellar medium - VI. Observations of two distant Type Ibn supernova candidates discovered by La Silla-QUEST

We present optical observations of the peculiar stripped-envelope supernovae (SNe) LSQ12btw and LSQ13ccw discovered by the La Silla-QUEST survey. LSQ12btw reaches an absolute peak magnitude of M-g = -19.3 +/- 0.2, and shows an asymmetric light curve. Stringent pre-discovery limits constrain its rise time to maximum light to less than 4 d, with a slower post-peak luminosity decline, similar to that experienced by the prototypical SN Ibn 2006jc. LSQ13ccw is somewhat different: while it also exhibits a very fast rise to maximum, it reaches a fainter absolute peak magnitude (M-g =-18.4 +/- 0.2), and experiences an extremely rapid post-peak decline similar to that observed in the peculiar SN Ib 2002bj. A stringent pre-discovery limit and an early marginal detection of LSQ13ccw allow us to determine the explosion time with an uncertainty of +/- 1 d. The spectra of LSQ12btw show the typical narrow He I emission lines characterizing Type Ibn SNe, suggesting that the SN ejecta are interacting with He-rich circumstellar material. The He I lines in the spectra of LSQ13ccw exhibit weak narrow emissions superposed on broad components. An unresolved H alpha line is also detected, suggesting a tentative Type Ibn/IIn classification. As for other SNe Ibn, we argue that LSQ12btw and LSQ13ccw likely result from the explosions of Wolf-Rayet stars that experienced instability phases prior to core collapse. We inspect the host galaxies of SNe Ibn, and we show that all of them but one are hosted in spiral galaxies, likely in environments spanning a wide metallicity range.

Deflagrations in hybrid CONe white dwarfs: a route to explain the faint Type Iax supernova 2008ha

Stellar evolution models predict the existence of hybrid white dwarfs (WDs) with a carbon-oxygen core surrounded by an oxygen-neon mantle. Being born with masses similar to 1.1 M-aS (TM), hybrid WDs in a binary system may easily approach the Chandrasekhar mass (M-Ch) by accretion and give rise to a thermonuclear explosion. Here, we investigate an off-centre deflagration in a near-M-Ch hybrid WD under the assumption that nuclear burning only occurs in carbon-rich material. Performing hydrodynamics simulations of the explosion and detailed nucleosynthesis post-processing calculations, we find that only 0.014 M-aS (TM) of material is ejected while the remainder of the mass stays bound. The ejecta consist predominantly of iron-group elements, O, C, Si and S. We also calculate synthetic observables for our model and find reasonable agreement with the faint Type Iax SN 2008ha. This shows for the first time that deflagrations in near-M-Ch WDs can in principle explain the observed diversity of Type Iax supernovae. Leaving behind a near-M-Ch bound remnant opens the possibility for recurrent explosions or a subsequent accretion-induced collapse in faint Type Iax SNe, if further accretion episodes occur. From binary population synthesis calculations, we find the rate of hybrid WDs approaching M-Ch to be of the order of 1 per cent of the Galactic SN Ia rate.

Hydrogen-Poor Superluminous Supernovae and Long-Duration Gamma-Ray Bursts Have Similar Host Galaxies

We present optical spectroscopy and optical/near-IR photometry of 31 host galaxies of hydrogen-poor superluminous supernovae (SLSNe), including 15 events from the Pan-STARRS1 Medium Deep Survey. Our sample spans the redshift range 0.1 ≲ z ≲ 1.6, and is the first comprehensive host galaxy study of this specific subclass of cosmic explosions. Combining the multi-band photometry and emission-line measurements, we determine the luminosities, stellar masses, star formation rates, and metallicities. We find that, as a whole, the hosts of SLSNe are a low-luminosity (〈M_B 〉 ≈ -17.3 mag), low stellar mass (〈M〉 ≈ 2 × 10⁸ M_⊙) population, with a high median specific star formation rate (〈sSFR〉 ≈ 2 Gyr^-1). The median metallicity of our spectroscopic sample is low, 12 + log (O/H) ≈ 8.35 ≈ 0.45 Z_⊙, although at least one host galaxy has solar metallicity. The host galaxies of H-poor SLSNe are statistically distinct from the hosts of GOODS core-collapse SNe (which cover a similar redshift range), but resemble the host galaxies of long-duration gamma-ray bursts (LGRBs) in terms of stellar mass, SFR, sSFR, and metallicity. This result indicates that the environmental causes leading to massive stars forming either SLSNe or LGRBs are similar, and in particular that SLSNe are more effectively formed in low metallicity environments. We speculate that the key ingredient is large core angular momentum, leading to a rapidly spinning magnetar in SLSNe and an accreting black hole in LGRBs.

Constraining the physical properties of Type II-Plateau supernovae using nebular phase spectra

We present a study of the nebular phase spectra of a sample of Type II-Plateau supernovae with identified progenitors or restrictive limits. The evolution of line fluxes, shapes and velocities is compared within the sample, and interpreted by the use of a spectral synthesis code. The small diversity within the data set can be explained by strong mixing occurring during the explosion, and by recognizing that most lines have significant contributions from primordial metals in the H envelope, which dominates the total ejecta mass in these types of objects. In particular, when using the [O I] 6300, 6364 Å doublet for estimating the core mass of the star, care has to be taken to account for emission from primordial O in the envelope. Finally, a correlation between the Hα line width and the mass of 56Ni is presented, suggesting that higher energy explosions are associated with higher 56Ni production.

The host galaxy and late-time evolution of the superluminous supernova PTF12dam

Superluminous supernovae (SLSNe) of Type Ic have a tendency to occur in faint host galaxies which are likely to have low mass and low metallicity. PTF12dam is one of the closest and best-studied superluminous explosions that has a broad and slowly fading light curve similar to SN 2007bi. Here we present new photometry and spectroscopy for PTF12dam from 200-500 d (rest frame) after peak and a detailed analysis of the host galaxy (SDSS J142446.21+461348.6 at z = 0.107). Using deep templates and image subtraction we show that the light curve can be fit with a magnetar model if escape of high-energy gamma rays is taken into account. The full bolometric light curve from -53 to +399 d (with respect to peak) cannot be fit satisfactorily with the pair-instability models. An alternative model of interaction with a dense circumstellar material (CSM) produces a good fit to the data although this requires a very large mass (˜13 M⊙) of hydrogen-free CSM. The host galaxy is a compact dwarf (physical size ˜1.9 kpc) and with Mg = -19.33 ± 0.10, it is the brightest nearby SLSN Ic host discovered so far. The host is a low-mass system (2.8 × 108 M⊙) with a star formation rate (5.0 M⊙ yr-1), which implies a very high specific star formation rate (17.9 Gyr-1). The remarkably strong nebular emission provide detections of the [O III] λ4363 and [O II] λλ7320, 7330auroral lines and an accurate oxygen abundance of 12 + log (O/H) = 8.05 ± 0.09. We show here that they are at the extreme end of the metallicity distribution of dwarf galaxies and propose that low metallicity is a requirement to produce these rare and peculiar SNe.

Type Ia supernovae from exploding oxygen-neon white dwarfs

Context. The progenitor problem of Type Ia supernovae (SNe Ia) is still unsolved. Most of these events are thought to be explosions of carbon-oxygen (CO) white dwarfs (WDs), but for many of the explosion scenarios, particularly those involving the externally triggered detonation of a sub-Chandrasekhar mass WD (sub-M-Ch, WD), there is also a possibility of having an oxygen-neon (ONe) WD as progenitor.

Aims. We simulate detonations of ONe WDs and calculate synthetic observables from these models. The results are compared with detonations in CO WDs of similar mass and observational data of SNe Ia.

Methods. We perform hydrodynamic explosion simulations of detonations in initially hydrostatic ONe WDs for a range of masses below the Chandrasekhar mass (M-Ch), followed by detailed nucleosynthetic postprocessing with a 384-isotope nuclear reaction network. The results are used to calculate synthetic spectra and light curves, which are then compared with observations of SNe Ia. We also perform binary evolution calculations to determine the number of SNe Ia involving ONe WDs relative to the number of other promising progenitor channels.

Results. The ejecta structures of our simulated detonations in sub-M-Ch, ONe WDs are similar to those from CO WDs. There are, however, small systematic deviations in the mass fractions and the ejecta velocities. These lead to spectral features that are systematically less blueshifted. Nevertheless, the synthetic observables of our ONe WD explosions are similar to those obtained from CO models.

Conclusions. Our binary evolution calculations show that a significant fraction (3-10%) of potential progenitor systems should contain an ONe WD. The comparison of our ONe models with our CO models of comparable mass (similar to 1.2 M-circle dot) shows that the less blueshifted spectral features fit the observations better, although they are too bright for normal SNe Ia.

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.

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.

A comparative study of Type II-P and II-L supernova rise times as exemplified by the case of LSQ13cuw

We report on our findings based on the analysis of observations of the Type II-L supernova LSQ13cuw within the framework of currently accepted physical predictions of core-collapse supernova explosions. LSQ13cuw was discovered within a day of explosion, hitherto unprecedented for Type II-L supernovae. This motivated a comparative study of Type II-P and II-L supernovae with relatively well-constrained explosion epochs and rise times to maximum (optical) light. From our sample of twenty such events, we find evidence of a positive correlation between the duration of the rise and the peak brightness. On average, SNe II-L tend to have brighter peak magnitudes and longer rise times than SNe II-P. However, this difference is clearest only at the extreme ends of the rise time versus peak brightness relation. Using two different analytical models, we performed a parameter study to investigate the physical parameters that control the rise time behaviour. In general, the models qualitatively reproduce aspects of the observed trends. We find that the brightness of the optical peak increases for larger progenitor radii and explosion energies, and decreases for larger masses. The dependence of the rise time on mass and explosion energy is smaller than the dependence on the progenitor radius. We find no evidence that the progenitors of SNe II-L have significantly smaller radii than those of SNe II-P.

SN 2011dh: DISCOVERY OF A TYPE IIb SUPERNOVA FROM A COMPACT PROGENITOR IN THE NEARBY GALAXY M51

On 2011 May 31 UT a supernova (SN) exploded in the nearby galaxy M51 (the Whirlpool Galaxy). We discovered this event using small telescopes equipped with CCD cameras and also detected it with the Palomar Transient Factory survey, rapidly confirming it to be a Type II SN. Here, we present multi-color ultraviolet through infrared photometry which is used to calculate the bolometric luminosity and a series of spectra. Our early-time observations indicate that SN 2011dh resulted from the explosion of a relatively compact progenitor star. Rapid shock-breakout cooling leads to relatively low temperatures in early-time spectra, compared to explosions of red supergiant stars, as well as a rapid early light curve decline. Optical spectra of SN 2011dh are dominated by H lines out to day 10 after explosion, after which He I lines develop. This SN is likely a member of the cIIb (compact IIb) class, with progenitor radius larger than that of SN 2008ax and smaller than the eIIb (extended IIb) SN 1993J progenitor. Our data imply that the object identified in pre-explosion Hubble Space Telescope images at the SN location is possibly a companion to the progenitor or a blended source, and not the progenitor star itself, as its radius (~1013 cm) would be highly inconsistent with constraints from our post-explosion spectra.