Skip to content

Fluctuation Theorems

Fluctuation theorems are exact statements about the probability of thermodynamic fluctuations in small systems driven away from equilibrium. They explain how the second law emerges from reversible microscopic dynamics even though individual runs can temporarily move entropy in the “wrong” direction.

The unifying idea is a comparison between a forward experiment and a suitably defined reverse experiment. If γ\gamma denotes a trajectory or measurement record and γ~\tilde\gamma denotes its time-reversed counterpart, many fluctuation relations have the form

PF[γ]PR[γ~]=eσ[γ],\frac{P_F[\gamma]} {P_R[\tilde\gamma]} = e^{\sigma[\gamma]},

where σ[γ]\sigma[\gamma] is dimensionless entropy production. Averaging this ratio in the right way gives an integral fluctuation theorem,

⟨e−σ⟩F=1,\left\langle e^{-\sigma}\right\rangle_F = 1,

and Jensen’s inequality gives

⟨σ⟩F≥0.\langle\sigma\rangle_F \ge 0.

Thus the ordinary second-law inequality is a consequence of a stronger probability statement. The theorem does not say negative entropy-production events are impossible. It says their probabilities are constrained relative to the corresponding reverse events.

“Fluctuation theorem” is a family name, not a single formula. The right theorem depends on the variables and dynamics:

SettingTypical random variableStandard relation
closed driven systemwork WWJarzynski equality and Crooks relation
open thermal trajectoryentropy production σ\sigmadetailed and integral fluctuation theorems
nonequilibrium steady stateexchanged heat or entropy flowsteady-state fluctuation theorem
transport junctiontransferred charge, particles, or energycounting-statistics fluctuation symmetry
feedback-controlled processwork plus informationfeedback-modified fluctuation relation

The shared structure is microscopic reversibility plus a clear operational definition of the measured quantity. In quantum mechanics the word “measured” is not cosmetic: projective energy measurements, weak measurements, quantum jumps, homodyne records, and full counting statistics define different probability distributions.

A compact classical-looking statement uses dimensionless entropy production,

σ=ΔSsyskB−∑αβαQα,\sigma = \frac{\Delta S_{\mathrm{sys}}}{k_B} - \sum_\alpha \beta_\alpha Q_\alpha,

where QαQ_\alpha is heat entering the system from bath α\alpha. The bath entropy change is therefore −βαQα-\beta_\alpha Q_\alpha.

When a forward protocol and a reverse protocol are related by time reversal, a detailed fluctuation theorem has the schematic form

PF(σ)PR(−σ)=eσ.\frac{P_F(\sigma)} {P_R(-\sigma)} = e^\sigma.

Multiplying by PR(−σ)e−σP_R(-\sigma)e^{-\sigma} and integrating over σ\sigma gives

⟨e−σ⟩F=1.\left\langle e^{-\sigma}\right\rangle_F=1.

Jensen’s inequality then implies

⟨σ⟩F≥0.\langle\sigma\rangle_F\ge0.

This is the mathematical core of the relation between fluctuation theorems and the second law.

For a closed system driven from Hamiltonian H0H_0 to HτH_\tau, assume the initial state is Gibbs:

ρ0=e−βH0Z0.\rho_0 = \frac{e^{-\beta H_0}}{Z_0}.

The equilibrium free-energy difference is

ΔF=Fτ−F0=−1βln⁡ZτZ0.\Delta F = F_\tau-F_0 = - \frac{1}{\beta} \ln \frac{Z_\tau}{Z_0}.

In the standard two-point measurement scheme, one measures the initial energy, evolves with the driven unitary, and measures the final energy. For outcomes nn and mm,

Wmn=Emτ−En0.W_{mn} = E_m^\tau-E_n^0.

The forward work distribution is

PF(W)=∑n,mpn0 pF(m∣n) δ(W−Wmn),P_F(W) = \sum_{n,m} p_n^0\, p_F(m\mid n)\, \delta(W-W_{mn}),

where

