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 F But 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)}
flowchart TD
A[19_Irreversible Thermodynamics And Fluctuations] --> B[Key Concepts]
A --> C[Core Principles]
A --> D[Practical Applications]
B --> E[Fundamental definitions]
C --> F[Design patterns]
D --> G[Real-world usage]

Irreversible thermodynamics describes systems that are driven away from equilibrium and eventually settle into a steady state. The Onsager relations reveal a deep symmetry: the way a temperature gradient drives particle flow is mathematically related to how a concentration gradient drives heat flow. Entropy production is always positive in irreversible processes, providing an arrow of time. Fluctuation theorems extend these ideas to small systems where thermal noise dominates, showing that entropy-decreasing fluctuations are possible but exponentially unlikely. These frameworks connect microscopic randomness to macroscopic irreversibility, explaining why heat flows from hot to cold and why perpetual motion machines are impossible.

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

This section provides detailed coverage of advanced concepts, including full derivations, proofs, and extended examples.

Complete mathematical derivations and proofs are provided where appropriate. Each step is explained to ensure understanding of the underlying reasoning.

Advanced examples demonstrate the application of concepts to complex problems. These examples go beyond standard exam questions to develop deeper understanding.

This material connects to current research and advanced applications in the field. Understanding these connections provides context for the study material.

Ensure you have mastered the prerequisite material before attempting this advanced content.