Ion Weibel Instability in the hybrid framework: the optimal resolution

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énic 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énic 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.

• The paper establishes resolution criteria in hybrid simulations to capture the ion Weibel instability, prescribing a minimum grid resolution that scales with Alfvénic Mach number.
• Systematic 1D and 2D hybrid and PIC simulations validate that under-resolving the ion scales suppresses instability growth while excessive resolution introduces non-physical whistler modes.
• The study prescribes a minimum resolution scaling law and a maximum limit of 30 cells per ion skin depth to ensure physical fidelity in modeling astrophysical shocks.

Optimal Resolution Criteria for the Ion Weibel Instability in Hybrid Simulations

Introduction

The paper "Ion Weibel Instability in the hybrid framework: the optimal resolution" (2604.05021) addresses a critical and previously unquantified issue in kinetic plasma theory: establishing resolution requirements for hybrid simulations to faithfully capture the ion Weibel instability, a key mechanism mediating dissipation and particle acceleration in high-Mach-number astrophysical shocks. Hybrid codes, in which ions are treated kinetically and electrons as a massless charge-neutralizing fluid, are central tools in modeling large-scale plasma environments such as those found at supernova remnants. However, the interplay between microinstability physics and numerical resolution has often been neglected, resulting in simulations that either under-resolve physical modes or are contaminated by numerically-induced artifacts.

This work develops a linear theory for the ion Weibel instability in the massless-electron (hybrid) limit, validates this framework through systematic one- and two-dimensional hybrid and PIC simulations, and prescribes both necessary and maximum spatial resolution as explicit functions of Alfvénic Mach number. These results are essential for accurate modeling of shock transition regions, particle injection, and acceleration processes in astrophysical and laboratory plasmas.

Linear Theory and Resolution Criteria

The theoretical analysis considers counterstreaming ion beams in the presence of a perpendicular background magnetic field. The hybrid model assumes electrons are a charge-neutralizing fluid with negligible inertia. The linear dispersion relation is analytically derived, and distinctions between the hybrid and full kinetic electron regimes are highlighted, especially regarding electron screening.

The framework demonstrates that, in weakly magnetized regimes, the ion Weibel instability’s fastest-growing mode has a wavelength that decreases with increasing Mach number. For hybrid models, neglecting electron mass is well-motivated only when the instability evolves faster than the ion gyroperiod but much more slowly than the electron gyroperiod. The theoretical spectrum yields key scaling laws for the dominant wavenumber:

Figure 1: Linear theory for $M_{\mathrm{A}}=30$ and $M_{\mathrm{A}}=100$, with growth rates, saturation, and polarization for the Weibel instability.

A precise fit ($k_{\mathrm{peak}} \propto M_{\mathrm{A}}^{0.51}$) to the fastest-growing wavenumber is found, establishing the minimum resolution requirement as:

$$N_{\mathrm{min}} = 9 \left(\frac{M_{\mathrm{A}}}{30}\right)^{0.51}$$

where $N$ is the number of grid cells per ion skin depth $d_i$. The required number of resolution elements thus grows sublinearly with Mach number. For instance, $M_{\mathrm{A}}=100$ demands $N \approx 17$, while $M_{\mathrm{A}}=30$ is consistent with $N=10$.

Figure 2: Dependence of peak mode wavenumber on $M_{\mathrm{A}}=100$0; orange denotes the linear theory, black the best fit $M_{\mathrm{A}}=100$1 scaling.

Polarization analysis reveals that the $M_{\mathrm{A}}=100$2 and $M_{\mathrm{A}}=100$3 nature of filamentary fluctuations also depends sensitively on resolution, with clear cross-over points for each Mach regime.

Limitations: Maximum Resolution and Whistler Mode Contamination

A pivotal insight of the analysis is that excessive spatial resolution in hybrid models introduces small-scale whistler modes that do not correspond to physical plasma behavior. The massless electron assumption causes an artificial divergence in whistler phase velocity at high wavenumbers, leading to numerically unstable or nonphysical solutions, in contrast to full kinetic models with finite electron mass.

Figure 3: Dispersion relation and phase velocity of whistler modes in hybrid and full kinetic models. Massless electrons cause artificial high phase velocities at short wavelengths.