pn0=e−βEn0Z0,pF(m∣n)=∣⟨mτ∣UF∣n0⟩∣2.p_n^0 = \frac{e^{-\beta E_n^0}}{Z_0}, \qquad p_F(m\mid n) = \left| \langle m^\tau|U_F|n^0\rangle \right|^2.

Under the usual microreversibility assumptions, the reverse protocol begins in the Gibbs state for HτH_\tau and runs the time-reversed drive. Crooks’ work relation is

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

Integrating over WW gives the Jarzynski equality:

⟨e−βW⟩F=e−βΔF.\left\langle e^{-\beta W}\right\rangle_F = e^{-\beta\Delta F}.

The equality is exact under its assumptions even for fast, irreversible driving. The average work bound follows from Jensen’s inequality:

⟨W⟩F≥ΔF.\langle W\rangle_F \ge \Delta F.

The two-point measurement scheme is standard because it gives a positive work distribution and a clean connection to energy conservation. It also has a cost: the first projective energy measurement removes initial coherences in the energy basis.

This is harmless when the initial state is Gibbs and therefore commutes with H0H_0. It becomes conceptually important when the initial state has energy coherence or when one wants to assign work without a strong initial measurement. Alternative quantum formulations include:

  • weak-measurement and quasiprobability approaches;
  • full counting statistics with counting fields;
  • interferometric measurements of characteristic functions;
  • quantum-jump and continuous-measurement trajectory definitions;
  • inclusive system-plus-bath energy measurements.

These formulations answer different operational questions. A theorem should always state which record defines the random variable.

For a weakly coupled open system described by a thermal quantum-jump unraveling, a trajectory is a sequence of no-jump evolutions and jumps. A jump associated with bath α\alpha and Bohr frequency ω\omega transfers energy

Qα=ℏωQ_\alpha = \hbar\omega

into the system when the jump raises the system energy, and Qα=−ℏωQ_\alpha=-\hbar\omega when it lowers the system energy. The sign convention must be stated.

Thermal rates satisfy local detailed balance. For a transition n→mn\to m caused by bath α\alpha,

Wm←n(α)Wn←m(α)=e−βα(Em−En)\frac{W_{m\leftarrow n}^{(\alpha)}} {W_{n\leftarrow m}^{(\alpha)}} = e^{-\beta_\alpha(E_m-E_n)}

when Em−EnE_m-E_n is the energy gained by the system. This microscopic ratio is the seed of the trajectory-level entropy production.

For a trajectory beginning in state nn and ending in state mm, one often writes

σ[γ]=ln⁡pn0pmτ−∑αβαQα[γ],\sigma[\gamma] = \ln \frac{p_n^0}{p_m^\tau} - \sum_\alpha \beta_\alpha Q_\alpha[\gamma],

where the first term is the stochastic system entropy change and the second is the bath entropy change. Under the right reverse process,

PF[γ]PR[γ~]=eσ[γ].\frac{P_F[\gamma]} {P_R[\tilde\gamma]} = e^{\sigma[\gamma]}.

This trajectory statement is the open-system analogue of the closed-system work relation. It relies on thermally consistent rates and a well-defined monitoring scheme. A master equation can be useful for state evolution while still failing to define a trustworthy thermodynamic trajectory ensemble if its jumps do not correspond to physical reservoir exchanges.

In mesoscopic transport and quantum optics, fluctuation relations are often written for generating functions. If P(N)P(N) is the probability of transferring NN particles, or P(Q)P(Q) the probability of transferring energy QQ, one introduces a counting field χ\chi and a generating function

Z(χ)=∑NP(N)eiχN.Z(\chi) = \sum_N P(N)e^{i\chi N}.

Thermodynamic consistency appears as a symmetry of Z(χ)Z(\chi) or of the long-time cumulant-generating function. For example, a two-reservoir particle current has an affinity built from chemical potentials and temperatures, and the fluctuation symmetry encodes the probability ratio between current flowing with and against the thermodynamic bias.

