Skip to content

Electromagnetic Waves

In free space (ρ=0\rho = 0, J=0\mathbf{J} = \mathbf{0}), take the curl of Faraday’s law:

×(×E)=t(×B)=μ0ε02Et2\nabla \times (\nabla \times \mathbf{E}) = -\frac{\partial}{\partial t}(\nabla \times \mathbf{B}) = -\mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2}

Using the identity ×(×E)=(E)2E\nabla \times (\nabla \times \mathbf{E}) = \nabla(\nabla \cdot \mathbf{E}) - \nabla^2 \mathbf{E} And E=0\nabla \cdot \mathbf{E} = 0:

2E=μ0ε02Et2\nabla^2 \mathbf{E} = \mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2}

Similarly: 2B=μ0ε02Bt2\nabla^2 \mathbf{B} = \mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{B}}{\partial t^2}.

These are wave equations with wave speed c=1/μ0ε03×108c = 1/\sqrt{\mu_0 \varepsilon_0} \approx 3 \times 10^8 m/s.

Theorem 5.1. Electromagnetic waves in free space are:

  1. Transverse: E\mathbf{E} and B\mathbf{B} are perpendicular to the direction of propagation.
  2. Mutually perpendicular: EB\mathbf{E} \perp \mathbf{B}.
  3. In phase: E=cBE = cB at every point.
  4. Linearly polarised (; other polarisations are superpositions).

Energy. The energy density of an EM wave is u=12(ε0E2+B2/μ0)u = \frac{1}{2}(\varepsilon_0 E^2 + B^2/\mu_0).

The Poynting vector S=1μ0E×B\mathbf{S} = \frac{1}{\mu_0}\mathbf{E} \times \mathbf{B} represents the energy Flux (power per unit area).

Problem. Show that E=E0cos(kzωt)x^\mathbf{E} = E_0 \cos(kz - \omega t)\,\hat{\mathbf{x}} satisfies the wave Equation and find the associated B\mathbf{B} field.

Solution. 2E=2Exz2x^=k2E0cos(kzωt)x^\nabla^2 \mathbf{E} = \frac{\partial^2 E_x}{\partial z^2}\hat{\mathbf{x}} = -k^2 E_0 \cos(kz - \omega t)\,\hat{\mathbf{x}}.

2Et2=ω2E0cos(kzωt)x^\frac{\partial^2 \mathbf{E}}{\partial t^2} = -\omega^2 E_0 \cos(kz - \omega t)\,\hat{\mathbf{x}}.

The wave equation requires k2=μ0ε0ω2k^2 = \mu_0 \varepsilon_0 \omega^2I.e., ω/k=c\omega/k = c.

From Faraday’s law: ×E=Bt\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}.

(×E)y=Exz=kE0sin(kzωt)(\nabla \times \mathbf{E})_y = -\frac{\partial E_x}{\partial z} = k E_0 \sin(kz - \omega t)

Byt=kE0sin(kzωt)    By=kωE0cos(kzωt)=E0ccos(kzωt)\frac{\partial B_y}{\partial t} = -k E_0 \sin(kz - \omega t) \implies B_y = \frac{k}{\omega} E_0 \cos(kz - \omega t) = \frac{E_0}{c}\cos(kz - \omega t)

So B=E0ccos(kzωt)y^\mathbf{B} = \frac{E_0}{c}\cos(kz - \omega t)\,\hat{\mathbf{y}}. \blacksquare

5.4 Poynting’s Theorem and Energy Conservation

Section titled “5.4 Poynting’s Theorem and Energy Conservation”

Poynting’s theorem is the statement of energy conservation for electromagnetic fields.

Derivation. Start with the two Maxwell equations containing time derivatives:

×E=Bt,×B=μ0J+μ0ε0Et\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}, \quad \nabla \times \mathbf{B} = \mu_0 \mathbf{J} + \mu_0 \varepsilon_0 \frac{\partial \mathbf{E}}{\partial t}

Compute B(×E)E(×B)\mathbf{B} \cdot (\nabla \times \mathbf{E}) - \mathbf{E} \cdot (\nabla \times \mathbf{B}):

