Ion Heating Research Articles

Earthward‐Tailward Asymmetry of Plasma Temperature in Reconnection Outflow in Earth's Magnetotail

AbstractTo explore the asymmetry in ion and electron heating at Earth's magnetotail at mid‐tail distances ( −30 ), we analyze near‐simultaneous observations of reconnection outflows from two opposite sides of reconnection sites at those distances using Magnetospheric Multiscale (MMS) and Acceleration, Reconnection, Turbulence and Electrodynamics of Moon's Interaction with the Sun (ARTEMIS) data. We report a pronounced temperature asymmetry between the earthward and tailward reconnection outflows. The asymmetry is more significant for electrons than for ions: Earthward moving ions are only three times hotter than tailward ones, but earthward moving electrons are 5–20 times hotter than tailward ones. The closed field‐line topology on the earthward side of the reconnection region, as opposed to the open topology on the tailward side, is likely a critical contributor to this asymmetry. These findings cast light on the underlying mechanisms of particle heating and energization in magnetotail reconnection, highlighting the significant role of Earth's dipolar magnetic field. This study offers insights for refining magnetic reconnection models, emphasizing the importance of incorporating realistic magnetic field topologies to accurately simulate the heating and energization processes observed in space plasma environments.

Multifluid Theory of Electrostatic Ion Cyclotron Waves in Partially Ionized Plasmas

Partially ionized plasmas universally exist in various astrophysical environments, such as the solar atmosphere and the E region of the ionosphere. In these contexts, the existence and propagation of waves in plasmas could be significantly influenced by effects of weakly ionized plasma (e.g., ion–neutral collisions). In this work, we investigate electrostatic ion cyclotron (EIC) waves in partially ionized plasmas based on the multifluid model with adiabatic electrons. Two distinct branches of EIC waves coexist in partially ionized plasmas: one branch is the conventional EIC waves; the other branch propagates around the “effective ion cyclotron frequency” which originates from self-consistent ion–neutral collisions. Furthermore, theoretical predictions in the new branch of EIC waves are qualitatively consistent with laboratory observations. In addition, a comparison between our theory and the previous work is also performed. This work can aid in understanding the acceleration and transverse heating of ions in partially ionized astrophysical plasmas where the ion–neutral collisions are frequent.

PIC simulation of a nonoscillatory perturbation on a subcritical fast magnetosonic shock wave

Abstract We use a two-dimensional particle-in-cell (PIC) simulation to study the propagation of subcritical fast magnetosonic shocks in electron-nitrogen plasma and their stability against an initial deformation. A slab of dense plasma launches two planar blast waves into a surrounding ambient plasma, which is permeated by a magnetic field that points out of the simulation box and is spatially uniform at the start of the simulation. One shock propagates into a spatially uniform ambient plasma. This reference shock has a Mach number of 1.75, and the heating of ions only along the shock normal compresses the ions that cross the shock to twice the upstream density. Drift instabilities lead to rapidly growing electron-cyclotron harmonic waves ahead of the location where the shock's density overshoot peaks, and to slowly growing lower-hybrid waves with a longer wavelength behind it. The second shock wave enters a perturbation layer that deforms it into a sine shape. Once the shock leaves the perturbation layer, the deformation is weakly damped and non-oscillatory, and the shock remains stable. Even without an external perturbation, and for the plasma parameters considered here, drift instabilities will cause ripples in the shock wave. These instabilities lead to a spatially and temporally varying compression of the plasma that crosses the shock.

Asymmetric thermo-electro-osmotic responses in charged conical nanochannels

Nanofluidic thermo-electric and thermo-osmotic responses have attracted increasing attention for their potential in low-grade heat energy recovery. However, the synergistic influence of geometric asymmetry and electrolyte physical properties on these responses is often overlooked. This study investigates asymmetric thermo-electro-osmotic responses by systematically varying parameters such as conicity, Debye length, and surface charge density within conical nanochannels. We employ the numerical model coupling extended Poisson-Nernst-Planck-Navier-Stokes and energy equations to simulate ion transport, fluid flow, and heat transfer. Results demonstrate pronounced asymmetries in thermo-electric and thermo-osmotic responses between forward and backward nanochannel configurations. The short-circuit current strongly depends on conicity, Debye length, and surface charge density, while the Seebeck coefficient is primarily influenced by Debye length, and surface charge density. Thermo-osmotic flow is significantly affected by all parameters, with flow direction reversals observed under specific conditions. The study highlights that Debye length and surface charge density significantly influence both thermo-electric and -osmotic responses, whereas geometric asymmetry predominantly affects thermo-osmotic response. This study provides a valuable framework for understanding fundamental thermal-driven transport phenomena at the nanoscale. The observed synergistic effects between various parameters offer insights for optimizing thermo-electric energy conversion and fluidic control within conical nanochannels, paving the way for innovative applications in nanofluidic technology and energy recovery systems.

