Coherence Functions and the Degree of Coherence

Coherence theory asks how much a light field resembles itself — one point compared to another, one instant compared to a later one. To make that precise it borrows three things from earlier in the chapter: from the quantization of the electromagnetic field, the mode operators a^\hat a and a^†\hat a^{\dagger}, so that "how bright" and "how correlated" become expectation values of normally ordered products; from photon statistics, the number distributions — Poissonian, Bose–Einstein, sub-Poissonian — that fix how photons cluster in time; and from the coherent states, the eigenstates of a^\hat a that behave most like a classical wave. Two notions of resemblance emerge, not one. First-order coherence compares the field amplitudes E(t)E(t) and E(t+τ)E(t+\tau); it governs interference and is carried by the complex degree of coherence g(1)(τ)g^{(1)}(\tau), whose modulus is literally the visibility of two-beam fringes. Second-order coherence compares the intensities I(t)∝∣E(t)∣2I(t)\propto|E(t)|^2; it governs photon-counting and is carried by g(2)(τ)g^{(2)}(\tau), the normalized rate of coincident clicks on two detectors. The decisive discovery — forced on physics by the 1956 Hanbury Brown–Twiss experiment — is that the two orders are logically independent. A thermal beam filtered to a razor-thin linewidth is as first-order coherent as a laser, yet it stays bunched at g(2)(0)=2g^{(2)}(0)=2 where the laser sits at 11; no amplitude measurement tells them apart, only intensity correlations. And when g(2)(0)g^{(2)}(0) drops below 11 — antibunching — no classical field can reproduce it, so the coherence hierarchy doubles as an operational witness of the quantized nature of light.

Key takeaways
  • First-order coherence g(1)(τ)=⟨E∗(t)E(t+τ)⟩/⟨∣E∣2⟩g^{(1)}(\tau)=\langle E^{*}(t)E(t+\tau)\rangle/\langle|E|^2\rangle measures amplitude correlation; its modulus equals fringe visibility, and via Wiener–Khinchin it is the Fourier partner of the spectrum — encoding the lineshape, nothing about photon statistics.
  • Second-order coherence g(2)(τ)=⟨I(t)I(t+τ)⟩/⟨I⟩2g^{(2)}(\tau)=\langle I(t)I(t+\tau)\rangle/\langle I\rangle^2 measures intensity (photon-number) correlation: g(2)(0)=1g^{(2)}(0)=1 for a laser, 22 for thermal light, 00 for a single photon.
  • For chaotic light the Siegert relation g(2)(τ)=1+∣g(1)(τ)∣2g^{(2)}(\tau)=1+|g^{(1)}(\tau)|^2 ties the two orders together, but only for that special Gaussian statistics — in general they are independent.
  • Classical fields must obey g(2)(0)≥1g^{(2)}(0)\ge1 and g(2)(0)≥g(2)(τ)g^{(2)}(0)\ge g^{(2)}(\tau); any violation is an unforgeable signature of quantized light.
  • Key idea: coherence is an infinite hierarchy of orders. First order g(1)g^{(1)} is amplitude correlation and equals visibility; second order g(2)g^{(2)} is intensity correlation and equals the normalized coincidence rate. The two are logically independent, and g(2)(0)<1g^{(2)}(0)<1 has no classical explanation.
The core insight

Coherence is an infinite hierarchy of orders. First order g(1)(τ)g^{(1)}(\tau) measures amplitude correlation and equals fringe visibility; second order g(2)(τ)g^{(2)}(\tau) measures intensity correlation and equals the normalized photon-coincidence rate. The two are logically independent, and g(2)(0)<1g^{(2)}(0)<1 has no classical explanation.

Where this lives today

The Hanbury Brown–Twiss layout of that 1956 experiment — a beam splitter, two detectors, a coincidence counter — is now the standard bench for measuring g(2)(τ)g^{(2)}(\tau), and it certifies every serious single-photon source. Quantum dots, nitrogen-vacancy centres in diamond, and defects in two-dimensional materials are graded by their measured g(2)(0)g^{(2)}(0), read out with superconducting-nanowire detectors whose picosecond timing resolves the antibunching dip directly. These sources feed quantum key distribution, where a low g(2)(0)g^{(2)}(0) closes the multiphoton loophole an eavesdropper exploits, and photonic quantum computers, where indistinguishable single photons are the raw material. The idea has travelled far beyond photons. Intensity interferometry, invented to size stars without phase stability, has been revived on arrays of Cherenkov telescopes to resolve stars at microarcsecond scales. And ultracold-atom experiments measure the atomic analogue of g(2)g^{(2)} and see the contrast quantum statistics demands — bunching for bosons, antibunching for a degenerate Fermi gas, the Pauli principle written into a coincidence histogram. In every case the same correlation function, born to describe fringes and clumped starlight, is doing the work.

