← ArticlesGeometric and Wave Optics: /20 SummaryEngineering · PhysicsLesson 4/13← PrevNext →
GuidePublished 14 Aug 202627 min readBy Kevin JoginPhysicsApplied Classical PhysicsGeometric and Wave OpticsContext and scope

Engineering · Physics · Applied Classical Physics

Geometric and Wave Optics: /20 Summary

Engineering handbook for geometric and wave optics, covering context and scope, /20 summary, geometric optics & the eikonal approximation.

Executive summary

This handbook section converts the supplied engineering material into a practical, source-controlled reference. It concentrates on the following learning outcomes.

Context and scope
/20 Summary
Geometric Optics & the Eikonal Approximation
From Wave Equation to Eikonal
Rays as Trajectories
Fermat's Principle

Context and scope

the supplied physics reference, Applications of Classical Physics — Chapters 7–10 Geometric Optics · Diffraction · Interferometry & Coherence · Nonlinear Optics



/20 Summary

  • Three regimes of wave propagation, distinguished by the ratio of wavelength λ\lambda to relevant length scale LL:
    • λL\lambda \ll LGeometric (ray) optics (eikonal limit, WKB)
    • λL\lambda \sim LDiffraction / wave optics (Helmholtz, Fresnel-Kirchhoff)
    • λL\lambda \gg L ⇒ Lumped/long-wavelength limit (not central to optics)
  • Three categories of phenomena the chapter develops:
    1. Propagation (geometric optics, eikonal, paraxial)
    2. Coherent superposition (diffraction, interference, holography)
    3. Statistical superposition (coherence theory, V-C-Z theorem, intensity interferometry)
  • Nonlinear optics = breakdown of linear superposition when EEatomicE \sim E_{\text{atomic}}; generates harmonics, mixes frequencies, enables solitons and parametric amplification.
  • Universal toolkit:
    • Eikonal |ψ|2=n2|\nabla\psi|^2 = n^2 for rays
    • Fresnel-Kirchhoff integral for diffraction
    • Wiener-Khinchin / Van Cittert-Zernike for coherence
    • χ(n)\chi^{(n)} susceptibility expansion for nonlinearity
  • Practical payoff: explains imaging resolution, telescopes & microscopes, lasers, fiber optics, holography, LIGO, parametric oscillators, frequency combs.


Master Map

mindmap
  root((Optics))
    Geometric Optics (Ch.7)
      Eikonal equation
      Ray equation
      Fermat principle
      Snell's law
      Hamiltonian rays
      Polarization transport
      Paraxial / Gaussian beams
      Caustics
      WKB
    Diffraction (Ch.8)
      Helmholtz equation
      Huygens-Fresnel
      Fresnel-Kirchhoff
      Fraunhofer
      Fresnel zones
      Airy, Rayleigh
      Fourier optics
      Babinet
    Coherence (Ch.9)
      Temporal coherence
      Spatial coherence
      Mutual coherence Γ
      Van Cittert-Zernike
      Michelson stellar
      Fabry-Perot
      HBT intensity interferometry
      Photon statistics
      Holography
    Nonlinear (Ch.10)
      χ⁽ⁿ⁾ expansion
      SHG, sum/diff freq
      Phase matching
      Parametric amp
      Four-wave mixing
      Kerr / Self-focusing
      Solitons
      Raman, Brillouin


Geometric Optics & the Eikonal Approximation


From Wave Equation to Eikonal

Start from scalar wave equation (or component of Maxwell) in slowly-varying medium n(𝐱)n(\mathbf{x}):

2ψn2c22ψt2=0\nabla^2 \psi - \frac{n^2}{c^2}\frac{\partial^2 \psi}{\partial t^2} = 0

Eikonal ansatz (WKB):

ψ(𝐱,t)=A(𝐱)ei[k0S(𝐱)ωt],k0=ω/c\psi(\mathbf{x}, t) = A(\mathbf{x})\, e^{i[k_0 S(\mathbf{x}) - \omega t]}, \quad k_0 = \omega/c

Substitute, separate large-k0k_0 and finite terms:

Order Equation Name
O(k02)O(k_0^2) $ \nabla S
O(k0)O(k_0) 2(S)A+A2S=02(\nabla S)\cdot\nabla A + A \nabla^2 S = 0 Transport equation for amplitude

Validity: |n|/nk0n|\nabla n|/n \ll k_0 n, i.e., nn varies over scales λ\gg \lambda.


Rays as Trajectories

Definition: ray = curve everywhere parallel to S\nabla S (gradient of phase).

Tangent vector 𝐭̂=S/|S|=S/n\hat{\mathbf{t}} = \nabla S/|\nabla S| = \nabla S/n.

Ray equation (from d/dsd/ds of eikonal):

dds(n𝐭̂)=n\boxed{\frac{d}{ds}\!\left(n\,\hat{\mathbf{t}}\right) = \nabla n}

Equivalent vector form: dds(nd𝐱ds)=n\frac{d}{ds}\!\left(n\frac{d\mathbf{x}}{ds}\right) = \nabla n.

Special cases:

  • Homogeneous medium (n=0\nabla n = 0): straight lines.
  • Snell's law (interface, nn jump): tangential component of n𝐭̂n\hat{\mathbf{t}} continuous → n1sinθ1=n2sinθ2n_1\sin\theta_1 = n_2\sin\theta_2.
  • Gradient-index (GRIN): ray bends toward higher nn (toward denser/cooler air → mirages, atmospheric refraction).

Fermat's Principle

The eikonal equation is equivalent to:

δnds=0\delta \int n\, ds = 0

— rays extremize the optical path length between two points.

flowchart LR
    A["Fermat\nδ∫n ds = 0"] <--> B["Eikonal\n∣∇S∣² = n²"]
    A <--> C["Ray eq\n(d/ds)(n t̂) = ∇n"]
    B <--> D["Snell's law\nn₁ sinθ₁ = n₂ sinθ₂"]
    A <--> D

