# Thermodynamic Derivation of the Fluctuation Theorem and Jarzynski Equality

Authors: Maarten H. P. Ambaum
Journal: ISRN Thermodynamics
 (2012)

---

## Abstract

A thermodynamic expression for the analog of the canonical ensemble for nonequilibrium systems is described based on a purely information theoretical interpretation of entropy. It is shown that this nonequilibrium canonical distribution implies some important results from nonequilibrium thermodynamics, specifically, the fluctuation theorem and the Jarzynski equality. Those results are therefore expected to be more widely applicable, for example, to macroscopic systems.


---

## Body

## 1. Introduction

The derivations of the fluctuation theorem [1, 2] and the Jarzynski equality [3] appear to depend on the underlying microscopic Hamiltonian dynamics. From this it would follow that these theorems are only relevant to microscopic systems, with their associated definitions of entropy and temperature. In contrast, a statistical mechanical description of macroscopic systems often depends on more general forms of entropy, primarily information entropy [4–6]. Two notable examples from fluid dynamics are the statistical mechanics of point vortices [7] and the statistical mechanics of two-dimensional incompressible flows [8]. In such cases, temperature is defined in terms of the change of entropy with the energy of the system [9] or, equivalently, in terms of the Lagrange multiplier for the energy under the maximization of entropy at a given expectation value of the energy [10].The question is whether for such macroscopic systems we can derive a fluctuation theorem or Jarzynski equality. This is of particular importance for climate science as there are strong indications that the global state of the climate system and, more generally, other components of the Earth system may be governed by thermodynamic constraints on entropy production [11–15]. The theoretical underpinning of those thermodynamic constraints is still lacking. The presence of a fluctuation theorem for such systems would be of great importance.Here we demonstrate that the information-theoretical definition of entropy implies the fluctuation theorem and the Jarzynski equality. It is shown that these results are due to the counting properties of entropy rather than the dynamics of the underlying system. As such, both these results are applicable to a much wider class of problems, specifically, macroscopic systems for which we can define an entropy and which are thermostated in some general sense.The central tenet is that for two statesA and B of a system, defined by two sets of macroscopic parameters, the ratio of the probabilities pB/pA for the system to be in either state is(1)pBpA=exp⁡(ΔABSk),
with ΔABS being the difference in entropy between the states B and A. This is essentially the Boltzmann definition of entropy: entropy is a counting property of the system. The theoretical background can be found in [10], where it is shown that this information theoretical interpretation reproduces the statistical mechanics based on Gibbs entropy but furthermore gives a justification of the Gibbs formulation as a statistical inference problem under limited knowledge of the system. Of note is that the entropy only has meaning in relation to the macroscopic constraints on the system (indicated by the subscripts A and B), constraints which can be arbitrarily complex and prescriptive, as may be needed for systems far from equilibrium. In an information-theoretical setting the previous definition of entropy is equivalent to the principle of indifference: the absence of any distinguishing information between microscopic states within any of the macroscopic states A or B is equivalent to equal prior (prior to obtaining additional macroscopic constraints) probabilities for the microscopic states [16]. Note also that we do not need to specify precisely at this point how the states are counted, or how an invariant measure can be defined on the phase space confined by A or B. The principle of indifference does not imply that all states are assumed equally probable; it is a statement that we cannot a priori assume a certain structure in phase space (such as a precisely defined invariant measure) in the absence of further information. The principle of indifference is not a statement about the structure of phase space; it is a principle of statistical inference and it is the only admissible starting point from an information theoretical point of view.


## 2. A General Form for the Canonical Ensemble