This language is especially useful when individual projective energy measurements on the system are not the most natural experimental record. It also connects fluctuation theorems to noise, response, and higher cumulants.

Fluctuation theorems are nonlinear nonequilibrium statements. Linear response and fluctuation-dissipation relations emerge when one expands near equilibrium.

For work, the Jarzynski equality gives a simple illustration. If WW is approximately Gaussian with mean μ\mu and variance vv, then

⟨e−βW⟩=e−βμ+β2v/2.\left\langle e^{-\beta W}\right\rangle = e^{-\beta\mu+\beta^2v/2}.

Jarzynski’s equality therefore implies

μ−ΔF=βv2.\mu-\Delta F = \frac{\beta v}{2}.

Near equilibrium, dissipated work is tied to work fluctuations. The full theorem is stronger than this Gaussian relation, but the relation shows why equilibrium fluctuations and irreversible response are deeply linked. See Fluctuation–Dissipation Relation for the linear-response side.

When a measurement outcome is used to choose the subsequent protocol, the reverse experiment and entropy budget must include information. A schematic feedback-modified Jarzynski relation is

⟨e−β(W−ΔF)−I⟩=1,\left\langle e^{-\beta(W-\Delta F)-I} \right\rangle = 1,

where II is an information term determined by the measurement record and feedback protocol. Jensen’s inequality then gives a bound of the form

⟨W⟩≥ΔF−kBT ⟨I⟩.\langle W\rangle \ge \Delta F - k_B T\,\langle I\rangle.

This does not mean information creates free energy from nowhere. It means the controller’s record is a thermodynamic resource. A complete cycle must include the physical memory, the measurement backaction, and the cost of resetting or maintaining the controller.

The open-system ingredients are developed in Quantum Instruments and Stochastic Master Equations.

Fluctuation relations are exact only after the process has been specified precisely. Common failure modes include:

  • using a work theorem for a heat distribution;
  • using the Jarzynski equality without an initial Gibbs state;
  • comparing a forward process to the wrong reverse process;
  • ignoring time-reversal parity of magnetic fields, spins, or driving parameters;
  • applying a two-point measurement relation to a state with important energy coherence without acknowledging the measurement backaction;
  • assigning heat to jumps in a phenomenological dissipator that is not thermally derived;
  • mixing local dissipators with strong internal interactions and then assuming local detailed balance;
  • treating feedback information as free without including memory and erasure.

The practical cure is to state the measured record, the initial ensemble, the reverse protocol, and the thermodynamic sign convention before writing the theorem.

Assume

PF(σ)PR(−σ)=eσ.\frac{P_F(\sigma)} {P_R(-\sigma)} = e^\sigma.

Derive ⟨e−σ⟩F=1\langle e^{-\sigma}\rangle_F=1.

Solution

Start from the forward average:

⟨e−σ⟩F=∫dσ PF(σ)e−σ.\left\langle e^{-\sigma}\right\rangle_F = \int d\sigma\, P_F(\sigma)e^{-\sigma}.

Using PF(σ)e−σ=PR(−σ)P_F(\sigma)e^{-\sigma}=P_R(-\sigma),

⟨e−σ⟩F=∫dσ PR(−σ).\left\langle e^{-\sigma}\right\rangle_F = \int d\sigma\,P_R(-\sigma).

Changing variables to σ′=−σ\sigma'=-\sigma gives

∫dσ′ PR(σ′)=1.\int d\sigma'\,P_R(\sigma')=1.

Use Jensen’s inequality to show that

⟨e−βW⟩=e−βΔF\left\langle e^{-\beta W}\right\rangle = e^{-\beta\Delta F}

implies ⟨W⟩≥ΔF\langle W\rangle\ge\Delta F.

Solution

For the convex function exe^x,

e⟨x⟩≤⟨ex⟩.e^{\langle x\rangle} \le \langle e^x\rangle.

Set x=−βWx=-\beta W. Then

e−β⟨W⟩≤e−βΔF.e^{-\beta\langle W\rangle} \le e^{-\beta\Delta F}.

