In incompressible magnetohydrodynamic (MHD) turbulence permeated by a strong uniform guide field B0B_0, the phenomenology of Bigot, Galtier and Politano (2008) predicts anisotropic self-similar energy decay — E⊥∼t−2/3E_\perp \sim t^{-2/3} for shear-Alfvén-wave energy and E∥∼t−1E_\parallel \sim t^{-1} for pseudo-Alfvén-wave energy — from an assumed low-wavenumber spectral exponent and a permanence-of-big-eddies argument. Those authors flagged as "a huge challenge for the future" the derivation of a von Kármán–Howarth relation for strongly anisotropic MHD, from which a Loitsianskii-type invariant could be found, settling whether the widely observed balance law k∥∼k⊥2/3k_\parallel \sim k_\perp^{2/3} persists at the largest scales (outside the inertial range) and fixing the relation between the perpendicular and parallel integral scales during decay. Later work has derived the anisotropic two-point Elsässer correlation equation and shown that self-similar decay requires a constant ratio of parallel to perpendicular similarity lengths, a condition verified in direct numerical simulations; however, no Loitsianskii- or Saffman-type invariant has been derived for the strong-guide-field regime, the anisotropic decay exponents are still not predicted from first principles, and the effect of finite cross-correlation and magnetic helicity on those exponents remains unestablished. The problem is to construct this missing invariant-based decay theory and its consequences, or to refute its existence.

Background

Magnetohydrodynamic (MHD) turbulence couples the velocity field v\mathbf{v} of a conducting fluid to its fluctuating magnetic field b\mathbf{b}; in many astrophysical plasmas (solar corona, solar wind, interstellar medium) a strong large-scale uniform magnetic field B0=B0e^∥\mathbf{B}_0 = B_0 \hat{e}_\parallel is present on top of the fluctuations. A strong guide field makes the turbulence strongly anisotropic: nonlinear energy transfer along B0\mathbf{B}_0 is suppressed, so the cascade proceeds mainly across the field and the spectrum extends to much larger perpendicular than parallel wavenumbers (k⊥≫k∥k_\perp \gg k_\parallel). The linear modes about B0\mathbf{B}_0 are shear-Alfvén waves (velocity and magnetic fluctuations perpendicular to both B0\mathbf{B}_0 and the wavevector) and pseudo-Alfvén waves (the incompressible limit of slow magnetosonic waves, polarized in the plane of B0\mathbf{B}_0 and the wavevector); under strong anisotropy the pseudo-waves are slaved to the shear waves. In the Elsässer variables z±=v±b\mathbf{z}^{\pm} = \mathbf{v} \pm \mathbf{b} the incompressible MHD equations take a symmetric form, and each Elsässer field carries an energy E±E^{\pm}.

For freely decaying turbulence, the classical theory of Navier–Stokes turbulence uses the von Kármán–Howarth equation — the exact evolution equation of the two-point velocity-correlation tensor — together with self-preservation of the correlation shape to derive decay laws; the prediction E(t)∼t−10/7E(t) \sim t^{-10/7} rests on the permanence of big eddies, i.e. on the low-wavenumber spectrum E(k→0)∼ksE(k \to 0) \sim k^s being frozen during decay so that E(t)∼ℓ−(s+1)E(t) \sim \ell^{-(s+1)}, with s=4s=4 tied to conservation of the Loitsianskii integral (a time-independent second moment of the velocity correlation; an analogous Saffman invariant arises for s=2s=2). For isotropic MHD turbulence, Galtier, Politano and Pouquet (1997) derived E(t)∼t−5/6E(t) \sim t^{-5/6} and ℓ(t)∼t1/6\ell(t) \sim t^{1/6} using the Iroshnikov–Kraichnan transfer time τtr=τNL2/τA\tau_{tr} = \tau_{NL}^2/\tau_A (the eddy turnover time squared over the Alfvén time), which accounts for Alfvén-wave propagation along the fluctuating field; magnetic helicity (the volume integral of a⋅b\mathbf{a} \cdot \mathbf{b}, with a\mathbf{a} the vector potential) modifies these laws in the helical case (Biskamp and Müller 1999).

