Plane Waves in Homogeneous Isotropic Media

Every calculation in the optics of continuous media begins by writing the light as a plane wave, so this chapter opens by pinning down exactly what that object is and what a medium does to it. We take the vacuum plane wave — its wavefront geometry, its transversality k⋅E=0\mathbf{k}\cdot\mathbf{E}=0, and its energy flow — from aq-em-monochromatic-plane-waves-vacuum, and we take the source-free optical wave equation ∇2E=με ∂t2E\nabla^2\mathbf{E}=\mu\varepsilon\,\partial_t^2\mathbf{E} from aq-optics-maxwell-equations-dielectric-media. What this article owns is the in-medium extension: the single scalar constraint that survives, the dispersion relation k=nω/ck=n\omega/c, and the way the refractive index rescales every field ratio the vacuum wave carried. A monochromatic plane wave is a disturbance whose phase is constant over every plane perpendicular to a fixed direction k^\hat{\mathbf{k}}, sliding forward rigidly at one speed. In a homogeneous medium the parameters ε,μ\varepsilon,\mu do not depend on position, so the wave keeps its shape; in an isotropic medium they are scalars, so D∥E\mathbf{D}\parallel\mathbf{E} and there is no preferred direction. Under these two assumptions the plane-wave substitutions ∇→ik\nabla\to i\mathbf{k}, ∂t→−iω\partial_t\to-i\omega turn Maxwell's coupled differential equations into algebra, and the whole of propagation collapses onto k=nω/ck=n\omega/c with n=εrμrn=\sqrt{\varepsilon_r\mu_r}. Three payoffs organize the rest of the chapter. First, the medium enters only through nn: the crests slow to vp=c/nv_p=c/n and the wave impedance drops to η=η0/n\eta=\eta_0/n, so the same field amplitude carries a different intensity in glass than in vacuum. Second, when the index is complex, n~=n+iκ\tilde n=n+i\kappa, the identical plane-wave solution describes an attenuated wave whose amplitude decays as e−κωz/ce^{-\kappa\omega z/c} — the seed of absorption. Third, because n=n(ω)n=n(\omega) in every real material, plane waves of different colours travel at different speeds, which is the origin of dispersion, group velocity, and pulse spreading. The plane wave is the vehicle that converts microscopic material response into observable propagation.

Key takeaways
  • A monochromatic plane wave E=E0 ei(k⋅r−ωt)\mathbf{E}=\mathbf{E}_0\,e^{i(\mathbf{k}\cdot\mathbf{r}-\omega t)} is the fundamental solution of Maxwell's equations in a linear, homogeneous, isotropic medium; every field is a Fourier superposition of these.
  • Substituting it into the optical wave equation collapses everything to one scalar constraint, the dispersion relation k=nω/ck=n\omega/c with n=εrμr≈εrn=\sqrt{\varepsilon_r\mu_r}\approx\sqrt{\varepsilon_r} optically.
  • Transversality is not an extra postulate: ∇ ⁣⋅ ⁣D=0\nabla\!\cdot\!\mathbf{D}=0 forces k⋅E=0\mathbf{k}\cdot\mathbf{E}=0, and k^,E,H\hat{\mathbf{k}},\mathbf{E},\mathbf{H} form a right-handed orthogonal triad.
  • The medium rescales the vacuum ratios by nn: phase velocity vp=c/nv_p=c/n, wavelength λ=λ0/n\lambda=\lambda_0/n, and wave impedance η=μ/ε=η0/n\eta=\sqrt{\mu/\varepsilon}=\eta_0/n with η0=376.730 Ω\eta_0=376.730\ \Omega.
  • Intensity is I=⟨S⟩=12cε0n ∣E0∣2=vp⟨u⟩I=\langle S\rangle=\tfrac12 c\varepsilon_0 n\,|\mathbf{E}_0|^2=v_p\langle u\rangle: the same amplitude carries nn times more power in a denser medium, and electric and magnetic energies are equal.
  • A complex index n~=n+iκ\tilde n=n+i\kappa turns the same solution into a decaying wave, ∣E∣2∝e−αz|\mathbf{E}|^2\propto e^{-\alpha z} with α=4πκ/λ0\alpha=4\pi\kappa/\lambda_0 — the plane-wave root of Beer–Lambert.
  • Key idea: in a linear homogeneous isotropic medium the substitutions ∇→ik, ∂t→−iω\nabla\to i\mathbf{k},\ \partial_t\to-i\omega reduce Maxwell's equations to algebra — transversality plus the single dispersion relation k=nω/ck=n\omega/c — with the refractive index nn encoding the material entirely.
The core insight

