Test particle simulation is a useful method for studying both linear and nonlinear wave-particle interactions in the magnetosphere. The gyro-averaged equations of particle motion for first-order and other cyclotron harmonic resonances with oblique whistler-mode waves were first derived by Bell [J. Geophys. Res. 89, 905 (1984)] and the most recent relativistic form was given by Ginet and Albert [Phys. Fluids B 3, 2994 (1991)], and Bortnik [Ph.D. thesis (Stanford University, 2004), p. 40]. However, recently we found there was a (- 1) l - 1 term difference between their formulas of perpendicular motion for the lth-order resonance. This article presents the detailed derivation process of the generalized resonance formulas, and suggests a check of the signs for self-consistency, which is independent of the choice of conventions, that is, the energy variation equation resulting from the momentum equations should not contain any wave magnetic components, simply because the magnetic field does not contribute to changes of particle energy. In addition, we show that the wave centripetal force, which was considered small and was neglect in previous studies of nonlinear interactions, has a profound time derivative and can significantly enhance electron phase trapping especially in high frequency waves. This force can also bounce the low pitch angle particles out of the loss cone. We justify both the sign problem and the missing wave centripetal force by demonstrating wave-particle interaction examples, and comparing the gyro-averaged particle motion to the full particle motion under the Lorentz force.

Radiation belt electrons can exhibit different types of pitch angle distributions in response to various magnetospheric processes. Butterfly distributions, characterized by flux minima at pitch angles around 90°, are broadly observed in both the outer and inner belts and the slot region. Butterfly distributions close to the outer magnetospheric boundary have been attributed to drift shell splitting and losses to the magnetopause. However, their occurrence in the inner belt and the slot region has hitherto not been resolved. By analyzing the particle and wave data collected by the Van Allen Probes during a geomagnetic storm, we combine test particle calculations and Fokker-Planck simulations to reveal that scattering by equatorial magnetosonic waves is a significant cause for the formation of energetic electron butterfly distributions in the inner magnetosphere. Another event shows that a large-amplitude magnetosonic wave in the outer belt can create electron butterfly distributions in just a few minutes.