physics2u
Depth
PU-102 / L12 · lecture

Energy of Charge Distributions and Field Energy

~7 min read · Tier toggle changes depth (T1 plain · T2 standard · T3 rigour)

BEFORE YOU START

You should be comfortable with the definition of electric potential \(\varphi\), Gauss’s law in integral form, and basic vector calculus identities, particularly the divergence theorem and the product rule for divergence.

Why this matters

Understanding where electrostatic energy resides is fundamental to designing capacitors and analyzing field interactions.

In previous lectures, we treated the electric field as a tool to calculate forces on test charges. However, in a system of multiple charges, the field itself carries energy. This lecture establishes that the total work required to assemble a static charge distribution is not merely a property of the charges themselves, but is stored locally within the electric field. This concept of field energy density is crucial for understanding how energy is transferred in electromagnetic systems and provides the foundation for the more general energy conservation laws in electrodynamics.

The Idea of Field Energy

Imagine you have a collection of point charges scattered infinitely far apart. To bring them together to form a specific configuration, you must do work against the electrostatic forces. This work is stored as potential energy in the system. Traditionally, we might think of this energy as residing in the interactions between pairs of charges. However, a more profound view, consistent with the field concept, is that the energy is distributed throughout space wherever the electric field exists. The stronger the field, the more energy is stored per unit volume. This localizes the energy, making it easier to analyze how energy flows and is conserved in dynamic situations.

Assembling a Charge Distribution

Consider a continuous charge distribution with volume charge density \(\rho(\vec{r})\). To calculate the total energy \(U\) required to assemble this distribution from infinity, we can imagine bringing infinitesimal amounts of charge \(dq\) from infinity to their final positions. The work done to bring a small charge \(dq\) to a point where the potential is already \(\varphi\) is \(dW = \varphi \, dq\). Since \(dq = \rho \, d^3r\), the total work is the integral of \(\varphi \rho \, d^3r\). However, because the potential \(\varphi\) at any point is generated by all the other charges in the distribution, we must be careful not to double-count the interaction energy. Each pair of interacting charges is counted twice in the direct integration of \(\rho \varphi\). Therefore, the total electrostatic energy is given by half the integral of the product of charge density and potential over all space.

\[ U = \frac{1}{2} \int \rho(\vec{r}) \, \varphi(\vec{r}) \, d^3r \]

This expression is known as the charge-potential form of the energy. It is useful for calculating the energy of discrete systems or when the charge distribution is known explicitly. However, it still relies on the source charges. We can eliminate \(\rho\) using Gauss’s law in differential form, \(\nabla \cdot \vec{E} = \rho / \varepsilon_0\), to express the energy purely in terms of the field.

Energy Density in the Field

Substituting \(\rho = \varepsilon_0 \nabla \cdot \vec{E}\) into the energy expression, we get \(U = \frac{\varepsilon_0}{2} \int (\nabla \cdot \vec{E}) \varphi \, d^3r\). Using the vector identity \(\nabla \cdot (\varphi \vec{E}) = \varphi (\nabla \cdot \vec{E}) + \vec{E} \cdot (\nabla \varphi)\), we can rewrite the integrand. Since \(\vec{E} = -\nabla \varphi\), the term \(\vec{E} \cdot (\nabla \varphi)\) becomes \(-E^2\). Applying the divergence theorem to the volume integral of \(\nabla \cdot (\varphi \vec{E})\), the surface integral vanishes at infinity because the potential falls off as \(1/r\) and the field as \(1/r^2\), making the product decay faster than \(1/r^2\) over a surface area growing as \(r^2\). This leaves us with an integral solely over the square of the electric field magnitude.

\[ U = \frac{\varepsilon_0}{2} \int |\vec{E}|^2 \, d^3r \]

This result is profound: the total energy is the sum of a local quantity, the energy density \(u\), integrated over all space. The energy density is defined as \(u = \frac{1}{2} \varepsilon_0 E^2\). This confirms that the electric field itself possesses energy, independent of the specific charges that created it. This local density is always non-negative, reflecting the fact that energy is a positive quantity stored in the field configuration.

Implications for Conductors and Capacitors

Although this lecture focuses on vacuum, the concept extends directly to conductors. Inside a conductor in electrostatic equilibrium, the electric field is zero, so the energy density is zero. All the energy is stored in the field outside the conductor. For a parallel plate capacitor, the field is uniform between the plates and zero outside. The energy density is constant between the plates, \(u = \frac{1}{2} \varepsilon_0 E^2\). Multiplying this density by the volume between the plates (\(A \times d\)) gives the total energy \(U = \frac{1}{2} \varepsilon_0 E^2 A d\). Since \(E = V/d\) and \(C = \varepsilon_0 A / d\), this reduces to the familiar formula \(U = \frac{1}{2} C V^2\). This consistency check validates the field energy interpretation and shows that the energy is physically located in the space between the plates, not in the metal plates themselves.

The derivation

◆ DERIVATION EMBED · D-045 · not re-derived here

Electrostatic Energy and Field Energy Density

