3.5: Electric Field Energy in a Dielectric
( \newcommand{\kernel}{\mathrm{null}\,}\)
In Chapter 1, we have obtained two key results for the electrostatic energy: Eq. (1.55) for a charge interaction with an independent (“external”) field, and a similarly structured formula (1.60), but with an additional factor 1⁄2, for the field induced by the charges under consideration. These relations are universal, i.e. valid for dielectrics as well, provided that the charge density includes all charges (including those bound into the elementary dipoles). However, it is convenient to recast them into a form depending on the density ρ(r) of only stand-alone charges.
If a field is created only by stand-alone charges under consideration and is proportional to ρ(r) (requiring that we deal with linear dielectrics only), we can repeat all the argumentation of the beginning of Sec. 1.3, and again arrive at Eq. (1.60), provided that ϕ is now the macroscopic field’s potential. Now we can recast this result in the terms of fields – essentially as this was done in Eqs. (1.62)-(1.64), but now making a clear difference between the macroscopic electric field E=−∇ϕ and the electric displacement field D that obeys the macroscopic Maxwell equation (32). Plugging ρ(r), expressed from that equation, into Eq. (1.60), we get
U=12∫(∇⋅D)ϕd3r.
Using the fact27 that for differentiable functions ϕ and D,
(∇⋅D)ϕ=∇⋅(ϕD)−(∇ϕ)⋅D,
we may rewrite Eq. (69) as
U=12∫∇⋅(ϕD)d3r−12∫(∇ϕ)⋅Dd3r.
The divergence theorem, applied to the first term on the right-hand side, reduces it to a surface integral of ϕDn. (As a reminder, in Eq. (1.63) the integral was of ϕ(∇ϕ)n∝ϕEn.) If the surface of the volume we are considering is sufficiently far, this surface integral vanishes. On the other hand, the gradient in the second term of Eq. (71) is just (minus) field E, so that it gives
U=12∫E⋅Dd3r=12∫E(r)ε(r)E(r)d3r≡ε02∫κ(r)E2(r)d3r.Field energy in a linear dielectric
U=12∫E⋅Dd3r=12∫E(r)ε(r)E(r)d3r≡ε02∫κ(r)E2(r)d3r.
This expression is a natural generalization of Eq. (1.65), and shows that we can, as we did in free space, represent the electrostatic energy in a local form:28
U=∫u(r)d3r, with u=12E⋅D=ε2E2=D22ε.Field energy in a linear dielectric
As a sanity check, in the trivial case ε=ε0( i.e. κ=1), this result is reduced to Eq. (1.65).
Of course, Eq. (73) is valid only for linear dielectrics, because our starting point, Eq. (1.60), is only valid if ϕ is proportional to ρ. To make our calculation more general, we should intercept the calculations of Sec. 1.3 at an earlier stage, at which we have not yet used this proportionality. For example, the first of Eqs. (1.56) may be rewritten, in the continuous limit, as
δU=∫ϕ(r)δρ(r)d3r,
where the symbol δ means a small variation of the function – e.g., its change in time, sufficiently slow to ignore the relativistic and magnetic-field effects. Applying such variation to Eq. (32), and plugging the resulting relation δρ=∇⋅δD into Eq. (74), we get
deltaU=∫(∇⋅δD)ϕd3r.
(Note that in contrast to Eq. (69), this expression does not have the front factor 1⁄2.) Now repeating the same calculations as in the linear case, for the energy density variation we get a remarkably simple (and general!) expression,
δu=E⋅δD≡3∑j=1EjδDj,Energy density’s variation
where the last expression uses the Cartesian components of the vectors E and D. This is as far as we can go for the general dependence D(E). If the dependence is linear and isotropic, as in Eq. (46), then δD=εδE and
δu=εE⋅δE≡εδ(E22).
The integration of this expression over the whole variation, from the field equal to zero to a certain final distribution E(r), brings us back to Eq. (73).
An important role of Eq. (76), in its last form, is to indicate that the Cartesian coordinates of E may be interpreted as generalized forces, and those of D as generalized coordinates of the field’s effect on a unit volume of the dielectric. This allows one, in particular, to form the proper Gibbs potential energy29 of a system with an electric field E(r) fixed, at every point, by some external source:
Gibbs potential energyUG=∫VuG(r)d3r,uG(r)=u(r)−E(r)⋅D(r).
The essence of this notion is that if the generalized external force (in our case, E) is fixed, the stable equilibrium of the system corresponds to the minimum of UG, rather than of the potential energy U as such – in our case, that of the field in our system.
As the simplest illustration of this concept, let us consider a very long cylinder (with an arbitrary cross-section’s shape), made of a uniform linear dielectric, placed into a uniform external electric field, parallel to the cylinder’s axis – see Fig. 13.