Following Boltzmann, we define the entropySA as the logarithm of the number of states accessible to a system under given macroscopic constraints A. For an isolated system, the entropy is related to the size ΦA of the accessible phase space:(2)SA=kln⁡ΦA.
For a classical gas system, A is defined by the energy U, volume V and molecule number N, the phase space size ΦA is the hyperarea of the energy shell, and it defines the usual microcanonical ensemble. For more complicated systems, where A may include several macroscopic order parameters, the energy shell becomes more confined; in the following we will still refer to the accessible phase space under constraints A as the energy shell. The hyperarea ΦA is nondimensionalised such that ΦA(U)dU is proportional to the number of states between energies U and U+dU. We will not consider other multiplicative factors which make the argument of the logarithm nondimensional; these contribute an additive entropy constant which will not be of interest to us here. Note also that the microcanonical ensemble does not include a notion of equilibrium: the system is assumed to be insulated, so it cannot equilibrate with an external system. It just moves around on the energy shell (defined by A) and the principle of indifference implies that all states, however improbable from a macroscopic point of view, are members of the ensemble. Of course, the number of unusual states (say, with nonuniform macroscopic properties not defined by A) is much lower than the number of regular states (say, with uniform macroscopic density) for macroscopic systems. Only for small systems, the distinction becomes important but it does not invalidate the previous formal definition of entropy. The previous definition of entropy also ensures that entropy is an extensive property such that for two independent systems considered together the total entropy is the sum of the individual entropies, S=S1+S2. The Boltzmann constant k ensures dimensional compatibility with the classical thermodynamic entropy when the usual equilibrium assumptions are made [10, 17].The hyperarea of the energy shell, and thus the entropy, can be a function of several variables which are set as external constraints, such as the total energyU, system volume, V, or particle number N for a simple gas system. For the canonical ensemble we consider a system that can exchange energy with some reservoir. We consider here only a theoretical canonical ensemble in that we consider the coupling between the two systems to be weak such that the interaction energy vanishes compared to the relevant energy fluctuations in the system.First, we need to define what a reservoir is. Following equilibrium thermodynamics, we formally define an inverse temperatureβ=(kT)-1 as(3)β=1k∂S∂U=1Φ∂Φ∂U.
We make no claim about the equality of β and the classical equilibrium inverse temperature; β is the expansivity of phase space with energy and as such can be defined for any system, whether it is in thermodynamic equilibrium or not. When an isolated system is prepared far from equilibrium (e.g., when it has a local equilibrium temperature which varies over the system), then β is still uniquely defined for the system as a nonlocal property of the energy shell that the system resides on. Because both energy and entropy in the weak coupling limit are extensive quantities, β must be an intensive quantity.Now consider a large isolated systemR with total (internal) energy UR. Let this system receive energy U′ from the environment. By expanding its entropy SR in powers of U, we can then write the entropy of this large system as(4)SR(UR+U′)=SR(UR)+kU′(β+12U′∂β∂U+O(U′2)).
We see that for finite U′, (∂β/∂U)-1 has to be an extensive quantity. But that means that for a very large system ∂β/∂U=𝒪(N-1), where N is a measure of the size of the system (such as particle number). For a classical thermodynamic system ∂β/∂U=-kβ2/CV with CV the heat capacity at constant volume. We conclude that for a very large system (N→∞), the entropy equals(5)SR(UR+U′)=SR(UR)+kβU′
for all relevant, finite energy exchanges U′. This expression for the entropy defines a reservoir. The size of the energy shell accessible to the reservoir is, for all relevant energy exchanges U′, exactly proportional to exp(βU′), with β an intensive and constant property of the reservoir. We do not require the reservoir to be in thermodynamic equilibrium. A change of energy in the reservoir pushes the reservoir to a different energy shell A′; the functional dependence of the size of the energy shell with energy defines the inverse temperature β, as in (3). However, it is not assured that a small and fast thermometer would measure an inverse temperature equal to β at some point in the reservoir; only if the reservoir is allowed to equilibrate, its inverse temperature is everywhere equal to β. Of course, this is how the temperature of a classical reservoir is determined in practice.Now suppose that a system of interest has energyU0. We then allow it to exchange heat U with a reservoir. If the system has energy U0+U, the reservoir must have given up energy U. We can write the hyperarea of the energy shell of the system Φ0 as a function of U. The total entropy of the system plus reservoir R can then be written as a function of the exchange energy, U, as(6)S=S0(U)+SR(UR)-kβU,
with S0=klnΦ0. The number of states at each level of exchange energy therefore is proportional to(7)Φ(U)∝Φ0(U)exp⁡(-βU),
where we omitted proportionality constants related to the additive entropy constants. Nowhere we assume that the system is in equilibrium with the reservoir. This means that Φ(U) is the relevant measure to construct an ensemble average for the system, even for far-from-equilibrium systems. Even the reservoir can be locally out of equilibrium, as discussed previously. We have also made no reference to the size of the system of interest, as long as it is much smaller than the reservoir. However, in contrast to systems in thermodynamic equilibrium, there is no guarantee that the extensive macroscopic variables, such as U, V, or N, define the state of the system in any reproducible sense. To fully define an out-of-equilibrium system we need to introduce order parameters that can describe the nonequilibrium aspects of the system.The previous density is an integrated version of the usual canonical distribution. The size of the energy shell of the system of interest,Φ0, can be written as an integral over states Γ such that(8)Φ0(U)=∫H0(Γ)=UdΓ,
with H0 being the Hamiltonian of the system of interest. With this definition, the density in (7) reduces to the usual canonical distribution exp(-βH0(Γ)) for states Γ. We will not make further use of this microscopic version of the density.


