• Double detonations in double white dwarf (WD) binaries undergoing unstable mass transfer have emerged in recent years as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. One potential outcome of this "dynamically driven double-degenerate double-detonation" (D^6) scenario is that the companion WD survives the explosion and is flung away with a velocity equal to its > 1000 km/s pre-SN orbital velocity. We perform a search for these hypervelocity runaway WDs using Gaia's second data release. In this paper, we discuss seven candidates followed up with ground-based instruments. Three sources are likely to be some of the fastest known stars in the Milky Way and are consistent with having previously been companion WDs in pre-SN Ia systems. The combined five-parameter astrometric solutions from Gaia and radial velocities from follow-up spectra yield tentative 6D confirmation of the D^6 scenario. The past position of one of these stars is also coincident with the center of a faint, old SN remnant, further strengthening the interpretation of these candidates as hypervelocity runaways from binary systems that underwent SNe Ia.
  • The recently discovered bright type II supernova OGLE14-073 evolved very slowly. The light curve rose to maximum for 90 days from discovery and then declined at a rate compatible with the radioactive decay of 56Co. In this study, we show that a pair-instability supernova is a plausible mechanism for this event. We calculate explosion models and light curves with the radiation hydrodynamics code STELLA starting from two M(ZAMS) = 150 solar masses, Z=0.001 progenitors. We obtain satisfactory fits to OGLE14-073 broadband light curves by including additional 56Ni in the centre of the models and mixing hydrogen down into the inner layers of the ejecta to a radial mass coordinate of 10 solar masses. The extra 56Ni required points to a slightly more massive progenitor star. The mixing of hydrogen could be due to large scale mixing during the explosion. We also present synthetic spectra for our models simulated with the Monte Carlo radiative transfer code ARTIS. The synthetic spectra reproduce the main features of the observed spectra of OGLE14-073. We conclude that OGLE14-073 is one of the most promising candidates for a pair-instability explosion.
  • We present a detailed spectral analysis of Type Iax SN 2011ay. Our spectra cover epochs between -3 and +19 days with respect to the maximum light in B-band. This time range allows us to employ a so-called abundance tomography technique. The synthetic spectral fitting was made with the 1D Monte Carlo radiative transfer code TARDIS. In this paper, we describe our method to fit multiple epochs with a self-consistent, stratified atmospheric model. We compare our results to previously published SYN++ models and the predictions of different explosion scenarios. Using a fixed density profile (exponential fit of W7), we find that a uniform abundance model cannot reproduce the spectral features before maximum light because of the emergence of excessively strong Fe lines. In our best-fit TARDIS model, we find an abundance profile that separated into two different regimes: a well-mixed region under 10,000 km s$^{-1}$ and a stratified region with decreasing IGE abundances above 10,000 km s$^{-1}$. Based on a detailed comparative analysis, our conclusion is that the available pure deflagration models cannot fully explain either the observed properties of SN 2011ay or the results of our TARDIS modeling. Further examinations are necessary to find an adequate explanation for the origin of this object.
  • Despite the importance of Type Ia supernovae (SNe Ia) for modern astrophysics, their detailed mechanism is still not fully understood. In this contribution, we present recent findings from numerical explosion models in the context of the observed diversity of SNe Ia and we discuss how these models can help to shed light on the explosion mechanism and the progenitor stars of SNe Ia. In addition, we introduce the Heidelberg Supernova Model Archive (HESMA), a new online data base where we provide integrated isotopic abundances and radially averaged ejecta profiles and synthetic observables for a wide range of state-of-the-art explosion models.
  • The gravitationally confined detonation (GCD) model has been proposed as a possible explosion mechanism for Type Ia supernovae in the single-degenerate evolution channel. Driven by buoyancy, a deflagration flame rises in a narrow cone towards the surface. For the most part, the flow of the expanding ashes remains radial, but upon reaching the outer, low-pressure layers of the white dwarf, an additional lateral component develops. This makes the deflagration ashes converge again at the opposite side, where the compression heats fuel and a detonation may be launched. To test the GCD explosion model, we perform a 3D simulation for a model with an ignition spot offset near the upper limit of what is still justifiable, 200 km. This simulation meets our deliberately optimistic detonation criteria and we initiate a detonation. The detonation burns through the white dwarf and leads to its complete disruption. We determine nucleosynthetic yields by post-processing 10^6 tracer particles with a 384 nuclide reaction network and we present multi-band light curves and time-dependent optical spectra. We find that our synthetic observables show a prominent viewing-angle sensitivity in UV and blue bands, which is in tension with observed SNe Ia. The strong dependence on viewing-angle is caused by the asymmetric distribution of the deflagration ashes in the outer ejecta layers. Finally, we perform a comparison of our model to SN 1991T. The overall flux-level of the model is slightly too low and the model predicts pre-maximum light spectral features due to Ca, S, and Si that are too strong. Furthermore, the model chemical abundance stratification qualitatively disagrees with recent abundance tomography results in two key areas: our model lacks low velocity stable Fe and instead has copious amounts of high-velocity 56Ni and stable Fe. We therefore do not find good agreement of the model with SN 1991T.
  • 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. 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. 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. 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. 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 (1.2 Msun) shows that the less blueshifted spectral features fit the observations better, although they are too bright for normal SNe Ia.
  • Binary stellar evolution calculations predict that Chandrasekhar-mass carbon/oxygen white dwarfs (WDs) show a radially varying profile for the composition with a carbon depleted core. Many recent multi-dimensional simulations of Type Ia supernovae (SNe Ia), however, assume the progenitor WD has a homogeneous chemical composition. In this work, we explore the impact of different initial carbon profiles of the progenitor WD on the explosion phase and on synthetic observables in the Chandrasekhar-mass delayed detonation model. Spectra and light curves are compared to observations to judge the validity of the model. The explosion phase is simulated using the finite volume supernova code LEAFS, which is extended to treat different compositions of the progenitor WD. The synthetic observables are computed with the Monte Carlo radiative transfer code ARTIS. Differences in binding energies of carbon and oxygen lead to a lower nuclear energy release for carbon depleted material; thus, the burning fronts that develop are weaker and the total nuclear energy release is smaller. For otherwise identical conditions, carbon depleted models produce less Ni-56. Comparing different models with similar Ni-56 yields shows lower kinetic energies in the ejecta for carbon depleted models, but only small differences in velocity distributions and line velocities in spectra. The light curve width-luminosity relation (WLR) obtained for models with differing carbon depletion is roughly perpendicular to the observed WLR, hence the carbon mass fraction is probably only a secondary parameter in the family of SNe Ia.
  • Type-Ia supernovae result from binary systems that include a carbon-oxygen white dwarf, and these thermonuclear explosions typically produce 0.5 M_solar of radioactive 56Ni. The 56Ni is commonly believed to be buried deeply in the expanding supernova cloud. Surprisingly, in SN2014J we detected the lines at 158 and 812 keV from 56Ni decay ({\tau}~8.8 days) earlier than the expected several-week time scale, only ~20 days after the explosion, and with flux levels corresponding to roughly 10% of the total expected amount of 56Ni. Some mechanism must break the spherical symmetry of the supernova, and at the same time create a major amount of 56Ni at the outskirts. A plausible explanation is that a belt of helium from the companion star is accreted by the white dwarf, where this material explodes and then triggers the supernova event.
  • In the version of the SD scenario of SNe Ia studied here, a CO WD explodes close to the Chandrasekhar limit after accreting material from a non-degenerate He companion. In the present study, we employ the Stellar GADGET code to perform 3D hydrodynamical simulations of the interaction of the SN Ia ejecta with the He companion taking into account its orbital motion and spin. It is found that only 2%--5% of the initial companion mass are stripped off from the outer layers of He companions due to the SN impact. The dependence of the unbound mass (or the kick velocity) on the orbital separation can be fitted in good approximation by a power law for a given companion model. After the SN impact, the outer layers of a He donor star are significantly enriched with heavy elements from the low-expansion-velocity tail of SN Ia ejecta. The total mass of accumulated SN-ejecta material on the companion surface reaches about > 10e-3 M_sun for different companion models. This enrichment with heavy elements provides a potential way to observationally identify the surviving companion star in SN remnants. Finally, by artificially adjusting the explosion energy of the W7 explosion model, we find that the total accumulation of SN ejecta on the companion surface is also dependent on the explosion energy with a power law relation in good approximation.
  • 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-200d. 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 characterised by weak or absent [Fe III] emission, pointing at a low ejecta ionisation state as a result of high densities. To constrain the ejecta and 56Ni 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 ~60d. Models with enough 56Ni 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 gamma-ray trapping. We instead propose a model with ~1 solar mass of 56Ni and ~2 solar masses of ejecta, which may be interpreted as the explosion of a Chandrasekhar-mass white dwarf (WD) enshrouded by 0.6-0.7 solar masses 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.
  • We present results for a suite of fourteen three-dimensional, high resolution hydrodynamical simulations of delayed-detonation modelsof 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 database for Chandrasekhar-mass delayed-detonation model nucleosynthetic yields and for deriving synthetic observables such as spectra and light curves. We employ a physically motivated, stochastic model based on turbulent velocity fluctuations and fuel density to calculate in situ the deflagration to detonation transition (DDT) 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 central density of 2.9 x 10^9 gcc, plus in addition one high central density (5.5 x 10^9 gcc), and one low central density (1.0 x 10^9 gcc) rendition of the 100 ignition kernel configuration. For each simulation we determined detailed nucleosynthetic yields by post-processing 10^6 tracer particles with a 384 nuclide reaction network. All delayed detonation models result in explosions unbinding the white dwarf, producing a range of 56Ni masses from 0.32 to 1.11 solar masses. As a general trend, the models predict that the stable neutron-rich iron group isotopes are not found at the lowest velocities, but rather at intermediate velocities (~3,000 - 10,000 km/s) in a shell surrounding a 56Ni-rich core. The models further predict relatively low velocity oxygen and carbon, with typical minimum velocities around 4,000 and 10,000 km/s, respectively.
  • 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 colour 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 SiII 6355A 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.
  • We present the first maximum-light ultraviolet (UV) through near-infrared (NIR) Type Ia supernova (SN Ia) spectrum. This spectrum of SN 2011iv was obtained nearly simultaneously by the Hubble Space Telescope at UV/optical wavelengths and the Magellan Baade telescope at NIR wavelengths. These data provide the opportunity to examine the entire maximum-light SN Ia spectral-energy distribution. Since the UV region of a SN Ia spectrum is extremely sensitive to the composition of the outer layers of the explosion, which are transparent at longer wavelengths, this unprecedented spectrum can provide strong constraints on the composition of the SN ejecta, and similarly the SN explosion and progenitor system. SN 2011iv is spectroscopically normal, but has a relatively fast decline (Delta m_15 (B) = 1.69 +/- 0.05 mag). We compare SN 2011iv to other SNe Ia with UV spectra near maximum light and examine trends between UV spectral properties, light-curve shape, and ejecta velocity. We tentatively find that SNe with similar light-curve shapes but different ejecta velocities have similar UV spectra, while those with similar ejecta velocities but different light-curve shapes have very different UV spectra. Through a comparison with explosion models, we find that both a solar-metallicity W7 and a zero-metallicity delayed-detonation model provide a reasonable fit to the spectrum of SN 2011iv from the UV to the NIR.
  • Xiaofeng Wang, Lifan Wang, Alexei V. Filippenko, Eddie Baron, Markus Kromer, Dennis Jack, Tianmeng Zhang, Greg Aldering, Pierre Antilogus, David Arnett, Dietrich Baade, Brian J. Barris, Stefano Benetti, Patrice Bouchet, Adam S. Burrows, Ramon Canal, Enrico Cappellaro, Raymond Carlberg, Elisa di Carlo, Peter Challis, Arlin Crotts, John I. Danziger, Massimo Della Valle, Michael Fink, Ryan J. Foley, Claes Fransson, Avishay Gal-Yam, Peter Garnavich, Chris L. Gerardy, Gerson Goldhaber, Mario Hamuy, Wolfgang Hillebrandt, Peter A. Hoeflich, Stephen T. Holland, Daniel E. Holz, John P. Hughes, David J. Jeffery, Saurabh W. Jha, Dan Kasen, Alexei M. Khokhlov, Robert P. Kirshner, Robert Knop, Cecilia Kozma, Kevin Krisciunas, Brian C. Lee, Bruno Leibundgut, Eric J. Lentz, Douglas C. Leonard, Walter H. G. Lewin, Weidong Li, Mario Livio, Peter Lundqvist, Dan Maoz, Thomas Matheson, Paolo Mazzali, Peter Meikle, Gajus Miknaitis, Peter Milne, Stefan Mochnacki, Ken'Ichi Nomoto, Peter E. Nugent, Elaine Oran, Nino Panagia, Saul Perlmutter, Mark M. Phillips, Philip Pinto, Dovi Poznanski, Christopher J. Pritchet, Martin Reinecke, Adam Riess, Pilar Ruiz-Lapuente, Richard Scalzo, Eric M. Schlegel, Brian Schmidt, James Siegrist, Alicia M. Soderberg, Jesper Sollerman, George Sonneborn, Anthony Spadafora, Jason Spyromilio, Richard A. Sramek, Sumner G. Starrfield, Louis G. Strolger, Nicholas B. Suntzeff, Rollin Thomas, John L. Tonry, Amedeo Tornambe, James W. Truran, Massimo Turatto, Michael Turner, Schuyler D. Van Dyk, Kurt Weiler, J. Craig Wheeler, Michael Wood-Vasey, Stan Woosley, Hitoshi Yamaoka
    Feb. 6, 2012 astro-ph.CO, astro-ph.HE
    We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects.
  • Type Ia supernovae (SNe Ia) are thought to result from thermonuclear explosions of carbon-oxygen white dwarf stars. Existing models generally explain the observed properties, with the exception of the sub-luminous 1991-bg-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 an underluminous explosion, though at the expense of requiring a single common-envelope phase, and component masses of ~0.9 M_sun. 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. While mass ratios can be slightly less than one and still produce an underluminous event, the masses have to be in the range 0.83-0.9 M_sun.