In a linear homogeneous isotropic medium the plane-wave substitutions ∇→ik\nabla\to i\mathbf{k} and ∂t→−iω\partial_t\to-i\omega reduce Maxwell's equations to algebra: transversality k⋅E=0\mathbf{k}\cdot\mathbf{E}=0 plus the scalar dispersion relation k=nω/ck=n\omega/c, n=εrμrn=\sqrt{\varepsilon_r\mu_r}. Everything a medium does to a wave — slow it, attenuate it, disperse it — is carried by that one index.

Where this lives today

The plane-wave decomposition is the computational backbone of modern optics and photonics. Every thin-film coating design, waveguide-mode solver, and metasurface simulation expands the field in plane waves (the angular spectrum) and propagates each with its own eikzze^{ik_z z}; a single line, k=nω/ck=n\omega/c, is evaluated billions of times inside these codes. The wave impedance η=η0/n\eta=\eta_0/n is the quantity that anti-reflection stacks and fibre-connector index-matching gels are engineered to match, and impedance mismatch is why a bare glass surface reflects 4% per face. In laser metrology the intensity relation I=12cε0n∣E0∣2I=\tfrac12 c\varepsilon_0 n|E_0|^2 converts a measured power and spot size into the peak field a sample sees — the number that decides whether an experiment is in the linear or the nonlinear regime. And the complex-index promotion n~=n+iκ\tilde n=n+i\kappa is the entry point to everything that follows: absorption spectroscopy, ellipsometry (which reads nn and κ\kappa off the polarization of a reflected plane wave), and the dispersion engineering that governs ultrafast pulses and long-haul fibre.

The physical picture

Stand in a field of parallel wavefronts sweeping past you. Each front is an infinite flat sheet on which the phase is the same everywhere; the whole comb of sheets slides rigidly along k^\hat{\mathbf{k}} at speed c/nc/n, no sheet changing its value as it moves. The medium's only job is to slow that glide by the factor nn and, if it is lossy, to bleed a little amplitude out of each successive sheet. Because the medium is isotropic, the wave feels no preferred axis: pick any k^\hat{\mathbf{k}} and the same dispersion relation holds — which is why a block of glass looks identical whichever way you turn it, unlike a crystal whose index depends on direction. The algebraic magic is that a plane wave diagonalises the derivative operators. Maxwell's coupled partial differential equations become a linear algebra problem in k\mathbf{k} and ω\omega; demanding a non-trivial solution collapses to one scalar equation, k2=n2ω2/c2k^2=n^2\omega^2/c^2. Everything about propagation in a simple medium hides inside that relation. Transversality then falls out for free: Gauss's law ∇⋅D=0\nabla\cdot\mathbf{D}=0 becomes k⋅E=0\mathbf{k}\cdot\mathbf{E}=0, so E\mathbf{E} must lie in the wavefront plane, and Faraday's law completes a right-handed triad with the Poynting vector pointing along k^\hat{\mathbf{k}}.

Figure
The transverse-EM triad in a medium: E ⟂ H ⟂ k with the fields in phase, wavefronts (planes of constant phase) stacked ⟂ k and spaced by λ = λ₀/n, the whole comb gliding at v_p = c/n
A transverse plane wave: (k̂, E, H) is a right-handed orthogonal triad k̂ E H (into/out of page) wavefronts (constant phase) λ = λ₀/n

The second figure is the quantitative payoff a photonics engineer actually reaches for: what changes when the same field amplitude crosses from vacuum into a medium. The frequency is fixed by the source, but the crest speed, the wavelength, the impedance, and the delivered intensity all rescale by the single number nn.

Figure
Vacuum → glass (n = 1.50): the same field amplitude E₀ = 1 kV/m gives v_p = c/n = 1.999×10⁸ m/s, impedance η = η₀/n = 251.15 Ω, and intensity I = ½cε₀n|E₀|² = 1.99 kW/m² — the n-factor the vacuum wave does not carry
The punchline: a medium rescales every plane-wave quantity by n vacuum n = 1 glass n = 1.50 v_p = c = 3.00×10⁸ m/s η₀ = 376.73 Ω I₀ = 1.33 kW/m² v_p = c/n = 1.999×10⁸ m/s η = η₀/n = 251.15 Ω I = ½cε₀n|E₀|² = 1.99 kW/m² Same E₀ = 1 kV/m, same ν — only n differs. Intensity scales ∝ n (more power in the denser medium). λ compresses to λ₀/n; the impedance drops, so the H-field amplitude H₀ = E₀/η rises.

The core equations

A monochromatic plane wave is written with a complex exponential; physical fields are the real part. The convention ei(k⋅r−ωt)e^{i(\mathbf{k}\cdot\mathbf{r}-\omega t)} (a wave moving toward +k^+\hat{\mathbf{k}} as tt grows) is fixed once and used throughout — it sets the sign of κ\kappa and of every susceptibility in the chapter.

