Electromagnetism & Maxwell’s Equations — A Guided Learning Path

Electromagnetism is the force responsible for virtually all chemistry, all light, and all electronics. This learning guide walks you from the basics of electric charge through to the full Maxwell equations, electromagnetic induction, and the prediction of light as an electromagnetic wave — building conceptual and mathematical understanding step by step, with interactive simulations at every stage.

Prerequisites and Learning Goals

This guide assumes you are comfortable with basic calculus (derivatives, integrals) and vectors (dot product, cross product, gradient, divergence, curl). A first course in mechanics is helpful but not required. By the end you should be able to:

Module 1: Electric Fields and Coulomb’s Law

Electric charge and the Coulomb force. All electromagnetic phenomena begin with electric charge. Charge is quantised (elementary charge e = 1.602×10−19 C) and conserved. Like charges repel; unlike charges attract. Coulomb’s law gives the force between two point charges q1, q2 separated by distance r:
F12 = k q1q2/r² &rcirc;12
k = 1/(4πε0) ≈ 8.99 × 10&sup9; N·m²/C²
ε0 = 8.854 × 10−12 F/m (permittivity of free space)
The force is along the line joining the charges (central force), falls as 1/r², and obeys the superposition principle: the total force on a charge equals the vector sum of all pairwise Coulomb forces.
The electric field concept. Instead of thinking about force between pairs, introduce the electric field E at a point P as the force that would be exerted on a unit positive test charge placed at P (in the absence of the test charge disturbing the source charges). For a point charge Q at the origin:
E(r) = kQ/r² &rcirc; = Q/(4πε0r²) &rcirc;

Force on charge q: F = qE
The electric field is a vector field defined everywhere in space. Field lines start on positive charges and end on negative charges (or go to infinity), and the field strength is proportional to line density.
Gauss’s law and symmetry. The total electric flux through any closed surface equals the enclosed charge divided by ε0:
&oiint; E·dA = Qenc0     (integral form)

∇·E = ρ/ε0     (differential form, ρ = charge density)
For a sphere of charge Q: choose a spherical Gaussian surface of radius r > R. By symmetry, E is uniform on the surface and radially outward, so: E × 4πr² = Q/ε0, giving E = Q/(4πε0r²) &rcirc;. Gauss’s law is most powerful for distributions with spherical, cylindrical, or planar symmetry.

Static Electricity Simulator: place positive and negative charges on a canvas and watch the electric field lines and equipotential surfaces update in real time. Try arranging a dipole (+q, −q) and verify that field lines start on + and end on −.

Module 2: Electric Potential and Energy

Electric potential. Because the electrostatic force is conservative (∇×E = 0 in static case), we can define a scalar potential V such that E = −∇V. The potential at a distance r from a point charge Q:
V(r) = kQ/r = Q/(4πε0r)

For a charge distribution: V = k ∫ dq/r'   (superposition)

Work to move charge q from A to B: W = q(VA − VB)
Equipotential surfaces are perpendicular to field lines. No work is done moving a charge along an equipotential.
Capacitors and energy storage. A capacitor stores charge Q on two conductors at potential difference V: Q = CV (capacitance C in farads). A parallel-plate capacitor with plate area A and separation d has C =&ε0A/d. Energy stored in the electric field:
UE = ½CV² = Q²/(2C) = ½ε0E² × (volume)

Energy density: uE = ½ε0E²   [J/m³]
The energy density in the field is a general result — it applies everywhere, not just inside capacitors.

Module 3: Magnetic Fields and Biot-Savart

The magnetic force on moving charges. A charged particle moving with velocity v in a magnetic field B experiences the Lorentz force:
F = qv × B

For a current element Idl: dF = Idl × B
Crucially, the magnetic force is always perpendicular to the velocity, so it does no work. It alters the direction of motion but not the kinetic energy. A charged particle in a uniform field moves in a circle at the cyclotron frequency ωc = qB/m (Larmor radius rL = mv/qB).
Sources of magnetic fields: Biot-Savart law. Moving charges create magnetic fields. The contribution to B from a current element Idl at position r from the element:
dB = (μ0/4π) Idl × &rcirc;/r²

