Skip to content

Electrostatics

2.1 Coulomb’s Law and the Electric Field

Coulomb’s Law: The force between two point charges q1q_1 and q2q_2 separated by distance rr:

F=14πε0q1q2r2r^\mathbf{F} = \frac{1}{4\pi\varepsilon_0} \frac{q_1 q_2}{r^2} \hat{\mathbf{r}}

The electric field due to a point charge qq at position r\mathbf{r}:

E(r)=14πε0qr2r^\mathbf{E}(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0} \frac{q}{|\mathbf{r}|^2} \hat{\mathbf{r}}

Superposition Principle: The field due to a collection of charges is the vector sum of individual Fields.

Electric Field: 1/r^2 Dependence

2.2 Gauss’s Law Applications

Example: Infinite plane of charge with surface charge density σ\sigma.

Choose a Gaussian “pillbox” of area AA straddling the plane. By symmetry, E\mathbf{E} is Perpendicular to the plane. Gauss’s law:

2EA=σAε0    E=σ2ε02EA = \frac{\sigma A}{\varepsilon_0} \implies E = \frac{\sigma}{2\varepsilon_0}

The field is uniform and perpendicular to the plane, pointing away from positive charge.

Example: Uniformly charged sphere of radius RR with total charge QQ.

For r>Rr \gt R: E=Q4πε0r2r^\mathbf{E} = \frac{Q}{4\pi\varepsilon_0 r^2} \hat{\mathbf{r}} (identical to a point charge).

For r<Rr \lt R: E=Qr4πε0R3E = \frac{Qr}{4\pi\varepsilon_0 R^3} (linear in rr).

2.3 Electric Potential

The electric potential is defined by E=V\mathbf{E} = -\nabla V (for electrostatics, where ×E=0\nabla \times \mathbf{E} = \mathbf{0}).

For a point charge: V(r)=14πε0qrV(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0} \frac{q}{r} (choosing V()=0V(\infty) = 0).

Theorem 2.1. ×E=0\nabla \times \mathbf{E} = \mathbf{0} in electrostatics implies E\mathbf{E} is Conservative, so the line integral ABEdl=V(A)V(B)\int_A^B \mathbf{E} \cdot d\mathbf{l} = V(A) - V(B) is Path-independent.

2.4 Poisson’s and Laplace’s Equations

Substituting E=V\mathbf{E} = -\nabla V into Gauss’s law:

(V)=2V=ρε0\nabla \cdot (-\nabla V) = -\nabla^2 V = \frac{\rho}{\varepsilon_0}

This is Poisson’s equation:

2V=ρε0\nabla^2 V = -\frac{\rho}{\varepsilon_0}

In regions with ρ=0\rho = 0This reduces to Laplace’s equation:

2V=0\nabla^2 V = 0

Theorem 2.2 (Uniqueness --- statement). The solution to Laplace’s (or Poisson’s) equation in a Region is unique given either Dirichlet boundary conditions (VV specified on the boundary) or Neumann boundary conditions (V/n\partial V / \partial n specified on the boundary).

2.5 Worked Example

Problem. Two infinite conducting plates at x=0x = 0 and x=dx = d are held at potentials V=0V = 0 and V=V0V = V_0 respectively. Find the potential and field between them.

Solution. Between the plates, ρ=0\rho = 0So 2V=0\nabla^2 V = 0. By symmetry, VV depends only on xx:

d2Vdx2=0    V(x)=Ax+B\frac{d^2V}{dx^2} = 0 \implies V(x) = Ax + B

Boundary conditions: V(0)=0    B=0V(0) = 0 \implies B = 0. V(d)=V0    A=V0/dV(d) = V_0 \implies A = V_0/d.

V(x)=V0dx,E=dVdxx^=V0dx^V(x) = \frac{V_0}{d} x, \quad \mathbf{E} = -\frac{dV}{dx}\hat{\mathbf{x}} = -\frac{V_0}{d}\hat{\mathbf{x}}

\blacksquare

2.6 Gauss’s Law: Cylindrical Symmetry

Example: Infinite line charge with linear charge density λ\lambda.

By cylindrical symmetry, E\mathbf{E} points radially outward and depends only on rr. Choose a Gaussian cylinder of radius rr and length LL:

EdA=E2πrL=λLε0\oint \mathbf{E} \cdot d\mathbf{A} = E \cdot 2\pi r L = \frac{\lambda L}{\varepsilon_0}

E=λ2πε0rr^\mathbf{E} = \frac{\lambda}{2\pi\varepsilon_0 r}\,\hat{\mathbf{r}}

Example: Coaxial cable. An inner conductor of radius aa carries linear charge density +λ+\lambdaAnd an outer conducting shell of radius bb carries λ-\lambda.

For r<ar \lt a: E=0\mathbf{E} = \mathbf{0} (conductor interior).

For a<r<ba \lt r \lt b: E=λ2πε0rr^\mathbf{E} = \frac{\lambda}{2\pi\varepsilon_0 r}\,\hat{\mathbf{r}}.

For r>br \gt b: E=0\mathbf{E} = \mathbf{0} (total enclosed charge is zero).

The potential difference between the conductors:

V(a)V(b)=abEdl=λ2πε0ln ⁣(ba)V(a) - V(b) = -\int_a^b \mathbf{E} \cdot d\mathbf{l} = \frac{\lambda}{2\pi\varepsilon_0}\ln\!\left(\frac{b}{a}\right)

2.7 The Uniqueness Theorem

Theorem 2.3 (Uniqueness for Dirichlet conditions). The solution to Poisson’s equation 2V=ρ/ε0\nabla^2 V = -\rho/\varepsilon_0 in a volume V\mathcal{V} is unique if VV is specified on the Boundary S\mathcal{S}.

Proof. Suppose V1V_1 and V2V_2 both satisfy Poisson’s equation with the same boundary Conditions. Define U=V1V2U = V_1 - V_2. Then 2U=0\nabla^2 U = 0 in V\mathcal{V} and U=0U = 0 on S\mathcal{S}.

Apply Green’s first identity with ϕ=ψ=U\phi = \psi = U:

V(U2U+U2)dV=SUUndA\int_{\mathcal{V}} \left(U\,\nabla^2 U + \lvert\nabla U\rvert^2\right) dV = \oint_{\mathcal{S}} U\,\frac{\partial U}{\partial n}\, dA

Since 2U=0\nabla^2 U = 0 and U=0U = 0 on S\mathcal{S}:

VU2dV=0\int_{\mathcal{V}} \lvert\nabla U\rvert^2\, dV = 0

Since the integrand is non-negative, U=0\nabla U = \mathbf{0} everywhere in V\mathcal{V}So UU is Constant. With U=0U = 0 on the boundary, U=0U = 0 throughout V\mathcal{V}. Hence V1=V2V_1 = V_2. \blacksquare

Theorem 2.4 (Uniqueness for Neumann conditions). The solution is unique up to an additive Constant when V/n\partial V/\partial n is specified on S\mathcal{S}.

Proof. The same argument applies, but now U/n=0\partial U/\partial n = 0 on S\mathcal{S} and the Right-hand side of Green’s identity vanishes for a different reason. We again conclude U=0\nabla U = \mathbf{0}So UU is constant. \blacksquare

2.8 Method of Images

The method of images replaces a problem with conductors by an equivalent problem with charges only, Exploiting the uniqueness theorem.

Point charge above a grounded plane. A charge qq is placed at distance dd above an Infinite grounded conducting plane (V=0V = 0 at z=0z = 0).

Replace the plane by an image charge q=qq' = -q at z=dz = -d. The potential for z>0z \gt 0 is:

V(x,y,z)=14πε0[qx2+y2+(zd)2qx2+y2+(z+d)2]V(x,y,z) = \frac{1}{4\pi\varepsilon_0}\left[\frac{q}{\sqrt{x^2 + y^2 + (z-d)^2}} - \frac{q}{\sqrt{x^2 + y^2 + (z+d)^2}}\right]

This satisfies 2V=0\nabla^2 V = 0 for z>0z \gt 0 (away from the charge), V=0V = 0 at z=0z = 0And V0V \to 0 as rr \to \infty. By the uniqueness theorem, this is the correct solution.

The force on qq is the force due to the image charge:

F=q24πε0(2d)2z^\mathbf{F} = -\frac{q^2}{4\pi\varepsilon_0 (2d)^2}\,\hat{\mathbf{z}}

