高斯定律與帕松方程式：從通量到場方程的嚴謹推導

從通量到場方程：高斯定律的再認識

直覺上，高斯定律告訴我們「穿過一個閉合曲面的電通量，正比於曲面內所封住的電荷」。這個積分敘述在高中已經出現，但它真正的威力在於：它其實是一條局部的、微分形式的場方程，是 Maxwell 方程組的第一條。本文要做的，是把這個積分敘述一路推到微分形式，再導出靜電學的核心工具——帕松方程式（Poisson equation），並說明它在數學物理上的地位。

考慮靜電場 $\mathbf{E}(\mathbf{r})$，高斯定律的積分形式為

$$\oint_{\partial V} \mathbf{E}\cdot d\mathbf{A} = \frac{Q_{\text{enc}}}{\varepsilon_0} = \frac{1}{\varepsilon_0}\int_V \rho\, dV,$$

其中 $\rho(\mathbf{r})$ 是電荷密度，$\partial V$ 是體積 $V$ 的封閉邊界。這條式子之所以成立，根本上來自 Coulomb 場的 $1/r^2$ 衰減：點電荷的場線通量與曲面半徑無關，因而通量只數「裡面有幾顆電荷」，與幾何形狀無關。

從積分形式到微分形式：散度定理的關鍵一步

要把積分式變成逐點成立的局部關係，工具是 Gauss–Ostrogradsky 散度定理：

$$\oint_{\partial V} \mathbf{E}\cdot d\mathbf{A} = \int_V (\nabla\cdot\mathbf{E})\, dV.$$

把它代回高斯定律：

$$\int_V (\nabla\cdot\mathbf{E})\, dV = \int_V \frac{\rho}{\varepsilon_0}\, dV.$$

由於 $V$ 是任意選取的體積，兩邊被積函數必須處處相等，於是得到微分形式的高斯定律：

$$\boxed{\;\nabla\cdot\mathbf{E} = \frac{\rho}{\varepsilon_0}\;}$$

這一步是整個論證的樞紐：「對任意體積積分相等」推出「被積函數逐點相等」，前提是被積函數連續。電荷密度若有面電荷這類奇異分布，就要回到積分形式或用 $\delta$ 函數處理。物理上，$\nabla\cdot\mathbf{E}$ 描述的是場線的「源」——正電荷是發散源，負電荷是匯。

引入電勢：帕松方程式的誕生

靜電場無旋，$\nabla\times\mathbf{E}=\mathbf{0}$，因此可寫成純量勢的梯度：

$$\mathbf{E} = -\nabla\phi.$$

負號是慣例，使電場由高勢指向低勢。把此式代入微分形式高斯定律：

$$\nabla\cdot(-\nabla\phi) = \frac{\rho}{\varepsilon_0}\quad\Longrightarrow\quad \boxed{\;\nabla^2\phi = -\frac{\rho}{\varepsilon_0}\;}$$

這就是帕松方程式，其中 $\nabla^2 = \nabla\cdot\nabla$ 為 Laplace 算子，在直角座標下 $\nabla^2\phi = \partial_x^2\phi + \partial_y^2\phi + \partial_z^2\phi$。在無電荷區域 $\rho=0$，它退化為 Laplace 方程式 $\nabla^2\phi = 0$，其解稱為調和函數（harmonic function）。

帕松方程式把靜電學問題從「處理向量場 $\mathbf{E}$」化簡為「解一條關於純量 $\phi$ 的二階線性偏微分方程」，再配合邊界條件即可完全決定場。這正是邊界值問題（boundary value problem）的典型形式：Dirichlet 條件指定邊界上的 $\phi$，Neumann 條件指定邊界上的法向導數 $\partial\phi/\partial n$（即面電荷或法向場）。

形式解與 Green 函數

帕松方程式的解可用 Coulomb 積分直接寫出。在無界空間且場在無窮遠衰減的邊界條件下，