The physical picture

Start with the older, first-order question. Split light in a Michelson interferometer and recombine it with a delay τ\tau. If the field keeps a steady phase across that delay, the beams add and subtract in step and bright-dark-bright fringes appear; if the phase wanders faster than τ\tau, the crest of one arm lands on the trough of the other as often as not and the fringes wash out to flat grey. The complex degree of coherence g(1)(τ)g^{(1)}(\tau) is the bookkeeping of that phase memory, and its magnitude is the fringe visibility. The delay at which the memory fades is the coherence time τc\tau_c, and cτcc\tau_c is the coherence length. Now the newer, second-order question, which no interferometer answers. Put two photon detectors on the outputs of a beam splitter and ask: given a click on A at time tt, is a click on B at t+τt+\tau more, less, or equally likely than chance? That normalized conditional rate is g(2)(τ)g^{(2)}(\tau). Thermal light says "more likely" — its photons arrive in clumps, so a first click flags a bright moment and a second follows: bunching, g(2)(0)=2g^{(2)}(0)=2. Laser light says "equally likely," g(2)=1g^{(2)}=1. A single-photon emitter says "impossible at zero delay" — one photon cannot fire both detectors at once, so g(2)(0)=0g^{(2)}(0)=0, antibunching.

Figure
Two orders of coherence, two instruments — an interferometer reads amplitude correlation as fringe visibility equal to the modulus of g-one, while a beam splitter with two detectors reads intensity correlation as the normalized coincidence rate g-two
Two questions, two instruments — the axes of coherence are independent First order — amplitude correlation interferometer · fringe visibility delay τ V = (I max − I min)/(I max + I min) = |g¹(τ)| fringe contrast reads the phase memory memory lasts a time τ_c; c·τ_c = coherence length Second order — intensity correlation beam splitter · coincidence rate D_A D_B × correlate g²(τ) = ⟨I(t)·I(t+τ)⟩ / ⟨I⟩² bunched 2 · random 1 · antibunched 0 no interferometer can measure this number

The crux is that these two axes do not talk to each other. First-order coherence is fixed by the spectrum — through Wiener–Khinchin, g(1)(τ)g^{(1)}(\tau) is the Fourier transform of the lineshape — while second-order coherence is fixed by the photon statistics. Filter a thermal lamp until its spectrum is as narrow as a laser's, matching g(1)g^{(1)} exactly, and you will not have removed one iota of its bunching. Quantum-mechanically, "perfectly coherent" means every correlation function factorizes, ⟨a^†⋯a^⟩=∏⟨a^†⟩⟨a^⟩\langle\hat a^{\dagger}\cdots\hat a\rangle=\prod\langle\hat a^{\dagger}\rangle\langle\hat a\rangle, which only coherent states achieve; every departure is a failure of that factorization, and below 11 no classical field can imitate it. The next figure plots the three canonical answers on one axis.

Figure
The g-two landscape versus delay for the three canonical sources: thermal light bunches at 2.0 and decays to 1 over the coherence time, laser light sits flat at 1.0, and a single photon dips to 0.0 at zero delay, with the Siegert value at one coherence time equal to 1.135
The g²(τ) landscape — one plot classifies every light source 2.0 1.0 0.0 g²(τ) delay τ (units of coherence time τ_c) 0 3 τ_c thermal · g²(0) = 2.0 (bunched) g²(τ_c) = 1 + e⁻² = 1.135 laser · g²(τ) = 1.0 (Poissonian) single photon · g²(0) = 0.0 (antibunched) classical light lives at or above 1; the blue dip below 1 is forbidden classically

The core equations

First-order coherence is built from the field autocorrelation. Under stationarity — statistics unchanged by shifting the time origin — every correlation depends only on the delay τ\tau, so we form the mutual coherence function and normalize by the mean intensity, which fixes g(1)(0)=1g^{(1)}(0)=1 and makes the modulus a number between 00 and 11.

Γ(1)(τ)=⟨E∗(t) E(t+τ)⟩,g(1)(τ)=⟨E∗(t) E(t+τ)⟩⟨∣E(t)∣2⟩(1)\Gamma^{(1)}(\tau)=\langle E^{*}(t)\,E(t+\tau)\rangle,\qquad g^{(1)}(\tau)=\frac{\langle E^{*}(t)\,E(t+\tau)\rangle}{\langle |E(t)|^2\rangle} \qquad (1)

