Generic Properties of Stochastic Entropy Production Simone Pigolotti1;2,∗ Izaak Neri1;3,y Edgar´ Rold´an1,z and Frank J¨ulicher1;4x 1 Max Planck Institute for the Physics of Complex Systems, N¨othnitzerstraße 38, 01187 Dresden, Germany 2 Biological Complexity Unit, Okinawa Institute for Science and Technology and Graduate University, Onna, Okinawa 904-0495, Japan 3Max Planck Institute of Molecular Cell Biology and Genetics, Pfotenhauerstraße 108, 01307 Dresden, Germany 4Center for Systems Biology Dresden, Pfotenhauerstraße 108, 01307 Dresden, Germany We derive an It^o stochastic differential equation for entropy production in nonequilibrium Langevin processes. Introducing a random-time transformation, entropy production obeys a one- dimensional drift-diffusion equation, independent of the underlying physical model. This transfor- mation allows us to identify generic properties of entropy production. It also leads to an exact uncertainty equality relating the Fano factor of entropy production and the Fano factor of the random time, which we also generalize to non steady-state conditions. PACS numbers: 05.70.Ln, 05.40.-a, 02.50.Le The laws of thermodynamics can be extended to meso- system is in contact with a thermostat at temperature T . scopic systems [1{5]. For such systems, energy changes The stochastic dynamics of the system can be described on the order of the thermal energy kBT are relevant. by the probability distribution P (X;~ t) to find the sys- Here, kB is the Boltzmann constant and T the tem- tem in a configuration X~ at time t . This probability perature. Therefore, thermodynamic observables asso- distribution satisfies the Smoluchowski equation ciated with mesoscopic degrees of freedom are stochas- tic. A key example of such thermodynamics observables @tP = −r~ · J;~ (1) is the stochastic entropy production in nonequilibrium processes. Recent experimental advances in micromanip- where the probability current is given by ulation techniques permit the measurement of stochastic entropy production in the laboratory [6{10]. J~ = µ · FP~ − D · r~ P: (2) Certain statistical properties of stochastic entropy pro- Here we have introduced the force at time t, F~ = duction are generic, i.e., they are independent of the ~ ~ ~ ~ physical details of a system. Examples of such generic −rU(X(t); t) + f(X(t); t), where U is a potential and ~ properties are the celebrated fluctuation theorems, for f is a non-conservative force. We always imply no flux reviews see [2,4,5]. Recently, it was shown that in- or periodic boundary conditions. The state-dependent ~ ~ fima and passage probabilities of entropy production mobility and diffusion tensors, µ(X(t)) and D(X(t)) re- are also generic [11]. Other statistical properties of spectively, are symmetric and obey the Einstein relation entropy production are system-dependent, such as the D = kBT µ. This system can also be represented by a mean value [12{15], the variance [16, 17], the first-passage Langevin equation with multiplicative noise as [21, 28] times of entropy production [18{20] and the large devia- dX~ p tion function [21, 22]. Nevertheless, these properties are = µ · F~ + r~ · D + 2 σ · ξ~ : (3) sometimes constrained by universal bounds [11, 14, 16, dt 17, 23{27]. It remains unclear which statistical properties Here ξ~(t) = (ξ (t); ξ (t); : : : ; ξ (t))T is a Gaussian of stochastic entropy production are generic, and why. 1 2 n white noise with mean hξi(t)i = 0 and autocorrelation 0 0 In this Letter, we introduce a theoretical frame- hξi(t)ξj(t )i = δijδ(t − t ) where h::: i denotes an ensem- work which addresses this question for nonequilibrium ble average. Here and throughout the paper the noise Langevin processes. We identify generic properties of terms are interpreted in the It^osense. The tensor σ obeys entropy production by their independence of a stochastic σσT = D and can be chosen as σ = D1=2. In the It^o arXiv:1704.04061v4 [cond-mat.stat-mech] 11 Aug 2017 variable τ which we call entropic time. We find that the interpretation, the term r~ · D is required for consistency evolution of steady-state entropy production as a func- with Eqs. (1) and (2) as it compensates a noise-induced tion of τ is governed by a simple one-dimensional drift- drift [29]. Examples of systems described by Eq. (31) diffusion process, independent of the underlying model. that we consider in this paper are represented in Fig.1: This allows us to identify a set of generic properties of en- a colloidal particle driven by a constant force along a tropy production and obtain exact results characterizing one-dimensional periodic potential (Fig.1A); a colloidal entropy production fluctuations. particle in a two-dimensional non-conservative force field We consider a mesoscopic system described by n slow pointing in the x direction (Fig.1B); and a chiral active T degrees of freedom X~ = (X1(t);X2(t);:::;Xn(t)) . The Brownian motion in two dimensions [30] (Fig.1C). 2 A) B) C) FIG. 1. Examples of nonequilibrium steady states. (A) Brow- Entropy production, nian particle driven by a constant non-conservative force in a pperiodic 1D sawtooth potential, dX=dt = µ[f − @X U(X)] + Time, ∗ ∗ 2Dξ, with the potential U(x) = (U0 x)=x for x 2 [0; x ] ∗ ∗ and U(x) = U0(1 − x)=(1 − x ) for x 2 [x ; 1]. (B)p 2D FIG. 2. Illustration of the decomposition of stochastic entropy production. In nonequilibrium steady states, the stochastic transport in a forcep field: dX=dt = µf cos(2πY ) + 2Dξx and dY=dt = 2Dξy. (C) Chiral active Brownian mo- entropy production Stot(t)=kB (green) is given by the sum of the monotonously increasing entropic time τ(t) (orange), and tion described by 3 degreesp of freedom: position coordinatesp the martingale process M(t) (blue), see Eq. (11). dX=dt = µf cos(φ) + 2Dξx, dY=dt =p µf sin(φ) + 2Dξy and orientation angle dφ/dt = µφ! + 2D!ξ!. In (B) and (C) U = 0 and f is an external non-conservative force. J~ · D−1 · J~ v = k ; (8) S B P 2 We now discuss the stochastic thermodynamics of the which on average equals the average rate of entropy pro- process described by Eq. (3). In It^o'scalculus, the rate duction, hvSi = hdStot=dti [5, 28]. Entropy fluctua- of change of the potential U(X~ (t); t) is given by It^o's ~ −1 tions are governed by the noise term ξS = ξ · σ · lemma [31]: p J=~ J~ · D−1 · J~ which is a one-dimensional Gaussian white noise with hξ (t)i = 0 and hξ (t)ξ (t0)i = δ(t−t0). dU dX~ h i S S S = @ U + r~ U(X~ (t); t) · + Tr D · r~ r~ U (4) The It^oEq. (7) is equivalent to the Langevin equation dt t dt for entropy production in the Stratonovich interpreta- where Tr denotes the trace and the dots denote tensor tion given in Ref. [33]. For each trajectory generated contractions. In stochastic thermodynamics, the first law by Eq. (3), Eq. (7) generates the corresponding entropy can be expressed as dU = dW + dQ, where dW is the production. From Eq. (7) we can derive several generic work performed on the system and dQ is the mesoscopic properties of stochastic entropy production in nonequi- heat exchanged with the thermostat during a time inter- librium processes. val dt [1]. In It^o'scalculus, the rates of change of work We first discuss properties of nonequilibrium steady and heat are given by [21] states for which @tP = 0. We now calculate the time derivative of e−Stot=kB in steady state. Using It^o'slemma, dW dX~ h i we obtain from Eq. (7) = @ U + f~ · + Tr D · r~ f~ (5) dt t dt −Stot=kB r de vS dQ dX~ h i = − 2 e−Stot=kB ξ ; (9) = −F~ · − Tr D · r~ F~ : (6) dt k S dt dt B The expressions (5) and (6) are the It^oversions of the which reveals that e−Stot=kB is a geometric Brownian stochastic work and mesoscopic heat originally defined by motion with zero drift and time-dependent diffusion Sekimoto using the Stratonovich interpretation [1, 32]. coefficient. The fact that e−Stot=kB has no drift im- −Stot=kB We define the stochastic entropy production Stot=kB plies that e is a martingale process [11, 31, 38]. as the logarithm of the ratio of probabilities of forward Using Stot(0) = 0 the integral fluctuation theorem and time-reversed stochastic trajectories [3, 21, 33]. This he−Stot(t)=kB i = 1 follows immediately from Eq. (9). definition is equivalent to dStot=dt = −(1=T )dQ=dt − In steady state, Eq. (7) can be simplified by introduc- ing the dimensionless entropic time kBd ln P (X~ (t); t)=dt, where the first term can be inter- preted as an exchange of entropy with the reservoir and 1 Z t the second term as a change of system entropy. Using τ = v (X~ (t0)) dt0 ; (10) k S Eq. (6) and It^o'slemma, as in Eq. (4) (see [34]), we ob- B 0 tain the following It^ostochastic differential equation for which is an example of a random time [31]. Note that, the entropy production rate in steady state, vS(X~ (t); t) = vS(X~ (t)) represents the expected rate of entropy production at a given point in dStot p = −2k @ ln P + v + 2k v ξ : (7) ~ dt B t S B S S phase space X(t) and τ thus represents the accumulated expected entropy production. In nonequilibrium situa- Here we define the entropic drift vS(X~ (t); t) ≥ 0 as tions with vS > 0, the entropic time τ(t) is monotonously 3 Entropy production at fixed entropic time Infimum Supremum Number of crossings 0.3 1 1 0.3 0.4 a) b) c) d) 0.25 0.2 0.8 0.8 0.2 0.2 0 0.6 0.6 0510 0.15 0.4 0.4 0.1 0.1 0.2 0.2 0.05 Probability density 0 0 0 0 -5 0501234 01234 012345678910 FIG.