反立方定律矢量场散度计算矛盾求解：直觉与公式结果不符

嘿，这个问题绝对是电磁学里经典的“奇点陷阱”，我来帮你理清逻辑～

首先，你的直觉完全没错：$\vec{\nabla}\cdot \vec{E} = 4\pi \delta^{3(\vec{r})$才是整个空间的正确结果。而你用球坐标散度公式得到的$-\frac{1}{r}4}$，问题出在公式的适用范围上。

1. 球坐标散度公式的隐藏前提 #

你用的球坐标散度公式：
$$\vec{\nabla}\cdot \vec{E} = \dfrac{1}{r^2}\dfrac{\partial (r^2 E_r)}{\partial r} + \dfrac{1}{r\sin\theta}\dfrac{\partial (\sin\theta E_\theta)}{\partial \theta} + \dfrac{1}{r\sin\theta}\dfrac{\partial E_\phi}{\partial \phi}$$
这个公式的推导基于场和所有偏导数在该区域连续可微的假设。但你的$\vec{E} = \dfrac{\hat{r}}{r^{3}$在$r=0$处是发散的，属于奇点，完全不满足“连续可微”的条件。所以用这个公式计算出来的$-\frac{1}{r}4}$，仅在$r≠0$的区域有效，但它无法描述原点处的奇异行为。

2. 用高斯定理推导正确结果（广义函数视角） #

要处理包含奇点的场，我们得从散度的定义本质——高斯定理出发：
$$\oint_{\partial V} \vec{E}·d\vec{S} = \int_V \vec{\nabla}\cdot \vec{E} , dV$$

• 当体积$V$不包含原点时：取任意不含原点的区域，比如两个半径分别为$a$和$b$（$a<b$）的同心球之间的体积。此时$\vec{E}$在整个区域内光滑，你用球坐标公式算的$-\frac{1}{r^4}$积分后，通量会完全抵消，左边的面积分结果为0，所以$\int_V \vec{\nabla}\cdot \vec{E} , dV=0$，即$r≠0$时$\vec{\nabla}\cdot \vec{E}=0$。
• 当体积$V$包含原点时：取一个半径为$r$的球面包围原点，计算面积分：
  $$\oint_{\partial V} \vec{E}·d\vec{S} = \oint \dfrac{\hat{r}}{r^3}·r2\hat{r}\sin\theta d\theta d\phi = \dfrac{1}{r}\int_0^{2\pi}\int_0\pi \sin\theta d\theta d\phi = 4\pi$$
  根据高斯定理，这个积分等于$\int_V \vec{\nabla}\cdot \vec{E} , dV=4\pi$。只有当$\vec{\nabla}\cdot \vec{E}$是一个在原点处取无穷大、积分后为$4\pi$的广义函数（也就是δ函数）时，才能满足这个结果，即$\vec{\nabla}\cdot \vec{E} = 4\pi \delta^3(\vec{r})$。

3. 总结一下 #

你用球坐标公式的计算本身没有错，但它只能描述$r≠0$区域的散度；而整个空间的散度需要用广义函数（δ函数）来涵盖原点的奇异行为，这时候高斯定理是更可靠的工具——毕竟散度的物理意义本来就和通量密切相关。

内容的提问来源于stack exchange，提问作者Stijn D'hondt