where: E(t)E(t) — complex analytic signal of the scalar field at a fixed point; ⟨⋅⟩\langle\cdot\rangle — time (ensemble) average, stationarity assumed; Γ(1)\Gamma^{(1)} — mutual coherence function; g(1)(τ)g^{(1)}(\tau) — complex degree of first-order coherence, g(1)(0)=1g^{(1)}(0)=1. Physically it records how well the field's phase at t+τt+\tau is predicted from its phase at tt; a slow decay means a long coherence time. Second-order coherence compares the intensity now with the intensity a delay τ\tau later. It is a fourth-order field moment, normalized so an uncorrelated stream gives exactly 11 — the quantity a two-detector experiment measures.

g(2)(τ)=⟨I(t) I(t+τ)⟩⟨I(t)⟩2=⟨E∗(t)E∗(t+τ)E(t+τ)E(t)⟩⟨∣E(t)∣2⟩2(2)g^{(2)}(\tau)=\frac{\langle I(t)\,I(t+\tau)\rangle}{\langle I(t)\rangle^{2}}=\frac{\langle E^{*}(t)E^{*}(t+\tau)E(t+\tau)E(t)\rangle}{\langle |E(t)|^{2}\rangle^{2}} \qquad (2)

where: I(t)∝∣E(t)∣2I(t)\propto|E(t)|^2 — instantaneous intensity; g(2)(τ)g^{(2)}(\tau) — normalized intensity correlation, dimensionless. Above 11, bright moments coincide (bunching); below 11 they avoid one another (antibunching); at 11, intensity now says nothing about later. The link between the abstract g(1)g^{(1)} and a laboratory reading is fringe visibility: when a field interferes with a delayed copy of itself, the depth of the fringes is governed entirely by the modulus of the degree of coherence.

V=Imax⁡−Imin⁡Imax⁡+Imin⁡=∣g(1)(τ)∣(3)\mathcal{V}=\frac{I_{\max}-I_{\min}}{I_{\max}+I_{\min}}=\bigl|g^{(1)}(\tau)\bigr| \qquad (3)

where: V\mathcal{V} — Michelson fringe visibility at path delay τ\tau; Imax⁡,Imin⁡I_{\max},I_{\min} — maximum and minimum fringe intensities. Because 0≤V≤10\le\mathcal V\le1, this is why ∣g(1)∣|g^{(1)}| never exceeds 11: full fringes (V=1\mathcal V=1) mean perfect first-order coherence, a flat field (V=0\mathcal V=0) means none. First-order coherence carries the spectrum and nothing else. The Wiener–Khinchin theorem states that g(1)(τ)g^{(1)}(\tau) and the power spectral density are a Fourier-transform pair, so measuring coherence versus delay reconstructs the lineshape.

S(ω)∝∫−∞∞g(1)(τ) eiωτ dτ(4)S(\omega)\propto\int_{-\infty}^{\infty} g^{(1)}(\tau)\,e^{i\omega\tau}\,\mathrm{d}\tau \qquad (4)

where: S(ω)S(\omega) — power spectral density; g(1)(τ)g^{(1)}(\tau) — first-order coherence, the Fourier conjugate of the spectrum. A broad spectrum gives a fast-decaying g(1)g^{(1)}, a narrow one a slow decay: the coherence time is roughly the inverse bandwidth, τc∼1/Δν\tau_c\sim1/\Delta\nu. For chaotic (thermal) light alone, the two orders are tied by the Siegert relation. It is not a general law — it holds precisely because thermal light has complex-Gaussian statistics — but where it holds it is exact and it drives intensity interferometry.

g(2)(τ)=1+∣g(1)(τ)∣2(5)g^{(2)}(\tau)=1+\bigl|g^{(1)}(\tau)\bigr|^{2} \qquad (5)

where: g(2)(τ)g^{(2)}(\tau) — intensity correlation of chaotic light. At zero delay ∣g(1)(0)∣=1|g^{(1)}(0)|=1 so g(2)(0)=2g^{(2)}(0)=2; at large delay g(1)→0g^{(1)}\to0 so g(2)→1g^{(2)}\to1. The bunching decays from 22 to 11 over exactly the coherence time set by the spectrum. To describe nonclassical light the correlation functions are promoted to quantum expectation values of the mode operators, written in normal order — every creation operator to the left of every annihilation operator. This is not a matter of taste; it is forced by the physics of detection, which absorbs photons.

g(n)({ti})=⟨a^†(t1)⋯a^†(tn) a^(tn)⋯a^(t1)⟩∏i⟨a^†(ti)a^(ti)⟩(6)g^{(n)}(\{t_i\})=\frac{\bigl\langle \hat a^{\dagger}(t_1)\cdots\hat a^{\dagger}(t_n)\,\hat a(t_n)\cdots\hat a(t_1)\bigr\rangle}{\prod_i \langle \hat a^{\dagger}(t_i)\hat a(t_i)\rangle} \qquad (6)

