# Dynamo theory

Illustration of the dynamo mechanism that creates the Earth's magnetic field: convection currents of fluid metal in the Earth's outer core, driven by heat flow from the inner core, organized into rolls by the Coriolis force, create circulating electric currents, which generate the magnetic field.[1]

In physics, the dynamo theory proposes a mechanism by which a celestial body such as Earth or a star generates a magnetic field. The dynamo theory describes the process through which a rotating, convecting, and electrically conducting fluid can maintain a magnetic field over astronomical time scales. A dynamo is thought to be the source of the Earth's magnetic field and the magnetic fields of Mercury and the Jovian planets.

## History of theory

When William Gilbert published de Magnete in 1600, he concluded that the Earth is magnetic and proposed the first hypothesis for the origin of this magnetism: permanent magnetism such as that found in lodestone. In 1919, Joseph Larmor proposed that a dynamo might be generating the field.[2][3] However, even after he advanced his hypothesis, some prominent scientists advanced alternative explanations. Einstein believed that there might be an asymmetry between the charges of the electron and proton so that the Earth's magnetic field would be produced by the entire Earth. The Nobel Prize winner Patrick Blackett did a series of experiments looking for a fundamental relation between angular momentum and magnetic moment, but found none.[4][5]

Walter M. Elsasser, considered a "father" of the presently accepted dynamo theory as an explanation of the Earth's magnetism, proposed that this magnetic field resulted from electric currents induced in the fluid outer core of the Earth. He revealed the history of the Earth's magnetic field through pioneering the study of the magnetic orientation of minerals in rocks.

In order to maintain the magnetic field against ohmic decay (which would occur for the dipole field in 20,000 years), the outer core must be convecting. The convection is likely some combination of thermal and compositional convection. The mantle controls the rate at which heat is extracted from the core. Heat sources include gravitational energy released by the compression of the core, gravitational energy released by the rejection of light elements (probably sulfur, oxygen, or silicon) at the inner core boundary as it grows, latent heat of crystallization at the inner core boundary, and radioactivity of potassium, uranium and thorium.[6]

At the dawn of the 21st century, numerical modeling of the Earth's magnetic field has not been successfully demonstrated, but appears to be in reach. Initial models are focused on field generation by convection in the planet's fluid outer core. It was possible to show the generation of a strong, Earth-like field when the model assumed a uniform core-surface temperature and exceptionally high viscosities for the core fluid. Computations which incorporated more realistic parameter values yielded magnetic fields that were less Earth-like, but also point the way to model refinements which may ultimately lead to an accurate analytic model. Slight variations in the core-surface temperature, in the range of a few millikelvins, result in significant increases in convective flow and produce more realistic magnetic fields.[7][8]

## Formal definition

Dynamo theory describes the process through which a rotating, convecting, and electrically conducting fluid acts to maintain a magnetic field. This theory is used to explain the presence of anomalously long-lived magnetic fields in astrophysical bodies. The conductive fluid in the geodynamo is liquid iron in the outer core, and in the solar dynamo is ionized gas at the tachocline. Dynamo theory of astrophysical bodies uses magnetohydrodynamic equations to investigate how the fluid can continuously regenerate the magnetic field.[9]

It was once believed that the dipole, which comprises much of the Earth's magnetic field and is misaligned along the rotation axis by 11.3 degrees, was caused by permanent magnetization of the materials in the earth. This means that dynamo theory was originally used to explain the Sun's magnetic field in its relationship with that of the Earth. However, this hypothesis, which was initially proposed by Joseph Larmor in 1919, has been modified due to extensive studies of magnetic secular variation, paleomagnetism (including polarity reversals), seismology, and the solar system's abundance of elements. Also, the application of the theories of Carl Friedrich Gauss to magnetic observations showed that Earth's magnetic field had an internal, rather than external, origin.

There are three requisites for a dynamo to operate:

• An electrically conductive fluid medium
• Kinetic energy provided by planetary rotation
• An internal energy source to drive convective motions within the fluid.[10]

