约束与权衡16 分钟

太阳外层为什么能维持高温,波动的能量怎样交给粒子?

日冕维持高温,需要太阳低层输入的机械与电磁能在合适位置交给粒子,并补偿辐射、热传导和外流带走的能量。Alfvén波能够传能,但波存在、能流足够、局部耗散与不可逆加热是不同环节。本文以DKIST的小尺度扭转波观测为入口,计算视线叠加和谱线宽度的关系、反射与级联的能量收支,以及有效速度扩散和间歇波动的加热含义;再比较近太阳探测器的粒子反演证据。重点限定安静日冕与开放磁场的相关过程,不把一种区域、一种读出或一套闭合公式升格为全部日冕加热的答案。

温度高,不等于能量只能从更热的地方传来

太阳可见表面的温度是几千开尔文,外侧稀薄日冕却能达到百万开尔文量级。这里并不要求热量自发从冷处沿温度梯度倒流到热处。太阳低层的运动可以先对磁场做功,把能量装进磁场应力、波动和流动,再在别处交给粒子。1

因此,要解释日冕高温,至少要讲清能量怎样进入、经过哪里、在哪里改变形式,以及怎样离开。日冕会辐射,热流可以沿磁场传导,太阳风还会带走能量。维持高温需要相应供给;只看见波,或者测到某处温度较高,还不能把这几笔账全部关上。4

本篇主要追踪安静日冕及开放磁场相关过程。入口观测位于日面外约0.1个太阳半径处,狭缝跨过开放区、安静环和其他结构;后面的粒子证据来自更远处近太阳风中的特定时段。它们约束能量链的不同环节,不能当成一次对同一团低日冕气体的完整追踪。

能量从低层运动进入磁场与波,再经小尺度场和粒子作用改变形式。各过程可能并存;图中箭头不表示已经逐段测得了全部功率。
能量从低层运动进入磁场与波,再经小尺度场和粒子作用改变形式。各过程可能并存;图中箭头不表示已经逐段测得了全部功率。

Alfvén波能搬运能量,却不自动产生热

磁化等离子体受到横向扰动时,磁场的恢复作用与流体惯性可以共同支持波传播。最简单的均匀介质中,Alfvén速度为:

$$ v_A=\frac{B_0}{\sqrt{\mu_0\rho}}. $$

对于一个线性行进Alfvén波,速度扰动与磁场扰动相联系。平均动能和磁能相等,总波能密度与$\rho\langle\delta v^2\rangle$有关,能流近似为:

$$ F_w=\rho\,v_{\rm rms}^2\,v_A. $$

这里使用均方根速度;若改用单一正弦波的峰振幅,公式需要相应的二分之一因子。对于宽频随机信号,“均方根乘以$\sqrt2$”是一种振幅约定,不能直接叫作实际最大峰值。

在理想、均匀、没有损失的线性模型里,波可以不断传过区域,而不把能量留在那里。真正与当地加热有关的,是波能通量的变化、其他能量形式的交换以及最后的损失项。还要检查一部分能量是否只是反射回去了、进入别的波模态,或被流动带走。34

望远镜看到的“小速度”,可能是叠加后的结果

Morton等利用DKIST的Cryo-NIRSP仪器观测1,074.7 nm的Fe XIII谱线,在细结构的整体摆动之外,识别出两侧常呈反相的残余多普勒信号。作者将其解释为小尺度扭转Alfvén运动,并用线性波模型、三维模拟及视线叠加计算检验这一解释。18

这不是简单地看见磁通管在照片上转圈。日冕对这条谱线近似光学薄,沿同一视线可能有许多发光结构。谱线中心给出的是按发射强度加权的平均速度,方向相反的运动会互相抵消。

观测处理还明显抑制了小于约2 Mm的空间变化和短于约30秒的时间变化。因此,这里的“小尺度波动”仍是日冕结构尺度上的信号,不能直接当成已经测到了粒子回旋尺度的耗散。1

在同一温度、相同本征线宽的高斯成分近似下,谱线质心和方差满足:

$$ \begin{aligned} \bar u&=\sum_i w_i u_i,\\ \sigma_{\rm line}^2 &=\sigma_{\rm thermal}^2+ \sum_iw_i u_i^2-\bar u^2, \qquad\sum_iw_i=1. \end{aligned} $$

已知的仪器展宽还要另外处理。这个关系说明:质心的小幅变化并不保证内部运动很小;被抵消的不同速度会出现在谱线宽度中。它也说明,谱线宽了不一定就是离子温度升高。

原研究估计非热线宽时,先把Fe XIII离子温度取为谱线形成峰值附近的约1.6 MK,再扣除热和仪器贡献。这样的线宽估计不能同时被当成完全独立的温度测量。1

本库把相位随机、振幅相同的发光成分叠在一起。在独立等权的设定下,质心均方速度的平均值约按$1/N$减小,成分自身的均方速度却没有随$N$减小。直接合成谱线并积分其矩,也验证了上面的方差关系。9

本库读出示意。左上:随机相位叠加抑制质心信号。右上:相同温度下,内部速度差让混合谱线变宽。左下:不同混合权重也能产生反相残差。右下:发射强度权重与质量权重不同,同一组速度可得到不同均方值。均非DKIST原始数据重处理。
本库读出示意。左上:随机相位叠加抑制质心信号。右上:相同温度下,内部速度差让混合谱线变宽。左下:不同混合权重也能产生反相残差。右下:发射强度权重与质量权重不同,同一组速度可得到不同均方值。均非DKIST原始数据重处理。

还有两层限制。第一,补回谱线宽度,得到的仍是发射强度加权的速度统计。在简单的同温碰撞辐射近似下,发射强度约随密度平方变化,而质量只随密度变化;两种权重不能随意互换。本库两个等体积成分的例子中,发射权重给出23.5 km²/s²,质量权重给出80 km²/s²,差别来自权重本身。

第二,反相残差也需要检验它是否能唯一指向扭转。设两个沿视线重叠成分的速度分别为$u_1(t),u_2(t)$,中心像素的权重各半,两侧分别是0.8/0.2和0.2/0.8。各侧减去中心后,就得到$\pm0.3(u_1-u_2)$:两边自动等幅、反相、同时过零。9

这个例子没有拟合真实日冕几何,不能据此断言论文信号是伪影。它说明这些关系需要与独立的结构、空间分辨率和混合条件共同解释。原研究也承认残差减除可能留下整体摆动成分,管外信号尤其要谨慎;现有固定狭缝数据还不能直接测出该信号的传播速度。12

振幅修正会显著改变能流估计

论文报告的视线方向整体摆动振幅约0.404 km/s,天空平面振幅约24.8 km/s;扭转相关视线振幅约0.321 km/s。以上振幅均按$\sqrt2$乘均方根的约定表示,实际运动是宽频的,通常只能找到少数相近周期的波包。1

若假定整体摆动的约61倍视线压低因子也适用于扭转信号,本库用这些中心值反算出约19.7 km/s的等效振幅,与作者约19.5 km/s的估计接近。但关键是假定能够共用这个修正因子,不能把修正后的数值称为直接读数。

再取自选密度$10^{-12}$ kg/m³和Alfvén速度1,000 km/s,按均方根约定计算,这个修正振幅对应约194 W/m²;未经修正的视线振幅只给出约0.052 W/m²。这个对照没有测定真实加热功率,它说明速度修正一旦进入平方,几何与权重条件就非常重要。9

作者结合其波模态推断和既往整体摆动波估计,给出约100—400 W/m²的合计能流范围,认为可能满足安静日冕和快太阳风的能量需求。1更稳妥的理解是:波动的能量规模值得认真追踪,但这些波最终在哪一层、以多大份额交给粒子,还需要下一段证据。

反射为什么可能让波动产生更小尺度

太阳大气的密度、磁场和Alfvén速度随位置变化,向外传播的波可能产生反射。在标准的非压缩MHD图景里,不同传播方向的扰动相互作用,使能量向更小的横向尺度转移。Chandran和Hollweg的分析、以及AWSoM模型,都使用了与返回波强度有关的级联闭合。34

