Oscillatory migrating magnetic fields in helical turbulence in spherical domains

We present direct numerical simulations of the equations of compressible magnetohydrodynamics in a wedge-shaped spherical shell, without shear, but with random helical forcing which has negative (positive) helicity in the northern (southern) hemisphere. We find a large-scale magnetic field that is nearly uniform in the azimuthal direction and approximately antisymmetric about the equator. Furthermore, the large-scale field in each hemisphere oscillates on nearly dynamical time scales with reversals of polarity and equatorward migration. Corresponding mean-field models also show similar migratory oscillations with a frequency that is nearly independent of the magnetic Reynolds number. This mechanism may be relevant for understanding equatorward migration seen in the solar dynamo.

💡 Research Summary

The paper investigates whether a large‑scale, oscillatory magnetic field with equatorward migration can be generated in a spherical domain without any imposed shear, relying solely on helical turbulence. Using direct numerical simulations (DNS) of the compressible magnetohydrodynamic (MHD) equations, the authors model a wedge‑shaped shell that spans a substantial portion of latitude (from roughly 45° to 135°) and a quarter of azimuth. The turbulence is driven by a random, statistically homogeneous forcing that injects kinetic helicity of opposite sign in the two hemispheres: negative helicity in the northern half and positive helicity in the southern half. This setup mimics the sign change of the α‑effect across the solar equator while deliberately omitting any differential rotation (Ω‑effect).

Key aspects of the numerical experiment include: (i) periodic boundary conditions in the azimuthal direction, (ii) perfectly conducting inner radial boundary and an open outer boundary, (iii) a forcing wavenumber that concentrates energy at scales a few times smaller than the shell depth, and (iv) magnetic Reynolds numbers (Re_M) ranging from 50 to 200. The authors monitor the azimuthally averaged magnetic field ⟨B⟩(θ,r,t), compute spherical‑harmonic spectra, and perform time‑frequency analysis to extract oscillation periods and migration speeds. In parallel, they construct a one‑dimensional mean‑field model based on the standard α‑effect and turbulent diffusivity η_t, with α measured directly from the DNS via the test‑field method.

The DNS results reveal a clear emergence of a dipolar‑like large‑scale field that is nearly uniform in azimuth and antisymmetric about the equator. Within each hemisphere the field exhibits quasi‑periodic reversals on a timescale of order the dynamical turnover time (≈0.8–1.2 τ_dyn). During each cycle the magnetic pattern drifts equatorward, with a migration speed of roughly 0.1 R_⊙ per dynamical time. Importantly, the oscillation frequency and migration speed show only a weak dependence on Re_M, indicating that the underlying wave dynamics are governed by the structural properties of the α‑effect rather than by the detailed turbulent intensity.

The mean‑field calculations, using the measured α(θ) profile (negative in the north, positive in the south) and a uniform η_t, reproduce the DNS behaviour remarkably well. The α²‑dynamo wave propagates toward the equator because the product α ∂α/∂θ is negative, a condition that selects equatorward propagation in the absence of shear. The frequency of the wave is set by the magnitude of α and η_t and remains essentially constant as Re_M varies, confirming that nonlinear quenching of α primarily affects the amplitude, not the phase speed, of the dynamo wave.

These findings have several important implications. First, they demonstrate that a pure α²‑dynamo, without any Ω‑effect, can generate a self‑sustained, oscillatory magnetic field that migrates equatorward—a feature traditionally attributed to the combined αΩ mechanism in solar dynamo theory. Second, the robustness of the oscillation period against changes in magnetic Reynolds number suggests that the solar cycle period might be insensitive to the exact level of turbulent diffusivity, provided the sign change of α across the equator is maintained. Third, the study offers a concrete numerical validation of the theoretical prediction that a sign reversal of the α‑effect alone can dictate the direction of dynamo wave propagation.

Nevertheless, the authors acknowledge limitations. The wedge geometry restricts the azimuthal extent, and the imposed forcing spectrum is idealised compared with realistic solar convection. Moreover, the simulations omit the strong radial shear present in the solar tachocline, which likely interacts with the α‑effect in the real Sun. Future work is proposed to gradually introduce differential rotation, explore broader azimuthal domains, and incorporate more realistic stratification and convection.

In summary, the paper provides compelling evidence that helical turbulence with hemispheric sign reversal can drive a large‑scale, oscillatory magnetic field that migrates toward the equator even in the complete absence of shear. This pure α²‑dynamo mechanism offers a fresh perspective on the origin of the solar equatorward migration and may be relevant for other astrophysical objects where shear is weak or absent.