等离子体的不稳定性

高温等离子体如果放着不管，就会自发地破坏自身的形状，试图发狂。这种「容易发狂」的性质就是等离子体不稳定性（plasma instability）。不稳定性会扰乱约束，是限制聚变堆性能的最大因素之一。本页从「不稳定性是什么」的直觉出发，逐一理解托卡马克中会造成问题的代表性模式，再层层递进到控制这些不稳定性的研究前沿。

先建立直觉（高中水平）

我们先用身边的例子来思考什么是不稳定。把铅笔平放在桌上，即使稍微推一下，它也会回到原位或停在原地。这就是稳定。然而，如果试图让铅笔用笔尖立起来，哪怕极其微小的倾斜都会不断变大，最终倒下。这就是不稳定。

关键在于，在不稳定的状态下，「微小的偏差会自己成长」。不需要有人用力去推。只要有一点点扰动，就会以此为契机而擅自变大。等离子体的不稳定性也是一样，被约束的等离子体中储存着大量的能量，这些能量试图转移到更低的状态，从而让微小的扰动一下子成长起来。

还有一个非常适合把握直觉的例子。在杯子里倒入水，再在上面轻轻浮上一层油，就会形成轻的油在上、重的水在下这样的稳定排列。然而反过来，如果试图把重的水放在轻的油之上，水就会落进油里，油则会穿透水浮上来，界面变得乱七八糟。重的东西在轻的东西之上的状态是不稳定的。这种现象被称为瑞利-泰勒不稳定性（Rayleigh-Taylor instability）。

在等离子体中，扮演这个「重的水」角色的有时是高压的等离子体，扮演「轻的油」角色的则是磁场。用磁场压住等离子体时，正好形成把水放在油上那样的配置，等离子体会试图冲破磁场向外喷出。把聚变堆中发生的许多不稳定性都看作这种瑞利-泰勒不稳定性的亲戚，就更容易把握整体图景。

在托卡马克这种装置中，人们在环形容器内对等离子体施加螺旋状的磁场来进行约束。只要这个螺旋的缠绕方式、流经等离子体的电流以及等离子体压力之间的平衡稍有不当，整个等离子体就会软塌塌地弯曲、局部膨胀或者被撕裂。从下一节开始，让我们给这些发狂的方式逐一命名并逐一观察。

用物理来理解（本科水平）

把不稳定性作为物理来处理时，首先要了解的是把等离子体视为一种流体、与磁场一起来描述其运动的电磁流体力学（magnetohydrodynamics, MHD）框架。MHD 处理的宏观不稳定性大致可分为电流驱动型和压力驱动型。前者以流经等离子体的电流为能量源，后者以等离子体的压力为能量源。

电流驱动型的代表是扭曲模（kink mode）。流过电流的细柱，由于自身产生的磁场，有时会陷入「弯成螺旋状比笔直更划算」的处境。就像电线扭转起来打成疙瘩一样，整个等离子体柱会变形成螺旋状。它的近亲有腊肠模（sausage mode），这种模式会让等离子体柱在一些地方产生收缩，形成像腊肠那样有粗有细的部分。在变细的部分，磁场增强而进一步收紧，收缩得以加速。

在这里登场的最重要参数就是安全因子（safety factor） $q$ 。 $q$ 表示磁力线在沿环的长方向（环向）绕几圈的期间，会沿短方向（极向）绕 1 圈。用式子来写，在圆形截面的近似下为

q(r) = \frac{r B_\phi}{R B_\theta}

其中 $r$ 是磁力线所经过的半径， $R$ 是环的大半径， $B_\phi$ 是环向磁场， $B_\theta$ 是极向磁场。 $q$ 越大，磁力线的扭转越平缓，等离子体越稳定。之所以叫「安全」因子，正是因为它是衡量稳定性余量的指标。

