Simulations of a Magnetic Fluctuation Driven Large Scale Dynamo and Comparison with a Two-scale Model

Models of large scale (magnetohydrodynamic) dynamos (LSD) which couple large scale field growth to total magnetic helicity evolution best predict the saturation of LSDs seen in simulations. For the simplest so called “{\alpha}2” LSDs in periodic boxes, the electromotive force driving LSD growth depends on the difference between the time-integrated kinetic and current helicity associated with fluctuations. When the system is helically kinetically forced (KF), the growth of the large scale helical field is accompanied by growth of small scale magnetic (and current) helicity which ultimately quench the LSD. Here, using both simulations and theory, we study the complementary magnetically forced(MF) case in which the system is forced with an electric field that supplies magnetic helicity. For this MF case, the kinetic helicity becomes the back-reactor that saturates the LSD. Simulations of both MF and KF cases can be approximately modeled with the same equations of magnetic helicity evolution, but with complementary initial conditions. A key difference between KF and MF cases is that the helical large scale field in the MF case grows with the same sign of injected magnetic helicity, whereas the large and small scale magnetic helicities grow with opposite sign for the KF case. The MF case can arise even when the thermal pressure is approximately smaller than the magnetic pressure, and requires only that helical small scale magnetic fluctuations dominate helical velocity fluctuations in LSD driving. We suggest that LSDs in accretion discs and Babcock models of the solar dynamo are actually MF LSDs.

💡 Research Summary

This paper investigates large‑scale magnetohydrodynamic (MHD) dynamos driven by magnetic fluctuations (MF) and compares them with the more familiar kinetically forced (KF) α² dynamos. In the traditional KF scenario, an externally imposed helical velocity field injects kinetic helicity ⟨v·ω⟩, which generates an electromotive force (EMF) ε = ⟨v×b⟩ that amplifies a large‑scale magnetic field. As the dynamo evolves, a small‑scale current helicity ⟨j·b⟩ of opposite sign builds up, reducing the α‑effect (α ∝ ⟨v·ω⟩ − ⟨j·b⟩) and eventually quenching the growth.

The authors introduce the complementary MF case, where a helical electric field is imposed, directly injecting magnetic helicity into the small‑scale magnetic field. In this configuration the small‑scale current helicity is the primary driver of the α‑effect, while kinetic helicity emerges later as a back‑reaction that limits the dynamo. Both cases are described by the same two‑scale helicity evolution equations: the total magnetic helicity H = ⟨A·B⟩ is split into a large‑scale component H₁ and a small‑scale component H₂, and their time derivatives contain the EMF term, resistive dissipation, and an external helicity source Sₕ (present only in the MF case). The α‑coefficient evolves as α = τ(⟨v·ω⟩ − ⟨j·b⟩), with τ the turbulent correlation time. The crucial difference lies in the initial conditions: KF starts with a positive kinetic helicity and zero current helicity, whereas MF starts with a positive current helicity and negligible kinetic helicity. Consequently, the sign relationship between large‑ and small‑scale magnetic helicities is opposite in KF (H₁ = −H₂) and identical in MF (H₁ = +H₂).

Numerical experiments were performed in a periodic cubic domain (128³ grid points) with a forcing wavenumber k_f = 5 and magnetic Prandtl number unity. For KF, a helical body force was applied to the momentum equation; for MF, a helical electric field was added to the induction equation. Both runs began from a quiescent, non‑magnetic state. In the KF run, kinetic helicity quickly drives exponential growth of the large‑scale magnetic energy, while the small‑scale current helicity grows with opposite sign and eventually saturates the dynamo. In the MF run, the imposed electric field creates a rapid increase of small‑scale current helicity, which immediately feeds a positive α‑effect; the large‑scale magnetic field then grows with the same helicity sign as the injected magnetic helicity. As the field strengthens, kinetic helicity is generated non‑linearly and provides the negative feedback that quenches the growth. The time histories of H₁, H₂, ⟨v·ω⟩, ⟨j·b⟩, and ε in both simulations are accurately reproduced by the two‑scale model when the appropriate initial conditions are used.

A key implication of the MF dynamo is that it can operate even when magnetic pressure exceeds thermal pressure (plasma β < 1), provided that the small‑scale magnetic fluctuations dominate over velocity fluctuations in the helicity budget. This situation is plausible in strongly magnetized astrophysical environments such as accretion disc coronae, magnetically dominated jets, or the solar tachocline where magnetic buoyancy and shear can generate helical magnetic structures. The authors therefore suggest that many large‑scale dynamos traditionally interpreted as kinetic‑driven may in fact be magnetically driven, including the Babcock–Leighton solar dynamo and the dynamo action associated with the magnetorotational instability in accretion discs.

In summary, the paper demonstrates that both kinetically forced and magnetically forced α² dynamos obey the same helicity evolution framework, differing only in the sign and magnitude of the initial helicity sources. The MF case reveals a “magnetic‑fluctuation‑driven” dynamo where kinetic helicity acts as the quenching agent, expanding the theoretical landscape of astrophysical dynamo mechanisms and offering new avenues for interpreting observations of large‑scale magnetic fields in high‑β and low‑β plasma regimes.