The induced surface charge density on the plane:

σ(x,y)=ε0Vzz=0=qd2π(x2+y2+d2)3/2\sigma(x,y) = -\varepsilon_0 \left.\frac{\partial V}{\partial z}\right|_{z=0} = -\frac{qd}{2\pi(x^2+y^2+d^2)^{3/2}}

Example: Point charge inside a grounded sphere. A charge qq is at distance aa from the centre Of a grounded conducting sphere of radius RR (a<Ra \lt R).

The image charge is q=qR/aq' = -qR/a located at distance b=R2/ab = R^2/a from the centre, along the same Radial line.

Solution: Verifying the image charge

We must verify that V=0V = 0 on the sphere. Place qq at distance aa from the origin along the zz-axis and qq' at distance bb along the zz-axis. At any point on the sphere at distance RR From the origin, the distances to qq and qq' are d1d_1 and d2d_2 where:

d12=R2+a22Racosθ,d22=R2+b22Rbcosθd_1^2 = R^2 + a^2 - 2Ra\cos\theta, \quad d_2^2 = R^2 + b^2 - 2Rb\cos\theta

For V=0V = 0 on the sphere, we need q/d1=q/d2q/d_1 = -q'/d_2 for all θ\theta. This requires the ratio d2/d1d_2/d_1 to be constant. Setting b=R2/ab = R^2/a:

d22d12=R2+R4/a22R3cosθ/aR2+a22Racosθ=R2a2\frac{d_2^2}{d_1^2} = \frac{R^2 + R^4/a^2 - 2R^3\cos\theta/a}{R^2 + a^2 - 2Ra\cos\theta} = \frac{R^2}{a^2}

The ratio is indeed constant. Choosing q=qR/aq' = -qR/a gives q/d1+q/d2=0q/d_1 + q'/d_2 = 0 on the sphere. \blacksquare

2.9 Multipole Expansion

For a localized charge distribution ρ(r)\rho(\mathbf{r}')The potential at large distance r=rr=rr = \lvert\mathbf{r}\rvert \gg r' = \lvert\mathbf{r}'\rvert is expanded using 1rr=n=0rnrn+1Pn(cosα)\frac{1}{\lvert\mathbf{r}-\mathbf{r}'\rvert} = \sum_{n=0}^{\infty} \frac{r'^n}{r^{n+1}} P_n(\cos\alpha) Where α\alpha is the angle between r\mathbf{r} and r\mathbf{r}':

V(r)=14πε0n=01rn+1rnPn(cosα)ρ(r)d3rV(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0}\sum_{n=0}^{\infty}\frac{1}{r^{n+1}}\int r'^n P_n(\cos\alpha)\,\rho(\mathbf{r}')\,d^3\mathbf{r}'

Monopole term (n=0n = 0):

V0=14πε0Qr,Q=ρ(r)d3rV_0 = \frac{1}{4\pi\varepsilon_0}\frac{Q}{r}, \quad Q = \int \rho(\mathbf{r}')\,d^3\mathbf{r}'

This is the potential of a point charge at the origin.

Dipole term (n=1n = 1):

V1=14πε0pr^r2,p=rρ(r)d3rV_1 = \frac{1}{4\pi\varepsilon_0}\frac{\mathbf{p} \cdot \hat{\mathbf{r}}}{r^2}, \quad \mathbf{p} = \int \mathbf{r}'\,\rho(\mathbf{r}')\,d^3\mathbf{r}'

Where p\mathbf{p} is the electric dipole moment.

Quadrupole term (n=2n = 2): Depends on the quadrupole moment tensor:

Qij=(3rirjr2δij)ρ(r)d3rQ_{ij} = \int (3r_i' r_j' - r'^2 \delta_{ij})\,\rho(\mathbf{r}')\,d^3\mathbf{r}'

V2=14πε012r3i,jQijr^ir^jV_2 = \frac{1}{4\pi\varepsilon_0}\frac{1}{2r^3}\sum_{i,j} Q_{ij}\,\hat{r}_i\,\hat{r}_j

For a neutral charge distribution (Q=0Q = 0), the dipole term dominates. If additionally p=0\mathbf{p} = \mathbf{0}The quadrupole term dominates.

Example: Dipole potential of two charges

