Field of a polarized object

Required math: calculus

Required physics: electrostatics

Reference: Griffiths, David J. (2007) Introduction to Electrodynamics, 3rd Edition; Prentice Hall – Chapter 4, Post 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}}+...\ \ \ \ \ (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}}+...\ \ \ \ \ (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}.

About these ads
Post a comment or leave a trackback: Trackback URL.


  • Jerry  On Friday, 2 August 2013 at 05:54

    I don’t get part a of this problem

    • georgianajackal  On Friday, 2 August 2013 at 06:12

      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.

    • georgianajackal  On Friday, 2 August 2013 at 06:18

      for volume bound charge, since div of position vector is 3 therefore it is -3k.

  • Pelos Bar  On Sunday, 6 October 2013 at 22:32

    when you transform the integral due to the standar theorem, why is there a minus sign? integration by parts?

    • growescience  On Monday, 7 October 2013 at 15:29

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

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

      \displaystyle  \mathbf{A}\cdot\nabla f=\nabla\left(f\mathbf{A}\right)-f\nabla\cdot\mathbf{A}

      We then make the substitutions {\mathbf{A}\rightarrow\mathbf{P}} and {f\rightarrow\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|}}.


Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s


Get every new post delivered to your Inbox.

Join 511 other followers

%d bloggers like this: