An Interactive Guide

The Kramers–Kronig Relations

Causality in time becomes analyticity in frequency. Absorption and dispersion are two sides of the same complex function — one cannot change without the other.

  χ'(ω)  =   (2/π) P ∫₀ⁿ ω' χ''(ω') / (ω'² − ω²) dω'
  χ''(ω) = −(2ω/π) P ∫₀ⁿ χ'(ω') / (ω'² − ω²)      dω'
Scroll to explore

Two halves of a complex response function are not independent

Any causal, linear, time-invariant response — an electric susceptibility, a refractive index, a loudspeaker transfer function, a microwave S-parameter — is a complex-valued function of frequency. The Kramers–Kronig relations state that the real and imaginary parts of such a response cannot be chosen independently. Each is the Hilbert transform of the other.

The same theorem appears in electromagnetism, dielectric spectroscopy, network analysis, rheology, seismic attenuation, filter theory, room acoustics, and fibre-optic communications. Every discipline dresses it differently, but the mathematics is one.

One Sentence Version

If the system cannot respond before it is poked, then measuring how much it absorbs determines how much it disperses — and vice versa — up to an integral over all frequencies.

Why the relations exist at all

Write a causal impulse response as h(t) with h(t)=0 for t<0. Its Fourier transform (physics convention) is

χ(ω) = ∫₀ⁿ h(t) eiωt dt

The frequency variable ω is real on physical grounds — it is the angular frequency of a steady-state Fourier component — but the integral on the right is a perfectly good function of ω regarded as a complex variable. Promoting ω = ωr + iωi with ωr, ωi ∈ ℝ turns χ(ω) into a candidate analytic function on the complex plane, and the K–K argument lives or dies on where that candidate is well-defined.

For Im ω > 0 — i.e. ωi > 0, the upper half of the complex frequency plane — the factor eiωt = eiωrt e−ωit decays exponentially as t → ∞, so the integral converges absolutely (causality, h(t)=0 for t<0, is what lets us drop the lower half of the time axis where this exponential would blow up). A convergent power-series definition in an open region is holomorphic. Therefore χ is analytic in the upper half complex-frequency plane. Apply Cauchy's integral formula on a closed contour along the real axis (with a small indentation at ω=ω₀) plus a large arc at infinity, and the two Kramers–Kronig integrals fall out.

The K–K integrals themselves only ever evaluate χ on the real axis (real ω), which is what experiments measure. The complex extension is scaffolding: it exists so that Cauchy's theorem has somewhere to live, and the conclusions are then read off on its real boundary.

The same logic in the engineering sign convention e+jωt swaps the role of upper and lower half planes. The content is identical; only the sign of an i flips.

Interactive · Causal vs Non-causal Impulse Response

Three columns × two rows. Top row shows real parts; bottom row shows imaginary parts. Left: complex impulse response h(t) = e−γ(t−τ) eiω₀(t−τ) for t ≥ τ, zero otherwise. Centre: real and imaginary parts of its spectrum χ(ω) — the “truth”. Right: the K–K prediction (gold) of each part computed from the other one via the Hilbert transform — Re χ = H[Im χ] on top, Im χ = −H[Re χ] on the bottom. The right column should match the centre column when h is causal (τ ≥ 0) and visibly disagree as soon as h leaks into negative time.

Four equivalent forms of the same relations

One-sided form, real h(t)

χ'(ω) = (2/π) P ∫₀ⁿ ω' χ''(ω') / (ω'² − ω²) dω' χ''(ω) = −(2ω/π) P ∫₀ⁿ χ'(ω') / (ω'² − ω²) dω'

Two-sided form, general complex spectrum

χ'(ω) = (1/π) P ∫−∞∞ χ''(ω') / (ω' − ω) dω' χ''(ω) = −(1/π) P ∫−∞∞ χ'(ω') / (ω' − ω) dω'

Hilbert-transform form

χ' = H[ χ'' ] χ'' = −H[ χ' ]

Once-subtracted form (for χ that does not vanish at ∞)

