arxiv: v4 [math.pr] 26 Apr 2017

Transcription

1 Cycle symmetry, limit theorems, and fluctuation theorems for diffusion processes on the circle Hao Ge,2, Chen Jia 3,4, Da-Quan Jiang 3,5 Biodynamic Optical Imaging Center, Peking University, Beijing 87, P.R. China 2 Beijing International Center for Mathematical Research, Peking University, Beijing 87, P.R. China 3 LMAM, School of Mathematical Sciences, Peking University, Beijing 87, P.R. China 4 Department of Mathematical Sciences, The University of Texas at Dallas, Richardson, Texas 758, U.S.A. 5 Center for Statistical Science, Peking University, Beijing 87, P.R. China arxiv:3.223v4 [math.pr] 26 Apr 27 haoge@pku.edu.cn (H. Ge), jiac@pku.edu.cn (C. Jia), jiangdq@math.pku.edu.cn (D.Q. Jiang) Abstract Cyclic structure and dynamics are of great interest in both the fields of stochastic processes and nonequilibrium statistical physics. In this paper, we find a new symmetry of the Brownian motion named as the quasi-time-reversal invariance. It turns out that such an invariance of the Brownian motion is the key to prove the cycle symmetry for diffusion processes on the circle, which says that the distributions of the forming times of the forward and backward cycles, given that the corresponding cycle is formed earlier than the other, are exactly the same. With the aid of the cycle symmetry, we prove the strong law of large numbers, functional central limit theorem, and large deviation principle for the sample circulations and net circulations of diffusion processes on the circle. The cycle symmetry is further applied to obtain various types of fluctuation theorems for the sample circulations, net circulation, and entropy production rate. Keywords: excursion theory, Bessel process, Haldane equality, cycle flux, nonequilibrium Mathematics Subject Classifications: 6J6, 58J65, 6J65, 6F, 82C3 Introduction Cyclic Markov processes are widely used to model stochastic systems in physics, chemistry, biology, meteorology, and other disciplines. By a cyclic Markov process, we mean a Markov process whose state space has a cyclic topological structure. In general, we mainly focus on two types of cyclic Markov processes: nearest-neighbor random walks on the discrete circle, also called nearestneighbor periodic walks, and diffusion processes on the continuous circle. The path of a cyclic Markov process constantly forms the forward and backward cycles and the cycle dynamics of cyclic Markov processes has been studied for a long time. Among these studies, many authors noticed that nearest-neighbor periodic walks have an interesting symmetry [ 4]. Let T + and T be the forming times of the forward and backward cycles by a nearest-neighbor periodic walk, respectively. Samuels [] and Dubins [2] proved that although the distributions of T + and T may be different, their distributions, given that the corresponding cycle is formed earlier than its reversed cycle, are the same: P (T + u T + < T ) = P (T u T < T + ), u. () This equality, which characterizes the symmetry of the forming times of the forward and backward cycles, is named as the cycle symmetry in this paper. Recently, Jia et al. [5] further generalized

2 the cycle symmetry to general discrete-time and continuous-time Markov chains based on some non-trivial equalities about taboo probabilities. In the previous studies, the cycle symmetry is mainly discussed for nearest-neighbor periodic walks. Besides, Kendall [6] and Pitman and Yor [7] have proved the cycle symmetry for Brownian motion with constant drift. It is then natural to ask whether the cycle symmetry also holds for the continuous version of nearest-neighbor periodic walks, namely, diffusion processes on the circle. Intuitively, the answer should be affirmative. A natural idea is to approximate diffusion processes by random walks. However, this idea seems to obscure the essence of the cycle symmetry to some extent. In this paper, we prove the cycle symmetry for diffusion processes on the circle using the intrinsic method of stochastic analysis, instead of approximation techniques. We find that the essence of the cycle symmetry lies in a new symmetry of the Brownian motion named as the quasi-time-reversal invariance. Based on Brownian excursion theory [8] and Williams Brownian paths decomposition theorem [9, ], we prove that the Brownian motion is invariant under a transformation called quasitime-reversal (see Definition 2.). In this paper, we use this new invariance of the Brownian motion to prove the cycle symmetry for diffusion processes on the circle. Let X be a diffusion process on the unit circle. Let T + and T be the forming times of the forward and backward cycles, respectively. Then the cycle symmetry for diffusion processes on the circle can also be formulated as (). Let T = T + T be the time needed for X to form a cycle for the first time. Then the cycle symmetry can be rewritten as P (T u T + < T ) = P (T u T < T + ), u. (2) This implies that the time needed for X to form a forward or backward cycle is independent of which one of these two cycles is formed. This independence result is important because it shows that the cycle dynamics for any diffusion process on the circle with any initial distribution is nothing but a renewal random walk (see Theorem 4.4). The cycle symmetry established in this paper has some interesting applications. In the literature, the most important functionals associated with a cyclic Markov process are the sample circulations. In fact, the circulation theory of nearest-neighbor random walks and general Markov chains have been well established [, 2]. In this paper, we develop the circulation theory for diffusion processes on the circle. Let N t + and Nt denote the numbers of the forward and backward cycles formed by the diffusion process X up to time t, respectively. Then the sample circulations J t + and Jt of X along the forward and backward cycles are defined as J t + = t N t +, J t = t N t, (3) respectively, and the sample net circulation is defined as J t = J t + Jt. In this paper, we prove the strong law of large numbers, functional central limit theorems, and large deviation principle for the sample circulations with the aid of the cycle symmetry. Over the past two decades, the studies on fluctuation theorems (FTs) for Markov processes have become a central topic in nonequilibrium statistical physics [3, 4]. So far, there has been a large amount of literature exploring various types of FTs [5 2]. In recent years, the FTs for the sample circulations and net circulation of cyclic Markov processes have been increasingly concerned [3, 4, 22]. Interestingly, the cycle symmetry turns out to be a useful tool to study the FTs for diffusion 2

