$$\require{cancel}$$

13.2: Chapter 2

Problem (2.1)

Given an electric field of the form $$\mathbf{E}=100 \text{x}$$ $$\hat{\mathbf{u}}_{\text{x}} \text{V} / \text{m}$$ find the total charge contained in the following volumes:

1) A cubical volume 1 cm on a side centered on the origin. The cube edges are parallel with the x, y, and z axes.

2) A cylindrical volume having a radius of 1 cm and a height of 2 cm centered at the origin. The axis of the cylinder is parallel with the z-axis.

Ex = 100x

Ey = 0

Ez = 0

$$\therefore \operatorname{div} \mathbf{E}=100=\rho_{\text{f}} / \varepsilon_{\text{o}}$$

This electric field corresponds to a uniform charge distribution $$\rho_{f}=100 \varepsilon_{0}$$ Coulombs/m3

∴ The total charge in

(1) The cube

\begin{aligned} Q=\left(100 \varepsilon_{0}\right)\left(10^{-6}\right)=10^{-4} \varepsilon_{0} & \\ =\frac{10^{-13}}{36 \pi}=8.84 \times 10^{-16} & \text {Coulombs. } \end{aligned}

(2) The cylinder

\begin{aligned} Q&=\left(100 \varepsilon_{0}\right)\left(2 \pi \times 10^{-6}\right)=2 \pi \varepsilon_{0} \times 10^{-4} \\ &=\frac{10^{-13}}{18}=5.56 \times 10^{-15} \text {Coulombs } \end{aligned}

since $$\varepsilon_{0}=\frac{1}{\mu_{0} C^{2}}=\frac{10^{-9}}{36 \pi}$$

Note that the above charge distribution though uniform must have planar symmetry (because Ey = Ez = 0).

Problem (2.2)

A free charge distribution is given by ρf = ar Coulombs/m3 for 0 ≤ r ≤ R and ρf = 0 for r > R. (The electric polarization P is everywhere zero).

a) Calculate the components of the electric field in and around this charge distribution. The problem has spherical symmetry so one can use Gauss' theorem (the divergence theorem).

b) Calculate the potential function corresponding to the electric field of part (a). Choose the arbitrary constants so that (1) V → 0 as r → ∞.

(2) V is continuous at r= R.

In this way show that the potential at r = 0 is given by $$V(0)=\frac{a R^{3}}{3 \varepsilon_{0}}$$ Volts.

(a) $$\text{div} \mathbf{E}=\frac{\rho_{\text{f}}}{\varepsilon_{\text{o}}}$$ since P ≡ 0

and therefore div P = 0.

So $$\int_{\text{V}} \text{div} \mathbf{E} \text{d} \tau=\int_{\text{s}} \mathbf{E} \cdot \mathbf{ds}=\frac{Q(\text{r})}{\varepsilon_{0}}$$

where Q(r) is the charge contained within a sphere of radius r.

$$Q(r)=a \int_{0}^{r}(r) 4 \pi r^{2} d r=\pi a r^{4}$$ Coulombs for r ≤ R.

But E has only a radial component by symmetry. Therefore for a spherical surface of radius r , $$\int_{S} \mathbf{E} \cdot \mathbf{d} \mathbf{s}=4 \pi r^{2} E_{r}$$

So for r≤R $$4 \pi r^{2} E_{r}=\frac{\pi a r^{4}}{\varepsilon_{0}}$$

or $$\text{E}_{\text{r}}=\frac{\text{ar}^{2}}{4 \varepsilon_{\text{O}}} \text{Volts} / \text{m}$$.

for r>R the charge is independent of r: Q = $$\pi$$aR4

