A deeper search for the progenitor of the Type Ic supernova 2002ap

Images of the site of the Type Ic supernova (SN) 2002ap taken before explosion were analysed previously by Smartt et al. We have uncovered new unpublished, archival pre-explosion images from the Canada-France-Hawaii Telescope (CFHT) that are vastly superior in depth and image quality. In this paper we present a further search for the progenitor star of this unusual Type Ic SN. Aligning high-resolution Hubble Space Telescope observations of the SN itself with the archival CFHT images allowed us to pinpoint the location of the progenitor site on the groundbased observations. We find that a source visible in the B- and R-band pre-explosion images close to the position of the SN is (1) not coincident with the SN position within the uncertainties of our relative astrometry and (2) is still visible similar to 4.7-yr post-explosion in late-time observations taken with the William Herschel Telescope. We therefore conclude that it is not the progenitor of SN 2002ap. We derived absolute limiting magnitudes for the progenitor of M-B >= -4.2 +/- 0.5 and M-R >= -5.1 +/- 0.5. These are the deepest limits yet placed on a Type Ic SN progenitor. We rule out all massive stars with initial masses greater than 7-8 M-circle dot (the lower mass limit for stars to undergo core collapse) that have not evolved to become Wolf-Rayet stars. This is consistent with the prediction that Type Ic SNe should result from the explosions of Wolf-Rayet stars. Comparing our luminosity limits with stellar models of single stars at appropriate metallicity (Z = 0.008) and with standard mass-loss rates, we find no model that produces a Wolf-Rayet star of low enough mass and luminosity to be classed as a viable progenitor. Models with twice the standard mass-loss rates provide possible single star progenitors but all are initially more massive than 30-40 M-circle dot. We conclude that any single star progenitor must have experienced at least twice the standard mass-loss rates, been initially more massive than 30-40 M-circle dot and exploded as a Wolf-Rayet star of final mass 10-12 M-circle dot. Alternatively a progenitor star of lower initial mass may have evolved in an interacting binary system. Mazzali et al. propose such a binary scenario for the progenitor of SN 2002ap in which a star of initial mass 15-20 M-circle dot is stripped by its binary companion, becoming a 5 M-circle dot Wolf-Rayet star prior to explosion. We constrain any possible binary companion to a main-sequence star of

The binary progenitor of Tycho Brahe's 1572 supernova

The brightness of type Ia supernovae, and their homogeneity as a class, makes them powerful tools in cosmology, yet little is known about the progenitor systems of these explosions. They are thought to arise when a white dwarf accretes matter from a companion star, is compressed and undergoes a thermonuclear explosion(1-3). Unless the companion star is another white dwarf ( in which case it should be destroyed by the mass-transfer process itself), it should survive and show distinguishing properties. Tycho's supernova(4,5) is one of only two type Ia supernovae observed in our Galaxy, and so provides an opportunity to address observationally the identification of the surviving companion. Here we report a survey of the central region of its remnant, around the position of the explosion, which excludes red giants as the mass donor of the exploding white dwarf. We found a type G0 - G2 star, similar to our Sun in surface temperature and luminosity ( but lower surface gravity), moving at more than three times the mean velocity of the stars at that distance, which appears to be the surviving companion of the supernova.

Detonations in sub-chandrasekhar-mass C+O white dwarfs

Explosions of sub-Chandrasekhar-mass white dwarfs (WDs) are one alternative to the standard Chandrasekhar-mass model of Type Ia supernovae (SNe Ia). They are interesting since binary systems with sub-Chandrasekhar-mass primary WDs should be common and this scenario would suggest a simple physical parameter which determines the explosion brightness, namely the mass of the exploding WD. Here we perform one-dimensional hydrodynamical simulations, associated post-processing nucleosynthesis, and multi-wavelength radiation transport calculations for pure detonations of carbon-oxygen WDs. The light curves and spectra we obtain from these simulations are in good agreement with observed properties of SNe Ia. In particular, for WD masses from 0.97 to 1.15 Msun we obtain 56Ni masses between 0.3 and 0.8 Msun, sufficient to capture almost the complete range of SN Ia brightnesses. Our optical light curve rise times, peak colors, and decline timescales display trends which are generally consistent with observed characteristics although the range of B-band decline timescales displayed by our current set of models is somewhat too narrow. In agreement with observations, the maximum light spectra of the models show clear features associated with intermediate-mass elements and reproduce the sense of the observed correlation between explosion luminosity and the ratio of the Si II lines at ?6355 and ?5972. We therefore suggest that sub-Chandrasekhar-mass explosions are a viable model for SNe Ia for any binary evolution scenario leading to explosions in which the optical display is dominated by the material produced in a detonation of the primary WD. © 2010. The American Astronomical Society.

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.

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.

2D simulations of the double-detonation model for thermonuclear transients from low-mass carbon-oxygen white dwarfs

