arxiv: v1 [cond-mat.stat-mech] 16 Nov 2018

Transcription

1 Journal of Statistical Physics manuscript No. will be inserted by the editor Statistical mechanical expressions of slip length Hiroyoshi Nakano Shin-ichi Sasa arxiv:8.6775v [cond-mat.stat-mech] 6 Nov 8 Received: date / Accepted: date Abstract We establish a microscopic foundation for the partial slip boundary condition provide two statistical mechanical expressions of the slip length based on microscopic dynamics linearized fluctuating hydrodynamics. We develop two linear response theories for microscopic systems assuming the separation of scales between the microscopic behavior, such as collision of particles, the macroscopic behavior, such as relaxation to global equilibrium. The first linear response theory describes the relaxation process of the fluid, the second describes the non-equilibrium steady state. From these linear response theories, we derive the partial slip boundary condition obtain two expressions for the slip length, involving the Green Kubo formulae on different time scales. These derivations suggest that the partial slip boundary condition is robustly realized when the slip length remains microscopic. We validate these expressions from the linearized fluctuating hydrodynamics the partial slip boundary condition. By calculating the microscopic expressions of the slip length exactly, we find that both microscopic expressions hold in the fluctuating hydrodynamic description giving the same value despite their difference in form. Keywords Hydrodynamics Boundary Condition Slip Length Green Kubo formula Linearized fluctuating hydrodynamics Introduction Over the past two decades, interest in applying hydrodynamics to nano- micro- scale system have increased, because of the remarkable developments in experimental techniques for small-scale systems [,,3]. From a large number of studies on this topic, it is confirmed that certain specific fluids confined in such systems slip on the solid surface but behave accurately as Newtonian liquids in the bulk region [4,5,6,7,8,9,,,,3,4,5,6,7,8,9,,,]. This result indicated the breakdown of macroscopic hydrodynamics that the behavior of fluids are described by the solutions of the Navier Stokes equation stick boundary condition. Here, the stick boundary condition means that a fluid has no relative velocity when in contact a solid surface [3]. Determining the correct boundary conditions to replace the stick boundary condition, is one of the important challenges in hydrodynamics. Slipping at the solid surface is often characterized by the partial slip boundary condition [4]. Specifically, the slip velocity v s of the fluid at the wall is taken to be linearly proportional to the shear rate γ at the wall as [5,6,7] v s = b γ, H. Nakano S. Sasa Department of Physics, Kyoto University, Kyoto 66-85, Japan h.takagi@scphys.kyoto-u.ac.jp, sasa@scphys.kyoto-u.ac.jp

2 Hiroyoshi Nakano, Shin-ichi Sasa where the slip length b represents the distance beyond the surface of the wall at which the fluid velocity extrapolates to zero. Much effort has been spent in investigating factors that affect the slip length such as surface roughness [8] wettability [9, 3, 3]. We shall try to describe the partial slip boundary condition based on the underlying microscopic theories. This microscopic foundation provides the limitations of the partial slip boundary condition the first step in considering the slip phenomena beyond the partial slip boundary condition. Despite its importance, this boundary condition is not fully understood compared the hydrodynamic equations whose microscopic foundation has been studied over the past century [3, 33, 34, 35, 36, 37, 38,39]. Pioneering work on the microscopic description of the partial slip boundary condition began Bocquet Barrat [4]. For a fluid in a semi-infinite container a stationary plane wall, they derived the Green Kubo formula for the friction force between fluid wall, γ BB, in the form [4,4, 4,43,44,45,46,47,48,49] γ BB = Sk B T dt F BB t F BB eq, where S denotes the surface area, k B the Boltzmann constant, T the temperature of the fluid, eq the canonical ensemble average at temperature T, F BB t the total force between wall fluid at time t. The Green Kubo formula does explain the dependence of slip length on microscopic properties such as the density near the wall the interaction between molecules of the fluid the wall [4, 4,43,44]. There remains an unsatisfactory point in their derivation. We recall that the derivation of hydrodynamics is focused towards the separation of scales a small parameter representing the ratio of the characteristic lengths of each scale. However, in these previous studies, the role of separation of scales is unclear. Bocquet Barrat derived as the linear response theory in the gap between the velocities of the fluid the wall. Because the hydrodynamic equations are derived as the linear response theory in a scale separation parameter [38], the approach of Bocquet Barrat does not correspond to that for the derivation of the hydrodynamic equations. Also, they emphasized that the order between the thermodynamic limit the infinite time limit is important for obtaining the finite slip length from the expression [4, 45, 46, 47, 48]. We shall show that, instead of introducing thermodynamic limit the infinite time limit, the correct value for the slip length can be obtained from the expression by introducing a cut-off length time specifying the separation of scale. Furthermore, our approach reveals that the boundary condition is a surface property, not a bulk property. We recall another approach to deriving the boundary condition that is related to Stokes law. Suppose a large sphere is moving in a fluid. From microscopic mechanics, Kirkwood derived the Fokker Planck equation that described the macroscopic motion of the sphere obtained a microscopic expression for the friction coefficient of the macroscopic sphere in the fluid, γ K, in the form [5] γ K = 3k B T dt F K t F K eq, 3 where k B denotes the Boltzmann constant, T the temperature of the fluid, eq a canonical ensemble average at temperature T, F K t the total force exerted on the sphere by the fluid at time t. The essence of Kirkwood s derivation is the separation of time scales between the macroscopic sphere the particles of the medium [36, 5, 5]. Therefore, Kirkwood s formula generally holds independently of the nature of the fluid e.g., dilute gas [53] fluctuating hydrodynamics [54]. In Kirkwood s approach, the macroscopic motion of the sphere is considered in contrast to the approach of Bocquet Barrat, which focuses solely on the macroscopic motion of fluid not wall. However, Kirkwood s formula also involves the boundary condition when the hydrodynamic equations hold in the bulk. From this point of view, several researchers have attempted to derive the boundary condition by directly calculating the friction coefficient of a macroscopic sphere γ K [55,56,57]. However, their calculations remained phenomenological. In addition, there are no discussions regarding the relationship between the boundary conditions derived in this approach the partial slip boundary condition.

