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.

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.

A Stubbornly Large Mass of Cold Dust in the Ejecta of Supernova 1987A

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 mu m data and improved imaging quality at 100 and 160 mu m compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 mu m [O-I] line flux, eliminating the possibility that line contaminations distort the previously estimated dustmass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 mu m flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5 +/- 0.1M(circle dot) of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3M(circle dot) of amorphous carbon and 0.5M(circle dot) of silicates, totalling 0.8M(circle dot) of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

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

Local star formation triggered by supernova shocks in magnetized diffuse neutral clouds

In this work, considering the impact of a supernova remnant (SNR) with a neutral magnetized cloud we derived analytically a set of conditions that are favourable for driving gravitational instability in the cloud and thus star formation. Using these conditions, we have built diagrams of the SNR radius, R(SNR), versus the initial cloud density, n(c), that constrain a domain in the parameter space where star formation is allowed. This work is an extension to previous study performed without considering magnetic fields (Melioli et al. 2006, hereafter Paper I). The diagrams are also tested with fully three-dimensional MHD radiative cooling simulations involving a SNR and a self-gravitating cloud and we find that the numerical analysis is consistent with the results predicted by the diagrams. While the inclusion of a homogeneous magnetic field approximately perpendicular to the impact velocity of the SNR with an intensity similar to 1 mu G within the cloud results only a small shrinking of the star formation zone in the diagram relative to that without magnetic field, a larger magnetic field (similar to 10 mu G) causes a significant shrinking, as expected. Though derived from simple analytical considerations these diagrams provide a useful tool for identifying sites where star formation could be triggered by the impact of a supernova blast wave. Applications of them to a few regions of our own Galaxy (e.g. the large CO shell in the direction of Cassiopeia, and the Edge Cloud 2 in the direction of the Scorpious constellation) have revealed that star formation in those sites could have been triggered by shock waves from SNRs for specific values of the initial neutral cloud density and the SNR radius. Finally, we have evaluated the effective star formation efficiency for this sort of interaction and found that it is generally smaller than the observed values in our own Galaxy (SFE similar to 0.01-0.3). This result is consistent with previous work in the literature and also suggests that the mechanism presently investigated, though very powerful to drive structure formation, supersonic turbulence and eventually, local star formation, does not seem to be sufficient to drive global star formation in normal star-forming galaxies, not even when the magnetic field in the neutral clouds is neglected.

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.

An Integral View of Fast Shocks Around Supernova 1006

Supernova remnants are among the most spectacular examples of astrophysical pistons in our cosmic neighborhood. The gas expelled by the supernova explosion is launched with velocities ~1000 kilometers per second into the ambient, tenuous interstellar medium, producing shocks that excite hydrogen lines. We have used an optical integral-field spectrograph to obtain high-resolution spatial-spectral maps that allow us to study in detail the shocks in the northwestern rim of supernova 1006. The two-component Hα line is detected at 133 sky locations. Variations in the broad line widths and the broad-to-narrow line intensity ratios across tens of atomic mean free paths suggest the presence of suprathermal protons, the potential seed particles for generating high-energy cosmic rays.

Investigating the cosmic-ray ionization rate in the galactic interstellar medium through observations of H3+

Observations of H3+ in the Galactic diffuse interstellar medium (ISM) have led to various surprising results, including the conclusion that the cosmic-ray ionization rate (zeta_2) is about 1 order of magnitude larger than previously thought. The present survey expands the sample of diffuse cloud sight lines with H3+ observations to 50, with detections in 21 of those. Ionization rates inferred from these detections are in the range (1.7+-1.0)x10^-16 s^-1 < zeta_2 < (10.6+-6.8)x10^-16 s^-1 with a mean value of zeta_2=(3.3+-0.4)x10^-16 s^-1. Upper limits (3sigma) derived from non-detections of H3+ are as low as zeta_2 < 0.4x10^-16 s^-1. These low upper-limits, in combination with the wide range of inferred cosmic-ray ionization rates, indicate variations in zeta_2 between different diffuse cloud sight lines. Calculations of the cosmic-ray ionization rate from theoretical cosmic-ray spectra require a large flux of low-energy (MeV) particles to reproduce values inferred from observations. Given the relatively short range of low-energy cosmic rays --- those most efficient at ionization --- the proximity of a cloud to a site of particle acceleration may set its ionization rate. Variations in zeta_2 are thus likely due to variations in the cosmic-ray spectrum at low energies resulting from the effects of particle propagation. To test this theory, H3+ was observed in sight lines passing through diffuse molecular clouds known to be interacting with the supernova remnant IC 443, a probable site of particle acceleration. Where H3+ is detected, ionization rates of zeta_2=(20+-10)x10^-16 s^-1 are inferred, higher than for any other diffuse cloud. These results support both the concept that supernova remnants act as particle accelerators, and the hypothesis that propagation effects are responsible for causing spatial variations in the cosmic-ray spectrum and ionization rate. Future observations of H3+ near other supernova remnants and in sight lines where complementary ionization tracers (OH+, H2O+, H3O+) have been observed will further our understanding of the subject.

Electron collisions with Fe-peak elements: Forbidden transitions between the low lying valence states 3d6, 3d54s, and 3d54p of Fe III:

Effective collision strengths are presented for the Fe-peak element Fe III at electron temperatures (Te in degrees Kelvin) in the range 2 × 103 to 1 × 106. Forbidden transitions results are given between the 3d6, 3d54s, and the 3d54p manifolds applicable to the modeling of laboratory and astrophysical plasmas.

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.

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.