Kinetic Plasma Research Articles

Effect of electron temperature on formation of travelling waves in plasma: Kinetic and hydrodynamic models

Abstract The kinetic formulation of the model problem of plasma waves excitation by a powerful short laser pulse is numerically studied for the first time. Kinetic and simplest hydrodynamic plasma models are also compared for the problem under consideration. It is shown that the considered hydrodynamic models do not provide good approximations to the solution to the Vlasov kinetic equation, namely, one leads to discontinuous solutions and the other has a significant qualitative distinction. At a low plasma temperature, the effect of non-isothermicity of the process is small, but it can lead to significant distortions of the solution during further heating. The results obtained here imply that the first two moments of the distribution function are not enough to describe the plasma hydrodynamics; higher-order moments should be used.

Computationally efficient high-fidelity plasma simulations by coupling multi-species kinetic and multi-fluid models on decomposed domains

A numerical method is developed for coupling a multi-species kinetic plasma model with a 5N-moment multi-fluid plasma model. The simulation domain is decomposed such that the local conditions satisfy the corresponding plasma model's region of validity. The method allows for hybrid simulations by formulating each model as a set of conservation laws and using a continuum numerical method to solve each model's governing equations in the subdomains of the decomposed domain. The models are coupled through fluxes across subdomain interfaces. Two methods are explored for the formulation of the fluxes that can be self-consistently represented by both plasma models. One method allows for flux calculations consistent with the 5N-moment multi-fluid plasma model and assumes thermodynamic equilibrium within each species of the kinetic plasma model. The second method ensures conservation of the distribution function as well as mass, momentum, and energy by formulating the fluxes using a composite underlying distribution function at the subdomain interfaces. The methods are compared in 1D1V simulations of a double rarefaction wave and a plasma sheath using the WARPXM framework, which solves each model using the discontinuous Galerkin finite element method. Both methods for formulating the fluxes perform well as the subdomain interface distribution function approaches a Maxwellian, with the consistent method being more robust to larger deviations. A simulation of the magnetized Kelvin-Helmholtz instability in 2D2V is also performed using the consistent method, which demonstrates the potential of the domain-decomposed hybrid method in facilitating speedup and reduction in required computational resources for high-fidelity plasma simulations, allowing for the investigation of problems that are beyond current capabilities.

Plasma‐enhanced evaporation and its impact on plasma properties and gaseous chemistry in a pin‐to‐water pulsed discharge

AbstractThe plasma in contact with liquids has led to various novel applications such as plasma biomedicine, material synthesis, and so on. However, the phenomenon of evaporation under plasma treatment and its impact on plasma–liquid interactions has a limited understanding. In this study, the spatially and temporally resolved behavior of water vapor production and its induced influences on plasma properties and gaseous chemistry were studied in detail in an atmospheric pressure pin‐to‐water pulsed He discharge. Diagnostic methods such as laser‐induced fluorescence (LIF) and high‐resolution optical emission spectroscopy (OES) were applied to determine the water vapor and OH radical densities, as well as key plasma parameters such as the gas temperature and electron density. It shows that the physicochemical properties of plasma vary among different discharge regions due to evaporation behavior stimulated during the pulsed discharge‐on phase. In addition, using simulation based on the experimental data, the mechanisms of how water vapor affects the observed spatiotemporal behaviors of OH radicals in different discharge regions are understood. Compared to the pin‐anode and liquid‐cathode sheath regions, proper electron parameters such as density and temperature, as well as water vapor density in the plasma‐positive column, significantly enhance the production of reactive OH radical through the dominant path of electron‐stimulated H2O dissociation. However, higher levels of electron parameters in the intense discharge region near the positive‐pin boundary enhance OH dissociation and finally result in the hollow distribution of OH density. From the global kinetic plasma simulation, the production of reactive hydroxide species playing key roles in plasma medicine treatments, such as O, H, HO2, H2O2, and hydrated ions including H+(H2O)4 and H+(H2O)5, are promoted noticeably as a result of the enhanced water evaporation process.

Phase-space Energization of Ions in Oblique Shocks