E(r,t)=E0 ei(k⋅r−ωt)(1)\mathbf{E}(\mathbf{r},t)=\mathbf{E}_0\,e^{i(\mathbf{k}\cdot\mathbf{r}-\omega t)} \qquad (1)

where: E0\mathbf{E}_0 — constant complex amplitude vector (V/m); k\mathbf{k} — wavevector, magnitude kk, along the propagation direction (rad/m); ω\omega — angular frequency (rad/s); r,t\mathbf{r},t — position and time. Surfaces of constant phase k⋅r−ωt=const\mathbf{k}\cdot\mathbf{r}-\omega t=\text{const} are planes ⊥k^\perp\hat{\mathbf{k}}. Inserting the ansatz into the optical wave equation leaves exactly one scalar constraint linking kk to ω\omega — the dispersion relation. All the material physics sits in the refractive index on its right-hand side.

k=n ωc,n=εrμr  → μr≈1   εr(2)k=n\,\frac{\omega}{c},\qquad n=\sqrt{\varepsilon_r\mu_r}\;\xrightarrow{\ \mu_r\approx1\ }\;\sqrt{\varepsilon_r} \qquad (2)

where: nn — refractive index; εr\varepsilon_r — relative permittivity at frequency ω\omega; μr≈1\mu_r\approx1 optically; c=1/ε0μ0c=1/\sqrt{\varepsilon_0\mu_0}. Because εr=εr(ω)\varepsilon_r=\varepsilon_r(\omega), equation (2) already contains dispersion: k(ω)k(\omega) is not simply proportional to ω\omega. Transversality is a consequence of the divergence equations, not an extra assumption, and it is where isotropy earns its keep.

k⋅E=0,k⋅H=0(3)\mathbf{k}\cdot\mathbf{E}=0,\qquad \mathbf{k}\cdot\mathbf{H}=0 \qquad (3)

where: H\mathbf{H} — magnetic field (A/m). ∇⋅D=0\nabla\cdot\mathbf{D}=0 gives ε k⋅E=0\varepsilon\,\mathbf{k}\cdot\mathbf{E}=0; a scalar ε\varepsilon (isotropy) then forces E⊥k\mathbf{E}\perp\mathbf{k}. In an anisotropic crystal only D⊥k\mathbf{D}\perp\mathbf{k} survives. Faraday's law ties H\mathbf{H} to E\mathbf{E} and fixes their ratio, the wave impedance, which is the vacuum value scaled down by nn.

H=nη0 k^×E,η=∣E∣∣H∣=με=η0n(3b)\mathbf{H}=\frac{n}{\eta_0}\,\hat{\mathbf{k}}\times\mathbf{E},\qquad \eta=\frac{|\mathbf{E}|}{|\mathbf{H}|}=\sqrt{\frac{\mu}{\varepsilon}}=\frac{\eta_0}{n} \qquad (3\text{b})

where: k^\hat{\mathbf{k}} — unit propagation vector; η0=μ0/ε0=376.730 Ω\eta_0=\sqrt{\mu_0/\varepsilon_0}=376.730\ \Omega — vacuum wave impedance; η\eta — impedance of the medium. (k^,E,H)(\hat{\mathbf{k}},\mathbf{E},\mathbf{H}) is right-handed. The phase velocity and wavelength follow directly from (2); the frequency is set at the source and is continuous across any interface.

vp=ωk=cn,λ=2πk=λ0n(4)v_p=\frac{\omega}{k}=\frac{c}{n},\qquad \lambda=\frac{2\pi}{k}=\frac{\lambda_0}{n} \qquad (4)

where: vpv_p — speed of surfaces of constant phase (m/s); λ0\lambda_0 — vacuum wavelength. Only the spatial scale and speed rescale by nn; ν=ω/2π\nu=\omega/2\pi does not. The time-averaged Poynting vector is the wave's intensity; carrying the nn from the impedance gives the working formula used throughout photonics.

I=⟨S⟩⋅k^=12 Re(E×H∗)⋅k^=∣E0∣22η=12 cε0n ∣E0∣2(5)I=\langle\mathbf{S}\rangle\cdot\hat{\mathbf{k}}=\tfrac12\,\mathrm{Re}(\mathbf{E}\times\mathbf{H}^*)\cdot\hat{\mathbf{k}}=\frac{|\mathbf{E}_0|^2}{2\eta}=\tfrac12\,c\varepsilon_0 n\,|\mathbf{E}_0|^2 \qquad (5)

where: II — intensity / irradiance (W/m²). The three forms are identical because 1/η=n/η0=cε0n1/\eta=n/\eta_0=c\varepsilon_0 n. Intensity scales ∝n\propto n for a fixed field amplitude. Finally, promoting the index to complex form turns the very same plane wave into an attenuated one — the single line that seeds the whole absorption story of this chapter.

