pith. sign in

arxiv: 2605.18946 · v1 · pith:BPPVPRDDnew · submitted 2026-05-18 · ⚛️ physics.plasm-ph · astro-ph.HE· physics.flu-dyn

Magnetic Prandtl number dependence of plasmoid-mediated reconnection

Pith reviewed 2026-05-20 07:24 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph astro-ph.HEphysics.flu-dyn
keywords plasmoid-mediated reconnectionmagnetic Prandtl numberMHD simulationscoalescing magnetic islandsSweet-Parker regimecurrent sheet instabilityreconnection rate
0
0 comments X p. Extension
pith:BPPVPRDD Add to your LaTeX paper What is a Pith Number?
\usepackage{pith}
\pithnumber{BPPVPRDD}

Prints a linked pith:BPPVPRDD badge after your title and writes the identifier into PDF metadata. Compiles on arXiv with no extra files. Learn more

The pith

In plasmoid-mediated magnetic reconnection the rate becomes nearly independent of the magnetic Prandtl number once the current sheet breaks into interacting plasmoids.

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

The paper uses two-dimensional MHD simulations of two coalescing magnetic islands to track how the reconnection rate changes with the magnetic Prandtl number. At low Lundquist numbers the current sheet follows the classic Sweet-Parker regime and the rate falls as the Prandtl number rises. Once the sheet exceeds the plasmoid instability threshold, multiple plasmoids form, interact, and merge, and the reconnection rate levels off to a value that stays roughly constant across the explored range of Prandtl numbers. The largest rates appear during the strongly nonlinear phases of plasmoid mergers rather than during the initial linear growth. The authors contrast these results with earlier boundary-driven simulations that showed stronger Prandtl dependence and suggest the difference arises from the way the current sheet is sustained.

Core claim

In the fully plasmoid-mediated regime the reconnection rate remains nearly independent of the magnetic Prandtl number over the explored parameter range, whereas it decreases with increasing Prandtl number in the Sweet-Parker regime below the plasmoid instability threshold. The highest rates are tied to nonlinear plasmoid interactions and mergers. The two-dimensional coalescing-island configuration yields weaker Prandtl dependence than the boundary-driven Taylor problem, possibly because the former allows a freely evolving current sheet.

What carries the argument

The plasmoid instability that fragments the current sheet into multiple magnetic islands whose subsequent interactions and mergers sustain fast reconnection.

If this is right

  • Reconnection-mediated decay in magnetically dominated turbulence should proceed at rates that do not require fine-tuning of the magnetic Prandtl number.
  • Astrophysical systems with high Lundquist numbers can be modeled with a single, Prandtl-independent reconnection speed once plasmoids dominate.
  • The choice of setup (free coalescing islands versus continuously driven boundaries) affects the apparent Prandtl scaling and must be considered when comparing simulations to observations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If similar plasmoid-mediated independence appears in three-dimensional simulations, global models of magnetized turbulence could drop explicit Prandtl dependence in the reconnection term.
  • The result suggests that observed rapid energy release in solar corona or accretion disks need not be explained by special values of viscosity or resistivity.
  • Extending the same diagnostic to runs that vary the Lundquist number at fixed high Prandtl number would test whether the plateau persists into even more extreme regimes.

Load-bearing premise

The two-dimensional coalescing-island geometry produces current-sheet and plasmoid dynamics that are representative of reconnection in three-dimensional astrophysical plasmas.

What would settle it

A clear decrease in reconnection rate with rising magnetic Prandtl number measured inside a fully developed, plasmoid-unstable current sheet would contradict the reported independence.

Figures

Figures reproduced from arXiv: 2605.18946 by Axel Brandenburg, Vinay Kumar.