A charge +q+q at z=+d/2z = +d/2 and q-q at z=d/2z = -d/2.

The dipole moment: p=q(d/2)z^+(q)(d/2)z^=qdz^\mathbf{p} = q(d/2)\,\hat{\mathbf{z}} + (-q)(-d/2)\,\hat{\mathbf{z}} = qd\,\hat{\mathbf{z}}.

On the zz-axis (θ=0\theta = 0): V1=qd4πε0r2V_1 = \frac{qd}{4\pi\varepsilon_0 r^2}.

In the equatorial plane (θ=π/2\theta = \pi/2): V1=0V_1 = 0.

The exact potential on the zz-axis is:

V=q4πε0(1rd/21r+d/2)=q4πε0dr2d2/4V = \frac{q}{4\pi\varepsilon_0}\left(\frac{1}{r-d/2} - \frac{1}{r+d/2}\right) = \frac{q}{4\pi\varepsilon_0}\frac{d}{r^2 - d^2/4}

For rdr \gg d: this reduces to V1=qd4πε0r2V_1 = \frac{qd}{4\pi\varepsilon_0 r^2}Confirming the Dipole approximation. \blacksquare

2.10 Dielectrics

Polarization. When an external field E\mathbf{E} is applied to a dielectric, the material Develops a polarization P\mathbf{P}The dipole moment per unit volume. This produces bound charges:

ρb=P,σb=Pn^\rho_b = -\nabla \cdot \mathbf{P}, \quad \sigma_b = \mathbf{P} \cdot \hat{\mathbf{n}}

The displacement field D\mathbf{D} is defined as:

D=ε0E+P\mathbf{D} = \varepsilon_0 \mathbf{E} + \mathbf{P}

Gauss’s law in terms of D\mathbf{D}:

D=ρf\nabla \cdot \mathbf{D} = \rho_f

Where ρf\rho_f is the free charge density. This form is useful because D\mathbf{D} depends Only on free charges, not bound charges.

Linear dielectrics. For an isotropic linear dielectric:

P=ε0χeE,D=εE\mathbf{P} = \varepsilon_0 \chi_e \mathbf{E}, \quad \mathbf{D} = \varepsilon \mathbf{E}

Where χe\chi_e is the electric susceptibility and ε=ε0(1+χe)\varepsilon = \varepsilon_0(1 + \chi_e) is the Permittivity. The relative permittivity (dielectric constant) is εr=ε/ε0=1+χe\varepsilon_r = \varepsilon/\varepsilon_0 = 1 + \chi_e.

Boundary conditions at dielectric interfaces (no free charges):

D1n=D2n    ε1E1n=ε2E2nD_{1n} = D_{2n} \implies \varepsilon_1 E_{1n} = \varepsilon_2 E_{2n}

E1t=E2tE_{1t} = E_{2t}

The tangential component of E\mathbf{E} is continuous, but the normal component changes. The angles of the field with respect to the normal satisfy ε1tanθ2=ε2tanθ1\varepsilon_1 \tan\theta_2 = \varepsilon_2 \tan\theta_1.

Example: Dielectric slab in a uniform field

A dielectric slab of permittivity ε\varepsilon and thickness dd is placed in a uniform External field E0\mathbf{E}_0 perpendicular to its faces.

Outside the slab: E=E0\mathbf{E} = \mathbf{E}_0.

Inside the slab: by continuity of DnD_n:

Din=Dout=ε0E0D_{\mathrm{in} = D_{\mathrm{out} = \varepsilon_0 E_0}}

E_{\mathrm{in} = \frac{D_{\mathrm{in}}{\varepsilon} = \frac{\varepsilon_0}{\varepsilon} E_0 = \frac{E_0}{\varepsilon_r}}}

The polarization: P=ε0χeEin=ε0(εr1)E0εrP = \varepsilon_0 \chi_e E_{\mathrm{in} = \varepsilon_0 (\varepsilon_r - 1) \frac{E_0}{\varepsilon_r}}.

The bound surface charge density on each face:

σb=±P=±ε0(11εr)E0\sigma_b = \pm P = \pm \varepsilon_0 \left(1 - \frac{1}{\varepsilon_r}\right) E_0

The bound charges produce a field opposing E0\mathbf{E}_0Reducing the net field inside the Dielectric. \blacksquare