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.

The potential of an ideal dipole is

\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. Pingback: Field of a polarized object – examples « Physics tutorials

  2. Pingback: Field of a polarized cylinder « Physics tutorials

  3. Pingback: Bound charges in a dielectric « Physics tutorials

  4. Pingback: Electric displacement « Physics tutorials

  5. Pingback: Point charge embedded in dielectric plane « Physics tutorials

  6. Pingback: Energy in a dielectric « Physics tutorials

  7. Pingback: Dielectric cube: bound charges « Physics tutorials

  8. Pingback: Point charge in dielectric sphere « Physics tutorials

  9. Pingback: Conducting sphere half-embedded in dielectric plane « Physics tutorials

  10. Pingback: Magnetization: bound currents | Physics tutorials

  11. Pingback: Electric field within a cavity inside a dielectric | Physics tutorials

  12. Pingback: Magnetic dipoles versus monopoles; an experiment | Physics tutorials

  13. Pingback: Uniform charge, polarization & magnetization | Physics tutorials

  14. Pingback: Magnetostatic formulas from electrostatic formulas | Physics tutorials

  15. Pingback: Maxwell’s equations in matter | Physics pages

Leave a Reply

Your email address will not be published. Required fields are marked *

You may use these HTML tags and attributes: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <s> <strike> <strong>