Normal type Ia supernovae from violent mergers of white dwarf binaries

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 M and 1.1 M 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 light curves of our merger, which produces about 0.62 M of 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. © 2012 The American Astronomical Society. All rights reserved.

Time-dependent three-dimensional spectrum synthesis for Type Ia supernovae

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 ?-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 comoving frame. Assuming a photoionization-dominated plasma, the equations of ionization 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. © 2009 RAS.

Off-center explosions of Chandrasekhar-mass white dwarfs:An explanation of super-bright type la supernovae?

Context. The recent discovery of a very bright type la supernova, SNLS-03D3bb (=SN 2003fg), in the Supernova Legacy Survey (SNLS) has raised the question of whether super-Chandrasekhar-mass white-dwarf stars are needed to explain such bright explosions. Progenitors of this sort could form by mergers of pairs of rather massive white dwarfs. Binary systems of two white dwarfs in close orbit, where their total mass significantly exceeds the Chandrasekhar mass, have not yet been found. Therefore SNLS-03D3bb could establish the first clear case of a double-degenerate progenitor of a (peculiar) type la supernovae. Moreover, if this interpretation is correct, it casts some doubt on the universality of the calibration relations used to make SNe la distance indicators for cosmology. Aims. We aim to evaluate the case for a super-Chandrasekhar-mass progenitor for SNLS-03D3bb in light of previous theoretical work on super-Chandrasekhar-mass explosions. Furthermore, we propose an alternative scenario involving only a Chandrasekhar-mass progenitor. Methods. We present a theoretically motivated critical discussion of the expected observational fingerprints of super-Chandrasekharmass explosions. As an alternative, we describe a simple class of aspherical Chandrasekhar-mass models in which the products of nuclear burning are displaced from the center. We then perform simple radiative transfer calculations to predict synthetic lightcurves for one such off-center explosion model. Results. In important respects, the expected observational consequences of super-Chandrasekhar-mass explosions are not consistent with the observations of SNLS-03D3bb. We demonstrate that the lopsided explosion of a Chandrasekhar-mass white dwarf could provide a better explanation. © ESO 2007.

3D deflagration simulations leaving bound remnants:A model for 2002cx-like Type Ia supernovae

2002cx-like supernovae are a sub-class of sub-luminous Type Ia supernovae (SNe). 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 prototypical 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 d 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.03Mconsisting 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. © 2013 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society.

Three-dimensional delayed-detonation models with nucleosynthesis for type ia supernovae

We present results for a suite of 14 three-dimensional, high-resolution hydrodynamical simulations of delayed-detonation models of 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 data base for Chandrasekhar-mass delayed-detonation model nucleosynthetic yields and for deriving synthetic observables such as spectra and light curves. We employ aphysically motivated, stochastic model based on turbulent velocity fluctuations and fuel density to calculate in situ the deflagration-to-detonation transition 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 a central density of 2.9 × 10 g cm, as well as one high central density (5.5 × 10 g cm) and one low central density (1.0 × 10 g cm) rendition of the 100 ignition kernel configuration. For each simulation, we determined detailed nucleosynthetic yields by postprocessing10 tracer particles with a 384 nuclide reaction network. All delayed-detonation models result in explosions unbinding thewhite dwarf, producing a range of 56Ni masses from 0.32 to 1.11M. As a general trend, the models predict that the stableneutron-rich iron-group isotopes are not found at the lowest velocities, but rather at intermediate velocities (~3000×10 000 km s) in a shell surrounding a Ni-rich core. The models further predict relatively low-velocity oxygen and carbon, with typical minimum velocities around 4000 and 10 000 km s, respectively. © 2012 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.

A filamentation instability for streaming cosmic rays

We demonstrate that cosmic rays form filamentary structures in the precursors of supernova remnant shocks due to their self-generated magnetic fields. The cosmic ray filamentation results in the growth of a long-wavelength instability, and naturally couples the rapid non-linear amplification on small scales to larger length-scales. Hybrid magnetohydrodynamics-particle simulations are performed to confirm the effect. The resulting large-scale magnetic field may facilitate the scattering of high-energy cosmic rays as required to accelerate protons beyond the knee in the cosmic ray spectrum at supernova remnant shocks. Filamentation far upstream of the shock may also assist in the escape of cosmic rays from the accelerator.

