Transport Of Energetic Particles Research Articles

Cross-scale interaction between microturbulence and meso-scale reversed shear Alfvén eigenmodes in DIII-D plasmas

This paper reports global nonlinear gyrokinetic simulations that couple meso-scale reversed shear Alfvén eigenmodes (RSAEs) driven by energetic particles (EPs) and ion temperature gradient (ITG) microturbulence driven by thermal plasma, using equilibrium and profiles from DIII-D discharge #159243. In simulations focusing only on the ITG, electrostatic ITG drives a huge thermal ion heat transport, which is reduced by a factor of 10 to a level close to the experimental value in electromagnetic simulation due to finite β effect. In the simulations coupling the RSAE and ITG, ITG can scatter the resonant EP nonlinearly trapped by the RSAE and damp the zonal flows generated by the RSAE. The regulation of the RSAE by the ITG greatly reduces the initial saturation amplitude of the RSAE but increases the RSAE amplitude and associated EP transport to experimental levels in the quasi-steady state. The RSAE effects on the ITG, specifically the stronger zonal flows generated by the RSAE and the RSAE frequency modulation of the ITG-induced thermal ion heat transport, in turn, leads to a reduction of the thermal ion heat transport by more than a factor of 2 . For a stronger background ITG, the regulation of the RSAE by the ITG is stronger, while the RSAE effects on the ITG are weaker. This work highlights the importance of cross-scale coupling in the dynamics of the AE turbulence and EP transport.

Effects of plasma nonuniformity on toroidal Alfvén eigenmode nonlinear decay

The parametric decay of toroidal Alfvén eigenmode (TAE) in nonuniform plasmas is investigated using nonlinear gyrokinetic equation. It is found that, the plasma nonuniformity not only significantly enhances the nonlinear coupling cross-section, but also qualitatively modifies the decay process. Specifically, the condition for spontaneous decay becomes the toroidal mode number of the sideband TAE being higher than that of the pump TAE, instead of the frequency of the sideband TAE being lower than the pump TAE in uniform plasmas. The consequences on TAE saturation and energetic particle transport are also discussed.

Transport of energetic particles in turbulent space plasmas: pitch-angle scattering, telegraph, and diffusion equations

Introduction: In this article, we revisit the pitch-angle scattering equation describing the propagation of energetic particles through magnetized plasma. In this case, solar energetic particles and cosmic rays interact with magnetohydrodynamic turbulence and experience stochastic changes in the pitch-angle. Since this happens over an extended period of time, a pitch-angle isotropization process occurs, leading to parallel spatial diffusion. This process is described well by the pitch-angle scattering equation. However, the latter equation is difficult to solve analytically even when considering special cases for the scattering coefficient.Methods: In the past, a so-called subspace approximation was proposed, which has important applications in the theory of perpendicular diffusion. Alternatively, an approach based on the telegraph equation (also known as telegrapher’s equation) has been developed. We show that two-dimensional subspace approximation and the description based on the telegraph equation are equivalent. However, it is also shown that the obtained distribution functions contain artifacts and inaccuracies that cannot be found in the numerical solution to the problem. Therefore, an N-dimensional subspace approximation is proposed corresponding to a semi-analytical/semi-numerical approach. This is a useful alternative compared to standard numerical solvers.Results and Discussion: Depending on the application, the N-dimensional subspace approximation can be orders of magnitude faster. Furthermore, the method can easily be modified so that it can be used for any pitch-angle scattering equation.

Wave–particle interactions in tokamaks

Transport consequences of the wave–particle interactions in the quasilinear plateau (QP) regime are presented. Eulerian approach is adopted to solve the drift kinetic equation that includes the physics of the nonlinear trapping (NT) and QP regimes. The localization of the perturbed distribution simplifies the test particle collision operator. It is shown that a mirror force like term responsible for the flattening of the distribution in the NT regime is subdominant in the QP regime, and controls the transition between these two regimes. Transport fluxes, flux-power relation, and nonlinear damping or growth rate are all calculated. There is no explicit collision frequency dependence in these quantities; however, the width of the resonance does. Formulas that join the asymptotic results of these two regimes to facilitate thermal and energetic particle transport, and nonlinear wave evolution of a single mode are presented.

