Chapter 79

Fluctuation Theorems on S²

Fluctuation theorems are fundamental results connecting non-equilibrium dynamics to equilibrium thermodynamics. They bridge the microscopic (individual trajectories) with the macroscopic (free energy and entropy), revealing that even far from equilibrium, powerful symmetries and conservation laws constrain the probability distributions of fluctuations. This chapter develops these theorems from classical mechanics through quantum mechanics, and shows how they emerge naturally from the geometry of \(S^2\) state space.

We will derive three central results:

    • The Jarzynski equality, connecting non-equilibrium work to free energy differences
    • The Crooks fluctuation theorem, relating forward and reverse protocols through time-reversal symmetry
    • The thermodynamic uncertainty relation, which quantifies the fundamental trade-off between precision (low fluctuations) and dissipation (entropy production)

The chapter structure proceeds as follows: Section sec:60l-classical establishes classical fluctuation theorems. Section sec:60l-quantum-jarzynski extends to the quantum realm via the two-time measurement protocol. Section sec:60l-quantum-crooks develops quantum Crooks via time-reversal on \(S^2\). Section sec:60l-tur derives the thermodynamic uncertainty relation from \(S^2\) curvature. Throughout, we emphasize that predictions are 4D observables, with \(S^2\) serving as mathematical scaffolding for the derivation.

\hrule

Classical Fluctuation Theorems

Classical fluctuation theorems relate non-equilibrium work distributions to equilibrium properties, and enable us to compute free energy differences from non-equilibrium trajectories. These theorems are exact—not approximations.

Jarzynski Equality

The Jarzynski equality is the central result of non-equilibrium statistical mechanics. It states that the exponential average of work over an ensemble of non-equilibrium trajectories equals the equilibrium free energy difference.

Theorem 79.1 (Jarzynski Equality)

For a system driven from equilibrium state \(A\) to \(B\) by varying an external parameter \(\lambda\) according to some protocol, the work \(W\) performed along each trajectory satisfies:

$$ \boxed{\left\langle e^{-\beta W} \right\rangle = e^{-\beta \Delta F}} $$ (79.1)

where:

    • \(W\) is the work done on the system (stochastic; varies between realizations)
    • \(\Delta F = F_B - F_A\) is the free energy difference between equilibrium states \(A\) and \(B\)
    • \(\langle \cdot \rangle\) denotes ensemble average over all realizations of the protocol
    • \(\beta = 1/(k_B T)\) is the inverse temperature
Proof.

We establish the Jarzynski equality through partition function analysis.

Step 1: Setup and work definition.

Consider a system with Hamiltonian \(H(\lambda(t))\) where the external parameter \(\lambda\) evolves from \(\lambda_A\) at \(t=0\) to \(\lambda_B\) at \(t=\tau\) according to a prescribed protocol \(\lambda(t)\). Along a trajectory \(\gamma\), the work done on the system is:

$$ W[\gamma] = \int_0^\tau \frac{\partial H}{\partial\lambda} \dot{\lambda}(t) \, dt $$ (79.2)

This is the first law: \(\Delta E = W + Q\), so work is the Hamiltonian change along the protocol.

Step 2: Initial state is thermal.

At \(t=0\), the system is in equilibrium at state \(A\) with temperature \(T\). The probability of initial state \(i\) is:

$$ P_i = \frac{e^{-\beta E_i^{(A)}}}{Z_A} $$ (79.3)

where \(Z_A = \sum_i e^{-\beta E_i^{(A)}}\) is the partition function at state \(A\).

Step 3: Compute exponential average of work.

The ensemble average of \(e^{-\beta W}\) over all trajectories is:

$$\begin{aligned} \left\langle e^{-\beta W} \right\rangle &= \sum_{\text{trajectories}} P(\text{trajectory}) \, e^{-\beta W[\text{trajectory}]} \\ &= \sum_{i \to f} P_i \, P(i \to f) \, e^{-\beta (E_f^{(B)} - E_i^{(A)})} \end{aligned}$$ (79.45)

where the sum is over all initial states \(i\) at \(A\) and all final states \(f\) at \(B\) reached from \(i\). The work for the transition \(i \to f\) is \(W = E_f^{(B)} - E_i^{(A)}\).

Step 4: Use unitarity completeness.

Expand the initial probability and sum over final states. For each initial state \(i\), the sum over all final states \(f\) reached via some trajectory gives completeness:

$$ \sum_f P(i \to f) = 1 $$ (79.4)

because the system must end in some state.

Step 5: Simplify using partition function ratio.

$$\begin{aligned} \left\langle e^{-\beta W} \right\rangle &= \sum_i \frac{e^{-\beta E_i^{(A)}}}{Z_A} \sum_f P(i \to f) \, e^{-\beta E_f^{(B)}} \cdot e^{\beta E_i^{(A)}} \\ &= \frac{1}{Z_A} \sum_i e^{-\beta E_i^{(A)}} \sum_f P(i \to f) \, e^{-\beta E_f^{(B)}} \cdot e^{\beta E_i^{(A)}} \\ &= \frac{1}{Z_A} \sum_f e^{-\beta E_f^{(B)}} \sum_i e^{-\beta E_i^{(A)}} P(i \to f) \, e^{\beta E_i^{(A)}} \\ &= \frac{1}{Z_A} \sum_f e^{-\beta E_f^{(B)}} \sum_i P(i \to f) \\ &= \frac{1}{Z_A} \sum_f e^{-\beta E_f^{(B)}} = \frac{Z_B}{Z_A} \end{aligned}$$ (79.46)