where: a^,a^†\hat a,\hat a^{\dagger} — annihilation/creation operators of the detected mode; the numerator is normally ordered. Because detection clicks by absorbing a photon, the coincidence rate is ⟨a^†a^†a^a^⟩\langle\hat a^{\dagger}\hat a^{\dagger}\hat a\hat a\rangle; normal ordering discards the vacuum term, so an empty mode gives zero counts. For a single mode at zero delay the numerator collapses to a photon-number expression, so g(2)(0)g^{(2)}(0) becomes a direct reporter of the number variance — the equation that connects coherence to photon statistics and exposes the classical boundary.

g(2)(0)=⟨a^†a^†a^a^⟩⟨a^†a^⟩2=1+⟨(Δn^)2⟩−⟨n^⟩⟨n^⟩2(7)g^{(2)}(0)=\frac{\langle \hat a^{\dagger}\hat a^{\dagger}\hat a\hat a\rangle}{\langle \hat a^{\dagger}\hat a\rangle^{2}}=1+\frac{\langle(\Delta \hat n)^{2}\rangle-\langle\hat n\rangle}{\langle\hat n\rangle^{2}} \qquad (7)

where: n^=a^†a^\hat n=\hat a^{\dagger}\hat a — number operator; ⟨(Δn^)2⟩\langle(\Delta\hat n)^2\rangle — number variance. Poissonian (⟨(Δn^)2⟩=⟨n^⟩\langle(\Delta\hat n)^2\rangle=\langle\hat n\rangle) gives g(2)(0)=1g^{(2)}(0)=1; super-Poissonian gives bunching above 11; sub-Poissonian (⟨(Δn^)2⟩<⟨n^⟩\langle(\Delta\hat n)^2\rangle<\langle\hat n\rangle) drives it below 11 — the exclusively quantum regime. Any classical field obeys two inequalities, and their violation is the operational definition of nonclassical light. Both follow from the fact that a classical intensity is real and non-negative.

g(2)(0)≥1,g(2)(0)≥g(2)(τ)(8)g^{(2)}(0)\ge1,\qquad g^{(2)}(0)\ge g^{(2)}(\tau) \qquad (8)

where: g(2)(0)g^{(2)}(0) — zero-delay intensity correlation. Classical correlations are strongest at zero delay, so a measured g(2)(0)<1g^{(2)}(0)<1, or a correlation that grows away from zero delay, cannot come from any classical field.

Derivation AFringe visibility equals the degree of coherence

The goal is to prove equation (3), that the modulus of the first-order degree of coherence is exactly the visibility of two-beam fringes — the bridge that turns an abstract correlation function into a number you read off an interferogram. We superpose the field with an equal-amplitude delayed copy, average the output intensity, identify the interference term with Γ(1)(τ)\Gamma^{(1)}(\tau), and extract the fringe extrema, assuming only stationarity and equal beam splitting.

Step 1. Write the output as the coherent sum of the two arms and time-average, I=⟨∣E(t)+E(t+τ)∣2⟩I=\langle|E(t)+E(t+\tau)|^{2}\rangle. Why: the device adds the field to a time-delayed copy, and the detector responds to the squared total.

Step 2. Expand into ⟨∣E(t)∣2⟩+⟨∣E(t+τ)∣2⟩+2 Re⟨E∗(t)E(t+τ)⟩\langle|E(t)|^2\rangle+\langle|E(t+\tau)|^2\rangle+2\,\mathrm{Re}\langle E^{*}(t)E(t+\tau)\rangle; by stationarity the first two are each ⟨∣E∣2⟩\langle|E|^2\rangle and the last is 2 Re Γ(1)(τ)2\,\mathrm{Re}\,\Gamma^{(1)}(\tau), so I=2⟨∣E∣2⟩(1+Re g(1)(τ))I=2\langle|E|^2\rangle\bigl(1+\mathrm{Re}\,g^{(1)}(\tau)\bigr). Why: the two arms carry equal mean intensity, and the cross term is the coherence function that normalizes into g(1)g^{(1)}.

Step 3. Write g(1)(τ)=∣g(1)∣ eiϕ(τ)g^{(1)}(\tau)=|g^{(1)}|\,e^{i\phi(\tau)}, with ϕ\phi sweeping through ω0τ\omega_0\tau as the delay scans, giving I∝1+∣g(1)∣cos⁡ϕI\propto1+|g^{(1)}|\cos\phi. Why: scanning the delay traces out a cosine fringe of amplitude ∣g(1)∣|g^{(1)}| about unit offset.

Step 4. Read the extrema at cos⁡ϕ=±1\cos\phi=\pm1: Imax⁡∝1+∣g(1)∣I_{\max}\propto1+|g^{(1)}|, Imin⁡∝1−∣g(1)∣I_{\min}\propto1-|g^{(1)}|, so V=(Imax⁡−Imin⁡)/(Imax⁡+Imin⁡)=∣g(1)(τ)∣\mathcal V=(I_{\max}-I_{\min})/(I_{\max}+I_{\min})=|g^{(1)}(\tau)|. Why: visibility strips the constant background and reports only the fringe amplitude — the coherence modulus itself.

