Search

Force-free current sheets, characterized by field-aligned electric currents and approximately uniform plasma pressures, have been widely observed in the planetary magnetosphere and throughout the heliosphere. Recent observations of force-free current sheets have clearly shown the presence of anisotropic electron distributions with different temperatures perpendicular and parallel to the local magnetic field. In most of the kinetic models for one-dimensional, force-free current sheets, however, the electron distributions are nearly isotropic, which necessitates the construction of new models accounting for the electron temperature anisotropy. In this paper, we develop a model for anisotropic force-free current sheets, by incorporating the magnetic moment as an additional invariant of motion into the nearly isotropic electron distribution function of a previous model. Despite the different electron distributions, the electromagnetic profiles of the new model are often close to those in the nearly isotropic model. The applicability of our model is then validated via a comparison to a typical force-free current sheet in the Jovian magnetodisk, which shows good agreement between the model and the observations. [ABSTRACT FROM AUTHOR]

In the negative hydrogen ion sources, the positive bias of the plasma grid with respect to the discharge chamber is often used to suppress co-extracted electrons. Because most of the co-extracted electrons are magnetically deflected onto the extraction grid, the difference between the extraction and acceleration current (i.e., Iext − Iacc) is commonly used to represent the number of co-extracted electrons. In the single-driver RF negative ion source at ASIPP, the variation of Iext − Iacc with the plasma grid bias voltage was accidentally found to have different trends with and without the Cs seeding. When the current density of Iacc (regarded as JH-) was larger than 120 A/m2 during Cs conditioning, the Iext − Iacc fell first and then rose with the increasing bias voltage. The previous research results from NIFS that the beam optics were changed by the bias voltage can be used to explain the phenomenon. In this paper, the electron temperature (Te) in front of the plasma grid was found to increase from the Langmuir probe measurement, when the bias voltage was higher than a threshold of 4.5 V to the float potential. Interestingly, the bias voltage corresponding to the trough of Iext − Iacc was higher than the float potential—a similar value of 4.1 V to the threshold of Te increasing. Since the negative hydrogen ions are easily destroyed by the high-energy electrons, the increasing of Te may also cause the change of beam optics and the increase in co-extracted electrons. The electron energy distribution function was used to explain the change of plasma parameters in front of the plasma grid. [ABSTRACT FROM AUTHOR]

When the hypersonic vehicle is flying, the plasma in the area near the stagnation point of the front end of the vehicle can be approximately seen as the fully ionized dusty plasma. Due to the existence of dust particles, dusty plasma affects the communication quality of the hypersonic vehicle. In this paper, the general Boltzmann equation applicable to dusty plasmas containing electrons and the Fokker–Planck–Landau collision model are combined to derive a general formula for the electron distribution function of fully ionized dusty plasmas. Considering the contribution of the collision effect and charging effect to the dispersion relationship of fully ionized dusty plasma, the dielectric constant of fully ionized dusty plasma under an external magnetic field is solved. The Wentzel–Kramers–Brillouin method is used to calculate the attenuation coefficient (α) of the THz wave in fully ionized dusty plasma, and the influence of the external magnetic field strength and other dusty plasma parameters on the attenuation characteristics of the THz circularly polarized wave is analyzed. The research results show that the α of the THz left-hand circularly polarized wave decreases with the increase in the external magnetic field strength, while the α of the THz right-hand circularly polarized wave increases. In addition, increasing the dust particle radius, dust particle density, and electron density in a certain frequency range can increase the α of the THz circularly polarized waves. These research results provide theoretical guidance for the exploration of the interaction mechanism between the THz waves and fully ionized dusty plasma. [ABSTRACT FROM AUTHOR]

When simulating a charged particle trajectory in plasmas with Monte Carlo methods, the traditional single-scattering model is accurate but computationally expensive. The traditional multiple-scattering model reduces the computation cost but assumes an exponential distribution of the scattering angle, which overestimates the cumulative scattering angle. In this paper, we develop a new Monte Carlo method for the charged particle transport simulations. The scattering angle is naturally sampled from a modified screened Rutherford scattering cross section, in which a greatly amplified minimum deflection angle is used. With this model, the simulations give the same scattering angle distribution with the single-scattering model but only take a comparable computation cost with the multiple-scattering model when the small scattering angle and Markovian approximations can be fulfilled. This model would be especially useful for the study of inertial-confinement-fusion, which is sensitive to the alpha-particle energy deposition. [ABSTRACT FROM AUTHOR]

This paper introduces a new model to find the velocity-space location of energetic ions generating ion cyclotron emission (ICE) in plasmas. ICE is thought to be generated due to inverted gradients in the v ⊥ direction of the velocity distribution function or due to anisotropies, i.e., strong gradients in the pitch direction. Here, we invert synthetic ICE spectra generated from first principles PIC-hybrid computations to find the locations of these ICE-generating ions in velocity space in terms of a probability distribution function. To this end, we compute 2D ICE weight functions based on the magnetoacoustic cyclotron instability, which reveals the velocity-space sensitivity of ICE measurements. As an example, we analyze the velocity-space sensitivity of synthetic ICE measurements near the first 15 harmonics for plasma parameters typical for the Large Helical Device. Furthermore, we investigate the applicability of a least-square subset search, Tikhonov regularization, and Lasso regularization to obtain the locations in velocity space of the ions generating the ICE. [ABSTRACT FROM AUTHOR]

In a recent paper [K. Min, Phys. Plasmas 30, 012904 (2023)], the formation of the banded chorus with a gap in intensity at half the electron cyclotron frequency ( Ω e / 2) is demonstrated by particle-in-cell simulations including an isotropic shell distribution at an intermediate energy. This follow-up study focuses on the phase space density (PSD) hill formation process and its role in the chorus wave damping at the gap. We first show that phase-trapped particles closely follow single wave characteristics in momentum space. This means that the formation of either PSD hole or hill is primarily determined by the temperature anisotropy, T ⊥ / T ‖ , of an initial distribution function. The critical value of T ⊥ / T ‖ increases (decreases) for a higher (lower) resonant frequency. We then revisit the recent banded chorus simulations to investigate how the presence of an isotropic shell distribution self-consistently affects chorus wave evolution at the gap. Initially, with an increasing wave frequency, more and more shell electrons get trapped and a PSD hill is formed. The enhanced PSD hill counteracts wave growth driven by phase-trapped anisotropic electrons and subsequently reduces wave amplitude. The weakened wave self-consistently feeds back to the particle trapping, ultimately suppressing both the PSD hole and hill. By the time the wave frequency reaches about 0.45 Ω e , the gyro-phase structure of the electron distribution becomes much less organized. In some cases, however, the wave growth at the upstream source region can be strong enough that waves still manage to go through the gap frequency, suggesting that additional process(es) should likely be accompanied. [ABSTRACT FROM AUTHOR]

