The First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique
|
|
- Eugenia Dorsey
- 6 years ago
- Views:
Transcription
1 Commun. Theor. Phys. (Beijing, China 35 (2 pp c International Academic Publishers Vol. 35, No. 4, April 5, 2 The First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique TIAN Chu-Shun, YU Hao, ZHANG Chi 2 LU Quan-Kang,3 Department of Physics, Fudan University, Shanghai 2433, China 2 Department of Physics, Maryl University, College Park, 2742, USA 3 T.D. Lee Physics Laboratory, Fudan University, Shanghai 2433, China (Received November 8, 999; Revised April 4, 2 Abstract For weak-coupling electronic plasmas in the presence of small fluctuations, we calculate Green-Kubo formula for transport coefficients employing the ab initio theory, i.e., the two-time ensemble technique. The microscopic interpretations of Onsager s hypothesis on fluctuations are presented in the framework. PACS numbers: 5.4.+j, 5.6.+w Key words: two-time ensemble technique, Green Kubo formula Introduction It is well known that Green Kubo formula is essential in understing linear transport processes. However, its generality also arises great difficulties in practice. Since 95 s there have been a few methods of calculating it, but they are all not the ab initio theory. Among these is the one imposed by Chrasekhar then further developed by Mclennan, ] which is based on Brownian motion, or equivalently, linearized Fokker Planck equation (FPE. Another important method is the field theory formalism. 2] A great advantage is that collective modes are incorporated, thus it is capable of treating anomalous transport. 3] Nevertheless, the problem still remains that the stochastic dynamical model is in prior, so the method is phenomenological cannot yield the first principle calculation via Green Kubo formula. In many realistic cases like plasmas, it is interesting to underst transport processes through time correlation of great value to explore calculating transport coefficients through Green Kubo formula based on the first principle. It is indicated through the exact Green Kubo formula derived from Liouville equation 4] that generally transport is in the close connection with both the time correlation of the same particle the correlation of different particles. So we naturally investigate the time correlation by the statistics of Γ space. In this paper we generalize the method adopted in Ref. 5] to dissipative cases. For simplicity, we discuss the correlation function of the weak-coupling plasma, taking into account neutral ion background small static electric fluctuations caused by electronic motions. Then we calculate the transport coefficients employing the twotime ensemble technique. 5] The better deeper understing of the previous theories are presented. 2 The Two-Time Ensemble Technique The similar idea was illustrated by Balescu. 6] However, no further research has been pursued to our knowledge. In Ref. 5], we introduced the so-called two-time ensemble technique for analyzing the time correlation. Introduce the two-time ensemble distribution, say f(γ, Γ; t + τ, t, Γ = (, 2,, N Γ = (, 2,, N. We use k to denote arguments (x k, p k k to denote arguments (x k, p k. The well-known continuing equation the measure preserving yield ( N τ + v k x + N k k= k= The initial condition is F k p k f(γ, Γ; t+τ, t =. ( f(γ, Γ; t, t = f(γ; tδ( δ(n N. (2 F k is the force exerted on particle k at time (t + τ. The current-current correlation is j e (tj e (t + τ = = N j= N p j p kf(γ, Γ; t + τ, td dn d dn k= p p f (, ; t + τ, td d + p 2 p f (, 2; t + τ, td d2. (3 The project supported by National Natural Science Foundation of China under Grant No High-Tech Program of China
2 No. 4 The First Principle Calculation of Green Kubo Formula with 43 Here the distributions f (, ; t + τ, t f (, 2; t + τ, t have been introduced f (, ; t + τ, t = N f(γ, Γ; t + τ, td2 dn d2 dn, (4 f (, 2; t + τ, t = N(N The Green Kubo formula for the electric conductivity is expressed as 4,6] σ = β j e = V f(γ, Γ; t + τ, td2 dn dd3 dn. (5 dτ j e j e (τ eq, (6 N n= (Ze/mp n is the microscopic electric current density with V being the volume of the system. The problem remains unsolved how to calculate both the f (, ; t+τ, t f (, 2; t+τ, t. This will be performed in the following section. For the purpose of illustrating the technique, we are restricted to the case of weak-coupling Coulombian plasmas. 3 Calculating f (, ; t + τ, t Due to the similarity of Eq. ( to ordinary Liouville s equation, the subdynamic theory also applies. 5,6] The great advantage of the theory is that it provides a systematic technique of obtaining kinetic equation based on perturbation theory. Adopting the notations of Ref. 6] simply following the technique, we can exp Eq. ( up to the second order (i.e., the weak-coupling approximation in the interaction parameter λ as τ V f(γ, Γ; t + τ, t = V Γ V f(γ, Γ; t + τ, t = V L f(γ, Γ; t + τ, t + λv L f(γ, Γ; t + τ, t + λ 2 dsv L CU (sl U ( sv f(γ, Γ; t + s, t. (7 Since Γ t play only the role of parameters, we drop them from now on for convenience. Then in the realization of specific correlation pattern, equation ( can be written as τ f ( ; t + τ = L f ( ; t + τ + λ dx 2L f 2 ( ; t + τ f (2 ; t + τ + λ 2 ds dx L 2U 2(sL 2U 2( s f ( ; t + τ f (2 ; t + τ. (8 Before further processing, some words are in order. It is well known that the regression of macroscopic small fluctuation can be described through linear hydrodynamics, which is the nucleus of Onsager s hypothesis. 7] The technique here presents its understing at the microscopic level. That is, in our opinion, the kinetic equation resulting in the evolution of macroscopic fluctuation preserves the same form of the one resulting in hydrodynamics. For the dispersive case, it has been put into practice in the previous work. 5] The present paper is dedicated to exploring how to generalize it to the dissipative case. Taking the above discussions into account, we reduce Eq. (8 into the formal Lau equation straightforward τ f ( ; t + τ + v f ( ; t + τ ee m v f ( ; t + τ = B dv 2 2 r g δ rs g r g s g 3 2 s f (x, v ; t + τ f (x, v 2; t + τ (9 with the self-consistent fluctuating static electric field E satisfying E = 4πe dv f ( ; t + τ Zn + ], (
3 44 TIAN Chu-Shun, YU Hao, ZHANG Chi LU Quan-Kang Vol. 35 B = 8π 5 dkk 3 V 2 k with V k being the Fourier component of the interaction potential V (r, we also take into account neutral ion background. Combined with the initial condition f (, ; t, t = f (; tδ(, equations (9 ( yield f ( ; t + τ. Consider the fundamental solution of Eqs (9 (: G(x, v, x, v; t+τ, t. For simplicity, we rewrite Eq. (9 as KG = v G + B ee εâg = m v { G = 4πe2 x m τ G + εâg = KG, ( dv 2 2 r g δ rs g r g s g 3 2G(x s, v, x, v; t + τg(x, v 2, x, v; t + τ (2 dx x x } dv G(x, v, x, v ; t + τ Zn + ] v G(x, v, x, v ; t + τ. (3 Here ε is a small parameter being characteristic of the perturbation by the mean field. In Eq. (, KG can be regarded as the counterpart of the regression hypothesis of fluctuations. While ÂG implies the relation between dispersion fluctuation. 5] Generally KG, or collision, dominates so that the perturbative technique can be made used of. Exp the solution G with respect to ε, with G ( satisfying G = G ( + ε n G (n (4 n= τ G( = KG ( (5 G ( τ= = δ(x xδ(v v. (6 Below we will find that the linear transport results from linearized K, hence we exp Eq. ( up to the first order in ε (that is, fluctuations are small after linearizing K, obtain G = G ( + ε τ ds e (τ s ˆK ( ÂG( τ=s. (7 The initial condition G ( τ= = implied by Eq. (6 has been employed. K is linearized as K. Equivalently, equation (7 can be expressed as G(x, v, x, v; t + τ = G ( (x, v, x, v; t + τ + ε τ ds e (τ s K { } dudwδ(x uδ(v w ÂG( (u, w, x, v; t + s] τ = G ( (x, v, x, v; t + τ + ε ds dudw e (τ s ˆK {δ(x uδ(v w ÂG( (u, w, x, v; t + s]}. (8 In the above equation, the linear property of e τ ˆK has been used. e (τ s ˆK δ(x uδ(v w is the fundamental solution of linearized Eq. (5. From Eq. (2 we notice that Lau collision operator is a nonlinear Fokker Planck operator, 6] KG(x, v, x, v; t + τ = Bm 2 r A r (v + ] 2 s B rs (v ] G(x, v, x, v; t + τ, (9 A r (v = dv 2 ( s s 2 g δ rs g r g s g 3 ] G(x, v 2, x, v; t + τ, (2
4 No. 4 The First Principle Calculation of Green Kubo Formula with 45 B rs (v = 2 dv g δ rs g r g s 2 g 3 G(x, v 2, x, v; t + τ, (2 the dummy indices r, s run from to 3. Here a key issue must be addressed. That is, although the time in Green Kubo formula runs over the positive axis, it is reasonable to assume that the integration is dominated by the time of hydrodynamic scale since in the smaller scale, the regular correlation is not formed, consequently leading to the counteraction of transient correlation. ] Moreover, we believe that Onsager s hypothesis about the evolution of small fluctuation implies the great similarity between the kinetic equations describing fluctuation hydrodynamics. Since the linear Lau equation yields hydrodynamic modes, naturally we linearize Eqs (9 (2 also. Thus G(v 2 in Eqs (2 (2 can be substituted by the stationary solution of δ(x x δ(v v, i.e., the well-known Maxwell Boltzmann distribution ( m 3/2 G(v 2 = n e mv 2 /2kBT. (22 2πk B T With the approximation of the small angle head-on collision introduced, A B are further simplified as with ( 4n γ(v 2 ( mv = ( 2n α(v 2 ( mv = A r (v = γ(v v r, (23 B rs (v = α(v δ rs + β(v v r v s rs (v (24 ( Ω v 3 dx e (x mv /2k BT = mv /2k B T /2 ( Ω v ( 2n β(v = γ(v 2 ( mv + ( 2 2n 2 dxx e (x mv /2k BT = mv /2k B T 3/2 ( Ω v 3 π ( mv 2k B T ( Ω v 3, (25 ( + 2π n 2( mv /2 ( Ω v, (26 dx mv /2k B T x e (x mv /2k BT (27 I (v = I 2. (28 I 3 Here I, I 2 I 3 can be explicitly expressed as ( 2n 2 ( mv /2 ( I = I 2 = Ω ] x 4 e ( Ω x2 v dx (x + mv /2k BT 3 4π ( 2n 2 ( mv /2 ( I 3 = Ω ] x 4 e x2 v dx (x + mv /2k, (3 BT 3 (29 Ω is the solid angle within which the head-on collision occurs. Ω/4π is a small parameter. In the limit of Ω/4π, I I 2 are dropped out. Moreover, the numerical calculation shows that in the regime of near equilibrium, A r (v the first two terms of Eq. (24 are much greater than rs (v. Hence it is also dropped out. Notice that the v r v s (r s terms involves higher-order hydrodynamic behavior (nondiagonal momentum flow or shear modes, we can drop them in the linear approximation. Finally, equation (9 reaches the stard form of linear FPE, ] α = τ G + v G = γ v (v G + 2 α 2 v G, (3 γ = 4 2π Ω ( m 3/2n 2 Bm 2 (32 4π 2k B T 4( + 2π 4 2π + 4 π /2 π /2 dxx e (x 3] 3 Ω ( m /2n 2 Bm 2. 4π 2k B T (33
5 46 TIAN Chu-Shun, YU Hao, ZHANG Chi LU Quan-Kang Vol. 35 Chrasekhar obtained the exact fundamental solution of Eq. (3, G = 8π(F I H 2 exp 3/2 IR2 2HR S+F S 2 2(F I H 2 ], (34 F = ( α 2 γ 3 2γ t e γ t e 2γ t ], I = ( α 2 γ 2 e 2γ t ], H = ( α 2 γ 2 2 e γ t ] 2, R = x x v e γ t γ, S = v v e γ t. It should be clarified that the starting point here substantially differs from that of Ref. ], although the exact solution of Ref. ] is employed. The latter is based on the stochastic model describing Brownian motion. 4 The Transport Coefficients For normal transport, combining Eqs (3 (6 yields σ ij = β dτ dd p i p j f (, ; t + τ, t = β dτ dd p i p j e ˆL τ δ( f (; t, (35 e ˆL τ δ( is nothing but the fundamental solution of the linearized Eqs ( (3 up to the first order in ε. Thus taking Eq. (34 into it, σ ij is written as σ ij = β dτ dd p i p j G(x, p, x, p ;τ τ f (x, p ; t + ɛ ds G(x u, v, w; τ s dudw ( ÂG(u x, w, p ; sf (x, p ; t] = 2 δij β(ze α /γ 2 + ɛβ dτ τ ds d d dudwp i p j (G(x u, v, w ; τ s. (36 In the last line we employ the integral performed by Mclennan ] with α γ expressed in terms of microscopic variables. We find that the electric conductivity consists of two parts. The leading part is proportional to diffusion coefficient, ] as one may expect. The second part arises from the first order self-consistent field correction, which can be neglected in many cases. 5 Conclusions In this paper, we calculate Green Kubo formula for transport coefficients through the first principle. The stochastic model has been recovered at the kinetic level. The technique gives us deep insight into the microscopic meaning of Onsager s hypothesis on fluctuations. Indeed, the present calculations are restricted to specific cases. We take into account both the weak-coupling approximation small static electric fluctuation approximation. Besides, the neutral ion background is employed, implying that only electronic motions account for fluctuations. It is essential that the calculations are valid for equilibrium systems, however, we hope to explore whether it is possible to generalize the present technique to other cases in future work. References ] S. Chrasechar, Rev. Mod. Phys. 5 (943 ; J.A. Mclennan, Phys. Fluids 9 ( ] B.C. Hohenberg B.I. Halperlin, Rev. Mod. Phys. 49 ( ] Y.P. HUO J.R. ZHENG, Theory of Non-equilibrium Statistical Mechanics, Chap. 5, Science Press, Beijing (987. 4] R. Kubo, M. Toda N. Hashitsume, Statistical Physics II, 2nd. ed., Springer-Verlag, Heidelberg (992. 5] C.S. TIAN Q.K. LU, Phys. Rev. E59 ( ] R. Balescu, Equilibrium Non-equilibrium Statistical Mechanics, Wiley-Interscience, New York (975. 7] L. Onsager, Phys. Rev. 37 (93 45; ibid. 38 (
The dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
More informationCHAPTER V. Brownian motion. V.1 Langevin dynamics
CHAPTER V Brownian motion In this chapter, we study the very general paradigm provided by Brownian motion. Originally, this motion is that a heavy particle, called Brownian particle, immersed in a fluid
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationThis is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or
Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects
More informationHIGH FRICTION LIMIT OF THE KRAMERS EQUATION : THE MULTIPLE TIME SCALE APPROACH. Lydéric Bocquet
HIGH FRICTION LIMIT OF THE KRAMERS EQUATION : THE MULTIPLE TIME SCALE APPROACH Lydéric Bocquet arxiv:cond-mat/9605186v1 30 May 1996 Laboratoire de Physique, Ecole Normale Supérieure de Lyon (URA CNRS 1325),
More informationNPTEL
NPTEL Syllabus Nonequilibrium Statistical Mechanics - Video course COURSE OUTLINE Thermal fluctuations, Langevin dynamics, Brownian motion and diffusion, Fokker-Planck equations, linear response theory,
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 5, April 14, 2006 1 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
1 Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 5, April 14, 2006 1 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationMajor Concepts Lecture #11 Rigoberto Hernandez. TST & Transport 1
Major Concepts Onsager s Regression Hypothesis Relaxation of a perturbation Regression of fluctuations Fluctuation-Dissipation Theorem Proof of FDT & relation to Onsager s Regression Hypothesis Response
More informationQuantum measurement theory and micro-macro consistency in nonequilibrium statistical mechanics
Nagoya Winter Workshop on Quantum Information, Measurement, and Quantum Foundations (Nagoya, February 18-23, 2010) Quantum measurement theory and micro-macro consistency in nonequilibrium statistical mechanics
More informationRational Form Solitary Wave Solutions and Doubly Periodic Wave Solutions to (1+1)-Dimensional Dispersive Long Wave Equation
Commun. Theor. Phys. (Beijing, China) 43 (005) pp. 975 98 c International Academic Publishers Vol. 43, No. 6, June 15, 005 Rational Form Solitary Wave Solutions and Doubly Periodic Wave Solutions to (1+1)-Dimensional
More informationVIII.B Equilibrium Dynamics of a Field
VIII.B Equilibrium Dynamics of a Field The next step is to generalize the Langevin formalism to a collection of degrees of freedom, most conveniently described by a continuous field. Let us consider the
More informationThe Kramers problem and first passage times.
Chapter 8 The Kramers problem and first passage times. The Kramers problem is to find the rate at which a Brownian particle escapes from a potential well over a potential barrier. One method of attack
More informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate Theory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
More informationPlasma waves in the fluid picture I
Plasma waves in the fluid picture I Langmuir oscillations and waves Ion-acoustic waves Debye length Ordinary electromagnetic waves General wave equation General dispersion equation Dielectric response
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More informationF r (t) = 0, (4.2) F (t) = r= r a (t)
Chapter 4 Stochastic Equations 4.1 Langevin equation To explain the idea of introduction of the Langevin equation, let us consider the concrete example taken from surface physics: motion of an atom (adatom)
More information08. Brownian Motion. University of Rhode Island. Gerhard Müller University of Rhode Island,
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 1-19-215 8. Brownian Motion Gerhard Müller University of Rhode Island, gmuller@uri.edu Follow this
More information1 Phase Spaces and the Liouville Equation
Phase Spaces and the Liouville Equation emphasize the change of language from deterministic to probablistic description. Under the dynamics: ½ m vi = F i ẋ i = v i with initial data given. What is the
More informationBrownian Motion: Fokker-Planck Equation
Chapter 7 Brownian Motion: Fokker-Planck Equation The Fokker-Planck equation is the equation governing the time evolution of the probability density of the Brownian particla. It is a second order differential
More informationSymmetry of the Dielectric Tensor
Symmetry of the Dielectric Tensor Curtis R. Menyuk June 11, 2010 In this note, I derive the symmetry of the dielectric tensor in two ways. The derivations are taken from Landau and Lifshitz s Statistical
More informationPartial factorization of wave functions for a quantum dissipative system
PHYSICAL REVIEW E VOLUME 57, NUMBER 4 APRIL 1998 Partial factorization of wave functions for a quantum dissipative system C. P. Sun Institute of Theoretical Physics, Academia Sinica, Beiing 100080, China
More information14. Energy transport.
Phys780: Plasma Physics Lecture 14. Energy transport. 1 14. Energy transport. Chapman-Enskog theory. ([8], p.51-75) We derive macroscopic properties of plasma by calculating moments of the kinetic equation
More informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 2 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationNon equilibrium thermodynamic transformations. Giovanni Jona-Lasinio
Non equilibrium thermodynamic transformations Giovanni Jona-Lasinio Kyoto, July 29, 2013 1. PRELIMINARIES 2. RARE FLUCTUATIONS 3. THERMODYNAMIC TRANSFORMATIONS 1. PRELIMINARIES Over the last ten years,
More informationOn the Boltzmann equation: global solutions in one spatial dimension
On the Boltzmann equation: global solutions in one spatial dimension Department of Mathematics & Statistics Colloque de mathématiques de Montréal Centre de Recherches Mathématiques November 11, 2005 Collaborators
More informationThe Boltzmann Equation and Its Applications
Carlo Cercignani The Boltzmann Equation and Its Applications With 42 Illustrations Springer-Verlag New York Berlin Heidelberg London Paris Tokyo CONTENTS PREFACE vii I. BASIC PRINCIPLES OF THE KINETIC
More informationCritical Phenomena under Shear Flow
Critical Phenomena under Shear Flow Pavlik Lettinga, Hao Wang, Jan K.G. Dhont Close to a gas-liquid critical point, effective interactions between particles become very long ranged, and the dynamics of
More information16. Working with the Langevin and Fokker-Planck equations
16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation
More informationFokker-Planck Equation with Detailed Balance
Appendix E Fokker-Planck Equation with Detailed Balance A stochastic process is simply a function of two variables, one is the time, the other is a stochastic variable X, defined by specifying: a: the
More informationNon equilibrium thermodynamics: foundations, scope, and extension to the meso scale. Miguel Rubi
Non equilibrium thermodynamics: foundations, scope, and extension to the meso scale Miguel Rubi References S.R. de Groot and P. Mazur, Non equilibrium Thermodynamics, Dover, New York, 1984 J.M. Vilar and
More informationLinear Response and Onsager Reciprocal Relations
Linear Response and Onsager Reciprocal Relations Amir Bar January 1, 013 Based on Kittel, Elementary statistical physics, chapters 33-34; Kubo,Toda and Hashitsume, Statistical Physics II, chapter 1; and
More informationNoise, AFMs, and Nanomechanical Biosensors
Noise, AFMs, and Nanomechanical Biosensors: Lancaster University, November, 2005 1 Noise, AFMs, and Nanomechanical Biosensors with Mark Paul (Virginia Tech), and the Caltech BioNEMS Collaboration Support:
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 15 Consider Hamilton s equations in the form I. CLASSICAL LINEAR RESPONSE THEORY q i = H p i ṗ i = H q i We noted early in the course that an ensemble
More informationarxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
More informationarxiv: v1 [cond-mat.stat-mech] 26 Jan 2008
Long Time Tail of the Velocity Autocorrelation Function in a Two-Dimensional Moderately Dense Hard Disk Fluid Masaharu Isobe Graduate School of Engineering, Nagoya Institute of Technology, Nagoya 466-8555,
More informationPlasma Astrophysics Chapter 1: Basic Concepts of Plasma. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University
Plasma Astrophysics Chapter 1: Basic Concepts of Plasma Yosuke Mizuno Institute of Astronomy National Tsing-Hua University What is a Plasma? A plasma is a quasi-neutral gas consisting of positive and negative
More informationMagnetic waves in a two-component model of galactic dynamo: metastability and stochastic generation
Center for Turbulence Research Annual Research Briefs 006 363 Magnetic waves in a two-component model of galactic dynamo: metastability and stochastic generation By S. Fedotov AND S. Abarzhi 1. Motivation
More informationA New Integrable Couplings of Classical-Boussinesq Hierarchy with Self-Consistent Sources
Commun. Theor. Phys. Beijing, China 54 21 pp. 1 6 c Chinese Physical Society and IOP Publishing Ltd Vol. 54, No. 1, July 15, 21 A New Integrable Couplings of Classical-Boussinesq Hierarchy with Self-Consistent
More information26. Non-linear effects in plasma
Phys780: Plasma Physics Lecture 26. Non-linear effects. Collisionless shocks.. 1 26. Non-linear effects in plasma Collisionless shocks ([1], p.405-421, [6], p.237-245, 249-254; [4], p.429-440) Collisionless
More informationONSAGER S RECIPROCAL RELATIONS AND SOME BASIC LAWS
Journal of Computational and Applied Mechanics, Vol. 5., No. 1., (2004), pp. 157 163 ONSAGER S RECIPROCAL RELATIONS AND SOME BASIC LAWS József Verhás Department of Chemical Physics, Budapest University
More informationElectron-Vibrational Energy Exchange in Nitrogen-Containing Plasma: a Comparison Between an Analytical Approach and a Kinetic Model
Electron-Vibrational Energy Exchange in Nitrogen-Containing Plasma: a Comparison Between an Analytical Approach and a Kinetic Model YANG Wei ( ) and DONG Zhiwei ( ) Institute of Applied Physics and Computational
More informationNon-equilibrium phenomena and fluctuation relations
Non-equilibrium phenomena and fluctuation relations Lamberto Rondoni Politecnico di Torino Beijing 16 March 2012 http://www.rarenoise.lnl.infn.it/ Outline 1 Background: Local Thermodyamic Equilibrium 2
More informationMicroscopic Deterministic Dynamics and Persistence Exponent arxiv:cond-mat/ v1 [cond-mat.stat-mech] 22 Sep 1999
Microscopic Deterministic Dynamics and Persistence Exponent arxiv:cond-mat/9909323v1 [cond-mat.stat-mech] 22 Sep 1999 B. Zheng FB Physik, Universität Halle, 06099 Halle, Germany Abstract Numerically we
More informationStochastic continuity equation and related processes
Stochastic continuity equation and related processes Gabriele Bassi c Armando Bazzani a Helmut Mais b Giorgio Turchetti a a Dept. of Physics Univ. of Bologna, INFN Sezione di Bologna, ITALY b DESY, Hamburg,
More informationFermi Liquid and BCS Phase Transition
Fermi Liquid and BCS Phase Transition Yu, Zhenhua November 2, 25 Abstract Landau fermi liquid theory is introduced as a successful theory describing the low energy properties of most fermi systems. Besides
More informationUNIVERSITY OF OSLO FACULTY OF MATHEMATICS AND NATURAL SCIENCES
UNIVERSITY OF OSLO FCULTY OF MTHEMTICS ND NTURL SCIENCES Exam in: FYS430, Statistical Mechanics Day of exam: Jun.6. 203 Problem :. The relative fluctuations in an extensive quantity, like the energy, depends
More informationLow Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model
Low Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model G. A. Radtke, N. G. Hadjiconstantinou and W. Wagner Massachusetts Institute of Technology,
More informationMathematical Structures of Statistical Mechanics: from equilibrium to nonequilibrium and beyond Hao Ge
Mathematical Structures of Statistical Mechanics: from equilibrium to nonequilibrium and beyond Hao Ge Beijing International Center for Mathematical Research and Biodynamic Optical Imaging Center Peking
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationV. Electrostatics Lecture 24: Diffuse Charge in Electrolytes
V. Electrostatics Lecture 24: Diffuse Charge in Electrolytes MIT Student 1. Poisson-Nernst-Planck Equations The Nernst-Planck Equation is a conservation of mass equation that describes the influence of
More informationTitle of communication, titles not fitting in one line will break automatically
Title of communication titles not fitting in one line will break automatically First Author Second Author 2 Department University City Country 2 Other Institute City Country Abstract If you want to add
More informationA Realization of Yangian and Its Applications to the Bi-spin System in an External Magnetic Field
Commun. Theor. Phys. Beijing, China) 39 003) pp. 1 5 c International Academic Publishers Vol. 39, No. 1, January 15, 003 A Realization of Yangian and Its Applications to the Bi-spin System in an External
More informationNonlinear Analysis: Modelling and Control, Vilnius, IMI, 1998, No 3 KINK-EXCITATION OF N-SYSTEM UNDER SPATIO -TEMPORAL NOISE. R.
Nonlinear Analysis: Modelling and Control Vilnius IMI 1998 No 3 KINK-EXCITATION OF N-SYSTEM UNDER SPATIO -TEMPORAL NOISE R. Bakanas Semiconductor Physics Institute Go štauto 11 6 Vilnius Lithuania Vilnius
More informationEvaluation of the rate constant and deposition velocity for the escape of Brownian particles over potential barriers
arxiv:1411.0692v1 [cond-mat.stat-mech] 22 Oct 2014 Evaluation of the rate constant and deposition velocity for the escape of Brownian particles over potential barriers Michael W Reeks School of Mechanical
More informationLectures on basic plasma physics: Kinetic approach
Lectures on basic plasma physics: Kinetic approach Department of applied physics, Aalto University April 30, 2014 Motivation Layout 1 Motivation 2 Boltzmann equation (a nasty bastard) 3 Vlasov equation
More informationNumerical Simulation of the Rarefied Gas Flow through a Short Channel into a Vacuum
Numerical Simulation of the Rarefied Gas Flow through a Short Channel into a Vacuum Oleg Sazhin Ural State University, Lenin av.5, 6283 Ekaterinburg, Russia E-mail: oleg.sazhin@uralmail.com Abstract. The
More informationBROWNIAN DYNAMICS SIMULATIONS WITH HYDRODYNAMICS. Abstract
BROWNIAN DYNAMICS SIMULATIONS WITH HYDRODYNAMICS Juan J. Cerdà 1 1 Institut für Computerphysik, Pfaffenwaldring 27, Universität Stuttgart, 70569 Stuttgart, Germany. (Dated: July 21, 2009) Abstract - 1
More informationBasic Concepts and Tools in Statistical Physics
Chapter 1 Basic Concepts and Tools in Statistical Physics 1.1 Introduction Statistical mechanics provides general methods to study properties of systems composed of a large number of particles. It establishes
More informationKinetic Models for Granular Flow
Journal of Statistical Physics, Vol. 97, Nos. 12, 1999 Kinetic Models for Granular Flow J. Javier Brey, 1 James W. Dufty, 2 and Andre s Santos 3 Received September 28, 1998; final March 29, 1999 The generalization
More informationSolution of time-dependent Boltzmann equation for electrons in non-thermal plasma
Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma Z. Bonaventura, D. Trunec Department of Physical Electronics Faculty of Science Masaryk University Kotlářská 2, Brno, CZ-61137,
More informationThis work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 4.0 License.
University of Rhode Island DigitalCommons@URI Equilibrium Statistical Physics Physics Course Materials 2015 07. Kinetic Theory I Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons
More informationHilbert Sixth Problem
Academia Sinica, Taiwan Stanford University Newton Institute, September 28, 2010 : Mathematical Treatment of the Axioms of Physics. The investigations on the foundations of geometry suggest the problem:
More informationF(t) equilibrium under H 0
Physics 17b: Statistical Mechanics Linear Response Theory Useful references are Callen and Greene [1], and Chandler [], chapter 16. Task To calculate the change in a measurement B t) due to the application
More information9.1 System in contact with a heat reservoir
Chapter 9 Canonical ensemble 9. System in contact with a heat reservoir We consider a small system A characterized by E, V and N in thermal interaction with a heat reservoir A 2 characterized by E 2, V
More informationLandau Bogolubov Energy Spectrum of Superconductors
Landau Bogolubov Energy Spectrum of Superconductors L.N. Tsintsadze 1 and N.L. Tsintsadze 1,2 1. Department of Plasma Physics, E. Andronikashvili Institute of Physics, Tbilisi 0128, Georgia 2. Faculty
More informationLarmor radius effects on impurity transport in turbulent plasmas
Larmor radius effects on impurity transport in turbulent plasmas arxiv:physics/0408029v1 [physics.plasm-ph] 5 Aug 2004 M. Vlad and F. Spineanu National Institute for Laser, Plasma and Radiation Physics,
More informationStochastic equations for thermodynamics
J. Chem. Soc., Faraday Trans. 93 (1997) 1751-1753 [arxiv 1503.09171] Stochastic equations for thermodynamics Roumen Tsekov Department of Physical Chemistry, University of Sofia, 1164 Sofia, ulgaria The
More informationActive Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: Hydrodynamics of SP Hard Rods
Active Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: of SP Hard Rods M. Cristina Marchetti Syracuse University Baskaran & MCM, PRE 77 (2008);
More informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
More informationLinear-response theory and the fluctuation-dissipation theorem: An executive summary
Notes prepared for the 27 Summer School at Søminestationen, Holbæk, July 1-8 Linear-response theory and the fluctuation-dissipation theorem: An executive summary Jeppe C. Dyre DNRF centre Glass and Time,
More informationLANGEVIN EQUATION AND THERMODYNAMICS
LANGEVIN EQUATION AND THERMODYNAMICS RELATING STOCHASTIC DYNAMICS WITH THERMODYNAMIC LAWS November 10, 2017 1 / 20 MOTIVATION There are at least three levels of description of classical dynamics: thermodynamic,
More informationHypocoercivity for kinetic equations with linear relaxation terms
Hypocoercivity for kinetic equations with linear relaxation terms Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (A JOINT
More informationEquilibration and decoupling of a relativistic gas in a Friedmann-Robertson-Walker spacetime
Equilibration and decoupling of a relativistic gas in a Friedmann-Robertson-Walker spacetime Juan M. Torres-Rincon (Frankfurt Institute for Advanced Studies) in collaboration with J. Tindall, J.-B. Rosé,
More informationHierarchical Modeling of Complicated Systems
Hierarchical Modeling of Complicated Systems C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park, MD lvrmr@math.umd.edu presented
More informationResearch of Micro-Rectangular-Channel Flow Based on Lattice Boltzmann Method
Research Journal of Applied Sciences, Engineering and Technology 6(14): 50-55, 013 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scientific Organization, 013 Submitted: November 08, 01 Accepted: December 8,
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More informationMonte Carlo Collisions in Particle in Cell simulations
Monte Carlo Collisions in Particle in Cell simulations Konstantin Matyash, Ralf Schneider HGF-Junior research group COMAS : Study of effects on materials in contact with plasma, either with fusion or low-temperature
More informationTopics in Nonequilibrium Physics. Nicolas Borghini
Topics in Nonequilibrium Physics Nicolas Borghini Version of September 8, 2016 Nicolas Borghini Universität Bielefeld, Fakultät für Physik Homepage: http://www.physik.uni-bielefeld.de/~borghini/ Email:
More informationDissipative nuclear dynamics
Dissipative nuclear dynamics Curso de Reacciones Nucleares Programa Inter universitario de Fisica Nuclear Universidad de Santiago de Compostela March 2009 Karl Heinz Schmidt Collective dynamical properties
More informationChapter 2. Deriving the Vlasov Equation From the Klimontovich Equation 19. Deriving the Vlasov Equation From the Klimontovich Equation
Chapter 2. Deriving the Vlasov Equation From the Klimontovich Equation 19 Chapter 2. Deriving the Vlasov Equation From the Klimontovich Equation Topics or concepts to learn in Chapter 2: 1. The microscopic
More informationCollaborators: Aleksas Mazeliauskas (Heidelberg) & Derek Teaney (Stony Brook) Refs: , /25
2017 8 28 30 @ Collaborators: Aleksas Mazeliauskas (Heidelberg) & Derek Teaney (Stony Brook) Refs: 1606.07742, 1708.05657 1/25 1. Introduction 2/25 Ultra-relativistic heavy-ion collisions and the Bjorken
More informationNew ideas in the non-equilibrium statistical physics and the micro approach to transportation flows
New ideas in the non-equilibrium statistical physics and the micro approach to transportation flows Plenary talk on the conference Stochastic and Analytic Methods in Mathematical Physics, Yerevan, Armenia,
More informationKinetic theory. Collective behaviour of large systems Statistical basis for the ideal gas equation Deviations from ideality
Kinetic theory Collective behaviour of large systems Statistical basis for the ideal gas equation Deviations from ideality Learning objectives Describe physical basis for the kinetic theory of gases Describe
More informationFractional Quantum Mechanics and Lévy Path Integrals
arxiv:hep-ph/9910419v2 22 Oct 1999 Fractional Quantum Mechanics and Lévy Path Integrals Nikolai Laskin Isotrace Laboratory, University of Toronto 60 St. George Street, Toronto, ON M5S 1A7 Canada Abstract
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 6 Feb 2004
A statistical model with a standard Gamma distribution Marco Patriarca and Kimmo Kaski Laboratory of Computational Engineering, Helsinki University of Technology, P.O.Box 923, 25 HUT, Finland arxiv:cond-mat/422v
More informationLifshitz Hydrodynamics
Lifshitz Hydrodynamics Yaron Oz (Tel-Aviv University) With Carlos Hoyos and Bom Soo Kim, arxiv:1304.7481 Outline Introduction and Summary Lifshitz Hydrodynamics Strange Metals Open Problems Strange Metals
More informationDynamique d un gaz de sphères dures et équation de Boltzmann
Dynamique d un gaz de sphères dures et équation de Boltzmann Thierry Bodineau Travaux avec Isabelle Gallagher, Laure Saint-Raymond Outline. Linear Boltzmann equation & Brownian motion Linearized Boltzmann
More informationUne approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck
Une approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (EN
More informationAn Improved F-Expansion Method and Its Application to Coupled Drinfel d Sokolov Wilson Equation
Commun. Theor. Phys. (Beijing, China) 50 (008) pp. 309 314 c Chinese Physical Society Vol. 50, No., August 15, 008 An Improved F-Expansion Method and Its Application to Coupled Drinfel d Sokolov Wilson
More informationLattice Boltzmann Method
3 Lattice Boltzmann Method 3.1 Introduction The lattice Boltzmann method is a discrete computational method based upon the lattice gas automata - a simplified, fictitious molecular model. It consists of
More informationStudy of Pre-equilibrium Fission Based on Diffusion Model
Commun. Theor. Phys. (Beijing, China) 45 (2006) pp. 325 331 c International Academic Publishers Vol. 45, No. 2, February 15, 2006 Study of Pre-equilibrium Fission Based on Diffusion Model SUN Xiao-Jun
More informationStatistical Properties of a Ring Laser with Injected Signal and Backscattering
Commun. Theor. Phys. (Beijing, China) 35 (2001) pp. 87 92 c International Academic Publishers Vol. 35, No. 1, January 15, 2001 Statistical Properties of a Ring Laser with Injected Signal and Backscattering
More informationDr. A A Mamun Professor of Physics Jahangirnagar University Dhaka, Bangladesh
SOLITARY AND SHOCK WAVES IN DUSTY PLASMAS Dr. A A Mamun Professor of Physics Jahangirnagar University Dhaka, Bangladesh OUTLINE Introduction Static Dust: DIA Waves DIA Solitary Waves DIA Shock Waves Mobile
More informationStatistical Mechanics
Franz Schwabl Statistical Mechanics Translated by William Brewer Second Edition With 202 Figures, 26 Tables, and 195 Problems 4u Springer Table of Contents 1. Basic Principles 1 1.1 Introduction 1 1.2
More information