3D Models for High Velocity Features in Type Ia Supernovae

Spectral synthesis in 3-dimensional (3D) space for the earliest spectra of Type Ia supernovae (SNe Ia) is presented. In particular, the high velocity absorption features that are commonly seen at the earliest epochs ($\sim 10$ days before maximum light) are investigated by means of a 3D Monte Carlo spectral synthesis code. The increasing number of early spectra available allows statistical study of the geometry of the ejecta. The observed diversity in strength of the high velocity features (HVFs) can be explained in terms of a ``covering factor'', which represents the fraction of the projected photosphere that is concealed by high velocity material. Various geometrical models involving high velocity material with a clumpy structure or a thick torus can naturally account for the observed statistics of HVFs. HVFs may be formed by a combination of density and abundance enhancements. Such enhancements may be produced in the explosion itself or may be the result of interaction with circumstellar material or an accretion disk. Models with 1 or 2 blobs, as well as a thin torus or disk-like enhancement are unlikely as a standard situation.Comment: 17 pages, 12 figures. Accepted for publication in the Astrophysical Journa

Multi-Dimensional Simulations for Early Phase Spectra of Aspherical Hypernovae: SN 1998bw and Off-Axis Hypernovae

Early phase optical spectra of aspherical jet-like supernovae (SNe) are presented. We focus on energetic core-collapse SNe, or hypernovae. Based on hydrodynamic and nucleosynthetic models, radiative transfer in SN atmosphere is solved with a multi-dimensional Monte-Carlo radiative transfer code, SAMURAI. Since the luminosity is boosted in the jet direction, the temperature there is higher than in the equatorial plane by ~ 2,000 K. This causes anisotropic ionization in the ejecta. Emergent spectra are different depending on viewing angle, reflecting both aspherical abundance distribution and anisotropic ionization. Spectra computed with an aspherical explosion model with kinetic energy 20 x 10^{51} ergs are compatible with those of the Type Ic SN 1998bw if ~ 10-20% of the synthesized metals are mixed out to higher velocities. The simulations enable us to predict the properties of off-axis hypernovae. Even if an aspherical hypernova explosion is observed from the side, it should show hypernova-like spectra but with some differences in the line velocity, the width of the Fe absorptions and the strength of the Na I line.Comment: 4 pages, 4 figures. Accepted for publication in The Astrophysical Journal Letter

Abundance stratification in Type Ia Supernovae - II: The rapidly declining, spectroscopically normal SN 2004eo

The variation of properties of Type Ia supernovae, the thermonuclear explosions of Chandrasekhar-mass carbon-oxygen white dwarfs, is caused by different nucleosynthetic outcomes of these explosions, which can be traced from the distribution of abundances in the ejecta. The composition stratification of the spectroscopically normal but rapidly declining SN2004eo is studied performing spectrum synthesis of a time-series of spectra obtained before and after maximum, and of one nebular spectrum obtained about eight months later. Early-time spectra indicate that the outer ejecta are dominated by oxygen and silicon, and contain other intermediate-mass elements (IME), implying that the outer part of the star was subject only to partial burning. In the inner part, nuclear statistical equilibrium (NSE) material dominates, but the production of 56Ni was limited to ~0.43 \pm 0.05 Msun. An innermost zone containing ~0.25 Msun of stable Fe-group material is also present. The relatively small amount of NSE material synthesised by SN2004eo explains both the dimness and the rapidly evolving light curve of this SN.Comment: 12 pages, 7 figures. Accepted for publication in MNRA

On the gamma-ray emission of Type Ia Supernovae

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

The Outermost Ejecta of Type Ia Supernovae

The properties of the highest velocity ejecta of normal Type Ia supernovae (SNe Ia) are studied via models of very early optical spectra of 6 SNe. At epochs earlier than 1 week before maximum, SNe with a rapidly evolving Si II 6355 line velocity (HVG) have a larger photospheric velocity than SNe with a slowly evolving Si II 6355 line velocity (LVG). Since the two groups have comparable luminosities, the temperature at the photosphere is higher in LVG SNe. This explains the different overall spectral appearance of HVG and LVG SNe. However, the variation of the Ca II and Si II absorptions at the highest velocities (v >~ 20,000 km/s) suggests that additional factors, such as asphericity or different abundances in the progenitor white dwarf, affect the outermost layers. The C II 6578 line is marginally detected in 3 LVG SNe, suggesting that LVG undergo less intense burning. The carbon mass fraction is small, only less than 0.01 near the photosphere, so that he mass of unburned C is only <~ 0.01 Msun. Radioactive 56Ni and stable Fe are detected in both LVG and HVG SNe. Different Fe-group abundances in the outer layers may be one of the reasons for spectral diversity among SNe Ia at the earliest times. The diversity among SNe Ia at the earliest phases could also indicate an intrinsic dispersion in the width-luminosity relation of the light curve.Comment: 13 pages, 10 figures, Accepted for publication in The Astrophysical Journa

Optical Emission from Aspherical Supernovae and the Hypernova SN 1998bw

A fully 3D Monte Carlo scheme is applied to compute optical bolometric light curves for aspherical (jet-like) supernova explosion models. Density and abundance distributions are taken from hydrodynamic explosion models, with the energy varied as a parameter to explore the dependence. Our models show initially a very large degree ($\sim 4$ depending on model parameters) of boosting luminosity toward the polar ($z$) direction relative to the equatorial ($r$) plane, which decreases as the time of peak is approached. After the peak, the factor of the luminosity boost remains almost constant ($\sim 1.2$) until the supernova enters the nebular phase. This behavior is due mostly to the aspherical $^{56}$Ni distribution in the earlier phase and to the disk-like inner low-velocity structure in the later phase. Also the aspherical models yield an earlier peak date than the spherical models, especially if viewed from near the z-axis. Aspherical models with ejecta mass \sim 10\Msun are examined, and one with the kinetic energy of the expansion $\sim 2 \pm 0.5 \times 10^{52}$ ergs and a mass of $^{56}$Ni \sim 0.4\Msun yields a light curve in agreement with the observed light curve of SN 1998bw (the prototypical hyper-energetic supernova). The aspherical model is also at least qualitatively consistent with evolution of photospheric velocities, showing large velocities near the z-axis, and with a late-phase nebular spectrum. The viewing angle is close to the z-axis, strengthening the case for the association of SN 1998bw with the gamma ray burst GRB980425.Comment: Accepted by the Astrophysical Journal. 28 pages, 14 figure

Abundance stratification in Type Ia supernovae - V. SN 1986G bridging the gap between normal and subluminous SNe Ia

A detailed spectroscopic analysis of SN 1986G has been performed. SN 1986G `bridges the gap' between normal and sub luminous type Ia supernova (SNe Ia). The abundance tomography technique is used to determine the abundance distribution of the elements in the ejecta. SN 1986G was found to be a low energy Chandrasekhar mass explosion. Its kinetic energy was 70% of the standard W7 model (0.9x10$^{51}$erg). Oxygen dominates the ejecta from the outermost layers down to $\sim$9000kms$^{-1}$ , intermediate mass elements (IME) dominate from $\sim$ 9000kms$^{-1}$ to $\sim$ 3500kms$^{-1}$ with Ni and Fe dominating the inner layers $<\sim$ 3500kms$^{-1}$. The final masses of the main elements in the ejecta were found to be, O=0.33M, IME=0.69M, stable NSE=0.21M, $^{56}$Ni=0.14M. An upper limit of the carbon mass is set at C=0.02M. The spectra of SN1986G consist of almost exclusively singly ionised species. SN1986G can be thought of as a low luminosity extension of the main population of SN Ia, with a large deflagration phase that produced more IMEs than a standard SN Ia.Comment: Accepted for publication in MNRAS, update

The Type Ic Hypernova SN 2003dh/GRB 030329

The spectra of SN 2003dh, identified in the afterglow of GRB030329, are modeled using radiation transport codes. It is shown that SN 2003dh had a high explosion kinetic energy ($\sim 4 \times 10^{52}$ erg in spherical symmetry), making it one of the most powerful hypernovae observed so far, and supporting the case for association between hypernovae and Gamma Ray Bursts. However, the light curve derived from fitting the spectra suggests that SN 2003dh was not as bright as SN 1998bw, ejecting only \sim 0.35\Msun of \Nifs. The spectra of SN 2003dh resemble those of SN 1998bw around maximum, but later they look more like those of the less energetic hypernova SN 1997ef. The spectra and the inferred light curve can be modeled adopting a density distribution similar to that used for SN 1998bw at $v > 25,000$\kms but more like that of SN 1997ef at lower velocities. The mass of the ejecta is \sim 8\Msun, somewhat less than in the other two hypernovae. The progenitor must have been a massive star (M \sim 35-40\Msun), as for other hypernovae. The need to combine different one-dimensional explosion models strongly indicates that SN 2003dh was an asymmetric explosion.Comment: 11 pages, 1 table and 5 figures. To appear in the Astrophysical Journal (Letters). Revised version taking referee's comments into account, minor change

Can differences in the nickel abundance in Chandrasekhar mass models explain the relation between brightness and decline rate of normal Type Ia Supernovae?

The use of Type Ia supernovae as distance indicators relies on the determination of their brightness. This is not constant, but it can be calibrated using an observed relation between the brightness and the properties of the optical light curve (decline rate, width, shape), which indicates that brighter SNe have broader, slower light curves. However, the physical basis for this relation is not yet fully understood. Among possible causes are different masses of the progenitor white dwarfs or different opacities in Chandrasekhar-mass explosions. We parametrise the Chandrasekhar-mass models presented by Iwamoto et al (1999), which synthesize different amounts of Ni, and compute bolometric light curves and spectra at various epochs. Since opacity in SNe Ia is due mostly to spectral lines, it should depend on the mass of Fe-peak elements synthesized in the explosion, and on the temperature in the ejecta. Bolometric light curves computed using these prescriptions for the optical opacity reproduce the relation between brightness and decline rate. Furthermore, when spectra are calculated, the change in colour between maximum and two weeks later allows the observed relation between M_B(Max) and Dm_{15}(B) to be reproduced quite nicely. Spectra computed at various epochs compare well with corresponding spectra of spectroscopically normal SNeIa selected to cover a similar range of Dm_{15}(B) values.Comment: 25 pages, including 6 figures. Accepted for publication in Ap