Conservation of magnetic-helicity fluctuations due to spatial decorrelation of fluxes in decaying MHD turbulence

This paper develops a theoretical framework to explain the conservation of magnetic-helicity fluctuations in decaying MHD turbulence by analyzing boundary terms associated with long-range spatial correlations, finding that such correlations are dynamically suppressed in local and most non-local gauges (including Coulomb) but can arise in specific non-local gauges, a result verified by high-resolution numerical simulations.

The Gist

The Big Picture: The "Ghost" of Magnetic Twists

Imagine you have a giant, invisible bowl of spaghetti (the universe) filled with magnetic noodles. These noodles are constantly tangling, twisting, and snapping apart. This is Magnetohydrodynamic (MHD) turbulence.

In physics, there's a rule called Magnetic Helicity. Think of this as a measure of how "knotted" or "twisted" the magnetic noodles are. If you have a lot of twists, you have high helicity.

For a long time, scientists knew that if you start with a bowl of already knotted noodles, the total amount of knots stays the same as the energy fades away. The noodles just get bigger and looser, but the total "twistiness" is preserved.

But what if you start with no knots at all? What if the magnetic field is a mess of random loops that cancel each other out perfectly? In 2021, a theory was proposed saying that even in this "knot-free" chaos, a specific type of "twistiness fluctuation" (let's call it the Hosking Integral) should be conserved.

The Problem:
The authors of this paper asked a tricky question: Is this conservation law actually real, or is it an illusion caused by how we measure things?

To measure these magnetic fields, physicists have to choose a "gauge." Think of a gauge like choosing a map projection.

• If you use a Mercator projection, Greenland looks huge.
• If you use a different projection, it looks small.
• The actual land doesn't change, but your measurement of its size does.

The authors were worried that maybe the "conservation" of this twistiness was just an artifact of the map (gauge) they were using, rather than a real physical law. They wanted to know: Does the "twistiness" stay constant no matter how we look at it, or does it leak away?

The Investigation: The "Long-Distance Phone Call"

To solve this, the authors had to look at how magnetic fields talk to each other over long distances.

The Analogy of the Party:
Imagine a crowded party (the turbulence).

1. Local Interaction: People talk to the person standing right next to them. This is easy to understand.
2. Long-Distance Interaction: In a normal room, if someone shouts, the sound travels. But in this magnetic "party," there are two weird ways people can communicate across the room:
   • The Pressure Wave (The Shout): If someone pushes a table, the whole room vibrates instantly. This is like the fluid pressure in space. It connects distant points.
   • The Magic Map (The Gauge): If you choose a specific way to draw the map (the gauge), it might force the "coordinates" of two people on opposite sides of the room to be mathematically linked, even if they aren't talking.

The authors wanted to know: Do these long-distance connections get strong enough to break the conservation law?

The Findings: The "No-Leak" Bucket

The team used a sophisticated mathematical method (borrowed from a 1956 paper on water turbulence) to calculate exactly how these long-distance connections behave.

1. The Good News (Most Maps are Safe):
They found that for almost all the standard "maps" (gauges) physicists use—including the most popular one called the Coulomb Gauge—the long-distance connections are too weak to break the rule.

• The Result: The "twistiness fluctuations" stay constant. The bucket doesn't leak.
• The Metaphor: Imagine you are trying to drain a bucket of water by poking a hole in the side. The authors proved that for most standard holes (gauges), the hole is so tiny that the water level (conservation) stays effectively the same.

2. The Bad News (One Weird Map):
However, they discovered a very strange, exotic "map" (a specific non-local gauge choice) where the long-distance connections do become strong enough to break the rule.

• The Result: If you use this weird map, the "twistiness" seems to leak away.
• The Metaphor: This is like finding a specific, bizarre angle to look at the bucket where it suddenly looks like it has a giant hole in it. The authors suspect this isn't a real physical leak, but rather a "ghost" created by the weird way this specific map forces distant points to talk to each other.

3. The Proof (The Simulation):
To be sure, they ran a massive computer simulation (a digital universe) with a resolution of over 12 billion grid points. They measured the "leakage" directly.

• The Result: In the real simulation (using the standard Coulomb gauge), there was no leakage. The correlation functions (the math describing the long-distance talk) died out fast enough to keep the conservation law intact.

Why Does This Matter?

This paper is like a quality control check for the laws of physics.

• For Astronomers: It confirms that we can trust the rules used to predict how magnetic fields in the early universe, or in stars, evolve over time. We know that even in a chaotic, knot-free environment, certain statistical properties remain stable.
• For Mathematicians: It settles a debate about whether these conservation laws depend on the "lens" (gauge) we use to view them. The answer is: No, not for the lenses we actually use.

The Takeaway

The authors effectively said: "We checked the math, we checked the physics, and we ran the simulation. The 'Hosking Integral' is a real, robust conservation law for decaying magnetic turbulence, provided you aren't using a bizarre, unnatural way of measuring it. The universe keeps its 'twistiness' fluctuations safe, even when the chaos is at its peak."

Technical Summary

1. Problem Statement

The paper addresses a fundamental question in the theory of decaying magnetohydrodynamic (MHD) turbulence: Is the integral of the two-point correlation function of magnetic helicity density ($I_H$) conserved in statistically isotropic, non-helical turbulence?

• Context: Previous work by Hosking & Schekochihin (2021) proposed that in non-helical MHD turbulence (where global magnetic helicity $H=0$), the quantity $I_H = \int \langle h(\mathbf{x})h(\mathbf{x}+\mathbf{r})\rangle d^3r$ (where $h = \mathbf{A} \cdot \mathbf{B}$) is conserved. This conservation leads to specific decay laws for magnetic energy ($E_M \propto t^{-10/9}$) and integral scale ($\xi_M \propto t^{4/9}$).
• The Core Issue: The formal conservation of $I_H$ depends on the vanishing of a boundary term, $C_\infty$, in the evolution equation for $I_H$. This term represents the correlation between magnetic helicity fluxes and helicity density at infinite separation.
  • If $C_\infty \neq 0$, long-range correlations could develop, violating conservation and altering decay laws.
  • Two potential mechanisms for such long-range correlations exist:
    1. Non-local pressure forces in incompressible flow.
    2. Non-local gauge choices for the vector potential $\mathbf{A}$ (specifically, gauges where the gauge function $\zeta$ depends on spatial integrals of the fields).
• Objective: To determine from first principles whether $C_\infty$ vanishes for physically relevant gauges and to verify this via high-resolution numerical simulations.

2. Methodology

The authors employ a dual approach combining analytical asymptotic theory and direct numerical simulation (DNS).

A. Analytical Theory (Batchelor & Proudman Approach)

The authors adapt the methodology of Batchelor & Proudman (1956), originally developed for hydrodynamic turbulence, to MHD.

• Initial Conditions: They assume an initial state where distant points are statistically independent (cumulants decay faster than any power law).
• Taylor Series Expansion: They analyze the time evolution of correlation functions (specifically the fifth-order correlations appearing in $C_\infty$) by expanding them in a Taylor series around $t=0$.
• Asymptotic Analysis: They identify the dominant terms at large separations ($r \to \infty$) by examining how non-local effects (pressure and gauge terms) propagate.
  • They evaluate the decay rates of correlation functions based on the Green's function solutions for pressure ($\nabla^2 P = \dots$) and gauge functions ($\nabla^2 \zeta = \phi$).
  • They rigorously check for symmetry cancellations (parity, isotropy) that might force specific terms to zero.

B. Numerical Simulation

• Code: Modified FLASH code solving 3D compressible isothermal MHD equations.
• Resolution: High-resolution simulation with $2304^3$ grid cells.
• Parameters:
  • Initial magnetic field: Non-helical Gaussian random field with spectrum $E_M(k) \propto k^4$.
  • Initial velocity: Zero.
  • Magnetic Prandtl number $P_m = 1$.
  • Initial Lundquist number $S_0 \approx 6 \times 10^4$.
• Gauge: The Coulomb gauge ($\nabla \cdot \mathbf{A} = 0$) is used, which is a Poisson-type non-local gauge.
• Measurement: Direct calculation of the three correlation functions contributing to $C_\infty$ at various times during the decay.

3. Key Contributions and Theoretical Findings

1. Conservation for Local and Poisson-Type Gauges

The authors prove that $C_\infty = 0$ (and thus $I_H$ is conserved) for a broad class of gauges:

• Local Gauges: Where $\zeta$ is a local function of $\mathbf{u}, \mathbf{B}, \mathbf{A}$.
• Poisson-Type Non-Local Gauges: Where $\zeta$ satisfies $\nabla^2 \zeta = \phi$, with $\phi$ being a local function of $\mathbf{u}$ and $\mathbf{B}$.
  • Mechanism: The analysis shows that terms involving a single appearance of the non-local pressure or gauge function decay as $r^{-4}$ or faster. However, due to symmetry arguments (specifically parity and reflection symmetry in isotropic turbulence), the leading-order terms that would decay as $r^{-4}$ vanish identically.
  • Result: The surviving terms decay as $r^{-8}$ or faster. Since $C_\infty$ requires a decay slower than $r^{-2}$ to be non-zero, these terms vanish, ensuring conservation.
  • Significance: This includes the Coulomb gauge, the most common choice in MHD simulations.

2. Identification of a Non-Conserving Gauge Class

The authors identify a specific, "exotic" class of non-local gauges where conservation may fail:

• Definition: Gauges defined by $\zeta = B_i Z_i$ where $\nabla^2 Z_i = \phi B_i$ (with $\phi$ a local function).
• Result: In this specific gauge, the correlation function $\langle \zeta B_i A'_l B'_l \rangle$ decays as $r^{-2}$. This implies $C_\infty$ is finite, and $I_H$ is not conserved.
• Interpretation: This suggests that $I_H$ conservation is not an absolute invariant of the physics but depends on the gauge choice if the gauge allows for specific non-local correlations that do not respect the symmetry cancellations found in standard gauges.

3. Resolution of the Pressure vs. Gauge Debate

The paper clarifies that while non-local pressure forces (hydrodynamic effect) and non-local gauge choices (electromagnetic effect) can theoretically induce long-range correlations, symmetry constraints in isotropic turbulence suppress these correlations for all standard gauges.

4. Numerical Results

The DNS results provide strong empirical support for the theoretical predictions:

• Decay Laws: The simulation confirms the decay laws $E_M \propto t^{-1.02}$ and $\xi_M \propto t^{0.43}$, which are consistent with the theoretical prediction of $I_H$ conservation ($E_M \propto t^{-10/9}$, $\xi_M \propto t^{4/9}$).
• Correlation Functions: The measured correlation functions contributing to $C_\infty$ in the Coulomb gauge were analyzed.
  • They decay rapidly in the "decorrelation range" ($r \gtrsim \xi_M$).
  • The decay is observed to be faster than $r^{-2}$ (consistent with $C_\infty = 0$).
  • Discrepancy Note: The observed decay ($r^{-3}$ to $r^{-4}$) is slightly slower than the theoretical prediction ($r^{-8}$). The authors attribute this to finite-box effects where parity-dependent quantities (like $\langle \mathbf{u} \cdot \mathbf{B} \rangle$) are not exactly zero in a periodic box, allowing sub-dominant terms to persist. However, the decay is still sufficiently fast to ensure $C_\infty \approx 0$.
• Self-Similarity: When normalized by characteristic scales, the correlation functions collapse onto a single curve, confirming self-similar decay.

5. Significance and Conclusion

• Validation of Selective Decay: The paper provides the first rigorous first-principles justification for the conservation of $I_H$ in non-helical MHD turbulence, validating the "Hosking integral" theory.
• Gauge Invariance Clarification: It establishes that while $I_H$ is formally gauge-dependent, its conservation is robust for all physically standard gauges (including Coulomb). The conservation breaks down only for highly specific, non-standard gauges that artificially correlate distant points in a way that violates symmetry constraints.
• Astrophysical Implications: The results reinforce the predictive power of the $t^{-10/9}$ decay law for non-helical magnetic fields in astrophysical contexts (e.g., the intergalactic medium, early universe magnetogenesis), where magnetic helicity is globally zero but fluctuations drive the dynamics.
• Methodological Advance: The successful application of Batchelor & Proudman's asymptotic expansion techniques to MHD turbulence resolves a long-standing ambiguity regarding the role of non-local interactions in turbulent decay.

In summary, the paper demonstrates that spatial decorrelation of fluxes ensures the conservation of magnetic helicity fluctuations in standard decaying MHD turbulence, provided the gauge choice does not introduce artificial long-range correlations that violate the symmetry of the system.