3 processes on the circle. In this paper, we prove various types of FTs for the sample circulations and net circulation, including the transient FT, the Kurchan-Lebowitz-Spohn-type FT, the integral FT, and the Gallavotti-Cohen-type FT. These FTs characterize the symmetry of X along the forward and backward cycles from different aspects. In particular, we prove that the sample circulations (J t +, J t ) satisfy a large deviation principle with rate t and good rate function I with the following highly non-obvious symmetry: I (x, x 2 ) = I (x 2, x ) γ(x x 2 ), x, x 2 R, (4) where γ is a constant characterizing whether X is a symmetric Markov process (see Theorem 6.). Moreover, we prove that the sample net circulation J t also satisfies a large deviation principle with rate t and good rate function I 2 with the following symmetry: I 2 (x) = I 2 ( x) γx, x R. (5) In nonequilibrium statistical physics, a central concept is entropy production rate. In this paper, we also prove that the sample entropy production rate of diffusion processes on the circle satisfies the Gallavotti-Cohen-type FT. The structure of this paper is organized as follows. In Section 2, we prove the quasi-time-reversal invariance of the Brownian motion using Brownian excursion theory. In Section 3, we apply such an invariance to prove the cycle symmetry for diffusion processes on the circle. In Section 4, we show that the cycle dynamics of any diffusion process on the circle can be represented as a renewal random walk. Sections 5 is devoted to the limit theorems and large deviations for the sample circulations. Sections 6 is devoted to the applications of the cycle symmetry in nonequilibrium statistical physics, including various types of FTs for the sample circulations, net circulation, and entropy production rate. Section 7 is an appendix including some detailed proofs. 2 Quasi-time-reversal invariance of Brownian motion In this section, we work with the canonical version of one-dimensional Brownian motion. Let W denote the Winner space, that is, the set of all continuous paths with values in R, let P denote the Winner measure, and let F denote the Borel σ-algebra on W completed with respect to P. Let W = {W t : t } be the coordinate process on W defined as W t (w) = w(t). Let {F t : t } be the natural filtration of W completed with respect to P. Then W is a standard Brownian motion defined on the filtered space (W, F, {F t }, P ) satisfying the usual conditions. For each w W, let τ(w) = inf{t : w(t) = } (6) be the hitting time of {, } by w and let g τ (w) = sup{ t < τ(w) : w(t) = } (7) be the last zero of w before τ(w). 3

4 Definition 2.. The quasi-time-reversal is a transformation φ on W defined as w(t), if t < g τ (w), w(g τ + τ t), if w(τ) = and g τ (w) t < τ(w), φ(w)(t) = w(g τ + τ t) +, if w(τ) = and g τ (w) t < τ(w), w(t) 2, if w(τ) = and t τ(w), w(t) + 2, if w(τ) = and t τ(w), w W. The definition of the quasi-time-reversal is somewhat complicated at first glance, but its intuitive implication is rather clear. Under the quasi-time-reversal, a continuous path w is taken time reversal between g τ (w) and τ(w) and is continuously spliced when needed (see Figure ) time reversal and translation translation g τ τ Figure. A Brownian path and its quasi-time-reversal. The blue curve represents a Brownian path and the red curve represents the transformed path under the quasi-time-reversal φ. During [, g τ ], the two curves coincide. During [g τ, τ], the blue curve is reversed in time and is translated to make the red curve continuous. During [τ, ), the blue curve is translated to make the red curve continuous. In fact, the expression of the quasi-time-reversal becomes much simpler if it is regarded as a transformation on W(S ). Specifically, let S = {z C : z = } denote the unit circle in the complex plane. Let W(S ) denote the Winner space on S, that is, the set of all continuous paths with values in S. Then the map w e 2πiw gives a one-to-one correspondence between W and W(S ). Let w be a path in W(S ). Under this one-to-one correspondence, τ(w) is the time needed for w to form a cycle for the first time and g τ (w) is the last exit time of by w before τ(w). Moreover, the quasi-time-reversal can be viewed as a transformation φ on W(S ) defined as w(t), if t < g τ (w), φ(w)(t) = w(g τ + τ t), if g τ (w) t < τ(w), w W(S ). w(t), if t τ(w), The following theorem shows that the Brownian motion is invariant under the quasi-time-reversal. Theorem 2.2. Let W be a standard one-dimensional Brownian motion. Let φ be the quasi-timereversal on W. Then φ(w ) is also a standard one-dimensional Brownian motion. 4

5 In order to prove the above theorem, we need some knowledge on Brownian excursion theory. We next introduce some notations. For each w W, let R(w) = inf{t > : w(t) = } be the first zero of w after time. Let U be the set of all continuous paths w such that < R(w) < and w(t) = for each t R(w). Let U be the Borel σ-algebra on U. In addition, let n be the Itô measure. Then n is a measure on the measurable space (U, U ). For the rigorous definition of the Itô measure n, please refer to [23, Chapter XII]. For each w W, let i (w) be the element u U such that w(t), if t < R(w), u(t) =, if t R(w). For each t, let i t (w) = i (θ t (w)), where θ t is the shift operator defined as θ t (w)(s) = w(t + s), and let g t (w) = sup{ s < t : w(s) = } be the last zero of w before t. It is easy to see that i gt is the excursion straddling the given time t. The next lemma gives the structure of the Brownian motion W between g τ and τ. In this paper, a three-dimensional Bessel process starting from x will be abbreviated as BES 3 (x). Lemma 2.3. Let ρ and ρ be two independent BES 3 (). Let ξ be a Bernoulli random variable independent of ρ and ρ with P (ξ = ) = P (ξ = ) = /2. Let H be a process defined as ρ t, if ξ =, H t = ρ t, if ξ =. Let β = inf{t : H t = }. Then the processes {W gτ +t : t τ g τ } and {H t : t β} have the same distribution. Proof. Let T and S be the hitting times of and by W, respectively. Obviously R, S, T, and τ can be viewed as defined on U. For each u U, let M(u) = max t R(u) u(t). It is a classical result [23, Chapter XII, Proposition 3.6] that for any x >, This shows that n(m x) = 2x. n(τ < R) = n(t < R) + n(s < R) = 2n(T < R) = 2n(M ) =. (8) Let C be a nonempty closed set in R that does not contain. Let T C be the hitting time of C by W and let g Tc (w) = sup{ s < T c : w(s) = } be the last zero of w before T C. It is a classical result [23, Chapter XII, Proposition 3.5] that for any nonnegative measurable function F on (U, U ), E[F (i gtc ) F gtc } = n(f {TC<R})/n(T C < R). (9) 5