3 Statistical mechanical expressions of slip length 3 The aim of this paper is to derive the partial slip boundary condition from the microscopic dynamics obtain the statistical mechanical expressions of the slip length. For this purpose, we first demonstrate that the two foregoing approaches are understood as two linear response theories to describe the macroscopic behavior on different time scales. Specifically, we consider the fluid enclosed between two parallel walls. By moving the walls parallel to each other, uniform shear flow is induced. In this system, we consider two characteristic time scales; one is a time scale of the microscopic motion of molecules, denoted by τ micro, the other is a time scale of the global equilibrium of fluid, denoted by τ macro. We assume the separation of scales, which is represented by τ micro τ macro. Given this separation, we introduce lower upper cut-off times concerning the macroscopic behavior of fluid as τ meso τ ss, which satisfy τ micro τ meso τ macro τ ss. We develop two linear response theories that describe the macroscopic behavior on the time scales τ meso τ ss. We obtain two statistical mechanical expressions of the slip length from these linear response theories. The first expression is obtained from the linear response theory on the τ meso scale, which coincides that of Bocquet Barrat, whereas the second one is obtained from the linear response theory on the τ ss scale, which is related to the application of Kirkwood s formula to our system. Our derivation of the partial slip boundary condition from the microscopic dynamics suggests that the necessary condition to realize the partial slip boundary condition is that the slip length remains microscopic. We validate these two expressions by analyzing the fluctuating hydrodynamics. Specifically, we consider the linearized fluctuating hydrodynamics the partial slip boundary condition. We explicitly calculate the time correlation functions contained in the two linear response theories e.g., the correlation function between the forces acting on the walls. From this calculation, we demonstrate that the two linear response theories derived from the underlying microscopic dynamics are consistent the linearized fluctuating hydrodynamics the partial slip boundary condition when the slip length remains microscopic. Then, we confirm that the two microscopic expressions of the slip length give the same value, despite these apparent difference in form. The remainder of this paper is organized as follows. In Sect., the microscopic setup of our model is introduced. In Sect. 3, we derive the exact formula to describe the time evolution of any ensembleaveraged quantities. This formula is the starting point for developing linear response theories. In Sect. 4, we introduce the separation of time scales explain the behavior of correlation functions of the momentum flux. In Sects. 5 6, we develop the linear response theories describing the macroscopic behavior on the time scales τ meso τ ss. Then, on the basis of these linear response theories, we derive the partial slip boundary condition obtain two statistical mechanical expressions of the slip length. In Sect. 7, we validate these expressions of the slip length from the linearized fluctuating hydrodynamics the partial slip boundary condition. The final section is devoted to a brief summary concluding remarks. Throughout this paper, the superscript a, b represents the indices in Cartesian coordinates x, y, z. Model We introduce a model for studying the boundary condition in fluid dynamics. A schematic illustration is shown in Fig.. The fluid consists of N particles that are confined to an L x L y L cuboid. The position momentum of the ith particle are denoted by r i, p i i =,,, N. We impose periodic boundary conditions along the x y axes introduce two parallel walls so as to confine particles in the z direction. We represent the two walls as collections of material points placed in planes z = z = L, which are referred to as walls A B, respectively. We move the walls in the x direction. Let q A i t, i =,,, N A, q B i t, i =,,, N B, be the position of the material points that consist of walls A B, respectively. The motion of the walls is given by the uniform motion of the material points: for i =,,, N A, q A i t = q A i + V A t 4 q B i t = q B i + V B t 5

4 4 Hiroyoshi Nakano, Shin-ichi Sasa Fig. Schematic illustration of our model. for i =,,, N B, where V A = V A,, V B = V B,,. The collection of positions momenta of all particles is denoted by Γ = r, p,, r N, p N. Similarly, the collections of positions of the material points that comprise walls A B are denoted by Γ A = q A,, q A N A Γ B = q B,, q B N B. The Hamiltonian of the system is given by HΓ ; Γ A, Γ B = N p i m + U{r i} N ; Γ A, Γ B 6 U{r i } N ; Γ A, Γ B N N U FF r i r j + U AF r i ; Γ A + U BF r i ; Γ B. 7 i<j U FF r describes the interaction potential between particles. U AF r; Γ A U BF r; Γ B denote the potentials of walls A B, respectively. U AF r; Γ A may be expressed as N A U AF r; Γ A U A r qj A, 8 j= where we introduce the interaction between the material point the fluid particle, U A r. U BF r; V B has a similar form. We employ the stochastic thermostat to thermalize the fluid. F bath r; T denotes the stochastic force applied to the particle at a point r by the thermostat temperature T. Then, the particles obey m d r a i dt = j i V FF r i r j r a i U AFr i ; Γ A r a i U BFr i ; Γ B r a i + F a bathr i ; T. 9 For this study, the configuration of the stochastic thermostat is: i the dynamics 9 satisfies the local detailed balance condition Sec. 3. gives a detailed account of this condition; ii the dynamics 9 are Galilean invariant in the x direction, implying that 9 is unchanged under the Galilean transformation rt r t = r vt v = v,, ; iii the characteristic strengths of interactions between particles, between a particle wall, between a particle the thermostat, denoted by F int, F wall F bath, respectively, obey F bath F int, F wall.

5 Statistical mechanical expressions of slip length 5. Microscopic expression of observed quantities Let Γ t denote the collection of positions momenta of all particles composing the fluid at time t. Similarly, let Γ A t Γ B t denote the collection of positions of all material points composing walls A B, respectively. We define the microscopic mass density field ˆρr; Γ the microscopic momentum density field ˆπ a r; Γ as ˆρr; Γ N mδ r r i, N ˆπ a r; Γ p a i δ r r i for a given microscopic configuration Γ. These microscopic densities satisfy the microscopic balance equations: ˆρr; Γ t t + ˆπa r; Γ t r a =, 3 ˆπ a r; Γ t t + Ĵ ab r; Γ t r b = + N ˆF a A r i t; Γt A + ˆF Br a i t; Γt B δ r r i t N Fbathr a i tδ r r i. 4 Here, ˆF a A r; Γ A ˆF a B r; Γ B are given by N A ˆF Ar; a Γ A j= N B ˆF Br; a Γ B j= U A r q A j r a, 5 U B r q B j r a, 6 which represent the total force acting on the fluid at a point r from walls A B, respectively. The microscopic momentum current Ĵ ab r; Γ is given by Ĵ ab r; Γ i p a i pb i m δr r i + i<j ˆF a ijγ r b i r b jdr; r i, r j 7 Dr; r i, r j dξδr r i r j r i ξ, 8 ˆF ijγ a U FF r i r j ri a. 9 Here, we recall the condition, which implies that the interaction between the particle the thermostat is sufficiently weak. From this condition, the last term on the right-h side in 4 may

