分层等离子体电流片在重力场中的电阻性不稳定性：撕裂模的引力修正与G模转变

《Journal of Plasma Physics》：Resistive instabilities of current sheets in stratified plasmas with a gravitational field

　　

磁重联可以通过撕裂模不稳定性自发发展，这一机制常被用来解释聚变装置中的破坏性不稳定性、太阳耀斑、盔状流罩顶端的周期性密度扰动以及地球日侧磁层顶的通量传输事件。然而，在许多此类环境中，重力、磁场曲率或其他形式的加速度常常导致重等离子体覆盖轻等离子体的有效重力场构型，其中嵌入了电流片。本文研究了在恒定重力加速度作用下存在密度梯度时，平板电流片对重联模的线性稳定性。研究表明，分层和重力的存在在有利和不利分层两种情况下都会改变撕裂模不稳定性的性质。有利分层抑制重联，而不利分层则强烈地使撕裂模失稳。此外，研究还表明，即使对于弱不利分层，在S ? 1条件下，经典的常ψ区域实际上也不存在。取而代之的是，重力修正的撕裂模逐渐转变为G模，这是一种重力驱动的重联模，其增长率标度为S-1/3。因此，不利分层只允许快速重联模存在。

论文解读

一、研究背景与科学问题

磁重联被认为是实验室、空间和天体等离子体中许多爆炸性现象背后的动力学机制（Yamada, Kulsrud & Ji 2010; Pucci et al. 2020b; Ji & Daughton 2022）。通过磁拓扑重构，重联将储存的磁能转化为等离子体加热、粒子加速和整体流动。这一机制是磁约束聚变装置（如tokamak锯齿崩溃）、太阳耀斑和日冕物质抛射、行星磁层内的地磁暴和亚暴以及吸积盘和相对论性喷流等高能天体物理现象的核心。

撕裂模不稳定性最早由Furth, Killeen & Rosenbluth（1963）分析，是电阻性磁流体力学（resistive magnetohydrodynamic, MHD）框架下电流片内磁重联发展的基本机制。过去几十年，撕裂模已在日益复杂的等离子体条件下被广泛研究，包括平衡剪切流和/或粘性、Hall和双流体效应、部分电离等离子体以及无碰撞或动力学机制。

然而，在许多自然和实验室环境中，重联可能发生在分层等离子体中，并伴随重力、磁场曲率或其他有效加速度的作用，有时导致重等离子体覆盖轻等离子体的构型嵌入电流片内。在地球磁层顶，太阳风动压的快速变化使边界穿过磁鞘致密等离子体与磁层稀薄等离子体之间的陡峭密度梯度加速，产生可驱动Rayleigh-Taylor模的有效重力，并可能在南向行星际磁场条件下影响重联。在日球层顶也存在类似条件。在实验室中，tokamak和stellarator等磁约束装置也容易出现交换型不稳定性和其他重力驱动的电阻性不稳定性。

然而，这些早期工作依赖于限制性假设——如常ψ、弱重力、弱磁剪切或渐近波长极限——从而限制了它们捕捉撕裂模如何被重力修正以及如何向重力驱动不稳定性转变的能力。

二、研究内容与结论

本研究在平面平板电流片中，通过引入垂直于电流层方向作用的重力加速度，研究了有利和不利分层对重联不稳定性线性阶段的影响。有利和不利分层分别指密度梯度与重力加速度平行或反平行的情形。虽然沿重力场方向具有负密度梯度的构型（不利分层）在理想情况下可能对交换模或Rayleigh-Taylor模不稳定，但研究将注意力限制在分层足够弱以满足Suydam稳定性判据的范围内。

在均匀电阻率条件下，已有研究表明平板片状箍缩在常ψ近似下对两种电阻性不稳定性不稳定：一种是当重力可忽略时占主导的长波长撕裂模，其增长率随Lundquist数标度为S^-3/5；另一种是短波长重力驱动模（G模），其增长率标度为S^-1/3。最近，Hopper, Wood & Bushby（2024）研究了有利分层对平板几何中撕裂模线性增长率的影响，发现增大的稳定化密度梯度抑制撕裂扰动的增长。本研究遵循其方法求解电阻层方程并假设线性密度剖面，但将其分析扩展到包含不利分层，允许重力驱动不稳定性的出现。

研究人员推导了描述从撕裂模主导区域到G模连续转变的一般色散关系，考察了增长率对分层强度（有利和不利）和Lundquist数的依赖关系。此外，还数值研究了具有两个等离子体区域之间平滑过渡的更现实平衡密度剖面，以评估增长率对密度剖面的敏感性。

三、关键技术方法

研究人员采用边界层方法（boundary-layer approach）对撕裂模方程进行渐近匹配分析。具体而言，在理想外区（outer region）中忽略电阻率和等离子体惯性，得到关于流函数（stream function）的二阶常微分方程；在电阻内层（inner layer）中求解包含电阻效应的薄层方程。通过Frobenius展开技术获得外区解在电阻层附近的行为，进而推导色散关系。数值方面，研究人员求解了更现实的tanh密度剖面下的本征值问题，以检验解析结果对密度剖面选择的敏感性。

四、研究结果

模型方程与平衡构型

研究人员从忽略等离子体可压缩性并包含均匀电阻率的电阻MHD方程组出发，采用Harris型电流片平衡，磁场为B? = B? tanh(x/a)?，密度剖面为线性形式ρ?(x) = ρ??(1 + εx/a)，其中ε ? 1以保证Boussinesq近似的适用性。通过线性化并引入流函数和磁通函数，得到两个耦合方程。分层参数定义为A = ε(ag/V_A2)，其正负分别对应不利和有利分层。

外区解与色散关系

在外区，理想MHD方程可简化为关于流函数的二阶方程。通过Frobenius展开，得到流函数在电阻层附近的渐近行为，其中指数由l = √(1/4 - A)决定。为保证幂律行为，需满足A < 1/4，这对应于Suydam稳定性条件。通过渐近匹配外区和内区解，推导出一般色散关系，描述了从撕裂模到G模的连续转变。

有利分层对撕裂模的影响

有利分层（A < 0）抑制撕裂模不稳定性。随着稳定化密度梯度的增强，撕裂扰动的增长被抑制，增长率降低。在有利分层渐近区域，增长率标度降至γτ_A ~ S^-1，此时不稳定性增长时间尺度与平衡电流片的电阻扩散时间尺度相当，固定平衡假设不再自洽。

不利分层与G模转变

不利分层（A > 0）强烈地使撕裂模失稳。研究表明，即使对于弱不利分层，在S ? 1条件下，经典的常ψ区域实际上也不存在。取而代之的是，重力修正的撕裂模逐渐转变为G模，这是一种重力驱动的重联模，其增长率标度为S^-1/3。因此，不利分层只允许快速重联模存在。

数值验证与密度剖面敏感性

研究人员数值求解了具有tanh密度剖面的更现实平衡构型，以评估增长率对密度剖面的敏感性。数值结果与解析色散关系在相应参数范围内一致，验证了理论分析的可靠性。

五、讨论与结论

研究表明，分层和重力在有利和不利两种情况下都会显著改变撕裂模不稳定性的性质。有利分层抑制重联，而不利分层强烈地使撕裂模失稳并使其转变为G模。这一发现对理解磁层顶、日球层顶等自然环境中密度梯度与有效重力共同作用下的磁重联过程具有重要意义。不利分层条件下快速重联模的存在意味着在这些环境中，磁重联可能以比传统撕裂模更快的速率发展，从而影响相关的能量释放和输运过程。