N-body charged particle simulation in two- and three-dimensional systems

A computational tool is developed to investigate charged particle transport in the presence of magnetic fields [D. Bohm, The characteristics of electrical discharges in magnetic fields, New York: McGraw-Hill (1949)] based on the first principle based N-body particle simulation. Instead of employing conventional softening parameters at the binary collisions, we have employed the analytical Kepler solutions when the kinetic energy and the Coulomb potential's energy became comparable, which are successfully connected to the numerical solutions of the N-body simulation. We assume that the interactions from the rest of the N-2 particles are small compared to the binary interactions (the effect of the central forces from the N-2 particles can be incorporated perturbatively). The work is extended to a three-dimensional system. At the binary collisions, reminding Kepler's first law, we find out local planes where the two orbits reside. In a boundary-less system, particles' diffusion coefficients are estimated by varying the strength of the background magnetic field.