REVIEW 3 major objections 4 minor 42 references
Magnetic model reveals why a 500 Mm solar filament erupted sideways
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
A neural-network magnetic field reconstruction shows that a 500 Mm filament's pre-eruption flux rope had an extended eastern footprint connected to the flare ribbon and coronal dimming, explaining the asymmetric eruption.
T0 review reviewed 2026-08-05 challenge →
load-bearing objection A careful single-event NLFFF study with a genuinely robust ~500 Mm flux-rope reconstruction; the dimming-expansion interpretation leans on the least stable connectivity in the ensemble, so it is plausible but not confirmed. the 3 major comments →
Magnetic structure and asymmetric eruption of a 500 Mm filament rooted in weak-field regions
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
Core claim
Using a PINN-based nonlinear force-free extrapolation of the pre-eruption photospheric field, the paper finds a channel of high current density whose length and shape match the observed 500 Mm filament. Field lines from the eastern portion of this channel form a magnetic flux rope with an extended, fanned-out footprint in a weak-field negative-polarity region; that footprint coincides with the area enclosed by the inverse J-shaped flare ribbon and with the initial coronal dimming. The same model shows overlying strapping field lines anchored in the region into which the dimming later expands, and a compact, multiply anchored western leg beneath strong sunspot-connected overlying fields. The
What carries the argument
The central object is a 500 Mm magnetic flux rope delivered by the NF2 physics-informed neural network nonlinear force-free extrapolation, in its vector-potential form. The neural network maps coordinates to a vector potential A, the field is B = curl A so divergence-free by construction, and training minimises a weighted sum of force-free, boundary, and potential-boundary losses; the lower boundary is the HMI vector magnetogram and the side and top boundaries are potential field. This mesh-free representation lets the method fill a roughly 730 x 550 x 300 Mm volume while allowing local departures from strict force-freeness in weak-field regions. The interpretation hinges on tracing field li
Load-bearing premise
The reconstruction stands or falls on the corona being close enough to force-free in the 5-G region that field lines traced from weak-field boundary pixels reflect real magnetic connections rather than numerical artifacts.
What would settle it
Take the same HMI boundary and seed regions and compute the field with a different NLFFF code or with chromospheric boundary data: the interpretation fails if the southern dimming region connects to the unconstrained twisted flux bundle in the preferred solution, or if the eastern MFR footprint does not fall inside the inverse J-shaped ribbon.
If this is right
- The pre-eruptive flux rope's eastern footprint, not the flare reconnection alone, sets where the core dimming appears and how large it can grow.
- Dimming growth into the strapping-field region is the observable signature of strapping-strapping reconnection, linking ribbon expansion to flux addition to the erupting rope.
- The western leg's multiple anchor points and overlying sunspot-connected field explain the suppressed dimming and partial confinement, so asymmetric eruptions can be diagnosed from pre-eruptive footpoint structure.
- PINN-based NLFFF extrapolation can be used on filaments that extend well beyond active-region cores, opening very large or quiet-Sun-rooted structures to quantitative coronal field modeling.
- The reconstructed free energy of about 3 x 10^32 erg and negative helicity are consistent with the inverse-S filament and its northern-hemisphere sign preference.
Where Pith is reading between the lines
- If footpoint geometry controls dimming, then for other whipping-like eruptions the pre-eruptive dimming region should be predictable from the fan of MFR field lines; a survey of similar events could test this without waiting for new instrumentation.
- The southern part of the dimming region is the fragile piece: in some sensitivity runs it connects to an unconstrained twisted flux bundle. A natural next test is to rerun the same event with chromospheric magnetogram constraints and see whether that connection disappears as the boundary information improves.
- This result suggests the NLFFF limit for filaments may lie beyond active-region cores; applying the same method to a fully quiescent filament, where one leg has no strong-field anchor at all, would test how far the force-free assumption can be pushed.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes the pre-eruptive magnetic configuration and asymmetric eruption of a roughly 500 Mm inverse S-shaped filament on 2023 February 24, partially rooted in AR 13229 and extending into weak-field regions. Using AIA, HMI, KSO, GONG, and STIX observations together with a physics-informed neural network (PINN) NLFFF extrapolation (NF2), the authors identify a large-scale magnetic flux rope whose high-current-density channel matches the observed filament. They report an extended eastern MFR footprint associated with the inverse J-shaped flare ribbon and the initial coronal dimming, while overlying strapping field lines connect to the region into which the dimming later expands. They interpret the dimming as stationary flux-rope and strapping-flux dimming, with later expansion driven by strapping-strapping and rope-strapping reconnection. The western leg shows multiple anchor points and stronger overlying fields, explaining the partial confinement and lack of dimming there. A 16-run sensitivity study over the force-free weighting factor lambda_ff is included in Appendix A.
Significance. If the interpretation holds, the paper provides a valuable demonstration that PINN-based NLFFF extrapolation can model large-scale filaments extending into weak-field regions, and it offers a physically coherent scenario connecting pre-eruptive MFR geometry to flare-ribbon and dimming asymmetries. The observational analysis is well documented, and the ensemble sensitivity study is a clear strength: the main MFR channel is robust for lambda_ff > 0.2, and the quantitative force-free metrics are reported for the full volume and for strong- and weak-field subregions. The central risk is that the load-bearing connectivity underlying the dimming-expansion scenario is the least robust part of the ensemble, and the model parameter lambda_ff was selected in part by matching the same observations used to validate the scenario. This does not undermine the existence of the MFR, but it weakens the specific claim about the strapping-field connectivity through the dimming region.
major comments (3)
- [Appendix A; Fig. A.2; Section 5] The interpretation in Section 5 that 'the area into which the dimming expanded is connected to strapping field lines that overlay the MFR's extended eastern leg' rests on exactly the part of the NLFFF ensemble that Appendix A shows to be unstable. The text states that for certain lambda_ff values the southernmost part of the dimming region is connected to an unconstrained twisted flux bundle, with no correlation to lambda_ff, and that lambda_ff = 0.60 gives a steep connectivity gradient that does not match observations. Since lambda_ff = 0.40 was selected partly because it produces the observed dimming boundary, the agreement of the selected model with the observed dimming is not an independent confirmation. I request that the authors either (i) reformulate the dimming-connectivity claim as one of several possible topologies with explicit uncertainty, or (ii) provide an independent test,
- [Sections 2.4, 4; Table A.1] The weak-field region that carries the dimming interpretation is also where the NLFFF assumptions are least secure. The mean unsigned flux in the final dimming mask is only ~5 G (Section 3.2), close to HMI noise levels, and Section 4 acknowledges that the real corona may deviate from a force-free state in such regions. Table A.1 shows that in the weak-field subregion E_div/E is about 4.9e-2 for the adopted lambda_ff = 0.40, i.e., near the 5% threshold used for reliable helicity computations, and the current-weighted angle theta_j does not improve systematically with lambda_ff there. The paper should quantify how the inferred connectivity changes under plausible perturbations of the weak-field boundary data (e.g., masking to the noise level) and should temper statements that the model 'effectively captures the essential large-scale connectivity' in this specific region.
- [Section 3.2; Fig. 6; Fig. A.2] The dimming mask is defined by a fixed LBR threshold of -0.5 (Eq. 1), and the field-line connectivity shown in Fig. 6 (right column) and Fig. A.2 (bottom row) uses that mask as seed regions. No sensitivity analysis is provided for this threshold. Because the connectivity conclusions concern a weak-field region where the dimensions and location of the mask are likely threshold-dependent, a modest change in the threshold could alter the seed region and hence the inferred strapping-field connectivity. The authors should show that their conclusions are robust to a reasonable range of LBR thresholds, or explicitly justify the chosen threshold physically.
minor comments (4)
- [Abstract and throughout] There are several typographical artifacts, e.g., 'e ffective' in the Abstract and 'whith' in Appendix A. These should be corrected in the final version.
- [Section 2.4, Eq. (7)] The total loss in Eq. (7) includes lambda_B0, which is decayed from 1000 to 1, but the text does not explain why this particular schedule is chosen or how sensitive the final solution is to the decay endpoint. A brief justification would improve reproducibility.
- [Section 2.3] The description of the CEA submap and its non-alignment with Carrington longitude/latitude is clear but would benefit from a small schematic or a more explicit statement of the projection-induced distortion in the top-right corner, since the authors note the 50-degree limitation there.
- [Section 4] The paper relies on the dimming classification of Veronig et al. (2025). It would be helpful to state explicitly which observational signatures, independent of the NLFFF model, support the stationary flux-rope versus moving flux-rope classification, so that the reader can separate the model-dependent and observation-based parts of the argument.
Circularity Check
Mild model-selection circularity for the southern dimming boundary, but the central MFR result is robust and independently corroborated.
specific steps
-
fitted input called prediction
[Appendix A (sensitivity analysis); selection used in Sect. 2.4; connectivity used in Sect. 3.3/Fig. 6 and Sect. 4]
"The value λff = 0.4 we selected for the results presented in this paper represents a balance between satisfying the quality metrics and aligning well with the observations."
The force-free weighting λff is a free parameter. The paper selects λff = 0.4 partly because it produces a connectivity boundary that aligns with the observed southern dimming boundary, and rejects λff = 0.60 because it 'does not match the observed evolution.' The same selected model is then used to infer the dimming-region connectivity (strapping fields) and to conclude that the area into which the dimming expanded is connected to strapping field lines. For this specific subregion, the model's agreement with the observation is therefore partly a selection criterion rather than an independent confirmation. However, the core MFR and the eastern footprint are robust across the ensemble, so the circularity is localized and not fatal to the main conclusion.
full rationale
The central claim—that a ~500 Mm MFR with an extended eastern footprint existed before the eruption—is an inference from a boundary-value NLFFF calculation and is robust across the 16-member λff ensemble. It is also corroborated by independent Hα/EUV observations of the filament shape and by the spatial correspondence of the MFR footprint with the inverse J-shaped ribbon. The self-citations to Jarolim et al. (2023) for the PINN method and to Veronig et al. (2025) for dimming categories are not load-bearing in a circular sense: the method is externally validated and applied to a new event, and the classification is a descriptive framework rather than a proof mechanism. The main circularity concern is the selection of λff = 0.4, which is chosen in part by visual agreement with the observed dimming boundary. Appendix A transparently shows that the southernmost part of the dimming region is connected to an unconstrained twisted flux bundle for some λff values, with no correlation to λff, indicating genuine difficulty finding a local force-free solution there. Because the paper uses the selected model to interpret exactly that dimming boundary, there is a mild validation circularity. But the broader MFR and strapping-flux interpretation remain supported by ensemble-robust features and by the qualitative agreement with the filament and flare-ribbon observations. Thus the paper is largely self-contained and not significantly circular, aside from this localized model-selection issue.
Axiom & Free-Parameter Ledger
free parameters (2)
- lambda_ff (force-free weighting factor) =
0.4
- Dimming threshold (LBR < -0.5) =
-0.5 log10 relative counts
axioms (5)
- domain assumption The coronal magnetic field is force-free (J x B = 0) over most of the volume, including weak-field regions.
- domain assumption The side and top boundaries of the computational volume can be approximated by a potential field solution.
- domain assumption The pre-eruption magnetic configuration at 19:34 UT is representative of the structure that erupted at about 20:03 UT.
- domain assumption HMI photospheric vector magnetogram, after error-masked clipping (Bdiff,clipped), is a reliable lower boundary in weak-field regions.
- domain assumption Cartesian geometry is adequate for this field of view because the filament lies within 50 degrees of disk center.
Cite this review
Pith. "Pith review of Magnetic structure and asymmetric eruption of a 500 Mm filament rooted in weak-field regions." pith.science (2026). https://pith.science/paper/S7IJI4QL
@misc{pith2026250818121,
author = {Pith},
title = {Pith review of: Magnetic structure and asymmetric eruption of a 500 Mm filament rooted in weak-field regions},
year = {2026},
howpublished = {\url{https://pith.science/paper/S7IJI4QL}},
note = {Machine review of arXiv:2508.18121}
}
read the original abstract
We performed a detailed analysis of the magnetic structure and asymmetric eruption of a large (about 500 Mm) inverse S-shaped filament partially located in AR 13229 on February 24, 2023. We linked the filament's pre-eruptive magnetic configuration to its highly asymmetric eruption dynamics and the formation of a large-scale coronal dimming in a weak-field region (mean unsigned flux of about 5 G). To reconstruct the coronal magnetic field, we applied a physics-informed neural network (PINN)-based nonlinear force-free field (NLFFF) extrapolation method to a pre-eruption HMI vector magnetogram. The NLFFF extrapolation reveals a large-scale magnetic flux rope (MFR) of about 500 Mm in length, consistent with the filament. We identified an extended MFR footprint to the east that connects to the J-shaped flare ribbon, outlining where the coronal dimming began. Overlying strapping fields connect to the area into which the dimming and flare ribbon later expand. This configuration explains the formation of the dimming as a stationary flux rope and strapping flux dimming, with subsequent expansion driven by the growth of the MFR footprint through strapping-strapping reconnection. Conversely, the western filament leg shows multiple anchor points and strong overlying magnetic fields, which suppressed the dimming and partially confined the eruption on that side. The reconstructed pre-eruptive NLFFF configuration offers a clear physical explanation for the asymmetries seen in the eruption, flare geometry, and coronal dimming. This demonstrates that PINN-based NLFFF extrapolation can effectively model large-scale filaments extending into weak-field regions, enhancing our understanding of complex solar eruptions.
Reference graph
Works this paper leans on
-
[1]
Afanasyev, A. N., Fan, Y ., Kazachenko, M. D., & Cheung, M. C. M. 2023, ApJ, 952, 136
work page 2023
-
[2]
Ahamed, A. A., Subramanian, S. P., Rahman, A. M., et al. 2025, New A, 114, 102312
work page 2025
-
[3]
Bobra, M. G., Sun, X., Hoeksema, J. T., et al. 2014, Sol. Phys., 289, 3549 Démoulin, P., Priest, E. R., & Lonie, D. P. 1996, J. Geophys. Res., 101, 7631
work page 2014
-
[4]
L., Wheatland, M
DeRosa, M. L., Wheatland, M. S., Leka, K. D., et al. 2015, ApJ, 811, 107
2015
-
[5]
Gary, G. A. 2001, Sol. Phys., 203, 71
work page 2001
-
[6]
Gary, G. A. & Hagyard, M. J. 1990, Sol. Phys., 126, 21
1990
-
[7]
Gibson, S. E. 2018, Living Reviews in Solar Physics, 15, 7
work page 2018
-
[8]
M., Török, T., Vršnak, B., Manchester, W., & Veronig, A
Green, L. M., Török, T., Vršnak, B., Manchester, W., & Veronig, A. 2018, Space Sci. Rev., 214, 46
work page 2018
-
[9]
W., Hill, F., Hubbard, R
Harvey, J. W., Hill, F., Hubbard, R. P., et al. 1996, Science, 272, 1284
1996
- [10]
-
[11]
T., Liu, Y ., Hayashi, K., et al
Hoeksema, J. T., Liu, Y ., Hayashi, K., et al. 2014, Sol. Phys., 289, 3483
2014
-
[12]
S., Acton, L
Hudson, H. S., Acton, L. W., & Freeland, S. L. 1996, ApJ, 470, 629
1996
- [13]
-
[14]
Janvier, M., Aulanier, G., Bommier, V ., et al. 2014, ApJ, 788, 60
work page 2014
-
[15]
Jarolim, R., Thalmann, J. K., Veronig, A. M., & Podladchikova, T. 2023, Nature Astronomy, 7, 1171
work page 2023
-
[16]
2022, The Innovation, 3, 100236
Jiang, C., Feng, X., Guo, Y ., & Hu, Q. 2022, The Innovation, 3, 100236
work page 2022
- [17]
-
[18]
J., Grimm, O., et al
Krucker, S., Hurford, G. J., Grimm, O., et al. 2020, A&A, 642, A15
2020
-
[19]
R., Title, A
Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17
2012
-
[20]
Liu, R., Alexander, D., & Gilbert, H. R. 2009, ApJ, 691, 1079
work page 2009
-
[21]
Martens, P. C. & Zwaan, C. 2001, ApJ, 558, 872
work page 2001
-
[22]
Mastrano, A., Yang, K. E., & Wheatland, M. S. 2020, Sol. Phys., 295, 97 Müller, D., St. Cyr, O. C., Zouganelis, I., et al. 2020, A&A, 642, A1
work page 2020
-
[23]
J., Freij, N., Stansby, D., et al
Mumford, S. J., Freij, N., Stansby, D., et al. 2023, SunPy
work page 2023
-
[24]
2014, Living Reviews in Solar Physics, 11, 1
Parenti, S. 2014, Living Reviews in Solar Physics, 11, 1
2014
-
[25]
D., Thompson, B
Pesnell, W. D., Thompson, B. J., & Chamberlin, P. C. 2012, Sol. Phys., 275, 3
2012
-
[26]
Peter, H., Warnecke, J., Chitta, L. P., & Cameron, R. H. 2015, A&A, 584, A68 Pötzi, W., Veronig, A., Jarolim, R., et al. 2021, Sol. Phys., 296, 164
work page 2015
-
[27]
A., et al
Purkhart, S., Collier, H., Hayes, L. A., et al. 2025, A&A, 698, A318
2025
-
[28]
Purkhart, S., Veronig, A. M., Dickson, E. C. M., et al. 2023, A&A, 679, A99
work page 2023
-
[29]
Purkhart, S., Veronig, A. M., Kliem, B., et al. 2024, A&A, 689, A259
work page 2024
-
[30]
Raissi, M., Perdikaris, P., & Karniadakis, G. E. 2019, Journal of Computational Physics, 378, 686 Rodríguez-Gómez, J. M., Kuckein, C., González Manrique, S. J., et al. 2024, ApJ, 964, 27
work page 2019
-
[31]
H., Bush, R
Schou, J., Scherrer, P. H., Bush, R. I., et al. 2012, Sol. Phys., 275, 229
2012
-
[32]
Schrijver, C. J., De Rosa, M. L., Metcalf, T. R., et al. 2006, Sol. Phys., 235, 161
work page 2006
- [33]
-
[34]
Sterling, A. C. & Hudson, H. S. 1997, ApJ, 491, L55
1997
-
[35]
K., Dumbovi´c, M., Dissauer, K., et al
Thalmann, J. K., Dumbovi´c, M., Dissauer, K., et al. 2023, A&A, 669, A72
work page 2023
-
[36]
K., Linan, L., Pariat, E., & Valori, G
Thalmann, J. K., Linan, L., Pariat, E., & Valori, G. 2019, ApJ, 880, L6 The SunPy Community, Barnes, W. T., Bobra, M. G., et al. 2020, The Astrophys- ical Journal, 890, 68 Trujillo Bueno, J. & del Pino Alemán, T. 2022, ARA&A, 60, 415
work page 2019
-
[37]
Valori, G., Démoulin, P., Pariat, E., & Masson, S. 2013, A&A, 553, A38
work page 2013
-
[38]
M., Dissauer, K., Kliem, B., et al
Veronig, A. M., Dissauer, K., Kliem, B., et al. 2025, Living Reviews in Solar Physics, 22, 2
work page 2025
-
[39]
Wheatland, M. S. & Régnier, S. 2009, ApJ, 700, L88
work page 2009
- [40]
-
[41]
Wiegelmann, T., Petrie, G. J. D., & Riley, P. 2017, Space Sci. Rev., 210, 249
2017
-
[42]
Wiegelmann, T., Thalmann, J. K., Inhester, B., et al. 2012, Sol. Phys., 281, 37 Article number, page 13 of 15 A&A proofs: manuscript no. aa55860-25corr Appendix A: Sensitivity analysis for PINN-based NLFFF modeling To assess the robustness of the NLFFF extrapolation results shown in Sect. 3.3, especially their sensitivity to the force-free weighting facto...
work page 2012
This paper was first reviewed by deepseek-v4-flash on August 5, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.