Runaway electrons exhibit kinetic instabilities with potentially beneficial consequences. The anomalous Doppler resonance between the electrons and whistler modes is a primary underlying mechanism. These instabilities require a first-principle nonlinear theoretical analysis, which would ultimately enable assessment of their impact on disruption mitigation and diagnostics, especially in ITER-relevant conditions. This paper presents recent progress in developing the required theoretical framework for isolated nonlinear resonances. By employing the generic bump-on-tail model, we predict the wave saturation levels in both near-threshold and strongly driven regimes. Remarkably, for the waves of interest, the parallel component of the wave vector dictates the parallel momentum of the resonant particles. This feature, together with the expected wave saturation level, provides a complete description of the resonance impact on the runaway electron distribution function. We show that the parametric dependence of the nonlinearly saturated mode provides a way to recover certain features of the runaway momentum distribution function. [ABSTRACT FROM AUTHOR]

It was recently shown that the use of Coulomb logarithms appropriate for bremsstrahlung radiation (rather than transport processes) along with corrections for the Langdon effect and ion screening reproduced measurements of collisional absorption in well-characterized underdense plasmas [D. Turnbull et al., Phys. Rev. Lett. 130, 145103 (2023)]. However, it was recognized at the time that the use of the standard absorption-reduction factor from Langdon's seminal paper was inconsistent with the use of Coulomb logarithms that are thermally averaged over a Maxwellian electron-velocity distribution function. A more accurate approach would be to average over the expected super-Gaussian distribution function while accounting for the Gaunt factor's velocity dependence, which somewhat mitigates the Langdon effect; however, at that time, this theory matched the data less well. This conflict is now eliminated with the additional insight that the ionization state of our mid-Z ion species (when present) was lower than had been assumed, as evidenced by the Thomson-scattering data and time-dependent Cretin simulations. We are now able to show that an improved treatment of the Langdon effect provides the best match to data. Otherwise, the prior conclusions remain unchanged. We also show an example of the substantial expected impact to the absorption rate in calculations of indirect-drive hohlraums. [ABSTRACT FROM AUTHOR]

Prompt redeposition is a process in which a neutral atom sputtered from a divertor target or the first wall of a tokamak returns on the surface right after ionization, before experiencing any collisions with the background plasma. This paper presents analytical solutions of kinetic equations for sputtered neutrals and resulting ions. Using obtained distribution functions, the redeposition coefficient and average impact angle are derived under assumptions of thin Debye sheath. Our expression for the prompt redeposition coefficient explicitly depends on the angular and energy distributions of sputtered atoms and, therefore, is more universal than often used Fussmann equation valid only for cosine distribution. The effect of the angular distributions of sputtered particles on the prompt redeposition efficiency and parameters of redeposited ions is analyzed. [ABSTRACT FROM AUTHOR]

Quasi electrostatic fluctuations in the upper-hybrid frequency range are commonly detected in the planetary magnetospheric environment. The origin of such phenomena may relate to the instability driven by a loss-cone feature associated with the electrons populating the dipole-like magnetic field. The present paper carries out a one-dimensional electrostatic particle-in-cell simulation accompanied by a reduced quasilinear kinetic theoretical analysis to investigate the dynamics of the upper-hybrid mode instability driven by an initial ring electron distribution function, which is a form of loss-cone distribution. A favorable comparison is found between the two approaches, which shows that the reduced quasilinear theory, which is grounded in the concept of a model of the particle distribution function that is assumed to maintain a fixed mathematical form except that the macroscopic parameters that define the distribution are allowed to evolve in time, can be an effective tool in the study of plasma instabilities, especially if it is guided by and validated against the more rigorous simulation result. [ABSTRACT FROM AUTHOR]

Electron cyclotron drift instability (ECDI) and ion–ion two stream instability (IITSI) are both kinetic instabilities that can be present in low-temperature, partially magnetized plasmas. The coupling of instabilities in a three-dimensional configuration leads to the existence of more than one unstable roots to the kinetic dispersion relation. In this paper, a generalized method has been developed for numerically evaluating solutions to the three-dimensional dispersion relation for coupled ECDI and IITSI, assuming cold singly and doubly charged ions and a Maxwellian velocity distribution function for the electrons. The present study demonstrates the coupling between ECDI and IITSI that affects the most unstable mode as a function of the wavenumbers in three dimensions and various plasma properties, including the applied electric field, magnetic field, electron temperature, ion velocities, and plasma density. One of the most notable results is that, while the most unstable mode with the largest growth rate is in the direction of the E × B drift in the two-dimensional cases, the most unstable mode for the three-dimensional configuration occurs in the oblique direction between the applied electric field and the E × B drift. This agrees with experimental observations in cross field plasma sources using coherent Thomson scattering. [ABSTRACT FROM AUTHOR]

