Quasi-perpendicular shocks of galaxy clusters in hybrid kinetic simulations

Abstract

Context. Cosmic ray acceleration in galaxy clusters is still an ongoing puzzle, with relativistic electrons forming radio relics at merger shocks and emitting synchrotron radiation. These shocks are also potential sources of ultra-high-energy cosmic rays, gamma rays, and neutrinos. Our recent work focuses on electron acceleration at low Mach number merger shocks in the hot intracluster medium which is characterized by high plasma beta. Using particle-in-cell (PIC) simulations, we previously showed that electrons are energized through the stochastic shock-drift acceleration process, which is facilitated by multi-scale turbulence, including ion-scale shock surface rippling. For the present work, we performed hybrid-kinetic simulations in a range of various quasi-perpendicular foreshock conditions, including plasma beta, magnetic obliquity, and the shock Mach number. Aims. We study the ion kinetic physics, which is responsible for the shock structure and wave turbulence, that in turn affects the particle acceleration processes. We cover the spatial and temporal scales, which allow the development of large-scale ion turbulence modes in the system. Methods. We applied a recently developed generalized fluid-particle hybrid numerical code that can combine fluid modeling for both electrons and ions with an arbitrary number of kinetic species. We limited this model to a standard hybrid simulation configuration with kinetic ions and fluid electrons. The model utilizes the exact form of the generalized Ohm’s law, allowing for an arbitrary choice of mass and energy densities, as well as the charge-to-mass ratio of the kinetic species. Results. We show that the properties of ion-driven multi-scale magnetic turbulence in merger shocks are in agreement with the ion structures observed in PIC simulations. In typical shocks with the sonic Mach number Ms = 3, the magnetic structures and shock front density ripples grow and saturate at wavelengths reaching approximately four ion Larmor radii. Only shocks with Ms ≳ 2.3 develop ripples. At very weak shocks with Ms ≲ 2.3, weak turbulence is formed downstream of the shock. We observed a moderate dependence of the strength of magnetic field fluctuations on the quasi-perpendicular magnetic field obliquity. However, as the field obliquity decreases, the shock front ripples exhibit longer wavelengths. Finally, we note that the steady-state structure of Ms = 3 shocks in high-beta plasmas shows evidence that there is little difference between 2D and 3D simulations. The turbulence near the shock front seems to be a 2D-like structure in 3D simulations.