Chaotic/turbulent cross-helical MHD dynamo: from laboratory to the Sun

Basic Information

• Paper ID: 2507.17621
• Title: Chaotic/turbulent cross-helical MHD dynamo: from laboratory to the Sun
• Author: A. Bershadskii (ICAR, Jerusalem, Israel)
• Classification: physics.flu-dyn
• Publication Date: October 15, 2025 (arXiv v4)
• Paper Link: https://arxiv.org/abs/2507.17621

Abstract

This study demonstrates through laboratory experiments, direct numerical simulations, and observations of solar full-disk magnetic fields and sunspot number dynamics that cross-helicity can dominate the magnetic field frequency power spectrum generated by magnetohydrodynamic (MHD) dynamos in chaotic/turbulent rotating flows. The theoretical analysis is based on Kolmogorov-type phenomenology within the framework of distributed chaos concepts. The research reveals that the solar full-disk magnetic field during periods of weak magnetic activity in the past two solar cycles exhibits deterministic chaotic behavior concentrated near the equator. Simultaneously, observational evidence indicates that during periods of strong solar magnetic activity, there exist two types of smooth chaotic (but non-deterministic) cross-helical dynamo states (high-frequency and low-frequency) separated by the 11-year phenomenon.

Research Background and Motivation

Problem Background

The primary source of magnetohydrodynamic (MHD) dynamos is the conversion of kinetic energy from the motion of conducting fluids (plasma) into magnetic energy. The cross-helicity density $h_{cr} = \langle \mathbf{u} \cdot \mathbf{b} \rangle$ (a fundamental invariant connecting the velocity field u and magnetic field b) should play an important role in this process, yet its importance is often underestimated.

Research Motivation

1. Theoretical Deficiency: Existing MHD dynamo theories often neglect the contribution of cross-helicity, particularly in astrophysical applications
2. Symmetry Breaking: The presence of cross-helicity is closely related to the breaking of global reflection symmetry, which presents a challenge in numerical simulations
3. Multi-scale Phenomena: A unified theoretical framework is needed to explain MHD dynamo phenomena at different scales from the laboratory to the Sun

Research Significance

Establishing a unified theoretical framework based on cross-helicity can explain various MHD dynamo phenomena from laboratory von Kármán flows to solar magnetic fields, which is of significant importance for understanding the mechanism of magnetic field generation in stellar interiors.

Core Contributions

1. Proposed a distributed chaos theory dominated by cross-helicity: Based on Kolmogorov phenomenology, established a theoretical framework describing magnetic field dynamics in MHD dynamos
2. Discovered characteristic power spectrum laws: Demonstrated that cross-helicity-dominated systems possess a stretched exponential spectrum of the form $E(f) \propto \exp^{-(f/f_\beta)^{1/2}}$
3. Revealed chaotic characteristics of solar magnetic fields: Through analysis of solar observational data from different periods, identified two types of states: deterministic chaos and distributed chaos
4. Established unified description from laboratory to Sun: Demonstrated universal laws from von Kármán flow experiments to solar magnetic field observations

Detailed Methodology

Theoretical Framework

Cross-Helicity Definition

For ideal MHD, the cross-helicity density is defined as: $h_{cr} = \langle \mathbf{u} \cdot \mathbf{b} \rangle$

In rotating systems, the generalized cross-helicity is: $\tilde{h}_{cr} = \frac{1}{2}\langle \mathbf{u} \cdot \mathbf{b} \rangle + \boldsymbol{\Omega} \cdot \mathbf{A}(t)$

where $\boldsymbol{\Omega}$ is the global angular velocity and $\mathbf{A}(t)$ is the time-varying component of the magnetic vector potential.

Kolmogorov-type Scaling Relations

Within the cross-helical inertial range, the relationship between the characteristic magnetic field value $b_c$ and characteristic frequency $f_c$ is: $b_c \propto h_{cr}^{1/2} \varepsilon^{-1/2} f_c^{1/2}$

where $\varepsilon$ is the total energy dissipation rate.

Local Symmetry Breaking

Even when the net cross-helicity is zero, the local cross-helicity density in chaotic/turbulent flows can be substantial. The locally sign-definite cross-helicity is: $H^\pm_{cr,j} = \int_{V_j} \mathbf{u}(\mathbf{x},t) \cdot \mathbf{b}(\mathbf{x},t) d\mathbf{r}$

The total sign-definite adiabatic invariant is: $I^\pm = \frac{1}{V} \sum_j H^\pm_{cr,j}$