Go deeper · Derivation B · 1 minThe Siegert relation for chaotic light
Derivation BThe Siegert relation for chaotic light

The aim is to derive equation (5), g(2)(τ)=1+∣g(1)(τ)∣2g^{(2)}(\tau)=1+|g^{(1)}(\tau)|^{2}, for thermal light and thereby explain why chaotic light always bunches to g(2)(0)=2g^{(2)}(0)=2. The engine is the Gaussian moment theorem: a chaotic field is a sum of many independent emitters, so by the central limit theorem its amplitude is a circular Gaussian variable whose higher moments factorize into second moments, which we apply to the fourth-order moment in g(2)g^{(2)}.

Step 1. For a complex Gaussian field the fourth-order moment splits as ⟨E1∗E2∗E2E1⟩=⟨E1∗E1⟩⟨E2∗E2⟩+⟨E1∗E2⟩⟨E2∗E1⟩\langle E_1^{*}E_2^{*}E_2E_1\rangle=\langle E_1^{*}E_1\rangle\langle E_2^{*}E_2\rangle+\langle E_1^{*}E_2\rangle\langle E_2^{*}E_1\rangle, with E1=E(t)E_1=E(t), E2=E(t+τ)E_2=E(t+\tau). Why: the Gaussian (Isserlis–Wick) theorem replaces a fourth-order average by all pairings of second-order averages, exact for thermal statistics.

Step 2. The first pairing is ⟨I⟩2\langle I\rangle^{2} (uncorrelated background), the second ∣Γ(1)(τ)∣2|\Gamma^{(1)}(\tau)|^{2} (the coherence term). Why: the first is what independent detectors register by chance, the second the real correlation surviving within τc\tau_c.

Step 3. Divide by ⟨I⟩2\langle I\rangle^{2} to normalize, giving g(2)(τ)=1+∣g(1)(τ)∣2g^{(2)}(\tau)=1+|g^{(1)}(\tau)|^{2}. Why: normalizing converts Γ(1)\Gamma^{(1)} into g(1)g^{(1)} and yields a number that starts at 22 and relaxes to 11.

Step 4. Limits: ∣g(1)(0)∣=1|g^{(1)}(0)|=1 gives g(2)(0)=2g^{(2)}(0)=2, and ∣g(1)(∞)∣=0|g^{(1)}(\infty)|=0 gives g(2)(∞)=1g^{(2)}(\infty)=1. Why: the value 22 is universal thermal bunching, independent of intensity or spectrum; the tail returns to 11 once the delay outruns τc\tau_c.

Go deeper · Derivation C · 1 minSecond-order coherence of Fock and coherent states
Derivation CSecond-order coherence of Fock and coherent states

The purpose is to compute g(2)(0)g^{(2)}(0) for the two extreme quantum states — the number state ∣n⟩|n\rangle and the coherent state ∣α⟩|\alpha\rangle — showing Fock states are antibunched at g(2)(0)=1−1/ng^{(2)}(0)=1-1/n while coherent states sit exactly at the classical boundary g(2)(0)=1g^{(2)}(0)=1. Both follow from reducing the normally ordered numerator with the single commutator [a^,a^†]=1[\hat a,\hat a^{\dagger}]=1, and together they bracket the whole coherence landscape.

Step 1. Reduce the numerator: a^†a^†a^a^=a^†(a^a^†−1)a^=n^2−n^=n^(n^−1)\hat a^{\dagger}\hat a^{\dagger}\hat a\hat a=\hat a^{\dagger}(\hat a\hat a^{\dagger}-1)\hat a=\hat n^{2}-\hat n=\hat n(\hat n-1), using [a^,a^†]=1[\hat a,\hat a^{\dagger}]=1. Why: the commutator rewrites the normally ordered product in terms of the number operator alone.

Step 2. For ∣n⟩|n\rangle the number is sharp, ⟨n^⟩=n\langle\hat n\rangle=n and ⟨n^(n^−1)⟩=n(n−1)\langle\hat n(\hat n-1)\rangle=n(n-1), giving g(2)(0)=n(n−1)/n2=1−1/ng^{(2)}(0)=n(n-1)/n^{2}=1-1/n. Why: a Fock state has zero number variance, so the factorial moment is exact.

Step 3. Read the single-photon case: n=1n=1 gives g(2)(0)=0g^{(2)}(0)=0, perfect antibunching; every finite nn gives g(2)(0)<1g^{(2)}(0)<1, and only n→∞n\to\infty recovers the classical 11. Why: this is the operational meaning of "single photon" — sub-Poissonian light no classical field can reach.

