No Arabic abstract
When two galaxy clusters encounter each other, the interaction results in a collisionless shock that is characterized by a low (1-4) sonic Mach number, and a high Alfv{e}nic Mach number. Our goal is to determine if, and to what extent, such shocks can accelerate particles to sufficient velocities that they can contribute to the cosmic ray spectrum. We combine two different computational methods, magnetohydrodynamics (MHD) and particle-in-cell (PIC) into a single code that allows us to take advantage of the high computational efficiency of MHD while maintaining the ability to model the behaviour of individual non-thermal particles. Using this method, we perform a series of simulations covering the expected parameter space of galaxy cluster collision shocks. Our results show that for shocks with a sonic Mach number below 2.25 no diffusive shock acceleration can take place because of a lack of instabilities in the magnetic field, whereas for shocks with a sonic Mach number $geq,3$ the acceleration is efficient and can accelerate particles to relativistic speeds. In the regime between these two extremes, diffusive shock acceleration can occur but is relatively inefficient because of the time- and space-dependent nature of the instabilities. For those shocks that show efficient acceleration, the instabilities in the upstream gas increase to the point where they change the nature of the shock, which, in turn, will influence the particle injection process.
Using large-scale fully-kinetic two-dimensional particle-in-cell simulations, we investigate the effects of shock rippling on electron acceleration at low-Mach-number shocks propagating in high-$beta$ plasmas, in application to merger shocks in galaxy clusters. We find that the electron acceleration rate increases considerably when the rippling modes appear. The main acceleration mechanism is stochastic shock-drift acceleration, in which electrons are confined at the shock by pitch-angle scattering off turbulence and gain energy from the motional electric field. The presence of multi-scale magnetic turbulence at the shock transition and the region immediately behind the main shock overshoot is essential for electron energization. Wide-energy non-thermal electron distributions are formed both upstream and downstream of the shock. The maximum energy of the electrons is sufficient for their injection into diffusive shock acceleration. We show for the first time that the downstream electron spectrum has a~power-law form with index $papprox 2.5$, in agreement with observations.
Giant radio halos are Mpc-scale diffuse radio sources associated with the central regions of galaxy clusters. The most promising scenario to explain the origin of these sources is that of turbulent re-acceleration, in which MeV electrons injected throughout the formation history of galaxy clusters are accelerated to higher energies by turbulent motions mostly induced by cluster mergers. In this Letter, we use the amplitude of density fluctuations in the intracluster medium as a proxy for the turbulent velocity and apply this technique to a sample of 51 clusters with available radio data. Our results indicate a segregation in the turbulent velocity of radio halo and radio quiet clusters, with the turbulent velocity of the former being on average higher by about a factor of two. The velocity dispersion recovered with this technique correlates with the measured radio power through the relation $P_{rm radio}proptosigma_v^{3.3pm0.7}$, which implies that the radio power is nearly proportional to the turbulent energy rate. Our results provide an observational confirmation of a key prediction of the turbulent re-acceleration model and possibly shed light on the origin of radio halos.
Understanding the thermodynamic state of the hot intracluster medium (ICM) in a galaxy cluster requires a knowledge of the plasma transport processes, especially thermal conduction. The basic physics of thermal conduction in plasmas with ICM-like conditions has yet to be elucidated, however. We use particle-in-cell simulations and analytic models to explore the dynamics of an ICM-like plasma (with small gyroradius, large mean-free-path, and strongly sub-dominant magnetic pressure) driven by the diffusive heat flux associated with thermal conduction. Lin- ear theory reveals that whistler waves are driven unstable electron heat flux, even when the heat flux is weak. The resonant interaction of electrons with these waves then plays a critical role in scattering electrons and suppressing the heat flux. In a 1D model where only whistler modes that are parallel to the magnetic field are captured, the only resonant electrons are moving in the opposite direction to the heat flux and the electron heat flux suppression is small. In 2D or more, oblique whistler modes also resonate with electrons moving in the direction of the heat flux. The overlap of resonances leads to effective symmetrization of the electron distribution function and a strong suppression of heat flux. The results suggest that thermal conduction in the ICM might be strongly suppressed, possibly to negligible levels.
We herein investigate shock formation and particle acceleration processes for both protons and electrons in a quasi-parallel high-Mach-number collisionless shock through a long-term, large-scale particle-in-cell simulation. We show that both protons and electrons are accelerated in the shock and that these accelerated particles generate large-amplitude Alfv{e}nic waves in the upstream region of the shock. After the upstream waves have grown sufficiently, the local structure of the collisionless shock becomes substantially similar to that of a quasi-perpendicular shock due to the large transverse magnetic field of the waves. A fraction of protons are accelerated in the shock with a power-law-like energy distribution. The rate of proton injection to the acceleration process is approximately constant, and in the injection process, the phase-trapping mechanism for the protons by the upstream waves can play an important role. The dominant acceleration process is a Fermi-like process through repeated shock crossings of the protons. This process is a `fast process in the sense that the time required for most of the accelerated protons to complete one cycle of the acceleration process is much shorter than the diffusion time. A fraction of the electrons is also accelerated by the same mechanism, and have a power-law-like energy distribution. However, the injection does not enter a steady state during the simulation, which may be related to the intermittent activity of the upstream waves. Upstream of the shock, a fraction of the electrons is pre-accelerated before reaching the shock, which may contribute to steady electron injection at a later time.
An extreme case of electron shock drift acceleration in low Mach number collisionless shocks is investigated as a plausible mechanism of initial acceleration of relativistic electrons in large-scale shocks in galaxy clusters where upstream plasma temperature is of the order of 10 keV and a degree of magnetization is not too small. One-dimensional electromagnetic full particle simulations reveal that, even though a shock is rather moderate, a part of thermal incoming electrons are accelerated and reflected through relativistic shock drift acceleration and form a local nonthermal population just upstream of the shock. The accelerated electrons can self-generate local coherent waves and further be back-scattered toward the shock by those waves. This may be a scenario for the first stage of the electron shock acceleration occurring at the large-scale shocks in galaxy clusters such as CIZA J2242.8+5301 which has well defined radio relics.