Numerical analysis of Debye shielding by N-body simulation

Basic properties of Debye shielding is investigated by directly calculating Coulomb force between the charged particles. Learning from N-body gravity simulation, numerical techniques are employed to avoid singularities when potential energy between two charged particles become comparable to the kinetic energy. Electron distribution function in Debye shielding is analyzed, which deviates from conventional Maxwellian distribution function.