Bigot, Galtier and Politano (2008) extended this phenomenology to the strongly anisotropic, strong-B0B_0 limit, distinguishing the perpendicular and parallel integral scales ℓ⊥\ell_\perp and ℓ∥\ell_\parallel (the correlation lengths of the energy-containing eddies across and along the guide field) and taking ℓ∥\ell_\parallel time-independent because parallel transfers are negligible. With a low-k⊥k_\perp spectrum E⊥(k⊥,k∥)∼k⊥DE_\perp(k_\perp,k_\parallel) \sim k_\perp^{D} (DD the space dimension) at scales beyond ℓ⊥\ell_\perp, permanence of large eddies gives αˉ=βˉ(D+1)\bar{\alpha} = \bar{\beta}(D+1) for exponents E⊥∼t−αˉE_\perp \sim t^{-\bar{\alpha}}, ℓ⊥∼tβˉ\ell_\perp \sim t^{\bar{\beta}}, and the anisotropic transfer-time argument gives 1=αˉ+2βˉ1 = \bar{\alpha} + 2\bar{\beta}, so αˉ=(D+1)/(D+3)\bar{\alpha} = (D+1)/(D+3), βˉ=1/(D+3)\bar{\beta} = 1/(D+3): in 3D, E⊥∼t−2/3E_\perp \sim t^{-2/3} (slower decay than the isotropic t−5/6t^{-5/6}) and ℓ⊥∼t1/6\ell_\perp \sim t^{1/6}. The incompressibility constraint k⊥2E⊥∼k∥2E∥k_\perp^2 E_\perp \sim k_\parallel^2 E_\parallel then yields E∥∼t−1E_\parallel \sim t^{-1} for the slaved pseudo-Alfvén energy. Integrations of the weak-turbulence kinetic equations of Alfvén-wave turbulence (the closed spectral equations governing the resonant interaction of weakly nonlinear Alfvén waves; Galtier, Nazarenko, Newell and Pouquet 2000) recovered these laws, and direct numerical simulations followed them only in the lowest-k∥k_\parallel planes, with finite-box and limited parallel resolution restricting the observable decay ranges. The model further assumes negligible magnetic helicity and weak cross-correlation between the velocity and magnetic fluctuations — i.e. small cross-helicity, the volume integral of v⋅b\mathbf{v} \cdot \mathbf{b}, a state the source paper calls balance turbulence.

A separate empirical fact is the balance law k∥∼k⊥2/3k_\parallel \sim k_\perp^{2/3}: within the inertial range (the range of scales between the large energy-containing eddies and the small dissipation scales, across which cascading energy is transferred), Goldreich–Sridhar critical balance (equality of the linear Alfvén time and the nonlinear turnover time) predicts this relation between the parallel and perpendicular wavenumbers at the edge of the cascade, and it is observed in many simulations (e.g. Cho, Lazarian and Vishniac 2002) and reproduced by heuristic models (Galtier, Pouquet and Mangeney 2005). The decay model above concerns only the largest scales, before the inertial range, so it neither implies nor excludes this law there; whether the 2/32/3 anisotropy persists at the energy-containing scales is precisely the question a Loitsianskii-type invariant for anisotropic MHD would settle. At the time of the source paper, no von Kármán–Howarth relation and no Loitsianskii-type invariant was known for strongly anisotropic MHD — the gap this problem addresses.

Problem Statement

