# Poynting vector

Dipole radiation of a dipole vertically in the page showing electric field strength (color) and Poynting vector (arrows) in the plane of the page.

In physics, the Poynting vector (or Umov–Poynting vector) represents the directional energy flux (the energy transfer per unit area per unit time) or power flow of an electromagnetic field. The SI unit of the Poynting vector is the watt per square metre (W/m2); kg/s3 in base SI units. It is named after its discoverer John Henry Poynting who first derived it in 1884.[1]: 132  Nikolay Umov is also credited with formulating the concept.[2] Oliver Heaviside also discovered it independently in the more general form that recognises the freedom of adding the curl of an arbitrary vector field to the definition.[3] The Poynting vector is used throughout electromagnetics in conjunction with Poynting's theorem, the continuity equation expressing conservation of electromagnetic energy, to calculate the power flow in electromagnetic fields.

## Definition

In Poynting's original paper and in most textbooks, the Poynting vector ${\displaystyle \mathbf {S} }$ is defined as the cross product[4][5][6]

${\displaystyle \mathbf {S} =\mathbf {E} \times \mathbf {H} ,}$
where bold letters represent vectors and

This expression is often called the Abraham form and is the most widely used.[7] The Poynting vector is usually denoted by S or N.

In simple terms, the Poynting vector S depicts the direction and rate of transfer of energy, that is power, due to electromagnetic fields in a region of space that may or may not be empty. More rigorously, it is the quantity that must be used to make Poynting's theorem valid. Poynting's theorem essentially says that the difference between the electromagnetic energy entering a region and the electromagnetic energy leaving a region must equal the energy converted or dissipated in that region, that is, turned into a different form of energy (often heat). So if one accepts the validity of the Poynting vector description of electromagnetic energy transfer, then Poynting's theorem is simply a statement of the conservation of energy.

If electromagnetic energy is not gained from or lost to other forms of energy within some region (e.g., mechanical energy, or heat), then electromagnetic energy is locally conserved within that region, yielding a continuity equation as a special case of Poynting's theorem:

${\displaystyle \nabla \cdot \mathbf {S} =-{\frac {\partial u}{\partial t}}}$
where ${\displaystyle u}$ is the energy density of the electromagnetic field. This frequent condition holds in the following simple example in which the Poynting vector is calculated and seen to be consistent with the usual computation of power in an electric circuit.

## Example: Power flow in a coaxial cable