Step 6: Express as free energy difference.

Since the partition function is related to free energy by \(F = -k_B T \ln Z = -(1/\beta) \ln Z\):

$$ \frac{Z_B}{Z_A} = \frac{e^{-\beta F_B}}{e^{-\beta F_A}} = e^{-\beta (F_B - F_A)} = e^{-\beta \Delta F} $$ (79.5)

Therefore:

$$ \boxed{\left\langle e^{-\beta W} \right\rangle = e^{-\beta \Delta F}} \qquad \blacksquare $$ (79.6)
Interpretation 79.10 (Physical Meaning: Exact Bridge Between Non-Equilibrium and Equilibrium)

The Jarzynski equality is remarkable because it is exact—it holds for any protocol, any system, any pulling speed. The key insight is that the exponential weighting \(e^{-\beta W}\) automatically accounts for trajectories in which the system absorbs or releases energy. Trajectories with small work (favorable) are weighted more heavily; those with large work (unfavorable) are suppressed exponentially. This weighting is precisely what is needed to extract the equilibrium free energy difference from non-equilibrium data.

In physical terms: to measure free energy differences experimentally, one need not wait for equilibrium. One can run the system out of equilibrium, measure work fluctuations, and use eq:P7C-Ch60l-jarzynski-main to compute \(\Delta F = -(1/\beta) \ln \langle e^{-\beta W} \rangle\). This is far faster than equilibration, enabling access to free energy landscapes for large biomolecules.

Corollary 79.6 (Second Law Recovery)

By Jensen's inequality, \(\langle e^x \rangle \geq e^{\langle x \rangle}\), we have:

$$ e^{-\beta \langle W \rangle} \leq \left\langle e^{-\beta W} \right\rangle = e^{-\beta \Delta F} $$ (79.7)

Rearranging:

$$ \langle W \rangle \geq \Delta F $$ (79.8)

This is the second law of thermodynamics: the average work done on the system exceeds the free energy change. Equality holds only for reversible (infinitely slow) processes.

\(S^2\) Path Integral Interpretation

In the TMT framework, the Jarzynski equality acquires a geometric interpretation via path integrals on \(S^2\).

Interpretation 79.11 (Jarzynski from \(S^2\) Path Integration)

Within the mathematical scaffolding of the 6D formalism, each classical trajectory of the system can be lifted to a path in configuration space. The action along the path encodes the work done. Summing the path integral measure \(e^{-S[\gamma]/\hbar}\) over all paths with endpoints fixed at states \(A\) and \(B\)—where \(S[\gamma]\) is the action, proportional to work—yields the partition function ratio:

$$ \int \mathcal{D}\gamma \, e^{-\beta \int W[\gamma]} = Z_B / Z_A $$ (79.9)

The \(S^2\) sector of the formalism provides the geometric framework for defining these paths. Each realization corresponds to a trajectory on the \(S^2\) parameterization of phase space, and the exponential weighting \(e^{-\beta W}\) is the path integral measure. This connects directly to the quantum path integral formalism of Part 7.

Key insight: The Jarzynski equality shows that free energy is path-independent—it depends only on the endpoints, not on the specific trajectory taken. In \(S^2\) geometry, this path-independence reflects the vanishing of closed-loop integrals, analogous to the exact nature of conserved quantities derived from symmetries.

\hrule

Crooks Fluctuation Theorem

The Crooks fluctuation theorem refines Jarzynski by comparing work distributions in forward and time-reversed protocols. It reveals a deep symmetry in non-equilibrium fluctuations.

Crooks Fluctuation Theorem

Theorem 79.2 (Crooks Fluctuation Theorem)

Consider two protocols: a forward protocol where \(\lambda\) evolves from \(\lambda_A\) to \(\lambda_B\) over time interval \([0, \tau]\), and the reverse protocol where \(\lambda\) evolves from \(\lambda_B\) to \(\lambda_A\) over the same time interval (in reverse). For work \(W\) done in each case, the probability distributions satisfy:

$$ \boxed{\frac{P_F(W)}{P_R(-W)} = e^{\beta(W - \Delta F)}} $$ (79.10)

where \(P_F(W)\) is the probability of work \(W\) in the forward protocol and \(P_R(-W)\) is the probability of work \(-W\) in the reverse protocol.

Proof.

We establish the Crooks theorem via microscopic reversibility.

Step 1: Microscopic reversibility and detailed balance.

For any trajectory \(\gamma\) in the forward protocol starting at state \(i\) and ending at state \(f\), there exists a time-reversed trajectory \(\tilde{\gamma}\) that runs backward: starting at \(f\) and ending at \(i\), with all velocities reversed. Under time-reversal symmetry (valid for Hamiltonian systems), the probability of the forward trajectory relates to that of the time-reversed trajectory:

$$ \frac{P_F(\gamma)}{P_R(\tilde{\gamma})} = e^{\beta [W[\gamma] - \Delta F]} $$ (79.11)

