Electron Velocity Distribution Research Articles

Effect of temperature anisotropy formed by fast electron beams moving in the flare loop on its excited electron-cyclotron maser instability

Context. The electron-cyclotron maser instability (ECMI) is a significant coherent radio emission mechanism widely utilized in various astrophysical radio phenomena. It is well known that the velocity anisotropic distribution of energetic electrons, which leads to an inverted perpendicular population in the vertical direction with ∂fb/∂v⊥ > 0, can provide the free energy necessary for the ECMI. Aims. The initial velocity distribution of energetic electrons leaving the acceleration region is typically isotropic or beam-like. However, as these energetic electrons travel along the magnetic field as fast electron beams (FEBs) in magnetic plasma, various velocity anisotropic distributions can emerge. In this paper, we examine the impact of temperature anisotropy formed by beam electrons traveling along a flare loop on the ECMI. Methods. By neglecting the energy loss of energetic electrons as they traverse the corona and invoking the conservation of energy and magnetic moments, we established the relationship between momentum dispersion and the magnetic field. Utilizing the magnetic field model of the flare loop, we calculated the evolution of momentum dispersion and the growth rates of the ECMI as FEBs precipitate along the flare loop. Results. The results demonstrate that the temperature anisotropy arising as FEBs descend along the flare loop significantly impacts the ECMI. The maximum growth rates of the excited modes exhibit a gradual increase initially and then decline rapidly after reaching a critical height for β0 = 0.2c and 0.15c. The results also show that the growth rates of the O2 mode are one order of magnitude smaller than those of the O1 and X2 modes. This indicates that the harmonic radiation is X-mode polarized. Notably, the temperature anisotropy of FEBs as they precipitate along the flare loop with different magnetic field models or at different heights has similar effects on the ECMI.

Whistler wave scattering of energetic electrons past 90°

The consequences of a 90° barrier in the scattering of energetic electrons by whistler waves are explored with self-consistent two-dimensional particle-in-cell simulations. In the presence of a 90° scattering barrier, a field-aligned heat flux of energetic electrons will rapidly scatter to form a uniform distribution with pitch angles 0<θ<90° but with a discontinuous jump at θ=90° to a lower energy distribution of electrons with 90°<θ<180°. However, simulations reveal that such a distribution contains a large reservoir of free energy that is released to drive large-amplitude, oblique-propagating whistler waves (δB/B0∼0.1). As a result, energetic electrons near a pitch angle 90° experience strong resonance scattering. Nearly half of the energetic electrons in the positive parallel velocity plane cross the 90° barrier and diffuse to negative parallel velocities. Thus, the late-time electron velocity distribution becomes nearly isotropic. This result has implications for understanding the regulation of energetic particle heat flux in space and astrophysical environments, including the solar corona, the solar wind, and the intracluster medium of galaxy clusters.

Importance of the electron plasma parameter for excitation of chorus and formation of magnetic field irregularity in the region of their excitation

A threshold condition has been established for excitation of VLF electromagnetic radiation with a chorus structure of the dynamic spectrum in the daytime magnetosphere using the BPA (Beam Pulse Amplifier) mechanism for amplifying short noise electromagnetic pulses. The kappa distribution was used as a model function of the electron velocity distribution in the magnetosphere. Calculations performed for this distribution have shown that the threshold for excitation of chorus largely depends on the electron plasma parameter equal to the ratio of gas-kinetic pressure of electrons to magnetic pressure. This pattern is not contradicted by the dependence of the probability of excitation of chorus on the degree of magnetic field irregularity, which we derived from observations made by the Van Allen Probe spacecraft. It is sharp fluctuations in the magnetic field strength near its local minima outside the plasmasphere where the radiation under study can be excited. If there is an irregularity, the probability of detecting chorus is >70 %; and if there is no or very low irregularity, the probability of the absence of any emissions is ~80 %. The results indicate a common reason for the excitation of chorus and the magnetic field irregularity — a small but finite value of the plasma parameter.

Существенность величины плазменного параметра электронов для возбуждения хоров и формирования нерегулярности магнитного поля в области их возбуждения

