Maxwell's equations

For thermodynamic relations, see Maxwell relations. For the history of the equations, see History of Maxwell's equations. For a general description of electromagnetism, see Electromagnetism.

Maxwell's equations (mid-left) as featured on a monument in front of Warsaw University's Center of New Technologies

Maxwell's equations are a set of partial differential equations that, together with the Lorentz force law, form the foundation of classical electromagnetism, classical optics, and electric circuits. They underpin all electric, optical and radio technologies, including power generation, electric motors, wireless communication, cameras, televisions, computers etc. Maxwell's equations describe how electric and magnetic fields are generated by charges, currents, and changes of each other. One important consequence of the equations is that they demonstrate how fluctuating electric and magnetic fields propagate at the speed of light. Known as electromagnetic radiation, these waves may occur at various wavelengths to produce a spectrum from radio waves to γ-rays. The equations are named after the physicist and mathematician James Clerk Maxwell, who between 1861 and 1862 published an early form of the equations, and first proposed that light is an electromagnetic phenomenon.

The equations have two major variants. The microscopic Maxwell equations have universal applicability, but are unwieldy for common calculations. They relate the electric and magnetic fields to total charge and total current, including the complicated charges and currents in materials at the atomic scale. The "macroscopic" Maxwell equations define two new auxiliary fields that describe the large-scale behaviour of matter without having to consider atomic scale details. However, their use requires experimentally determining parameters for a phenomenological description of the electromagnetic response of materials.

The term "Maxwell's equations" is often used for equivalent alternative formulations. Versions of Maxwell's equations based on the electric and magnetic potentials are preferred for explicitly solving the equations as a boundary value problem, analytical mechanics, or for use in quantum mechanics. The spacetime formulations (i.e., on spacetime rather than space and time separately), are commonly used in high energy and gravitational physics because they make the compatibility of the equations with special and general relativity manifest. In fact, Einstein developed special and general relativity to accommodate the absolute speed of light that drops out of the Maxwell equations with the principle that only relative movement has physical consequences.

Since the mid-20th century, it has been understood that Maxwell's equations are not exact, but a classical field theory approximation of some aspects of the fundamental theory of quantum electrodynamics, although some quantum features, such as quantum entanglement, are completely absent and in no way approximated. (For example, quantum cryptography has no approximate version in Maxwell theory.) In many situations, though, deviations from Maxwell's equations are immeasurably small. Exceptions include nonclassical light, photon–photon scattering, quantum optics, and many other phenomena related to photons or virtual photons.

Formulation in terms of electric and magnetic fields (microscopic or in vacuum version)

In the electric and magnetic field formulation there are four equations. The two inhomogeneous equations describe how the fields vary in space due to sources. Gauss's law describes how electric fields emanate from electric charges. Gauss's law for magnetism describes magnetic fields as closed field lines not due to magnetic monopoles. The two homogeneous equations describe how the fields "circulate" around their respective sources. Ampère's law with Maxwell's addition describes how the magnetic field "circulates" around electric currents and time varying electric fields, while Faraday's law describes how the electric field "circulates" around time varying magnetic fields.

A separate law of nature, the Lorentz force law, describes how the electric and magnetic field act on charged particles and currents. A version of this law was included in the original equations by Maxwell but, by convention, is no longer.

The precise formulation of Maxwell's equations depends on the precise definition of the quantities involved. Conventions differ with the unit systems, because various definitions and dimensions are changed by absorbing dimensionful factors like the speed of light c. This makes constants come out differently. The most common form is based on conventions used when quantities measured using SI units, but other commonly used conventions are used with other units including Gaussian units based on the cgs system, Lorentz–Heaviside units (used mainly in particle physics), and Planck units (used in theoretical physics).

The vector calculus formulation below has become standard. It is mathematically much more convenient than Maxwell's original 20 equations and is due to Oliver Heaviside. The differential and integral equations formulations are mathematically equivalent and are both useful. The integral formulation relates fields within a region of space to fields on the boundary and can often be used to simplify and directly calculate fields from symmetric distributions of charges and currents. On the other hand, the differential equations are purely local and are a more natural starting point for calculating the fields in more complicated (less symmetric) situations, for example using finite element analysis. For formulations using tensor calculus or differential forms, see alternative formulations. For relativistically invariant formulations, see relativistic formulations.

Formulation in SI units convention

Name	Integral equations	Differential equations	Meaning
Gauss's law	${\scriptstyle \partial \Omega }$ $\mathbf {E} \cdot \mathrm {d} \mathbf {S} ={\frac {1}{\varepsilon _{0}}}\iiint _{\Omega }\rho \,\mathrm {d} V$	$\nabla \cdot \mathbf {E} ={\frac {\rho }{\varepsilon _{0}}}$	The electric flux leaving a volume is proportional to the charge inside.
Gauss's law for magnetism	${\scriptstyle \partial \Omega }$ $\mathbf {B} \cdot \mathrm {d} \mathbf {S} =0$	$\nabla \cdot \mathbf {B} =0$	There are no magnetic monopoles; the total magnetic flux through a closed surface is zero.
Maxwell–Faraday equation (Faraday's law of induction)	$\oint _{\partial \Sigma }\mathbf {E} \cdot \mathrm {d} {\boldsymbol {\ell }}=-{\frac {\mathrm {d} }{\mathrm {d} t}}\iint _{\Sigma }\mathbf {B} \cdot \mathrm {d} \mathbf {S}$	$\nabla \times \mathbf {E} =-{\frac {\partial \mathbf {B} }{\partial t}}$	The voltage induced in a closed loop is proportional to the rate of change of the magnetic flux that the loop encloses.
Ampère's circuital law (with Maxwell's addition)	$\oint _{\partial \Sigma }\mathbf {B} \cdot \mathrm {d} {\boldsymbol {\ell }}=\mu _{0}\iint _{\Sigma }\mathbf {J} \cdot \mathrm {d} \mathbf {S} +\mu _{0}\varepsilon _{0}{\frac {\mathrm {d} }{\mathrm {d} t}}\iint _{\Sigma }\mathbf {E} \cdot \mathrm {d} \mathbf {S}$	$\nabla \times \mathbf {B} =\mu _{0}\left(\mathbf {J} +\varepsilon _{0}{\frac {\partial \mathbf {E} }{\partial t}}\right)$	The magnetic field induced around a closed loop is proportional to the electric current plus displacement current (rate of change of electric field) that the loop encloses.

Formulation in Gaussian units convention

Main article: Gaussian units

