Plasmas electron transport properties can theoretically be obtained according to the electric field in two distinct ways. The first method directly solves the electron Boltzmann Equation (BE) by some means of approximation, such as Two-Term Approximation (TTA) [
4,
18,
26,
27]. The second method tracks the random motion of electron bunches and collision with neutral particles by Monte Carlo Code (MCC) [
21,
23,
28]. The TTA technique is used in BOLSIG+ software that is widely used by researchers in the Low-Temperature Plasma (LTP) area, several physical models allow its use in many different working conditions, namely: electron drift and diffusion under the influence of electric and magnetic fields at any angle with each other; density-gradient expansion allowing calculating the transverse and longitudinal bulk/flux swarm parameters; The Influence of Electron-Electron Collisions on the first anisotropy; HF excitement; Pulsed Townsend (PT) and Steady State Townsend (SST) electron density growth due to non-conservative electron scattering mechanisms.
For electronic collisions, EEDF represented by
f, is obtained from solving the Boltzmann Equation (
1):
where
is the electronic speed vector,
is the electric field,
accounts for changes in
f due to the collisions with the species [
4],
e electron charge and
m electron mass.
Since solving the Boltzmann equation directly is a highly challenging task, it can be expanded in spherical coordinates, and the binomial approximation method can be employed to simplify the complexity. In this context, the distribution function, denoted as “
f,” can be expanded as follows:
Here,
represents the angle between the electron velocity vector and the direction of the electric field, while
and
correspond, respectively, to the homogeneity and heterogeneity of the electron energy distribution function, your solution provides valuable foundational data for obtaining various electron transport parameters [
29].
Swarm parameters can be comprehended by considering the spatiotemporal progression of a singular electron avalanche within an unchanging, stable electric field E. In the BOLSIG+ software, these parameters are obtained via a set of mathematical equations. The mean energy in Equation (
4) is determined by integrating the electron energy
with the isotropic portion of the electron distribution function
, which is the zero th-order term in the spherical harmonics expansion of the velocity space. The quantity
, where
e and
represent the electron charge and mass, respectively, and N is the gas density. The Gaussian probability density function characterizing the avalanche is propelled by a center-of-mass velocity Equation (
3) that is determined by the electron mobility Equation (
5), which includes the effective total momentum-transfer cross-section
, including superelastic and electron-ion Coulomb collisions. The avalanche undergoes both transverse and longitudinal diffusion processes, with diffusion coefficients
Equation (
6) and
Equation (
7), respectively, that depend on the first-order component of the density-gradient expansion of the distribution function
, which is a perturbation of
due to the electron density gradient. The reduced Townsend coefficient Equation (
8) determines the first Townsend ionization coefficient and the total number of electrons created per unit length, which is described by the cross-section
, which is the reduced excitation coefficient for process x [
1,
4].
The MCC Method is employed in the METHES program that calculates transport parameters, such as velocity w and diffusion coefficient D for flux Equations (
10) and (
11) and bulk Equations (
12) and (
13) respectively, reaction rates, and EEDF, mean energy Equation (
8) is derived by time-averaging the energies of all electron trajectories for arbitrary gas mixtures where positions
r, velocities
v, through the code written object-oriented in Matlab.