Analogous to Maupertuis principle in mechanics (δpdq=0\delta\int p\,dq = 0). Provides the bridge "geometric optics ↔︎ classical mechanics," with npn \leftrightarrow p.


Hamiltonian Formulation

Treat phase SS as Hamilton's action; rays as trajectories in phase space (𝐱,𝐤)(\mathbf{x}, \mathbf{k}) with 𝐤=k0S\mathbf{k} = k_0\nabla S playing role of momentum.

Optical Hamiltonian:

H(𝐱,𝐤)=12k0(|𝐤|2n2n2)k0(one form)orH=12(k2k02n2)H(\mathbf{x}, \mathbf{k}) = \frac{1}{2k_0}\!\left(\frac{|\mathbf{k}|^2}{n^2} - n^2\right) k_0 \;\;\text{(one form)}\quad \text{or}\quad H = \tfrac{1}{2}(k^2 - k_0^2 n^2)

Hamilton's equations:

d𝐱dσ=𝐤H,d𝐤dσ=𝐱H\frac{d\mathbf{x}}{d\sigma} = \partial_{\mathbf{k}} H, \qquad \frac{d\mathbf{k}}{d\sigma} = -\partial_{\mathbf{x}} H

with σ\sigma a parameter along the ray.

Liouville's theorem holds: phase-space volume in (𝐱,𝐤)(\mathbf{x},\mathbf{k}) space is conserved along rays. This is the basis of conservation of étendue (area × solid angle × n2n^2) in optical systems.


Polarization Transport (Parallel Transport along a Ray)

For Maxwell (not scalar) waves, polarization vector 𝐞̂\hat{\mathbf{e}} obeys:

d𝐞̂ds=(𝐞̂lnn)𝐭̂\frac{d\hat{\mathbf{e}}}{ds} = -(\hat{\mathbf{e}}\cdot\nabla\ln n)\,\hat{\mathbf{t}}

i.e., the projection onto the wavefront is parallel-transported.

Consequence: in a twisted ray (e.g., optical fiber with torsion), polarization picks up a geometric (Berry) phase = solid angle subtended by 𝐭̂\hat{\mathbf{t}} on the unit sphere. Verified in helical fibers (Tomita & Chiao 1986).


Caustics

A caustic is where neighboring rays converge → amplitude diverges → eikonal breaks down. Locally, two rays merge; the geometric-optics field has a square-root singularity.

Examples: focal point of a lens; bright curve on bottom of a swimming pool; rainbow caustic (Airy function structure across it).

Topology: Arnold's classification — fold, cusp, swallowtail, butterfly, umbilic. Locally diffraction smooths the singularity into Airy- and Pearcey-function patterns.


Paraxial Approximation & Gaussian Beams

Paraxial = rays make small angles with optical axis.

Paraxial wave equation: write ψ=u(𝐫,z)eikz\psi = u(\mathbf{r}_\perp, z)\,e^{ikz}, uu slowly varying:

2ikuz+2u=0\boxed{2ik\frac{\partial u}{\partial z} + \nabla_\perp^2 u = 0}

— a Schrödinger equation in (𝐫,z)(\mathbf{r}_\perp, z) with ztz\leftrightarrow t.

Fundamental Gaussian beam (TEM₀₀ mode):

u(r,z)=w0w(z)exp[r2w(z)2]exp[ikz+ikr22R(z)iη(z)]u(r, z) = \frac{w_0}{w(z)}\exp\!\left[-\frac{r^2}{w(z)^2}\right]\exp\!\left[ikz + ik\frac{r^2}{2R(z)} - i\eta(z)\right]

Parameter Formula Meaning
Beam waist w0w_0 Minimum radius
Rayleigh range zR=πw02/λz_R = \pi w_0^2/\lambda Distance over which beam doubles in area
Beam radius w(z)=w01+(z/zR)2w(z) = w_0\sqrt{1 + (z/z_R)^2} 1/e1/e radius
Radius of curvature R(z)=z[1+(zR/z)2]R(z) = z[1 + (z_R/z)^2] Of wavefront
Gouy phase η(z)=arctan(z/zR)\eta(z) = \arctan(z/z_R) π/2\pi/2 accumulated through focus
Far-field divergence θ=λ/(πw0)\theta = \lambda/(\pi w_0) Half-angle of cone

Key invariant (uncertainty-style):

w0×θ=λ/πsmall waist ⇔ large divergencew_0 \times \theta = \lambda/\pi \quad \Rightarrow \quad \text{small waist ⇔ large divergence}


ABCD Matrix Optics (Paraxial Ray Tracing)