A threshold condition has been established for excitation of VLF electromagnetic radiation with a chorus structure of the dynamic spectrum in the daytime magnetosphere using the BPA (Beam Pulse Amplifier) mechanism for amplifying short noise electromagnetic pulses. The kappa distribution was used as a model function of the electron velocity distribution in the magnetosphere. Calculations performed for this distribution have shown that the threshold for excitation of chorus largely depends on the electron plasma parameter equal to the ratio of gas-kinetic pressure of electrons to magnetic pressure. This pattern is not contradicted by the dependence of the probability of excitation of chorus on the degree of magnetic field irregularity, which we derived from observations made by the Van Allen Probe spacecraft. It is sharp fluctuations in the magnetic field strength near its local minima outside the plasmasphere where the radiation under study can be excited. If there is an irregularity, the probability of detecting chorus is >70 %; and if there is no or very low irregularity, the probability of the absence of any emissions is ~80 %. The results indicate a common reason for the excitation of chorus and the magnetic field irregularity — a small but finite value of the plasma parameter.

ThunderBoltz: an open-source direct simulation Monte Carlo Boltzmann solver for plasma transport, chemical kinetics, and 0D modeling

Abstract Plasma-neutral interactions, including reactive kinetics, are often either studied in 0D using ODE-based descriptions, or in multi-dimensional fluid or particle-based plasma codes. The latter case involves a complex assembly of procedures that are not always necessary to test effects of underlying physical models and mechanisms for particle-based descriptions. Here we present ThunderBoltz, a lightweight, publicly available 0D direct simulation Monte Carlo code designed to accommodate a generalized combination of species and arbitrary cross sections without the overhead of expensive field solves. It can produce electron, ion, and neutral velocity distributions in applied AC/DC E-field and/or static B-field scenarios. The code is built in the C++ standard library and includes a convenient Python interface that allows for input file generation (including compatibility with cross section data from the LXCat database), electron transport and reaction rate constant post-processing, input parameter constraint satisfaction, calculation scheduling, and diagnostic plotting. These codes can be accessed at the repository: https://github.com/lanl/ThunderBoltz. In this work we compare ThunderBoltz transport calculations against Bolsig+ calculations, benchmark test problems, and swarm experiment data, finding good agreement with all three in the appropriate field regimes. In addition, we present example use cases where the electron, ion, and background neutral particle species are self-consistently evolved to model the background kinetics, a feature that is absent in fixed-background Monte Carlo and n-term Boltzmann solvers. The latter functionality allows for the possibility of particle-based chemical kinetics simulations of the plasma and neutral species as a new alternative to ODE-based approaches.

Particle Simulation Study of Obliquely Propagating Whistler Waves in Low‐Beta Plasmas

AbstractWhistler mode waves, which are electromagnetic emissions commonly observed in various space plasma environments, play critical roles in electron dynamics. While most whistler waves are driven by temperature anisotropy and propagate parallel to the background magnetic field, these waves can also be excited in the oblique direction when electron plasma beta is very low (<0.025). Although linear theory accounts for the excitation processes of obliquely propagating whistler waves, the subsequent evolution and saturation processes remain inadequately understood. This study utilizes two‐dimensional self‐consistent simulations to investigate the complete wave‐particle interaction process. By scanning a broad parameter range, we derive scaling laws for the wave intensity, linking the wave properties to initial and final plasma conditions. The saturated temperature anisotropy from simulations can explain the upper bound anisotropy constraint observed by spacecraft well. Additional phase space analysis shows that both Landau and cyclotron resonance play critical roles in the evolution of electron velocity distribution, albeit in different energy ranges. Oblique whistler waves can effectively heat electrons through Landau resonance, creating a plateau distribution in the parallel direction. This research advances our understanding of the mechanisms behind obliquely propagating whistler waves in low‐beta plasmas and their impact on electron dynamics.

Analysis of plasma wakefield generation by a Laguerre-Gaussian laser beam in laser-plasma accelerators with external injection

This study presents a comprehensive modeling of wakefield generation through external injection utilizing a Laguerre–Gaussian (LG) laser beam in a bubble/blowout regime. The wakefield dynamics are simulated in two dimensions using the particle-in-cell (2D-PIC) method via Wake-T tool, aiming to investigate the underlying mechanisms and characteristics of this process. The simulation results provide insights into the behavior of electrons within the wakefield, their acceleration, phase spaces of the electron beam, velocity distribution, and longitudinal and transversal profiles of the laser electric field in the plasma. The presented model serves as a valuable tool for further investigations into wakefield generation with external injection using LG laser beams, facilitating advancements in this field of study.

