Development of gravitational N-body simulation algorithm for charged particle dynamics

Charged particles'  transport process is investigated employing a gravitational N-body numerical simulation algorithm. To circumvent severe computational restrictions at the close encounters of two charged particles, we employ analytical Kepler solutions which are smoothly connected to the numerical solutions. The connection is at the order ot Landau length. We also discuss how to incorporate the magnetic field perturbatively. Having alike charges in the system and thus Debye shielding, plasma N-body simulation is expected to be less demanding in constructing hierarchies compared to the gravitational N-body simulation whose central forces are all attractive. After examining the mechanism of the like-charged-particles' collisions and alike-charged-particles collsions, the statistical behavior of the system (particle's diffusion rate) is examined in comparison with the Bohm's analysis which is given by "D ~T/B".