Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOn magnetic field amplification and particle acceleration near non-relativistic astrophysical shocks: Particles in MHD Cells simulations
252
0
0.0
(
0
)
Added by
Allard Jan van Marle
Publication date
2017
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We present simulations of magnetized astrophysical shocks taking into account the interplay between the thermal plasma of the shock and supra-thermal particles. Such interaction is depicted by combining a grid-based magneto-hydrodynamics description of the thermal fluid with particle in cell techniques devoted to the dynamics of supra-thermal particles. This approach, which incorporates the use of adaptive mesh refinement features, is potentially a key to simulate astrophysical systems on spatial scales that are beyond the reach of pure particle-in-cell simulations. We consider in this study non-relativistic shocks with various Alfvenic Mach numbers and magnetic field obliquity. We recover all the features of both magnetic field amplification and particle acceleration from previous studies when the magnetic field is parallel to the normal to the shock. In contrast with previous particle-in-cell-hybrid simulations, we find that particle acceleration and magnetic field amplification also occur when the magnetic field is oblique to the normal to the shock but on larger timescales than in the parallel case. We show that in our simulations, the supra-thermal particles are experiencing acceleration thanks to a pre-heating process of the particle similar to a shock drift acceleration leading to the corrugation of the shock front. Such oscillations of the shock front and the magnetic field locally help the particles to enter the upstream region and to initiate a non-resonant streaming instability and finally to induce diffuse particle acceleration.
rate research
Read More
We use particle-in-magnetohydrodynamics-cells to model particle acceleration and magnetic field amplification in a high Mach, parallel shock in three dimensions and compare the result to 2-D models. This allows us to determine whether 2-D simulations can be relied upon to yield accurate results in terms of particle acceleration, magnetic field amplification and the growth rate of instabilities. Our simulations show that the behaviour of the gas and the evolution of the instabilities are qualitatively similar for both the 2-D and 3-D models, with only minor quantitative differences that relate primarily to the growth speed of the instabilities. The main difference between 2-D and 3-D models can be found in the spectral energy distributions (SEDs) of the non-thermal particles. The 2-D simulations prove to be more efficient, accelerating a larger fraction of the particles and achieving higher velocities. We conclude that, while 2-D models are sufficient to investigate the instabilities in the gas, their results have to be treated with some caution when predicting the expected SED of a given shock.
Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic electron-positron jet front propagating into an ambient electron-positron plasma with and without initial magnetic fields. We find small differences in the results for no ambient and modest ambient magnetic fields. New simulations show that the Weibel instability created in the collisionless shock front accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. Furthermore, the non-linear fluctuation amplitudes of densities, currents, electric, and magnetic fields in the electron-positron shock are larger than those found in the electron-ion shock studied in a previous paper at the comparable simulation time. This comes from the fact that both electrons and positrons contribute to generation of the Weibel instability. Additionally, we have performed simulations with different electron skin depths. We find that growth times scale inversely with the plasma frequency, and the sizes of structures created by the Weibel instability scale proportional to the electron skin depth. This is the expected result and indicates that the simulations have sufficient grid resolution. The simulation results show that the Weibel instability is responsible for generating and amplifying nonuniform, small-scale magnetic fields which contribute to the electrons (positrons) transverse deflection behind the jet head.
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.
The Fermi LAT discovery that classical novae produce >100 MeV gamma-rays establishes that shocks and relativistic particle acceleration are key features of these events. These shocks are likely to be radiative due to the high densities of the nova ejecta at early times coincident with the gamma-ray emission. Thermal X-rays radiated behind the shock are absorbed by neutral gas and reprocessed into optical emission, similar to Type IIn (interacting) supernovae. Gamma-rays are produced by collisions between relativistic protons with the nova ejecta (hadronic scenario) or Inverse Compton/bremsstrahlung emission from relativistic electrons (leptonic scenario), where in both scenarios the efficiency for converting relativistic particle energy into LAT gamma-rays is at most a few tens of per cent. The ratio of gamma-ray and optical luminosities, L_gam/L_opt, thus sets a lower limit on the fraction of the shock power used to accelerate relativistic particles, e_nth. The measured values of L_gam/L_opt for two classical novae, V1324 Sco and V339 Del, constrains e_nth > 1e-2 and > 1e-3, respectively. Inverse Compton models for the gamma-ray emission are disfavored given the low electron acceleration efficiency, e_nth ~ 1e-4-1e-3, inferred from observations of Galactic cosmic rays and particle-in-cell (PIC) numerical simulations. A fraction > 100(0.01/e_nth) and > 10(0.01/e_nth) per cent of the optical luminosity is powered by shocks in V1324 Sco and V339 Del, respectively. Such high fractions challenge standard models that instead attribute all nova optical emission to the direct outwards transport of thermal energy released near the white dwarf surface.
Synchrotron radio emission from non-relativistic jets powered by massive protostars has been reported, indicating the presence of relativistic electrons and magnetic fields of strength ~0.3-5 mG. We study diffusive shock acceleration and magnetic field amplification in protostellar jets with speeds between 300 and 1500 km/s. We show that the magnetic field in the synchrotron emitter can be amplified by the non-resonant hybrid (Bell) instability excited by the cosmic-ray streaming. By combining the synchrotron data with basic theory of Bell instability we estimate the magnetic field in the synchrotron emitter and the maximum energy of protons. Protons can achieve maximum energies in the range 0.04-0.65 TeV and emit gamma rays in their interaction with matter fields. We predict detectable levels of gamma rays in IRAS 16547-5247 and IRAS 16848-4603. The gamma ray flux can be significantly enhanced by the gas mixing due to Rayleigh-Taylor instability. The detection of this radiation by the Fermi satellite in the GeV domain and the forthcoming Cherenkov Telescope Array at higher energies may open a new window to study the formation of massive stars, as well as diffusive acceleration and magnetic field amplification in shocks with velocities of about 1000 km/s.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
