Entropy prouction of iffusion in spatially perioic eterministic systems arxiv:nlin/0203046v [nlin.cd] 2 Mar 2002 J. R. Dorfman, P. Gaspar, 2 an T. Gilbert 3 Department of Physics an Institute for Physical Science an Technology, University of Marylan, College Park, MD 20742, USA 2 Center for Nonlinear Phenomena an Complex Systems, Université Libre e ruxelles, Coe Postal 23, Campus Plaine, -050 russels, elgium 3 Department of Chemical Physics, The Weizmann Institute of Science, Rehovot 7600, Israel (Date: December 30, 203) This paper presents an ab initio erivation of the expression given by irreversible thermoynamics for the rate of entropy prouction for ifferent classes of iffusive processes. The first class are Lorentz gases, where non-interacting particles move on a spatially perioic lattice, an collie elastically with fixe scatterers. The secon class are perioic systems where N particles interact with each other, an one of them is a tracer particle which iffuses among the cells of the lattice. We assume that, in either case, the ynamics of the system is eterministic an hyperbolic, with positive Lyapunov exponents. This work extens methos originally evelope for a chaotic two-imensional moel of iffusion, the multi-baker map, to higher imensional, continuous time ynamical systems appropriate for systems with one or more moving particles. Here we express the rate of entropy prouction in terms of hyroynamic measures that are etermine by the fractal properties of microscopic hyroynamic moes that escribe the slowest ecay of the system to an equilibrium state. PACS numbers: 05.45.Df, 05.60.-k,05.70.Ln I. INTRODUCTION In 902, Gibbs escribe a mechanism by which the entropy coul increase towar its equilibrium value in mechanical systems obeying Newton s equations []. Gibbs mechanism is base on the assumption that the microscopic ynamics is mixing. The mixing woul allow coarse-graine probabilities to reach their equilibrium values after a long time, a result that has receive a rigorous meaning in the moern efinition of mixing [2]. The secon ingreient of Gibbs mechanism is the assumption that the entropy of a physical system shoul be efine as a quantity which is now referre to as the coarse-graine entropy. The use of the coarse-graine entropy coul be ustifie by the fact that, if the entropy shoul be given accoring to oltzmann by the logarithm of the number of complexion of a system, then it can only be efine by introucing cells of non-vanishing size in systems escribe by continuous coorinates.[20] The aim of the present paper is to apply the program set up by Gibbs to hyperbolic, eterministic ynamical systems sustaining a transport process of iffusion. We suppose that the systems obey Liouville s theorem, namely that phase space volumes are preserve by the ynamics, which is a maor assumption use by Gibbs. Specifically, the systems we consier are either perioic Lorentz gases, in which moving particles iffuse through a lattice interacting only with fixe scatterers, or are pe- Electronic aress: r@ipst.um.eu Electronic aress: gaspar@ulb.ac.be Electronic aress: thomas.gilbert@weizmann.ac.il rioic repetitions of interacting N-particle systems such that a tagge particle is followe as it unergoes iffusion among the unit cells. The present work is an extension to continuous time, an interacting systems, of previous works [3, 4, 5, 6, 7] concerne with the multi-baker map, as a chaotic moel of iffusion. In the previous work, we use the fact that an initial non-equilibrium istribution function, in the multi-baker phase space, rapily evelops a fractal structure ue to the chaotic nature of the ynamics. This structure is such that variations of the istribution function on arbitrarily fine scales evelop as the system evolves in time. The final stages of the approach to equilibrium are then controlle by the ecay of fractal, microscopic hyroynamic moes of the system, in this case, iffusive moes, which ecay with time as exp( Dk 2 t), where k is a wave number characterizing a particular moe associate with a macroscopic ensity variation over a istance of orer k, D is the iffusion coefficient, an t is the time. For the multi-baker system it is possible to express the rate of entropy prouction in this final stage in terms of measures of sets which are etermine by the non-equilibrium phase space istribution in the set, in particular, by the values of the fractal hyroynamic moes in the set. In this paper we show that it is possible to apply the same methos to calculate the rate of entropy prouction for iffusive flows in perioic Lorentz gases an for tracer iffusion in perioic, interacting N-particle systems, as long as the microscopic ynamics is eterministic, mixing, an chaotic. Our metho is base on the explicit construction of the microscopic fractal hyroynamic moes of iffusion, which characterize the long time relaxation of the system towar thermoynamic equilibrium. Our main result is that we obtain by this
2 metho exactly the expression for the rate of entropy prouction as given by irreversible thermoynamics for these systems[8]. The source of this agreement can be trace to the role playe by the fractal hyroynamic moes, both for requiring a coarse graining of the phase space to properly incorporate the effects of their fractal properties on entropy prouction in the system, as well as for escribing the slowest ecay of the system as it relaxes to equilibrium. The plan of the paper is the following. The mathematical methos neee to escribe spatially perioic systems are provie in Sec. II. Then the non-equilibrium istribution is efine in Sec. III. Once the istribution function has been constructe, we can ientify the microscopic hyroynamic moes of iffusion, an they are then constructe in Sec. IV. Here we use the properties of these moes to introuce the notion of a hyroynamic measure of a set in phase space, an we also ientify the sets that are use in the coarse graining of the phase space. These sets are not arbitrary, but must have some specific properties in orer to be useful for the calculation of the rate of entropy prouction, which is carrie out in Sec. V. We conclue with a iscussion of the metho an results obtaine here, an with an outline of irections for future work in Sec. VI. II. SPATIALLY PERIODIC SYSTEMS A. The phase space We consier a eterministic ynamical system of phase space imension M which is spatially perioic in the form of a -imensional lattice L. We will label the positions of the perioic cells M on the lattice by the vector l L, an the phase space coorinates within an elementary cell by X M. The lattice L is isomorphic to Z. The total phase space of the system is the irect prouct M L of imension M. The time isplacement operator over a time interval t, acting on points ( l, X), is enote Φ t, which is also calle the flow. On the other han, the time isplacement operator acting on points X insie the basic unit cell is enote φ t. Examples of such systems are the following.. The two-imensional Lorentz gas In this system a point particle moves in free flight an unergoes elastic collisions on har isks forming a triangular lattice. The phase space is efine by the spatial an velocity coorinates ( r, v) of the moving particle. The spatial coorinates vary over the plane with the exclusion of the area occupie by the isks: r Q with im Q = 2. The velocity coorinates form another twoimensional vector v R 2. Energy is conserve uring the motion so that each energy shell is preserve by the ynamics. In each energy shell, we can thus reuce the coorinates to the two positions r = (x, y) Q an the velocity angle ϕ [0, 2π[. The position space ecomposes into a triangular lattice of hexagonal elementary cells, each containing a single isk: Q = l L T l C, where T l enotes the translation by the lattice vector l, an C is the elementary cell of the position space Q. The elementary cell of the phase space has thus the coorinates X = (x, y, ϕ) M = C [0, 2π[. In the twoimensional Lorentz gas, the phase space has the imension im M = 3 while the lattice has the imension = im L = 2. The flow of the har-isk Lorentz gas preserves the Lebesgue measure X = x y ϕ. We notice that the horizon of the har-isk Lorentz gas must be finite in orer for the iffusion coefficient to be finite [9], which we assume in the following. Similar consierations apply to the Lorentz gas in which a point particle moves in a perioic lattice of attractive Yukawa potentials. In this system the iffusion coefficient is positive an finite for a large enough energy [0]. Here, too, the phase space has 3 imensions an is the union of the constant energy surfaces for each cell of the lattice. 2. The three-imensional Lorentz gas This Lorentz gas is the irect generalization of the twoimensional one. The spatial coorinates vary over the space with the exclusion of the volume occupie by the spheres: r Q with im Q = 3. The velocity coorinates form another three-imensional vector v R 3. In each energy shell, we can reuce the coorinates to the three positions r = (x, y, z) Q an the two velocity spherical angles θ [0, π] an ϕ [0, 2π[. The position space ecomposes into a lattice of elementary cells containing one or several isks: Q = l L T l C. The elementary cell of the phase space has thus the coorinates X = (x, y, z, cosθ, ϕ) M = C [, +] [0, 2π[, where, again, C is the elementary cell of the position space Q. In the three-imensional Lorentz gas, the phase space has the imension im M = 5 while the lattice has the imension = im L = 3. The flow of the har-isk Lorentz gas preserves the Lebesgue measure X = xy z cosθ ϕ. Here also, we suppose that the horizon is finite in orer for the iffusion coefficient to be finite. 3. Diffusion of a tracer in a system on a torus The molecular ynamics simulation of the iffusion of a tracer particle moving in a flui can be performe by consiering a finite number of particles moeling the flui an the tracer particle, all of them moving with interactions, in a omain elimite by perioic bounary conitions. The total number of particles is equal to N. The center of mass can be taken at rest. The vector l
3 can be use to locate the position of the cell containing the tracer particle as it moves on the checkerboar lattice mae of infinitely many images of the system, which tiles the -imensional space of the system. Then the iffusion coefficient of the tracer particle can be compute by aing the appropriate lattice vector to l each time the tracer particle crosses a bounary. The energy an total momentum are to be conserve so that here X enotes the phase space coorinates of an elementary cell of the phase space, which is of imension M = im M = 2N 2 after elimination of the egrees of freeom of the center of mass, while the lattice is of imension = im L = 2 or 3. In summary, we will enote by X the coorinates on an energy-momentum shell of a micro-canonical ensemble for the perioic system. We suppose that the flow Φ t preserves the Lebesgue measure X. Moreover, we assume that the iffusion coefficient of the system is finite, that the mean rift vanishes, an that the microscopic ynamics is chaotic an mixing. In the sequel l will enote the lattice position vector an T the translation operator on the lattice. Subsets of the unit cell M will be enote by capital Roman letters. The notation A will be use for a set belonging to the elementary phase space cell at the origin l = 0 of the lattice. The notation A l = T l A will be use when we want to refer to a set in a specific lattice cell at position l.. Lattice Fourier transforms We will also nee to efine lattice Fourier transforms []. We will nee the preliminary result : Definition Consier a function G( l, X), which is a function of the lattice coorinate l, an the unit cell coorinate X, then this function can be expresse in terms of a lattice Fourier transform as G( l, X) = k e i k l G( k, X), () where enotes the first rillouin zone of the reciprocal lattice, is its volume, an G is the lattice Fourier transform of G. It is important to note that : Lemma. If G is only a function of the lattice vector l, then G oes not epen upon the unit cell coorinate X, but only upon k. III. THE NON-EQUILIRIUM DISTRIUTION We consier a perioic eterministic ynamical system with finite iffusion coefficient. A. The non-equilibrium measure We now construct the statistical ensemble we will use for the rest of this paper. We suppose that the coorinates of the ensemble are istribute on the lattice in such a way that the istribution can be escribe by an initial ensemble ensity ρ( l, X, 0), where ρ enotes the number of systems per unit phase space volume. We take this ensity to be close to that of total equilibrium an write it in the form ρ( l, X, 0) = ρ eq [ + R( ] l, X, 0), (2) where the equilibrium istribution ρ eq is, for Lorentz gases, uniform with respect to the cells l, an with respect to the phase variables X in agreement with the assumption that the Lebesgue measure X is preserve by the flow Φ t. For the case of tracer iffusion, ρ eq is the equilibrium micro-canonical istribution for the N particles in a cell on the lattice. The initial eviation from equilibrium in the cell locate at l is enote by R( l, X, 0), which we assume to be Lebesgue integrable when weighte with the equilibrium istribution. Using the lattice Fourier transform, Eq. (), we can express this eviation in the form R( l, X, 0) = k e i k l R( k, X). (3) The phase space ensity (2) leas us to the efinition of the non-equilibrium measure of a set A l belonging to the phase space cell T l M corresponing to the lattice vector l : Definition 2 The equilibrium measure ν(a l ) of a set A l is efine by ν(a l ) Xρ eq, (4) A l an the non-equilibrium measure, µ t (A l ), of the same set by µ t (A l ) Xρ( l, X, t) = ν(a l ) + δµ t (A l ), (5) A l where δµ t (A l ) = X ρ eq R( l, X, t). A l (6) We will simplify matters a bit by assuming that the initial eviation from equilibrium R( l, X, 0) epens only upon the cell l but not on the initial phase X of the system within the cell. In this case the Fourier transform R( k, X) oes not epen upon X, either, an will henceforth be enote by F k.
4 The time epenent istribution function ρ( l, X, t) is the solution of Liouville s equation an is given by ρ( l, X, t) = ρ eq [ + R( ] l, X, t), (7) { = ρ eq + } k F k e i k L[Φ t ( l,x,0)]. Here L enotes the proection on the lattice coorinate, thus L [ Φ t ( ] l, X, 0) is the lattice vector of the cell in which a moving particle woul be locate at time t if it were in cell l at time t = 0, with phase X. We can express the time epenent eviation from total equilibrium in the form R( l, X, t) = k F k e i k [ l+ (X,t)], (8) where the backwar isplacement of the lattice vector l over a time interval t is efine by (X, t) = L [ Φ t ( ] l, X, 0) l. (9) We will make heavy use of the fact that (X, t) epens upon the time interval ( t, 0) an upon the phase point X, at the initial time, but not upon the initial cell l. In other wors, the perioicity of the lattice an the ynamics prouce a wining number (X, t) that oes not epen upon the cell in which the traectory is locate at the initial time. We remark here that the long time limit of the isplacement vector (X, t) may be a willy varying function of the phase coorinate X. Thus we expect that the ecaying moes of the time epenent istribution (7) are singular functions of the phase coorinates. We efine the time epenent ensity of the tracer particle by integrating the phase space ensity over the coorinates X of an elementary phase space cell of the lattice [2] : n( l, t) X ρ( l, X, t) = µ t (M l ). (0) M Using Eq. (8), we obtain the ensity as n( l, t) = ν(m) () [ + k F k e i k l M Xρ ] eq e i k (X,t) M X ρ. eq We consier times t that are long compare to the mean time between collisions of the moving particles, but short compare to the time neee for the system to relax to total equilibrium. For such times, we expect the time epenent eviation from total equilibrium to ecay exponentially with a rate s k given by the van Hove relation [2] s k lim t t ln ei k (X,t) M, = lim t t ln M Xρ eq e i k (X,t) M Xρ eq, (2) which gives the ecay rate of a hyroynamic moe of iffusion of wave number k. An expansion in powers of the wave number gives s k = D k 2 + O( k 4 ), (3) with iffusion coefficient D. We have here assume that the iffusive motion of the tracer particle is invariant uner space inversion so that all the o powers of the wave number vanish. We notice that the existence of the successive terms of the expansion in powers of the wave number epens on the existence of the super-urnett an higher iffusion coefficients[3], which has been recently prove for the har-isk perioic Lorentz gas with a finite horizon [4]. The efinition (2) shows that M Xρ eq e i k (X,t) M Xρ = C( k, t) e skt, (4) eq where C( k, t) is a function of the wave number k with a sub-exponential epenence on time, i. e. lim t t lnc( k, t) = 0. (5) Accoringly, the tracer ensity can be written as [ n( l, t) = ν(m) + ] k F k e i k l C( k, t) e s kt (6) Notice that n( l, t) obeys a form of the iffusion equation, appropriate for our lattice system, given by n( l, t) t = ν(m) = D ν(m) k F k e i k l C( k, t)s k e skt + k F k k 2 e i k l C( k, t)e s kt +, (7) which, in the scaling limit where the size of the unit cell becomes small an for large times an small wave numbers, is the iffusion equation: IV. n( l, t) t D 2 n( l, t) l 2. (8) THE HYDRODYNAMIC MODES A. The hyroynamic measures Spatially perioic eviations from total equilibrium characterize by the wave number k relax exponentially at the rate given by van Hove s relation [2]. Our purpose is here to etermine the non-equilibrium state corresponing to this moe of exponential relaxation. This
5 state can be efine as a measure, which we call a hyroynamic measure, associate with the hyroynamic moe of iffusion. It is the microscopic analog of the solutions exp(i k l Dk 2 t) of wave number k for the macroscopic iffusion equation (8). We introuce the hyroynamic measures by consiering the eviations from the equilibrium measure for a set A l : δµ t (A l ) = X ρ eq R( l, X, t), A l = Xρ eq A l k F k e i k [ l+ (X,t)], = ν(m) k F k e i k l M X ρ eq e i k (X,t) M Xρ eq A Xρ eq e i k (X,t) M Xρ eq e i k (X,t), (9) where we have use the property that the backwar isplacement (X, t) is inepenent of the initial lattice vector l so that the integral over the set A l is equal to the integral over the set A of the elementary cell at the origin of the lattice. In the last line, we have factorize the exponential ecay accoring to Eq. (4), which is inepenent of the set A, from a further factor, which epens on the set A but which is expecte to have a well-efine limit for t because both its numerator an enominator are expecte to ecay exponentially as exp(s k t). This observation motivates the Definition 3 The hyroynamic measure χ k (A, t) is efine by A χ k (A, t) ν(m) Xρ eq e i k (X,t) M Xρ eq e i k (X,t). (20) We emphasize that the hyroynamic measures are inepenent of the cell location l. We notice that the hyroynamic measures are complex measures because of the lattice Fourier transform. One important property of the hyroynamic measures is that the total hyroynamic measure of a unit cell is constant in time, as follows from the efinition 3. That is, Lemma 3. For the set A = M, i. e., the phase space region associate with an entire unit cell, the hyroynamic measure is Another observation is that χ k (M, t) = ν(m). (2) Lemma 3.2 If we make a k-expansion of the hyroynamic measure of the form χ k (A, t) = ν(a) +i k T(A, t)+ k k : g (A, t) +, (22) it follows from Eq. (2), that T(M, t) = 0 an that g(m, t) = 0, etc. Properties like these have alreay been use in the various symmetric, multi-baker moels [3, 4, 5, 6]. For a system with vanishing mean rift (X, t) M = 0, the two first coefficients of the k-expansion of the hyroynamic measures can be expresse as T(A, t) = ν(a) (X, t) A (23) g(a, t) = ν(a) (24) 2 [ (X, t) (X, t) M (X, t) (X, ] t) A with the efinition A A Xρ eq( ) A Xρ eq (25) Thanks to the hyroynamic measures Eq. (20) an Eq. (4), we finally erive from Eq. (9) an expression for the measure of a set A in cell M l : µ t (A l ) = ν(a) + δµ t (A l ), = ν(a) (26) + k F k e i k l C( k, t) e skt χ k (A, t).. Conservation of measure an e Rham-type equation Since the time evolution is a measure preserving Liouville operator, the measure of any set A remains constant as the set follows the motion of the system in phase space. Therefore we may express this conservation of measure as µ t+τ (A) = µ t (Φ τ A). (27) Here τ enotes some time interval, an Φ τ A is the preimage of the set A uner the flow, obtaine by following the backwar evolution of the points of A over a time interval τ. This simple result has some important consequences, among them, a e Rham-type equation for the hyroynamic measures. We suppose that the sets A are sufficiently small that all the points in them will flow through the same sequence of cells over some large time interval ( T τ T ). In such a case the set of points Φ τ A are all in the same cell with location enote by l+ (X A, τ) which is etermine by the backwar evolution of an arbitrary phase point X A in the set A. Using Eq. (26), we can express the application of Eq. (27) to a set A l as
6 k F k e i k l C( k, t + τ) e (t+τ)s k χ k (A, t + τ) = k F k e i k [ l+ (X A,τ)] C( k, t) e skt χ k (Φ τ A, t). (28) Since this equation must be true for all cells l for all F k, an for all sets A satisfying the above conition, the only way it can be satisfie is if the integrans are equal almost everywhere. Equating the integrans leas to the equation C( k, t + τ) e s kτ χ k (A, t + τ) = e i k (X A,τ) C( k, t) χ k (Φ τ A, t). (29) Since the hyroynamic measures o not epen on the lattice vector, the pre-image Φ τ A uner the full flow over the lattice can be reuce to the pre-image φ τ A uner the flow efine with perioic bounary conitions insie the elementary cell M at l = 0. Now we suppose that the ynamics is hyperbolic in orer to assert that the hyroynamic measures χ k (A, t) reach asymptotic forms exponentially rapily, that is, on a time scale of the orer of the inverse of the positive Lyapunov exponent for the system. In this case, we can replace the hyroynamic measures in Eq. (29) by their asymptotic forms, enote by χ k (A), in the long time limit t. Moreover, we have the property that lim t C( k, t + τ) C( = (30) k; t) as a consequence of the sub-exponential behavior (5) of the functions C( k, t). In the long time limit t, combining Eqs. (29) an (30), we obtain the Lemma 3.3 The hyroynamic measures satisfy a e Rham-type equation : e s kτ χ k (A) = e i k (X A,τ) χ k (φ τ A). (3) Explicit solutions of this equation have been foun for multi-baker maps [7, 5, 6]. For the har-isk Lorentz gas, these solutions lea to the cumulative functions constructe in reference [7] with one-imensional sets A an in reference [] with two-imensional sets A. An alternative form of Eq. (3) is e s kτ χ k (φ τ A) = e i k (X φ τ A,τ) χ k (A). (32) Equation (32) has an expansion in powers of the wavenumber k which will be useful in the calculation of the rate of entropy prouction. In obtaining these expansions we will make use of the k-expansions of the χ k functions given in Eq. (22). The wave number expansion of Eq. (32) leas to the following equation for terms of orer k: k T(φ τ A) = k T(A) + ν(a) k (X φτ A, τ). (33) C. Partition of phase space an sum rules We consier a partition {A } of the elementary cell M at l = 0 of the phase space into isoint sets A : A M A = M, A i A =, i,, i. (34) We notice that the images φ τ A also form a partition of the elementary cell of phase space: M = A M φ τ A. (35) We can apply the e Rham-type equation (32) to one set A of the partition (34) an sum both members of Eq. (32) over all the sets A M to obtain: e s kτ with the notation χ k (φ τ A ) = e i k χ k (A ), (36) (X φτ A, τ). (37) Since the sets φ τ A form a partition of M into isoint sets, we infer from Eq. (2) that χ k (φ τ A ) = χ k (M) = ν(m). (38) so that Eq. (36) becomes ν(m) e s kτ = e i k χ k (A ). (39) Now, we perform a wavenumber expansion of both members of Eq. (39) using Eqs. (3) an (22). Using the properties ν(a ) = ν(m), (40) T(A ) = T(M) = 0, (4) g(a ) = g(m) = 0, (42) the ientification of the terms which are of the first an secon powers of the wave number k gives us the following two sum rules: ν(a ) = 0, (43) [ T(A ) + T(A ) + ν(a )] = 2 D τ ν(m). (44)
7 Equation (44) is funamental for the following evelopment because it constitutes a sum rule relating the iffusion coefficient to the first coefficients T(A ) of the wavenumber expansion of the hyroynamic measures which is linear in the wave number k. The measures T(A ) have been interprete elsewhere as the stationary non-equilibrium measures associate with a graient of concentration of tracer particles across the system. In the case of the multi-baker maps, T(A ) is given by the ifference of the Takagi function at both ens of the oneimensional sets A [6]. The sum rule (44) thus relates the iffusion coefficient to the generalization of the Takagi function for the present system. V. ENTROPY PRODUCTION A. Definitions In this section we are going to calculate of the rate of irreversible entropy prouction over a time τ, assuming that t τ. For our calculation of the rate of entropy prouction in a unit cell of the perioic lattice, we use a partition of the total phase space into the small isoint sets A efine above Eq. (28). We suppose that the partition is invariant uner the spatial translations T l. The phase space cell locate at the lattice vector l is ecompose by this partition as: M l = A M l A. (45) This partition can be seen as a translationally invariant gri extening over the whole phase space. We begin by efining the entropy of the lattice cell M l at time t as the coarse-graine entropy of this cell with respect to the partition (45): S t (M l {A }) µ t (A ) ln µ t(a ) ν(a ) A M l +S eq (M l {A }), (46) where we have set oltzmann s constant equal to unity, k =. The first term on the right han sie of Eq. (46) is the non-equilibrium relative entropy with respect to the equilibrium entropy for this partition. The equilibrium entropy is given by S eq (M l {A }) = ν(a ) ln ν(a ), (47) c A M l where c is a constant which fixes the absolute value of the equilibrium entropy. The time variation of the entropy over a time interval τ is of course only ue to the change in the relative entropy, an is efine as the ifference τ S(M l ) S t (M l {A }) S t τ (M l {A }), = S t (M l {A }) S t (Φ τ M l {Φ τ A }). (48) On the other han, the entropy flow is efine as the ifference between the entropy which enters the cell M l an the entropy which exits that cell: τ es(m l ) S t τ (Φ τ M l {A }) S t τ (M l {A }), = S t (M l {Φ τ A }) S t (Φ τ M l {Φ τ A }). (49) Accoringly, the entropy prouction over a time τ, assuming that t τ, is efine as τ i S(M l ) τ S(M l ) τ es(m l ), = S t (M l {A }) S t (M l {Φ τ A }). (50). Calculation of the entropy prouction Equation (50) gives the expression of the entropy prouction as the ifference between the entropy with respect to the original partition into sets A an the entropy with respect to a partition into sets which are the images Φ τ A of the sets A after time τ. We notice that each set Φ τ A belongs to a single unit cell M l by a previous assumption. Moreover, since the partition is invariant uner translation from cell to cell, the partition {Φ τ A } is ientical to the partition {φ τ A } obtaine by using the flow on the torus. Written out in full, this entropy prouction is τ i S(M l ) = µ t (A ) ln µ t(a ) ν(a ) A M l + µ t (φ τ A ) ln µ t(φτ A ) (5) ν(a ) φ τ A M l where we have use ν(φ τ A ) = ν(a ). Next we expan in powers of the eviations of the measures from their
8 equilibrium values an fin τ i S(M l ) = [δµ t (φτ A )] 2 ν(a ) φ τ A M l [δµ t (A )] 2 ν(a ) A M l 2 2 + O ( δµ 3 t).(52) We now use the explicit forms for the measures δµ t (A) given by the secon term on the right han sie of Eq. (26). After some algebra an the use of the conservation of measures, as well as the summation formulas (4) an (42), we fin that the right han sie of Eq. (52) becomes τ i S(M l ) = 2 2 e (s k +s k2 )t k F k ν(a ) k k2 : k 2 F k2 e i l ( k + k 2) C( k, t) C( k 2, t) [ T(A ) T(A ) T(φ τ A ) T(φ τ A )]. (53) Here the summation is over the sets A that form a partition of the unit cell M l. Now we use the ientity (33) an the sum rule (44) to obtain our central result τ i S(M l ) = D τ ν(m) 2 k F k k 2 F k2 k k 2 e i l ( k + k 2) C( k, t) C( k 2, t) e (s k +s k2 )t [ n( D τ ] 2 l, t) n eq, (54) l since the tracer ensity is expresse accoring Eq. (6) an n eq = ν(m). Here we have use the isotropy of the motion, in orer to eliminate correlations between the isplacements in orthogonal irections. This results in the factor k k 2 appearing in the integran in Eq. (54). We have also implie a scaling limit in orer to write the last line of this equation. VI. CONCLUSIONS In this paper, the irreversible entropy prouction has been erive from statistical mechanics for a process of iffusion in perioic, eterministic ynamical systems. The erivation starts from Gibbs coarse-graine entropy, assumes that the system is spatially perioic, an that the ynamics satisfies Liouville s theorem, is chaotic an mixing. We chose a coarse graining partition of phase space which has the property that any two traectories starting in the same set of the partition will remain close together over some specifie time interval. The central quantities appearing in our erivation are the hyroynamic measures. They efine, at the microscopic level, the hyroynamic moes of iffusion which are exponentially ampe at a rate given by the van Hove ispersion relation for iffusion. These hyroynamic measures escribe the approach to the thermoynamic equilibrium uner the iffusion process. In eterministic systems, the hyroynamic measures turn out to be singular. Inee, the quantity T(A) is a measure escribing a nonequilibrium stationary state corresponing to a graient of concentration across the system an, in the multibaker, the cumulative function of this measure is known to be the continuous but non-ifferentiable Takagi function [5]. On the other han, it has been shown elsewhere that the hyroynamic measures are singular in perioic Lorentz gases [7]. It shoul be emphasize here that the hyperbolicity of the system is use to argue that the time epenent hyroynamic measures, efine by Eq. (20) approach their asymptotic forms, with variations on arbitrarily fines scales in phase space, on time scales etermine by the positive Lyapunov exponents, which are very short compare to the time scales necessary for the relaxation of the ensity istribution to equilibrium. Of course, not all systems with macroscopic iffusion an a positive rate of entropy prouction are chaotic, e.g wintree moels with non-overlapping trees[8]. In such systems the mechanism responsible for positive entropy prouction is expecte to evelop on algebraic rather than exponential time scales, as is the case for chaotic systems. Moreover, if the hyroynamic measures were regular, expression (54) for the entropy prouction woul vanish, as shown in the case of the multi-baker map [3, 6]. Accoringly, the singular character of the hyroynamic measures plays an essential role in the positiveness of the entropy prouction expecte from irreversible thermoynamics. A further point we want to emphasize is that the present erivation leas naturally to a positive entropy prouction in agreement with the secon law of thermoynamics. Finally we point out that our erivation of the expression given by irreversible thermoynamics
9 for the rate of entropy prouction applies not only to Lorentz gases, but also to tracer iffusion taking place in a spatially perioic N particle system, where the mechanism for tracer iffusion is provie by the particles interacting with each other. Without this interaction, the motion of the tracer particle woul be ballistic, an the mean square isplacement woul grow quaratically with time. The next steps to be taken in the evelopment of this approach to the theory of entropy prouction in the relaxation of flui systems to thermal equilibrium is to generalize the metho given here for tracer iffusion in a perioic N particle system to viscous an heat flows. We shall report on the erivation of entropy prouction for these other transport processes in future publications. Finally one woul like to remove the restriction to perioic systems, an to consier the entropy prouction for a general, isolate N particle system relaxing to thermal equilibrium. This remains open for future work. It is necessary mention that our approach to entropy prouction in fluis has been criticize by Rononi an Cohen in a series of papers[9]. This is not the place to provie a etaile response to their criticisms, which we will o in separate publications. However, it is appropriate here to mention two issues that are of some importance for our response to their comments. () Rononi an Cohen fin it troublesome that the time neee for our metho to be applicable epens on the nature of the partition chosen for the calculation of the relative entropy an its change with time. While this observation is inee correct, the consequences are not problematic. The essential point to note is that for a non-zero rate of irreversible entropy prouction to result from our metho, we require a partition which is coarser than the scale of variation of the non-equilibrium istribution function. As we have trie to make clear, the non-equilibrium istribution function evelops a fractal structure on a time scale set by the magnitue of the largest positive Lyapunov exponent. For the systems escribe here, this is a time scale which is very short compare to the time scale over which the system relaxes to equilibrium. In fact, this time scale is on the orer of the mean free time between collisions, a scale sufficiently short that more traitional methos for computing the rate of entropy prouction, base upon the oltzmann equation, for example, o not apply, either. (2) Rononi an Cohen have also obecte to our use of multi-baker maps an Lorentz gases as examples of systems to which non-equilibrium thermoynamics might be applie. While one might obect to their view of the utility of simple moel systems for the evelopment of physical intuition, our calculation of the rate of entropy prouction for tracer iffusion in an N particle system, shows that the metho can be applie to much more general systems than multi-baker maps an Lorentz gases. Our future work will be evote to seeing exactly how general this metho might be. Acknowlegments: The authors woul like to thank the Max Planck Institute for Physics of Complex Systems, Dresen, for its hospitality while part of this work was carrie out. JRD thanks the National Science Founation for support uner grant PHY 98-20824; PG thanks the National Fun for Scientific Research (FNRS elgium) for financial support; TG thanks the Israeli Council for Higher Eucation an the Feinberg Postoctoral Fellowships Program at the Weizmann Institute of Science for financial support. [] J. W. Gibbs, Elementary Principles in Statistical Mechanics Yale U. Press, New Haven, (902); reprinte by Dover Publ. Co, New York, (960). [2] A non-technical account may be foun in J. R. Dorfman, An Introuction to Chaos in Nonequilibrium Statistical Mechanics, (Cambrige U. Press, Cambrige, 999); A more mathematical iscussion is given by A. Lasota an M. C. Mackey, Chaos, fractals, an noise, (Springer- Verlag, New York, 994). [3] P. Gaspar, J. Stat. Phys. 88 25, (997). [4] S. Tasaki an P. Gaspar Theoretical Chemistry Accounts 02, 385, (999). [5] S. Tasaki an P. Gaspar, J. Stat. Phys. 0, 25, (2000). [6] T. Gilbert, J. R. Dorfman, an P. Gaspar, Phys. Rev. Lett. 85, 606, (2000). [7] T. Gilbert, an J. R. Dorfman, Physica A 282, 427, (2000). [8] S. e Groot an P. Mazur Non-equilibrium Thermoynamics, (North-Hollan, Amsteram, 962); reprinte by Dover Publ. Co., New York, (984). [9] L. A. unimovich, an Ya. G. Sinai, Commun. Math. Phys. 78, 247, 479, (980). [0] A. Knauf, Commun. Math. Phys. 0, 89, (987); Ann. Phys. (N. Y.) 9, 205, (989). [] P. Gaspar, Chaos, Scattering an Statistical Mechanics, (Cambrige U. Press, Cambrige, 998). [2] L. Van Hove, Phys. Rev. 95, 249, (954). [3] See, for example, H. van eieren, Rev. Mo. Phys. 54, 95, (982). [4] N. I. Chernov an C. P. Dettmann, Physica A279 37, (2000); C. P. Dettmann preprint nlin.cd/0003038. [5] S. Tasaki, an P. Gaspar, J. Stat. Phys. 8, 935, (995). [6] T. Gilbert, J. R. Dorfman, an P. Gaspar, Nonlinearity 4, 339, (200). [7] P. Gaspar, I. Claus, T. Gilbert, an J. R. Dorfman, Phys. Rev. Lett. 86, 506, (200). [8] E. H. Hauge an E. G. D. Cohen, J. Math. Phys. 0, 397, (969). [9] L. Rononi an E. G. D. Cohen, Nonlinearity, 3, 905, (2000); an (to be publishe). [20] With the avent of quantum mechanics this size turne out to be equal to h f where h is the Planck constant an f the number of egrees of freeom. In classical systems, the size of the cells remains arbitrary, which is relate to the fact that classical physics oes not fix the constant of entropy.
0 [2] We will use the term ensity of the tracer particle to refer both to the ensity of moving particles in a Lorentz gas as well as to the ensity of the tracer particle in the system of N moving, interacting particles.