Although problems in electromagnetics with arbitrary geometries are notoriously difficult to solve, we can find a relatively simple solution in the case of power transmission through a section of coaxial cable analyzed in cylindrical coordinates as depicted in the accompanying diagram. We can take advantage of the model's symmetry: no dependence on θ (circular symmetry) nor on Z (position along the cable). The model (and solution) can be considered simply as a DC circuit with no time dependence, but the following solution applies equally well to the transmission of radio frequency power, as long as we are considering an instant of time (during which the voltage and current don't change), and over a sufficiently short segment of cable (much smaller than a wavelength, so that these quantities are not dependent on Z). The coaxial cable is specified as having an inner conductor of radius R1 and an outer conductor whose inner radius is R2 (its thickness beyond R2 doesn't affect the following analysis). In between R1 and R2 the cable contains an ideal dielectric material of relative permittivity εr and we assume conductors that are non-magnetic (so μ = μ0) and lossless (perfect conductors), all of which are good approximations to real-world coaxial cable in typical situations.

Illustration of electromagnetic power flow inside a coaxial cable according to the Poynting vector S, calculated using the electric field E (due to the voltage V) and the magnetic field H (due to current I).
DC power transmission through a coaxial cable showing relative strength of electric (${\displaystyle E_{r}}$) and magnetic (${\displaystyle H_{\theta }}$) fields and resulting Poynting vector (${\displaystyle S_{z}=E_{r}\cdot H_{\theta }}$) at a radius r from the center of the coaxial cable. The broken magenta line shows the cumulative power transmission within radius r, half of which flows inside the geometric mean of R1 and R2.

The center conductor is held at voltage V and draws a current I toward the right, so we expect a total power flow of P = V · I according to basic laws of electricity. By evaluating the Poynting vector, however, we are able to identify the profile of power flow in terms of the electric and magnetic fields inside the coaxial cable. The electric fields are of course zero inside of each conductor, but in between the conductors (${\displaystyle R_{1}) symmetry dictates that they are strictly in the radial direction and it can be shown (using Gauss's law) that they must obey the following form:

${\displaystyle E_{r}(r)={\frac {W}{r}}}$
W can be evaluated by integrating the electric field from ${\displaystyle r=R_{2}}$ to ${\displaystyle R_{1}}$ which must be the negative of the voltage V:
${\displaystyle -V=\int _{R_{2}}^{R_{1}}{\frac {W}{r}}dr=-W\ln \left({\frac {R_{2}}{R_{1}}}\right)}$
so that:
${\displaystyle W={\frac {V}{\ln(R_{2}/R_{1})}}}$

The magnetic field, again by symmetry, can only be non-zero in the θ direction, that is, a vector field looping around the center conductor at every radius between R1 and R2. Inside the conductors themselves the magnetic field may or may not be zero, but this is of no concern since the Poynting vector in these regions is zero due to the electric field's being zero. Outside the entire coaxial cable, the magnetic field is identically zero since paths in this region enclose a net current of zero (+I in the center conductor and −I in the outer conductor), and again the electric field is zero there anyway. Using Ampère's law in the region from R1 to R2, which encloses the current +I in the center conductor but with no contribution from the current in the outer conductor, we find at radius r:

{\displaystyle {\begin{aligned}I=\oint _{C}\mathbf {H} \cdot ds&=2\pi rH_{\theta }(r)\\H_{\theta }(r)&={\frac {I}{2\pi r}}\end{aligned}}}
Now, from an electric field in the radial direction, and a tangential magnetic field, the Poynting vector, given by the cross-product of these, is only non-zero in the Z direction, along the direction of the coaxial cable itself, as we would expect. Again only a function of r, we can evaluate S(r):
${\displaystyle S_{z}(r)=E_{r}(r)H_{\theta }(r)={\frac {W}{r}}{\frac {I}{2\pi r}}={\frac {W\,I}{2\pi r^{2}}}}$
where W is given above in terms of the center conductor voltage V. The total power flowing down the coaxial cable can be computed by integrating over the entire cross section A of the cable in between the conductors:
{\displaystyle {\begin{aligned}P_{\text{tot}}&=\iint _{\mathbf {A} }S_{z}(r,\theta )\,dA=\int _{R_{2}}^{R_{1}}2\pi rdrS_{z}(r)\\&=\int _{R_{2}}^{R_{1}}{\frac {W\,I}{r}}dr=W\,I\,\ln \left({\frac {R_{2}}{R_{1}}}\right).\end{aligned}}}

Substituting the earlier solution for the constant W we find:

${\displaystyle P_{\mathrm {tot} }=I\ln \left({\frac {R_{2}}{R_{1}}}\right){\frac {V}{\ln(R_{2}/R_{1})}}=V\,I}$
that is, the power given by integrating the Poynting vector over a cross section of the coaxial cable is exactly equal to the product of voltage and current as one would have computed for the power delivered using basic laws of electricity.

## Other forms

In the "microscopic" version of Maxwell's equations, this definition must be replaced by a definition in terms of the electric field E and the magnetic flux density B (described later in the article).

It is also possible to combine the electric displacement field D with the magnetic flux B to get the Minkowski form of the Poynting vector, or use D and H to construct yet another version. The choice has been controversial: Pfeifer et al.[8] summarize and to a certain extent resolve the century-long dispute between proponents of the Abraham and Minkowski forms (see Abraham–Minkowski controversy).

The Poynting vector represents the particular case of an energy flux vector for electromagnetic energy. However, any type of energy has its direction of movement in space, as well as its density, so energy flux vectors can be defined for other types of energy as well, e.g., for mechanical energy. The Umov–Poynting vector[9] discovered by Nikolay Umov in 1874 describes energy flux in liquid and elastic media in a completely generalized view.

## Interpretation

The Poynting vector appears in Poynting's theorem (see that article for the derivation), an energy-conservation law:

${\displaystyle {\frac {\partial u}{\partial t}}=-\mathbf {\nabla } \cdot \mathbf {S} -\mathbf {J_{\mathrm {f} }} \cdot \mathbf {E} ,}$

where Jf is the current density of free charges and u is the electromagnetic energy density for linear, nondispersive materials, given by

${\displaystyle u={\frac {1}{2}}\!\left(\mathbf {E} \cdot \mathbf {D} +\mathbf {B} \cdot \mathbf {H} \right)\!,}$

where

• E is the electric field;
• D is the electric displacement field;
• B is the magnetic flux density;
• H is the magnetizing field.[10]: 258–260

The first term in the right-hand side represents the electromagnetic energy flow into a small volume, while the second term subtracts the work done by the field on free electrical currents, which thereby exits from electromagnetic energy as dissipation, heat, etc. In this definition, bound electrical currents are not included in this term and instead contribute to S and u.

For linear, nondispersive and isotropic (for simplicity) materials, the constitutive relations can be written as

${\displaystyle \mathbf {D} =\varepsilon \mathbf {E} ,\quad \mathbf {B} =\mu \mathbf {H} ,}$

where

Here ε and μ are scalar, real-valued constants independent of position, direction, and frequency.

In principle, this limits Poynting's theorem in this form to fields in vacuum and nondispersive[clarification needed] linear materials. A generalization to dispersive materials is possible under certain circumstances at the cost of additional terms.[10]: 262–264

One consequence of the Poynting formula is that for the electromagnetic field to do work, both magnetic and electric fields must be present. The magnetic field alone or the electric field alone cannot do any work.[11]

## Plane waves

In a propagating electromagnetic plane wave in an isotropic lossless medium, the instantaneous Poynting vector always points in the direction of propagation while rapidly oscillating in magnitude. This can be simply seen given that in a plane wave, the magnitude of the magnetic field H(r,t) is given by the magnitude of the electric field vector E(r,t) divided by η, the intrinsic impedance of the transmission medium:

${\displaystyle |\mathbf {H} |={\frac {|\mathbf {E} |}{\eta }}}$

where |A| represents the vector norm of A. Since E and H are at right angles to each other, the magnitude of their cross product is the product of their magnitudes. Without loss of generality let us take X to be the direction of the electric field and Y to be the direction of the magnetic field. The instantaneous Poynting vector, given by the cross product of E and H will then be in the positive Z direction:

${\displaystyle {\mathsf {S_{z}}}={\mathsf {E_{x}}}\cdot {\mathsf {H_{y}}}={\frac {\left|{\mathsf {E_{x}}}\right|^{2}}{\eta }}}$.

Finding the time-averaged power in the plane wave then requires averaging over time periods large compared to the frequency:

${\displaystyle \left\langle {\mathsf {S_{z}}}\right\rangle ={\frac {\left\langle \left|{\mathsf {E_{x}}}\right|^{2}\right\rangle }{\eta }}={\frac {\mathsf {E_{rms}^{2}}}{\eta }}}$

where Erms is the root mean square electric field amplitude. In the important case that E(t) is sinusoidally varying at some frequency with peak amplitude Epeak, its rms voltage is given by ${\displaystyle {\mathsf {E_{peak}}}/{\sqrt {2}}}$, with the average Poynting vector then given by:

${\displaystyle \left\langle {\mathsf {S_{z}}}\right\rangle ={\frac {\mathsf {E_{peak}^{2}}}{2\eta }}}$

This is the most common form for the energy flux of a plane wave, since sinusoidal field amplitudes are most often expressed in terms of their peak values, and complicated problems are typically solved considering only one frequency at a time. However, the expression using Erms is totally general, applying, for instance, in the case of noise whose RMS amplitude can be measured but where the "peak" amplitude is meaningless. In free space the intrinsic impedance η is simply given by the impedance of free space η0 377 Ω. In non-magnetic dielectrics (such as all transparent materials at optical frequencies) with a specified dielectric constant εr, or in optics with a material whose refractive index ${\displaystyle {\mathsf {n}}={\sqrt {\epsilon _{r}}}}$, the intrinsic impedance is found as:

${\displaystyle \eta ={\frac {\eta _{0}}{\sqrt {\epsilon _{r}}}}}$.

In optics, the value of radiated flux crossing a surface, thus the average Poynting vector component in the direction normal to that surface, is technically known as the irradiance, more often simply referred to as the intensity (a somewhat ambiguous term).

## Formulation in terms of microscopic fields

The "microscopic" (differential) version of Maxwell's equations admits only the fundamental fields E and B, without a built-in model of material media. Only the vacuum permittivity and permeability are used, and there is no D or H. When this model is used, the Poynting vector is defined as

${\displaystyle \mathbf {S} ={\frac {1}{\mu _{0}}}\mathbf {E} \times \mathbf {B} ,}$

where

This is actually the general expression of the Poynting vector[dubious ].[12] The corresponding form of Poynting's theorem is

${\displaystyle {\frac {\partial u}{\partial t}}=-\nabla \cdot \mathbf {S} -\mathbf {J} \cdot \mathbf {E} ,}$

where J is the total current density and the energy density u is given by

${\displaystyle u={\frac {1}{2}}\!\left(\varepsilon _{0}|\mathbf {E} |^{2}+{\frac {1}{\mu _{0}}}|\mathbf {B} |^{2}\right)\!,}$

where ε0 is the vacuum permittivity. It can be derived directly from Maxwell's equations in terms of total charge and current and the Lorentz force law only.

The two alternative definitions of the Poynting vector are equal in vacuum or in non-magnetic materials, where B = μ0H. In all other cases, they differ in that S = (1/μ0) E × B and the corresponding u are purely radiative, since the dissipation term JE covers the total current, while the E × H definition has contributions from bound currents which are then excluded from the dissipation term.[13]

Since only the microscopic fields E and B occur in the derivation of S = (1/μ0) E × B and the energy density, assumptions about any material present are avoided. The Poynting vector and theorem and expression for energy density are universally valid in vacuum and all materials.[13]

## Time-averaged Poynting vector

The above form for the Poynting vector represents the instantaneous power flow due to instantaneous electric and magnetic fields. More commonly, problems in electromagnetics are solved in terms of sinusoidally varying fields at a specified frequency. The results can then be applied more generally, for instance, by representing incoherent radiation as a superposition of such waves at different frequencies and with fluctuating amplitudes.

We would thus not be considering the instantaneous E(t) and H(t) used above, but rather a complex (vector) amplitude for each which describes a coherent wave's phase (as well as amplitude) using phasor notation. These complex amplitude vectors are not functions of time, as they are understood to refer to oscillations over all time. A phasor such as Em is understood to signify a sinusoidally varying field whose instantaneous amplitude E(t) follows the real part of Emejωt where ω is the (radian) frequency of the sinusoidal wave being considered.

In the time domain, it will be seen that the instantaneous power flow will be fluctuating at a frequency of 2ω. But what is normally of interest is the average power flow in which those fluctuations are not considered. In the math below, this is accomplished by integrating over a full cycle T = 2π / ω. The following quantity, still referred to as a "Poynting vector", is expressed directly in terms of the phasors as:

${\displaystyle \mathbf {S} _{\mathrm {m} }={\tfrac {1}{2}}\mathbf {E} _{\mathrm {m} }\times \mathbf {H} _{\mathrm {m} }^{*},}$

where denotes the complex conjugate. The time-averaged power flow (according to the instantaneous Poynting vector averaged over a full cycle, for instance) is then given by the real part of Sm. The imaginary part is usually ignored, however, it signifies "reactive power" such as the interference due to a standing wave or the near field of an antenna. In a single electromagnetic plane wave (rather than a standing wave which can be described as two such waves travelling in opposite directions), E and H are exactly in phase, so Sm is simply a real number according to the above definition.

The equivalence of Re(Sm) to the time-average of the instantaneous Poynting vector S can be shown as follows.

{\displaystyle {\begin{aligned}\mathbf {S} (t)&=\mathbf {E} (t)\times \mathbf {H} (t)\\&=\operatorname {Re} \!\left(\mathbf {E} _{\mathrm {m} }e^{j\omega t}\right)\times \operatorname {Re} \!\left(\mathbf {H} _{\mathrm {m} }e^{j\omega t}\right)\\&={\tfrac {1}{2}}\!\left(\mathbf {E} _{\mathrm {m} }e^{j\omega t}+\mathbf {E} _{\mathrm {m} }^{*}e^{-j\omega t}\right)\times {\tfrac {1}{2}}\!\left(\mathbf {H} _{\mathrm {m} }e^{j\omega t}+\mathbf {H} _{\mathrm {m} }^{*}e^{-j\omega t}\right)\\&={\tfrac {1}{4}}\!\left(\mathbf {E} _{\mathrm {m} }\times \mathbf {H} _{\mathrm {m} }^{*}+\mathbf {E} _{\mathrm {m} }^{*}\times \mathbf {H} _{\mathrm {m} }+\mathbf {E} _{\mathrm {m} }\times \mathbf {H} _{\mathrm {m} }e^{2j\omega t}+\mathbf {E} _{\mathrm {m} }^{*}\times \mathbf {H} _{\mathrm {m} }^{*}e^{-2j\omega t}\right)\\&={\tfrac {1}{2}}\operatorname {Re} \!\left(\mathbf {E} _{\mathrm {m} }\times \mathbf {H} _{\mathrm {m} }^{*}\right)+{\tfrac {1}{2}}\operatorname {Re} \!\left(\mathbf {E} _{\mathrm {m} }\times \mathbf {H} _{\mathrm {m} }e^{2j\omega t}\right)\!.\end{aligned}}}

The average of the instantaneous Poynting vector S over time is given by:

${\displaystyle \langle \mathbf {S} \rangle ={\frac {1}{T}}\int _{0}^{T}\mathbf {S} (t)\,dt={\frac {1}{T}}\int _{0}^{T}\!\left[{\tfrac {1}{2}}\operatorname {Re} \!\left(\mathbf {E} _{\mathrm {m} }\times \mathbf {H} _{\mathrm {m} }^{*}\right)+{\tfrac {1}{2}}\operatorname {Re} \!\left({\mathbf {E} _{\mathrm {m} }}\times {\mathbf {H} _{\mathrm {m} }}e^{2j\omega t}\right)\right]dt.}$

The second term is the double-frequency component having an average value of zero, so we find:

${\displaystyle \langle \mathbf {S} \rangle =\operatorname {Re} \!\left({\tfrac {1}{2}}{\mathbf {E} _{\mathrm {m} }}\times \mathbf {H} _{\mathrm {m} }^{*}\right)=\operatorname {Re} \!\left(\mathbf {S} _{\mathrm {m} }\right)}$

According to some conventions, the factor of 1/2 in the above definition may be left out. Multiplication by 1/2 is required to properly describe the power flow since the magnitudes of Em and Hm refer to the peak fields of the oscillating quantities. If rather the fields are described in terms of their root mean square (RMS) values (which are each smaller by the factor ${\displaystyle {\sqrt {2}}/2}$), then the correct average power flow is obtained without multiplication by 1/2.

## Resistive dissipation

If a conductor has significant resistance, then, near the surface of that conductor, the Poynting vector would be tilted toward and impinge upon the conductor. Once the Poynting vector enters the conductor, it is bent to a direction that is almost perpendicular to the surface.[14]: 61  This is a consequence of Snell's law and the very slow speed of light inside a conductor. The definition and computation of the speed of light in a conductor can be given.[15]: 402  Inside the conductor, the Poynting vector represents energy flow from the electromagnetic field into the wire, producing resistive Joule heating in the wire. For a derivation that starts with Snell's law see Reitz page 454.[16]: 454

The density of the linear momentum of the electromagnetic field is S/c2 where S is the magnitude of the Poynting vector and c is the speed of light in free space. The radiation pressure exerted by an electromagnetic wave on the surface of a target is given by

${\displaystyle P_{\mathrm {rad} }={\frac {\langle S\rangle }{\mathrm {c} }}.}$

## Uniqueness of the Poynting vector

The Poynting vector occurs in Poynting's theorem only through its divergence ∇ ⋅ S, that is, it is only required that the surface integral of the Poynting vector around a closed surface describe the net flow of electromagnetic energy into or out of the enclosed volume. This means that adding a solenoidal vector field (one with zero divergence) to S will result in another field that satisfies this required property of a Poynting vector field according to Poynting's theorem. Since the divergence of any curl is zero, one can add the curl of any vector field to the Poynting vector and the resulting vector field S′ will still satisfy Poynting's theorem.

However even though the Poynting vector was originally formulated only for the sake of Poynting's theorem in which only its divergence appears, it turns out that the above choice of its form is unique.[10]: 258–260, 605–612  The following section gives an example which illustrates why it is not acceptable to add an arbitrary solenoidal field to E × H.

## Static fields

Poynting vector in a static field, where E is the electric field, H the magnetic field, and S the Poynting vector.

The consideration of the Poynting vector in static fields shows the relativistic nature of the Maxwell equations and allows a better understanding of the magnetic component of the Lorentz force, q(v × B). To illustrate, the accompanying picture is considered, which describes the Poynting vector in a cylindrical capacitor, which is located in an H field (pointing into the page) generated by a permanent magnet. Although there are only static electric and magnetic fields, the calculation of the Poynting vector produces a clockwise circular flow of electromagnetic energy, with no beginning or end.

While the circulating energy flow may seem unphysical, its existence is necessary to maintain conservation of angular momentum. The momentum of an electromagnetic wave in free space is equal to its power divided by c, the speed of light. Therefore the circular flow of electromagnetic energy implies an angular momentum.[17] If one were to connect a wire between the two plates of the charged capacitor, then there would be a Lorentz force on that wire while the capacitor is discharging due to the discharge current and the crossed magnetic field; that force would be tangential to the central axis and thus add angular momentum to the system. That angular momentum would match the "hidden" angular momentum, revealed by the Poynting vector, circulating before the capacitor was discharged.

## References

1. ^ Stratton, Julius Adams (1941). Electromagnetic Theory (1st ed.). New York: McGraw-Hill. ISBN 978-0-470-13153-4.
2. ^ "Пойнтинга вектор". Физическая энциклопедия (in Russian). Retrieved 2022-02-21.
3. ^ Nahin, Paul J. (2002). Oliver Heaviside: The Life, Work, and Times of an Electrical Genius of the Victorian Age. p. 131. ISBN 9780801869099.
4. ^ Poynting, John Henry (1884). "On the Transfer of Energy in the Electromagnetic Field". Philosophical Transactions of the Royal Society of London. 175: 343–361. doi:10.1098/rstl.1884.0016.
5. ^ Grant, Ian S.; Phillips, William R. (1990). Electromagnetism (2nd ed.). New York: John Wiley & Sons. ISBN 978-0-471-92712-9.
6. ^ Griffiths, David J. (2012). Introduction to Electrodynamics (3rd ed.). Boston: Addison-Wesley. ISBN 978-0-321-85656-2.
7. ^ Kinsler, Paul; Favaro, Alberto; McCall, Martin W. (2009). "Four Poynting Theorems". European Journal of Physics. 30 (5): 983. arXiv:0908.1721. Bibcode:2009EJPh...30..983K. doi:10.1088/0143-0807/30/5/007. S2CID 118508886.
8. ^ Pfeifer, Robert N. C.; Nieminen, Timo A.; Heckenberg, Norman R.; Rubinsztein-Dunlop, Halina (2007). "Momentum of an Electromagnetic Wave in Dielectric Media". Reviews of Modern Physics. 79 (4): 1197. arXiv:0710.0461. Bibcode:2007RvMP...79.1197P. doi:10.1103/RevModPhys.79.1197.
9. ^ Umov, Nikolay Alekseevich (1874). "Ein Theorem über die Wechselwirkungen in Endlichen Entfernungen". Zeitschrift für Mathematik und Physik. 19: 97–114.
10. ^ a b c d Jackson, John David (1998). Classical Electrodynamics (3rd ed.). New York: John Wiley & Sons. ISBN 978-0-471-30932-1.
11. ^ "K. McDonald's Physics Examples - Railgun" (PDF). puhep1.princeton.edu. Retrieved 2021-02-14.
12. ^ Zangwill, Andrew (2013). Modern Electrodynamics. Cambridge University Press. p. 508. ISBN 9780521896979.
13. ^ a b Richter, Felix; Florian, Matthias; Henneberger, Klaus (2008). "Poynting's Theorem and Energy Conservation in the Propagation of Light in Bounded Media". EPL. 81 (6): 67005. arXiv:0710.0515. Bibcode:2008EL.....8167005R. doi:10.1209/0295-5075/81/67005. S2CID 119243693.
14. ^ Harrington, Roger F. (2001). Time-Harmonic Electromagnetic Fields (2nd ed.). McGraw-Hill. ISBN 978-0-471-20806-8.
15. ^ Hayt, William (2011). Engineering Electromagnetics (4th ed.). New York: McGraw-Hill. ISBN 978-0-07-338066-7.
16. ^ Reitz, John R.; Milford, Frederick J.; Christy, Robert W. (2008). Foundations of Electromagnetic Theory (4th ed.). Boston: Addison-Wesley. ISBN 978-0-321-58174-7.
17. ^ Feynman, Richard Phillips (2011). The Feynman Lectures on Physics. Vol. II: Mainly Electromagnetism and Matter (The New Millennium ed.). New York: Basic Books. ISBN 978-0-465-02494-0.