Determination of turbulent heating rate and relaxed states in finite Larmor radius magnetohydrodynamic turbulence with helicity barrier

This paper derives exact cascade laws for energy and generalized helicity in finite Larmor radius magnetohydrodynamic turbulence to quantify the ion heating rate arising from the helicity barrier and establishes the corresponding relaxed states characterized by velocity-magnetic field alignment under strong anisotropy.

The Gist

Imagine the Sun's atmosphere (the corona) and the stream of particles blowing away from it (the solar wind) as a giant, chaotic kitchen where the ingredients are not flour and sugar, but super-hot plasma. This plasma is a mix of electrons and ions (heavier charged particles like protons).

For a long time, scientists had a puzzle: Why are the heavy ions in this solar soup getting so hot, specifically in a direction sideways to the Sun's magnetic field? Standard theories of fluid turbulence were like trying to explain a tornado using only a flat map; they couldn't account for the specific "spin" and size of the ions that made them heat up.

This paper introduces a new, more detailed "recipe" called Finite Larmor Radius Magnetohydrodynamics (FLR-MHD). Think of this as upgrading from a blurry, low-resolution photo of the solar wind to a high-definition 3D model that accounts for the actual size of the ions as they spin.

Here is a breakdown of what the authors discovered, using simple analogies:

1. The "Helicity Barrier": A Traffic Jam in Space

In normal fluid turbulence, energy usually flows like water down a waterfall, cascading from big swirls to tiny ripples until it disappears as heat.

However, in this specific solar plasma, the authors found a "traffic jam" caused by something called helicity (a measure of how twisted or knotted the magnetic and velocity fields are).

• The Analogy: Imagine a highway where cars (energy) are trying to drive from a wide open road (large scales) into a narrow tunnel (tiny scales). Suddenly, a massive construction zone (the Helicity Barrier) appears at a specific size.
• The Result: Most of the cars can't get through the construction zone. They pile up right before it. Only a tiny, trickle of cars manages to squeeze through to the other side.

2. The Heating Mechanism: The Pile-Up

Why does this matter for heating?

• Because the energy piles up at this "barrier," the pressure builds up.
• Eventually, this buildup forces the energy to change direction. Instead of just getting smaller, the energy gets squeezed into a very specific, narrow channel that allows it to interact with the ions in a way that heats them up sideways.
• The Paper's Claim: The authors derived a mathematical "receipt" (an exact law) that allows scientists to calculate exactly how much energy is stuck at the barrier versus how much gets through. The difference between these two amounts is the heating rate of the ions. It's like calculating how much fuel is wasted in traffic versus how much actually reaches the destination.

3. No "Steady State": The Unbalanced Scale

In many physics problems, scientists assume a "steady state" where things flow smoothly and evenly.

• The Discovery: The authors found that in this solar plasma, if the flow is unbalanced (one type of wave is much stronger than the other), a steady state is impossible.
• The Analogy: Imagine a seesaw that is heavily weighted on one side. You cannot get it to balance perfectly in the middle. The "Helicity Barrier" prevents the system from ever reaching a calm, steady flow. Instead, the system is constantly shifting, with energy accumulating at the barrier and then releasing in bursts.

4. The "Relaxed" State: When the Chaos Settles

The paper also asks: "If we stop stirring the pot (stop adding energy), how does the plasma finally settle down?"

• The Finding: The plasma doesn't just stop moving. It settles into a specific, organized pattern where the velocity of the particles and the magnetic field lines align with each other.
• The Catch: Because the Sun's magnetic field is so strong and directional (like a long, straight river), the particles cannot twist into a perfect spiral (a "Beltrami" state). Instead, they align in a way that respects the strong magnetic "river," creating a state with a specific pressure gradient.

5. Connecting the Dots: From Big to Small

The authors showed that their new, complex model acts like a universal adapter:

• At large scales (far away from the ions' size), their math simplifies to match the old, well-known theories of solar turbulence.
• At very small scales (inside the ion's spin), it simplifies to match theories about electron behavior.
• In the middle (where the ions live), their new model explains the "missing link" that previous theories couldn't solve.

Summary

This paper provides the mathematical tools to measure exactly how much the Sun's ions are heated by turbulence. It explains that a "traffic jam" of magnetic energy (the helicity barrier) forces energy to pile up and then release in a way that selectively heats heavy ions sideways. This helps solve the mystery of why the solar corona is so hot and why the solar wind accelerates the way it does.

Technical Summary

1. Problem Statement

The paper addresses a critical unsolved problem in space plasma physics: the anomalous heating of the solar corona and the preferential perpendicular heating of heavy ions in the solar wind.

• The Limitation of Existing Models: Standard Magnetohydrodynamics (MHD) and Reduced MHD (RMHD) models are valid only at scales much larger than the ion gyroradius ($\rho_i$). These models predict energy cascades primarily to smaller perpendicular scales, failing to explain the generation of small parallel scales required for Ion Cyclotron Resonance (ICR), which is the mechanism believed to heat ions perpendicularly.
• The Kinetic Gap: While kinetic models (like gyrokinetics) can explain ICR, they are computationally expensive. A hybrid model, Finite Larmor Radius Magnetohydrodynamics (FLR-MHD), bridges this gap by treating electrons as a fluid and ions kinetically. However, a systematic analytical framework to quantify the energy transfer and heating rates within FLR-MHD turbulence was previously missing.
• The Helicity Barrier: Recent numerical studies suggested that in FLR-MHD, a "helicity barrier" forms near the ion gyroradius scale ($k_\perp \rho_i \sim 1$), potentially blocking energy cascades. The paper aims to analytically derive the exact laws governing this process to determine the ion heating rate.

2. Methodology

The authors employ a rigorous analytical approach based on the exact relation theory (Yaglom-type relations) adapted for FLR-MHD.

• Governing Equations: They utilize the FLR-MHD equations derived from gyrokinetics for low-beta plasmas ($\beta \ll 1$) with a strong mean magnetic field ($\mathbf{B}_0$). The system is characterized by two quadratic invariants in the inviscid limit: Total Energy ($E$) and Generalized Helicity ($H$).
• Derivation of Exact Laws:
  • The authors derive exact relations for the two-point correlation functions of energy and generalized helicity in homogeneous turbulence.
  • They derive these laws in two forms:
    1. Divergence Form: Expressing the time evolution of correlations as the divergence of flux terms (analogous to the Kolmogorov 4/5 law).
    2. Alternative (Banerjee-Galtier) Form: A non-divergence form useful for analyzing relaxed states and nonlinear transfer rates.
• Limit Analysis: The derived exact laws are analyzed in two asymptotic limits:
  • Large Scales ($k_\perp \rho_i \ll 1$): Recovering the RMHD limit.
  • Small Scales ($k_\perp \rho_i \gg 1$): Recovering the Electron Reduced MHD (ERMHD) limit.
• Relaxed States: Using the Principle of Vanishing Nonlinear Transfer (PVNLT), they investigate the statistical states the turbulence relaxes to when external forcing is removed.

3. Key Contributions and Results

A. Exact Laws for Energy and Helicity Cascades

The paper provides the first analytical derivation of exact laws for FLR-MHD turbulence:

• Energy Cascade Law: An exact equation relating the time derivative of the energy correlator to the divergence of a flux involving velocity, magnetic field, density, and current fluctuations.
• Generalized Helicity Cascade Law: An exact equation for the generalized helicity correlator.
• Recovery of Known Limits:
  • In the large-scale limit, the derived laws reduce exactly to the known exact laws for RMHD (involving cross-helicity).
  • In the small-scale limit, they reduce to new exact laws for ERMHD (involving magnetic helicity).

B. The Helicity Barrier and Absence of Global Stationarity

A major theoretical finding is the impossibility of a global stationary cascade in strongly imbalanced FLR-MHD turbulence.

• The Contradiction: If a global stationary state existed, the ratio of helicity injection to energy injection ($\sigma = \varepsilon_H / \varepsilon_E$) would need to be constant across all scales. However, the derived exact laws show that this ratio is scale-dependent (constant only in the small-scale limit, not the large-scale limit).
• The Barrier Mechanism: This contradiction confirms the existence of a helicity barrier at $k_\perp \rho_i \sim 1$. The counter-cascading nature of cross-helicity (forward) and magnetic helicity (inverse) creates a bottleneck.
  • Most energy accumulates at the barrier.
  • Only a small, balanced fraction of energy ($\varepsilon_1 - \varepsilon_2$) passes through to smaller scales.
  • This accumulation drives the parallel wavenumber ($k_\parallel$) up until it satisfies the ICR condition ($k_\parallel v_A \sim \Omega_i$).

C. Determination of Ion Heating Rate

The authors propose a method to calculate the ion heating rate ($\varepsilon_{heat}$) based on the segment-wise stationary states:

• Assuming a dissipation mechanism (like ICR) exists at the barrier, the system can be treated as having two separate stationary cascades: one for $k_\perp \rho_i < 1$ (rate $\varepsilon_1$) and one for $k_\perp \rho_i > 1$ (rate $\varepsilon_2$).
• The heating rate is the difference: $\varepsilon_{heat} = \varepsilon_1 - \varepsilon_2$.
• This allows the heating rate to be calculated directly from two-point fluctuations of plasma variables (density, potential, magnetic field) using the derived exact laws, without needing full kinetic simulations.

D. Relaxed States

Using the PVNLT framework, the authors determined the relaxed states of FLR-MHD turbulence:

• Alignment: The relaxed states exhibit an alignment between velocity fluctuations ($\mathbf{u}_\perp$) and magnetic field fluctuations ($\mathbf{b}_\perp$).
• No Beltrami Alignment: Unlike standard MHD, Beltrami alignment (where $\nabla \times \mathbf{u} \propto \mathbf{u}$) is impossible due to the strong anisotropy of the system. The relaxed state is characterized by a finite pressure gradient and specific alignment conditions involving the vorticity and current density, rather than a simple linear relationship between fields.
• Limit Consistency: These relaxed states correctly reduce to the known relaxed states of RMHD (large scales) and ERMHD (small scales).

4. Significance and Impact

• Analytical Validation: The paper provides the first rigorous analytical proof of the "helicity barrier" mechanism, moving beyond previous numerical approximations.
• Quantifying Solar Heating: It offers a practical tool for space physicists to estimate ion heating rates in the solar wind and corona using in-situ spacecraft data (e.g., Parker Solar Probe, Solar Orbiter) by applying the derived exact relations to measured two-point correlations.
• New Theoretical Framework: The derivation of exact laws for ERMHD and the relaxed states of FLR-MHD fills a gap in plasma turbulence theory, providing a bridge between fluid and kinetic descriptions.
• Future Applications: The authors suggest these relations can be used to study intermittency effects and extend the model to include finite electron inertia (KREHM framework).

In summary, this work establishes the theoretical foundation for understanding how turbulent energy is transferred across the ion gyroscale in space plasmas, explaining the mechanism behind perpendicular ion heating via the helicity barrier, and providing exact formulas to quantify this heating.