Quasineutral particle simulation technique for whistlers

We present a new hybrid fluid/PIC simulation scheme for whistlers, which eliminates both the speed-of-light time scale and the electron plasma oscillation time scale, and concentrates simulation resources on the resonant parts of electron phase space that control whistler evolution. The code runs with time steps on the order of the electron gyrofrequency, with extremely accurate energy conservation and numerical stability. Examples are shown of application to whistler instability growth and saturation, and ducting of whistlers in density channels.