Thermonuclear explosions may arise in binary star systems in which a carbon-oxygen (CO) white dwarf (WD) accretes helium-rich material from a companion star. If the accretion rate allows a sufficiently large mass of helium to accumulate prior to ignition of nuclear burning, the helium surface layer may detonate, giving rise to an astrophysical transient. Detonation of the accreted helium layer generates shock waves that propagate into the underlying CO WD. This might directly ignite a detonation of the CO WD at its surface (an edge-lit secondary detonation) or compress the core of the WD sufficiently to trigger a CO detonation near the centre. If either of these ignition mechanisms works, the two detonations (helium and CO) can then release sufficient energy to completely unbind the WD. These 'double-detonation' scenarios for thermonuclear explosion of WDs have previously been investigated as a potential channel for the production of Type Ia supernovae from WDs of ~ 1 M . Here we extend our 2D studies of the double-detonation model to significantly less massive CO WDs, the explosion of which could produce fainter, more rapidly evolving transients. We investigate the feasibility of triggering a secondary core detonation by shock convergence in low-mass CO WDs and the observable consequences of such a detonation. Our results suggest that core detonation is probable, even for the lowest CO core masses that are likely to be realized in nature. To quantify the observable signatures of core detonation, we compute spectra and light curves for models in which either an edge-lit or compression-triggered CO detonation is assumed to occur. We compare these to synthetic observables for models in which no CO detonation was allowed to occur. If significant shock compression of the CO WD occurs prior to detonation, explosion of the CO WD can produce a sufficiently large mass of radioactive iron-group nuclei to significantly affect the light curves. In particular, this can lead to relatively slow post-maximum decline. If the secondary detonation is edge-lit, however, the CO WD explosion primarily yields intermediate-mass elements that affect the observables more subtly. In this case, near-infrared observations and detailed spectroscopic analysis would be needed to determine whether a core detonation occurred. We comment on the implications of our results for understanding peculiar astrophysical transients including SN 2002bj, SN 2010X and SN 2005E. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 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.

Monte Carlo radiation hydrodynamics:Methods, tests and application to Type Ia supernova ejecta

In astrophysical systems, radiation-matter interactions are important in transferring energy and momentum between the radiation field and the surrounding material. This coupling often makes it necessary to consider the role of radiation when modelling the dynamics of astrophysical fluids. During the last few years, there have been rapid developments in the use of Monte Carlo methods for numerical radiative transfer simulations. Here, we present an approach to radiation hydrodynamics that is based on coupling Monte Carlo radiative transfer techniques with finite-volume hydrodynamical methods in an operator-split manner. In particular, we adopt an indivisible packet formalism to discretize the radiation field into an ensemble of Monte Carlo packets and employ volume-based estimators to reconstruct the radiation field characteristics. In this paper the numerical tools of this method are presented and their accuracy is verified in a series of test calculations. Finally, as a practical example, we use our approach to study the influence of the radiation-matter coupling on the homologous expansion phase and the bolometric light curve of Type Ia supernova explosions. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 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.

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.

The host galaxy of the super-luminous SN 2010gx and limits on explosive nickel-56 production

Super-luminous supernovae have a tendency to occur in faint host galaxies which are likely to have low mass and low metallicity. While these extremely luminous explosions have been observed from z=0.1 to 1.55, the closest explosions allow more detailed investigations of their host galaxies. We present a detailed analysis of the host galaxy of SN 2010gx (z=0.23), one of the best studied super-luminous type Ic supernovae. The host is a dwarf galaxy (M_g=-17.42+/-0.17) with a high specific star formation rate. It has a remarkably low metallicity of 12+log(O/H)=7.5+/-0.1 dex as determined from the detection of the [OIII] 4363 Angs line. This is the first reliable metallicity determination of a super-luminous stripped-envelope supernova host. We collected deep multi-epoch imaging with Gemini + GMOS between 240-560 days after explosion to search for any sign of radioactive nickel-56, which might provide further insights on the explosion mechanism and the progenitor's nature. We reach griz magnitudes of m_AB~26, but do not detect SN 2010gx at these epochs. The limit implies that any nickel-56 production was similar to or below that of SN 1998bw (a luminous type Ic SN that produced around 0.4 M_sun of nickel-56). The low volumetric rates of these supernovae (~10^-4 of the core-collapse population) could be qualitatively matched if the explosion mechanism requires a combination of low-metallicity (below 0.2 Z_sun), high progenitor mass (>60 M_sun) and high rotation rate (fastest 10% of rotators).

The superluminous supernova PS1-11ap: bridging the gap between low and high redshift

We present optical photometric and spectroscopic coverage of the superluminous supernova (SLSN) PS1-11ap, discovered with the Pan-STARRS1 Medium Deep Survey at z = 0.524. This intrinsically blue transient rose slowly to reach a peak magnitude of Mu = −21.4 mag and bolometric luminosity of 8 × 1043 erg s−1 before settling on to a relatively shallow gradient of decline. The observed decline is significantly slower than those of the SLSNe-Ic which have been the focus of much recent attention. Spectroscopic similarities with the lower redshift SN2007bi and a decline rate similar to 56Co decay time-scale initially indicated that this transient could be a candidate for a pair instability supernova (PISN) explosion. Overall the transient appears quite similar to SN2007bi and the lower redshift object PTF12dam. The extensive data set, from 30 d before peak to 230 d after, allows a detailed and quantitative comparison with published models of PISN explosions. We find that the PS1-11ap data do not match these model explosion parameters well, supporting the recent claim that these SNe are not pair instability explosions. We show that PS1-11ap has many features in common with the faster declining SLSNe-Ic, and the light-curve evolution can also be quantitatively explained by the magnetar spin-down model. At a redshift of z = 0.524, the observer-frame optical coverage provides comprehensive rest-frame UV data and allows us to compare it with the SLSNe recently found at high redshifts between z = 2 and 4. While these high-z explosions are still plausible PISN candidates, they match the photometric evolution of PS1-11ap and hence could be counterparts to this lower redshift transient.

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.