推导电磁波波动方程的过程
电磁波波动方程的推导是经典电磁学的巅峰之作，它将电场与磁场统一在同一个数学框架下，并预言了光本质上是一种电磁波。要完成这一推导，我们需要从麦克斯韦方程组出发，在真空中（无电荷、无电流）通过矢量分析将其转化为波动方程的形式。
在推导波动方程之前，我们必须定义研究的物理环境。在真空中，电荷密度 $\rho = 0$ 且电流密度 $\mathbf{J} = 0$。此时，麦克斯韦方程组简化为以下形式：
- 高斯电场定律：$\nabla \cdot \mathbf{E} = 0$
- 高斯磁场定律：$\nabla \cdot \mathbf{B} = 0$
- 法拉第电磁感应定律：$\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}$
- 安培-麦克斯韦定律：$\nabla \times \mathbf{B} = \mu_0 \epsilon_0 \frac{\partial \mathbf{E}}{\partial t}$
其中，$\mathbf{E}$ 为电场强度，$\mathbf{B}$ 为磁感应强度，$\epsilon_0$ 为真空介电常数，$\mu_0$ 为真空磁导率。
数学工具：矢量恒等式
在推导过程中，我们需要用到一个关键的矢量分析恒等式，即“旋度的旋度”公式：
$$\nabla \times (\nabla \times \mathbf{A}) = \nabla(\nabla \cdot \mathbf{A}) - \nabla^2 \mathbf{A}$$
这个公式将一个矢量的二阶旋度运算转化为散度的梯度与拉普拉斯算子（$\nabla^2$）的差。
电场波动方程的推导步骤
我们将通过对法拉第定律进行空间求旋的操作，逐步导出电场的波动方程。
第一步：对法拉第定律求旋
取麦克斯韦方程组中的第三式 $\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}$，在等式两边同时作用旋度算子 $\nabla \times$：
$$\nabla \times (\nabla \times \mathbf{E}) = \nabla \times \left( -\frac{\partial \mathbf{B}}{\partial t} \right)$$
第二步：交换空间导数与时间导数
由于空间坐标和时间是独立变量，我们可以将旋度算子 $\nabla \times$ 与时间偏导数 $\frac{\partial}{\partial t}$ 互换位置：
$$\nabla \times (\nabla \times \mathbf{E}) = -\frac{\partial}{\partial t} (\nabla \times \mathbf{B})$$
第三步：代入安培-麦克斯韦定律
此时，等式右边的 $\nabla \times \mathbf{B}$ 可以用麦克斯韦方程组的第四式 $\nabla \times \mathbf{B} = \mu_0 \epsilon_0 \frac{\partial \mathbf{E}}{\partial t}$ 来替换：
$$\nabla \times (\nabla \times \mathbf{E}) = -\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}$$
第四步：利用矢量恒等式简化左侧
利用前文提到的恒等式 $\nabla \times (\nabla \times \mathbf{E}) = \nabla(\nabla \cdot \mathbf{E}) - \nabla^2 \mathbf{E}$，将左侧展开：
$$\nabla(\nabla \cdot \mathbf{E}) - \nabla^2 \mathbf{E} = -\mu_0 \epsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2}$$
第五步：应用真空条件
根据真空中的高斯电场定律 $\nabla \cdot \mathbf{E} = 0$，左侧的第一项 $\nabla(\nabla \cdot \mathbf{E})$ 为零。于是方程简化为：
$$-\nabla^2 \mathbf{E} = -\mu_0 \epsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2}$$
整理得：
$$\nabla^2 \mathbf{E} - \mu_0 \epsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2} = 0$$
磁场波动方程的推导
磁场 $\mathbf{B}$ 的推导过程与电场完全对称。我们从安培-麦克斯韦定律出发：
- 对 $\nabla \times \mathbf{B} = \mu_0 \epsilon_0 \frac{\partial \mathbf{E}}{\partial t}$ 两边求旋。
- 利用 $\nabla \times (\nabla \times \mathbf{B}) = \nabla(\nabla \cdot \mathbf{B}) - \nabla^2 \mathbf{B}$。
- 代入法拉第定律 $\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}$。
- 利用 $\nabla \cdot \mathbf{B} = 0$。
最终得到磁场的波动方程：
$$\nabla^2 \mathbf{B} - \mu_0 \epsilon_0 \frac{\partial^2 \mathbf{B}}{\partial t^2} = 0$$
物理意义与结论
1. 波速的确定
上述两个方程在数学形式上完全一致，都是标准的三维波动方程。其通用形式为 $\nabla^2 \Psi - \frac{1}{v^2} \frac{\partial^2 \Psi}{\partial t^2} = 0$，其中 $v$ 为波速。通过对比可以得出：
$$\frac{1}{v^2} = \mu_0 \epsilon_0 \implies v = \frac{1}{\sqrt{\mu_0 \epsilon_0}}$$
将 $\mu_0$ 和 $\epsilon_0$ 的数值代入，计算结果约为 $2.9979 \times 10^8 \text{ m/s}$，这恰好与实验测得的光速 $c$ 一致。
2. 电磁波的特性
从波动方程及其推导过程中，我们可以得出电磁波的几个核心物理图像：
- 自持性：变化的磁场产生电场，变化的电场产生磁场，两者在空间中交替产生，无需介质即可传播。
- 横波性质：由于 $\nabla \cdot \mathbf{E} = 0$ 和 $\nabla \cdot \mathbf{B} = 0$，电场 $\mathbf{E}$ 和磁场 $\mathbf{B}$ 始终垂直于波的传播方向 $\mathbf{k}$。
- 正交关系：$\mathbf{E}$ 和 $\mathbf{B}$ 彼此正交，且相位相同。
总结
推导电磁波波动方程的过程实际上是将麦克斯韦方程组中的“耦合关系”通过二阶微分运算转化为“独立演化关系”的过程。这一推导不仅在数学上证明了电磁场的波动性，更在物理上揭示了光作为电磁波的本质，为后来的相对论和量子力学奠定了基础。