介质球极化电荷分布计算详解

极化电荷分布可以通过计算电场和介电常数的关系来得到。在球内部，由于电场被介质极化了，所以相当于在球内部引入了一个等效的自由电荷密度$\rho_p=- abla\cdot\vec{P}$，其中$\vec{P}$是极化矢量。根据高斯定理，球内部的极化电荷密度为：

$$\rho_p = - abla\cdot\vec{P}= -rac{1}{4\pi}\nabla\cdot\int_V\frac{\vec{P}}{r^2}\mathrm{d}V$$

其中，$V$是球的体积，$r$是距离球心的距离。根据电势能的定义，$\vec{P}$可以表示为$\vec{P}=\frac{1}{4\pi e}\nabla\cdot(\vec{E}\cdot\vec{E})$，其中$\vec{E}$是电场强度。因此，可以将上式中的$\vec{P}$用$\vec{E}$表示，得到：

$$\rho_p = -rac{1}{4\pi}\nabla\cdot\int_V\frac{1}{4\pi e}\nabla\cdot(\vec{E}\cdot\vec{E})\frac{\mathrm{d}V}{r^2}$$

根据矢量恒等式$\nabla\cdot(\vec{E}\cdot\vec{E})=\vec{E}\cdot\nabla\times\vec{E}+(\nabla\cdot\vec{E})\vec{E}$，将其代入上式，得到：

$$\rho_p = -rac{1}{4\pi}\int_V\frac{1}{4\pi e}\left[\vec{E}\cdot\nabla\times\vec{E}+(\nabla\cdot\vec{E})\vec{E}\right]\frac{\mathrm{d}V}{r^2}$$

由于球外为真空，即$\nabla\times\vec{E}=0$，所以上式中第一项为零。而在球内部，由于球对电场的影响相当于在球心引入了一个等效的自由电荷$q=\frac{4}{3}\pi R^3\rho$，其中$\rho$是球内均匀分布的电荷密度。因此，在球内部有$\nabla\cdot\vec{E}=\frac{\rho}{\varepsilon_0}$。代入上式，得到：

$$\rho_p = -rac{1}{4\pi}\int_V\frac{1}{4\pi e}\left(\frac{\rho}{\varepsilon_0}\right)\vec{E}\frac{\mathrm{d}V}{r^2}-rac{1}{4\pi}\int_{\partial V}\frac{1}{4\pi e}\vec{E}\cdot\mathrm{d}\vec{S}$$

其中，第二项是在球的表面上的积分。由于球面的法向量和电场强度的方向相同，所以上式中第二项可以化简为：

$$\rho_p = -rac{1}{4\pi}\int_V\frac{1}{4\pi e}\left(\frac{\rho}{\varepsilon_0}\right)\vec{E}\frac{\mathrm{d}V}{r^2}-rac{1}{4\pi e}\frac{1}{\varepsilon_0}\vec{E}(R)$$

其中，$\vec{E}(R)$是球面上的电场强度。由于球内部的电场强度可以用球内自由电荷密度计算得到，即$\vec{E}=\frac{1}{4\pi\varepsilon_0}\frac{q}{r^2}\hat{r}$，其中，$\hat{r}$是距离球心的单位矢量。因此，可以将上式中的$\vec{E}$代入，得到：

$$\rho_p = -rac{1}{4\pi}\int_V\frac{1}{4\pi e}\left(\frac{\rho}{\varepsilon_0}\right)\frac{1}{4\pi\varepsilon_0}\frac{q}{r^4}\hat{r}\mathrm{d}V-rac{1}{4\pi e}\frac{1}{\varepsilon_0}\frac{1}{4\pi\varepsilon_0}\frac{q}{R^2}\hat{r}$$

将$q=\frac{4}{3}\pi R^3\rho$代入，得到：

$$\rho_p = -rac{1}{4\pi}\int_V\frac{1}{4\pi e}\left(\frac{\rho}{\varepsilon_0}\right)\frac{1}{4\pi\varepsilon_0}\frac{4}{3}\pi R^3\rho\frac{1}{r^4}\hat{r}\mathrm{d}V-rac{1}{4\pi e}\frac{1}{\varepsilon_0}\frac{1}{4\pi\varepsilon_0}\frac{4}{3}\pi R^3\rho\frac{1}{R^2}\hat{r}$$

化简得到：

$$\rho_p = -rac{1}{3e}\frac{\rho^2}{\varepsilon_0}\left(\frac{R}{r}\right)^3-rac{1}{4\pi e\varepsilon_0}\frac{1}{r^2}\frac{4}{3}\pi R^3\rho$$

因此，极化电荷密度分布为：

$$\rho_p = \begin{cases}-rac{1}{3e}\frac{\rho^2}{\varepsilon_0}\left(\frac{R}{r}\right)^3 & r<R\ -rac{1}{4\pi e\varepsilon_0}\frac{1}{r^2}\frac{4}{3}\pi R^3\rho & r\geq R\end{cases}$$

其中，$r$是距离球心的距离，$\rho$是球内的电荷密度，$e$是介质的介电常数。