Numerical scheme for evaluating the collision integrals for triple interactions in relativistic plasma

TL;DR

This paper introduces a new numerical scheme for accurately computing collision integrals for triple interactions in relativistic plasma, improving the modeling of plasma kinetics in astrophysical and laboratory contexts.

Contribution

A novel conservative numerical scheme based on direct integration of QED matrix elements for triple interactions in relativistic plasma.

Findings

Reaction rates agree with analytic expressions where available.

The scheme effectively models relativistic plasma out of equilibrium.

Results are applicable to astrophysical and laboratory plasma studies.

Abstract

Binary interactions in relativistic plasma, such as Coulomb and Compton scattering as well as pair creation and annihilation are well known and studied in details. Triple interactions, namely relativistic bremsstrahlung, double Compton scattering, radiative pair production, triple pair production/annihilation and their inverse processes, are usually considered as emission processes in astrophysical problems, as well as in laboratory plasmas. Their role in plasma kinetics is fundamental [1]. We present a new conservative scheme for computation of Uehling-Uhlenbeck collision integral for all triple interactions in relativistic plasma based on direct integration of exact QED matrix elements. Reaction rates for thermal distributions are compared, where possible, with the corresponding analytic expressions, showing good agreement. Our results are relevant for quantitative description of relativistic plasmas out of equilibrium, both in astrophysical and laboratory conditions.