Mass limits for the progenitor star of supernova 2001du and other Type II-P supernovae

The supernova SN 2001du was discovered in the galaxy NGC 1365 at a distance of 19 +/- 2 Mpc, and is a core-collapse event of Type II-P. Images of this galaxy, of moderate depth, have been taken with the Hubble Space Telescope approximately 6.6 yr before discovery and include the supernova position on the WFPC2 field of view. We have observed the supernova with the WFPC2 to allow accurate differential astrometry of SN 2001du on the pre-explosion frames. As a core-collapse event it is expected that the progenitor was a massive, luminous star. There is a marginal detection (3sigma) of a source close to the supernova position on the pre-discovery V -band frame, but it is not precisely coincident and we do not believe it to be a robust detection of a point source. We conclude that there is no stellar progenitor at the supernova position and derive sensitivity limits of the pre-discovery images that provide an upper mass limit for the progenitor star. We estimate that the progenitor had a mass of less than 15 M-circle dot . We revisit two other nearby Type II-P supernovae that have high-quality pre-explosion images, and refine the upper mass limits for the progenitor stars. Using a new distance determination for SN 1999gi from the expanding photosphere method, we revise the upper mass limit to 12 M-circle dot . We present new HST images of the site of SN 1999em, which validate the use of lower spatial resolution ground-based images in the progenitor studies and use a new Cepheid distance to the galaxy to measure an upper mass limit of 15 M-circle dot for that progenitor. Finally we compile all the direct information available for the progenitors of eight nearby core-collapse supernovae and compare their mass estimates. These are compared with the latest stellar evolutionary models of pre-supernova evolution which have attempted to relate metallicity and mass to the supernovae type. Although this is statistically limited at present, reasonable agreement is already found for the lower-mass events (generally the II-P), but some discrepancies appear at higher masses.

A comparison between star formation rate diagnostics and rate of core collapse supernovae within 11 Mpc

Aims. The core collapse supernova rate provides a strong lower limit for the star formation rate (SFR). Progress in using it as a cosmic SFR tracer requires some confidence that it is consistent with more conventional SFR diagnostics in the nearby Universe. This paper compares standard SFR measurements based on H alpha, far ultraviolet (FUV) and total infrared (TIR) galaxy luminosities with the observed core collapse supernova rate in the same galaxy sample. The comparison can be viewed from two perspectives. Firstly, by adopting an estimate of the minimum stellar mass to produce a core collapse supernova one can determine a SFR from supernova numbers. Secondly, the radiative SFR can be assumed to be robust and then the supernova statistics provide a constrain on the minimum stellar mass for core collapse supernova progenitors.

The Unification of Asymmetry Signatures of Type Ia Supernovae

We present a compilation of the geometry measures acquired using optical and IR spectroscopy and optical spectropolarimetry to probe the explosion geometry of Type Ia supernovae (SNe Ia). Polarization measurements are sensitive to asymmetries in the plane of the sky, whereas line profiles in nebular phase spectra are expected to trace asymmetries perpendicular to the plane of the sky. The combination of these two measures can overcome their respective projection effects, completely probing the structures of these events. For nine normal SNe Ia, we find that the polarization of Si II ?6355 at 5 days before maximum (p Si II ) is well correlated with its velocity evolution (\dot{v}_Si II), implying that \dot{v}_Si II is predominantly due to the asymmetry of the SNe. We find only a weak correlation between the polarization of Si II and the reported velocities (v neb) for peak emission of optical Fe II and Ni II lines in nebular spectra. Our sample is biased, with polarization measurements being only available for normal SNe that subsequently exhibited positive (i.e., redshifted) v neb. In unison these indicators are consistent with an explosion in which the outer layers are dominated by a spherical oxygen layer, mixed with an asymmetric distribution of intermediate-mass elements. The combination of spectroscopic and spectropolarimetric indicators suggests a single geometric configuration for normal SNe Ia, with some of the diversity of observed properties arising from orientation effects.

Double-detonation sub-Chandrasekhar supernovae: Synthetic observables for minimum helium shell mass models

