假如只能看到一群粒子造成的加热,却不能直接看到粒子本身,能否从“哪里被加热了”推断它们的能量?2017 年 8 月 21 日那场 A 级微耀斑,为这个问题提供了一个难得的检验机会。

在上篇中,我们介绍了如何用 DEM 约束热辐射,再次确认 NuSTAR 能谱中的非热成分。拟合给出的电子低能截止约为 7.1±0.57.1\pm0.5 keV。与此同时,图像显示两个超过 10 MK 的热源最初出现在日冕中,明显高于色球足点。

这两条线索指向一个共同的问题:由能谱得到的电子参数,是否能解释热源出现的位置及其随时间的演化? 在 Li 等(2022) 的研究中,我们据此首次给出了非热电子低能截止的自洽成像证据:能谱、热源位置和电子输运计算相互吻合。这种约束仍依赖物理模型,并不是对电子截止能量的直接成像测量。

1. 低能截止:一个小参数,牵动很大的能量估算

太阳耀斑中被加速的非热电子,常用幂律分布描述:能量越高,电子越少。如果把这条分布一直延伸到任意低的能量,电子总数和总能量就可能变得不合理。因此,模型需要规定它在低能端如何结束或过渡。最常用的一种近似就是引入低能截止 EcE_c(Holman 等,2011)。

在这里采用的截断幂律模型中,EcE_c 是非热幂律分布的下限。它不是说太阳大气里不存在低于这个能量的电子——低能热电子大量存在;它也不是探测器的最低能道,或热、非热光子谱相交的位置。真实的非热分布向热背景过渡,可能比理想化的“突然截断”更复杂。

用 F(E)=AE−δF(E)=AE^{-\delta} 表示单位电子能量区间内的注入率,且 E≥EcE\geq E_c。当 δ>2\delta>2 时,非热电子功率满足

Pnth=∫Ec∞EF(E) dE=Aδ−2Ec 2−δ.P_{\rm nth}=\int_{E_c}^{\infty}E F(E)\,dE =\frac{A}{\delta-2}E_c^{\,2-\delta}.

这个式子要求能量单位一致;若 EE 使用 keV,积分结果转为 erg/s 时还需要乘以 1.602×10−91.602\times10^{-9} erg/keV。更关键的是它揭示的趋势:在高能端归一化 AA 和谱指数固定时,低能截止越低,推算的非热能量越大。 对很陡的电子谱,这种变化尤其显著。

本事件的电子谱指数约为 6.8,因此电子数和能量输入都主要集中在低能端附近。这不是一束所有电子都具有 7.1 keV 的单能电子,而是一个以低能端占主导、仍包含较高能电子的连续分布。确定 EcE_c,也就成为确定能量输入的重要环节。

2. 为什么不能只靠能谱

困难在于,低能截止影响最明显的能段,往往也是热辐射很强的地方。改变热等离子体的温度分布、背景,或非热分布的低能端,都可能改变同一段观测能谱(Holman 等,2011)。同一事件的早期 NuSTAR 能谱分析见 Glesener 等(2020)。

上篇使用的 DEM 方法,先从 AIA 多通道图像约束热成分,再通过 NuSTAR 的响应计算热计数谱,最后拟合剩余的非热成分。这样得到约 7.1 keV 的截止值,比完全依靠任意热模型的拟合有了更强的观测约束。

但研究没有停在这里。既然这些电子被认为参与了加热,就应该进一步检查:它们在观测到的日冕环境中,会在哪里损失能量?计算的位置,是否与真实的热源相符?

3. 两个日冕热源,对电子输运提出了什么问题

在通常的图景中,电子从日冕中的加速区域沿磁场运动,较高能电子可以穿过稀薄的日冕,到达较致密的下层大气并沉积能量。这能够形成足点辐射和色球蒸发。

