Subgrid modelling of ion pitch-angle scattering for magnetosheath waves in a global hybrid-Vlasov simulation

Numerical simulations are widely used in modern space physics and are an essential tool to understand or discover new phenomena which cannot be observed using spacecraft measurements. However, numerical simulations are limited by the space grid resolution of the system and the computational costs of having a high spatial resolution. Therefore, some physics may be unresolved in part of the system due to its low spatial resolution. We have previously identified, using Vlasiator, that the proton cyclotron instability is not resolved for grid cell sizes larger than four times the inertial length in the solar wind, for waves in the downstream of the quasi-perpendicular shock in the magnetosheath of a global hybrid-Vlasov simulation. This leads to unphysically high perpendicular temperature and a dominance of the mirror mode waves. In this study, we use high-resolution simulations to measure and quantify how the proton cyclotron instability diffuses and isotropizes the velocity distribution functions. We investigate the process of pitch-angle scattering during the development of the instability and propose a method for the sub-grid modelling of the diffusion process of the instability at low resolution. This allows us to model the isotropization of the velocity distribution functions and to reduce the temperature anisotropy in the plasma while saving computational resources.