03 · Intermediate → advanced

Maxwell equations

Assemble all four vacuum Maxwell equations in local and fixed-surface integral form, keeping enclosed charge and signed current distinct; derive source-free wave equations and transverse plane-wave constraints, show the finite vacuum wave speed and correctly orient E, B and propagation. State boundary, material and quasistatic limits before applying simplified field laws.

270–330 minIntermediate → advanced
In this lesson
  • Write all four vacuum Maxwell equations in SI, distinguishing ρ and J from fields E and B.
  • Relate differential divergence equations to fixed-volume flux equations by Gauss's theorem.
  • Relate differential curl equations to fixed-surface circulation equations by Stokes with oriented surfaces.
  • Explain the difference between ∇·B=0 (closed magnetic flux) and ∇×E=−∂B/∂t (induced electric circulation).
  • State Maxwell–Ampère's electric-flux correction and why a changing electric field can generate B even where J=0.
  • Derive source-free vacuum wave equations for E and B using double-curl identities and the transverse divergence constraints.
  • Deduce c=1/√(μ₀ε₀) and check its dimensions; distinguish vacuum c from phase velocity in an arbitrary dispersive material.
  • Derive k·E₀=k·B₀=0 and B₀=(1/ω)k×E₀ for a monochromatic plane-wave convention.
  • State the standard limiting interface conditions with surface charge/current and specify absence of singular sheet time derivatives and proper orientations.
  • Solve six worked problems on flux/circulation signs, wave propagation, material-law scope and boundary conditions.

The electrostatic and magnetostatic source laws are not independent recipes that remain unchanged at high frequency. Faraday induction and displacement current complete a coupled four-equation description. We write Maxwell's equations in vacuum with SI units and a FIXED observation frame, derive both integral laws and wave propagation, and separate exact equations from approximate device formulas.

0103.1 — Four LOCAL Maxwell equations in vacuum

Let ρ(r,t) denote total electric charge density and J(r,t) total electric current density in the chosen vacuum-source formulation. The fields E and B obey four equations, the first two specifying divergence and the latter two specifying curl. An electric charge produces electric flux; magnetic flux has no monopole source in the classical theory; varying B produces electric circulation; J and varying E produce magnetic circulation.

A pointwise source-free region means ρ=0 and J=0 THERE, not that E and B must vanish. Waves pass through empty space because the time derivatives couple the fields. Under ordinary vacuum conditions the four equations are local laws, whereas simple static Coulomb/Biot–Savart formulas have additional source/time assumptions.

∇·E=ρ/ε0; ∇·B=0
∇×E=−∂B/∂t
∇×B=μ0J+μ0ε0∂E/∂t

0203.2 — Electric and magnetic GAUSS laws in integral form

Gauss: a closed surface and its enclosed chargePerspective projection of a closed spherical Gaussian surface S around an enclosed positive charge +q. Outward normals point away from the surface. A negative charge outside the sphere is not enclosed and does not contribute to net flux.+qS
n̂
−q
qext
Qenc = +q
Gauss: a closed surface and its enclosed charge

The dashed equator is the BACK half of a closed spherical surface, not an open circular contour. Only +q is enclosed: ∯s E·n dA = +q/ε0. An outside charge can alter E locally but not the enclosed charge.

Choose a fixed volume V bounded by a closed outward-oriented surface ∂V. Apply the divergence theorem ∮∂V F·n̂dS=∫V∇·F d³r. Electric Gauss gives the total signed electric flux equal to Qinside/ε₀, where Qinside=∫Vρd³r.

For B the closed flux is exactly zero. This does NOT imply that B vanishes at the surface; inward and outward contributions can cancel. The surface must be closed to use these divergence statements, and its normal is outward by convention.

∮∂V E·n̂dS=Qinside/ε0; Qinside=∫Vρd3r
∮∂V B·n̂dS=0
[∮E·dS]=V·m; [∮B·dS]=Wb

0303.3 — Faraday and Ampère–Maxwell as oriented circulations

Ampère: a current creates a circulating magnetic fieldPositive steady current I along positive z, toward the viewer, marked by a dot in a circle. Viewed from positive z, the magnetic field B and H circulate counterclockwise. A concentric dashed Ampere contour C is oriented counterclockwise with normal positive z.I (+z) ⊙BφC
Bφ, Hφ
B = μ0H
∮C H·dℓ = Ienc
+z ⊙
Ampère: a current creates a circulating magnetic field

For a long straight steady current I>0 out of the page, the right-hand rule gives B counterclockwise (upward at the right-hand side of the wire). In vacuum, B=μ0H, and the oriented contour obeys ∮C H·dℓ=Ienc. This is Ampère's law, NOT Faraday induction.