然而,本事件最早出现的两个 >10 MK 热源,明显离开了色球足点。研究考虑了两种解释:

  1. 局部高密度厚靶。 两处局部密度显著增强,像两个集中阻挡电子的区域,使一段较宽能量范围内的电子都主要停在相近的位置。
  2. 柱密度与低能截止共同决定热化位置。 电子沿途不断碰撞、损失能量,不同初始能量的电子停在不同深度;由于谱很陡,大多数电子靠近低能截止,它们集中在相近高度热化,形成两个高温源。

要区分这两种解释,需要由沿环的密度分布计算电子停止位置,再与高温源的位置比较。

4. 从 DEM 到柱密度,再到电子停止位置

DEM 提供了不同温区的 EM。若视线方向的有效深度为 ll,采用均匀密度与单位填充因子的近似,电子密度可以估算为

ne≃EMcoll.n_e\simeq\sqrt{\frac{\mathrm{EM}_{\rm col}}{l}}.

其中 EMcol\mathrm{EM}_{\rm col} 是沿视线积分的柱发射度,单位为 cm−5\mathrm{cm^{-5}}。论文在估算耀斑环的密度前,扣除了宁静日冕背景的 EM,并将环宽约 5.4′′5.4'' 作为视线深度的近似。

电子沿磁环前进时,决定碰撞损失的并不只是某一点的密度,而是它沿途遇到的累积柱密度:

N(p)=∫injectionpne(s) ds.N(p)=\int_{\mathrm{injection}}^{p} n_e(s)\,ds.

它可以理解为:从注入位置走到 pp,电子总共穿过了多少靶粒子。注意它的单位是 cm−2\mathrm{cm^{-2}},与 DEM 所给出的 EM 不是同一个量。

在非相对论碰撞损失模型中,电子的能量沉积与沿途累积柱密度有关(Emslie,1978)。采用 Li 等(2022) 文中的靶区条件,使初始能量为 EsE_s 的电子停止所需的柱密度约为

Ns≃1.5×1017(EskeV)2cm−2.N_s\simeq1.5\times10^{17} \left(\frac{E_s}{\mathrm{keV}}\right)^2\mathrm{cm^{-2}}.

这里“停止”指电子损失大部分定向非热能量、逐渐热化,不是电子消失。这个关系也不是适用于任何太阳大气环境的万能公式;它依赖所采用的碰撞输运和靶区条件,不能替代对返回电流等效应的完整处理。

代入便可看出差别:5 keV 电子需要约 3.75×1018 cm−23.75\times10^{18}\ \mathrm{cm^{-2}},7.1 keV 电子需要约 7.56×1018 cm−27.56\times10^{18}\ \mathrm{cm^{-2}}。后者大约是前者的两倍。这两项是按上式计算的示例值。给定沿环的密度分布,就能寻找累积柱密度达到这些数值的位置。

5. 关键检验:7.1 keV 的位置与热源吻合吗

研究先假定电子从环的中部注入,利用各个时刻的 DEM 密度分布,分别计算 5 keV 和 7.1 keV 电子向两侧传播的停止位置。

沿环密度分布与时间—距离图,7.1 keV电子停止位置的黑线与蓝色高温EM源接近,5 keV的灰线则明显偏离
Li 等(2022)图 5(b、c):左图显示沿环密度随时间演化;右图将计算的电子停止位置与 >10 MK 的 EM 时间—距离图叠加。黑线对应 7.1 keV,灰线对应 5 keV;蓝色背景是高温 EM。右图中其他彩色曲线表示所标热源位置的密度演化,不能与停止位置混为一谈。原始论文,CC BY 4.0。

检验有三个关键结果。

不同初始能量的电子,确实停在不同位置。 5 keV 和 7.1 keV 的计算轨迹明显分开。这不支持“只由两个局部密度峰把不同能量的电子都挡在同一个高度”的简单解释。

7.1 keV 的停止位置与最初两个热源相符。 在 18:50:42 UT 两个 >10 MK 热源出现时,计算位置与其位置吻合,而 5 keV 的结果更靠近注入区域,与观测不符。这个 7.1 keV 并不是为了让图中的黑线经过蓝色源才随意挑选的,它来自前面的非热能谱拟合。