Taking logarithms and multiplying by −1/β-1/\beta reverses the inequality because −1/β<0-1/\beta\lt0:

⟨W⟩≥ΔF.\langle W\rangle \ge \Delta F.

Consider a two-level Hamiltonian with energies 00 and ϵ0\epsilon_0 initially, changed suddenly to 00 and ϵτ\epsilon_\tau with the same eigenvectors. The initial state is Gibbs. Verify the Jarzynski equality.

Solution

There are no transitions because the eigenvectors are unchanged. The ground-state work is 00, and the excited-state work is ϵτ−ϵ0\epsilon_\tau-\epsilon_0.

The initial probabilities are

pg=11+e−βϵ0,pe=e−βϵ01+e−βϵ0.p_g = \frac{1}{1+e^{-\beta\epsilon_0}}, \qquad p_e = \frac{e^{-\beta\epsilon_0}} {1+e^{-\beta\epsilon_0}}.

Therefore

⟨e−βW⟩=pg+pe e−β(ϵτ−ϵ0).\left\langle e^{-\beta W}\right\rangle = p_g + p_e\,e^{-\beta(\epsilon_\tau-\epsilon_0)}.

Substituting pgp_g and pep_e gives

⟨e−βW⟩=1+e−βϵτ1+e−βϵ0=ZτZ0.\left\langle e^{-\beta W}\right\rangle = \frac{1+e^{-\beta\epsilon_\tau}} {1+e^{-\beta\epsilon_0}} = \frac{Z_\tau}{Z_0}.

Since

e−βΔF=ZτZ0,e^{-\beta\Delta F} = \frac{Z_\tau}{Z_0},

the Jarzynski equality holds.

Assume the feedback-modified relation

⟨e−β(W−ΔF)−I⟩=1.\left\langle e^{-\beta(W-\Delta F)-I} \right\rangle = 1.

Use Jensen’s inequality to derive the corresponding average-work bound.

Solution

Let

x=−β(W−ΔF)−I.x = -\beta(W-\Delta F)-I.

Jensen’s inequality gives

e⟨x⟩≤⟨ex⟩=1.e^{\langle x\rangle} \le \left\langle e^x\right\rangle = 1.

Thus ⟨x⟩≤0\langle x\rangle\le0, or

−β(⟨W⟩−ΔF)−⟨I⟩≤0.-\beta(\langle W\rangle-\Delta F) - \langle I\rangle \le 0.

Rearranging,

⟨W⟩≥ΔF−kBT ⟨I⟩.\langle W\rangle \ge \Delta F - k_B T\,\langle I\rangle.
  • C. Jarzynski, “Nonequilibrium equality for free energy differences,” Physical Review Letters 78, 2690-2693 (1997).
  • G. E. Crooks, “Entropy production fluctuation theorem and the nonequilibrium work relation for free energy differences,” Physical Review E 60, 2721-2726 (1999).
  • J. Kurchan, “A quantum fluctuation theorem,” arXiv:cond-mat/0007360 (2000).
  • H. Tasaki, “Jarzynski relations for quantum systems and some applications,” arXiv:cond-mat/0009244 (2000).
  • P. Talkner, E. Lutz, and P. Hanggi, “Fluctuation theorems: Work is not an observable,” Physical Review E 75, 050102(R) (2007).
  • M. Esposito, U. Harbola, and S. Mukamel, “Nonequilibrium fluctuations, fluctuation theorems, and counting statistics in quantum systems,” Reviews of Modern Physics 81, 1665-1702 (2009).
  • M. Campisi, P. Hanggi, and P. Talkner, “Colloquium: Quantum fluctuation relations: Foundations and applications,” Reviews of Modern Physics 83, 771-791 (2011).
  • T. Sagawa and M. Ueda, “Generalized Jarzynski equality under nonequilibrium feedback control,” Physical Review Letters 104, 090602 (2010).
  • U. Seifert, “Stochastic thermodynamics, fluctuation theorems and molecular machines,” Reports on Progress in Physics 75, 126001 (2012).