In the double-detonation scenario for Type Ia supernovae, it is suggested that a detonation initiates in a shell of helium-rich material accreted from a companion star by a sub-Chandrasekhar-mass white dwarf. This shell detonation drives a shock front into the carbon-oxygen white dwarf that triggers a secondary detonation in the core. The core detonation results in a complete disruption of the white dwarf. Earlier studies concluded that this scenario has difficulties in accounting for the observed properties of Type Ia supernovae since the explosion ejecta are surrounded by the products of explosive helium burning in the shell. Recently, however, it was proposed that detonations might be possible for much less massive helium shells than previously assumed (Bildsten et al.). Moreover, it was shown that even detonations of these minimum helium shell masses robustly trigger detonations of the carbon-oxygen core (Fink et al.). Therefore, it is possible that the impact of the helium layer on observables is less than previously thought. Here, we present time-dependent multi-wavelength radiative transfer calculations for models with minimum helium shell mass and derive synthetic observables for both the optical and ? -ray spectral regions. These differ strongly from those found in earlier simulations of sub-Chandrasekhar-mass explosions in which more massive helium shells were considered. Our models predict light curves that cover both the range of brightnesses and the rise and decline times of observed Type Ia supernovae. However, their colors and spectra do not match the observations. In particular, their B - V colors are generally too red. We show that this discrepancy is mainly due to the composition of the burning products of the helium shell of the Fink et al. models which contain significant amounts of titanium and chromium. Using a toy model, we also show that the burning products of the helium shell depend crucially on its initial composition. This leads us to conclude that good agreement between sub-Chandrasekhar-mass explosions and observed Type Ia supernovae may still be feasible but further study of the shell properties is required.

Delay times and rates for Type Ia supernovae and thermonuclear explosions from double-detonation sub-Chandrasekhar mass models

We present theoretical delay times and rates of thermonuclear explosions that are thought to produce Type Ia supernovae (SNe Ia), including the double-detonation sub-Chandrasekhar mass model, using the population synthesis binary evolution code startrack. If detonations of sub-Chandrasekhar mass carbon-oxygen white dwarfs following a detonation in an accumulated layer of helium on the white dwarf's surface ('double-detonation' models) are able to produce thermonuclear explosions which are characteristically similar to those of SNe Ia, then these sub-Chandrasekhar mass explosions may account for at least some substantial fraction of the observed SN Ia rate. Regardless of whether all double-detonations look like 'normal' SNe Ia, in any case the explosions are expected to be bright and thus potentially detectable. Additionally, we find that the delay time distribution of double-detonation sub-Chandrasekhar mass SNe Ia can be divided into two distinct formation channels: the 'prompt' helium-star channel with delay times

Double-detonation sub-Chandrasekhar supernovae: Can minimum helium shell masses detonate the core?

The explosion of sub-Chandrasekhar mass white dwarfs via the double detonation scenario is a potential explanation for type Ia supernovae. In this scenario, a surface detonation in a helium layer initiates a detonation in the underlying carbon/oxygen core leading to an explosion. For a given core mass, a lower bound has been determined on the mass of the helium shell required for dynamical burning during a helium flash, which is a necessary prerequisite for detonation. For a range of core and corresponding minimum helium shell masses, we investigate whether an assumed surface helium detonation is capable of triggering a subsequent detonation in the core even for this limiting case. We carried out hydrodynamic simulations on a co-expanding Eulerian grid in two dimensions assuming rotational symmetry. The detonations are propagated using the level-set approach and a simplified scheme for nuclear reactions that has been calibrated with a large nuclear network. The same network is used to determine detailed nucleosynthetic abundances in a post-processing step. Based on approximate detonation initiation criteria in the literature, we find that secondary core detonations are triggered for all of the simulated models, ranging in core mass from 0.810 up to 1.385 M? with corresponding shell masses from 0.126 down to 0.0035 M?. This implies that, as soon as a detonation triggers in a helium shell covering a carbon/oxygen white dwarf, a subsequent core detonation is virtually inevitable.

Light curves for off-centre ignition models of Type Ia supernovae

Motivated by recent models involving off-centre ignition of Type Ia supernova explosions, we undertake three-dimensional time-dependent radiation transport simulations to investigate the range of bolometric light-curve properties that could be observed from supernovae in which there is a lop-sided distribution of the products from nuclear burning. We consider both a grid of artificial toy models which illustrate the conceivable range of effects and a recent three-dimensional hydrodynamical explosion model. We find that observationally significant viewing angle effects are likely to arise in such supernovae and that these may have important ramifications for the interpretation of the observed diversity of Type Ia supernova and the systematic uncertainties which relate to their use as standard candles in contemporary cosmology. © 2007 RAS.

A theoretical color-velocity correlation for supernovae associated with gamma-ray bursts

