第九章统计热力学初步物理化学
Statistical Thermodynamics
学习要求:
明确统计热力学的基本假设,理解最概然分布与平衡分布及摘取最大项原理掌握 Boltzmann 分布律及其各物理量的意义与适用条件;理解粒子配分函数、体系配分函数的意义与表达式,配分函数的析因子性质。
理解不同独立子体系的配分函数,q 及 Θ与热力学函数间的关系。
重点掌握平动能与平动配分函数,转动能与转动配分函数,振动能与振动配分函数的计算。
理解系统的热容、熵及其他热力学函数与配分函数的关系。
第九章 统计热力学初步引言
§ 9.5 粒子配分函数的计算
§ 9.3 最概然分布与平衡分布
§ 9.4 玻耳兹曼分布
§ 9.2 能级分布的微态数及系统的总微态数
§ 9.6 系统的热力学能与配分函数的关系
§ 9.10 理想气体反应的标准平衡常数
§ 9.11 系综理论简介
§ 9.1 粒子各运动形式的能级及能级的简并度
§ 9.7 系统的摩尔定容热容与配分函数的关系
§ 9.8 系统的熵与配分函数的关系
§ 9.9 其它热力学函数与配分函数的关系引 言统计热力学的研究方法物质的宏观性质本质上是微观粒子不停地运动的客观反应。虽然每个粒子都遵守力学定律,
但是无法用力学中的微分方程去描述整个体系的运动状态,所以必须用统计学的方法。
根据统计单位的力学性质,经过统计平均推求体系的热力学性质,将体系的微观性质与宏观性质联系起来,这就是统计热力学的研究方法。
① 为系统的热力学量及热力学量之间的关系提供微观解释,反过来也使得从系统宏观热力学性质推测系统的微观结构成为可能;
② 可以直接从系统内部粒子的微观运动性质及结构数据计算得出平衡系统各种宏观性质的具体数据。
统计热力学是联系微观与宏观的桥梁统计热力学的基本任务根据对物质结构的某些基本假定,以及实验所得的光谱数据,求得物质结构的一些基本常数,如核间距、键角、振动频率等,从而 计算分子配分函数。再根据配分函数求出物质的热力学性质,这就是统计热力学的基本任务。
定域子系统和离域子系统粒子(子) ( particles) —— 聚集在气体、液体、固体中的分子、原子、离子等。
按粒子运动情况不同,可分为,
定域子系统 ( system of localized particles)
(或称为可别粒子系统 system of distinguishable
particles)
离域子系统 ( system of non-localized particles)
(或称为等同粒子系统 system of indistinguishable
particles)
离域子系统,粒子是不可区分的(气体、液体)。
例如,气体的分子,总是处于混乱运动之中,彼此无法分辨,所以气体是非定位体系,它的微观状态数在粒子数相同的情况下要比定位体系少得多。
定域子系统,粒子是可以区分的(固体),例如,在晶体中,粒子在固定的晶格位置上作振动,每个位置可以想象给予编号而加以区分,
所以定位体系的微观态数是很大的。
定域子系统和离域子系统独立子系统和相依 (倚 )子系统按粒子间相互作用情况不同,可分为:
独立子系统 ( system of independent particles)
—— 粒子之间除弹性碰撞之外,无其它相互作用
(理想气体)。
相依 (倚 )子系统 ( system of interacting particles)
—— 粒子之间存在相互作用(实际气体、液体、
固体)。
本章只讨论 独立子系统 。
量子态的确定在总粒子数为 N、总能量为 U、体积为 V的独立子系统中,有定态薛定谔方程
( 1)根据测量原理,系统的总能量 U为上式的本征值;所有系统所允许的量子态均为对应本征值 U
的简并态。
1 2 1 2?H ( r,r,) E ( r,r,)NN?,r?,r
i i i i i i?H ( r ) ( r )
而由单个粒子的定态薛定谔方程可得
1
N
i
i
E?

12
1
N
ii
i
( r,r,) ( r )
N?,r
( 2)对于独立子系统,由于粒子间无作用力,
各粒子相互独立,因此系统的哈密顿算符可分离为各粒子哈密顿算符之和。
1
N
i
i
HH

( 3)对于全同粒子系统,每个粒子的哈密顿算符形式等价,因而具有完全相同的本征值集合,合并得到
Od
O
d


原则上,对于给定的独立子系统,只要知道单粒子定态薛定谔方程的解,求得分布数,
就可以得到系统的波函数。
系统处在该量子态时任意可观测物理量的平均值由下式给出,
这几个能级的大小次序是:
§ 9.1粒子各运动形式的能级及能级的简并度
t r v e n
一个分子的能量内部运动的能量平动能 (εt)
转动能 (εr )
振动能 (εV )
电子的能量 (εe )
核运动能量 (εn )
42-420J·mol-1
4.2-42KJ·mol-1
更高
4.2× 10-21J·mol-1
若分子中各运动形式可近似认为彼此独立,则分子的能量等于各独立运动形式所具有的能量之和:
同时,其简并度等于各独立运动形式的简并度之积:
t r v e n
t r v e ng g g g g g
运动自由度对于一个具有 n个原子的分子,通常有 3n个自由度,分别为:
3个平动自由度( xyz轴方向的平动)
3个转动自由度(围绕三个轴的旋转)
3n-6个振动自由度对于线型分子,转动自由度为 2(围绕线轴的旋转可忽略),振动自由度为 3n-5
2
2 2 2
t 328 x y z/
h ( n n n )
mV