本库用一个有限区域的能量盒说明这一层关系。$w_+$、$w_-$分别表示两个方向的波能;区域接受固定输入$P$,两种波都能离开区域,内部的反射交换项为$R(w_+-w_-)$。另以一个自选系数$C$表示向未解析小尺度的转移:

$$ \begin{aligned} Q_+&=Cw_+\sqrt{w_-},\qquad Q_-=Cw_-\sqrt{w_+},\\ \frac{\mathrm d(w_++w_-)}{\mathrm dt} &=P-\frac{w_++w_-}{\tau}-Q_+-Q_-. \end{aligned} $$

反射交换在总账中相消,自己不创造能量。这里将波的相位平均掉了,并没有复现真实反射方程或整个太阳的几何;$Q$也只是这个闭合中离开大尺度波的功率。

令$P=1,\tau=1,C=4$,在稳态比较三种区域内交换率,得到:

区内交换率 转入小尺度的功率份额 仍离开该区域的份额
0 约0 约1
0.1 0.354 0.646
1 0.575 0.425

没有持续产生返回成分时,最初很小的逆向扰动逐渐消失,这一特定级联通道便趋于零。加入内部交换后,一部分能量更容易转到小尺度。9

本库能量盒。左图虚线为逆向成分,中图为转入小尺度的功率,右图仅比较稳态份额。模型起初已有储存波能,所以过渡阶段的转移功率可以短时超过输入;累计能量仍闭合。
本库能量盒。左图虚线为逆向成分,中图为转入小尺度的功率,右图仅比较稳态份额。模型起初已有储存波能,所以过渡阶段的转移功率可以短时超过输入;累计能量仍闭合。

这不意味着真实日冕中所有加热都必须靠这种反射通道。非均匀结构、模态转换、局部磁能释放和粒子共振还可能参与。例子固定的是已经进入目标区域的能量;更低层的反射也可能让更少能量抵达该区域,不能由这张表推出“反射越强,整个日冕必然越热”。

到了粒子尺度,磁场怎样把能量交给粒子

单个带电粒子的运动满足Lorentz力。对它的动能求导:

$$ \frac{\mathrm d}{\mathrm dt} \left(\frac12mv^2\right)=q\,\mathbf E\cdot\mathbf v. $$

磁力改变运动方向,本身不直接改变粒子的动能;磁场结构和波动产生、组织的电场却可以做功。不同加热机制的差别,落在电场与粒子运动怎样保持关系、以及这种关系是否能长期产生净能量转移。

回旋共振是一种相干路线。粒子绕磁场回旋,当它沿磁场看到的波频率与回旋运动满足相应共振关系时,连续几次作用可以在合适相位上积累。净交换方向还取决于粒子分布;有些不稳定分布会把粒子能量交还给波。因此,探测到回旋波不自动等于探测到净加热。

随机加热是另一种路线。回旋轨道尺度上的场变化足够强时,粒子原先近似保持的磁矩可能被破坏,多次非相干作用改变垂直随机运动。常用的低$\beta$随机加热模型含有类似:

$$ \mathcal H\propto \delta u^3\exp\!\left(-\frac{c_2v_\perp}{\delta u}\right) $$

的非线性因子,具体还包含回旋尺度、频率和分布函数等项。这里$\delta u$代表相关尺度的扰动速度,$v_\perp$是粒子垂直于背景磁场的速度。5

这使“先把振幅平均,再算加热”与“先算每次作用,再平均”可能很不一样。本库保持$\delta u$的均方根为$0.1v_\perp$:一种情况始终取该振幅;另一种95%的时间为零,5%的时间为$0.447v_\perp$。取$c_2=0.34$,仅比较上述模型因子,后者平均值约为前者的62.7倍。9

它不说明任何尖峰都必然高效加热,而是说明在带有阈值抑制的机制里,均方振幅不足以指定结果。事件的大小分布、所处尺度、相位及模型适用条件都重要。

一个能够关上能量账的速度扩散例子

经过适当统计化,有些波粒作用可以用速度空间扩散描述。注意这里扩散的是速度分布,不是气体在空间中从一处流向另一处。对于已经积分掉平行速度、具有回旋对称性的简化分布$F(v_\perp)$,可写成:

$$ \frac{\partial F}{\partial t} =\frac1{v_\perp}\frac{\partial}{\partial v_\perp} \left(v_\perp D\,\frac{\partial F}{\partial v_\perp}\right). $$

若边界没有粒子通量,粒子数积分中的权重是$2\pi v_\perp$。相应垂直随机运动的能量变化率与:

$$ Q_\perp=-2\pi m \int_0^\infty v_\perp^2D \frac{\partial F}{\partial v_\perp}\,\mathrm dv_\perp $$

有关。由此可见,只有扩散系数还不够,还要看它与分布梯度在什么速度位置重叠。5

本库指定一个在典型垂直速度附近起作用的正扩散系数,并让它随一个有限波能库的剩余能量减弱;粒子获得多少随机运动能,就从波能库扣除多少。计算中,初始约0.500单位粒子能量和0.5单位波能,逐渐变为约1.000单位粒子能量与接近零的波能;粒子数与总能量守恒。9

本库有效模型。左上显示速度分布变化,右上核对波能和粒子能;左下的熵增属于采用扩散闭合后的分布。右下是相同均方振幅下,非零大振幅事件的时间份额改变时,随机加热模型因子的变化。
本库有效模型。左上显示速度分布变化,右上核对波能和粒子能;左下的熵增属于采用扩散闭合后的分布。右下是相同均方振幅下,非零大振幅事件的时间份额改变时,随机加热模型因子的变化。

随机运动能增加,不保证分布已经变成新的麦克斯韦分布。这个扩散模型还预先包含了相位平均和统计化,不能拿它的熵增去证明精细、无碰撞、封闭的Vlasov动力学也必然产生同样熵增。热流、分布的非平衡结构和所用熵定义都需要分别说明。7

扩散描述本身也不是随机加热的专属标签;共振的准线性理论同样可以得到扩散方程。要区分机制,还需比较扩散和能量交换在速度空间的形状、强度及其对等离子体条件的依赖。

近太阳探测器正在怎样检验这一段机制

2025年和2026年的两项工作研究了Parker Solar Probe的一段亚Alfvénic太阳风,比较平行传播的离子回旋波共振与随机加热。2026年研究明确采用2021年4月28日11—13时的窗口,估计其磁连接相对稳定,来自同一负极性日冕洞;这段流的热压相对磁压较低,两个方向的Alfvénic波动能量又明显不平衡。56

2025年的比较发现,纳入大振幅间歇事件后,某种随机加热估计可达到局部级联率的量级,而平行回旋波的净作用不足。2026年的工作进一步把质子分布拟合到速度网格,利用沿途变化和导引中心方程反推扩散系数,再比较不同理论的速度空间加热形状。

Ervin等2026年论文图3的原始比较图,未经本库重拟。粉色为作者的经验反演,实蓝线为计入间歇性的模型;蓝虚线与紫线均在图中放大了1000倍,所有曲线又按经验峰值归一化。因此不能把图上相近的线高当成相近的原始加热率。
Ervin等2026年论文图3的原始比较图,未经本库重拟。粉色为作者的经验反演,实蓝线为计入间歇性的模型;蓝虚线与紫线均在图中放大了1000倍,所有曲线又按经验峰值归一化。因此不能把图上相近的线高当成相近的原始加热率。

图源:Ervin等,Physical Review Letters 137, 065201,Fig. 3;CC BY 4.0。从正式PDF提取原图,未裁切或改动数据。6

这是有价值的机制检验路线,但“经验反演”仍有条件:需要近似稳定、回旋对称的分布及指定的磁通管面积变化,需处理仪器视场与插值,还需把不同位置测量间的变化解释为空间演化。116个滚动窗口估计也不能当成116次独立实验。这两篇研究共享相关流事件,属于互补分析。

