Plasma Dynamics Research Articles

Radio Instrument Package for Lunar Ionospheric Observation: A Concept Study

AbstractThe lunar ionosphere is a ∼100 km thick layer of electrically charged plasma surrounding the moon. Despite knowledge of its existence for decades, the structure and dynamics of the lunar plasma remain a mystery due to lack of consistent observational capacity. An enhanced observational picture of the lunar ionosphere and improved understanding of its formation/loss mechanisms is critical for understanding the lunar environment as a whole and assessing potential safety and economic hazards associated with lunar exploration and habitation. To address the high priority need for observations of the electrically charged constituents nikear the lunar surface, we introduce a concept study for the Radio Instrument Package for Lunar Ionospheric Observation (RIPLIO). RIPLIO would consist of a multi‐CubeSat constellation (at least two satellites) in lunar orbit for the purpose of conducting “crosslink” radio occultation measurements of the lunar ionosphere, with at least one satellite carrying a very high frequency (VHF) transmitter broadcasting at multiple frequencies, and at least one satellite flying a broadband receiver to monitor transmitting satellites. Radio occultations intermittently occur when satellite‐to‐satellite signals cross through the lunar ionosphere, and the resulting phase perturbations of VikHF signals may be analyzed to infer the ionosphere electron content and high‐ resolution vertical electron density profiles. As demonstrated in this study, RIPLIO would provide a novel means for lunar observation, with the potential to provide long‐term, high‐resolution observations of the lunar ionosphere with unprecedented pan‐lunar detail.

Role of electron temperature at extremely low density in negative positive ion plasma

It is pointed out that the role of electrons in the dynamics of pair ion and negative positive ion plasmas cannot be neglected even at extremely low density of electrons, i.e., ne0n+0≪mem+ (where nj0 is the background density of jth species and mj is the mass of the particles of the j-species while j = e, +, −) because electron thermal velocity is almost always larger than thermal velocities of ions, i.e., vT± ≪ vTe. An analysis of electrostatic waves in unmagnetized negative positive ion electron (NPIE) and pair ion electron (PIE) plasmas is presented for the case ωpe ≪ ωp+ (where ωpj=(4πnj0e2mj)1/2 is the plasma oscillation frequency corresponding to j-species). The electron dynamics contribute to electrostatic perturbations at ion plasma oscillation time scale at longer wavelengths for λDe2k2<1 where λDe=(Te4πne0)1/2 is the electron Debye length. On the other hand, the electron plasma wave turns into thermal wave when the conditions ωpe ≪ ωp± and 1≪λDe2k2 hold simultaneously and ion acoustic wave approaches the sum of ion plasma oscillation frequencies of positive and negative ions. The only electrostatic normal mode of such a plasma is the ion plasma wave corresponding to longer wavelengths, which also includes the contribution of electrons. The electron thermal wave is separated from plasma oscillations and electron time scale disappears with respect to electrostatic perturbations. Similar situation also occurs in plasmas having negatively charged dust particles. To elaborate these points, the analytical results are applied to the two experiments with NPIE and PIE plasmas.

Nonlinear Interaction of Low-frequency Alfvén Waves and Ions

The interaction between Alfvén waves and particles is a critical phenomenon in space and laboratory plasmas, and it has been observed that low-frequency Alfvén waves can accelerate and heat ions through subharmonic resonant interactions. In this study, we use test particle simulations to verify the nonlinear heating of parallelly propagating low-frequency Alfvén waves on ions and analyze the underlying process in terms of the Poincaré section. Our results demonstrate that low-frequency Alfvén waves can periodically pick up ions, leading to oscillations of average parallel velocity and temperature of the plasma by phase mixing, ultimately resulting in the stabilization of acceleration. Furthermore, we have developed an analysis that can estimate the time required for heating and accelerating. In the case of multiple waves, our findings indicate that the presence of more chaotic modes does not necessarily result in better wave heating. We have also discussed the effect of random phases on the heating process. Overall, this research sheds light on the crucial role played by the interaction between Alfvén waves and particles in astrophysics and provides new insights into the mechanisms underlying the heating and acceleration of ions through subharmonic resonant interactions with low-frequency Alfvén waves. These findings may have significant implications for the understanding of plasma dynamics in a range of astrophysical environments.

First measurement of the 7Li(D, n) astrophysical S-factor in laser-induced full plasma

Nuclear fusion reactions which produce energies and most of the elements in the Universe are the keys to driving the process of the Big Bang nucleosynthesis, as well as the evolution of stars. For fusions in real astrophysical conditions, it is crucial to directly measure them in the plasma environments. However, it is challenging to obtain reaction cross sections quantitatively because of the high instabilities of plasma. Here, an innovative self-calibration method is proposed to eliminate the large uncertainties caused by the plasma instabilities. And for the first time, it is applied to investigating the deuteron-lithium nuclear-fusion reaction 7Li(D,n)8Be in the plasma, which is related to the Big Bang lithium abundance puzzle. Further, the reaction's astrophysical S-factor in the Gamow window around 173 keV in the plasma is measured to be S=(24±8) MeV⋅barn, which demonstrates from an experimental point of view that the coulomb screening effect does not play an essential role in this 7Li consuming reaction related to the Big Bang lithium abundance puzzle. This self-calibration method for nuclear reactions in the plasma can also be applied to exploring other nuclear-fusion reactions, such as 7Li(p,α)4He and 6Li(p,α)3He, as well as providing a new probe for detecting plasma dynamics.

Statistical Analysis of Plasma Dynamics in Gyrokinetic Simulations of Stellarator Turbulence.

A geometrical method for assessing stochastic processes in plasma turbulence is investigated in this study. The thermodynamic length methodology allows using a Riemannian metric on the phase space; thus, distances between thermodynamic states can be computed. It constitutes a geometric methodology to understand stochastic processes involved in, e.g., order-disorder transitions, where a sudden increase in distance is expected. We consider gyrokinetic simulations of ion-temperature-gradient (ITG)-mode-driven turbulence in the core region of the stellarator W7-X with realistic quasi-isodynamic topologies. In gyrokinetic plasma turbulence simulations, avalanches, e.g., of heat and particles, are often found, and in this work, a novel method for detection is investigated. This new method combines the singular spectrum analysis algorithm with a hierarchical clustering method such that the time series is decomposed into two parts: useful physical information and noise. The informative component of the time series is used for the calculation of the Hurst exponent, the information length, and the dynamic time. Based on these measures, the physical properties of the time series are revealed.

Origin and trends in NH4+ observed in the coma of 67P

ABSTRACT The European Space Agency/Rosetta mission escorted comet 67P/Churyumov–Gerasimenko and witnessed the evolution of its coma from low activity (∼2.5–3.8 au) to rich ion-neutral chemistry (∼1.2–2.0 au). We present an analysis of the ion composition in the coma, focusing on the presence of protonated high proton affinity (HPA) species, in particular $\mathrm{{ NH}_{4}}^{+}$. This ion is produced through the protonation of NH3 and is an indicator of the level of ion-neutral chemistry in the coma. We aim to assess the importance of this process compared with other $\mathrm{{ NH}_{4}}^{+}$ sources, such as the dissociation of ammonium salts embedded in dust grains. The analysis of $\mathrm{{ NH}_{4}}^{+}$ has been possible thanks to the high mass resolution of the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis/Double Focusing Mass Spectrometer (ROSINA/DFMS). In this work, we examine the $\mathrm{{ NH}_{4}}^{+}$ data set alongside data from the Rosetta Plasma Consortium instruments, and against outputs from our in-house ionospheric model. We show that increased comet outgassing around perihelion yields more detections of $\mathrm{{ NH}_{4}}^{+}$ and other protonated HPA species, which results from more complex ion-neutral chemistry occurring in the coma. We also reveal a link between the low magnetic field strength associated with the diamagnetic cavity and higher $\mathrm{{ NH}_{4}}^{+}$ counts. This suggests that transport inside and outside the diamagnetic cavity is very different, which is consistent with 3D hybrid simulations of the coma: non-radial plasma dynamics outside the diamagnetic cavity is an important factor affecting the ion composition.

Quantum algorithm for the linear Vlasov equation with collisions

The Vlasov equation is a nonlinear partial differential equation that provides a first-principles description of the dynamics of plasmas. Its linear limit is routinely used in plasma physics to investigate plasma oscillations and stability. In this paper, we present a quantum algorithm that simulates the linearized Vlasov equation with and without collisions, in the one-dimensional electrostatic limit. Rather than solving this equation in its native spatial and velocity phase space, we adopt an efficient representation in the dual space yielded by a Fourier-Hermite expansion. For a given simulation time, the Fourier-Hermite representation is exponentially more compact, thus yielding a classical algorithm that can match the performance of a previously proposed quantum algorithm for this problem. This representation results in a system of linear ordinary differential equations (ODEs) which can be solved with well-developed quantum algorithms: a Hamiltonian simulation in the collisionless case, and quantum ODE solvers in the collisional case. In particular, we demonstrate that a quadratic speedup in system size is attainable.

