- Email: [email protected]
THE GIBBS DISTRIBUTION
THE GIBBS DISTRIBUTION § 28. The Gibbs distribution L e t u s now turn to the problem stated in Chapter I of finding the distribution function for a subsystem, i.e. any macroscopic body which is a small part of some large closed system. T h e most convenient a n d general method of approaching the solution of this problem is based on the application of the microcanonical distribution t o the whole system. Distinguishing the body in question from the rest of the closed system, we may consider the system as consisting of these two parts. T h e rest of the system will be called the " m e d i u m " in relation to the body. T h e microcanonical distribution (6.6) can be written in the form 0
dw = c o n s t a n t X o ( £ + £ ' - # > ) d r d r ' ,
(28.1)
where E, d T and E\ dT' relate t o the body a n d the medium respectively, a n d ( 0) 2? is the given value of the energy of the closed system, which must be equal to the sum E+E' of the energies of the body a n d the medium. O u r object is t o find the probability wn of a state of the whole system such that the body concerned is in some definite quantum state (with energy Ε J, i.e. a microscopically defined state. T h e microscopic state of the medium is of n o interest, so that we shall suppose this t o be in a macroscopically defined state. L e t A T " be the statistical weight of the macroscopic state of the medium a n d let ΔΕ' be the range of values of the energy of the medium corresponding t o the range ΔΓ' of quantum states in the sense discussed in §7. T h e required probability wn can be found by taking άΓ = 1 in (28.1), putting Ε = En a n d integrating with respect t o Γ ' : wn =
œnstantxJo^+E'-J^dr'.
Let-T^isO be the total number of quantum states of the medium with energy n o t exceeding E'. Since the integrand depends only on E\ we can change t o integration with respect t o E\ putting dF = (dr'(EydE')dE'. T h e derivative dF/dE' is replaced (cf. § 7) by
drVd£ = r
sp
- d
eS'W/AE', 79
80
The Gibbs
Distribution
where .S"(JE") is the entropy of the medium as a function of its energy; ΔΕ' is, of course, also a function of E\ Thus we have wn = constant Χ | ~
<5(£' +En-Ε™)
àE'.
Owing to the presence of the delta function, the result of the integration is m simply to replace E' by E -En: w„ = constantX (-^=Λ
.
(28.2)
We now use the fact that, since the body is small, its energy En is small in ( 0 ) comparison with £ . The quantity ΔΕ' undergoes only a very small relative (0) change when E' varies slightly, and so in ΔΕ' we can simply p u t E' = 2 s ; it then becomes a constant independent of En. In the exponential factor e ^ , {0) we must expand S'(E —En) in powers of En as far as the linear term : S'(EW-En)
=
sr(EW)-EHâSr(EW)làE<®.
The derivative of the entropy S' with respect t o energy is just 1/Γ, where Τ is the temperature of the system; the temperatures of the body and the mediu m are the same, since the system is assumed to be in equilibrium. Thus we have finally for wn the expression w„ = Ae-W,
(28.3)
where A is a normalisation constant independent of En. This is one of the most important formulae in statistical physics. It gives the statistical distribution of any macroscopic body which is a comparatively small part of a large closed system. The distribution (28.3) is called the Gibbs distribution or canonical distribution*, it was discovered by J. W . Gibbs for classical statistics in 1901. The normalisation constant A is given by the condition Ewn = 1, whence
7 =Σ ^ ·
( ·) 28
4
The mean value of any physical quantity / pertaining to the body can b e calculated by means of the Gibbs distribution, using the formula η
= ΣΛ^-^/Σ*-^·
(28-5)
§28
The Gibbs
81
Distribution
In classical statistics an expression exactly corresponding to (28.3) is obtained for the distribution function in phase space : Q(p,q)=
(28.6)
Ae-****,
where E(p9 q) is the energy of the body as a function of its coordinates and 1. momenta. The normalisation constant A is given by the condition E
T
J ρ àp dq = A j e~ ^> M dp dq = 1.
(28.7)
In practice, cases are frequently encountered where it is n o t the entire microscopic motion of the particles which is quasi-classical, but only the motion corresponding to some of the degrees of freedom, whereas the motion with respect to the remaining degrees of freedom is quantised (for example, the translational motion of the molecules may be quasi-classical while the motion of the atoms within the molecules is quantised). Then the energy levels of the body may be written as functions of the quasi-classical coordinates and momenta : En = En(p9 q\ where η denotes the set of q u a n t u m numbers defining the "quantised p a r t " of the motion, for which ρ and q are parameters. T h e Gibbs distribution formula then becomes dw„(p, q) = Ae-stP-Wdpcl
dqcl,
(28.8)
where dpcldqcl is the product of differentials of the "quasi-classical" coordinates and momenta. Finally, the following comment is necessary concerning the group of p r o b lems which may be solved by means of the Gibbs distribution. We have spoken of the latter throughout as the statistical distribution for a subsystem, as in fact it is. It is very important to note, however, that this same distribution can quite succesfully be used also to determine the fundamental statistical properties of bodies forming closed systems, since such properties of a body as the values of the thermodynamic quantities or the probability distributions for the coordinates and velocities of its individual particles are clearly independent of whether we regard the body as a closed system or as being placed in an imaginary thermostat (§ 7). But in the latter case the body becomes a "subsystem" and the Gibbs distribution is immediately applicable to it. T h e difference between bodies forming closed and non-closed systems when the Gibbs distribution is used appears essentially only in the treatment of t h e fairly unimportant problem of fluctuations in the total energy of the body. The Gibbs distribution gives for the mean fluctuation of this quantity a n o n t To avoid misunderstanding, let us mention once more that the wn (or ρ) are monotonie functions of energy and need not have maxima for E= E. It is the distribution function with respect to energy, obtained by multiplying wn by άΓ(Ε)/άΕ, wich has a sharp maximum at is = E.
82
The Gibbs
Distribution
zero value, which is meaningful for a body in a medium but is entirely spurious for a closed system, since the energy of such a body is by definition constant and does n o t fluctuate. The possibility of applying the G i b b s distribution (in the manner described) t o closed systems is also seen from the fact that this distribution hardly differs from the microcanonical distribution, while being very much more convenient for practical calculations. F o r the microcanonical distribution is, roughly speaking, equivalent to regarding as equally probable all microstates of the body which correspond to a given value of its energy. T h e canonical distribution is " s p r e a d " over a certain range of energy values, b u t the width of this range (of the order of the mean fluctuation of energy) is negligible for a macroscopic body.
§ 29. The Maxwellian distribution T h e energy E(p9 q) in the Gibbs distribution formula of classical statistics can always be written as the sum of two p a r t s : the kinetic energy a n d the potential energy. T h e first of these is a quadratic function of the momenta of 1 the a t o m s , and the second is a function of their coordinates, the form of which depends on the law of interaction between the particles within the body (and o n the external field, if any). If the kinetic a n d potential energies are denoted by K(p) and U(q) respectively, then E(p, q) = K(p) + U(q\ and the probability dw = ρ(ρ, q) dp dq becomes dw =
Ae-vuWe-xWdpdq,
i.e. is the product of two factors, one of which depends only on the coordinates and the other only on the momenta. This means that the probabilities for momenta and coordinates are independent, in the sense that any particular values of the momenta d o n o t influence the probabilities of the various values of the coordinates, and vice versa. Thus the probability of the various values of the momenta can be written K
T
dwp = ae- W dp,
(29.1)
a n d the probability distribution for the coordinates is u T
dwq = be- ^ dq.
(29.2)
Since the sum of the probabilities of all possible values of the momenta must be unity (and the same applies to the coordinates), each of the probabilities dwp and dwq must be normalised, i.e. their integrals over all possible t It is assumed that Cartesian coordinates are used.
§ 29
The Maxwellian
83
Distribution
values of the momenta and coordinates respectively for the body concerned must be equal to unity. F r o m these conditions we can determine the constants a and b in (29.1) and (29.2). Let us consider the probability distribution for the momenta, and once again emphasise the very important fact that in classical statistics this distribution does n o t depend on the nature of the interaction of particles within the system or on the nature of the external field, and so can be expressed in a form applicable to all bodies.* The kinetic energy of the whole body is equal to the sum of the kinetic energies of each of the atoms composing it, and the probability again falls into a product of factors, each depending on the momenta of only one atom. This again shows that the m o m e n t u m probabilities of different atoms are independent, i.e. the momentum of one does not affect the probabilities o f various momenta of any other. We can therefore write the probability distribution for the momenta of each atom separately. F o r an atom of mass m the kinetic energy is (pl+pl+p*)/2m, where px> Py, pz are the Cartesian components of its momentum, and the probability distribution is +p
p
dwp = cw-tâ l+ W
2mT
dpx dpy
dp2.
The constant a is given by the normalisation condition. The integrations over px, ρy and pz are separate; by means of the formula oo
x
J e-« * dx
= V(^/a)
— oo
s,2
we find a = (2nmT)~ the final form *"·
9
=
and the m o m e n t u m probability distribution takes
e
( M ? * ~
(
P
l
+
P
l
^ ' *
m
T
Changing from momenta to velocities (p = ms\ sponding velocity distribution as
/ m\ d Wy
=
U!L.\
à x P à yP
ά
ρ
-
·
3)
( 2 9
we can write the corre-
s/2
-m(vlWi+v\)/2T e
X D VÀ YV D ZV
4)
( 2 9 e
This is the Maxwellian distribution (J. C. Maxwell, 1860). It again consists of a product of three independent factors
t In quantum statistics this statement is not true in general.
84
The Gibbs
Distribution
each of which gives the probability distribution for a single velocity c o m p o nent. If the body consists of molecules (e.g. a polyatomic gas), then together with the Maxwellian distribution for the individual atoms there is a similar distribution for the translational motion of each molecule as a whole: from the kinetic energy of the molecule we can separate a term which gives the energy of the translational motion, and so the required distribution separates in the form (29.4), where m must now be taken as the total mass of the molecule, a n d vx9 Vyy vz as the velocity components of its centre of mass. It should be emphasised that the Maxwellian distribution for the translational motion of molecules can be valid quite independently of the nature of the motion of the atoms within the molecule (and the rotation of the molecule), and in partic1 ular when a quantised description of the latter is necessary. The expression (29.4) is written in terms of Cartesian coordinates in "velocity space". If we change from Cartesian to spherical polar coordinates, the result is \
-mvmT e v2
g 0 φ ά
ns i d
d yV
(29.6)
where ν is the absolute magnitude of the velocity, and 0 and φ the polar angle and azimuthal angle which determine the direction of the velocity. Integration with respect to angle gives the probability distribution for the absolute magnitude of the velocity: dwv =
2^)
3 /2
*-™ ν dv. β/23
(29.7)
It is sometimes convenient to use cylindrical coordinates in velocity space. Then d
Wr =
emiOlv ),2Tv
( ί τ ) ' ~ +' r
dv, dv2 άφ,
(29.8)
where vx is the velocity component along the z-axis, vr the component perpendicular to that axis, and φ the angle which gives the direction of this component. Let us calculate the mean kinetic energy of an atom. According to the definition of the mean, and using (29.5), we find for any Cartesian velocity com-
t The Maxwellian distribution clearly applies also to the Brownian motion of particles suspended in a liquid.
§29
The Maxwellian
85
Distribution
ponenf
— oo
(29.9)
= T/m.
The mean value of the kinetic energy of the atom is therefore 3772. We can thus say that the mean kinetic energy of all the particles in the body in classical statistics is always 3NT/2, where Ν is the total number of particles. PROBLEMS
PROBLEM. 1. Find the mean value of the nth power of the absolute magnitude of the velocity. SOLUTION. Using (29.7), we find
^
=
^J
{^τ)
4π
4
e-^V "
2
dv
ο 2 /27γ/*
/n + 3\
t For reference we shall give the values of the integrals of the form eo e
/„ = J é>- *Vd*, ο which often occur in applications of the Maxwellian distribution. The substitution 2 (xx = y gives oo
ο =
ia-<»+i"T(i«+|),
where Γ(χ) is the gamma function. In particular, if η = 2r with r > 0, then
(2r-11)!!
+i 2'+i
-ι/
π
^ αΐ τ + ι ·
where ( 2 r - 1 ) ! ! = 1-3-5- . . . ( 2 r - l ) . If r = 0, then
/ 0= l V W a ) . If Λ = 2 r + l , then The same integral from — oo to -f oo is zero if η = 2r-f-1 and twice the integral from 0 to oo if η = 2r.
86
The Gibbs
Distribution
In particular, if η is even ( = 2r), then r
IF = ( 7 > i ) ( 2 r + l ) ! ! ; if η = 2 r + l , then ^
+ 1
=
2 /2T\ -γ/π \ m /
( 2 r + 12 ) /
(r+l)L
- 4 - | — )
PROBLEM 2. Find the mean square fluctuation of the velocity. SOLUTION. The result of Problem 1 with η = 1 and η = 2 gives ((Avf)
= (Γ//π)(3-8/π).
= v*-v~
PROBLEM 3. Find the mean energy, the mean square energy, and the meiui square fluctuation of the kinetic energy of an atom. SOLUTION. From the results of Problem 1 we find έ=
= 3Γ/2,
\mtf ε
2
=
2
2
15Γ /4, 2
{(Δε) ) = 3Γ /2. PROBLEM 4 . Find the probability distribution for the kinetic energy of an atom. SOLUTION. d H ,
< v(!^ =
/ V e d e
-
PROBLEM 5 . Find the probability distribution for the angular velocities of rotation of molecules. SOLUTION. Just as for translational motion, we can write the probability distribution for the rotation of each molecule separately (in classical statistics). The kinetic energy of rotation of a molecule regarded as a rigid body (which is permissible, owing to the smallness of the atomic vibrations within the molecule) is
^, = | ^ (
+
«
+
« = i(f f f). +
+
where Il9 7 2, 7 3 are the principal moments of inertia, Ωΐ9Ω2, Ω3 are the components of the angular velocity along the principal axes of inertia, and Ml = ΙχΩΐ9 M2 = Ι2Ω2, M 3 = / 3£2 3 are the components of the angular momentum, which act as generalised momenta with respect to the velocities Ωΐ9 Ω2, Ωζ. The normalised probability distribution for the angular-momentum components is dwM = ( 2 * 7 T W ^ W - W exp [ - -L
+ *L + ^
jd
Ml
M 2 άΜΆ
d
id for the angular velocity dw0 = ( 2 π 7 Τ
3/ 2
( W /
exp [ - JL { ψ \ + ψ \ + ά Ω
ι
d i 3 2 dÛs.
§ 30
The Probability Distribution for an Oscillator
87
PROBLEM 6. Find the mean squares of the absolute magnitudes of the angular velocity and angular momentum of a molecule. SOLUTION. The above distributions give
M
2
= 7X7! + / , + / 3) .
§ 30. The probability distribution for an oscillator Let us consider a body whose atoms are executing small oscillations about some equilibrium positions. They may be atoms in a crystal or in a gas molecule; in the latter case the motion of the molecule as a whole does not affect the oscillations of the atoms within it and so does not influence the results. As we know from mechanics, the Hamiltonian (the energy) of a system consisting of an arbitrary number of particles executing small oscillations can be written as a sum : a
where qa are what are called the normal coordinates of the oscillations (equal to zero at points of equilibrium), pa = ςΛ are the corresponding generalised momenta, and ω α are the oscillation frequencies. In other words, E(p, q) is a sum of independent terms, each corresponding to a separate normal oscillation (or, as we say, to an oscillator). In q u a n t u m mechanics the same is true of the Hamiltonian operator of the system, so that each oscillator is independently quantised and the energy levels of the system are given by the sums
ΣΑω(Ηα + -ΐ), β
α
the na being integers. As a result of these facts the Gibbs distribution for the whole system is a product of independent factors, each giving the statistical distribution for a separate oscillator. In consequence we shall consider a single oscillator in what follows. Let us determine the probability distribution for the coordinate q of an oscillator;* the suffix α which gives the number of the oscillator will be omitted henceforward. In classical statistics the solution to this problem would be 2 2 very simple: since the potential energy of the oscillator is \co q , the probability distribution is dwq =
2T
Ae-<*W dq9 1 2
t The normal coordinate has the dimensions cm-g ^ S P - D*
.
88
The Gibbs
Distribution
or, determining A from the normalisation condition,
^=vfe "" e
W d
*
;
(3(U)
the integration with respect to q may be taken from — OO to -F <», since the integral is rapidly convergent. Let us now consider the solution of this problem in the q u a n t u m case. Let y>J(q) be the wave functions of the stationary states of the oscillator, corresponding to the energy levels εη = tw(n + \). If the oscillator is in the nth state, the quantum probability distribution for its coordinate is given by y% (in the present case the functions ψη are real, 2 and so we write simply y% instead of the squared modulus l y j ) . T h e required statistical probability distribution is obtained by multiplying yfc b y the probability wn of finding the oscillator in the wth state, and then summing over all possible states. According to the Gibbs distribution, e T
w„ = ae~ ^
9
where a is a constant. Thus we have the formula f
dwq=adq
V7
(30.2)
n«0
which is, of course, entirely in agreement with the general formula (5.8). T o calculate the sum, we can proceed as follows. With the notation dwq = qq dq, we form the derivative
dq
n~o
r
dq
Using the momentum operator ρ = ih d/dq and the fact that the oscillator m o m e n t u m has non-zero matrix elements (see Quantum Mechanics, § 23) only for transitions with η n± 1, we can write άψη
i
Λ
(Pn-l, ηψη-Ι 4"Ρ,Ι+1, ηψη+l)
=
ω
f
λ
§ 30
89
The Probability Distribution for an Oscillator
Here we have used the relations
between the m o m e n t u m and coordinate matrix elements. Thus
In the first sum we change the summation suffix from η to η +1 and use the relations ε,,+ι = εη + Λω,
tf-i,o
q»+i,„ = qn,n+u
= 0,
obtaining
In an exactly similar manner we can prove that OO T
qqq = a(\ +e-**' )
£ Π«=0
ίη.^ιψηψη^^.
A comparison of the two equations gives dp9
/2ω
whence
Λω\
{
0 CO
ÄCÜ ]
— ^2— tanh -^τ?· Determining the constant from the normalisation condition, we finally obtain the formula dw, = ί ^
tanh — )
exp J _
-
tanh
I df
() 3
( F . Bloch, 1932). Thus in the quantum case also the probabilities of various values of the coordinate of an oscillator are distributed according to a law 2 of the form exp ( — a t ? ) , but the coefficient a differs from that in the classical case. In the limit hco « Γ , where the quantisation is unimportant, formula (30.3) becomes (30.1), as we should expect. In the opposite limiting case hco » 7 \ formula (30.3) becomes d Wq
=
][nh
CP X
q2(
^~ °^
>
dq
a
3
90
The Gibbs
Distribution
i.e. the purely quantum probability distribution for the coordinate in the ground state of the oscillator.* This corresponds to the fact that when Τ
5}*
(300
In the limit of classical mechanics (fao <
2Τ
dp.
(30.5)
PROBLEM
Determine the coordinate density matrix for a harmonic oscillator. SOLUTION. The coordinate density matrix for an oscillator in statistical equilibrium is q(q,q')=a
£
T
e-°»t y)n(q')y)n(q);
see the footnote following (5.4). We put q = r-f s, q' = r—s9 and calculate the derivative (dqlds)r. As in the similar calculation above, we find 3ρ 3$
3ρ dq
3ρ dq'
= - ^ ( 1 + *-*^) £
c7w
[^
t e ) r f ) - ^ ) ^ +l ( ^ ) ] .
> w + 1 + 1
On calculating sq = \(q—q*)q in the same way and comparing it with the derivative, we obtain (dq/ds), = - sq(2colfi) coth (fao/lT), whence 2,ω . fao) q(qt «7Ο = A(r) exp ^ - J ^ - coth .
t This is the squared modulus of the wave function of the ground state of the oscillator.
§ 31
91
The Free Energy in the Gibbs Distribution
The function A(r) is determined by the condition that for s = 0, i.e. q = q' = r, the "diagonal elements" q(q, q) of the density matrix agree with (30.3). The final result is ^ , , ) = ( - t a n h - )
exp{
^ t a n h -
j j - c o t h - j .
§ 31. The free energy in the Gibbs distribution According to formula (7.9) the entropy of a body can be calculated as the mean logarithm of its distribution function :
S = -
Wn
(31.1)
a n d this is the form most frequently used. The same method gives in the classical case, using (7.12), the expression
rf*-*(*«>i/7.
ρ = (2π*)-*
(31.2)
The normalisation condition for the distribution (31.1) is Ywn
ΟΓ
= ^ Y e - ^ r =
η
1,
η
(31.3)
F^-Tlog^e-VT. η
This formula is fundamental in thermodynamic applications of the G i b b s distribution. It affords, in principle, the possibility of calculating the thermodynamic functions for any body whose energy spectrum is known. T h e sum in the logarithm in (31.3) is usually called the partition function (or sum over states). It is just the trace of the operator exp (—A/T) 9 where Ά f is the Hamiltonian of the b o d y : E
Ζ ΞΞ £ e- »!
T
= tr exp (-ß/T).
(31.4)
n * In accordance with the general rules, exp (—6IT) denotes an operator whose eigenfunctions are the same as those of the operator Û and whose eigenvalues are E T e- *> .
The Gibbs
92
Distribution
This notation has the advantage that any complete set of wave functions may be used in order to calculate the trace. A similar formula in classical statistics is obtained from the normalisation condition for the distribution (31.2). First of all, however, we must take account of the following fact, which was unimportant so long as we were discussing the distribution function as such and not relating the normalisation coefficient to a particular quantitative property of the body, viz. its free energy. If, for example, two identical atoms change places, then afterwards the microstate of the body is represented by a different phase point, obtained from the original one by replacing the coordinates and momenta of one atom by those of the other. On the other hand, since the interchanged atoms are identical, the two states of the body are physically identical. Thus a number of points in phase space correspond to one physical microstate of the body. In integrating the distribution (31.2), however, each state must of course be 1 taken only once. " In other words, we must integrate only over those regions of phase space which correspond to physically different states of the body. This will be denoted by a prime to the integral sign. Thus we have the formula F =
-riog f
é?-*M>/rdT;
(31.5)
here and in all similar cases below, d T denotes the volume element in phase space divided by (2nhy: dT= dpdq/(2nhy. (31.6) T h u s the partition function in the quantum formula (31.3) becomes an integral over states. As already mentioned in § 29, the classical energy E(p, q) can always be written as the sum of the kinetic energy K(p) and the potential energy U(q). The kinetic energy is a quadratic function of the momenta, and the integration with respect to the latter can be effected in a general form. The problem of calculating the partition function therefore actually reduces Uiq),T with respect to the coordinates. to that of integrating the function e~ In the practical calculation of the partition function it is usually convenient t o extend the region of integration and include an appropriate correction t This becomes particularly evident if we consider the classical partition function (integral over states) as the limit of the quantum partition function. In the latter the summation is over all the different quantum states, and there is no problem (remembering that, because of the principle of symmetry of wave functions in quantum mechanics, the quantum state is unaffected by interchanges of identical particles). From the purely classical viewpoint the need for this interpretation of the statistical integration arises because otherwise the statistical weight would no longer be multiplicative, and so the entropy and the other thermodynamic quantities would no longer be additive.
§31
The Free Energy in the Gibbs
Distribution
93
factor. F o r example, let us consider a gas of Ν identical atoms. Then we can integrate with respect to the coordinates of each atom separately, extending the integration over the whole volume occupied by the gas; but the result must be divided by the number of possible permutations of Ν atoms, which is N\. In other words, the integral J' can be replaced by the integral over all phase space, divided by N\ : (31.7) Similarly, it is convenient to extend the region of integration for a gas consisting of Ν identical molecules: the integration with respect to the coordinates of each molecule as a whole (i.e. the coordinates of its centre of mass) is carried out independently over the whole volume, whilst that with respect to the coordinates of the atoms within the molecule is carried out over the " v o l u m e " belonging to each molecule (i.e. over a small region in which there is an appreciable probability of finding the atoms forming the molecule). Then the integral must again be divided by N\.
PROBLEMS
PROBLEM 1. The potential energy of the interaction between the particles in a body is a homogeneous function of degree η in their coordinates. Using similarity arguments, determine the form of the free energy of such a body in classical statistics. SOLUTION. In the partition function Z=
tftrtl/TdT,
j' η/2
we replace each q by Xq and eachη ρ by λ ρ, where λ is an arbitrary constant. If at the same time we replace Τ by λ Τ, the integrand is unchanged, but the limits of integration with respect to the coordinates are altered: the linear size of the region 3 of integration is multiplied by l/λ, and so the volume is multiplied by l/λ . In order to restore the limits of integration, we must therefore at the same time replace a N ( 1 + n2 ^> V by 7?V. The result of these changes is to multiply the integral by A because of the change of variables in άΓ (s = 3N coordinates and the same number of momenta, Ν being the number of particles in the body). Thus we conclude that the substitutions V 7?Vy Τ ATgive 3ΛΓ 1+η
Ζ -+ λ < ^>Ζ. The most general form of function Z(Vt T) having this property is
where / i s an arbitrary function of one variable.
94
The Gibbs
Distribution
Hence we find for the free energy an expression of the form 3/η
F=-3(\+lln)NTlog
(1)
Τ+ΝΤφ(νΤ- /Ν),
which involves only one unknown function of one variable; the number Ν is included in the second term in (1) so that F shall have the necessary property of additivity. PROBLEM 2. Derive the virial theorem for a macroscopic body for which the potential energy of interaction of the particles is a homogeneous function of degree η in their coordinates. SOLUTION. Following the derivation of the virial theorem in mechanics (see Mechanics, § 10), we calculate the time derivative of the sum27r^p, where r and ρ are the radius vectors and momenta of the particles in the body. Since r = dK(p)/dp and K(p) is a homogeneous function of degree two in the momenta, we have ^-27r.p = 27p. at
*Β?ί+Στ·ρ = dp
2Κ(ρ)+Στ.ρ.
The particles in the body execute a motion in a finite region of space with velocities which do not become infinite. The quantity 27r«p is therefore bounded and the mean value of its time derivative is zero, so that 2Κ+(Στ·ρ)
= 0,
where Κ = K(p). The derivatives ρ are determined by the forces acting on the particles in the body. In summing over all particles we must take into account not only the forces of interaction between the particles but also the forces exerted on the surface of the body by surrounding bodies: <2r.p> = - ^ Σ τ .
-P j> r-df =
-nU-3PV;
the surface integral is transformed to a volume integral and we use the fact that div r = 3. Thus we have 2K—nU—3PV = 0 or, in terms of the total energy Ε = U+K. (n+2)K = nE+3PV. (2) This is the required theorem. It is valid in both classical and quantum theory. In the classical case, the mean kinetic energy Κ = 3NT/2, and (2) gives E+OlnyPV = 3 ( | + l/n)NT.
(3)
This formula could also be derived from the expression (1) for the free energy (Problem 1). When the particles interact by Coulomb's law (n = — 1), we have from (2) K=
-E+3PV.
This is the limiting case of the relativistic relation 2
2
E—3PV— Emc ^ ( 1 - z^/c ), in which the energy Ε includes the rest energy of the particles in the body; see Fields,§ 35.
§32
Thermodynamic
Perturbation
Theory
95
§ 32. Thermodynamic perturbation theory In the actual calculation of thermodynamic quantities there occur cases where the energy E(p9 q) of a body contains relatively small terms which may be neglected to a first approximation. These may be, for instance, the potential energy of the particles of the body in an external field. T h e conditions under which such terms may be regarded as small are discussed below. In these cases a kind of "perturbation t h e o r y " may be employed for the calculation of the thermodynamic quantities (R. E. Peierls, 1932). W e shall first show how this is to be done when the classical Gibbs distribution is applicable. We write the energy E(p9 q) in the form E(p9 q)
=
(32.1)
q) + V(p9 q)9
£ 0( / 7 ,
where V represents the small terms. T o calculate the free energy of the b o d y , we p u t g—F/T
= J*
-lE0(ptq)+V(ptq))/Tdp
e
s j e - i w ^ ^ + ^ d ^ ;
(32.2)
in the expansion in powers of V we shall always omit terms above the second order, in order to calculate the corrections only to the first and second orders of approximation. Taking logarithms and again expanding in series, we h a v e to the same accuracy F=
F
Fo + J
E
/T
^V-^yi °- ^M
^—L-jJ*
dT
άΓ
ye[F0-Eo(ptq)VT
where F 0 denotes the " u n p e r t u r b e d " free energy, calculated for V = 0. The resulting integrals are the mean values of the corresponding quantities over the " u n p e r t u r b e d " Gibbs distribution. Taking the averaging in 2
this sense and noticing that V —V
2
2
we have finally
= ((V—V) )9
F = Fo + P—-L
((V-
2
V) ).
(32.3)
Thus the first-order correction to the free energy is just the mean value of the energy perturbation V. The second-order correction is always negative, and is determined by the mean square of the deviation of V from its mean value. In particular, if the mean value V is zero, the perturbation reduces the free energy.
96
The Gibbs
Distribution
A comparison of the terms of the second and first orders in (32.3) enables us to ascertain the condition for this perturbation method to be applicable. Here it must be remembered that both the mean value V and the mean 2 square ((V— V) ) are roughly proportional to the number of particles; cf. the discussion in § 2 concerning r.m.s. fluctuations of the thermodynamic quantities for macroscopic bodies. W e can therefore formulate the desired condition by requiring that the perturbation energy per particle should be 1 small in comparison with Γ. " Let us now carry out the corresponding calculations for the q u a n t u m case. Instead of (32.1) we must now use the analogous expression for the Hamiltonian operator: + Ϋ. Η=Η0 According to the q u a n t u m perturbation theory (see Quantum Mechanics, § 38), the energy levels of the perturbed system are given, correct t o the second-order terms, by En = Ef^Vnn
(32.4)
Y^y, +
0)
where t h e £ * are the unperturbed energy levels (assumed non-degenerate); t h e prime to the sum signifies that the term with m = η must be omitted. This expression is to be substituted in the formula e-FfT
£=
ET
e' ^
η a n d expanded in the same way as above. We thus easily obtain
η
η m
- " J T Σ y*nJ»n + i
&η 1
(Σ P « * . ) >
3 2
( ·5)
where wn = exp [(F0—Ej®)/T] is the unperturbed Gibbs distribution. The diagonal matrix element Vm is just the mean value of the perturbation energy V in the given (/rth) q u a n t u m state. The sum
Ση Vnn Wn —
Vnn
is therefore the value of V averaged both over the q u a n t u m state of the body a n d over the (unperturbed) statistical distribution with respect t o the varit In expanding the integrand in (32.2) we have, strictly speaking, expanded in terms of a quantity V/T, which is proportional to the number of particles and is therefore certainly not small, but the further expansion of the logarithm causes the large terms to cancel, and so a series in powers of a small quantity is obtained.
§32
Thermodynamic
Perturbation
Theory
97
ous q u a n t u m states. This value gives the correction to the free energy in the first-order approximation, formally the same as the classical result above. F o r m u l a (32.5) may be rewritten as F=
1_ -^((Vnn-Vnnf).
.» \ Vnm?{wm-w„)
Fo + F,nn-4-ΣΣ 1
(32.6)
All the second-order terms in this expression are negative, since wm — wn has the same sign as E^—E^. Thus the correction to the free energy in the second-order approximation is negative in the q u a n t u m case also. As in the classical case, the condition for this method to be applicable is that the perturbation energy per particle should be small compared with T. O n the other hand, the condition for the applicability of the ordinary quant u m perturbation theory (leading to the expression (32.4) for En) is, as we know, that the matrix elements of the perturbation should be small compared with the separations of the corresponding energy levels; roughly speaking, the perturbation energy must be small compared with the separations of the energy levels between which allowed transitions can take placed These two conditions are not the same, since the temperature is unrelated to the energy levels of the body. It may happen that the perturbation energy is small compared with Γ, b u t is not small, or indeed is even large, compared with the significant separations between energy levels. In such cases the "perturbation t h e o r y " for thermodynamic quantities, i.e. formula (32.6), will be applicable while the perturbation theory for the energy levels themselves, i.e. formula (32.4), is n o t ; that is, the limits of convergence of the expansion r e p resented by formula (32.6) may be wider than those of (32.4), from which the former expansion has been derived. T h e converse case is, of course, also possible (at sufficiently low temperatures). F o r m u l a (32.6) is considerably simplified if n o t only the perturbation energy but also the differences between energy levels are small in comparison (E^-E^/T, with T. Expanding the difference wm - wn in (32.6) in powers of we find in this case 2
E ' < l ^ m| > + < ( F „ B- ^ ) m
2 m
m
m
t These are in general the transitions in which the states of only a small number of particles in the body are changed.
98
The Gibbs
Distribution
and we obtain an expression which is formally exactly the same as formula (32.3). Thus in this case the q u a n t u m formula is in formal agreement with the 1 classical formula. "
§ 33. Expansion in powers of h F o r m u l a (31.5) is essentially the first and principal term in an expansion of the quantum formula (31.3) for the free energy in powers of h in the quasiclassical case. It is of considerable interest to derive the next non-vanishing term in this expansion (E. P . Wigner, G . E. Uhlenbeck and L. Gropper, 1932). The problem of calculating the free energy amounts to that of calculating the partition function. F o r this purpose we use the fact that the latter is the trace of the operator exp (—ßß) (see (31.4)), with the notation β = IJT in order to simplify the writing of the involved expressions. The trace of an operator may be calculated by means of any complete set of orthonormal wave functions. F o r these it is convenient to use the wave functions of free motion of a system of Ν non-interacting particles in a large b u t finite volume V. These functions are (33.1) where the qt are the Cartesian coordinates of the particles and the p{ the corresponding momenta, labelled by the suffix /, which takes the values 1,2, . . . , s, where s = 3N is the number of degrees of freedom of the system of Ν particles. The subsequent calculations apply equally to systems containing identical particles (atoms) and to those where the particles are different. In order to allow in a general manner for a possible difference between the particles, we shall add to the particle mass a suffix indicating the degree of freedom: m{. Of course the three mi corresponding to any one particle are always equal. The existence of identical particles in a body means that, in the quantum theory, exchange effects must be taken into account. This means, first of all, that the wave functions (33.1) must be made symmetrical or antisymmetrical in the particle coordinates, depending on the statistics obeyed by the particles. It is found, however, that this effect leads only to exponentially small terms in the free energy, and so is of n o interest. Secondly, the identity of particles in q u a n t u m mechanics affects the manner in which the summation t The more powerful methods of the "diagram technique", allowing the whole perturbation theory series to be considered for thermodynamic quantities, will be described in Part 2.
§33
Expansion in Powers of h
99
over different values of the particle momenta must be carried out. We shall meet this later, for example in calculating partition functions for an ideal quantum gas. The effect produces a term of the third order in h in the free 2 energy (as shown later) and so again does n o t affect the terms of order h which we shall calculate here. Thus the exchange effects can be ignored in the calculation. In each of the wave functions (33.1) the momenta p, have definite constant values. The possible values of each pt form a dense discrete set (the distances between neighbouring values being inversely proportional to the linear dimensions of the volume occupied by the system). The summation of the matrix elements exp (—ßä)pp with respect to all possible values of the momenta may therefore be replaced by integration with respect to ρ (dp = dpx dp2... dps), bearing in mind t h a t the number of q u a n t u m N states "belonging" to the volume V dp of phase space (all values of the coordinates of each particle in the volume V and values of the momenta in dp) is N
V dp!(2nhy. We shall use the notation / = exp
[-(//i)E j Ml
exp(-^)exp
[ ( / / * ) Ç M, J .
(33.2)
The required matrix elements are obtained by integrating with respect to all the coordinates : exp ( - / ? # ) , , = ^
J/d?.
(33.3)
The partition function is then obtained by integration with respect to the momenta. Altogether, therefore, we must integrate I over all phase space, or more precisely over those of its regions which correspond to physically different states of the body, as explained in § 31. This is again denoted by a prime to the integral sign :
Z = J > - Ä = J'/dr. η
(33.4)
Let us first calculate / by means of the following procedure. W e take the derivative | g - = - e x p [ - ( / / * ) T M i] Ä { « p
KWJ»*,]/},
100
The Gibbs
Distribution
the operator Ü acting on all factors to its right, and expand the right-hand side, using the explicit expression for the Hamiltonian of the b o d y :
where U = U(ql9 q2, ...» q5) is the potential energy of interaction between all particles in the body. By means of (33.5) we obtain after a straightforward calculation the following equation for / :
ν & l
2i
8 7
Γ·/
, \
87
927
where
'
E(p 9)
= %-êr
*
+u
(33
i
6)
is the usual classical expression for the energy of the body. This equation is to be solved with the obvious condition that I — I when 0 = 0. The substitution / =
(33.7)
-ßE(Pe,g)x
gives
II. op
Σ
» . r_ ψ
τ 2m, [ -
2
β
^
*_ g Ä
η ^
+
aqt Μ
ft
öqi
+ f t
oqf
\
(»z)' \
8)
with the boundary condition χ = 1 for β = 0. In order to obtain an expansion in powers of ή, we solve equation (33.8) by successive approximations, putting X = l + h
X l
+ h*x2+
.•·,
(33.9)
with χ ι = 0, χ 2 = 0, . . . for β = 0. Substituting this expansion in equation (33.8) and separating terms in different powers of h, we obtain the equations
§33
Expansion in Powers of h
101
The first equation gives χ „ and then the second equation gives χ 2. A simple calculation leads to the results
The required partition function (33.4) is Ζ = J' (1 +
2
+ Α χ 2> - ^ < * . *> dr.
(33.11)
T h e term of the first order in A in this integral is easily seen to be zero, since ßEOfq) the integrand %ie~ ' in that term is an odd function of the m o m e n t a (E(p, q) being quadratic in the momenta and χι, by (33.10), linear), a n d s o the result on integrating with respect to momenta is zero. Thus we can write (33.11) as
Z= where
(l+*%)J'e-«*«dr,
is the value of χ 2 averaged over the classical Gibbs distribution : y~ —
Substituting this expression for the partition function in formula (31.3), w e have for the free energy F=
F c l- - l o g ( l + A * 5 Q > ß
or, to the same accuracy, F=Fcl-h%/ß.
(33.12)
H e r e F d is the free energy in classical statistics (formula (31.5)). T h u s the next term after the classical expression in the expansion of t h e free energy is of the second order in h. This is n o t accidental: in equation (33.8), solved here by the method of successive approximations, the q u a n t u m constant appears only as /A, and so the resulting expansion is one in powers of ih\ b u t the free energy, being a real quantity, can contain only powers of ih which are real. Thus this expansion of the free energy (ignoring exchange effects) is an expansion in even powers of A. It remains to calculate the mean value χ^. W e have seen in § 29 t h a t in classical statistics the probability distributions for the coordinates a n d
102
The Gibbs
Distribution
momenta are independent. The averaging over momenta and over coordinates can therefore be made separately. The mean value of the product of two different momenta is clearly zero. The mean value of the square pf is mjß. W e can therefore write (rnilß)bik,
(PiPk) =
where àik = 1 for i — k a n d 0 for i ^ k. Having averaged with respect to momenta by means of this formula, we obtain x
r
-
Z
v
J
/
^
-
V
^
Y
-
^
(3313)
T h e two terms here may be combined, since the mean values are related by the formula
dqf This is easily seen by noticing that
T h e first term on the right-hand side gives only a surface effect, a n d since the body is macroscopic this effect may be neglected in comparison with the second term. Substituting the resulting expression for %2 in formula (33.12), and replacing β by 1/Γ, we find the following final expression for the free energy: F
=
1 F
c *2
24Γ
Yi
±
+
(
m
\
r "!/V 3?i7
(33.15)
We see that the correction to the classical value is always positive, a n d is determined by the mean squares of the forces acting on the particles. This correction decreases with increasing particle mass a n d increasing temperature. According to the above discussion, the next term in the expansion given here would be of the fourth order. This enables us to calculate quite independently the term of order #* which occurs in the free energy because of the peculiarities of the summation over momenta resulting from the identity of particles in quantum mechanics. T h e term in question is formally the same as the correction term which appears in a similar calculation for an ideal gas, a n d is given by formula (56.14): = + —- 2g
— VTWrrfiK
f33 16) ÔA0 V >
103
Expansion in Powers of h
§33
for a body consisting of Ν identical particles. The upper sign applies for Fermi statistics and the lower sign for Bose statistics; g is the total degree of degeneracy with respect to the directions of the electron and nuclear angular momenta. F r o m these formulae we can also obtain the correction terms in the p r o b a bility distribution functions for the coordinates and momenta of the atoms of the body. According to the general results i n § 5, the m o m e n t u m probability distribution is given by the integral of / with respect to q (see (5.10)): dwp = constant Xàp j I àq. ßE{pq)
The term %\e~~ ' in / contains a total derivative with respect to the coordinates, and the integral of it gives a surface effect which can be neglected. Thus we have 2m
βυd
à
P JO
dwp = c o n s t a n t X e x p ^-ß^P^ ^
+&%2)*~ <2-
The third and fourth terms in the expression (33.10) for χ 2 give a small constant (not involving the momenta) o n integration, and this can be neglected ßu in the same approximation. Taking o u t also the factor j e~ dq and including it in the constant coefficient, we have
6
mmk
dqidqk\
The mean values which appear here are related by
~dHT ^ 3?i
dqk ~
~W~W dqi
dqk '
similarly to (33.14). Hence
= „„s „ xexp M
t
Σ
ξ
Α (33.17)
This expression can be conveniently rewritten in the following final form : «hr, = c o n s.t a n^ t x e x p |ί -
1
Γ ν T
A
&
V
V
PiPk
W — W _ J Ij I .ap, (33.18)
the bracket in (33.17) being replaced by an exponential function to the same degree of accuracy.
104
The Gibbs
Distribution
Thus we see that the correction to the classical distribution function for the momenta is equivalent to adding to the kinetic energy in the exponent an expression quadratic in the momenta, with coefficients depending on the law of interaction between the particles in the body. If it is desired to find the probability distribution for any one m o m e n t u m pi9 then (33.17) must be integrated with respect to all the other momenta. All the terms involving the squares pi (k ?± i) will then give constants negligible compared with unity, while the terms containing products of different momenta give zero. The result is, again in exponential form, d » A = constantxexp { - - ^
[l-
} öp,.
(33.19)
T h u s the distribution obtained differs from the Maxwellian only in t h a t the true temperature Τ is replaced by a somewhat higher "effective temperature":
Similarly we can calculate the corrected coordinate distribution function b y integrating / with respect to the m o m e n t a : dwq = constant X dq j I dp. The same calculations as led to (33.13) give .
. v
1
f
l"rr
1
** ^
& V
/MY
dwg = c o n s t a n t x e x p | - -
+
1
2
9 ^ll A ^j}
dq. (33.20)
§ 34. The Gibbs distribution for rotating bodies The problem of the thermodynamic relations for rotating bodies has already been considered in § 26. Let us now see how the Gibbs distribution is to be formulated for rotating bodies. This will complete the investigation of their statistical properties. As regards the uniform translational motion, Galileo's relativity principle shows that, as already mentioned in § 26, this motion has only a trivial effect on the statistical properties and so needs n o special consideration. In a system of coordinates rotating with the body, the usual Gibbs distribution is valid; in classical statistics, E
T
ρ = (27dt)-°eV'- '<
<»V 9
(34.1)
Ç 34
105
The Gibbs Distribution for Rotating Bodies
where £"(/;, q) is the energy of the body in this system, as a function of the coordinates and momenta of its particles, and F the free energy in the same system (which, of course, is not the same as the free energy of the body when a t rest). T h e energy E'(p, q) is related to the energy E(p, q) in a fixed system by q), (34.2) E'(p, q) = E(p, q)-Q.M{p9 where Ω is the angular velocity of rotation and M(p, q) the angular m o m e n t u m of t h e body (see § 26). Substituting (34.2) in (34.1), we find the G i b b s f distribution for a rotating body in the f o r m ρ = (2nh)~
s
eir-*U>. *)+ß-M
(343)
In classical statistics the Gibbs distribution for a rotating body can also be represented in another form. T o obtain this, we use the following expression for the energy of the body in the rotating coordinate system: 2
E' = E\mvf -\Em(S}Xtf
(34.4)
+ U,
where the v' are the velocities of the particles relative to the rotating system, a n d the r their radius vectors (see Mechanics, § 39). Denoting by E0(y\
r) == E\mtf
2
(34.5)
+U
the p a r t of the energy which is independent of Ω, we obtain the G i b b s distribution in t h e f o / m ρ = ( 2 a * ) - * e x p Ji
[f'-^,
T)+^Em<ßXtf
]
T h e function ρ determines the probability corresponding to the element of phase space dx1 dy± dz1 . . . dplx dp'ly dp'lz . . . , where p ' = w y ' + m f l X r (see Mechanics, § 39). Since, in obtaining the differentials of the momenta, we must regard the coordinates as constant, d p ' = m dv', a n d the probability distribution expressed in terms of the coordinates and velocities of the particle is dw = C e x p J £ - I Xdxxdy1dzx
| £ 0( v ' ,
r)-Ejm(QXr)
. . . dv'lxdv'^dv'lz
(34.6)
5
where C denotes for brevity the factor (2nh)~ together with the p r o d u c t of the particle masses which appears when we go from the m o m e n t u m differentials to the velocity differentials. * The distribution (34.3), like the ordinary Gibbs distribution, is fully in agreement with the result (4.2) derived in § 4 from Liouville's theorem: the logarithm of the distribution function is a linear function of the energy and angular momentum of the body.
106
The Gibbs
Distribution
F o r a body at rest we have dw =
C J T - w . d v t
d
xi
di
d Z l . . . évlx
y
r
dvly
dvlz
. . ( 3 4 . 7 )
n ow a
function of the velocities in with the same expression (34.5) for 2s 0(v, ) > the fixed coordinate system. Thus we see that the Gibbs distribution for the coordinates and velocities for a rotating body differs from that for a body at 2 rest only by the additional potential energy — y 2 m ( ß X r ) . In other words, as regards the statistical properties of the body, the rotation is equivalent to the existence of an external field corresponding to the centrifugal force. T h e statistical properties are not affected by the Coriolis force. It should be emphasised, however, that this last result applies only to classical statistics. In the quantum case the expression w = exp
Q.U)IT]
(34.8)
gives the statistical operator corresponding to (34.3) for a rotating body. Formally we can reduce this operator to a form analogous to (34.6), the velocities v' being replaced by the operators v' = p'/m—Ωχτ, but the components of this vector operator d o not commute, unlike those of the operator ? in the fixed system. The statistical operators corresponding to the expressions (34.6) and (34.7) will therefore in general be markedly different from each other, quite apart from the fact that one of them contains the centrifugal energy.
§ 35. The Gibbs distribution for a variable number of particles So far we have always tacitly assumed that the number of particles in a body is some given constant, and have deliberately passed over the fact that in reality particles may be exchanged between different subsystems. In other words, the number Ν of particles in a subsystem will necessarily fluctuate about its mean value. In order to formulate precisely what we mean by the number of particles, we shall use the term subsystem to refer to a part of the system which is enclosed in a fixed volume. Then Ν will denote the number of particles within that volume.* Thus the problem arises of generalising the Gibbs distribution t o bodies with a variable number of particles. Here we shall write the formulae for bodies consisting of identical particles; the further generalisation t o systems containing different particles is obvious (§ 85).
t In deriving the Gibbs distribution in§ 28 we have in essence already understood subsystems in this sense; in going from (28.2) to (28.3) we differentiated the entropy whilst regarding the volume of the body (and therefore of the medium) as constant.
§ 35
The Gibbs Distribution for a Variable Number of Particles
107
The distribution function n o w depends n o t only on the energy of the quantum state b u t also on the number Ν of particles in the body, a n d the energy levels EnN are of course themselves different for different Ν (as indicated by the suffix N). The probability that the body contains Ν particles a n d is in the nth state will be denoted by w^. The form of this function can be determined in exactly the same way as the function wn in § 28. T h e only difference is that the entropy of the medium is now a function n o t only of its energy E' b u t also of the number N' of particles (0) in i t : S' = S\E\ N"). Writing Ε = E^-EnN a n d N' = N —TV (where Ν is (0) the number of particles in the body, a n d N the given total number of particles in the entire closed system, which is large compared with N), we have in accordance with (28.2) wnN =
œnstàntXQx^{S\E^-EnN,N^-N)}\
the quantity ΔΕ' is regarded as constant, as in § 28. Next, we expand S" in powers of E^ a n d N, again taking only the linear terms. Equation (24,5), in the form + *
shows that (dS/dE)yN S'(E«»-EnN,
ON,
= 1/Γ, (dS/dN)EtV
-
NM-N)
= -μ/Τ.
tf<°>)-^
Hence
+^ ,
the chemical potential μ (and the temperature) being the same for the body and the medium, from the conditions of equilibrium. Thus we obtain for the distribution function the expression wnN
N E
= Aei« - »*VT.
(3
)
5 > 1
The normalisation constant A can be expressed in terms of the thermodynamic quantities in the same way as in § 3 1 . T h e entropy of the body is S = - ( l o g wflN)
= -\ogA-^r
+γ>
and so TlogA
=
E-TS-μΝ.
But E—TS= F , and the difference F— μΝ is the thermodynamic potential Ω. Thus riog A = .Q, a n d (35.1) may be rewritten as (35.2)
108
The Gibbs
Distribution
This is the final formula for the Gibbs distribution for a variable number of particles.* The normalisation condition for the distribution (35.2) requires that the result of summing the wnN first over all quantum states (for a given N) and then over all values of Ν should be equal to unity:
£Σ>«* = ^ Ç ( ^ Σ ' / r
/ r
e
EnNlT
)=
L
Hence we obtain the following expression for the thermodynamic potential Ω : Ω = - Τ log Ç
Ç
W
J.
(35.3)
This formula together with (31.3) can be used to calculate the thermodynamic quantities for specific bodies. Formula (31.3) gives the free energy of the body as a function of Γ , Ν and V9 and (35.3) gives the potential Ω as a function of Γ , μ and V. In classical statistics the probability distribution has the form dwN =
QNdpWâtfM9
where = {2nh)~* 4Ρ+μΝ-Ε*ίΡ. 4W m (354) The variable Ν is written as a subscript to the distribution function, and the same letter is written as a superscript to the element of phase volume in order t o emphasise that a different phase space (of 2s dimensions) corresponds t o each value of N. The formula for Ω correspondingly becomes ρΝ
fl = - Γ log
evWTf
E
e~ »<
τ
™ άΓΝJ.
(35.5)
Finally, we may say a few words concerning the relation between the G i b b s distribution (35.2) for a variable number of particles derived here and the previous distribution (31.1). First of all, it is clear that, for the determination of all the statistical properties of the body except the fluctuations of the total number of particles in it, these two distributions are entirely equivalent. On neglecting the fluctuations of the number N9 we obtain Ω -h μΝ = F , and the distribution (35.2) is identical with (31.1). The relation between the distributions (31.1) and (35.2) is to a certain extent analogous to that between the microcanonical and canonical distributions. The description of a subsystem by means of the microcanonical distribution is equivalent to neglecting the fluctuations of its total energy; the canonical distribution in its usual form (31.1) takes into account these fluct This distribution is sometimes called the grand canonical ensemble.
§36
The Derivation of the Thermodynamic
Relations
109
tuations. The latter form in turn neglects the fluctuations in the number o f particles, and may be said to be "microcanonical with respect to the number of particles'*; the distribution (35.2) is "canonical" with respect to both the energy and the number of particles. T h u s all three distributions, the microcanonical and the two forms of the Gibbs distribution, are in principle suitable for determining the thermodynamic properties of the body. The only difference from this point of view lies in the degree of mathematical convenience. In practice the microcanonical distribution is the least convenient and is never used for this purpose. The Gibbs distribution for a variable number of particles is usually the most convenient.
§ 36. The derivation of the thermodynamic relations from the Gibbs distribution The Gibbs distribution plays a fundamental part throughout statistical physics. We shall therefore give here another justification of it. This distribution has essentially been derived in §§ 4 and 6 directly from Liouville's theorem. We have seen that the application of Liouville's theorem (together with considerations of the multiplicativity of distribution functions for subsystems) enables us to deduce that the logarithm of the distribution function of a subsystem must be a linear function of its energy : (36.1)
log wn = *+ßEn9
the coefficients β being the same for all subsystems in a given closed system (see (6.4), and the corresponding relation (4.5) for the classical case). Hence wn =
a+ßE
e *;
using the purely formal notation β = — 1/Γ, α = F/T9 we have an expression of the same form as the Gibbs distribution (31.1). It remains to show that t h e fundamental thermodynamic relations can be derived from the Gibbs distribution itself, i.e. in a purely statistical manner. We have already seen that the quantity β9 and therefore T9 must be the same for all parts of a system in equilibrium. It is also evident that β < 0, i.e. Τ > 0, since otherwise the normalisation sum Σχνη must diverge: owing to the presence of the kinetic energy of the particles, the energy En can take arbitrarily large values. All these properties agree with the fundamental properties of the thermodynamic temperature. T o derive a quantitative relation, we start from the normalisation c o n d i tion £^-*n)/r=
1.
110
The Gibbs
Distribution
We differentiate this equation, regarding the left-hand side as a function of Γ a n d of various quantities λΐ9 λ 2, . . . which represent the external conditions to which the body considered is subject; these quantities may, for example, determine the shape and size of the volume occupied by the body. T h e energy levels En depend on λΐ9 . . . as parameters. Differentiation gives
where for simplicity only one external parameter is used. Hence
dFÇ wm = άλÇ wn
+
(f- ς wEy n
n
O n the left-hand side Ewn = 1, and on the right-hand side
Σ
w
n
_
dE„
dEn
» » = > Ç^-gr-'ar. E
E
Using also the formulae F-E
=—TS dËJdX
and
f
= 9Ä/8A,
(36.2)
we have finally dF =
-SdT+dß/dX.M.
This is the general form for the differential of the free energy. In the same way we can derive the Gibbs distribution for a variable number of particles. If the number of particles is regarded as a dynamical variable, it is clear that it will be an "integral of the m o t i o n " , and additive, for a closed system. W e must therefore write log wnN = oc+ßEn + γΝ,
(36.3)
where λ, like β, must be the same for all parts of a system in equilibrium. Putting α = Ω/Τ9 β = — l/T9 γ = μ/Τ9 we obtain a distribution of the form (35.2), and then by the same method as above we can deduce an expression for the differential of the potential Ω.
t If the Hamiltonian Û (and therefore its eigenvalues E^ depends on a parameter A, then dEjdX^idDidX)^ see Quantum Mechanics, (11.16). On statistical averaging this gives (36.2).