pith. sign in

arxiv: 2604.16195 · v1 · submitted 2026-04-17 · 🌌 astro-ph.SR

FastQSL 2: A Comprehensive Toolkit for Magnetic Connectivity Analysis

Jun Chen , Thomas Wiegelmann , Li Feng , Chaowei Jiang , Rui Liu This is my paper

Pith reviewed 2026-05-10 07:14 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords quasi-separatrix layersQSLmagnetic field line tracingspherical coordinatessolar magnetic connectivitymagnetic reconnectionsolar wind modelingfield line footpoints
0
0 comments X

The pith

FastQSL 2 traces magnetic field lines over the full sphere by switching to a second coordinate system near the poles.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces an updated FastQSL toolkit for locating quasi-separatrix layers in magnetic fields, where strong connectivity gradients favor current buildup and reconnection. It extends the original code to spherical coordinates by employing a second spherical coordinate system for field line tracing around the polar regions. This change eliminates the singularity that blocked accurate calculations at the north and south poles. The toolkit further supports arbitrary output meshes, supplies magnetic field and current density values, and exports traced field lines. These additions expand the tool's reach to global solar magnetic connectivity studies and allow derivation of footpoint-based parameters for solar wind speed modeling.

Core claim

FastQSL 2 supports spherical coordinates by utilizing a second spherical coordinate system for tracing magnetic field lines around the polar regions. This approach completely resolves the singularity problem at the two poles while accommodating arbitrary mesh shapes for output, providing both magnetic field and electric current density on the mesh, and saving the traced magnetic field lines.

What carries the argument

A second spherical coordinate system that switches in for field line tracing near the poles to avoid the coordinate singularity.

Load-bearing premise

Numerical field line integration stays accurate and stable when the code switches between the two spherical coordinate systems across the full sphere.

What would settle it

A direct comparison of field lines traced by FastQSL 2 against known analytic solutions for a global dipole field on a sphere, checking for mismatches or instabilities near the poles.

Figures

Figures reproduced from arXiv: 2604.16195 by Chaowei Jiang, Jun Chen, Li Feng, Rui Liu, Thomas Wiegelmann.

