3 Induced and remanent magnetism

Chapter 3
Induced and remanent magnetism

BACKGROUND: For a review of basic quantum mechanics and statistical mechanics, read relevant chapters from an introductory Chemistry text book.

Scientists in the late 19th century thought that it might be possible to exploit the magnetic record retained in accidental records to study the geomagnetic field in the past. Work in the mid 20th century provided the theoretical and experimental basis for presuming that such materials might retain a record of past geomagnetic fields. There are several books and articles that describe the subject in detail (see e.g., the supplemental readings). We present here a brief overview of theories on how rocks get and stay magnetized. We will begin with magnetism at the atomic level caused by electronic orbits and spins giving rise to induced magnetizations. Then we will see how electronic spins working in concert give rise to permanently magnetized substances (like magnetic minerals) making remanent magnetization possible.

3.1 Magnetism at the atomic level

We learned in Chapter 1 that magnetic fields are generated by electric currents. Given that there are no wires leading into or out of permanent magnets, you may well ask, “Where are the currents?” At the atomic level, the electric currents come from the motions of the electrons. From here quantum mechanics quickly gets esoteric, but some rudimentary understanding is helpful. In this chapter we will cover the bare minimum necessary to grasp the essentials of rock magnetism.

In Chapter 1 we took the classical (pre-quantum mechanics) approach and suggested that the orbit of an electron about the nucleus could be considered a tiny electric current with a correspondingly tiny magnetic moment. But quantum physics tells us that this “planetary” view of the atom cannot be true. An electron zipping around a nucleus would generate radio waves, losing energy and eventually would crash into the nucleus.


Figure 3.1: Plot of radial distribution and “dot-density” for the 1s electron shell.

Apparently, this does not happen, so the classical approach is fatally flawed and we must turn to quantum mechanics.

In quantum mechanics, electronic motion is stabilized by the fact that electrons can only have certain energy states; they are quantized. The energy of a given electron can be described in terms of solutions, Ψ, to something called Schr�dinger’s wave equation. The function Ψ(r,θ,ϕ) gives the probability of finding an electron at a given position. [Remember from Chapter 2 that r,θ,ϕ are the three spherical coordinates.] It depend on three special quantum numbers (n,l,m):
          l    m
Ψr,θ,ϕ = R n(r)Yl (ϕ,θ),

The number n is the so-called “principal” quantum number. The Rnl(r) are functions specific to the element in question and the energy state of the electron n. It is evaluated at an effective radius r in atomic units. The Y lm are a fully normalized complex representation of the spherical harmonics introduced in Section 2.2. For each level n, the number l ranges from 0 to n-1 and m from l backwards to -l.

The lowest energy of the quantum wave equations is found by setting n equal to unity and both l and m to zero. Under these conditions, the solution to the wave equation is given by:

          3 -ρ∕2
R1,0 = 2Z 2e   ,
        1  1
Y0,0 = (--)2

where Z is the atomic number and ρ is 2Zr∕n. Note that at this energy level, there is no dependence of Y on ϕ or θ. Substituting these two equations into Equation 3.1 gives the probability density Ψ for an electron as a function of radius of r. This is sketched as the line in Figure 3.1. Another representation of the same idea is shown in the inset, whereby the density of dots at a given radius reflects the probability distribution shown by the solid curve. The highest dot density is found at a radius of about one atomic unit, tapering off the farther away from the center of the atom. Because there is no dependence on θ or ϕ the probability distribution is a spherical shell. All the l,m = 0 shells are spherical and are often referred to as the 1s, 2s, 3s shells, where the numbers are the energy levels n. A surface with equal probability is a sphere and example of one such shell is shown in Figure 3.2a.

For l = 1, m will have values of -1, 0 and 1 and the Y lm(ϕ,θ)s are given by:

        ∘ ---                  ∘ --               ∘ ---
Y-11 = 1-  -3-sin θe-iϕ,  Y01 = 1-  3-cosθ,  Y11= --1   3--sinθeiϕ.
      2   2π                 2   π              2    2π
Shells with l = 1 depend not only on radial distance but also on the angles ϕ and θ, so they are not spheres, but more complicated shapes. A surface of equal probability for one such shell (the m = 1 shell) is shown in Figure 3.2b. Shells with l = 1 are called the “p” shells.