Distributed Chaos Theory

Basic Concepts

When the characteristic frequency $f_c$ exhibits random fluctuations, through ensemble averaging: $E(f) \propto \int_0^\infty P(f_c) \exp^{-(f/f_c)} df_c$

where $P(f_c)$ is the probability distribution of $f_c$.

Probability Distribution Derivation

Assuming the magnetic field intensity $b_c$ follows a half-normal distribution $P(b_c) \propto \exp^{-(b_c^2/2\sigma^2)}$, combined with the scaling relations, we obtain: $P(f_c) \propto f_c^{-1/2} \exp^{-(f_c/4f_\beta)}$

This is a chi-squared distribution function.

Stretched Exponential Spectrum

Substituting the above probability distribution into the ensemble average yields the characteristic stretched exponential spectrum: $E(f) \propto \exp^{-(f/f_\beta)^{1/2}}$

Experimental Setup

Laboratory Experiments

von Kármán Flow Experiment:

• Working medium: liquid sodium
• Magnetic Reynolds number: $Rm \approx 50$
• Configuration: two coaxial counter-rotating impellers
• Measurement: fluctuations of the axial magnetic field component

Numerical Simulations

Direct Numerical Simulation (DNS):

• Governing equations: standard MHD equations
• Boundary conditions: three-dimensional periodic domain with side length $2\pi$
• Forcing: two Taylor-Green vortices simulating impeller rotation
• Non-dimensional parameters: magnetic diffusivity $\eta$, viscosity $\nu$

Solar Observational Data

Data Sources:

1. Full-disk magnetic field: SOLIS Vector Spectromagnetograph (2003-2017)
2. Historical data: Wilcox Solar Observatory (1975-2003)
3. Sunspot number: SILSO database (1700-2025)

Data Processing:

• Short time series: maximum entropy method for power spectrum calculation
• Long time series: fast Fourier transform (FFT)
• Moving average: smoothing using Carrington rotation period (27 days)

Experimental Results

Laboratory and Numerical Simulation Results

Figure 3 displays the power spectra of magnetic field fluctuations from the von Kármán flow experiment and corresponding DNS. Despite significant differences in boundary conditions and physical parameters (the experiment produces an axial dipole while the simulation produces an equatorial dipole), the properly normalized power spectra are in complete agreement, demonstrating the universality of the stretched exponential spectrum $E(f) \propto \exp^{-(f/f_\beta)^{1/2}}$.

Solar Magnetic Field Observational Results

Weak Magnetic Activity Period (2003-2017)

Figure 4 shows the power spectrum of the solar full-disk magnetic field, exhibiting deterministic chaotic characteristics:

• Exponential spectrum: $E(f) \propto \exp^{-(f/f_c)}$
• Characteristic frequency: $f_c$ corresponds to the stellar rotation period $T_c = 24.47$ days
• Physical mechanism: coherent structures produced by equatorial solar rotation dominate the deterministic chaos of the global magnetic field

Strong Magnetic Activity Period (1975-2003)

Figure 5 displays the power spectrum after Carrington rotation smoothing, presenting distributed chaos characteristics:

• Stretched exponential spectrum: $E(f) \propto \exp^{-(f/f_\beta)^{1/2}}$, $\beta = 1/2$
• Parameters: $f_\beta = 0.0015 \text{ d}^{-1}$

Sunspot Number Analysis

Daily Data Analysis (1880-2003)

Figure 6 shows the power spectrum of sunspot numbers smoothed by Carrington rotation:

• Prominent 11-year period peak
• Stretched exponential spectrum at timescales longer than the rotation period
• Parameters: $f_\beta = 0.0017 \text{ d}^{-1}$

Annual Average Data Analysis (1700-2025)

Figure 7 presents analysis at longer timescales:

• Stretched exponential characteristics at frequencies below the 11-year period
• Parameters: $f_\beta = 0.0033 \text{ y}^{-1}$
• Indicates that the degree of randomization in magnetic field dynamics is imprinted in sunspot number dynamics

Reconstructed Magnetic Field Analysis

Based on phase space reconstruction theory, using the covering system method to reconstruct the solar magnetic field (Figures 8-9):

• Embedding dimension: $D = 3$
• Reconstruction time span: 1750-2005
• Results: similarly exhibit stretched exponential spectrum characteristics

Related Work

MHD Dynamo Theory