2 222
t 2 2 28
y zx n nnh ()
m a b c

量子力学中把能级可能有的微观状态数称为该能级的 简并度 (degeneration),用符号 g表示。
简并度亦称为退化度或 统计权重 。
对立方容器 a=b=c,V=a3
1.三维平动子三维平动子例如,气体分子平动能的公式为:
2
2 2 2
t 328 x y z/
h
( n n n )
mV

式中 分别是在 轴方向的平动量子数,当,则只有一种可能的状态,则,是非简并的。
x y zn,n n和 zyx 和,
2
t 32 38 /
h
mV

1ig?
1 1 1x y zn,n,n
x y zn n n
这时,在 相同的情况下,有三种不同的微观状态,则 。
i?
3ig?
t
2
3 / 2 68
h
mV当 时
2 1 1
1 2 1
1 1 2
J,转动量子数,J= 0,1,2,…
I—— 转动惯量 ( moment of inertia),与结构有关,
数值可由光谱数据获得。
对于双原子分子,有式中,R0 = r1 + r2,μ —— 折合质量 ( reduced mass)
20IR 12
12
mm
mm?

2.刚性转子
2
r 2( 1 ) 0 1 28
hJ J J
I
,,,
r 21,JgJ
简并度只考虑双原子分子刚性转子
3.一维谐振子
v ( ) h,,,
1 0 1 2
2

ν—— 粒子的振动频率,与结构有关,数值可由光谱数据获得。
υ—— 振动量子数 υ= 0,1,2,
V 1,g
一维谐振子电子运动及核运动的能级差一般都很大,一般的温度变化难以产生能级的跃迁或激发,所以本章只讨论最简单的情况,即一般认为系统中各粒子的这两种运动处于 基态 。
不同物质电子运动基态能级的简并度 ge,0和核运动基态能级的简并度 gn,0可能有区别,但对指定物质而言,都为常数。
4.电子和原子核
§ 9.2 能级分布的微态数及系统的总微态数分布数 —— 任一能级 i上的粒子数目 ni称为能级 i上的分布数。
能级分布 —— N个粒子在各个能级上的分布,
称为能级分布,简称分布。
要说明一种能级分布,就要阐明各能级上的粒子分布数。
1.能级分布 (energy level distribution)
例,一个定域子系统中,只有 3个一维谐振子,它们分别在 ABC三个定点上振动,总能量为 9hν/2,
即,N = 3,U = 9hν/2。由 知,系统中粒子可能处于的能级有,ε0= hν/2,ε1=3 hν/2,ε2=
5hν/2,ε3= 7hν/2,不可能有粒子处于能量大于 ε3的能级上,否则,系统的总能量会超过 9hν/2。
v
1
2( )h
系统的可能的能级分布方式有:
能级分布 能级分布数 Σni Σniεi =9hν /2n
0 n1 n2 n3
Ⅰ 0 3 0 0 3 3× 3 hν /2=9hν /2
Ⅱ 2 0 0 1 3 2× hν /2+1× 7hν /2=9hν /2
Ⅲ 1 1 1 0 3 1× hν /2+1× 3hν /2
+1× 5hν /2=9hν /2
2.状态分布 (state distribution)
能级分布只说明在各个能级上分布的粒子数,但在能级有简并度或粒子可以区别的情况下,同一能级分布还可以对应不同的状态分布。
状态分布 —— 粒子在各量子态上的分布。
显然,要描述一种状态分布,就需要知道各个量子上的粒子数 —— 状态分布数 。
那么,一种能级分布要用一定数目的几套状态分布数来描述。
υ= 3
υ= 2 ● ● ● ● ● ● ● ● ●
υ= 1 ●●● ● ●● ● ● ●
υ= 0 ●● ● ● ●●● ● ● ● ● ●
微态 1 2 3 4 5 6 7 8 9 10
能级 Ⅰ Ⅱ Ⅲ
还是上例上例中,Ω = Σ WD = 1+3+6 = 10
D
D
W
一般又将粒子的量子态称为 微观状态,简称 微态 。
显然,一种能级分布 D有一定的微态数 WD,全部能级分布的微态数之和即为系统的总微态数 Ω。
计算一种能级分布的微态数的本质
—— 排列组合问题 。
由于定域子系统和离域子系统中,粒子存在是否能区分的问题,其 WD的计算也有所不同。
3.定域子系统能级分布微态数的计算
i
i
ii
i
NN
NU?
全部能级分布的微态数之和即为系统的总微态 Ω。
D
D
W
一个由 N 个可区分的独立粒子组成的宏观体系,
在量子化的能级上可以有多种不同的分配方式。
无论哪种分配都必须满足如下两个条件:
能级的简并度是 1,任何能级的分布数是 1,
即 N个粒子分布在 ε 1~ ε n 共 n个不同能级上,
则 这种分配的微态数 即 N个粒子的全排列,

WD=N!( 能级 Ⅲ )
12
12
i
iN N N


