Ion Cyclotron Emission Research Articles

Experimental identification of Ion Cyclotron Emission on HL-2A using YOLO neural network algorithm

Abstract Identification of magnetohydrodynamics (MHD) instabilities with neural networks has been extensively applied in the research of magnetically controlled fusion plasmas. Ion Cyclotron Emission (ICE) is a potential fast ion diagnostic method in burning plasmas. To assess ICE as a fast ion diagnostic for ITER, real-time identification of ICE is required in the fast ion diagnostic flow. In the present work, we employed YOLO (You Only Look Once) to identify core and edge ICE in a large labeled database of HL-2A discharges, achieving a precision of 85.4% and a recall rate of 77.3%. Subsequent improvements to the YOLO model resulted in a noteworthy 8.3% increment in the recall rate. The developed identification method demonstrates significant potential for real-time application in identifying MHD instabilities.

The consequences of tritium mix for simulated ion cyclotron emission spectra from deuterium-tritium plasmas

Abstract Measurements of ion cyclotron emission (ICE) are obtained from most large magnetically confined fusion (MCF) plasma experiments, and may be used in future to quantify properties of the fusion-born alpha-particle population in deuterium-tritium (DT) plasmas in ITER. ICE is driven by spatially localised, strongly non-Maxwellian, minority energetic ion populations which relax collectively under the magnetoacoustic cyclotron instability (MCI). ICE spectral peaks are typically observed at, near, or separated by, integer harmonics of the energetic ion cyclotron frequency. Here, for the first time, we study how simulated ICE spectra from DT plasmas vary with tritium concentration and compare with an observed JET DT ICE spectrum. We incorporate a population of thermal tritons, in addition to thermal deuterons and an energetic minority population of fusion-born alpha-particles, in simulations with the kinetic particle-in-cell (PIC) code EPOCH. This code has previously been used extensively for interpretation of ICE observations in terms of the self-consistent gyro-resolved collective Maxwell-Lorentz dynamics of tens of millions of simulation particles. Our simulation parameters are relevant to the ICE-generative outer midplane edge region of JET DT plasma 26148, which included 11% tritium. Quantifying the variation of simulated ICE power spectra with tritium concentration reveals that our simulation with 11% tritium concentration most accurately represents the observed ICE spectrum from this plasma. This outcome is encouraging for the diagnostic application of ICE to fusion plasmas containing two thermal ion species, including future DT plasmas.

Detection of ion cyclotron emission by using an ion cyclotron range of frequency antennas-based diagnostic system in experimental advanced superconducting tokamak

In the experimental advanced superconducting tokamak (EAST), a novel ion cyclotron range of frequency (ICRF) antenna-based diagnostic system is designed to measure ion cyclotron emission (ICE) driven by high-energy ions. The diagnostic system includes ICRF antenna straps, a three-tune impedance matching system, a coaxial switching system, a direct current block, and a data acquisition and storage system. Using the coaxial switching system, the ICRF antenna can be switched from the heating mode to the coupling mode between two discharges. In the 2023 EAST experiment campaign, core ICE was observed using the ICRF antenna-based diagnostic system during neutron beam injection heating, and the obtained results agreed well with the signal detected by the previous high-frequency B-dot probe-based diagnostic system.

Development and implementation of ion cyclotron emission diagnostic system on J-TEXT tokamak

An ion cyclotron emission (ICE) diagnostic based on a magnetic probe has been designed and installed on the J-TEXT tokamak. This ICE diagnostic is used to study the high-frequency magnetic fluctuations driven by energetic particles. The signals detected by the ICE probe are divided into two equal power signals, which enter the two systems respectively. The first system converts high-frequency signals into power signals in different frequency bands using bandpass filters and power detectors. This reduces the amount of data to be processed and allows the data to be uploaded in time after a discharge. Thus, the approximate time resolution, frequency range (10 MHz-150 MHz), and power of the high-frequency signals can be quickly judged. The second system uses low-pass filters and an oscilloscope to collect data by Nyquist sampling law, so that the frequency spectrum of signals with frequencies from 7 kHz to 200 MHz can be analyzed. Therefore, if significant signals are detected in the first system on the data server, corresponding stored data can be extracted from the second system for large-scale data processing. This structure enhances the efficiency of data storage and analysis. The high-frequency magnetic fluctuations driven by fast electrons or runaway electrons under ohmic heating have been measured on J-TEXT, which verified the feasibility of the system.

Predicting ion cyclotron emission from neutral beam heated plasmas in Wendelstein7-X stellarator

