Early observations of Type Ia supernovae (SNe$\,$Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. We present exquisite observations of SN$\,$2019yvq, the second observed SN$\,$Ia, after iPTF$\,$14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN$\,$2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for an SN$\,$Ia ($M_g \approx -18.5\,$mag at peak) yet featured very high absorption velocities ($v \approx 15,000\,\mathrm{km\,s}^{-1}$ for Si II $\lambda$6355 at peak). We find that many of the observational features of SN$\,$2019yvq, aside from the flash, can be explained if the explosive yield of radioactive $^{56}\mathrm{Ni}$ is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN$\,$2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of $^{56}\mathrm{Ni}$ in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN$\,$2019yvq. In closing, we predict that the nebular spectra of SN$\,$2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger.

Growing evidence suggests that Type Iax supernovae might be the result of thermonuclear deflagrations of Chandrasekhar-mass white dwarfs in binary systems. We carry out Monte Carlo radiative transfer simulations and predict spectropolarimetric features originating from the supernova explosion and subsequent ejecta interaction with the companion star. Specifically, we calculate viewing-angle dependent flux and polarisation spectra for a 3D model simulating the deflagration of a Chandrasekhar-mass white dwarf and, for a second model, simulating the ejecta interaction with a main-sequence star. We find that the intrinsic signal is weakly polarised and only mildly viewing-angle dependent, owing to the overall spherical symmetry of the explosion and the depolarising contribution of iron-group elements dominating the ejecta composition. The interaction with the companion star carves out a cavity in the ejecta and produces a detectable, but modest signal that is significant only at relatively blue wavelengths ($\lesssim$ 5000 $\unicode{x212B}$). In particular, increasingly fainter and redder spectra are predicted for observer orientations further from the cavity, while a modest polarisation signal $P\sim0.2$ per cent is found at blue wavelengths for orientations 30$^\circ$ and 45$^\circ$ away from the cavity. We find a reasonable agreement between the interaction model viewed from these orientations and spectropolarimetric data of SN 2005hk and interpret the maximum-light polarisation signal seen at blue wavelengths for this event as a possible signature of the ejecta-companion interaction. We encourage further polarimetric observations of SNe Iax to test whether our results can be extended and generalised to the whole SN Iax class.

We present the HOLISMOKES project on strong gravitational lensing of supernovae (SNe) as a probe of SN physics and cosmology. We investigate the effects of microlensing on early-phase SN Ia spectra, and find that within 10 rest-frame days after SN explosion, distortions of SN Ia spectra due to microlensing are typically negligible (<1% distortion within the 68% credible region, and ~10% distortion within the 95% credible region). This shows great prospects of using lensed SNe Ia to obtain intrinsic early-phase SN spectra for deciphering SN Ia progenitors. As a demonstration of the usefulness of lensed SNe Ia for cosmology, we simulate a sample of mock lensed SN Ia systems that are expected to have accurate and precise time-delay measurements in the era of the Legacy Survey of Space and Time (LSST). Adopting realistic yet conservative uncertainties on their time-delay distances and lens angular diameter distances (of 6.6% and 5%, respectively), we find that a sample of 20 lensed SNe Ia would allow a constraint on the Hubble constant ($H_0$) with 1.3% uncertainty in the flat $\Lambda$CDM cosmology. We find a similar constraint on $H_0$ in an open $\Lambda$CDM cosmology, while the constraint degrades to 3% in a flat wCDM cosmology. We anticipate lensed SNe to be an independent and powerful probe of SN physics and cosmology in the upcoming LSST era.

Sub-Chandrasekhar mass white dwarfs accreting a helium shell on a carbon-oxygen core are potential progenitors of normal Type Ia supernovae. This work focuses on the details of the onset of the carbon detonation in the double detonation sub-Chandrasekhar model. In order to simulate the influence of core-shell mixing on the carbon ignition mechanism, the helium shell and its detonation are followed with an increased resolution compared to the rest of the star treating the propagation of the detonation wave more accurately. This significantly improves the predictions of the nucleosynthetic yields from the helium burning. The simulations were carried out with the AREPO code. A carbon-oxygen core with a helium shell was set up in one dimension and mapped to three dimensions. We ensured the stability of the white dwarf with a relaxation step before the hydrodynamic detonation simulation started. Synthetic observables were calculated with the radiative transfer code ARTIS. An ignition mechanism of the carbon detonation was observed, which received little attention before. In this "scissors mechanism", the impact the helium detonation wave has on unburnt material when converging opposite to its ignition spot is strong enough to ignite a carbon detonation. This is possible in a carbon enriched transition region between the core and shell. The detonation mechanism is found to be sensitive to details of the core-shell transition and our models illustrate the need to consider core-shell mixing taking place during the accretion process. Even though the detonation ignition mechanism differs form the converging shock mechanism, the differences in the synthetic observables are not significant. Though they do not fit observations better than previous simulations, they illustrate the need for multi-dimensional simulations.

We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium (non-LTE) population and ionisation solver, a new multi-frequency radiation field model, and a new atomic dataset with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionisation, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf in which the possible effects of gravitational settling of Ne22 prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M$_\odot$ white dwarf model obtained when 5.5 Gyr of very-efficient settling is assumed to a similar model without settling. We find that this degree of Ne22 settling has only a modest effect on the resulting nebular spectra due to increased Ni58 abundance. Due to the high ionisation in sub-Chandrasekhar models, the nebular [Ni II] emission remains negligible, while the [Ni III] line strengths are increased and the overall ionisation balance is slightly lowered in the model with Ne22 settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe III] emission relative to [Fe II] in comparison to observations of normal Type Ia supernovae.