Radio-frequency capacitively coupled plasmas (CCPs) are one of the key technologies enabling the latest etching processes in 3D NAND and FinFET manufacturing. These processes rely crucially on the precise control of the uniformity of ion/radical fluxes and ion angular and energy distribution function (IAEDF) in CCPs. The plasma behavior and scaling properties are dependent on the plasma chemistry in these processes, e.g., electro-positive Ar plasmas vs highly electro-negative O2 plasmas. With the large number of process and design parameters influencing the plasma properties, computational modeling has become an important tool in conjunction with experimental diagnostics in understanding the intricate physical mechanisms in CCPs. In this paper, a 2D particle-in-cell plasma model is used to study the kinetic behavior of low-pressure (<5 Pa) CCPs in two different representative chemistries: Ar and O2. The low-frequency RF source is at 1.356 MHz while 27.12 MHz is used for the high frequency. Simulations show a shift of the peak in the plasma density from the center of the chamber to the edge as the pressure increases from 0.3 to 2.6 Pa. The computed magnitude and spatial profile of electron density compare reasonably with experimental measurements over a range of pressure. Comparison between electro-positive and electro-negative plasmas are discussed. Modeling results for the dual frequency CCP highlight the effect of plasma uniformity on the IAEDF, especially near the outer edge of the electrodes. Collisions in the sheath increase the population of low-energy ions as the pressure is increased to 2.6 Pa. [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]

In regions where a plasma is not fully ionized, such as the edge and scrape-off layer (SOL) regions in a tokamak, the charged particles may be subject to strong sources from interactions with neutral atoms and molecules. Such sources, e.g., from electron impact ionization, can introduce kinetic effects, as the ionized particles may have flow velocity and temperature different from that of the main species. If treated in the conventional fluid picture, this kinetic effect emerges as a frictional heating term. In this paper, the physics of this term is discussed, both for un-magnetized and magnetized plasmas. The fluid source terms are mapped back to the kinetic sources to provide a consistent picture for future model comparison. In the limits of low and high ratios between the rates of thermalization and ionization, a multi-ion species drift-fluid model is applied to assess the impact of this kinetic effect on SOL drift-plane plasma transport. This is done by modeling a seeded blob where the ions follow either a single- or double-Maxwellian velocity distribution function (VDF). It is found that the robustness of the magnetized plasma VDF in the drift-plane and the limited effect on the vorticity source ensure that the impact of kinetic effects on the perpendicular blob evolution is small, even in the limit of high ionization to thermalization rate ratio, where kinetic effects to the ion VDF are significant. [ABSTRACT FROM AUTHOR]

In this paper, we propose a convenient method for calculating the electron energy distribution function (EEDF) in a gas-discharge plasma, where the presence of high-energy electrons results in the formation of a partial EEDF anisotropy. A detailed description of this method is provided. The performance of the proposed method is demonstrated to calculate the EEDF of plasma column of a high-voltage nanosecond discharge in argon with an extended hollow cathode. The applicability assessment of the proposed method calculating the EEDF in the plasma column of pulsed discharge sustained by the 50 Hz pulse and gas pressure from 5 to 40 Torr is given. It was shown that high-energy electrons were effectively formed in this type of discharge with an increase in the EEDF in the high-energy region. This technique can also be usefully employed for other types of discharges where fast electrons are formed. [ABSTRACT FROM AUTHOR]

Strongly magnetizing a plasma adds a range of waves that do not exist in unmagnetized plasmas and enlarges the laser-plasma interaction (LPI) landscape. In this paper, we use particle-in-cell simulations to investigate strongly magnetized LPI in one dimension under conditions relevant for magneto-inertial fusion experiments, focusing on a regime where the electron-cyclotron frequency is greater than the plasma frequency and the magnetic field is at an oblique angle with respect to the wave vectors. We show that when electron-cyclotron-like hybrid wave frequency is about half the laser frequency, the laser light resonantly decays to magnetized plasma waves via primary and secondary instabilities with large growth rates. These distinct magnetic-field-controlled instabilities, which we collectively call two-magnon decays, are analogous to two-plasmon decays in unmagnetized plasmas. Since additional phase mixing mechanisms are introduced by the oblique magnetic field, collisionless damping of large-amplitude magnetized waves substantially broadens the electron distribution function, especially along the direction of the magnetic field. During this process, energy is transferred efficiently from the laser to plasma waves and then to electrons, leading to a large overall absorptivity when strong resonances are present. The enhanced laser energy absorption may explain hotter-than-expected temperatures observed in magnetized laser implosion experiments and may also be exploited to develop more efficient laser-driven x-ray sources. [ABSTRACT FROM AUTHOR]

This paper develops an approximate steady state Fokker Planck theory for energetic electron transport in a spherical laser fusion target. First, we apply the theory to only a small population of electrons, which is specified at the outset. Then, one determines the nonlocal electron energy flux of these electrons at a particular position. These energetic electrons may either come from the tail of the thermal distribution function or be generated by an instability. This paper develops two approximate methods of solution, which we call the "characteristic method" and "sparse eigenfunction." The former works only in planar geometry and the latter in both planar and spherical geometry. Comparison of the two methods in planar geometry shows that even though the approximations are very different, they give about the same result, increasing their credibility. For the example we have chosen, spherical effects are not important for electrons from the tail of the distribution function but may well be for instability generated electrons, which have much higher energy. Comparing planar to spherical, one finds an additional spherical barrier protecting the fuel. It turns out that the associated fuel preheat a Fokker Planck model predicts is considerably less than that predicted by the Krook models as developed at both NRL and other places. [ABSTRACT FROM AUTHOR]

Analytical solutions of the Vlasov equation for periodically driven systems are of importance in several areas of plasma physics and dynamical systems and are usually approximated using ponderomotive theory. In this paper, we derive the plasma distribution function predicted by ponderomotive theory using Hamiltonian averaging theory and compare it with solutions obtained by the method of characteristics. Our results show that though ponderomotive theory is relatively much easier to use, its predictions are very restrictive and are likely to be very different from the actual distribution function of the system. We also analyse all possible initial conditions which lead to periodic solutions of the Vlasov equation for periodically driven systems and conjecture that the irreducible polynomial corresponding to the initial condition must only have squares of the spatial and momentum coordinate. The resulting distribution function for other initial conditions is aperiodic and can lead to complex relaxation processes within the plasma. [ABSTRACT FROM AUTHOR]

