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.

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.

Multi-Dimensional Simulations of the Expanding Supernova Remnant of SN 1987A

The expanding remnant from SN 1987A is an excellent laboratory for investigating the physics of supernovae explosions. There is still a large number of outstanding questions, such as the reason for the asymmetric radio morphology, the structure of the pre-supernova environment, and the efficiency of particle acceleration at the supernova shock. We explore these questions using three-dimensional simulations of the expanding remnant between days 820 and 10,000 after the supernova. We combine a hydrodynamical simulation with semi-analytic treatments of diffusive shock acceleration and magnetic field amplification to derive radio emission as part of an inverse problem. Simulations show that an asymmetric explosion, combined with magnetic field amplification at the expanding shock, is able to replicate the persistent one-sided radio morphology of the remnant. We use an asymmetric Truelove & McKee progenitor with an envelope mass of 10 M-circle dot and an energy of 1.5 x 10(44) J. A termination shock in the progenitor's stellar wind at a distance of 0 ''.43-0 ''.51 provides a good fit to the turn on of radio emission around day 1200. For the H II region, a minimum distance of 0 ''.63 +/- 0 ''.01 and maximum particle number density of (7.11 +/- 1.78) x 10(7) m(-3) produces a good fit to the evolving average radius and velocity of the expanding shocks from day 2000 to day 7000 after explosion. The model predicts a noticeable reduction, and possibly a temporary reversal, in the asymmetric radio morphology of the remnant after day 7000, when the forward shock left the eastern lobe of the equatorial ring.

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.

Galactic Edge Clouds I. Molecular Line Observations and Chemical Modelling of Edge Cloud 2

Edge Cloud 2 (EC2) is a molecular cloud, about 35 pc in size, with one of the largest galactocentric distances known to exist in the Milky Way. We present observations of a peak CO emission region in the cloud and use these to determine its physical characteristics. We calculate a gas temperature of 20 K and a density of n(H2)~10^4 cm-3. Based on our CO maps, we estimate the mass of EC2 at around 10^4 Msolar and continuum observations suggest a dust-to-gas mass ratio as low as 0.001. Chemical models have been developed to reproduce the abundances in EC2, and they indicate that heavy element abundances may be reduced by a factor of 5 relative to the solar neighborhood (similar to dwarf irregular galaxies and damped Lya systems), very low extinction (A_V <4 mag) due to a very low dust-to-gas mass ratio, an enhanced cosmic-ray ionization rate, and a higher UV field compared to local interstellar values. The reduced abundances may be attributed to the low level of star formation in this region and are probably also related to the continuing infall of primordial (or low-metallicity) halo gas since the Milky Way formed. Finally, we note that shocks from the old supernova remnant GSH 138-01-94 may have determined the morphology and dynamics of EC2.

Magnetic instability in a dilute circular rarefaction wave

The growth of magnetic fields in the density gradient of a rarefaction wave has been observed in simulations and in laboratory experiments. The thermal anisotropy of the electrons, which gives rise to the magnetic instability, is maintained by the ambipolar electric field. This simple mechanism could be important for the magnetic field amplification in astrophysical jets or in the interstellar medium ahead of supernova remnant shocks. The acceleration of protons and the generation of a magnetic field by the rarefaction wave, which is fed by an expanding circular plasma cloud, is examined here in form of a 2D particle-in-cell simulation. The core of the plasma cloud is modeled by immobile charges, and the mobile protons form a small ring close to the cloud's surface. The number density of mobile protons is thus less than that of the electrons. The protons of the rarefaction wave are accelerated to 1/10 of the electron thermal speed, and the acceleration results in a thermal anisotropy of the electron distribution in the entire plasma cloud. The instability in the rarefaction wave is outrun by a TM wave, which grows in the dense core distribution, and its magnetic field expands into the rarefaction wave. This expansion drives a secondary TE wave. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4769128]

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. 

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.

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.

Evolution of slow electrostatic shock into a plasma shock mediated by electrostatic turbulence