The upcoming Large Synoptic Survey Telescope (LSST) will detect many strongly lensed Type Ia supernovae (LSNe Ia) for time-delay cosmography. This will provide an independent and direct way for measuring the Hubble constant $H_0$, which is necessary to address the current $4.4 \sigma$ tension in $H_0$ between the local distance ladder and the early Universe measurements. We present a detailed analysis of different observing strategies for the LSST, and quantify their impact on time-delay measurement between multiple images of LSNe Ia. For this, we produced microlensed mock-LSST light curves for which we estimated the time delay between different images. We find that using only LSST data for time-delay cosmography is not ideal. Instead, we advocate using LSST as a discovery machine for LSNe Ia, enabling time delay measurements from follow-up observations from other instruments in order to increase the number of systems by a factor of 2 to 16 depending on the observing strategy. Furthermore, we find that LSST observing strategies, which provide a good sampling frequency (the mean inter-night gap is around two days) and high cumulative season length (ten seasons with a season length of around 170 days per season), are favored. Rolling cadences subdivide the survey and focus on different parts in different years; these observing strategies trade the number of seasons for better sampling frequency. In our investigation, this leads to half the number of systems in comparison to the best observing strategy. Therefore rolling cadences are disfavored because the gain from the increased sampling frequency cannot compensate for the shortened cumulative season length. We anticipate that the sample of lensed SNe Ia from our preferred LSST cadence strategies with rapid follow-up observations would yield an independent percent-level constraint on $H_0$.

We present a comprehensive dataset of optical and near-infrared photometry and spectroscopy of type~Ia supernova (SN) 2016hnk, combined with integral field spectroscopy (IFS) of its host galaxy, MCG -01-06-070, and nearby environment. Properties of the SN local environment are characterized by means of single stellar population synthesis applied to IFS observations taken two years after the SN exploded. SN 2016hnk spectra are compared to other 1991bg-like SNe Ia, 2002es-like SNe Ia, and Ca-rich transients. In addition, abundance stratification modelling is used to identify the various spectral features in the early phase spectral sequence and the dataset is also compared to a modified non-LTE model previously produced for the sublumnious SN 1999by. SN 2016hnk is consistent with being a sub-luminous (M$_{\rm B}=-16.7$ mag, s$_{\rm BV}$=0.43$\pm$0.03), highly reddened object. IFS of its host galaxy reveals both a significant amount of dust at the SN location, as well as residual star formation and a high proportion of old stellar populations in the local environment compared to other locations in the galaxy, which favours an old progenitor for SN 2016hnk. Inspection of a nebular spectrum obtained one year after maximum contains two narrow emission lines attributed to the forbidden [Ca II] $\lambda\lambda$7291,7324 doublet with a Doppler shift of 700 km s$^{-1}$. Based on various observational diagnostics, we argue that the progenitor of SN 2016hnk was likely a near Chandrasekhar-mass ($M_{\rm Ch}$) carbon-oxygen white dwarf that produced 0.108 $M_\odot$ of $^{56}$Ni. Our modeling suggests that the narrow [Ca II] features observed in the nebular spectrum are associated with $^{48}$Ca from electron capture during the explosion, which is expected to occur only in white dwarfs that explode near or at the $M_{\rm Ch}$ limit.

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.

