纳米结构体系物理化学性质的理论研究方法与实例-固体表面物理化学课件.ppt

1、State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室Statistical Thermodynamics and Chemical KineticsState Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室Lecture 6State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室第六章

2、 相倚子体系（1）对实际的宏观体系，粒子间相互作用一般不可忽略。采用系综统计仍能从原则上建立其正则或巨正则配分函数。理论处理中，主要问题都集中在如何求解因粒子间相互作用而产生的构型配分函数项，但其计算又牵涉到势能函数的具体表现形式。由于影响粒子间相互作用的因素极端复杂，迄今有关相倚子体系的统计力学处理仍停留在不同程度的近似阶段。State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室6.1 晶体的频率分布晶体的频率分布 爱因斯坦模型中假定了晶体中所有原子的振动都是各自独立，且振动频率相同，这

3、显然过于理想化，不合事实。光谱实验表明：即使是简单晶体，其红外光谱亦呈现宽广即使是简单晶体，其红外光谱亦呈现宽广的吸收谱带，暗示晶体中各简正振动方式的频率变化必存的吸收谱带，暗示晶体中各简正振动方式的频率变化必存在一定的分布规律在一定的分布规律。如何确定晶体中原子振动的频率分布，就成为了以统计力学方法来揭示晶体的宏观热力学性质的一个极为关键问题。State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室6.1.1 晶体的能量函数晶体的能量函数晶体中原子间的距离很近，存在强烈的原子间相互作用。任

4、一原子均被其周围近邻原子包围，即中心原子处于邻近原子所产生的力场约束之中，其运动只能围绕着自身的平衡位置微弱振动，其周围势场决定了其振动的频率和振幅。右图给出了一维晶体中原子振动的势能曲线。aa线和cc线分别表示BA和B-C原子间的相互作用势，虚线为二者之迭加。推及三维晶体，其情形大致相仿。State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室现考虑一个由N个原子构成的晶体，每个点阵原子均有C个配位数，其运动自由度为s。由坐标的简正变换可将晶体中N个原子的振动化解为sN个简正振动。近似地，