Track a ray by (y,y)(y, y') where yy = height from axis, y=dy/dzy' = dy/dz:

(youtyout)=(ABCD)(yinyin)\begin{pmatrix} y_{\text{out}} \\ y'_{\text{out}}\end{pmatrix} = \begin{pmatrix} A & B \\ C & D \end{pmatrix}\begin{pmatrix} y_{\text{in}} \\ y'_{\text{in}}\end{pmatrix}

Element Matrix
Free space length dd (1d01)\begin{pmatrix} 1 & d \\ 0 & 1\end{pmatrix}
Thin lens focal ff (101/f1)\begin{pmatrix} 1 & 0 \\ -1/f & 1\end{pmatrix}
Flat interface n1n2n_1\to n_2 (100n1/n2)\begin{pmatrix} 1 & 0 \\ 0 & n_1/n_2\end{pmatrix}
Curved interface RR (10(n1n2)/(n2R)n1/n2)\begin{pmatrix} 1 & 0 \\ (n_1-n_2)/(n_2R) & n_1/n_2\end{pmatrix}
Mirror radius RR (102/R1)\begin{pmatrix} 1 & 0 \\ -2/R & 1\end{pmatrix}

System matrix = product of element matrices in reverse order of propagation. Always det=nin/nout\det = n_{\text{in}}/n_{\text{out}}.

For Gaussian beam propagation: complex parameter qq satisfies q=(Aq+B)/(Cq+D)q' = (A q + B)/(C q + D), where 1/q=1/Riλ/(πw2)1/q = 1/R - i\lambda/(\pi w^2).



Diffraction


Helmholtz Equation & Green's Function

For monochromatic field ψ(𝐱,t)=U(𝐱)eiωt\psi(\mathbf{x},t) = U(\mathbf{x}) e^{-i\omega t}:

(2+k2)U=0(Helmholtz)(\nabla^2 + k^2) U = 0 \quad\text{(Helmholtz)}

Spherical-wave Green's function:

G(r)=eikrrG(r) = \frac{e^{ikr}}{r}

— solution to (2+k2)G=4πδ3(𝐱)(\nabla^2 + k^2)G = -4\pi\delta^3(\mathbf{x}).


Huygens-Fresnel Principle (Modernized)

Every point on a wavefront acts as a source of a secondary spherical wave; the total field is the superposition.

Kirchhoff's integral theorem:

U(P)=14πV[UnGGnU]dSU(P) = \frac{1}{4\pi}\oint_{\partial V}\!\left[U\,\partial_n G - G\,\partial_n U\right]\, dS

For diffraction from an aperture Σ\Sigma (illuminated by spherical/plane wave), Fresnel-Kirchhoff formula:

U(P)=iλΣU0(P0)eikrrK(θ)dA\boxed{U(P) = \frac{-i}{\lambda}\int_\Sigma U_0(P_0)\,\frac{e^{ikr}}{r}\,K(\theta)\, dA}

with obliquity factor K(θ)=12(1+cosθ)K(\theta) = \tfrac{1}{2}(1 + \cos\theta) (Kirchhoff) or cosθ\cos\theta (Rayleigh-Sommerfeld).


Fresnel Number — Regime Classification

For aperture of size aa, viewing at distance zz, wavelength λ\lambda:

NF=a2λz\boxed{N_F = \frac{a^2}{\lambda z}}

Regime NFN_F Approximation
Near field (Fresnel) NF1N_F \gg 1 Quadratic-phase approximation; complicated patterns inside geometric shadow
Far field (Fraunhofer) NF1N_F \ll 1 Drop quadratic phase; pattern = Fourier transform of aperture
Geometric optics NFN_F \to \infty Sharp shadow

Fraunhofer Diffraction = Fourier Optics

For an aperture in plane z=0z=0 with field U0(𝐫)U_0(\mathbf{r}_\perp), observed at distance za2/λz \gg a^2/\lambda:

U(𝛒)=eikziλzeikρ2/2zU0(𝐫)ei𝐤𝐫d2r\boxed{U(\boldsymbol\rho) = \frac{e^{ikz}}{i\lambda z} e^{ik\rho^2/2z} \int U_0(\mathbf{r}_\perp)\, e^{-i\mathbf{k}_\perp\cdot\mathbf{r}_\perp}\, d^2 r_\perp}

with 𝐤=k𝛒/z\mathbf{k}_\perp = k\boldsymbol\rho/z. The far field is the 2-D spatial Fourier transform of the aperture amplitude.

Standard apertures:

Aperture Field U0U_0 Far-field intensity I(θ)I(\theta)
Slit width aa rect I0sinc2(kasinθ/2)I_0\,\text{sinc}^2(ka\sin\theta/2)
Two slits, sep. dd, width aa comb × rect sinc² × cos² (envelope × fringes)
NN-slit grating, period dd comb (N teeth) × rect I0sinc2[sin(Nkdsinθ/2)/sin(kdsinθ/2)]2I_0\,\text{sinc}^2 \cdot [\sin(Nkd\sin\theta/2)/\sin(kd\sin\theta/2)]^2
Circular aperture, radius aa disk I0[2J1(kasinθ)/(kasinθ)]2I_0\,[2J_1(ka\sin\theta)/(ka\sin\theta)]^2 (Airy)
Gaussian Gaussian Gaussian

Airy disk:

  • First zero at kasinθ=3.8317sinθ=1.22λ/(2a)=1.22λ/Dka\sin\theta = 3.8317 \Rightarrow \sin\theta = 1.22\lambda/(2a) = 1.22\lambda/D for diameter DD.
  • 84% of power within first zero.

Rayleigh criterion: two point sources resolved when one's center sits on other's first dark ring:

θmin=1.22λ/D(circular aperture)\boxed{\theta_{\min} = 1.22\,\lambda/D \quad \text{(circular aperture)}}


Diffraction Grating

NN slits, period dd, width aa:

I(θ)=I0sinc2(πasinθλ)single-slit envelope[sin(Nπdsinθ/λ)sin(πdsinθ/λ)]2principal maxima at dsinθ=mλI(\theta) = I_0\,\underbrace{\text{sinc}^2\!\left(\tfrac{\pi a \sin\theta}{\lambda}\right)}_{\text{single-slit envelope}}\cdot\underbrace{\left[\frac{\sin(N\pi d\sin\theta/\lambda)}{\sin(\pi d\sin\theta/\lambda)}\right]^2}_{\text{principal maxima at } d\sin\theta = m\lambda}

Quantity Formula Meaning
Principal maxima dsinθm=mλd\sin\theta_m = m\lambda mm = diffraction order
Width of max Δθλ/(Ndcosθm)\Delta\theta \approx \lambda/(Nd\cos\theta_m) Sharper with more slits
Resolving power R=λ/Δλ=mNR = \lambda/\Delta\lambda = mN Higher order + more slits = better
Free spectral range ΔλFSR=λ/m\Delta\lambda_{FSR} = \lambda/m Beyond this, orders overlap

Fresnel Diffraction & Zones

For near-field, must keep quadratic phase. The aperture is divided into Fresnel zones — annular regions whose contributions to U(P)U(P) alternate in sign.

Fresnel zone radii (axis-to-zone-edge distance, for source/observation at distance zz): rn=nλzr_n = \sqrt{n\lambda z}

Field at axis = approximate sum of alternating contributions:

U(axis)12(U1+Ulast)U(\text{axis}) \approx \tfrac{1}{2}(U_1 + U_{\text{last}})

Zone plate = alternating opaque/clear zones → focuses like a lens (without refraction). Focal length f=a2/(λ)f = a^2/(\lambda) for first-zone radius aa.

Fresnel integrals (used for slit edge, knife edge):

C(u)=0ucos(πt2/2)dt,S(u)=0usin(πt2/2)dtC(u) = \int_0^u \cos(\pi t^2/2)\, dt, \quad S(u) = \int_0^u \sin(\pi t^2/2)\, dt

→ Cornu spiral.


Babinet's Principle

Field of an aperture + field of its complement = field of unobstructed wave:

Uaperture(P)+Ucomplement(P)=Uunobstructed(P)U_{\text{aperture}}(P) + U_{\text{complement}}(P) = U_{\text{unobstructed}}(P)

Consequence: away from the geometric shadow (where Uunobstructed=0U_{\text{unobstructed}} = 0), an opaque disk and a circular hole produce the same diffraction pattern in intensity.


Diffraction from the Information-Theory View

Optical imaging = bandlimited spatial Fourier reconstruction.

Object plane Lens Image plane
f(𝐫)f(\mathbf{r}) spatial filter H(𝐤)H(\mathbf{k}) g=f*hg = f * h

Abbe / Rayleigh resolution limit: maximum spatial frequency the lens admits is kmax=(2π/λ)sinα=2πNA/λk_{\max} = (2\pi/\lambda)\sin\alpha = 2\pi\,\text{NA}/\lambda. The image-plane resolution is

Δxλ/(2NA)\Delta x \approx \lambda/(2\,\text{NA})

NA = numerical aperture = nsinαn\sin\alpha at object side.

Super-resolution techniques (STED, PALM, structured illumination, near-field) break this classical limit by exploiting nonlinearity or proximity.



Interferometry & Coherence


Two-Beam Interference (Coherent Case)

Two fields E1,E2E_1, E_2 with intensities I1,I2I_1, I_2 and phase difference δ\delta:

Itot=I1+I2+2I1I2cosδI_{\text{tot}} = I_1 + I_2 + 2\sqrt{I_1 I_2}\cos\delta

Visibility (fringe contrast):

V=ImaxIminImax+Imin=2I1I2I1+I2(coherent)\boxed{V = \frac{I_{\max} - I_{\min}}{I_{\max} + I_{\min}} = \frac{2\sqrt{I_1 I_2}}{I_1 + I_2} \quad \text{(coherent)}}

Equal intensities I1=I2I_1=I_2V=1V=1. Unequal or partially coherent ⇒ V<1V<1.


Multi-Beam Interference — Fabry-Perot

Two parallel mirrors, reflectivity RR, separation LL:

ItI0=11+Fsin2(δ/2),F=4R(1R)2,δ=2kLcosθ\frac{I_t}{I_0} = \frac{1}{1 + F\sin^2(\delta/2)}, \quad F = \frac{4R}{(1-R)^2}, \quad \delta = 2kL\cos\theta

— the Airy function (transmission, not Airy disk).

Quantity Formula Meaning
FSR ΔνFSR=c/(2L)\Delta\nu_{FSR} = c/(2L) Spacing between transmission peaks
Finesse =πF/2=πR/(1R)\mathcal{F} = \pi\sqrt{F}/2 = \pi\sqrt{R}/(1-R) FSR/FWHM
Resolution Δν=c/(2L)\Delta\nu = c/(2L\mathcal{F}) Linewidth
Q factor Q=ν/Δν=2L/λQ = \nu/\Delta\nu = 2L\mathcal{F}/\lambda Photon lifetime

Used in: laser cavities, optical filters, gravitational wave detectors (LIGO arm cavities, 450\mathcal{F}\sim 450), atomic clocks.


Coherence — Statistical Description of Light

Real light has finite spectral width and angular extent. Coherence theory quantifies "how well-defined is the phase" in time + space.

Mutual coherence function (two points 𝐫1,𝐫2\mathbf{r}_1, \mathbf{r}_2, time lag τ\tau):

Γ12(τ)=E*(𝐫1,t)E(𝐫2,t+τ)\Gamma_{12}(\tau) = \langle E^*(\mathbf{r}_1, t)\, E(\mathbf{r}_2, t+\tau)\rangle

Normalized:

γ12(τ)=Γ12(τ)Γ11(0)Γ22(0),|γ12|1\gamma_{12}(\tau) = \frac{\Gamma_{12}(\tau)}{\sqrt{\Gamma_{11}(0)\,\Gamma_{22}(0)}}, \quad |\gamma_{12}| \le 1

General visibility:

V=2I1I2I1+I2|γ12|V = \frac{2\sqrt{I_1 I_2}}{I_1 + I_2}\, |\gamma_{12}|


Temporal Coherence (Self-Coherence)

For single point (𝐫1=𝐫2\mathbf{r}_1 = \mathbf{r}_2), γ(τ)\gamma(\tau) measures how the field correlates with delayed copy of itself.

Coherence time: τc1/Δν\tau_c \sim 1/\Delta\nu (width of |γ(τ)||\gamma(\tau)|). Coherence length: Lc=cτcL_c = c\tau_c.

Wiener-Khinchin (again): Γ11(τ)\Gamma_{11}(\tau) ↔︎ power spectrum S(ω)S(\omega). So:

  • Spectrally narrow source ⇒ long coherence time/length.
  • White light: τc1/Δω\tau_c \sim 1/\Delta\omega \sim optical cycle.

Typical numbers:

Source Δν\Delta\nu LcL_c
White light (vis) 3×1014\sim 3\times 10^{14} Hz 1μ\sim 1\,\mum
LED 1013\sim 10^{13} Hz 30μ\sim 30\,\mum
HeNe laser (multimode) 109\sim 10^9 Hz 30\sim 30 cm
Stabilized HeNe (single mode) 105\sim 10^5 Hz 3\sim 3 km
Best optical clock laser <1< 1 Hz >108> 10^8 m

Michelson interferometer sweeps τ\tau; envelope of fringes = |γ11(τ)||\gamma_{11}(\tau)|Fourier-transform spectroscopy (FTIR).


Spatial Coherence

For lag τ=0\tau = 0 but two separated points 𝐫1,𝐫2\mathbf{r}_1, \mathbf{r}_2: γ12(0)\gamma_{12}(0) measures spatial coherence.

Van Cittert-Zernike theorem: for a distant, incoherent, quasi-monochromatic source with intensity distribution Is(𝛔)I_s(\boldsymbol\sigma) (angular):

γ12(0)=Is(𝛔)ei𝐤(𝐫2𝐫1)𝛔/Rd2σIsd2σ\boxed{\gamma_{12}(0) = \frac{\int I_s(\boldsymbol\sigma)\, e^{-i\mathbf{k}\cdot(\mathbf{r}_2-\mathbf{r}_1)\cdot\boldsymbol\sigma/R}\, d^2\sigma}{\int I_s\, d^2\sigma}}

— the spatial coherence function is the Fourier transform of the source brightness distribution.

Implications:

  • Coherence area on Earth ≈ λR/d\lambda R/d where dd = source size.
  • For Sun (d1.4×109d\approx 1.4\times 10^9 m at R1.5×1011R\approx 1.5\times 10^{11} m, λ=500\lambda = 500 nm): 50μ\sim 50\,\mum. Why pinhole experiments need pinholes.
  • For star (apparent diameter θ*\theta_*): coherence baseline = λ/θ*\lambda/\theta_*.

Stellar Interferometry

Michelson stellar interferometer (1920, measured Betelgeuse): two-aperture interference; baseline DDV=|J1(πθ*D/λ)|V = |J_1(\pi\theta_*D/\lambda)| for uniform disk. First null gives θ*=1.22λ/D\theta_* = 1.22\lambda/D.

Modern long-baseline interferometers (CHARA, VLTI, GRAVITY): kilometer baselines → milli-arcsec resolution. Used for stellar diameters, exoplanet imaging, AGN black holes (Event Horizon Telescope = VLBI at radio wavelengths).


Hanbury Brown–Twiss (Intensity Interferometry)

Correlate intensities at two detectors instead of amplitudes:

g(2)(τ)=I(t)I(t+τ)I2g^{(2)}(\tau) = \frac{\langle I(t)I(t+\tau)\rangle}{\langle I\rangle^2}

For thermal/chaotic light: g(2)(0)=2g^{(2)}(0) = 2 ("bunching"), decaying to 1 at ττc\tau\gg\tau_c. For coherent (laser): g(2)(τ)=1g^{(2)}(\tau) = 1 (uncorrelated arrivals). For single-photon source: g(2)(0)=0g^{(2)}(0) = 0 ("antibunching" — purely quantum, no classical wave can do this).

HBT-Siegert relation (for chaotic Gaussian light): g(2)(τ)=1+|g(1)(τ)|2g^{(2)}(\tau) = 1 + |g^{(1)}(\tau)|^2, with g(1)=γ11(τ)g^{(1)} = \gamma_{11}(\tau).

Advantage of intensity interferometry: insensitive to atmospheric phase fluctuations → very long baselines. HBT measured Sirius stellar diameter (1956). Modern revivals (VERITAS, etc.) probe stellar surfaces.


Photon Statistics

For monochromatic light, intensity ⟶ photon counts \to:

Light n\langle n\rangle (Δn)2\langle(\Delta n)^2\rangle Mandel QQ
Thermal (chaotic) n\bar n n+n2\bar n + \bar n^2 (super-Poisson) n>0\bar n > 0
Coherent (laser) n\bar n n\bar n (Poisson) 00
Number/Fock state nn 00 1-1
Squeezed light n\bar n <n< \bar n (sub-Poisson) <0<0

Mandel Q=(Δn2n)/nQ = (\langle\Delta n^2\rangle - \langle n\rangle)/\langle n\rangle.

Standard quantum limit (SQL) for phase: ΔϕSQL=1/n\Delta\phi_{SQL} = 1/\sqrt{\bar n} — Poisson counting noise. Squeezed light beats SQL (used in upgraded LIGO).


Holography

Record interference pattern between object wave UOU_O and reference wave URU_R on a photographic plate; intensity recorded:

I=|UO+UR|2=|UO|2+|UR|2+UR*UO+URUO*I = |U_O + U_R|^2 = |U_O|^2 + |U_R|^2 + U_R^* U_O + U_R U_O^*

Readout: illuminate hologram with reference URU_R → third term gives |UR|2UO|U_R|^2 U_O → reconstructs object wave (virtual 3-D image).

Key feature: stores both amplitude and phase of object wave. Originally Gabor (1948, Nobel 1971) for electron microscopy; revolutionized with lasers (Leith-Upatnieks 1962). Modern variants: rainbow holograms, digital holography, holographic optical elements.



Nonlinear Optics


Nonlinear Polarization

In linear media: 𝐏=ϵ0χ(1)𝐄\mathbf{P} = \epsilon_0\chi^{(1)}\mathbf{E}.

In nonlinear: expand in powers of 𝐄\mathbf{E}:

Pi=ϵ0[χij(1)Ej+χijk(2)EjEk+χijkl(3)EjEkEl+]\boxed{P_i = \epsilon_0\!\left[\chi^{(1)}_{ij} E_j + \chi^{(2)}_{ijk} E_j E_k + \chi^{(3)}_{ijkl} E_j E_k E_l + \cdots\right]}

Symmetries:

  • χ(2)0\chi^{(2)} \ne 0 only in non-centrosymmetric media (BBO, KDP, LiNbO₃, GaAs).
  • χ(3)\chi^{(3)} allowed in all media (incl. glass, air).

Magnitudes: χ(2)1012\chi^{(2)} \sim 10^{-12} m/V, χ(3)1022\chi^{(3)} \sim 10^{-22} m²/V². Nonlinear effects need intense fields → lasers, especially pulsed.


Three-Wave Mixing (from χ(2)\chi^{(2)})

Two input waves at ω1,ω2\omega_1, \omega_2 generate output components at ω1±ω2\omega_1 \pm \omega_2 (and harmonics 2ω1,2ω22\omega_1, 2\omega_2, plus DC). Energy and momentum conservation:

ω3=ω1+ω2(or differences)\omega_3 = \omega_1 + \omega_2 \quad\text{(or differences)} 𝐤3=𝐤1+𝐤2(phase matching)\mathbf{k}_3 = \mathbf{k}_1 + \mathbf{k}_2 \quad \text{(phase matching)}

Coupled-wave equations (slowly varying envelope approximation):

dA3dz=iω3n3cdeffA1A2eiΔkz,Δk=k1+k2k3\frac{dA_3}{dz} = \frac{i\omega_3}{n_3 c} d_{\text{eff}} A_1 A_2\, e^{i\Delta k z}, \quad \Delta k = k_1 + k_2 - k_3

Conversion efficiency sinc2(ΔkL/2)\propto \text{sinc}^2(\Delta k\, L/2)phase mismatch kills it.


Phase Matching

Phase matching = ensuring Δk=0\Delta k = 0 over the interaction length. Hard because of dispersion ($n(\omega) \ne $ const).

Methods:

  1. Birefringent phase matching: use ordinary/extraordinary rays in anisotropic crystal. Tune ne(θ)n_e(\theta) via crystal orientation.
  2. Temperature tuning: n(T)n(T) varies; lock to phase-matching condition.
  3. Quasi-phase matching (QPM): periodically poled crystal (e.g., PPLN) flips sign of χ(2)\chi^{(2)} every coherence length Lc=π/ΔkL_c = \pi/\Delta k. Equivalent to grating momentum kick.

Second-Harmonic Generation (SHG)

Special case of three-wave mixing with ω1=ω2=ω\omega_1 = \omega_2 = \omega, ω3=2ω\omega_3 = 2\omega.

Conversion efficiency (low-depletion, phase-matched):

η=P2ωPω=8π2deff2L2Iωn3ϵ0cλ2\eta = \frac{P_{2\omega}}{P_\omega} = \frac{8\pi^2 d_{\text{eff}}^2 L^2 I_\omega}{n^3\epsilon_0 c\lambda^2}

— quadratic in length, linear in input intensity. Doubles 1064 nm Nd:YAG to 532 nm green; the basis of green laser pointers.


Parametric Processes

Parametric amplification: pump at ωp\omega_p, signal at ωs\omega_s → idler at ωi=ωpωs\omega_i = \omega_p - \omega_s (with ωs+ωi=ωp\omega_s + \omega_i = \omega_p, 𝐤s+𝐤i=𝐤p\mathbf{k}_s + \mathbf{k}_i = \mathbf{k}_p). Signal grows exponentially when phase-matched:

As(z)=As(0)cosh(γz),γIpA_s(z) = A_s(0)\cosh(\gamma z), \quad \gamma \propto \sqrt{I_p}

Optical parametric oscillator (OPO): place in cavity → coherent tunable source. Widely used for mid-IR / UV generation.

Spontaneous parametric down-conversion (SPDC): quantum process where pump photon spontaneously splits into entangled signal-idler pair. Foundation of quantum optics (Bell tests, quantum key distribution, two-photon interferometry).


Four-Wave Mixing & Kerr Effect (from χ(3)\chi^{(3)})

Generic: ω4=ω1+ω2ω3\omega_4 = \omega_1 + \omega_2 - \omega_3.

Kerr effect (degenerate FWM, ω1=ω2=ω3=ω\omega_1 = \omega_2 = \omega_3 = \omega):

n(I)=n0+n2I,n2=3χ(3)4ϵ0cn02n(I) = n_0 + n_2 I, \quad n_2 = \frac{3\chi^{(3)}}{4\epsilon_0 c n_0^2}

— refractive index depends on intensity.

Consequences:

Phenomenon Mechanism
Self-focusing Beam intensity ⇒ lens-like n2In_2 I; collapses above critical power Pcrλ2/n0n2P_{cr} \sim \lambda^2/n_0 n_2
Self-phase modulation (SPM) Time-dependent I(t)I(t) ⇒ time-dependent phase ⇒ spectral broadening
Cross-phase modulation (XPM) One beam's intensity modulates another's phase
Optical Kerr shutter Intense pulse rotates polarization of weak probe
Solitons Balance of SPM (spectral broadening) and GVD (group-velocity dispersion) in fiber

Solitons

Nonlinear Schrödinger equation (envelope AA in fiber):

iAzβ222At2+γ|A|2A=0i\frac{\partial A}{\partial z} - \frac{\beta_2}{2}\frac{\partial^2 A}{\partial t^2} + \gamma |A|^2 A = 0

with β2<0\beta_2 < 0 (anomalous dispersion) and γ>0\gamma > 0 (Kerr).

Fundamental soliton: A(t,z)=A0sech(t/T0)eiγA02z/2A(t, z) = A_0\,\text{sech}(t/T_0)\,e^{i\gamma A_0^2 z/2}. Propagates without distortion — dispersion and nonlinearity cancel.

Used in fiber-optic communications, mode-locked lasers (passive soliton mode-locking).


Stimulated Raman & Brillouin Scattering

SRS: photon scatters off molecular vibration; emits Stokes photon at ωΩv\omega - \Omega_v (with Ωv\Omega_v \sim THz for molecules). Gain Ip\propto I_p.

SBS: photon scatters off acoustic phonon; backward-propagating Stokes at ωΩB\omega - \Omega_B with ΩB\Omega_B \sim GHz. Limits power in long fibers (couples forward to backward).


Frequency Combs (a Modern Triumph)

Mode-locked laser → train of pulses with rep rate frf_r → spectrum is a comb of teeth at fn=nfr+f0f_n = n f_r + f_0 (f0f_0 = carrier-envelope offset).

Self-referencing (f-2f): broaden spectrum to octave; SHG of low-frequency tooth interferes with high-frequency tooth → measure f0f_0. Locking both frf_r and f0f_0 to atomic clocks ⇒ all comb teeth become frequency standards.

Impact: revolutionized optical frequency metrology, optical clocks, exoplanet RV spectroscopy, attosecond pulse generation. Hänsch & Hall Nobel 2005.



Workflow / Process

flowchart TD
    A[Optical problem] --> B{Length scales}
    B -->|λ ≪ L| C[Geometric optics]
    B -->|λ ~ L| D[Wave optics]
    B -->|λ ≫ L| E[Lumped/quasi-static]
    C --> F[Eikonal / ray equation]
    C --> G[Paraxial / Gaussian beam]
    D --> H{Aperture vs distance}
    H -->|N_F ≫ 1| I[Fresnel diffraction]
    H -->|N_F ≪ 1| J[Fraunhofer = Fourier transform]
    D --> K{Coherent or partially coherent?}
    K -->|Coherent| L[Add amplitudes, then intensity]
    K -->|Partial| M[Mutual coherence Γ]
    M --> N[Van Cittert-Zernike]
    A --> O{Intensity high?}
    O -->|Yes, &#124;E&#124; ~ &#124;E_atomic&#124;| P[Nonlinear χⁿ]
    P --> Q[Phase matching condition]


Comparison Tables


Wave/Ray Regimes

Regime Validity Key equation What's lost
Geometric optics λL\lambda \ll L Eikonal $ \nabla S
Paraxial / Gaussian beams θ1\theta \ll 1 2ikzu+2u=02ik\partial_z u + \nabla_\perp^2 u = 0 Wide-angle accuracy
Fresnel diffraction NF1N_F \gg 1 Quadratic-phase integral Edge waves (already in)
Fraunhofer diffraction NF1N_F \ll 1 Fourier transform Near-field detail
Full Maxwell always (2c2t2)E=0(\nabla^2 - c^{-2}\partial_t^2)E = 0 + sources Nothing

Imaging Resolution

System Resolution limit Comment
Eye (D≈3 mm) θ1\theta \sim 1' Atmospheric and biological limit ~ this
Telescope (D = 2.4 m, Hubble) θ=0.05\theta = 0.05'' At 550 nm
Microscope (NA = 1.4 oil) Δx=0.2μ\Delta x = 0.2\,\mum At 550 nm
Stellar interferometer (D=200 m baseline) θ=0.5\theta = 0.5 mas At 500 nm
EHT VLBI (Earth-size, mm wave) θ=20μ\theta = 20\,\muas M87, Sgr A* black holes
STED, PALM (super-resolution) <50< 50 nm Beats diffraction via nonlinearity / statistics

Coherence Concepts

Type Quantifies Function Limit
Temporal Spectral purity γ11(τ)\gamma_{11}(\tau) τc1/Δν\tau_c \sim 1/\Delta\nu
Spatial Angular purity γ12(0)\gamma_{12}(0) Area (λR/dsrc)2\sim (\lambda R/d_{\text{src}})^2
Second-order Photon statistics g(2)(τ)g^{(2)}(\tau) Classical 1\ge 1; quantum <1<1 possible

Nonlinear Processes (χ(2)\chi^{(2)})

Process Inputs Output Phase match
SHG ω,ω\omega, \omega 2ω2\omega 2k1=k22k_1 = k_2
Sum freq. ω1,ω2\omega_1, \omega_2 ω1+ω2\omega_1 + \omega_2 𝐤1+𝐤2=𝐤3\mathbf{k}_1+\mathbf{k}_2 = \mathbf{k}_3
Diff. freq. ω1,ω2\omega_1, \omega_2 ω1ω2\omega_1 - \omega_2 similar
Parametric amp ωp,ωs\omega_p, \omega_s ωi=ωpωs\omega_i = \omega_p-\omega_s 𝐤p=𝐤s+𝐤i\mathbf{k}_p = \mathbf{k}_s + \mathbf{k}_i
SPDC ωp\omega_p ωs,ωi\omega_s, \omega_i entangled similar; quantum

Photon Statistics Summary

Light g(2)(0)g^{(2)}(0) (Δn)2/n\langle(\Delta n)^2\rangle/\langle n\rangle Where
Thermal/chaotic 2 1+n1+\bar n Stars, blackbody
Coherent (laser) 1 1 Above threshold laser
Single-photon source 0 < 1 NV centers, quantum dots
Squeezed <1<1 or >1>1 variable Parametric amplifiers, GW detectors


Common Mistakes

  • Using geometric optics where λL\lambda \sim L. Always check NFN_F or eikonal validity.
  • Confusing Fresnel and Fraunhofer regimes. Plug in numbers for NFN_F before choosing approximation.
  • Adding intensities instead of amplitudes for coherent superposition. Coherent: |U1+U2|2|U_1 + U_2|^2. Incoherent: |U1|2+|U2|2|U_1|^2 + |U_2|^2.
  • Forgetting the obliquity factor / Kirchhoff sign. Different texts use different forms; verify before plugging into a problem.
  • Using paraxial Gaussian beam formulas at large angles. Breaks down when w0λw_0 \sim \lambda.
  • Treating Airy disk first zero as the resolution by Sparrow criterion. Sparrow uses where pattern flattens, 0.95λ/D\sim 0.95 \lambda/D. Rayleigh = 1.22λ/D1.22\lambda/D. Pick one consistently.
  • Confusing visibility VV with mutual coherence |γ||\gamma|. They coincide only for equal intensities.
  • Forgetting that VCZ gives Fourier transform of source intensity, not amplitude. Source is incoherent; coherence is built up by propagation distance.
  • Using g(2)(0)>1g^{(2)}(0) > 1 as proof of quantum effect. That's classical bunching (thermal light). Only g(2)(0)<1g^{(2)}(0) < 1 is non-classical.
  • Ignoring polarization in χ(2)\chi^{(2)} processes. It's a 3-tensor; ordinary/extraordinary mixing essential for phase matching.
  • Computing SHG efficiency without phase matching condition. A factor of 10410^410610^6 comes from sin2(ΔkL/2)/(ΔkL/2)2\sin^2(\Delta k L/2)/(\Delta k L/2)^2.
  • Confusing Kerr nonlinearity n2n_2 with electro-optic Pockels effect. Pockels is χ(2)\chi^{(2)} (linear in EE); Kerr is χ(3)\chi^{(3)} (quadratic in EE).
  • Treating coherence as binary (coherent/incoherent). It's continuous: 0|γ|10 \le |\gamma| \le 1.
  • Assuming all light is monochromatic. Real light has finite Δν\Delta\nu; affects coherence, interferometer fringe envelope.


Expert Insights

Geometric optics is the classical limit of wave optics, exactly analogous to classical mechanics being the limit of QM. Eikonal ↔︎ Hamilton-Jacobi; ray equation ↔︎ Newton's equation; λ\hbar \leftrightarrow \lambda.

Caustics are where geometric optics commits suicide. Diffraction integrals there yield Airy / Pearcey / higher catastrophe functions. The colors of rainbows live exactly in this regime.

The lens does a Fourier transform. This is the central insight of Fourier optics — image plane = FT of object plane. Spatial filtering in the Fourier plane = phase contrast, dark field, edge enhancement, modern phase microscopy.

Resolution ≠ sampling. A telescope's 1.22λ/D1.22\lambda/D is the angular resolution; a CCD pixel array can be coarser (under-sampling, aliasing) or finer (over-sampling). Match them at Nyquist.

Coherence is geometric. It's not a property of "the light" alone but of "the light + the experimental setup." A pinhole increases spatial coherence (at cost of flux). A monochromator increases temporal coherence (same trade-off).

The fluctuation-dissipation theorem hides in HBT. The bunching of thermal photons is a stat-mech statement about Δn\Delta n for Bose distribution — directly inherited from Part II.

Holography is a metaphor for many things. Storing phase = storing 3-D info. Quantum holography, gravitational holography (AdS/CFT) borrow the name.

The Van Cittert-Zernike theorem and the Wiener-Khinchin theorem are sibling FTs. One in space (incoherent source), one in time (random process). Both unmask randomness as Fourier dual of correlations.

Quasi-phase matching is the right answer. Most useful crystals can't be birefringently phase-matched at desired wavelengths. PPLN (periodically poled LiNbO₃) made nonlinear optics industrial.

Squeezed light is not "more photons in less space." It's a redistribution of the quantum uncertainty between amplitude and phase quadratures. Anti-squeezing the unmeasured quadrature is the cost.

Solitons are particles to fluid dynamicists, but pulses to optical engineers. Their stability comes from balance of two opposing effects (dispersion vs nonlinearity in optics; depth-dispersion vs amplitude in shallow water).

Frequency combs make light into a ruler. Phase-locking a comb to a microwave standard transfers atomic-clock precision to optical frequencies — billion-fold gain in precision.

Atmospheric phase fluctuations kill amplitude interferometry but not intensity interferometry. That's the deep reason HBT survives at intercontinental baselines while Michelson interferometers struggle past 100 m.

Numerical aperture is the same as sinθ\sin\theta for the marginal ray. Higher NA = wider acceptance angle = better resolution AND smaller depth of focus. Trade-off baked in.

Diffraction limit is not fundamental. It's the linear far-field limit. Nonlinear / near-field / statistical tricks (STED, PALM, NSOM) all break it routinely now.


Engineering use and verification

State the model, coordinate system, assumptions, boundary conditions and validity range before using an equation. Track dimensions and sign conventions through each derivation, test limiting cases, and distinguish mathematical possibility from physical realisability. Where a model informs engineering design, compare it with measurement or a second method and quantify the effect of idealisations rather than hiding them inside numerical precision.

  • Confirm scope, assumptions, interfaces and required outcome.
  • Check dimensions, sign conventions, boundary conditions and limiting cases.
  • Identify current project, customer and regulatory requirements.
  • Separate source examples from mandatory acceptance criteria.
  • Check calculations, tables and selections by an independent method.
  • Verify safety, maintainability and credible failure modes.
  • Record evidence, revisions, approvals and unresolved limitations.
  • Validate the result under representative operating conditions.

Continue learning

Statistical Physics and Kinetic Theory: Expert InsightsGuide · PhysicsNEXT LESSON →Geometric and Wave Optics: Final TakeawaysGuide · PhysicsStatistical Physics and Kinetic Theory: /20 SummaryGuide · PhysicsElasticity, Stress and Continuum Mechanics: /20 SummaryGuide · Physics