The Nearby Supernova Factory has made spectrophotometric observations of Type Ia supernovae since $2004$. This work presents an interim version of the data produced, including $210$ supernovae observed between $2004$ and $2013$.

There are now hundreds of publicly available supernova spectral time series. Radiative transfer modeling of this data gives insights into the physical properties of these explosions such as the composition, the density structure, or the intrinsic luminosity---this is invaluable for understanding the supernova progenitors, the explosion mechanism, or for constraining the supernova distance. However, a detailed parameter study of the available data has been out of reach due to the high dimensionality of the problem coupled with the still significant computational expense. We tackle this issue through the use of machine-learning emulators, which are algorithms for high-dimensional interpolation. These use a pre-calculated training dataset to mimic the output of a complex code but with run times orders of magnitude shorter. We present the application of such an emulator to synthetic type II supernova spectra generated with the TARDIS radiative transfer code. The results show that with a relatively small training set of 780 spectra we can generate emulated spectra with interpolation uncertainties of less than one percent. We demonstrate the utility of this method by automatic spectral fitting of two well-known type IIP supernovae; as an exemplary application, we determine the supernova distances from the spectral fits using the tailored-expanding-photosphere method. We compare our results to previous studies and find good agreement. This suggests that emulation of TARDIS spectra can likely be used to perform automatic and detailed analysis of many transient classes putting the analysis of large data repositories within reach.

Type Ia Supernovae (SNe Ia) are widely used to measure the expansion of the Universe. Improving distance measurements of SNe Ia is one technique to better constrain the acceleration of expansion and determine its physical nature. This document develops a new SNe Ia spectral energy distribution (SED) model, called the SUpernova Generator And Reconstructor (SUGAR), which improves the spectral description of SNe Ia, and consequently could improve the distance measurements. This model is constructed from SNe Ia spectral properties and spectrophotometric data from The Nearby Supernova Factory collaboration. In a first step, a PCA-like method is used on spectral features measured at maximum light, which allows us to extract the intrinsic properties of SNe Ia. Next, the intrinsic properties are used to extract the average extinction curve. Third, an interpolation using Gaussian Processes facilitates using data taken at different epochs during the lifetime of a SN Ia and then projecting the data on a fixed time grid. Finally, the three steps are combined to build the SED model as a function of time and wavelength. This is the SUGAR model. The main advancement in SUGAR is the addition of two additional parameters to characterize SNe Ia variability. The first is tied to the properties of SNe Ia ejecta velocity, the second is correlated with their calcium lines. The addition of these parameters, as well as the high quality the Nearby Supernova Factory data, makes SUGAR an accurate and efficient model for describing the spectra of normal SNe Ia as they brighten and fade. The performance of this model makes it an excellent SED model for experiments like ZTF, LSST or WFIRST.

A non-local-thermodynamic-equilibrium (NLTE) level population model of the first and second ionisation stages of iron, nickel and cobalt is used to fit a sample of XShooter optical + near-infrared (NIR) spectra of Type Ia supernovae (SNe Ia). From the ratio of the NIR lines to the optical lines limits can be placed on the temperature and density of the emission region. We find a similar evolution of these parameters across our sample. Using the evolution of the Fe II 12$\,$570$\,\mathring{A}\,$to 7$\,$155$\,\mathring{A}\,$line as a prior in fits of spectra covering only the optical wavelengths we show that the 7200$\,\mathring{A}\,$feature is fully explained by [Fe II] and [Ni II] alone. This approach allows us to determine the abundance of Ni II$\,$/$\,$Fe II for a large sample of 130 optical spectra of 58 SNe Ia with uncertainties small enough to distinguish between Chandrasekhar mass (M$_{\text{Ch}}$) and sub-Chandrasekhar mass (sub-M$_{\text{Ch}}$) explosion models. We conclude that the majority (85$\%$) of normal SNe Ia have a Ni/Fe abundance that is in agreement with predictions of sub-M$_{\text{Ch}}$ explosion simulations of $\sim Z_\odot$ progenitors. Only a small fraction (11$\%$) of objects in the sample have a Ni/Fe abundance in agreement with M$_{\text{Ch}}$ explosion models.

As a candidate 'super-Chandrasekhar' or 09dc-like Type Ia supernova (SN Ia), SN 2012dn shares many characteristics with other members of this remarkable class of objects but lacks their extraordinary luminosity. Here, we present and discuss the most comprehensive optical data set of this SN to date, comprised of a densely sampled series of early-time spectra obtained within the Nearby Supernova Factory project, plus photometry and spectroscopy obtained at the VLT about 1 yr after the explosion. The light curves, colour curves, spectral time series and ejecta velocities of SN 2012dn are compared with those of other 09dc-like and normal SNe Ia, the overall variety within the class of 09dc-like SNe Ia is discussed, and new criteria for 09dc-likeness are proposed. Particular attention is directed to additional insight that the late-phase data provide. The nebular spectra show forbidden lines of oxygen and calcium, elements that are usually not seen in late-time spectra of SNe Ia, while the ionisation state of the emitting iron plasma is low, pointing to low ejecta temperatures and high densities. The optical light curves are characterised by an enhanced fading starting ~60 d after maximum and very low luminosities in the nebular phase, which is most readily explained by unusually early formation of clumpy dust in the ejecta. Taken together, these effects suggest a strongly perturbed ejecta density profile, which might lend support to the idea that 09dc-like characteristics arise from a brief episode of interaction with a hydrogen-deficient envelope during the first hours or days after the explosion.

There is now solid experimental evidence of at least one supernova explosion within 100 pc of Earth within the last few million years, from measurements of the short-lived isotope 60Fe in widespread deep-ocean samples, as well as in the lunar regolith and cosmic rays. This is the first established example of a specific dated astrophysical event outside the Solar System having a measurable impact on the Earth, offering new probes of stellar evolution, nuclear astrophysics, the astrophysics of the solar neighborhood, cosmic-ray sources and acceleration, multi-messenger astronomy, and astrobiology. Interdisciplinary connections reach broadly to include heliophysics, geology, and evolutionary biology. Objectives for the future include pinning down the nature and location of the established near-Earth supernova explosions, seeking evidence for others, and searching for other short-lived isotopes such as 26Al and 244Pu. The unique information provided by geological and lunar detections of radioactive 60Fe to assess nearby supernova explosions make now a compelling time for the astronomy community to advocate for supporting multi-disciplinary, cross-cutting research programs.

We present substantial extensions to the Monte Carlo radiative transfer code TARDIS to perform spectral synthesis for type II supernovae. By incorporating a non-LTE ionization and excitation treatment for hydrogen, a full account of free-free and bound-free processes, a self-consistent determination of the thermal state and by improving the handling of relativistic effects, the improved code version includes the necessary physics to perform spectral synthesis for type II supernovae to high precision as required for the reliable inference of supernova properties. We demonstrate the capabilities of the extended version of TARDIS by calculating synthetic spectra for the prototypical type II supernova SN1999em and by deriving a new and independent set of dilution factors for the expanding photosphere method. We have investigated in detail the dependence of the dilution factors on photospheric properties and, for the first time, on changes in metallicity. We also compare our results with two previously published sets of dilution factors by Eastman et al. (1996) and by Dessart & Hillier (2005), and discuss the potential sources of the discrepancies between studies.

