Maxwell's Equations and Gauge Theory: Electromagnetism as a Principal Bundle

A compass needle aligns with a magnetic field threading through empty air. A lightning bolt drives a current through the atmosphere via the same electric force that pulls charged particles in a wire. Radio waves cross a continent at the speed of light, which turns out to be the same constant that emerges from the equations governing both electricity and magnetism. All of this follows from four equations that James Clerk Maxwell assembled in 1865.

The classical formulation in vector calculus is powerful and concrete, but it hides something. The equations have a symmetry—gauge invariance—that suggests the fields E and B are not the most fundamental objects. Behind them sits the vector potential A, and behind A sits the geometry of a principal fibre bundle over spacetime. The same geometric structure that appears in electromagnetism reappears, with a non-abelian gauge group, in the quantum chromodynamics of quarks and the electroweak theory of the Standard Model. In 2000, the Clay Mathematics Institute designated the question of whether the non-abelian (Yang–Mills) version of this structure always has a positive mass gap as one of the seven Millennium Prize Problems, with a million-dollar prize still unclaimed.

This post develops four parallel languages for Maxwell's equations: the classical vector calculus form, the coordinate-free language of differential forms on Minkowski spacetime, the spacetime algebra (geometric algebra) encoding all four equations in a single line, and the principal fibre bundle construction that reveals gauge theory as differential geometry. At each level the physical content is identical; what changes is how much of the underlying structure becomes visible.

1. The Classical Laws

The four Maxwell equations rest on four independent experimental observations, each discovered in a 46-year window between 1785 and 1831. Coulomb measured the force between static charges. Biot and Savart quantified the force between current-carrying wires. Faraday discovered that a moving magnet induces a voltage in a nearby coil. Ampere related closed line integrals of $\mathbf{B}$ to enclosed currents. In 1865, Maxwell added one term to Ampere's law to repair an inconsistency with charge conservation. As a consequence, he found that light is an electromagnetic wave.

Before stating the laws, two pieces of notation need grounding.

Charge is an intrinsic scalar property of matter, analogous to mass but with one fundamental difference: it can be positive, negative, or zero. Two particles with the same sign of charge repel each other; two with opposite signs attract. Coulomb quantified this force in 1785, and that measurement is the foundational observation that justifies introducing charge as a concept at all. Classical theory offers no deeper explanation for why matter carries charge; it is a label that organizes the forces we observe. Field theory adds structural content to that label: charge is the source term that couples matter to the electromagnetic field. Without charge, Gauss's law reduces to the homogeneous equation $\nabla \cdot \mathbf{E} = 0$, and the field has no reason to be nonzero anywhere. Charge is what forces the field to respond to the presence of matter.

Charge density $\rho(\mathbf{x}, t)$ extends this to continuous distributions. When charge is spread through a volume (as in a plasma or the conduction band of a metal), we cannot assign a charge to each geometric point; instead we ask how much charge occupies an infinitesimally small region near $\mathbf{x}$. Let $B_\varepsilon(\mathbf{x})$ be a ball of radius $\varepsilon$ centered at $\mathbf{x}$. Then

where $Q(B_\varepsilon, t)$ is the total charge inside that ball at time $t$ and $\mathrm{Vol}(B_\varepsilon) = \tfrac{4}{3}\pi\varepsilon^3$. This is the same construction as mass density in continuum mechanics: charge per unit volume, measured in coulombs per cubic meter. The total charge in any region $V$ is recovered by integration:

Point charges and the Dirac measure. An electron carries all of its charge at a single spatial point, to every precision ever measured. A density function for this situation must be zero everywhere except at the electron's location, yet integrate over all of space to the total charge $e$. No function in the ordinary (Lebesgue) sense can achieve this: a function that is zero almost everywhere has integral zero, regardless of the value at any isolated point.

The $k$-dimensional Dirac delta $\delta^{(k)}(\mathbf{x})$ resolves this by working in the broader class of distributions: continuous linear functionals on smooth test functions, rather than pointwise-valued maps. The three-dimensional Dirac delta centered at the origin is defined by its action on any smooth compactly supported test function $\varphi$:

It is the evaluation functional at the origin. A concrete way to approach it: $\delta^{(3)}(\mathbf{x})$ is the distributional limit of the family of normalized Gaussians

