Magnetic Prandtl number dependence of plasmoid-mediated reconnection

This study demonstrates that while the magnetic Prandtl number significantly influences reconnection rates in the Sweet-Parker regime, this dependence weakens considerably in the fully plasmoid-mediated regime where rates become nearly independent of the Prandtl number, a finding that helps reconcile discrepancies with boundary-driven Taylor problem simulations.

The Gist

Imagine the universe is filled with a super-hot, electrically charged soup called plasma. In this soup, invisible magnetic field lines act like giant rubber bands. Sometimes, these rubber bands get tangled, stretched, and then suddenly snap and reconnect. This snapping process is called magnetic reconnection, and it's the reason for explosive events like solar flares or the aurora borealis. It's how the universe quickly turns stored magnetic energy into heat and motion.

For a long time, scientists thought this snapping happened slowly, like a slow leak in a tire. But we know from looking at the sky that these events happen incredibly fast. To explain this speed, scientists discovered that the "rubber bands" don't just snap in one spot; they break into a chaotic chain reaction of smaller loops and islands, a process called the plasmoid instability. Think of it like a long, thin rope that, when pulled too tight, doesn't just break once, but shatters into a thousand tiny, snapping pieces all at once.

The Big Question: Does "Thickness" Matter?

In this study, the researchers wanted to know if the speed of this snapping depends on how "thick" or "sticky" the plasma is. They used a specific measurement called the Magnetic Prandtl number to describe this stickiness.

• Low Stickiness (Low Prandtl): Imagine the plasma is like water.
• High Stickiness (High Prandtl): Imagine the plasma is like thick honey.

Previous studies suggested that if you make the plasma thicker (more honey-like), the snapping should slow down significantly. It was like saying, "If you try to snap a thick rubber band, it takes much longer than a thin one."

The Experiment: Two Swirling Islands

To test this, the researchers didn't just push on a magnetic field from the outside (which is how previous studies did it). Instead, they set up a simulation where two giant magnetic "islands" naturally swirled together and merged.

Think of it like two whirlpools in a bathtub slowly spinning toward each other. As they merge, the space between them gets squeezed into a thin, stretched-out sheet. This is where the reconnection happens. Because the islands are moving on their own, the "snap" happens spontaneously, just like it does in real space storms, rather than being forced by a human hand.

What They Found

The results were surprising and changed the rules of the game:

1. Before the Snap (The Slow Phase): When the magnetic field wasn't stretched enough to break into pieces, the old rules held true. The thicker the plasma (higher stickiness), the slower the reconnection. It behaved exactly like the "thick rubber band" theory.
2. After the Snap (The Fast Phase): Once the field stretched enough to trigger the "plasmoid instability" (the chain reaction of snapping), the rules changed completely. The speed of the snap stopped caring about the stickiness. Whether the plasma was like water or honey, the reconnection happened at nearly the same fast speed.

The Secret Sauce: The Party of Plasmoids

Why did the stickiness stop mattering? The researchers found that in their "swirling islands" setup, the snapping didn't happen just once. It created a chaotic party of many small magnetic islands (plasmoids) that crashed into each other, merged, and bounced around.

• The Old View: Previous studies looked at the moment just before the chaos really started. They saw the first few snaps and thought, "Okay, stickiness matters here."
• The New View: The researchers watched the full chaos. They saw that the fastest speeds happened when these little islands were crashing into each other and merging. In this wild, non-linear dance, the "stickiness" of the fluid became irrelevant. The sheer violence of the collisions drove the speed, not the fluid's thickness.

Why This Matters

The paper suggests that previous studies might have been looking at the "calm before the storm" rather than the storm itself. In real astrophysical systems (like the space around stars or galaxies), the magnetic fields are constantly swirling and merging on their own, creating this chaotic, high-speed environment.

So, if you want to know how fast energy is released in the universe, you shouldn't worry about how "thick" the plasma is. Once the chaos of merging magnetic islands begins, the universe snaps its magnetic rubber bands at a blistering, consistent speed, regardless of the fluid's texture.

In short: When magnetic fields get really tangled and start breaking into pieces, the speed of the explosion is determined by the chaos of the crash, not by how thick the fluid is.

Technical Summary

Technical Summary: Magnetic Prandtl Number Dependence of Plasmoid-Mediated Reconnection

Problem Statement
Magnetic reconnection is a fundamental mechanism for energy dissipation in plasmas, yet its rate in high-Lundquist number ($S$) astrophysical environments remains a subject of debate. Classical Sweet–Parker theory predicts a reconnection rate scaling as $S^{-1/2}$, which is too slow to explain observed explosive phenomena. While the plasmoid instability is known to enhance reconnection rates, the dependence of this enhanced rate on the magnetic Prandtl number ($Pr_M$) is inconsistent across different simulation setups. Specifically, studies of the boundary-driven Taylor problem (e.g., Comisso et al. 2015, hereafter CGW15) report a scaling of $V_{rec} \propto Pr_M^{-1/2}$ in the plasmoid-mediated regime. Conversely, recent observations of decaying magnetohydrodynamic (MHD) turbulence suggest decay timescales that are largely independent of $Pr_M$. This paper investigates whether the discrepancy arises from the specific boundary conditions and dynamical evolution of the current sheets, specifically comparing externally driven configurations with self-consistently evolving, coalescing magnetic structures.

