# Field of a polarized object

Required math: calculus

Required physics: electrostatics

Reference: Griffiths, David J. (2007) Introduction to Electrodynamics, 3rd Edition; Prentice Hall – Section 4.2.1; Problem 4.10.

$\displaystyle V(\mathbf{r})=\frac{1}{4\pi\epsilon_{0}r^{2}}\hat{\mathbf{r}}\cdot\mathbf{p} \ \ \ \ \ (1)$

where ${\mathbf{r}}$ is the vector from the dipole to the observation point.

If we now consider an object consisting of polarized dielectric, we can define the polarization density, or polarization per unit volume as ${\mathbf{P}}$. The potential due to such an object is then

$\displaystyle V(\mathbf{r})=\frac{1}{4\pi\epsilon_{0}}\int\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|^{2}}\frac{\left(\mathbf{r}-\mathbf{r}'\right)\cdot\mathbf{P}(\mathbf{r}')}{\left|\mathbf{r}-\mathbf{r}'\right|}d^{3}\mathbf{r}' \ \ \ \ \ (2)$

Note that this is a volume integral over the primed coordinates ${\mathbf{r}'}$, that is, over the location of the volume element containing the polarized material.

We can transform this integral by doing a calculation of a gradient:

 $\displaystyle \nabla'\left(\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|}\right)$ $\displaystyle =$ $\displaystyle \frac{\partial}{\partial x'}\left[\frac{1}{\left[\left(x-x'\right)^{2}+\left(y-y'\right)^{2}+\left(z-z'\right)^{2}\right]^{1/2}}\right]\hat{\mathbf{x}}+\ldots\ \ \ \ \ (3)$ $\displaystyle$ $\displaystyle =$ $\displaystyle \frac{\left(x-x'\right)}{\left[\left(x-x'\right)^{2}+\left(y-y'\right)^{2}+\left(z-z'\right)^{2}\right]^{3/2}}\hat{\mathbf{x}}+\ldots\ \ \ \ \ (4)$ $\displaystyle$ $\displaystyle =$ $\displaystyle \frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|^{2}}\frac{\left(\mathbf{r}-\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|} \ \ \ \ \ (5)$

where we’ve omitted the terms in ${\hat{\mathbf{y}}}$ and ${\hat{\mathbf{z}}}$ in the first two lines since they have the same form. We can therefore write the original integral as

$\displaystyle V(\mathbf{r})=\frac{1}{4\pi\epsilon_{0}}\int\mathbf{P}\cdot\nabla'\left(\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|}\right)d^{3}\mathbf{r'} \ \ \ \ \ (6)$

A standard theorem from vector calculus says, for a function ${f}$ and vector field ${\mathbf{A}}$:

$\displaystyle \nabla\left(f\mathbf{A}\right)=\mathbf{A}\cdot\nabla f+f\nabla\cdot\mathbf{A} \ \ \ \ \ (7)$

so we can transform the integral to get

$\displaystyle V(\mathbf{r})=\frac{1}{4\pi\epsilon_{0}}\int\nabla'\cdot\left(\frac{\mathbf{P}}{\left|\mathbf{r}-\mathbf{r}'\right|}\right)d^{3}\mathbf{r'}-\frac{1}{4\pi\epsilon_{0}}\int\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|}\nabla'\cdot\mathbf{P}d^{3}\mathbf{r}' \ \ \ \ \ (8)$

The first integral is a volume integral of a divergence so we can apply the divergence theorem to transform this to a surface integral, so we finally get

$\displaystyle V(\mathbf{r})=\frac{1}{4\pi\epsilon_{0}}\int\left(\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|}\right)\mathbf{P}\cdot d\mathbf{a}'-\frac{1}{4\pi\epsilon_{0}}\int\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|}\nabla'\cdot\mathbf{P}d^{3}\mathbf{r}' \ \ \ \ \ (9)$

where the first integral is over the surface of the polarized object, and the second integral is over its volume. That is, the potential of a polarized object can be expressed as the sum of the potential of a surface charge density and a volume charge density, where we have

 $\displaystyle \sigma_{b}$ $\displaystyle \equiv$ $\displaystyle \mathbf{P}\cdot\hat{\mathbf{n}}\ \ \ \ \ (10)$ $\displaystyle \rho_{b}$ $\displaystyle \equiv$ $\displaystyle -\nabla'\cdot\mathbf{P} \ \ \ \ \ (11)$

These charge distributions are known as bound charges, which is why we’ve used a subscript ${b}$. If we’re interested in the electric field of a dipole distribution, we can work out these integrals and then take the negative gradient to get ${\mathbf{E}=-\nabla V}$, or, if the problem has the right symmetry, we can use Gauss’s law to work out the fields of the two charge distributions and then add them together.

As an example, suppose we have a sphere of radius ${R}$ with a polarization density given by

$\displaystyle \mathbf{P}(\mathbf{r}')=k\mathbf{r}' \ \ \ \ \ (12)$

where ${k}$ is a constant. Then

 $\displaystyle \sigma_{b}$ $\displaystyle =$ $\displaystyle kR\ \ \ \ \ (13)$ $\displaystyle \rho_{b}$ $\displaystyle =$ $\displaystyle -k\nabla'\cdot\mathbf{r}\ \ \ \ \ (14)$ $\displaystyle$ $\displaystyle =$ $\displaystyle -3k \ \ \ \ \ (15)$

where the last line can be found by expressing ${\mathbf{r}}$ in rectangular coordinates. Note that the total charge in the shell ${q_{s}=4\pi R^{2}kR}$ is equal and opposite to the total volume charge ${\frac{4}{3}\pi R^{3}\left(-3k\right)}$.

Rather than working out the integrals above, we can use Gauss’s law to find the electric field from a spherical shell and a uniformly charged sphere. We’ve already solved this problem in examples 1 and 2 here, so we can just quote the results. Inside the sphere, the spherical shell contributes nothing, and the volume charge gives a radial field of

 $\displaystyle E$ $\displaystyle =$ $\displaystyle -\frac{kr}{\epsilon_{0}} \ \ \ \ \ (16)$

Outside the sphere, the shell contributes

$\displaystyle E_{s}=\frac{kR^{3}}{\epsilon_{0}r^{2}} \ \ \ \ \ (17)$

and the volume charge contributes

$\displaystyle E_{v}=-\frac{kR^{3}}{\epsilon_{0}r^{2}} \ \ \ \ \ (18)$

so the total field is ${E_{s}+E_{v}=0}$.

## 20 thoughts on “Field of a polarized object”

1. georgianajackal

for volume bound charge, since div of position vector is 3 so it’s directly -3k, you can also find the div of r by expressing it as xi+yj+zk and then find their derivatives independently.

1. growescience

Not integration by parts, just subtraction. The standard theorem says

$displaystyle nablaleft(fmathbf{A}right)=mathbf{A}cdotnabla f+fnablacdotmathbf{A}$
so

$displaystyle mathbf{A}cdotnabla f=nablaleft(fmathbf{A}right)-fnablacdotmathbf{A}$

We then make the substitutions ${mathbf{A}rightarrowmathbf{P}}$ and ${frightarrowfrac{1}{left|mathbf{r}-mathbf{r}'right|}}$.