Establish a first-principles theory of self-similar energy decay for incompressible, homogeneous MHD turbulence with a strong uniform guide field B0B_0 — the regime of Bigot, Galtier and Politano (2008) — by determining whether a Loitsianskii-type (or Saffman-type) invariant exists in this regime: a time-independent functional of the low-wavenumber perpendicular and parallel spectra, derivable from the von Kármán–Howarth-type equation for the anisotropic two-point Elsässer correlations (the anisotropic form of which has been derived by Wan et al. 2012 and sequels, so the open challenge is the invariant and its consequences, not the two-point equation itself). The theory must determine: (i) the relation between the perpendicular and parallel integral scales ℓ⊥(t)\ell_\perp(t) and ℓ∥(t)\ell_\parallel(t) during self-similar decay, and in particular whether the balance law k∥∼k⊥2/3k_\parallel \sim k_\perp^{2/3}, well observed inside the inertial range, persists at the largest scales outside the inertial range; and (ii) the anisotropic decay exponents for shear- and pseudo-Alfvén-wave energy, derived rather than assumed, thereby testing the heuristic predictions E⊥∼t−2/3E_\perp \sim t^{-2/3} and E∥∼t−1E_\parallel \sim t^{-1}. How these exponents would be modified under finite cross-correlation (non-balance turbulence) and nonzero magnetic helicity, and how domain size, initial k∥k_\parallel-extent of the spectrum, and dissipation coefficients limit the time ranges over which the asymptotic decay exponents are observable in wave-kinetic and direct numerical calculations, are open questions adjacent to this challenge but are not required outputs of the theory. Extending the decay analysis beyond the incompressible, homogeneous, strong-guide-field limit (compressibility, inhomogeneity, intermediate B0B_0) is an open-ended program rather than part of the completion boundary of this problem.