本库还保留一处未闭合的方法问题:一般二维情况下,要先保留扩散系数与分布梯度的联合结构,再对未显示的速度方向积分;不能无条件把两者分别积分后再相乘。该研究从二维系数到一维加热率的书面步骤需要进一步的定义或闭合说明。完整反演实现和处理后系数未取得,因此本文没有独立确认其归一化曲线或绝对加热率,也不根据这一点断言作者实际程序一定算错。

更重要的是,这个近太阳风时段与DKIST的低日冕狭缝,并不是同一次沿流测量。即使某种机制能解释前者,也不能直接拿来填写后者的全部加热预算。

怎样理解现阶段已经知道的部分

已有观测让我们看到,低日冕存在丰富的横向和扭转相关运动;理论与全局模型说明,反射、级联和磁能释放能够把低层输入的能量重新分配。粒子理论和原位研究则给出了共振、随机作用及其他通道的具体检验方式。

全局模型仍需要闭合这些环节。AWSoM把波传播、反射、耗散、电子与离子能量分配、热传导和辐射放在同一系统中计算,能够与多波段辐射和太阳风性质比较;但其中的级联长度和粒子分配处方仍有模型输入,合成图像吻合不等于每一笔微观加热都被独立测量。4

间歇磁重联也不必与波动形成二选一。它可以直接改变粒子和流动的能量,也可以先激发波;波动产生的小尺度结构又可能促进局部磁能释放。区分它们应看具体区域、时间尺度与能量收支,而不是给所有安静区、活动区和太阳风时段指定同一种主导过程。

对本篇主问题,能够清楚回答的是:太阳低层的机械功和电磁能可以跨越温度梯度,把能量送到稀薄外层;在合适尺度和等离子体条件下,场与粒子的相互作用使随机运动能增加。真正尚未完全确定的,是不同区域中有多少输入经过哪条路线、在哪里沉积,以及这些沉积怎样共同抵偿当地的辐射、热流和外流损失。