A zero-dimensional global model is developed with an aim to study the influence of air and water vapor on an electron energy distribution function, electron density, and temperature as well as on the chemical composition of atmospheric pressure helium plasmas, with special focus on the main reactive oxygen and nitrogen species important for applications of low temperature plasmas. The main channels for the production and consumption of electrons and electron energy gain and loss are examined by a global model based on a parametric study with variations of the mole fractions of air and water vapor. The calculations are done for 100, 1000, and 10 000 ppm of air in plasma, and for each of these values, the content of water vapor is taken to be 100, 1000, 2000, ..., to 10 000 ppm. The variations of the most important production and consumption processes for reactive oxygen and nitrogen species are analyzed in detail. According to the results presented in this paper, water vapor highly contributes to electronegativity of the plasma through pronounced attachment of electrons, which then leads to an increase in electron temperature. For high water vapor content, vibrational excitations of water molecules are one of the main electron energy loss processes, leading to a further effect on the plasma composition. Water decreases concentrations of atomic oxygen, nitrogen, and ozone, while there is an increase in nitrogen oxide, OH, H2O2, and HO2 concentrations. Cross-sectional data for electron scattering on ground and excited state neutrals are mainly taken from the Quantemol-DB database. [ABSTRACT FROM AUTHOR]

The present paper is dedicated to the study of the discrepancies encountered in electron temperature (Te) measurements carried out with electron cyclotron emission (ECE) and Thomson scattering (TS) diagnostics in the core of the JET tokamak. A large database of discharges has been collected, including high-performance scenarios performed with deuterium only and deuterium–tritium mixtures. Discrepancies have been found between core Te measurements taken with an X-mode ECE interferometer (TECE) and a LIDAR TS system (TLID) for T e > 5 keV. Depending on the plasma scenario, TECE has been found to be systematically higher or lower than TLID. Discrepancies have also been observed between the peaks of the ECE spectrum in the second (X2) and third (X3) harmonic domains, even in high optical thickness conditions. These discrepancies can be interpreted as evidence of the presence of non-Maxwellian features in the electron energy distribution function (EEDF). In order to investigate the relation between the shape of the EEDF and the measured discrepancies, a model for bipolar perturbations of Maxwellian EEDF has been developed. The model allows analytical calculations of ECE absorption and emission coefficients; hence, the comparison of modeled ECE spectra with experimental data. The different experimental results observed for the various JET scenarios have been found to be qualitatively reproducible by adapting the model parameters, suggesting that bipolar distortions of the bulk EEDF could play a role in giving rise to the reported discrepancies between ECE and TS measurements. [ABSTRACT FROM AUTHOR]

Kinetic plasma instabilities are important for regulating the temperature anisotropies of electrons and ions in solar wind. For the low beta regime, it is known that electromagnetic ion/electron cyclotron instabilities are important, but in the literature these unstable modes are discussed under the assumption of parallel propagation. The present paper extends the analysis to two (or with cylindrical symmetry, three) dimensions. The analysis is further extended to include quasilinear description with the assumption of the bi-Maxwellian velocity distribution function. Such an analysis lays the foundation for an eventual study in which cyclotron instabilities as well as obliquely propagating unstable modes such as the mirror instability are simultaneously taken into account. The present paper first lays down the basis for such future efforts in which the two- or three dimensional linear and quasilinear theories of cyclotron instabilities in the low beta regime are formulated. [ABSTRACT FROM AUTHOR]

The aim of the present paper is to examine the influence of assumption on the electron energy distribution function on the relation between the plasma potential and the electron temperature for both electropositive (argon) and electronegative (chlorine) plasmas. A one-dimensional fluid model is used for simplicity although similar results were obtained using a self-consistent two-dimensional fluid model coupled with the Maxwell's equations for inductively coupled plasmas. We find that for electropositive plasma only a bi-Maxwellian electron energy distribution function provides reasonable results compared to measurements in low-pressure inductively coupled plasmas, namely, the increasing plasma potential for increasing electron temperature. For electronegative plasma, the plasma potential is an increasing function of the electron temperature for all electron distributions considered in the model. However, the scaling factors do not agree with the conventional plasma theory. We explain these results by the deviation of electrons from a Boltzmann distribution, which is due to non-equilibrium and non-local nature of plasma at the low-pressure conditions. [ABSTRACT FROM AUTHOR]

Variable-width technology applied to the folded waveguide traveling wave tube (FW-TWT) is proposed in this paper to suppress the lower band-edge oscillation and to expand the operating bandwidth when the operating frequency shifts toward the lower band for obtaining larger output power. Changes in the width of the slow wave structure (SWS) determine the variation of the lower-cutoff-frequency (LCF), which makes the LCF's distribution a frequency range rather than a specific frequency-point. The lower band-edge oscillation will not be focused on a special position but distributed in a frequency-band range in which multiple oscillation modes compete with each other near the LCF. The lower band-edge oscillation can be suppressed by mode competition to some extent. In addition, the width variation of SWS also results in the phase-velocity variation, and thus, the operating bandwidth of FW-TWT can be further improved by the optimized design of variable-width-SWS parameters. Beam-wave interaction simulations predicted that power and bandwidth improvements have been achieved without the lower band-edge oscillation. Meanwhile, four solutions are given for further certifying the effectiveness of the variable-width method in lower band-edge oscillation suppression and bandwidth widening. An output power level of over 350 W in a frequency range of 85–90 GHz (5 GHz) is obtained, with a maximum power of ∼648 W and an electron efficiency of ∼12.1% at 88 GHz with a beam current of 260 mA and a beam voltage of 20.5 kV. [ABSTRACT FROM AUTHOR]

Extensive information, such as ion temperature and wave-particle interactions, can be obtained by the direct measurement of ion distribution functions using laser-induced fluorescence (LIF). This nonintrusive plasma diagnostics provides an important window into the ion motions in phase-space. Previous simulation results suggest that LIF measurements, which are performed on metastable ions produced directly from neutral gas particles and also from ions in other electronic states, place restrictions on the metastable lifetime. In the case where metastable population is produced from direct ionization of neutral atoms, the velocity distribution measured using LIF will only faithfully represent processes which act on the ion dynamics in a time shorter than the metastable lifetime. In this paper, the metastable lifetime effects are explored experimentally for the first time using wave-particle interaction in an Argon multidipole plasma. We demonstrate that this systematic effect can be corrected using the theory addressed in this paper based on the metastable lifetime and relative fraction of metastables produced from pre-existing ions. [ABSTRACT FROM AUTHOR]