n~=n+iκ  ⇒  E∝e−κωz/c ei(nωz/c−ωt),α=2κωc=4πκλ0(6)\tilde n=n+i\kappa\;\Rightarrow\; \mathbf{E}\propto e^{-\kappa\omega z/c}\,e^{i(n\omega z/c-\omega t)},\qquad \alpha=\frac{2\kappa\omega}{c}=\frac{4\pi\kappa}{\lambda_0} \qquad (6)

where: κ\kappa — extinction coefficient; α\alpha — intensity absorption coefficient (Beer–Lambert). The real part of n~\tilde n advances the phase; the imaginary part decays the amplitude. Both n~\tilde n (article 1487) and α\alpha (article 1488) are developed next.

Derivation AFrom the wave equation to the dispersion relation and transversality

The goal is to show that in a source-free, linear, homogeneous, isotropic medium the plane-wave ansatz reduces the four Maxwell equations to two statements — a single scalar dispersion relation k=nω/ck=n\omega/c and the transverse orthogonal triad — and that isotropy is exactly what makes E\mathbf{E} (not merely D\mathbf{D}) transverse. We insert the ansatz, read off the algebraic replacements, and assemble the field geometry, recalling the vacuum plane-wave object and the wave equation by name from their owners.

Step 1. Recall from aq-optics-maxwell-equations-dielectric-media the source-free optical wave equation ∇2E=με ∂t2E\nabla^2\mathbf{E}=\mu\varepsilon\,\partial_t^2\mathbf{E}, obtained by taking the curl of Faraday's law and using ∇⋅E=0\nabla\cdot\mathbf{E}=0. Why: this is the object the plane wave must solve; the material enters only through the product με\mu\varepsilon.

Step 2. For the ansatz (1), ∇→ik\nabla\to i\mathbf{k} and ∂t→−iω\partial_t\to-i\omega, so the wave equation becomes −k2E+με ω2E=0-k^2\mathbf{E}+\mu\varepsilon\,\omega^2\mathbf{E}=0, i.e. k2=με ω2=εrμr ω2/c2k^2=\mu\varepsilon\,\omega^2=\varepsilon_r\mu_r\,\omega^2/c^2. Taking the positive root defines n=εrμrn=\sqrt{\varepsilon_r\mu_r} and gives k=nω/ck=n\omega/c. Why: a plane wave diagonalises the derivatives, so a differential equation becomes one algebraic constraint — the dispersion relation carries all of propagation.

Step 3. Gauss's law ∇⋅D=0\nabla\cdot\mathbf{D}=0 becomes iε k⋅E=0i\varepsilon\,\mathbf{k}\cdot\mathbf{E}=0. Because ε\varepsilon is a scalar (isotropy) and non-zero, k⋅E=0\mathbf{k}\cdot\mathbf{E}=0: the electric field is transverse. Likewise ∇⋅B=0\nabla\cdot\mathbf{B}=0 gives k⋅H=0\mathbf{k}\cdot\mathbf{H}=0. Why: this is exactly the transversality recalled from aq-em-monochromatic-plane-waves-vacuum, but now the step k⋅D=0⇒k⋅E=0\mathbf{k}\cdot\mathbf{D}=0\Rightarrow\mathbf{k}\cdot\mathbf{E}=0 needs isotropy — in a crystal it fails, and the wave acquires a longitudinal E\mathbf{E} component (birefringence).

Step 4. Faraday's law ∇×E=−μ ∂tH\nabla\times\mathbf{E}=-\mu\,\partial_t\mathbf{H} becomes ik×E=iωμ Hi\mathbf{k}\times\mathbf{E}=i\omega\mu\,\mathbf{H}, so H=(k×E)/(ωμ)=(n/η0) k^×E\mathbf{H}=(\mathbf{k}\times\mathbf{E})/(\omega\mu)=(n/\eta_0)\,\hat{\mathbf{k}}\times\mathbf{E} using k=nω/ck=n\omega/c and η0=cμ0\eta_0=c\mu_0. Hence H⊥E\mathbf{H}\perp\mathbf{E}, the triad (k^,E,H)(\hat{\mathbf{k}},\mathbf{E},\mathbf{H}) is right-handed, and ∣E∣/∣H∣=η0/n=η|\mathbf{E}|/|\mathbf{H}|=\eta_0/n=\eta. Why: this fixes the magnetic amplitude and the impedance, the quantity whose mismatch at a boundary sets the Fresnel reflection.