Choose a FIXED smooth open surface S with oriented closed boundary C. Its positive contour direction follows the right-hand rule about n̂. Stokes' theorem converts the local curl equations to two signed circulation equations; the time derivatives act on flux through that fixed S.

The minus sign in Faraday makes the induced E circulation oppose the INCREASE of a positive B flux for a passive loop. In Ampère–Maxwell the electric-flux term permits nonzero magnetic circulation with no conduction current through S. For a moving material contour, the emf contains u×B as developed separately in EM-03.01; do not apply the fixed-loop electric-circulation equation alone to a moving wire.

∮C E·dℓ=−d/dt∫S B·n̂dS
∮C B·dℓ=μ0∫S J·n̂dS+μ0ε0d/dt∫S E·n̂dS
C=∂S; C direction ↔ n̂ via right-hand rule

0403.4 — A consistency check: Maxwell implies charge continuity

Apply the divergence operator to Maxwell–Ampère. Since ∇·(∇×B)=0, one obtains 0=μ₀∇·J+μ₀ε₀∂(∇·E)/∂t. Insert Gauss's law ∇·E=ρ/ε₀ to derive ∂ρ/∂t+∇·J=0.

This shows the four local field equations are compatible with local charge conservation; it does not replace the physical statement that charge is conserved. The detailed interpretation of integrated current at a moving boundary belongs in the following continuity chapter.

0=∇·(∇×B)=μ0[∇·J+ε0∂(∇·E)/∂t]
∂ρ/∂t+∇·J=0

0503.5 — Source-free regions: derive the ELECTRIC wave equation

In a vacuum region with ρ=0 and J=0, Maxwell gives ∇·E=0 and ∇×B=μ₀ε₀∂E/∂t. Take the curl of ∇×E=−∂B/∂t. Commute the spatial curl with the time derivative for sufficiently regular fields, and substitute Maxwell–Ampère.

The double-curl identity ∇×(∇×E)=∇(∇·E)−∇²E reduces to −∇²E because there is no local charge. Thus E obeys the three-component homogeneous wave equation with coefficient μ₀ε₀. The elimination of ∇(∇·E) is NOT valid as written where ρ varies spatially.

∇×(∇×E)=−∂(∇×B)/∂t=−μ0ε0∂²E/∂t2
∇×(∇×E)=∇(∇·E)−∇²E=−∇²E (ρ=0)
∇²E−μ0ε0∂²E/∂t2=0

0603.6 — Magnetic wave equation and speed in vacuum

Take the curl of ∇×B=μ₀ε₀∂E/∂t in the same source-free region, substitute ∇×E=−∂B/∂t, and use ∇·B=0 in the double-curl identity. The magnetic field satisfies the same wave equation.

Comparing with ∇²F−(1/v²)∂²F/∂t²=0 yields v=c=1/√(μ₀ε₀). In modern SI the vacuum light speed is fixed at c=299 792 458 m/s, while the measured μ₀ and ε₀ have consistent values. This c is the vacuum wave speed, not an automatic phase or group velocity in a dispersive material.

∇²B−μ0ε0∂²B/∂t2=0
c=1/√(μ0ε0)=299 792 458 m/s (vacuum)
[μ0ε0]=s2/m2

0703.7 — Plane-wave geometry: transverse E and B

Traveling wave: E, B and propagation form a right-handed triadAt a single position and time in a source-free vacuum plane wave, propagation vector k and energy flux S point right along positive x. Electric field E points up along positive y and magnetic field B points out of page along positive z, marked with a dot. The three field vectors are mutually perpendicular.B : +z ⊙
E : +y ↑
k, S : +x →
B = k×E/ω
S = E×B/μ0
E ⊥ B ⊥ k
Traveling wave: E, B and propagation form a right-handed triad

For the convention cos(kx−ωt), choose k∥+x, E∥+y and B∥+z (out of page). Maxwell–Faraday gives B=(1/ω)k×E, so E×B∥+x and S=E×B/μ0 travels with the wave. The arrows depict one instant when E>0, not a longitudinal electric field or a spatial wave profile. In vacuum |B0|=|E0|/c.

Use the real part of complex fields E(r,t)=E₀exp[i(k·r−ωt)] and B(r,t)=B₀exp[i(k·r−ωt)] with positive ω, source-free vacuum and nonzero k. The divergence laws give k·E₀=0 and k·B₀=0. Faraday gives ik×E₀=iωB₀, so B₀=(1/ω)k×E₀.

For |k|=ω/c, the magnetic amplitude magnitude is |B₀|=|E₀|/c. If a wave propagates toward +z with E toward +x, then B points +y: x̂×ŷ=ẑ. Reversing k reverses B for the same E orientation. This geometry follows from the SIGNED complex convention, not an unsigned arrow mnemonic.