Electron-only Reconnection and Ion Heating in 3D3V Hybrid-Vlasov Plasma Turbulence

We perform 3D3V hybrid-Vlasov simulations of turbulence with quasi-isotropic, compressible injection near ion scales to mimic the Earth’s magnetosheath plasma, and investigate the novel electron-only reconnection, recently observed by NASA’s Magnetospheric Multiscale mission, and its impact on ion heating. Retaining electron inertia in the generalized Ohm's law enables collisionless magnetic reconnection. Spectral analysis shows a shift from kinetic Alfvén waves to inertial kinetic Alfvén and inertial whistler waves near electron scales. To distinguish the roles of inertial scale and gyroradius (d i and ρ i), three ion beta (β i = 0.25, 1, 4) values are studied. Ion-electron decoupling increases with β i, as ions become less mobile when the injection scale is closer to ρ i than d i, highlighting the role of ρ i in achieving an electron magnetohydrodynamic regime at sub-ion scales. This regime promotes electron-only reconnection in turbulence with small-scale injection at β i ≳ 1. We observe significant ion heating even at large β i, with Q i/ϵ ≈ 69%, 91%, and 96% at β i = 0.25, 1, and 4, respectively. While ion heating is anisotropic at β i ≤ 1 (T i,⊥ > T i,∥), it is marginally anisotropic at β i > 1 (T i,⊥ ≳ T i,∥). Our results show ion turbulent heating in collisionless plasmas is sensitive to the separation between injection scales (λ inj) and ρ i, β i, and finite-k ∥ effects, necessitating further investigation for accurate modeling. These findings have implications for other collisionless astrophysical environments, like high-β plasmas in intracluster medium, where processes such as microinstabilities or shocks may inject energy near ion-kinetic scales.

The Effects of Turbulence on Heliosheath Ions and Implications for Energetic Neutral Atoms

The distribution of ions in the heliosheath—the region between the heliospheric termination shock and the heliopause—is important for understanding remote observations of energetic neutral atoms (ENAs). The ion distributions were estimated previously based on hybrid simulations of the heating and evolution of solar wind and interstellar pickup ions across the solar wind termination shock, but these estimates only provide the distributions near the shock. In this work, we use self-consistent hybrid kinetic simulations to investigate the effects of turbulence on ion distributions in the heliosheath. The simulations are compared against Voyager observations, constraining the feasible amplitude and compressibility of turbulence. We find that the heating due to turbulent dissipation can lead to a significant increase in the temperature of thermal solar wind ions. Both turbulent velocity fluctuations and the heating of solar wind ions increase the charge-exchange source for ENAs at low energies (around 100 eV), where current ENA models underpredict observations by more than an order of magnitude. However, the effects of turbulence are likely not strong enough to fully explain these discrepancies.

Observations of Significant Ion Energy Outflows Associated With Cusp Ion Outflows and the Role of Poynting Flux as an Energy Source

AbstractThe cusp ion outflow on Earth has been extensively studied for decades. However, the energy flux associated with the ion outflow, which is of equivalent importance to the number flux or mass flux, has been rarely studied. Here, we present the first systematic study on the cusp ion energy flux and the energy budget along the cusp magnetic flux tube using quasi‐conjunction observations from the Polar and Fast satellites. Significant ion energy fluxes (several 10s mW/ mapped to 100 km altitude) away from the Earth are frequently observed in the mid‐altitude cusp (3–6 Re). Observations at low altitudes (2 Re) show that the upward ion energy flux is associated with ion outflows originating from the ionosphere. In addition, we show that the ion outflows experience intense energization well above the ionosphere. The only possible energy source for this energization is the earthward Poynting flux in mid‐altitude. The electrons are accelerated downward and are another energy sink (not source). The altitude profile of the energy fluxes suggests a transition region between 2 and 4 Re where both ion heating and electron acceleration primarily occur and where significant Poynting flux is dissipated. Analysis of the E/B ratio shows that the Poynting flux is carried by both the Alfven waves and quasi‐static structures. The above results place important constraints on possible local ion energization mechanisms.