We carry out the first multi-dimensional radiative transfer calculations to simultaneously compute synthetic spectra and light curves for models of supernovae driven by fast bipolar outflows. These allow us to make self-consistent predictions for the orientation dependence of both color evolution and spectral features. We compare models with different degrees of asphericity and metallicity and find significant observable consequences of both. In aspherical models, we find spectral and light curve features that vary systematically with observer orientation. In particular, we find that the early-phase light curves are brighter and bluer when viewed close to the polar axis but that the peak flux is highest for equatorial (off-axis) inclinations. Spectral line features also depend systematically on observer orientation, including the velocity of the Si II 6355 Å line. Consequently, our models predict a correlation between line velocity and color that could assist the identification of supernovae associated with off-axis jet-driven explosions. The amplitude and range of this correlation depends on the degree of asphericity, the metallicity, and the epoch of observation but we find that it is always present and acts in the same direction. © 2012. The American Astronomical Society. All rights reserved..

Normal type Ia supernovae from violent mergers of white dwarf binaries

One of the most important questions regarding the progenitor systems of Type Ia supernovae (SNe Ia) is whether mergers of two white dwarfs can lead to explosions that reproduce observations of normal events. Here we present a fully three-dimensional simulation of a violent merger of two carbon-oxygen white dwarfs with masses of 0.9 M and 1.1 M combining very high resolution and exact initial conditions. A well-tested combination of codes is used to study the system. We start with the dynamical inspiral phase and follow the subsequent thermonuclear explosion under the plausible assumption that a detonation forms in the process of merging. We then perform detailed nucleosynthesis calculations and radiative transfer simulations to predict synthetic observables from the homologously expanding supernova ejecta. We find that synthetic color light curves of our merger, which produces about 0.62 M of Ni, show good agreement with those observed for normal SNe Ia in all wave bands from U to K. Line velocities in synthetic spectra around maximum light also agree well with observations. We conclude that violent mergers of massive white dwarfs can closely resemble normal SNe Ia. Therefore, depending on the number of such massive systems available these mergers may contribute at least a small fraction to the observed population of normal SNe Ia. © 2012 The American Astronomical Society. All rights reserved.

On the γ-ray emission of Type Ia supernovae

A multidimension, time-dependent Monte Carlo code is used to compute sample ?-ray spectra to explore whether unambiguous constraints could be obtained from ?-ray observations of Type Ia supernovae. Both spherical and aspherical geometries are considered and it is shown that moderate departures from sphericity can produce viewing-angle effects that are at least as significant as those caused by the variation of key parameters in 1D models. Thus, ?-ray data could, in principle, carry some geometrical information, and caution should be applied when discussing the value of ?-ray data based only on 1D explosion models. In light of the limited sensitivity of current ?-ray observatories, the computed theoretical spectra are studied to revisit the issue of whether useful constraints could be obtained for moderately nearby objects. The most useful ?-ray measurements are likely to be of the light curve and time-dependent hardness ratios, but sensitivity higher than currently available, particularly at relatively hard energies (~2-3 MeV), is desirable. © 2008 The Authors. Journal compilation © 2008 RAS.

Time-dependent three-dimensional spectrum synthesis for Type Ia supernovae

A Monte Carlo code (artis) for modelling time-dependent three-dimensional spectral synthesis in chemically inhomogeneous models of Type Ia supernova ejecta is presented. Following the propagation of ?-ray photons, emitted by the radioactive decay of the nucleosynthesis products, energy is deposited in the supernova ejecta and the radiative transfer problem is solved self-consistently, enforcing the constraint of energy conservation in the comoving frame. Assuming a photoionization-dominated plasma, the equations of ionization equilibrium are solved together with the thermal balance equation adopting an approximate treatment of excitation. Since we implement a fully general treatment of line formation, there are no free parameters to adjust. Thus, a direct comparison between synthetic spectra and light curves, calculated from hydrodynamic explosion models, and observations is feasible. The code is applied to the well-known W7 explosion model and the results tested against other studies. Finally, the effect of asymmetric ejecta on broad-band light curves and spectra is illustrated using an elliptical toy model. © 2009 RAS.

Off-center explosions of Chandrasekhar-mass white dwarfs:An explanation of super-bright type la supernovae?