As might be expected, the shells for l = 2 are even more complicated that for l = 1. These shells are called “d” shells and two examples are shown in Figure 3.2c and d.


Figure 3.2: Examples of surfaces of equal probability of the first three shells (l = 1,2,3). Surfaces created with Orbital Viewer.

Returning to the tiny circuit idea, somehow the motion of the electrons in their shells acts like an electronic circuit and creates a magnetic moment. In quantum mechanics, the angular momentum vector of the electron L is quantized, for example as integer multiples of , the “reduced” Planck’s constant (or h _ 2π where h = 6.63 x 10-34 Js). The magnetic moment arising from the orbital angular momentum is given by:

|m | = --e-|L|,
where μe is the mass of an electron (9.11 x 10-31 kg), qe= -1.69 x 10-19C. The smallest value of L is so the fundamental unit of magnetic moment arising from the oribit of electrons is given by:
       -ℏqe           -24kg-m2-  C--          -24   2
|mb | = 2μe = 9.27 × 10     s   ⋅ kg = 9.27 × 10   Am  .

This is known as the Bohr magneton.


Figure 3.3: Electronic structure of elements from Na to Zn.

So far we have not mentioned one last quantum number, s. This is the “spin” of the electron and has a value of 1 2. The spin itself produces a magnetic moment which is given by 2smb, hence is numerically identical to that produced by the orbit.

Atoms have the same number of electrons as protons in order to preserve charge balance. Hydrogen has but one lonely electron which in its lowest energy state sits in the 1s electronic shell. Helium has a happy pair, so where does the second electron go? To fill in their electronic shells, atoms follow three rules:

  1. No two electrons may have the same set of quantum numbers. This is Pauli’s exclusion principle. Because spin (s) can be 1 2, two electrons fit in one orbital. When a single electron occupies a given orbital, it is called “unpaired” and has a magnetic moment of 1 mb.
  2. Orbitals are filled in order of increasing energy. The energy state of a given orbital is dependent on the context (whether the atom is bound to other atoms or not), but in general they will be filled according to the scheme shown in Figure 3.3.
  3. Electrons are added so that the spins remain as parallel as possible (Hund’s rule). Notice in Figure 3.3 that when filling the third energy level (n = 3), all five d shells are filled with one kind of spin (say, all up, or +1 2), before the electrons begin to pair up. Also, because the energies of the shells change somewhat according to the context they are in, the 4s shell will actually give up an electron to a d shell, before the d shells begin to pair up. Hund’s rule gives the atoms with some d shell electrons (the so-called “transition elements”, e.g., Cr, Mn, Fe, Co and Ni) the possibility of large magnetic moments.

Each unpaired spin has a moment of one Bohr magneton mb. The elements with the most unpaired spins are the transition elements which are responsible for most of the paramagnetic behavior observed in rocks. For example, in Figure 3.3 we see that Mn has a structure of: (1s22s22p63s23p6)3d54s2, hence has five unpaired spins and a net moment of 5 mb. Fe has a structure of (1s22s22p63s23p6)3d64s2 with a net moment of 4 mb, In minerals, the transition elements are in a variety of oxidation states. Fe commonly occurs as Fe2+ and Fe3+. When losing electrons to form ions, transition metals lose the 4s electrons first, so we have for example, Fe3+ with a structure of (1s22s22p63s23p6)3d5, or 5 mb. Similarly Fe2+ has 4 mb and Ti4+ has no unpaired spins. Iron is the main magnetic species in geological materials, but Mn2+ (5 mb) and Cr3+ (3 mb) occur in trace amounts.

3.2 Induced magnetization

We have learned that there are two sources of magnetic moments in electronic motions: the orbits and the (unpaired) spins. These moments respond to external magnetic fields giving rise to an induced magnetization, a phenomenon alluded to briefly in Chapter 1. We will consider first the contribution of the electronic orbits.


Figure 3.4: Larmor precession. The orbit of the electron has an angular momentum vector L which creates a magnetic moment. In the presence of a magnetic field H, the moment experiences a torque which causes a change in angular momentum ΔL. The precession of the electronic orbit about H creates an induced magnetic moment Δm in a sense opposite to the applied field H.

