2.8: Appendix 2A

Last updated
Save as PDF

Page ID: 23150

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \) \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)\(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\) \(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\)\(\newcommand{\AA}{\unicode[.8,0]{x212B}}\)

2.8 Appendix 2A.

It was pointed out in sections 2.2.1 and 2.2.2 above that the potential function, V(\(\vec R\)) generated by a distribution of electric dipoles, P(\(\vec r\)), can be calculated in two ways:

\[\text { (1) } \quad \mathrm{V}(\overrightarrow{\mathrm{R}})=\frac{1}{4 \pi \epsilon_{0}} \int \int \int_{S_{\text {pace }} d V_{\text {vol }}} \frac{(-\operatorname{div}(\overrightarrow{\mathrm{P}}))}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|}. \label{2.36}\]

This equation for the potential function is calculated from the distribution of bound charges, ρ_b = −div(\(\vec P\)). The second equation for the potential function can be written as the potential due to point dipoles \(\vec P\)dV_vol summed over the entire distribution of dipoles:

\[\text { (2) } \quad \mathrm{V}(\overrightarrow{\mathrm{R}})=\frac{1}{4 \pi \epsilon_{0}} \int \int \int_{S_{p a c e}} d V_{\text {vol }} \frac{\overrightarrow{\mathrm{P}} \cdot(\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}})}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|^{3}}. \label{2.37}\]

These two formulae, Equations (\ref{2.36} and \ref{2.37}), give the same potential function apart from a possible constant that has no effect on the resulting electric field. This statement can be proved by applying Gauss’ Theorem to the function

\[\operatorname{div}\left(\frac{\overrightarrow{\mathrm{P}}}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|}\right)=\operatorname{div}\left(\frac{\overrightarrow{\mathrm{P}}}{\sqrt{[X-x]^{2}+[Y-y]^{2}+[Z-z]^{2}}}\right). \nonumber\]

The divergence is calculated with respect to the co-ordinates of the source point, (x,y,z):

\[\operatorname{div}\left(\frac{\overrightarrow{\mathrm{P}}}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|}\right)=\frac{\partial}{\partial x}\left(\frac{\mathrm{P}_{x}}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|}\right)+\frac{\partial}{\partial y}\left(\frac{\mathrm{P}_{y}}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|}\right)+\frac{\partial}{\partial z}\left(\frac{\mathrm{P}_{z}}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|}\right).\nonumber\]

By direct differentiation one can readily show that

\[\operatorname{div}\left(\frac{\overrightarrow{\mathrm{P}}}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|}\right)=\frac{\operatorname{div}(\overrightarrow{\mathrm{P}})}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|}+\frac{\overrightarrow{\mathrm{P}} \cdot(\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}})}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|^{3}}.\nonumber\]

Remember that the differentiations are with respect to the co-ordinates of \(\vec r\), (x,y,z), and not with respect to the observer co-ordinates \(\vec R\), (X,Y,Z). Integrate the above equation over a volume, V_vol, bounded by a surface S and apply Gauss’ Theorem, section 1.3.3, to the term on the left. The result is

\[\int \int_{S} \frac{d S(\overrightarrow{\mathrm{P}} \cdot \hat{\mathbf{n}})}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|}=\int \int \int_{V_{\text {rot }}} \frac{d V_{\text {vol }} d i v(\overrightarrow{\mathrm{P}})}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|}+\int \int \int_{V_{\text {rot }}} \frac{d V_{\text {vol }} \overrightarrow{\mathrm{P}} \cdot(\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}})}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|^{3}}.\nonumber\]

Now let the volume V_vol become very large so that the surface S recedes to infinity. If the polarization distribution is limited to a finite region of space, as we shall assume, the surface integral must vanish because the polarization density on the surface, S, is zero. We are left with the identity

\[-\int \int \int_{V_{\text {vol }}} \frac{d V_{\text {vol }} d i v(\overrightarrow{\mathrm{P}})}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|}=\int \int \int_{V_{\text {red }}} \frac{d V_{\text {vol }}}{|\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}|^{3}} \cdot(\overrightarrow{\mathrm{R}}-\overrightarrow{\mathrm{r}}) \label{2.38}\]

Upon multiplying both sides of Equation (\ref{2.38}) by \(1 /\left(4 \pi \epsilon_{0}\right)\) one obtains the integral of Equation (\ref{2.36}) on the left and the integral of Equation (\ref{2.37}) on the right. It follows that the same value for the potential will be obtained, aside from a possible unimportant constant, whether one uses the formulation based upon the potential for a point charge or the formulation based upon the potential function for a point dipole.