Nonlinear dynamics of the reversed shear Alfvén eigenmode in burning plasmas

In a tokamak fusion reactor operated at steady state, the equilibrium magnetic field is likely to have reversed shear in the core region, as the noninductive bootstrap current profile generally peaks off-axis. The reversed shear Alfvén eigenmode (RSAE) as a unique branch of the shear Alfvén wave in this equilibrium, can exist with a broad spectrum in wavenumber and frequency, and be resonantly driven unstable by energetic particles (EP). After briefly discussing the RSAE linear properties in burning plasma condition, we review several key topics of the nonlinear dynamics for the RSAE through both wave-EP resonance and wave-wave coupling channels, and illustrate their potentially important role in reactor-scale fusion plasmas. By means of simplified hybrid MHD-kinetic simulations, the RSAEs are shown to have typically broad phase space resonance structure with both circulating and trapped EP, as results of weak/vanishing magnetic shear and relatively low frequency. Through the route of wave-EP nonlinearity, the dominant saturation mechanism is mainly due to the transported resonant EP radially decoupling with the localized RSAE mode structure, and the resultant EP transport generally has a convective feature. The saturated RSAEs also undergo various nonlinear couplings with other collective oscillations. Two typical routes as parametric decay and modulational instability are studied using nonlinear gyrokinetic theory, and applied to the scenario of spontaneous excitation by a finite amplitude pump RSAE. Multiple RSAEs could naturally couple and induce the spectral energy cascade into a low frequency Alfvénic mode, which may effectively transfer the EP energy to fuel ions via collisionless Landau damping. Moreover, zero frequency zonal field structure could be spontaneously excited by modulation of the pump RSAE envelope, and may also lead to saturation of the pump RSAE by both scattering into stable domain and local distortion of the continuum structure.

Solar Energetic Particle and the Heliospheric Current Sheet

The effect of the heliospheric current sheet (HCS) on the propagation of solar energetic particles (SEPs) remains poorly known. In this study we address this question by surveying energetic (∼2.0–9.6 MeV nucleon–1) helium data acquired by the energetic particle acceleration, composition, and transport (EPACT) sensor on board the Wind spacecraft. A superposed epoch analysis of 319 HCS crossings made by Wind reveals a sharp drop in the SEP fluxes at the HCS for the low-energy channels and little change across the HCS for the high-energy channels. To help understand the statistical result, we studied a total of 15 SEP flux dropout (a decrease of ∼50% or more) events that coincided with the crossing of the HCS. One of the common features of these SEP events is that they were initiated in the western hemisphere but far away from the longitude of HCS crossings, suggesting that the source of SEPs was well connected initially but was cut off later after Wind moved to the opposite hemisphere (e.g., HCS crossing). Further analysis of the events suggests that the percentage of flux dropouts decreases with increasing energy. It is suggested that a strong scattering of MeV helium may have occurred as the particle gyroradius is comparable to the thickness of the current sheet. This study clearly provides solid evidence for the HCS as a barrier to suppressing SEP flux of MeV energies from the onset hemisphere to the other.

Identification of the weak-to-strong transition in Alfvénic turbulence from space plasma

AbstractPlasma turbulence is a ubiquitous dynamical process that transfers energy across many spatial and temporal scales in astrophysical and space plasma systems. Although the theory of anisotropic magnetohydrodynamic (MHD) turbulence has successfully described natural phenomena, its core prediction of an Alfvénic transition from weak to strong MHD turbulence when energy cascades from large to small scales has not been observationally confirmed. Here we report evidence for the Alfvénic weak-to-strong transition in small-amplitude, turbulent MHD fluctuations in Earth’s magnetosheath using data from four Cluster spacecraft. Our observations demonstrate the universal existence of strong turbulence accompanied by weak turbulent fluctuations on large scales. Moreover, we find that the nonlinear interactions of MHD turbulence are crucial to the energy cascade, as they broaden the cascade directions and fluctuating frequencies. The observed connection between weak and strong MHD turbulence systems may be present in many astrophysical environments, such as star formation, energetic particle transport, turbulent dynamos, and solar corona or solar wind heating.