Step 5. Ampère's law ∇×H=ε ∂tE\nabla\times\mathbf{H}=\varepsilon\,\partial_t\mathbf{E} reproduces k2=με ω2k^2=\mu\varepsilon\,\omega^2, confirming the solution is complete and self-consistent — no independent constraint is left. Why: the fourth equation must not over-determine the wave; that it returns the dispersion relation is the consistency check. *Verified: k=nω/ck=n\omega/c from με=n2/c2\mu\varepsilon=n^2/c^2; k⋅E=0\mathbf{k}\cdot\mathbf{E}=0 symbolically; ∣E∣/∣H∣=η0/n|\mathbf{E}|/|\mathbf{H}|=\eta_0/n; η0=376.730 Ω\eta_0=376.730\ \Omega in verify_aq-optics-plane-waves-isotropic-media.py.*

Go deeper · Derivation B · 1 minIntensity and the equal sharing of energy: I = ½cε₀n|E₀|² = v_p⟨u⟩
Derivation BIntensity and the equal sharing of energy: I = ½cε₀n|E₀|² = v_p⟨u⟩

The aim is to compute the time-averaged intensity of the plane wave from its fields, show that the electric and magnetic energy densities are equal, and prove that energy is transported along k^\hat{\mathbf{k}} at exactly the phase velocity in a lossless non-dispersive medium. We form the mean Poynting vector, add the two energy densities, and take their ratio — the origin of the working intensity formula I=12cε0n∣E0∣2I=\tfrac12 c\varepsilon_0 n|\mathbf{E}_0|^2.

Step 1. With H=(n/η0)k^×E\mathbf{H}=(n/\eta_0)\hat{\mathbf{k}}\times\mathbf{E} and k^×(k^×E)=−E\hat{\mathbf{k}}\times(\hat{\mathbf{k}}\times\mathbf{E})=-\mathbf{E} for transverse E\mathbf{E}, the mean flux is ⟨S⟩=12Re(E×H∗)=12(n/η0)∣E0∣2 k^\langle\mathbf{S}\rangle=\tfrac12\mathrm{Re}(\mathbf{E}\times\mathbf{H}^*)=\tfrac12(n/\eta_0)|\mathbf{E}_0|^2\,\hat{\mathbf{k}}. Since n/η0=cε0nn/\eta_0=c\varepsilon_0 n, the intensity is I=12cε0n∣E0∣2I=\tfrac12 c\varepsilon_0 n|\mathbf{E}_0|^2. Why: this is the number a power meter reads; the factor nn is what the vacuum owner's I=12cε0∣E0∣2I=\tfrac12 c\varepsilon_0|\mathbf{E}_0|^2 does not carry.

Step 2. The time-averaged energy density is ⟨u⟩=14ε∣E0∣2+14μ∣H0∣2\langle u\rangle=\tfrac14\varepsilon|\mathbf{E}_0|^2+\tfrac14\mu|\mathbf{H}_0|^2. Using ∣H0∣=(n/η0)∣E0∣|\mathbf{H}_0|=(n/\eta_0)|\mathbf{E}_0| and μ(n/η0)2=μεrμr(ε0/μ0)=ε\mu (n/\eta_0)^2=\mu\varepsilon_r\mu_r(\varepsilon_0/\mu_0)=\varepsilon (for μr≈1\mu_r\approx1), the magnetic term equals the electric term. Why: equipartition of electric and magnetic energy is a defining feature of a freely propagating wave and makes the bookkeeping collapse to one term.

Step 3. Hence ⟨u⟩=12ε∣E0∣2\langle u\rangle=\tfrac12\varepsilon|\mathbf{E}_0|^2, and the energy-transport velocity is ⟨S⟩/⟨u⟩=(12cε0n∣E0∣2)/(12ε0εr∣E0∣2)=cn/εr=c/n=vp\langle S\rangle/\langle u\rangle=(\tfrac12 c\varepsilon_0 n|\mathbf{E}_0|^2)/(\tfrac12\varepsilon_0\varepsilon_r|\mathbf{E}_0|^2)=cn/\varepsilon_r=c/n=v_p. Why: in a lossless non-dispersive medium the energy rides with the crests; only when nn disperses does the energy velocity depart from vpv_p toward the group velocity (article 1491). *Verified: I=12cε0n∣E0∣2I=\tfrac12 c\varepsilon_0 n|\mathbf{E}_0|^2 (glass n=1.5, E₀=1 kV/m → 1.99 kW/m²), electric = magnetic energy, and ⟨S⟩/⟨u⟩=c/n\langle S\rangle/\langle u\rangle=c/n in verify_aq-optics-plane-waves-isotropic-media.py.*

Limiting cases