Abiotic Stress in Cotton: Insights into Plant Responses and Biotechnological Solutions

Climate change has rapidly increased incidences of frequent extreme abiotic stresses, such as heat, drought, salinity, and waterlogging. Each of these stressors negatively affects the cotton crop (Gossypium spp.) and results in significant yield decreases. Every stressful event causes specific changes in the metabolism and physiology of plants, which are linked to complex molecular alterations. Understanding the molecular mechanisms that regulate a plant’s response to stress is essential to developing stress-resistant cotton varieties that can withstand various stress factors. Gene expressions in response to multiple stresses have been studied and mapped. These genes include ion transporters and heat shock proteins, which are vital to allowing adaptive responses. These approaches showed the ability to employ advanced genome sequencing and multi-omics techniques to identify dynamic gene expression patterns and elucidate intricate regulatory networks. Using genetic variation in combination with molecular techniques, it would be possible to generate stress-resilient cotton varieties that would enable sustainable cotton output in the face of abiotic stresses. Here, we reviewed the effects of major abiotic stressors on cotton plants, such as heat, salinity, drought, heavy metals, and waterlogging. We also examine the vast network of proteins, genes, and stress-sensitive signaling pathways that help cotton tolerate abiotic stress.

Summary of the 11th Conference on Magnetic Confined Fusion Theory and Simulation

This conference report summarizes recent progress in plasma theory and simulation that was presented in contributed papers and discussions at the 11th Conference on Magnetic Confined Fusion Theory and Simulation (CMCFTS) held in Chengdu, China, 27–30 October, 2023. Progress in various fields has been achieved. For example, results on zonal flow generation by mode coupling, simulations of the key physics of divertor detachment, energetic particle effects on magnetohydrodynamic (MHD) modes in addition to ion- and electron-scale turbulence, physics of edge coherent modes and edge-localized modes, and the optimization of ion heating schemes as well as confinement scenarios using advanced integrated modeling are presented at the conference. In this conference, the scientific research groups were organized into six categories: (a) edge and divertor physics; (b) impurity, heating, and current drive; (c) energetic particle physics; (d) turbulent transport; (e) MHD instability; and (f) integrated modeling and code development. A summary of the highlighted progress in these working groups is presented.

Three-dimensional transient modeling and simulation of high-current multi-component vacuum arc under transverse magnetic field

This paper established a three-dimensional transient MHD model of a high-current multi-component vacuum arc under a transverse magnetic field (TMF) to simulate the dynamic characteristics of plasma parameters in the uniform motion mode of a contracted arc under lower TMF. The model considered the ionization–recombination process among the multi-components and used the P1 model to calculate the thermal radiation losses. The simulation results show that the ion number density of The TMF arc is on the order 1024 m−3, and when the plasma temperature reaches 4.3 eV, the effects of thermal radiation cannot be ignored. The majority of the plasma from the cathode and anode undergo acceleration and then deceleration before finally intersecting at the center of the arc, creating an extreme value of the density and the pressure, where a large amount of kinetic energy is converted into internal energy. The arc will be bent and deflected under the action of the ampere force, and the deflection of Cu3+ is especially obvious. The TMF affects the collision strength of the jet in the intersection area and, thus, affects the ion density, in which the change of Cu2+ is dominant. The ionization–recombination process of ions is mainly determined by the electron temperature, which is affected by the arc current. As the current decreases, Cu1+ increases and Cu2+ and Cu3+ decrease, and the change in ionization rate is the main reason for the change in the proportion of each ionic component; the heat flux density to the anode and cathode is on the order of 1010 W/m2, which heats the front of the arc and ensures the stable movement of the arc. Meanwhile, due to the high ion number density, the ion heat flux accounts for the main part of the heat flux, and the anode exhibits a higher proportion of ion heat flux.

