EM Wave Theory and Vector Calculus
EM Wave Theory and Vector Calculus
Prof. J. K. Goswamy
UIET, Panjab University
Chandigarh.
Overview of Presentation
❖ Vector Calculus
❖ Maxwell’s Equations
❖ EM Waves Propagation
❖ Laws of Optics
Vector Calculus
Introduction
The development of the physics has been aided
by mathematical tools. These tools generalize
physical relations thereby offering deeper
interpretation.
The understanding and applications of calculus
and vectors have played an important role in
treatise of physics.
These two, when merged together, yield vector
calculus, which is one of the potent tools.
We plan to introduce concepts of vector calculus.
Concept of Field
The field can be defined as the region of space in which a physical
quantity has unique value at each point.
U U U
dU = dx + dy + dz
x y z
E E E
dE = dx + dy + dz
x y z
Gradient
Consider a function f(x,y,z) dependent on three variables. Let’s give an
increment of dx, dy and dz to these variables, then change in f(x,y,z) is
given to be:
df =
f
x
dx +
f
y
dy +
f
z
f
dz = iˆ +
f ˆ f ˆ ˆ
j+ (
k i dx + ˆjdy + kˆdz )
x y z
( )( )
df = f dl
f ˆ f ˆ f ˆ
where f = i+ j+ k = Gradient of f ( x, y, z )
x y z
( )( )
df = f dl = f dl cos
(df )max = f dl ( when cos = 1)
It does not depend anyway on the path taken between the end
points. This is referred to as fundamental theorem for gradients.
( )
b
f dl = f (b) − f (a)
a
▪ abstract in nature
The 𝛻 operator acts on the scalar or vector field to generate the following:
• Curl gives vector resultant through cross product 𝛻 × 𝐴Ԧ with vector field.
Divergence of a Vector Field
The scalar product of the vector field function 𝐸 𝑥, 𝑦, 𝑧 , defined at point P,
= E dS = ( E cos )dS = En dS
S S S
The scalar product gives the product of normal component of vector field
function with the area of the surface. It is also referred to as normal flux.
The term flux was first introduced to describe velocity of flow of liquids as it
would yield volume of flow per unit time. However later, it was also adopted
for the quantities which actually do not flow.
The divergence of a vector field is the ratio of
B 0 Field is divergent from source. B 0 Field is convergent from source. B = 0 Field is uniform.
Gauss’ Divergence Theorem
E dS = DivE dV
V
orientation of the area bound by the closed curve with respect to the direction
of the vector field. It gets maximized for a particular orientation of the vector
area enclosed by the closed curve and is called maximum circulation.
into large number of smaller loops, then the net circulation is the sum of
circulation over the constituting loops. Hence we can write mathematically as:
n
B dl = B dl
C i =1 Ci
Curl of Vector Field: The maximum value of circulation per unit area drawn
around a point, when the area enclosed by the curve shrinks to vanishing
value, gives the curl of the vector field. It can be expressed mathematically as:
1
CurlB = Lim
S →0 S
C
B dl
Interpretation: The curl is associated with the whirling motion of a vector field
about the point. It is a pseudo vector along the direction of vector area.
Components: The three components of Curl of vector field are given as:
(curlB ) ( )
1
x = curlB iˆ = Lim
S x →0 S
x Sx
B d l
(curlB ) ( )
1
y = curlB ˆj = Lim
S y →0 S
y Sy
B d l
(curlB ) ( )
1
z = curlB kˆ = Lim B dl
S z →0 S
z Sz
Curl of a Vector
The curl of a vector field is given as :
iˆ ˆj kˆ
ˆ Bz B y B
ˆj Bx − Bz + kˆ y − Bx
B = = i − +
x y z y z z x x y
Bx B y Bz
B dl =
C
curlB dS
S
Stoke’s Theorem
Irrotational Fields Solenoidal Fields
Vector fields, that obey the Vector fields, that obey the
relation 𝛻 × 𝐸=0, are irrotational relation 𝛻. 𝐸 =0, are solenoidal
fields. fields.
d
( f + g) = +
df dg
dx dx dx
In the similar spirit we can have the following relations related to gradient,
divergence and curl of sum of fields:
( f + g ) = f + g
( A + B) = A + B
( A + B) = A + B
Rule II(Multiplication by a constant)
d
(kf ) = k df
dx dx
( kf ) = kf
( kA) = k A
( kA) = k A
Rule III (Product Rule)
d
( fg ) = f + g
dg df
dx dx dx
( fg) = fg + gf
( ) ( ) (
)( ) ( )
A B = A B + B A + A B + B A
( )
fA = f A + A f
( )
(
( A B) = B A − A B)
( )
fA = f A + A f
( ) ( ) ( ) ( )
( A B) = B A − A B − A B − B A
Rule IV(Quotient Rule)
df dg
g − f
d f dx dx
=
dx g g2
f gf − fg
=
g g2
A g A − A g
=
g g2
A g A − A g
= 2
g g
Second Derivatives
Divergence of Gradient of scalar field
( )
2 2 2
= + 2 + 2 = 2
x 2
y z
We define Laplacian operator as :
2 2 2
= 2+ 2+ 2
2
x y z
Calculate the surface integral of the vector field 𝑣Ԧ = 2𝑥𝑧𝑖Ƹ + 𝑥 + 2 𝑗Ƹ + 𝑦 𝑧 2 − 3 𝑘 over five
surfaces of a cube of side=2 and having one of it vertices at origin of coordinate system.
Exclude the surface of cube lying in the xy plane of the coordinate system.
Check the divergence theorem using the function 𝑣Ԧ = 𝑦 2𝑖Ƹ + 2𝑥𝑦 + 𝑧 2 𝑗Ƹ + 2𝑦𝑧𝑘 and the unit
cube having of its vertices situated at the origin.
Check the Stoke’s theorem for the vector field 𝑣Ԧ = 2𝑥𝑧 + 𝑦 2 𝑖Ƹ + 4𝑦𝑧 2𝑘 and the unit side
square lying in yz plane with one of its vertices situated at origin.
Integration of Fields
B dr = Bx dx + B y dy + Bz dz
• The closed line integral represents the value of function over a
closed loop and are expressed as:
B dl
✓ Usually the line integral depends upon the path
taken from a to b.
✓ Vector field whose line integral is independent of
the path is known as conservative field.
Surface Integrals
• The open surface integral is cumulative value of field over an open
surface and is represented as
B dS = Bx dS x + B y dS y + Bz dS z
• The surface area elements dSx is parallel to YZ plane. Similar is the
case for other surface area elements.
• The closed surface integral is cumulative value of field over a
closed surface and is represented as
B dS
✓ For open surface, the direction of the surface
integral can be inwards or outwards.
✓ For closed surface, the direction of surface
integral is always considered outwards.
Volume Integral
1 n
qi conjunction with the
E=
4 0 i =1 ri
rˆ
2 i principle of superposition,
constitutes the basic
If a charge Q suffers force (F) due to distribution of
framework to deal with the
discrete charges
problems in electrostatics.
F = QE
Electric Field due to Continuous Charge Distribution
If we have a continuous distribution of charges over a line or surface or
volume, then electric field due to it, in general, is given as:
1 dq
E=
4 0 r 2
rˆ
The charge is integrated over a line or surface or volume. Hence the electric
1 (r )dS
E (r ) =
4 0
surface
r 2
rˆ Surface Ch arg e Distribution
1 (r )dV
E (r ) =
4 0
volume
r 2
rˆ Volume Ch arg e Distribution
= E dS = EdS cos = ( E cos )dS
= En dS
The electric flux through an open surface is expressed as:
= E dS
surface
1
EdV = 0 dV
V
Maxwell’s
First
Equation
E =
0
𝑏
The line integral of the electric field along path 1 is 𝐸 𝑎. 𝑑𝑙
𝑎
The line integral of the electric field along path 1 is 𝐸 𝑏. 𝑑𝑙
𝑎 𝑞 𝑞 1 1
𝐸 𝑏. 𝑑𝑙 = 4𝜋𝜀
2 𝑟.Ƹ 𝑑𝑟𝑟Ƹ + 𝑟𝑑𝜃 𝜃 + 𝑟𝑠𝑖𝑛𝜃𝑑𝜙 𝜙 = −
0𝑟 4𝜋𝜀0 𝑟𝑏 𝑟𝑎
𝑏 𝑎
𝐸 ׯ. 𝑑𝑙 = 𝐸 × ∇ װ. 𝑑𝑎 = 𝐸 𝑎. 𝑑𝑙 + 𝐸 𝑏. 𝑑𝑙 = 0
⇒∇ ×𝐸 =0
𝜌
The equations ∇. 𝐸 = 𝜀 𝑎𝑛𝑑 ∇ × 𝐸 = 0 uniquely define the electric
0
field revealing that its source is charge density and is irrotational in nature.
Maxwell’s Second Equation
Biot – Savert’s Law & Divergence of Magnetic Field
The magnetic field intensity due to a steady state current (I), flowing in a
B = (
0 I dl
4
) ( )
r I − I r 2 = 0
r
through a closed surface is always zero.
• It implies that each magnet or source of
We know that :
( )
magnetic field enclosed contributes equal
I = 0 (current independent of r ) outward and inward magnetic flux.
( )
r = 0 (Cross product of two parallel vectors) • This is possible if a magnet or source of
magnetic field always exists as a dipole.
B = 0
Maxwell’s 2nd Equation : Qualitative Basis
Every magnet comprises of a north and a south pole.
If the magnet is divided into two pieces, it results in formation of two new
o outside the magnet, the field lines emerge out from the N-pole and
o inside the magnet, the field lines traverse from south to north pole.
will be equal to the flux penetrating into the surface thereby resulting in
overall zero magnetic flux.
The situation remains same even when many magnets are enclosed by
a common surface.
B dS = 0 Maxwell’s
Second
BdV = 0 Equation
B = 0
Maxwell’s Third Equation
Faraday’s Experiments
Michael Faraday in 1831 performed a series of experiments to
relate time variation in the magnetic flux through a circuit with the
emf developed across its ends.
enunciated:
d
e=−
dt
The Faraday’s law states that time variation of magnetic flux through a
d d
e=− = − B dS
dt dt
d
e = E dl = − B dS • It is time dependent equation.
dt
• This indicates that time variation
d
E dS = − dt B dS in magnetic field yields electric
field varying in space.
• This equation gives the desired
dB
E = − transformation of magnetic field
dt to electric field.
( ) ( )
E = −B = 0
t
From the above expression, the magnetic field can be expressed in
terms of another quantity called vector magnetic potential as 𝑩 = 𝜵 × 𝑨.
Hence Ampere’s circuital rule is valid for magnetostatics (i.e. where the
magnetic field due to current is uniform everywhere) only.
▪ Energetics in EM Waves.
Propagation of EM Waves in Dielectric Medium
Maxwell Equations in Perfect Dielectric Medium
A perfect dielectric medium is characterized by:
1
medium with the same velocity i.e 𝑣 = .
𝜀𝜇
1
If the medium is air, then speed is 𝑐 =
𝜀0 𝜇0
Given 𝜺𝟎 = 𝟖. 𝟖𝟓 × 𝟏𝟎−𝟏𝟐𝑭/𝒎
𝝁𝟎 = 𝟒𝝅 × 𝟏𝟎−𝟕 𝑯/𝒎
The electric and magnetic waves are inseparable from each other and
propagate together in space as EM wave.
➢ The time variation of magnetic field generates electric field
varying in space as per Faraday’s law of EM induction.
➢ Similarly time varying electric field produces space varying
magnetic field in accordance with Ampere-Maxwell law.
Transverse Nature of EM Waves
A plane wave is one in which the wave variables are constant in a
The x-components of both electric and magnetic fields do not show any
variation along x-direction (i.e. along the direction of propagation of EM wave)
and hence can be considered to be a constant quantity.
Using third Maxwell’s
equation and plane wave conditions
B H
E = − = −
t t
E z E y H x H x
− = − =0 (13)
y z t t
E x E z H y E z H y
− = − = (14)
z x t x t
E y E x H z E y H z
− = − = − (15)
x y t x t
y-direction.
direction and magnetic field vector oscillating along Two fields are interwoven
as well as inseparable as
z-direction.
the time variation in one
The solutions of differential equations (5) and (6) field generates the other
i (t − kx )
are:
E y = E0 e
field with space variation.
i (t − kx )
(19)
z
H = H 0 e
Relation of E and H Fields in EM Waves
When the EM waves propagate in free space or perfect
Ey E0 0
Z= = = 377
Hz H0 0
✓ When EM waves propagate in free
Transverse nature of light can also space or air, the electric field vector is
about 377 times more intense than
be expressed by generalized field
the magnetic field vector.
transformation equation as: ✓ In the wave theory of light, only
electric vector oscillations are
Bz = (E y )
k
considered due to this reason.
✓
( )
Even retina of our eye is more
1
ˆ
Bz = k E y (21) sensitive to electric vector oscillations
v
than those of magnetic field.
Propagation of EM Waves in Conducting Medium
EM Waves in Conducting Medium
The conducting medium is characterized by 𝜎 ≠ 0 𝑎𝑛𝑑 𝜌 = 0.
𝜕𝐸
but for the magnetic diffusivity term 𝜇𝜎 𝜕𝑡 , which arises as a result of the
For good conductors we have
The equation (30) is the progressive but damped electric wave in a conducting
medium. The amplitude of the wave decays as it penetrates deeper into the
conducting medium due to energy loss caused by Joule’s heating. This
phenomenon is termed as the skin effect.
The skin depth () is defined as the distance traversed by the EM waves,
at which its electric field is attenuated to 1/e of its value at the surface.
1 2
= = (31)
k
The skin depth is frequency dependent and the skin effect is severe at
the high frequencies of EM waves.
Anomalous Dispersion of EM Waves in a Conductor
The dispersion of EM waves occurs in a medium when their phase and
The electric field vector oscillations in the propagating electric waves in the
dv d dv v d r
r 2v + v 2 r = 0 =−
d d d 2 r d
dv d r
vg = v − = v 1 + (32)
d 2 r d
𝜕𝐸𝑦 𝜕𝐻𝑧
Two electric fields are related to each other as = −𝜇
𝜕𝑥 𝜕𝑡
The oscillating electric and magnetic vectors for EM wave propagation in the
H z
t
=
t
( )
H 0 e ( it −x ) = iH z
− E y = −iH z
Ey i
Z = =
Hz
Putting = (1 + i )
, we get :
2
2 i 2 1 + i 1 + i 1 i
Z= = = = + = cos + i sin
1+ i 2 2 2 2 4 4
i 4
Z= e (33) Impedance of Conductor
When EM waves propagate through a conductor, the electric and magnetic field oscillations
differ by a phase of π/4.
Response of Medium to EM Waves
The ratio of current density to displacement current density is:
J E y i2
= =− = e (34)
JD E y i
t
Material
Remarks
Conductor ~100 Conduction current dominates.
Semiconductor 0.01 to 100 -
Dielectric ~0.01 Displacement current dominates.
Energetics in EM Waves
Work Done (Energy) in the Electric Field
Total electric potential energy associated with discrete distribution of charges
1
W = 0 E dv
2
(35)
2 v
Work Done (Energy) in the Magnetic Field
The energy associated with flow of current (I) in the inductor (L) is given as:
1 2
W = LI
2
= LI = B da = ( A) da = A dl
1 1
W = ( LI ) I = I
2 2
1 1
W = I A dl = A ( Idl )
2 2
We know that Idl = Jdadl = Jdv
1
W = ( A J )dv
2
U sin g Maxwell fourth equation, we get :
1 1
W =
20 A ( B )dv =
20 ( A) Bdv
1
W = dv
2
B (36)
20
Concept of Poynting Vector
The cross-product of electric and magnetic field vectors is called
Poynting vector. Mathematically it is expressed as:
S = EH (37)
If we consider an EM wave propagating in x-direction with its electric
Sx = Ey H z
The Poynting vector denotes the time rate of flow of energy per unit area
of the medium through which the EM wave is propagating. It is often
S = c 0 E y2 = c 0 H z2 (39)
Average Value of Poynting Vector
The time average of the Poynting vector is expressed as:
1
T
E H
T
S av = iˆ E y H z dt = iˆ 0 0 (t − kx)dt
sin 2
T 0 T 0
E H
S av = iˆ 0 0
2
E H
S av = ˆj 0 kˆ 0
2 2
S av = Erms H rms
1 1
ˆ ˆ
S av = i c 0 E y + i c 0 H z2
2
(40)
2 2
P = S dA = S dV
V
S = ( E H ) = H ( E ) + E ( H )
S = H ( E ) − E ( H )
Using Maxwell’s third
and fourth equations,
we get:
B D
S = (− ) H − E (J + )
t t
B D
S = − H − J E − E
t t
we
Total power transmitted through closed surface is expressed as:
1
S dA = −
V
2 t
( D E + B H ) − J E dV
(41)
1 1 1
p = 2 S = 0 E02 cos 2 (kx − t + )iˆ p = uemiˆ (42)
c c c
1
p = 0 E02iˆ (43)
2c
1
I = S = c 0 E02 (44)
2
Laws of Optics
Reflection & Transmission of EM Waves : Normal Incidence
Let an EM wave be incident normally at the boundary of two media of
denoted by Ei, Er and Et and that for magnetic field vectors be denoted
by Hi, Hr and Ht respectively.
Ei Er
=− = Z1 (46)
Hi Hr
Et
= Z2 (47)
Ht
At normal incidence, the normal components of electric and magnetic field
Ei + Er = Et (48) Hi + Hr = Ht (49)
Using equations (46, 47) in equation (49), we get:
Ei Er Et Z
− = Ei − Er = 1 Et (50)
Z1 Z1 Z 2 Z2
From (48) and (50), we get:
Z1 + Z 2 Et Z − Z1 Et
Ei = (51) Er = 2 (52)
Z2 2 Z2 2
E r Z 2 − Z1
r= = (53)
Ei Z1 + Z 2
Et 2Z 2
t = = (54)
Ei Z1 + Z 2
Case I: EM waves are incident from air to conducting medium, then Z 1>>Z2.
Et Z
t= =2 2
Ei Z1
When EM waves are incident on conductor from air :
Et E
=2 i H t = 2H i (55) ✓ Incident electric vector suffers phase reversal on
Z2 Z1 reflection.
Er ✓ Transmitted magnetic field vector is twice in
r= = −1 Er = − Ei (56)
Ei magnitude as compared to the incident one.
Case II: EM waves propagate from conducting medium into air then Z1<< Z2.
Et 2Z 2
t= = =2 Et = 2 Ei (57)
Ei Z 2 When EM waves propagate from conductor to air :
Er Z 2 ✓ Reflected components of electric and magnetic
r= = =1 Er = Ei (58) fields remain same in magnitude without any
Ei Z 2
phase reversal.
E E
Also r = i Hr = Hi (59) ✓ Transmitted electric vector is twice in magnitude
Z1 Z1 as compared to incident one.
Relation Between Impedance & Refractive Index
The refractive index of the medium is expressed as:
c 1
n= = . = . = r r
v 0 0 0 0
r = 1 for non − magnetic subs tan ces
n = r (60)
If Z0 and Z are impedance of free space and medium, then we get:
0 Zo
Z0 = Z= = r = n (61)
0 Z
The reflection and transmission coefficients of amplitude are given as:
Z −Z 2 Z1
r= 2 1 t=
Z 2 + Z1 Z 2 + Z1
Let Z 2 = Z and Z1 = Z 0 then
1− n
r= (62)
1+ n
2
t= (63)
1+ n
The reflection coefficients of energy is given as:
Si EH
R= = i i
Sr Er H r
0 0
H i = Ei H r = Er
0 0
2
E 1− n
2
R = r = r 2 =
Ei 1 + n
2
Er 1− n
2
R = = r 2 = (64)
i
E 1 + n
The transmission coefficients of energy is given as:
St EH
T= = t t
Si Ei H i
1− n
2
4n
T = 1− R = 1− =
1+ n (1 + n )2
4n
T = (65)
(1 + n) 2
Reflection & Transmission of EM Waves: Oblique Incidence
Let’s consider two media of refractive indices n1 and n2 having interface at XY
Let the angles made by reflected and transmitted waves with z-axis be θr and
θt respectively.
ki v1 = k r v1 = kt v2 = (70)
More explicitly, the condition given by equation (69) can be expressed as:
1 E1⊥ = 2 E2⊥
II
1
E = E II
2
⊥ ⊥
B1 = B2 (75)
1 B II = 1 B II
1
2
1 2
Applying boundary conditions given in equations (76) using the equations
(66-68), we get:
1 (E0i + E0 r )z = 2 (E0t )z
( E
0i + E )
0r x, y = ( E )
0t x , y
(B0i + B0 r )z = (B0t )z (76)
(B + B ) =
1 1
(B0t )x , y
1
0i 0r x, y
1
Let’s consider that the polarization of incident wave is parallel to the plane of
incidence (i.e. XZ plane) and reflected and transmitted waves are also
polarized in the same plane.
form as E0r=E0i.
Case II: Let’s consider intermediate angle θb where the reflected wave is
completely extinguished.
n2
tan b = (93) Brewester ' s Law
n1
Numerical Problems
Show that in EM waves, electrostatic energy density is equal to magnetic
energy density.
Show that the energy flux in a plane EM polarized wave propagating in free
(0.24oC/s)
2
Show that the skin depth in poor conductor is and hence independent
of frequency of EM waves.
If alternating electric field 𝐸 = 𝐸0 𝑐𝑜𝑠𝜔𝑡 is applied to a conductor, show that
the displacement current is negligible as compared to conduction current at
any frequency lower than the optical frequencies.
Calculate the skin depth for silver at afrequency of 100GHz. Given that
σ=2x107S/m and μ=4πx107.
Show that Sav = 1.327 x 10-3Emax Jm2s-1 if given that Emax = 27.45√Sav V/m
Calculate the poynting vector at the surface of sun. Given that it radiates
3.8 x 1026W and radius of sun is 0.7Gm. Also evaluate the electric and
magnetic field vector amplitudes. (1340W/m2, 10004V/m, 2.27A/m).
Starting with continuity equation and assuming ohm’s law, show that the charge
𝜎𝜚 𝜕𝜚
density obeys the equation + = 0. Solve this equation and discuss the
𝜀 𝜕𝑡
The EM wave has the electric and magnetic field vectors in the conductor
𝜎
expressed as 𝐸𝑥 = 𝐴𝑒 −𝑘𝑧 𝑒 𝑖(𝜔𝑡−𝑘𝑧) , 𝐻𝑦 = 𝐴 𝑒 −𝑘𝑧 𝑒 𝑖(𝜔𝑡−𝑘𝑧) 𝑒 𝑖𝜋/4 𝑤ℎ𝑒𝑟𝑒 𝑘 =
𝜔𝜇
𝜔𝜇𝜎
. Show that the average value of Poynting vector in the conductor is 𝑆𝑎𝑣 =
2
𝐴2 𝜎
𝑒 −2𝑘𝑧 .
2 2𝜇𝜔