Three-dimensional delayed-detonation models with nucleosynthesis for type ia supernovae

We present results for a suite of 14 three-dimensional, high-resolution hydrodynamical simulations of delayed-detonation models of Type Ia supernova (SN Ia) explosions. This model suite comprises the first set of three-dimensional SN Ia simulations with detailed isotopic yield information. As such, it may serve as a data base for Chandrasekhar-mass delayed-detonation model nucleosynthetic yields and for deriving synthetic observables such as spectra and light curves. We employ aphysically motivated, stochastic model based on turbulent velocity fluctuations and fuel density to calculate in situ the deflagration-to-detonation transition probabilities. To obtain different strengths of the deflagration phase and thereby different degrees of pre-expansion, we have chosen a sequence of initial models with 1, 3, 5, 10, 20, 40, 100, 150, 200, 300 and 1600 (two different realizations) ignition kernels in a hydrostatic white dwarf with a central density of 2.9 × 10 g cm, as well as one high central density (5.5 × 10 g cm) and one low central density (1.0 × 10 g cm) rendition of the 100 ignition kernel configuration. For each simulation, we determined detailed nucleosynthetic yields by postprocessing10 tracer particles with a 384 nuclide reaction network. All delayed-detonation models result in explosions unbinding thewhite dwarf, producing a range of 56Ni masses from 0.32 to 1.11M. As a general trend, the models predict that the stableneutron-rich iron-group isotopes are not found at the lowest velocities, but rather at intermediate velocities (~3000×10 000 km s) in a shell surrounding a Ni-rich core. The models further predict relatively low-velocity oxygen and carbon, with typical minimum velocities around 4000 and 10 000 km s, respectively. © 2012 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.

Sub-luminous type Ia supernovae from the mergers of equal-mass white dwarfs with mass 0.9M

Type Ia supernovae are thought to result from thermonuclear explosions of carbong'oxygen white dwarf stars. Existing models generally explain the observed properties, with the exception of the sub-luminous 1991bg-like supernovae. It has long been suspected that the merger of two white dwarfs could give rise to a type Ia event, but hitherto simulations have failed to produce an explosion. Here we report a simulation of the merger of two equal-mass white dwarfs that leads to a sub-luminous explosion, although at the expense of requiring a single common-envelope phase, and component masses of 0.9M. The light curve is too broad, but the synthesized spectra, red colour and low expansion velocities are all close to what is observed for sub-luminous 1991bg-like events. Although the mass ratios can be slightly less than one and still produce a sub-luminous event, the masses have to be in the range 0.83M to 0.9M. © 2010 Macmillan Publishers Limited. All rights reserved.

Nebular emission-line profiles of Type Ib/c supernovae - Probing the ejecta asphericity

In order to assess qualitatively the ejecta geometry of stripped-envelope core-collapse supernovae (SNe), we investigate 98 late-time spectra of 39 objects, many of them previously unpublished. We perform a Gauss-fitting of the [O ] ??6300, 6364 feature in all spectra, with the position, full width at half maximum and intensity of the ?6300 Gaussian as free parameters, and the ?6364 Gaussian added appropriately to account for the doublet nature of the [O ] feature. On the basis of the best-fitting parameters, the objects are organized into morphological classes, and we conclude that at least half of all Type Ib/c SNe must be aspherical. Bipolar jet models do not seem to be universally applicable, as we find too few symmetric double-peaked [O ] profiles. In some objects, the [O ] line exhibits a variety of shifted secondary peaks or shoulders, interpreted as blobs of matter ejected at high velocity and possibly accompanied by neutron-star kicks to assure momentum conservation. At phases earlier than ~200 d, a systematic blueshift of the [O ] ??6300, 6364 line centroids can be discerned. Residual opacity provides the most convincing explanation of this phenomenon, photons emitted on the rear side of the SN being scattered or absorbed on their way through the ejecta. Once modified to account for the doublet nature of the oxygen feature, the profile of Mg i] ?4571 at sufficiently late phases generally resembles that of [O ] ??6300, 6364, suggesting negligible contamination from other lines and confirming that O and Mg are similarly distributed within the ejecta. © 2009 RAS.

Super-luminous Type Ic Supernovae: Catching a Magnetar by the Tail

