Physics

In the white paper we describe a national program in Low Temperature Plasma (LTP). The program should take advantage of the research opportunities of 3 rapidly growing areas (nanomaterial plasma synthesis, plasma medicine, microelectronics). The main theme is to achieve a fundamental understanding of Low Temperature Plasmas as they are applied to these different applications. This understanding will allow U.S. industry to meet the challenges of international competition.

This manuscript presents verification cases that are developed to study the electrothermal instability (ETI). Specific verification cases are included to ensure that the unit physics components necessary to model the ETI are accurate, providing a path for fluid-based codes to effectively simulate ETI in the linear and nonlinear growth regimes. Two software frameworks with different algorithmic approaches are compared for accuracy in their ability to simulate diffusion of a magnetic field, linear growth of the ETI, and a fully nonlinear ETI evolution. The nonlinear ETI simulations show early time agreement, with some differences emerging, as noted in the wavenumber spectrum, late into the nonlinear development of ETI. A sensitivity study explores the role of equation-of-state (EOS), vacuum density, and vacuum resistivity. EOS and vacuum resistivity are found to be the most critical factors in the modeling of nonlinear ETI development.

A multi-species magnetized collisional nonthermal plasma system containing inertial ion species, non-inertial electron species following nonthermal kappa-distribution, and immobile dust particles are considered to model the dissipative dust-ion-acoustic (DIA) soliton modes, both theoretically and numerically. The electrostatic solitary modes are found to be associated with the low frequency dissipative dust-ion-acoustic solitary waves (DIASWs). The ion-neutral collision is taken into account, and the influence of ion-neutral collisional effects on the dynamics of dissipative DIASWs is investigated. It is reported that most of the plasma medium in space and laboratory are far from thermal equilibrium, and the particles in such plasma system are well fitted via the kappa-nonthermal distribution than via the thermal Maxwellian distribution. The reductive perturbation approach is adopted to derive the damped KdV equation, and the solitary wave solution of the damped KdV equation is derived via the tangent hyperbolic method to analyze the basic features (amplitude, width, speed, time evolution, etc.) of dissipative DIASWs. The propagation nature and also the basic features of dissipative DIASWs are seen to influence significantly due to the variation of the plasma configuration parameters and also due to the variation of the supethermality index kappa in the considered plasma system. The implication of the results of this study could be useful for better understanding the electrostatic localized disturbances, in the ion length and time scale, in space and experimental dusty plasmas, where the presence of excess energetic electrons and ion-neutral collisional damping are accountable.

The JET baseline scenario is being developed to achieve high fusion performance and sustained fusion power. However, with higher plasma current and higher input power, an increase in pulse disruptivity is being observed. Although there is a wide range of possible disruption causes, the present disruptions seem to be closely related to radiative phenomena such as impurity accumulation, core radiation, and radiative collapse. In this work, we focus on bolometer tomography to reconstruct the plasma radiation profile and, on top of it, we apply anomaly detection to identify the radiation patterns that precede major disruptions. The approach makes extensive use of machine learning. First, we train a surrogate model for plasma tomography based on matrix multiplication, which provides a fast method to compute the plasma radiation profiles across the full extent of any given pulse. Then, we train a variational autoencoder to reproduce the radiation profiles by encoding them into a latent distribution and subsequently decoding them. As an anomaly detector, the variational autoencoder struggles to reproduce unusual behaviors, which includes not only the actual disruptions but their precursors as well. These precursors are identified based on an analysis of the anomaly score across all baseline pulses in two recent campaigns at JET.

A unified theoretical description is developed for the formation of an ionized filament channel in a dense-gas medium and the evolution of electronic degrees of freedom in this channel in the laser pulse wake, as illustrated on an example of high-pressure argon. During the laser pulse, the emerging free electrons gain energy via inverse Bremsstrahlung on neutral atoms, enabling impact ionization and extensive collisional excitation of the atoms. A kinetic model of these processes produces the radial density distributions in the immediate wake of the laser pulse. After the pulse, the thermalized electron gas drives the system evolution via impact ionization (from the ground and excited states) and collisional excitation of the residual neutral atoms, while the excited atoms are engaged in Penning ionization. The interplay of these three processes determines the electron gas cooling dynamics. The local imbalance of the free-electron and ion densities induces a transient radial electric field, which depends critically on the electron temperature. The evolving radial profiles of the electron, ion, and excited-atom densities, as well as the profiles of electron temperature and induced electric field are obtained by solving the system of diffusion-reaction equations numerically. All these characteristics evolve with two characteristic timescales, and allow for measuring the electronic stage of the wake channel evolution via linear and nonlinear light-scattering experiments.