For a long straight wire at distance R:
B = μ0I / (2πR)     (circling the wire)
Magnetic field lines form closed loops — they have no beginning or end.
Ampère’s law. For a closed Ampèrian loop C enclosing current Ienc:
B·dl = μ0Ienc     (integral form)

∇×B = μ0J     (differential form, J = current density)

Solenoid (n turns/m): B = μ0nI   (inside, uniform)
Just as Gauss’s law exploits symmetry for E, Ampère’s law is most powerful for symmetric current distributions: infinite wire, solenoid, toroid.
Gauss’s law for magnetism. Unlike electric field lines which start/end on charges, magnetic field lines always close on themselves:
&oiint; B·dA = 0     (no magnetic monopoles)

∇·B = 0     (differential form)
Magnetic flux through any closed surface is zero. The absence of magnetic monopoles is one of the fundamental asymmetries between electricity and magnetism.

Module 4: Faraday Induction and Lenz’s Law

Magnetic flux and Faraday’s law. Define the magnetic flux through a surface S:
ΦB = ∫∫S B·dA

Faraday’s law: ϵ = −dΦB/dt

N-turn coil: ϵ = −N dΦB/dt
The induced EMF drives a current in the coil. Note: you do not need to have a physical circuit — a changing B creates a circulating E field in empty space (from the integral form: ∮E·dl = −dΦB/dt).
Lenz’s law — the sign matters. The minus sign in Faraday’s law enforces Lenz’s law: the induced current always flows in the direction that opposes the change in flux. If flux is increasing through a loop, the induced current creates a field opposing the increase. This is a direct consequence of energy conservation: if the induced current aided the flux change, you could create energy from nothing.

Physical examples: a magnet dropped through a copper tube falls slowly (eddy currents brake it); regenerative braking in electric cars (the motor acts as a generator, opposing wheel rotation); transformer heating in iron cores.

Motional EMF and generators. When a conductor moves in a magnetic field, the free electrons experience F = qv×B, creating a charge separation and an EMF. For a rod of length L moving at velocity v perpendicular to field B:
ϵ = BLv     (motional EMF)

Rotating coil (N turns, area A, angular frequency ω):
ϵ(t) = NBAω sin(ωt)     (AC generator)
🧲

Electromagnetic Induction Simulator: choose the bar-magnet mode and drag the magnet toward and away from the coil. Notice the galvanometer deflects only while the magnet is moving, and reverses direction when the motion reverses. Switch to AC Generator mode to see ϵ(t) = NBAω sin(ωt) traced in real time as the coil spins.

Module 5: Maxwell’s Displacement Current

The missing term in Ampère’s law. Ampère’s law (∇×B = μ0J) is inconsistent when applied to a charging capacitor. Maxwell introduced the displacement current density:
JD = ε0E/∂t

Corrected Ampère-Maxwell: ∇×B = μ0J + μ0ε0E/∂t
Physical meaning: a time-varying electric field acts as a source of magnetic field, even with no real current present. This is the key insight that leads to self-sustaining electromagnetic waves.

Module 6: The Full Maxwell Equations

The four Maxwell equations in differential form in SI units:

(I) ∇·E = ρ/ε0      [Gauss electric]

(II) ∇·B = 0      [Gauss magnetic]

(III) ∇×E = −∂B/∂t      [Faraday]

(IV) ∇×B = μ0J + μ0ε0E/∂t      [Ampère-Maxwell]

Together with the Lorentz force law (F = q(E + v×B)), these equations completely describe all classical electromagnetic phenomena for given charge and current distributions.

In a linear, isotropic, homogeneous material with relative permittivity εr and relative permeability μr, replace ε0 by ε = εrε0 and μ0 by μ =μrμ0. The wave speed in the medium becomes v = c/n where the refractive index n = √(εrμr).

Module 7: Electromagnetic Waves

Deriving the wave equation. In vacuum, take the curl of Faraday’s law (III) and substitute (IV):
∇×(∇×E) = −∇×(∂B/∂t)

Using identity ∇×(∇×E) = ∇(∇·E) − ∇²E and ∇·E = 0 (in vacuum):

∇²E = μ0ε0 ∂²E/∂t²

Wave speed: c = 1/√(μ0ε0) = 2.9979 × 10&sup8; m/s
Plane wave solution. The simplest solution is a monochromatic plane wave (angular frequency ω, wave vector k = ω/c) propagating in the +x direction:
Ey(x,t) = E0 cos(kx − ωt + φ)
Bz(x,t) = (E0/c) cos(kx − ωt + φ)