5、每个简正振动又被模拟为单维谐振子振动。依经典力学，谐振子的能量包括了动能和势能两项：22212121xkpmxUxmxv)()()()(,sNKKKxsNKKvxUpmoVoVE12121(6.1)(6.2)m和k 分别为振子的有效质量和力常数。则整块晶体的总能量为V(o)表示晶体中所有原子皆处于平衡位置时的结合能，亦称为点阵能(晶格能)或构型能。State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室晶体点阵能的估算：以 AA表示拆散晶体中一对最近邻原子对必须耗费的能量。若不计次近邻和更远

6、的原子，因N个原子共可形成NC/2个最近邻原子对，即有sNKiKvioVE1)(,)(2/)(AACNoVKKvhV)(21,iEieNV),(6.3)(6.4)(6.5)(6.6)而谐振子的振动能还可由量子力学来求解其能级公式，对第K个简正振动有：可设正则系综中与量子态i相对应的成员能量为故可通过量子态求和来建立晶体的正则配分函数：(V为量子数)State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室isNKoViKvee1)(,)(sNKioViKvee1)()()(,若指定V(o)为零

7、，则有sNKhhKKeeNV12/1),(iEieNV),(6.7)(6.8)sNKiKvioVE1)(,)(振子的振动量子态系综成员量子态sNKhhoVsNKKvoVKKeeeqe1211/)(,)(ioVKiKvee)(,)(然而，真实的情况是3N个谐振子的特征频率 K未必全同，因此进一步求解晶体的配分函数就必须设法确定频率分布。State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室6.1.2 德拜连续介质模型与德拜热容德拜连续介质模型与德拜热容 对于晶体中原子振动的频率分布，迄今最实

8、用的处理方法为德拜连续介质模型：设想晶体是一块各向同性的连续介质，原子在晶体中振动宛若传播于这块介质的弹性波；其中每一简正振动方式的波动又相当于一个驻波。dNdnD329)(6.9)3233DNnP)()(即为分布于 +d 间隔内的振动方式数，D则为最大的特征频率，即德拜频率德拜频率,与晶体中原子间的力常数及原子约化质量有关，一般由实验确定。故可应用驻波条件导出频率分布公式为：亦有频率的几率分布函数：State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室dendnhhDD)ln()()(l

9、n12100将正则配分函数(指定V(o)为零)取其对数，即有NKhNKKKeh31311-21lnlnDDdeNdhNhDD0230331929)ln()(DDhNdNhD83329033(6.10)(6.11)(6.12)sNKhhKKee121/4331D00DDdnNdP)()(若各简正振动的频率变化是近似连续的，则上式中求和可由积分替代，U0,V为晶体中3N个简正振动的基态能量之和。(Why?)DdeNUhDV023019)ln(,其中VU,03233/DNnP/)()()(State Key Laboratory for Physical Chemistry of Solid Sur

10、faces厦门大学固体表面物理化学国家重点实验室令hkThxhkThTxkhDDDDDD-,33020291ln31ln)1ln(DDDxxxxxxxdxxdxexDDNxNUDVln3ln,01)1ln(nnxxnee45122114143110202nnnnxxxnnndxxexdxexD!)(ln(6.13)(6.14)(6.15)(6.16)(6.17)(6.18)dxexxNUxxDVD)ln(ln,190230)ln,340/(5DVxNU(1)高温状态(2)低温极限状态因 x xD 1,有State Key Laboratory for Physical Chemistry of

11、 Solid Surfaces厦门大学固体表面物理化学国家重点实验室进而可分别导出高、低温极限下晶体的热力学性质：(1)高温状态kNTUCTNkNkxNkkUkSVNVDD0333443)lnln(lnln,3433043434443512512545445159DDDDVDDDRTTkNCTNkxNkxNkS;(6.19)德拜低温热容德拜低温热容立方定律立方定律(6.20)(6.21)159430DVVNxNkTUU,ln(6.22)(6.23)(6.24)NkTUUVVN30,ln(2)低温极限State Key Laboratory for Physical Chemistry of S

12、olid Surfaces厦门大学固体表面物理化学国家重点实验室6.2 实际气体的梅逸理论实际气体的梅逸理论 实际气体的非理想性完全来自分子间的相互吸引或排斥。在众多已建立的描述实际气体的经验或半经验的状态方程中，以范德华方程和维里方程应用最广：NkTNbVVaNp)(22.)()(12VTCVTBVNkTp),.,(21),(212NpKKKKqqqVpmqpE 气体能量函数可表示为：Vp为气体中全部分子相互作用势能总和;实际上解析不开，一般都将它简化为两两分子对间吸引与排斥的势能和。(6.27)State Key Laboratory for Physical Chemistry of S

13、olid Surfaces厦门大学固体表面物理化学国家重点实验室如以Uij(r)表示相距为r的i、j分子间的势能,那么 jiijNiNijijNprUrUqqqV)()(),.,(11121于是求解Vp就变成了如何确立分子对的势能作用形式！因而，应用统计力学处理实际气体，其主要任务就是选择一个合适的势能函数去解剖状态方程中一系列依赖于分子性质的宏观函数，并由此导出气体的热力学性质。(6.28)State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室6.2.1 分子间势能函数分子间势能函数 分

14、子间相互作用可以形象地描述为“长程吸引、短程排斥”。分子间的吸引包括取向力、诱导力和色散力等不同表现形式。)(,)()(,)(0000rrrrUrrrUmnmrBrArU)(2)Lennard-Jones势：设分子非硬球，势能函数为(m3或m=6)(一般取 m=6,n=12)(6.29)(6.30)实际处理中，可采用经验或半经验势能函数，如：(1)基松(keesom)势:假定分子为硬球，r0为最接近距离，势能函数表示为：State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室平衡距离r0时有

