Physics 212: Statistical mechanics II Lecture IV
|
|
- Dulcie Burke
- 6 years ago
- Views:
Transcription
1 Physics 22: Statistical mechanics II Lecture IV Our program for kinetic theory in the last lecture and this lecture can be expressed in the following series of steps, from most exact and general to most specific and approximate (but also useful!): Liouville s theorem BBGKY hierarchy Boltzmann eqn. Navier-Stokes eqn. () The first step, from Liouville s theorem to the BBGKY hierarchy, will be completed in this lecture. The only assumption is that the system being dealt with is a gas with binary collisions, and the BBGKY hierarchy is still exact for this specific system. We will say a bit about the approximations required to get from the BBGKY hierarchy to the closed Boltzmann equation for f toward the end of this lecture, and start deriving the Navier-Stokes equation in the next lecture. Let s review what we learned last time. Suppose we have a problem of N particles moving in the Hamiltonian (switching back to the usual variables (r, p) rather than the generalized conjugate variables (q, p)) p 2 H = 2m + N U i + v ij i<j U i = U(r i ), v ij = v ji = v( r i r j ). (2) One of the main results from last lecture is Liouville s equation for the time evolution of g, the full N-particle distribution normalized to unity: Nd g t + ( ) g g ṗ i + ṙ i = 0. (3) p i r i For the specific case of the dilute gas with binary collisions, Liouville s equation becomes (where F = U, K ij = ri v ij ) g N t = F i K ij pi g p i g. (4) j=,j i Here g is the distribution function in N-particle phase space. For compactness, in most of this lecture we will suppress the time and sometimes other variables appearing in g, when the meaning is clear. Now rewrite Liouville s equation in the form [ t + h N(r, p,..., r N, p N )] g = 0, (5) where the differential operator h N is defined as h N (r, p,..., r N, p N ) = [ ] pi + F i pi
2 + 2 i,j= K ij ( pi pj ). (6) Our goal for most of this lecture will be to find an exact equation for the evolution of the n-body distribution functions f n by doing a partial integration of Liouville s equation. The definition of the -body distribution function f (z ), with z = (p, r), is f(z ) = δ(z z i ) = N g(z,..., z n )δ(z z ) dz... dz N = N g(z,..., z n ) dz 2... dz N. (7) Here the normalization factor N arose because the sum over i in the second form contains N identical terms. That is, the one-body distribution function f is a particle density (normalized to N), while the function g is normalized to over all of phase space. g differs from f N by a normalization factor N!. This choice of normalization may seem strange, but it will make the BBGKY equation derived below simpler. is Now the general definition of the n-body distribution function f n (z,..., z n ), with z = (p, r), f n (z,..., z n) = = j=,j i N! (N n)!... δ(z z i )δ(z 2 z j )... g(z,..., z n, z n+,..., z N ) dz n+... dz N, (8) where, as before for f, the combinatorial factor N!/(N n)! arises because each term in the sums over i, j,... contributes an equal result, given by the integral in the final form. A simple way to remember this normalization is that, just as f is normalized to give a particle density, f 2 is normalized to give a density of ordered pairs, f 3 of ordered triples, etc. as Now we show a surprising property of the differential operator h n, which is defined for n N Recall that h n (r, p,..., r n, p n ) = [ ] pi + F i pi + K ij ( pi pj ). (9) 2 i,j= ( ) t + h N g = 0. (0) Since the normalization of g is constant, integrating g t over all the space and momentum variables should give zero. Then, from (??), we also have h N (z,..., z N ) dz... dz N = 0. () This actually holds for all the h n for n < N as well: h n (z,..., z n ) dz... dz n = 0. (2) 2
3 To see this, go back to the definition of h n. It contains one set of terms which are spatial gradients of g multipled by space-independent terms, and another set of terms which are momentum gradients of g multipled by momentum-independent terms. To see that the first set vanish, consider holding the momentum variables constant and integrating the spatial variables. This integral over the gradient of g evaluates to g on the boundary of spatial volume, which for a bounded system we can take to be zero. Similarly the momentum-gradient terms can be argued to integrate to zero, under the assumption of a bounded system so boundary terms vanish. The operator h N can be broken down into terms containing only the first n particles, terms containing only the later N n particles, and terms which contain one particle from the first group and one particle from the second group: h N (z,..., z N ) = h n (z,..., z n ) + h N n (z n+,..., z N ) + j=n+ K ij ( pi pj ). (3) The factor of 2 in (??) disappears in the sum here because we no longer have both possibilities i < j and i > j. Now integrate out the last N n variables in the evolution of g: then the h N n part integrates to zero, by the lemma (??) shown above. Finally, t + h n(z,..., z n ) + j=n+ K ij ( pi pj ) g dz n+... dz N = 0. (4) The p j -gradient terms in the above integrate to zero for a bounded system as their coefficient is momentum-independent (same argument as above). Multiplying by N!/(N n)! then gives t + h n(z,..., z n )) f n (z,..., z n ) = N! (N n)! n j=n+ n K ij pi g dz n+... dz N N! = K i,n+ pi g dz n+... dz N (N n + )! n = K i,n+ pi f n+ (z,..., z n+ ) dz n+. (5) Here in passing from the first line to the second we used the fact that the sum over j contains N n equal terms, because of the integration. (We could not similarly get rid of the sum over i because the z i variables aren t integrated.) It is essential that the Hamiltonian contains only binary collisions (no three-body or higher interactions) so that the evolution of f at time t is determined by f and f 2 at that t, not by any higher f i. The expression (??) for the time evolution of f n in terms of f n+ is known as the BBGKY hierarchy. The key assumption that allowed us to write the hierarchy in this simple way was that collisions involve only two particles. Now we can say a bit about the physical meaning of the hierarchy. The first equation of the hierarchy is + F p f (t, z ) = K 2 p f 2 (t, z, z 2 ) dz 2. (6) 3
4 Here we have restored time as a variable in order to connect to the Boltzmann equation derived previously. The left side of the first hierarchy equation is just that of the Boltzmann equation, so somehow the right side must contain information about collisional processes that increase the entropy S of f. The second equation of the hierarchy is t + p + F p + p 2 m x 2 + F 2 p2 + ) 2 K 2 ( p p2 ) f 2 (t, z, z 2 ) = (K 3 p + K 23 p2 )f 3 (t, z, z 2, z 3 ) dz 3. (7) There is an important difference between the left sides of the first and second hierarchy equations. In the first hierarchy equation, there are no terms involving collisions, only streaming terms resulting from single-particle motion. These streaming terms, as seen before, cannot change the entropy, so important physics is lost if we neglect the collisional term on the right side of the first equation. In the second equation, however, there is already a collisional term on the left side. In fact, a simple estimate shows that for the dilute gas with λ = na 3, the left side collisional term K 2 dominates the right side collisional terms. Hence it is a good approximation for the dilute gas to truncate the hierarchy at the second equation, t + p + F p + p 2 m x 2 + F 2 p2 + ) 2 K 2 ( p p2 ) f 2 (t, z, z 2 ) = 0. (8) This is the first of several approximations involved in deriving the Boltzmann equation from the BBGKY hierarchy: it just amounts to neglect of any correlations beyond f 2. Some of the correlations in f 2 will also have to be neglected in order to obtain a closed equation for f. The second assumption will be that, on a relatively rapid time scale, f 2 equilibrates to be nearly a product distribution f 2 (t, z, z 2 ) f (t, z )f (t, z 2 ) except for rare collision events. Then f 2 continues to evolve in time, but its evolution is due to changes in f, rather than to changes in the pair correlations. Under this assumption f 2 t 0 and the second hierarchy equation can be used to simplify the first to + F p f (t, z ) = (p x + p 2 x2 )f 2 (t, z, z 2 ) dz 2. (9) m The remaining part of obtaining the Boltzmann equation requires working through the kinematics of collisions in some detail, and can be found in Huang and other textbooks (part of this may be an exercise on problem set 2). Instead we will focus in this lecture and the next on understanding more physical consequences of the Boltzmann equation. By asking the question of what distributions make zero the collisional term of the Boltzmann equation, we were led to the equation f f(t, r, p ) = exp(c 0 + c p + c 2 p 2 ). (20) Here and for the rest of this lecture we revert to the Boltzmann equation notation, where f, f 2, f, f 2 are shorthand for the one-particle distribution function evaluated at four different points in momentum space and one position. The above form for the distribution function came about from the 4
5 requirement that log f be a sum of quantities conserved in collisions: the collision term will vanish if for all incoming and outgoing momenta with we have f f 2 = f f 2, or p + p 2 = p + p 2, (2) log f(p ) + log f(p 2 ) = log f(p ) + log f(p 2). (22) Understanding conservation laws more generally will lead to the Navier-Stokes equation and other fundamental equations of hydrodynamics. We will then study a few examples of hydrodynamics such as Stokes law (used later for Brownian motion) and turbulence. Let be a conserved quantity like energy or momentum that moves with a particle, and define χ = χ(r, p ) as the value of χ for particle, etc. Then for any collision, χ + χ 2 = χ + χ 2. (23) We now show a local conservation law for such quantities: the collisional term in the Boltzmann equation does not change them, so ( ) f dp = 0. (24) t coll To show this, we make use of the same variable-switching tricks used in the proof of the H theorem in Lecture II to write the desired expression in a more symmetric form: ( ) f dp = t coll 4 w(f f 2 f f 2 )(χ + χ 2 χ χ 2) dp dp 2 dp dp 2 = 0. (25) Note that we did not need to require that the distribution f satisfied the Boltzmann transport equation. (In the derivation of the H theorem, we had log f in place of χ and so on.) Now integrating both sides of the Boltzmann equation over momentum (but not position) and using the above gives m x + F p f(t, r, p) dp = 0. (26) In the next lecture we will discuss physical interpretation of this law for conserved quantities such as energy and momentum. 5
G : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 1 Defining statistical mechanics: Statistical Mechanics provies the connection between microscopic motion of individual atoms of matter and macroscopically
More informationFrom Hamiltonian particle systems to Kinetic equations
From Hamiltonian particle systems to Kinetic equations Università di Roma, La Sapienza WIAS, Berlin, February 2012 Plan of the lectures 1 Particle systems and BBKGY hierarchy (the paradigm of the Kinetic
More informationTime-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics
Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle
More informationCHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,
More informationCHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation:
CHAPTER 9 Microscopic Approach: from Boltzmann to Navier-Stokes In the previous chapter we derived the closed Boltzmann equation: df dt = f +{f,h} = I[f] where I[f] is the collision integral, and we have
More informationp 2 V ( r i )+ i<j i + i=1
2. Kinetic Theory The purpose of this section is to lay down the foundations of kinetic theory, starting from the Hamiltonian description of 0 23 particles, and ending with the Navier-Stokes equation of
More informationLecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions
Lecture Models for heavy-ion collisions (Part III: transport models SS06: Dynamical models for relativistic heavy-ion collisions Quantum mechanical description of the many-body system Dynamics of heavy-ion
More informationAST1100 Lecture Notes
AST11 Lecture Notes Part 1G Quantum gases Questions to ponder before the lecture 1. If hydrostatic equilibrium in a star is lost and gravitational forces are stronger than pressure, what will happen with
More information7 Relic particles from the early universe
7 Relic particles from the early universe 7.1 Neutrino density today (14 December 2009) We have now collected the ingredients required to calculate the density of relic particles surviving from the early
More informationElectrical Transport in Nanoscale Systems
Electrical Transport in Nanoscale Systems Description This book provides an in-depth description of transport phenomena relevant to systems of nanoscale dimensions. The different viewpoints and theoretical
More informationUniversità di Pisa Facoltà di Scienze Matematiche Fisiche e Naturali
Università di Pisa Facoltà di Scienze Matematiche Fisiche e Naturali Corso di Laurea Specialistica in Fisica della Materia Anno Accademico 2005/2006 Tesi di Laurea Contributions to new generation ultra-cold
More information8.333: Statistical Mechanics I Fall 2007 Test 2 Review Problems
8.333: Statistical Mechanics I Fall 007 Test Review Problems The second in-class test will take place on Wednesday 10/4/07 from :30 to 4:00 pm. There will be a recitation with test review on Monday 10//07.
More informationChapter 10 Momentum, System of Particles, and Conservation of Momentum
Chapter 10 Momentum, System of Particles, and Conservation of Momentum 10.1 Introduction... 1 10.2 Momentum (Quantity of Motion) and Average Impulse... 1 Example 10.1 Impulse for a Non-Constant Force...
More information11.D.2. Collision Operators
11.D.. Collision Operators (11.94) can be written as + p t m h+ r +p h+ p = C + h + (11.96) where C + is the Boltzmann collision operator defined by [see (11.86a) for convention of notations] C + g(p)
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 informationTangent Planes, Linear Approximations and Differentiability
Jim Lambers MAT 80 Spring Semester 009-10 Lecture 5 Notes These notes correspond to Section 114 in Stewart and Section 3 in Marsden and Tromba Tangent Planes, Linear Approximations and Differentiability
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More informationContact Geometry of Mesoscopic Thermodynamics and Dynamics
Entropy 2014, 16, 1652-1686; doi:10.3390/e16031652 Article OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Contact Geometry of Mesoscopic Thermodynamics and Dynamics Miroslav Grmela École
More informationThe Physics of Fluids and Plasmas
The Physics of Fluids and Plasmas An Introduction for Astrophysicists ARNAB RAI CHOUDHURI CAMBRIDGE UNIVERSITY PRESS Preface Acknowledgements xiii xvii Introduction 1 1. 3 1.1 Fluids and plasmas in the
More informationA Deeper Look into Phase Space: The Liouville and Boltzmann Equations
A Deeper Look into Phase Space: The Liouville and Boltzmann Equations Brian Petko Statistical Thermodynamics Colorado School of Mines December 7, 2009 INTRODUCTION Phase Space, in a sense, can be a useful
More informationShear instabilities. Chapter Energetics of shear instabilities
Chapter 7 Shear instabilities In this final Chapter, we continue our study of the stability of fluid flows by looking at another very common source of instability, shear. By definition, shear occurs whenever
More informationChapter 10 Momentum, System of Particles, and Conservation of Momentum
Chapter 10 Momentum, System of Particles, and Conservation of Momentum 10.1 Introduction... 1 10. Momentum (Quantity of Motion) and Impulse... 1 10..1 Average Force, Momentum, and Impulse... 10.. Non-Constant
More informationSTATISTICAL PHYSICS II
. STATISTICAL PHYSICS II (Statistical Physics of Non-Equilibrium Systems) Lecture notes 2009 V.S. Shumeiko 1 . Part I CLASSICAL KINETIC THEORY 2 1 Lecture Outline of equilibrium statistical physics 1.1
More informationN-body simulations. Phys 750 Lecture 10
N-body simulations Phys 750 Lecture 10 N-body simulations Forward integration in time of strongly coupled classical particles or (rigid-body composite objects) dissipation Brownian motion Very general
More informationPlasmas as fluids. S.M.Lea. January 2007
Plasmas as fluids S.M.Lea January 2007 So far we have considered a plasma as a set of non intereacting particles, each following its own path in the electric and magnetic fields. Now we want to consider
More informationLecture 7 - Separable Equations
Lecture 7 - Separable Equations Separable equations is a very special type of differential equations where you can separate the terms involving only y on one side of the equation and terms involving only
More informationThe Accelerator Hamiltonian in a Straight Coordinate System
Hamiltoninan Dynamics for Particle Accelerators, Lecture 2 The Accelerator Hamiltonian in a Straight Coordinate System Andy Wolski University of Liverpool, and the Cockcroft Institute, Daresbury, UK. Given
More informationUnder evolution for a small time δt the area A(t) = q p evolves into an area
Physics 106a, Caltech 6 November, 2018 Lecture 11: Hamiltonian Mechanics II Towards statistical mechanics Phase space volumes are conserved by Hamiltonian dynamics We can use many nearby initial conditions
More informationMITOCW watch?v=ybcsmyk5xmg
MITOCW watch?v=ybcsmyk5xmg The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources for free. To
More informationStructure formation. Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München
Structure formation Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München Structure formation... Random density fluctuations, grow via gravitational instability galaxies, clusters, etc. Initial perturbations
More informationFrom the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1)
Lecture 9 9. Administration None. 9. Continuation of energy equation From the last time, we ended with an expression for the energy equation ρ D (e + ) u = ρg u + (τ u) q (9.) Where ρg u changes in potential
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamics in High Energy Accelerators Part 2: Basic tools and concepts Nonlinear Single-Particle Dynamics in High Energy Accelerators This course consists of eight lectures: 1.
More informationSpecial Theory of Relativity Prof. Shiva Prasad Department of Physics Indian Institute of Technology, Bombay. Lecture - 15 Momentum Energy Four Vector
Special Theory of Relativity Prof. Shiva Prasad Department of Physics Indian Institute of Technology, Bombay Lecture - 15 Momentum Energy Four Vector We had started discussing the concept of four vectors.
More informationQFT. Unit 11: Cross Sections and Decay Rates
QFT Unit 11: Cross Sections and Decay Rates Decays and Collisions n When it comes to elementary particles, there are only two things that ever really happen: One particle decays into stuff Two particles
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 information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationLecture 5. Alexey Boyarsky. October 21, Legendre transformation and the Hamilton equations of motion
Lecture 5 Alexey Boyarsky October 1, 015 1 The Hamilton equations of motion 1.1 Legendre transformation and the Hamilton equations of motion First-order equations of motion. In the Lagrangian formulation,
More informationPHYSICS 219 Homework 2 Due in class, Wednesday May 3. Makeup lectures on Friday May 12 and 19, usual time. Location will be ISB 231 or 235.
PHYSICS 219 Homework 2 Due in class, Wednesday May 3 Note: Makeup lectures on Friday May 12 and 19, usual time. Location will be ISB 231 or 235. No lecture: May 8 (I m away at a meeting) and May 29 (holiday).
More informationHydrodynamic Limits for the Boltzmann Equation
Hydrodynamic Limits for the Boltzmann Equation François Golse Université Paris 7 & Laboratoire J.-L. Lions golse@math.jussieu.fr Academia Sinica, Taipei, December 2004 LECTURE 2 FORMAL INCOMPRESSIBLE HYDRODYNAMIC
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationParticle-Simulation Methods for Fluid Dynamics
Particle-Simulation Methods for Fluid Dynamics X. Y. Hu and Marco Ellero E-mail: Xiangyu.Hu and Marco.Ellero at mw.tum.de, WS 2012/2013: Lectures for Mechanical Engineering Institute of Aerodynamics Technical
More informationmeans is a subset of. So we say A B for sets A and B if x A we have x B holds. BY CONTRAST, a S means that a is a member of S.
1 Notation For those unfamiliar, we have := means equal by definition, N := {0, 1,... } or {1, 2,... } depending on context. (i.e. N is the set or collection of counting numbers.) In addition, means for
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 informationFrom Newton s law to the linear Boltzmann equation without cut-off
From Newton s law to the linear Boltzmann equation without cut-off Nathalie Ayi (1) (1) Université Pierre et Marie Curie 27 Octobre 217 Nathalie AYI Séminaire LJLL 27 Octobre 217 1 / 35 Organisation of
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 informationFokker-Planck collision operator
DRAFT 1 Fokker-Planck collision operator Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3NP, UK (This version is of 16 April 18) 1. Introduction In these
More informationCourse Okt 28, Statistische Mechanik plus. Course Hartmut Ruhl, LMU, Munich. Personen. Hinweis. Rationale
Okt 28, 2016 ASC, room A239, phone 089-21804210, email hartmut.ruhl@lmu.de. Patrick Böhl, ASC, room A205, phone 089-21804640, email patrick.boehl@physik.uni-muenchen.de. Das folgende Material ist nicht
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationKinetic theory of gases
Kinetic theory of gases Toan T. Nguyen Penn State University http://toannguyen.org http://blog.toannguyen.org Graduate Student seminar, PSU Jan 19th, 2017 Fall 2017, I teach a graduate topics course: same
More informationA = , A 32 = n ( 1) i +j a i j det(a i j). (1) j=1
Lecture Notes: Determinant of a Square Matrix Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk 1 Determinant Definition Let A [a ij ] be an
More informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationLecture 7. 1 Wavepackets and Uncertainty 1. 2 Wavepacket Shape Changes 4. 3 Time evolution of a free wave packet 6. 1 Φ(k)e ikx dk. (1.
Lecture 7 B. Zwiebach February 8, 06 Contents Wavepackets and Uncertainty Wavepacket Shape Changes 4 3 Time evolution of a free wave packet 6 Wavepackets and Uncertainty A wavepacket is a superposition
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
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 informationThe kinetic equation (Lecture 11)
The kinetic equation (Lecture 11) January 29, 2016 190/441 Lecture outline In the preceding lectures we focused our attention on a single particle motion. In this lecture, we will introduce formalism for
More informationHamiltonian aspects of fluid dynamics
Hamiltonian aspects of fluid dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS 01/29/08, 01/31/08 Joris Vankerschaver (CDS) Hamiltonian aspects of fluid dynamics 01/29/08, 01/31/08 1 / 34 Outline
More informationIn this section again we shall assume that the matrix A is m m, real and symmetric.
84 3. The QR algorithm without shifts See Chapter 28 of the textbook In this section again we shall assume that the matrix A is m m, real and symmetric. 3.1. Simultaneous Iterations algorithm Suppose we
More informationLecture 6. Four postulates of quantum mechanics. The eigenvalue equation. Momentum and energy operators. Dirac delta function. Expectation values
Lecture 6 Four postulates of quantum mechanics The eigenvalue equation Momentum and energy operators Dirac delta function Expectation values Objectives Learn about eigenvalue equations and operators. Learn
More informationLiouville Equation. q s = H p s
Liouville Equation In this section we will build a bridge from Classical Mechanics to Statistical Physics. The bridge is Liouville equation. We start with the Hamiltonian formalism of the Classical Mechanics,
More informationHamiltonian Solution I
Physics 4 Lecture 5 Hamiltonian Solution I Lecture 5 Physics 4 Classical Mechanics II September 7th 2007 Here we continue with the Hamiltonian formulation of the central body problem we will uncover the
More informationThe Vlasov equation for cold dark matter and gravity
The Vlasov equation for cold dark matter and gravity Alaric Erschfeld Ruprecht-Karls-Universität Heidelberg Master seminar, 25.11.2016 1 / 50 Table of Contents 1 Introduction 2 The Vlasov equation 3 Virial
More informationKinetic Transport Theory
Lecture notes on Kinetic Transport Theory Christian Schmeiser 1 Preface Kinetic transport equations are mathematical descriptions of the dynamics of large particle ensembles in terms of a phase space (i.e.
More informationANALYSIS QUALIFYING EXAM FALL 2016: SOLUTIONS. = lim. F n
ANALYSIS QUALIFYING EXAM FALL 206: SOLUTIONS Problem. Let m be Lebesgue measure on R. For a subset E R and r (0, ), define E r = { x R: dist(x, E) < r}. Let E R be compact. Prove that m(e) = lim m(e /n).
More informationQuantization of a Scalar Field
Quantization of a Scalar Field Required reading: Zwiebach 0.-4,.4 Suggested reading: Your favorite quantum text Any quantum field theory text Quantizing a harmonic oscillator: Let s start by reviewing
More informationTTK Coherent quantum Boltzmann equations from cqpa Matti Herranen a, Kimmo Kainulainen b;c and Pyry Matti Rahkila b;c a Institut f ur Theoretisc
TTK-10-34 Coherent quantum Boltzmann equations from cqpa Matti Herranen a, Kimmo Kainulainen b;c and Pyry Matti Rahkila b;c a Institut f ur Theoretische Teilchenphysik und Kosmologie, RWTH Aachen University,
More informationIntroductory Analysis I Fall 2014 Homework #9 Due: Wednesday, November 19
Introductory Analysis I Fall 204 Homework #9 Due: Wednesday, November 9 Here is an easy one, to serve as warmup Assume M is a compact metric space and N is a metric space Assume that f n : M N for each
More informationSecond quantization: where quantization and particles come from?
110 Phys460.nb 7 Second quantization: where quantization and particles come from? 7.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system? 7.1.1.Lagrangian Lagrangian
More information4 Divergence theorem and its consequences
Tel Aviv University, 205/6 Analysis-IV 65 4 Divergence theorem and its consequences 4a Divergence and flux................. 65 4b Piecewise smooth case............... 67 4c Divergence of gradient: Laplacian........
More informationDerivation of Kinetic Equations
CHAPTER 2 Derivation of Kinetic Equations As we said, the mathematical object that we consider in Kinetic Theory is the distribution function 0 apple f(t, x, v). We will now be a bit more precise about
More informationPoint vortex mean field theories multi-intensities and kinetic model
Point vortex mean field theories multi-intensities and kinetic model Takashi Suzuki A. Static Theory ~ what should it be? Hamiltonian statistical mechanics (Gibbs) micro-canonical statistics H total energy
More informationThe Distribution Function
The Distribution Function As we have seen before the distribution function (or phase-space density) f( x, v, t) d 3 x d 3 v gives a full description of the state of any collisionless system. Here f( x,
More informationPhysics 212: Statistical mechanics II Lecture XI
Physics 212: Statistical mechanics II Lecture XI The main result of the last lecture was a calculation of the averaged magnetization in mean-field theory in Fourier space when the spin at the origin is
More information[ ] ( ) L = L = r p. We can deal with this more conveniently by writing: so that L becomes: α α α. α α N N N N
Lecture 5: Angular Momentum and Energy of a System The total angular momentum is given by the sum of the angular momenta of all particles in the system: We can deal with this more conveniently by writing:
More informationFluid Mechanics for Astrophysicists
Fluid Mechanics for Astrophysicists MAUCA Fundamental Course 2018 Course Notes Instructor: Oliver HAHN oliver.hahn@oca.eu Contents 1 From Microscopic Particles to the Fluid Equations 3 1.1 Describing many
More informationModule 17: Systems, Conservation of Momentum and Center of Mass
Module 17: Systems, Conservation of Momentum and Center of Mass 17.1 External and Internal Forces and the Change in Momentum of a System So far we have restricted ourselves to considering how the momentum
More informationPhysics 342 Lecture 26. Angular Momentum. Lecture 26. Physics 342 Quantum Mechanics I
Physics 342 Lecture 26 Angular Momentum Lecture 26 Physics 342 Quantum Mechanics I Friday, April 2nd, 2010 We know how to obtain the energy of Hydrogen using the Hamiltonian operator but given a particular
More informationLecture 7: More Arithmetic and Fun With Primes
IAS/PCMI Summer Session 2000 Clay Mathematics Undergraduate Program Advanced Course on Computational Complexity Lecture 7: More Arithmetic and Fun With Primes David Mix Barrington and Alexis Maciel July
More informationPhysics 6303 Lecture 2 August 22, 2018
Physics 6303 Lecture 2 August 22, 2018 LAST TIME: Coordinate system construction, covariant and contravariant vector components, basics vector review, gradient, divergence, curl, and Laplacian operators
More informationEntropy and irreversibility in gas dynamics. Joint work with T. Bodineau, I. Gallagher and S. Simonella
Entropy and irreversibility in gas dynamics Joint work with T. Bodineau, I. Gallagher and S. Simonella Kinetic description for a gas of hard spheres Hard sphere dynamics The system evolves under the combined
More informationBeyond the Second Law of Thermodynamics
Beyond the Second Law of Thermodynamics C. Van den Broeck R. Kawai J. M. R. Parrondo Colloquium at University of Alabama, September 9, 2007 The Second Law of Thermodynamics There exists no thermodynamic
More informationPDE Methods for Mean Field Games with Non-Separable Hamiltonian: Data in Sobolev Spaces (Continued) David Ambrose June 29, 2018
PDE Methods for Mean Field Games with Non-Separable Hamiltonian: Data in Sobolev Spaces Continued David Ambrose June 29, 218 Steps of the energy method Introduce an approximate problem. Prove existence
More informationA dual form of the sharp Nash inequality and its weighted generalization
A dual form of the sharp Nash inequality and its weighted generalization Elliott Lieb Princeton University Joint work with Eric Carlen, arxiv: 1704.08720 Kato conference, University of Tokyo September
More informationSUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, 2011.
SUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, 2011. C. Ringhofer ringhofer@asu.edu, math.la.asu.edu/ chris OVERVIEW Quantum and classical description of many particle
More informationElastic Scattering. R = m 1r 1 + m 2 r 2 m 1 + m 2. is the center of mass which is known to move with a constant velocity (see previous lectures):
Elastic Scattering In this section we will consider a problem of scattering of two particles obeying Newtonian mechanics. The problem of scattering can be viewed as a truncated version of dynamic problem
More informationClassical Statistical Mechanics: Part 1
Classical Statistical Mechanics: Part 1 January 16, 2013 Classical Mechanics 1-Dimensional system with 1 particle of mass m Newton s equations of motion for position x(t) and momentum p(t): ẋ(t) dx p =
More informationComputing the Universe, Dynamics IV: Liouville, Boltzmann, Jeans [DRAFT]
Computing the Universe, Dynamics IV: Liouville, Boltzmann, Jeans [DRAFT] Martin Weinberg June 24, 1998 1 Liouville s Theorem Liouville s theorem is a fairly straightforward consequence of Hamiltonian dynamics.
More informationSelected Topics in Mathematical Physics Prof. Balakrishnan Department of Physics Indian Institute of Technology, Madras
Selected Topics in Mathematical Physics Prof. Balakrishnan Department of Physics Indian Institute of Technology, Madras Module - 11 Lecture - 29 Green Function for (Del Squared plus K Squared): Nonrelativistic
More informationLecture 12: Hydrodynamics in heavy ion collisions. Elliptic flow Last lecture we learned:
Lecture 12: Hydrodynamics in heavy ion collisions. Elliptic flow Last lecture we learned: Particle spectral shapes in thermal model ( static medium) are exponential in m T with common slope for all particles.
More informationL(q, q) = m 2 q2 V (q) 2 m + V (q)
Lecture 7 Phase Space, Part 1 MATH-GA 71.1 Mechanics 1 Phase portraits 1.1 One dimensional system Consider the generic one dimensional case of a point mass m described by a generalized coordinate q and
More informationChapter Five Schneider s Solution to Hilbert s Seventh Problem (and Beyond) (Revised January 2, 2011)
Chapter Five Schneider s Solution to Hilbert s Seventh Problem (and Beyond) (Revised January 2, 2011) In this lecture we will briefly examine Schneider s solution to Hilbert s Seventh Problem and then
More informationAdvanced Topics in Equilibrium Statistical Mechanics
Advanced Topics in Equilibrium Statistical Mechanics Glenn Fredrickson 2. Classical Fluids A. Coarse-graining and the classical limit For concreteness, let s now focus on a fluid phase of a simple monatomic
More informationLecture: Wave-induced Momentum Fluxes: Radiation Stresses
Chapter 4 Lecture: Wave-induced Momentum Fluxes: Radiation Stresses Here we derive the wave-induced depth-integrated momentum fluxes, otherwise known as the radiation stress tensor S. These are the 2nd-order
More informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 1 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
More informationQuantum Mechanics - I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras
Quantum Mechanics - I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 14 Exercises on Quantum Expectation Values (Refer Slide Time: 00:07) In the last couple
More informationBeta and gamma decays
Beta and gamma decays April 9, 2002 Simple Fermi theory of beta decay ² Beta decay is one of the most easily found kinds of radioactivity. As we have seen, this result of the weak interaction leads to
More informationBasic Quantum Mechanics Prof. Ajoy Ghatak Department of Physics Indian Institute of Technology, Delhi
Basic Quantum Mechanics Prof. Ajoy Ghatak Department of Physics Indian Institute of Technology, Delhi Module No. # 05 The Angular Momentum I Lecture No. # 02 The Angular Momentum Problem (Contd.) In the
More informationLecture 10. Central potential
Lecture 10 Central potential 89 90 LECTURE 10. CENTRAL POTENTIAL 10.1 Introduction We are now ready to study a generic class of three-dimensional physical systems. They are the systems that have a central
More informationECE521 Lectures 9 Fully Connected Neural Networks
ECE521 Lectures 9 Fully Connected Neural Networks Outline Multi-class classification Learning multi-layer neural networks 2 Measuring distance in probability space We learnt that the squared L2 distance
More informationPotential/density pairs and Gauss s law
Potential/density pairs and Gauss s law We showed last time that the motion of a particle in a cluster will evolve gradually, on the relaxation time scale. This time, however, is much longer than the typical
More informationMITOCW watch?v=wr88_vzfcx4
MITOCW watch?v=wr88_vzfcx4 PROFESSOR: So we're building this story. We had the photoelectric effect. But at this moment, Einstein, in the same year that he was talking about general relativity, he came
More information