Ion Weibel Instability in the hybrid framework: the optimal resolution

TL;DR

This paper develops a linear theory and simulation validation for ion Weibel instability in hybrid models, providing guidelines for optimal resolution to accurately capture shock microphysics.

Contribution

It introduces a resolution criterion for hybrid simulations of ion Weibel instability, ensuring accurate modeling of shock microinstabilities in collisionless plasma.

Findings

Hybrid simulations can reliably reproduce Weibel instability features with proper resolution.

A minimum spatial resolution depends on the Alfvénic Mach number.

Excessive resolution can produce unphysical small-scale modes.

Abstract

The study of collisionless shocks and their role in cosmic-ray acceleration has gained increasing importance through both observations and simulations. Accurately modeling the shock transition region, where particle injection occurs, requires a proper description of the microinstabilities governing its structure. In high-Mach-number shocks, such as those associated with supernova remnants, the ion Weibel instability is believed to provide the dominant dissipation mechanism. In this work, we investigate the ion Weibel instability driven by counterstreaming beams in the presence of an external perpendicular magnetic field. We employ hybrid simulations, in which ions are treated kinetically while electrons are modeled as a charge-neutralizing fluid. Although hybrid models are widely employed to study collisionless shocks, the resolution requirements needed to accurately capture ion-scale instabilities remain poorly understood. We address this issue by developing a linear theory of the ion Weibel instability tailored to the massless electron assumption of hybrid models and validating it with one- and two-dimensional simulations over a wide range of Alfv\'enic Mach numbers. We show that hybrid simulations can reliably reproduce the growth, saturation, and polarization of Weibel-generated magnetic fields in weakly magnetized regimes, provided that the relevant ion-scale modes are properly resolved. From the scaling of the dominant mode, we derive a minimum spatial resolution required as a function of Alfv\'enic Mach number. We also demonstrate that excessive resolution introduces unphysical small-scale whistler modes inherent to the massless-electron approximation. We validate the analysis by comparing the results with full particle-in-cell simulations. Together, these results provide practical guidance for hybrid simulations of collisionless shocks and beam-driven plasma systems.