Context. The recent discovery of a very bright type la supernova, SNLS-03D3bb (=SN 2003fg), in the Supernova Legacy Survey (SNLS) has raised the question of whether super-Chandrasekhar-mass white-dwarf stars are needed to explain such bright explosions. Progenitors of this sort could form by mergers of pairs of rather massive white dwarfs. Binary systems of two white dwarfs in close orbit, where their total mass significantly exceeds the Chandrasekhar mass, have not yet been found. Therefore SNLS-03D3bb could establish the first clear case of a double-degenerate progenitor of a (peculiar) type la supernovae. Moreover, if this interpretation is correct, it casts some doubt on the universality of the calibration relations used to make SNe la distance indicators for cosmology. Aims. We aim to evaluate the case for a super-Chandrasekhar-mass progenitor for SNLS-03D3bb in light of previous theoretical work on super-Chandrasekhar-mass explosions. Furthermore, we propose an alternative scenario involving only a Chandrasekhar-mass progenitor. Methods. We present a theoretically motivated critical discussion of the expected observational fingerprints of super-Chandrasekharmass explosions. As an alternative, we describe a simple class of aspherical Chandrasekhar-mass models in which the products of nuclear burning are displaced from the center. We then perform simple radiative transfer calculations to predict synthetic lightcurves for one such off-center explosion model. Results. In important respects, the expected observational consequences of super-Chandrasekhar-mass explosions are not consistent with the observations of SNLS-03D3bb. We demonstrate that the lopsided explosion of a Chandrasekhar-mass white dwarf could provide a better explanation. © ESO 2007.

Multidimensional simulations of radiative transfer in Type Ia supernovae

A three-dimensional Monte Carlo code for modelling radiation transport in Type Ia supernovae is described. In addition to tracking Monte Carlo quanta to follow the emission, scattering and deposition of radiative energy, a scheme involving volume-based Monte Carlo estimators is used to allow properties of the emergent radiation field to be extracted for specific viewing angles in a multidimensional structure. This eliminates the need to compute spectra or light curves by angular binning of emergent quanta. The code is applied to two test problems to illustrate consequences of multidimensional structure on the modelling of light curves. First, elliptical models are used to quantify how large-scale asphericity can introduce angular dependence to light curves. Secondly, a model which incorporates complex structural inhomogeneity, as predicted by modern explosion models, is used to investigate how such structure may affect light-curve properties. © 2006 RAS.

On the brightness distribution of type ia supernovae from violent white dwarf mergers

We investigate the brightness distribution expected for thermonuclear explosions that might result from the ignition of a detonation during the violent merger of white dwarf (WD) binaries. Violent WD mergers are a subclass of the canonical double degenerate scenario where two carbon-oxygen (CO) WDs merge when the larger WD fills its Roche lobe. Determining their brightness distribution is critical for evaluating whether such an explosion model could be responsible for a significant fraction of the observed population of Type Ia supernovae (SNe Ia). We argue that the brightness of an explosion realized via the violent merger model is mainly determined by the mass of Ni produced in the detonation of the primary COWD. To quantify this link, we use a set of sub-Chandrasekhar mass WD detonation models to derive a relationship between primary WD mass (m) and expected peak bolometric brightness (M). We use this m-M relationship to convert the masses of merging primary WDs from binary population models to a predicted distribution of explosion brightness. We also investigate the sensitivity of our results to assumptions about the conditions required to realize a detonation during violent mergers ofWDs. We find a striking similarity between the shape of our theoretical peak-magnitude distribution and that observed for SNe Ia: our model produces a M distribution that roughly covers the range and matches the shape of the one observed for SNe Ia. However, this agreement hinges on a particular phase of mass accretion during binary evolution: the primary WD gains ~0.15-0.35M? from a slightly evolved helium star companion. In our standard binary evolution model, such an accretion phase is predicted to occur for about 43 per cent of all binary systems that ultimately give rise to binary CO WD mergers. We also find that with high probability, violent WD mergers involving the most massive primaries (?1.3M?, which should produce bright SNe) have delay times ?500 Myr. © 2012 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.

3D deflagration simulations leaving bound remnants:A model for 2002cx-like Type Ia supernovae

2002cx-like supernovae are a sub-class of sub-luminous Type Ia supernovae (SNe). Their light curves and spectra are characterized by distinct features that indicate strong mixing of the explosion ejecta. Pure turbulent deflagrations have been shown to produce such mixed ejecta. Here, we present hydrodynamics, nucleosynthesis and radiative-transfer calculations for a 3D full-star deflagration of a Chandrasekhar-mass white dwarf. Our model is able to reproduce the characteristic observational features of SN 2005hk (a prototypical 2002cx-like supernova), not only in the optical, but also in the near-infrared. For that purpose we present, for the first time, five near-infrared spectra of SN 2005hk from -0.2 to 26.6 d with respect to B-band maximum. Since our model burns only small parts of the initial white dwarf, it fails to completely unbind the white dwarf and leaves behind a bound remnant of ~1.03Mconsisting mainly of unburned carbon and oxygen, but also enriched by some amount of intermediate-mass and iron-group elements from the explosion products that fall back on the remnant.We discuss possibilities for detecting this bound remnant and how it might influence the late-time observables of 2002cx-like SNe. © 2013 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society.