← 返回研究总览

已发表于《天体物理学杂志》· 2026

剪切驱动的磁化湍流中的粒子加速

在湍流剪切层中,粒子一次次受到微小的随机扰动。它们怎样逐渐获得能量,形成高能粒子群?

三幅并排的模拟图,从左至右分别显示剪切层达到统计稳态后的气体密度、速率和磁场强度。
剪切驱动的湍流层达到统计稳态后,气体密度(左)、速率(中)和磁场强度(右)呈现不同的空间结构。在剪切最强的区域,磁场受到拉伸和折叠,在右图形成明亮的细丝。

当相邻的等离子体以不同速度流过彼此时,交界面上的速度剪切可引发开尔文–亥姆霍兹不稳定性。界面先卷成涡旋,随后涡旋破碎,发展成湍流,将磁场搅得错综复杂。喷流的边缘、从落入星系团的星系中剥离出来的气体,以及星系风的边界,都能出现这样的剪切层。我想知道的是,穿过这片磁场的带电粒子会发生什么变化。

我们模拟了仅由速度剪切持续驱动的亚声速、非相对论性磁化湍流,并计入粒子对流体的反作用。是否持续向湍流供能,对粒子加速至关重要:只要驱动持续,粒子就能不断获得能量;一旦停止驱动,湍流中储存的能量会在几个大涡翻转时间内耗散,粒子加速也随之停止。

原因藏在粒子回旋轨道的形状里。湍流电场会扭曲轨道:在电场力做正功的半个回旋周期中,粒子沿力的方向走得更远,获得的能量更多;在电场力做负功的另半个周期中,粒子沿反方向走得较短,损失的能量较少。经过多次随机扰动,这种不对称使粒子平均获得净能量,且平均增益与剪切速度的平方成正比,正是二阶费米加速的特征。

初始能量相同的粒子逐渐形成明显的非热高能尾部。反复穿越剪切层的粒子以近似几何布朗运动的方式积累能量,动量因而呈对数正态分布。当回旋半径超过湍流层厚度时,粒子的运动方向趋于与背景磁场垂直。

发表于 《天体物理学杂志》 1009,247(2026)。
二维 MHD-PIC 模拟采用 Sun 与 Bai(2023) 介绍的 Athena++ MHD-PIC 模块。

不同时刻的粒子动量分布:最初的尖峰逐渐展宽,形成向高能端延伸的尾部。
粒子动量分布随时间的变化,颜色表示不同时刻。初始单能分布的尖峰逐渐变平,并形成宽广的非热尾部。只要湍流持续受到驱动,高能端的截断就会继续向更高的动量移动。
湍流动能和磁能的功率谱,横轴为波数,参考线的幂律指数为 −2。
统计稳态下的湍流动能谱和磁能谱。虚线标示惯性区的 k−2 参考标度。在小尺度上,磁能超过动能。
上图是一颗粒子穿越剪切层 600 次的累计能量变化,数值从 0 起伏上升到约 7。下图显示每次穿越的能量得失,红色表示增加,黑色表示减少,两者的幅度相近。
一颗粒子在反复穿越剪切层的过程中如何积累能量?下图记录每次穿越时的能量变化,红色的增益和黑色的损失看起来几乎对称。把它们逐次加起来,就得到上图的曲线。粒子每次获得的能量平均只比损失的略多一点,经过 600 次穿越,能量却已增至初始值的数倍。
双对数图展示粒子的平均相对能量增益随剪切速度的变化。拟合得到的幂律指数为 2.01 ± 0.02。
对不同的剪切速度重复测量,发现粒子的平均相对能量增益与剪切速度呈幂律关系,指数为 2.01 ± 0.02,与二阶费米加速预言的平方关系吻合。
数百条粒子动量轨迹从同一初值向外散开,颜色由蓝到红表示末态动量由低到高;右侧是由此形成的对数正态分布。
图中各条曲线对应反复穿越剪切层的不同粒子。它们的初始动量相同,随后逐渐分散,形成右侧的对数正态分布。这种分布源于许多微小随机因子的连乘。虚线给出了测得的均值和理论预期的均值。