## 3. Fluctuation Theorems

The canonical density in (7) can be expanded by parametrizing each energy shell with some continuous coordinate υ so that every part of phase space has coordinates (U,υ). The coordinate v is again a macroscopic coordinate so that any combination (U,υ) can correspond to many microscopic states. At each value of υ the differential ϕ(U,υ)dUdυ is proportional to the number of states between coordinate values U and U+dU, and υ and υ+dυ, and it is normalised such that(9)∫ϕ(U,υ)dυ=Φ0(U).
The parametrisation is arbitrary at this point and can be chosen such as to divide the phase space in as fine a structure as desired for a given application. We can define an entropy S0(U,υ) again as the logarithm of the number of available states for the system of interest corresponding to the subset of phase space defined by (U,υ):(10)S0(U,υ)=kln⁡ϕ(U,υ).Now consider a process that occurs on the energy shellU where some variable changes from A→B. On the parametrized energy shell this corresponds to a coordinate shift from υ(A)→υ(B). The number of corresponding states changes from ϕ(U,υ(A))→ϕ(U,υ(B)). We can use detailed balance to express the ratio of the probability of making this transition to the probability of making the reverse transition as the ratio of the number of states at (U,υ(A)) to the number of states at (U,υ(B)):(11)pA→BpB→A=ϕ(U,υ(B))ϕ(U,υ(A))=exp⁡(ΔABSk),
where ΔABS/k=S0(U,υ(B))-S0(U,υ(A)). If, in addition, during the process A→B the energy of the system of interest changes from UA→UB through exchange with the reservoir, then the previous ratio of probabilities can still be expressed as exp(ΔABS/k) but now with(12)ΔABS=S0(UB,υ(B))-S0(UA,υ(A))-kβ(UB-UA).
We can always write the entropy change of the system of interest as the sum of the entropy change due to heat exchange with the reservoir and an irreversible entropy change associated with uncompensated heat [14, 18], namely, S0(UB,υ(B))-S0(UA,υ(A))=kβ(UB-UA)+ΔiS0. We thus conclude that ΔABS=ΔiS0; that is, the relevant entropy change in (11) equals the irreversible entropy change of the system of interest. So for processes that occur either on or across energy shells, we have(13)pA→BpB→A=exp⁡(ΔiS0k),
with ΔiS0 being the irreversible entropy change of the system in a process A→B. The right-hand side of this equation is only dependent on the irreversible entropy change ΔiS0 between the two states of the system of interest. So this equation must be true for any pair of states (A,B) that are related by the same irreversible entropy change. We thus arrive at the fluctuation theorem [1, 2]:(14)p(ΔiS)p(-ΔiS)=exp⁡(ΔiSk),
with p(ΔiS) being the probability that the system of interest makes a transition with irreversible entropy change of ΔiS and p(-ΔiS) being the probability for the opposite change.The fluctuation theorem applies to spontaneous processes that occur in thermostated but otherwise isolated systems. We next consider processes that occur when we modify the system of interest by changing some external macroscopic parameters. The entropy of the energy shellU is then also a function of some parameter λ, namely, S=Sλ(U,υ). Without loss of generality we set λ=0 at A and λ=1 at B. In this case the irreversible entropy change in (13) is(15)ΔiSk=S1(UB,υ(B))-S0(UA,υ(A))-kβ(UB-UA).
Apart from this, there is no change in the considerations leading to the fluctuation theorem. By definition, thermostated systems that receive work WAB from their environment have an irreversible entropy change equal to(16)ΔiSk=β(WAB-ΔABF),
with ΔABF being the change in free energy going from A to B. Recognising that the right-hand side is again only a function of the difference between the two states, we arrive at the Crooks fluctuation theorem [19]:(17)p01(W)p10(-W)=exp⁡(β(W-Δ01F)),
with p01(W) being the probability that the system absorbs work W when λ changes from 0 to 1, and p10(-W) being the probability that the system performs work W when λ changes in reverse from 1 to 0. Because the transition probabilities can be normalised with respect to the exchanged work, it is straightforward to use this equation to show that the expectation value of exp(-β(W-Δ01F)) equals unity, or equivalently,(18)〈exp⁡(-βW)〉=exp⁡(-βΔ01F).
This is the Jarzynski equality [3].The consistency of the previous argument is strengthened by the following independent route to calculate free energy changes. The phase space measureΦ(U) can be normalised with the partition function Zλ:(19)Zλ=∫Φλ(U)exp⁡(-βU)dU,
where Φλ(U) is proportional to the number of accessible states of the isolated the system of interest when the external parameter is set to λ. The equilibrium free energy for the thermostated system is(20)Fλ=-β-1ln⁡Zλ.
Next we consider what happens to the equilibrium free energy of the system when we vary λ from 0 to 1. The partition function at λ=1 satisfies(21)Z1=∫Φ1(U)exp⁡(-βU)dU=∫Φ0(U)exp⁡(ΔSk)exp⁡(-βU)dU=Z0〈exp⁡(ΔSk)〉,
where 〈.〉 denotes an ensemble average over the initial ensemble, and ΔS=kln(Φ1(U)/Φ0(U)). As before, the entropy change can be written as the sum of the entropy change due to heat exchange with the reservoir and the irreversible entropy change due to uncompensated heat. Because the system plus the reservoir is thermally insulated, any heat given to the reservoir must be compensated by work performed by the external parameter change. The entropy change can therefore be written as 〈exp(ΔS/k)〉=〈exp(ΔiS/k-βW)〉 so that we find(22)Z1Z0=〈exp⁡(ΔiSk-βW)〉.
Because (16) is true for any microscopic realisation of the process, we find that the right-hand side of the previous equation is the same for every realisation and it is equal to exp(-βΔF). This is consistent with the equilibrium expression for the free energy, (20), from which follows that exp(-βΔF)=Z1/Z0. The previous equation is only apparently in contradiction to the Jarzynski equality, (18). To arrive at the Jarzynski equality we recognise that (16) implies that 〈exp(β(ΔF-W))〉=〈exp(-ΔiS/k)〉=1, where the last equality follows from integrating the fluctuation theorem over all values of ΔiS.


