No Arabic abstract
Galactic cosmic rays are believed to be accelerated at supernova remnants. However, whether supernova remnants can be Pevatrons is still very unclear. In this work we argue that PeV cosmic rays can be accelerated during the early phase of a supernova blast wave expansion in dense red supergiant winds. We solve in spherical geometry a system combining a diffusive-convection equation which treats cosmic-ray dynamics coupled to magnetohydrodynamics to follow gas dynamics. The fast shock expanding in a dense ionized wind is able to trigger the fast non-resonant streaming instability over day timescales, and energizes cosmic-rays even under the effect of p-p losses. We find that such environments make the blast wave a Pevatron, although the maximum energy depends on various parameters such as the injection rate and mass-loss rate of the winds. Multi-PeV energies can be reached if the progenitor mass loss rates are of the order of $10^{-3}$ Msun yr$^{-1}$. It has been recently invoked that, prior to the explosion, hydrogen rich massive stars can produce enhanced mass loss rates. These enhanced rates would then favor the production of a Pevatron phase in early times after the shock breakout.
We have performed magnetohydrodynamical simulations to study the amplification of magnetic fields in the precursors of shock waves. Strong magnetic fields are required in the precursors of the strong shocks that occur in supernova remnants. Observations also suggest that magnetic field amplification takes place in the weak shocks that occur in galaxy clusters and that produce so-called radio relics. Here, we extend the study of magnetic field amplification by cosmic-ray driven turbulence to weak shocks. The amplification is driven by turbulence that is produced by the cosmic-ray pressure acting on the density inhomogeneities in the upstream fluid. The clumping that has been inferred from X-ray data for the outskirts of galaxy clusters could provide some of the seed inhomogeneities. Magnetic field power spectra and Faraday maps are produced. Furthermore, we investigate how the synchrotron emission in the shock precursor can be used to verify the existence of this instability and constrain essential plasma parameters.
We present a nonlinear Monte Carlo model of efficient diffusive shock acceleration (DSA) where the magnetic turbulence responsible for particle diffusion is calculated self-consistently from the resonant cosmic-ray (CR) streaming instability, together with non-resonant short- and long-wavelength CR-current-driven instabilities. We include the backpressure from CRs interacting with the strongly amplified magnetic turbulence which decelerates and heats the super-alfvenic flow in the extended shock precursor. Uniquely, in our plane-parallel, steady-state, multi-scale model, the full range of particles, from thermal (~eV) injected at the viscous subshock, to the escape of the highest energy CRs (~PeV) from the shock precursor, are calculated consistently with the shock structure, precursor heating, magnetic field amplification (MFA), and scattering center drift relative to the background plasma. In addition, we show how the cascade of turbulence to shorter wavelengths influences the total shock compression, the downstream proton temperature, the magnetic fluctuation spectra, and accelerated particle spectra. A parameter survey is included where we vary shock parameters, the mode of magnetic turbulence generation, and turbulence cascading. From our survey results, we obtain scaling relations for the maximum particle momentum and amplified magnetic field as functions of shock speed, ambient density, and shock size.
Cosmic ray (CR) currents through magnetised plasma drive strong instabilities producing amplification of the magnetic field. This amplification helps explain the CR energy spectrum as well as observations of supernova remnants and radio galaxy hot spots. Using magnetohydrodynamic (MHD) simulations, we study the behaviour of the non-resonant hybrid (NRH) instability (also known as the Bell instability) in the case of CR currents perpendicular and parallel to the initial magnetic field. We demonstrate that extending simulations of the perpendicular case to 3D reveals a different character to the turbulence from that observed in 2D. Despite these differences, in 3D the perpendicular NRH instability still grows exponentially far into the non-linear regime with a similar growth rate to both the 2D perpendicular and 3D parallel situations. We introduce some simple analytical models to elucidate the physical behaviour, using them to demonstrate that the transition to the non-linear regime is governed by the growth of thermal pressure inside dense filaments at the edges of the expanding loops. We discuss our results in the context of supernova remnants and jets in radio galaxies. Our work shows that the NRH instability can amplify magnetic fields to many times their initial value in parallel and perpendicular shocks.
Collisionless shocks are ubiquitous in the Universe and often associated with strong magnetic field. Here we use large-scale particle-in-cell simulations of non-relativistic perpendicular shocks in the high-Mach-number regime to study the amplification of magnetic field within shocks. The magnetic field is amplified at the shock transition due to the ion-ion two-stream Weibel instability. The normalized magnetic-field strength strongly correlates with the Alfvenic Mach number. Mock spacecraft measurements derived from PIC simulations are fully consistent with those taken in-situ at Saturns bow shock by the Cassini spacecraft.
We model the multi-wavelength emission in the southern hotspot of the radio quasar 4C74.26. The synchrotron radio emission is resolved near the shock with the MERLIN radio-interferometer, and the rapid decay of this emission behind the shock is interpreted as the decay of the amplified downstream magnetic field as expected for small scale turbulence. Electrons are accelerated to only 0.3 TeV, consistent with a diffusion coefficient many orders of magnitude greater than in the Bohm regime. If the same diffusion coefficient applies to the protons, their maximum energy is only ~100 TeV.