6 6 Hiroyoshi Nakano, Shin-ichi Sasa be neglected in the discussion of this paper. Therefore, we employ as the microscopic balance equation for ˆπ a r; Γ ˆπ a r; Γ t t + Ĵ ab r; Γ t r b = N ˆF a A r i t; Γ A t + ˆF a Br i t; Γ B t δ r r i t, instead of 4. [Γ ] = Γ t denotes the history of the system over the time interval [, τ]. For a given initial <t<τ state Γ tot Γ, Γ A, Γ B, T [Γ ] Γ tot, V A, V B represents the probability density that [Γ ] is realized wall velocities V A V B. We assume that the initial state Γ is chosen according to a given distribution P ini Γ. We now focus on any two microscopic quantities ˆfΓ ĝγ given as functions of microscopic state Γ. The ensemble average of ˆfΓ τ at time t = τ, fτ V A,V B, is expressed in terms of P ini Γ T [Γ ] Γ tot, V A, V B as fτ V A,V B = d[γ ] ˆfΓ τ P ini Γ T [Γ ] Γ tot, V A, V B. Similarly, the ensemble-averaged correlation function between ˆfΓ τ at time t = τ ĝγ at time t = is expressed as fτg V A,V B = d[γ ] ˆfΓ τ ĝγ P ini Γ T [Γ ] Γ tot, V A, V B. In this paper, the motion of the fluid is characterized by macroscopic mass density field ρr, t; V A, V B, macroscopic velocity field v a r, t; V A, V B, stress tensor field σ ab r, t; V A, V B. These quantities are defined as ρr, t; V A, V B ˆρr, t V A,V B, 3 π a r, t; V A, V B ˆπ a r, t V A,V B, 4 v a r, t; V A, V B πa r, t; V A, V B ρr, t; V A, V B, 5 σ ab r, t; V A, V B J ab r, t; V A, V B ρr, t; V A, V B v a r, t; V A, V B v b r, t; V A, V B 6 J ab r, t; V A, V B = Ĵ ab r, t V A,V B. 7 Furthermore, we focus on the force acting on the fluid from each wall. The microscopic definition is given by ˆF a αγ t ; Γ α t = N ˆF αr a i t; Γt α, 8 where α = A, B. The ensemble average is denoted by F a αt; V A, V B = ˆF a αγ t ; Γ α t V A,V B.

7 Statistical mechanical expressions of slip length 7 3 Generalized Green Kubo formula In this section, we give the exact formula describing the time evolution of any ensemble-averaged quantities. We prepare initial states as follows. First, we assume that the system is in equilibrium at temperature T. Next, we apply an impulse force murδt to this system where we set u = u x z,,. Then, the probability density just after t =, P leq Γ Γ A, Γ B, is represented by P leq Γ Γ A, Γ B = [ Z k B T, ur exp HΓ ; Γ A, Γ B k B T N ] u x z i p x i. 9 We consider that the initial microscopic state Γ is chosen according to 9. For any physical quantity ˆfΓ any velocity field u x z, we have fτ V A,V B = f V A,V B + τ k B T + τ ds k B T d 3 r k B T ds N VA u x z i fsfar x i, V A,V B N VB u x z i fsfbr x i, V A,V B τ ds ux z fsj xz r, V A,V B, 3 z which is the exact formula to describe the time evolution of fτ V A,V B. We refer to 3 as the generalized Green Kubo formula. This formula is the starting point for developing the linear response theories in the subsequent sections. In this section, we show a proof of the generalized Green Kubo formula 3. The key point in our derivation is an efficient use of a non-equilibrium identity, which is known as the local detailed balance condition. 3. Local detailed balance condition By considering the energetics for a given history [Γ ], we define the heat transferred into the fluid, Q [Γ ] Γ tot, V A, V B, the work applied to the fluid from the walls, W [Γ ] Γ tot, V A, V B as Q [Γ ] Γ tot N, V A, V B = W [Γ ] Γ tot, V A, V B = ˆΓ ˆΓ ds py i m The first law of the thermodynamics is expressed in the form ṗ y i + HΓ s; Γ A s, Γ B s r y i, 3 ds x V A ˆF A Γ s ; Γs A x + V B ˆF B Γ s ; Γs B. 3 HΓ t, t HΓ, = W [Γ ] Γ tot, V A, V B + Q [Γ ] Γ tot, V A, V B. 33 We use the convention that work applied to the fluid is positive heat dissipated from the fluid is negative. We write the time reversal of the microscopic state Γ as Γ = r, r,, r N, p, p,, p N, that of the history of the system [Γ ] as [Γ ] = Γτ t. The local detailed balance condition is expressed as T [Γ ] Γ tot, V A, V B T [Γ ] = exp Γ tot τ, V A, V B <t<τ Q [Γ ] Γ tot, V A, V B k B T, 34

8 8 Hiroyoshi Nakano, Shin-ichi Sasa where Γ tot = Γ, Γ A, Γ B. As mentioned in Sec., we assume that the dynamics 9 satisfies the condition 34 by choosing the appropriate stochastic thermostat. We note that the local detailed balance condition is derived from the microscopic dynamics under some assumptions e.g., Hamiltonian particle systems [58]. In this sense, the local detailed balance condition is often referred to as the local fluctuation theorem. 3. Derivation of 3 By substituting 9 34 into, we obtain fτ V A,V B = Wleq [Γ ] Γ tot d[γ ]fγ τ exp, V A, V B, u k B T P leq Γ τ Γ A τ, Γ B τ T [Γ ] Γ tot τ, V A, V B W leq [Γ ] Γ tot, V A, V B, u = W [Γ ] Γ tot N N, V A, V B u x z i tp x i t u x z i p x i. 36 We rewrite the right-h side of 35 by using the ensemble average along the time-reversal history. By noting we obtain 35 W [Γ ] Γτ tot, V A, V B = W [Γ ] Γ tot, V A, V B, 37 fτ V A,V B = d[γ ]f Γτ exp W leq [Γ ] Γτ tot, V A, V B, u k B T P leq Γτ Γτ A, Γτ B T [Γ ] Γτ tot, V A, V B = f exp W leq [Γ ] Γ tot, V A, V B, u VA, V B, 38 k B T where we have defined f Γ fγ. Similarly, we rewrite as fτg V A,V B = f g τ exp W leq [Γ ] Γ tot, V A, V B, u VA, V B. 39 k B T Next, using the microscopic force-balance equation, we rewrite W leq [Γ ] Γ, Γ A, Γ B, V A, V B, u. Multiplying by ur integrating over [, τ] all space, we obtain = N N u x z i tp x i t u x z i p x i τ ds N u x z i s ˆF a A r i s; Γs A + ˆF Br a i s; Γs B + Substituting 3 4 into 36, we obtain W leq [Γ ] Γ tot, V A, V B, u = τ τ ds τ ds d 3 r ux z Ĵ xz r; Γ s. 4 z N VA u x z i s x ˆF A r i s; Γs A + V B u x z i s x ˆF B r i s; Γs B ds d 3 r ux z Ĵ xz r; Γ s. 4 z

