第三章 量子统计物理学基础 热力学和统计物理： 经典统计物理和量子统计物理：粒子遵从经典（量子）力学规律。 热力学：从若干（宏观）经验定律出发，通过数学上的推导获得系统的宏观性质； 统计物理：从单个微观粒子的力学运动规律出发，加上统计的假设，来描述宏观物理量的行为。宏观量是相应微观物理量的统计平均值。 经典统计物理和量子统计物理：粒子遵从经典（量子）力学规律。 平衡态统计物理和非平衡态统计物理：研究系统与时间无关的性质或系统的时间演化行为（如前两章我们已讲述的）从现在开始我们的讨论仅限于平衡态统计物理。

3.1 经典统计系综 先从经典统计出发： 给定系统的动力学状态可用系统的广义坐标q和与之共轭的广义动量p来确定。 对由N个粒子组成的系统，可记为 其中 为第i个粒子的坐标和动量。(q,p)是6N维的相空间(Γ空间)的一个代表点，称为相点，它代表系统的一个微观状态。代表点在Γ空间的运动反映系统微观状态的演化。 系统的动力学函数或力学量：表征系统的状态，并能加以观测的量，它是q,p的函数，可记为b(q,p)。其中，表征系统能量的动力学函数H(q,p)非常重要，称为哈密顿量(Hamiltonian)。 系统的运动方程（哈密顿正则方程）： 任意力学量b(q,p)的运动方程: 上面后两式称为力学量b和H的泊松符号(Poisson bracket)。

统计系综： 统计物理认为系统的动力学状态遵从统计规律性（对比牛顿力学的确定性）。即在一定的宏观条件下，某一时刻系统以一定的概率处于某一状态或某种状态范围内。并假设，宏观量是相应微观量对系统可能处的所有动力学状态的统计平均值。 如何获得统计平均值？ 大量的重复测量！ 统计系综：由大量处于相同宏观条件下，性质完全相同而各处于某一微观状态、并各自独立的系统的集合。系综在相空间里的几何表示是无数多个相点的集合。 密度函数D(q,p,t)：相点(q,p)附近单位相体积元内相点的数目。 特别地，概率密度函数ρ(q,p,t)满足归一化条件（D=Nρ,N为总粒子数）：

刘维尔定理 系综的概率密度函数在运动中不变，即 在体积元dΩ=dqdp里，经过时间dt后，代表点的增加为 而通过平面qi(对应的面积为dA=dq1…dqi-1dqi+1dqfdp1…dpf)进入的代表点为 通过平面qi+dqi走出的代表点为： 因此净进入的代表点数为： 考虑所有qi,pi我们发现 利用正则方程及其推论： 我们发现（刘维尔定理）：

3.2 量子统计系综 由N个粒子组成的系统的状态用波函数来描写：ψ(q1,q2,…,qN,t)，时刻t在(q1,q2,…,qN)找到该N个粒子的概率为 纯粹系综和混合系综： 纯粹系综：每次测量，系综中N个粒子都处于同一态 ，可以用单一态矢量来描写： 这里 是纯态态矢量。 混合系综：每次测量，系统以一定的概率可处于多个态上。混合系综是由若干纯态混合来描写，即 参加混合的态： 各态混合的概率： P1, P2, ,…,Pi, … 且 几个例子： 1. 考虑位置空间x，找到粒子处于x的概率密度为： 纯粹系综： 各 间有干涉。 混合系综： 各 间没有干涉。

统计算符 2. 算符的平均值： 考虑算符 ，其平均值为： 纯粹系综： 各 间有干涉。 混合系综： 各 间没有干涉。 考虑算符 ，其平均值为： 纯粹系综： 各 间有干涉。 混合系综： 各 间没有干涉。 统计算符 统计算符对应于经典统计物理里的概率密度函数，其在任意表象中的矩阵形式称为密度矩阵。 对混合系综，我们定义统计算符为： 若 为完全正 交归一的基矢（ ），我们有：