In the case of the Earth, the magnetic field is induced and constantly maintained by the convection of liquid iron in the outer core. A requirement for the induction of field is a rotating fluid. Rotation in the outer core is supplied by the Coriolis effect caused by the rotation of the Earth. The Coriolis force tends to organize fluid motions and electric currents into columns (also see Taylor columns) aligned with the rotation axis. Induction or creation of magnetic field is described by the induction equation:

${\displaystyle {\frac {\partial \mathbf {B} }{\partial t}}=\eta \nabla ^{2}\mathbf {B} +\nabla \times (\mathbf {u} \times \mathbf {B} )}$

where u is velocity, B is magnetic field, t is time, and ${\displaystyle \eta =1/(\sigma \mu )}$ is the magnetic diffusivity with ${\displaystyle \sigma }$ electrical conductivity and ${\displaystyle \mu }$ permeability. The ratio of the second term on the right hand side to the first term gives the magnetic Reynolds number, a dimensionless ratio of advection of magnetic field to diffusion.

### Tidal heating supporting a dynamo

Tidal forces between celestial orbiting bodies cause friction that heats up their interiors. This is known as tidal heating, and it helps keep the interior in a liquid state. A liquid interior that can conduct electricity is required to produce a dynamo. Saturn's Enceladus and Jupiter's Io have enough tidal heating to liquify their inner cores, but they may not create a dynamo because they cannot conduct electricity.[11][12] Mercury, despite its small size, has a magnetic field, because it has a conductive liquid core created by its iron composition and friction resulting from its highly elliptical orbit.[13] It is theorized that the Moon once had a magnetic field, based on evidence from magnetized lunar rocks, due to its short-lived closer distance to Earth creating tidal heating.[14] An orbit and rotation of a planet helps provide a liquid core, and supplements kinetic energy that supports a dynamo action.

## Kinematic dynamo theory

In kinematic dynamo theory the velocity field is prescribed, instead of being a dynamic variable. This method cannot provide the time variable behavior of a fully nonlinear chaotic dynamo but is useful in studying how magnetic field strength varies with the flow structure and speed.

Using Maxwell's equations simultaneously with the curl of Ohm's law, one can derive what is basically the linear eigenvalue equation for magnetic fields (B) which can be done when assuming that the magnetic field is independent from the velocity field. One arrives at a critical magnetic Reynolds number above which the flow strength is sufficient to amplify the imposed magnetic field, and below which it decays.

The most functional feature of kinematic dynamo theory is that it can be used to test whether a velocity field is or is not capable of dynamo action. By applying a certain velocity field to a small magnetic field, it can be determined through observation whether the magnetic field tends to grow or not in reaction to the applied flow. If the magnetic field does grow, then the system is either capable of dynamo action or is a dynamo, but if the magnetic field does not grow, then it is simply referred to as non-dynamo.

The membrane paradigm is a way of looking at black holes that allows for the material near their surfaces to be expressed in the language of dynamo theory.

### As a spontaneous breakdown of topological supersymmetry

Kinematic dynamo can be also viewed as the phenomenon of the spontaneous breakdown of the topological supersymmetry of the associated stochastic differential equation related to the flow of the background matter.[15] Within supersymmetric theory of stochastics, this supersymmetry is an intrinsic property of all stochastic differential equations, its meaning is the preservation of the continuity of the phase space of the model by continuous time flows, and its spontaneous breakdown is the stochastic generalization of the concept of deterministic chaos.[16] In other words, kinematic dynamo is a manifestation of the chaoticity of the underlying flow of the background matter.

## Nonlinear dynamo theory

The kinematic approximation becomes invalid when the magnetic field becomes strong enough to affect the fluid motions. In that case the velocity field becomes affected by the Lorentz force, and so the induction equation is no longer linear in the magnetic field. In most cases this leads to a quenching of the amplitude of the dynamo. Such dynamos are sometimes also referred to as hydromagnetic dynamos.[17] Virtually all dynamos in astrophysics and geophysics are hydromagnetic dynamos.