能级:,,,
一种分配方式:,,,
12
i
i
N!
N ! N !
N!
N!


12
1
NN
D N N NW C C 1
1 1 2 1 2
( N N ) !N!
N ! ( N N ) ! N ! ( N N N ) !


(能级 Ⅰ,Ⅱ )
能级的简并度是 1,能级的分布数是
n1,n2,…,由于同一能级上各粒子的状态数相同,则 这种分配的微态数为:
12
12
12
i
i
i
,,,
g,g,,g
N,N,,N



能级各能级简并度一种分配方式各能级的简并度是 g1,g2,…,能级的分布数是 n1,n2,…,由于同一能级的粒子可处于不同量子态,则将 N1个粒子放在 能级上,共有 种微态数。
111 NNN Cg?1?
先从 N个分子中选出 N1个粒子放在 能级上,
有 种取法;
1?
1NNC
但 能级上有 个不同状态,每个分子在能级上都有 种放法,所以共有 种放法;
1? 1g 1?
1g 11Ng
1 1 2 2
112
N N N N
D N N NW ( g C ) ( g C )
12 1
12
1 2 1 2
NN ( N N ) !N!ggN ! ( N N ) ! N ! ( N N N ) !
12
12
1 2 i
NN N!gg
N ! N ! N !
iN
i
i i
gN!
N!
4.离域子系统能级分布微态数的计算假设每个量子态所容纳的粒子数没有限制任一能级非简并,gi=1,由于每个能级的简并度为 1,即只有一种排列方式,而且粒子是不可区分的,所以,在任一能级上粒子的分布方式只有一种。这样,系统某一分布 D
的微态数 WD = 1。
能级,ε1,ε2,ε3,…,εk
粒子数,n1,n2,n3,…,nk
简并度,g1,g2,g3,…,gk
对第一个能级而言,就是 n1个粒子分布在 g1个量子态上,则该能级上的微态数就是 n1在 g1个量子态上的分布方式数,若将 g1个量子态看成 g1个相连的盒子,那么,该能级上的微态数就是将 n1个粒子放在
g1个盒子中有多少可能的方式的问题,该问题可转化为 n1个粒子与 g1-1个盒壁混合在一起的排列数问题,而 n1个粒子与 g1-1个盒壁混合一起的全排列数为 [n1+( g1-1) ]!
11
11
1
1
( n g )!
( g )! n !


1 1 1 1 1 1 1 1
11
1 2 1
1
( n g ) ( n g ) ( n g n ) ( g ) !
( g ) ! n !

11()gn若
1
1
1
ng
n!?
1
1
ii
D
i ii
( n g ) !W
( g ) ! n !


in
i
i i
g
n!
所以,该分布的微态数为,
由于粒子是不可区分的,盒壁互换也不会影响分配方式,所以能出现的方式数为:
5.系统的总微态数
D
D
W
总微态数 Ω也有定值,即
( N,U,V )
N,U,V确定的系统能有哪些分布方式是确定的,各分布方式的微态数 WD是可计算的,
§ 9.3 最概然分布与平衡分布
1.概率 ( probability)
概率指某一件事或某一种状态出现的机会大小。
热力学概率体系在一定的宏观状态下,可能出现的微观总数,通常用 Ω 表示。
2.等概率定理对于 U,V和 N确定的某一宏观体系,任何一个可能出现的微观状态,都 有相同的数学概率,
这称为 等概率原理 。例如,某宏观体系的总微态数为 Ω,则每一种微观状态 P出现的数学概率都相等,即:
1P
3.最概然分布最概然分布 ( the most probable distribution)
若某种分布的微态数是 WD,则该分布出现的概率是,PD = WD /Ω
那么,在指定 N,U,V条件下,微态数最大的分布出现的概率最大。
所以,微态数最大的分布 —— 最概然分布。
4.最概然分布与平衡分布
D
N!W
M ! ( N M ) !
设某独立子系统中有 N个粒子分布于同一能级的 A、
B两个量子态上,当 A上粒子数为 M时,B上的粒子数为 N-M。因粒子可区别,则上述分布方式的微态数为:
不同的 M值代表着不同的分布方式,其中,
M=N/2时 WD最大,为最概然分布 (WB),
系统总的微态数为 Ω,它们分别为:
22B
N!W
( N / ) ! ( N / ) !
00
2
NN
N
D
MM
N!W
M ! ( N M ) !



随着 N的增大,系统微态数增加,但因总微态数增加的程度更大,所以,最概然分布的数学概率反而下降。在偏离最可几分布同样范围内,各种分布的数学几率之和随着 N的增大而增大。
取 N=10及 N=20两种情况在 N,U,V确定的系统达平衡时可以认为粒子的分布方式几乎不随时间而变,这种分布被称为 平衡分布 。反过来说,平衡分布就是最概然分布所能代表的那些分布。
DB/pp
1
0 105.
/MN
N=10,N=20时独立定域子在同一能级 A、
B两个量子态上分布的 PD/PB-M/N图
§ 9.4 玻耳兹曼分布
e ej i/ kT / kTj i i
j j i i
N n N n g
ee j i/ k T / k Ti
ij
NN
g


