One Forth

P06 · 第 6 章 · 第三编 相变与综合复习

输运、涨落与统计物理综合复习

以系统与环境交换约束选择微正则、正则或巨正则系综，由配分函数连接热力学势、平均量和涨落，比较量子占据统计，并把 Boltzmann 输运、涨落响应和临界有限尺寸纳入统一框架。

难度 5/5约 200 分钟最近修订 2026-07-14
报告页面错误
预备知识相变、序参量与临界现象微正则系综与熵的统计解释正则系综、配分函数与涨落巨正则系综与化学势Bose–Einstein 与 Fermi–Dirac 统计

本章目标

  1. 按照能量、粒子和体积交换边界选择微正则、正则或巨正则系综。
  2. 从 Z 或 Ξ 的导数计算热力学势、平均能量、粒子数与响应。
  3. 推导能量和粒子数涨落关系，并说明相对涨落随系统尺度的条件性衰减。
  4. 比较 Maxwell–Boltzmann、Bose–Einstein 和 Fermi–Dirac 占据分布。
  5. 从 Boltzmann 方程与弛豫时间近似连接微观碰撞、扩散和输运系数。
  6. 在临界区和有限数值样本中正确处理相关长度、自相关及热力学极限。
页面阅读位置0% · 仅保存在此浏览器
章节未开始
本册完成进度0/6 章 · 0%
本页目录

系统边界决定概率空间

统计物理先规定宏观约束，再在允许微观态上建立概率。微观态可由经典相空间点或量子本征态标记；宏观态只给能量、体积、粒子数等少量变量。三种基本边界是：

系综与环境交换固定控制量归一化对象
微正则不交换能量和粒子E,V,NE,V,N能壳内态数 Ω(E,V,N)\Omega(E,V,N)
正则交换能量，不交换粒子T,V,NT,V,NZ(T,V,N)Z(T,V,N)
巨正则交换能量和粒子T,V,μT,V,\muΞ(T,V,μ)\Xi(T,V,\mu)

温度 TT 用 K，能量用 J，体积用 m3\mathrm{m^3}，化学势 μ\mu 用 J/粒子或 J mol−1\mathrm{J\,mol^{-1}}。下文以单粒子能量使用 J，因此 β=(kBT)−1\beta=(k_BT)^{-1} 的单位为 J−1\mathrm{J^{-1}}，指数 βE\beta E 无量纲。Boltzmann 常数 kBk_B 的单位为 J K−1\mathrm{J\,K^{-1}}。

系综不是三套互斥的物理定律，而是不同环境约束下的概率描述。在短程相互作用、稳定单相和热力学极限中，它们对体密度平均量常给出相同结果；小系统、长程相互作用、相共存或临界有限尺寸处，涨落和大偏差可能保留明显系综差异。

微正则：从态数到温度

孤立系统在窄能壳内的可及微观态按等概率先验取样。熵定义为

S(E,V,N)=kBln⁡Ω(E,V,N),S(E,V,N)=k_B\ln\Omega(E,V,N),

其中 Ω\Omega 是无量纲态数；连续经典相空间需要用适当相空间单元和能壳宽度使对数自变量无量纲。热力学导数为

1T=(∂S∂E)V,N,pT=(∂S∂V)E,N,−μT=(∂S∂N)E,V.\frac1T=\left(\frac{\partial S}{\partial E}\right)_{V,N}, \qquad \frac pT=\left(\frac{\partial S}{\partial V}\right)_{E,N}, \qquad -\frac\mu T=\left(\frac{\partial S}{\partial N}\right)_{E,V}.

两个弱耦合子系统交换能量时，总态数为乘积，总熵为和。固定总能量下最大化 SA(EA)+SB(E−EA)S_A(E_A)+S_B(E-E_A)，得到 1/TA=1/TB1/T_A=1/T_B。二阶导数为负才对应稳定最大值。这个推导不把温度当作粒子的平均动能定义；平均动能关系只在特定经典模型中成立。

