Plasma Physics (physics.plasm-ph)

The ability to accurately model and predict neutral transport in the boundary plasma is important for tokamak operation. Nevertheless, validation of neutral transport models can be challenging due to the difficulty in measuring neutral particle distributions. Taking advantage of the localized neutral gas puff associated with the Gas Puff Imaging (GPI) diagnostic on the Alcator C-Mod, a validation study of the neutral transport code DEGAS 2 has been performed for helium and deuterium neutrals. Absolutely calibrated measurements of helium and hydrogen line emission are compared with simulated emission from DEGAS 2, accounting for the measured gas flow rates and employing a realistic geometry. Good agreement in peak brightness and profile shape is found for a deuterium puff case. However, helium line emission measurements are found to be lower by a factor of three than that predicted in the steady state DEGAS 2 simulations for a wide range in plasma density. Discrepancies in the light emission profile shape are evident as well. As possible causes for this discrepancy, two effects are discussed. First is the effect of local cooling due to gas puff. Second is the role of time-dependent turbulence effects on neutral penetration.

Atmospheric pressure plasma jets (APPJs) are versatile devices with numerous applications. This work focuses on APPJs generated at the tip of long, flexible tubes using the jet transfer technique. The plasma source consists of a primary discharge and a secondary discharge forming the plasma jet. Discharge power measurements were carried out in a way that it was possible to separate the contribution of the primary discharge from the total power dissipated by the plasma source. Both power and effective current were analyzed under different operating conditions. The results show that the variation of the primary discharge power is much lower than the power dissipated by the plasma jet. Additionally, the electrical characteristics of the plasma device were analyzed. Notable differences were observed between the negative and positive phases of the discharge, with a more resistive load in the negative one, which suggests that the electrical equivalent circuit model changes according to the voltage polarity.

The study investigated the effects of atmospheric plasma on various dielectric materials (FR-4, PTFE, Al$_2$O$_3$) and electrode materials (copper, silver, gold-plated copper) used in surface micro-discharge plasma sources. XPS and laser microscopy were used to analyze changes in surface properties and chemical composition after 10 hours of plasma exposure. Al$_2$O$_3$ showed the highest resistance, with no significant changes in Al/O ratio or oxidation state. PTFE underwent oxidation, with fluorine substituted by oxygen, forming carbonyl, hydroxyl and aldehyde groups. FR-4 showed the most substantial changes, with etching of the epoxy matrix exposing glass fibers. This was in line with the XPS results, which show a higher O/C-ratio, less epoxy groups and higher nitrogen signals. Silver electrodes were most resilient, with sintered particles redistributed across the Al$_2$O$_3$ surface. Copper and gold-plated copper electrodes were more susceptible to oxidation and degradation, especially in areas exposed to plasma filaments. The results highlight the importance of carefully selecting materials based on resistance to plasma exposure and compatibility with the application.

In this paper, we present three-dimensional (3D) Particle-In-Cell (PIC) simulations to study the stability of 2D Bernstein-Greene-Kruskal (BGK) modes in a magnetized plasma with a finite background magnetic field. The simulations were performed using the Plasma Simulation Code (PSC) [Germaschewski et al., J. of Comp. Phys. 318, 305 (2016)], as in our recent study using 2D PIC simulations [McClung et al., Phys. Plasmas 31, 042302 (2024)], in order to see if and how the previous results would change with 3D effects. We found that solutions that are stable (unstable) in 2D simulations are still stable (unstable) in the new 3D simulations. However, the instability develops slower in 3D than in 2D and forms an unstable spiral wave structure that is in-phase along the axial direction. We have also simulated cases with an electron density bump (EDB) at the center, in addition to cases with an electron density hole (EDH) considered in our previous study, and found differences in the unstable spiral wave structures between the two cases. Additionally, we have generalized our simulations to have an increased electron thermal velocity, as well as using initial conditions solved from the complete Vlasov-Maxwell system of equations. We found that these generalizations did not change the overall behavior of the simulations and the instability that evolves.

We report observations of direct evidence of energetic protons being accelerated above ~400 keV within the reconnection exhaust of a heliospheric current sheet (HCS) crossing by NASA's Parker Solar Probe (PSP) at a distance of ~16.25 solar radii (Rs) from the Sun. Inside the extended exhaust, both the reconnection-generated plasma jets and the accelerated protons propagated toward the Sun, unambiguously establishing their origin from HCS reconnection sites located beyond PSP. Within the core of the exhaust, PSP detected stably trapped energetic protons up to ~400 keV, which is approximately 1000 times greater than the available magnetic energy per particle. The differential energy spectrum of the accelerated protons behaved as a pure power-law with spectral index of about -5. Supporting simulations using the kglobal model suggest that the trapping and acceleration of protons up to ~400 keV in the reconnection exhaust is likely facilitated by merging magnetic islands with a guide field between ~0.2-0.3 of the reconnecting magnetic field, consistent with the observations. These new results, enabled by PSP's proximity to the Sun, demonstrate that magnetic reconnection in the HCS is a significant new source of energetic particles in the near-Sun solar wind. The discovery of in-situ particle acceleration via magnetic reconnection at the HCS provides valuable insights into this fundamental process which frequently converts the large magnetic field energy density in the near-Sun plasma environment and may be responsible for heating the sun's atmosphere, accelerating the solar wind, and energizing charged particles to extremely high energies in solar flares.

We consider the issue of strict, fully discrete \emphlocal energy conservation for a whole class of fully implicit local-charge- and global-energy-conserving particle-in-cell (PIC) algorithms. Earlier studies demonstrated these algorithms feature strict global energy conservation. However, whether a local energy conservation theorem exists (in which the local energy update is governed by a flux balance equation at every mesh cell) for these schemes is unclear. In this study, we show that a local energy conservation theorem indeed exists. We begin our analysis with the 1D electrostatic PIC model without orbit-averaging, and then generalize our conclusions to account for orbit averaging, multiple dimensions, and electromagnetic models (Darwin). In all cases, a temporally, spatially, and particle-discrete local energy conservation theorem is shown to exist, proving that these formulations (as originally proposed in the literature), in addition to being locally charge conserving, are strictly locally energy conserving as well. In contrast to earlier proofs of local conservation in the literature \citepxiao2017local, which only considered continuum time, our result is valid for the fully implicit time-discrete version of all models, including important features such as orbit averaging. We demonstrate the local-energy-conservation property numerically with a paradigmatic numerical example.

The turbulence characteristics of the scrape-off-layer (SOL) plasma in the W7-X stellarator are investigated using a gas-puff-imaging (GPI) diagnostic, newly installed and operated during the OP 2.1 campaign. The SOL plasma on W7-X features a chain of magnetic islands intersected by discrete divertor plates at five separate toroidal locations, forming a set of island divertors for heat and particle exhaust. GPI measures 2D emission fluctuations in the magnetic island plasma. A dedicated experiment was conducted to characterize the broadband turbulence and filament dynamics in the island divertor SOL in the standard magnetic configuration. This involves shifting the poloidal location of the magnetic island with respect to the GPI field of view. Overall, initial analyses confirm minimal intermittent type fluctuations, compared to tokamaks. A relationship between the radial profile of the connection length and the sheared poloidal flow structure is discussed. Turbulence statistical analyses reveal poloidal asymmetry of the fluctuation properties and a birth region of intermittent fluctuations indicated by skewness reversal. Their birth locations appear sensitive to plasma conditions. Finally, 2D cross-correlation analysis is presented to study turbulence propagation in the presence of multiple shear layers in the poloidal flows, showing a dominant poloidal motion and evidencing a shear effect.

Collisionless tearing instability with a power-law distribution function in a relativistic pair plasma with a guide field is studied. When the current sheet is supported by plasma pressure, the tearing mode is suppressed as the particle spectrum hardens. In the force-free limit, the instability growth rate becomes independent of the particle spectrum. We apply these results to relativistic MHD turbulence, where magnetic energy greatly exceeds plasma rest energy, and derive an expression for the transverse size of turbulent eddies unstable to tearing mode. We also establish the critical plasma magnetization parameter above which charge starvation prevents the tearing instability. These results might be useful for developing more accurate models of particle acceleration in relativistic astrophysical sources.

Pair creation by a spatially homogeneous, time-dependent electric field is studied within the framework of scalar quantum electrodynamics. We employ the standard Bogoliubov transformation approach to compute the single-particle distribution function in an adiabatic basis. We analyzed the distribution function of created particles in two different adiabatic bases (related by a unitary transformation). A novel dynamical scaling is observed while analyzing the oscillatory momentum spectrum of the pairs created by the Sauter pulsed field at intermediate times, calculated using the two adiabatic bases. In these bases, the same oscillatory momentum spectra are observed but at different times. However, when we scale the time by the point marking the end of the transient stage of dynamical evolution for each case of central momentum, the respective momentum spectra overlap. Furthermore, we study the time evolution of the momentum spectrum in the multi-photon regime and find that the spectra show a multi-modal profile structure at finite times for both choices of basis.