从Maxwell方程组到波动方程

Maxwell方程组是十九世纪最伟大的公式，代表了传统物理学人对公式美学的孜孜追求，也影响了无数后来者的物理美学品味。

回顾历史，当1864年，Maxwell发出那篇著名的《电磁场的动力学理论》时，实则列出了二十个公式，以总结前人的物理学成果，我们将分量公式合并为矢量，可以得到八个式子，即

名称	公式	物理意义
总电流定律	${\mathbf{J}}_{tot}=\mathbf{J}+\frac{\partial \mathbf{D}}{\partial t}$	总电流等于传导电流加位移电流
磁场方程	$\mu \mathbf{H}=\nabla \times \mathbf{A}$	磁场强度为矢势的旋度
安培环路定理	$\nabla \times \mathbf{H}={\mathbf{J}}_{tot}$	变化磁场产生电流
洛伦兹力	$\mathbf{E}=\mu \mathbf{v}\times \mathbf{H}-\frac{\partial \mathbf{A}}{\partial t}-\nabla \varphi$	磁场中运动电荷产生指向法线方向的力
电弹性方程	$\mathbf{E}=\frac{1}{ϵ}\mathbf{D}$	场强和电位移矢量的关系
欧姆定律	$\mathbf{E}=\frac{1}{\sigma }\mathbf{J}$	电流和电势的关系
高斯定律	$\nabla \cdot \mathbf{D}=\rho$	电荷产生电场
连续性方程	$\nabla \cdot \mathbf{J}=-\frac{\partial \rho }{\partial t}$	电荷变化产生传导电流

以上符号分别表示

符号	物理意义	符号	物理意义
$\mathbf{H}$	磁场强度	$\mathbf{D}$	电位移矢量
		$\mathbf{E}$	电场强度
$\mathbf{J}$	传导电流密度	$\rho$	自由电荷密度
$\mathbf{J}tot$	总电流密度
$\mathbf{A}$	矢量势	$\varphi$	标势
$\mu$	磁导率	$\sigma$	电导率
$ϵ$	电容率

二十年后，Heaviside对这二十个公式进行重新编排，得到了我们熟悉的形式，并将其命名为麦克斯韦方程组：

$$\begin{array}{rl}\nabla \times \mathbf{E}& =-\frac{\partial \mathbf{B}}{\partial t}\\ \nabla \times \mathbf{H}& =\mathbf{J}+\frac{\partial \mathbf{D}}{\partial t}\\ \nabla \cdot \mathbf{D}& =\rho \\ \nabla \cdot \mathbf{B}& =0\end{array}$$

而$\mathbf{D},\mathbf{B},\mathbf{J}$的定义，则作为物质方程而存在

$$\begin{array}{rl}& \mathbf{D}={\epsilon }_0\mathbf{E}+\mathbf{P}={\epsilon }_0\epsilon \mathbf{E}=\epsilon E\\ & \mathbf{B}=\mu H={\mu }_0\mathbf{H}={\mu }_0(\mathbf{H}+\mathbf{M})\\ & \mathbf{J}=\sigma \mathbf{E}\end{array}$$

其中，$\mathbf{P}$为介质在外场作用下的电极化强度矢量；$\mathbf{M}$为磁化强度矢量。

对于真空而言，$\epsilon ={\epsilon }_0,\mu ={\mu }_0$皆为标量，且真空中点电荷与电流均为0，所以可以进一步写为

$$\begin{array}{rlr}\nabla \times \mathbf{E}& =-\frac{\partial \mathbf{B}}{\partial t}& \nabla \cdot \mathbf{D}=0\\ \nabla \times \mathbf{B}& ={\epsilon }_0{\mu }_0\frac{\partial \mathbf{E}}{\partial t}& \nabla \cdot \mathbf{B}=0\end{array}$$

对上式中左侧两个旋度公式再取旋度，得到

