Ece329lecture Notes
Ece329lecture Notes
Fundamental building blocks of matter — electrons and protons at atomic produced without permission from the author.
1
which is Newton’s 2nd law of motion 1 for a particle under the influence of
Lorentz force
F = q(E + v × B). (2)
In view of (1), the operational definitions of fields E(r, t) and B(r, t) arise
from particle acceleration a = dv dt
that can be measured in the lab frame:
the electric field E is evidently force per unit stationary charge (i.e., v = 0)
whereas field B describes an additional force per charge in transport (i.e.,
qv) that acts in a direction perpendicular to v.
There are important differences between gravitational and electromag-
netic interactions: Gravitational interactions are always attractive indicating
that particle masses m that generate the gravitational field g(r, t) must all
have the same algebraic sign (taken to be positive by convention). Electro-
magnetic interactions, on the other hand, are attractive or repulsive depend-
ing on particle charges q which can be positive or negative. By convention
a positive charge q = e ≈ 1.6 × 10−19 C is attributed to the fundamental
particle know as proton, while, again by convention, q = −e for an electron,
the sole companion of the proton within a hydrogen atom2 . Protons and
electrons are charged elementary building blocks3 of all atoms (hydrogen as
1
Valid so long as |v| ≪ c where c is the speed of light in vacuum.
2
Hydrogen atom exists as a consequence of mutual attraction between proton and electron counterbal-
anced by quantum mechanical constraints on allowed energy states — the constraints include the influence
of short-lived virtual particle/anti-particle pairs interacting with the proton and electron in a sporadic
manner.
3
Atoms can also contain in their nuclei varying numbers of an uncharged particle known as the neutron
which is responsible for different isotopes of chemical elements (e.g., the hydrogen isotope known as
deuterium contains a neutron in addition to a proton and an electron). While neutrons have no net
2
well as atoms of heavier elements) that constitute the matter around us. In a
collection of fundamental particles the total mass is always a monotonically
increasing function of the number of particles. However, that is not the case
with total charge since individual particle charges can be positive or nega-
tive. In fact, the net charge density ρ(r, t) found in macroscopic amounts of
matter is typically close to zero as a result of having nearly equal numbers of
protons and electrons in ordinary matter composed of charge-neutral atoms
and molecules4 .
charge, they consist of charged sub-nuclear particles known as up ( 23 e) and down (− 31 e) quarks whose
motions within the neutron establish currents and a magnetic moment.
4
The reason why intrinsically weaker gravity becomes dominant in the macro world.
3
Fields are relative:
Physical laws that we use today to describe our surroundings have been
developed to have identical forms in all reference frames in uniform motion
with respect to one another. For instance, Lorentz force law on a charge q is
F = q(E + v × B) and F′ = q(E′ + v′ × B′ ) (3)
in terms of unprimed and primed variables measured in two reference frames.
Moreover, particle charge q and the speed of light c are assigned invariant5
values in reference frames in relative motion (thus q ′ and c′ are unnecessary
to invoke in physical models). The ramifications of these restrictions con-
stituting the special theory of relativity (first described by Einstein in
1905 and covered at UIUC in PHYS 325) are in full accord with experimen-
tal measurements. They are also well matched by Newtonian relations
(approximate but more intuitive laws of dynamics covered in PHYS 211) if
and when the relative speed of primed and unprimed frames is negligible
compared to the speed of light c.
Since in Newtonian descriptions mass m and acceleration dv dt have invari-
ant values in all reference frames, it follows that if and when |v′ − v| ≪ c,
then F′ = F, in which case (3) implies
E′ + v′ × B′ = E + v × B. (4)
5
Other “relativistic invariants”
p between different
p reference frames include particle (rest) masses and the
so-called “spacetime interval” t2 − L2 /c2 = t′2 − L′2 /c2 between two events ocurring at two locations
and two times separated by a distance L and time-delay t, respectively. Relativistic invariants are the
most prized physical quantitites to focus on in relativistic models (simply because they remain fixed in all
reference frames). Note that distances L 6= L′ and time-delays t 6= t′ are not relativistic invariants!
4
Then, for a stationary charge in the primed frame, we have v′ = 0 and
E′ = E + v × B, (5)
which indicates that force per unit stationary charge in the primed frame
— i.e., the electric field in the primed frame — is a linear combination of
electrical and magnetic forces exerted on the same charge as seen from an-
other reference frame (unprimed) where the charge appears to have a vector
velocity v.
Thus, electric and magnetic fields needed in the formulation of charged
particle interactions are not unrelated to one another — rather, they intermix
in a manner that depends on the reference frame6 being used for analysis
purposes. Note that charges q which are stationary in one reference frame
(and therefore carry no electrical current) will appear to be in motion in
another frame and thus carry electrical currents I. It must therefore be
evident that the equations for E and B in any reference frame must be cross-
coupled and depend on both charge and current densities that are measured
in the same frame.
6
Given E and B measured in the lab, E′ and B′ measured by an observer moving through the lab
with a constant velocity v are well approximated by E′ ≈ E + v × B and B′ ≈ B − v×E c2
so long as
8
|v| ≪ c = 3 × 10 m/s, the speed of light in free space (shown by relativistic analysis discussed in PHYS
225 — exact transformation formulae are E′k = Ek , B′k = Bk , E′⊥ = γ(E⊥ + v × B⊥ ), B′⊥ = γ(B⊥ − v×Ec2
⊥
),
1
where γ = q v2 ).
1−
c2
5
Maxwell’s field equations:
The required set of coupled equations governing E and B was “discovered” in
1864 by James Clerk Maxwell to be (first introduced in PHYS 212 in integral
form and discussed throughout this course)
ρ
∇·E= ǫo Divergence eqn’s ∇·B=0
∇ × E = − ∂B
∂t
Curl eqn’s ∇ × B = µoJ + µoǫo ∂E
∂t
where
H 1 1 F
µo ≡ 4π × 10−7 and ǫo = ≈
m µoc2 36π × 109 m
in mksA units and
1 m
c=√ ≈ 3 × 108
µ o ǫo s
is the speed of light in free space. Furthermore ρ = ρ(r, t) refers to the net
charge density and J = J(r, t) to the current density in the measurement
frame, whereas ∇ · E and ∇ × E refer to the divergence and curl of vector
field E generated by partial differentiation of the orthogonal components of
E (concepts introduced in MATH 241 and reviewed in Lecture 4).
6
Solutions of Maxwell’s equations — waves and static fields (AC/DC):
Maxwell’s partial differential equations shown above, describing the coupled
dynamics of electric and magnetic fields E and B in response to space and
time varying source fields ρ and J, require an extended study to appreciate
their full ramifications and predictions. All predictions of these equations
have been experimentally verified and it has been found out that everything
that is known and observed about electricity and magnetism can be explained
in terms of these equations and their quantized forms. t=0 c
One of their predictions, derived specifically in Lecture 18, is that they
support traveling wave solutions of the form (a)
L
z
E(r, t) ∝ B(r, t) ∝ cos(2πf (t − )) (6)
c
in regions where J = ρ = 0. These are co-sinusoidal field perturbations hav- t = L/2c
c
ing oscillation frequencies f , oscillation periods T = f1 , wavelengths λ = fc ,
(b)
and they travel in 3D space with the speed of light c in free space. Since L
perturbation will travel across a region of size L during a time interval L/c
(c)
as illustrated in the margin. L
7
Electrostatics Magnetostatics
ρ
∇·E= ǫo
Divergence eqn’s ∇·B=0
shown in the left and right columns above — these were obtained by sim-
ply setting the terms ∂E/∂t and ∂B/∂t in the curl equations to zero. In-
dependent “curl-free” static electric fields E(r) and “divergence-free” static
magnetic fields B(r) satisfying these simplified equations are naturally far
easier to determine than the coupled dynamic fields E(r, t) and B(r, t) to be
encountered in response to time-varying sources ρ(r, t) and J(r, t).
Quasi-static fields:
Even though in practical cases of interest (in physics and engineering) time-
varying sources are the “rule” and static sources an “exception”, learning to
solve the simplified set of electrostatics and magnetostatics equations turns
out to be invaluable. The reason is, static solutions often provide accurate
approximations — known as quasi-static approximation — for time-varying
field problems involving slowly-varying sources ρ(r, t) and J(r, t).
More specifically, if the source variation period T is much longer than the
travel time L/c of field perturbations across a region of size L, that is, if
L
T ≫ , (7)
c
8
then field calculations for the entire region can be done statically using the t=0 c
instantaneous (as opposed to retarded or previous) values of field sources
ρ and J. This is true because under the given condition source strengths (a)
L
will remain nearly constant over time intervals needed to communicate the
new fields to the most distant corners of the region of interest. We can also
re-state the same inequality (7) as t = L/2c
c
c
L ≪ cT = =λ (8) (b)
L
f
using the definition of wavelength λ introduced earlier. The indication is
then, any system with a physical size L that is small in terms of wavelength t = L/c
c
λ of the applied field variations can be analyzed quasi-statically by starting
from Maxwell’s static equations. (c)
L
9
computing7 . While the physical dimensions of electronic chips and micro-
circuits are generally very small, such elements can still be electrically large
in the sense that L ∼ λ because of reduced wavelengths λ at high clock
speeds f = 1/T . Thus, even the computer engineers (CompE’s) amongst us
need to understand and learn how to mitigate (and take advantage of) the
ramifications of Maxwell’s equations.
7
E.g., Taflove, “Why study electromagnetics”, IEEE APM, 44, 132, 2002.
10
Details and study plan:
So much for the big picture about fields and waves encountered in nature
and engineering systems and circuits. Working details of how fields and wave
effects can be computed and characterized will be provided in the remaining
parts of these notes.
Over the course of 39 lectures we will develop and study, in succession,
the equations and applications of electrostatics (Lectures 1-11), magnetism
(Lectures 12-15), and electromagnetics (Lectures 16-39) with a focus on time
varying (quasi-static as well as wave-like) phenomena.
ECE 329:
We start by finding out how the equations of electrostatics arise from the
familiar Coulomb’s law (like charges repel, unlike charges attract) and the
idea of field superpositions. We learn to solve electrostatic problems using the
notion of electrostatic potential (voltage) and develop the notions of polariza-
tion, conduction, charge continuity, and capacitance in quasi-static settings
of practical importance.
Next we learn how magnetic fields arise from charges in motion (a relative
concept depending on the reference frame of the observer) and develop the
governing laws of magnetostatics (also an extension of Coulomb’s law seen
from different reference frames). The vector potential is introduced for mag-
netic field calculations from prescribed current configurations, and notions
of magnetization and inductance are subsequently developed and applied in
11
quasi-static settings.
Just like time-varying electric fields imply time-varying charge densities
(or vice versa) in electro-quasi-statics (EQS), time-varying currents imply
time-varying magnetic flux in magneto-quasi-statics (MQS). We also learn
that time-varying magnetic-flux is accompanied by time varying electric fields
— a key finding of Faraday’s called induced field with paradigm shifting
ramifications and applications — and requires the modification of curl-free
electric field condition into a dynamic equation known as Faraday’s law.
Finally, the full set of Maxwell’s equations is reached after adding a
time-varying electric field term to the curl equation of magnetostatics. This
change, first introduced by Maxwell in order to make sure that the govern-
ing equations of electricity and magnetism are consistent with conservation
of charge, acknowledges the two-way coupling and feedback between electric
and magnetic fields: time-varying magnetic fields induce time-varying elec-
tric fields — Faraday effect — and time-varying electric fields in turn induce
time-varying magnetic fields (call it the “Maxwell effect”) in order to sustain
electromagnetic field variations in regions far away from charges and current
loops — that is the way nature seems to work (and here we are to observe
all that thanks to Maxwell effect allowing us to be here).
A study of wave solutions of Maxwell’s equations follows, including plane
TEM waves in free space, linear and circular polarized waves, waves in con-
ducting media, normal incidence of waves on planar interfaces of homoge-
neous regions, energy and momentum transfer, guided waves in two-wire
transmission-line (TL) systems, transient response on TL circuits, resonant
12
oscillations in TL cavities, sinusoidal steady-state analysis of TL’s and dis-
tributed circuits, Smith Chart applications, and finally losses in TL systems.
That is the full scope of the 39 lectures of ECE 329 — the course ends
with an intensive study of distributed circuit concepts based on transmission
lines, a study that complements the lumped circuit techniques examined and
mastered in earlier courses.
ECE 329 is only the first half of our first-pass study of the fields and waves
topics essential in electrical engineering education. Important topics such as
radiation and antennas (generation details of electromagnetic waves by time-
varying currents) and dispersion (frequency dependence of wave propagation
speeds in material media) are barely mentioned or not at all in ECE 329.
These constitute the main topics of the follow-on course, ECE 350.
ECE 350:
ECE 350 starts with the discussion of electromagnetic radiation theory and
transmission antennas, continues with propagation and wave guidance ef-
fects (including dispersion, phase and group velocities, Doppler shifts, oblique
incidence, evanescence and tunneling effects, guided modes in parallel-plate,
rectangular, and dielectric slab waveguides), treats cavity fluctuations (in-
cluding resonant modes, blackbody radiation in 3D cavities, thermal noise),
and concludes with a discussion of antenna reception (including effective area,
available power, link equations).
13
Beyond ECE 329 and 350:
Students having gone through ECE 329 and 350 will find themselves ready to
encounter higher level courses in our curriculum focusing on different applica-
tion areas and frequency regimes of the implications of Maxwell’s equations.
It is a life-long endeavor to master these relationships which have precipitated
the scientific upheavals of the 20th century (relativity and quantum mechan-
ics) and have remained intact and essential despite the upheavals unlike most
aspects of classical physics. Our high speed electronics and communication
networks and devices are intrinsically and fundamentally based on fields and
wave concepts. Progress and innovation in these areas will require a deep
understanding of fields and waves and how they interact with novel materials
and structures.
Learn the basics and then go and invent the next thing!
14
ECE 329 Lecture Notes — Summer 09/11/21, Erhan Kudeki Copyright ©2021 Reserved — no parts of this
1 Vector fields and Lorentz force produced without permission from the author.
15
and
B(r, t) = (Bx(r, t), By (r, t), Bz (r, t))
that we refer to as E and B for brevity (dependence on position
r and time t is implied ). Maxwell’s equations:
ρ
• Field vectors E and B and electric charge and current densities ρ and ∇·E =
ǫo
∇·B = 0
J — describing the distribution and motions of charge carriers — are ∂B
∇×E = −
related by (i.e., satisfy) a coupled set of linear constraints known as ∂t
∂E
∇×B = µo J + µo ǫ o .
Maxwell’s equations, shown in the margin. ∂t
such that
– Maxwell’s equations are expressed in terms of divergence and curl
F = q(E + v × B),
of field vectors — recall MATH 241 — or, equivalently, in terms
of closed surface and line integrals of the fields enclosing arbitrary with
tists including Gauss, Ampere, Faraday, and Maxwell. in mksA units, where
1 m
They remain intact and essential despite the scientific upheavals c= √
µo ǫ o
≈ 3 × 108
s
(paradigm shifts) of 20th century: relativity and quantum physics10 .
is the speed of light in free space.
16
Given the charge and current densities ρ and J, Maxwell’s equations can be
solved for the fields E and B.
17
F
• Charge carrier positions r, velocities ṙ, and accelerations r̈ = m , as well
as forces F, fields E and B, and current density J are all described, in
general, in terms of 3D vectors.
where
p q
– |r| ≡ x + y + z and |E| ≡ Ex2 + Ey2 + Ez2 etc., are vector
2 2 2
z
r = (x, y, z)
magnitudes,
= xx̂ + y ŷ + z ẑ
r E
– r̂ ≡ |r| and Ê ≡ |E| etc., are associated unit vectors, y
– with dot products ẑ
ŷ
◦ r̂ · r̂ = 1, Ê · Ê = 1, x̂ · x̂ = 1, etc., but x
x̂
◦ x̂ · ŷ = x̂ · ẑ = ŷ · ẑ = 0 UNIT VECTORS AND A POSITION
VECTOR IN RIGHT-HANDED
– and cross products CARETESIAN COORDINATES
◦ x̂ × ŷ = ẑ,
◦ ŷ × ẑ = x̂,
◦ ẑ × x̂ = ŷ,
adopting a right-handed convention (see the margin note in
the next page).
18
• Recall that Right handed con-
vention: cross product vec-
tor points in the direction indi-
– Dot product A · B is defined as |A| times |B| times the cosine of cated by the thumb of your right
angle θ between A and B. hand when you rotate your fin-
gers from vector A toward vector
◦ Thus dot product is zero when angle θ is 90◦, as in the case B through angle θ you decide to
use.
of x̂ and ŷ, etc.
– Cross product A × B is defined as a vector with a magnitude B = |B|b̂
|A| times |B| times the sine of angle θ between A and B and a
|B| sin θ
direction orthogonal to both A and B in a right-handed sense n̂
θ â
(see margin note) . |B| cos θ A = |A|â
◦ Thus cross product is zero when the vectors cross multiplied A · B = |A||B| cos θ
DOT PRODUCT:product of
are collinear (θ = 0◦) or anti-linear (θ = 180◦). projected vector lengths
19
Also, v · B = (5, 10, 0) · (0, 0, 2) = 5 · 0 + 10 · 0 + 0 · 2 = 0.
F = q(E + v × B) v1 = 0
x
to determine the fields E and B at the origin. F1 = 2x̂
Solution: Using the Lorentz force formula first with F = F1 and v =v1, we note that z
Wb v2 = 2ŷ
Bx = 3 , and By = Bz = 0. x
m2
Hence, vector F2 = 2x̂ − 6ẑ
Wb z
B = 3x̂ 2 .
m
y
v3 = 3ẑ F3 = 2x̂ + 9ŷ
x
• In your first homework you will be asked to do a sequence of vector
exercises, including problems on volume, surface, and line integrals of
vector or scalar functions of space (i.e., “fields”). These problems should
be worked out with the help of your PHYS 212 and/or MATH 241 texts
and notes.
21
as ECE 210 (lumped circuits and linear systems con-
cepts including time- and frequency-domain approaches
and phasors).
22
2 Static electric fields — Coulomb’s and Gauss’s
laws
Static electric fields E(r) are produced by static (non-time-varying) distri-
bution of charges and obey and the electrostatic laws shown in the margin
where ρ(r) denotes the net charge density in 3D volume. Over the next few Laws of
lectures we will find out how these laws emerge from Coulomb’s law. electrostatics:
At the most elementary level, each stationary point charge (electron or
proton) Q is surrounded by its radially directed electrostatic field E given ∇ · E = ρ/ǫo
by Coulomb’s law, and in the presence of multiple charges the field vectors ∇×E=0
of all the charges are added vectorially (linear superposition holds) to obtain
a superposition field E.
r̂
• Coulomb’s law specifies the electric field of a stationary charge Q at q
the origin as
Q r = |r|r̂
z
E(r) = r̂ Q
4πǫor2 y
x
1
as a function of position vector r = (x, y, z), where ǫo ≈ 36π×10 9 F/m Force exerted by Q on q:
is a scaling constant known as permittivity of free space, F = qE
p with electric field
r = |r| = x2 + y 2 + z 2 Q
E= r̂
r 4πǫo|r|2
is radial distance from the charge, and r̂ = r radial unit vector pointing With multiple Q’s superpose
away from the charge. multiple E’s
1
stationary “test charge” q brought within distance r of Q (see
figure in the margin). If qQ > 0, force F is
repulsive (directed along
The existence of a Coulomb field accompanying each charge carrier in its rest r̂), if qQ < 0 it is at-
frame1 is taken to be a fundamental property of charge carriers (established tractive — like charges
by measurements). repel, unlike charges at-
tract.
• When multiple static charges Qn are present in a region, the force
on a stationary test charge q can be described as qE in terms of a
z
superposition field r − rn
X Qn q Qn
E= r̂
2 n rn
n
4πǫ o rn
r
2
Example 1: Charges Q1 = 4πǫo and Q2 = −2Q1 are located at coordinates r1 =
(1, 0, 0) = x̂ and r2 = (0, 1, 0) = ŷ, respectively. What is the expression for E(r)
and what is the explicit value of vector E(0)?
y
Q1(r − r1) Q2 (r − r2 )
E(r) = +
4πǫo |r − r1 |3 4πǫo |r − r2|3 1
E(r) 10
tors Ê(r) ≡ |E(r)| for the field E(r) obtained in Example 1 on a suit-
able grid established on xy-plane — such plots are useful or visualiza-
tion purposes. Note that arrows emanate out of the positive charge at 5
(x, y) = (1, 0) and converge upon the negative charge at (x, y) = (0, 1).
x
-10 -5 5 10
out the flux lines over the vector map shown in the margin!
-10
3
E, dS I Q
• According to Coulomb’s law, electrostatic field of a charge Q placed at S E · dS =
ǫo
the origin points out in the radial direction r̂ away from the origin and
S
has a magnitude Q
Q (a)
Er =
4πǫor2
E, dS I Q
E · dS =
that depends on radial distance r, but it does not depend on direction S′
ǫo
r̂. The product of Er with ǫo and the surface area of a sphere at radius (b) S′
Q
r, namely, S = 4πr2, yields I Q
E, dS S ′′ E · dS =
ǫo
ǫ o Er S = Q Q
(c) S ′′
independent of the radius of the sphere. Let’s re-write the same result
as I E, dS I
S ′′′ E · dS = 0
ǫo E · dS = Q, (d)
S S ′′′
Q
where Surface integral depends
only on the net amount of
charge contained within
H the surface --- charges
outside the surface don’t
– the “closed surface integral” S E · dS is called the flux of E over matter; surface shape doesn’t
matter; also charge motion
surface S bounding the volume V = 4π 3
r3 , within the surface does not
matter.
– which in turn denotes the limiting value of the sum of dot products
Ej · ∆Sj computed over all surface elements of S having incremen-
tal areas |∆Sj | and unit vectors ∆Sj /|∆Sj | pointing away from
volume V — the limiting value is obtained as all |∆Sj | approach
zero (i.e., with increasingly finer subdivision of S into |∆Sj | ele-
ments).
4
H
Although we obtained the equality ǫo S E·dS = Q above only for a spherical
surface S centered about charge Q, we can easily convince ourselves — see the
sketches on the right — that the equality should hold even when we distort
the shape of surface S and/or displace Q away from the center so long as
we do not move Q outside of S. All such variations are permitted because
of inverse r-square dependence of the Coulomb’s law and additive nature of
fields, and if Q is moved outside the surface then the surface integral (flux)
simply goes to zero.
• Hence, given an arbitrary shaped volume V enclosed by an arbitrary
shaped surface S and including a net electrical charge QV , and defining
a displacement field Displacement
D ≡ ǫoE, FV = C
D = ǫoE[=] m m m2
we obtain I
D · dS = QV , Gauss’s law
S
a constraint known as Gauss’s law. At this stage, the introduction of Gauss’s law
D is simply a notational convenience.
Gauss’s law offers an alternative to implementing an explicit sum of Coulomb
fields for calculating static field distributions E or D = ǫoE — the alternative
method can be used when charge distributions have simplifying symmetry
properties as will be illustrated in the next set of examples.
Also, later on we will learn that Gauss’s law is valid even when charges
QV within volume V are non-static (i.e., in motion), a condition under which
Coulomb’s law is no longer valid.
5
z
Q
∆z =
Example 2: Charged particles Q are located uniformly along the z-axis with an λ
average line density of λ C/m extending from z = −∞ to +∞. We will compute
y
the electrostatic field E of this charge distribution at a distance r from z-axis.
r = r(x̂ cos φ + ŷ sin φ)
Having an average charge density of λ C/m implies that individual charges Q φ
x
are spaced from one another by a distance ∆z = Qλ along the z-axis. Assuming
that charge locations are z = n∆z, where n is any integer, and using Coulomb’s Q point charge in C
law, we find that
∞ ∞
X Q r − ẑn∆z X λ∆z(r − ẑn∆z) λ charge density in C/m
E(r) = 2 |r − ẑn∆z|
= 3
,
n=−∞
4πǫo |r − ẑn∆z| n=−∞
4πǫo |r − ẑn∆z|
6
• The result
λ
E = r̂
2πǫor
obtained above, valid for r ≫ ∆z, and labelled as macroscopic
field , also represents at any r (and z) the space average of the mi-
croscopic field taken over small volumes having dimensions of many
∆z’s (inter-particle separations).
7
z
λ
Solution using Gauss’s law: We first notice that macroscopic electric field of a S r
charge distribution along the z-axis having an average charge density of λ C/m
λ
Er =
should be pointing in radial direction r̂ away from the z-axis (why?). L 2πǫor
y
Also its magnitude Er should be independent of azimuth angle φ by symmetry.
φ
x
As a consequence, we can apply Gauss’s law
I
D · dS = QV
S
as
ǫo Er 2πrL = λL
over the surface S of a cylindrical volume V of some length L and radius r
centered about the z-axis as shown in the margin — notice that our “clever”
choice of surface S in this problem resulted in the evaluation of the flux integral
in Gauss’s law without doing any calculus.
8
3 Gauss’s law and static charge densities
z
S
We continue with examples illustrating the use of Gauss’s law in macroscopic
y
field calculations: A ρs
Ex(x) =
2ǫo
Example 1: Point charges Q are distributed over x = 0 plane with an average surface ρS x
charge density of ρs C/m2 . Determine the macroscopic electric field E of this
charge distribution using Gauss’s law.
Ex(x)
Solution: First, invoking Coulomb’s law, we convince ourselves that the field produced
ρs
by surface charge density ρs C/m2 on x = 0 plane will be of the form E = x̂Ex (x) sgn(x)
where Ex (x) is an odd function of x because y- and z-components of the field will 2ǫo
cancel out due to the symmetry of the charge distribution. In that case we can x
apply Gauss’s law over a cylindrical integration surface S having circular caps of
area A parallel to x = 0, and obtain
I
D · dS = QV ⇒ ǫo Ex (x)A − ǫo Ex (−x)A = Aρs ,
S
which leads, with Ex (−x) = −Ex (x), to
ρs
Ex (x) = for x > 0.
2ǫo
Hence, in vector form
ρs
E = x̂ sgn(x),
2ǫo
where sgn(x) is the signum function, equal to ±1 for x ≷ 0.
1
z
y
Solution: Symmetry arguments based on Coulomb’s law once again indicates that we W W
expect a solution of the form E = x̂Ex (x) where Ex (x) is an odd function of x. −
2 2
Ex(x) ρW
In that case, applying Gauss’s law with a cylindrical surface S having circular caps 2ǫo
of area A parallel to x = 0 extending between −x and x < W2 , we obtain
I W W x
−
D · dS = QV ⇒ ǫo Ex (x)A − ǫo Ex (−x)A = ρ2xA, 2 2
S
2
Note that the field solution depicted in the margin in terms of Ex (x) plot is a con-
tinuous function of x as opposed to the discontinuous Ex (x) solution obtained in
Example 1 for the macroscopic field of a surface charge.
– Coulomb field
Q
E = r̂ of a point charge Q,
4πǫor2
– Field
λ
E = r̂ of constant line density λ,
2πǫor
– Field
ρs
E = x̂ sgn(x) of constant surface density ρs ,
2ǫo
– Field
ρx
E = x̂ of constant volume density ρ
ǫo
as building blocks — that is, the above field equations can be super-
posed to determine the field structure of charge distributions ρ(x, y, z)
that can be expressed as superpositions of simpler charge distributions
with known field structures. Some examples...
3
Example 3: Consider a pair of surface charges ρs > 0 and −ρs C/m2 of equal mag-
nitudes placed on x = − W2 and x = W2 surfaces. Determine the electric field of
this charge distribution depicted in the margin.
E = E+ + E− = = x̂ rect( )
0, otherwise, ǫo W
x
as depicted in the margin.
ρs > 0 −ρs
Note that the field lines of our solution point from positive charges on one surface to
the negative charges resting on the other surface — this field has the structure Ex(x)
ρs
of fields encountered in parallel plate capacitors that we will be studying soon. ǫo
W W x
−
2 2
4
z
−ρ1 < 0 ρ2 > 0
Example 4: An infinite charged slab of width W1, located over −W1 < x < 0, has
E
a negative volumetric charge density of −ρ1 C/m3 , ρ1 > 0. A second slab of - +
width W2 and positive charge density ρ2 is located over 0 < x < W2 as shown x
in the margin. Compute the electric field of this static charge configuration if
W1ρ1 = W2ρ2 , implying that the entire system is charge neutral (i.e., a net charge
−W1 W2
of zero).
E1x(x)
Solution: We note that the field of slab W1 can be written as ρ1W1
2ǫo
ρ1 W1
x̂ 2ǫo , for x < −W1
W1 −W1 x
ρ (x+ 2 )
E1 = −x̂ 1 ǫo , for − W1 < x < 0
−x̂ ρ1 W1 ,
for x > 0
2ǫo
E2x(x)
as depicted in the margin. Likewise, the field of slab W2 is
ρ2W2
ρ2 W2
−x̂ 2ǫoW, for x < 0
2ǫo
2
E2 = x̂ ρ2 (x−ǫo
2 )
, for 0 < x < W2 W2 x
ρ W
x̂ 2 2 ,
for x > W2.
2ǫo
5
Gauss’ Law in terms of
• Charge density formalism which we find convenient to use for macro- charge density:
scopic field calculations can also be “adjusted” to describe the distri-
butions of isolated point charges via the use of impulses or delta I Z
functions in space. D · dS = ρdV
S V
– For example
ρ(x, y, z) = Qδ(x − xo)δ(y − yo )δ(z − zo)
can be regarded as a 3D volumetric charge density function rep- ρ(x, y, z) = Qδ(x)δ(y)δ(z)
resenting a point charge Q located at a coordinate z
r = (x, y, z) = (xo, yo, zo) ≡ ro.
y
◦ This is justified because we can regard δ(x − xo) to be zero Q
x
everywhere except at x = xo. By extension, the product
δ(x − xo)δ(y − yo )δ(z − zo)
3D impulse here
is zero everywhere except at r = ro = (xo, yo, zo ) — therefore where point charge
Q is localized over
the density function ρ(x, y, z) defined above behaves correctly a region of zero
volume
to indicate the absence of charges everywhere with the ex-
ception of ro. Furthermore, the area property of the impulse
implies that the volume integral of the charge density yields
Z Z Z Z
ρdV = Qδ(x − xo)δ(y − yo )δ(z − zo)dxdydz = Q
as it should.
6
◦ Notice that the shifted impulses δ(x − xo), etc., must have ρ(x, y, z) = ρs (y, z)δ(x − xo)
m−1 units in order to maintain dimensional consistency in the
above expression. ρs > 0
– Another example is z
ρ(x, y, z) = ρs(y, z)δ(x − xo) y
representing a surface charge density of ρs(y, z) C/m2 on x = xo
plane.
xo x
Example 5: Figure in the margin depicts (for the d = 1) the Ê-field of a pair of
charges ±Q located at (0, 0, ± d2 ) derived from
1.0
Q(r − d2 ẑ) −Q(r + d2 ẑ)
E(r) = +
4πǫo |r − d2 ẑ|3 4πǫo |r + d2 ẑ|3
0.5
Q (x, y, z − d2 ) (x, y, z + d2 )
= [ − ] V/m.
4πǫo |(x, y, z − d2 )|3 |(x, y, z + d2 )|3
z
0.0
R
Determine the electric flux xy E · dS across the entire xy-plane using dS =
−ẑdxdy. -0.5
Solution: Because of linearity, the flux we want to calculate equals the sum of the flux -1.0
due to charge Q at (0, 0, d2 ) above xy-plane and the flux due to charge −Q at -1.0 -0.5 0.0
x
0.5 1.0
7
E · dS = ǫQo for any S surrounding Q, we can, by symmetry,
H
Since by Gauss’s law S
infer that
Q
Z
E · (−ẑdxdy) =
xy 2ǫo
when only charge Q is considered — the logic here is, half of flux S E · dS = ǫQo
H
emanating from charge Q should go up and the remaining half H should go−Q
down
crossing the xy-plane in downward direction. Likewise, since S E · dS = ǫo for
any S surrounding −Q, again by symmetry, we can infer
Q
Z
E · (−ẑdxdy) =
xy 2ǫo
due to charge −Q only — the logic in this case is, half of flux ǫQo “entering” charge
−Q is “coming from” above crossing the xy-plane in downward direction.
8
4 Divergence and curl
Expressing the total charge QV contained in a volume V as a 3D volume
integral of charge density ρ(r), we can write Gauss’s law examined during
the last few lectures in the general form
I Z
D · dS = ρdV.
S V
This equation asserts that the flux of displacement D = ǫoE over any closed
surface S equals the net electrical charge contained in the enclosed volume
V — only the charges included within V affect the flux of D over surface
S, with charges
H outside surface S making no net contribution to the surface
integral S D · dS.
z
• Gauss’s law stated above holds true everywhere in space over all sur- (x, y, z + ∆z)
faces S and their enclosed volumes V , large and small. y (x, y + ∆y, z)
(x, y, z) (x + ∆x, y, z)
• Application of Gauss’s law to a small volume ∆V = ∆x∆y∆z sur-
rounded by a cubic surface ∆S of six faces, leads, in the limit of van-
ishing ∆x, ∆y, and ∆z, to the differential form of Gauss’s law expressed
x
in terms of a divergence operation to be reviewed next:
1
– Again under the same assumption
I
D·dS ≈ (Dx|2−Dx|1)∆y∆z+(Dy|4−Dy|3)∆x∆z+(Dz|6−Dz|5)∆x∆y
S
x
∇ · D = ρ,
D = (Dx , Dy , Dz )
2
∂
following the usual dot product rules, except that the product of ∂x and Dx,
for instance, is treated as a partial derivative ∂D
∂x . In the left side above
x
∂x ∂y ∂z
∂ ∂ ∂ 2
C
= 5 + 0 + 0 = 5 3.
m
z
0
Note that the divergence of vector D is a scalar quantity which is the volumetric
-1
-3
-3 -2 -1 0 1 2 3
x
3
z
−ρ1 < 0 ρ2 > 0
Ex(x)
W2
• Summarizing the results so far, Gauss’s law can be expressed in integral
−W1
x
as well as differential forms given by
I Z
ρ1W1
D · dS = ρdV ⇔ ∇ · D = ρ. −
ǫo
S V
4
over volume V on both sides)
I Z
D · dS = ∇ · D dV
S V
Given their curl and divergences, vector fields can be uniquely recon-
structed in regions V of 3D space provided they are known at the
bounding surface S of region V , however large (even infinite) S and V
may be — this is known as Helmholtz theorem (proof outlined in
Lecture 7).
5
operator ∇, like a cross product
x̂ ŷ ẑ
∂ ∂ ∂ ∂ ∂ ∂
∇ × E ≡ ( , , ) × (Ex , Ey , Ez ) = ∂x ∂y ∂z (curl of E)
∂x ∂y ∂z
Ex Ey Ez
∂Ez ∂Ey ∂Ez ∂Ex ∂Ey ∂Ex
= x̂( − ) − ŷ( − ) + ẑ( − ).
∂y ∂z ∂x ∂z ∂x ∂y
E = x̂ cos y + ŷ1
x̂ ŷ ẑ
∂ ∂ ∂
∇×E = ∂x ∂y ∂z
2
cosy 1 0
∂ ∂ ∂ ∂ ∂ ∂ 0
y
= x̂( 0 − 1) − ŷ( 0 − cos y) + ẑ( 1 − cos y)
∂y ∂z ∂x ∂z ∂x ∂y
= x̂0 − ŷ0 + ẑ(0 + sin y) = ẑ sin y -2
6
It is apparent that curl ∇ × E is stronger in those regions where E is rapidly
varying in directions orthogonal to the direction of E itself.
7
◦ The calculation is slightly more complicated than the following z
example (although similar in many ways) where we show that λ
the static electric field of an infinite line charge is curl-free. S r
λ
Er =
L 2πǫor
y
Example 4: Recall that the static field of a line charge λ distributed on the z-axis is φ
x
λ
E(x, y, z) = r̂ ,
2πǫo r
where
x y
r2 = x2 + y 2 and r̂ = x̂ cos φ + ŷ sin φ = ( , , 0).
r r
Show that field E satisfies the condition ∇ × E = 0.
8
5 Curl-free fields and electrostatic potential
• Mathematically, we can generate a curl-free vector field E(x, y, z) as
∂V ∂V ∂V
E = −( , , ),
∂x ∂y ∂z
by taking the gradient of any scalar function V (r) = V (x, y, z). The
gradient of V (x, y, z) is defined to be the vector
∂V ∂V ∂V
∇V ≡ ( , , ),
∂x ∂y ∂z
pointing in the direction of increasing V ; in abbreviated notation, curl-
free fields E can be indicated as
E = −∇V.
1
– The prescription E = −∇V , including the minus sign (optional,
but taken by convention in electrostatics), ensures that electro-
static field E points from regions of “high potential” to “low po-
tential” as illustrated in the next example. Electrostatic fields E
point from regions of
“high V ” to “low V ”
Example 1: Given an electrostatic potential
V (x, y, z) = x2 − 6y V
4
Note that this field is directed from regions of high potential to low potential. Also note
that electric field vectors are perpendicular everywhere to “equipotential” contours. -4
-4 -2 0 2 4
2
in 3D space, since, as shown below, such integrals are “path independent” for
z p = (xp, yp, zp)
curl-free fields E = −∇V .
Ej
• The vector line integral ∆lj
y
Z o C
E · dl C′
p
o = (xo, yo, zo)
over an integration path C extending from a point p = (xp, yp, zp ) in
3D space to some other point o = (xo, yo, zo ) is defined to be x
– the limiting value of the sum of dot products Ej ·∆lj computed over
all sub-elements of path C having incremental lengths |∆lj | and
unit vectors ∆lj /|∆lj | directed from p towards o — the limiting
value is obtained as all |∆lj | approach zero (i.e., with increasingly
finer subdivision of C into |∆lj | elements).
3
will in general be path dependent except for when E is curl-free. Curl-free: path-independent
line integrals
Cu
1.0
Example 2: The field E = x̂y ± ŷx is curl-free with the + sign, but not with − as
0.8
verified below by computing ∇ × E. Calculate the line integral of E (for both
signs, ±) from a point o = (0, 0, 0) to point p = (1, 1, 0) for two different paths
0.6
C going through points u = (0, 1, 0) and l = (1, 0, 0), respectively (see margin). Cl
z
Solution: First we note that 0.4
x̂ ŷ ẑ 0.2
∂ ∂ ∂
∇ × (x̂y ± ŷx) = ∂x ∂y ∂z = ẑ(±1 − 1)
y ±x 0 0.0
0.0 0.2 0.4 0.6 0.8 1.0
x
which confirms that E = x̂y ± ŷx is curl-free with with + sign, but not with −. “Curly”: path-dependent line
In either case, the integral to be performed is integrals
Z p Z p Z p Cu
1.0
E · dl = (Exdx + Ey dy + Ez dz) = (y dx ± x dy).
o o o
0.8
z
o y=0 x=0 0.4
Clearly, the result shows that the line integral o E · dl is path independent for x
E = x̂y + ŷx which is curl-free, and path dependent for E = x̂y − ŷx in which
case ∇ × E 6= 0.
4
• The mathematical reason why curl-free fields have path-independent
line integrals is because in those occasions the integrals can be written
in terms of exact differentials:
∂V ∂V ∂V
dV ≡ dx+ dy+ dz, is precisely −Ex dx−Ey dy−Ez dz = −E·dl. Vo = 0
∂x ∂y ∂z
– In that case x
Z o Z o Z p
E · dl = − dV = dV = Vp − Vo
p p o
5
• The physical reason why this integral formula for potential Vp works z Vp =
Z
o
· dl
p E
with any integration path is the principle of energy conservation: E(r)
y
Ro dl
– integral p E · dl, namely the “voltage drop” from p to o, repre-
sents the work done per unit charge by the field E in moving Vo = 0
charges from location p to location o1, so if the line integral were
x
path-dependent (in reaching from p to o) there would be ways As long as E is curl-free, line
of creating net energy by making a charge q follow special closed integral is path-independent and
produces the voltage drop from
point p to "ground" o.
paths within the electrostatic field E, in violation of the general
principle of energy conservation (that permits energy conversion
but not creation or destruction).
1
Either to increase the kinetic energy of the charge if charge transport from p to o is unimpeded (as
for a test charge accelarating between a pair of capacitor plates) or else in pushing the charge against
frictional forces (as through a resistive wire) both at the expense of the energy stored in the field. On the
other hand, work done (i.e., the voltage drop) by the field would be negative if charges q > 0 were moved
from p to o against the local electric field (as within a battery), in which case there would be a positive
voltage rise from p to o representing energy gain for the field per unit charge transported from p to o.
6
z Z
p
Vp = − o E · dl
y
Example 3: Given that Vo = V (0, 0, 0) = 0 and Z
V
E = 2xx̂ + 3z ŷ + 3(y + 1)ẑ ,
m Y
o
determine the electrostatic potential Vp = V (X, Y, Z) at point p = (X, Y, Z) in X x
volts.
A voltmeter with its (+/red)
Solution: Assuming that the field is curl-free (it is), so that any integration path can probe contacting point p and
be used, we find that its (-/black) probe contact-
Z o Z p Z p ing point o would display (by
Vp = E · dl = − E · dl = − (2x dx + 3z dy + 3(y + 1) dz) analog or digital means) the
p o o numerical value of
Z X Z Y Z Z Z o
= − 2x dx|y,z=0 − 3z dy|x=X,z=0 − 3(y + 1) dz|x=X,y=Y Vp = E · dl
p
0 0 0
2
= −X − 0 − 3(Y + 1)Z. with the integration path
consisting of the path defined
This implies by the probe wires. The volt-
V (x, y, x) = −x2 − 3(y + 1)z V. meter reading would be inde-
pendent of the path config-
Note that uration when the field E is
electrostatic.
−∇(−x2 − 3(y + 1)z) = ∇(x2 + 3(y + 1)z)
For the voltmeter not to per-
= x̂2x + ŷ3z + ẑ3(y + 1)
turb the field it is probing,
its input impedance need to
yields the original field E, which is an indication that E is indeed curl-free. be much greater than the
impedance between points p
and o.
7
Alternate Solution — Exact Differential Method: Note that
Therefore
V (x, y, z) = −x2 − 3yz − 3z + C,
where the integration constant C should chosen so that V (0, 0, 0) = 0. The result
is
V (x, y, z) = −x2 − 3(y + 1)z
as before.
8
Example 5: According to Coulomb’s law electrostatic field of a proton with charge
Q = e (where −e is electronic charge) located at the origin is given as z
e
E= r̂,
4πǫo r2
where z=r
p (x, y, z)
r= x2 + y 2 + z 2 and r̂ = .
r
Determine the electrostatic potential field V established by charge Q = e with
the provision that V → 0 as r → ∞ (i.e., ground at infinity).
Solution: Field E and its potential V will exhibit spherical symmetry in this problem. e y
Therefore, with no loss of generality, we can calculate the line integral from a x
point p at a distance r from the origin to a point o at ∞ (the specified ground)
along, say, the z-axis. Approaching the problem that way, the potential drop
Or else (exact differential
from r to ∞ is method):
Z ∞
e e −e
V (r) = 2
ẑ · ẑdz E·dl = ( r̂)·r̂dr = d( ) = −dV
z=r 4πǫo z 4πǫo r 2 4πǫo r
e ∞ e
= − |r = . leading to
4πǫo z 4πǫo r
e
V (r) =
4πǫo r
(using an integration con-
• To convert electrostatic potential Vp (in volts) at any point p to poten- stant of zero).
tial energy of a charge q brought to the same point, it is sufficient to
multiply Vp with q (or just the sign of q, depending on which energy
units we want to use — see the next example).
9
Example 6: In view of Example 5, what are the potential energies of a proton e and
an electron −e placed at distance r = a away from the proton at the origin,
where distance
4πǫo ~2
a≡ 2 = 0.529 × 10−10 m
e me
stands for Bohr radius — it is the mean distance of the ground state electron in
a hydrogen atom from the center of the atom. Recall that e = 1.602 × 10−19 C
and ǫo ≈ 10−9/36π F/m.
For the proton, potential energy in Joules is calculated by multiplying V (a) = 27.2
V with q = e = 1.602 × 10−19 C. However, by referring to 1.602 × 10−19 J of
energy as 1 eV (electron-volt), it is more convenient to refer to potential energy
eV (a) of the proton at r = a as
Likewise, for a particle with charge q = −e, i.e., an electron, potential energy at the
same location is
−eV (a) = −27.2 eV.
10
6 Circulation and boundary conditions z
y
Since curl-free static electric fields have path-independent line integrals, it dS o=p
follows that over closed paths C (when points p and o coincide)
C
S
I
E · dl = 0, ∆lj
C Ej
x
H
where the C E · dl is called the circulation of field E over closed path C Closed loop integral over path
bounding a surface S (see margin). C enclosing surface S.
Note that the area increment
dS of surface S is taken by
convention to point in the
right-hand-rule direction
with respect to "circulation"
Example 1: Consider the static electric field variation direction C.
ρx
E(x, y, z) = x̂ z y
ǫo
that will be encountered within a uniformly charged slab of an infinite extent in
y and z directions and a finite width in xH direction centered about x = 0. Show (-3,4,0) (3,4,0)
Solution: Integration path C is shown in Hthe figure in the margin. With the help of
the figure we expand the circulation C E · dl as
Z 3 Z 4 Z −3 Z 0
ρx ρ3 ρx ρ(−3)
E = x̂ · x̂dx + x̂ · ŷdy + x̂ · x̂dx + x̂ · ŷdy
x=−3 ǫo y=0 ǫo x=3 ǫo y=4 ǫo
Z 3 Z −3
ρx ρx
= dx + 0 + dx + 0 = 0.
x=−3 ǫo x=3 ǫo
1
H
Note that in expanding C E · dl above for the given path C, we took dl as x̂dx
and ŷdy in turns (along horizontal and vertical edges of C, respectively) and
ordered the integration limits in x and y to traverse C in a counter-clockwise
direction as indicated in the diagram.
z
y
• Vector fields E having zero circulations over all closed paths C are dS
known as conservative fields (for obvious reasons having to do with C
S
their use in modeling static fields compatible with conservation theo- dl
rems). E
x
– The concepts of curl-free and conservative fields overlap, that is STOKE’S THM:
Circulation of E around close
I path C equals the flux over
enclosed surface S of the curl
E · dl = 0 ⇔ ∇ × E = 0 of E taken in direction of dS.
C dS points in right-hand-rule
direction with respect to
over all closed paths C and at each r. "circulation" direction C.
2
– the incremental area element dS on the right points across area S z
in the direction indicated by a right-hand rule as follows: y
dS
Point your right thumb in chosen circulation direction C; then your
right fingers point through surface S in the direction that should be C
S
adopted for dS. dl
H E
– Given Stoke’s theorem, C E · dl = 0 follows immediately for all x
STOKE’S THM:
C, if ∇ × E = 0 is true over all r. Circulation of E around close
path C equals the flux over
enclosed surface S of the curl
• Stoke’s theorem clearly implies that curl is circulation per unit of E taken in direction of dS.
area, just as the divergence theorem showed that divergence is flux dS points in right-hand-rule
direction with respect to
per unit volume. "circulation" direction C.
– The only difference is, curl also has a direction, which is the normal
C
unit of the plane that contains the maximal value of circulation z 4
2
per unit area found at that location (over all possible orientations (x, y, z + ∆z) dS = ∆y∆z x̂
of dS). y
1 (x, y + ∆y, z)
3
(x, y, z)
• We will verify Stoke’s thm after explaining the circulation per unit
area notion in steps:
3
C
– For a small rectangular contour “Cx ” on a constant x plane with z 4
2
sufficiently small ∆y and ∆z dimensions parallel to y and z axes (x, y, z + ∆z) dS = ∆y∆z x̂
(see figure in the margin), we have y
1 (x, y + ∆y, z)
I 3
(x, y, z)
E · dl ≈ Ez|2∆z − Ey|4∆y − Ez|1 ∆z + Ey|3∆y
Cx
= (Ez|2 − Ez|1 )∆z − (Ey|4 − Ey|3)∆y.
It follows that x
1 Ez|2 − Ez|1 Ey|4 − Ey|3
I
E · dl ≈ ( − )
∆y∆z Cx ∆y ∆z
and
1 ∂Ez ∂Ey
I
lim E · dl = ( − ) = x̂ · ∇ × E,
∆y,∆z→0 ∆y∆z Cx ∂y ∂z
meaning that the x component of ∇ × E is the circulation of E
per unit area on a constant x surface.
– Likewise, y and z components of ∇ × E are circulations of E per
unit area on constant y and z surfaces, and in general
1 1 1
I I I
∇×E = lim ( E·dl, E·dl, E·dl).
∆x,∆y,∆z→0 ∆y∆z C ∆x∆z Cy ∆x∆y Cz
x
4
magnitudes corresponding to the maximum values of circulations C
per unit area at every point. S
pertinent for a closed path C and its enclosed surface S, we will make Sum of circulations over small
squares cancel in the interior
use of the diagram shown in the margin. edges and only survive around the
exterior path C. This way,
circulation around C matches
the sum of the fluxes of curl E
– The circulations over small squares shown in the diagram are ap- calculated over the small squares.
proximately equal to the products of their areas and the normal Laws of
components of ∇ × E calculated at the center points (basd on electrostatics:
what we learned above).
– When allR such circulations covering surface S are added up, the
∇×E = 0
result is S ∇ × E · dS in the limit of vanishing size for the squares,
H ∇ · ǫo E = ρ
– as well as C E · dl because in the grand sum of the circulations
They also apply “quasi-statically”
over all the squares, all the contributions mutually cancel out (like over a region of dimension L
the overlapping edges of red and blue squares) except for those when a time-varying field source
ρ(r, t) has a time-constant τ much
calculated along the periphery C! longer than the propagation time
delay L/c of E(r, t) field varia-
tions across the region (c is the
We can now summarize the general constraints governing static electric fields speed of light).
as
In electro-quasistatics (EQS)
∇ × E(r) = 0, ∇ · D(r) = ρ(r), where D(r) = ǫoE(r). E(r, t) will be accompanied by
a slowly varying magnetic field
B(r, t) (to be studied starting in
Lecture 12).
5
• Vector fields E(r) and D(r) governed by these equations will in general
be continuous functions of position coordinates r = (x, y, z) except at
boundary surfaces where charge density function ρ(r) requires a repre-
sentation in terms of a surface charge density ρs (r).
ρ(r) = ρs δ(z)
1
would be
ρs ρs 0
z
E(r) = ẑ sgn(z) ⇒ D(r) = ẑ sgn(z).
2ǫo 2
-1
ρs ρs
E(r) = ẑ sgn(z)+Eo(r) and D(r) = ẑ sgn(z)+ǫoEo(r).
2ǫo 2
Since fields Eo(r) and Do(r) vary continuously, these field expressions
must satisfy
where
E+ ≡ E(x, y, 0+) and E− ≡ E(x, y, 0− )
6
D+
refer to limiting values of E at z = 0 plane from above and below,
ẑ D−
respectively, and likewise for z=0
– The equations can be further simplified as around the dotted path yields
Et+ = Et−
Dn+ − Dn− = ρs and Et+ = Et−
in w → 0 limit.
Having D = 0 in the region x < 0 requires that the field due to surface charge ρso
on x = 5 m plane must cancel the field due ρs = 2 C/m2 on x = 0 plane — this
requires that ρso be −2 C/m2 .
8
D = 3ŷ for x < 0.
ρs = 2C/m2 ρs = −6C/m2
z
Example 3: In the region x < 0 measurements indicate a constant displacement field
D = 3ŷ C/m2. Also, x = 0 and x = 5 m planes contain surface charge densities
of ρs = 2 C/m2 and ρs = −6 C/m2 respectively. Determine D for x > 0 if D is
known to be uniform in the intervals 0 < x < 5 m and x > 5 m. y
Solution: First we note that E = D ǫo
= ŷ ǫ3o V/m is tangential to x = 0 and x = 5 m
surfaces. Since the tangential component of E cannot change at any boundary, x
we will have a uniform Ey = ǫ3o in all regions, −∞ < x < ∞, implying that x=5m
Dy = 3 C/m2 throughout (caused by charges at |y| → ∞).
In summary,
ŷ3,
for x < 0,
D = x̂2 + ŷ3, C.
for 0 < x < 5 m m2
−x̂4 + ŷ3, for x > 5 m
9
7 Poisson’s and Laplace’s equations
Summarizing the properties of electrostatic fields we have learned so far, they
satisfy the laws of electrostatics shown in the margin and, in addition, Laws of
electrostatics:
E = −∇V as a consequence of ∇ × E = 0.
∇ · E = ρ/ǫo
• Using these relations, we can re-write Gauss’s law as
∇×E=0
ρ
∇ · E = −∇ · (∇V ) = ,
ǫo
from which it follows that
ρ
∇2V = − , (Poisson’s eqn)
ǫo
where
2 ∂ 2V ∂ 2V ∂ 2V
∇V ≡ + +
∂x2 ∂y 2 ∂z 2
is known as Laplacian of V . Poisson’s eqn:
ρ
– A special case of Poisson’s equation corresponding to having ∇2 V = −
ǫo
ρ(x, y, z) = 0
1
Focusing our attention first on Laplace’s equation, we note that the equation z V (d) = −3 V
densities (to be determined in Example 2), which are known to be responsible for
an electrostatic field E = ẑEz measured in between the plates. Each plate has V (z) = Az + B
some unique and constant electrostatic potential V since neither E(r) nor V (r)
can dependent the coordinates x or y given the geometry of the problem. V (z)
Using Laplace’s equation determine V (z) and E(z) between the plates if the potential
of the plate at z = 0 is 0 (the ground), while the potential of the plate at z = d
is −3 V.
Solution: Since the potential function V = V (z) between the plates is only dependent
on z, it follows that Laplace’s equation simplifies as
2 ∂ 2V ∂ 2V ∂ 2V ∂ 2V
∇V = + + = = 0.
∂x2 ∂y 2 ∂z 2 ∂z 2
This equation can be satisfied by
V (z) = Az + B
where A and B are constants to be determined. Now applying the given boundary
conditions, we first notice that (at the lower plate)
V (0) = (Az + B)|z=0 = B = 0.
2
z V (d) = −3 V
Applying the second boundary condition (at the top plate) we find
3V y
V (2) = (Az + 0)|z=2 = 2A = −3 V ⇒ A = − . V (0) = 0
2m V (z) =?
z=0
x
The upshot is, potential function
3
V (z) = − z, for 0 < z < 2 m.
2 z
Finally, we determine the electric field between the plates as z=d=2m
−3 3
3 ∂ 3 3V V (z) = − z
2
E = −∇V = −∇(− z) = ẑ ( z) = ẑ . 3 3
2 ∂z 2 2m E = −∇(− z) = ẑ
2 2
0
V (z)
Example 2: In Example 1 what are the surface charge densities of the metallic plates
located at z = 0 and z = 2 m surfaces?
3
Notice that our solution with equal and opposite charge densities on the parallel
surfaces implies that electrostatic fields are zero within the conducting plates
where the fields due to two charged surfaces are canceling out. This conclusion is
consistent with having constant electrostatic potentials within conducting regions
as will be discussed in the next lecture.
z
Example 3: A pair of copper blocks separated by a distance d = 3 m in x direction
Vo = 0
hold surface charge densities of ρs = ±2 C/m2 on surfaces facing one another as +
+
-
-
shown in the margin. The blocks are assigned constant potentials Vo = 0 and Vp + -
+ -
(see figure). What is the potential difference Vp ? + - Vp =?
− + +
E =0 + E -
+ -
Solution: Let D = x̂ǫo Ex denote the displacement vector in between the blocks, and + -
+
let D− = 0 denote the displacement vector within the block with a surface at +
-
- x
+
x = 0. Then the boundary condition equation used at x = 0 implies that 3 m
-
C 2
x̂ · (D+ − D− ) = ǫo Ex = 2 ⇒ Ex = .
m2 ǫo
In that case, potential difference between the blocks is
2 6
V = Ex d = 3= .
ǫo ǫo
Since the block on the left is at a higher potential (electric field vectors point
from high to low potential) assigned as Vo = 0, we must have
6
Vp = − .
ǫo
4
z
Poisson’s equation −ρ1 < 0 ρ2 > 0
ρ
∇2 V = −
ǫo
E
is used in regions where the charge density ρ(r) is non-zero. The following - +
x
example illustrates a possible use of Poisson’s equation.
−W1 W2
E1x(x)
Example 4: An infinite charged slab of width W1, located over −W1 < x < 0, has a ρ1W1
negative volumetric charge density of −ρ1 C/m3 , ρ1 > 0. A second slab of width 2ǫo
W2 and positive charge density ρ2 is located over 0 < x < W2 as shown in the
−W1 x
margin. The electric field of this static charge configuration under the constraint
W1ρ1 = W2ρ2 was computed in an earlier section as
(
−x̂ ρ1 (x+Wǫo
1)
, for − W1 < x < 0 E2x(x)
E=
x̂ ρ2 (x−W
ǫo
2)
, for 0 < x < W2 ρ2W2
2ǫo
and is depicted in the margin. Determine the electrostatic potential in the re- W2 x
gion and the potential difference V21 ≡ V (W2) − V (−W1) satisfying Poisson’s
equation.
Solution: This is a one dimensional geometry where E and potential V depend only on Ex(x)
−W1 W2
coordinate x. Therefore, Poisson’s equation ∇2V = −ρ/ǫo takes the simplified
form x
d2 V ρ(x)
= − .
dx2 ǫo
ρ1W1
−
Integral of this equation over x yields in the left dV
dx = −Ex , which implies, given ǫo
the electric field result from above,
(
ρ1 (x+W1 )
dV ǫo , for − W1 < x < 0
= ρ2 (x−W2 )
dx − ǫo , for 0 < x < W2
5
Integrating dV
dx
once more (i.e., finding suitable anti-derivatives with integration
constants), we find
ρ1 (x+W1 )2
(
2ǫo + V1 , for − W1 < x < 0
V (x) = 2
− ρ2 (x−W2ǫo
2)
+ V2 , for 0 < x < W2 z
−ρ1 < 0 ρ2 > 0
where the integration constants included on each line have been selected so that
V2 = V (W2), V1 = V (−W1). E
- +
Requiring a unique potential value at x = 0 (we can only associate a single potential x
energy level with each position in space) compatible with this expression for V (x),
we obtain −W1
ρ1 (0 + W1)2 ρ2 (0 − W2)2 Ex(x)
W2
+ V1 = − + V2 , −W1 W2
2ǫo 2ǫo
x
from which
ρ2 W22 + ρ1 W12 ρ2 W2 (W1 + W2 ) ρ1 W1(W1 + W2)
V21 = V2 − V1 = = = . ρ1W1
2ǫo 2ǫo 2ǫo −
ǫo
Note that the equation above can be solved for W1, W2, and W2 + W1 in terms of
V12 , ρ2 , and ρ1 , providing useful formulas for diode design (see ECE 440). We V (x)
V2
can also get useful specific formulae for V1 and V2 by imposing V (0) = 0, i.e.,
choosing x = 0 to be the reference point.
−W1
W2 x
V1
6
with an arbitrary ρ existing over a finite region in space can be obtained z ρ(r′ )
′
as r−r
ρ(r′)
Z
V (r) = ′
d3 r′ r r′
4πǫo|r − r |
where d3r′ ≡ dx′ dy ′ dz ′ and the 3D integral on the right over the primed
coordinates is performed over the entire region where the charge density y
is non-zero. O
x
– Verification: The solution above can be verified by combining a
number of results we have seen earlier on:
1. In Lecture 5 we learned that the electric potential V (r) of a
point charge e at the origin is
e
V (r) = .
4πǫo|r|
Clearly, this singular result is a solution of Poisson’s equa-
tion above (and the stated boundary condition) for a charge
density input of
ρ(r) = eδ(r).
2. Using ECE 210-like terminology and notation, the above re-
sult can be represented as
1
δ(r) → Poisson’s Eqn →
4πǫo|r|
7
identifying the output on the right as 3D “impulse response”
of the linear and shift-invariant (LSI) system represented
by Poisson’s equation.
3. Because of shift-invariance, we have
1
δ(r − r′) → Poisson’s Eqn → ,
4πǫo|r − r′|
meaning that a shifted impulse causes a shifted impulse re-
sponse.
The shifted impulse response is usually called “Green’s
function” G(r, r′) in EM theory.
4. Because of linearity, we are allowed to use superpositioning
arguments like
1
Z Z
ρ(r′ )δ(r−r′)d3r′ = ρ(r) → Poisson’s Eqn → ρ(r′ ) ′
d3r′ = V (r),
4πǫo|r − r |
which concludes our verification of the electrostatic1 potential
solution. Note how we made use of the sifting property of the
impulse (from ECE 210) in above calculation.
1
Also, in quasi-statics we use ρ(r′ , t) to obtain V (r, t) over regions small compared to λ = c/f , with f
the highest frequency in ρ(r′ , t).
8
• As an application of the general solution of Poisson’s equation, namely
ρ ρ(r′ )
Z
2 3 ′
∇ V =− ⇒ V (r) = d r,
ǫo 4πǫo|r − r′ |
we next provide an outline of the proof of Helmholtz theorem (see
p states that any vector field F(x, y, z) that vanishes in
Lecture 4) which
the limit r = x2 + y 2 + z 2 → ∞ can be reconstructed uniquely from
its divergence and curl:
F = −∇V + ∇ × A
9
– Likewise, the curl of F (with ∇ × ∇V = 0) leads us to, with a
divergence-free 3 A, to
∇′ × F(r′ ) 3 ′
Z
2 2
∇×F = ∇×∇×A = ∇(∇·A)−∇ A = −∇ A ⇒ A(r) = dr
4π|r − r′|
once again in analogy with Poisson’s equation4 .
3
To confirm ∇ · A = 0 directly, use the identity ∇ · (αG) = α∇ · G + G · ∇α to expand A(r) as
after also using integration by parts for an integrand that vanishes as |r| → ∞ and a symmetry relation
∇|r − r′ |−1 = −∇′ |r − r′ |−1 which is easy to confirm.
4
While the vector field A identified above is divergence-free, ∇×A in the F = −∇V +∇×A expansion
can also be replaced with ∇×A′ so long as A′ = A+∇Ψ since ∇×∇Ψ is unconditionally zero independent
of the choice of Ψ. Note it is possible to specify Ψ so that ∇ · A′ = ∇ · ∇Ψ = ∇2 Ψ 6= 0 in which case A′
will be a divergent solution of the ∇ × F equation above! The additive term ∇Ψ in A′ is analogous to
allowing a constant number to be added to V ! The freedom to specify Ψ and thus ∇ · A′ at will is known
as gauge freedom and any choice of Ψ making ∇ · A′ = 0 is known as Coulomb’s gauge.
10
−ρs
8 Conductors, dielectrics, and polarization (a)
- - - - - - - - - -
ρs
So far in this course we have examined static field configurations of charge Eo = ẑ
ǫo
distributions assumed to be fixed in free space in the absence of nearby
materials (solid, liquid, or gas) composed of neutral atoms and molecules. + + + + + + + + + + ρs
In the presence of material bodies composed of large number of charge-
neutral atoms (in fluid or solid states) static charge distributions giving rise (b)- - - - - - - - - - −ρs
to electrostatic fields can be typically1 found: ρs + + + + + + + + +
Eo
+
σ>0 E=0
1. On exterior surfaces of conductors in “steady-state”, - - - - - - - - - -
−ρs
Eo
2. In crystal lattices occupied by ionized atoms, as in depletion regions of + + + + + + + + + + ρs
semiconductor junctions in diodes and transistors.
A conducting slab inserted into a
region with field E_o (as shown
In this lecture we will examine these configurations and response of materials in b)develops surface charge which
cancels out E_o within the slab.
to applied electric fields. E_o relates to surface charge
as dictated by Gauss’s law and
superposition principle.
1
relation2
J = σE. (Ohm’s Law)
(a)
- - - - - - - - - - −ρs
• Simple physics-based models for σ will be discussed later in Lecture 11. ρs
Eo = ẑ
For now it is sufficient to note that: ǫo
within the conductor after applying an external electric field Eo (see E_o relates to surface charge
as dictated by Gauss’s law and
margin). superposition principle.
– The reason is, mobile free charges (e.g., electrons in metallic con-
ductors) within the conductor will be pulled or pushed by the
applied field Eo to pile up on exterior surfaces of the conductor
2
Linear behavior is possible provided charge carriers suffer occasional collisions within the medium.
3
PEC is an “idealization” that has no real counterpart, even though it is convenient to treat high
conductivity materials such as copper as PEC in certain approximate models and calculations. For “su-
perconducting materials” σ → ∞ only in the low frequency limit.
2
until a surface charge density ρs that is generated produces a sec- (a)
- - - - - - - - - - −ρs
ondary field −Eo that exactly cancels out the applied Eo within ρs
Eo = ẑ
the interior of the conductor. ǫo
• The transient “time-constant” τ for the decay of charge density ρ (and A conducting slab inserted into a
region with field E_o (as shown
hence E, as claimed above) in a homogeneous4 conductor (constant σ) in b)develops surface charge which
cancels out E_o within the slab.
can be obtained using the continuity equation E_o relates to surface charge
as dictated by Gauss’s law and
superposition principle.
∂ρ
+∇·J=0
∂t
representing the mathematical statement of charge conservation (de-
rived in Lecture 16). Using J = σE and ∇ · E = ρ/ǫo, we have
σ
∇ · J = σ∇ · E = ρ
ǫo
4
See Fisher and Varney, Am. J. Phys., 44, 464 (1976), for a discussion of contact potential between
different metals.
3
above, from which it follows that
∂ρ σ − ǫσo t
+ ρ = 0 with a damped solution ρ(t) = ρ(0)e .
∂t ǫo
The decay time-constant
ǫo
τ=
σ
is typically very short (∼ 10−18 s) in metallic conductors, which is why
such conductors are usually considered to be in steady-state (and have
zero interior fields).
5
Also applicable quasi-statically when externally applied fields Eo (t) change slowly with time-constants
much longer than ǫo /σ. The way conductors are treated in high frequency electromagnetic problems will
be described later on.
4
Dielectric materials and polarization:
• Dielectric materials consist of a large number of charge-neutral atoms - - - - - - - - - -
or molecules and ideally contain no mobile charge carriers (i.e., σ = 0). Eo
Dielectric P
• Electric fields produced by charges located outside or within a dielectric slab E = Eo −
ǫo
material will polarize the dielectric — meaning that its constituent Eo
atoms or molecules will be “stretched out” to expose their internal or + + + + + + + + + +
A dielectric slab inserted into a
“bound” charges, electrons and protons — which will in turn cause the region with an initial field E_o
will become polarized.
electric field inside the dielectric to become weaker than (but not zero,
Inside the polarized dielectric the
as in conductors) what the field would have been in the absence of field will be weaker than E_o, but
not reduced to zero as in a
polarization effect. conductor.
We will next examine this polarization process and see how Gauss’s law can
be re-stated to facilitate field calculations in dielectric materials containing
bound charge carriers, i.e., atomic/molecular electrons and protons which
are not free to drift away from one another indefinitely (neglecting possible
ionization events).
5
ρs + + + + + + + + + + ρs
−ẑ ≡ E1
• In the presence of an electric field E = ẑEz in the dielectric each neutral −ρs - -
d
- - - - - - - - ǫo
atom of the medium will be in a distorted (but not ripped apart) state ∆z 0 = E2
forming a ẑ oriented electric dipole, which can be visualized as a ρs + + + + + + + + + +
−ẑ
ρs
−ρs - ǫo
proton-electron pair with a small proton displacement d in z direction - - - - - - - - -
0
with respect to the electron. ρs + + + + + + + + + + ρs
−ẑ
−ρs - - - - - - - - - - ǫo
– Consider a regular array of such dipoles 0
ρs + + + + + + + + + + ρs
p ≡ edẑ, −ẑ
ǫo
−ρs - - - - - - - - - -
with ∆x, ∆y, and ∆z spacings between the dipoles (see margin), e
0
ρs =
so that the volumetric dipole density is ∆x∆y
1
Nd ≡ m−3,
∆x∆y∆z
within the array, and, furthermore,
e C
ρs =
∆x∆y m2
is the magnitude of charge density of the adjacent proton and
electron layers (see margin again) formed by arrays of adjacent
dipoles displaced in z by intervals ∆z.
– Assuming that the array is infinite in extent in x and y directions,
the proton and electron layers with surface charge densities ±ρs
will produce interior electric fields
ρs e/ǫo
E1 = −ẑ = −ẑ
ǫo ∆x∆y
6
ρs + + + + + + + + + + ρs
(pointing in opposite direction to E = ẑEz ), and exterior fields d −ẑ
ǫo
≡ E1
−ρs - - - - - - - - - -
∆z 0 = E2
E2 = 0 ρs + + + + + + + + + + ρs
−ẑ
−ρs - - - - - - - - - - ǫo
in between the dipole layers. Space averaged macroscopic electric
0
field within the array (with a spatial weighting proportional to ρs + + + + + + + + + + ρs
−ẑ
the size of regions with the fields E1 and E2) produced by the −ρs - - - - - - - - - - ǫo
7
• Let’s re-arrange the expression for E from above as
ǫ o E + P = ǫ o Eo
after multiplying it with ǫo and moving P to the left. Now, the term on
the right is ǫoEo = Do respresenting the displacement vector outside the
dielectric slab, and if we were to “adopt” the left hand side expression
as the “displacement vector” for the interior, i.e, take
D = ǫo E + P - - - - - - - - - -
Eo
in regions with non-zero P, then we would see that Dielectric P
slab E = Eo −
ǫo
– D = Do, i.e., the displacement is the same inside and outside the
Eo
slab, while electric fields E and Eo inside and outside differ by a + + + + + + + + + +
8
– outside dielectrics, Gauss’s law is ∇ · (ǫoE) = ρf , since ρb = 0 in
that case, and this can be expressed as ∇ · D = ρf , with D = ǫoE
as usual;
– using ∇ · D = ρf also within dielectrics, with D = ǫoE + P, leads
to ∇ · (ǫoE) = ρf − ∇ · P, which is the expected form for Gauss’s
law within dielectrics with ρb = −∇ · P;
– in conclusion, we can write Gauss’s law in a general form applicable - - - - - - - - - -
everywere as Eo
∇ · D = ρf , Dielectric P
slab E = Eo −
ǫo
with the understanding that D ≡ ǫoE + P and “includes” the
Eo
effect of bound charge density ρb = −∇ · P which may be non- + + + + + + + + + +
• The differential form of Gauss’s law that we will now write (without
the subscript f on ρf ) as
9
integral of the differential form) as
I Z
D · dS = ρdV,
S V
where the right side denotes the net free charge inside volume V .
P = ǫoχeE,
where
ǫ = ǫo(1 + χe) ≡ ǫr ǫo
is known as the permittivity of the dielectric, and
ǫ r = 1 + χe
10
◦ for glass 4 − 10,
◦ dry-to-wet earth 5 − 10, silicon 11 − 12, distilled water 81.
In certain materials χe and ǫ are found to be tensors — mean-
ing that P and D are no longer aligned with E. Such materials
are said to be anisotropic, but they will not be studied in this
course. Also, there is an exception to the condition χe ≥ 0 — in
collisionless plasmas χe < 0, as discussed in ECE 450.
11
9 Static fields in dielectric media
• Summarizing important results from last lecture:
D = ǫE = ǫoE + P,
P = D − ǫo E
if needed.
1
z
Example 1: A perfect dielectric slab having a finite thickness W in the x direction E = 18x̂ E = 3x̂ E = 18x̂
is surrounded by free space and has a constant electric field E = 18x̂ V/m in
its exterior. Induced polarization of bound charges inside dielectric reduces the
electric field strength inside the slab from 18x̂ V/m to E = 3x̂ V/m. What are ǫ = ǫo ǫ = ǫr ǫo ǫ = ǫo
the displacement field D and polarization P outside and inside the slab, and
what are the dielectric constant ǫr and electric susceptibility χe of the slab? x
2
• Our revised definition of displacement D = ǫE, where ǫ = ǫr ǫo, implies,
when combined with E = −∇V and ∇ · D = ρ, a revised form of
Poisson’s equation
ρ
∇2 V = − ,
ǫ
– provided that dielectric constant ǫr is independent of position so
that ∇ · D = ∇ · (ǫE) = ǫ∇ · E is a valid intermediate step in the
derivation of Poisson’s equation.
– Under the same condition Laplace’s equation ∇2V = 0 also re-
mains valid.
– Dielectrics where ǫr is independent of position are said to be ho-
mogeneous.
◦ In inhomogeneous dielectrics where ǫ varies with position
neither equation is valid, and one has to resort to the full
form of Gauss’s law in field and potential calculations.
In other words, don’t use Laplace’s/Poisson’s equations
in inhomogeneous media.
In the next example we have two homogeneous slabs side-by-side
making up an inhomogeneous configuration. In that case we can
use Laplace/Poisson within the slabs one at a time and then match
the results at the boundary using boundary condition equations
as shown.
3
ρs = 2ǫo V (2) =?
z
4
Based on above calculations of constants A and B, the potential solution for the
region is (
2z V, 0<z<1
V (z) =
2 + (z − 1) V, 1 < z < 2.
It follows that V (2) = 3 V.
Note that electric fields −2ẑ V/m and −ẑ V/m in the bottom and top layers point
from high to low potential regions. Electric field E is discontinuous at the bound-
ary at z = 1 m while displacement D is continuous — the continuity of normally
directed D is demanded by boundary condition equations in the absence of sur-
face charge. z
E=0 V =0
d
ρsd = ?
d1 ρs
Example 3: A pair of infinite conducting plates at z = 0 and z = d are grounded
and have equal potentials, say, V = 0. The region 0 < z < d is occupied by E=?
free space (i.e., ǫ = ǫo ) except that an infinite charge sheet with a static surface ρs0 = ?
charge density ρs is located at z = d1 < d. Determine (a) the electrostatic field 0
E(z) in regions 0 < z < d1 and d1 < z < d, and (b) the surface charge densities E=0 V =0
ρs0 and ρsd at z = 0 and z = d on conductor surfaces if d1 = d/2.
If ρs in Example 3 is a slowly-
Solution: (a) Laplace’s equation for the given geometry requires a linear (in z) poten- varying function of time, then
tial solution in regions 0 < z < d1 and d1 < z < d. Since electrostatic E = −∇V , slowly varying E, ρs0 , and ρsd cal-
we can therefore represent the electric field in these regions as culated with instantaneous values
( of ρs would constitute quasi-static
solutions which are valid so long
−ẑVo /d1, 0 < z < d1
E= as d ≪ c/f , with f the highest
+ẑVo /d2, d1 < z < d frequency in ρs (t).
5
where Vo ≡ V (d1) and d2 ≡ d − d1. Hence,
(
−ẑǫo Vo /d1 , 0 < z < d1
D = ǫo E = ,
+ẑǫo Vo /d2, d1 < z < d
6
z
V (2) =?
ρs = −2ǫo
2ǫo
7
To determine V (2), we integrate E = ẑ2(1 − 4z ) V/m from top to bottom plate
(grounded), obtaining
Z 0 Z 0
z
V (2) = E · dl = 2(1 − )dz
z=2 z=2 4
2
z 4 3
= 2(z − )|02 = −2(2 − ) = −2 · = −3 V.
8 8 2
8
10 Capacitance and conductance
z A A
Parallel-plate capacitor: Consider a pair of conducting plates with surface
areas A separated by some distance d in free space (see margin).
The plates are initially charge neutral, but then some amount of electrons
Q y −Q
are transferred from one plate to the other so that the plates acquire equal
and opposite charges Q and −Q, distributed with surface densities of ± Q A on
plate surfaces facing one another (as shown in the margin). x
√ d
• That way, in steady state and for d ≪ A, a field configuration con-
fined mainly to the region between the plates is acquired, satisfying the
condition that static field inside a conductor
√ should be zero. A weak
“fringing field” can be ignored if d ≪ A and thus the geometry well
approximates the case with infinite plates.
1
charged) right plate is
Z (d,0,0) z A A
d
Q d
Z
V = E · dl = dx = Q.
(0,0,0) x=0 ǫ o A ǫ o A
Q y−Q
The last result can be expressed as a linear charge-voltage relation
Q = CV x
with d
A
C ≡ ǫo
d
representing the capacitance of the parallel conducting plate ar-
dV
rangement that we call parallel plate capacitor. I = C
dt
• By differentiating Q = CV we obtain the charging rate of the capacitor
as + V (t) -
dQ dV
I= =C
dt dt
which is only possible, for ideal capacitors, if the capacitor plates are
externally connected to a circuit supplying a current as shown on the
right where the direction of I = dQ
dt
is in the direction of voltage drop V
across the capacitor, from the positively charged plate to the negatively
charged plate as shown.
2
given by Q −Q
dQ dV +
+
-
=C = I − GV, + E
-
-
dt dt +
+
-
-
J = σE
where G stands for the conductance of the capacitor (derived later +
+
-
-
x
in this lecture) and I the external current flowing into the non- +
+
-
-
ideal capacitor. +
+
-
-
– Therefore, for non-ideal capacitors the external current d
dV V (t)
I=C + GV, + -
dt I = GV (t) + C
dV
dt C
meaning that part of I goes into changing stored charge Q = CV
of capacitor plates and the rest to conduct a GV amount of leakage
1
current of the capacitor plates, and the equivalent circuit of the R=
G
non-ideal capacitor there contains a “shunt resistance” R = G1
accompanying C as shown in the margin.
A capacitor with a perfect dielectric between its plates will hold its charge
and stored energy indefinitely. However, if the dielectric is imperfect and
has a finite conductivity σ, charge will be transported from the positive to
negative plate by a volumetric current density
J = σE,
4
which will result in a quasi-static discharge of the capacitor and the loss of Q −Q
+ -
the stored energy W to Ohmic dissipation in the imperfect dielectric. +
+
-
E -
Just as capacitance C characterizes the energy and charge storage “ca- +
+
-
-
+ J = σE
pacity” of the capacitor, we can define a conductance G that relates the +
-
-
x
+ -
quasi-static discharge current I in between the plates of a capacitor to po- +
+
-
-
tential drop V : + -
d
5
resistance, but for a given σ, R increases with length d and de-
r
creases with increasing cross-sectional area A. Simple conductivity b a
models and J will be discussed next lecture.
z
Coaxial Cable: When we study guided wave propagation later in the course ℓ
we will learn about coaxial cables.
6
– It follows that
Q
Er = ,
2πǫℓr
and voltage drop from inner to outer conductor is
Z b Z b Z b
Q Q dr Q b
V = Er dr = dr = = ln .
r=a r=a 2πǫℓr 2πℓǫ r=a r 2πℓǫ a
Clearly, once again Q = CV , with
2π
C = b ℓǫ
ln a
representing the capacitance of the coax segment of length l.
7
z
Diode junctions: In Example 4 in Lecture 7 we derived the expression
−ρ1 < 0 ρ2 > 0
for potential drop V across a charged region of a total width of W1 + W2,
such that in region 1 where −W1 < x < 0 the charge density ρ = −ρ1 is
negative, while in region 2 where 0 < x < W2 the charge density ρ = ρ2 -
E +
is positive, with the additional constraint that the entire region is charge
x
neutral, meaning that ρ1W1 = ρ2W2.
By solving Poisson’s equation for this charge density configuration (see
margin) encountered in junction regions of semiconductor diodes (described −W1
Ex(x)
W2
in detail ECE 440) we had established that the voltage drop from x = W2 −W1 W2
to x = −W1 across the junction is given by x
Using the expressions above for junction voltage V and width W1 + W2,
−W1
we will next derive an expression for small signal capacitance of the diode
junction:
W2 x
V1
• In region 2 where x > 0, the junction holds a total positive charge of
Q = ρ2 W2A per cross-sectional area A.
Q
• Therefore, substituting A
for ρ2W2 in the expression for V above, and
8
√ z
also using the W1 + W2 ∝ V expression derived above, we obtain
−ρ1 < 0 ρ2 > 0
q
ρ1 +ρ2
ρ2W2(W1 + W2) Q 2ǫoV ρ1 ρ2
V = = , E
2ǫo 2ǫoA - +
2ǫoρ1 ρ2 √
r
Q=A V W2
ρ1 + ρ2 −W1
Ex(x)
−W1 W2
representing a non-linear charge-voltage relation (for a given charge x
profiles satisfying ρ1W1 = ρ2 W2).
dQ
The slope of any Q vs V curve is given by the derivative dV , whether or
dQ
not the curve is linear. The slope dV of a non-linear charge-voltage curve V (x)
V2
can be interpreted as a small signal capacitance C. For a diode junction,
differentiating the above equation, we find that −W1
dQ
r
ǫo ρ1ρ2 W2 x
C= =A . V1
dV 2V (ρ1 + ρ2 )
Small changes dV in junction voltage will accompany small changes dQ =
CdV in stored charge Q of the junction, but the amount CdV will itself
depend on V because C ∝ V −1/2.
9
11 Lorentz-Drude models for conductivity and (a)
In this lecture we will describe simple microscopic models for conductivity σ q>0
and electric susceptibility χe of material media composed of free and bound In the absence of an applied electric
field E, free charge q exhibits a
charge carriers. The models were first developed by Lorentz and Drude prior "random walk" between collisions such
that its average velocity v is zero.
to the establishment of quantum mechanics. In these models free charge By collisions we refer to the inreactions
of q with zero-mean miscroscopic electric
fields within the conductor due to charges
carriers motions are described using Newtonian dynamics and atoms are entrapped in the lattice.
1
carriers and −m vτ denotes a macroscopic friction force proportional to 1 m
no macroscopic current density J will be found in the absence of E in In 1 second all N charges in the
volume cross the shaded surface,
regions with homogeneous charge carrier densities. transporting Nq Coulombs of charge
per second per unit area.
mν m2 I= S J · dS
2
mass, collision frequency, and number density of qs, ms, νs, and
Ns, respectively, then current density can be expressed as Ohm’s Law
and
J = σE, DC
with conductivity
X Ns qs2
σ= σs and σs =
s
msνs
denoting the medium and species conductivities, respectively, un-
der DC conditions.
• With a time varying field E the corresponding current density will also
be time varying, in which case conductivity σ should be defined in the
frequency domain using phasor techniques (remember ECE 210).
J̃ = σ Ẽ,
3
with
X Nsqs2
σ= σs and σs = .
s
ms(νs + jω)
AC
– Note that the AC conductivity just derived can be approximated conductivity
by the DC conductivity derived earlier for all AC frequencies ω
much smaller than species collision frequencies νs .
◦ In many cases of practical interest, this condition can be easily
met, and we are often well justified to ignore the frequency
dependence and complex character of conductivity σ revealed
in above derivation.
– More advanced quantum mechanical derivations of σs give the
same results except with effective masses specified by quantum
theory replacing the particle masses ms used in classical models.
• Typical DC conductivities:
4
when the DC conductivity vanishes as a consequence of correlated
charge carrier motions which are ignored in the Lorentz-Drude
model.
– Possible values of ωo for a given atom can be obtained from the en-
ergy levels of bound states of the atom (calculated using standard
quantum2 models like in PHYS 214) and time constants τd = α1
(which are finite because energies ~ωo radiated away are also finite)
are related to observed line widths (2α) in the emission spectra of
excited atoms.
2
For a quantum mechanical derivation of susceptibility, see, e.g., Mott, “Elements of Wave Mechanics”,
Chapt 4, Sect 10 (1958); Miller, “Quantum Mechanics for Scientists and Engineers”, Sect 7.3 (2008).
5
Electron displacement having the inferred damped co-sinusoid form
−t/τd
cos( ωo2 − α2 t) ≈ roe−t/τd cos(ωo t)
p
r(t) = roe
+ E -
+
is “zero-input response” (remember ECE 210) of a linear second-order +
-
-
ODE that can be constructed using Newton’s second law of classical +
+ −e r e -
-
+
mechanics: +
-
+ −e r e -
-
d2 r + −e r e
– If we assume that mass m times acceleration dt2
of a displaced +
-
-
+
electron equals the sum of a -
we get
d2 r dr
m 2 = −eE − mωo2r − m2α ,
dt dt
for which r(t) given above is the zero-input solution in the absence
of E.
6
Consequently, dipole moment of a single polarized atom is
e2
p = −er = E,
mωo2 + E -
+ -
+
and polarization field in a dielectric with a dipole density of Nd is + r e
-
-
+ −e -
+
Nd e2 +
-
P = Nd p = E. + −e r e -
mωo2 + −e r e -
-
+ -
+
This result can also be written as -
P = ǫoχeE,
where DC
Nd e2/mǫo susceptibility
χe ≡ 2
ωo
is DC susceptibility. AC susceptibility can be derived using phasor
techniques, but at frequencies ω ≪ ωo , AC susceptibility is well ap-
proximated by the DC susceptibility derived above.
7
• Time variation of E will imply the time variation of electron displace-
ment r, so that there will be in effect a non-zero electron velocity
dr e dE
v= =−
dt mωo2 dt
capable of producing a current.
8
12 Magnetic force and fields and Ampere’s law
Pairs of wires carrying currents I running in the same (opposite) direction
are known to attract (repel) one another. In this lecture we will explain the
mechanism — the phenomenon is a relativistic1 consequence of electrostatic
charge interactions, but it is more commonly described in terms of magnetic
fields. This will be our introduction to magnetic field effects in this course. I
1
F
Brief summary of special relativity: Observations indicate that light (EM) waves can be
“counted” like particles and yet travel at one and the same speed c = 3 × 108 m/s in all reference
frames in relative motion. As first recognized by Albert Einstein, these facts preclude the possibility that
F
a particle velocity u could appear as I
u′ = u − v (Newtonian)
to an observer approaching the particle with a velocity v; instead, u must transform to the observer’s
frame as
u−v
u′ = , (relativistic) “Things should be made as simple
1 − uv
c2 as possible – but no simpler.”
so that if u = c, then u′ = c also. This “relativistic” velocity transformation in turn requires that positions — Albert Einstein
x and times t of physical events transform (between the frames) as
v
x′ = γ(x − vt) and t′ = γ(t − x), (relativistic)
c2
where γ ≡ √ 1
, rather than as
1−v2 /c2
x′ = x − vt and t′ = t, (Newtonian)
′
so that dxdt
= u and dx
dt′
= u′ are related by the relativistic formula for u′ given above.
Relativistic transformations imply a number of “counter-intuitive” effects ordinarily not noticed unless
|v| is very close to c. One of them is Lorentz contraction, implied by dx = dx′ /γ at a fixed t: since γ > 1,
dx < dx′ , and moving objects having velocities v appear shorter then they are when viewed from other
reference frames where v is determined. A second one is time dilation, implied by dt′ = dt/γ at a fixed x′ :
since γ > 1, dt′ < dt, and moving clocks having velocities v and fixed x′ run slower than clocks in other
reference frames where v is determined. Consider taking PHYS 325 to learn more about special relativity.
1
• Consider a current carrying stationary wire in the lab frame:
2
– In the electron frame the wire is found to have a positive charge
density λ′, and thus it has a radial electrostatic field
λ′ ′
E = r̂
2πǫor
implying an electrostatic force F′ = qE′ on a stationary charge q.
a positive density for non-zero |v| ≪ c. (e.g., articles in Am. J. Phys.: Webster, 29, 262, 1961; Matzek
and Russel, 36, 905, 1968; Arista and Lopez, 43, 525, 1975; Zapolsky, 56, 1137, 1988). ′ v 2 Iv
λ ≈ λ+ 2 = 2 = Ivµoǫo
c c
3
and force F′ = qE′ can be transformed back to the lab frame, (a) In the "electron frame" the wire
appears positively charged and
repelsa test charge q with
where q appears to be moving with velocity v, as (with no ap- force F’=qE’
While we assumed q to be stationary in the reference frame of the electrons Magnetic field lines close upon
themselves unlike electric field
in the above discussion (for the sake of simplicity), the results obtained above lines which start and stop on
point charges.
are found to be valid for all particle velocities v measured in the lab frame. Right hand rule: point your
right thumb in the direction
4
We also get the same result using the approximation F = F′ that can be justified when |v| ≪ c, which of current flow; your fingers
is typically true by a large margin for electron speeds in current carrying conducting metals — see HW. will point in direction φ̂.
4
Also, if there are multiple current filaments In , each generating its own field
Bn, force F on q can be calculated using a superposition method as with
electrostatic fields.
Magnetic field B of the infinite current filament I obtained above can
also be obtained by superposing the magnetic field increments
µoIdl × r̂
dB ≡ (Biot-Savart law)
4πr2
of directed current increments Idl, where r = rr̂ is a position vector extend- z
ing from the location of the current increment to the field position where dB
is being specified — this formula, known as Biot-Savart law, is only valid
when used in terms of infinitesimal segments Idl of time-invarying current I
loops.
µoI
r B= φ̂
• Magnetic field B of the infinite line current I “derived” above satisfies 2πr
a circulation relation I
B · dl = µoIC ,
C
with IC = I.
This integral for the circulation of static magnetic field B is found to be
valid (experimentally) for all closed circulation paths C, and is known
as Ampere’s law (for static magnetic fields). In Ampere’s law
– IC stands for the net sum of all filament currents In crossing any
surface S bounded by path C,
◦ flowing in the direction given by the “right-hand-rule”:
5
when the right thumb is pointed in the direction of dl along path
C, the direction of filament current In is specified as the direction
of the fingers of your right hand through surface S bounded by
contour C.
◦ Filament currents not crossing S — i.e., current filaments not
“linked” to path C — should not be included on the right hand
side of Ampere’s law.
where
I
6
• Stoke’s theorem re-stated for a vector field H as
I Z
H · dl = ∇ × H · dS
C S
7
– Likewise,
J = ẑ y rect(y − 0.5) δ(x)
J(x, y, z) = ẑI(z)δ(x − xo)δ(y − yo )
representats a line current I(z) measured in A units flowing in
z-direction along a filament defined by the intersections of x = xo
z y
and y = yo surfaces.
– As a most extreme case,
J(x, y, z, t) = Qvδ(x − xo)δ(y − yo )δ(z − zo)
represents the time-varying current density of a point charge Q xo x
at coordinates (x, y, z) = (xo(t), yo(t), zo(t)) moving with velocity
v = (ẋo(t), ẏo(t), żo(t)).
flowing on x = 0 plane (as shown in the margin). What is the total current I
flowing on the same plane measured in A units?
8
13 Current sheet, solenoid, vector potential and
current loops Js = ẑJs
x
Example 1: Consider a uniform surface current density Js = Js ẑ A/m flowing on µoJs
B(x) =
x = 0 plane (see figure in the margin) — the current sheet extends infinitely in 2
y and z directions. Determine B and H.
Solution: Since the current sheet extends infinitely in y and z directions we expect B As shown in Example 1 mag-
netic field of a current sheet
to depend only on coordinate x. Also, the field should be the superposition of the is independent of distance
fields of an infinite number of current filaments, which suggests, by right-hand- |x| from the current sheet.
rule, B = ŷB(x), where B(x) is an odd function of x. To determine B(x), such Also H changes discontinu-
that B(−x) = −B(x), we apply Ampere’s law by computing the circulation of ously across the current sheet
by an amount Js .
B around the rectangular path C shown in the figure in the margin. We expand
I
B · dl = µo IC
C
as
B(x)L + 0 − B(−x)L + 0 = µo Js L,
from which we obtain
µo Js µo Js Js
B(x) = ⇒ B = ŷ sgn(x) and H = ŷ sgn(x).
2 2 2
1
Example 2: Consider a slab of thickness W over − W2 < x < W2 which extends in-
finitely in y and z directions and conducts a uniform current density of J = ẑJo
A/m2 . Determine H if the current density is zero outside the slab.
Solution: Given the geometric similarities between this problem and Example 1, we
postulate that B = ŷB(x), where B(x) is an odd function of x, that is B(−x) =
−B(x). To determine B(x) we apply Ampere’s law by computing the circulation
of B around the rectangular path C shown in the figure in the margin. For
x < W2 , we expand I
B · dl = µo IC
C
as
z
B(x)L + 0 − B(−x)L + 0 = µo Jo 2xL ⇒ B(x) = µo Jo x.
y
W
For x ≥ 2 , the expansion gives
Joẑ
W
B(x)L + 0 − B(−x)L + 0 = µo Jo W L ⇒ B(x) = µo Jo . B L B = ŷB(x)
2 C
W x
Hence, we find that
( W W
ŷJo x, |x| < W2 −
2 2
H= By (x)
ŷJo W2 sgn(x), otherwise. µoJoW
2
Note that the solution plotted in the margin shows no discontinuity at x = ± W2
or elsewhere. −
W W x
2 2
2
B = ẑB
multi loop coil as depicted. A solenoid with its loop carrying a current I
in φ̂ direction (as shown), produces effectively a surface current density of
Js = IN φ̂ A/m, where N is the number density (1/m) of current loops in
C
the solenoid. In Example 3 we compute the magnetic field of the infinite
solenoid using Ampere’s law.
L I
The assumption of zero magnetic flux density B = 0 for the exterior region is justified
because:
(a) if the exterior field is non-zero, then it must be independent of x and y (follows
from Ampere’s law applied to any exterior path C with IC = 0), and
(b) the finite interior flux Ψ = µo IN πa2 can only be matched with the flux of
the infinitely extended exterior region when the constant exterior flux density
(because of (a)) is vanishingly small.
3
• Static electric fields: Curl-free and are governed by
∇ × E = 0, ∇ · D = ρ where D = ǫE
with ǫ = ǫr ǫo.
4
– Vector potential A can be used in magnetostatics in similar ways
to how electrostatic potential V is used in electrostatics.
◦ In electrostatics we can assign V = 0 to any point in space
that is convenient in a given problem.
◦ In magnetostatics we can assign ∇ · A to any scalar that is
convenient in a given problem.
– For example, if we make the assignment1
∇ · A = 0,
∇ × H = J ⇒ ∇ × B = µoJ,
it produces
∇2A = −µoJ,
which is the magnetostatic version of Poisson’s equation
ρ
∇2 V = − .
ǫo
1
With this assignment — known as Coulomb gauge — A acquires the physical meaning of “potential
momentum per unit charge”, just as scalar potential V is “potential enegy per unit charge” (see Konopinski,
Am. J. Phys., 46, 499, 1978).
5
– In analogy with solution
ρ(r′ )
Z
3 ′
V (r) = d r z J(r′)
4πǫo|r − r′ | r−r ′
∇ · A = 0 and ∇ × A = B.
Once A is available, obtaining B = ∇ × A is then just a matter of taking a
curl.
• Magnetic flux density B of a single current loop I can be calculated
after determining its vector potential as follows: z
y
– For a loop of radius a on z = 0 plane, we can express the corresponding current
density as
(−y ′ , x′, 0)
a
p
′ ′ ′2 ′2
J(r ) = Iδ(z )δ( x + y − a) p
x′2 + y ′2
x
I
where the ratio on the right is the unit vector φ̂′ .
r
– Inserting this into the general solution for vector potential, and performing
the integration over z ′ , we obtain Jφ = Iδ(z)δ( x2 + y 2 − a)
2
Also, in quasi-statics we use J(r′ , t) to obtain A(r, t) and B = ∇ × A over regions small compared to
λ = c/f , with f the highest frequency in J(r′, t).
6
µo I (−y ′ , x′, 0)
Z p
A(r) = ′2 ′2
δ( x + y − a) p p dx′ dy ′
4π (x − x′)2 + (y − y ′ )2 + z 2 x′2 + y ′2
µo I (−y ′, x′, 0)
Z
′
= δ(r − a) p r′ dr′ dφ′
4π (x − x′ )2 + (y − y ′ )2 + z 2 r′
Z π
µo I (−a sin φ′ , a cos φ′ , 0)
= p dφ′ ≡ x̂Ax (r) + ŷAy (r).
4π −π (x − a cos φ ) + (y − a sin φ ) + z
′ 2 ′ 2 2
7
◦ However, it can be easily seen that if x = 0 (i.e., along the z-axis), Bx = 0
(as symmetry would dictate) and
µo aI π a µo Ia2
Z
Bz = dφ = .
4π −π (a2 + z 2 )3/2 2(a2 + z 2 )3/2 3
For |z| ≫ a, 2
µo Ia2
Bz ≈ , 1
2|z|3
z
0
which is positive and varies with the inverse third power of distance |z|.
-1 I
– Also, Bx and Bz integrals can be performed numerically. Figure -2
600
H
• Note that circulation C B · dl around each closed field line (“linking” 400
the current loop) equals a fixed value of µoI — this dictates that the
200
z
0
lines closer to the current loop than on longer field lines linking the
loop further out. As a result |B| can be shown to vary as r−3 for large
-200
r. -400
-600
-600 -400 -200 0 200 400 600
• It can be shown that the equations for magnetic field lines of a current x
r = L sin2 θ
8
radial distance of the field line on θ = 90o plane, and the field line B∝I
formula is accurate only for r ≫ a. The Earth’s magnetic field had
qv × B
such a magnetic dipole topology as shown. qv × B
I
• Lorentz force due to the magnetic fields of a pair of current loops — also
known as magnetic dipoles — turns out to be “attractive” when the cur- I
rent directions agree (see margin). Bar magnets carrying “equivalent”
Loops with parallel
current loops of atomic origins interact with one another in exactly the currents attract one
another
9
14 Faraday’s law and induced emf
Michael Faraday discovered (in 1831, less than 200 years ago) that a chang-
ing current in a wire loop induces current flows in nearby wires — today
we describe this phenomenon as electromagnetic induction: the current
change in the first loop causes the magnetic field produced by the current to
change, and magnetic field change, in turn, is said to induce 1 (i.e., produce)
electric fields which drive the currents in nearby wires. Definitions of E and
• While static electric fields produced by static charge distributions are B have not changed:
unconditionally curl-free, time-varying electric fields produced by cur-
rent distributions with time-varying components are found to have, in recall that
accordance with Faraday’s observations, non-zero curls specified by • E is force per unit sta-
tionary charge
∂B
∇×E=− Faraday’s law • B gives an additional
∂t force v × B per unit
at all positions r in all reference frames of measurement. Using Stoke’s charge in motion with
velocity v in the mea-
theorem, the same constraint can also be expressed in integral form as surement frame.
∂B
I Z
E · dl = − · dS Faraday’s law
C S ∂t
for all surfaces S bounded by all closed and directed paths C (with
the direction of C, indicated by an arrow, and direction of vector dS
related by right hand rule).
1
Relativistic derivation of static B given in Lecture 12 can be extended to show that Coulomb interactions of charges
in time-varying motions require a description in terms of time-varying B and E — see, e.g., Am. J. Phys.: Tessman, 34,
1048 (1966); Tessman and Finnel, 35, 523 (1967); Kobe, 54, 631 (1986). Thus, the cause of induced E is not really the
time-varying B, but rather the time-varying current J that is also producing the variation in B.
1
• The right hand side of the integral form equation above includes the
flux of rate of change of magnetic field B over surface S. C(∆t) δS
E · dl = − B · dS, Ψ(0) = S
B(r, 0) · dS, and
C dt S Ψ(∆t) =
Z
B(r, ∆t) · dS +
Z
B(r, ∆t) · dS.
S δS
where the right hand side is now expressed in terms of the rate of change Ψ(∆t) − Ψ(0)
Thus, =
of magnetic flux Z
∆t
Z B(r, ∆t) − B(r, 0) Z dS
· dS + δS B(r, ∆t) · .
S ∆t ∆t
Ψ≡ B · dS
S Hence in limit ∆t → 0
linking contour C over any surface S bounded by C.
dΨ Z ∂B Z
= S · dS − C v × B · dl,
dt ∂t
• This modification (the exchange of the order of integration and time
since
derivative on the right side) would not be permissible if path C were Z dS Z ∆tv × dl
B(r, ∆t) · = C B(r, ∆t) ·
moving within the measurement frame or being deformed in time — δS ∆t
Z
∆t
=− (v × B) · dl
but in such cases we could still express Faraday’s integral form equation C
because
with − dΨ
dt
on the right side, provided that we also modify the left side B · (v × dl) = dl · (B × v),
as in both representing the volume of
d
I Z
a parallelepiped formed by the vectors
(E + v × B) · dl = − B · dS dl, v, and B.
C dt S Note that velocity v does not have to be
where v denotes the velocity of motion or deformation of path C. constant around contour C.
2
• A physical interpretation of the final equation
d
I Z
(E + v × B) · dl = − B · dS Integral form Faraday’s Law
C dt
| {z S }
dΨ
E =− Faraday’s eqn.
dt
is as follows: Magnetic field lines con-
tributing to Ψ form links
– the circulation integral on the left is the “voltage drop” once with path C (bounding S)
around the directed closed path C, representing the work done like the links in an ordinary
per unit charge (by the Lorentz force ∝ E + v × B) taken a chain — hence, Ψ is said to
full circle around C, which was denoted by Michael Faraday with a be the flux linking path C.
symbol E and called the emf (short for electro-motive force, which
B
is a bad name since E is work, and not force, per unit charge) for
the closed path, equaling the decay rate − dΨdt
of its linked mag-
v
netic flux Ψ (due to all sources of magnetic flux density B in the
region). dS
S
– if/when path C is occupied by a conducting wire loop of some
total conductance G = R1 , and a resistance R = G1 , a current C
3
driven by the non-zero field E + v × B within the wire accounting
for the non-zero E if/when − dΨ
dt is non-zero.
– in equivalent circuit models of conducting wire loops, Faraday’s
equation, re-written as RI = − dΨdt
, is effectively Kirchfoff’s voltage
law (KVL) applied to the loop, with RI on the left denoting the
(sum of all) voltage drops in the direction of C, while − dΨdt
on the
right denoting a voltage rise also in the direction of C.
◦ note that the emf E describes both the voltage drop RI and
voltage rise − dΨ
dt
appearing in the circuit model for the con-
ducting wire loop since E = RI and E = − dΨ dt are both true.
– in modern parlance (since Maxwell) the term emf and its symbol
E are used to refer to and denote sources of energy, e.g., battery
voltages and magnetic flux rate − dΨdt that drive currents I = R
E
B
around closed circuits3 .
v
• If path C is fixed in the measurement frame, then v = 0, and KVL for
such a stationary loop reads as dS
S
dΨ
I
E · dl = − ; C
C dt
4
dΨ
I
(E + v × B) · dl = −
C dt
because in that case force per unit charge advected with path C S
will be E + v × B according to Lorentz force (note: any additional
velocity vq of a moving charge along C does not contribute because C I=
E
R
(vq × B) · dl = 0 if dl and vq are parallel).
Think of EMF as the sum of all the "voltage
rises" around the loop traversed in the
– In either case, if C is a physical conducting path with a total direction of loop current I that needs to
match the total "voltage drop" RI around
I = dt is fulfilled.
R
around C in the circulation direction (determined by dl and dS
directions used in accordance with the right-hand-rule).
5
1. Fixed C, but time-varying B,
2. B =const. (in space and time), but time-varying C (rotating or
changing size),
3. An inhomogeneous static B = B(r) in the measurement frame
and C in motion.
B
• Note that even in the absence of any electric field E in the measurement
frame, a non-zero emf
v
dΨ
I
(v × B) · dl = −
C dt S
E
can exist because of the motion of C through an inhomogeneous mag- C I=
R
netic field (if the field is homogeneous then dΨ
dt will be zero, implying
zero E), which will of course appear as an emf
dΨ′
I
′ ′
E · dl = − ′
C dt
for a second observer moving with C who sees a time varying electric
field E′ = v × B in her own frame (in addition to the inhomogeneous
but constant magnetic field B of the first frame appearing as a time-
varying magnetic field B′)5 .
E′
See Scanlon et. al., Am. J. Phys., 37, 698 (1969), for a discussion of I ′ =
5
R
for rigid C with resistance
R observed from different reference frames.
6
– Thus, having non-zero electric field circulations
I
E′ · dl′
C
7
z
Boe−t/τ ẑ
y
Example 1: If
B = Bo e−t/τ ẑ,
what is the emf E taken over a stationary circular loop C of radius r = 10 m on x
z = 0 plane in counter-clockwise direction (looking down on z = 0 plane)? What
is current I if the loop resistance is R? C
8
Example 2: Consider the magnetic flux density
µo I
B= φ̂
2πr
y
produced by current I flowing along the x axis. What is the emf E of a square
loop C of area 4 m2 moving on xy-plane with edges parallel to x- and y-axes, if 2 m/s
its center is located at
H y = 2t m as a function of time? Compute the emf E first
dΨ
as − dt and then as C (v × B) · dl to verify that the same values are obtained. 2m
2t
Solution: Given the described geometry, we have C
Z 1 Z 2t+1
µo I µo I 2t + 1 I x
Ψ(t) = dx dy = ln( ).
−1 2t−1 2πy π 2t − 1
Thus, the emf E is
dΨ µo I 2t − 1 ∂ 2t + 1 µo I 4 µo I
− =− ( ) ( )= = .
dt π 2t + 1 ∂t 2t − 1 π (2t + 1)(2t − 1) π(t2 − 41 )
µo I
Alternatively, since v = 2ŷ m/s, and v × B = 2 2πr x̂, we find, using dl = ±x̂dx
and ±ŷdy in turns,
µo I µo I µo I
I
E = (v × B) · dl = 2 2−2 2=
C 2π(2t − 1) 2π(2t + 1) π(t2 − 41 )
in consistency with the above result.
9
Example 3: A conducting loop of a radius r = 0.1 m (see figure in the margin) is
ω
being rotated about the x axis with frequency of f = 2π = 60 Hz in a region
with a DC magnetic field of B = 10ẑ T. Determine the induced current in the
loop if the loop resistance is 12 Ω.
Solution: The maximum value of the magnetic flux linking the loop should be
10
z
C ~v = 3x̂ m/s
X X X X X X X
Example 4: A conducting bar of resistance R1 = 1 Ω ohms is moved in the x-direction R2 2m
R1 = 1 Ω
with a velocity v = 3x̂ m/s on a pair of perfect conducting (R = 0) stationary X X X X X X X
x
3t
rails 2 m apart terminated with a load resistance R2 at x = 0, all constituting a
rectangular contour C to be taken counterclockwise. A constant magnetic field
of B = 1ŷ T is linked throught contour C such that the flux Ψ = −1 × 2 × 3t
and the emf E = −dΨ/dt = 6 V. Hence, Faraday’s law demands that
I Z t Z b
(E + v × B) · dl = (E + v × B)1 · dl + (E)2 · dl = 6 Moving bar in the presence
C b t of a constant magnetic field
produces an emf and electric
where the two integrals (with b and t referring to bottom and top rail contact fields in the lab frame that
points) correspond to voltage drops across resistors R1 and R2 , respectively. But drive a loop current I.
since Z t
(v × B)1 · dl = 3 × 1 × 2 = 6, Example
H 4 illustrates
H how the
b E · dl part of emf (E + v ×
it follows that B) · dl caused by a motion v =
Z t Z b 3x̂ m/s is zero (with non-zero
(E)1 · dl + (E)2 · dl = 0 ⇒ Ez1 − Ez2 = 0 ⇒ Ez2 = Ez1, static Ez components)!!
b t
i.e., identical static fields within the moving and stationary bars across the perfect
conducting rails. This may be a surprising claim/result — let’s give two examples
to illustrate how this happens:
1. Let R2 = 2 Ω ohms. Then I = 6/3 = 2 A. It follows that voltage drops (E + v × B)1 · 2ẑ = 2
V across R1 and (E)2 · (−2ẑ) = 4 V across R2 , yielding Ez1 = Ez2 = −2 V/m.
2. Let R2 = ∞ — open ckt load to the moving conductor. Then I = 6/∞ = 0 A. It follows that
(E + v × B)1 · 2ẑ = 0 V across R1 and (E)2 · (−2ẑ) = 6 V across R2 , yielding Ez1 = Ez2 = −3
V/m. Note that in this
R t case the entire emf appears across the open termination (gap in the
loop C and the emf b (E + v × B)1 · dl = 0 across resistor R1 ).
11
Example 5: An infinite solenoid producing a constant − dΨ dt = 8 V, passes through
small a loop consisting of a 1 Ω resistor on the right and a 3 Ω resistor on the
left, connected in series — see margin plot. What is the current Ic through this
resistor loop, and what voltages would be measured (by a voltmeter) across the
individual resistors?
Solution: The magnetic flux produced by the solenoid will be confined to its interior
Transformers which operate
as long as dI/dt (and thus dΨ/dt, as specified) is constant and emf E = −dΨ/dt based on an inductive cou-
is non-time varying (see below). In that case, with constant emf E = − dΨ
dt = 8
pling principle, and electric
dynamos (and motors) which
V of the encircling resistor loop in the setup, the loop current Ic is the ratio of produce motion induced emfs
E (and rotating coils) are stud-
E and the total loop resistance 4 Ω, i.e., Ic = R = 2 A. Consequently, 1 and 3
ied in depth in power courses
Ω resistors will develop 2 and 6 V drops, respectively, in the direction of the 2A starting with ECE 330.
current!! Note that:
• the loop has no battery to support this current flow — it has instead been excited
“inductively”.
• with constant dI/dt, there is zero magnetic field at the locations of the loop wire and
resistors (static E in the solenoid exterior is curl-free!) — thus, the emf of the loop is
not being produced by a time varying local magnetic field; it is rather a consequence of
the time-varing current I(t) in the solenoid (which is also responsible for time-varying
Ψ), with the relation E = −dΨ/dt being “incidental”!
• what a voltmeter measures across the resistors — whether 2 or 6 V — depends on
whether its probes contacting points A and B are placed to the right or to the left of
the solenoid!! That’s because the field E produced by the time-varying current I(t) is
no longer conservative across the system and consequently the line integral of E is path
dependent — we have to be more careful about what we mean by voltage in these new
situations!
12
y
A
x
B = (12 − 3t)ẑ T.
• Voltage drop V2 = −1 V across R2 from point A to point B. A voltmeter reads and dis-
plays the voltage of its
• A voltmeter connected from A (positive lead) to B will read 2 V if and only if its own path constituted by the
leads form a path identical to the path defined by R1 (from A to B). placement of its own probe
wires contacting the mea-
• A voltmeter connected from A (positive lead) to B will read -1 V if and only if surement nodes A and B.
its leads form a path identical to the path defined by R2.
• A voltmeter connected from A (postive lead) to B will read 0.5 V if its leads form
a diagonal path from A to B.
– To see this, notice that Faraday’s law applied for the triangular loop in-
cluding the voltmeter and R2 would have an emf of 1.5 V equaling the sum
of voltmeter reading VR and 1 V drop across resistor R2 .
13
(a) A one turn coil with current I
generates its own linked magnetic
15 Inductance — coil, solenoid, shorted coax flux LI as shown, where a non-
negative L is the inductance of
the coil.
• Given a circular coil with some resistance R and conducting some cur- 3
rent I, the magnetic flux Ψ produced by I and “linking” the coil itself
2
Ψ = LI Ψ
z
0
L=
I -3
-3 -2 -1 0 1 2 3
1
(a) A one turn coil with current I
– The current I and self-emf E are then the solutions of differential generates its own linked magnetic
flux LI as shown, where a non-
equations negative L is the inductance of
the coil.
dI dE
RI = −L and RE = −L ,
dt dt 3
z
0
in RC circuits).
-1 I
◦ Note that τ = L/R implies that when the inducance L is
large, so is time constant τ , and current decay in the induc-
-2
time-varying current sources (just like capacitors behaving like (b) An equivalent circuit model
for the coil expressed in terms of
time-varying voltage sources) as they relase their stored energy lumped resistor R and inductor L
forming a loop carrying the loop
(while maintaining a voltage rise −L dI dt determined by other
current I R
elements in their connected circuits).
dI
• For an inductor consisting of n-loops, the emf E measured across all I, E = −L
dt
n-loops is naturally (since n emf’s add up)
+ dI -
V (t) = L
d d dI dt
E = n(− Ψ) = − nΨ ≡ −L The emf RI=-LdI/dt of the coil
dt dt dt appears as a voltage rise across
the inductor in the ckt model,
as well as a voltage drop across
implying an inductance the resistor, both taken in the
direction of current I. Voltage
nΨ drop V across the inductor in the
L≡ . current direction is LdI/dt, as
I we learned in our circuit courses.
2
Example 1: An n-turn coil has a resistance R = 1 Ω and inductance of 1 µH. If it is
conducting 3 A current at t = 0, determine I(t) for t > 0.
Solution: Current flow in the resistive n-turn coil will be driven by self-emf E = −L dI
dt
matching a voltage drop RI. Hence
dI dI R R 6
RI = −L ↔ + I = 0 ⇒ I(t) = I(0)e− L t = 3e−10 t A.
dt dt L
3
Inductance of long solenoid: Consider a long solenoid with length ℓ, B = ẑB
cross-sectional area A, and a density of N loops per unit length as examined
in Example 3 of Lecture 12 (see figure in the margin). As determined in
Example 3, the magnetic flux density in the interior of the solenoid is
B = µoIN ẑ
4
• Notice that the stored energy of the solenoid is
1 1
µo|Hz |2 = µoH · H
2 2
times its volume Aℓ occupied by the field H inside the solenoid. That
suggests that
1
w = µo H · H
2
can be interpreted as stored magnetostatic energy per unit volume in
general.
µ = (1 + χm )µo
5
Inductance of shorted coax: Consider a coaxial cable of some length ℓ I
which is “shorted” at one end (with a wire connecting the inner and outer
r Short
conductors), so that a steady current I can flow on the inner conductor of b a I
radius a to return on the interior surface of the outer conductor at radius z
b after having circulated through the short. We will next determine the B
inductance L of such an inductor after first computing the magnetic flux
density Bφ that will be produced by the inner conductor current I. In Bφ
calculation we will assume ℓ ≫ b so that an “infinite coax” approximation ℓ
can be invoked. Shorted coax circulates
a current I linking a
• Expanding the integral form of Ampere’s law magnetic flux Ψ
confined to a region
I bounded by the outer
conductor of the coax.
B · dl = µoIC
C
as
Bφ2πr = µoI
over a circular integration contour C of a radius r > a, we find that
the magnetic flux density in the interior of the coax cable is
µo I
Bφ = .
2πr
• Therefore, the magnetic flux linked by the closed current path I (see
figure in the margin) is
Z b
µo dr µo b
Z
Ψ = B · dS = ℓ I = ℓ ln I.
S 2π a r 2π a
6
Clearly, we have a linear relation Ψ = LI, with
ln ab
L≡ ℓµo,
2π
which is the inductance of a shorted coax of a finite length ℓ.
ln ab
L= µo ,
2π
which should be contrasted with capacitance per unit length
2π
C= ǫo
ln ab
of the same coax configuration.
Notice how L and C are proportional to ǫo and µo, respectively,
having proportionality constants which are inverses of one another.
7
Inductance of shorted parallel plates: If a pair of parallel plates of
areas A = W ℓ and separation d were shorted at one end, we would obtain
effectively an inductor with a per length inductance
d
L= µo
W
that accompanies per length capacitance
W
C= ǫo
d
of the same parallel plate configuration. This follows from a generalization of
our finding above that the proportionality constants of L and C are arithmetic
inverses of one another.
8
16 Charge conservation, continuity eqn, displace- (a) At t=0 volume V contains
a neutral atom but no net charge
1
• Consider two distinct surfaces S1 and S2 bounded by the same closed
loop C (as shown in the margin) such that a volume V is contained
between the two surfaces.
– Let dS2
Z S2
I1 = J · dS1 V dS
S1
and Z dS1
I2 = J · dS2 C
S2
S1
denote currents flowing through surfaces S1 and S2, respectively.
– Note that current I1 through surface S1 enters volume V , while
current I2 through surface S2 exits volume V (with the directions
assigned to dS1 and dS2).
– If I1 6= I2, then current out is not matched by the current in,
and as a result, the net charge QV contained in volume V increases
with time at a rate I1 − I2 provided that charge is conserved in
the sense discussed above. In that case, we have
dQV
= I1 − I2 .
dt
This relationship can be expressed as
d
Z Z Z
ρdV = J · dS1 − J · dS2
dt V S1 S2
2
since, in terms of charge density ρ, charge in volume V is
Z
QV = ρdV.
V
3
• While Faraday’s law
∂B
∇×E=−
∂t
indicates that time-varying B induces time-varying electric fields E,
Ampere’s law, written as
∇ × H = J,
B = µoH.
4
Verification of Maxwell’s claim: Since ∇ × H is divergence-free
(just like the curl of vector potential A, namely B), it follows that
the divergence of Maxwell’s modified Ampere’s law — often called
Ampere-Maxwell equation — is
∂
∇ · (∇ × H) = ∇ · J + ∇ · D = 0.
∂t
– In the absence of the second term due to displacement current,
this results would be inconsistent with the continuity equation
∂ρ
+ ∇ · J = 0,
∂t
∂ρ
unless ∂t = 0 (the static case).
– By, contrast, including the second term, the result above is rec-
ognized as the continuity equation per se, since by Gauss’s law
— assuming that it applies with no change under time varying
situations —
∂ ∂ρ
∇·D= .
∂t ∂t
• The modified Ampere’s law
∂D
∇×H=J+
∂t
postulated by Maxwell under the assumption that Gauss’s law is also
valid under time-varying conditions, leads to some specific predictions
about how time-varying fields should behave.
5
• These predictions — concerning the propagation of electromagnetic
waves — were validated experimentally by Heinrich Hertz around 1888.
6
in terms of suitably defined permittivities and permeabilities ǫ and
µ — see next Lecture.
∇·B=0
7
magnetic fields” in nature, ∇ · B would have equaled the density of those
charges, and magnetic field lines would have started and stopped on them
(rather than looping into themselves). But no one has observed of any evi-
dence for such magnetic charges anywhere, even going back to the very early D+
n̂
times in the history of the universe (accessible by making observations of w D−
very far astronomical objects). So, ∇ · B = 0.
– We had already seen how the first and third boundary condition
Constraint
equations arise. I
B · dS = 0
– The second boundary condition equation concerning the normal S
component of B is another consequence of the absence of magnetic applied over the dotted volume (seen in
profile) yields
charges (see margin).
Bn+ − Bn− = 0
– A detailed justification of the last boundary condition concerning
tangential H will be given explicitly during Lecture 19. This equa- in w → 0 limit.
8
17 Magnetization current, Maxwell’s equations
in material media
• Consider the microscopic-form Maxwell’s equations
∇·D = ρ Gauss’s law
∇·B = 0
∂B
∇×E = − Faraday’s law
∂t
∂D
∇×H = J+ , Ampere’s law
∂t
where
D = ǫo E
B = µoH.
1
– subscripts f indicate charge and current density contributions due
to free charge carriers,
– the term −∇ · P denotes the bound charge density,
– the term ∂P∂t denotes the polarization current density due to
oscillating dipoles (already discussed in Lecture 11), and
– ∇×M is a “magnetization current density” also due to bound
charges, an effect that we will discuss and clarify later in this
section.
∇ · ǫo E = ρ Gauss’s law
∂ǫoE
∇ × µ−1
o B = J + , Ampere’s law
∂t
we obtain
D = ǫeE + P = ǫE
and
H = µ−1 −1
o B − M = µ B,
2
respectively, and drop the subscripts f which will no longer be needed.
Using these new definitions, the full set of Maxwell’s equations now
read as (the same form as before)
∇·D = ρ Gauss’s law
∇·B = 0
∂B
∇×E = − Faraday’s law
∂t
∂D
∇×H = J+ , Ampere’s law
∂t
with
D = ǫE
B = µH,
where ρ and J are understood to be due to free charge carriers only.
• We had already seen many aspects of the above procedure for obtaining
the macroscopic form field equations earlier on (e.g., in Lectures 8 and
11).
3
• Just like “free charge” density and currents, “bound charge” densities
and currents also have to satisfy the continuity equation
∂ρ
+ ∇ · J = 0.
∂t
– This equation is automatically satisfied if we substitute
ρ = ρb = −∇ · P
and
∂P
J = Jb =
∂t
in it.
Verification:
∂ρb ∂ ∂P
+ ∇ · Jb = (−∇ · P) + ∇ · =0
∂t ∂t ∂t
since the order of time derivative and divergence can be exchanged
on the right.
– But the same equation is also satisfied if we take
∂P
Jb = +∇×M
∂t
for any vector field M simply because vector ∇ × M is divergence
free.
4
Consequently, it is not sufficient to represent bound currents in mate-
rial media as simply ∂P
∂t , if bound carriers can also conduct divergence-
free currents due to closed-loop orbits.
H = µ−1
o B−M
1
Physical models of nucleons involve bound charge carriers known as quarks which cannot be observed
in a free state.
2
All elementary charge carriers carry an intrinsic magnetization proportional to charge-to-mass ratio
q
m
and a “spin angular momentum” having quantized values of ± ~2 N.m.s in any measurement direction.
Using Heisenberg’s uncertainty principle, ∆p∆r ≥ ~2 , we can interpret the spin angular momentum of
an elementary particle as the lower bound of ∆p∆r, the product of quantum uncertainties in particle
momentum and position. There is no classical interpretation of spin angular momentum for point particles.
5
• To get a physical picture about magnetization M and the physical
origin of H = µ−1 o B − M consider a solenoid wound around some
cylindrical shaped material as shown in the margin. We know that
with a solenoid current Io, we would have Ho = N Io ẑ in the interior of
a solenoid with N loops per unit legth aligned with the z-axis, and a
corresponding magnetic flux density Bo = µoN Ioẑ when the solenoid
Io
core is free space. This will be modified to some B = Bo + µoM when
a material core is introduced into the same space, where µoM stands Within the core, stacks of
atomic loop currents are
effective solenoids giving
for the (additional) macroscopic (space averaged) magnetic flux den- rise to an additional
magnetic flux density of
an average magnitude
sity produced by microscopic current loops localized within the atoms Bo = µoIoN
inside the solenoid
µo Al Il Nl = µoM
but outside the cylindrical that adds to Bo
constituting the core. magnetic core
1
– If there are Na = ∆x∆y∆z atoms per unit volume in the core, with
∆x separations in x direction and so forth, loop currents Il of a
stack of atoms with ∆z separations in z would produce an effective
Il
solenoid an internal z-directed magnetic flux density of µo ∆z ẑ and
zero exterior field.
– Since one such atomic stack solenoid with a loop area of Al will be
found for every ∆x∆y cross-sectional area of the core, a macro-
scopic average magnetic flux density produced by these atomic
solenoids would be calculated as (this calculation is similar to
finding the average polarization field in a dielectric as discussed in
Al Il
Lecture 8) ∆x∆y × µo ∆z ẑ = µoNaIl Al ẑ ≡ µoM, with M = Na m,
m ≡ Il Al ẑ.
6
◦ Here m is the magnetic dipole moment of each current loop
(analogous to electric dipole p = qr), M is the magnetization
field vector (analogous to P = Nap), which is a simple product
of m per magnetized atom and the atomic number density Na
in the core.
– Superposing the magnetic flux densities of µoM andBo, we obtain
B = Bo +µoM for the core region, or for any region of space having
a non-zero magnetization M, which then leads to the general result
H = µ−1o B − M, which is further discussed below.
– Notice, whether the flux density B = Bo +µoM inside the material
medium is stronger or weaker in magnitude than Bo depends on
the direction of M, which, in turn, depends on the algebraic sign
of microscopic loop currents Il introduced above.
◦ Negative Il is found in diamagnetic materials where |B| <
|Bo|, while positive Il in paramagnetic and ferromagnetic
materials where |B| > |Bo|, as discussed below.
– Also, the expression H = µ−1 −1
o B − M leads to H = µo B = Ho
in the exterior region where M = 0, indicating that while fields B
and Boif the interior and exterior are different, H is the same in
both regions (analogous with D in dielectrics).
7
• Lab measurements — e.g., inductances L measured for coils wound
around magnetic materials3 — show that for a large class of materials
M ≡ µ−1
o B−H
3
Recall from Lecture 15 that L ∝ µ when inductors are wound around materials with permeability µ.
8
• For a large class of materials with M ∝ H, it is observed that |χm| ≪ 1.
In that case, the material is called
– Diamagnetic if χm < 0:
◦ Diamagnetism occurs when an applied magnetic field induces
electron orbital angular momentum in a collection of atoms
having no net permanent magnetization M — in such materi-
als electron clouds around atomic nuclei spin up in accordance
with Lenz’s to generate magnetic fields opposing the applied
magnetic field so as to keep B = µH smaller than µoH. This
happens in materials that we ordinarily think of being non-
magnetic (wood, glass, water, etc.). Diamagnetic materials
are in fact very weakly repulsed by permanent magnets since
µ ≈ µo in all diamagnetic materials.
– Paramagnetic if χm > 0:
◦ Paramagnetism occurs in materials composed of atoms hav-
ing permanent magnetic dipole moments due to electron spin
angular momentum — magnetic dipoles of such atoms co-
align with the applied magnetic field due to v × B related
torques, leading to M pointing in the applied B direction4 .
This happens for atoms with unfilled inner electron shells, be-
cause in filled shells electron spins are opposite (due to Pauli
4
In these materials the described paramagnetism overcomes the diamagnetic tendency of the material
caused by the orbital angular momenta of its electrons around atomic nuclei.
9
exclusion principle) and cancel one another. Unfilled outer
shells do not usually give rise to paramagnetism because in-
teractions between adjacent atoms in that case give rise to
opposite spins of their outer shell electrons. Paramagnetic
materials are very weakly attracted to permanent magnets
(e.g., aluminum, lithium, tungsten).
B = B(H),
10
els of χe examined in Lecture 11) relevant for applications including
various magnetic imaging techniques.
11
18 Wave equation and plane TEM waves in source-
free media
With this lecture we start our study of the full set of Maxwell’s equations
shown in the margin by first restricting our attention to homogeneous and
non-conducting media with constant ǫ and µ and zero σ.
• Our first objective is to show that non-trivial (i.e., non-zero) time- ∇·D = ρ
varying field solutions of these equations can be obtained even in the ∇·B = 0
absence of ρ and J. ∂B
∇×E = −
∂t
– We already know static ρ and J to be the source of static electric ∂D
∇×H = J+ .
and magnetic fields. ∂t
– We will come to understand that time varying ρ and J, which
necessarily obey the continuity equation
∂ρ
+ ∇ · J = 0,
∂t
constitute the source of time-varying electromagnetic fields.
1
• Such field solutions in fact represent electromagnetic waves, a familiar
example of which is light.
– radiowaves,
– microwaves,
– infrared,
– light,
– ultraviolet,
– X-rays, and gamma rays,
2
• In source-free and homogeneous regions where ρ = J = 0 and ǫ and
µ are constant, we can simplify Maxwell’s equations as shown in the
margin.
∇ × (∇ × E) = ∇(∇ · E) − ∇2E
which should be familiar from an earlier homework problem.
4
satisfies the 1D wave equation and represents an x-polarized time-
periodic field solution with an oscillation frequency ω.
– 1D wave equation can also be satisfied by
√
Ex = cos(ω(t + µǫz)).
• Let us next find out the magnetic field intensity H that accompanies
the x-polarized electric field wave solution
z
E = x̂ cos(ω(t ∓ )).
v
– Since the curl of E is
x̂ ŷ ẑ
∂ ∂ ∂ ∂Ex z ω
∇×E= ∂x ∂y ∂z = ŷ = ±ŷ sin(ω(t ∓ )) ,
∂z v v
Ex 0 0
5
Faraday’s law z
∂H y
∇ × E = −µ
∂t E×H
requires that H should satisfy
∂H z ω H
−µ = ±ŷ sin(ω(t ∓ )) . z
∂t v v E = x̂f (t − )
v
Finding the time-dependent anti-derivative (and remembering v = x
√
1/ µǫ), we obtain
r
ǫ z z y
H = ±ŷ cos(ω(t ∓ )).
µ v
6
• Since Maxwell’s equations with constant µ and ǫ are linear and time-
invariant (LTI), the field solutions above can be further generalized by
using their weighted and time-shifted superpositions such as
X
f (t) = An cos(ωn t + θn )
n
and ∞
1
Z
f (t) = F (ω)ejωtdω
2π −∞
having frequency dependent weighting factors An and F (ω). And since
according to Fourier analysis all practical signals f (t) can be synthe-
sized in these forms, it follows that the field solutions above are valid
with arbitrary waveforms f (t). d’Alembert
wave
Solutions solutions
z
E, H ∝ f (t ∓ )
v
of the 1D scalar wave equation with arbitrary f (t) are known as d’Alembert
wave solutions.
• d’Alembert solution
z
E, H ∝ f (t − )
v
describes electromagnetic waves traveling in +z direction, whereas so-
lution
z
E, H ∝ f (t + )
v
7
Traveling wave in +z direction with speed v=c:
v η pulse
τ 0 τ t 0 τ to τ t
− −
Question: What about z-polarized waves 2 2 2 2
t t − to
z △( )
τ
△(
τ
)
E = ẑf (t ∓ ), 1
Triangle
1
v pulse
τ 0 τ t 0 τ to τ t
can they exist? −
2 2
−
2 2
t[u(t) − u(t − to)]
8
z
19 d’Alembert wave solutions, radiation from y
current sheets E×H
these wave-fields are constant (have the same vector value) at planes
of constant phase, e.g., on planes defined by
z
t − = const.,
v
1
which are planes transverse to the propagation direction (direction of
vector E × H).
Not all waves solutions of Maxwell’s equations are uniform plane — for in-
stance non-uniform TEM waves with spherical surfaces of constant phase are
ubiquitous, but they will be examined later on (in ECE 450, mainly).
After the next set of examples we will examine how uniform plane waves
can be radiated by infinite planes of surface currents. By contrast, spherical
waves are produced by compact antennas having finite dimensions.
Example 1: Let
t − y/c
E = x̂△( )
τ
be a wave solution in free space where △( τt ) is a triangular waveform of duration
τ peaking at t = 0 (defined in ECE 210). We will next provide two different
solutions demonstrating how the wave field B accompanying E can be found.
Solution 1: We recognize the given wave field E as a TEM uniform plane wave travel-
ing in y-direction given the t−y/c dependence of phase. Consequently, we obtain
H by dividing E with η = ηo and rotating it by 90◦ from x̂-direction to co-align
it with E × H vector. As a result,
△( t−y/c
τ ) △( t−y/c )
H = −ẑ = −ẑ p τ .
ηo µo /ǫo
Hence,
√ t − y/c △( t−y/c
τ
)
B = µo H = −ẑ µo ǫo △( ) = −ẑ .
τ c
2
Solution 2: According to Faraday’s law,
x̂ ŷ ẑ
∂B ∂ ∂ ∂ ∂Ex
= −∇ × E = − ∂x ∂y ∂z = ẑ
∂t ∂y
Ex 0 0
t − y/c ∂ t − y/c −1 t − y/c
= ẑ△′ ( ) ( ) = ẑ △′ ( )
τ ∂y τ cτ τ
d
with the help of chain rule of differentiation, where △’(u) ≡ du △(u). Finding
the time-dependent anti derivative, we directly obtain (as before)
△( t−y/c
τ )
B = −ẑ .
c
Solution: Since
t − y/c △( t−y/c
τ )
E = x̂△( ) and B = −ẑ ,
τ c
it follows that Lorentz force
t − y/c ẑ
F = q(E + v × B) = q△( )(x̂ − v × ).
τ c
3
Clearly, the first term of F proportional to x̂ is dominant, unless v = |v| is close
to c.
Example 3: Consider an x̂-polarized plane TEM wave field in free space propagating
in +z direction such that
z t
E(z, t) = x̂f (t − ), with f (t) = At rect( ),
c τ
where c = 3 × 108 m/s = 300 m/µs is the speed of light in free space, τ = 1 µs,
V/m
and A = 2 µs . A plot of f (t) vs t (labelled in µs units is shown in the margin.
Determine the corresponding H(z, t) and make the following plots:
Solution: (a) t-plots at fixed z’s: Since z/c = 2 µs for z = 600 m, it follows that
t − 2µs V
Ex (600 m, t) = 2 (t − 2µs) rect( )
1µs m
is a shifted version of
t V
Ex (0, t) = 2 t rect(
)
1µs m
already plotted above. A graph showing both waveforms (black for z = 0 and
red for z = 600 m) is in the margin.
4
(b) z-plots at fixed t’s: In this case we wish to deptict
z 0 − z/c V
Ex (z, 0) = 2 (0 − ) rect( )
c 1µ m
and
z 2µ − z/c V
Ex (z, 2 µs) = 2 (2µ − ) rect( ) .
c 1µ m
The minus sign in front of z in the first term on the right indicates that the slopes
of the curves to be plotted are negative. Hence, we end up with the descending
ramp waveforms (black for t = 0 and red for t = 2µs) shown in the margin.
• The answer must be, far away ρ and J variations (linked by continuity
equation) that we have not considered in our equations so far.
• We will next describe how plane TEM waves can be produced — radi-
ated — by time-varying infinite current sheets by starting from familiar
static and quasi-static solutions:
5
z
• Consider first a static and constant surface current density
Js = x̂Jx A/m y
+
H
flowing on z = 0 surface as shown graphically in the margin. This
infinite surface current will produce a static magnetic field Js = x̂Jx
Jx x
A/m for z ≷ 0
H(z) = ∓ŷ
2
also shown in the margin as we learned in Lecture 13. H−
– Note that the fields point in opposing directions above and below
the surface current in compliance with the right hand rule and
obey the boundary condition equation for tangential H.
– Also, H is not accompanied by an electric field E since static
currents produce only static magnetic fields.
6
• But the exact field solution of Maxwell’s equtions valid for all z is z
equally easy to obtain: just replace Jx (t) above with Jx(t ∓ vz ) and E+ × H+
E+
replace ≈ with = so that y
z H+
Jx (t ∓ v)
H(z, t) = ∓ŷ A/m for z ≷ 0
2 Js = x̂Jx (t)
2
We use Jx (t ∓ vz ) rather than Jx (t ± vz ) for z ≷ 0 because we assume that Jx (t) on z = 0 surface is the
only field source — in that case causality principle dictates that we use only the solutions propagating
away from the source (just like when a pebble drops in a pond, ripples propagate away).
7
Conclusion: Evidently, a time varying surface current
± ηf (t ∓ vz ) ± f (t ∓ vz )
E = −x̂ and H = ∓ŷ in regions z ≷ 0
2 2
propagating away from the z = 0 plane.
Note that: z
E+ × H+
1. Ex and Hy waveforms are proportional to delayed versions of surface E+
y
current Jx(t) at each location z above and below the current sheet, with +
H
the reference directions of E and Js opposing one another.
Js = x̂Jx (t)
±
2. fields E are continuous on z = 0 surface in compliance with tangential
x
boundary condition equations.
H−
3. fields H± exhibit a discontinuity on z = 0 surface that matches the E−
current density of the same surface, once again in compliance with
E− × H−
tangential boundary condition equations.
8
Example 4: A current sheet on z = 0 surface is described by
t
Js (t) = x̂f (t), with f (t) = At rect( ),
τ
A/m
where τ = 1 µs and A = 1 µs . A plot of the current waveform f (t) is plotted in
the margin. Assuming that the current sheet is embedded in free space, construct
the following plots:
(b) We have, multiplying Hy with ηo = 120π Ω, and adjusting the signs so that E
and Js are pointing in opposite directions,
z 2µ ∓ zc V
Ex (z, 2µs) = −60π(2µ ∓ ) rect( ) for z ≷ 0.
c 1µ m
Plots are shown in the margin.
9
20 Poynting theorem and monochromatic waves
• The magnitude of Poynting vector
S=E×H
1
∂ 1 1
( ǫE · E + µH · H) + ∇ · (E × H) + J · E = 0.
∂t 2 2
2
is energy flux per unit area measured in
VA W J/s
= 2= 2
mm m m
3
z
y
Example 1: On z = 0 plane we have a time-harmonic surface current specified as S+
+ η z
E = −x̂ f (t − )
2 v Js = x̂f (t)
A
Js = x̂f (t) = x̂2 cos(ωt)
m
where ω is some frequency of oscillation.
1 z
η z H− = ŷ f (t + )
(a) Determine the radiated TEM wave fields E(z, t) and H(z, t) in the regions z ≷ 0, E = −x̂ f (t + )
−
2 v
2 v
S−
(b) The associated Poynting vectors E × H, and
Solution: (a) With reference to the solution of the current sheet radiator depicted
in the margin (from last lecture), we that an x-polarized surface current f (t)
produces the wave fields
η z 1 z
Ex = − f (t ∓ ) and Hy = ∓ f (t ∓ )
2 v 2 v
in the surrounding regions propagating away from the current sheet on both sides.
Given that f (t) = 2 cos(ωt), this implies that
Ex = −η cos(ωt ∓ βz) and Hy = ∓ cos(ωt ∓ βz)
where
ω
β= and η = ηo ≈ 120π Ω
c
since the current sheet is surrounded by vacuum. Hence in vector form we have
V A,
E(z, t) = −η cos(ωt ∓ βz)x̂ and H(z, t) = ∓ cos(ωt ∓ βz)ŷ m
m
where the upper signs are for z > 0, and lower signs for z < 0.
4
(b) The associated Poynting vectors are
W
S = E × H = ±η cos2(ωt ∓ βz)ẑ .
m2
Note that the time-average value of vector S points in the direction of wave
propagation on both sides of the current sheet.
(c) Since on z = 0 surface of the current sheet the electric field vector is
V
E(0, t) = −η cos(ωt)x̂ ,
m
it follows that Js · E on the same surface is
A V W
Js (t) · E(0, t) = (x̂2 cos(ωt) ) · (−η cos(ωt)x̂ ) = −2η cos2 (ωt) 2 .
m m m
5
can express the time-average of this identity as z
y
1 1 S+
hcos2(ωt + φ)i = h [1 + cos(2ωt + 2φ)]i = , + η z
E = −x̂ f (t − )
2 2 2 v Js = x̂f (t)
6
W
−Js · E = 2η cos2(ωt) .
m2
Clearly, its time-aveage works out as
W W
h−Js · Ei = η = 120π .
m2 m2
– Note that h−Js · Ei exactly matches the sum of |hE × Hi| calcu-
lated on both sides of the current sheet, in conformity with energy
conservation principle (Poynting theorem).
7
21 Monochromatic waves and phasor notation
• Recall that we reached the traveling-wave d’Alembert solutions
z
E, H ∝ f (t ∓ )
v cos(ωt)
Period
1 2π
via the superposition of time-shifted and amplitude-scaled versions of T =
ω
f (t) = cos(ωt),
t
namely the monochromatic waves
-1
z
A cos[ω(t ∓ )] = A cos(ωt ∓ βz), cos(βz)
Wavelength
v 1 2π
λ=
with amplitudes A where β
ω √
β≡ = ω µǫ z
v
can be called wave-number in analogy with wave-frequency ω. -1
1
Period
– In general, monochromatic solutions of 1D wave-equations ob- 1
cos(ωt)
2π
tained in various branches of science and engineering can all be rep- T =
ω
resented in the same format as above in terms of wave-frequency
/ wave-wavenumber pairs ω and β having a ratio t
ω
v≡ -1
β
Wavelength
cos(βz)
recognized as the wave-speed and specific dispersion relations 1 2π
λ=
such as: β
√
β = ω µǫ, -1
Dispersion relations
2. Acoustic waves in monoatomic gases with temperature T (K) between
and atomic mass m (kg): wavefrequency ω
r
m and
β=ω 5 , wavenumber β
3 KT
determine the
3. TEM waves in collisionless
q plasmas (ionized gases) with plasma propagation veloc-
frequency ω = Ne :
2 ity
p mǫo
ω
1q 2 v= = λf
β= ω − ωp2 . β
c
for all types of
wave motions.
2
– For any type of wave solution — TEM, acoustic, plasma wave
— once the dispersion relation is available (meaning that it has
been derived from fundamental physical laws governing the specific
wave type), wave propagation velocity is always obtained as
ω
v=
β
or, equivalently, as
λ
v= = λf
T
where
2π
λ≡ Wavelength
β
and
2π 1
T = ≡ Waveperiod.
ω f
propagatingWaveCos.eps
3
4
• Monochromatic x-polarized waves
V
E = Eo cos(ωt ∓ βz) x̂
m
can also be expressed in phasor form as
V
Ẽ = Eoe∓jβz x̂
m
such that
Re{Ẽejωt} = Eo cos(ωt ∓ βz) x̂ = E
in view of Euler’s identity.
5
Example 2: Given that
representing the sum of wave fields propagating in opposite directions, the corre-
sponding phasor
H̃ = x̂H +e−jβz − j ŷH − ejβz .
The corresponding E-field phasor is
from which
Make sure to check that all the signs make sense, and if you think you have
caught an error, let us know.
6
wave. z
y
S+
+ η z
E = −x̂ f (t − )
2 v Js = x̂f (t)
where
ω
β= and η = ηo ≈ 120π Ω
c
since the current sheet is surrounded by vacuum. Converting these into phasors,
we find
Ẽ± = −ηe∓jβz x̂ and H̃± = ∓e∓jβz ŷ.
7
z
• In the last lecture we calculated the time-average E × H and Js · E of y
the fields examined in Example 3 using a time-domain approach. The + η z
S+
E = −x̂ f (t − )
same calculations can be carried out in terms of phasors Ẽ, H̃, and J̃s 2 v Js = x̂f (t)
as follows:
1 1 1 z
hE × Hi = Re{Ẽ × H̃∗} and hJs · Ei = Re{J̃s · Ẽ∗} η z
E = −x̂ f (t + )
− H− = ŷ f (t + )
2 v
2 2 2 v
S−
where Ẽ × H̃∗ ≡ S̃ is called complex Poynting vector. Instantaneous power
V I ∗ + cc V Iej2ωt + cc
A V v(t)i(t) = + .
J̃s = 2x̂ and Ẽ±(z) = −ηe∓jβz x̂ 4 4
m m The second term has a zero time aver-
age. It follows that time-average power
in Example 3, it follows that
V I ∗ + cc 1
hv(t)i(t)i = = Re{V I ∗ }
1 W 4 2
h−Js(t) · E(0, t)i = Re{−J̃s · Ẽ∗(0)} = η ≈ 120π 2 , since
2 m
V I ∗ + cc = V I ∗ + V ∗ I = 2Re{V I ∗ }.
in conformity with the result from last lecture.
(Also see ECE 210 text.)
8
22 Phasor form of Maxwell’s equations and
damped waves in conducting media
• When the fields and the sources in Maxwell’s equations are all monochro-
matic functions of time expressed in terms of their phasors, Maxwell’s
equations can be transformed into the phasor domain. ∇·D = ρ
∇·B = 0
– In the phasor domain all ∂B
∇×E = −
∂t
∂ ∂D
→ jω ∇×H = J+ .
∂t ∂t
and all variables D, ρ, etc. are replaced by their phasors D̃, ρ̃, and
so on. With those changes Maxwell’s equations take the form shown
in the margin.
– Also in these equations it is implied that
∇ · D̃ = ρ̃
D̃ = ϵẼ ∇ · B̃ = 0
B̃ = µH̃ ∇ × Ẽ = −jω B̃
J̃ = σ Ẽ ∇ × H̃ = J̃ + jω D̃
1
We will do that, and and after that relax the requirement J̃ = 0 with
J̃ = σ Ẽ to examine how TEM waves propagate in conducting media.
– Using ∇ · Ẽ = 0
∇ · H̃ = 0
∇ × [∇ × Ẽ = −jωµH̃] ⇒ −∇2Ẽ = −jωµ∇ × H̃
∇ × Ẽ = −jωµH̃
which combines with the Ampere’s law to produce ∇ × H̃ = jωϵẼ
∇2Ẽ + ω 2µϵẼ = 0.
Ẽ = x̂Ẽx(z),
where γ is to be determined.
2
– Upon substitution into wave equation both of these lead to
(γ 2 + ω 2µϵ)Ẽx = 0,
which yields
γ 2 + ω 2 µϵ = 0 ⇒ γ 2 = −ω 2µϵ
as we already knew.
– Furthermore, using the phasor form Faraday’s law it is easy to show
that
e∓jβz
!
µ
H̃y = ± with η = .
η ϵ
3
• With ρ̃ = 0 but J̃ = σ Ẽ ̸= 0, implying non-zero conductivity σ, the
pertinent phasor form equations are as shown in the margin.
1.
2 2 ⇒⇒ "
γ = −ω µϵ = (jωµ)(jωϵ) γ = (jωµ)(σ + jωϵ)
σ ̸= 0
2. ! # #
µ jωµ ⇒⇒ jωµ
η= = η= .
ϵ jωϵ σ ̸= 0 σ + jωϵ
Note that the modified γ and η satisfy
γ γη
γη = jωµ and = σ + jωϵ µ =
η jω
γ
leading to useful relations shown in the margin (assuming real valued σ = Re{ }
η
σ and ϵ). 1 γ
ϵ = Im{ }
ω η
4
(a) Damped wave snapshot at t=0
• In terms of γ and η above, we can express an x-polarized plane wave together with exponential envelope
Eo ∓γz
Ẽ = x̂Eoe∓γz and H̃ = ±ŷ e z
η
e−αz cos(ωt − βz)|t=0
where Eo is an arbitrary complex constant (complex wave amplitude).
• In expanded forms γ and η appear as: (b) Snaphot at t>0, with t=0 waveform
for comparison
"
γ = (jωµ)(σ + jωϵ) ≡ α+jβ, so that α = Re{γ} and β = Im{γ},
and z
# # #
e−αz cos(ωt − βz)
jωµ jωµ jωµ
η= ≡ |η|ejτ so that |η| = | | and τ = ∠ .
σ + jωϵ σ + jωϵ σ + jωϵ
– β appears within cosine
argument and deter-
1. In the special case of a perfect dielectric with σ = 0, we find mines the wavelength
! 2π
√ µ λ=
β
γ = jω µϵ ≡ jβ and η = ,
ϵ and propagation speed
and, therefore, ω
vp = .
β
5
2. Another case of imperfect dielectric (or “lousy” conductor) occurs
when σ is not zero, but it is so small that are justified in using
! ! e−αz
ωµσ " ωµ
α≈β≈ = πf µσ while |η| = and τ = ∠η = 45o.
2 σ z
7
(a) Damped wave snapshot at t=0
• Wavelength together with exponential envelope
2π vp e−αz
λ= =
β f
now depends on frequency f via both the numerator and the denomi- z
nator, and measures twice the distance between successive nodes of the
e−αz cos(ωt − βz)|t=0
waveform.
• Penetration depth (also called skin depth if very small) (b) Snaphot at t>0, with t=0 waveform
for comparison
1 1
δ≡ = "
α Re{ (jωµ)(σ + jωϵ)}
is the distance for the field strength to be reduced by e−1 factor in its z
8
23 Imperfect dielectrics, good conductors
2π 1
Condition β α |η| τ λ= β
δ= α
Perfect √ pµ 2π
σ=0 ω ǫµ 0 ǫ 0 √
ω ǫµ ∞
dielectric
Imperfect σ √ pµ pµ q
ωǫ
≪1 ∼ ω ǫµ β 12 ωǫ
σ
= σ
2 ǫ
∼ ǫ
∼ σ
2ωǫ
∼ 2π
√
ω ǫµ
2
σ
ǫ
µ
dielectric
Good σ
√ √ p ωµ
≫1 ∼ πf µσ ∼ πf µσ 45◦ ∼ √ 2π ∼ √ 1
conductor ωǫ σ πf µσ πf µσ
Perfect
σ=∞ ∞ ∞ 0 - 0 0
conductor
• The table above summarizes TEM wave parameters in homogeneous
conducting media where the propagation velocity x-polarized phasor
ω
vp = Ẽ = x̂Eoe∓αz e∓jβz
β
(note that it can be frequency dependent) and field phasors can be accompanied by
expressed in formats similar to that shown in the margin, keeping in
mind that propagation direction coincides with vector H̃ = ±ŷ Eηo e∓αz e∓jβz .
S̃ ≡ Ẽ × H̃∗
such that
1
hSi = hE × Hi = Re{S̃}
2
is the average energy flux per unit area (time-average Poynting vector).
1
Example 1: Consider the plane TEM wave
V
Ẽ = ŷ2e−αz e−jβz
,
m
in an imperfect dielectric. Determine H̃ and time-average Poynting vector
hSi. Compute hSi at z = 0 and z = 10 m, if ǫ = 4ǫo , µ = µo , σ = 10−3 S/m,
and ω = 2π · 109 rad/s
Solution: Using right hand rule, so that E × H points in propagation direction ẑ, we
find that
2 2 A
H̃ = −x̂ e−αz e−jβz ≈ −x̂ p e−αz e−jβz e−jτ
η µ/ǫ m
pµ
using |η| = ǫ from the table above for a perfect dielectric.
2
Hence, at z = 0,
2 2 1 W
hSi = ẑ p cos τ ≈ ẑ = ẑ ,
µ/ǫ 60π 30π m2
whereas, at z = 10 m,
2 2 −6π/10 0.15 W
e−2·30π·10 ·10 cos τ ≈ ẑ
−3
hSi = ẑ p e ≈ ẑ .
µ/ǫ 60π 30π m2
• Note that in above example power transmitted per unit area has dropped
to 15% of its value upon propapagating over a relatively short distance
of 10 m.
– In the physical terms, the lost power of the wave is gained by the
propagation medium in the form of heat — average Joule heating ⇐This is what we
hJ · Ei in the medium will be positive and account for the loss of want to happen in a
the wave power (as seen in a HW problem). microwave oven.
3
Example 2: Repeat Example 1 for ω = 2π · 103 rad/s and σ = 4 S/m (sea water) in
which case the propagation medium becomes a good conductor.
Solution: Using right hand rule, so that E × H points in propagation direction ẑ, we
have
2 2 A
H̃ = −x̂ e−αz e−jβz ≈ −x̂ e−αz e−jβz e−jτ
η |η| m
as well as
1 1 2
hSi = Re{Ẽ × H̃∗ } = Re{ŷ2e−αz e−jβz × (−x̂ e−αz e−jβz e−jτ )∗}
2 2 |η|
1 2 2
= − Re{ŷ2e−αz × x̂ e−αz ejτ } = ẑ e−2αz cos τ.
2 |η| |η|
4
whereas, at z = 10 m,
100 −2· π ·10 100 W
hSi = ẑ e 25 ≈ ẑ 0.081 2 .
π π m
The next example identifies the penetration depth in sea water at 1 kHz.
5
Example 3: What is the penetration depth δ = α−1 in a medium with σ = 4 S/m,
ǫ = 81ǫo , and µ = µo for ω = 2π · 103 rad/s.
6
24 Signal transmission, circular polarization
Since in perfect dielectrics the propagation velocity vp = v and the intrinsic
impedance η are frequency independent (i.e., propagation is non-dispersive),
d’Alembert plane wave solutions of the form
z f (t − zv )
E = x̂f (t − ) and H = ŷ
v η
are valid in such media.
6
m(t)
• Consider a waveform 4
-6
2π
– ω is some specific frequency having a corresponding period T = ω, 6
m(t) cos(ωt)
– m(t) is some arbitrary signal (e.g., a voice signal, a message) 4
-4 -2 2 4 6
-2
In that case, -4
t
-6
– f (t) specified above can be called narrowband AM, and
– ω the carrier frequency of modulating cosine of the message
signal m(t).
1
The corresponding x-polarized wave fields propagating in z direction
can then be represented as Field 1
z m(t − vz )
E = m(t − ) cos(ωt − βz)x̂ and H = cos(ωt − βz)ŷ
v η
√
where β = ω µǫ as usual1.
• By contrast, Field 2
z
E = m(t − ) cos(ωt − βz)ŷ
v
represents an AM wave field with a y-polarized carrier, and so does Field 3
z
E = m(t − ) sin(ωt − βz)ŷ
v
but with a carrier that has been time-delayed by a quarter period.
2
CIRCULAR POLARIZATION:
which has a circular polarized carrier. Since this is just a superpo- Field vector rotates instead
sition of two d’Alembert waves, the accompanying H is easily found to of oscillating.
z
H = m(t − )[cos(ωt − βz)ŷ − sin(ωt − βz)x̂]/η. cos(ωt − βz)x̂ + sin(ωt − βz)ŷ
v z
y
or Left-circular
3
x̂ − j ŷ
• Also, y
z
cos(ωt)x̂ + sin(ωt)ŷ ⇔ x̂ − j ŷ
RIGHT CIRCULAR
rotates in the direction that your right-hand fingers curl when x-comp leads y-comp because of -j
vector
cos(ωt)x̂ − sin(ωt)ŷ ⇔ x̂ + j ŷ t=0
E x
rotates in the direction that your left-hand fingers curl when the t>0
thumb is directed in propagation direction z.
LEFT CIRCULAR
4
Example 1: A circular polarized wave field vector is given as
Ẽ = (ẑ + j ŷ)ejβx .
Determine the propagation and rotation directions of the field vector as well as
its helicity.
Solution: The propagation direction is −x since the exponent in ejβx lacks a minus
sign.
When I direct my right thumb in −x direction, my fingers curl from z- toward y-axis,
LEFT CIRCULAR
which is curling in the wrong direction. Hence this wave is not right-circular! It
(i) z-comp lags y-comp because of +j.
is left-circular.
(ii)E-vector rotates from y towards z.
5
polarizations.
• EXAMPLE: Right- and left circular waves propagating in z directions
are weighted superpositions of orthogonal x- and y-polarized fields
as in (expressed in terms of phasors): basis functions Circulars
in terms of
x̂e−jβz and ŷe−jβz
linears
superpose to form right- and left-circular waves
(x̂ − j ŷ)e−jβz and (x̂ + j ŷ)e−jβz
using the weights
1, −j and 1, j
respectively.
6
• It can be argued that right- and left-circular wave pair forms an in-
trinsically more fundamental set of basis functions than, say, x̂- and
ŷ-polarized waves, because while the selection of which direction is x
and which direction is y can be arbitrary, there is no arbitrariness in
how helicity is assigned to circular polarized modes propagating in a
given direction2 .
7
Note that this figure
only shows one linear
component of the sur-
Example 2: On z = 0 plane we have a time-varying surface current density
face current on
A z = 0 plane. One linear
Js (t) = m(t)[cos(ωt)x̂ + sin(ωt)ŷ]
m
component causes a lin-
with a carrier frequency of ω. Determine the radiated wave fields E± and the ear polarized radiation.
polarization (and the helicity if appropriate) of the carrier.
An orthogonal pair of
Solution: We have already learned that a surface current Js (t) on z = 0 plane will linear components will
produce TEM wave fields
η z conspire to radiate a cir-
E± = − Js (t ∓ )
2 v cular polarized wave as
in surrounding regions. With the given Js (t) , this implies in Example 2 when they
η z V are 90◦ out of phase.
E± = − m(t ∓ )[cos(ωt ∓ βz)x̂ + sin(ωt ∓ βz)ŷ] , z
2 v m
y
which has a circular-polarized carrier
S+
+ η z
cos(ωt ∓ βz)x̂ + sin(ωt ∓ βz)ŷ E = −x̂ f (t − )
2 v Js = x̂f (t)
cos(ωt)x̂ + sin(ωt)ŷ. 1 z
η z H− = ŷ f (t + )
E = −x̂ f (t + )
−
2 v
This vector rotates from x- toward y-axis, and therefore the carrier of E+ is 2 v
S−
right-circular and the carrier of E− is left-circular.
8
Example 3: In Example 2, what is the average power density of the circular polarized
carrier signal
V
Ec = cos(ωt − βz)x̂ + sin(ωt − βz)ŷ
m
in the region z > 0, assumed to be vacuum?
Make sure you check and follow all the sign changes that take place in
Example 3.
9
25 Wave reflection and transmission
In this lecture we will examine the phenomenon of plane-wave reflections at
an interface separating two homogeneous regions where Maxwell’s equations
allow for traveling TEM wave solutions. The solutions will also need to
satisfy the boundary condition equations repeated in the margin. We will
consider a propagation scenario in which (see margin): n̂ · (D+ − D−) = ρs
n̂ · (B+ − B−) = 0
1. Region 1 where z < 0 is occupied by a perfect dielectric with medium n̂ × (E+ − E−) = 0
parameters µ1, ǫ1, and σ1 = 0,
n̂ × (H+ − H− ) = Js
2. Region 2 where z > 0 is homogeneous with medium parameters µ2, ǫ2,
y
and σ2,
Hi Ht
z
• In Region 1 we envision an incident plane-wave with linear-polarized field phasors
Er
Eo −jβ1 z Hr
Ẽi = x̂Eoe−jβ1 z and H̃i = ŷ e ,
η1
where
1
Fields above satisfy Maxwell’s equations in Region 1, but if there were no
other fields in Regions 1 and 2 boundary condition equations requiring
continuous tangential E and H at the z = 0 interface would be violated.
In order to comply with the boundary condition equations we postulate
a set of reflected and transmitted wave fields in Regions 1 and 2 as follows: Incident:
• In Region 1 we postulate a reflected plane-wave with linear-polarized field phasors
jβ1 z ΓEo jβ1 z Ẽi = x̂Eoe−jβ1 z ,
Ẽr = x̂ΓEoe and H̃r = −ŷ e
η1 Eo −jβ1 z
including an unknown Γ that we will refer to as reflection coefficient.
H̃i = ŷ e ,
η1
– Note that the reflected wave propagates in −z direction (direction of H̃r and
the exponential terms have been adjusted accordingly). Reflected:
2
• To determine the unknowns Γ and τ we enforce the following boundary
conditions at z = 0 where the fields simplify as shown in the margin: Incident at z = 0:
3
Special cases: Reflection coeff.:
η2 − η1
1. Region 2 is a perfect conductor with σ2 → ∞: In that case Γ= ,
η2 + η1
η2 → 0, and consequently
4
y
Hi Ht
Ei Et
Example 1: An plane-wave in vacuum,
x Region 2
√ Region 1
Ẽi = x̂ 120πe−jβ1 z m
V
,
z
9
is incident at z = 0 on a dielectric medium with µ = µo and ǫ = ǫ.
4 o
Determine Er
the average Poynting vectors hSi i, hSr i, and hSt i of the incident, reflected, and Hr
transmitted fields.
5
The transmitted wave, likewise, has the field phasors
4 √ 4 √
Ẽt = x̂ 120πe−jβ2 z and H̃t = 2 ŷ 120πe−jβ2z
5 5 3 ηo
and
1 1 4 3 120π 1 4 3 W
hSt i = Re{Ẽt × H̃∗t } = ẑ ( )2 ≈ ẑ ( )2 .
2 2 5 2 ηo 2 5 2 m2
As for the incident wave
√ 1 √
Ẽi = x̂ 120πe−jβ1 z and H̃i = ŷ 120πe−jβ1 z
ηo
and
1 1 120π 1 W
hSi i = Re{Ẽi × H̃∗i } = ẑ ≈ ẑ .
2 2 ηo 2 m2
Note: We have
1 1 16 3 1 1 24 1
|hSr i| + |hSt i| = ( + ) = ( + ) = = |hSi i|
2 25 25 2 2 25 25 2
in compliance with energy conservation (as expected) — energy flux per unit
area of the transmitted and reflected waves add up the that of the
incident wave!
6
26 Standing waves, radiation pressure y
Hi Ht = 0
We continue in this lecture with our studies of wave reflection and transmis-
Ei Et = 0
sion at a plane boundary between two homogeneous media.
Region 1
x Region 2
• In case of total reflection from a perfectly conducting mirror placed at
z
z = 0 surface, Γ = −1, and the incident and reflected waves in z < 0
Er
region combine to produce standing waves of electric and magnetic Hr
field:
to the time-domain (by multiplying with ejωt and taking the real
time as usual) we obtain:
2Eo z
E(z, t) = x̂2Eo sin(β1z) sin(ωt) and H(z, t) = ŷ cos(β1z) cos(ωt);
η1
– Standing waves carry no net energy, that is, with standing wave
fields we have
z
hE × Hi = 0,
because of the cancellation of the power transported by its travel-
ing wave components in opposite directions. λ
Verification: Using the phasors 2
Note: Nulls in Ex and Hy are
2Eo separated by half wavelength.
Ẽ = −j x̂2Eo sin(β1z) and H̃ = ŷ cos(β1 z),
η1 Adjacent nulls of Ex and Hy
are separated by quarter
we have wavelength.
2
– Note that E(0, t) = 0 on z = 0 surface satisfying the tangential Ex(z, t) ∝ sin(βz) sin(ωt)
2Eo A,
Js = x̂ cos(ωt) m
η1
z
satisfying
−ẑ × H(0, t) = Js.
λ
Js on mirror surfaces is really a convenient idealization of volume cur-
2
rents flowing in thin layers — just a few skin depths — near good- Note: Nulls in Ex and Hy are
separated by half wavelength.
conductor surfaces (real-life mirrors are good but not perfect conduc-
tors!). Radiation due to Js in effect causes the reflected wave and also Adjacent nulls of Ex and Hy
are separated by quarter
cancels out the incident wave field in z > 0. wavelength.
the limiting case of a perfect conductor is naturally reached. Conductor shorts Ex on its
surface where a current flows.
3
• Going back to the partial reflection case, consider the transmitted fields
Ẽt = x̂τ Eoe−γ2 z
Ẽt and H̃t in Region 2 shown in the margin. Also shown in the margin
are the phasors for current density Jt and magnetic flux density Bt. τ Eo −γ2 z
H̃t = ŷ e
η2
In the box below we integrate the volumetric current density Jt of a good
−γ2 z
conductor from z = 0 to ∞ and find out that this “depth integral” matches J̃t = σ2Ẽt = x̂σ2τ Eoe
the surface current density found above for the case of perfect conductor. In B̃ = µ H̃ = ŷ µ2τ Eo e−γ2 z
t 2 t
this calculation we assume that Region 1 is vacuum, and also take µ2 = µo: η2
and therefore
√
2η2 2η2 2σ2η2 2 jωµ2 σ2 2γ2
τ= ≈ and σ2 τ ≈ = = .
ηo + η2 ηo ηo ηo ηo
The depth integral of the volumetric current density in Region 2, that is, the effective
surface current of the region is then
Z ∞ Z ∞
1 2Eo
J̃t dz = x̂ σ2τ Eo e−γ2 z dz = x̂Eo (σ2τ ) = x̂
0 0 γ2 ηo
in phasor form, matching the phasor of the time-domain result from above, namely
2Eo A
Js = x̂ cos(ωt) m
ηo
representing the surface current on an idealized perfect conductor surface.
4
Surface resistance: Let J̃s stand for the effective surface current of a good con-
ductor with a propagation constant
p
γ ≈ jωµσ = α + jβ = α + jα
5
• Let’s finally calculate the magnetic component of the Lorentz force on
charge carriers of a good conductor due to the penetrating fields: Radiation pressure propor-
tional to
Radiation pressure: If there are N free charge carriers per unit volume inside a hSi i/c
reflecting mirror, then
N F = N qv × Bt = Jt × Bt shows that plane-TEM waves
not only carry and transport
will be the force per unit volume of the mirror, expressed in terms of current density energy, but also momentum.
Jt = N qv and the magnetic flux density Bt .
Its integral over all z can be interpreted as the total force per unit area of the
mirror, Z ∞ TEM waves not only
Prad = Jt × Bt dz, heat, but also push!
0
having a magnitude known as radiation pressure of the reflecting wave. This is a
time-varying quantity, with a time-average
Z ∞ ((It can also be shown that the momen-
1 tum density of the wave is
hPrad i = Re{J̃ × B̃∗}dz
0Z 2
∞ hSi i/c2 N.s/m3
1 µ2τ Eo −2α2z
= ẑ Re{(σ2τ Eo )( )}e dz
0 2 η2 and (spin) angular momentum density
6
27 Guided TEM waves on TL systems
• An x̂ polarized plane TEM wave propagating in z direction is depicted in the mar-
x
gin.
E z
– A pair of conducting plates placed at x = 0 and x = d would not perturb E×H
the fields except that charge and current density variations would be induced H
Unguided uniform plane
wave propagation in a
on plate surfaces at x = 0 and x = d (on both sides) to satisfy Maxwell’s homogeneous medium
y
boundary condition equations.
• If charge and currents were confined only to interior surfaces of the plates facing one
x
another, fields E and H accompanying them would be restricted to the region in d
I
between the plates, constituting what we would call guided waves. E z
E×H Plate 1
– Such a guided wave field confined to the region between the plates will sat- W
H
I
Plate 2
isfy Maxwell’s equations including a minor fringing component that can be
y
neglected when the plate width W is much larger than plate separation d.
– Guided waves produce wavelike surface charge and current variations on plate
surfaces.
– Conversely, wavelike charge and current variations on plate surfaces would
produce guided wave fields.
1
x
How such excitations propagate away from their “source points” on TL systems will d
be our main subject of study for the rest of the semester. I
E z
E×H Plate 1
• In a parallel-plate TL we ignore any fringing fields and assume that H
I
W Plate 2
TEM wave fields y
E = x̂Ex (z, t) and H = ŷHy (z, t)
occupy the region between the plates. For these fields uniform in x and
y, Faraday’s and Ampere’s laws reduce to scalar expressions
∂H ∂Ex ∂Hy
∇ × E = −µ ⇒ = −µ
∂t ∂z ∂t
and
∂E ∂Hy ∂Ex
∇ × H = σE + ǫ ⇒ − = σEx + ǫ .
∂t ∂z ∂t
• Now, multiply both equations by d and let Note that voltage drop
Z 1
V ≡ Ex d voltage drop from plate 2 to plate 1 V = E · dl = Ex d
2
2
• We can re-write these equations as
∂V ∂I ∂I ∂V
− =L and − =C + GV
∂z ∂t ∂z ∂t
utilizing
d W W
L=µ , C=ǫ , G=σ
W d d
appropriate for the parallel-plate TL — we recognize these parameters
as inductance, capacitance, and conductance of the parallel plate TL. Telegrapher’s
equations:
– In the equations above the GV term accounts for Ohmic losses of
wave fields having to do with currents leaking between the wires
(plates) of the TL.
∂V ∂I
– Another possible loss term that we have not picked up — because − = L
∂z ∂t
we assumed infinite conducting plates — is a missing RI term in ∂I ∂V
− = C
the right-hand-side of the first equation. ∂z ∂t
Rather than correcting for that at this stage, we will drop the GV term
from the second equation, and focus our attention for a while (until
the last day of the semester, in fact) on ideal lossless transmission lines
governed by the equations shown in the margin — they are known as
known as telegrapher’s equations1.
1
Telegrapher’s equations were first compiled by Oliver Heaviside (of close-up method, unit-step, and
countless other contributions) in 1880’s. Telegraphy was being used worldwide by 1850’s as a means of
rapid communications.
3
∂
• Except for − ∂z on the left, the telegrapher’s equations look like the
V − I relations of inductors and capacitors (which is the best way of
remembering them). Telegrapher’s
equations:
• The equations can be readily combined to obtain a 1D scalar wave
equation
∂ 2V ∂ 2V
= LC 2 . ∂V ∂I
∂z 2 ∂t − = L
In analogy to ∂z ∂t
∂ 2 Ex ∂ 2 Ex ∂I ∂V
= µǫ 2 , − = C
∂z 2 ∂t ∂z ∂t
the wave equation for V has d’Alembert wave solutions where
z 1 1 µ
V (z, t) = f (t ∓ ) where v ≡ √ =√ . C = ǫGF, L =
GF
,
v LC µǫ
with “geometrical factor”
• In that case the second telegrapher’s equation demands
W
GF = parallel-plate
∂I ∂V z d
− =C = Cf ′ (t ∓ )
∂z ∂t v =
2π
coax
ln ab
implying an anti-derivative π
= twin-lead
cosh−1 2aD
z f (t ∓ vz )
I(z, t) = ±Cv f (t ∓ ) = ±
v Zo
with √ r r
1 LC L 1 µ
Zo ≡ = = = .
Cv C C GF ǫ
4
• In summary, d’Alembert wave solutions of telegrapher’s equations are
z f (t ∓ vz )
V (z, t) = f (t ∓ ) and I(z, t) = ±
v Zo
with a propagation speed
1 1
v=√ =√ a) t <0
LC µǫ ++ Wire 2
that equals the wave speed of the associated electric and magnetic fields, 3V
Wire 1
and voltage-to-current ratio representing a characteristic impedance --
r 0 z
r
L 1 µ
Zo = = .
C GF ǫ
b) t = 0+
++++I Wire 2
Telegrapher’s equations and their d’Alembert solutions provide us with a
3V V
“handle” on the following physics: Wire 1
----
I
• Suppose that + and - terminals of a 3 V battery makes contact with
the terminals of a charge neutral TL at t = 0 as depicted in the margin.
We will assume that V (z, t) = I(z, t) = 0 on the TL for t < 0. c) t > 0+
+++++++I Wire 2
As soon as contact is made between the terminals of the battery and the
3V V z = vt
TL, the excess + and - charges on battery terminals will “spill onto” the TL Wire 1
5
a) t <0
– electrons move from the - terminal of the battery onto the bottom ++ Wire 2
wire of the TL,
3V
– replenished by an equal amount of electrons moving from the top Wire 1
--
wire into the battery via its + terminal, 0 z
I
the rest of the TL!!!
– because, for instance, electrons at some z > 0 on the top wire will
start moving towards the battery terminal only after the neighbor-
ing electrons at z − deplete the region leaving some excess positive
charge.
6
analogy with ripples spreading out on a pond surface when perturbed
by a falling pebble.
7
Evidently, a short section ∆z of the TL has an equivalent T-network with I(z) I(z + ∆z)
1. a series inductance ∆zL carrying a current I(z, t) ≈ I(z + ∆z, t), and Wire 2 + ∆zL
V (z) ∆zC
2. a shunt capacitance ∆zC carrying a voltage V (z + ∆z, t) ≈ V (z, t)
Wire 1 -
as shown in the margin. z z + ∆z
This lumped-circuit equivalent is only accurate for ∆z so small that
I(z, t) ≈ I(z + ∆z, t) and V (z + ∆z, t) ≈ V (z, t)
are both true, which is of course possible only if ∆z ≪ λ, λ being the shortest
wavelength in I(z, t) ∝ H(z, t) and V (z, t) ∝ E(z, t) waveforms. TL’s can also support non-TEM
a+b < λ
π (coax). At high frequen-
the familiar formula from circuit theory. cies when these modes cannot be
8
28 Distributed circuits and bounce diagrams
Last lecture we learned that voltage and current variations on TL’s are gov-
erned by telegrapher’s equations and their d’Alembert solutions — the latter
can be expressed as
z z f (t − vz ) g(t + zv )
V (z, t) = f (t − ) + g(t + ) and I(z, t) = −
v v Zo Zo
in terms of r
1 L Source
ckt
Transmission line
v=√ and Zo = Rg I(z, t) Wire 2
IL
LC C + + Load
+ V (z, t) Zo VL RL
and functions f (t) and g(t) corresponding to signal waveforms propagated in -
fi (t) I(z, t) - Wire 1 -
+z and −z directions, respectively. 0 l z
1
– We want to determine voltage and current signals V (z, t) and
I(z, t) on the TL and the load RL for time t > 0 in terms of
source signal fi(t) assuming that fi (t) = 0 for t < 0.
l RL − Zo l 2l
g(t + ) = f (t − ) ⇒ g(t) = ΓLf (t − )
v |RL {z
+ Zo} v v
ΓL
where
RL − Zo
ΓL =
RL + Zo
is the load reflection coefficient in the TL circuit. We can re-write
the d’Alembert solution for V (z, t) and I(z, t) in terms of only f (t) as
z z 2l f (t − zv ) ΓLf (t + zv − 2lv )
V (z, t) = f (t− )+ΓLf (t+ − ) and I(z, t) = − .
v v v Zo Zo
• Assuming that fi(t) = 0 = f (t) for t < 0, we can relate f (t) to fi (t) in
t > 0 interval using the KVL equation at z = 0 that states
fi(t) = Rg I(0, t) + V (0, t),
2
which is, using V (z, t) and I(z, t) at z = 0,
f (t) Zo
fi (t) = Rg + f (t) ⇒ f (t) = fi (t) Source Transmission line
Zo Rg + Zo ckt
Rg I(z, t) IL
| {z } Wire 2
+ + Load
τg + V (z, t) Zo VL RL
-
fi (t) I(z, t) - Wire 1 -
where 0 l z
Zo
τg =
Rg + Zo
is the injection coefficient of the TL circuit1 .
2l
• Thus, for the epoch 0 < t < v, we have the voltage and current
solutions
z z 2l τg fi(t − vz ) ΓLτg fi (t + vz − 2lv )
V (z, t) = τg fi(t− )+ΓL τg fi (t+ − ) and I(z, t) = −
v v v Zo Zo
on the line.
1
Note how f (t) appears to be related to fi (t) according to a voltage division rule with Zo representing
the resistance across which voltage f (t) is measured.
3
– So far fi(t) function is arbitrary and the above results would also
be valid for fi(t) = δ(t), Dirac’s impulse, in which case
| Zo {z Zo v} z
} | {z τg
l
I(0, t) V (0, t) τ g ΓL v
2l
v τ g ΓL Γg
which can be re-arranged as 3l
τg Γ2L Γg v
Rg Rg 2l 4l
δ(t) = (1 + )f (t) + (1 − )ΓL f (t − ), v τg Γ2L Γ2g
Zo Zo | {z v }
τg Γ3L Γ2g
2l
τg δ(t − )
v t
4
– Hence, solving this for f (t), we find, for this epoch,
Rg − Zo 2l
f (t) = τg δ(t) + ΓLτg δ(t − ),
R +Z v
| g {z o}
Γg
where
Rg − Zo
Γg =
Rg + Zo
is the source reflection coefficient of the TL circuit.
2l 4l
– Substituting f (t) for the epoch v
< t < v
within voltage and
current formulae
z z 2l f (t − vz ) ΓLf (t + vz − 2lv )
V (z, t) = f (t− )+ΓL f (t+ − ) and I(z, t) = −
v v v Zo Zo
we obtain the “extended” voltage and current impulse response
functions
z z 2l z 2l z 4l
V (z, t) = τg δ(t− )+ΓLτg δ(t+ − )+Γg ΓLτg δ(t− − )+Γg Γ2Lτg δ(t+ − )
v v v v v v v
and
z z 2l z 2l z 4l
I(z, t) = Zo−1 [τg δ(t− )−ΓLτg δ(t+ − )+Γg ΓL τg δ(t− − )−Γg Γ2L τg δ(t+ − )]
v v v v v v v
respectively.
5
◦ At this point the algebra is pretty messy, but a straightforward Rg I(z, t) IL
puts. l
τ g ΓL v
– The horizontal axis represents position z of the traveling impulses 2l
v τ g ΓL Γg
while time t is represented by a downward pointing axis. 3l
τg Γ2L Γg v
– The first slanted line on the top of the diagram, representing the 4l
v τg Γ2L Γ2g
traveling impulse
z
τg δ(t − ), τg Γ3L Γ2g
v
(first term of hz (t) = V (z, t)) is “reflected” at time t = vℓ from load t
6
to become a forward propagating impulse
z 2ℓ
τg ΓLΓg δ(t −
− )
v v
represented by the third line of the diagram.
◦ Reflection at Rg is in effect the same physical process as re-
flection at RL and therefore its coefficient Γg is identical with Rg IL
I(z, t)
ΓL except for the replacement of RL by Rg . + +
+ V (z, t) Zo VL RL
– The bounce diagram is advanced in time with further reflections δ(t)
-
I(z, t) - -
occurring at both ends. 0 l z
7
• Although these series formulae look daunting, only the lower order
terms usually matter — that is true because |ΓL| ≤ 1 and |Γg | ≤ 1 and Rg I(z, t) IL
• We typically rely on the bounce diagram technique more so than the Bounce diagram
l
z
series expressions developed above. This will be illustrated by several τg
l
examples in the next lecture. τ g ΓL v
2l
v τ g ΓL Γg
– The main idea is to combine delayed versions of the circuit input 3l
fi(t) with the impulse weights indicated on the bounce diagram, τg Γ2L Γg v
4l
since, in general, the convolution δ(t − Tz ) ∗ fi (t) = fi (t − Tz ) for v τg Γ2L Γ2g
8
29 Bounce diagrams
• Last lecture we obtained the implulse-response functions Source matched to
line:
z z 2l Zo I(z, t)
V (z, t) = τg [δ(t − ) + ΓLδ(t + − )]
v v v +
+ V (z, t) Zo RL
-
and f (t) I(z, t) -
τg z z 2l 0 l z
I(z, t) = [δ(t − ) − ΓLδ(t + − )]
Zo v v v Bounce diagram l
for the voltage and current in the TL circuit shown in the margin where z
τg
l
the source is matched to the line so that τg = 21 — circuit response with τg Γ L v
an arbitrary input f (t) is obtained by convolving these with f (t) (as 2l
v
shown in Example 1 in last lecture).
1
– the next line down is the depiction of backward propagating im-
pulse τg ΓLδ(t + zv − 2lv ).
Zo Bounce diagram
Rg and τg = , l
Rg + Zo τg
z
l
2l
for which the bounce diagram is not terminated at t = because the
v τ g ΓL v
2l
backward propagating impulse on the line arriving at z = 0 at time v τ g ΓL Γg
l
∞
X z 2l τ g ΓL v
V (z, t) = τg (ΓLΓg )n δ(t − −n ) 2l
n=0
v v v τ g ΓL Γg
∞ 3l
X
n z 2l v
+τg ΓL (ΓLΓg ) δ(t + − (n + 1) ). τg Γ2L Γg
n=0
v v 4l
v τg Γ2L Γ2g
– It can also be shown that the first term of V (z, t) above is derived
from the formal solution of the equation
2l
V +(t) = τg δ(t) + ΓLΓg V +(t − )
v
which is obtained from
δ(t) − V (0, t)
I(0, t) =
Rg
enforced at z = 0. We have effectively by-passed such a formal
approach to the problem by using the bounce diagram technique.
3
• These awful series formulae above are hardly needed in most appli-
cations when only the first few terms of the series are sufficient for
reasonably accurate results (like in the next example).
4
z = 1200 m
I(z, t)
+
+
- V (z, t) RL = 100 Ω
u(t) I(z, t) - Zo = 50 Ω
Example 1: Consider a TL circuit where Zo = 50 Ω, v = c, l = 2400 m, Rg = 0, and 0 l z
RL = 100 Ω. Determine and plot V (1200, t) if f (t) = u(t).
Bounce diagram
l
z
Solution: For this circuit Γg = −1 1 ΓL =
1
3 l
Zo Rg − Zo RL − Zo 1 1
v
τg = = 1, Γg = = −1, and ΓL = = . 3
Rg + Zo Rg + Zo RL + Zo 3 2l 1
v −
3
Also, the transit time across the TL is 1
3l
− v
9
l 2400 m 4l 1
= = 8 µs. v
v 300 × 106 m/s 9
1
From the bounce diagram shown in the margin, the impulse response for z = 1200 27
Replacing the δ(t) in this expression with the unit-step u(t), the specified source 1.5
5
Bounce diagram
l
1 z
• Note that as t → ∞, V (1200, t) → 1 V in Example 1, as if DC condi- Γg = −1 1 ΓL =
3 l
1
tions prevail and the TL becomes a pair of wires in the lumped circuit 3
v
2l
sense. v −
1
3
3l
1
− v
– DC steady-state corresponds to ω = 0 and signal wavelength λ → 4l
9
1
∞. In that limit l ≪ λ is always valid and TL can be treated like v
9
1
an ordinary lumped circuit. 27
1
– Of course this simplification can only occur with f (t) ∝ u(t), or t
−
27
example. 1.5
1.0
0.5
t us
0 10 20 30 40 50 60
6
Example 2: In the TL circuit described in Example 1, determine V (z, t) and I(z, t)
for a new source signal f (t) = rect( Tt )+2rect( t−T
T ), T = 1 µs. Plot V (z, t) versus
z at t = 3 µs and t = 11 µs.
n=0
3 c 3 n=0 3 c
0.5
z
where c is to be entered in µs units. Also,
t us
-2 -1 1 2 3 4
∞ ∞
1 X 1 n z 1 1X 1 n z VHz,3L
I(z, t) = (− ) f (t − − n16) − (− ) f (t + − (n + 1)16). 2.0
50 n=0 3 c 50 3 n=0 3 c
1.5
∞ ∞
X 1 z 1X 1 n z
V (z, 3) = (− )nf (3 − − n16) + (− ) f (3 + − (n + 1)16), 0.5
n=0
3 c 3 n=0 3 c
zm
500 1000 1500 2000
which is plotted in the margin using f (t) = rect(t) + 2rect(t − 1). Likewise, at VHz,11L
2.0
t = 11 µs,
1.5
∞ ∞
X 1 n z 1X 1 n z
V (z, 11) = (− ) f (11 − − n16) + (− ) f (11 + − (n + 1)16). 1.0
n=0
3 c 3 n=0 3 c
0.5
zm
0 500 1000 1500 2000
7
30 Multi-line circuits Rg = 2Z1
Z1 Z2 = 2Z1
Load
+ RL = Z2
-
v1 v2 = 2v1
f (t)
• In this lecture we will extend the bounce diagram technique to solve 0 l/2 l
margin. t
1
Rg = 2Z1
– Zj is the characteristic impedance of the line of the incident pulse, Load
Z1 Z2 = 2Z1
while +
-
f (t)
v1 v2 = 2v1
RL = Z2
0 l/2 l
– Zk is the impedance of the cascaded line into which the transmitted
pulse is injected. τg = Γg = Γ12 =
1
3
τ12 =
4
3
ΓL = 0
0 l z
Verification:
1
3
4
9
1
– Let 9
81
– let
Vj+τjk and Vj+τjk /Zk
denote the voltage and current on line Zk adjacent to line Zj .
This notation identifies Γjk and τjk as reflection and transmission coef-
ficients at the junction.
– Taking
Vj+(1 + Γjk ) = Vj+τjk
and
Vj+(1 − Γjk )/Zj = Vj+τjk /Zk
in order to enforce voltage and current continuity, we can solve for
Γjk and τjk expressions stated above.
2
Rg = 2Z1
Load
+ Z1 Z2 = 2Z1
- RL = Z2
v1 v2 = 2v1
f (t)
0 l/2 l
Example 1: In the circuit shown in the margin with two TL segments, line 2 has twice
the characteristic impedance and propagation velocity of line 1, i.e., 1 4
ΓL = 0
τg = Γg = Γ12 = τ12 =
0 3 3 l z
Z2 = 2Z1 and v2 = 2v1. 1
3
Solution: We have
L2 L1 1
Z2 = 2Z1 ⇒ =4 27
C2 C1 4
81
and 1
81
1 1
v2 = 2v1 ⇒ =4 . t
L2 C2 L1C1
The product of the two equations gives
1 1 1
= 16 ⇒ C2 = C1 ,
C22 C12 4
while their ratio leads to
L2 = L1 .
3
Rg = 2Z1
Load
+ Z1 Z2 = 2Z1
- RL = Z2
v1 v2 = 2v1
f (t)
Example 2: In the circuit of Example 1, determine V (z, t) and I(z, t) if 0 l/2 l
Solution: From the bounce diagram we infer the following impulse-response for the 3
4
voltage variable: 1
9
9
∞
1 X 1 2n z l 1 z l
V (z, t) = ( ) [δ(t − − n ) + δ(t + − (n + 1) )] 1
3 n=0 3 v1 v1 3 v1 v1 27
4
81
l
for z < 2
, and 1
81
∞
1 X 1 2n 4 z l/2 t
V (z, t) = ( ) δ(t − − (4n + 1) )
3 n=0 3 3 v2 v2
l
for 2 < z < l. The impulse response for the current is t
∞
1 X 1 2n z l 1 z l VHz,3 usL
-0.5
∞
1 X 1 4 z l/2
( )2n δ(t − − (4n + 1) )
-1.0
I(z, t) =
3Z2 n=0
3 3 v2 v2
l
for 2
< z < l. Using
l 2400
= = 16 µs
v1 150
and replacing δ(t) with f (t) = sin(2πt)u(t) the plot depicted in the margin was
obtained.
4
V+−V− V ++
Z1 Z2
+ +
Example 3: Two TL’s with characteristic impedances Z1 and Z2 are joined at a junc- V + + V −- -V ++
Z1 R Z2
tion that also includes a “shunt” resistance R as shown in the diagram in the
margin. Determine the reflection coefficient Γ12 and transmission coefficient τ12 z
at the junction.
5
is the parallel combination of R and Z2 . Solving these equations, we find that
V− Zeq − Z1
Γ12 ≡ =
V+ Zeq + Z1
and
V ++ 2Zeq
τ12 = + = .
V Zeq + Z1
z
Hint: in this ckt Γ12 has the
usual form in terms of Zeq ≡
R+Z2 . For τ12 , we need 1+Γ12
multiplied by a voltage divi-
Exercise: Two TL’s with characteristic impedances Z1 and Z2 are joined at a junction sion factor Z2 /(R + Z2 ).
that also includes a series resistance R as shown in the margin. Determine the
reflection coefficient Γ12 and transmission coefficient τ12 at the junction.
6
31 Periodic oscillations in lossless TL ckts +
C I(t) V (t)
• Lossless LC circuits (see margin) can support source-free and co-sinusoidal L -
voltage and current oscillations at a frequency of
1
ω=√ |V (z, t)|
LC |I|
1
• Consider first a TL segment of some length ℓ having no electrical con-
nection to any elements at either end, as shown in the margin.
|V (z, t)|
– Effectively, both ends of the TL have been “open circuited”, and |I|
thus TL current I(z, t) needs to vanish at z = 0 and z = ℓ. Since
f (t − vz ) g(t + vz )
I(z, t) = −
Zo Zo Open λ Open
4
2
Since no other constraint is imposed, any waveform with the spec-
ified period is admissible, and the most general such expression is
given by the Fourier series
∞
X
f (t) = Fo + Fn cos(nωot + θn )
n=1
3
Also1 ∞ 1.0
X
V (z, t) = 2Fn cos(nωot + θn ) cos(nβoz) 0.5
n=1
z HmL
jθn −jnβo z jnβo z 50 100 150 200 250 300
P
from the phasor Ṽ (z) = n Fn e [e +e ].
-0.5
In summary: -1.0
4
implying that
λ
ℓ=n ,
2
that is, the line length is an integer multiple of half-
wavelength at each resonance.
5
|I|
• For TL’s of length ℓ open at one end shorted at the other end, |V |
resonant wavelengths and frequencies can be identified by requiring ℓ
to be an odd multiple of λ4
– the reason for this is, the nulls of waveforms ∝ cos(βz) and sin(βz) Open λ
4
Short
6
1.0
0.5
Example 1: A lossless TL of 600 m length is left open at z = 0 and shorted at
z = l = 600 m. Determine (a) resonant frequencies of the line, (b) resonant 100 200 300 400 500 600
z HmL
current modes, (c) resonant voltage modes obtained from the current modes
using the telegrapher’s equations. The line has a characteristic impedance of -0.5
7
(c) Let’s find the voltage modes Vn (z, t) from the current modes above using one of
the telegrapher’s equations,
∂V ∂I
− =L .
∂z ∂t
Substituting In (z, t) into this equation and differentiating we find
∂V π π π
= −(2n + 1) L sin((2n + 1) z) cos((2n + 1) t).
∂z 4 1200 4
Next finding the anti-derivative of the above, we conclude
π π
Vn (z, t) = 300L cos((2n + 1) z) cos((2n + 1) t).
1200 4
A snapshot of the animation of resonant voltage modes is shown in the margin. 1.0
0.5
Note that the amplitudes of current and voltage modes cannot be assigned indepedent
of one another — above we set the current amplitudes to unity and obtained z HmL
100 200 300 400 500 600
300L for voltage amplitudes.
-0.5
-1.0
8
t = 1 µs
• How can one get source free oscillations in a TL?
t=0
One answer is, the TL might have been connected to a source in the 3V V (z, 1)
past before being disconnected from it. z
0 600 m
will be allowed on the source free line. (per mode) where K is the
Boltzmann constant. Lossy
Therefore, V (z, t) and I(z, t) for t > 1 µs can be expressed as a weighted
lines with finite conductivity
superposition of the resonant modes of the line with resonant frequen-
also produce thermal noise.
cies nωo, subject to the initial conditions V (z, 1) and I(z, 1) given
Thermal noise is easy to
above.
detect and routinely inter-
feres with weak communi-
cation signals that we care
about!
9
32 Input impedance and microwave resonators Series Parallel
• In the last lecture we learned that when a shorted stub is open circuited
at its input port, it shows resonance if the stub length ℓ is an odd
multiple of λ4 .
1
and the input port of the stub coincides with a voltage max and
a current null, i.e., I(0, t) = 0 — thus the input impedance Zin of
the stub is infinite at these resonances, just like the impedance of
the parallel LC-circuit depicted above.
– Thus this set of resonant frequencies are referred to as parallel
resonances of the shorted stub.
0
v f
Series Resonances
2ℓ
short
The diagram in the margin marks the locations of parallel and series
resonance frequencies of the shorted stub associated with infinite and zero l
input impedance Zin . +
Thus, a shorted stub, included in a circuit such as the one shown in the + RL VL(t)
margin, will exhibit extreme behavior at these special frequencies — namely fi(t) - -
it will appear as a
2
• short at its series resonances, causing the entire input signal fi (t) to
appear as VL (t) across the load RL , and v
Parallel Resonances
4ℓ
We next focus our attention on how Zin of the stub appears at other l
frequencies not coinciding with any of the resonances discussed
+
above.
+ RL VL(t)
- -
• In the following we will assume that the TL stub, as well as the circuit fi(t)
it is connected to, are all in sinusoidal steady state at a frequency
determined by the frequency of the sinusoidal source fi (t).
and
z z
Re{V +ejω(t− v )} − Re{V − ejω(t+ v )} ˜ = V + −jβz
e − V − jβz
e
I(z, t) = ⇔ I(z) ,
Zo Zo
where
ω
√
– β= v = ω LC is the wavenumber at frequency ω, and
3
– V + and V − are phasors of forward and backward propagating
voltage waves on the line evaluated at z = 0.
the TL.
growing to the left from the short termination toward the input terminal
of the line. d l 0
In that case, the input impedance of the shorted stub can be denoted
as
Ṽ (d = l)
Z(l) = ,
˜ = l)
I(d
where
+ jβd − −jβd ˜ V +ejβd − V − e−jβd
Ṽ (d) ≡ V e +V e and I(d) ≡ .
Zo
• An immediate benefit of our new notation comes when we apply the
voltage boundary condition at the short termination.
4
GENERATOR, LOAD
– We apply it as Circuit
Input
V (0, t) = 0 ⇔ Ṽ (0) = V + + V − = 0 port Zo Short
and thus d l 0
Ṽ (l)
Z(l) = = jZo tan(βl). Π Π 3Π 5Π 3Π 7Π
2Π
Βl
˜
Π
4 2 4 4 2 4
I(l)
-5
Note that: input impedance Z(l) = 0 + jX(l) is (see margin for X(l))
1. purely reactive for all l,
λ
• a shorted stub of length l = 4 appears at its input terminals like λ
l=
an open (see margin). 4
6
6. at l = λ2 the input impedance of the shorted stub returns back to
zero, which in turn indicates, in view of (5), that
short
λ
• an open ended stub of length l = 4
must appear at its input l
terminals like a short (see margin).
+
RL VL(t)
Next set of examples illustrate the uses of shorted/opened TL stubs as circuit +
- -
elements. f i (t)
Amplitude Resp.
1.0
Example 1: A shorted TL stub of length l = 3 m is connected in series with a resistor
RL = 50 Ω as shown in the diagram in the margin. Plot the magnitude of the 0.8
and v = c on the stub. Interpret the amplitude response curve |H(ω)| in terms 0.4
f HMHzL
ω
Solution: Using β = c and voltage division, we find that frequency response 50 100 150 200 250 300
ṼL RL 1
H(ω) = = = .
F̃i RL + jZo tan(βl) 1 + j tan( ωc l)
The plot of |H(ω)| with the given parameters is shown in the margin. The peaks
of the amplitude response occur at the series resonance frequencies of the shorted
stub when its input impedance is zero (an effective short). The nulls of the
amplitude response correspond to parallel resonances of the stub when it appears
like an open at its input terminals.
7
Series network:
short
Example 2: Consider a shorted TL connected at d = l to an inductor L. Determine
the resonances of the combined network. L l
Solution: The input impedance of the shorted line is
√
Z(l) = jZo tan(βl) = jZo tan(ω LCl)
whereas inductor L has an impedance ZL = jωL. If the inductor and shorted stub
are connected in series (see margin), then the series resonances of the network
will be observed when the network input impedance
√ Parallel network:
ZL + Z(l) = jωL + jZo tan(ω LCl)
short
equals zero. The parallel resonances of the network will be observed when the
impedance is infinite. While the series resonance frequencies of the stub will be
shifted because of the inductor, parallel resonances will not shift (infinities due l
to tan function cannot be shifted by the finite additive term due to the inductor).
The shifted series resonance frequencies ωn can be found graphically by plotting
ZL + Z(l) and looking for the zero crossings.
L
If the inductor and shorted stub are connected in parallel, then the parallel resonances
of the network will be observed when the network input admittance
1 1
YL + Y (l) = + √
jωL jZo tan(ω LCl)
equals zero (same as infinite input impedance). The series resonances, on the
other hand, will be observed when the admittance is infinite (same as zero input
impedance). Series resonances of the stub will not be shifted with, unlike its
parallel resonances. The shifted parallel resonance frequencies ωn will equal the
series resonance frequencies of the series connected network described above.
8
short
l
Example 3: A shorted TL stub of length l = 3 m is connected in series with a a +
C
capacitor C = 10 pf and a resistor RL = 50 Ω as shown in the diagram in the RL VL(t)
+
margin. Plot the magnitude of the frequency response H(ω) = ṼF̃L as a function - -
ω fi(t)
of frequency f = 2π if Zo = 50 Ω and v = c on the stub. Interpret the amplitude
Can change the short
response curve |H(ω)| in terms of resonance frequencies of the shorted line. location on the stub
to tune the filter
ω
Solution: Using β = c and voltage division, we find that frequency response
ṼL RL 1
H(ω) = = = . Amplitude Resp.
F̃i RL + 1
jωC + jZo tan(βl) 1 + 1
jωRL C + j tan( ωc l) 1.0
0.8
The plot of |H(ω)| with the given parameters is shown in the margin. The peaks 0.6
of the amplitude response occur at the shifted series resonance frequencies of the
0.4
shorted T.L. stub. The nulls of the amplitude response correspond to parallel
resonances of the stub when it appears open. 0.2
f HMHzL
50 100 150 200 250 300
9
Example 4: (a) If in the TL circuit shown in the margin IR = 2∠0o A, what is the
line length l in terms of wavelength λ of the given source frequency on the line?
Solution: (a) If IR = 2∠0o A, then KCL application at the source terminal implies
that I(l) = 0.
(b) If IR = 0, then V (l) = (50 Ω)IR = 0, implying that the T.L. has a short at d = l. I(l) I(0) = 0
IR + +
26 V (l)
Since d = 0 is an open, we need to have l some odd multiple of λ4 .
0A V (0)
50 Ω - -
d
l 0
10
33 TL circuits with half- and quarter-wave trans-
formers GENERATOR, LOAD
Circuit I(d)
1
• In the TL circuit shown in the margin an arbitrary load ZL is connected
to a TL of length l = λ2 at the source frequency.
Given that
λ 2π λ
e±jβ 2 = e±j λ 2 = e±jπ = −1, Half-wave transformer:
– inverts the algebraic sign of its voltage and current inputs at the
load end (and vice versa), and
– has an input impedance identical with the load impedance since
Vin −VL
Zin ≡ = = ZL.
Iin −IL
These very simple results are easy to remember and use.
2
• In the TL circuit shown in the margin an arbitrary load ZL is connected
to a TL of length l = λ4 at the source frequency. Quarter-wave transformer:
+
+ 50 + j50 Ω
Vg = 1006 0 V -
Example 2: The load and the transformer of Example 1 are connected to a source Zo = 50 Ω -
with voltage phasor Vg = 100∠0o V at the input port. What is the load current λ 0
IL and what is the average power absorbed by the load? 4
Solution: Since Vin = Vg = 100∠0o V, the current-forcing formula for the quarter-wave
transformer implies
Vin 100
IL = −j = −j = −j2 A.
Zo 50
To find the average power absorbed, we first note that load voltage
VL = ZLIL = (50 + j50)(−j2) = 100 − j100 V.
Thus,
1 1
PL = Re{VL IL∗ } = Re{(100 − j100)(j2)} = 100 W.
2 2
4
I(d)
Zg = 25 Ω
Vg = +
ZL =
V (d)
j10 V +- Zo = 50 Ω
100 Ω
-
Solution: First we determine input impedance Zin by noting that ZL = 100 Ω trans-
forms to itself, namely 100 Ω at d = 0.5λ, but then it transforms from d = 0.5λ
to 0.75λ as
Zo2 502
Zin = = = 25 Ω.
Z(0.5λ) 100
Hence, using voltage division, we find,
Zin 25
Vin = Vg = j10 = j5 V.
Zg + Zin 25 + 25
Next, using half-wave transformer rule, we notice that
and finally applying the quarter-wave current forcing equation with V (0.25λ) we
get
V (0.25λ) −j5
IL = −j = −j = −0.1 A.
Zo 50
Clearly, then, the load voltage is
5
Vin
IL2 = −j
Zo
Example 4: In the circuit shown in the margin, ZL1 = 50 Ω, ZL2 = 100 Ω, and Zo1 = ZL2
Zo2 Vin
Zo2 = 50 Ω. Determine IL1 and IL2 if Vin = 5 V. Both T.L. sections are quarter- IL1 = −j
Zo
wave transformers.
+
Vin Zo1 ZL1
Solution: Using the current-forcing equation, we have -
Vin 5 λ 0
IL1 = IL2 = −j = −j = −j0.1 A. 4
Zo 50
Consequently,
VL1 = IL1ZL1 = −j0.1 A × 50 Ω = −j5 V
and
VL2 = IL2ZL2 = −j0.1 A × 100 Ω = −j10 V.
Thus, total avg power absorbed is
1 ∗ 1
P = Re{VL1IL1 } + Re{VL2I ∗L2}
2 2
1 1
= = Re{−j5 × j0.1} + Re{−j10 × j0.1} = 0.75 W.
2 2
6
34 Line impedance, generalized reflection coef- LOAD
GENERATOR,
I(d)
ficient, Smith Chart Circuit
+
Input
port Zo V (d)
• Consider a TL of an arbitrary length l terminated by an arbitrary load - ZL = RL + jXL
ZL = RL + jXL. −l l 0 z
• Our analysis starts at the load end of the TL where V (0) and I(0)
stand for the load voltage and current, obeying Ohm’s law
V (0) = ZL I(0).
+ − V+−V− ZL − Zo +
V +V = ZL ⇒ V− = V .
Zo ZL + Zo
1
– Define a load reflection coefficient “Load reflection coefficient”
ZL − Zo is a well justified name for
ΓL ≡
ZL + Zo ΓL since the forward travel-
ing wave with phasor V + ejβd
and re-write the voltage and current phasors as
gets reflected from the load.
+ jβd −j2βd V +ejβd [1 − ΓLe−j2βd]
V (d) = V e [1 + ΓLe ] and I(d) = .
Zo
– Define a generalized reflection coefficient The term “generalized reflec-
tion coefficient” is also well
Γ(d) ≡ ΓLe−j2βd
justified even if there is no
and re-write the voltage and current phasors as reflection taking place at ar-
bitrary d — the reason is, if
+ jβd V +ejβd[1 − Γ(d)]
V (d) = V e [1 + Γ(d)] and I(d) = . the line were cut at location
Zo
d and the stub with the load
– Line impedance is then defined as were replaced by a lumped
V (d) 1 + Γ(d) load having a reflection co-
Z(d) = = Zo
I(d) 1 − Γ(d) efficient equal to Γ(d), then
there would be no modifica-
for all values of d on the line extending from the load point d = 0
tion of the voltage and cur-
all the way to the input port at d = l.
rent variations on the line to-
wards the generator.
With the dependence on d of Z(d) as well as Γ(d) tacitly implied,
we can re-write this important relation and its inverse as Each location d on the line
2. Since
|Γ| = |ΓL| and ∠Γ(d) = ∠ΓL − 2βd
property (1) implies that Γ(d) is a complex number which is constrained
to be on or within the unit-circle on the complex plane.
3. Relationships
Z 1+Γ Z − Zo
= ⇔ Γ=
Zo 1 − Γ Z + Zo
between Γ and Z are known as bilinear transformations — here the
term bilinear refers to the numerator as well as the denominator of
these transformations being linear in the variable being transformed
(from right to left).
3
Bilinear (or Möbius) transformations are known to have the general
property of mapping straight lines into circles on the complex num-
ber plane.
and
Z Unit
≡ z ≡ r + jx, Radius
Circle
Zo
known as normalized impedance, lead to an ingenious graphical aid Re
4
• SC can be constructed by first noting that
z − 1 r + jx − 1 [(r − 1) + jx][(r + 1) − jx] (r2 + x2 − 1) + j2x
Γ= = = = ≡ Γr +jΓi ;
z + 1 r + jx + 1 (r + 1)2 + x2 (r + 1)2 + x2
thus Γi
2 2
(r + x − 1) 2x SmithChart
Γr = and Γ i = , 1
(r + 1)2 + x2 (r + 1)2 + x2 .5 2
r 2 1 2 1 1 2. + 1. ä
r+1 r+1 x x 0
-1
.2 .5
0.4 - 0.2 ä
1 2 2.61803 r=5
1
Γr
describing r and x dependent circles, respectively, on complex plane -.2 x=-5
• Typical SC usage: -1
5
SmithChart
– Rotation by an angle of 2βd amounts to 1
3. Also, 0 .2 .5 1
80.447214, 0.0368959<
2 2.61803 r=5
1 0.4 - 0.2 ä
y(d) ≡
z(d) -.2 x=-5
z(d). -1
Verification: Since
1+Γ 1 1 − Γ 1 + (−Γ)
z= ⇒ y= = = ;
1−Γ z 1 + Γ 1 − (−Γ)
hence whereas z is the transform of Γ, y is the transform of −Γ,
having the same magnitude as Γ but an angle off by ±180◦.
– Therefore, “reflect” on the SC across the origin to jump
from z(d) to y(d) if you need the value of the normalized
admittance.
6
I(d)
+
Zo = 50 Ω V (d) ZL = 50 + j100 Ω
-
Example 1: A transmission line is terminated by an inductive load of
d 0.125λ 0
ZL = 50 + j100 Ω.
ZL 50 + j100
z(0) = = = 1 + j2.
Zo 50 -.2 0.2 - 0.4 ä x=-5
λ
Enter z(0) on the SC and then rotate clockwise by 8 ⇔(quarter circle) to obtain
the normalized input impedance
-.5 -2
-1
.2 x=5
-.2 1. - 2. ä x=-5
80.707107, -0.0625<
-.5 -2
7 -1
Blow up of the SC’s used in Example 1:
SmithChart SmithChart
1 1
.5 2 .5 2
1. + 2. ä 0.2 + 0.4 ä
0 .2 .5 1 2 r=5
5.82843 0 .2 .5 1 2 r=5
5.82843
80.707107, -0.0625<
-.5 -2 -.5 -2
8
35 Smith Chart examples .5
SmithChart
1
.2 x!5
2. # 1. $
!0.447214, 0.0368959"
Example 1: A load ZL = 100 + j50 Ω is connected across a TL with Zo = 50 Ω and 0 .2 .5 1 2 2.61803 r!5
l = 0.4λ. At the generator end, d = l, the line is shunted by an impedance 0.4 " 0.2 $
Zs = 100 Ω. What are the input impedance Zin and admittance Yin of the line, ".2 x!"5
z(0) = = = 2 + j1
Zo 50 .5 2
is entered in the SC shown in the margin on the top. Clockwise rotation (from 0.600056 # 0.663401 $
!0.447214, 0.136896"
load toward generator) at fixed |Γ| (red circle) by
.2 x!5
takes us to ".2
0.749913 " 0.829077 $
x!"5
for the overall normalized input admittance of the shunted line as shown on the
SC in the bottom — the corresponding normalized input impedance is .2 0.555603 # 0.368536 $ x!5
!0.361132, 0.176393"
1
zin = ≈ 0.56 + j0.37. 0 .2 .5 1 2
2.13054 r!5
Hence, the unnormalized input impedance and admittance are ".2 x!"5
"1
1
I(d)
Zg = 25 Ω
Vg = +
Zs = 100 Ω ZL =
100! 0 V +- V (d)
100 + j50 Ω
Example 2: The TL network described in Example 1 is connected to a generator with
Zo = 50 Ω
-
Solution:
(b) The shunt element Zs = 100 Ω sees the same voltage Vin and conducts a current
Vin /Zs . Therefore it absorbs an average power of
1 Vin ∗ |Vin |2 |Vg Zin|2
P = Re{Vin ( ) } = =
2 Zs 2Zs 2Zs|Zg + Zin |2
|100 · (27.8 + j18.4)|2
≈ ≈ 17.78 W.
2 · 100 · |25 + 27.8 + j18.4|2
The remainder of 44.44 W will be absorbed in ZL .
2
I(d) SmithChart
1
+ .5 2
Zin = 50 + j50 Ω V (d)
Zo = 50 Ω
ZL =?
- 1. # 1. $
.2 x!5
!0.447214, 0.0881041"
"1
Zin 50 + j50 .5 2
zin = = =1+j
Zo 50 0.594079 # 0.655216 $
in the SC as shown in the margin on the top. Counter-clockwise rotation (from .2 x!5
!0.447214, "0.111896"
z(0) ≈ 0.76 − j0.84 and y(0) ≈ 0.59 + j0.66
as shown on the next SC at the load point. Hence, we find ".5 "2
and
1
YL = Yo y(0) ≈ (0.59 + j0.66) = 0.012 + j0.013 S.
50
3
SmithChart
ΓLs = −0.5 1
.5 2
Zo = 50 Ω
.2 x!5
I(d)
0.333333 3.
+ .2 .5 1 2 3. r!5
!0.5, 0"
0
V (d)
Zin =? Zo = 50 Ω
ΓL = 0.5
- ".2 x!"5
"1
SmithChart
1
.5 2
Solution: Recall that the SC covers the unit circle of the complex plane and therefore ".2 !0.5, "0.2" x!"5
can be entered directly in the SC as shown on the top SC in the margin. Clockwise "1
SmithChart
rotation (from load toward generator) at fixed |Γ| (red circle) by 1
.5 2
◦ ◦
0.2λ ⇔ 0.4 × 360 = 144
takes us to
.2 x!5
0.333333 3.
!0.5, 0.25"
.5 1 2 3. r!5
"1
4
SmithChart
1
.5 2
for the shunt connected stub in the third SC and rotating clockwise by
.2 .5 1 2 3. r!5
0.3λ ⇔ 0.6 × 360◦ = 216◦
0
!0.5, "0.05"
we obtain ".2 x!"5
"1
.5 2
and entering it in the next SC. Finally rotating clockwise once again by .2
2.065 # 1.61837 $
x!5
"1
SmithChart
1
.5 2
.2 x!5
0.283936 # 0.0577194 $
!0.558955, "0.00997561"
".2 x!"5
".5 "2
"1
5
Example 5: What is the load impedance ZLs terminating the shunt connected stub
in Example 4?
ΓLs = −0.5,
it follows from the bilinear transformation linking zLs and ΓLs that
1 + ΓLs 1 − 0.5 1
zLs = = = .
1 − ΓLs 1 + 0.5 3
Hence, the impedance is
50
ZLs = Zo zLs = Ω.
3
Solution: This is similar to Example 5. Given that the load reflection coefficient is
ΓL = 0.5,
it follows from the bilinear transformation linking zL and ΓL that
1 + ΓL 1 + 0.5
zL = = = 3.
1 − ΓL 1 − 0.5
Hence, the impedance is
ZL = Zo zL = 150 Ω.
6
36 Smith Chart and VSWR Generator
I(d)
Transmission line
Wire 2
+ Load
Zg
• Consider the general phasor expressions
+ V (d) Zo
- ZL
F = Vg - Wire 1
0
V +ejβd(1 − ΓLe−j2βd )
l d
+ jβd
V (d) = V e (1 + ΓLe −j2βd
) and I(d) = |V (d)|max
|V (d)|
Zo
describing the voltage and current variations on TL’s in sinusoidal |V (d)|min dmax
steady-state. dmin
– Unless ΓL = 0, these phasors contain reflected components, which Complex addition displayed
graphically superposed on a
means that voltage and current variations on the line “contain” Smith Chart
.5 2
"1
at locations d = dmax and dmin such that |1 + Γ(d)| maximizes for d = dmax
Γ(dmax) = ΓLe−j2βdmax = |ΓL| and Γ(dmin) = ΓLe−j2βdmin = −|ΓL|,
|1 + Γ(d)| minimizes for d = dmin
and such that Γ(dmin) = −Γ(dmax)
λ
dmax − dmin is an odd multiple of .
4
1
– These results can be most easily understood and verified graphi- Generator Transmission line
I(d) Wire 2
- ZL
F = Vg - Wire 1
0
• We define a parameter known as voltage standing wave ratio, or
l d
|V (d)|max
VSWR − 1
|V (d)|min
|V (dmax )| 1 + |ΓL| dmax
VSWR ≡ = ⇔ |ΓL| = .
VSWR + 1
dmin
|V (dmin)| 1 − |ΓL|
Notice that the VSWR and |ΓL| form a bilinear transform pair just Complex addition displayed
graphically superposed on a
like Smith Chart
1+Γ z−1
z= ⇔ Γ= . SmithChart
1
1−Γ z+1 .5 2
Since
1 + Γ(dmax)
Γ(dmax ) = |ΓL| ⇒ VSWR = , .2
) Γ(d) z(dmax )
x!5
1 − Γ(dmax ) 1 + Γ (d =VSWR
this analogy between the transform pairs also implies that 0 .2 .5 1 2 VSWR r!5
1 Γ(dmax ) = |ΓL|
z(dmax ) = VSWR, ".2 x!"5
as explicitly marked on the the SC shown in the margin . Consequently, ".5 "2
"1
– the VSWR of any TL can be directly read off from its SC plot |1 + Γ(d)| maximizes for d = dmax
2
• The extreme values the VSWR can take are:
SmithChart
1
.5 2
.2
z(d x!5)
|1 + Γ(d)| maximizes for d = dmax
Γ ( d) Γ(d) max
1+ =VSWR
0 .2 .5 1 2 VSWR r!5
1 Γ(dmax ) = |ΓL|
x!"5
".2
|1 + Γ(d)| minimizes for d = dmin
such that Γ(dmin) = −Γ(dmax)
".5 "2
"1
3
• In the lab it is easy and useful to determine the VSWR and dmax or Generator Transmission line
I(d) Wire 2
- ZL
F = Vg - Wire 1
VSWR − 1 |V (d)|max
|V (d)|
|ΓL| =
VSWR + 1
|V (d)|min
is easily determined, and
dmax
dmin
SmithChart
1
1 Γ(dmax ) = |ΓL|
it follows that
".2 x!"5
1 + ΓL
ΓL = |ΓL|ej2βdmax ⇒ zL = .
1 − ΓL ".5 "2
"1
• alternatively, zL can be obtained directly on the SC by rotating counter- |1 + Γ(d)| maximizes for d = dmax
clockwise by dmax from the location of
|1 + Γ(d)| minimizes for d = dmin
z(dmax ) = VSWR. such that Γ(dmin) = −Γ(dmax)
4
SmithChart
1
.5 2
by the load.
".2 x!"5
Solution:
(a) As shown in the top SC in the margin, VSWR=4 is entered in the SC as ".5 "2
z(dmax ) = 4 + j0, and constant |ΓL | circle is then drawn (red circle) passing "1
SmithChart
through z(dmax ) = 4. 1
.5 2
A counter-clockwise rotation from z(dmin) = 0.25 by one fourth of a full circle corre- .2 x!5
zL ≈ 0.4706 − j0.8823
".2 x!"5
as shown in the second SC. Hence, this gives 0.470588 " 0.882353 #
!0.6, "0.125"
"1
5
Therefore the voltage and current phasors at the voltage minimum location are
20 V
V (dmin) = 20 V and I(dmin) = .
12.5 Ω
Generator Transmission line
Average power transported toward the load at d − dmin is, therefore, I(d) Wire 2
+ Load
1 1 20 400 Zg
P (dmin ) = Re{V (dmin)I(dmin)∗} = Re{20 }= W = 16 W. +
- V (d) Zo ZL
2 2 12.5 25 F = Vg - Wire 1
l d 0
Since the TL is assumed to be lossless we should have |V (d)|max
|V (d)|
PL = P (dmin ) = 16 W.
|V (d)|min dmax
dmin
Solution: Given that l = 0.625λ and dmin = 0.125λ, we note that there is just one
half-wave transformer between l = 0.625λ and dmin = 0.125λ. Therefore
Vin = −V (dmin) = −20 V and Zin = Z(dmin) = 12.5 Ω.
But also
Zin
Vin = Vg .
Zg + Zin
Consequently,
Zg + Zin 50 + 12.5 62.5
Vg = Vin = −20 = −20 = −100 V.
Zin 12.5 12.5
6
SmithChart
1
.5 2
Example 3: Determine V + and V − in the circuit of Examples 1 and 2 above such that .2 x!5
the voltage phasor on the line is given by
0.25 4.
V (d) = V + ejβd + V − e−jβd . 0 .2 .5 1 2 r!5
4.
!0.6, 0.25"
Solution: Looking back to Example 1 (also see the SC’s in the margin), we first note ".2 x!"5
that
VSWR − 1 4 − 1
|ΓL | = = = 0.6 = Γ(dmax) = −Γ(dmin).
VSWR + 1 4+1 ".5 "2
+ −j π4
V = 50e V.
.2 .5 1 2 r!5
Since 0 4.
j π2
ΓL = Γ(0) = Γ(dmin)ej2βdmin = −0.6e ,
it follows that ".2
0.470588 " 0.882353 #
x!"5
π π π
V − = ΓL V + = −0.6ej 2 × 50e−j 4 = −30ej 4 V. !0.6, "0.125"
".5 "2
"1
7
SmithChart
1
.5 2
Example 4: Determine the load voltage and current VL = V (0) and IL = I(0) in the .2 x!5
circuit of Examples 1-3 above.
0.25 4.
Solution: In general, 0 .2 .5 1 2 r!5
4.
!0.6, 0.25"
V + ejβd − V − e−jβd
V (d) = V + ejβd − V − e−jβd and I(d) = . ".2 x!"5
Zo
Therefore,
".5 "2
+ −
V −V "1
VL = V (0) = V + + V − and IL = I(0) = . SmithChart
Zo 1
.5 2
Using Zo = 50 Ω and 0.470588 $ 0.882353 #
−j π4 j π4
V + = 50e V and V − = −30e V .2 x!5
!0.6, "0.125"
".5 "2
"1
8
37 Smith Chart and impedance matching
• In lossless TL circuits the average power input Pin at the generator end I(d)
+ Load
precisely matches the average power delivered to the load, PL. +
-
Zg
V (d) Zo ZL
F = Vg -
In fact, Pin and PL also match the average power P (d) transported on
the line at an arbitrary d. d l 0
• We have in general
1
P (d) = Re{V (d)I ∗(d)}
2
+ jβd
1 + jβd − −jβd V e − V − e−jβd ∗ Power tx’ed toward the load:
= Re{(V e + V e )( )}
2 Zo |V +|2
1 |V + |2 |V −|2 V − V +∗e−j2βd − (V −V +∗e−j2βd )∗ .
= Re{ − + } 2Zo
2 Zo Zo Zo
|V +|2 |V − |2
= − . Power tx’ed toward the gen-
2Zo 2Zo
erator:
– Note that P (d) is the difference of power transported |V − |2
|V + |2 .
2Zo
toward the load by the “forward-going” wave, and 2Zo
|V − |2
2Zo
toward the generator by the reflected wave.
– Also note that Power reflection coefficient:
Impedance matching achieves VSWR=1 between the generator and the match-
ing network inserted at a location between the load and the generator.
• The inserted network should be designed to yield an input impedance
equal Zo at its input terminals.
The following examples illustrate different ways of achieving an
impedance match.
2
Quarter-wave matching
λ
4
q
Zq = Zo Z(d1 )
Zo ZL
Example 1: Quarter-wave matching of resistive loads:
To reduce the VSWR on the line connected to the generator to unity, we can insert
a quarter-wave transformer right after ZL — i.e., at d1 = 0 in the circuit
shown in the margin — with a characteristic impedance
√ √
Zq = 25 × 50 = 1250 = 35.35 Ω.
The impedance at the input terminals of the quarter-wave transformer (on the left)
is then Zo , i.e., 50 Ω, implying a perfect impedance match.
3
Quarter-wave matching
λ
4
q
Zq = Zo Z(d1 )
Example 2: Quarter-wave matching of reactive loads: Zo ZL
yields a real valued Zq as needed. This procedure leads to having d1 = dmax or 0.5 - 0.5 ä
x=-5
d1 = dmin corresponding to the positions of voltage maxima and minima on the -.2
line.
-.5 -2
As shown in the margin,
-1
for
d1 ≈ 0.250λ − 0.162λ = 0.088λ z(d1 ) = z(dmax ) = VSWR ≈ 2.62
is a suitable choice for quarter-wave matching. In that case we need as marked on the SC.
√ √ Also
Zq = 131 × 50 = 50 × 2.62 Ω
dmax ≈ 0.088λ
for the quarter wave transformer in order match to load to a line with Zo = 50 Ω. since, as marked on the SC,
the angle of ΓL is 0.088λ.
4
Single-stub
tuning Shorted
stub
ls
y(d1) = 1 + jb
so that a stub, with a normalized input admittance .2 x=5
0.4 + 0.2 ä
ystub = −jb,
.2 .5 1 2 2.61803 r=5
can be added in parallel to have a combined admittance of 0
2. - 1. ä
In specific -.5 -2
ZL 1
zL = = 2 − j1 and yL = = 0.4 + j0.2
Zo zL
-1
SmithChart
1
as shown on the SC on the top in the margin. We rotate clockwise on the SC by
an amount corresponding to d1 to obtain .5 2
y(d1 ) = 1 + j1 1. + 1. ä
.2 x=5
on the “g = 1” or “y = 1 + jb” circle as shown in the bottom SC. From the
amount of rotation we determine
0 .2 .5 1 2 2.61803 r=5
d1 ≈ 0.162λ − 0.037λ = 0.125λ.
0.5 - 0.5 ä
x=-5
5
-.2
80.447214, -0.161896<
The required input impedance of the shorted stub to achieve
y(d1) + ystub = 1 + j0
is
ystub = −1j. SmithChart
1
.5 2
To achieve this input admittance the required stub length is
λ
ls = = 0.125λ .2 x=5
8 0.4 + 0.2 ä
-.5 -2
.5 2
3 λ8 .
.2 x=5
– d2 = d1 +
Vary l1 and l2 until VSWR is reduced to 1 near the generator end. 0 .2 .5 1 2 2.61803 r=5
The advantage of double-stub tuning is avoiding changes of stub loca- 0.5 - 0.5 ä
-.2 x=-5
tions when ZL is changed. It’s implementation on a SC is considerably 80.447214, -0.161896<
-1
6
38 Distribution networks
• A corporate ladder network that combines 4
CORPORATE LADDER NETWORK
identical loads ZL into a single equivalent input ZL ZL ZL ZL
Loads
impedance ZL is shown in the margin.
1
• A hybrid combiner network shown in the mar- HYBRID COMBINER --- SUM and DIFFERENCE OF ISOLATED
gin can be used to excite two identical TL loads GENERATOR SIGNALS V_A AND V_B ARE APPLIED TO LOADS R
√ R √
nal generators VA and VB having equal internal R
Zo = 2R
+
Zo = 2R
R
LOAD - VA LOAD
resistances R matched to the load resistance. INPUT A Zo =
√
2R
√
Zo = 2R
Zo = 2R, INPUT B
• This very special situation allows one to calculate the various terminal voltages on the hybrid
due to VA and VB one-at-a-time as if loads R were isolated from generator-B and -A (by “virtual
shorts” existing across generator terminals when VB and VA are suppressed) in turns, and then
superpose the results.
2
• Terminal voltages obtained with that procedure (those shown on the diagram) turn
out to be valid when both generators are active as can easily be checked for self-
consistency by using the current-forcing equations introduced earlier. For instance,
the total current into generator-A terminal (flowing from both sides) is
j VA − VB j VA + VB VA
IA = − √ (−j √ ) − √ (−j √ ) = − ,
R 2 2 2 R 2 2 2 2R
and hence the voltage drop from the same terminal to the ground is
VA VA
IA R + V A = − R + VA =
2R 2
as marked explicitly on the diagram. All self-consistency tests that can be applied
with the given expressions are passed, and so the results given are valid.
• The input and output ports of a HYBRID COMBINER --- SUM and DIFFERENCE OF ISOLATED
hybrid combiner can be swapped GENERATOR SIGNALS V_A AND V_B ARE APPLIED TO LOADS R
√ R √
impedance R, and output signals the R
Zo = 2R
+
Zo = 2R
R
LOAD - VA LOAD
sum and difference of generator volt- INPUT A Zo =
√
2R
√
ages. Zo = 2R
VB
2
√ √ R
Zo = 2R Zo = 2R
λ +
- VB
4
INPUT B
3
39 Lossy lines IDEAL LOSSLESS T.L.:
∂V ∂I
− = (jωL + R)I and − = (jωC + G)V.
∂z ∂z
1
• We will next show that
2
While waves governed by lossy γ and Zo (see margin) can exhibit substan-
tially different beavior than the lossless waves (examined in the previous
sections), at high frequencies the wave properties are reasonably similar as
alluded in item (2) above. We next examine this simplified high-frequency
limit.
p
• At high frequencies ω, such that ωL ≫ R and ωC ≫ G, we have γ= (jωL + R)(jωC + G)
– characteristic impedance s
jωL + R
Zo =
s r
jωL + R L jωC + G
Zo = ≈
jωC + G C
with
√ R G 1 R
β ≈ ω LC and α ≈ β( + ) = ( + GZo).
q
2ωL 2ωC 2 Zo
L L
1
In fact Zo = C
is exact even for a lossy line if R
= GC .
3
2π
• Note that β = λ
is the same as in the lossless case, and since
R G
α ≈ β( + ) ≪ β,
2ωL 2ωC
the “penetration depth” δ ≡ α1 of voltage and current waves on the TL
is much longer than a wavelength λ = 2π β
in this regime.
In summary, in the high-frequency regime, characteristic impedance Zo
and wavenumber β are (practically) the same as they are on lossless lines,
but signals do attenuate by a factor e±αd which should not be (and cannot
be) neglected over long distances d exceeding many wavelengths λ.
• At lower frequencies where the above approximations cannot be justi-
fied, a more careful analysis of lossy line equations is warranted.
• Finally, for an air-filled coax with inner and outer radii a and b, it can
be shown that the attenuation constant
q
f µo 1 b
1R 1 πσ b (1 + a )
α= = ηo ,
2 Zo 2 2π ln( ab )
which minimizes, at a fixed outer radius b, for ab ≈ 3.6, which in turn
results in an “optimal” characteristic impedance of
ηo b b
Zo = ln( ) = 60 ln( ) Ω ≈ 75 Ω
2π a a
for the same coax. Note that this result is independent of σ, the con-
ductivity of inner and outer conductors of the coax.
4
– For a dielectric filled coax having ǫ = 94 ǫo — implying vp = 23 c =
2×108 m/s — the same ratio ab ≈ 3.6 of outer and inner conductor
radii leads to Zo ≈ 50 Ω, the most common Zo encountered in
practical applications.
– The above result should also explain why having a thicker coax —
larger b — is better when losses are a concern.