The main idea of the theory is that any small magnetic field existing in the outer core, creates currents in the moving fluid there due to Lorenz force. These currents create further magnetic field due to Ampere's law. With the fluid motion, the currents are carried in a way that the magnetic field gets stronger (as long as ${\displaystyle \mathbf {u} \cdot (\mathbf {J} \times \mathbf {B} )}$ is negative[18]). Thus a "seed" magnetic field can get stronger and stronger until it reaches some value that is related to existing non-magnetic forces.

Numerical models are used to simulate fully nonlinear dynamos. The following equations are used:

• The induction equation, presented above.
• Maxwell's equations for negligible electric field:
${\displaystyle \nabla \cdot \mathbf {B} =0}$
${\displaystyle \nabla \times \mathbf {B} =\mu _{0}\mathbf {J} }$
${\displaystyle \nabla \cdot \mathbf {u} =0,}$
• The Navier-Stokes equation for conservation of momentum, again in the same approximation, with the magnetic force and gravitation force as the external forces:
${\displaystyle {\frac {D\mathbf {u} }{Dt}}=-{\frac {1}{\rho _{0}}}\nabla p+\nu \nabla ^{2}\mathbf {u} +\rho '\mathbf {g} +2\mathbf {\Omega } \times \mathbf {u} +\mathbf {\Omega } \times \mathbf {\Omega } \times \mathbf {R} +{\frac {1}{\rho _{0}}}\mathbf {J} \times \mathbf {B} ,}$
where ${\displaystyle \nu }$ is the kinematic viscosity, ${\displaystyle \rho _{0}}$is the mean density and ${\displaystyle \rho '}$ is the relative density perturbation that provides buoyancy (for thermal convection ${\displaystyle \rho '=\alpha \Delta T}$ where ${\displaystyle \alpha }$ is coefficient of thermal expansion), ${\displaystyle \Omega }$ is the rotation rate of the Earth, and ${\displaystyle \mathbf {J} }$ is the electric current density.
• A transport equation, usually of heat (sometimes of light element concentration):
${\displaystyle {\frac {\partial T}{\partial t}}=\kappa \nabla ^{2}T+\epsilon }$
where T is temperature, ${\displaystyle \kappa =k/\rho c_{p}}$ is the thermal diffusivity with k thermal conductivity, ${\displaystyle c_{p}}$ heat capacity, and ${\displaystyle \rho }$ density, and ${\displaystyle \epsilon }$ is an optional heat source. Often the pressure is the dynamic pressure, with the hydrostatic pressure and centripetal potential removed.

These equations are then non-dimensionalized, introducing the non-dimensional parameters,

${\displaystyle Ra={\frac {g\alpha TD^{3}}{\nu \kappa }},E={\frac {\nu }{\Omega D^{2}}},Pr={\frac {\nu }{\kappa }},Pm={\frac {\nu }{\eta }}}$

where Ra is the Rayleigh number, E the Ekman number, Pr and Pm the Prandtl and magnetic Prandtl number. Magnetic field scaling is often in Elsasser number units ${\displaystyle B=(\rho \Omega /\sigma )^{1/2}}$.

### Energy conversion between magnetic and kinematic energy

The scalar product of the above form of Navier-Stokes equation with ${\displaystyle \rho _{0}\mathbf {u} }$ gives the rate of increase of kinetic energy density, ${\displaystyle (1/2)\rho _{0}u^{2}}$, on the left-hand side. The last term on the right-hand side is then ${\displaystyle \mathbf {u} \cdot (\mathbf {J} \times \mathbf {B} )}$, the local contribution to the kinetic energy due to Lorentz force.