Conceptual Design of the Tail Research EXperiment in Space Plasma Environment Research Facility (SPERF-TREX)

Abstract The Space Plasma Environment Research Facility (SPERF) for ground simulation of space plasma environment is a key component of Space Environment Simulation Research Infrastructure (SESRI), a major national science and technology infrastructure for fundamental researches. It is designed to investigate outstanding issues in space plasma environment, such as energetic particles acceleration, transport, and interaction with electromagnetic waves, as well as magnetic reconnection processes, in the magnetospheric plasmas, etc. Tail-Research EXperiment (TREX) is part of the SPERF for laboratory studies of space physics relevant to magnetic reconnection, dipolarization and hydromagnetic waves excitation in the magnetotail. SPERF-TREX is designed to carry out three types of experiments: the tail plasmoid for magnetic reconnection, dipolarization front formation, and magnetohydrodynamic waves excited by high speed plasma jet. In this paper, the scientific goals and three scenarios of SPERF-TREX for typical processes in space plasmas are presented, and experimental plans for SPERF-TREX are also reviewed, together with plasma sources applied to generate the plasma with desired parameters and various magnetic configurations.

Parallel Diffusion Coefficient of Energetic Charged Particles in the Inner Heliosphere from the Turbulent Magnetic Fields Measured by Parker Solar Probe

Diffusion coefficients of energetic charged particles in turbulent magnetic fields are a fundamental aspect of diffusive transport theory but remain incompletely understood. In this work, we use quasi-linear theory to evaluate the spatial variation of the parallel diffusion coefficient κ ∥ from the measured magnetic turbulence power spectra in the inner heliosphere. We consider the magnetic field and plasma velocity measurements from Parker Solar Probe made during Orbits 5–13. The parallel diffusion coefficient is calculated as a function of radial distance from 0.062 to 0.8 au, and the particle energy from 100 keV to 1 GeV. We find that κ ∥ increases exponentially with both heliocentric distance and energy of particles. The fluctuations in κ ∥ are related to the episodes of large-scale magnetic structures in the solar wind. By fitting the results, we also provide an empirical formula of κ ∥ = (5.16 ± 1.22) × 1018 r 1.17±0.08 E 0.71±0.02 (cm2 s−1) in the inner heliosphere, which can be used as a reference in studying the transport and acceleration of solar energetic particles as well as the modulation of cosmic rays.

Energetic particle acceleration and transport with the novel Icarus+PARADISE model

With the rise of satellites and mankind's growing dependence on technology, there is an increasing awareness of space weather phenomena related to high-energy particles. Shock waves driven by coronal mass ejections (CMEs) and corotating interaction regions (CIRs) occasionally act as potent particle accelerators, generating hazardous solar energetic particles (SEPs) that pose risks to satellite electronics and astronauts. Numerical simulation tools capable of modelling and predicting large SEP events are thus highly demanded. We introduce the new Icarus+PARADISE model as an advancement of the previous EUHFORIA+PARADISE model. Icarus, based on the MPI-AMRVAC framework, is a three-dimensional magnetohydrodynamic code that models solar wind configurations from 0.1 au onwards, encompassing transient structures like CMEs or CIRs. Differing from EUHFORIA's uniform-only grid, Icarus incorporates solution adaptive mesh refinement (AMR) and grid stretching. The particle transport code PARADISE propagates energetic particles as test particles through these solar wind configurations by solving the focused transport equation in a stochastic manner. We validate our new model by reproducing EUHFORIA+PARADISE results. This is done by modelling the acceleration and transport of energetic particles in a synthetic solar wind configuration containing an embedded CIR. Subsequently, we illustrate how the simulation results vary with grid resolution by employing different levels of AMR.. The resulting intensity profiles illustrate increased particle acceleration with higher levels of AMR in the shock region, better capturing the effects of the shock.

