1、1,化学原理第九章 统计热力学简介,2,主 线,等同性修正,非定位体系,3,第九章 统计热力学简介,9.1 概论,9.5 各配分函数的求法及其对热力学函数的贡献,*9.3 Bose-Einstein统计和Fermi-Dirac统计,9.4 配分函数,9.2 Boltzmann 统计,*9.6 晶体的热容问题,9.7 分子的全配分函数,9.8 用配分函数计算 和反应的平衡常数,4,9.1 概 论,统计热力学的研究方法和目的,统计体系的分类,统计热力学的基本假定,5,统计热力学的研究方法和目的,物质的宏观性质本质上是微观粒子不停地运 动的客观反应。虽然每个粒子都遵守力学定律， 但是无法用力学中的微
2、分方程去描述整个体系的 运动状态，所以必须用统计学的方法。,根据统计单位的力学性质（例如速度、动量、位置、振动、转动等），经过统计平均推求体系的热力学性质，将体系的微观性质与宏观性质联系起来，这就是统计热力学的研究方法。,6,统计热力学的基本任务,根据对物质结构的某些基本假定，以及实验所得的光谱数据，求得物质结构的一些基本常数，如核间距、键角、振动频率等。,利用这些数据可以计算分子配分函数。再根据配分函数求出物质的热力学性质，这就是统计热力学的基本任务。,7,统计热力学的基本任务,该方法的局限性：计算时必须假定结构的模型，而人们对物质结构的认识也在不断深化，这势必引入一定的近似性。另外，对大的
3、复杂分子以及凝聚体系，计算尚有困难。,该方法的优点： 将体系的微观性质与宏观性质联系起来，对于简单分子计算结果常是令人满意的。不需要进行复杂的低温量热实验，就能求得相当准确的熵值。,8,统计体系的分类,目前，统计主要有三种：,一种是Maxwell-Boltzmann统计，通常称为Boltzmann统计。,1900年Planck提出了量子论，引入了能量量子化的概念，发展成为初期的量子统计。,在这时期中，Boltzmann有很多贡献，开始是用经典的统计方法，而后来又有发展，加以改进，形成了目前的Boltzmann统计。,9,统计体系的分类,1924年以后有了量子力学，使统计力学中力学的基础发生改变
4、，随之统计的方法也有改进，从而形成了Bose-Einstein统计和Fermi-Dirac统计，分别适用于不同体系。,但这两种统计在一定条件下通过适当的近似，可与Boltzmann统计得到相同结果。,10,统计体系的分类,定位体系（localized system）,定位体系又称为定域子体系，这种体系中的粒子彼此可以分辨。例如，在晶体中，粒子在固定的晶格位置上作振动，每个位置可以想象给予编号而加以区分，所以定位体系的微观态数是很大的。,根据统计单位是否可以分辨，把体系分为定位体系和非定位体系,11,统计体系的分类,非定位体系（non-localized system）,非定位体系又称为离域子体
5、系，基本粒子之间不可区分。例如，气体的分子，总是处于混乱运动之中，彼此无法分辨，所以气体是非定位体系，它的微观状态数在粒子数相同的情况下要比定位体系少得多。,12,统计体系的分类,根据统计单位之间有无相互作用，又可把统计体系分为近独立粒子体系和非独立粒子体系,独立粒子体系（assembly of independent particles）,独立粒子体系是本章主要的研究对象,粒子之间的相互作用非常微弱，因此可以忽略不计，所以独立粒子体系严格讲应称为近独立粒子体系。这种体系的总能量应等于各个粒子能量之和，即：,13,统计体系的分类,非独立粒子体系（assembly of interacting
6、particles）,非独立粒子体系又称为相依粒子体系，体系中粒子之间的相互作用不能忽略，体系的总能量除了包括各个粒子的能量之和外，还包括粒子之间的相互作用的位能，即：,非理想气体就是非独立粒子体系,14,统计热力学的基本假定,概率（probability）,热力学概率,指某一件事或某一种状态出现的机会大小,体系在一定的宏观状态下，可能出现的微观态的总数，通常用 表示。,15,统计热力学的基本假定,等概率假定,例如，某宏观体系的总微态数为 ，则每一种微观状态 P出现的数学概率都相等，即：,对于U, V 和 N 确定的某一宏观体系，任何一个可能出现的微观状态，都有相同的数学概率，所以这假定又称为
7、等概率原理。,16,9.2 Boltzmann 统计,定位体系的最概然分布,Boltzmann公式的讨论非定位体系的最概然分布,撷取最大项法及其原理,值的推导,Boltzmann公式的其他形式,17,定位体系的最概然分布,一个由 N 个可区分的独立粒子组成的宏观体系（U，V，N为定值），在量子化的能级上可以有多种不同的分配方式。,设其分配方式为：,18,定位体系的最概然分布,但无论哪一种分布方式，都必须满足如下两个条件,或,或,这种分布的微态数相当于将N个不同的球在两个限制条件下分成若干不同的堆，根据排列组合公式，有：,19,定位体系的最概然分布,这是一种分布，在满足这两个条件下，可以有各种不
8、同的分布，则总的微观状态数为：,20,设有n个项进行求和，每一项都取最大值，则有,每种分配的 值各不相同，但Boltzmann认为其中有一项的值最大，即 ，在粒子数足够多的宏观体系中，可以近似用 来代表所有的微观数，这就是最概然分布。,21,由于,所以,问题在于如何在两个限制条件下，找出一种合适的分布 ，才能使 有极大值，在数学上就是求条件极值的问题。即：,22,将上式取对数，并用Stirling公式展开,再用Lagrange乘因子法，求得最概然的分布为：,式中 和 是Lagrange乘因子法中引进的待定因子。,23,先求,值的推导,已知,所以,或,或,最概然分布公式中已消去了,24,值的推导
9、,已知,代入得,再求,25,值的推导,根据复合函数的性质,可以证明上式中的方括号等于零，故而得,26,值的推导,因为,所以,这就是Boltzmann最概然分布公式,27,值的推导,已知,所以,又因为,所以,这就是定位体系的熵和Helmholtz自由能的计算公式,28,Boltzmann 公式的讨论,简并度（degeneration）,量子力学中把能级可能有的微观状态数称为该能级的简并度，用符号 表示。简并度亦称为退化度或统计权重。,能量是量子化的，但每一个能级上可能有若干个不同的量子状态存在，反映在光谱上就是代表某一能级的谱线常常是由好几条非常接近的精细谱线所构成。,29,例如，气体分子平动能
10、的公式为：,式中 分别是在 轴方向的平动量子数,当,只有一种可能的状态，是非简并的，,30,例如，气体分子平动能的公式为：,当,可分别为：,体系具有三种可能的状态，是简并的,31,有简并度时定位体系的微态数,设有 N 个粒子的某定位体系的一种分布为：,先从N个分子中选出N1个粒子放在 能级上，有 种取法；,但 能级上有 个不同状态，每个分子在 能级上都有 种放法，所以共有 种放法；,32,这样将N1个粒子放在 能极上，共有 种微态数。依次类推，这种分配方式的微态数为：,33,由于分配方式很多，所以在U、V、N一定的条件下，所有的总微态数为：,求和的限制条件仍为：,再采用最概然分布概念，令：,用
11、 Stiring 公式和 Lagrange 乘因子法求条件极值，得到微态数为极大值时的分布方式 为：,34,与不考虑简并度的公式相比，只多了 项,35,非定位体系的最概然分布,非定位体系由于粒子不能区分，它在能级上分布的微态数一定少于定位体系，所以对定位体系微态数的计算式进行等同粒子的修正，即将计算公式除以 。,则非定位体系在U、V、N一定的条件下，所有的总微态数为：,36,同样采用最概然分布的概念，用Stiring公式和Lagrange乘因子法求条件极值，得到微态数为极大值时的分布方式 （非定位）为：,由此可见，定位体系与非定位体系，最概然的分布公式是相同的。,但熵和Helmholtz自由能
12、计算式差一些常数项，但在计算变化值时可以消去。,37,Boltzmann 公式的其他形式,（1）将 i 能级和 j 能级上粒子数进行比较，用最概然分布公式相比，消去相同项，得：,（2）在经典力学中不考虑简并度，则上式成为,38,设最低能级为,在 能级上的粒子数为 ，略去 标号，则上式可写作：,这公式使用方便，例如讨论压力在重力场中的分布，设各个高度温度相同，即得：,39,撷取最大项法及其原理,设为定位体系，其中一种分布方式的微态数为,取对数，得：,将上面两式相减，得：,40,在上式中，,若是最概然分布，t 有极大值,因是最概然分布，将t 换作,41,因为,引用级数公式,略去 及更高次项，,又因
13、,42,在一个等分为二的长方形盒子中，均匀分布时，,由于分子运动，发生1%偏离，即,即,这个数值很小，表示 是“尖锐的极大”,43,(2) 能否用最概然分布的微观状态数代替总的微观状态数?,在粒子数足够大时，设,可以用数学方法证明,若某一能态的粒子数处于的间隔为,则所有可能分布的微态数为,即,说明了最概然分布足以代表体系的一切分布。,44,9.4 配分函数,配分函数的定义,配分函数的分离,配分函数与热力学函数的关系,45,配分函数的定义,根据Boltzmann最概然分布公式（略去标号 ）,令分母的求和项为：,q 称为分子配分函数，或配分函数（partition function),配分函数是量
14、纲一的量，单位为1,求和项中 称为Boltzmann因子。,46,配分函数q是对体系中一个粒子的所有可能状态的Boltzmann因子求和，因此q又称为状态和。,配分函数q是属于一个粒子的，与其余粒子无关，故称为粒子的配分函数。,将q代入最概然分布公式，得：,q中的任何一项与q之比，等于分配在该能级上粒子的分数，q中任两项之比等于这两个能级上最概然分布的粒子数之比，这正是q被称为配分函数的由来。,47,配分函数与热力学函数的关系,先讨论粒子数为N的非定位体系的热力学函数,（1）Helmholz自由能A,48,配分函数与热力学函数的关系,（2）熵 S,或根据以前得到的熵的表达式直接得到下式：,49
15、,配分函数与热力学函数的关系,（3）热力学能U,或从 两个表达式一比较就可得上式。,50,配分函数与热力学函数的关系,（4）Gibbs自由能G,将 A，p 代入，得：,51,配分函数与热力学函数的关系,(5)焓H,(6)定容热容CV,根据以上各个表达式，只要知道配分函数，就能求出非定位体系的各热力学函数值。,52,定位体系配分函数与热力学函数的关系,根据非定位体系求配分函数与热力学函数关系相同的方法，得：,53,定位体系配分函数与热力学函数的关系,54,定位体系配分函数与热力学函数的关系,由上列公式可见，U，H 和CV的表达式在定位和非定位体系中是一样的；,而A，S 和 G的表达式中，定位体系
16、少了与 有关的常数项，而这些在计算函数的变化值时是可以消去的。,本章主要讨论非定位体系,55,配分函数的分离,一个分子的能量可以认为是由分子的整体运动能量即平动能，以及分子内部运动的能量之和。,分子内部的能量包括转动能( )、振动能( )、电子的能量( )和核运动能量( )，各能量可看作独立无关。,这几个能级的大小次序是：,56,配分函数的分离,平动能约为,分子的总能量等于各种能量之和，即：,电子和核的能量 则更高,转动能约为,振动能约为,57,配分函数的分离,各不同的能量有相应的简并度,当总能量为 时，总简并度等于各种能量简并度的乘积，即：,根据配分函数的定义将 和 的表达式代入，得,58,
17、从数学上可以证明，几个独立变数乘积之和等于各自求和的乘积，于是上式可写作：,59,比较定位体系和非定位体系Helmholtz自由能的表达式，,两者仅在平动项上差了,今后的问题是如何计算各种运动的贡献,60,9.5 各配分函数的求法及其对热力学函数的贡献,原子核配分函数,电子配分函数,平动配分函数,转动配分函数,振动配分函数,单原子理想气体的热力学函数,61,原子核配分函数,式中 分别代表原子核在基态和第一激发态的能量， 分别代表相应能级的简并度。,62,由于化学反应中，核总是处于基态，另外基态与第一激发态之间的能级间隔很大，所以一般把方括号中第二项及以后的所有项都忽略不计，则：,如将核基态能级
18、能量选为零，则上式可简化为：,即原子核的配分函数等于基态的简并度，它来源于核的自旋作用。式中 sn 是核的自旋量子数。,63,对于多原子分子，核的总配分函数等于各原子的核配分函数的乘积,由于核自旋配分函数与温度、体积无关，所以对热力学能、焓和等容热容没有贡献。,但对熵、Helmholtz自由能和Gibbs自由能有相应的贡献。,从化学反应的角度看，一般忽略核自旋配分函数的贡献，仅在计算规定熵时会计算它的贡献。,64,电子配分函数,电子能级间隔也很大， 除F, Cl 少数元素外，方括号中第二项也可略去。虽然温度很高时，电子也可能被激发，但往往电子尚未激发，分子就分解了。所以通常电子总是处于基态，则
19、：,65,电子配分函数,若将 视为零，则,式中 j 是电子总的角动量量子数。电子绕核运动总动量矩也是量子化的，沿某一选定轴上的分量可能有 2j+1个取向。,某些自由原子和稳定离子的 是非简并的。如有一个未配对电子，可能有两种不同的自旋，如 它的,66,电子配分函数,电子配分函数对热力学函数的贡献为,67,平动配分函数,设质量为m的粒子在体积为 的立方体内运动，根据波动方程解得平动能表示式为：,式中h是普朗克常数， 分别是 轴上的平动量子数，其数值为 的正整数。,68,将 的表示式代入：,因为对所有量子数从 求和，包括了所有状态，所以公式中不出现 项。,在三个轴上的平动配分函数是类似的，只解其中
20、一个 ，其余类推。,69,因为 是一个很小的数值，所以求和号用积分号代替，得：,引用积分公式： 则上式得：,70,和 有相同的表示式，只是把a换成 b或 c，故,71,这就是Sackur-Tetrode公式，用来计算理想气体的平动熵,平动配分函数对热力学函数的贡献,72,对于1 mol 理想气体，Sackur-Tetrode 公式为,根据,根据热力学函数之间的关系，可以得到,73,根据热力学函数之间的关系，可以得到,代入相应的 表示式即得。,74,单原子理想气体热力学函数,由于单原子分子内部运动没有转动和振动，所以只有原子核、电子和外部的平动对热力学函数有贡献。,理想气体是非定位体系，所以它的
21、一系列热力学函数用配分函数的计算式分别分列如下：,（1）Helmholtz自由能 A,75,第一项是核和电子处于基态时的能量，第二项是与简并度有关的项。在计算热力学函数变量时，这些都可以消去。,76,这公式也称为Sachur-Tetrode公式。,（2）熵 S,可用来计算但原子理想气体的熵,77,因为 对热力学能没有贡献，只有平动能有贡献，所以：,（3）热力学能 U,（4）定容热容 CV,这个结论与经典的能量均分原理的结果是一致的，单原子分子只有三个平动自由度，每个自由度贡献 ，则N个粒子共有,78,对于理想气体， ，代入 A 的表示式，得：,（5）化学势,对1 mol气体分子而言，各项均乘以
22、阿伏伽德罗常数 , , 则1 mol气体化学势为,79,（5）化学势,当处于标准态时， ，则：,从该式可看出， 一定时， 只是T的函数。两式相减得：,80,将A的表示式代入，由于其它项均与体积无关，只有平动项中有一项与V有关，代入即得理想气体状态方程。,用统计热力学的方法可以导出理想气体状态方程，这是经典热力学无法办到的。,（6）状态方程式,81,转动配分函数,单原子分子的转动配分函数等于零，异核双原子分子、同核双原子分子和线性多原子分子的 有类似的形式，而非线性多原子分子的 表示式较为复杂。,（1）异核双原子分子的 ，设其为刚性转子绕质心转动，能级公式为：,式中J是转动能级量子数，I是转动惯
23、量，设双原子质量分别为 ，r为核间距，则：,82,转动配分函数,转动角动量在空间取向也是量子化的，所以能级简并度为：,称为转动特征温度，因等式右边项具有温度的量纲。将 代入 表达式，得：,83,令：,代入上式后，得：,从转动惯量 I 求 得,除H2外，大多数分子的 很小,在常温下， ，因此用积分号代替求和号,84,对于转动特征温度较高的分子，应该使用下式,85,对于同核双原子和线性多原子分子，还要除以对称数,对于非线性多原子分子，转动配分函数为,分别为三个轴上的转动惯量。,86,振动配分函数,双原子分子的,设分子作只有一种频率 的简谐振动，振动是非简并的， ，其振动能为：,式中v为振动量子数，
24、当v=0时， 称为零点振动能,87,令,称为振动特征温度，也具有温度量纲，则上式为,88,振动特征温度是物质的重要性质之一，Qv 越高，处于激发态的百分数越小，qV 表示式中第二项及其以后项可略去不计。,也有的分子 QV 较低，如碘的QV310K ，则第一激发态项就不能忽略。,在低温时， ，则 ，引用数学近似公式：,89,则 的表示式为：,将零点振动能视为零, 即 则：,90,多原子分子振动自由度,多原子分子的,为平动自由度, 为转动自由度,因此，线型多原子分子的 为：,n 为分子中原子总数,91,非线型多原子分子的 为：,92,a,b,x,z,y,图 7.1 双原子分子在空间中的取向,93,
25、9.7 分子的全配分函数,根据配分函数的定义及可分离的性质，分子的全配分函数应该由5个部分组成，即：,94,9.7 分子的全配分函数,对于单原子分子,对于双原子分子,95,9.7 分子的全配分函数,对于线型多原子分子,96,9.7 分子的全配分函数,对于非线型多原子分子,97,9.8 用配分函数计算 和反应的平衡常数,化学平衡体系的公共能量标度,从自由能函数计算平衡常数,热函函数,从配分函数求平衡常数,98,化学平衡体系的公共能量标度,粒子的能量零点,对于同一物质粒子的能量零点，无论怎样选取，都不会影响其能量变化值的求算。通常粒子的能量零点是这样规定的：,当转动和振动量子数都等于零时 的能级定
26、为能量坐标原点，这时粒子的能量等于零。,公共能量标度,化学平衡体系中有多种物质，而各物质的能量零点又各不相同，所以要定义一个公共零点,99,通常选取0 K作为最低能级，从粒子的能量零点到公共零点的能量差为,粒子的能量零点和 公共能量零点的关系,100,按公共的能量零点计算的分子能量为,按公共能量标度计算的配分函数为,101,按公共能量零点用非定位体系的配分函数计算的热力学函数的表示式为,式中,102,103,采用公共零点后， A，G，H，U的配分函数表达式中多了 项,而 和 p 的表达式不变。,在统计热力学中常选择0 K作为最低能级，因此 就是N个分子在0 K时的能量,当分子混合并且发生了化学
27、变化时，必须使用公共的能量表度。,104,从自由能函数计算平衡常数,自由能函数（free energy function）,称 为自由能函数,因为,所以,在0K时 ，所以,也是自由能函数,105,当 ，又设在标准状态下,自由能函数可以从配分函数求得。各种物质在不同温度时的自由能函数值有表可查。,求平衡常数,设任意反应,106,等式右边第一项是反应前后各物质自由能函数的差值,第二项中的 是0 K时该反应热力学能的变化值,在0 K时,107,1已知 值和各物质的自由能函数值，倒算 值。,求算 值的方法,2从Gibbs自由能的定义式求,108,两边同时加一个、减一个 ，移项整理得：,等式左方易于求得
28、，从而可得,3根据热化学中的Kirchhoff公式求,等式右方的值都可以从热化学中求得,109,4由分子解离能D来计算,设反应为：,解离产物的基态,反应物的基态,生成物的基态,用分子解离能来计算0K时的 反应热力学能的变化值,110,5由热函函数求,已知反应焓变和热函函数值，可求得 值。,对于1 mol物质，在标准状态下，有,等式左方称为热函函数。其数值可以通过配分函数求得。当T为298.15 K时， 值有表可查。,111,从配分函数求平衡常数,设反应,分子之间的关系为,112,是用分子数目表示的平衡常数，q是将零点能分出以后的总配分函数。,如果将平动配分函数中的V再分出，则配分函数用 f 表示,在平衡时，A有极小值。利用Lagrange乘因子法，求出能使A有极小值的,113,求出各配分函数 f 值，可得到平衡常数 值,对于理想气体，,114,若气体反应为,不难证明在平衡后有如下关系若气体反应为,在配分函数中，浓度C的单位是：,若单位用,平衡常数值必须作相应的换算,115,N-tm 关系,将2N个球分布在两个盒子中(不考虑顺序)的均匀分布，tm、等变化情况：,当N越大， lntm/ln就越可以互换，而tm/越小。所以说lntm=ln是对的，而说tm=是欠妥的。,