正则：配分函数是生成器

小系统与大热库弱耦合时，系统处于能量本征态 ii 的概率为

pi=e−βEiZ,Z=∑ie−βEi.p_i=\frac{e^{-\beta E_i}}{Z}, \qquad Z=\sum_i e^{-\beta E_i}.

配分函数无量纲。Helmholtz 自由能、平均能量、熵和热容由

F=−kBTln⁡Z,U=−∂ln⁡Z∂β,F=-k_BT\ln Z, \qquad U=-\frac{\partial\ln Z}{\partial\beta},
S=kB(ln⁡Z+βU),CV=(∂U∂T)V,N=Var⁡(E)kBT2.S=k_B(\ln Z+\beta U), \qquad C_V=\left(\frac{\partial U}{\partial T}\right)_{V,N} =\frac{\operatorname{Var}(E)}{k_BT^2}.

最后一式是涨落—响应关系：平衡能量涨落决定热容。Var⁡(E)\operatorname{Var}(E) 单位为 J2\mathrm{J^2}，除以 kBT2k_BT^2 后是 J K−1\mathrm{J\,K^{-1}}。它要求正则平衡和固定 V,NV,N，不能直接用于任意非平稳时间序列。

例 1：两能级系统的概率、能量与热容

单个两能级单元的能量为 00 和 ϵ\epsilon，且 ϵ=2kBT\epsilon=2k_BT。配分函数

Z=1+e−2.Z=1+e^{-2}.

激发态概率为

p1=e−21+e−2=0.119,p_1=\frac{e^{-2}}{1+e^{-2}}=0.119,

平均能量 U=ϵp1=0.238kBTU=\epsilon p_1=0.238k_BT。热容可由方差得到

CkB=(βϵ)2p1(1−p1)=4(0.119)(0.881)=0.420.\frac C{k_B}=(\beta\epsilon)^2p_1(1-p_1) =4(0.119)(0.881)=0.420.

低温时激发概率趋零，高温时两态概率趋于 1/21/2，两端热容都趋零，因此中间出现有限峰。该峰是有限自由度的平滑交叉，不是热力学相变。

独立子系统、可加性与经典理想气体

若 Hamilton 算符可写成互不作用部分之和 H=HA+HBH=H_A+H_B，状态空间也作张量积或经典相空间直积，则正则配分函数分解为

ZA+B=ZAZB,FA+B=FA+FB.Z_{A+B}=Z_AZ_B, \qquad F_{A+B}=F_A+F_B.

这解释了短程系统自由能的可加性。存在界面相互作用时还会有表面修正；长程作用可能使交叉项与体积项同阶，普通热力学极限需要重新定义。

经典单原子理想气体的平动单粒子配分函数含热德布罗意波长

λT=h2πmkBT,\lambda_T=\frac{h}{\sqrt{2\pi mk_BT}},

单位为 m。对 NN 个全同、稀薄粒子，

ZN=1N!(VλT3)N.Z_N=\frac1{N!}\left(\frac{V}{\lambda_T^3}\right)^N.

1/N!1/N! 去除经典标号造成的重复计数并恢复广延熵。经典低简并条件是 nλT3≪1n\lambda_T^3\ll1，其中数密度 n=N/Vn=N/V 的单位为 m−3\mathrm{m^{-3}}，乘积无量纲。温度升高、质量增大或密度降低都会减小量子交换效应；单说“宏观样品”不足以保证经典统计。

巨正则：粒子数也成为随机量

与热粒子库接触时，微观态按 Ei−μNiE_i-\mu N_i 加权：

Pi=e−β(Ei−μNi)Ξ,Ξ=∑ie−β(Ei−μNi).\mathbb P_i=\frac{e^{-\beta(E_i-\mu N_i)}}{\Xi}, \qquad \Xi=\sum_i e^{-\beta(E_i-\mu N_i)}.

巨势记作

