Ballistic annihilation kinetics: The case of discrete velocity distributions

The kinetics of the annihilation process A+A→O/ with ballistic particle motion is investigated when the distribution of particle velocities is discrete. This discreteness is the source of many intriguing phenomena. In the mean field limit, the densities of different velocity species decay in time with different power law rates for many initial conditions. For a one-dimensional symmetric system containing particles with velocity 0 and ±1, there is a particular initial state for which the concentrations of all three species decay as ${\mathit{t}}^{\mathrm{-}2/3}$, where t is the time. For the case of a fast ‘‘impurity’’ in a symmetric background of + and - particles, the impurity survival probability decays as exp(-const×${\ln }^2$t). In a symmetric four-velocity system in which there are particles with velocities ±${\mathit{v}}_1$ and ±${\mathit{v}}_2$, there again is a special initial condition where the two species decay at the same rate ${\mathit{t}}^{\mathrm{-}\mathrm{\alpha }}$, with α≃0.72. Efficient algorithms are introduced to perform the large-scale simulations necessary to observe these unusual phenomena clearly.