as $\sigma \to 0$. For any fixed smooth $\varphi$, one has $\int_{\mathbb{R}^3} g_\sigma\,\varphi \to \varphi(\mathbf{0})$ as $\sigma \to 0$. The Dirac delta is what that Gaussian family converges to in the distributional sense: a spike of unit mass that has collapsed to a single point. The $k$-dimensional version $\delta^{(k)}(\mathbf{x})$ is the same construction in $\mathbb{R}^k$, satisfying $\int_{\mathbb{R}^k} \delta^{(k)}(\mathbf{x})\,\varphi(\mathbf{x})\,\mathrm{d}^k x = \varphi(\mathbf{0})$.

The charge density of a point charge $Q$ at the origin is then $\rho = Q\,\delta^{(3)}(\mathbf{x})$, a distributional statement with full physical weight: substituting into Gauss's law and integrating over any sphere enclosing the origin recovers Coulomb's inverse-square field exactly. A point particle moving along a worldline $\boldsymbol{\gamma}(t)$ generalizes this to $\rho = Q\,\delta^{(3)}(\mathbf{x} - \boldsymbol{\gamma}(t))$, with current density $\mathbf{J} = Q\,\dot{\boldsymbol{\gamma}}(t)\,\delta^{(3)}(\mathbf{x} - \boldsymbol{\gamma}(t))$. In field theory, the interaction between a charged particle and the electromagnetic potential is the spacetime integral $\int J^\mu A_\mu\,\mathrm{d}^4x$. For a point particle, the delta function collapses this volume integral to a line integral $\oint_{\boldsymbol{\gamma}} A$ along the worldline. That line integral is precisely the geometric object whose ambiguity encodes gauge freedom and whose measurability produces the Aharonov–Bohm phase in Section 8.

In a conductor or plasma, $\rho$ is instead a smooth function spread continuously through space. Both cases, the point-supported delta and the smooth continuous distribution, are valid instances of the same distributional framework; only the degree of spatial concentration differs.

On quantum mechanics. The treatment above assigns charge density via Dirac deltas for point particles or smooth functions for continuous distributions, and in both cases the electron is treated as if it occupies a definite location or a definite extended region. Quantum mechanics alters this. An electron is not at a definite point; it occupies a quantum state described by a wavefunction $\psi(\mathbf{x},t)$, and $|\psi(\mathbf{x},t)|^2$ is the probability density for finding the electron near $\mathbf{x}$. The charge density it contributes is accordingly $\rho = -e|\psi|^2$ (negative because the electron carries charge $-e$). This is not charge spread spatially in the classical sense; it is the expectation value of the charge density operator in the state $\psi$. In an atom, $\psi$ is an orbital, and $-e|\psi|^2$ is a smooth function spread over the orbital volume that sources the electric field through Gauss's law exactly as a classical continuous distribution does. Section 10 develops this connection: how quantum charge densities enter Maxwell's equations, how the resulting electrostatic potential feeds back into the Schrödinger equation, and how the coupled system is solved computationally.

The closed surface integral $\oiint_S \mathbf{F} \cdot \mathrm{d}\mathbf{S}$ adds up the outward flux of a vector field $\mathbf{F}$ through every infinitesimal patch $\mathrm{d}\mathbf{S} = \hat{\mathbf{n}}\,\mathrm{d}A$ of a closed surface $S$ (one with no boundary, like a sphere or a box). It measures how much of $\mathbf{F}$ is leaving versus entering the enclosed region. The divergence theorem relates this to a volume integral: $\oiint_S \mathbf{F}\cdot\mathrm{d}\mathbf{S} = \iiint_V \nabla\cdot\mathbf{F}\,\mathrm{d}V$, which is what converts the integral form of each law into the differential (pointwise) form below.

Each law is stated as a proposition and derived from the corresponding experiment below.

Derivation. Coulomb's law gives the field of a point charge $Q$ at the origin as $\mathbf{E} = \frac{Q}{4\pi\varepsilon_0}\frac{\hat{\mathbf{r}}}{r^2}$. The flux through a sphere $S_r$ of radius $r$ is