Nonlinear evolution of driven ion acoustic waves in plasmas with super-Gaussian electron velocity distributions: Fluid vs particle-in-cell simulations

Nonlinear evolution of driven ion acoustic waves in plasmas with super-Gaussian electron velocity distributions: Fluid vs particle-in-cell simulations

Kinetic and thermodynamic approach to precisely solve the unsteady Rayleigh flow problem of a rarefied homogeneous charged gas under external force influence

Abstract An extension and further development of our previous article [J. Non-equilibrium Thermodyne. 37 (2012), 119–141] is presented. We study the irreversible non-equilibrium thermodynamics (INT) properties of the exact solution to the dilute homogeneously charged gas problem with unsteady Rayleigh flow. In contrast to previous research, the charged gas flows under the influence of an external force, the flat plate oscillates, and the displacement current term is considered, leading to significant advancements in understanding natural plasma dynamics. We are solving the Boltzmann kinetic equation (BKE) Krook model supplemented by Maxwell’s equations. We used a travelling wave and moments method with an electron velocity distribution function (EVDF). To the best of our knowledge, as three new scientific achievements, we introduced a new mathematical model for calculating the thermodynamic forces, kinetic coefficients, and fluxes variables, Equations (28–40) and (50–54). Second, we determined, with reasonable accuracy, the thermodynamic equilibrium time of electrons, t equ = 26.7955, under an external force. We clarify the difference between equilibrium EVDF and perturbed EVDF and take advantage of BKE to account for non-equilibrium thermodynamic principles. For diamagnetic and paramagnetic plasmas, the extended Gibbs equation predicts ratios between various contributions to the internal energy change (IEC) is presented. A standard laboratory argon plasma model is used to apply the results.

Role of spontaneous thermal emissions in inflationary laser Raman instability

The role of thermal fluctuations on the stimulated Raman backscattering instability is investigated by means of Vlasov and particle-in-cell (PIC) simulations in a regime of strong linear Landau damping of the Langmuir wave. The instability is initially convective and amplifies thermal noise, leading to a low-amplitude back-scattered laser sideband. Linear Landau damping of the Langmuir sideband modifies and flattens the electron velocity distribution function at the resonant velocity, leading to a gradual decrease in the Landau damping rate and an increase in the convective amplification. The Langmuir wave traps electrons resulting in a rapid nonlinear absolute instability and large amplitude flashes of backscattered light off large amplitude Langmuir waves with trapped electrons, leading to the production of hot electrons. Conditions for simulating realistic thermal noise with Vlasov and PIC simulations are discussed and defined.

Properties of collisional plasma sheath with ionization source term and two-temperature electrons in an oblique magnetic field

A collisional magnetized plasma sheath with two groups of electrons has been studied using a fluid model including the effects of the ionization source term and the collisional force between ions and neutral atoms. Two kinds of non-Maxwellian descriptions of electron velocity distribution, non-extensive distribution and truncated distribution, are applied in the model, and the ionization effects of both kinds are considered. By applying Sagdeev potential, the modified Bohm sheath criterion is derived. The effects of ionization, magnetic field, and high-temperature electron concentration ratio on plasma sheath density, potential, sheath thickness, and ion kinetic energy are studied. In cases with high background gas density, ion density accumulates at the sheath edge position, forming a peak and manifesting as a rapid drop in the potential profile. The distribution characteristics of electrons have a significant impact on the transport properties of ions. Oscillations and non-monotonic characteristics of net charge near the sheath edge occur as the magnetic field angle increases, leading to an increase in the sheath layer width. It can be seen that in the case of a collisional sheath structure with high-temperature electrons, it is essential to consider the sheath changes induced by the ionization and the collisional force. Compared to a symmetric electron velocity distribution, the actual thickness of the sheath layer in a truncated electron distribution assumption could be significantly reduced.

Intermittency, bursty turbulence, and ion and electron phase-space holes formation in collisionless current-carrying plasmas