CaseResultMeaning
Vacuumn=1, k=ω/c, vp=c, η=376.73 Ωn=1,\ k=\omega/c,\ v_p=c,\ \eta=376.73\ \Omegano medium to slow or load the wave; the vacuum impedance is universal
Transparent dielectricn=εr∈R, vp=c/nn=\sqrt{\varepsilon_r}\in\mathbb{R},\ v_p=c/nreal index, no attenuation; E\mathbf{E} and H\mathbf{H} oscillate exactly in phase (glass, n≈1.5n\approx1.5)
Non-magnetic optical mediumμr≈1⇒n≈εr, η=η0/n\mu_r\approx1\Rightarrow n\approx\sqrt{\varepsilon_r},\ \eta=\eta_0/nmagnetic response vanishes optically, so the index comes entirely from εr\varepsilon_r
Absorbing mediumn~=n+iκ, ∣E∣2∝e−αz\tilde n=n+i\kappa,\ |\mathbf{E}|^2\propto e^{-\alpha z}complex index: phase advances with nn, amplitude decays with κ\kappa, α=4πκ/λ0\alpha=4\pi\kappa/\lambda_0
Anisotropic crystalk⋅D=0\mathbf{k}\cdot\mathbf{D}=0 but k⋅E≠0\mathbf{k}\cdot\mathbf{E}\neq0isotropy lost: E\mathbf{E} is no longer purely transverse; two indices per direction (birefringence)

Worked examples

1 · Intensity carried by a plane wave in glass (confidence builder)

A monochromatic beam has electric-field amplitude E0=1.0 kV/mE_0=1.0\ \mathrm{kV/m} inside crown glass of index n=1.50n=1.50. Find the intensity it carries, and compare with the amplitude that would deliver the solar constant 1361 W/m21361\ \mathrm{W/m^2} in vacuum.

Step 1. Use I=12cε0n∣E0∣2I=\tfrac12 c\varepsilon_0 n|E_0|^2. Insert c=2.998×108c=2.998\times10^8, ε0=8.854×10−12\varepsilon_0=8.854\times10^{-12}, n=1.50n=1.50, E0=103E_0=10^3: I=12(2.998×108)(8.854×10−12)(1.50)(106)I=\tfrac12(2.998\times10^8)(8.854\times10^{-12})(1.50)(10^6). Grouping, 12cε0=1.327×10−3\tfrac12 c\varepsilon_0=1.327\times10^{-3}, times n∣E0∣2=1.50×106n|E_0|^2=1.50\times10^6 gives I=1.99 kW/m2I=\mathbf{1.99\ kW/m^2}.

Step 2. The same field in vacuum (n=1n=1) would carry I0=1.99/1.50=1.33 kW/m2I_0=1.99/1.50=1.33\ \mathrm{kW/m^2} — the intensity is larger in glass by the factor nn, because the medium loads the wave (lower impedance ⇒\Rightarrow larger H0H_0 for the same E0E_0).

Step 3. For the solar constant in vacuum, invert I=12cε0∣E0∣2I=\tfrac12 c\varepsilon_0|E_0|^2: E0=2I/(cε0)=2(1361)/(2.998×108⋅8.854×10−12)=1013 V/mE_0=\sqrt{2I/(c\varepsilon_0)}=\sqrt{2(1361)/(2.998\times10^8\cdot8.854\times10^{-12})}=\mathbf{1013\ V/m}. So bright sunlight is an optical field of only about a kilovolt per metre — a useful yardstick. I = 1.99 kW/m² in glass; sunlight ↔ E₀ = 1013 V/m.

2 · Impedance and the magnetic-field amplitude

For the same wave (E0=1.0 kV/mE_0=1.0\ \mathrm{kV/m}, n=1.50n=1.50), find the wave impedance of the glass and the magnetic-field amplitude H0H_0, and confirm the intensity via I=12ηH02I=\tfrac12\eta H_0^2.

Step 1. Impedance: η=η0/n=376.730/1.50=251.15 Ω\eta=\eta_0/n=376.730/1.50=\mathbf{251.15\ \Omega}. The medium presents a lower impedance than vacuum, exactly as a denser transmission line does.

Step 2. Magnetic amplitude: H0=E0/η=103/251.15=3.982 A/mH_0=E_0/\eta=10^3/251.15=\mathbf{3.982\ A/m} (equivalently B0=μ0H0=5.00 μTB_0=\mu_0 H_0=5.00\ \mu\mathrm{T}). Because η\eta dropped, H0H_0 is larger than the vacuum value E0/η0=2.654 A/mE_0/\eta_0=2.654\ \mathrm{A/m}.

Step 3. Cross-check the intensity: I=12ηH02=12(251.15)(3.982)2=1.99×103 W/m2I=\tfrac12\eta H_0^2=\tfrac12(251.15)(3.982)^2=1.99\times10^3\ \mathrm{W/m^2}, matching Example 1. The three equivalent forms I=∣E0∣2/2η=12ηH02=12cε0n∣E0∣2I=|E_0|^2/2\eta=\tfrac12\eta H_0^2=\tfrac12 c\varepsilon_0 n|E_0|^2 are one statement. η = 251.15 Ω, H₀ = 3.98 A/m, I = 1.99 kW/m² (consistent).