We report extensive observational data for five of the lowest redshift Super-Luminous Type Ic Supernovae (SL-SNe Ic) discovered to date, namely, PTF10hgi, SN2011ke, PTF11rks, SN2011kf, and SN2012il. Photometric imaging of the transients at +50 to +230 days after peak combined with host galaxy subtraction reveals a luminous tail phase for four of these SL-SNe. A high-resolution, optical, and near-infrared spectrum from xshooter provides detection of a broad He I ?10830 emission line in the spectrum (+50 days) of SN2012il, revealing that at least some SL-SNe Ic are not completely helium-free. At first sight, the tail luminosity decline rates that we measure are consistent with the radioactive decay of 56Co, and would require 1-4 M ? of 56Ni to produce the luminosity. These 56Ni masses cannot be made consistent with the short diffusion times at peak, and indeed are insufficient to power the peak luminosity. We instead favor energy deposition by newborn magnetars as the power source for these objects. A semi-analytical diffusion model with energy input from the spin-down of a magnetar reproduces the extensive light curve data well. The model predictions of ejecta velocities and temperatures which are required are in reasonable agreement with those determined from our observations. We derive magnetar energies of 0.4 <~ E(1051 erg) lsim 6.9 and ejecta masses of 2.3 <~ M ej(M ?) lsim 8.6. The sample of five SL-SNe Ic presented here, combined with SN 2010gx—the best sampled SL-SNe Ic so far—points toward an explosion driven by a magnetar as a viable explanation for all SL-SNe Ic.

Slowly fading super-luminous supernovae that are not pair-instability explosions

Super-luminous supernovae that radiate more than 10 44 ergs per second at their peak luminosity have recently been discovered in faint galaxies at redshifts of 0.1-4. Some evolve slowly, resembling models of 'pair-instability' supernovae. Such models involve stars with original masses 140-260 times that of the Sun that now have carbon-oxygen cores of 65-130 solar masses. In these stars, the photons that prevent gravitational collapse are converted to electron-positron pairs, causing rapid contraction and thermonuclear explosions. Many solar masses of 56 Ni are synthesized; this isotope decays to 56 Fe via 56 Co, powering bright light curves. Such massive progenitors are expected to have formed from metal-poor gas in the early Universe. Recently, supernova 2007bi in a galaxy at redshift 0.127 (about 12 billion years after the Big Bang) with a metallicity one-third that of the Sun was observed to look like a fading pair-instability supernova. Here we report observations of two slow-to-fade super-luminous supernovae that show relatively fast rise times and blue colours, which are incompatible with pair-instability models. Their late-time light-curve and spectral similarities to supernova 2007bi call the nature of that event into question. Our early spectra closely resemble typical fast-declining super-luminous supernovae, which are not powered by radioactivity. Modelling our observations with 10-16 solar masses of magnetar-energized ejecta demonstrates the possibility of a common explosion mechanism. The lack of unambiguous nearby pair-instability events suggests that their local rate of occurrence is less than 6 × 10 -6 times that of the core-collapse rate. © 2013 Macmillan Publishers Limited. All rights reserved.

A statistical analysis of circumstellar material in Type Ia supernovae

A key tracer of the elusive progenitor systems of Type Ia supernovae (SNe Ia) is the detection of narrow blueshifted time-varying Na I D absorption lines, interpreted as evidence of circumstellar material surrounding the progenitor system. The origin of this material is controversial, but the simplest explanation is that it results from previous mass-loss in a system containing a white dwarf and a non-degenerate companion star. We present new single-epoch intermediate-resolution spectra of 17 low-redshift SNe Ia taken with XShooter on the European Southern Observatory Very Large Telescope. Combining this sample with events from the literature, we confirm an excess (∼20 per cent) of SNe Ia displaying blueshifted narrow Na I D absorption features compared to redshifted Na I D features. The host galaxies of SNe Ia displaying blueshifted absorption profiles are skewed towards later-type galaxies, compared to SNe Ia that show no Na I D absorption and SNe Ia displaying blueshifted narrow Na I D absorption features have broader light curves. The strength of the Na I D absorption is stronger in SNe Ia displaying blueshifted Na I D absorption features than those without blueshifted features, and the strength of the blueshifted Na I D is correlated with the B − V colour of the SN at maximum light. This strongly suggests the absorbing material is local to the SN. In the context of the progenitor systems of SNe Ia, we discuss the significance of these findings and other recent observational evidence on the nature of SN Ia progenitors. We present a summary that suggests that there are at least two distinct populations of normal, cosmologically useful SNe Ia.

Three-dimensional pure deflagration models with nucleosynthesis and synthetic observables for type ia supernovae