Three-dimensional analysis of THz radiation mechanisms in interaction of arbitrarily polarized multi-color laser pulses and obliquely magnetized plasma

Three-dimensional THz radiation generation via wave-mixing and photo-current mechanisms of multi-color femtosecond laser fields in obliquely magnetized plasma for various laser polarization states is considered. The main motive behind present study has been to investigate the development of more efficient and versatile THz radiation generation for a range of applications through a comprehensive mechanism describing the entire THz frequency band. Two and three-color laser waves interactions with plasma electrons results in a transverse nonlinear current density which causes THz radiation generation. Variation of THz electric field under different laser polarizations is examined analytically and numerically. The impacts of the applied oblique magnetic field and the field oblique angle on the THz electric field for both wave mixing and photo-current schemes are studied and result are compared. A method to implement THz radiation for special applications was developed showing that the external magnetic field is crucial in controlling plasma dynamics. The waveforms of THz angular distribution in forward direction are examined and the impact of different harmonics, various polarization states, external magnetic field, oblique angle of the applied external magnetic field, and the length of plasma on generated waveforms is considered. Based on the results of study, laser harmonics and applied external magnetic field increase electrons current density, leading to enhancement of radiation power without change in directionality diagram. These findings are in close agreement with experimental results.

An overview of new measurements of flow, chirality and vorticity from STAR experiment

In relativistic heavy-ion collisions, the properties of quark–gluon plasma (QGP) and complex dynamics of multi-scale processes in Quantum Chromodynamics (QCD) are studied by analyzing the final state produced particles in a variety of different ways. In these proceedings, we present an overview of new detailed measurements of flow, chirality and vorticity by the STAR experiment at Relativistic Heavy Ion Collider (RHIC). Furthermore, STAR’s future opportunities for the precision measurements on small systems, fixed-target (FXT) mode, and Beam Energy Scan (BES-II) program are discussed.

Quantum three-wave instability

For the three-wave interaction, the lowest-order nonlinear interaction in plasma dynamics, we describe how the quantum system, which is time-independent, nonchaotic, finite-dimensional, and Hermitian, can give rise to a linear instability corresponding to that in the classical system. We show that the instability is realized in the quantum regime as a cascade of the wave function in the space of occupation number states, and the unstable quantum system has a richer spectrum and a much longer recurrence time than the stable quantum system. The conditions for instability of the quantum three-wave interaction are described.

Magnetic nozzle performance in a cluster of helicon plasma thrusters

A numerical study of the plasma dynamics in a Helicon Plasma Thrusters’ (HPT) cluster is presented. For the first time in the literature, the three-dimensional (3D) plasma dynamics occurring in the plume of a HPTs’ cluster is analyzed. The physical investigation relies on ProPic, a 3D particle-in-cell (PIC) code specifically designed to simulate the plasma dynamics in magnetic nozzles and in a non-axi-symmetric domain. The code has been validated against experiments reported in the literature and cross-validated with Starfish, an open-source two-dimensional PIC software. The physical investigation has revealed an interesting mutual influence between the thrusters that constitute the cluster. Three significant phenomena that affect the cluster’s performance have been identified. The first phenomenon is related to the effect that clustering has on the shape of the magnetic field lines and, in turn, on the divergence angle of the plume. The second phenomenon is related to electron currents flowing among different thrusters, which affect the potential drop across the plume. The third phenomenon is related to the effect that neighboring thrusters have on the plasma potential map and, in turn, on the expansion of the ions.

The Merging of a Coronal Dimming and the Southern Polar Coronal Hole

We report on the merging between the southern polar coronal hole and an adjacent coronal dimming induced by a coronal mass ejection on 2022 March 18, resulting in the merged region persisting for at least 72 hr. We use remote sensing data from multiple co-observing spacecraft to understand the physical processes during this merging event. The evolution of the merger is examined using Extreme-UltraViolet (EUV) images obtained from the Atmospheric Imaging Assembly on board the Solar Dynamic Observatory and Extreme Ultraviolet Imager, which is on board the Solar Orbiter spacecraft. The plasma dynamics are quantified using spectroscopic data obtained from the EUV Imaging Spectrometer on board Hinode. The photospheric magnetograms from the Helioseismic and Magnetic Imager are used to derive the magnetic field properties. To our knowledge, this work is the first spectroscopical analysis of the merging of two open-field structures. We find that the coronal hole and the coronal dimming become indistinguishable after the merging. The upflow speeds inside the coronal dimming become more similar to that of a coronal hole, with a mixture of plasma upflows and downflows observable after the merging. The brightening of the bright points and the appearance of coronal jets inside the merged region further imply ongoing reconnection processes. We propose that component reconnection between the coronal hole and coronal dimming fields plays an important role during this merging event because the footpoint switching resulting from the reconnection allows the coronal dimming to intrude onto the boundary of the southern polar coronal hole.

