in order to insure that the Liouville equation for f(?; t) is still valid. These equations of motion will give rise to a distribution function f(?; t)
|
|
- Janice Stone
- 5 years ago
- Views:
Transcription
1 G : Statistical Mechanics Notes for Lecture 21 Consider Hamilton's equations in the form I. CLASSICAL LINEAR RESPONSE THEORY _q i _p i We noted early in the course that an ensemble of systems evolving according to these equations of motion would generate an equilibrium ensemble (in this case, microcanonical). Recall that the phase space distribution function f(?; t) satised a + ilf = where il = f:::; Hg. We noted that =, then f = f(h) is a pure function of the Hamiltonian which dened the general class of distribution functions valid for equilibrium ensembles. What does it mean, however, 6=? To answer this, consider the problem of a simple harmonic oscillator. In an equilibrium ensemble of simple harmonic oscillators at temperature T, the members of the ensemble will undergo oscillatory motion about the potential minimum, with the amplitude of this motion determined by the temperature. Now, however, consider driving each oscillator with a time-dependent driving force F (t). Depending on how complicated the forcing function F (t) is, the motion of each member of the ensemble will, no longer, be simple oscillatory motion about the potential minimum, but could be a very complex kind of motion that explores large regions of the potential energy surface. In other words, the ensemble of harmonic oscillators has been driven away from equilibrium by the time-dependent force F (t). Because of this nonequilibrium behavior of the ensemble, averages over the ensemble could become time-dependent quantities rather than static quantities. Indeed, the distribution function f(?; t), itself, could be time-dependent. This can most easily be seen by considering the equation of motion for a forced oscillator mx =?m! 2 x + F (t) The solution now depends on the entire history of the forcing function F (t), which can introduce explicit timedependence into the ensemble distribution function. A. Generalized equations of motion The most general way a system can be driven away from equilibrium by a forcing function F e (t) is according to the equations of motion: _q i + C i (?)F e (t) _p i + D i (?)F e (t) where the 3N functions C i and D i are required to satisfy the incompressibility condition i i 1
2 in order to insure that the Liouville equation for f(?; t) is still valid. These equations of motion will give rise to a distribution function f(?; + ilf = 6=. (We assume that f is normalized so that R d?f(?; t) = 1.) What does the Liouville equation say about the nature of f(?; t) in the limit that C i and D i are small, so that the displacement away from equilibrius, itself, small? To examine this question, we propose to solve the Liouville equation perturbatively. Thus, let us assume a solution of the form f(?; t) = f (H(?)) + f(?; t) Note, also, that the equations of motion _? take a perturbative form and as a result, the Liouville operator contains two pieces: where il = f:::; Hg and f (H) is assumed to satisfy _?(t) = _? + _?(t) il = _? r? = _? r? + _? r? = il + il il f (H(?)) = _? means the Hamiltonian part of the equations of motion _q i _p i For an observable A(?), the ensemble average of A is a time-dependent quantity: ha(t)i = d?a(?)f(?; t) which, when the assumed form for f(?; t) is substituted in, gives ha(t)i = d?a(?)f (?) + d?a(?)f(?; t) = hai + d?a(?)f(?; t) where hi means average with respect to f (?). B. Perturbative solution of the Liouville equation Substituting the perturbative form for f(?; t) into the Liouville equation, (f (?) + f(?; t)) + (il + il(t))(f (?) + f(?; t)) = =@t =. Thus, working to linear order in small quantities, one obtains the following equation for + il f(?; t) =?ilf (?) which is just a rst-order inhomogeneous dierential equation. This can easily be solved using an integrating factor, and one obtains the result f(?; t) =? dse?il(t?s) il(s)f (?) 2
3 Note that But, using the chain rule, we have ilf (?) = ilf (?)? il f (?) = ilf (?) = _? r? f (?) _? r? f (?) = = = i + _q pi where F i (q 1 ; :::; q N ) is the force on the ith particle. Dene called the dissipative ux. Then j(?) i (F i + D i F e (t))? F i D i (?) p i? C i (?) F i F e (t) Now, suppose f (?) is a canonical distribution function then C i (?) F i? D i (?) p i _? r? f j(?)f e(t) f (H(?)) = 1 Q(N; V; T ) e?h(?) pi + C i F e (t) =?f (H) _? r? f (?) = f (?)j(?)f e (t) Thus, the solution for f(?; t) is f(?; t) =? The ensemble average of the observable A(?) now becomes ha(t)i = hai? = hai? = hai? d?a(?) ds ds dse?il(t?s) f (?)j(?)f e (s) dse?il(t?s) f (?)j(?)f e (s) d?a(?)e?il(t?s) f (?)j(?)f e (s) d?f (?)A(?)e?iL(t?s) j(?)f e (s) Recall that the classical propagator is exp(ilt). Thus the operator appearing in the above expression is a classical propagator of the unperturbed system for propagating backwards in time to?(t? s). An observable A(?) evolves in time according to 3
4 da dt = ila A(t) = e ilt A() A(?t) = e?ilt A() Now, if we take the complex conjugate of both sides, we nd A (t) = A ()e?ilt where now the operator acts to the left on A (). However, since observables are real, we have A(t) = A()e?iLt which implies that forward evolution in time can be achieved by acting to the left on an observable with the time reversed classical propagator. Thus, the ensemble average of A becomes ha(t)i = hai? = hai? dsf e (s) d? f (? )A(? t?s (? ))j(? ) dsf e (s)hj()a(t? s)i where the quantity on the last line is an object we have not encountered yet before. It is known as an equilibrium time correlation function. An equilibrium time correlation function is an ensemble average over the unperturbed (canonical) ensemble of the product of the dissipative ux at t = with an observable A evolved to a time t? s. Several things are worth noting: 1. The nonequilibrium average ha(t)i, in the linear response regime, can be expressed solely in terms of equilibrium averages. 2. The propagator used to evolve A(?) to A(?; t? s) is the operator exp(il (t? s)), which is the propagator for the unperturbed, Hamiltonian dynamics with C i = D i =. That is, it is just the dynamics determined by H. 3. Since A(?; t? s) = A(?(t? s)) is a function of the phase space R variables evolved to a time t? s, we must now specify over which set of phase space variables the integration d? is taken. The choice is actually arbitrary, and for convenience, we choose the initial conditions. Since?(t) is a function of the initial conditions?(), we can write the time correlation function as hj()a(t? s)i = 1 Q d? e?h(?) j(? )A(? t?s (? )) C. General properties of time correlation functions Dene a time correlation function between two quantities A(?) and B(?) by C AB (t) = ha()b(t)i = d?f(?)a(?)e ilt B(?) The following properties follow immediately from the above denition: ha()b(t)i = ha(?t)b()i 1. C AB () = ha(?)b(?)i 2. Thus, if A = B, then 4
5 known as the autocorrelation function of A, and If we dene A = A? hai, then which just measures the uctuations in the quantity A. C AA (t) = ha()a(t)i C AA () = ha 2 i C AA () = h(a) 2 i = h(a? hai) 2 i = ha 2 i? hai 2 3. A time correlation function may be evaluated as a time average, assuming the systes ergodic. In this case, the phase space average may be equated to a time average, and we have C AB (t) = lim T!1 1 T? t T?t dsa(?(s))b(?(t + s)) which is valid for t << T. In molecular dynamics simulations, where the phase space trajectory is determined at discrete time steps, the integral is expressed as a sum X C AB (kt) = 1 N?k A(? k )B(? k+j ) N? k where N is the total number of time steps, t is the time step and N c << N. j=1 k = ; 1; 2; :::; N c 4. Onsager regression hypothesis: In the long time limit, A and B eventually become uncorrelated from each other so that the time correlation function becomes For the autocorrelation function of A, this becomes C AB (t) = ha()b(t)i! haihbi C AA (t)! hai 2 Thus, C AA (t) decays from ha 2 i at t = to hai 2 as t! 1. An example of a signal and its time correlation function appears in the gure below. In this case, the signal is the magnitude of the velocity along the bond of a diatomic molecule interacting with a Lennard-Jones bath. Its time correlation function is shown beneath the signal: 5
6 FIG. 1. Over time, it can be seen that the property being autocorrelated eventually becomes uncorrelated with itself. 6
7 II. TIME CORRELATION FUNCTIONS AND TRANSPORT COEFFICIENTS A. The shear viscosity The shear viscosity of a system measures is resistance to ow. A simple ow eld can be established in a system by placing it between two plates and then pulling the plates apart in opposite directions. Such a force is called a shear force, and the rate at which the plates are pulled apart is the shear rate. A set of microscopic equations of motion for generating shear ow is _q i = p i + y i^x _p i = F i? p yi ^x where is a parameter known as the shear rate. These equations have the conserved quantity H = (p i + y i^x) 2 + U(q 1 ; ::; q N ) The physical picture of this dynamical system corresponds to the presence of a velocity ow eld v(y) = y^x shown in the gure. The ow eld points in the ^x direction and increases with increasing y-value. Thus, layers of a uid, for example, will slow past each other, creating an anisotropy in the system. From the conserved quantity, one can see that the momentum of a particle is the value of p i plus the contribution from the eld evaluated at the position of the particle p i! p i + v(y i ) 7
8 FIG. 2. Such an applied external shearing force will create an asymmetry in the internal pressure. In order to describe this asymmetry, we need an analog of the internal pressure that contains a dependence on specic spatial directions. Such a quantity is known as the pressure tensor and can be dened analogously to the isotropic pressure P that we encountered earlier in the course. Recall that an estimator for the pressure was p = 1 3V p 2 i + q i F i and P = hpi in equilibrium. Here, V is the volume of the system. By analogy, one can write down an estimator for the pressure tensor p : and p = 1 V (pi ^e )(p i ^e ) + (q i ^e )(F i ^e ) P = hp i where ^e is a unit vector in the direction, = x; y; z. This (nine-component) pressure tensor gives information about spatial anisotropies in the system that give rise to o-diagonal pressure tensor components. The isotropic pressure can be recovered from 8
9 P = 1 3 which is just 1/3 of the trace of the pressure tensor. While most systems have diagonal pressure tensors due to spatial isotropy, the application of a shear force according to the above scheme gives rise to a nonzero value for the xy component of the pressure tensor P xy. In fact, P xy is related to the velocity ow eld by a relation of the form X P P xy =? where the coecient is known as the shear viscosity and is an example of a transport coecient. Solving for we nd =? P xy =? lim hp xy (t)i t!1 where hp xy (t)i is the nonequilibrium average of the pressure tensor estimator using the above dynamical equations of motion. Let us apply the linear response formula to the calculation of the nonequilibrium average of the xy component of the pressure tensor. We make the following identications: Thus, the dissipative ux j(?) becomes F e (t) = 1 C i (?) = y i^x D i (?) =?p yi ^x j(?) = According to the linear response formula, so that the shear viscosity becomes = = = V p xy C i F i? D i p i p i ^x y i (F i ^x) + p yi (pi ^y)(p i ^x) hp xy (t)i = hp xy i? V = lim? hp xyi + V t!1 + (q i ^y)(f i ^x) dshp xy ()p xy (t? s)i dshp xy ()p xy (t)i Recall that h i means average of a canonical distribution with =. It is straightforward to show that hp xy i = for an equilibrium canonical distribution function. Finally, taking the limit that t! 1 in the above expression gives the result = V kt 1 dthp xy ()p xy (t)i which is a relation between a transport coecient, in this case, the shear viscosity coecient, and the integral of an equilibrium time correlation function. Relations of this type are known as Green-Kubo relations. Thus, we have expressed a new kind of thermodynamic quantity to an equilibrium time correlation function, which, in this case, is an autocorrelation function of the xy component of the pressure tensor. 9
10 B. The diusion constant The diusive ow of particles can be studied by applying a constant force f to a system using the microscopic equations of motion _q i = p i _p i = F i (q 1 ; ::; q N ) + f ^x which have the conserved energy H = p 2 i 2 + U(q 1 ; :::; q N )? f x i Since the force is applied in the ^x direction, there will be a net ow of particles in this direction, i.e., a current J x. Since this current is a thermodynamic quantity, there is an estimator for it: u x =? and J x = hu x i. The constant force can be considered as arising from a potential eld (x) =?xf The potential will give rise to a concentration which is opposite to the potential gradient and related to it by =? However, Fick's law tells how to relate the particle current J x to the concentration gradient where D is the diusion constant. Solving for D gives J = =? D kt f D =?kt J x f hu x (t)i =?kt lim t!1 f Let us apply the linear response formula again to the above nonequilibrium average. Again, we make the identication: F e (t) = 1 D i = f ^x C i = Thus, hu x (t)i = hu x i? = hu x i? f dsfh ds X i;j! _x i () _x i (t? s) h _x i () _x j (t? s)i! i In equilibrium, it can be shown that there are no cross correlations between dierent particles. Consider the initial value of the correlation function. From the virial theorem, we have which vanishes for i 6= j. In general, h _x i _x j i = ij h _x 2 i i h _x i () _x j (t)i = ij h _x i () _x i (t? s)i 1
11 Thus, hu x (t)i = hu x i? f ds _x i () _x i (t? s)i In equilibrium, hu x i = being linear in the velocities (hence momenta). Thus, the diusion constant is given by, when the limit t! 1 is taken, 1 D = h _x i () _x i (t)i However, since no spatial direction is preferred, we could also choose to apply the external force in the y or z directions and average the result over the these three. This would give a diusion constant D = dt h _q i () _q i (t)i The quantity h _q i () _q i (t)i is known as the velocity autocorrelation function, a quantity we will encounter again in other contexts. 11
G : 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 informationwhere (E) is the partition function of the uniform ensemble. Recalling that we have (E) = E (E) (E) i = ij x (E) j E = ij ln (E) E = k ij ~ S E = kt i
G25.265: Statistical Mechanics Notes for Lecture 4 I. THE CLASSICAL VIRIAL THEOREM (MICROCANONICAL DERIVATION) Consider a system with Hamiltonian H(x). Let x i and x j be specic components of the phase
More informationG : Statistical Mechanics Notes for Lecture 3 I. MICROCANONICAL ENSEMBLE: CONDITIONS FOR THERMAL EQUILIBRIUM Consider bringing two systems into
G25.2651: Statistical Mechanics Notes for Lecture 3 I. MICROCANONICAL ENSEMBLE: CONDITIONS FOR THERMAL EQUILIBRIUM Consider bringing two systems into thermal contact. By thermal contact, we mean that the
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 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 informationG : 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 informationB. Physical Observables Physical observables are represented by linear, hermitian operators that act on the vectors of the Hilbert space. If A is such
G25.2651: Statistical Mechanics Notes for Lecture 12 I. THE FUNDAMENTAL POSTULATES OF QUANTUM MECHANICS The fundamental postulates of quantum mechanics concern the following questions: 1. How is the physical
More informationJavier Junquera. Statistical mechanics
Javier Junquera Statistical mechanics From the microscopic to the macroscopic level: the realm of statistical mechanics Computer simulations Thermodynamic state Generates information at the microscopic
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 informationBrownian Motion and Langevin Equations
1 Brownian Motion and Langevin Equations 1.1 Langevin Equation and the Fluctuation- Dissipation Theorem The theory of Brownian motion is perhaps the simplest approximate way to treat the dynamics of nonequilibrium
More informationUnstable K=K T=T. Unstable K=K T=T
G252651: Statistical Mechanics Notes for Lecture 28 I FIXED POINTS OF THE RG EQUATIONS IN GREATER THAN ONE DIMENSION In the last lecture, we noted that interactions between block of spins in a spin lattice
More informationThus, the volume element remains the same as required. With this transformation, the amiltonian becomes = p i m i + U(r 1 ; :::; r N ) = and the canon
G5.651: Statistical Mechanics Notes for Lecture 5 From the classical virial theorem I. TEMPERATURE AND PRESSURE ESTIMATORS hx i x j i = kt ij we arrived at the equipartition theorem: * + p i = m i NkT
More informationDerivation of the GENERIC form of nonequilibrium thermodynamics from a statistical optimization principle
Derivation of the GENERIC form of nonequilibrium thermodynamics from a statistical optimization principle Bruce Turkington Univ. of Massachusetts Amherst An optimization principle for deriving nonequilibrium
More informationFrom quantum to classical statistical mechanics. Polyatomic ideal gas.
From quantum to classical statistical mechanics. Polyatomic ideal gas. Peter Košovan peter.kosovan@natur.cuni.cz Dept. of Physical and Macromolecular Chemistry Lecture 5, Statistical Thermodynamics, MC260P105,
More informationViscous Fluids. Amanda Meier. December 14th, 2011
Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass density, velocity and pressure. An Eulerian description of uids focuses on the transport
More information6 VORTICITY DYNAMICS 41
6 VORTICITY DYNAMICS 41 6 VORTICITY DYNAMICS As mentioned in the introduction, turbulence is rotational and characterized by large uctuations in vorticity. In this section we would like to identify some
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechanics Rajdeep Sensarma sensarma@theory.tifr.res.in Quantum Dynamics Lecture #2 Recap of Last Class Schrodinger and Heisenberg Picture Time Evolution operator/ Propagator : Retarded
More informationn i,j+1/2 q i,j * qi+1,j * S i+1/2,j
Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/TERMO-5-97 DATE: December 9,997 TITLE A comparison of complete vs. simplied viscous terms in boundary layer
More informationChapter 5. Sound Waves and Vortices. 5.1 Sound waves
Chapter 5 Sound Waves and Vortices In this chapter we explore a set of characteristic solutions to the uid equations with the goal of familiarizing the reader with typical behaviors in uid dynamics. Sound
More information2 GOVERNING EQUATIONS
2 GOVERNING EQUATIONS 9 2 GOVERNING EQUATIONS For completeness we will take a brief moment to review the governing equations for a turbulent uid. We will present them both in physical space coordinates
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 informationSimple Harmonic Motion
Simple Harmonic Motion (FIZ 101E - Summer 2018) July 29, 2018 Contents 1 Introduction 2 2 The Spring-Mass System 2 3 The Energy in SHM 5 4 The Simple Pendulum 6 5 The Physical Pendulum 8 6 The Damped Oscillations
More informationNanoscale simulation lectures Statistical Mechanics
Nanoscale simulation lectures 2008 Lectures: Thursdays 4 to 6 PM Course contents: - Thermodynamics and statistical mechanics - Structure and scattering - Mean-field approaches - Inhomogeneous systems -
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 informationAbsolute Equilibrium Entropy
To appear in the Journal of Plasma Physics Absolute Equilibrium Entropy By JOHN V. SHEBALIN 1 National Aeronautics and Space Administration, Langley Research Center, Hampton, Virginia 3681, USA The entropy
More informationOnsager theory: overview
Onsager theory: overview Pearu Peterson December 18, 2006 1 Introduction Our aim is to study matter that consists of large number of molecules. A complete mechanical description of such a system is practically
More informationIntroduction. Statistical physics: microscopic foundation of thermodynamics degrees of freedom 2 3 state variables!
Introduction Thermodynamics: phenomenological description of equilibrium bulk properties of matter in terms of only a few state variables and thermodynamical laws. Statistical physics: microscopic foundation
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 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 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 informationDamped Oscillations *
OpenStax-CNX module: m58365 1 Damped Oscillations * OpenStax This wor is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 4.0 By the end of this section, you will be
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 informationEnsemble averaged dynamic modeling. By D. Carati 1,A.Wray 2 AND W. Cabot 3
Center for Turbulence Research Proceedings of the Summer Program 1996 237 Ensemble averaged dynamic modeling By D. Carati 1,A.Wray 2 AND W. Cabot 3 The possibility of using the information from simultaneous
More informationHeat and Mass Transfer over Cooled Horizontal Tubes 333 x-component of the velocity: y2 u = g sin x y : (4) r 2 The y-component of the velocity eld is
Scientia Iranica, Vol., No. 4, pp 332{338 c Sharif University of Technology, October 2004 Vapor Absorption into Liquid Films Flowing over a Column of Cooled Horizontal Tubes B. Farhanieh and F. Babadi
More informationStatistical Mechanics in a Nutshell
Chapter 2 Statistical Mechanics in a Nutshell Adapted from: Understanding Molecular Simulation Daan Frenkel and Berend Smit Academic Press (2001) pp. 9-22 11 2.1 Introduction In this course, we will treat
More informationPhysics 70007, Fall 2009 Answers to Final Exam
Physics 70007, Fall 009 Answers to Final Exam December 17, 009 1. Quantum mechanical pictures a Demonstrate that if the commutation relation [A, B] ic is valid in any of the three Schrodinger, Heisenberg,
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationHolographic entropy production
1 1 School of Physics, University of Chinese Academy of Sciences ( 中国科学院大学物理学院 ) (Based on the joint work [arxiv:1204.2029] with Xiaoning Wu and Hongbao Zhang, which received an honorable mention in the
More informationStatistical mechanics of classical systems
Statistical mechanics of classical systems States and ensembles A microstate of a statistical system is specied by the complete information about the states of all microscopic degrees of freedom of the
More informationCHAPTER 8 ENTROPY GENERATION AND TRANSPORT
CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.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 informationIII. Kinetic Theory of Gases
III. Kinetic Theory of Gases III.A General Definitions Kinetic theory studies the macroscopic properties of large numbers of particles, starting from their (classical) equations of motion. Thermodynamics
More informationLECTURE 15 + C+F. = A 11 x 1x1 +2A 12 x 1x2 + A 22 x 2x2 + B 1 x 1 + B 2 x 2. xi y 2 = ~y 2 (x 1 ;x 2 ) x 2 = ~x 2 (y 1 ;y 2 1
LECTURE 5 Characteristics and the Classication of Second Order Linear PDEs Let us now consider the case of a general second order linear PDE in two variables; (5.) where (5.) 0 P i;j A ij xix j + P i,
More informationHydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010
Hydrodynamics Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 What is Hydrodynamics? Describes the evolution of physical systems (classical or quantum particles, fluids or fields) close to thermal
More informationMathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts 1 and Stephan Matthai 2 3rd Febr
HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts and Stephan Matthai Mathematics Research Report No. MRR 003{96, Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL
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 informationC. Show your answer in part B agrees with your answer in part A in the limit that the constant c 0.
Problem #1 A. A projectile of mass m is shot vertically in the gravitational field. Its initial velocity is v o. Assuming there is no air resistance, how high does m go? B. Now assume the projectile is
More informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationSpin Relaxation and NOEs BCMB/CHEM 8190
Spin Relaxation and NOEs BCMB/CHEM 8190 T 1, T 2 (reminder), NOE T 1 is the time constant for longitudinal relaxation - the process of re-establishing the Boltzmann distribution of the energy level populations
More informationDecoherence and Thermalization of Quantum Spin Systems
Copyright 2011 American Scientific Publishers All rights reserved Printed in the United States of America Journal of Computational and Theoretical Nanoscience Vol. 8, 1 23, 2011 Decoherence and Thermalization
More informationFirst Law of Thermodyamics U = Q + W. We can rewrite this by introducing two physical. Enthalpy, H, is the quantity U + pv, as this
First Law of Thermodyamics U = Q + W where U is the increase in internal energy of the system, Q is the heat supplied to the system and W is the work done on the system. We can rewrite this by introducing
More informationUNIVERSITY OF MARYLAND Department of Physics College Park, Maryland. PHYSICS Ph.D. QUALIFYING EXAMINATION PART II
UNIVERSITY OF MARYLAND Department of Physics College Park, Maryland PHYSICS Ph.D. QUALIFYING EXAMINATION PART II August 23, 208 9:00 a.m. :00 p.m. Do any four problems. Each problem is worth 25 points.
More information4 Evolution of density perturbations
Spring term 2014: Dark Matter lecture 3/9 Torsten Bringmann (torsten.bringmann@fys.uio.no) reading: Weinberg, chapters 5-8 4 Evolution of density perturbations 4.1 Statistical description The cosmological
More information4.1 Eigenvalues, Eigenvectors, and The Characteristic Polynomial
Linear Algebra (part 4): Eigenvalues, Diagonalization, and the Jordan Form (by Evan Dummit, 27, v ) Contents 4 Eigenvalues, Diagonalization, and the Jordan Canonical Form 4 Eigenvalues, Eigenvectors, and
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 notes for /12.586, Modeling Environmental Complexity. D. H. Rothman, MIT October 22, 2014
Lecture notes for 12.086/12.586, Modeling Environmental Complexity D. H. Rothman, MIT October 22, 2014 Contents 1 Origin of biogeochemical cycles 1 1.1 The carbon cycle......................... 1 1.1.1
More informationOn modeling pressure diusion. in non-homogeneous shear ows. By A. O. Demuren, 1 M. M. Rogers, 2 P. Durbin 3 AND S. K. Lele 3
Center for Turbulence Research Proceedings of the Summer Program 1996 63 On modeling pressure diusion in non-homogeneous shear ows By A. O. Demuren, 1 M. M. Rogers, 2 P. Durbin 3 AND S. K. Lele 3 New models
More informationG : Statistical Mechanics
G5.651: Statistical Mechanics Notes for Lecture 1 I. DERIVATION OF THE DISCRETIZED PATH INTEGRAL We begin our discussion of the Feynman path integral with the canonical ensemble. The epressions for the
More informationContents. 6 Systems of First-Order Linear Dierential Equations. 6.1 General Theory of (First-Order) Linear Systems
Dierential Equations (part 3): Systems of First-Order Dierential Equations (by Evan Dummit, 26, v 2) Contents 6 Systems of First-Order Linear Dierential Equations 6 General Theory of (First-Order) Linear
More informationHAMILTON S PRINCIPLE
HAMILTON S PRINCIPLE In our previous derivation of Lagrange s equations we started from the Newtonian vector equations of motion and via D Alembert s Principle changed coordinates to generalised coordinates
More informationEconomics 472. Lecture 10. where we will refer to y t as a m-vector of endogenous variables, x t as a q-vector of exogenous variables,
University of Illinois Fall 998 Department of Economics Roger Koenker Economics 472 Lecture Introduction to Dynamic Simultaneous Equation Models In this lecture we will introduce some simple dynamic simultaneous
More informationA. Bers, A. K. Ram, and S. D. Schultz. Plasma Science and Fusion Center,
COUPLING TO ELECTRON BERNSTEIN WAVES IN TOKAMAKS* A. Bers, A. K. Ram, and S. D. Schultz Plasma Science and Fusion Center, Massachusetts Institute of Technology, Cambridge, MA 02139. U.S.A. Abstract The
More informationBoundary Dissipation in a Driven Hard Disk System
Boundary Dissipation in a Driven Hard Disk System P.L. Garrido () and G. Gallavotti (2) () Institute Carlos I for Theoretical and Computational Physics, and Departamento de lectromagnetismo y Física de
More informationA Chaotic Phenomenon in the Power Swing Equation Umesh G. Vaidya R. N. Banavar y N. M. Singh March 22, 2000 Abstract Existence of chaotic dynamics in
A Chaotic Phenomenon in the Power Swing Equation Umesh G. Vaidya R. N. Banavar y N. M. Singh March, Abstract Existence of chaotic dynamics in the classical swing equations of a power system of three interconnected
More informationMapping Closure Approximation to Conditional Dissipation Rate for Turbulent Scalar Mixing
NASA/CR--1631 ICASE Report No. -48 Mapping Closure Approximation to Conditional Dissipation Rate for Turbulent Scalar Mixing Guowei He and R. Rubinstein ICASE, Hampton, Virginia ICASE NASA Langley Research
More informationPhysics 106b: Lecture 7 25 January, 2018
Physics 106b: Lecture 7 25 January, 2018 Hamiltonian Chaos: Introduction Integrable Systems We start with systems that do not exhibit chaos, but instead have simple periodic motion (like the SHO) with
More informationwhere the viscosity coecients are functions of the invariants of, A = 1 + rv, and are
Appendix: Final form of the algebraic mass ux model The pertinent facts concerning the computational implementation of the mass ux representation are presented here. The general form of the model, valid
More information424 Index. Eigenvalue in quantum mechanics, 174 eigenvector in quantum mechanics, 174 Einstein equation, 334, 342, 393
Index After-effect function, 368, 369 anthropic principle, 232 assumptions nature of, 242 autocorrelation function, 292 average, 18 definition of, 17 ensemble, see ensemble average ideal,23 operational,
More information10 Shallow Water Models
10 Shallow Water Models So far, we have studied the effects due to rotation and stratification in isolation. We then looked at the effects of rotation in a barotropic model, but what about if we add stratification
More informationHandout 10. Applications to Solids
ME346A Introduction to Statistical Mechanics Wei Cai Stanford University Win 2011 Handout 10. Applications to Solids February 23, 2011 Contents 1 Average kinetic and potential energy 2 2 Virial theorem
More informationgr-qc/ Sep 1995
DAMTP R95/45 On the Evolution of Scalar Metric Perturbations in an Inationary Cosmology R. R. Caldwell University of Cambridge, D.A.M.T.P. Silver Street, Cambridge CB3 9EW, U.K. email: R.R.Caldwell@amtp.cam.ac.uk
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More informationAnalysis of the simulation
Analysis of the simulation Marcus Elstner and Tomáš Kubař January 7, 2014 Thermodynamic properties time averages of thermodynamic quantites correspond to ensemble averages (ergodic theorem) some quantities
More information5.80 Small-Molecule Spectroscopy and Dynamics
MIT OpenCourseWare http://ocw.mit.edu 5.80 Small-Molecule Spectroscopy and Dynamics Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 5.80 Lecture
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 informationP 1 P * 1 T P * 1 T 1 T * 1. s 1 P 1
ME 131B Fluid Mechanics Solutions to Week Three Problem Session: Isentropic Flow II (1/26/98) 1. From an energy view point, (a) a nozzle is a device that converts static enthalpy into kinetic energy. (b)
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 informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationto satisfy the large number approximations, W W sys can be small.
Chapter 12. The canonical ensemble To discuss systems at constant T, we need to embed them with a diathermal wall in a heat bath. Note that only the system and bath need to be large for W tot and W bath
More information5 Eigenvalues and Diagonalization
Linear Algebra (part 5): Eigenvalues and Diagonalization (by Evan Dummit, 27, v 5) Contents 5 Eigenvalues and Diagonalization 5 Eigenvalues, Eigenvectors, and The Characteristic Polynomial 5 Eigenvalues
More informationTwo recent works on molecular systems out of equilibrium
Two recent works on molecular systems out of equilibrium Frédéric Legoll ENPC and INRIA joint work with M. Dobson, T. Lelièvre, G. Stoltz (ENPC and INRIA), A. Iacobucci and S. Olla (Dauphine). CECAM workshop:
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 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 informationTime Evolution of Modeled Reynolds Stresses in Planar Homogeneous Flows
NASA/TM-97-206265 Time Evolution of Modeled Reynolds Stresses in Planar Homogeneous Flows T. Jongen Swiss Federal Institute of Technology, Lausanne, Switzerland T. B. Gatski Langley Research Center, Hampton,
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 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 informationspring mass equilibrium position +v max
Lecture 20 Oscillations (Chapter 11) Review of Simple Harmonic Motion Parameters Graphical Representation of SHM Review of mass-spring pendulum periods Let s review Simple Harmonic Motion. Recall we used
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 9 I. OVERVIEW Our treatment of the classical ensembles makes clear that the free energy is a quantity of particular importance in statistical mechanics.
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 informationContinuum Mechanics Fundamentals
Continuum Mechanics Fundamentals James R. Rice, notes for ES 220, 12 November 2009; corrections 9 December 2009 Conserved Quantities Let a conseved quantity have amount F per unit volume. Examples are
More information1.11 Some Higher-Order Differential Equations
page 99. Some Higher-Order Differential Equations 99. Some Higher-Order Differential Equations So far we have developed analytical techniques only for solving special types of firstorder differential equations.
More informationChapter 9: Statistical Mechanics
Chapter 9: Statistical Mechanics Chapter 9: Statistical Mechanics...111 9.1 Introduction...111 9.2 Statistical Mechanics...113 9.2.1 The Hamiltonian...113 9.2.2 Phase Space...114 9.2.3 Trajectories and
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 informationSpurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics
UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Spurious Chaotic Solutions of Dierential Equations Sigitas Keras DAMTP 994/NA6 September 994 Department of Applied Mathematics and Theoretical Physics
More informationLecture 8: Ch Sept 2018
University of Alabama Department of Physics and Astronomy PH 301 / LeClair Fall 2018 Lecture 8: Ch. 4.5-7 12 Sept 2018 1 Time-dependent potential energy Sometimes we need to consider time-dependent forces,
More informationPath integral in quantum mechanics based on S-6 Consider nonrelativistic quantum mechanics of one particle in one dimension with the hamiltonian:
Path integral in quantum mechanics based on S-6 Consider nonrelativistic quantum mechanics of one particle in one dimension with the hamiltonian: let s look at one piece first: P and Q obey: Probability
More informationStatistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany
Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals From Micro to Macro Statistical Mechanics (Statistical
More informationOur goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always
October 5 Relevant reading: Section 2.1, 2.2, 2.3 and 2.4 Our goal is to solve a general constant coecient linear second order ODE a d2 y dt + bdy + cy = g (t) 2 dt where a, b, c are constants and a 0.
More information7 The Quasi-Biennial Oscillation (QBO)
7 The Quasi-Biennial Oscillation (QBO) (Reviewed by Baldwin et al., Rev. Geophys., 001) Previously we noted the remarkable quasi-periodic reversal of zonal winds in the tropical stratosphere, the quasi-biennial
More informationContents. Contents. Contents
Physics 121 for Majors Class 18 Linear Harmonic Last Class We saw how motion in a circle is mathematically similar to motion in a straight line. We learned that there is a centripetal acceleration (and
More informationAsymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks
Physica A 275 (2000) 483 504 www.elsevier.com/locate/physa Asymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks V. Kumaran Department of Chemical Engineering, Indian
More information