Search

A simple analytical model of a hollow cathode effect in a direct-current (DC) discharge with parallel-plate cathode configuration is presented. The model is based on the empirical formulation of non-local ionization rate produced by fast electrons within the hollow cathode and takes into account the loss of ions due to ambipolar diffusion. The model reproduces a transition from a DC discharge with two separate near-cathode regions to a single hollow-cathode type discharge with the accompanying dramatic increase in current density as intercathode separation L decreases. It is shown that for typical conditions of a planar hollow-cathode discharge, ambipolar losses are relevant and contribute to deviation from pL scaling laws. Comparison with classical experiments show good agreement up until the smallest intercathode separations. Application of the model for assessment of experimental data on properties of discharges with parallel-plate cathode configuration and for studies of secondary electron emission from cathode surfaces is discussed. [ABSTRACT FROM AUTHOR]

The main objective of the present paper is the detailed analysis of the electron kinetics near the Paschen's curve minimum of pulsed breakdown of nitrogen gas. Three main questions are to be answered. First, whether the breakdown curve minimum corresponds to the threshold electric field necessary for the runaway electron generation. Second, what is the role of these electrons in the vicinity of the Paschen's curve minimum. Third, what is the ionization cost near the minimum. To answer these questions, the one-dimensional model is used in which electrons are modeled using the Particle-in-Cell approach, while ions are modeled in the drift-diffusion approximation. [ABSTRACT FROM AUTHOR]

The conditions for the generation of runaway electrons in an air gap are compared at different degrees of inhomogeneity of the electric field distribution provided by varying the opening angle of the conical cathode: in the range 40 ° – 120 ° in experiments and 0 ° – 180 ° in calculations. It is demonstrated that, in a weakly inhomogeneous electric field (according to the proposed classification, this corresponds to cones with angles greater than the Taylor angle of 98.6 °), the runaway condition has a local character. The transition of free electrons into the runaway mode is determined by the local distribution of the electric field near their starting point—the tip of the cone. The local electric field strength must exceed a threshold value comparable to the strength critical for the runaway of electrons in a uniform field. In a strongly inhomogeneous field (cones with angles less than 98.6 °), this condition is not sufficient for electrons to run away throughout the gap. Electrons accelerating in the near-cathode region may begin to slow down in a weak field at a distance from the cathode. In this case, the runaway condition becomes nonlocal. It is determined by the dynamics of electrons in the entire gap, primarily in the near-anode region, and reduces to the requirement that the potential difference applied to the gap exceeds a certain threshold value. [ABSTRACT FROM AUTHOR]

A model is presented for simulating the electron beam within the entire gyrotron geometry using an electron optics code. The model incorporates the influence of the electromagnetic fields generated by the cavity's operational mode through a simple self-consistent approach. This simulation method proves highly beneficial for investigating the effects occurring within the complete gyrotron structure. One such application includes the examination of ions and electrons, resulting from the ionization of residual gas molecules upon collision with the electron beam. By calculating the trajectories of these particles, interesting conclusions can be drawn. [ABSTRACT FROM AUTHOR]

We investigate the electron acceleration dynamics in spatially stochastic wave fields. For a bounded system, the spectra (frequency and wave number) of the stochastic wave fields are discrete so that they can form spatiotemporal "singular" structures once their phases are in synchronization. As a consequence, the electrons will experience significant scattering by these singular structures, which is an in-phase effect of an ensemble of non-resonant wave–particle interactions. In the presence of parallel symmetry breaking, it is found that the governing Fokker–Planck equation has a structure similar to that of the resonant wave–particle interaction but with a broader parameter regime in the velocity/energy space of the electron distribution function. [ABSTRACT FROM AUTHOR]

Effect of kinetic electrons on negative triangularity plasmas has been investigated and compared against the corresponding positive triangularity plasmas, using the global gyrokinetic code X-point Gyrokinetic Code with scale-separated delta-f option without Coulomb collisions. Our model magnetic equilibria have strong positive and negative triangularities and weak magnetic shear. However, unusually large ρ i / a and low density plasmas are chosen to maximize the nonlocal effect to investigate the finite ρi effect and to be clearly away from kinetic ballooning modes. Similar conclusions to previous flux tube and global simulations have been obtained in this highly nonlocal model plasma: it is essential to include kinetic electrons in the micro-instability study of negative triangularity plasmas. Most physics findings agree with existing reports, with some disagreement. We offer a new "effective trapping fraction" concept that can add to the explanation of the growth rate difference between NT and PT plasmas, pointing to the significant variation in trapped particle fractions that have turning points in the mode growth regions. [ABSTRACT FROM AUTHOR]

This work demonstrates the feasibility of creating a sub-millimeter, subsonic nitrogen gas jet using a 178-μm diameter orifice to conduct laser wakefield acceleration (LWFA) with 1-TW, 40-fs laser pulses. More importantly, our findings reveal that using a blade to impede part of the gas flow and create an asymmetric density profile with a shortened down-ramp leads to a notable reduction in pointing fluctuations and an increase in the total charge of the output electron beams. As evidenced by the corresponding particle-in-cell simulation, the laser intensity is more effectively sustained toward the downstream end of the shaped gas jet, allowing for effective excitation of low-amplitude plasma waves that help preserve the accelerated electrons over the target rear side. In contrast, the pulse intensity drops significantly within the rear side of the unshaped gas jet, resulting in continuously diminishing plasma waves and decreased beam charge. The steeper gradient of the density down-ramp in the shaped gas jet also leads to a more rapid increase in the plasma wavelength over a reduced propagation distance, which helps mitigate the dephasing of accelerated electrons and increase the charge at the high-energy side of the spectrum. Our study paves the way for the future development of few-TW LWFA using a subsonic gas jet with sharp edges to further enhance the properties of output electron beams. [ABSTRACT FROM AUTHOR]