Parametric study of non-relativistic electrostatic shocks and the structure of their transition layer

Nonrelativistic electrostatic unmagnetized shocks are frequently observed in laboratory plasmas and they are likely to exist in astrophysical plasmas. Their maximum speed, expressed in units of the ion acoustic speed far upstream of the shock, depends only on the electron-to-ion temperature ratio if binary collisions are absent. The formation and evolution of such shocks is examined here for a wide range of shock speeds with particle-in-cell simulations. The initial temperatures of the electrons and the 400 times heavier ions are equal. Shocks form on electron time scales at Mach numbers between 1.7 and 2.2. Shocks with Mach numbers up to 2.5 form after tens of inverse ion plasma frequencies. The density of the shock-reflected ion beam increases and the number of ions crossing the shock thus decreases with an increasing Mach number, causing a slower expansion of the downstream region in its rest frame. The interval occupied by this ion beam is on a positive potential relative to the far upstream. This potential pre-heats the electrons ahead of the shock even in the absence of beam instabilities and decouples the electron temperature in the foreshock ahead of the shock from the one in the far upstream plasma. The effective Mach number of the shock is reduced by this electron heating. This effect can potentially stabilize nonrelativistic electrostatic shocks moving as fast as supernova remnant shocks. 

Time-Resolved Characterization of the Formation of a Collisionless Shock

We report on the temporally and spatially resolved detection of the precursory stages that lead to the formation of an unmagnetized, supercritical collisionless shock in a laser-driven laboratory experiment. The measured evolution of the electrostatic potential associated with the shock unveils the transition from a current free double layer into a symmetric shock structure, stabilized by ion reflection at the shock front. Supported by a matching particle-in-cell simulation and theoretical considerations, we suggest that this process is analogous to ion reflection at supercritical collisionless shocks in supernova remnants.

Universal behaviour of shock precursors in the presence of efficient cosmic ray acceleration

The self-consistent interaction between energetic particles and self-generated hydromagnetic waves in a cosmic ray pressure dominated plasma is considered. Using a three-dimensional hybrid magnetohydrodynamics (MHD)-kinetic code, which utilizes a spherical harmonic expansion of the Vlasov-Fokker-Planck equation, high-resolution simulations of the magnetic field growth including feedback on the cosmic rays are carried out. It is found that for shocks with high cosmic ray acceleration efficiency, the magnetic fields become highly disorganized, resulting in near isotropic diffusion, independent of the initial orientation of the ambient magnetic field. The possibility of sub-Bohm diffusion is demonstrated for parallel shocks, while the diffusion coefficient approaches the Bohm limit from below for oblique shocks. This universal behaviour suggests that Bohm diffusion in the root-mean-squared field inferred from observation may provide a realistic estimate for the maximum energy acceleration time-scale in young supernova remnants. Although disordered, the magnetic field is not self-similar suggesting a non-uniform energy-dependent behaviour of the energetic particle transport in the precursor. Possible indirect radiative signatures of cosmic ray driven magnetic field amplification are discussed.

Self-organized electromagnetic field structures in laser-produced counter-streaming plasmas

Self-organization(1,2) occurs in plasmas when energy progressively transfers from smaller to larger scales in an inverse cascade(3). Global structures that emerge from turbulent plasmas can be found in the laboratory(4) and in astrophysical settings; for example, the cosmic magnetic field(5,6,) collisionless shocks in supernova remnants(7) and the internal structures of newly formed stars known as Herbig-Haro objects(8). Here we show that large, stable electromagnetic field structures can also arise within counter-streaming supersonic plasmas in the laboratory. These surprising structures, formed by a yet unexplained mechanism, are predominantly oriented transverse to the primary flow direction, extend for much larger distances than the intrinsic plasma spatial scales and persist for much longer than the plasma kinetic timescales. Our results challenge existing models of counter-streaming plasmas and can be used to better understand large-scale and long-time plasma self-organization.

Cosmic ray acceleration at oblique shocks

We show that the diffusion approximation breaks down for particle acceleration at oblique shocks with velocities typical of young supernova remnants. Higher order anisotropies flatten the spectral index at quasi-parallel shocks and steepen the spectral index at quasi-perpendicular shocks. We compare the theory with observed spectral indices.

