[Physics] Feynman’s calculation of surface charge density in a dielectric: is the surface polarization charge a true an electric charge

capacitancedielectricelectrostatics

This question pertains to The Feynman Lectures on Physics, Vol II 10-3

Feynman says the following about the surface polarization charge on a dielectric:

This charge can be calculated as follows. If $A$ is the area of the plate, the number of electrons that appear at the surface is the product of $A$ and $N$, the number per unit volume, and the displacement $\delta$, which we assume here is perpendicular to the surface. The total charge is obtained by multiplying by the electronic charge $q_{e}$. To get the surface density of the polarization charge induced on the surface, we divide by $A$. The magnitude of the surface charge density is

$$\sigma_{\text{pol}}=Nq_{e}\delta.$$

What this says to me is that the effective charge of one dipole moment aligned with the “line of sight” to the field point is equal to the dipole moment. It makes sense that the dipole moment would be used to calculate the charge contribution, but I don't understand why it is considered equal to the charge contribution.

In general an electric charge produces a $r^{-1}$ Coulombic potential. A dipole moment aligned with the line of sight to the field point will produce a potential proportional to $r^{-2}$. That is, the dipole potential is $$\phi\left[\mathfrak{r}\right]=\frac{1}{4\pi\epsilon_{o}}\frac{\mathfrak{p}\cdot\mathfrak{r}}{r^{3}}=\frac{1}{4\pi\epsilon_{o}}\frac{q_{e}\vec{\delta}\cdot\mathfrak{r}}{r^{3}},$$

with the angle between $\vec{\delta}$ and $\mathfrak{r}$ equal to zero. So

$$\phi\left[r\right]=\frac{1}{4\pi\epsilon_{o}}\frac{q_{e}\delta}{r^{2}}.$$

If we represent the effective charge per dipole by $q_{\text{eff}}$, then if it's really an electric charge, the potential would be

$$\phi\left[r\right]=\frac{1}{4\pi\epsilon_{o}}\frac{q_{\text{eff}}}{r}.$$

So, I ask: is the surface polarization charge really an electric charge?


Edit to add:

This part I understood when I posted my question:

This image depicts one idealized model I concocted to explain the mathematics of calculating the surface charge.
enter image description here

The cells of the dielectric (atoms, molecules, whatever) are imagined as having overlapping positive and negative rectangles of uniform charge distribution of density $\rho$. In the absence of an electric field (right) their center points are aligned. So there is no net charge at any point in the diagram. When a field is applied, the centers are displaced as shown (left). The distance between points which originally overlapped is what Feynman calls $\delta$. The positive and negative charge distributions are treated as rigid squares in the illustration.

Since the applied field and the dielectric properties of the material are assumed to be uniform there is no strictly interior point of the dielectric where the displaced charge distributions do not cancel. Within the surface layer, the negative charge distribution is displaced inward from the original surface boundary (solid black line), while the positive charge distribution is displaced outward by the same amount, $\frac{\delta}{2}.$

Integrating the electric flux over the Gaussian shell, depicted as a dashed line, results in no contribution from points strictly in the interior of the material, and no contribution from the region strictly beyond the surface layer. Since the surface layer of area $\Delta{A}$ is missing an amount of negative charge $\Delta{A}\delta\rho$, the net charge determined will be $+\Delta{A}\delta\rho$.

If the total positive movable charge per cell is $+q$, and the number of cells per volume is $N$, then $\rho=qN$. Thus the amount of exposed positive charge over the entire surface area $A$ is $ANq\delta$. Setting $q=q_e$, and dividing by the volume gives $\sigma_\text{pol}=Nq_e\delta$, as advertised by Feynman.

This is why I question whether surface charge is a true electric charge

Now consider the arrangement shown here:
dielectric wafers

The gray material is dielectric, and the black is conductive. The two parts are connected by a charged capacitor. This will induce a "surface charge" on the outer surfaces. Nowhere will there be a Coulombic potential due to this assembly.

