Dynamics of ballistic annihilation

Jaroslaw Piasecki, Emmanuel Trizac 1, 2, Michel Droz

Physical Review E: Statistical, Nonlinear, and Soft Matter Physics 66 (2002) 066111

The problem of ballistically controlled annihilation is revisited for general initial velocity distributions and arbitrary dimension. An analytical derivation of the hierarchy equations obeyed by the reduced distributions is given, and a scaling analysis of the corresponding spatially homogeneous system is performed. This approach points to the relevance of the non-linear Boltzmann equation for dimensions larger than one and provides expressions for the exponents describing the decay of the particle density n(t) ~ t^{-\xi} and the root mean-square velocity ${\bar v} ~ t^{-\gamma}$ in term of a parameter related to the dissipation of kinetic energy. The Boltzmann equation is then solved perturbatively within a systematic expansion in Sonine polynomials. Analytical expressions for the exponents $\xi$ and $\gamma$ are obtained in arbitrary dimension as a function of the parameter $\mu$ characterizing the small velocity behavior of the initial velocity distribution. Moreover, the leading non-Gaussian corrections to the scaled velocity distribution are computed. These expressions for the scaling exponents are in good agreement with the values reported in the literature for continuous velocity distributions in $d=1$. For the two dimensional case, we implement Monte-Carlo and molecular dynamics simulations that turn out to be in excellent agreement with the analytical predictions.

  • 1. Laboratoire de Physique Théorique d'Orsay (LPT),
    CNRS : UMR8627 – Université Paris XI - Paris Sud
  • 2. Laboratoire de Physique Théorique et Modèles Statistiques (LPTMS),
    CNRS : UMR8626 – Université Paris XI - Paris Sud