Key relations: EB ⊥ propagation direction
Amplitude ratio: E0/B0 = c
B = (1/c) kˆ × E
The electric and magnetic energy densities are equal (½ε0E² = B²/(2μ0)).
Energy transport: the Poynting vector. The electromagnetic energy flux (power per unit area) is:
S = (1/μ0) E × B     [W/m²]

Time-averaged Poynting flux: ⟨S⟩ = E0B0/(2μ0) = E0²/(2μ0c)

Radiation pressure: Prad = ⟨S⟩/c   (absorbed)   or   2⟨S⟩/c (reflected)
Radiation pressure is tiny for everyday light intensities but significant in stellar interiors (supporting the Sun against gravitational collapse), near-critical laser pulses, and proposed solar sails.
📡

Maxwell Waves Simulator renders the oscillating E and B field vectors as an animated transverse wave. Verify that they are perpendicular and in phase, and that the Poynting vector (cross product) points along the propagation direction.

Module 8: Electromagnetic Induction Applications

Self-inductance and RL circuits. A coil opposes changes in current through itself via self-EMF ϵ = −L dI/dt. In a series RL circuit (resistance R, inductance L) driven by EMF ϵ0:
I(t) = (ϵ0/R)(1 − e−t/τ)

Time constant: τ = L/R     (inductive time constant)

Stored magnetic energy: UL = ½LI²
Solenoid: L = μ0N²A/ℓ
LC oscillators and RLC circuits. In an ideal LC circuit (no resistance), energy oscillates between the inductor (magnetic) and capacitor (electric) at the resonant frequency:
ω0 = 1/√(LC)     (resonance)

Q(t) = Q0 cos(ω0t + φ)
I(t) = −Q0ω0 sin(ω0t + φ)

With resistance R: ωd = √(ω0² − (R/2L)²)
The LC resonator is the electromagnetic analogue of the mechanical harmonic oscillator (spring-mass). It is fundamental to radio receivers (tuning circuits), oscillators, and filters.
Transformers and impedance matching. An ideal transformer with primary turns Np and secondary turns Ns:
Vs/Vp = Ns/Np
Is/Ip = Np/Ns
Zin = (Np/Ns)² Zload     (impedance transformation)
The impedance transformation property is used in audio amplifiers to match the high output impedance of an amplifier to a low-impedance speaker (8 Ω), and in RF circuits to match antenna impedance (often 50 Ω) to transmission line impedance.

Module 9: Boundary Conditions and Waveguides

At an interface between two media, Maxwell’s equations imply specific boundary conditions on the field components. For the tangential E and normal B (from Gauss and Faraday):

E1t = E2t     (tangential E continuous)
B1n = B2n     (normal B continuous)
D1n − D2n = σfree     (normal D discontinuous by surface charge)
H1t − H2t = Kfree     (tangential H discontinuous by surface current)

At a perfect conductor surface: Etangential = 0 and Bnormal = 0. These conditions confine electromagnetic waves inside metallic pipes — waveguides. A rectangular waveguide of width a and height b supports TEmn and TMmn modes with cutoff frequencies fc = c/(2) √[(m/a)² + (n/b)²]. Below the cutoff frequency, waves are evanescent (exponentially decaying, no net energy transport). Waveguides are used in microwave engineering: radar, satellite communications, lab microwave measurements.

Module 10: Relativity and Electromagnetism

Maxwell’s equations are Lorentz-covariant: they have the same form in all inertial reference frames. This is why Einstein’s 1905 paper was titled “On the Electrodynamics of Moving Bodies” — special relativity was born from the inconsistency of classical mechanics with Maxwell’s equations.

Key consequences for electromagnetism:

Summary: the hierarchy of classical EM theory

Electrostatics (static charges: Gauss E, Coulomb) → Magnetostatics (steady currents: Gauss B, Biot-Savart, Ampère) → Faraday induction (changing B → E) → Maxwell displacement current (changing E → B) → Self-sustaining EM waves (light, radio, X-rays) → Lorentz covariance (4D unification, tensor Fμν) → QED (quantisation of the EM field, photons, virtual particles)

Recommended Simulations for This Module

Further Reading