$$\begin{array}{r}\nabla \times (\nabla \times \mathbf{E})=-\frac{\partial }{\partial t}(\nabla \times \mathbf{B})=-{\epsilon }_0{\mu }_0\frac{{\partial }^2\mathbf{E}}{\partial t^2}\\ \nabla \times (\nabla \times \mathbf{B})={\epsilon }_0{\mu }_0\frac{\partial }{\partial t}(\nabla \times \mathbf{E})=-{\epsilon }_0{\mu }_0\frac{{\partial }^2\mathbf{B}}{\partial t^2}\end{array}$$

很多人对矢量的计算方法不太熟悉，所以下面着重推导一下$\nabla \times (\nabla \times \mathbf{E})$。通过上标来表示空间分量，下标表示对空间坐标的导数，例如记空间矢量$\stackrel{\rightharpoonup }{r}=(r^x,r^y,r^z)$，场强$\mathbf{E}=(E^x,E^y,E^z)$，记$\frac{\partial E^x}{\partial y}=E_y^x$，则

$$\begin{array}{rl}\nabla \times (\nabla \times \mathbf{E})& =\nabla \times \left[\begin{array}{ccc}\stackrel{\rightharpoonup }{n^x}& \stackrel{\rightharpoonup }{n^y}& \stackrel{\rightharpoonup }{n^z}\\ \frac{\partial }{\partial x}& \frac{\partial }{\partial y}& \frac{\partial }{\partial z}\\ E^x& E^y& E^z\end{array}\right]\\ & =\nabla \times [E_y^z-E_z^y,-(E_x^z-E_z^x),E_x^y-E_y^x]\\ & =\left[\begin{array}{ccc}\stackrel{\rightharpoonup }{n^x}& \stackrel{\rightharpoonup }{n^y}& \stackrel{\rightharpoonup }{n^z}\\ \frac{\partial }{\partial x}& \frac{\partial }{\partial y}& \frac{\partial }{\partial z}\\ E_y^z-E_z^y& E_z^x-E_x^z& E_x^y-E_y^x\end{array}\right]\\ & ={\left[\begin{array}{c}(E_{xy}^y-E_{yy}^x)-(E_{zz}^x-E_{xz}^z)\\ -[(E_{xx}^y-E_{yx}^x)-(E_{yz}^z-E_{zz}^y)]\\ (E_{zx}^x-E_{xx}^z)-(E_{yy}^z-E_{zy}^y)\end{array}\right]}^T\\ & =-{\left[\begin{array}{c}{E}_{zz}^x+E_{yy}^x+E_{xx}^x\\ E_{zz}^y+E_{xx}^y+E_{yy}^y\\ E_{zz}^z+E_{yy}^z+E_{xx}^z\end{array}\right]}^T+{\left[\begin{array}{c}{E}_{xx}^x+E_{yx}^y+E_{xz}^z\\ E_{yy}^y+E_{xy}^x+E_{yz}^z\\ E_{xz}^x+E_{yz}^y+E_{zz}^z\end{array}\right]}^T\\ & =-{\nabla }^2\mathbf{E}+\nabla (\nabla \cdot \mathbf{E})\end{array}$$

其中，$\nabla \mathbf{E}=0$，所以可得到波动方程

$$\begin{array}{r}({\nabla }^2-{\epsilon }_0{\mu }_0\frac{{\partial }^2}{\partial t^2})\mathbf{E}=0\\ ({\nabla }^2-{\epsilon }_0{\mu }_0\frac{{\partial }^2}{\partial t^2})\mathbf{B}=0\end{array}$$

令$c=\frac{1}{\sqrt{{\epsilon }_0{\mu }_0}}$，并定义达朗贝尔算子$\square =-\frac{1}{c^2}\frac{{\partial }^2}{\partial t^2}+{\nabla }^2$，则波动方程可以写为

$$\square \mathbf{B}=0,\square \mathbf{E}=0$$