Energetic Particle Physics Research Articles

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.

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.

Gamma ray diagnostics for high time resolution measurement in large helical device

A 1-inch LaBr3:Ce gamma ray scintillation detector, characterized by high time and energy resolution, was installed to advance energetic particle physics studies in Large Helical Device (LHD). We reduced the size of the scintillator from 3 inches to 1 inch to reinforce the radiation shielding to reduce the unwanted signal in the detector induced by fast neutrons and stray gamma rays according to the commissioning results. The radiation shielding composed of 10% borated polyethylene and lead was redesigned to suppress the gamma-ray induced signal based on the Monte Carlo three-dimensional neutron and gamma-ray transport code MCNP6. We increased the lead thickness from 50 mm to 77.5 mm to suppress the stray gamma ray effect. The gamma ray spectrum was measured in the hydrogen neutral beam heated deuterium plasma with 6LiF pellet injection. We might obtain a gamma ray peak near 0.48 MeV due to the 6Li(d,p'γ)7Li reaction.

Full f and gyrokinetic particle simulations of Alfvén waves and energetic particle physics

In this work, we focus on the development of the particle-in-cell scheme and its application to the studies of Alfvén waves and energetic particle (EP) physics in tokamak plasmas. The and full f schemes are formulated on the same footing adopting mixed variables and the pullback scheme for electromagnetic problems. The TRIMEG-GKX code (Lu et al 2021 J. Comput. Phys. 440 110384) has been upgraded using cubic spline finite elements and full f and schemes. The toroidal Alfvén eigenmode (TAE) driven by EPs has been simulated for the International Tokamak Physics activity (ITPA)-TAE case featured by a small electron skin depth , which is a challenging parameter regime for electromagnetic simulations, especially for the full f model. The simulation results using the scheme are in good agreement with previous work. Excellent performance of the mixed variable/pullback scheme has been observed for both full f and schemes. Simulations with mixed full f EPs and electrons and thermal ions demonstrate the good features of this novel scheme in mitigating the noise level. The full f scheme is a natural choice for EP physics studies, which allows for large variations of EP profiles and distributions in velocity space, providing a powerful tool for kinetic studies using realistic experimental distributions related to intermittent and transient plasma activities.

An E and B gyrokinetic simulation model for kinetic Alfvén waves in tokamak plasmas

The gyrokinetic particle simulation is a powerful tool for studies of transport, nonlinear phenomenon, and energetic particle physics in tokamak plasmas. While most gyrokinetic simulations make use of the scalar and vector potentials, a new model (GK-E&B) has been developed by using the E and B field in a general form and has been implemented in simulating kinetic Alfvén waves in uniform plasma [Chen et al., Sci. China: Phys., Mech. Astron. 64, 245211 (2021)]. In our work, the Chen et al. GK-E&B model has been expressed, in general, tokamak geometry using the local orthogonal coordinates and general tokamak coordinates. Its reduction for uniform plasma is verified, and the numerical results show good agreement with the original work. The theoretical dispersion relation and numerical results in the local model in screw pinch geometry are also in excellent agreement. Numerical results show excellent performance in a realistic parameter regime of burning plasmas with high values of β/(Mek⊥2ρi2), which is a challenge for traditional methods due to the “cancellation” problem. As one application, the GK-E&B model is implemented with kinetic electrons in the local single flux surface limit. With the matched International Tokamak Physics Activity-Toroidicity-induced Alfvén Eigenmodes parameters adopted, numerical results show the capability of the GK-E&B in treating the parallel electron Landau damping for realistic tokamak plasma parameters. As another application, the global GK-E&B model has been implemented with the dominant electron contribution in the cold electron limit. Its capability in simulating the finite E|| due to the finite electron mass is demonstrated.

The Evolution of Research on Abundances of Solar Energetic Particles

Sixty years of study of energetic particle abundances have made a major contribution to our understanding of the physics of solar energetic particles (SEPs) or solar cosmic rays. An early surprise was the observation in small SEP events of huge enhancements in the isotope 3He from resonant wave–particle interactions, and the subsequent observation of accompanying enhancements of heavy ions, later found to increase 1000-fold as a steep power of the mass-to-charge ratio A/Q, right across the elements from H to Pb. These “impulsive” SEP events have been related to magnetic reconnection on open field lines in solar jets; similar processes occur on closed loops in flares, but those SEPs are trapped and dissipate their energy in heat and light. After early controversy, it was established that particles in the large “gradual” SEP events are accelerated at shock waves driven by wide, fast coronal mass ejections (CMEs) that expand broadly. On average, gradual SEP events give us a measure of element abundances in the solar corona, which differ from those in the photosphere as a classic function of the first ionization potential (FIP) of the elements, distinguishing ions and neutrals. Departures from the average in gradual SEPs are also power laws in A/Q, and fits of this dependence can determine Q values and thus estimate source plasma temperatures. Complications arise when shock waves reaccelerate residual ions from the impulsive events, but excess protons and the extent of abundance variations help to resolve these processes. Yet, specific questions about SEP abundances remain.

The development of an implicit full f method for electromagnetic particle simulations of Alfvén waves and energetic particle physics

In this work, an implicit scheme for particle-in-cell/Fourier electromagnetic simulations is developed and applied to studies of Alfvén waves in one dimension and three dimensional tokamak plasmas. An analytical treatment is introduced to achieve efficient convergence of the iterative solution of the implicit field-particle system. First, its application to the one-dimensional uniform plasma demonstrates its applicability in a broad range of β/me values. Second, toroidicity induced Alfvén eigenmodes (TAE) are simulated in a three dimensional axisymmetric tokamak plasma, using the widely studied case defined by the International Tokamak Physics Activity (ITPA) Energetic Particle (EP) Topical Group. The real frequency and the growth (or damping) rate of the TAE with (or without) EPs agree with previous results reasonably well. The full f electromagnetic particle scheme established in this work provides a possible natural choice for EP transport studies where large profile variation and arbitrary particle distribution functions need to be treated in kinetic simulations.

MHD-kinetic hybrid code based on structure-preserving finite elements with particles-in-cell

We present a STRUcture-Preserving HYbrid code - STRUPHY - for the simulation of magneto-hydrodynamic (MHD) waves interacting with a population of energetic particles far from thermal equilibrium (kinetic species). The implemented model features linear, ideal MHD equations in curved, three-dimensional space, coupled nonlinearly to the full-orbit Vlasov equations via a current coupling scheme. The algorithm is based on finite element exterior calculus (FEEC) for MHD and particle-in-cell (PIC) methods for the kinetic part; it provably conserves mass, energy, and the divergence-free magnetic field, irrespective of metric (= space curvature), mesh parameters and chosen order of the scheme. These properties enable reliable long-time simulations of energetic particle physics in complex geometries, covering the whole range of MHD waves. In STRUPHY, the finite element spaces are built from tensor products of univariate B-splines on the logical cuboid and can be made high-order by increasing the polynomial degree. Time-stepping is based on a skew-symmetric splitting with implicit sub-steps, mitigating CFL conditions from fast magneto-acoustic waves. High-order time splitting schemes can be used in this regard.

Progress of Experimental Studies in the HL-2A Tokamak

During the last several years, the HL-2A experiment has made significant progress in the following areas: (1) lower-hybrid wave (LHW) heating and current drive, (2) plasma confinement and turbulent transport, (3) magnetohydrodynamic (MHD) instabilities and energetic particle physics and (4) H-mode and edge localized mode (ELM) control. The results show that the LHW system working in the co-current mode can reach higher driving efficiency and full non-inductive lower-hybrid current drive (LHCD) has been achieved. The intrinsic poloidal torque characterized by the divergence of the residual stress is deduced from synthesis for the first time. The dynamics of spectral symmetry breaking in drift wave turbulence is in good agreement with the development of the poloidal torque to drive the edge poloidal flow. The influence of the cross-phase dynamics on turbulent stress was also investigated. The ion internal transport barrier has been observed in the NBI-heated plasma, and inside the barrier the ion thermal transport is reduced to the neoclassical level. Besides, micro-turbulence is modulated by the rotation frequency of the magnetic island, and this modulation effect is related to a critical island width. Strong E × B shear is found at the island boundary. Three kinds of axisymmetric modes, beta-induced Alfven eigenmode (BAE), toroidal Alfven eigenmode (TAE) and the ellipticity-induced Alfven eigenmode (EAE), are found to be driven unstable by nonlinear mode coupling between Alfven eigenmodes and tearing mode which is well explained by the nonlinear gyrokinetic theory. The fishbone and tearing modes were actively controlled by the electron cyclotron resonance heating (ECRH). The dynamics of the edge plasma flows and turbulence during the L–I–H transition have been dedicatedly investigated. The geodesic acoustic mode (GAM) and limit cycle oscillation (LCO) coexist for a short time and disappear in the H-mode plasma with the increasing of E × B shear flow before the I–H transition, which plays an important role in the turbulence suppression. Different techniques, such as LHW, ECRH, resonant magnetic perturbation (RMP), and impurity seeding by the laser blow-off (LBO) and supersonic molecular beam injection (SMBI), have been successfully applied to control the large ELMs. It has been found that pedestal turbulence enhancement might be responsible for the observed mitigation effect.

Fast-ion D alpha diagnostic with 3D-supporting FIDASIM in the Large Helical Device

For understanding the physics of energetic particles, the deuterium experimental campaigns started in Large Helical Device (LHD) from March 2017. To investigate the behavior of energetic particles, a Fast-ion D Alpha (FIDA) diagnostic was installed on the LHD. In the FIDA diagnostic, the Doppler-shifted D alpha light from fast-neutrals are utilized as signals of energetic particles, where these fast-neutrals are produced by the charge exchange process between fast-ions in the plasma and actively injected neutrals by neutral beam (NB). The advantages of the FIDA diagnostic are the velocity and the spatially resolved measurement of fast-ions at the crossing point between its line of sight (LOS) and the incident line of NB. The most recent FIDASIM is enhanced to simulate signals produced in three-dimensional magnetic configurations. The new version of FIDASIM uses the fast-ion distribution function produced by GNET as input to simulate FIDA signals at LHD. In order to validate the new version of the code, measurements of radial profiles of fast-ions using the FIDA diagnostic are performed in magnetohydrodynamic (MHD)-quiescent plasmas. The measured spectra are in good agreement with the theoretical prediction by 3D-supporting FIDASIM at the center of the plasma (R= 3.5 m∼3.7 m, reff/a99= −0.28∼0.05) on the LHD when the line averaged electron density is ne_avg < 1.23 × 1019 m−3. On the other hand, the measured spectra are in disagreement with the theoretical prediction by 3D-supporting FIDASIM at even the center of the plasma when the line averaged electron density is ne_avg ≥ 1.23 × 1019 m−3.

Summary of the 16th IAEA technical meeting on energetic particles in magnetic confinement systems—theory of plasma instabilities

This report summarizes the contributions presented at the 16th IAEA Technical Meeting on Energetic Particles in Magnetic Confinement Systems—Theory of Plasma Instabilities, held in Shizuoka, Japan, 3–6 September 2019. The meeting brought together about 100 experts from nuclear fusion research sites worldwide to discuss the physics of energetic particles and plasma instabilities, at the first joint meeting from the two scientific disciplines. The main topics of the meeting were: alpha particles physics; transport of energetic particles; effects of energetic particles in magnetic confinement fusion devices; collective phenomena such as Alfvén eigenmodes, energetic particle modes and others; runaway electrons and disruptions; diagnostics for energetic particles; control of energetic particles confinement; multiscale physics and instabilities in burning plasmas.

Small, Low-energy, Dispersive Solar Energetic Particle Events Observed by Parker Solar Probe

The Energetic Particle Instrument–Low Energy (EPI-Lo) experiment has detected several weak, low-energy (~30–300 keV nucleon⁻¹) solar energetic particle (SEP) events during its first two closest approaches to the Sun, providing a unique opportunity to explore the sources of low-energy particle acceleration. As part of the Parker Solar Probe (PSP) Integrated Science Investigation of the Sun (IS⊙IS) suite, EPI-Lo was designed to investigate the physics of energetic particles; however, in the special lowest-energy time-of-flight only product used in this study, it also responds to solar photons in a subset of approximately sunward-looking apertures lacking special light-attenuating foils. During the first three perihelia, in a frame rotating with the Sun, PSP undergoes retrograde motion, covering a 17° heliographic longitudinal range three times during the course of the ~11-day perihelion passes, permitting a unique spatial and temporal study into the location, correlation, and persistence of previously unmeasurable SEPs. We examine the signatures of these SEPs (during the first PSP perihelion pass only) and the connection to possible solar sources using remote observations from the Solar Dynamics Observatory (SDO), the Solar TErrestrial RElations Observatory (STEREO), and the ground-based Global Oscillation Network Group (GONG). The orientation of the Sun relative to STEREO, SDO, and GONG makes such identifications challenging, but we do have several candidates, including an equatorial coronal hole at a Carrington longitude of ~335°. To analyze observations from EPI-Lo, which is a new type of particle instrument, we examine instrumental effects and provide a preliminary separation of the ion signal from the photon background.

Preface

G P ZankEditor, Proceedings of the 18th Annual International Astrophysics Conference,Center for Space Plasma and Aeronomic Research (CSPAR) and Department of SpacePhysics, University of Alabama in Huntsville, Huntsville AL 35805, USAE-mail: garyp.zank@gmail.comThe 18th Annual International Astrophysics Conference was held at the Sheraton Pasadena, Pasadena, California, USA from February 18 to February 22, 2019. The theme of the meeting was "The Physics of Energetic Particles: Universal Processes from the Solar Corona to the Very Local Interstellar Medium and the Physics they Enable.” Energetic particles are ubiquitous throughout the solar wind, being integral to the solar corona, accelerated by interplanetary shock waves, present in the quiet solar wind, possibly accelerated by magnetic flux ropes or magnetic islands in the vicinity of the heliospheric current sheet, at planetary bow shocks, present in the distant heliosphere and the heliospheric termination shock, in the inner heliosheath, and now observed by Voyager 1 in the very local interstellar medium. Although numerous mechanisms have been identified that accelerate particles, no consensus exists as to which is primary, and indeed several acceleration processes can operate simultaneously in close proximity. The acceleration of particles to high energies represents a fundamental dissipative process for the plasma and can therefore modify the underlying plasma physical processes in important ways. The meeting explored observations, theories, and looked forward to anticipated new missions that will shine a light on this theme, the Parker Solar Probe and IMAP. The 18th AIAC explored these topics via a combination of invited 25 minute talks and a few 40 minute plenary talks.The 18th AIAC also celebrated Prof. J.R. (Randy) Jokipii’s 80th birthday. Randy has been a centralfigure in the physics of energetic particles and transport, including particle acceleration, particularly at shock waves, and the modulation of cosmic rays in the heliosphere. His seminal early paper, “Cosmicray propagation. I. Charged particles in a random magnetic field” defined the field for decades after its publication in 1966. This meeting was an opportunity to celebrate Randy’s accomplishments with some selected plenary presentations.Finally, there are several people to thank for their kindness and considerable help with the 18th AIAC. Adele Corona suffered a horrific car accident on the second day of the meeting. She has thankfully largely recovered. However, we are profoundly grateful to Prof.’s Dave McComas and Nathan Schwadron for selflessly stepping forward to taking over the running and management of the meeting after Adele was rushed to hospital. Their kindness and organizational skills ensured that the meeting continued with virtually no interruption. I should like to thank them both and of course Adele Corona for another outstanding meeting, and her help in providing the logistical support for this volume of papers.

Overview of HL-2A recent experiments

The mission of HL-2A is to explore the key physical topics relevant to ITER and to the advanced tokamak operation (e.g. the operation of future HL-2M), such as the access of H-mode, energetic particle physics, edge-localized mode (ELM) mitigation/suppression and disruption mitigation. Since the 2016 Fusion Energy Conference, the HL-2A team has focused on the investigations on the following areas: (i) pedestal dynamics and L–H transition, (ii) techniques of ELM control, (iii) turbulence and transport, (iv) energetic particle physics. The HL-2A results demonstrated that the increase of mean shear flow plays a key role in triggering L–I and I–H transitions, while the change of flow is mainly induced by the ion pressure gradient. Both mitigation and suppression of ELMs were realized by laser blow-off seeded impurity (Al, Fe, W). The 30% Ne mixture supersonic molecular beam injection seeding also robustly induced ELM mitigation. The ELMs were mitigated by low-hybrid current drive. The stabilization of m/n = 1/1 ion fishbone activities by electron cyclotron resonance heating was found on the HL-2A. A new m/n = 2/1 ion fishbone activity was observed recently, and the modelling indicated that passing fast ions dominantly contribute to the driving of 2/1 fishbone. The non-linear coupling between the toroidal Alfven eigenmode and tearing mode (TM) led to the generation of a high frequency mode with the toroidal mode number n = 0. The turbulence was modulated by TM when the island width exceeds a threshold and the modulation is localized merely in the inner area of the islands. Meanwhile, turbulence radial spread took place across the island region.

Summary of the fundamental plasma physics session in the first AAPPS-DPP conference

Understanding of fundamental plasma physics contributes to progress in plasma science in many fields, such as magnetic confinement fusion, solar/space physics and application of plasma techniques. Intersecting common fundamental problems stimulate cross-fertilization of different topical areas. This paper summarizes the highlights of progress reported in the Fundamental Plasma Physics Session of the first conference of the Association of Asia Pacific Physical Societies-Division of plasma physics. Progress in the following fields is reviewed: (1) basic plasma theory; (2) self-organization; (3) reconnection and particle acceleration; (4) magneto-hydrodynamic instabilities and energetic particle physics; (5) magnetic confinement physics fundamentals; and (6) exploiting fundamentals for performance.

Velocity-space sensitivity of time-of-flight neutron spectrometer at EAST in deuterium plasma

The Time-Of-Flight Enhanced Diagnostics (TOFED) neutron spectrometer with a double-ring structure has been installed at the Experimental Advanced Superconducting Tokamak (EAST) to perform advanced neutron emission spectroscopy diagnosis for deuterium plasma. In order to reduce the random coincidence from the background neutrons and gamma-rays, TOFED was moved outside the experimental hall and placed in the newly-built nuclear diagnostics laboratory in 2017. In this paper, the instrument-specific weight functions of TOFED are derived by taking the instrument response matrix and the radial line of sight in this new layout into consideration. The results show that the instrument is predominantly sensitive to counter-passing particles in the region where time-of-flights < 69.4 ns, while events at higher time-of-flights (corresponding lower neutron energies) are mostly representative of co-passing ions. The instrument-specific weight functions express the relationship between data in a given channel of the spectrum and the velocity space region that contributes to that. The results can be applied for energetic particle physics studies at EAST, in particular to compare data from different diagnostic techniques.

Progress of microwave diagnostics development on the HL-2A tokamak

Sets of powerful microwave diagnostic systems have been developed on the HL-2A tokamak, including 18 subsystems with high spatial and temporal resolutions. These systems have been applied in the HL-2A tokamak experiments to investigate the core plasma transport, turbulence, MHD, energetic particle physics and so on. In this paper, the microwave diagnostics and their applications are reviewed. Some new technologies, including GHz sampling digital correlation electron cyclotron emission (DCECE), multi-channel correlation reflectometers, and solid state terahertz interferometers, are also presented in the paper.

Review of the experiments for energetic particle physics on HL-2A

This paper reviews the energetic particle (EP) experiments during electron cyclotron resonance heating (ECRH) and neutral beam injection in the HL-2A tokamak. A number of important results are summarized, which relate to ITER physics, including the behavior of the multi-mode instability, the nonlinear interaction between wave–wave and wave–particles, the losses of EP induced by the instabilities, the effect of the EP instabilities on the thermal plasma confinement and the control of the EP instabilities by means of ECRH. Systematic experiments indicate that when the drive is great enough, the nonlinear effects and the multi-mode coexistence may play an important role, which affect the transport both of the EPs and the background plasma confinement, and these instabilities could be controlled. Some new phenomena about the EP induced instabilities discovered recently on the device, such as high frequency reversed shear Alfvén eigenmodes, Alfvénic ion temperature gradient modes, the geodesic acoustic mode induced by energetic electrons excited by interaction between tearing mode and beta induced Alfvén eigenmode and double e-fishbone in negative magnetic shear discharges etc, have also been presented in the paper.

The HESPERIA HORIZON 2020 Project and Book on Solar Particle Radiation Storms Forecasting and Analysis

AbstractThis article presents the High Energy Solar Particle Events forecasting and Analysis (HESPERIA) project, supported by the HORIZON 2020 programme of the European Union (Project 637324) as well as the resultant recently published book entitled Solar Particle Radiation Storms Forecasting and Analysis, The HESPERIA HORIZON 2020 Project and Beyond, edited by Malandraki and Crosby, Springer, Astrophysics and Space Sciences Library, 2018, ISBN 978‐3‐319‐60051‐2. The book reviews the results of the HESPERIA project as well as our current understanding of solar energetic particle physics.

Study of energetic particle physics with advanced ECEI system on the HL-2A tokamak

Understanding the physics of energetic particles (EP) is crucial for the burning plasmas in next generation fusion devices such as ITER. In this work, three types of internal kink modes (a saturated internal kink mode (SK), a resonant internal kink mode (RK), and a double e-fishbone) excited by energetic particles in the low density discharges during ECRH/ECCD heating have been studied by the newly developed 24(poloidal) × 16(radial) = 384 channel ECEI system on the HL-2A tokamak. The SK and RK rotate in the electron diamagnetic direction poloidally and are destabilized by the energetic trapped electrons. The SK is destabilized in the case of q min > 1, while the RK is destabilized in the case of q min = 1/1 modes propagating in the opposite directions poloidally, has been observed during plasma current ramp-up with counter-ECCD. Strong thermal transfer and mode coupling between the two m/n = 1/1 modes have been studied.