电子压强等向化控制的磁发电机与魏贝尔不稳定性在无碰撞10阶矩等离子体模型中的研究

《Journal of Plasma Physics》：From Weibel seeds to dynamo beyond pair-plasmas

　　

摘要：研究人员采用无碰撞10阶矩流体模型（10-moment collisionless plasma model）并结合Gkeyll求解器，研究电子热流闭合（electron heat-flux closure）所调控的电子压强等向化（electron-pressure isotropisation）对磁发电机（magnetic dynamo）与魏贝尔不稳定性（Weibel instability）的控制作用。通过在弱碰撞星团内介质（intracluster medium, ICM）相关参数下开展三维受迫湍流模拟，研究人员以Galloway–Proctor（GP）流驱动亚声速剪切流，初始磁场设为无磁化状态，由电子魏贝尔不稳定性产生种子磁场，随后过渡到发电机放大阶段。研究发现，电子闭合参数k0,e控制磁场阻尼率γB?∝k4/k0,e（等效于调控有效磁雷诺数Rm），离子闭合参数k0,i控制流动阻尼率γU?∝k2/k0,i（等效于流体雷诺数Re）。当质量比mi/me=100时，电子与离子动力学充分解耦：电子魏贝尔不稳定性由电子压强各向异性Δe驱动，饱和后磁场能谱峰值向低波数迁移，进入与外流翻转时间tturn相当的发电机指数增长阶段，最终磁场动能与流动动能趋于近均分。增大k0,e（减弱电子等向化强度）可降低有效磁场阻尼、提高等效Rm，使发电机增长率增大且魏贝尔阶段相对缩短，系统更趋近磁流体力学（MHD）行为；反之较小k0,e时魏贝尔阶段显著。研究还表明10阶矩模型中闭合参数可灵活复现不同微观耗散特性，为跨越种子场生成与发电机放大的多尺度模拟提供了介于全动理学（fully kinetic）与MHD之间的有效途径。

研究背景方面，宇宙中磁场起源与放大是天体物理长期未解难题，星系团内介质（intracluster medium, ICM）等弱碰撞等离子体（weakly collisional plasma）中，传统磁流体力学（magnetohydrodynamics, MHD）无法描述无碰撞动理学效应如压强张量各向异性、微尺度不稳定性等；全动理学粒子模拟（如PIC）虽能捕捉动理学但受限于质量比、尺度分离与计算成本，难以同时解析种子磁场生成的微观魏贝尔不稳定性（Weibel instability）与宏观发电机（dynamo）放大过程。现有成对等离子体（pair-plasma）研究结论对电子–质子真实质量比适用性不明，因此有必要在较大离子–电子质量比下采用兼顾动理学特征与计算效率的模型开展研究。研究人员利用无碰撞10阶矩流体模型（10-moment collisionless fluid model，每物种演化密度、动量、全压强张量共10个矩，热流由闭合关系给定）结合Gkeyll框架，在质量比m_i/m_e=100、Prandtl数P_m?1（R_m?1，R_e<1，对标ICM参数）下模拟亚声速剪切湍流中的磁生成全过程，分析电子闭合参数k_0,e（控制电子压强等向化速率与等效磁雷诺数R_m）与离子闭合参数k_0,i（控制流动阻尼与等效流体雷诺数R_e）的影响，明确电子/离子动力学解耦条件，并检验该模型对魏贝尔不稳定性与发电机过程的刻画能力，论文发表于《Journal of Plasma Physics》。

关键技术方法：研究人员采用无碰撞10阶矩流体模型（10-moment collisionless fluid model）在Gkeyll求解器中实现，每物种演化密度、三动量分量、六个独立压强张量分量（共10矩），热流采用压强等向化闭合（pressure isotropisation closure）引入物种依赖的闭合参数k_0,α（α为离子i、电子e）来近似动理学等向化效应并确定有效阻尼率；模拟为三维周期性盒，采用有限体积法离散，外力驱动为Galloway–Proctor（GP）时变流形式，初始磁场按多模余弦叠加赋予微小种子（β_e^?1~10^?6），初始等离子体参数为ICM对标值（T_e=T_i=1 keV，n_e=2.3×10²⁸m^?3，m_i/m_e=100，u₀=0.35c_s），盒尺度L=70δ_i（δ_i为离子惯性长度），分辨率N=140³；另开展一维线性阻尼测试扫描k_0,α与波数k以标定γ_U^?、γ_B^?标度，二维魏贝尔基准测试对照动理学结果调参，质量比扫描（m_i/m_e从1到100）检验电子–离子解耦，电子闭合参数扫描k_0,e/k₀={1,2,4,8,32}分析对魏贝尔阶段与发电机增长率的影响。