玻耳兹曼对独立子系统的平衡分布做了定量的描述。
e j / kTjn (按量子态计算)
e i / k Ti i j in g n g (按能级计算)
1.玻耳兹曼分布( Boltzmann distribution)
粒子配分函数 (partition function)q
de f
e j kT
j
q (按量子态计算)
d e f
e i kTi
i
qg (按能级计算)
因此得玻耳兹曼公式:
e ej i/ k T / k Tj i iNNn n gqq
符合此式的分布方式称为玻耳兹曼分布。
实质上 玻耳兹曼分布 即是最概然分布,故它可以代表平衡分布 。
e
e
i
k
kT
ii
kT
kk
ng
ng

由上式可以得出任何两个能级 i,k上粒子分布数 ni、
nk之比为,
任一能级 i上分布的粒子数 ni与系统的总粒子数 N之比为
ee
e
ii
i
/ k T / k T
i i i
/ k T
i
i
n g g
N g q



式中 称为能级 i的有效状态数或有效容量。e i / kTig
2.玻耳兹曼分布的推导
D l n l n l n l n
in
i
D i i i
iii
gW N ! W N ! ( n g n ! )
n!
由斯特林近似式 l n l nN ! N N N
以定域子系统为例玻耳兹曼分布需满足的三条件:
1 0i
i
nN

2 0ii
i
nU
DDd d l n
dd00ii
WW
nn
l n l n l n l nD i i i i i
i
W N N N ( n g n n n )
则:

1
2
0
0
i
i
ii
i
( n N )
( n U )




ln
l n l n
12D
i i i i i
i
i i i
ii
ZW
N ln N N ( n g n n n )
( n N ) ( n U )






则:
Z为极值时,则:iZ / n 0
l n l n e e 10 ii i i i ig n n g kT

DDl n l nl n l n
2
2
1 0
ii
i i i
WWgn
n n n


则,lnWD是极大值,即 WD是极大值,其对应分布为最概然分布或平衡分布、玻耳兹曼分布。
e e e
e
i
i
/ k T
ii / k T
ii i
i
NNN n g
gq


i / k T
ii
Nn g e
q

独立子系统中粒子的任一能级 i的能量值
t r v e ni,i,i,i,i,i
该能级的简并度
t r v e ni,i,i,i,i,ig g g g g g?
§ 9.5粒子配分函数的计算
1.配分函数的析因子性质
t r v e n
t,r,v e n e
,i,i,i,i,i / k T
i i,i,i,i
i
q g g g g g
t,r,
v,e,n
tr
v e n
ee
e e e
ii
ii,i
/ kT / kT
,i,i
ii
/ kT / kT / kT
,i,i,i
i i i
gg
g g g












代入 e
i kT
i
i
qg 得:
q= qt qr qv qe qn lnq=lnqt+lnqr+lnqv+lnqe+lnqn
上两式称为配分函数的 析因子性质 。
t
r
v,
e
n
tt
rr
vv
ee
nn
e
e
e
e
e
,i
,i
i
,i
,i
kT
,i
i
kT
,i
i
kT
,i
i
kT
,i
i
kT
,i
i
qg
qg
qg
qg
qg


平动配分函数转动配分函数振动配分函数电子运动配分函数核运动配分函数设粒子的能级为 ε 0,ε 1,ε 2,…,则
0i1
01e e e
kTk T / k T
i
i
q g g g
若选择粒子的基态能级作为能量的零点,且某独立运动形式基态能级的能量为 ε0,能级 i的能量为 εi,则以基态作为能量零点时能级 i的能量 εi0为:
2.能量零点的选择对配分函数 q值的影响
0 0ii
则上述各能级的能量分别为
00 0 0,02 2 0,01 1 0,
上式表明,选择不同的能量零点会影响配分函数的值。
但能量零点的选择对计算玻耳兹曼分布中任一能级上粒子的分布数 ni是没有影响的。因为令
00 e
i kT
i
i
qg
000 e e eii( ) kT / kT / kT
ii
ii
q g ( g )
00 e kTqq
则:
或 00 e kTqq
e i kTiiNng
q
0 0
00
e
e
i kT
ikT
N g
q



0
0 e
i kT
i
N g
q

t,0
r,0
v,0
e,0
n,0
0
tt
0
rr
0
vv
0
ee
0
nn
e
e
e
e
e
kT
kT
kT
kT
kT
qq
qq
qq
qq
qq



平动配分函数转动配分函数振动配分函数电子运动配分函数核运动配分函数因 εt,0≈0,εr,0=0,故常温下 qt0≈qt,qr0=qr。但 qv0和 qv、
qe0和 qe,qn0和 qn的差别不能忽略。
t e
j kT
j
q
对三维平动子 2 222
t 2 2 28
y zx n nnh
m a b c



3,平动配分函数的计算
2 222
t 2 2 2e x p 8
x y z
y zx
( n,n,n )
n nnh
q
m k T a b c



22 2222
2 2 2
1 1 1
e x p e x p e x p
8 8 8
x y z
y zx
n n n
hn hnhn
m k T a m k T b m k T c