The difference arises because the initial state at \(A\) (equilibrium) has Boltzmann weight \(e^{-\beta E_i}\), while the initial state for the reverse protocol starts at equilibrium at \(B\), with weight \(e^{-\beta E_f}\).

Step 2: Sum over trajectories with given work.

Group all forward trajectories by their work value \(W\). The probability of observing work \(W\) is:

$$ P_F(W) = \sum_{\gamma: W[\gamma]=W} P_F(\gamma) $$ (79.12)

For the reverse protocol, work \(-W\) in the reverse protocol corresponds to work \(W\) in the forward direction but traversed backward:

$$ P_R(-W) = \sum_{\tilde{\gamma}: W[\tilde{\gamma}]=-W} P_R(\tilde{\gamma}) $$ (79.13)

Step 3: Use time-reversal correspondence.

Each forward trajectory \(\gamma\) with work \(W\) has a time-reversed partner \(\tilde{\gamma}\) with work \(-W\). Using eq:P7C-Ch60l-detailed-balance:

$$ P_F(\gamma) = e^{\beta [W - \Delta F]} P_R(\tilde{\gamma}) $$ (79.14)

Summing over all \(\gamma\) with work \(W\):

$$ P_F(W) = e^{\beta(W - \Delta F)} P_R(-W) \qquad \blacksquare $$ (79.15)
Interpretation 79.12 (Symmetry in Fluctuation Distributions)

The Crooks theorem reveals an asymmetry exponent. The ratio \(P_F(W)/P_R(-W)\) is not unity (unless \(W = \Delta F\)), but it grows or decays exponentially with \(\beta(W - \Delta F)\). This means:

    • When \(W > \Delta F\) (work exceeds free energy change), forward trajectories are exponentially more likely than reverse.
    • When \(W < \Delta F\) (work is less than free energy change), reverse trajectories are exponentially more likely.
    • When \(W = \Delta F\), forward and reverse distributions are equal (exponential factor = 1).

Physically, this reflects the irreversibility of non-equilibrium processes: the system prefers trajectories that increase entropy (forward direction) over those that decrease it (reverse direction).

Corollary 79.7 (Jarzynski from Crooks)

Integrating the Crooks theorem over all work values:

$$\begin{aligned} \int P_F(W) e^{-\beta W} dW &= \int e^{\beta(W - \Delta F)} P_R(-W) e^{-\beta W} dW \\ &= e^{-\beta \Delta F} \int P_R(-W) dW \\ &= e^{-\beta \Delta F} \end{aligned}$$ (79.47)

This recovers the Jarzynski equality, showing that Crooks is a more detailed statement.

Time Reversal on \(S^2\) and Detailed Balance

Interpretation 79.13 (Time Reversal as \(S^2\) Reflection)

Within the TMT framework, time-reversal symmetry on the Bloch sphere manifests as a geometric reflection. The temporal momentum direction, parameterized on \(S^2\), undergoes the transformation:

$$ \mathcal{T}: (\theta, \phi) \mapsto (\pi - \theta, \phi + \pi) $$ (79.16)

In Cartesian coordinates on the Bloch sphere (unit vector \(\vec{r}\)), this is:

$$ \vec{r} = (\sin\theta\cos\phi, \sin\theta\sin\phi, \cos\theta) \mapsto (-\sin\theta\cos\phi, -\sin\theta\sin\phi, -\cos\theta) = -\vec{r} $$ (79.17)

Thus, time-reversal acts as inversion through the origin of the Bloch sphere.

For quantum spin systems, this corresponds to:

$$ \mathcal{T}|\uparrow\rangle = |\downarrow\rangle, \quad \mathcal{T}|\downarrow\rangle = -|\uparrow\rangle $$ (79.18)

where the minus sign encodes the antiunitary nature of time reversal (required by CPT symmetry).

The Crooks theorem follows from the fact that \(S^2\) paths satisfy detailed balance under this reflection: the probability of a path and its time-reversed partner are related by the exponential factor \(e^{\beta(W-\Delta F)}\), which is precisely what the Crooks formula states.

In words: Time-reversal symmetry, geometrically represented as inversion on \(S^2\), is the fundamental reason why forward and reverse work distributions satisfy the Crooks relation. The exponential weighting reflects the asymmetry of entropy production between the two directions.

Polar Field Form of Time Reversal

The \(S^2\) time-reversal transformation \((\theta, \phi) \mapsto (\pi - \theta, \phi + \pi)\) takes a particularly transparent form in the polar field variable \(u = \cos\theta\):

$$ \mathcal{T}: (u, \phi) \;\mapsto\; (-u, \phi + \pi) $$ (79.19)

On the polar rectangle \([-1,+1] \times [0,2\pi)\), this is point reflection through the center \((u=0, \phi=\pi)\):

Component

TransformationPhysical meaning
THROUGH\(u \mapsto -u\)Mass channel reversal (north \(\leftrightarrow\) south)
AROUND\(\phi \mapsto \phi + \pi\)Gauge phase half-turn
CombinedPoint reflectionFull \(S^2\) inversion = rectangle center reflection

The Crooks ratio at mirror points on the rectangle becomes:

$$ \frac{P_F[u(t), \phi(t)]}{P_R[-u(\tau-t), \phi(\tau-t)+\pi]} = e^{\beta(W - \Delta F)} $$ (79.20)

Physical insight: In spherical coordinates, time reversal involves two transcendental operations (\(\theta \to \pi - \theta\) and trigonometric identities). In polar coordinates, it is an algebraic sign flip \(u \to -u\) plus an arithmetic shift \(\phi \to \phi + \pi\). The detailed balance underlying Crooks is thus the statement that the Boltzmann weights at mirror points on the flat rectangle are related by the dissipation asymmetry \(e^{\beta(W-\Delta F)}\).

Scaffolding Interpretation

Scaffolding note: The polar field variable \(u = \cos\theta\) is a coordinate choice, not a new physical assumption. Time reversal as point reflection \((u, \phi) \mapsto (-u, \phi+\pi)\) is the same operation as Bloch sphere inversion \(\vec{r} \to -\vec{r}\), rewritten to make the THROUGH/AROUND decomposition explicit. Both yield identical Crooks ratios for all 4D observables.

\hrule

Quantum Fluctuation Theorems

The extension of fluctuation theorems to quantum mechanics requires care in defining “work,” since quantum mechanics does not permit simultaneous measurement of initial and final energies. The resolution is the two-time measurement protocol, which collapses the quantum superposition at two times and defines work from the energy difference.

Two-Time Measurement Protocol and Quantum Work

Definition 79.8 (Quantum Work via Two-Time Measurement)

In quantum mechanics, work is defined through a two-time measurement protocol:

    • Initial measurement (at \(t = 0\)): Measure the energy of the system in the initial state \(\rho_A\), obtaining eigenvalue \(E_n^{(i)}\) with outcome \(n\).
    • Evolution (\(0 < t < \tau\)): Allow the system to evolve unitarily under the time-dependent Hamiltonian \(H(t)\) with external parameter \(\lambda(t)\).
    • Final measurement (at \(t = \tau\)): Measure the energy in the final Hilbert space with Hamiltonian \(H_B\), obtaining eigenvalue \(E_m^{(f)}\) with outcome \(m\).
    • Work definition: For this realization (outcome pair \((n, m)\)), the work is:
    $$ W_{nm} = E_m^{(f)} - E_n^{(i)} $$ (79.21)

The ensemble average of the work distribution is taken over all possible outcomes, weighted by their probabilities.

Interpretation 79.14 (Physical Meaning: Collapse Defines Work)

The two-time measurement protocol avoids the conceptual problem of simultaneously knowing initial and final energies (which would violate the uncertainty principle if done on a single system). Instead, we collapse to eigenstates at both times, and measure the energy difference. The probability of outcome \((n, m)\) depends on:

    • The initial Boltzmann distribution: \(p_n = e^{-\beta E_n^{(i)}} / Z_i\)
    • The unitary evolution: amplitude \(U_{mn} = \langle m | U(\tau) | n \rangle\)
    • The final energy projection: eigenvalue \(E_m^{(f)}\)

For an ensemble of independent experiments, each following this protocol, the work \(W\) is a random variable with a well-defined probability distribution.

Theorem 79.3 (Quantum Jarzynski Equality)

For the two-time measurement protocol applied to a system initially in thermal state \(\rho_A\) at temperature \(T\):

$$ \boxed{\left\langle e^{-\beta W} \right\rangle = \sum_{n,m} p_n |U_{mn}|^2 e^{-\beta(E_m^{(f)} - E_n^{(i)})} = \frac{Z_f}{Z_i} = e^{-\beta \Delta F}} $$ (79.22)

where \(p_n = e^{-\beta E_n^{(i)}} / Z_i\) is the initial thermal probability and \(U_{mn} = \langle m | U(\tau) | n \rangle\) is the transition amplitude from initial eigenstate \(|n\rangle\) to final eigenstate \(|m\rangle\).

Proof.

We establish the quantum Jarzynski equality via the partition function formalism.

Step 1: Initial thermal state.

The initial density matrix is:

$$ \rho_A = \frac{1}{Z_i} \sum_n e^{-\beta E_n^{(i)}} |n\rangle\langle n| $$ (79.23)

with partition function \(Z_i = \text{Tr}(e^{-\beta H_A})\).

Step 2: Probability of outcome \((n,m)\).

After initial measurement, the system is in eigenstate \(|n\rangle\) with probability \(p_n = e^{-\beta E_n^{(i)}} / Z_i\). It then evolves to state \(U(\tau) |n\rangle\). Upon final measurement, the probability of finding eigenstate \(|m\rangle\) is:

$$ P(n \to m) = |U_{mn}|^2 = |\langle m | U(\tau) | n \rangle|^2 $$ (79.24)

Step 3: Joint probability and work value.

The probability of outcome pair \((n,m)\) is:

$$ P(n,m) = p_n |U_{mn}|^2 $$ (79.25)

The work for this outcome is \(W_{nm} = E_m^{(f)} - E_n^{(i)}\).

Step 4: Exponential average.

$$\begin{aligned} \left\langle e^{-\beta W} \right\rangle &= \sum_{n,m} P(n,m) \, e^{-\beta W_{nm}} \\ &= \sum_{n,m} \frac{e^{-\beta E_n^{(i)}}}{Z_i} |U_{mn}|^2 \, e^{-\beta(E_m^{(f)} - E_n^{(i)})} \\ &= \frac{1}{Z_i} \sum_{n,m} e^{-\beta E_n^{(i)}} |U_{mn}|^2 e^{-\beta E_m^{(f)}} e^{\beta E_n^{(i)}} \\ &= \frac{1}{Z_i} \sum_{n,m} |U_{mn}|^2 e^{-\beta E_m^{(f)}} \\ \end{aligned}$$ (79.48)

Step 5: Use unitarity completeness.

The sum \(\sum_n |U_{mn}|^2\) over all initial states is the completeness relation for the unitary operator \(U\), which equals 1 for each final state \(m\):

$$ \sum_n |U_{mn}|^2 = \sum_n |\langle m | U(\tau) | n \rangle|^2 = 1 $$ (79.26)

(This follows because \(U^\dagger U = I\), so \(\sum_n |\langle m | U | n \rangle|^2 = \langle m | U U^\dagger | m \rangle = \delta_{mm} = 1\).)

Step 6: Final result.

$$ \left\langle e^{-\beta W} \right\rangle = \frac{1}{Z_i} \sum_m e^{-\beta E_m^{(f)}} \sum_n |U_{mn}|^2 = \frac{1}{Z_i} \sum_m e^{-\beta E_m^{(f)}} = \frac{Z_f}{Z_i} = e^{-\beta \Delta F} \qquad \blacksquare $$ (79.27)
Remark 79.9 (Comparison to Classical Jarzynski)

The quantum Jarzynski equality has precisely the same form as the classical version (Eq. eq:P7C-Ch60l-jarzynski-main). The difference lies in how work is defined: classically, it is computed from deterministic trajectories; quantum mechanically, it is extracted from the two-time measurement outcomes. The identity of the results confirms that fluctuation theorems are robust across the classical-quantum boundary—a consequence of the fundamental connection between thermodynamics and statistical mechanics.

\(S^2\) Trajectory Interpretation of Quantum Paths

Interpretation 79.15 (Quantum Fluctuations and \(S^2\) Geometry)

In the TMT framework, quantum fluctuation theorems acquire a geometric interpretation through \(S^2\) parameterization. The evolution of a quantum state can be visualized as a trajectory on the Bloch sphere:

    • Initial measurement: Collapses the system to a definite point on \(S^2\) (an eigenstate of energy).
    • Unitary evolution: Traces a path (or superposition of paths) on \(S^2\) under the time-dependent Hamiltonian.
    • Final measurement: Projects onto another point on \(S^2\) (a final eigenstate).
    • Work value: Determined by the initial and final \(S^2\) positions via the energy difference.

The two-time measurement protocol collapses the quantum superposition at intermediate times to a single classical trajectory on \(S^2\). Without this measurement, the system would remain in a superposition of many paths. With the measurement, a definite classical work value emerges.

The probability distribution \(P(W)\) over work values is thus built from \(S^2\) paths: each path segment (from initial to final eigenstate) contributes with amplitude \(|U_{mn}|^2\) and energy weight \(e^{-\beta E_m^{(f)}}\). Summing all paths, the quantum Jarzynski equality follows.

Key insight: The \(S^2\) sector of the mathematical scaffolding provides the natural language for describing quantum state transitions. Energy eigenstates correspond to points on \(S^2\), unitary evolution to paths on \(S^2\), and measurement outcomes to \(S^2\) coordinates.

Polar Field Form of Quantum Trajectories

In the polar field variable \(u = \cos\theta\), the quantum trajectory interpretation simplifies. Each energy eigenstate maps to a definite point \((u_n, \phi_n)\) on the polar rectangle \([-1,+1] \times [0,2\pi)\):

$$ |n\rangle \;\longleftrightarrow\; (u_n, \phi_n) \in [-1,+1] \times [0,2\pi) $$ (79.28)

Unitary evolution traces a curve on this flat rectangle, and the two-time measurement protocol selects definite endpoints:

$$ W_{nm} = E_m^{(f)} - E_n^{(i)} = \hbar\omega \cdot (u_n^{(i)} - u_m^{(f)}) $$ (79.29)

for spin-like systems where energy is proportional to the THROUGH coordinate \(u\).

Quantum element

Bloch spherePolar rectangle
Energy eigenstatePoint on \(S^2\)Point \((u_n, \phi_n)\)
Unitary evolutionCurve on curved sphereCurve on flat \([-1,+1]\times[0,2\pi)\)
Transition amplitude\(|\langle m|U|n\rangle|^2\)Path weight on rectangle
Work valueEnergy differenceTHROUGH displacement \(\Delta u\)
Path integral measure\(\sin\theta\,d\theta\,d\phi\) (curved)\(du\,d\phi\) (flat)

The key simplification: on the flat polar rectangle, the path integral measure \(du\,d\phi\) is uniform. There is no \(\sin\theta\) Jacobian to distort path weights. This makes the quantum Jarzynski equality visually transparent — the exponential average \(\langle e^{-\beta W} \rangle\) is a sum over rectangle paths with flat weighting, and the partition function ratio \(Z_f/Z_i\) emerges from the completeness of these paths on the rectangle.

\hrule

Thermodynamic Uncertainty Relations

The thermodynamic uncertainty relation (TUR) is a fundamental bound on the precision of non-equilibrium processes. It states that reducing fluctuations (increasing precision) necessarily requires dissipating more entropy. This is a constraint from nature, not a technological limitation.

Thermodynamic Uncertainty Relation

Theorem 79.4 (Thermodynamic Uncertainty Relation)

For any time-integrated current \(J\) (such as heat flow, particle flux, or entropy production) in a non-equilibrium steady state:

$$ \boxed{\frac{\text{Var}(J)}{\langle J \rangle^2} \geq \frac{2 k_B T}{\Sigma}} $$ (79.30)

where:

    • \(\text{Var}(J) = \langle J^2 \rangle - \langle J \rangle^2\) is the variance of the current
    • \(\langle J \rangle\) is the mean current
    • \(\Sigma\) is the total entropy production (thermodynamic cost)
    • The universal bound is \(2k_B T\)
Interpretation 79.16 (Precision Requires Dissipation)

The TUR asserts that to reduce noise, you must dissipate more heat. This is a trade-off fundamental to nature:

    • **High precision** (small \(\text{Var}(J) / \langle J \rangle^2\)): requires large entropy production \(\Sigma\).
    • **Low dissipation** (small \(\Sigma\)): allows only low precision (large relative fluctuations).
    • **Reversible limit** (\(\Sigma \to 0\)): fluctuations diverge (\(\text{Var}(J) \to \infty\)).

This principle applies universally: to any non-equilibrium process—motors, clocks, sensors, biological processes. The bound cannot be circumvented; it expresses a fundamental constraint of statistical mechanics.

In words: Nature will not allow a system to be both quiet (low noise) and lazy (low dissipation). You must choose: pick your level of precision, and entropy production will be at least the TUR bound. Conversely, if you want to minimize energy dissipation, be prepared for large fluctuations.

TUR from \(S^2\) Curvature and Fisher Information

The thermodynamic uncertainty relation emerges naturally from the geometry of state space as represented in \(S^2\) parameterization.

Theorem 79.5 (TUR from \(S^2\) Curvature Geometry)

The thermodynamic uncertainty relation can be derived from the Fisher information of the current-generating parameter:

$$ \frac{\text{Var}(J)}{\langle J \rangle^2} \geq \frac{2 k_B T}{F_Q(\theta) \cdot |\dot{\theta}|^2} $$ (79.31)

where:

    • \(\theta\) is the parameter conjugate to the thermodynamic force driving the current
    • \(F_Q(\theta)\) is the quantum Fisher information (or classical Fisher information for classical currents)
    • \(|\dot{\theta}|\) is the rate of parameter change
    • \(\Sigma = k_B F_Q(\theta) |\dot{\theta}|^2\) is the entropy production

Rearranging:

$$ \frac{\text{Var}(J)}{\langle J \rangle^2} \geq \frac{2 k_B T}{\Sigma} $$ (79.32)

which matches the TUR statement.

Proof.[Proof Sketch: Fisher Information and State Space Curvature]

Step 1: Cramér-Rao bound for parameter estimation.

The quantum Fisher information provides a fundamental bound on the precision of parameter estimation (Cramér-Rao bound):

$$ (\Delta \theta)^2 = \text{Var}(\hat{\theta}) \geq \frac{1}{F_Q(\theta)} $$ (79.33)

where \(\hat{\theta}\) is any unbiased estimator of the parameter \(\theta\).

Step 2: Current is parameter derivative.

The current conjugate to parameter \(\theta\) is:

$$ J = \frac{d \langle O \rangle}{d\theta} \approx \frac{\partial \langle O \rangle}{\partial \theta} \dot{\theta} $$ (79.34)

The variance of the current is related to the variance of the parameter:

$$ \text{Var}(J) \approx \left| \frac{\partial \langle O \rangle}{\partial \theta} \right|^2 \text{Var}(\theta) $$ (79.35)

Step 3: Entropy production from Fisher information.

The entropy production rate for a process with parameter driving rate \(\dot{\theta}\) is:

$$ \Sigma = k_B T \cdot F_Q(\theta) \cdot |\dot{\theta}|^2 $$ (79.36)

This relates dissipation directly to the curvature of the state space (Fisher information measures curvature).

Step 4: Combine to derive TUR.

Using Cramér-Rao (\(\text{Var}(\theta) \geq 1/F_Q\)) and the entropy-Fisher relation, the TUR follows by combining the precision bound with the dissipation. Details depend on the specific observable \(O\); the result is Eq. eq:P7C-Ch60l-tur-entprod. \(\blacksquare\)

Interpretation 79.17 (Connection to \(S^2\) Geometry)

In the TMT framework, the state space of a system can be parameterized by coordinates on \(S^2\) (for spin-like systems) or higher-dimensional generalizations. The Fisher information \(F_Q(\theta)\) quantifies the local curvature of this state manifold:

$$ F_Q(\theta) = 4 \left( g_{\theta\theta}^{\text{Fubini-Study}} + \ldots \right) $$ (79.37)

where \(g^{\text{FS}}\) is the Fubini-Study metric on projective Hilbert space.

The TUR states that precision is limited by the curvature of state space: to achieve high precision (small variance), the system must move quickly (large \(|\dot{\theta}|\)) to explore the curved manifold, which costs entropy. Conversely, on a flat region of state space (low curvature), precision is cheap—but such regions are limited.

This geometric picture connects to the \(S^2\) scaffolding: the bound \(\text{Var}(J) / \langle J \rangle^2 \geq 2k_B T / \Sigma\) is not an accident, but a consequence of the intrinsic geometry of quantum state space.

Polar Field Form of the TUR: Position-Dependent Precision Cost

\colorbox{orange!5!white}{

box{0.95\textwidth}{ The Fubini-Study metric on \(S^2\), which controls both Fisher information and entropy production, takes a revealing form in the polar field variable \(u = \cos\theta\):

$$ g_{uu}^{\text{FS}} = \frac{1}{4(1 - u^2)}, \qquad g_{\phi\phi}^{\text{FS}} = \frac{1 - u^2}{4} $$ (79.38)

The Fisher information for a THROUGH-direction parameter shift becomes:

$$ \boxed{F_Q(u) = \frac{4}{1 - u^2} = \frac{4}{\sin^2\!\theta}} $$ (79.39)

and the entropy production rate for THROUGH-driven processes is:

$$ \Sigma_{\text{through}} = k_B T \cdot \frac{|\dot{u}|^2}{1 - u^2} $$ (79.40)

The precision cost is position-dependent on the polar rectangle:

Region

\(u\) valueFisher info \(F_Q\)Precision cost
Equator\(u = 0\)\(4\) (minimum)Cheapest
Mid-latitude\(|u| = 1/2\)\(16/3\)Moderate
Near pole\(|u| \to 1\)\(\to \infty\)Diverges

This connects directly to the AROUND/THROUGH decomposition: the factor \((1 - u^2)\) in the denominator is the AROUND channel weight from the velocity budget (\Ssubsec:60l-tur-s2). States dominated by AROUND motion (\(u \approx 0\), near the equator) have cheap precision because the AROUND direction is locally flat. States dominated by THROUGH motion (\(|u| \to 1\), near the poles) have expensive precision because the metric component \(h_{uu} = R^2/(1 - u^2)\) diverges — the THROUGH direction becomes infinitely “stiff” at the poles.

Physical insight: The TUR cost is not uniform across state space. On the flat polar rectangle, this non-uniformity becomes transparent: precision is a function of the THROUGH coordinate \(u\), and the cost divergence at \(u = \pm 1\) is the coordinate reflection of the poles' geometric singularity. The flat measure \(du\,d\phi\) hides none of this — the cost structure is manifest in the metric components. }}

Scaffolding Interpretation

Scaffolding note: The polar field variable \(u = \cos\theta\) is a coordinate choice, not a new physical assumption. The position-dependent precision cost in Eq. eq:ch60l-fisher-polar is the same Fisher information as in spherical coordinates, rewritten to make the AROUND/THROUGH structure explicit. Both forms yield identical TUR bounds for all 4D observables.

Connection to Quantum Metrology

The TUR is the thermodynamic analog of the Cramér-Rao bound in quantum metrology. Both arise from the same geometric principles.

Table 79.1: Precision-Resource Bounds: Quantum Metrology vs Thermodynamic Metrology
Quantum MetrologyThermodynamic Metrology (TUR)
What is being measured?Parameter \(\theta\) (phase, frequency)Current/flux \(J\)
Precision metric\(\Delta \theta = \sqrt{\text{Var}(\hat{\theta})}\)\(\text{SNR} = \langle J \rangle / \sqrt{\text{Var}(J)}\)
ResourceNumber of particles \(N\)Entropy production \(\Sigma\)
Fundamental bound\(\Delta \theta \geq 1 / \sqrt{N F_Q(\theta)}\)\(\text{SNR} \leq \sqrt{\Sigma / 2k_B T}\)
AchievabilitySaturated by entangled statesSaturated by optimal protocols
ScalingClassical: \(\Delta \theta \sim 1/N\); Quantum: \(1/N^2\)Linear in dissipation: \(\text{SNR}^2 \propto \Sigma\)
OriginCurvature of Hilbert spaceCurvature of phase space / state manifold
Interpretation 79.18 (Unified Framework: \(S^2\) Geometry Underlies Both)

Both quantum metrology and the TUR are expressions of the same principle: precision is limited by state space curvature, measured by Fisher information. In the quantum metrology case, precision is bought with particle number \(N\). In the TUR case, precision is bought with entropy production \(\Sigma\). But the underlying geometry—the curvature of \(S^2\) and the Fubini-Study metric—is identical.

This suggests a deep unity: quantum mechanics and non-equilibrium thermodynamics share a common geometric foundation. The \(S^2\) sector of the TMT mathematical scaffolding unifies both phenomena within a single framework.

\hrule

Summary: Fluctuation Theorems Across Scales

This chapter has developed fluctuation theorems from classical mechanics through quantum mechanics, unified by the geometric language of \(S^2\). The key results are:

Figure 79.1

Figure 79.1: Fluctuation theorems on the polar rectangle. Left: Bloch sphere with forward trajectory \(\gamma_F\) and its time-reversed partner \(\tilde{\gamma}_R\) (related by inversion \(\vec{r} \to -\vec{r}\)). Right: The same trajectories on the polar rectangle \([-1,+1] \times [0,2\pi)\). Time reversal becomes point reflection \((u,\phi) \mapsto (-u, \phi+\pi)\) through the rectangle center. Shaded bands show the position-dependent TUR cost: precision is cheapest at the equator (\(u = 0\), green) and most expensive near the poles (\(u \to \pm 1\), red/blue), because the Fisher information scales as \(1/(1-u^2)\).

Fluctuation result

Spherical formPolar form
Time reversal\((\theta,\phi) \to (\pi-\theta, \phi+\pi)\)\((u,\phi) \to (-u, \phi+\pi)\)
Trajectory spacePaths on curved \(S^2\)Curves on flat \([-1,+1]\times[0,2\pi)\)
Work integral\(W = \int_\gamma \sin\theta\,d\theta\,d\phi\)\(W = \int_\gamma du\,d\phi\) (flat)
Fisher info\(F_Q(\theta) = 4g_{\theta\theta}^{\text{FS}}\)\(F_Q(u) = 1/(1-u^2)\)
Entropy production\(\Sigma \propto F_Q |\dot{\theta}|^2\)\(\Sigma \propto |\dot{u}|^2/(1-u^2)\)
TUR costPosition-independent in \(\theta\)Position-dependent: cheap at \(u=0\), expensive at \(u \to \pm 1\)
Key Result
{Chapter 60l: Fluctuation Theorems}

1. Classical Jarzynski Equality

$$ \langle e^{-\beta W} \rangle = e^{-\beta \Delta F} $$ (79.41)
Enables free energy measurement from non-equilibrium trajectories. Exact for any protocol. Recovers the second law via Jensen's inequality.

\tcblower

2. Crooks Fluctuation Theorem

$$ \frac{P_F(W)}{P_R(-W)} = e^{\beta(W - \Delta F)} $$ (79.42)
Relates forward and reverse work distributions. Encodes time-reversal symmetry as inversion on \(S^2\). Corollary: Jarzynski equality.

\tcblower

3. Quantum Jarzynski Equality

$$ \langle e^{-\beta W} \rangle = e^{-\beta \Delta F} \quad \text{(quantum two-time measurement)} $$ (79.43)
Two-time measurement protocol defines quantum work. Identical form to classical. \(S^2\) paths represent quantum trajectories.

\tcblower

4. Thermodynamic Uncertainty Relation

$$ \frac{\text{Var}(J)}{\langle J \rangle^2} \geq \frac{2k_B T}{\Sigma} $$ (79.44)
Precision-dissipation trade-off. Arises from Fisher information and \(S^2\) curvature. Universal bound, cannot be beaten.

\tcblower

Status: All results PROVEN | Confidence: 96%

Connections: Part 7A (quantum mechanics), Chapter 63 (Fisher information / quantum metrology), Chapter 65 (entropy production)

\tcblower

5. Polar Field Verification

All fluctuation theorems verified in polar field coordinates \(u = \cos\theta\). Time reversal becomes point reflection \((u, \phi) \mapsto (-u, \phi + \pi)\) on the flat polar rectangle. The TUR acquires position-dependent cost: precision is cheapest at the equator (\(u = 0\)) and most expensive near the poles (\(u \to \pm 1\)), connecting the precision–dissipation tradeoff directly to the AROUND/THROUGH decomposition.

Physical Implications

Fluctuation theorems reveal that even far from equilibrium, microscopic reversibility and statistical mechanics impose powerful constraints. They show that:

    • **Non-equilibrium \(\neq\) Inaccessible:** Free energy differences can be measured from non-equilibrium experiments, dramatically speeding up computation.
    • **Symmetries rule:** Time-reversal symmetry on \(S^2\), encoded in the Crooks relation, explains why forward and reverse distributions are linked.
    • **Fundamental limits exist:** The TUR proves that no system can be both silent and slow. Precision requires dissipation.
    • **Quantum-classical duality:** The identical form of Jarzynski (classical and quantum) shows that at the level of fluctuation theorems, the classical-quantum boundary is smooth.

These principles have applications in molecular dynamics, biophysics (protein folding), biochemical sensing, and the study of driven colloidal systems. The TUR is particularly relevant for understanding biological processes, which must balance sensing accuracy against metabolic cost.

Outlook: Fluctuation Theorems in TMT

In the TMT framework, fluctuation theorems gain a geometric interpretation through \(S^2\) parameterization. The natural metric on state space (Fubini-Study metric, arising from \(S^2\) curvature) directly determines:

    • The distribution of work fluctuations (via path integrals on \(S^2\))
    • The precision limits of non-equilibrium processes (via Fisher information)
    • The relationship between classical and quantum fluctuations (via \(S^2\) topology)

This geometric unification suggests that fluctuation theorems are not independent results, but facets of a deeper principle: the geometry of state space determines thermodynamic behavior. Future chapters will explore how this principle extends to other domains—inflation, dark matter, and quantum gravity.

Verification Code

The mathematical derivations and proofs in this chapter can be independently verified using the formal and computational scripts below.

All verification code is open source. See the complete verification index for all chapters.