Optimization Algorithm for Determining the Source Surface Radius Based on Parker Solar Probe in situ Measurements from Encounters 1 to 19

技术摘要：基于帕克太阳探测器原位观测结果确定源面半径的优化算法

问题陈述
势场源面（PFSS）外推法是利用光球磁图估算大规模日冕磁场的标准方法。该模型中的一个关键参数是源面半径（$R_{ss}$），即磁场被假定为径向的球面边界。虽然 $R_{ss} = 2.5 R_s$ 是一个常规参考值，但先前的研究表明，单一数值无法普遍满足所有观测约束（例如：日冕空洞范围、开磁通量和磁极性）。已知最优的 $R_{ss}$ 会随太阳活动而变化，然而，如何确定一个严谨的、受观测约束的数值仍然具有挑战性。本文通过最小化 PFSS 外推结果与原位磁场测量值之间的差异，来解决确定最优 $R_{ss}$ 的逆问题，特别利用了帕克太阳探测器（PSP）在第 1 至 19 次遭遇期间的数据。

方法论
作者将 $R_{ss}$ 的确定视为自由边界框架下的形状优化问题。

1. 数学公式化：

   • 正向 PFSS 问题： 被公式化为一个带有混合狄利克雷（Dirichlet）和诺伊曼（Neumann）条件的二阶椭圆边界值问题。作者证明了正向解在 Sobolev 空间 $H^1(\Omega)$ 中的适定性（存在性与唯一性）。
   • 逆问题： 被构建为在可容许定义域集合上最小化目标泛函。作者证明了可容许集合的紧致性和目标泛函的连续性，从而确立了逆问题全局极小值的存在性。这些证明同时适用于球形和非球形的源面几何结构。

2. 数据与预处理：

   • 输入： 来自全球振荡网络组（GONG）的综合光球磁图作为内边界条件。
   • 观测： PSP（FIELDS 和 SWEAP 仪器）在第 1 至 19 次遭遇期间测得的径向磁场（$B_r$）和径向太阳风速度（$v_r$）。
   • 处理： 对 PSP 数据进行异常值剔除（采用四分位距准则）和时间平均处理（20 分钟步长）。通过线性插值填充太阳风速度间隙，以实现连续的帕克螺旋线反向映射。
   • 验证： 使用来自 1 AU 附近先进太阳观测者号（ACE）的独立原位测量数据进行交叉验证。

3. 优化算法：

   • 目标泛函： 算法旨在最小化观测到的 PSP $B_r$ 与 PFSS 预测的 $B_r$ 之间的均方误差（MSE）。PFSS 预测是通过沿帕克螺旋线将航天器位置回溯至源面，采样 PFSS 场，并将其径向缩放到航天器距离（$r^2 B_r = \text{const}$）来获得的。
   • 归一化： 为了消除不同遭遇期间的幅度差异，在计算 MSE 之前，观测场和模拟场均经过 Z-score 归一化处理。
   • 求解器： 采用基于梯度的迭代方案（算法 1）。梯度使用中心有限差分进行近似。算法迭代直至收敛或达到最大迭代次数。
   • 验证： 该算法通过针对解析势场解和添加了高斯噪声的合成数据进行了验证，证明了其恢复预设 $R_{ss}$ 的能力。

4. 多目标分析：

   • 为了理解权衡关系，作者引入了两个额外的指标：磁场缩放因子（$\alpha_k$，与开磁通量相关）和极性预测准确度（$P^k_a$）。
   • 进行帕累托分析（算法 2）以识别 MSE 目标与极性预测误差之间的帕累托前沿（Pareto frontier）。这量化了优化主要是由幅度一致性还是极性准确性所驱动。

主要结果

• 最优 $R_{ss}$ 趋势： PSP 第 1 至 19 次遭遇的最优 $R_{ss}$ 值总体上呈现从太阳极小期向第 25 太阳周上升期的增加趋势。这一趋势得到了独立 ACE 数据的证实。
• 与参考值的比较： 与固定的参考值 $R_{ss} = 2.5 R_s$ 相比，优化后的解通常表现出与观测磁场幅度更好的符合度（减少了开磁通量差异），同时保持或提高了磁极性预测的准确度。
• 太阳极小期 vs. 上升期：
  • 在太阳极小期（早期遭遇），优化主要由幅度一致性（开磁通量）驱动。选择较低的 $R_{ss}$ 值以增加模拟的开磁通量，从而匹配观测值。
  • 在上升期，帕累托前沿变得更加聚集，表明极性预测准确度成为了主导指标。最优 $R_{ss}$ 随之增加，以更好地捕捉观测到的复杂扇区结构和极性反转。
• 局部偏差： 特定遭遇（如第 9 次和第 13 次）显示出偏离总体趋势的局部现象。对于第 13 次遭遇，由于穿越了复杂的行星际日冕物质抛射（ICME），优化选择了极低的 $R_{ss}$（$1.27 R_s$）以捕捉复杂的极性反转，这凸显了在扰动环境下帕克螺旋线假设的局限性。

意义与主张
本文声称为 PFSS 逆问题提供了数学上严谨的基础，证明了源面半径估计中唯一极小值的存在性。通过构建一个经过解析解和原位数据验证的特定优化算法，本研究表明 $R_{ss}$ 不是一个静态常数，而是一个随太阳活动演化的动态参数。

作者强调，不存在单一的“全局最优”$R_{ss}$；相反，最优值取决于特定的观测约束和优先考虑的评估指标（例如：开磁通量 vs. 极性）。研究指出，主导指标从太阳极小期的幅度一致性转向了上升期的极性一致性。此外，论文承认了“开磁通量问题”，指出虽然优化改善了与原位磁通量的符合度，但并未解决日冕空洞磁通量估算与原位测量之间的根本差异，这表明未来的改进可能需要更好的极区磁场观测或更先进的日冕模型。这项工作为未来利用来自太阳轨道器（Solar Orbiter）等任务的非球形源面形状优化奠定了框架。