Infernal Modes in Reversed Shear

• Infernal modes are pressure-driven MHD instabilities in tokamak plasmas with low or reversed magnetic shear, where weakened field-line bending permits rapid mode growth.
• The topic explains how non-monotonic safety-factor profiles and toroidal sideband coupling, influenced by shaping effects like elongation and triangularity, modulate mode structure and dispersion relations.
• Experimental and simulation data, such as from the CFETR 1GW SSO equilibrium, highlight strategies to optimize q-profiles and control rotation for managing infernal instabilities in advanced tokamaks.

Infernal modes in reversed shear refer to a class of pressure-driven long-wavelength MHD instabilities that manifest in toroidal fusion plasmas with extended regions of very low or reversed magnetic shear. These modes are closely associated with non-monotonic safety-factor ($q$) profiles, particularly those exhibiting flat or locally reversed shear, such that $q'(r) \approx 0$ over a finite interval. In these configurations, the stabilizing field-line bending effect is significantly weakened, allowing the destabilizing influence of pressure gradients and toroidal coupling to dominate. Infernal modes are characterized by broad radial structures spanning the low-shear region, discrete spectra of fast-growing modes for given $(m,n)$, and, under sufficient resistivity, formation of a cascade of subdominant oscillatory branches. Their relevance is amplified in advanced tokamak scenarios and reactor-scale plasmas (e.g., CFETR 1GW SSO), where the interplay between equilibrium shaping, kinetic effects, and current profiles can both promote and modulate infernal activity (Coste-Sarguet et al., 3 Sep 2025, Li et al., 2024, Han et al., 2021).

1. Overview of Governing Physics and Equilibrium Features

Infernal modes arise within the framework of linear resistive MHD, typically around static axisymmetric equilibria described by $\{\mathbf{B}_0, p_0, \rho_0\}$. The key ingredient is an extended region where the local magnetic shear $s(r) = r\,dq/dr / q$ is small or even negative. In such settings, rational surfaces (where $q(r_s) = m/n$) are densely packed or nearly degenerate, causing the coupling between neighboring poloidal harmonics ($m$, $m\pm1$) to intensify.