The scalar product of the induction equation with ${\displaystyle (1/\mu _{0})\mathbf {B} }$ gives the rate of increase of the magnetic energy density, ${\displaystyle (1/2\mu _{0})B^{2}}$, on the left-hand side. The last term on the right-hand side is then ${\displaystyle (1/\mu _{0})\mathbf {B} \cdot \left(\nabla \times (\mathbf {u} \times \mathbf {B} )\right)}$. Since the equation is volume-integrated, this term is equivalent up to a boundary term (and with the double use of the scalar triple product identity) to ${\displaystyle -\mathbf {u} \cdot \left((1/\mu _{0})(\nabla \times \mathbf {B} )\times \mathbf {B} )\right)=-\mathbf {u} \cdot (\mathbf {J} \times \mathbf {B} )}$ (where one of Maxwell's equations was used). This is the local contribution to the magnetic energy due to fluid motion.

Thus the term ${\displaystyle -\mathbf {u} \cdot (\mathbf {J} \times \mathbf {B} )}$ is the rate of transformation of kinetic energy to magnetic energy. This has to be non-negative at least in part of the volume, for the dynamo to produce magnetic field.[18]

### Order of magnitude of the magnetic field created by Earth's dynamo

The above formula for the rate of conversion of kinetic energy to magnetic energy, is equivalent to a rate of work done by a force of ${\displaystyle \mathbf {J} \times \mathbf {B} }$ on the outer core matter, whose velocity is ${\displaystyle \mathbf {u} }$. This work is the result of non-magnetic forces acting on the fluid.

Of those, the gravitational force and the centrifugal force are conservative and therefore have no overall contribution to fluid moving in closed loops. Ekman number (defined above), which is the ratio between the two remaining forces, namely the viscosity and Coriolis force, is very low inside Earth's outer core, because its viscosity is low (1.2-1.5 x10−2 pascal-second [19]) due to its liquidity.

Thus the main time-averaged contribution to the work is from Coriolis force, whose size is ${\displaystyle -2\rho \,\mathbf {\Omega } \times \mathbf {u} }$, though this quantity and ${\displaystyle \mathbf {J} \times \mathbf {B} }$ are related only indirectly and are not in general equal locally (thus they affect each other but not in the same place and time).

The current density J is itself the result of the magnetic field according to Ohm's law. Again, due to matter motion and current flow, this is not necessarily the field at the same place and time. However these relations can still be used to deduce orders of magnitude of the quantities in question.

In terms of order of magnitude, ${\displaystyle J\,B\sim \rho \,\Omega \,u}$ and ${\displaystyle J\sim \sigma uB}$, giving ${\displaystyle \sigma \,u\,B^{2}\sim \rho \,\Omega \,u}$, or:

${\displaystyle B\sim {\sqrt {\frac {\rho \,\Omega }{\sigma }}}}$

The exact ratio between both sides is the square root of Elsasser number.

Note that the magnetic field direction cannot be inferred from this approximation (at least not its sign) as it appears squared, and is, indeed, sometimes reversed, though in general it lies on a similar axis to that of ${\displaystyle \mathbf {\Omega } }$.

For earth outer core, ρ is approximately 104 kg/m3,[19] Ω=2π/day = 7.3x10−5 seconds and σ is approximately 107Ω−1m−1.[20] This gives 2.7x10−4 Tesla.

The magnetic field of a magnetic dipole has an inverse cubic dependence in distance, so its order of magnitude at the earth surface can be approximated by multiplying the above result with (Router core/REarth)3 = (2890/6370)3 = 0.093, giving 2.5x10−5 Tesla, not far from to the measured value of 3x10−5 Tesla at the equator.

## Numerical models

A visual representation of the Glatzmaier model before dipole reversal

Broadly, models of the geodynamo attempt to produce magnetic fields consistent with observed data given certain conditions and equations as mentioned in the sections above. Implementing the magnetohydrodynamic equations successfully was of particular significance because they pushed dynamo models to self-consistency. Though geodynamo models are especially prevalent, dynamo models are not necessarily restricted to the geodynamo; solar and general dynamo models are also of interest. Studying dynamo models has utility in the field of geophysics as doing so can identify how various mechanisms form magnetic fields like those produced by astrophysical bodies like Earth and how they cause magnetic fields to exhibit certain features, such as pole reversals.

The equations used in numerical models of dynamo are highly complex. For decades, theorists were confined to two dimensional kinematic dynamo models described above, in which the fluid motion is chosen in advance and the effect on the magnetic field calculated. The progression from linear to nonlinear, three dimensional models of dynamo was largely hindered by the search for solutions to magnetohydrodynamic equations, which eliminate the need for many of the assumptions made in kinematic models and allow self-consistency.

A visual representation of the Glatzmaier model during dipole reversal

The first self-consistent dynamo models, ones that determine both the fluid motions and the magnetic field, were developed by two groups in 1995, one in Japan[21] and one in the United States.[22][23] The latter was made as a model with regards to the geodynamo and received significant attention because it successfully reproduced some of the characteristics of the Earth's field.[18] Following this breakthrough, there was a large swell in development of reasonable, three dimensional dynamo models.[18]

Though many self-consistent models now exist, there are significant differences among the models, both in the results they produce and the way they were developed.[18] Given the complexity of developing a geodynamo model, there are many places where discrepancies can occur such as when making assumptions involving the mechanisms that provide energy for the dynamo, when choosing values for parameters used in equations, or when normalizing equations. In spite of the many differences that may occur, most models have shared features like clear axial dipoles. In many of these models, phenomena like secular variation and geomagnetic polarity reversals have also been successfully recreated.[18]

#### Observations

A visual representation of the Glatzmaier model after dipole reversal

Many observations can be made from dynamo models. Models can be used to estimate how magnetic fields vary with time and can be compared to observed paleomagnetic data to find similarities between the model and the Earth. Due to the uncertainty of paleomagnetic observations, however, comparisons may not be entirely valid or useful[18]. Simplified geodynamo models have shown relationships between the dynamo number (determined by variance in rotational rates in the outer core and mirror-asymmetric convection (e.g. when convection favors one direction in the north and the other in the south)) and magnetic pole reversals as well as found similarities between the geodynamo and the Sun's dynamo.[18] In many models, it appears that magnetic fields have somewhat random magnitudes that follow a normal trend that average to zero.[18] In addition to these observations, general observations about the mechanisms powering the geodynamo can be made based on how accurately the model reflects actual data collected from Earth.

#### Modern Modelling

The complexity of dynamo modelling is so great that models of the geodynamo are limited by the current power of supercomputers, particularly because calculating the Ekman and Rayleigh number of the outer core is extremely difficult and requires a vast number of computations.

Many improvements have been proposed in dynamo modelling since the self-consistent breakthrough in 1995. One suggestion in studying the complex magnetic field changes is applying spectral methods to simplify computations.[24] Ultimately, until considerable improvements in computer power are made, the methods for computing realistic dynamo models will have to be made more efficient, so making improvements in methods for computing the model is of high importance for the advancement of numerical dynamo modelling.

## References

1. ^ "How does the Earth's core generate a magnetic field?". USGS FAQs. United States Geological Survey. Archived from the original on 18 January 2015. Retrieved 21 October 2013.
2. ^ Larmor, J. (1919). "How could a rotating body such as the Sun become a magnet?". Reports of the British Association. 87: 159–160.
3. ^ Larmor, J. (1919). "Possible rotational origin of magnetic fields of sun and earth". Electrical Review. 85: 412ff. Reprinted in Engineering, vol. 108, pages 461ff (3 October 1919).
4. ^ Nye, Mary Jo (1 March 1999). "Temptations of theory, strategies of evidence: P. M. S. Blackett and the earth's magnetism, 1947–52". The British Journal for the History of Science. 32 (1): 69–92. doi:10.1017/S0007087498003495.
5. ^ Merrill, McElhinny & McFadden 1996, page 17 claim that in 1905, shortly after composing his special relativity paper, Albert Einstein described the origin of the Earth's magnetic field as being one of the great unsolved problems facing modern physicists. However, they do not provide details on where he made this statement.
6. ^ Sanders, Robert (2003-12-10). "Radioactive potassium may be major heat source in Earth's core". UC Berkeley News. Retrieved 2007-02-28.
7. ^ Sakuraba, Ataru; Paul H. Roberts (4 October 2009). "Generation of a strong magnetic field using uniform heat flux at the surface of the core". Nature Geoscience. 2 (11): 802–805. Bibcode:2009NatGe...2..802S. doi:10.1038/ngeo643.
8. ^ Buffett, Bruce (2009). "Geodynamo: A matter of boundaries". Nature Geoscience. 2 (2): 741–742. Bibcode:2009NatGe...2..741B. doi:10.1038/ngeo673.
9. ^ Brandenburg, Axel (2007). "Hydromagnetic dynamo theory". Scholarpedia. Retrieved October 7, 2016.
10. ^ E. Pallé (2010). The Earth as a Distant Planet: A Rosetta Stone for the Search of Earth-Like Worlds (Astronomy and Astrophysics Library). Berlin: Springer. pp. 316–317. ISBN 978-1-4419-1683-9. Retrieved 17 July 2010.
11. ^ Steigerwald, Bill (October 6, 2010). "Saturn's Icy Moon May Keep Oceans Liquid with Wobble". NASA. Retrieved August 14, 2012.
12. ^ Cassis, Nikki (March 19, 2012). "Geologic map of Jupiter's moon Io details an otherworldly volcanic surface". Astrogeology Science Center. Retrieved August 14, 2012.
13. ^ "Mercury's Surprising Core and Landscape Curiosities". MESSENGER. Carnegie Institution for Science. March 21, 2012. Retrieved August 14, 2012.
14. ^ Stevens, Tim (November 9, 2011). "Ancient lunar dynamo may explain magnetized moon rocks". University of California. Retrieved August 14, 2012.
15. ^ Ovchinnikov, I.V.; Ensslin, T. A. (April 2016). "Kinematic dynamo, supersymmetry breaking, and chaos". Physical Review D. 93 (8): 085023. arXiv:1512.01651. Bibcode:2016PhRvD..93h5023O. doi:10.1103/PhysRevD.93.085023.
16. ^ Ovchinnikov, I.V. (March 2016). "Introduction to Supersymmetric Theory of Stochastics". Entropy. 18 (4): 108. arXiv:1511.03393. Bibcode:2016Entrp..18..108O. doi:10.3390/e18040108.
17. ^ Parker, Eugene N. (September 1955). "Hydromagnetic Dynamo Models". The Astrophysical Journal. 122: 293–314. Bibcode:1955ApJ...122..293P. doi:10.1086/146087.
18. Kono, Masaru; Roberts, Paul H. (2002). "Recent geodynamo simulations and observations of the geomagnetic field". Reviews of Geophysics. 40 (4): 1–53. Bibcode:2002RvGeo..40.1013K. doi:10.1029/2000RG000102.
19. ^ a b de Wijs, G. A., Kresse, G., Vočadlo, L., Dobson, D., Alfe, D., Gillan, M. J., & Price, G. D. (1998). The viscosity of liquid iron at the physical conditions of the Earth's core. Nature, 392(6678), 805.
20. ^ Ohta, K., Kuwayama, Y., Hirose, K., Shimizu, K., & Ohishi, Y. (2016). Experimental determination of the electrical resistivity of iron at Earth’s core conditions. Nature, 534(7605), 95. Link to a summary
21. ^ Kageyama, Akira; Sato, Tetsuya (1 January 1995). "Computer simulation of a magnetohydrodynamic dynamo. II". Physics of Plasmas. 2 (5): 1421–1431. Bibcode:1995PhPl....2.1421K. doi:10.1063/1.871485.
22. ^ Glatzmaier, Gary A.; Roberts, Paul H. (1995). "A three-dimensional self-consistent computer simulation of a geomagnetic field reversal". Nature. 377 (6546): 203–209. Bibcode:1995Natur.377..203G. doi:10.1038/377203a0.
23. ^ Glatzmaier, G; Roberts, Paul H. (1995). "A three-dimensional convective dynamo solution with rotating and finitely conducting inner core and mantle". Physics of the Earth and Planetary Interiors. 91 (1–3): 63–75. Bibcode:1995PEPI...91...63G. doi:10.1016/0031-9201(95)03049-3.
24. ^ Avery, Margaret S.; Constable, Catherine G.; Davies, Christopher J.; Gubbins, David (2019-01-01). "Spectral methods for analyzing energy balances in geodynamo simulations". Physics of the Earth and Planetary Interiors. 286: 127–137. doi:10.1016/j.pepi.2018.10.002. ISSN 0031-9201.