The force balance of electrons during kinetic anti-parallel magnetic reconnection

Fully kinetic simulations are applied to the study of 2D anti-parallel reconnection, elucidating the dynamics by which the electron fluid maintains force balance within both the ion diffusion region (IDR) and the electron diffusion region (EDR). Inside the IDR, magnetic field-aligned electron pressure anisotropy (pe∥≫pe⊥) develops upstream of the EDR. Compared to previous investigations, the use of modern computer facilities allows for simulations at the natural proton to electron mass ratio mi/me=1836. In this high-mi/me-limit, the electron dynamics change qualitatively, as the electron inflow to the EDR is enhanced and mainly driven by the anisotropic pressure. Using a coordinate system with the x-direction aligned with the reconnecting magnetic field and the y-direction aligned with the central current layer, it is well known that for the much studied 2D laminar anti-parallel and symmetric scenario the reconnection electric field at the X-line must be balanced by the ∂pexy/∂x and ∂peyz/∂z off-diagonal electron pressure stress components. We find that the electron anisotropy upstream of the EDR imposes large values of ∂pexy/∂x within the EDR, and along the direction of the reconnection X-line, this stress cancels with the stress of a previously determined theoretical form for ∂peyz/∂z. The electron frozen-in law is instead broken by pressure tensor gradients related to the direct heating of the electrons by the reconnection electric field. The reconnection rate is free to adjust to the value imposed externally by the plasma dynamics at larger scales.

Investigation of the plasma column evolution and dynamics in the PW-7 plasma accelerator

Investigation of the plasma column evolution and dynamics in the PW-7 plasma accelerator

Data-driven Radiative Magnetohydrodynamics Simulations with the MURaM Code

We present a method of conducting data-driven simulations of solar active regions and flux emergence with the MURaM radiative magnetohydrodynamics (MHD) code. The horizontal electric field that is derived from the full velocity and magnetic vectors is implemented at the photospheric (bottom) boundary to drive the induction equation. The energy equation accounts for thermal conduction along magnetic fields, optically thin radiative loss, and heating of coronal plasma by viscous and resistive dissipation, which allows for a realistic presentation of the thermodynamic properties of coronal plasma that are key to predicting the observational features of solar active regions and eruptions. To validate this method, the photospheric data from a comprehensive radiative MHD simulation of solar eruption (the ground truth) are used to drive a series of numerical experiments. The data-driven simulation reproduces the accumulation of free magnetic energy over the course of flux emergence in the ground truth with an error of 3%. The onset time is approximately 8 minutes delayed compared to the ground truth. However, a precursor-like signature can be identified at the correct onset time. The data-driven simulation captures key eruption-related emission features and plasma dynamics of the ground truth flare over a wide temperature span, from to . The evolution of the flare and coronal mass ejection as seen in synthetic extreme ultraviolet images is also reproduced with high fidelity. This method helps to understand the evolution of magnetic field in a more realistic coronal environment and to link the magnetic structures to observable diagnostics.

Extending the Spectral Difference Method with Divergence Cleaning (SDDC) to the Hall MHD Equations

The Hall Magnetohydrodynamic (MHD) equations are an extension of the standard MHD equations that include the “Hall” term from the general Ohm’s law. The Hall term decouples ion and electron motion physically on the ion inertial length scales. Implementing the Hall MHD equations in a numerical solver allows more physical simulations for plasma dynamics on length scales less than the ion inertial scale length but greater than the electron inertial length. The present effort is an important step towards producing physically correct results to important problems, such as the Geospace Environmental Modeling (GEM) Magnetic Reconnection problem. The solver that is being modified is currently capable of solving the resistive MHD equations on unstructured grids using the spectral difference scheme which is an arbitrarily high-order method that is relatively simple to parallelize. The GEM Magnetic Reconnection problem is used to evaluate whether the Hall MHD equations have been correctly implemented in the solver using the spectral difference method with divergence cleaning (SDDC) algorithm by comparing against the reconnection rates reported in the literature.