9 Statistical mechanical expressions of slip length 9 Using 38, 39 4, we write the derivative of fs V A,V B in s as fs V A,V B s = k B T N VA u x z i fsfar x i, V A,V B we obtain 3 by integrating 4 over [, τ]. + N VB u x z i fsf k B T Br x i, V A,V B d 3 r ux z fsj xz r, V A,V B. 4 z 4 Separation of length time scales 4. Assumption In our model, there are two length time scales. The first scale is characterized by the length time appearing in the molecular description. The maximum of the microscopic length scales, such as the molecular size, the interaction length, or the mean free path is denoted as ξ micro the corresponding time scale, which may be the collision time or the mean free time, is denoted as τ micro. The second scale is related to the global equilibration. The characteristic length time scales are given by the system size L the time τ macro necessary for the relaxation of the system to the global equilibrium, respectively. We assume the separation of length time scales, which is represented by τ micro τ macro, 43 ξ micro L Poperties of the correlation functions of the momentum fluxes We now focus on the correlation functions of the momentum fluxes in the equilibrium state. We prepare the initial states according to P leq Γ Γ A, Γ B u x z = consider the time evolution V A = V B =. eq denotes the ensemble average in this system, which corresponds to the average in equilibrium. The correlation functions of the momentum fluxes are, for example, FA xτf A x eq, J xz r, τfa x eq J xz r, τj xz r, eq. First, we focus on the τ-dependence of these correlation functions. We naturally conjecture that these correlation functions behave in the same manner respect to time τ. Therefore, we choose FA xτf A x eq as an example. By keeping the separation of time scales in mind, we assume that the τ-dependence of FA xτf A x eq is divided into four stages. First, when < τ τ micro, FA xτf A x eq decays to zero. After that, when τ micro τ τ macro, FA xτf A x eq remains zero. Eventually, when τ τ macro, FA xτf A x eq couples the slow macroscopic motion become non-zero. Finally, when τ τ macro, the relaxation is completed FA xτf A x eq again remains zero. In summary, we express this behavior as >, for < τ τ micro, FAτF x A x, for τ micro τ τ macro, eq 45, for τ τ macro,, for τ macro τ. The other correlation functions such as J xz r, τf x A eq J xz r, τj xz r, eq behave in the similar manner.

10 Hiroyoshi Nakano, Shin-ichi Sasa The macroscopic behavior of the fluid is related to the motion on the τ macro scale. By taking this into account, we introduce the lower upper cut-off times τ meso τ ss for the macroscopic behavior as τ micro τ meso τ macro τ ss. 46 Recalling the separation of length scales, we may choose the length scale corresponding to τ meso, i.e., ξ meso, to satisfy ξ micro ξ meso ξ macro. Because the correlation length of the momentum fluxes is reasonably estimated as the order of ξ micro when τ = τ meso, we assume for r r ξ meso. Similarly, J xz r, τ meso J xz r, eq 47 J xz r, τ meso F x A eq 48 is assumed for z ξ meso. Also, because there are no correlations when τ = τ ss, we have where r is any position. J xz r, τ ss J xz r, eq, 49 J xz r, τ ss F x A eq, Small parameters We introduce two small parameters representing the extent of the separation of scales, ɛ ξ micro /ξ macro, 5 ζ ξ meso /ξ macro. 5 In this study, we are mainly concerned a macroscopic system of system size m a shear rate sec. As the molecular size is typically m, ɛ is estimated to be of order 8. Recalling that ξ meso satisfies the relation ξ micro ξ meso ξ macro, we then assume that ξ meso is of order 6 m, which leads that ζ 4. As a result, the relationship between ɛ ζ is ζ = ɛ /. 53 In the reminder of this paper, we develop two linear response theories to describe the macroscopic motions of the fluid at time τ = τ meso τ = τ ss. 5 Linear response theory for the relaxation process of fluid We now develop the linear response theory that describes the macroscopic behavior at the time τ = τ meso. We prepare an initial state chosen according to 9, consider the time evolution during the time interval [, τ meso ]. In this interval, the macroscopic behavior of the fluid is determined by the velocity fields u x z given as the initial condition. We derive the constitutive equations of the fluid the velocity u x z. These constitutive equations lead to the partial slip boundary condition the expression of the slip length proposed by Bocquet Barrat [4]. In this section, we convert all the quantities to dimensionless forms setting ξ micro = m = V A V B =. We express dimensionless quantities a subscripted asterisk.

11 Statistical mechanical expressions of slip length 5. Scale separation parameter Let us assume that u x z depends on the small parameter ɛ, du x z dz = Oɛ, 54 d u x z dz = Oɛ, 55 u x V A = Oɛ, 56 V B u x L = Oɛ. 57 For example, u x z is given by u x z V A V B = K K k= πk sin L z + u x, 58 where K > is any integer u x is chosen to satisfy Because the dimensionless form of 58 is given by u x z = K K k= sin ɛπkz + u x, satisfies the assumptions 54, 55, We note that conditions confine us to a slip of order ɛ. These conditions are essential to the derivation of the partial slip boundary condition. Previously [59], we discussed what kind of walls satisfies when the fluid is in a steady state. In this section, we derive the linear response formula focusing on the first order approximation of 3 in ɛ. 5. Constitutive equation in bulk We focus on the behavior of the shear stress σ xz r, t ; V A, V B. Our starting point is 3 ˆfΓ = Ĵ xz r ; Γ τ = τ meso, J xz r, τ meso ; V A, V B = k B T τmeso ds + τmeso ds k B T d 3 r k B T N N τmeso VA u x z i J xz r, s FA r x i, V A,V B VB u x z i J xz r, s FB r x i, V A,V B ds du x z dz J xz r, s J xz r, V A,V B, 6 where we have used the property J xz r, ; V A, V B =. We now derive a simpler form of 6 by the separation of scales. First, we remark on the local approximation. As mentioned in Sect. 4, during time interval < τ < τ meso, the velocity field v x z, τ ; V A, V B remains u x z given by the initial condition,

