New analytic approach to multivelocity annihilation in the kinetic theory of reactions

Probabilistic ballistic annihilation with continuous velocity distributions

We investigate the problem of ballistically controlled reactions where particles either annihilate upon collision with probability p, or undergo an elastic shock with probability 1-p. Restricting to homogeneous systems, we provide in the scaling regime that emerges in the long time limit, analytical expressions for the exponents describing the time decay of the density and the root-mean-square velocity, as continuous functions of the probability p and of a parameter related to the dissipation of energy. We work at the level of molecular chaos (nonlinear Boltzmann equation), and using a systematic Sonine polynomials expansion of the velocity distribution, we obtain in arbitrary dimension the first non-Gaussian correction and the corresponding expressions for the decay exponents. We implement Monte Carlo simulations in two dimensions, which are in excellent agreement with our analytical predictions. For p<1, numerical simulations lead to the conjecture that unlike for pure annihilation (p=1), the velocity distribution becomes universal, i.e., does not depend on the initial conditions.

Some exact results for Boltzmann's annihilation dynamics

The problem of ballistic annihilation for a spatially homogeneous system is revisited within Boltzmann's kinetic theory in two and three dimensions. Analytical results are derived for the time evolution of the particle density for some isotropic discrete bimodal velocity modulus distributions. According to the allowed values of the velocity modulus, different behaviors are obtained: power law decay with nonuniversal exponents depending continuously upon the ratio of the two velocities, or exponential decay. When one of the two velocities is equal to zero, the model describes the problem of ballistic annihilation in the presence of static traps. The analytical predictions are shown to be in agreement with the results of two-dimensional molecular dynamics simulations.

Dynamics of ballistic annihilation

The problem of ballistically controlled annihilation is revisited for general initial velocity distributions and an 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 nonlinear Boltzmann equation for dimensions larger than 1 and provides expressions for the exponents describing the decay of the particle density n(t) proportional, variant t(-xi) and the root-mean-square velocity v proportional, variant t(-gamma) in terms 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.

Kinetics and scaling in ballistic annihilation

We study the simplest irreversible ballistically controlled reaction, whereby particles having an initial continuous velocity distribution annihilate upon colliding. In the framework of the Boltzmann equation, expressions for the exponents characterizing the density and typical velocity decay are explicitly worked out in arbitrary dimension. These predictions are in excellent agreement with the complementary results of extensive Monte Carlo and molecular dynamics simulations. We finally discuss the definition of universality classes indexed by a continuous parameter for this far from equilibrium dynamics with no conservation laws.

Stochastic ballistic annihilation and coalescence

We study a class of stochastic ballistic annihilation and coalescence models with a binary velocity distribution in one dimension. We obtain an exact solution for the density which reveals a universal phase diagram for the asymptotic density decay. By universal we mean that all models in the class are described by a single phase diagram spanned by two reduced parameters. The phase diagram reveals four regimes, two of which contain the previously studied cases of ballistic annihilation. The two new phases are a direct consequence of the stochasticity. The solution is obtained through a matrix product approach and builds on properties of a q-deformed harmonic oscillator algebra.