Numerical Simulations of Plasmas Generated by Air Explosion

The explosion process of explosives is very complex, often accompanied by phenomena such as luminescence, heat, and electromagnetic radiation. Many studies have shown that the electromagnetic radiation of explosive explosions is related to plasma generated by explosion. In this paper, through the thermal equilibrium ionization model, the plasma model of explosive explosion is established, the calculation formula of plasma state parameters is obtained, and the plasma dynamics program is written. The numerical simulation analysis of the air explosion conditions is carried out by using the plasma dynamics program, and the relationship between the electron density and conductivity of the air explosion plasma and the mass of explosive and distance are summarized.

Coherent microwave scattering for diagnostics of small plasma objects: A review

Measurements of parameters of small-size plasmas are very challenging because many traditional diagnostic approaches cannot be used. Constructive coherent microwave scattering (CMS) offers a convenient diagnostic solution for such small plasmas. This work reviews the development and applications of constructive coherent microwave scattering by the Electric Propulsion and Plasma Laboratory at Purdue University. It presents fundamentals of CMS with an emphasis on Thomson, collisional, and Rayleigh scattering in short, thin, unmagnetized plasma media. Additionally, we review examples of CMS application for diagnostics of temporally resolved plasma dynamics and electron decay, photoionization rates, electron momentum-transfer collision frequencies, and number densities of selective species in gaseous mixtures. These applications are relevant for various research fields including strong field and femtosecond filamentation physics, plasma-assisted ignition and combustion, and combustion and spacecraft electric propulsion diagnostics.

Early dynamics of laser-induced plasma and cavitation bubble in water

The knowledge of early laser-induced plasma and cavitation bubble dynamics are of great importance for a better understanding of underwater laser-induced breakdown spectroscopy. In this work, we investigated the temporal evolution of the plasma and bubble at an early stage before 2 μs, by using fast imaging and shadowgraph techniques. It showed that the initial plasma forms at the laser focal point and then grows up toward the incident laser. As time evolves, the plasma expands simultaneously forward and backward to its peripheral region and a sub-plasma structure can be observed at low laser energies. The plasma has a good pulse-to-pulse repeatability during the laser pulse, while soon after the laser pulse (from 20 to 50 ns), it undergoes a severe pulse-to-pulse fluctuation with an inhomogeneous emission distribution. At longer times, the plasma becomes much stable again. To clarify the interrelationship between the bubble dynamics and the plasma characteristics, shadowgraph images of the bubble were shown and compared with the plasma images. It showed that the morphology of the bubble and its evolution have a good consistency with the morphology of plasma. The bubble radius as a function of time can be well expressed by the equation R=at0.4, and a is a constant related to the laser energy. By inspecting the evolution curve of the bubble wall and the bubble pressure calculated based on the Gilmore model, a transition phase from 20 to 50 ns can be identified between the moving breakdown phase and the thermal expansion phase. In this transition phase, both the bubble expansion speed and the bubble pressure drop drastically that may account for the severe plasma fluctuations as observed in this period. These results of early bubble dynamics provide insights into the critical role of bubble-plasma interaction on the characterization of laser-induced plasma in water.

New Challenges in Nuclear Fusion Reactors: From Data Analysis to Materials and Manufacturing

The construction and operation of the first generation of magnetically controlled nuclear fusion power plants require the development of proper physics and the engineering bases. The analysis of data, recently collected by the actual largest and most important tokamak in the world JET, that has successfully completed his second deuterium and tritium campaign in 2021 (DTE2) with a full ITER like wall main chamber, has provided an important consolidation of the ITER physics basis. Thermonuclear plasmas are highly nonlinear systems characterized by the need of numerous diagnostics to measure physical quantities to guide, through proper control schemes, external actuators. Both modelling and machine learning approaches are required to maximize the physical understanding of plasma dynamics and at the same time, engineering challenges have to be faced. Fusion experiments are indeed extremely hostile environments for plasma facing materials (PFM) and plasma-facing components (PFC), both in terms of neutron, thermal loads and mechanical stresses that the components have to face during either steady operation or off-normal events. Efforts are therefore spent by the community to reach the ultimate goal ahead: turning on the first nuclear fusion power plant, DEMO, by 2050. This editorial is dedicated at reviewing some aspects touched in recent studies developed in this dynamic, challenging project, collected by the special issue titled “New Challenges in Nuclear Fusion Reactors: From Data Analysis to Materials and Manufacturing”.