In this paper, the single-momentum path integral Monte Carlo method, previously developed for simple quantum systems and hydrogen plasma, is adapted to simulations of the uniform electron gas. The developed method is based on the combination of Wigner formalism and the path integral approach and is able to calculate various thermodynamic values and distribution functions without differentiation of the partition function. Since the exchange interaction between electrons is taken into account by the Gram determinants of the exchange matrix, the fermionic sign problem is reduced significantly, and in the case of coordinate-depending variables, is completely eliminated. The method was applied to study thermodynamic properties of the uniform electron gas in warm dense matter regime. Average kinetic, potential, and exchange-correlation energy were calculated in a wide range of states. [ABSTRACT FROM AUTHOR]

Tailoring the ion energy distribution function (IEDF) is vital for advanced plasma processing applications. Capacitively coupled plasma (CCP) discharges excited using a non-sinusoidal waveform have shown its capability to control IEDF through the generation of plasma asymmetry and DC self-bias. In this paper, we performed a particle-in-cell simulation study to investigate the IEDF in a symmetric capacitive discharge excited by a saw-tooth-like current waveform at a very high frequency. At a constant driving frequency of 27.12 MHz, the simulation results predict that the ion energy asymmetry in the discharge scales with the discharge current amplitude. A transition from a single narrow ion energy peak to a bi-modal type IEDF is observed with an increase in the current density amplitude. Further studies at a constant current density and varying the fundamental excitation frequency show that the ion energy asymmetry enhances with a reduction in the driving frequency. Increase in the plasma asymmetry and significant DC self-bias at a lower driving frequency is observed to be one of the principal factors responsible for the observed asymmetry in the ion energy peaks. An investigation of DC self-bias and plasma potential confirms that the powered electrode energy peak corresponds to the DC self-bias with respect to the plasma potential, and the grounded electrode peak corresponds to the plasma potential. These results suggest that although lower driving frequency is beneficial for generating the discharge asymmetry and large DC self-bias, a narrow low energy IEDF is plausible in very high frequency driven CCP systems. [ABSTRACT FROM AUTHOR]

Non-adiabatic electron response leads to significant changes in ion temperature gradient (ITG) eigenmodes, leading, in particular, to fine-structures that are significantly extended along the magnetic field lines at corresponding mode rational surfaces (MRSs). These eigenmodes can nonlinearly interact with themselves to drive zonal flows via the so-called self-interaction mechanism. In this paper, the effect of collisions on these processes are studied. In the presence of non-adiabatic electrons, the linear growth rate of ITG eigenmodes decreases with the increasing collisionality. Detailed velocity space analysis of the distribution function shows that this results from collisions leading to a more adiabatic-like response of electrons away from MRSs. In linear simulations, collisions are furthermore found to broaden the radial width of the fine-structures, which translates to narrower tails of the eigenmode in extended ballooning space. The characteristic parallel scale length associated with these tails is shown to scale with the mean free path of electrons. In nonlinear turbulence simulations accounting for physically relevant values of collisionality, the fine-structures located at MRSs, together with the associated drive of zonal flows via self-interaction, are shown to persist and play a significant role. [ABSTRACT FROM AUTHOR]

This paper presents a study of the interaction between Alfvén modes and zonal structures, considering a realistic ASDEX Upgrade equilibrium. The results of gyrokinetic simulations with the global, electromagnetic, particle-in-cell code ORB5 are presented, where the modes are driven unstable by energetic particles with a bump-on-tail equilibrium distribution function, with radial density gradient. Two regimes have been observed. At low energetic particle concentration, the zonal structure (identified as an energetic particle-driven geodesic acoustic mode) is more unstable than the Alfvén mode. In the regime at high energetic particle concentration, the Alfvén mode is more unstable than the zonal structure. The interplay between the modes leads to a modification of their growth rates as well as to a modification of their saturation levels. The theoretical explanation of the mode interaction is given in terms of three-wave coupling of the energetic particle-driven geodesic acoustic mode and Alfvén mode, mediated by the curvature–pressure coupling term of the energetic particles. [ABSTRACT FROM AUTHOR]

Velocity-space anisotropy can significantly modify fusion reactivity. The nature and magnitude of this modification depends on the plasma temperature, as well as the details of how the anisotropy is introduced. For plasmas that are sufficiently cold compared to the peak of the fusion cross section, anisotropic distributions tend to have higher yields than isotropic distributions with the same thermal energy. At higher temperatures, it is instead isotropic distributions that have the highest yields. However, the details of this behavior depend on exactly how the distribution differs from an isotropic Maxwellian. This paper describes the effects of anisotropy on fusion yield for the class of anisotropic distribution functions with the same energy distribution as a 3D isotropic Maxwellian and compares those results with the yields from bi-Maxwellian distributions. In many cases, especially for plasmas somewhat below reactor-regime temperatures, the effects of anisotropy can be substantial. [ABSTRACT FROM AUTHOR]

This paper discusses the spatial resolution that is required to study inhomogeneous, low-temperature plasmas and is based on a review of low-temperature plasma electron kinetics and methods for probe measurements of electron energy distribution functions (EEDFs). It is stated that EEDFs can be extracted from probe measurements by applying an appropriate probe theory. The Druyvesteyn formula is most commonly used for this extraction and has been used in numerous publications, but more general theory can be used for a wider range of gas pressures. It is demonstrated that the Druyvesteyn formula can be obtained from the general theory as a limiting case. This paper justifies the application of wall probes in plasma studies of an energetic part of EEDFs. This justification is made for an idealized probe. We briefly review the methods for studying anisotropic plasmas and their usefulness in plasma research. It is demonstrated that to determine anisotropic electron energy distribution functions, a planar, one-sided probe is most convenient. [ABSTRACT FROM AUTHOR]