15、(U/r)=0,即得：1200600-2)()()(rrUrrUrU(6.31)/260ArB U(r)r-U0U(r)r-U00r0r*BAKeesomLennard-Jones故有：这叫作L-J(6-12)势能函数.2600/ArU12006002 rUBrUA，State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室6.2.2 配分函数构型积分项配分函数构型积分项 对实际气体，采用相空间来求解其正则配分函数，若假定分子间作用力不影响分子的内部运动，则体系的 可表示如下:dehNNVEN

16、31!1),(KKKKKdqdpdddxdydzdqdpdpdpdpzyx;2222zyxppppNNKpKNdqdqdpdpVmphNNV.exp.!),(1123211(6.32)(6.33)暂设气体分子仅有平动运动，则NpNKKNdqdqVdpdpmkTphN.)exp(.exp.!1123211State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室其中 即平动子单位体积配分函数，而一般称为构型积分。在理想气体场合下，Vp=0,则若将分子内部运动有关的配分函数项qint(=qrqv)

17、引入，即得即非理想气体正则配分函数。其中 K的计算是个紧要问题。KNtKNNqhmkTNNV!)(2!1),(*2322/32*2hmkTqtNpKdqdqV.)exp(.1NNKVdqdq.1KNKNtNqNqqNV!)(!)(),(*int*(6.34)(6.35)(6.36)(6.37)State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室现定义则有由此，K变为一个多项式积分。上式的展开几乎是无穷系列，其求解依旧是个难题。Mayer提出了著名的集团积分理论，从原则上建立起 K的可解模

18、式。ijrUfeij1)(NjitsstijjiijNjiijNjirUNjiijKdqdqfffdqdqfdqdqedqdqrUij.)1(.)1(.)(exp(.111)(1(6.38)(6.39)State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室此式右端共含64个独立项，除首项1外，含6项，有15项等等。Mayer理论的基本构思：将 展开后的各项按分子集团重新归类组合，并以集团分布来囊括 的全部内容。6.2.3 集团积分集团积分 为便于阐明Mayer集团积分理论，先考虑一个由四分

19、子组成的“气体”体系，并以1、2、3、4分别标记分子。据此有342423141312.1)1(ffffffffffjitsstijjiijjiij jiijfjitsstijffjiijf)1(jiijf)1(6.40)State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室(1)集团排列集团排列 展开后任何一个fij都关联着相互作用的一对分子对，从而每个独立项都可看成是一定数目的分子通过fij连结起来的排列单元。因而可以利用图形来表示这样的排列单元，譬如，在i,j分子对之间连线来代替对应的

20、fij。这种图形为集团排列图。jiijf)(1State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室四分子体系集团排列图ml 表示l分子集团分布数，括号内数为同类结构型式的排列图数。m1=21234am1=4(1)b1234m2=1(6)1234m2=2c(3)1234d(12)m1=1m3=1f12f13f12f34f121234e m1=1m3=1(12)f12f13f2312341234f12f13f14f12341234m4=1(4)m4=1(12)f12f13f14f23f24f

21、3412341234gm4=1m4=1(12)(3)f12f13f14f23hif12f13f14f23f34m4=1(6)m4=1(1)1State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室(2)集团分布和集团组合集团分布和集团组合 对一个含N个分子的实际气体，其可能产生的集团排列图将是一个天文数字，只能根据分子的集团组合来说明其分布。如前图b,其单分子集团有2，双分子集团为1，构成了一种集团分布。所谓集团分布，即对一个给定的集团排列，分别指明每一类分子集团的数目各为多少，并以ml表示

22、l分子集团的分布数。如下图，其集团分布数ml为：12324252627281617181920219101112131423456781522m1=7;m2=3;m3=2m4=1;m5=1;m6=0而实现该种集团分布的组合样式数相当大！State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室 对一个指定ml（l=1,2,l)的集团排列，各分子集团的连结肯定存在许许多多的组合式样，其可能实现的组合方式数即称该集团分布ml的组合方式数。同理，对一个含有N个分子的气体体系，其可能产生的集团分布式样