The kinetics of electrons and negative ions was numerically studied under the conditions, which are typical for multi-pulsed guided ionization waves in He flows ejected into ambient air. It was shown that, in He with a small admixture of O2 or air, O2− ions formed due to three-body electron attachment in the discharge afterglow are rapidly converted to O4− cluster ions in pulse off time. We suggested that seed electrons are generated at the leading edge of a pulse in two steps, and these steps are dissociation of O4− ions to form O2− ions followed by electron detachment from the O2− ions. The rate of O4− dissociation was calculated in He–O2 mixtures for various reduced electric fields on the basis of a Monte Carlo simulation. This rate, as well as the rate of electron detachment from O2− ions calculated previously, was used to simulate seed electron generation in a high-voltage nanosecond pulse in He–O2 mixtures for different gas pressures. It was shown that the formation of O4− ions drastically hinders the generation of seed electrons in multi-pulsed guided ionization waves in He–O2 and He–air mixtures. [ABSTRACT FROM AUTHOR]

Hollow cathodes for plasma switch applications are investigated via 2D3V particle-in-cell simulations of the channel and plume region. The kinetic nature of the plasma within the channel is dependent on the thermalization rate of electrons, emitted from the insert. When Coulomb collisions occur at a much greater rate than ionization or excitation collisions, the electron energy distribution function rapidly relaxes to a Maxwellian and the plasma within the channel can be described accurately via a fluid model. In contrast, if inelastic processes are much faster than Coulomb collisions, then the electron energy distribution function in the channel exhibits a notable high-energy tail, and a kinetic treatment is required. This criterion is applied to hollow cathodes from the literature, revealing that a fluid approach is suitable for most electric propulsion applications, whereas a kinetic treatment can be more critical to accurate modeling of plasma switches. [ABSTRACT FROM AUTHOR]

Collisionless absorption of laser energy by an electron via laser–cluster interaction in an ambient magnetic field (B0) has recently renewed interest. Previously, using a rigid sphere model (RSM) and an extensive particle-in-cell (PIC) simulation with linearly polarized (LP) laser light, we have shown that an auxiliary field B0 in a transverse direction to the laser polarization significantly enhances the laser absorption [Swain et al., Sci. Rep. 12, 11256 (2022)]. In this LP case, the average energy ( E A ) of an electron rises nearly 30 − 70 times of its ponderomotive energy (Up). The two-stage laser absorption by cluster electrons has been attributed via anharmonic resonance (AHR) followed by electron-cyclotron resonance (ECR) satisfying the improved phase-matching and frequency-matching conditions simultaneously. In the present work, we study the effect of circularly polarized (CP) laser fields on the cluster-electron dynamics considering left/right circular polarizations with an ambient B0. In typical conditions, without B0, we show that both LP and CP light yield almost the same level of laser absorption (about 3 U p or less) by an electron. However, with B0, CP light enhances the electron's energy further by ≈ 10 − 20 U p beyond the previously reported values ≈ 30 − 70 U p by the LP light. These ejected electrons from cluster show narrow cone-like propagation as a weakly relativistic electron beam with an angular spread Δ θ < 5 ° and the beam quality improves in CP than LP. In all cases, RSM and PIC results show good agreement. [ABSTRACT FROM AUTHOR]

Phase-space holes are well-known Bernstein–Greene–Kruskal (B.G.K.) modes and are formed by particle-trapping in solitary potential waveforms. They exhibit orbital particle trajectories in the phase-space, due to which they are also referred to as phase-space vortices. In this article, we develop the theory of phase-space hydrodynamics for electron and ion phase-space in collisionless plasmas. The analogy between ordinary two-dimensional fluids and 1 D − 1 V phase-space has been explored by introducing a momentum equation and a phase-space vorticity field, which enable the fluid-like analyses of the plasma phase-space. The developed kinetic-hydrodynamic equations are then employed to address the vortical nature of phase-space holes by exploring their fluid-analogous vortex-like characteristics, an identification technique of phase-space vortices, an exact derivation of the Schamel-df equations, and a measurable definition of the particle-trapping β parameter. This article introduces a new technique to the study of phase-space holes which focuses on the fluid-analogous vortical nature of the phase-space holes and prevents the need for an initial assumption of the trapped and free particle phase-space densities, thus presenting itself as a precursor to the Schamel-pseudopotential method. [ABSTRACT FROM AUTHOR]

This short report describes the phenomenon of space-charge-limited (SCL) current transport between ballistic and collision-dominated regimes for an inhomogeneous medium with a finite emission area. This intermediate regime can be considered as a weakly collisional regime where the collisional mean free path is comparable to the length of the medium. The SCL current density is calculated as a function of the degree of collision, inhomogeneity of the medium, and the geometrical properties of the emitting area. The inhomogeneity of the medium is characterized by a parameter of (0 < α ≤ 1) , where α = 1 denotes a perfect homogenous medium. The calculated SCL current density is enhanced by finite emission area effects by a factor of 1 + F × G , where F measures the mean position of the electrons in the medium, and G is a geometrical correction factor due to finite emission area. The enhancement is found to be higher in the collisional regime as compared to the ballistic regime. A higher inhomogeneity (with smaller α) also increases the enhancement. Smooth transition between the fully ballistic SCL transport (Child–Langmuir model) and the collision-dominated SCL transport (Mott–Gurney model) is demonstrated and verified, respectively, by the particle-in-cell code and the device simulator. This model is useful for the characterization of high-current SCL transport where the non-ideal conditions (such as inhomogenous medium and weakly collisional regime) cannot be described by the existing SCL models. [ABSTRACT FROM AUTHOR]

