第六章-近独立粒子的最概然分布(习题课)

第六章 近独立粒子的最概然分布(习题课)本章题型一、基本概念：1、粒子相空间、自由度；广义坐标、广义动量；粒子微观状态、系统微观状态；经典相格与粒子微观状态；系统宏观态与系统微观态。

2、等概率原理（统计物理学的基本假设）：平衡态孤立系统的各个微观态出现的概率相等。

最概然分布作为平衡态下的分布近似。

3、近独立粒子孤立系统的粒子分布和与一个分布相对应的系统的微观状态数及各分布出现的几率、最概然分布。

,,,,21l τττ∆∆∆,,,,21l εεε}{l a,,,,21l ωωω,,,,21l a a a与分布}{l a 对应的微观状态数为()l a Ω分布{}l a 要满足的条件是：N a ll =∑E =∑ll l a ε系统总的微观状态数()()lm man a l a a lΩΩ=Ω∑~总系统某时刻的微观状态只是其中的一个。

在宏观短，微观长时间内（一瞬间）系统经历了所有的微观状态()()lm man a l a a lΩΩ∑~----各态历经假说。

且各微观态出现的概率相等()()lmman a l a a lΩ≈Ω=∑11ρ()le a a l lm l βεαωδ--=⇒=Ω0ln ---玻耳慈曼分布。

此分布（宏观态）的概率为()()()()()()1=ΩΩ≈ΩΩ=Ω=∑lmman lm man a l lm man lm man lm a a a a a a p lρ 即：最概然分布几乎就是孤立系统的平衡态分布。

4、热力学第一定律的统计解释：Q d W d dU +=l ll l ll l l da d a dU a U ∑∑∑+=⇒=εεε比较可知：l ll d a W d ε∑= l ll da Q d ∑=ε即：从统计热力学观点看，做功：通过改变粒子能级引起内能变化； 传热：通过改变粒子分布引起内能变化。

二、相关公式 1、分布与微观状态数①、 ()l a l lll l B M a a ω∏=Ω∏!N!.. ②、 ()∏--+=Ωl ll l l E B a a a )!1(!)!1(..ωω ③、 ()∏-=Ωll ll l D F a a a )!(!!..ωω④、 ()l a r l l ll l cl h a N a ) ( ! !ω∆∏∏=Ω 2、最概然分布玻耳兹曼分布le a l l βεαω--=玻色-爱因斯坦分布1-=+l e a ll βεαω费米-狄拉克分布1+=+l e a ll βεαω本章题型※、第一类是求粒子运动状态在μ空间的相轨迹:关键是由已知条件写出广义坐标q 和广义动量p 满足的函数关系()0,=p q f 。

※、第二类是求粒子能态密度()εD ；已知粒子的哈密顿量H 与广义坐标q 和广义动量p 满足的函数关系()p q H H ,=，求粒子能态密度()εD 。

不同方法有不同步骤，方法有：方法一：量子力学方法。

第一步，解薛定谔方程()()p q ,p q,H ψ=ψε，求能量本证值i ε 第二步，求出粒子能量小于ε的量子态数()εω第三步，求出粒子能量在ε到εεd +范围的量子态数()εεd D 。

方法二：半经典近似法。

该方法的依据是：对自由度为r 的一个粒子，对每一个可能的状态对于μ空间中大小为r h 的一个相体积元，因此，粒子能量小于ε的量子态数为()()⎰⎰<=εεωp q H r hdqdp,由此求得粒子能量在到范围的量子态数()()εεεωεεd d d d D =。

计算步骤：第一步、写出粒子自由度r 和粒子哈密顿()p q H H ,=。

第二步、由()()⎰⎰<=εεωp q H r hdqdp,求出粒子能量小于的状态数。

第三步、求出粒子态密度()()εεωεd d D =。

[例1]、对于二维自由粒子，在长度L 2内，求粒子在ε到εεd +的能量范围内量子态数()εεd D 。

方法一：解，量子力学方法:边长为L 的正方形平面内，粒子哈密顿算符的能量本征方程为()εϕϕϕ=+=22ˆ21H Y X P P m设：()()()y Y x X y x =,ϕ 则22222222222112 εεm dy Y d Y dx X d X XY XY y x m -=+⇒=⎪⎪⎭⎫ ⎝⎛∂∂+∂∂- 2222222222;1;1εm k k k dy Y d Y k dx X d X y x y x =+-=-=其中 解得：()()()()()y p x p iy k x k i y x y x e e y Y x X y x ++=== A1A 1,ϕ 利用周期性边界条件：⎪⎭⎫⎝⎛=⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛=⎪⎭⎫ ⎝⎛-2L ,2L ,;,2L ,2L x x y y ϕϕϕϕ得：2,1,0;,2;2±±===y x y y x x n n n Lp n L p ππ 由上式可知，量子数y x n n ，完全决定了粒子的量子状态。

以y x n n ，为直角坐标轴，构成二维量子数空间，每一组数()y x n n ，对应一个点，它代表一个量子态，这种点成为代表点，此空间中边长为1的一个正方形（面积为1）内有1个代表点，即相应于1个量子态。

由()()2222222221y x y x n n mLp p m +=+= πε可知，在数空间中能量ε的等能线为半径()2122221222R ⎪⎪⎭⎫ ⎝⎛=+= πεmL nn yx 的圆，它所包围的面积为2222R πεπmL =，而单位面积对应1个量子态，所以粒子能量小于ε的量子态数为()222 πεεωmL =，所以粒子在ε到εεd +的能量范围内的量子态数()()επεεεωεεmd h L d d d d D 222== 其中：()m hL D 222πε=为态密度，显然此情况在数空间态密度是均匀的。

方法二： 解，半经典方法:由()2221yx p p m+=ε可知，在二维动量空间中，等能线满足εm p p y x 222=+，等能线为半径等于εm 2的圆，由此求得粒子能量小于ε的量子态数：()επεωεm hL h dp dxdydp Am p p yx y x 2222222==⎰⎰≤+所以粒子在ε到εεd +的能量范围内的量子态数()()επεεεωεεmd hL d d d d D 222== ※、第三类确定孤立系统的粒子分布和与一个分布相对应的系统的微观状态数及各分布出现的几率或求最概然分布。

[例2]：（1）假设某种类型分子的许可能级为0、ω、ω2、ω3、……，而且都是非简并的，如果体系含有6个分子，问与总能量ω3相联系的是什么样的分布？并根据公式∏∏=Ωla l ll la ω!N!M.B 计算每种分布的微观态数D Ω,并由此确定各种分布的几率（设各种微观态出现的几率相等）。

（2）、在题（1）中，如0和ω两能级是非简并的，而ω2和ω3两个能级分别是6度和10度简并。

试重复上面的计算。

解：（1）粒子的在各能级的分布可以描述如下：能 级 ,,,4321εεεε， 能量值 ωωω,32,,0简并度 ,11,1,1 ，分布数 ,421,,a a a分布{}l a 要满足的条件是：6==∑N a ll ， ωε3E ==∑ll l a满足上述限制条件的分布可以有： {}{} 0,1,0,0,5a :D l 1={}{} 0,0,1,1,4a :D l 2={}{} 0,0,0,3,3a :D l 3=则各分布所对应的微观态数为：615!6!1D =⨯=Ω 3014!6!2D =⨯=Ω 2013!3!6!3D =⨯=Ω 所以此种情况下体系的总的微观状态数为56321=Ω+Ω+Ω=Ω总 各分布的几率为：107.056611D D ==ΩΩ=总P 536.0563022D D ==ΩΩ=总P 357.0562033D D ==ΩΩ=总P （2）粒子的在各能级的分布可以描述如下：能 级 ,,,4321εεεε， 能量值 ωωω,32,,0 简并度 ,106,1,1 ， 分布数 ,421,,a a a分布{}l a 要满足的条件是：6==∑N a ll ， ωε3E ==∑ll l a满足上述限制条件的分布可以有： {}{} 0,1,0,0,5a :D l 1={}{} 0,0,1,1,4a :D l 2={}{} 0,0,0,3,3a :D l 3=则各分布所对应的微观态数为：60015!6!1D =⨯=Ω 08164!6!2D =⨯=Ω 2013!3!6!3D =⨯=Ω 所以此种情况下体系的总的微观状态数为260321=Ω+Ω+Ω=Ω总 各分布的几率为：230.02606011D D ==ΩΩ=总P 692.026018022D D ==ΩΩ=总P 077.02602033D D ==ΩΩ=总P [例3]：设系统含有两种粒子，其粒子数分别为N 和N ’.粒子间的相互作用很弱，可看作是近独立的。

假设粒子可分辨，处在一个个体量子态的粒子数不受限制。

试证明，在平衡态下两种粒子的最概然分布分别为：l e a l l βεαω--=和'--'='le a l l βεαω其中l ε和'l ε是两种粒子的能级，l ω和'l ω是能级简并度。

证： 粒子A 能级，粒子数分布：l ε——{a l }——简并度l ω 粒子B 能级，粒子数分布：'l ε——{a ’l }——简并度'l ω 体系两种粒子分布要满足的条件为： N a ll =∑，N a ll '='∑ E =''+∑∑ll l ll l a a εε分布{}l a ，对应的微观状态数为∏∏=Ωla l ll la ω!N!1分布{}l a '，对应的微观状态数为∏∏''''=Ωla lll l a ω!!N 2 则系统的微观态数为21Ω⋅Ω=Ω上式表明：当第一类粒子的分布为{a l }，而同时第二类粒子的分布为{a ’l }时系统的微观态数。

在平衡下两种粒子的最可几分布是对应于在限制条件N a ll =∑，N a ll '='∑E =''+∑∑lll lll a a εε下使21ln ln Ω⋅Ω=Ω为极大的分布。

利用斯特林公式可得：l ll l ll l ll l ll a a a N a a a N ωω''+''-''++-=Ω⋅Ω=Ω∑∑∑∑ln ln ln N ln ln ln N ln ln 21由0ln 21=Ω⋅Ωδ，得 0ln ln ln 21='⎪⎪⎭⎫⎝⎛''-⎪⎪⎭⎫ ⎝⎛-=Ω⋅Ω∑∑l l l l l l lla a a a δωδωδ 而由限制条件可得：0=∑llaδ，0='∑ll a δ0=''+∑∑lllllla a δεδε引入拉氏不定乘子βαα，，'，得0ln ln ln 21='⎪⎪⎭⎫ ⎝⎛'+'+''-⎪⎪⎭⎫ ⎝⎛++-=⎪⎭⎫⎝⎛''+-''--Ω⋅Ω∑∑∑∑∑∑l l l l l l l l l l l l l l l l l l l l a a a a a a a a δεβαωδβεαωδεδεβδαδαδ根据拉格朗日未定乘子法原理，每个l a δ及l a 'δ的系数都等于零，所以得：[][]⎩⎨⎧'-'-'='--=⇒⎪⎪⎭⎪⎪⎬⎫='+'+''=++l l ll l l l l l l lla a a a εβαωβεαωεβαωβεαωexp exp 0ln 0ln 讨论：（1）、上面的推导表明，两种粒子各自遵从玻耳兹曼分布，两分布的α，α'不同，但有共同的β，原因在于开始就假设两种粒子的粒子数和能量具有确定值，这意味着在相互作用中两粒子可以交换能量，但不会相互转化。