[论文解读] An iterative method in a probabilistic approach to the spectral inverse problem: Differential emission measure from line spectra and broadband data

F. F. Goryaev, S. Parenti+4 more|arXiv (Cornell University)|Oct 25, 2010

Solar and Space Plasma Dynamics参考文献 26被引用 11

一句话总结

本文提出了一种贝叶斯迭代方法（BIM），用于从太阳极紫外和X射线光谱中重建微分发射度（DEM）分布，采用基于贝叶斯定理的概率方法避免正则化。该方法成功恢复了活动区中的多温态等离子体结构，确认在4–6 MK温度范围内存在高达2–4 MK温度范围30%的发射度，支持纳米耀斑加温范式。

ABSTRACT

Inverse problems are of great importance in astrophysics for deriving information about the physical characteristics of hot optically thin plasma sources from their EUV and X-ray spectra. We describe and test an iterative method developed within the framework of a probabilistic approach to the spectral inverse problem for determining the thermal structures of the emitting plasma. We also demonstrate applications of this method to both high resolution line spectra and broadband imaging data. Our so-called Bayesian iterative method (BIM) is an iterative procedure based on Bayes' theorem and is used to reconstruct differential emission measure (DEM) distributions. To demonstrate the abilities of the BIM, we performed various numerical tests and model simulations establishing its robustness and usefulness. We then applied the BIM to observable data for several active regions (AR) previously analyzed with other DEM diagnostic techniques: both SUMER/SOHO (Landi and Feldman, 2008) and SPIRIT/CORONAS-F (Shestov et al., 2010) line spectra data, and XRT/Hinode (Reale et al., 2009) broadband imaging data. The BIM results show that this method is an effective tool for determining the thermal structure of emitting plasma and can be successfully used for the DEM analysis of both line spectra and broadband imaging data. The BIM calculations correlate with recent studies confirming the existence of hot plasma in solar ARs. The BIM results also indicate that the coronal plasma may have the continuous distributions predicted by the nanoflare paradigm.

研究动机与目标

• 开发一种稳健的、无需正则化的光谱反问题求解方法，用于太阳物理学。
• 从高分辨率谱线光谱和宽带成像数据中重建微分发射度（DEM）分布。
• 将该方法与现有DEM分析进行对比验证，评估其检测太阳活动区中高温等离子体的能力。
• 研究与纳米耀斑加温范式一致的连续温度分布的存在性。
• 将BIM结果与Solarsoft中既有的滤波比法和迭代方法进行比较。

提出的方法

• 该方法通过迭代应用贝叶斯定理，更新温度和发射度的后验概率分布。
• 将观测通量和响应函数视为似然函数，结合先验知识以约束解的范围。
• 通过将反问题表述为概率形式而非弗雷德霍姆积分形式，避免了病态性问题。
• 通过迭代优化实现稳定DEM分布的收敛，无需引入正则化。
• 该方法在合成数据上进行了测试，并应用于SUMER/SOHO、SPIRIT/CORONAS-F和XRT/Hinode的真实观测数据。
• X射线滤波器和谱线的响应函数被用作贝叶斯更新过程中的似然核。

实验结果

研究问题

• RQ1 贝叶斯迭代方法是否能在光谱反问题中无需正则化即可重建可靠的DEM分布？
• RQ2 BIM在从谱线光谱和宽带成像数据中恢复多温态等离子体结构方面表现如何？
• RQ3 BIM是否能在非耀斑活动区中检测到T > 4 MK的高温等离子体，与纳米耀斑范式一致？
• RQ4 BIM结果与以往使用滤波比法或其他迭代技术的DEM分析在定量上如何比较？
• RQ5 当应用于具有宽响应函数的宽带数据时，BIM存在哪些局限性？

主要发现

• BIM与Landi & Feldman（2008）对SUMER/SOHO数据的结果在log T ≈ 6.5处表现出极佳的定量一致性，冷却峰值处差异仅为30–50%。
• 对于SPIRIT/CORONAS-F数据，BIM证实了在log T ≈ 6.8–7.2 MK温度范围内存在高温等离子体，其发射度与2–4 MK范围相当。
• 在XRT/Hinode宽带数据中，BIM在非耀斑活动区的一个子区域检测到4–6 MK温度范围的高温等离子体，其发射度最高可达2–4 MK范围的30%。
• 该方法揭示了活动区等离子体在log T ≈ 6.5–6.6处存在更宽的温度成分，与多温态结构一致。
• 仅在使用XRT滤波器时出现伪低温峰的人工伪影，这是由其宽温度响应函数引起，而非BIM本身固有的问题。
• BIM结果支持日冕中存在连续温度分布，与纳米耀斑加温模型一致。