Gaussian units are a popular system of units, that are part of the centimetre–gram–second system of units (cgs). When using cgs units it is conventional to use a slightly different definition of electric field E_cgs = c⁻¹ E_SI. This implies that the modified electric and magnetic field have the same units (in the SI convention this is not the case making dimensional analysis of the equations different: e.g. for an electromagnetic wave in vacuum $|\mathbf {E} |_{\mathrm {SI} }=|c|_{\mathrm {SI} }|\mathbf {B} |_{\mathrm {SI} }$, ). The CGS system uses a unit of charge defined in such a way that the permittivity of the vacuum ε₀ = ⁠1/4πc⁠, hence μ₀ = ⁠4π/c⁠. These units are sometimes preferred over SI units in the context of special relativity, in which the components of the electromagnetic tensor, the Lorentz covariant object describing the electromagnetic field, have the same unit without constant factors. Using these different conventions, the Maxwell equations become:

Name	Integral equations	Differential equations	Meaning
Gauss's law	${\scriptstyle \partial \Omega }$ $\mathbf {E} \cdot \mathrm {d} \mathbf {S} =4\pi \iiint _{\Omega }\rho \,\mathrm {d} V$	$\nabla \cdot \mathbf {E} =4\pi \rho$	The electric flux leaving a volume is proportional to the charge inside.
Gauss's law for magnetism	${\scriptstyle \partial \Omega }$ $\mathbf {B} \cdot \mathrm {d} \mathbf {S} =0$	$\nabla \cdot \mathbf {B} =0$	There are no magnetic monopoles; the total magnetic flux through a closed surface is zero.
Maxwell–Faraday equation (Faraday's law of induction)	$\oint _{\partial \Sigma }\mathbf {E} \cdot \mathrm {d} {\boldsymbol {\ell }}=-{\frac {1}{c}}{\frac {d}{dt}}\iint _{\Sigma }\mathbf {B} \cdot \mathrm {d} \mathbf {S}$	$\nabla \times \mathbf {E} =-{\frac {1}{c}}{\frac {\partial \mathbf {B} }{\partial t}}$	The voltage induced in a closed loop is proportional to the rate of change of the magnetic flux that the loop encloses.
Ampère's circuital law (with Maxwell's addition)	$\oint _{\partial \Sigma }\mathbf {B} \cdot \mathrm {d} {\boldsymbol {\ell }}={\frac {1}{c}}\left(4\pi \iint _{\Sigma }\mathbf {J} \cdot \mathrm {d} \mathbf {S} +{\frac {d}{dt}}\iint _{\Sigma }\mathbf {E} \cdot \mathrm {d} \mathbf {S} \right)$	$\nabla \times \mathbf {B} ={\frac {1}{c}}\left(4\pi \mathbf {J} +{\frac {\partial \mathbf {E} }{\partial t}}\right)$	The magnetic field integrated around a closed loop is proportional to the electric current plus displacement current (rate of change of electric field) that the loop encloses.

Key to the notation

Symbols in bold represent vector quantities, and symbols in italics represent scalar quantities, unless otherwise indicated.

The equations introduce the electric field, E, a vector field, and the magnetic field, B, a pseudovector field, each generally having a time and location dependence. The sources are

• the total electric charge density (total charge per unit volume), ρ, and
• the total electric current density (total current per unit area), J.

The universal constants appearing in the equations are

• the permittivity of free space, ε₀, and
• the permeability of free space, μ₀.

Differential equations

In the differential equations,

• the nabla symbol, ∇, denotes the three-dimensional gradient operator,
• the ∇⋅ symbol denotes the divergence operator,
• the ∇× symbol denotes the curl operator.

Integral equations

In the integral equations,

• Ω is any fixed volume with closed boundary surface ∂Ω, and
• Σ is any fixed surface with closed boundary curve ∂Σ,

Here a fixed volume or surface means that it does not change over time. The equations are correct, complete and a little easier to interpret with time-independent surfaces. For example, since the surface is time-independent, we can bring the differentiation under the integral sign in Faraday's law:

${\frac {d}{dt}}\iint _{\Sigma }\mathbf {B} \cdot \mathrm {d} \mathbf {S} =\iint _{\Sigma }{\frac {\partial \mathbf {B} }{\partial t}}\cdot \mathrm {d} \mathbf {S} \,,$

Maxwell's equations can be formulated with possibly time-dependent surfaces and volumes by substituting the left hand side with the right hand side in the integral equation version of the Maxwell equations.

• ${\scriptstyle \partial \Omega }$ is a surface integral over the boundary surface ∂Ω, with the loop indicating the surface is closed
• $\iiint _{\Omega }$ is a volume integral over the volume Ω,
• $\oint _{\partial \Sigma }$ is a line integral around the boundary curve ∂Σ, with the loop indicating the curve is closed.
• $\iint _{\Sigma }$ is a surface integral over the surface Σ,
• The total electric charge Q enclosed in Ω is the volume integral over Ω of the charge density ρ (see the "macroscopic formulation" section below):

$Q=\iiint _{\Omega }\rho \ \mathrm {d} V,$

where dV is the volume element.

• The net electric current I is the surface integral of the electric current density J passing through a fixed surface, Σ:

$I=\iint _{\Sigma }\mathbf {J} \cdot \mathrm {d} \mathbf {S} ,$

where dS denotes the vector element of surface area S, normal to surface Σ. (Vector area is sometimes denoted by A rather than S, but this conflicts with the notation for magnetic potential).

Relationship between differential and integral formulations

The equivalence of the differential and integral formulations are a consequence of the Gauss divergence theorem and the Kelvin–Stokes theorem.

Flux and divergence

Volume Ω and its closed boundary ∂Ω, containing (respectively enclosing) a source (+) and sink (−) of a vector field F. Here, F could be the E field with source electric charges, but not the B field which has no magnetic charges as shown. The outward unit normal is n.

The "sources of the fields" (i.e. their divergence) can be determined from the surface integrals of the fields through the closed surface ∂Ω. I.e. the electric flux is

${\scriptstyle \partial \Omega }$ $\mathbf {E} \cdot \mathrm {d} \mathbf {S} =\iiint _{\Omega }\nabla \cdot \mathbf {E} \,\mathrm {d} V$

where the last equality uses the Gauss divergence theorem. Using the integral version of Gauss's equation we can rewrite this to

$\iiint _{\Omega }\left(\nabla \cdot \mathbf {E} -{\frac {\rho }{\epsilon _{0}}}\right)\,\mathrm {d} V=0$

Since Ω can be chosen arbitrarily, e.g. as an arbitrary small ball with arbitrary center, this implies that the integrand must be zero, which is the differential equations formulation of Gauss equation up to a trivial rearrangement. Gauss's law for magnetism in differential equations form follows likewise from the integral form by rewriting the magnetic flux

${\scriptstyle \partial \Omega }$ $\mathbf {B} \cdot \mathrm {d} \mathbf {S} =\iiint _{\Omega }\nabla \cdot \mathbf {B} \,\mathrm {d} V=0$.

