波动方程的推导与解的物理意义
波动光学是光学的一个重要分支，它将光视为在空间中传播的电磁波，用以解释几何光学无法涵盖的干涉、衍射和偏振等复杂现象。在深入探讨惠更斯原理、光的干涉、光的衍射、光的色散与光的偏振等具体子主题之前，我们需要建立一个统一的数学与物理框架。本文将从麦克斯韦方程组出发，推导描述光波传播的通用波动方程，并深入剖析其数学解法及对应的物理意义。
波动光学宏观理论的基石是经典电动力学中的麦克斯韦方程组。在没有自由电荷和自由电流的各向同性、均匀介质中，麦克斯韦方程组的微分形式可以写为：
- $\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 \epsilon \frac{\partial \mathbf{E}}{\partial t}$ （安培-麦克斯韦定律）
其中，$\mathbf{E}$ 为电场强度，$\mathbf{B}$ 为磁感应强度，$\epsilon$ 为介质的电容率（介电常数），$\mu$ 为介质的磁导率。
为了推导电场的波动方程，我们可以对法拉第电磁感应定律（方程3）两边取旋度：
$$\nabla \times (\nabla \times \mathbf{E}) = \nabla \times \left(-\frac{\partial \mathbf{B}}{\partial t}\right)$$
利用向量恒等式 $\nabla \times (\nabla \times \mathbf{E}) = \nabla(\nabla \cdot \mathbf{E}) - \nabla^2 \mathbf{E}$，并结合电场高斯定律 $\nabla \cdot \mathbf{E} = 0$，左边简化为 $-\nabla^2 \mathbf{E}$。
对于右边，交换空间微商与时间微商的顺序，并代入安培-麦克斯韦定律（方程4）：
$$\nabla \times \mathbf{B} = \mu \epsilon \frac{\partial \mathbf{E}}{\partial t}$$
最终我们可以得到电场强度的矢量波动方程：
$$\nabla^2 \mathbf{E} - \mu \epsilon \frac{\partial^2 \mathbf{E}}{\partial t^2} = 0$$
同理，也可以推导出磁场强度的波动方程：
$$\nabla^2 \mathbf{B} - \mu \epsilon \frac{\partial^2 \mathbf{B}}{\partial t^2} = 0$$
波动方程的标准形式与波速
为了更直观地理解波动方程，我们可以将其写成标量形式（针对电场的某一分量 $u$）：
$$\nabla^2 u - \frac{1}{v^2} \frac{\partial^2 u}{\partial t^2} = 0$$
通过对比系数，我们可以定义光在介质中的传播速度 $v$：
$$v = \frac{1}{\sqrt{\mu \epsilon}}$$
在真空中，介电常数和磁导率分别为 $\epsilon_0$ 和 $\mu_0$，此时光速为：
$$c = \frac{1}{\sqrt{\mu_0 \epsilon_0}}$$
介质中的光速 $v$ 与真空中光速 $c$ 的比值定义为介质的绝对折射率 $n$：
$$n = \frac{c}{v} = \sqrt{\frac{\mu \epsilon}{\mu_0 \epsilon_0}}$$
对于大多数光学透明介质，其相对磁导率 $\mu_r \approx 1$，因此折射率主要由相对电容率 $\epsilon_r$ 决定，即 $n \approx \sqrt{\epsilon_r}$。
波动方程的解及其物理意义
波动方程是一个二阶线性偏微分方程，在不同的坐标系和边界条件下有多种形式的解。最基础且最重要的解是简谐平面波解。
1. 简谐平面波的数学表达
假设光波沿任意方向 $\mathbf{k}$ 传播，其复振幅形式的解可以表示为：
$$\mathbf{E}(\mathbf{r}, t) = \mathbf{E}_0 \cos(\mathbf{k} \cdot \mathbf{r} - \omega t + \phi_0)$$
或者使用复数形式（在物理计算中更为常用）：
$$\mathbf{E}(\mathbf{r}, t) = \text{Re} \left[ \mathbf{E}_0 e^{i(\mathbf{k} \cdot \mathbf{r} - \omega t)} \right]$$
2. 核心物理参数解析
- 振幅 $\mathbf{E}_0$：决定了光波的振动强弱。在物理上，光波的强度（辐照度）与电场振幅的平方成正比，$I \propto |\mathbf{E}_0|^2$。
- 波矢 $\mathbf{k}$：其方向代表光波的传播方向，其大小（波数）为 $k = |\mathbf{k}| = \frac{2\pi}{\lambda}$（其中 $\lambda$ 为介质中的波长）。
- 角频率 $\omega$：由光源的振动频率决定，$\omega = 2\pi \nu$。频率在光波穿过不同介质时保持不变。
- 相位 $\Phi = \mathbf{k} \cdot \mathbf{r} - \omega t + \phi_0$：决定了空间中某一点在某一时刻的振动状态。
3. 其他常见类型的解
除了平面波，波动方程在数学处理上还允许其他重要的特解，以适应不同的物理场景：
- 球面波：适用于点光源发出的光，振幅随传播距离 $r$ 的增加而衰减（与 $1/r$ 成正比）。
- 高斯光束：激光器输出的光束形态，其能量在横截面上呈高斯分布，是描述实际激光传播的最优近似解。
波动光学建模的全局视角
在掌握了通用波动方程及其解的物理意义后，我们可以俯瞰整个波动光学的内容体系。不同的光学现象本质上都是波动方程在特定条件下的表现：
- 惠更斯原理：提供了从已知波前构建未来波前的几何方法，是理解衍射现象的直观基础。
- 光的干涉：两列或多列满足相干条件的波叠加，导致波动方程解的振幅重新分布（极大与极小）。
- 光的衍射：波在传播过程中遇到障碍物时偏离直线传播路径的现象，本质上是波的连续次级子波叠加（数学上表现为基尔霍夫衍射积分，源于波动方程的边界值问题）。
- 光的色散：由于介质参数 $\epsilon$（从而折射率 $n$）随频率 $\omega$ 变化，导致不同频率成分的波速不同。
- 光的偏振：体现了电场矢量 $\mathbf{E}$ 振动方向相对于传播方向的横波特性。
通过建立严谨的波动方程数学模型，我们得以将上述复杂的物理现象统一在麦克斯韦电磁理论的大厦之下，为后续的定量计算和实验设计奠定了坚实的理论基础。