The avalanche of runaway electrons is thought to pose a significant obstacle to the success of reactor scale devices such as ITER. As a result, a significant effort has been devoted toward quantifying both the threshold for the initiation of the avalanche of runaway electrons and the efficiency of the avalanche mechanism. In this work, these two quantities are computed utilizing a guiding-center formulation with large-angle collision operators of varying physics fidelity. The use of a guiding-center formulation, while computationally more costly compared to bounce-averaged approaches, provides a conceptually straightforward means of incorporating tokamak geometry. It is found that while the avalanche threshold is only weakly impacted by toroidal geometry for fully ionized low-Z plasmas, it can be significantly impacted if high-Z impurities are present. Furthermore, it is shown that the efficiency of the avalanche mechanism depends sensitively on the impurity content, the charge state of the underlying impurities, and the radial profile of the seed electron population. Finally, the commonly employed Møller secondary source term used to model the generation of secondary electrons is shown to yield avalanche growth rates and thresholds in good agreement with a more complete conservative large-angle collision operator.
The impact of collisionality on the runaway electron avalanche during a tokamak disruption