Effects Of Energetic Particles 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.

Simulation of non-resonant high-order harmonics energetic particle modes in tokamak plasmas

Based on the parameters of the HL-2A experiment, the effect of energetic particles (EPs) on non-resonant high-order harmonics energetic particle modes (EPMs) with q min>1 is investigated in the present work. Hybrid kinetic-magnetohydrodynamic nonlinear code M3D-K is performed to simulate the linear properties and the nonlinear evolution of the non-resonant EPM during neutral beam injection (NBI). To deeply understand the physical mechanism of interaction resonant between energetic-ions and non-resonant EPM, this work compares the effects of passing energetic particles and trapped energetic particles on the non-resonant EPM instabilities. It is numerically identified that EPs’ effects on high n harmonics (m/n = 2/2, 3/3, 4/4) instability are more obvious than the m/n = 1/1 mode. Furthermore, the effects of energetic particles injection energy, the minimum safety factor q min , toroidal rotation and beam ion distribution on the features of high n harmonics are also investigated specifically. Toroidal rotation is found to suppress high n harmonics, which is more obvious for the modes driven by trapped particles. Nonlinear simulation results show that these non-resonant high n harmonics can induce larger energetic ion transport, which may affect the plasma confinement performance.

Predictive modeling of Alfvén eigenmode stability in inductive scenarios in JT-60SA

The JT-60SA device offers unique conditions before ITER for the study of the interaction of energetic particles with plasma waves. With similar dimensions to JET, e.g., a major radius but with a slightly more elongated plasma volume, JT-60SA is used as a high-power device where additional heating power (including 10 MW of the 500 keV Neutral Beam Injection) of up to 41 MW and the potential for high non-inductive plasma current operation pave the path for numerous challenges in physics on MHD stability, in particular, when considering the effects of energetic particles. Several operational scenarios with ITER and DEMO-relevant plasma regimes, in terms of non-dimensional plasma parameters, are anticipated. In this work, the stability of Alfvén eigenmodes (AEs) in variants of two of the most relevant operational scenarios with single null is analyzed: a full Ip inductive scenario at high density (1.1 × 1020 m−3 on-axis electron density) and 5.48MA/2.05T toroidal plasma current and magnetic field, and an advanced (hybrid) scenario with an ion energy transport barrier (ITB) and 3.5MA/2.28T toroidal plasma current and magnetic field. The workflow included the CRONOS code to establish the scenario, the ASCOT code to calculate the slowing-down energetic particle distributions for a positive/negative ion source-based neutral beam, and the MISHKA/CASTOR-K suite to calculate the MHD spectra of AEs and the associated drive/damping contributions from the NBI energetic ions, as well as the thermal ion landau damping. The systematic analysis, over a large Fourier space of the toroidal mode number/mode frequency, provides evidence that although a significant fraction of supra-Alfvénic particles stemming from the negative ion source-based neutral beam (500 keV) can, in some cases, drive to AEs in both scenarios, it is not enough to overcome the thermal ion landau damping. In addition, the advanced scenario with ITB is shown to be stable against AEs localized in the vicinity of the barrier as well, offering good prospects of sustainability of the plasma performance and of ITB. Finally, some sensitivity scan results are shown on the influence of fast ion density and q-profile on the AE mode spectra and stability.

How accurate are flux-tube (local) gyrokinetic codes in modeling energetic particle effects on core turbulence?

Flux-tube (local) gyrokinetic codes are widely used to simulate drift-wave turbulence in magnetic confinement devices. While a large number of studies show that flux-tube codes provide an excellent approximation for turbulent transport in medium-large devices, it still needs to be determined whether they are sufficient for modeling supra-thermal particle effects on core turbulence. This is called into question given the large temperature of energetic particles (EPs), which makes them hardly confined on a single flux-surface, but also due to the radially broad mode structure of EP-driven modes. The primary focus of this manuscript is to assess the range of validity of flux-tube codes in modeling fast ion effects by comparing radially global turbulence simulations with flux-tube results at different radial locations for realistic JET parameters using the gyrokinetic code GENE. To extend our study to a broad range of different plasma scenarios, this comparison is made for four different plasma regimes, which differ only by the profile of the ratio between the plasma kinetic and magnetic pressure. The latter is artificially rescaled to address the (i) electrostatic limit and regimes with (ii) marginally stable, (iii) weakly unstable and (iv) strongly unstable fast ion modes. These EP-driven modes are identified as Alfvénic ion temperature gradient modes (AITG)/kinetic beta-induced Alfvén eigenmodes (KBAE) via linear ORB5 and LIGKA simulations. It is found that the local flux-tube simulations can recover well the global results only in the electrostatic and marginally stable cases. When the AITG/KBAE becomes linearly unstable, the local approximation fails to correctly model the radially broad fast ion mode structure and the consequent global zonal patterns. According to this study, global turbulence simulations are likely required in regimes with linearly unstable AITG/KBAEs. In conditions with different fast ion-driven modes, these results might change.

