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

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

Particle simulation study of electron heating by counter-streaming ion beams ahead of supernova remnant shocks

The growth and saturation of Buneman-type instabilities is examined with a particle-in-cell (PIC) simulation for parameters that are representative for the foreshock region of fast supernova remnant shocks. A dense ion beam and the electrons correspond to the upstream plasma and a fast ion beam to the shock-reflected ions. The purpose of the 2D simulation is to identify the nonlinear saturation mechanisms, the electron heating and potential secondary instabilities that arise from anisotropic electron heating and result in the growth of magnetic fields. We confirm that the instabilities between both ion beams and the electrons saturate by the formation of phase space holes by the beam-aligned modes. The slower oblique modes accelerate some electrons, but they cannot heat up the electrons significantly before they are trapped by the faster beam-aligned modes. Two circular electron velocity distributions develop, which are centred around the velocity of each ion beam. They develop due to the scattering of the electrons by the electrostatic wave potentials. The growth of magnetic fields is observed, but their amplitude remains low.

The Early Asymmetries of Supernova 2008D/XRF 080109

Spectropolarimetry of the Type Ib SN 2008D, associated with the X-ray Flash (XRF) 080109, at two separate epochs, is presented. The epochs of these observations correspond to V-band light curve maximum and 15 days after light curve maximum (or 21 and 36 days after the XRF). We find SN 2008D to be significantly polarized, although the largest contribution is due to the interstellar polarization component of Q ISP = 0% ± 0.1% and U ISP = -1.2% ± 0.1%. At the two epochs, the spectropolarimetry of SN 2008D is classified as being D1+L(He I)+L(Ca II). The intrinsic polarization of continuum wavelength regions is <0.4%, at both epochs, implying an asymmetry of the photosphere of <10%. Similar to other Type Ibc SNe, such as 2005bf, 2006aj, and 2007gr, we observed significant polarization corresponding to the spectral features of Ca II, He I, Mg I, Fe II and, possibly, O I ?7774, about a close-to-spherically symmetric photosphere. We introduce a new plot showing the chemically distinct line-forming regions in the ejecta and comment on the apparent ubiquity of highly polarized high-velocity Ca II features in Type Ibc SNe. The polarization angle of Ca II IR triplet was significantly different, at both epochs, to those of the other species, suggesting high-velocity Ca II forms in a separate part of the ejecta. The apparent structure in the outer layers of SN 2008D has implications for the interpretation of the early-time X-ray emission associated with shock breakout. We present two scenarios, within the jet-torus paradigm, which explain the lack of an apparent geometry discontinuity between the two observations: (1) a jet which punched a hole straight through the progenitor and deposited Ni outside the ejecta and (2) a jet which stalled inside the radius of the photosphere as observed at the second epoch. The lack of a peculiar polarization signature, suggesting strongly asymmetric excitation of the ejecta, and the reported properties of the shock-breakout favor the second scenario.

Interpreting the near-infrared spectra of the 'golden standard' Type Ia supernova 2005cf

We present nine near-infrared (NIR) spectra of supernova (SN) 2005cf at epochs from -10 to +42d 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 ~8d before and ~10d 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 Mgii at early epochs and a series of relatively unblended Coii lines at late phases allow us to constrain the regions of the ejecta in which the respective elements are abundant. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

Monte Carlo radiation hydrodynamics:Methods, tests and application to Type Ia supernova ejecta

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. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

Constraining type Ia supernova models:SN 2011fe as A test case

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 (WDs) - 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 Co production. SN 2011fe is close enough to be followed sufficiently long to study this effect. © © 2012 The American Astronomical Society. All rights reserved.

