Modeling Type Ia supernova explosions

Despite their astrophysical significanceas a major contributor to cosmic nucleosynthesis and as distance indicators in observational cosmologyType Ia supernovae lack theoretical explanation. Not only is the explosion mechanism complex due to the interaction of (potentially turbulent) hydrodynamics and nuclear reactions, but even the initial conditions for the explosion are unknown. Various progenitor scenarios have been proposed. After summarizing some general aspects of Type Ia supernova modeling, recent simulations of our group are discussed. With a sequence of modeling starting (in some cases) from the progenitor evolution and following the explosion hydrodynamics and nucleosynthesis we connect to the formation of the observables through radiation transport in the ejecta cloud. This allows us to analyze several models and to compare their outcomes with observations. While pure deflagrations of Chandrasekhar-mass white dwarfs and violent mergers of two white dwarfs lead to peculiar events (that may, however, find their correspondence in the observed sample of SNe Ia), only delayed detonations in Chandrasekhar-mass white dwarfs or sub-Chandrasekhar-mass explosions remain promising candidates for explaining normal Type Ia supernovae. © 2011 Elsevier B.V. All rights reserved.

CI Aql: a Type Ia supernova progenitor?

If recurrent novae are progenitors of Type Ia supernovae, their white dwarfs must have masses close to the Chandrasekhar limit. The most reliable means of determining white dwarf masses in recurrent novae is dynamically, via radial-velocity and rotational-broadening measurements of the companion star. Such measurements require the system to be both eclipsing and to show absorption features from the secondary star. Prior to the work reported here, the only dynamical mass estimate of a recurrent nova was for U Sco, which has a white dwarf mass of 1.55 +/- 0.24 Msolar (Thoroughgood et al. 2001). We present new time-resolved, intermediate-resolution spectroscopy of the eclipsing recurrent nova CI Aquilae (CI Aql) during quiescence. We find the mass of the white dwarf to be 1.00 +/- 0.14 Msolar and the mass of the secondary star to be 2.32 +/- 0.19 Msolar. We estimate the radius of the secondary to be 2.07 +/- 0.06 Rsolar, implying that it is a slightly-evolved early A-type star. The high mass ratio of q = 2.35 +/- 0.24 and the high secondary-star mass implies that the mass transfer occurs on a thermal timescale. We suggest that CI Aql is rapidly evolving into a supersoft X-ray source, and ultimately may explode as a Type Ia supernova within 10 Myr.

Spectroscopic Observations of SN 2012fr: A Luminous Normal Type Ia Supernova with Early High Velocity Features and Late Velocity Plateau

We present 65 optical spectra of the Type Ia supernova SN 2012fr, of which 33 were obtained before maximum light. At early times SN 2012fr shows clear evidence of a high-velocity feature (HVF) in the Si II 6355 line which can be cleanly decoupled from the lower velocity "photospheric" component. This Si II 6355 HVF fades by phase -5; subsequently, the photospheric component exhibits a very narrow velocity width and remains at a nearly constant velocity of v~12,000 km/s until at least 5 weeks after maximum brightness. The Ca II infrared (IR) triplet exhibits similar evidence for both a photospheric component at v~12,000 km/s with narrow line width and long velocity plateau, as well as a high-velocity component beginning at v~31,000 km/s two weeks before maximum. SN 2012fr resides on the border between the "shallow silicon" and "core-normal" subclasses in the Branch et al. (2009) classification scheme, and on the border between normal and "high-velocity" SNe Ia in the Wang et al. (2009a) system. Though it is a clear member of the "low velocity gradient" (LVG; Benetii et al., 2005) group of SNe Ia and exhibits a very slow light-curve decline, it shows key dissimilarities with the overluminous SN 1991T or SN 1999aa subclasses of SNe Ia. SN 2012fr represents a well-observed SN Ia at the luminous end of the normal SN Ia distribution, and a key transitional event between nominal spectroscopic subclasses of SNe Ia.

Type Ia supernova bolometric light curves and ejected mass estimates from the nearby supernova factory

We present a sample of normal Type Ia supernovae (SNe Ia) from the Nearby Supernova Factory data set with spectrophotometry at sufficiently late phases to estimate the ejected mass using the bolometric light curve.Wemeasure Ni masses from the peak bolometric luminosity, then compare the luminosity in the Co-decay tail to the expected rate of radioactive energy release from ejecta of a given mass. We infer the ejected mass in a Bayesian context using a semi-analytic model of the ejecta, incorporating constraints from contemporary numerical models as priors on the density structure and distribution of Ni throughout the ejecta. We find a strong correlation between ejected mass and light-curve decline rate, and consequently Ni mass, with ejected masses in our data ranging from 0.9 to 1.4 M. Most fast-declining (SALT2 x <-1) normal SNe Ia have significantly sub-Chandrasekhar ejected masses in our fiducial analysis.