Measurements of ion cyclotron emission (ICE) are planned for magnetically confined fusion plasmas heated by neutral beam injection (NBI) in the Wendelstein 7-X stellarator (W7-X). Freshly injected NBI ions in the edge region, whose velocity-space distribution function approximates a delta-function, are potentially unstable against the magnetoacoustic cyclotron instability (MCI), which could drive a detectable ICE signal. Prediction of ICE from NBI protons in W7-X hydrogen plasmas is challenging, owing to the low ratio of the ions’ perpendicular velocity to the local Alfvén speed, v⊥(NBI)/VA≃0.14 . We address this from first principles, using the particle-in-cell kinetic code EPOCH. This self-consistently solves the Lorentz force equation and Maxwell’s equations for tens of millions of computational ions (both thermal majority and energetic NBI minority) and electrons, fully resolving gyromotion and hence capturing the cyclotron resonant phenomenology which gives rise to ICE. Our simulations predict an ICE signal which is predominantly electrostatic while incorporating a significant electromagnetic component. Its frequency power spectrum reflects novel MCI physics, reported here for the first time. The NBI ions relaxing under the MCI first drive broadband field energy at frequencies a little below the lower hybrid frequency ωLH, across the wavenumber range kωc/VA=40 –60, where ω c and VA denote ion cyclotron frequency and Alfvén velocity. Nonlinear coupling between these waves then excites spectrally structured ICE with narrow peaks, at much lower frequencies, typically the proton cyclotron frequency and its lower harmonics. The relative strength of these peaks depends on the specifics of the NBI ion velocity-space distribution and of the local plasma conditions, implying diagnostic potential for the predicted ICE signal from W7-X.

ICE modes with high wave numbers

This work deals with high-k modes of the ion cyclotron emission (ICE)—the modes with frequencies close to ion cyclotron harmonics and wave numbers well exceeding those of fast magnetoacoustic modes (FMM). These modes exist due to finite Larmor radius of the ions. They are responsible for the ICE recently observed in the DIII-D tokamak [N. A. Crocker et al., Nucl. Fusion 62, 026023 (2022)]. The work is aimed at knowing the radial structure of high-k eigenmodes. Its analysis is relied on a comprehensive study carried out by employing a general dispersion equation for traveling waves, which allowed to group modes and determine conditions of their existence. Differential equations for various eigenmodes (standing waves) were obtained and analyzed. In particular, electrostatic modes and electromagnetic modes with various ratios of the longitudinal phase velocity to the electron thermal velocity were considered. A numerical code solving these equations was developed. It is found that high-k modes with sufficiently large poloidal mode numbers are located at the plasma periphery, and their amplitudes in the plasma core are very small. It is concluded that high-k modes are not an exceptional case, although presumably FMMs can explain most ICE experiments. They exist in plasmas with various magnitudes of β, representing different wave branches associated with finite Larmor radius of the ions. The mentioned experiment on DIII-D is discussed.

Interpretation of ion cyclotron emission from sub-Alfvénic beam-injected ions heated plasmas soon after L-H mode transition in EAST

Ion cyclotron emission (ICE) at deuterium ion cyclotron harmonics, driven by sub-Alfvénic beam-injected deuterium ions, has been observed by the high-frequency B-dot probe in the EAST tokamak. The origin of ICE shifts from the plasma core to the plasma edge soon after an L-H mode transition, where the beam-injected deuterium ions have a relatively peak bump-on tail structure in the energy direction and a very intense pitch angle anisotropy. Based on the fast ion distribution function obtained from the TRANSP/NUBEAM code, together with a linear analysis theory of magnetoacoustic cyclotron instability (MCI), the growth rates of MCI could be calculated. It is shown that MCI, resulting in the generation of obliquely propagating fast Alfvén waves at deuterium ion cyclotron harmonics, can occur under such conditions. And the temporal evolution of the MCI growth rate closely follows that of the observed ICE amplitude in the EAST.

A linear parameters study of ion cyclotron emission using drift ring beam distribution