In collision-poor space plasmas, protons with an excess of kinetic energy or temperature in the direction perpendicular to the background magnetic field can excite the electromagnetic ion cyclotron (EMIC) instability. This instability is expected to be highly sensitive to suprathermal protons, which enhance the high-energy tails of the observed velocity distributions and are well reproduced by the (bi-)Kappa distribution functions. In this paper, we present the results of a refined quasi-linear approach, able to describe the effects of suprathermal protons on the extended temporal evolution of EMIC instability. It is, thus, shown that suprathermals have a systematic stimulating effect on the EMIC instability, enhancing not only the growth rates and the range of unstable wavenumbers but also the magnetic fluctuating energy density reached at the saturation. In effect, the relaxation of anisotropic temperature also becomes more efficient, i.e., faster in time and closer to isotropy. [ABSTRACT FROM AUTHOR]

A numerical algorithm is developed and tested that implements the kinetic treatment of electromagnetic radiation propagating through plasmas whose properties have small scale fluctuations, which was developed in a companion paper. This method incorporates the effects of refraction, damping, mode structure, and other aspects of large-scale propagation of electromagnetic waves on the distribution function of quanta in position and wave vector, with small-scale effects of nonuniformities, including scattering and mode conversion approximated as causing drift and diffusion in wave vector. Numerical solution of the kinetic equation yields the distribution function of radiation quanta in space, time, and wave vector. Simulations verify the convergence, accuracy, and speed of the methods used to treat each term in the equation. The simulations also illustrate the main physical effects and place the results in a form that can be used in future applications. [ABSTRACT FROM AUTHOR]

Synchrotron radiation has markedly different behavior in ∼ 10 keV and in ∼ 100 keV plasma. We show that high-energy electrons that occupy the tail of velocity distribution function have disproportionate impact on power loss of ∼ 100 keV plasma. If electrons with energy more than cutoff energy are redistributed while keeping the Maxwellian distribution function below cutoff energy intact, both emission and absorption of synchrotron radiation act to decrease the lost power. These novel radiation transport effects in non-equilibrium plasma suggest large utility in the deconfinement of high-energy electrons to reduce synchrotron radiation in applications where the radiation is deleterious. [ABSTRACT FROM AUTHOR]

The comment questions the formulation of the κ velocity distribution function used in our paper as compared to a slightly different form used by the authors. The difference in the distribution function necessarily leads to different number densities, thermal pressures, etc. We show that the restriction with their distribution function is that the macroscopic temperature (or average kinetic energy) is the same for all spatially uniform systems with a family of κ distributions including the Maxwellian case. The distribution function used in our paper and widely adopted in various studies of nonthermal systems, however, does not impose such a constraint; in particular, the temperature has κ dependence reflecting the kinetic nature of different statistical systems. The points made in the comment are trivial and misleading. [ABSTRACT FROM AUTHOR]

In a recent pair of papers, the present author discussed a self-consistent theory of asymptotically steady-state electron distribution function and Langmuir turbulence intensity in one [P. H. Yoon, Phys. Plasmas 18, 122303 (2011)] and three [P. H. Yoon, Phys. Plasmas 19, 012304 (2012)] dimensions. The resulting electron distribution function is a type of kappa distribution that features a non-Maxwellian energetic tail component. However, while the one-dimensional solution is rigorously correct, the three-dimensional solution, which was obtained using the cylindrical coordinate representation, contains two features that may be inconsistent for field-free plasmas. One is the assumption that the resonance condition can be approximated by ω-k·v≈ω-k∥v∥. Needless to say, this is not the most general condition. The second inconsistency is that while the electron distribution is isotropic in velocity, the Langmuir turbulence intensity depends on the wave propagation direction. While these features may not be too unrealistic in the presence of an implicit ambient magnetic field, they certainly cannot be correct if the plasma is genuinely unmagnetized. In the present paper, we rectify such shortcomings by properly reformulating the problem using a spherical coordinate system in a truly free-field plasma. [ABSTRACT FROM AUTHOR]

The nonlinear regime of electromagnetic oblique instabilities is investigated by means of a "noiseless" semi-Lagrangian Vlasov–Maxwell solver. Starting from an initial equilibrium configuration with two counterstreaming electron beams, qualitatively different nonlinear regimes are shown to exist depending on the nature of the solutions of the linear dispersion relation, whose properties have been discussed in the companion paper I [Ghizzo et al., Phys. Plasmas 27, 072103 (2020)]. This behavior is in contrast with existing theories of the oblique instability, which are based on the excitation of a single eigenmode at a time: nonlinear transitions toward regimes dominated by low-frequency modes are generally shown to be possible. The emphasis here is on gaining a better understanding of the multiplicity of electromagnetic oblique unstable modes and on modeling their back-reaction on plasma wave-particle interactions and energy conversion mechanisms. The latter are shown to depend on the saturation scenario of the different regimes of the oblique instability. A new regime is discussed, in which a stochastic heating occurs at the expenses of the magnetic energy first amplified by the oblique modes and in which a (reversible) violation of entropy conservation is made possible by large amplitude phase-space fluctuations of the distribution function. [ABSTRACT FROM AUTHOR]

The thin current sheet of the Earth's magnetotail is the key element of the entire magnetosphere dynamics that significantly depends on a current sheet configuration. Spacecraft measurements have shown that the ion pressure tensor in thin current sheets is essentially isotropic, but widely used two-dimensional current sheet models with isotropic pressure tensor underestimate the current density in the current sheet. In this paper, we develop a set of one-dimensional kinetic models, which demonstrate that although the plasma can be on average isotropic, due to complex kinetic features of the ion distribution function, the current density can be larger than that predicted by isotropic models. These models demonstrate that the current density in thin current sheets crucially depends on the kinetic features of the ion distribution function. [ABSTRACT FROM AUTHOR]