Recent progress in the interaction between energetic particles and tearing modes.

The dynamics of energetic particles and tearing modes and the interactions between them are of great significance for magnetically confined fusion plasmas. In this review, we focus on these issues in the context of tokamak plasmas. The interaction between energetic particles and tearing modes is considered from two perspectives: (i) the influence of energetic particles on tearing modes and (ii) the transport of energetic particles by tearing modes. The influence of energetic particles on tearing modes is described on the basis of a general dispersion relation for tearing modes. The effects of energetic particles are considered separately in the outer region and the island region of a tearing mode. The physics mainly results from the modification of the perturbed parallel current by energetic particles without wave-particle resonance. In addition, the resonance between energetic particles and tearing modes is also reviewed. For the transport of energetic particles, transport of both circulating and trapped energetic particles by tearing mode is reviewed. Our descriptions of physical phenomena here are based on an analytical approach, while the experiments and simulations are used to illustrate and confirm our results. Finally, a number of open issues are discussed.

Effects of anisotropic energetic particles on zonal flow residual level

In tokamak plasmas, the interaction among the micro-turbulence, zonal flows (ZFs) and energetic particles (EPs) can affect the turbulence saturation level and the consequent confinement quality and thus, is important for future burning plasmas. In this work, the EP anisotropy effects on the ZF residual level are studied by using anisotropic EP distributions with dependence on pitch. Significant effects on the long wavelength ZFs have been found when small to moderate width around the dominant pitch in the EP distribution function is assumed. In addition, it is found that ZF residual level is enhanced by barely passing/trapped and/or deeply trapped EPs, but it is suppressed by well passing and/or intermediate trapped EPs. Numerical calculation shows that for ASDEX Upgrade plasmas, typical EP distribution functions can bring in −3% to +5.5% mitigation/enhancement in ZF residual level, depending on the EP distribution functions.

Effects of energetic particles on the density-gradient-driven collisionless trapped electron mode instability in tokamak plasmas

We investigate analytically the effects of energetic particles (EPs) on the instability of the density-gradient-driven collisionless trapped electron mode (CTEM) through linear gyrokinetic theory and bounce kinetic theory in tokamak plasmas. The effects of EPs, including fusion-born alpha particles and neutral-beam-injection-driven beam ions, on the CTEM instability are compared for the dynamic model with slowing-down (SD) and equivalent Maxwellian (EM) equilibrium EP distribution functions and dilution model. It is found that the density-gradient-driven CTEM instability in the long wavelength regime can be further destabilized by EPs mainly due to the downshift in the real frequency of the mode by dilution effects. This is attributed to more resonant electrons around the smaller phase velocity of the drift wave and the consequent stronger excitation of CTEM instability. The growth rate is slightly higher for the dilution model as compared to that for the dynamic model since the Landau damping effects of EPs are neglected in the dilution model. Moreover, there is no significant difference in the growth rate between the cases of SD and EM equilibrium EP distribution functions, except for the case of alpha particles and with relatively higher electron temperature.

Influence of Off-Axis Neutral Beam Injection on Resistive Wall Mode Stability

Investigated are the effects of energetic particles (EPs), arising from off-axis neutral beam injection, on the resistive wall mode (RWM) instability in tokamaks. The analytic results for a cylindrical plasma show that EPs contribute the strongest kinetic damping to the RWM, as the deposition radius approaches the peak location of the mode eigenfunction. Toroidal computations are also carried out utilizing the MARS-K code, which includes EPs with anisotropic equilibrium distribution in the particle pitch angle space. These self-consistent numerical results confirm the analytic finding, that the RWM growth rate decreases as the peak of the EPs’ radial density profile approaches the maximum of the mode displacement. This study thus demonstrates the possibility of the RWM stabilization by the off-axis neutral beam injection in tokamak experiments.