t t t,x,y,zq q q?
qt,x,,qt,y,qt,z是一维平动子的配分函数 。
22
t 2
1
22
t 2
1
22
t 2
1
e x p
8
e x p
8
e x p
8
x
y
z
x
,x
n
y
,y
n
z
,z
n
hn
q
m kTa
hn
q
m kTb
hn
q
m kTc










以 qt,x为例求 qt:
令 2
2
2A 8
h
kT a
A是常数,对通常温度和体积条件下的气体,
A2<<1。则:
2 2 2 2
22
t 2 10
1
e x p e d e d8 xx
x
A n A nx
,x x x
n
hnq ( ) n n
m k T a






1
2
t
1
2
t
1
2
t
2
2
2
,x
,y
,z
m kT
qa
h
m kT
qb
h
m kT
qc
h


h:普朗克常数,h=
6.626× 10-34J·S;
玻耳兹曼常数 k=R/L
=1.38× 10-23J.K-1

3 3
2 2
t t t t 32
2 2
,x,y,z
m k T m k T
q q q q a b c V
hh



各项分别积分后得:
如果用 ft表示立方容器中平动子一个平动自由度的配分函数,则
3
ttqf?
1
12
3
t 2
2 m k T
fV
h



ft也是量纲为一的量 。 ( q是 量纲为一的量 )
对理想气体,整理后的平动配分函数计算式为
qt=8.2052× 107N(M/kg·mol-1)3/2(T/K)5/2/(P/pa)
对于 双原子分子 2rr
21 2 18
hJ ( J ) g J
I
式中 Θr =h2/8π2Ik 其单位为 K,由光谱得到的 I值算出,称为 转动特征温度 (rotational characteristic
temperature)。
4.转动配分函数的计算
r,
r r,e
i kT
i
i
qg
π
2
r 2
0
2 1 e x p 1 8
J
hq ( J ) J ( J )
Ik T

则在通常温度下,T>>Θr,可用积分代替加和,得:
r1r
r 0
0
2 1 e x p 1 2 1 e dJ ( J ) / T
J
q ( J ) J ( J ) ( J ) JT




π 2
r 2
r
8
T I k Tq
Θh
考虑到线性分子旋转一周会有 σ次不可分辨的几何位置,则:
π 2
r 2
r
8
T I k Tq
Θh

式中 σ称为 对称数
σ为 对称数,
对同核双原子分子,σ= 2;异核双原子分子,σ= 1。
因转动基态能级 εr,o=0,故 qr0=qr 。
将式中各常数的数值代入上式,最终可得:
qr=2.483× 1045(I/㎏ ·m 2)(T/K)/?
双原子分子的转动自由度数为 2,以 fr表示每个转动自由度配分函数的几何平均值,则
1
1 2
2
rr
r
Tfq



5.振动配分函数的计算
v
2
vv
00
1e e x p e e
2
,i hh
k T k T k T
,i
i
q g ( ) h / k T







对一维谐振子,gv,i=1,εv=(υ+1/2)hν,则其单位为 K(可由光谱得到的 ν 值算出 ),故称为 振动特征温度 (vibrational characteristic temperature) 。
一般情况下 Θv>>T(即 hν/k>>T)。
v
d e f
Θ h / k?令令 e -Θ v /T=x,则:
v v v22 2
v
0
e e e 1/ T / T / Tq ( x x

)
式中 0<x<1,故( 1+x+x2+…)=1 /(1-x),即:
v
v
v
vv
2
2
v
22 22
1e
e
1 1e
11
eeee
/T
/T
/T
/ T / T h / k T h / k T
q
x






以基态能级的能量为能量零点时的振动配分函数 qv0为
v0
v
02
vv 22
1
ee
ee
11
1e 1e
,/ k T h / k T
h / k T h / k T
/Th / k T
qq







当 hν<<kT时 0v kTq hv? 2e1
2!
x xx(利用公式 )
0v 1q?当 hν>>kT时
vvqf?
因一维谐振子的振动自由度数为 1,故
:一个振动自由度的配分函数
vf
e 0 e 1 e 2 e 0e e,0 e,1 e,2 e,0e e e e,,,,k T k T k T k Tq g g g g
只讨论粒子的电子运动全部处于基态,求和项自第二项起均可被忽略,则对大多数双原子分子
e0 1,g?
e00e e e 0e,kT,q q g =
常数
6.电子运动配分函数
qt,qr,qv,qe的数量级与 T,V的关系配分函数 数量级 ( T? 3 0 0 K )

32
t
2
2 π m k Tq
Vh
10
30
~ 10
32
2
r 2
8 π I k T
q
h?
10 ~ 10
2
2
v
e
1e
h v k T
h v k T
q
10
- 2
~ 10
0
e e,0
qg?
1 ~ 10
q0e取一项是常数,取两项以上是 T的函数;
qv,qr都是 Τ 的函数;
qt是 T和 V的函数,但
qt/V就只是 T的函数 。
qt>>qr,qv,qe。
7.核运动的配分函数只考虑核运动全部处于基态
n0
n n 0 e
,
kT
,qg

0
n n,0qg
常数独立子系统的热力学能
ii
i
Un
e i kTiiNngq

则,e ikT
ii
i
NUg
q

