Gamma-ray diagnostics of Type Ia supernovae. Predictions of observables from three-dimensional modeling

Context. Although the question of progenitor systems and detailed explosion mechanisms still remains a matter of discussion, it is commonly believed that Type Ia supernovae (SNe Ia) are production sites of large amounts of radioactive nuclei. Even though the gamma-ray emission due to radioactive decays is responsible for powering the light curves of SNe Ia, gamma rays themselves are of particular interest as a diagnostic tool because they directly lead to deeper insight into the nucleosynthesis and the kinematics of these explosion events. Aims: We study the evolution of gamma-ray line and continuum emission of SNe Ia with the objective of analyzing the relevance of observations in this energy range. We seek to investigate the chances for the success of future MeV missions regarding their capabilities for constraining the intrinsic properties and the physical processes of SNe Ia. Methods: Focusing on two of the most broadly discussed SN Ia progenitor scenarios - a delayed detonation in a Chandrasekhar-mass white dwarf (WD) and a violent merger of two WDs - we used three-dimensional explosion models and performed radiative transfer simulations to obtain synthetic gamma-ray spectra. Both chosen models produce the same mass of 56Ni and have similar optical properties that are in reasonable agreement with the recently observed supernova SN 2011fe. We examine the gamma-ray spectra with respect to their distinct features and draw connections to certain characteristics of the explosion models. Applying diagnostics, such as line and hardness ratios, the detection prospects for future gamma-ray missions with higher sensitivities in the MeV energy range are discussed. Results: In contrast to the optical regime, the gamma-ray emission of our two chosen models proves to be quite different. The almost direct connection of the emission of gamma rays to fundamental physical processes occurring in SNe Ia permits additional constraints concerning several explosion model properties that are not easily accessible within other wavelength ranges. Proposed future MeV missions such as GRIPS will resolve all spectral details only for nearby SNe Ia, but hardness ratio and light curve measurements still allow for a distinction of the two different models at 10 Mpc and 16 Mpc for an exposure time of 106 s. The possibility of detecting the strongest line features up to the Virgo distance will offer the opportunity to build up a first sample of SN Ia detections in the gamma-ray energy range and underlines the importance of future space observatories for MeV gamma rays.

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.

Type Ia Supernovae : Homogeneity and Diversity

The use of type Ia supernovae as distance estimators has shown that about 75% of the energy content of the universe has a negative equation of state parameter and thus, drives the acceleration of the universe. Constraining the exact nature of this energy is one of the main goals in cosmology. As the statistics of observed high-redshift supernovae increases, systematic effects become the limiting factor to pursue such investigations, thus deeper understanding of the physical properties of SNe is of great importance. In this thesis we investigate spectral homogeneity and diversity of local and high redshift supernovae. Special emphasis has been given to the analysis of optical spectra of local peculiar supernovae 1999aa and 1999ac. The study of the spectra of SN 1999aa pointed out that this SN could be a link between the extreme peculiar SN 1991T and normal SNe. Moreover, the identification of a high velocity component of Ca II and possibly of a low velocity component of C III suggests some degree of asphericity in the ejecta of this supernova. Evidence for a deflagration of a C+O white dwarf was found in the early spectra of SN 1999ac. The spectral proprieties of a vast sample of local SNe are also studied by means of newly introduced spectral indicators. These were used to possibly improve the intrinsic spread of SN peak magnitudes to 0.15 mag, independently of light curve parameters. The first quantitative comparison between local and high redshift supernova is carried out. No evidence for extreme peculiar sub-luminous SNe was found in our data set including 13 SNe with redshift range z=0.279-0.912. Furthermore, SN2002fd (z=0.279) was found to show spectral characteristics similar to SN 1991T/SN 1999aa-like supernovae. We also present a feasibility study of the Hubble diagram in rest frame I-band up to z~0.5, and show the possibility to probe the presence of intergalactic dust, which could possibly mimic the effect of dark energy in the Hubble diagram.

Type Ia supernovae from violent mergers of carbon-oxygen white dwarfs: polarization signatures