$$\phi(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0}\int \frac{\rho(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\, d^3r'.$$

要驗證它確實滿足帕松方程式，關鍵恆等式是

$$\nabla^2\!\left(\frac{1}{|\mathbf{r}-\mathbf{r}'|}\right) = -4\pi\,\delta^3(\mathbf{r}-\mathbf{r}').$$

換言之，$G(\mathbf{r},\mathbf{r}') = \dfrac{1}{4\pi|\mathbf{r}-\mathbf{r}'|}$ 是 Laplace 算子的 Green 函數，滿足 $\nabla^2 G = -\delta^3(\mathbf{r}-\mathbf{r}')$。把點源的響應對所有源疊加，就得到上面的積分解——這完全是線性疊加原理的體現。對有邊界的問題，則需構造滿足對應邊界條件的 Green 函數（例如鏡像電荷法本質上就是在求半空間 Dirichlet Green 函數）。

定量小範例：均勻帶電球的勢與場

考慮一個總電荷 $Q$、半徑 $R$、均勻分布的固體球，電荷密度 $\rho = \dfrac{Q}{\frac{4}{3}\pi R^3}$。利用球對稱性，球內帕松方程式在球座標下只剩徑向項：

$$\frac{1}{r^2}\frac{d}{dr}\!\left(r^2 \frac{d\phi}{dr}\right) = -\frac{\rho}{\varepsilon_0}.$$

步驟一：積分一次。兩邊乘 $r^2$ 並對 $r$ 積分：

$$r^2\frac{d\phi}{dr} = -\frac{\rho}{3\varepsilon_0}r^3 + C_1.$$

要求 $r=0$ 處場有限，故 $C_1 = 0$，得徑向場

$$E_r = -\frac{d\phi}{dr} = \frac{\rho\, r}{3\varepsilon_0} = \frac{Q\,r}{4\pi\varepsilon_0 R^3}\quad (r\le R).$$

這與直接用高斯定律對半徑 $r$ 球面所得 $E_r = \dfrac{Q_{\text{enc}}}{4\pi\varepsilon_0 r^2}$、其中 $Q_{\text{enc}} = Q (r/R)^3$ 完全一致，互為印證。

步驟二：再積分得勢。由 $d\phi/dr = -\rho r/(3\varepsilon_0)$，

$$\phi(r) = -\frac{\rho\, r^2}{6\varepsilon_0} + C_2 \quad (r\le R).$$

步驟三：套邊界條件定常數。球外為 $\phi_{\text{out}}(r) = \dfrac{Q}{4\pi\varepsilon_0 r}$。在 $r=R$ 要求 $\phi$ 連續：

$$-\frac{\rho R^2}{6\varepsilon_0} + C_2 = \frac{Q}{4\pi\varepsilon_0 R}.$$

代入 $\rho = \dfrac{3Q}{4\pi R^3}$，解得 $C_2 = \dfrac{3Q}{8\pi\varepsilon_0 R}$，最終球內電勢為

$$\phi(r) = \frac{Q}{8\pi\varepsilon_0 R}\left(3 - \frac{r^2}{R^2}\right)\quad (r\le R).$$

球心電勢 $\phi(0) = \dfrac{3Q}{8\pi\varepsilon_0 R}$，恰為球面電勢 $\dfrac{Q}{4\pi\varepsilon_0 R}$ 的 $1.5$ 倍。整個過程展示了帕松方程式的標準解法：求通解、用對稱性與有限性砍掉發散項、再以邊界連續性鎖定積分常數。

深入探討（研究所視角）

帕松方程式在更高層的理論架構中只是冰山一角。從變分原理看，靜電場滿足的方程可由場能泛函極小化導出。定義拉格朗日密度

$$\mathcal{L} = -\frac{\varepsilon_0}{2}(\nabla\phi)^2 + \rho\,\phi,$$

對應的作用量 $S=\int \mathcal{L}\, d^3r$。對 $\phi$ 做變分，Euler–Lagrange 方程 $\dfrac{\partial\mathcal{L}}{\partial\phi} - \nabla\cdot\dfrac{\partial\mathcal{L}}{\partial(\nabla\phi)} = 0$ 恰好給出 $\varepsilon_0\nabla^2\phi + \rho = 0$，即帕松方程式。這說明靜電平衡態是「場能最小、源耦合最大」的折衷，與 Thomson 定理（導體上電荷自我安排以使能量最小）同源。靜電能 $U = \frac{\varepsilon_0}{2}\int (\nabla\phi)^2 d^3r$ 正是泛函的「動能」項。

進入相對論場論，靜電的帕松方程式是協變方程的一個特例。引入四維勢 $A^\mu=(\phi/c,\mathbf{A})$ 與四維電流 $J^\mu=(c\rho,\mathbf{J})$，在 Lorenz 規範 $\partial_\mu A^\mu = 0$ 下，Maxwell 方程合併為單一波動方程

$$\Box A^\mu = \mu_0 J^\mu,\qquad \Box \equiv \frac{1}{c^2}\partial_t^2 - \nabla^2.$$

當場與源不隨時間變化、且只看純量部分時，$\partial_t\to 0$，$\Box\to -\nabla^2$，立刻回到 $\nabla^2\phi = -\rho/\varepsilon_0$。換言之，帕松方程式是達朗貝爾（d'Alembertian）波動算子在靜態極限下的退化形式；電磁波、推遲勢（retarded potential）都是它的動態推廣。

帕松方程式的結構也跨越了電磁學的疆界。Newton 重力的勢滿足 $\nabla^2\Phi = 4\pi G\rho_m$，形式完全平行，只差耦合常數與符號——這反映兩種長程作用力都源自 $1/r$ 勢、$1/r^2$ 力。在凝態與半導體物理中，Poisson–Boltzmann 方程把右側電荷密度寫成隨勢變化的 Boltzmann 分布 $\rho\propto e^{-q\phi/k_BT}$，描述電解質中的 Debye 屏蔽；其線性化版本 $\nabla^2\phi = \kappa^2\phi$（Debye–Hückel 方程）引入屏蔽長度 $\kappa^{-1}$，是電化學雙電層與半導體 PN 接面空乏區分析的基石。在密度泛函理論（DFT）的自洽場計算中，每一步都要解一次帕松方程式以更新 Hartree 勢，使它成為計算材料科學中被呼叫最頻繁的偏微分方程之一。

最後，從位勢論（potential theory） 的數學視角，Laplace 方程的解（調和函數）具有平均值性質與極值原理：調和函數在區域內不會出現嚴格的內部極大或極小，這正是 Earnshaw 定理的數學根據——純靜電力無法使帶電粒子穩定懸浮，因為穩定平衡需要勢的局部極小，而電荷所在的真空區域 $\phi$ 是調和的，不存在這種極小點。這條看似抽象的定理，深刻限制了從離子阱設計到磁懸浮的整個工程方向，也說明了為何 Paul trap 必須引入時變場才能困住離子。一條二階線性偏微分方程，竟同時牽動著古典場論、相對論、凝態物理與精密量子實驗的設計哲學。