统计热力学基本方法

1、 第五章统计热力学基本方法在第四章我们论证了最概然分布的微观状态数lntm可以代替平衡系统的总微观状态数lnO,而最概然分布的微观状态数又可以用粒子配分函数来表示。在此基础上，为了达到从粒子的微观性质计算系统的宏观热力学性质之目的，本章还需重点解决以下两个问题（：1）导出系统的热力学量与分子配分函数之间的定量关系；（2）解决分子配分函数的计算问题。5.1热力学量与配分函数的关系本节的主要目的是推导出系统的热力学函数与表征分子微观性质的分子配分函数间的定量关系。在此之前先证明0=-1/（kT）一求待定乘子0对独立可别粒子系统：lnQ将Stirling近似公式代入、展开得lnt=ln(N!n単)=

2、lnN!+工NlngmN!.11.lnQ=NlnN+YNlng-XnlnNXlnN!.代入Boltzmann关系式(46)得iiiiiS=k(NInN+XNlng-工nlnN).iii按Boltzmann分布律公式Ni=igiexp（0i），代入上式的lnN.中,利用粒子数与能量守qiii恒关系得独立可别粒子系统：独立不可别粒子系统：上式表明S是（U,N，0丿的函数,偏微分法则对（51a,b）式微分结果均为S=k(Nlnq-0U)(51a)S=k(Nlnq-0U-lnN!)(51b)而0是U，N，V的函数，当N一定时，根据复合函数的V,N而丿V,52)又q=Xgexp(Ps)所以.代入（52）

3、式得=1xq.V,N.=-k0对照热力学中的特征偏微商关系亜E丿V,N=丄qV,N空JdU丿V,N=丄便可以得到p_Tgexp(Ps)=U.53)iN1kTl二热力学函数U,S,F与粒子配分函数q的关系1热力学能U由（53）式得U=N但nq,将B一丄代入得*押丿kTV,NU=NkT2plnq（54a）V,N系统的摩尔热力学能U=RT2（沁）（54b）mI8T丿V,N由于（53）式对独立可别与独立不可别粒子系统具有相同的形式,所以（54）式适用与整个独立粒子系统。上面是从Blotzmann关系式出发导出热力学能与粒子配分函数间的关系，此关系也可以直接由Blotzmann分布律推出。独立粒子系统的

4、热力学能是所有粒子运动能量的总和,且平衡时粒子在各种形式的运动能级上的分布均服从Blotzmann分布律，所以U=艺N恪=艺Ngi列yJ灯)N*iii=0i=0qi=0&gexp(-/kT)i55)因为gi和込均与温度T无关，则(dq/dT)=d芳V,N*gexp(-/kT)iii=0V,NkT2iii=0gexp(-/kT)i?gexp(-s/kT)=kT2Gq；dT)iiiV,Ni=0代入(55)式得U=NkT2(5q-dT)=NkT(dnq：6T)V,NV,N2熵S将0=1/kT代入(5la)与(5lb)式得S（可别粒子系）=kNlnq+U(56a)S（不可别粒子系）=kln必+U（56

5、b）*N!丿T将（54a）式代入得S与q的关系S（可别粒子系）=kNlnq+NkTdlnq（57a）IdT丿V,NS（不可别粒子系）=klnfqN丿+NkTf池（57b）IN!丿IdT丿V,N再利用关系式S=klnQ=klnt可以得到最概然分布微观状态数与粒子配分函数间的关系mtm（可别粒子系）=qNexp（U/kT（58a）tm（不可别粒子系）=qNexp（U/kT）/N!（58b）显然最概然分布的微观状态数可以用粒子配分函数来表示，由此可见粒子配分函数在统计力学中占有极其重要的地位。3Helmholtz自由能F将（54）和（56）式代入定义式F=UTS，贝V F=NkT2(dInq/8T)