$\therefore \quad E_{r}=\frac{\pi a R^{4}}{4 \pi \varepsilon_{O} r^{2}} \text { Volts/m (Coulomb's Law) }\nonumber.$

(b) Since E has only a radial component, the potential function will depend only upon r:

$E_{r}=-\left(\frac{\partial V}{\partial r}\right)\nonumber$

∴ in the region r ≤ R $$V=\frac{-a r^{3}}{12 \varepsilon_{0}}+V_{0}$$

in the region r > R $$V=\frac{\pi a R^{4}}{4 \pi \varepsilon_{0} r}$$

(The constant is zero so that V → 0 as r → ∞).

At r = R we require V to be continuous. Therefore

$V_{0}-\frac{a R^{3}}{12 \varepsilon_{0}}=\frac{a R^{3}}{4 \varepsilon_{0}}=\frac{3 a R^{3}}{12 \varepsilon_{0}}\nonumber$

So $$V_{0}=\frac{4 a R^{3}}{12 \varepsilon_{0}}=a R^{3} / 3 \varepsilon_{0}$$

The potential at the center of the charge distribution is therefore (aR3/3 εo) Volts

Problem (2.3)

A cube of side length L m is centered on the origin and its edges are parallel with the x, y, and z axes. The electric dipole vector per unit volume, P, is given by $$\mathbf{P}=\text{P}_{0} \left(\text{x} \hat{\mathbf{u}}_{\text{x}}+\text{y} \hat{\mathbf{u}}_{\text{y}}+\text{z} \hat{\mathbf{u}}_{\text{z}}\right) \ \text { Coulombs } / \text{m}^{2}$$

a) Calculate the bound charge density ρb = - div P.

b) Calculate the surface bound charge density on each face of the cube.

c) Show that the total bound charge on the cube is zero.

(a) $$\mathbf{P}=\text { Po }\left(\text{x} \hat{\mathbf{u}}_{\mathbf{x}}+\text{y} \hat{\mathbf{u}}_{\text{y}}+\text{z} \hat{\mathbf{u}}_{\text{z}}\right)$$ inside cube

P ≡ 0 Outside cube

∴ ρb = - div P = - 3Po inside the cube

= 0 outside the cube.

The total bound charge inside the cube is therefore Qv = - 3PoL3 Coulombs

(b) The discontinuity in the normal component of P gives an effective surface charge density on each face of the cube.

For example, on the face f1 there is a discontinuity in Px which is illustrated in the sketch below.

∴ The surface bound charge density is $$\sigma_{\text{b}}=+\frac{\text{P}_{\text{O}} \text{L}}{2}$$ Coulombs/m2 The total surface charge on f1 is Qs = σbL2 = PoL3/2 Coulombs. There is a similar charge on each of the other faces. Therefore the total surface charge = 6Qs = 3PoL3. The same as the volume charge.

Problem (2.4)

A disc of charge whose diameter is R meters is centered on the origin with its plane normal to the z-axis as shown in the sketch.

(a) Calculate the potential function V(z) on the axis of the disc. Sketch V(z).

(b) Make a sketch of Ez(z). Show that as z → 0 $$\mathbf{E}_{\text{z}}=\frac{\sigma}{2 \varepsilon_{0}}$$ for z > 0 and $$E_{z}=-\frac{\sigma}{2 \varepsilon_{0}}$$ for z < 0.

