The Unified Foundation of Classical Electrodynamics
An electromagnetic wave (EM wave) is a disturbance that propagates through space, carrying energy. Crucially, it is not a mechanical wave requiring a medium; it is a self-propagating vibration of electric ($\mathbf{E}$) and magnetic ($\mathbf{B}$) fields. The fundamental insight is that a changing magnetic field creates an electric field (Faraday's Law), and a changing electric field creates a magnetic field (Ampère-Maxwell Law, including displacement current). This interdependence allows the two fields to continuously generate each other, leading to sustained propagation across vast distances.
The EM wave represents the ultimate synthesis of classical physics, unifying the previously separate phenomena of electricity, magnetism, and light into a single, cohesive framework.
In a vacuum, where there are no free charges ($\rho = 0$) and no free currents ($\mathbf{J} = \mathbf{0}$), the four fundamental equations simplify:
The full set of equations includes charge density ($\rho$) and current density ($\mathbf{J}$), which act as sources for the fields.
The wave equation is derived directly from the coupled Maxwell's equations in free space (Eqs. 1-4). The key to decoupling the fields is the vector identity: $$ \nabla \times (\nabla \times \mathbf{F}) = \nabla (\nabla \cdot \mathbf{F}) - \nabla^2 \mathbf{F} \quad \quad (9) $$
Proof for $\mathbf{E}$:
1. Start with Faraday's Law (3):
$$ \nabla \times \mathbf{E} = - \frac{\partial \mathbf{B}}{\partial t} \quad \quad (3) $$2. Take the curl ($\nabla \times$) of both sides:
$$ \nabla \times (\nabla \times \mathbf{E}) = - \nabla \times \left( \frac{\partial \mathbf{B}}{\partial t} \right) = - \frac{\partial}{\partial t} \left( \nabla \times \mathbf{B} \right) \quad \quad (10) $$3. Simplify the LHS using the vector identity (9). Since $\nabla \cdot \mathbf{E} = 0$ (Eq. 1):
$$ \nabla \times (\nabla \times \mathbf{E}) = \nabla (\nabla \cdot \mathbf{E}) - \nabla^2 \mathbf{E} = 0 - \nabla^2 \mathbf{E} = - \nabla^2 \mathbf{E} \quad \quad (11) $$4. Substitute the Ampère-Maxwell Law (4) into the RHS of (10):
$$ - \frac{\partial}{\partial t} \left( \nabla \times \mathbf{B} \right) = - \frac{\partial}{\partial t} \left( \mu_0 \epsilon_0 \frac{\partial \mathbf{E}}{\partial t} \right) = - \mu_0 \epsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2} \quad \quad (12) $$5. Equating (11) and (12) yields the 3D Wave Equation for the Electric Field:
$$ \nabla^2 \mathbf{E} = \mu_0 \epsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2} \quad \quad (13) $$Proof for $\mathbf{B}$:
1. Start with the Ampère-Maxwell Law (4):
$$ \nabla \times \mathbf{B} = \mu_0 \epsilon_0 \frac{\partial \mathbf{E}}{\partial t} \quad \quad (4) $$2. Take the curl ($\nabla \times$) of both sides:
$$ \nabla \times (\nabla \times \mathbf{B}) = \mu_0 \epsilon_0 \nabla \times \left( \frac{\partial \mathbf{E}}{\partial t} \right) = \mu_0 \epsilon_0 \frac{\partial}{\partial t} \left( \nabla \times \mathbf{E} \right) \quad \quad (14) $$3. Simplify the LHS using the vector identity (9). Since $\nabla \cdot \mathbf{B} = 0$ (Eq. 2):
$$ \nabla \times (\nabla \times \mathbf{B}) = \nabla (\nabla \cdot \mathbf{B}) - \nabla^2 \mathbf{B} = 0 - \nabla^2 \mathbf{B} = - \nabla^2 \mathbf{B} \quad \quad (15) $$4. Substitute Faraday's Law (3) into the RHS of (14):
$$ \mu_0 \epsilon_0 \frac{\partial}{\partial t} \left( - \frac{\partial \mathbf{B}}{\partial t} \right) = - \mu_0 \epsilon_0 \frac{\partial^2 \mathbf{B}}{\partial t^2} \quad \quad (16) $$5. Equating (15) and (16) yields the 3D Wave Equation for the Magnetic Field:
$$ \nabla^2 \mathbf{B} = \mu_0 \epsilon_0 \frac{\partial^2 \mathbf{B}}{\partial t^2} \quad \quad (17) $$Proof of Propagation Speed ($c$):
The general form of a linear, homogeneous, 3D wave equation for any field $f$ (which could be displacement, pressure, or electric field component) propagating with velocity $v$ is:
$$ \nabla^2 f = \frac{1}{v^2} \frac{\partial^2 f}{\partial t^2} \quad \quad (18) $$This equation dictates that the spatial curvature ($\nabla^2 f$, representing how the field changes over distance) must be proportional to the temporal acceleration ($\partial^2 f / \partial t^2$, representing how the field changes over time), with the proportionality constant being $1/v^2$.
We compare this general form (18) to the derived Electromagnetic Wave Equation for the electric field (Eq. 13):
$$ \nabla^2 \mathbf{E} = \mu_0 \epsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2} \quad \quad (13) $$By comparing the constants multiplying the time-derivative terms, we can establish the relationship:
$$ \frac{1}{v^2} = \mu_0 \epsilon_0 $$Solving for the speed of propagation, $v$:
$$ v = \frac{1}{\sqrt{\mu_0 \epsilon_0}} \quad \quad (19) $$This result shows that EM disturbances propagate at a constant speed determined solely by fundamental electromagnetic constants: the permeability of free space ($\mu_0$) and the permittivity of free space ($\epsilon_0$). When Maxwell calculated this speed using the experimentally known values for $\mu_0$ and $\epsilon_0$:
$$ v = \frac{1}{\sqrt{(4\pi \times 10^{-7} \, \text{N/A}^2) \times (8.854 \times 10^{-12} \, \text{F/m})}} \approx 2.9979 \times 10^8 \, \text{m/s} $$This value was precisely the known speed of light, $c$. This mathematical identity led Maxwell to the profound conclusion that **light itself is an electromagnetic wave**.
$$ v = c = \frac{1}{\sqrt{\mu_0 \epsilon_0}} \quad \quad (20) $$To simplify the four coupled Maxwell equations (5)-(8) into two independent wave equations for the sources ($\rho$ and $\mathbf{J}$), we introduce the scalar potential ($\phi$) and the magnetic vector potential ($\mathbf{A}$).
The fields are defined in terms of the potentials as:
$$ \mathbf{B} = \nabla \times \mathbf{A} \quad \quad (21) $$ $$ \mathbf{E} = - \nabla \phi - \frac{\partial \mathbf{A}}{\partial t} \quad \quad (22) $$The fields $\mathbf{E}$ and $\mathbf{B}$ are invariant under certain transformations of the potentials (Gauge freedom). We impose a **Gauge Condition** to decouple the resulting complex equations. The **Lorenz Gauge** condition is chosen for this purpose:
$$ \nabla \cdot \mathbf{A} + \mu_0 \epsilon_0 \frac{\partial \phi}{\partial t} = 0 \quad \quad (23) $$Applying (21) and (22) into the General Maxwell's Equations (5) and (8), and utilizing the Lorenz Gauge (23), results in the Inhomogeneous Wave Equations:
$$ \nabla^2 \mathbf{A} - \mu_0 \epsilon_0 \frac{\partial^2 \mathbf{A}}{\partial t^2} = - \mu_0 \mathbf{J} \quad \quad (24) $$ $$ \nabla^2 \phi - \mu_0 \epsilon_0 \frac{\partial^2 \phi}{\partial t^2} = - \frac{\rho}{\epsilon_0} \quad \quad (25) $$The operator $\nabla^2 - \mu_0 \epsilon_0 \frac{\partial^2}{\partial t^2}$ is the **D'Alembertian operator** ($\Box^2$). These equations show that potentials, and thus the fields they create, are sourced by $\mathbf{J}$ and $\rho$, and propagate away from the sources at speed $c$.
The Poynting theorem describes the conservation of electromagnetic energy. It relates the work done on charges ($\mathbf{E} \cdot \mathbf{J}$), the change in stored field energy ($u$), and the flow of energy across the boundaries ($\mathbf{S}$).
Proof:
1. Start with the rate of work done by the fields per unit volume, $\mathbf{E} \cdot \mathbf{J}$. Substitute $\mathbf{J}$ from the General Ampère-Maxwell Law (8):
$$ \mathbf{E} \cdot \mathbf{J} = \frac{1}{\mu_0} \mathbf{E} \cdot (\nabla \times \mathbf{B}) - \epsilon_0 \mathbf{E} \cdot \frac{\partial \mathbf{E}}{\partial t} \quad \quad (26) $$2. Use the vector identity $\nabla \cdot (\mathbf{E} \times \mathbf{B}) = \mathbf{B} \cdot (\nabla \times \mathbf{E}) - \mathbf{E} \cdot (\nabla \times \mathbf{B})$ to isolate $\mathbf{E} \cdot (\nabla \times \mathbf{B})$.
3. Substitute the resulting expression and Faraday's Law (7), $\nabla \times \mathbf{E} = - \partial \mathbf{B} / \partial t$, into (26). Use the rule $\mathbf{F} \cdot \frac{\partial \mathbf{F}}{\partial t} = \frac{1}{2} \frac{\partial}{\partial t} (\mathbf{F}^2)$ for both fields:
$$ \mathbf{E} \cdot \mathbf{J} = - \frac{\partial}{\partial t} \left( \frac{1}{2}\epsilon_0 E^2 + \frac{1}{2\mu_0} B^2 \right) - \nabla \cdot \left( \frac{1}{\mu_0} \mathbf{E} \times \mathbf{B} \right) \quad \quad (27) $$This is the differential form of the Poynting Theorem. The terms are defined as: $$ \mathbf{E} \cdot \mathbf{J} = \text{Power dissipated per unit volume (Joule heating).} $$ $$ u = \frac{1}{2}\epsilon_0 E^2 + \frac{1}{2\mu_0} B^2 = \text{Total energy density of the EM field.} \quad \quad (28) $$ $$ \mathbf{S} = \frac{1}{\mu_0} \mathbf{E} \times \mathbf{B} = \text{Poynting Vector (EM power flow per unit area).} \quad \quad (29) $$
The theorem (27) can be rewritten as: $$ \nabla \cdot \mathbf{S} + \frac{\partial u}{\partial t} = - \mathbf{E} \cdot \mathbf{J} \quad \quad (30) $$
An EM wave is a **transverse wave**, meaning the fields oscillate perpendicular to the direction of propagation ($\mathbf{k}$). The fields $\mathbf{E}$ and $\mathbf{B}$ are also perpendicular to each other, forming a right-handed system $(\mathbf{E}, \mathbf{B}, \mathbf{k})$.
By substituting the general plane wave solution $\mathbf{E}(\mathbf{r}, t) = \mathbf{E}_0 e^{i(\mathbf{k} \cdot \mathbf{r} - \omega t)}$ into Gauss's Law ($\nabla \cdot \mathbf{E} = 0$, Eq. 1): $$ \nabla \cdot \mathbf{E} = i (\mathbf{k} \cdot \mathbf{E}_0) e^{i(\mathbf{k} \cdot \mathbf{r} - \omega t)} = 0 \quad \quad (31) $$ Since $i$ and the exponential term are non-zero, this necessitates $\mathbf{k} \cdot \mathbf{E}_0 = 0$. This confirms $\mathbf{E}$ is perpendicular to $\mathbf{k}$. The same process for $\mathbf{B}$ (Eq. 2) proves $\mathbf{k} \cdot \mathbf{B} = 0$, thus $\mathbf{B} \perp \mathbf{k}$.
Substituting the plane wave solutions into Faraday's Law (3), $\nabla \times \mathbf{E} = - \frac{\partial \mathbf{B}}{\partial t}$, yields:
$$ i \mathbf{k} \times \mathbf{E} = i \omega \mathbf{B} \quad \quad (32) $$Which simplifies to:
$$ \mathbf{B} = \frac{1}{\omega} (\mathbf{k} \times \mathbf{E}) \quad \quad (33) $$Equation (33) proves that $\mathbf{B}$ is perpendicular to both $\mathbf{k}$ and $\mathbf{E}$. Since $\omega/k = c$ (the phase speed), the magnitude relation is:
$$ B = \frac{k}{\omega} E = \frac{1}{c} E \quad \quad (34) $$ $$ \frac{E}{B} = c \quad \quad (35) $$In free space, the energy is equally distributed between the electric and magnetic fields. Since $u_E = \frac{1}{2}\epsilon_0 E^2$ and $u_B = \frac{1}{2\mu_0} B^2$, substituting $E = cB$ and $c^2 = 1/(\mu_0 \epsilon_0)$ into $u_E$:
$$ u_E = \frac{1}{2} \epsilon_0 (cB)^2 = \frac{1}{2} \epsilon_0 \left( \frac{1}{\mu_0 \epsilon_0} \right) B^2 = \frac{1}{2\mu_0} B^2 = u_B \quad \quad (36) $$The ratio of the electric field magnitude to the magnetic field magnitude, scaled by $\mu_0$, is the **Characteristic Impedance of Free Space** ($Z_0$):
$$ Z_0 = \mu_0 \frac{E}{B} = \mu_0 c = \sqrt{\frac{\mu_0}{\epsilon_0}} \approx 377 \, \Omega \quad \quad (37) $$The flow of energy associated with an EM wave is quantified by its **Intensity ($I$)**, which is the time-averaged magnitude of the Poynting Vector $\mathbf{S}$.
Derivation of $I_{avg}$:
For a plane wave, $\mathbf{E}$ and $\mathbf{B}$ are in phase, perpendicular, and the fields vary as $\cos(kx-\omega t)$. The instantaneous Poynting magnitude is:
$$ S(t) = \frac{1}{\mu_0} |\mathbf{E} \times \mathbf{B}| = \frac{E B}{\mu_0} = \frac{E (E/c)}{\mu_0} = \frac{E^2}{\mu_0 c} $$If $E = E_0 \cos(kx-\omega t)$, then $E^2 = E_0^2 \cos^2(kx-\omega t)$. The time average of $\cos^2(\theta)$ over one cycle is $\langle \cos^2(\theta) \rangle = 1/2$.
$$ I_{avg} = \langle S \rangle = \frac{\langle E^2 \rangle}{\mu_0 c} = \frac{E_0^2 / 2}{\mu_0 c} $$Using $c^2 = 1/(\mu_0 \epsilon_0)$, so $\frac{1}{\mu_0 c} = c \epsilon_0$:
$$ I_{avg} = \frac{1}{2} E_0^2 (\epsilon_0 c) = \frac{1}{2} c \epsilon_0 E_0^2 \quad \quad (38) $$Because EM waves carry momentum, they exert a pressure on any surface they strike. The pressure ($P_{rad}$) is the average momentum flux and is directly related to the wave intensity ($I_{avg}$):
This effect, though tiny, confirms the momentum-carrying nature of light and is critical in astrophysical processes.
The EM spectrum covers a vast range of waves, all propagating at speed $c$ in a vacuum and satisfying the derived wave equations. They are classified by frequency ($f$) or wavelength ($\lambda$), related by $c = f \lambda$.
[Image of the complete Electromagnetic Spectrum]Light is simply the narrow band of the spectrum visible to the human eye. The mathematical framework developed here applies universally across the entire spectrum, from low-frequency radio waves to high-frequency gamma rays.
Formulas: $E_y = E_0 \sin(kx - \omega t)$, $k = \omega/c$, $f = \omega/(2\pi)$.
Solution: From the equation, $k = 1.0 \times 10^7 \, \text{m}^{-1}$.
$$ \omega = k c = (1.0 \times 10^7 \, \text{m}^{-1}) (3.0 \times 10^8 \, \text{m/s}) = 3.0 \times 10^{15} \, \text{rad/s} $$ $$ f = \frac{\omega}{2\pi} = \frac{3.0 \times 10^{15} \, \text{rad/s}}{2\pi} \approx 4.77 \times 10^{14} \, \text{Hz} $$Formula: $E_0 / B_0 = c \quad \text{(Eq. 35)}$.
Solution: The ratio of the field magnitudes equals the speed of light.
$$ E_0 = B_0 c = (5.0 \times 10^{-11} \, \text{T}) (3.0 \times 10^8 \, \text{m/s}) = 15.0 \times 10^{-3} \, \text{V/m} = 15.0 \, \text{mV/m} $$Formula: $\lambda = c/f$.
Solution: Convert frequency to SI units: $f = 98.1 \, \text{MHz} = 98.1 \times 10^6 \, \text{Hz}$.
$$ \lambda = \frac{3.0 \times 10^8 \, \text{m/s}}{98.1 \times 10^6 \, \text{Hz}} \approx 3.06 \, \text{m} $$Formula: $I_{avg} = \frac{1}{2} c \epsilon_0 E_0^2 \quad \text{(Eq. 38)}$.
Solution: Substitute values: $\epsilon_0 \approx 8.854 \times 10^{-12} \, \text{F/m}$.
$$ I_{avg} = \frac{1}{2} (3.0 \times 10^8) (8.854 \times 10^{-12}) (10.0)^2 $$ $$ I_{avg} \approx 13.28 \, \text{W/m}^2 $$Formulas: $\mathbf{S} \propto \mathbf{E} \times \mathbf{B} \quad \text{(Eq. 29)}$, $\frac{E}{B} = c \quad \text{(Eq. 35)}$.
Solution: Propagation direction is $+z$, so $\mathbf{\hat{k}} = \mathbf{\hat{z}}$. Magnetic field is $\mathbf{\hat{B}} = \mathbf{\hat{x}}$. We must have $\mathbf{\hat{z}} \propto \mathbf{\hat{E}} \times \mathbf{\hat{x}}$. This requires $\mathbf{\hat{E}} = -\mathbf{\hat{y}}$ (since $(-\mathbf{\hat{y}}) \times (\mathbf{\hat{x}}) = \mathbf{\hat{z}}$). Thus, $\mathbf{E}$ is in the $-\mathbf{\hat{y}}$ direction.
$$ E = B c = (2.0 \times 10^{-7} \, \text{T}) (3.0 \times 10^8 \, \text{m/s}) = 60.0 \, \text{V/m} $$ $$ \mathbf{E} = -60.0 \, \mathbf{\hat{y}} \, \text{V/m} $$Formula: $v = \frac{1}{\sqrt{\mu \epsilon}} \quad \text{(Generalization of Eq. 20)}$.
Solution: In the medium, $\mu = \mu_0$ and $\epsilon = 4\epsilon_0$.
$$ v = \frac{1}{\sqrt{\mu_0 (4\epsilon_0)}} = \frac{1}{2 \sqrt{\mu_0 \epsilon_0}} = \frac{c}{2} $$The wave speed is $v = c/2$. The ratio of fields $E/B$ equals the wave speed $v$ in any non-conducting medium. Thus:
$$ \frac{E}{B} = v = \frac{c}{2} $$Formulas: $u_E = \frac{1}{2}\epsilon_0 E^2$, $u_B = \frac{1}{2\mu_0} B^2$, $E/B = c \quad \text{(Eq. 35)}$, $c^2 = 1/(\mu_0 \epsilon_0) \quad \text{(From Eq. 20)}$.
Proof: We wish to show $u_E = u_B$. Substitute $E = cB$ into the expression for $u_E$:
$$ u_E = \frac{1}{2} \epsilon_0 E^2 = \frac{1}{2} \epsilon_0 (cB)^2 = \frac{1}{2} \epsilon_0 c^2 B^2 $$Now substitute the expression for $c^2$:
$$ u_E = \frac{1}{2} \epsilon_0 \left( \frac{1}{\mu_0 \epsilon_0} \right) B^2 = \frac{1}{2} \frac{1}{\mu_0} B^2 = u_B $$Formula: $F = P_{rad} A$, $P_{rad} = I_{avg}/c$ (absorbing surface) $\quad \text{(Eq. 39)}$.
Solution: The total force exerted on the absorbing surface is $F = P_{rad} A = (I_{avg}/c) A$.
$$ F = \frac{I_{avg} A}{c} = \frac{(1.0 \, \text{W/m}^2) (0.01 \, \text{m}^2)}{3.0 \times 10^8 \, \text{m/s}} $$ $$ F \approx 3.33 \times 10^{-11} \, \text{N} $$Formulas: $\mathbf{B} = \frac{1}{\omega} (\mathbf{k} \times \mathbf{E}) \quad \text{(Eq. 33)}$, $\mathbf{E} \perp \mathbf{B}$, $E/B = c \quad \text{(Eq. 35)}$.
Solution: The argument $\cos(\omega t + k x)$ means the wave propagates in the **$-x$ direction**. Thus, $\mathbf{\hat{k}} = -\mathbf{\hat{x}}$. The magnetic field direction is $\mathbf{\hat{B}} = \mathbf{\hat{z}}$.
The energy flow (Poynting) relation is $\mathbf{\hat{k}} \propto \mathbf{\hat{E}} \times \mathbf{\hat{B}}$.
$$ -\mathbf{\hat{x}} \propto \mathbf{\hat{E}} \times \mathbf{\hat{z}} $$To satisfy the right-hand rule, $\mathbf{\hat{E}}$ must be in the $-\mathbf{\hat{y}}$ direction, since $(-\mathbf{\hat{y}}) \times (\mathbf{\hat{z}}) = -\mathbf{\hat{x}}$. The magnitude is $E = c B_0$.
$$ \mathbf{E} = - c B_0 \mathbf{\hat{y}} \cos(\omega t + k x) $$Formulas: $P_{rad} = I_{avg}/c$ (absorbing surface) $\quad \text{(Eq. 39)}$, $F = P_{rad} A$.
Solution:
Radiation Pressure:
$$ P_{rad} = \frac{I_{avg}}{c} = \frac{1.0 \times 10^3 \, \text{W/m}^2}{3.0 \times 10^8 \, \text{m/s}} \approx 3.33 \times 10^{-6} \, \text{N/m}^2 \text{ (or Pa)} $$Total Force:
$$ F = P_{rad} A = (3.33 \times 10^{-6} \, \text{N/m}^2) (0.50 \, \text{m}^2) \approx 1.67 \times 10^{-6} \, \text{N} $$