6 It follows from (8) and (9) that for any Γ U, P (i gτ and that Γ) = n(γ, τ < R) = n(γ, T < R) + n(γ, S < R) = n(γ, T < R) + n( Γ, T < R). n(γ, T < R) = n(t < R)P (i gt Γ) = 2 P (i g T Γ). Thus we have P (i gτ Γ) = 2 [P (i g T Γ) + P (i gt Γ)]. Let B and B be two independent standard Brownian motions defined on some probability space (Ω, F, P ) and let T and T be their hitting times of, respectively. Let e and ẽ be the excursions of B and B straddling T and T, respectively. Let η be a Bernoulli random variable independent of B and B with P (η = ) = P (η = ) = /2. Let J be the process defined as e, if η =, J = ẽ, if η =. Then we have P (J Γ) = 2 [P (e Γ) + P (ẽ Γ)] = 2 [P (i g T Γ) + P (i gt Γ)] = P (i gτ Γ). This shows that i gτ and J have the same distribution. Thus the processes {W gτ +t : t τ g τ } and {J t : t γ} have the same distribution, where γ is the hitting time of {, } by J. By Williams Brownian paths decomposition theorem [9, ], the excursion e before reaching is a BES 3 () before reaching. Thus the processes {J t : t γ} and {H t : t β} have the same distribution. Lemma 2.4. The processes {W t : t g τ }, {W gτ +t : t τ g τ }, and {W τ+t W τ : t } are independent. Proof. In view of (9), it is easy to see that i gτ is independent of F gτ. This shows that the excursion straddling the hitting time τ is independent of the past of the Brownian motion up to time g τ (see also [23, Page 492, Lines -3]). Thus the first and second processes are independent. By the strong regenerative property of the Brownian motion, the third process is independent of the first two processes. This completes the proof of this lemma. We are now in a position to prove Theorem 2.2. Proof of Theorem 2.2. By the definition of φ, the processes {W t : t g τ } and {φ(w ) t : t g τ } are the same, and the processes {W τ+t W τ : t } and {φ(w ) τ+t φ(w ) τ : t } are also the same. By Lemma 2.4, in order to prove that φ(w ) is a Brownian motion, we only need to prove that the processes {W gτ +t : t τ g τ } and {φ(w ) gτ +t : t τ g τ } have the same distribution. We continue to use the notations in Lemma 2.3. For any t τ g τ, W τ t if W τ =, φ(w ) gτ +t = W τ t + if W τ =. 6

7 Thus the processes {φ(w ) gτ +t : t τ g τ } and {(ρ β t ) {ξ=} +( ρ β t ) {ξ=} : t β} have the same distribution. Let α be the hitting time of by ρ. It is a classical result [23, Chapter VII, Proposition 4.8] that the processes {ρ t : t α} and { ρ α t : t α} have the same distribution. This suggests that the processes {(ρ β t ) {ξ=} + ( ρ β t ) {ξ=} : t β} and { ρ t {ξ=} + ρ t {ξ=} : t β} have the same distribution. By Lemma 2.3, the processes {W gτ +t : t τ g τ } and { ρ t {ξ=} + ρ t {ξ=} : t β} have the same distribution. This shows that the processes {W gτ +t : t τ g τ } and {φ(w ) gτ +t : t τ g τ } have the same distribution. Remark 2.5. For any given x >, if we replace τ by the hitting time of {x, x} and change the definition of the quasi-time-reversal φ correspondingly, then the result in Theorem 2.2 still holds. 3 Cycle symmetry for diffusion processes on the circle In this section, we shall use the quasi-time-reversal invariance of the Brownian motion to prove the cycle symmetry for diffusion processes on the circle. Let X = {X t : t } be a one-dimensional diffusion process with diffusion coefficient a : R (, ) and drift b : R R. The initial distribution of X can be arbitrary. We assume that a and b are continuous periodic functions: a(x + ) = a(x), b(x + ) = b(x), x R. Due to the periodicity, a and b can be viewed as defined on S and X can be viewed as a diffusion process on S. In the sequel, we shall use the symbol X to denote both the diffusion process on S and its lifted process on R. The specific implication of X should be clear from the context. We shall construct the diffusion process X as the weak solution to the stochastic differential equation dx t = b(x t )dt + σ(x t )dw t, () where σ = a /2 and W is a standard Brownian motion defined on some filtered space (Ω, G, G t, P ). Since a and b are continuous periodic functions and a >, the Stroock-Varadhan uniqueness theorem [24, Theorem 7.2.] ensures that the martingale problem for (a, b) is well-posed. By the equivalence between the martingale-problem formulation and the weak-solution formulation [25, Theorem 2.], the weak solution to () always exists and is unique in law. Recall that the potential function U and scale function s of X are defined as U(x) = 2 x b(y) σ 2 dy, s(x) = (y) x e U(y) dy, respectively. In addition, recall that the distribution µ Y of any real-valued continuous process Y = {Y t : t } is a probability measure on W defined as µ Y (A) = P (Y A), A F. Definition 3.. The cycle forming time T of X is defined as T = inf{t : X t X = }. If X T X =, then we say that X forms the forward cycle at time T. If X T X =, then we say that X forms the backward cycle at time T. 7

8 Definition 3.2. For each n, the nth cycle forming time T n of X is defined as T n = inf{t T n : X t X Tn = }, where we assume that T =. If X Tn X Tn =, then we say that X forms the forward cycle at time T n. If X Tn X Tn =, then we say that X forms the backward cycle at time T n. Intuitively, T n is the time needed for X to form a cycle for the nth time and T = T. Under our assumptions, all these cycle forming times T n are finite a.s. [23, Chapter VII, Exercise 3.2]. Definition 3.3. Let ν + = inf{n : X forms the forward cycle at time T n } and let ν = inf{n : X forms the backward cycle at time T n }. Then the forming time T + of the forward cycle is defined as T + = T ν + and the forming time T of the backward cycle is defined as T = T ν. Intuitively, T + is the time needed for X to form a forward cycle for the first time and T is the time needed for X to form a backward cycle for the first time. It is clear that T = T + T. Theorem 3.4. Let X be a diffusion process solving the stochastic differential equation () with any initial distribution, where b : R R and σ : R (, ) are continuous functions with period. Then (i) for any u >, P (T + u, T + < T ) P (T u, T < T + ) = P (T + < T ( ) ) P (T < T + ) = exp b(y) 2 σ 2 (y) dy ; (ii) for any u, P (T + u T + < T ) = P (T u T < T + ). () Remark 3.5. The above theorem shows that although the distributions of T + and T may not be the same, their distributions, given that the corresponding cycle is formed earlier than its reversed cycle, are the same. The equality (), which characterizes the symmetry of the forming times of the forward and backward cycles for diffusion processes on the circle, will be named as the cycle symmetry in this paper. The proof of the cycle symmetry is divided into three steps. First, using the quasi-time-reversal invariance of the Brownian motion, we shall prove the result when σ and b is smooth. Second, using some transformation techniques, we shall prove the result when σ and b are both smooth. Third, using some approximation techniques, we shall prove the result when σ and b are both continuous. Lemma 3.6. Let X be a diffusion process solving the stochastic differential equation dx t = b(x t )dt + dw t, X =, (2) where b : R R is a C function with period. Then for any u and A F u, ( µ X (A) = exp [ u ]) U(w u ) + (b 2 (w t ) + b (w t ))dt dµ W (w), 2 where U is the potential function of X. A 8

