Fourier-Transform Spectroscopy
Fourier-transform spectroscopy turns the ordinary spectrometer inside out. A grating or prism fans the colours across space and reads them channel by channel; a Fourier-transform spectrometer refuses to separate the colours at all. Instead it splits a single beam in a scanning Michelson interferometer, lets the two copies interfere, and records one number — the detected brightness — as the path difference is swept. The whole spectrum is folded into that single scan, the interferogram, and a Fourier transform unfolds it back into colour. No slit, no ruling, no moving grating: just one mirror sliding through a few centimetres while a detector watches the fringes bloom and fade. The reason this works is the two-beam interference term. When a field is split, delayed by a time , and recombined, the modulated part of the output is the field's own autocorrelation — the very quantity that measures temporal coherence. The Wiener–Khinchin theorem then supplies the punchline: the Fourier transform of that autocorrelation is the power spectrum. So the interferogram is not a picture of the spectrum but its Fourier transform, and reading one off the other is a single numerical step. This one idea makes the infrared spectrum of almost every molecule on Earth accessible in seconds, and it carries two structural bonuses — the multiplex and throughput advantages — that a slit spectrometer can never match.
- Fourier-transform spectroscopy (FTS) recovers a spectrum not by dispersing light in space but by recording an interferogram: the output of a scanning Michelson interferometer versus optical path difference .
- The modulated interferogram is the real part of the field autocorrelation , — the same degree of temporal coherence that governs fringe visibility.
- By the Wiener–Khinchin theorem, the Fourier transform of that autocorrelation is the power spectrum, so — a delay scan replaces a grating and every colour is measured at once.
- The maximum path difference sets the resolution, ; the sampling step sets the Nyquist limit against aliasing.
- The multiplex (Fellgett) and throughput (Jacquinot) advantages give FTS a large signal-to-noise edge in the infrared.
- Key idea: the interferogram of a scanning Michelson is the field's autocorrelation; its Fourier transform with respect to path delay is the source's power spectrum, and the scan length alone fixes the resolution at δσ = 1/(2Δmax).
A two-beam interferometer measures the autocorrelation of the light, not the spectrum. Sweep the path difference to record the interferogram ; subtract its mean and Fourier-transform it, and the Wiener–Khinchin theorem returns the spectrum . Because the transform pair is autocorrelation ↔ spectrum, narrow lines make long interferograms and broad bands make short bursts — and the longest delay you reach, , is exactly what limits how finely you can resolve.
Where this lives today
Fourier-transform spectroscopy is the default way to take an infrared spectrum. Virtually every chemistry lab runs an FTIR: a broadband source, a Michelson with a potassium-bromide or diamond beamsplitter, and a helium–neon reference laser deliver a molecule's vibrational fingerprint across roughly 400–4000 cm⁻¹ in seconds, with the Fellgett and Jacquinot advantages supplying the sensitivity to see thin films, trace gases, and reactions as they run. The same architecture, scaled up, flies in space: atmospheric sounders such as IASI and climate instruments measuring CO₂, CH₄, and H₂O rely on large maximum path differences for the fine resolution needed to pick apart rotational–vibrational lines. The idea also generalises far beyond optics. Fourier-transform nuclear magnetic resonance and Fourier-transform ion-cyclotron-resonance mass spectrometry record a time-domain signal and transform it to a spectrum by exactly the same multiplex logic, inheriting the same window-length and Nyquist rules. And because the interferogram envelope is the degree of temporal coherence, an FTS scan doubles as a coherence-length meter: a Lorentzian decay reveals lifetime or collision broadening, a Gaussian decay reveals Doppler broadening — the line-broadening mechanism read straight off the fading of the fringes.
The physical picture
Imagine building a spectrum out of pure tones. Each optical frequency , passing through a two-beam interferometer, writes a single cosine into the path-difference axis. A real source is a chorus of such tones with weights ; the detector adds all their cosines together, and the result is the interferogram. Reverse the operation — a Fourier transform — and the weights fall back out. The interferometer is a physical Fourier synthesiser, and the computer that transforms its output is the analyser that reads the score back. At zero path difference every tone is in step and they pile up into a tall central burst. As the mirror moves, the tones drift out of phase and the fringes wash out — quickly for a broad spectrum, slowly for a narrow one. That decay is the loss of temporal coherence, and its width in delay is the coherence length. Narrow in frequency means wide in delay, and vice versa: the reciprocity that runs through all of Fourier optics.
A reference laser makes all of this metrology rather than art. A helium–neon beam runs through the same interferometer and marks off equal steps of path difference by its own single fringe; every sample of the interferogram is therefore taken at a path difference known to a fraction of a wavelength. FTS inherits its wavenumber accuracy directly from that reference wavelength — which is why the method is prized not only for sensitivity but for the sheer trustworthiness of its frequency axis. The single most important number in the whole method is how far the mirror travels. The finite scan truncates the interferogram, and truncation blurs the recovered spectrum. The next figure makes the rule concrete: resolution and scan length are simple reciprocals.
The core equations
Start from what the detector actually sees. A 50:50 beamsplitter makes two copies of the field, one delayed by , and the detector time-averages the squared sum. The fast optical oscillation averages away, leaving a constant plus one modulated term — and that term is the field talking to a delayed copy of itself.
where: — detected intensity at optical path difference ; — the field autocorrelation (self-coherence) function; — the mean intensity (the DC level); — the time delay between the two arms. The measurement is an autocorrelation of the light, full stop; the constant carries no colour information and is subtracted before transforming. The bridge from that autocorrelation to the spectrum is the Wiener–Khinchin theorem. For a stationary field the autocorrelation and the power spectral density are a Fourier-transform pair — the spectrum is not measured, it is computed.
where: — power spectral density (the source spectrum); — optical frequency. Reading left to right builds the interferogram from the spectrum; reading right to left is what the spectrometer's computer does. Everything else is bookkeeping on this one pair. In practice the field is real, the spectrum one-sided, and the interferogram symmetric, so the complex transform collapses to a cosine transform in the directly measured path difference . Working in wavenumber (the spectroscopist's unit, cm) removes the factor of entirely.
where: — wavenumber; — the subtracted mean level; — optical path difference, twice the mirror displacement. The cosine transform is its own inverse up to a constant, so the same operation runs in both directions. The scan cannot be infinite, and the truncation is what limits resolution. Cutting the interferogram off at convolves the true spectrum with a sinc instrument line shape whose width sets the smallest resolvable interval.
where: — wavenumber resolution (cm for in cm); — frequency resolution; — maximum optical path difference. Resolution is pure geometry — how far the mirror travels — and owes nothing to the detector. Sampling at steps then sets the Nyquist ceiling above which frequencies alias. Finally, the reason FTS conquered the infrared: because every spectral element is measured throughout the whole scan rather than one at a time, a detector-noise-limited FTS beats a scanning monochromator by the square root of the number of elements.
where: — number of resolved spectral elements measured simultaneously; the equality is the Fellgett (multiplex) advantage and holds when additive detector noise dominates.
The goal is to show, from scratch, why a scanning two-beam interferometer measures a Fourier transform of the spectrum rather than the spectrum itself, and therefore why a single numerical transform of the recorded interferogram returns the colours. We build the detected intensity from the superposed fields, isolate the modulated part as the field autocorrelation, and invoke the Wiener–Khinchin theorem to identify its transform as the power spectrum — the exact chain equations (1)–(3) assert.
Step 1. Split the field into two equal copies and delay one by , so the recombined field is ; the detector reports the time average . Why: a square-law detector integrates over millions of optical cycles, so every term oscillating at the optical frequency averages to zero except the correlation between the field and its delayed self.
Step 2. Recognise the surviving term as the autocorrelation , so the modulated interferogram is exactly . Why: this is the same self-coherence function that sets fringe visibility, so the interferogram literally traces the light's degree of temporal coherence as the delay grows.
Step 3. Expand the field in its spectral components, , and use stationarity, , to collapse the double integral to . Why: a stationary source cannot correlate different frequencies, so the delta function reduces the correlation of the whole field to a single integral weighted by the spectrum — this is the Wiener–Khinchin theorem in the making.
Step 4. Take the real part and invert the cosine transform: , so with . Why: the cosine transform is self-inverse up to a constant, so the very operation that built the interferogram from the spectrum, run backwards on the measured interferogram, hands back the spectrum. *Verified: for a monochromatic input the interferogram reduces to , and a numerical cosine transform of a two-line (doublet) spectrum reproduces the beating interferogram and recovers both line positions, in verify_aq-optics-fourier-transform-spectroscopy.py.*
Go deeper · Derivation B · 1 minWhy the scan length fixes the resolution
The aim here is to pin down the single most consequential specification of a Fourier-transform spectrometer — its resolving power — and to show it is set entirely by how far the mirror travels, not by the detector or the source. We treat the finite scan as a rectangular window on the interferogram and track what that truncation does to the recovered spectrum, arriving at and the Nyquist sampling condition.
Step 1. Write the finite scan as the ideal (infinite) interferogram multiplied by a boxcar of half-width : . Why: no instrument scans to infinite delay, and the boxcar captures exactly what is thrown away — everything beyond the last mirror position.
Step 2. Transforming a product gives a convolution, so the recovered spectrum is the true spectrum smeared by the transform of the boxcar, with . Why: truncation in delay must blur in frequency; the sharper you want the spectrum, the longer you must scan — the two are Fourier reciprocals.
Step 3. Read the width of that sinc: its first zero sits at , which is the smallest interval over which two lines stay distinct (the Rayleigh-like criterion), so and . Why: two lines closer than the first zero merge under the instrument line shape, so the sinc half-width is the resolution — a pure consequence of the scan length.
Step 4. Sample the interferogram at steps ; the recovered spectrum then repeats with period , so to avoid overlap the spectrum must lie below (equivalently ). Why: discrete sampling is periodic in frequency, so anything above the Nyquist wavenumber folds back into the band — a reference laser fixes precisely to keep that ceiling under control. *Verified: cm⁻¹ (equivalently GHz) at cm, the cm needed to split lines 1 cm⁻¹ apart, and the Nyquist relation , in verify_aq-optics-fourier-transform-spectroscopy.py.*
Limiting cases
| Case | Result | Meaning |
|---|---|---|
| Monochromatic source | a single line gives an undamped cosine; its transform is a δ-spike at | |
| Lorentzian line, width | exponential envelope; coherence length read straight off the decay | |
| Gaussian (Doppler) line | Gaussian envelope; the transform of a Gaussian is a Gaussian, recovering the Doppler profile | |
| Finite scan | rectangular truncation → sinc instrument line shape with sidelobes | |
| Sampling step | frequencies above the Nyquist limit alias; a bandpass or fine sampling prevents fold-over | |
| Multiplex, elements | Fellgett gain in the detector-noise limit; vanishes (or reverses) when photon shot noise dominates |
Worked examples
1 · Resolution of an FTIR spectrometer (confidence builder)
A benchtop FTIR scans to a maximum optical path difference of cm. What wavenumber resolution does it achieve, what is that in frequency, and how far must it scan to separate two lines 1 cm⁻¹ apart?
Step 1. Resolution is the reciprocal of twice the scan length: . That is fine enough to resolve the rotational fine structure on a vibrational band of a small gas-phase molecule.
Step 2. Convert to frequency: . The same instrument therefore places optical features to better than a part in of a typical Hz optical frequency.
Step 3. Invert the resolution rule for the required scan: to reach cm⁻¹ you need . Two lines 1 cm⁻¹ apart are just resolved by a mere half-centimetre of travel.
Step 4. The scaling is the design lever: quadrupling the wavenumber resolution to 0.0125 cm⁻¹ demands a cm scan — a metre-scale optical bench and correspondingly stringent mirror alignment. Resolution in FTS is bought with length, linearly. δσ = 0.05 cm⁻¹ (1.5 GHz) at 10 cm; 0.5 cm scan splits 1 cm⁻¹.
2 · The Fellgett multiplex advantage (exam level)
A mid-infrared spectrum is to be recorded across resolution elements. Compare the signal-to-noise of a Fourier-transform spectrometer with a scanning grating monochromator that visits one element at a time, assuming the same total measurement time and detector-noise-limited operation.
Step 1. The monochromator spends only of the run on each element, gathering signal from one channel while the other go unmeasured. The FTS detector sees light from all elements throughout the entire scan — the multiplex principle.
Step 2. With additive detector noise independent of signal, the noise is the same for both instruments over equal time, while the FTS accumulates every element's photons the whole time. The signal-to-noise ratios then differ by the square root of the number of elements: .
Step 3. Put in the number: . The Fourier-transform spectrometer delivers about a signal-to-noise improvement — or reaches the same quality in of the time.
Step 4. The advantage is conditional. It is largest in the infrared, where cheap thermal detectors are detector-noise-limited; it shrinks toward the visible, and under photon shot noise from a bright source it can vanish or even reverse (the multiplex disadvantage), because piling all the light onto one detector piles up its shot noise too. √1000 ≈ 32× SNR, detector-noise-limited only.
Physics of the portrait
The whole method rests on one mathematical fact — that the autocorrelation of a signal and its power spectrum are a Fourier-transform pair — and that fact has a name.
Origins · 1 minWiener–Khinchin theorem and how this was discovered
The Wiener–Khinchin theorem states that for a wide-sense-stationary random signal the power spectral density is the Fourier transform of the autocorrelation function. Norbert Wiener, the American founder of cybernetics, established it in his 1930 memoir Generalized Harmonic Analysis, giving a rigorous spectral theory for signals — like light from a thermal source or noise in a circuit — that have no ordinary Fourier transform because they never decay. The Soviet probabilist Aleksandr Khinchin gave the theorem its clean form for stationary stochastic processes in 1934, and the joint name has stuck. Its reach is enormous: it is why a noisy resistor has a definable spectrum (Johnson–Nyquist noise), why radio astronomers can build spectra from correlator outputs, and — the case here — why a two-beam interferometer, which physically measures the autocorrelation of an optical field, delivers the source's spectrum after one Fourier transform. The theorem is the mathematical licence for the entire enterprise of Fourier-transform spectroscopy: it guarantees that the interferogram and the spectrum carry identical information, just in conjugate variables.
Going further
References. Hecht, Optics, 5th ed., §7.4 and §12.1 (coherence, Fourier methods) · Loudon, The Quantum Theory of Light, 3rd ed., Ch. 3 (coherence and correlation functions) · Saleh & Teich, Fundamentals of Photonics, §11.3 (Fourier spectroscopy) · Demtröder, Laser Spectroscopy, Vol. 1, §4.4 (Fourier-transform spectroscopy) · Wiener (1930) and Khinchin (1934), the autocorrelation–spectrum theorem.
Provenance. Equations follow from the two-wave interference term and the Wiener–Khinchin theorem as in Hecht §12.1, Loudon Ch. 3, Saleh & Teich §11.3, and Demtröder §4.4. Recomputed and verified here: (i) the interferogram from superposition; (ii) the cosine-transform inversion , checked numerically on a doublet; (iii) resolution cm⁻¹ (1.5 GHz) at cm; (iv) the Nyquist limit ; (v) the Fellgett scaling in the detector-noise limit. Line-shape envelopes (Lorentzian , Gaussian) are standard Fourier pairs; no non-standard numerical constants are asserted.
Builds on. Two-Wave Interference and the Interference Term · Temporal Coherence and Coherence Length · Mutual Coherence Function and Degree of Coherence · Optical Path Difference and Phase Retardation · Fringe Visibility and Contrast. Leads to. Wavefront-Splitting and Amplitude-Splitting Interferometers · Spatial Coherence and the Van Cittert–Zernike Theorem · Conditions for Observable Interference.
Is the interferogram the same thing as the spectrum, just drawn differently?
No — they are Fourier transforms of each other, which means they carry the same information in conjugate variables but look nothing alike. The horizontal axis of the interferogram is optical path difference (a length), while the spectrum's axis is frequency or wavenumber. A single spectral line, sharp in frequency, spreads into a cosine that fills the whole interferogram; a broad band, wide in frequency, collapses into a narrow burst near zero path difference. You cannot read a spectral feature off the interferogram by eye; the transform is essential.
Why does scanning the mirror farther improve resolution?
Because truncating the interferogram at a maximum path difference is mathematically a rectangular window, and windowing in delay convolves the recovered spectrum with a sinc whose width is . Two lines closer than that width blur into one. The farther you scan, the narrower the sinc, and the finer the detail you can separate — so resolution and scan length are exact reciprocals, . Nothing about the detector or the source changes this; only geometry does.
What is the reference laser for?
It supplies the ruler for the path-difference axis. A helium–neon laser passes through the same interferometer, and its single sharp fringe marks off equal, precisely known increments of path difference; the interferogram is sampled at each of those marks. This does two jobs: it fixes the sampling step (and hence the Nyquist limit) to a fraction of the laser wavelength, and it pins the absolute wavenumber scale of the final spectrum to the reference wavelength. FTS owes its exceptional frequency accuracy almost entirely to this trick.
Why is real data apodized, and what does that cost?
A finite rectangular scan gives a sinc instrument line shape whose sidelobes reach about of the main peak — ripples that can be mistaken for weak spectral lines. Apodization multiplies the interferogram by a smooth taper (triangular, Happ–Genzel, Norton–Beer) that goes gently to zero at the scan ends, which suppresses those sidelobes. The price is a broader main lobe: apodization trades a modest loss of resolution for a much cleaner, artifact-free line shape, and the choice of taper is a deliberate compromise for the measurement at hand.
When does the multiplex advantage disappear?
The Fellgett advantage assumes the noise is additive and independent of the signal — true for thermal infrared detectors, where it delivers its full gain. If instead the dominant noise is photon shot noise from a bright source, the noise grows with the total light on the detector; piling all elements onto one detector then piles up their shot noise too, and the advantage shrinks to nothing or reverses into a multiplex disadvantage. This is why FTS dominates the infrared but competes on more even terms with dispersive and, especially, modern array-detector instruments in the visible and ultraviolet.
Why can an FTS interferogram be asymmetric, and how is that fixed?
Ideally the interferogram is an even function of path difference and a pure cosine transform returns the spectrum. But wavelength-dependent phase errors — from a dispersive beamsplitter, imperfect compensation, or electronic delays — make the autocorrelation complex, so the interferogram is no longer symmetric. Transforming it naively distorts the intensities. The remedy is to record a double-sided interferogram (or a short double-sided plus a long single-sided one) and apply a phase correction, such as the Mertz or Forman method, taking the magnitude of the complex transform so the recovered line intensities are undistorted.