Ion Cyclotron Emission (ICE) holds great potential as a diagnostic tool for fast ions in fusion devices. The theory of Magnetoacoustic Cyclotron Instability (MCI), as an emission mechanism for ICE, states that MCI is driven by a velocity distribution of fast ions that approximates to a drift ring beam. In this study, the influence of key parameters (velocity spread of the fast ions, number density ratio, and instability propagation angle) on the linear MCI is systematically investigated using the linear kinetic dispersion relation solver BO (Xie 2019 Comput. Phys. Commun. 244 343). The computational spectra region considered extends up to 40 times the ion cyclotron frequency. By examining the influence of these key parameters on MCI, several novel results have been obtained. In the case of MCI excited by super-Alfvénic fast ions (where the unique perpendicular speed of fast ion is greater than the perpendicular phase velocity of the fast Alfvén waves), the parallel velocity spread significantly affects the bandwidth of harmonics and the continuous spectrum, while the perpendicular velocity spread has a decisive effect on the MCI growth rate. As the velocity spread increases, the linear relationship between the MCI growth rate and the square root of the number density ratio transitions to a linear relationship between the MCI growth rate and the number density ratio. This finding provides a linear perspective explanation for the observed linear relation between fast ion number density and ICE intensity in JET. Furthermore, high harmonics are more sensitive to changes in propagation angle than low harmonics because a decrease in the propagation angle alters the dispersion relation of the fast Alfvén wave. In the case of MCI excited by sub-Alfvénic fast ions (where the unique perpendicular speed of fast ion is less than the perpendicular phase velocity of the fast Alfvén waves), a significant growth rate increase occurs at high harmonics due to the transition of sub-Alfvénic fast ions to super-Alfvénic fast ions. Similarly, for MCI excited by greatly sub-Alfvénic fast ions (where the unique perpendicular speed of fast ion is far less than the perpendicular phase velocity of the fast Alfvén waves), the growth rate at high harmonics also experiences a drastic increase compared to the low harmonic, thereby expanding the parameter range of the velocity spread.

Identification of core ion cyclotron instabilities on HL-2A tokamak

Instabilities in multiplies of ion cyclotron frequency range are identified and termed as core ion cyclotron emission (ICE) in recent HL-2A neutral beam injection heated experiments. Characteristics of the core ICE are presented, including frequency dependence and harmonics features. The detected frequencies are found to agree well with the multiplies of the deuterium cyclotron frequency around the magnetic axis. Additionally, the core ICE exhibits a predominantly compressional property. Observations of distinct spectrum features and individual excitation of each harmonic have demonstrated that the core ICE harmonics are independent multiple modes. Notably, the variation of plasma current is a necessary condition for exciting the 4th harmonic ICE individually. The results suggest that the drive mechanism of core ICE varies between the different frequency ranges.

Velocity-space sensitivity and inversions of synthetic ion cyclotron emission

This paper introduces a new model to find the velocity-space location of energetic ions generating ion cyclotron emission (ICE) in plasmas. ICE is thought to be generated due to inverted gradients in the v⊥ direction of the velocity distribution function or due to anisotropies, i.e., strong gradients in the pitch direction. Here, we invert synthetic ICE spectra generated from first principles PIC-hybrid computations to find the locations of these ICE-generating ions in velocity space in terms of a probability distribution function. To this end, we compute 2D ICE weight functions based on the magnetoacoustic cyclotron instability, which reveals the velocity-space sensitivity of ICE measurements. As an example, we analyze the velocity-space sensitivity of synthetic ICE measurements near the first 15 harmonics for plasma parameters typical for the Large Helical Device. Furthermore, we investigate the applicability of a least-square subset search, Tikhonov regularization, and Lasso regularization to obtain the locations in velocity space of the ions generating the ICE.

On equations for ion cyclotron modes in ‘warm’ bounded plasmas

An equation describing eigenmodes in the ion cyclotron frequency range in ‘warm’ bounded plasmas, i.e. eigenmodes which are absent in the two-fluid model but exist in kinetic theory due to finite Larmor radius of the ions, is derived for the first time. It is valid for electrostatic modes but the developed approach is generic. Calculations are carried out for two cases: first, for a homogeneous magnetic field; second, taking into account the effects of toroidicity in tokamaks. It is found that, in general, equations for eigenmodes in warm plasmas do not reduce to second-order differential equations (in contrast to those which are usually used to describe the radial structure of eigenmodes in fusion devices). The study of modes in warm plasmas is of interest, in particular, in connection with the recent observations of superthermal ion cyclotron emission in the NSTX-U spherical torus and DIII-D tokamak, which can be hardly explained by conventional theories employing fast magnetoacoustic modes.

Mechanism for Collective Energy Transfer from Neutral Beam-Injected Ions to Fusion-Born Alpha Particles on Cyclotron Timescales in a Plasma