6、VnTkln(q/N!)+NkT(dlnq/ST)J得F(不可别粒子系)=kTln(qN/N!)(59a)或F(可别粒子系)=NkTlnq(59b)根据Stirling公式，(59a)式也可以写成：F(不可别粒子系)=NkTln(q/N)+l(59c)其它热力学性质与粒子配分函数？的关系1压力p将(59)式代入热力学关系p=_(6FI叫N,贝yp=F/叫N=-6-kTln(qN/N!)I叫Np=NkTlnqISV)TN(510)上，v2焓H将(54)和(510)式代入定义式H=U+pV,得H=NkT2lnqIdT)yN+NkTVlnqISV)T,N(511)3Gibbs自由能G将(59)和(5

7、10)式代入热力学关系G=F+pV,得G(不可别粒子系)=-kTln(qN/N!)+NkTVlnq1叫,N(512a)G(可别粒子系)=kTln(qN)+Nk7Vlnq叫,N(512b)上面以粒子配分函数表示出了独立粒子系统的五个主要热力学状态函数U，S,H,F和G。这些公式是联系物质的微观结构与宏观热力学性质的基本关系式当知道了q的具体形式后，就可以求得这些热力学函数。另外,由此出发,利用其他热力学关系式如CV=U18T)V；Cp=(8UIQT)p；卩=(8FI叫V等即可求得任何需要的热力学性质。从以上这些结果可以看出，可别粒子系统和不可别粒子系统的内能U和恰H的表达式完全相同，只是热力学函

8、数S,F,G相差一些常数项。这是由于两种系统的微观状态数不同，导致S不同，当然与S有关的F和G也就有所不同。但是，在求这些热力学函数的差值时，这些常数项即可相互消去。零点能选择所产生的影响各种能级的能量值都与零点能的选择有关。关于零点能的选择一般有两种方式：绝对零点标度，即选择共同的零点。这样，粒子的各种运动形式的基态能量就有一定的数值町。例如，振动基态能量0=加/2；相对零点标度，即选择各种运动形式自身的基态能量为能量标度的零点。这样，粒子基态的能量值就规定为零。例如，在这种零点标度下，振动基态的能量0=0。1对配分函数的影响设选择绝对零点标度时，i能级的能值为&,而选择相对零点标度i时i能

9、级的能值为显然=&“(513)iii0根据定义，当规定基态的能量为S0时(采用绝对零点标度)配分函数q为q二工gexp-sI(kT)二exp-eI(kT)工gexp-(e-e)I(kT)(514)ii0ii0ii当规定基态的能量为零时(采用相对零点标度)配分函数q为q二工gexp-s/(kT)二iii515)工giexp-(si-s。)/(kT)i比较上述两式得q=qexP-so/(kT)(516a)也可写作：Inq=Inq-&0/(kT)或Inq=InqU0J(RT)(516b)其中U0,m=Ls0是绝对零度时(各种运动形式均处于基态)系统的摩尔热力学能。0,m02对热力学函数表达式的影响为

10、了方便，统计力学常把基态的能量规定为零，即采用相对零点标度。当用q代替q时，热力学函数的表示式应进行相应的修正。U二NkT2(dlnq.cT)V,N(1)热力学能U将(516b)式代入(54)得=NkT2(C(lnq-s/kT)CTTOC o 1-5 h z0V,N HYPERLINK l bookmark74 o Current Document =NkT2(ClnqCT)+NsV,N0所以U=NkT2(Clnq；CT)+U(517) HYPERLINK l bookmark76 o Current Document V,N0U0表示全部分子都处在基态时系统的能量。可见，零点能的选择对热力学

11、能的表达式产生影响，二者相差一个U0项，但是不管零点能如何选择都不会影响XU。由于H，F，G均与U有关，所以零点能的选择对这三种热力学函数的统计表达式都会产生影响。若用q代替q，在H，F，G的表示式中应增加一个U0项。(2)熵S将(516b)式代入(56a)式得S=Nklnq+NkTClnqCT)V,N=Nklnq-s/(kT)+NkTClnq-s/(kT).CT00jV,N=Nklnq+NkT(Clnq.CT)V,N可见，零点能的选取对熵s没有任何影响。用类似的方法可以证明,p和CV也不受零点能选择的影响。5.2分子配分函数的计算系统各热力学函数均可用粒子的配分函数来表达。正是粒子配分函数将

