Recent analytical studies and particle-in-cell simulations suggested that the electron velocity distribution function in a Hall thruster plasma is non-Maxwellian and anisotropic. 1,2 The electron average kinetic energy in the direction parallel to walls is several times larger than the electron average kinetic energy in direction normal to the walls. Electrons are stratified into several groups depending on their origin (e.g., plasma discharge or thruster channel walls) and confinement (e.g., lost on the walls or trapped in the plasma). Practical analytical formulas are derived for wall fluxes, secondary electron fluxes, plasma parameters, and conductivity. The calculations based on analytical formulas agree well with the results of numerical simulations. The self-consistent analysis demonstrates that elastic electron scattering on collisions with atoms and ions plays a key role in formation of the electron energy distribution function and plasma-wall interaction. The fluxes of electrons from the plasma bulk are shown to be proportional to the rate of scattering to loss cone, thus collision frequency determines the wall potential and secondary electron fluxes. Secondary electron emission from the walls is shown to enhance the electron conductivity across the magnetic field, while having almost no effect on insulating properties of the near-wall sheaths. Such a self-consistent decoupling between secondary electron emission effects on electron energy losses and electron crossed-field transport is currently not captured by the existing fluid and hybrid models of the Hall thrusters.

The increasing need to demonstrate the correctness of computer simulations has highlighted the importance of benchmarks. We define in this paper a representative simulation case to study low-temperature partially-magnetized plasmas. Seven independently developed Particle-In-Cell codes have simulated this benchmark case, with the same specified conditions. The characteristics of the codes used, such as implementation details or computing times and resources, are given. First, we compare at steady-state the time-averaged axial profiles of three main discharge parameters (axial electric field, ion density and electron temperature). We show that the results obtained exhibit a very good agreement within 5% between all the codes. As ExB discharges are known to cause instabilities propagating in the direction of electron drift, an analysis of these instabilities is then performed and a similar behaviour is retrieved between all the codes. A particular attention has been paid to the numerical convergence by varying the number of macroparticles per cell and we show that the chosen benchmark case displays a good convergence. Detailed outputs are given in the supplementary data, to be used by other similar codes in the perspective of code verification. 2D axial-azimuthal Particle-In-Cell benchmark for low-temperature partially ...

The Electron Cyclotron Drift Instability (ECDI) driven by the electron E × B drift in partially magnetized plasmas is investigated with highly resolved particle-in-cell simulations. The emphasis is on two-dimensional effects involving the parallel dynamics along the magnetic field in a finite length plasma with dielectric walls. It is found that the instability develops as a sequence of growing cyclotron harmonics demonstrating wave breaking and complex nonlinear interactions, being particularly pronounced in ion density fluctuations at short wavelengths. At the same time, nonlinear evolution of fluctuations of the ion and electron density, as well as the anomalous electron current, shows cascade toward long wavelengths. Tendency to generate long wavelength components is most clearly observed in the spectra of the electron density and the anomalous current fluctuations. An intense but slowly growing mode with a distinct eigen-mode structure along the magnetic field develops at a later nonlinear stage enhancing the tendency toward long wavelength condensation. The latter mode having a finite wavelength along the magnetic field is identified as the Modified Two-Stream Instability (MTSI). It is shown that the MTSI mode results in strong parallel heating of electrons.

A one-dimensional kinetic theory of sheaths surrounding planar, electron-emitting surfaces is presented which accounts for plasma electrons lost to the surface and the temperature of the emitted electrons. It is shown that ratio of plasma electron temperature to emitted electron temperature significantly affects the sheath potential when the plasma electron temperature is within an order of magnitude of the emitted electron temperature. The sheath potential goes to zero as the plasma electron temperature equals the emitted electron temperature, which can occur in the afterglow of an rf plasma and some low-temperature plasma sources. These results were validated by particle in cell simulations. The theory was tested by making measurements of the sheath surrounding a thermionically emitting cathode in the afterglow of an rf plasma. The measured sheath potential shrunk to zero as the plasma electron temperature cooled to the emitted electron temperature, as predicted by the theory.