式中某一能级 εi的粒子数为 ni,由玻耳兹曼分布公式得到:
§ 9.6系统的热力学能与配分函数的关系
1.热力学能与配分函数的关系
22
e
11
ee
i
ii
kT
i
iV
V
ik T k T
i i i
ii
q
g
TT
gg
k T k T









移项,2 e ikT
ii
iV
qk T g
T


e
i
kT
i
i
qg


由于
22 ln
VV
N q qU k T N k T
q T T


则:
上式即独立子系统的热力学能与配分函数的关系式。
将配分函数的析因子性质代入,得
2 t r v e nln
V
q q q q qU N k T
T


式中仅 qt与 V有关,故
2 2 2tv r
22 e n
ln d lnd ln
dd
d ln d ln
dd
V
qq q
U N k T N k T N k T
T T T
q q
N k T N k T
TT




t r v e nU U U U U U
所以若以基态能量为起点 ( 各运动形式基态能量规定为零时 ),系统热力学能 U0用以上同样的方法可导出 。
0
0 ln
V
qU N k T
T


22t r
tr
22ve
ve
2 n
n
ln dln
d
dl n dl n
dd
dln
d
V
q q
U N k T U N k T
TT
qq
U N k T U N k T
TT
q
U N k T
T




上式说明系统的热力学能与能量零点的选择有关。
式中,Nε0是系统中全部粒子均处于基态时的能量,可以认为是系统于 0k时的热力学能 U0,则 U0=U-U0。
同样 0 0 0 0 0 0
t r v e nU U U U U U
0 0U U N00 e kTqq
因 则:
而且有:
00
t t r r
0
vv
00
en
2
0 0
U U U U
N h v
UU
UU


电子及核运动处于基态
(1) Ut0的计算 02 t
tt
3
2
2
2
ln
2
ln
V
V
q
U U N kT
T
m kT
V
h
N kT
T










2.Ut0,Ur0及 Uv0的计算当系统物质量为 1mol,即 N=1molL,摩尔平动热力学能为,又因平动自由度数为 3,故每个平动自由度的摩尔能量为 1/2 RT。
0t 32U / R T?
0
t
3
2U N k T?

3 3
1
2 2
2
2 2
3
2
2
2 32 3
ln
2 2
2
V
m kT mk
V VT
h h
TT
m kT
V
h











双原子等线型分子的转动自由度为 2,所以
1mol物质每个转动自由度对热力学能的贡献同样是 。1
2 RT
对 1mol物质 0
rU R T?
( 2) Ur0的计算

0 2 2rr
rr
d ln
d ln
dd
T
q n !
U U N k T N k T N k T
T T r ! n r !



v
vv
v
v
vv
0
0 2 2v
v
2
2
v
2
1
vv
1
d ln
d ln
1e
dd
1 e e
1e
e1
1 e e 1
T
TT
T
/T
/ T / T
q
U NkT NkT
TT
T
NkT
Nk Nk


















通常情况下 Θv/T>>1,Uv0≈0 (振动能级不充分开放 )
(3)Uv0 的计算
0vm 0 m72,U R T U U R T② 振动能级充分开放,
0vm 0 m5 02,U R T U U① 振动能级不充分开放:
对 1mol气体,Uv0=RT
当 Θv/T<<1 (振动能级充分开放 )
0
vv
v
1
11
U N k N k T
T


综上所述,在粒子的电子运动与核运动均处于基态时,
单原子气体的摩尔热力学能为,( 转动及振动运动均可不予考虑 )
m t e n 0 m
3
2,U U U U R T U
对双原子气体:
t r v e nU U U U U U
m
ln


2
V,
V V
qC R T
TT

1m o lNL且已知:
V
ln2 qU N k T
T

mm

V,V
UC
T
§ 9.7系统的摩尔定容热容与配分函数的关系
1.摩尔定容热容与配分函数的关系又 且 ε0是定值e 00 kTqq
m
ln


0
2
V,
V V
qC RT
TT
对比上两式可知物质的 Cv,m不受能量零点选择的影响。
当电子运动和核运动均处于基态时,将 q0的析因子性质代入
t r v0 0 0 0q q q q
tv r
m
l n d l nd d l n d
d d d d



00 0
2 2 2
V,
V V
qq qC RT RT RT
T T T T T T
表明:物质的摩尔定容热容是 1mol物质平动,转动,
振动三种独立运动对热容贡献之和 。

t
t
r
r
v
v
ln
d dl n
dd
dl nd
dd













0
2
0
2
0
2
V,
V V
V,
V,
q
C RT
TT
q
C RT
TT
q
C RT
TT
m t r vV,V,V,V,C C C C

2.CV,t,CV,r,CV,v的计算
t?
3
2V,CR
r?V,CR
双原子分子等线型分子,能级充分开放?
vv
v
v ee



22
1TTV,CR
T
当 >>T,即 >>1 (通常情况下 )v? v
T
v? 0V,C

<<T,即 <<1
(能级得到充分开放时 )
v? v
T
v?V,CR

