动量守恒方程的积分形式推导
在经典电动力学中，动量守恒定律是描述物质与电磁场相互作用的核心原理之一。与纯力学系统不同，电磁系统中的动量不仅存在于运动的粒子（机械动量）中，还存在于电磁场本身（场动量）中。为了在宏观尺度上描述这种复杂的相互作用，我们需要从微分形式的力密度出发，通过数学变换推导出动量守恒的积分形式。
在进行推导之前，我们需要明确以下几个关键物理量的定义：
机械动量密度 ($\mathbf{g}_{mech}$)：
对于电荷密度为 $\rho_m$、速度为 $\mathbf{v}$ 的物质，其机械动量密度定义为：
$$\mathbf{g}_{mech} = \rho_m \mathbf{v}$$电磁场动量密度 ($\mathbf{g}_{em}$)：
电磁场携带的动量密度由坡印廷矢量 $\mathbf{S} = \frac{1}{\mu_0} (\mathbf{E} \times \mathbf{B})$ 与光速的关系决定，表达式为：
$$\mathbf{g}_{em} = \epsilon_0 (\mathbf{E} \times \mathbf{B})$$洛伦兹力密度 ($\mathbf{f}_{em}$)：
电磁场对电荷和电流施加的单位体积力为：
$$\mathbf{f}_{em} = \rho \mathbf{E} + \mathbf{J} \times \mathbf{B}$$
其中 $\rho$ 是电荷密度，$\mathbf{J}$ 是电流密度。麦克斯韦应力张量 ($\mathbf{T}$)：
这是一个二阶对称张量，用于描述电磁场在空间中传递动量的压力和剪切力。其分量 $T_{ij}$ 定义为：
$$T_{ij} = \epsilon_0 \left( E_i E_j - \frac{1}{2} \delta_{ij} E^2 \right) + \frac{1}{\mu_0} \left( B_i B_j - \frac{1}{2} \delta_{ij} B^2 \right)$$
麦克斯韦应力张量的微分关系
动量守恒推导的关键在于建立力密度 $\mathbf{f}_{em}$ 与应力张量 $\mathbf{T}$ 之间的数学联系。通过结合麦克斯韦方程组（高斯定律、安培-麦克斯韦定律等），可以证明洛伦兹力密度可以改写为：
$$\mathbf{f}_{em} = \nabla \cdot \mathbf{T} - \epsilon_0 \mu_0 \frac{\partial \mathbf{S}}{\partial t}$$
这个等式具有深刻的物理意义：
- $\nabla \cdot \mathbf{T}$ 项代表了电磁场通过表面施加的力（类似于力学中的应力梯度）。
- $\epsilon_0 \mu_0 \frac{\partial \mathbf{S}}{\partial t}$ 项代表了电磁场动量密度的随时间的变化率（注意 $\epsilon_0 \mu_0 = 1/c^2$）。
积分形式的推导步骤
我们将上述微分关系应用于一个闭合体积 $V$ 内的电磁相互作用。
第一步：建立力与动量变化的关系
根据牛顿第二定律的连续体形式，体积 $V$ 内物质受到的电磁力等于其机械动量随时间的变化率：
$$\int_V \mathbf{f}{em} , dV = \frac{d}{dt} \int_V \mathbf{g}{mech} , dV$$
即：
$$\int_V \mathbf{f}{em} , dV = \frac{d\mathbf{P}{mech}}{dt}$$
其中 $\mathbf{P}_{mech}$ 是体积 $V$ 内的总机械动量。
第二步：代入应力张量表达式
将微分形式的 $\mathbf{f}{em}$ 代入上述积分式中：
$$\frac{d\mathbf{P}{mech}}{dt} = \int_V \left( \nabla \cdot \mathbf{T} - \epsilon_0 \mu_0 \frac{\partial \mathbf{S}}{\partial t} \right) dV$$
利用积分的线性性质，将其拆分为两部分：
$$\frac{d\mathbf{P}_{mech}}{dt} = \int_V (\nabla \cdot \mathbf{T}) , dV - \int_V \epsilon_0 \mu_0 \frac{\partial \mathbf{S}}{\partial t} , dV$$
第三步：应用散度定理
对第一项应用高斯散度定理（Divergence Theorem），将体积积分转化为穿过体积边界 $\partial V$ 的面积分：
$$\int_V (\nabla \cdot \mathbf{T}) , dV = \oint_{\partial V} \mathbf{T} \cdot d\mathbf{a}$$
其中 $d\mathbf{a}$ 是指向体积外部的单位面积矢量。
对于第二项，如果体积 $V$ 不随时间变化，我们可以将时间导数符号移到积分号外面：
$$\int_V \epsilon_0 \mu_0 \frac{\partial \mathbf{S}}{\partial t} , dV = \frac{d}{dt} \int_V \epsilon_0 \mu_0 \mathbf{S} , dV$$
由于 $\epsilon_0 \mu_0 \mathbf{S} = \epsilon_0 (\mathbf{E} \times \mathbf{B}) = \mathbf{g}{em}$，该项即为场动量的变化率：
$$\int_V \epsilon_0 \mu_0 \frac{\partial \mathbf{S}}{\partial t} , dV = \frac{d\mathbf{P}{em}}{dt}$$
第四步：整理得到最终形式
将上述结果合并，得到：
$$\frac{d\mathbf{P}{mech}}{dt} = \oint{\partial V} \mathbf{T} \cdot d\mathbf{a} - \frac{d\mathbf{P}_{em}}{dt}$$
移项后，我们得到动量守恒的积分形式：
$$\frac{d}{dt} (\mathbf{P}{mech} + \mathbf{P}{em}) = \oint_{\partial V} \mathbf{T} \cdot d\mathbf{a}$$
物理意义解析
最终得到的方程 $\frac{d\mathbf{P}{total}}{dt} = \oint{\partial V} \mathbf{T} \cdot d\mathbf{a}$ 揭示了电磁系统中动量守恒的本质：
- 总动量守恒：方程左侧是系统总动量（机械动量 + 场动量）的变化率。
- 动量通量：方程右侧的面积分表示通过封闭曲面 $\partial V$ 的动量通量。如果系统是一个孤立系统（即边界处没有电磁应力传递，$\mathbf{T} \cdot d\mathbf{a} = 0$），那么总动量 $\mathbf{P}_{total}$ 将保持不变。
- 动量交换机制：该方程说明，机械动量的变化不仅取决于电磁场对物质的直接作用，还必须考虑场本身动量的变化。这种“场与物质交换动量”的过程，是理解辐射压力、电磁推进等现象的基础。
示例应用：辐射压力
考虑一束强激光照射在平面镜上。激光携带了显著的电磁动量 $\mathbf{P}{em}$。当激光被反射时，其动量发生改变 $\Delta \mathbf{P}{em}$。根据上述积分形式，为了满足动量守恒，镜子必须获得一个反向的机械动量变化 $\Delta \mathbf{P}_{mech}$。通过计算麦克斯韦应力张量在镜面边界上的积分，我们可以精确地计算出光子对物质产生的压力，即辐射压力。