Circulation and curl

Surface Σ with closed boundary ∂Σ. F could be the E or B fields. Again, n is the unit normal. (The curl of a vector field doesn't literally look like the "circulations", this is a heuristic depiction.)

The "circulation of the fields" (i.e. their curls) can be determined from the line integrals of the fields around the closed curve ∂Σ. E.g. for the magnetic field

$\oint _{\partial \Sigma }\mathbf {B} \cdot \mathrm {d} {\boldsymbol {\ell }}=\iint _{\Sigma }(\nabla \times \mathbf {B} )\cdot \mathrm {d} \mathbf {S}$,

where we used the Kelvin–Stokes theorem. Using the modified Ampere law in integral form and the writing the time derivative of the flux as the surface integral of the partial time derivative of E we conclude that

$\iint _{\Sigma }\left(\nabla \times \mathbf {B} -\mu _{0}\left(\mathbf {J} +\epsilon _{0}{\frac {\partial \mathbf {E} }{\partial t}}\right)\right)\cdot \mathrm {d} \mathbf {S} =0$.

Since Σ can be chosen arbitrarily, e.g. as an arbitrary small, arbitrary oriented, and arbitrary centered disk, we conclude that the integrand must be zero. This is Ampere's modified law in differential equations form up to a trivial rearrangement. Likewise, the Faraday law in differential equations form follows from rewriting the integral form using the Kelvin–Stokes theorem.

The line integrals and curls are analogous to quantities in classical fluid dynamics: the circulation of a fluid is the line integral of the fluid's flow velocity field around a closed loop, and the vorticity of the fluid is the curl of the velocity field.

Conceptual descriptions

Gauss's law

Gauss's law describes the relationship between a static electric field and the electric charges that cause it: The static electric field points away from positive charges and towards negative charges. In the field line description, electric field lines begin only at positive electric charges and end only at negative electric charges. 'Counting' the number of field lines passing through a closed surface, therefore, yields the total charge (including bound charge due to polarization of material) enclosed by that surface divided by dielectricity of free space (the vacuum permittivity). More technically, it relates the electric flux through any hypothetical closed "Gaussian surface" to the enclosed electric charge.

Gauss's law for magnetism: magnetic field lines never begin nor end but form loops or extend to infinity as shown here with the magnetic field due to a ring of current.

Gauss's law for magnetism

Gauss's law for magnetism states that there are no "magnetic charges" (also called magnetic monopoles), analogous to electric charges. Instead, the magnetic field due to materials is generated by a configuration called a dipole. Magnetic dipoles are best represented as loops of current but resemble positive and negative 'magnetic charges', inseparably bound together, having no net 'magnetic charge'. In terms of field lines, this equation states that magnetic field lines neither begin nor end but make loops or extend to infinity and back. In other words, any magnetic field line that enters a given volume must somewhere exit that volume. Equivalent technical statements are that the sum total magnetic flux through any Gaussian surface is zero, or that the magnetic field is a solenoidal vector field.

Faraday's law

In a geomagnetic storm, a surge in the flux of charged particles temporarily alters Earth's magnetic field, which induces electric fields in Earth's atmosphere, thus causing surges in electrical power grids. (Not to scale.)

The Maxwell–Faraday version of Faraday's law of induction describes how a time varying magnetic field creates ("induces") an electric field. This dynamically induced electric field has closed field lines similar to a magnetic field, unless superposed by a static (charge induced) electric field. This aspect of electromagnetic induction is the operating principle behind many electric generators: for example, a rotating bar magnet creates a changing magnetic field, which in turn generates an electric field in a nearby wire.

Ampère's law with Maxwell's addition

Magnetic core memory (1954) is an application of Ampère's law. Each core stores one bit of data.

Ampère's law with Maxwell's addition states that magnetic fields can be generated in two ways: by electric current (this was the original "Ampère's law") and by changing electric fields (this was "Maxwell's addition").

Maxwell's addition to Ampère's law is particularly important: it makes the set of equations mathematically consistent for non static fields, without changing the laws of Ampere and Gauss for static fields. However, as a consequence, it predicts that a changing magnetic field induces an electric field and vice versa. Therefore, these equations allow self-sustaining "electromagnetic waves" to travel through empty space (see electromagnetic wave equation).

The speed calculated for electromagnetic waves, which could be predicted from experiments on charges and currents, exactly matches the speed of light; indeed, light is one form of electromagnetic radiation (as are X-rays, radio waves, and others). Maxwell understood the connection between electromagnetic waves and light in 1861, thereby unifying the theories of electromagnetism and optics.

Charge conservation

The invariance of charge can be derived as a corollary of Maxwell's equations. The left hand side of the modified Ampere's Law has zero divergence by the div–curl identity. Combining the right hand side, Gauss's law, and interchange of derivatives gives:

$0=\nabla \cdot \nabla \times \mathbf {B} =\mu _{0}\left(\nabla \cdot \mathbf {J} +\varepsilon _{0}{\frac {\partial }{\partial t}}\nabla \cdot \mathbf {E} \right)=\mu _{0}\left(\nabla \cdot \mathbf {J} +{\frac {\partial \rho }{\partial t}}\right)$

i.e.

${\frac {\partial \rho }{\partial t}}+\nabla \cdot \mathbf {J} =0$.

By the Gauss Divergence Theorem, this means the rate of change of charge in a fixed volume equals the net current flowing through the boundary:

${\frac {d}{dt}}Q_{\Omega }={\frac {d}{dt}}\iiint _{\Omega }\rho \mathrm {d} V=-$ ${\scriptstyle \partial \Omega }$ ${\mathbf {J}}\cdot {\rm {d}}{\mathbf {S}}=-I_{\partial \Omega }.$

In particular, in an isolated system the total charge is conserved.

Vacuum equations, electromagnetic waves and speed of light

Further information: Electromagnetic wave equation, Inhomogeneous electromagnetic wave equation, and Sinusoidal plane-wave solutions of the electromagnetic wave equation

This 3D diagram shows a plane linearly polarized wave propagating from left to right with the same wave equations where E = E₀ sin(−ωt + k ⋅ r) and B = B₀ sin(−ωt + k ⋅ r)

In a region with no charges (ρ = 0) and no currents (J = 0), such as in a vacuum, Maxwell's equations reduce to:

${\begin{aligned}\nabla \cdot \mathbf {E} &=0\quad &\nabla \times \mathbf {E} &=-{\frac {\partial \mathbf {B} }{\partial t}},\\\nabla \cdot \mathbf {B} &=0\quad &\nabla \times \mathbf {B} &={\frac {1}{c^{2}}}{\frac {\partial \mathbf {E} }{\partial t}}.\end{aligned}}$

Taking the curl (∇×) of the curl equations, and using the curl of the curl identity ∇ × (∇ × X) = ∇(∇·X) − ∇²X we obtain the wave equations

${\begin{aligned}{\frac {1}{c^{2}}}{\frac {\partial ^{2}\mathbf {E} }{\partial t^{2}}}-\nabla ^{2}\mathbf {E} =0\\{\frac {1}{c^{2}}}{\frac {\partial ^{2}\mathbf {B} }{\partial t^{2}}}-\nabla ^{2}\mathbf {B} =0\end{aligned}}$

which identify

$c={\frac {1}{\sqrt {\mu _{0}\varepsilon _{0}}}}=2.99792458\times 10^{8}\,{\text{m}}\,{\text{s}}^{-1}$

with the speed of light in free space. In materials with relative permittivity, ε_r, and relative permeability, μ_r, the phase velocity of light becomes

$v_{\text{p}}={\frac {1}{\sqrt {\mu _{0}\mu _{\text{r}}\varepsilon _{0}\varepsilon _{\text{r}}}}}$

which is usually less than c.

In addition, E and B are mutually perpendicular to each other and the direction of wave propagation, and are in phase with each other. A sinusoidal plane wave is one special solution of these equations. Maxwell's equations explain how these waves can physically propagate through space. The changing magnetic field creates a changing electric field through Faraday's law. In turn, that electric field creates a changing magnetic field through Maxwell's addition to Ampère's law. This perpetual cycle allows these waves, now known as electromagnetic radiation, to move through space at velocity c.

Macroscopic formulation

The above equations are the "microscopic" version of Maxwell's equations, expressing the electric and the magnetic fields in terms of the (possibly atomic-level) charges and currents present. This is sometimes called the "general" form, but the macroscopic version below is equally general, the difference being one of bookkeeping.

The microscopic version is sometimes called "Maxwell's equations in a vacuum": this refers to the fact that the material medium is not built into the structure of the equations, but appears only in the charge and current terms. The microscopic version was introduced by Lorentz, who tried to use it to derive the macroscopic properties of bulk matter from its microscopic constituents.

"Maxwell's macroscopic equations", also known as Maxwell's equations in matter, are more similar to those that Maxwell introduced himself.

Name	Integral equations (SI convention)	Differential equations (SI convention)	Differential equations (Gaussian convention)
Gauss's law	${\scriptstyle \partial \Omega }$ $\mathbf {D} \cdot \mathrm {d} \mathbf {S} =\iiint _{\Omega }\rho _{\text{f}}\,\mathrm {d} V$	$\nabla \cdot \mathbf {D} =\rho _{\text{f}}$	$\nabla \cdot \mathbf {D} =4\pi \rho _{\text{f}}$
Gauss's law for magnetism	${\scriptstyle \partial \Omega }$ $\mathbf {B} \cdot \mathrm {d} \mathbf {S} =0$	$\nabla \cdot \mathbf {B} =0$	$\nabla \cdot \mathbf {B} =0$
Maxwell–Faraday equation (Faraday's law of induction)	$\oint _{\partial \Sigma }\mathbf {E} \cdot \mathrm {d} {\boldsymbol {\ell }}=-{\frac {d}{dt}}\iint _{\Sigma }\mathbf {B} \cdot \mathrm {d} \mathbf {S}$	$\nabla \times \mathbf {E} =-{\frac {\partial \mathbf {B} }{\partial t}}$	$\nabla \times \mathbf {E} =-{\frac {1}{c}}{\frac {\partial \mathbf {B} }{\partial t}}$
Ampère's circuital law (with Maxwell's addition)	$\oint _{\partial \Sigma }\mathbf {H} \cdot \mathrm {d} {\boldsymbol {\ell }}=\iint _{\Sigma }\mathbf {J} _{\text{f}}\cdot \mathrm {d} \mathbf {S} +{\frac {d}{dt}}\iint _{\Sigma }\mathbf {D} \cdot \mathrm {d} \mathbf {S}$	$\nabla \times \mathbf {H} =\mathbf {J} _{\text{f}}+{\frac {\partial \mathbf {D} }{\partial t}}$	$\nabla \times \mathbf {H} ={\frac {1}{c}}\left(4\pi \mathbf {J} _{\text{f}}+{\frac {\partial \mathbf {D} }{\partial t}}\right)$

Unlike the "microscopic" equations, the "macroscopic" equations separate out the bound charge Q_b and bound current I_b to obtain equations that depend only on the free charges Q_f and currents I_f. This factorization can be made by splitting the total electric charge and current as follows:

${\begin{aligned}Q&=Q_{\text{f}}+Q_{\text{b}}=\iiint _{\Omega }\left(\rho _{\text{f}}+\rho _{\text{b}}\right)\,\mathrm {d} V=\iiint _{\Omega }\rho \,\mathrm {d} V\\I&=I_{\text{f}}+I_{\text{b}}=\iint _{\Sigma }\left(\mathbf {J} _{\text{f}}+\mathbf {J} _{\text{b}}\right)\cdot \mathrm {d} \mathbf {S} =\iint _{\Sigma }\mathbf {J} \cdot \mathrm {d} \mathbf {S} \end{aligned}}$

Correspondingly, the total current density J splits into free J_f and bound J_b components, and similarly the total charge density ρ splits into free ρ_f and bound ρ_b parts.

The cost of this factorization is that additional fields, the displacement field D and the magnetizing field H, are defined and need to be determined. Phenomenological constituent equations relate the additional fields to the electric field E and the magnetic B-field, often through a simple linear relation.

For a detailed description of the differences between the microscopic (total charge and current including material contributes or in air/vacuum) and macroscopic (free charge and current; practical to use on materials) variants of Maxwell's equations, see below.

Bound charge and current

Main articles: Current density, Bound charge, and Bound current

Left: A schematic view of how an assembly of microscopic dipoles produces opposite surface charges as shown at top and bottom. Right: How an assembly of microscopic current loops add together to produce a macroscopically circulating current loop. Inside the boundaries, the individual contributions tend to cancel, but at the boundaries no cancelation occurs.

When an electric field is applied to a dielectric material its molecules respond by forming microscopic electric dipoles – their atomic nuclei move a tiny distance in the direction of the field, while their electrons move a tiny distance in the opposite direction. This produces a macroscopic bound charge in the material even though all of the charges involved are bound to individual molecules. For example, if every molecule responds the same, similar to that shown in the figure, these tiny movements of charge combine to produce a layer of positive bound charge on one side of the material and a layer of negative charge on the other side. The bound charge is most conveniently described in terms of the polarization P of the material, its dipole moment per unit volume. If P is uniform, a macroscopic separation of charge is produced only at the surfaces where P enters and leaves the material. For non-uniform P, a charge is also produced in the bulk.

Somewhat similarly, in all materials the constituent atoms exhibit magnetic moments that are intrinsically linked to the angular momentum of the components of the atoms, most notably their electrons. The connection to angular momentum suggests the picture of an assembly of microscopic current loops. Outside the material, an assembly of such microscopic current loops is not different from a macroscopic current circulating around the material's surface, despite the fact that no individual charge is traveling a large distance. These bound currents can be described using the magnetization M.

The very complicated and granular bound charges and bound currents, therefore, can be represented on the macroscopic scale in terms of P and M which average these charges and currents on a sufficiently large scale so as not to see the granularity of individual atoms, but also sufficiently small that they vary with location in the material. As such, Maxwell's macroscopic equations ignore many details on a fine scale that can be unimportant to understanding matters on a gross scale by calculating fields that are averaged over some suitable volume.

Auxiliary fields, polarization and magnetization

The definitions (not constitutive relations) of the auxiliary fields are:

${\begin{aligned}\mathbf {D} (\mathbf {r} ,t)&=\varepsilon _{0}\mathbf {E} (\mathbf {r} ,t)+\mathbf {P} (\mathbf {r} ,t)\\\mathbf {H} (\mathbf {r} ,t)&={\frac {1}{\mu _{0}}}\mathbf {B} (\mathbf {r} ,t)-\mathbf {M} (\mathbf {r} ,t)\end{aligned}}$

where P is the polarization field and M is the magnetization field which are defined in terms of microscopic bound charges and bound currents respectively. The macroscopic bound charge density ρ_b and bound current density J_b in terms of polarization P and magnetization M are then defined as

${\begin{aligned}\rho _{\text{b}}&=-\nabla \cdot \mathbf {P} \\\mathbf {J} _{\text{b}}&=\nabla \times \mathbf {M} +{\frac {\partial \mathbf {P} }{\partial t}}\end{aligned}}$

If we define the total, bound, and free charge and current density by

${\begin{aligned}\rho &=\rho _{\text{b}}+\rho _{\text{f}},\\\mathbf {J} &=\mathbf {J} _{\text{b}}+\mathbf {J} _{\text{f}},\end{aligned}}$

and use the defining relations above to eliminate D, and H, the "macroscopic" Maxwell's equations reproduce the "microscopic" equations.

Constitutive relations

Main article: Constitutive equation § Electromagnetism

In order to apply 'Maxwell's macroscopic equations', it is necessary to specify the relations between displacement field D and the electric field E, as well as the magnetizing field H and the magnetic field B. Equivalently, we have to specify the dependence of the polarisation P (hence the bound charge) and the magnetisation M (hence the bound current) on the applied electric and magnetic field. The equations specifying this response are called constitutive relations. For real-world materials, the constitutive relations are rarely simple, except approximately, and usually determined by experiment. See the main article on constitutive relations for a fuller description.

For materials without polarisation and magnetisation, the constitutive relations are (by definition)

$\mathbf {D} =\varepsilon _{0}\mathbf {E} ,\quad \mathbf {H} ={\frac {1}{\mu _{0}}}\mathbf {B}$

where ε₀ is the permittivity of free space and μ₀ the permeability of free space. Since there is no bound charge, the total and the free charge and current are equal.

An alternative viewpoint on the microscopic equations is that they are the macroscopic equations together with the statement that vacuum behaves like a perfect linear "material" without additional polarisation and magnetisation. More generally, for linear materials the constitutive relations are

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

where ε is the permittivity and μ the permeability of the material. For the displacement field D the linear approximation is usually excellent because for all but the most extreme electric fields or temperatures obtainable in the laboratory (high power pulsed lasers) the interatomic electric fields of materials of the order of 10¹¹ V/m are much higher than the external field. For the magnetizing field $\mathbf {H}$, however, the linear approximation can break down in common materials like iron leading to phenomena like hysteresis. Even the linear case can have various complications, however.

• For homogeneous materials, ε and μ are constant throughout the material, while for inhomogeneous materials they depend on location within the material (and perhaps time).
• For isotropic materials, ε and μ are scalars, while for anisotropic materials (e.g. due to crystal structure) they are tensors.
• Materials are generally dispersive, so ε and μ depend on the frequency of any incident EM waves.

Even more generally, in the case of non-linear materials (see for example nonlinear optics), D and P are not necessarily proportional to E, similarly H or M is not necessarily proportional to B. In general D and H depend on both E and B, on location and time, and possibly other physical quantities.

In applications one also has to describe how the free currents and charge density behave in terms of E and B possibly coupled to other physical quantities like pressure, and the mass, number density, and velocity of charge-carrying particles. E.g., the original equations given by Maxwell (see History of Maxwell's equations) included Ohms law in the form

$\mathbf {J} _{\text{f}}=\sigma \mathbf {E} \,.$

Alternative formulations

For an overview, see Mathematical descriptions of the electromagnetic field.

For the equations in quantum field theory, see quantum electrodynamics.

Following is a summary of some of the numerous other ways to write the microscopic Maxwell's equations, showing they can be formulated using different mathematical formalisms. In addition, we formulate the equations using "potentials". Originally they were introduced as a convenient way to solve the homogeneous equations, but it was originally thought that all the observable physics was contained in the electric and magnetic fields (or relativistically, the Faraday tensor). The potentials play a central role in quantum mechanics, however, and act quantum mechanically with observable consequences even when the electric and magnetic fields vanish (Aharonov–Bohm effect). See the main articles for the details of each formulation. SI units are used throughout.

Formalism	Formulation	Homogeneous equations	Inhomogeneous equations
Vector calculus	Fields 3D Euclidean space + time	${\begin{aligned}\nabla \cdot \mathbf {B} &=0\\\nabla \times \mathbf {E} +{\frac {\partial \mathbf {B} }{\partial t}}&=0\end{aligned}}$	${\begin{aligned}\nabla \cdot \mathbf {E} &={\frac {\rho }{\varepsilon _{0}}}\\\nabla \times \mathbf {B} -{\frac {1}{c^{2}}}{\frac {\partial \mathbf {E} }{\partial t}}&=\mu _{0}\mathbf {J} \end{aligned}}$
	Potentials (any gauge) 3D Euclidean space + time	${\begin{aligned}\mathbf {B} &=\mathbf {\nabla } \times \mathbf {A} \\\mathbf {E} &=-\mathbf {\nabla } \varphi -{\frac {\partial \mathbf {A} }{\partial t}}\end{aligned}}$	${\begin{aligned}-\nabla ^{2}\varphi -{\frac {\partial }{\partial t}}\left(\mathbf {\nabla } \cdot \mathbf {A} \right)&={\frac {\rho }{\varepsilon _{0}}}\\(-\nabla ^{2}+{\frac {1}{c^{2}}}{\frac {\partial ^{2}}{\partial t^{2}}})\mathbf {A} +\mathbf {\nabla } \left(\mathbf {\nabla } \cdot \mathbf {A} +{\frac {1}{c^{2}}}{\frac {\partial \varphi }{\partial t}}\right)&=\mu _{0}\mathbf {J} \end{aligned}}$
	Potentials (Lorenz gauge) 3D Euclidean space + time	${\begin{aligned}\mathbf {B} &=\mathbf {\nabla } \times \mathbf {A} \\\mathbf {E} &=-\mathbf {\nabla } \varphi -{\frac {\partial \mathbf {A} }{\partial t}}\\\mathbf {\nabla } \cdot \mathbf {A} &=-{\frac {1}{c^{2}}}{\frac {\partial \varphi }{\partial t}}\\\end{aligned}}$	${\begin{aligned}(-\nabla ^{2}+{\frac {1}{c^{2}}}{\frac {\partial ^{2}}{\partial t^{2}}})\varphi &={\frac {\rho }{\varepsilon _{0}}}\\(-\nabla ^{2}+{\frac {1}{c^{2}}}{\frac {\partial ^{2}}{\partial t^{2}}})\mathbf {A} &=\mu _{0}\mathbf {J} \end{aligned}}$
Tensor calculus	Fields space + time spatial metric independent of time	${\begin{aligned}\partial _{[i}B_{jk]}&=\\\nabla _{[i}B_{jk]}&=0\\\partial _{[i}E_{j]}+{\frac {\partial B_{ij}}{\partial t}}&=\\\nabla _{[i}E_{j]}+{\frac {\partial B_{ij}}{\partial t}}&=0\end{aligned}}$	${\begin{aligned}{\frac {1}{\sqrt {h}}}\partial _{i}{\sqrt {h}}E^{i}&=\\\nabla _{i}E^{i}&={\frac {\rho }{\epsilon _{0}}}\\-{\frac {1}{\sqrt {h}}}\partial _{i}{\sqrt {h}}B^{ij}-{\frac {1}{c^{2}}}{\frac {\partial }{\partial t}}E^{j}&=&\\-\nabla _{i}B^{ij}-{\frac {1}{c^{2}}}{\frac {\partial E^{j}}{\partial t}}&=\mu _{0}J^{j}\\\end{aligned}}$
	Potentials space (with topological restrictions) + time spatial metric independent of time	${\begin{aligned}B_{ij}&=\partial _{[i}A_{j]}\\&=\nabla _{[i}A_{j]}\\E_{i}&=-{\frac {\partial A_{i}}{\partial t}}-\partial _{i}\phi \\&=-{\frac {\partial A_{i}}{\partial t}}-\nabla _{i}\phi \\\end{aligned}}$	${\begin{aligned}-{\frac {1}{\sqrt {h}}}\partial _{i}{\sqrt {h}}(\partial ^{i}\phi +{\frac {\partial A^{i}}{\partial t}})&=\\-\nabla _{i}\nabla ^{i}\phi -{\frac {\partial }{\partial t}}\nabla _{i}A^{i}&={\frac {\rho }{\epsilon _{0}}}\\-{\frac {1}{\sqrt {h}}}\partial _{i}({\sqrt {h}}h^{im}h^{jn}\partial _{[m}A_{n]})+{\frac {1}{c^{2}}}{\frac {\partial }{\partial t}}({\frac {\partial A^{j}}{\partial t}}+\partial ^{j}\phi )&=\\-\nabla _{i}\nabla ^{i}A^{j}+{\frac {1}{c^{2}}}{\frac {\partial ^{2}A^{j}}{\partial t^{2}}}+R_{i}^{j}A^{i}+\nabla ^{j}(\nabla _{i}A^{i}+{\frac {1}{c^{2}}}{\frac {\partial \phi }{\partial t}})&=\mu _{0}J^{j}\\\end{aligned}}$
	Potentials (Lorenz gauge) space (with topological restrictions) + time spatial metric independent of time	${\begin{aligned}B_{ij}&=\partial _{[i}A_{j]}\\&=\nabla _{[i}A_{j]}\\E_{i}&=-{\frac {\partial A_{i}}{\partial t}}-\partial _{i}\phi \\&=-{\frac {\partial A_{i}}{\partial t}}-\nabla _{i}\phi \\\nabla _{i}A^{i}&=-{\frac {1}{c^{2}}}{\frac {\partial \phi }{\partial t}}\\\end{aligned}}$	${\begin{aligned}-\nabla _{i}\nabla ^{i}\phi +{\frac {1}{c^{2}}}{\frac {\partial ^{2}\phi }{\partial t^{2}}}&={\frac {\rho }{\epsilon _{0}}}\\-\nabla _{i}\nabla ^{i}A^{j}+{\frac {1}{c^{2}}}{\frac {\partial ^{2}A^{j}}{\partial t^{2}}}+R_{i}^{j}A^{i}&=\mu _{0}J^{j}\\\end{aligned}}$
Differential forms	Fields Any space + time	${\begin{aligned}dB&=0\\dE+{\frac {\partial B}{\partial t}}&=0\\\end{aligned}}$	${\begin{aligned}d*E={\frac {\rho }{\epsilon _{0}}}\\d*B-{\frac {1}{c^{2}}}{\frac {\partial *E}{\partial t}}={\mu _{0}}J\\\end{aligned}}$
	Potentials (any gauge) Any space (with topological restrictions) + time	${\begin{aligned}B&=dA\\E&=-d\phi -{\frac {\partial A}{\partial t}}\\\end{aligned}}$	${\begin{aligned}-d*(d\phi +{\frac {\partial A}{\partial t}})&={\frac {\rho }{\epsilon _{0}}}\\d*dA+{\frac {1}{c^{2}}}{\frac {\partial }{\partial t}}*(d\phi -{\frac {\partial A}{\partial t}})&=\mu _{0}J\\\end{aligned}}$
	Potential (Lorenz Gauge) Any space (with topological restrictions) + time spatial metric independent of time	${\begin{aligned}B&=dA\\E&=-d\phi -{\frac {\partial A}{\partial t}}\\d*A&=-*{\frac {1}{c^{2}}}{\frac {\partial \phi }{\partial t}}\\\end{aligned}}$	${\begin{aligned}*(-\Delta \phi +{\frac {1}{c^{2}}}{\frac {\partial ^{2}}{\partial t^{2}}}\phi )&={\frac {\rho }{\epsilon _{0}}}\\*(-\Delta A+{\frac {1}{c^{2}}}{\frac {\partial ^{2}A}{\partial ^{2}t}})&=\mu _{0}J\\\end{aligned}}$

where

• In the vector formulation on Euclidean space + time, φ is the electrical potential, and A is the vector potential.

Relativistic formulations

For the equations in special relativity, see classical electromagnetism and special relativity and covariant formulation of classical electromagnetism.

For the equations in general relativity, see Maxwell's equations in curved spacetime.

The Maxwell equations can also be formulated on a spacetime-like Minkowski space where space and time are treated on equal footing. The direct spacetime formulations make manifest that the Maxwell equations are relativistically invariant. Because of this symmetry electric and magnetic field are treated on equal footing and are recognised as components of the Faraday tensor. This reduces the four Maxwell equations to two, which simplifies the equations, although we can no longer use the familiar vector formulation. In fact the Maxwell equations in the space + time formulation are not Galileo invariant and have Lorenz invariance a hidden symmetry. This was a major source of inspiration for the development of relativity theory. The space + time formulation is not a non-relativistic approximation, however, they describe the same physics by simply renaming variables. For this reason the relativistic invariant equations are usually simply called the Maxwell equations as well.

Formalism	Formulation	Homogeneous equations	Inhomogeneous equations
Tensor calculus	Fields Minkowski space	$\partial _{[\alpha }F_{\beta \gamma ]}=0$	$\partial _{\alpha }F^{\alpha \beta }=\mu _{0}J^{\beta }$
	Potentials (any gauge) Minkowski space	$F_{\alpha \beta }=2\partial _{[\alpha }A_{\beta ]}$	$2\partial _{\alpha }\partial ^{[\alpha }A^{\beta ]}=\mu _{0}J^{\beta }$
	Potentials (Lorenz gauge) Minkowski space	${\begin{aligned}F_{\alpha \beta }&=2\partial _{[\alpha }A_{\beta ]}\\\partial _{\alpha }A^{\alpha }&=0\end{aligned}}$	$\partial _{\alpha }\partial ^{\alpha }A^{\beta }=\mu _{0}J^{\beta }$
	Fields Any spacetime	${\begin{aligned}\partial _{[\alpha }F_{\beta \gamma ]}&=\\\nabla _{[\alpha }F_{\beta \gamma ]}&=0\end{aligned}}$	${\begin{aligned}{\frac {1}{\sqrt {-g}}}\partial _{\alpha }({\sqrt {-g}}F^{\alpha \beta })&=\\\nabla _{\alpha }F^{\alpha \beta }&=\mu _{0}J^{\beta }\end{aligned}}$
	Potentials (any gauge) Any spacetime (with topological restrictions)	${\begin{aligned}F_{\alpha \beta }&=2\partial _{[\alpha }A_{\beta ]}\\&=2\nabla _{[\alpha }A_{\beta ]}\end{aligned}}$	${\begin{aligned}{\frac {2}{\sqrt {-g}}}\partial _{\alpha }({\sqrt {-g}}g^{\alpha \mu }g^{\beta \nu }\partial _{[\mu }A_{\nu ]})&=\\2\nabla _{\alpha }(\nabla ^{[\alpha }A^{\beta ]})&=\mu _{0}J^{\beta }\end{aligned}}$
	Potentials (Lorenz gauge) Any spacetime (with topological restrictions)	${\begin{aligned}F_{\alpha \beta }&=2\partial _{[\alpha }A_{\beta ]}\\&=2\nabla _{[\alpha }A_{\beta ]}\\\nabla _{\alpha }A^{\alpha }&=0\end{aligned}}$	$\nabla _{\alpha }\nabla ^{\alpha }A^{\beta }-R^{\beta }{}_{\alpha }A^{\alpha }=\mu _{0}J^{\beta }$
Differential forms	Fields Any spacetime	$\mathrm {d} F=0$	$\mathrm {d} {\star }F=\mu _{0}J$
	Potentials (any gauge) Any spacetime (with topological restrictions)	$F=\mathrm {d} A$	$\mathrm {d} {\star }\mathrm {d} A=\mu _{0}J$
	Potentials (Lorenz gauge) Any spacetime (with topological restrictions)	${\begin{aligned}F&=\mathrm {d} A\\\mathrm {d} {\star }A&=0\end{aligned}}$	${\star }\Box A=\mu _{0}J$

• In the tensor calculus formulation, the electromagnetic tensor F_αβ is an antisymmetric covariant rank 2 tensor; the four-potential, A_α, is a covariant vector; the current, J^α, is a vector; the square brackets, [ ], denote antisymmetrization of indices; ∂_α is the derivative with respect to the coordinate, x^α. In Minkowski space coordinates are chosen with respect to an inertial frame; (x^α) = (ct,x,y,z), so that the metric tensor used to raise and lower indices is η_αβ = diag(1,−1,−1,−1). The d'Alembert operator on Minkowski space is ◻ = ∂_α∂^α as in the vector formulation. In general spacetimes, the coordinate system x^α is arbitrary, the covariant derivative ∇_α, the Ricci tensor, R_αβ and raising and lowering of indices are defined by the Lorentzian metric, g_αβ and the d'Alembert operator is defined as ◻ = ∇_α∇^α. The topological restriction is that the second real cohomology group of the space vanishes (see the differential form formulation for an explanation). Note that this is violated for Minkowski space with a line removed, which can model a (flat) spacetime with a point-like monopole on the complement of the line.
• In the differential form formulation on arbitrary space times, F = F_αβdx^α ∧ dx^β is the electromagnetic tensor considered as a 2-form, A = A_αdx^α is the potential 1-form, J is the current 3-form, d is the exterior derivative, and ${\star }$ is the Hodge star on forms defined (up to its an orientation, i.e. its sign) by the Lorentzian metric of spacetime. Note that in the special case of 2-forms such as F, the Hodge star ${\star }$ depends on the metric tensor only for its local scale. This means that, as formulated, the differential form field equations are conformally invariant, but the Lorenz gauge condition breaks conformal invariance. The operator $\Box =(-{\star }\mathrm {d} {\star }\mathrm {d} -\mathrm {d} {\star }\mathrm {d} {\star })$ is the d'Alembert–Laplace–Beltrami operator on 1-forms on an arbitrary Lorentzian spacetime. The topological condition is again that the second real cohomology group is trivial. By the isomorphism with the second de Rham cohomology this condition means that every closed 2-form is exact.

Other formalisms include the geometric algebra formulation and a matrix representation of Maxwell's equations. Historically, a quaternionic formulation was used.

Solutions

Maxwell's equations are partial differential equations that relate the electric and magnetic fields to each other and to the electric charges and currents. Often, the charges and currents are themselves dependent on the electric and magnetic fields via the Lorentz force equation and the constitutive relations. These all form a set of coupled partial differential equations, which are often very difficult to solve. In fact, the solutions of these equations encompass all the diverse phenomena in the entire field of classical electromagnetism. A thorough discussion is far beyond the scope of the article, but some general notes follow.

Like any differential equation, boundary conditions and initial conditions are necessary for a unique solution. For example, even with no charges and no currents anywhere in spacetime, many solutions to Maxwell's equations are possible, not just the obvious solution E = B = 0. Another solution is E = constant, B = constant, while yet other solutions have electromagnetic waves filling spacetime. In some cases, Maxwell's equations are solved through infinite space, and boundary conditions are given as asymptotic limits at infinity. In other cases, Maxwell's equations are solved in just a finite region of space, with appropriate boundary conditions on that region: For example, the boundary could be an artificial absorbing boundary representing the rest of the universe, or periodic boundary conditions, or (as with a waveguide or cavity resonator) the boundary conditions may describe the walls that isolate a small region from the outside world.

Jefimenko's equations (or the closely related Liénard–Wiechert potentials) are the explicit solution to Maxwell's equations for the electric and magnetic fields created by any given distribution of charges and currents. It assumes specific initial conditions to obtain the so-called "retarded solution", where the only fields present are the ones created by the charges. Jefimenko's equations are not so helpful in situations when the charges and currents are themselves affected by the fields they create.

Numerical methods for differential equations can be used to approximately solve Maxwell's equations when an exact solution is impossible. These methods usually require a computer, and include the finite element method and finite-difference time-domain method. For more details, see Computational electromagnetics.

Overdetermination of Maxwell's equations

Maxwell's equations seem overdetermined, in that they involve six unknowns (the three components of E and B) but eight equations (one for each of the two Gauss's laws, three vector components each for Faraday's and Ampere's laws). (The currents and charges are not unknowns, being freely specifiable subject to charge conservation.) This is related to a certain limited kind of redundancy in Maxwell's equations: It can be proven that any system satisfying Faraday's law and Ampere's law automatically also satisfies the two Gauss's laws, as long as the system's initial condition does. This explanation was first introduced by Julius Adams Stratton in 1941. Although it is possible to simply ignore the two Gauss's laws in a numerical algorithm (apart from the initial conditions), the imperfect precision of the calculations can lead to ever-increasing violations of those laws. By introducing dummy variables characterizing these violations, the four equations become not overdetermined after all. The resulting formulation can lead to more accurate algorithms that take all four laws into account.

Limitations of the Maxwell equations as a theory of electromagnetism

While Maxwell's equations (along with the rest of classical electromagnetism) are extraordinarily successful at explaining and predicting a variety of phenomena, they are not exact, but approximations. In some special situations, they can be noticeably inaccurate. Examples include extremely strong fields (see Euler–Heisenberg Lagrangian) and extremely short distances (see vacuum polarization). Moreover, various phenomena occur in the world even though Maxwell's equations predict them to be impossible, such as "nonclassical light" and quantum entanglement of electromagnetic fields (see quantum optics). Finally, any phenomenon involving individual photons, such as the photoelectric effect, Planck's law, the Duane–Hunt law, single-photon light detectors, etc., would be difficult or impossible to explain if Maxwell's equations were exactly true, as Maxwell's equations do not involve photons. For the most accurate predictions in all situations, Maxwell's equations have been superseded by quantum electrodynamics.

Variations

Popular variations on the Maxwell equations as a classical theory of electromagnetic fields are relatively scarce because the standard equations have stood the test of time remarkably well.

Magnetic monopoles

Main article: Magnetic monopole

Maxwell's equations posit that there is electric charge, but no magnetic charge (also called magnetic monopoles), in the universe. Indeed, magnetic charge has never been observed (despite extensive searches) and may not exist. If they did exist, both Gauss's law for magnetism and Faraday's law would need to be modified, and the resulting four equations would be fully symmetric under the interchange of electric and magnetic fields.

See also

2

References

Further reading can be found in list of textbooks in electromagnetism

Historical publications

• On Faraday's Lines of Force – 1855/56 Maxwell's first paper (Part 1 & 2) – Compiled by Blaze Labs Research (PDF)
• On Physical Lines of Force – 1861 Maxwell's 1861 paper describing magnetic lines of Force – Predecessor to 1873 Treatise
• James Clerk Maxwell, "A Dynamical Theory of the Electromagnetic Field", Philosophical Transactions of the Royal Society of London 155, 459–512 (1865). (This article accompanied a December 8, 1864 presentation by Maxwell to the Royal Society.)
  • A Dynamical Theory Of The Electromagnetic Field – 1865 Maxwell's 1865 paper describing his 20 Equations, link from Google Books.
• J. Clerk Maxwell (1873) A Treatise on Electricity and Magnetism
  • Maxwell, J.C., A Treatise on Electricity And Magnetism – Volume 1 – 1873 – Posner Memorial Collection – Carnegie Mellon University
  • Maxwell, J.C., A Treatise on Electricity And Magnetism – Volume 2 – 1873 – Posner Memorial Collection – Carnegie Mellon University

The developments before relativity:

• Joseph Larmor (1897) "On a dynamical theory of the electric and luminiferous medium", Phil. Trans. Roy. Soc. 190, 205–300 (third and last in a series of papers with the same name).
• Hendrik Lorentz (1899) "Simplified theory of electrical and optical phenomena in moving systems", Proc. Acad. Science Amsterdam, I, 427–43.
• Hendrik Lorentz (1904) "Electromagnetic phenomena in a system moving with any velocity less than that of light", Proc. Acad. Science Amsterdam, IV, 669–78.
• Henri Poincaré (1900) "La théorie de Lorentz et le Principe de Réaction", Archives Néerlandaises, V, 253–78.
• Henri Poincaré (1902) La Science et l'Hypothèse
• Henri Poincaré (1905) "Sur la dynamique de l'électron", Comptes rendus de l'Académie des Sciences, 140, 1504–8.
• Catt, Walton and Davidson. "The History of Displacement Current". Wireless World, March 1979.

External links

• "Maxwell equations", Encyclopedia of Mathematics, EMS Press, 2001 [1994]
• maxwells-equations.com — An intuitive tutorial of Maxwell's equations.
• Mathematical aspects of Maxwell's equation are discussed on the Dispersive PDE Wiki.

Modern treatments

• Electromagnetism (ch. 11), B. Crowell, Fullerton College
• Lecture series: Relativity and electromagnetism, R. Fitzpatrick, University of Texas at Austin
• Electromagnetic waves from Maxwell's equations on Project PHYSNET.
• MIT Video Lecture Series (36 x 50 minute lectures) (in .mp4 format) – Electricity and Magnetism Taught by Professor Walter Lewin.

Other

• Feynman's derivation of Maxwell equations and extra dimensions
• Nature Milestones: Photons – Milestone 2 (1861) Maxwell's equations