In this paper, the breaking process of the vacuum circuit breakers is briefly described, and the main factors that influence the post-arc process are analyzed. Based on the detachable vacuum interrupter, the experimental platform of AC breaking is built, and the experiments of the post-arc current measurement and cathode spot observation are completed. The voltage and current during the period of current zero are measured, and the images of cathode spots during the arcing process are obtained. The experimental results show that at the same frequency, the post-arc current Ipac and the post-arc charge Qpac increase with the increase in d i d t cz ( d i d t at current zero). Under the same conditions of d i d t cz , Ipac and Qpac decrease with the increase in frequency, but both are unchanged with frequency when the frequency is high enough. With the increase in current frequency, the distribution diameter of cathode spots decreases under the condition of the same d i d t cz value. According to the characteristics of cathode spot distribution, the metal vapor density under different conditions is calculated. Comparing the calculation results of metal vapor with the experimental results of post-arc current, the concept of metal vapor ionization is put forward. An ionization process of metal vapor after the current zero crossing occurs. The ionization of metal vapor will cause an increase in the number of electrons after the current zero crossing, which affects Ipac, Qpac, and sheath growth. [ABSTRACT FROM AUTHOR]

We analyze a new mechanism for the creation and confinement of energetic electrons in a mirror-configuration plasma. A Fermi–Ulam-type process, driven by end-localized coherent electrostatic oscillations, provides axial acceleration, while a natural non-adiabaticity of μ provides phase decorrelation and energy isotropization. This novel 2D combination causes the electron energy distribution function, calculated with a diffusive-loss model, to assume a Maxwellian shape with the μ non-adiabaticity, reducing loss-cone escape and annulling the absolute-barrier energy-limiting Chirikov criterion of lower dimensional models. The theoretical predictions are compared with data from an experiment. [ABSTRACT FROM AUTHOR]

Electron-beam-sustained plasmas are of vital importance for separately controlling ion flux and ion energy. In this paper, we use an implicit particle-in-cell Monte Carlo method to study plasma kinetics in an electron-beam-sustained plasma under operating conditions relevant to the use of such plasmas for polymer processing. The results indicate that the electron and ion densities are uniformly distributed because of the uniform ionization rate and heating rate. The electron-energy distribution function is Druyvesteyn-like with an ultrahigh concentration of low-energy electrons and a high-energy tail. Low-energy electrons are beneficial for protecting the substrate in material processing and a high-energy tail is useful for the precise control of plasma-gas chemistry. For ion-energy distribution functions at the electrode surface, the low-energy (< 5 eV) ion occupation rate increases with decreasing beam current or beam energy. The proportion of low-energy ions bombarding the electrode exceeds 99%, which indicates the superiority of electron-beam-generated plasma compared with a voltage- or current-driven discharge to obtain independent control of ion flux and ion energy. The results obtained herein are important for nondestructive etching in plasma processing because of the unique plasma characteristics provided by electron-beam injection. [ABSTRACT FROM AUTHOR]

Measuring the equations of state of a compressed magnetized plasma is important for both advancing fusion experiments and understanding natural systems such as stellar winds. In this paper, we present results from our experiments on the thermodynamics of compressed magnetized plasmas; we call these studies "magnetothermodynamics." In these experiments, we generate parcels of relaxed, magnetized plasma at one end of the linear Swarthmore Spheromak eXperimental device and observe their compression in a closed conducting boundary installed at the other end. Plasma parameters are measured during compression. Instances of ion heating during compression are identified by constructing a pressure-volume diagram using the measured density, temperature, and volume of the magnetized plasma. An axial scan of the ion temperature at upstream locations suggests that the increase in ion temperature arises due to the compression of the magnetized plasma in the conducting boundary. The theoretically predicted magnetohydrodynamic (MHD) and double adiabatic equations of state are compared with experimental measurements to estimate the adiabatic nature of the compressed plasma. The equilibrium of our magnetized plasmas is well-described by magnetohydrodynamics; however, we find that the MHD equation of state is not supported by our data. Our results are more consistent with the parallel Chew-Goldberger-Low equation of state, suggesting that there is significant anisotropy in the ion distribution function. [ABSTRACT FROM AUTHOR]

Chorus waves are widely investigated in space physics due to their significant role played in the radiation belts. Multiple structures of chorus frequency chirping elements, such as rising-tone structure, falling-tone structure, and hook-like structure, are observed in both space plasmas and the ground simulation experiments. In this paper, the numerical simulation on the generation of descending elements of chorus waves with parameters of a launching device, Dipole Research EXperiment, is conducted by a hybrid code named DAWN. The descending element of chorus waves is self-consistently generated without additional requirements except for a hot electron component as the free energy source. A case of ascending element under the magnetosphere parameters is also carried out to verify our simulations. Phase space distributions of electrons for the two chirping structures are analyzed, and the simulation results are consistent with previous theoretical predictions. It is obtained from the simulation results that the ascending element of chorus is observed together with an electron phase space hole (a depleted distribution function), while the descending elemet is observed together with an electron phase space clump. [ABSTRACT FROM AUTHOR]

In this paper, we present a study of anomalous diffusion using a Fokker-Planck description with fractional velocity derivatives. The distribution functions are found using numerical means for varying degree of fractionality of the stable Lévy distribution. The statistical properties of the distribution functions are assessed by a generalized normalized expectation measure and entropy in terms of Tsallis statistical mechanics. We find that the ratio of the generalized entropy and expectation is increasing with decreasing fractionality towards the well known so-called sub-diffusive domain, indicating a self-organising behavior. [ABSTRACT FROM AUTHOR]

Test-particle simulations of electrons interacting with finite-amplitude, obliquely propagating whistler waves are carried out in order to investigate the acceleration of relativistic electrons by these waves. According to the present findings, an efficient acceleration of relativistic electrons requires a narrow range of oblique propagation angles, close to the whistler resonance cone angle, when the wave amplitude is held constant at relatively low value. For a constant wave propagation angle, it is found that a range of oblique whistler wave amplitudes permits the acceleration of relativistic electrons to O(MeV) energies. An initial distribution of test electrons is shown to form a power-law distribution when plotted in energy space. It is also found that the acceleration is largely uniform in electron pitch-angle space. [ABSTRACT FROM AUTHOR]