特点： 若 是正交归一的态矢量，则 是统计算符的本征矢，这时密度矩阵为ρij=Piδij. 统计算符的求和中若只有一项i不为零，我们回到了纯粹系综。因此我们上面的定义对两种系综都成立。 统计算符的迹为1，与表象无关。即： 统计算符平方的迹对混合系综小于1，对纯粹系综等于1。 统计算符是厄密算符，故其本征值为实数。 求统计算符的两个例子：见杨展如书第8-9页。 系综的熵算符：熵算符的定义为 而系综的熵由此为： 上式最后一个等式我们已取 为 的正交归一的本征态矢量。这称作von Neumann熵。

量子统计里的刘维尔定理 我们可以采用两种绘景： 薛定鄂绘景：态矢量显含时间，而算符不显含时间； 海森堡绘景：态矢量不显含时间，而算符显含时间。 用哪个？ 注意到 我们采用薛定鄂绘景。此时： 由薛定鄂方程 为系统的哈密顿算符，可得 所以 这就是量子刘维尔方程，其中 为量子泊松符号。

刘维尔方程的形式解： 我们可以定义演变算符： ，则 代入到刘维尔方程中我们发现 若H不显含t，则 密度算符的形式解为： 如用能量表象的完备正交矢展开，我们发现： 统计平衡时（定态），统计算符不随时间变化，这时统计算符和系统的哈密顿算符对易。若无简并，则统计算符是哈密顿算符的任意函数；若有简并，密度算符是哈密顿算符和所有与哈密顿算符对易的算符的函数。反之，若统计算符是哈密顿算符的任意函数，则其不随时间变化！

3.3 几种平衡态量子统计系综 3.3.1 微正则系综 考虑一个由封闭的、能量孤立的系统组成的系综。系统体积为V，总粒子数为N。而能量有微小变化（在E到E+ΔE之间）。 等概率假设:孤立系达到平衡态时，系统处于任一可能状态的概率相等。 平衡时统计算符和哈密顿算符对易，因此在能量表象里ρnm=Pnδnm。若总状态数为Ω(E)，由等概率假设有 任何一个物理量的平均值为： 特别地，系统的熵为：

3.3.2 正则系综 微正则系综的极值性质：对由孤立系组成的系综中，系统状态在ΔE内的一切可能分布里，微正则分布对应的熵最大（熵增加原理）！ 证：由于对所有x>0有：ln(x)>=1-1/x，设 是任一个可能的统计算符， 是微正则分布对应的统计算符，令 ，我们发现： 上式两边取迹后易知： 3.3.2 正则系综 考虑一个封闭系统，它可以与外界交换能量，但不能交换粒子。可设想为与外界大热源接触而达到统计平衡的系统。平衡时有确定的粒子数N，确定的温度T和确定的体积V。 设系统和热源组成的（孤立）复合系统的总能量为E0，系统处于能量Es（E0>>Es）。这时热源可处于能量为Er=E0-Es的任何一个状态，由等概率假设得： 但 因此

归一化后我们发现： 这里Z是配分函数 其中s 对所有粒子数为N和体积为V的微观状态求和。 考虑能量的本征矢 ，统计算符可表示为： 配分函数可写为： 任一动力学量的平均值为： 自由能定义为： 正则系综的极值性质：在具有相同平均能量的所有可能的分布里，正则分布的熵最大（熵增加原理）。 证：设 是任一个可能的统计算符， 是正则分布对应的统计算符。故有： 利用此式即容易得证：

3.3.3 巨正则系综 考虑一个开放系统，它可以与外界交换能量和交换粒子。可设想为与外界大热源和大粒子源接触而达到统计平衡的系统。平衡时有确定的化学势μ，确定的温度T和确定的体积V。 设系统和热源及粒子源组成的复合系统的总粒子数为N，总能量为E，系统的粒子数为Nn（N>>Nn），处于能量En（E>>En）。这时热源可处于粒子数为Nr=N-Nn，能量为Er=E-En的任何一个状态，由等概率假设得： 但 因此 归一化后有： 这里Ξ是巨配分函数，它常写为： 这里 是粒子数为N的正则配分函数， 是易逸度。

考虑能量为E，粒子数为N的完备本征矢，类似前面的情形容易发现巨正则系综的统计算符为： 巨配分函数可写为： 物理量的平均值为： 热力学势定义为： 巨正则系综的极值性质：在具有相同平均能量和平均粒子数的所有可能的分布里，巨正则分布的熵最大（熵增加原理）。 证：设 是任一个可能的统计算符， 是巨正则分布对应的统计算符。故有： 利用此式即容易得证：

