3 The Effect of Dielectric and Magnetic Media On Electric and Magnetic Fields
3 The Effect of Dielectric and Magnetic Media On Electric and Magnetic Fields
3 The Effect of Dielectric and Magnetic Media On Electric and Magnetic Fields
3.1 Polarization
where hpi i is the average dipole moment of the ith type of molecule in the vicinity
of point r, and Ni is the average number of such molecules per unit volume at r.
It is easily demonstrated that any divergence of the polarization field P (r)
gives rise to an effective charge density ρb in the medium. In fact,
ρb = −∇·P . (3.2)
This charge density is attributable to bound charges (i.e., charges which arise
from the polarization of neutral atoms), and is usually distinguished from the
charge density ρf due to free charges, which represents a net surplus or deficit of
electrons in the medium. Thus, the total charge density ρ in the medium is
ρ = ρ f + ρb . (3.3)
It must be emphasized that both terms in this equation represent real physical
charge. Nevertheless, it is useful to make the distinction between bound and
free charges, especially when it comes to working out the energy associated with
electric fields in dielectric media.
84
Gauss’ law takes the differential form
ρ ρf + ρ b
∇·E = = . (3.4)
0 0
This expression can be rearranged to give
∇·D = ρf , (3.5)
where
D = 0 E + P (3.6)
is termed the electric displacement, and has the same dimensions as P (dipole
moment per unit volume). The divergence theorem tells us that
I Z
D·dS = ρf dV. (3.7)
S V
In other words, the flux of D out of some closed surface S is equal to the total
free charge enclosed within that surface. Unlike the electric field E (which is the
force acting on unit charge) or the polarization P (the dipole moment per unit
volume), the electric displacement D has no clear physical meaning. The only
reason for introducing it is that it enables us to calculate fields in the presence
of dielectric materials without first having to know the distribution of polarized
charges. However, this is only possible if we have a constitutive relation connecting
E and D. It is conventional to assume that the induced polarization P is directly
proportional to the electric field E, so that
P = 0 χe E, (3.8)
where χe is termed the electric susceptibility of the medium. It follows that
D = 0 E, (3.9)
where
= 1 + χe (3.10)
is termed the dielectric constant or relative permittivity of the medium. (Likewise,
0 is termed the permittivity of free space.) It follows from Eqs. (3.5) and (3.9)
that
ρf
∇·E = . (3.11)
0
85
Thus, the electric fields produced by free charges in a dielectric medium are anal-
ogous to those produced by the same charges in a vacuum, except that they are
reduced by a factor . This reduction can be understood in terms of a polarization
of the atoms or molecules of the dielectric medium that produces electric fields
in opposition to that of given charge. One immediate consequence is that the
capacitance of a capacitor is increased by a factor if the empty space between
the electrodes is filled with a dielectric medium of dielectric constant (assuming
that fringing fields can be neglected).
It must be understood that Eqs. (3.8)–(3.11) are just an approximation which
is generally found to hold under terrestrial conditions (provided that the fields are
not too large) for isotropic media. For anisotropic media (e.g., crystals) Eq. (3.9)
generalizes to
D = 0 ·E, (3.12)
where is a second rank tensor known as the dielectric tensor. For strong elec-
tric fields D ceases to vary linearly with E. Indeed, for sufficiently strong elec-
tric fields neutral molecules are disrupted and the whole concept of a dielectric
medium becomes meaningless.
When the space near a set of charges contains dielectric material of non-uniform
dielectric constant then the electric field no longer has the same form as in vac-
uum. Suppose, for example, that the space is occupied by two dielectric media
whose uniform dielectric constants are 1 and 2 . What are the matching condi-
tions on E and D at the boundary between the two media?
Imagine a Gaussian pill-box enclosing part of the boundary surface between
the two media. The thickness of the pill-box is allowed to tend towards zero, so
that the only contribution to the outward flux of D comes from its flat faces.
These faces are parallel to the bounding surface and lie in each of the two media.
Their outward normals are dS1 (in medium 1) and dS2 , where dS1 = −dS2 .
Assuming that there is no free charge inside the disk (which is reasonable in the
86
limit where the volume of the disk tends towards zero), then Eq. (3.7) yields
where n21 is the normal to the boundary surface, directed from medium 1 to
medium 2. If the fields and charges are non time varying then Maxwell’s equa-
tions yield ∇∧E = 0, which give the familiar boundary condition (obtained by
integrating around a small loop which straddles the boundary surface)
In other word, the normal component of the electric displacement and the tan-
gential component of the electric field are both continuous across any boundary
between two dielectric materials.
87
z=0
ε2 P
ε1
R2 R1
/
A A
z ->
q / d d q
1 Ez (z = 0+ ) = 2 Ez (z = 0− ), (3.19a)
Ex (z = 0+ ) = Ex (z = 0− ), (3.19b)
Ey (z = 0+ ) = Ey (z = 0− ). (3.19c)
In order to solve this problem we will employ a slightly modified form of the
well known method of images. Since ∇ ∧ E = 0 everywhere, the electric field can
be written in terms of a scalar potential. So, E = −∇φ. Consider the region
z > 0. Let us assume that the scalar potential in this region is the same as that
obtained when the whole of space is filled with the dielectric 1 and, in addition to
the real charge q at position A, there is a second charge q 0 at the image position
A0 (see diagram). If this is the case then the potential at some point P in the
region z > 0 is given by
1 q q0
φ(z > 0) = + , (3.20)
4π0 1 R1 R2
p p
where R1 = ρ2 + (d − z)2 and R2 = ρ2 + (d + z)2 , when written in terms of
cylindrical polar coordinates (ρ, ϕ, z). Note that the potential (3.20) clearly is a
88
solution of Eq. (3.16) in the region z > 0. It gives ∇·E = 0, with the appropriate
singularity at the position of the point charge q.
Consider the region z < 0. Let us assume that the scalar potential in this
region is the same as that obtained when the whole of space is filled with the
dielectric 2 and a charge q 00 is located at the point A. If this is the case then the
potential in this region is given by
1 q 00
φ(z < 0) = . (3.21)
4π0 2 R1
The above potential is clearly a solution of Eq. (3.17) in the region z < 0. It gives
∇·E = 0, with no singularities.
It now remains to choose q 0 and q 00 in such a manner that the boundary
conditions (3.19) are satisfied. The boundary conditions (3.19b) and (3.19c) are
obviously satisfied if the scalar potential is continuous at the interface between
the two dielectric media:
The boundary condition (3.19a) implies a jump in the normal derivative of the
scalar potential across the interface:
∂φ(z = 0+ ) ∂φ(z = 0− )
1 = 2 . (3.23)
∂z ∂z
The first matching condition yields
q + q0 q 00
= , (3.24)
1 2
whereas the second yields
q − q 0 = q 00 . (3.25)
Here, use has been made of
∂ 1 ∂ 1 d
=− = 2 . (3.26)
∂z R1 z=0 ∂z R2 z=0 (ρ + d2 )3/2
89
Equations (3.24) and (3.25) imply that
0 2 − 1
q = − q, (3.27a)
2 + 1
22
q 00 = q. (3.27b)
2 + 1
∆χe = 1 − 2 . (3.28)
This implies that there is a polarization charge sheet on the interface between
the two dielectric media. In fact,
where n21 is a unit normal to the interface pointing from medium 1 to medium
2 (i.e., along the positive z-axis). Since
90
E0 E1 E0
z=0 z=b
91
Let us now generalize this result. Consider a dielectric medium whose dielec-
tric constant varies with z. The medium is assumed to be of finite extent and
is surrounded by a vacuum, so that (z) → 1 as |z| → ∞. Suppose that this
dielectric is placed in a uniform z-directed electric field E 0 . What is the field
E(z) inside the dielectric?
We know that the electric displacement inside the dielectric is given by D(z) =
0 (z) E(z). We also know from Eq. (3.5) that, since there are no free charges
and this is a one-dimensional problem,
dD(z) d[(z)E(z)]
= 0 = 0. (3.34)
dz dz
Furthermore, E(z) → E0 as |z| → ∞. It follows that
E0
E(z) = . (3.35)
(z)
Thus, the electric field is inversely proportional to the dielectric constant inside
the dielectric medium. The polarization charge density inside the dielectric is
given by
dE(z) d 1
ρb = 0 = 0 E0 . (3.36)
dz dz (z)
92
∂φ ∂φ
= . (3.38b)
∂θ r=a− ∂θ r=a+
Furthermore,
φ(r, θ, ϕ) → −E0 r cos θ (3.39)
as r → 0. Here, (r, θ, ϕ) are spherical polar coordinates centred on the sphere.
Let us search for an axisymmetric solution, φ = φ(r, θ). Since the solutions
to Poisson’s equation are unique, we know that if we can find such a solution
which satisfies all of the boundary conditions then we can be sure that this is the
correct solution. Equation (3.37) reduces to
1 ∂ 2 (rφ)
1 ∂ ∂φ
+ 2 sin θ = 0. (3.40)
r ∂r2 r sin θ ∂θ ∂θ
where l is a non-negative integer, the Al and Bl are arbitrary constants, and Pl (x)
is a solution to Legendre’s equation,
d dP l
(1 − x2 ) + l(l + 1) Pl = 0, (3.42)
dx dx
1 dl 2
Pl (x) = l l
(x − 1)l . (3.43)
2 l! dx
Since Eq. (3.42) is a Sturm-Liouville type equation, and the Legendre polynomi-
als satisfy Sturm-Liouville type boundary conditions at x = ±1, it immediately
93
follows that the Pl (cos θ) are orthogonal functions which form a complete set in
θ-space. The orthogonality relation can be written
Z 1
2
Pl0 (x)Pl (x) dx = δll0 . (3.44)
−1 2l + 1
where we have rejected the r −(l+1) radial solutions because they diverge unphys-
ically as r → 0. In the latter region
∞
X
φ(r, θ) = (Bl rl + Cl r−(l+1) )Pl (cos θ). (3.46)
l=0
However, it is clear from the boundary condition (3.39), and Eq. (3.43), that the
only non-vanishing Bl is B1 = −E0 . This follows since P1 (cos θ) = cos θ. The
boundary condition (3.38b) (which integrates to give φ(r = a − ) = φ(r = a+ ) for
a potential which is well behaved in θ-space) gives
C1
A1 = −E0 + , (3.47)
a3
and
Cl
Al = (3.48)
a2l+1
for l 6= 1. Note that it is appropriate to match the coefficients of the P l (cos θ)
since these functions are orthogonal. The boundary condition (3.38a) yields
C1
A1 = −E0 − 2 , (3.49)
a3
94
and
Cl
l Al = −(l + 1) (3.50)
a2l+1
for l 6= 1. Equations (3.48) and (3.50) give A l = Cl = 0 for l 6= 1. Equations (3.47)
and (3.49) reduce to
3
A1 = − E0 , (3.51a)
2+
−1
C1 = a 3 E0 . (3.51b)
+2
The solution is therefore
3
φ=− E0 r cos θ (3.52)
2+
for r ≤ a, and
a3
−1
φ = −E0 r cos θ + E0 cos θ (3.53)
+2 r2
for r > a.
Equation (3.52) is the potential of a uniform z-directed electric field of strength
3
E1 = E0 . (3.54)
2+
Note that E1 < E0 provided that > 1. Thus, the electric field strength is re-
duced inside the dielectric sphere due to partial shielding by polarization charges.
Outside the sphere the potential is equivalent to that of the applied field E 0 plus
the field of a point electric dipole, located at the origin and pointing in the z-
direction, whose dipole moment is
−1
p = 4π0 a 3 E0 . (3.55)
+2
This dipole moment can be interpreted as the volume integral of the polarization
P over the sphere. The polarization is
−1
P = 0 ( − 1) E1 ẑ = 30 E0 ẑ. (3.56)
+2
95
Since the polarization is uniform there is zero polarization charge density inside
the sphere. However, there is a polarization charge sheet on the surface of the
sphere whose density is given by σpol = P ·r̂ (see Eq. (3.29) ). It follows that
−1
σpol = 30 E0 cos θ. (3.57)
+2
where the integral is taken over all space and φ(r) is the electrostatic potential.
Here, it is assumed that the original charges and the dielectric are held fixed, so
96
that no mechanical work is performed. It follows from Eq. (3.5) that
Z
δU = φ ∇·δD d3 r, (3.61)
where δD is the change in the electric displacement associated with the charge
increment. Now the above equation can also be written
Z Z
δU = ∇·(φ δD) d r − ∇φ·δD d3 r,
3
(3.62)
giving Z Z
δU = φ δD·dS − ∇φ·δD d3 r, (3.63)
where use has been made of Gauss’s theorem. If the dielectric medium is of finite
spatial extent then we can neglect the surface term to give
Z Z
δU = − ∇φ·δD d r = E ·δD d3 r.
3
(3.64)
97
3.7 The force density within a dielectric medium
Equation (3.67) was derived by considering a virtual process in which true charges
are added to a system of charges and dielectrics which are held fixed, so that no
mechanical work is done against physical displacements. Let us now consider
a different virtual process in which the physical coordinates of the charges and
dielectric are given a virtual displacement δr at each point in space, but no free
charges are added to the system. Since we are dealing with a conservative system,
the energy expression (3.67) can still be employed, despite the fact that it was
derived in terms of another virtual process. The variation in the total electrostatic
energy δU when the system undergoes a virtual displacement δr is related to the
electrostatic force density f acting within the dielectric medium via
Z
δU = − f ·δr d3 r. (3.69)
If the medium is moving with a velocity field u then the rate at which electrostatic
energy is drained from the E and D fields is given by
dU
Z
= − f ·u d3 r. (3.70)
dt
Let us now consider the energy increment due to both a change δρ f in the
free charge distribution and a change δ in the dielectric constant, caused by the
displacements. From Eq. (3.67)
1
Z
2
D δ(1/) + 2D·δD/ d3 r,
δU = (3.71)
20
or
0
Z Z
δU = − E δ d r + E ·δD d3 r.
2 3
(3.72)
2
Here, the first term represents the energy increment due to the change in dielec-
tric constant associated with the virtual displacements, whereas the second term
corresponds to the energy increment caused by displacements of the free charges.
The second term can be written
Z Z Z Z
E ·δD d r = − ∇φ·δD d r = φ ∇·δD d r = φ δρf d3 r,
3 3 3
(3.73)
98
where surface terms have been neglected. Thus, Eq. (3.72) implies that
Z
dU ∂ρf 0 2 ∂
= φ − E d3 r. (3.74)
dt ∂t 2 ∂t
99
Integrating the first term by parts and neglecting any surface contributions, we
obtain Z Z
− φ ∇·(ρf u) d r = ρf ∇φ·u d3 r.
3
(3.82)
Likewise,
0 2 d 0 d
Z Z
E ρm ∇·u d3 r = − ∇ E2 ρm ·u d3 r. (3.83)
2 dρm 2 dρm
Thus, Eq. (3.81) becomes
Z
dU 0 2 0 d
= −ρf E + E ∇ − ∇ E 2 ρm · u d3 r. (3.84)
dt 2 2 dρm
Comparing with Eq. (3.70), we see that the force density inside the dielectric is
given by
0 2 0 d
f = ρf E − E ∇ + ∇ E 2 ρm . (3.85)
2 2 dρm
The first term in the above equation is the standard electrostatic force density.
The second term represents a force which appears whenever an inhomogeneous
dielectric is placed in an electric field. The last term, known as the electrostriction
term, gives a force acting on a dielectric in an inhomogeneous electric field. Note
that the magnitude of the electrostriction force depends explicitly on the dielectric
equation of state of the material, through d/dρ m . The electrostriction term gives
zero net force acting on any finite region of dielectric if we can integrate over a
large enough portion of the dielectric that its extremities lie in a field free region.
For this reason the term is frequently omitted, since in the calculation of the total
forces acting on dielectric bodies it usually does not contribute. Note, however,
that if the electrostriction term is omitted an incorrect pressure variation within
the dielectric is obtained, even though the total force is given correctly.
Let us now investigate what a dielectric equation of state actually looks like.
Suppose that a dielectric medium is made up of identical molecules which develop
100
a dipole moment
p = α0 E (3.86)
when placed in an electric field E. The constant α is called the molecular polar-
izability. If N is the number density of such molecules then the polarization of
the medium is
P = N p = N α0 E, (3.87)
or
NA ρ m α
P = 0 E, (3.88)
M
where ρm is the mass density, NA is Avogadro’s number, and M is the molecular
weight. But, how does the electric field experienced by an individual molecule
relate to the average electric field in the medium? This is not a trivial question
since we expect the electric field to vary strongly (on atomic length-scales) inside
the dielectric.
Suppose that the dielectric is polarized with a mean electric field E 0 which is
uniform (on macroscopic length-scales) and directed along the z-axis. Consider
one of the molecules which constitute the dielectric. Let us draw a sphere of ra-
dius a about this particular molecule. This is intended to represent the boundary
between the microscopic and the macroscopic range of phenomena affecting the
molecule. We shall treat the dielectric outside the sphere as a continuous medium
and the dielectric inside the sphere as a collection of polarized molecules. Ac-
cording to Eq. (3.29) there is a polarization surface charge of magnitude
on the inside of the sphere, where (r, θ, ϕ) are spherical polar coordinates, and
P = P ẑ = 0 ( − 1)E0 ẑ is the uniform polarization of the dielectric. The
magnitude of Ez at the molecule due to the surface charge is
1 σpol cos θ
Z
Ez = − dS, (3.90)
4π0 a2
where dS = 2πa2 sin θ dθ is a surface element of the sphere. It follows that
Z π
P P
Ez = cos2 θ sin θ dθ = . (3.91)
20 0 30
101
It is easily demonstrated that Eθ = Eϕ = 0 at the molecule. Thus, the field at
the molecule due to the surface charges on the sphere is
P
E= . (3.92)
30
The field due to the individual molecules within the sphere is obtained by
summing over the dipole fields of these molecules. The electric field at a distance
r from a dipole p is
1 p 3(p·r)r
E=− − . (3.93)
4π0 r3 r5
It is assumed that the dipole moment of each molecule within the sphere is the
same, and also that the molecules are evenly distributed throughout the sphere.
This being the case, the value of Ez at the molecule due to all of the other
molecules within in the sphere,
3(px xz + py yz + pz z 2 )
1 X pz
Ez = − − , (3.94)
4π0 r3 r5
mols
is zero, since
X X X 1X 2
x2 = y2 = z2 = r (3.95)
3
mols mols mols mols
and X X X
xy = yz = zx = 0. (3.96)
mols mols mols
It is easily seen that Eθ = Eϕ = 0. Hence, the electric field at the molecule due
to the other molecules within the sphere vanishes.
It is clear that the net electric field seen by an individual molecule is
P
E = E0 + . (3.97)
30
This is larger than the average electric field E0 in the dielectric. The above
analysis indicates that this effect is ascribable to the long range (rather than the
102
short range) interactions of the molecule with the other molecules in the medium.
Making use of Eq. (3.88) and the definition P = 0 ( − 1)E0 , we obtain
−1 NA ρ m α
= . (3.98)
+2 3M
This is called the Clausius-Mossotti relation. This formula is found to work pretty
well for a wide class of dielectric liquids and gases. The Clausius-Mossotti relation
yields
d ( − 1)( + 2)
= . (3.99)
dρm 3ρm
∇p = f . (3.100)
where 1 and 2 refer to two general points in the liquid. Here, it is assumed that
the liquid possesses an equation of state, so that p = p(ρ m ). If the liquid is
essentially incompressible (ρm ' constant) then
2
0 ρm d
p2 − p 1 = E2 . (3.103)
2 dρm 1
103
Finally, if the liquid obeys the Clausius-Mossotti relation then
2
0 E 2 ( − 1)( + 2)
p2 − p 1 = . (3.104)
2 3 1
104
that the electrostriction term makes no contribution to the line integral (3.106).
It follows that
0 B 2
Z
pA − p B = E ∇·dl. (3.107)
2 A
The only contribution to this integral comes from the vacuum/dielectric interface
in the vicinity of point A (since is constant inside the liquid, and E = 0 in the
vicinity of point B). Suppose that the electric field at point A has normal and
tangential (to the surface) components En and Et , respectively. Making use of
the boundary conditions that Et and En are constant across a vacuum/dielectric
interface, we obtain
Z
0 2 2 2 d
pA − p B = Et ( − 1) + En () 2
, (3.108)
2 1
giving
2
0 ( − 1) E
pA − p B = Et 2 + n . (3.109)
2
This electrostatic pressure difference can be equated to the hydrostatic pressure
difference ρm g h to determine the height h that the liquid rises between the plates.
At first sight, the above analysis appears to suggest that the dielectric liquid
is drawn upward by a surface force acting on the vacuum/dielectric interface
in the region between the plates. In fact, this is far from being the case. A
brief examination of Eq. (3.104) shows that this surface force is actually directed
downwards. According to Eq. (3.85), the force which causes the liquid to rise
between the plates is a volume force which develops in the region of non-uniform
electric field at the base of the capacitor, where the field splays out between
the plates. Thus, although we can determine the height to which the fluid rises
between the plates without reference to the electrostriction force, it is, somewhat
paradoxically, this force which is actually responsible for supporting the liquid
against gravity.
Let us consider another paradox concerning the electrostatic forces exerted in
a dielectric medium. Suppose that we have two charges embedded in a uniform
dielectric . The electric field generated by each charge is the same as that in
vacuum, except that it is reduced by a factor . Therefore, we expect that the
force exerted by one charge on another is the same as that in vacuum, except
105
that it is also reduced by a factor . Let us examine how this reduction in
force comes about. Consider a simple example. Suppose that we take a parallel
plate capacitor and insert a block of solid dielectric between the plates. Suppose,
further, that there is a small vacuum gap between the faces of the block and each
of the capacitor plates. Let ±σ be the surface charge densities on each of the
capacitor plates, and let ±σp be the polarization charge densities which develop on
the outer faces of the intervening dielectric block. The two layers of polarization
charge produce equal and opposite electric fields on each plate, and their effects
therefore cancel each other. Thus, from the point of view of electrical interaction
alone there would appear to be no change in the force exerted by one capacitor
plate on the other when a dielectric slab is placed between them (assuming that
σ remains constant during this process). That is, the force per unit area (which
is attractive) remains
σ2
fs = . (3.110)
20
However, in experiments in which a capacitor is submerged in a dielectric liquid
the force per unit area exerted by one plate on another is observed to decrease to
σ2
fs = . (3.111)
20
This apparent paradox can be explained by taking into account the difference
in liquid pressure in the field filled space between the plates and the field free
region outside the capacitor. This pressure difference is balanced by internal
elastic forces in the case of the solid dielectric discussed earlier, but is transmitted
to the plates in the case of the liquid. We can compute the pressure difference
between a point A on the inside surface of one of the capacitor plates and a point
B on the outside surface of the same plate using Eq. (3.107). If we neglect end
effects then the electric field is normal to the plates in the region between the
plates and is zero everywhere else. Thus, the only contribution to the line integral
(3.107) comes from the plate/dielectric interface in the vicinity of point A. Using
Eq. (3.109), we find that
2
0 1 σ 1
pA − p B = 1− E2 = 1− , (3.112)
2 20
106
where E is the normal field strength between the plates in the absence of dielectric.
The sum of this pressure force and the purely electrical force (3.110) yields a net
attractive force per unit area
σ2
fs = (3.113)
20
acting between the plates. Thus, any decrease in the forces exerted by charges
on one another when they are immersed or embedded in some dielectric medium
can only be understood in terms of mechanical forces transmitted between these
charges by the medium itself.
3.10 Magnetization
All matter is built up out of atoms, and each atom consists of electrons in mo-
tion. The currents associated with this motion are termed atomic currents. Each
atomic current is a tiny closed circuit of atomic dimensions, and may therefore
be appropriately described as a magnetic dipole. If the atomic currents of a given
atom all flow in the same plane then the atomic dipole moment is directed normal
to the plane (in the sense given by the right-hand rule) and its magnitude is the
product of the total circulating current and the area of the current loop. More
generally, if j(r) is the atomic current density at the point r then the magnetic
moment of the atom is
1
Z
m= r ∧ j d3 r, (3.114)
2
where the integral is over the volume of the atom. If there are N such atoms or
molecules per unit volume then the magnetization M (i.e., the magnetic dipole
moment per unit volume) is given by M = N m. More generally,
X
M (r) = Ni hmi i, (3.115)
i
where hmi i is the average magnetic dipole moment of the ith type of molecule in
the vicinity of point r, and Ni is the average number of such molecules per unit
volume at r.
Consider a general medium which is made up of molecules which are polar-
izable and possess a net magnetic moment. It is easily demonstrated that any
107
circulation in the magnetization field M (r) gives rise to an effective current den-
sity jm in the medium. In fact,
jm = ∇ ∧ M . (3.116)
This current density is called the magnetization current density, and is usually
distinguished from the true current density jt , which represents the convection
of free charges in the medium. In fact, there is a third type of current called a
polarization current, which is due to the apparent convection of bound charges.
It is easily demonstrated that the polarization current density j p is given by
∂P
jp = . (3.117)
∂t
Thus, the total current density j in the medium is given by
∂P
j = jt + ∇ ∧ M + . (3.118)
∂t
It must be emphasized that all terms on the right-hand side of this equation
represent real physical currents, although only the first term is due to the motion
of real charges (over more than atomic dimensions).
The Ampère-Maxwell equation takes the form
∂E
∇ ∧ B = µ 0 j + µ 0 0 , (3.119)
∂t
which can also be written
∂D
∇ ∧ B = µ 0 jt + µ 0 ∇ ∧ M + µ 0 , (3.120)
∂t
where use has been made of the definition D = 0 E + P . The above expression
can be rearranged to give
∂D
∇ ∧ H = jt + , (3.121)
∂t
where
B
H= −M (3.122)
µ0
108
is termed the magnetic intensity, and has the same dimensions as M (i.e., mag-
netic dipole moment per unit volume). In a steady-state situation, Stokes’s the-
orem tell us that I Z
H ·dl = jt ·dS. (3.123)
C S
In other words, the line integral of H around some closed curve is equal to
the flux of true current through any surface attached to that curve. Unlike the
magnetic field B (which specifies the force e v ∧ B acting on a charge e moving
with velocity v) or the magnetization M (the magnetic dipole moment per unit
volume), the magnetic intensity H has no clear physical meaning. The only
reason for introducing it is that it enables us to calculate fields in the presence of
magnetic materials without first having to know the distribution of magnetization
currents. However, this is only possible if we possess a constitutive relation
connecting B and H.
M = χm H, (3.124)
B = µH, (3.125)
10 Data obtained from the Handbook of Chemistry and Physics, Chemical Rubber Company
Press, Baca Raton, FL
109
Material χm
Aluminium 2.3 × 10−5
Copper −0.98 × 10−5
Diamond −2.2 × 10−5
Tungsten 6.8 × 10−5
Hydrogen (1 atm) −0.21 × 10−8
Oxygen (1 atm) 209.0 × 10−8
Nitrogen (1 atm) −0.50 × 10−8
where
µ = µ0 (1 + χm ) (3.126)
is termed the magnetic permeability of the material in question. (Likewise, µ 0
is termed the permeability of free space.) It is clear from Table 1 that the per-
meabilities of common diamagnetic and paramagnetic materials do not differ
substantially from that of free space. In fact, to all intents and purposes the
magnetic properties of such materials can be safely neglected (i.e., µ = µ 0 ).
3.12 Ferromagnetism
There is, however, a third class of magnetic materials called ferromagnetic ma-
terials. Such materials are characterized by a possible permanent magnetization,
and generally have a profound effect on magnetic fields (i.e., µ/µ 0 1). Unfor-
tunately, ferromagnetic materials do not exhibit a linear dependence between M
and H or B and H, so that we cannot employ Eqs. (3.124) and (3.125) with
constant values of χm and µ. It is still expedient to use Eq. (3.125) as the defi-
nition of µ, with µ = µ(H), however this practice can lead to difficulties under
certain circumstances. The permeability of a ferromagnetic material, as defined
by Eq. (3.125), can vary through the entire range of possible values from zero
to infinity, and may be either positive or negative. The most sensible approach
is to consider each problem involving ferromagnetic materials separately, try to
110
Figure 2: Magnetization curve and relative permeability of commercial iron (an-
nealed)
determine which region of the B-H diagram is important for the particular case
in hand, and then make approximations appropriate to this region.
First, let us consider an unmagnetized sample of ferromagnetic material. If
the magnetic intensity, which is initially zero, is increased monotonically, then
the B-H relationship traces out a curve such as that shown in Fig. 2. This is
called a magnetization curve. It is evident that the permeabilities µ derived from
the curve (according to the rule µ = B/H) are always positive, and show a wide
range of values. The maximum permeability occurs at the “knee” of the curve.
In some materials this maximum permeability is as large as 10 5 µ0 . The reason
for the knee in the curve is that the magnetization M reaches a maximum value
in the material, so that
B = µ0 (H + M ) (3.127)
111
B
-H C H
112
second and fourth quadrants of the diagram and positive in the first and third
quadrants. The shape of the hysteresis loop depends not only on the nature of the
ferromagnetic material but also on the maximum value of H to which the material
is subjected. However, once this maximum value, Hmax , becomes sufficient to
produce saturation in the material the hysteresis loop does not change shape with
any further increase in Hmax .
Ferromagnetic materials are used either to channel magnetic flux (e.g., around
transformer circuits) or as sources of magnetic field (permanent magnets). For
use as a permanent magnet, the material is first magnetized by placing it in a
strong magnetic field. However, once the magnet is removed from the external
field it is subject to a demagnetizing H. Thus, it is vitally important that a
permanent magnet should possess both a large remanence and a large coercivity.
As will become clear later on, it is generally a good idea for the ferromagnetic
materials used to channel magnetic flux around transformer circuits to possess
small remanences and small coercivities.
What are the matching conditions for B and H at the boundary between two
media? The governing equations for a steady state situation are
∇·B = 0, (3.128)
and
∇ ∧ H = jt . (3.129)
Integrating Eq. (3.128) over a Gaussian pill-box enclosing part of the boundary
surface between the two media gives
where n21 is the unit normal to this surface directed from medium 1 to medium
2. Integrating Eq. (3.129) around a small loop which straddles the boundary
surface yields
(H2 − H1 ) ∧ n21 = 0, (3.131)
113
assuming that there is no true current sheet flowing in this surface. In general,
there is a magnetization current sheet flowing in the boundary surface whose
density is given by
Jm = n21 ∧ (M2 − M1 ), (3.132)
where M1 is the magnetization in medium 1 at the boundary, etc. It is clear that
the normal component of the magnetic field and the tangential component of the
magnetic intensity are both continuous across any boundary between magnetic
materials.
∇ ∧ H = 0. (3.133)
∇2 φm = −ρm , (3.136)
where
ρm = −∇·M . (3.137)
114
Thus, the magnetostatic field H is determined by Poisson’s equation. We can
think of ρm as an effective magnetic charge density. Of course, this magnetic
charge has no physical reality. We have only introduced it in order to make
the problem of the steady magnetic field generated by a set of permanent mag-
nets look formally the same as that of the steady electric field generated by a
distribution of charges.
The unique solution of Poisson’s equation, subject to sensible boundary con-
ditions at infinity, is well known:
1 ρm (r 0 ) 3 0
Z
φm (r) = d r. (3.138)
4π |r − r 0 |
This yields
1 ∇0 ·M (r 0 ) 3 0
Z
φm (r) = − d r. (3.139)
4π |r − r 0 |
If the magnetization field M (r) is well behaved and localized we can integrate
by parts to obtain
1 1
Z
φm (r) = M (r 0 )·∇0 d3 r 0 . (3.140)
4π |r − r |
0
Now
0 1 1
∇ = −∇ , (3.141)
|r − r 0 | |r − r 0 |
so our expression for the magnetic potential can be written
1 M (r 0 ) 3 0
Z
φm (r) = − ∇· d r. (3.142)
4π |r − r 0 |
115
magnetization asymptotically produces a dipole magnetic field whose strength is
determined by the net magnetic moment of the distribution.
It is often a good approximation to treat the magnetization field M (r) as a
discontinuous quantity. In other words, M (r) is specified throughout the “hard”
ferromagnets in question, and suddenly falls to zero at the boundaries of these
magnets. Integrating Eq. (3.137) over a Gaussian pill-box which straddles one
of these boundaries leads to the conclusion that there is an effective magnetic
surface charge density,
σm = n·M , (3.144)
on the surface of the ferromagnets, where M is the surface magnetization, and
n is a unit outward directed normal to the surface. Under these circumstances
Eq. (3.139) yields
1 ∇0 ·M (r 0 ) 3 0 1 M (r 0 )·dS 0
Z Z
φm (r) = − d r + , (3.145)
4π V |r − r 0 | 4π S |r − r 0 |
where V represents the volume occupied by the magnets and S is the bounding
surface to V . Here, dS is an outward directed volume element to S. It is clear
that Eq. (3.145) consists of a volume integral involving the volume magnetic
charges ρm = −∇·M and a surface integral involving the surface magnetic charges
σm = n·M . If the magnetization is uniform throughout the volume V then the
first term in the above expression vanishes and only the surface integral makes a
contribution.
We can also write B = ∇ ∧ A in order to satisfy ∇·B = 0 automatically. It
follows from Eqs. (3.121) and (3.122) that
∇ ∧ H = ∇ ∧ (B/µ0 − M ) = 0, (3.146)
which gives
∇2 A = −µ0 jm , (3.147)
since jm = ∇ ∧ M . The unique solution to Eq. (3.147), subject to sensible
boundary conditions at infinity, is very well known:
µ0 jm (r 0 ) 3 0
Z
A(r) = d r. (3.148)
4π |r − r 0 |
116
Thus,
µ0 ∇0 ∧ M (r 0 ) 3 0
Z
A(r) = d r. (3.149)
4π |r − r 0 |
∇2 φm = 0, (3.151)
on the surface of the sphere. One of the matching conditions at the surface of
the sphere is that the tangential component of H must be continuous. It follows
from Eq. (3.134) that the scalar magnetic potential must be continuous at r = a,
so that
φm (r = a+ ) = φm (r = a− ). (3.153)
117
Integrating Eq. (3.136) over a Gaussian pill-box straddling the surface of the
sphere yields
r=a+
∂φm
= −σm = −M0 cos θ. (3.154)
∂r r=a−
In other words, the magnetic charge sheet on the surface of the sphere gives rise
to a discontinuity in the radial gradient of the magnetic scalar potential at r = a.
The most general axisymmetric solution to Eq. (3.151) which satisfies physical
boundary conditions at r = a and r = ∞ is
∞
X
φm (r, θ) = Al rl Pl (cos θ) (3.155)
l=0
Bl = Al a2l+1 (3.157)
(l + 1) Bl
− l+2
− lAl al−1 = −M0 δl1 (3.158)
a
for all l, since Pl (cos θ) = cos θ. It follows that
Al = B l = 0 (3.159)
for l 6= 1, and
M0
A1 = , (3.160a)
3
M0 a3
B1 = . (3.160b)
3
118
Thus,
M0 a2 r
φm (r, θ) = cos θ (3.161)
3 a2
for r < a, and
M0 a2 a
φm (r, θ) = cos θ (3.162)
3 r2
for r ≥ a. Since there is a uniqueness theorem associated with Poisson’s equation,
we can be sure that this axisymmetric potential is the only solution to the problem
which satisfies physical boundary conditions at r = 0 and infinity.
In the vacuum region outside the sphere
where
4 3
m= πa M . (3.165)
3
This, of course, is the magnetic field of a magnetic dipole m. Not surprisingly,
the net dipole moment of the sphere is equal to the integral of the magnetization
M (which is the dipole moment per unit volume) over the volume of the sphere.
Inside the sphere we have H = −∇φm and B = µ0 (H + M ), giving
M
H=− , (3.166)
3
and
2
B= µ0 M . (3.167)
3
Thus, both the H and B fields are uniform inside the sphere. Note that the
magnetic intensity is oppositely directed to the magnetization. In other words,
the H field acts to demagnetize the sphere. How successful it is at achieving
119
µ = -2 µ 0
B
operating point
BR
B-H curve
0
-H
−µ0 HC 0
this depends on the shape of the hysteresis curve in the negative H and posi-
tive B quadrant. This curve is sometimes called the demagnetization curve of
the magnetic material which makes up the sphere. Figure 4 shows a schematic
demagnetization curve. The curve is characterized by two quantities: the reten-
tivity BR (i.e., the residual magnetic field strength at zero magnetic intensity)
and the coercivity µ0 Hc (i.e., the negative magnetic intensity required to demag-
netize the material: this quantity is conventionally multiplied by µ 0 to give it the
units of magnetic field strength). The operating point (i.e., the values of B and
µ0 H inside the sphere) is obtained from the intersection of the demagnetization
curve and the curve B = µH. It is clear from Eqs. (3.166) and (3.167) that
µ = −2 µ0 (3.168)
for a uniformly magnetized sphere in the absence of external fields. The magne-
tization inside the sphere is easily calculated once the operating point has been
determined. In fact, M0 = B − µ0 H. It is clear from Fig. 4 that for a magnetic
material to be a good permanent magnet it must possess both a large retentivity
and a large coercivity. A material with a large retentivity but a small coercivity
120
is unable to retain a significant magnetization in the absence of a strong external
magnetizing field.
with
3µ
B= B0 , (3.172)
µ + 2µ0
121
where, in general, µ = µ(B). Clearly, for a highly permeable material (i.e.,
µ/µ0 1, which is certainly the case for soft iron) the magnetic field strength
inside the sphere is approximately three times that of the externally applied field.
In other words, the magnetic field is amplified inside the sphere.
The amplification of the magnetic field by a factor three in the high perme-
ability limit is specific to a sphere. It can be shown that for elongated objects
(e.g., rods), aligned along the direction of the external field, the amplification
factor can be considerably larger than this.
It is important to realize that the magnetization inside a ferromagnetic mate-
rial cannot increase without limit. The maximum possible value of M is called the
saturation magnetization, and is usually denoted M s . Most ferromagnetic mate-
rials saturate when they are placed in external magnetic fields whose strengths
are greater than, or of order, one tesla. Suppose that our soft iron sphere first
attains the saturation magnetization when the unperturbed external magnetic
field strength is Bs . It follows from (3.170b) and (3.171) (with µ µ 0 ) that
B = B0 + 2Bs (3.173)
inside the sphere, for B0 > Bs . In this case, the field amplification factor is
B Bs
=1+2 . (3.174)
B0 B0
Thus, for B0 Bs the amplification factor approaches unity. We conclude that
if a ferromagnetic material is placed in an external field which greatly exceeds
that required to cause saturation then the material effectively loses its magnetic
properties, so that µ ' µ0 . Clearly, it is very important to avoid saturating the
soft magnets used to channel magnetic flux around transformer circuits. This sets
an upper limit on the magnetic field strengths which can occur in such circuits.
122
surrounding the region in question by a material of high permeability. It is vitally
important that a material used as a magnetic shield does not develop a permanent
magnetization in the presence of external fields, otherwise the material itself may
become a source of magnetic fields. The most effective commercially available
magnetic shielding material is called Mumetal, and is an alloy of 5% Copper, 2%
Chromium, 77% Nickel, and 16% Iron. The maximum permeability of Mumetal
is about 105 µ0 . This material also possesses a particularly low retentivity and
coercivity. Unfortunately, Mumetal is extremely expensive. Let us investigate
how much of this material is actually required to shield a given region from an
external magnetic field.
Consider a spherical shell of magnetic shielding, made up of material of per-
meability µ, placed in a formerly uniform magnetic field B 0 = B0 ẑ. Suppose
that the inner radius of the shell is a and the outer radius is b. Since there are
no free currents in the problem, we can write H = −∇φm . Furthermore, since
B = µH and ∇ · B = 0, it is clear that the magnetic scalar potential satisfies
Laplace’s equation, ∇2 φm = 0, throughout all space. The boundary conditions
are that the potential must be well behaved at r = 0 and r → ∞, and also that
the tangential and the normal components of H and B, respectively, must be
continuous at r = a and r = b. The boundary conditions on H merely imply that
the scalar potential φm must be continuous at r = a and r = b. The boundary
conditions on B yield
Let us try the following test solution for the magnetic potential:
B0 α
φm = − r cos θ + 2 cos θ (3.176)
µ0 r
for r > b, γ
φm = β r + 2 cos θ (3.177)
r
123
for b ≥ r ≥ a, and
φm = δ r cos θ (3.178)
for r < a. This potential is certainly a solution of Laplace’s equation throughout
space. It yields the uniform magnetic field B0 as r → ∞, and satisfies physical
boundary conditions at r = 0 and infinity. Since there is a uniqueness theorem
associated with Poisson’s equation, we can be certain that this potential is the
correct solution to the problem provided that the arbitrary constants α, β, etc.
can be adjusted in such a manner that the boundary conditions at r = a and
r = b are also satisfied.
The continuity of φm at r = a and r = b requires that
γ
β a + 2 = δ a, (3.179)
a
and
γ B0 α
βb+ 2 =− b + 2. (3.180)
b µ0 b
The boundary conditions (3.175) yield
2γ
µ0 δ = µ β − 3 , (3.181)
a
and
B0 2α 2γ
µ0 − − 3 =µ β− 3 . (3.182)
µ0 b b
It follows that
(2µ + µ0 )(µ − µ0 )
µ0 α = 3 3 2
(b3 − a3 )B0 ,
(2µ + µ0 )(µ + 2µ0 ) − 2(a /b )(µ − µ0 )
(3.183a)
3(2µ + µ0 )µ0
µ0 β = − B0 , (3.183b)
(2µ + µ0 )(µ + 2µ0 ) − 2(a3 /b3 )(µ − µ0 )2
3(µ − µ0 )µ0
µ0 γ = − 3 3 2
a 3 B0 , (3.183c)
(2µ + µ0 )(µ + 2µ0 ) − 2(a /b )(µ − µ0 )
9µµ0
µ0 δ = − B0 . (3.183d)
(2µ + µ0 )(µ + 2µ0 ) − 2(a3 /b3 )(µ − µ0 )2
124
Consider the limit of a thin, high permeability shell for which b = a + d,
d/a 1, and µ/µ0 1. In this limit, the field inside the shell is given by
3 µ0 a
B' B0 . (3.184)
2 µ d
Thus, if µ ' 105 µ0 for Mumetal, then we can reduce the magnetic field strength
inside the shell by almost a factor of 1000 using a shell whose thickness is only
1/100 th of its radius. Clearly, a little Mumetal goes a long way! Note, however,
that as the external field strength, B0 , is increased, the Mumetal shell eventually
saturates, and µ/µ0 gradually falls to unity. Thus, extremely strong magnetic
fields (typically, B0 ∼
>
1 tesla) are hardly shielded at all by Mumetal, or similar
magnetic materials.
jt = σ(E + E 0 ), (3.185)
where jt is the true current density, σ is the conductivity, and E is the inductive
electric field. Taking the scalar product with jt , we obtain
0 jt 2
E ·jt = − E ·jt . (3.186)
σ
The left-hand side of this equation represents the rate at which the battery does
work on the conductor. The first term on the right-hand side is the rate of Joule
heating inside the conductor. We tentatively identify the remaining term with
the rate at which energy is fed into the magnetic field. If all fields are quasi-
stationary (i.e., slowly varying) then the displacement current can be neglected,
and the Ampère-Maxwell equation reduces to ∇ ∧ H = jt . Substituting this
125
expression into Eq. (3.186) and integrating over all space, we get
2
(∇ ∧
Z Z Z
H)
E 0 ·(∇ ∧ H) d3 r = d3 r − E ·(∇ ∧ H) d3 r. (3.187)
σ
The last term can be integrated by parts using the relation
∇·(E ∧ H) = H ·(∇ ∧ E) − E ·(∇ ∧ H). (3.188)
The divergence theorem plus the Faraday-Maxwell equation yield
∂B
Z Z Z
E ·(∇ ∧ H) d3 r = − H · d3 r − (E ∧ H)·dS. (3.189)
∂t
Since E ∧ H falls off at least as fast as 1/r 5 in electrostatic and quasi-stationary
magnetic fields (1/r 2 comes from electric monopole fields, and 1/r 3 from magnetic
dipole fields), the surface integral in the above expression can be neglected. Of
course, this is not the case for radiation fields, for which E and H fall off like
1/r. Thus, the constraint of “quasi-stationarity” effectively means that the fields
vary sufficiently slowly that any radiation fields can be neglected.
The total power expended by the battery can now be written
2
(∇ ∧ ∂B 3
Z Z Z
H)
E 0 ·(∇ ∧ H) d3 r = d3 r + H · d r. (3.190)
σ ∂t
The first term on the right-hand side has already been identified as the energy
loss rate due to Joule heating. The last term is obviously the rate at which energy
is fed into the magnetic field. The variation δU in the magnetic field energy can
therefore be written Z
δU = H ·δB d3 r. (3.191)
This result is analogous to the result (3.64) for the variation in the energy of an
electrostatic field.
In order to make Eq. (3.191) integrable, we must assume a functional relation-
ship between H and B. For a medium which magnetizes linearly the integration
can be carried out in much the same manner as Eq. (3.67), to give
1
Z
U= H ·B d3 r. (3.192)
2
126
Thus, the magnetostatic energy density inside a linear magnetic material is given
by
H· B
W = . (3.193)
2
Unfortunately, most interesting magnetic materials, such as ferromagnets, exhibit
a nonlinear relationship between H and B. For such materials, Eq. (3.191) can
only be integrated between definite states, and the result, in general, depends on
the past history of the sample. For ferromagnets, the integral of Eq. (3.191) has
a finite, non-zero value when B is integrated around a complete magnetization
cycle. This cyclic energy loss is given by
Z I
∆U = H · dB d3 r. (3.194)
In other words, the energy expended per unit volume when a magnetic material
is carried through a magnetization cycle is equal to the area of its hysteresis loop
as plotted in a graph of B against H. Thus, it is particularly important to ensure
that the magnetic materials used to form transformer cores possess hysteresis
loops with comparatively small areas, otherwise the transformers are likely to be
extremely inefficient.
127