Magnetosonic waves in the Martian ionosphere driven by upstream proton cyclotron waves: Two-point observations by MAVEN and Mars Express

Recent observations from the Mars Atmosphere and Volatile EvolutioN (MAVEN) and Mars Express (MEX) spacecraft have suggested that pressure pulses originating from upstream proton cyclotron waves (PCWs) can “ring” the Martian magnetopause at the same frequency and drive magnetosonic waves in the upper ionosphere of Mars, thereby transporting energy from the solar wind into the ionosphere. However, the limitation of single-spacecraft measurements prevents simultaneous observations of the driver and response of this “ringing” process of the Martian magnetosphere. Here we utilize two-point measurements from MAVEN and MEX to characterize the ringing probability at which upstream PCWs drive compressional fluctuations in the ionospheric magnetic field. We develop an algorithm to identify PCW-driven magnetosonic waves in the upper ionosphere of Mars from the two-point magnetic field data. The derived ringing probability is higher on the dayside, outside strong crustal magnetic fields, and under high solar wind density conditions. We also show that the median power of dayside ionospheric magnetic field fluctuations is enhanced by a factor of ∼2 at corresponding frequencies in the presence of upstream PCWs compared to the median power in the absence of upstream PCWs. These results demonstrate the prevalence of energy deposits into the dayside Martian ionosphere from the solar wind mediated by the PCW-driven ringing of the magnetosphere. Future studies, possibly with new multi-point observations, should address the detailed processes of wave propagation and energy transport through the system and the long-term impact of this chain of processes on the planetary ion heating in the ionosphere and atmospheric loss from Mars.

Reduction or Enhancement of Stellarator Turbulence by Impurities

A systematic study of the impact of impurities on the turbulent heat fluxes is presented for the stellarator Wendelstein 7-X (W7-X) and, for comparison, the Large Helical Device and ITER. By means of nonlinear multispecies gyrokinetic simulations, it is shown that impurities, depending on the sign of their density gradient, can significantly enhance or reduce turbulent ion heat losses. For the relevant scenario of turbulence reduction, an optimal impurity concentration that minimizes the ion heat diffusivity emerges as a universal feature. This result demonstrates the potential of impurities for controlling turbulence and accessing enhanced confinement regimes in fusion plasmas and, in particular, in W7-X. Published by the American Physical Society 2024

The transient thermomechanical behavior of solid oxide fuel cell during the operation condition variations

Solid oxide fuel cell (SOFC) is a promising energy conversion technology with high efficiency, low emission and fuel flexibility. In practical and dynamic operation, the thermal stress at the interface of different layers can cause cracks or even thermal mechanical failure. This study investigates the transient thermomechanical behavior of SOFC under changing operating conditions. A three-dimensional dynamical model is developed by considering the electrochemical reaction, electron and ion transport, gas diffusion, heat and momentum transfer, and thermal-mechanical interaction. When the operating conditions vary, the time it takes for the thermal stress to reach a steady state is consistent with the time required for heat transfer to stabilize. The electrolyte layer is subjected to the maximum gradient of thermal stress change after the voltage is switched. With an increase in the inlet gas temperature, the change in the cell temperature can be divided into two stages: a rapid rise and a slow rise. Appropriate prolongation of the temperature rise time is beneficial to the thermal stability of SOFC. The increase in inlet gas temperature changes the distribution trend of cell temperature and thermal stress, while the 50 % reduction in cathode and anode inlet velocity changes its magnitude.

Modeling of ion cyclotron resonance frequency heating of proton-boron plasmas in EHL-2 spherical tokamak