23、当然也是一个大数。这里亟需解决的问题是，怎样计算每一种集团分布ml的组合方式数(A)以及全部可能分布ml的组合方式数。根据排列组合公式，对一给定的集团分布ml，其对应的组合方式数当如下：llmmlNAl!)!(!(6.41)State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室再将一切可能实现的集团分布所产生的A总加起来，即得该分子体系全部集团分布的组合方式数：此处，对ml的限制条件是：即任一集团分布所包含分子数均必恒等于体系分子总数。注意：这里给出的组合方式数并未概括各分子集团具体实现的

24、排列连结式样数。因而上式的 A实际上并不等于全部的集团排列图数。这种不足将由集团函数定义得到补充。!)!(!lllmllmmmlNANmlll(6.42)(6.43)State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室(3)集团函数与集团积分集团函数与集团积分 为便于按分子集团归类，定义集团函数如下：记S1,2l为l分子集团的集团函数;并令对单分子集团，因其不与任何分子发生关联,故S1=1;对双分子集团，S12=f12。三分子集团(l=3)，有四种不同的连结排列样式，分别对应f12f13

25、、f12f23、f13f23、和f12f13 f23等四个独立项，其中三个fijfis是等价的，于是令 S123=f12f13+f12f23+f13f23+f12f13f23=3 f12f13+f12f13f23 (6.43)四分子集团的连结排列样式更多，S1234有38个独立项：S1234=16 f12f13 f14+12 f12f13f14f23+3 f13f14f23f24 +6 f12f13f14f23f34+f12f13f23f14f24f34 (6.44)如此类推，一个l分子集团的集团函数所包含的独立项数必与该分子集团的排列连结样式数相当！State Key Laboratory

26、for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室 显然，各分子集团所包括的独立项正是 K式中待积分的内容。为此，定义集团积分：于是有其中 称为1积分。又llldqdqSVlb.!11.2,11111111dqVdqSVb12112211222121!21dqdqfVdqdqSVbdrrrf021)(422132112333121!31dqdqdqSVb3212313122!21dqdqdqfffV(6.45)(6.46)(6.47)(6.48)(6.49)(6.50)State Key Laboratory for Phys

27、ical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室四分子集团积分更复杂，可确定为(4)K的集团积分式的集团积分式 将任一集团组合样式之全部S12l连乘,称做集团组合函数集团组合函数，而与其对应的集团组合积分亦被确定为因此，一个集团组合积分代表 K积分的一部分。进一步，一种集团分布ml对 K的贡献无疑将等于 4321342423141312242314131224231413332131432112344)63(.!314132.!41dqdqdqdqfffffffffffffffVdqdqdqdqSVb ),.2,1()!(llVblIiml

28、lllmllmlllmmVbNVblmlNAIlll!)(!)!(!)!(!1(6.51)(6.52)(6.53)State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室则这样一来，实际气体的构型积分被整顿成集团积分的表达式子。进一步可导出实际气体的正则及巨正则配分函数，分别为：（因 )再令 (即单位体积下l分子集团配分函数)故!)(!llmllmlKmVbN!)*)(!*)(),(llmllmllKNmVbqNqNVlllbqq*)(!)(),(llmllmlmVqNV1)(exp),()