ΦG=−kBTln⁡Ξ=−pV\Phi_G=-k_BT\ln\Xi=-pV

（均匀宏观系统）。平均粒子数与涨落满足

⟨N⟩=1β(∂ln⁡Ξ∂μ)T,V,\langle N\rangle=\frac1\beta \left(\frac{\partial\ln\Xi}{\partial\mu}\right)_{T,V},
Var⁡(N)=kBT(∂⟨N⟩∂μ)T,V.\operatorname{Var}(N)=k_BT \left(\frac{\partial\langle N\rangle}{\partial\mu}\right)_{T,V}.

在普通单相、短程系统中，方差常与 NN 同阶，因此相对标准差按 N−1/2N^{-1/2} 衰减。临界点相关体积变大，或接近相共存时分布呈多峰，这个独立块估计会失效。

例 2：由粒子数涨落反求响应

某小系统在 T=300 KT=300\,\mathrm K 时测得 ⟨N⟩=1000\langle N\rangle=1000、 Var⁡(N)=2500\operatorname{Var}(N)=2500。相对标准差为

Var⁡(N)⟨N⟩=501000=5.0%.\frac{\sqrt{\operatorname{Var}(N)}}{\langle N\rangle} =\frac{50}{1000}=5.0\%.

由涨落关系

(∂⟨N⟩∂μ)T,V=2500kB(300).\left(\frac{\partial\langle N\rangle}{\partial\mu}\right)_{T,V} =\frac{2500}{k_B(300)}.

右侧单位为 J−1\mathrm{J^{-1}}，因为这里的 μ\mu 是每粒子能量。若使用摩尔化学势，粒子数变量也必须改成物质的量并同步使用气体常数，不能混合 kBk_B 与 mol。

量子统计与经典极限

全同非相互作用粒子的单粒子能级平均占据数为

nˉBE(ϵ)=1eβ(ϵ−μ)−1,nˉFD(ϵ)=1eβ(ϵ−μ)+1.\bar n_{BE}(\epsilon)=\frac1{e^{\beta(\epsilon-\mu)}-1}, \qquad \bar n_{FD}(\epsilon)=\frac1{e^{\beta(\epsilon-\mu)}+1}.

Bose 粒子允许同一单粒子态有任意占据；Fermi 粒子的每个单粒子态受 Pauli 限制，若另有自旋简并度，应逐态或乘简并数计算。光子数不守恒，平衡黑体辐射的化学势为零；守恒原子数的 Bose 气体一般有由粒子数约束确定的 μ\mu。

当 eβ(ϵ−μ)≫1e^{\beta(\epsilon-\mu)}\gg1，两者都趋于 Maxwell–Boltzmann 形式 nˉ≈e−β(ϵ−μ)\bar n\approx e^{-\beta(\epsilon-\mu)}。经典极限由低相空间占据决定，不等于只说“温度高”。密度、质量和量子热波长共同决定简并程度。

例 3：同一能级的 Bose、Fermi 与经典占据

若某能级满足 ϵ−μ=2kBT\epsilon-\mu=2k_BT，则

nˉBE=1e2−1=0.157,nˉFD=1e2+1=0.119,\bar n_{BE}=\frac1{e^2-1}=0.157, \quad \bar n_{FD}=\frac1{e^2+1}=0.119,

而经典近似为 e−2=0.135e^{-2}=0.135。Bose 占据高于经典值，Fermi 占据低于经典值。若差值增至 6kBT6k_BT，三者都接近 e−6e^{-6}，统计差异减小。该比较针对单粒子态；总粒子数还需对态密度积分并计简并度。

从平衡分布到 Boltzmann 输运

平衡系综不给出空间流和恢复时间。稀薄经典粒子的单粒子分布 f(r,p,t)f(\mathbf r,\mathbf p,t) 可满足