We present JEKYLL, a new code for modelling of supernova (SN) spectra and lightcurves based on Monte-Carlo (MC) techniques for the radiative transfer. The code assumes spherical symmetry, homologous expansion and steady state for the matter, but is otherwise capable of solving the time-dependent radiative transfer problem in non-local-thermodynamic-equilibrium (NLTE). The method used was introduced in a series of papers by Lucy, but the full time-dependent NLTE capabilities of it have never been tested. Here, we have extended the method to include non-thermal excitation and ionization as well as charge-transfer and two-photon processes. Based on earlier work, the non-thermal rates are calculated by solving the Spencer-Fano equation. Using a method previously developed for the SUMO code, macroscopic mixing of the material is taken into account in a statistical sense. In addition, a statistical Markov-chain model is used to sample the emission frequency, and we introduce a method to control the sampling of the radiation field. Except for a description of JEKYLL, we provide comparisons with the ARTIS, SUMO and CMFGEN codes, which show good agreement in the calculated spectra as well as the state of the gas. In particular, the comparison with CMFGEN, which is similar in terms of physics but uses a different technique, shows that the Lucy method does indeed converge in the time-dependent NLTE case. Finally, as an example of the time-dependent NLTE capabilities of JEKYLL, we present a model of a Type IIb SN, taken from a set of models presented and discussed in detail in an accompanying paper. Based on this model we investigate the effects of NLTE, in particular those arising from non-thermal excitation and ionization, and find strong effects even on the bolometric lightcurve. This highlights the need for full NLTE calculations when simulating the spectra and lightcurves of SNe.

We present direct spectroscopic modeling of five Type Iax supernovae (SNe) with the one dimensional Monte Carlo radiative transfer code TARDIS. The abundance tomography technique is used to map the chemical structure and physical properties of the SN atmosphere. Through via fitting of multiple spectral epochs with self-consistent ejecta models, we can then constrain the location of some elements within the ejecta. The synthetic spectra of the best-fit models are able to reproduce the flux continuum and the main absorption features in the whole sample. We find that the mass fractions of IGEs and IMEs show a decreasing trend toward the outer regions of the atmospheres using density profiles similar to those of deflagration models in the literature. Oxygen is the only element, which could be dominant at higher velocities. The stratified abundance structure contradicts the well-mixed chemical profiles predicted by pure deflagration models. Based on the derived densities and abundances, a template model atmosphere is created for the SN Iax class and compared to the observed spectra. Free parameters are the scaling of the density profile, the velocity shift of the abundance template, and the peak luminosity. The results of this test support the idea that all SNe Iax can be described by a similar internal structure, which argues for a common origin of this class of explosions.

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.

Double detonations in double white dwarf (WD) binaries undergoing unstable mass transfer have emerged in recent years as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. One potential outcome of this "dynamically driven double-degenerate double-detonation" (D^6) scenario is that the companion WD survives the explosion and is flung away with a velocity equal to its > 1000 km/s pre-SN orbital velocity. We perform a search for these hypervelocity runaway WDs using Gaia's second data release. In this paper, we discuss seven candidates followed up with ground-based instruments. Three sources are likely to be some of the fastest known stars in the Milky Way, with total Galactocentric velocities between 1000 and 3000 km/s, and are consistent with having previously been companion WDs in pre-SN Ia systems. However, although the radial velocity of one of the stars is > 1000 km/s, the radial velocities of the other two stars are puzzlingly consistent with 0. The combined five-parameter astrometric solutions from Gaia and radial velocities from follow-up spectra yield tentative 6D confirmation of the D^6 scenario. The past position of one of these stars places it within a faint, old SN remnant, further strengthening the interpretation of these candidates as hypervelocity runaways from binary systems that underwent SNe Ia.