In the previous studies of nonlinear saturation of the Buneman instability caused by high electron drift velocity relative to ions, the phase-space holes and the plateau on the electron velocity distribution function were identified as features of the saturation stage of instability [notably in the paper by Omura et al., J. Geophys. Res. 108, 1197 (2003)]. We have performed a much longer simulation of the Buneman instability and observed a secondary instability. This secondary instability generates fast electron-acoustic waves. By analyzing the phase-space plot of ions and electrons, we show that the fast electron heating and the formation of the plateau of electron velocity distribution function are not due to the quasi-linear diffusion but due to the nonlinear interaction of ion- and electron-acoustic solitary waves (phase-space holes) by exchange of trapped electrons in each wave. We also report the details on the intermittent and bursty nature of turbulence driven by this instability.

Non-Thermal Solar Wind Electron Velocity Distribution Function.

The quiet-time solar wind electrons feature non-thermal characteristics when viewed from the perspective of their velocity distribution functions. They typically have an appearance of being composed of a denser thermal "core" population plus a tenuous energetic "halo" population. At first, such a feature was empirically fitted with the kappa velocity space distribution function, but ever since the ground-breaking work by Tsallis, the space physics community has embraced the potential implication of the kappa distribution as reflecting the non-extensive nature of the space plasma. From the viewpoint of microscopic plasma theory, the formation of the non-thermal electron velocity distribution function can be interpreted in terms of the plasma being in a state of turbulent quasi-equilibrium. Such a finding brings forth the possible existence of a profound inter-relationship between the non-extensive statistical state and the turbulent quasi-equilibrium state. The present paper further develops the idea of solar wind electrons being in the turbulent equilibrium, but, unlike the previous model, which involves the electrostatic turbulence near the plasma oscillation frequency (i.e., Langmuir turbulence), the present paper considers the impact of transverse electromagnetic turbulence, particularly, the turbulence in the whistler-mode frequency range. It is found that the coupling of spontaneously emitted thermal fluctuations and the background turbulence leads to the formation of a non-thermal electron velocity distribution function of the type observed in the solar wind during quiet times. This demonstrates that the whistler-range turbulence represents an alternative mechanism for producing the kappa-like non-thermal distribution, especially close to the Sun and in the near-Earth space environment.

Suprathermal Electron Transport in the Solar Wind: Effects of Coulomb Collisions and Whistler Turbulence

The nature and radial evolution of solar wind electrons in the suprathermal energy range are studied. A wave–particle interaction tensor and a Fokker–Planck Coulomb collision operator are introduced into the kinetic transport equation describing electron collisions and resonant interactions with whistler waves. The diffusion tensor includes diagonal and off-diagonal terms, and the Coulomb collision operator applies to arbitrary electron velocities describing collisions with both background protons and electrons. The background proton and electron densities and temperatures are based on previous turbulence models that mediate the supersonic solar wind. The electron velocity distribution functions and electron heat flux are calculated. Comparison and analysis of the numerical results with analytical solutions and observations in the near-Sun region are made. The numerical results reproduce well the creation of the sunward electron deficit observed in the near-Sun region. The deficit of the electron velocity distribution function below the core Maxwellian fit at low velocities results from Coulomb collisions, and the excess part above the core Maxwellian fit at high velocities is determined by strong wave–particle interactions.

The Generation of 150 km Echoes Through Nonlinear Wave Mode Coupling

AbstractA fundamental problem in plasma turbulence is understanding how energy cascades across multiple scales. In this paper, a new weak turbulence theory is developed to explain how energy can be transferred from Langmuir and Upper‐Hybrid waves to ion‐acoustic waves. A kinetic approach is used where the Boltzmann equation is Fourier‐Laplace transformed, and the nonlinear term is retained. A unique feature of this approach is the ability to calculate power spectra at low frequencies, for any wavelength or magnetic aspect angle. The results of this theory explain how the predominant type of 150‐km radar echoes are generated in the ionosphere. First, peaks in the suprathermal electron velocity distribution drive a bump‐on‐tail like instability that excites the Upper‐Hybrid mode. This excited wave then couples nonlinearly to the ion‐acoustic mode, generating the ∼10 dB enhancement observed by radars. This theory also explains why higher frequency radars like ALTAIR do not observe these echoes.

Loss cone effects and monotonic sheath conditions of a partially magnetized plasma sheath

