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)=1kBTddt⟨A(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)⟩=kBTme−t/τ\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 f∼kBT/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−βW dW=∫PR(−W) e−βΔF dW=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=−kBTln⁡Z=−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⟩=p1⋅e0+p2⋅e−β(ϵ′−ϵ)=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−β(F′−F)=e−βF′eβF=(1+e−βϵ′)1Z=⟨e−βW⟩✓e^{-\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)=1kBTddt⟨A(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 equality⟨e−β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 f∼kBT/h≈6f \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 protocol⟨e−β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×10−23)(300)(104)≈1.66×10−16S_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 Δf≈1.66×10−16×106=1.66×10−10  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×10−10≈1.29×10−5  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 f≪fcf \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/4≈8.7×\sqrt{300/4} \approx 8.7\times.

■\blacksquare