Canonical Ensemble
Total Page:16
File Type:pdf, Size:1020Kb
ME346A Introduction to Statistical Mechanics { Wei Cai { Stanford University { Win 2011 Handout 8. Canonical Ensemble January 26, 2011 Contents Outline • In this chapter, we will establish the equilibrium statistical distribution for systems maintained at a constant temperature T , through thermal contact with a heat bath. • The resulting distribution is also called Boltzmann's distribution. • The canonical distribution also leads to definition of the partition function and an expression for Helmholtz free energy, analogous to Boltzmann's Entropy formula. • We will study energy fluctuation at constant temperature, and witness another fluctuation- dissipation theorem (FDT) and finally establish the equivalence of micro canonical ensemble and canonical ensemble in the thermodynamic limit. (We first met a mani- festation of FDT in diffusion as Einstein's relation.) Reading Assignment: Reif x6.1-6.7, x6.10 1 1 Temperature For an isolated system, with fixed N { number of particles, V { volume, E { total energy, it is most conveniently described by the microcanonical (NVE) ensemble, which is a uniform distribution between two constant energy surfaces. const E ≤ H(fq g; fp g) ≤ E + ∆E ρ (fq g; fp g) = i i (1) mc i i 0 otherwise Statistical mechanics also provides the expression for entropy S(N; V; E) = kB ln Ω. In thermodynamics, S(N; V; E) can be transformed to a more convenient form (by Legendre transform) of Helmholtz free energy A(N; V; T ), which correspond to a system with constant N; V and temperature T . Q: Does the transformation from N; V; E to N; V; T have a meaning in statistical mechanics? A: The ensemble of systems all at constant N; V; T is called the canonical NVT ensemble. Q: What is the probability distribution function for the canonical ensemble? Is it uniform? A: The distribution is not uniform. It is called Boltzmann's distribution, which we will develop below. Q: How do we define temperature in statistical mechanics? A: In thermodynamics, temperature is established through the zeroth law | transitivity of thermal equilibrium. This means that temperature is a property that emerges in the thermodynamic limit, N ! 1, where N is the number of particles. | The existence of temperature is independent of the type of the material (e.g. gas, liquid, or solid). On the other hand, a system with say, N < 10, particles does not have a well defined temperature (in the usual thermodynamic sense), unless it is in thermal contact with another system, i.e., a thermostat, which has a large number of particles and a well- defined temperature. The thermostat is considered to be many many times larger than the system of inter- est, so that the heat exchange with the system has a negligible effect on the energy per particle and temperature of the thermostat. Through thermal contact with the thermostat, the system will eventually settle down to an equilibrium state with temperature equal to that of the thermostat. 2 The only requirement for the thermostat is that it is much larger than the system of interest. | You could imagine the rest of the universe as the thermostat. The only problem there is that the universe is not in thermal equilibrium. For convenience, the ideal gas is usually used as a model for the thermostat (because its analytic expression for S(N; V; E) is known). But this is not necessary. The property of the system of interest will not change if we use a different type of thermostat. 2 Boltzmann's distribution Because the system can exchange heat with the thermostat, its total energy is no longer conserved. Hence, we should no longer expect its energy to be confined as E ≤ H(fqig; fpig) ≤ E + ∆E Instead, any energy is allowed. Q: How can we obtain the equilibrium density distribution ρc(fqig; fpig) for the canonical ensemble? 3 A: The key idea is that, the (system + thermostat) can be considered an isolated system, and can be described by the microcanonical ensemble. Let H be the Hamiltonian of the system and HB be the Hamiltonian of the heat bath. ^ H + HB = E is conserved. Hence the system and heat bath together can be described by the microcanonical ensemble when they have reached thermal equilibrium. const E^ ≤ H(fq g; fp g) + H (fqBg; fpBg) ≤ E^ + ∆E ρ (fq g; fp g; fqBg; fpBg) = i i B i i mc i i i i 0 otherwise (2) B B where fqig; fpig represent the system's degrees of freedom, and fqi g; fpi g represent the heat bath's degrees of freedom. If we only want to know the statistical distribution concerning the system's degrees of freedom, all we need to do is to integrate out the degrees of freedom corresponding to the heat bath. Z 3NB Y B B B B ρc(fqig; fpig) = dqj dpj ρmc(fqig; fpig; fqi g; fpi g) (3) j=1 This is very similar to the study of momentum distribution of one gas molecule in a gas tank of N molecules (see \Microcanonical ensemble" notes). ) ρ(fqig; fpig) is proportional to the number of ways the thermostat molecules can rearrange ^ B B ^ themselves such that E ≤ H(fqig; fpig) + HB(fqi g; fpi g) ≤ E + ∆E. ^ ρc(fqig; fpig) = const · ΩB(E − H(fqig; fpig)) (4) 4 Recall that the entropy SB of the heat bath is a function of its energy EB. ^ SB(EB) = kB ln ΩB(EB) EB = E − H(fqig; fpig) (5) ^ Also notice that E, EB E = H(fqig; fpig). Hence we can do a Taylor expansion of SB around E = 0. ^ @SB SB(EB) ≈ SB(E) − · H(fqig; fpig) (6) @EB The higher order terms are neglected. We recognize that 1 ≡ @SB , where T is the temperature of the thermostat. T must also T @EB be the temperature of the system at thermal equilibrium. Therefore, H(fq g; fp g) S (E ) = const − i i (7) B B T SB(EB) H(fqig; fpig) ΩB(EB) = exp = const · exp − (8) kB kBT H(fqig; fpig) ρc(fqig; fpig) = const · exp − (9) kBT This is the canonical distribution, also called Boltzmann's distribution or Boltzmann's law. 1 H(fqig; fpig) ρc(fqig; fpig) = exp − (10) Z~ kBT The term in the exponential is called Boltzmann's factor. Z~ is a normalization factor Z 3N Y H(fqig; fpig) Z~ ≡ dq dp exp − (11) i i k T i=1 B We can compare this result with the momentum distribution of one gas molecule in the microcanonical ensemble (see \Microcanonical ensemble" notes). 3 Helmholtz free energy Do not say \it's just a normalization factor". We will see that Z~ is actually very important. In the miocrocanonical ensemble, we also have a nomalization constant Ω~ (volume of phase space) 1 1 E ≤ H(fqig; fpig) ≤ E + ∆E ρmc(fqig; fpig) = (12) Ω~ 0 otherwise 5 Ω~ is eventually connected to entropy Ω~ S(N; V; E) = k ln (13) B N!h3N In other words, Ω~ S = k ln Ω; Ω ≡ (14) B N!h3N Q: What's the meaning of the normalization constant Z in the canonical ensemble? Q: What is entropy in the canonical ensemble? A: Use Shannon's entropy formula, X S = −kB pi ln pi i=1 Z 1 H H ~ = −kB dq dp exp − · − − ln Z Γ Z kBT kBT Z 1 ~ = dq dp ρc(q; p) H(q; p) + kB ln Z (15) T Γ R where Γ means integration over the entire phase space. Notice that the first term simply contains the ensemble average of energy. Z E ≡ hHi ≡ dq dp ρc(q; p) H(q; p) (16) Γ E S = + k ln Z;~ k T ln Z~ = TS − E (17) ) T B B Recall in thermodynamics, the Helmholtz free energy is defined as A ≡ E − TS. Therefore, ~ A = −kBT ln Z (18) Adding quantum corrections Z~ A(N; V; T ) = −k T ln (19) B N! h3N Z~! Z ≡ N!h3N is called the partition function. Notice that energy E is the ensemble average of the Hamiltonian. But the free energy A cannot be written as an ensemble average. A is proportional to the log of the normalization constant A = −kB ln Z. (Analogous to Boltzmann's entropy expression for the canonical ensemble S = kB ln Ω.) 6 4 Fluctuation-dissipation theorem on energy In the canonical ensemble, the system acquire a temperature by having a thermal contact with a thermostat (heat bath) with temperature T . Thus the system is no longer isolated any more. Its total energy, i.e., Hamiltonian H(fqig; fpig) is no longer conserved. In other words, we should expect some fluctuation of total energy in the canonical ensemble. On the other hand, fluctuations are not considered in thermodynamics. At constant N; V; T , the appropriate thermodynamics potential is A(N; V; T ), from which we can compute a definite value for energy E = A + TS, with S ≡ − (@A=@T )N;V . Hence, in thermodynamics, we expect the system to simultaneously have a definite tem- perature T and total energy E. In statistical mechanics, if the system have a well defined temperature, its total energy E must fluctuate. Q: How do we reconcile this difference between statistical mechanics and thermodynamics? A: (1) The total energy in thermodynamics should be identified as the ensemble average of the Hamiltonian in statistical mechanics. 3N 1 Z Y 1 E = hHi ≡ dq dp exp(−βH(fq g; fp g)) · H(fq g; fp g); β ≡ (20) ~ i i i i i i k T Z i=1 B (2) We expect the statistical fluctuation of the Hamiltonian to diminish in the thermody- namic limit (N ! 1).