Examining energization of kinetic plasmas in phase space is a growing topic of interest, owing to the wealth of data in phase space compared to traditional bulk energization diagnostics. Via the field-particle correlation (FPC) technique and using multiple means of numerically integrating the plasma kinetic equation, we have studied the energization of ions in phase space within oblique collisionless shocks. The perspective afforded to us with this analysis in phase space allows us to characterize distinct populations of energized ions. In particular, we focus on ions that reflect multiple times off the shock front through shock-drift acceleration, and how to distinguish these different reflected populations in phase space using the FPC technique. We further extend our analysis to simulations of three-dimensional shocks undergoing more complicated dynamics, such as shock ripple, to demonstrate the ability to recover the phase-space signatures of this energization process in a more general system. This work thus extends previous applications of the FPC technique to more realistic collisionless shock environments, providing stronger evidence of the technique’s utility for simulation, laboratory, and spacecraft analysis.

Ponderomotive recoil for electromagnetic waves

When waves damp or amplify on resonant particles in a plasma, nonresonant particles experience a recoil force that conserves the total momentum between particles and electromagnetic fields. This force is important to understand, as it can completely negate current drive and rotation drive mechanisms that are predicted on the basis of only resonant particles. Here, the existing electrostatic theory of this recoil force is extended to electromagnetic waves. While the result bears close similarity to historical fluid theories of laser–plasma interactions, it now incorporates both resonant and nonresonant particles, allowing momentum conservation to be self-consistently proven. Furthermore, the result is shown to be generally valid for kinetic plasmas, which is verified through single-particle hot-plasma simulations. The new form of the force provides physical insight into the nature of the generalized Minkowski (plasmon) momentum of geometrical optics, which is shown to correspond to the momentum gained by the field and nonresonant particles as the wave is self-consistently ramped up from vanishing amplitude.

Quantifying the Agyrotropy of Proton and Electron Heating in Turbulent Plasmas

An important aspect of energy dissipation in weakly collisional plasmas is that of energy partitioning between different species (e.g., protons and electrons) and between different energy channels. Here we analyse pressure–strain interaction to quantify the fractions of isotropic compressive, gyrotropic, and nongyrotropic heating for each species. An analysis of kinetic turbulence simulations is compared and contrasted with corresponding observational results from Magnetospheric Multiscale Mission data in the magnetosheath. In assessing how protons and electrons respond to different ingredients of the pressure–strain interaction, we find that compressive heating is stronger than incompressive heating in the magnetosheath for both electrons and protons, while incompressive heating is stronger in kinetic plasma turbulence simulations. Concerning incompressive heating, the gyrotropic contribution for electrons is dominant over the nongyrotropic contribution, while for protons nongyrotropic heating is enhanced in both simulations and observations. Variations with plasma β are also discussed, and protons tend to gain more heating with increasing β.

Computer-aided identification, synthesis, and biological evaluation of DNA polymerase η inhibitors for the treatment of cancer

In cancer cells, Pol η allows DNA replication and cell proliferation even in the presence of chemotherapeutic drug-induced damages, like in the case of platinum-containing drugs. Inhibition of Pol η thus represents a promising strategy to overcome drug resistance and preserve the effectiveness of chemotherapeutic drugs. Here, we report the discovery of a novel class of Pol ƞ inhibitors, with 35 active close analogs. Compound 21 (ARN24964) stands out as the best inhibitor, with an IC50 value of 14.7 μM against Pol η and a good antiproliferative activity when used in combination with cisplatin – with a synergistic effect in three different cancer cell lines (A375, A549, OVCAR3). Moreover, it is characterized by a favorable drug-like profile in terms of its aqueous kinetic solubility, plasma and metabolic stability. Thus, ARN24964 is a promising compound for further structure-based drug design efforts toward developing drugs to solve or limit the issue of drug resistance to platinum-containing drugs in cancer patients.

Optimal lameness induction model development using amphotericin B in meat goats.

