# Electric dipole moment: Wikis

Note: Many of our articles have direct quotes from sources you can cite, within the Wikipedia article! This article doesn't yet, but we're working on it! See more info or our list of citable articles.

# Encyclopedia

A water molecule. A molecule of water is polar because of the unequal sharing of its electrons in a "bent" structure. A separation of charge is present with negative charge in the middle (red shade), and positive charge at the ends (blue shade).

In physics, the electric dipole moment is a measure of the separation of positive and negative electrical charges in a system of charges, that is, a measure of the charge system's overall polarity.

In the simple case of two point charges, one with charge + q and one with charge q, the electric dipole moment p is:

$\boldsymbol{p} = q \, \boldsymbol{d}$

where d is the displacement vector pointing from the negative charge to the positive charge. Thus, the electric dipole moment vector p points from the negative charge to the positive charge. There is no inconsistency here, because the electric dipole moment has to do with orientation of the dipole, that is, the positions of the charges, and does not indicate the direction of the field originating in these charges.

An idealization of this two-charge system is the electrical point dipole consisting of two (infinite) charges only infinitesimally separated, but with a finite p = q d.

## General case

More generally, for a continuous distribution of charge confined to a volume V, the corresponding expression for the dipole moment is:

$\boldsymbol{p}(\boldsymbol{r}) = \int_{V} \rho(\boldsymbol{r_0})\, (\boldsymbol{r_0}-\boldsymbol{r}) \ d^3 \boldsymbol{r_0},$

where r locates the point of observation and d3r0 denotes an elementary volume in V. For an array of point charges, the charge density becomes a sum of Dirac delta functions:

$\rho (\boldsymbol{r}) = \sum_{i=1}^N \, q_i \, \delta (\mathbf{r} - \mathbf{r}_i ) \ ,$

where each $\mathbf{r}_i$ is a vector from some reference point to the charge qi. Substitution into the above integration formula provides:

$\boldsymbol{p}(\boldsymbol{r}) = \sum_{i=1}^N \, q_i \int\delta (\mathbf{r_0} - \mathbf{r}_i )\, (\boldsymbol{r_0}-\boldsymbol{r}) \ d^3 \boldsymbol{r_0}$$= \sum_{i=1}^N \, q_i (\boldsymbol{r_i}-\boldsymbol{r}),$

This expression is equivalent to the previous expression in the case of charge neutrality and N = 2. For two opposite charges, denoting the location of the positive charge of the pair as $\boldsymbol {r_+}$ and the location of the negative charge as $\boldsymbol {r_-}$ :

$\boldsymbol{p}(\boldsymbol{r})$$=q_1(\boldsymbol{r_1}-\boldsymbol{r})+q_2(\boldsymbol{r_2}-\boldsymbol{r})$$= q(\boldsymbol{r_+}-\boldsymbol{r})-q(\boldsymbol{r_-}-\boldsymbol{r})$$=q (\boldsymbol{r_+} - \boldsymbol{r_-})=q\boldsymbol d \ ,$

showing that the dipole moment vector is directed from the negative charge to the positive charge because the position vector of a point is directed outward from the origin to that point.

The dipole moment is most easily understood when the system has an overall neutral charge; for example, a pair of opposite charges, or a neutral conductor in a uniform electric field. For a system of charges with no net charge, visualized as an array of paired opposite charges, the relation for electric dipole moment is:

$\boldsymbol{p}(\boldsymbol{r}) = \sum_{i=1}^{N} \, \int q_i \left( \delta (\mathbf{r_0} - (\mathbf{r}_i + \boldsymbol{d_i}) )- \delta ( \mathbf{r_0} - \mathbf{r}_i ) \right)\, (\boldsymbol{r_0}-\boldsymbol{r}) \ d^3 \boldsymbol{r_0}$$= \sum_{i=1}^{N} \, q_i \left( \boldsymbol{r_i +d_i}-\boldsymbol{r} -(\boldsymbol{r_i }-\boldsymbol{r}) \right) = \sum_{i=1}^{N} q_i\boldsymbol{d}_i \, = \sum_{i=1}^{N} \boldsymbol{p}_i \ ,$