MOVES – IV. Modelling the influence of stellar XUV-flux, cosmic rays, and stellar energetic particles on the atmospheric composition of the hot Jupiter HD 189733b

ABSTRACT Hot Jupiters provide valuable natural laboratories for studying potential contributions of high-energy radiation to pre-biotic synthesis in the atmospheres of exoplanets. In this fourth paper of the Multiwavelength Observations of an eVaporating Exoplanet and its Star (MOVES) programme, we study the effect of different types of high-energy radiation on the production of organic and pre-biotic molecules in the atmosphere of the hot Jupiter HD 189733b. Our model combines X-ray and UV observations from the MOVES programme and 3D climate simulations from the 3D Met Office Unified Model to simulate the atmospheric composition and kinetic chemistry with the STAND2019 network. Also, the effects of galactic cosmic rays and stellar energetic particles are included. We find that the differences in the radiation field between the irradiated dayside and the shadowed nightside lead to stronger changes in the chemical abundances than the variability of the host star’s XUV emission. We identify ammonium (NH4+) and oxonium (H3O+) as fingerprint ions for the ionization of the atmosphere by both galactic cosmic rays and stellar particles. All considered types of high-energy radiation have an enhancing effect on the abundance of key organic molecules such as hydrogen cyanide (HCN), formaldehyde (CH2O), and ethylene (C2H4). The latter two are intermediates in the production pathway of the amino acid glycine (C2H5NO2) and abundant enough to be potentially detectable by JWST.

Analytic investigation of combined effects of anisotropic thermal transport and energetic particles on stability of resistive plasma resistive wall mode

The combined effects of anisotropic thermal transport and trapped energetic particles (EPs) on the stability of the resistive plasma resistive wall mode (RPRWM) are investigated by an energy-principle based analytical model. The results qualitatively confirm that of a recent toroidal modeling study [Bai et al., Phys. Plasmas 27, 072502 (2020)], in which the thermal transport can stabilize the RPPWM depending on the parameters of both the plasma equilibrium and energetic particles. The analytical model predicts a complete stabilization of the RPRWM in highly resistive plasmas, at sufficiently high EPs' pressure, and finite plasma flow. The stabilizing effect of thermal transport originates from its enhancement of energy dissipations associated with both the resistive layer and the trapped energetic particles.

Progress in extending high poloidal beta scenarios on DIII-D towards a steady-state fusion reactor and impact of energetic particles

To prepare for steady-state operation of future fusion reactors (e.g. the International Thermonuclear Experimental Reactor and China Fusion Engineering Test Reactor (CFETR)), experiments on DIII-D have extended the high poloidal beta (βP) scenario to reactor-relevant edge safety factor q95 ∼ 6.0, while maintaining a large-radius internal transport barrier (ITB) using negative magnetic shear. Excellent energy confinement quality (H98y2 > 1.5) is sustained at high normalized beta (βN ∼ 3.5). This high-performance ITB state with Greenwald density fraction near 100% and qmin ≥ 3 is achieved with toroidal plasma rotation Vtor ∼ 0 at ρ ≥ 0.6. This is a key result for reactors expected to have low Vtor. At high βP (≥1.9), large Shafranov shift can stabilize turbulence leading to a high confinement state with a low pedestal and an ITB. At lower βP (<1.9), negative magnetic shear in the plasma core contributes to turbulence suppression and can compensate for reduced Shafranov shift to continue to access a large-radius ITB and excellent confinement with low Vtor, consistent with the results of gyrofluid transport simulations. These high-βP cases are characterized by weak/no Alfvén eigenmodes (a.e.) and classical fast-ion transport. At high density, the fast-ion deceleration time decreases and Δβfast is lower; these reduce a.e. drive. The reverse-shear Alfvén eigenmodes are weaker or stable because the negative magnetic shear region is located at higher radius, away from the peaked fast-ion profile. Resistive wall modes can be a limitation at simultaneous high βN, low internal inductance, and low rotation. Analysis suggests that additional off-axis external current drive could provide a more stable path at reduced q95. Based on a DIII-D high-βP plasma with large-radius ITB, two scenarios are proposed for CFETR Q = 5 steady-state operation with ∼1 GW fusion power: a lower- ( ∼ 0.66) and a higher- ( ∼ 0.75) case. Using a Landau closure model, multiple energetic particle (EP) effects on the a.e. stability are analyzed modifying the growth rate of the a.e.s triggered by the neutral-beam-injection EPs and alpha particles, although the stabilizing/destabilizing effect is weak for the cases analyzed. The stabilizing effects of the combined EP species β, energy, and density profile in CFETR need further investigation.