3.2.1 Orbital contribution and diamagnetism

The angular momentum of electrons is quantized in magnitude but also has direction (see L in Figure 3.4). The angular momentum vector has an associated magnetic moment vector mb. A magnetic field H exerts a torque on the moment, which nudges it (and the momentum vector associated with it) to the side (ΔL). L therefore will precess around the magnetic field direction, much like a spinning top precesses around the direction of gravity. The precession of L is called Larmor precession.

The changed momentum vector from Larmor precession in turn results in a changed magnetic moment vector Δm. The sense of the change in net moment is always to oppose the applied field. Therefore, the response of the magnetic moments of electronic orbitals creates an induced magnetization MI that is observable outside the substance; it is related to the applied field by:

MI  = χdH.

We learned in Chapter 1 that the proportionality between induced magnetization and the applied field is known as the magnetic susceptibility. The ratio MIH for the response of the electronic orbitals is termed the diamagnetic susceptibility χd; it is negative, essentially temperature independent and quite small. This diamagnetic response is a property of all matter, but for substances whose atoms possess atomic magnetic moments, diamagnetism is swamped by effects of magnetic fields on the atomic magnetic moments. In the absence of unpaired electronic spins, diamagnetic susceptibility dominates the magnetic response. Common diamagnetic substances include quartz (SiO2), calcite (CaCO3) and water (H2O). The mass normalized susceptibility of quartz is -0.62 x 10-9 m3kg-1 to give you an idea of the magnitudes of these things.


Figure 3.5: a) Paramagnetic magnetization (obtained from the Langevin function L(a) versus a = mB∕kT.) b) Paramagnetic magnetization as a function of temperature (Curie Law).

3.2.2 Role of electronic spins and paramagnetism

In many geological materials, the orbital contributions cancel out because they are randomly oriented with respect to one another and the magnetization arises from the electronic spins. We mentioned that unpaired electronic spins behave as magnetic dipoles with a moment of one Bohr magneton. In the absence of an applied field, or in the absence of the ordering influence of neighboring spins which are known as exchange interactions, the electronic spins are essentially randomly oriented. An applied field acts to align the spins which creates a net magnetization equal to χpH where χp is the paramagnetic susceptibility. For any geologically relevant conditions, the induced magnetization is linearly dependent on the applied field. In paramagnetic solids, atomic magnetic moments react independently to applied magnetic fields and to thermal energy. At any temperature above absolute zero, thermal energy vibrates the crystal lattice, causing atomic magnetic moments to oscillate rapidly in random in orientations. In the absence of an applied magnetic field, atomic moments are equally distributed in all directions with a resultant magnetization of zero.

A useful first order model for paramagnetism was worked out by P. Langevin in 1905. (Of course in messy reality things are a bit more complicated, but Langevin theory will work well enough for us at this stage.) Langevin theory is based on a few simple premises:

  1. Each unpaired spin contributes a dipole moment.
  2. In the absence of an applied field, the moments are essentially randomly oriented, i.e., all directions are equally likely to occur.
  3. Application of a magnetic field exerts an aligning torque on the atomic magnetic moments. The magnetic energy Em (see also Section 1.4 in Chapter 1) of a magnetic moment m at an angle θ with an external magnetic field B = μoH is given by:
    Em  = - m ⋅B  = - mB cos θ.

    Magnetic energy is at a minimum when the magnetic moment is lined up with the magnetic field.

  4. There is competition between the magnetic energy Em and the thermal energy kT where k is Boltzmann’s constant (1.38 x 10-23 m2 kg s-2 K-1) and T is temperature in kelvin).

Consider an atomic magnetic moment, (m = 2mb = 1.85×10-23 Am2), in a magnetic field of 10-2 T, (for reference, the largest geomagnetic field at the surface is about 65 μT – see Chapter 2). The aligning energy is therefore mB = 1.85 × 10-25 J). However, thermal energy at 300K (traditionally chosen as a temperature close to room temperature providing easy arithmetic) is Boltzmann’s constant times the temperature, or about 4 x 10-21 J. So thermal energy is several orders of magnitude larger than the aligning energy and the net magnetization is small even in this rather large (compared to the Earth’s field) magnetizing field.