9 Proof. Since b is a continuous periodic function, b(x t ) is a bounded continuous adapted process. By Novikov s condition, the process ( M t = exp t b(x s )dw s 2 t ) b(x s ) 2 ds is a martingale. Let Q be a probability measure defined as dq = M u dp. By Girsanov s theorem, {X t : t u} is a standard Brownian motion under Q. Thus for any bounded measurable function f on C[, u], E P f(w ) = E P f(x) dq ( dp = EP f(x) exp 2 u U (X t )dx t + 2 Since b is a C function, U is a C 2 function. By Itô s formula, we obtain that This shows that u du(x t ) = U (X t )dx t + 2 U (X t )dt = U (X t )dx t b (X t )dt. E P f(w ) = E P f(x) exp This implies the result of this lemma. ( [ U(X u ) + 2 u ]) (b 2 (X t ) + b (X t ))dt. ) b 2 (X t )dt. We are now in a position to prove the cycle symmetry when σ and b is smooth. Lemma 3.7. Let X be a diffusion process solving the stochastic differential equation (2), where b : R R is a C function with period. Then for any u >, P (T u, X T = ) P (T u, X T = ) = P (X ( T = ) P (X T = ) = exp 2 ) b(y)dy. Proof. Let τ and g τ be defined as in (6) and (7), respectively. For an arbitrarily fixed u >, let A and B be two subsets of W defined as A = {w W : τ(w) u, w τ = }, B = {w W : τ(w) u, w τ = }. (3) Let φ be the quasi-time-reversal on W. It follows from Theorem 2.2 that φ µ W = µ W, where φ µ W is the push-forward measure defined as φ µ W ( ) = µ W (φ ( )). Since φ is a one-to-one map from B onto A, it follows from Lemma 3.6 that ( µ X (A) = exp A 2 ( = exp 2 B [ U(w u ) + [ U(φ(w) u ) + u ]) p(w t )dt dφ µ W (w) ]) p(φ(w) t )dt dµ W (w). where p = b 2 + b is a continuous function with period. Since b is a C function with period, for any w B, wu U(φ(w) u ) = U(w u + 2) = 2 wu = 2 b(y)dy 4 u wu+ b(y)dy 2 b(y)dy 2 w u b(y)dy = U(w u ) + 2U(). 9 wu+2 w u+ b(y)dy

10 Since p is a continuous function with period, for any w B, u = p(φ(w) t )dt = gτ p(w t )dt + gτ τ p(w t )dt + g τ p(w t )dt + τ u The above calculations show that ( µ X (A) = exp [ U(w u ) + 2U() + 2 This shows that B This completes the proof of this lemma τ g τ p(w gτ +τ t + )dt + p(w t )dt = u u u τ p(w t )dt. P (T u, X T = ) P (T u, X T = ) = µ X(A) µ X (B) = e U(). p(w t + 2)dt ]) p(w t )dt dµ W (w) = e U() µ X (B). We are now in a position to prove the cycle symmetry when σ and b are both smooth. Lemma 3.8. Let X be a diffusion process solving the stochastic differential equation dx t = b(x t )dt + σ(x t )dw t, X =, (4) where b : R R and σ : R (, ) are smooth functions with period. Then for any u >, P (T u, X T = ) P (T u, X T = ) = P (X ( T = ) ) P (X T = ) = exp b(x) 2 σ 2 (x) dx. Proof. Let f be the function on R defined as f(x) = x σ(y) dy. It is easy to see that f is a smooth diffeomorphism on R. By Itô s formula, it is easy to check that f(x) = {f(x t ) : t } is a diffusion process solving the stochastic differential equation dy t = c(y t )dt + dw t, Y =, where c = ( f b + 2 f σ 2) f is a function satisfying c(f( + )) = c(f( )). Since σ is a smooth function with period, for any x R, f(x + ) f(x) = x+ x σ(y) dy = dy = f(). σ(y) This shows that c( + f()) = c( ) and thus c is a smooth function with period f(). Since f is a smooth diffeomorphism, the cycle forming time T of the diffusion process X is exactly the cycle forming time of the diffusion process f(x): By Lemma 3.7, we have P (T u, f(x T ) = f()) P (T u, f(x T ) = f()) = T = inf{t : f(x t ) = f()}. ( P (f(x T ) = f()) P (f(x T ) = f()) = exp 2 f() c(y)dy ).

11 Since σ is a smooth function with period, it is easy to check that f() = f( ) and f() Thus we obtain that c(y)dy = P (T u, f(x T ) = f()) P (T u, f(x T ) = f( )) = which gives the desired result. ( f (x)b(x) + ) 2 f (x)σ 2 (x) f b(x) (x)dx = σ 2 (x) dx. P (f(x ( T ) = f()) P (f(x T ) = f( )) = exp 2 ) b(x) σ 2 (x) dx, Lemma 3.9. Let X be a diffusion process solving the stochastic differential equation (4), where b : R R and σ : R (, ) are continuous functions with period. For an arbitrarily fixed u >, let A and B be two subsets of W as defined in (3). Then µ X ( A) = µ X ( B) =. Proof. Let Ā be the closure of A and let A be the interior of A. It is easy to check that Ā A {τ < u, w(τ) =, h >, such that t [τ, τ + h], w(t) }, A {τ < u, w(τ) =, h >, t [τ, τ + h], such that w(t) > }. The above two relations show that A {τ = u, w(τ) = } {τ < u, w(τ) =, h >, such that t [τ, τ + h], w(t) } {τ < u, w(τ) =, h >, such that t [τ, τ + h], w(t) }. (5) Let M be the diffusion process solving the stochastic differential equation dm t = σ(m t )dw t, M =. By Girsanov s theorem, µ X and µ M are equivalent on the measurable space (W, F u ). Since M is a strong Markov process, it follows from (5) that µ M ( A) P (M u = ) + P ( h >, such that t [, h], M t } + P ( h >, such that t [, h], M t }, where P x ( ) = P ( M = x). It is a classical result that for each u > and x R, the transition function P u (x, dy) of M has a density p(u, x, y) with respect to the Lebesgue measure [24, Lemma 9.2.2]. This shows that P (M u = ) = for each u >. Since M is a continuous martingale, there exists a Brownian motion B, such that M t = B [M,M]t [23, Chapter V, Theorem.6], where [M, M] t = t σ2 (M s )ds is a strictly increasing process. Thus P ( h >, such that t [, h], M t } = P ( h >, such that t [, h], B t } =. Similarly, we can prove that P ( h >, such that t [, h], M t } =. This shows that µ X ( A) = µ M ( A) =. Similarly, we can prove that µ X ( B) =. The following result can be found in [24, Theorem..4].