A current-driven instability in parallel, relativistic shocks

Recently, Bell ( 2004 Mon. Not. R. Astron. Soc. 353 550) has reanalysed the problem of wave excitation by cosmic rays propagating in the pre-cursor region of a supernova remnant shock front. He pointed out a strong, non-resonant, current-driven instability that had been overlooked in the kinetic treatments by Achterberg ( 1983 Astron. Astrophys. 119 274) and McKenzie and Volk ( 1982 Astron. Astrophys. 116 191), and suggested that it is responsible for substantial amplification of the ambient magnetic field. Magnetic field amplification is also an important issue in the problem of the formation and structure of relativistic shock fronts, particularly in relation to models of gamma-ray bursts. We have therefore generalized the linear analysis to apply to this case, assuming a relativistic background plasma and a monoenergetic, unidirectional incoming proton beam. We find essentially the same non-resonant instability observed by Bell and show that also, under GRB conditions, it grows much faster than the resonant waves. We quantify the extent to which thermal effects in the background plasma limit the maximum growth rate.

The transport of cosmic rays in self-excited magnetic turbulence

The process of diffusive shock acceleration relies on the efficacy with which hydromagnetic waves can scatter charged particles in the precursor of a shock. The growth of self-generated waves is driven by both resonant and non-resonant processes. We perform high-resolution magnetohydrodynamic simulations of the non-resonant cosmic ray driven instability, in which the unstable waves are excited beyond the linear regime. In a snapshot of the resultant field, particle transport simulations are carried out. The use of a static snapshot of the field is reasonable given that the Larmor period for particles is typically very short relative to the instability growth time. The diffusion rate is found to be close to, or below, the Bohm limit for a range of energies. This provides the first explicit demonstration that self-excited turbulence reduces the diffusion coefficient and has important implications for cosmic-ray transport and acceleration in supernova remnants.

A cosmic ray current-driven instability in partially ionised media

Context. We investigate the growth of hydromagnetic waves driven by streaming cosmic rays in the precursor environment of a supernova remnant shock.

Aims. It is known that transverse waves propagating parallel to the mean magnetic field are unstable to anisotropies in the cosmic ray distribution, and may provide a mechanism to substantially amplify the ambient magnetic field. We quantify the extent to which temperature and ionisation fractions modify this picture.

Methods. Using a kinetic description of the plasma we derive the dispersion relation for a collisionless thermal plasma with a streaming cosmic ray current. Fluid equations are then used to discuss the effects of neutral-ion collisions.

Results. We calculate the extent to which the environment into which the cosmic rays propagate influences the growth of the magnetic field, and determines the range of possible growth rates.

Conclusions. If the cosmic ray acceleration is efficient, we find that very large neutral fractions are required to stabilise the growth of the non-resonant mode. For typical supernova parameters in our Galaxy, thermal effects do not significantly alter the growth rates. For weakly driven modes, ion-neutral damping can dominate over the instability at more modest ionisation fractions. In the case of a supernova shock interacting with a molecular clouds, such as in RX J1713.7-3946, with high density and low ionisation, the modes can be rapidly damped.

STEADY-STATE SOLUTIONS IN NONLINEAR DIFFUSIVE SHOCK ACCELERATION

Stationary solutions to the equations of nonlinear diffusive shock acceleration play a fundamental role in the theory of cosmic-ray acceleration. Their existence usually requires that a fraction of the accelerated particles be allowed to escape from the system. Because the scattering mean free path is thought to be an increasing function of energy, this condition is conventionally implemented as an upper cutoff in energy space-particles are then permitted to escape from any part of the system, once their energy exceeds this limit. However, because accelerated particles are responsible for the substantial amplification of the ambient magnetic field in a region upstream of the shock front, we examine an alternative approach in which particles escape over a spatial boundary. We use a simple iterative scheme that constructs stationary numerical solutions to the coupled kinetic and hydrodynamic equations. For parameters appropriate for supernova remnants, we find stationary solutions with efficient acceleration when the escape boundary is placed at the point where growth and advection of strongly driven nonresonant waves are in balance. We also present the energy dependence of the distribution function close to the energy where it cuts off-a diagnostic that is in principle accessible to observation.