Lameness continues to be a critical health and welfare concern associated with goat production. Amphotericin B (amp B) is an antimicrobial successful in inducing transient lameness for research purposes previously in livestock animals. The objectives of this study were to (1) identify which of three varying doses of amp B would be most effective in inducing lameness in meat type goats and (2) develop a facial grimace scale for goats. Lameness was produced by an intra-articular injection of amphotericin B into the left hind lateral claw distal interphalangeal joint with either a 5 mg/0.25 mL (high-low, 5 mg of amphotericin B in a volume of 0.25 mL), 5 mg/0.5 mL (high-high, 5 mg of amphotericin B in a volume of 0.5 mL), or a 2.5 mg/0.25 mL (low-low, 2.5 mg of amphotericin B in a volume of 0.25 mL). A saline treatment of 0.5 mL was used as control (0.9% sterile saline solution). Lameness response was analyzed by infrared thermography (IRT) at the induced joint, mechanical-nociception threshold (MNT), visual lameness scoring (VLS), a visual analogue scale (VAS), kinetic gait analysis (KGA), plasma cortisol (CORT), substance P (Sub P), and behavior scoring. The IRT and MNT values differed by timepoint (P ≤ 0.0001). Results from VLS showed the HL treatment was the most effective at inducing lameness (6/6 goats became lame compared to HH 4/6 and LL 2/6). At 24, 48, and 72 h, VAS scores were significantly higher when comparing HL to all other treatment groups (P = 0.0003). Both behavior observers (1 and 2) reported a significant time effect (P = 0.05), with goats exhibiting more facial grimacing at 24 h post-lameness induction. From these data, an optimal dose for a repeatable lameness induction model in goats was aquired. An effective Goat Grimace Scale (GGS) was also developed to evaluate pain responses in goats.

Vlasov equation, waves and dispersion relations in fractal dimensions: Landau damping and the toroidal ion temperature gradient instability problem

In this study, we have used the concept of ‘product-like fractal measure’ introduced by Li and Ostoja-Starzewski in their formulation of anisotropic and elastic media to introduce Vlasov equation in fractal dimension. We have been motivated by this new concept since it is safely applied to several systems with length scales bounded by lower and upper cutoffs via the dimensional regularization technique. We have considered through this study collisionless plasma in fractal dimension and we have analyzed the impacts of the Vlasov equation in some points in plasma physics starting from the Hamiltonian dynamics. Several results have been obtained and discussed accordingly. The main outcomes of the present study concern first the facts that waves in cold plasma are associated with low fractal dimensions and hence by low complexity and that the effective damping is affected by the sign of the fractal dimension and differs from the conventional Landau damping in magnitude and in sign for negative dimension. Besides, for a very low fractal dimension, the Landau damping is negligible or removed. This result was observed in kinetic turbulent plasmas and may solve the paradox problem emerging in the application of critical balance to a kinetic turbulence cascade.

Multispecies structure-preserving particle discretization of the Landau collision operator

This paper proposes a novel numerical integrator for modeling multispecies Coulomb collisions in kinetic plasmas. The proposed scheme provides an energy-, momentum-, and positivity-preserving particle discretization of the nonlinear Landau collision operator, extending the works of Carrillo et al. [J. Comput. Phys. 7, 100066 (2020)] and Hirvijoki [Plasma Phys. Controlled Fusion 63, 044003 (2021)]. The discrete-time conservation properties are analyzed both algebraically and numerically, and an efficient, graphics processing unit-parallelized implementation is validated against inhomogeneous temperature relaxation, isotropization, and thermalization examples. The results agree with analytical estimates, confirming the method capable of reproducing physics.

Proton and Helium Heating by Cascading Turbulence in a Low-beta Plasma

How ions are energized and heated is a fundamental problem in the study of energy dissipation in magnetized plasmas. In particular, the heating of heavy ions (including 4He2+, 3He2+, and others) has been a constant concern for understanding the microphysics of impulsive solar flares. In this article, via two-dimensional hybrid-kinetic particle-in-cell simulations, we study the heating of helium ions (4He2+) by turbulence driven by cascading waves launched at large scales from the left-handed polarized helium ion cyclotron wave branch of a multi-ion plasma composed of electrons, protons, and helium ions. We find significant parallel (to the background magnetic field) heating for both helium ions and protons due to the formation of beams and plateaus in their velocity distribution functions along the background magnetic field. The heating of helium ions in the direction perpendicular to the magnetic field starts with a lower rate than that in the parallel direction, but overtakes the parallel heating after a few hundreds of the proton gyro-periods due to cyclotron resonances with mainly obliquely propagating waves induced by the cascade of injected helium ion cyclotron waves at large scales. There is, however, little evidence for proton heating in the perpendicular direction due to the absence of left-handed polarized cyclotron waves near the proton cyclotron frequency. Our results are useful for understanding the preferential heating of 3He and other heavy ions in the 3He-rich solar energetic particle events, in which helium ions play a crucial role as a species of background ions regulating the kinetic plasma behavior.