对于扭曲模，已知圆柱等离子体的 $m = 1$ 外部扭曲模会在等离子体表面的安全因子 $q_a > 1$ 时被稳定化，这被称为克鲁斯卡尔-沙夫拉诺夫极限（Kruskal-Shafranov limit）。在实际的托卡马克中，为了避免更高次的模式等，通常以 $q_a \geq 3$ 左右作为稳定运行的标准。至于等离子体内部，会出现与 $q = 1$ 面共振的内部扭曲模，这会成为锯齿振荡（sawtooth oscillation）的原因。在锯齿振荡中，中心区域的温度缓慢上升之后又急剧跌落，如此循环反复，因其波形类似锯齿而得名。周期大约为 10〜100 ms。

压力驱动型的代表是气球模（ballooning mode）。在环的外侧（环外侧），磁场的弯曲方向对等离子体不利，即所谓坏曲率的区域。在这里，高压的等离子体会像气球在薄弱处鼓起那样，试图局部膨胀。正如其名，就是气球（balloon）的形象。等离子体能把压力提高到什么程度，用归一化贝塔（normalized beta） $\beta_N$ 这个量来衡量，经验上已知存在一个被称为特罗荣极限（Troyon limit）的上限（在壁较远的情况下 $\beta_N$ 约为 2.8，在实际的托卡马克中可达 3〜3.5 左右）。所谓贝塔，是等离子体压力除以磁场压力的比值，是表示约束效率的重要指标。

还有一种因电阻起作用而发生的不稳定性。在理想的等离子体中，磁力线与等离子体一起运动，磁力线不会重新连接。然而现实的等离子体具有微小的电阻，由此会发生磁力线断开并重新连接的磁重联（magnetic reconnection）。当这种重联在共振面上推进时，那里就会发展出撕裂模（tearing mode）。撕裂（tear）是「扯裂」的意思，表示磁力线被扯裂并重新连接的样子。

深入理论（研究生水平）

作为撕裂模的结果而出现的是磁岛（magnetic island）。原本呈同心圆状层层嵌套的磁面，在共振面周围局部地重新连接，从截面看会形成像岛屿一样封闭的磁力线结构。由于磁岛内部的磁力线短路，温度和密度会横穿岛屿迅速均匀化。所谓约束就是在中心和边缘之间保持很大的梯度，因此岛内梯度的消失就意味着约束性能相应地下降。

尤其棘手的是新经典撕裂模（neoclassical tearing mode, NTM）。托卡马克等离子体中存在由压力梯度而自发流动的自举电流（bootstrap current）。由于磁岛内部压力梯度消失，这个电流也会缺失。缺失的电流会朝着进一步扩大磁岛的方向起作用，因此一旦形成小岛（种子磁岛）并超过临界尺寸，岛就会开始自我成长。NTM 是限制高贝塔运行的主要因素，尤其是立在 $q = 3/2$ 面和 $q = 2$ 面上的岛会成为问题。长周期化的锯齿（怪兽锯齿）有时也会制造出 NTM 的种子磁岛。

扭曲类中也有涉及电阻效应的重要模式。电阻壁模（resistive wall mode, RWM）是指本应被周围的完全导体壁稳定化的外部扭曲模，由于壁的有限电阻而随时间渗出，缓慢地成长起来。已知 RWM 可以通过等离子体的旋转或主动的反馈控制来稳定化。

作为局限于边界的不稳定性，ELM（edge localized mode, 边缘局域模）很重要。在高约束状态即 H 模下，等离子体边缘会立起一道陡峭的压力壁，即台基（pedestal）。这个陡峭的梯度和流经其中的大自举电流，会引发同时兼具压力驱动和电流驱动的剥离-气球不稳定性，周期性地吐出台基的能量。这就是 ELM。尤其是 Type I ELM，频率大约为 10〜100 Hz，一次会释放台基能量的百分之几到百分之十几。释放出的能量会作为瞬间热负荷集中到偏滤器上，在大型装置中会导致壁材料损伤，成为严重的问题。