The Sagdeev pseudo-potential approach has been employed extensively in theoretical studies to determine large-amplitude (fully) nonlinear solutions in a variety of multi-species plasmas. Although these solutions are repeatedly considered as solitary waves (and even solitons), their temporal stability has never been proven. In this paper, a numerical study of the Vlasov-Poisson system is made to follow their temporal evolution in the presence of numerical noise and thereby test their long-time propagation stability. Considering the ion-acoustic regime, both constituents of the plasma, i.e., electrons and ions are treated following their distribution functions in these sets of fully-kinetic simulations. The findings reveal that the stability of the Sagdeev solution depends on a combination of two parameters, i.e., velocity and trapping parameter. It is shown that there exists a critical value of trapping parameter for both fast and slow solutions which separates stable from unstable solutions. In the case of stable solutions, it is shown that these nonlinear structures can propagate for long periods, which confirms their status as solitary waves. Stable solutions are reported for both Maxwellian and Kappa distribution functions. For unstable solutions, it is demonstrated that the instability causes the Sagdeev solution to decay by emitting ion-acoustic wave-packets on its propagation trail. The instability is shown to take place in a large range of velocities and even for Sagdeev solutions with a velocity much higher than the ion-sound speed. Besides, in order to validate our simulation code, two precautionary measures are taken. First, the well-known effect of the ion dynamics on a stationary electron hole solution is presented as a benchmarking test of the approach. Second, In order to verify the numerical accuracy of the simulations, the conservation of energy and entropy is presented. [ABSTRACT FROM AUTHOR]

In this paper, we have studied linear and nonlinear propagation of electron acoustic waves (EAWs) comprising cold and hot populations in which the ions form the neutralizing background. The hot electrons have been assumed to follow the generalized ( r , q ) distribution which has the advantage that it mimics most of the distribution functions observed in space plasmas. Interestingly, it has been found that unlike Maxwellian and kappa distributions, the electron acoustic waves admit not only rarefactive structures but also allow the formation of compressive solitary structures for generalized ( r , q ) distribution. It has been found that the flatness parameter r , tail parameter q , and the nonlinear propagation velocity u affect the propagation characteristics of nonlinear EAWs. Using the plasmas parameters, typically found in Saturn's magnetosphere and the Earth's auroral region, where two populations of electrons and electron acoustic solitary waves (EASWs) have been observed, we have given an estimate of the scale lengths over which these nonlinear waves are expected to form and how the size of these structures would vary with the change in the shape of the distribution function and with the change of the plasma parameters. [ABSTRACT FROM AUTHOR]

The paper investigates theoretically the sensitivities of ignition capsule implosion performance on the hard x-ray spectral distribution of hohlraum. In the simulation, the hohlraum radiation is represented by a Planckian spectrum for the main drive plus a gaussian bump centered at energy Ec for preheating x-rays. Simulation results show that with the increasing of center energy Ec, the Atwood number at the fuel-ablator interface increases rapidly due to the preheating and expanding of the inner undoped CH layer. The growing of Atwood number indicates the hydrodynamic instability (HI) growth and mixing at this interface. On the other hand, the increasing of Ec results in a large density gradient scale length of ablation front and stabilizes the HI growth at ablation front. The changes of the hard x-ray spectrum have significant influences on other important implosion parameters including the ablator mass remaining, shock timing, implosion velocity, and yield as well. High-precision results on the hard x-ray spectral distribution of hohlraum are thus critical for optimizing the ignition capsule design to limit the HI growth. [ABSTRACT FROM AUTHOR]

In order to model the non-equilibrium plasma within the discharge region of a Hall thruster, the velocity distribution functions (VDFs) must be obtained accurately. A direct kinetic (DK) simulation method that directly solves the plasma Boltzmann equation can achieve better resolution of VDFs in comparison to particle simulations, such as the particle-in-cell (PIC) method that inherently include statistical noise. In this paper, a one-dimensional hybrid-DK simulation, which uses a DK simulation for heavy species and a fluid model for electrons, is developed and compared to a hybrid-PIC simulation. Time-averaged results obtained from the hybrid-DK simulation are in good agreement with hybrid-PIC results and experimental data. It is shown from a comparison of using a kinetic simulation and solving the continuity equation that modeling of the neutral atoms plays an important role for simulations of the Hall thruster discharge plasma. In addition, low and high frequency plasma oscillations are observed. Although the kinetic nature of electrons is not resolved due to the use of a fluid model, the hybrid-DK model provides spatially and temporally well-resolved plasma properties and an improved resolution of VDFs for heavy species with less statistical noise in comparison to the hybrid-PIC method. [ABSTRACT FROM AUTHOR]

The propagation and evolution of an ultrashort laser pulse in a large amplitude plasma wave are investigated based on the photon kinetic theory. The photon number distribution function for a laser pulse in the phase space is analytically obtained by solving the photon kinetic equation in the background plasma wave. And then, the behavior of the laser pulse can be described by combining the single photon dynamics and the photon number distribution function. The evolutions of the photon number density in the coordinate and frequency domain space are discussed, and broadening or compressing of the laser pulse is also displayed in this paper. In particular, the frequency shift of the entire laser pulse is analyzed, which reflects a way of energy transformation between the laser pulse and the plasma wave. [ABSTRACT FROM AUTHOR]

Analytical studies and numerical simulations show that the electron velocity distribution function in a Hall thruster discharge with crossed electric and magnetic fields is not Maxwellian. This is due to the fact that the mean free path between collisions is greater than both the Larmor radius and the characteristic dimensions of the discharge channel. However in numerical models of Hall thrusters, a hydrodynamic approach is often used to describe the electron dynamics, because discharge simulation in a fully kinetic approach requires large computing resources and is time consuming. A more accurate modeling of the electron flow in the hydrodynamic approximation requires taking into account the non-Maxwellian character of the distribution function and finding its moments, an approach that reflects the properties of electrons drifting in crossed electric and magnetic fields better than the commonly used Euler or Navier-Stokes approximations. In the present paper, an expression for the electron velocity distribution function in rarefied spatially homogeneous stationary plasma with crossed electric and magnetic fields and predominance of collisions with heavy particles is derived in the relaxation approximation. The main moments of the distribution function including longitudinal and transversal temperatures, the components of the viscous stress tensor, and of the heat flux vector are calculated. Distinctive features of the hydrodynamic description of electrons with a strongly non-equilibrium distribution function and the prospects for further development of the proposed approach for calculating the distribution function in spatially inhomogeneous plasma are discussed. [ABSTRACT FROM AUTHOR]