Turbulence Suppression by Energetic Particle Effects in Modern Optimized Stellarators.

Turbulent transport is known to limit the plasma confinement of present-day optimized stellarators. To address this issue, a novel method to strongly suppress turbulence in such devices is proposed, namely the resonant wave-particle interaction of suprathermal particles-e.g., from ion-cyclotron-resonance-frequency heating-with turbulence-driving microinstabilities like ion-temperature-gradient modes. The effectiveness of this mechanism is demonstrated via large-scale gyrokinetic simulations, revealing an overall turbulence reduction by up to 65% in the case under consideration. Comparisons with a tokamak configuration highlight the critical role played by the magnetic geometry and the first steps into the optimization of fast particle effects in stellarator devices are discussed. These results hold the promise of new and still unexplored stellarator scenarios with reduced turbulent transport, essential for achieving burning plasmas in future devices.

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.

Toroidal modeling of anisotropic thermal transport and energetic particle effects on stability of resistive plasma resistive wall mode

Effects of anisotropic thermal transport on the linear stability of the resistive plasma resistive wall mode (RPRWM) are investigated by the magnetohydrodynamic-kinetic hybrid code MARS-K [Liu et al., Phys. Plasmas 15, 112503 (2008)], including the kinetic contribution from energetic particles (EPs). It is found that thermal transport can further stabilize the RPRWM in the presence of drift kinetic contributions from EPs. This is different from the fluid model, which always predicts destabilization of the mode by thermal transport. Furthermore, the thermal transport effect is found to amplify the role played by an adiabatic term, associated with the radial distribution of EPs' birth energy, in modifying the mode stability as well as the mode eigenfunction. The shape of the equilibrium profiles of EPs, in particular that of the temperature, is also found to strongly affect the mode stabilization. This profile effect is more pronounced in the absence of thermal transport. MARS-K computations show that the stabilizing effect by thermal transport is more likely to occur at slower plasma rotation and lower EP energy.

Toroidal models of magnetic field with twisted structure

We present and discuss properties of the following magnetic field models in a magnetic cloud: Miller and Turner solution, modified Miller–Turner solution, Romashets–Vandas toroidal and integral models, and Krittinatham–Ruffolo model. Helicity of the magnetic field in all the models is the main feature of magnetic clouds. The first three models describe the magnetic field inside an ideal torus. In the integral model, parameters of a generating torus ambiguously determine the volume and form of the magnetic field region. In the Krittinatham–Ruffolo model, the cross-section radius of the torus is variable, thereby it corresponds more closely to the real form of magnetic clouds in the inner heliosphere. These models can be used to interpret in-situ observations of the magnetic flux rope, to study a Forbush decrease in magnetic clouds and transport effects of solar energetic particles injected into a coronal mass ejection.

Тороидальные модели магнитного поля с винтовой структурой

We present and discuss properties of the following magnetic field models in a magnetic cloud: Miller and Turner solution, modified Miller–Turner solution, Romashets–Vandas toroidal and integral models, and Krittinatham–Ruffolo model. Helicity of the magnetic field in all the models is the main feature of magnetic clouds. The first three models describe the magnetic field inside an ideal torus. In the integral model, parameters of a generating torus ambiguously determine the volume and form of the magnetic field region. In the Krittinatham–Ruffolo model, the cross-section radius of the torus is variable, thereby it corresponds more closely to the real form of magnetic clouds in the inner heliosphere. These models can be used to interpret in-situ observations of the magnetic flux rope, to study a Forbush decrease in magnetic clouds and transport effects of solar energetic particles injected into a coronal mass ejection.

Effect of energetic particles on pulsed magnetron sputtering of hard nanocrystalline MBCN (M = Ti, Zr, Hf) films with high electrical conductivity

