Skip to content

Irreversible Thermodynamics and Fluctuations

19.1 Fluctuation-Dissipation in the Time Domain

Section titled “19.1 Fluctuation-Dissipation in the Time Domain”

The classical fluctuation-dissipation theorem relates the autocorrelation function of a fluctuating variable to the linear response function:

χ(t)=1kBTddtA(t)A(0)\chi(t) = \frac{1}{k_BT}\frac{d}{dt}\langle A(t)A(0)\rangle

For example, the velocity autocorrelation function of a Brownian particle:

v(t)v(0)=kBTmet/τ\langle v(t)v(0)\rangle = \frac{k_BT}{m}e^{-t/\tau}

Gives the mobility μ=eτ/m\mu = e\tau/m (Einstein relation).

The voltage noise spectrum across a resistor RR at temperature TT:

SV(f)=4kBTRS_V(f) = 4k_BTR

This is white noise (frequency-independent up to fkBT/hf \sim k_BT/h).

The voltage fluctuation in bandwidth Δf\Delta f:

V2=4kBTRΔf\langle V^2 \rangle = 4k_BTR\,\Delta f

The Jarzynski equality (1997) connects non-equilibrium work to equilibrium free energy differences:

eβW=eβΔF\langle e^{-\beta W}\rangle = e^{-\beta\Delta F}

Where the average is over many realisations of a process that drives the system from equilibrium state AA to equilibrium state BB in time τ\tau.

Consequences:

  • By Jensen”s inequality: WΔF\langle W \rangle \geq \Delta F (the average work is never less than the free energy change).
  • For quasi-static processes: W=ΔF\langle W \rangle = \Delta F and the distribution of WW is a delta function.
  • For fast (far-from-equilibrium) processes: W>ΔF\langle W \rangle > \Delta FBut the exponential average still equals eβΔFe^{-\beta\Delta F}.

This remarkable result has been verified experimentally in single-molecule pulling experiments (RNA, DNA hairpins) using optical tweezers.

The Crooks theorem (1999) relates the work distributions for forward and reverse processes:

PF(W)PR(W)=eβ(WΔF)\frac{P_F(W)}{P_R(-W)} = e^{\beta(W - \Delta F)}

Where PF(W)P_F(W) is the probability distribution of work for the forward process and PR(W)P_R(W) for the reverse process.

This implies the Jarzynski equality as a special case:

PF(W)eβWdW=PR(W)eβΔFdW=eβΔF\int P_F(W)\,e^{-\beta W}\,dW = \int P_R(-W)\,e^{-\beta\Delta F}\,dW = e^{-\beta\Delta F}

Worked Example 19.1: Jarzynski Equality for a Two-Level System

Consider a two-level system with ϵ1=0\epsilon_1 = 0 and ϵ2=ϵ\epsilon_2 = \epsilonInitially in equilibrium at inverse temperature β\beta.

The free energy: F=kBTlnZ=kBTln(1+eβϵ)F = -k_BT\ln Z = -k_BT\ln(1 + e^{-\beta\epsilon}).

Now the energy gap is suddenly changed from ϵ\epsilon to ϵ"\epsilon". The work done is:

