Dielectrics
Polarization, bound charges, the Clausius-Mossotti relation, and the dielectric sphere in a uniform field.
4.1 Polarization & Bound Charges
When an external electric field is applied to a dielectric material, the constituent atoms or molecules develop induced dipole moments. The macroscopic polarization $\mathbf{P}$ (dipole moment per unit volume) creates effective bound charges:
$\rho_b$ is the volume bound charge density and $\sigma_b$ is the surface bound charge density.
Derivation 1: Bound Charges from Polarization
Step 1: Potential of a polarized dielectric
A polarized volume produces a potential equivalent to a collection of dipoles. The potential at $\mathbf{r}$ is:
$$\Phi(\mathbf{r}) = \frac{1}{4\pi\epsilon_0}\int_\mathcal{V}\frac{\mathbf{P}(\mathbf{r}')\cdot(\mathbf{r}-\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|^3}\,d^3r'$$
Step 2: Use the identity $\frac{\mathbf{r}-\mathbf{r}'}{|\mathbf{r}-\mathbf{r}'|^3} = \nabla'\frac{1}{|\mathbf{r}-\mathbf{r}'|}$
$$\Phi = \frac{1}{4\pi\epsilon_0}\int\mathbf{P}\cdot\nabla'\frac{1}{|\mathbf{r}-\mathbf{r}'|}\,d^3r'$$
Step 3: Integrate by parts
Using $\nabla'\cdot(\mathbf{P}f) = f\nabla'\cdot\mathbf{P} + \mathbf{P}\cdot\nabla'f$ and the divergence theorem:
$$\Phi = \frac{1}{4\pi\epsilon_0}\left[\oint_\mathcal{S}\frac{\mathbf{P}\cdot\hat{n}'}{|\mathbf{r}-\mathbf{r}'|}\,da' - \int_\mathcal{V}\frac{\nabla'\cdot\mathbf{P}}{|\mathbf{r}-\mathbf{r}'|}\,d^3r'\right]$$
Step 4: Identify bound charges
Comparing with $\Phi = \frac{1}{4\pi\epsilon_0}\int\frac{\rho_{\text{eff}}}{|\mathbf{r}-\mathbf{r}'|}\,d^3r'$, we identify:
$$\sigma_b = \mathbf{P}\cdot\hat{n}, \qquad \rho_b = -\nabla\cdot\mathbf{P}$$
Step 5: Total bound charge is zero
$$Q_b = \oint\sigma_b\,da + \int\rho_b\,d^3r = \oint\mathbf{P}\cdot d\mathbf{a} - \int\nabla\cdot\mathbf{P}\,d^3r = 0$$
by the divergence theorem. Polarization merely redistributes charge; it cannot create net charge.
4.2 The Electric Displacement $\mathbf{D}$
Including both free and bound charges in Gauss's law: $\nabla\cdot\mathbf{E} = (\rho_f + \rho_b)/\epsilon_0 = (\rho_f - \nabla\cdot\mathbf{P})/\epsilon_0$. Define:
For linear, isotropic dielectrics: $\mathbf{P} = \epsilon_0\chi_e\mathbf{E}$ and $\mathbf{D} = \epsilon\mathbf{E}$ where $\epsilon = \epsilon_0(1 + \chi_e) = \epsilon_0\epsilon_r$.
Derivation 2: Boundary Conditions at Dielectric Interfaces
Step 1: Normal component of $\mathbf{D}$
Apply $\oint\mathbf{D}\cdot d\mathbf{a} = Q_{f,\text{enc}}$ to a thin pillbox straddling the interface:
$$D_1^{\perp} - D_2^{\perp} = \sigma_f$$
Step 2: Tangential component of $\mathbf{E}$
Apply $\oint\mathbf{E}\cdot d\boldsymbol{\ell} = 0$ around a thin rectangular loop:
$$E_1^{\parallel} = E_2^{\parallel}$$
Step 3: Refraction of field lines
At an interface with no free surface charge, $\epsilon_1 E_1^{\perp} = \epsilon_2 E_2^{\perp}$ and $E_1^{\parallel} = E_2^{\parallel}$. The field lines bend according to:
$$\frac{\tan\theta_1}{\tan\theta_2} = \frac{\epsilon_1}{\epsilon_2}$$
Step 4: Energy density in a dielectric
The energy density in a linear dielectric is:
$$u = \frac{1}{2}\mathbf{D}\cdot\mathbf{E} = \frac{1}{2}\epsilon|\mathbf{E}|^2 = \frac{|\mathbf{D}|^2}{2\epsilon}$$
Step 5: Stress tensor in a dielectric
The Maxwell stress tensor generalizes to:
$$T_{ij} = \epsilon E_i E_j - \frac{1}{2}\delta_{ij}\epsilon|\mathbf{E}|^2$$
At dielectric interfaces, the discontinuity in $T_{ij}$ produces electrostriction forces.
4.3 Clausius-Mossotti Relation
The Clausius-Mossotti relation connects the macroscopic dielectric constant $\epsilon_r$ to the microscopic atomic polarizability $\alpha$:
Derivation 3: Clausius-Mossotti Relation
Step 1: Define atomic polarizability
Each atom develops a dipole moment $\mathbf{p} = \alpha\mathbf{E}_{\text{local}}$ where $\mathbf{E}_{\text{local}}$ is the field actually experienced by the atom (not the macroscopic average field).
Step 2: Compute the local field (Lorentz cavity)
Carve a small spherical cavity around the atom. The local field has three contributions: the external field, the field from bound charges on the cavity surface, and the field from nearby atoms (zero for cubic symmetry):
$$\mathbf{E}_{\text{local}} = \mathbf{E} + \frac{\mathbf{P}}{3\epsilon_0}$$
Step 3: Relate $\mathbf{P}$ to $\mathbf{E}_{\text{local}}$
$$\mathbf{P} = N\alpha\mathbf{E}_{\text{local}} = N\alpha\left(\mathbf{E} + \frac{\mathbf{P}}{3\epsilon_0}\right)$$
Step 4: Solve for $\mathbf{P}$ in terms of $\mathbf{E}$
$$\mathbf{P} = \frac{N\alpha}{1 - N\alpha/(3\epsilon_0)}\mathbf{E} = \epsilon_0\chi_e\mathbf{E}$$
Step 5: Express in terms of $\epsilon_r$
Since $\epsilon_r = 1 + \chi_e$, rearranging gives the Clausius-Mossotti relation:
$$\frac{\epsilon_r - 1}{\epsilon_r + 2} = \frac{N\alpha}{3\epsilon_0}$$
Note the catastrophe: as $N\alpha/(3\epsilon_0) \to 1$, $\epsilon_r \to \infty$. This signals a ferroelectric transition where spontaneous polarization develops.
4.4 Dielectric Sphere in Uniform Field
Derivation 4: Dielectric Sphere in Uniform External Field
Step 1: Set up the problem
A sphere of radius $R$ and permittivity $\epsilon$ is placed in a uniform field $\mathbf{E}_0 = E_0\hat{z}$. By azimuthal symmetry, the solution involves only $P_1(\cos\theta) = \cos\theta$.
Step 2: Write general solutions inside and outside
Inside ($r < R$): $\Phi_{\text{in}} = -E_{\text{in}}r\cos\theta$ (uniform field, regular at origin).
Outside ($r > R$): $\Phi_{\text{out}} = -E_0 r\cos\theta + A\frac{\cos\theta}{r^2}$ (uniform + dipole).
Step 3: Apply boundary conditions at $r = R$
Continuity of $E_\theta$ (tangential): $E_{\text{in}} = E_0 + A/R^3$.
Continuity of $D_r$ (normal, no free charge): $\epsilon E_{\text{in}} = \epsilon_0(E_0 - 2A/R^3)$.
Step 4: Solve for $E_{\text{in}}$ and $A$
$$E_{\text{in}} = \frac{3\epsilon_0}{\epsilon + 2\epsilon_0}E_0 = \frac{3}{\epsilon_r + 2}E_0$$
$$A = \frac{\epsilon_r - 1}{\epsilon_r + 2}R^3 E_0$$
Step 5: Physical interpretation
The field inside the sphere is uniform but reduced by the factor $3/(\epsilon_r + 2)$. As $\epsilon_r \to \infty$, $E_{\text{in}} \to 0$ (conducting limit). The external field is that of the applied field plus a dipole with moment $\mathbf{p} = 4\pi\epsilon_0\frac{\epsilon_r-1}{\epsilon_r+2}R^3\mathbf{E}_0$.
4.5 Electrostatic Energy in Dielectrics
Derivation 5: Energy of a Dielectric in an External Field
Step 1: Energy to polarize a dielectric element
The work to increase the polarization by $d\mathbf{P}$ in a volume element is $dW = -\mathbf{E}\cdot d\mathbf{P}\,d^3r$.
Step 2: Integrate for a linear dielectric
For $\mathbf{P} = \epsilon_0\chi_e\mathbf{E}$, integrating from 0 to the final field:
$$W = \frac{1}{2}\int\mathbf{D}\cdot\mathbf{E}\,d^3r = \frac{1}{2}\int\epsilon|\mathbf{E}|^2\,d^3r$$
Step 3: Force on a dielectric
The force on a dielectric body in a non-uniform field can be found from the energy: $\mathbf{F} = -\nabla W$ at constant charge, or $\mathbf{F} = +\nabla W$ at constant potential.
Step 4: Dielectrophoresis
A dielectric sphere in a non-uniform field experiences a force:
$$\mathbf{F} = 2\pi\epsilon_0 R^3\frac{\epsilon_r - 1}{\epsilon_r + 2}\nabla|\mathbf{E}_0|^2$$
Step 5: Capacitance with dielectrics
For a parallel plate capacitor with dielectric: $C = \epsilon_r\epsilon_0 A/d$. The energy stored is $W = Q^2/(2C)$. Inserting a dielectric at constant charge decreases the energy (the dielectric is pulled in); at constant voltage it increases (external work is done).
Historical Context
Michael Faraday introduced the concept of dielectric materials in the 1830s, showing that insulating materials between capacitor plates increase the capacitance. The molecular theory of dielectrics was developed by Ottaviano-Fabrizio Mossotti (1850) and Rudolf Clausius (1879) independently. Peter Debye extended the theory to include permanent dipole moments in 1912, explaining the temperature dependence of the dielectric constant of polar molecules. The local field correction (Lorentz cavity) was introduced by H.A. Lorentz. Lars Onsager improved the theory in 1936 by considering the reaction field of the cavity on the molecule.
Applications
High-k Dielectrics
Modern transistor gates use hafnium oxide ($\epsilon_r \approx 25$) instead of silicon dioxide ($\epsilon_r \approx 3.9$), allowing thicker films while maintaining capacitance.
Optical Tweezers
Dielectrophoretic forces from focused laser beams trap and manipulate biological cells and nanoparticles, exploiting the gradient force on polarizable objects.
Ferroelectric Memory
Ferroelectric materials exhibit spontaneous polarization (the Clausius-Mossotti catastrophe), enabling non-volatile memory storage (FeRAM).
Metamaterials
Engineered dielectric structures can produce effective negative permittivity, enabling phenomena like perfect lensing and cloaking.
Simulation: Dielectric Response
This simulation shows the field distribution inside and outside a dielectric sphere in a uniform field, the Clausius-Mossotti relation, bound surface charge, and field reduction as a function of permittivity.
Dielectrics: Sphere in Field, Clausius-Mossotti & Bound Charges
PythonClick Run to execute the Python code
Code will be executed with Python 3 on the server
Video Lectures: Physics of Dielectrics
Electric Dipole: Dipole in External Electric Field
Electric Dipole: Potential & Field Due to a Dipole