The neoclassical tearing mode (NTM) is expected to form a magnetic island structure which locally flattens the plasma pressure, hence losing the bootstrap current    . Low mode number NTMs, especially the (m and n are the poloidal and toroidal mode numbers)and 2/1 modes, are observed to cause the confinement degradation and/or disruptions, limiting the achievable plasma β value well below the predictions of ideal magnetohydrodynamic (MHD) calculations for positive magnetic shear -.
In general, NTMs can be stabilized by compensating the perturbed bootstrap current or tailoring the plasma current density profile using local radio frequency (rf) current. The driven current should be in the direction along the equilibrium plasma current, to compensate the missing bootstrap current inside the island       . The stabilization of NTMs by rf current drive has been successfully carried out in tokamak experiments, e.g. on ASDEX-U, DIII-D, JT-60U     . The modified Rutherford equation (MRE), governing the time evolution of the nonlinear island, has been calculated by several authors    . However, the predictions of the required rf current for mode stabilization are still in question and this is the purpose of this paper.
2. Island Evolution Equation
The MRE which describes the temporal evolution of the full width of a magnetic island
in which , is the rational surface and is the plasma resistivity.
is the tearing mode stability index, defined as the logarithmic discontinuity in the radial derivative of the perturbed magnetic flux function .In the vicinity of the rational surface the well-known “constant approximation” is utilized. The bootstrap current perturbation term can be written as
where is the poloidal magnetic field, the plasma pressure,
the magnetic shear length at the rational surface,
the safety factor, and
Equation (1) is obtained through averaging of the parallel projection of the Ohm’s law over helical flux surfaces in the vicinity of the island
It is convenient to introduce the normalized perturbed flux function   
with at the island O-point and at the X-point, where and is the helical angle. Following Ref. , the flux surface average operator is defined as
for , and
for . Here, , , and is a general function.
In this way, the driven current term is given by
where is the magnitude of the rf driven current and is the half width of the driven current. With a Gaussian profile for the driven current by rf waves  
the stabilization efficiency is given by
where is the peak amplitude of the driven current and the square box function is equal to 1 for and to 0 otherwise, defining the helical angle of the rf current deposition, where , is the mode rotation frequency, and is the instantaneous wave deposition width along the helical angle.
In this paper, we neglect the modification of the equilibrium plasma current profile by the current drive and the reduced island width, which can modify the .
3. Required rf Current for Mode Stabilization
In this paper the inverse aspect ratio is utilized, where is the plasma minor radius and the major radius. Results obtained from the MRE utilizing Equation (1) for the magnetic island and its stabilization with rf current are shown in Figure 1, where , at and are taken. is assumed. The solid curve is the case without the driven current term. The dashed and the dot-dashed curves are for modulated current drive (MCD, 50% on-time) with and for non-modulated current drive (NMCD) with , respectively. For these values of the maximal of is zero, i.e., the required driven current for mode stabilization.
The required for mode stabilization is shown as a function of for in Figure 2. The solid (dashed) curve is for MCD (NMCD). It is seen the required depends linearly on for both MCD and NMCD.
Figure 1. Evaluation of from MRE for the magnetic island with , at and . is utilized. The solid curve is the case without the driven current. The dashed and dot-dashed curves are for MCD with and for NMCD with , respectively.
Figure 2. Required for mode stabilization versus for . The solid (dashed) curve is for MCD (NMCD).
Figure 3 shows the ratio as a function of the half width of driven current for . There is an approximate linear and a quadratic relationship for MCD and NMCD, respectively. According to Equations (1) and
(7) the required driven current is proportional to y, where . For
different island width there is always a linear (quadratic) relationship between y and with MCD (NMCD).
For above results the localized current drive is applied right at the resonant surface. The radial misalignment between the rf waves deposition and the resonant surface significantly affects the NTM stabilization. For and , the required for mode stabilization is shown in Figure 4 as a function of the normalized radial deviation of rf current deposition from the
rational surface, , for MCD (solid curve) and NMCD (dashed). It is
Figure 3. Required for mode stabilization versus for . The solid (dashed) curve is for MCD (NMCD).
Figure 4. Required for mode stabilization versus the radial deviation of the rf current , where . The solid (dashed) curve is for MCD (NMCD).
seen between the ratio and the radial deviation xcd there is always an approximate quadratic relationship for both MCD and NMCD. A shift of rf current deposition away from the rational surface greatly decreases the stabilizing efficiency, since the applied driven current is not around the island’s o-point. The required driven current for NMCD is much more than that for MCD when . The reason is that NMCD deposited around the island separatrix and the island X-point is destabilizing.
The magnetic island width changes with the pressure profile, which in turn affects the island growth. The modification of the equilibrium plasma current profile by the current drive can affect the . These are not considered in our calculations. And the numerical simulation on the required driven current for mode stabilization which consistently calculates the pressure profile and the plasma current density profile will be the future work.
In summary, the required driven current for NTMs’ stabilization by rf current is investigated by the MRE. The simple linear and/or quadratic relations are found between the required rf current for mode stabilization and the bootstrap current density, the radial width and radial deviation of the rf current from the rational surface for both MCD and NMCD.
This work was supported by the National Magnetic Confinement Fusion Science Program of China (Contract No. 2012GB103000 and 2012GB103003) and the National Natural Science Foundation of China (Grant No. 11475225).