Ion cyclotron resonance heating (ICRH) stands out as a widely utilized and cost-effective auxiliary method for plasma heating, bearing significant importance in achieving high-performance discharges in p-11B plasmas. In light of the specific context of p-11B plasma in the EHL-2 device, we conducted a comprehensive scan of the fundamental physical parameters of the antenna using the full-wave simulation program TORIC. Our preliminary result indicated that for p-11B plasma, optimal ion heating parameters include a frequency of 40 MHz, with a high toroidal mode number like to heat the majority H ions. In addition, we discussed the impact of concentration of minority ion species on ion cyclotron resonance heating when 11B serves as the heavy minority species. The significant difference in charge-to-mass ratio between boron and hydrogen ions results in a considerable distance between the hybrid resonance layer and the tow inverted cyclotron resonance layer, necessitating a quite low boron ion concentration to achieve effective minority heating. We also considered another method of direct heating of hydrogen ions in the presence of boron ion minority. It is found that at appropriate boron ion concentrations ( ), the position of the hybrid resonance layer approaches that of the hydrogen ion cyclotron resonance layer, thereby altering the polarization at this position and significantly enhancing hydrogen ion fundamental absorption.

Recent advance progress of HL-3 experiments

Abstract Since the first plasma realized in 2020, a series of key systems on HL-3 (known as HL-2M before) tokamak have been equipped/upgraded, including in-vessel components (the first wall, lower divertor, and toroidal cryogenic/water-cooling/baking/glow discharge systems, etc.), auxiliary heating system of 11 MW, and 28 diagnostic systems (to measure the plasma density, electron temperature, radiation, magnetic field, etc.). Magnet field systems were commissioned firstly for divertor plasma discharges. During the 2nd experimental campaign of HL-3 tokamak, several great progresses have been achieved. Firstly, the successful operation with plasma current larger than 1 MA was achieved under a divertor configuration. Secondly, the advanced divertor concept with two distinct snowflake configurations was realized. It is found that the distribution of ion saturation current and heat flux on bottom plate becomes wide due to magnetic surface expansion, demonstrating the advantage of such configuration in the heat flux mitigation. In addition, using the combination of NBI, ECRH and LHCD, the standard sawtoothing high confinement mode of megampere plasma was firstly accessed on the HL-3. The successful commissioning of HL-3 is beneficial for the initial operation of ITER.

Neoclassical transport computations in non-isothermal tokamak plasmas

Neoclassical transport in a non-isothermal plasma in which each plasma species has different equilibrium temperature is investigated by solving the drift kinetic equation with a Fokker–Planck (FP) collision operator in a circular tokamak model. Since it is known that a linearized FP operator does not have a self-adjoint property in a non-isothermal plasma, approximate models are developed for comparison to intend to have the self-adjoint property in the non-isothermal case. To achieve this, we set a common temperature that the system should reach after a long time, and the individual temperature of each particle species is expressed by a parameter to measure a shift of the individual temperature from the common one. Then, both the Vlasov part and the collision term of the kinetic equation are expanded around the common temperature, taking the temperature shift parameter up to the first order. It is found that the lowest order collision term of expansion preserves the self-adjointness while the first-order, nonlinear FP term does not. A large difference of the ion heat neoclassical transport is found in comparison between the developed models with and without the self-adjointness and the original FP collision term in the non-isothermal plasma, especially in a strong temperature equilibration regime, showing that a contribution of the collision term without the self-adjointness seems to be significant. Furthermore, when an impurity species is included, the result is complicated where the usual enhancement in the main ion particle transport coefficient, due to the impurity effect, is rather suppressed with the increase in the ion heat transport coefficients by the non-isothermal effects.

Overview of Large Helical Device experiments of basic plasma physics for solving crucial issues in reaching burning plasma conditions

Abstract Recently, experiments on basic plasma physics issues for solving future problems in fusion energy have been performed on a Large Helical Device. There are several problems to be solved in future devices for fusion energy. Emerging issues in burning plasma are: alpha-channeling (ion heating by alpha particles), turbulence and transport in electron dominant heating helium ash exhaust, reduction of the divertor heat load. To solve these problems, understanding the basic plasma physics of (1) wave–particle interaction through (inverse) Landau damping, (2) characteristics of electron-scale (high-k) turbulence, (3) ion mixing and the isotope effect, and (4) turbulence spreading and detachment, is necessary. This overview discusses the experimental studies on these issues and turbulent transport in multi-ion plasma and other issues in the appendix.

Scaling of Ion Bulk Heating in Magnetic Reconnection Outflows for the High-Alfvén-speed and Low-β Regime in Earth’s Magnetotail