12、系统中粒子的微观运动形态与系统的宏观性质联系起来。下面的任务就是依据量子力学所提供的有关粒子结构及运动形态的结果，计算出微观粒子的配分函数。由于前面的讨论局限于独立粒子系统，这类系统的最好例子就是理想气体，所以下面将只讨论如何计算理想气体分子的配分函数。一配分函数的析因子性质作为近似，我们可以认为分子的运动形式由彼此独立的平动、转动、振动、电子运动和核自i=1 旋运动所组成。因此任一能级的能量.应等于各种运动形式的能量之和：i&=,+e,+e,+e.+&.(518)it,ir,iv,ie,in,i同时，任一能级的简并度gi等于各种运动形式的简并度之乘积：gi=gt,gr,igv,ige,ign

13、,i(519)上两式中右端各因子依次分别表示平动、转动、振动、电子运动、核自旋运动的能量和相应的简并度。将(518)式和(519)式代入配分函数的定义式，则q二工gexp-&/(kT)iii+&+&+&+&)/(kTJt,ir,iv,ie,in,ie,in,i二工gggggexpLGt,ir,iv,ii从数学上可以证明，几个独立变数的乘积之和等于各自求和的乘积。于是上式可写作q=q-q-q-q-qtrven其中q二工gexp-&/(kT)称为平动配分函数tt,it,iiq=工gexp-&/(kT)称为转动配分函数TOC o 1-5 h zrr,ir,iiq二工gexp-&/(kT)称为振动配分

14、函数vv,iv,iiq二工gexp-&/(kT)称为电子配分函数ee,ie,iiq=工gexp-&/(kT)称为核自旋配分函数nn,in,ii520)(521a)(521b)(521c)(521d)(521e)由于平动是分子的外部运动，转动、振动、电子运动和核自旋运动是分子的内部运动，所以q=qtqi(522)其中q=qqqq(523)irven称为分子的内配分函数。可见，分子配分函数等于各种运动形式的配分函数的乘积。这一规律叫作配分函数的析因子性质。二各种运动形式配分函数的计算1平动配分函数根据(35)式，三维平动子的能量&t,izh2n2、n2n2+rb2所以q=工gexp-&/(kT)=

15、艺艺艺exptt,it,ih2n2n=1n=1n=1xyz8mkTa2n2n2、+ub2C2丿=另expn=1xh28mkTa2丿-艺exph28mkTb2-另exp丿n=1zh28mkTc2丿面分别介绍如何计算理想气体分子各种运动形式的配分函数。=qt,xqt,yqt,z 其中qt,qt和qt分别表示在x,y,z三个坐标方向运动的一维平动配分函数。需要指出，TOC o 1-5 h z匚x，y，z当能级为&时，由于n,n,n不同，应有不同的微观状态数，因此上式第一等式中的g是ixyzt,i该能级的简并度，其求和是对所有能级求和。但在第二等式和第三等式中，求和是对所有的量子态n，n，n求和，它已

16、经包括了全部可能的微观状态，因此就不再出现g.项了。xyz/、1为运算方便起见，令h=a2，贝yqexpHa2n2/ HYPERLINK l bookmark114 o Current Document 8mkTa2t,xxnx=1对通常温度和体积的理想气体a21（例如300K,a=1.0cm条件下的氢分子a2=3.96X10-17）,即平动能可作为连续处理，所以上式中的求和可用积分代替，即qt,xsexp(-a2n2)dnJexp(-a2n2)dn1xxxx根据积分公式Jex1(p-(ax2)dx二一21a丿2可得，q1t,x2(a2丿2a将a的值代入得qt,x(2兀mkT)1524)用类似的方法可以得到qt,y(2兀mkT)12-bqt,z(2兀mkT)剧=一-c所以q二q-qtt,x-qt,yt,zQm空.abc=皿仏.vh3h3(525)把式中的常数值以及m=Mr/（1000L）代入后得qt=1.88x1026（MrT）3/2V（526）式中Mr为相对分子质量，L为Avogadro常数，V的单位是m3。因此只要知道相对分子质量M及系统的温度T和体积V,便可以很方便地由（526）式计算分子的平动配分函数。r另外，对平动运动，处于基态时t00（例如298.15K时1.0dm3体积内的氮分子t0=10-26J）,所以q=qt=（加mkT）32.v（527）t,otth32.转