12 Hiroyoshi Nakano, Shin-ichi Sasa the correlation functions of the momentum fluxes decay, e.g., F x A τ meso F x A eq. Therefore, the first term of 6 may be rewritten as k B T k B T τmeso ds N N V τ meso A u x z i eq VA u x z i J xz r, s FA r x i, V A,V B ds J xz r, s F x A r i, V A,V B. 6 With this approximation, we neglect non-local effects. Next, we consider the linear response approximation. We focus on 6. From 56, we note that 6 is at least of first order in ɛ. Because J xz r, τ meso FA x r i, V A,V B is estimated as the value in the steady state V A, V B u x z, we hypothesize that J xz r, τ meso FA x r i, V A,V B may be rewritten as J xz r, τ meso FA r x i, V A,V B J xz r, τ meso FA r x i, eq + O V A u x z i, ux z i, 6 z where the definition of eq is given in Sec. 4.. In the linear response theory, we neglect the second term. Collecting the local linear response approximations, we rewrite 6 as J xz r, τ meso ; V A, V B N V τ meso A u x k B T z i eq + k B T N V τ meso B u x z i eq k B T du x z dz d 3 r τmeso ds J xz r, s FA r x i, eq ds J xz r, s FB r x i, eq ds J xz r, s J xz r,. 63 eq Here, we define the bulk as the region ξ meso z L ξ meso. From 48, we find that the first second terms of 63 are zero in the bulk. Furthermore, 47 implies that the integral of the third term is independent of the walls when r is chosen in the bulk. Hence, we obtain J xz r, τ meso ; V A, V B = η du x z dz 64 η = k B T d 3 r τmeso ds J xz r, s J xz r, eq, 65 where η is independent of the walls. By substituting 64 into 6, we obtain where we have used v z z, τ meso ; V A, V B =. σ xz r, τ meso ; V A, V B = η du x z dz, 66

13 Statistical mechanical expressions of slip length Constitutive equation at the wall We next focus on the behavior of the force per unit area acting on the fluid from wall A, which is written as F x A t ; V A, V B /L x L y. By substituting ˆfΓ = ˆF x A Γ ; Γ A /L x L y τ = τ meso into 3, we obtain F x A τ meso ; V A, V B L x L y = + k B T L x L y k B T L x L y k B T L x L y τmeso τmeso d 3 r ds ds N N τmeso VA u x z i F x A s F x A r i, V A,V B VB u x z i F x A s F x B r i, V A,V B du x ds z F x dz A s J xz r, V A,V B. 67 We apply again the procedure given in Sect. 5.. We focus on the first term. Recalling ˆF x A r, Γ A in z > ξ micro, we note that only the particles lying in < z < ξ micro at time t = contribute to the first term. Therefore, the first term is rewritten as k B T L x L y V A u x τmeso ds k B T L x L y N VA u x z i F x A s F x A r i, V A,V B τmeso In the same manner, we rewrite 67 in the approximate form ds F x A s F x A V A,V B. 68 FA x τ meso ; V A, V B = γ AA τ meso V A u x L x L + γ AB τ meso V B u x L y γ AJ τ meso dux dz z = 69 γ αβ τ = k B T L x L y τ ds F x α s F x β eq, 7 γ αj τ = d 3 r k B T L x L y τ ds F x α s J xz r, eq, 7 where α, β = A, B. As the distance between walls A B is of order ξ macro, we have Then, 69 is rewritten as γ AB τ meso. 7 FA x τ meso ; V A, V B γ AA τ meso u x L x L V A γ AJ τ meso dux. 73 y dz z = Similarly, by focusing on the force acting on the fluid from wall B, we have F x B τ meso ; V A, V B L x L y γ BB τ meso V B u x L γ BJ τ meso dux dz z =L. 74

14 4 Hiroyoshi Nakano, Shin-ichi Sasa 5.4 Expression of the slip length To connect the constitutive equation in the bulk 66 at wall A 73, we assume the force-balance at time τ meso near the wall, dx dy σ xz r, τ meso ; V A, V B F A x τ meso ; V A, V B. 75 L x L y z =ξ meso L x L y Substituting into 75, we rewrite 75 in terms of u x z as du x η z γ AA τ meso u dz z =ξ x V A meso Furthermore, by exping the left-h side around z =, we obtain u x V A d A z A dux z dz z = d A = γ AJ τ meso dux z. 76 dz z = + d A ξ meso d u x z dz z = 77 η γ AA τ meso, 78 z A = γ AJ τ meso γ AA τ meso. 79 We now consider the order of each term of 77. From 53 55, the second term on the right-h side is estimated as Oɛ 3/. Because we are interested in the first order approximation in ɛ, we may neglect this term. Then, 77 is rewritten as u x V A = b meso A du x z dz z = 8 b meso A = d A z A, 8 which represents the partial slip boundary condition slip length b meso A. d A z A are independent of a specific value of τ meso because of the separation of scales, e.g., 45. This expression is the same as that proposed by Bocquet Barrat [4]. We note that the assumption 56 leads to b meso A Oɛ, which implies that d A z A remain microscopic, b meso A d A z A ξ micro. 8 Because the typical value of the slip length reported in recent experiments is about nm [,3], the result 8 is reasonable. We may rewrite 8 in the form u x z A V A = d A du x z dz z =z A. 83 In Ref. [4], taking the expression 83 into account, Bocquet Barrat called z A the hydrodynamic position of the wall defined d A as the slip length. However, we shall define the slip length as the length measured from the real position of the wall, z =. That is, we refer to b meso A as the slip length. We finally note that in the same manner we may derive the partial slip boundary condition for wall B V B u x L = d B z B dux z dz z =L 84

15 Statistical mechanical expressions of slip length 5 d B = η γ BB τ meso, 85 z B = γ BJ τ meso γ BB τ meso Linear response theory for the steady state We next develop the linear response theory that describes the macroscopic behaviors at time τ = τ ss. The fluid at time τ = τ ss is assumed to be in the non-equilibrium steady state. We note that the steady state is independent of the velocity u x z given as the initial condition. Instead, the steady state is characterized only by the wall velocities V A V B. Then, we focus on the linear response regime for V A V B. In this section, calculations are performed quantities retaining their dimensional units. 6. Linear response regime First, we consider the momentum density field π x r, τ ss ; V A, V B in the steady state. Substituting ˆfΓ = ˆπ x r; Γ τ = τ ss into 3 yields π x r, τ ss ; V A, V B = τss k B T ds N VA u x z i π x r, sfar x i, V A,V B + τss N ds k B T τss d 3 r k B T VB u x z i π x r, sf x Br i, V A,V B ds dux z dz π x r, sj xz r, V A,V B. 87 As the steady state is independent of the initial condition, the right-h side of 87 is independent of u x z. Then, by using the initial condition u x z =, we simplify 87 as π x r, τ ss ; V A, V B = V τss A k B T + V B k B T τss ds π x r, sf x A V A,V B ds π x r, sf x B V A,V B. 88 We apply a linear approximation to 88 in V A V B. We assume that ˆπ x r, τ ss V A,V B may be exped in V A V B, which leads to the expansion of the right-h side of 88 in V A V B. Within the linear response regime, 88 may be rewritten as where α = A, B. π x r, τ ss ; V A, V B = γ P A r, τ ss V A + γ P B r, τ ss V B 89 γ P α r, τ ss = τss ds π x r, sf k B T α x eq, 9