By the divergence theorem, $\oiint_{S_r} \mathbf{E}\cdot\mathrm{d}\mathbf{S} = \iiint_V \nabla\cdot\mathbf{E}\,\mathrm{d}V$. Replacing $Q$ by a continuous distribution and equating the integrands gives $\nabla\cdot\mathbf{E} = \rho/\varepsilon_0$ pointwise. $\square$

Derivation. The Biot–Savart law gives $\mathbf{B}(\mathbf{x}) = \frac{\mu_0}{4\pi}\int \mathbf{J}(\mathbf{x}')\times\nabla_x\!\left(\frac{-1}{|\mathbf{x}-\mathbf{x}'|}\right)\mathrm{d}^3x'$, which can be written as [eq:b-potential] with $\mathbf{A} = \frac{\mu_0}{4\pi}\int\frac{\mathbf{J}(\mathbf{x}')}{|\mathbf{x}-\mathbf{x}'|}\,\mathrm{d}^3x'$. Then $\nabla\cdot\mathbf{B} = \nabla\cdot(\nabla\times\mathbf{A}) = 0$ by the identity $\nabla \cdot (\nabla \times \mathbf{A}) = 0$, which holds identically for any smooth vector field. $\square$

Derivation. Faraday observed in 1831 that a changing magnetic flux $\Phi_B = \iint_\Sigma \mathbf{B}\cdot\mathrm{d}\mathbf{S}$ through a loop drives an EMF:

For a stationary loop, $-\frac{\mathrm{d}\Phi_B}{\mathrm{d}t} = -\iint_\Sigma \frac{\partial\mathbf{B}}{\partial t}\cdot\mathrm{d}\mathbf{S}$. By Stokes' theorem, $\oint_C \mathbf{E}\cdot\mathrm{d}\boldsymbol{\ell} = \iint_\Sigma(\nabla\times\mathbf{E})\cdot\mathrm{d}\mathbf{S}$. Since this holds for every $\Sigma$, the integrands must match: $\nabla\times\mathbf{E} = -\partial_t\mathbf{B}$. $\square$

Derivation. Ampere's part. The Biot–Savart law for a long straight wire gives $\mathbf{B} = \frac{\mu_0 I}{2\pi s}\hat{\boldsymbol{\phi}}$ at distance $s$. Integrating around a circular Amperian loop: $\oint \mathbf{B}\cdot\mathrm{d}\boldsymbol{\ell} = \mu_0 I_{\mathrm{enc}}$. By Stokes, this becomes $\nabla\times\mathbf{B} = \mu_0\mathbf{J}$ in the static case.

Maxwell's correction. Taking the divergence: $0 = \nabla\cdot(\nabla\times\mathbf{B}) = \mu_0\nabla\cdot\mathbf{J}$, but charge conservation requires $\nabla\cdot\mathbf{J} = -\partial_t\rho$, which is nonzero during transients. By Gauss's law, $\partial_t\rho = \varepsilon_0\partial_t(\nabla\cdot\mathbf{E}) = \nabla\cdot(\varepsilon_0\partial_t\mathbf{E})$. Adding the displacement current $\varepsilon_0\partial_t\mathbf{E}$ to the right side restores consistency: $\nabla\cdot(\nabla\times\mathbf{B}) = \mu_0(\nabla\cdot\mathbf{J} + \partial_t\rho) = 0$ by continuity. $\square$

Proof. Apply $\nabla\times$ to [eq:faraday-eq] (Faraday's law): $\nabla\times(\nabla\times\mathbf{E}) = -\partial_t(\nabla\times\mathbf{B}) = -\mu_0\varepsilon_0\partial_t^2\mathbf{E}$ using [eq:ampere-eq] (vacuum Ampere–Maxwell law). The identity $\nabla\times(\nabla\times\mathbf{E}) = \nabla(\nabla\cdot\mathbf{E}) - \nabla^2\mathbf{E} = -\nabla^2\mathbf{E}$ (with $\nabla\cdot\mathbf{E} = 0$ in vacuum) gives $\nabla^2\mathbf{E} = \mu_0\varepsilon_0\partial_t^2\mathbf{E}$. Maxwell identified $c = 1/\sqrt{\mu_0\varepsilon_0} \approx 3\times10^8$ m/s with the known speed of light and concluded that light is an electromagnetic wave. The same derivation for $\mathbf{B}$ gives an identical equation. $\square$

2. The Vector Potential and Gauge Freedom

The equation $\nabla \cdot \mathbf{B} = 0$ says that the magnetic field has no divergence. A standard result from vector calculus (the Poincaré lemma in disguise) says that any divergence-free vector field on a simply connected domain can be written as the curl of another vector field. This defines the magnetic vector potential.

The fields E and B are determined by the pair $(\varphi, \mathbf{A})$, but not uniquely: different potentials can produce the same fields.

Proof. Direct substitution: $\mathbf{B}' = \nabla \times (\mathbf{A} + \nabla\chi) = \nabla \times \mathbf{A} + \nabla \times \nabla\chi = \mathbf{B} + 0 = \mathbf{B}$ since any gradient is curl-free. For E: $\mathbf{E}' = -\nabla(\varphi - \partial_t\chi) - \partial_t(\mathbf{A} + \nabla\chi) = -\nabla\varphi + \nabla\partial_t\chi - \partial_t\mathbf{A} - \partial_t\nabla\chi = \mathbf{E}$ since partial derivatives commute. $\square$

The freedom in choosing $\chi$ is called gauge freedom. Two gauge choices that often appear in practice:

The potentials $(\varphi, \mathbf{A})$ were originally treated as mathematical auxiliaries with no independent physical meaning. Section 8 shows this view is wrong: the Aharonov–Bohm effect demonstrates that the vector potential has measurable physical consequences in regions where $\mathbf{B} = 0$.

3. Differential Forms and the Faraday Tensor

Vector calculus works beautifully in flat three-dimensional space, but it cannot be easily extended to curved spacetime or to general dimensions. Differential forms provide the coordinate-free language that makes this extension natural. The key payoff: all four Maxwell equations reduce to two compact statements.

A $k$-form on a smooth manifold $M$ is a field of antisymmetric multilinear maps on the tangent space, eating $k$ tangent vectors and returning a number. These spaces are organized into the de Rham complex: a sequence of vector spaces connected by the exterior derivative $\mathrm{d}$, which maps $k$-forms to $(k+1)$-forms, subject to the central identity $\mathrm{d}^2 = 0$. This single algebraic fact encodes conservation laws, Bianchi identities, and the existence of gauge potentials all at once.

The de Rham complex is the continuum structure that discrete exterior calculus (DEC) faithfully discretizes: smooth manifolds become simplicial complexes, smooth $k$-forms become $k$-cochains, and the exterior derivative $\mathrm{d}$ becomes a combinatorial coboundary operator that inherits $\mathrm{d}^2 = 0$ exactly from the boundary identity $\partial^2 = 0$. Section 10 develops this construction and uses it to formulate Maxwell's equations on a discrete mesh, where $\nabla \cdot \mathbf{B} = 0$ and gauge invariance hold exactly at the discrete level rather than approximately. For Maxwell's equations, the relevant complex runs over four-dimensional Minkowski spacetime $\mathbb{R}^{1,3}$ up to $\Omega^4$.

Proof. We work in coordinates $(t,x,y,z)$ with $\eta = \mathrm{diag}(+1,-1,-1,-1)$. We use the convention in which the component form of $F$ is $F_{0i} = -E_i$, $F_{ij} = \varepsilon_{ijk}B_k$ (consistent with the STA bivector in Definition 4.2). In this convention the field 2-form reads

Part 1: $\mathrm{d}F = 0$. Because $\mathrm{d}F$ is a 3-form in four dimensions it has exactly four independent components. The general formula for the antisymmetrised derivative of a 2-form is

The $(x,y,z)$ component — the only purely spatial triple. Only the magnetic 2-forms contribute (the electric terms involve $\mathrm{d}t$ and vanish on a purely spatial triple):

Setting this to zero gives $\nabla\cdot\mathbf{B} = 0$, the no-monopole condition.

The $(t,x,y)$ component. Only $E_x$, $E_y$, and $B_z$ contribute, since these are the only field components whose basis 2-forms share two of the three indices $\{t,x,y\}$:

Setting this to zero gives $\partial_t B_z = \partial_y E_x - \partial_x E_y = -(\nabla\times\mathbf{E})_z$, the $z$-component of Faraday's law $\partial_t\mathbf{B} + \nabla\times\mathbf{E} = 0$. The remaining two components follow by cycling $(x,y,z)$:

which are the $x$- and $y$-components of $\partial_t\mathbf{B} + \nabla\times\mathbf{E} = 0$.

Part 2: $\mathrm{d}{\star}F = {\star}J$. The Hodge star in $(\mathbb{R}^{1,3}, \eta)$ acts on basis 2-forms by ${\star}(\mathrm{d}x^\mu\wedge \mathrm{d}x^\nu) = \tfrac{1}{2}\eta^{\mu\alpha}\eta^{\nu\beta}\varepsilon_{\alpha\beta\rho\sigma}\,\mathrm{d}x^\rho\wedge\mathrm{d}x^\sigma$ with $\varepsilon_{0123}=+1$. Applying this to $F$, using $\star\star F = -F$ to verify signs:

The Hodge star has swapped the electric and magnetic sectors. Now applying $\mathrm{d}$ to ${\star}F$:

The $(x,y,z)$ component. Only the spatial 2-forms in ${\star}F$ contribute a purely spatial 3-form:

The source 3-form ${\star}J$ (the Hodge dual of the current 1-form $J = \rho\,\mathrm{d}t - J_x\,\mathrm{d}x - J_y\,\mathrm{d}y - J_z\,\mathrm{d}z$) has $({\star}J)_{xyz} = -\rho$ in this convention. Setting $(\mathrm{d}{\star}F)_{xyz} = ({\star}J)_{xyz}$ gives $\nabla\cdot\mathbf{E} = \rho$, Gauss's law.

The $(t,x,y)$ component. The magnetic 2-forms $-B_x\,\mathrm{d}t\wedge\mathrm{d}x$, $-B_y\,\mathrm{d}t\wedge\mathrm{d}y$ and the spatial form $-E_z\,\mathrm{d}x\wedge\mathrm{d}y$ contribute:

Setting equal to $({\star}J)_{txy} = -J_z$ gives $(\nabla\times\mathbf{B})_z - \partial_t E_z = J_z$, the $z$-component of Ampere's law $\nabla\times\mathbf{B} - \partial_t\mathbf{E} = \mathbf{J}$. The $(t,y,z)$ and $(t,x,z)$ components give the remaining components by the same calculation.

All eight equations — four from [eq:maxwell-forms-eq] — are precisely Maxwell's four equations in classical form. $\square$

The forms language also makes the charge conservation law automatic: applying $\mathrm{d}$ to $\mathrm{d}{\star}F = {\star}J$ and using $\mathrm{d}^2 = 0$ gives $\mathrm{d}{\star}J = 0$, which in coordinates is $\partial_t \rho + \nabla \cdot \mathbf{J} = 0$, the continuity equation.

4. Spacetime Algebra: $\nabla F = J$

The differential forms approach breaks Maxwell's equations into two statements. Geometric algebra (Clifford algebra applied to spacetime) collapses them into one. The gain is more than notational: the algebraic structure makes the spinor formulation of relativistic quantum mechanics and the Yang–Mills generalization (Section 7) more transparent.

Proof. Write $\nabla = \gamma^\mu\partial_\mu$ with $\gamma^\mu = \eta^{\mu\nu}\gamma_\nu$, so $\gamma^0 = \gamma_0$ and $\gamma^i = -\gamma_i$. The Clifford product of a grade-1 vector $v$ with a grade-2 bivector $F$ decomposes by grade:

where $\langle\cdot\rangle_1$ (the dot part) decreases grade by 1 and $\langle\cdot\rangle_3$ (the wedge part) increases grade by 1. Since $F$ is a bivector, $\nabla F$ is the sum of a grade-1 vector and a grade-3 trivector.

Grade-1 part: inhomogeneous equations. Expanding $\gamma^\mu \cdot (\partial_\mu F)$ in the canonical blade basis (computed explicitly via the Clifford multiplication rules $\gamma_\mu\gamma_\nu = -\gamma_\nu\gamma_\mu$ for $\mu\neq\nu$ and $\gamma_\mu^2 = \eta_{\mu\mu}$):

The $\gamma_0$ coefficient is $\nabla\cdot\mathbf{E}$; equating to $J^0 = \rho$ gives Gauss's law. The $\gamma_k$ coefficient is $(\nabla\times\mathbf{B})_k - \partial_t E_k$; equating to $J^k$ gives Ampere's law $\nabla\times\mathbf{B} - \partial_t\mathbf{E} = \mathbf{J}$.

Grade-3 part: homogeneous equations. Expanding $\gamma^\mu\wedge(\partial_\mu F)$ in the trivector basis:

The $\gamma_0\gamma_k\gamma_l$ coefficients are the components of $\partial_t\mathbf{B} + \nabla\times\mathbf{E}$; each set to zero gives a component of Faraday's law. The purely spatial trivector coefficient is

which set to zero gives $\nabla\cdot\mathbf{B}=0$. Since $J = J^\mu\gamma_\mu$ is a grade-1 object, equating $\nabla F = J$ forces the grade-3 part to vanish identically, recovering all four homogeneous equations simultaneously. The grade-1 equation then yields the four inhomogeneous equations of [eq:maxwell-forms-eq]. $\square$

5. Principal Fibre Bundles

The gauge potential A of electromagnetism is a 1-form on spacetime with values in a Lie algebra. This structure has a precise geometric home: it is a connection on a principal fibre bundle. Understanding that home requires building up the definition carefully, but the physical picture is always a family of gauge-equivalent copies of the structure group, one sitting over each point of spacetime.

Think of a helix wrapping around a cylinder. The cylinder is the base space (the spacetime you navigate), the circle at each height is the fibre (the gauge group), and the helix tells you how to move consistently from one fibre to the next. That consistent prescription is a connection.

6. Connections and Curvature

A connection on a principal bundle answers the question: given a path in the base manifold, how do we lift it to a path in the total space in a consistent, $G$-equivariant way? The answer defines parallel transport, and the failure of parallel transport to close around loops is curvature.

Concretely, a connection splits the tangent space of $P$ at each point into a vertical part (tangent to the fibre, "staying at the same base point") and a horizontal part (transverse to the fibre, "moving in the base direction"). The choice of horizontal subspace, required to be G-equivariant and smooth, is the connection.

7. Yang–Mills Theory

In 1954, Chen-Ning Yang and Robert Mills asked: what happens if the gauge group is non-abelian? The same geometric structure (principal bundle, connection, curvature) applies, but the non-commutativity of the group means the gauge bosons themselves carry charge and interact. In electromagnetism, photons are electrically neutral and do not interact with each other. In a non-abelian gauge theory, the gauge bosons have charges under their own gauge group and self-interact.

8. The Aharonov–Bohm Effect

Consider a solenoid: a tightly wound coil of wire carrying current. Inside the solenoid there is a strong magnetic field. Outside, by symmetry and Ampere's law, the field is nearly zero. Classical electromagnetism predicts that electrons traveling only in the field-free region outside the solenoid should be unaffected by whatever happens inside.

In 1959, Yakir Aharonov and David Bohm predicted that this classical intuition is wrong. Even though $B = 0$ outside, the vector potential $A$ is nonzero (it must be, since $B = \nabla\times A$ is nonzero inside). And the vector potential has a directly measurable effect on quantum mechanics: it shifts the phase of the electron wavefunction along any path. If two electron paths encircle the solenoid from opposite sides and recombine, the phase difference between the paths depends on the total enclosed magnetic flux.

The prediction was confirmed experimentally by Chambers (1960) and definitively by Tonomura and collaborators (1986) using electron holography with a toroidal magnet fully enclosed in a superconducting shield, ensuring $\mathbf{B} = 0$ everywhere on the electron paths to very high precision.

In the language of fibre bundles, the AB phase is the holonomy of the connection A around the closed loop formed by the two paths. Holonomy measures how much parallel transport around a loop fails to return to the starting point in the fibre. For an abelian bundle, holonomy is a phase in $U(1)$.

9. Dirac Monopoles and Chern Classes

No magnetic monopole has ever been observed. But Dirac showed in 1931 that if a single monopole exists anywhere in the universe, electric charge must be quantized in discrete units: every particle's charge is an integer multiple of the electron charge. This is a striking example of topology constraining physics.

The argument uses the same two-patch construction that classifies line bundles over spheres, now called the Wu–Yang construction.