12 Lemma 3.. Let X be a diffusion process solving the stochastic differential equation (4), where b : R R and σ : R (, ) are bounded continuous functions. For each n, let X n be a diffusion process solving the stochastic differential equation dx t = b n (X t )dt + σ n (X t )dw t, X =, where b n : R R and σ n : R (, ) are bounded continuous functions. Assume that for any M >, the following two conditions hold: Then µ Xn sup sup n x M lim sup n x M ( σ 2 n (x) + b n (x) ) <, ( σ 2 n (x) σ 2 (x) + b n (x) b(x) ) =. µ X as n, where stands for weak convergence. We are now in a position to prove the cycle symmetry when σ and b are both continuous. Proof of Theorem 3.4. We only need to prove the result when X starts from a given point x R. Without loss of generality, we assume that x =. It is easy to see that (ii) is a direct corollary of (i). Thus we only need to prove (i). Since the subset of smooth periodic functions are dense in the space of continuous periodic functions, we can find a sequence of smooth functions b n with period such that b n converges to b uniformly. Similarly, we can find a sequence of smooth functions σ n > with period such that σ n converges to σ uniformly. Let X n be the diffusion process solving the stochastic differential equation dx t = b n (X t )dt + σ n (X t )dw t, X =. For any u >, let A and B be two subsets of W as defined in (3). By Lemma 3.8, we have ( µ Xn (A) ) µ Xn (B) = exp b n (x) 2 σn(x) 2 dx. (6) Since b, σ, b n, σ n are all continuous periodic functions, it is easy to check that the conditions in Lemma 3. are satisfied. Thus µ Xn µ X. Moreover, it follows from Lemma 3.9 that µ X ( A) = µ X ( B) =. This shows that µ Xn (A) µ X (A) and µ Xn (B) µ X (B). Then we can obtain the desired result by letting n in (6). 4 Renewal random walk representation for diffusion processes on the circle In this section, we shall use the cycle symmetry to prove that the cycle dynamics of any diffusion process on the circle with any initial distribution is nothing but a renewal random walk. Here a renewal random walk stands for a nearest-neighbor random walk whose interarrival times are i.i.d. positive random variables. Proposition 4.. T and X T X are independent. Proof. It is easy to see that T + < T is equivalent to X T X = and T < T + is equivalent to X T X =. By Theorem 3.4, we obtain that P (T u X T X = ) = P (T u X T X = ). This implies that T and X T X are independent. 2

13 Remark 4.2. Intuitively, T is the time needed for X to form a forward or backward cycle for the first time and X T X characterizes which one of these two cycles is formed. Thus the above corollary shows that the forming time of a forward or backward cycle for a diffusion process on the circle is independent of which one of these two cycles is formed. The following result is interesting in its own right. Theorem 4.3. The distributions of T and X T X are independent of the initial distribution of X. Proof. Let µ be the initial distribution of X. It follows from Theorem 3.4 that ( P µ (X T X = ) ) P µ (X T X = ) = exp b(y) 2 σ 2 (y) dy. This shows that the distribution of X T X is independent of the initial distribution of X. We shall now view X as a diffusion process on S. Without loss of generality, we assume that X is the coordinate process on W (S ). Let T 2 be the second cycle forming time of X and let S = T 2 T. By the strong Markov property of X, for any t, P µ (S t G T ) = P (T θ T t G T ) = P XT (T t) = P µ (T t), (7) where θ T be the shift operator on W (S ) and the last equality holds because X T and X represent the same point on S. This shows that T and S are independent and have the same distribution. Since (7) holds for any initial distribution µ, we have for any x S, P µ (S t, T t) = P µ (T t) 2, P x (S t, T t) = P x (T t) 2. The above two equations suggest that [ ] 2 P x (T t)µ(dx) = P x (T t) 2 µ(dx). S S Since the above equation holds for any distribution µ, the function g(x) = P x (T t) must be a constant. This shows that the distribution of T is independent of the initial distribution of X. The above results suggest that the cycle dynamics of any diffusion process on the circle with any initial distribution is simply a renewal random walk: we only need to wait a random time for the process to form a cycle, toss a coin to decide whether the forward or the backward cycle is formed, and repeat the above procedures independently. This fact is stated rigorously in the following theorem. Theorem 4.4. For any t and n, let N t = inf{n : T n+ > t} and let ξ n = X Tn X Tn, where T =. Let L = {L t : t } be a process defined as L t = N t n= ξ n. Then L is a renewal random walk. Proof. For any n, let τ n = T n+ T n. By the strong Markov property of X and Proposition 4., for any x R, t,, t m, and s,, s m {, }, m m P x (τ t,, τ m t m, ξ = s,, ξ m = s m ) = P x (T t i ) P x (X T X = s i ). By Theorem 4.3, it is easy to see that m m P (τ t,, τ m t m, ξ = s,, ξ m = s m ) = P (T t i ) P (X T X = s i ). This shows that {τ n : n } are i.i.d. positive random variables and {ξ n : n } are i.i.d. Bernoulli random variables independent of {τ n : n }. Thus L is a renewal random walk. 3 i= i= i= i=

14 5 Limit theorems and large deviations for sample circulations In this section, we shall use the renewal random walk representation for diffusion processes on the circle to study the limit theorems and large deviations for the sample circulations. 5. Some results on Markov renewal processes To make the paper more self-contained, we recall the following definition. Definition 5.. Let {ξ n : n } be an irreducible discrete-time Markov chain with finite state space E. Assume that each i E is associated with a Borel probability measure φ i on (, ). Let {τ n : n } be a sequence of positive random variables such that given {ξ n : n }, the random variables {τ n : n } are independent and have the distribution P (τ n ξ,, ξ n ) = φ ξn ( ). Then {(ξ n, τ n ) : n } is called a Markov renewal process. Remark 5.2. It is easy to see that a renewal random walk can be represented as a Markov renewal process, where {τ n : n } are i.i.d positive random variables and {ξ n : n } are i.i.d. Bernoulli random variables independent of {τ n : n }. In fact, the limit theorems and large deviations for Markov renewal processes are well-established. Before we state these results, we introduce several notations. Let {(ξ n, τ n ) : n } be a Markov renewal process, where ξ = {ξ n : n } is an irreducible discrete-time Markov chain on a finite state space E = {,, N} with transition probability matrix P = (p ij ) and invariant distribution π = (π,, π N ). For any n and t, let T n = n k= τ k be the nth jump time of the Markov renewal process and let N t = inf{n : T n+ > t} be the number of jumps up to time t. Let r = (r ij ) be a given matrix and define the reward process R = {R t : t } by N t R t = r ξn ξ n. (8) n= Let g = (g,, g N ) be a column vector defined as g i = j E p ij r ij, and let f be a solution to the following Poisson equation: (I P )f = g πg, where is a column vector whose components are all. Since P is irreducible, it is easy to check that the rank of I P is N. Thus the solution to the Poisson equation is unique up to an additive constant. For any i, j E, set m i = tφ i (dt), v i = Moreover, set m = π i m i, i E [, ) [, ) t 2 φ i (dt), H ij = r ij + f j f i. u = m π i p ij r ij, σ 2 = m π i p ij (Hij 2 2um i H ij + u 2 v i ). i,j E i,j E 4

