Quasi-linear theory of perpendicular ion heating by critically balanced turbulence

This paper employs quasi-linear theory to analytically derive a unified ion heating rate that smoothly transitions between stochastic and cyclotron-resonant mechanisms depending on the imbalance of Alfvénic turbulence, while explaining the suppression of heating at small amplitudes due to magnetic moment conservation.

The Gist

The Big Picture: Heating Up Space Soup

Imagine the space around our Sun (the solar wind) and the atmosphere above it (the solar corona) as a giant, invisible pot of "plasma soup." This soup is made of charged particles (ions and electrons) and magnetic fields.

Usually, when you heat a pot of soup on a stove, the heat spreads out evenly. But in space, things are different. The "stove" is turbulence—chaotic, swirling motions in the magnetic field. The paper asks a specific question: How does this turbulence heat the ions (the heavy particles in the soup), and why do they get hotter on their sides (perpendicular to the magnetic field) rather than just getting hotter overall?

The authors found that the answer depends on how "balanced" the turbulence is.

The Two Ways to Heat the Soup

The paper describes two main mechanisms for how the magnetic turbulence kicks the ions, making them spin faster and hotter. Think of these as two different ways to push a child on a swing:

1. The "Stochastic" Push (Balanced Turbulence):
   Imagine the swing is being pushed by a crowd of people from both sides (left and right) with equal strength. The pushes are random and chaotic. Sometimes you get a push from the left, sometimes the right. The child doesn't move in a perfect rhythm; they just get jostled around, gaining energy through a "random walk."

   • In the paper: This happens when the turbulence is balanced (equal energy moving with and against the magnetic field). The ions get kicked by random fluctuations, breaking their smooth spinning motion and heating them up.

2. The "Resonant" Push (Imbalanced Turbulence):
   Now imagine the swing is only being pushed by a crowd from one side. The pushes are rhythmic and perfectly timed. If the pusher hits the swing at exactly the right moment in its arc, the swing goes higher and higher very efficiently.

   • In the paper: This happens when the turbulence is imbalanced (mostly energy moving in one direction). The ions "resonate" with the waves, like a swing matching the rhythm of a pusher. This is called cyclotron-resonant heating.

The "Goldilocks" Discovery

The most important finding of this paper is that these two methods aren't actually separate worlds. They are part of a smooth spectrum.

The authors created a mathematical model (a "recipe") that describes the turbulence in space. They found that as you change the balance of the turbulence (from equal pushes to one-sided pushes), the heating mechanism smoothly transitions from the "random jostling" style to the "perfect rhythm" style.

The Universal Formula:
Regardless of whether the turbulence is balanced or imbalanced, the heating rate follows a specific, predictable pattern.

• The Analogy: Think of the turbulence amplitude (how strong the waves are) as the "volume" of the music.
  • If the volume is too low (small waves), the ions don't heat up much because they hold onto their "magnetic moment" (a rule that says they keep spinning smoothly unless the wave is strong enough to break that rule). It's like trying to push a heavy swing with a gentle breeze; nothing happens.
  • Once the volume gets loud enough, the heating kicks in.
  • The paper proves that the heating rate always looks like a specific mathematical curve: it starts very low (suppressed) and then rises sharply as the turbulence gets stronger.

Why This Matters

Before this paper, scientists had different theories for balanced turbulence (stochastic) and imbalanced turbulence (resonant). They treated them as separate problems.

This paper shows that it's all the same physics, just viewed through different lenses.

• The "Imbalance" Knob: The authors show that the "imbalance" of the turbulence (how much more energy is flowing one way than the other) changes the shape of the turbulence's "frequency spectrum" (the range of wave speeds).
• The Result: This change in shape is what switches the heating mechanism from "random jostling" to "perfect rhythm."

The "Suppression" Effect

The paper also explains why ions don't heat up instantly when turbulence is weak.

• The Analogy: Imagine a spinning top. If you tap it gently, it keeps spinning smoothly. It resists the tap. This is the conservation of magnetic moment.
• The paper mathematically proves that for small waves, this "resistance" is very strong, and heating is almost zero. But once the waves get strong enough to overcome this resistance, the heating explodes. The paper provides a precise formula for exactly how this "resistance" fades away as the waves get stronger.

Summary

In short, the authors used advanced math (quasi-linear theory) to show that:

1. Ions in space are heated by magnetic turbulence.
2. Whether the turbulence is balanced or unbalanced, the heating follows a single, universal rule.
3. The mechanism shifts smoothly from "random kicking" to "rhythmic pushing" as the turbulence becomes more one-sided.
4. There is a "threshold" where weak turbulence fails to heat ions because the ions are too "stubborn" (conserving their magnetic moment), but once the turbulence gets loud enough, the heating kicks in efficiently.

This helps scientists understand how the Sun's corona gets so hot and how the solar wind accelerates, providing a single mathematical framework to explain observations that previously seemed contradictory.

Technical Summary

Technical Summary: Quasi-linear theory of perpendicular ion heating by critically balanced turbulence