We investigate whether pure deflagration models ofChandrasekhar-mass carbon-oxygen white dwarf stars can account for one or more subclass of the observed population of Type Ia supernova (SN Ia) explosions. We compute a set of 3D full-star hydrodynamic explosion models, in which the deflagration strength is parametrized using the multispot ignition approach. For each model, we calculate detailed nucleosynthesis yields in a post-processing step with a 384 nuclide nuclear network. We also compute synthetic observables with our 3D Monte Carlo radiative transfer code for comparison with observations. For weak and intermediate deflagration strengths (energy release E {less-than or approximate} 1.1 × 10 erg), we find that the explosion leaves behind a bound remnant enriched with 3 to 10 per cent (by mass) of deflagration ashes. However, we do not obtain the large kick velocities recently reported in the literature. We find that weak deflagrations with E ~ 0.5 × 10 erg fit well both the light curves and spectra of 2002cx-like SNe Ia, and models with even lower explosion energies could explain some of the fainter members of this subclass. By comparing our synthetic observables with the properties of SNe Ia, we can exclude the brightest, most vigorously ignited models as candidates for any observed class of SN Ia: their B-V colours deviate significantly from both normal and 2002cx-like SNe Ia and they are too bright to be candidates for other subclasses.

A spectral synthesis code for rapid modelling of supernovae

We present TARDIS-an open-source code for rapid spectral modelling of supernovae (SNe). Our goal is to develop a tool that is sufficiently fast to allow exploration of the complex parameter spaces of models for SN ejecta. This can be used to analyse the growing number of highquality SN spectra being obtained by transient surveys. The code uses Monte Carlo methods to obtain a self-consistent description of the plasma state and to compute a synthetic spectrum. It has a modular design to facilitate the implementation of a range of physical approximations that can be compared to assess both accuracy and computational expediency. This will allow users to choose a level of sophistication appropriate for their application. Here, we describe the operation of the code and make comparisons with alternative radiative transfer codes of differing levels of complexity (SYN++, PYTHON and ARTIS). We then explore the consequence of adopting simple prescriptions for the calculation of atomic excitation, focusing on four species of relevance to Type Ia SN spectra-Si II, SII, MgII and Ca II. We also investigate the influence of three methods for treating line interactions on our synthetic spectra and the need for accurate radiative rate estimates in our scheme.

The effect of helium accretion efficiency on rates of Type Ia supernovae: Double detonations in accreting binaries

he double-detonation explosion scenario of Type Ia supernovae (SNe Ia) has gained increased support from the SN Ia community as a viable progenitor model, making it a promising candidate alongside the well-known single degenerate and double degenerate scenarios. We present delay times of double-detonation SNe, in which a sub-Chandrasekhar mass carbon–oxygen white dwarf (WD) accretes non-dynamically from a helium-rich companion. One of the main uncertainties in quantifying SN rates from double detonations is the (assumed) retention efficiency of He-rich matter. Therefore, we implement a new prescription for the treatment of accretion/accumulation of He-rich matter on WDs. In addition, we test how the results change depending on which criteria are assumed to lead to a detonation in the helium shell. In comparing the results to our standard case (Ruiter et al.), we find that regardless of the adopted He accretion prescription, the SN rates are reduced by only ∼25 per cent if low-mass He shells (≲0.05 M⊙) are sufficient to trigger the detonations. If more massive (0.1 M⊙) shells are needed, the rates decrease by 85 per cent and the delay time distribution is significantly changed in the new accretion model – only SNe with prompt (<500 Myr) delay times are produced. Since theoretical arguments favour low-mass He shells for normal double-detonation SNe, we conclude that the rates from double detonations are likely to be high, and should not critically depend on the adopted prescription for accretion of He.

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.

'Super-Chandrasekhar' Type Ia Supernovae at nebular epochs

We present a first systematic comparison of superluminous Type Ia supernovae (SNe Ia) at late epochs, including previously unpublished photometric and spectroscopic observations of SN 2007if, SN 2009dc and SNF20080723-012. Photometrically, the objects of our sample show a diverse late-time behaviour, some of them fading quite rapidly after a light-curve break at ∼ 150-200 d. The latter is likely the result of flux redistribution into the infrared, possibly caused by dust formation, rather than a true bolometric effect. Nebular spectra of superluminous SNe Ia are characterized by weak or absent [Fe III] emission, pointing at a low ejecta ionization state as a result of high densities. To constrain the ejecta and Ni masses of superluminous SNe Ia, we compare the observed bolometric light curve of SN 2009dc with synthetic model light curves, focusing on the radioactive tail after ∼60 d. Models with enough Ni to explain the light-curve peak by radioactive decay, and at the same time sufficient mass to keep the ejecta velocities low, fail to reproduce the observed light-curve tail of SN 2009dc because of too much γ -ray trapping.We instead propose a model with ∼1M of Ni and ∼2 M of ejecta, which may be interpreted as the explosion of a Chandrasekhar-mass white dwarf (WD) enshrouded by 0.6-0.7 M of C/O-rich material, as it could result from a merger of two massive C/O WDs. This model reproduces the late light curve of SN 2009dc well. A flux deficit at peak may be compensated by light from the interaction of the ejecta with the surrounding material.