For a static charge distribution of density \(\rho(\vec{r})\) held in vacuum, the total work done to assemble it from infinitely dispersed charge equals \(U=\tfrac{1}{2}\int \rho\,\varphi\,d^3r\), and this same energy can be rewritten purely in terms of the electric field as \(U=\dfrac{\varepsilon_0}{2}\int |\vec{E}|^2\,d^3r\), so that the electrostatic energy is expressible as the volume integral of a strictly non-negative local field-energy density \(u=\tfrac{1}{2}\varepsilon_0 E^2\).

\[ U=\frac{1}{2}\int \rho\,\varphi\,d^3r=\frac{\varepsilon_0}{2}\int |\vec{E}|^2\,d^3r,\qquad u=\frac{1}{2}\varepsilon_0 E^2 \]
Open full derivation →
Tier 3 · watch out

The formula \(u = \tfrac12 \varepsilon_0 E^2\) is not the same as \(\tfrac12\rho\varphi\) point by point: the two integrals agree over all space, but \(\tfrac12\rho\varphi\) is zero wherever there is no charge, while \(u\) is not. Also, \(\tfrac12\varepsilon_0\int E^2\) is always positive and includes the (infinite) self-energy of point charges, whereas the sum over pairs \(\tfrac12\sum_{i\ne j} q_i q_j/(4\pi\varepsilon_0 r_{ij})\) leaves it out and can be negative. In time-dependent problems energy also flows, described by the Poynting vector \(\vec S = \vec E\times\vec B/\mu_0\) and \(\partial u/\partial t + \nabla\cdot\vec S = -\vec J\cdot\vec E\). C

Worked example

Problem. Calculate the electrostatic energy stored in the electric field of a spherical capacitor consisting of two concentric conducting spheres with radii \(a = 0.05 \ \mathrm{m}\) and \(b = 0.10 \ \mathrm{m}\), carrying a charge \(Q = 1.0 \times 10^{-9} \ \mathrm{C}\) on the inner sphere and \(-Q\) on the outer sphere. Assume the space between them is vacuum.

  1. Determine the electric field \(\vec{E}(r)\) in the region \(a < r < b\) using Gauss’s law. By symmetry, \(E(r) = \frac{Q}{4 \pi \varepsilon_0 r^2}\).
  2. Write the expression for the energy density \(u(r) = \frac{1}{2} \varepsilon_0 E(r)^2 = \frac{1}{2} \varepsilon_0 \left( \frac{Q}{4 \pi \varepsilon_0 r^2} \right)^2 = \frac{Q^2}{32 \pi^2 \varepsilon_0 r^4}\).
  3. Integrate the energy density over the volume between the spheres. The volume element is \(d^3r = 4 \pi r^2 \, dr\). The total energy is \(U = \int_a^b u(r) 4 \pi r^2 \, dr\).
  4. Substitute \(u(r)\): \(U = \int_a^b \frac{Q^2}{32 \pi^2 \varepsilon_0 r^4} 4 \pi r^2 \, dr = \frac{Q^2}{8 \pi \varepsilon_0} \int_a^b \frac{dr}{r^2}\).
  5. Evaluate the integral: \(\int_a^b r^{-2} \, dr = \left[ -\frac{1}{r} \right]_a^b = \frac{1}{a} - \frac{1}{b}\).
  6. Calculate the numerical value: \(U = \frac{Q^2}{8 \pi \varepsilon_0} \left( \frac{1}{a} - \frac{1}{b} \right)\). Using \(\varepsilon_0 \approx 8.854 \times 10^{-12} \ \mathrm{F/m}\), \(Q = 10^{-9} \ \mathrm{C}\), \(a = 0.05 \ \mathrm{m}\), \(b = 0.10 \ \mathrm{m}\). Note that \(\frac{1}{8 \pi \varepsilon_0} = \frac{k}{2}\) where \(k \approx 9 \times 10^9\). Computing directly: \(\frac{1}{8 \pi \varepsilon_0} \approx 4.494 \times 10^{9}\ \mathrm{m/F}\). The term \(\left( \frac{1}{0.05} - \frac{1}{0.10} \right) = 20 - 10 = 10 \ \mathrm{m}^{-1}\). \(Q^2 = 10^{-18} \ \mathrm{C}^2\). \(U \approx 4.494 \times 10^{9} \times 10^{-18} \times 10 = 4.49 \times 10^{-8} \ \mathrm{J}\).
\[ U \approx 4.49 \times 10^{-8} \ \mathrm{J} \]

Units check. \(\frac{\mathrm{C}^2}{\mathrm{F/m} \cdot \mathrm{m}} = \frac{\mathrm{C}^2}{\mathrm{F}} = \mathrm{J}\). The units are consistent.

CHECK YOURSELF

Three quick tests

  1. Is the electrostatic energy density always positive?
    AnswerYes, because it is proportional to \(E^2\), which is non-negative.
  2. Where is the energy stored in a parallel plate capacitor?
    AnswerIn the electric field in the space between the plates.
  3. Why does the surface integral term vanish when converting \(\int \rho \varphi\) to \(\int E^2\)?
    AnswerBecause \(\varphi\) and \(E\) decay sufficiently fast at infinity such that their product integrated over a large sphere approaches zero.
PROBLEMS

Problem set for PU-102 →

Each derivation in this unit ends with five graded problems with full solutions.

NOTES

Lecture notes drafted with a local model and checked line by line by the editor; every number in the worked example was recomputed.

Derivations cited: D-045 · Threads: force, energy, fields