The Physics of Electromagnetism
From medieval lodestones to Maxwell's wave equations and modern MRI systems — a complete technical survey of one of the four fundamental forces of nature.
Historical Evolution
Maxwell's Equations
Oliver Heaviside — a self-taught physicist and former telegrapher — reduced Maxwell's original 20 scalar equations into the four elegant partial differential equations recognized today, by introducing vector calculus and elevating the E and B fields to fundamental primacy.
Gauss's Law for Electricity
The divergence of the electric field is proportional to local charge density. Electric charges act as absolute sources or sinks of field lines.
Gauss's Law for Magnetism
The divergence of the magnetic field is always zero — strictly forbidding magnetic monopoles. Field lines form continuous unbroken loops.
Faraday's Law of Induction
A time-varying magnetic field induces a circulating, non-conservative electric field. The negative sign (Lenz's Law) ensures the induced field opposes the change.
Ampère–Maxwell Law
Circulating magnetic fields arise from steady currents J and Maxwell's crucial displacement current — enabling self-sustaining wave propagation.
Coupled with the Lorentz force law \(\mathbf{F} = q(\mathbf{E} + \mathbf{v} \times \mathbf{B})\), these four equations entirely dictate the macroscopic behavior of classical electromagnetism.
Biot–Savart & Ampère's Circuital Law
In static scenarios where fields do not vary over time, the generation of magnetic fields by steady direct currents is modeled by the Biot–Savart Law (1820):
Differential Field Element
Relates the differential magnetic field contribution from a current element \(Id\mathbf{l}\) at distance \(r\).
Circuital Law
The closed line integral of B around any loop equals \(\mu_0\) times the enclosed current. Most useful for symmetric geometries like solenoids.
Pivotal Experiments
Faraday's Induction Ring
Built in his basement at the Royal Institution, Faraday constructed the world's first electrical transformer using a soft iron ring wound with two separate copper coils insulated with cotton. The galvanometer deflected — but only transiently — revealing that EMF is induced not by the presence of a field, but by its change over time.
Hertz's Spark-Gap Transmitter
Under Helmholtz's encouragement, Hertz built a spark-gap oscillator using a Ruhmkorff induction coil. The transmitter's cylindrical brass dipole, 26 cm long, bridged a gap with a plasma arc and emitted invisible EM radiation. His resonant copper loop receiver detected waves up to 50 feet away, confirming Maxwell's prediction and proving waves travel at exactly \(c\).
EM Wave Visualizer
The electromagnetic wave equation \(\nabla^2 \mathbf{E} - \mu_0\epsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2} = 0\) predicts transverse waves propagating at \(c = \frac{1}{\sqrt{\mu_0\epsilon_0}} = 299{,}792{,}458\text{ m/s}\). Adjust frequency and amplitude below.
Relativistic Formulation
Classical electromagnetism is recognized as the first inherently relativistic field theory. Galilean transformations fail because light must travel at exactly \(c\) for all observers — Einstein resolved this with Lorentz transformations and Minkowski spacetime.
Space and time are inextricably bound in a four-displacement vector \(x^\alpha = (ct, x, y, z)\). The geometry is quantified by the Minkowski Metric Tensor \(\eta^{\mu\nu}\) with signature \((+\,-\,-\,-)\).
Four-current: J^ν = (cρ, J)
Four-potential: A^α = (φ/c, A)
The electric scalar potential and magnetic vector potential are unified into a single covariant electromagnetic four-potential, making gauge invariance manifest.
The E and B fields are not invariant standalone vectors. Under a Lorentz boost, a purely electric field in one frame appears as a mixture of E and B in a moving frame. They are unified into the antisymmetric rank-2 Faraday tensor:
The six independent off-diagonal entries precisely map to the three E and three B field components, with antisymmetry \(F_{\alpha\beta} = -F_{\beta\alpha}\) ensuring zero diagonal entries.
The entire structure of classical electrodynamics reduces to two tensor equations. The inhomogeneous Maxwell equations (Gauss + Ampère) become:
The homogeneous equations (Faraday + Gauss magnetism) become the Bianchi identity \(\partial_{[\alpha} F_{\beta\gamma]} = 0\). These two equations encode all four of Maxwell's laws in a fully Lorentz-covariant form.
The four-force on a charged particle is \(dp^\alpha/d\tau = qF^{\alpha\beta} u_\beta\), valid at any velocity up to \(c\). Energy density and the Poynting vector (energy flux) are:
S = (1/μ₀)(E × B)
When an EM wave strikes a surface, momentum density \(\mathbf{S}/c^2\) transfers as radiation pressure — crucial in astrophysics and laser mechanics.
MRI: Applied Electromagnetism
Magnetic Resonance Imaging represents the most sophisticated practical implementation of classical electrodynamics — manipulating static fields, RF waves, and spatially varying linear gradients to extract anatomical images from hydrogen protons in biological tissue.
Larmor Precession
Protons precess at \(f_0 = \gamma B_0\). For hydrogen, \(\gamma \approx 42.58\) MHz/T. At 3T, protons precess at ~127 MHz — a radio frequency.
RF Excitation
A short RF pulse at precisely the Larmor frequency tips the net magnetization into the transverse plane via resonant energy absorption.
Slice Selection
A gradient applied during the RF pulse restricts resonance to a single 2D plane by making precession frequency a linear function of position.
Phase Encoding
A transient gradient introduces proportional phase shifts along one spatial axis, encoding positional information into the phase of the signal.
k-Space & Fourier
Acquired data maps into a 2D spatial frequency domain (k-space), reconstructed into the final image via a 2D discrete Fourier transform.
Eddy Currents
Rapidly switching gradients induce Faraday eddy currents in conductive hardware. Active shielding — a secondary shield coil set — cancels external flux by an order of magnitude.
Anti-Helmholtz Spacing
Two coaxial loops spaced at exactly \(\sqrt{3}\times a\) cancel the third-order spherical harmonic \(z^3\) error, achieving <5% linearity deviation within a central spherical volume 0.5–0.6a from the isocenter.
Optimal Linearity Geometry
Golay (1958) proved that optimal transverse gradient linearity requires coil length 2.18r, 120° arc subtension, and inner margin separation of 0.78r.
Gradient Performance Specifications
| Metric | Definition | Typical Range (1.5T–3T) |
|---|---|---|
| Peak Gradient Strength | Maximum deliverable gradient amplitude (mT/m). Dictates ultra-high spatial resolution and diffusion weighting capability. | 45–100 mT/m (Specialty: up to 300 mT/m) |
| Rise Time | Time for gradient hardware to ramp from zero to peak amplitude. Shorter rise times enable faster pulse sequences. | 0.1–0.3 ms |
| Slew Rate | Peak Strength ÷ Rise Time (T/m/s). Governs minimum echo times and echo spacing in fast-spin sequences. | 100–200 T/m/s (Specialty: up to 750 T/m/s) |
| Effective Capability | When all three axes driven equally, vectorial combination boosts effective gradient by factor of √3 (~1.73×). | 1.73× single-axis spec |
| Duty Cycle | Fraction of time at maximum amplitude without exceeding thermal limits. Constrained by water-cooling efficiency. | Sequence-dependent |
Modern scanners have abandoned purely geometric "classic" Golay designs in favor of computational optimization algorithms that generate complex, asymmetrical "fingerprint" winding patterns — maximizing gradient linearity while simultaneously minimizing inductance, resistance, acoustic vibration, and thermal heating.