3.4 计算密度矩阵举例 置于边长为L的立方体里的自由单粒子：杨展如书第20－21页。 波函数满足周期性边界条件，其能量本征值为 正则系综的统计算符为： 这里 是哈密顿量的本征矢。 在坐标表象下，统计算符的矩阵元为： 其中我们利用了： 而哈密顿量的系综平均值为：

2. 磁场中的单粒子：杨展如书第21-22页。 其哈密顿量为： 这里 是单粒子的自旋， 为泡利自旋算符，我们并取磁场为沿z方向。 在 为对角化的表象里，其泡利矩阵为 ，于是正则系综的密度矩阵为： 而 的系综平均值为：

3.5 三种独立粒子系统的最概然统计分布 自然界最基本的粒子分为两种： 波色粒子：粒子不可分辨，有整数自旋，系统每个状态的粒子数不受限制，多粒子系的波函数用对称波函数描述； 费米粒子：粒子不可分辨，有半整数自旋，系统每个状态的粒子数最多为一个（泡利不相容原理），多粒子系的波函数用反对称波函数描述。 经典近似：玻尔兹曼粒子，粒子可分辨，系统每个状态的粒子数不受限制。 考虑一个由大量全同独立粒子组成的孤立系统（总粒子数N，体积V，总能量E）。设εi为单个粒子的能级，gi为该能级的简并度，α表示该能级的简并态，则 与填布数(n1,n2,…,ni,…)对应的微观状态数为： 波色分布： 费米分布： 由等概率原理可知，每个可能的微观状态出现的概率相等，因此使微观状态数为极大的分布出现的概率最大（最概然分布），这时有δlnΩ=0。因总能量和总粒子数恒定，我们有限制条件：δN=0及δE=0。用拉格朗日法求极值即可得出分布。 波尔兹曼分布： （ 极限）

以波色粒子为例，利用斯特林公式 ln m!=m(ln m-1)，近似有： 因此我们有 及 因所有δni独立，上式中的所有系数必为零。由此即得波色分布： 类似可得： 对费米粒子系统： 对玻尔兹曼粒子系统：

3.6 配分函数和统计热力学 所有表征系统平衡热力学性质的热力学函数都可以用配分函数表示出来（对经典和量子统计都适用）。对量子统计来说，配分函数可以通过对统计算符求迹获得。 这里我们将简述不同系综中统计热力学的基本公式。 正则系综： 基本热 力学量为自由能F，由定义有: 两边对β求微商，即得系统的内能U： 利用 ，可以发现系统的熵为： 以β,V,N为独立变量，可导出自由能的增量： 这里 压强 化学势

定容热容量为： 能量涨落为： 能量相对涨落为： 对理想气体， 因此 ，由此可见对宏观的大系统来说能量的涨落是极低的。 巨正则系综： 基本热 力学量为热力学势Ω，由定义有: 两边分别对β和μ微分，易得平均能量和平均粒子数的表达式： 利用 ，可以发现系统的熵为： 以β,V,μ为独立变量，可导出热力学势的增量(这里记 ， )：

粒子数涨落： 由 因此粒子数涨落为： 其相对涨落为： 其中ρ=N/V为系统的密度，应用热力学关系：VdP=SdT+Ndμ 可得 因此 这里κ为系统的等温压缩系数。我们最后有： 对理想气体，可以验证粒子数相对涨落约为1/N。但在气液相变的临界点附近，压缩系数κ趋于无穷，这时粒子数的涨落很重要。如临界乳光现象，就是由大密度涨落导致的光散射效应的增强。 对费米和波色气体，我们有： 其中“＋”为费米气体。根据上面的公式可以发现相对粒子数涨落为： 其中“－”号为费米气体。 由此可知费米气体的涨落很小，而波色气体当 取任何值时涨落都不为零，相对涨落的数量级为