(This problem is not as trivial as it looks. Remember V(z) must be an even function of z. It must also go to zero as | z | → ∞, and it must be continuous at z = 0. The answer is $$V(z)=\frac{\sigma}{2 \varepsilon_{0}}[\sqrt{R^{2}+z^{2}}-|z|]$$.

$V(z)=\frac{1}{4 \pi \varepsilon_{0}} \quad \int_{0}^{R} \frac{2 \pi r d r \sigma}{\sqrt{r^{2}+z^{2}}}=\frac{\sigma}{4 \varepsilon_{0}} \int_{z^{2}}^{R+z^{2}} \frac{d u}{\sqrt{u}}\nonumber$

where u = r2 + z2 du = 2r dr

$\therefore V(z)=\frac{\sigma}{2 \varepsilon_{0}}[\sqrt{R^{2}+z^{2}}-\sqrt{z^{2}}]\nonumber$

(There is a temptation to write $$\sqrt{z^{2}}=z$$ but this would be wrong because one must use the +'ve root of z2 even when z is negative. Hence $$\sqrt{z^{2}}=|z|$$.)

For z > 0 but z small $$V \rightarrow-\frac{\sigma_{z}}{2 \varepsilon_{0}}+\frac{\sigma_{R}}{2 \varepsilon_{0}}$$

For z < 0 but z small $$V \rightarrow \frac{\sigma_{z}}{2 \varepsilon_{0}}+\frac{\sigma R}{2 \varepsilon_{0}}$$

Therefore for z > 0 $$\lim _{z \rightarrow 0}-\frac{\partial V}{\partial z}=E_{z}=\frac{\sigma}{2 \varepsilon_{0}}$$.

for z < 0 $$\lim _{z \rightarrow 0}-\frac{\partial V}{\partial z}=E_{z}=-\frac{\sigma}{2 \varepsilon_{0}}$$.

for z > 0 but |z| >> R $$V \rightarrow \frac{\pi \sigma R^{2}}{4 \pi \varepsilon_{0}}\left(\frac{1}{z}\right)$$,

for z < 0 but |z| >> R $$\text{V} \rightarrow-\frac{\pi \sigma \text{R}^{2}}{4 \pi \varepsilon_{0}}\left(\frac{1}{\text{z}}\right)$$,

ie. the potential looks like a point charge q = $$\pi$$R2σ Coulombs.

(b) $$E_{z}=-\frac{\partial V}{\partial z}$$

for z > 0 $$E_{z}=\frac{\sigma}{2 \varepsilon_{0}}\left(1-\frac{z}{\sqrt{R^{2}+z^{2}}}\right)$$

∴ @ z = 0 $$E_{z}=+\frac{\sigma}{2 \varepsilon_{0}}$$.

for z < 0 $$\text{E}_{\text{Z}}=\frac{\sigma}{2 \varepsilon_{0}}\left(-1+\frac{\text{z}}{\sqrt{\text{R}^{2}+\text{z}^{2}}}\right)$$

∴ @ z = 0 $$E_{z}=-\frac{\sigma}{2 \varepsilon_{0}}$$

Problem (2.5)

Two thin infinite plane metal plates are parallel and separated by a gap d meters as shown in the sketch. Plate #1 carries a surface charge density of +σ Coulombs/m2. Plate #2 carries a surface charge of -σ Coulombs/m2. In the metal E = 0, otherwise the charges in the metal would move and one would not have an electrostatic problem. Let the direction normal to the plates be the z direction.

(a) Use Gauss' theorem to calculate the electric field strength in the gap between the plates. Let this value be Eo.

(b) What is the value of Dz between the plates?

(c) What is the potential difference between the two metal plates?

(d) Suppose that a slab of matter whose thickness was (d/2) meters was slipped between the two metal plates. Suppose further that this slab were polarized such that Pz = Po. What would now be the potential difference between the two plates?

(e) Show that Dz is continuous across the faces of the polarized slab.

(a) Since we have two infinite sheets of charge the electric field is uniform and parallel with z (normal to the plates). Use a pill box which penetrates the metal surface on the left

$\int_{S}(\mathbf{E} \cdot \mathbf{n}) \text{d} S=Q / \varepsilon_{0}\) \quad \quad \quad $$\therefore \text{E}_{\text{o}} \text{A}=\frac{\sigma_{\text{A}}}{\varepsilon_{\text{o}}}\nonumber$ or Eo = σ/εo Volts/meter. (b) Dz = ε0Ez + Pz Here Pz = 0 ∴ Dz = ε0Eo = σ Coulombs/m2. Notice that for σ = 1 Coulomb/m2 the electric field would be 1.13 x 1011 Volts/meter. This is huge: air breaks down in a field of ~ 3 x 106 Volts/meter. Therefore 1 Coulomb/m2 is a huge charge density. (c) ∆V = Eod = σd/ε0 Volts. (d) Outside the slab Ez = 0 Inside the slab \(E_{z}=-\frac{P_{0}}{\varepsilon_{0}}$$

When placed between the two metal plates the field distributions add. Therefore in the gap Ez = Eo but in the slab $$\text{E}_{\text{z}}=\text{E}_{\text{O}}-\frac{\text{P}_{\text{O}}}{\varepsilon_{0}}$$.

The total potential drop between the plates will be

$\Delta \text{V}=\text{E}_{\text{O}}\left(\frac{\text{d}}{2}\right)+\left[\text{E}_{\text{O}}-\frac{\text{P}_{\text{O}}}{\varepsilon_{0}}\right]\left(\frac{\text{d}}{2}\right)\nonumber$

$\therefore \Delta \text{V}=\text{E}_{\text{O}} \text{d}-\frac{\text{P}_{\text{o}} \text{d}}{2 \varepsilon_{0}}=\frac{\sigma \text{d}}{\varepsilon_{0}}-\frac{\text{P}_{\text{O}} \text{d}}{2 \varepsilon_{0}}\nonumber$

The potential drop is decreased by the presence of the slab.

(e) In the gap between the slab and the plates Dz = ε0Eo = σ Coulombs/m2.

In the slab

\begin{aligned} &\begin{array}{l} D_{z}=\varepsilon_{0} E_{z}+P_{0} \\ \quad \ =\varepsilon_{0} E_{0}-P_{0}+P_{0} \end{array}\\ &\quad \ =\varepsilon_{0} \text{E}_{0}=\sigma \text { Coulombs } / \text{m}^{2} \end{aligned}\nonumber

Dz is continuous across the slab faces!

Problem (2.6)

An ellipsoid of revolution has the shape of a cigar with its axis oriented along z. The length of the cigar is 1 cm and its diameter is $$\frac{1}{2}$$ cm. The cigar is uniformly polarized: The polarization is given by

$\mathbf{P}=\text{P}_{\text{X}} \hat{\mathbf{u}}_{\text{X}}+\text{P}_{\text{Y}} \hat{\mathbf{u}}_{\text{Y}}+\text{P}_{\text{Z}} \hat{\mathbf{u}}_{\text{Z}}\nonumber$

where

Px = 0.1 Coulombs/m2

Py = 0.2 Coulombs/m2

Pz = 0.3 Coulombs/m2.

Calculate the electric field components in the ellipsoid. (They turn out to be huge ~ 1010 V/m. Air breaks down in a field of ~ 106 V/m).

For the cigar whose length is 2d and whose diameter is 2R the depolarizing coefficient is given by ( where $$\frac{R}{d}<1$$)

$$\text{N}_{\text{z}}=\left(\frac{1-\varepsilon^{2}}{\varepsilon^{3}}\right)\left(\frac{1}{2} \ln \left(\frac{1+\varepsilon}{1-\varepsilon}\right)-\varepsilon\right) \quad \text { where } \varepsilon=\sqrt{1-\left(\frac{\text{R}}{\text{d}}\right)^{2}}$$.

For this problem the ratio $$\frac{R}{d}=\frac{1}{2}$$ and therefore $$\varepsilon=\sqrt{3 / 4}$$.

according to my calculations, Nz = 0.1736.

But the sum rule states that Nx + Ny + Nz = 1 and therefore the sum Nx + Ny = 0.826. By symmetry Nx = Ny, therefore Nx = Ny = 0.413.

$$\mathbf{E}=-\frac{\text{N}_{\text{X}} \text{P}_{\text{X}}}{\varepsilon_{0}} \hat{\mathbf{u}}_{\text{x}}-\frac{\text{N}_{\text{Y}} \text{P}_{\text{Y}}}{\varepsilon_{0}} \hat{\mathbf{u}}_{\text{Y}}-\frac{\text{N}_{\text{Z}} \text{P}_{\text{z}}}{\varepsilon_{0}} \hat{\mathbf{u}}_{\text{z}}, \text { where } \varepsilon_{0}=8.85 \times 10^{-12} \text{MKS}$$,

therefore Ex = -0.47x1010 Volts/meter,

Ey = -0.94x1010 Volts/meter

Ez = -0.59x1010 Volts/meter.

Problem (2.7)

An uncharged uniformly polarized disc of radius R meters and thickness 2D meters is shown in the figure. The polarization, P0 Coulombs/m2, is directed along the axis of the disc.

(a) Calculate the effective bound charge density, ρb= - div P, everywhere.

(b) Use the bound charge density of part (a) to calculate the potential function along the axis of the disc.

(c) Calculate the electric field along the axis of the disc. Check your answer by looking at three limits: (1) the limit (D/R)<<1; (2) the limit z>0 and (z/R)>>1; and (3) the limit z<0 and (|z|/R)>>1.

(d) Calculate the displacement vector D= ε0E + P for all points along the axis of the disc.

(a) P= 0 everywhere outside the disc and therefore ρb= - divP = 0 everywhere outside the disc. Everywhere inside the disc P is constant and so its divergence is zero; ρb= 0 inside the disc. The discontinuity in the normal component of P on the surfaces at z=-D and at z=+D produces surface bound charge densities. The surface charge density carried by the surface at z= -D is σb= -P0 Coulombs/m2; the surface charge density carried by the surface at z= +D is σb= +P0 Coulombs/m2.

(b) In order to calculate the potential function along the axis of the disc that is generated by the two surface charge distributions, it is useful to begin by considering just one plane surface charge distribution; see the sketch below.

$\text{d} \text{V}_{\text{p}}=\frac{1}{4 \pi \varepsilon_{0}} \frac{\sigma 2 \pi \text{rd} \text{r}}{\sqrt{\text{r}^{2}+\text{z}^{2}}}\nonumber,$

and

$V_{p}=\frac{\sigma}{4 \varepsilon_{0}} \int_{0}^{R} \frac{2 r d r}{\sqrt{r^{2}+z^{2}}}=\frac{\sigma}{2 \varepsilon_{0}}(\sqrt{R^{2}+z^{2}}-\sqrt{z^{2}})\nonumber.$

Notice that the potential function must be symmetric in z. There is a temptation to replace $$\sqrt{z^{2}}$$ in Equation (1) by z but that would be quite wrong because z is an odd function. One must replace $$\sqrt{z^{2}}$$ by |z|.

For z>0, but (z/R)<<1, $$V \rightarrow \frac{\sigma_{R}}{2 \varepsilon_{0}}-\frac{\sigma_{z}}{2 \varepsilon_{0}}$$.

For z<0, but (|z|/R)<<1, $$V \rightarrow \frac{\sigma_{R}}{2 \varepsilon_{0}}+\frac{\sigma_{z}}{2 \varepsilon_{0}}$$.

Therefore near the charged disc the electric field has the value Ez= +σ/2ε0 on the right, and Ez= -σ/2ε0 on the left; this is the expected result based upon an analysis of an infinite uniformly charged plane. Far from the charged disc, $$\frac{R}{|z|}<<1$$, one finds

For z>0, $$V \rightarrow \frac{\pi R^{2} \sigma}{4 \pi \varepsilon_{0}}\left(\frac{1}{z}\right)$$.

For z<0, $$V \rightarrow-\frac{\pi R^{2} \sigma}{4 \pi \varepsilon_{0}}\left(\frac{1}{z}\right)$$.

From a great distance the disc of charge looks like a point charge, where Q= $$\pi$$R2σ Coulombs.

Returning to the problem of the polarized disc, the potential function along the axis of the disc can be written by inspection using Equation (1).

For z≥ D:

$V(z)=\frac{P_{0}}{2 \varepsilon_{0}} \quad(\sqrt{R^{2}+(z-D)^{2}}-\sqrt{R^{2}+(z+D)^{2}}+2 D)\nonumber.$

For -D≤ z ≤ +D:

$V(z)=\frac{P_{0}}{2 \varepsilon_{0}}(\sqrt{R^{2}+(z-D)^{2}}-\sqrt{R^{2}+(z+D)^{2}}+2 z)\nonumber.$

For z≤ -D:

$V(z)=\frac{P_{0}}{2 \varepsilon_{0}}(\sqrt{R^{2}+(z-D)^{2}}-\sqrt{R^{2}+(z+D)^{2}}-2 D)\nonumber.$

(c) The electric field along the axis is given by $$\text{E}_{\text{z}}=-\frac{\partial \text{V}}{\partial \text{z}}$$.

For z≥ +D:

$E_{z}(z)=\frac{P_{0}}{2 \varepsilon_{0}}\left(\frac{(z+D)}{\sqrt{R^{2}+(z+D)^{2}}}-\frac{(z-D)}{\sqrt{R^{2}+(z-D)^{2}}}\right)\nonumber.$

For -D ≤ z ≤ +D:

$E_{z}(z)=\frac{P_{0}}{2 \varepsilon_{0}}\left(\frac{(z+D)}{\sqrt{R^{2}+(z+D)^{2}}}-\frac{(z-D)}{\sqrt{R^{2}+(z-D)^{2}}}-2\right)\nonumber.$

For z ≤ -D:

$E_{z}(z)=\frac{P_{0}}{2 \varepsilon_{0}}\left(\frac{(z+D)}{\sqrt{R^{2}+(z+D)^{2}}}-\frac{(z-D)}{\sqrt{R^{2}+(z-D)^{2}}}\right)\nonumber.$

In the limit as (D/R)→0 the field outside the disc goes to zero like (2D/R); notice that the electric field is symmetric in z. The field inside the disc approaches the value $$\text{E}_{\text{z}}=\frac{2 \text{D}}{\text{R}}-\frac{\text{P}_{0}}{\varepsilon_{0}}$$; i.e. the field approaches the value expected for an infinite pair of oppositely charged planes in the limit D→0.

In the limit |z|/R→∞, the electric field approaches the limit $$\text{E}_{\text{z}}=\frac{1}{4 \pi \varepsilon_{0}} 4 \pi \text{R}^{2} \text{DP}_{0} \frac{1}{|\text{z}|^{3}}$$; i.e. the field due to a point dipole of moment p= 2D$$\pi$$R2P0.

(d) By definition D= ε0E+P. Outside the slab D= ε0E since P=0. Inside the slab the term P0 just cancels the constant term in Ez. The displacement vector is continuous through the surfaces of the slab. It is given for all points along the axis by $$D_{z}(z)=\frac{P_{0}}{2}\left(\frac{(z+D)}{\sqrt{R^{2}+(z+D)^{2}}}-\frac{(z-D)}{\sqrt{R^{2}+(z-D)^{2}}}\right)$$.

Problem (2.8)

An uncharged uniformly polarized disc of radius R meters and thickness 2D meters is shown in the figure. The polarization, P0 Coulombs/m2, is directed along the axis of the disc.

Calculate the potential function along the axis of the disc for z≥ D by summing the potential contributions from a collection of point dipoles. Show that

$V(z)=\frac{P_{0}}{2 \varepsilon_{0}}(\sqrt{R^{2}+(z-D)^{2}}-\sqrt{R^{2}+(z+D)^{2}}+2 D)\nonumber,$

and therefore that the field for z≥D is given by

$E_{z}(z)=\frac{P_{0}}{2 \varepsilon_{0}}\left(\frac{(z+D)}{\sqrt{R^{2}+(z+D)^{2}}}-\frac{(z-D)}{\sqrt{R^{2}+(z-D)^{2}}}\right)\nonumber.$

(Compare with the results of Problem(2.7)).

The contribution from the illustrated volume element to the potential at P can be written

$\text{d} \text{V}=\frac{1}{4 \pi \varepsilon_{0}} \frac{\text{P}_{0} \text{d} \tau \cos \theta}{\left(r^{2}+(z-\xi)^{2}\right)}\nonumber.$

The total potential at z is given by

$V(z)=\frac{1}{4 \pi \varepsilon_{0}} \int_{0}^{R} \int_{-D}^{D} \frac{P_{0}(z-\xi) d \tau}{\left(r^{2}+(z-\xi)^{2}\right)^{3 / 2}}\nonumber,$

where $$d \tau=2 \pi r d r d \xi$$. The integrations can be readily carried out. The result is

$V(z)=\frac{P_{0}}{2 \varepsilon_{0}}(\sqrt{R^{2}+(z-D)^{2}}-\sqrt{R^{2}+(z+D)^{2}}+2 D) \ \text { Volts. }\nonumber.$

$E_{z}(z)=-\frac{\partial V}{\partial z}=\frac{P_{0}}{2 \varepsilon_{0}}\left(\frac{(z+D)}{\sqrt{R^{2}+(z+D)^{2}}}-\frac{(z-D)}{\sqrt{R^{2}+(z-D)^{2}}}\right)\nonumber.$

Problem (2.9)

An uncharged uniformly polarized disc of radius R meters and thickness 2D meters is shown in the figure. The polarization, P0 Coulombs/m2, is directed along the axis of the disc.

The electric field at the center of the disc is, by direct calculation (see Problem(2.7)),

$\text{E}_{\text{Z}}(0)=-\frac{\text{P}_{0}}{\varepsilon_{0}}\left(1-\frac{\text{D}}{\sqrt{\text{R}^{2}+\text{D}^{2}}}\right)\nonumber.$

Compare this value for the electric field with that obtained using the depolarizing coefficient for an ellipsoid of revolution having the same ratio of (D/R) as the disc. Carry out the calculation for (a) $$(\text{D} / \text{R})=\frac{1}{10}$$, and for (b) $$(\text{D} / \text{R})=\frac{1}{100}$$.

The appropriate depolarizing coefficient for a pancake shaped ellipsoid is stated in the E&M notes, Figure (2.15):

$\text{N}_{\text{z}}=\frac{\text{R}^{2} \text{D}}{\left(\text{R}^{2}-\text{D}^{2}\right)^{3 / 2}}\left(\frac{\sqrt{\text{R}^{2}-\text{D}^{2}}}{\text{D}}-\operatorname{Tan}^{-1}\left(\frac{\sqrt{\text{R}^{2}-\text{D}^{2}}}{\text{D}}\right)\right) \text { for }(\text{D} / \text{R})<1\nonumber.$

If (D/R)<<1, $$\text{N}_{\text{z}} \cong 1-\frac{\pi}{2}\left(\frac{\text{D}}{\text{R}}\right)$$. Using the exact expression for Nz

(a) For R= 10D one finds Nz= 0.861, and this gives

$\text{E}_{\text{z}}(0)=-0.861\left(\frac{\text{P}_{0}}{\varepsilon_{0}}\right)\nonumber.$

By direct calculation, the exact value is $$-0.8995\left(\frac{\text{P}_{0}}{\varepsilon_{0}}\right)$$.

(b) For R= 100D one finds Nz= 0.9845, and therefore

$E_{z}(0)=-0.9845\left(\frac{P_{0}}{\varepsilon_{0}}\right)\nonumber.$

By direct calculation, the value is $$-0.9900\left(\frac{\text{P}_{0}}{\varepsilon_{0}}\right)$$.

Cylindrical discs are often approximated as ellipsoids of revolution, especially in magnetic problems, for purposes of estimating the first order correction to the field at the center of a disc having an infinite radius, $$E_{z}(0)=-\frac{P_{0}}{\varepsilon_{0}}$$.

Problem (2.10)

An uncharged sphere of radius R is polarized in such a way that the polarization vector P is radial, and its magnitude is given by $$\text{P}_{\text{r}}(\text{r})=\text{P}_{0}\left(\frac{\text{r}}{\text{R}}\right)^{2}$$.

(a) Calculate the electric field at all points inside the sphere.

(b) Calculate the electric field at all points outside the sphere.

The polarization vector possesses only a radial component, therefore $$\operatorname{div} \mathbf{P}=\frac{1}{r^{2}} \frac{\partial}{\partial r}\left(r^{2} P_{r}\right)$$. The bound charge density is $$\rho_{b}=-\operatorname{div} \mathbf{P}=-\left(\frac{4 P_{0}}{R^{2}}\right) r$$.

(a) The electric field inside the sphere can be calculated from Gauss' theorem because the field must be radial by symmetry. Thus

$4 \pi r^{2} E_{r}=\frac{1}{\varepsilon_{0}} \int_{0}^{r} \rho_{b} 4 \pi r^{2} d r=-\frac{4 \pi P_{0} r^{4}}{\varepsilon_{0} R^{2}}\nonumber$

or

$\text{E}_{\text{r}}=-\frac{\text{P}_{0}}{\varepsilon_{0}}\left(\frac{\text{r}}{\text{R}}\right)^{2}\nonumber.$

(b) There is a surface charge density on the sphere, ρs= P0 Coulombs/m2 because of the discontinuity in the normal component of the polarization vector. The total charge contained within a sphere whose radius is slightly larger than the radius R is zero. Therefore the electric field is zero everywhere outside the sphere.

Problem (2.11)

Consider an uncharged sphere having a very large radius R which is uniformly polarized along the z direction. The polarization is Pz= P0.

(a) What is the direction and strength of the electric field inside the sphere? How does this field depend upon the radius of the sphere?

(b) A tiny spherical cavity of radius b, b/R<<1, is cut out of the sphere at some point not too far from its center. The polarization in the remainder of the big sphere remains unchanged. Use the principle of superposition to calculate the electric field strength inside the small cavity of radius b.

(a) The depolarizing factors obey the sum rule Nx + Ny + Nz= 1. But for a sphere Nx=Ny=Nz, therefore each is equal to (1/3). Inside the sphere $$E_{z}=-\frac{P_{0}}{3 \varepsilon_{0}}$$. This fieldstrength does not depend at all on the radius of the sphere.

(b) The field inside the tiny sphere of polarized material which has been cut out of the big sphere is $$\text{E}_{\text{z}}=-\frac{\text{P}_{0}}{3 \varepsilon_{0}}$$. When this is added to the field in the cavity of radius b it must give a total field equal to the field strength before the tiny sphere was removed. It can therefore be concluded that the field in the cavity is zero!

Problem (2.12)

Consider an uncharged cylinder of radius R and length L. The axis of the cylinder lies along the zaxis. Let both R and L become infinitely large, but in such a way that the ratio (R/L)→0.

(a) Let the material of the cylinder be polarized along its length, i.e. Pz= P0. What is the direction and strength of the electric field inside the cylinder?

(b) Let the material of the cylinder be polarized transverse to its axis, along the x-axis say; i.e. Px= P0. What is the direction and magnitude of the electric field inside the cylinder?

(a) The depolarizing coefficient for a very long needle in the direction of its length is zero. Therefore, when the cylinder is polarized along its axis there is no electric field inside it.

(b) By symmetry the transverse depolarizing coefficients must be equal: Nx=Ny. But from the sum rule Nx+Ny= 1, since Nz=0. It follows that Nx=Ny=(1/2). The electric field inside the cylinder is given by $$E_{x}=-\frac{P_{0}}{2 \varepsilon_{0}}$$.

This problem demonstrates that the dipole field has such a long range that the electric field inside an infinitely large body depends upon its shape.

Problem (2.13)

Consider two charges q1=Q and q2=- $$\alpha$$Q. The charge q1 is located at (-b,0,0); the charge q2 is located at (b,0,0).

(a) Let $$\alpha$$=2. Show that the equipotential V=0 is a sphere of radius $$R=\frac{4 b}{3}$$ centered at $$x_{0}=-\frac{5 b}{3}$$.

(b) Let $$\alpha$$= 1/2. Show that the equipotential V=0 is again a sphere of radius $$R=\frac{4 b}{3}$$ but centered at $$x_{0}=+\frac{5 b}{3}$$.

The equipotential V=0 can be replaced by a grounded metal sphere without disturbing the potential distribution outside the sphere. This construction therefore provides the solution of the problem of a point charge brought up to a grounded conducting sphere.

$$V_{p}=\frac{1}{4 \pi \varepsilon_{0}}\left(\frac{q_{1}}{r_{1}}+\frac{q_{2}}{r_{2}}\right), \text { where } q_{1}=Q \text { and } q_{2}=-\alpha_{Q}$$

$$r_{1}=\sqrt{(x+b)^{2}+y^{2}+z^{2}}$$, and $$r_{2}=\sqrt{(x-b)^{2}+y^{2}+z^{2}}$$. Therefore V=0 when $$r_{2}=\alpha_{r_{1},} \text { or }\left(r_{2}\right)^{2}=\alpha^{2}\left(r_{1}\right)^{2}$$. A bit of algebra gives

$x^{2}+y^{2}+z^{2}+2 b x\left(\frac{\alpha^{2}+1}{\alpha^{2}-1}\right)+b^{2}=0\quad \quad \quad \quad (1) \nonumber$

By adding $$\text{b}^{2}\left(\frac{\alpha^{2}+1}{\alpha^{2}-1}\right)^{2}$$ to both sides of (1) this equation can be written

$\left(x+b\left(\frac{\alpha^{2}+1}{\alpha^{2}-1}\right)\right)^{2}+y^{2}+z^{2}=\left(\frac{2 \alpha b}{\left|\alpha^{2}-1\right|}\right)^{2}\quad \quad \quad \quad (2)\nonumber$

This is the equation of a sphere centered at $$x_{0}=-b \frac{\alpha^{2}+1}{\alpha^{2}-1}$$, and having a radius $$\text{R}=\frac{2 \alpha \text{b}}{\left|\alpha^{2}-1\right|}$$.

Problem (2.14)

Consider two charges q1=Q and q2= -Q. The charge q1 is located at (-b,0,0); the charge q2 is located at (b,0,0).

(a) Show that V=0 on the plane x=0. The region to the right of x=0 can be replaced by a conducting metal without disturbing the potential in the region x<0. This construction provides the solution of the problem of a point charge brought up to a grounded conducting plane.

(b) Show that the charge q1=Q is attracted to a grounded metal plane with a force

$\text{F}_{\text{x}}=\frac{\text{Q}^{2}}{4 \pi \varepsilon_{0}} \frac{1}{4 \text{b}^{2}} \ \text { Newtons }\nonumber.$

The electric field at x=-b is just that due to a point charge -Q located at x=+b. Therefore $$E_{x}=\frac{Q}{4 \pi \varepsilon_{0}} \frac{1}{4 b^{2}}$$, and Ey= Ez= 0. Thus the force on the charge pulling it towards the metal surface is just
$\text{F}_{\text{x}}=\frac{\text{Q}^{2}}{4 \pi \varepsilon_{0}} \frac{1}{4 \text{b}^{2}} \text { Newtons. }\nonumber.$