RAPID COMMUNICATIONS PHYSICAL REVIEW E 82, 030104共R兲 共2010兲 Unifying approach for fluctuation theorems from joint probability distributions Reinaldo García-García,1 Daniel Domínguez,1 Vivien Lecomte,2,3 and Alejandro B. Kolton1 1 Centro Atómico Bariloche and Instituto Balseiro, 8400 S. C. de Bariloche, Argentina Laboratoire de Probabilités et Modèles Aléatoires (CNRS UMR 7599), Université Pierre et Marie Curie–Paris VI and Université Paris-Diderot–Paris VII, 2 Place Jussieu, 75005 Paris, France 3 DPMC-MaNEP, Université de Genève, 24 Quai Ernest-Ansermet, 1211 Geneva 4, Switzerland 共Received 12 July 2010; published 22 September 2010兲 2 Any decomposition of the total trajectory entropy production for Markovian systems has a joint probability distribution satisfying a generalized detailed fluctuation theorem, when all the contributing terms are odd with respect to time reversal. The expression of the result does not bring into play dual probability distributions, hence easing potential applications. We show that several fluctuation theorems for perturbed nonequilibrium steady states are unified and arise as particular cases of this general result. In particular, we show that the joint probability distribution of the system and reservoir trajectory entropies satisfy a detailed fluctuation theorem valid for all times. DOI: 10.1103/PhysRevE.82.030104 PACS number共s兲: 05.70.Ln, 05.40.⫺a The nonequilibrium stochastic thermodynamics of small systems has attracted a lot of attention. From the experimental side the development of techniques for microscopic manipulation has allowed one to study fluctuations in small systems 关1,2兴. From the theoretical side a group of exact relations known as fluctuation theorems 共FTs兲 关3–8兴 has shed light into the principles governing dissipation and fluctuations in nonequilibrium phenomena, as in driven systems in contact with thermal bath. Formally, the generality of the FTs can be attributed to the way probability distribution functions of different observables behave under time-reversal symmetry-breaking perturbations 共see 关9,10兴 for reviews on FT兲. Oono and Paniconi 关11兴 proposed a phenomenological framework for a “nonequilibrium steady-state 共NESS兲 thermodynamics” aimed at describing fluctuating systems subjected to an external protocol. In this approach, the total exchanged heat during a time interval by a system initially prepared in a NESS is expressed as the sum of two contributions: Qtot = Qex + Qhk. The “excess heat” Qex is, in average, associated with the energy exchange during transitions between steady states while the “housekeeping heat” Qhk corresponds, in average, to the energy we need to supply to maintain the system in a NESS. Hatano and Sasa 关12兴 introduced a formal framework for these phenomenological ideas and derived a FT which extends the second law of thermodynamics for transitions between NESSs controlled by external parameters 共t兲. The Hatano and Sasa FT is applicable to trajectories x共t兲 evolved with a Langevin or more generally a Markovian dynamics starting from an initial condition sampled from a NESS compatible with the initial values 共0兲 of the control parameters. Under identical conditions different FTs were subsequently proposed for NESS, involving the above decomposition of the exchanged heat. We can distinguish between the so-called integral and detailed FTs for systems initially prepared as described above. The integral fluctuation theorems 共IFTs兲 are exact relations for the average over histories of different stochastic trajectory functionals W关x兴, such as 具e−W典 = 1. Examples are the Jarzynski FT for the total work 关7兴, the Hatano-Sasa FT 关12兴, and the Speck-Seifert 1539-3755/2010/82共3兲/030104共4兲 FT 关13兴. The so-called detailed FTs 共DFTs兲 are, on the other hand, exact relations for the probability distribution functions 共PDFs兲 of different observables W, such as P共W兲 / PR共−W兲 = eW, where PR共W兲 corresponds to the timereversed protocol R共t兲 = 共 − t兲, and a NESS initial condition compatible with 共兲. Examples are given by the Crooks relation 关8兴 or Seifert fluctuation theorem 关14兴. While observables satisfying a DFT trivially satisfy an IFT, the opposite is not always true. In many recently formulated DFTs, a modified “dual” PDF P†R共W兲 enters into play 关15,16兴, which corresponds to trajectories in a system with same stationary PDF but with reversed steady probability current 关i.e., such that P†R共W兲 = PR共W兲 when detailed balance holds兴. In general the dual dynamics is different from the original dynamics of the system. Thus, the presence of dual distributions is a strong limitation to the experimental use of such a DFT. A central result of our work is that generalized DFTs can be established without relying on dual probabilities by considering joint probability distributions for different complementary observables, instead of single PDFs. The joint probability distributions arise naturally from the above mentioned separation of two contributions to the total heat exchanged in a NESS. From this joint DFT all the known DFTs and IFTs follow in a straightforward way. Among the fluctuation theorems formulated for Markovian dynamics, the total trajectory entropy production S关x ; 兴 = ln共P关x ; 兴 / P关xR ; R兴兲 plays a fundamental role 关5,6,14兴. Here, P关x ; 兴 共P关xR ; R兴兲 denotes the probability density of trajectory x 共time-reversed trajectory xR兲 in the forward 共backward兲 protocol. We include in P the initial distribution of x. We omit hereafter the final time in all trajectory functionals and use calligraphic symbols to denote functionals and normal symbols to denote their values. The total trajectory entropy production is odd upon time reversal: S关xR ; R兴 = −S关x ; 兴. We show that any decomposition of S in M distinct contributions, S关x ; 兴 = 兺iM Ai关x ; 兴, each of them being odd Ai关xR ; R兴 = −Ai关x ; 兴, has a generating function satisfying the symmetry 030104-1 ©2010 The American Physical Society RAPID COMMUNICATIONS GARCÍA-GARCÍA et al. 冓 冉 M exp − 兺 iAi关x; 兴 = 冓冉 i PHYSICAL REVIEW E 82, 030104共R兲 共2010兲 冊冔 M exp − 兺 共1 − i兲Ai关x; R兴 i 冊冔 Qtot关x; 兴 = −  共1兲 , R with S = 兺 Ai . = DxP关xR ; R兴O关x; 兴e−S关x DxP关x; R兴O关xR ; 兴e dtẋ兵x共x; 兲 − 关xU共x; ␣兲 − f兴其, 共6兲 for the housekeeping heat, and Qex关x; 兴 = − = 具exp兵− 兺Ni 共1 共3兲 − S关x; 兴兲典R R − i兲Ai关x; 兴其典R . R Equation 共2兲 is proved in a similar way, or also follows from Eq. 共1兲 since it is a symmetry for the generating function of the joint distribution P共A1 , . . . , A M 兲. Before considering their particular application for explicit decompositions of S we note that symmetries 共1兲 and 共2兲 are valid for all times for systems prepared in any initial distribution Pi(x共0兲). In particular, we see that the total trajectory entropy production FTs 具e−S典 = 1 and P共S兲 / PR共−S兲 = eS hold without further assumption. We will consider two generic frameworks where our result applies: systems described 共i兲 by Langevin dynamics and 共ii兲 by a Markovian dynamics on discrete configurations, exemplifying their parallel features. First we consider a particle driven by a constant force f in a potential U, ẋ = − xU„x; ␣共t兲… + f + , dtẋx共x; 兲, 共7兲 for the excess heat. The Hatano-Sasa functional 关12兴 Y关x ; 兴 is then defined as Y关x; 兴 ⬅ R;R兴 where we have used S关x ; 兴 = −S关xR ; R兴 together with the change of variable x → xR. Considering S关x ; 兴 = 兺iM Ai关x ; 兴 the proof comes around: = 具exp共− 冕 冕 dt˙ 共x; 兲 = Qex关x; 兴 + ⌬共x; 兲, 共8兲 0 R 具exp共− 冕 0 −S关x;R兴 兺iM iAi关xR ; 兴 共5兲 0 i=1 = 具O关xR ; 兴e−S关x; 兴典R , 兺iM iAi关x; 兴兲典 Qhk关x; 兴 = 共2兲 Note that the result involves no use of dual PDFs. To prove Eq. 共1兲 we start by noting that the average of any observable O关x兴 over trajectories satisfies the relation 冕 冕 dtẋ关xU共x; ␣兲 − f兴. The total exchanged heat in a trajectory can be split as Qtot = Qhk + Qex 关11兴. Defining 共x ; 兲 = −ln SS共x , 兲 from the steady-state probability density SS共x , 兲 at fixed values of = 共␣ , f兲 Hatano and Sasa 关12兴 proposed M P共A1,A2, . . . ,A M 兲 = eS R P 共− A1,− A2, . . . ,− A M 兲 0 where i are arbitrary parameters and 具 . . . 典R denotes average over trajectories in the reversed protocol 关17兴. This symmetry is equivalent to the following generalized DFT for the joint probability of Ai关x ; 兴’s: 具O关x; 兴典 = 冕 共4兲 where ␣共t兲 represents a set of control parameters and 共t兲 is the Gaussian white noise, with 具共t兲典 = 0 and 具共t兲共t⬘兲典 = 2T␦共t − t⬘兲 due to a thermal bath at temperature T. For example, this system can be out of equilibrium when U has periodic boundary conditions. We consider for simplicity a single degree of freedom x, but our results are easily generalized, e.g., in larger dimensions and/or with more particles. For a stochastic trajectory generated by Eq. 共4兲 we define the total exchanged heat as 关18兴 where ⌬共x ; 兲 = (x共兲 ; 共兲) − (x共0兲 ; 共0兲) is a time boundary term. We now assume that the system is initially prepared with a distribution corresponding to a NESS compatible with 共0兲, Pi(x共0兲) = SS(x共0兲 , 共0兲). With the previous definitions of Eqs. 共5兲–共8兲 it is known that the total trajectory entropy production S can be decomposed as the sum of two contributions, in two different ways 关10,11,13,16兴: S = Y + Qhk = ⌬ + Qtot . 共9兲 Similar decompositions also exist for Markovian dynamics: we consider now discrete configurations 兵C其 with transition rates W共C → C⬘ ; 兲 between configurations. They depend on , an external control parameter which may vary in time. The probability density at time t obeys the Markovian dynamics t P共C , t兲 = 兺C⬘W(C⬘ → C ; 共t兲)P共C⬘ , t兲 − r(C ; 共t兲)P共C , t兲, where r共C ; 兲 = 兺C⬘W共C → C⬘ ; 兲 is the escape rate from configuration C. A history of the system is described by the succession of configuration 共C0 , . . . , CK兲 visited by the system, with Ck being visited between tk and tk+1. Starting from initial distribution Pi(C共0兲 , 共0兲), the probability of a history is P关C ; 兴 K W共Ck−1 → Ck , tk兲Pi(C共0兲 , 共0兲), = exp关−兰0 dtr(C共t兲 ; 共t兲)兴兿k=1 meaning that the mean value of an historydependent observable O is given by 具O典 We obtain = 兺Kⱖ0兺C0¯CK兰t0dtK¯兰t02dt1O关C , 兴P关C , 兴. that the total trajectory entropy production S关C ; 兴 = ln共P关C ; 兴 / P关CR ; R兴兲 has a first decomposition S = ⌬ + Qtot with ⌬ = log关Pi(C共0兲 , 共0兲) / Pi(C共兲 , 共兲)兴 and K Qtot = 兺 log k=1 W共Ck−1 → Ck, tk兲 W共Ck → Ck−1, tk兲 . 共10兲 Although there is no natural definition of  we write Qtot to exemplify the parallel with Langevin dynamics. Turning to the first decomposition, let us now assume that the initial distribution is steady state: Pi = SS = e−. One directly checks that the Hatano-Sasa functional Y关C , 兴 = 兰0 dt˙ can be written as 030104-2 RAPID COMMUNICATIONS PHYSICAL REVIEW E 82, 030104共R兲 共2010兲 UNIFYING APPROACH FOR FLUCTUATION THEOREMS… K Y = 关共C, 兲兴0 − 兺 关共Ck, tk兲 − 共Ck−1, tk兲兴. 共11兲 k=1 Let us now derive from a unified view the known FTs. As an intermediate step, it is useful to define a dual trajectory weight P†关x兴 as 关15,16兴 P†关x; 兴 = P关x; 兴e−Q Besides, defining the housekeeping work as K Qhk关C, 兴 = 兺 log k=1 W共Ck−1 → Ck, tk兲 W共Ck → Ck−1, tk兲 K + 兺 共Ck, tk兲 k=1 − 共Ck−1, tk兲, 共12兲 we check that the decomposition S = Y + Qhk holds 关19兴. The parallel between Markovian and Langevin frameworks also appears by specializing to rates W共k → k ⫾ 1兲 = e−共/2兲共Vk⫾1−Vk兲 of jumping on a one-dimensional lattice from site k to k ⫾ 1 in a tilted potential Vk = Uk − kf: in the continuum limit, one recovers the Langevin observables 关19兴. In the first decomposition, Y can be identified with the so-called nonadiabatic contribution 共since it vanishes for quasistatic protocols兲 to the trajectory entropy Sna ⬅ Y, while Qhk 共which is continuously produced in the steady state兲 can be identified with the so-called adiabatic part Sa ⬅ Qhk 关16兴. On the other hand, in the second decomposition of S, ⌬ can be identified with the system contribution Ss ⬅ ⌬, while Qtot can be identified with the reservoir contribution Sr ⬅ Qtot to the total trajectory entropy production. The entropy decompositions of Eq. 共9兲 satisfy the conditions for the application of identity 共1兲 or 共2兲 since each term is odd with respect to time reversal in both decompositions. We can thus write DFTs 共valid for all times 兲 for the joint probabilities as P共Y, Qhk兲 Y+Qhk , R hk = e P 共− Y,− Q 兲 共13兲 P共⌬, Qtot兲 tot = e ⌬+Q . P 共− ⌬,− Qtot兲 共14兲 R The corresponding IFTs can be written as hk hk 具e−Y−Q 典 = 具e−共1−兲YR−共1−兲QR 典R , tot tot 具e−⌬−Q 典 = 具e−共1−兲⌬R−共1−兲QR 典R . 共15兲 共16兲 Here, XR denotes X关x ; R兴, with X representing Y, Qhk, Qtot, or ⌬. From Eqs. 共13兲 and 共14兲 we have, in terms of Ss, Sr, Sa, and Sna, that P共Ss , Sr兲 / PR共−Ss , −Sr兲 = eSs+Sr and P共Sa , Sna兲 / PR共−Sa , −Sna兲 = eSa+Sna. It is worth noting that these relations do not involve dual PDFs, and thus they can be tested for a physical system with a given dynamics. We also note that while one can show that Sa and Sna satisfy each one separately a DFT by using dual PDFs 关16兴, Ss and Sr satisfy a joint DFT although they do not satisfy separately a DFT. 共In the asymptotic long-time limit a DFT for Sr can be formulated for systems with bounded energy under steadystate conditions 关3–6,14兴, while for systems with unbounded energies an “extended fluctuation relation” is necessary 关20兴.兲 hk关x;兴 共17兲 . For Markovian dynamics the dual probability P† corresponds to the dynamics in the so-called dual rates W†共C → C⬘ , 兲 ⬅ e−关共C⬘,兲−共C,兲兴W共C⬘ → C , 兲 which share the same steady state as the original dynamics. In the case of the Langevin dynamics of Eq. 共4兲, it corresponds to trajectories generated by the equation ẋ = −xU†(x ; ␣共t兲) + f † + , with U† = 2 /  − U and f † = −f. This equation also has the same steady state as the original one. Let us stress that the dual dynamics corresponds to trajectories in a different physical system. We now follow Eq. 共17兲 to write the dual joint PDF related to Eq. 共13兲 as hk P†共Y, Qhk兲 = P共Y,− Qhk兲eQ , 共18兲 which is normalized. Integrating this relation over Y, we first hk obtain the DFT 关16兴 P共Qhk兲 = eQ P†共−Qhk兲, and hence hk the IFT 具e−Q 典 = 1 关13兴. Using now successively Eqs. 共13兲 and 共18兲, P共Y兲 = eY = eY 冕 冕 hk d共Qhk兲eQ PR共− Y,− Qhk兲 d共Qhk兲P†R共− Y, Qhk兲, 共19兲 we see that the DFT P共Y兲 = eY P†R共−Y兲 关15兴 holds. This implies the corresponding IFT 具e−Y典 = 1 关12兴 关also derived from setting = 1 and = 0 in Eq. 共15兲 and using the Speck-Seifert IFT兴. Finally, we note that although we have so far assumed steady-state initial conditions, relations 共14兲 and 共16兲 remain valid for an arbitrary initial distribution Pi(C共0兲 , 共0兲) by replacing ⌬ → log关Pi(C共兲 , 共兲) / Pi(C共0兲 , 共0兲)兴. By taking = = 1 in Eq. 共16兲 we thus reobtain the IFT noted by Seifert 关14兴. As an illustration of our approach, let us show how joint FTs provide insights on the experimental error in the evaluation of entropy productions. Consider an experiment where the steady state can be evaluated for different values of the control parameter , e.g., microspheres optically driven in a liquid 关21兴. Having in hand an experimental evaluation exp of , we write S = Ye + ␦Y + Qhk , 共20兲 where Ye关x ; 兴 = 兰0 dt˙ exp and ␦Y = Y − Ye is the difference between exact and experimental Hatano-Sasa functionals. Starting from the NESS associated with 共0兲 and assuming that each of the terms in Eq. 共20兲 is odd upon time reversal, we can use Eq. 共2兲 for M = 3, obtaining P共Y e, ␦Y, Qhk兲 = PR共− Y e,− ␦Y,− Qhk兲eY e+␦Y+Q . hk Using Eq. 共18兲, this gives P共Y e , ␦Y兲 = P†R共−Y e , −␦Y兲eY e+␦Y , and hence also the IFT 030104-3 RAPID COMMUNICATIONS PHYSICAL REVIEW E 82, 030104共R兲 共2010兲 GARCÍA-GARCÍA et al. 具e−Ye典 = 具e−␦YR典R† . 共21兲 From the analysis of 具e−␦YR典R† for specific cases of experimental errors, one can estimate the dispersion of the experimentally obtained 具e−Ye典 around 具e−Y典 = 1 关19兴. As a second example let us consider a system initially prepared in a NESS with a variation of its parameters i共t兲 = 0i + ␦i共t兲 in such a way that 兩␦i共t兲 / 0i 兩 Ⰶ 1, with 0i = i共0兲 and ␦i共0兲 = 0. In this case a modified fluctuationdissipation theorem has been derived in 关22兴, which relates dissipation under small perturbations around a NESS with fluctuations in the corresponding steady state. Expanding the exponentials in Eq. 共15兲 up to second order in ␦ and then to second order in and order zero in , we arrive at 共see 关19兴 for details兲 具Bij共0, 兲典0 = 冓 冋 Bij共,0兲exp −  冕 0 dtẋ共t兲vs„x共t兲; 0… 册冔 , 0 共22兲 where Bij共t , t⬘兲 = 关(x共t兲 ; ) / i兴关(x共t⬘兲 ; ) / j兴 and vs = −1x − 共xU − f兲 corresponds to the average velocity in the NESS associated with . For a Boltzmann-Gibbs steady state, this result reduces to the symmetry Cij共兲 = Cij共−兲, with Cij共兲 = 具Bij共0 , 兲典0. Equation 共22兲 can also be derived from Eq. 共19兲 in Ref. 关15兴. With the use of the joint PDF we can obtain further new results. Let us introduce a weighted correlation function as 0 CW ij 共,0兲 = 冓 冋 0  Bij共0, 兲exp − 冓冋 exp −  2 2 冕 0 冕 0 dtẋ共t兲vs„x共t兲; 0… dtẋ共t兲vs„x共t兲; … 0 册冔 册冔 0 . 0 共23兲 关1兴 C. Bustamante, J. Liphardt, and F. Ritort, Phys. Today 58共7兲, 43 共2005兲. 关2兴 F. Ritort, Adv. Chem. Phys. 137, 31 共2008兲. 关3兴 Denis J. Evans, E. G. D. Cohen, and G. P. Morriss, Phys. Rev. Lett. 71, 2401 共1993兲. 关4兴 G. Gallavotti and E. G. D. Cohen, Phys. Rev. Lett. 74, 2694 共1995兲. 关5兴 J. Kurchan, J. Phys. A 31, 3719 共1998兲. 关6兴 J. L. Lebowitz and H. Spohn, J. Stat. Phys. 95, 333 共1999兲. 关7兴 C. Jarzynski, Phys. Rev. Lett. 78, 2690 共1997兲; Phys. Rev. E 56, 5018 共1997兲. 关8兴 G. E. Crooks, J. Stat. Phys. 90, 1481 共1998兲; Phys. Rev. E 61, 2361 共2000兲. 关9兴 J. Kurchan, J. Stat. Mech.: Theory Exp. 共2007兲, P07005. 关10兴 R. J. Harris and G. M. Schütz, J. Stat. Mech.: Theory Exp. 共2007兲, P07020. 关11兴 Y. Oono and M. Paniconi, Prog. Theor. Phys. Suppl. 130, 29 共1998兲. 关12兴 Takahiro Hatano and Shin-ichi Sasa, Phys. Rev. Lett. 86, 3463 This correlation function carries explicit information about the lack of detailed balance and reduces to the usual one when the system is able to equilibrate. Using Eq. 共15兲 with = 21 and repeating the same procedure used to obtain Eq. W 共22兲 we arrive at the symmetry CW ij 共 , 0兲 = Cij 共0 , 兲, which reduces to the known result for equilibrium dynamics when detailed balance holds. In conclusion, identities 共1兲 and 共2兲 and their immediate consequences for Markovian systems are the main results of this work. Equation 共1兲 or Eq. 共2兲 indeed contains, as particular cases, several known FTs such as the ones previously derived by Hatano and Sasa 关12兴, Speck and Seifert 关13兴, Seifert 关14兴, Chernyak et al. 关15兴, and Esposito and Van den Broeck 关16兴. In addition, an exact DFT, valid for all times , holds for the joint distribution of the reservoir and system entropy contributions to the total trajectory entropy production, although each contribution does not do it separately, as given by Eq. 共14兲. Also a similar DFT holds for the joint distribution of the adiabatic and nonadiabatic entropy contributions to the total trajectory entropy, as given by Eq. 共13兲. For the type of NESS discussed here, two-variable joint PDFs are all that is needed for the corresponding DFTs since there are M = 2 decompositions of the total trajectory entropy production, like Eq. 共9兲, for this case. We have shown an example with M = 3 for handling experimental errors in the Hatano-Sasa FT. In any case, in the light of Eqs. 共1兲 and 共2兲, obtaining an adequate minimal M decomposition of the total trajectory entropy production constitutes the cornerstone toward the derivation of generalized FTs for nonequilibrium systems. This work was supported by CNEA, CONICET 共Grant No. PIP11220090100051兲, and ANPCYT 共Grant No. PICT2007886兲. V.L. was supported in part by the Swiss NSF under MaNEP and Division II, and thanks CNEA for hospitality. 共2001兲. 关13兴 T. Speck and U. Seifert, J. Phys. A 38, L581 共2005兲. 关14兴 U. Seifert, Phys. Rev. Lett. 95, 040602 共2005兲. 关15兴 V. Y. Chernyak, M. Chertkov, and C. Jarzynski, J. Stat. Mech.: Theory Exp. 共2006兲, P08001. 关16兴 M. Esposito and C. Van den Broeck, Phys. Rev. Lett. 104, 090601 共2010兲. 关17兴 A similar splitting of the entropy has been considered in 关6兴 but applied for the case of constant . 关18兴 K. Sekimoto, Prog. Theor. Phys. Suppl. 130, 17 共1998兲. 关19兴 R. García-García, V. Lecomte, A. B. Kolton, and D. Domínguez 共unpublished兲. 关20兴 R. van Zon and E. G. D. Cohen, Phys. Rev. Lett. 91, 110601 共2003兲. 关21兴 E. H. Trepagnier, C. Jarzynski, F. Ritort, G. E. Crooks, C. J. Bustamante, and J. Liphardt, Proc. Natl. Acad. Sci. U.S.A. 101, 15038 共2004兲. 关22兴 J. Prost, J.-F. Joanny, and J. M. R. Parrondo, Phys. Rev. Lett. 103, 090601 共2009兲. 030104-4
© Copyright 2025 Paperzz