15 The following lemma, which gives the strong law of large numbers and functional central limit theorem for Markov renewal processes, is due to Glynn and Hass [26]. Lemma 5.3. Assume that i E m i <. Then R t lim t t = u, a.s. For any η, let U η = {U η t : t } be the process defined as U η t = η /2 (R ηt uηt). If v i < for any i E, then U η σw on D[, ) as η, where D[, ) is the Skorokhod space and W is a standard Brownian motion. The following lemma, which gives the large deviations for Markov renewal processes, is due to Mariani and Zambotti [27]. Lemma 5.4. For any t, let Q t C(E E, [, )) be the empirical flow defined as Q t (i, j) = t N t n= {ξn =i,ξ n=j}. (9) Then the law of Q t satisfies a large deviation principle with rate t and good rate function I : C(E E, [, )) [, ]. Moreover, the rate function I is convex. The explicit expression of the rate function I is rather complicated. Readers who are interested in this expression may refer to Equations (5)-(7) in [27]. 5.2 Strong law of large numbers for sample circulations Let X be a diffusion process solving the stochastic differential equation (), where b : R R and σ : R (, ) are continuous functions with period. For further references, we introduce some notations. For any t, let N t = inf{n : T n+ > t} be the number of cycles formed by X up to time t and let N t N t N t + = {XTn X Tn =}, Nt = n= n= {XTn X Tn = } be the numbers of the forward and backward cycles formed by X up to time t, respectively. Definition 5.5. The sample circulations J t + and Jt along the forward and backward cycles up to time t are defined as J t + = t N t +, J t = t N t, (2) respectively. The sample net circulation J t of X up to time t is defined as J t = J t + Jt. The following theorem gives the strong law of large numbers for the sample circulations. Recall that if we regard X as a diffusion process on S, then X is always ergodic with respect to its unique invariant distribution [28]. 5

16 Theorem 5.6. Let ρ be the invariant distribution of X. Then for any initial distribution of X, lim J + t t = J +, lim J t t = J, a.s., (2) where J + and J are two positive constants satisfying J + ( ) J = exp b(x) 2 σ 2 (x) dx, J + J = b(x)ρ(dx), S with b viewed as a function on S in the second equation. Proof. By Theorem 4.4, the cycle dynamics of X is a renewal random walk and thus N t is a renewal process. By the strong law of large numbers and elementary renewal theorem, we obtain that J + = lim t t N t + By Theorem 3.4, we have N t = lim t t N t N t n= J + J = P (X T X = ) P (X T X = ) = P (T + < T ) P (T < T + ) = exp {XTn X Tn =} = P (X T X = ), a.s. (22) ET ( 2 Since X is the solution to the stochastic differential equation (), we have X t t = X t + t t b(x s )ds + t t ) b(x) σ 2 (x) dx. σ(x s )dw s. (23) It follows from Birkhoff s ergodic theorem that t lim b(x s )ds = b(x)ρ(dx), t t S a.s. (24) Let M t = t σ(x s)dw s. Since M is a continuous martingale, there exists a Brownian motion B such that M t = B [M,M]t. By Birkhoff s ergodic theorem again, we obtain that [M, M] t lim t t = lim t t By Khinchin s law of the iterated logarithm, it is easy to check that t σ 2 (X s )ds = σ 2 (x)ρ(dx) >, S a.s. (25) B t lim t t This fact, together with (24) and (25), shows that lim t t t σ(x s )dw s = lim t B [M,M]t t =, a.s. B [M,M]t = lim [M, M] t t [M, M] t t =, a.s. (26) In addition, it is easy to see that N t + Nt = X TNt X and X TNt X t. This implies that J t X t t t. (27) Combining (23), (24), and (26), we finally obtain that X t J = lim J t = lim = b(x)ρ(dx), t t t S a.s.. (28) This completes the proof of this theorem. 6