B(×E)=BBt=t(B22)\mathbf{B} \cdot (\nabla \times \mathbf{E}) = -\mathbf{B} \cdot \frac{\partial \mathbf{B}}{\partial t} = -\frac{\partial}{\partial t}\left(\frac{B^2}{2}\right)

E(×B)=μ0EJμ0ε0EEt=μ0EJt(ε0E22)-\mathbf{E} \cdot (\nabla \times \mathbf{B}) = -\mu_0 \mathbf{E} \cdot \mathbf{J} - \mu_0 \varepsilon_0 \mathbf{E} \cdot \frac{\partial \mathbf{E}}{\partial t} = -\mu_0 \mathbf{E} \cdot \mathbf{J} - \frac{\partial}{\partial t}\left(\frac{\varepsilon_0 E^2}{2}\right)

Using the vector identity (E×B)=B(×E)E(×B)\nabla \cdot (\mathbf{E} \times \mathbf{B}) = \mathbf{B} \cdot (\nabla \times \mathbf{E}) - \mathbf{E} \cdot (\nabla \times \mathbf{B}):

(E×B)=μ0JEμ0ε0t(E22)t(B22)\nabla \cdot (\mathbf{E} \times \mathbf{B}) = -\mu_0 \mathbf{J} \cdot \mathbf{E} - \mu_0 \varepsilon_0 \frac{\partial}{\partial t}\left(\frac{E^2}{2}\right) - \frac{\partial}{\partial t}\left(\frac{B^2}{2}\right)

Dividing by μ0\mu_0 and rearranging:

S=JE+ut\boxed{-\nabla \cdot \mathbf{S} = \mathbf{J} \cdot \mathbf{E} + \frac{\partial u}{\partial t}}

Where S=1μ0E×B\mathbf{S} = \frac{1}{\mu_0}\mathbf{E} \times \mathbf{B} is the Poynting vector and u=12(ε0E2+B2μ0)u = \frac{1}{2}\left(\varepsilon_0 E^2 + \frac{B^2}{\mu_0}\right) is the energy density.

Interpretation: The rate of energy leaving a volume equals the work done on charges plus The rate of increase of field energy. In integral form:

SSdA=ddtVudV+VJEdV-\oint_S \mathbf{S} \cdot d\mathbf{A} = \frac{d}{dt}\int_V u\,dV + \int_V \mathbf{J} \cdot \mathbf{E}\,dV

Intensity. For a plane wave, the time-averaged Poynting vector is:

S=E022μ0ck^=12ε0cE02k^\langle\mathbf{S}\rangle = \frac{E_0^2}{2\mu_0 c}\,\hat{\mathbf{k}} = \frac{1}{2}\varepsilon_0 c E_0^2\,\hat{\mathbf{k}}

Example: Radiation pressure

A plane wave normally incident on a perfectly absorbing surface exerts a radiation pressure. The momentum flux of the wave is S/c\langle S \rangle/c per unit area, so:

Pabs=Sc=ε0E022P_{\mathrm{abs} = \frac{\langle S \rangle}{c} = \frac{\varepsilon_0 E_0^2}{2}}

For a perfectly reflecting surface, the momentum transfer is doubled:

Pref=2Sc=ε0E02P_{\mathrm{ref} = \frac{2\langle S \rangle}{c} = \varepsilon_0 E_0^2}

A 1 kW/m2^2 beam (like sunlight near Earth) exerts a pressure of about 3.3 μ3.3\ \muPa on a Perfect absorber. \blacksquare

Example: Polarization of EM waves

Linear polarization. E=E0cos(kzωt)x^\mathbf{E} = E_0\cos(kz - \omega t)\,\hat{\mathbf{x}}. The field Oscillates in a fixed direction.

Circular polarization. Two orthogonal linear polarizations with a phase difference of π/2\pi/2:

E=E0cos(kzωt)x^±E0sin(kzωt)y^\mathbf{E} = E_0\cos(kz - \omega t)\,\hat{\mathbf{x}} \pm E_0\sin(kz - \omega t)\,\hat{\mathbf{y}}