16 6 Hiroyoshi Nakano, Shin-ichi Sasa We can obtain linear approximations of other quantities by a similar procedure. In particular, we employ in the reminder of this section the following approximate expressions: F x A τ ss; V A, V B L x L y = γ AA τ ss V A + γ AB τ ss V B, 9 F x B τ ss; V A, V B L x L y = γ BA τ ss V A + γ BB τ ss V B, 9 where γ αβ τ α = A, B is defined by Force balance in steady state Because of the nature of the steady state, the microscopic balance equation for the momentum density,, leads to dxdyσ xz r, τ ss ; V A, V B = F A xτ ss; V A, V B = F B xτ ss; V A, V B 93 L x L y L x L y L x L y z for ξ micro z L ξ micro. This is physically reasonable, but in this subsection, we derive 93 from the microscopic description. First, we extend the microscopic momentum current Ĵ ab r; Γ to the system walls. Using 5 8, we rewrite the first term on the right-h side of as N ˆF Ar a i ; Γ A ab ĴA δr r i = r; Γ, Γ A r a Ĵ ab A r; Γ, Γ A = N N A U A r qj A j= r a N N A j= U A r q A j r a r=ri δr q A j 94 r=ri r b i q Ab j Dr; r i, q A j. 95 Similarly, we introduce Ĵ ab B r; Γ, Γ B rewrite the second term on the right-h side of. Substituting these equations into, we obtain ˆπ x r; Γ t + Ĵ xb r; Γ t t r b N N A U A r qj A = r a j= + xb ĴA r; Γ t r b + δr qj r=rit A xb ĴB r; Γ t r b N N B U B r qj B j= r a δr qj r=rit B. 96 As our concern is the fluid in the region < z < L, we may neglect the right-h side of 96. Equation 96 implies that Ĵ A abr; Γ, Γ A Ĵ B abr; Γ, Γ B correspond to the microscopic momentum current induced by walls A B, respectively. We now focus on the force balance equation in the steady state. Because ˆπ x r, t V A,V B is independent of time t in the steady state, 96 is rewritten as r b J xb r, τ ss ; V A, V B + JA xb r, τ ss ; V A, V B + JB xb r, τ ss ; V A, V B = 97 for < z < L, where J ab α r, t; V A, V B α = A, B is defined by J ab α r, t; V A, V B = Ĵ ab α r, t V A,V B. 98

17 Statistical mechanical expressions of slip length 7 Spatially averaging 97 in the xy plane at a position z gives dxdy J xz r, τ ss ; V A, V B + JA xz r, τ ss ; V A, V B + JB xz r, τ ss ; V A, V B = const 99 L x L y z in z, where we have used the periodic boundary conditions along the x y axes. By definition 95, ĴA abr; Γ, Γ A is zero in the region z ξ micro where the potential force of wall A cannot reach. Similarly, ĴB abr; Γ, Γ A is zero in the region z L ξ micro. Using these properties, 99 is rewritten as dxdyj xz r, τ ss ; V A, V B = const L x L y z in ξ micro z L ξ micro. Next, we introduce ξ A so that < ξ A ξ micro holds. Because fluid particles are seldom present in the thin layer < z < ξ A due to the repulsive force of wall A, it follows that Ĵ ab r; Γ in this layer < z < ξ A. Using Gauss s divergence theorem relation 95, we have dxdyĵ A xz r; Γ, Γ A L x L y z=ξ A = dxdydz L x L y z ξ A z Ĵ A xz r; Γ, Γ A = N N N A dxdydz ˆF a L x L Ar i ; Γ, Γ A δr r i + y z ξ A j= U A r qj A r=ri r a δr qj A. Similarly, as particles are seldom present in the region z < ξ A, the first term of is zero. Therefore, we obtain dxdyĵ A xz r; Γ, Γ A ˆF A xγ t; Γt A, 3 L x L y z=ξ A L x L y where we have used 8. From 99,,, 3, we find J xz z, τ ss ; V A, V B F x A τ ss; V A, V B L x L y 4 in ξ micro z L ξ micro. Finally, taking into account that v z z, τ ss ; V A, V B = in the steady state, we obtain the first equality of 93. By a similar procedure for wall B, we obtain the second equality of Galilean invariance of steady state The transport coefficients in 89, 9 9 are dependent on each other because of the Galilean invariance of the steady state. In this subsection, we derive relations between these transport coefficients. Let the probability density in the non-equilibrium steady state be P ss ri N, p i N ; V A, V B. The Galilean invariance in the steady state is expressed in terms of P ss ri N, p i N ; V A, V B as P ss ri N, p i N ; V A + V, V B + V = P ss ri N, p i N ; V A, V B, 5 where p i = px i + V, py i, pz i for any V. As the microscopic force acting on the fluid from wall A, 8, is independent of the momentum of the particles composing the fluid, F x A τ ss; V A + V, V B + V

18 8 Hiroyoshi Nakano, Shin-ichi Sasa is independent of V. In the linear response regime, this implies that F x A τ ss; V A, V B is proportional to the relative velocity between the walls, V A V B. Using this property 9, we find In the same manner, we obtain γ AA τ ss + γ AB τ ss =. 6 γ BA τ ss + γ BB τ ss =. 7 Next, because the microscopic momentum density ˆπ a r; Γ depends on the the momentum of the particles composing the fluid, ˆπ x r, τ ss V A+V,V B +V depends on V. Using 5, we obtain ˆπ x r, τ ss V A,V B = dγ = dγ N p x i δ r r i Pss ri N, p i N ; V A, V B N p x i V δ r r i Pss ri N, p i N ; V A + V, V B + V = ˆπ x r, τ ss V A+V,V B +V V ρr, τ ss ; V A + V, V B + V. 8 In particular, by focusing on the linear response regime, 8 is rewritten as From 89 9, we have ˆπ x r, τ ss V A+V,V B +V = ˆπ x r, τ ss V A,V B + V ρr, τ ss ; V A, V B. 9 γ P A r, τ ss + γ P B r, τ ss = ρr, τ ss ;,. Finally, by applying the Onsager reciprocity theorem to 9 9, we obtain γ AB τ ss = γ BA τ ss. In summary, we find that the linear response relations 89, 9 9 may be expressed using only two dependent transport coefficients. 6.4 Expression of the slip length The linear response theory derived in this section is consistent the stard hydrodynamics. Then, we assume that the macroscopic behavior of the fluid in the bulk is described as an incompressible Newtonian liquid. That is, we have ρr, t; V A, V B = ρ, σ xz r, t; V A, V B = η vx r, t; V A, V B z 3 in the bulk, where ρ is a constant value η is known. In this subsection, we derive the most general boundary condition under this assumption. From 93 3, the velocity field in the steady state is expressed as L x L y z dxdyv x r, τ ss ; V A, V B = F x A τ ss; V A, V B ηl x L y z z + L x L y z dxdyv x r, τ ss ; V A, V B,4 in the bulk, where z is any position in the bulk. Here, it is plausible that the velocity fields in the steady state, v x r, τ ss ; V A, V B, are independent of x, y. Then, we focus on the z-dependence of