Shaping effects—elongation ($\kappa$), triangularity ($\delta$)—modify the effective pressure–field-line-bending competition via the generalized Mercier index $DI(r)$: $$DI(r) = D_{i0}(r) + \Delta D_i(\text{shaping}), \quad D_{i0}(r) = -\left( \frac{a}{R_0} \right)\frac{r q'}{q^2}, \quad \Delta D_i(\text{shaping}) = \left( \frac{a}{R_0} \right)\frac{\kappa-1}{2}(1-2\delta)$$ where $a(r) = -(R_0/q(r)) dp_0/dr$ defines the local pressure drive (Coste-Sarguet et al., 3 Sep 2025).

Advanced scenarios, exemplified by the CFETR 1GW SSO equilibrium, feature deeply reversed $q$ in the core ($q_\text{min}\approx 1.12$) and pronounced flat-$q$ plateaus at the edge ($q_\text{flat}\approx 7.4$), predisposing both core and edge to infernal activity (Han et al., 2021).

2. Modal Structure and Dispersion Relations

In low- or reversed-shear regions, the eigenmode structure is governed by an integro-differential equation for the radial plasma displacement $\xi_m(r)$, with coupling to immediate sidebands (resistive layer physics entering via Ohm’s law and rational layer matching).

Ideal Shearless Spectrum

For $|s|\ll1$ and monotonic $q>q_s$, neglecting $q'$ yields the “shearless” infernal spectrum. The generalized dispersion relation reduces to: $$D_0(y^2/\omega_A^2) \equiv \frac{x J_{m+2}(x) - 2(m+1)J_{m+1}(x)}{x J_{m+1}(x) - 2m J_{m}(x)} - \frac{2m+2}{x} = 0$$ with $x^2 = Q y^2/\omega_A^2$, $Q = 4 + 2 r_1^2 D_I^* + m^2$, and $\omega_A$ the Alfvén frequency. Each $x_j$ root corresponds to a discrete mode, forming an unstable spectrum with the fundamental branch ($j=0$) broad across the low-shear region and higher $j$ branches more radially oscillatory (Coste-Sarguet et al., 3 Sep 2025).

Reversed Shear and Generalized $q$ Profiles

When $q(r)$ exhibits a minimum ($q'(r_s)=0,\,q''>0$), field-line bending diminishes, and analytic Bessel-type solutions break down. An energy-minimization plus singular-layer approach yields the full dispersion: $$D(y^2/\omega_A^2) \equiv 2(m+1)^2 \frac{J_m(k r_s)}{k r_s J_{m+1}(k r_s)} + \frac{s_s^2+2(m-1)r_s^2 q_s'' q(r_s)}{q_s^2} - \frac{q_s^2 k^2}{(m+1)^2} = 0$$ which reduces continuously to the ideal case as $s_s,A_q\to 0$ (Coste-Sarguet et al., 3 Sep 2025).

Analytic Instability Thresholds

• Ideal interchange criterion: $D_I^* > [4m(m+2)]/[(m+1)(m+3)]$
• Resistive infernal (reversed-shear): $$|D_I^*| < C_m \sqrt{q'' r_s^2}/(q_s' r_s q_s^{3/2})$$, where $C_m$ depends on $(m,n)$—explicit thresholds are given in (Coste-Sarguet et al., 3 Sep 2025).

3. Physical Drivers and Discrete Mode Cascade

Several mechanisms coalesce in reversed-shear configurations to drive infernal modes:

• Field-line bending suppression: As $s\to 0$, the stabilization term drops, allowing pressure drive to dominate.
• Toroidal sideband coupling: Geometric effects couple $(m,n)$ modes to $m\pm1$ harmonics, enhancing instability even when formally higher order in aspect ratio expansion.
• Cascade spectrum: The eigenmode sequence ($j=0,1,2\dots$) corresponds to increasingly oscillatory radial structures. As long as the resistive growth $y_j \tau_R \gg 1$, subdominant branches remain unstable (Coste-Sarguet et al., 3 Sep 2025).

Kinetic generalizations (KIM) bridge the infernal and kinetic ballooning mode (KBM) branches, especially evident in global gyrokinetic simulations. The transition arises due to the changing proximity of rational surfaces and local shear: KIM occurs where the rational surface is nearest the minimum shear, saturating as $\gamma \sim \beta$ and peaking near optimal $|s| \approx 0.2$. FLR and FOW effects modify thresholds and radial localization (Li et al., 2024).

4. Role of Resistivity, Kinetics, and Numerical Modeling

Resistivity, compressibility, and kinetic terms influence both mode structure and instability thresholds:

• Resistive layer: Inner-layer width $\delta_R \propto (\Lambda/y)^{1/2}$, with $\Lambda$ normalized resistivity. Sufficiently large $y\tau_R$ admits ideal-like infernal roots; otherwise, “resistive interchange” dominates (Coste-Sarguet et al., 3 Sep 2025).
• Finite Larmor radius (FLR)/orbit width (FOW): In kinetic treatment, FLR and FOW set mode envelope width $\ell_e$ and stabilize the small-scale MHD drive above threshold $\beta$. For $|\beta|$ above threshold and $|s|$ small, kinetic infernal modes can reach twice the KBM growth rate (Li et al., 2024).
• Numerical limitations: Conventional MHD codes frequently reduce the perturbed vector potential to $\delta A_\parallel \mathbf{b}$, suppressing parallel magnetic perturbations and altering the Mercier parameter. This omission artificially stabilizes tearing/interchange and masks the infernal branch except at vanishingly small $\beta$. Modular linear solvers that preserve full $\delta\mathbf{A}$, include sidebands up to $m\pm2$, and use adaptive grids reliably capture the entire infernal spectrum in shaped reversed-shear equilibria (Coste-Sarguet et al., 3 Sep 2025).

Characteristic Kinetic and Growth-Rate Data

β (%)	ω<sub>r</sub> (c<sub>s</sub>/R₀)	γ (c<sub>s</sub>/R₀)	ω<sub>r</sub> (c<sub>s</sub>/R₀, rev. shear)	γ (c<sub>s</sub>/R₀, rev. shear)
1.0	0 (ITG)	0.005	0 (ITG)	0.007
1.4	−1.3 (KBM)	0.008	−1.2 (ITG)	0.009
2.0	−1.4	0.01	−1.5 (KIM)	0.018

For $n=10$, $R/L_{Ti}=6.92$, $R/L_n=2.22$ (Li et al., 2024).

5. Experimental Manifestation and Practical Consequences

Analysis of the CFETR 1GW SSO reversed-shear equilibrium provides a clear practical demonstration:

• Profile Features: Deep $q$ reversal in the core, broad edge $q$ plateau.
• Dominant Eigenstructure: Edge-localized infernal harmonics ($m=8$ at $\psi\approx0.98$) possessing peak radial displacement 10–20$\times$ larger than core infernal components ($m=2,3$).
• Stability Boundaries: Edge-localized infernal mode lowers both no-wall and ideal-wall $\beta_N$ limits; the operational target $\beta_N=2.86$ exceeds the computed no-wall limit ($\beta_N=2.76$), necessitating additional stabilization (Han et al., 2021).

Rotational stabilization—particularly edge-localized rotation—is found to be critical. Complete suppression of the dominant edge infernal component at the design wall location requires

$$\Omega_\text{edge} \geq 1.5\%\,\Omega_{A0}, \quad \Omega_{A0} = v_A / R_0$$

where $v_A$ is the Alfvén speed and $R_0$ the major radius. Both uniform and edge-peaked rotation profiles suffice if this condition is met precisely at the edge plateau (Han et al., 2021).

6. Implications for Design and Control in Advanced Tokamaks

The infernal mode spectrum fundamentally constrains the operational space in reversed- and low-shear scenarios:

• Scenario Design: Avoid positioning rational surfaces in flat or reversed-shear regions coincident with strong pressure gradients unless deliberate rapid reconnection (e.g., hybrid sawtooth pacing) is desired.
• Profile Optimization: Raised minimum shear, radial displacement of $q_{min}$, reduction of $R/L_{Ti}$ near weak-shear regions, and controlled shaping (including negative triangularity) serve as suppression strategies (Coste-Sarguet et al., 3 Sep 2025, Li et al., 2024).
• Modeling Requirements: Neglecting parallel field perturbations and limit-sideband physics in stability codes results in overestimating achievable $\beta$ or safety-factor limits, masking critical instability branches.

A plausible implication is that restrained use of advanced bootstrap-current shaping and $q$-profile tailoring is required to balance energy confinement gains against the destabilizing persistence of infernal spectra in such advanced scenarios.

7. Connection to Kinetic and Global Non-MHD Effects

Global gyrokinetic simulations exhibit a seamless transition between kinetic ballooning and infernal branches. In the low-shear/well regimes, electromagnetic inertial instabilities driven by ion pressure gradient ($d\beta/dr$) combine features of both classical infernal and kinetic ballooning modes. As $|s|$ decreases, the mode envelope localizes about the shear minimum, with the growth rate maximized for moderate but finite $|s|$ before declining as stabilization from residual field-line bending (and FLR) overcomes the drive (Li et al., 2024).

This continuous spectrum, dependence on rational surface separation $\Delta r_{m,m+1}$, and the associated physics point to the necessity of adopting global (as opposed to ballooning local) analyses and kinetic-MHD models in reversed-shear stability studies for reactor-relevant plasmas.