The dipole moment of an array of charges, \mathbf p = \sum_{i=1}^N q_i \mathbf {d}_i \, , determines the degree of polarity of the array, but for a neutral array it is simply a vector property of the array with no information about the array's absolute location. The dipole moment
density of the array
p(
r) contains both the location of the array and its dipole moment. When it comes time to calculate the electric field in some region containing the array, Maxwell's equations are solved, and the information about the charge array is contained in the
polarization density P(
r) of Maxwell's equations. Depending upon how fine-grained an assessment of the electric field is required, more or less information about the charge array will have to be expressed by
P(
r). As explained below, sometimes it is sufficiently accurate to take . Sometimes a more detailed description is needed (for example, supplementing the dipole moment density with an additional quadrupole density) and sometimes even more elaborate versions of
P(
r) are necessary. It now is explored just in what way the polarization density
P(
r) that enters
Maxwell's equations is related to the dipole moment
p of an overall neutral array of charges, and also to the
dipole moment density p(
r) (which describes not only the dipole moment, but also the array location). Only static situations are considered in what follows, so
P(
r) has no time dependence, and there is no
displacement current. First is some discussion of the polarization density
P(
r). That discussion is followed with several particular examples. A formulation of
Maxwell's equations based upon division of charges and currents into "free" and "bound" charges and currents leads to introduction of the
D- and
P-fields: \mathbf{D} = \varepsilon _0 \mathbf{E} + \mathbf{P}\, , where
P is called the
polarization density. In this formulation, the divergence of this equation yields: \nabla \cdot \mathbf{D} = \rho_\text{f} = \varepsilon _0 \nabla \cdot \mathbf{E} +\nabla \cdot \mathbf{P}\, , and as the divergence term in
E is the
total charge, and
ρf is "free charge", we are left with the relation: \nabla \cdot \mathbf{P} = -\rho_\text{b} \, , with
ρb as the bound charge, by which is meant the difference between the total and the free charge densities. As an aside, in the absence of magnetic effects, Maxwell's equations specify that \nabla \times \mathbf{E} = \boldsymbol{0}\, , which implies \nabla \times \left( \mathbf{D} - \mathbf{P} \right) = \boldsymbol{0}\, , Applying
Helmholtz decomposition: \mathbf{D} - \mathbf{P} = -\nabla \Phi \, , for some scalar potential
φ, and: \nabla \cdot (\mathbf{D} - \mathbf{P}) = \varepsilon_0 \nabla \cdot \mathbf{E} = \rho_\text{f} + \rho_\text{b} = - \nabla^2 \Phi\, . Suppose the charges are divided into free and bound, and the potential is divided into \Phi = \Phi_\text{f} + \Phi_\text{b}\, . Satisfaction of the boundary conditions upon
φ may be divided arbitrarily between
φf and
φb because only the sum
φ must satisfy these conditions. It follows that
P is simply proportional to the electric field due to the charges selected as bound, with boundary conditions that prove convenient. In particular, when
no free charge is present, one possible choice is . Next is discussed how several different dipole moment descriptions of a medium relate to the polarization entering Maxwell's equations.
Medium with charge and dipole densities As described next, a model for polarization moment density
p(
r) results in a polarization \mathbf{P}(\mathbf{r}) = \mathbf{p}(\mathbf{r}) restricted to the same model. For a smoothly varying dipole moment distribution
p(
r), the corresponding bound charge density is simply \nabla \cdot \mathbf{p} (\mathbf{r}) = -\rho_\text{b}, as we will establish shortly via
integration by parts. However, if
p(
r) exhibits an abrupt step in dipole moment at a boundary between two regions, ∇·
p(
r) results in a surface charge component of bound charge. This surface charge can be treated through a
surface integral, or by using discontinuity conditions at the boundary, as illustrated in the various examples below. As a first example relating dipole moment to polarization, consider a medium made up of a continuous charge density
ρ(
r) and a continuous dipole moment distribution
p(
r). The potential at a position
r is: \Phi (\mathbf{r}) = \frac{1}{4 \pi \varepsilon_0} \int \frac{\rho \left(\mathbf{r}_0\right)}{\left|\mathbf{r} - \mathbf{r}_0\right|} d^3 \mathbf{r}_0 \ + \frac {1}{4 \pi \varepsilon_0}\int \frac{\mathbf{p} \left(\mathbf{r}_0\right) \cdot \left(\mathbf{r} - \mathbf{r}_0\right)} = \frac{\mathbf{r} - \mathbf{r}_0}{\left|\mathbf{r} - \mathbf{r}_0\right|^3} the polarization integral can be transformed: \begin{align} \int \frac{\mathbf{p} \left(\mathbf{r}_0\right) \cdot (\mathbf{r} - \mathbf{r}_0)}{\left|\mathbf{r} - \mathbf{r}_0\right|^3 } d^3 \mathbf{ r}_0 = {} & \int \mathbf{p} \left(\mathbf{r}_0\right) \cdot \nabla_{\mathbf{r}_0} \frac{1}{\left|\mathbf{r} - \mathbf{r}_0\right|} d^3 \mathbf{r}_0 , \\ ={} & \int \nabla_{\mathbf{r}_0} \cdot \frac {\mathbf{p} \left(\mathbf{r}_0\right)}{\left|\mathbf{r} - \mathbf{r}_0\right|} d^3 \mathbf{r}_0 - \int \frac{\nabla_{\mathbf{r}_0} \cdot \mathbf{p} \left(\mathbf{r}_0\right)}{\left|\mathbf{r} - \mathbf{r}_0\right|} d^3 \mathbf{r}_0 , \end{align} where the vector identity \nabla\cdot(\mathbf{A}{B}) = (\nabla\cdot\mathbf{A}){B} + \mathbf{A}\cdot(\nabla{B}) \implies \mathbf{A}\cdot(\nabla{B}) = \nabla\cdot(\mathbf{A}{B}) - (\nabla\cdot\mathbf{A}){B} was used in the last steps. The first term can be transformed to an integral over the surface bounding the volume of integration, and contributes a surface charge density, discussed later. Putting this result back into the potential, and ignoring the surface charge for now: \Phi (\mathbf{r}) = \frac{1}{4 \pi \varepsilon_0} \int \frac{\rho{\left(\mathbf{r}_0\right)} - \nabla_{\mathbf{r}_0} \cdot \mathbf{p}{\left(\mathbf{r}_0\right)}}{\left|\mathbf{r} - \mathbf{r}_0\right|} d^3 \mathbf{r}_0 , where the volume integration extends only up to the bounding surface, and does not include this surface. The potential is determined by the total charge, which the above shows consists of: \rho_\text{total} \left(\mathbf{r}_0\right) = \rho\left(\mathbf{r}_0\right) - \nabla_{\mathbf{r}_0} \cdot \mathbf{p} \left(\mathbf{r}_0\right)\, , showing that: -\nabla_{\mathbf{r}_0} \cdot \mathbf{p} \left(\mathbf{r}_0\right) = \rho_\text{b}\, . In short, the dipole moment density
p(
r) plays the role of the polarization density
P for this medium. Notice,
p(
r) has a non-zero divergence equal to the bound charge density (as modeled in this approximation). It may be noted that this approach can be extended to include all the multipoles: dipole, quadrupole, etc. Using the relation: \nabla \cdot \mathbf{D} = \rho_\text{f} \, , the polarization density is found to be: \mathbf{P}(\mathbf{r}) = \mathbf{p}_\text{dip} - \nabla \cdot \mathbf{p}_\text{quad} + \cdots\, , where the added terms are meant to indicate contributions from higher multipoles. Evidently, inclusion of higher multipoles signifies that the polarization density
P no longer is determined by a dipole moment density
p alone. For example, in considering scattering from a charge array, different multipoles scatter an electromagnetic wave differently and independently, requiring a representation of the charges that goes beyond the dipole approximation.
Surface charge Above, discussion was deferred for the first term in the expression for the potential due to the dipoles. Integrating the divergence results in a surface charge. The figure at the right provides an intuitive idea of why a surface charge arises. The figure shows a uniform array of identical dipoles between two surfaces. Internally, the heads and tails of dipoles are adjacent and cancel. At the bounding surfaces, however, no cancellation occurs. Instead, on one surface the dipole heads create a positive surface charge, while at the opposite surface the dipole tails create a negative surface charge. These two opposite surface charges create a net electric field in a direction opposite to the direction of the dipoles. This idea is given mathematical form using the potential expression above. Ignoring the free charge, the potential is: \Phi\left(\mathbf{r}\right) = \frac{1}{4 \pi \varepsilon_0} \int \nabla_{\mathbf{r}_0} \cdot \left(\mathbf{p} \left(\mathbf{r}_0\right) \frac{1}{\left|\mathbf{r} - \mathbf{r}_0\right|} \right) d^3 \mathbf{r}_0 - \frac{1}{4 \pi \varepsilon_0} \int \frac{\nabla_{\mathbf{r}_0} \cdot \mathbf{p} \left(\mathbf{r}_0\right)}{\left|\mathbf{r} - \mathbf{r}_0\right|} d^3 \mathbf{r}_0\, . Using the
divergence theorem, the divergence term transforms into the surface integral: \int \nabla_{\mathbf{r}_0} \cdot \left(\mathbf{p} \left(\mathbf{r}_0\right) \frac{1}{\left|\mathbf{r} - \mathbf{r}_0\right|}\right) d^3\mathbf{r}_0 = \int \frac{\mathbf{p} \left(\mathbf{r}_0\right) \cdot d \mathbf{A}_0}\left|\mathbf{r} - \mathbf{r}_0\right| \, , with d
A0 an element of surface area of the volume. In the event that
p(
r) is a constant, only the surface term survives: \Phi(\mathbf{r}) = \frac{1}{4 \pi \varepsilon_0} \int \frac{1}{\left|\mathbf{r} - \mathbf{r}_0\right|}\ \mathbf{p} \cdot d\mathbf{A}_0 \, , with d
A0 an elementary area of the surface bounding the charges. In words, the potential due to a constant
p inside the surface is equivalent to that of a
surface charge \sigma = \mathbf{p} \cdot d \mathbf{A} which is positive for surface elements with a component in the direction of
p and negative for surface elements pointed oppositely. (Usually the direction of a surface element is taken to be that of the outward normal to the surface at the location of the element.) If the bounding surface is a sphere, and the point of observation is at the center of this sphere, the integration over the surface of the sphere is zero: the positive and negative surface charge contributions to the potential cancel. If the point of observation is off-center, however, a net potential can result (depending upon the situation) because the positive and negative charges are at different distances from the point of observation.{{refn|group=notes|name="multipole1"|refn=A brute force evaluation of the integral can be done using a multipole expansion: \frac{1}{\left|\mathbf{r} - \mathbf{r}_0\right|} = \sum_{\ell,\ m} \frac{4\pi}{2\ell + 1} \frac{1}{r} \left(\frac{r_0}{r}\right)^\ell {Y^*}_{\ell}^m \left(\theta_0,\ \varphi_0\right) Y_{\ell}^m \left(\theta,\ \varphi\right). }} The field due to the surface charge is: \mathbf{E}\left(\mathbf{r}\right) = -\frac{1}{4 \pi \varepsilon_0} \nabla_\mathbf{r} \int \frac{1}{\left|\mathbf{r} - \mathbf{r}_0\right|}\ \mathbf{p} \cdot d\mathbf{A}_0\, , which, at the center of a spherical bounding surface is not zero (the
fields of negative and positive charges on opposite sides of the center add because both fields point the same way) but is instead: In the case when the polarization is
outside a spherical cavity, the field in the cavity due to the surrounding dipoles is in the
same direction as the polarization. In particular, if the
electric susceptibility is introduced through the approximation: \mathbf{p}(\mathbf{r}) = \varepsilon_0 \chi(\mathbf{r}) \mathbf{E}(\mathbf{r})\, , where , in this case and in the following, represent the
external field which induces the polarization. Then: \nabla \cdot \mathbf{p}(\mathbf{r}) = \nabla \cdot \left(\chi(\mathbf{r}) \varepsilon_0 \mathbf{E}(\mathbf{r})\right) = -\rho_\text{b}\, . Whenever
χ(
r) is used to model a step discontinuity at the boundary between two regions, the step produces a surface charge layer. For example, integrating along a normal to the bounding surface from a point just interior to one surface to another point just exterior: \varepsilon_0 \hat{\mathbf{n}} \cdot \left[\chi\left(\mathbf{r}_+\right) \mathbf{E}\left(\mathbf{r}_+\right) - \chi\left(\mathbf{r}_-\right) \mathbf{E}\left(\mathbf{r}_-\right)\right] = \frac{1}{A_\text{n}} \int d \Omega_\text{n}\ \rho_\text{b} = 0 \, , where
An, Ωn indicate the area and volume of an elementary region straddling the boundary between the regions, and \hat{\mathbf{n}} a unit normal to the surface. The right side vanishes as the volume shrinks, inasmuch as
ρb is finite, indicating a discontinuity in
E, and therefore a surface charge. That is, where the modeled medium includes a step in permittivity, the polarization density corresponding to the dipole moment density \mathbf{p}(\mathbf{r}) = \varepsilon_0 \chi(\mathbf{r}) \mathbf{E}(\mathbf{r}) necessarily includes the contribution of a surface charge. A physically more realistic modeling of
p(
r) would have the dipole moment density drop off rapidly, but smoothly to zero at the boundary of the confining region, rather than making a sudden step to zero density. Then the surface charge will not concentrate in an infinitely thin surface, but instead, being the divergence of a smoothly varying dipole moment density, will distribute itself throughout a thin, but finite transition layer.
Dielectric sphere in uniform external electric field s of the
D-field in a dielectric sphere with greater susceptibility than its surroundings, placed in a previously uniform field. The
field lines of the
E-field (not shown) coincide everywhere with those of the
D-field, but inside the sphere, their density is lower, corresponding to the fact that the
E-field is weaker inside the sphere than outside. Many of the external
E-field lines terminate on the surface of the sphere, where there is a bound charge. The above general remarks about surface charge are made more concrete by considering the example of a dielectric sphere in a uniform electric field. The sphere is found to adopt a surface charge related to the dipole moment of its interior. A uniform external electric field is supposed to point in the
z-direction, and spherical polar coordinates are introduced so the potential created by this field is: \Phi_\infty = -E_\infty z = -E_\infty r \cos\theta \, . The sphere is assumed to be described by a
dielectric constant κ, that is, \mathbf{D} = \kappa \varepsilon_0 \mathbf{E} \, , and inside the sphere the potential satisfies Laplace's equation. Skipping a few details, the solution inside the sphere is: \Phi_ while outside the sphere: \Phi_> = \left(Br + \frac{C}{r^2} \right) \cos\theta \, . At large distances, so . Continuity of potential and of the radial component of displacement determine the other two constants. Supposing the radius of the sphere is
R, A = -\frac{3}{\kappa + 2} E_\infty\ ;\ C = \frac{\kappa - 1}{\kappa + 2} E_\infty R^3\, , As a consequence, the potential is: \Phi_> = \left(-r + \frac{\kappa - 1}{\kappa + 2} \frac{R^3}{r^2}\right) E_\infty \cos\theta\, , which is the potential due to applied field and, in addition, a dipole in the direction of the applied field (the
z-direction) of dipole moment: \mathbf{p} = 4 \pi \varepsilon_0 \left(\frac{\kappa - 1}{\kappa + 2} R^3\right) \mathbf{E}_\infty\, , or, per unit volume: \frac{\mathbf{p}}{V} = 3 \varepsilon_0 \left(\frac{\kappa - 1}{\kappa + 2}\right) \mathbf{E}_\infty\, . The factor is called the
Clausius–Mossotti factor and shows that the induced polarization flips sign if . Of course, this cannot happen in this example, but in an example with two different dielectrics
κ is replaced by the ratio of the inner to outer region dielectric constants, which can be greater or smaller than one. The potential inside the sphere is: \Phi_ leading to the field inside the sphere: -\nabla \Phi_ showing the depolarizing effect of the dipole. Notice that the field inside the sphere is
uniform and parallel to the applied field. The dipole moment is uniform throughout the interior of the sphere. The surface charge density on the sphere is the difference between the radial field components: \sigma = 3 \varepsilon_0 \frac{\kappa - 1}{\kappa + 2} E_\infty \cos\theta = \frac{1}{V} \mathbf{p} \cdot \hat{\mathbf{R}}\, . This linear dielectric example shows that the dielectric constant treatment is equivalent to the uniform dipole moment model and leads to zero charge everywhere except for the surface charge at the boundary of the sphere.
General media If observation is confined to regions sufficiently remote from a system of charges, a multipole expansion of the exact polarization density can be made. By truncating this expansion (for example, retaining only the dipole terms, or only the dipole and quadrupole terms, or
etc.), the results of the previous section are regained. In particular, truncating the expansion at the dipole term, the result is indistinguishable from the polarization density generated by a uniform dipole moment confined to the charge region. To the accuracy of this dipole approximation, as shown in the previous section, the dipole moment
density p(
r) (which includes not only
p but the location of
p) serves as
P(
r). At locations
inside the charge array, to connect an array of paired charges to an approximation involving only a dipole moment density
p(
r) requires additional considerations. The simplest approximation is to replace the charge array with a model of ideal (infinitesimally spaced) dipoles. In particular, as in the example above that uses a constant dipole moment density confined to a finite region, a surface charge and depolarization field results. A more general version of this model (which allows the polarization to vary with position) is the customary approach using
electric susceptibility or
electrical permittivity. A more complex model of the point charge array introduces an
effective medium by averaging the microscopic charges; A related approach is to divide the charges into those nearby the point of observation, and those far enough away to allow a multipole expansion. The nearby charges then give rise to
local field effects. In a common model of this type, the distant charges are treated as a homogeneous medium using a dielectric constant, and the nearby charges are treated only in a dipole approximation. The approximation of a medium or an array of charges by only dipoles and their associated dipole moment density is sometimes called the
point dipole approximation, the
discrete dipole approximation, or simply the
dipole approximation. == Electric dipole moments of fundamental particles ==