19 Statistical mechanical expressions of slip length 9 v x r, τ ss ; V A, V B omit the x, y-dependence. By substituting the linear response relations 89, 9 9 into 4, we obtain v x z, τ ss ; V A, V B = γ AAτ ss V A + γ AB τ ss V B z z η + γ P Az, τ ss V A + γ P B z, τ ss V B. 5 ρ Using 6, 7,,, 5 is expressed in terms of two transport coefficients as v x z, τ ss ; V A, V B = γ AAτ ss V A V B z + η γ P Bz, τ ss ρ + γ AAτ ss z V A V B + V A. 6 η Here, we extrapolate the velocity field in the bulk, 6, to the whole region. Then, we obtain the boundary condition for the extrapolated velocity field: v x z =, τ ss ; V A, V B V A = b ss A v x z, τ ss ; V A, V B z 7 z= b ss A = ηγ P Bz, τ ss z. 8 ργ AA τ ss As long as we are interested in the velocity field in the bulk, we may use the boundary condition for the extrapolated velocity field 7. Because 8 is independent of v x z =, τ ss ; V A, V B, V A, V B, we identify 7 as the partial slip boundary condition constant slip length b ss A. Therefore, we confirm that the partial slip boundary condition is the most general boundary condition consistent the linear response theory in this section. We here show that the slip length b ss A is independent of z. Note that z in 8 comes from the reference point for the velocity field v x z, τ ss ; V A, V B, 4. Because v x z, τ ss ; V A, V B is independent of the reference point, we find that the slip length b ss A is also independent of z : db ss A dz =. 9 Finally, we similarly derive the partial slip boundary condition for wall B: v x z = L, τ ss ; V A, V B V B = b ss v x z, τ ss ; V A, V B B z z=l b ss B = ηγ P Az, τ ss L z. ργ BB τ ss 7 Connection between two linear response theories In Sects. 5 6, we developed two linear response theories that describe the macroscopic behavior belonging to different time scales. From these two theories, we obtained two microscopic expressions of the slip length: 8 8. In this section, we establish a connection between these linear response formulae using a linearized fluctuating hydrodynamics the partial slip boundary condition. By explicitly calculating the correlation function F x A sf x A eq from τ meso -scale to τ ss -scale, we demonstrate that the linearized fluctuating hydrodynamics the partial slip boundary condition is consistent the two linear response theories derived in Sect That is, we obtain two microscopic expressions of the slip length premised on linearized fluctuating hydrodynamics. Through

20 Hiroyoshi Nakano, Shin-ichi Sasa this calculation, we confirm that the two microscopic expressions, 8 8, provide different expressions of the slip length. The method of calculation in this section is essentially the same as that in [6], which discusses that the cross correlations between the rom stress tensor other fluctuating fields vanish. However, by directly calculating these correlations, we shall show that these correlations do not vanish the previous results are modified. 7. Linearized fluctuating hydrodynamics We study a linearized fluctuating hydrodynamics the partial slip boundary condition. The geometry of the system is the same as that studied in Sect.. The system is defined in the region [, L x ] [, L y ] [, L]. Let ṽ a r, t be a fluctuating velocity field in an incompressible fluid. The time evolution of ṽ a r, t is described by the incompressible Navier Stokes equation stochastic fluxes, ρ ṽa t + J ab r b =, where twice repeated indices are assumed to be summed over. The momentum flux tensor J ab r, t is given by where s ab r, t is the Gaussian rom stress tensor satisfying ṽ a J ab = pδ ab η r b + ṽb r a + s ab, 3 s ab r, t s cd r, t = k B T η [δ ac δ bd + δ ad δ bc ] 3 δ abδ cd δ 3 r r δt t. 4 The pressure term in 3 is used to enforce the incompressibility condition, ṽ a r, t r a =. 5 In this study, we ignore the nonlinear effect induced by the advection term. We impose periodic boundary conditions along the x y axes the partial slip boundary condition slip length B A B B at z = z = L respectively, specifically, ṽ x ṽ x r r V A = B A, 6 z= z z= V B ṽ x ṽ x r r = B B, 7 z=l z z=l where V A V B are the wall velocities. In this section, we focus on stationary walls. That is, we set to V A = V B =. Because we measure all the quantities on the macroscopic scales, we convert all the quantities to the dimensionless form by setting L = τ marco = k B T =. Here, τ macro is given by τ macro = L ρ η, 8 which represents the relaxation time of the incompressible fluid. These dimensionless quantities are expressed by putting the superscript asterisk. We note that the units in this section are different from those in Sec. 5.