nero- a post-maximum supernova radiation transport code

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 non-local thermodynamic equilibrium (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 500days 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, Fransson & Kozma and to the nebular code of Mazzali et al. 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' Type Ia supernova (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. © 2011 The Authors Monthly Notices of the Royal Astronomical Society © 2011 RAS.

Modeling Type Ia supernova explosions

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

Late-time supernova light curves:The effect of internal conversion and auger electrons

Energy release from radioactive decays contributes significantly to supernova light curves. Previous works, which considered the energy deposited by ?-rays and positrons produced by Ni, Co, Ni, Co, Ti and Sc, have been quite successful in explaining the light curves of both core collapse and thermonuclear supernovae. We point out that Auger and internal conversion electrons, together with the associated X-ray cascade, constitute an additional heat source. When a supernova is transparent to ?-rays, these electrons can contribute significantly to light curves for reasonable nucleosynthetic yields. In particular, the electrons emitted in the decay of Co, which are largely due to internal conversion from a fortuitously low-lying 3/2 state in the daughter Fe, constitute an additional significant energy-deposition channel. We show that when the heating by these electrons is accounted for, a slow-down in the light curve of SN 1998bw is naturally obtained for typical hypernova nucleosynthetic yields. Additionally, we show that for generic Type Ia supernova yields, the Auger electrons emitted in the ground-state to ground-state electron capture decay of Fe exceed the energy released by the Ti decay chain for many years after the explosion. © 2009 RAS.

Cosmic-ray acceleration and escape from supernova remnants

Galactic cosmic-ray (CR) acceleration to the knee in the spectrum at a few PeV is only possible if the magnetic field ahead of a supernova remnant (SNR) shock is strongly amplified by CRs escaping the SNR. A model formulated in terms of the electric charge carried by escaping CRs predicts the maximum CR energy and the energy spectrum of CRs released into the surrounding medium. We find that historical SNRs such as Cas A, Tycho and Kepler may be expanding too slowly to accelerate CRs to the knee at the present time.

TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF THE NONRESONANT, COSMIC-RAY-DRIVEN INSTABILITY IN SUPERNOVA REMNANT SHOCKS

In supernova remnants, the nonlinear amplification of magnetic fields upstream of collisionless shocks is essential for the acceleration of cosmic rays to the energy of the "knee" at 10(15.5) eV. A nonresonant instability driven by the cosmic ray current is thought to be responsible for this effect. We perform two-dimensional, particle-in-cell simulations of this instability. We observe an initial growth of circularly polarized nonpropagating magnetic waves as predicted in linear theory. It is demonstrated that in some cases the magnetic energy density in the growing waves can grow to at least 10 times its initial value. We find no evidence of competing modes, nor of significant modification by thermal effects. At late times, we observe saturation of the instability in the simulation, but the mechanism responsible is an artifact of the periodic boundary conditions and has no counterpart in the supernova-shock scenario.

CI Aql: a Type Ia supernova progenitor?

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

A spectroscopically normal type Ic supernova from a very massive progenitor

We present observations of the Type Ic supernova (SN Ic) 2011bm spanning a period of about one year. The data establish that SN 2011bm is a spectroscopically normal SN Ic with moderately low ejecta velocities and with a very slow spectroscopic and photometric evolution (more than twice as slow as SN 1998bw). The Pan-STARRS1 retrospective detection shows that the rise time from explosion to peak was 40 days in the R band. Through an analysis of the light curve and the spectral sequence, we estimate a kinetic energy of 7-17 foe and a total ejected mass of 7-17 Mo, 5-10 Mo of which is oxygen and 0.6-0.7 Mo is 56Ni. The physical parameters obtained for SN 2011bm suggest that its progenitor was a massive star of initial mass 30-50 Mo. The profile of the forbidden oxygen lines in the nebular spectra show no evidence of a bi-polar geometry in the ejected material.

The Type IIb Supernova 2011dh from a Supergiant Progenitor

A set of hydrodynamical models based on stellar evolutionary progenitors is used to study the nature of SN 2011dh. Our modeling suggests that a large progenitor star ---with R ~200 Rsun---, is needed to reproduce the early light curve of SN 2011dh. This is consistent with the suggestion that the yellow super-giant star detected at the location of the SN in deep pre-explosion images is the progenitor star. From the main peak of the bolometric light curve and expansion velocities we constrain the mass of the ejecta to be ~2 Msun, the explosion energy to be E= 6-10 x 10^50 erg, and the 56Ni mass to be approximately 0.06 Msun. The progenitor star was composed of a helium core of 3 to 4 Msun and a thin hydrogen-rich envelope of ~0.1 M_sun with a main sequence mass estimated to be in the range of 12--15 Msun. Our models rule out progenitors with helium-core masses larger than 8 Msun, which correspond to M_ZAMS > 25 Msun. This suggests that a single star evolutionary scenario for SN 2011dh is unlikely.

PS1-12sk is a Peculiar Supernova From a He-rich Progenitor System in a Brightest Cluster Galaxy Environment

We report on our discovery and observations of the Pan-STARRS1 supernova (SN) PS1-12sk, a transient with properties that indicate atypical star formation in its host galaxy cluster or pose a challenge to popular progenitor system models for this class of explosion. The optical spectra of PS1-12sk classify it as a Type Ibn SN (c.f. SN 2006jc), dominated by intermediate-width (3x10^3 km/s) and time variable He I emission. Our multi-wavelength monitoring establishes the rise time dt = 9-23 days and shows an NUV-NIR SED with temperature > 17x10^3 K and a peak rise magnitude of Mz = -18.9 mag. SN Ibn spectroscopic properties are commonly interpreted as the signature of a massive star (17 - 100 M_sun) explosion within a He-enriched circumstellar medium. However, unlike previous Type Ibn supernovae, PS1-12sk is associated with an elliptical brightest cluster galaxy, CGCG 208-042 (z = 0.054) in cluster RXC J0844.9+4258. The expected probability of an event like PS1-12sk in such environments is low given the measured infrequency of core-collapse SNe in red sequence galaxies compounded by the low volumetric rate of SN Ibn. Furthermore, we find no evidence of star formation at the explosion site to sensitive limits (Sigma Halpha