Editorial: Numerical simulations of plasma thrusters and/or related technologies

In the last few decades, electric propulsion has become in high demand on modern spacecrafts, from 8 very small nano-satellites to large communications orbital platforms. The difficulty of experimentally 9 characterizing the plasma discharge channels and the operative in-space conditions of plasma thrusters make 10 simulations a real game changer in this advanced technology/research field. In this respect, the rapid growth 11 of high-performance computing and esa-scale flops of modern architecture supercomputers has enabled 12 tackling more and more complex and computationally costly scenarios. With the use of different numerical 13 approaches, typically based on kinetic, hybrid or fluid plasma models, a large number of applications 14 can be covered: thruster physics and performance estimation, plasma thruster plume interaction with the 15 spacecraft, synthetic simulations of plasma diagnostic tools, etc.This research topic aimed to investigate, from a numerical point of view, plasma propulsion applications 17 or physical phenomena. When possible, in order to validate and/or tune the numerical tools, simulation

Alternative formulation of weak magnetohydrodynamic turbulence theory

In a recent paper [P. H. Yoon and G. Choe, Phys. Plasmas 28, 082306 (2021)], the weak turbulence theory for incompressible magnetohydrodynamics is formulated by employing the method customarily applied in the context of kinetic weak plasma turbulence theory. Such an approach simplified certain mathematical procedures including achieving the closure relationship. The formulation in the above-cited paper starts from the equations of incompressible magnetohydrodynamic (MHD) theory expressed via Elsasser variables. The derivation of nonlinear wave kinetic equation therein is obtained via a truncated solution at the second-order of iteration following the standard practice. In the present paper, the weak MHD turbulence theory is alternatively formulated by employing the pristine form of incompressible MHD equation rather than that expressed in terms of Elsasser fields. The perturbative expansion of the nonlinear momentum equation is carried out up to the third-order iteration rather than imposing the truncation at the second order. It is found that while the resulting wave kinetic equation is identical to that obtained in the previous paper cited above, the third-order nonlinear correction plays an essential role for properly calculating derived quantities such as the total and residual energies.

Kinetic plasma dynamics in a radial model of a Hall thruster with a curved magnetic field

A 1D particle-in-cell model of a Hall thruster discharge is used to analyze the effect of a curved magnetic topology in the radial plasma response and the plasma fluxes to dielectric walls. The kinetic solution shows a significant replenishment of the velocity distribution function tail and temperature isotropization for both negative (i.e. anode pointing) and positive curvatures. The new radial magnetic force is electron confining or expanding for, respectively, negative and positive curvatures, and this modifies significantly the electric and pressure radial forces. As a consequence, the plasma density near the wall and the degree of radial ion defocusing are affected: they are highly reduced for negative curvatures, the case of higher interest. For positive curvatures, the kinetic solution shows that the radial ion flow becomes supersonic within the plasma bulk, away from the Debye sheaths. An ancillary quasineutral fluid model is presented to explain this feature and other aspects of the kinetic solution. Some kinetic studies on additional phenomena complete the work.

Axisymmetric Pulsar Magnetosphere Revisited

We present a global kinetic plasma simulation of an axisymmetric pulsar magnetosphere with self-consistent e ± pair production. We use the particle-in-cell method and log-spherical coordinates with a grid size 4096 × 4096. This allows us to achieve a high voltage induced by the pulsar rotation and investigate pair creation in a young pulsar far from the death line. We find the following: (1) The energy release and e ± creation are strongly concentrated in the thin, Y-shaped current sheet, with a peak localized in a small volume at the Y-point. (2) The Y-point is shifted inward from the light cylinder by ∼15% and “breathes” with a small amplitude. (3) The dense e ± cloud at the Y-point is in ultrarelativistic rotation, which we call superrotation, because it exceeds corotation with the star. The cloud receives angular momentum flowing from the star along the poloidal magnetic field lines. (4) Gamma-ray emission peaks at the Y-point and is collimated in the azimuthal direction, tangent to the Y-point circle. (5) The separatrix current sheet between the closed magnetosphere and the open magnetic field lines is sustained by the electron backflow from the Y-point cloud. Its thickness is self-regulated to marginal charge starvation. (6) Only a small fraction of dissipation occurs in the separatrix inward of the Y-point. A much higher power is released in the equatorial plane, including the Y-point where the created dense e ± plasma is spun up and intermittently ejected through the nozzle between the two open magnetic fluxes.