21 Statistical mechanical expressions of slip length We remark on the ultraviolet cutoff of the fluctuating hydrodynamics. The ultraviolet cutoff is given by ξ meso, τ meso. In particular, τ meso in the hydrodynamic description is assumed to be τ meso = ξ mesoρ. 9 η Using an ultraviolet cutoff, we rewrite the properties of the stochastic flux, 4, in a more precise form s ab r, t s cd r, t = η [ δ ac δ bd + δ ad δ bc 3 δ abδ cd ]D 3 r r T t t. 3 where D 3d T s are given by D 3d { = /ζ 3, for d < ζ, =, for d ζ, 3 T s { = /ζ, for s < ζ, =, for s ζ. 3 Here, ζ is defined by 5. The functions D 3d T s correspond to regularized functions of δ 3 r δt. We assume that the slip lengths B A B B are Oξ micro. Through later discussions, this assumption is found to be inevitable when we compare the results obtained from the linearized fluctuating hydrodynamics those obtained from the underlying microscopic theory. 7. Calculated quantities Let F A t F B t be forces that the fluid exerts on walls A B at time t, respectively. F A t is expressed in terms of the momentum flux J ab r, t as F At = dx dy J xz r, t z =ζ = η z Ṽ x ζ, t + S xz ζ, t, 33 where Ṽ x z, t = dx dy ṽ x r, t, z 34 S ab z, t = dx dy s ab r, t. z 35 F B t is also written as the similar form. We note that F A t is defined by the surface integral at z = ζ, not z =. By the definition 4, Sxz - S xz correlation function is calculated as S xz z, t S xz z, t = η L xl yδ z z δ t t, 36 where δ z z is interpreted as the same form as 3. Furthermore, we define the time correlation functions: C αβt, t F αt F βt, 37

22 Hiroyoshi Nakano, Shin-ichi Sasa where α, β = A, B. Using the explicit form of F A t given in 33, 37 is rewritten as the sum of four terms, C AAt, t = η z Ṽ x ζ, t z Ṽ x ζ, t + S xz ζ, t S xz ζ, t η z Ṽ x ζ, t S xz ζ, t η z Ṽ x ζ, t S xz ζ, t. 38 In this section, we calculate CAA t, t in equilibrium using the Green function. For a function ft, we define its Fourier transform as fω = dte iωt ft, 39 where we use the notation fω for the Fourier transform of ft. Due to the time translational invariance in equilibrium, CAA t, t is expressed in the form dω CAAt, t = π C AAω e iω t t 4 C AAω = η z Ṽ x ζ, ω z Ṽ x ζ, ω + S xz ζ, ω S xz ζ, ω η z Ṽ x ζ, ω S xz ζ, ω η z Ṽ x ζ, ω S xz ζ, ω. 4 Hereafter, we perform the calculation of C AA ω in the frequency domain. 7.3 Expression of C AA ω in terms of Green function Integrating 3 over the xy plane yields the equation of motion for Ṽx z, t, In the frequency domain, 4 is expressed as z q Ṽ x = η z Sxz ρ Ṽx t = η z Ṽx + z Sxz q = iω ρ η. 44 The stochastic properties of S xz in the frequency domain are given by S xz z, ω S xz z, ω = η L xl yδ z z δ ω + ω. 45 The Green function of 43 is defined as the solution of the equation z q G z, z ; ω = η δ z z 46 boundary conditions G z =, z; ω = BA G z, z; ω z z, 47 = G z =, z; ω = BB G z, z; ω z z. 48 =

23 Statistical mechanical expressions of slip length 3 As discussed in Refs. [6,6], the general solutions of 46 are given by G z, z; ω = ge q z + ge q z + η q e q z z, 49 where g g are arbitrary constants. Using the boundary conditions, g g are determined to be g = η q q BA + q BBe q z q BA q BBe q z, 5 g = η q + q BA q BBe q z q BA q BBe q z 5 = + q B A + q B Be q q B A q B Be q. 5 Therefore, the Green function is expressed as The solution of 43 the partial slip boundary condition is expressed in terms of the Green function as Ṽ x z, ω = dz G z, z, ω z Sxz z, ω. 53 We note that the integration limits on the right-h side are z = z =, not z = ζ z = ζ. Using this expression for the fluctuating velocity fields 53, we may rewrite 4 in terms of the Green function. We sketch this calculation below; see Appendix A for technical aspects of the calculation. We consider the velocity-velocity correlation function. By using 53, we obtain = Ṽx z, ω Ṽx z, ω dz dz G z, z, ω Gz, z, ω z Sxz z, ω z Sxz z, ω. 54 Using the properties of the rom stress 45 performing an integration by parts yield Ṽx z, ω Ṽx z, ω = L xl yη dz dz Gz, z, ω Gz, z, ω z z δz z = L xl yη dz z G z, z; ω z G z, z; ω + L xl yη G z, ; ω z G z, ; ω Gz, ; ω z G z, ; ω + L xl yη δ G z, ; ω G z, ; ω + G z, ; ω G z, ; ω. 55 Here, by recalling the regularization of the delta function, we find that δ in the last line equals ζ. Then, the first term of 4 is expressed as z Ṽ x ζ, ω z Ṽ x ζ, ω = L xl yη dz z z G ζ, z; ω z z G ζ, z; ω + L xl yη z G ζ, ; ω z z G ζ, ; ω z G ζ, ; ω z z G ζ, ; ω + L xl yη ζ z G ζ, ; ω z G ζ, ; ω + z G ζ, ; ω z G ζ, ; ω. 56

24 4 Hiroyoshi Nakano, Shin-ichi Sasa Similar calculations apply to the remaining terms of 4 lead to z Ṽ x ζ, ω S xz ζ, ω + z Ṽ x ζ, ω S xz ζ, ω = L xl yη z z G ζ, ζ; ω + z z G ζ, ζ; ω. 57 Finally, by recalling 45, we obtain the full expression of C AA ω in terms of the Green function, C AA ω η L xl y = δ η z z G ζ, ζ; ω + z z G ζ, ζ; ω + η dz z z G ζ, z; ω z z G ζ, z; ω + η z G ζ, ; ω z z G ζ, ; ω z G ζ, ; ω z z G ζ, L; ω + η δ z G ζ, ; ω z G ζ, ; ω + z G ζ, ; ω z G ζ, ; ω,58 where δ = ζ. We note that the third fourth lines of 58 stem from the surface terms of the integration by parts. When we consider the infinite system or the stick boundary condition, these surface terms are zero. Having imposed the partial slip boundary condition on the walls, we cannot neglect these terms. 7.4 Behavior of C AA ω in low- high- frequency regions Using the explicit form of the Green function , we may obtain the full expression of C AA ω. However, this expression is rather complicated. Therefore, we focus on C AA ω in the low- high- frequency regions, which are defined as ω τmacro =, 59 ω τmeso = ζ, 6 respectively. These regions are simply expressed as ω ω. Before proceeding the calculation of C AA ω in the low- high- frequency regions, we remark on the full expression of C AA ω. We show that the first second lines of 58 do not contain terms proportional to δ. The explicit expression for the derivative of the Green function in these lines is z z Gζ, z ; ω = q e q η + q BBe q z q e q η q B B e q z + η δ ζ z. 6 From this expression, we find that the second third terms in the first line of 58 contain terms proportional to δ, η z z G ζ, ζ; ω + z z G ζ, ζ; ω = δ + other terms, 6 the second line of 58 contains terms proportional to δ, η dz z z G ζ, z; ω z z G ζ, z; ω = δ + other terms. 63 Substituting 6 63 into 58, we find that there is no term proportional to δ in the first second lines of 58. Here, we note that the second third terms in the first line of 58