Suggested Supplemental Reading
In this lecture we will review the basic physical principles behind magnetism. We will be using primarily the Système International (SI) units that are based on meters-kilograms-seconds. There are other systems of units that are important in magnetism and the most prevalent of these (electromagnetic units of the cgs system) will be covered later in the lecture.
Magnetic fields, like gravitational fields, cannot be seen or touched. We can feel the pull of the Earth’s gravitational field on ourselves and the objects around us, but we do not experience magnetic fields in such a direct way. We know of the existence of magnetic fields by their effect on objects such as magnetized pieces of metal, naturally magnetic rocks such as lodestone, or temporary magnets such as copper coils that carry an electrical current. If we place a magnetized needle on a cork in a bucket of water, it will slowly align itself with the local magnetic field. Turning on the current in a copper wire can make a nearby compass needle jump. Observations like these led to the development of the concept of magnetic fields.
Electric currents make magnetic fields, so we can define what is meant by a “magnetic field” in terms of the electric current that generates it. Figure 1.1a is a picture of what happens when we pierce a flat sheet with a wire carrying a current i. When iron filings are sprinkled on the sheet, the filings line up with the magnetic field produced by the current in the wire. A circle tangential to the field is shown in Figure 1.1b, which illustrates the right-handrule (see inset to Figure 1.1b). If your right thumb points in the direction of (positive) current flow (the direction opposite the flow of the electrons), your fingers will curl in the direction of the magnetic field.
The magnetic field H points at right angles to both the direction of current flow and to the radial unit vector r in Figure 1.1b. The magnitude of H is proportional to the strength of the current i. In the simple case illustrated in Figure 1.1b the magnitude of H is given by Ampère’s law:
So, now we know the units of H: Am-1.
Ampère’s Law in its most general form is one of Maxwell’s equations of electromagnetism: In a steady electrical field,
× H = Jf, where Jf is the electric current density. In other words: The curl (or circulation) of the magnetic field is equal to the current density. The origin of the term “curl” for the cross product of the gradient operator with a vector field is suggested in Figure 1.1a in which the iron filings seem to curl around the wire.
We have seen that an electrical current in a wire produces a magnetic field that curls around the wire. If we bend the wire into a loop with an area
r2 that carries a current i, as shown in Figure 1.2a, the current loop creates the magnetic field shown by pattern of the the iron filings. This magnetic field is that same as the field that would be produced by a magnet with a magnetic moment m shown in Figure 1.2b. This moment is created by the current i and also depends on
the area of the current loop (the bigger the loop, the bigger the moment), hence m = i
r2. The moment created by a set of loops (as shown in Figure 1.2c is the sum of the n individual loops, i.e.:
| (1.1) |
So, now we know the units of m: Am2.
|
|
The magnetic field is a vector field because at any point the field has both direction and magnitude. Consider the field of a bar magnet made visible by iron filings as shown in Figure 1.3. The direction of the field at any point is given by the arrows while the strength depends on how close the field lines are to one another. The magnetic field lines are known as “magnetic flux”. The density of flux lines is one measure of the strength of the magnetic field: the magnetic induction B.
|
|
Magnetic flux density (i.e., magnetic induction) can therefore be quantified when a conductor moves through it. Magnetic induction can be thought of as something that creates a potential difference with voltage V in a conductor of length l when the conductor moves relative to the magnetic induction B with velocity v (see Figure 1.3): V = vlB. From this we can derive the units of magnetic induction: the tesla (T). One tesla is the magnetic induction that generates a potential of one volt in a conductor of length 1 meter when moving 1 meter per second. Ergo, 1 T = 1 V . s . m-2.
Another way of looking at this is that if magnetic induction (B) is the flux density, this must be the flux
per unit area. So an increment of flux d
is the field B times the increment of area dA. The area here is the length of the wire l times its displacement ds in time dt. The instantaneous velocity is dv = ds/dt so or d
= BdA and the rate of change of flux is:
| (1.2) |
Equation 1.2 is known as Faraday’s law and in its most general form is the fourth of Maxwell’s equations. We see from this equation that the units of magnetic flux must be a volt-second which is unit in its own right, the weber (Wb). The weber is defined as the amount of magnetic flux which, when passed through a one-turn coil of a conductor carrying a current of one ampere, produces an electric potential of one volt. This definition suggests a means to measure the strength of magnetic induction and is the basis of the “fluxgate” magnetometer.
A magnetic moment m has a magnetostatic energy (Em) associated with it. This is the energy that tends to align compass needles with the magnetic field (see Figure 1.4. This energy is given by m . B or mB cos
where m and B are the magnitudes of
m and B, respectively. Magnetic energy has of joules.
|
|
Magnetization M is a moment per unit volume (units of Am-1) or per unit mass (Am2kg-1). Sub-atomic charges such as protons and electrons can be thought of as tracing out tiny circuits and behaving as tiny magnetic moments. They respond to external magnetic fields and give rise to an induced magnetization. The relationship between the magnetization induced in a material MI and the external field H is defined as:
| (1.3) |
The parameter
b is known as the bulk magnetic susceptibility of the material; it can be a complicated function of orientation, temperature, state of stress, time scale of observation and applied field but is often treated as a scalar.
Certain materials can produce magnetic fields in the absence of external magnetic fields (i.e., they are permanent magnets). As we shall see later in the course, these so-called “spontaneous” magnetic moments are also the result of spins of electrons which, in some crystals, act in a coordinated fashion, thereby producing a net magnetic field. The resulting magnetization can be fixed by various mechanisms and can preserve records of ancient magnetic fields. This remanent magnetization forms the basis of the field of paleomagnetism and will be discussed at length in the rest of this class.
From the foregoing discussion, we see that B and H are closely related. In paleomagnetic practice, both B and H are referred to as the “magnetic field”. Strictly speaking, B is the induction and H is the field, but the distinction is often blurred. The relationship between B and H is given by:
| (1.4) |
where
o is a parameter known as “the permeability of free space”. In the SI system,
o has dimensions of henries per meter and is given by
o = 4
× 10-7H . m-1.
So far, we have derived magnetic units in terms of the Système International (SI). In practice, you will notice that in many laboratories and in the literature people frequently use what are known as cgs units, based on centimeters, grams and seconds. You may wonder why any fuss would be made over using meters as opposed to centimeters because the conversion is trivial. With magnetic units, however, the conversion is far from trivial and has been the source of confusion and many errors. So, in the interest of clearing things up, we will briefly outline the cgs approach to magnetic units.
The derivation of magnetic units in cgs is entirely different from SI. The approach we will take here (see Cullity, 1972) starts with the concept of a magnetic pole with strength p. By analogy to Coulomb’s law, the force between two poles p1,p2 instead of with current loops as we did for SI units. Coulomb’s Law states that the force between two charges (q1,q2) is:
| (1.5) |
where r is the distance between the two charges. In cgs units, the proportionality constant k is simply unity, whereas in SI units it is 1 __ 4![]()
0 where
0 = 107 _ 4
c2
and c is the speed of light in a vacuum (hence
0 = 8.859 . 10-12 AsV-1m-1). [You can see why many people really prefer cgs but we are not allowed to publish in cgs in AGU journals so we just must grin and bear it!]
For magnetic units, we use pole strength p1,p1 in units of “electrostatic units” or esu, so Equation 1.5 becomes
Force in cgs is in units of dynes (dyn) so,
A magnetic pole, as an isolated electric charge, will create a magnetic induction
oH in the space around it. One unit of field strength (defined as one “oersted” or Oe) is the unit of field strength that exerts a force of one dyne on a unit of pole strength. The relationship between force, pole and field is written as:
Returning to the lines of force idea developed for magnetic fields earlier, let us define the oersted to be 1 line of force per square centimeter. Imagine a sphere with a radius r surrounding the magnetic monopole. The surface area of such a sphere is 4
r2. The sphere is a unit sphere (r = 1), the field strength at the surface is 1 Oe, then there must be 4
lines of force passing through it.
Proceeding to the notion of magnetic moment, from a cgs point of view, we start with a magnet of length l with two poles of strength p at each end. Placing the magnet in a field
oH, we find that it experiences a torque
proportional to p,l,H such that
| (1.6) |
Recalling our earlier discussion of magnetic moment, you will realize that pl is simply the magnetic moment m. The units of torque are energy, which are ergs in cgs, so the units of magnetic moment are ergs/oersted. We therefore define the “electromagnetic unit” (emu) as being one erg/oersted. [Some use emu to refer to the magnetization (volume normalized moment, see above), but this is incorrect.]
You will have noticed the use of the parameter
0 in the above treatment - a parameter missing in Cullity (1972) and in many books and articles using the cgs units. The reason for this is that
0 is unity in cgs units and simply converts from oersteds (H) and gauss (B) which are therefore used interchangeably. It was inserted in this derivation to remind us that there IS a difference and that the difference becomes very
important when we convert to SI because
0 is not unity, but 4
x 10-7 ! For conversion between commonly used cgs and SI paramters, please refer to Table 1. 1.
1 H = kg m2A-2s-2, 1 emu = 1 G cm3, B = |
An isolated electrical charge produces electrical fields that begin at the source (the charge) and spread (diverge) outward (see Figure 1.5a). Because there is no return flux to an oppositely charged “sink”, there is a net flux out of the dashed box shown in the figure. The “divergence” of the electrical field is defined as
.E which quantifies the net flux (see supplement to Chapter 1 for more). In the case of the field around an electric charge, the divergence is non-zero.
Magnetic fields are different from electrical fields in that there is no equivalent to an isolated electrical charge; there are only pairs of “opposite charges”, or magnetic dipoles. Therefore, any line of flux starting at one magnetic pole, returns to its sister pole and there is no net flux out of the box shown in Figure 1.5b; the magnetic field has no divergence (Figure 1.5b). This property of magnetic fields is another of Maxwell’s equations:
. B = 0.
|
|
We have already seen that the curl of the magnetic field (
×H) depends on the current density which is not always zero. Therefore, magnetic fields cannot generally be represented as the gradient of a scalar field. However, in the special case away from electric currents, the magnetic field can be written as the gradient of a scalar field that is known as the magnetic potential
m, i.e.,
|
|
The presence of a magnetic moment m creates a magnetic field which is the gradient of a scalar field. We also know that the divergence of the magnetic field is zero, hence
2
m = 0. This is LaPlace’s equation which is the starting point for spherical harmonic analysis discussed briefly in Lecture 2.
The magnetic potential
m is a function the vector r with radial distance r and angle
from the moment. Given a dipole moment m, the solution to LaPlace’s equation for the simple case of a magnetic field produced by m is:
| (1.7) |
The radial and tangential components of H at P (Figure 1.6) are:
Maxwell’s equations tell us that electric and changing magnetic fields are closely linked and can effect each other. Moving an electrical conductor through a magnetic field will cause electrons to flow, generating an electrical current. This is the principal of electric motors. In Figure 1.7 we see a design for a machine that will turn mechanical energy into magnetic field energy. The rotating disk is made of metal. As the disk turns in the presence of an initial magnetic field, the electrons scurry at right angles to the field, generating an electric potential (Figure 1.7b). The brush connection allows a current to flow through the wire wound into a coil, in turn generating a magnetic field. If the rotating disk is spun in the right direction, the magnetic field will be in the same sense as the initial field, amplifying the effect and generating a much larger magnetic field. More complicated setups using two disks whose fields interact with one another generate chaotic magnetic behavior that can switch polarities even if the mechanical motion remains steady. While a very poor analogue for the Earth’s magnetic field, it demonstrates that moving electrical conductors can generate a magnetic field. In the Earth of course the moving electrical conductor is the molten iron outer core.
|
|
In this appendix we will review the basic math concepts necessary to understand the chapter on magnetism. We will start with basic vector math and then review useful operators grad, div and curl.
|
|
To add the two vectors (see Figure A1) A and B, we break each vector into components Ax,y and Bx,y. For example, Ax = |A|cos
,Ay = |A|sin
where |A| is the length of the vector A. The components of the resultant vector C are: Cx = Ax + Bx,Cy = Ay + By. These can be converted back to polar coordinates of magnitude and angles if desired.
To subtract two vectors, compute the components as in addition, but the components of the vector difference C are: Cx = Ax - Bx,Cy = Ay - By.
There are two ways to multiply vectors. The first is the dot product whereby A.B = AxBx + AyBy. This is a scalar and is actually the cosine of the angle between the two vectors if the A and B are taken as unit vectors (assume a magnitude of unity in the component calculation.
|
|
The other way to perform vector multiplication is the cross product (see Figure A2), which produces a vector orthogonal to both A and B and whose components are given by:

To calculate the determinant, we follow these rules:
or
In paleomagnetism, we often have to change coordinate systems, from say sample coordinates to geographic, or to correct for tilting of the geological units. The way to do this in a simple 2-D case is illustrated in Figure A3. Given the vector shown in Figure A3a, that is oriented at an angle
from the X1 axis. To change to a second set of axes X'1,X'2,
we first have to define a set of coefficients called “direction cosines”. For example, the direction cosine a12 is the cosine of the angle between the old X1 and the new X'2,
12. We can define four of these direction cosines to fully describe the relationship between the two coordinate systems:
The first subscript always refers to the old system and the second refers to the new.
To find the new coordinates x'1,x'2 from the old, we just have:
In three dimensions we have:

A short cut notation to this is: x'i = aijxj. This just means that for each axis i, just sum through the j’s for all the dimensions.
|
|
We often wish to differentiate a function along three orthogonal axes. For example, imagine we know the topography of a ski area (see Figure B1). For every location (in say, X and Y coordinates), we know the height above sea level. This is a scalar function. Now imagine we want to build a ski resort, so we need to know the direction of steepest descent and the slope (red arrows in Figure B1).
|
|
To convert the scalar field (height versus position) to a vector field (direction and magnitude of greatest slope) mathematically, we would simply differentiate the topography function. Let’s say we had a very weird two dimensional, sinusoidal topography such that z = f(x) = sinx with z the height and x is the distance from some marker. The slope in the x direction (
), then would be
d _ dxf(x).
If f(x,y,z) were a three dimentional topography then the gradient of the topography function would be:
For short hand, we define a “vector differential operator” to be a vector whose components are
This can also be written in polar coordinates:
Just as the direction and magnitude of maximum slope of the topography is a the gradient of the scalar function of height, the magnetic field is the gradient of a scalar function of something we will define as the magnetic potential. In Lecture 1, we said that the magnetic field H is the gradient of a scalar potential field
m, or
The divergence of a vector function (e.g. H) is written as:
The trick here is to treat
as a vector and use the rules for dot products described in the section A of this appendix. In cartesian coordinates, this is:
Like all dot products, the divergence of a vector function is a scalar.
The name divergence is well chosen because
.H is a measure of how much the vector field “spreads out” (diverges) from the point in question. In fact, what divergence quantifies is the balance between vectors coming in to a particular region versus those that go out. The example in Figure B2 depicts a vector function whereby the magnitude of the vector increases linearly with distance away from the central point. An example of such a function would be v(r) = r. The divergence of this function
is:
Now consider Figure B3, which depicts a vector function that is constant over space, i.e. v(r) = k. The divergence of this function is:
The zero divergence means that for every vector leaving the box, there is an equal and opposite vector coming in. Put another way, no net flux results in a zero divergence. The fact that the divergence of the magnetic field is zero means that there are no point sources (monopoles), as opposed to electrical fields that have divergence related to the presence of electrons or protons.
The curl of the vector function B is defined as
× B. In cartesian coordinates we have
Curl is a measure of how much the vector function “curls” around a given point. The function describing the velocity of water in a whirlpool has a significant curl, while that of a smoothly flowing stream does not.
Consider Figure B4 which depicts a vector function v = -y
+ x
. The curl of this function is:

or
The magnetic field has a non-zero curl in the presence of currents or changing electric fields. In free space, away from currents (lightning!!), the magnetic field has zero curl.
Butler, R. F. (1992), Paleomagnetism: Magnetic Domains to Geologic Terranes, Blackwell Scientific Publications.
Cullity, B. (1972), Introduction to Magnetic Materials, Addison-Wesley Publishing Company.
Jiles, D. C. (1992), ‘Numerical determination of hysteresis parameters for the modeling of magnetic properties using the theory of ferromagnetic hysteresis’, IEEE Transactions on Magnetics 28, 27-35.
|
[ Home ]
[ Databases ]
[ Events ]
[ Tools ]
[ Publications ]
[ Links ]
[ Google Site Search ]
[ Metadata ] [ Who's Who ] [ Browsers and Plugins ] [ Database Indexes ] | |
| This page was last updated on 04-Apr-2008 Sponsored by NSF EAR 0000998 |
Supported by the San Diego Supercomputer Center and the Scripps Institution of Oceanography |