Skip to content

Fluctuation-Dissipation Theorem

The fluctuation-dissipation theorem (FDT) connects the response of a system to a small perturbation with the spontaneous fluctuations of the system at equilibrium.

Consider a Hamiltonian H0\mathcal{H}_0 perturbed by a time-dependent field:

H(t)=H0f(t)A\mathcal{H}(t) = \mathcal{H}_0 - f(t)A

Where AA is an observable conjugate to the field f(t)f(t). The change in A(t)\langle A(t) \rangle to first order in ff is:

A(t)A0=tχAA(tt")f(t)dt\langle A(t) \rangle - \langle A \rangle_0 = \int_{-\infty}^{t} \chi_{AA}(t - t")\, f(t')\, dt'

Where the response function is:

χAA(t)=iθ(t)[A(t),A(0)]0\chi_{AA}(t) = \frac{i}{\hbar}\theta(t)\langle[A(t), A(0)]\rangle_0

In the classical limit, the FDT takes a simpler form. The dynamic susceptibility χ(ω)=χ(ω)+iχ(ω)\chi(\omega) = \chi'(\omega) + i\chi''(\omega) relates to the power spectrum S(ω)S(\omega) of fluctuations:

S(ω)=2kBTωχ(ω)S(\omega) = \frac{2k_B T}{\omega}\,\chi''(\omega)

For a harmonic oscillator with damping γ\gamma and natural frequency ω0\omega_0:

χ(ω)=γω(ω02ω2)2+γ2ω2\chi''(\omega) = \frac{\gamma\omega}{(\omega_0^2 - \omega^2)^2 + \gamma^2\omega^2}

The fluctuation spectrum is Lorentzian, peaked at ω0\omega_0.

The FDT predicts thermal (Johnson—Nyquist) noise in a resistor:

V2=4kBTRΔf\langle V^2 \rangle = 4k_B T R \Delta f

Where RR is the resistance and Δf\Delta f is the bandwidth. This noise is fundamental — it arises from thermal fluctuations of charge carriers and cannot be eliminated.

Worked Example 13.1: Johnson--Nyquist Noise Calculation

A 1010 kΩ\Omega resistor at room temperature (T=300T = 300 K) measured with bandwidth Δf=1\Delta f = 1 MHz:

V2=4×1.38×1023×300×104×106\langle V^2 \rangle = 4 \times 1.38 \times 10^{-23} \times 300 \times 10^4 \times 10^6

=4×1.38×1023×3×1012= 4 \times 1.38 \times 10^{-23} \times 3 \times 10^{12}

=1.66×1010 V2= 1.66 \times 10^{-10} \text{ V}^2

Vrms=1.66×10101.29×105 V=12.9 \muVV_{\text{rms} = \sqrt{1.66 \times 10^{-10}} \approx 1.29 \times 10^{-5} \text{ V} = 12.9 \text{ \mu\text{V}}}

This sets a fundamental limit on the sensitivity of electrical measurements.

Worked Example 13.2: Brownian Motion and Einstein Relation

The Einstein relation is a special case of the FDT for Brownian motion. The diffusion constant DD relates to the mobility μ\mu:

D=μkBTD = \mu k_B T

For a spherical particle of radius rr in a fluid with viscosity η\eta:

μ=16πηr(Stokes drag)\mu = \frac{1}{6\pi\eta r} \quad \text{(Stokes drag)}

So D=kBT/(6πηr)D = k_B T/(6\pi\eta r).

For a 11 μ\muM diameter sphere in water (η=103\eta = 10^{-3} Pa\cdotS) at T=300T = 300 K:

D=1.38×1023×3006π×103×0.5×106=4.14×10219.42×1094.39×1013 m2/sD = \frac{1.38 \times 10^{-23} \times 300}{6\pi \times 10^{-3} \times 0.5 \times 10^{-6}} = \frac{4.14 \times 10^{-21}}{9.42 \times 10^{-9}} \approx 4.39 \times 10^{-13} \text{ m}^2/\text{s}

The mean squared displacement in time tt is x2=2Dt\langle x^2 \rangle = 2Dt. In 1 second: x20.94\sqrt{\langle x^2 \rangle} \approx 0.94 μ\muM.

QuantityExpressionPhysical Content
Response functionχAA(t)=iθ(t)[A(t),A(0)]0\chi_{AA}(t) = \frac{i}{\hbar}\theta(t)\langle[A(t), A(0)]\rangle_0Causal linear response
Classical FDTS(ω)=2kBTωχ(ω)S(\omega) = \frac{2k_B T}{\omega}\chi''(\omega)Fluctuations \leftrightarrow dissipation
Johnson—NyquistV2=4kBTRΔf\langle V^2 \rangle = 4k_B T R \Delta fVoltage noise in resistor
Einstein relationD=μkBTD = \mu k_B TDiffusion \leftrightarrow mobility
Nyquist formulaSV(ω)=2kBTR(ω)S_V(\omega) = 2k_B T R(\omega)Generalised impedance noise
  1. FDT assumes thermal equilibrium: The system must be in equilibrium at temperature TT. For out-of-equilibrium systems (driven, ageing, glasses), generalised fluctuation-dissipation relations apply with an effective temperature.
  2. Quantum corrections at low temperature: The classical FDT S(ω)=2kBTχ(ω)/ωS(\omega) = 2k_B T \chi''(\omega)/\omega is valid only for ωkBT\hbar\omega \ll k_B T. At low temperatures, the quantum FDT gives S(ω)=coth(ω/2kBT)χ(ω)S(\omega) = \hbar\coth(\hbar\omega/2k_B T)\,\chi''(\omega) which includes zero-point fluctuations.
  3. Response functions must be causal: The response function χAA(t)\chi_{AA}(t) must vanish for t<0t < 0 (causality). The Kramers—Kronig relations follow from this, connecting the real and imaginary parts of χ(ω)\chi(\omega).
  4. Linearity requirement: FDT is a result of linear response theory. For large perturbations, nonlinear effects break the simple relation between fluctuations and dissipation.
  • Electrical engineering: Johnson—Nyquist noise sets the fundamental noise floor in amplifiers, sensors, and communication systems. Cryogenic cooling reduces thermal noise for sensitive measurements.
  • Brownian motion: The Einstein relation enables determining Boltzmann’s constant via tracking colloidal particles, or measuring Avogadro’s number from diffusion measurements.
  • Optical trapping: Fluctuations of a trapped bead in optical tweezers obey the FDT, allowing calibration of trap stiffness from the power spectrum of position fluctuations.
  • Gravitational wave detectors: Thermal noise in mirror suspensions and test masses limits the sensitivity of LIGO at intermediate frequencies.
  • Biophysics: Single-molecule force spectroscopy (optical tweezers, AFM) uses fluctuation analysis to extract spring constants and dissipation in biomolecules.
  • Kramers—Kronig relations: Causality imposes integral relations between χ(ω)\chi'(\omega) and χ(ω)\chi''(\omega), which are intimately related to the FDT.
  • Onsager regression hypothesis: The relaxation of macroscopic nonequilibrium fluctuations follows the same laws as spontaneous fluctuations at equilibrium — a precursor to the FDT.
  • Nonequilibrium statistical mechanics: The fluctuation theorem (Evans—Searles, Crooks) extends fluctuation—dissipation ideas to far-from-equilibrium regimes, relating work distributions to free energy differences.
  • Quantum optics: The FDT applied to the electromagnetic field yields the Planck spectrum, connecting blackbody radiation to vacuum fluctuations and dissipation.
SystemFluctuationDissipationFDT Relation
ResistorVoltage noise V2\langle V^2 \rangleResistance RRV2=4kBTRΔf\langle V^2 \rangle = 4k_B T R \Delta f
Brownian particlePosition fluctuations x2\langle x^2 \rangleDrag coefficient γ\gammaD=kBT/γD = k_B T / \gamma
Harmonic oscillatorAmplitude fluctuationsDamping rate Γ\GammaSx(ω)=2kBTωχ(ω)S_x(\omega) = \frac{2k_B T}{\omega}\,\chi''(\omega)
Blackbody radiationField fluctuationsAbsorption cross sectionPlanck spectrum S(ω)=ω3π2c2(eω/kBT1)1S(\omega) = \frac{\hbar\omega^3}{\pi^2 c^2}(e^{\hbar\omega/k_BT} - 1)^{-1}
Magnetic systemMagnetisation noiseMagnetic susceptibilityχ(ω)=ω2kBTSM(ω)\chi''(\omega) = \frac{\omega}{2k_B T} S_M(\omega)

Additional Worked Example: Fluctuation-Dissipation in an RLC Circuit

Section titled “Additional Worked Example: Fluctuation-Dissipation in an RLC Circuit”

Problem. An RLC circuit at temperature TT has resistance RR, inductance LL, and capacitance CC. Find the power spectrum of voltage fluctuations across the capacitor using the FDT.

Solution. The impedance of the RLC circuit is Z(ω)=R+i(ωL1/(ωC))Z(\omega) = R + i(\omega L - 1/(\omega C)). The FDT in impedance form states: V2ω=2kBTRe[Z(ω)]\langle V^2 \rangle_\omega = 2k_B T \, \text{Re}[Z(\omega)]. For the RLC circuit:

VC2ω=2kBTR1/(iωC)R+i(ωL1/(ωC))2\langle V_C^2 \rangle_\omega = 2k_B T R \left| \frac{1/(i\omega C)}{R + i(\omega L - 1/(\omega C))} \right|^2

The total mean-square voltage across the capacitor is VC2=kBT/C\langle V_C^2 \rangle = k_B T / C (equipartition). This is independent of RR, illustrating that the fluctuation-dissipation relation always yields the correct thermal equilibrium result regardless of the dissipation mechanism.