Interacting supernovae from photoionization-confined shells around red supergiant stars

Betelgeuse, a nearby red supergiant, is a runaway star with a powerful stellar wind that drives a bow shock into its surroundings. This picture has been challenged by the discovery of a dense and almost static shell that is three times closer to the star than the bow shock and has been decelerated by some external force. The two physically distinct structures cannot both be formed by the hydrodynamic interaction of the wind with the interstellar medium. Here we report that a model in which Betelgeuse's wind is photoionized by radiation from external sources can explain the static shell without requiring a new understanding of the bow shock. Pressure from the photoionized wind generates a standing shock in the neutral part of the wind and forms an almost static, photoionization-confined shell. Other red supergiants should have significantly more massive shells than Betelgeuse, because the photoionization-confined shell traps up to 35 per cent of all mass lost during the red supergiant phase, confining this gas close to the star until it explodes. After the supernova explosion, massive shells dramatically affect the supernova lightcurve, providing a natural explanation for the many supernovae that have signatures of circumstellar interaction.

Superluminous supernovae from PESSTO

We present optical spectra and light curves for three hydrogen-poor superluminous supernovae followed by the Public ESO Spectroscopic Survey of Transient Objects (PESSTO). Time series spectroscopy from a fewdays aftermaximum light to 100 d later shows them to be fairly typical of this class, with spectra dominated by Ca II, MgII, FeII, and Si II, which evolve slowly over most of the post-peak photospheric phase. We determine bolometric light curves and apply simple fitting tools, based on the diffusion of energy input by magnetar spin-down, Ni-56 decay, and collision of the ejecta with an opaque circumstellar shell. We investigate how the heterogeneous light curves of our sample (combined with others from the literature) can help to constrain the possible mechanisms behind these events. We have followed these events to beyond 100-200 d after peak, to disentangle host galaxy light from fading supernova flux and to differentiate between the models, which predict diverse behaviour at this phase. Models powered by radioactivity require unrealistic parameters to reproduce the observed light curves, as found by previous studies. Both magnetar heating and circumstellar interaction still appear to be viable candidates. A large diversity is emerging in observed tail-phase luminosities, with magnetar models failing in some cases to predict the rapid drop in flux. This would suggest either that magnetars are not responsible, or that the X-ray flux from the magnetar wind is not fully trapped. The light curve of one object shows a distinct rebrightening at around 100 d after maximum light. We argue that this could result either from multiple shells of circumstellar material, or from a magnetar ionization front breaking out of the ejecta.

Interacting supernovae and supernova impostors. SN 2007sv: the major eruption of a massive star in UGC 5979

We report the results of the photometric and spectroscopic monitoring campaign of the transient SN 2007sv. The observables are similar to those of Type IIn supernovae, a well-known class of objects whose ejecta interact with pre-existing circumstellar material (CSM). The spectra show a blue continuum at early phases and prominent Balmer lines in emission; however, the absolute magnitude at the discovery of SN 2007sv (M-R=-14.25 +/- 0.38) indicate it to be most likely a supernova impostor. This classification is also supported by the lack of evidence in the spectra of very high velocity material as expected in supernova ejecta. In addition, we find no unequivocal evidence of broad lines of alpha- and/or Fe-peak elements. The comparison with the absolute light curves of other interacting objects (including Type IIn supernovae) highlights the overall similarity with the prototypical impostor SN 1997bs. This supports our claim that SN 2007sv was not a genuine supernova, and was instead a supernova impostor, most likely similar to the major eruption of a luminous blue variable.

Detecting the progenitors of core collapse supernovae

The masses and the evolutionary states of the progenitors of core-collapse supernovae are not well constrained by direct observations. Stellar evolution theory generally predicts that massive stars with initial masses less than about 30M_sol should undergo core-collapse when they are cool M-type supergiants. However the only two detections of a SN progenitor before explosion are SN1987A and SN1993J, and neither of these was an M-type supergiant. Attempting to identify the progenitors of supernovae is a difficult task, as precisely predicting the time of explosion of a massive star is impossible for obvious reasons. There are several different types of supernovae which have different spectral and photometric evolution, and how exactly these are related to the evolutionary states of the progenitor stars is not currently known. I will describe a novel project which may allow the direct identification of core-collapse supernovae progenitors on pre-explosion images of resolved, nearby galaxies. This project is now possible with the excellent image archives maintained by several facilities and will be enhanced by the new initiatives to create Virtual Observatories, the earliest of which ASTROVIRTEL is already producing results.