∂f∂t+v⋅∇rf+F⋅∇pf=C[f],\frac{\partial f}{\partial t} +\mathbf v\cdot\nabla_{\mathbf r}f +\mathbf F\cdot\nabla_{\mathbf p}f =C[f],

其中 C[f]C[f] 是碰撞项。f d3r d3pf\,\mathrm d^3r\,\mathrm d^3p 的归一化按具体约定给出粒子数或概率。弛豫时间近似写成

C[f]≈−f−flocτ,C[f]\approx-\frac{f-f_{\mathrm{loc}}}{\tau},

τ\tau 用 s。它假定偏离局域平衡不大，并把复杂碰撞压缩为单一时间尺度；强非平衡、守恒投影或多尺度散射需要更完整碰撞算符。

随机游走的扩散系数数量级为

D∼ℓ22dτ,D\sim\frac{\ell^2}{2d\tau},

ℓ\ell 为平均步长（m），DD 单位 m2 s−1\mathrm{m^2\,s^{-1}}。线性响应区中通量与梯度成正比，例如 j=−D∇n\mathbf j=-D\nabla n。负号表示粒子从高密度流向低密度；它来自熵增方向和构成关系，不是矢量单位的结果。

例 4：由微观步长估算扩散尺度

三维随机运动的有效步长 ℓ=2.0 nm\ell=2.0\,\mathrm{nm}，步间隔 τ=0.50 ps\tau=0.50\,\mathrm{ps}。估算

D=(2.0×10−9)26(0.50×10−12)=1.33×10−6 m2 s−1.D=\frac{(2.0\times10^{-9})^2}{6(0.50\times10^{-12})} =1.33\times10^{-6}\,\mathrm{m^2\,s^{-1}}.

在 t=1.0 mst=1.0\,\mathrm{ms} 内，三维均方位移满足

⟨r2⟩=6Dt=8.94×10−5 m=89.4 μm.\sqrt{\langle r^2\rangle}=\sqrt{6Dt} =8.94\times10^{-5}\,\mathrm m=89.4\,\mu\mathrm m.

这一结果要求步之间近似独立且长时间已进入扩散区。短于碰撞时间时运动可能近似弹道，⟨r2⟩∝t2\langle r^2\rangle\propto t^2，不能继续使用扩散公式。

局域平衡、Knudsen 数与输运边界

Boltzmann 输运把宏观场看成局域平衡参数的缓慢变化。令系统宏观变化尺度为 LgL_g，平均自由程为 ℓ\ell，Knudsen 数

Kn=ℓLg\mathrm{Kn}=\frac\ell{L_g}

无量纲。Kn≪1\mathrm{Kn}\ll1 时一个宏观梯度尺度内发生许多碰撞，密度、温度和流速可作为平滑局域变量，扩散和流体构成关系较可靠；Kn≳1\mathrm{Kn}\gtrsim1 时进入稀薄或弹道区，边界层和非局域效应显著。

弛豫时间还给出平均自由程估计 ℓ∼vthτ\ell\sim v_{\mathrm{th}}\tau。若热速度为 500 m s−1500\,\mathrm{m\,s^{-1}}、τ=2.0×10−10 s\tau=2.0\times10^{-10}\,\mathrm s，则 ℓ∼1.0×10−7 m\ell\sim1.0\times10^{-7}\,\mathrm m。在 Lg=1.0 mmL_g=1.0\,\mathrm{mm} 的装置中 Kn∼10−4\mathrm{Kn}\sim10^{-4}，连续扩散近似合理；在 Lg=0.20 μmL_g=0.20\,\mu\mathrm m 的通道中 Kn∼0.5\mathrm{Kn}\sim0.5，同一扩散系数和无滑移边界已不可直接沿用。

碰撞项必须保持相应守恒矩：对只含弹性二体碰撞，∫C[f] d3p=0\int C[f]\,\mathrm d^3p=0、动量和能量加权积分也应为零。Boltzmann 的 HH 定理在分子混沌等假设下给出趋向平衡的熵方向，但它不是对任意强关联动力学的普遍证明。局域平衡、Markov 碰撞和稀薄性都应列入适用范围。

