[论文解读] A New Multi-Constraint Potential Field Source Surface (PFSS) Extrapolation Model

C. Antonio, I. Chifu+2 more|arXiv (Cornell University)|Mar 20, 2026

Solar and Space Plasma Dynamics被引用 0

一句话总结

本文提出一种多约束PFSS外推方法，将3D日冕环几何引入PFSS框架，在保持高效性的同时提高与环结构的拟合并控制散度与（可选的）力自由性。

ABSTRACT

The Potential Field Source Surface (PFSS) model is the most used approach for extrapolating the global coronal magnetic field, offering efficiency and strong performance at large scales. However, PFSS assumes a potential coronal field, so it cannot account for distortions from electric currents. More advanced methods, such as nonlinear force-free field (NLFFF) models, can represent these effects but are much more computationally intensive. Recent observational techniques also allow 3D reconstruction of coronal loops, which trace solar magnetic field geometry. This work develops a new approach that constrains the PFSS model using 3D coronal loop information, improving agreement with observations while keeping efficiency. The model is based on PFSS field constraints from photospheric data but allows magnetic field deviations from the potential state within loop-influenced regions, maintaining control over divergence and force-freeness. We adapted NLFFF optimization to the PFSS framework, enabling multiple physical constraints. Our functional includes up to three terms: divergence-free, loop geometry, and force-free. The resulting Python algorithm was tested with synthetic loops, using Carrington rotation 2284 as the lower boundary. This method yields magnetic field solutions that better match the geometry of included loops and controls divergence and force-freeness. Our results show that 3D coronal loop information can be incorporated into PFSS, largely preserving computational efficiency even with many loops. This approach lets PFSS better reflect observed coronal structures without significant computational cost.

研究动机与目标

• 将3D日冕环几何引入PFSS外推。
• 在保持PFSS边界条件的同时，允许环附近偏离势场状态。
• 确保磁场保持散度为零并评估力自由性的潜在提升。
• 提供一个灵活的优化框架，可以开关约束以提高效率与探索性。

提出的方法

• 通过定义总泛函Lt = L1 + L2 + ... + Ln，将散度、环切触性和（可选的）力自由项纳入NLFFF优化形式，改造为PFSS。
• 使用二阶有限差分的Python算法，通过Landweber型更新以自适应步长μn迭代最小化Lt。
• 将磁场插值到3D日冕环以计算环约束L2，并再插回网格，确保一致性和最小散度。
• 约束下边界匹配光球面PFSS输入，同时保持源表面和底边界固定。
• 允许开关L2（环切触）和L3（力自由）约束，以研究它们对收敛、切触和全局可压性的影响。

Figure 1: GONG radial magnetic field CR 2284 synoptic map. The map is used as a lower boundary condition for the pfsspy .

实验结果

研究问题

• RQ1 3D日冕环几何是否能够在不牺牲全局计算效率的前提下，稳定地整合进PFSS外推？
• RQ2 将环切触约束引入后，结果磁场的散度和力自由性有何变化？
• RQ3 在包含NLFFF-inspired项时，环对齐与传统PFSS约束（散度与力自由性）之间的权衡是什么？
• RQ4 切换环切触约束开启/关闭如何影响收敛性和最终场拓扑？
• RQ5 方法对正则化参数ξ2和环数量的敏感性如何？

主要发现

• 多约束PFSS可以产生与3D日冕环更好对齐的磁场解，同时在大体上保持计算效率。
• 引入环切触约束在提高环几何拟合度的同时，可能带来更高的散度/较少的力自由性（在环包含阶段）。
• 在初始收敛后关闭环约束可降低附加散度，同时在一定程度上保留改进的环对齐。
• 两种优化设置表明，包含力自由约束（L3）可引导解向更高的力自由性和更低的全局散度，但对环切触的提升存在代价。
• 最优ξ2取值取决于约束数量与配置：设置I约为10，设置II约为约100。
• 在现代硬件上，方法以适度的迭代次数实现近似分钟级收敛。

Figure 2: Contour plot of the pfsspy output radial magnetic field at some height reference equal to r=1.15 R ⊙ .