Step 4. For ∣α⟩|\alpha\rangle use a^∣α⟩=α∣α⟩\hat a|\alpha\rangle=\alpha|\alpha\rangle: numerator ∣α∣4|\alpha|^{4}, denominator ⟨n^⟩2=∣α∣4\langle\hat n\rangle^{2}=|\alpha|^{4}, so g(2)(0)=1g^{(2)}(0)=1 and indeed g(n)=1g^{(n)}=1 to every order. Why: the coherent state factorizes all correlation functions — the quantum embodiment of a perfectly coherent classical wave.

Limiting cases

The value of g(2)(0)g^{(2)}(0) classifies the source at a glance; the standard cases are:

CaseResultMeaning
Coherent (laser)g(1)(τ)=e−iω0τ, g(2)(τ)=1g^{(1)}(\tau)=e^{-i\omega_0\tau},\ g^{(2)}(\tau)=1First-order coherent for all delays and Poissonian; the reference point of the whole hierarchy.
Chaotic / thermalg(2)(0)=2, g(2)(τ)=1+∣g(1)(τ)∣2g^{(2)}(0)=2,\ g^{(2)}(\tau)=1+\vert g^{(1)}(\tau)\vert^{2}Siegert relation; photon bunching decaying from 22 to 11 over the coherence time τc\tau_c.
Single photon (Fock ∣1⟩\vert1\rangle)g(2)(0)=0g^{(2)}(0)=0Perfect antibunching — one photon cannot split between two detectors; impossible classically.
Fock state ∣n⟩\vert n\rangleg(2)(0)=1−1ng^{(2)}(0)=1-\tfrac{1}{n}Sub-Poissonian for every finite nn; approaches the classical value 11 only as n→∞n\to\infty.
Any classical fieldg(2)(0)≥1g^{(2)}(0)\ge1The classical floor — the value 11 is reached from above, never from below.

Worked examples

1 · Coherence length of a filtered spectral lamp, and its bunching (confidence builder)

A low-pressure mercury lamp is filtered to a Lorentzian lineshape of full width Δν=1.0 GHz\Delta\nu=1.0\ \mathrm{GHz} around the green 546.1 nm546.1\ \mathrm{nm} line. Find its coherence time and length, the fringe visibility at a Michelson delay equal to one coherence time, and the second-order coherence there.

Step 1. For a Lorentzian spectrum the Wiener–Khinchin partner is a two-sided exponential, ∣g(1)(τ)∣=e−πΔν∣τ∣|g^{(1)}(\tau)|=e^{-\pi\Delta\nu|\tau|}, so the 1/e1/e coherence time is τc=1/(πΔν)=1/(π×1.0×109)=0.318 ns\tau_c=1/(\pi\Delta\nu)=1/(\pi\times1.0\times10^{9})=\mathbf{0.318\ ns}. The narrower the line, the longer the phase memory.

Step 2. The coherence length is the distance light travels in that time, ℓc=c τc=(2.998×108)(3.18×10−10)=9.5 cm\ell_c=c\,\tau_c=(2.998\times10^{8})(3.18\times10^{-10})=\mathbf{9.5\ cm} — the largest Michelson path imbalance the fringes survive.

Step 3. At delay τ=τc\tau=\tau_c the modulus has fallen to ∣g(1)(τc)∣=e−1=0.368|g^{(1)}(\tau_c)|=e^{-1}=0.368, and by equation (3) that is the fringe visibility, V=0.37\mathcal V=\mathbf{0.37} — fringes washed to about a third of full contrast.

Step 4. Because the lamp is chaotic, Siegert gives the second-order coherence at that delay: g(2)(τc)=1+∣g(1)(τc)∣2=1+e−2=1.135g^{(2)}(\tau_c)=1+|g^{(1)}(\tau_c)|^{2}=1+e^{-2}=\mathbf{1.135}. The bunching, which starts at g(2)(0)=2g^{(2)}(0)=2, has decayed most of the way to 11 over one coherence time. τ_c = 0.318 ns, ℓ_c = 9.5 cm, 𝒱 = 0.37, g²(τ_c) = 1.135 — one bandwidth fixes all four.

2 · Certifying a single-photon source

A quantum-dot emitter is characterized in a Hanbury Brown–Twiss setup. Compare the ideal Fock-state coincidence values, explain why simply attenuating a laser can never mimic them, and interpret a measured g(2)(0)=0.08g^{(2)}(0)=0.08.

