Exact Timestep for a Pairwise Coulomb Collision

Abstract

Standard numerical integrators work well for many-body Coulomb repulsion problems when the timestep is much shorter than the timescale of relative position changes. However, for ‘hard’ collisions in which two particles have a near miss and exchange a lot of momentum within one timestep, they understandably struggle. This note proposes using the exact solution of Keplerian two-body orbits (usually hyperbolic) to calculate the momentum exchange with other particles: either a selection of the ‘closest’ ones or all of them.