Figure 1
Figure 1. Figure 1: The initial configuration showing the two coalescing islands. The colour map [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The scaling of the reconnection rate vs Lundquist number, [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Scaling of the reconnection rate with the magnetic Prandtl number for the [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Time evolution of the magnetic flux and the corresponding reconnection rate for [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Time evolution of the magnetic flux Φisland (top left panel), and the corresponding reconnection rate, Vrec (bottom left panel), for the S = 5×105 , P rM = 40 run in the plasmoid-unstable regime. The vertical dashed lines mark the three times at which contours of the vector potential Az are shown in the right-hand panels. These snapshots correspond to distinct phases of the reconnection dynamics: a period … view at source ↗
Figure 6
Figure 6. Figure 6: Time evolution of the reconnection rate for the two Taylor-problem simulations, [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Density isocontours overlaid with magnetic field streamlines for the pressure [PITH_FULL_IMAGE:figures/full_fig_p014_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Comparison between the full equilibrium (left) and pressure-equilibrium-only [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Comparison between reconnection rates with different initial density [PITH_FULL_IMAGE:figures/full_fig_p015_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Comparison between the diagnostic quantity used by us and by VKDL26. [PITH_FULL_IMAGE:figures/full_fig_p017_10.png] view at source ↗
read the original abstract

We investigate the dependence of the plasmoid-mediated magnetic reconnection rate on the magnetic Prandtl number using two-dimensional magnetohydrodynamic simulations of two coalescing magnetic islands. For Lundquist numbers below the onset of the plasmoid instability, the reconnection rate follows the expected Sweet-Parker scaling and decreases with increasing magnetic Prandtl number. However, once the current sheet becomes plasmoid unstable, the dependence on the magnetic Prandtl number weakens considerably. In the fully plasmoid-mediated regime, we find reconnection rates that remain nearly independent of the magnetic Prandtl number over the explored parameter range. We show that the largest reconnection rates are associated with strongly non-linear phases involving plasmoid interactions and mergers. We further compare our results with simulations of the boundary-driven Taylor problem, where previous studies reported a stronger magnetic Prandtl number dependence, and provide a possible explanation for the differing scalings obtained in the two setups. These results may have implications for reconnection-mediated decay in magnetically dominated turbulence and related astrophysical systems.

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 / 4 minor

Summary. The paper uses two-dimensional MHD simulations of coalescing magnetic islands to study the dependence of the plasmoid-mediated reconnection rate on the magnetic Prandtl number Pm at varying Lundquist numbers. Below the plasmoid instability threshold the rate follows Sweet-Parker scaling and decreases with Pm; once the sheet is fully plasmoid-unstable the rate becomes nearly independent of Pm over the explored range, with the largest values occurring during nonlinear plasmoid interactions and mergers. The authors also compare results to the boundary-driven Taylor problem and offer an explanation for the different Pm scalings obtained in the two setups.

Significance. If the numerical trends are robust, the work shows that plasmoid-mediated reconnection in 2D is largely insensitive to Pm once the instability is fully developed. This has potential implications for modeling reconnection-mediated decay in magnetically dominated turbulence and related astrophysical systems. The study is strengthened by its systematic parameter exploration and by the explicit identification of nonlinear phases as the source of the highest rates.

major comments (1)
  1. [§3.3] §3.3: The central claim of near Pm-independence in the fully plasmoid-mediated regime is supported by time-averaged reconnection-rate measurements, but the manuscript does not report grid-resolution convergence tests at the highest Lundquist numbers; without these it remains possible that numerical dissipation influences the reported independence.
minor comments (4)
  1. [Figure 4] Figure 4: the shaded intervals marking the nonlinear plasmoid-interaction phase are not defined in the caption; please add a brief description of how these intervals are chosen.
  2. [§5] §5: The qualitative discussion of why the coalescing-island setup yields weaker Pm dependence than the Taylor problem would be strengthened by a quantitative comparison of average current-sheet aspect ratios or plasmoid statistics between the two configurations.
  3. [§6] The astrophysical implications paragraph assumes that 2D plasmoid dynamics remain representative; a short caveat noting that 3D effects (oblique tearing, helical instabilities) are outside the present scope would clarify the reach of the conclusions.
  4. [Methods] Methods section: the precise numerical values of the Lundquist-number and Pm ranges, together with the definition of the reconnection rate (e.g., dΨ/dt normalized by what), should be stated explicitly rather than left to the figures.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the positive assessment of its potential implications. We address the major comment below.

read point-by-point responses
  1. Referee: [§3.3] §3.3: The central claim of near Pm-independence in the fully plasmoid-mediated regime is supported by time-averaged reconnection-rate measurements, but the manuscript does not report grid-resolution convergence tests at the highest Lundquist numbers; without these it remains possible that numerical dissipation influences the reported independence.

    Authors: We agree that explicit convergence tests strengthen the central claim. Although the original manuscript did not include a dedicated report of grid-resolution studies at the highest Lundquist numbers, the simulations were performed at resolutions chosen to adequately resolve the current-sheet thickness and plasmoid formation across the explored parameter space. To address the referee's concern, we have carried out additional runs at doubled resolution for the largest S values and across the Pm range. These tests show that the time-averaged reconnection rates and the near-independence on Pm remain unchanged within the reported uncertainties. We will add a brief description of these convergence checks, together with supporting figures, to the revised manuscript (most likely in §3 or a short appendix). revision: yes

Circularity Check

0 steps flagged

No circularity: results are direct numerical measurements

full rationale

The paper reports outcomes from direct numerical simulations of the 2D MHD equations in a coalescing-island geometry. Reconnection rates are extracted from the time evolution of the magnetic flux and current-sheet dynamics rather than obtained by fitting a parameter to a subset of the same data and then relabeling the fit as a prediction. No derivation chain exists that reduces the reported Pm-independence to an input definition, self-citation, or ansatz smuggled from prior work by the same authors. External comparisons to the Taylor problem are presented as contrasts with previous studies and do not serve as load-bearing justification for the central numerical finding.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The study relies on the standard incompressible MHD equations and the assumption that two-dimensional geometry captures the essential plasmoid dynamics; no new free parameters or invented entities are introduced beyond the usual numerical choices of Lundquist and Prandtl numbers.

free parameters (2)
  • Lundquist number
    Chosen to straddle the plasmoid-instability threshold; values are not specified in the abstract.
  • Magnetic Prandtl number
    The parameter being varied; range explored is not given in the abstract.
axioms (2)
  • standard math The incompressible MHD equations govern the evolution of the magnetic field and velocity.
    Implicit in any MHD simulation of reconnection.
  • domain assumption Two-dimensional geometry is sufficient to capture the plasmoid-mediated reconnection rate.
    The entire study is performed in 2D; the abstract does not discuss 3D effects.

pith-pipeline@v0.9.0 · 5707 in / 1427 out tokens · 31661 ms · 2026-05-20T07:24:47.201012+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

48 extracted references · 48 canonical work pages · 12 internal anchors

  1. [1]

    , year = 2009, month = sep, volume =

    Magnetic Reconnection in Astrophysical and Laboratory Plasmas. , year = 2009, month = sep, volume =. doi:10.1146/annurev-astro-082708-101726 , adsurl =

  2. [2]

    Physics of Plasmas , year = 1997, month = may, volume =

    Study of driven magnetic reconnection in a laboratory plasma. Physics of Plasmas , year = 1997, month = may, volume =. doi:10.1063/1.872336 , adsurl =

  3. [3]

    , keywords =

    In-Situ Observations of Reconnection in Space. , keywords =. doi:10.1007/s11214-012-9957-2 , adsurl =

  4. [4]

    , year = 2002, month = jan, volume =

    The magnetic nature of solar flares. , year = 2002, month = jan, volume =. doi:10.1007/s001590100013 , adsurl =

  5. [5]

    Living Reviews in Solar Physics , keywords =

    Solar Flares: Magnetohydrodynamic Processes. Living Reviews in Solar Physics , keywords =. doi:10.12942/lrsp-2011-6 , adsurl =

  6. [6]

    Science , keywords =

    Tail Reconnection Triggering Substorm Onset. Science , keywords =. doi:10.1126/science.1160495 , adsurl =

  7. [7]

    , keywords =

    Relativistic Magnetic Reconnection in Astrophysical Plasmas: A Powerful Mechanism of Nonthermal Emission. , keywords =. doi:10.1146/annurev-astro-020325-115713 , archivePrefix =. 2506.02101 , primaryClass =

  8. [8]

    , keywords =

    Magnetic Reconnection: An Alternative Explanation of Radio Emission in Galaxy Clusters. , keywords =. doi:10.3847/2041-8213/ad9f2d , archivePrefix =. 2407.11156 , primaryClass =

  9. [9]

    Electromagnetic Phenomena in Cosmical Physics , year = 1958, editor =

    The Neutral Point Theory of Solar Flares. Electromagnetic Phenomena in Cosmical Physics , year = 1958, editor =. doi:10.1017/S0074180900237704 , adsurl =

  10. [10]

    , keywords =

    Sweet's Mechanism for Merging Magnetic Fields in Conducting Fluids. , keywords =. doi:10.1029/JZ062i004p00509 , adsurl =

  11. [11]

    Instability of current sheets and formation of plasmoid chains

    Instability of current sheets and formation of plasmoid chains. Physics of Plasmas , keywords =. doi:10.1063/1.2783986 , archivePrefix =. astro-ph/0703631 , primaryClass =

  12. [12]

    AGU Fall Meeting Abstracts , year = 2009, series =

    Fast Reconnection in Extended High-Lundquist-Number Current Sheets Due to the Plasmoid Instability. AGU Fall Meeting Abstracts , year = 2009, series =

  13. [13]

    Scaling laws of resistive magnetohydrodynamic reconnection in the high-Lundquist-number, plasmoid-unstable regime

    Scaling laws of resistive magnetohydrodynamic reconnection in the high-Lundquist-number, plasmoid-unstable regime. Physics of Plasmas , keywords =. doi:10.1063/1.3420208 , archivePrefix =. 1003.5951 , primaryClass =

  14. [14]

    Fast magnetic reconnection in the plasmoid-dominated regime

    Fast Magnetic Reconnection in the Plasmoid-Dominated Regime. , keywords =. doi:10.1103/PhysRevLett.105.235002 , archivePrefix =. 1008.3330 , primaryClass =

  15. [15]

    Magnetic reconnection and stochastic plasmoid chains in high-Lundquist-number plasmas

    Magnetic reconnection and stochastic plasmoid chains in high-Lundquist-number plasmas. Physics of Plasmas , keywords =. doi:10.1063/1.3703318 , archivePrefix =. 1108.4040 , primaryClass =

  16. [16]

    Reconnection in a Weakly Stochastic Field

    Reconnection in a Weakly Stochastic Field. , keywords =. doi:10.1086/307233 , archivePrefix =. astro-ph/9811037 , primaryClass =

  17. [17]

    Numerical Tests of Fast Reconnection in Weakly Stochastic Magnetic Fields

    Numerical Tests of Fast Reconnection in Weakly Stochastic Magnetic Fields. , keywords =. doi:10.1088/0004-637X/700/1/63 , archivePrefix =. 0903.2052 , primaryClass =

  18. [18]

    Reconnection Studies Under Different Types of Turbulence Driving

    Reconnection studies under different types of turbulence driving. Nonlinear Processes in Geophysics , keywords =. doi:10.5194/npg-19-297-2012 , archivePrefix =. 1203.2971 , primaryClass =

  19. [19]

    , keywords =

    Investigating Turbulence Effects on Magnetic Reconnection Rates through 3D Resistive Magnetohydrodynamic Simulations. , keywords =. doi:10.3847/1538-4357/addc62 , archivePrefix =. 2405.15909 , primaryClass =

  20. [20]

    , keywords =

    The scaling of collisionless, magnetic reconnection for large systems. , keywords =. doi:10.1029/1999GL900481 , adsurl =

  21. [21]

    Physics of Plasmas , keywords =

    The Hall fields and fast magnetic reconnection. Physics of Plasmas , keywords =. doi:10.1063/1.2901194 , adsurl =

  22. [22]

    Plasmoid instability in the semi-collisional regime

    Plasmoid instability in the semi-collisional regime. Journal of Plasma Physics , keywords =. doi:10.1017/S002237781800106X , archivePrefix =. 1804.05145 , primaryClass =

  23. [23]

    , keywords =

    Collisionless Reconnection in an Electron-Positron Plasma. , keywords =. doi:10.1103/PhysRevLett.95.245001 , adsurl =

  24. [24]

    , keywords =

    Reconnection of Quasi-singular Current Sheets: The ``Ideal'' Tearing Mode. , keywords =. doi:10.1088/2041-8205/780/2/L19 , adsurl =

  25. [25]

    General Theory of the Plasmoid Instability

    General theory of the plasmoid instability. Physics of Plasmas , keywords =. doi:10.1063/1.4964481 , archivePrefix =. 1608.04692 , primaryClass =

  26. [26]

    Physics of Fluids , keywords =

    Reconnection rates of magnetic fields including the effects of viscosity. Physics of Fluids , keywords =. doi:10.1063/1.864502 , adsurl =

  27. [27]

    JOSS , keywords =

    The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained. JOSS , keywords =. doi:10.21105/joss.02807 , archivePrefix =. 2009.08231 , primaryClass =

  28. [28]

    Parallel Computing , keywords =

    Scalable communication for high-order stencil computations using CUDA-aware MPI. Parallel Computing , keywords =. doi:10.1016/J.PARCO.2022.102904 , archivePrefix =. 2103.01597 , primaryClass =

  29. [29]

    Extended theory of the Taylor problem in the plasmoid-unstable regime

    Extended theory of the Taylor problem in the plasmoid-unstable regime. Physics of Plasmas , keywords =. doi:10.1063/1.4918331 , archivePrefix =. 1502.06408 , primaryClass =

  30. [30]

    Physical Review X , keywords =

    Tearing Instability and Current-Sheet Disruption in the Turbulent Dynamo. Physical Review X , keywords =. doi:10.1103/PhysRevX.12.041027 , archivePrefix =. 2201.07757 , primaryClass =

  31. [31]

    Physical Review X , keywords =

    Reconnection-Controlled Decay of Magnetohydrodynamic Turbulence and the Role of Invariants. Physical Review X , keywords =. doi:10.1103/PhysRevX.11.041005 , archivePrefix =. 2012.01393 , primaryClass =

  32. [32]

    Nature Communications , keywords =

    Cosmic-void observations reconciled with primordial magnetogenesis. Nature Communications , keywords =. doi:10.1038/s41467-023-43258-3 , archivePrefix =. 2203.03573 , primaryClass =

  33. [33]

    The Open Journal of Astrophysics , keywords =

    Quasi-two-dimensionality of three-dimensional, magnetically dominated, decaying turbulence. The Open Journal of Astrophysics , keywords =. doi:10.33232/001c.122855 , archivePrefix =. 2401.01965 , primaryClass =

  34. [34]

    Physical Review Research , keywords =

    Magnetic island merger as a mechanism for inverse magnetic energy transfer. Physical Review Research , keywords =. doi:10.1103/PhysRevResearch.1.012004 , archivePrefix =. 1901.02448 , primaryClass =

  35. [35]

    , keywords =

    Inverse energy transfer in decaying, three-dimensional, non-helical magnetic turbulence due to magnetic reconnection. , keywords =. doi:10.1093/mnras/staa3849 , archivePrefix =. 2007.07325 , primaryClass =

  36. [36]

    arXiv e-prints , keywords =

    Spectral Evolution and Current Sheet Analysis as Probes of Reconnection-Mediated Decay in Magnetically Dominated Turbulence. arXiv e-prints , keywords =. doi:10.48550/arXiv.2601.17955 , archivePrefix =. 2601.17955 , primaryClass =

  37. [37]

    , keywords =

    Resistively controlled primordial magnetic turbulence decay. , keywords =. doi:10.1051/0004-6361/202449267 , archivePrefix =. 2401.08569 , primaryClass =

  38. [38]

    Journal of Plasma Physics , keywords =

    Hosking integral in non-helical Hall cascade. Journal of Plasma Physics , keywords =. doi:10.1017/S0022377823000028 , archivePrefix =. 2211.14197 , primaryClass =

  39. [39]

    Physical Review Research , keywords =

    Decay law of magnetic turbulence with helicity balanced by chiral fermions. Physical Review Research , keywords =. doi:10.1103/PhysRevResearch.5.L022028 , archivePrefix =. 2302.00512 , primaryClass =

  40. [40]

    Loureiro, N. F. and Schekochihin, A. A. and Uzdensky, D. A. , title =. Physical Review E , publisher =. 2013 , doi =

  41. [41]

    and Kowal, Grzegorz and de Gouveia Dal Pino, Elisabete M

    Vicentin, Giovani H. and Kowal, Grzegorz and de Gouveia Dal Pino, Elisabete M. and Lazarian, Alex , title =. The Astrophysical Journal , publisher =. 2026 , doi =

  42. [42]

    Journal of Plasma Physics , publisher =

    Kumar, Vinay and Bhat, Pallavi , title =. Journal of Plasma Physics , publisher =. 2025 , doi =

  43. [43]

    and Mac Low, Mordecai-Mark and Collins, David C

    Oishi, Jeffrey S. and Mac Low, Mordecai-Mark and Collins, David C. and Tamura, Moeko , title =. The Astrophysical Journal , publisher =. 2015 , doi =

  44. [44]

    Boundary shearing of an initially homogeneous magnetic field

    Heating and activity of the solar corona 1. Boundary shearing of an initially homogeneous magnetic field. , keywords =. doi:10.1029/96JA00428 , adsurl =

  45. [45]

    Journal of Plasma Physics , keywords =

    Scaling of the Hosking integral in decaying magnetically dominated turbulence. Journal of Plasma Physics , keywords =. doi:10.1017/S002237782200109X , archivePrefix =. 2206.07513 , primaryClass =

  46. [46]

    Nonhelical inverse transfer of a decaying turbulent magnetic field

    Nonhelical Inverse Transfer of a Decaying Turbulent Magnetic Field. , keywords =. doi:10.1103/PhysRevLett.114.075001 , archivePrefix =. 1404.2238 , primaryClass =

  47. [47]

    arXiv e-prints , keywords =

    Revision of upper bound on volume-filling intergalactic magnetic fields with LOFAR. arXiv e-prints , keywords =. doi:10.48550/arXiv.2412.14825 , archivePrefix =. 2412.14825 , primaryClass =

  48. [48]

    Formation of Plasmoid Chains in Fusion Relevant Plasmas

    Comisso, L. and Grasso, D. and Waelbroeck, F.L. , title =. J.Phys.Conf.Ser. , year =. doi:10.1088/1742-6596/561/1/012004 , url =. 1412.0933 , archivePrefix =