The recently discovered bright type II supernova OGLE14-073 evolved very slowly. The light curve rose to maximum for 90 days from discovery and then declined at a rate compatible with the radioactive decay of 56Co. In this study, we show that a pair-instability supernova is a plausible mechanism for this event. We calculate explosion models and light curves with the radiation hydrodynamics code STELLA starting from two M(ZAMS) = 150 solar masses, Z=0.001 progenitors. We obtain satisfactory fits to OGLE14-073 broadband light curves by including additional 56Ni in the centre of the models and mixing hydrogen down into the inner layers of the ejecta to a radial mass coordinate of 10 solar masses. The extra 56Ni required points to a slightly more massive progenitor star. The mixing of hydrogen could be due to large scale mixing during the explosion. We also present synthetic spectra for our models simulated with the Monte Carlo radiative transfer code ARTIS. The synthetic spectra reproduce the main features of the observed spectra of OGLE14-073. We conclude that OGLE14-073 is one of the most promising candidates for a pair-instability explosion.

We present 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.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

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.

We present a detailed spectral analysis of Type Iax SN 2011ay. Our spectra cover epochs between -3 and +19 days with respect to the maximum light in B-band. This time range allows us to employ a so-called abundance tomography technique. The synthetic spectral fitting was made with the 1D Monte Carlo radiative transfer code TARDIS. In this paper, we describe our method to fit multiple epochs with a self-consistent, stratified atmospheric model. We compare our results to previously published SYN++ models and the predictions of different explosion scenarios. Using a fixed density profile (exponential fit of W7), we find that a uniform abundance model cannot reproduce the spectral features before maximum light because of the emergence of excessively strong Fe lines. In our best-fit TARDIS model, we find an abundance profile that separated into two different regimes: a well-mixed region under 10,000 km s$^{-1}$ and a stratified region with decreasing IGE abundances above 10,000 km s$^{-1}$. Based on a detailed comparative analysis, our conclusion is that the available pure deflagration models cannot fully explain either the observed properties of SN 2011ay or the results of our TARDIS modeling. Further examinations are necessary to find an adequate explanation for the origin of this object.

Despite the importance of Type Ia supernovae (SNe Ia) for modern astrophysics, their detailed mechanism is still not fully understood. In this contribution, we present recent findings from numerical explosion models in the context of the observed diversity of SNe Ia and we discuss how these models can help to shed light on the explosion mechanism and the progenitor stars of SNe Ia. In addition, we introduce the Heidelberg Supernova Model Archive (HESMA), a new online data base where we provide integrated isotopic abundances and radially averaged ejecta profiles and synthetic observables for a wide range of state-of-the-art explosion models.

The light from distant supernovae (SNe) can be magnified through gravitational lensing when a foreground galaxy is located along the line of sight. This line-up allows for detailed studies of SNe at high redshift that otherwise would not be possible. Spectroscopic observations of lensed high-redshift Type Ia supernovae (SNe Ia) are of particular interest since they can be used to test for evolution of their intrinsic properties. The use of SNe Ia for probing the cosmic expansion history has proven to be an extremely powerful method for measuring cosmological parameters. However, if systematic redshift-dependent properties are found, their usefulness for future surveys could be challenged. We investigate whether the spectroscopic properties of the strongly lensed and very distant SN Ia PS1-10afx at $z=1.4$ deviates from the well-studied populations of normal SNe Ia at nearby or intermediate distance. We created median spectra from nearby and intermediate-redshift spectroscopically normal SNe Ia from the literature at -5 and +1 days from light-curve maximum. We then compared these median spectra to those of PS1-10afx. We do not find signs of spectral evolution in PS1-10afx. The observed deviation between PS1-10afx and the median templates are within what is found for SNe at low- and intermediate-redshift. There is a noticeable broad feature centred at $\rm \lambda\sim 3500$~\AA{}, which is present only to a lesser extent in individual low and intermediate redshift SN Ia spectra. From a comparison with a recently developed explosion model, we find this feature to be dominated by iron peak elements, in particular, singly ionized cobalt and chromium.

