Coupled spin models for magnetic variation of planets and stars

Geomagnetism is characterized by intermittent polarity reversals and rapid fluctuations. We have recently proposed a coupled macro-spin model to describe these dynamics based on the idea that the whole dynamo mechanism is described by the coherent interactions of many small dynamo elements. In this paper, we further develop this idea and construct a minimal model for magnetic variations. This simple model naturally yields many of the observed features of geomagnetism: its time evolution, the power spectrum, the frequency distribution of stable polarity periods, etc. This model has coexistent two phases; i.e. the cluster phase which determines the global dipole magnetic moment and the expanded phase which gives random perpetual perturbations that yield intermittent polarity flip of the dipole moment. This model can also describe the synchronization of the spin oscillation. This corresponds to the case of sun and the model well describes the quasi-regular cycles of the solar magnetism. Furthermore, by analyzing the relevant terms of MHD equation based on our model, we have obtained a scaling relation for the magnetism for planets, satellites, sun, and stars. Comparing it with various observations, we can estimate the scale of the macro-spins.

💡 Research Summary

The paper presents a unified framework for describing magnetic field variations in the Earth, other planets and satellites, and the Sun, based on a coupled macro‑spin model. The authors start by identifying the minimal dynamo element in the Earth’s outer core as an inward‑winding Taylor‑cell (vortex column) that can generate magnetic field through its vorticity. Each such element is abstracted as a vector “macro‑spin” representing the direction and strength of the associated magnetic moment.

Two extreme coupling schemes are considered: short‑range coupling (SCS), where only neighboring spins interact, and long‑range coupling (LCS), where every spin interacts with every other with equal strength. The authors argue that electric currents and magnetic fields in the core provide a natural mechanism for global (long‑range) coupling, and therefore focus on the LCS model. The interaction energy is taken as λ s_i·s_j, with λ < 0, analogous to ferromagnetic coupling, implying that the macro‑spins tend to align. This “macroscopic ferromagnetism” is distinct from the microscopic iron spins, which are demagnetized at core temperatures.

The LCS model exhibits coexistence of two phases. In the “cluster phase” the spins are strongly aligned, producing a global dipole moment that reproduces the observed average geomagnetic field. In the “expanded phase” the spins are essentially random, providing a background of stochastic perturbations. When fluctuations in the expanded phase become sufficiently large, they can trigger a reversal of the global dipole, reproducing the intermittent polarity flips seen in the geological record. Numerical simulations of the model generate time series that match key geomagnetic observables: a power‑law (≈ 1/f) power spectrum, a distribution of stable polarity intervals ranging from thousands to millions of years, and realistic reversal rates.

By varying the coupling strength and damping parameters, the model also displays a synchronized oscillatory regime. In this regime all spins oscillate with nearly the same phase, leading to quasi‑periodic behavior. The authors map this regime onto the solar magnetic cycle, showing that the same macro‑spin dynamics can produce a regular ≈ 11‑year cycle when the system is near a synchronization transition. The synchronization mechanism is compared to the Kuramoto model and to the mean‑field Hamiltonian (HMF) model, highlighting that long‑range interactions naturally generate collective oscillations without external forcing.

To connect the abstract spin model with physical quantities, the authors analyze the governing magnetohydrodynamic (MHD) equations for an incompressible, rotating, electrically conducting fluid. They identify the relevant dimensionless numbers (Reynolds, Rossby, Rayleigh, magnetic Reynolds) for the Earth’s core and use them to estimate the characteristic size, magnetic moment, and coupling strength of a macro‑spin. Extending this scaling analysis to other celestial bodies, they derive a simple relation for the surface magnetic field strength:

 B ∝ Ω R (ρ ν)^{‑1/2},

where Ω is the rotation rate, R the radius, ρ the density, and ν the kinematic viscosity. When applied to planets, moons, the Sun, and various stars, the relation reproduces observed magnetic field strengths within an order of magnitude, suggesting that the macro‑spin size is of order tens of kilometers in planetary cores and larger in stellar convection zones.

Overall, the paper demonstrates that a minimal, analytically tractable model based on globally coupled macro‑spins can capture the essential statistical and dynamical features of planetary and stellar magnetism. It provides a bridge between detailed numerical MHD simulations and phenomenological descriptions, offering a new perspective on how large‑scale magnetic fields emerge, fluctuate, and reverse across a wide range of astrophysical objects. The work opens avenues for further quantitative comparison with paleomagnetic records, solar observations, and future high‑resolution dynamo simulations.