Element abundance and the physics of solar energetic particles

The acceleration and transport of solar energetic particles (SEPs) cause their abundance, measured at a constant velocity, to be enhanced or suppressed as a function of the magnetic rigidity of each ion, and hence, of its atomic mass-to-charge ratio of A/Q. Ion charges, in turn, depend upon the source electron temperature. In small “impulsive” SEP events, arising from solar jets, acceleration during magnetic reconnection causes steep power-law abundance enhancements. These impulsive SEP events can have 1,000-fold enhancements of heavy elements from sources at ∼2.5 MK and similar enhancements of 3He/4He and of streaming electrons that drive type-III radio bursts. Gamma-ray lines show that solar flares also accelerate 3He-rich ions, but their electrons and ions remain trapped in magnetic loops, so they dissipate their energy as X-rays, γ-rays, heat, and light. “Gradual” SEPs accelerated at shock waves, driven by fast coronal mass ejections (CMEs), can show power-law abundance enhancements or depressions, even with seed ions from the ambient solar corona. In addition, shocks can reaccelerate seed particles from residual impulsive SEPs with their pre-existing signature heavy-ion enhancements. Different patterns of abundance often show that heavy elements are dominated by a source different from that of H and He. Nevertheless, the SEP abundance, averaged over many large events, defines the abundance of the corona itself, which differs from the solar photosphere as a function of the first ionization potential (FIP) since ions, with FIP <10 eV, are driven upward by forces of electromagnetic waves, which neutral atoms, with FIP >10 eV, cannot feel. Thus, SEPs provide a measurement of element abundance in the solar corona, distinct from solar wind, and may even better define the photosphere for some elements.

Saturation of Fishbone Instability by Self-Generated Zonal Flows in Tokamak Plasmas.

Gyrokinetic simulations of the fishbone instability in DIII-D tokamak plasmas find that self-generated zonal flows can dominate the nonlinear saturation by preventing coherent structures from persisting or drifting in the energetic particle phase space when the mode frequency down-chirps. Results from the simulation with zonal flows agree quantitatively, for the first time, with experimental measurements of the fishbone saturation amplitude and energetic particle transport. Moreover, the fishbone-induced zonal flows are likely responsible for the formation of an internal transport barrier that was observed after fishbone bursts in this DIII-D experiment. Finally, gyrokinetic simulations of a related ITER baseline scenario show that the fishbone induces insignificant energetic particle redistribution and may enable high performance scenarios in ITER burning plasma experiments.

Small-amplitude Compressible Magnetohydrodynamic Turbulence Modulated by Collisionless Damping in Earth’s Magnetosheath: Observation Matches Theory

Plasma turbulence is a ubiquitous dynamical process that transfers energy across many spatial and temporal scales and affects energetic particle transport. Recent advances in the understanding of compressible magnetohydrodynamic (MHD) turbulence demonstrate the important role of damping in shaping energy distributions on small scales, yet its observational evidence is still lacking. This study provides the first observational evidence of substantial collisionless damping (CD) modulation on the small-amplitude compressible MHD turbulence cascade in Earth’s magnetosheath using four Cluster spacecraft. Based on an improved compressible MHD decomposition algorithm, turbulence is decomposed into three eigenmodes: incompressible Alfvén modes and compressible slow and fast (magnetosonic) modes. Our observations demonstrate that CD enhances the anisotropy of compressible MHD modes because CD has a strong dependence on wave propagation angle. The wavenumber distributions of slow modes are mainly stretched perpendicular to the background magnetic field ( B 0) and weakly modulated by CD. In contrast, fast modes are subjected to a more significant CD modulation. Fast modes exhibit a weak, scale-independent anisotropy above the CD truncation scale. Below the CD truncation scale, the anisotropy of fast modes enhances as wavenumbers increase. As a result, fast-mode fractions in the total energy of compressible modes decrease with the increase of perpendicular wavenumber (to B 0) or wave propagation angle. Our findings reveal how the turbulence cascade is shaped by CD and its consequences for anisotropies in the space environment.