来源与依据

  1. S1 原始DKIST观测与波模态建模

    Evidence for small-scale torsional Alfvén waves in the solar corona

    R. J. Morton; Y. Gao; E. Tajfirouze; H. Tian; T. Van Doorsselaere; T. A. Schad · 2025

    在线2025-10-24,NatAstron10,42–53(2026),正式OAHTML · 相关全文已核

    正式Main/Results、Methods数据处理/非热线宽/残差/MonteCarlo与其他来源段

    狭缝位于约.1Rsun日面外,跨开放与安静环及其他结构;无法直接测传播速度。低通削弱<2Mm/<30s。振幅按sqrt2*RMS定义,宽频信号不是实测峰值;61倍修正与100–400W/m²合计是条件性推断,不当局部热测量。

  2. S2 原始误差、视线权重与波模型补充

    Supplementary information: small-scale torsional Alfvén waves

    R. J. Morton and coauthors · 2025

    正式MOESM1;实际读取上述相关段落 · 相关全文已核

    §1–5误差/滤波/局域性/发射率加权,§6波剖面与MC相关说明

    Gaussian谱线一阶矩为发射强度加权;不自动成为质量权重。误差估计有像素差与单谱Bayes不同口径,条纹主导部分残差;不自称本库复做原始谱线拟合/EMD/PLUTO。

  3. S3 原始反射与级联近似理论

    Alfvén Wave Reflection and Turbulent Heating in the Solar Wind from 1 Solar Radius to 1 AU: An Analytical Treatment

    Benjamin D. G. Chandran; Joseph V. Hollweg · 2009

    arXiv0911.1068v1,2009-11-05;PDF编译页眉2018不作发表年 · 相关全文已核

    导言、Heinemann–Olbert框架及Eq20–33附近非线性闭合/长周期假设

    强不平衡、指定背景与横向相关长度下的近似,不能直接套全部观测频段。本文能量盒改用守恒的非相干交换,只说明记账与逆向成分作用,不冒充其解析解复现。

  4. S4 原始全局日冕—太阳风模型作者稿

    Alfvén Wave Solar Model (AWSoM): Coronal Heating

    B. van der Holst; I. V. Sokolov; X. Meng; M. Jin; W. B. Manchester IV; G. Tóth; T. I. Gombosi · 2013

    实际arXiv1311.4093v3,2013-12-03;对应后续AWSoM研究 · 相关全文已核

    摘要、热传导/电子离子能量分配及Eq11–13波能反射/耗散

    反射项总能量抵消,Gamma±依返回波与Lperp;粒子分配和传导有闭合选择。未复跑全3DAWSoM或合成EUV,匹配观测不当每项微观机制独立验证。

  5. S5 原始PSP加热估计与模型比较

    Stochastic Heating in the Sub-Alfvénic Solar Wind

    Trevor A. Bowen; Tamar Ervin; Alfred Mallet; Benjamin D. G. Chandran; Nikos Sioulas; Philip A. Isenberg; Stuart D. Bale; Jonathan Squire; Kristopher G. Klein; Oreste Pezzi · 2025

    PRL135,255201,2025-12-18正式PDF;CCBY4.0 · 相关全文已核

    Main、Eq1–3随机加热/分布导数与间歇性、EndMatter相关假设

    常用δu³exp(-c2/epsilon)参数依模型与标定;RMS后计算不等于逐时刻计算后平均。局部LET与加热比较不是全日冕预算。本文未重算SCaM时间序列或相干波功率;原文April8与图注April28冲突记录,不编造独立事件。

  6. S6 原始PSP速度分布反演研究

    Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas

    Tamar Ervin; Trevor A. Bowen; Alfred Mallet; Philip A. Isenberg; Kristopher G. Klein; Stuart D. Bale; Benjamin D. G. Chandran; Roberto Livi; Ali Rahmati; Davin E. Larson · 2026

    PRL137,065201,2026-08-06正式PDF;arXiv2607.11880PDF及源归档另存;CCBY4.0 · 相关全文已核

    Main Eq1–3、Fig1–3、Discussion/EndMatter;作者源main.tex中的二维到一维写法

    Apr28,2021,11–13UTC,低β/强不平衡/稳定连接假设;9点约束、116滚动窗口,不当独立样本。完整反演SI/处理后D和实现未取得,书面D积分再乘独立积分梯度需要额外定义;不据此判实际程序必错,也不独立确认归一化曲线或绝对率。原图两条线×1000已回图核。

  7. S7 原始无碰撞能量与熵理论

    Energy dissipation and entropy in collisionless plasma

    Senbei Du; Gary P. Zank; Xiaocan Li; Fan Guo · 2020

    PRE101,033208,2020-03-24正式PDF · 相关全文已核

    摘要/Introduction及Kinetic and Fluid Entropy相关起始定义

    分清温度、随机能量、流体熵代理与细粒度动力学熵;本文不重跑其PIC。自有FokkerPlanck已作统计闭合,不能用其熵增证明封闭Vlasov细粒度熵增。

  8. S8 观测目录与读取范围

    DKIST data catalogue and linked analysis-code provenance

    NSO DKIST Data Center; Tom Schad; Richard J. Morton · 2026

    目录核对2026-09-13;代码树cryonirsp0951eba与NUWT24a2f7c存档 · 题录已核

    2023-10-30同次计划的公开目录及链接GitHub树/Readme

    当前JAXBIO为CI图像,无光谱轴;同期LHWWPP有光谱轴,5.914GiB/4024frames,均公开可下载。所链示例notebook为BZJOM另一数据集。只取得目录与工具信息,未取科学帧;不把目录绿灯或别的数据当本例复现,不称数据永久不可访问。

  9. S9 本库简化模型与数值核对

    T84 — coronal readout, wave-energy budgets and effective velocity diffusion

    Codex · 2026

    2026-09-13本库计算;实际自审,不称独立审阅 · 相关全文已核

    readout/wave_box/diffusion_model/particles_and_intermittency,results.json

    发射权重与谱矩、残差混合反例、RMS振幅能流换算、守恒能量盒、回旋对称速度扩散与间歇模型;不称原始DKIST/PSP复算。网格核全时间能量差和概率L1,不仅核预定终态总能量。

    本地原件
最后更新 2026-09-13

来源