Helium ash alpha particles at ∼100 keV in magnetically confined fusion plasmas may have the same Larmor radius, as well as cyclotron frequency, as the energetic beam-injected deuterons that heat the plasma. While the velocity-space distribution of the helium ash is monotonically decreasing, that of the energetic deuterons is a delta function in the edge plasma. Here we identify, by means of first principles particle-in-cell computations, a new physical process by which Larmor radius matching enables collective gyroresonant energy transfer between these two colocated minority energetic ion populations, embedded in majority thermal plasma. This newly identified nonlinear phenomenon rests on similar underlying physics to widely observed ion cyclotron emission from suprathermal minority ion populations.

Doublet splitting of fusion alpha particle driven ion cyclotron emission

Ion cyclotron emission (ICE) originating from confined populations of fast ions in toroidal fusion plasmas is an important non-invasive, passive diagnostic for current and next-generation devices. The ability to model ICE signals accurately is an essential step towards inferring the characteristics of confined energetic alpha particles or fast neutral beam ions and background plasma. In this paper, the linear growth rates of the magnetoacoustic cyclotron instability, which is the leading explanation for ICE, are calculated in high-resolution 2D space for parameters corresponding to JET Pulse No. 26148. These calculations further highlight the origin of doublet splitting of the cyclotron harmonics and the perceived need to invoke nonlinear interactions with 1D3V particle-in-cell studies to account for lower harmonics with substantial amplitudes. The inclusion of signals of ICE from a well-resolved range of propagation angles can also account for these effects.

Mode structure measurements of ion cyclotron emission and sub-cyclotron modes on DIII-D

We report mode structure measurements of coherent ion cyclotron emission (ICE) and sub-cyclotron modes on DIII-D. Through a dedicated experiment, we aimed to characterize a variety of modes via the upgraded ICE diagnostic in both L- and H-mode plasmas. In the L-mode plasmas, autopower spectrum peaks at harmonics of the ion cyclotron frequency f ci were observed and are localized in the core of the plasma. Sub-cyclotron modes (f ∼ 0.5f ci) were also observed in L-mode plasmas when the toroidal magnetic field strength was lowered from B T = 2.17 T. In H-mode plasmas, many ICE harmonics localized to the plasma edge were observed, with some exceeding the nominal ICE diagnostic bandwidth of f ∈ [0, 100] MHz. Polarization estimates made using an orthogonal pair of vertical and horizontal magnetic pickup loops on the outer wall of the machine indicate that ICE and sub-cyclotron modes have compressional polarization at the plasma edge, the latter being consistent with simulation efforts for comparable DIII-D plasmas. For all modes, the same harmonics are observed on both the centerpost and outer wall loops, indicating that ICE is poloidally extended. Finally, toroidal mode numbers were calculated using three outer wall loops for both core ICE and sub-cyclotron modes in L-mode plasmas. The sub-cyclotron case served as a benchmark for our calculation method, with measured numbers of n ∈ [−24, −18] roughly agreeing with heuristic estimates of n ∈ [−20, −13]. Core ICE mode numbers were measured to be n ∈ [−10, 5], which is comparable to measurements made on other machines.

A Pockels cell enabled heterodyne phase contrast imaging diagnostic for detection of ion cyclotron emission.

This work describes a novel optical heterodyne detection scheme that significantly extends the frequency response of the phase contrast imaging method to detect electron density fluctuations in tens of megahertz frequency range. The system employs a variable frequency electro-optic modulator to allow operation at any frequency in the range 10-40MHz without the need to realign the system. The frequency coverage of the system makes it suitable to measure the radial structure of the electron density component of ion cyclotron emission on devices having confining magnetic field of a few tesla, thus extending the purely temporal measurements provided so far by magnetic probes.

Verification of a fully kinetic ion model for electromagnetic simulations of high-frequency waves in toroidal geometry

For the study of high-frequency electromagnetic waves in tokamaks, an electromagnetic simulation model, in which the ion dynamics is described by a six-dimensional Vlasov equation and the electron dynamics is described by a drift kinetic equation, is formulated and implemented in the global gyrokinetic toroidal code (GTC). Analytic dispersion relations are derived in reduced systems and compared with various theories to verify the model. Linear simulations of a generalized ion Bernstein wave and ion cyclotron emission are verified by comparing the GTC simulation results with analytic dispersion relation theory and magnetoacoustic cyclotron instability theory, respectively, in cylindrical geometry.

First observation and interpretation of spontaneous collective radiation from fusion-born ions in a stellarator plasma