Nanocrystalline MBCN (M = Ti, Zr, Hf) films were deposited onto Si substrates using pulsed magnetron co-sputtering of a single B4C-M target (at a fixed 45% fraction of M in the target erosion area) with a repetition frequency of 10 kHz, voltage pulse duration of 85 μs and a fixed target power. The substrate temperature was adjusted to 450 °C during the depositions on the substrates at a floating potential. The total pressure of 95% Ar + 5% N2 gas mixture was 0.5 Pa for TiBCN and ZrBCN, but it had to be increased up to 1.7 Pa for HfBCN to decrease very high compressive stress in the films resulting in their delamination from the substrate. Energy-resolved mass spectroscopy was used to correlate the energy of Ar+ ions bombarding the growing films with high positive voltage overshoots after the negative voltage pulses. Monte-Carlo simulations were carried out to estimate the energy and flux of sputtered M atoms and backscattered Ar atoms at the substrate. It was found that the energy and flux of the backscattered Ar atoms increase significantly with the mass of the M atoms in the target. The MBCN films with 41–46 at.% of M, 25–31 at.% of B, 7–10 at.% of C and 14–22 at.% of N exhibit a high hardness (18–37 GPa), high elastic recovery (69–85%), high H/E⁎ ratio (0.10–0.14) and low electrical resistivity (1.7–2.7 μΩm) at a low internal stress (less than 0.8 GPa).

Investigation of ion fishbone stability on HL-2A using NIMROD

Numerical calculations of resistive internal kink mode with effects of energetic particles (EPs) on HL-2A have been performed using the hybrid kinetic-MHD model inplemented in the NIMROD code. The m/n = 1/1 resistive internal kink mode is unstable in MHD limit. However, with kinetic effects of energetic ions, a fishbone mode is excited with mode frequency around 10 kHz. We calculate the impact of resistivity on the growth rate and frequency of ion fishbone mode, and the results are in good agreement with the analytic solutions, which are obtained by solving the fishbone dispertion relation including resistivity effect. The effects of βfrac and cut off velocity of EP on fishbone mode are calculated in detail, where βfrac is the ratio of EP pressure to background plasma pressure. This work presents a clear explanation of the stabilizing effect of ECRH on ion fishbone, which is first observed on HL-2A.

Heavy ion ranges from first-principles electron dynamics

The effects of incident energetic particles, and the modification of materials under irradiation, are governed by the mechanisms of energy losses of ions in matter. The complex processes affecting projectiles spanning many orders of magnitude in energy depend on both ion and electron interactions. Developing multi-scale modeling methods that correctly capture the relevant processes is crucial for predicting radiation effects in diverse conditions. In this work, we obtain channeling ion ranges for tungsten, a prototypical heavy ion, by explicitly simulating ion trajectories with a method that takes into account both the nuclear and the electronic stopping power. The electronic stopping power of self-ion irradiated tungsten is obtained from first-principles time-dependent density functional theory (TDDFT). Although the TDDFT calculations predict a lower stopping power than SRIM by a factor of three, our result shows very good agreement in a direct comparison with ion range experiments. These results demonstrate the validity of the TDDFT method for determining electronic energy losses of heavy projectiles, and in turn its viability for the study of radiation damage.

Analysis of the MHD stability and energetic particles effects on EIC events in LHD plasma using a Landau-closure model

The aim of this study is to perform a theoretical analysis of the magnetohydrodynamic (MHD) stability and energetic particle effects on a LHD equilibria, calculated during a discharge where energetic-ion-driven resistive interchange mode (EIC) events were triggered. We use the reduced MHD equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particles species, including the effect of the acoustic modes, multiple energetic particles (EP) species, helical couplings and helically trapped EP. We add the Landau damping and resonant destabilization effects using a closure relation. The simulations suggest that the helically trapped EP driven by the perpendicular neutral beam injector (NBI) further destabilizes the 1/1 MHD-like mode located at the plasma periphery (r/a = 0.88). If the of the EP driven by the perpendicular NBI is larger than 0.0025 a 1/1 EIC with a frequency around 3 kHz is destabilized. If the effect of the passing EP driven by the tangential NBI is included on the model, any enhancement of the injection intensity of the tangential NBI below leads to a decrease of the instability growth rate. The simulations indicate that the perpendicular NBI EP is the main driver of the EIC events, as it was observed in the experiment. If the effect of the helical couplings are added in the model, an 11/13 EIC is destabilized with a frequency around 9 kHz, inward shifted (r/a = 0.81) compared to the 1/1 EIC. Thus, one possible explanation for the EIC frequency chirping down from 9 to 3 kHz is a transition between the 11/13 to the 1/1 EIC due to a weakening of the destabilizing effect of the high n modes, caused by a decrease of the EP drive due to a loss of helically trapped EP or a change in the EP distribution function after the EIC burst. The experimental data during the EIC bursting phase shows a complex mode structure and an inward shift of the instability, although no direct evidence of the proposed transition has been observed yet.