The tip of E\mathbf{E} traces a circle. The ++ sign gives left-circular polarization (LCP) and the - sign gives right-circular polarization (RCP).

Elliptical polarization. The general case with arbitrary amplitudes and phase:

E=E0xcos(kzωt)x^+E0ycos(kzωt+δ)y^\mathbf{E} = E_{0x}\cos(kz - \omega t)\,\hat{\mathbf{x}} + E_{0y}\cos(kz - \omega t + \delta)\,\hat{\mathbf{y}}

\blacksquare

In a conductor with conductivity σ\sigmaOhm’s law gives J=σE\mathbf{J} = \sigma\mathbf{E}. Substituting into the Ampere-Maxwell law:

×B=μ0σE+μ0ε0Et\nabla \times \mathbf{B} = \mu_0\sigma\mathbf{E} + \mu_0\varepsilon_0\frac{\partial \mathbf{E}}{\partial t}

For a monochromatic wave E=E0eiωt\mathbf{E} = \mathbf{E}_0\,e^{-i\omega t}This leads to a complex Wave number:

k~2=μ0ε0ω2+iμ0σω\tilde{k}^2 = \mu_0\varepsilon_0\omega^2 + i\mu_0\sigma\omega

Writing k~=k+iκ\tilde{k} = k + i\kappa where kk is the real part (wave number) and κ\kappa is the Imaginary part (attenuation constant):

E(z,t)=E0eκzcos(kzωt)\mathbf{E}(z,t) = \mathbf{E}_0\,e^{-\kappa z}\cos(kz - \omega t)

The field decays exponentially. The skin depth is the distance over which the amplitude Falls by a factor of 1/e1/e:

δ=1κ\delta = \frac{1}{\kappa}

For a good conductor (σε0ω\sigma \gg \varepsilon_0\omega):

δ=2μ0σω\delta = \sqrt{\frac{2}{\mu_0\sigma\omega}}

Example: Skin depth in copper at 60 Hz and 1 MHz

Copper: σ=5.96×107\sigma = 5.96 \times 10^7 S/m, μr1\mu_r \approx 1.

At f=60f = 60 Hz (ω=2π×60\omega = 2\pi \times 60 rad/s):

δ=24π×107×5.96×107×2π×608.5 mm\delta = \sqrt{\frac{2}{4\pi \times 10^{-7} \times 5.96 \times 10^7 \times 2\pi \times 60}} \approx 8.5\ \mathrm{mm}

At f=1f = 1 MHz (ω=2π×106\omega = 2\pi \times 10^6 rad/s):

δ=24π×107×5.96×107×2π×10665 μm\delta = \sqrt{\frac{2}{4\pi \times 10^{-7} \times 5.96 \times 10^7 \times 2\pi \times 10^6}} \approx 65\ \mu\mathrm{m}

The skin depth decreases as 1/f1/\sqrt{f}So higher-frequency signals are confined to thinner Surface layers. \blacksquare

Electromagnetic waves can be guided by hollow conducting pipes (waveguides). Consider a Rectangular waveguide with dimensions aa (width) and bb (height).

TE modes (transverse electric, Ez=0E_z = 0, Bz0B_z \neq 0). The lowest-order mode is \mathrm{TE_}{10}With fields:

Ey=E0sin ⁣(πxa)cos(kgzωt)E_y = E_0 \sin\!\left(\frac{\pi x}{a}\right)\cos(k_g z - \omega t)

Bx=kgωE0sin ⁣(πxa)cos(kgzωt)B_x = -\frac{k_g}{\omega}E_0 \sin\!\left(\frac{\pi x}{a}\right)\cos(k_g z - \omega t)

Bz=πωaE0cos ⁣(πxa)sin(kgzωt)B_z = \frac{\pi}{\omega a}E_0 \cos\!\left(\frac{\pi x}{a}\right)\sin(k_g z - \omega t)

Where the guide wave number is kg=(ω/c)2(π/a)2k_g = \sqrt{(\omega/c)^2 - (\pi/a)^2}.

Cutoff frequency. Waves propagate only when ω>ωc\omega \gt \omega_c where:

ωc,mn=cπ(ma)2+(nb)2\omega_{c,mn} = c\pi\sqrt{\left(\frac{m}{a}\right)^2 + \left(\frac{n}{b}\right)^2}

For the \mathrm{TE_}{10} mode: fc=c2af_c = \frac{c}{2a}.

Phase and group velocities. In a waveguide, the phase velocity exceeds cc:

vp=ωkg=c1(ωc/ω)2>cv_p = \frac{\omega}{k_g} = \frac{c}{\sqrt{1 - (\omega_c/\omega)^2}} \gt c

The group velocity (signal velocity) is less than cc:

vg=dωdkg=c1(ωc/ω)2<cv_g = \frac{d\omega}{dk_g} = c\sqrt{1 - (\omega_c/\omega)^2} \lt c

They satisfy vpvg=c2v_p\,v_g = c^2.

:::caution Common Pitfall The phase velocity in a waveguide exceeds ccBut this does not violate special relativity. No information or energy travels faster than cc; the signal velocity is the group velocity vg<cv_g \lt c. The phase velocity is the speed of the wave crests, which is a purely Kinematic quantity.

An oscillating electric dipole is the simplest source of electromagnetic radiation.

Consider a dipole p(t)=p0cos(ωt)z^\mathbf{p}(t) = p_0\cos(\omega t)\,\hat{\mathbf{z}}. In the radiation zone (rλr \gg \lambda), the fields are:

E=μ0p0ω24πsinθrcos[ω(tr/c)]θ^\mathbf{E} = -\frac{\mu_0 p_0 \omega^2}{4\pi}\frac{\sin\theta}{r}\cos[\omega(t - r/c)]\,\hat{\boldsymbol{\theta}}

B=μ0p0ω24πcsinθrcos[ω(tr/c)]ϕ^\mathbf{B} = -\frac{\mu_0 p_0 \omega^2}{4\pi c}\frac{\sin\theta}{r}\cos[\omega(t - r/c)]\,\hat{\boldsymbol{\phi}}

The fields fall off as 1/r1/r (not 1/r21/r^2 as for static fields), which is characteristic of Radiation.

Radiation pattern. The intensity varies as sin2θ\sin^2\thetaWith maximum radiation in the Equatorial plane (θ=π/2\theta = \pi/2) and zero along the dipole axis (θ=0,π\theta = 0, \pi).

Total radiated power. Integrating the Poynting vector over a sphere:

P=μ0p02ω412πcP = \frac{\mu_0 p_0^2 \omega^4}{12\pi c}

Larmor formula. For a point charge qq undergoing acceleration aa:

P=q2a26πε0c3P = \frac{q^2 a^2}{6\pi\varepsilon_0 c^3}

This is the non-relativistic limit and is valid whenever vcv \ll c.

Derivation: Power radiated by an oscillating dipole

The time-averaged Poynting vector magnitude in the radiation zone:

S=12μ0EθBϕ=μ0p02ω432π2csin2θr2\langle S \rangle = \frac{1}{2\mu_0}\lvert E_\theta\rvert\,\lvert B_\phi\rvert = \frac{\mu_0 p_0^2\omega^4}{32\pi^2 c}\frac{\sin^2\theta}{r^2}

The total power through a sphere of radius rr:

P=02π ⁣ ⁣0πSr2sinθdθdϕ=μ0p02ω432π2c2π0πsin3θdθP = \int_0^{2\pi}\!\!\int_0^\pi \langle S \rangle\, r^2\sin\theta\,d\theta\,d\phi = \frac{\mu_0 p_0^2\omega^4}{32\pi^2 c} \cdot 2\pi \int_0^\pi \sin^3\theta\,d\theta

Using 0πsin3θdθ=4/3\int_0^\pi \sin^3\theta\,d\theta = 4/3:

P=μ0p02ω432π2c2π43=μ0p02ω412πcP = \frac{\mu_0 p_0^2\omega^4}{32\pi^2 c} \cdot 2\pi \cdot \frac{4}{3} = \frac{\mu_0 p_0^2\omega^4}{12\pi c}

\blacksquare

:::