Using the principles of statistical mechanics, we find that the probability density of a particular magnetic moment having a magnetic energy of Em is given by:
P (E ) ∝ exp (- Em ∕kT ).

From this we see that the degree of alignment depends exponentially on the ratio of magnetic energy to thermal energy. The degree of alignment with the magnetic field controls the net magnetization M. When spins are completely aligned, the substance has a saturation magnetization Ms. The probability density function leads directly to the following relation (derived in Appendix A.2.1):
M              1
--- = [coth a - -] = L(a).
Ms             a

where a = mB∕kT. The function enclosed in square brackets is known as the Langevin function (L).

Equation 3.6 is plotted in Figure 3.5a and predicts several intuitive results: 1) M = 0 when B = 0 and 2) M∕Ms = 1 when the applied magnetic field is infinite. Furthermore, M is some 90% of Ms when mB is some 10-20 times kT. When kT >> mB,L(a) is approximately linear with a slope of ~ 13. At room temperature and fields up to many tesla, L(a) is approximately mB∕3kT. If the moments are unpaired spins (m = mb), then the maximum magnetization possible (Ms) is given by the number of moments N, their magnitude (mb) normalized by the volume of the material v or Ms = Nmb∕v, and

M     mbμo
M--≃  3kT--H.

Please note that we have neglected all deviations from isotropy including quantum mechanical effects as well as crystal shape, lattice defects, and state of stress. These complicate things a little, but to first order the treatment followed here provides a good approximation. We can rewrite the above equation as:
M     m μ         N m2 μ   1
---=  --b-o⋅Ms  = ----b-o ⋅--= χp.
H     3kT           3kv    T

To first order, paramagnetic susceptibility χp is positive, larger than diamagnetism and inversely proportional to temperature. This inverse T dependence (see Figure 3.5b) is known as Curie’s law of paramagnetism. The paramagnetic susceptibility of, for example, biotite is 790 x 10-9 m3 kg-1, or about three orders of magnitude larger than quartz (and of the opposite sign!).

We have considered the simplest case here in which χ can be treated as a scalar and is referred to as the bulk magnetic susceptibility χb. In detail, magnetic susceptibility can be quite complicated. The relationship between induced magnetization and applied field can be affected by crystal shape, lattice structure, dislocation density, state of stress, etc., which give rise to possible anisotropy of the susceptibility. Furthermore, there are only a finite number of electronic moments within a given volume. When these are fully aligned, the magnetization reaches saturation. Thus, magnetic susceptibility is both anisotropic and non-linear with applied field.


Figure 3.6: Exchange energy associated with overlapping orbitals. Example of super-exchange between the 3d orbitals of two iron cations through the 2p orbitals of the intervening oxygen anion. The two electrons in the 2p shells are, by necessity antiparallel. These are shared by the 3d shells, hence the two cations have anti-parallel spins. [Figure redrawn from O’Reilly, 1984.]

3.3 Ferromagnetism

Some substances give rise to a magnetic field in the absence of an applied field. This magnetization is called remanent or spontaneous magnetization, also loosely known as ferromagnetism (sensu lato). Magnetic remanence is caused by strong interactions between neighboring spins that occur in certain crystals.

The so-called exchange energy is minimized when the spins are aligned parallel or anti-parallel depending on the details of the crystal structure. Exchange energy is a consequence of the Pauli exclusion principle (no two electrons can have the same set of quantum numbers). In the transition elements, the 3d orbital is particularly susceptible to exchange interactions because of its shape and the prevalence of unpaired spins, so remanence is characteristic of certain crystals containing transition elements with unfilled 3d orbitals.

In oxides, oxygen can form a bridge between neighboring cations which are otherwise too far apart for direct overlap of the 3d orbitals in a phenomenon known as superexchange. In Figure 3.6 the 2p electrons of the oxygen are shared with the neighboring 3d shells of the iron ions. Pauli’s exclusion principle means that the shared electrons must be antiparallel to each of the electrons in the 3d shells. The result is that the two cations are coupled. In the case shown in Figure 3.6 there is an Fe2+ ion coupled antiparallel to an Fe3+ ion. For two ions with the same charge, the coupling will be parallel. Exchange energies are huge, equivalent to the energy associated with the same moment in a field of the order of 1000 T. [The largest field available in the Scripps paleomagnetic laboratory is about 2.5 T, and that only fleetingly.]