吻合还体现在随时间的演化上。 随着环内密度增加,电子走更短的距离就能遇到足够的柱密度,因此停止位置逐渐向环中部靠近。计算出的趋势与两个高温源向中部延伸、汇合的过程一致。

因此,证据不仅是“两点在某一刻重合”,而是能谱参数、沿环密度、热源的空间位置及其时间演化,共同支持了同一套电子输运解释。

6. “自洽的成像证据”究竟说明了什么

这一物理图景可以这样理解:陡峭的非热电子谱使大量电子集中在低能截止附近;由于初始能量较低,它们不必到达色球,就可以在日冕中经过碰撞热化,形成两处高温源。较高能的电子仍可能继续进入下层大气并驱动蒸发。热源的后续扩展还可能涉及热传导和蒸发等过程。

这里必须区分电子能量与等离子体温度。7.1 keV 描述的是注入非热电子的低能截止,>10 MK 描述的是被加热等离子体的温区。 图像不是直接拍到了“一束 7.1 keV 电子”;它拍到的是加热后的辐射结构,再由密度和输运模型建立两者的联系。

“自洽”也不表示每一项测量都完全独立。能谱分析中的热辐射约束和输运计算中的密度都使用了 DEM,因此共享部分仪器响应、温度反演和几何假设。但它们检验的是不同观测特征:一个检验 X 射线谱,另一个检验热源的位置与演化。两种约束吻合,为低能截止提供了单纯能谱拟合之外的空间证据。

7. 换一个几何假设,结论还站得住吗

DEM 的误差会传到密度,视线深度会改变密度的绝对值,电子注入位置也会影响沿途柱密度。论文附录专门检查了这些问题:在观测强度上加入噪声、重复 100 次 DEM 计算估计误差,并传播至柱密度和停止位置;还考虑了相对于约 5.4′′5.4'' 视线深度而言已经较大的深度偏差。

在所检验的环中部注入情形下,即使取 3′′3'' 的视线深度不确定性,计算停止位置与热源的吻合仍然保持。作为对照,改变为明显偏离环中部的注入位置,所得轨迹难以同时匹配两侧高温结构。这说明几何假设会受到观测约束,并非任意位置都可以解释数据。

查看论文中的误差与注入位置检验
论文图6比较环中部与偏离中部的注入位置,展示DEM、视线深度和低能截止误差对停止位置的影响
Li 等(2022)图 6。上排采用环中部注入,并估计视线深度不确定性为 3″;下排采用另一注入位置及 2″ 的深度不确定性。灰色误差区及误差棒显示传播后的不确定性。两排用于检验具体假设,不代表已穷尽所有可能的输运模型。原始论文,CC BY 4.0。

仍需保留边界:AIA 图像的饱和、配准与响应误差,视线叠加、填充因子、环的三维几何,以及返回电流或其他输运过程,都可能影响定量结果。论文展示的是在本事件中,经多种观测和误差检验支持的一套自洽解释,而不是彻底消除了低能截止的模型依赖性。

8. 从一个参数,走向加热发生在哪里

这项研究的意义不止是把低能截止写成一个带误差的数字。它把一个难以从能谱单独确定的参数,转化成了可以用图像检查的问题:哪些电子能够走到色球,哪些会在日冕中热化,最热的物质应该出现在哪里?

对于这场微耀斑,热源的位置提示低能电子的热化可以直接参与日冕加热;较高能电子及色球蒸发仍然发挥作用,但不需要把所有高温物质都解释为先在色球被加热、再整体上流的结果。不同能量的电子,其能量沉积位置可能不同。

这为研究其他微耀斑提供了一条值得检验的方法:先尽可能准确地约束热成分,再结合空间结构和密度演化检验电子输运。至于其他事件是否具有同样的低能截止和热源结构,需要各自的观测证据,而不是直接套用本事件的数值。

回到上篇:非热成分、能量分配与色球蒸发。方法基础:DEM 简介与使用方法。

参考文献