Analysis and simulation for a model of electron impact excitation/deexcitation and ionization/recombination

This paper describes a kinetic model and a corresponding Monte Carlo simulation method for excitation/deexcitation and ionization/recombination by electron impact in a plasma free of external fields. The atoms and ions in the plasma are represented by continuum densities and the electrons by a particle distribution. A Boltzmann-type equation is formulated and a corresponding H-theorem is formally derived. An efficient Monte Carlo method is developed for an idealized analytic model of the excitation and ionization collision cross sections. To accelerate the simulation, the reduced rejection method and binary search method are used to overcome the singular rate in the recombination process. Numerical results are presented to demonstrate the efficiency of the method on spatially homogeneous problems. The evolution of the electron distribution function and atomic states is studied, revealing the possibility under certain circumstances of system relaxation towards stationary states that are not the equilibrium states, a potential non-ergodic behavior.