Skip to main content

Electromagnetism

1. Maxwell's Equations

1.1 The Four Equations

Maxwell's equations are the foundation of classical electromagnetism. In SI units:

Integral Form:

SEdA=Qencε0(GausssLaw)\oint_S \mathbf{E} \cdot d\mathbf{A} = \frac{Q_{\mathrm{enc}}}{\varepsilon_0} \quad \mathrm{(Gauss's Law)}

SBdA=0(GausssLawforMagnetism)\oint_S \mathbf{B} \cdot d\mathbf{A} = 0 \quad \mathrm{(Gauss's Law for Magnetism)}

CEdl=dΦBdt(FaradaysLaw)\oint_C \mathbf{E} \cdot d\mathbf{l} = -\frac{d\Phi_B}{dt} \quad \mathrm{(Faraday's Law)}

CBdl=μ0Ienc+μ0ε0dΦEdt(AmpereMaxwellLaw)\oint_C \mathbf{B} \cdot d\mathbf{l} = \mu_0 I_{\mathrm{enc}} + \mu_0 \varepsilon_0 \frac{d\Phi_E}{dt} \quad \mathrm{(Ampere-Maxwell Law)}

Differential Form:

E=ρε0(GausssLaw)\nabla \cdot \mathbf{E} = \frac{\rho}{\varepsilon_0} \quad \mathrm{(Gauss's Law)}

B=0(GausssLawforMagnetism)\nabla \cdot \mathbf{B} = 0 \quad \mathrm{(Gauss's Law for Magnetism)}

×E=Bt(FaradaysLaw)\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t} \quad \mathrm{(Faraday's Law)}

×B=μ0J+μ0ε0Et(AmpereMaxwellLaw)\nabla \times \mathbf{B} = \mu_0 \mathbf{J} + \mu_0 \varepsilon_0 \frac{\partial \mathbf{E}}{\partial t} \quad \mathrm{(Ampere-Maxwell Law)}

where ρ\rho is the charge density, J\mathbf{J} is the current density, ε0\varepsilon_0 is the permittivity of free space, and μ0\mu_0 is the permeability of free space.

1.2 Derivation from Integral to Differential Form

Gauss's Law. Apply the divergence theorem to the integral form:

SEdA=V(E)dV=1ε0VρdV\oint_S \mathbf{E} \cdot d\mathbf{A} = \int_V (\nabla \cdot \mathbf{E})\, dV = \frac{1}{\varepsilon_0}\int_V \rho\, dV

Since this holds for any volume VV: E=ρ/ε0\nabla \cdot \mathbf{E} = \rho / \varepsilon_0.

Faraday's Law. Apply Stokes' theorem:

CEdl=S(×E)dA=SBtdA\oint_C \mathbf{E} \cdot d\mathbf{l} = \int_S (\nabla \times \mathbf{E}) \cdot d\mathbf{A} = -\int_S \frac{\partial \mathbf{B}}{\partial t} \cdot d\mathbf{A}

Since this holds for any surface SS: ×E=B/t\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t.

1.3 Continuity Equation

Taking the divergence of the Ampere-Maxwell law:

(×B)=0=μ0J+μ0ε0t(E)\nabla \cdot (\nabla \times \mathbf{B}) = 0 = \mu_0 \nabla \cdot \mathbf{J} + \mu_0 \varepsilon_0 \frac{\partial}{\partial t}(\nabla \cdot \mathbf{E})

Using Gauss's law: J+ρt=0\nabla \cdot \mathbf{J} + \frac{\partial \rho}{\partial t} = 0.

This is the continuity equation, expressing conservation of charge.

2. Electrostatics

2.1 Coulomb's Law and the Electric Field

Coulomb's Law: The force between two point charges q1q_1 and q2q_2 separated by distance rr:

F=14πε0q1q2r2r^\mathbf{F} = \frac{1}{4\pi\varepsilon_0} \frac{q_1 q_2}{r^2} \hat{\mathbf{r}}

The electric field due to a point charge qq at position r\mathbf{r}:

E(r)=14πε0qr2r^\mathbf{E}(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0} \frac{q}{|\mathbf{r}|^2} \hat{\mathbf{r}}

Superposition Principle: The field due to a collection of charges is the vector sum of individual fields.

2.2 Gauss's Law Applications

Example: Infinite plane of charge with surface charge density σ\sigma.

Choose a Gaussian "pillbox" of area AA straddling the plane. By symmetry, E\mathbf{E} is perpendicular to the plane. Gauss's law:

2EA=σAε0    E=σ2ε02EA = \frac{\sigma A}{\varepsilon_0} \implies E = \frac{\sigma}{2\varepsilon_0}

The field is uniform and perpendicular to the plane, pointing away from positive charge.

Example: Uniformly charged sphere of radius RR with total charge QQ.

For r>Rr > R: E=Q4πε0r2r^\mathbf{E} = \frac{Q}{4\pi\varepsilon_0 r^2} \hat{\mathbf{r}} (identical to a point charge).

For r<Rr < R: E=Qr4πε0R3E = \frac{Qr}{4\pi\varepsilon_0 R^3} (linear in rr).

2.3 Electric Potential

The electric potential is defined by E=V\mathbf{E} = -\nabla V (for electrostatics, where ×E=0\nabla \times \mathbf{E} = \mathbf{0}).

For a point charge: V(r)=14πε0qrV(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0} \frac{q}{r} (choosing V()=0V(\infty) = 0).

Theorem 2.1. ×E=0\nabla \times \mathbf{E} = \mathbf{0} in electrostatics implies E\mathbf{E} is conservative, so the line integral ABEdl=V(A)V(B)\int_A^B \mathbf{E} \cdot d\mathbf{l} = V(A) - V(B) is path-independent.

2.4 Poisson's and Laplace's Equations

Substituting E=V\mathbf{E} = -\nabla V into Gauss's law:

(V)=2V=ρε0\nabla \cdot (-\nabla V) = -\nabla^2 V = \frac{\rho}{\varepsilon_0}

This is Poisson's equation:

2V=ρε0\nabla^2 V = -\frac{\rho}{\varepsilon_0}

In regions with ρ=0\rho = 0, this reduces to Laplace's equation:

2V=0\nabla^2 V = 0

Theorem 2.2 (Uniqueness). The solution to Laplace's (or Poisson's) equation in a region is unique given either Dirichlet boundary conditions (VV specified on the boundary) or Neumann boundary conditions (V/n\partial V / \partial n specified on the boundary).

2.5 Worked Example

Problem. Two infinite conducting plates at x=0x = 0 and x=dx = d are held at potentials V=0V = 0 and V=V0V = V_0 respectively. Find the potential and field between them.

Solution. Between the plates, ρ=0\rho = 0, so 2V=0\nabla^2 V = 0. By symmetry, VV depends only on xx:

d2Vdx2=0    V(x)=Ax+B\frac{d^2V}{dx^2} = 0 \implies V(x) = Ax + B

Boundary conditions: V(0)=0    B=0V(0) = 0 \implies B = 0. V(d)=V0    A=V0/dV(d) = V_0 \implies A = V_0/d.

V(x)=V0dx,E=dVdxx^=V0dx^V(x) = \frac{V_0}{d} x, \quad \mathbf{E} = -\frac{dV}{dx}\hat{\mathbf{x}} = -\frac{V_0}{d}\hat{\mathbf{x}}

\blacksquare

3. Magnetostatics

3.1 The Biot-Savart Law

The magnetic field due to a steady current II in a wire element dld\mathbf{l}:

dB=μ0I4πdl×r^r2d\mathbf{B} = \frac{\mu_0 I}{4\pi} \frac{d\mathbf{l} \times \hat{\mathbf{r}}}{r^2}

For a complete circuit:

B(r)=μ0I4πdl×r^rr2\mathbf{B}(\mathbf{r}) = \frac{\mu_0 I}{4\pi} \oint \frac{d\mathbf{l} \times \hat{\mathbf{r}}'}{|\mathbf{r} - \mathbf{r}'|^2}

3.2 Ampere's Law

For steady currents (E/t=0\partial \mathbf{E} / \partial t = 0):

CBdl=μ0Ienc\oint_C \mathbf{B} \cdot d\mathbf{l} = \mu_0 I_{\mathrm{enc}}

Example: Infinite straight wire carrying current II.

By cylindrical symmetry, BB is constant on circles centred on the wire. Choose an Amperian loop of radius rr:

Bdl=B2πr=μ0I    B=μ0I2πr\oint \mathbf{B} \cdot d\mathbf{l} = B \cdot 2\pi r = \mu_0 I \implies B = \frac{\mu_0 I}{2\pi r}

Example: Solenoid. For a long solenoid with nn turns per unit length carrying current II:

B=μ0nI(inside),B=0(outside)B = \mu_0 n I \quad \mathrm{(inside)}, \quad B = 0 \quad \mathrm{(outside)}

3.3 Magnetic Vector Potential

Since B=0\nabla \cdot \mathbf{B} = 0, we can write B=×A\mathbf{B} = \nabla \times \mathbf{A}, where A\mathbf{A} is the magnetic vector potential.

In the Coulomb gauge (A=0\nabla \cdot \mathbf{A} = 0), the vector potential satisfies

2A=μ0J\nabla^2 \mathbf{A} = -\mu_0 \mathbf{J}

This is Poisson's equation for each component of A\mathbf{A}.

For a current loop, the solution is:

A(r)=μ04πJ(r)rrd3r\mathbf{A}(\mathbf{r}) = \frac{\mu_0}{4\pi} \int \frac{\mathbf{J}(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|}\, d^3\mathbf{r}'

4. Electrodynamics

4.1 Faraday's Law of Induction

A changing magnetic field induces an electric field:

×E=Bt\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}

Lenz's Law: The induced EMF opposes the change in flux that produced it.

Example. A circular loop of radius RR in a uniform magnetic field B(t)=B0cos(ωt)z^\mathbf{B}(t) = B_0 \cos(\omega t)\,\hat{\mathbf{z}}.

The flux: ΦB=πR2B0cos(ωt)\Phi_B = \pi R^2 B_0 \cos(\omega t).

The induced EMF: E=dΦBdt=πR2B0ωsin(ωt)\mathcal{E} = -\frac{d\Phi_B}{dt} = \pi R^2 B_0 \omega \sin(\omega t).

4.2 Displacement Current

Maxwell's key insight: Ampere's law ×B=μ0J\nabla \times \mathbf{B} = \mu_0 \mathbf{J} is inconsistent with the continuity equation. Adding the displacement current term μ0ε0E/t\mu_0 \varepsilon_0 \partial \mathbf{E}/\partial t resolves this:

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

4.3 Worked Example

Problem. A parallel-plate capacitor with circular plates of radius RR is being charged by a current II. Find the magnetic field between the plates at distance rr from the axis.

Solution. Between the plates, J=0\mathbf{J} = 0, but there is a changing electric field. The displacement current density is Jd=ε0EtJ_d = \varepsilon_0 \frac{\partial E}{\partial t}.

E=σε0=QπR2ε0E = \frac{\sigma}{\varepsilon_0} = \frac{Q}{\pi R^2 \varepsilon_0}, so Et=IπR2ε0\frac{\partial E}{\partial t} = \frac{I}{\pi R^2 \varepsilon_0}.

By symmetry, use an Amperian loop of radius r<Rr < R:

Bdl=μ0ε0tEdA\oint \mathbf{B} \cdot d\mathbf{l} = \mu_0 \varepsilon_0 \frac{\partial}{\partial t}\int \mathbf{E} \cdot d\mathbf{A}

B2πr=μ0ε0IπR2ε0πr2=μ0Ir2R2B \cdot 2\pi r = \mu_0 \varepsilon_0 \cdot \frac{I}{\pi R^2 \varepsilon_0} \cdot \pi r^2 = \frac{\mu_0 I r^2}{R^2}

B=μ0Ir2πR2B = \frac{\mu_0 I r}{2\pi R^2}

\blacksquare

5. Electromagnetic Waves

5.1 The Wave Equation

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.

5.2 Properties of EM Waves

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 (in general; 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).

5.3 Worked Example

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^2, i.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

6. Potentials and Gauge Transformations

6.1 Scalar and Vector Potentials

We can express the fields in terms of potentials:

E=VAt,B=×A\mathbf{E} = -\nabla V - \frac{\partial \mathbf{A}}{\partial t}, \quad \mathbf{B} = \nabla \times \mathbf{A}

In electrostatics, A=0\mathbf{A} = \mathbf{0} and E=V\mathbf{E} = -\nabla V.

6.2 Gauge Transformations

The potentials are not unique. The transformation

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

for any scalar function χ(r,t)\chi(\mathbf{r}, t) leaves E\mathbf{E} and B\mathbf{B} unchanged. This is a gauge transformation.

Common gauges:

  • Coulomb gauge: A=0\nabla \cdot \mathbf{A} = 0. Useful in magnetostatics.
  • Lorenz gauge: A+μ0ε0Vt=0\nabla \cdot \mathbf{A} + \mu_0 \varepsilon_0 \frac{\partial V}{\partial t} = 0. Simplifies the wave equations for VV and A\mathbf{A}:

2Vμ0ε02Vt2=ρε0\nabla^2 V - \mu_0 \varepsilon_0 \frac{\partial^2 V}{\partial t^2} = -\frac{\rho}{\varepsilon_0}

2Aμ0ε02At2=μ0J\nabla^2 \mathbf{A} - \mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{A}}{\partial t^2} = -\mu_0 \mathbf{J}

Common Pitfall

The Lorenz gauge (with one "r") is named after Ludvig Lorenz, not Hendrik Lorentz. It is frequently misspelled "Lorentz gauge." The two are different people, and the correct spelling is "Lorenz gauge."