In this Letter, we propose the conditions for monotonic plasma sheaths adjacent to a floating wall in the presence of an applied, oblique magnetic field. The electron velocity distribution function (VDF) at the sheath edge obtained from a kinetic model exhibits a loss cone shaped truncation. Using an approximation of the truncated VDF, we derive an analytical framework of the sheath edge condition (i.e., Bohm condition), namely, the relation of ion injection velocity and sheath potential drop as a function of the magnetic field angle. The results show that the sheath edge velocity and total potential drop decrease for a steady-state sheath, which eventually collapses when the magnetic field lines become parallel to the wall.

Solar Energetic Particle Charge States and Abundances with Nonthermal Electrons

An important aspect of solar energetic particle (SEP) events is their source populations. Elemental abundance enhancements of impulsive SEP events, originating in presumed coronal reconnection episodes, can be fitted to steep power laws of A/Q, where A and Q are the atomic mass and ionic charge. Since thermal electron energies are enhanced and nonthermal electron distributions arise in the reconnection process, we might expect that ionic charge states Q would be increased through ionization interactions with those electron populations during the acceleration process. The temperature estimated from the SEPs corresponds to the charge state during the acceleration process, while the actual charge state measured in situ may be modified as the SEPs pass through the corona. We examine whether the temperature estimation from the A/Q would differ with various κ values in a κ function representing high-energy tail deviating from a Maxwellian velocity distribution. We find that the differences in the A/Q between a Maxwellian and an extreme κ distribution are about 10%–30%. We fit power-law enhancement of element abundances as a function of their A/Q with various κ values. Then, we find that the derived source region temperature is not significantly affected by whether or not the electron velocity distribution deviates from a Maxwellian, i.e., thermal, distribution. Assuming that electrons are heated in the acceleration region, the agreement of the SEP charge state during acceleration with typical active region temperatures suggests that SEPs are accelerated and leave the acceleration region in a shorter time than the ionization timescale.

Study of the Electron Velocity Distribution Function in Weakly Ionized Radiofrequency Plasma

Study of the Electron Velocity Distribution Function in Weakly Ionized Radiofrequency Plasma

Excitation characterization of whistler waves and electrostatic oscillation in asymmetric magnetic reconnection

Utilizing the 2.5D Particle-in-Cell (PIC) method, we observe that large-amplitude whistler waves are present near the separatrix on the magnetospheric side and in the outflow region during asymmetric magnetic reconnection. Minimum variance analyses confirm their quasi-parallel and anti-parallel propagation to the background magnetic field, respectively, and their accuracy is validated through parallel Poynting fluxes (S∥). The excited whistler waves have intermittent property due to the evolution and movement of the reconnection layer. Wave-Particle Interaction Analysis (WPIA) reveals a pronounced energy conversion between the waves and electrons. Our simulation results demonstrate that whistler waves near the separatrix are excited by the cyclotron resonance duo to the instability of electron beam by analyzing the electron velocity distribution functions (VDFs). Meanwhile, the electrostatic oscillation is generated by the positive parallel drift of electron beam. However, the whistler waves in the outflow region are excited by the Landau resonance duo to anisotropic instability distributions of electron beam-mode.

On the possibility of mode coupling for low frequency electromagnetic waves in plasmas with anisotropy of temperature

We study the dispersion relation for low frequency electromagnetic waves propagating along the ambient magnetic field and investigate the possibility of occurrence of coupling between waves in the ion cyclotron branch and waves in the whistler branch. The results obtained show that the coupling may occur in the case of conditions leading to the ion cyclotron instability, for sufficiently high value of the ratio between perpendicular and parallel ion temperature, and does not occur in the case of conditions leading to the ion firehose instability. The results also show that the decrease in the value of the plasma beta may lead to the disappearance of the mode coupling conditions. Regarding the effect of the electron population, it is shown that the change in the shape of the electron velocity distribution, from Maxwellian to bi-Kappa form, does not change the results obtained, as long as the electron temperatures are isotropic, but the increase in anisotropy in the electron temperatures may lead to the disappearance of the coupling between the different waves. The consequences of the frequency dependency of the mode coupling conditions are discussed considering wave propagation in an inhomogeneous medium, leading to the conclusion that the energy of a packet of waves of a given mode can be absorbed or mode converted over an extended region of space. These findings can be of relevance for the analysis and understanding of processes related to the conversion between ion cyclotron waves and whistler waves.