The violent merger of two carbon-oxygen white dwarfs has been proposed as a viable progenitor for some Type Ia supernovae. However, it has been argued that the strong ejecta asymmetries produced by this model might be inconsistent with the low degree of polarization typically observed in Type Ia supernova explosions. Here, we test this claim by carrying out a spectropolarimetric analysis for the model proposed by Pakmor et al. for an explosion triggered during the merger of a 1.1 and 0.9 M⊙ carbon-oxygen white dwarf binary system. Owing to the asymmetries of the ejecta, the polarization signal varies significantly with viewing angle. We find that polarization levels for observers in the equatorial plane are modest (≲1 per cent) and show clear evidence for a dominant axis, as a consequence of the ejecta symmetry about the orbital plane. In contrast, orientations out of the plane are associated with higher degrees of polarization and departures from a dominant axis. While the particular model studied here gives a good match to highly polarized events such as SN 2004dt, it has difficulties in reproducing the low polarization levels commonly observed in normal Type Ia supernovae. Specifically, we find that significant asymmetries in the element distribution result in a wealth of strong polarization features that are not observed in the majority of currently available spectropolarimetric data of Type Ia supernovae. Future studies will map out the parameter space of the merger scenario to investigate if alternative models can provide better agreement with observations.

Searching for swept-up hydrogen and helium in the late-time spectra of 11 nearby Type Ia supernovae

The direct detection of a stellar system that explodes as a Type Ia supernova (SN Ia) has not yet been successful. Various indirect methods have been used to investigate SN Ia progenitor systems but none have produced conclusive results. A prediction of single-degenerate models is that H- (or He-) rich material from the envelope of the companion star should be swept up by the SN ejecta in the explosion. Seven SNe Ia have been analysed to date looking for signs of H-rich material in their late-time spectra and none were detected. We present results from new late-time spectra of 11 SNe Ia obtained at the Very Large Telescope using XShooter and FORS2. We present the tentative detection of Hα emission for SN 2013ct, corresponding to ∼0.007 M_⊙ of stripped/ablated companion star material (under the assumptions of the spectral modelling). This mass is significantly lower than expected for single-degenerate scenarios, suggesting that >0.1 M_⊙ of H-rich is present but not observed. We do not detect Hα emission in the other 10 SNe Ia. This brings the total sample of normal SNe Ia with non-detections (<0.001–0.058 M_⊙) of H-rich material to 17 events. The simplest explanation for these non-detections is that these objects did not result from the explosion of a CO white dwarf accreting matter from a H-rich companion star via Roche lobe overflow or symbiotic channels. However, further spectral modelling is needed to confirm this. We also find no evidence of He-emission features, but models with He-rich companion stars are not available to place mass limits.

Predicting polarization signatures for double-detonation and delayed-detonation models of Type Ia supernovae

Calculations of synthetic spectropolarimetry are one means to test multidimensional explosion models for Type Ia supernovae. In a recent paper, we demonstrated that the violent merger of a 1.1 and 0.9 M⊙ white dwarf binary system is too asymmetric to explain the low polarization levels commonly observed in normal Type Ia supernovae. Here, we present polarization simulations for two alternative scenarios: the sub-Chandrasekhar mass double-detonation and the Chandrasekhar mass delayed-detonation model. Specifically, we study a 2D double-detonation model and a 3D delayed-detonation model, and calculate polarization spectra for multiple observer orientations in both cases. We find modest polarization levels (<1 per cent) for both explosion models. Polarization in the continuum peaks at ∼0.1–0.3 per cent and decreases after maximum light, in excellent agreement with spectropolarimetric data of normal Type Ia supernovae. Higher degrees of polarization are found across individual spectral lines. In particular, the synthetic Si II λ6355 profiles are polarized at levels that match remarkably well the values observed in normal Type Ia supernovae, while the low degrees of polarization predicted across the O I λ7774 region are consistent with the non-detection of this feature in current data. We conclude that our models can reproduce many of the characteristics of both flux and polarization spectra for well-studied Type Ia supernovae, such as SN 2001el and SN 2012fr. However, the two models considered here cannot account for the unusually high level of polarization observed in extreme cases such as SN 2004dt.

New Mass Limit for White Dwarfs: Super-Chandrasekhar Type Ia Supernova as a New Standard Candle

