Lecture 18.Pdf
Total Page:16
File Type:pdf, Size:1020Kb
Lecture 18: Ensembles, free energies, max/min principles: A Summary 18.1 Reading Assignment for Lectures 15--18: PKT Chapter 7 Midterm coming up Friday, Oct 29. Problem Set 3 due Oct. 22 (accepted up thru Oct. 25) now available at the website. Summary: (schema) Ensemble/distribution: What you know: Description: Associated free energy: Derivatives: Max/Min principles: 1. Microcanonical ensemble, “equal a priori probability”: Note: From my point of view, this ensemble is the fundamental postulate of equilibrium statistical mechanics and the basis of the derivation of thermodynamics from stat mech. You know values of extensive variables. e.g., E,V,N. 1 1 All states between E and E+ΔE weighted equally, Pn = = . W(E,V,N,("E)) #(E,V,N)"E S(E,V,N) = kB lnW, kB ln " “entropy” 1 P µ dS = dE + dV " dN T T T Isolated system maximizes entropy. ! ! 2. Canonical ensemble: Appropriate for system in contact with heat bath at temperature T. ! E " n kBT Microstates weighted according to Pn ~ e , E E " n " n 1 kBT kBT i.e., Pn = e , Z(T,V,N) = #e , “partition function” Z n ! F(T,V,N) = "kT ln Z = E "TS “Helmholtz free energy” dF = "SdT " PdV + µdN System minimizes F at fixed T. ! 3. Grand (canonical) ensemble: ! Appropriate for system in contact with heat bath at T and particle reservoir(s) at chemical potential(s) µ. ! 1 " (En "µNn ) kBT Microstates weighted according to Pn ~ e , 1 µN # (En #µNn ) 1 kBT kBT i.e., Pn = e , "(T,V,µ) = $e Z(T,V,N). “Grand partition function” " N ! "(T,V,µ) = #kBT ln$ = E #TS + µN = #PV “Grand free energy” d" = #SdT # PdV # Ndµ System minimizes " at fixed T and µ. ! ! 4. Constant pressure ensemble: Appropriate for system in contact with heat bath at T and subject to constant pressure P. ! ! 1 " (En +PVn ) kBT Microstates weighted according to Pn ~ e , i.e., 18.2 1 PV " (En +PVn ) " 1 kBT kBT Pn = e , Z(T,P,N) = # dVe Z(T,V,N) “constant P partition function” Z V G T,P,N k T ln! E TS PV N “Gibbs free energy” ( ) = " B Z = - + = µ dG = "SdT + VdP + µdN System minimizes G at fixed T and P. ! All these are just different fundamental representations of the thermodynamic information. ! For systems of a few degrees of freedom, they are inequivalent (use appropriate one!). ! For thermodynamically large systems, they are equivalent, provided variables are chosen to correspond. ! I want you to be able to use them, as required. I will not expect you to be able to derive them all from the microcanonical distribution, although I have indicated the steps. Let me tie up a few loose ends: The Grand Canonical distribution: Suppose small system s is in thermal+particle contact with the large bath B. i.e., they are in contact across a “barrier” which particles pass freely but which can, for example, withstand pressure differences. Q: How are states of s distributed? Now, both the total energy and the total number of particles of the small system are no longer fixed. (s) (s) Claim: Probability of each microstate n of s with En , N n is 1 " E s"µBN s B ( ) e kBT P = , “Grand Canonical Ensemble” n B B #s(T ,V,µ ) ! ! 1 "S µ "S where T and µ are the temperature and chemical potential of the bath = B ;# = B and the T "EB T "N B 1 % (s) (s)( µNs ! $ ' En $µNn * kBT& ) kBT normalization factor, "s(T,V,µ) # +e = +e Zs(T,V,N s) is called the n N ! s “grand partition function” of the system and the sum is over all states n of s (any number and energy). Proof: s+B is distributed microcanonically, so, for each particular state s,B, the probability is 1 1 Ps+B (E,N) = = . ! total number of states of (s + B)(E to E + "E,N) #s+B (E,N)"E ! Consider now the particular state n of s with energy Es and number Ns. 18.3 " s s s % total number of states of B$ E - E( ) to E - E( ) + (E,N - N ( ) ' " (s) (s) % # n n n & Pn$ En ,N n ' = # & total number of states of (s + B)(E to E + (E,N) . B" (s) (s) % ) $ E * En ,N * N n '( E # & B" (s) (s) % = + ) $ E * En ,N * N n ' )s+B (E,N)(E # & 1 B$ s s ' S & E#E( ),N#N( )) B$ (s) (s) ' k % n n ( But, " & E # E ,N # N ) = e B , so, expanding, % n n ( ! B B B# (s) (s) & (s) )S (s) )S S % E " E ,N " N n ( = SB (E) " E " N + small... $ n ' n )E n )N E,N E,N ! Second-derivative terms are small in that they go as the inverse of extensive quantities. "SB 1 "SB µ B But, = and = # ; hence, B B B ! "E T "N T 1 % s s ( $ ' E( )$µ N( )* B% (s) (s) ( k T & n B n ) P " # ' E $ E ,N $ N * " e B B . QED n & n n ) ! This probability distribution is called the “grand (canonical) ensemble.” Note the “spread” of energy and number in the system! Comment: When there are several types of particles in contact, then, ! 1 $ (s) (s) (s) ' # & En #µ1,BN1,n #µ2,BN2,n #...) kBTB % ( Pn " e . When the “system” is thermodynamically large, then choice of T and µ fixes very precisely the energy * * E and number N that you find in the grand canonical ensemble at values E (T,V,µ),N (T,V,µ), so ! that averages of (almost) all quantities are the same in microcanonical, canonical, and grand canonical ensembles, provided that T,µ,...are chosen to make E=E*, N=N*, etc., i.e., µ is just a “dial” for adjusting N just as T was a dial for adjusting E. Proof: (not required) sketched below: ! You can go all the way back to S; but, I’ll use the previous change to canonical. µN µN F(T,V,N) # # % k T k T k T Thus, for example, "(T,V,µ) = $e B Z(T,V,N) = $e B e B . N=0 N=0 E S(E,V,N) E # % # % k T k T k T Compare: Z(T,V,N) = $ dE "(E,V,N)e B = $ dE e B e B . ! E0 E0 The summand is a product of an increasing function of N and a decreasing function of N. "F * There is a maximum at 0 = µ ! , which fixes N (T, V ,µ) "N T ,V !Can show that this is a maximum and that it is sharp by looking at second and higher derivatives. Upshot: * Sum is dominated by term with (near) N = N . 18.4 µN*#F(T,V,N*) This dominance is strong enough so that "(T,V,µ) = e kT $ O(N)(1+ small); ! * * hence, "(T,V,µ) # $kBT ln% = F(T,V,N (T,V,µ)) $ µN (T,V,µ) + small i.e., " = #kT ln$ = F # µN = E #TS # µN = #PV (Euler) “Grand (canonical) free energy”! (see Lect. 14.1) with, dF(T,V,N) = "SdT " PdV + µdN , so ! d"(T,V,µ) = dF # d(µN) = (#SdT # PdV + µdN) # (µdN + Ndµ) = #SdT # PdV # Ndµ = d" ! Example: Classical Ideal Gas in Grand Ensemble ! N p2 # & % $ 1 V 1 dp 2mk T mk T ! Recall: Z(T,V,N) = % ( , where = & e B = B . N!% d ( " 2# 2 $ "th ' th $% h 2#h µN , µ / µN N # # k T & ) . k T 1 k T e B V V e B Calculate Grand partition function: " = $e B Z T,V,N = $ ( + = exp. 1 . ! ( ) d d ! N 0 N 0 N! ( % + . % 1 = = ' th * . th 1 - 0 µ kTVe kBT Evaluate Grand free energy: " T,V,µ = #k T ln$ = # . ( ) B d ! %th µ µ #$ k T ekBT #$ Ve kBT Calculate derivatives P = " = B ; N = " = . #V d #µ d ! %th %th Eliminating the chemical potential gives the ideal-gas equation of state, as usual. For future reference: µ ! N e kBT ) " N % , A. Note that n , which is equivalent to k T ln dln as above. = = d µ = B + $ ' + (th. V "th * # V & - µ! k T N! e B B. With several different types of (non-interacting classical) particles α, n! " = d . ! ! V #! r 2 r 2 r p p 0 C. If particle has additional internal energy, then "# (p) = $ + "# and 2m# 2m# µ" & 0 N e kBT % " n " e kT . " # = d V $" ! ! .