17 Definition 5.7. The limits J + and J in (2) are called the circulations of X along the forward and backward cycles, respectively. The net circulation J of X is defined as J = J + J. Intuitively, J + and J represent the numbers of the forward and backward cycles formed by X per unit time, respectively, and J represents the net number of cycles formed by X per unit time. Due to this reason, the net circulation J is also called the rotation number [28]. The following definition originates from nonequilibrium statistical physics. Definition 5.8. The affinity γ of X is defined as γ = log J + J = 2 b(y) σ 2 (y) dy. 5.3 Functional central limit theorem for sample circulations Lemma 5.9. For any p >, ET p <. Proof. Without loss of generality, we assume that X starts from. It is easy to see that the scale function s of X is an strictly increasing C 2 function with s() =. By Itô s formula, it is easy to check that M = s(x) is a diffusion process solving the stochastic differential equation dm t = f(m t )dw t, Y =. where f = (s σ) s. Let σ z be the hitting time of z by s(x). Then T = σ s( ) σ s(). Since the distribution of T only depends on the value of f on [s( ), s()], we can assume that f a for some constant a >. Since M is a continuous local martingale, there exists a Brownian motion B such that M t = B [M,M]t, where [M, M] t = t f 2 (M s )ds a 2 t. Let τ z be the hitting time of z by B. It is easy to check that T τ s( ) τ s() /a 2. This shows that This completes the proof of this lemma. ET p a 2p E[τ s( ) τ s() ] p <. The functional central limit theorem for the sample circulations is stated in the next theorem. Theorem 5.. For any η, let J η,+ = {J η,+ t } be three processes defined as : t }, J η, = {J η,+ t : t }, J η = {J η t : t J η,+ t = η /2 (N ηt + J + ηt), J η, t = η /2 (Nηt J ηt), J η t = J η,+ t J η, t, respectively. Then (J η,+, J η,, J η ) ΣW on D([, ), R 3 ) as η, where D([, ), R 3 ) is the Skorokhod space, W is a standard three-dimensional Brownian motion, and ρ + + αρ 2 + 2αρ ρ + 2ρ + + 2αρ + (ρ + ρ ) ΣΣ = m 2αρ ρ + ρ + αρ 2 2ρ + 2αρ (ρ + ρ ) 2ρ + + 2αρ + (ρ + ρ ) 2ρ + 2αρ (ρ + ρ ) + α(ρ + ρ ) 2 with m = ET, α = ET 2 /(ET ) 2 2, ρ + = P (T + < T ), and ρ = P (T < T + ). 7

18 Proof. Let K = (J η,+, J η,, J η ). By [29, Chapter 29, Theorem 29.6], the desired result holds if and only if λ K λ ΣW for any row vector λ = (a, b, c) R 3. Note that λ K t = aj η,+ t where u = aj + + bj + cj and + bj η, t + cj η t = η /2 (R ηt uηt), N t R t = an t + + bnt + c(n t + Nt ) = (a + c) {ξn=} + (b c) {ξn= } is a reward process with r, = r, = b c and r, = r, = a + c. By Lemma 5.9, it is easy to see that m = m = ET < and v = v = ET 2 <. By Lemma 5.3, it is straightforward to check that λ K λ ΣW. 5.4 Large deviations for sample circulations n= The large deviation principle for the sample circulations is stated in the following theorem. For the precise definition of the large deviation principle for a family of probability measures, please refer to [3, Page 3]. Theorem 5.. The law of (J + t, J t, J t) satisfies a large deviation principle with rate t and good rate function I : R 3 [, ]. Moreover, I is convex and I (J +, J, J) =. Proof. By Theorem 4.4, L is a renewal random walk and thus can be viewed as a Markov renewal process. Let {Q t : t } be the empirical flow of L defined in (9). Then J + t N t N t = t {ξn=} = Q t (i, ), Jt = t {ξn= } = Q t (i, ). n= i E n= i E We define a continuous map F : C(E E, [, )) R 3 as ( ) F (Q) = (Q(i, ) Q(i, )). (29) i E Q(i, ), i E Q(i, ), i E Then we have (J + t, J t, J t) = F (Q t ). By Lemma 5.4, the law of Q t satisfies a large deviation principle with rate t and good rate function I. By the contraction principle, the law of (J + t, J t, J t) satisfies a large deviation principle with rate t and good rate function I : R 3 [, ] defined as I (x) = inf I(Q). (3) Q F (x) We next prove that I is convex. We make a crucial observation that the map F is linear. Moreover, it follows from Lemma 5.4 that I is convex. These two facts suggest that for any λ, λ 2, λ 3 satisfying λ + λ 2 + λ 3 = and x, x 2, x 3 R 3, I (λ x + λ 2 x 2 + λ 3 x 3 ) inf I(Q) Q λ F (x )+λ 2F (x 2)+λ 3F (x 3) = inf Q F (x ),Q 2 F (x 2),Q 3 F (x 3) I(λ Q + λ 2 Q 2 + λ 3 Q 3 ) = λ inf I(Q) + λ 2 inf I(Q) + λ 3 inf I(Q) Q F (x ) Q F (x 2) Q F (x 3) = λ I (x ) + λ 2 I (x 2 ) + λ 3 I (x 3 ), 8

19 which shows that I is also convex. We finally prove that I (J +, J, J) =. Let B ɛ = {y R 3 : y (J +, J, J) ɛ}. Since I is lower semi-continuous and (J + t, J t, J t) (J +, J, J), a.s., we obtain that I (J +, J, J) lim ɛ inf y B ɛ I (y) lim sup This clearly shows that I (J +, J, J) =. ɛ lim sup t t log P ((J + t, J t, J t) B ɛ ) =. Since the cycle dynamics of X is simply a renewal random walk, where the interarrival times and jump probability are both independent of the initial distribution of X. This suggests that the rate function I is also independent of the initial distribution of X. 6 Fluctuation theorems for diffusion processes on the circle In this section, we shall use the cycle symmetry to prove various types of FTs for diffusion processes on the circle. 6. Fluctuation theorems for sample circulations The next proposition characterizes the symmetry of the distribution of the sample circulations. Results of the following type are called transient FTs in nonequilibrium statistical physics. Proposition 6.. Let N = {,, 2, }. For each t and any m, n N, where γ is the affinity of X. P (N + t = n, N t = m) P (N + t = m, N t = n) = eγ(n m), Proof. With the notations in Theorem 4.4, we have P (N + t = n, N t = m) = s,,s n+m P (T n+m t < T n+m+, ξ = s,, ξ n+m = s n+m ), where s,, s n+m ranges over all choices such that n of them equal and m of them equal. By Theorem 4.4, it is easy to see that P (T n+m t < T n+m+, ξ = s,, ξ n+m = s n+m ) = P (T n+m t < T n+m+ )ρ n +ρ m, where ρ + = P (T + < T ) and ρ = P (T < T + ). By Theorem 3.4, we obtain that P (N + t = n, N t = m) = s,,s n+m P (T n+m t < T n+m+ )ρ n +ρ m = P (T n+m t < T n+m+ )C n n+mρ m + ρ n e γ(n m) = P (N + t = m, N t = n)e γ(n m), which gives the desired result. The following lemma gives a lower bound for the survival function of the cycle forming time T. In order not to interrupt things, we defer the proof of this lemma to the final section of this paper. 9

20 Lemma 6.2. Let T be the cycle forming time of X. Then there exists ρ >, such that for any t, P (T > t) e ρt. The next result follows directly from the above lemma. Lemma 6.3. There exists ρ >, such that for any β and t, Ee βnt e (eβ )ρt. Proof. By Lemma 6.2, there exists ρ > such that P (T > t) e ρt. This shows that the cycle forming time T is stochastically larger than an exponential random variable with rate ρ. Thus it is easy to see that the nth cycle forming time T n is stochastically larger than the sum of n independent exponential random variables with rate ρ. This further suggests that N t is stochastically dominated by a Poisson random variable R t with parameter ρt. Thus we obtain that Ee βnt = which gives the desired result. n= βe βx P (N t x)dx βe βx P (R t x)dx ( ) = Ee βrt βn (ρt)n = e e ρt = exp (e β )ρt, n! The next proposition characterizes the symmetry of the moment generating function of the sample circulations. nonequilibrium statistical physics. Proposition 6.4. Let Results of the following type are called Kurchan-Lebowitz-Spohn-type FTs in g t (λ, λ 2 ) = Ee λn + t +λ 2N t = Ee t(λj + t +λ 2J t ). Then for each t and any λ, λ 2 R, we have g t (λ, λ 2 ) < and Proof. By Lemma 6.3, we obtain that g t (λ, λ 2 ) = g t (λ 2 γ, λ + γ). g t (λ, λ 2 ) Ee λ N + t + λ 2 N t Ee α(n + t +N t ) = Ee αnt <, (3) where α = max{ λ, λ 2 }. Moreover, it follows from Proposition 6. that g t (λ, λ 2 ) = e λn+λ2m P (N t + = n, Nt = m) = n,m N n,m N which gives the desired result. e (λ+γ)n+(λ2 γ)m P (N + t = m, N t = n) = g t (λ 2 γ, λ + γ), The next theorem characterizes the symmetry of the rate function of the sample circulations. Theorems of the following type are called Gallavotti-Cohen-type FTs in nonequilibrium statistical physics. Recall that the Legendre-Fenchel transform of a function f : R n [, ] is defined as f (λ) = sup x R n {λ x f(x)}. and a function f : R n (, ] is called proper if f(x) < for at least one x. 2

21 Theorem 6.5. The law of (J t +, J t ) satisfies a large deviation principle with rate t and good rate function I : R 2 [, ], which is convex and satisfies I (J +, J ) =. Moreover, the rate function I has the following symmetry: for any x, x 2 R, I (x, x 2 ) = I (x 2, x ) γ(x x 2 ). Proof. By Theorem 5. and the contraction principle, the law of (J t +, J t ) satisfies a large deviation principle with rate t and good rate function I : R 2 [, ]. The proof of the facts that I is convex and I (J +, J ) = follows the same line as that of Theorem 5.. By a strengthened version of Varadhan s lemma [3, Theorem 4.3.], if there exists β > such that the following moment condition is satisfied: then for any λ, λ 2 R, By Lemma 6.3, we obtain that lim sup t + t +λ log Eeβt(λJ 2J t ) <, t lim t t log g t(λ, λ 2 ) = I(x, x 2 ). (32) Ee βt(λj + t +λ 2J t ) Ee β λ N + t +β λ 2 N t Ee βαnt exp where α = max{ λ, λ 2 }. This shows that for any β, lim sup t ( ) (e βα )ρt, + t +λ log Eeβt(λJ 2J t ) lim sup t t t log (e βα )ρ <. EeβαNt Thus we have proved (32). It thus follows from (32) and Proposition 6.4 that I (λ, λ 2 ) = I (λ 2 γ, λ + γ). This implies that I (x, x 2 ) = sup {λ x + λ 2 x 2 I(λ 2 γ, λ + γ)} λ,λ 2 R = sup {(λ γ)x + (λ 2 + γ)x 2 I(λ 2, λ )} = I (x 2, x ) γ(x x 2 ). λ,λ 2 R (33) On the other hand, the Fenchel-Moreau theorem [32, Theorem 4.2.] shows that if a function f is proper, then f = f if and only if f is convex and lower semi-continuous. By Theorem 5., I is a good rate function which is also convex. This shows that I is proper, convex, and lower semicontinuous. It thus follows from the Fenchel-Moreau theorem that I = I. This fact, together with (33), gives the desired result. 6.2 Fluctuation theorems for sample net circulation The following transient FT characterizes the symmetry of the distribution of the sample net circulation. Proposition 6.6. For each integer k, we have P (N t = k) P (N t = k) = eγk. (34) 2

22 Proof. By Proposition 6., we have P (N t = k) = P (N t + = n, Nt = m) n m=k = P (N t + = m, Nt = n)e γ(n m) = P (N t = k)e γk, n m=k which gives the desired result. The following Kurchan-Lebowitz-Spohn-type FT characterizes the symmetry of the moment generating function of the sample net circulation. Proposition 6.7. For each t and λ R, we have Ee λnt < and Ee λnt = Ee (λ+γ)nt. Proof. If we take λ = λ and λ 2 = λ in Proposition 6.4, then we have Ee λnt < and which is just the desired result. Ee λnt = g t (λ, λ) = g t ( (λ + γ), λ + γ) = Ee (λ+γ)nt, Results of the following type are called integral FTs in nonequilibrium statistical physics. Corollary 6.8. For each t, Ee γnt =. Proof. If we take λ = γ in Proposition 6.7, then we obtain the desired result. The following Gallavotti-Cohen-type FT characterizes the symmetry of the rate function of the sample net circulation. Theorem 6.9. The law of J t satisfies a large deviation principle with rate t and good rate function I 2 : R [, ], which is convex and satisfies I 2 (J) =. Moreover, the rate function I 2 has the following symmetry: for each x R, I 2 (x) = I 2 ( x) γx. Proof. By Theorem 6.5 and the contraction principle, the law of J t satisfies a large deviation principle with rate t and good rate function I 2 (x) = inf I (x, x 2 ). x x 2=x The proof of the facts that I 2 is convex and I 2 (J) = follows the same line as that of Theorem 5.. By Theorem 6.5, we have I 2 (x) = inf [I (x 2, x ) γ(x x 2 )] = inf I (x 2, x ) γx = I 2 ( x) γx. x x 2=x x 2 x = x This completes the proof of this theorem. Remark 6.. It is easy to see that the FTs for the sample circulations and net circulation can be easily extended to general renewal random walks whenever the distribution of the interarrival times satisfies the inequality in Lemma

23 6.3 Relationship with reversibility The FTs discussed above are closely related to the reversibility for diffusion processes on the circle. The following theorem gives several equivalent conditions for X to be a symmetric Markov process, where X should be viewed as a diffusion process on S rather than on R. Theorem 6.. Let ρ be the invariant distribution of X. Then the following statements are equivalent: (i) X is a symmetric Markov process with respect to ρ; (ii) γ = ; (iii) J = ; (iv) I (x, x 2 ) = I (x 2, x ) for any x, x 2 R; (v) I 2 (x) = I 2 ( x) for any x R. Proof. It is a classical result that (i) and (ii) are equivalent [, 28]. By the definition of the affinity γ, it is easy to see that γ = if and only if J + = J. This shows that (ii) and (iii) are equivalent. By Theorem 6.5 and Theorem 6.9, it is easy to see that (ii),(iv), and (v) are equivalent. 6.4 Fluctuation theorems for sample entropy production rate Entropy production rate is a central concept in nonequilibrium statistical physics. In the previous work, the transient and integral FTs for the sample entropy production rate have been established for pretty general stochastic processes [33]. However, it turns out to be rather difficult to establish the large deviation principle and Gallavotti-Cohen-type FT for the sample entropy production rate. Here, we shall give a simple proof of the Gallavotti-Cohen-type FT for the sample entropy production rate of diffusion processes on the circle. Let X = {X t : t } be a diffusion process on S with smooth diffusion coefficient a : S (, ) and smooth drift b : S R. In fact, X is a Brownian motion with drift on the Riemannian manifold S with some changed Riemannian metric [, Chapter 5, Page 2]. Thus X has a unique invariant distribution ρ, which has a strictly positive smooth density η(x) with respect to the volume element of S [34, Chapter V, Proposition 4.5]. Definition 6.2. Let Y be the lifted process of X on R with diffusion coefficient a : R (, ) and drift b : R R. The sample entropy production rate E t of X up to time t is defined as E t = t [ ] 2b t a (log(aη)) (Y s ) dy s, where η is viewed as a periodic function on R and the symbol means that the stochastic integral is taken in the Stratonovich sense. Remark 6.3. If X is a stationary diffusion process on S, there is still another way to defined the sample entropy production rate E t. Let µ + t denote the distribution of the process {X s : s t} and let µ t denote the distribution of its timed-reversal {X t s : s t}. Let C([, t], S ) denote the space of all continuous functions on [, t] with values in S. Then µ + t and µ t are both probability measures on C([, t], S ). With these notations, the sample entropy production rate E t 23