That was the essence of my original question. The surface charge doesn't act like an isolated electric charge.

But that was an observation made because something else Feynman said caused me a problem. Here's the real reason I balked at Feynman's discussion:

This charge can be calculated as follows. If $A$ is the area of the plate, the number of electrons that appear at the surface is the product of $A$ and $N$, the number per unit volume, and the displacement $\delta$, which we assume here is perpendicular to the surface.

Referring back to the first graphic. Rather than treating the charges as being distributed over a volume, we instead consider them to be point charges arranged in the lattice, as shown. In this case any displacement will move all the positive charge points above the line connecting the black dots. In that case, there would be a Gaussian shell containing a net charge of $q_{e}N\Delta{A},$ which is huge. I'm sure that's not what Feynman had intended, but he really didn't explain very clearly what he did mean.

The problem may be that I'm pushing the idealized models past there limitations.

Best Answer

Yes, it exists as a what is known as a "bound charge". This answer might look lengthy, but do read through it carefully. I have attempted to explain it as simply as possible.

Consider a closed cubical surface within an unpolarized material. Since the material is neutral and the surface encloses no net charge, the free charge density is zero: $$Q_{f} = 0$$ Now, imagine that the material becomes polarized parallel to the cube, such that each atom can be represented as a dipole of length $d$ with center at the original center of the atom. Consider the face of the cube where the dipoles are pointing straight at the surface.

Some positive charge will have moved "out" of the surface, such that it will no longer be enclosed by the surface. Convince yourself that the positive part of the dipole $+q$ will only exit the surface if the center of the dipole is within a distance $d/2$ below the surface. Similarly, dipoles whose centers lie within a distance $d/2$ above the surface will have their negative part $-q$ pushed into the surface. In other words, we only need to consider the dipoles within a distance $d/2$ of the surface as they are the only ones that contribute to this effect. This effect is known as "bound charge", which we will denote by $Q_{b}$.

Now, consider a small area element $da$ at the surface. Let there be $n$ dipoles per unit volume, each with a dipole moment $p$ and charge $q$. The small volume element containing the affected dipoles is $$2 \left(\frac{d}{2}\right) \text{d}a$$ as dipoles both above and below by a distance $d/2$ have to be considered as mentioned above. The number of dipoles is $$(nd)\text{d}a$$

The charge is therefore $$\text{d}Q_{b} = (-qnd)\text{d}a = -np\text{d}a = -P \text{d}a$$ where $P$ is the polarization per unit volume. We have used the relations $p = qd$ and $P = np$. This is because positive charge has left the volume and negative charge has entered, so there is a net negative charge. Now, let us generalize it by writing it in vector form.

If the dipoles are inclined at an angle $\theta$ to the surface, then we just have to take the dot product, so our expression becomes: $$\text{d}Q_{b} = -\vec P\cdot \text{d}\vec a$$ If the surface is parallel to the dipoles, $\theta = 0$, so there will be no effect, as expected.

We have now derived our relationship $$Q_{b} = - \int \vec P \cdot \text{d}\vec a$$ which is a surface integral. In polarized materials, the free electric field $\vec{D}$ and the polarization field $\vec{P}$ combine to give the total electric field $\vec{E}$. Using Gauss' law: $$\int \vec{E} \cdot \text{d} \vec a = \frac{Q_{total}}{\epsilon_{0}} = \frac{Q_{f} + Q_{b}}{\epsilon_{0}} = \frac{1}{\epsilon_{0}}\left(\int \vec{D} \cdot \text{d} \vec a - \int \vec{P} \cdot \text{d} \vec a\right)$$ Since this is true for any surface, the integrands must be equal: $$\epsilon_{0} \vec{E} = \vec{D} - \vec{P}$$ This can be further expressed in terms of the electric susceptibility. Further elaboration on the potential can be found here.

Related Question