Type Ia supernova cosmology depends on the ability to fit and standardize observations of supernova magnitudes with an empirical model. We present here a series of new models of Type Ia Supernova spectral time series that capture a greater amount of supernova diversity than possible with the models that are currently customary. These are entitled SuperNova Empirical MOdels (\textsc{SNEMO}\footnote{https://snfactory.lbl.gov/snemo}). The models are constructed using spectrophotometric time series from $172$ individual supernovae from the Nearby Supernova Factory, comprising more than $2000$ spectra. Using the available observations, Gaussian Processes are used to predict a full spectral time series for each supernova. A matrix is constructed from the spectral time series of all the supernovae, and Expectation Maximization Factor Analysis is used to calculate the principal components of the data. K-fold cross-validation then determines the selection of model parameters and accounts for color variation in the data. Based on this process, the final models are trained on supernovae that have been dereddened using the Fitzpatrick and Massa extinction relation. Three final models are presented here: \textsc{SNEMO2}, a two-component model for comparison with current Type~Ia models; \textsc{SNEMO7}, a seven component model chosen for standardizing supernova magnitudes which results in a total dispersion of $0.100$~mag for a validation set of supernovae, of which $0.087$~mag is unexplained (a total dispersion of $0.113$~mag with unexplained dispersion of $0.097$~mag is found for the total set of training and validation supernovae); and \textsc{SNEMO15}, a comprehensive $15$ component model that maximizes the amount of spectral time series behavior captured.

The observed sub-class of "superluminous" Type Ia supernovae lacks a convincing theoretical explanation. If the emission of such objects were powered exclusively by radioactive decay of 56Ni formed in the explosion, a progenitor mass close to or even above the Chandrasekhar limit for a non-rotating white dwarf star would be required. Masses significantly exceeding this limit can be supported by differential rotation. We, therefore, explore explosions and predict observables for various scenarios resulting from differentially rotating carbon-oxygen white dwarfs close to their respective limit of stability. Specifically, we have investigated a prompt detonation model, detonations following an initial deflagration phase ("delayed detonation" models), and a pure deflagration model. In postprocessing steps, we performed nucleosynthesis and three-dimensional radiative transfer calculations, that allow us, for the first time, to consistently derive synthetic observables from our models. We find that all explosion scenarios involving detonations produce very bright events. The observables predicted for them, however, are inconsistent with any known subclass of Type Ia supernovae. Pure deflagrations resemble 2002cx-like supernovae and may contribute to this class. We discuss implications of our findings for the explosion mechanism and for the existence of differentially rotating white dwarfs as supernova progenitors.

As part of an on-going effort to identify, understand and correct for astrophysics biases in the standardization of Type Ia supernovae (SNIa) for cosmology, we have statistically classified a large sample of nearby SNeIa into those located in predominantly younger or older environments. This classification is based on the specific star formation rate measured within a projected distance of 1kpc from each SN location (LsSFR). This is an important refinement compared to using the local star formation rate directly as it provides a normalization for relative numbers of available SN progenitors and is more robust against extinction by dust. We find that the SNeIa in predominantly younger environments are DY=0.163\pm0.029 mag (5.7 sigma) fainter than those in predominantly older environments after conventional light-curve standardization. This is the strongest standardized SN Ia brightness systematic connected to host-galaxy environment measured to date. The well-established step in standardized brightnesses between SNeIa in hosts with lower or higher total stellar masses is smaller at DM=0.119\pm0.032 mag (4.5 sigma), for the same set of SNeIa. When fit simultaneously, the environment age offset remains very significant, with DY=0.129\pm0.032 mag (4.0 sigma), while the global stellar mass step is reduced to DM=0.064\pm0.029 mag (2.2 sigma). Thus, approximately 70% of the variance from the stellar mass step is due to an underlying dependence on environment-based progenitor age. Standardization using only the SNeIa in younger environments reduces the total dispersion from 0.142\pm0.008 mag to 0.120\pm0.010 mag. We show that as environment ages evolve with redshift a strong bias on measurement of the dark energy equation of state parameters can develop. Fortunately, data to measure and correct for this effect is likely to be available for many next-generation experiments. [abstract shorten]

Type Ia Supernovae (SNe Ia) are widely used to measure the expansion of the Universe. To perform such measurements the luminosity and cosmological redshift ($z$) of the SNe Ia have to be determined. The uncertainty on $z$ includes an unknown peculiar velocity, which can be very large for SNe Ia in the virialized cores of massive clusters. We determine which SNe Ia exploded in galaxy clusters. We then study how the correction for peculiar velocities of host galaxies inside the clusters improves the Hubble residuals. Using 145 SNe Ia from the Nearby Supernova Factory we found 11 candidates for membership in clusters. To estimate the redshift of a cluster we applied the bi-weight technique. Then, we use the galaxy cluster redshift instead of the host galaxy redshift to construct the Hubble diagram. For SNe Ia inside galaxy clusters the dispersion around the Hubble diagram when peculiar velocities are taken into account is smaller in comparison with a case without peculiar velocity correction, with a $wRMS=0.130\pm0.038$ mag instead of $wRMS=0.137\pm0.036$ mag. The significance of this improvement is 3.58 $\sigma$. If we remove the very nearby Virgo cluster member SN2006X ($z<0.01$) from the analysis, the significance decreases to 1.34 $\sigma$. The peculiar velocity correction is found to be highest for the SNe Ia hosted by blue spiral galaxies, with high local specific star formation rate and smaller stellar mass, seemingly counter to what might be expected given the heavy concentration of old, massive elliptical galaxies in clusters. As expected, the Hubble residuals of SNe Ia associated with massive galaxy clusters improve when the cluster redshift is taken as the cosmological redshift of the SN. This fact has to be taken into account in future cosmological analyses in order to achieve higher accuracy for cosmological redshift measurements. Here we provide an approach to do so.

Context. Observations of Type Ia supernovae (SNe Ia) can be used to derive accurate cosmological distances through empirical standardization techniques. Despite this success neither the progenitors of SNe Ia nor the explosion process are fully understood. The U-band region has been less well observed for nearby SNe, due to technical challenges, but is the most readily accessible band for high-redshift SNe. Aims. Using spectrophotometry from the Nearby Supernova Factory, we study the origin and extent of U-band spectroscopic variations in SNe Ia and explore consequences for their standardization and the potential for providing new insights into the explosion process. Methods. We divide the U-band spectrum into four wavelength regions {\lambda}(uNi), {\lambda}(uTi), {\lambda}(uSi) and {\lambda}(uCa). Two of these span the Ca H&K {\lambda}{\lambda} 3934, 3969 complex. We employ spectral synthesis using SYNAPPS to associate the two bluer regions with Ni/Co and Ti. Results. (1) The flux of the uTi feature is an extremely sensitive temperature/luminosity indicator, standardizing the SN peak luminosity to 0.116 $\pm$ 0.011 mag RMS. A traditional SALT2.4 fit on the same sample yields a 0.135 mag RMS. Standardization using uTi also reduces the difference in corrected magnitude between SNe originating from different host galaxy environments. (2) Early U-band spectra can be used to probe the Ni+Co distribution in the ejecta, thus offering a rare window into the source of lightcurve power. (3) The uCa flux further improves standardization, yielding a 0.086 $\pm$ 0.010 mag RMS without the need to include an additional intrinsic dispersion to reach {\chi}$^2$/dof $\sim$ 1. This reduction in RMS is partially driven by an improved standardization of Shallow Silicon and 91T-like SNe.

We present photometry and spectroscopy of nine Type II-P/L supernovae (SNe) with redshifts in the 0.045 < z < 0.335 range, with a view to re-examining their utility as distance indicators. Specifically, we apply the expanding photosphere method (EPM) and the standardized candle method (SCM) to each target, and find that both methods yield distances that are in reasonable agreement with each other. The current record-holder for the highest-redshift spectroscopically confirmed SN II-P is PS1-13bni (z = 0.335 +0.009 -0.012), and illustrates the promise of Type II SNe as cosmological tools. We updated existing EPM and SCM Hubble diagrams by adding our sample to those previously published. Within the context of Type II SN distance measuring techniques, we investigated two related questions. First, we explored the possibility of utilising spectral lines other than the traditionally used Fe II 5169 to infer the photospheric velocity of SN ejecta. Using local well-observed objects, we derive an epoch-dependent relation between the strong Balmer line and Fe II 5169 velocities that is applicable 30 to 40 days post-explosion. Motivated in part by the continuum of key observables such as rise time and decline rates exhibited from II-P to II-L SNe, we assessed the possibility of using Hubble-flow Type II-L SNe as distance indicators. These yield similar distances as the Type II-P SNe. Although these initial results are encouraging, a significantly larger sample of SNe II-L would be required to draw definitive conclusions.

Upcoming high-cadence transient survey programmes will produce a wealth of observational data for Type Ia supernovae. These data sets will contain numerous events detected very early in their evolution, shortly after explosion. Here, we present synthetic light curves, calculated with the radiation hydrodynamical approach Stella for a number of different explosion models, specifically focusing on these first few days after explosion. We show that overall the early light curve evolution is similar for most of the investigated models. Characteristic imprints are induced by radioactive material located close to the surface. However, these are very similar to the signatures expected from ejecta-CSM or ejecta-companion interaction. Apart from the pure deflagration explosion models, none of our synthetic light curves exhibit the commonly assumed power-law rise. We demonstrate that this can lead to substantial errors in the determination of the time of explosion. In summary, we illustrate with our calculations that even with very early data an identification of specific explosion scenarios is challenging, if only photometric observations are available.

The scientific yield of current and future optical surveys is increasingly limited by systematic uncertainties in the flux calibration. This is the case for Type Ia supernova (SN Ia) cosmology programs, where an improved calibration directly translates into improved cosmological constraints. Current methodology rests on models of stars. Here we aim to obtain flux calibration that is traceable to state-of-the-art detector-based calibration. We present the SNIFS Calibration Apparatus (SCALA), a color (relative) flux calibration system developed for the SuperNova Integral Field Spectrograph (SNIFS), operating at the University of Hawaii 2.2 m (UH 88) telescope. By comparing the color trend of the illumination generated by SCALA during two commissioning runs, and to previous laboratory measurements, we show that we can determine the light emitted by SCALA with a long-term repeatability better than 1%. We describe the calibration procedure necessary to control for system aging. We present measurements of the SNIFS throughput as estimated by SCALA observations. The SCALA calibration unit is now fully deployed at the UH\,88 telescope, and with it color-calibration between 4000 {\AA} and 9000 {\AA} is stable at the percent level over a one-year baseline.

Correction of Type Ia Supernova brightnesses for extinction by dust has proven to be a vexing problem. Here we study the dust foreground to the highly reddened SN 2012cu, which is projected onto a dust lane in the galaxy NGC 4772. The analysis is based on multi-epoch, spectrophotometric observations spanning 3,300 - 9,200 {\AA}, obtained by the Nearby Supernova Factory. Phase-matched comparison of the spectroscopically twinned SN 2012cu and SN 2011fe across 10 epochs results in the best-fit color excess of (E(B-V), RMS) = (1.00, 0.03) and total-to-selective extinction ratio of (RV , RMS) = (2.95, 0.08) toward SN 2012cu within its host galaxy. We further identify several diffuse interstellar bands, and compare the 5780 {\AA} band with the dust-to-band ratio for the Milky Way. Overall, we find the foreground dust-extinction properties for SN 2012cu to be consistent with those of the Milky Way. Furthermore we find no evidence for significant time variation in any of these extinction tracers. We also compare the dust extinction curve models of Cardelli et al. (1989), O'Donnell (1994), and Fitzpatrick (1999), and find the predictions of Fitzpatrick (1999) fit SN 2012cu the best. Finally, the distance to NGC4772, the host of SN 2012cu, at a redshift of z = 0.0035, often assigned to the Virgo Southern Extension, is determined to be 16.6$\pm$1.1 Mpc. We compare this result with distance measurements in the literature.

Over the past years type Ia supernovae (SNe Ia) have become a major tool to determine the expansion history of the Universe, and considerable attention has been given to, both, observations and models of these events. However, until now, their progenitors are not known. The observed diversity of light curves and spectra seems to point at different progenitor channels and explosion mechanisms. Here, we present a new way to compare model predictions with observations in a systematic way. Our method is based on the construction of a metric space for SN Ia spectra by means of linear Principal Component Analysis (PCA), taking care of missing and/or noisy data, and making use of Partial Least Square regression (PLS) to find correlations between spectral properties and photometric data. We investigate realizations of the three major classes of explosion models that are presently discussed: delayed-detonation Chandrasekhar-mass explosions, sub-Chandrasekhar-mass detonations, and double-degenerate mergers, and compare them with data. We show that in the PC space all scenarios have observed counterparts, supporting the idea that different progenitors are likely. However, all classes of models face problems in reproducing the observed correlations between spectral properties and light curves and colors. Possible reasons are briefly discussed.

We investigate the consequences of fairly normal Type Ia supernovae being embedded in compact and dense envelopes of carbon and oxygen rich circumstellar material by means of detailed radiation hydrodynamic simulations. Our main focus rests on exploring the effects of the interaction between ejecta and circumstellar material on the ejecta evolution and the broad-band light curve. In our calculations, we find that a strong reverse shock efficiently decelerates and compresses the ejecta material. This leads to a significant broadening of the optical light curve, a longer rise to maximum and a slower decline in the tail phase. During the interaction, substantial radiative energy is generated, which mostly emerges in the extreme ultraviolet and X-ray regime. Only if reprocessing due to radiation-matter interactions is very efficient, a significant boost in the optical light curve is observed. We discuss these findings in particular in the context of the super-luminous event SN 2009dc. As our calculations are able to reproduce a number of its peculiar properties, we conclude that the flavour of the interaction scenario investigated in this work constitutes a promising candidate to explain such 'Super-Chandrasekhar' supernovae.

Based on optical imaging and spectroscopy of the Type II-Plateau SN 2013eq, we present a comparative study of commonly used distance determination methods based on Type II supernovae. The occurrence of SN 2013eq in the Hubble flow (z = 0.041 +/- 0.001) prompted us to investigate the implications of the difference between "angular" and "luminosity" distances within the framework of the expanding photosphere method (EPM) that relies upon a relation between flux and angular size to yield a distance. Following a re-derivation of the basic equations of the EPM for SNe at non-negligible redshifts, we conclude that the EPM results in an angular distance. The observed flux should be converted into the SN rest frame and the angular size, theta, has to be corrected by a factor of (1+z)^2. Alternatively, the EPM angular distance can be converted to a luminosity distance by implementing a modification of the angular size. For SN 2013eq, we find EPM luminosity distances of D_L = 151 +/- 18 Mpc and D_L = 164 +/- 20 Mpc by making use of different sets of dilution factors taken from the literature. Application of the standardized candle method for Type II-P SNe results in an independent luminosity distance estimate (D_L = 168 +/- 16 Mpc) that is consistent with the EPM estimate.

Calculations of synthetic spectropolarimetry are one means to test multi-dimensional explosion models for Type Ia supernovae. In a recent paper, we demonstrated that the violent merger of a 1.1 and 0.9 M$_{\odot}$ 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 two-dimensional double-detonation model and a three-dimensional 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 $\sim$ 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 {\lambda}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 {\lambda}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.

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.

A new method to study the intrinsic color and luminosity of type Ia supernovae (SNe Ia) is presented. A metric space built using principal component analysis (PCA) on spectral series SNe Ia between -12.5 and +17.5 days from B maximum is used as a set of predictors. This metric space is built to be insensitive to reddening. Hence, it does not predict the part of color excess due to dust-extinction. At the same time, the rich variability of SN Ia spectra is a good predictor of a large fraction of the intrinsic color variability. Such metric space is a good predictor of the epoch when the maximum in the B-V color curve is reached. Multivariate Partial Least Square (PLS) regression predicts the intrinsic B band light-curve and the intrinsic B-V color curve up to a month after maximum. This allows to study the relation between the light curves of SNe Ia and their spectra. The total-to-selective extinction ratio RV in the host-galaxy of SNe Ia is found, on average, to be consistent with typical Milky-Way values. This analysis shows the importance of collecting spectra to study SNe Ia, even with large sample publicly available. Future automated surveys as LSST will provide a large number of light curves. The analysis shows that observing accompaning spectra for a significative number of SNe will be important even in the case of "normal" SNe Ia.

Ultraviolet (UV) observations of Type Ia supernovae (SNe Ia) probe the outermost layers of the explosion, and UV spectra of SNe Ia are expected to be extremely sensitive to differences in progenitor composition and the details of the explosion. Here we present the first study of a sample of high signal-to-noise ratio SN Ia spectra that extend blueward of 2900 A. We focus on spectra taken within 5 days of maximum brightness. Our sample of ten SNe Ia spans the majority of the parameter space of SN Ia optical diversity. We find that SNe Ia have significantly more diversity in the UV than in the optical, with the spectral variance continuing to increase with decreasing wavelengths until at least 1800 A (the limit of our data). The majority of the UV variance correlates with optical light-curve shape, while there are no obvious and unique correlations between spectral shape and either ejecta velocity or host-galaxy morphology. Using light-curve shape as the primary variable, we create a UV spectral model for SNe Ia at peak brightness. With the model, we can examine how individual SNe vary relative to expectations based on only their light-curve shape. Doing this, we confirm an excess of flux for SN 2011fe at short wavelengths, consistent with its progenitor having a subsolar metallicity. While most other SNe Ia do not show large deviations from the model, ASASSN-14lp has a deficit of flux at short wavelengths, suggesting that its progenitor was relatively metal rich.

We present results based on observations of SN 2015H which belongs to the small group of objects similar to SN 2002cx, otherwise known as type Iax supernovae. The availability of deep pre-explosion imaging allowed us to place tight constraints on the explosion epoch. Our observational campaign began approximately one day post-explosion, and extended over a period of about 150 days post maximum light, making it one of the best observed objects of this class to date. We find a peak magnitude of M$_r$ = -17.27 $\pm$ 0.07, and a ($\Delta m_{15})_r$ = 0.69 $\pm$ 0.04. Comparing our observations to synthetic spectra generated from simulations of deflagrations of Chandrasekhar mass carbon-oxygen white dwarfs, we find reasonable agreement with models of weak deflagrations that result in the ejection of ~0.2 M$_{\odot}$ of material containing ~0.07 M$_{\odot}$ of 56Ni. The model light curve however, evolves more rapidly than observations, suggesting that a higher ejecta mass is to be favoured. Nevertheless, empirical modelling of the pseudo-bolometric light curve suggests that $\lesssim$0.6 M_sun of material was ejected, implying that the white dwarf is not completely disrupted, and that a bound remnant is a likely outcome.

There has been much debate about the origin of the diffuse $\gamma$--ray background in the MeV range. At lower energies, AGNs and Seyfert galaxies can explain the background, but not above $\simeq$0.3 MeV. Beyond $\sim$10 MeV blazars appear to account for the flux observed. That leaves an unexplained gap for which different candidates have been proposed, including annihilations of WIMPS. One candidate are Type Ia supernovae (SNe Ia). Early studies concluded that they were able to account for the $\gamma$--ray background in the gap, while later work attributed a significantly lower contribution to them. All those estimates were based on SN Ia explosion models which did not reflect the full 3D hydrodynamics of SNe Ia explosions. In addition, new measurements obtained since 2010 have provided new, direct estimates of high-z SNe Ia rates beyond $z\sim$2. We take into account these new advances to see the predicted contribution to the gamma--ray background. We use here a wide variety of explosion models and a plethora of new measurements of SNe Ia rates. SNe Ia still fall short of the observed background. Only for a fit, which would imply $\sim$150\% systematic error in detecting SNe Ia events, do the theoretical predictions approach the observed fluxes. This fit is, however, at odds at the highest redshifts with recent SN Ia rates estimates. Other astrophysical sources such as FSRQs do match the observed flux levels in the MeV regime, while SNe Ia make up to 30--50\% of the observed flux.

HD 188112 is a bright (V = 10.2 mag) hot subdwarf B (sdB) star with a mass too low to ignite core helium burning and is therefore considered as a pre-extremely low mass (ELM) white dwarf (WD). ELM WDs (M $\le$ 0.3 Msun) are He-core objects produced by the evolution of compact binary systems. We present in this paper a detailed abundance analysis of HD 188112 based on high-resolution Hubble Space Telescope (HST) near and far-ultraviolet spectroscopy. We also constrain the mass of the star's companion. We use hybrid non-LTE model atmospheres to fit the observed spectral lines and derive the abundances of more than a dozen elements as well as the rotational broadening of metallic lines. We confirm the previous binary system parameters by combining radial velocities measured in our UV spectra with the already published ones. The system has a period of 0.60658584 days and a WD companion with M $\geq$ 0.70 Msun. By assuming a tidally locked rotation, combined with the projected rotational velocity (v sin i = 7.9 $\pm$ 0.3 km s$^{-1}$) we constrain the companion mass to be between 0.9 and 1.3 Msun. We further discuss the future evolution of the system as a potential progenitor of a (underluminous) type Ia supernova. We measure abundances for Mg, Al, Si, P, S, Ca, Ti, Cr, Mn, Fe, Ni, and Zn, as well as for the trans-iron elements Ga, Sn, and Pb. In addition, we derive upper limits for the C, N, O elements and find HD 188112 to be strongly depleted in carbon. We find evidence of non-LTE effects on the line strength of some ionic species such as Si II and Ni II. The metallic abundances indicate that the star is metal-poor, with an abundance pattern most likely produced by diffusion effects.

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 polarisation 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. (2012) for an explosion triggered during the merger of a 1.1 M$_{\odot}$ and 0.9 M$_{\odot}$ carbon-oxygen white dwarf binary system. Owing to the asymmetries of the ejecta, the polarisation signal varies significantly with viewing angle. We find that polarisation levels for observers in the equatorial plane are modest ($\lesssim$ 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 polarisation and departures from a dominant axis. While the particular model studied here gives a good match to highly-polarised events such as SN 2004dt, it has difficulties in reproducing the low polarisation levels commonly observed in normal Type Ia supernovae. Specifically, we find that significant asymmetries in the element distribution result in a wealth of strong polarisation 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.

The fortuitous occurrence of a type II-Plateau (IIP) supernova, SN 2014bc, in a galaxy for which distance estimates from a number of primary distance indicators are available provides a means with which to cross-calibrate the standardised candle method (SCM) for type IIP SNe. By applying calibrations from the literature we find distance estimates in line with the most precise measurement to NGC 4258 based on the Keplerian motion of masers (7.6$\pm$0.23\,Mpc), albeit with significant scatter. We provide an alternative local SCM calibration by only considering type IIP SNe that have occurred in galaxies for which a Cepheid distance estimate is available. We find a considerable reduction in scatter ($\sigma_I = 0.16$\, mag.), but note that the current sample size is limited. Applying this calibration, we estimate a distance to NGC 4258 of $7.08\pm0.86$ Mpc.

SN 2013dy is a Type Ia supernova for which we have compiled an extraordinary dataset spanning from 0.1 to ~ 500 days after explosion. We present 10 epochs of ultraviolet (UV) through near-infrared (NIR) spectra with HST/STIS, 47 epochs of optical spectra (15 of them having high resolution), and more than 500 photometric observations in the BVrRiIZYJH bands. SN 2013dy has a broad and slowly declining light curve (delta m(B) = 0.92 mag), shallow Si II 6355 absorption, and a low velocity gradient. We detect strong C II in our earliest spectra, probing unburned progenitor material in the outermost layers of the SN ejecta, but this feature fades within a few days. The UV continuum of SN 2013dy, which is strongly affected by the metal abundance of the progenitor star, suggests that SN 2013dy had a relatively high-metallicity progenitor. Examining one of the largest single set of high-resolution spectra for a SN Ia, we find no evidence of variable absorption from circumstellar material. Combining our UV spectra, NIR photometry, and high-cadence optical photometry, we construct a bolometric light curve, showing that SN 2013dy had a maximum luminosity of 10.0^{+4.8}_{-3.8} * 10^{42} erg/s. We compare the synthetic light curves and spectra of several models to SN 2013dy, finding that SN 2013dy is in good agreement with a solar-metallicity W7 model.

The Public European Southern Observatory Spectroscopic Survey of Transient Objects (PESSTO) began as a public spectroscopic survey in April 2012. We describe the data reduction strategy and data products which are publicly available through the ESO archive as the Spectroscopic Survey Data Release 1 (SSDR1). PESSTO uses the New Technology Telescope with EFOSC2 and SOFI to provide optical and NIR spectroscopy and imaging. We target supernovae and optical transients brighter than 20.5mag for classification. Science targets are then selected for follow-up based on the PESSTO science goal of extending knowledge of the extremes of the supernova population. The EFOSC2 spectra cover 3345-9995A (at resolutions of 13-18 Angs) and SOFI spectra cover 0.935-2.53 micron (resolutions 23-33 Angs) along with JHK imaging. This data release contains spectra from the first year (April 2012 - 2013), consisting of all 814 EFOSC2 spectra and 95 SOFI spectra (covering 298 distinct objects), in standard ESO Phase 3 format. We estimate the accuracy of the absolute flux calibrations for EFOSC2 to be typically 15%, and the relative flux calibration accuracy to be about 5%. The PESSTO standard NIR reduction process does not yet produce high accuracy absolute spectrophotometry but the SOFI JHK imaging will improve this. Future data releases will focus on improving the automated flux calibration of the data products.

Stellar evolution models predict the existence of hybrid white dwarfs (WDs) with a carbon-oxygen core surrounded by an oxygen-neon mantle. Being born with masses ~1.1 Msun, hybrid WDs in a binary system may easily approach the Chandrasekhar mass (MCh) by accretion and give rise to a thermonuclear explosion. Here, we investigate an off-centre deflagration in a near-MCh hybrid WD under the assumption that nuclear burning only occurs in carbon-rich material. Performing hydrodynamics simulations of the explosion and detailed nucleosynthesis post-processing calculations, we find that only 0.014 Msun of material is ejected while the remainder of the mass stays bound. The ejecta consist predominantly of iron-group elements, O, C, Si and S. We also calculate synthetic observables for our model and find reasonable agreement with the faint Type Iax SN 2008ha. This shows for the first time that deflagrations in near-MCh WDs can in principle explain the observed diversity of Type Iax supernovae. Leaving behind a near-MCh bound remnant opens the possibility for recurrent explosions or a subsequent accretion-induced collapse in faint Type Iax SNe, if further accretion episodes occur. From binary population synthesis calculations, we find the rate of hybrid WDs approaching MCh to be on the order of 1 percent of the Galactic SN Ia rate.

Recent progress in numerical simulations of thermonuclear supernova explosions brings up a unique opportunity in studying the progenitors of Type Ia supernovae. Coupling state-of-the-art explosion models with detailed hydrodynamical simulations of the supernova remnant evolution and the most up-to-date atomic data for X-ray emission calculations makes it possible to create realistic synthetic X-ray spectra for the supernova remnant phase. Comparing such spectra with high quality observations of supernova remnants could allow to constrain the explosion mechanism and the progenitor of the supernova. The present study focuses in particular on the oxygen emission line properties in young supernova remnants, since different explosion scenarios predict a different amount and distribution of this element. Analysis of the soft X-ray spectra from supernova remnants in the Large Magellanic Cloud and confrontation with remnant models for different explosion scenarios suggests that SNR 0509-67.5 could originate from a delayed detonation explosion and SNR 0519-69.0 from an oxygen-rich merger.

We present observational data for a peculiar supernova discovered by the OGLE-IV survey and followed by the Public ESO Spectroscopic Survey for Transient Objects. The inferred redshift of $z=0.07$ implies an absolute magnitude in the rest-frame $I$-band of M$_{I}\sim-17.6$ mag. This places it in the luminosity range between normal Type Ia SNe and novae. Optical and near infrared spectroscopy reveal mostly Ti and Ca lines, and an unusually red color arising from strong depression of flux at rest wavelengths $<5000$ \AA. To date, this is the only reported SN showing Ti-dominated spectra. The data are broadly consistent with existing models for the pure detonation of a helium shell around a low-mass CO white dwarf and "double-detonation" models that include a secondary detonation of a CO core following a primary detonation in an overlying helium shell.

The measurement of gamma-ray lines from the decay chain of 56Ni provides unique information about the explosion in supernovae. The 56Ni freshly-produced in the supernova powers the optical light curve, as it emits gamma-rays upon its radioactive decay to 56Co and then 56Fe. Gamma-ray lines from 56Co decay are expected to become directly visible through the overlying white dwarf material several weeks after the explosion, as they progressively penetrate the overlying material of the supernova envelope, diluted as it expands. The lines are expected to be Doppler-shifted or broadened from the kinematics of the 56Ni ejecta. With the SPI spectrometer on INTEGRAL and using an improved instrumental background method, we detect the two main lines from 56Co decay at 847 and 1238 keV from SN2014J at 3.3 Mpc, significantly Doppler-broadened, and at intensities (3.65+/-1.21) 10^-4 and (2.27+/-0.69) 10^-4 ph cm^-2s^-1, respectively, at brightness maximum. We measure their rise towards a maximum after about 60-100 days and decline thereafter. The intensity ratio of the two lines is found consistent with expectations from 56Co decay (0.62+/-0.28 at brightness maximum, expected is 0.68). We find that the broad lines seen in the late, gamma-ray transparent phase are not representative for the early gamma-ray emission, and rather notice the emission spectrum to be complex and irregular until the supernova is fully transparent to gamma-rays, with progressive uncovering of the bulk of 56Ni. We infer that the explosion morphology is not spherically symmetric, both in the distribution of 56Ni and of the unburnt material which occults the 56Co emission. Comparing light curves from different plausible models, the resulting 56Ni mass is determined as 0.49+/-0.09 M_sol.

We show that the X-ray line flux of the Mn Kalpha 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 3D explosion models: a near-Chandrasekhar mass delayed detonation and a violent merger of two 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 synthesises 3.5 times more radioactive 55Fe than the merger model. As a result, we find that the peak Mn Kalpha 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. For the delayed-detonation scenario, a line detection is feasible with Chandra up to 3 Mpc for an exposure time of 10^6 s. We find that it should be possible with currently existing X-ray instruments (with exposure times 5x10^5 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).

We develop a new framework for use in exploring Type Ia Supernova (SN Ia) spectra. Combining Principal Component Analysis (PCA) and Partial Least Square analysis (PLS) we are able to establish correlations between the Principal Components (PCs) and spectroscopic/photometric SNe Ia features. The technique was applied to ~120 supernova and ~800 spectra from the Nearby Supernova Factory. The ability of PCA to group together SNe Ia with similar spectral features, already explored in previous studies, is greatly enhanced by two important modifications: (1) the initial data matrix is built using derivatives of spectra over the wavelength, which increases the weight of weak lines and discards extinction, and (2) we extract time evolution information through the use of entire spectral sequences concatenated in each line of the input data matrix. These allow us to define a stable PC parameter space which can be used to characterize synthetic SN Ia spectra by means of real SN features. Using PLS, we demonstrate that the information from important previously known spectral indicators (namely the pseudo-equivalent width (pEW) of Si II 5972 / Si II 6355 and the line velocity of S II 5640 / Si II 6355) at a given epoch, is contained within the PC space and can be determined through a linear combination of the most important PCs. We also show that the PC space encompasses photometric features like B or V magnitudes, B-V color and SALT2 parameters c and x1. The observed colors and magnitudes, that are heavily affected by extinction, cannot be reconstructed using this technique alone. All the above mentioned applications allowed us to construct a metric space for comparing synthetic SN Ia spectra with observations.

The bulk of p isotopes is created in the 'gamma processes' mainly by sequences of photodisintegrations and beta decays in explosive conditions in Type Ia supernovae (SNIa) or in core collapse supernovae (ccSN). The contribution of different stellar sources to the observed distribution of p-nuclei in the Solar System is still under debate. We explore single degenerate Type Ia supernovae in the framework of two-dimensional SNIa delayed-detonation explosion models. Travaglio et al. (2011, TRV11) discussed the sensitivity of p-nuclei production to different SNIa models, i.e. delayed detonations of different strength, deflagrations, and the dependence on selected s-process seed distributions. Here we present a detailed study of p-process nucleosynthesis occuring in SNIa with s-process seeds at different metallicities. Based on the delayed-detonation model DDT-a of TRV11, we analyze the dependence of p-nucleosynthesis on the s-seed distribution obtained from different strengths of the 13C-pocket. Using a Galactic chemical evolution code (see Travaglio et al. 2004) we study the contribution of SNIa to the solar stable p-nuclei. We find that explosions of Chandrasekhar-mass single degenerate systems produce a large amount of p-nuclei in our Galaxy, both in the range of light (A < 120) and heavy p-nuclei, at almost flat average production factors (within a factor of about 3). We discussed in details p-isotopes such as 94Mo with a behavior diverging from the average, which we attribute to uncertainties in the nuclear data or in SNIa modelling.

The nucleosynthesis of proton-rich isotopes is calculated for multi-dimensional Chandrasekhar-mass models of Type Ia supernovae with different metallicities. The predicted abundances of the short-lived radioactive isotopes 92Nb, 97Tc, 98Tc and 146Sm are given in this framework. The abundance seeds are obtained by calculating s-process nucleosynthesis in the material accreted onto a carbon-oxygen white dwarf from a binary companion. A fine grid of s-seeds at different metallicities and 13C-pocket efficiencies is considered. A galactic chemical evolution model is used to predict the contribution of SNIa to the solar system p-nuclei composition measured in meteorites. Nuclear physics uncertainties are critical to determine the role of SNeIa in the production of 92Nb and 146Sm. We find that, if standard Chandrasekhar-mass SNeIa are at least 50% of all SNIa, they are strong candidates for reproducing the radiogenic p-process signature observed in meteorites.

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.

SN 2014J in M82 is the closest detected Type Ia supernova (SN Ia) in at least 28 years and perhaps in 410 years. Despite its small distance of 3.3 Mpc, SN 2014J is surprisingly faint, peaking at V = 10.6 mag, and assuming a typical SN Ia luminosity, we infer an observed visual extinction of A_V = 2.0 +/- 0.1 mag. But this picture, with R_V = 1.6 +/- 0.2, is too simple to account for all observations. We combine 10 epochs (spanning a month) of HST/STIS ultraviolet through near-infrared spectroscopy with HST/WFC3, KAIT, and FanCam photometry from the optical to the infrared and 9 epochs of high-resolution TRES spectroscopy to investigate the sources of extinction and reddening for SN 2014J. We argue that the wide range of observed properties for SN 2014J is caused by a combination of dust reddening, likely originating in the interstellar medium of M82, and scattering off circumstellar material. For this model, roughly half of the extinction is caused by reddening from typical dust (E(B-V ) = 0.45 mag and R_V = 2.6) and roughly half by scattering off LMC-like dust in the circumstellar environment of SN 2014J.

We investigate whether pure deflagration models of Chandrasekhar-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_nuc <~ 1.1 x 10^51 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_nuc ~ 0.5 x 10^51 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.

One of the main questions concerning Type Ia supernovae is the nature of the binary companion of the exploding white dwarf. A major discriminant between different suggested models is the presence and physical properties of circumstellar material at the time of explosion. If present, this material will be ionized by the ultra-violet radiation of the explosion and later recombine. This ionization-recombination should manifest itself as time-variable absorption features that can be detected via multi-epoch high-spectral-resolution observations. Previous studies have shown that the strongest effect is seen in the neutral sodium D lines. We report on observations of neutral sodium absorption features observed in multi-epoch high-resolution spectra of 14 Type Ia supernova events. This is the first multi-epoch high-resolution study to include multiple SNe. No variability in line strength that can be associated with circumstellar material is detected. We find that ~18% of the events in the extended sample exhibit time-variable sodium features associated with circumstellar material. We explore the implication of this study on our understanding of the progenitor systems of Type Ia supernovae via the current Type Ia supernova multi-epoch high-spectral-resolution sample.

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 OI in late-time spectra, indicating unburned oxygen close to the centre of the ejecta. Most explosion models for SNe Ia cannot explain the narrow OI 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 solar masses, respectively, adequately reproduces the early-time observables of SN 2010lp. Moreover, our model predicts oxygen close to the centre of the explosion ejecta, a pre-requisite for narrow OI emission in nebular spectra as observed in SN 2010lp.

We perform a model independent reconstruction of the cosmic expansion rate based on type Ia supernova data. Using the Union 2.1 data set, we show that the Hubble parameter behaviour allowed by the data without making any hypothesis about cosmological model or underlying gravity theory is consistent with a flat LCDM universe having H_0 = 70.43 +- 0.33 and Omega_m=0.297 +- 0.020, weakly dependent on the choice of initial scatter matrix. This is in closer agreement with the recently released Planck results (H_0 = 67.3 +- 1.2, Omega_m = 0.314 +- 0.020) than other standard analyses based on type Ia supernova data. We argue this might be an indication that, in order to tackle subtle deviations from the standard cosmological model present in type Ia supernova data, it is mandatory to go beyond parametrized approaches.

In a companion paper, Seitenzahl et al. (2013) presented a set of three-dimensional delayed detonation models for thermonuclear explosions of near-Chandrasekhar mass white dwarfs (WDs). Here, we present multi-dimensional radiative transfer simulations that provide synthetic light curves and spectra for those models. The model sequence explores both changes in the strength of the deflagration phase (controlled by the ignition configuration) and the WD central density. In agreement with previous studies, we find that the strength of the deflagration significantly affects the explosion and the observables. Variations in the central density also have an influence on both brightness and colour, but overall it is a secondary parameter in our set of models. In many respects, the models yield a good match to normal Type Ia supernovae (SNe Ia): peak brightness, rise/decline time scales and synthetic spectra are all in reasonable agreement. There are, however, several differences. In particular, the models are too red around maximum light, manifest spectral line velocities that are a little too high and yield I-band light curves that do not match observations. Although some of these discrepancies may simply relate to approximations made in the modelling, some pose real challenges to the models. If viewed as a complete sequence, our models do not reproduce the observed light-curve width-luminosity relation (WLR) of SNe Ia: all our models show similar B-band decline rates, irrespective of peak brightness. This suggests that simple variations in the strength of the deflagration phase in Chandrasekhar-mass deflagration-to-detonation models do not readily explain the observed diversity of normal SNe Ia. This may imply that some other parameter within the Chandrasekhar-mass paradigm is key to the WLR, or that a substantial fraction of normal SNe Ia arise from an alternative explosion scenario.

In this letter a late-phase spectrum of SN 2010lp, a subluminous Type Ia supernova (SN Ia), is presented and analysed. As in 1991bg-like SNe Ia at comparable epochs, the spectrum is characterised by relatively broad [FeII] and [CaII] emission lines. However, instead of narrow [FeIII] and [CoIII] lines that dominate the emission from the innermost regions of 1991bg-like SNe, SN 2010lp shows [OI] 6300,6364 emission, usually associated with core-collapse SNe and never observed in a subluminous thermonuclear explosion before. The [OI] feature has a complex profile with two strong, narrow emission peaks. This suggests oxygen to be distributed in a non-spherical region close to the centre of the ejecta, severely challenging most thermonuclear explosion models discussed in the literature. We conclude that given these constraints violent mergers are presently the most promising scenario to explain SN 2010lp.

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.

Besides the fact that the gamma-ray emission due to radioactive decays is responsible for powering the light curves of Type Ia supernovae (SNe Ia), gamma rays themselves are of particular interest as a diagnostic tool because they provide a direct way to obtain deeper insights into the nucleosynthesis and the kinematics of these explosion events. 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 use three-dimensional explosion models and perform 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. In contrast to the optical regime, the gamma-ray emission of our two chosen models proves to be rather different. The almost direct connection of the emission of gamma rays to fundamental physical processes occuring 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 and 16 Mpc for an exposure time of 10^6 s, respectively. The possibility to detect 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.

In the SD scenario of SNe Ia the companion survives the SN explosion and thus should be visible near the center of the SN remnant and may show some unusual features. A promising approach to test progenitor models of SNe Ia is to search for the companion in SNRs. Here we present the results of 3D hydrodynamics simulations of the interaction between the SN Ia blast wave and a MS companion taking into consideration its orbital motion and spin. The primary goal of this work is to investigate the rotation of surviving companions after SN Ia explosions in the WD+MS scenario. We use Eggleton's code including the optically thick accretion wind model to obtain realistic models of companions. The impact of the SN blast wave on these companions is followed in 3D hydrodynamic simulations employing the SPH code GADGET3. We find that the rotation of the companion does not significantly affect the amount of stripped mass and the kick velocity caused by the SN impact. However, in our simulations, the rotational velocity of the companion is significantly reduced to about 14% to 32% of its pre-explosion value due to the expansion of the companion and the fact that 55%-89% of the initial angular momentum is carried away by the stripped matter. Compared with the observed rotational velocity of the presumed companion star of Tycho's SN, Tycho G, of 6 km/s the final rotational velocity we obtain is still higher by at least a factor of two. Whether this difference is significant, and may cast doubts on the suggestion that Tycho G is the companion of SN 1572, has to be investigated in future studies. Based on binary population synthesis results we present, for the first time, the expected distribution of rotational velocities of companions after the explosion which may provide useful information for the identification of the surviving companion in observational searches in other historical SNRs.

Motivated by the fact that calibrated light curves of Type Ia supernovae (SNe Ia) have become a major tool to determine the expansion history of the Universe, considerable attention has been given to, both, observations and models of these events over the past 15 years. Here, we summarize new observational constraints, address recent progress in modeling Type Ia supernovae by means of three-dimensional hydrodynamic simulations, and discuss several of the still open questions. It will be be shown that the new models have considerable predictive power which allows us to study observable properties such as light curves and spectra without adjustable non-physical parameters. This is a necessary requisite to improve our understanding of the explosion mechanism and to settle the question of the applicability of SNe Ia as distance indicators for cosmology. We explore the capabilities of the models by comparing them with observations and we show how such models can be applied to study the origin of the diversity of SNe Ia.

We explore SNIa as p-process sources in the framework of two-dimensional SNIa models using enhanced s-seed distributions as directly obtained from a sequence of thermal pulse instabilities. The SNIa WD precursor is assumed to have reached the Chandrasekhar mass limit in a binary system by mass accretion from a giant/main sequence companion. We apply the tracer-particle method to reconstruct the nucleosynthesis from the thermal histories of Lagrangian particles, passively advected in the hydrodynamic calculations. For each particle we follow the explosive nucleosynthesis with a detailed nuclear reaction network. We select tracers within the typical temperature range for p-process production, 1.5-3.7 109K, and analyse in detail their behaviour, exploring the influence of different s-process distributions on the p-process nucleosynthesis. We find that SNIa contribute to a large fraction of p-nuclei, both the light p-nuclei and the heavy-p nuclei at a quite flat average production factor. For the first time, the very abundant Ru and Mo p-isotopes are reproduced at the same level as the heavy p-nuclei. We investigate the metallicity effect on the p-process production. Starting with a range of s-seeds distributions obtained for different metallicities, running SNIa two-dimensional models and using a simple chemical evolution code, we give estimates of the SNIa contribution to the solar p-process composition. We find that SNIa contribute for at least 50% at the solar p-nuclei composition, in a primary way.

We present 32 epochs of optical (3300-9700 \AA) spectrophotometric observations of the nearby quintessential "normal" type Ia supernova (SN Ia) SN 2011fe in the galaxy M101, extending from -15 to +97 d with respect to B-band maximum, obtained by the Nearby Supernova Factory collaboration. SN 2011fe is the closest (\mu = 29.04) and brightest (Bmax = 9.94 mag) SN Ia observed since the advent of modern large scale programs for the intensive periodic followup of supernovae. Both synthetic light curve measurements and spectral feature analysis attest to the normality of SN 2011fe. There is very little evidence for reddening in its host galaxy. The homogeneous calibration, intensive time sampling, and high signal-to-noise ratio of the data set make it unique. Thus it is ideal for studying the physics of SN Ia explosions in detail, and for furthering the use of SNe Ia as standardizable candles for cosmology. Several such applications are shown, from the creation of a bolometric light curve and measurement of the 56Ni mass, to the simulation of detection thresholds for unburned carbon, direct comparisons with other SNe Ia, and existing spectral templates.

H and He features in photospheric spectra have seldom been used to infer quantitatively the properties of Type IIb, Ib and Ic supernovae (SNe IIb, Ib and Ic) and their progenitor stars. Most radiative transfer models ignored NLTE effects, which are extremely strong especially in the He-dominated zones. In this paper, a comprehensive set of model atmospheres for low-mass SNe IIb/Ib/Ic is presented. Long-standing questions such as how much He can be contained in SNe Ic, where He lines are not seen, can thus be addressed. The state of H and He is computed in full NLTE, including the effect of heating by fast electrons. The models are constructed to represent iso-energetic explosions of the same stellar core with differently massive H/He envelopes on top. The synthetic spectra suggest that 0.06 - 0.14 M_sun of He and even smaller amounts of H suffice for optical lines to be present, unless ejecta asymmetries play a major role. This strongly supports the conjecture that low-mass SNe Ic originate from binaries where progenitor mass loss can be extremely efficient.

Extremely luminous, super-Chandrasekhar (SC) Type Ia Supernovae (SNe Ia) are as yet an unexplained phenomenon. We analyse a well-observed SN of this class, SN 2009dc, by modelling its photospheric spectra with a spectral synthesis code, using the technique of 'Abundance Tomography'. We present spectral models based on different density profiles, corresponding to different explosion scenarios, and discuss their consistency. First, we use a density structure of a simulated explosion of a 2 M_sun rotating C-O white dwarf (WD), which is often proposed as a possibility to explain SC SNe Ia. Then, we test a density profile empirically inferred from the evolution of line velocities (blueshifts). This model may be interpreted as a core-collapse SN with an ejecta mass ~ 3 M_sun. Finally, we calculate spectra assuming an interaction scenario. In such a scenario, SN 2009dc would be a standard WD explosion with a normal intrinsic luminosity, and this luminosity would be augmented by interaction of the ejecta with a H-/He-poor circumstellar medium. We find that no model tested easily explains SN 2009dc. With the 2 M_sun WD model, our abundance analysis predicts small amounts of burning products in the intermediate-/high-velocity ejecta (v > 9000 km/s). However, in the original explosion simulations, where the nuclear energy release per unit mass is large, burned material is present at high v. This contradiction can only be resolved if asymmetries strongly affect the radiative transfer or if C-O WDs with masses significantly above 2 M_sun exist. In a core-collapse scenario, low velocities of Fe-group elements are expected, but the abundance stratification in SN 2009dc seems 'SN Ia-like'. The interaction-based model looks promising, and we have some speculations on possible progenitor configurations. However, radiation-hydro simulations will be needed to judge whether this scenario is realistic at all.

The properties of underluminous type Ia supernovae (SNe Ia) of the 91bg subclass have yet to be theoretically understood. Here, we take a closer look at the structure of the dim SN Ia 2005bl. We infer the abundance and density profiles needed to reproduce the observed spectral evolution between -6 d and +12.9 d with respect to B maximum. Initially, we assume the density structure of the standard explosion model W7; then we test whether better fits to the observed spectra can be obtained using modified density profiles with different total masses and kinetic energies. Compared to normal SNe Ia, we find a lack of burning products especially in the rapidly-expanding outer layers (v>~15000 km/s). The zone between ~8500 and 15000 km/s is dominated by oxygen and includes some amount of intermediate mass elements. At lower velocities, intermediate mass elements dominate. This holds down to the lowest zones investigated in this work. This fact, together with negligible-to-moderate abundances of Fe-group elements, indicates large-scale incomplete Si burning or explosive O burning, possibly in a detonation at low densities. Consistently with the reduced nucleosynthesis, we find hints of a kinetic energy lower than that of a canonical SN Ia: The spectra strongly favour reduced densities at >~13000 km/s compared to W7, and are very well fitted using a rescaled W7 model with original mass (1.38 M_sun), but a kinetic energy reduced by ~30 % (i.e. from 1.33e51 erg to 0.93e51 erg).

2002cx-like supernovae are a sub-class of sub-luminous Type Ia supernovae. Their light curves and spectra are characterized by distinct features that indicate strong mixing of the explosion ejecta. Pure turbulent deflagrations have been shown to produce such mixed ejecta. Here, we present hydrodynamics, nucleosynthesis and radiative transfer calculations for a 3D full-star deflagration of a Chandrasekhar-mass white dwarf. Our model is able to reproduce the characteristic observational features of SN 2005hk (a proto-typical 2002cx-like supernova), not only in the optical, but also in the near-infrared. For that purpose we present, for the first time, five near-infrared spectra of SN 2005hk from -0.2 to 26.6 days with respect to B-band maximum. Since our model burns only small parts of the initial white dwarf, it fails to completely unbind the white dwarf and leaves behind a bound remnant of ~1.03 solar masses -- consisting mainly of unburned carbon and oxygen, but also enriched by some amount of intermediate-mass and iron-group elements from the explosion products that fall back on the remnant. We discuss possibilities for detecting this bound remnant and how it might influence the late-time observables of 2002cx-like SNe.

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. 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 SNe Ia. We argue that the brightness of an explosion realized via the violent merger model is determined by the mass of nickel 56 produced in the detonation of the primary WD. We use a set of sub-Chandrasekhar mass WD detonation models to derive a relationship between primary WD mass and expected peak bolometric brightness. We use this relationship to convert the masses of merging primary WDs from binary population models to a predicted distribution of explosion brightness. We find a striking similarity between the shape of our theoretical peak-magnitude distribution and that observed for SNe Ia: our model produces a brightness 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.35 Msun from a slightly-evolved helium star companion. We also find that with high probability, violent WD mergers involving the most massive primaries (which should produce bright SNe) have delay times <500 Myr (Abridged).

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 perform multi-dimensional, time-dependent radiation transfer simulations for hard X-ray and gamma-ray emissions, following radioactive decays of 56Ni and 56Co, for two-dimensional delayed detonation models of Type Ia supernovae (SNe Ia). The synthetic spectra and light curves are compared with the sensitivities of current and future observatories for an exposure time of 10^6 seconds. The non-detection of the gamma-ray signal from SN 2011fe at 6.4 Mpc by SPI on board INTEGRAL places an upper limit for the mass of 56Ni of \lesssim 1.0 Msun, independently from observations in any other wavelengths. Signals from the newly formed radioactive species have not been convincingly measured yet from any SN Ia, but the future X-ray and gamma-ray missions are expected to deepen the observable horizon to provide the high energy emission data for a significant SN Ia sample. We predict that the hard X-ray detectors on board NuStar (launched in 2012) or ASTRO-H (scheduled for launch in 2014) will reach to SNe Ia at \sim15 Mpc, i.e., one SN every few years. Furthermore, according to the present results, the soft gamma-ray detector on board ASTRO-H will be able to detect the 158 keV line emission up to \sim25 Mpc, i.e., a few SNe Ia per year. Proposed next generation gamma-ray missions, e.g., GRIPS, could reach to SNe Ia at \sim20 - 35 Mpc by MeV observations. Those would provide new diagnostics and strong constraints on explosion models, detecting rather directly the main energy source of supernova light.

The identity of the progenitor systems of SNe Ia is still uncertain. In the single-degenerate (SD) scenario, the interaction between the SN blast wave and the outer layers of a main sequence (MS) companion star strips off H-rich material which is then mixed into the ejecta. Strong contamination of the SN ejecta with stripped material could lead to a conflict with observations of SNe Ia. This constrains the SD progenitor model. In this work, our previous simulations based on simplified progenitor donor stars have been updated by adopting more realistic progenitor-system models that result from fully detailed, state-of-the-art binary evolution calculations. We use Eggleton's stellar evolution code including the optically thick accretion wind model and the possibility of the effects of accretion disk instabilities to obtain realistic models of companions for different progenitor systems. The impact of the SN blast wave on these companion stars is followed in three-dimensional hydrodynamic simulations employing the SPH code GADGET3. We find that the stripped masses range from 0.11 to 0.18 M_sun. The kick velocity is between 51 and 105 km/s. We find that the stripped mass and kick velocity depend on the ratio of the orbital separation to the radius of a companion. They can be fitted by a power law for a given companion model. However, the structure of the companion star is also important for the amount of stripped material. With more realistic companion star models than in previous studies, our simulations show that the H masses stripped from companions are inconsistent with the best observational limits (< 0.01 M_sun) derived from nebular spectra. However, a rigorous forward modeling based on impact simulations with radiation transfer is required to reliably predict observable signatures of the stripped H and to conclusively assess the viability of the considered SN Ia progenitor scenario.

We present nine near-infrared (NIR) spectra of supernova (SN) 2005cf at epochs from -10 d to +42 d with respect to B-band maximum, complementing the existing excellent data sets available for this prototypical Type Ia SN at other wavelengths. The spectra show a time evolution and spectral features characteristic of normal Type Ia SNe, as illustrated by a comparison with SNe 1999ee, 2002bo and 2003du. The broad-band spectral energy distribution (SED) of SN 2005cf is studied in combined ultraviolet (UV), optical and NIR spectra at five epochs between ~ 8 d before and ~ 10 d after maximum light. We also present synthetic spectra of the hydrodynamic explosion model W7, which reproduce the key properties of SN 2005cf not only at UV-optical as previously reported, but also at NIR wavelengths. From the radiative-transfer calculations we infer that fluorescence is the driving mechanism that shapes the SED of SNe Ia. In particular, the NIR part of the spectrum is almost devoid of absorption features, and instead dominated by fluorescent emission of both iron-group material and intermediate-mass elements at pre-maximum epochs, and pure iron-group material after maximum light. A single P-Cygni feature of Mg II at early epochs and a series of relatively unblended Co II lines at late phases allow us to constrain the regions of the ejecta in which the respective elements are abundant.

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.

Mergers of two carbon-oxygen white dwarfs have long been suspected to be progenitors of Type Ia Supernovae. Here we present our modifications to the cosmological smoothed particle hydrodynamics code Gadget to apply it to stellar physics including but not limited to mergers of white dwarfs. We demonstrate a new method to map a one-dimensional profile of an object in hydrostatic equilibrium to a stable particle distribution. We use the code to study the effect of initial conditions and resolution on the properties of the merger of two white dwarfs. We compare mergers with approximate and exact binary initial conditions and find that exact binary initial conditions lead to a much more stable binary system but there is no difference in the properties of the actual merger. In contrast, we find that resolution is a critical issue for simulations of white dwarf mergers. Carbon burning hotspots which may lead to a detonation in the so-called violent merger scenario emerge only in simulations with sufficient resolution but independent of the type of binary initial conditions. We conclude that simulations of white dwarf mergers which attempt to investigate their potential for Type Ia supernovae should be carried out with at least 10^6 particles.

The nearby supernova SN 2011fe can be observed in unprecedented detail. Therefore, it is an important test case for Type Ia supernova (SN Ia) models, which may bring us closer to understanding the physical nature of these objects. Here, we explore how available and expected future observations of SN 2011fe can be used to constrain SN Ia explosion scenarios. We base our discussion on three-dimensional simulations of a delayed detonation in a Chandrasekhar-mass white dwarf and of a violent merger of two white dwarfs-realizations of explosion models appropriate for two of the most widely-discussed progenitor channels that may give rise to SNe Ia. Although both models have their shortcomings in reproducing details of the early and near-maximum spectra of SN 2011fe obtained by the Nearby Supernova Factory (SNfactory), the overall match with the observations is reasonable. The level of agreement is slightly better for the merger, in particular around maximum, but a clear preference for one model over the other is still not justified. Observations at late epochs, however, hold promise for discriminating the explosion scenarios in a straightforward way, as a nucleosynthesis effect leads to differences in the 55Co production. SN 2011fe is close enough to be followed sufficiently long to study this effect.

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.

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.

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 \mathrm{M_\odot}$ and $1.1 \mathrm{M_\odot}$ 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 lightcurves of our merger, which produces about $0.62 \mathrm{M_\odot}$ of $^{56}\mathrm{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.

Thermonuclear explosions may arise in binaries in which a CO white dwarf (WD) accretes He from a companion. If the accretion rate allows a sufficiently large mass of He to accumulate prior to ignition of nuclear burning, the He surface layer may detonate, giving rise to an astrophysical transient. Detonation of the He layer generates shock waves that propagate into the underlying CO WD. This might directly ignite a detonation at the edge of the CO WD or compress the core of the WD sufficiently to trigger a CO detonation near the centre. If either ignition mechanism works, the two detonations can release sufficient energy to completely unbind the WD. Here we extend our 2D studies of this double-detonation model to low-mass CO WDs. 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 realized in nature. We compute spectra and light curves for models in which either an edge-lit or compression-triggered CO detonation is assumed to occur and compare these to 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 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, NIR observations and detailed spectroscopic analysis would be needed to determine whether core detonation occurred. We comment on the implications of our results for understanding peculiar astrophysical transients including SN 2002bj, SN 2010X and SN 2005E.

Based on the largest homogeneously reduced set of Type Ia supernova luminosity data currently available -- the Union2 sample -- we reconstruct the expansion history of the Universe in a model-independent approach. Our method tests the geometry of the Universe directly without reverting to any assumptions made on its energy content. This allows us to constrain Dark Energy models and non-standard cosmologies in a straightforward way. The applicability of the presented method is not restricted to testing cosmological models. It can be a valuable tool for pointing out systematic errors hidden in the supernova data and planning future Type Ia supernova cosmology campaigns.

We present theoretical delay times and rates of thermonuclear explosions that are thought to produce Type Ia supernovae, including the double-detonation sub-Chandrasekhar mass model, using the population synthesis binary evolution code StarTrack. If detonations of sub-Chandrasekhar mass carbon-oxygen white dwarfs following a detonation in an accumulated layer of helium on the white dwarf's surface ("double-detonation" models) are able to produce thermonuclear explosions which are characteristically similar to those of SNe Ia, then these sub-Chandrasekhar mass explosions may account for at least some substantial fraction of the observed SN Ia rate. Regardless of whether all double-detonations look like 'normal' SNe Ia, in any case the explosions are expected to be bright and thus potentially detectable. Additionally, we find that the delay time distribution of double-detonation sub-Chandrasekhar mass SNe Ia can be divided into two distinct formation channels: the 'prompt' helium-star channel with delay times <500 Myr (~10% of all sub-Chandras), and the 'delayed' double white dwarf channel, with delay times >800 Myr spanning up to a Hubble time (~90%). These findings coincide with recent observationally-derived delay time distributions which have revealed that a large number of SNe Ia are prompt with delay times <500 Myr, while a significant fraction also have delay times spanning ~1 Gyr to a Hubble time.

We explore SNe Ia as p-process sites in the framework of two-dimensional SN Ia delayed detonation and pure deflagration models. The WD precursor is assumed to have reached the Chandrasekhar mass in a binary system by mass accretion from a giant/main sequence companion. We use enhanced s-seed distributions, obtained from a sequence of thermal pulse instabilities both in the AGB phase and in the accreted material. We apply the tracer-particle method to reconstruct the nucleosynthesis by the thermal histories of Lagrangian particles, passively advected in the hydrodynamic calculations. For each particle we follow the explosive nucleosynthesis with a detailed network for all isotopes up to 209Bi. We find that SNe Ia can produce a large amount of p-nuclei, both the light p-nuclei below A=120 and the heavy-p nuclei, at quite flat average production factors, tightly related to the s-process seed distribution. For the first time, we find a stellar source able to produce both, light and heavy p-nuclei almost at the same level as 56Fe, including the very debated neutron magic 92,94Mo and 96,98Ru. We also find that there is an important contribution from p-process nucleosynthesis to the s-only nuclei 80Kr, 86Sr, to the neutron magic 90Zr, and to the neutron-rich 96Zr. Finally, we investigate the metallicity effect on p-process. Starting with different s-process seed distributions, for two metallicities Z = 0.02 and Z = 0.001, running SNe Ia models with different initial composition, we estimate that SNe Ia can contribute to, at least, 50% of the solar p-process composition.

The interpretation of supernova (SN) spectra is essential for deriving SN ejecta properties such as density and composition, which in turn can tell us about their progenitors and the explosion mechanism. A very large number of atomic processes are important for spectrum formation. Several tools for calculating SN spectra exist, but they mainly focus on the very early or late epochs. The intermediate phase, which requires a NLTE treatment of radiation transport has rarely been studied. In this paper we present a new SN radiation transport code, NERO, which can look at those epochs. All the atomic processes are treated in full NLTE, under a steady-state assumption. This is a valid approach between roughly 50 and 500 days after the explosion depending on SN type. This covers the post-maximum photospheric and the early and the intermediate nebular phase. As a test, we compare NERO to the radiation transport code of Jerkstrand et al. (2011) and to the nebular code of Mazzali et al. (2001). All three codes have been developed independently and a comparison provides a valuable opportunity to investigate their reliability. Currently, NERO is one-dimensional and can be used for predicting spectra of synthetic explosion models or for deriving SN properties by spectral modelling. To demonstrate this, we study the spectra of the 'normal' SN Ia 2005cf between 50 and 350 days after the explosion and identify most of the common SN Ia line features at post maximum epochs.

We propose to perform a continuously scanning all-sky survey from 200 keV to 80 MeV achieving a sensitivity which is better by a factor of 40 or more compared to the previous missions in this energy range. The Gamma-Ray Imaging, Polarimetry and Spectroscopy (GRIPS) mission addresses fundamental questions in ESA's Cosmic Vision plan. Among the major themes of the strategic plan, GRIPS has its focus on the evolving, violent Universe, exploring a unique energy window. We propose to investigate $\gamma$-ray bursts and blazars, the mechanisms behind supernova explosions, nucleosynthesis and spallation, the enigmatic origin of positrons in our Galaxy, and the nature of radiation processes and particle acceleration in extreme cosmic sources including pulsars and magnetars. The natural energy scale for these non-thermal processes is of the order of MeV. Although they can be partially and indirectly studied using other methods, only the proposed GRIPS measurements will provide direct access to their primary photons. GRIPS will be a driver for the study of transient sources in the era of neutrino and gravitational wave observatories such as IceCUBE and LISA, establishing a new type of diagnostics in relativistic and nuclear astrophysics. This will support extrapolations to investigate star formation, galaxy evolution, and black hole formation at high redshifts.

The origin of subluminous Type Ia supernovae (SNe Ia) has long eluded any explanation, as all Chandrasekhar-mass models have severe problems reproducing them. Recently, it has been proposed that violent mergers of two white dwarfs of 0.9 M_sun could lead to subluminous SNe Ia events that resemble 1991bg-like SNe~Ia. Here we investigate whether this scenario still works for mergers of two white dwarfs with a mass ratio smaller than one. We aim to determine the range of mass ratios for which a detonation still forms during the merger, as only those events will lead to a SN Ia. This range is an important ingredient for population synthesis and one decisive point to judge the viability of the scenario. In addition, we perform a resolution study of one of the models. Finally we discuss the connection between violent white dwarf mergers with a primary mass of 0.9 M_sun and 1991bg-like SNe Ia. The latest version of the smoothed particle hydrodynamics code Gadget3 is used to evolve binary systems with different mass ratios until they merge. We analyze the result and look for hot spots in which detonations can form. We show that mergers of two white dwarfs with a primary white dwarf mass of ~0.9 M_sun and a mass ratio larger than about $0.8$ robustly reach the conditions we require to ignite a detonation and thus produce thermonuclear explosions during the merger itself. We also find that while our simulations do not yet completely resolve the hot spots, increasing the resolution leads to conditions that are even more likely to ignite detonations. (abridged)

SN 2009dc shares similarities with normal Type Ia supernovae, but is clearly overluminous, with a (pseudo-bolometric) peak luminosity of log(L) = 43.47 [erg/s]. Its light curves decline slowly over half a year after maximum light, and the early-time near-IR light curves show secondary maxima, although the minima between the first and second peaks are not very pronounced. Bluer bands exhibit an enhanced fading after ~200 d, which might be caused by dust formation or an unexpectedly early IR catastrophe. The spectra of SN 2009dc are dominated by intermediate-mass elements and unburned material at early times, and by iron-group elements at late phases. Strong C II lines are present until ~2 weeks past maximum, which is unprecedented in thermonuclear SNe. The ejecta velocities are significantly lower than in normal and even subluminous SNe Ia. No signatures of CSM interaction are found in the spectra. Assuming that the light curves are powered by radioactive decay, analytic modelling suggests that SN 2009dc produced ~1.8 solar masses of 56Ni assuming the smallest possible rise time of 22 d. Together with a derived total ejecta mass of ~2.8 solar masses, this confirms that SN 2009dc is a member of the class of possible super-Chandrasekhar-mass SNe Ia similar to SNe 2003fg, 2006gz and 2007if. A study of the hosts of SN 2009dc and other superluminous SNe Ia reveals a tendency of these SNe to explode in low-mass galaxies. A low metallicity of the progenitor may therefore be an important pre-requisite for producing superluminous SNe Ia. We discuss a number of explosion scenarios, ranging from super-Chandrasekhar-mass white-dwarf progenitors over dynamical white-dwarf mergers and Type I 1/2 SNe to a core-collapse origin of the explosion. None of the models seem capable of explaining all properties of SN 2009dc, so that the true nature of this SN and its peers remains nebulous.

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_solar with corresponding shell masses from 0.126 down to 0.0035 M_solar. 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.

Using the population synthesis binary evolution code StarTrack, we present theoretical rates and delay times of Type Ia supernovae arising from various formation channels. These channels include binaries in which the exploding white dwarf reaches the Chandrasekhar mass limit (DDS, SDS, and helium-rich donor scenario) as well as the sub-Chandrasekhar mass scenario, in which a white dwarf accretes from a helium-rich companion and explodes as a SN Ia before reaching the Chandrasekhar mass limit. We find that using a common envelope parameterization employing energy balance with alpha=1 and lambda=1, the supernova rates per unit mass (born in stars) of sub-Chandrasekhar mass SNe Ia exceed those of all other progenitor channels at epochs t=0.7 - 4 Gyr for a burst of star formation at t=0. Additionally, the delay time distribution of the sub-Chandrasekhar model can be divided in to two distinct evolutionary channels: the `prompt' helium-star channel with delay times < 500 Myr, and the `delayed' double white dwarf channel with delay times > 800 Myr spanning up to a Hubble time. These findings are in agreement with recent observationally-derived delay time distributions which predict that a large number of SNe Ia have delay times < 1 Gyr, with a significant fraction having delay times < 500 Myr. We find that the DDS channel is also able to account for the observed rates of SNe Ia. However, detailed simulations of white dwarf mergers have shown that most of these mergers will not lead to SNe Ia but rather to the formation of a neutron star via accretion-induced collapse. If this is true, our standard population synthesis model predicts that the only progenitor channel which can account for the rates of SNe Ia is the sub-Chandrasekhar mass scenario, and none of the other progenitors considered can fully account for the observed rates.

Abridged. In the double detonation scenario for Type Ia supernovae (SNe Ia) a detonation initiates in a shell of He-rich material accreted from a companion star by a sub-Chandrasekhar-mass White Dwarf (WD). This shell detonation drives a shock front into the carbon-oxygen (C/O) WD that triggers a secondary detonation in the core. The core detonation results in a complete disruption of the WD. Earlier studies concluded that this scenario has difficulties in accounting for the observed properties of SNe Ia since the explosion ejecta are surrounded by the products of explosive He burning in the shell. Recently, it was proposed that detonations might be possible for much less massive He shells than previously assumed. Moreover, it was shown that even detonations of these minimum He shell masses robustly trigger detonations of the C/O core. Here we present time-dependent multi-wavelength radiative transfer calculations for models with minimum He shell mass and derive synthetic observables for both the optical and {\gamma}-ray spectral regions. These differ strongly from those found in earlier simulations of sub-Chandrasekhar-mass explosions in which more massive He shells were considered. Our models predict light curves which cover both the range of brightnesses and the rise and decline times of observed SNe Ia. However, their colours and spectra do not match the observations. In particular, their B-V colours are generally too red. We show that this discrepancy is mainly due to the composition of the burning products of the He shell of our models which contain significant amounts of Ti and Cr. Using a toy model, we also show that the burning products of the He shell depend crucially on its initial composition. This leads us to conclude that good agreement between sub-Chandrasekhar-mass explosions and observed SNe Ia may still be feasible but further study of the shell properties is required.

Explosions of sub-Chandrasekhar-mass white dwarfs are one alternative to the standard Chandrasekhar-mass model of Type Ia supernovae. They are interesting since binary systems with sub-Chandrasekhar-mass primary white dwarfs should be common and this scenario would suggest a simple physical parameter which determines the explosion brightness, namely the mass of the exploding white dwarf. Here we perform one-dimensional hydrodynamical simulations, associated post-processing nucleosynthesis and multi-wavelength radiation transport calculations for pure detonations of carbon-oxygen white dwarfs. The light curves and spectra we obtain from these simulations are in good agreement with observed properties of Type Ia supernovae. In particular, for white dwarf masses from 0.97 - 1.15 Msun we obtain 56Ni masses between 0.3 and 0.8 Msun, sufficient to capture almost the complete range of Type Ia supernova brightnesses. Our optical light curve rise times, peak colours 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 Angstroms. We therefore suggest that sub-Chandrasekhar mass explosions are a viable model for Type Ia supernovae 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 white dwarf.

Type Ia supernovae (SNe Ia) are a prime tool in observational cosmology. A relation between their peak luminosities and the shapes of their light curves allows to infer their intrinsic luminosities and to use them as distance indicators. This relation has been established empirically. However, a theoretical understanding is necessary in order to get a handle on the systematics in SN Ia cosmology. Here, a model reproducing the observed diversity of normal SNe Ia is presented. The challenge in the numerical implementation arises from the vast range of scales involved in the physical mechanism. Simulating the supernova on scales of the exploding white dwarf requires specific models of the microphysics involved in the thermonuclear combustion process. Such techniques are discussed and results of simulations are presented.

The nucleosynthetic characteristics of various explosion mechanisms of Type Ia supernovae (SNe Ia) is explored based on three two-dimensional explosion simulations representing extreme cases: a pure turbulent deflagration, a delayed detonation following an approximately spherical ignition of the initial deflagration, and a delayed detonation arising from a highly asymmetric deflagration ignition. Apart from this initial condition, the deflagration stage is treated in a parameter-free approach. The detonation is initiated when the turbulent burning enters the distributed burning regime. This occurs at densities around $10^{7}$ g cm$^{-3}$ -- relatively low as compared to existing nucleosynthesis studies for one-dimensional spherically symmetric models. The burning in these multidimensional models is different from that in one-dimensional simulations as the detonation wave propagates both into unburned material in the high density region near the center of a white dwarf and into the low density region near the surface. Thus, the resulting yield is a mixture of different explosive burning products, from carbon-burning products at low densities to complete silicon-burning products at the highest densities, as well as electron-capture products synthesized at the deflagration stage. In contrast to the deflagration model, the delayed detonations produce a characteristic layered structure and the yields largely satisfy constraints from Galactic chemical evolution. In the asymmetric delayed detonation model, the region filled with electron capture species (e.g., $^{58}$Ni, $^{54}$Fe) is within a shell, showing a large off-set, above the bulk of $^{56}$Ni distribution, while species produced by the detonation are distributed more spherically (abridged).

The delayed detonation model describes the observational properties of the majority of type Ia supernovae very well. Using numerical data from a three-dimensional deflagration model for type Ia supernovae, the intermittency of the turbulent velocity field and its implications on the probability of a deflagration-to-detonation (DDT) transition are investigated. From structure functions of the turbulent velocity fluctuations, we determine intermittency parameters based on the log-normal and the log-Poisson models. On the other hand, the analysis of the turbulent velocity fluctuations in the vicinity of the flame front by Roepke suggests a much higher probability of large velocity fluctuations on the grid scale in comparison to the log-normal intermittency model. Following Pan et al., we computed probability density functions for a DDT for the different distributions. Assuming that a DDT can occur in the stirred flame regime, as proposed by Woosley et al., the log-normal model would imply a delayed detonation between 0.7 and 0.8 seconds after the beginning of the deflagration phase for the multi-spot ignition scenario used in the simulation. However, the probability drops to virtually zero if a DDT is further constrained by the requirement that the turbulent velocity fluctuations reach about 500 km/s. Under this condition, delayed detonations are only possible if the distribution of the velocity fluctuations is not log-normal. From our calculations follows that the distribution obtained by Roepke allow for multiple DDTs around 0.8 seconds after ignition at a transition density close to 1x10^7 g/cm^3.

The theory of radiative transfer provides the link between the physical conditions in an astrophysical object and the observable radiation which it emits. Thus accurately modelling radiative transfer is often a necessary part of testing theoretical models by comparison with observations. We describe a new radiative transfer code which employs Monte Carlo methods for the numerical simulation of radiation transport in expanding media. We discuss the application of this code to the calculation of synthetic spectra and light curves for a Type Ia supernova explosion model and describe the sensitivity of the results to certain approximations made in the simulations.

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.

Despite the importance of Type Ia supernovae as standard candles for cosmology and to the chemical evolution of the Universe, we still have no consistent picture of the nature of these events. Much progress has been made in the hydrodynamical explosion modelling of supernovae Ia in the last few years and fully 3-D explosion models are now available. However those simulations are not directly comparable to observations: to constrain explosion models, radiative transfer calculations must be carried out. We present a new 3-D Monte Carlo radiative transfer code which allows forward modelling of the spectral evolution of Type Ia supernovae from first principles, using hydrodynamical explosion models as input. Here, as a first application, we calculate line-of-sight dependent colour light curves for a toy model of an off-centre explosion.

In order to assess qualitatively the ejecta geometry of stripped-envelope core-collapse supernovae, we investigate 98 late-time spectra of 39 objects, many of them previously unpublished. We perform a Gauss-fitting of the [O I] 6300, 6364 feature in all spectra, with the position, full width at half maximum (FWHM) and intensity of the 6300 Gaussian as free parameters, and the 6364 Gaussian added appropriately to account for the doublet nature of the [O I] feature. On the basis of the best-fit parameters, the objects are organised into morphological classes, and we conclude that at least half of all Type Ib/c supernovae must be aspherical. Bipolar jet-models do not seem to be universally applicable, as we find too few symmetric double-peaked [O I] profiles. In some objects the [O I] line exhibits a variety of shifted secondary peaks or shoulders, interpreted as blobs of matter ejected at high velocity and possibly accompanied by neutron-star kicks to assure momentum conservation. At phases earlier than ~200d, a systematic blueshift of the [O I] 6300, 6364 line centroids can be discerned. Residual opacity provides the most convincing explanation of this phenomenon, photons emitted on the rear side of the SN being scattered or absorbed on their way through the ejecta. Once modified to account for the doublet nature of the oxygen feature, the profile of Mg I] 4571 at sufficiently late phases generally resembles that of [O I] 6300, 6364, suggesting negligible contamination from other lines and confirming that O and Mg are similarly distributed within the ejecta.

We analyze the statistical properties of the turbulent velocity field in the deflagration model for Type Ia supernovae. In particular, we consider the question of whether turbulence is isotropic and consistent with the Kolmogorov theory at small length scales. Using numerical data from a high-resolution simulation of a thermonuclear supernova explosion, spectra of the turbulence energy and velocity structure functions are computed. We show that the turbulent velocity field is isotropic at small length scales and follows a scaling law that is consistent with the Kolmogorov theory until most of the nuclear fuel is burned. At length scales greater than a certain characteristic scale, turbulence becomes anisotropic. Here, the radial velocity fluctuations follow the scaling law of the Rayleigh-Taylor instability, whereas the angular component still obeys Kolmogorov scaling. In the late phase of the explosion, this characteristic scale drops below the numerical resolution of the simulation. The analysis confirms that a subgrid-scale model for the unresolved turbulence energy is required for the consistent calculation of the flame speed in deflagration models of Type Ia supernovae, and that the assumption of isotropy on these scales is appropriate.

The nature of Type Ia supernova progenitors is still unclear. The outstanding characteristic of the single-degenerate scenario is that it contains hydrogen in the binary companion of the exploding white dwarf star, which, if mixed into the ejecta of the supernova in large amounts may lead to conflicts with the observations thus ruling out the scenario. We investigate the effect of the impact of Type Ia supernova ejecta on a main sequence companion star of the progenitor system. With a series of simulations we investigate how different parameters of this system affect the amount of hydrogen stripped from the companion by the impact. The stellar evolution code GARSTEC is used to set up the structure of the companion stars mimicking the effect of a binary evolution phase. The impact itself is simulated with the smoothed particle hydrodynamics code GADGET2. We reproduce and confirm the results of earlier grid-based hydrodynamical simulation. Parameter studies of the progenitor system are extended to include the results of recent binary evolution studies. The more compact structure of the companion star found here significantly reduces the stripped hydrogen mass. The low hydrogen masses resulting from a more realistic companion structure are consistent with current observational constraints. Therefore, the single-degenerate scenario remains a valid possibility for Type Ia supernova progenitors. These new results are not a numerical effect, but the outcome of different initial conditions.

SNe Ia are good distance indicators because the shape of their light curves, which can be measured independently of distance, varies smoothly with luminosity. This suggests that SNe Ia are a single family of events. Similar correlations are observed between luminosity and spectral properties. In particular, the ratio of the strengths of the SiII \lambda 5972 and \lambda 6355 lines, known as R(SiII), was suggested as a potential luminosity indicator. Here, the physical reasons for the observed correlation are investigated. A Monte-Carlo code is used to construct a sequence of synthetic spectra resembling those of SNe with different luminosities near B maximum. The influence of abundances and of ionisation and excitation conditions on the synthetic spectral features is investigated. The ratio R(SiII) depends ssentially on the strength of SiII \lambda 5972, because SiII \lambda 6355 is saturated. In less luminous objects, SiII \lambda 5972 is stronger because of a rapidly increasing SiII/SiIII ratio. Thus, the correlation between R(SiII) and luminosity is the effect of ionisation balance. The SiII \lambda 5972 line itself may be the best spectroscopic luminosity indicator for SNe Ia, but all indicators discussed show scatter which may be related to abundance distributions.

As part of the European Supernova Collaboration we obtained extensive photometry and spectroscopy of the type Ia SN 2002dj covering epochs from 11 days 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-IR colours and a reduced spectral emission beyond 6500 A. For high-velocity gradient Type Ia supernovae, 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 objects. 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 Type Ia supernovae 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.

We present the spectra of 36 Supernovae (SNe) of various types, obtained by the European Supernova Collaboration. Because of the spectral classification and the phase determination at their discovery the SNe did not warrant further study, and the spectra we present are the only available for the respective objects. In this paper we present and discuss this material using a new software for the automated classification of SNe spectra. As a validation of the software, we verify the classification and phase estimate reported for these objects in their discovery / classification circulars. For the comparison, the software uses the library of template spectra of Padova-Asiago Supernova Archive (ASA). For each spectrum of our sample we present a brief, individual discussion, highlighting the main characteristics and possible peculiarities. The comparison with ASA spectra confirms the previous classification of all objects and refines the age estimates. For our software we determine numerical limits of "safe" spectral classification and the uncertainties of the phase determination.

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.

Optical observations of the Type Ia supernova (SN Ia) 2005bl in NGC 4070, obtained from -6 to +66 d with respect to the B-band maximum, are presented. The photometric evolution is characterised by rapidly-declining light curves and red colours at peak and soon thereafter. With M_B,max = -17.24 the SN is an underluminous SN Ia, similar to the peculiar SNe 1991bg and 1999by. This similarity also holds for the spectroscopic appearance, the only remarkable difference being the likely presence of carbon in pre-maximum spectra of SN 2005bl. A comparison study among underluminous SNe Ia is performed, based on a number of spectrophotometric parameters. Previously reported correlations of the light-curve decline rate with peak luminosity and R(Si) are confirmed, and a large range of post-maximum Si II lambda6355 velocity gradients is encountered. 1D synthetic spectra for SN 2005bl are presented, which confirm the presence of carbon and suggest an overall low burning efficiency with a significant amount of leftover unburned material. Also, the Fe content in pre-maximum spectra is very low, which may point to a low metallicity of the precursor. Implications for possible progenitor scenarios of underluminous SNe Ia are briefly discussed.

The onset of the thermonuclear runaway in a Chandrasekhar-mass white dwarf, leading to the explosion as a type Ia supernova, is studied with hydrodynamical simulations. We investigate the evolution of temperature fluctuations (``bubbles'') in the WD's convective core by means of 2D numerical simulations. We show how the occurrence of the thermonuclear runaway depends on various bubble parameters. The relevance of the progenitor's composition for the ignition process is also discussed.

In the "double-detonation sub-Chandrasekhar" model for type Ia supernovae, a carbon-oxygen (C + O) white dwarf accumulates sufficient amounts of helium such that a detonation ignites in that layer before the Chandrasekhar mass is reached. This detonation is thought to trigger a secondary detonation in the C + O core. By means of one- and two-dimensional hydrodynamic simulations, we investigate the robustness of this explosion mechanism for generic 1-M_sun models and analyze its observable predictions. Also a resolution dependence in numerical simulations is analyzed. The propagation of thermonuclear detonation fronts, both in helium and in the carbon-oxygen mixture, is computed by means of both a level-set function and a simplified description for nuclear reactions. The decision whether a secondary detonation is triggered in the white dwarf's core or not is made based on criteria given in the literature. In a parameter study involving different initial flame geometries for He-shell masses of 0.2 and 0.1 M_sun, we find that a secondary detonation ignition is a very robust process. Converging shock waves originating from the detonation in the He shell generate the conditions for a detonation near the center of the white dwarf in most of the cases considered. Finally, we follow the complete evolution of three selected models with 0.2 M_sun of He through the C/O-detonation phase and obtain nickel-masses of about 0.40 to 0.45 M_sun. Although we have not done a complete scan of the possible parameter space, our results show that sub-Chandrasekhar models are not good candidates for normal or sub-luminous type Ia supernovae. The chemical composition of the ejecta features significant amounts of nickel in the outer layers at high expansion velocities, which is inconsistent with near-maximum spectra. (abbreviated)

We present VRIJHK photometry, and optical and near-infrared spectroscopy, of the heavily extinguished Type Ia supernova (SN) 2002cv, located in NGC 3190, which is also the parent galaxy of the Type Ia SN 2002bo. SN 2002cv, not visible in the blue, has a total visual extinction of 8.74 +- 0.21 mag. In spite of this we were able to obtain the light curves between -10 and +207 days from the maximum in the I band, and also to follow the spectral evolution, deriving its key parameters. We found the peak I-band brightness to be Imax = 16.57 +- 0.10 mag, the maximum absolute I magnitude to be MmaxI = -18.79 +- 0.20, and the parameter dm15(B) specifying the width of the B-band light curve to be 1.46 +- 0.17 mag. The latter was derived using the relations between this parameter and dm40(I) and the time interval dtmax(I) between the two maxima in the I-band light curve. As has been found for previously observed, highly extinguished SNe Ia, a small value of 1.59 +- 0.07 was obtained here for the ratio Rv of the total-to-selective extinction ratio for SN 2002cv, which implies a small mean size for the grains along the line of sight toward us. Since it was found for SN 2002bo a canonical value of 3.1, here we present a clear evidence of different dust properties inside NGC 3190.

Aims. We study the optical spectroscopic properties of Type Ia Supernova (SN Ia) 2004dt, focusing our attention on the early epochs. Methods. Observation triggered soon after the SN 2004dt discovery allowed us to obtain a spectrophotometric coverage from day -10 to almost one year (~353 days) after the B band maximum. Observations carried out on an almost daily basis allowed us a good sampling of the fast spectroscopic evolution of SN 2004dt in the early stages. To obtain this result, low-resolution, long-slit spectroscopy was obtained using a number of facilities. Results. This supernova, which in some absorption lines of its early spectra showed the highest degree of polarization ever measured in any SN Ia, has a complex velocity structure in the outer layers of its ejecta. Unburnt oxygen is present, moving at velocities as high as ~16,700 km/s, with some intermediate-mass elements (Mg, Si, Ca) moving equally fast. Modeling of the spectra based on standard density profiles of the ejecta fails to reproduce the observed features, whereas enhancing the density of outer layers significantly improves the fit. Our analysis indicates the presence of clumps of high-velocity, intermediate-mass elements in the outermost layers, which is also suggested by the spectropolarimetric data.

Type Ia supernovae are thought to be thermonuclear explosions of accreting white dwarfs that reach a critical mass limit. Despite their importance as cosmological distance indicators, the nature of their progenitors has remained controversial. Here we report the detection of circumstellar material in a normal Type Ia supernova. The expansion velocities, densities and dimensions of the circumstellar envelope indicate that this material was ejected from the progenitor system. The relatively low expansion velocities appear to favor a progenitor system where a white dwarf accretes material from a companion star which is in the red-giant phase at the time of explosion.

A simulation of the thermonuclear explosion of a Chandrasekhar-mass C+O white dwarf, the most popular scenario of a type Ia supernova (SN Ia), is presented. The underlying modeling is pursued in a self-consistent way, treating the combustion wave as a turbulent deflagration using well tested methods developed for laboratory combustion and based on the concept of `large eddy simulations' (LES). Such consistency requires to capture the onset of the turbulent cascade on resolved scales. This is achieved by computing the dynamical evolution on a 1024$^3$ moving grid, which resulted in the best-resolved three-dimensional SN Ia simulation carried out thus far, reaching the limits of what can be done on present supercomputers. Consequently, the model has no free parameters other than the initial conditions at the onset of the explosion, and therefore it has considerable predictive power. Our main objective is to determine to which extent such a simulation can account for the observations of normal SNe Ia. Guided by previous simulations with less resolution and a less sophisticated flame model, initial conditions were chosen that yield a reasonably strong explosion and a sufficient amount of radioactive nickel for a bright display. We show that observables are indeed matched to a reasonable degree. In particular, good agreement is found with the light curves of normal SNe Ia. Moreover, the model reproduces the general features of the abundance stratification as inferred from the analysis of spectra. This indicates that it captures the main features of the explosion mechanism of SNe Ia. However, we also show that even a seemingly best-choice pure deflagration model has shortcomings that indicate the need for a different mode of nuclear burning at late times, perhaps the transition to a detonation at low density.

An extensive set of optical and NIR photometry and low-resolution spectra the Type Ia supernova (SN Ia) 2003du was obtained using a number of facilities. The observations started 13 days before B-band maximum light and continued for 480 days with exceptionally good time sampling. The optical photometry was calibrated through the S-correction technique. The UBVRIJHK light curves and the color indices of SN 2003du closely resemble those of normal SNe Ia. SN 2003du reached a B-band maximum of 13.49 (+/-0.02) mag on JD2452766.38 (+/-0.5). We derive a B-band stretch parameter of 0.988 (+/-0.003), which corresponds to dM15=1.02 (+/-0.05), indicative of a SN Ia of standard luminosity. The reddening in the host galaxy was estimated by three methods, and was consistently found to be negligible. We estimate a distance modulus of 32.79 (+/-0.15) mag to the host galaxy, UGC 9391. The peak UVOIR bolometric luminosity of 1.35(+/-0.20) 10^43 erg/s and Arnett's rule implies that M(Ni56)=0.68 (+/-0.14)M_sun of Ni56 was synthesized during the explosion. Modeling of the UVOIR bolometric light curve also indicates M(Ni56) in the range 0.6-0.8 M_sun. The spectral evolution of SN 2003du at both optical and NIR wavelengths also closely resembles normal SNe Ia. In particular, the Si II ratio at maximum R(Si II)=0.22 (+/-0.02) and the time evolution of the blueshift velocities of the absorption line minima are typical. The pre-maximum spectra of SN 2003du showed conspicuous high-velocity features in the Ca II H&K doublet and infrared triplet, and possibly in Si II 6355, lines. We compare the time evolution of the profiles of these lines with other well-observed SNe Ia and we suggest that the peculiar pre-maximum evolution of Si II 6355 line in many SNe Ia is due to the presence of two blended absorption components.

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.

The ESC-RTN optical spectroscopy data-set for SN 2005cf is presented and analyzed. The observations range from -11.6 and +77.3 days with respect to B-band maximum light. The evolution of the spectral energy distribution of SN 2005cf is characterized by the presence of high velocity SiII and CaII features. SYNOW synthetic spectra are used to investigate the ejecta geometry of silicon. Based on the synthetic spectra the SiII high velocity feature appears detached at 19500 km/s. We also securely establish the presence of such feature in SN 1990N, SN 1994D, SN 2002er and SN 2003du. On a morphological study both the CaII IR Triplet and H&K absorption lines of SN 2005cf show high velocity features centered around 24000 km/s. When compared with other Type Ia SNe based on the scheme presented in Benetti et al. 2005 SN 2005cf definitely belongs to the LVG group.

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 SN Ia. However, the ejected 56Ni mass derived by modelling the bolometric light curve (about 0.45 solar masses) lies near the lower limit of the 56Ni mass distribution observed in normal SNe Ia. Accordingly, SN 2004eo shows a relatively rapid post-maximum decline in the light curve (Delta m_(B) = 1.46), small expansion velocities in the ejecta, and a depth ratio Si II 5972 / Si II 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. (2005). Similar behaviour is seen in a few other SNe Ia. Thus, there may in fact exist a few SNe Ia with intermediate physical properties.

We present early-time optical and near-infrared photometry of supernova (SN) 2005cf. The observations, spanning a period from about 12 days before to 3 months after maximum, have been obtained through the coordination of observational efforts of various nodes of the European Supernova Collaboration and including data obtained at the 2m Himalayan Chandra Telescope. From the observed light curve we deduce that SN 2005cf is a fairly typical SN Ia with a post-maximum decline (Delta m_15(B) = 1.12) close to the average value and a normal luminosity of M_B,max = -19.39+/-0.33. Models of the bolometric light curve suggest a synthesised 56Ni mass of about 0.7 solar masses. The negligible host galaxy interstellar extinction and its proximity make SN 2005cf a good Type Ia supernova template.

The recent discovery of a very bright type Ia supernova, SNLS-03D3bb ($\equiv$ SN 2003fg), in the Supernova Legacy Survey (SNLS) has raised the question of whether super-Chandrasekhar-mass white dwarf stars are needed to explain such explosions. In principle, such a progenitor could form by the mergers of two rather massive white dwarfs. Binary systems of two white dwarfs in close orbit where their total mass significantly exceeds the Chandrasekhar mass, have not yet been found. Therefore SNLS-03D3bb could establish the first clear case of a double-degenerate progenitor of a (peculiar) type Ia supernovae. Moreover, if this interpretation is correct, it casts some doubt on the universality of the calibration relations used to make SNe Ia distance indicators for cosmology. Here we present a critical discussion of the expected observational finger prints of super-Chandrasekhar-mass explosions; in important respects, these are not consistent with the observations of SNLS-03D3bb. We demonstrate that the lop-sided explosion of a Chandrasekhar-mass white dwarf could provide a better explanation.

Type Ia supernovae, the thermonuclear explosions of white dwarf stars composed of carbon and oxygen, were instrumental as distance indicators in establishing the acceleration of the universe's expansion. However, the physics of the explosion are debated. Here we report a systematic spectral analysis of a large sample of well observed type Ia supernovae. Mapping the velocity distribution of the main products of nuclear burning, we constrain theoretical scenarios. We find that all supernovae have low-velocity cores of stable iron-group elements. Outside this core, nickel-56 dominates the supernova ejecta. The outer extent of the iron-group material depends on the amount of nickel-56 and coincides with the inner extent of silicon, the principal product of incomplete burning. The outer extent of the bulk of silicon is similar in all SNe, having an expansion velocity of ~11000 km/s and corresponding to a mass of slightly over one solar mass. This indicates that all the supernovae considered here burned similar masses, and suggests that their progenitors had the same mass. Synthetic light curve parameters and three-dimensional explosion simulations support this interpretation. A single explosion scenario, possibly a delayed detonation, may thus explain most type Ia supernovae.

We present a preliminary analysis of an extensive set of optical observations of the Type Ia SN 2005hk. We show that the evolution of SN 2005hk closely follows that of the peculiar SN 2002cx. SN 2005hk is more luminous than SN 2002cx, while still under-luminous compared to normal Type Ia supernovae. The spectrum at 9 days before maximum is dominated by conspicuous Fe III and Ni III lines, and the Si II 6355 line is also clearly visible. All these features have low velocity (~6000 km/s). The near maximum spectra show lines of Si II, S II, Fe II, Fe III, as well as other intermediate mass and iron group elements. Analysis with the code for synthetic spectra SYNOW indicates that all these spectral lines have similar velocities.

Modeling type Ia supernova (SN Ia) explosions in three dimensions allows to eliminate any undetermined parameters and provides predictive power to simulations. This is necessary to improve the understanding of the explosion mechanism and to settle the question of the applicability of SNe Ia in cosmological distance measurements. Since the models contain no tunable parameters, it is also possible to directly assess their validity on the basis of a comparison with observations. Here, we describe the modeling of SNe Ia as thermonuclear explosions in which the flame after ignition near the center of the progenitor white dwarf star propagates outward in the sub-sonic deflagration mode accelerated by the interaction with turbulence. We explore the capabilities of this model by comparison with observations and show in a preliminary approach, how such a model can be applied to study the origin of the diversity of SNe Ia.

We study intermediate--z SNe Ia using the empirical physical diagrams which enable to learn about those SNe explosions. This information can be very useful to reduce systematic uncertainties of the Hubble diagram of SNe Ia up to high z. The study of the expansion velocities and the measurement of the ratio $\mathcal{R}$(\SiII) allow to subtype those SNe Ia as done for nearby samples. The evolution of this ratio as seen in the diagram $\mathcal{R}$(\SiII)--(t) together with $\mathcal{R}$(\SiII)$_{max}$ versus (B-V)$_{0}$ indicate consistency of the properties at intermediate z compared with local SNe. At intermediate--z, the expansion velocities of Ca II and Si II are similar to the nearby counterparts. This is found in a sample of 6 SNe Ia in the range 0.033$\leq z \leq$0.329 discovered within the {\it International Time Programme} (ITP) of {\it Cosmology and Physics with SNe Ia} during the spring of 2002. Those supernovae were identified using the 4.2m William Herschel Telescope. Two SNe Ia at intermediate z were of the cool FAINT type, one being a SN1986G--like object highly reddened. The $\mathcal{R}$(\SiII) ratio as well as subclassification of the SNe Ia beyond templates help to place SNe Ia in their sequence of brightness and to distinguish between reddened and intrinsically red supernovae. This test can be done with very high z SNe Ia and it will help to reduce systematic uncertainties due to extinction by dust. It should allow to map the high-z sample into the nearby one.

Type Ia supernovae (SNe Ia) are one of the major tools to determine the cosmological parameters. Utilizing them as distance indicators, it is possible to geometrically survey the universe. To this end, the intrinsic scatter in the luminosities of these events needs to be calibrated using empirical relations between observables. A theoretical explanation for these relations is still lacking and can only be provided by a sound understanding of the mechanism of SNe Ia. Recently there has been significant progress in modeling SNe Ia. We report on numerical simulations of the explosion process, compare the results with observations of nearby SNe Ia, and discuss current uncertainties of the models. The presented simulations shed some light on the origin of the diversity of SNe Ia. Such simulations will pave the way towards an understanding of SN Ia diversities and correlations of their properties and ultimately provide a tool to validate the cosmological implications of SN Ia distance measurements.

We present results from a systematic study of the effects of initial parameters on three-dimensional thermonuclear supernova models.

We review the astrophysical modeling of type Ia supernova explosions and describe numerical methods to implement numerical simulations of these events. Some results of such simulations are discussed.

We present a systematic study of the diversity of three-dimensional deflagration simulations of Type Ia supernova explosions arising from variations of the initial parameters. By altering the carbon mass fraction, the central density, and the metallicity of the progenitor white dwarf star, we set up a grid of numerical explosion models. While changing the central density has the largest impact on the explosion energy, the largest variation in the 56Ni production is found by changing the metallicity of the models. Varying the carbon mass fraction hardly affects the 56Ni synthesized although it alters the energetics of the explosion. Possible consequences for the shape of light curves of Type Ia supernovae are discussed.

The explosion of a carbon-oxygen white dwarf as a Type Ia supernova is known to be sensitive to the manner in which the burning is ignited. Studies of the pre-supernova evolution suggest asymmetric, off-center ignition, and here we explore its consequences in two- and three-dimensional simulations. Compared with centrally ignited models, one-sided ignitions initially burn less and release less energy. For the distributions of ignition points studied, ignition within two hemispheres typically leads to the unbinding of the white dwarf, while ignition within a small fraction of one hemisphere does not. We also examine the spreading of the blast over the surface of the white dwarf that occurs as the first plumes of burning erupt from the star. In particular, our studies test whether the collision of strong compressional waves can trigger a detonation on the far side of the star as has been suggested by Plewa et al. (2004). The maximum temperature reached in these collisions is sensitive to how much burning and expansion has already gone on, and to the dimensionality of the calculation. Though detonations are sometimes observed in 2D models, none ever happens in the corresponding 3D calculations. Collisions between the expansion fronts of multiple bubbles also seem, in the usual case, unable to ignite a detonation. "Gravitationally confined detonation" is therefore not a robust mechanism for the explosion. Detonation may still be possible in these models however, either following a pulsation or by spontaneous detonation if the turbulent energy is high enough.

Optical and near-infrared observations of the Type Ic supernova (SN) 2004aw are presented, obtained from day -3 to day +413 with respect to the B-band maximum. The photometric evolution is characterised by a comparatively slow post-maximum decline of the light curves. The peaks in redder bands are significantly delayed relative to the bluer bands, the I-band maximum occurring 8.4 days later than that in B. With an absolute peak magnitude of -18.02 in the V band the SN can be considered fairly bright, but not exceptional. This also holds for the U through I bolometric light curve, where SN 2004aw has a position intermediate between SNe 2002ap and 1998bw. Spectroscopically SN 2004aw provides a link between a normal Type Ic supernova like SN 1994I and the group of broad-lined SNe Ic. The spectral evolution is rather slow, with a spectrum at day +64 being still predominantly photospheric. The shape of the nebular [O I] 6300,6364 line indicates a highly aspherical explosion. Helium cannot be unambiguously identified in the spectra, even in the near-infrared. Using an analytical description of the light curve peak we find that the total mass of the ejecta in SN 2004aw is 3.5-8.0 M_Sun, significantly larger than in SN 1994I, although not as large as in SN 1998bw. The same model suggests that about 0.3 M_Sun of {56}Ni has been synthesised in the explosion. No connection to a GRB can be firmly established.

We present optical and near-infrared photometry and spectroscopy of the Type Ia SN 2003cg, which exploded in the nearby galaxy NGC 3169. The observations cover a period between -8.5 and +414 days post-maximum. SN 2003cg is a normal but highly-reddened Type Ia event. Its B magnitude at maximum B_max = 15.94+/-0.04 and Delta m_15(B)_obs = 1.12+/-0.04 (Delta m_15(B)_intrinsic = 1.25+/-0.05). Allowing R_v to become a free parameter within the Cardelli et al. (1989) extinction law, simultaneous matches to a range of colour curves of normal SNe Ia yielded E(B-V) = 1.33+/-0.11, and R_v = 1.80+/-0.19. While the value obtained for R_v is small, such values have been invoked in the past, and may imply a grain size which is small compared with the average value for the local ISM.

We present synthetic bolometric and broad-band UBVRI light curves of SNe Ia, for four selected 3-D deflagration models of thermonuclear supernovae. The light curves are computed with the 1-D hydro code STELLA, which models (multi-group time-dependent) non-equilibrium radiative transfer inside SN ejecta. Angle-averaged results from 3-D hydrodynamical explosion simulations with the composition determined in a nucleosynthetic postprocessing step served as the input to the radiative transfer model. The predicted model UBV light curves do agree reasonably well with the observed ones for SNe Ia in the range of low to normal luminosities, although the underlying hydrodynamical explosion models produced only a modest amount of radioactive Ni56 and relatively low kinetic energy in the explosion. The evolution of predicted B and V fluxes in the model with a Ni56 mass of 0.42 M_sun follows the observed decline rate after the maximum very well, although the behavior of fluxes in other filters somewhat deviates from observations, and the bolometric decline rate is a bit slow. Using our models, we check the validity of Arnett's rule and the accuracy of the procedure for extracting the Ni56 mass from the observed light curves. We find that the comparison between theoretical light curves and observations provides a useful tool to validate SN Ia models. The steps necessary to improve the agreement between theory and observations are set out.

We present a one-equation subgrid scale model that evolves the turbulence energy corresponding to unresolved velocity fluctuations in large eddy simulations. The model is derived in the context of the Germano consistent decomposition of the hydrodynamical equations. The eddy-viscosity closure for the rate of energy transfer from resolved toward subgrid scales is localised by means of a dynamical procedure for the computation of the closure parameter. Therefore, the subgrid scale model applies to arbitrary flow geometry and evolution. For the treatment of microscopic viscous dissipation a semi-statistical approach is used, and the gradient-diffusion hypothesis is adopted for turbulent transport. A priori tests of the localised eddy-viscosity closure and the gradient-diffusion closure are made by analysing data from direct numerical simulations. As an a posteriori testing case, the large eddy simulation of thermonuclear combustion in forced isotropic turbulence is discussed. We intend the formulation of the subgrid scale model in this paper as a basis for more advanced applications in numerical simulations of complex astrophysical phenomena involving turbulence.

The dynamics of the explosive burning process is highly sensitive to the flame speed model in numerical simulations of type Ia supernovae. Based upon the hypothesis that the effective flame speed is determined by the unresolved turbulent velocity fluctuations, we employ a new subgrid scale model which includes a localised treatment of the energy transfer through the turbulence cascade in combination with semi-statistical closures for the dissipation and non-local transport of turbulence energy. In addition, subgrid scale buoyancy effects are included. In the limit of negligible energy transfer and transport, the dynamical model reduces to the Sharp-Wheeler relation. According to our findings, the Sharp-Wheeler relation is insuffcient to account for the complicated turbulent dynamics of flames in thermonuclear supernovae. The application of a co-moving grid technique enables us to achieve very high spatial resolution in the burning region. Turbulence is produced mostly at the flame surface and in the interior ash regions. Consequently, there is a pronounced anisotropy in the vicinity of the flame fronts. The localised subgrid scale model predicts significantly enhanced energy generation and less unburnt carbon and oxygen at low velocities compared to earlier simulations.

We present results of simulations of stellar collapse and explosions in spherical symmetry for progenitor stars in the 8-10 solar mass range with an O-Ne-Mg core. The simulations were continued until nearly one second after core bounce and were performed with the Prometheus/Vertex code with a variable Eddington factor solver for the neutrino transport, including a state-of-the-art treatment of neutrino-matter interactions. Particular effort was made to implement nuclear burning and electron capture rates with sufficient accuracy to ensure a smooth continuation, without transients, from the progenitor evolution to core collapse. Using two different nuclear equations of state (EoSs), a soft version of the Lattimer & Swesty EoS and the significantly stiffer Wolff & Hillebrandt EoS, we found no prompt explosions, but instead delayed explosions, powered by neutrino heating and the neutrino-driven baryonic wind which sets in about 200 ms after bounce. The models eject little nickel (< 0.015 solar masses), explode with an energy of about or slightly more than 10**50 erg, and leave behind neutron stars (NSs) with a baryonic mass near 1.36 solar masses. Different from previous models of such explosions, the ejecta during the first second have a proton-to-baryon ratio of Ye > 0.46, which suggests a chemical composition that is not in conflict with galactic abundances. No low-entropy matter with Ye << 0.5 is ejected. This excludes such explosions as sites of a low-entropy r-process. The low explosion energy and nucleosynthetic implications are compatible with the observed properties of the Crab supernova, and the small nickel mass supports the possibility that our models explain some subluminous Type II-P supernovae.

In the framework of the Chandrasekhar-mass deflagration model for Type Ia supernovae (SNe Ia), a persisting free parameter is the initial morphology of the flame front, which is linked to the ignition process in the progenitor white dwarf. Previous analytical models indicate that the thermal runaway is driven by temperature perturbations (''bubbles'') that develop in the white dwarf's convective core. In order to probe the conditions at ignition (diameters, temperatures and evolutionary timescales), we have performed hydrodynamical 2D simulations of buoyant bubbles in white dwarf interiors. Our results show that fragmentation occurring during the bubble rise affects the outcome of the bubble evolution. Possible implications for the ignition process of SNe Ia are discussed.

We present a systematic survey of the capabilities of type Ia supernova explosion models starting from a number of flame seeds distributed around the center of the white dwarf star. To this end we greatly improved the resolution of the numerical simulations in the initial stages. This novel numerical approach facilitates a detailed study of multi-spot ignition scenarios with up to hundreds of ignition sparks. Two-dimensional simulations are shown to be inappropriate to study the effects of initial flame configurations. Based on a set of three-dimensional models, we conclude that multi-spot ignition scenarios may improve type Ia supernova models towards better agreement with observations. The achievable effect reaches a maximum at a limited number of flame ignition kernels as shown by the numerical models and corroborated by a simple dimensional analysis.

The abundance stratification in the ejecta of the normal Type Ia supernova 2002bo is derived by fitting a series of spectra obtained at close time intervals. A Monte Carlo code, modified to include abundance stratification, is used to compute synthetic spectra at 13 epochs in the photospheric phase, starting 13d before B maximum. A description of the abundance distribution above 7600 km/s is thus obtained. Abundances in deeper layers, down to zero velocity, are derived from models of two nebular-phase spectra. Elements synthesized in different stages of burning are significantly, but not completely, mixed in the ejecta. A total 56Ni mass of 0.52M_sun is derived. Evidence for intermediate-mass elements at high velocities (> 18 000 km/s) is found, most clearly in Si II 6355A, in Ca II H&K, and in the Ca II infrared triplet. Carbon lines are not seen at any velocity, with possible implications for the progenitor/explosion scenario. A synthetic bolometric light curve computed using the inferred abundance distribution is in very good agreement with the observed one, providing an independent check. In particular, the fast rise of the light curve is reproduced very well. This is due to outward mixing of 56Ni, which is clearly well determined by the spectroscopic modelling.

We study the dynamics of thermonuclear flames propagating in fuel stirred by stochastic forcing. The fuel consists of carbon and oxygen in a state which is encountered in white dwarfs close to the Chandrasekhar limit. The level set method is applied to represent the flame fronts numerically. The computational domain for the numerical simulations is cubic, and periodic boundary conditions are imposed. The goal is the development of a suitable flame speed model for the small-scale dynamics of turbulent deflagration in thermonuclear supernovae. Because the burning process in a supernova explosion is transient and spatially inhomogeneous, the localised determination of subgrid scale closure parameters is essential. We formulate a semi-localised model based on the dynamical equation for the subgrid scale turbulence energy $k_{\mathrm{sgs}}$. The turbulent flame speed $s_{\mathrm{t}}$ is of the order $\sqrt{2k_{\mathrm{sgs}}}$. In particular, the subgrid scale model features a dynamic procedure for the calculation of the turbulent energy transfer from resolved toward subgrid scales, which has been successfully applied to combustion problems in engineering. The options of either including or suppressing inverse energy transfer in the turbulence production term are compared. In combination with the piece-wise parabolic method for the hydrodynamics, our results favour the latter option. Moreover, different choices for the constant of proportionality in the asymptotic flame speed relation, $s_{\mathrm{t}}\propto\sqrt{2k_{\mathrm{sgs}}}$, are investigated.

We investigate the metallicity effect (measured by the original 22Ne content) on the detailed nucleosynthetic yields for 3D hydrodynamical simulations of the thermonuclear burning phase in SNe Ia. Calculations are based on post-processes of the ejecta, using passively advected tracer particles, as explained in details by Travaglio et al.(2004). The nuclear reaction network employed in computing the explosive nucleosynthesis contains 383 nuclear species. For this work we use the high resolution multi-point ignition (bubbles) model b30_3d_768 (Travaglio et al.2004 for the solar metallicity case), and we cover a metallicity range between 0.1xZ_sun up to 3xZ_sun. We find a linear dependence of the 56Ni mass ejected on the progenitor's metallicity, with a variation in the 56Ni mass of ~25% in the metallicity range explored. Moreover, the largest variation in 56Ni occurs at metallicity greater than solar. Almost no variations are shown in the unburned material 12C and 16O. The largest metallicity effect is seen in the alpha-elements. Finally, implications for the observed scatter in the peak luminosities of SNe Ia are also discussed.

The use of type Ia supernovae as distance indicators for cosmology has initiated a search for theoretical arguments supporting the empirical calibration methods applied. To this end, as a first step, a sound understanding of the origin of the observed diversity in type Ia supernova properties is needed. Here we present a first systematic study of effects resulting from changing some physical parameters of three-dimensional deflagration models of thermonuclear supernovae. In our study we vary the progenitor's carbon-to-oxygen ratio and its central density prior to ignition because both properties are not well determined by stellar evolution theory and they may change from supernova to supernova. Next we compute for these explosion models the nucleosynthesis yields in a post-processing step. This, in addition, allows us to study variations in the progenitor's metallicity by means of different $^{22}$Ne mass fractions in the initial composition. We find that the progenitor's carbon-to-oxygen ratio and its central density affect the energy release of the models and thus the expansion velocity of the supernova. Moreover, we find that changing the metallicity and the central density changes the production of radioactive $^{56}$Ni and thus affects the luminosity. In contrast, the carbon-to-oxygen ratio has little effect on the $^{56}$Ni production. Implications of the found variations of the explosion energy and the produced $^{56}$Ni mass for the type Ia supernova diversity are discussed.

Late time synthetic spectra of Type Ia supernovae, based on three-dimensional deflagration models, are presented. We mainly focus on one model,"c3_3d_256_10s", for which the hydrodynamics (Roepke 2005) and nucleosynthesis (Travaglio et al. 2004) was calculated up to the homologous phase of the explosion. Other models with different ignition conditions and different resolution are also briefly discussed. The synthetic spectra are compared to observed late time spectra. We find that while the model spectra after 300 to 500 days show a good agreement with the observed Fe II-III features, they also show too strong O I and C I lines compared to the observed late time spectra. The oxygen and carbon emission originates from the low-velocity unburned material in the central regions of these models. To get agreement between the models and observations we find that only a small mass of unburned material may be left in the center after the explosion. This may be a problem for pure deflagration models, although improved initial conditions, as well as higher resolution decrease the discrepancy. The relative intensity from the different ionization stages of iron is sensitive to the density of the emitting iron-rich material. We find that clumping, with the presence of low density regions, is needed to reproduce the observed iron emission, especially in the range between 4000 and 6000 AA. Both temperature and ionization depend sensitively on density, abundances and radioactive content. This work therefore illustrates the importance of including the inhomogeneous nature of realistic three-dimensional explosion models. We briefly discuss the implications of the spectral modeling for the nature of the explosion.

We present an extensive set of optical spectroscopy of the nearby type Ia supernova, SN 2002er, with 24 epochs spanning -11 to +34 days. Its spectral evolution is fairly typical of a type Ia supernova although it suffers high extinction. Nevertheless, there are differences in the spectral evolution when compared to coeval spectra of other normal type Ia supernova with comparable decline-rate parameters. Modelling of the photospheric phase spectra using a homogeneous abundance distribution in the atmosphere provides a fair match to the observations, but only by pushing the adopted distance and risetime close to the observational limits. Future improvements here will require models with a more realtistic stratified abundance distribution. From simple modelling of a nebular spectrum obtained at +215d, we infer a 56Ni mass of 0.69Msun, consistent with that derived from the light curve.

The piece-wise parabolic method (PPM) is applied to simulations of forced isotropic turbulence with Mach numbers $\sim 0.1... 1$. The equation of state is dominated by the Fermi pressure of an electron-degenerate fluid. The dissipation in these simulations is of purely numerical origin. For the dimensionless mean rate of dissipation, we find values in agreement with known results from mostly incompressible turbulence simulations. The calculation of a Smagorinsky length corresponding to the rate of numerical dissipation supports the notion of the PPM supplying an implicit subgrid scale model. In the turbulence energy spectra of various flow realisations, we find the so-called bottleneck phenomenon, i.e., a flattening of the spectrum function near the wavenumber of maximal dissipation. The shape of the bottleneck peak in the compensated spectrum functions is comparable to what is found in turbulence simulations with hyperviscosity. Although the bottleneck effect reduces the range of nearly inertial length scales considerably, we are able to estimate the value of the Kolmogorov constant. For steady turbulence with a balance between energy injection and dissipation, it appears that $C\approx 1.7$. However, a smaller value is found in the case of transonic turbulence with a large fraction of compressive components in the driving force. Moreover, we discuss length scales related to the dissipation, in particular, an effective numerical length scale $\Delta_{\mathrm{eff}}$, which can be regarded as the characteristic smoothing length of the implicit filter associated with the PPM.

Evidence of high-velocity features such as those seen in the near-maximum spectra of some Type Ia Supernovae (eg SN 2000cx) has been searched for in the available SNIa spectra observed earlier than one week before B maximum. Recent observational efforts have doubled the number of SNeIa with very early spectra. Remarkably, all SNeIa with early data (7 in our RTN sample and 10 from other programmes) show signs of such features, to a greater or lesser degree, in CaII IR, and some also in SiII 6255A line. High-velocity features may be interpreted as abundance or density enhancements. Abundance enhancements would imply an outer region dominated by Si and Ca. Density enhancements may result from the sweeping up of circumstellar material by the highest velocity SN ejecta. In this scenario, the high incidence of HVFs suggests that a thick disc and/or a high-density companion wind surrounds the exploding white dwarf, as may be the case in Single Degenerate systems. Large-scale angular fluctuations in the radial density and abundance distribution may also be responsible: this could originate in the explosion, and would suggest a deflagration as the more likely explosion mechanism. CSM-interaction and surface fluctuations may coexist, possibly leaving different signatures on the spectrum. In some SNe the HVFs are narrowly confined in velocity, suggesting the ejection of blobs of burned material.

The deflagration mode of flame propagation in Type Ia supernova (SN Ia) models requires a correct description of the interaction of the flame with turbulent motions. It is well-known that turbulent combustion proceeds in different regimes. For most parts of the deflagration in SNe Ia the flamelet regime applies. This has been modeled in previous multi-dimensional simulations. However, at late stages of the explosion, the flame will propagate in the so-called distributed regime. We investigate the effects of this regime on SN Ia models in a first and simplified approach and show that the trend of effects seems capable of curing some problems of current pure deflagration models.

We present the current status of our construction of synthetic spectra for type Ia supernovae. These properly take into account the effects of NLTE and an adequate representation of line blocking and blanketing. The models are based on a sophisticated atomic database. We show that the synthetic spectrum reproduces the observed spectrum of 'normal' SN-Ia near maximum light from the UV to the near-IR. However, further improvements are necessary before truly quantitative analyses of observed SN-Ia spectra can be performed. In particular, the inner boundary condition has to be fundamentally modified. This is due to the dominance of electron scattering over true absorption processes coupled with the flat density structure in these objects

We present full-star simulations of Type Ia supernova explosions on the basis of the standard Chandrasekhar-mass deflagration model. Most simulations so far considered only one spatial octant and assumed mirror symmetry to the other octants. Two full-star models are evolved to homologous expansion and compared with previous single-octant simulations. Therefrom we analyze the effect of abolishing the artificial symmetry constraint on the evolution of the flame surface. It turns out that the development of asymmetries depends on the chosen initial flame configuration. Such asymmetries of the explosion process could possibly contribute to the observed polarization of some Type Ia supernova spectra.

Extensive light and colour curves for the Type Ia supernova SN 2002er are presented as part of the European Supernova Collaboration. We have collected UBVRI photometry from ten different telescopes covering the phases from 7 days before until 619 days after maximum light. Corrections for the different instrumental systems and the non-thermal spectrum of the supernova (S-corrections) have been applied. With the densely sampled light curves we can make detailed comparisons to other well-observed objects. SN 2002er most closely resembles SN 1996X after maximum, but clearly shows a different colour evolution before peak light and a stronger shoulder in V and R bands compared to other well-observed SNe Ia. In particular, the rise time appears to be longer than what is expected from rise-time vs.decline-rate relation. We use several methods to determine the reddening towards SN 2002er based on the colour evolution at near peak and at late phases. The uvoir (bolometric) light curve shows great similarity with SN 1996X, but also indications of a higher luminosity, longer rise time and a more pronounced shoulder 25 days past maximum. The interpretation of the light curves was done with two independent light curve codes. Both find that given the luminosity of SN 2002er the 56Ni mass exceeds 0.6 Msun with prefered values near 0.7 Msun. Uncertainties in the exact distance to SN 2002er are the most serious limitation of this measurement. The light curve modelling also indicates a high level of mixing of the nickel in the explosion of SN 2002er.

We present the results of nucleosynthesis calculations based on multidimensional (2D and 3D) hydrodynamical simulations of the thermonuclear burning phase in SNIa. The detailed nucleosynthetic yields of our explosion models are calculated by post-processing the ejecta, using passively advected tracer particles. The nuclear reaction network employed in computing the explosive nucleosynthesis contains 383 nuclear species. We analyzed two different choices of ignition conditions (centrally ignited, in which the spherical initial flame geometry is perturbated with toroidal rings, and bubbles, in which multi-point ignition conditions are simulated). We show that unburned C and O varies typically from ~40% to ~50% of the total ejected material.The main differences between all our models and standard 1D computations are, besides the higher mass fraction of unburned C and O, the C/O ratio (in our case is typically a factor of 2.5 higher than in 1D computations), and somewhat lower abundances of certain intermediate mass nuclei such as S, Cl, Ar, K, and Ca, and of 56Ni. Because explosive C and O burning may produce the iron-group elements and their isotopes in rather different proportions one can get different 56Ni-fractions (and thus supernova luminosities) without changing the kinetic energy of the explosion. Finally, we show that we need the high resolution multi-point ignition (bubbles) model to burn most of the material in the center (demonstrating that high resolution coupled with a large number of ignition spots is crucial to get rid of unburned material in a pure deflagration SNIa model).

We present numerical calculations of the parameters $C_{\nu}$, $C_{\epsilon}$ and $C_{\kappa}$ associated with the common closures for turbulence production, dissipation and diffusion. In the case of homogeneous and isotropic turbulence, these parameters are expected to be statistically scale-invariant within the inertial subrange. In order to scrutinise this conjecture, we utilised a generalisation of the Germano filtering formalism, which is applicable to compressible flows as well. The filtering of data obtained from three-dimensional direct numerical simulations of forced isotropic turbulence with Mach numbers in the range $\sim 0.1...1$ then yielded values of the closure parameters associated with different length scales. The results indicate that the closure parameters are nearly universal for subsonic or moderately transonic flows, although the resolution of $432^{3}$ grid cells in our simulations is not quite sufficient to clearly establish scale invariance. In addition, it was found that the customary assumption of a kinetic Prandtl number of about unity for the gradient-diffusion closure is flawed due to the misalignment between turbulent flux and the gradient of the turbulence energy. Nevertheless, sound correlation can be achieved if the flux magnitude rather than the flux vector is locally matched. This conclusion is particularly useful for the family of subgrid scale models based on the turbulence energy equation. Furthermore, the parameter of production $C_{\nu}$ was computed in the fashion of dynamical procedures. Thereby, superior agreement between modelled and explicitly evaluated turbulence stresses in comparison to the eddy-viscosity closure with constant $C_{\nu}$ was verified.

We present preliminary results of the observational campaign on SN 2002er carried out by the European Supernova Collaboration (ESC). For this supernova the optical coverage from epochs -11 days to +50 days is excellent, both in imaging and spectroscopy. The extinction deduced from the photometry is E(B-V)=0.35, while Deltam15=1.31. In general SN 2002er behaves like a "normal" type Ia supernova, very similar to SN 1996X.

Type Ia supernovae have recently received considerable attention because it appears that they can be used as "standard candles" to measure cosmic distances out to billions of light years away from us. Observations of type Ia supernovae seem to indicate that we are living in a universe that started to accelerate its expansion when it was about half its present age. These conclusions rest primarily on phenomenological models which, however, lack proper theoretical understanding, mainly because the explosion process, initiated by thermonuclear fusion of carbon and oxygen into heavier elements, is difficult to simulate even on supercomputers. Here, we investigate a new way of modeling turbulent thermonuclear deflagration fronts in white dwarfs undergoing a type Ia supernova explosion. Our approach is based on a level set method which treats the front as a mathematical discontinuity and allows for full coupling between the front geometry and the flow field. New results of the method applied to the problem of type Ia supernovae are obtained. It is shown that in 2-D with high spatial resolution and a physically motivated subgrid scale model for the nuclear flames numerically "converged" results can be obtained, but for most initial conditions the stars do not explode. In contrast, simulations in 3-D, do give the desired explosions and many of their properties, such as the explosion energies, lightcurves and nucleosynthesis products, are in very good agreement with observed type Ia supernovae.

The propagation of cellularly stabilized thermonuclear flames is investigated by means of numerical simulations. In Type Ia supernova explosions the corresponding burning regime establishes at scales below the Gibson length. The cellular flame stabilization - which is a result of an interplay between the Landau-Darrieus instability and a nonlinear stabilization mechanism - is studied for the case of propagation into quiescent fuel as well as interaction with vortical fuel flows. Our simulations indicate that in thermonuclear supernova explosions stable cellular flames develop around the Gibson scale and that deflagration-to-detonation transition is unlikely to be triggered from flame evolution effects here.

One of the major challenges for theoretical modeling of Type Ia supernova explosions is to explain the diversity of these events and the empirically established correlation between their peak luminosity and light curve shape. In the framework of the so-called Chandrasekhar mass models, the progenitor's carbon-to-oxygen ratio has been suggested to be a principal source of peak luminosity variations due to a variation in the production of radioactive $^{56}$Ni during the explosion. The underlying idea is that an enhanced carbon mass fraction should result in a more vigorous explosion since here the energy release from nuclear reactions is increased. It was suspected that this would produce a higher amount of $^{56}$Ni in the ejecta. In this letter we describe a mechanism resulting from an interplay between nucleosynthesis and turbulent flame evolution which counteracts such an effect. Based on three-dimensional simulations we argue that it is nearly balanced and only minor differences in the amount of synthesized $^{56}$Ni with varying carbon mass fraction in the progenitor can be expected. Therefore this progenitor parameter is unlikely to account for the observed variations in Type Ia supernova luminosity. We discuss possible effects on the calibration of cosmological measurements.

We investigate the interaction of thermonuclear flames in Type Ia supernova explosions with vortical flows by means of numerical simulations. In our study, we focus on small scales, where the flame propagation is no longer dominated by the turbulent cascade originating from large-scale effects. Here, the flame propagation proceeds in the cellular burning regime, resulting from a balance between the Landau-Darrieus instability and its nonlinear stabilization. The interaction of a cellularly stabilized flame front with a vortical fuel flow is explored applying a variety of fuel densities and strengths of the velocity fluctuations. We find that the vortical flow can break up the cellular flame structure if it is sufficiently strong. In this case the flame structure adapts to the imprinted flow field. The transition from the cellularly stabilized front to the flame structure dominated by vortices of the flow proceeds in a smooth way. The implications of the results of our simulations for Type Ia Supernova explosion models are discussed.

We present a numerical investigation of the cellular burning regime in Type Ia supernova explosions. This regime holds at small scales (i.e. below the Gibson scale), which are unresolved in large-scale Type Ia supernova simulations. The fundamental effects that dominate the flame evolution here are the Landau-Darrieus instability and its nonlinear stabilization, leading to a stabilization of the flame in a cellular shape. The flame propagation into quiescent fuel is investigated addressing the dependence of the simulation results on the specific parameters of the numerical setup. Furthermore, we investigate the flame stability at a range of fuel densities. This is directly connected to the questions of active turbulent combustion (a mechanism of flame destabilization and subsequent self-turbulization) and a deflagration-to-detonation transition of the flame. In our simulations we find no substantial destabilization of the flame when propagating into quiescent fuels of densities down to ~10^7 g/cm^3, corroborating fundamental assumptions of large-scale SN Ia explosion models. For these models, however, we suggest an increased lower cutoff for the flame propagation velocity to take the cellular burning regime into account.

We present optical/near-infrared photometry and spectra of the type Ia SN 2002bo spanning epochs from -13 days before maximum B-band light to +102 days after. The pre-maximum optical coverage is particularly complete. In some respects, SN 2002bo behaves as a typical "Branch normal" type Ia supernova (SN Ia) at optical and IR wavelengths. We find a B-band risetime of 17.9+-0.5 days, a Dm_{15}(B) of 1.13+-0.05, and a M_B=-19.41+-0.42. However, comparison with other type Ia supernovae having similar Delta m_{15}(B) values indicates that in other respects SN 2002bo is unusual. While the optical spectra of SN 2002bo are very similar to those of SN 1984A, lower velocities and a generally more structured appearance are found in SNe 1990N, 1994D and 1998bu. For supernovae having Dm_(15)(B) > 1.2, we confirm the variation of R(SiII) (Nugent et al. 1995) with Dm_(15)(B). However, for supernovae such as SN2002bo, with lower values of Dm_(15)(B) the relation breaks down. Moreover, the evolution of R(SiII) for SN 2002bo is strikingly different from that shown by other type Ia supernovae. The velocities of SN 2002bo and 1984A derived from SII 5640A, SiII 6355A and CaII H&K lines are either much higher and/or evolve differently from those seen in other normal SNe Ia events. We suggest that the unusually low temperature, the presence of high-velocity intermediate-mass elements and the low abundance of carbon at early times indicates that burning to Si penetrated to much higher layers than in more normal type Ia supernovae. This may be indicative of a delayed-detonation explosion.

We present a numerical model which allows us to investigate thermonuclear flames in Type Ia supernova explosions. The model is based on a finite-volume explicit hydrodynamics solver employing PPM. Using the level-set technique combined with in-cell reconstruction and flux-splitting schemes we are able to describe the flame in the discontinuity approximation. We apply our implementation to flame propagation in Chandrasekhar-mass Type Ia supernova models. In particular we concentrate on intermediate scales between the flame width and the Gibson-scale, where the burning front is subject to the Landau-Darrieus instability. We are able to reproduce the theoretical prediction on the growth rates of perturbations in the linear regime and observe the stabilization of the flame in a cellular shape. The increase of the mean burning velocity due to the enlarged flame surface is measured. Results of our simulation are in agreement with semianalytical studies.

We present the results of three-dimensional hydrodynamical simulations of the subsonic thermonuclear burning phase in type Ia supernovae. The burning front model contains no adjustable parameters so that variations of the explosion outcome can be linked directly to changes in the initial conditions. In particular, we investigate the influence of the initial flame geometry on the explosion energy and find that it appears to be weaker than in 2D. Most importantly, our models predict global properties such as the produced nickel masses and ejecta velocities within their observed ranges without any fine tuning.

We present a numerical model to study the behavior of thermonuclear flames in the discontinuity approximation. This model is applied to investigate the Landau-Darrieus instability under conditions found in Type Ia supernova explosions of Chandrasekhar mass white dwarfs. This is a first step to explore the flame microstructure in these events. The model reproduces Landau's linearized stability analysis in early stages of the flame evolution and the stabilization in a cellular flame structure in the nonlinear stage.

Type Ia supernovae have become an indispensable tool for studying the expansion history of the universe, yet our understanding of the explosion mechanism is still incomplete. We describe the variety of discussed scenarios, sketch the most relevant physics, and report recent advances in multidimensional simulations of Chandrasekhar mass white dwarf explosions.

Following up on earlier work on this topic (Reinecke et al. 1999, A&A 347, pp. 724 and 739), we present an improved set of numerical models for simulations of white dwarfs exploding as Type Ia supernovae (SNe Ia). Two-dimensional simulations were used to test the reliability and numerical robustness of these algorithms; the results indicate that integral quantities like the total energy release are insensitive to changes of the grid resolution (above a certain threshold), which was not the case for our former code. The models were further enhanced to allow fully three-dimensional simulations of SNe Ia. A direct comparison of a 2D and a 3D calculation with identical initial conditions shows that the explosion is considerably more energetic in three dimensions; this is most likely caused by the assumption of axisymmetry in 2D, which inhibits the growth of flame instabilities in the azimuthal direction and thereby decreases the flame surface.

In this article we discuss the first simulations of two- and three-dimensional Type Ia supernovae with an improved hydrodynamics code. After describing the various enhancements, the obtained results are compared to those of earlier code versions, observational data and the findings of other researchers in this field.

Using the results of recent work in shear instabilities in stratified fluids, we show that the resonant interaction between large-scale flows in the accreted H/He envelope of white dwarf stars and interfacial gravity waves can mix with the star's envelope with the white dwarf's surface material, leading to the enhancement of the envelope's C/O abundance to levels required by extant models for nova outbursts.

Type Ia supernovae (SN Ia) are generally believed to be the result of the thermonuclear disruption of Chandrasekhar-mass carbon-oxygen white dwarfs, mainly because such thermonuclear explosions can account for the right amount of nickel, which is needed to explain the light curves and the late-time spectra, and the abundances of intermediate-mass nuclei which dominate the spectra near maximum light. Because of their enormous brightness and apparent homogeneity SN Ia have become an important tool to measure cosmological parameters. In this article the present understanding of the physics of thermonuclear explosions is reviewed. In particular, we focus our attention on subsonic ("deflagration") fronts, i.e. we investigate fronts propagating by heat diffusion and convection rather than by compression. Models based upon this mode of nuclear burning have been applied very successfully to the SN Ia problem, and are able to reproduce many of their observed features remarkably well. However, the models also indicate that SN Ia may differ considerably from each other, which is of importance if they are to be used as standard candles.

We present results of three-dimensional (3D) simulations of the magnetohydrodynamic Kelvin-Helmholtz instability in a stratified shear layer. The magnetic field is taken to be uniform and parallel to the shear flow. We describe the evolution of the fluid flow and the magnetic field for a range of initial conditions. In particular, we investigate how the mixing rate of the fluid depends on the Richardson number and the magnetic field strength. It was found that the magnetic field can enhance as well as suppress mixing. Moreover, we have performed two-dimensional (2D) simulations and discuss some interesting differences between the 2D and 3D results.

In this paper we present the results of numerical simulations of the Kelvin-Helmholtz instability in a stratified shear layer. This shear instability is believed to be responsible for extra mixing in differentially rotating stellar interiors and is the prime candidate to explain the abundance anomalies observed in many rotating stars. All mixing prescriptions currently in use are based on phenomenological and heuristic estimates whose validity is often unclear. Using three-dimensional numerical simulations, we study the mixing efficiency as a function of the Richardson number and compare our results with some semi-analytical formalisms of mixing.

Because calibrated light curves of Type Ia supernovae have become a major tool to determine the local expansion rate of the Universe and also its geometrical structure, considerable attention has been given to models of these events over the past couple of years. There are good reasons to believe that perhaps most Type Ia supernovae are the explosions of white dwarfs that have approached the Chandrasekhar mass, M_ch ~ 1.39 M_sun, and are disrupted by thermonuclear fusion of carbon and oxygen. However, the mechanism whereby such accreting carbon-oxygen white dwarfs explode continues to be uncertain. Recent progress in modeling Type Ia supernovae as well as several of the still open questions are addressed in this review. Although the main emphasis will be on studies of the explosion mechanism itself and on the related physical processes, including the physics of turbulent nuclear combustion in degenerate stars, we also discuss observational constraints.

Type Ia supernovae (SN Ia) are generally believed to be the result of the thermonuclear disruption of Chandrasekhar-mass carbon-oxygen white dwarfs, mainly because such thermonuclear explosions can account for the right amount of Ni-56, which is needed to explain the light curves and the late-time spectra, and the abundances of intermediate-mass nuclei which dominate the spectra near maximum light. Because of their enormous brightness and apparent homogeneity SN Ia have become an important tool to measure cosmological parameters. In this article the present understanding of the physics of thermonuclear explosions is reviewed. In particular, we focus our attention on subsonic (``deflagration'') fronts, i.e. we investigate fronts propagating by heat diffusion and convection rather than by compression. Models based upon this mode of nuclear burning have been applied very successfully to the SN Ia problem, and are able to reproduce many of their observed features remarkably well. However, the models also indicate that SN Ia may differ considerably from each other, which is of importance if they are to be used as standard candles.

We present a new way of modeling turbulent thermonuclear deflagration fronts in Chandrasekhar-mass white dwarfs, consisting of carbon and oxygen, undergoing a type Ia supernova explosion. Our approach is a front capturing/tracking hybrid scheme, based on a level set method, which treats the front as a mathematical discontinuity and allows for full coupling between the front geometry and the flow field. First results of the method applied to the problem of type Ia supernovae are discussed. It will be shown that even in 2-D and even with a physically motivated sub-grid model numerically ``converged'' results are difficult to obtain.

We present a statistical model which shows the influence of turbulence on a thermonuclear flame propagating in C+O white dwarf matter. Based on a Monte Carlo description of turbulence, it provides a method for investigating the physics in the so-called distributed burning regime. Using this method we perform numerical simulations of turbulent flames and show that in this particular regime the flamelet model for the turbulent flame velocity loses its validity. In fact, at high turbulent intensities burning in the distributed regime can lead to a deceleration of the turbulent flame and thus induces a competing process to turbulent effects that cause a higher flame speed. It is also shown that in dense C+O matter turbulent heat transport is described adequately by the Peclet number, Pe, rather than by the Reynolds number, which means that flame propagation is decoupled from small-scale turbulence. Finally, at the onset of our results we argue that the available turbulent energy in an exploding C+O white dwarf is probably too low in order to make a deflagration to detonation transition possible.

We investigate the possibility of a delayed detonation in a type Ia supernova under the assumption that the transition to detonation is triggered by turbulence only. Our discussion is based on the Zeldovich mechanism and suggests that typical turbulent velocities present during the explosion are not strong enough to allow this transition to occur. Although we are able to show that in carbon-rich matter (e.g., $X(^{12}$C$) = 0.75$) the possibility of a deflagration to detonation transition (DDT) is enhanced, even in this case the turbulent velocities needed are larger than the expected value of $u'(L) \approx 10^7 {cm s}^{-1}$ on a length-scale of $L \approx 10^6$ cm. Thus we conclude that a DDT may not be a common event during a thermonuclear explosion of a Chandrasekhar-mass white dwarf.

Radioactive isotopes such as 44Ti offer probably the most direct probe into nucleosynthesis environments in the interior of exploding stars, when the associated gamma-ray activities in the explosion remnant are detected and translated back to the isotopic abundances at the time of the explosion. In this paper, we assert that the procedure may not necessarily be straightforward, at least in the case of 44Ti, an orbital-electron capture decay isotope. Using the analytic model of McKee and Truelove (1995) for young supernova remnants, and assuming the existence of overdense 56Fe-dominated clumps that contain also 44Ti, we show that a high degree of ionization may be caused by the reverse shock so that the electron-capture rate of 44Ti could be significantly reduced from its laboratory value. When applied to Cas A, this increases under certain conditions the current 44Ti-activity by a factor 1.5 - 2.5, which yields a better compatibility between the COMPTEL observation of the 1.16 MeV line activity associated with the 44Ti decay and the SN model predictions of the initial 44Ti abundance. This possibility is, however, subject to various uncertainties, and in particular to the unknown properties and radial distribution of the the clumps in the ejecta.

First results of simulations are presented which compute the dynamical evolution of a Chandrasekhar-mass white dwarf, consisting of equal amounts of carbon and oxygen, from the onset of violent thermonuclear burning, by means of a new two-dimensional numerical code. Since in the interior of such a massive white dwarf nuclear burning progresses on microscopic scales as a sharp discontinuity, a so-called flamelet, which cannot be resolved by any numerical scheme, and since on macroscopic scales the burning front propagates due to turbulence, we make an attempt to model both effects explicitly in the framework of a finite-volume hydrodynamics code. Turbulence is included by a sub-grid model, following the spirit of large eddy simulations, and the well-localized burning front is treated by means of a level set, which allows us to compute the geometrical structure of the front more accurately than with previous methods. The only free parameters of our simulations are the location and the amount of nuclear fuel that is ignited as an initial perturbation. We find that models in which explosive carbon burning is ignited at the center remain bound by the time the front reaches low densities, where we stopped the computations because our description of combustion is no longer applicable. In contrast, off-center ignition models give rise to explosions which, however, are still too weak for typical Type Ia supernovae. Possible reasons for this rather disappointing result are discussed.

We present a new way of modeling deflagration fronts in reactive fluids, the main emphasis being on turbulent thermonuclear deflagration fronts in white dwarfs undergoing a Type Ia supernova explosion. Our approach is based on a level set method which treats the front as a mathematical discontinuity and allows full coupling between the front geometry and the flow field. With only minor modifications, this method can also be applied to describe contact discontinuities. Two different implementations are described and their physically correct behaviour for simple testcases is shown. First results of the method applied to the concrete problems of Type Ia supernovae and chemical hydrogen combustion are briefly discussed; a more extensive analysis of our astrophysical simulations is given in (Reinecke et al. 1998, MPA Green Report 1122b).

We present first results of three-dimensional (3D-) calculations of turbulent and degenerate hydrogen-burning on top of a C+O white dwarf of one solar mass. The simulations are carried out by means of a code which solves Euler's equation for an arbitrary equation of state together with a nuclear reaction network and the energy input from nuclear reactions on a Cartesian grid covering a fraction of the white dwarf's surface and accreted atmosphere. The flow patterns we obtain are very different from those of earlier 2D simulations using the same initial conditions and the same numerical resolution. The possibility of self-enrichment of the accreted hydrogen-rich atmosphere with carbon and oxygen from the surface layers of the white dwarf during the violent phase of the burning is investigated, and it is demonstrated that self-enrichment proceeds too slowly if the accreted gas has near-solar CNO-abundances at the onset of the thermonuclear runaway. As a result, we do not find a fast nova outburst. This conclusion remains valid if the initial metallicity of the accreted gas is raised by a factor of five. Therefore we conclude that fast nova outbursts indeed require huge enrichments of C and O, as postulated from spherically symmetric models, and that the mechanism which leads to such enhancements must operate prior to the outburst.

Collection of "Extended Abstracts" of talks presented at the SFB-375 workshop NEUTRINO ASTROPHYSICS, Ringberg Castle, Tegernsee, Germany, 20-24 Oct. 1997. Topics include Solar Neutrinos, Supernovae, Gamma-Ray Bursts, High-Energy Neutrinos (atmospheric, astrophysical), Cosmology, and Future Prospects. The focus is on astrophysical and experimental/observational aspects of astrophysical neutrinos while particle theory and neutrino laboratory experiments are not covered.

We present the results of two-dimensional calculations of turbulent nuclear burning of hydrogen-rich material accreted onto a white dwarf of 1 solar mass. The main aim of the present paper is to investigate the question as to whether and how the general properties of the burning are affected by numerical resolution effects. In particular, we want to see whether or not convective overshooting into the surface layers of the C+O white dwarf can lead to self-enrichment of the initially solar composition of the hydrogen-rich envelope with carbon and oxygen from the underlying white dwarf core. Our explicit hydrodynamic code is based on the PPM-method and computes the onset of the thermonuclear runaway on a Cartesian grid. In contrast to previous works we do not observe fast mixing of carbon and oxygen from the white dwarf's surface into the envelope by violent overshooting of large eddies. The main features of the flow fields in our simulations are the appearance of small persistent coherent structures of very high vorticity (and velocity) compared to the background flow. Their typical linear scales are about 10 to 20 grid zones and thus their physical size depends on the numerical resolution, i.e, their size decreases with increasing resolution. The two simulations (low and high resolution) which are presented here show only moderate differences in spatially integrated quantities such as laterally averaged temperature, energy generation rate, and chemical composition. We have not expanded both simulations equally long, but for the physical time under consideration the major difference seems to be that the highly resolved simulation is a bit less violent. In conclusion, we do find some self-enrichment, but on time-scales much longer than in previous calculations.

Several circumstantial arguments point to the formation of the third r-process peak at A about 190, near platinum, in stars of mass of about 8-10 solar masses: 1) The delayed production of europium with respect to iron imposes a time scale that restricts the progenitor stars to less than about 10 solar masses; 2) the r-process demands a dominant robust mechanism at least for barium and above, since the relative abundance pattern of those r-process elements in low-metallicity stars is consistent with the solar pattern; 3) stars of about 8-10 solar masses produce nearly identical degenerate O/Ne/Mg cores that collapse due to electron capture; and 4) the resulting low-mass cores may produce both an r-process in a prompt explosion and a subsequent r-process in a neutrino driven wind. The prompt explosion of an O/Ne/Mg core yields low entropy and low electron fraction, and hence may produce a reasonable r-process peak at A about 190 as well as all of the r-process elements with Z greater than 56. The possible differences in the neutrino-driven wind and associated r-process due to the low-mass neutron stars expected in this mass range are also discussed.

First results of numerical simulations are presented which compute the dynamical evolution of a neutron star with a mass slightly below the minimum stable mass by means of a new implicit (general relativistic) hydrodynamic code. We show that such a star first undergoes a phase of quasi-static expansion, caused by slow nuclear $\beta$-decays, lasting for about 20 seconds, but then explodes violently. The kinetic energy of the explosion is around $10^{49} erg$, the peak luminosity in electron anti-neutrinos is of order $10^{52} erg/s$, and the thermodynamic conditions of the expanding matter are favorable for r-process nucleosynthesis. These results are obtained for the Harrison-Wheeler equation of state and a simple and, possibly, unrealistic treatment of $\beta$-decay rates and nuclear fission, which were adopted for comparison with previous works. However, we do not expect that the outcome will change qualitatively if more recent nuclear input physics used. Although our study does not rely on a specific scenario ofhow a neutron star starting from a bigger (and stable) masscan reach the dynamical phase, we assume that the final mass-loss event happens on a very short time scale, i.e., on a time scale shorter than a sound-crossing time, by removing a certain amount of mass as an initial perturbation. This assumption implies that the star has no time to adjust its nuclear composition to the new mass mass through a sequence of quasi-equilibria. In the latter case, however, there exists no stable configuration below the minimum mass, because the equation of state of fully catalyzed matter is too soft. Therefore, the dynamics of the explosion will not be too different from what have obtained if different initial perturbations are assumed.

A series of two dimensional numerical simulations of explosive nuclear burning is presented for white dwarfs near the Chandraskhar mass. We assume that the burning begins as a slow deflagration front at or near the center of the star, and continues until the density in the burning regions has declined to about 10$^7$ g cm$^{-3}$, where the flame is essentially extinguished. We employ a novel numerical representation of the turbulent flame brush based upon ideas previously developed for modelling laboratory combustion and explore in some detail the sensitivity of the outcome to the manner in which burning is initiated. In particular, we simulate 1) a centrally ignited deflagration, 2) off-center ignition at a single ``point", and 3) simultaneous off-center ignition at five ``points". We find that the amount of $^{56}$Ni that is produced and other observable properties depend sensitively upon how the fuel is ignited.