Problem Statement
In collisionless astrophysical plasmas, such as the solar corona and wind, turbulent energy dissipation is a primary mechanism for heating ions and electrons. Observations indicate that ion heating is often anisotropic, with ions preferentially heated perpendicular to the local magnetic field. While two primary mechanisms have been proposed to explain this—stochastic heating (driven by the breaking of magnetic moment conservation in low-frequency, broad-spectrum turbulence) and cyclotron-resonant heating (driven by resonant interactions with waves)—the transition between these regimes and their behavior under varying levels of turbulence imbalance (the difference in energy between Alfvénic fluctuations propagating parallel and antiparallel to the magnetic field) remains a subject of theoretical investigation. Existing empirical models for stochastic heating lack a formal derivation that accounts for general levels of imbalance, and the connection between stochastic and resonant heating in critically balanced turbulence requires further analytical clarification.

Methodology
The authors employ the framework of quasi-linear theory to analytically compute the perpendicular heating rates of ions interacting with large-scale Alfvénic fluctuations. The approach involves:

1. Quasi-linear Diffusion: Starting from the Vlasov equation, the authors derive diffusion coefficients in velocity space for ions interacting with small-amplitude electromagnetic fluctuations. They simplify these coefficients by assuming purely Alfvénic fluctuations and the large-scale limit ($k_\perp \rho_i \ll 1$), where $\rho_i$ is the ion thermal gyroradius.
2. Turbulence Modeling: To close the system, a novel model for the wavevector-frequency spectrum of Reduced Magnetohydrodynamics (RMHD) turbulence is developed. This model incorporates the concept of critical balance (CB) and explicitly accounts for turbulence imbalance, quantified by the normalized cross-helicity $\sigma_c$. The model describes how the spectrum broadens in frequency due to nonlinear interactions in balanced turbulence ($\sigma_c = 0$) and sharpens along the linear dispersion relation in imbalanced turbulence ($\sigma_c \to 1$).
3. Analytical and Numerical Calculation: The derived diffusion coefficients are integrated over a Maxwellian ion distribution to calculate the perpendicular heating rate $Q_\perp$. The authors perform both analytical derivations for limiting cases (balanced and fully imbalanced) and numerical integrations for general imbalance levels. They also numerically evolve the distribution function to observe the structural changes in velocity space.

Key Contributions

• Unified Heating Framework: The paper provides a single analytical framework that describes ion heating across the entire spectrum of turbulence imbalance, bridging the gap between stochastic heating and cyclotron-resonant heating.
• Model Spectrum: A new model for the RMHD wavevector-frequency spectrum is introduced, which captures the dependence of spectral broadening on the imbalance parameter $\sigma_c$. This model is validated against numerical simulations of RMHD turbulence.
• Analytical Derivation of Suppression: The authors analytically derive the suppression factor for heating at small turbulent amplitudes. In the balanced limit, this derivation recovers the functional form of the empirical stochastic heating formula (Chandran et al. 2010a), providing a theoretical basis for the exponential suppression related to magnetic moment conservation.
• Mechanism Transition: The work demonstrates how the heating mechanism transitions smoothly: in balanced turbulence, the broadened spectrum allows for stochastic-like heating across a wide range of frequencies; in imbalanced turbulence, the spectrum narrows, shifting the mechanism toward cyclotron-resonant heating where ions interact with specific resonant frequencies.

Results

• General Heating Rate Form: The perpendicular heating rate is found to follow a general form regardless of the imbalance level:
  $$Q_\perp \propto \xi_{\rho,th}^3 F(\xi_{\rho,th}; \sigma_c)$$
  where $\xi_{\rho,th}$ is the amplitude of $\rho_i$-scale velocity fluctuations and $F$ is an imbalance-dependent suppression factor that vanishes as $\xi_{\rho,th} \to 0$.
• Balanced Limit ($\sigma_c = 0$): The analytical result for balanced turbulence yields a suppression factor closely resembling the empirical exponential form $e^{-c_2/\xi_{\rho,th}}$. This confirms that the suppression arises fundamentally from the temporal correlations of the turbulence and the broadening of the spectrum, rather than solely from the ion distribution function.
• Imbalanced Limit ($\sigma_c \to 1$): In the fully imbalanced limit, the heating rate is governed by resonant interactions. The suppression at small amplitudes is sharper due to the requirement that ions resonate with a single frequency mode.
• Scaling Behavior: The heating rate exhibits a transition in scaling with turbulent amplitude $\xi_{\rho,th}$. At large amplitudes, $Q_\perp \propto \xi_{\rho,th}^3$. At small amplitudes, the rate scales as $\xi_{\rho,th}^6$ due to interactions with weak fluctuations above the critical balance cone, though this contribution is suppressed when considering only modes below the cone.
• Velocity Space Evolution: Numerical evolution of the distribution function confirms that balanced turbulence produces a "flat-topped" distribution characteristic of stochastic heating, while imbalanced turbulence leads to diffusion along contours of constant energy in the wave frame, characteristic of resonant heating.

Significance
The paper claims to consolidate the qualitative understanding of ion heating in astrophysical plasmas by unifying stochastic and cyclotron-resonant heating mechanisms within a single quasi-linear formalism. By recovering the empirical stochastic heating formula from first principles in the balanced limit and extending the theory to general imbalance, the work yields specific quantitative predictions. These predictions are intended to aid in the analysis of numerical simulations and spacecraft observations (e.g., Parker Solar Probe data) to better constrain the partitioning of turbulent energy and the evolution of plasma thermal structure in the heliosphere and beyond. The authors note that while their model neglects certain small-scale kinetic effects (such as the transition to ion-cyclotron waves at $k_\parallel d_i \sim 1$ and the helicity barrier), it provides a robust baseline for understanding large-scale Alfvénic heating.