Jovian Electrons in the Inner Heliosphere: Opportunities for Multi-spacecraft Observations and Modeling

In this paper we explore the idea of using multi-spacecraft observations of Jovian electrons to measure the 3D distribution of these particles in the inner heliosphere. We present simulations of Jovian electron intensities along selected spacecraft trajectories for 2021 and compare these, admittedly qualitatively, to these measurements. Using the data–model comparison we emphasize how such a study can be used to constrain the transport parameters in the inner heliosphere, and how this can lead to additional insight into energetic particle transport. Model results are also shown along the expected trajectories of selected spacecraft, including the off-ecliptic phase of the Solar Orbiter mission from 2025 onward. Lastly, we revisit the use of historical data and discuss upcoming missions that may contribute to Jovian electron measurements.

The Effect of Solar Wind on Charged Particles’ Diffusion Coefficients

The transport of energetic charged particles through magnetized plasmas is ubiquitous in interplanetary space and astrophysics, and the important physical quantities are the parallel and perpendicular diffusion coefficients of energetic charged particles. In this paper, the influence of solar wind on particle transport is investigated. Using the focusing equation, we obtain parallel and perpendicular diffusion coefficients, accounting for the solar wind effect. For different conditions, the relative importance of the solar wind effect to diffusion is investigated. It is shown that, when energetic charged particles are close to the Sun, for parallel diffusion, the solar wind effect needs to be taken into account. These results are important for studying energetic charged particle transport processes in the vicinity of the Sun.

Effects of the resonance modification by electron cyclotron current drive on the linear and nonlinear dynamics of energetic particle driven magnetohydrodynamics modes in Heliotron J

This study investigates the effects of electron cyclotron current drive (ECCD) on the linear and nonlinear dynamics of energetic-particle (EP)-driven magnetohydrodynamics (MHD) instability (mode) in Heliotron J (HJ), using MEGA, a nonlinear hybrid EP-MHD simulation code. The effects of ECCD are included through the modification of the equilibrium magnetic field. Five MHD equilibria with plasma currents () of −2.00, −1.00, 0.00, 0.50, and 1.00 kA are considered, where negative and positive represent counter and co-ECCD, respectively. To account for the EP finite charge-exchange loss observed in HJ, a bump-on-tail velocity distribution is assumed for EP. The study finds that the stabilization trends in terms of the linear growth rate () are qualitatively similar to the FAR3D linear simulation and experimental results. In addition to the enhancements of the local magnetic shear and continuum damping by ECCD as proposed in previous studies, the study suggests that changes in the linear EP and shear Alfvén wave resonance are another important factor, especially for the mode. The study shows that in the currentless equilibrium, the high-velocity toroidicity-induced resonance between the mode and the co-passing EP () has the highest EP energy transfer rate, but is located slightly below the NBI injection velocity. The co-ECCD and cntr-ECCD can shift this resonance into higher and lower velocity regions, respectively. This can lead to resonance removal for co-ECCD, while cntr-ECCD can initially increase the number of resonant EPs for a small value of shift. If the shift of the resonance by cntr-ECCD is sufficiently large, the study observes stabilization due to the reduction in the EP population in the low-velocity region. The study also finds that changes in resonance affect nonlinear dynamics such as EP transport and frequency chirping. The convective transport of the co- and counter-passing EPs can be controlled, which can only be achieved for energetic-particle-mode due to the asymmetry in the resonance velocity between co- and cntr-passing EPs. The changes in the frequency chirping from the downward chirping branch to the upward chirping branch can also be explained by the changes in the resonance velocity caused by ECCD.

Simulations of Field Line Random Walk in Noisy Slab Turbulence