Step 1. From g(2)(0)=1−1/ng^{(2)}(0)=1-1/n, the ideal number states give ∣1⟩→0|1\rangle\to0, ∣2⟩→0.5|2\rangle\to0.5, ∣5⟩→0.8|5\rangle\to0.8, ∣10⟩→0.9|10\rangle\to0.9. A perfect single-photon source reads exactly 0\mathbf{0}; even a two-photon state already sits at 0.50.5, on the classical boundary. Antibunching weakens fast with photon number.

Step 2. Attenuating a bright laser does not help. A coherent state has g(2)(0)=1g^{(2)}(0)=1 regardless of ∣α∣|\alpha|, so dimming it to one photon per pulse still gives g(2)(0)=1g^{(2)}(0)=1 — the weak beam is mostly vacuum with an occasional Poissonian burst, including two-photon events. Weakness is not antibunching; a faint laser is still classical.

Step 3. The certification threshold is g(2)(0)<0.5g^{(2)}(0)<0.5, which rules out any classical mixture faking single photons by emitting two-photon states half the time. A measurement of g(2)(0)=0.08g^{(2)}(0)=0.08 is comfortably below it.

Step 4. As a purity, g(2)(0)≈0.08g^{(2)}(0)\approx0.08 means two-photon emission is suppressed to roughly 8%8\% of the Poissonian level — good enough for quantum key distribution, where multiphoton pulses are the security leak eavesdroppers exploit. Ideal |1⟩ gives 0; an attenuated laser stays at 1; the measured 0.08 certifies a near-ideal single-photon source.

Physics of the portrait

Coherence functions are where classical optics, invented to describe fringes, was forced to become quantum — and the person who forced it, giving the degree of coherence its modern operator definition, is this article's natural face.

Origins · 1 minRoy Glauber and how this was discovered
RG
Roy Glauber
1925–2018 · founder of quantum optics; defined the coherence functions and the degree of coherence

Roy J. Glauber built the quantum theory of optical coherence in a pair of 1963 papers that answered a crisis the Hanbury Brown–Twiss experiment had opened. When Hanbury Brown and Twiss reported in 1956 that photons from a thermal source arrive in correlated clumps, many physicists doubted the result could survive a photon picture at all. Glauber resolved it by defining the coherence functions G(n)G^{(n)} as normally ordered expectation values of the field operators, arguing that because photodetection absorbs photons the theory must be written in normal order. In the same work he introduced the coherent states ∣α⟩|\alpha\rangle as eigenstates of the annihilation operator and showed they factorize every G(n)G^{(n)} — the precise quantum meaning of "fully coherent." His framework explained thermal bunching, predicted the antibunching that only quantized light can show, and turned the vague adjective "coherent" into a measurable ladder of numbers. He received the 2005 Nobel Prize in Physics for this theory, half a century after the experiment that provoked it.

History
1934
— van Cittert computes the coherence of light from an extended source
1938
— Zernike formalizes the degree of coherence and partial coherence
1956
— Hanbury Brown and Twiss measure intensity correlations in starlight from Sirius
1963
— Glauber builds the quantum theory of optical coherence from normally ordered field operators
1977
— Kimble, Dagenais and Mandel observe photon antibunching, the first g² below one

Going further

Texts. Loudon, The Quantum Theory of Light, 3rd ed., Ch. 3 and 6 (coherence functions, photon statistics) · Mandel and Wolf, Optical Coherence and Quantum Optics (1995), the comprehensive reference · Foot, Atomic Physics, Ch. 8 (coherence and the HBT effect) · Saleh and Teich, Fundamentals of Photonics, 3rd ed., Ch. 11–12 · Hecht, Optics, 5th ed., Ch. 12 (classical coherence) · Fox, Quantum Optics, Ch. 5–6. The key results to keep in view are: (i) V=∣g(1)∣\mathcal V=|g^{(1)}| from two-beam interference; (ii) the Siegert relation g(2)(τ)=1+∣g(1)(τ)∣2g^{(2)}(\tau)=1+|g^{(1)}(\tau)|^{2} from the Gaussian moment theorem, giving g(2)(0)=2g^{(2)}(0)=2; (iii) g(2)(0)=1−1/ng^{(2)}(0)=1-1/n for Fock states, yielding 0, 0.5, 0.8, 0.90,\,0.5,\,0.8,\,0.9 for n=1,2,5,10n=1,2,5,10; (iv) thermal g(2)(0)=2g^{(2)}(0)=2 from ⟨n^(n^−1)⟩=2nˉ2\langle\hat n(\hat n-1)\rangle=2\bar n^{2}; (v) coherent-state g(n)=1g^{(n)}=1 to all orders; (vi) the filtered-lamp numbers τc=0.318 ns\tau_c=0.318\ \mathrm{ns}, ℓc=9.5 cm\ell_c=9.5\ \mathrm{cm}, V=0.37\mathcal V=0.37, g(2)(τc)=1.135g^{(2)}(\tau_c)=1.135. The standard values g(2)(0)=1,2,0g^{(2)}(0)=1,2,0 for laser, thermal, and single-photon light are consistent across the treatments listed above.

