Three-dimensional N-body charged particle simulation employing analytical Kepler solutions in the presence of background magnetic field.

Charged particle diffusion in the presence of magnetic field is studied numerically by a N-body simulation. Central forces (Coulomb forces) between charged particles are computated directly without employing field quantities. When the distance between the two charged particles become comparable to the Landau length (where the kinetic and potential energy becomes comparable), numerical solutions are replaced by analytical Kepler solutions (by regarding the problem as a binary collision problem). The computational method with the analytical solutions are applied in three dimensional computation (by finding out the two-dimensional planes). The electron diffusion coefficient is computed in boundary-less opened domains whose scaling lies between the the classical estimation D_e ~ T_e/ B^2 and the Bohm diffusion D_e ~ T_e/ B, depending on the ratio between the Larmor radius and the Landau length.