Skip to content

Lines That End Inside the Filling

Irodov 3.104DifficultyDifficulty 3 of 5UniquenessUniqueness 3 of 5

Problem Irodov 3.104

A parallel-plate capacitor has plates of area SS a distance dd apart. The gap is filled with an isotropic dielectric whose permittivity changes linearly across the gap, perpendicular to the plates, from ε1\varepsilon_1 at one plate to ε2\varepsilon_2 at the other, with ε2>ε1\varepsilon_2 > \varepsilon_1. Find:

  1. the capacitance;

  2. the volume density of bound charge as a function of ε\varepsilon, when the capacitor carries the charge qq and its field E\vect E points toward increasing ε\varepsilon.

The capacitor edge on. The permittivity of the filling grows linearly with the distance x from the lower plate, from ε_1 there to ε_2 at the upper plate, a distance d away.

The capacitor edge on. The permittivity of the filling grows linearly with the distance xx from the lower plate, from ε1\varepsilon_1 there to ε2\varepsilon_2 at the upper plate, a distance dd away.

Given
  • plate area SS
  • gap dd
  • ε\varepsilon linear from ε1\varepsilon_1 to ε2>ε1\varepsilon_2 > \varepsilon_1
  • charge qq
  • E\vect E toward growing ε\varepsilon
Find
(a) CC · (b) ρ′(ε)\rho'(\varepsilon)

Solution

Step 1The permittivity profile and the displacement

Measure xx across the gap from the plate where the permittivity is ε1\varepsilon_1. A linear profile with the given end values is

ε(x)=ε1+ε2−ε1d x(1) \varepsilon(x) = \varepsilon_1 + \frac{\varepsilon_2 - \varepsilon_1}{d}\,x \tag{1}

Let that plate carry qq and the other −q-q. Edge effects aside, the fields are uniform along the plates and normal to them, and Gauss’s theorem for D\vect D sees only the free charge, so DD is the same at every xx:

D=qSD = \frac{q}{S}
Step 2The field and the voltage

Dividing by ε0ε\varepsilon_0\varepsilon gives the field at each depth:

E(x)=qε0S ε(x)(2) E(x) = \frac{q}{\varepsilon_0 S\,\varepsilon(x)} \tag{2}

The voltage is its integral across the gap. With (1), change the variable to ε\varepsilon, using dx=d dε/(ε2−ε1)dx = d\,d\varepsilon/(\varepsilon_2 - \varepsilon_1):

V=∫0dE dx=qε0S∫0ddxε(x)=qε0S⋅dε2−ε1∫ε1ε2dεε=q dε0S(ε2−ε1) ln⁡ε2ε1\begin{aligned} V &= \int_0^d E\,dx = \frac{q}{\varepsilon_0 S}\int_0^d \frac{dx}{\varepsilon(x)} = \frac{q}{\varepsilon_0 S}\cdot\frac{d}{\varepsilon_2 - \varepsilon_1}\int_{\varepsilon_1}^{\varepsilon_2}\frac{d\varepsilon}{\varepsilon}\\ &= \frac{q\,d}{\varepsilon_0 S(\varepsilon_2 - \varepsilon_1)}\,\ln\frac{\varepsilon_2}{\varepsilon_1} \end{aligned}
Step 3The capacitance

The capacitance is C=q/VC = q/V, and qq cancels:

C=ε0(ε2−ε1)Sd ln⁡(ε2/ε1)C = \dfrac{\varepsilon_0(\varepsilon_2 - \varepsilon_1)S}{d\,\ln(\varepsilon_2/\varepsilon_1)}
(3)
Step 4The bound volume charge

The field points toward growing ε\varepsilon, so it runs along +x+x and the plate at x=0x = 0 is the positive one, as assumed. The polarization is P=D−ε0E\vect P = \vect D - \varepsilon_0\vect E, along +x+x; with (2):

P=qS−qS ε=qS(1−1ε)P = \frac{q}{S} - \frac{q}{S\,\varepsilon} = \frac{q}{S}\left(1 - \frac{1}{\varepsilon}\right)

The bound volume density is ρ′=−∇⋅P\rho' = -\nabla\cdot\vect P, which here is −dP/dx-dP/dx. Differentiating through ε(x)\varepsilon(x), whose slope is (ε2−ε1)/d(\varepsilon_2 - \varepsilon_1)/d by (1):

dPdx=qS⋅1ε2 dεdx=q(ε2−ε1)S d ε2\frac{dP}{dx} = \frac{q}{S}\cdot\frac{1}{\varepsilon^2}\,\frac{d\varepsilon}{dx} = \frac{q(\varepsilon_2 - \varepsilon_1)}{S\,d\,\varepsilon^2}

so that

ρ′=−q(ε2−ε1)S d ε2\rho' = -\dfrac{q(\varepsilon_2 - \varepsilon_1)}{S\,d\,\varepsilon^2}
(4)
The field for a positive lower plate, drawn for ε growing from 1 to 4. The number of lines falls as 1/ε, and each line that stops ends on a negative bound charge inside the filling; the positive bound charge on the upper face balances them.

The field for a positive lower plate, drawn for ε\varepsilon growing from 1 to 4. The number of lines falls as 1/ε1/\varepsilon, and each line that stops ends on a negative bound charge inside the filling; the positive bound charge on the upper face balances them.

The bound charge is negative throughout the filling, densest where ε\varepsilon is smallest, next to the positive plate.

Check

The bound charge must total zero. The volume charge is S∫0dρ′ dxS\int_0^d \rho'\,dx; changing the variable to ε\varepsilon as before,

S∫0dρ′ dx=−q∫ε1ε2dεε2=−q(1ε1−1ε2).S\int_0^d \rho'\,dx = -q\int_{\varepsilon_1}^{\varepsilon_2}\frac{d\varepsilon}{\varepsilon^2} = -q\left(\frac{1}{\varepsilon_1} - \frac{1}{\varepsilon_2}\right).

The surface bound charges are −P(0)S=−q(1−1/ε1)-P(0)S = -q(1 - 1/\varepsilon_1) at the positive plate and +P(d)S=q(1−1/ε2)+P(d)S = q(1 - 1/\varepsilon_2) at the negative one, and the three sum to zero. Uniform limit: as ε2→ε1\varepsilon_2 \to \varepsilon_1, ln⁡(ε2/ε1)≈(ε2−ε1)/ε1\ln(\varepsilon_2/\varepsilon_1) \approx (\varepsilon_2 - \varepsilon_1)/\varepsilon_1, so (3) tends to ε0ε1S/d\varepsilon_0\varepsilon_1 S/d, and (4) vanishes.

Remark

The capacitance is the plain parallel-plate value ε0εˉS/d\varepsilon_0\bar\varepsilon S/d with an effective permittivity εˉ=(ε2−ε1)/ln⁡(ε2/ε1)\bar\varepsilon = (\varepsilon_2 - \varepsilon_1)/\ln(\varepsilon_2/\varepsilon_1), the logarithmic mean of the two end values. A series stack averages 1/ε1/\varepsilon, not ε\varepsilon, so this mean always falls below the arithmetic mean (ε1+ε2)/2(\varepsilon_1 + \varepsilon_2)/2: the weak end of the filling costs more than the strong end gains.