SN 2010LP - A type IA supernova from a violent merger of two carbon-oxygen white dwarfs

SN 2010lp is a subluminous Type Ia supernova (SN Ia) with slowly evolving lightcurves. Moreover, it is the only subluminous SN Ia observed so far that shows narrow emission lines of [O I] in late-time spectra, indicating unburned oxygen close to the center of the ejecta. Most explosion models for SNe Ia cannot explain the narrow [O I] emission. Here, we present hydrodynamic explosion and radiative transfer calculations showing that the violent merger of two carbon-oxygen white dwarfs of 0.9 and 0.76 M adequately reproduces the early-time observables of SN 2010lp. Moreover, our model predicts oxygen close to the center of the explosion ejecta, a pre-requisite for narrow [O I] emission in nebular spectra as observed in SN 2010lp. © 2013. The American Astronomical Society. All rights reserved.

Polarization spectral synthesis for Type Ia supernova explosion models

We present a Monte Carlo radiative transfer technique for calculating synthetic spectropolarimetry for multidimensional supernova explosion models. The approach utilizes 'virtual-packets' that are generated during the propagation of the Monte Carlo quanta and used to compute synthetic observables for specific observer orientations. Compared to extracting synthetic observables by direct binning of emergent Monte Carlo quanta, this virtual-packet approach leads to a substantial reduction in the Monte Carlo noise. This is not only vital for calculating synthetic spectropolarimetry (since the degree of polarization is typically very small) but also useful for calculations of light curves and spectra. We first validate our approach via application of an idealized test code to simple geometries. We then describe its implementation in the Monte Carlo radiative transfer code ARTIS and present test calculations for simple models for Type Ia supernovae. Specifically, we use the well-known one-dimensional W7 model to verify that our scheme can accurately recover zero polarization from a spherical model, and to demonstrate the reduction in Monte Carlo noise compared to a simple packet-binning approach. To investigate the impact of aspherical ejecta on the polarization spectra, we then use ARTIS to calculate synthetic observables for prolate and oblate ellipsoidal models with Type Ia supernova compositions.

Systematic Uncertainties Associated with the Cosmological Analysis of the First PAN-STARRS1 Type Ia Supernova Sample

We probe the systematic uncertainties from the 113 Type Ia supernovae (SN Ia) in the Pan-STARRS1 (PS1) sample along with 197 SN Ia from a combination of low-redshift surveys. The companion paper by Rest et al. describes the photometric measurements and cosmological inferences from the PS1 sample. The largest systematic uncertainty stems from the photometric calibration of the PS1 and low-z samples. We increase the sample of observed Calspec standards from 7 to 10 used to define the PS1 calibration system. The PS1 and SDSS-II calibration systems are compared and discrepancies up to ∼0.02 mag are recovered. We find uncertainties in the proper way to treat intrinsic colors and reddening produce differences in the recovered value of w up to 3%. We estimate masses of host galaxies of PS1 supernovae and detect an insignificant difference in distance residuals of the full sample of 0.037 ± 0.031 mag for host galaxies with high and low masses. Assuming flatness and including systematic uncertainties in our analysis of only SNe measurements, we find w = -1.120^+0.360_-0.206(Stat)^+0.269_-0.291(Sys). With additional constraints from Baryon acoustic oscillation, cosmic microwave background (CMB) (Planck) and H₀ measurements, we find w = -1.166^+0.072_-0.069 and Ω_m = 0.280^+0.013_-0.012 (statistical and systematic errors added in quadrature). The significance of the inconsistency with w = -1 depends on whether we use Planck or Wilkinson Microwave Anisotropy Probe measurements of the CMB: w_{BAO+H0+SN+WMAP} = -1.124^+0.083_-0.065.

Early ultraviolet emission in the Type Ia supernova LSQ12gdj: No evidence for ongoing shock interaction

We present photospheric-phase observations of LSQ12gdj, a slowly declining, UV-bright Type Ia supernova. Classified well before maximum light, LSQ12gdj has extinction-corrected absolute magnitude MB = -19.8, and pre-maximum spectroscopic evolution similar to SN 1991T and the super-Chandrasekhar-mass SN 2007if. We use ultraviolet photometry from Swift, ground-based optical photometry, and corrections from a near-infrared photometric template to construct the bolometric (1600-23 800 Å) light curve out to 45 d past B-band maximum light. We estimate that LSQ12gdj produced 0.96 ± 0.07 M· of 56Ni, with an ejected mass near or slightly above the Chandrasekhar mass. As much as 27 per cent of the flux at the earliest observed phases, and 17 per cent at maximum light, is emitted bluewards of 3300 Å. The absence of excess luminosity at late times, the cutoff of the spectral energy distribution bluewards of 3000 Å and the absence of narrow line emission and strong Na I D absorption all argue against a significant contribution from ongoing shock interaction. However, ~10 per cent of LSQ12gdj's luminosity near maximum light could be produced by the release of trapped radiation, including kinetic energy thermalized during a brief interaction with a compact, hydrogen-poor envelope (radius <1013 cm) shortly after explosion; such an envelope arises generically in double-degenerate merger scenarios.

5.9-keV Mn K-shell X-ray luminosity from the decay of 55Fe in Type Ia supernova models

We show that the X-ray line flux of the Mn Kα line at 5.9 keV from the decay of 55Fe is a promising diagnostic to distinguish between Type Ia supernova (SN Ia) explosion models. Using radiation transport calculations, we compute the line flux for two three-dimensional explosion models: a near-Chandrasekhar mass delayed detonation and a violent merger of two (1.1 and 0.9 M⊙) white dwarfs. Both models are based on solar metallicity zero-age main-sequence progenitors. Due to explosive nuclear burning at higher density, the delayed-detonation model synthesizes ˜3.5 times more radioactive 55Fe than the merger model. As a result, we find that the peak Mn Kα line flux of the delayed-detonation model exceeds that of the merger model by a factor of ˜4.5. Since in both models the 5.9-keV X-ray flux peaks five to six years after the explosion, a single measurement of the X-ray line emission at this time can place a constraint on the explosion physics that is complementary to those derived from earlier phase optical spectra or light curves. We perform detector simulations of current and future X-ray telescopes to investigate the possibilities of detecting the X-ray line at 5.9 keV. Of the currently existing telescopes, XMM-Newton/pn is the best instrument for close (≲1-2 Mpc), non-background limited SNe Ia because of its large effective area. Due to its low instrumental background, Chandra/ACIS is currently the best choice for SNe Ia at distances above ˜2 Mpc. For the delayed-detonation scenario, a line detection is feasible with Chandra up to ˜3 Mpc for an exposure time of 106 s. We find that it should be possible with currently existing X-ray instruments (with exposure times ≲5 × 105 s) to detect both of our models at sufficiently high S/N to distinguish between them for hypothetical events within the Local Group. The prospects for detection will be better with future missions. For example, the proposed Athena/X-IFU instrument could detect our delayed-detonation model out to a distance of ˜5 Mpc. This would make it possible to study future events occurring during its operational life at distances comparable to those of the recent supernovae SN 2011fe (˜6.4 Mpc) and SN 2014J (˜3.5 Mpc).

EARLY RADIO AND X-RAY OBSERVATIONS OF THE YOUNGEST NEARBY TYPE Ia SUPERNOVA PTF 11kly (SN 2011fe)

On 2011 August 24 (UT) the Palomar Transient Factory (PTF) discovered PTF11kly (SN 2011fe), the youngest and most nearby Type Ia supernova (SN Ia) in decades. We followed this event up in the radio (centimeter and millimeter bands) and X-ray bands, starting about a day after the estimated explosion time. We present our analysis of the radio and X-ray observations, yielding the tightest constraints yet placed on the pre-explosion mass-loss rate from the progenitor system of this supernova. We find a robust limit of from sensitive X-ray non-detections, as well as a similar limit from radio data, which depends, however, on assumptions about microphysical parameters. We discuss our results in the context of single-degenerate models for SNe Ia and find that our observations modestly disfavor symbiotic progenitor models involving a red giant donor, but cannot constrain systems accreting from main-sequence or sub-giant stars, including the popular supersoft channel. In view of the proximity of PTF11kly and the sensitivity of our prompt observations, we would have to wait for a long time (a decade or longer) in order to more meaningfully probe the circumstellar matter of SNe Ia.

Consequences of dark matter-dark energy interaction on cosmological parameters derived from type Ia supernova data

Models where the dark matter component of the Universe interacts with the dark energy field have been proposed as a solution to the cosmic coincidence problem, since in the attractor regime both dark energy and dark matter scale in the same way. In these models the mass of the cold dark matter particles is a function of the dark energy field responsible for the present acceleration of the Universe, and different scenarios can be parametrized by how the mass of the cold dark matter particles evolves with time. In this article we study the impact of a constant coupling delta between dark energy and dark matter on the determination of a redshift dependent dark energy equation of state w(DE)(z) and on the dark matter density today from SNIa data. We derive an analytical expression for the luminosity distance in this case. In particular, we show that the presence of such a coupling increases the tension between the cosmic microwave background data from the analysis of the shift parameter in models with constant w(DE) and SNIa data for realistic values of the present dark matter density fraction. Thus, an independent measurement of the present dark matter density can place constraints on models with interacting dark energy.

Constraining the dark energy and smoothness parameter with type Ia supernovae and gamma-ray bursts