研究结果：首先为基线情形（baseline case，k_0,e/k₀=1，k_0,i/k₀=13.6），磁场能量演化分三阶段：极短时（t/t_turn<0.1）为电子魏贝尔不稳定性主导的快速指数增长（峰值瞬时γ_B/ω_pe对应k~Δ_e^1/2/δ_e），磁场呈丝状结构位于剪切流双曲固定点强应变区；随后进入与外翻转时间相当的发电机指数增长阶段（γ_Bt_turn≈0.76），能谱峰值由高k移向盒尺度k₀，磁结构失去丝状走向大尺度，最终饱和于磁场与流动动能近均分；流动动能先略降后恢复。压力各向异性显示Δ_e在魏贝尔阶段先升后因磁场反作用于各向异性而稳定在较低值，Δ_i持续受剪切流产生并由离子闭合等向化控制，电子与离子Δ解耦。形态上魏贝尔阶段B集中在剪切最强区且沿不稳定分隔线排列，发电机阶段大尺度结构显现，饱和态近各向同性但有沿GP流优势方向（x向）的排列。

其次为增强电子等向化（increased electron pressure isotropisation），即增大k_0,e/k₀={2,4,8,32}的扫描结果：随k_0,e增大，魏贝尔阶段减弱至消失（k_0,e/k₀≥4时无清晰魏贝尔峰），磁场阻尼率γ_B^?∝1/k_0,e减小即等效R_m增大，发电机指数增长率γ_Bt_turn从~0.76（k_0,e/k₀=1）升至更高（k_0,e/k₀=32时更大），与高R_m快发电机行为一致；饱和能谱仍呈两段以上斜率变化，低k（≤4k₀）较平（~k^?1.6），更高k更陡（~k^?5.5）；流动动能演化与基线类似，x向分量仍占优。魏贝尔特征波数k_EB,max与最大生长率γ_B,max随k_0,e增大下降，Δ_e,max亦降。

附录部分：附录A为10阶矩系统中的魏贝尔不稳定性，对照零热流闭合色散关系与动理学Z函数形式，发现本文采用的等向化闭合在Δ_e小时与原闭合（k_0,α=0）差～ζ_α^?4；用成对等离子体、加速度a_y=a₀sin(k₀x)驱动剪切流测试表明，取k_0,e/k₀≈1.3可使Δ_e,max与k_EB,max匹配动理学估计，但γ_B,max需k_0,e/k₀≈3才匹配；能谱饱和时出现奇谐波对应电流片结构，与Mart’yanov等结论一致。附录B为压强各向异性与归一化压强演化，在p_⊥,e/p_∥,e–β_e,∥图上，魏贝尔阶段电子仅盒尺度附近微弱磁化（β_e,∥~10³），饱和机制属Alfvén电流 regime而非磁化饱和；动态期非磁化剪切仍可产生一定Δ_e（尤其低β_e,∥处p_⊥,e/p_∥,e>1似绝热行为），饱和阶段Δ_e普遍压低；离子Δ_i～k₀u_i/(k_0,iv_th,i)平衡剪切产生与闭合等向化，与大k_0,i基线符合好；极小k_0,e时Δ_e在动态期可达较高值但饱和仍被磁结构散射等限制。

讨论与结论部分：研究人员总结认为无碰撞10阶矩流体模型结合物种依赖的等向化闭合参数k_0,α，可分别调控流动阻尼（∝k²/k_0,i，等效R_e）与磁场阻尼（∝k⁴/k_0,e，等效R_m），在质量比m_i/m_e=100时电子与离子动力学充分解耦，电子主导磁场阻尼与魏贝尔不稳定性，离子主导流动阻尼；模拟复现了由电子魏贝尔不稳定性生成种子磁场（Δ_e驱动，丝状结构，生长率与波数标度与理论一致）、过渡到外包翻转时间尺度的发电机指数增长、最终近均分的全过程，与成对等离子体动理学结果定性相似但增加了电子–离子解耦特征；增大k_0,e（弱电子等向化）降低磁场阻尼、升高等效R_m，使发电机增长率增大、魏贝尔阶段缩短、系统更趋近MHD行为。该模型计算成本远低于全动理学，允许质量比扫描与较真实参数设置，为弱碰撞等离子体中跨种子场微观生成到宏观发电机放大的研究提供了有效途径；但闭合固定形式难以完全适应从无磁化到磁化过渡中动理学效应（相混合、热输运）的定性变化，未来可结合理论或数据驱动改进闭合，也可扩展多组分、引入自引力等提升真实性。论文发表在《Journal of Plasma Physics》。