k·E0=0; k·B0=0; B0=(1/ω)k×E0
|k|=ω/c ⇒ |B0|=|E0|/c
k∥+z, E∥+x ⇒ B∥+y

0803.8 — Interface conditions: identify sheets and limits

Let n̂ point from medium/region 1 to region 2. For ordinary finite fields near an interface, a thin Gaussian pillbox and shrinking Stokes loop give the jumps of normal E, normal B, tangential E and tangential B. In a vacuum formulation with TOTAL surface charge σ and surface current K (with no singular sheet time-derivative fields), the jumps below hold with the chosen n̂.

The tangential B jump is a VECTOR cross-product condition, not a freely assignable scalar sign. Across material interfaces one often instead uses D and H with FREE σ and K and constitutive assumptions. When a singular sheet has explicit surface polarization/magnetization or time-dependent generalized sheet responses, additional effective interface terms may enter; do not blindly apply ordinary textbook jumps.

n̂·(E2−E1)=σtotal/ε0; n̂·(B2−B1)=0
n̂×(E2−E1)=0; n̂×(B2−B1)=μ0Ktotal
materials (ordinary): n̂·(D2−D1)=σfree; n̂×(H2−H1)=Kfree

0903.9 — Worked vacuum-wave and regime check

At f=1.00 THz, λ₀=c/f≈2.998×10⁻⁴ m=300 μm in vacuum. For a monochromatic plane wave with electric amplitude E₀=300 V/m, the magnetic amplitude is B₀=E₀/c≈1.00×10⁻⁶ T=1.00 μT. These are paired field amplitudes of the SAME wave, not separately chosen source fields.

This vacuum calculation does not predict a wavelength inside an arbitrary conductor, graphene device or dispersive dielectric. The corresponding material propagation depends on the frequency-dependent constitutive response and boundaries. An apparently small physical object can still violate an instantaneous quasistatic approximation when relevant delay, transport or resonance times matter.

f=1.00 THz ⇒ λ0=c/f≈300 μm (vacuum)
E0=300 V/m ⇒ B0=E0/c≈1.00 μT

03.10 — Six mastery exercises with worked solutions

  1. State all four vacuum Maxwell equations in local form, identify the two divergence constraints and their electric/magnetic surface flux counterparts.
  2. For a FIXED +z-oriented loop enclosing a uniform increasing +z magnetic flux, determine the sign of ∮E·dℓ relative to counterclockwise traversal; explain whether a global single-valued electrostatic potential alone can represent E.
  3. Starting from source-free Maxwell, derive ∇²E−μ₀ε₀∂²E/∂t²=0, clearly identifying the assumption needed to eliminate ∇(∇·E).
  4. For a vacuum plane wave traveling in +z with E₀=300 x̂ V/m and f=1.00 THz, determine B₀ vector, wavelength and units.
  5. At a regular interface with normal from 1 to 2 and surface current K=K x̂, compute n̂×(B₂−B₁) when n̂=ẑ; describe why a scalar B₂−B₁=μ₀K is not the general boundary equation.
  6. Derive continuity from vacuum Ampère–Maxwell and Gauss, and explain why J=0 in an empty region does not force ∇×B=0 during wave propagation.
Worked solutions

1. ∇·E=ρ/ε₀, ∇·B=0, ∇×E=−∂tB, ∇×B=μ₀J+μ₀ε₀∂tE. Their closed-surface forms are ∮E·dS=Qinside/ε₀ and ∮B·dS=0, with outward normal.

2. For +z flux increasing, Faraday gives ∮E·dℓ=−dΦB/dt<0: clockwise electric circulation viewed from +z. A global single-valued static V would give ∮(−∇V)·dℓ=0, so it is insufficient.

3. Curl Faraday: ∇×∇×E=−μ₀ε₀∂t²E. With ρ=0, ∇·E=0, so the left side is −∇²E. Rearranging gives the vacuum electric wave equation.

4. B₀=(1/c) ẑ×(300 x̂ V/m)=+(300/c)ŷ T≈+1.00 μT ŷ. λ₀=c/f≈2.998×10⁻⁴ m≈300 μm.

5. With n̂=ẑ and K=Kx̂, the jump condition reads ẑ×(B₂−B₁)=μ₀Kx̂. This requires a tangential jump component ΔB_y=−μ₀K (other components depend on boundary/source constraints), not a scalar unsigned jump.

6. Divergence of ∇×B=μ₀J+μ₀ε₀∂tE gives 0=μ₀(∇·J+∂tρ). In vacuum J=0 yet ∇×B=μ₀ε₀∂tE can remain nonzero in a propagating wave.

EM-03.04 will develop charge continuity for fixed and moving control volumes, derive its integral flux form and connect interface-charge accumulation to measurable circuit and displacement-current balances.