1、Statistical Thermodynamics and Chemical Kinetics,State Key Laboratory for Physical Chemistry of Solid Surfaces,厦门大学固体表面物理化学国家重点实验室,Lecture 3,4.1 理想气体4.1.1 化学位及其标准态 设有给定T、V、N条件下的理想气体，其分子的化学位当为：,Chapter 4 简单体系,（不考虑粒子间相互作用）,表示指定分子基态能量为零的配分函数。,（4.1）,（4.2）,若已知0、 1、 2为分子各能级能量，则分子配分函数为：,再定义分子的标准配分函数为：,其中qi
2、nt,0 为分子(基态能量为零)的内配分函数；,（4.4）,理想气体分子的平动配分函数由三维平动子公式：,又V=NkT/p, 故可得,（4.3）,分子的内配分函数：,即为标准状态下，1mol气体中一个分子的配分函数。,分子基态能量为零的配分函数可表示为：,由于基态0绝对值未知，实际上o的绝对值仍不可得。,因热力学中某组分化学位以j/mol为单位，,j/mol为单位,其中,（4.5）,（4.6）,（4.7）,因而分子的化学位可表示为：,4.1.2 气体的热容,利用配分函数的析因子性质，体系热容可分解为各单一运动形态的贡献之和。则有平动、转动和振动等三种运动形式的热容，分别为：,理想气体热容的统计
3、表达式如下：,（4.8）,（4.9）,（4.10）,例如，对双原子分子而言，其摩尔热容当为：,4.1.3 多组份气体的 及,： 设气体中仅含A、B两组份，分子数分别为NA和NB, 则其正则配分函数可定义为：,（4.11）,对更多组分气体有：,（4.12）,（4.13）,Ei(A,B) 等于全部A分子和B分子各自单独存在时的能量之和:,Ei(A,B): 体系量子态i的能量,（2） : 对双组份理想气体的开放体系有,（4.14）,（4.15）,同理，多组分理想气体的开放体系有,4.1.4 混合熵,不同种类理想气体等温混合，体系内能不发生变化，但体系总熵增加，其净增量即为过程混合熵Smix:,（4.
4、17）,（4.16）,（4.18）,回顾玻尔兹曼公式，体系熵的增加归因于其微观状态数的增加。因而，理想气体混合过程体系内部必然发生了的变化。,等温混合前，气体A、B的正则配分函数分别为：,等温混合后, 混合气体总体积V=VA+VB, 则有：,（4.21）,（4.19）,（4.20）,由道尔顿分压定律当知：,(混合后）,(混合前）,离域子体系等温混合熵产生的唯一根源： 不同组分气体混合后，分子活动的空间范围都比原来纯态时扩大了。微观上，体系空间扩大将导致分子平动能的变小及能级间距变窄，在等温混合时为维护体系内能不变，必然导致能级分布数的变化，其结果是体系微观状态数增多，总熵随之增加。,代入前式，
5、并令N = NA+NB，即有,如为多组分气体理想混合，则有,（4.22）,（4.23）,思考题：当A、B为同种气体等温混合时，重复上述推导过程，证明混合熵为零！,4.1.5 吉布斯佯谬,设使上述混合熵推导过程中的气体分子A、B为同类分子，那么合二为一之后体系总熵是否增加了呢？答案否定的，此过程并不导致熵变。但是在历史上有关此类过程的思索却出现了吉布斯佯谬。,（4.24）,若直接应用上式来计算混合熵，后果如何呢？,早期经典统计缺乏量子论概念，对离域子体系和定域子体系在微观态的实现方式上未能加以严格区别，以致在应用相空间求算体系配分函数时潜伏巨大缺陷。,例如，对单原子分子气体，其经典微正则配分函数
6、为,（4.26）,混合前后体系熵变为,A、B气体理想等温混合前后各体系的配分函数分别为,且等温混合过程满足：,（4.25）,依经典统计，体系熵为,与正则系综统计法所得结果一致。,不论A、B是否同类分子，上述推算所得混合熵总大于0。与熵的加和性相抵触。,（4.26）,这一似是而非的命题是Gibbs提出的，故称Gibbs佯谬。,其根源在于经典微正则配分函数表达式中未能体现微观粒子不可分辨性的修正。,考虑不可分辨性后，有,（4.27）,（4.28）,当A、B为同类分子时，必有：,借助斯特林公式对上式展开最终就有,思考题：若A、B粒子不同，试重复上述推导过程，得到混合熵的表达式是否与式（4.26）相同
7、？,4.2 金属晶体,4.2.1 晶体的Einstein Model. (理想定域子体系),一维谐振子振动配分函数,一维谐振子能级：,（4.29）,金属中N个原子的振动各自独立；原子的每一振动自由度都服从虎克定律；所有的振动自由度都有相同的特征频率 () 。,要确定金属晶体的正则或巨正则配分函数，首要问题是确定晶体中原子振动的频率分布，迄今这仍是一个相当复杂的问题。先后提出了各种各样的理论模型。,据此，金属的正则配分函数为,由此可进一步推得Einstein晶体的主要热力学性质：,所有原子振动基态的能量和,（4.30）,（4.31）,（4.32）,ii）温度趋近0K时，,一般定义E=h/k, 为
8、振动的Einstein特征温度。,i）在足够高的温度下，,即所有原子均处振动基态。,当T0K时，S 0。因为此时所有原子均处于振动基态，与其对应的微观态只有一个即=1.,（4.33）,（4.34）,（4.35）,如从巨正则系综出发，则有,补充： 一维谐振子,一维谐振子能级：,其配分函数：,（4.36）,（4.37）,（4.38）,（4.39）,与式（4.30）相比，仅以N 代替 N。,4.2.2 金属热容,根据Einstein模型，体系内能的正则平均为：,（4.40）,（4.41）,则摩尔热容为：,则热容平均值为,且,（4.42）,爱因斯坦晶体热容曲线,1,2,3,在足够高温下，金属的热容趋近
9、于3R, 与Dulong-Petit经验定律一致。在T 0K时，热容趋近于0, 大体能反映实验事实。严格地说，在较低温度范围，由上式计算的热容仍与实际情形有一定偏差，在于Einstein model过于理想化。,4.2.3 晶体与其蒸汽的平衡（气固平衡）,设晶体A(s)与其蒸汽A(g)达成平衡： A(s) A(g)其平衡条件当为：(s)= (g),（4.43）,（4.44）,若A(s)为Einstein晶体，则,如A(g)为理想气体，有,（4.39）,连立后即有,（4.44）,（4.45）,0= 0(g) - 0(s)为气态A与固态A的原子基态能量之差。,形式上，上式与克拉贝龙公式 lnp =
10、 -H/RT +B 十分相似，且有,若以摩尔数计量，则成,为晶体在0K下的升华能。,两相平衡为高温状态，近似有,4.3 理想固溶体,固溶体是一类多组分混合物，其主要特征是不同组分的粒子共同构成一个晶格。在晶格的点阵点上，任一组分粒子的排布均为随机的。形成理想固溶体的基本条件是：1) 相互混合的组分的晶格同型，且粒子大小尺寸相近；2) 任意交换不同点阵位置的两个粒子，体系能量不变（或几乎不变）。,空间构型方式数（或称该能级的简并度）依排列组合公式，有,（4.46）,（4.47）,（4.48）,4.3.1 二元固溶体的正则配分函数当A、B两种组分(粒子数分别为NA、NB)形成固溶体时，其微观态的实
11、现方式除了粒子的能级(或量子态) 分布外，还包括异类粒子在同一晶格点阵上的空间构型分布。,体系某一能级i的能量可表示为：Ei(A,B) = Ei(A) + Ei(B),（4.49）,（4.50）,（4.51）,（4.52）,定义构型配分函数,即理想固溶体，混合过程内能变化为零。,二元固溶体的某些状态性质平均：,式中k*=-kT lnqk为固溶体中k组分粒子的标准化学位。,（4.53）,（4.54）,（4.55）,（4.56）,故体系混合自由能为：,4.3.2 固溶体混合熵,其结果表观上与离域子体系的混合熵公式一样； 从推导过程看，两类混合体系的混合熵增加的起因均为混合后体系微观态数增加所致，但
12、微观态变化的图象却完全不同-定域子体系: 源自于排列方式数增加 离域子体系: 源自于粒子运动空间增大所致的平动能级变窄,4.3.3 固溶体的巨正则配分函数,对二元固溶体开放体系，其巨正则配分函数可导出为：,（4.57）,借助下列二项式展开,二元固溶体状态性质的巨正则平均(1) 粒子数:,（4.58）,（4.59）,（4.60）,(2) U:,（4.61）,（4.62）,(3) S: 依熵的巨正则平均定义当得：,效仿撷取最大项原理, 因,因而有,（4.63）,4.4 气固吸附,气体在固体表面的吸附即气体分子被定域在固体表面的吸附位。倘若发生的吸附是单分子层的，则可表示如下：,基本假设： 固体表面
13、上吸附分子的铺展是单分子层的；吸附分子间不存在任何相互作用力。 分子吸附或脱附因之是完全随机的！,4.4.1单层吸附的理想模型与Langmuir吸附方程,系综统计法研究：设固体表面有NS个全同吸附位，气固两相平衡时，有NA个吸附分子被吸附在固体表面，每一分子占据一个吸附位。,（4.64）,指定固体表面原子的配分函数qs为1，qA(S)为吸附分子A的配分函数。,吸附相的正则配分函数可导出为：,吸附分子A的化学位,吸附相自由能为,0(S)为吸附分子A基态能量，=NA/NS为表面覆盖度。,（4.65）,（4.66）,平衡时两相化学位必相等：,与吸附相平衡的另一相是气相A(g)，假定系理想气体，其化学
14、位为,故可得0= 0(g) - 0(S)为气态A与吸附态A的基态能量之差。,（4.67）,思考题：请由巨正则系综导出上式。,这实际上是Langmuir方程的变体，可令,4.4.2 吸附熵,据熵的统计表达式，对吸附相,（4.68）,暂确认吸附分子只有振动运动，故前一式称为吸附分子的振动熵, 而后一式称为吸附分子的构型熵。,若定义,对构型熵展开：,（4.69）,故脱附过程的熵变为：,而气相的熵为,前者即为通常定义的吸附熵，实质上是吸附过程中的热熵。后者(构型熵)与温度无关，即使T0K 仍然存在，因而无法直接通过量热法测定。,4.4.3 多层吸附 -B.E.T公式,在吸附剂表面上，吸附质亦可发生多分
15、子层吸附。,独立吸附位模型,Hill于1960年提出了独立吸附位模型 ，并应用巨正则系综导出多层吸附方程。该方法具有普遍性，可用于处理各种不同类型的吸附体系。,任何一吸附位都可有限制或无限制地沿着垂直于表面的方向叠加吸附分子。或者说每个吸附分子本身构成一个吸附位。但从第二层开始，吸附力性质转属物理吸附，故多吸附层的叠加类似于气体分子的凝聚液化。,以NS 和NA表示固体表面吸附位数和吸附分子总数；am表示吸附相中已形成m个吸附分子层的吸附位数；若将任一吸附位连同叠加在它上面的m个吸附分子看成是构成固体表面的一个独立单元，其配分函数q(m)由析因子性质有：,独立吸附位模型,其中q(0)=1为空白吸
16、附位的配分函数，q1为第一层吸附分子的配分函数，q2为高层吸附分子的配分函数。,在理想情况下，固体表面各个吸附位互不干扰，任何吸附位如何叠加均为随机。,设吸附相中各独立吸附位分布的一套组合样式am为：吸附层数 m: 0, 1, 2, , l 吸附位数 a: a0, a1, a2, al 任何这样的一套分布必须满足：,（4.70）,（4.71）,以上X表示一切可能组合样式;,当NA给定不变时，该多层吸附相的正则配分函数为,X表示某一套分布amX在NS个吸附位上可能实现的构型方式数：,令E(m)K为指定量子态i下，第K个独立吸附位上m个吸附分子的能量和，即有,（4.72）,第K个吸附位(层数为m)
17、的配分函数：所有可及量子态的波尔兹曼因子求和！,则多层吸附体系(T、NS、NA给定)的正则配分函数为,吸附过程中气相与吸附相时刻在发生分子的交换，因此确立吸附相的巨正则配分函数更具实际意义。 可得：,lNs 吸附体系吸附分子总数的极大值,l 单个吸附位上吸附分子层数的极大值,上式可借助多项式展开来进一步简化！,（4.73）,A(S) 吸附分子的因子 (假定所有吸附层的分子化学位相同!),X 已包含NA的一切可能变化，且恒有,借助多项试展开，式(4.71)可表示为：,（4.74）,（4.75）,显然，代表了单个独立吸附位的巨正则配分函数。,1) 在单层吸附情况下，m的最大值为1，且有,令,（4.
18、76）,2) 对多层吸附，若吸附分子层数(l)不受限制，则有：,代入前式即有：,（4.77）,（4.78）,注：平衡压强p、饱和蒸气压p0、平衡气体吸附量V 、单层饱和吸附量Vm,知x1, 有,故得此即多层吸附B.E.T公式。,与通常表达式比较,4.5 化学平衡常数,据热力学原理，反应达平衡时，满足,(4.6),(4.79),(4.80),(4.81),qo0为标准状态下1mol气体中分子配分函数。,4.5.1 Kp的统计表达式 设有气相反应: aA + bB cC + dD,(4.82),产物与反应物的分子基态能量差。,公共能标原点,反应物分子基态,产物分子基态,气态原子基态,D(R),D(
19、P),实验可测,量子化学计算,运用统计力学研究化学平衡时，无需获得单个分子基态能量的绝对值！,亦可推出以浓度或粒子数表示的平衡常数(KC或KN)关系式:,(4.83),(4.84),为粒子单位体积的配分函数,4.5.2 化学反应的涨落,化学反应中当然会存在涨落现象，各组分的涨落程度明显取决于该组分在平衡状态下的平均粒子数。,(4.85),暂考虑气相反应： A B （单分子异构反应）采用巨正则系综，可合理假定反应体系总分子数N固定不变。,则体系巨正则配分函数为：,(4.86),(4.87),(4.88),(4.89),(4.90),因平衡时有 A= B或 A= B ，于是有,而上述反应A、B分子
20、数的散差分别为,则反应物分子与产物分子的相对涨落分别为,可见KN越大，产物分子相对于反应物分子的涨落就越小。,二者的相对涨落比为,4.5.3 同位素交换反应,属化学反应中的一类特殊反应，可利用同位素交换来富集同位素。此类反应的几个典型示例有：H2 + D2 2HD; H2 + 2DI 2HI + D2; (3) O216 + 2H2O18 O218 + 2H2O16,标准态下分子配分函数（基态能量为零),(4.91),0= 20(HD) - 0(H2) - 0(D2),其显著特点为：反应物的分子数总是与产物的分子数相等，且两者的分子又存在一一对应关系。,例如反应(1), 由平衡常数的定义有：,
21、同位素交换反应中，同一交换分子对的两个分子，其核质量数虽略有不同，但核电荷及核外电子云分布实际上无甚差别，二者的基态能量差可近似地归结为分子的振动基态之差:,(4.92),(4.93),(4.94),反应物及产物均为电子数相同的双原子分子，暂不计较分子中电子及核自旋运动，故上式中q00 仅含平动、转动和振动三种运动形态的配分函数项，即,双原子分子的解离能可由光谱数据获得。,一般地，反应的产物分子与反应物分子的基态能量之差0可看作是反应物分子与产物分子的解离能之差D，即,分子平动配分函数： (标准态),双原子分子转动配分函数（刚性哑铃转子近似）：,双原子分子振动配分函数（谐振子近似）：,对称因子
22、：,双原子分子转动惯量：,（转动基态能量为0）,（平动基态能量为0）,利用分子的平动、转动、振动配分函数，可得:,(4.95),(4.96),(4.97),又知0=657 kjmol-1(实验值)。并将分子的m、I代入得：,若交换反应的温度比分子振动的特征温度低得多，即 T v=h/k, qv,01。 即有,H2 + D2 2HD反应的Kp值：实验 vs.计算,显然对称因子Re对同位素交换反应的KP起关键作用：,在量子力学中，同核双原子分子的振动基态能量0=h/2, 同位素交换反应的分子基态能量差主要源于分子的振动基态能量差(电子及核自旋运动基态能差)，这就是零点能效应。量子力学结合统计力学研究化学反应时，必须注意零点能！,