χ(ω) − χ(ωa) = (ω − ωa)/(iπ) · P ∫ χ(ω') / [(ω' − ω)(ω' − ωa)] dω'

The P is the Cauchy principal value — it handles the pole on the real axis. The subtracted form is essential for optical constants where n(∞)→1, not zero.

The Lorentzian oscillator — K–K in one line

A single damped harmonic oscillator driven by an electric field gives a susceptibility with one pair of complex-conjugate poles:

χ(ω) = F / (ω₀² − ω² − iγω)

Splitting χ = χ' + iχ'' gives

χ'(ω) = F (ω₀² − ω²) / [ (ω₀² − ω²)² + γ²ω² ] χ''(ω) = F γω / [ (ω₀² − ω²)² + γ²ω² ]

You can verify by direct integration that these satisfy K–K. Because the poles sit in the lower half plane (physics convention), χ(ω) is analytic in the upper half plane — just as causality demands.

Interactive · Lorentzian χ(ω) in Different Plot Styles

Same complex function, three conventions. The electrical-engineer's Bode plot, the dielectric-spectroscopist's Cole–Cole semicircle, and the physicist's χ'/χ'' spectra all carry identical information.

The contour that unlocks the proof

The Kramers–Kronig derivation is a single contour integral. Walk along the real axis from −∞ to +∞, indent with a small semicircle above the pole at ω₀, then close the loop with a giant semicircle in the upper half plane. The enclosed region contains no singularities, so Cauchy says the total integral is zero.

Diagram · The Kramers–Kronig Contour

Large arc at R→∞ vanishes by Jordan's lemma. Small indent at ω₀ contributes −iπ χ(ω₀) by the Sokhotski–Plemelj formula (half of a full residue). The rest is the principal-value integral along ℝ.

Rearranging:

χ(ω₀) = (1/iπ) P ∫−∞∞ χ(ω') / (ω' − ω₀) dω'

Split into real and imaginary parts and you have the relations.

Same theorem, different skins

One of the reasons the Kramers–Kronig literature is so hard to read across fields is that every community chose its own notation. Here is the translation table.

DisciplineQuantityDecompositionTypical plot
Optics / condensed matterε(ω), ñ(ω) = n + iκdispersion + absorptionn, κ vs ℏω
Dielectric spectroscopyε* = ε' − jε'' (eng) or ε' + iε'' (phys)storage + lossCole–Cole semicircle
Microwave / RFSij(ω), Z(ω)real & imag, or mag & phaseSmith chart, Bode
Controls / signalsH(jω) = |H|ejφgain & phaseBode, Nyquist
RheologyG*(ω) = G' + iG''elastic + loss modulilog–log vs ω
Particle physicsf(ω) scattering amplitudereal + i·(σtot/4π)forward amplitude vs energy
Audio / DSPH(ejω)log-magnitude & phasedB / linear phase
Sign-convention trap

Physics writes ε = ε' + iε'' with ε'' > 0 for absorption. Some engineering texts write ε = ε' − jε'' with the same ε'' > 0. Both describe a lossy medium — the sign flip lives in e−iωt vs e+jωt. Always check the kernel before copying a formula across a disciplinary boundary.

Where engineers actually use it

DSP
Minimum-phase reconstruction
Given |H(ejω)|, the K–K relations on log H recover phase uniquely — the discrete-time Bode theorem. Core of the complex cepstrum.
RF / SI
Causality enforcement on VNA data
Measured S-parameters are band-limited and noisy. K–K post-processing yields a physically realisable macromodel for SPICE or IBIS-AMI.
OPTICS
Reflectance to refractive index
Bulk samples give only R(ω); K–K on log r recovers the phase, and Fresnel inverts to n + iκ. Workhorse of solid-state optics.
EM SIM
Causal dispersive materials in FDTD
Fit measured ε(ω) to Lorentzian / Debye / Drude poles — poles in the lower half plane are automatically K–K-consistent, and yield stable ADE or PLRC updates.
COMMS
Group delay — phase bundles
Channel magnitude is not free of phase. Equalisers exploiting K–K give minimum group delay; non-minimum-phase channels demand latency.
AUDIO
Loudspeaker & room EQ
Minimum-phase drivers respond to IIR magnitude EQ without added delay; rooms decompose as minimum-phase × allpass, only the first is K–K-invertible.
MATERIALS
Broadband dielectric consistency
Spectra measured from mHz to GHz are cross-checked: if ε'(ω) and ε''(ω) don't obey K–K, something in the calibration is wrong.
FIBRE OPTICS
K–K coherent receivers (2016–)
A single-ended photodiode detects intensity only; a minimum-phase field plus K–K recovers the optical phase, eliminating one coherent-detector branch.

Chapter index

Each chapter is a standalone Markdown page in the chapters/ directory. The five parts build from mathematical foundations, through the theorem itself, to the applications above.

Part I · Mathematical Foundations
Part II · The Kramers–Kronig Relations
Part III · Conventions
Part IV · Applications
Part V · Practice & Further Reading