# Main Challenges of Heating Plasma with Waves at the Ion Cyclotron Resonance Frequency (ICRF)

Guillaume Urbanczyk

► To cite this version:

Guillaume Urbanczyk. Main Challenges of Heating Plasma with Waves at the Ion Cyclotron Resonance Frequency (ICRF). Advances in Fusion Energy Research - From Theory to Models, Algorithms, and Applications, IntechOpen, 2022, 10.5772/intechopen.105394. [https://doi.org/10.5772/intechopen.105394](https://doi.org/10.5772/intechopen.105394)

HAL Id: hal-04242138

https://hal.science/hal-04242138v1

Submitted on 4 Apr 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Main Challenges of Heating Plasma with Waves at the Ion Cyclotron Resonance Frequency (ICRF)

Guillaume Urbanczyk

## Abstract

Of all the techniques used for heating plasmas in fusion devices, waves in the Ion Cyclotron Resonance Frequency (ICRF ~ MHz) continue to be exceptionally advantageous and unique insofar as it enables to deposit of power directly on ions in the core, significantly enhancing fast ion population together with fusion reaction products. However, because of the multiple inherent challenges—such as matching robustness, antenna design, wave coupling efficiency, wave propagation, wave absorption, and plasma surface interactions due to radiofrequency (RF) sheath excitation—ICRF is often one of the most complex heating systems to implement successfully. This chapter provides a brief introduction of these challenges and their respective underlying physics, together with examples of both simulations and experimental results from various tokamaks around the world. Finally, ICRF advantages and applications on present and future devices and perspectives of technological solutions are discussed and summarized.

**Keywords:** ICRF, wave coupling, RF sheath, plasma surface interactions

## 1. Introduction

The ICRF electrostatic wave is first produced by a generator called a tetrode, which is essentially a multistage amplifier of the power composed of a cascade of electron tubes. Each stage of which increases the power by approximately an order of magnitude, from watts to megawatts [1].

High-power waves enter the coaxial line of characteristic impedance $Z_0$, that matches the output impedance of the generator; typically, 50Ω or 30Ω as in the case of the Full-Tungsten Environment Superconducting Tokamak (WEST) [2] which **Figure 1** shows an overview of the ICRF system with all key elements.

1. Electrostatic waves then propagate towards the antenna and cross at some point the so-called *matching* unit which ensures a smooth transition between the impedance of the coaxial line $Z_0$ and the one of the antenna straps $Z_S$ to avoid
<figure><img src="image_2.png" /><figcaption>**Figure 1.** ICRF system seen in 3D (a) and from profile (b), with a schematic of the key components of the RF circuit (c) from the generator to the antenna (case of WEST tokamak) adapted from [3].</figcaption></figure>

power reflection. If successfully matched, the wave induces a current on the strap at the extremity of the transmission line.

2. This current oscillating at the ion cyclotron frequency (one the order of MHz), results in the excitation of electromagnetic ICRF waves and their *coupling* to the low-density edge plasma.

3. These *waves* must then propagate from the edge up to the resonance layer in the core, where they can be *absorbed by the ions*. Energy is then spread to all charged particles through collisions, resulting in *plasma heating*.
These three aspects, namely the *matching*, wave *coupling*, and *absorption*, are the key points mandatory for the successful operation of an ICRF system. Section 2 provides the physics background for each of these steps, while section 3 explains what happens in practice due to all of the imperfect conditions and introduces the challenges inherent in ICRF operation; power reflection, arcing, Radiofrequency (RF) sheath excitation, erosion, impurity generation, and plasma contamination. Despite all of these challenges, the fourth and final sections present the key assets that still make ICRF an attractive auxiliary heating system and some perspectives on solutions worthy of further investigation.

## 2. The key 3 steps of ICRH (Ion Cyclotron Resonance Heating)

### 2.1 RF circuit matching

Matching the RF circuit consists of building smooth transitions between the purely real impedance of the transmission line ($Z_0 = R_0 \approx 10 \rightarrow 50\Omega$) and the complex impedance of the antenna straps ($Z_S = R_C + jX_S$) which depends on the plasma. When the system is perfectly matched, the electrostatic wave flows from the generator to the straps without seeing any obstacle or being subject to any reflection. This can be achieved in many ways, relying on various components each with a variable impedance such as:

* *Single stub* also called “trombones” [4] in the Joint European Torus (JET) ITER-Like Antenna [5, 6] (ITER=International Tokamak Experimental Reactor), Alcator C-Mod [7] and the National Spherical Torus Experiment (NSTX) [8]

* *Double stubs loaded by shorted tunable capacitors* on TEXTOR [9] and analyzed for ASDEX-Upgrade [10]

* *Liquid stubs* arranged in a triple stub configuration in the Large Helical Device (LHD) [11] and the Experimental Advanced Superconducting Tokamak (EAST) [12]. The reactance of the short-circuited stubs of fixed geometrical length is varied by changing the height of the dielectric fluid (oil) inside the stubs

* *Fast Ferrite Tuners*: the reactance of a matching stub of fixed geometrical length is here modified by changing the magnetic characteristics of the ferrites which are disposed inside the stubs which were used in Alcator C-Mod [7] and ASDEX Upgrade [13]

* *Sliding Impedances* [14]

* *Single stubs in the Conjugate-T configuration* [15]

* *Tunable line-stretchers* (or *phase shifters*) in JET [6]

* *Tunable vacuum capacitors* in JET ITER-like antenna [16] and WEST ICRF antennas [17]. These are complex components ($Z_C = X_L - jX_C$), which self-inductance ($X_L$), and capacitance ($X_C$) can be tuned to match any load at the strap (cf. purple elements in **Figure 1**)
If the antenna is facing vacuum, the load is stable as it basically only corresponds to ohmic losses in the antenna. In this case, almost perfect matching can be achieved with the resistive part of the load remaining small ($R_c \approx 0.1\Omega$). In front of a plasma though, the load is larger than vacuum, so that $R_c$ can increase by up to two orders of magnitude, in correlation with the efficiency by which the power is coupled.

In this case, however, the load (plasma) can change rapidly. This means $Z_s$ is seldom if ever perfectly matched to $Z_0$, and that the power launched by the generator ($P_{in}$) is hardly exactly equal to the one coupled by the antenna to the plasma:

$$P_{coupled} = P_{in}(1 - |\Gamma|^2) \quad (1)$$

where $\Gamma$ if the reflection coefficient defined as $\Gamma = (Z_s-Z_0) / (Z_s + Z_0)$.

The fraction of uncoupled power is therefore reflected in the transmission line, which can damage the generator if it is too large. This is avoided by introducing a quarter wavelength phase-shift in the matching circuit—for instance by tuning the height of the oil in EAST stub tuner or by using an impedance transformer in WEST— so that in the section between the generator and the matching circuit, only forward waves exist.

In the section between the matching circuit and the antenna though, coexist waves travelling both forward and backward, which superimpose on one another and give rise to a standing wave pattern. The maximum and minimum voltages of such waves are defined as follows:

$$V_{max} = (1 + |\Gamma|) \sqrt{2Z_0 P_{in}} = \sqrt{2Z_0 P_{coupled} V_{SWR}} \quad (2)$$

$$V_{min} = (1 - |\Gamma|) \sqrt{2Z_0 P_{in}} = \sqrt{2Z_0 P_{coupled} / V_{SWR}} \quad (3)$$

And the ratio of these voltages is the so-called Voltage Standing Wave Ratio:

$$V_{SWR} = \frac{1 + |\Gamma|}{1 - |\Gamma|} \quad (4)$$

It follows that VSWR grows from 1 in the ideal perfectly matched “flat-line” case ($V_{max} = V_{min}$), to infinity when almost all the power is reflected. If the VSWR is too large, the difference of electrical potentials in the lines becomes so high that there is a significant risk of arcing. To prevent such deleterious events, operational safety limits are set on the voltage ($V_{max}$) and the current ($I_{max}$), above which the power shuts down.

$$I_{max} = (1 + |\Gamma|) \sqrt{2Y_0 P_{in}} = \sqrt{2Y_0 P_{coupled} V_{SWR}} \quad (5)$$

with $Y_0=1/Z_0$. We now understand how crucial it is to transport the power along well-matched lines. It naturally follows that minimizing length of the lines with these undesirable effects by placing the matching system as close as possible to the antenna, will be a key aspect of the system's efficiency (cf. section 2.5 of [18]).

In addition, decouplers are used to prevent the mutual coupling of adjacent lines which influence each other due to the proximity of the straps to one another. This effect can make a source behave as a receiver and endanger its generator. Decouplers are used in DIII-D [19], EAST [20], Alcator C-Mod [7], NSTX [8] and will be used on ITER [21]. This is also important to equalize the voltages at the inputs of the array and
guarantee a homogeneous excitation over the whole surface of the antenna's front face, which can otherwise have deleterious effects on local fields and RF sheath excitation. Strap-decoupling can also be improved by separating them with a septum in the antenna box to reduce mutual influence.

Excellent explanations of the topic can further be found in the fourth chapter of [18] with more detailed calculations and specific applications to the case of WEST. Once well matched, the system and the antenna are then ready to launch waves which coupling to the plasma remains to be optimized.

## 2.2 Wave coupling

The front of the antenna faces the plasma inside the vessel. As the power emitted from the generator reaches the antenna, it induces an oscillating current $I_{RF}$ along the strap as highlighted in red in **Figures 1b and 2**, which excites ICRF waves at the edge of the plasma, as efficiently as the load or the coupling resistance ($R_c$) is high.

<figure><img src="image_3.png" /><figcaption>**Figure 2.**
Poloidal cut of a tokamak vacuum vessel, with ICRF waves excited by the antenna, propagating in the well-confined region of the plasma, partially absorbed and mode converted near the resonance layer.</figcaption></figure>
$$P_{\text{coupled}} = \frac{1}{2} R_c I_{\text{RF}}^2 \qquad (6)$$

The problem in coupling ICRF waves at the edge, lies in the fact that to avoid exposure to dense and hot plasma—which would induce unacceptably large heat loads—antennas are retracted away and sit in the Scrape-Off Layer (SOL), a low-density plasma (pink regions in **Figures 2 and 3**) where the Fast Wave (FW) is generally not propagative but evanescent ($k_{\perp}^2 < 0$). A wave is evanescent when the density is smaller than its cutoff density ($n_{\text{co}}$). Therefore, the larger the distance from the strap to the cutoff layer ($D_{\text{Strap-co}}$), the lower the coupling efficiency and vice versa. To give a rough idea, 1 cm increase of $D_{\text{Strap-co}}$ generally results in about 20% drop of the coupling efficiency, meaning that about 1/5 of the power vanishes every centimeter until reaching densities above $n_{\text{co}}$. In addition, the nature of the spectrum excited by the antenna, according to its geometry (straps width, height and spacing), and in particular, the value of the parallel wave vector ($k_{\parallel}$) corresponding to the main components of the spectrum, plays a key role in the coupling process:

$$R_c = \frac{Z_0}{\text{VSWR}} \sim R_0 \exp(-2 \langle k_{\parallel} \rangle D_{\text{Strap-co}}) \qquad (7)$$

One can observe, in **Figure 3**, how this parameter typically influences the fast-wave cutoff density which drops by over an order of magnitude as $k_{\parallel}$ decreases from 18 down to $6 \text{ m}^{-1}$.

ICRF waves comport two modes called *fast* and *slow* waves in reference to their respective group velocity (details of mathematical calculations can be found in the second chapter of [22]):

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Properties of 120MHz ICRF fast (FW—dashed lines) and slow (SW—solid line) wave modes in a realistic fusion plasma with similar proportions of deuterium and tritium confined by a magnetic field of 10 tesla. Waves are propagative when the sign of $k_{\perp}^2$ is positive, and evanescent elsewhere. Cutoff densities of each branch are written explicitly on the graph.</figcaption></figure>
* Fast Waves (FW) are evanescent at low densities typical of the plasma edge and propagative for densities characteristic of the plasma core (Figure 3). FW have wavelengths on the order of few meters and an electric field perpendicular to the magnetic field. Therefore, when talking about coupling maximization, it is common to implicitly refer to this mode, because it is the one suited to heating the ions which rotate perpendicularly around the magnetic field lines.

* Slow Waves (SW) are propagative at very low densities characteristic of the outer most “far” SOL ($n_e < 5.2 \ 10^{17} \ m^{-3}$ in Figure 3) and evanescent above. SW have millimetric wavelengths and an electric field parallel to the magnetic field. Because the current straps are not perfectly perpendicular to the tilted magnetic field lines, SW are parasitically excited and are often responsible for exciting the RF-sheath, leading to deleterious plasma surface interactions. The role the faraday screen (Figure 1a) with bars in principle parallel to the magnetic field, is intended to locally increase the conductivity in the parallel direction to screen the SW in front of the straps. Unfortunately, we will see in the next section that this is not enough.

At this point, one should remember when hearing people say that they are trying to “optimize the coupling of an ICRF antenna”, that they are basically trying to maximize its coupling resistance ($R_c$), either by reducing $D_{\text{Strap-co}}$ or $k_{\parallel}$ (Eq. (7)).

$D_{\text{Strap-co}}$ can be minimized in different ways [23]:

* Put the antenna closer to the plasma, but this will come with stronger heat loads and usually impurity production

* Globally increase the density of the plasma, but this is limited and associated with an increase of the power lost in radiations

* Locally increase the density by fueling the plasma from valves magnetically connected to the antenna and ideally as close as possible to the middle plane [24, 25]

* Increase the frequency of the waves, which must remain compatible with an efficient scenario for heating the plasma (discussed in the next section)

$k_{\parallel}$ can be minimized at two different stages:

* While designing the antenna, basically by increasing the space between adjacent straps, but this often limited by the available space in the torus

* During the experiments, by reducing the phase difference of the currents on the adjacent straps ($\phi$). Most of the time however, this trick is first limited by the currents induced on the surrounding passive structures, which play an important role in plasma surface interactions [26]. These currents basically tend to add up when currents are in phase (monopole $\phi=0°$) instead of compensating each other for largest phasing (dipole $\phi=180°$) [27]. Secondly, the larger $k_{\parallel}$ the more efficient the wave absorption [28].

Once efficiently coupled to the plasma, wave absorption remains to be optimized.
## 2.3 Wave absorption

Waves that propagate towards the center of the plasma finally reach the resonance layer where ions frequency ($w_0$) corresponds to a multiple of the wave frequency ($n_{wci}$), enlightened by the vertical yellow line in Figure 2. In this region, the frequency of the wave coincides with the fundamental or harmonics of the cyclotron frequency of ion species to be heated. One can tune the wave frequency, the magnetic field profile, and the plasma composition to pick the desired heating scheme among a variety of techniques [29] by matching the Doppler-shifted wave frequency with the ion cyclotron harmonics:

$$ \omega = n\omega_{ci} + k_{//}v_{//},I \quad (8) $$

* Heating a minority at its fundamental resonance (n=0): such as hydrogen in a deuterium plasma (H-D) with a small concentration of hydrogen, typically $n_H/n_H+n_D<15\%$. So far, minority heating is the most routinely used in most devices. However, while it has long been applied on two-ion species plasmas, it has been found that it can also be applied with great efficiency in plasmas with three (or more) ion species.

* The three-ions scheme also consists in heating the fundamental of the minority species such as a little helium in a deuterium tritium fusion plasma (ex: ($He^3$)-DT), but here the resonance location of the minority must fall in between the resonances of both others majority species [30], where the amplitude of the left-handed polarization ($E_+$) can reach very high values and lead to strong diffusion (D) of the wave power:

$$ D\alpha|J_0|^2|E_+|^2 \quad (9) $$

* Heating a majority at its harmonic resonance (n>0): this can be an alternative to the minority heating typically if it proves challenging to accurately controlling the isotopic ratio. However, this scenario becomes truly effective only when a population of fast ions is created to improve power absorption. This process usually relies on harmonic heating. While generating a population of fast ions large enough to boost up the absorption, it can also be directly injected by a Neutral Beam Injector (NBI), hence the so-called synergy between NBI and ICRH (cf. experimental [31] and simulation [32] results in JET tokamak). Unfortunately, not only NBI is not available in all devices, but if the confinement is not good enough, for instance in case of low plasma current, this method is counteracted by fast-ions losses.

# 3. Challenges behind each step

## 3.1 Real-time matching

In practice, before main experiments start, ICRF operators usually “prepare the matching” and pre-set the parameters (ex: Oil level in the stub tuner in EAST, capacitance of the capacitors in WEST ...) with the antenna facing vacuum. Later, the matching is adjusted with the antenna facing the plasma, one of the many reasons why
campaigns always include a period for commissioning systems at the beginning. The system starts with low power in case of sources of mismatch or poor coupling efficiency to avoid the risk of arcing due to large fractions of power reflected, and gradually raise it up to power levels relevant for the experiments (MW order). One of the key goals during the experiments is to further improve plasma confinement, often leading to operate in H-mode [33]. H-mode plasmas represent a great challenge for all RF circuit matching, first because during the L-H transition, a pedestal forms at the edge where the density typically drops with sharper gradient compared to L-mode plasmas, lowering the coupling efficiency ($R_c$) sometimes by half. But also, the so-called Edge-Localized Modes (ELMs) induce important load variations through transient bursts, basically as large as their frequency is low (from kHz range for the smallest type-III ELMs down to tens of Hz for the largest type-I ELMs with $R_c$ suddenly increasing by factors up to 5 within less than a millisecond. One can further show for a strap matched by a 2-port matching unit, if the matching is not reconfigured, the VSWR will raise in similar proportions as $R_c$ (cf. calculation details in section 2.7 of [18]), risking provoking arcs either in the transmission lines [34] or in the antenna box [35]. All this motivates real-time control of the impedance matching that can be achieved with:

* Capacitors, like in Tore Supra Classical antennas [36] and JET-ILA [37]. Amplitudes of the incident and reflected voltages and their phase shift are typically measured either on the transmission lines or ideally as close as possible to the antenna, and mixed to provide error signals driving the capacitors

* Conjugate-T and similar principle as with capacitors [13]

* Fast Ferrite Tuners, like in Alcator C-mod antennas [7]

* Impedance matching by real-time controlling the *generator's frequency* [4]

* Decouplers and *double stub* in ITER [21]

Various algorithms can then be used to automatically tune parameters of interest for each component. However so far, no controller allows tackling sub millisecond variations typically induced by ELMs, hence the need for so-called *load resilient* matching schemes [38].

## 3.2 Load resilience

Without load resilience, an ICRF system would stop delivering power every time its VSWR would exceed some safety value, which would basically make it incompatible with steady operation during a discharge with ELMs. A very convincing comparison between resilient and non-resilient system can be found in Figure 13 of [13]. Among existing load-resilient schemes, one can quote:

* 3dB *hybrid coupler* diverting the reflected power to a dummy load, such as the lossy network in ASDEX-Upgrade [39]

* Conjugate-T concept, connecting pairs of straps, with various designs:
* *External* conjugate-T relying on line-stretchers in JET [13] and EAST [40]

* *Internal* conjugate-T relying on capacitors in JET-ILA [5], Tore Supra [41] and WEST [17]

* *Stub* conjugate-T in Alcator C-Mod [15]

More details can be found on the conjugate-T under section 4 of the fourth chapter of [18] with detailed mathematical calculations and applications to the internal conjugate-T of WEST and experimental results in [42].

* *Travelling wave antenna* concept [43] which in case of load drop, takes profit from the mutual coupling between straps to let the wave propagate along the array and recirculate, making it inherently load resilient. This concept has been used in DIIID [44, 45] for driving current and envisaged for heating plasmas in WEST [46]. While never tested for heating tokamak plasmas, the concept somehow lost its appeal, essentially due to the integration challenges it raises, that become even more problematic in the perspective of fusion reactors which must maximize the surface of the tritium breading blanket and therefore focus on heating systems with large power densities.

Note that in most of the cases, load resilient solutions are combined with matching units, so that the problem is most often tackled from the generator perspective, with first goal to deliver constant power despite load fluctuations at the antenna.

After having explained the importance of operating with a well-matched system, we will now place ourselves from the inner vessel perspective and further discuss the challenges that concern an ICRF antenna facing a fluctuating plasma. There the main goal will rather be to prevent the coupling resistance from dropping under some value, not only critical for the VSWR and arcing in the transmission lines, but also for the excitation of near fields and RF sheath.

### 3.3 RF-sheath and potential rectification

The formation of a sheath on any component facing the plasma is a natural phenomenon [47]. The sheath is a thin layer (usually few millimeters wide) which forms on the surface of all materials in contact with the plasma. Across this layer, a separation of the charges (therefore an electric field) makes it possible to preserve the ambipolarity of the fluxes of charged particles from the plasma towards a surface (cf. **Figure 4**). Since the mass of electrons is much smaller than that of ions, their speed is much greater. Therefore, when exposed to a plasma, the surface of materials becomes negatively charged, repelling electrons, and attracting ions approaching below a Debye length ($\lambda_{De} = \sqrt{\epsilon_0 k_B T_e / n_e e^2}$) with $\epsilon_0$ and $k_B$ respectively electric and Boltzmann constants, $T_e$ and $n_e$ respectively the electron temperature and density and the charge of the electron $e = 1.6 \times 10^{-19}$ coulombs.

The properties of the sheath can undergo substantial changes depending on several parameters that can typically influence its electric field and width $\delta_{sheath}$:

$$ \delta_{sheath} = \lambda_{De} \left( \frac{e \cdot V_{DC}}{k_B T_e} \right)^m \quad (10) $$
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Sketch of a sheath layer forming on a material facing the plasma.</figcaption></figure>

with $V_{DC}$ the DC potential and $m$ a parameter that typically changes between 2/3 and 3/4 depending on the incidence angle of magnetic field lines on materials and the collisional properties of the sheath as explored through Particle in Cell simulations in [48]. In general, the electric field accelerates the ions towards the wall, leading to an increase of the sputtering and the production of impurities (Figure 4) in proportions usually comparable to the increase in electron temperature:

$$V_{\text{Thermal sheath}} = 3k_{B}T_{e}/2 \qquad (11)$$

On the other hand, in the presence of ICRF heating, this phenomenon is further aggravated by the so-called rectification of the oscillating potentials associated with the slow wave carrying electric field in the direction parallel to the magnetic field lines (cf. red curve in Figure 3). Let us now express the current instantly conducted to any object exposed to the plasma as:

$$\tilde{I}(t) = I_{\text{sat}}^{+} \left( 1 - \exp \left( \frac{e}{k_{B}T_{e}} (\phi_{f} - \tilde{V}(t)) \right) \right) \qquad (12)$$

with $I_{sat}^{+}$ the ion saturated current, $\phi_{f}$ the floating potential and $\tilde{V}(t) = V_{DC} + V_{RF} \cos(\omega_{0}t)$ the sheath potential in presence of an RF wave.

Then the RF → DC rectification, specific to waves at the ion cyclotron frequency, results from the fact that the electrons (light compared to the ions) react instantaneously to the oscillating electric field of the wave ($V_{RF}$), while the ions react only to the average value.

$$\langle \tilde{I}(t) \rangle = \bar{I} = I_{\text{sat}}^{+} \left[ 1 - \exp \left( \frac{e}{k_{B}T_{e}} (\phi_{f} - V_{\text{DC}}) \right) I_{0} \left( \frac{eV_{\text{RF}}}{k_{B}T_{e}} \right) \right] \qquad (13)$$
with $I_0$ a Bessel function. The consequence at the level of the components facing the plasma on the scale of an RF period, is the appearance of a DC current due to the privileged drift of the electrons. Thus, to compensate this current, the DC potential of the sheath must adjust and get biased by an extra potential $V_b$:

$$V_b = \frac{k_B T_e}{e} I_0 \left( \frac{eV_{RF}}{k_B T_e} \right) \qquad (14)$$

The RF → DC rectification of the sheath therefore has the effect of increasing the potential drop (electric field) through the skin between the plasma and the wall, i.e., the acceleration of the ions towards the materials. This acceleration often boosts ions incident energies up to critical levels for exposed components (i.e., above sputtering yields threshold values). For deeper understanding of the topic, the reader is highly encouraged to refer to excellent tutorial [49] and reviews of key experimental [26–50] and modelling [51] results.

## 3.4 Physical sputtering yield

The other key process by which ICRF operation often enhances impurity production is the physical sputtering yield $Y_{Eff}$. This parameter represents the probability for an ion to sputter an atom from a target, given its charge, its incident energy, and their respective mass:

$$\text{Impurity outflux} = \Gamma_{\text{target atoms out}} = \Gamma_{\text{ion in } Y_{\text{Eff}}} \qquad (15)$$

For example, the graph in **Figure 5** shows how the sputtering of a graphite (~carbon) target evolves for different bombarding species. The first important aspect of the sputtering yield is that it varies non-linearly with the energy of the incident ion. The second one is that only above a threshold energy of about 30eV, do the incident species start to sputter carbon atoms. The red curve also represents the self-sputtering of lithium which is sometimes used for material coating to prevent from important plasma contamination by high-Z impurities. One can see that a lithium coating will typically get eroded at much lower energy, while plasma facing components in graphite will require larger energies. We have here taken the case of carbon components, sometimes coated with lithium, exposed to various species, which is the case in EAST tokamak. Carbon however won't be used in future devices plasma facing components due to tritium retention. Tungsten has now become the most suitable candidate. In WEST for instance, antenna limiters are coated with tungsten, which sputtering follows similar trends (cf. **Figure 3** of [23]).

## 3.5 Poor coupling efficiency and near-field effects

RF-sheath excitation is the phenomenon proper to ICRF which is responsible for its peculiar complexity, due to the non-linear trade-off relation between the maximization of the coupling and the minimization of the interactions with the plasma. For any object in principle, moving away from the plasma results in a strict decrease of their interactions. However, for a classic ICRF antenna, moving away first induces a drop in its coupling efficiency, therefore an increase of excited fields (consistent with an increase of the VSWR in the transmission lines), which in turn can in critical cases result in a local but exponential increase of the RF sheath potential (cf. Figure 3 of [52]
<figure><img src="image_6.png" /><figcaption>**Figure 5.**
Sputtering yield of various ion species on a target in graphite.</figcaption></figure>

and section 3.3 of this chapter). In addition, the sputtering yield not only changes non-linearly, but also reaches its largest values for energies typically in the range of potential rectified by RF sheaths (cf. red region between 100eV and 1keV in **Figure 4**). It follows that in some cases, increasing the distance between an ICRF antenna and the plasma to reduce their interactions based on linear intuition, can unexpectedly cause an increase of impurity production from local sources, reacting to sharp increase of RF sheath potentials. These phenomena can be modeled in different ways as summarized in [53]. This behavior was also experimentally observed in WEST during a scan of the radial position of an ICRF antenna, where the floating potential measured by a reciprocating probe did not directly decrease as the antenna was retracted but first passed by a maximum for intermediate position (cf. Figure 4 in [50] or Figure 6.23 of [22] in open access). A wide variety of RF-sheath induced plasma surface interactions have also been observed in devices such as:

*   EAST with the partial melting of metallic plates at the corners of an ICRF antenna, which precisely did not melt at the closest point to the plasma, but at the locations where simulations predict RF sheath excitation to be maximal (cf. Figure 14 in [54] or Figure 6.27 of [22] in open access)

*   Tore Supra with large, enhanced potentials measured with retarded field analyzer [55] causing large heat loads on antennas [56]
* A LineAr Plasma Device (LAPD) using diagnostic magnetically connected to ICRF antenna along large number of reproducible discharges to provide great evidence of ICRF interaction process in the edge plasma [57]

* Alcator C-Mod with strong potentials measured with several different probes [58] caused by RF sheath excitation [59], leading to the increase of impurity sources and deleterious influence on plasma performance [60]

* ASDEX Upgrade with similar experiments [61] that have inspired a series of successful technological upgrades [62, 63] discussed hereafter in section 4

* JET is also subject to these effects [26, 64]

Interactions occurring at vicinity or in regions magnetically connected to the active antenna frame (red start with most branches around the antenna in **Figure 2**) are often classified as *near field effects*. Hence, interactions respectively taking place elsewhere will be conveniently referred to as *far field effects*<fnref n="1" />.

## 3.6 Abnormal wave absorption and far-field effects

Following the same logic, we will categorize interactions occurring in regions without magnetic connection to the active antenna as *far field effects*. These can be observed in unexpected locations when abnormal propagation and absorption of the wave take place. For instance, if the absorption efficiency is low, the power launched in the plasma is not fully absorbed at the first pass, and a non-negligible fraction of unabsorbed power propagates and reflects on farther objects such as the divertor or the inner wall as represented in **Figure 2**. In these cases, not only heating performance drop, but RF-sheath can also potentially become excited in global fashion. Such effects were well-observed in

* EAST *high field side wall* (cf. red star with 4 branches in **Figure 2**) when the antenna operated in monopole phasing, causing a drop in the absorption efficiency and an increase of fields excitation at the inner wall [54]

* WEST *divertor* if too low hydrogen ratio leads to poor wave absorption [65]

* NSTX *divertor* where operating at high harmonics led to the excitation of surface waves [66] and substantial fraction of their power to be coupled to the edge filamentary structures parallel to the magnetic field lines [67], resulting in potential rectification at their extremity and deleterious interactions [68].

<footnote n="1">¹ Here, *near* and *far field effects* should not be confused with the near and far electromagnetic fields diffusion patterns around an object, but simply refer to the location where the effect occurs with respect to the ICRF antenna.</footnote>
# 4. The beauty of ICRH: assets and technological perspectives

## 4.1 Assets

Despite the many challenges, ICRF has remained a very important tool for heating plasmas, with constant progress over the past fifty years [69]. ICRF is the only heating system that can directly deposit power onto the ions, allowing to generate fast ions and significantly boost fusion performance. While minority heating has long remained the most widely used method, this past two decades have been marked by the emergence of several modelling tools allowing to predict wave absorption efficiency [70] and power thermalization on the various species [71, 72]. These tools have allowed to explore much wider varieties of heating schemes [29]. Most promising ones were then tried experimentally, and their efficiency confirmed. One can typically quote scenarios based on the synergy with NBI [31, 32, 73] and these relying on three-ion scheme [74]. In particular, the three-ion scheme has gained increasing interest not only due to the efficiency in generating fast ions, but also for the operation flexibility it offers. Indeed, instead of relying on specific fraction of a minority specie—which can be difficult to maintain over operation—it proves much easier to keep control of rough proportions of two (or more) majority species and deposit power on a third very minor one which fraction is easier to control. Most promising scenario for fusion reactor D-T plasmas would therefore consist in depositing power most likely on ³He, which has been successfully achieved in JET and Alcator C-Mod [30] and foreseen for future devices like ITER [74] and SPARC [75].

Another application of ICRH is the control of impurity population in the core. It has for instance been observed in JET, that for ICRF power above 4MW, heating H minority can result in high-Z impurity screening out of the core [76]. Furthermore, the three-ion scheme even allows to directly pump out an impurity by depositing part of the power on it. This was observed with beryllium in JET [29], with Argon on TFR [77] and Alcator C-Mod [78]. The same trick could further allow to pump out tungsten impurity by targeting the second harmonic of W⁵⁶⁺ in SPARC's high temperature plasmas [78], a powerful tool to prevent tungsten accumulation in the core.

## 4.2 Technological perspectives

A few research and development ideas for improving current design of antennas yet remain to be tested. One can quote:

* *Travelling wave antennas* [43] could improve the reliability of ICRF systems thanks to its inherent load resilience and much larger coupling capacities, offering the possibility to place the antenna 10cm further away from the plasma than a classic antenna and still be able to couple similar amount of power. This could have tremendous impact not only on the coupling efficiency but also on the interactions with the plasma and its contamination by metallic impurity. The disadvantage of this concept is its low power density due to the large number of straps (>5). This makes any design hardly compatible with fusion reactor worried of maximizing the surface of the tritium breeding blanket. But maybe some tricks have not yet been thought about ...
* *Active limiters* could allow compensating actively the image currents induced on limiters by nearby current straps. Relying on the concept of proximity effect [27], these currents can be cancelled by powering straps with appropriate phasing and power ratio, as was shown in ASDEX-Upgrade 3-strap antennas [63] and Alcator C-Mod 4-strap antennas [79]. However, operating the ICRF system with optimal settings is not only challenging in terms of control, but mostly limiting in terms of power density. *Active limiters* would then consist in applying a fraction of the power sent on the straps, directly on the limiters. This could allow reducing the local impurity sources, by canceling off the local currents that tend to be responsible for RF sheath excitation. In addition, such solution would offer two key advantages:

○ *more flexibility* during ICRF operation by opening the possibility to routinely tune the antenna *phasing* without aggravating interactions with the plasma. In other words, this would allow operating with different antenna spectra to cope with plasma changes regardless of near-field effects

<figure><img src="image_7.png" /><figcaption>**Figure 6.** Summary of the key components of an ICRF system, and a selection of the main phenomena that can occur in various location. Black boxes enlighten phenomena inherent to ICRF powering. Blue boxes show the physical parameters that are most directly influenced. Red boxes show the potentially deleterious consequences of parameters changes. Green boxes with green and orange borders show the main solutions respectively existing and to be explored.</figcaption></figure>
* *maximize power density*. While this is already quite important in nowadays experimental tokamaks, it will be a critical aspect on future fusion reactors, which aim at maximizing the surface of the tritium breeding blanket for producing fuel and extracting the energy radiated by the plasma to generate electricity. Simultaneously satisfying requirements on plasma heating, fuel production and energy exhaust, pushes towards more compact designs with larger power density currently hardly compatible with low level of impurity production. Active limiters could therefore be a major technological solution.

**5. Conclusion**

Despite all challenges, heating fusion plasmas with waves in the ion cyclotron
range of frequency is still very attractive. While several technological solutions remain
to be explored, a variety of existing ones are already routinely used and allow suc-
cessful operation in many devices around the globe. A summary of the key compo-
nents of an ICRF system, and a selection of the main phenomena that can occur in
various location is provided in **Figure 6**. To conclude on the perspective of fusion
reactors, while ICRF may not be the optimal solution in terms of power density, note
that as the intensity of the magnetic field due to progress made on low temperature
superconductors and the size of the devices will increase, ICRF is the only heating
system which difficulty to implement will not be affected.

For further reading of on this topic, the reader is highly encourage to refer to a
review of the “Recent Progress in ICRF in Magnetic Confinement Fusion” [80] and
references therein.

Acknowledgements

Laurent Colas, Julien Hillairet, Karl Vulliez, Annika Ekedahl, and more globally WEST and EAST ICRF Teams are warmly acknowledged for their efforts in making experiments possible and the numerous fruitful discussions. I also want to thank a lot Wouter Tierens, Mari Usoltceva, Vladimir Bobkov, and the rest of ASDEX ICRF team, Ernesto Lerche from ERM and colleagues from the 107 team of Institut Jean Lamour for all the things we collaborated over the past few years.

Author details

Guillaume Urbanczyk
Institute of Plasma Physics, Chinese Academy of Sciences, Hefei, China

*Address all correspondence to: guiguiurban@hotmail.com

IntechOpen

© 2022 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

<figure><img src="image_8.png" /></figure>
**References**

[1] Faugel H, Bobkov V, Fünfgelder H. The prospects of ICRF generators at ASDEX upgrade. AIP Conference Proceedings. 2020;2254:070006. DOI: 10.1063/5.0013529

[2] Bourdelle C et al. WEST Physics Basis. Nuclear Fusion. 2015;55:063017. DOI: 10.1088/0029-5515/55/6/063017

[3] Vulliez K, Chen Z, Ferlay F, Winkler K, Helou W, Hillairet J. The mechanical structure of the WEST ion cyclotron resonant heating launchers. Fusion Engineering and Design. 2015; 96-97:611-615. DOI: 10.1016/j.fusengdes.2015.06.050

[4] Wade TJ, Jacquinot J, Bosia G, Sibley A, Schmid M, et al. Development of the JET ICRH plant. Fusion Engineering and Design. 1994;24:23-46. DOI: 10.1016/0920-3796(94)90035-3

[5] Durodié F, Dumortier P, Blackman T, Wooldridge E, Lerche E, Helou W, et al. ITER-like antenna for JET first results of the advanced matching control algorithms. Fusion Engineering and Design. 2017;123 (November):253-258. DOI: 10.1016/j.fusengdes.2017.05.113

[6] Monakhov I, Graham M, Blackman T, Dowson S, Durodie F, Jacquet P, et al. Design and operations of a load-tolerant external conjugate-T matching system for the A2 ICRH antennas at JET. Nuclear Fusion. 2013;53(8):083013. DOI: 10.1088/0029-5515/53/8/083013

[7] Lin Y et al. ICRF antenna matching systems with ferrite tuners for the Alcator C-Mod tokamak. Fusion Engineering and Design. 2015;100: 239-248. DOI: 10.1016/j.fusengdes. 2015.06.022

[8] Swain DW et al. Loading and asymmetry measurements and modeling

for the national spherical torus experiment ion cyclotron range of frequencies system. Fusion Science and Technology. 2003;43:503-513. DOI: 10.13182/FST03-A297

[9] Durodié F, Vervier M. First results of the automatic matching device for textor's ICRH system, Fusion Technology. 1992. pp. 477-480. DOI: 10.1016/B978-0-444-89995-8.50088-6. ISBN 9780444899958

[10] Hangan D-M et al. Analysis of dynamic matching networks for the ICRF system at ASDEX upgrade. Fusion Engineering and Design. 2011;86: 736-741. DOI: 10.1016/j.fusengdes. 2011.02.001

[11] Saito K et al. Real-time impedance matching system for ICRF heating in LHD. Fusion Engineering and Design. 2008;83:245-248. DOI: 10.1016/j.fusengdes.2008.01.017

[12] Zhao Y-P et al. EAST ion cyclotron resonance heating system for long pulse operation. Fusion Engineering and Design. 2014;89:2642-2646. DOI: 10.1016/j.fusengdes.2014.06.017

[13] Braun F et al. ICRF system enhancements at ASDEX upgrade. Fusion Engineering and Design. 2001; 56-57:551-555. DOI: 10.1016/S0920-3796 (01)00274-5

[14] Lamalle PU et al. Commissioning of the wideband matching system for ICRH of ELMy JET Plasmas. AIP Conference Proceedings. 2001;595:118

[15] Parisot A. Design of an ICRF fast matching system for the Alcator C-Mod tokamak. PhD Thesis. 2004. Available from: https://dspace.mit.edu/bitstream/handle/1721.1/93763/
04RR002_full.pdf?sequence=1&
isAllowed=y

[16] Durodié F et al. Physics and engineering results obtained with the ion cyclotron range of frequencies ITER-like antenna on JET. Plasma Physics and Controlled Fusion. 2012;54:074012. DOI: 10.1088/0741-3335/54/7/074012

[17] Helou W, Colas L, Hillairet J, Milanesio D, Mollard P, Argouarch A, et al. Radio-frequency electrical design of the WEST long pulse and load-resilient ICRH launchers. Fusion Engineering and Design. 2015;96-97 (October):473-476. DOI: 10.1016/j.fusengdes.2015.01.005

[18] Helou W. Design and operation of antennas at the ion cyclotron and lower hybrid range of frequencies for nuclear fusion reactors. PhD Thesis. 2017. Available from: http://www.theses.fr/2017AIXM0395

[19] Goulding RH, Hoffman DJ, Ryan PM. Power compensators for phased operation of antenna arrays on JET and DIII-D. AIP Conference Proceedings. 1994;289:351. DOI: 10.1063/1.44958

[20] Gen C, Mao YZ, Zhao YP, Yuan S, Zhang XJ, Qing CM. Development and test of decoupler for ICRF antenna in EAST. Fusion Engineering and Design. 2016;107:32-37. DOI: 10.1016/j.fusengdes.2016.04.010

[21] Grine D et al. Summary and results of the study of the hybrid matching option implementation of the ITER ICRH system. Fusion Engineering and Design. 2012;87:167-178. DOI: 10.1016/j.fusengdes.2011.12.006

[22] Urbanczyk G. Interaction of High-Power waves with the plasma periphery of WEST & EAST tokamaks. PhD

Thesis. 2019. Available from: https://hal.univ-lorraine.fr/tel-02517768/document

[23] Urbanczyk G, Colas L, Hillairet J, Lerche E, Fedorczak N, et al. RF wave coupling, plasma heating and characterization of induced plasma-material interactions in WEST L-mode discharges. Nuclear Fusion. 2021;61: 086027. DOI: 10.1088/1741-4326/ac0d11

[24] Urbanczyk G, Zhang XZ, Colas L, et al. Optimization of discharges with ion cyclotron range of frequencies using local gas injection in EAST. Nuclear Fusion. 2019;59:066023. DOI: 10.1088/1741-4326/ab1610

[25] Zhang W, Jacquet P, Lerche E, Bilato R, Bobkov V, et al. 3D simulations of gas puff effects on edge plasma and ICRF coupling in JET. Nuclear Fusion. 2017;57:056042. DOI: 10.1088/1741-4326/aa6817

[26] Bobkov V, Aguiam D, Bilato R, Brezinsek S, Colas L, et al. Impact of ICRF on the scrape-off layer and on plasma wall interactions: From present experiments to fusion reactor. Nuclear Materials and Energy. 2019;18:131-140. DOI: 10.1016/j.nme.2018.11.017

[27] Colas L, Lu L-F, Křivská A, et al. Spatial proximity effects on the excitation of Sheath RF voltages by evanescent slow waves in the ion cyclotron range of frequencies. Plasma Physics and Controlled Fusion. 2017;59: 025014. DOI: 10.1088/1361-6587/59/2/025014

[28] Lerche E, Bobkov V, Jacquet P, Mayoral M-L, Messiaen A, et al. Recent experiments on alternative dipole phasing with the JET A2 ICRF antennas. AIP Conference Proceedings. 2009;1187: 93. DOI: 10.1063/1.3273845
[29] Lerche E, Van Eester D, Jacquet P, Casson F, et al. ICRH options for JET-ILW DTE2 operation. AIP Conference Proceedings. 2020;2254:030007. DOI: 10.1063/5.0013530

[30] Kazakov Y, Ongena J, Wright J, et al. Efficient generation of energetic ions in multi-ion plasmas by radio-frequency heating. Nature Physics. 2017;13: 973-978. DOI: 10.1038/nphys4167

[31] Kirov KK, Kazakov Y, Nocente M, Ongena J, et al. Synergistic ICRH and NBI heating for fast ion generation and maximising fusion rate in mixed plasmas at JET. AIP Conference Proceedings. 2020;2254:030011. DOI: 10.1063/5.0014235

[32] Huynh P, Lerche EA, Van Eester D, Bilato R, et al. Modeling ICRH and ICRH-NBI synergy in high power JET scenarios using European transport simulator (ETS). AIP Conference Proceedings. 2020;2254:060003. DOI: 10.1063/5.0014240

[33] Keilhacker M et al. Confinement studies in L and H-type Asdex discharges. Plasma Physics and Controlled Fusion. 1984;26(1A):49-63. DOI: 10.1088/0741-3335/26/1A/305

[34] Bobkov V. Studies of high voltage breakdown phenomena on ICRF antennas. 2003. Available from: http://mediatum.ub.tum.de/doc/602957/document.pdf

[35] Bobkov V, Braun F, Hartmann DA, et al. Influence of ELMs on operation of ICRF antennas in ASDEX Upgrade. Journal of Nuclear Materials. 2005; 337-339:776-780. DOI: 10.1016/j.jnucmat.2004.10.109

[36] Ladurelle L, Agarici G, Beaumont B, et al. First results of automatic matching system on Tore Supra ICRH antennas:

fast matching network for ICRH systems. In: Proceedings of 19th SOFT. Lisbonne, Portugal; 1996. DOI: 10.1016/B978-0-444-82762-3.50118-X

[37] Durodié F, Dumortier P, Helou W, et al. Circuit model of the ITER-like antenna for JET and simulation of its control algorithms. AIP Conference Proceedings. 2015;1689(070013). DOI: 10.1063/1.4936520

[38] Noterdaeme J-M, Bobkov V, Brémond S, Parisot A, et al. Matching to ELMy plasmas in the ICRF domain. Fusion Engineering and Design. 2005;74: 191-198. DOI: 10.1016/j.fusengdes.2005.06.071

[39] Faugel H, Angene P, Becker W, Braun F, Bobkov V, et al. The ASDEX upgrade ICRF system: Operational experience and developments. Fusion Engineering and Design. 2005;74: 319-324. DOI: 10.1016/j.fusengdes.2005.06.268

[40] Yang H, Zhang XJ, Yuan S, Qin CM, et al. Experimental validation of the newly designed ICRF antenna on EAST. Private Communication

[41] Vulliez K, Argouarch A, Bosia G, et al. Validation of the load-resilient ion cyclotron resonance frequency antenna concept on Tore Supra plasmas. Nuclear Fusion. 2008;48:065007. DOI: 10.1088/0029-5515/48/6/065007

[42] Hillairet J, Mollard P, Colas L, Helou W, Urbanczyk G, et al. WEST actively cooled load resilient ion cyclotron resonance heating system results. Nuclear Fusion. 2021;61:096030. DOI: 10.1088/1741-4326/ac1759

[43] Ragona R, Messiaen A. Conceptual study of an ICRH traveling-wave antenna system for low-coupling conditions as expected in DEMO.
Nuclear Fusion. 2016;56:076009.
DOI: 10.1088/0029-5515/56/7/076009

[44] Pinsker RI, Prater R, Moeller CP, deGrassie JS, et al. Experiments on helicons in DIII-D-investigation of the physics of a reactor-relevant non-inductive current drive technology. Nuclear Fusion. 2018;58:106007.
DOI: 10.1088/1741-4326/aad1f8

[45] Van Compernolle B, Brookman MW, Moeller CP, Pinsker RI, Garofalo AM, O'Neill R, et al. The high-power helicon program at DIII-D: Gearing up for first experiments. Nuclear Fusion. 2021; 61(11):116034. DOI: 10.1088/1741-4326/ac25c0

[46] Ragona R, Batal T, Hillairet J, Messiaen A, et al. Progress on the design of a DEMO high power ICRH travelling wave antenna mock-up to be tested on WEST. AIP Conference Proceedings. 2020;2254:070014. DOI: 10.1063/5.0014215

[47] Stangeby PC. The Plasma Boundary of Magnetic Fusion Devices (1st ed.). Boca Raton: CRC Press. 2000.
DOI: 10.1201/9780367801489

[48] Moritz J, Heuraux S, Gravier E, Lesur M, Brochard F, et al. Sheath size and Child-Langmuir law in one dimensional bounded plasma system in the presence of an oblique magnetic field: PIC results. Physics of Plasmas. 2021;28:083501. DOI: 10.1063/5.0055790

[49] Myra JR. A tutorial on radio frequency sheath physics for magnetically confined fusion devices. Journal of Plasma Physics. 2021;87(5): 905870504. DOI: 10.1017/S0022377821000878

[50] Colas L, Urbanczyk G, Goniche M, Hillairet J, Bernard J-M, Bourdelle C, et al. The geometry of the ICRF-induced

wave-SOL interaction. A multi-machine experimental review in view of the ITER operation. Nuclear Fusion. 2022; 62(1):016014. DOI: 10.1088/1741-4326/ac35f9

[51] Colas L, Lu LF, Jacquot J, Tierens W, et al. SOL RF physics modelling in Europe, in support of ICRF experiments. EPJ Web of Conferences. 2017;157: 01001. DOI: 10.1051/epjconf/201715701001

[52] Myra JR, D'Ippolito DA. Resonance cone interaction with a self-consistent radio-frequency sheath. Physical Review Letters. 2008;101:195004. DOI: 10.1103/PhysRevLett.101.195004

[53] Zhang W, Bilato R, Bobkov VV, et al. Recent progress in modeling ICRF-edge plasma interactions with application to ASDEX upgrade. Nuclear Fusion. 2021. DOI: 10.1088/1741-4326/ac38c8

[54] Urbanczyk G, Zhang XJ, Colas L, Dumont R, Tierens W, et al. Metallic impurity content behavior during ICRH-heated L-mode discharges in EAST. Nuclear Fusion. 2020;60:126003. DOI: 10.1088/1741-4326/abae82

[55] Martin K, Gunn JP, Heuraux S, Colas L, Faudot E, Jacquot J. Measurement of sheath potential in RF-biased flux tubes using a retarding field analyzer in Tore Supra tokamak. Journal of Nuclear Materials. 2013;438: S509-S512. DOI: 10.1016/j.jnucmat.2013.01.105

[56] Corre Y, Firdaouss M, Colas L, Argouarch A, Guilhem D, Gunn J, et al. Characterization of heat flux generated by ICRH heating with cantilevered bars and a slotted box Faraday screen. Nuclear Fusion. 2012;52:103010. DOI: 10.1088/0029-5515/52/10/103010
[57] Martin MJ, Gekelman W, Van Compernolle B, Pribyl P, Carter T. Experimental observation of convective cell formation due to a fast wave antenna in the large plasma device. Physical Review Letters. 2017;119:205002. DOI: 10.1103/PhysRevLett.119.205002

[58] Ochoukov R, Whyte DG, Brunner D, D'Ippolito DA, et al. ICRF-enhanced plasma potentials in the SOL of Alcator C-Mod. Plasma Physics and Controlled Fusion. 2014;56:015004. DOI: 10.1088/0741-3335/56/1/015004

[59] Wukitch SJ, LaBombard B, Lin Y, Lipschultz B, Marmar E, et al. ICRF specific impurity sources and plasma sheaths in Alcator C-Mod. Journal of Nuclear Materials. 2009;390-391: 951-954. DOI: 10.1016/j.jnucmat.2009.01.245

[60] Wukitch SJ, Lipschultz B, Marmar E, Lin Y, et al. RF plasma edge interactions and their impact on ICRF antenna performance in Alcator C-Mod. Journal of Nuclear Materials. 2007; 363-365:491-497. DOI: 10.1016/j.jnucmat.2007.01.273

[61] Dux R, Bobkov V, Herrmann A, Janzer A, Kallenbach A, et al. Plasma-wall interaction and plasma behaviour in the non-boronised all tungsten ASDEX Upgrade. Journal of Nuclear Materials. 2009;390-391:858-863. DOI: 10.1016/j.jnucmat.2009.01.225

[62] Bobkov V, Balden M, Bilato R, Braun F, et al. ICRF operation with improved antennas in ASDEX Upgrade with W wall. Nuclear Fusion. 2013;53: 093018. DOI: 10.1088/0029-5515/53/9/093018

[63] Bobkov V, Aguiam D, Bilato R, Brezinsek S, Colas L, et al. Making ICRF power compatible with a high-Z wall in ASDEX upgrade. Plasma Physics and

Controlled Fusion. 2017;59:014022. DOI: 10.1088/0741-3335/59/1/014022

[64] Colas L, Jacquet P, Bobkov V, Brix M, Meneses L, et al. 2D mappings of ICRF-induced SOL density modifications on JET. In: 45th EPS Conference on Plasma Physics Proceedings. http://ocs.ciemat.es/ EPS2018PAP/pdf/O4.101.pdf

[65] Urbanczyk G, Colas L, Zhang XJ, Helou W, Hillairet J, et al. ICRH coupling optimization and impurity behavior in EAST and WEST. AIP Conference Proceedings. 2020;2254:030012. DOI: 10.1063/5.0018453

[66] Tierens W, Colas L, EUROfusion MST1 Team. Slab-geometry surface waves on steep gradients and the origin of related numerical issues in a variety of ICRF codes. Journal of Plasma Physics. 2021;87(4):905870405. DOI: 10.1017/S002237782100074X

[67] Tierens W, Zhang W, Manz P, et al. The importance of realistic plasma filament waveforms for the study of resonant wave-filament interactions in tokamak edge plasmas editors-pick Physics of Plasmas. 2020;27:052102. DOI: 10.1063/5.0007098

[68] Perkins RJ, Hosea JC, et al. High-harmonic fast-wave power flow along magnetic field lines in the scrape-off layer of NSTX. Physical Review Letters. 2012;109:045001. DOI: 10.1103/PhysRevLett.109.045001

[69] Noterdaeme JM. Fifty years of progress in ICRF, from first experiments on the model C stellarator to the design of an ICRF system for DEMO. AIP Conference Proceedings. 2020;2254: 020001. DOI: 10.1063/5.0014254

[70] Van Eester D, Lerche E. Integro-differential modeling of ICRH wave
propagation and damping at arbitrary
cyclotron harmonics and wavelengths in
tokamaks. AIP Conference Proceedings.
2014;1580:298. DOI: 10.1063/1.4864547

[71] Dumont RJ. Variational approach to radiofrequency waves in magnetic fusion devices. Nuclear Fusion. 2009;49:075033. DOI: 10.1088/0029-5515/49/7/075033

[72] Brambilla M. Numerical simulation of ion cyclotron waves in tokamak plasmas. Plasma Physics and Controlled Fusion. 1999;41:1. DOI: 10.1088/0741-3335/41/1/002

[73] Bilato R, Volpe F, Poli E, Köhn A, et al. Electron Cyclotron Heating in RFP plasmas. AIP Conference Proceedings. 2009;1187:81. DOI: 10.1063/1.3273812

[74] Kazakov YO, Ongena J, Wright JC, et al. Physics and applications of three-ion ICRF scenarios for fusion research. Physics of Plasmas. 2021;28:020501. DOI: 10.1063/5.0021818

[75] Lin Y, Wright J, Wukitch S. Physics basis for the ICRF system of the SPARC tokamak. Journal of Plasma Physics. 2020;86(5):865860506. DOI: 10.1017/S0022377820001269

[76] Lerche E, Goniche M, Jacquet P, Van Eester D, Bobkov V, Colas L, et al. Optimization of ICRH for core impurity control in JET-ILW. Nuclear Fusion. 2016;56:036022. DOI: 10.1088/0029-5515/56/3/036022

[77] TFR Group. Impurity pump-out at the two-ion hybrid resonance during ICRF experiments in TFR tokamak plasmas. Nuclear Fusion. 1982;22:956

[78] Rice JE, Lin Y, Perks CJ, Reinke ML, et al. Argon pumpout by ICRF waves in C-Mod L- and I-mode plasmas. Nuclear Fusion. 62 086009; 2022. DOI: 10.1088/1741-4326/ac6ef0

[79] Wukitch S et al. Towards ICRF antennas compatible with high performance plasmas: characterization and mitigation of ICRF antenna – plasma edge interaction. In: 22nd RFPPC Conference. France: Aix-en-Provence; 2017

[80] Hillairet J. Recent progress in ion cyclotron range of frequency in magnetic confinement fusion. From modelling technology and experiments. Reviews of Modern Physics. 2022. (Under reviewing process)