Prerequisites. Quantization of the Electromagnetic Field · Photon Statistics and the Mandel Q Parameter · Coherent States. Leads to. Hanbury Brown–Twiss Experiment · Photon Antibunching · Single-Photon Sources · Phase-Space Quasiprobability Distributions.

Ask the tutor: Ask me to walk through why a filtered thermal lamp and a laser can share the same fringes yet split apart the instant you count coincidences.

What is the difference between first-order and second-order coherence?

First-order coherence, g(1)(τ)g^{(1)}(\tau), measures how well the field's amplitude and phase now predict its amplitude and phase a delay τ\tau later; its modulus is fringe visibility and its decay time is the coherence time. Second-order coherence, g(2)(τ)g^{(2)}(\tau), measures how well the intensity now predicts the intensity later; it is the normalized rate at which two detectors click together, telling you whether photons arrive bunched, random, or antibunched. The first is about waves and interference, the second about particles and photon statistics — and the point of the article is that they are independent.

Why can a laser and filtered thermal light have the same g¹ but different g²?

Because g(1)g^{(1)} is set by the spectrum and g(2)g^{(2)} by the photon statistics — separate properties. Wiener–Khinchin makes g(1)(τ)g^{(1)}(\tau) the Fourier transform of the lineshape, so filtering a thermal lamp to a laser's narrow bandwidth gives an identical g(1)g^{(1)} and identical fringes. But filtering does not touch how photons are spread in time: the thermal beam stays bunched at g(2)(0)=2g^{(2)}(0)=2, the laser Poissonian at g(2)(0)=1g^{(2)}(0)=1. An interferometer cannot see the difference; only a Hanbury Brown–Twiss coincidence measurement can.

Is the Siegert relation always true?

No — this is a common trap. The Siegert relation g(2)(τ)=1+∣g(1)(τ)∣2g^{(2)}(\tau)=1+|g^{(1)}(\tau)|^2 holds only for chaotic (thermal) light, whose complex-Gaussian statistics let the Gaussian moment theorem factorize the fourth-order moment and force g(2)(0)=2g^{(2)}(0)=2. For a laser it fails — coherent light has ∣g(1)∣=1|g^{(1)}|=1 yet g(2)=1g^{(2)}=1, not 22 — and for nonclassical light it fails even harder, since g(2)(0)g^{(2)}(0) can be 00. It is a property of one statistical class of light, not a universal law.

What does g²(0) < 1 actually mean physically?

It means the photons actively avoid arriving together — antibunching — which requires the number variance to fall below its Poissonian value, ⟨(Δn^)2⟩<⟨n^⟩\langle(\Delta\hat n)^2\rangle<\langle\hat n\rangle. No classical field can do this: a classical intensity is real and non-negative, so its variance is non-negative and forces g(2)(0)≥1g^{(2)}(0)\ge1. A sub-mean variance needs a quantized field with a nearly definite photon number — a Fock state, or a single atom that emits one photon and must re-excite before it can emit again. A measured g(2)(0)<1g^{(2)}(0)<1 is therefore an unforgeable certificate that light is quantized.

Why is normal ordering forced, rather than a convention?

Because a photodetector registers a count by absorbing a photon, the action of a^\hat a. A coincidence between two detectors is two absorptions, an expectation value ⟨a^†a^†a^a^⟩\langle\hat a^{\dagger}\hat a^{\dagger}\hat a\hat a\rangle with all creation operators to the left — normal order. Any other ordering mixes in the field's zero-point fluctuations and predicts spurious counts from the vacuum, which clicks nothing. Normal ordering removes that vacuum term, so an empty mode gives zero coincidences and g(n)g^{(n)} reports only real photon correlations — imposed by the physics of detection, not chosen for convenience.

How does the coherence time relate to the linewidth of a source?

They are Fourier reciprocals. Through Wiener–Khinchin, g(1)(τ)g^{(1)}(\tau) is the Fourier transform of the spectrum, so a bandwidth Δν\Delta\nu gives a coherence that decays over τc∼1/Δν\tau_c\sim1/\Delta\nu. A broad white source loses phase memory almost instantly (femtosecond τc\tau_c), while a stabilized laser of hertz-level linewidth stays coherent for a fraction of a second. For a Lorentzian line, τc=1/(πΔν)\tau_c=1/(\pi\Delta\nu): the worked example takes a 1 GHz1\ \mathrm{GHz} line to τc=0.318 ns\tau_c=0.318\ \mathrm{ns} and ℓc=9.5 cm\ell_c=9.5\ \mathrm{cm}. Narrowing a spectrum lengthens its coherence.

Before this

Where this leads