mt
3
2V,V,C C R
对单原子分子:不必考虑粒子的转动和振动
m t r v
5
2V,V,V,V,C C C C R
对双原子分子:平动、转动能级充分开放,
振动基本处于基态时
§ 9.8 系统的熵与配分函数的关系系统的 N,U,V确定后,各状态函数均已确定,所以
( ) = (,,)S S N,U,V N U V
1 1 1 1 2 2 2 2
1 1 1 1 2 2 2 2
1 1 1 1 2 2 2 2
( ) = ( ) + ( )
( ) = ( ) ( )
ln ( ) = ln ( ) + ln ( )
S N,U,V S N,U,V S N,U,V
N,U,V N,U,V N,U,V
N,U,V N,U,V N,U,V


则,S 与 Ω 应存在一定的关系。因
1.玻耳兹曼熵定理此式即玻耳兹曼熵定理,它表明了独立子系统熵 S
与系统总微观状态数 Ω的关系。
lnS则:
lnSk
利用理想气体的平衡分布是玻耳兹曼分布,
可得引入比例常数
lnSc
2.摘取最大项原理若以 WB代表最可几分布的微态数,根据摘取最大项原理,当粒子数 N很大时,,于是得
ln lnBW
ln BS k W?
由玻耳兹曼熵定理表明,隔离系统的熵说明其总微态数的多少,这就是熵的统计意义。 Ω↗,S↗ 。
0K时纯物质完美晶体中粒子均处于基态,Ω=1,S0=0
3.熵的统计意义离域子系统在 N,U,V确定的条件下,当 ni<<gi时,
l n l n l n !!
in
i
B B i i i
ii i
gW W ( n g n )
n
由斯特林近似式,
及:玻耳兹曼分布式 代入上式
NNNN ln!ln
e i kTiin ( N / q ) g?
l n l n l n
l n l n l n
l n l n
B i i i i i
i
ii
i i i i i i
i
ii
ii
i
W n g n n n
nN
n g n n g n
q kT
nq q U
n n N N
N kT N kT







4.熵与配分函数的关系同理,对定域子系统,! ini
B
ii
gWN
n!


又 ln
BS k W?
l n l n
00q U q U
S N k N k N k N kN T N T则
l n l n
0
0UUS N k q N k q
TT
则:
由上两式可知,系统的熵与能量零点选择无关。
以离域子系统为例,以上各项为:
t t t
r r r
v v v
e e e
n n n
ln
ln
ln
ln
ln





00
00
00
00
00
S N k ( q / N ) U / T N k
S N k q U / T
S N k q U / T
S N k q U / T
S N k q U / T
由 q=qtqrqvqeqn和 U=Ut+Ur+Uv+Ue+Un可推导得:系统的熵是粒子各种独立运动形式对熵的贡献之和 。
t r v e nS S S S S S
5.统计熵的计算由于目前还不能求得电子运动和核运动对熵的贡献,而且在一般的物理化学过程中其值不变,故通常把统计热力学方法计算出系统的 St,Sr,Sv之和称为统计熵,以 S表示 。
t r vS S S S
对理想气体,每摩尔粒子数 Nmol-1=L,m=M/L,
V=nRT/p,n=1mol,则,
-1m,t l n l n l nk g m o l35 2 0 7 2 322
a
M T PS R,kp
上式称为 萨克尔 — 泰特洛德 (Sackur-Tetrode)方程,
是计算理想气体摩尔平动熵常用公式。
(1) St的计算
t t t
32
00
3
2 32/( m k T )q q V U N k T /
h
tt
t l n l n
00 32
3
25
2
/qU ( m k T ) V N k
S N k N k N kN T N h?
上式适用于通常转动能级充分开放的线性分子 。
m,r
r
ln TS R R


( 2) Sr的计算
00r r r r q q T / U N k T
00t t t rl n l nS N k q U / T N k ( T / ) N k
(3)Sv的计算
00t t t rl n l nS N k q U / T N k ( T / ) N k
vv
11
0 0 1
v v v vl n l n 1 e e 1
TTS N k q U / T N k N k T




vv0 1 0 1v v v1 e e 1/ T / Tq ( ) U N k ( )
某些物质 298.15K的 Sm?统计熵与 Sm?量热熵统计熵与量热熵的差别称为 残余熵 。
残余熵产生的原因:低温下量热实验中系统未能达到真正的平衡态。
146.6
205.14
186.3
206.59
223.07
146.34
205.15
186.88
206.80
223.16
Ne
O2
HCl
HI
Cl2
统计熵 /J·mol-1·K-1物质 量热熵 /J·mol-1·K-1
6.统计熵与量热熵的简单比较
(1)亥姆兹霍函数与配分函数的关系将 U和 S代入亥姆霍兹函数的定义式 A=U-TS,可得:
ln

Nq
A k T N!
(离域子系统 )
ln NA k T q (定域子系统 )
§ 9.9其它热力学函数与配分函数的关系
1.A,G,H与配分函数的关系
lnlnN TqqG k T N k T VN ! V(离域子系统 )
lnlnN TqG kT q NkT V V(定域子系统 )
(2)吉布斯函数与配分函数的关系
ln TT qAp NkT G A pVVV由 得
(3) 焓与配分函数的关系将 p与 U代入焓的定义式 H=U+pV,得:
l n l n2
T T
qqH N k T N k T V
VV