As temperature increases, crystals expand and exchange becomes weaker. Above a temperature characteristic of each crystal type (known as the Curie temperature Tc), cooperative spin behavior disappears entirely and the material becomes paramagnetic.


Figure 3.7: Behavior of magnetization versus temperature of a ferromagnetic substance. Below Tc, the magnetization follows Equation 3.9 and is the ferromagnetic magnetization. Above Tc the magnetization follows Equation 3.8 and is the induced magnetization. [Redrawn from Tauxe, 1998.]

While the phenomenon of ferromagnetism results from complicated interactions of neighboring spins, it is useful to think of the ferromagnetic moment as resulting from a quasi-paramagnetic response to a huge internal field. This imaginary field is termed the Weiss molecular field Hw. In Weiss theory, Hw is proportional to the magnetization of the substance, i.e.,

Hw  = βM,

where β is the constant of proportionality. The total magnetic field that the substance experiences is:

Htot = H + Hw  = H  + βM,

where H is the external field. By analogy to paramagnetism, we can substitute a = μomb(Htot)∕kT) for H in Langevin function:
 M      (μ m  (H + βM  ))
--- = L  --o-b---------- .
Ms             kT

For temperatures above the Curie temperature Tc (i.e. T - Tc > 0) there is by definition no internal field, hence βM is zero. Substituting Nmb∕v for Ms, and using the low-field approximation for L(a), Equation 3.8 can be rearranged to get:
M       μoN m2
---=  -------b---≡  χf.
H     v3k(T - Tc)

Equation 3.9 is known as the Curie-Weiss law and governs ferromagnetic susceptibility above the Curie temperature (dashed line in Figure 3.7).


Figure 3.8: Various data sets for the behavior of Ms(T) for magnetite.

Below the Curie temperature Hw >> H; we can neglect the external field H and get:

M--     μomb-βM--
M   = L(   kT    ).

Substituting again for Ms and rearranging, we get:
M-- = L(N-m-bβ ⋅-M-) = L( Tc⋅ M--),
Ms        vkT   Ms        T   Ms

where Tc is the Curie temperature and is given by:

     N m2 β
Tc = ----b-.
Equation 3.10 can be solved graphically or numerically and is sketched (solid line) in Figure 3.7. Below the Curie temperature, exchange interactions are strong relative to the external field and the magnetization is governed by Equation 3.10. Above the Curie temperature, it follows the Curie-Weiss law (Equation 3.9).

We have treated ferromagnetism from a classical point of view and this is strictly incorrect because ferromagnetism results primarily from quantum mechanical phenomena. The primary difference between the classical derivation and the quantum mechanical one lies in the fact that in quantum mechanics, only certain angles of the magnetic moments are allowed, as opposed to all directions in Langevin theory. In the end, the predictions of magnetization as a function of temperature are different in detail. The end product of the quantum mechanical treatment (see Dunlop and �zdemir, 1997) is that the variation of saturation magnetization as a function of temperature can be reasonably well approximated (near the Curie Temperature, Tc) by a normalized power law variation:
 Ms(T )    Tc - T γ
M--(T-) = [T---T--],
  s  o      c    o

where γ is 0.5 from simple molecular field theory and To is absolute zero (in kelvin). Dunlop and �zdemir (1997) cite a value of around 0.43 for γ, but the data sets cited by Dunlop and �zdemir (1997; e.g., Figure 3.5 on page 52) are actually best-fit with values for γ of about 0.36 – 0.39 (see Figure 3.8). These curves have been normalized by their inferred curie Temperatures which are around 565C (data of B. Moskowitz, cited in Banerjee, 1991).


Figure 3.9: Types of spin alignment in ferromagnetism (sensu lato): a) ferromagnetism (sensu stricto), b) antiferromagnetism, c) spin-canted antiferromagnetism, d) defect anti-ferromagnetism, e) ferrimagnetism.