Type Ia supernovae, sparked off by exploding white dwarfs of mass close to the Chandrasekhar limit, play the key role in understanding the expansion rate of the Universe. However, recent observations of several peculiar type Ia supernovae argue for its progenitor mass to be significantly super-Chandrasekhar. We show that strongly magnetized white dwarfs not only can violate the Chandrasekhar mass limit significantly, but exhibit a different mass limit. We establish from a foundational level that the generic mass limit of white dwarfs is 2.58 solar mass. This explains the origin of overluminous peculiar type Ia supernovae. Our finding further argues for a possible second standard candle, which has many far reaching implications, including a possible reconsideration of the expansion history of the Universe. DOI: 10.1103/PhysRevLett.110.071102

On the nature of the circumstellar medium of the remarkable Type Ia/IIn supernova SN 2002ic

We present results from the first high-resolution, high signal-to-noise ratio spectrum of SN 2002ic. The resolved Ha line has a P Cygni-type profile, clearly demonstrating the presence of a dense, slow-moving (~100 km s-1) outflow. We have additionally found a huge near-infrared excess, hitherto unseen in Type Ia supernovae. We argue that this is due to an infrared light-echo arising from the pre-existing dusty circumstellar medium. We deduce a circumstellar medium mass probably exceeding 0.3 Msolar produced by a mass-loss rate greater than several times 10-4 Msolar yr-1. For the progenitor, we favour a single-degenerate system where the companion is a post-asymptotic giant branch star. As a by-product of our optical data, we are able to provide a firm identification of the host galaxy of SN 2002ic.

Signatures of delayed detonation, asymmetry, and electron capture in the mid-infrared spectra of supernovae 2003hv and 2005df

We present mid-infrared (5.2-15.2 mu m) spectra of the Type Ia supernovae (SNe Ia) 2003hv and 2005df observed with the Spitzer Space Telescope. These are the first observed mid-infrared spectra of thermonuclear supernovae, and show strong emission from fine-structure lines of Ni, Co, S, and Ar. The detection of Ni emission in SN 2005df 135 days after the explosion provides direct observational evidence of high-density nuclear burning forming a significant amount of stable Ni in a SN Ia. The SN 2005df Ar lines also exhibit a two-pronged emission profile, implying that the Ar emission deviates significantly from spherical symmetry. The spectrum of SN 2003hv also shows signs of asymmetry, exhibiting blueshifted [Co (III)], which matches the blueshift of [Fe (II)] lines in nearly coeval near-infrared spectra. Finally, local thermodynamic equilibrium abundance estimates for the yield of radioactive Ni-56 give M-56Ni approximate to 0.5 M-circle dot, for SN 2003hv, but only M-56Ni approximate to 0.13-0.22 M-circle dot for the apparently subluminous SN 2005df, supporting the notion that the luminosity of SNe Ia is primarily a function of the radioactive 56Ni yield. The observed emission-line profiles in the SN 2005df spectrum indicate a chemically stratified ejecta structure, which matches the predictions of delayed detonation (DD) models, but is entirely incompatible with current three-dimensional deflagration models. Furthermore, the degree that this layering persists to the innermost regions of the supernova is difficult to explain even in a DD scenario, where the innermost ejecta are still the product of deflagration burning. Thus, while these results are roughly consistent with a delayed detonation, it is clear that a key piece of physics is still missing from our understanding of the earliest phases of SN Ia explosions.

An early-time infrared and optical study of the Type Ia Supernova 1998bu in M96

We present first-season infrared (IR) and optical photometry and spectroscopy of the Type Ia Supernova 1998bu in M96. We also report optical polarimetry of this event. SN 1998bu is one of the closest type Ia supernovae of modern times, and the distance of its host galaxy is well determined. We find that SN 1998bu is both photometrically and spectroscopically normal. However, the extinction to this event is unusually high, with A(V) = 1.0 +/- 0.11. We find that SN 1998bu peaked at an intrinsic M-V = -19.37 +/- 0.23. Adopting a distance modulus of 30.25 (Tanvir et al.) and using Phillips et al.'s relations for the Hubble constant, we obtain H-0 = 70.4 +/- 4.3 km s(-1) Mpc(-1). Combination of our IR photometry with those of Jha et al. provides one of the most complete early-phase IR light curves for a SN Ia published so far. In particular, SN 1998bu is the first normal SN Ia for which good pre-t(Bmax) IR coverage has been obtained. It reveals that the J, H and K light curves peak about 5 days earlier than the flux in the B-band curve.

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.

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.

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.

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.