29、,(lllNNVqeNVVe(6.54)(6.55)lllmN(6.56)(6.57)(6.58)State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室6.2.4 粗放处理及其验证粗放处理及其验证利用以上给定的巨正则配分函数，可求出实际气体各状态函数平均值如下：(1)(2)(3)1,)(lnlllVllVlVqVqNTbVqkTNkTbqVUllllVllV2,23)*(lnllllllqkTpVqkTkTpVln(6.59)(6.60)(6.61)(6.62)State Key Labo

30、ratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室粗放处理：粗放处理：实际气体密度公式当为：对一般气体，都甚小，在极端条件下，可取并以 2=(/q1)2代入上式，给出因q1=q*,q2=(q*)1/2,所以进一步有.32)(332211qqqlqVNlll2212qq2212121qqq)1(*1q)ln(ln*ln)(ln*ln10111kTVkTNkTNRTqkTkTkTqkT(6.63)(6.64)(6.65)(6.66)(6.67)(6.68)ln(lnln),ln*(ln*lnkTkTVkTNkTkTq

31、kTNRTqkT0State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室与实际气体通常的化学位表达式比较：)1(*1q0000NqkTNRTqkTln*ln(6.69)(6.70)(6.71)fkT ln0可知气体的逸度：pVkTNf)(11p为气体压强，为气体的逸度系数;因非理想气体，pNkT/V,不能直接由上式确定！对式(6.62)做粗放处理：)(221qqkTqkTplll)()(1121-1.21-1VNkTkTp)/()(/2-1-111pf(6.72)(6.73)State K

32、ey Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室 进一步将(6.72)与维里方程对照：)()(1121-1.21-1VNkTkTp.)()(12VTCVTBVNkTp 可知第二维里系数满足：2-2-/11/)(/)(NTBVTB(6.74)(6.72)若将范德华方程改写为：)kTa(bVNVNkTp1 则有：2-1)kTa(b(6.75)(6.76)因此，首要问题变为 1积分的求解。State Key Laboratory for Physical Chemistry of Solid Surfa

33、ces厦门大学固体表面物理化学国家重点实验室 根据定义有：若分子间势能函数取Keesom势 则有：(6.77)若假定U(r)kT,近似地取 则有：(6.75)(6.76)drredrrrfrU02021144)()()()()()()(06000 rrrrUrrrU)()()()()(drredrrdrredrrerrUrrrUr0000202202114114(6.78)()(rUerU1)()()(KTUrrdrrrrkTUrdrrrUdrrrrr03030260030202134344344000State Key Laboratory for Physical Chemistry of

34、 Solid Surfaces厦门大学固体表面物理化学国家重点实验室以v0表示分子体积有：进而有 则气体的第二维里系数为:(6.81)类似地，亦可利用L-J势来求1积分，进而获得实际气体的各状态性质。则有：(6.84)(6.82)000044bUUab,)/(/)(kTUNNNTB0001442(6.83)()(KTUKTUrr000030301883434632430300/)/(rr与范德华方程中的参数比较：)(kTU)kTa(b0001442State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国

35、家重点实验室6.3 液体的自由体积液体的自由体积 液体的聚集状态介于气、固体之间。其密度高于气体，分子间相互作用力较大。统计力学理论处理液态体系的主要困难在于描述其运动情态主要困难在于描述其运动情态的精确模型不容易建立的精确模型不容易建立，它不存在像气体或固体那样的理想极限。处理液体体系两类近似方法：一是将液体看作是高度稠密的气体，并按照离域子体系的特性去建立描述体系微观态的分布函数，如径向分布函数法；二是基于液体分子排布的“长程无序、短程有序”特征，设想局域的液态分子呈现类似固体晶体的点阵排列，即似晶模型。State Key Laboratory for Physical Chemistry

36、 of Solid Surfaces厦门大学固体表面物理化学国家重点实验室6.3.1 似晶模型似晶模型-囚胞法囚胞法 假设液体中所有分子均置于一个周期性的点阵“格子”里，每一格子容纳一个分子，由一定数量的邻近分子包围。任一分子无不处于周围分子产生的势能场的约束之中。故可认为分子被囚禁在一个似晶格子的笼子内，其活动空间受限，这种格子称为囚胞。所谓的似晶点阵，充其量不过是一种平均的统计行为。State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室 对含N个粒子的纯组分液体，在不考虑分子内部运动的

37、情况下，其能量函数可设想为：AA为拆散液体中最近邻的一对分子对必须消耗的能量，而UK(qK)表示第K个分子自其平衡位置偏离qK距离时周围分子对它产生的势能。配位数C是一假想的动态平均值。上式同时确定全体分子的相互作用能为：类似固态晶体，定义)(21(21).(2121121KKKNKAANPKNKqUpmNCqqVpmENKKKNKKKAANPqUqUNCqqV1110V21)()()().(AANCV21)0(6.86)(6.87)(6.88)State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家

38、重点实验室体系正则配分函数为：其中为液体自由体积，意即分子在其似晶格子里的自由活动空间。引入分子的内配分函数，则有NftVNzyxUNtVNNKKKNVNNKKKNKNViEqedxdydzeqedqdqqUhmkTedqdqdpdpqUpmheeNVi)(.)(.)(exp.2.)(21(exp.1),(*)0(),(*)0(112/32)0(11213)0(dxdydzezyxUf),(.NfVNftVqeqqeNV)()(),(*)0(int*)0(6.89)(6.90)(6.91)State Key Laboratory for Physical Chemistry of Solid

39、Surfaces厦门大学固体表面物理化学国家重点实验室 在极端场合下，设分子仅能紧贴着囚胞中心来回晃动，且x、y、z三个轴向的势能变化均符合谐振子模式，则即有这就回到了经典晶体的状态。倘若认真计较液体中分子的非定域属性，则前述处理必然发生问题。实际情况下，还应考虑分子交换格子。由此，可将N个格子看成是N个分子公有，任何时刻，每个格子照样接纳一个分子，但不同时刻可能为不同分子所居留。2/332/),(22 kkTdxedxdydzexkzyxUfNVNNVhkTekkThmkTeNV3)0(2/32/32)0(22),(mk21(6.92)(6.93)State Key Laboratory f

40、or Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室因此，对与某一分子，其自由体积扩大N倍，即因而有 经过修正后，设分子间作用力不复存在，则分子的自由体积遍及整个体系空间，VfV,则对应的 亦还原为理想离域子体系的正则配分函数:无论如何，似晶模型下液体的统计分析最终归结为如何求解 f？其关键是液体分子间势能函数设计,例如球势阱模型球势阱模型。!)(),(*)0(NVqeNVNfVffNV!)(!)(),(*)0(NqNVqNVqeNVNNNfV(6.94)(6.95)(6.96)State Key Laboratory for P

41、hysical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室如右图，球笼的体积为给定T、V下每个分子分摊到的空间体积，即则球笼半径为又以r0和b分别表示两分子间最接近距离和分子的不可入体积，则该分子可及的活动空间（即自由体积，图中阴影部分）为其中r0可由L-J势能函数近似取值，而b相当于范德华方程的体积因子.arNV/3/143a33/13/130)2()(34braf3/230rb(6.97)(6.98)(6.99)State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学

42、固体表面物理化学国家重点实验室6.3.2 沟通熵沟通熵 体系的内能和熵平均分别为：若将正则配分函数表示为则有最后一项klneN=Nk,即称为沟通熵沟通熵。fAANVVqNNCU*ln21ln,!)*(ln)*ln(lnNVqkVqNkkUkSNffNNfVNfNVeqeNNqe)(!)()(*)0(*)0(NffekqNkqNkkUkSln)*ln()*ln(lnState Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室是由于分子从“准定域”到非定域状态而产生的熵。体系其它热力学性质可确定为：

43、NkTqNkTNCkTFfAA)*ln(21ln)()*ln(210,TkTqkTCNFfAAVNfNAANfNAANkTCVNkTVCNVFP,ln)(21ln)(21)2(ln3/13/13/2,bNfState Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室又NC AA/2可与液体汽化能LV相当，并利用Hildebrand近似：LV N2a/V,则有:NC AA/2 N2a/V于是即or此为由囚胞法导出的液体的状态方程。其适用范围为靠近临界状态下的液相。VLVaNCVNAA22,)(21

44、)2(3/13/13/222bkTVaNPkTbVaNP)2()(3/23/122RTVbNVVaNP)2()(3/23/1022State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室6.3.3 蒸汽压公式蒸汽压公式 设在给定温度下，液体A与其蒸汽达平衡。再假定蒸汽为理想气体，则可导出液体蒸汽压公式为近似地有 (摩尔汽化热)这简直成了通常的克拉贝龙方程。在正常沸点下，得标准汽化熵为glgAlA);()(fAAglANRTkTgCkTp000ln)(2/lnfAVRTRTHplnlnRTLH

45、V)ln(0fbVRTRTHSState Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室作业：作业：pp.311-312 习题1,5 State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室6.1.3 一维晶体频谱公式一维晶体频谱公式可以根据点阵结构来剖析晶体中原子的振动，其数学处理步骤为：首先建立点阵原子的运动方程，该点阵原子乃系处在周围原子产生的势能场中振动。然后求解原子的振动方程并导出

46、频率分布函数。此处以一维晶体为例来说明这一方法。按点阵结构，原子围绕其平衡位置的振动可描绘成传播于晶体中的行波。若该点阵为有限排列，则因处两端原子固定而形成驻波。可以证明，每一个驻波相当于晶体中的一种简正振动。驻波总数即等于晶体的振动自由度。其次，能够传播的驻波其频率限制在某一确定范围内。State Key Laboratory for Physical Chemistry of Solid Surfaces厦门大学固体表面物理化学国家重点实验室对一维晶体，其许可的频率变化区间为mkmkm1;10dNdnm222)(dNdnm2)(而分布在+d之间的驻波数被确定为：当振动的波长远大于点阵中原子对的距离时，其对应的频率分布亦变为此为一维晶体的德拜频谱一维晶体的德拜频谱。三维晶体的频率分布情况更为复杂-德拜的连续介质模型及频率分布公式。(6.25)(6.26)