The collision of two plasma clouds at a speed that exceeds the ion acoustic speed can result in the formation of shocks. This phenomenon is observed not only in astrophysical scenarios, such as the propagation of supernova remnant (SNR) blast shells into the interstellar medium, but also in laboratory-based laser-plasma experiments. These experiments and supporting simulations are thus seen as an attractive platform for small-scale reproduction and study of astrophysical shocks in the laboratory. We model two plasma clouds, which consist of electrons and ions, with a 2D particle-in-cell simulation. The ion temperatures of both clouds differ by a factor of ten. Both clouds collide at a speed that is realistic for laboratory studies and for SNR shocks in their late evolution phase, like that of RCW86. A magnetic field, which is orthogonal to the simulation plane, has a strength that is comparable to that of SNR shocks. A forward shock forms between the overlap layer of both plasma clouds and the cloud with cooler ions. A large-amplitude ion acoustic wave is observed between the overlap layer and the cloud with hotter ions. It does not steepen into a reverse shock because its speed is below the ion acoustic speed. A gradient of the magnetic field amplitude builds up close to the forward shock as it compresses the magnetic field. This gradient gives rise to an electron drift that is fast enough to trigger an instability. Electrostatic ion acoustic wave turbulence develops ahead of the shock, widens its transition layer, and thermalizes the ions, but the forward shock remains intact. © 2014 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.

Turbulent amplification of magnetic fields in laboratory laser-produced shock waves

X-ray and radio observations of the supernova remnant Cassiopeia A reveal the presence of magnetic fields about 100 times stronger than those in the surrounding interstellar medium. Field coincident with the outer shock probably arises through a nonlinear feedback process involving cosmic rays. The origin of the large magnetic field in the interior of the remnant is less clear but it is presumably stretched and amplified by turbulent motions. Turbulence may be generated by hydrodynamic instability at the contact discontinuity between the supernova ejecta and the circumstellar gas9. However, optical observations of Cassiopeia A indicate that the ejecta are interacting with a highly inhomogeneous, dense circumstellar cloud bank formed before the supernova explosion. Here we investigate the possibility that turbulent amplification is induced when the outer shock overtakes dense clumps in the ambient medium. We report laboratory experiments that indicate the magnetic field is amplified when the shock interacts with a plastic grid. We show that our experimental results can explain the observed synchrotron emission in the interior of the remnant. The experiment also provides a laboratory example of magnetic field amplification by turbulence in plasmas, a physical process thought to occur in many astrophysical phenomena.

The binary progenitor of Tycho Brahe's 1572 supernova

The brightness of type Ia supernovae, and their homogeneity as a class, makes them powerful tools in cosmology, yet little is known about the progenitor systems of these explosions. They are thought to arise when a white dwarf accretes matter from a companion star, is compressed and undergoes a thermonuclear explosion(1-3). Unless the companion star is another white dwarf ( in which case it should be destroyed by the mass-transfer process itself), it should survive and show distinguishing properties. Tycho's supernova(4,5) is one of only two type Ia supernovae observed in our Galaxy, and so provides an opportunity to address observationally the identification of the surviving companion. Here we report a survey of the central region of its remnant, around the position of the explosion, which excludes red giants as the mass donor of the exploding white dwarf. We found a type G0 - G2 star, similar to our Sun in surface temperature and luminosity ( but lower surface gravity), moving at more than three times the mean velocity of the stars at that distance, which appears to be the surviving companion of the supernova.

Deflagrations in hybrid CONe white dwarfs: a route to explain the faint Type Iax supernova 2008ha

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 similar to 1.1 M-aS (TM), hybrid WDs in a binary system may easily approach the Chandrasekhar mass (M-Ch) by accretion and give rise to a thermonuclear explosion. Here, we investigate an off-centre deflagration in a near-M-Ch 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 M-aS (TM) 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-M-Ch WDs can in principle explain the observed diversity of Type Iax supernovae. Leaving behind a near-M-Ch 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 M-Ch to be of the order of 1 per cent of the Galactic SN Ia rate.