For this simple problem, the equilibrium value of D inside the cylinder may be, of course, readily found without any appeal to energies. Indeed, the solution of the Laplace equation inside the cylinder, with the boundary condition (37) is evident: E(r)=Eext, and so that Eq. (46) immediately yields D(r)=εEext. However, one may wonder why does the minimum of the potential energy U, given by Eq. (73) in its last form,
UV=D22ε,
correspond to a different (zero) value of D. The Gibbs potential energy (78) immediately removes the contradiction. Indeed, for our uniform case, this energy per unit volume of the cylinder is
UGV=UV−E⋅D=D22ε−E⋅D≡3∑j=1(D2j2ε−EjDj),
and its minimum as a function of every Cartesian component of D corresponds to the correct value of the displacement: Dj=εEj, i.e. D=εE=εEext . So, the minimum of the Gibbs potential energy indeed corresponds to the systems’ equilibrium, and it may be very useful for analyses of the polarization dynamics.
Note also that Eq. (80), at this equilibrium point (only!), may be rewritten as
UGV=UV−E⋅D=D22ε−Dε⋅D≡−D22ε,
i.e. formally coincides with Eq. (79), besides the opposite sign. Another useful general relation (not limited to linear dielectrics) may be obtained by taking the variation of the uG expressed by Eq. (78), and then using Eq. (76):
δuG=δu−δ(E⋅D)=E⋅δD−(δE⋅D+E⋅δD)≡−D⋅δE.
In order to see how these expressions (with their perhaps counter-intuitive negative signs30) work, let us plug mathbfD from Eq. (33):
δuG=−(ε0E+P)⋅δE≡−δ(ε0E22)−P⋅δE.
So far, this relation is general. In the particular case when the polarization P is field-independent, we may integrate Eq. (83) over the full electric field’s variation, say from 0 to some finite value E, getting
uG=−ε0E22−P⋅E.
Again, the Gibbs energy is relevant only if E is dominated by an external field Eext, independent of the orientation of the polarization P. If, in addition, P(r)≠0 only in some finite volume V, we may integrate Eq. (84) over the volume, getting
UG=−p⋅Eext + const, with p≡∫VP(r)d3r,
where “const” means the terms independent of p. In this expression, we may readily recognize Eq. (15a) for an electric dipole p of a fixed magnitude, which was obtained in Sec. 1 in a different way.
This comparison shows again that UG is nothing extraordinary; it is just the relevant part of the potential energy of the system in a fixed external field, including the energy of its interaction with the field. Still, I would strongly recommend the reader to get a better gut feeling of the relation between the two potential energies, U and UG – for example, by using them to solve a very simple problem: calculate the force of attraction between the plates of a plane capacitor.
Reference
27 See, e.g., MA Eq. (11.4a).
28 In the Gaussian units, each of the last three expressions should be divided by 4π.
29 See, e.g., CM Sec. 1.4, in particular Eq. (1.41), and Sec. 2.1. Note that as Eq. (78) clearly illustrates, once again, that the difference between the potnetial energies UG and U, usually discussed in courses of statistical physics and/or thermodynamics as the difference between the Gibbs and Helmholtz free energies (see, e.g., SM 1.4), is more general than the effects of random thermal motion, addressed by these disciplines.
30 Some psychological relief may be provided that the fact that you may add to UG (and U) any constant – positive if you like.