2.理想气体的标准摩尔吉布斯函数 Gm,TΘ
理想气体的标准摩尔吉布斯函数是在标准压力
pΘ =100 kPa下,温度为 T时理想气体的标准摩尔吉布斯函数:
tl n l 1TTq / V ) n q / V ) V( (
lnln
l n l n
ln
N
TT
qq
G kT N kTV
N ! V
N kT q N kT N N kT N kT
N kT ( q / N )





=-
式中 q0为以基态能量为起点 1mol某理想气体在标准压力及温度 T下的配分函数,U0,m是 1mol理想气体温度降至 0K时的热力学能 。
m
0
m m 0 m
ln
ln
ln ln
T
,T T
,T,T,
N n L G n L k T ( q / N )
G G / n L k T ( q / N )
G R T ( q / N ) G R T ( q / N ) U+



$$
3.理想气体的标准摩尔吉布斯自由能函数
mm ln
0
0,T,GU qR
TN
(Gm,T?-U0,m)/T称为标准摩尔吉布斯自由能函数,它是温度的函数 。 由于 0K时 U0,m H0,m,
故它也可用 (Gm,T?-H0,m)/T表示 。
与标准摩尔吉布斯自由能函数类似,定义标准焓函数为 ( Hm,T?-U0,m)/T。 2
m
2
ln ln
ln
VTT
,T
V
qq
N kT N kT V
H TV
H
nn
q
RT RT
T








4.理想气体的标准摩尔焓函数
0
2
m 0 m
0
m 0 m
ln
ln
,T,
V
,T,
V
q
H RT U RT
T
HU q
RT R
TT






又因 U0,m≈H0,m,故标准摩尔焓函数可近似表示为 。它也是温度的函数。?
mm? 0,T,H H T
一些热力学数据表中给出了某些理想气体在不同温度下的 值及 值或,应用这些数据及气体在某温度下的标准摩尔生成焓,可以计算理想气体化学反应在某温度下的 值,进而可以计算该温度下反应的标准平衡常数 。
mm? 0,T,G U Tmm? 0,T,H U T
mm? 0,T,HU
rm? G
K
§ 9.10 理想气体反应的标准平衡常数对理气间任一反应
B
B
0B
r m B m,B
B
lnG ( T ) R T K G

B m,B 0 m,B B 0 m,B
BB
ln,,R T K G U U

m,B m,B
B B m,B
BB
mm
r r m
ln






0
0
0
0
11
11
,
,
,
,
GU
KU
R T R T
GU
U
R T R T
式中 Δr(Gm?-U0,m)/T 为反应的标准摩尔吉布斯自由能函数变,可查表。 是单位反应进度在 0K
时的热力学能变化 。
m? 0,U
r 0 m B 0 m,B
B
,,UU

r 0 m r m 2 9 8 K r m 2 9 8 K 0 m,,,,U H H U
m,B 0 m,Bm 0 m
rB
B
,,GUGU
TT?

r m 2 9 8 K 0 m B m 2 9 8 K 0 m,B
B
,,,,H U H U
式中:
§ 9.11 系综理论简介力学量的长时间平均与系综平均如,N,V,T一定的系统,其能量 E
宏观,一定微观:有波动( E1,E2,E3,?)
系统的力学量的观测值都是在给定条件下和观测时间内系统所经历的一切微观态的力学性质的统计平均。
则平衡态的 E是系统随时间变动所经历的系统态的能量 Ei的统计平均值 。
,1 1,2 2
,1 1,2 2
NN
N
tt
t E E
EEE P E P E
式中? 很长时间内,此系统经历的微观状态数 (显然是很大的数 ),( 注意 与 N的区别 ) 。
而 Nt,1,Nt,2,Nt,3? 其中出现在微观状态 (
即系统态 )1,2,3,… 的次数 。
N
N
N
N
N
N 2,
2,
1,
1,
t
E
t
E PP,
分别为出现系统态 1,2,… 的概率个与所研究的具体系统的 N,V,T相同的系统组成的系综:
N
… 从微观状态看,个系统中分别在系统态 1,2,3,… 的个数; E1,E2,E3,… 各系统态 1,2,3,… 的能量。
,,,321 NNN N
这称为力学量 E的 系综平均 。
则这 个系统能量的平均值是
N
NN 2S,211,S EEE
N
⑵ 系综的分类
(ii)微正则系综 —— N个隔离系统 (N,V,E一定 )
的集合
(i)正则系综 —— N个封闭的定温系统 (N,V,T一定 )的集合
(iii)巨正则系综 —— N个敞开的定温系统 (μ,V,T
一定 )的集合系综? 构想的具有与所研究的具体系统同样条件的 (如
N,V,E一定,等 )大量的 (N个 )系统的集合 。
⑴ 系综的定义系综的定义及分类
(ii):孤立系统 (即 N,V,E一定 )的全部可能达到的量子态是等概率的。用于封闭系统 (即 N,V,T一定
)时可表述为:能量相同的各量子态有相同的概率。
统计热力学基本假设
(i):当 N→∞ 的极限情况下,在实际系统中任何力学量的长时间平均值 〈 Gt〉 等于该力学量的系综平均值 〈 G〉,只要统计系综和实际系统的热力学状态及环境完全相同。
N,V
刚性绝缘壁刚性导热壁
N个粒子数为 N,体积为 V
的系统组成的,超,隔离系统 (系综 )