最剧烈的现象是破裂（disruption）。这是等离子体放电突然崩溃而终止的大事件，与其说是单一的不稳定性，不如说是多个不稳定性连锁发生的灾变。当密度提得过高而超过格林沃尔德极限（Greenwald limit），或贝塔提得过高而超过贝塔极限时，大的磁岛会成长起来使磁面破坏，最终导致破裂。过程分两个阶段推进。首先在热淬灭（thermal quench, 约 1〜10 ms）中储存的热能一下子丧失，接着在电流淬灭（current quench, 约 10〜100 ms）中等离子体电流急剧减小。电流的急剧变化会在容器中产生巨大的电磁力，晕电流和逃逸电子（runaway electron）会对设备造成严重损伤。如何避免破裂，以及在无法避免时如何减轻损害，是聚变堆设计的核心课题之一。

另一方面，一点一点地侵蚀约束的是微观不稳定性（microinstability）。这不是破坏整个装置那样的大变形，而是在等离子体中以离子或电子的尺度发生的细微扰动，会产生湍流（turbulence）。湍流会横穿梯度把热和粒子输运出去，成为远远超出经典理论预测的反常输运（anomalous transport）的真凶。代表性的有由离子温度梯度驱动的离子温度梯度模（ion temperature gradient mode, ITG）、捕获电子参与的捕获电子模（trapped electron mode, TEM）、以及电子尺度的电子温度梯度模（electron temperature gradient mode, ETG）。它们与 MHD 不稳定性不同，是在把等离子体连粒子的速度分布也一并处理的动理学理论（kinetic theory）框架下描述的，基于回旋动理学（gyrokinetics）的大规模模拟成为理论研究的主角。与输运的关系将在输运页面进一步详细讲解。

研究前沿（博士水平）

在当前的研究中，重心正从「发生之后再抑制」这些不稳定性，逐渐转向「在发生之前就预测并避免」的方向。

在 NTM 控制方面，电子回旋电流驱动（electron cyclotron current drive, ECCD）已确立为标准手段。瞄准磁岛的 O 点（岛的中心）局部地驱动电流，补上缺失的自举电流，从而阻止岛的成长。为了持续把束流准确打在瞄准的位置上，实时控制成为研究的焦点。

在 ELM 控制方面，共振磁扰动（resonant magnetic perturbation, RMP）被寄予厚望。这是一种从外部线圈施加小的螺旋磁场、稍微扰乱边缘磁面，把大的 Type I ELM 变成细碎而温和的释放，或者完全抑制的方法。通过弹丸注入进行的 ELM 调步（人为诱发小的 ELM 以把能量释放分成小份的手法）也在并行研究。在 ITER 这样的大型装置中，由于不受控的 Type I ELM 的热负荷无法承受，因此完全抑制或充分缓解 ELM 成为必须实现的要求。

破裂预测是近年来机器学习被格外积极应用的领域。人们正在推进从大量诊断信号中实时检测破裂的前兆、趁还来得及时转入缓解动作的研究。作为缓解手段，人们开发了用大量气体注入（massive gas injection）或碎裂弹丸注入（shattered pellet injection）一下子注入杂质，把储存的能量作为辐射温和地释放掉，同时抑制逃逸电子产生的手法。逃逸电子本身的产生机制及其抑制，也是活跃的研究课题。

在微观不稳定性和湍流领域，人们持续开展把回旋动理学模拟与实验测量相互对照、试图从第一性原理预测湍流输运的验证（validation）工作。湍流自身产生的带状流（zonal flow）削弱湍流这一自我控制机制，以及约束突然变好的输运垒（transport barrier）的形成机制，至今仍是核心问题被研究着。论文中频繁出现 ITG, TEM, ETG, zonal flow, gyrokinetics, ExB shear 等关键词。

理解自测

第 1 题「不稳定」并不是被外力猛推而损坏。那么不稳定状态的本质特征是什么呢？

第 2 题安全因子 q 大的时候，为什么等离子体容易变得稳定呢？

第 3 题磁岛一旦形成，为什么约束性能会下降呢？

第 4 题新经典撕裂模（NTM）为什么会「岛一旦形成就自己成长」呢？

第 5 题微观不稳定性（如 ITG 和 ETG 等）与扭曲、气球这类 MHD 不稳定性在本质上有什么不同呢？