Figure 1
Figure 1. Figure 1: Left: two-spherical-coordinate system. Right: the mesh around a polar region. The grids of φ, ϑ are shown by the crosses of the gray mesh; the grids of φ2, ϑ2 are shown by the crosses of the purple mesh. In this spherical coordinate system, ∇ × B = eφ r  ∂ Br ∂ϑ − ∂ (r Bϑ) ∂r  + eϑ r  ∂ (r Bφ) ∂r − 1 cos ϑ ∂ Br ∂φ  + er r cos ϑ  ∂ Bϑ ∂φ − ∂ (cos ϑ Bφ) ∂ϑ  . (12) A magnetic field line can be traced by… view at source ↗
Figure 2
Figure 2. Figure 2: Left: Br at r = 1.04 Rsun with selected magnetic field lines from both the inserted flux rope and the background field overlaid; the distribution of the twist number is displayed on a planar cross-section. Right: the distributions of Q on the surface of r = 1.04 Rsun, the same planar cross-section, and the surface of a hollowed-out dome; the purple region marks one segment of a separator detected with Qloc… view at source ↗
Figure 3
Figure 3. Figure 3: The Q-map at r = 1.04 Rsun of the same field used in [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The distributions of Qlocal for rlocal = 0.05, 0.2, 1, 2 in the plane of y = 0 with Bquadrupole. Arge et al. (2003) found the solar wind speed at the first Lagrangian point from December 1994 to the end of 1995 can be roughly modeled by vsw = 265 + 25 f 2/7 s [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Parameters used for modeling solar wind speed. The magnetic field is extrapolated from the synoptic magnetogram of FMG using PFSS Model; the field of Quadrupole2 is not embedded here. R0 = 1.04 Rsun, R1 = 2.5 Rsun. (a) Br at r = R0. (b) Target point type. (c) fs at r = R0. For a closed field line, its fs is set to 1000. (d) θb at r = R0. For a closed field line, its θb is set to 0. 3.3.2. Slip-Squashing Fa… view at source ↗
Figure 6
Figure 6. Figure 6: successfully reproduces the plots of [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
read the original abstract

We present a new version of FastQSL for locating quasi-separatrix layers (QSLs) -- regions characterized by strong magnetic connectivity gradients, preferential current buildup, and subsequent magnetic reconnection. This version now supports spherical coordinates, utilizing a second spherical coordinate system for tracing magnetic field lines around the polar regions. This approach completely resolves the singularity problem at the two poles. Furthermore, our code accommodates arbitrary mesh shapes for output, can provide both magnetic field and electric current density on the mesh, and can save the traced magnetic field lines. We suggest using $Q_\mathrm{local}$ calculated through a localized mapping to locate (quasi-)separators. By quickly and accurately outputting the footpoint coordinates of magnetic field lines, FastQSL can be used to derive the two key parameters used for modeling solar wind speed and slip-squashing factors for the case of zero boundary flow. Compared with the first version, FastQSL 2 achieves significant improvements in terms of application scope.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 0 minor

Summary. The manuscript presents FastQSL 2, an updated toolkit for locating quasi-separatrix layers (QSLs) in magnetic fields. It claims support for spherical coordinates via a second spherical coordinate system that completely resolves polar singularities during field-line tracing, accommodation of arbitrary output meshes, output of magnetic field and current density, saving of traced field lines, and use of Q_local for (quasi-)separators. The code enables derivation of solar wind speed parameters and slip-squashing factors, representing a significant expansion in scope over the first version.

Significance. If the numerical claims hold, the toolkit would be a useful addition for solar physics, enabling full-sphere QSL analysis without polar artifacts and supporting large-scale connectivity studies relevant to reconnection and solar wind modeling.

major comments (1)
  1. [Abstract] Abstract: The assertion that the second spherical coordinate system 'completely resolves' the polar singularity lacks any supporting error metrics, convergence tests, stability analysis during coordinate switches, or validation against analytic solutions (e.g., dipole field). This is load-bearing for the central improvement claimed over FastQSL 1.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive feedback on our manuscript. We address the single major comment point by point below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The assertion that the second spherical coordinate system 'completely resolves' the polar singularity lacks any supporting error metrics, convergence tests, stability analysis during coordinate switches, or validation against analytic solutions (e.g., dipole field). This is load-bearing for the central improvement claimed over FastQSL 1.

    Authors: We agree that the abstract would benefit from explicit references to supporting validation. The method relies on a mathematically constructed secondary spherical coordinate system whose poles are offset from those of the primary system, thereby removing the coordinate singularity by design during field-line tracing. In the revised manuscript we will update the abstract to reference the numerical validation against an analytic dipole field (including quantitative error metrics near the poles) and a brief stability check on the coordinate transformation, both of which appear in Section 3 of the current text. These additions will be incorporated without lengthening the abstract substantially. revision: yes

Circularity Check

0 steps flagged

No circularity: computational toolkit paper with no derivations or fitted predictions

full rationale

The manuscript describes an updated software implementation (FastQSL 2) for QSL mapping, including a dual spherical coordinate scheme for polar regions. No mathematical derivation chain, first-principles predictions, parameter fitting, or output quantities that reduce to inputs by construction are present. References to the prior version serve only for comparison of scope; they do not supply load-bearing premises or uniqueness theorems. The work is self-contained as a methods description and contains none of the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper relies on standard numerical methods for field-line integration and coordinate transformations without introducing new free parameters, physical axioms beyond conventional MHD, or invented entities.

axioms (1)
  • standard math Magnetic field lines are traced by numerical integration of the vector field
    Standard technique invoked for connectivity mapping; no specific section quoted as abstract-only.

pith-pipeline@v0.9.0 · 5477 in / 1201 out tokens · 54335 ms · 2026-05-10T07:14:44.834376+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

34 extracted references · 34 canonical work pages

  1. [1]

    N., Odstrcil, D., Pizzo, V

    Arge, C. N., Odstrcil, D., Pizzo, V. J., & Mayer, L. R. 2003, in American Institute of Physics Conference Series, Vol. 679, Solar Wind Ten, ed. M. Velli, R. Bruno, F. Malara, & B. Bucci (AIP), 190–193, doi: 10.1063/1.1618574

    work page doi:10.1063/1.1618574 2003

  2. [2]

    B., Wilkins, C

    Aslanyan, V., Scott, R. B., Wilkins, C. P., et al. 2024, The Astrophysical Journal, 971, 137, doi: 10.3847/1538-4357/ad55ca

    work page doi:10.3847/1538-4357/ad55ca 2024

  3. [3]

    2005, A&A, 444, 961, doi: 10.1051/0004-6361:20053600

    Aulanier, G., Pariat, E., & D´ emoulin, P. 2005, A&A, 444, 961, doi: 10.1051/0004-6361:20053600

    work page doi:10.1051/0004-6361:20053600 2005

  4. [4]

    A., & Prior, C

    Berger, M. A., & Prior, C. 2006, Journal of Physics A: Mathematical and General, 39, 8321

    work page 2006

  5. [5]

    M., Downs, C., Linker, J

    Caplan, R. M., Downs, C., Linker, J. A., & Mikic, Z. 2021, ApJ, 915, 44, doi: 10.3847/1538-4357/abfd2f 18Chen et al

    work page doi:10.3847/1538-4357/abfd2f 2021

  6. [6]

    2023, The Astrophysical Journal Letters, 951, L35

    Chen, J., Cheng, X., Kliem, B., & Ding, M. 2023, The Astrophysical Journal Letters, 951, L35

    work page 2023

  7. [7]

    2020, The Astrophysical Journal, 890, 158 D´ emoulin, P., Bagala, L

    Chen, J., Liu, R., Liu, K., et al. 2020, The Astrophysical Journal, 890, 158 D´ emoulin, P., Bagala, L. G., Mandrini, C. H.,

    work page 2020

  8. [8]

    C., & Rovira, M

    Henoux, J. C., & Rovira, M. G. 1997, A&A, 325, 305 D´ emoulin, P., Priest, E., & Lonie, D. 1996, Journal of Geophysical Research: Space Physics, 101, 7631

    work page 1997

  9. [9]

    2025, Chinese Journal of Space Science, 45, 913, doi: 10.11728/cjss2025.04.2025-0054

    Deng, Y., Tian, H., Jiang, J., et al. 2025, Chinese Journal of Space Science, 45, 913, doi: 10.11728/cjss2025.04.2025-0054

    work page doi:10.11728/cjss2025.04.2025-0054 2025

  10. [10]

    2019, Research in Astronomy and Astrophysics, 19, 157, doi: 10.1088/1674-4527/19/11/157

    Deng, Y.-Y., Zhang, H.-Y., Yang, J.-F., et al. 2019, Research in Astronomy and Astrophysics, 19, 157, doi: 10.1088/1674-4527/19/11/157

    work page doi:10.1088/1674-4527/19/11/157 2019

  11. [11]

    1969, Low-order classical Runge-Kutta formulas with stepsize control and their application to some heat transfer problems, Vol

    Fehlberg, E. 1969, Low-order classical Runge-Kutta formulas with stepsize control and their application to some heat transfer problems, Vol. 315 (National aeronautics and space administration)

    work page 1969

  12. [12]

    2023, Solar Physics, 298, 68

    Gan, W., Zhu, C., Deng, Y., et al. 2023, Solar Physics, 298, 68

    work page 2023

  13. [13]

    2020, in APS Meeting Abstracts, Vol

    Gekelman, W., Dehaas, T., Prior, C., & Yeates, A. 2020, in APS Meeting Abstracts, Vol. 2020, APS Division of Plasma Physics Meeting

    work page 2020

  14. [14]

    Hoeksema, J. T. 1984, PhD thesis, Stanford

    work page 1984

  15. [15]

    Kutta, W. 1901, Z. Math. Phys., 46

    work page 1901

  16. [16]

    S., et al

    Liu, R., Kliem, B., Titov, V. S., et al. 2016, The Astrophysical Journal, 818, 148

    work page 2016

  17. [17]

    W., & Strauss, H

    Longcope, D. W., & Strauss, H. R. 1994a, ApJ, 426, 742, doi: 10.1086/174111 —. 1994b, ApJ, 437, 851, doi: 10.1086/175045

    work page doi:10.1086/174111

  18. [18]

    1976, American Journal of Physics, 44, 1210, doi: 10.1119/1.10264

    Mathews, J. 1976, American Journal of Physics, 44, 1210, doi: 10.1119/1.10264

    work page doi:10.1119/1.10264 1976

  19. [19]

    2019, A&A, 624, A51, doi: 10.1051/0004-6361/201834668

    Moraitis, K., Pariat, E., Valori, G., & Dalmasse, K. 2019, A&A, 624, A51, doi: 10.1051/0004-6361/201834668

    work page doi:10.1051/0004-6361/201834668 2019

  20. [20]

    2012, A&A, 541, A78, doi: 10.1051/0004-6361/201118515

    Pariat, E., & D´ emoulin, P. 2012, A&A, 541, A78, doi: 10.1051/0004-6361/201118515

    work page doi:10.1051/0004-6361/201118515 2012

  21. [21]

    Pontin, D. I. 2011, Advances in Space Research, 47, 1508, doi: 10.1016/j.asr.2010.12.022

    work page doi:10.1016/j.asr.2010.12.022 2011

  22. [22]

    2000, Magnetic reconnection : MHD theory and applications

    Priest, E., ed. 2000, Magnetic reconnection : MHD theory and applications

    work page 2000

  23. [23]

    1995, Journal of Geophysical Research: Space Physics, 100, 23443

    Priest, E., & D´ emoulin, P. 1995, Journal of Geophysical Research: Space Physics, 100, 23443

    work page 1995

  24. [24]

    B., Pontin, D

    Scott, R. B., Pontin, D. I., & Hornig, G. 2017, ApJ, 848, 117, doi: 10.3847/1538-4357/aa8a64

    work page doi:10.3847/1538-4357/aa8a64 2017

  25. [25]

    2017, The Astrophysical Journal, 840, 89

    Tassev, S., & Savcheva, A. 2017, The Astrophysical Journal, 840, 89

    work page 2017

  26. [26]

    Titov, V. S. 2007, ApJ, 660, 863, doi: 10.1086/512671

    work page doi:10.1086/512671 2007

  27. [27]

    S., Forbes, T

    Titov, V. S., Forbes, T. G., Priest, E. R., Miki´ c, Z., & Linker, J. A. 2009, ApJ, 693, 1029, doi: 10.1088/0004-637X/693/1/1029

    work page doi:10.1088/0004-637x/693/1/1029 2009

  28. [28]

    S., Hornig, G., & Démoulin, P

    Titov, V. S., Hornig, G., & D´ emoulin, P. 2002, Journal of Geophysical Research (Space Physics), 107, 1164, doi: 10.1029/2001JA000278

    work page doi:10.1029/2001ja000278 2002

  29. [29]

    M., & Sheeley, Jr., N

    Wang, Y. M., & Sheeley, Jr., N. R. 1992, ApJ, 392, 310, doi: 10.1086/171430

    work page doi:10.1086/171430 1992

  30. [30]

    N., Hock, R., Eparvier, F., et al

    Woods, T. N., Hock, R., Eparvier, F., et al. 2011, ApJ, 739, 59, doi: 10.1088/0004-637X/739/2/59

    work page doi:10.1088/0004-637x/739/2/59 2011

  31. [31]

    E., Xia, C., & Guo, Y

    Yang, K. E., Xia, C., & Guo, Y. 2024, Kai-E-Yang/QSL: QSL v1, Zenodo, doi: 10.5281/ZENODO.13831036

    work page doi:10.5281/zenodo.13831036 2024

  32. [32]

    R., & Page, M

    Yeates, A. R., & Page, M. H. 2018, Journal of Plasma Physics, 84, 775840602, doi: 10.1017/S0022377818001204

    work page doi:10.1017/s0022377818001204 2018

  33. [33]

    2022, ApJ, 937, 26, doi: 10.3847/1538-4357/ac8d61

    Zhang, P., Chen, J., Liu, R., & Wang, C. 2022, ApJ, 937, 26, doi: 10.3847/1538-4357/ac8d61

    work page doi:10.3847/1538-4357/ac8d61 2022

  34. [34]

    2019, Annales Geophysicae, 37, 325, doi: 10.5194/angeo-37-325-2019

    Zhu, P., Wang, Z., Chen, J., Yan, X., & Liu, R. 2019, Annales Geophysicae, 37, 325, doi: 10.5194/angeo-37-325-2019

    work page doi:10.5194/angeo-37-325-2019 2019