SN 2014J in M82 is the closest Type Ia supernova (SN Ia) in decades. The proximity allows for detailed studies of supernova physics and provides insights into the circumstellar and interstellar environment. In this work we analyze Spitzer mid-IR data of SN 2014J in the 3.6 and 4.5 {\mu}m wavelength range, together with several other nearby and well-studied SNe Ia. We compile the first composite mid-IR light-curve templates from our sample of SNe~Ia, spanning the range from before peak brightness well into the nebular phase. Our observations indicate that SNe Ia form a very homogeneous class of objects at these wavelengths. Using the low-reddening supernovae for comparison, we constrain possible thermal emission from circumstellar dust around the highly reddened SN 2014J. We also study SNe 2006X and 2007le, where the presence of matter in the circumstellar environment has been suggested. No significant mid-IR excess is detected, allowing us to place upper limits on the amount of pre-existing dust in the circumstellar environment. For SN 2014J, $M_{dust} < 10^{-5}$ M$_{\odot}$ within $r_{dust} \sim 10^{17}$ cm, which is insufficient to account for the observed extinction. Similar limits are obtained for SNe 2006X and 2007le.

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 present a search for an electromagnetic counterpart of the gravitational wave source GW151226. Using the Pan-STARRS1 telescope we mapped out 290 square degrees in the optical i_ps filter starting 11.5hr after the LIGO information release and lasting for a further 28 days. The first observations started 49.5hr after the time of the GW151226 detection. We typically reached sensitivity limits of i_ps = 20.3-20.8 and covered 26.5% of the LIGO probability skymap. We supplemented this with ATLAS survey data, reaching 31% of the probability region to shallower depths of m~19. We found 49 extragalactic transients (that are not obviously AGN), including a faint transient in a galaxy at 7Mpc (a luminous blue variable outburst) plus a rapidly decaying M-dwarf flare. Spectral classification of 20 other transient events showed them all to be supernovae. We found an unusual transient, PS15dpn, with an explosion date temporally coincident with GW151226 which evolved into a type Ibn supernova. The redshift of the transient is secure at z=0.1747 +/- 0.0001 and we find it unlikely to be linked, since the luminosity distance has a negligible probability of being consistent with that of GW151226. In the 290 square degrees surveyed we therefore do not find a likely counterpart. However we show that our survey strategy would be sensitive to NS-NS mergers producing kilonovae at D < 100 Mpc which is promising for future LIGO/Virgo searches.

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.

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

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.

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 progenitor problem of Type Ia supernovae (SNe Ia) is still unsolved. Most of these events are thought to be explosions of carbon-oxygen (CO) white dwarfs (WDs), but for many of the explosion scenarios, particularly those involving the externally triggered detonation of a sub-Chandrasekhar mass WD (sub-M Ch WD), there is also a possibility of having an oxygen-neon (ONe) WD as progenitor. We simulate detonations of ONe WDs and calculate synthetic observables from these models. The results are compared with detonations in CO WDs of similar mass and observational data of SNe Ia. We perform hydrodynamic explosion simulations of detonations in initially hydrostatic ONe WDs for a range of masses below the Chandrasekhar mass (M Ch), followed by detailed nucleosynthetic postprocessing with a 384-isotope nuclear reaction network. The results are used to calculate synthetic spectra and light curves, which are then compared with observations of SNe Ia. We also perform binary evolution calculations to determine the number of SNe Ia involving ONe WDs relative to the number of other promising progenitor channels. The ejecta structures of our simulated detonations in sub-M Ch ONe WDs are similar to those from CO WDs. There are, however, small systematic deviations in the mass fractions and the ejecta velocities. These lead to spectral features that are systematically less blueshifted. Nevertheless, the synthetic observables of our ONe WD explosions are similar to those obtained from CO models. Our binary evolution calculations show that a significant fraction (3-10%) of potential progenitor systems should contain an ONe WD. The comparison of our ONe models with our CO models of comparable mass (1.2 Msun) shows that the less blueshifted spectral features fit the observations better, although they are too bright for normal SNe Ia.

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