涨落—耗散与时间尺度

平衡涨落描述系统自发偏离平均值的大小；线性响应描述弱外场下平均量的改变。两者由同一微观动力学约束。经典扩散中，随机碰撞造成均方位移，阻力又让速度扰动衰减；Einstein 关系把扩散系数和迁移率连接。具体形式依赖变量、守恒律和时间反演性质，不能只凭“都有涨落”写出常数。

时间关联函数

CA(t)=⟨δA(t)δA(0)⟩C_A(t)=\langle\delta A(t)\delta A(0)\rangle

给出记忆衰减。输运系数常可写成平衡流关联函数的时间积分，这类 Green–Kubo 关系要求系统平衡、积分收敛和变量归一化明确。临界附近相关长度和弛豫时间同时增长，出现临界慢化；静态指数不足以决定动力学指数。

临界区、有限尺寸与系综边界

连续相变附近 ξ=ξ0∣t∣−ν\xi=\xi_0|t|^{-\nu} 增大，响应按幂律增强。有限边长 LL 截断发散，伪临界峰满足尺寸标度。计算时应报告 LL 的物理单位或格距单位、边界条件、热化区间、随机种子和自相关时间。

系综等价依赖热力学极限与稳定性。相共存的能量或粒子数分布可能双峰；长程相互作用可导致非加性；小系统的表面项与体积项同阶。此时微正则热容、正则能量分布和巨正则粒子涨落可能显示不同现象，不能用“大系统都一样”跳过边界检查。

数值抽样的最小证据

Metropolis 型抽样生成的是相关 Markov 链。样本均值误差取决于有效样本数，而非保存构型数。至少应记录接受率、热化丢弃长度、积分自相关时间、独立链或分块误差，并检查能量、粒子数或序参量直方图。重加权和数据塌缩能复用附近参数数据，但不能外推到没有分布重叠的区域。

解析配分函数、直接枚举、Monte Carlo 和分子动力学回答的问题不同。Monte Carlo 步通常没有真实秒单位；分子动力学时间有物理单位，但受积分步长、势函数和恒温方案影响。不得把两者横轴都写成“时间”后直接比较输运系数。

常见误区

配分函数只是归一化常数

它确实归一化概率，同时其对温度、体积和化学势的导数生成自由能、平均量与涨落。

相对涨落总按 N 的负二分之一次方消失

该估计要求可分成大量弱相关块。临界关联、相共存、长程作用或小系统可使标度改变。

弛豫时间近似保持所有守恒律

单一 −(f−floc)/τ-(f-f_{loc})/\tau 需要恰当选择局域平衡和守恒投影；任意目标分布可能错误耗散粒子数、动量或能量。

练习

练习 1：按边界选择系综
所属知识
系综
难度
2/5

为三类实验装置选择系综，并说明为什么“温度已知”通常不属于严格微正则约束。

查看提示
逐项判断能量和粒子能否穿过边界，再列固定控制量。
查看解答
孤立刚性封闭系统选微正则 E,V,N；与热库交换能量但不交换粒子选正则 T,V,N；与热粒子库同时交换能量和粒子选巨正则 T,V,μT,V,\mu。若边界还允许体积变化，应加入压强约束并选择相应扩展系综。
练习 2：从 Z 求能量方差
所属知识
涨落响应
难度
3/5

推导 CV=Var⁡(E)/(kBT2)C_V=\operatorname{Var}(E)/(k_BT^2) 并核对单位。