## 4. Discussion

We have shown that the fluctuation theorem (14) and Jarzynski equality (18) follow from general counting properties of entropy and not from the underlying dynamics. As such we expect both results to be widely applicable to systems that are in some sense thermostated, that is, systems that are able to settle on a given expectation value for the total energy by interaction with a reservoir.The climate system is potentially a nontrivial example of such a system: the incoming short-wave radiation from the Sun is depleted by long-wave (thermal infrared) radiation from the Earth to space. The corresponding equilibrium temperature is the bolometric temperature of the planet (about 255 K in case of the Earth [14]). (The bolometric radiation temperature of the Earth is substantially lower than the observed average surface temperature of about 288 K, because of the greenhouse effect of the atmosphere.) It is not obvious how to apply the fluctuation theorems to the climate system and how the entropy production in the climate system is related to the actual climate on Earth. For example, most of the entropy production in the climate system is due to degradation of radiation (e.g., [20]); namely, short wavelength visible sunlight is thermalized by molecular absorption into molecular thermal energy corresponding to long wavelength infrared radiation. This degradation of radiative energy is the main source of entropy production in the climate system, but as this entropy production only resides in the photon field, its relation to, for example, kinetic energy dissipation in the atmosphere is not clear. So from this example it appears that we need to select the relevant forms of entropy production before we can use it to make inferences about the climate system.It remains to be seen whether the fluctuation theorems can be usefully applied to complex systems such as the climate, but we believe that the derivation presented here can pave the way for attempts in that direction.



---

*Source: 528737-2012-05-17.xml*