Quadrature-based moment methods for kinetic plasma simulations

Quadrature-based moment methods (QBMM) are applied to rarefied gas flow and plasma physics problems represented by the Boltzmann equation or the Vlasov-Poisson system of equations. QBMM are a computationally advantageous alternative to both direct Eulerian and Lagrangian particle solvers, able to provide a noise-free solution to the Boltzmann equation and capture non-equilibrium velocity distribution functions (VDF) without excessive computational cost. To provide a closure of the moment equations, the VDF is assumed to be represented by the sum of weighted kernel functions that are placed at given velocity abscissas. Three QBMM approaches, QMOM (Dirac delta kernels), HyQMOM (hyperbolic QMOM), and EQMOM (Gaussian kernels), are applied to canonical one-dimensional rarefied gas flow and plasma physics problems. A new algorithm for EQMOM that is able to treat pathological VDFs and provide control over large abscissas for more than two nodes is developed. QMOM and EQMOM are suitable for problems where the flow does not have shock-like features (for QMOM) or a discontinuity in the VDF (for EQMOM) and where the VDF is far from non-equilibrium and features strong phase-mixing. HyQMOM is generally applicable to all test problems and yields good representations of both the moments and the VDF. HyQMOM also provides significant wall-time improvement (factor of 10–100) compared to other direct Eulerian solvers.

Ion and Electron Acceleration in Fully Kinetic Plasma Turbulence

Turbulence is often invoked to explain the origin of nonthermal particles in space and astrophysical plasmas. By means of 3D fully kinetic particle-in-cell simulations, we demonstrate that turbulence in low-β plasmas (β is the ratio of plasma pressure to magnetic pressure) accelerates ions and electrons into a nonthermal energy distribution with a power-law energy range. The ion spectrum is harder than the electron one, and both distributions get steeper for higher β. We show that the energization of electrons is accompanied by a significant energy-dependent pitch-angle anisotropy, with most electrons moving parallel to the local magnetic field, while ions stay roughly isotropic. We demonstrate that particle injection from the thermal pool occurs in regions of high current density. Parallel electric fields associated with magnetic reconnection are responsible for the initial energy gain of electrons, whereas perpendicular electric fields control the overall energization of ions. Our findings have important implications for the origin of nonthermal particles in space and astrophysical plasmas.

Quasilinear theory of general electromagnetic fluctuations including discrete particle effects for magnetized plasmas: General analysis

The general quasilinear Fokker–Planck kinetic equation for the gyrophase-averaged plasma particle distribution functions in magnetized plasmas is derived, making no restrictions on the energy of the particles and on the frequency of the electromagnetic fluctuations and avoiding the often made Coulomb approximation of the electromagnetic interactions. The inclusion of discrete particle effects breaks the dichotomy of nonlinear kinetic plasma theory divided into the test particle and the test fluctuation approximation because it provides expression of both the non-collective and collective electromagnetic fluctuation spectra in terms of the plasma particle distribution functions. Within the validity of the quasilinear approach, the resulting full quasilinear transport equation can be regarded as a determining nonlinear equation for the time evolution of the plasma particle distribution functions.

A parallel low-rank solver for the six-dimensional Vlasov–Maxwell equations

Continuum Vlasov simulations can be utilized for highly accurate modelling of fully kinetic plasmas. Great progress has been made recently regarding the applicability of the method in realistic plasma configurations. However, a reduction of the high computational cost that is inherent to fully kinetic simulations would be desirable, especially at high velocity space resolutions. For this purpose, low-rank approximations can be employed. The so far available low-rank solvers are restricted to either electrostatic systems or low dimensionality and can therefore not be applied to most space, astrophysical and fusion plasmas. In this paper we present a new parallel low-rank solver for the full six-dimensional electromagnetic Vlasov–Maxwell equations that can utilize distributed memory architectures. Special care is taken to ensure the conservation of mass and a good representation of Gauss's law. The low-rank Vlasov solver is applied to standard benchmark problems of plasma turbulence and magnetic reconnection and compared to the full grid method. It yields accurate results at significantly reduced computational cost.

Low-frequency Alfvén waves and parametric instabilities in fluid and kinetic plasmas

Low-frequency Alfvén waves and parametric instabilities in fluid and kinetic plasmas