As we have seen, below the Curie temperature, certain crystals have a permanent (remanent) magnetization resulting from the alignment of unpaired electronic spins over a large area within the crystal. Spins may be either parallel or anti-parallel; the sense of spin alignment is controlled entirely by crystal structure. The energy term associated with this phenomenon is the exchange energy. There are three categories of spin alignment: ferromagnetism (sensu stricto), ferrimagnetism and antiferromagnetism (see Figure 3.9).


Figure 3.10: a) Response of a magnetic moment to the torque of an applied field for isolated moments. b) Response of coupled moments to a perturbation. Neighboring spins produce an effect known as “spin waves”.

In ferromagnetism (sensu stricto, Figure 3.9a), the exchange energy is minimized when all the spins are parallel, as occurs in pure iron. When spins are perfectly antiparallel (antiferromagnetism, Figure 3.9b), there is no net magnetic moment, as occurs in ilmenite. Occasionally, the antiferromagnetic spins are not perfectly aligned in an antiparallel orientation, but are canted by a few degrees. This spin-canting (Figure 3.9c) gives rise to a weak net moment, as occurs in hematite, a common magnetic mineral (see Chapter 6). Also, antiferromagnetic materials can have a net moment if spins are not perfectly compensated owing to defects in the crystal structure, as occurs in fine-grained hematite. The uncompensated spins result in a so-called defect moment (Figure 3.9d). We note in passing that the temperature at which spins become disordered in antiferromagnetic substances is termed the N�el temperature. In ferrimagnetism, spins are also aligned antiparallel, but the magnitudes of the moments in each direction are unequal, resulting in a net moment (Figure 3.9e).

In figures like Figure 3.9, electronic spins are depicted as being simply aligned with some minimum energy direction (aligned with the field, or along some easy axis). Yet we already know about the paramagnetic effect of misalignment through random thermal fluctuations. We learned that an external magnetic field generates a torque on the electronic spins, and in isolation, a magnetic moment will respond to the torque in a manner similar in some respects to the way a spinning top responds to gravity: the magnetic moment will precess about the applied field direction, spiraling in and come to a rest parallel to it (Figure 3.10a). Because of the strong exchange coupling in ferromagnetic phases, spins tend to be aligned parallel (or antiparallel) to one another and the spiralling is done in a coordinated fashion, with neighboring spins as parallel as possible to one another (Figure 3.10b). This phenomenon is known as a spin wave.

SUPPLEMENTAL READINGS: O’Reilly (1984), Chapter 3.1; Dunlop and �zdemir (1997), Chapter 2.1 to 2.7.

3.4 Problems

Problem 1

a) Given one Bohr magneton (mb) in the Earth’s field (40 μT), write a program using Python that calcuates magnetostatic interaction energy (-mbB cosθ) for angles 0180. Make a plot of this with the matplotlib module in Python.

b) Calculate the thermal energy at room temperature (300K). How does this compare with the interaction energy?

Problem 2

Fayalite (Fe2SiO4) is a paramagnetic solid with magnetic susceptibility χ = 4.4 x 10-4 (cgs units) at 0C (= 273K). A single crystal of fayalite has a volume of 2 cm3. This crystal is placed in a magnetic field, H = 10 oe at 0C. What is the resulting induced magnetic moment m of this crystal?

a) Do this problem first in cgs units. Then convert your answer to SI using the conversion factors in Table 1.1 in Chapter 1.

b) Do the problem again by first converting all the parameters into SI units. Check your answer by converting the SI answer that you get back to cgs. You should get the same answer (but you would be surprised how many people do this wrong).

Problem 3

If fayalite is placed in a magnetic field H= 100 oe at a temperature of 500C (= 773K), what is the resulting magnetization, M?

Problem 4

MnS is a paramagnetic solid. At 300K there are 4 x 1028 molecules of MnS per m3. Look up the number of unpaired spins for the cationic magnetic moment of Mn2+ in the text and find the paramagnetic susceptibility, χ, of MnS at 300K?

Problem 5

a) Read into a Pandas DataFrame the datafile Chapter_3/BMoskinBan91.txt provided. Make a plot of magnetization versus temperature. What is the Curie temperature of the material?

b) Using this Equation 3.11 from the chapter, find the value for γ between 0.35 and 0.43 at intervals of 0.01 that fits the best. Plot the data as in Figure 3.8 in the chapter, i.e. Ms(T)∕Ms(To) against T∕Tc.