We survey 20 reconnection outflow events observed by Magnetospheric MultiScale in the low-β and high-Alfvén-speed regime of the Earth’s magnetotail to investigate the scaling of ion bulk heating produced by reconnection. The range of inflow Alfvén speeds (800–4000 km s−1) and inflow ion β (0.002–1) covered by this study is in a plasma regime that could be applicable to the solar corona and flare environments. We find that the observed ion heating increases with increasing inflow (upstream) Alfvén speed, V A, based on the reconnecting magnetic field and the upstream plasma density. However, ion heating does not increase linearly as a function of available magnetic energy per particle, miVA2 . Instead, the heating increases progressively less as miVA2 rises. This is in contrast to a previous study using the same data set, which found that electron heating in this high-Alfvén-speed and low-β regime scales linearly with miVA2 , with a scaling factor nearly identical to that found for the low-V A and high-β magnetopause. Consequently, the ion-to-electron heating ratio in reconnection exhausts decreases with increasing upstream V A, suggesting that the energy partition between ions and electrons in reconnection exhausts could be a function of the available magnetic energy per particle. Finally, we find that the observed difference in ion and electron heating scaling may be consistent with the predicted effects of a trapping potential in the exhaust, which enhances electron heating, but reduces ion heating.

The effects of finite electron inertia on helicity-barrier-mediated turbulence

Understanding the partitioning of turbulent energy between ions and electrons in weakly collisional plasmas is crucial for the accurate interpretation of observations and modelling of various astrophysical phenomena. Many such plasmas are ‘imbalanced’, wherein the large-scale energy input is dominated by Alfvénic fluctuations propagating in a single direction. In this paper, we demonstrate that when strongly-magnetised plasma turbulence is imbalanced, nonlinear conservation laws imply the existence of a critical value of the electron plasma beta (the ratio of the thermal to magnetic pressures) that separates two dramatically different types of turbulence in parameter space. For betas below the critical value, the free energy injected on the largest scales is able to undergo a familiar Kolmogorov-type cascade to small scales where it is dissipated, heating electrons. For betas above the critical value, the system forms a ‘helicity barrier’ that prevents the cascade from proceeding past the ion Larmor radius, causing the majority of the injected free energy to be deposited into ion heating. Physically, the helicity barrier results from the inability of the system to adjust to the disparity between the perpendicular-wavenumber scalings of the free energy and generalised helicity below the ion Larmor radius; restoring finite electron inertia can annul, or even reverse, this disparity, giving rise to the aforementioned critical beta. We relate this physics to the ‘dynamic phase alignment’ mechanism (that operates under yet lower beta conditions and in pair plasmas), and characterise various other important features of the helicity barrier, including the nature of the nonlinear wavenumber-space fluxes, dissipation rates, and energy spectra. The existence of such a critical beta has important implications for heating, as it suggests that the dominant recipient of the turbulent energy, ions or electrons, can depend sensitively on the characteristics of the plasma at large scales.

Azimuthal ion dynamics at the inner pole of an axisymmetric Hall thruster

The azimuthal dynamics of ions along the inner pole of a Hall thruster with a centrally mounted cathode and a magnetic shielding topography are experimentally investigated. A time-averaged laser-induced fluorescence diagnostic is implemented to characterize the azimuthal ion velocity distribution, and its moments are computed numerically to infer bulk rotation speed and ion temperature. It is found that the time-averaged ion swirl velocity grows to 2 km/s in the near-pole region, and the cathode ions exhibit ion temperatures in the azimuthal direction approaching 8 eV. Both of these quantities exceed the speeds and temperatures anticipated from classical acceleration and heating. Time-resolved laser-induced fluorescence is then employed to investigate the role of plasma fluctuations in driving the time-averaged ion properties. Semicoherent fluctuations at 90 kHz are observed in the ion velocity distribution and its associated moments. These oscillations are correlated with the gradient-driven anti-drift wave, which propagates azimuthally in the near-field cathode plume. Quasilinear theory is used to construct a 1D model for acceleration and heating of the ion population as a result of the anti-drift mode. This approach demonstrates qualitative agreement with the time-averaged ion velocity and temperature, suggesting that the anti-drift mode may be a dominant driver of azimuthal ion acceleration and heating in front of the cathode keeper and the inner half of the inner front pole cover. These results are discussed in terms of their relevance to the erosion of thruster surfaces in the near-field cathode plume.