We develop a universal framework which allows quickly solve a wide class of problems for longitudinal space charge effects in beams and plasmas in cylindrical geometry. We introduce the longitudinal dielectric permittivity for the beam of charged particles, which describes its collective space charge response. The analyis yields an effective plasma frequency, which depends on the transverse geometry of the system. This dielectric permittivity mirrors the dielectric permittivity of plasma and matches the one dimensional (1D) expression once the transverse size of the beam is large. Several particle species can be included as additive terms describing susceptibility of each specie. The developed approach allows to study stability criteria for collective beam-beam and beam-plasma instabilities for arbitrary transverse distributions in particle densities.

We apply field-particle correlations -- a technique that tracks the time-averaged velocity-space structure of the energy density transfer rate between electromagnetic fields and plasma particles -- to data drawn from a hybrid Vlasov-Maxwell simulation of Alfvén Ion-Cyclotron turbulence. Energy transfer in this system is expected to include both Landau and cyclotron wave-particle resonances, unlike previous systems to which the field-particle correlation technique has been applied. In this simulation, the energy transfer rate mediated by the parallel electric field E ∥ comprises approximately 60% of the total rate, with the remainder mediated by the perpendicular electric field E ⊥ . The parallel electric field resonantly couples to protons, with the canonical bipolar velocity-space signature of Landau damping identified at many points throughout the simulation. The energy transfer mediated by E ⊥ preferentially couples to particles with v tp ≲ v ⊥ ≲3 v tp in agreement with the expected formation of a cyclotron diffusion plateau. Our results demonstrate clearly that the field-particle correlation technique can distinguish distinct channels of energy transfer using single-point measurements, even at points in which multiple channels act simultaneously, and can be used to determine quantitatively the rates of particle energization in each channel.

We show that the recently introduced two-field flux-balanced Hasegawa-Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa-Wakatani models, which confirms the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multiscale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.

Tunable diode laser induced fluorescence (TDLIF) measurements are discussed and quantitatively evaluated for nitrogen admixtures in argon plasma under the influence of a strong magnetic field. TDLIF measurements were used to evaluate light-transport properties in a strongly magnetized optically thick argon/nitrogen plasma under different pressure conditions. Therefore, a coupled system of rate balance equations was constructed to describe laser pumping of individual magnetic sub-levels of 2 p 8 state through frequency-separated sub-transitions originating from 1 s 4 magnetic sub-levels. The density distribution (alignment) of 2 p 8 multiplet was described by balancing laser pumping with losses including radiative decay, transfer of excitation between the neighboring sub-levels in the 2 p 8 multiplet driven by neutral collisions (argon and nitrogen) and quenching due to electron and neutral collisions. Resulting 2 p 8 magnetic sub-level densities were then used to model polarization dependent fluorescence, considering self-absorption, which could be directly compared with polarization-resolved TDLIF measurements. This enables to estimate the disalignment rate constant for the 2 p 8 state due to collisions by molecular nitrogen. A comparison to molecular theory description is given providing satisfactory agreement. The presented measurement method and model can help to describe optical emission of argon and argon-nitrogen admixtures in magnetized conditions and provides a basis for further description of optical emission spectra in magnetized plasmas.

Gyrokinetic theory is arguably the most important tool for numerical studies of transport physics in magnetized plasmas. However, exact local energy-momentum conservation law for the electromagnetic gyrokinetic system has not been found despite continuous effort. Without such a local conservation law, energy-momentum can be instantaneously transported across spacetime, which is unphysical and casts doubt on the validity of numerical simulations based on the gyrokinetic theory. Standard Noether's procedure for deriving conservation laws from corresponding symmetries does not apply to gyrokinetic systems because the gyrocenters and electromagnetic field reside on different manifolds. To overcome this difficulty, we developed a high-order field theory on heterogeneous manifolds for classical particle-field systems and apply it to derive exact local conservation laws, in particular the energy-momentum conservation law, for the electromagnetic gyrokinetic system. A weak Euler-Lagrange equation is established to replace the standard Euler-Lagrange equation for the particles. It is discovered that an induced weak Euler-Lagrange current enters the local conservation laws. And it is the new physics captured by the high-order field theory on heterogeneous manifolds.