Physics of the portrait

The plane wave is the atom of wave optics: not a real beam (which is finite and therefore a superposition of plane waves) but the exact eigenmode from which every real field is built by Fourier synthesis. Its power is that it turns partial differential equations into algebra, so that a medium's entire influence is compressed into one number, nn.

Origins · 1 minplane wave and how this was discovered
PW
plane wave
the monochromatic eigenmode E₀e^{i(k·r−ωt)} — the building block of wave optics · electromagnetic optics of media

The plane wave is less a discovery than a language. Its roots lie in the wave equation that d'Alembert wrote for a vibrating string in 1747, whose solutions are rigidly translating profiles f(x∓vt)f(x\mp vt); the sinusoidal, single-frequency member of that family became the workhorse once Fourier showed in 1822 that any disturbance is a superposition of sinusoids. When Maxwell identified light as an electromagnetic wave in 1865 and computed its speed from ε0,μ0\varepsilon_0,\mu_0, the transverse plane wave became the concrete object that carried the identification: E\mathbf{E} and H\mathbf{H} perpendicular to each other and to the propagation, in phase, with the ratio η0=μ0/ε0\eta_0=\sqrt{\mu_0/\varepsilon_0}. John Henry Poynting (1884) and Oliver Heaviside gave the energy flux S=E×H\mathbf{S}=\mathbf{E}\times\mathbf{H} its meaning, so that the plane wave's intensity became a measurable transport of energy rather than a formal quantity. The idealization is deliberately unphysical — an infinite sheet of infinite extent and duration — yet it is controlled: a real beam of waist w0w_0 behaves as a plane wave over its Rayleigh range, and a laser's coherence time can exceed a second, so within those windows every formula here applies with negligible error. Everything in this chapter — refractive index, absorption, phase and group velocity, dispersion, Kramers–Kronig — is bookkeeping performed on this one object as it meets matter.

History
1747
— d'Alembert's wave equation for the vibrating string: rigidly translating solutions
1822
— Fourier establishes that any disturbance is a superposition of sinusoids (plane-wave basis)
1865
— Maxwell identifies light as an electromagnetic wave and computes its speed from ε₀, μ₀
1884
— Poynting formulates the energy-flux vector S = E × H
1885
— Heaviside independently arrives at the E × H energy flux, confirming the transport law

Going further

Texts. Hecht, Optics, 5th ed., §3.2–3.5 · Saleh & Teich, Fundamentals of Photonics, 3rd ed., Ch. 5 (electromagnetic optics) · Loudon, The Quantum Theory of Light, 3rd ed., Ch. 1 · Foot, Atomic Physics, §7.1 · Jackson, Classical Electrodynamics, 3rd ed., §7.1.

Provenance. The vacuum plane wave, its transversality k⋅E=0\mathbf{k}\cdot\mathbf{E}=0, and its energy flux are recalled by name from aq-em-monochromatic-plane-waves-vacuum; the source-free optical wave equation ∇2E=με ∂t2E\nabla^2\mathbf{E}=\mu\varepsilon\,\partial_t^2\mathbf{E} from aq-optics-maxwell-equations-dielectric-media. A verifier should confirm: (1) k=nω/ck=n\omega/c from με=n2/c2\mu\varepsilon=n^2/c^2 (Jackson §7.1, Hecht §3.2); (2) the field ratio ∣E∣/∣H∣=η=η0/n|E|/|H|=\eta=\eta_0/n with η0=376.730 Ω\eta_0=376.730\ \Omega; (3) the intensity I=12cε0n∣E0∣2I=\tfrac12 c\varepsilon_0 n|E_0|^2 and its equality with vp⟨u⟩v_p\langle u\rangle via electric–magnetic equipartition (Saleh & Teich Ch. 5); (4) α=4πκ/λ0\alpha=4\pi\kappa/\lambda_0 from the complex-index propagation factor (Loudon Ch. 1). The worked numbers (glass n=1.5n=1.5: vp=1.999×108v_p=1.999\times10^8 m/s, η=251.15 Ω\eta=251.15\ \Omega, I=1.99I=1.99 kW/m²; sunlight E0=1013E_0=1013 V/m) are recomputable from the formulas.

Prerequisites. Maxwell's Equations in Dielectric Media · Dielectric Susceptibility and Permittivity · Electric Polarization and the Displacement Field. Leads to. The Refractive Index · Phase Velocity · Complex Refractive Index and the Extinction Coefficient · Group Velocity and Group Index.

Is a plane wave a real thing, or just a mathematical convenience?

