Jarzynski equality

Nonequilibrium work identity: ⟨e^{-βW}⟩ = e^{-βΔF} for protocols starting in equilibrium, yielding ΔF from nonequilibrium experiments and implying ⟨W⟩≥ΔF.
Jarzynski equality

Statement

Let a system at inverse temperature β=1/(kBT)\beta=1/(k_B T) (with TT the ) start in equilibrium at control parameter λ0\lambda_0, typically the .

Drive the system with an arbitrary (possibly fast) protocol λt\lambda_t over t[0,τ]t\in[0,\tau]. For each realization, measure the work WW performed on the system (see ). Let

ΔF  =  F(λτ)F(λ0) \Delta F \;=\; F(\lambda_\tau)-F(\lambda_0)

be the equilibrium free-energy difference at the same β\beta (see ).

Then the Jarzynski equality is

eβW  =  eβΔF, \left\langle e^{-\beta W}\right\rangle \;=\; e^{-\beta \Delta F},

where \langle\cdot\rangle denotes the ensemble average over repeated realizations (an with respect to the induced work distribution).

Conditions (minimal checklist)

  • Initial state is equilibrium at λ0\lambda_0 (canonical, or appropriate equilibrium ensemble).
  • Dynamics are microreversible (time-reversal symmetric Hamiltonian dynamics, or suitable stochastic dynamics coupled to a heat bath).
  • Work is defined consistently with the driven parameter (protocol work).

These assumptions are also those used for the , from which Jarzynski can be derived.

Consequences

Second-law bound (via Jensen)

Since exp\exp is convex,

eβW    eβW  =  eβΔF, e^{-\beta\langle W\rangle} \;\le\; \left\langle e^{-\beta W}\right\rangle \;=\; e^{-\beta\Delta F},

hence

W    ΔF, \langle W\rangle \;\ge\; \Delta F,

consistent with the .

Dissipated work

Define the dissipated work

Wdiss  =  WΔF. W_{\mathrm{diss}} \;=\; W-\Delta F.

Jarzynski becomes eβWdiss=1\langle e^{-\beta W_{\mathrm{diss}}}\rangle = 1, emphasizing that fluctuations with unusually small WdissW_{\mathrm{diss}} control the exponential average.

Near-equilibrium expansion

If WW has small fluctuations around its mean, cumulant expansion gives

ΔFWβ2Var(W)+, \Delta F \approx \langle W\rangle - \frac{\beta}{2}\,\mathrm{Var}(W) + \cdots,

linking nonequilibrium work fluctuations to equilibrium free-energy differences (compare with linear-response identities such as ).

Practical note (convergence)

Estimating ΔF\Delta F via β1lneβW-\beta^{-1}\ln\langle e^{-\beta W}\rangle can be statistically challenging because rare low-work events dominate the exponential average. Using both forward and reverse processes with often improves efficiency.