Skip to content

Advanced Electrodynamics

11.1 Multipole Expansion

The scalar potential of a localised charge distribution at large distances (rdr \gg dWhere dd is the size of the distribution):

ϕ(r)=14πε0[Qr+pr^r2+12ijQijr^ir^jr3+]\phi(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0}\left[\frac{Q}{r} + \frac{\mathbf{p}\cdot\hat{\mathbf{r}}}{r^2} + \frac{1}{2}\sum_{ij}Q_{ij}\frac{\hat{r}_i\hat{r}_j}{r^3} + \cdots\right]

Monopole term: Q=ρdVQ = \int \rho\, dV (total charge).

Dipole term: p=r"ρ(r)dV\mathbf{p} = \int \mathbf{r}"\,\rho(\mathbf{r}')\,dV' (electric dipole moment).

Quadrupole term: Qij=(3rirjr2δij)ρ(r)dVQ_{ij} = \int (3r'_ir'_j - r'^2\delta_{ij})\,\rho(\mathbf{r}')\,dV' (traceless quadrupole tensor).

The quadrupole term is important for nuclei with spin I1I \geq 1 and for non-spherical charge distributions. The quadrupole moment Q=QzzQ = Q_{zz} (in the principal axis frame) characterises the deviation from spherical symmetry.

11.2 Gauge Transformations and Potentials

The scalar and vector potentials are not unique. The gauge transformation:

A=A+χ,ϕ=ϕχt\mathbf{A}' = \mathbf{A} + \nabla\chi, \quad \phi' = \phi - \frac{\partial\chi}{\partial t}

Leaves E\mathbf{E} and B\mathbf{B} unchanged for any scalar function χ(r,t)\chi(\mathbf{r}, t).

Common gauges:

GaugeConditionUse
CoulombA=0\nabla \cdot \mathbf{A} = 0Static problems, quantum mechanics
LorenzA+1c2ϕt=0\nabla \cdot \mathbf{A} + \frac{1}{c^2}\frac{\partial\phi}{\partial t} = 0Relativistic problems, radiation
Temporalϕ=0\phi = 0Some scattering problems

In the Lorenz gauge, both A\mathbf{A} and ϕ\phi satisfy wave equations with sources:

2A=μ0J,2ϕ=ρε0\Box^2\mathbf{A} = -\mu_0\mathbf{J}, \quad \Box^2\phi = -\frac{\rho}{\varepsilon_0}

Where 2=21c22t2\Box^2 = \nabla^2 - \frac{1}{c^2}\frac{\partial^2}{\partial t^2} is the d’Alembertian.

11.3 Electromagnetic Stress-Energy Tensor

The electromagnetic stress-energy tensor TμνT^{\mu\nu} encodes the energy density, momentum density, and stress:

T00=12(ε0E2+B2μ0)(energy density)T^{00} = \frac{1}{2}\left(\varepsilon_0 E^2 + \frac{B^2}{\mu_0}\right) \quad \text{(energy density)}

T0i=1c(E×B)i=Sic(momentum density)T^{0i} = \frac{1}{c}(\mathbf{E} \times \mathbf{B})_i = \frac{S_i}{c} \quad \text{(momentum density)}

Tij=ε0EiEj1μ0BiBj+12δij(ε0E2+B2μ0)(Maxwell stress tensor)T^{ij} = -\varepsilon_0 E_i E_j - \frac{1}{\mu_0}B_i B_j + \frac{1}{2}\delta_{ij}\left(\varepsilon_0 E^2 + \frac{B^2}{\mu_0}\right) \quad \text{(Maxwell stress tensor)}

Conservation law: μTμν=fν\partial_\mu T^{\mu\nu} = -f^\nu where fνf^\nu is the Lorentz force density on charges.

Radiation pressure: For a normally incident plane wave with intensity II:

Prad=Ic=12ε0E02P_{\text{rad} = \frac{I}{c} = \frac{1}{2}\varepsilon_0 E_0^2}

For a perfect reflector, the radiation pressure is 2I/c2I/c (momentum transfer is doubled).

Worked Example 11.1: Radiation Pressure from Sunlight

Solar constant at Earth: I=1361I = 1361 W/m2^2.

Radiation pressure on a perfectly absorbing surface:

P=Ic=13613×108=4.54×106 N/m2=4.54 μPaP = \frac{I}{c} = \frac{1361}{3 \times 10^8} = 4.54 \times 10^{-6}\ \text{N}/m^2 = 4.54\ \mu\text{Pa}

For a perfect reflector: P=9.07μPaP = 9.07\,\mu\text{Pa}.

This is tiny compared to atmospheric pressure (10510^5 Pa), but is significant for:

  • Solar sails: A 100 m ×\times 100 m sail with 90% reflectivity experiences F0.12F \approx 0.12 N, producing acceleration a0.6a \approx 0.6 mm/s2^2 for a 100 kg sail.
  • Asteroid deflection: Sustained radiation pressure can perturb asteroid orbits over years.
  • Atom optics: Laser cooling uses radiation pressure to slow atoms to microkelvin temperatures.

Common Pitfalls (Additional)

  1. vp>cv_p > c does not violate relativity: The phase velocity in a waveguide exceeds ccBut no information or energy travels faster than cc. The group velocity (signal velocity) is always vg<cv_g < c. Similarly, the refracted phase front in a prism can appear to move faster than ccBut the actual signal does not.

  2. Gauge choice matters for potentials, not fields: Different gauges give different A\mathbf{A} and ϕ\phi for the same E\mathbf{E} and B\mathbf{B}. In quantum mechanics, the Hamiltonian depends on the gauge, but all physical observables are gauge-invariant. The Aharonov—Bohm effect shows that even in regions where E=B=0\mathbf{E} = \mathbf{B} = 0The vector potential A\mathbf{A} has measurable physical effects.

  3. Multipole expansion convergence: The multipole expansion converges only outside a sphere that encloses all charges. Inside the charge distribution, the expansion diverges and must not be used. The expansion parameter is d/rd/r where dd is the source size and rr is the observation distance.

  4. Radiation fields vs. Near fields: At distances rλr \ll \lambda (near field), the fields are dominated by 1/r21/r^2 (induction) and 1/r31/r^3 (electrostatic/magnetostatic) terms. The radiation fields (1/r\propto 1/r) dominate only in the far field (rλr \gg \lambda). Do not apply the Larmor formula or radiation resistance in the near field.

  5. Poynting vector is not unique: The Poynting vector S=E×H\mathbf{S} = \mathbf{E} \times \mathbf{H} is gauge-dependent and can be nonzero even in static situations (e.g., a charged capacitor in a constant magnetic field). Only the surface integral Sda\oint \mathbf{S} \cdot d\mathbf{a} (total power flow) is physically meaningful.

Problems (Additional)

Problem 19: TE$_{10}$ Mode Field Patterns

For a rectangular waveguide (a×ba \times b) operating in TE10_{10} mode at frequency ff:

(a) Write the complete expressions for all six field components (Ex,Ey,Ez,Bx,By,BzE_x, E_y, E_z, B_x, B_y, B_z).

(b) Sketch the field pattern: show the direction and relative magnitude of E\mathbf{E} and B\mathbf{B} in the xyxy-plane at z=0z = 0.

(c) Find the positions of maximum surface current density on the walls and explain why the waveguide loss is minimised by making the broad wall dimension aa as large as possible (for a given ff).

Solution:

(a) For TE10_{10}: Bz=B0cos(πx/a)ei(βzωt)B_z = B_0\cos(\pi x/a)\,e^{i(\beta z - \omega t)}.

Ex=0,Ey=iωμ0aπB0sin ⁣(πxa)ei(βzωt)E_x = 0, \quad E_y = \frac{i\omega\mu_0 a}{\pi}B_0\sin\!\left(\frac{\pi x}{a}\right)e^{i(\beta z - \omega t)}

Ez=0E_z = 0

Bx=iβaπB0sin ⁣(πxa)ei(βzωt),By=0B_x = -\frac{i\beta a}{\pi}B_0\sin\!\left(\frac{\pi x}{a}\right)e^{i(\beta z - \omega t)}, \quad B_y = 0

Bz=B0cos ⁣(πxa)ei(βzωt)B_z = B_0\cos\!\left(\frac{\pi x}{a}\right)e^{i(\beta z - \omega t)}

(b) The electric field EyE_y is purely vertical, with a sin(πx/a)\sin(\pi x/a) profile: zero at the side walls (x=0,ax = 0, a) and maximum at the centre (x=a/2x = a/2). The magnetic field forms closed loops in the xzxz-plane.

(c) Surface current K=n^×H\mathbf{K} = \hat{\mathbf{n}} \times \mathbf{H}. On the broad walls (y=0,by = 0, b): K\mathbf{K} has components from BxB_x and BzB_zWith maximum at x=a/2x = a/2 (where sin(πx/a)=1\sin(\pi x/a) = 1). The power loss per unit length is:

Ploss=Rs2K2dlP_{\text{loss} = \frac{R_s}{2}\oint |\mathbf{K}|^2\, dl}

Where Rs=ωμ0/(2σ)R_s = \sqrt{\omega\mu_0/(2\sigma)} is the surface resistance. For fixed ffIncreasing aa reduces the current density on the broad walls and increases the power-handling capacity.

Problem 20: Antenna Radiation Pattern

A half-wave dipole antenna of length =λ/2\ell = \lambda/2 carries a sinusoidal current distribution:

I(z)=I0cos(kz),λ/4zλ/4I(z) = I_0\cos(kz), \quad -\lambda/4 \leq z \leq \lambda/4

(a) Calculate the radiation fields E\mathbf{E} and B\mathbf{B} in the far field.

(b) Find the angular distribution of radiated power dP/dΩdP/d\Omega.

(c) Calculate the total radiated power and the radiation resistance. Compare with the short-dipole result Rrad=197(/λ)2ΩR_{\text{rad} = 197(\ell/\lambda)^2\,\Omega}.

Solution:

(a) The vector potential in the far field:

Az=μ04πeikrrλ/4λ/4I0cos(kz)eikzcosθdzA_z = \frac{\mu_0}{4\pi}\frac{e^{ikr}}{r}\int_{-\lambda/4}^{\lambda/4}I_0\cos(kz')\,e^{-ikz'\cos\theta}\,dz'

The integral evaluates to:

Az=μ0I04πeikrr2cos ⁣(π2cosθ)ksin2θA_z = \frac{\mu_0 I_0}{4\pi}\frac{e^{ikr}}{r}\frac{2\cos\!\left(\frac{\pi}{2}\cos\theta\right)}{k\sin^2\theta}

The radiation fields:

Eθ=ikAzsinθ=iμ0cI04πeikrrcos ⁣(π2cosθ)sinθE_\theta = ikA_z\sin\theta = \frac{i\mu_0 c I_0}{4\pi}\frac{e^{ikr}}{r}\frac{\cos\!\left(\frac{\pi}{2}\cos\theta\right)}{\sin\theta}

Bϕ=Eθ/cB_\phi = E_\theta/c

(b) The angular distribution:

dPdΩ=r22Z0Eθ2=Z0I0232π2cos2 ⁣(π2cosθ)sin2θ\frac{dP}{d\Omega} = \frac{r^2}{2Z_0}|E_\theta|^2 = \frac{Z_0 I_0^2}{32\pi^2}\frac{\cos^2\!\left(\frac{\pi}{2}\cos\theta\right)}{\sin^2\theta}

(c) Total power:

P=Z0I0232π202π ⁣ ⁣0πcos2 ⁣(π2cosθ)sin2θsinθdθdϕP = \frac{Z_0 I_0^2}{32\pi^2}\int_0^{2\pi}\!\!\int_0^\pi \frac{\cos^2\!\left(\frac{\pi}{2}\cos\theta\right)}{\sin^2\theta}\sin\theta\,d\theta\,d\phi

=Z0I0216π0πcos2 ⁣(π2cosθ)sinθdθ= \frac{Z_0 I_0^2}{16\pi}\int_0^\pi \frac{\cos^2\!\left(\frac{\pi}{2}\cos\theta\right)}{\sin\theta}\,d\theta

With the substitution u=cosθu = \cos\theta: 11cos2(πu/2)1u2du=1.2188\int_{-1}^{1}\frac{\cos^2(\pi u/2)}{1-u^2}\,du = 1.2188 (the Siegel integral).

P=377×1.218816πI02=9.16I02P = \frac{377 \times 1.2188}{16\pi}I_0^2 = 9.16\,I_0^2

Radiation resistance: Rrad=2P/I02=18.3ΩR_{\text{rad} = 2P/I_0^2 = 18.3\,\Omega}.

For comparison, a short dipole (λ\ell \ll \lambda) of length λ/2\lambda/2 would give Rrad=197×0.25=49.3ΩR_{\text{rad} = 197 \times 0.25 = 49.3\,\Omega}. The half-wave dipole has lower radiation resistance because the current distribution (cosine) has less total effective acceleration than a uniform current.

The directivity of the half-wave dipole is D=1.64D = 1.64 (2.15 dBi), slightly higher than the short dipole (D=1.5D = 1.5).