Left-over, ablated material from a possible non-degenerate companion can reveal itself after about one year in spectra of Type Ia SNe (SNe Ia). We have searched for such material in spectra of SN 2011fe (at 294 days after the explosion) and for SN 2014J (315 days past explosion). The observations are compared with numerical models simulating the expected line emission. The spectral lines sought for are H-alpha, [O I] 6300 and [Ca II] 7291,7324, and the expected width of these lines is about 1000 km/s. No signs of these lines can be traced in any of the two supernovae. When systematic uncertainties are included, the limits on hydrogen-rich ablated gas in SNe 2011fe and 2014J are 0.003 M_sun and 0.0085 M_sun, respectively, where the limit for SN 2014J is the second lowest ever, and the limit for SN 2011fe is a revision of a previous limit. Limits are also put on helium-rich ablated gas. These limits are used, in conjunction with other data, to argue that these supernovae can stem from double-degenerate systems, or from single-degenerate systems with a spun up/spun down super-Chandrasekhar white dwarf. For SN 2011fe, other types of hydrogen-rich donors can likely be ruled out, whereas for SN 2014J a main-sequence donor system with large intrinsic separation is still possible. Helium-rich donor systems cannot be ruled out for any of the two supernovae, but the expected short delay time for such progenitors makes this possibility less likely, especially for SN 2011fe. The broad [Ni II] 7378 emission in SN 2014J is redshifted by about +1300 km/s, as opposed to the known blueshift of roughly -1100 km/s for SN 2011fe. [Fe II] 7155 is also redshifted in SN 2014J. SN 2014J belongs to a minority of SNe Ia that both have a nebular redshift of [Fe II] 7155 and [Ni II] 7378, and a slow decline of the Si II 6355 absorption trough just after B-band maximum.

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.

In this letter we present an optical spectrum of SN 2011fe taken 1034 d after the explosion, several hundred days later than any other spectrum of a Type Ia supernova (disregarding light-echo spectra and local-group remnants). The spectrum is still dominated by broad emission features, with no trace of a light echo or interaction of the supernova ejecta with surrounding interstellar material. Comparing this extremely late spectrum to an earlier one taken 331 d after the explosion, we find that the most prominent feature at 331 d - [Fe III] emission around 4700 A - has entirely faded away, suggesting a significant change in the ionisation state. Instead, [Fe II] lines are probably responsible for most of the emission at 1034 d. An emission feature at 6300-6400 A has newly developed at 1034 d, which we tentatively identify with Fe I {\lambda}6359, [Fe I] {\lambda}{\lambda}6231, 6394 or [O I] {\lambda}{\lambda}6300, 6364. Interestingly, the features in the 1034-d spectrum seem to be collectively redshifted, a phenomenon that we currently have no convincing explanation for. We discuss the implications of our findings for explosion models, but conclude that sophisticated spectral modelling is required for any firm statement.

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 study the very early lightcurve of supernova 2014J (SN 2014J) using the high-cadence broad-band imaging data obtained by the Kilodegree Extremely Little Telescope (KELT), which fortuitously observed M 82 around the time of the explosion, starting more than two months prior to detection, with up to 20 observations per night. These observations are complemented by observations in two narrow-band filters used in an H$\alpha$ survey of nearby galaxies by the intermediate Palomar Transient Factory (iPTF) that also captured the first days of the brightening of the \sn. The evolution of the lightcurves is consistent with the expected signal from the cooling of shock heated material of large scale dimensions, $\gsim 1 R_{\odot}$. This could be due to heated material of the progenitor, a companion star or pre-existing circumstellar environment, e.g., in the form of an accretion disk. Structure seen in the lightcurves during the first days after explosion could also originate from radioactive material in the outer parts of an exploding white dwarf, as suggested from the early detection of gamma-rays. The model degeneracy translates into a systematic uncertainty of $\pm 0.3$ days on the estimate of the first light from SN 2014J.

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.