The existence of inhomogeneities in the observed Universe modifies the distance-redshift relations thereby affecting the results of cosmological tests in comparison to the ones derived assuming spatially uniform models. By modeling the inhomogeneities through a Zeldovich-Kantowski-Dyer-Roeder approach which is phenomenologically characterized by a smoothness parameter alpha, we rediscuss the constraints on the cosmic parameters based on type Ia supernovae (SNe Ia) and gamma-ray bursts (GRBs) data. The present analysis is restricted to a flat Lambda CDM model with the reasonable assumption that Lambda does not clump. A chi(2) analysis using 557 SNe Ia data from the Union2 compilation data (R. Amanullah et al., Astrophys. J. 716, 712 (2010).) constrains the pair of parameters (Omega(m), alpha) to Omega(m) = 0.27(-0.03)(+0.08) (2 sigma) and alpha >= 0.25. A similar analysis based only on 59 Hymnium GRBs (H. Wei, J. Cosmol. Astropart. Phys. 08 (2010) 020.) constrains the matter density parameter to be Omega(m) = 0.35(-0.24)(+0.62) (2 sigma) while all values for the smoothness parameter are allowed. By performing a joint analysis, it is found that Omega(m) = 0.27(-0.06)(+0.06) and alpha >= 0.52. As a general result, although considering that current GRB data alone cannot constrain the smoothness alpha parameter, our analysis provides an interesting cosmological probe for dark energy even in the presence of inhomogeneities.

The peculiar Type Ia supernova iPTF14atg: Chandrasekhar-mass explosion or violent merger?

iPTF14atg, a subluminous peculiar Type Ia supernova (SN Ia) similar to SN 2002es, is the first SN Ia for which a strong UV flash was observed in the early-time light curves. This has been interpreted as evidence for a single-degenerate (SD) progenitor system, where such a signal is expected from interactions between the SN ejecta and the non-degenerate companion star. Here, we compare synthetic observables of multidimensional state-of-the-art explosion models for different progenitor scenarios to the light curves and spectra of iPTF14atg. From our models, we have difficulties explaining the spectral evolution of iPTF14atg within the SD progenitor channel. In contrast, we find that a violent merger of two carbon-oxygen white dwarfs with 0.9 and 0.76 M⊙, respectively, provides an excellent match to the spectral evolution of iPTF14atg from 10 d before to several weeks after maximum light. Our merger model does not naturally explain the initial UV flash of iPTF14atg. We discuss several possibilities like interactions of the SN ejecta with the circumstellar medium and surface radioactivity from an He-ignited merger that may be able to account for the early UV emission in violent merger models.

Optical and infrared observations of SN 2002dj: some possible common properties of fast-expanding Type Ia supernovae

As part of the European Supernova Collaboration, we obtained extensive photometry and spectroscopy of the Type Ia supernova (SN Ia) SN 2002dj covering epochs from 11 d before to nearly two years after maximum. Detailed optical and near-infrared observations show that this object belongs to the class of the high-velocity gradient events as indicated by Si, S and Ca lines. The light curve shape and velocity evolution of SN 2002dj appear to be nearly identical to SN 2002bo. The only significant difference is observed in the optical to near-infrared colours and a reduced spectral ernission beyond 6500 A. For high-velocity gradient SNe Ia, we tentatively identify a faster rise to maximum, a more pronounced inflection in the V and R light curves after maximum and a brighter, slower declining late-time B light curve as common photometric properties of this class of object. They also seem to be characterized by a different colour and colour evolution with respect to 'normal' SNe Ia. The usual light Curve shape parameters do not distinguish these events. Stronger, more blueshifted absorption features of intermediate-mass elements and lower temperatures are the most prominent spectroscopic features of SNe Ia displaying high-velocity gradients. It appears that these events burn more intermediate-mass elements in the outer layers. Possible connections to the metallicity of the progenitor star are explored.

ESC and KAIT observations of the transitional type Ia SN 2004eo

We present optical and infrared observations of the unusual Type Ia supernova (SN) 2004eo. The light curves and spectra closely resemble those of the prototypical SN 1992A, and the luminosity at maximum (M-B = -19.08) is close to the average for a Type Ia supernova (SN Ia). However, the ejected Ni-56 mass derived by modelling the bolometric light curve (about 0.45M(circle dot)) lies near the lower limit of the Ni-56 mass distribution observed in normal SNe Ia. Accordingly, SN 2004eo shows a relatively rapid post-maximum decline in the light curve [Delta m(15)(B)(true) = 1.46], small expansion velocities in the ejecta and a depth ratio Si II lambda 5972/ Si II lambda 6355 similar to that of SN 1992A. The physical properties of SN 2004eo cause it to fall very close to the boundary between the faint, low-velocity gradient and high-velocity gradient subgroups proposed by Benetti et al. Similar behaviour is seen in a few other SNe Ia. Thus, there may in fact exist a few SNe Ia with intermediate physical properties.