Magnetic field lines in magnetohydrodynamic turbulence are stochastic curves and, therefore, understanding their structure is a difficult problem of theoretical physics. One way to study the properties of such field lines is to perform simulations in which field line equations are solved numerically. In most previous articles studies have been performed in the context of the two-component turbulence model where one assumes a superposition of slab and two-dimensional modes. In the current paper we perform simulations for so-called noisy slab turbulence which is an example for three-dimensional turbulence with small Kubo numbers. We discuss in detail the challenges of the numerical magnetic field creation for this type of turbulence, perform the simulations, and compare our findings with analytical theory. It is demonstrated that a large number of wave modes is needed in order to obtain accurate results. Furthermore, we also show that analytical theories are very accurate for this type of turbulence model. The obtained results will also be useful for the analytical description of energetic particle transport in the considered turbulence model.

A Surrogate Model for Studying Solar Energetic Particle Transport and the Seed Population

AbstractThe high energy particles originating from the Sun, known as solar energetic particles (SEPs), contribute significantly to the space radiation environment, posing serious threats to astronauts and scientific instruments on board spacecraft. The mechanism that accelerates the SEPs to the observed energy ranges, their transport in the inner heliosphere, and the influence of suprathermal seed particle spectrum are open questions in heliophysics. Accurate predictions of the occurrences of SEP events well in advance are necessary to mitigate their adverse effects but prediction based on first principle models still remains a challenge. In this scenario, adopting a machine learning approach to SEP modeling and prediction is desirable. However, the lack of a balanced database of SEP events restrains this approach. We addressed this limitation by generating large data sets of synthetic SEP events sampled from the physics‐based model, Energetic Particle Radiation Environment Module (EPREM). Using this data, we developed neural networks‐based surrogate models to study the seed population parameter space. Our models, EPREM‐S, run thousands to millions of times faster (depending on computer hardware), making simulation‐based inference workflows practicable in SEP studies while providing predictive uncertainty estimates using a deep ensemble approach.

Investigation of alpha-particle transport by Alfvén eigenmodes in CFETR using a simplified diffusion model

The aim of this work is to analyze the alpha particle transport induced by Alfvén eigenmodes (AEs) based on the latest design of China Fusion Engineering Test Reactor (CFETR). Firstly, the stability of AEs is analyzed in order to understand the physical characteristics of AEs in CFETR with latest design parameters. AEs driven by alpha particles are analyzed using the gyrokinetic ion/fluid electron hybrid code GEM. The transport of alpha particles is investigated using a reduced energetic particle transport model based on the resonance broadened quasi-linear model. GEM results show that AEs with many toroidal mode numbers can be destabilized simultaneously by alpha particles in CFETR, and the toroidal mode number of the most unstable mode is . It is found that the excitation threshold of the alpha particle central beta for the mode is substantially below the expected alpha particle beta value of CFETR. All these results indicate that AEs can be driven strongly in CFETR plasmas. Furthermore, simulation results using the reduced transport model show that, due to multiple unstable AEs, the alpha particle density decreases in the core and alpha particles are transported from the core region () to the outer region (). In addition, we find that the minimum value of the safety factor and the central beta value of alpha particle have significant effects on the alpha particle transport, and the AE-induced alpha particle transport level becomes lower with smaller linear growth rates. Finally, it is found that the radial range of the redistributed alpha particle density profile depend on the radial location.

Energetic particle marginal stability profile for HL-2M integrated simulation based on neural network module

A critical gradient model is employed to develop a module of energetic particle (EP) marginal stability profiles in OMFIT integrated simulations for studying EP transport. Currently, each iteration of transport evolution is approximately 10 min in the integrated simulation, whereas, the EP marginal stability profile, which serves as an input in the integrated simulation could take much longer; the reason being a combination of the TGLFEP and EPtran codes is employed in our previous investigation. To reduce the simulation time, the critical gradient is predicted by a neural network instead of the TGLFEP code, and the EPtran code is revised with parallel computing, so that the running time of this module can be controlled to within 5 min. The predictions are in good agreement with previous approaches. The integrated simulation of HL-2M with Alfven eigenmodes transported by neutral beam EP profiles indicates that EP transport reduces the total pressure and current as expected, but could also under some conditions raise the safety factor in the core, which is favorable for reversed magnetic shear and high-performance plasmas.