Magneto Statics
Magneto Statics
Magneto Statics
5.1 The Lorentz Force Law 5.1.1 Magnetic Fields Outside a magnetic the lines emerge from the north pole
and enter the south pole; within the magnet they are
By analogy with electrostatics, why don’t we study directed from the south pole to the north pole. The dots
magnetostatics first? Due to lack of magnetic monopole. represents the tip of an arrow coming toward you. The
cross represents the tail of an arrow moving away.
F = qvB sin θ = qv × B
The magnetic field, B, at a point along the tangent to a field
line. The direction of B is that of the force on the north pole of Since F is always perpendicular to v, a magnetic force does
a bar magnet, or the direction in which a compass needle no work on a particle and cannot be used to change its
points. The strength of the field is proportional to the number kinetic energy.
of lines passing through a unit area normal to the field (flux
density). The SI unit of magnetic field is the Tesla (T). 1 T=104 G
3 4
The Lorentz Force Law Force on a Current-Carrying Conductor
When a particle is subject to both electric and magnetic When a current flows in a magnetic field, the electrons as a
fields in the same region, what is the total force on it? whole acquire a slow drift speed, vd, and experience a
magnetic force, which is then transmitted to the wire.
F = q (E + v × B )
F = qvB sin θ = (nA e)vd B
This is called the Lorentz force law. This axiom is found in
experiments. = (nAevd ) B
=I B
The magnetic force do no work.
F = I ×B
dWmag = Fmag ⋅ dl = q ( v × B) ⋅ vdt = 0 F = I B sin θ
Really? But, how do you explain a magnetic crane lifts a n: is the number of the conductor per unit volume.
container? : is defined to be in the direction in which the current is flowing.
5 6
Example:
Force on a Current-Carrying Conductor The Magnetic Force on a Semicircular Loop
The force on an infinitesimal current element is A wire is bent into a semicircular loop of radius R. It carries a
current I, and its plane is perpendicular to a uniform magnetic
dF = Id × B field B, as shown below. Find the force on the loop.
Solution:
dF = Id × B
dFy = IRB sin θ dθ
π
Fy = IRB ∫ sin θ dθ
0
= 2 IRB = I (2 R) B
The force on a wire is the vector sum (integral) of the forces The x-components of the forces on such elements will
on all current elements. cancel in pairs.
The net force on any close current-carrying loop is zero.
7 8
The Motion of Charged Particles in Magnetic Fields Cyclotron Motion
How does a charged particle move with an initial velocity v What are the frequency and the period? Are they independent
perpendicular to a uniform magnetic field B? of the speed of the particle? Yes.
Since v and B are perpendicular, the particle experiences a The period of the orbit is
force F=qvB of constant magnitude directed perpendicular.
2πr 2πm m 2π
Under the action of such a force, the particle will move in a T= = =
circular path at constant speed. From Newton’s second law, v qB q B
F=ma, we have 1 qB q B
fc = = =
mv 2 mv T 2πm m 2π
qvB = ⇒ r=
r qB f c B = 2.8 MHz Gauss
The radius of the orbit is directly proportional to the linear The frequency is called the cyclotron frequency.
momentum of the particle and inversely proportional to the
magnetic field strength. All particles with the same charge-to-mass ratio, q/m, have
9 the same period and cyclotron frequency. 10
E
v=
E = − Ej E B1
⇒ v = i = vi
B = − Bk B mv mE
r= =
qB2 qB1 B2
Only those particles with speed v=E/B pass through the
crossed fields undeflected. This provides a convenient way m B1 B2
= r
of either measuring or selecting the velocities of charged q E
particles.
15 16
Example: Mass Spectrometer Current and Surface Current
In a mass spectrometer shown below, two isotopes of an The current in a wire is the charge per unit time passing a
elements with mass m1 and m2 are accelerated from rest by given point.
a potential difference V. They then enter a uniform B normal Current is measured in coulombs-per-second, or amperes (A).
to the magnetic field lines. What is the ratio of the radii of
their paths? 1 A=1 C/s
Solution:
The surface current density, K, is defined as follows:
2qV consider a ”ribbon” of infinitesimal width d ⊥, running parallel
v= to the flow. Then, dI
m K=
d ⊥
mv 2mV
r= = then r1 / r2 = (m1 / m2 ) In words, K is the current per unit width-perpendicular-to-flow.
qB qB 2
d (σ ⊥ // )
d( )
dI dt d //
Note1: How particle is accelerated by a potential difference? K= ≈ =σ =σv
17 d ⊥ d ⊥ dt 18
d ( ρ a⊥ dQ d
dt dt V∫
// ) where = ρ dτ .
d( )
dI dt d //
J= ≈ =ρ = ρv
da⊥ da⊥ dt d dρ
∫V (∇ ⋅ J )dτ = − dt V∫ ρ dτ ⇒ ∇ ⋅ J = − dt
19 continuity equation 20
5.2 The Biot-Savart Law 5.2.1 Steady Currents 5.2.2 The Magnetic Field of a Steady Current
Stationary charges produce electric fields that are constant The Biot-Sarvart law:
in time. Steady currents produce magnetic fields that are µ0 I × rˆ µ I dl ′ × rˆ
also constant in time. B(r ) = ∫
4π r 2
dl ′ = 0 ∫
4π r2
Stationary charges ⇒ constant electric fields; electrostatics. The integration is along the current
Steady currents ⇒ constant magnetic fields; magnetostatics. path, in the direction of the flow.
µ0: the permeability of free space.
Steady current means that a continuous flow that goes on Definition of magnetic field B: newtons per ampere-meter
forever without change and without charge piling up
anywhere. They represent suitable approximations as long
or tesla (T). 1 T=1 N/(A⋅m)
as the fluctuations are reasonably slow.
The Biot-Sarvart law plays a role analogous to Coulomb’s
∇⋅J = 0 law in electrostatics.
21 22
π
Con't
Example 5.5 Find the magnetic field a distance s from a long µ I dl′ × rˆ µ I 2
cos 2 θ s
B(r ) = 0 ∫ 2 = 0 ∫π cos θ dθ
straight wire carrying a steady current I. 4π r 4π s 2 cos 2 θ
−
2
π
µ0 I µI I
= sin θ 2π = 0 (= 2 ×10−7 Tesla)
4π s −
2 2π s s
µ0 I dl′ × rˆ µ0 K × rˆ
4π ∫ r 2
B(r ) = da ′
4π ∫ r 2
Sol : B(r ) = For surface current:
Solution:
Sine the solenoid is a series of closely
packed loops, we may divided into current
loops of width dz, each of which contains
ndz turns, where n =N/L is the number of
turns per unit length.
The current within such a loop is (ndz)I.
27 28
Solenoid (II)
Con’t Homework #9
z = a tan θ ⇒ dz = a sec θ dθ 2
1
= µ0 nI (sin θ 2 − sin θ1 )
2
B = µ 0 nI (infinite long solenoid)
29 30
∇ × B = µ0 J
∫ B(r) ⋅ dl = ∫ 2π0 s φˆ ⋅ φˆ sdφ = µ0 I
In fact for any loop that encloses the wire would give the Does this differential equation apply to any shape of the
same answer. Really? current loop? Yes, to be prove next.
31 32
5.3.2 The Divergence and Curl of B The Divergence of B
The divergence of B:
The Biot-Sarvart Law for the
µ0 J (r′) × rˆ µ J (r′) × rˆ
general case of a volume
current:
∇ ⋅ B(r ) = ∇ ⋅ (
4π ∫ r 2
dτ ′) = 0 ∫ ∇ ⋅ (
4π r2
)dτ ′
J (r′) × rˆ rˆ rˆ
∇⋅( 2
) = 2 ⋅ (∇ × J ) − J ⋅ (∇ × 2 )
r r r
µ0 J (r′) × rˆ
4π ∫ r 2
B(r ) = dτ ′ ∇ ⋅ ( A × B) = B ⋅ (∇ × A) − A ⋅ (∇ × B)
rˆ 1
2
= −∇( ) (Prob. 1.13)
r r
0 0
J (r′) × rˆ rˆ rˆ
The integration is over the primed coordinates. ∇⋅( 2
) = 2 ⋅ (∇ × J ) − J ⋅ (∇ × 2 )
r r r
The divergence and the curl are to be taken with ∴ ∇ ⋅ B = 0 The divergence of a magnetic field is zero.
respective to the unprimed coordinates.
33 34
µ 0I
B= r (r < R )
2πR 2
39 40
Solenoid Toroid
Example 5.9 An ideal infinite solenoid has n turns per unit Example 5.10 A toroidal coil (shaped like a doughnut) is
length and carries a current I. Find its magnetic field inside. tightly wound with N turns and carries a current I. We
assume that it has a rectangular cross section, as shown
Solution:
below. Find the field strength within the toroid.
∫B⋅d Solution:
b c d a
= ∫ B⋅d + ∫ B⋅d + ∫ B⋅d + ∫ B⋅d
a
b
b c d
∫B⋅d = B ∫ d = µ 0 NI
= ∫ B⋅d µ 0NI
B=
a
2πr
BLab = µ 0nLab I
The field is not uniform; it varies
as 1/r. The toroidal fields are
B = µ 0nI
used in research on fusion power.
41 42
∇ λ = −∇ ⋅ A 0 λ = 1 ∇ ⋅ A0
∫ dτ ′
2
4π r
43 It is always possible to make the vector potential divergenceless.44
The Vector Potential and Scalar Potential Example 5.11 A spherical shell, of radius R, carrying a
uniform surface charge σ, is set spinning at angular velocity
Using the Coulomb gauge, we obtain : ∇ 2 A = − µ0 J ω. Find the vector potential it produce at point r.
µ J (r′)
A= 0 ∫ dτ ′
4π r
For line and surface current,
µ I µ0 K (r′)
4π ∫ r
A = 0 ∫ dl ′ A= da′
4π r Sol :First, let the observer is in the z axis and ω is tilted at an angle ψ
µ K (r′)
What happens when the curl of B vanishes? Vector potential is A(r ) = 0 ∫ da′
4π r
Magnetostatic scalar potential. The surface current density K (r′) = σ v′
∇ × B = 0 ⇒ B = −∇U
⇒ ∇ 2U = 0 (similiar to Laplace's equation)
45 46
− µ0 R 3σω sinψ yˆ
1
u
=
2 ∫
−1 r + R − 2rRu
2 2
du
Sol : A cute method that does the job. We have derived five formulas
interrelating three fundamental
∫ B ⋅ da = Φ = ∫ (∇ × A) ⋅ da = ∫ A ⋅ dl quantities: J, A and B.
where Φ is the flux of Β through the loop in question.
∫ B ⋅ dl = µ I
0 enc ⇒ ∫ A ⋅ dl = Φ
Using a circular " amperian loop" at a radius inside the solenoid. Comments:
µ nI •There is one “missing link” in the diagram.
∫ A ⋅ dl = A2πs = ∫ B ⋅ da = µ0 nI (πs ) ⇒ A = 02 sφˆ for s < R
2
∫B⋅d = µ 0 I enc
Consider a wafer-thin pillbox. Gauss’s law states that Consider a thin rectangular loop. The
curl of the Ampere’s law states that P
∫ B ⋅ da = 0
S
The ends gives nothing (as ε→0), and the sides give
The sides of the pillbox contribute nothing to the flux, in the
limit as the thickness ε goes to zero.
//
( Babove − Bbelow
//
) = µ0 K ⇒ Babove
//
− Bbelow
//
= µ0 K
⊥ ⊥ ⊥ ⊥ In short, B above − B below = µ 0 K × nˆ , where nˆ points " upward."
( Babove − Bbelow ) A = 0 ⇒ Babove = Bbelow
How about the vector potential A?
51 52
Boundary Conditions in Terms of Vector Potential 5.4.3 Multipole Expansion of the Vector Potential
1 r′ r′
A above = A below = (1 + ( ) cos θ ′ + ( ) 2 ((3 cos 2 θ ′ − 1) / 2) + …)
r r r
1 ∞
r ′
= ∑ ( ) P (cosθ ′)
r =0 r
⊥ ⊥
∇ ⋅ A = 0 ⇒ Aabove = Abelow The vector potential of a current loop
µI 1 µI ∞ 1
A= 0 ∫ dl ′ = 0 ∑ n +1 ∫ (r ′) n Pn (cos θ ′)dl ′
∇ × A = B ⇒ ∫ A ⋅ dl = ∫ B ⋅ da = Φ ⇒ Aabove = Abelow
// //
4π r- r′ 4π n =0 r
µ I 1 1 1
= 0 ∫ dl ′ + 2 ∫ r ′ cos θ ′dl ′ + 3 ∫ r ′2 P2 (cos θ ′)dl ′ +
4π r r r
µ I 1 1 1
Part I Recalling Stokes’ theorem ∫ (∇ × v) ⋅ da = ∫ v ⋅ dl
S
A = 0 ∫ dl′ + 2 ∫ r′ cos θ ′dl′ + r 3 ∫ r′ P (cosθ ′)dl′ +
2
Let v=cT ∇ × ( fA) = ∇f × A + f (∇ × A) P
4π r
2
r
magnetic monopole term is always zero. ∫ (∇ × cT ) ⋅ da = ∫ (∇T × c + T (∇ × c)) ⋅ da = −c ⋅ ∫ ∇T × da
S S S
µ0 I 1 µI 1 ∫ cT ⋅ dl = c ⋅ ∫ Tdl ∫ ∇T × da = − ∫ Tdl
2 ∫
A dip = r ′ cosθ ′dl′ = 0 2 ∫ (rˆ ⋅ r′)dl′ S
4π r 4π r
P P P
µ0 m
B dip = ∇ × A = (2 cos θrˆ + sin θθˆ )
4πr 3
57 58
Chapter 6 Magnetic Fields in Matter 6.1.2. Torques and Forces on Magnetic Dipoles
6.1 Magnetization
6.1.1 Diamagnets, Paramagnets, Ferromagnets A magnetic dipole experiences a torque in a magnetic field,
just as an electric dipole does in an electric field.
All the magnetic phenomena are due to electric charges in
motion: Any current loop could be built up from infinitesimal
Electrons orbiting around nuclei rectangles, with all the “internal” side canceling. There is
magnetic dipoles
Electrons spinning about their axes no actual loss of generality in using the shape.
When a magnetic field is applied, a net alignment of these
magnetic dipoles occurs, and the medium becomes
magnetically polarized, or magnetized.
The magnetic polarization M, unlike electrical polarization
P, might be parallel to B (paramagnets) or opposite to B
(diamagnets).
Let’s calculate the torque on a rectangular current loop in
A few substances (ferromagnets) retain their magnetization a uniform magnetic field.
even after the external field has been removed. 1 2
5 6
6.1.3. Effect of a Magnetic Field on Atomic Orbits Electron Speeds Up or Slows Down
1 e2 v2
Electrons not only spin; they also revolve around the nucleus. = me without the magnetic field.
4πε 0 R 2 R
Let’s assume the orbit is a circle of radius R. The current
looks like steady (really?) The centripetal force comes from two sources:
the electric force and the magnetic force.
e e ev
Current I = = =
T 2πR 2πR The negative charge
v of the electron 1 e2 v′ 2
+ ev ′B = m
4πε 0 R′2 R′
e
evR
Orbital dipole moment m = − zˆ (m = IπR 2 )
2 Assume R′ ≅ R
m (v′2 − v 2 ) me (v′ + v)(v′ − v)
ev′B = e =
R R
eRB
∴ ∆v = (v′ − v) ≈ When B is turn on, the electron speeds up.
7 2me 8
The Dipole Moment and The Diamagnetism 6.1.4 Magnetization
A change in the orbital speed means a change in the In the present of a magnetic field, matter becomes magnetized.
dipole moment. Upon microscopic examination, it contains many tiny dipoles,
1 e2 R 2
∆m = − e∆vR = − B with a net alignment along some direction.
2 2me
Two mechanisms account for this magnetic polarization:
The change in m is opposite to the direction of B.
1.Paramagnetism: the dipoles associated with the spins of
In the presence of a magnetic field, each atom picks up a unpaired electrons experience a torque tending to line them up
little “extra” dipole moment, and the increments are all parallel to the field.
antiparallel to the field. This is the mechanism responsible
2. Diamagnetism: the orbital speed of the electrons is altered
for diamagnetism.
in such a way as to change the orbital dipole moment in the
This is a universal phenomenon, affecting all atoms, but it is direction opposite to the field.
typically much weaker than paramagnetism.
We describe the state of magnetic polarization by the vector
quantity:
9 M ≡ magnetic dipole moment per unit volume. 10
6.2 The Field of a Magnetized Object Vector potential and Bound Currents
6.2.1 Bound Currents
Can the equation be expressed in a more illuminating form,
Suppose we have a piece of magnetized as in the electrical case? Yes!
material (i.e. M is given). What field does By exploiting the identity, ∂ ∂ ∂ 1
this object produce? (xˆ ′ + yˆ ′ + zˆ ′ )
1 rˆ ∂x′ ∂y′ ∂z′ ( x − x′) 2 + ( y − y′) 2 + ( z − z′) 2
The vector potential of a single dipole m is ∇′ = xˆ ′( x − x′) + yˆ ′( y − y′) + zˆ ′( z − z ′) rˆ
r r2 = =
(( x − x′) 2 + ( y − y′) 2 + ( z − z ′) 2 )3/ 2 r 2
µ0 m × rˆ µ0
A(r ) = 1
4π r 2 The vector potential is A(r ) =
4π ∫ M (r′) × (∇′ )dτ ′
r
In the magnetized object, each volume element carries a Using the product rule ∇ × ( fA) = ∇f × A + f (∇ × A)
dipole moment Mdτ’, so the total vector potential is and integrating by part, we have
µ0 1 M (r′)
4π ∫ r
µ M (r′) × rˆ A(r ) = [∇′ × M (r′)]dτ ′ − ∫ ∇′ × [ ]dτ ′
A(r ) = 0 ∫ dτ ′ r
4π r2 µ0 1 µ 1
how? Prob. 1.60
4π ∫ r ∫ [M(r′) × nˆ ′]da′
= [∇′ × M (r′)]dτ ′ +
0
11 4π r 12
Vector potential and Bound Currents
µ0 1 µ 1
A(r ) = ∫
4π r
[∇′ × M (r′)]dτ ′ + 0 ∫ [M (r′) × nˆ ′]da′
4π r
J b = ∇′ × M (r′) K b = M (r′) × nˆ ′
volume current surface current
∫ (∇ ⋅ E)dτ = v∫ E ⋅ da µ0 Jb µ Kb
∫ ∫
Gauss's law A(r ) = dτ ′ + 0 da′
v S 4π v r 4π S r
(∇ ⋅ ( v × c))dτ = c ⋅ (∇ × v )dτ
∫ ∫v
Let E = v × c, v The electrical analogy
v∫ ( v × c) ⋅ da = −c ⋅ v∫ v × da volume charge density ρ b = −∇ ⋅ P
S S
4π r − µ0 R 3σω sinψ yˆ
1
u
The surface current density K (r′) = σ v′ =
2 ∫
−1 r + R − 2rRu
2 2
du
15 16
− µ0 R 3σω sinψyˆ ( R 2 + r 2 + Rr ) | R − r | −( R 2 + r 2 − Rr )( R + r ) Example 6.1 Find the magnetic field of a
A(r ) = (− )
2 3R 2 r 2 uniformly magnetized sphere of radius R.
µ 0 Rσ
(ω × r ) inside
A(r ) = 24 Sol : Choosing the z axis along the direction of M,
µ 0 R 3σ (ω × r ) outside J ′b = ∇ × M = 0
2r we have
K ′b = M × nˆ ′ = M sin θφˆ
Reverting to the “natural” coordinate, we have The surface current density is analogous to that of a spinning
spherical shell with uniform surface current density.
19 20
6.3 The Auxiliary Field H
The Role of H in Magnetostatics
6.3.1 Ampere’s Law in Magnetized Materials
What is the difference between bound current and free H plays a role in magnetostatics analogous to D in the
current? J =J +J electrostotics.
b f
Ampere’s law can be written: D allows us to write Gauss’s law in terms of free change alone.
1
(∇ × B) = J = J f + J b = J f + ∇ × M D = ε 0 E + P, ∇ ⋅ D=ρ f
µ0
H permits us to express Ampere’s law in terms of free current
1
⇒ ∇×( B − M) = J f alone.
H=
1
B − M, ∇ × H = J f
µ0 µ0
H What we can control directly.
In terms of H, then the Ampere’s law reads
∇ × H = J f (differential form) Why can’t we turn the bound currents on or off
independently?
∫ H ⋅ dl = I f (integral form)
21 22
25 26
∇×H = J f ⇒ H above
//
− H below
//
= K f × nˆ
Problems: 6.4, 6.10, 6.13, 6.15
∇ ⋅ H = −∇ ⋅ M ⇒ H above
⊥ ⊥
− H below ⊥
= −( M above ⊥
− M below )
∫ H ⋅ dl = I f (integral form)
∫ M ⋅ da ≠ 0
HA = nAI ∴ H = nIzˆ Gaussian pillbox
⇒ ∇⋅M ≠ 0
B = µ 0 (1 + χ m )nIzˆ
The enhancement of the magnetic field strength depends on
the susceptibility of the material.
and J b = ∇ × M = ∇ × χ m H = χ m J f
Is there a material that the field is significantly enhanced? 31 32
6.4.2 Ferromagnetism Ferromagnetic Domains
37 38
39 40
Classical picture of the magnetization process
Properties of ferrite materials (II)
--- By treating the spinning electron as a gyroscopic top.
Ferrites are ceramiclike materials with specific resistivities
that may be as much as 1014 greater than that of metals and If an electron is located in a uniform static magnetic field
with dielectric constants around 10 to 15 or greater. B0, a torque is given by
m e
T = m × B0 = − p × B0 = B0 × P
p mo
Ferrites are made by sintering a mixture of metal oxides dP e
and have the general chemical composition MO·Fe2O3, =T= B 0 × P = ω0 × P
dt mo
where M is a divalent metal such as Mn, Mg, Fe, Zn, Ni, Cd,
etc. eB0
where ω0 = is called the Larmor fequency;
mo
=
Relative permeabilities of several thousand are common. P= is angular momentum; and
The magnetic properties of ferrites arise mainly from the 2
magnetic dipole moment associated with the electron e=
m= is magnetic dipole moment.
spin. 41 2mo 42
At a high enough temperature a zero net magnetization The unit of B is gauss; the unit of H is Oersted. They have same
results. This temperature is called the Curie temperature, Tc. dimension.
What does ∆H and Mr mean?
45 46
Ferrite linewidth and remanent magnetization
49 50
Forced precession of spinning electron (II) Real and imaginary permeability &
propagation and attenuation constant.
propagation constant
β ± = ω εµ ±
RHCP LHCP
ωm
M x+ = H + ,
ω0 − ω ωm ωm
+
⇒M = H+ ⇒ µ + = µ0 (1 + ) RHCP
− jω ω − ω ω0 − ω
M y+ = m
H + . 0
ω0 − ω ω
µ − = µ0 (1 + m ) LHCP
ω0 + ω 51 52
Faraday rotation --- a nonreciprocal effect Microwave gyrator
Consider linearly polarized
electric field at z=0, represented Gyrator with a
as the sum of a RHCP and a twist section.
LHCP wave:
E0 E
E ( z =0) = xˆE0 = ( xˆ − jyˆ ) + 0 ( xˆ + jyˆ )
2 2
These two polarized waves propagate with different
propagation constants.
E0 E
E ( z =A) = ( xˆ − jyˆ )e− jβ+A + 0 ( xˆ + jyˆ )e− jβ−A Gyrator without
2 2 a twist section.
β −β β −β
= E0 xˆ cos( + − )A − yˆ sin( + − )Ae− j ( β+ +β− )A / 2
2 2
E β −β
θ = tan−1 y = −( + − )A. This effect is called Faraday rotation.
Ex 2 53 54
55 56
Four-port circulator Three-port circulator
57 58