Monte Carlo implementation of a guiding-center Fokker-Planck kinetic equation

A Monte Carlo method for the collisional guiding-center Fokker-Planck kinetic equation is derived to include the effects of background magnetic-field nonuniformity. It is shown that, in the limit of a homogeneous magnetic field, the guiding-center Monte Carlo collision operators reduce to the standard Coulomb operators. The coefficients required for the implementation of the method are also given.