This necessitates a practical maximum resolution in hybrid simulations, identified as $M_{\mathrm{A}}=100$4 cells per $M_{\mathrm{A}}=100$5, independent of $M_{\mathrm{A}}=100$6, to avoid unphysical small-scale whistler excitation. Therefore, for a given Mach number:

$M_{\mathrm{A}}=100$7

Beyond $M_{\mathrm{A}}=100$8, no physically meaningful spatial resolution exists within hybrid models, and fully kinetic treatments become necessary.

Validation with 1D and 2D Simulations

Extensive 1D and 2D hybrid and PIC simulations using the dHybridR code corroborate the theoretical predictions. In the 1D case, for $M_{\mathrm{A}}=100$9, a resolution of $k_{\mathrm{peak}} \propto M_{\mathrm{A}}^{0.51}$0 accurately resolves the dominant $k_{\mathrm{peak}} \propto M_{\mathrm{A}}^{0.51}$1-mode; under-resolving ($k_{\mathrm{peak}} \propto M_{\mathrm{A}}^{0.51}$2) suppresses growth rates, shifts spectral power, and distorts magnetic field polarization.

Figure 4: 1D hybrid simulation results showing power spectra and polarization for $k_{\mathrm{peak}} \propto M_{\mathrm{A}}^{0.51}$3 and $k_{\mathrm{peak}} \propto M_{\mathrm{A}}^{0.51}$4. Convergence is reached as predicted by theory.

For $k_{\mathrm{peak}} \propto M_{\mathrm{A}}^{0.51}$5, $k_{\mathrm{peak}} \propto M_{\mathrm{A}}^{0.51}$6 is sufficient, with accurate agreement in growth, saturation, and polarization characteristics with the theoretical model.

2D simulations further demonstrate that, for $k_{\mathrm{peak}} \propto M_{\mathrm{A}}^{0.51}$7, hybrid results match the full PIC solutions at ion scales, while $k_{\mathrm{peak}} \propto M_{\mathrm{A}}^{0.51}$8 exhibits spurious small-scale and oblique fluctuations incompatible with accurate modeling.

Figure 5: 2D hybrid and full PIC results for $k_{\mathrm{peak}} \propto M_{\mathrm{A}}^{0.51}$9. Hybrid at optimal ($$N_{\mathrm{min}} = 9 \left(\frac{M_{\mathrm{A}}}{30}\right)^{0.51}$$0) and excessive ($$N_{\mathrm{min}} = 9 \left(\frac{M_{\mathrm{A}}}{30}\right)^{0.51}$$1) resolutions compared to PIC with realistic mass ratio.

Implications and Future Directions

The implications for shock acceleration studies and plasma turbulence modeling are immediate. Simulations that under-resolve the critical scales will misrepresent the growth, saturation, and topology of Weibel-induced fields, which directly modulates particle injection, cosmic-ray acceleration, and shock reformation. Conversely, over-resolving, or using hybrid codes outside validated Mach number ranges, can contaminate transition-region structure and field statistics via artificial whistler mode power.

The practical guidelines derived here provide parameterization for future large-scale hybrid shock and beam instabilities studies, ensuring physical fidelity while maximizing computational efficiency. For higher Mach regimes or studies probing sub-ion scales, more advanced models should be adopted. Recent developments in hybrid-kinetic and extended fluid models, including finite electron mass [Munoz2018, Jain2023], and closure approaches incorporating kinetic corrections [Hunana2022, Jikei2021], may enable hybrid studies at higher $$N_{\mathrm{min}} = 9 \left(\frac{M_{\mathrm{A}}}{30}\right)^{0.51}$$2 and offer improved fidelity for small-scale phenomena.

Conclusion

This work solidifies the theoretical and numerical foundation for simulating ion Weibel-driven field evolution in hybrid codes, explicitly quantifying both the minimum and maximum spatial resolution as a function of plasma parameters. These results establish best practices for hybrid modeling of high-Mach-number shocks, providing prescriptive boundaries that safeguard against both physical under-resolution and numerical pathology. The extension of hybrid methods through improved electron models is a promising research avenue for expanding the physically meaningful parameter space of hybrid simulations in collisionless shock and astrophysical plasma research.