Binary stellar evolution calculations predict that Chandrasekhar-mass carbon/oxygen white dwarfs (WDs) show a radially varying profile for the composition with a carbon depleted core. Many recent multi-dimensional simulations of Type Ia supernovae (SNe Ia), however, assume the progenitor WD has a homogeneous chemical composition. In this work, we explore the impact of different initial carbon profiles of the progenitor WD on the explosion phase and on synthetic observables in the Chandrasekhar-mass delayed detonation model. Spectra and light curves are compared to observations to judge the validity of the model. The explosion phase is simulated using the finite volume supernova code LEAFS, which is extended to treat different compositions of the progenitor WD. The synthetic observables are computed with the Monte Carlo radiative transfer code ARTIS. Differences in binding energies of carbon and oxygen lead to a lower nuclear energy release for carbon depleted material; thus, the burning fronts that develop are weaker and the total nuclear energy release is smaller. For otherwise identical conditions, carbon depleted models produce less Ni-56. Comparing different models with similar Ni-56 yields shows lower kinetic energies in the ejecta for carbon depleted models, but only small differences in velocity distributions and line velocities in spectra. The light curve width-luminosity relation (WLR) obtained for models with differing carbon depletion is roughly perpendicular to the observed WLR, hence the carbon mass fraction is probably only a secondary parameter in the family of SNe Ia.

Type-Ia supernovae result from binary systems that include a carbon-oxygen white dwarf, and these thermonuclear explosions typically produce 0.5 M_solar of radioactive 56Ni. The 56Ni is commonly believed to be buried deeply in the expanding supernova cloud. Surprisingly, in SN2014J we detected the lines at 158 and 812 keV from 56Ni decay ({\tau}~8.8 days) earlier than the expected several-week time scale, only ~20 days after the explosion, and with flux levels corresponding to roughly 10% of the total expected amount of 56Ni. Some mechanism must break the spherical symmetry of the supernova, and at the same time create a major amount of 56Ni at the outskirts. A plausible explanation is that a belt of helium from the companion star is accreted by the white dwarf, where this material explodes and then triggers the supernova event.

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.

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

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.

The most favored progenitor scenarios for Type Ia supernovae (SNe Ia) involve the single-degenerate (SD) scenario and the double-degenerate scenario. The absence of stripped hydrogen (H) in the nebular spectra of SNe Ia challenges the SD progenitor models. Recently, it was shown that pure deflagration explosion models of Chandrasekhar-mass white dwarfs ignited off-center reproduce the characteristic observational features of 2002cx-like SNe Ia very well. In this work we predict, for the first time, the amount of stripped H for the off-center pure deflagration explosions. We find that their low kinetic energies lead to inefficient H mass stripping (less than 0.01 M_sun), indicating that the stripped H may be hidden in (observed) late-time spectra of SN 2002cx-like SNe Ia.

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.

The progenitors of Type Ia Supernovae (SNe Ia) are still unknown, despite significant progress during the last years in theory and observations. Violent mergers of two carbon--oxygen (CO) white dwarfs (WDs) are one candidate suggested to be responsible for at least a significant fraction of normal SNe Ia. Here, we simulate the merger of two CO WDs using a moving-mesh code that allows for the inclusion of thin helium (He) shells (0.01\,\msun) on top of the WDs, at an unprecedented numerical resolution. The accretion of He onto the primary WD leads to the formation of a detonation in its He shell. This detonation propagates around the CO WD and sends a converging shock wave into its core, known to robustly trigger a second detonation, as in the well-known double-detonation scenario for He-accreting CO WDs. However, in contrast to that scenario where a massive He shell is required to form a detonation through thermal instability, here the He detonation is ignited dynamically. Accordingly the required He-shell mass is significantly smaller, and hence its burning products are unlikely to affect the optical display of the explosion. We show that this scenario, which works for CO primary WDs with CO- as well as He-WD companions, has the potential to explain the different brightness distributions, delay times and relative rates of normal and fast declining SNe Ia. Finally, we discuss extensions to our unified merger model needed to obtain a comprehensive picture of the full observed diversity of SNe Ia.

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.

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 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.