which is the vector sum of the individual dipole moments of the neutral charge pairs. (Because of overall charge neutrality, the dipole moment is independent of the observer's position r.) Thus, the value of p is independent of the choice of reference point, provided the overall charge of the system is zero.

When discussing the dipole moment of a non-neutral system, such as the dipole moment of the proton, a dependence on the choice of reference point arises. In such cases it is conventional to choose the reference point to be the center of mass of the system or the center of charge, not some arbitrary origin.[1] This convention ensures that the dipole moment is an intrinsic property of the system.

## Potential and field of an electric dipole

An ideal dipole consists of two charges with infinitesimal separation. The potential and field of such an ideal dipole are found next as a limiting case of an example of two opposite charges at non-zero separation.

Two closely spaced opposite charges have a potential of the form:

$\phi ( \boldsymbol{r} )=\frac {q}{4 \pi \varepsilon _0 | \boldsymbol{ r}- \boldsymbol{r}_1 |} -\frac {q}{4 \pi \varepsilon _0 | \boldsymbol{ r}- \boldsymbol{r}_2 | } \ ,$

with

$\boldsymbol{r}_2 - \boldsymbol{r}_1 =\boldsymbol d \ ,$

with d the charge separation. Introducing the radius to the center of charge, say R, and the unit vector in the direction of R:

${\boldsymbol {R}} = \frac{\boldsymbol{r}_1 + \boldsymbol{r}_2}{2}\ ; \ \boldsymbol {\hat{R}} = \frac {\boldsymbol {R}}{R} \ ,$

some mathematical manipulation (see multipole expansion and quadrupole) allows this potential to be expressed as a series in d/R as:[2][3]

$\phi ( \boldsymbol{R} )=\frac {1}{4 \pi \varepsilon _0} \frac {q\boldsymbol {d \cdot \hat{R}}}{R^2} + \mathrm{ other \ terms } \approx \frac {1}{4 \pi \varepsilon _0} \frac {\boldsymbol {p\cdot \hat{R}}}{R^2} \ ,$

where the other terms in the series are small at distances R large enough to make d/R small.[4] Here, the electric dipole moment p is, as above:

$\boldsymbol p = q \boldsymbol d \ .$

The result for the dipole potential also can be expressed as:[5]

$\phi ( \boldsymbol{R} )=- \boldsymbol {p\cdot \nabla}\frac {1}{4 \pi \varepsilon _0 R}\ ,$

which relates the dipole potential to that of a point charge. A key point is that the potential of the dipole falls off faster with distance R than that of the point charge.

The field of the dipole is the gradient of the potential, leading to:[5]

$\boldsymbol E = \frac {3 \boldsymbol {p \cdot \hat{R}}}{4 \pi \varepsilon_0 R^3} \boldsymbol { \hat{R}}-\frac {\boldsymbol{p}}{4 \pi \varepsilon_0 R^3} \ .$

Thus, although two closely spaced opposite charges are not an ideal electric dipole (because their potential at close approach is not that of a dipole), at distances much larger than their separation, their dipole moment p appears directly in their potential and field.

As the two charges are brought closer together (d is made smaller), the dipole term in the multipole expansion based on the ratio d/R becomes the only significant term at ever closer distances R, and in the limit of infinitesimal separation the dipole term in this expansion is all that matters. As d is made infinitesimal, however, the dipole charge must be made to increase to hold p constant. This limiting process results in a "point dipole".

## Dipole moment density and polarization density

The dipole moment of an array of charge,

$\boldsymbol p = \sum_{i=1}^N \ q_i \boldsymbol {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 directions about where the array happens to be located. 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 P(r) = p(r). 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 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:

$\boldsymbol{D} = \varepsilon _0 \boldsymbol{E} + \boldsymbol{P}\ ,$

where P is called the polarization density. In this formulation, the divergence of this equation yields:

$\nabla \cdot \boldsymbol{D} = \rho_f = \varepsilon _0 \nabla \cdot \boldsymbol{E} +\nabla \cdot \boldsymbol{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 \boldsymbol{P} = -\rho_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 curl E = 0, which implies curl (D − P) = 0. Applying Helmholtz decomposition:[6]

$\boldsymbol{ (D-P) = -\nabla } \varphi \ ,$

for some scalar potential φ, and:

$\boldsymbol {\nabla \cdot (D-P)} =\varepsilon_0 \boldsymbol {\nabla \cdot E}=\rho_f +\rho_b = -\nabla ^2 \varphi \ .$

Suppose the charges are divided into free and bound, and the potential is divided into φ = φf + φ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 the electric field of the charges selected as bound, with boundary conditions that prove convenient.[7][8] In particular, when no free charge is present, one possible choice is P = ε0 E.

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 P(r) = p(r) restricted to the same model. For a smoothly varying dipole moment distribution p(r), the corresponding bound charge density is simply ∇·p(r) = −ρb. However, in the case of a p(r) that exhibits an abrupt step in dipole moment at a boundary between two regions, ∇·p(r) exhibits 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).[9] The potential at a position r is:[10][11]

$\phi ( \boldsymbol{r} ) = \frac {1}{4 \pi \varepsilon_0}\int \frac { \rho ( \boldsymbol{ r}_0 )} {| \boldsymbol{ r}- \boldsymbol{r}_0 | } d^3 \boldsymbol{ r}_0 \ + \frac {1}{4 \pi \varepsilon_0}\int \frac { \boldsymbol{p} ( \boldsymbol{ r}_0 )\boldsymbol{\cdot (r - r_0)}} {| \boldsymbol{ r}- \boldsymbol{r}_0 |^3 } d^3 \boldsymbol{ r}_0 ,$

where ρ(r) is the unpaired charge density, and p(r) is the dipole moment density.[12] Using an identity:

$\nabla_{\boldsymbol {r}_0} \frac {1}{|\boldsymbol r - \boldsymbol{r}_0|} = \frac {\boldsymbol r - \boldsymbol{r}_0}{|\boldsymbol r - \boldsymbol{r}_0|^3}$

the polarization integral can be transformed:

$\frac {1}{4 \pi \varepsilon_0}\int \frac { \boldsymbol{p} ( \boldsymbol{ r}_0 )\boldsymbol{\cdot (r - r_0)}} {| \boldsymbol{ r}- \boldsymbol{r}_0 |^3 } d^3 \boldsymbol{ r}_0 =\frac {1}{4 \pi \varepsilon_0}\int \boldsymbol{p} ( \boldsymbol{ r}_0 )\boldsymbol{\cdot \nabla}_{\boldsymbol {r}_0} \frac {1}{|\boldsymbol r - \boldsymbol{r}_0|} d^3 \boldsymbol{ r}_0 ,$
$=\frac {1}{4 \pi \varepsilon_0}\int \boldsymbol{\nabla_{\boldsymbol {r_0}}\cdot} \left( \boldsymbol{p} ( \boldsymbol{ r}_0 ) \frac {1}{|\boldsymbol r - \boldsymbol{r}_0|} \right) d^3 \boldsymbol{ r}_0 -\frac {1}{4 \pi \varepsilon_0}\int \frac {\boldsymbol{\nabla_{\boldsymbol {r_0}}\cdot} \boldsymbol{p} ( \boldsymbol{ r}_0 )}{|\boldsymbol r - \boldsymbol{r}_0|} d^3 \boldsymbol{ r}_0 ,$

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 ( \boldsymbol{r} ) = \frac {1}{4 \pi \varepsilon_0}\int \frac { \rho ( \boldsymbol{ r}_0 )-\boldsymbol{\nabla_{\boldsymbol {r_0}}\cdot} \boldsymbol{p} ( \boldsymbol{ r}_0 )} {| \boldsymbol{ r}- \boldsymbol{r}_0 | } d^3 \boldsymbol{ 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_{total} (\boldsymbol{ r}_0)= \rho ( \boldsymbol{ r}_0 )-\boldsymbol{\nabla_{\boldsymbol {r_0}}\cdot} \boldsymbol{p} ( \boldsymbol{ r}_0 ) \ ,$

showing that:

$-\boldsymbol{\nabla_{\boldsymbol {r_0}}\cdot} \boldsymbol{p} ( \boldsymbol{ r}_0 ) = \rho_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.[13][14] Using the relation:

$\nabla \cdot \boldsymbol{D} = \rho_f \ ,$

the polarization density is found to be:

$\boldsymbol{P( r )} = -\boldsymbol{\nabla \cdot p_{Dip}} -\boldsymbol{\nabla \cdot p_{Quad}} + ... \ ,$

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. 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.[15]

#### Surface charge

A uniform array of identical dipoles is equivalent to a surface charge.

Above, discussion was deferred for the leading divergence term in the expression for the potential due to the dipoles. This term 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. The potential is:

$\phi ( \boldsymbol{r} ) =\frac {1}{4 \pi \varepsilon_0}\int \boldsymbol{\nabla_{\boldsymbol {r_0}}\cdot} \left( \boldsymbol{p} ( \boldsymbol{ r}_0 ) \frac {1}{|\boldsymbol r - \boldsymbol{r}_0|} \right) d^3 \boldsymbol{ r}_0$$-\frac {1}{4 \pi \varepsilon_0}\int \frac {\boldsymbol{\nabla_{\boldsymbol {r_0}}\cdot} \boldsymbol{p} ( \boldsymbol{ r}_0 )}{|\boldsymbol r - \boldsymbol{r}_0|} d^3 \boldsymbol{ r}_0 \ .$

Using the divergence theorem, the divergence term transforms into the surface integral:

$\frac {1}{4 \pi \varepsilon_0}\int \boldsymbol{\nabla_{\boldsymbol {r_0}}\cdot} \left( \boldsymbol{p} ( \boldsymbol{ r}_0 ) \frac {1}{|\boldsymbol r - \boldsymbol{r}_0|} \right) d^3 \boldsymbol{ r}_0$
$=\frac {1}{4 \pi \varepsilon_0}\int \frac {\boldsymbol{p} ( \boldsymbol{ r}_0 )\boldsymbol{\cdot } d \boldsymbol {A_0 } } {|\boldsymbol r - \boldsymbol{r}_0|} \ ,$

with dAo an element of surface area of the volume. In the event that p(r) is a constant, only the surface term survives:

$\phi ( \boldsymbol{r} )$$=\frac {1}{4 \pi \varepsilon_0}\int \frac {1}{|\boldsymbol r - \boldsymbol{r}_0|}\ \boldsymbol{p} \cdot d\boldsymbol{A_0} \ ,$

with dAo 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 σ = dA, 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.[16] The field due to the surface charge is:

$\boldsymbol E ( \boldsymbol{r} ) =-\frac {1}{4 \pi \varepsilon_0} \boldsymbol{\nabla}_{\boldsymbol {r}}\int \frac {1}{|\boldsymbol r - \boldsymbol{r}_0|}\ \boldsymbol{p} \cdot d\boldsymbol{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 :[17]

$\boldsymbol E =-\frac {\boldsymbol p}{3 \varepsilon_0} \ .$

If we suppose the polarization of the dipoles was induced by an external field, the polarization field opposes the applied field and sometimes is called a depolarization field.[18][19] 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.[20]

In particular, if the electric susceptibility is introduced through the approximation:

$\boldsymbol{p(r)} = \varepsilon_0 \chi(\boldsymbol r ) \boldsymbol {E(r)} \ ,$

then:

$\boldsymbol { \nabla \cdot p(r)}=\boldsymbol { \nabla \cdot} \left( \chi \boldsymbol{ (r)}\varepsilon_0 \boldsymbol {E(r)}\right) =-\rho_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{\boldsymbol n} \cdot \left( \chi \boldsymbol{ (r_+)}\boldsymbol {E(r_+)}-\chi \boldsymbol{ (r_-)}\boldsymbol {E(r_-)}\right) =\frac{1}{A_n} \int d \Omega_n \ \rho_b = 0 \ ,$

where An, Ωn indicate the area and volume of an elementary region straddling the boundary between the regions, and $\hat{\boldsymbol 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 p(r) = χ(r)E(r) necessarily includes the contribution of a surface charge.[21][22][23]

It may be noted that a physically more realistic modeling of p(r) would cause the dipole moment density to taper off continuously to zero at the boundary of the confining region, rather than making a sudden step to zero density. Then the surface charge becomes zero at the boundary, and the surface charge is replaced by the divergence of a continuously varying dipole-moment density.

#### Dielectric sphere in uniform external electric field

Electric field lines in dielectric sphere with greater susceptibility than its surroundings.[24]

The above general remarks about surface charge are made more concrete by considering the example of a dielectric sphere in a uniform electric field.[25][26] 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, D = κε0E, and inside the sphere the potential satisfies Laplace's equation. Skipping a few details, the solution inside the sphere is:

$\phi_< = A r \cos \theta \ ,$

while outside the sphere:

$\phi_> = \left(Br + \frac {C}{r^2} \right ) \cos \theta \ .$

At large distances, φ> → φ so B = -E . Continuity of potential and of the radial component of displacement D = κε0E 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:

$\boldsymbol p = 4 \pi \varepsilon_0 \left(\frac {\kappa-1}{\kappa+2}{R^3} \right) \boldsymbol{E_{\infty}} \ ,$

or, per unit volume:

$\frac {\boldsymbol p}{V} = {3}\varepsilon_0 \left(\frac {\kappa-1}{\kappa+2}\right) \boldsymbol{E_{\infty}} \ .$

The factor (κ-1)/(κ+2) is called the Clausius-Mossotti factor and shows that the induced polarization flips sign if κ < 1. 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_< = -\frac{3}{\kappa +2} E_{\infty}r \cos \theta \ ,$

leading to the field inside the sphere:

$\boldsymbol {-\nabla} \phi_< = \frac{3}{\kappa +2} \boldsymbol{ E_{\infty}} =\left( 1-\frac {\kappa-1}{\kappa+2} \right)\boldsymbol{ E_{\infty}} \ ,$

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} \boldsymbol{ p \cdot \hat{R}}\ .$

his example shows that, for this example, 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 a electric susceptibility or electrical permittivity.

A more complex model of the point charge array introduces an effective medium by averaging the microscopic charges;[19] for example, the averaging can arrange that only dipole fields play a role.[27][28] 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.[17][29] 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.[30] 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.[31][32][33]

## Dipole moments of fundamental particles

Much experimental work is continuing on measuring the electric dipole moments (EDM) of fundamental and composite particles, namely those of the neutron and electron. As EDMs violate both the Parity (P) and Time (T) symmetries, their values yield a mostly model-independent measure (assuming CPT symmetry is valid) of CP-violation in nature. Therefore, values for these EDMs place strong constraints upon the scale of CP-violation that extensions to the standard model of particle physics may allow.

Indeed, many theories are inconsistent with the current limits and have effectively been ruled out, and established theory permits a much larger value than these limits, leading to the strong CP problem and prompting searches for new particles such as the axion.

Current generations of experiments are designed to be sensitive to the supersymmetry range of EDMs, providing complementary experiments to those done at the LHC.

## References and in-line notes

1. ^ Christopher J. Cramer (2004). Essentials of computational chemistry (2 ed.). Wiley. p. 307. ISBN 0470091827.
2. ^ David E Dugdale (1993). Essentials of Electromagnetism. Springer. pp. 80–81. ISBN 1563962535.
3. ^ Kikuji Hirose, Tomoya Ono, Yoshitaka Fujimoto (2005). First-principles calculations in real-space formalism. Imperial College Press. p. 18. ISBN 1860945120.
4. ^ Each succeeding term provides a more detailed view of the distribution of charge, and falls off more rapidly with distance. For example, the quadrupole moment is the basis for the next term: $Q_{ij} = \int d^3 \boldsymbol{r_0} \left( 3x_i x_j -r_0^2 \delta_{ij} \right) \rho( \boldsymbol{r_0}) \ ,$ with ro = (x1, x2, x3). See HW Wyld (1999). Mathematical Methods for Physics. Westview Press. p. 106. ISBN 0738201251.
5. ^ a b BB Laud (1987). Electromagnetics (2 ed.). New Age International. p. 25. ISBN 0852264992.
6. ^ Jie-Zhi Wu, Hui-Yang Ma, Ming-De Zhou (2006). "§2.3.1 Functionally Orthogonal Decomposition". Vorticity and vortex dynamics. Springer. pp. 36 ff. ISBN 3540290273.
7. ^ For example, one could place the boundary around the bound charges at infinity. Then φb falls off with distance from the bound charges. If an external field is present, and zero free charge, the field can be accounted for in the contribution of φf, which would arrange to satisfy the boundary conditions and Laplace's equation2φf = 0.
8. ^ In principle, one could add the same arbitrary curl to both D and P, which would cancel out of the difference D − P. However, assuming D and P originate in in a simple division of charges into free and bound, they are at bottom electric fields and so have zero curl.
9. ^ This medium can be seen as an idealization growing from the multipole expansion of the potential of an arbitrarily complex charge distribution, truncation of the expansion, and the forcing of the truncated form to apply everywhere. The result is a hypothetical medium. See Jack Vanderlinde (2004). "§7.1 The electric field due to a polarized dielectric". Classical Electromagnetic Theory. Springer. ISBN 1402026994.
10. ^ Uwe Krey, Anthony Owen (2007). Basic Theoretical Physics: A Concise Overview. Springer. pp. 138–143. ISBN 3540368043.
11. ^ T Tsang (1997). Classical Electrodynamics. World Scientific. p. 59. ISBN 9810230419.
12. ^ For example, for a system of ideal dipoles with dipole moment p confined within some closed surface, the dipole density p(r) is equal to p inside the surface, but is zero outside. That is, the dipole density includes a Heaviside step function locating the dipoles inside the surface.
13. ^ George E Owen (2003). Introduction to Electromagnetic Theory (republication of the 1963 Allyn & Bacon ed.). Courier Dover Publications. p. 80. ISBN 0486428303.
14. ^ Pierre-François Brevet (1997). Surface second harmonic generation. Presses polytechniques et universitaires romandes. p. 24. ISBN 2880743451.
15. ^ See Daniel A. Jelski, Thomas F. George (1999). Computational studies of new materials. World Scientific. p. 219. ISBN 9810233256.  and EM Purcell & CR Pennypacker (1973). "Scattering and Absorption of Light by Nonspherical Dielectric Grains". Astrophysical Journal 186: 705–714. doi:10.1086/152538.
16. ^ A brute force evaluation of the integral can be done using a multipole expansion: $\frac{1}{|\boldsymbol{r-r_0}|}$ = $\sum_{\ell,\ m } \frac{4\pi}{2 \ell +1}$·$\frac {1}{r} \left({\frac {r_0}{r}}\right)^{\ell}$·$\ {Y^*}_{\ell}^m (\theta_0 , \ \phi_0) Y_{\ell}^m (\theta, \ \phi)$. See HW Wyld (1999). Mathematical Methods for Physics. Westview Press. p. 104. ISBN 0738201251.
17. ^ a b H. Ibach, Hans Lüth (2003). Solid-state Physics: an introduction to principles of materials science (3 ed.). Springer. p. 361. ISBN 354043870X.
18. ^ Yasuaki Masumoto, Toshihide Takagahara (2002). Semiconductor quantum dots: physics, spectroscopy, and applications. Springer. p. 72. ISBN 3540428054.
19. ^ a b Yutaka Toyozawa (2003). Optical processes in solids. Cambridge University Press. p. 96. ISBN 0521556058.
20. ^ For example, a droplet in a surrounding medium experiences a higher or a lower internal field depending upon whether the medium has a higher or a lower dielectric constant than that of the droplet. See Paul S. Drzaic (1995). Liquid crystal dispersions. World Scientific. p. 246. ISBN 9810217455.
21. ^ Wai-Kai Chen (2005). The electrical engineering handbook. Academic Press. p. 502. ISBN 0121709604.
22. ^ Julius Adams Stratton (2007). Electromagnetic theory (reprint of 1941 ed.). Wiley-IEEE. p. 184. ISBN 0470131535.
23. ^ Edward J. Rothwell, Michael J. Cloud (2001). Electromagnetics. CRC Press. p. 68. ISBN 084931397X.
24. ^ Based upon equations from Andrew Gray (1888). The theory and practice of absolute measurements in electricity and magnetism. Macmillan & Co.. pp. 126–127. , which refers to papers by Sir W. Thomson.
25. ^ HW Wyld (1999). Mathematical Methods for Physics (2 ed.). Westview Press. pp. 233 ff. ISBN 0738201251.
26. ^ Julius Adams Stratton (2007). Electromagnetic theory (Wiley-IEEE reissue ed.). Piscataway, NJ: IEEE Press. p. 205 ff. ISBN 0470131535.
27. ^ John E Swipe & RW Boyd (2002). "Nanocomposite materials for nonlinear optics based upon local field effects". in Vladimir M. Shalaev. Optical properties of nanostructured random media. Springer. p. 3. ISBN 3540420312.
28. ^ Emil Wolf (1977). Progress in Optics. Elsevier. p. 288. ISBN 0720415152.
29. ^ Mark Fox (2006). Optical Properties of Solids. Oxford University Press. p. 39. ISBN 0198506120.
30. ^ Lev Kantorovich (2004). "§8.2.1 The local field". Quantum theory of the solid state. Springer. p. 426. ISBN 1402021534.
31. ^ Pierre Meystre (2001). Atom Optics. Springer. p. 5. ISBN 0387952748.
32. ^ Bruce T Draine (2001). "The discrete dipole approximation for light scattering by irregular targets". in Michael I. Mishchenko. Light scattering by nonspherical particles. Academic Press. p. 132. ISBN 0124986609.
33. ^ MA Yurkin & AG Hoekstra (2007). "The discrete dipole approximation: an overview and recent developments". Journal of Quantitative Spectroscopy and Radiative Transfer 106: 558–589. doi:10.1016/j.jqsrt.2007.01.034.

Melvin Schwartz (1987). "Electrical dipole moment". Principles of Electrodynamics (reprint of 1972 ed.). Courier Dover Publications. p. 49ff. ISBN 0486654931.

. A molecule of water is polar because of the unequal sharing of its electrons in a "bent" structure. A separation of charge is present with negative charge in the middle (red shade), and positive charge at the ends (blue shade).]] In physics, the electric dipole moment is a measure of the separation of positive and negative electrical charges in a system of charges, that is, a measure of the charge system's overall polarity.

In the simple case of two point charges, one with charge $\left\{+\right\}q$ and one with charge $\left\{-\right\}q$, the electric dipole moment p is:


 \boldsymbol{p} = q \, \boldsymbol{d}


where d is the displacement vector pointing from the negative charge to the positive charge. Thus, the electric dipole moment vector p points from the negative charge to the positive charge. There is no inconsistency here, because the electric dipole moment has to do with orientation of the dipole, that is, the positions of the charges, and does not indicate the direction of the field originating in these charges.

An idealization of this two-charge system is the electrical point dipole consisting of two (infinite) charges only infinitesimally separated, but with a finite p = q d.

## General case

More generally, for a continuous distribution of charge confined to a volume V, the corresponding expression for the dipole moment is:

$\boldsymbol\left\{p\right\}\left(\boldsymbol\left\{r\right\}\right) = \int_\left\{V\right\} \rho\left(\boldsymbol\left\{r_0\right\}\right)\, \left(\boldsymbol\left\{r_0\right\}-\boldsymbol\left\{r\right\}\right) \ d^3 \boldsymbol\left\{r_0\right\},$

where r locates the point of observation and d3r0 denotes an elementary volume in V. For an array of point charges, the charge density becomes a sum of Dirac delta functions:

$\rho \left(\boldsymbol\left\{r\right\}\right) = \sum_\left\{i=1\right\}^N \, q_i \, \delta \left(\mathbf\left\{r\right\} - \mathbf\left\{r\right\}_i \right) \ ,$

where each $\mathbf\left\{r\right\}_i$ is a vector from some reference point to the charge $q_i$. Substitution into the above integration formula provides:

$\boldsymbol\left\{p\right\}\left(\boldsymbol\left\{r\right\}\right) = \sum_\left\{i=1\right\}^N \, q_i \int\delta \left(\mathbf\left\{r_0\right\} - \mathbf\left\{r\right\}_i \right)\, \left(\boldsymbol\left\{r_0\right\}-\boldsymbol\left\{r\right\}\right) \ d^3 \boldsymbol\left\{r_0\right\}$$= \sum_\left\{i=1\right\}^N \, q_i \left(\boldsymbol\left\{r_i\right\}-\boldsymbol\left\{r\right\}\right),$

This expression is equivalent to the previous expression in the case of charge neutrality and $N = 2$. For two opposite charges, denoting the location of the positive charge of the pair as $\boldsymbol \left\{r_+\right\}$ and the location of the negative charge as $\boldsymbol \left\{r_-\right\}$ :

$\boldsymbol\left\{p\right\}\left(\boldsymbol\left\{r\right\}\right)$$=q_1\left(\boldsymbol\left\{r_1\right\}-\boldsymbol\left\{r\right\}\right)+q_2\left(\boldsymbol\left\{r_2\right\}-\boldsymbol\left\{r\right\}\right)$$= q\left(\boldsymbol\left\{r_+\right\}-\boldsymbol\left\{r\right\}\right)-q\left(\boldsymbol\left\{r_-\right\}-\boldsymbol\left\{r\right\}\right)$$=q \left(\boldsymbol\left\{r_+\right\} - \boldsymbol\left\{r_-\right\}\right)=q\boldsymbol d \ ,$

showing that the dipole moment vector is directed from the negative charge to the positive charge because the position vector of a point is directed outward from the origin to that point.

The dipole moment is most easily understood when the system has an overall neutral charge; for example, a pair of opposite charges, or a neutral conductor in a uniform electric field. For a system of charges with no net charge, visualized as an array of paired opposite charges, the relation for electric dipole moment is:

$\boldsymbol\left\{p\right\}\left(\boldsymbol\left\{r\right\}\right) = \sum_\left\{i=1\right\}^\left\{N\right\} \, \int q_i \left\left( \delta \left(\mathbf\left\{r_0\right\} - \left(\mathbf\left\{r\right\}_i + \boldsymbol\left\{d_i\right\}\right) \right)- \delta \left( \mathbf\left\{r_0\right\} - \mathbf\left\{r\right\}_i \right) \right\right)\, \left(\boldsymbol\left\{r_0\right\}-\boldsymbol\left\{r\right\}\right) \ d^3 \boldsymbol\left\{r_0\right\}$

= \sum_{i=1}^{N} \, q_i \left( \boldsymbol{r_i +d_i}-\boldsymbol{r} -(\boldsymbol{r_i }-\boldsymbol{r}) \right) = \sum_{i=1}^{N} q_i\boldsymbol{d}_i \, = \sum_{i=1}^{N} \boldsymbol{p}_i \ ,

which is the vector sum of the individual dipole moments of the neutral charge pairs. (Because of overall charge neutrality, the dipole moment is independent of the observer's position r.) Thus, the value of p is independent of the choice of reference point, provided the overall charge of the system is zero.

When discussing the dipole moment of a non-neutral system, such as the dipole moment of the proton, a dependence on the choice of reference point arises. In such cases it is conventional to choose the reference point to be the center of mass of the system or the center of charge, not some arbitrary origin.[1] This convention ensures that the dipole moment is an intrinsic property of the system.

## Potential and field of an electric dipole

An ideal dipole consists of two opposite charges with infinitesimal separation. The potential and field of such an ideal dipole are found next as a limiting case of an example of two opposite charges at non-zero separation.

Two closely spaced opposite charges have a potential of the form:

$\phi \left( \boldsymbol\left\{r\right\} \right)=\frac \left\{q\right\}\left\{4 \pi \varepsilon _0 | \boldsymbol\left\{ r\right\}- \boldsymbol\left\{r\right\}_+ |\right\} -\frac \left\{q\right\}\left\{4 \pi \varepsilon _0 | \boldsymbol\left\{ r\right\}- \boldsymbol\left\{r\right\}_- | \right\} \ ,$

with charge separation, d, defined as

$\boldsymbol d = \boldsymbol\left\{r\right\}_+ - \boldsymbol\left\{r\right\}_- \ ,$

The radius to the center of charge, R, and the unit vector in the direction of R are given by:

Taylor expansion in d/r (see multipole expansion and quadrupole) allows this potential to be expressed as a series.[2][3]

$\phi \left( \boldsymbol\left\{R\right\} \right)=\frac \left\{1\right\}\left\{4 \pi \varepsilon _0\right\} \frac \left\{q\boldsymbol \left\{d \cdot \hat\left\{R\right\}\right\}\right\}\left\{R^2\right\} + O\left(\frac\left\{d^2\right\}\left\{R^2\right\}\right) \approx \frac \left\{1\right\}\left\{4 \pi \varepsilon _0\right\} \frac \left\{\boldsymbol \left\{p\cdot \hat\left\{R\right\}\right\}\right\}\left\{R^2\right\} \ ,$

where higher order terms in the series are vanishing at large distances, R, compared to d.[4] Here, the electric dipole moment p is, as above:

$\boldsymbol p = q \boldsymbol d \ .$

The result for the dipole potential also can be expressed as:[5]

$\phi \left( \boldsymbol\left\{R\right\} \right)=- \boldsymbol \left\{p\cdot \nabla\right\}\frac \left\{1\right\}\left\{4 \pi \varepsilon _0 R\right\}\ ,$

which relates the dipole potential to that of a point charge. A key point is that the potential of the dipole falls off faster with distance R than that of the point charge.

The field of the dipole is the gradient of the potential, leading to:[5]

$\boldsymbol E = \frac \left\{3 \boldsymbol \left\{p \cdot \hat\left\{R\right\}\right\}\right\}\left\{4 \pi \varepsilon_0 R^3\right\} \boldsymbol \left\{ \hat\left\{R\right\}\right\}-\frac \left\{\boldsymbol\left\{p\right\}\right\}\left\{4 \pi \varepsilon_0 R^3\right\} \ .$

Thus, although two closely spaced opposite charges are not an ideal electric dipole (because their potential at close approach is not that of a dipole), at distances much larger than their separation, their dipole moment p appears directly in their potential and field.

As the two charges are brought closer together (d is made smaller), the dipole term in the multipole expansion based on the ratio d/R becomes the only significant term at ever closer distances R, and in the limit of infinitesimal separation the dipole term in this expansion is all that matters. As d is made infinitesimal, however, the dipole charge must be made to increase to hold p constant. This limiting process results in a "point dipole".

## Dipole moment density and polarization density

The dipole moment of an array of charge,

$\boldsymbol p = \sum_\left\{i=1\right\}^N \ q_i \boldsymbol \left\{d_i\right\} \ ,$

determines the degree of polarity of the array, but for a neutral array it is simply a vector property of the array with no directions about where the array happens to be located. 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 P(r) = p(r). 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 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:

$\boldsymbol\left\{D\right\} = \varepsilon _0 \boldsymbol\left\{E\right\} + \boldsymbol\left\{P\right\}\ ,$

where P is called the polarization density. In this formulation, the divergence of this equation yields:

$\nabla \cdot \boldsymbol\left\{D\right\} = \rho_f = \varepsilon _0 \nabla \cdot \boldsymbol\left\{E\right\} +\nabla \cdot \boldsymbol\left\{P\right\}\ ,$

and as the divergence term in E is the total charge, and ρf is "free charge", we are left with the relation:

$\nabla \cdot \boldsymbol\left\{P\right\} = -\rho_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 curl E = 0, which implies curl (D − P) = 0. Applying Helmholtz decomposition:[6]

$\boldsymbol\left\{ \left(D-P\right) = -\nabla \right\} \varphi \ ,$

for some scalar potential φ, and:

$\boldsymbol \left\{\nabla \cdot \left(D-P\right)\right\} =\varepsilon_0 \boldsymbol \left\{\nabla \cdot E\right\}=\rho_f +\rho_b = -\nabla ^2 \varphi \ .$

Suppose the charges are divided into free and bound, and the potential is divided into φ = φf + φ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.[7][8] In particular, when no free charge is present, one possible choice is P = ε0 E.

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 P(r) = p(r) restricted to the same model. For a smoothly varying dipole moment distribution p(r), the corresponding bound charge density is simply ∇·p(r) = −ρb. However, in the case of a p(r) that exhibits an abrupt step in dipole moment at a boundary between two regions, ∇·p(r) exhibits 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).[9] The potential at a position r is:[10][11]

$\phi \left( \boldsymbol\left\{r\right\} \right) = \frac \left\{1\right\}\left\{4 \pi \varepsilon_0\right\}\int \frac \left\{ \rho \left( \boldsymbol\left\{ r\right\}_0 \right)\right\} \left\{| \boldsymbol\left\{ r\right\}- \boldsymbol\left\{r\right\}_0 | \right\} d^3 \boldsymbol\left\{ r\right\}_0 \ + \frac \left\{1\right\}\left\{4 \pi \varepsilon_0\right\}\int \frac \left\{ \boldsymbol\left\{p\right\} \left( \boldsymbol\left\{ r\right\}_0 \right)\boldsymbol\left\{\cdot \left(r - r_0\right)\right\}\right\} \left\{| \boldsymbol\left\{ r\right\}- \boldsymbol\left\{r\right\}_0 |^3 \right\} d^3 \boldsymbol\left\{ r\right\}_0 ,$

where ρ(r) is the unpaired charge density, and p(r) is the dipole moment density.[12] Using an identity:

$\nabla_\left\{\boldsymbol \left\{r\right\}_0\right\} \frac \left\{1\right\}\left\{|\boldsymbol r - \boldsymbol\left\{r\right\}_0|\right\} = \frac \left\{\boldsymbol r - \boldsymbol\left\{r\right\}_0\right\}\left\{|\boldsymbol r - \boldsymbol\left\{r\right\}_0|^3\right\}$

the polarization integral can be transformed:

$\frac \left\{1\right\}\left\{4 \pi \varepsilon_0\right\}\int \frac \left\{ \boldsymbol\left\{p\right\} \left( \boldsymbol\left\{ r\right\}_0 \right)\boldsymbol\left\{\cdot \left(r - r_0\right)\right\}\right\} \left\{| \boldsymbol\left\{ r\right\}- \boldsymbol\left\{r\right\}_0 |^3 \right\} d^3 \boldsymbol\left\{ r\right\}_0 =\frac \left\{1\right\}\left\{4 \pi \varepsilon_0\right\}\int \boldsymbol\left\{p\right\} \left( \boldsymbol\left\{ r\right\}_0 \right)\boldsymbol\left\{\cdot \nabla\right\}_\left\{\boldsymbol \left\{r\right\}_0\right\} \frac \left\{1\right\}\left\{|\boldsymbol r - \boldsymbol\left\{r\right\}_0|\right\} d^3 \boldsymbol\left\{ r\right\}_0 ,$
$=\frac \left\{1\right\}\left\{4 \pi \varepsilon_0\right\}\int \boldsymbol\left\{\nabla_\left\{\boldsymbol \left\{r_0\right\}\right\}\cdot\right\} \left\left( \boldsymbol\left\{p\right\} \left( \boldsymbol\left\{ r\right\}_0 \right) \frac \left\{1\right\}\left\{|\boldsymbol r - \boldsymbol\left\{r\right\}_0|\right\} \right\right) d^3 \boldsymbol\left\{ r\right\}_0 -\frac \left\{1\right\}\left\{4 \pi \varepsilon_0\right\}\int \frac \left\{\boldsymbol\left\{\nabla_\left\{\boldsymbol \left\{r_0\right\}\right\}\cdot\right\} \boldsymbol\left\{p\right\} \left( \boldsymbol\left\{ r\right\}_0 \right)\right\}\left\{|\boldsymbol r - \boldsymbol\left\{r\right\}_0|\right\} d^3 \boldsymbol\left\{ r\right\}_0 ,$

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 \left( \boldsymbol\left\{r\right\} \right) = \frac \left\{1\right\}\left\{4 \pi \varepsilon_0\right\}\int \frac \left\{ \rho \left( \boldsymbol\left\{ r\right\}_0 \right)-\boldsymbol\left\{\nabla_\left\{\boldsymbol \left\{r_0\right\}\right\}\cdot\right\} \boldsymbol\left\{p\right\} \left( \boldsymbol\left\{ r\right\}_0 \right)\right\} \left\{| \boldsymbol\left\{ r\right\}- \boldsymbol\left\{r\right\}_0 | \right\} d^3 \boldsymbol\left\{ r\right\}_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_\left\{total\right\} \left(\boldsymbol\left\{ r\right\}_0\right)= \rho \left( \boldsymbol\left\{ r\right\}_0 \right)-\boldsymbol\left\{\nabla_\left\{\boldsymbol \left\{r_0\right\}\right\}\cdot\right\} \boldsymbol\left\{p\right\} \left( \boldsymbol\left\{ r\right\}_0 \right) \ ,$

showing that:

$-\boldsymbol\left\{\nabla_\left\{\boldsymbol \left\{r_0\right\}\right\}\cdot\right\} \boldsymbol\left\{p\right\} \left( \boldsymbol\left\{ r\right\}_0 \right) = \rho_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.[13][14] Using the relation:

$\nabla \cdot \boldsymbol\left\{D\right\} = \rho_f \ ,$

the polarization density is found to be:

$\boldsymbol\left\{P\left( r \right)\right\} = -\boldsymbol\left\{\nabla \cdot p_\left\{Dip\right\}\right\} -\boldsymbol\left\{\nabla \cdot p_\left\{Quad\right\}\right\} + \ldots \ ,$

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. 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.[15]

#### Surface charge

File:Dipole
A uniform array of identical dipoles is equivalent to a surface charge.

Above, discussion was deferred for the leading divergence term in the expression for the potential due to the dipoles. This term 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. The potential is:

$\phi \left( \boldsymbol\left\{r\right\} \right) =\frac \left\{1\right\}\left\{4 \pi \varepsilon_0\right\}\int \boldsymbol\left\{\nabla_\left\{\boldsymbol \left\{r_0\right\}\right\}\cdot\right\} \left\left( \boldsymbol\left\{p\right\} \left( \boldsymbol\left\{ r\right\}_0 \right) \frac \left\{1\right\}\left\{|\boldsymbol r - \boldsymbol\left\{r\right\}_0|\right\} \right\right) d^3 \boldsymbol\left\{ r\right\}_0$$-\frac \left\{1\right\}\left\{4 \pi \varepsilon_0\right\}\int \frac \left\{\boldsymbol\left\{\nabla_\left\{\boldsymbol \left\{r_0\right\}\right\}\cdot\right\} \boldsymbol\left\{p\right\} \left( \boldsymbol\left\{ r\right\}_0 \right)\right\}\left\{|\boldsymbol r - \boldsymbol\left\{r\right\}_0|\right\} d^3 \boldsymbol\left\{ r\right\}_0 \ .$

Using the divergence theorem, the divergence term transforms into the surface integral:

$\frac \left\{1\right\}\left\{4 \pi \varepsilon_0\right\}\int \boldsymbol\left\{\nabla_\left\{\boldsymbol \left\{r_0\right\}\right\}\cdot\right\} \left\left( \boldsymbol\left\{p\right\} \left( \boldsymbol\left\{ r\right\}_0 \right) \frac \left\{1\right\}\left\{|\boldsymbol r - \boldsymbol\left\{r\right\}_0|\right\} \right\right) d^3 \boldsymbol\left\{ r\right\}_0$
$=\frac \left\{1\right\}\left\{4 \pi \varepsilon_0\right\}\int \frac \left\{\boldsymbol\left\{p\right\} \left( \boldsymbol\left\{ r\right\}_0 \right)\boldsymbol\left\{\cdot \right\} d \boldsymbol \left\{A_0 \right\} \right\} \left\{|\boldsymbol r - \boldsymbol\left\{r\right\}_0|\right\} \ ,$

with dAo an element of surface area of the volume. In the event that p(r) is a constant, only the surface term survives:

$\phi \left( \boldsymbol\left\{r\right\} \right)$$=\frac \left\{1\right\}\left\{4 \pi \varepsilon_0\right\}\int \frac \left\{1\right\}\left\{|\boldsymbol r - \boldsymbol\left\{r\right\}_0|\right\}\ \boldsymbol\left\{p\right\} \cdot d\boldsymbol\left\{A_0\right\} \ ,$

with dAo 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 σ = dA, 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.[16] The field due to the surface charge is:

$\boldsymbol E \left( \boldsymbol\left\{r\right\} \right) =-\frac \left\{1\right\}\left\{4 \pi \varepsilon_0\right\} \boldsymbol\left\{\nabla\right\}_\left\{\boldsymbol \left\{r\right\}\right\}\int \frac \left\{1\right\}\left\{|\boldsymbol r - \boldsymbol\left\{r\right\}_0|\right\}\ \boldsymbol\left\{p\right\} \cdot d\boldsymbol\left\{A_0\right\} \ ,$

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 :[17]

$\boldsymbol E =-\frac \left\{\boldsymbol p\right\}\left\{3 \varepsilon_0\right\} \ .$

If we suppose the polarization of the dipoles was induced by an external field, the polarization field opposes the applied field and sometimes is called a depolarization field.[18][19] 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.[20]

In particular, if the electric susceptibility is introduced through the approximation:

$\boldsymbol\left\{p\left(r\right)\right\} = \varepsilon_0 \chi\left(\boldsymbol r \right) \boldsymbol \left\{E\left(r\right)\right\} \ ,$

then:

$\boldsymbol \left\{ \nabla \cdot p\left(r\right)\right\}=\boldsymbol \left\{ \nabla \cdot\right\} \left\left( \chi \boldsymbol\left\{ \left(r\right)\right\}\varepsilon_0 \boldsymbol \left\{E\left(r\right)\right\}\right\right) =-\rho_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\left\{\boldsymbol n\right\} \cdot \left\left( \chi \boldsymbol\left\{ \left(r_+\right)\right\}\boldsymbol \left\{E\left(r_+\right)\right\}-\chi \boldsymbol\left\{ \left(r_-\right)\right\}\boldsymbol \left\{E\left(r_-\right)\right\}\right\right) =\frac\left\{1\right\}\left\{A_n\right\} \int d \Omega_n \ \rho_b = 0 \ ,$

where An, Ωn indicate the area and volume of an elementary region straddling the boundary between the regions, and $\hat\left\{\boldsymbol n\right\}$ 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 p(r) = χ(r)E(r) necessarily includes the contribution of a surface charge.[21][22][23]

It may be noted that a physically more realistic modeling of p(r) would cause the dipole moment density to taper off continuously to zero at the boundary of the confining region, rather than making a sudden step to zero density. Then the surface charge becomes zero at the boundary, and the surface charge is replaced by the divergence of a continuously varying dipole-moment density.

#### Dielectric sphere in uniform external electric field

File:Dielectric
Electric field lines in dielectric sphere with greater susceptibility than its surroundings.[24]

The above general remarks about surface charge are made more concrete by considering the example of a dielectric sphere in a uniform electric field.[25][26] 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_\left\{\infty\right\} = - E_\left\{\infty\right\}z = -E_\left\{\infty\right\} r \cos \theta \ .$

The sphere is assumed to be described by a dielectric constant κ, that is, D = κε0E, and inside the sphere the potential satisfies Laplace's equation. Skipping a few details, the solution inside the sphere is:

$\phi_< = A r \cos \theta \ ,$

while outside the sphere:

$\phi_> = \left\left(Br + \frac \left\{C\right\}\left\{r^2\right\} \right \right) \cos \theta \ .$

At large distances, φ> → φ so B = -E. Continuity of potential and of the radial component of displacement D = κε0E determine the other two constants. Supposing the radius of the sphere is R,

As a consequence, the potential is:

$\phi_> = \left\left( \left\{-r\right\}+\frac \left\{\kappa-1\right\}\left\{\kappa+2\right\}\frac Template:R^3\left\{r^2\right\} \right\right)E_\left\{\infty\right\} \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:

$\boldsymbol p = 4 \pi \varepsilon_0 \left\left(\frac \left\{\kappa-1\right\}\left\{\kappa+2\right\}\left\{R^3\right\} \right\right) \boldsymbol\left\{E_\left\{\infty\right\}\right\} \ ,$

or, per unit volume:

$\frac \left\{\boldsymbol p\right\}\left\{V\right\} = \left\{3\right\}\varepsilon_0 \left\left(\frac \left\{\kappa-1\right\}\left\{\kappa+2\right\}\right\right) \boldsymbol\left\{E_\left\{\infty\right\}\right\} \ .$

The factor (κ-1)/(κ+2) is called the Clausius-Mossotti factor and shows that the induced polarization flips sign if κ < 1. 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_< = -\frac\left\{3\right\}\left\{\kappa +2\right\} E_\left\{\infty\right\}r \cos \theta \ ,$

leading to the field inside the sphere:

$\boldsymbol \left\{-\nabla\right\} \phi_< = \frac\left\{3\right\}\left\{\kappa +2\right\} \boldsymbol\left\{ E_\left\{\infty\right\}\right\} =\left\left( 1-\frac \left\{\kappa-1\right\}\left\{\kappa+2\right\} \right\right)\boldsymbol\left\{ E_\left\{\infty\right\}\right\} \ ,$

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 = \left\{3\right\}\varepsilon_0\frac \left\{\kappa-1\right\}\left\{\kappa+2\right\} E_\left\{\infty\right\} \cos \theta =\frac\left\{1\right\}\left\{V\right\} \boldsymbol\left\{ p \cdot \hat\left\{R\right\}\right\}\ .$

his example shows that, for this example, 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 a electric susceptibility or electrical permittivity.

A more complex model of the point charge array introduces an effective medium by averaging the microscopic charges;[19] for example, the averaging can arrange that only dipole fields play a role.[27][28] 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.[17][29] 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.[30] 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.[31][32][33]

## Dipole moments of fundamental particles

Much experimental work is continuing on measuring the electric dipole moments (EDM) of fundamental and composite particles, namely those of the neutron and electron. As EDMs violate both the Parity (P) and Time (T) symmetries, their values yield a mostly model-independent measure (assuming CPT symmetry is valid) of CP-violation in nature. Therefore, values for these EDMs place strong constraints upon the scale of CP-violation that extensions to the standard model of particle physics may allow.

Indeed, many theories are inconsistent with the current limits and have effectively been ruled out, and established theory permits a much larger value than these limits, leading to the strong CP problem and prompting searches for new particles such as the axion.

Current generations of experiments are designed to be sensitive to the supersymmetry range of EDMs, providing complementary experiments to those done at the LHC.

## References and in-line notes

1. ^ Christopher J. Cramer (2004). Essentials of computational chemistry (2 ed.). Wiley. p. 307. ISBN 0-470-09182-7.
2. ^ David E Dugdale (1993). Essentials of Electromagnetism. Springer. pp. 80–81. ISBN 1-56396-253-5.
3. ^ Kikuji Hirose, Tomoya Ono, Yoshitaka Fujimoto (2005). First-principles calculations in real-space formalism. Imperial College Press. p. 18. ISBN 1-86094-512-0.
4. ^ Each succeeding term provides a more detailed view of the distribution of charge, and falls off more rapidly with distance. For example, the quadrupole moment is the basis for the next term: $Q_\left\{ij\right\} = \int d^3 \boldsymbol\left\{r_0\right\} \left\left( 3x_i x_j -r_0^2 \delta_\left\{ij\right\} \right\right) \rho\left( \boldsymbol\left\{r_0\right\}\right) \ ,$ with ro = (x1, x2, x3). See HW Wyld (1999). Mathematical Methods for Physics. Westview Press. p. 106. ISBN 0-7382-0125-1.
5. ^ a b BB Laud (1987). Electromagnetics (2 ed.). New Age International. p. 25. ISBN 0-85226-499-2.
6. ^ Jie-Zhi Wu, Hui-Yang Ma, Ming-De Zhou (2006). "§2.3.1 Functionally Orthogonal Decomposition". Vorticity and vortex dynamics. Springer. pp. 36 ff. ISBN 3-540-29027-3.
7. ^ For example, one could place the boundary around the bound charges at infinity. Then φb falls off with distance from the bound charges. If an external field is present, and zero free charge, the field can be accounted for in the contribution of φf, which would arrange to satisfy the boundary conditions and Laplace's equation2φf = 0.
8. ^ In principle, one could add the same arbitrary curl to both D and P, which would cancel out of the difference D − P. However, assuming D and P originate in in a simple division of charges into free and bound, they are at bottom electric fields and so have zero curl.
9. ^ This medium can be seen as an idealization growing from the multipole expansion of the potential of an arbitrarily complex charge distribution, truncation of the expansion, and the forcing of the truncated form to apply everywhere. The result is a hypothetical medium. See Jack Vanderlinde (2004). "§7.1 The electric field due to a polarized dielectric". Classical Electromagnetic Theory. Springer. ISBN 1-4020-2699-4.
10. ^ Uwe Krey, Anthony Owen (2007). Basic Theoretical Physics: A Concise Overview. Springer. pp. 138–143. ISBN 3-540-36804-3.
11. ^ T Tsang (1997). Classical Electrodynamics. World Scientific. p. 59. ISBN 981-02-3041-9.
12. ^ For example, for a system of ideal dipoles with dipole moment p confined within some closed surface, the dipole density p(r) is equal to p inside the surface, but is zero outside. That is, the dipole density includes a Heaviside step function locating the dipoles inside the surface.
13. ^ George E Owen (2003). Introduction to Electromagnetic Theory (republication of the 1963 Allyn & Bacon ed.). Courier Dover Publications. p. 80. ISBN 0-486-42830-3.
14. ^ Pierre-François Brevet (1997). Surface second harmonic generation. Presses polytechniques et universitaires romandes. p. 24. ISBN 2-88074-345-1.
15. ^ See Daniel A. Jelski, Thomas F. George (1999). Computational studies of new materials. World Scientific. p. 219. ISBN 981-02-3325-6.  and EM Purcell & CR Pennypacker (1973). "Scattering and Absorption of Light by Nonspherical Dielectric Grains". Astrophysical Journal 186: 705–714. doi:10.1086/152538.
16. ^ A brute force evaluation of the integral can be done using a multipole expansion: $\frac\left\{1\right\}\left\{|\boldsymbol\left\{r-r_0\right\}|\right\}$ = $\sum_\left\{\ell,\ m \right\} \frac\left\{4\pi\right\}\left\{2 \ell +1\right\}$·$\frac \left\{1\right\}\left\{r\right\} \left\left(\left\{\frac \left\{r_0\right\}\left\{r\right\}\right\}\right\right)^\left\{\ell\right\}$·$\ \left\{Y^*\right\}_\left\{\ell\right\}^m \left(\theta_0 , \ \phi_0\right) Y_\left\{\ell\right\}^m \left(\theta, \ \phi\right)$. See HW Wyld (1999). Mathematical Methods for Physics. Westview Press. p. 104. ISBN 0-7382-0125-1.
17. ^ a b H. Ibach, Hans Lüth (2003). Solid-state Physics: an introduction to principles of materials science (3 ed.). Springer. p. 361. ISBN 3-540-43870-X.
18. ^ Yasuaki Masumoto, Toshihide Takagahara (2002). Semiconductor quantum dots: physics, spectroscopy, and applications. Springer. p. 72. ISBN 3-540-42805-4.
19. ^ a b Yutaka Toyozawa (2003). Optical processes in solids. Cambridge University Press. p. 96. ISBN 0-521-55605-8.
20. ^ For example, a droplet in a surrounding medium experiences a higher or a lower internal field depending upon whether the medium has a higher or a lower dielectric constant than that of the droplet. See Paul S. Drzaic (1995). Liquid crystal dispersions. World Scientific. p. 246. ISBN 981-02-1745-5.
21. ^ Wai-Kai Chen (2005). The electrical engineering handbook. Academic Press. p. 502. ISBN 0-12-170960-4.
22. ^ Julius Adams Stratton (2007). Electromagnetic theory (reprint of 1941 ed.). Wiley-IEEE. p. 184. ISBN 0-470-13153-5.
23. ^ Edward J. Rothwell, Michael J. Cloud (2001). Electromagnetics. CRC Press. p. 68. ISBN 0-8493-1397-X.
24. ^ Based upon equations from Andrew Gray (1888). The theory and practice of absolute measurements in electricity and magnetism. Macmillan & Co.. pp. 126–127. , which refers to papers by Sir W. Thomson.
25. ^ HW Wyld (1999). Mathematical Methods for Physics (2 ed.). Westview Press. pp. 233 ff. ISBN 0-7382-0125-1.
26. ^ Julius Adams Stratton (2007). Electromagnetic theory (Wiley-IEEE reissue ed.). Piscataway, NJ: IEEE Press. p. 205 ff. ISBN 0-470-13153-5.
27. ^ John E Swipe & RW Boyd (2002). "Nanocomposite materials for nonlinear optics based upon local field effects". In Vladimir M. Shalaev. Optical properties of nanostructured random media. Springer. p. 3. ISBN 3-540-42031-2.
28. ^ Emil Wolf (1977). Progress in Optics. Elsevier. p. 288. ISBN 0-7204-1515-2.
29. ^ Mark Fox (2006). Optical Properties of Solids. Oxford University Press. p. 39. ISBN 0-19-850612-0.
30. ^ Lev Kantorovich (2004). "§8.2.1 The local field". Quantum theory of the solid state. Springer. p. 426. ISBN 1-4020-2153-4.
31. ^ Pierre Meystre (2001). Atom Optics. Springer. p. 5. ISBN 0-387-95274-8.
32. ^ Bruce T Draine (2001). "The discrete dipole approximation for light scattering by irregular targets". In Michael I. Mishchenko. Light scattering by nonspherical particles. Academic Press. p. 132. ISBN 0-12-498660-9.
33. ^ MA Yurkin & AG Hoekstra (2007). "The discrete dipole approximation: an overview and recent developments". Journal of Quantitative Spectroscopy and Radiative Transfer 106: 558–589. doi:10.1016/j.jqsrt.2007.01.034.

Melvin Schwartz (1987). "Electrical dipole moment". Principles of Electrodynamics (reprint of 1972 ed.). Courier Dover Publications. p. 49ff. ISBN 0-486-65493-1.