Methodology
The authors perform two-dimensional direct numerical simulations (DNS) using the compressible isothermal visco-resistive MHD equations solved with the Pencil Code. The study focuses on two distinct configurations:

1. Coalescing Magnetic Islands: The primary setup involves two initially separated magnetic islands that coalesce, forming a current sheet at their interface. This configuration mimics the spontaneous generation of current sheets in decaying turbulence. The initial magnetic vector potential is smoothed to avoid discontinuities, and the density profile is initialized to balance the full Lorentz force (magnetic tension and pressure), ensuring a force-balanced equilibrium.
2. Boundary-Driven Taylor Problem: To facilitate comparison with prior literature, the authors reproduce the setup of CGW15. This involves a stable equilibrium current sheet ($B_y = x$) where reconnection is triggered externally via boundary perturbations.

The simulations span a wide range of Lundquist numbers ($S$) and magnetic Prandtl numbers ($Pr_M$). The reconnection rate ($V_{rec}$) is diagnosed by measuring the rate of magnetic flux depletion within the islands, defined via the time derivative of the vector potential extrema.

Key Contributions and Results

• Sweet–Parker Regime ($S < S_c$): For Lundquist numbers below the critical threshold for plasmoid instability ($S_c \approx 10^5$), the simulations recover the expected Sweet–Parker scaling. The reconnection rate decreases with increasing $Pr_M$, following a scaling of approximately $V_{rec} \propto Pr_M^{-1/4}$, consistent with theoretical predictions by Park et al. (1984) and previous literature.

• Plasmoid-Mediated Regime ($S \gtrsim 10^5$): Once the current sheet becomes unstable to the plasmoid instability, the dependence on $Pr_M$ changes significantly.

  • Coalescing Islands: In the self-consistent coalescing setup, the reconnection rate becomes nearly independent of $Pr_M$ over the explored parameter range ($1 \le Pr_M \le 50$). The authors observe that the highest reconnection rates occur during strongly non-linear phases characterized by the interaction and merger of multiple plasmoids. Quiescent phases, dominated by single small plasmoids, exhibit slower rates. The authors suggest that at astrophysically large $S$, the current sheet would host a dense population of interacting plasmoids, sustaining a persistent, fast, and $Pr_M$-independent reconnection state.
  • Taylor Problem Reproduction: Simulations reproducing the CGW15 boundary-driven setup confirm the previously reported $V_{rec} \propto Pr_M^{-1/2}$ scaling during the plasmoid-dominated phase. However, the authors note that these simulations, like the original CGW15 work, terminate shortly after the onset of plasmoid coalescence due to the development of sharp gradients. Consequently, the measured rates reflect an early plasmoid phase before a sustained non-linear interaction regime is established.

• Resolution of Discrepancy: The paper attributes the differing scalings to the dynamical state of the plasmoid population. The boundary-driven Taylor problem measures reconnection during the linear or early non-linear phase where plasmoids are isolated. In contrast, the coalescing-islands setup naturally evolves into a state of continuous plasmoid interaction and merger. The authors argue that the $Pr_M$-independent scaling is a feature of this strongly non-linear, interaction-dominated regime, whereas the $Pr_M^{-1/2}$ scaling applies to regimes where plasmoid interactions are limited or transient.

• Initial Condition Sensitivity: The authors also address a discrepancy with VKDL26 regarding an intermediate $S^{-1/3}$ scaling regime. They demonstrate that this regime appears only when the initial density profile balances magnetic pressure but neglects magnetic tension (pressure equilibrium only). By using a full force-balanced equilibrium (balancing both tension and pressure), their simulations recover the standard $S^{-1/2}$ scaling up to the onset of plasmoid instability, eliminating the intermediate regime.

Significance and Implications
The authors claim that their results suggest plasmoid-mediated reconnection in self-consistently evolving current sheets exhibits qualitatively different behavior from that predicted in externally driven configurations. Specifically, the strongly non-linear regime dominated by plasmoid interactions and mergers may lead to reconnection rates that are only weakly sensitive to the magnetic Prandtl number.

This finding has potential implications for:

1. Decaying MHD Turbulence: It offers a possible explanation for the $Pr_M$-independent decay timescales observed in recent studies of magnetically dominated turbulence (Brandenburg et al. 2024), suggesting that such systems operate in the strongly non-linear, interaction-dominated regime.
2. Primordial Magnetic Fields: The results may necessitate a re-evaluation of estimates regarding the survivability of primordial magnetic fields in cosmic voids, which have previously relied on the $Pr_M^{-1/2}$ scaling assumption.

The authors maintain modesty regarding these implications, noting that all simulations are two-dimensional and that the parameter range, while extending into the plasmoid regime, remains limited compared to many astrophysical systems. They identify the verification of these scalings in fully non-linear 3D systems as a critical direction for future work.