3.7 巨正则系综：理想气体的统计分布和物态方程 这里我们将通过巨正则系综对理想气体的统计分布和物态方程作较严格的推导（对比3.5节通过微正则系综获得的最概然分布）。 理想气体的能量和粒子数可写为： 这里 是动量为p的单个粒子的能量， 是动量为p的粒子数。因此巨配分函数可写为： 对 的求和里可能有简并。由3.5节我们知道对波色气体，粒子数为 ，简并为 的微观状态数为 ；对费米气体，粒子数为 ，简并为 的微观状态为 。 因此，考虑到简并后，我们有 动量为p的态的平均粒子数为： 其中“-”号是对波色统计，这与前面得到的相同。

物态方程： 对压强，有： 对系统总粒子数，有： 令V趋于无穷大，上面式子里的求和可以化为积分: 可以把物态方程显式地表达出来。对波色气体，注意到p=0项当 时发散,因此我们必须对p=0项单独考虑。 结果为（无简并情形）： 理想费米气体： 这里 是比容， 理想波色气体： 对波色气体，p=0的态的平均粒子占据数为 ，若 不是一个小量，它会对比容公式右边的第二项产生明显影响，它意味着系统中有限部分的粒子占据了p=0态，这与波色-爱因斯坦凝聚有关。

3.8 热力学函数的奇异性 李－杨定理 我们知道热力学函数可以用来描述统计平衡时宏观系统的热力学行为，那么自然界中经常出现的相变现象，能否用同一热力学量（的单一数学表达式）来描述？ 相变时某些热力学量有奇异性（发散），如果上面的描述可行的话，这个奇异性必然与配分函数的行为有关（因为这些热力学量可以用配分函数表达出来）！ 李－杨在理论上严格地解决了这个问题。虽然在有限体积时，不会出现热力学量的奇异性，但当体积V趋于无穷大，同时保持粒子密度N/V恒定时（这个条件称为热力学极限），奇异性可能出现（李－杨定理）。 原因：解析函数序列的极限并不一定解析（留意巨配分函数定义中的求和）。

有限体积的一个例子： 考虑N个粒子的系统，其体积为V，粒子间相互作用为短程吸引势包围的钢球势（如左下图），因此有限体积内最多只能容纳有限个粒子M(V)。如N>M(V)，必有两粒子接触，从而能量为无穷大，这时配分函数 巨配分函数可写为易逸度z的M阶多项式，因所有系数大于0，巨配分函数没有正实根，且大于1： 因此压强 是z的单调解析函数； 再考虑单位体积的倒数1/v=N/V，由于N和V均有限，N不能超过M(V)，故1/v大于0且有限。由 知，奇点 都不在正实轴上，因此1/v在物理区域也是z的解析函数。

因此P和1/v的物态方程也是解析的。 但在热力学极限： 上面右边等式的求极限和微分次序一般不能任意交换，这时可能出现热力学函数的奇异性。 李、杨证明了下面两个定理： 对所有z>0，极限 存在，且是z的连续非递减函数。如V的表面积的增加不快于 ，则该极限也与V的形状无关。 注释：当V趋于无穷大时，巨配分函数的某些根 可能会任意地趋近正实轴。此时设 是正实轴上的点，在该点的任意小的邻域内都存在巨配分函数的零点， 在 出现奇异性（其某阶导数可能不连续）。 设在复z平面上包含一段正实轴的区域R内，不含巨配分函数的零点，则在区域R内，当V趋于无穷大时， 均匀地收敛到它的极限，该极限对区域R内所有z是解析的。 注释：区域R可看作一个热力学相的定义域。分两种情况（见下页）：

(b) 多相： 以 为界，有两个区域（相），它们各自使定理2成立。 (a) 单一相（无相变）： R包含了整个正实轴。 (c) 一级相变（ 不连续） ： (d) 二级相变（ 不连续）：

一个例子 设巨配分函数为： ，这里体积V为整数。故巨配分 函数在正实轴上无零点。其根都在复z平面的单位圆上，当V趋于无穷大时，一些零点无限趋近于正实轴。 该巨配分函数在z<1和z>1有不同的极限： z<1时: z>1时: Ξ的零点分布 由上两式解出物态方程可发现该系统在z=1处存在一级相变。

前三章练习 1。考虑一个n个放射性粒子组成的系统。设在单位时间内粒子衰变的概率为γn，试写下衰变过程的主方程，并用特征函数法求解。 2。试从巨配分函数推导理想波尔兹曼气体的物态方程。