陕西师范大学精品课程 ……,物理化学,
第 1 页 共 40 页 2004-7-15
第三章 统计热力学基础
第一节 绪论
一、统计热力学研究的对象、任务和方法
1,对象
研究的是大量微观粒子所构成的宏观物质处于平衡态时的性质。从研究对象来讲，
统计热力学和热力学一样都是研究宏观物质处于平衡态时的性质。热力学是根据从经验归纳得到的四条基本定律，而不管物质的微观结构和微观运动形态，因此只能得到联系各种宏观性质的一般规律，而不能给出微观性质与宏观性质之间的联系。而统计热力学则是从物质的微观结构出发来了解其宏观性质。两者研究的深度不同，所以有必要讨论后者。
2,任务
研究粒子所构成的体系的宏观行为，从粒子的微观性质来寻求体系的平均的宏观性质，这就是统计热力学的任务和研究内容。由此可见，统计热力学是从微观到宏观过渡的理论。它具有统计平均的性质，是联系物质的宏观性质与微观结构、沟通热力学与量子力学的一座桥梁。
3,方法
统计力学的研究方法是微观的方法，它根据统计单位的力学性质如速度、动量、位置、振动、转动等，用统计的方法来推求体系的热力学性质，例如压力、热容、熵等热力学函数。统计力学建立了体系的微观性质和宏观性质之间的联系。从这个意义上，统计力学又可称为统计热力学。
对于简单分子，应用统计热力学的方法进行处理，其结果是令人满意的。当然统计热力学也有自身的局限性，由于人们对于物质结构的认识不断深化，不断地修改充实物质结构的模型，同时模型本身也有近似性，所以由此得到的结论也有近似性。
从历史的发展来看，最早是由玻兹曼（ Boltzmann）以经典力学为基础建立的统计方法,称为经典统计热力学。 1900 年普朗克 （ Planck） 提出了量子论,麦克斯韦 （ Maxwell）
将能量量子化的概念引入统计热力学，对经典统计进行某些修正，发展成为麦克斯韦－
玻兹曼统计热力学方法。 1924 年量子力学建立后，在统计力学中不但所依赖的力学基础要改变，而且所用的统计方法也需要改变。由此产生了玻色－爱因斯坦（ Bose-Einstein）
统计和费米－狄拉克（ Fermi-Dirac）统计，分别适用于不同的体系。这两种统计方法都陕西师范大学精品课程 ……,物理化学,
第 2 页 共 40 页 2004-7-15
可以在一定的条件下通过适当的近似而得到玻兹曼统计。本章的内容就是简要介绍麦克斯韦－玻兹曼统计热力学的基本原理和应用。
二、统计体系的分类
在统计热力学中，按照构成体系的微观粒子（称为,统计单位,）的不同特性，可以将体系分为不同的类型。
按照粒子是否可以分辨，把体系分为定位体系（ localized system） （或称为定域子体系）和非定位体系（ non-localized system） （离域子体系）,前者的粒子可以彼此分辨，
而后者的粒子彼此不能分辨。例如气体分子处于无序运动之中，彼此无法区别，因此是离域子体系。而晶体，由于粒子是束缚在晶格位置上作振动运动，每个位置可以想象给予编号而加以区别，所以晶体是定域子体系。
按照统计单位之间有无相互作用,又可以把体系分为近独立粒子体系 (assembly of
independent particles)和非独立粒子体系（ assembly of interacting particles） 。前者或简称为独立粒子体系，其粒子之间的相互作用非常微弱，可以忽略不计，如理想气体，这种体系的总能量等于各个粒子的能量之和，即 =
i
i
U ε
∑；后者或称为相依粒子体系，其粒子之间其的相互作用不容忽略，如高圧下的实际气体等，这种体系的总能量除了各个粒子的能量之和外，还存在粒子之间相互作用的位能，即
I111
= (x,y,z,......x,y,z )
iNN
i
UUε +
∑
。 显然,粒子之间绝对无相互作用的体系是不存在的，
但可以把那些粒子之间的相互作用非常微弱可以忽略不计的体系，如低圧气体，作为独立粒子体系进行处理。本章中仅限于讨论独立粒子体系。
统计力学可分为两大阶段：经典统计力学和量子统计力学。前者是在 19 世纪末发展并成熟起来。在许多场合能给出满意的结果，但某些情况下它无法解释一些实验结果。
后者在二十世纪二十年代（ 1926 年）量子力学建立后发展起来的。它比经典统计力学能解释更广泛的宏观现象。本章着重讨论经典统计力学，只对量子统计力学稍加介绍。
三．数学知识
1．排列与组合
在统计热力学中，需要讨论粒子在不同能级上的分布，这在数学上相当于排列组合问 题。因此，先扼要介绍一些排列组合的有关知识。
(1)在 N 个不同的物体中，每次取出 m 个按照一定的顺序排成一列，称为从 N 个物体中每次取 m 个物体的排列；其排列的方式数为
陕西师范大学精品课程 ……,物理化学,
第 3 页 共 40 页 2004-7-15
A
N
m
= N!/(N - m)! (3.1?1)
当 m = N 时，上述排列称为全排列，全排列的方式总数为
A
N
N
= N!/(N – N)! = N ！ ~ (3.1?2)
其中规定 0！ = 1。
(2)若在 N 个物体中有 n
1
个相同，另外 n
2
个也彼此相同，其余的各不相同，则这 N
个物体的全排列方式数为
n!/N
1
!N
2
! (3.1?3)
(3) 将 N 个相同的物体放入 M 个相同的容器中（每个容器的容量不限），则放置的方式数为
（ N+M-1） !/N!(M-1)! (3.1?4)
(4) 将 N 个不同的物体放入 M 个容器中（每个容器容量不限），则放置的方式数为
M
N
(3.1?5)
(5) 在 N 个不同的物体中，每次提取 m 个，不管排列顺序编为一组，称为从 N 个不同物体中每次取出 m 个物体的组合，其组合数为
C
N
m
= N!/m!(N-m)! (3.1?6)
2,斯特林（ stirling）公式
在统计热力学中，常常要计算 N!。
阶乘 N！可展开如下式
23
1 1 139
!2 1
e 12 288 51840
N
N
NN
NN N
π
=++?
" (3.1?7)
(1) 当 N 是不太小的整数时，可近似为
1n ! 1n 2 1n 2 1n
e
N
N
NN NNNππ
==+?
(3.1?8)
误差约为
N12
1
。
(2) 当 N 是很大的整数时，上式进一步简化为
1n ! 1n 1n
e
N
N
NNN-
==
(3.1?9)
陕西师范大学精品课程 ……,物理化学,
第 4 页 共 40 页 2004-7-15
误差约为 ()Nπ2n1
2
1
。式（ 3.1?8）和式（ 3.1?9）称为斯特林公式。即,
当 N 很大时，有
lnN! = NlnN -N
当 N 很小时,有
!2()
e
=π
N
N
NN
N 越大，所得结果越精确。
第二节 粒子体系的统计分布
一,体系的微观状态及微观状态数,
2,体系的微观状态
用微观性质来描述的体系运动状态称体系的微观状态。
体系的微观状态是由各个粒子的微观状态所决定。粒子的微观状态是瞬息万变的，
由大量粒子组成的体系的微观运动状态也是千变万化的，如何描述粒子及体系的微观运动状态呢？经典力学与量子力学有不同的描述方法。
经典力学：粒子运动遵守牛顿运动方程，常用空间坐标 (q
x
,q
y
,q
z)
、瞬时速度或动量
(p
x
,p
y
,p
z)
来描述粒子的运动状态。在经典力学中，可根据粒子的空间坐标识别它们，故在经典力学中认为粒子是可别的。
量子力学：认为微观粒子是等同的、不可区别的。同时粒子具有波粒二相性，根据测不准原理,粒子不可能有确定的坐标和动量数值,所以不能用经典力学的方法来描述。
量子力学用本征函数 φ（波函数）和相应的本征值（能量值）来描述，微观粒子的运动服从薛定谔方程，即
22 2 2
22 2 2
()(,)
8
++ + =
h φφφ
Vxyz φE
mx y zπ
通过解此方程得到粒子的波函数 φ
i
相对应的能量 E
i
，具有一定 φ
i
,E
i
的状态称一种量子状态或量子态。从能级表达式得出能量是不连续的、量子化的。粒子处于不同的能级就表现不同的状态。
3,体系的微观状态数?
在体系的体积、总能量一定的情况下，含有 N 个粒子的体系中各种分布类型的样数之和，称体系的总微观状态数。 也称为体系的热力学几率。
陕西师范大学精品课程 ……,物理化学,
第 5 页 共 40 页 2004-7-15
设每种分布的微观状态数为 t
j
,则体系的总微观状态数就 等于各种分布的微观状态数之和，即?= Σt
j
。
我们知道,密闭体系的热力学能是熵和体积的特征函数,则体系的熵可表示为 S（ U、
V,N） 。由玻尔兹曼（ Boltzmann）公式（ S = k
B
ln?）知，该体系的总微观状态数? 也可以表示为?（ U,V,N） 。这就是说当体系的热 力学性质（ U,V,N）确定时，体系的能级和能级的简并度一定，体系的总微观状态数也一定，而且巨大数目的不同的微观状态对应于一个给定的宏观状态。
熵 S（体系的状态性质）与热力学几率?（体系的状态性质）之间存在某种函数关系，这个关系可表示为,
S = f(?)
玻兹曼定理为,
S = k
B
ln?
它揭示了体系的熵函数与其热力学 几率之间的关系。可以证明
k
B
=R/L=1.3806×10
-23
J·K
1
,即 k
B
是玻兹曼常数。 这是联系微观性质与宏观性质的关系式。
只要知道体系?，由此式可求出热力学函数 S，进而还可求出其它热力学函数。所以?
是重要的热力学函数。
二,统计热力学的基本假定
统计热力学认为：,对于宏观处于一定平衡状态的体系而言，任何一个可能出现的微观状态都具有相同的数学几率。,
统计热力学的这个基本假设，就是认为在所有 可能出现的微观状态中，任何一种状态都没有明显理由比其它微观状态出现的可能性更大些，这称为,等可几假设,。
上述假定的出发点是认为体系的热力学性质是所有可能出现的微观状态的统计平均。当我们对体系进行宏观测量时,需要一定的时间，在此时间内，体系将经历所有可能的微观状态。因此，宏观测得的某个物理量实际上是相应微观量的平均值，其中每个微观状态对平均值的贡献是相同的。这个假设的合理性已经由其引出的结论与实验事实相一致而得到证明。
由上述假定，对于拥有? 个微观状态的热力学体系，每一个微观状态出现的几率应为 1/?，而某一分布 x 出现的几率则为
P
x
= t
x
/?
式中 t
x
是该分布所拥有的微观状态数。此式表明，虽然各微观状态出现的几率相同，但陕西师范大学精品课程 ……,物理化学,
第 6 页 共 40 页 2004-7-15
各种分布出现的几率是不相同的。
为了说明统计热力学的基本假定，设某系统有 4 个可辨粒子 a,b,c,d，分配于两个相连的、容积相等的空间Ⅰ和Ⅱ之中，所有可能的分配形式如表 3.1 所列。
表 3.1 可分辨粒子在两个等容积空间的分配方式
分布方式 空间Ⅰ 空间Ⅱ 微观状态数 数学几率
（ 4,0）分布 a b c d 0 C
4
4
=1 1/16
（ 3,1）分布
a b c
a b d
a c d
b c d
d
c
b
a
C
4
3
=4 4/16
（ 2,2）分布
a b
a c
a d
b c
b d
c d
c d
b d
b c
a d
a c
a b
C
4
2
=6 6/16
（ 1,3）分布
a
b
c
d
b c d
a c d
a b d
a b c
C
4
1
=4 4/16
（ 0,4）分布 0 a b c d C
4
0
=1 1/16
在表 3.1 中列出的每一种可能的分配形式称为一个微观状态，所有可能的分配形式总数称为体系的总微观状态数，用? 表示。由表 3.1 可见，所研究体系的总微观状态数
=16。
根据等可几假设，上例中每个微观状态出现的数学几率都是 1/16。
三．经典统计 —麦克斯韦－玻兹曼统计
运用麦克斯韦－玻兹曼统计法所研究的体 系，应该具有两个特性,(1) 宏观状态确定的封闭体系； (2) 近独立粒子体系。
以下讨论的麦克斯韦－玻兹曼统计体系是按粒子的量子状态的分布来处理的，即能量是量子化的。这样做的目的是使推导过程更清晰和容易理解、接受。
陕西师范大学精品课程 ……,物理化学,
第 7 页 共 40 页 2004-7-15
1,非间并度、间并度
非间并度：一种能级只有一种量子状态即一种能量对应一种波函数。
n
1
n
2 ……….
n
i
ε
1
ε
2 ………,
ε
i
φ
1
φ
2 ………
φ
i
间并度：一种能级有多种量子状态即一种能量对应多个波函数。
n
1
n
2 ……………
n
i
ε
1
ε
2
………
,
ε
i
11 12 1g
...
i
φφ φ
21 22 2g
...
i
φφ φ ………
12
...
i
ii ig
φφ φ
注,g
i
是能级 ε
i
具有的量子状态数,称该能级的间并度或者统计权重。
2．麦克斯韦－玻兹曼统计,求体系的微观状态数?
(1) 对非间并情况,首先体系 U,N,V 固定,N 个粒子可区分
例,N 个可区分粒子的一种特殊分布,
有 n
1
个粒子的能量同为 ε
1
即 n
1
个粒子分配在 ε
1
能级上
有 n
2
个粒子的能量同为 ε
2
即 n
2
个粒子分配在 ε
2
能级上
┋ ┋
有 n
i
个粒子的能量同为 ε
i
即 n
i
个粒子分配在 ε
i
能级上
示意为,
能级
ε
1
ε
2
……,ε
i
一种分布方式 n
1
n
2 ………,
n
i
另一种分布方式 n
1
′
n
2
′
……….
n
i
′
由于粒子有弹性碰撞，在运动中相互交换能量，所以 N 个粒子在满足 U,N,V 一定的条件下可以有许多种不同的分布方式。
先设其中一分布方式的微观状态数，即花样数为 t
这个问题相当于将 N 个不同的物体分成若干堆，按排列组合公式
!
!
=
∏
i
i
i
N
t
n
(3.2?1)
那么对各种分布方式，由以上公式我们都可算出它的 t。
则体系的总微观状态数? 等于各个 t 的加和,即
陕西师范大学精品课程 ……,物理化学,
第 8 页 共 40 页 2004-7-15
!
!
==
∑∑
∏
i
ii
i
i
N
t
n
(3.2?2)
(2) 间并态的情况,在 U,N,V 一定的条件下,N 个可区分粒子在各能级上的一种分布
一种分布方式,
ε
1
ε
2
………ε
i
n
1
n
2 ……….
n
i
g
1
g
2
………g
i
即,能级是间并的，在 ε
i
能级上有 g
i
个间并度，即拥有 g
i
个不同的量子状态，而且假设每个量子状态容纳的粒子数不受限制，则 n
i
个处于能级 ε
i
的粒子中,每个粒子都可能有 g
i
个分布方式，故 n
i
个粒子在 g
i
个间并度内有
i
n
i
g 个微观状态数。
那么,ε
1
能级上有 n
1
个粒子，由于有 g
1
个间并度，产生
1
1
n
g 个微观状态
ε
2
能级上有 n
2
个粒子，由于有 g
2
个间并度，产生
2
2
n
g 个微观状态
┋ ┋ ┋
ε
i
能级上有 n
i
个粒子，由于有 g
i
个间并度，产生
i
n
i
g 个微观状态
所以一种分布方式，由于全部间并分布，所产生的微观状态数为
12
=
∏
ii
nnnn
i
gg g g
故一种分布方式所有的微观状态数 t 应为非间并时的微观状态数乘以间并分布时的微观状态数，即
i
!
!
!!
==
∏∏
∏
i
i
n
n
i
ii
i
gN
tgN
nn
(3.2?3)
而总的微观状态数? 为各种分布方式的 t 之和，即
!
!
==
∑∑
i
n
i
i
ii
i
g
tN
n
(3.2?4)
这样由 (3.2?2) (3.2?4)式就可求得非间并及间并态下的体系的总微观状态数。
注意：前提是在 U,N,V 一定的条件下。
四．最可几分布
对于指定状态的宏观体系，它的各种分布所拥有的微观状态数大小不一，其中必有一种分布所包含的微观状态数最多或出现的几率最大，称为最可几分布。
由表 3.1 可以看出，尽管各微观状态具有相同的数学几率，但各种分布的数学几率是不相同的。数学几率最大的分布称为,最可几分布,，上例中（ 2,2）分布就是最可几分陕西师范大学精品课程 ……,物理化学,
第 9 页 共 40 页 2004-7-15
布，其数学几率为 6/16。统计热力学认为，最可几分布可以代表体系的平衡分布。它可以代表体系处于热力学平衡状态时的一切分布状态。也就是说，对于一个粒子数目众多的实际平衡体系而言，其微观状态虽然千变万化，但基本上都是辗转于最可几分布以及与最可几分布没有实质差别的那些分布之中。因此，最可几分布是统计热力学最关注的分布。
最可几分布的微观状态数以 t
max
表示，则
≈ t
max
即 N 足够大时，在体系的各种分布方式中，求?只需考虑 t 中的最大项 t
m
对?的贡献就可以了，其它项可略去不计。实际上,N 很大时，体系的分布方式数目大的无法精确计算，也没必要计算。 （注：最可几分布方式中各能级的分布数 n
i
都以 n
i
*
表示。 ）
也就是只要求出与 t
m
相应的的那一套分布 n
1
*
n
2
*
…n
i
*
的微观状态数,就可代表体系的一切微观状态数。
S=k
B
lnt
max
=k
B
ln?
目的：求 t
m
时的那一套分布 n
1
*
n
2
*
…n
i
*
的表达式。求最可几分布是统计热力学的核心问题。
1,两个限制条件：首先 t
m
时的那一套分布 n
1
*
n
2
*
…n
i
*
必须满足 U,N 一定，以式子表示为,
∑ n
i
*
= N （ 3.2?5）
∑ n
i
*
ε=U （3.2?6）
2．非间并情况 n
i
*
的表达式
因为对任一分布
!
!
i
i
N
t
n
=
∏
现在问题就是在以上两式限制条件下，如何选择 n
i
才能使 t 的数值最大。数学上就是求 t 有极大值时的 n
i
。 t 有极大值时 lnt 必有极大值（为了运算方便采用 lnt）
即对
!
!
i
i
N
t
n
=
∏
取对数,
lnt=lnN!－∑ln n
i
!
（ 1）
由斯特林公式,lnN!=NlnN－ N
所以 lnt= NlnN－ N－∑ n
i
lnn
i
+∑ n
i
= NlnN－∑ n
i
lnn
i
（ 2）
陕西师范大学精品课程 ……,物理化学,
第 10 页 共 40 页 2004-7-15
对
ln t
取极值，即作微分
ln
ln (ln 1) 0
iii
i
t
dt dn n dn
n
==?+=
∑
即∑(ln n
i
+1)dn
i
=0 （ 3）
同时上式应满足两个限制条件，即
∑ dn
i
=d N=0 （ 4）
∑ ε
i
dn
i
=dU =0 （ 5）
现在就是求既满足函数 lnt 有极值，又满足两个限制条件的变数值 n
1
n
2
…n
i
。在数字上采用拉格朗奇未定乘子法。此法就是求多元函数具有极值条件的方法。
即 在（ 4） （ 5）式分别乘以待定因子 α,β,再和（ 3）式相加得,
∑ (lnn
i
+1+α+βε
i
)dn
i
=0
1 并入 α内
∑ (lnn
i
+α+βε
i
)dn
i
=0 （ 6）
也就是最可几分布方式中的一套分布数 n
1
*
n
2
*
…n
i
*
必满足（ 6） 。
（注：最可几分布方式中各能级的分布数 n
i
都以 n
i
*
表示。 ）
因为 微变量 dn
i
≠ 0，则 lnn
i
+α+βε
i
=0
选择 α,β 值得
*
e
α?β
=
i
ε
i
n (3.2?7)
此式就是独立可别粒子体系的最可几分布的表达式。
即当 n
i
适合于（ 3.2?7）的那一种分配就是微观状态数最多的一种分配，这种分布就叫最可几分布，也叫玻兹曼分布。
同理对间并情况得到,
*
e
α?β
=
i
ε
ii
ng （ 3.2?8）
α,β 值的推导：因上式中含有两个待定因子 α和 β，需求出。
1.求 e
α
*
ee
β?α
= =
∑∑
i
ε
i
nN
e
e
α
β
∴ =
∑
i
ε
N
①
陕西师范大学精品课程 ……,物理化学,
第 11 页 共 40 页 2004-7-15
对间并情况：
*
ee
β?α
= =
∑∑
i
ε
ii
ngN
e
e
α
β
∴ =
∑
i
ε
i
N
g
②
①②分别代入 (3.2?7) （ 3.2?8）得,
*
e
e
β
β
=
∑
i
i
ε
i ε
N
n （3.2?9）
*
e
e
β
β
=
∑
i
i
ε
i
i ε
i
Ng
n
g
（ 3.2?10)
2.求 β 值
β 与体系的热力学温度 T 有关，即 β=1/k
B
T 则
B
B
/
*
/
e
e
=
∑
i
i
ε kT
i ε kT
N
n （ 3.2?11）
B
B
/
*
/
e
e
=
∑
i
i
ε kT
i
i ε kT
i
Ng
n
g
(3.2?12)
这就是定域子和经典离域子体系的最可几分布所遵从的规律,称为 Boltzmann 分布定律。
五．粒子的配分函数 q
1.分子配分函数的物理意义
分子配分函数的定义为
B
=
∑
i
ε kT
i
qge (3.2?13)
其中 g
i
为 i 能级的简并度，即 i 能级所有的量子状态数。由于体 系总能量的限制，并不是所有的能级及其量子状态都能被粒子所占据，玻兹曼因子
B
e
i
ε kT
小于 1 就是与 i 能级的能量 ε
i
有关的有效分数，因此,g
i
B
i
ε kT
e 是表示 i 能级的有效量子状态数，或称为有效状态数。 Σg
i
iB
ε kT
e 则表示所有能级的有效状态数之和，通常简称为,状态和,，这就是分子配分函数所表示的物理意义。
下面对代表的物理意义又作讨论,
由玻兹曼分布定律可知
BiB
B
*
ee
e
==
∑
i
i
ε kTε kT
ii i
ε kT
i
ng g
Nq g
（ 3.2?14）
该式表明分配在 i 能级上的粒子数 n
i
与体系总粒子数 N 之比等于该能级的有效状态数与陕西师范大学精品课程 ……,物理化学,
第 12 页 共 40 页 2004-7-15
所有能级的有效状态数总和之比。若考虑 i 能级和 j 能级，则由（ 3.9）式得到
B
B
/*
/*
e
e
=
i
j
ε kT
ii
ε kT
j
j
ng
n
g
（ 3.2?15）
这个关系式表明分配在两个不同的能级 i 和 j 上的粒子数之比等于两能级的有效状态数之比。以上两式也都称为玻兹曼公式。因此，人们常说体系中 N 个粒子是按玻兹曼公式分配于各个能级。
对非间并情况 g
i
=1,ε
0
表示最低能级，其粒子数为 n
0
*
则得,
0B
()/**
0
e
=
i
εε kT
i
nn （ 3.2?16）
若
0
ε =0，则
B
/**
0
e
=
i
ε kT
i
nn
六．最可几分布与平衡分布,
有一种分布所包含的微观状态数最多或出现的几率最大，称为最可几分布。
它有两个特点,
(1) 当 N 很大时，其它任何分布方式的微观状态数与最可几分布的比较可忽略不计，
即有一种分布所包含的微观状态数最多或出现的几率最大，称为最可几分布。
(2) 最可几分布代表平衡分布，代表体系平衡时的一切分布。
举例说明,N 个独立定域子在一个二维简并能级的两个量子态上的分布，设一个量子态上的粒子为 N? m 个，另一个为 m,则
m
!
()!!
=
N
t
Nmm-
()
m
00
!
!!
==
==
∑∑
NN
mm
N
t
Nmm
-
由二项式定理,
2
0
m1mm
!!
2!( 2)! !( )!
!
()!!
(1 ) 1
=
++ ++
=
+=+
∑
N
m
x
NN
xxNx
NmN
N
Nmm
+x
--
-
令 x=1 则
n
=0
!
2
()!!
= =
∑
N
m
N
Nmm-
在二项式中最大系数,
陕西师范大学精品课程 ……,物理化学,
第 13 页 共 40 页 2004-7-15
2
!
()!()!
22
=
N
N
t
NN
或
1
2
+
=
N
m 代替体系的一切分布
运用 Stirling 公式近似求得,
2
2
2
π
=
N
N
t
N
取对数
2
ln 2
2
ln +ln
π
=
N
N
t
N
（第二项可忽略）
2
ln ln≈
N
t?
当 N 足够大时,ln? 可用 lnt
max
来代替，即体系的总微观状态数的对数可由最可几分布所拥有的微观状态数的对数来代替。此即颉取最大项法，它使统计热力学的推导大为简化。
最可几分布的热力学几率为
1
2
max
(2/ ) 2
N
tNπ=?
设粒子数 N=10
24
，则最可几分布的几率为
11
13max
22
24
(,) ( ) ( ) 810
22 ππ10
tNN
p
N
== = ≈×
×
上述结果表明，最可几分布出现的几率只是一个很小的数，而且随着 N 的增大,P
（ N/2,N/2）的值反而减小。
当某一分布的分布数 M 与最可几分布 N/2 有一微小偏差时，这一分布的几率为：
!
!+!
222
()=
22
±
±=
N
N
NN
tm m m
N
Pm
-
(1)
应用斯特林公式 != 2
e
N
N
NN
π
后，最终得到
陕西师范大学精品课程 ……,物理化学,
第 14 页 共 40 页 2004-7-15
22
11
()= (2)
2
2
()()22
11+
22
,
22 2
2
(3)
()()
22
NN
m+m
NN
Pm
NN
m+m mm
NN
NN N
mm
N
NN
N
m+m
±?
π
±≈
=
-
-
-
由对数级数公式,
2
2
1
ln(1 )=
2
1
ln(1 )=
2
xx x
xxx
++?
-
--- -
取前两项
2
2
22
2
22
22
ln 1 1+
1
e(4)
22
11+
=
∴ =
NN
m+m
m
N
NN
m+m
mmm
NNN
mm
NN
-
-
(3),(4)代入 (2)中得,
22
12 2
()= e e
2
2
mm
NN
N
Pm
N
N
±=?
π
π
所以在此 m 区间范围内各个分布的总几率为,
2
2
2
+2
+2 +2 2
y
22
2
d=
2
21
( ) e d e dy 0.99993
2
=?
=? = =
π
π
∑
∫∫
m
N
N
N
N
mN
N
m
N
Pm m
N
--
m 偏离最可几分布数 N/2 如此小时,热力学几率接近 1,也就是说在
2
N
m- 到
2
N
+m
的狭小区间内，各个分布所具有的几率已非常接近体系的全部分布所具有的几率 1。因此最可几分布可以代表体系平衡时的一切分布。
七．经典统计的不可辨性修正,
经典统计认为粒子可辨，独立可辨粒子体系某一分布方式的微观状态数为,
陕西师范大学精品课程 ……,物理化学,
第 15 页 共 40 页 2004-7-15
!
!
=
∏
i
n
i
i
i
g
tN
n
若对独立不可辨粒子体系可对上式作一修正,因为可别粒子体系的 N 个可别粒子在不同能级上的排列方式数为 !N,而对不可别粒子体系,N 个粒子只有一种排列方式。
所以在 (1)中扣除 !N,将除以 !N 称为不可辨性修正。
故对独立不可辨粒子体系某一分布方式的微观状态数为,
!
!
!!
=
∏
∏
i
i
n
i
n
ii
i
g
N
ng
t=
Nn
这就是经典统计中不可别粒子体系的表达式。再用处理可别粒子体系类似的方法，
可得到玻尔兹曼分布定律公式,
B
B
/
*
/
e
e
=
∑
i
i
kT
i
i kT
i
Ng
n
g
ε
ε
（ 3.2?17）
即不可别与可别粒子体系的玻尔兹曼分布定律公式是一样的。 这种处理方法称为修正的玻尔兹曼统计，只是近似的。
八．粒子体系的量子统计
据每一量子状态所容纳粒子数的多少分,
玻色 -爱因斯坦统计：每个量子态容纳粒子数没限制。
费米 -狄拉克统计：每个量子态只容纳一个粒子。
玻色粒子是光子和由偶数个基本粒子（如质子、电子、中子）构成的粒子，不受泡里原理制约，具有整数自旋量子数的粒子。
费米粒子是基本粒子和由奇数个基本粒子构成的粒子，遵守泡里不相容原理，具有半整数自旋量子数的粒子。
1,玻色（ Bose） —爱因斯坦（ Einstein）统计
对热力学性质（ U,V,N）确定，自旋量子数为整数的粒子所组成的不可别独立粒子体系，体系中的每一个量子态能容纳的粒子数不限。因为量子态是可别的，可用排列组合的知识进行数学分析。
在第 i 个能级上（简并度为 g
i
） N
i
个粒子的分配样式，相当于 N
i
个相同的小球在 g
i
个箱中分配，每个箱子中能容纳的球数不限，其分配样式数即微观状态数 f
i
为
陕西师范大学精品课程 ……,物理化学,
第 16 页 共 40 页 2004-7-15
()
()!1!
!1
+
=
ii
ii
i
gN
gN
f （ 3.2?18）
体系中某分布 j 的微观状态数为各能级的微观状态数的乘积
()
()!1!
!1
11
+
==
ΠΠ
==
ii
ii
k
i
i
k
i
j
gN
gN
ft （ 3.2?19）
体系的总微观状态数? 仍为满足
∑
=
i
i
NN 和
ii
i
NUε =
∑
的各种分布微观状态数的加和
( )
()1
1!
!1!
i
i
k
ii
j
ij
NN ii
NU
Ng
t
Ng
ε
=
=
=
+?
==
∑
∑
∑∑
Π
i
（ 3.2?20）
在
∑
=
i
i
NN 和
ii
i
NUε =
∑
的限制条件下，用拉格朗日未定因子法得到最可几分布中，
第 I 能级上分配的粒子数
1
i
i
i
g
N
e
αβε
=
（ 3.2?21）
2,费米（ Fermi） —狄拉克（ Dirac）统计
对热力学性质（ U,V,N）确定，自旋量子数为半整数的粒子所组成的不可别独立粒子体系，体系中的一个量子态只能容纳一个粒子。因为量子态可别，可用排列组合的知识进行数学分析。
在能级 ε
i
上,N
i
个粒子在 g
i
个量子态中分配，每个量子态只能容纳一个粒子，且
N
i
< g
i
，其分配样式数 f
i
相当于从 g
i
个小球进行组合时的样式数。
()!!
!
iii
iN
gi
NgN
g
Cf
i
i
== （ 3.2?22）
体系中某种分布 j 的微观状态数 t
j
为各个能级微观状态数的乘积
()!!
!
11
iii
i
k
i
i
k
i
j
NgN
g
ft
Π=Π=
==
（ 3.2?23）
体系的总微观状态数
()
i
1
!
!!
=
=
=
== Π
∑
∑
∑∑
i
i
k
i
j
i
j NN
ii i
NU
g
t
NgN
ε
（ 3.2?24）
在
∑
=
i
i
NN 和
∑
=∈
i
ii
UN 的限制条件下，用拉格朗日未定因了法得到最概然分布中，
陕西师范大学精品课程 ……,物理化学,
第 17 页 共 40 页 2004-7-15
第 I 能级上分配的粒子数
1
i
i
i
g
N
e
αβε
=
+
（ 3.2?25）
式（ 3.2?18）至式（ 3.2?25）中各种符号所代表的意义与它们在经典统计中的意义相同。
同样，式（ 3.2?21）和式（ 3.2?25）中的
B
1
kT
β =?,即对于任何两个达成热平衡的独立粒子体系，它们的 β 值一定趋于一致。
九．三种统计的关系
三种统计的不同,
经典统计和量子统计的力学基础不同，且同属于量子统计的费米 —狄拉克统计和玻色 —爱因斯坦统计处理的对象不同，这些导致了三种统计关于某一分布方式微观状态数的表达式不同，最可几分布中某能级上分配的粒子数的表达式也不同。
1,某一分布的微观状态数表达式不同,
经典统计,
!
i
n
i
i
g
t
n
=
∏
玻色 -爱因斯坦统计,
(1)!
!( 1)!
ii
ii
ng
t
ng
+?
=
∏
费米 -狄拉克统计,
!
!( )!
i
ii i
g
t
ng n
=
∏
2,最可几分布定律的表达式不同,
经典统计,
*
i
i
i
ii
g
nge
e
αβε
αβε
+
==
玻色 -爱因斯坦统计：
*
1
i
i
i
g
n
e
αβε+
=
费米 -狄拉克统计,
*
1
i
i
i
g
n
e
αβε+
=
+
联系：当压力不太高、温度不太低时对等同粒子体系（气体分子体系）两种量子统计都可还原为玻兹曼统计。
故对气体这样的等同粒子体系可以用经典统计的理论加以描述，不会引起太多误差。但也有例外,
空腔辐射频率分布问题（光子气体）只遵守玻色 -爱因斯坦统计。
金属和半导体中的电子分布（电子气体）只遵守费米 -狄拉克统计。
陕西师范大学精品课程 ……,物理化学,
第 18 页 共 40 页 2004-7-15
第三节 配分函数与热力学函数的关系
统计热力学的主要任务之一就是要建立微观性质与宏观性质之间的关系。玻兹曼定理 S = k
B
ln?= k
B
ln t
max
的重要意义就在于此。然而在统计热力学中往往并不是直接通过具体计算最可几分布的微观状态数 t
max
来沟通微观和宏观，而是通过配分函数来建立二者之间的联系，配分函数虽然是一个无因次的量，但其值的大小却与体系的温度、体积、
分子的质量、核间距离、分子的平动、转动、振动以及电子和核的各种内部运动都有密切的联系，并且体系的各种热力学性质都可以表示为配分函数的函数 f(q)，因此配分函数是联系微观性质和宏观性质 的桥梁。统计热力学的首要 任务就是求算分子的配分函数。只要得到配分函数，就可以计算体系的热力学量。
一．分子配分函数与热力学性质的关系
1．热力学能
独立粒子体系的热力学能等于各粒子能量的总和，即 U = Σn
i
ε
i
。根据玻兹曼分布定律
B
=
i
kT
ii
N
nge
q
ε
,故
B
=?
∑
i
kT
ii
N
Uge
q
ε
ε
如果将 q 对 T 求偏微商，可得
B
2
,B
=?
∑
i
kT
i
i
VN
q
ge
TkT
ε
ε
即
B
2
B
,
=
∑
i
kT
ii
VN
q
ge k T
T
ε
ε
所以
22
BB
,,
ln
==
VN VN
Nq q
UkT NkT
qT T
（ 3.3?1）
显然这个关系对于定域子和离域子体系都是适用的。
若将式
0/
B
0
=
kT
qqe
ε
代入（ 3.3?1）式，则
2 0
B0
,
ln?
=+
VN
q
UNkT N
T
ε
令 U
0
= Nε
0
，它表示体系中的 N 个粒子全部处于基态时的总能量，即体系为 0 K 时的热力学能。于是上式可表示为
陕西师范大学精品课程 ……,物理化学,
第 19 页 共 40 页 2004-7-15
2 0
0B
,
ln?
=
VN
q
UU NkT
T
（ 3.3?2）
若规定 ε
0
= 0，则 U
0
= 0，也就是规定体系在 0 K 时的热力学能为零。显然，能量标度零点的选择不同，计算得到的体系热力学能值亦不同。
2．熵
不同体系的 t
max
不同，根据玻兹曼定理 S = k
B
ln t
max
，体系的熵值也不相同。对于定域子体系，由式
!
!
max
i
n
i
i
n
g
Nt
i
Π
=,所以
BBB
ln ! ln ln
!
==?
∑
Π
i
n
ii
i
i
gn
Sk N Nk Nk n
ng
将玻兹曼分布定律代入上式，可得
B
BB
,
(localized) ln
ln
ln
N
N
VN
U
Skq
T
q
kqNkT
T
=+
=+
（ 3.3?3）
而对离域子体系，
!
max
i
n
i
i
n
g
t
i
Π
=,所以
BB
ln ln
!
==?
∑
Π
i
n
ii
i
i
gn
Sk Nk k n
ng
将玻兹曼分布定律代入，即得
B
BB
,
(non - loc) ln
!
ln
ln
!
N
N
VN
qU
Sk
NT
qq
kNkT
NT
=+
=+
（ 3.3?4）
比较定域子体系和离域子体系的结果，由于在相同情况下二者的 t
max
相差 N！倍，故引起相应的熵值相差 k ln N! 。
如果将（ 3.3?1）式代入（ 3.3?3）或（ 3.3?4）式，可以证明，能量标度零点的不同选择对计算体系的熵值没有影响。
3．其它热力学量
根据 U,S 与 q的关系式，借助热力学函数之间的关系，容易得到其它热力学量的统计热力学表达式。
例如，体系的亥姆霍兹自由能 F，吉布斯自由能 G，热力学能 U，焓 H，定容热容陕西师范大学精品课程 ……,物理化学,
第 20 页 共 40 页 2004-7-15
C
V
等，在定域子体系的情况下，分别有以下关系式
BB
,
ln
() ln
=+
VN
q
Sloc kN q NkT
T
B
BB
,
2
BB
,,
2
B
.
B
,
() ln
ln
() ln
ln ln q
()
ln
()
ln
=?
=? +
=+
=
=
N
N
TN
VN TN
V
VN
V
TN
Floc kT q
q
Gloc kT q NkTV
V
q
Hloc NkT NkTV
TV
q
Cloc NkT
TT
q
pNkT
V
2
B
,
()
=
VN
q
Uloc NkT
T
类似地，也可以导出离域子体系的热力学量的统计表达式。无论是离域子体系还是定域子体系,U,H,C
V
的表达式是一样的，只是在函数 F,S,G 上相差一些常数项。
在通常计算热力学量的改变量? 时，这些常数项可以互相消去。
对可别粒子体系和不可别粒子体系函数 S,F,G 表示的比较,
BB
,
ln
ln
VN
q
SNkqNkT
T
=+
可别
BB
,
ln
ln
!
N
VN
qq
SNk NkT
NT
=+
不可别
B
lnFNkTq=?
可别
B
ln
!
N
q
FkT
N
=?
不可别
BB
,
ln
ln
TN
q
G NkT q NkTV
V
=? +
可别
BB
,
ln
ln
!
N
TN
qq
GNkT NkTV
NV
=+
不可别
理想气体有 TNkPV
B
=,则吉布斯自由能的定义可写成 TNkFG
B
+=,在此式的基础上，引用亥姆霍兹自由能统计热力学表达式得
陕西师范大学精品课程 ……,物理化学,
第 21 页 共 40 页 2004-7-15
()
B
1lnGNkT q=?
可别
B
ln
q
GNkT
N
=?
不可别
二．能量标度零点的选择
分子配分函数是分子各个能级的有效状态数之 和，而各个能级的有效状态数又与该能级的能量 ε
i
有关，欲确定能级能量必须首先选定标度零点。通常，分子能级能量标度零点有两种选择方案。一种是选取能量的绝对零点为起点，而确定基态能量为某一 ε
0
，
于是分子配分函数
B0B1B 2B
////
01
eeee
==+++
∑
"
i
kT kT kT kT
i
qg g g g
εεεε
另一种选择方案是规定基态的能量 ε
0
= 0，称为相对零点。这样得到的分子配分函数以
q
0
表示，即
B 1B 2B
/ //
001
eee
==+++
∑
"
i
kT kT kT
i
qg gg g
ε εε
式中?ε
i
= ε
i
-ε
0
，表示 i 能级相对于基态的能量值。能量标度零点的选择不同，求得的分子配分函数值也不同，显然二者的关系为
0B
/
0
e
=
kT
qq
ε
（ 3.3?5）
但对于玻兹曼分布定律
BB
0
ee
==
ii
kT kT
ii
NN
ng g
qq
εε
（ 3.3?6）
可见能量标度零点的选择不同，不会影响玻兹曼分布关系，而对某些热力学量的数值却会产生一定影响。具体情况看,分子配分函数与热力学性质的关系,部分。
三．对热力学三大定律本质的统计解释
1．热力学第一定律的本质
在组成不变的封闭体系中,若发生了一个微小的可逆变化,则根据热力学第一定律，
体系热力学能的变化为
dU = δQ + δW
由统计热力学原理可知，独立粒子体系的热力学能为 U = ∑n
i
ε
i
,当封闭体系经历了一个可逆变化后，热力学能的变化为
ddd
ii ii
Unnε ε=+
∑ ∑
(1)
陕西师范大学精品课程 ……,物理化学,
第 22 页 共 40 页 2004-7-15
上式右边的第一项表示能级固定时，由于能级分布数发生改变所引起的热力学能变化值，第二项则表示能级分布数固定时，由于能级改变所引起的热力学能增量。从经典力学原理可知，对于组成不变的封闭体系，热力学能的改变只能是体系与环境之间通过热和功的交换来体现。
2．热力学第二定律的本质
由熵的热力学定义式及式 (1)，得
1
R
ii
i
dQ
dS dn
TT
ε==
∑
(2)
上式就是热力学第二定律的表达式，它表明可逆过程的熵变与能级分布数的改变有关。而能级分布数的改变以为意味着体系的微观状态数发生了改变。
熵变是与体系微观状态数或热力学几率? 的变化相联系的。有公式,
S = k
B
ln?+ C (3)
式中 C 是积分常数。若?=1 时,S=0，则上式变成
S = k
B
ln?
此即为 Boltzmann 定理的数学表达式。由式可见，熵是体系微观状态数的一种量度。微观状态数? 较少的状态对应于较有序的状态,反之,? 值大的状态对应于较无序的状态。
因此，微观状态数? 的大小反映了体系有序程度的大小，亦即熵是体系有序程度或混乱程度的量度。当?=1 时，只有一个微观状态，体系最为有序，混乱程度为零，熵值为零。
基于以上讨论，我们可以作如下表述：在孤立体系中，自发变化的方向总是从较有序的状态向较无序的状态变化，即从微观状态数少的状态向微观状态数多的状态变化，从熵值小的状态向熵值大的状态变化，这就是热力学第二定律的本质。
3．热力学第三定律的本质
当 T→ 0 时，所有粒子都处于基态能级，此时?
0
=1，即把所有粒子放在一个能级上只有一个放法，体系只有一个微观状态，因此从玻兹曼定理，即式（ 6-25）可以得出结论：在 0 K 时物质的熵值为零，即
S
0
= k
B
ln?
0
= k
B
ln1 = 0 (4)
上式可以看作是热力学第三定律的统计表达式，这与热力学第三定律的表述,在 0 K 时任何纯物质的完美晶体的熵值为零,的结论是一致的。
第四节 配分函数的计算
一．平动、转动、振动的能级公式和能级的简并度
陕西师范大学精品课程 ……,物理化学,
第 23 页 共 40 页 2004-7-15
1,平动
质量为 m 的粒子，在体积为 a×b×c的方盒中运动，该三维平动子平动的能级公式为
()
2
222
y
xz
txyz 222
,,
8
=++
n
nnh
nnn
ma b c
ε （ 3.4?1）
其中 n
x
,n
y
,n
z
为 x,y,z 轴上的平动量子数，其取值为 1,2,3,…（正整数） 。
一组 n
x
,n
y
,n
z
的数值就决定一个平动量子态。如果运动空间是体积为 V 的立方盒，则该三维平动子平动的能级公式可写成
() ()
2
222
txyz x y z
23
,,
8
h
nnn n n n
mV
ε =++ （ 3.4?2）
三维平动子的基态能级是非简并的，其余以能级多数是简并的。如对在立方盒中运动的三维平动子，若某能级的能量
2
t
23
6
8
h
mV
ε = ×,则该能级有（ 2,1,1）,（ 1,2,1）,）
（ 1,1,2）三个量子态，即该能级的简并度为 3。
2、转动
线性刚性转子的能级公式为
() ( )
2
r
2
1
8
h
JJJ
I
ε
π
=+ （ 3.4?3）
其中 J 为转动量子数，其取值为 0,1,2,…； I 为转动惯量，对双原子分子
22
21
21
rr
mm
mm
I μ=
+
= （ 3.4?4）
其中 m
1
,m
2
为构成分子的两个原子的质量,μ 称为折合质量,r 为平衡时两原子的质心（核间）距离。
每一个由 J 决定的转动能级的简并度为
g
i,r
= 2J + 1 （ 3.4?5）
3、振动
单维简谐振子的能级公式为
()
v
11
22
hv hcvευ υ υ
=+ =+
（ 3.4?6）
其中 υ 为振动量子数，其取值为 0,1,2,…； v 为简谐振子的振动频率； C 为光陕西师范大学精品课程 ……,物理化学,
第 24 页 共 40 页 2004-7-15
速； v 为波数。
单维简谐振子的能级都是非简并的，即 g
i
,
v
=1 。三维谐振子的能量为
hv
zyx
+++
2
3
υυυ,能级的简并度为
g
i
,
v
=(S+1)(S+2)/2 （ 3.4?7）
其中
xyz
S υ υυ=++。
二．配分函数的析因性质
分子的运动可以分解为平动（ t）,转动（ r）,振动（ v）,电子运动（ e）和核运动（ n）
等。分子的各种运动形式可以看作是相对独立的，因此分子的能量可以认为是各种形式运动能量的加和，即
ε
i
= ε
i,t
+ ε
i,r
+ ε
i,v
+ ε
i,e
+ ε
i,n
各种运动形式分别具有其相应的简并度，分子总的简并度等于各个能级上简并度的乘积
g
i
= g
i,t
·g
i,r
·g
i,v
·g
i,e
·g
i,n
根据配分函数的定义，得
,t,r,v,e,n
trven
B
exp exp
++++
=?=
∑∑
i i i i i
i
i i,i,i,i,i,
ii
qg ggggg
kT k T
εεεεε
ε
数学上可以证明，几个独立变量乘积之和等于各自求和的乘积，于是
,t,r,v
,t,r,v
BB B
,e,n
,e,n
BB
trven
exp exp exp
exp e xp
=
∑∑∑
∑∑
i i i
i i i
iii
i i
i i
ii
qg g g
kT kT kT
gg
kT kT
qqqqq
εεε
εε
（ 3.4?8）
式中 q
t
,q
r
,q
v
,q
e
,q
n
分别称为平动配分函数、转动配分函数、振动配分函数、电子配分函数、核配分函数，分别代表分子的各种运动形式对分子配分函数的贡献。 （ 3.4?8）
式表明分子配分函数能够分解为各种形式运动贡献的乘积，这称为分子配分函数的析因子性质。应该指出，分子配分 函数的析因子性质只是对独 立粒子体系才是正确的。由
（ 3.4?8）式我们得到,
平动配分函数为
B
,t
t,t
=
∑
i
kT
i
i
qge
ε
（ 3.4?9）
转动配分函数为
陕西师范大学精品课程 ……,物理化学,
第 25 页 共 40 页 2004-7-15
B
,r
r,r
i
kT
i
i
qge
ε?
=
∑
（ 3.4?10）
振动配分函数为
vB
v,v
kT
i
i
qge
ε?
=
∑
（ 3.4?11）
电子配分函数为
B
,e
e,e
i
kT
i
i
qg
ε?
=
∑
（ 3.4?12）
核配分函数为
B
,n
n,n
i
kT
i
i
qge
ε?
=
∑
（ 3.4?13）
式（ 3.4?9）至式（ 3.4?13）中 g
i,t
,g
i,r
,g
i,v
,g
i,e
,g
i,n
和 ε
i,t
,ε
i,r
,ε
i,v
,ε
i,e
,ε
i,n
分别为平动、转动、振动、电子运动、核运动这五种运动形式能级的简并度和能级能量。
将（ 3.4?8 ）式应用于各热力学量的表达式，可得
2
B
,
222tvr
BBB
,,,
22e n
BB
,,
ln
ln lnln
ln ln
=
=++
++
VN
VNVN VN
VNVN
q
UNkT
T
qqq
Nk T Nk T Nk T
TTT
q q
Nk T Nk T
TT
(3.4?14)
=U
t
+U
r
+U
v
+U
e
+U
n
B
tvr
Bt Br Bv
e n
Be Bn
(loc) ln
ln ln ln
ln ln
N
NNN
NN
U
Skq
T
UUU
kq kq kq
TTT
U U
kq kq
TT
=+
=+++++
++++
(3.4?15)
=S
t
+S
r
+S
v
+S
e
+S
n
B
tt vr
BBB
e n
Be Bn
(non) ln
!
ln ln ln
!
ln ln
N
N
NN
NN
qU
Sk
NT
qU UU
kkqkq
NT T T
U U
kq kq
TT
=+
=+++++
++++
(3.4?16)
=S
t
+S
r
+S
v
+S
e
+S
n
使用类似的处理，其它热力学量均可以分解为 各种运动的贡献之和。应注意的是陕西师范大学精品课程 ……,物理化学,
第 26 页 共 40 页 2004-7-15
比较 （ 3.4?15） 和 （ 3.4?16） 式,可以看出定域子体系和离域子体系的平动熵 S
t
相差 k lnN!,
其余各项都是相同的。这是因为这两种粒子的区别就在于分子整体的平动，而分子的内部运动则是一样的。
三．各种运动形式的分子配分函数的计算及应用
根据分子配分函数的析因性质，只要求算出各 种运动形式的分子配分函数，就可以得到总配分函数，从而计算体系各热力学量。
1．平动配分函数
分子的平动，就是分子作为一个整体在容许的 空间范围内的运动，这相当于一个粒子了三维势箱中的运动，分子可以简化为三维平动子。根据量子力学原理，在长、宽、
高分别为 a,b,c 的三维势箱中，质量为 m 的三维平动子的能级公式为
2
22 2
y
x
txyz
222
(,,)
8
z
n
hn n
nnn
ma b c
ε
=++
式中 h 是普朗克常数,n
x
,n
y
,n
z
分别是 x,y,z 三个方向上的平动量子数，都只能取 1 到无穷大的任意正整数。当 n
x
,n
y
,n
z
一定时，就确定了一个平动量子态，量子态是非简并的。所以分子的平动配分函数为
xyz
xyz
2
222
y
xz
t
222
B
2
2 222
y
x z
BBB
xyz
exp
8
exp exp exp
888
=?+
=
=
∑∑∑
∑∑∑
nnn
nnn
n
nnh
q
mk T a b c
n
hhh
mk T a mk T b mk T c
qqq
式中 q
x
,q
y
,q
z
分别是三个方向上单维平动子的配分函数，各与一个平动自由度相对应。
如果令
2
2
2
B
8
=
h
mk Ta
α,则
xx
22
22x
2
11
B
exp exp( )
8
∞∞
==
==?
∑∑
nn
nh
qn
mk T a
α
一般来说,α
2
的值是很小的。当 α
2
<<1 时，一系列连续相差很小的数值求和，在数学上可以看作是连续的，因而可以用积分代替求和，即
x
22 22
xx x
0
1
exp( ) exp( )
2
n
qnnd
π
αα
α
∞
∞
=
=?=? =
∑
∫
陕西师范大学精品课程 ……,物理化学,
第 27 页 共 40 页 2004-7-15
对于 q
y
和 q
z
我们可作类似的处理。于是有
1/2 1/2 1/2
BBB
xyz
22
2
,,
πππ
=?=?=?
mk T mk T mk T
qaqbqc
hhh
所以
3/2 3/2
BB
t
22
=?=?
mk T mk T
q abc V
hh
ππ
（ 3.4?17）
其中 V= abc，是体系的体积。
将平动配分函数应用于理想气体（独立离域子体系）,有
22t
tB B B
,,
ln 3ln 3
22
==?=
VNVN
q T
UNkT NkT NkT
TT
（ 3.4?18）
t
,t B
,
3
2
==
v
VN
U
CNk
T
（ 3.4?19）
tt
tB B t B B
t
BB
3
ln ln ln !
!2
5
ln
2
=+=?+
=+
N
qU
Sk NkqkN Nk
NT
q
Nk Nk
N
（ 3.4?20）
（ 3.4?20）式称为沙克尔 (Sackur)-特鲁德 (Tetrode)方程，可直接用于计算理想气体的平动熵。平动对于其它热力学量的贡献亦可用类似的方法求得。
[例题 3.2] 计算 298.15 K,10
5
Pa 下,1 mol N
2
的平动配分函数。
解答,N
2
的摩尔质量为 14.008×2×10
－ 3
kg·mol
－ 1
，所以
326
23
14.008 2
10 kg 4.6515 10 kg
6.023 10
×
=×=×
×
m
h = 6.626×10
－ 34
J·s k
B
= 1.38×10
－ 23
J·K
-1
而 V = RT／ p = 0.02445 m
3
于是
3/2
29B
t 2
2
1.01 10
=?=×
mk T
qV
h
π
2．转动配分函数
对于双原子分子，除平动外，分子内部还有转 动和振动。一般说来，转动和振动互相有关，但为了简化，可近似认为二者相互独立，并将转动看作是刚性转子绕质心的运动。
根据量子力学处理的结果，线型刚性转子的能级公式为
2
r
2
() ( 1)
8
=+
π
h
JJJ
I
ε
陕西师范大学精品课程 ……,物理化学,
第 28 页 共 40 页 2004-7-15
式中 J 是转动量子数，可取 0,1,2,…等非负整数； I 是转动惯量，对于双原子分子
2
21
212
r
mm
mm
rI
+
== μ （ 3.4?21）
式中 m
1
,m
2
分别是两个原子的质量,μ 称为约化质量,r 是两个原子的质心距离。
分子的转动角动量的空间取向是量子化的，转动能级的简并度 g
r
= 2J+1，所以
2
r
2
0
B
(2 1) exp ( 1)
8
∞
=
=+?+
π
∑
J
h
qJ J
IkT
令 Θ
r
= h
2
/ 8π
2
I k
B
,称为分子转动的特征温度。通常温度下，对于大多数气体 Θ
r
<<T，
因此也可以用积分代替加和，即
r
r
0
r
2
B
2
(1)
(2 1)exp
8
∞
+
=+? =
π
=
∫
JJ T
qJ dJ
T
Ik T
h
Θ
Θ
（ 3.4?22）
（ 3.4－ 22）适用于异核双原子分子，也适用于非对称的线型多原子分子，如 HCN 等。
对于同核双原子分子或对称的线型多原子分子，由于分子对称性的原因
2
B
r
2
r
8π
==
IkTT
q
hσΘ σ
（ 3.4?23）
式中 σ 称为分子的对称数，即分子围绕对称由旋转 360°时具有相同位置的次数。对于同核双原子分子 σ =2，异核双原子分子 σ =1。
非线型多原子分子的情况比较复杂，可以导出
23/2 1/2
Bxyz
r
3
(8 ) ( )ππ
=
kT III
q
hσ
（ 3.4－ 24）
式中 I
x
,I
y
,I
z
分别是 x,y,z 三个轴向上的转动惯量。根据（ 3.4－ 22）,（ 3.4－ 23）,
（ 3.4－ 24）式计算分子的转动配分函数，其能量标度零点均选在转动的基态能级。
将分子的转动配分函数应用于双原子分子或线型多原子分子体系，有
22
r
rB B B
,
r
,r B
,
r
rBr Br B
ln ln
ln ln
===
==
=+= +
NV
v
VN
N
q T
UNkT NkT NkT
TT
U
CNk
T
U
Sk q Nk qNk
T
对于非线型多原子分子，有
陕西师范大学精品课程 ……,物理化学,
第 29 页 共 40 页 2004-7-15
rB
,r B
rBr B
3
2
3
2
3
ln
2
=
=
=+
v
UNkT
CNk
SNk q Nk
[例题 3.3] 已知 NO 的转动惯量 I = 16.4×10
－ 47
kg·m
2
，计算 298 K 时 NO 分子的转动配分函数
r
q 和该气体的摩尔转动热力学能 U
r,m
,摩尔转动熵 S
r,m
及摩尔定容热容 C
v,r,m
。
解答 NO 是异核双原子分子，对称数 σ = 1，所以
r
q = (8π
2
Ik
B
T) / h
2
= 8×3.14
2
×16.4×10
－ 47
×1.38×10
－ 23
×298 / (6.626×10
－ 34
)
2
= 121.2
U
r,m
= Lk
B
T = RT = 8.314 ×298 = 2478 (J·mol
－ 1
)
S
r,m
= Lk
B
ln
r
q + Lk
B
= Rln
r
q + R = 8.314 ln121.2 + 8.314 =48.2 (J·K
－ 1
·mol
－ 1
)
C
v,r,m
= Lk
B
= R = 8.314 (J·K
－ 1
·mol
－ 1
)
3．振动配分函数
分子的振动可以近似简化为若干个简谐振动,每个振动自由度相当于一个单维简谐振子。双原子分子只有一个振动自由度，可作为一维简谐振子处理。
根据量子力学的处理结果，单维简谐振子的能级公式为
v
1
()
2
hvευ υ
=+
式中 υ 为振动量子数，其取值为 0,1,2,…非负整数； v 为简谐振子的振动频率。当量子数
υ = 0 时 ε
v
（ 0） = 1/2 hv，称为零点振动能，亦即振动基态能级的能量。单维简谐振子的能级都是非简并的，于是
v
0
B
BBB
B
B
1
2
exp
2
exp 1 exp exp
2
1
exp
2
1exp
∞
=
+
=?
=+?+?+
=
∑
"
hv
q
kT
hv hv hv
kT kT kT
hv
kT
hv
kT
υ
υ
（ 3.4?25）
令 Θ
v
= hv/ k
B
,称为分子特征振动温度，则
陕西师范大学精品课程 ……,物理化学,
第 30 页 共 40 页 2004-7-15
v
v
v
1
exp
2
1exp
q
T
T
Θ
Θ
=
（ 3.4?26）
绝大多数气体的 Θ
v
都很高，通常温度下 Θ
v
>>T，因此上式可近似为
11
vB22
//
v
ee
==
ThkT
q
Θν
（ 3.4?27）
说明在通常温度下，气体分子几乎总是处于振动基态。只有当温度 T 接近 Θ
v
时，其它能级对配分函数才有实际贡献。当温度升高,T>>Θ
v
时，上述求和可用积分代替而得到
B
v
v
==
kTT
q
hvΘ
（ 3.4?28）
表明此时各振动能级对配分函数均有实际贡献。
若将能量零点选在振动基态能级，即令 ε
v
(0) = 0，则
v B
0v / /
11
()
1e 1e
==
T hkT
q
Θ ν
（ 3.4?29）
当 Θ
v
>>T 时,q(
0
)
v
= 1，于是可使计算大大简化。
对于由 n 个原子组成的多原子分子来说，分子具有多个振动自由度，总的振动配分函数应等于各单维简谐振子配分函数的乘积。
如果是线型多原子分子，总自由度为 3n，其中平动自由度为 3，转动自由度为 2，
振动自由度为 3n?5； 而对于非线型多原子分子,转动自由度为 3,振动自由度应为 3n?6。
于是振动配分函数可写为
B
0v /
1
()
1e
=
∏
i
f
hkT
i
q
ν
（ 3.4?30）
式中 f 为分子的振动自由度，其值为 3n?5 或 3n?6 。
即,
（ⅰ）线型分子
()
v
v
2
35
v
1
e
1e
=
′
=Π
T
n
T
i
q
Θ
Θ
（ 3.4?31）
或者
v
35
v
1
1
1e
=
=Π
n
T
i
q
Θ
（ 3.4?32）
（ⅱ）非线型分子
陕西师范大学精品课程 ……,物理化学,
第 31 页 共 40 页 2004-7-15
()
v
v
2
36
v
1
e
1e
=
′
=Π
T
n
T
i
q
Θ
Θ
（ 3.4?33）
或者
v
36
v
1
1
1e
=
=Π
n
T
i
q
Θ
（ 3.4?34）
式（ 3.4?31）和式（ 3.4?33）的能量标度零点与式（ 3.4－ 25）相同，式（ 3.4?32）和式
（ 3.4?34）的能量标度零点与式（ 3.4－ 29）相同。
[例题 3.4] 已知 NO 分子的特征振动波数 /ννc= =1907 cm
－ 1
，求 298 K 时该分子的振动配分函数和 NO 气体的摩尔转动能及摩尔振动熵。
解答,Θ
v
= hν/ k
B
= 6.626×10
－ 34
×(3×10
8
×1907×10
2
) / (1.38×10
－ 23
)=2747 K
因 Θ
v
(2747K) >> T(298K)
故 （ q
0
）
v
≈ 1
0v
m0,m B
,
ln( )
() 0
= ≈
NV
q
UU LkT
T
m0,m
v,m B 0 v
()
ln( ) 0
=+ ≈
UU
SLkq
T
4,电子配分函数
如果将能量零点选在基态，则电子的配分函数
,e 1,e
0e,e 0,e 1,e
BB
() exp exp
=?=+?+
∑
"
i
i
i
qg gg
kT kT
ε
ε
电子的能级间隔较大,?ε
i,e
一般几百 kJ·mok
－ 1
，如基态到第一激发态大约有 400 kJ·mok
－ 1
，粗略的估计，如果?ε
i,e
/ k
B
T > 5，则上式中第二项以后就可以忽略不计。除非在一个相当高的温度下，电子总是处在基态，而当温度增高时经常是在电子尚未激发之前分子就已经分解，因此上式第二项及其以后各项一般均可忽略，即
( q
0
)
e
= g
0,e
亦即电子配分函数等于电子运动基态的简并度。若电子运动的总角动量量子数为 j，对应 j 值有 (2j+1)个简并度。所以电子配分函数也可以写为
( q
0
)
e
= 2j + 1
应该注意的是，有少数原子如卤素原子和少数分子如 NO 等，它们的电子基态与第一激发态之间的能量间隔并不是很大，这时就需要考虑到第二项的贡献，即
陕西师范大学精品课程 ……,物理化学,
第 32 页 共 40 页 2004-7-15
1,e
0e 0,e 1,e
B
() exp
=+?
qgg
kT
ε
5,核配分函数
原子核运动的能级间隔更大，在通常的物理和化学变化中，原子核总是处在基态。
核配分函数的表达式为
B
,n
n,n
e
=
∑
i
kT
i
i
qg
ε
若将核运动基态能级的能量选作零，则上式变成
q
n
’ = g
0,n
核能级的简并度来源于原子核的自旋作用。若核自旋量子数为 s
n
，则原子或原子分子核运动基态能级的简并度为
g
0,n
= 2S
n
+ 1
对于多原子分子，其核运动基态能级的简并度是组成分子的各原子核运动基态能级的简并度之积，即
g
0,n
= ∏(2S
n
+ 1)
i
第五节 体系热力学性质的计算
上面我们对各种运动形式的配分函数求算作了讨论，有了 q 值就可具体计算体系的热力学函数值。这里主要讨论理想气体体系，即独立等同粒子体系。
一．分子的全配分函数
对于独立粒子体系而言，由于分子配分函数的析因子性质,分子的配分函数可以写为
trven
qqqqqq=
一般情况下，原子核的运动状态不会发生变化。而在某些化学反应中，电子转移会导致价电子运动状态发生变化，于是分子的全配分函数可以写为
trve
qqqqq=
通常在一般化学问题中，电子和核的运动状态都不会发生变化，电子配分函数和核配分函数没有必要进行计算，因此分子的全配分函数可以简化为
trv
qqqq=
综合以上对各种运动形式配分函数的讨论，在 不考虑电子和核运动的情况下，若采用基态能量为能量零点的规定，则各种不同分子的全配分函数可表示如下
陕西师范大学精品课程 ……,物理化学,
第 33 页 共 40 页 2004-7-15
(1) 单个原子
由于单个原子既无转动亦无振动，只需要考虑平动，故
3/2
B
0t
2
2π
==?
mk T
qq V
h
（ 3.5－ 1）
(2) 双原子分子
双原子分子有三个平动自由度，两个转动自由度和一个振动自由度，故有
B
3/2
2
BB
0 /
22
281
1e
ππ
=
hkT
mk T Ik T
qV
hh
ν
σ
（ 3.5－ 2）
由以上关系式知，只要知道了分子的质量 m、分子的转动惯量 I（ μr
2
）,对称数 σ
和振动频率 ν 等微观性质，就可以计算分子的全配分函数，从而求得宏观体系的热力学量。
二、体系热力学性质的计算
1．单原子分子
将单原子分子的配分函数表达式（ 3.5－ 1）代入不可别粒子体系的热力学性质的配分函数表达式中（纯气 体是不可别粒子体系）,可得到单原子分子气体体系相应的热力学性质。
(1) 亥姆霍兹自由能
e,0
e,0
t
net B n B eB
n,0
Bn,0 B
3/2
BB
3
n,0 e,0 B n.0 e,0
3/2
BBB
3
ln ln ln
!
[ exp( )] [ exp( )]
2
ln[( ) ] ln !
()l
2
ln[( ) ] ln
=++=
=
+
=+
+?
N
q
FFFF NkTq NkTqkT
N
Nk T g Nk T g
kT kT
mkT
Nk T V k T N
h
NN NkTgg
mkT
Nk T V Nk T N Nk T
h
ε
ε
π
εε
π
(3.5－ 3)
式中第一项是核和电子处于基态时的能量，第二项表示相应的间并度。在讨论热力学变量时，这些量是常量，都可以消去。
(2) 压力
B
Nk T
P
V
= (3.5－ 4)
(3) 摩尔热力学能 RTU
2
3
m
= (3.5－ 5)
(4) 摩尔焓 RTH
2
5
m
= (3.5－ 6)
陕西师范大学精品课程 ……,物理化学,
第 34 页 共 40 页 2004-7-15
(5) 摩尔恒容热容和摩尔恒压热容
RC
V
2
3
m,
= (3.5－ 7)
RC
P
2
3
m,
= (3.5－ 8)
(6) 摩尔熵
m
35
ln ln ln 9.685
22
r
p
SRMRTR
p
=+
\
(3.5－ 9)
(7) 摩尔吉布斯自由能
35
ln ln ln 30.47
22
+
r
p
GRT RM RTT T
p
m
=
\
(3.5－ 10)
式（ 3.5－ 10）称为沙克尔（ Sackur） ——特鲁德（ Tetrode）方程式。在应用式（ 3.5－ 9）
和式 （ 3.5－ 10） 计算惰性气体的 S
m
和 G
m
时,其中普适气体常数 R 必须取 8.314 J·K
－ 1
·mol
－ 1
，相对摩尔质量 M
r
的量纲为 g·mol
－ 1
，温度 T 的量纲为 K,压力 p 和标准压力 p
的量纲一定要相同。这样计算出来的 S
m
的单位为 J·K
－ 1
·mol
－ 1
,G
m
的单位为 J·mol
－ 1
。
(8) 化学势
对理想气体
B
=
Nk T
V
p
，代入 F 的表示式，得
n,0 e,0 B n,0 e,0
3/2
B
BBBB
3
()l
2
ln( ) ln ln
=+
π
+
kT g g
mk T
kT kT kT kT kT p
h
μ εε
对理想气体：
B
=
Nk T
V
p
对 1 摩尔气体而言，各项均乘以阿伏加德罗常数 L,Lk
B
=R,则 1 摩尔气体化学势为,
,0,0,0,0
2/3
B
B
2
()ln
2
ln( ) ln ln
=+
π
+
ne ne
LRTg
mk T
RTRkTRTp
h
μ εε
若气体的压力为标准压力 p
，则标准态的化学势 μ
可写为,
,
()
=
TV
F
N
μ
陕西师范大学精品课程 ……,物理化学,
第 35 页 共 40 页 2004-7-15
3/2
0,n 0,e 0 n 0 e
3
(2 )
( ) ln ln
ln ln
,,
mkT
LRTgRT
h
RT kT RT RT p
μεε
π
=+
+
\
\
对一定体系来说，上式右方都只是温度的函数，所以标准态的化学势可写作 μ
(T)，将上面两式相减，则得,
() ln( / )TRTppμμ=+
\ \
这就是气体化学势的表达式，这里也可以具体地看出 (T)中包含着哪些项。
2．双原子分子气体
将双原子分子的配分函数表达式（ 3.5－ 2）代入不可别粒子体系的各种热力学性质的配分函数表达式中，得到双原子分子气体体系相应的热力学性质。
(1) 热力学能
B
B
m
5
21
hv k T
hvkT
URT
e
=+
(2) 摩尔恒容热容
()
BB
B
2
2
,m
B
2
B
5
1
2
5
2
=+?
=+
hv k T hv k T
V
hv k T
hv
CRR ee
kT
hv
RR e
kT
显然，在室温时大多数双原子分子气体的 RCRC
PV
2
7
,
2
5
m,m,
== 。
(3) 熵尔熵
其中
t
S,
r
S,
v
S 为平动、转动、振动运动形式对摩尔熵的贡献。并且
()
()()
BB
32
B
B
t 3
2
B
r 2
1
v
B
2
5
ln
2
8
ln
ln 1 e e 1
π
=?+
π
=+
= +?
hv k T hv k T
mk T
kT
SR R
hP
Ik T
SR R
h
hv
SR R
kT
σ
(4) 吉布斯自由能
陕西师范大学精品课程 ……,物理化学,
第 36 页 共 40 页 2004-7-15
()
()
B
32
2
B
B
mm
2
8
ln 1 1
π
π
=
hv k T
mk T
Ik T
GRT V e
Lh hσ
其中 L 为阿伏加德罗（ Avogadro）常数。在计算熵时要注意，经典统计对粒子的不可别性只在平动运动形式的贡献项内进行修正。
[例题 3.5] 求氩（ Ar）气在其正常沸点 87.3 K 和标准压力时的摩尔热力学能
m
U,摩尔熵
m
S 及定圧摩尔热容
,mp
C 。
解答 氩为单个原子，可作为理想气体
3/2
20 3/ 2B
t 2
2
1.879 10 ( )
π
==?= ×
mk T
qq V MT V
h
Ar 的摩尔质量 M = 39.9 g·mol
－ 1
，标准压力 p
时的摩尔体积
m
V = RT / p
= 8.314×87.3 / 101325 = 7.16×10
－ 3
(m
3
)
故 q = 1.879×10
20
×(39.9×87.3)
3/2
×(7.16×10
－ 3
×10
6
) = 2.766×10
29
m
U = (3/2)RT = (3/2)×8.314×87.3 = 1089 (J·mol
－ 1
)
,mv
C = (3/2)R = (3/2)×8.314 = 12.47 (J·K
－ 1
·mol
－ 1
)
,mp
C = C
v,m
+ R = 20.87 (J·K
－ 1
·mol
－ 1
)
m
S = R ln( q /L) + (5/2)R = 8.314×[ln(2.766×10
29
／ 6.023×10
23
) +5/2]
= 129.2 (J·K
－ 1
·mol
－ 1
)
根据热力学第三定律，采用量热法得到的结果是 129.1J·K
－ 1
·mol
－ 1
。
三、理想气体状态方程的导出
将
B
ln
!
=?
N
q
FkT
N
代入 ()
T
F
p
V
=?
中，得
B
ln
()
=?
T
q
pNkT
V
又
trven
3/2
B
rven2
ln ln ln( )
2
ln ln ln( )
=+
π
=++
qq qq
mk T
Vqqq
h
式中，右端首尾两项与体系的体积 V 无关，故有,
BB
,
ln 1?
==?
TN
q
pNkT nkT
VV
陕西师范大学精品课程 ……,物理化学,
第 37 页 共 40 页 2004-7-15
而 Nk
B
=(N/k
B
)·Lk
B
=nR,L 为阿佛加德罗常数。所以
pV=nRT
四,统计熵的计算
在统计热力学中，对于非定位体系，有,
S = S
t
+ S
r
+ S
v
+ S
e
+ S
n
其中
平动熵
tt
tB B,
() ln
ln ( )
!
=+
N
VN
qq
Sk NkT
NT
转动熵
rr
rB B,
() ln
ln ( )
!
=+
N
VN
qq
Sk NkT
NT
振动熵
vv
vB B,
() ln
ln ( )
!
=+
N
VN
qq
Sk NkT
NT
电子熵
ee
eB B,
() ln
ln ( )
!
=+
N
VN
qq
Sk NkT
NT
核 熵
nn
nB B,
() ln
ln ( )
!
=+
N
VN
qq
Sk NkT
NT
式（ 6-110）表明，体系的熵值是物质各种运动形式对熵的贡献之和。
对于一般的物理及化学体系，可以忽略 S
n
，对少数物质需考虑 S
e
的贡献，所以统计熵主要来自 S
t
,S
r
和 S
v
的贡献。
1,平动熵
tt
tB B,
() ln
ln ( )
!
=+
N
VN
qq
Sk NkT
NT
应用平动配分函数，可将上式化简为,
3/2B
tB
2
25
[ln( ) ]
2
π
=+?
mk T V
SNk
hN
上式称为 Sackur-Tetrode 方程，表示平动运动对熵的贡献。
当理想气体的物质的量为 1 mol 时，有
3/2Bm
tB 2
3/2BB
B 2
25
[ln( ) ]
2
25
[ln( ) ]
2
π
=+?
π
=+?
mk T V
SNk
hL
mk T k T
Nk
hp
陕西师范大学精品课程 ……,物理化学,
第 38 页 共 40 页 2004-7-15
应用 SI 单位,k
B
=1.381×10
－23
J·K
-1
,h = 6.626×10
－34
J·s，若压力为标准压力 p
，则上式成为
S
m,t
= R(3/2 ln M + 5/2 lnT – 1.165)
式中 M 为摩尔质量，单位为 g·mol
－1
。
2,转动熵
对线型分子
当 Θr/T≤0.01 时，将
r
q 代入 S
r
的表达式，得
rB
r
(1 )=+
T
SNk
σΘ
3,振动熵
对双原子分子
将振动配分函数代入 S
v
的表达式，得
v
v
/v
vB/
/
[ln(1e)]
e1
=?
T
T
T
SNk
Θ
Θ
Θ
4,量热熵和光谱熵
在分子的各种运动形式中，平动、转动和振动等称为热运动形式，它们对体系熵的贡献称为热熵。体系的热能和热熵都能在体系升降温度时吞吐。在一般的温度范围内不能吞吐能量，而能使体系内粒子的位置或构型的无序运动形式所贡献的熵称为构型熵。
其中包括分子的电子运动（但对 NO
2
O
2
等可以是热熵）,核自旋运动及一系列甚至未被认识的不可穷尽的运动形式。构型熵的来源是作为晶体前身的气体的热熵，在气体凝聚成晶体的过程中，有那么一部分热熵变成了晶体的构型熵。在量热实验中，我们可以测定一个体系在降温过程中倾吐出来的热熵。包括平动熵、转动熵、振动熵及构型熵转化的热熵的总和，这样得出熵值称为体系的量热熵。就气体来说，分子的平动、转动和振动等热运动都会贡献出相应的热熵；对晶体物质来说，再加上源于热熵的构型熵，它们的总和称为气体或晶体的光谱熵。简单地说,
(1) 量热熵（或规定熵）
在热力学第三定律的基础上，根据热力学第二定律，利用量热测出的各种 H
p
和 △ H
数据而算得的 TK 时的熵叫量热熵。
陕西师范大学精品课程 ……,物理化学,
第 39 页 共 40 页 2004-7-15
mb
mb
0
0
0
d
d d d
=?=
=+++
∫
∑
∫∫∫
T
P
TT
TTT
PPP
TT
C
SSS T
T
CCC H
TTT
T
(2) 光谱熵（或统计熵）
根据统计热力学的方法，利用光谱实验数据先算出转动、振动、电子等运动形式的配分函数，从而得出的各种熵值，再加上平动熵而得到 TK 时的总熵值叫光谱熵。
在一般温度范围内吞吐的即是热熵,则我们在量热实验中测量的就只能是热熵了。 虽然分子内部的全部运动形式不可穷尽，但在一定温度范围内表现出来的热运动形式实际上已经明确。这也是我们计算统计熵的根据，即我们在实验中测量和在理论中计算以及使实际上已经明确。这也正是我们计算统计熵的根据，即我们在实验中测量和在理论中计算以及使实验和理论相互核对的正是热熵。而热熵与构型熵是可以相互转化的，当我们给一个热力学体系降温至复写纸的凝固点（
m
T ）时，体系会凝为固体，这是我们所能观测到的。在体系的内部实际上不存在着两个我们观测不到的温度：转化温度（
t
T ）和冻结温度（
f
T ） 。晶体中粒子从无序的位置或取向构型 可逆地转化为有序的位置或取向构型的温度称为转化温度；在某温度以下，由于动力学的原因（向有序的转化需较高的活化能）,体系中粒子的位置或构型将受到冻结，该温度称为冻结温度（一般
f
T ＜
m
T ＝。
对于转化温度高于冻结温度的晶体；可逆降温至温度达到转化温度和以下时，体系内部完成从无序向有序的可逆转化。这种情况下，来自气体热熵的构型熵又重新转化为能在降温过程中倾吐出来的热熵，即气体在降温过程中将它原有的全部热熵（包括源于气体熵的晶体的构型熵）作为量热熵倾吐出来。在这种情况下，光谱熵和量热熵在实验的误差范围内是一致的。对转化温度低于冻结温度的晶体，在降温过程中，由于动力学的原因，它的位置或取向构型的无序状态（该无序的位置或取向构型的总数? 可由统计的方法计算得到）被冻结，即被冻结为? 个无序状态中的任意一个状态。晶体没有机会将这种无序的位置或取向的构型熵（ S=k
B
ln?）重新以热熵的形式倾吐出来，来自它的前身
（气休） 热熵的这部分构型熵被永远保留在晶体中 （尽管被冻结的只是? 个状态中一个，
但被冻结的构型熵仍应以 S=k
B
ln?计算,即将被冻结的那个状态视为?个无序状态的,全权代表,） 。在这种情况下，我们由量热实验得到的量热熵只是体系全部热熵的一部分，
而通过计算得到的光谱熵更准确。这就是我们所说的统计热力学的一大成功之处：统计陕西师范大学精品课程 ……,物理化学,
第 40 页 共 40 页 2004-7-15
热力学关于低温熵的计算值比量热实验的测值更准确。简单地说,
一般情况下，量热熵和光谱熵能较好地吻合。个别物质在低温下，光谱熵大于量热熵，光谱熵更准确。
第 1 页 共 40 页 2004-7-15
第三章 统计热力学基础
第一节 绪论
一、统计热力学研究的对象、任务和方法
1,对象
研究的是大量微观粒子所构成的宏观物质处于平衡态时的性质。从研究对象来讲，
统计热力学和热力学一样都是研究宏观物质处于平衡态时的性质。热力学是根据从经验归纳得到的四条基本定律，而不管物质的微观结构和微观运动形态，因此只能得到联系各种宏观性质的一般规律，而不能给出微观性质与宏观性质之间的联系。而统计热力学则是从物质的微观结构出发来了解其宏观性质。两者研究的深度不同，所以有必要讨论后者。
2,任务
研究粒子所构成的体系的宏观行为，从粒子的微观性质来寻求体系的平均的宏观性质，这就是统计热力学的任务和研究内容。由此可见，统计热力学是从微观到宏观过渡的理论。它具有统计平均的性质，是联系物质的宏观性质与微观结构、沟通热力学与量子力学的一座桥梁。
3,方法
统计力学的研究方法是微观的方法，它根据统计单位的力学性质如速度、动量、位置、振动、转动等，用统计的方法来推求体系的热力学性质，例如压力、热容、熵等热力学函数。统计力学建立了体系的微观性质和宏观性质之间的联系。从这个意义上，统计力学又可称为统计热力学。
对于简单分子，应用统计热力学的方法进行处理，其结果是令人满意的。当然统计热力学也有自身的局限性，由于人们对于物质结构的认识不断深化，不断地修改充实物质结构的模型，同时模型本身也有近似性，所以由此得到的结论也有近似性。
从历史的发展来看，最早是由玻兹曼（ Boltzmann）以经典力学为基础建立的统计方法,称为经典统计热力学。 1900 年普朗克 （ Planck） 提出了量子论,麦克斯韦 （ Maxwell）
将能量量子化的概念引入统计热力学，对经典统计进行某些修正，发展成为麦克斯韦－
玻兹曼统计热力学方法。 1924 年量子力学建立后，在统计力学中不但所依赖的力学基础要改变，而且所用的统计方法也需要改变。由此产生了玻色－爱因斯坦（ Bose-Einstein）
统计和费米－狄拉克（ Fermi-Dirac）统计，分别适用于不同的体系。这两种统计方法都陕西师范大学精品课程 ……,物理化学,
第 2 页 共 40 页 2004-7-15
可以在一定的条件下通过适当的近似而得到玻兹曼统计。本章的内容就是简要介绍麦克斯韦－玻兹曼统计热力学的基本原理和应用。
二、统计体系的分类
在统计热力学中，按照构成体系的微观粒子（称为,统计单位,）的不同特性，可以将体系分为不同的类型。
按照粒子是否可以分辨，把体系分为定位体系（ localized system） （或称为定域子体系）和非定位体系（ non-localized system） （离域子体系）,前者的粒子可以彼此分辨，
而后者的粒子彼此不能分辨。例如气体分子处于无序运动之中，彼此无法区别，因此是离域子体系。而晶体，由于粒子是束缚在晶格位置上作振动运动，每个位置可以想象给予编号而加以区别，所以晶体是定域子体系。
按照统计单位之间有无相互作用,又可以把体系分为近独立粒子体系 (assembly of
independent particles)和非独立粒子体系（ assembly of interacting particles） 。前者或简称为独立粒子体系，其粒子之间的相互作用非常微弱，可以忽略不计，如理想气体，这种体系的总能量等于各个粒子的能量之和，即 =
i
i
U ε
∑；后者或称为相依粒子体系，其粒子之间其的相互作用不容忽略，如高圧下的实际气体等，这种体系的总能量除了各个粒子的能量之和外，还存在粒子之间相互作用的位能，即
I111
= (x,y,z,......x,y,z )
iNN
i
UUε +
∑
。 显然,粒子之间绝对无相互作用的体系是不存在的，
但可以把那些粒子之间的相互作用非常微弱可以忽略不计的体系，如低圧气体，作为独立粒子体系进行处理。本章中仅限于讨论独立粒子体系。
统计力学可分为两大阶段：经典统计力学和量子统计力学。前者是在 19 世纪末发展并成熟起来。在许多场合能给出满意的结果，但某些情况下它无法解释一些实验结果。
后者在二十世纪二十年代（ 1926 年）量子力学建立后发展起来的。它比经典统计力学能解释更广泛的宏观现象。本章着重讨论经典统计力学，只对量子统计力学稍加介绍。
三．数学知识
1．排列与组合
在统计热力学中，需要讨论粒子在不同能级上的分布，这在数学上相当于排列组合问 题。因此，先扼要介绍一些排列组合的有关知识。
(1)在 N 个不同的物体中，每次取出 m 个按照一定的顺序排成一列，称为从 N 个物体中每次取 m 个物体的排列；其排列的方式数为
陕西师范大学精品课程 ……,物理化学,
第 3 页 共 40 页 2004-7-15
A
N
m
= N!/(N - m)! (3.1?1)
当 m = N 时，上述排列称为全排列，全排列的方式总数为
A
N
N
= N!/(N – N)! = N ！ ~ (3.1?2)
其中规定 0！ = 1。
(2)若在 N 个物体中有 n
1
个相同，另外 n
2
个也彼此相同，其余的各不相同，则这 N
个物体的全排列方式数为
n!/N
1
!N
2
! (3.1?3)
(3) 将 N 个相同的物体放入 M 个相同的容器中（每个容器的容量不限），则放置的方式数为
（ N+M-1） !/N!(M-1)! (3.1?4)
(4) 将 N 个不同的物体放入 M 个容器中（每个容器容量不限），则放置的方式数为
M
N
(3.1?5)
(5) 在 N 个不同的物体中，每次提取 m 个，不管排列顺序编为一组，称为从 N 个不同物体中每次取出 m 个物体的组合，其组合数为
C
N
m
= N!/m!(N-m)! (3.1?6)
2,斯特林（ stirling）公式
在统计热力学中，常常要计算 N!。
阶乘 N！可展开如下式
23
1 1 139
!2 1
e 12 288 51840
N
N
NN
NN N
π
=++?
" (3.1?7)
(1) 当 N 是不太小的整数时，可近似为
1n ! 1n 2 1n 2 1n
e
N
N
NN NNNππ
==+?
(3.1?8)
误差约为
N12
1
。
(2) 当 N 是很大的整数时，上式进一步简化为
1n ! 1n 1n
e
N
N
NNN-
==
(3.1?9)
陕西师范大学精品课程 ……,物理化学,
第 4 页 共 40 页 2004-7-15
误差约为 ()Nπ2n1
2
1
。式（ 3.1?8）和式（ 3.1?9）称为斯特林公式。即,
当 N 很大时，有
lnN! = NlnN -N
当 N 很小时,有
!2()
e
=π
N
N
NN
N 越大，所得结果越精确。
第二节 粒子体系的统计分布
一,体系的微观状态及微观状态数,
2,体系的微观状态
用微观性质来描述的体系运动状态称体系的微观状态。
体系的微观状态是由各个粒子的微观状态所决定。粒子的微观状态是瞬息万变的，
由大量粒子组成的体系的微观运动状态也是千变万化的，如何描述粒子及体系的微观运动状态呢？经典力学与量子力学有不同的描述方法。
经典力学：粒子运动遵守牛顿运动方程，常用空间坐标 (q
x
,q
y
,q
z)
、瞬时速度或动量
(p
x
,p
y
,p
z)
来描述粒子的运动状态。在经典力学中，可根据粒子的空间坐标识别它们，故在经典力学中认为粒子是可别的。
量子力学：认为微观粒子是等同的、不可区别的。同时粒子具有波粒二相性，根据测不准原理,粒子不可能有确定的坐标和动量数值,所以不能用经典力学的方法来描述。
量子力学用本征函数 φ（波函数）和相应的本征值（能量值）来描述，微观粒子的运动服从薛定谔方程，即
22 2 2
22 2 2
()(,)
8
++ + =
h φφφ
Vxyz φE
mx y zπ
通过解此方程得到粒子的波函数 φ
i
相对应的能量 E
i
，具有一定 φ
i
,E
i
的状态称一种量子状态或量子态。从能级表达式得出能量是不连续的、量子化的。粒子处于不同的能级就表现不同的状态。
3,体系的微观状态数?
在体系的体积、总能量一定的情况下，含有 N 个粒子的体系中各种分布类型的样数之和，称体系的总微观状态数。 也称为体系的热力学几率。
陕西师范大学精品课程 ……,物理化学,
第 5 页 共 40 页 2004-7-15
设每种分布的微观状态数为 t
j
,则体系的总微观状态数就 等于各种分布的微观状态数之和，即?= Σt
j
。
我们知道,密闭体系的热力学能是熵和体积的特征函数,则体系的熵可表示为 S（ U、
V,N） 。由玻尔兹曼（ Boltzmann）公式（ S = k
B
ln?）知，该体系的总微观状态数? 也可以表示为?（ U,V,N） 。这就是说当体系的热 力学性质（ U,V,N）确定时，体系的能级和能级的简并度一定，体系的总微观状态数也一定，而且巨大数目的不同的微观状态对应于一个给定的宏观状态。
熵 S（体系的状态性质）与热力学几率?（体系的状态性质）之间存在某种函数关系，这个关系可表示为,
S = f(?)
玻兹曼定理为,
S = k
B
ln?
它揭示了体系的熵函数与其热力学 几率之间的关系。可以证明
k
B
=R/L=1.3806×10
-23
J·K
1
,即 k
B
是玻兹曼常数。 这是联系微观性质与宏观性质的关系式。
只要知道体系?，由此式可求出热力学函数 S，进而还可求出其它热力学函数。所以?
是重要的热力学函数。
二,统计热力学的基本假定
统计热力学认为：,对于宏观处于一定平衡状态的体系而言，任何一个可能出现的微观状态都具有相同的数学几率。,
统计热力学的这个基本假设，就是认为在所有 可能出现的微观状态中，任何一种状态都没有明显理由比其它微观状态出现的可能性更大些，这称为,等可几假设,。
上述假定的出发点是认为体系的热力学性质是所有可能出现的微观状态的统计平均。当我们对体系进行宏观测量时,需要一定的时间，在此时间内，体系将经历所有可能的微观状态。因此，宏观测得的某个物理量实际上是相应微观量的平均值，其中每个微观状态对平均值的贡献是相同的。这个假设的合理性已经由其引出的结论与实验事实相一致而得到证明。
由上述假定，对于拥有? 个微观状态的热力学体系，每一个微观状态出现的几率应为 1/?，而某一分布 x 出现的几率则为
P
x
= t
x
/?
式中 t
x
是该分布所拥有的微观状态数。此式表明，虽然各微观状态出现的几率相同，但陕西师范大学精品课程 ……,物理化学,
第 6 页 共 40 页 2004-7-15
各种分布出现的几率是不相同的。
为了说明统计热力学的基本假定，设某系统有 4 个可辨粒子 a,b,c,d，分配于两个相连的、容积相等的空间Ⅰ和Ⅱ之中，所有可能的分配形式如表 3.1 所列。
表 3.1 可分辨粒子在两个等容积空间的分配方式
分布方式 空间Ⅰ 空间Ⅱ 微观状态数 数学几率
（ 4,0）分布 a b c d 0 C
4
4
=1 1/16
（ 3,1）分布
a b c
a b d
a c d
b c d
d
c
b
a
C
4
3
=4 4/16
（ 2,2）分布
a b
a c
a d
b c
b d
c d
c d
b d
b c
a d
a c
a b
C
4
2
=6 6/16
（ 1,3）分布
a
b
c
d
b c d
a c d
a b d
a b c
C
4
1
=4 4/16
（ 0,4）分布 0 a b c d C
4
0
=1 1/16
在表 3.1 中列出的每一种可能的分配形式称为一个微观状态，所有可能的分配形式总数称为体系的总微观状态数，用? 表示。由表 3.1 可见，所研究体系的总微观状态数
=16。
根据等可几假设，上例中每个微观状态出现的数学几率都是 1/16。
三．经典统计 —麦克斯韦－玻兹曼统计
运用麦克斯韦－玻兹曼统计法所研究的体 系，应该具有两个特性,(1) 宏观状态确定的封闭体系； (2) 近独立粒子体系。
以下讨论的麦克斯韦－玻兹曼统计体系是按粒子的量子状态的分布来处理的，即能量是量子化的。这样做的目的是使推导过程更清晰和容易理解、接受。
陕西师范大学精品课程 ……,物理化学,
第 7 页 共 40 页 2004-7-15
1,非间并度、间并度
非间并度：一种能级只有一种量子状态即一种能量对应一种波函数。
n
1
n
2 ……….
n
i
ε
1
ε
2 ………,
ε
i
φ
1
φ
2 ………
φ
i
间并度：一种能级有多种量子状态即一种能量对应多个波函数。
n
1
n
2 ……………
n
i
ε
1
ε
2
………
,
ε
i
11 12 1g
...
i
φφ φ
21 22 2g
...
i
φφ φ ………
12
...
i
ii ig
φφ φ
注,g
i
是能级 ε
i
具有的量子状态数,称该能级的间并度或者统计权重。
2．麦克斯韦－玻兹曼统计,求体系的微观状态数?
(1) 对非间并情况,首先体系 U,N,V 固定,N 个粒子可区分
例,N 个可区分粒子的一种特殊分布,
有 n
1
个粒子的能量同为 ε
1
即 n
1
个粒子分配在 ε
1
能级上
有 n
2
个粒子的能量同为 ε
2
即 n
2
个粒子分配在 ε
2
能级上
┋ ┋
有 n
i
个粒子的能量同为 ε
i
即 n
i
个粒子分配在 ε
i
能级上
示意为,
能级
ε
1
ε
2
……,ε
i
一种分布方式 n
1
n
2 ………,
n
i
另一种分布方式 n
1
′
n
2
′
……….
n
i
′
由于粒子有弹性碰撞，在运动中相互交换能量，所以 N 个粒子在满足 U,N,V 一定的条件下可以有许多种不同的分布方式。
先设其中一分布方式的微观状态数，即花样数为 t
这个问题相当于将 N 个不同的物体分成若干堆，按排列组合公式
!
!
=
∏
i
i
i
N
t
n
(3.2?1)
那么对各种分布方式，由以上公式我们都可算出它的 t。
则体系的总微观状态数? 等于各个 t 的加和,即
陕西师范大学精品课程 ……,物理化学,
第 8 页 共 40 页 2004-7-15
!
!
==
∑∑
∏
i
ii
i
i
N
t
n
(3.2?2)
(2) 间并态的情况,在 U,N,V 一定的条件下,N 个可区分粒子在各能级上的一种分布
一种分布方式,
ε
1
ε
2
………ε
i
n
1
n
2 ……….
n
i
g
1
g
2
………g
i
即,能级是间并的，在 ε
i
能级上有 g
i
个间并度，即拥有 g
i
个不同的量子状态，而且假设每个量子状态容纳的粒子数不受限制，则 n
i
个处于能级 ε
i
的粒子中,每个粒子都可能有 g
i
个分布方式，故 n
i
个粒子在 g
i
个间并度内有
i
n
i
g 个微观状态数。
那么,ε
1
能级上有 n
1
个粒子，由于有 g
1
个间并度，产生
1
1
n
g 个微观状态
ε
2
能级上有 n
2
个粒子，由于有 g
2
个间并度，产生
2
2
n
g 个微观状态
┋ ┋ ┋
ε
i
能级上有 n
i
个粒子，由于有 g
i
个间并度，产生
i
n
i
g 个微观状态
所以一种分布方式，由于全部间并分布，所产生的微观状态数为
12
=
∏
ii
nnnn
i
gg g g
故一种分布方式所有的微观状态数 t 应为非间并时的微观状态数乘以间并分布时的微观状态数，即
i
!
!
!!
==
∏∏
∏
i
i
n
n
i
ii
i
gN
tgN
nn
(3.2?3)
而总的微观状态数? 为各种分布方式的 t 之和，即
!
!
==
∑∑
i
n
i
i
ii
i
g
tN
n
(3.2?4)
这样由 (3.2?2) (3.2?4)式就可求得非间并及间并态下的体系的总微观状态数。
注意：前提是在 U,N,V 一定的条件下。
四．最可几分布
对于指定状态的宏观体系，它的各种分布所拥有的微观状态数大小不一，其中必有一种分布所包含的微观状态数最多或出现的几率最大，称为最可几分布。
由表 3.1 可以看出，尽管各微观状态具有相同的数学几率，但各种分布的数学几率是不相同的。数学几率最大的分布称为,最可几分布,，上例中（ 2,2）分布就是最可几分陕西师范大学精品课程 ……,物理化学,
第 9 页 共 40 页 2004-7-15
布，其数学几率为 6/16。统计热力学认为，最可几分布可以代表体系的平衡分布。它可以代表体系处于热力学平衡状态时的一切分布状态。也就是说，对于一个粒子数目众多的实际平衡体系而言，其微观状态虽然千变万化，但基本上都是辗转于最可几分布以及与最可几分布没有实质差别的那些分布之中。因此，最可几分布是统计热力学最关注的分布。
最可几分布的微观状态数以 t
max
表示，则
≈ t
max
即 N 足够大时，在体系的各种分布方式中，求?只需考虑 t 中的最大项 t
m
对?的贡献就可以了，其它项可略去不计。实际上,N 很大时，体系的分布方式数目大的无法精确计算，也没必要计算。 （注：最可几分布方式中各能级的分布数 n
i
都以 n
i
*
表示。 ）
也就是只要求出与 t
m
相应的的那一套分布 n
1
*
n
2
*
…n
i
*
的微观状态数,就可代表体系的一切微观状态数。
S=k
B
lnt
max
=k
B
ln?
目的：求 t
m
时的那一套分布 n
1
*
n
2
*
…n
i
*
的表达式。求最可几分布是统计热力学的核心问题。
1,两个限制条件：首先 t
m
时的那一套分布 n
1
*
n
2
*
…n
i
*
必须满足 U,N 一定，以式子表示为,
∑ n
i
*
= N （ 3.2?5）
∑ n
i
*
ε=U （3.2?6）
2．非间并情况 n
i
*
的表达式
因为对任一分布
!
!
i
i
N
t
n
=
∏
现在问题就是在以上两式限制条件下，如何选择 n
i
才能使 t 的数值最大。数学上就是求 t 有极大值时的 n
i
。 t 有极大值时 lnt 必有极大值（为了运算方便采用 lnt）
即对
!
!
i
i
N
t
n
=
∏
取对数,
lnt=lnN!－∑ln n
i
!
（ 1）
由斯特林公式,lnN!=NlnN－ N
所以 lnt= NlnN－ N－∑ n
i
lnn
i
+∑ n
i
= NlnN－∑ n
i
lnn
i
（ 2）
陕西师范大学精品课程 ……,物理化学,
第 10 页 共 40 页 2004-7-15
对
ln t
取极值，即作微分
ln
ln (ln 1) 0
iii
i
t
dt dn n dn
n
==?+=
∑
即∑(ln n
i
+1)dn
i
=0 （ 3）
同时上式应满足两个限制条件，即
∑ dn
i
=d N=0 （ 4）
∑ ε
i
dn
i
=dU =0 （ 5）
现在就是求既满足函数 lnt 有极值，又满足两个限制条件的变数值 n
1
n
2
…n
i
。在数字上采用拉格朗奇未定乘子法。此法就是求多元函数具有极值条件的方法。
即 在（ 4） （ 5）式分别乘以待定因子 α,β,再和（ 3）式相加得,
∑ (lnn
i
+1+α+βε
i
)dn
i
=0
1 并入 α内
∑ (lnn
i
+α+βε
i
)dn
i
=0 （ 6）
也就是最可几分布方式中的一套分布数 n
1
*
n
2
*
…n
i
*
必满足（ 6） 。
（注：最可几分布方式中各能级的分布数 n
i
都以 n
i
*
表示。 ）
因为 微变量 dn
i
≠ 0，则 lnn
i
+α+βε
i
=0
选择 α,β 值得
*
e
α?β
=
i
ε
i
n (3.2?7)
此式就是独立可别粒子体系的最可几分布的表达式。
即当 n
i
适合于（ 3.2?7）的那一种分配就是微观状态数最多的一种分配，这种分布就叫最可几分布，也叫玻兹曼分布。
同理对间并情况得到,
*
e
α?β
=
i
ε
ii
ng （ 3.2?8）
α,β 值的推导：因上式中含有两个待定因子 α和 β，需求出。
1.求 e
α
*
ee
β?α
= =
∑∑
i
ε
i
nN
e
e
α
β
∴ =
∑
i
ε
N
①
陕西师范大学精品课程 ……,物理化学,
第 11 页 共 40 页 2004-7-15
对间并情况：
*
ee
β?α
= =
∑∑
i
ε
ii
ngN
e
e
α
β
∴ =
∑
i
ε
i
N
g
②
①②分别代入 (3.2?7) （ 3.2?8）得,
*
e
e
β
β
=
∑
i
i
ε
i ε
N
n （3.2?9）
*
e
e
β
β
=
∑
i
i
ε
i
i ε
i
Ng
n
g
（ 3.2?10)
2.求 β 值
β 与体系的热力学温度 T 有关，即 β=1/k
B
T 则
B
B
/
*
/
e
e
=
∑
i
i
ε kT
i ε kT
N
n （ 3.2?11）
B
B
/
*
/
e
e
=
∑
i
i
ε kT
i
i ε kT
i
Ng
n
g
(3.2?12)
这就是定域子和经典离域子体系的最可几分布所遵从的规律,称为 Boltzmann 分布定律。
五．粒子的配分函数 q
1.分子配分函数的物理意义
分子配分函数的定义为
B
=
∑
i
ε kT
i
qge (3.2?13)
其中 g
i
为 i 能级的简并度，即 i 能级所有的量子状态数。由于体 系总能量的限制，并不是所有的能级及其量子状态都能被粒子所占据，玻兹曼因子
B
e
i
ε kT
小于 1 就是与 i 能级的能量 ε
i
有关的有效分数，因此,g
i
B
i
ε kT
e 是表示 i 能级的有效量子状态数，或称为有效状态数。 Σg
i
iB
ε kT
e 则表示所有能级的有效状态数之和，通常简称为,状态和,，这就是分子配分函数所表示的物理意义。
下面对代表的物理意义又作讨论,
由玻兹曼分布定律可知
BiB
B
*
ee
e
==
∑
i
i
ε kTε kT
ii i
ε kT
i
ng g
Nq g
（ 3.2?14）
该式表明分配在 i 能级上的粒子数 n
i
与体系总粒子数 N 之比等于该能级的有效状态数与陕西师范大学精品课程 ……,物理化学,
第 12 页 共 40 页 2004-7-15
所有能级的有效状态数总和之比。若考虑 i 能级和 j 能级，则由（ 3.9）式得到
B
B
/*
/*
e
e
=
i
j
ε kT
ii
ε kT
j
j
ng
n
g
（ 3.2?15）
这个关系式表明分配在两个不同的能级 i 和 j 上的粒子数之比等于两能级的有效状态数之比。以上两式也都称为玻兹曼公式。因此，人们常说体系中 N 个粒子是按玻兹曼公式分配于各个能级。
对非间并情况 g
i
=1,ε
0
表示最低能级，其粒子数为 n
0
*
则得,
0B
()/**
0
e
=
i
εε kT
i
nn （ 3.2?16）
若
0
ε =0，则
B
/**
0
e
=
i
ε kT
i
nn
六．最可几分布与平衡分布,
有一种分布所包含的微观状态数最多或出现的几率最大，称为最可几分布。
它有两个特点,
(1) 当 N 很大时，其它任何分布方式的微观状态数与最可几分布的比较可忽略不计，
即有一种分布所包含的微观状态数最多或出现的几率最大，称为最可几分布。
(2) 最可几分布代表平衡分布，代表体系平衡时的一切分布。
举例说明,N 个独立定域子在一个二维简并能级的两个量子态上的分布，设一个量子态上的粒子为 N? m 个，另一个为 m,则
m
!
()!!
=
N
t
Nmm-
()
m
00
!
!!
==
==
∑∑
NN
mm
N
t
Nmm
-
由二项式定理,
2
0
m1mm
!!
2!( 2)! !( )!
!
()!!
(1 ) 1
=
++ ++
=
+=+
∑
N
m
x
NN
xxNx
NmN
N
Nmm
+x
--
-
令 x=1 则
n
=0
!
2
()!!
= =
∑
N
m
N
Nmm-
在二项式中最大系数,
陕西师范大学精品课程 ……,物理化学,
第 13 页 共 40 页 2004-7-15
2
!
()!()!
22
=
N
N
t
NN
或
1
2
+
=
N
m 代替体系的一切分布
运用 Stirling 公式近似求得,
2
2
2
π
=
N
N
t
N
取对数
2
ln 2
2
ln +ln
π
=
N
N
t
N
（第二项可忽略）
2
ln ln≈
N
t?
当 N 足够大时,ln? 可用 lnt
max
来代替，即体系的总微观状态数的对数可由最可几分布所拥有的微观状态数的对数来代替。此即颉取最大项法，它使统计热力学的推导大为简化。
最可几分布的热力学几率为
1
2
max
(2/ ) 2
N
tNπ=?
设粒子数 N=10
24
，则最可几分布的几率为
11
13max
22
24
(,) ( ) ( ) 810
22 ππ10
tNN
p
N
== = ≈×
×
上述结果表明，最可几分布出现的几率只是一个很小的数，而且随着 N 的增大,P
（ N/2,N/2）的值反而减小。
当某一分布的分布数 M 与最可几分布 N/2 有一微小偏差时，这一分布的几率为：
!
!+!
222
()=
22
±
±=
N
N
NN
tm m m
N
Pm
-
(1)
应用斯特林公式 != 2
e
N
N
NN
π
后，最终得到
陕西师范大学精品课程 ……,物理化学,
第 14 页 共 40 页 2004-7-15
22
11
()= (2)
2
2
()()22
11+
22
,
22 2
2
(3)
()()
22
NN
m+m
NN
Pm
NN
m+m mm
NN
NN N
mm
N
NN
N
m+m
±?
π
±≈
=
-
-
-
由对数级数公式,
2
2
1
ln(1 )=
2
1
ln(1 )=
2
xx x
xxx
++?
-
--- -
取前两项
2
2
22
2
22
22
ln 1 1+
1
e(4)
22
11+
=
∴ =
NN
m+m
m
N
NN
m+m
mmm
NNN
mm
NN
-
-
(3),(4)代入 (2)中得,
22
12 2
()= e e
2
2
mm
NN
N
Pm
N
N
±=?
π
π
所以在此 m 区间范围内各个分布的总几率为,
2
2
2
+2
+2 +2 2
y
22
2
d=
2
21
( ) e d e dy 0.99993
2
=?
=? = =
π
π
∑
∫∫
m
N
N
N
N
mN
N
m
N
Pm m
N
--
m 偏离最可几分布数 N/2 如此小时,热力学几率接近 1,也就是说在
2
N
m- 到
2
N
+m
的狭小区间内，各个分布所具有的几率已非常接近体系的全部分布所具有的几率 1。因此最可几分布可以代表体系平衡时的一切分布。
七．经典统计的不可辨性修正,
经典统计认为粒子可辨，独立可辨粒子体系某一分布方式的微观状态数为,
陕西师范大学精品课程 ……,物理化学,
第 15 页 共 40 页 2004-7-15
!
!
=
∏
i
n
i
i
i
g
tN
n
若对独立不可辨粒子体系可对上式作一修正,因为可别粒子体系的 N 个可别粒子在不同能级上的排列方式数为 !N,而对不可别粒子体系,N 个粒子只有一种排列方式。
所以在 (1)中扣除 !N,将除以 !N 称为不可辨性修正。
故对独立不可辨粒子体系某一分布方式的微观状态数为,
!
!
!!
=
∏
∏
i
i
n
i
n
ii
i
g
N
ng
t=
Nn
这就是经典统计中不可别粒子体系的表达式。再用处理可别粒子体系类似的方法，
可得到玻尔兹曼分布定律公式,
B
B
/
*
/
e
e
=
∑
i
i
kT
i
i kT
i
Ng
n
g
ε
ε
（ 3.2?17）
即不可别与可别粒子体系的玻尔兹曼分布定律公式是一样的。 这种处理方法称为修正的玻尔兹曼统计，只是近似的。
八．粒子体系的量子统计
据每一量子状态所容纳粒子数的多少分,
玻色 -爱因斯坦统计：每个量子态容纳粒子数没限制。
费米 -狄拉克统计：每个量子态只容纳一个粒子。
玻色粒子是光子和由偶数个基本粒子（如质子、电子、中子）构成的粒子，不受泡里原理制约，具有整数自旋量子数的粒子。
费米粒子是基本粒子和由奇数个基本粒子构成的粒子，遵守泡里不相容原理，具有半整数自旋量子数的粒子。
1,玻色（ Bose） —爱因斯坦（ Einstein）统计
对热力学性质（ U,V,N）确定，自旋量子数为整数的粒子所组成的不可别独立粒子体系，体系中的每一个量子态能容纳的粒子数不限。因为量子态是可别的，可用排列组合的知识进行数学分析。
在第 i 个能级上（简并度为 g
i
） N
i
个粒子的分配样式，相当于 N
i
个相同的小球在 g
i
个箱中分配，每个箱子中能容纳的球数不限，其分配样式数即微观状态数 f
i
为
陕西师范大学精品课程 ……,物理化学,
第 16 页 共 40 页 2004-7-15
()
()!1!
!1
+
=
ii
ii
i
gN
gN
f （ 3.2?18）
体系中某分布 j 的微观状态数为各能级的微观状态数的乘积
()
()!1!
!1
11
+
==
ΠΠ
==
ii
ii
k
i
i
k
i
j
gN
gN
ft （ 3.2?19）
体系的总微观状态数? 仍为满足
∑
=
i
i
NN 和
ii
i
NUε =
∑
的各种分布微观状态数的加和
( )
()1
1!
!1!
i
i
k
ii
j
ij
NN ii
NU
Ng
t
Ng
ε
=
=
=
+?
==
∑
∑
∑∑
Π
i
（ 3.2?20）
在
∑
=
i
i
NN 和
ii
i
NUε =
∑
的限制条件下，用拉格朗日未定因子法得到最可几分布中，
第 I 能级上分配的粒子数
1
i
i
i
g
N
e
αβε
=
（ 3.2?21）
2,费米（ Fermi） —狄拉克（ Dirac）统计
对热力学性质（ U,V,N）确定，自旋量子数为半整数的粒子所组成的不可别独立粒子体系，体系中的一个量子态只能容纳一个粒子。因为量子态可别，可用排列组合的知识进行数学分析。
在能级 ε
i
上,N
i
个粒子在 g
i
个量子态中分配，每个量子态只能容纳一个粒子，且
N
i
< g
i
，其分配样式数 f
i
相当于从 g
i
个小球进行组合时的样式数。
()!!
!
iii
iN
gi
NgN
g
Cf
i
i
== （ 3.2?22）
体系中某种分布 j 的微观状态数 t
j
为各个能级微观状态数的乘积
()!!
!
11
iii
i
k
i
i
k
i
j
NgN
g
ft
Π=Π=
==
（ 3.2?23）
体系的总微观状态数
()
i
1
!
!!
=
=
=
== Π
∑
∑
∑∑
i
i
k
i
j
i
j NN
ii i
NU
g
t
NgN
ε
（ 3.2?24）
在
∑
=
i
i
NN 和
∑
=∈
i
ii
UN 的限制条件下，用拉格朗日未定因了法得到最概然分布中，
陕西师范大学精品课程 ……,物理化学,
第 17 页 共 40 页 2004-7-15
第 I 能级上分配的粒子数
1
i
i
i
g
N
e
αβε
=
+
（ 3.2?25）
式（ 3.2?18）至式（ 3.2?25）中各种符号所代表的意义与它们在经典统计中的意义相同。
同样，式（ 3.2?21）和式（ 3.2?25）中的
B
1
kT
β =?,即对于任何两个达成热平衡的独立粒子体系，它们的 β 值一定趋于一致。
九．三种统计的关系
三种统计的不同,
经典统计和量子统计的力学基础不同，且同属于量子统计的费米 —狄拉克统计和玻色 —爱因斯坦统计处理的对象不同，这些导致了三种统计关于某一分布方式微观状态数的表达式不同，最可几分布中某能级上分配的粒子数的表达式也不同。
1,某一分布的微观状态数表达式不同,
经典统计,
!
i
n
i
i
g
t
n
=
∏
玻色 -爱因斯坦统计,
(1)!
!( 1)!
ii
ii
ng
t
ng
+?
=
∏
费米 -狄拉克统计,
!
!( )!
i
ii i
g
t
ng n
=
∏
2,最可几分布定律的表达式不同,
经典统计,
*
i
i
i
ii
g
nge
e
αβε
αβε
+
==
玻色 -爱因斯坦统计：
*
1
i
i
i
g
n
e
αβε+
=
费米 -狄拉克统计,
*
1
i
i
i
g
n
e
αβε+
=
+
联系：当压力不太高、温度不太低时对等同粒子体系（气体分子体系）两种量子统计都可还原为玻兹曼统计。
故对气体这样的等同粒子体系可以用经典统计的理论加以描述，不会引起太多误差。但也有例外,
空腔辐射频率分布问题（光子气体）只遵守玻色 -爱因斯坦统计。
金属和半导体中的电子分布（电子气体）只遵守费米 -狄拉克统计。
陕西师范大学精品课程 ……,物理化学,
第 18 页 共 40 页 2004-7-15
第三节 配分函数与热力学函数的关系
统计热力学的主要任务之一就是要建立微观性质与宏观性质之间的关系。玻兹曼定理 S = k
B
ln?= k
B
ln t
max
的重要意义就在于此。然而在统计热力学中往往并不是直接通过具体计算最可几分布的微观状态数 t
max
来沟通微观和宏观，而是通过配分函数来建立二者之间的联系，配分函数虽然是一个无因次的量，但其值的大小却与体系的温度、体积、
分子的质量、核间距离、分子的平动、转动、振动以及电子和核的各种内部运动都有密切的联系，并且体系的各种热力学性质都可以表示为配分函数的函数 f(q)，因此配分函数是联系微观性质和宏观性质 的桥梁。统计热力学的首要 任务就是求算分子的配分函数。只要得到配分函数，就可以计算体系的热力学量。
一．分子配分函数与热力学性质的关系
1．热力学能
独立粒子体系的热力学能等于各粒子能量的总和，即 U = Σn
i
ε
i
。根据玻兹曼分布定律
B
=
i
kT
ii
N
nge
q
ε
,故
B
=?
∑
i
kT
ii
N
Uge
q
ε
ε
如果将 q 对 T 求偏微商，可得
B
2
,B
=?
∑
i
kT
i
i
VN
q
ge
TkT
ε
ε
即
B
2
B
,
=
∑
i
kT
ii
VN
q
ge k T
T
ε
ε
所以
22
BB
,,
ln
==
VN VN
Nq q
UkT NkT
qT T
（ 3.3?1）
显然这个关系对于定域子和离域子体系都是适用的。
若将式
0/
B
0
=
kT
qqe
ε
代入（ 3.3?1）式，则
2 0
B0
,
ln?
=+
VN
q
UNkT N
T
ε
令 U
0
= Nε
0
，它表示体系中的 N 个粒子全部处于基态时的总能量，即体系为 0 K 时的热力学能。于是上式可表示为
陕西师范大学精品课程 ……,物理化学,
第 19 页 共 40 页 2004-7-15
2 0
0B
,
ln?
=
VN
q
UU NkT
T
（ 3.3?2）
若规定 ε
0
= 0，则 U
0
= 0，也就是规定体系在 0 K 时的热力学能为零。显然，能量标度零点的选择不同，计算得到的体系热力学能值亦不同。
2．熵
不同体系的 t
max
不同，根据玻兹曼定理 S = k
B
ln t
max
，体系的熵值也不相同。对于定域子体系，由式
!
!
max
i
n
i
i
n
g
Nt
i
Π
=,所以
BBB
ln ! ln ln
!
==?
∑
Π
i
n
ii
i
i
gn
Sk N Nk Nk n
ng
将玻兹曼分布定律代入上式，可得
B
BB
,
(localized) ln
ln
ln
N
N
VN
U
Skq
T
q
kqNkT
T
=+
=+
（ 3.3?3）
而对离域子体系，
!
max
i
n
i
i
n
g
t
i
Π
=,所以
BB
ln ln
!
==?
∑
Π
i
n
ii
i
i
gn
Sk Nk k n
ng
将玻兹曼分布定律代入，即得
B
BB
,
(non - loc) ln
!
ln
ln
!
N
N
VN
qU
Sk
NT
kNkT
NT
=+
=+
（ 3.3?4）
比较定域子体系和离域子体系的结果，由于在相同情况下二者的 t
max
相差 N！倍，故引起相应的熵值相差 k ln N! 。
如果将（ 3.3?1）式代入（ 3.3?3）或（ 3.3?4）式，可以证明，能量标度零点的不同选择对计算体系的熵值没有影响。
3．其它热力学量
根据 U,S 与 q的关系式，借助热力学函数之间的关系，容易得到其它热力学量的统计热力学表达式。
例如，体系的亥姆霍兹自由能 F，吉布斯自由能 G，热力学能 U，焓 H，定容热容陕西师范大学精品课程 ……,物理化学,
第 20 页 共 40 页 2004-7-15
C
V
等，在定域子体系的情况下，分别有以下关系式
BB
,
ln
() ln
=+
VN
q
Sloc kN q NkT
T
B
BB
,
2
BB
,,
2
B
.
B
,
() ln
ln
() ln
ln ln q
()
ln
()
ln
=?
=? +
=+
=
=
N
N
TN
VN TN
V
VN
V
TN
Floc kT q
q
Gloc kT q NkTV
V
q
Hloc NkT NkTV
TV
q
Cloc NkT
TT
q
pNkT
V
2
B
,
()
=
VN
q
Uloc NkT
T
类似地，也可以导出离域子体系的热力学量的统计表达式。无论是离域子体系还是定域子体系,U,H,C
V
的表达式是一样的，只是在函数 F,S,G 上相差一些常数项。
在通常计算热力学量的改变量? 时，这些常数项可以互相消去。
对可别粒子体系和不可别粒子体系函数 S,F,G 表示的比较,
BB
,
ln
ln
VN
q
SNkqNkT
T
=+
可别
BB
,
ln
ln
!
N
VN
SNk NkT
NT
=+
不可别
B
lnFNkTq=?
可别
B
ln
!
N
q
FkT
N
=?
不可别
BB
,
ln
ln
TN
q
G NkT q NkTV
V
=? +
可别
BB
,
ln
ln
!
N
TN
GNkT NkTV
NV
=+
不可别
理想气体有 TNkPV
B
=,则吉布斯自由能的定义可写成 TNkFG
B
+=,在此式的基础上，引用亥姆霍兹自由能统计热力学表达式得
陕西师范大学精品课程 ……,物理化学,
第 21 页 共 40 页 2004-7-15
()
B
1lnGNkT q=?
可别
B
ln
q
GNkT
N
=?
不可别
二．能量标度零点的选择
分子配分函数是分子各个能级的有效状态数之 和，而各个能级的有效状态数又与该能级的能量 ε
i
有关，欲确定能级能量必须首先选定标度零点。通常，分子能级能量标度零点有两种选择方案。一种是选取能量的绝对零点为起点，而确定基态能量为某一 ε
0
，
于是分子配分函数
B0B1B 2B
////
01
eeee
==+++
∑
"
i
kT kT kT kT
i
qg g g g
εεεε
另一种选择方案是规定基态的能量 ε
0
= 0，称为相对零点。这样得到的分子配分函数以
q
0
表示，即
B 1B 2B
/ //
001
eee
==+++
∑
"
i
kT kT kT
i
qg gg g
ε εε
式中?ε
i
= ε
i
-ε
0
，表示 i 能级相对于基态的能量值。能量标度零点的选择不同，求得的分子配分函数值也不同，显然二者的关系为
0B
/
0
e
=
kT
ε
（ 3.3?5）
但对于玻兹曼分布定律
BB
0
ee
==
ii
kT kT
ii
NN
ng g
εε
（ 3.3?6）
可见能量标度零点的选择不同，不会影响玻兹曼分布关系，而对某些热力学量的数值却会产生一定影响。具体情况看,分子配分函数与热力学性质的关系,部分。
三．对热力学三大定律本质的统计解释
1．热力学第一定律的本质
在组成不变的封闭体系中,若发生了一个微小的可逆变化,则根据热力学第一定律，
体系热力学能的变化为
dU = δQ + δW
由统计热力学原理可知，独立粒子体系的热力学能为 U = ∑n
i
ε
i
,当封闭体系经历了一个可逆变化后，热力学能的变化为
ddd
ii ii
Unnε ε=+
∑ ∑
(1)
陕西师范大学精品课程 ……,物理化学,
第 22 页 共 40 页 2004-7-15
上式右边的第一项表示能级固定时，由于能级分布数发生改变所引起的热力学能变化值，第二项则表示能级分布数固定时，由于能级改变所引起的热力学能增量。从经典力学原理可知，对于组成不变的封闭体系，热力学能的改变只能是体系与环境之间通过热和功的交换来体现。
2．热力学第二定律的本质
由熵的热力学定义式及式 (1)，得
1
R
ii
i
dQ
dS dn
TT
ε==
∑
(2)
上式就是热力学第二定律的表达式，它表明可逆过程的熵变与能级分布数的改变有关。而能级分布数的改变以为意味着体系的微观状态数发生了改变。
熵变是与体系微观状态数或热力学几率? 的变化相联系的。有公式,
S = k
B
ln?+ C (3)
式中 C 是积分常数。若?=1 时,S=0，则上式变成
S = k
B
ln?
此即为 Boltzmann 定理的数学表达式。由式可见，熵是体系微观状态数的一种量度。微观状态数? 较少的状态对应于较有序的状态,反之,? 值大的状态对应于较无序的状态。
因此，微观状态数? 的大小反映了体系有序程度的大小，亦即熵是体系有序程度或混乱程度的量度。当?=1 时，只有一个微观状态，体系最为有序，混乱程度为零，熵值为零。
基于以上讨论，我们可以作如下表述：在孤立体系中，自发变化的方向总是从较有序的状态向较无序的状态变化，即从微观状态数少的状态向微观状态数多的状态变化，从熵值小的状态向熵值大的状态变化，这就是热力学第二定律的本质。
3．热力学第三定律的本质
当 T→ 0 时，所有粒子都处于基态能级，此时?
0
=1，即把所有粒子放在一个能级上只有一个放法，体系只有一个微观状态，因此从玻兹曼定理，即式（ 6-25）可以得出结论：在 0 K 时物质的熵值为零，即
S
0
= k
B
ln?
0
= k
B
ln1 = 0 (4)
上式可以看作是热力学第三定律的统计表达式，这与热力学第三定律的表述,在 0 K 时任何纯物质的完美晶体的熵值为零,的结论是一致的。
第四节 配分函数的计算
一．平动、转动、振动的能级公式和能级的简并度
陕西师范大学精品课程 ……,物理化学,
第 23 页 共 40 页 2004-7-15
1,平动
质量为 m 的粒子，在体积为 a×b×c的方盒中运动，该三维平动子平动的能级公式为
()
2
222
y
xz
txyz 222
,,
8
=++
n
nnh
nnn
ma b c
ε （ 3.4?1）
其中 n
x
,n
y
,n
z
为 x,y,z 轴上的平动量子数，其取值为 1,2,3,…（正整数） 。
一组 n
x
,n
y
,n
z
的数值就决定一个平动量子态。如果运动空间是体积为 V 的立方盒，则该三维平动子平动的能级公式可写成
() ()
2
222
txyz x y z
23
,,
8
h
nnn n n n
mV
ε =++ （ 3.4?2）
三维平动子的基态能级是非简并的，其余以能级多数是简并的。如对在立方盒中运动的三维平动子，若某能级的能量
2
t
23
6
8
h
mV
ε = ×,则该能级有（ 2,1,1）,（ 1,2,1）,）
（ 1,1,2）三个量子态，即该能级的简并度为 3。
2、转动
线性刚性转子的能级公式为
() ( )
2
r
2
1
8
h
JJJ
I
ε
π
=+ （ 3.4?3）
其中 J 为转动量子数，其取值为 0,1,2,…； I 为转动惯量，对双原子分子
22
21
21
rr
mm
mm
I μ=
+
= （ 3.4?4）
其中 m
1
,m
2
为构成分子的两个原子的质量,μ 称为折合质量,r 为平衡时两原子的质心（核间）距离。
每一个由 J 决定的转动能级的简并度为
g
i,r
= 2J + 1 （ 3.4?5）
3、振动
单维简谐振子的能级公式为
()
v
11
22
hv hcvευ υ υ
=+ =+
（ 3.4?6）
其中 υ 为振动量子数，其取值为 0,1,2,…； v 为简谐振子的振动频率； C 为光陕西师范大学精品课程 ……,物理化学,
第 24 页 共 40 页 2004-7-15
速； v 为波数。
单维简谐振子的能级都是非简并的，即 g
i
,
v
=1 。三维谐振子的能量为
hv
zyx
+++
2
3
υυυ,能级的简并度为
g
i
,
v
=(S+1)(S+2)/2 （ 3.4?7）
其中
xyz
S υ υυ=++。
二．配分函数的析因性质
分子的运动可以分解为平动（ t）,转动（ r）,振动（ v）,电子运动（ e）和核运动（ n）
等。分子的各种运动形式可以看作是相对独立的，因此分子的能量可以认为是各种形式运动能量的加和，即
ε
i
= ε
i,t
+ ε
i,r
+ ε
i,v
+ ε
i,e
+ ε
i,n
各种运动形式分别具有其相应的简并度，分子总的简并度等于各个能级上简并度的乘积
g
i
= g
i,t
·g
i,r
·g
i,v
·g
i,e
·g
i,n
根据配分函数的定义，得
,t,r,v,e,n
trven
B
exp exp
++++
=?=
∑∑
i i i i i
i
i i,i,i,i,i,
ii
qg ggggg
kT k T
εεεεε
ε
数学上可以证明，几个独立变量乘积之和等于各自求和的乘积，于是
,t,r,v
,t,r,v
BB B
,e,n
,e,n
BB
trven
exp exp exp
exp e xp
=
∑∑∑
∑∑
i i i
i i i
iii
i i
i i
ii
qg g g
kT kT kT
gg
kT kT
qqqqq
εεε
εε
（ 3.4?8）
式中 q
t
,q
r
,q
v
,q
e
,q
n
分别称为平动配分函数、转动配分函数、振动配分函数、电子配分函数、核配分函数，分别代表分子的各种运动形式对分子配分函数的贡献。 （ 3.4?8）
式表明分子配分函数能够分解为各种形式运动贡献的乘积，这称为分子配分函数的析因子性质。应该指出，分子配分 函数的析因子性质只是对独 立粒子体系才是正确的。由
（ 3.4?8）式我们得到,
平动配分函数为
B
,t
t,t
=
∑
i
kT
i
i
qge
ε
（ 3.4?9）
转动配分函数为
陕西师范大学精品课程 ……,物理化学,
第 25 页 共 40 页 2004-7-15
B
,r
r,r
i
kT
i
i
qge
ε?
=
∑
（ 3.4?10）
振动配分函数为
vB
v,v
kT
i
i
qge
ε?
=
∑
（ 3.4?11）
电子配分函数为
B
,e
e,e
i
kT
i
i
qg
ε?
=
∑
（ 3.4?12）
核配分函数为
B
,n
n,n
i
kT
i
i
qge
ε?
=
∑
（ 3.4?13）
式（ 3.4?9）至式（ 3.4?13）中 g
i,t
,g
i,r
,g
i,v
,g
i,e
,g
i,n
和 ε
i,t
,ε
i,r
,ε
i,v
,ε
i,e
,ε
i,n
分别为平动、转动、振动、电子运动、核运动这五种运动形式能级的简并度和能级能量。
将（ 3.4?8 ）式应用于各热力学量的表达式，可得
2
B
,
222tvr
BBB
,,,
22e n
BB
,,
ln
ln lnln
ln ln
=
=++
++
VN
VNVN VN
VNVN
q
UNkT
T
qqq
Nk T Nk T Nk T
TTT
q q
Nk T Nk T
TT
(3.4?14)
=U
t
+U
r
+U
v
+U
e
+U
n
B
tvr
Bt Br Bv
e n
Be Bn
(loc) ln
ln ln ln
ln ln
N
NNN
NN
U
Skq
T
UUU
kq kq kq
TTT
U U
kq kq
TT
=+
=+++++
++++
(3.4?15)
=S
t
+S
r
+S
v
+S
e
+S
n
B
tt vr
BBB
e n
Be Bn
(non) ln
!
ln ln ln
!
ln ln
N
N
NN
NN
qU
Sk
NT
qU UU
kkqkq
NT T T
U U
kq kq
TT
=+
=+++++
++++
(3.4?16)
=S
t
+S
r
+S
v
+S
e
+S
n
使用类似的处理，其它热力学量均可以分解为 各种运动的贡献之和。应注意的是陕西师范大学精品课程 ……,物理化学,
第 26 页 共 40 页 2004-7-15
比较 （ 3.4?15） 和 （ 3.4?16） 式,可以看出定域子体系和离域子体系的平动熵 S
t
相差 k lnN!,
其余各项都是相同的。这是因为这两种粒子的区别就在于分子整体的平动，而分子的内部运动则是一样的。
三．各种运动形式的分子配分函数的计算及应用
根据分子配分函数的析因性质，只要求算出各 种运动形式的分子配分函数，就可以得到总配分函数，从而计算体系各热力学量。
1．平动配分函数
分子的平动，就是分子作为一个整体在容许的 空间范围内的运动，这相当于一个粒子了三维势箱中的运动，分子可以简化为三维平动子。根据量子力学原理，在长、宽、
高分别为 a,b,c 的三维势箱中，质量为 m 的三维平动子的能级公式为
2
22 2
y
x
txyz
222
(,,)
8
z
n
hn n
nnn
ma b c
ε
=++
式中 h 是普朗克常数,n
x
,n
y
,n
z
分别是 x,y,z 三个方向上的平动量子数，都只能取 1 到无穷大的任意正整数。当 n
x
,n
y
,n
z
一定时，就确定了一个平动量子态，量子态是非简并的。所以分子的平动配分函数为
xyz
xyz
2
222
y
xz
t
222
B
2
2 222
y
x z
BBB
xyz
exp
8
exp exp exp
888
=?+
=
=
∑∑∑
∑∑∑
nnn
nnn
n
nnh
q
mk T a b c
n
hhh
mk T a mk T b mk T c
qqq
式中 q
x
,q
y
,q
z
分别是三个方向上单维平动子的配分函数，各与一个平动自由度相对应。
如果令
2
2
2
B
8
=
h
mk Ta
α,则
xx
22
22x
2
11
B
exp exp( )
8
∞∞
==
==?
∑∑
nn
nh
qn
mk T a
α
一般来说,α
2
的值是很小的。当 α
2
<<1 时，一系列连续相差很小的数值求和，在数学上可以看作是连续的，因而可以用积分代替求和，即
x
22 22
xx x
0
1
exp( ) exp( )
2
n
qnnd
π
αα
α
∞
∞
=
=?=? =
∑
∫
陕西师范大学精品课程 ……,物理化学,
第 27 页 共 40 页 2004-7-15
对于 q
y
和 q
z
我们可作类似的处理。于是有
1/2 1/2 1/2
BBB
xyz
22
2
,,
πππ
=?=?=?
mk T mk T mk T
qaqbqc
hhh
所以
3/2 3/2
BB
t
22
=?=?
mk T mk T
q abc V
hh
ππ
（ 3.4?17）
其中 V= abc，是体系的体积。
将平动配分函数应用于理想气体（独立离域子体系）,有
22t
tB B B
,,
ln 3ln 3
22
==?=
VNVN
q T
UNkT NkT NkT
TT
（ 3.4?18）
t
,t B
,
3
2
==
v
VN
U
CNk
T
（ 3.4?19）
tt
tB B t B B
t
BB
3
ln ln ln !
!2
5
ln
2
=+=?+
=+
N
qU
Sk NkqkN Nk
NT
q
Nk Nk
N
（ 3.4?20）
（ 3.4?20）式称为沙克尔 (Sackur)-特鲁德 (Tetrode)方程，可直接用于计算理想气体的平动熵。平动对于其它热力学量的贡献亦可用类似的方法求得。
[例题 3.2] 计算 298.15 K,10
5
Pa 下,1 mol N
2
的平动配分函数。
解答,N
2
的摩尔质量为 14.008×2×10
－ 3
kg·mol
－ 1
，所以
326
23
14.008 2
10 kg 4.6515 10 kg
6.023 10
×
=×=×
×
m
h = 6.626×10
－ 34
J·s k
B
= 1.38×10
－ 23
J·K
-1
而 V = RT／ p = 0.02445 m
3
于是
3/2
29B
t 2
2
1.01 10
=?=×
mk T
qV
h
π
2．转动配分函数
对于双原子分子，除平动外，分子内部还有转 动和振动。一般说来，转动和振动互相有关，但为了简化，可近似认为二者相互独立，并将转动看作是刚性转子绕质心的运动。
根据量子力学处理的结果，线型刚性转子的能级公式为
2
r
2
() ( 1)
8
=+
π
h
JJJ
I
ε
陕西师范大学精品课程 ……,物理化学,
第 28 页 共 40 页 2004-7-15
式中 J 是转动量子数，可取 0,1,2,…等非负整数； I 是转动惯量，对于双原子分子
2
21
212
r
mm
mm
rI
+
== μ （ 3.4?21）
式中 m
1
,m
2
分别是两个原子的质量,μ 称为约化质量,r 是两个原子的质心距离。
分子的转动角动量的空间取向是量子化的，转动能级的简并度 g
r
= 2J+1，所以
2
r
2
0
B
(2 1) exp ( 1)
8
∞
=
=+?+
π
∑
J
h
qJ J
IkT
令 Θ
r
= h
2
/ 8π
2
I k
B
,称为分子转动的特征温度。通常温度下，对于大多数气体 Θ
r
<<T，
因此也可以用积分代替加和，即
r
r
0
r
2
B
2
(1)
(2 1)exp
8
∞
+
=+? =
π
=
∫
JJ T
qJ dJ
T
Ik T
h
Θ
Θ
（ 3.4?22）
（ 3.4－ 22）适用于异核双原子分子，也适用于非对称的线型多原子分子，如 HCN 等。
对于同核双原子分子或对称的线型多原子分子，由于分子对称性的原因
2
B
r
2
r
8π
==
IkTT
q
hσΘ σ
（ 3.4?23）
式中 σ 称为分子的对称数，即分子围绕对称由旋转 360°时具有相同位置的次数。对于同核双原子分子 σ =2，异核双原子分子 σ =1。
非线型多原子分子的情况比较复杂，可以导出
23/2 1/2
Bxyz
r
3
(8 ) ( )ππ
=
kT III
q
hσ
（ 3.4－ 24）
式中 I
x
,I
y
,I
z
分别是 x,y,z 三个轴向上的转动惯量。根据（ 3.4－ 22）,（ 3.4－ 23）,
（ 3.4－ 24）式计算分子的转动配分函数，其能量标度零点均选在转动的基态能级。
将分子的转动配分函数应用于双原子分子或线型多原子分子体系，有
22
r
rB B B
,
r
,r B
,
r
rBr Br B
ln ln
ln ln
===
==
=+= +
NV
v
VN
N
q T
UNkT NkT NkT
TT
U
CNk
T
U
Sk q Nk qNk
T
对于非线型多原子分子，有
陕西师范大学精品课程 ……,物理化学,
第 29 页 共 40 页 2004-7-15
rB
,r B
rBr B
3
2
3
2
3
ln
2
=
=
=+
v
UNkT
CNk
SNk q Nk
[例题 3.3] 已知 NO 的转动惯量 I = 16.4×10
－ 47
kg·m
2
，计算 298 K 时 NO 分子的转动配分函数
r
q 和该气体的摩尔转动热力学能 U
r,m
,摩尔转动熵 S
r,m
及摩尔定容热容 C
v,r,m
。
解答 NO 是异核双原子分子，对称数 σ = 1，所以
r
q = (8π
2
Ik
B
T) / h
2
= 8×3.14
2
×16.4×10
－ 47
×1.38×10
－ 23
×298 / (6.626×10
－ 34
)
2
= 121.2
U
r,m
= Lk
B
T = RT = 8.314 ×298 = 2478 (J·mol
－ 1
)
S
r,m
= Lk
B
ln
r
q + Lk
B
= Rln
r
q + R = 8.314 ln121.2 + 8.314 =48.2 (J·K
－ 1
·mol
－ 1
)
C
v,r,m
= Lk
B
= R = 8.314 (J·K
－ 1
·mol
－ 1
)
3．振动配分函数
分子的振动可以近似简化为若干个简谐振动,每个振动自由度相当于一个单维简谐振子。双原子分子只有一个振动自由度，可作为一维简谐振子处理。
根据量子力学的处理结果，单维简谐振子的能级公式为
v
1
()
2
hvευ υ
=+
式中 υ 为振动量子数，其取值为 0,1,2,…非负整数； v 为简谐振子的振动频率。当量子数
υ = 0 时 ε
v
（ 0） = 1/2 hv，称为零点振动能，亦即振动基态能级的能量。单维简谐振子的能级都是非简并的，于是
v
0
B
BBB
B
B
1
2
exp
2
exp 1 exp exp
2
1
exp
2
1exp
∞
=
+
=?
=+?+?+
=
∑
"
hv
q
kT
hv hv hv
kT kT kT
hv
kT
hv
kT
υ
υ
（ 3.4?25）
令 Θ
v
= hv/ k
B
,称为分子特征振动温度，则
陕西师范大学精品课程 ……,物理化学,
第 30 页 共 40 页 2004-7-15
v
v
v
1
exp
2
1exp
q
T
T
Θ
Θ
=
（ 3.4?26）
绝大多数气体的 Θ
v
都很高，通常温度下 Θ
v
>>T，因此上式可近似为
11
vB22
//
v
ee
==
ThkT
q
Θν
（ 3.4?27）
说明在通常温度下，气体分子几乎总是处于振动基态。只有当温度 T 接近 Θ
v
时，其它能级对配分函数才有实际贡献。当温度升高,T>>Θ
v
时，上述求和可用积分代替而得到
B
v
v
==
kTT
q
hvΘ
（ 3.4?28）
表明此时各振动能级对配分函数均有实际贡献。
若将能量零点选在振动基态能级，即令 ε
v
(0) = 0，则
v B
0v / /
11
()
1e 1e
==
T hkT
q
Θ ν
（ 3.4?29）
当 Θ
v
>>T 时,q(
0
)
v
= 1，于是可使计算大大简化。
对于由 n 个原子组成的多原子分子来说，分子具有多个振动自由度，总的振动配分函数应等于各单维简谐振子配分函数的乘积。
如果是线型多原子分子，总自由度为 3n，其中平动自由度为 3，转动自由度为 2，
振动自由度为 3n?5； 而对于非线型多原子分子,转动自由度为 3,振动自由度应为 3n?6。
于是振动配分函数可写为
B
0v /
1
()
1e
=
∏
i
f
hkT
i
q
ν
（ 3.4?30）
式中 f 为分子的振动自由度，其值为 3n?5 或 3n?6 。
即,
（ⅰ）线型分子
()
v
v
2
35
v
1
e
1e
=
′
=Π
T
n
T
i
q
Θ
Θ
（ 3.4?31）
或者
v
35
v
1
1
1e
=
=Π
n
T
i
q
Θ
（ 3.4?32）
（ⅱ）非线型分子
陕西师范大学精品课程 ……,物理化学,
第 31 页 共 40 页 2004-7-15
()
v
v
2
36
v
1
e
1e
=
′
=Π
T
n
T
i
q
Θ
Θ
（ 3.4?33）
或者
v
36
v
1
1
1e
=
=Π
n
T
i
q
Θ
（ 3.4?34）
式（ 3.4?31）和式（ 3.4?33）的能量标度零点与式（ 3.4－ 25）相同，式（ 3.4?32）和式
（ 3.4?34）的能量标度零点与式（ 3.4－ 29）相同。
[例题 3.4] 已知 NO 分子的特征振动波数 /ννc= =1907 cm
－ 1
，求 298 K 时该分子的振动配分函数和 NO 气体的摩尔转动能及摩尔振动熵。
解答,Θ
v
= hν/ k
B
= 6.626×10
－ 34
×(3×10
8
×1907×10
2
) / (1.38×10
－ 23
)=2747 K
因 Θ
v
(2747K) >> T(298K)
故 （ q
0
）
v
≈ 1
0v
m0,m B
,
ln( )
() 0
= ≈
NV
q
UU LkT
T
m0,m
v,m B 0 v
()
ln( ) 0
=+ ≈
UU
SLkq
T
4,电子配分函数
如果将能量零点选在基态，则电子的配分函数
,e 1,e
0e,e 0,e 1,e
BB
() exp exp
=?=+?+
∑
"
i
i
i
qg gg
kT kT
ε
ε
电子的能级间隔较大,?ε
i,e
一般几百 kJ·mok
－ 1
，如基态到第一激发态大约有 400 kJ·mok
－ 1
，粗略的估计，如果?ε
i,e
/ k
B
T > 5，则上式中第二项以后就可以忽略不计。除非在一个相当高的温度下，电子总是处在基态，而当温度增高时经常是在电子尚未激发之前分子就已经分解，因此上式第二项及其以后各项一般均可忽略，即
( q
0
)
e
= g
0,e
亦即电子配分函数等于电子运动基态的简并度。若电子运动的总角动量量子数为 j，对应 j 值有 (2j+1)个简并度。所以电子配分函数也可以写为
( q
0
)
e
= 2j + 1
应该注意的是，有少数原子如卤素原子和少数分子如 NO 等，它们的电子基态与第一激发态之间的能量间隔并不是很大，这时就需要考虑到第二项的贡献，即
陕西师范大学精品课程 ……,物理化学,
第 32 页 共 40 页 2004-7-15
1,e
0e 0,e 1,e
B
() exp
=+?
qgg
kT
ε
5,核配分函数
原子核运动的能级间隔更大，在通常的物理和化学变化中，原子核总是处在基态。
核配分函数的表达式为
B
,n
n,n
e
=
∑
i
kT
i
i
qg
ε
若将核运动基态能级的能量选作零，则上式变成
q
n
’ = g
0,n
核能级的简并度来源于原子核的自旋作用。若核自旋量子数为 s
n
，则原子或原子分子核运动基态能级的简并度为
g
0,n
= 2S
n
+ 1
对于多原子分子，其核运动基态能级的简并度是组成分子的各原子核运动基态能级的简并度之积，即
g
0,n
= ∏(2S
n
+ 1)
i
第五节 体系热力学性质的计算
上面我们对各种运动形式的配分函数求算作了讨论，有了 q 值就可具体计算体系的热力学函数值。这里主要讨论理想气体体系，即独立等同粒子体系。
一．分子的全配分函数
对于独立粒子体系而言，由于分子配分函数的析因子性质,分子的配分函数可以写为
trven
qqqqqq=
一般情况下，原子核的运动状态不会发生变化。而在某些化学反应中，电子转移会导致价电子运动状态发生变化，于是分子的全配分函数可以写为
trve
qqqqq=
通常在一般化学问题中，电子和核的运动状态都不会发生变化，电子配分函数和核配分函数没有必要进行计算，因此分子的全配分函数可以简化为
trv
qqqq=
综合以上对各种运动形式配分函数的讨论，在 不考虑电子和核运动的情况下，若采用基态能量为能量零点的规定，则各种不同分子的全配分函数可表示如下
陕西师范大学精品课程 ……,物理化学,
第 33 页 共 40 页 2004-7-15
(1) 单个原子
由于单个原子既无转动亦无振动，只需要考虑平动，故
3/2
B
0t
2
2π
==?
mk T
qq V
h
（ 3.5－ 1）
(2) 双原子分子
双原子分子有三个平动自由度，两个转动自由度和一个振动自由度，故有
B
3/2
2
BB
0 /
22
281
1e
ππ
=
hkT
mk T Ik T
qV
hh
ν
σ
（ 3.5－ 2）
由以上关系式知，只要知道了分子的质量 m、分子的转动惯量 I（ μr
2
）,对称数 σ
和振动频率 ν 等微观性质，就可以计算分子的全配分函数，从而求得宏观体系的热力学量。
二、体系热力学性质的计算
1．单原子分子
将单原子分子的配分函数表达式（ 3.5－ 1）代入不可别粒子体系的热力学性质的配分函数表达式中（纯气 体是不可别粒子体系）,可得到单原子分子气体体系相应的热力学性质。
(1) 亥姆霍兹自由能
e,0
e,0
t
net B n B eB
n,0
Bn,0 B
3/2
BB
3
n,0 e,0 B n.0 e,0
3/2
BBB
3
ln ln ln
!
[ exp( )] [ exp( )]
2
ln[( ) ] ln !
()l
2
ln[( ) ] ln
=++=
=
+
=+
+?
N
q
FFFF NkTq NkTqkT
N
Nk T g Nk T g
kT kT
mkT
Nk T V k T N
h
NN NkTgg
mkT
Nk T V Nk T N Nk T
h
ε
ε
π
εε
π
(3.5－ 3)
式中第一项是核和电子处于基态时的能量，第二项表示相应的间并度。在讨论热力学变量时，这些量是常量，都可以消去。
(2) 压力
B
Nk T
P
V
= (3.5－ 4)
(3) 摩尔热力学能 RTU
2
3
m
= (3.5－ 5)
(4) 摩尔焓 RTH
2
5
m
= (3.5－ 6)
陕西师范大学精品课程 ……,物理化学,
第 34 页 共 40 页 2004-7-15
(5) 摩尔恒容热容和摩尔恒压热容
RC
V
2
3
m,
= (3.5－ 7)
RC
P
2
3
m,
= (3.5－ 8)
(6) 摩尔熵
m
35
ln ln ln 9.685
22
r
p
SRMRTR
p
=+
\
(3.5－ 9)
(7) 摩尔吉布斯自由能
35
ln ln ln 30.47
22
+
r
p
GRT RM RTT T
p
m
=
\
(3.5－ 10)
式（ 3.5－ 10）称为沙克尔（ Sackur） ——特鲁德（ Tetrode）方程式。在应用式（ 3.5－ 9）
和式 （ 3.5－ 10） 计算惰性气体的 S
m
和 G
m
时,其中普适气体常数 R 必须取 8.314 J·K
－ 1
·mol
－ 1
，相对摩尔质量 M
r
的量纲为 g·mol
－ 1
，温度 T 的量纲为 K,压力 p 和标准压力 p
的量纲一定要相同。这样计算出来的 S
m
的单位为 J·K
－ 1
·mol
－ 1
,G
m
的单位为 J·mol
－ 1
。
(8) 化学势
对理想气体
B
=
Nk T
V
p
，代入 F 的表示式，得
n,0 e,0 B n,0 e,0
3/2
B
BBBB
3
()l
2
ln( ) ln ln
=+
π
+
kT g g
mk T
kT kT kT kT kT p
h
μ εε
对理想气体：
B
=
Nk T
V
p
对 1 摩尔气体而言，各项均乘以阿伏加德罗常数 L,Lk
B
=R,则 1 摩尔气体化学势为,
,0,0,0,0
2/3
B
B
2
()ln
2
ln( ) ln ln
=+
π
+
ne ne
LRTg
mk T
RTRkTRTp
h
μ εε
若气体的压力为标准压力 p
，则标准态的化学势 μ
可写为,
,
()
=
TV
F
N
μ
陕西师范大学精品课程 ……,物理化学,
第 35 页 共 40 页 2004-7-15
3/2
0,n 0,e 0 n 0 e
3
(2 )
( ) ln ln
ln ln
,,
mkT
LRTgRT
h
RT kT RT RT p
μεε
π
=+
+
\
\
对一定体系来说，上式右方都只是温度的函数，所以标准态的化学势可写作 μ
(T)，将上面两式相减，则得,
() ln( / )TRTppμμ=+
\ \
这就是气体化学势的表达式，这里也可以具体地看出 (T)中包含着哪些项。
2．双原子分子气体
将双原子分子的配分函数表达式（ 3.5－ 2）代入不可别粒子体系的各种热力学性质的配分函数表达式中，得到双原子分子气体体系相应的热力学性质。
(1) 热力学能
B
B
m
5
21
hv k T
hvkT
URT
e
=+
(2) 摩尔恒容热容
()
BB
B
2
2
,m
B
2
B
5
1
2
5
2
=+?
=+
hv k T hv k T
V
hv k T
hv
CRR ee
kT
hv
RR e
kT
显然，在室温时大多数双原子分子气体的 RCRC
PV
2
7
,
2
5
m,m,
== 。
(3) 熵尔熵
其中
t
S,
r
S,
v
S 为平动、转动、振动运动形式对摩尔熵的贡献。并且
()
()()
BB
32
B
B
t 3
2
B
r 2
1
v
B
2
5
ln
2
8
ln
ln 1 e e 1
π
=?+
π
=+
= +?
hv k T hv k T
mk T
kT
SR R
hP
Ik T
SR R
h
hv
SR R
kT
σ
(4) 吉布斯自由能
陕西师范大学精品课程 ……,物理化学,
第 36 页 共 40 页 2004-7-15
()
()
B
32
2
B
B
mm
2
8
ln 1 1
π
π
=
hv k T
mk T
Ik T
GRT V e
Lh hσ
其中 L 为阿伏加德罗（ Avogadro）常数。在计算熵时要注意，经典统计对粒子的不可别性只在平动运动形式的贡献项内进行修正。
[例题 3.5] 求氩（ Ar）气在其正常沸点 87.3 K 和标准压力时的摩尔热力学能
m
U,摩尔熵
m
S 及定圧摩尔热容
,mp
C 。
解答 氩为单个原子，可作为理想气体
3/2
20 3/ 2B
t 2
2
1.879 10 ( )
π
==?= ×
mk T
qq V MT V
h
Ar 的摩尔质量 M = 39.9 g·mol
－ 1
，标准压力 p
时的摩尔体积
m
V = RT / p
= 8.314×87.3 / 101325 = 7.16×10
－ 3
(m
3
)
故 q = 1.879×10
20
×(39.9×87.3)
3/2
×(7.16×10
－ 3
×10
6
) = 2.766×10
29
m
U = (3/2)RT = (3/2)×8.314×87.3 = 1089 (J·mol
－ 1
)
,mv
C = (3/2)R = (3/2)×8.314 = 12.47 (J·K
－ 1
·mol
－ 1
)
,mp
C = C
v,m
+ R = 20.87 (J·K
－ 1
·mol
－ 1
)
m
S = R ln( q /L) + (5/2)R = 8.314×[ln(2.766×10
29
／ 6.023×10
23
) +5/2]
= 129.2 (J·K
－ 1
·mol
－ 1
)
根据热力学第三定律，采用量热法得到的结果是 129.1J·K
－ 1
·mol
－ 1
。
三、理想气体状态方程的导出
将
B
ln
!
=?
N
q
FkT
N
代入 ()
T
F
p
V
=?
中，得
B
ln
()
=?
T
q
pNkT
V
又
trven
3/2
B
rven2
ln ln ln( )
2
ln ln ln( )
=+
π
=++
qq qq
mk T
Vqqq
h
式中，右端首尾两项与体系的体积 V 无关，故有,
BB
,
ln 1?
==?
TN
q
pNkT nkT
VV
陕西师范大学精品课程 ……,物理化学,
第 37 页 共 40 页 2004-7-15
而 Nk
B
=(N/k
B
)·Lk
B
=nR,L 为阿佛加德罗常数。所以
pV=nRT
四,统计熵的计算
在统计热力学中，对于非定位体系，有,
S = S
t
+ S
r
+ S
v
+ S
e
+ S
n
其中
平动熵
tt
tB B,
() ln
ln ( )
!
=+
N
VN
Sk NkT
NT
转动熵
rr
rB B,
() ln
ln ( )
!
=+
N
VN
Sk NkT
NT
振动熵
vv
vB B,
() ln
ln ( )
!
=+
N
VN
Sk NkT
NT
电子熵
ee
eB B,
() ln
ln ( )
!
=+
N
VN
Sk NkT
NT
核 熵
nn
nB B,
() ln
ln ( )
!
=+
N
VN
Sk NkT
NT
式（ 6-110）表明，体系的熵值是物质各种运动形式对熵的贡献之和。
对于一般的物理及化学体系，可以忽略 S
n
，对少数物质需考虑 S
e
的贡献，所以统计熵主要来自 S
t
,S
r
和 S
v
的贡献。
1,平动熵
tt
tB B,
() ln
ln ( )
!
=+
N
VN
Sk NkT
NT
应用平动配分函数，可将上式化简为,
3/2B
tB
2
25
[ln( ) ]
2
π
=+?
mk T V
SNk
hN
上式称为 Sackur-Tetrode 方程，表示平动运动对熵的贡献。
当理想气体的物质的量为 1 mol 时，有
3/2Bm
tB 2
3/2BB
B 2
25
[ln( ) ]
2
25
[ln( ) ]
2
π
=+?
π
=+?
mk T V
SNk
hL
mk T k T
Nk
hp
陕西师范大学精品课程 ……,物理化学,
第 38 页 共 40 页 2004-7-15
应用 SI 单位,k
B
=1.381×10
－23
J·K
-1
,h = 6.626×10
－34
J·s，若压力为标准压力 p
，则上式成为
S
m,t
= R(3/2 ln M + 5/2 lnT – 1.165)
式中 M 为摩尔质量，单位为 g·mol
－1
。
2,转动熵
对线型分子
当 Θr/T≤0.01 时，将
r
q 代入 S
r
的表达式，得
rB
r
(1 )=+
T
SNk
σΘ
3,振动熵
对双原子分子
将振动配分函数代入 S
v
的表达式，得
v
v
/v
vB/
/
[ln(1e)]
e1
=?
T
T
T
SNk
Θ
Θ
Θ
4,量热熵和光谱熵
在分子的各种运动形式中，平动、转动和振动等称为热运动形式，它们对体系熵的贡献称为热熵。体系的热能和热熵都能在体系升降温度时吞吐。在一般的温度范围内不能吞吐能量，而能使体系内粒子的位置或构型的无序运动形式所贡献的熵称为构型熵。
其中包括分子的电子运动（但对 NO
2
O
2
等可以是热熵）,核自旋运动及一系列甚至未被认识的不可穷尽的运动形式。构型熵的来源是作为晶体前身的气体的热熵，在气体凝聚成晶体的过程中，有那么一部分热熵变成了晶体的构型熵。在量热实验中，我们可以测定一个体系在降温过程中倾吐出来的热熵。包括平动熵、转动熵、振动熵及构型熵转化的热熵的总和，这样得出熵值称为体系的量热熵。就气体来说，分子的平动、转动和振动等热运动都会贡献出相应的热熵；对晶体物质来说，再加上源于热熵的构型熵，它们的总和称为气体或晶体的光谱熵。简单地说,
(1) 量热熵（或规定熵）
在热力学第三定律的基础上，根据热力学第二定律，利用量热测出的各种 H
p
和 △ H
数据而算得的 TK 时的熵叫量热熵。
陕西师范大学精品课程 ……,物理化学,
第 39 页 共 40 页 2004-7-15
mb
mb
0
0
0
d
d d d
=?=
=+++
∫
∑
∫∫∫
T
P
TT
TTT
PPP
TT
C
SSS T
T
CCC H
TTT
T
(2) 光谱熵（或统计熵）
根据统计热力学的方法，利用光谱实验数据先算出转动、振动、电子等运动形式的配分函数，从而得出的各种熵值，再加上平动熵而得到 TK 时的总熵值叫光谱熵。
在一般温度范围内吞吐的即是热熵,则我们在量热实验中测量的就只能是热熵了。 虽然分子内部的全部运动形式不可穷尽，但在一定温度范围内表现出来的热运动形式实际上已经明确。这也是我们计算统计熵的根据，即我们在实验中测量和在理论中计算以及使实际上已经明确。这也正是我们计算统计熵的根据，即我们在实验中测量和在理论中计算以及使实验和理论相互核对的正是热熵。而热熵与构型熵是可以相互转化的，当我们给一个热力学体系降温至复写纸的凝固点（
m
T ）时，体系会凝为固体，这是我们所能观测到的。在体系的内部实际上不存在着两个我们观测不到的温度：转化温度（
t
T ）和冻结温度（
f
T ） 。晶体中粒子从无序的位置或取向构型 可逆地转化为有序的位置或取向构型的温度称为转化温度；在某温度以下，由于动力学的原因（向有序的转化需较高的活化能）,体系中粒子的位置或构型将受到冻结，该温度称为冻结温度（一般
f
T ＜
m
T ＝。
对于转化温度高于冻结温度的晶体；可逆降温至温度达到转化温度和以下时，体系内部完成从无序向有序的可逆转化。这种情况下，来自气体热熵的构型熵又重新转化为能在降温过程中倾吐出来的热熵，即气体在降温过程中将它原有的全部热熵（包括源于气体熵的晶体的构型熵）作为量热熵倾吐出来。在这种情况下，光谱熵和量热熵在实验的误差范围内是一致的。对转化温度低于冻结温度的晶体，在降温过程中，由于动力学的原因，它的位置或取向构型的无序状态（该无序的位置或取向构型的总数? 可由统计的方法计算得到）被冻结，即被冻结为? 个无序状态中的任意一个状态。晶体没有机会将这种无序的位置或取向的构型熵（ S=k
B
ln?）重新以热熵的形式倾吐出来，来自它的前身
（气休） 热熵的这部分构型熵被永远保留在晶体中 （尽管被冻结的只是? 个状态中一个，
但被冻结的构型熵仍应以 S=k
B
ln?计算,即将被冻结的那个状态视为?个无序状态的,全权代表,） 。在这种情况下，我们由量热实验得到的量热熵只是体系全部热熵的一部分，
而通过计算得到的光谱熵更准确。这就是我们所说的统计热力学的一大成功之处：统计陕西师范大学精品课程 ……,物理化学,
第 40 页 共 40 页 2004-7-15
热力学关于低温熵的计算值比量热实验的测值更准确。简单地说,
一般情况下，量热熵和光谱熵能较好地吻合。个别物质在低温下，光谱熵大于量热熵，光谱熵更准确。