During bursty MHD events, transient ion cyclotron emission (ICE) is observed from deuterium plasmas in the large helical device (LHD) heliotron-stellarator. Unusually, the frequencies of the successive ICE spectral peaks are not close to integer multiples of the local cyclotron frequency of an energetic ion population in the likely emitting region. We show that this ICE is probably driven by a subset of the fusion-born protons near their birth energy MeV. This subset has a kinetic energy component parallel to the magnetic field, , significantly greater than its perpendicular energy , for which , the Alfvén speed. First principles computations of the collective relaxation of this proton population, within a majority thermal deuterium plasma, are carried out using a particle-in-cell approach. This captures the full gyro-orbit kinetics of all ions which, together with an electron fluid, evolve self-consistently with the electric and magnetic fields under the Maxwell–Lorentz equations. The simulated ICE spectra are derived from the Fourier transform of the fields which are excited. We find substantial frequency shifts in the peaks of the simulated ICE spectra, which correspond closely to the measured ICE spectra following the resonance condition for nth proton harmonic. This suggests that the transient ICE in LHD is generated by the identified subset of the fusion-born protons, relaxing under the magnetoacoustic cyclotron instability. So far as is known, this is the first report of a collective radiation signal from fusion-born ions in anon-tokamak magnetically confined plasma. Disambiguation between two or more energetic ion species that could potentially generate complex observed ICE spectra is an increasing challenge, and the results and methodology developed here will assist this. Our approach is also expected to be relevant to ICE driven by ion beams with lower parallel velocities, for example in cylindrical plasma experiments.

Novel internal measurements of ion cyclotron frequency range fast-ion driven modes

Novel internal measurements and analysis of ion cyclotron frequency range fast-ion driven modes in DIII-D are presented. Observations, including internal density fluctuation () measurements obtained via Doppler backscattering, are presented for modes at low harmonics of the ion cyclotron frequency localized in the edge. The measurements indicate that these waves, identified as coherent ion cyclotron emission (ICE), have high wave number, k ⊥ ρ fast ≳ 1, consistent with the cyclotron harmonic wave branch of the magnetoacoustic cyclotron instability, or electrostatic instability mechanisms. Measurements show extended spatial structure (at least ∼1/6 the minor radius). These edge ICE modes undergo amplitude modulation correlated with edge localized modes (ELM) that is qualitatively consistent with expectations for ELM-induced fast-ion transport.

Development of the ion cyclotron emission diagnostic on the HL-2A tokamak

An ion cyclotron emission (ICE) diagnostic, which is based on a B-dot probe, has been recently designed and installed on HL-2A tokamak. The diagnostic is used to study various high-frequency magnetic field fluctuations which can be excited by energetic ions and runaway electrons in the plasma. The ICE diagnostic on HL-2A includes a high-frequency B-dot probe, direct current (DC) block, radio frequency splitters, filter bank and power detectors. The filter bank is composed of 16 channels filters, with the center frequency covering from 10 to 160 MHz, 10 MHz step length and 8 MHz bandwidth. The log detectors with a large dynamic range (from −80 dBm to −20 dBm) are used to detect the bandpass power. Test results of the B-dot probe, filters and power detectors are shown. The signals can also be sampled with a fast analog-to-digital converter with a 14-bit depth, 100 MHz bandwidth and 250 MSample/s sampling rate.

Identification of slow-wave ion cyclotron emission on JT-60U

Ion cyclotron emissions (ICEs) are generally considered to be emissions of fast waves driven by velocity-space instabilities due to fast ions. However, on JT-60U, measured properties of ICEs observed at frequencies lower than the bulk ion cyclotron frequency (L-ICEs) do not match the fast-wave dispersion relation, indicating the possibility of L-ICEs being slow-wave emissions. In this paper, we show that L-ICE observation on JT-60U can be explained in terms of a slow-wave emission driven by fast ions, i.e. a slow-wave ICE. To investigate whether the slow wave can be driven by fast ions, its linear growth rate has been calculated in a typical discharge using a wave dispersion code. This wave dispersion code is capable of performing calculations with an arbitrary spatially-localized velocity distribution. For the growth rate calculation, we have directly used fast ion velocity distributions evaluated by an orbit following Monte-Carlo simulation. It is found that negative-ion-source neutral beam (N-NB) injected fast ions can destabilize the slow wave. Its frequency and wavenumber at high linear growth rates are close to experimental observations of L-ICE. In addition, observed L-ICE frequencies agree reasonably well with frequencies based on the slow-wave dispersion relation at measured toroidal wavenumbers in several discharges. Moreover, the observed frequencies also agree with the Doppler-shifted cyclotron resonance frequencies for the N-NB injected ion. Therefore, L-ICE in JT-60U is identified as a slow-wave ICE driven by the N-NB injected fast ions and distinguished from other ICEs through a measurement of its frequency and wavenumber.