We solve the equations of motion of a one-dimensional planar Heisenberg (or Vaks-Larkin) model consisting of a system of interacting macro-spins aligned along a ring. Each spin has unit length and is described by its angle with respect to the rotational axis. The orientation of the spins can vary in time due to random forcing and spin-spin interaction. We statistically describe the behaviour of the sum of all spins for different parameters. The term "domino model" in the title refers to the interaction among the spins. We compare the model results with geomagnetic field reversals and find strikingly similar behaviour. The aggregate of all spins keeps the same direction for a long time and, once in a while, begins flipping to change the orientation by almost 180 degrees (mimicking a geomagnetic reversal) or to move back to the original direction (mimicking an excursion). Most of the time the spins are aligned or anti-aligned and deviate only slightly with respect to the rotational axis (mimicking the secular variation of the geomagnetic pole with respect to the geographic pole). Reversals are fast compared to the times in between and they occur at random times, both in the model and in the case of the Earth's magnetic field.

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.