Thermodynamic Uncertainty Relations

We take the entropy production of the previous chapters — the Schnakenberg current-times-affinity sum that entropy production makes strictly positive out of equilibrium — and ask a sharper question than "is it positive?" We ask what that dissipation buys. The thermodynamic uncertainty relation (TUR) answers with a hard trade-off: the precision of any fluctuating nonequilibrium current is bounded by the entropy produced to sustain it. Extend the fluctuation–dissipation theorem past linear response, stand it next to the fluctuation theorems of Jarzynski and Crooks, and the TUR is the inequality that emerges — precision is never free. Take a system relaxed into a nonequilibrium steady state: a Markov jump process, an overdamped colloid in a force field, a molecular motor stepping along a track. Any current — particle flux, heat flow, motor cycles — integrated over a time tt has a mean ⟨Jt⟩\langle J_t\rangle and a variance Var(Jt)\mathrm{Var}(J_t). Barato and Seifert conjectured in 2015, and Gingrich, Horowitz, Perunov and England proved in 2016, that Var(Jt) Σt≥2kB⟨Jt⟩2\mathrm{Var}(J_t)\,\Sigma_t\ge2k_B\langle J_t\rangle^2, where Σt\Sigma_t is the total entropy produced. The relation is astonishing in its reach: it holds for every current, in every steady state, whatever the network or rates, and saturates only in a near-equilibrium Gaussian limit. Two proofs illuminate it. A large-deviation argument bounds the current's rate function by a universal parabola set by the entropy production; a Cramér–Rao argument treats the current as an estimator and caps the Fisher information by the dissipation. Both routes reveal the same content — a signal-to-noise ratio ceiling. And because only the mean and variance enter, the TUR doubles as a thermodynamic thermometer: measure a current's fluctuations and read off a model-free lower bound on how much entropy the hidden machinery must be burning.

Key takeaways
  • The thermodynamic uncertainty relation bounds the squared relative fluctuation of any nonequilibrium current by the entropy production: Var(J)/⟨J⟩2≥2kB/Σ\mathrm{Var}(J)/\langle J\rangle^2\ge 2k_B/\Sigma.
  • Equivalently, the squared signal-to-noise ratio is capped by half the dissipation: (⟨J⟩/σJ)2≤Σ/2kB(\langle J\rangle/\sigma_J)^2\le\Sigma/2k_B. Precision costs dissipation.
  • It is universal and model-free — true for every current in every continuous-time Markov / overdamped-Langevin NESS, without knowing the rates.
  • It extends the fluctuation–dissipation theorem far from equilibrium: near equilibrium it saturates and reduces to the FDT; far away it has slack but never reverses.
  • Two proofs: a large-deviation parabolic bound on the rate function, and a Cramér–Rao inequality bounding Fisher information by entropy production.
  • As inference: Σ≥2kB⟨J⟩2/Var(J)\Sigma\ge2k_B\langle J\rangle^2/\mathrm{Var}(J) turns measured fluctuations into a guaranteed floor on dissipation — the basis of thermodynamic inference for molecular motors and clocks.
  • Key idea: in a nonequilibrium steady state the squared signal-to-noise ratio of any current is bounded by half the total entropy production in units of kBk_B; a precise clock or motor must burn free energy at a minimum rate set by the TUR.
The core insight

In a nonequilibrium steady state, the squared signal-to-noise ratio of any current is bounded by half the total entropy production in units of kBk_B: (⟨J⟩/σJ)2≤Σ/(2kB)\big(\langle J\rangle/\sigma_J\big)^2\le\Sigma/(2k_B). Precision costs dissipation, and the exchange rate is exactly 2kB2k_B.

Where this lives today

The TUR has become a working instrument, not just a theorem. Its inference form, σ^min⁡=2kB⟨J⟩2/(Var(J) t)\hat\sigma_{\min}=2k_B\langle J\rangle^2/(\mathrm{Var}(J)\,t), lets experimentalists put a guaranteed floor on the entropy production of a system whose microscopic rates are unknown — a stalled kinesin, an active colloid, a beating flagellum — from trajectory data alone. This "thermodynamic inference" is now a standard estimator in single-molecule biophysics and active matter. In simulation, checking ⟨J⟩\langle J\rangle, Var(J)\mathrm{Var}(J) and an independent Σ\Sigma is a consistency test: an apparent violation flags a sampling error, hidden non-Markovian memory, or a discrete-time update that falls outside the standard TUR class. And in the design of nanoscale heat engines and sensors, the TUR draws a hard boundary on the constancy axis: maximizing power at fixed reliability forces a minimum dissipation, so the achievable efficiency-at-given-precision is bounded. A steady stream of generalizations — finite-time, discrete-time, periodically driven, multidimensional, and quantum — keeps extending the reach of the same core idea.

The physical picture

Imagine a molecular motor stepping forward and back along a track — a noisy random walk with a bias. To make its progress reliable (small fractional spread in position after many steps) you must bias each step strongly, and a strong bias means burning more fuel per step. The TUR makes this folk wisdom exact: to halve the relative noise you must quadruple the minimum dissipation. Entropy production is the budget a steady state spends to swim against thermal randomness, and a sharp directed signal cannot be extracted from omnidirectional noise for free. Near equilibrium the bound is tight — there the current is linear-response-Gaussian and the TUR is the fluctuation–dissipation theorem. Push far from equilibrium and the inequality develops slack: you dissipate more than the minimum the bound would allow, but you never beat it. The gap between what you spend and the TUR floor is a clean measure of how far from equilibrium — and how inefficiently precise — a machine runs.

Figure
The thermodynamic uncertainty relation carves out a forbidden region: precision can never exceed half the entropy production, with a near-equilibrium point sitting on the boundary and a driven point strictly inside
Precision (⟨J⟩/σ_J)² vs dissipation Σ/k_B — the TUR ceiling Σ / k_B (⟨J⟩/σ_J)² boundary: precision = Σ/2k_B forbidden by the TUR allowed near-equilibrium ring: ON the line (saturation = FDT) driven far from equilibrium: strictly inside (slack)

Every steady-state current lands on or below the line precision =Σ/2kB=\Sigma/2k_B. The blue point is an overdamped driven ring, which saturates the bound exactly; the green point is a strongly driven network, which dissipates more than the minimum and sits inside the allowed wedge.

The information-geometric picture sharpens this. A tiny virtual perturbation that speeds up the dynamics changes the current's statistics; how distinguishable that change is — the Fisher information — is capped by the entropy production. Cramér–Rao then says any current used as an estimator has variance at least the inverse Fisher information. That is the TUR again, now read as a statement about how much a steady state can "know" about its own driving.

Figure
A precise clock must dissipate: the number of coherent cycles a molecular motor can complete is bounded by the free energy it burns per cycle, and the bound is tight only near equilibrium
TUR quality factor Q = Var(J)·Σ / (2k_B⟨J⟩²) ≥ 1 affinity A = ln(k₊/k₋) (drive) Q Q = 1: TUR saturated (tightest) equilibrium (A→0): Q → 1 far from equilibrium: Q exceeds 1, extra dissipation biased hopper: Q(A) = A·(x+1)/[2(x−1)], x = e^A

The quality factor Q=Var(J)Σ/(2kB⟨J⟩2)Q=\mathrm{Var}(J)\Sigma/(2k_B\langle J\rangle^2) is the TUR written as Q≥1Q\ge1. For a biased hopper it equals 11 only as the drive A→0A\to0 (equilibrium, where the TUR meets the FDT) and climbs above 11 as the drive grows — the machine pays more dissipation than the precision strictly demands.

The core equations

The central statement is a single inequality tying three measured quantities together. Read it as: to shrink the relative noise of a current you must spend entropy, at an exchange rate of exactly 2kB2k_B.

Var(Jt)⟨Jt⟩2 Σt  ≥  2kB⟺Var(Jt) Σt  ≥  2kB ⟨Jt⟩2(1)\frac{\mathrm{Var}(J_t)}{\langle J_t\rangle^2}\,\Sigma_t\;\ge\;2k_B\qquad\Longleftrightarrow\qquad \mathrm{Var}(J_t)\,\Sigma_t\;\ge\;2k_B\,\langle J_t\rangle^2 \qquad (1)

where: JtJ_t — time-integrated current (net particles, cycles, or charge counted up to tt); ⟨Jt⟩\langle J_t\rangle — steady-state mean; Var(Jt)\mathrm{Var}(J_t) — its variance; Σt\Sigma_t — total entropy produced (system + reservoirs) over tt, in energy/temperature units; kBk_B — Boltzmann constant. Reading it physically: the left side is (relative noise)−1^{-1} paid for by dissipation; equality is reachable only in a Gaussian near-equilibrium limit. The entropy production rate is the Schnakenberg current–affinity sum — a result entropy production (the previous article) derives and proves nonnegative; here we recall it by name and use it.

Σt=σ t,σ=kB∑i<j(pjsswij−pisswji)ln⁡pjsswijpisswji  ≥  0(2)\Sigma_t=\sigma\,t,\qquad \sigma=k_B\sum_{i<j}\big(p_j^{ss}w_{ij}-p_i^{ss}w_{ji}\big)\ln\frac{p_j^{ss}w_{ij}}{p_i^{ss}w_{ji}}\;\ge\;0 \qquad (2)

where: σ\sigma — entropy production rate; wijw_{ij} — transition rate from state jj to state ii; pissp_i^{ss} — stationary probability of state ii. Each term is a current times its affinity; the product (a−b)ln⁡(a/b)≥0(a-b)\ln(a/b)\ge0 makes σ≥0\sigma\ge0 — the second law. It vanishes iff every net current vanishes (detailed balance). At long times any current becomes drift-plus-diffusion, which is what lets us define a precision rate.

⟨Jt⟩≃v t,Var(Jt)≃2D t(t→∞)(3)\langle J_t\rangle\simeq v\,t,\qquad \mathrm{Var}(J_t)\simeq2D\,t\quad(t\to\infty) \qquad (3)

where: vv — mean current rate (drift); DD — diffusivity of the integrated current (half its variance rate). With these, the TUR sheds its explicit tt-dependence and becomes a statement about rates.

v2D  ≤  σkB⟺σ  ≥  kBv2D(4)\frac{v^2}{D}\;\le\;\frac{\sigma}{k_B}\qquad\Longleftrightarrow\qquad \sigma\;\ge\;\frac{k_B v^2}{D} \qquad (4)

where: v2/Dv^2/D — the precision rate, large when a current is fast and quiet; σ/kB\sigma/k_B — dissipation rate in natural units. This is the form to reach for: it says the "sharpness" of a current is capped by the entropy rate. Rearranged, σ^min⁡=kBv2/D\hat\sigma_{\min}=k_B v^2/D is a model-free lower bound on dissipation from measured vv and DD alone. Finally, the large-deviation route needs the current's rate function and the universal parabola that bounds it.

I(j)=sup⁡s[sj−λ(s)],I(j)  ≤  (j−v)24v2 σkB(5)I(j)=\sup_s\big[sj-\lambda(s)\big],\qquad I(j)\;\le\;\frac{(j-v)^2}{4v^2}\,\frac{\sigma}{k_B} \qquad (5)

where: λ(s)\lambda(s) — scaled cumulant generating function of the current (a construct large deviation theory owns and we recall); I(j)I(j) — rate function, with P(Jt/t≈j)∼e−tI(j)P(J_t/t\approx j)\sim e^{-tI(j)}. The parabola passes through I(v)=0I(v)=0 and carries the entropy production as its curvature scale — the engine of the Gingrich–Horowitz–Perunov–England proof.

Derivation AThe TUR from the large-deviation parabolic bound

This is the original rigorous route: bound the current's large-deviation rate function by a universal parabola whose curvature is set by the entropy production, then read the TUR off the curvature at the minimum. We work in a Markov-jump nonequilibrium steady state and recall from large deviation theory that the scaled cumulant generating function λ(s)\lambda(s) encodes all current statistics, with λ′(0)=v\lambda'(0)=v and λ′′(0)=2D\lambda''(0)=2D. The one subtle number is the curvature of the Gaussian rate function at its minimum — get it right and the factor of 2kB2k_B falls out cleanly.

Step 1. Recall the Gaussian core of any current's rate function. From ⟨Jt⟩=vt\langle J_t\rangle=vt and Var(Jt)=2Dt\mathrm{Var}(J_t)=2Dt, the central-limit form is P(Jt/t≈j)∼e−t(j−v)2/4DP(J_t/t\approx j)\sim e^{-t(j-v)^2/4D}, so

I(j)=(j−v)24D,I′′(v)=12D.I(j)=\frac{(j-v)^2}{4D},\qquad I''(v)=\frac{1}{2D}.

Why: Jt/tJ_t/t has mean vv and variance 2D/t2D/t; the Gaussian exponent is (j−v)2/(2⋅2D/t)=t(j−v)2/4D(j-v)^2/(2\cdot2D/t)=t(j-v)^2/4D. Differentiating twice gives I′′(v)=1/(2D)I''(v)=1/(2D), not 1/(4D)1/(4D) — a slip the dry base carried in this step (see the note below).

Step 2. Invoke the Gingrich–Horowitz–Perunov–England universal parabola. Tilting the jump rates to scale the current while accounting for the entropy cost bounds the rate function above by

I(j)≤P(j)=(j−v)24v2 σkB,P′′(v)=σ2kBv2.I(j)\le P(j)=\frac{(j-v)^2}{4v^2}\,\frac{\sigma}{k_B},\qquad P''(v)=\frac{\sigma}{2k_B v^2}.

Why: the KL divergence between the tilted and original path measures is quadratic in the tilt and controlled by Σt\Sigma_t; optimizing the simplest tilt yields this parabola through I(v)=0I(v)=0.

Step 3. Compare curvatures at the shared minimum j=vj=v. Since I(j)≤P(j)I(j)\le P(j) with both zero and minimized at vv,

I′′(v)≤P′′(v)⟹12D≤σ2kBv2.I''(v)\le P''(v)\quad\Longrightarrow\quad \frac{1}{2D}\le\frac{\sigma}{2k_B v^2}.

Why: two nonnegative functions sharing a zero-minimum must obey the same order in their curvatures there.

Step 4. Rearrange to the rate-form TUR and restore time:

v2D≤σkB→ ⟨Jt⟩=vt, Var=2Dt, Σt=σt Var(Jt) Σt≥2kB⟨Jt⟩2.\frac{v^2}{D}\le\frac{\sigma}{k_B}\quad\xrightarrow{\ \langle J_t\rangle=vt,\ \mathrm{Var}=2Dt,\ \Sigma_t=\sigma t\ }\quad \mathrm{Var}(J_t)\,\Sigma_t\ge2k_B\langle J_t\rangle^2.

Why: multiplying v2/D≤σ/kBv^2/D\le\sigma/k_B by 2Dt22Dt^2 and substituting the long-time forms reproduces equation (1) exactly. Verified: `verify_aq-sm-thermodynamic-uncertainty-relations.py` confirms I(j)=(j−v)2/4D⇒I′′(v)=1/2DI(j)=(j-v)^2/4D\Rightarrow I''(v)=1/2D, the parabola curvature σ/2kBv2\sigma/2k_Bv^2, and the implication I′′≤P′′⇒v2/D≤σ/kBI''\le P''\Rightarrow v^2/D\le\sigma/k_B. Base-error note: the dry base's Step 6 wrote I′′(v)=1/(4D)I''(v)=1/(4D); the correct value is 1/(2D)1/(2D) (the final TUR is unchanged because the base's last step used the correct v2/D≤σ/kBv^2/D\le\sigma/k_B).

Go deeper · Derivation B · 1 minExact saturation for the overdamped driven ring
Derivation BExact saturation for the overdamped driven ring

Here is the structural result this article owns: a case where the TUR is not merely an inequality but an exact equality for all times. A single overdamped colloid driven around a ring has Gaussian current statistics and purely dissipative entropy production, and the two conspire — through the Einstein relation — to hit the bound identically. This is the cleanest demonstration that 2kB2k_B is the right constant, and it is the near-equilibrium limit where the TUR reduces to the fluctuation–dissipation theorem.

Step 1. Model and current. A particle on a ring, mobility μ\mu, diffusion D0=μkBTD_0=\mu k_BT (Einstein relation), driven by a constant tangential force ff: x˙=μf+2D0 ξ(t)\dot x=\mu f+\sqrt{2D_0}\,\xi(t). The integrated current is the displacement,

⟨Jt⟩=μf t,Var(Jt)=2D0 t.\langle J_t\rangle=\mu f\,t,\qquad \mathrm{Var}(J_t)=2D_0\,t.

Why: free drift–diffusion gives an exactly Gaussian displacement with drift μf\mu f and diffusion D0D_0.

Step 2. Entropy production. The force does work at rate f⋅⟨x˙⟩=μf2f\cdot\langle\dot x\rangle=\mu f^2, fully dissipated as heat at temperature TT:

σ=μf2T,Σt=μf2T t.\sigma=\frac{\mu f^2}{T},\qquad \Sigma_t=\frac{\mu f^2}{T}\,t.

Why: in a steady state the system's own entropy is constant, so all the injected work is exported as heat; divide by TT for the entropy rate.

Step 3. Assemble the product and use the Einstein relation:

Var(Jt) Σt=(2D0t)(μf2Tt)=2μkBT⋅μf2T t2=2kB(μft)2.\mathrm{Var}(J_t)\,\Sigma_t=(2D_0 t)\Big(\tfrac{\mu f^2}{T}t\Big)=2\mu k_BT\cdot\frac{\mu f^2}{T}\,t^2=2k_B(\mu f t)^2.

Why: substituting D0=μkBTD_0=\mu k_BT cancels the temperature and leaves a perfect square in the mean current.

Step 4. Recognize the bound, exactly:

Var(Jt) Σt=2kB⟨Jt⟩2for all t.\mathrm{Var}(J_t)\,\Sigma_t=2k_B\langle J_t\rangle^2\quad\text{for all }t.

Why: ⟨Jt⟩=μft\langle J_t\rangle=\mu f t, so the right side is 2kB(μft)22k_B(\mu f t)^2 — the TUR is saturated identically, not just asymptotically. Verified: `verify_aq-sm-thermodynamic-uncertainty-relations.py` proves Var(Jt)Σt−2kB⟨Jt⟩2=0\mathrm{Var}(J_t)\Sigma_t-2k_B\langle J_t\rangle^2=0 symbolically using D0=μkBTD_0=\mu k_BT, and that the rate form v2/D0=σ/kBv^2/D_0=\sigma/k_B holds with equality.

Go deeper · Derivation C · 1 minThe Cramér–Rao route: precision bounded by Fisher information
Derivation CThe Cramér–Rao route: precision bounded by Fisher information

The second proof recasts the TUR as an instance of the oldest inequality in statistics — Cramér–Rao — and explains why the bound feels like an uncertainty principle. The idea is to treat a current as an estimator of a fictitious perturbation parameter, then bound its variance below by the inverse Fisher information of the trajectory ensemble, and finally cap that Fisher information by the entropy production. This information-geometric picture (Hasegawa–Van Vu; Dechant–Sasa) unifies the TUR with speed limits and connects it to the Fisher information that also governs quantum metrology.

Step 1. Perturb the dynamics by a parameter θ\theta that uniformly rescales time (speeds up the process). For any current viewed as an estimator, Cramér–Rao gives

Var(J)≥(∂θ⟨J⟩)2I(θ).\mathrm{Var}(J)\ge\frac{(\partial_\theta\langle J\rangle)^2}{\mathcal I(\theta)}.

Why: Cramér–Rao lower-bounds the variance of any unbiased estimator by the inverse Fisher information I(θ)\mathcal I(\theta) of the underlying distribution.

Step 2. Choose the time-rescaling so that ∂θ⟨J⟩=⟨J⟩\partial_\theta\langle J\rangle=\langle J\rangle:

Var(J)≥⟨J⟩2I(θ).\mathrm{Var}(J)\ge\frac{\langle J\rangle^2}{\mathcal I(\theta)}.

Why: a uniform speed-up multiplies every integrated current by the same factor, so its derivative equals the current itself.

Step 3. Cap the path Fisher information by the entropy production,

I(θ)≤ΣkB⟹Var(J)≥kB⟨J⟩2Σ,\mathcal I(\theta)\le\frac{\Sigma}{k_B}\quad\Longrightarrow\quad \mathrm{Var}(J)\ge\frac{k_B\langle J\rangle^2}{\Sigma},

which rearranges to Var(J)Σ≥kB⟨J⟩2\mathrm{Var}(J)\Sigma\ge k_B\langle J\rangle^2, and with the tighter time-rescaling estimate to the factor 2kB2k_B of equation (1). Why: the Fisher information of the time-rescaling perturbation is exactly (or bounded by) the dissipation — distinguishing "fast" from "slow" dynamics costs entropy. Verified: the Cramér–Rao structure and the biased-hopper example (which satisfies v2/D≤σ/kBv^2/D\le\sigma/k_B with slack, saturating near equilibrium) are checked in `verify_aq-sm-thermodynamic-uncertainty-relations.py`.

Limiting cases

CaseResultMeaning
Equilibrium⟨J⟩=0, Σ=0\langle J\rangle=0,\ \Sigma=0Detailed balance: every current mean and the dissipation vanish; the TUR is the trivially true 0≥00\ge0.
Near-equilibrium / linear responseVar(J) Σ→2kB⟨J⟩2\mathrm{Var}(J)\,\Sigma\to2k_B\langle J\rangle^2Gaussian fluctuations and Onsager symmetry saturate the bound — the tight limit, recovering the FDT.
Far from equilibriumVar(J) Σ>2kB⟨J⟩2\mathrm{Var}(J)\,\Sigma>2k_B\langle J\rangle^2Strict inequality with growing slack; the bound holds but is no longer tight.
Finite timeΣt Var(Jt)≥2kB⟨Jt⟩2\Sigma_t\,\mathrm{Var}(J_t)\ge2k_B\langle J_t\rangle^2Valid for all tt, not only asymptotically (Horowitz–Gingrich), with Σt\Sigma_t the entropy up to tt.
Periodic / time-dependent drivingstandard TUR can be violatedA time-homogeneous NESS is required; explicit driving needs modified bounds (Barato–Chetrite, Koyuk–Seifert).

Worked examples

Example 1 — Reading dissipation off a trajectory. A colloid on a force-driven ring is watched for a time tt, giving mean displacement ⟨Jt⟩=10\langle J_t\rangle=10 (length units) and variance Var(Jt)=5\mathrm{Var}(J_t)=5. What is the least entropy it can have produced?

Step 1. Rearrange the TUR to bound Σt\Sigma_t from below:

Σt≥2kB⟨Jt⟩2Var(Jt)=2kB(10)25=40 kB.\Sigma_t\ge\frac{2k_B\langle J_t\rangle^2}{\mathrm{Var}(J_t)}=\frac{2k_B(10)^2}{5}=\mathbf{40\,k_B}.

Step 2. Interpret. The system must have produced at least 40 kB40\,k_B of entropy to reach this precision. If it is in the Gaussian saturating regime, 40 kB40\,k_B is also the actual value; otherwise the true Σt\Sigma_t is larger — the measured fluctuations set a floor, never a ceiling.

Step 3. This is thermodynamic inference in one line: no rates, no model, just ⟨J⟩\langle J\rangle and Var(J)\mathrm{Var}(J). The bound is model-free.

Example 2 — A driven hopper never beats the bound. A particle hops right at rate k+k_+ and left at rate k−k_- on a lattice, with x≡k+/k−x\equiv k_+/k_- the drive. Check the rate-form TUR and where it saturates.

Step 1. Compound-Poisson statistics give v=k+−k−v=k_+-k_-, variance rate 2D=k++k−2D=k_++k_-, and (recalling the Schnakenberg form) dissipation rate σ/kB=(k+−k−)ln⁡(k+/k−)\sigma/k_B=(k_+-k_-)\ln(k_+/k_-).

Step 2. The rate-form TUR v2/D≤σ/kBv^2/D\le\sigma/k_B reduces to the pure inequality

2(x−1)x+1≤ln⁡x,\frac{2(x-1)}{x+1}\le\ln x,

whose difference has derivative (x−1)2/[x(x+1)2]≥0(x-1)^2/[x(x+1)^2]\ge0 — so it holds for all x>1x>1, with equality only at x=1x=1 (equilibrium).

Step 3. Numbers: at x=2x=2 the quality factor is Q=1.040Q=1.040 (4% slack); at x=20x=20, Q=1.66Q=1.66 (66% slack). Far from equilibrium the hopper dissipates far more than the minimum, while at x→1x\to1 it saturates and the TUR becomes the FDT.

Physics of the portrait

Origins · 1 minmolecular motors and how this was discovered
MM
molecular motors
the arena of the TUR, 1990s–present · Where precision is paid for in ATP

Molecular motors — kinesin and dynein walking cytoskeletal tracks, the rotary F1-ATPase, the bacterial flagellar motor — are the biological machines that dramatized the TUR. Each takes discrete, noisy steps fueled by ATP hydrolysis, and each must be reliable despite ferocious thermal buffeting. The TUR quantifies the price of that reliability: the number of coherent steps or oscillations a motor can complete is bounded by the free energy it dissipates per cycle, roughly N≲Δμ/2kBT\mathcal N\lesssim\Delta\mu/2k_BT, where Δμ≈20 kBT\Delta\mu\approx20\,k_BT is the chemical-potential drop of one ATP. André Barato and Udo Seifert conjectured the TUR in 2015 while analyzing exactly such biochemical clocks; Todd Gingrich, Jordan Horowitz, Nikolay Perunov and Jeremy England proved it in 2016. The relation reframed molecular-motor biophysics as a three-way trade-off among power, efficiency, and precision — with dissipation as the currency that buys the last one.

History
2015
— Barato and Seifert conjecture the thermodynamic uncertainty relation for biomolecular processes
2016
— Gingrich, Horowitz, Perunov and England prove the TUR from a large-deviation bound
2016
— Pietzonka, Barato and Seifert extend the bound to finite times and affinities
2019
— Hasegawa and Van Vu derive the TUR from the Cramér–Rao inequality
2020
— Horowitz and Gingrich review thermodynamic uncertainty relations across the field

Going further

References. Barato & Seifert, Phys. Rev. Lett. 114, 158101 (2015) · Gingrich, Horowitz, Perunov & England, Phys. Rev. Lett. 116, 120601 (2016) · Pietzonka, Barato & Seifert, Phys. Rev. E 93, 052145 (2016) · Horowitz & Gingrich, Nature Physics 16, 15 (2020) — review · Seifert, Rep. Prog. Phys. 75, 126001 (2012) — stochastic thermodynamics · Hasegawa & Van Vu, 'Uncertainty relations in stochastic processes: an information inequality approach', Phys. Rev. E 99, 062126 (2019) — the Cramér–Rao route · Sethna, Statistical Mechanics: Entropy, Order Parameters, and Complexity (2nd ed., 2021)

Provenance. All equations are standard results of stochastic thermodynamics. A verifier should: (1) confirm the master-equation entropy production formula and its positivity (Seifert 2012); (2) reproduce the parabolic large-deviation bound and the TUR (Gingrich, Horowitz, Perunov, England, 2016; conjecture Barato & Seifert, 2015); (3) check the Langevin saturation using the Einstein relation D0=μkBTD_0=\mu k_BT; (4) verify the Cramér–Rao route (Hasegawa & Van Vu, 2019). Note (base-error repair): the dry base's large-deviation derivation stated the Gaussian rate-function curvature as I′′(v)=1/(4D)I''(v)=1/(4D); the correct value is 1/(2D)1/(2D) (since I(j)=(j−v)2/4DI(j)=(j-v)^2/4D), which is what the parabola-matching step actually requires — corrected in-article and in `verify_aq-sm-thermodynamic-uncertainty-relations.py`. All numbers recomputed there.

Builds on: Markov Processes and the Markov Property · The Master Equation · Detailed Balance and Microscopic Reversibility · Fluctuation–Dissipation Theorem · Fokker–Planck Equation Leads to: Langevin Equation and Stochastic Differential Equations · Fluctuation Theorems (Jarzynski–Crooks) · Stochastic Thermodynamics of Molecular Motors

Is the TUR a kind of uncertainty principle like Heisenberg's?

It has the same flavor — a product of two quantities bounded below — but the physics is thermodynamic, not quantum. Heisenberg's relation bounds the product of conjugate-observable spreads by ℏ\hbar; the TUR bounds the product of a current's relative variance and its entropy production by 2kB2k_B. The shared structure is real: both can be derived from Cramér–Rao/Fisher-information arguments, and there are genuinely quantum TURs. But the TUR's floor is set by dissipation, a classical nonequilibrium quantity, and it holds for ordinary thermal noise with no ℏ\hbar in sight.

Where does the factor of exactly 2 come from?

From the Gaussian core of the current statistics. The rate function near the mean is I(j)=(j−v)2/4DI(j)=(j-v)^2/4D, whose curvature is 1/2D1/2D; the universal parabola set by the entropy production has curvature σ/2kBv2\sigma/2k_Bv^2; matching them gives v2/D≤σ/kBv^2/D\le\sigma/k_B, i.e. Var⋅Σ≥2kB⟨J⟩2\mathrm{Var}\cdot\Sigma\ge2k_B\langle J\rangle^2. The 2 is the same 2 that appears in Var(Jt)=2Dt\mathrm{Var}(J_t)=2Dt (the diffusive variance) and in D0=μkBTD_0=\mu k_BT via the Einstein relation. Miscount the curvature as 1/4D1/4D and you would spuriously loosen the bound by a factor of two — a slip worth guarding against.

Can the TUR really tell me the dissipation without knowing the model?

Yes, but only a lower bound. From a trajectory you measure any convenient current's mean and variance and compute σ^min⁡=2kB⟨J⟩2/(Var(J) t)\hat\sigma_{\min}=2k_B\langle J\rangle^2/(\mathrm{Var}(J)\,t). This is guaranteed ≤\le the true dissipation rate, with equality only in the Gaussian near-equilibrium regime. Different currents give different bounds; the tightest is the best estimate. The power is that it needs no knowledge of the hidden network — hence "thermodynamic inference." The limitation is that far from equilibrium the bound can badly underestimate the real dissipation.

Why does a more precise clock have to burn more fuel?

Because precision is a directed signal pulled out of thermal noise, and directing anything against noise dissipates. A clock that ticks N\mathcal N times coherently is a current with small relative fluctuations; the TUR says N≲Δμ/2kBT\mathcal N\lesssim\Delta\mu/2k_BT, so more coherent ticks demand more free energy per cycle. This is why biological clocks like the KaiABC oscillator consume ATP: accuracy is purchased with dissipation, and the TUR sets the exchange rate. A free (undissipating) clock would have Δμ=0\Delta\mu=0 and hence no coherence at all.

Does the TUR contradict the fluctuation theorems, or complement them?

They complement each other. The fluctuation theorems (Jarzynski, Crooks, Gallavotti–Cohen) are equalities about the full distribution of work or entropy production; the TUR is an inequality about just the first two moments of a current. In fact the TUR can be derived as a corollary of the large-deviation structure that also yields the Gallavotti–Cohen symmetry — both live inside the same rate-function framework. The fluctuation theorems tell you the exact odds of a rare entropy-consuming event; the TUR tells you the unavoidable price of a precise typical current.

What changes for periodically driven systems?

The dynamics is no longer a time-homogeneous steady state, so the standard proof does not apply and the bound can be beaten. The driving protocol injects temporal structure that can correlate fluctuations and sharpen a current below the static 2kB2k_B floor. The generalized relations (Barato–Chetrite; Koyuk–Seifert) restore a valid inequality by including the driving period and a driving-induced contribution — the bound becomes weaker or differently shaped, but the second-law spirit survives. This is why one must always check whether the rates are truly time-independent before invoking the standard TUR.

Before this

Where this leads