It is an idealization — an infinite sheet of infinite duration that no source can produce — but a controlled one. Any real field is an exact superposition of plane waves (its angular spectrum in space, its frequency spectrum in time), so nothing is lost by working with one plane wave at a time and adding them up. A collimated laser beam of waist w0w_0 behaves as a single plane wave over its Rayleigh range zR=πw02/λz_R=\pi w_0^2/\lambda — metres for a millimetric visible beam — and over the coherence length set by its linewidth. Within that window every formula in this article applies to better than a part in 10310^3. The corrections (finite width →\to diffraction, finite duration →\to linewidth) are themselves organized as small departures from the plane wave, which is why it is the right starting point rather than a crude one.

The frequency doesn't change entering a medium but the wavelength does — why?

The frequency is imposed by the source and by the boundary: at an interface the fields on both sides must oscillate in step at every instant, which is only possible if they share the same ω\omega. What the medium changes is the speed of the crests, vp=c/nv_p=c/n. Since vp=νλv_p=\nu\lambda with ν\nu fixed, the wavelength must shrink to λ=λ0/n\lambda=\lambda_0/n — the crests pack closer together while ticking past at the same rate. This is exactly why colour (a frequency) is preserved when light enters water or glass, even though the wave physically slows and compresses.

How can the intensity be larger in glass for the same E₀ — where does the extra power come from?

No extra power is created; the comparison just holds a different quantity fixed. Intensity is I=12cε0n∣E0∣2I=\tfrac12 c\varepsilon_0 n|E_0|^2, so for a fixed field amplitude it rises with nn. The physical reason is the impedance: η=η0/n\eta=\eta_0/n is lower in glass, so the same E0E_0 is accompanied by a larger magnetic amplitude H0=E0/ηH_0=E_0/\eta, and the energy flux 12E0H0\tfrac12 E_0 H_0 grows. If instead you held the power fixed (as when a beam refracts into glass), the field amplitude would adjust and there is no paradox. In practice, what is continuous across a normal interface is the tangential EE and HH, and matching them gives the Fresnel transmission — a bookkeeping this article's impedance sets up.

Why does isotropy matter so much for transversality?

Because the fundamental divergence law is about D\mathbf{D}, not E\mathbf{E}: ∇⋅D=0\nabla\cdot\mathbf{D}=0 gives k⋅D=0\mathbf{k}\cdot\mathbf{D}=0 always. Whether this makes E\mathbf{E} transverse depends on how D\mathbf{D} relates to E\mathbf{E}. In an isotropic medium D=εE\mathbf{D}=\varepsilon\mathbf{E} with a scalar ε\varepsilon, so D∥E\mathbf{D}\parallel\mathbf{E} and k⋅E=0\mathbf{k}\cdot\mathbf{E}=0 follows. In a crystal ε\varepsilon is a tensor, D\mathbf{D} and E\mathbf{E} point in different directions, and k⋅E\mathbf{k}\cdot\mathbf{E} is generally non-zero: the wave has a longitudinal electric component, the energy flow (Poynting vector) walks off from k\mathbf{k}, and each propagation direction supports two waves with two indices. Birefringence, wave plates, and phase-matching in nonlinear crystals all live in that broken isotropy.

What is the wave impedance really, and why do engineers care about it?

It is the ratio of the electric to the magnetic field amplitude in the wave, η=∣E∣/∣H∣=μ/ε=η0/n\eta=|E|/|H|=\sqrt{\mu/\varepsilon}=\eta_0/n, and it plays exactly the role that characteristic impedance plays on a transmission line. When a wave meets a boundary between two media, the reflection is governed by the mismatch of their impedances: equal impedances ⇒\Rightarrow no reflection, large mismatch ⇒\Rightarrow strong reflection. That is why a single-layer anti-reflection coating uses an index nf=n0nsn_f=\sqrt{n_0 n_s} (a geometric-mean impedance match) and why index-matching fluid on a fibre connector kills the back-reflection. The 4% reflection from each face of a bare glass plate is just the impedance step from η0\eta_0 to η0/1.5\eta_0/1.5.

Does energy always travel at the phase velocity c/n?

In a lossless, non-dispersive medium, yes: the derivation above gives ⟨S⟩/⟨u⟩=vp=c/n\langle S\rangle/\langle u\rangle=v_p=c/n. But that relied on nn being constant. When n=n(ω)n=n(\omega) — which is always true in a real material — a pulse (a superposition of frequencies) has its envelope, and hence its energy and information, move at the group velocity vg=c/ngv_g=c/n_g, generally different from vpv_p. Near an absorption line the two can differ dramatically, and vpv_p can even exceed cc without anything physical exceeding cc. So "energy travels at c/nc/n" is the correct non-dispersive statement; the dispersive correction is the subject of the phase-velocity and group-velocity articles that follow.

Before this

Where this leads