“…Early simulations were one dimensional in space and used hybrid particle‐in‐cell (hybrid‐PIC) algorithms [e.g., Leroy et al , ]. They soon moved to full particle simulations [e.g., Lembège and Dawson , ] and increased in dimensionality to 2‐D hybrid‐ and full‐PIC simulations [e.g., Lembège and Savoini , ; Savoini et al , ], gaining in complexity and detail of the physical description at large scales [e.g., Lin , ; Omidi et al , ; Blanco‐Cano et al , ] as well as for local kinetic processes [e.g., Ofman and Gedalin , ] and even reaching three spatial dimensions for hybrid‐PIC simulations of the magnetosphere [e.g., Lin and Wang , ]. PIC‐based methods suffer from the statistical noise inherent to the random sampling of the particle distribution injected in the simulation and the low number of particles that can feasibly be simulated.…”