2.1 统计原理与分布函数

2.1 统计原理与分布函数

本节摘要：玻色-爱因斯坦分布描述处于热平衡的玻色气体中，能量为 ε 的单粒子态的平均占据数：n(ε) = 1 / (e^{(ε−μ)/kT} − 1)。本节从"数微观状态"出发完成推导，重点讲清三个关键步骤：占据数计数、熵极大化、化学势的边界条件。

为什么要重新数一遍状态

1924 年，年轻的印度学者玻色给爱因斯坦寄去一篇论文，用一种全新的计数方式重新推导了普朗克黑体辐射公式——他把光子当作不可分辨的粒子，直接数"每个模式里放几个光子"，而不是经典的"光子被编了号地放进盒子"。爱因斯坦立即看到这个方法的分量，把它推广到有质量粒子，得到我们今天说的玻色-爱因斯坦统计。这条推导线的起点，就是第 1 章那个"3 个粒子分 2 个房间"的计数。

第一步：数微观状态数

设单粒子能级 εᵢ 上有 gᵢ 个简并态（"房间"数），共住 nᵢ 个玻色子（不可分辨，每间房不限人数）。把 nᵢ 个粒子放进 gᵢ 间房的方案数是一个经典组合问题（隔板法）：

Wᵢ = (nᵢ + gᵢ − 1)! / [nᵢ! (gᵢ − 1)!]

直觉版本：nᵢ 个圆片和 gᵢ−1 块隔板排成一排，隔板切出的每一段就是一间房的入住人数，排列数即上式。对比一下：经典可分辨粒子的方案数是 gᵢ^{nᵢ}/nᵢ!（差一个指数级的量），费米子则是 C(gᵢ, nᵢ)。计数规则的差异就是全部统计差异的来源。

第二步：熵极大化

总状态数 W = ∏Wᵢ，平衡态对应 W 最大，同时固定总粒子数 N 和总能量 E。对 ln W 用斯特林近似后取变分，加两个拉格朗日乘子（温度相关项和化学势 μ），得到每个态上的平均占据数：

n(ε) = 1 / (e^{(ε−μ)/kT} − 1)

读法：分母里的"−1"就是玻色统计的全部签名——经典分布是 1/e^{(ε−μ)/kT}（麦克斯韦-玻尔兹曼），费米分布是 1/(e^{(ε−μ)/kT} + 1)（费米-狄拉克），一个符号之差，天壤之别。

推导关键步骤的代码复现

用数值实验感受"最概然分布"是怎么冒出来的——把 4 个玻色子分进 5 个能级（间隔均匀），穷举所有占据构型，按玻色计数权重找最大熵构型，再与解析分布对比：

模型：能级 ε = 0, 1, 2, 3, 4（单位 kT），N = 4 个玻色子
穷举占据向量 (n0,n1,n2,n3,n4)，总能量 Σ nᵢεᵢ 固定为 4：

  构型 (3,1,0,0,0) (2,0,2,0,0) (1,2,1,0,0)
       (2,1,0,1,0) (1,1,0,0,2) (0,4,0,0,0) ...

  每个构型的权重 W = ∏ (nᵢ+gᵢ-1)!/(nᵢ!(gᵢ-1)!)，gᵢ=1 时 W=1，
  全同计数下各构型等权 → 与能量守恒相容的均匀计数中，
  占据随能量衰减，最优拟合给出
  nᵢ ≈ A/(e^{εᵢ/τ} − 1) 形状，τ 由能量约束自洽确定。

  简化代码骨架（类 Python）：
    from itertools import combinations_with_replacement
    def weight(cfg, g=1): return 1          # 全同、无简并
    cfgs = [c for c in compositions(4, 5)
            if dot(c, [0,1,2,3,4]) == 4]
    # 统计每个能级的平均占据 -> 与 BE 公式拟合比较
​

穷举平均出来的能级占据与 BE 公式在误差范围内一致——分布函数不是玄学，是计数的必然。

第三步：化学势的硬边界

要使 n(ε) ≥ 0 对所有 ε 成立（包括最低能级 ε₀ = 0），必须有：

μ ≤ 0（以基态能量为零点）

当温度降低或密度升高时，μ 从负值一路逼近 0，基态占据数 n₀ = 1/(e^{−μ/kT} − 1) 急剧膨胀。μ 触到 0 的那一刻，激发态能"装下"的粒子数达到上限，剩余粒子别无去处只能跳进基态——这正是 BEC 的临界条件，第 3 章将把这个"装不下"算成定量判据。

占据数随化学势逼近零的爆发

μ 越接近 0，低能级占据越高，曲线整体"抬升"——这就是降温或加压时系统走向凝聚的统计图像。

化学势的角色问答

光子为什么没有化学势？

光子数不守恒（壁面吸收和发射随时改变总数），没有"固定 N"的约束，μ 退化为 0。这让 BE 分布直接变成普朗克分布——2.3 节的入口。

什么时候 μ 严格为负？

任何粒子数守恒且未凝聚的玻色气体。实验上超冷原子气体的 μ 可以用密度和温度反推出来，是表征体系离凝聚多远的核心参数（第 3 章的临界判据就用它表达）。

本节要点回顾

• 计数即统计：玻色计数（隔板法）与经典计数在权重上差指数倍，这是一切差异的源头。
• 分布函数：n(ε) = 1/(e^{(ε−μ)/kT} − 1)，签名是分母的"−1"。
• 推导路线：状态计数 → 斯特林近似 → 拉格朗日乘子极大化熵，三步完成。
• 化学势边界：μ ≤ 0 是硬约束；μ → 0 对应基态占据宏观化，即凝聚前夜。
• 数值可验：小系统穷举计数能复现分布形状，说明公式不神秘。

有了公式，下一节把它放进三种统计的坐标系，看它何时退化成经典、何时与费米统计分道扬镳。

代入数字：分布函数的三个工作区间

玻色-爱因斯坦分布 n(ε) = 1/(exp((ε−μ)/kT) − 1) 的行为取决于 (ε−μ)/kT 的大小。以光子气体（μ=0）为例，取 kT = 0.025 eV（室温），看几个能量档的占据数，能直观看到"低能模式被疯抢、高能模式近经典"的过渡。

import math
kT = 0.025  # eV，室温
mu = 0.0    # 光子气化学势为零
for eps in [0.001, 0.005, 0.025, 0.25, 2.5]:
    x = (eps - mu) / kT
    n_be = 1 / (math.exp(x) - 1)
    n_mb = math.exp(-x)
    print(f"eps={eps:6.3f} eV  x={x:5.1f}  BE占据={n_be:9.2f}  MB占据={n_mb:8.4f}")
# 低能端 BE 占据远超 MB（红外群聚），高能端两者趋于一致
​

化学势的物理边界

μ 趋近 ε₀ 的过程值得盯紧：温度高于临界温度时，逸出基态的粒子由激发态容纳；温度降到 Tc 以下，激发态"满员"，多余的粒子只能全部灌进基态——这就是凝聚分数 N₀/N = 1 − (T/Tc)^1.5 的来源。化学势被基态能量卡死，是玻色系统独有的饱和机制，费米子没有类似限制，因为它的每个态本来就只容一个粒子，永远不会出现"无处可去只能集体沉底"的剧情。