Traditional MHD dynamo theory primarily focuses on magnetic helicity and the $\alpha$-effect, while the role of cross-helicity is often overlooked. Recent reviews by Yokoi and others have begun to emphasize the anti-diffusion effects of cross-helicity in plasma physics.

Chaotic Dynamics

The identification of deterministic chaos is typically based on exponential power spectra, while the concept of distributed chaos extends traditional chaos theory by describing smooth non-deterministic systems with random characteristic frequencies through stretched exponential spectra.

Solar Magnetic Field Observations

Recent observations indicate that the global solar magnetic field has continuously weakened over the past several solar cycles, providing an opportunity to study deterministic chaos. Research by Hanasoge and others also supports the view that solar turbulence is suppressed and confined to the equatorial region.

Conclusions and Discussion

Main Conclusions

1. Unified Theoretical Framework: The cross-helicity-dominated MHD dynamo theory can provide a unified description of phenomena at different scales from the laboratory to the Sun
2. Classification of Chaotic States:
   • Weak magnetic activity period: deterministic chaos concentrated in the equatorial region
   • Strong magnetic activity period: distributed chaos with two types of states (high-frequency and low-frequency) separated by the 11-year phenomenon
3. Universal Laws: The universal appearance of the stretched exponential spectrum $E(f) \propto \exp^{-(f/f_\beta)^{1/2}}$ across different systems validates the theory

Physical Mechanisms

The small-scale fluctuation dynamo mechanism within the cross-helical inertial range can explain the observed characteristics of magnetic field dynamics. Spontaneous breaking of local reflection symmetry leads to the formation of clumps with non-zero cross-helicity and magnetic helicity, whose long lifetimes sustain the dynamo process.

Limitations

1. Theoretical Assumptions: Scaling relations based on Kolmogorov phenomenology may require modification in strong magnetic field regimes
2. Observational Constraints: Solar interior magnetic fields cannot be directly observed and must be inferred indirectly from surface magnetic fields and sunspot numbers
3. Numerical Simulations: Limited by computational resources, unable to fully simulate the complex magnetohydrodynamic processes in the solar interior

Future Directions

1. Multi-scale Modeling: Develop MHD dynamo models capable of simultaneously handling global and local effects
2. Observational Verification: Validate theoretical predictions using higher-precision solar magnetic field observational data
3. Experimental Extension: Conduct laboratory verification under higher magnetic Reynolds number conditions
4. Astrophysical Applications: Extend the theory to magnetic field generation mechanisms in other stars and planets

In-Depth Evaluation

Strengths

1. High Originality: First systematic treatment of cross-helicity as the dominant factor in MHD dynamos, providing a new theoretical perspective
2. Cross-scale Unification: Successfully establishes a unified theoretical framework spanning from millimeter scales in the laboratory to million-kilometer scales on the Sun
3. Sufficient Experimental Verification: Combines evidence from laboratory experiments, numerical simulations, and solar observations at three levels with compelling arguments
4. Mathematical Rigor: The theoretical derivation based on Kolmogorov phenomenology is logically clear with accurate mathematical formulation

Shortcomings

1. Insufficient Mechanism Explanation: The physical mechanism explaining why cross-helicity can dominate the MHD dynamo process requires further deepening
2. Parameter Dependence: The specific relationships between key parameters in the theory (such as $f_\beta$) and system physical parameters remain unclear
3. Limited Predictive Capability: While capable of explaining observed phenomena, the theory's ability to predict future solar activity requires further verification

Impact

1. Theoretical Contribution: Provides a new research direction for MHD dynamo theory, potentially triggering a paradigm shift in the field
2. Application Value: Offers important guidance for solar physics, stellar physics, and planetary magnetic field research
3. Interdisciplinary Significance: Connects multiple disciplines including plasma physics, astrophysics, and nonlinear dynamics

Applicable Scenarios

1. Astrophysics: Stellar interior magnetic field generation, planetary magnetosphere dynamics
2. Plasma Physics: Self-organization phenomena of magnetic fields in laboratory plasmas
3. Geophysics: Improvements and refinements to Earth dynamo theory
4. Engineering Applications: Optimization of liquid metal MHD generator design

References

The paper cites 34 important references covering multiple fields including MHD dynamo theory, chaotic dynamics, and solar physics observations, providing a solid theoretical foundation and observational support for the research.

Overall Assessment: This is an excellent paper of significant theoretical importance and practical value. Through an innovative theoretical framework, it successfully explains MHD dynamo phenomena across multiple scales, opening new directions for research in related fields. Although certain details require further refinement, its scientific contributions and impact are undeniable.