We present ten medium-resolution, high signal-to-noise ratio near-infrared (NIR) spectra of SN 2011fe from SpeX on the NASA Infrared Telescope Facility (IRTF) and Gemini Near-Infrared Spectrograph (GNIRS) on Gemini North, obtained as part of the Carnegie Supernova Project. This data set constitutes the earliest time-series NIR spectroscopy of a Type Ia supernova (SN Ia), with the first spectrum obtained at 2.58 days past the explosion and covering -14.6 to +17.3 days relative to B-band maximum. C I {\lambda}1.0693 {\mu}m is detected in SN 2011fe with increasing strength up to maximum light. The delay in the onset of the NIR C I line demonstrates its potential to be an effective tracer of unprocessed material. For the first time in a SN Ia, the early rapid decline of the Mg II {\lambda}1.0927 {\mu}m velocity was observed, and the subsequent velocity is remarkably constant. The Mg II velocity during this constant phase locates the inner edge of carbon burning and probes the conditions under which the transition from deflagration to detonation occurs. We show that the Mg II velocity does not correlate with the optical light-curve decline rate {\Delta}m15. The prominent break at ~1.5 {\mu}m is the main source of concern for NIR k-correction calculations. We demonstrate here that the feature has a uniform time evolution among SNe Ia, with the flux ratio across the break strongly correlated with {\Delta}m15. The predictability of the strength and the onset of this feature suggests that the associated k-correction uncertainties can be minimized with improved spectral templates.

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.

We carry out the first multi-dimensional radiative transfer calculations to simultaneously compute synthetic spectra and light curves for models of supernovae driven by fast bipolar outflows. These allow us to make self-consistent predictions for the orientation dependence of both colour evolution and spectral features. We compare models with different degrees of asphericity and metallicity and find significant observable consequences of both. In aspherical models, we find spectral and light curve features that vary systematically with observer orientation. In particular, we find that the early phase light curves are brighter and bluer when viewed close to the polar axis but that the peak flux is highest for equatorial (off-axis) inclinations. Spectral line features also depend systematically on observer orientation, including the velocity of the SiII 6355A line. Consequently, our models predict a correlation between line velocity and color that could assist the identification of supernovae associated with off-axis jet-driven explosions. The amplitude and range of this correlation depends on the degree of asphericity, the metallicity and the epoch of observation but we find that it is always present and acts in the same direction.

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.

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.

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.

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.

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.

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.

A Monte Carlo code (ARTIS) for modelling time-dependent three-dimensional spectral synthesis in chemically inhomogeneous models of Type Ia supernova ejecta is presented. Following the propagation of gamma-ray photons, emitted by the radioactive decay of the nucleosynthesis products, energy is deposited in the supernova ejecta and the radiative transfer problem is solved self-consistently, enforcing the constraint of energy conservation in the co-moving frame. Assuming a photoionisation dominated plasma, the equations of ionisation equilibrium are solved together with the thermal balance equation adopting an approximate treatment of excitation. Since we implement a fully general treatment of line formation, there are no free parameters to adjust. Thus a direct comparison between synthetic spectra and light curves, calculated from hydrodynamic explosion models, and observations is feasible. The code is applied to the well known W7 explosion model and the results tested against other studies. Finally the effect of asymmetric ejecta on broad band light curves and spectra is illustrated using an elliptical toy model.

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.

Dwarf nova outbursts result from enhanced mass transport through the accretion disc of a cataclysmic variable system. We assess the question of whether these outbursts are caused by an enhanced mass transfer from the late-type main sequence star onto the white dwarf (so-called mass transfer instability model, MTI) or by a thermal instability in the accretion disc (disc instability model, DIM). We compute non-LTE models and spectra of accretion discs in quiescence and outburst and construct spectral time sequences for discs over a complete outburst cycle. We then compare our spectra to published optical spectroscopy of the dwarf nova SS Cygni. In particular, we investigate the hydrogen and helium line profiles that are turning from emission into absorption during the rise to outburst. The evolution of the hydrogen and helium line profiles during the rise to outburst and decline clearly favour the disc-instability model. Our spectral model sequences allow us to distinguish inside-out and outside-in moving heating waves in the disc of SS Cygni, which can be related to symmetric and asymmetric outburst light curves, respectively.