W={0withprob.p1=1/Zϵϵwithprob.p2=eβϵ/ZW = \begin{cases} 0 & \text{with} prob. p_1 = 1/Z \\ \epsilon' - \epsilon & \text{with} prob. p_2 = e^{-\beta\epsilon}/Z \end{cases}

The Jarzynski average:

eβW=p1e0+p2eβ(ϵϵ)=1Z+eβϵZ=1+eβϵZ\langle e^{-\beta W}\rangle = p_1 \cdot e^0 + p_2 \cdot e^{-\beta(\epsilon' - \epsilon)} = \frac{1}{Z} + \frac{e^{-\beta\epsilon'}}{Z} = \frac{1 + e^{-\beta\epsilon'}}{Z}

The new free energy: F=kBTln(1+eβϵ)F' = -k_BT\ln(1 + e^{-\beta\epsilon'}).

eβΔF=eβ(FF)=eβFeβF=(1+eβϵ)1Z=eβWe^{-\beta\Delta F} = e^{-\beta(F' - F)} = e^{-\beta F'}e^{\beta F} = (1 + e^{-\beta\epsilon'})\frac{1}{Z} = \langle e^{-\beta W}\rangle \quad \checkmark

The Jarzynski equality is verified exactly for this two-level system, even though the process is far from equilibrium (sudden quench).

TheoremStatementConnection to equilibrium
Fluctuation-dissip.χ(t)=1kBTddtA(t)A(0)\chi(t) = \frac{1}{k_BT}\frac{d}{dt}\langle A(t)A(0)\rangleResponse \leftrightarrow equilibrium fluctuations
Johnson-NyquistSV(f)=4kBTRS_V(f) = 4k_BTRVoltage noise \leftrightarrow resistance
Einstein relationD=μkBTD = \mu k_BTDiffusion \leftrightarrow mobility
Jarzynski equalityeβW=eβΔF\langle e^{-\beta W}\rangle = e^{-\beta\Delta F}Non-equilibrium work \leftrightarrow free energy
Crooks theoremPF(W)/PR(W)=eβ(WΔF)P_F(W)/P_R(-W) = e^{\beta(W - \Delta F)}Forward/reverse work distributions

The fluctuation-dissipation theorem unifies these: the Einstein relation and Johnson-Nyquist formula are special cases of the FDT applied to Brownian motion and electrical circuits, respectively.

  • Applying the FDT only to equilibrium systems. The standard FDT assumes the system is in thermal equilibrium. Fix: For non-equilibrium steady states, use generalised fluctuation-dissipation relations that include additional correlation terms.
  • Confusing white noise with infinite power. Johnson-Nyquist noise is white only up to fkBT/h6f \sim k_BT/h \approx 6 THz at 300 K; above this, quantum effects cut off the spectrum. Fix: Use SV(f)=4kBTR[f/(kBT)]/[exp(f/kBT)1]S_V(f) = 4k_BTR \cdot [\hbar f/(k_BT)]/[\exp(\hbar f/k_BT) - 1] for the quantum-corrected spectrum.
  • Assuming Jarzynski equality only applies to slow processes. The equality holds for arbitrarily fast (even instantaneous) processes. Fix: The work distribution for a fast process has large tails, but the exponential average still equals eβΔFe^{-\beta\Delta F} — verify with the two-level example.
  • Forgetting to take the exponential average in experiments. The average eβW\langle e^{-\beta W}\rangle is dominated by rare trajectories with negative work, requiring many samples to converge. Fix: Use Bennett’s acceptance ratio or Hummer-Szabo estimator for better convergence.
  • Single-molecule biophysics: Optical tweezers measure the work needed to unfold RNA/DNA hairpins; the Jarzynski equality extracts the folding free energy without requiring reversible pulling.
  • Nanoscale heat transfer: The FDT predicts thermal noise in nanomechanical resonators (cantilevers, membranes), limiting force sensitivity in AFM and gravitational-wave detectors.
  • Circuit design: Johnson-Nyquist noise sets the fundamental noise floor in amplifiers and receivers; cryogenic cooling reduces SVS_V linearly with TT.
  • Molecular dynamics: The Crooks theorem is used to compute free energy differences from non-equilibrium pulling simulations, avoiding expensive equilibrium sampling.
  • Brownian ratchets: Fluctuation theorems constrain the efficiency of molecular motors and information-driven devices (Maxwell’s demon).
ConceptTypeDomainKey formula
Fluctuation-dissipationGeneral theoremNear equilibriumχ(ω)=ω2kBTSA(ω)\chi''(\omega) = \frac{\omega}{2k_BT}S_A(\omega)
Johnson-Nyquist noiseSpecific applicationElectrical circuitsSV=4kBTRS_V = 4k_BTR
Einstein relationSpecific applicationBrownian motionD=μkBTD = \mu k_BT
Jarzynski equalityNon-equilibriumAny driving protocoleβW=eβΔF\langle e^{-\beta W}\rangle = e^{-\beta\Delta F}
Crooks theoremNon-equilibriumForward/reverse pairsPF(W)/PR(W)=eβ(WΔF)P_F(W)/P_R(-W) = e^{\beta(W - \Delta F)}

19.9 Worked Example: Johnson-Nyquist Noise in an RC Circuit

Section titled “19.9 Worked Example: Johnson-Nyquist Noise in an RC Circuit”

Problem. A 1010 kΩ\Omega resistor at T=300T = 300 K is connected to a 11 nF capacitor. Compute the RMS voltage fluctuation across the capacitor and the noise power in a 11 MHz bandwidth.

Solution. The voltage noise spectral density is SV(f)=4kBTR=4(1.38×1023)(300)(104)1.66×1016S_V(f) = 4k_BTR = 4(1.38\times10^{-23})(300)(10^4) \approx 1.66\times10^{-16} V2^2/Hz. In a Δf=1\Delta f = 1 MHz bandwidth:

V2=SVΔf1.66×1016×106=1.66×1010  V2\langle V^2\rangle = S_V\,\Delta f \approx 1.66\times10^{-16} \times 10^6 = 1.66\times10^{-10}\;\mathrm{V}^2

Vrms=1.66×10101.29×105  V=12.9  μVV_{\rm rms} = \sqrt{1.66\times10^{-10}} \approx 1.29\times10^{-5}\;\mathrm{V} = 12.9\;\mu\mathrm{V}

The RC low-pass filter (fc=1/(2πRC)16f_c = 1/(2\pi RC) \approx 16 kHz) limits the effective bandwidth if the capacitor is considered, but at ffcf \ll f_c the full 11 MHz bandwidth applies. This illustrates why sensitive electronics are cryogenically cooled: reducing TT from 300 K to 4 K reduces VrmsV_{\rm rms} by 300/48.7×\sqrt{300/4} \approx 8.7\times.

\blacksquare