查看提示
先用 U=−∂βln⁡ZU=-\partial \beta \ln Z，再对 β\beta 求导；利用 dβ/dT=−1/(kBT2)d\beta/dT=-1/(k_BT^{2})。
查看解答
∂2βln⁡Z=Var⁡(E)\partial^{2}\beta \ln Z=\operatorname{Var}(E)。又 CV=∂U/∂T=(∂U/∂β)(dβ/dT)=Var⁡(E)/(kBT2)C_V=\partial U/\partial T=(\partial U/\partial \beta)(d\beta/dT)=\operatorname{Var}(E)/(k_BT^{2})。方差非负，所以稳定正则系统的 CVC_V 非负；等式要求 Z 收敛且导数可交换。
练习 3：量子占据的经典极限
所属知识
量子统计
难度
3/5

证明 Bose 与 Fermi 分布在低占据极限趋于 Maxwell–Boltzmann 形式。

查看提示
令 x=eβ(ϵ−μ)x=e^{\beta(\epsilon-\mu)} 很大，比较 1/(x−1)1/(x-1)、1/(x+1) 与 1/x。
查看解答
x≫1x\gg 1 时，两分母都由 x 主导，故 nBE≈nFD≈x−1=e−β(ϵ−μ)n_{BE}\approx n_{FD}\approx x^{-1}=e^{-\beta(\epsilon-\mu)}。下一阶 Bose 修正为正，Fermi 修正为负。经典极限还要求各态低占据，不能只由单个高能级推断整个气体。
练习 4：扩散方程的量纲
所属知识
输运
难度
3/5

由连续性方程和 Fick 定律推导扩散方程，并核对 DD 的 SI 单位。

查看提示
在 ∂tn=D∇2n\partial_t n=D\nabla^{2}n 中比较左右单位，并从连续性方程识别通量。
查看解答
∂tn\partial_t n 单位为 n⋅s−1n\cdot s^{-1}，∇2n\nabla^{2}n 为 n⋅m−2n\cdot m^{-2}，所以 D 单位 m2⋅s−1m^{2}\cdot s^{-1}。Fick 定律 j=−D∇nj=-D\nabla n 给通量单位为粒子数⋅m−2⋅s−1\cdot m^{-2}\cdot s^{-1}。负号表示沿密度降低方向。
练习 5：有效样本数
所属知识
数值抽样
难度
3/5

说明为什么一百万个 Monte Carlo 构型不等于一百万个独立样本，并写出有效样本数估计。

查看提示
相关链的均值方差约放大 2τint2\tau_{\mathrm{int}}；先统一 τint\tau_{\mathrm{int}} 与采样间隔的单位。
查看解答
若保存 M 个等间隔样本，常用近似 Meff≈M/(2τint)M_{eff}\approx M/(2\tau_{\mathrm{int}})，其中 τint\tau_{\mathrm{int}} 以采样间隔为单位。标准误按独立样本估计再乘 2τint\sqrt{2\tau_{\mathrm{int}}}。应配合分块或多链检查；临界慢化会使 τint\tau_{\mathrm{int}} 随 L 增长。
练习 6：系综等价失效检查
所属知识
热力学极限
难度
4/5

列出至少四项判断微正则、正则和巨正则能否视为等价的条件。

查看提示
检查相互作用是否短程可加、系统是否单相、关联长度是否远小于 L，以及分布是否单峰。
查看解答
普通等价需要稳定、可加、短程和适当热力学极限。相共存双峰、ξ≈L\xi \approx L、长程非加性或表面占比显著时应分别比较各系综概率分布和大偏差函数。平均值接近也不保证涨落相同。

关系与资源

课程 · 2013

Statistical Mechanics I: Statistical Mechanics of Particles

Mehran Kardar

用于核对 P06 的系综推导、配分函数、涨落关系、量子占据分布和多粒子例题。

打开官方来源
课程 · 2014

Statistical Mechanics II: Statistical Physics of Fields

Mehran Kardar

用于核对 P06 临界现象章节的对称性、临界指数、尺度关系和重整化群适用边界。

打开官方来源

MIT OpenCourseWare 8.333 用于系综、量子统计与涨落主线，8.334 用于场、尺度与临界现象。两项均为已登记课程；正文中的概率归一化、单位、极限和算例可独立复核。