References

  1. B. Bigot, S. Galtier, H. Politano, 'Energy Decay Laws in Strongly Anisotropic Magnetohydrodynamic Turbulence', Phys. Rev. Lett. 100, 074502 (2008). DOI: 10.1103/PhysRevLett.100.074502, https://doi.org/10.1103/PhysRevLett.100.074502 (arXiv:0712.2024, https://arxiv.org/abs/0712.2024)
  2. B. Bigot, S. Galtier, H. Politano, 'Development of Anisotropy in Incompressible Magnetohydrodynamic Turbulence', Phys. Rev. E 78, 066301 (2008). DOI: 10.1103/PhysRevE.78.066301, https://doi.org/10.1103/PhysRevE.78.066301 (arXiv:0808.3061, https://arxiv.org/abs/0808.3061)
  3. T. de Kármán, L. Howarth, 'On the Statistical Theory of Isotropic Turbulence', Proc. R. Soc. Lond. A 164, 192 (1938). DOI: 10.1098/rspa.1938.0013, https://doi.org/10.1098/rspa.1938.0013
  4. S. Galtier, H. Politano, A. Pouquet, 'Self-Similar Energy Decay in Magnetohydrodynamic Turbulence', Phys. Rev. Lett. 79, 2807 (1997). DOI: 10.1103/PhysRevLett.79.2807, https://doi.org/10.1103/PhysRevLett.79.2807
  5. S. Galtier, S. Nazarenko, A. C. Newell, A. Pouquet, 'A weak turbulence theory for incompressible magnetohydrodynamics', J. Plasma Phys. 63, 447 (2000). DOI: 10.1017/s0022377899008284, https://doi.org/10.1017/s0022377899008284
  6. S. Galtier, A. Pouquet, A. Mangeney, 'On spectral scaling laws for incompressible anisotropic magnetohydrodynamic turbulence', Phys. Plasmas 12, 092310 (2005). DOI: 10.1063/1.2052507, https://doi.org/10.1063/1.2052507
  7. P. Goldreich, S. Sridhar, 'Toward a theory of interstellar turbulence. 2: Strong alfvenic turbulence', Astrophys. J. 438, 763 (1995). DOI: 10.1086/175121, https://doi.org/10.1086/175121
  8. J. Cho, A. Lazarian, E. T. Vishniac, 'Simulations of Magnetohydrodynamic Turbulence in a Strongly Magnetized Medium', Astrophys. J. 564, 291 (2002). DOI: 10.1086/324186, https://doi.org/10.1086/324186
  9. D. Biskamp, W.-C. Müller, 'Decay Laws for Three-Dimensional Magnetohydrodynamic Turbulence', Phys. Rev. Lett. 83, 2195 (1999). DOI: 10.1103/PhysRevLett.83.2195, https://doi.org/10.1103/PhysRevLett.83.2195
  10. J. J. Podesta, 'Laws for third-order moments in homogeneous anisotropic incompressible magnetohydrodynamic turbulence', J. Fluid Mech. 609, 171 (2008). DOI: 10.1017/S0022112008002280, https://doi.org/10.1017/S0022112008002280
  11. S. Galtier, 'Exact vectorial law for axisymmetric magnetohydrodynamics turbulence', Astrophys. J. 704, 1371 (2009). DOI: 10.1088/0004-637X/704/2/1371, https://doi.org/10.1088/0004-637X/704/2/1371
  12. P. A. Davidson, 'On the decay of Saffman turbulence subject to rotation, stratification or an imposed magnetic field', J. Fluid Mech. 663, 268 (2010). DOI: 10.1017/S0022112010003496, https://doi.org/10.1017/S0022112010003496
  13. N. Okamoto, P. A. Davidson, Y. Kaneda, 'On the decay of low-magnetic-Reynolds-number turbulence in an imposed magnetic field', J. Fluid Mech. 651, 295 (2010). DOI: 10.1017/S0022112009994137, https://doi.org/10.1017/S0022112009994137
  14. M. Wan, S. Oughton, S. Servidio, W. H. Matthaeus, 'von Kármán self-preservation hypothesis for magnetohydrodynamic turbulence and its consequences for universality', J. Fluid Mech. 697, 296 (2012). DOI: 10.1017/jfm.2012.61, https://doi.org/10.1017/jfm.2012.61
  15. R. Bandyopadhyay, S. Oughton, M. Wan, W. H. Matthaeus, R. Chhiber, T. N. Parashar, 'Finite Dissipation in Anisotropic Magnetohydrodynamic Turbulence', Phys. Rev. X 8, 041052 (2018). DOI: 10.1103/PhysRevX.8.041052, https://doi.org/10.1103/PhysRevX.8.041052 (arXiv:1812.02379, https://arxiv.org/abs/1812.02379)
  16. R. Bandyopadhyay, W. H. Matthaeus, S. Oughton, M. Wan, 'Evolution of similarity lengths in anisotropic magnetohydrodynamic turbulence', J. Fluid Mech. 876, 5 (2019). DOI: 10.1017/jfm.2019.513, https://doi.org/10.1017/jfm.2019.513 (arXiv:1906.07764, https://arxiv.org/abs/1906.07764)
  17. A. Briard, T. Gomez, 'The decay of isotropic magnetohydrodynamics turbulence and the effects of cross-helicity', J. Plasma Phys. 84, 905840110 (2018). DOI: 10.1017/s0022377818000120, https://doi.org/10.1017/s0022377818000120
  18. D. N. Hosking, A. A. Schekochihin, 'Reconnection-Controlled Decay of Magnetohydrodynamic Turbulence and the Role of Invariants', Phys. Rev. X 11, 041005 (2021). DOI: 10.1103/PhysRevX.11.041005, https://doi.org/10.1103/PhysRevX.11.041005
  19. Y. Yang, M. Wan, W. H. Matthaeus, S. Chen, 'Energy budget in decaying compressible MHD turbulence', J. Fluid Mech. 916, A4 (2021). DOI: 10.1017/jfm.2021.199, https://doi.org/10.1017/jfm.2021.199
  20. A. Salhi, A. Khlifi, R. Marino, F. Feraco, R. Foldes, C. Cambon, 'Waves and non-propagating modes in stratified MHD turbulence subject to a weak mean magnetic field', J. Fluid Mech. 1001, A28 (2024). DOI: 10.1017/jfm.2024.1100, https://doi.org/10.1017/jfm.2024.1100
  21. R. Sasmal, S. Banerjee, 'Universal relations between parallel and perpendicular spectral power law exponents in non-axisymmetric magnetohydrodynamic turbulence' (2025). arXiv:2504.18109, https://arxiv.org/abs/2504.18109
  22. N. P. Müller, S. Galtier, 'Dynamics of fast magnetosonic wave turbulence' (2026). arXiv:2605.22710, https://arxiv.org/abs/2605.22710