We present details of recent extensions of the runaway electron (RE) fluid model implemented in the fusion magnetohydrodynamics code JOREK [M. Hoelzl et al., Nucl. Fusion 61, 065001 (2021)] to include the effects of partially ionized impurity species and deuterium neutrals. The model treats the interaction of runaway electrons with the background plasma via current-coupling. The code is separately benchmarked using ITER (https://www.iter.org/) relevant scenarios, with the GO [G. Papp et al., Nucl. Fusion 53, 123017 (2013)] code in relation to runaway electron beam formation and with the DINA [Khayrutdinov and Lukash, J. Comp. Phys. 109(2), 193–201 (1993)] code in relation to simultaneous runaway beam formation and vertical plasma motion. Benchmark results show a decent agreement in both the cases, which are also discussed. [ABSTRACT FROM AUTHOR]

Quasi-parallel magnetized collisionless shocks are believed to be one of the most efficient accelerators in the universe. Compared to quasi-perpendicular shocks, quasi-parallel shocks are more difficult to form in the laboratory and to simulate because of their large spatial scales and long formation times. Our two-dimensional particle-in-cell simulations show that the early stages of quasi-parallel shock formation are achievable in experiments planned for the National Ignition Facility and that particles accelerated by diffusive shock acceleration (DSA) are expected to be observable in the experiment. Repetitive ion acceleration by crossings of the shock front, a key feature of DSA, is seen in the simulations. Other characteristic features of quasi-parallel shocks such as upstream wave excitation by energetic ions are also observed, and energy partition between the ions and the electrons in the downstream of the shock is briefly discussed. [ABSTRACT FROM AUTHOR]

Plasma-based acceleration (PBA) is being considered for a next generation linear collider (LC). In some PBA-LC designs for the electron arm, the extreme beam parameters are expected to trigger background ion motion within the witness beam, which can lead to longitudinally varying nonlinear focusing forces and result in an unacceptable emittance growth of the beam. To mitigate this, we propose to use quasi-adiabatic plasma density ramps as matching sections at the entrance and exit of each stage. We match the witness electron beam to the low density plasma entrance, where the beam initially has a large matched spot size so the ion motion effects are relatively small. As the beam propagates in the plasma density upramp, it is quasi-adiabatically focused, and its distribution maintains a non-Gaussian equilibrium distribution in each longitudinal slice throughout the process, even when severe ion collapse has occurred. This only causes small amounts of slice emittance growth. The phase mixing between slices with different betatron frequencies leads to additional projected emittance growth within the acceleration stage. A density downramp at the exit of an acceleration section can eliminate much of the slice and projected emittance growth as the beam and ion motion adiabatically defocuses and decreases, respectively. Simulation results from QuickPIC with Azimuthal Decomposition show that within a single acceleration stage with a 25 GeV energy gain, this concept can limit the projected emittance growth to only ∼2% for a 25 GeV, 100 nm emittance witness beam and ∼20% for a 100 GeV, 100 nm normalized emittance witness beam. The trade-off between the adiabaticity of the plasma density ramp and the initial ion motion at the entrance for a given length of the plasma density ramp is also discussed. [ABSTRACT FROM AUTHOR]

A set of equations is developed that extends the macroscale magnetic reconnection simulation model kglobal to include particle ions. The extension from earlier versions of kglobal, which included only particle electrons, requires the inclusion of the inertia of particle ions in the fluid momentum equation. The new equations will facilitate the exploration of the simultaneous non-thermal energization of ions and electrons during magnetic reconnection in macroscale systems. Numerical tests of the propagation of Alfvén waves and the growth of firehose modes in a plasma with anisotropic electron and ion pressure are presented to benchmark the new model. [ABSTRACT FROM AUTHOR]

The exponential growth (avalanching) of runaway electrons (REs) during a tokamak disruption continues to be a large uncertainty in RE modeling. The present work investigates the impact of tokamak geometry on the efficiency of the avalanche mechanism across a broad range of disruption scenarios. It is found that the parameter ν * , crit , describing the collisionality at the critical energy to run away, delineates how toroidal geometry impacts RE formation. In particular, utilizing a reduced but self-consistent description of plasma power balance, it is shown that for a high-density deuterium-dominated plasma, ν * , crit is robustly less than one, resulting in a substantial decrease in the efficiency of the RE avalanche compared to predictions from slab geometry. In contrast, for plasmas containing a substantial quantity of neon or argon, ν * , crit ≳ 1 , no reduction of the avalanche is observed due to toroidal geometry. This sharp contrast in the impact of low- vs high-Z material results primarily from the relatively strong radiative cooling from high-Z impurities, enabling the plasma to be radiatively pinned at low temperatures and thus large electric fields, even for modest quantities of high-Z material. [ABSTRACT FROM AUTHOR]

Microturbulence in magnetic confined plasmas contributes to energy exchange between particles of different species as well as the particle and heat fluxes. Although the effect of turbulent energy exchange has not been considered significant in previous studies, it is anticipated to have a greater impact than collisional energy exchange in low collisional plasmas such as those in future fusion reactors. In this study, gyrokinetic simulations are performed to evaluate the energy exchange due to ion temperature gradient (ITG) turbulence in a tokamak configuration. The energy exchange due to the ITG turbulence mainly consists of the cooling of ions in the ∇ B -curvature drift motion and the heating of electrons streaming along a field line. It is found that the ITG turbulence transfers energy from ions to electrons regardless of whether the ions or electrons are hotter, which is in marked contrast to the energy transfer by Coulomb collisions. This implies that the ITG turbulence should be suppressed from the viewpoint of sustaining the high ion temperature required for fusion reactions since it prevents energy transfer from alpha-heated electrons to ions as well as enhancing ion heat transport toward the outside of the reactor. Furthermore, linear and nonlinear simulation analyses confirm the feasibility of quasilinear modeling for predicting the turbulent energy exchange in addition to the particle and heat fluxes. [ABSTRACT FROM AUTHOR]

Understanding laser–target coupling is of the utmost importance for achieving high performance in laser-direct-drive (LDD) inertial confinement fusion (ICF) experiments. Thus, accurate modeling of electron transport and deposition through ICF-relevant materials and conditions is necessary to quantify the total thermal conduction and ablation. The stopping range is a key transport quantity used in thermal conduction models; in this work, we review the overall role that the electron mean free path (MFP) plays in thermal conduction and hydrodynamic simulations. The currently used modified Lee–More model employs various physics approximations. We discuss a recent model that uses time-dependent density functional theory (TD-DFT) to eliminate these approximations in both the calculation of the electron stopping power and corresponding MFP in conduction zone polystyrene (CH) plasma. In general, the TD-DFT calculations showed a larger MFP (lower stopping power) than the standard modified Lee–More model. Using the TD-DFT results, an analytical model for the electron deposition range, λ T D − DFT (ρ , T , K) , was devised for CH plasmas between ρ = [ 0.05 − 1.05 ] g / cm 3 , k B T = [ 100 − 1000 ] eV. We implemented this model into LILAC, for simulations of a National Ignition Facility-scale LDD implosion and compared key physics quantities to ones obtained by simulations using the standard model. The implications of the obtained results and the path moving forward to calculate this same quantity in conduction-zone deuterium–tritium plasmas are further discussed, to hopefully close the understanding gap for laser target coupling in LDD-ICF simulations. [ABSTRACT FROM AUTHOR]

This work characterizes the core transport physics of SPARC early-campaign plasmas using the PORTALS-CGYRO framework. Empirical modeling of SPARC plasmas with L-mode confinement indicates an ample window of breakeven (Q > 1) without the need of H-mode operation. Extensive modeling of multi-channel (electron energy, ion energy, and electron particle) flux-matched conditions with the nonlinear CGYRO code for turbulent transport coupled to the macroscopic plasma evolution using PORTALS reveals that the maximum fusion performance to be attained will be highly dependent on the near-edge pressure. Stiff core transport conditions are found, particularly when fusion gain approaches unity, and predicted density peaking is found to be in line with empirical databases of particle source-free H-modes. Impurity optimization is identified as a potential avenue to increase fusion performance while enabling core-edge integration. Extensive validation of the quasilinear TGLF model builds confidence in reduced-model predictions. The implications of projecting L-mode performance to high-performance and burning-plasma devices is discussed, together with the importance of predicting edge conditions. [ABSTRACT FROM AUTHOR]

Disorder-induced heating (DIH) prevents ultracold neutral plasma into the electron strong coupling regime. Here, we propose a scheme to suppress electronic DIH via optical lattice. We simulate the evolution dynamics of ultracold neutral plasma constrained by three-dimensional optical lattice using the classical molecular dynamics method. The results show that for experimentally achievable conditions, electronic DIH is suppressed by a factor of 1.3, and the Coulomb coupling strength can reach 0.8, which is approaching the strong coupling regime. Suppressing electronic DIH via optical lattice may pave a way for the research of electronic strongly coupled plasma. [ABSTRACT FROM AUTHOR]

A novel method of spatial separation of electrons with different initial energies in the helical electron beams of gyro-traveling wave tube (gyro-TWT), configured by two orthogonal external transverse magnetic fields, is presented. To demonstrate the method, an advanced asymmetric two-stage depressed collector has been designed. The structure of the collector and the electronic classification method are analyzed in detail in this article. The simulation shows that the application of this method can realize electronic classification effectively. The overall efficiency of the gyro-TWT improved from 32.6% to 55.7% and increased by 70.8%. The maximum power dissipation and the maximum thermal loading on the electrode were 2.12 kW/cm2 and 299 °C, respectively. The secondary electron analysis shows that the rectangular groove on the electrode surface can effectively suppress the influence of secondary electrons on the collector's efficiency. [ABSTRACT FROM AUTHOR]

A theory for presheath-like structures near probes biased at and above the plasma potential is developed for collisionless plasmas with an electron-neutral mean free path on the order of the chamber scale. The theory predicts presheath-like perturbations to the plasma that result from the free streaming of electrons and an effusion loss process from the chamber at the electrode. For these situations, a loss-cone-like velocity distribution function for electrons is predicted where the loss angle of the depletion region corresponds to the angular size of the electrode at a specified distance. The angle of the loss cone becomes 180° at the sheath edge. In comparison with a previous collisional electron presheath model that required electrons satisfy a Bohm criterion at the sheath edge [Scheiner et al., Phys. Plasmas 22, 123520 (2015)], the present work suggests that no such condition is needed for collisionless low pressure plasmas in the ≲ 10 mTorr range. The theory predicts the generation of a density depletion of roughly 0.5ne and an electron velocity moment of tens of percent of the electron thermal speed by the sheath edge in a presheath with a potential drop of less than T i / e. The range of this presheath perturbation is determined by the electrode geometry instead of the collisional mean free path. These predictions are tested against previously published particle in cell simulations and are found to be in good agreement. [ABSTRACT FROM AUTHOR]

We perform a numerical study of the linear dynamics of tearing modes in slab incompressible electron-magnetohydrodynamics (EMHD) by considering some parameter ranges, which can be of interest for laboratory plasmas (e.g., helicon devices) or for astrophysics (e.g., solar-wind turbulence). To this purpose, several non-ideal effects are simultaneously retained (finite electron inertia, resistivity, and electron viscosity), and we make distinction between the dissipation coefficients in the direction parallel and perpendicular to the guide field. We thus identify some new reconnection regimes, characterized by a departure from the customary monotonic power-law scalings of the growth rates with respect to the non-ideal parameters. The results here presented can provide a useful indication for future studies of EMHD regimes relevant to experiments and for extensions of the EMHD tearing mode modeling to more complete regimes including kinetic effects (e.g., "electron-only" reconnection in kinetic regimes). [ABSTRACT FROM AUTHOR]

A one-dimensional model, in which the ionization and collision are omitted in the Debye sheath region, is applied to investigate the effect of the superthermal electrons on the heat flux through a magnetized sheath. For different temperatures and concentrations of the superthermal electron, and magnitudes and directions of the magnetic field, the profiles of plasma parameters including ion density and flow velocity perpendicular to the wall, the background and superthermal electron densities, and sheath potential in the presheath region are calculated. The variation of the plasma density and sheath potential drop at the Debye sheath entrance with the superthermal electrons and magnetic field modifies the particle and heat fluxes across the Debye sheath to the material surface. The sheath heat transmission coefficient can increase significantly even for a very small superthermal electron population. The dependence of the sheath heat transmission coefficient on the magnetic field angle decreases with the contribution of the superthermal electron in a strong magnetized sheath. When investigation of the heat flux including the superthermal electrons to a water-cooled W/Cu monoblack for the tokamak divertor, compared to the case of without superthermal electrons, it is found that the increase in both heat flux to the material surface and surface temperature of the material is mainly due to the enhancement of the sheath potential drop caused by the superthermal electrons, but the increase in the latter is not as pronounced as the former. [ABSTRACT FROM AUTHOR]

The problem of the trigger of the neoclassical tearing mode is addressed in this work by evaluating the non-linear terms of the generalized Rutherford equation (GRE) for a set of JET disruptions. The linear stability index Δ ′ 0 , even for positive values, is not enough to describe the trigger of the mode as the stabilizing non-linear effects tend to prevent the growth of a mode below a certain width. First, an analysis on the contribution of the stabilizing effect of the curvature and the destabilizing effect of the bootstrap is done. Second, the work focuses on the role of the ion polarization current, a return current that flows parallel to the magnetic field due to the difference in the drift motion of the electrons and ions. This contribution is thought to play an important role in the onset of an island of width W, because it scales as 1 / W 3 , making it a dominant term in the GRE when W is small. The assessment is carried out over a subset of pulses, producing interesting observations that were then generalized across the entire database, obtaining consistent results. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The propagation features of dust acoustic waves in a three-component plasma system composed of regularized Kappa distributed electrons, Maxwellian ions, and dust grains carrying positive charges are investigated. The reductive perturbation technique is employed to derive the KdV equation. A generalized expression for the polarization force is derived and the effect of the polarization force is taken into consideration as well. The bifurcation analysis is used, and the solitary wave solution was investigated. The critical value of the superthermal spectral index κ is introduced at which the solitonic structure turns up from rarefactive to compressive. It is found that in the range 0 < κ < 2.23 , a rarefactive structure is obtained while the compressive structure appears for κ > 2.23. In addition, it is found that by increasing the value of cutoff parameter α, the polarization strength increases too. All the obtained results are helpful to investigate the characteristics of the nonlinear wave propagating in the mesosphere region. [ABSTRACT FROM AUTHOR]

In negative ion sources for the ITER Neutral Beam Injection system, the co-extraction of electrons is one of the main limiting factors. The current of co-extracted electrons can be decreased by applying a positive bias voltage to the Plasma Grid (PG) with respect to its source walls. Simulations using three-dimensional Particle-in-Cell Monte Carlo Collision (3D-PIC MCC) model are a powerful tool for studying the extraction region of such ion sources. However, the inclusion of both PG and source walls in the simulation domain is difficult due to numerical constraints. This study uses the 3D-PIC MCC code ONIX to explore the effects of particle injection models on plasma characteristics, using a flux injection model to regulate particle influx for a flat transition in potential from the bulk plasma to the simulation domain. Biasing of the PG above floating potential is possible using the flux injection scheme and results in a notable reduction in co-extracted electrons, corroborating with established experimental observations. [ABSTRACT FROM AUTHOR]

The purpose of this article is to study turbulent transport for laboratory plasmas in toroidal devices by gyrokinetic analyses. Linear analysis is performed to clarify the dominant mode for tokamak plasmas. The dissipative trapped electron mode (d-TEM) and the ion temperature gradient (ITG) mode are predicted using the Sugama collision model operator [Sugama et al., Phys. Plasmas 16, 112503 (2009)]. Nonlinear gyrokinetic analysis is used to quantify turbulent transport. The nonlinear simulation results show the levels of particle and energy transport, where the d-TEM and ITG mode are unstable. The effect of zonal flows is studied by the linear and nonlinear simulation results. The results of the analysis are compared when two types of model collision operator, which are the Sugama and Lenard–Bernstein [Phys. Rev. 112, 1456 (1958)] collision model operators, are used. In this study, the simulation results using the Sugama collision operator show a stronger effect of the zonal flows on the turbulent transport than those using the Lenard–Bernstein collision operator, as predicted by the linear simulation result such as the zonal flow decay time. [ABSTRACT FROM AUTHOR]

Nonlinear kinetic structures, called electron phase-space holes (EHs), are regularly observed in space and experimental magnetized plasmas. The existence of EHs is conditioned and varies according to the ambient magnetic field and the parameters of the electron beam(s) that may generate them. The objective of this paper is to extend the 3D Bernstein–Greene–Kruskal model with cylindrical geometry developed by L.-J. Chen et al. ["Bernstein–Greene–Kruskal solitary waves in three-dimensional magnetized plasma," Phys. Rev. E 69, 055401 (2004)] and L.-J. Chen et al., ["On the width-amplitude inequality of electron phase space holes," J. Geophys. Res. 110, A09211 (2005)] to include simultaneously finite effects due to (i) the strength of the ambient magnetic field B 0 , by modifying the Poisson equation with a term derived from the electron polarization current, and (ii) the drift velocity ue of the background plasma electrons with respect to the EH, by considering velocity-shifted Maxwellian distributions for the boundary conditions. This allows us to more realistically determine the distributions of trapped and passing particles forming the EHs, as well as the width-amplitude relationships for their existence. [ABSTRACT FROM AUTHOR]

The whistlers in space plasmas and in magnetic fusion experiments are destabilized by beams of fast electrons. While the linear regime of instability is analytically tractable, in most practical cases, the instability operates at the saturated level during the stages of observation and measurement. The saturated states, however, involve nonlinear whistlers, which remain best accessible for analysis by kinetic simulations. Results of electromagnetic Vlasov simulations are presented, analyzing an anisotropic electron beam driven whistler instability. The simulations cover the initially unstable regime followed by a saturated or marginally stable regime. Both regimes are separated by an intermediate nonlinear regime during which the electron distribution undergoes a kinetically self-consistent modification. A linearly obtained generalized marginal stability condition is applied to the stabilized state. The condition obtained in its dispersive version shows the β | | at threshold and, in turn, the residual anisotropy, to be a function of the whistler mode number k. [ABSTRACT FROM AUTHOR]

This comment discusses the feasibility of hot ion mode T i / T e = 4 for proton–boron fusion, which is critical for the roadmap proposed in Liu et al. [Phys. Plasmas 31, 062507 (2024)]. The hot ion mode T i / T e = 4 has been calculated to be far from accessible ( T i / T e < 1.5 for T i = 150 keV) under the most optimal conditions if fusion provides the heating [Xie, Introduction to Fusion Ignition Principles: Zeroth Order Factors of Fusion Energy Research (USTC Press, Hefei, 2023)], i.e., that all fusion power serves to heat the ions and that electrons acquire energy only through interactions with ions. We also discuss if hot ion mode of T i / T e = 4 could be achieved by an ideal heating method, which is much more efficient than fusion itself (near 20 times fusion power for T i = 150 keV) and only heats the ions, whether it makes sense economically. [ABSTRACT FROM AUTHOR]

We present a sensitivity analysis of the multipactor susceptibility zones to variations in the secondary electron yield (SEY) of materials, specifically focusing on the first and second unity crossover points of SEYs. In conducting this research, we leveraged our semi-analytic approach, which allows for the rapid prediction of the full multipactor zones with enhanced accuracy. Using this approach, we unveil several unique features of multipactor susceptibility zones, including the infinite extension of different-order multipactor zones and the overlap between them. Building upon this prediction capability, our results complement previous findings on the same topic and reveal that the multipactor zones depend not only on the first crossover point but also on the second crossover point of SEY, which this latter predominantly impacts multipactor susceptibility zones for low SEY materials. To validate our predictions, we present two distinct sets of multipactor experiments, providing empirical support for our results. [ABSTRACT FROM AUTHOR]

Particle energization due to magnetic reconnection is an important unsolved problem for myriad space and astrophysical plasmas. Electron energization in magnetic reconnection has traditionally been examined from a particle, or Lagrangian, perspective using particle-in-cell (PIC) simulations. Guiding-center analyses of ensembles of PIC particles have suggested that Fermi (curvature drift) acceleration and direct acceleration via the reconnection electric field are the primary electron energization mechanisms. However, both PIC guiding-center ensemble analyses and spacecraft observations are performed in an Eulerian perspective. For this work, we employ the continuum Vlasov–Maxwell solver within the Gkeyll simulation framework to reexamine electron energization from a kinetic continuum, Eulerian, perspective. We separately examine the contribution of each drift energization component to determine the dominant electron energization mechanisms in a moderate guide-field Gkeyll reconnection simulation. In the Eulerian perspective, we find that the diamagnetic and agyrotropic drifts are the primary electron energization mechanisms away from the reconnection x-point, where direct acceleration dominates. We compare the Eulerian (Vlasov Gkeyll) results with the wisdom gained from Lagrangian (PIC) analyses. [ABSTRACT FROM AUTHOR]

Direct laser acceleration of electrons is an important energy deposition mechanism for laser-irradiated plasmas that is particularly effective at relativistic laser intensities in the presence of quasi-static laser-driven plasma electric and magnetic fields. These radial electric and azimuthal magnetic fields provide transverse electron confinement by inducing betatron oscillations of forward-moving electrons undergoing laser acceleration. Electrons are said to experience a betatron resonance when the frequency of betatron oscillations matches the average frequency of the laser field oscillations at the electron position. In this paper, we show that the modulation of the laser frequency as seen by an electron performing betatron oscillations can be another important mechanism for net energy gain that is qualitatively different from the betatron resonance. Specifically, we show that the frequency modulation experienced by the electron can lead to net energy gain in the regime where the laser field performs three oscillations per betatron oscillation. There is no net energy gain in this regime without the modulation because the energy gain is fully compensated by the energy loss. The modulation slows down the laser oscillation near transverse stopping points, increasing the time interval during which the electron gains energy and making it possible to achieve net energy gain. [ABSTRACT FROM AUTHOR]

We extend the bounce-averaged kinetic (BK) electron model to be applicable in general tokamak magnetic geometries and implement it on the global δf particle-in-cell gyrokinetic code gKPSP. We perform a benchmark study of the updated BK model against the gyrokinetic electron model in flux-tube codes, CGYRO and GENE. From the comparisons among the simulations based on the local parameters of a KSTAR L-mode plasma, we confirm a reasonable agreement among the linear results from the different codes. In the nonlinear gKPSP simulation with a narrow plasma gradient region whose width comparable to the mode correlation length, ion and electron heat fluxes are compatible with those calculated by CGYRO. However, with an unstable region sufficiently wider than the mode correlation length, gKPSP predicts 2–3 times larger turbulent heat fluxes. Taking into account the differences between the flux-tube and global simulations, the overall agreement is encouraging for further validation and development of the BK electron model. In global simulations using a wide range of the experimental plasma profiles, we find an intricate coupling of turbulence spreading and a zonal flow in determining the radial profiles of turbulent heat fluxes, which has not been reported to date. [ABSTRACT FROM AUTHOR]

The resonant field penetration to magnetic islands is the central MHD mechanism of non-axisymmetric plasma responses in a tokamak such as disruptive locking or favorable ELM stabilization. The resonant field penetration can be induced by any non-ideal processes as manifested in the delicate balance under the generalized Ohm's law. Here, we show that the viscous effects by electrons are not ignorable in the field penetration unlike previous presumption, even if the electron viscosity is as small as the square root of its mass compared to the ions. It is clear that its effects become only bigger if the electron viscosity becomes anomalously large. The work strictly follows the three-field model in the linear regime targeting the prediction of the onset of the field penetration and successfully extending it with electron viscosity and identifying new regimes. The results also indicate that the error field thresholds become more strongly dependent on plasma density than ones predicted in the linear regimes without the electron viscosity, which is consistent with experimental observations and thus a significant implication. [ABSTRACT FROM AUTHOR]

We consider the steady-state limiting current that can be carried by an infinite periodic array of thin electron sheets spaced by period p in a planar diode of gap voltage V and gap separation d. Our primary assumptions are as follows: (1) electron motion is restricted by an infinite magnetic field to the direction normal to the electrode surfaces, (2) all electrons are emitted from the cathode with initial kinetic energy Ein, and (3) electron motion is non-relativistic. The limiting current density, averaged over a period and normalized to the classical 1D Child–Langmuir (CL) current density (including a factor that accounts for nonzero Ein), is found to depend only on the two dimensionless parameters p/d and Ein/eV. This average limiting current density is computed from the maximum current density for which the iterative solution of a non-linear integral equation converges. Numerical results and empirical curve fits for the limiting current are presented, together with an analysis as p/d and Ein/eV approach zero or infinity, in which cases previously published results are recovered. Our main finding is that, while the local anode current density within each electron sheet is infinite in our model (that is, it exceeds the classical 1D CL value by an "infinite" factor), the period average anode current density is in fact still bounded by the classical 1D CL value. This study therefore provides further evidence that the classical 1D Child–Langmuir current density is truly a fundamental limit that cannot be circumvented. [ABSTRACT FROM AUTHOR]

Starting from the beginning of their research in the early 2000s, the ultracold plasmas were considered as a promising tool to achieve considerable values of the Coulomb coupling parameter for electrons. Unfortunately, this was found to be precluded by a sharp spontaneous increase in temperature, which was often attributed to the so-called disorder-induced heating. It is the aim of the present paper to quantify the effect of spontaneous heating as a function of the initial ionic disorder and, thereby, to estimate the efficiency of its mitigation, e.g., by the Rydberg blockade. As a result of the performed simulations, we found that the dynamics of electrons exhibited a well-expressed transition from the case of the quasi-regular arrangement of ions to the disordered one; the magnitude of the effect being about 30%. Thereby, we can conclude that the two-step formation of ultracold plasmas—involving the intermediate stage of the blockaded Rydberg gas—can really serve as a tool to increase the degree of Coulomb coupling, but the efficiency of this method is moderate. [ABSTRACT FROM AUTHOR]

We present a new analytical model for the expansion of a thin slab of plasma into vacuum. By considering the rising plasma scale length during the initial heating phase, we were able to give the plasma a smooth quadratic behavior at the origin while describing its exponentially falling density at a large distance. We show this functional form to be a solution to the plasma equations and validate its predictions against numerical simulations and experimental measurements. We demonstrate the applicability of the model to experimental scenarios in which solid foils are turned into tens of micrometer-scale plasmas, to serve as targets for direct laser acceleration of electrons. [ABSTRACT FROM AUTHOR]

Electron resonant scattering by high-frequency electromagnetic whistler-mode waves has been proposed as a mechanism for solar wind electron scattering and pre-acceleration to energies that enable them to participate in shock drift acceleration around the Earth's bow shock. However, observed whistler-mode waves are often sufficiently intense to resonate with electrons nonlinearly, which prohibits the application of quasi-linear diffusion theory. This is the second of two accompanying papers devoted to developing a new theoretical approach for quantifying the electron distribution evolution subject to multiple resonant interactions with intense whistler-mode wave-packets. In the first paper, we described a probabilistic approach, applicable to systems with short wave-packets. For such systems, nonlinear resonant effects can be treated by diffusion theory, but with diffusion rates different from those of quasi-linear diffusion. In this paper, we generalize this approach by merging it with a mapping technique. This technique can be used to model the electron distribution evolution in the presence of significantly non-diffusive resonant scattering by intense long wave-packets. We verify our technique by comparing its predictions with results from a numerical integration approach. [ABSTRACT FROM AUTHOR]

Adiabatic heating of solar wind electrons at the Earth's bow shock and its foreshock region produces transversely anisotropic hot electrons that, in turn, generate intense high-frequency whistler-mode waves. These waves are often detected by spacecraft as narrow-band, electromagnetic emissions in the frequency range of [0.1, 0.5] of the local electron gyrofrequency. Resonant interactions between these waves and electrons may cause electron acceleration and pitch-angle scattering, which can be important for creating the electron population that seeds shock drift acceleration. The high intensity and coherence of the observed whistler-mode waves prohibit the use of quasi-linear theory to describe their interaction with electrons. In this paper, we aim to develop a new theoretical approach to describe this interaction, which incorporates nonlinear resonant interactions, gradients of the background density and magnetic field, and the fine structure of the waveforms that usually consist of short, intense wave-packet trains. This is the first of two accompanying papers. It outlines a probabilistic approach to describe the wave–particle interaction. We demonstrate how the wave-packet size affects electron nonlinear resonance at the bow shock and foreshock regions, and how to evaluate electron distribution dynamics in such a system that is frequented by short, intense whistler-mode wave-packets. In the Paper II, this probabilistic approach is merged with a mapping technique, which allows us to model systems containing short and long wave-packets. [ABSTRACT FROM AUTHOR]

In this paper, the effect of appearance of fast electrons behind the grid cathode in the direction reverse to the anode is studied computationally. Fast electrons are observed within 0.5 ns after application of a short voltage pulse. The results obtained confirm the possibility of generation of backward fast electrons (some of them are in a runaway mode) and explain the main trends of this process. It is shown that backward fast electrons are supported by the ionization wave (IW). The IW evolution proceeds via two phases. During the first phase, first fast electrons are observed moving toward the anode. Then, multiple individual IWs starting from each individual cathode wire are formed in the anode–cathode gap. The duration of this stage is 0.3 ns and corresponds to the pulse rise time. At the second phase, the separate individual IWs merge in a single flat ionization wave. The IW penetrates through the cathode wires and propagates in the direction reverse to the anode. The preferred direction of fast electrons propagation also reverses. Now, the trajectories of fast electrons are mainly directed away from the anode. The duration of the recorded flux of fast electrons is of the order of a few picoseconds. This time interval correlates with the available experimental data. [ABSTRACT FROM AUTHOR]

We analyze and model fast-electron-induced Kα emission from an experiment in which a high-intensity ultra-short laser irradiated foil and bulk titanium targets. The motion of electrons inside the targets is calculated allowing for multiple scattering and collisional energy loss, while outside the target, electric fields of arbitrary configurations are assumed. It is shown that both the radial Kα-intensity distributions and the somewhat non-intuitive dependence of the absolute Kα emission on the target thickness can be reproduced by taking into account the fast-electron refluxing with an electric field configuration based on the target normal sheath acceleration model. We infer the presence of a sheath electric field on the order of TV/m, extending to about 100 μ m in the radial direction. In addition, we obtain a temporal profile of the Kα radiation. [ABSTRACT FROM AUTHOR]

Ring-beam energetic electrons can drive two different kinetic instabilities, including the bump-on-tail instability by the beam component and the electron cyclotron maser instability by the ring component. It is critical to understand how the two instabilities develop and interact, and how their competition for free energy affects the wave modes and further plasma emission process. Here, we present fully kinetic particle-in-cell simulations of the ring-beam interaction with weakly magnetized plasmas of coronal conditions. We found that both the beam-Langmuir (BL) mode and the upper-hybrid (UH) mode can be excited efficiently, along with significant fundamental and harmonic plasma emissions. In comparison with the corresponding pure-beam case, in the ring-beam case, the BL mode has a suppressed range of k ⊥ , slower damping and thus stronger intensity, and the fundamental emission is considerably stronger, while the harmonic emission has comparable intensity yet being more directional. We further investigated the effect of the pitch angle (α) of the ring-beam distribution, found that the (1) BL mode is more limited in the range of k ⊥ with larger α and (2) the UH mode is excited in localized blobs that can be clustered with a straight line whose slope is ∼ cot α. Most of these intriguing results can be interpreted with the linear kinetic theory. [ABSTRACT FROM AUTHOR]

The characteristic behaviors of the sheath in a collisional dusty plasma composed of nonthermal electrons, ions, neutral atoms, and negatively charged dusts are investigated. To suit the realistic environment, the ionization effect is considered. The result reveals that the peculiarities of the sheath relied on ionization frequency σ , non-thermality parameter b, and dust grains concentration μ. At the sheath edge, the requirement of ion-entering-sheath-velocity enhances with increased b and μ , conversely, reduces with the increment of σ. Nevertheless, when the ionization is pronounced, the impact of b on the Bohm velocity is almost negligible. In the sheath, the increased σ leads to the reduction of the sheath thickness, which results in the redistribution of particles densities. It is noted that ion accumulation is present near the sheath edge. As expected, the present results can give more insight into the interaction processes that happened in the plasma–wall transition region. [ABSTRACT FROM AUTHOR]

Resonant interactions with electromagnetic whistler-mode waves are a primary driver of energetic electron dynamics in the Earth's radiation belts. The most intense waves can resonate with electrons nonlinearly, and effects of such nonlinear resonant interactions significantly differ from the classical quasi-linear diffusion. There have been continuous efforts on the theoretical investigation and implementation of these effects into radiation belt models, but not all nonlinear effects have been revealed yet. The two most investigated effects are phase trapping and phase bunching, which are responsible for electron acceleration and precipitation into the Earth's atmosphere, respectively, i.e., for the first cyclotron resonance with waves generated at the equator and propagating to higher latitudes, phase trapping increases electrons' energy, whereas phase bunching decreases the electron pitch-angle (and magnetic moment). However, recent studies reported a new effect called positive phase bunching, which may increase the electron pitch-angle and move them away from the loss-cone. This paper aims to characterize possible contributions of this effect to long-term electron dynamics, including multiple resonant interactions. Using an iterated mapping technique, we show that although the positive phase bunching effect can modify electron trajectories, it does not change the average rate of electron mixing in phase space. Thus, this effect may be safely neglected in long-term simulations of radiation belt dynamics. We also discuss possible verification of the positive phase bunching effect using short (single resonance), bursty electron precipitation events. [ABSTRACT FROM AUTHOR]

Among various non-linear regimes of electron cyclotron instabilities that may be driven by a strong electron cyclotron resonance (ECR) heating in open magnetic traps, the continuous generation of a quasi-monochromatic wave, characterized by a smooth release of a free energy from a system, is especially important for applications, e.g., aimed at performance optimization of ECR ion sources. We study an universal model based on the quasilinear description of plasma–wave interactions that explains the condition of a continuous-wave regime of the cyclotron instability observed experimentally including, in particular, a strong dependence on a magnetic configuration. [ABSTRACT FROM AUTHOR]

During transient instabilities in a 2 eV, highly collisional MHD-driven plasma jet experiment, evidence of a 6 keV electron tail was observed via x-ray measurements. The cause for this unexpected high energy tail is explored using numerical simulations of the Rutherford scattering of a large number of electrons and ions in the presence of a uniform electric field that is abruptly turned on as in the experiment. When the only active processes are Rutherford scattering and acceleration by the electric field, contrary to the classical Fokker–Planck theory of plasma resistivity, it is found that no steady state develops, and instead, the particle kinetic energy increases continuously. However, when a power loss mechanism is introduced mimicking atomic line radiation, a near steady state can develop and, in this case, an energetic electron tail similar to that observed in the experiment can develop. The reasons underlying this behavior are analyzed, and it is shown that an important consideration is that Rutherford scattering is dominated by the cumulative effect of grazing collisions, whereas atomic line radiation requires an approximately direct rather than a grazing collision. [ABSTRACT FROM AUTHOR]