Learning objectives / 学习目标
- 用 occupation-number description 解释不可分辨性为何改变微观态计数，而不是给经典粒子补一个标签修正。
- 从单粒子能级的 grand partition function 推导 Bose–Einstein 与 Fermi–Dirac 平均占据数。
- 解释 Pauli exclusion 对费米占据、费米面与低温激发的约束。
- 用 fugacity \(z\)、thermal de Broglie wavelength \(\lambda_{th}\) 和 \(n\lambda_{th}^3\) 判断经典极限。
- 正确陈述 photon、phonon、守恒 massive boson 与 fermion 的 chemical-potential 条件。
- 对量子统计计算执行能量零点、简并度、粒子数守恒、量纲与极限检查。
Prerequisites / 先修知识
- 巨正则系综、配分函数与 Boltzmann factor。
- 经典理想气体与能量均分的失效。
- 单粒子量子态、波函数交换对称性、spin 与 Pauli exclusion 的基础概念。
- 指数函数、几何级数、组合数和三维自由粒子能量 \(\varepsilon=p^2/(2m)\)。
Scope and conventions / 范围与约定
- 讨论 thermodynamic equilibrium 下的理想或弱相互作用量子气体；分布式本身不包含强关联、超流动力学或相对论修正。
- 单粒子态标为 \(i\)，能量为 \(\varepsilon_i\)，occupation number 为 \(n_i\)；总粒子数 \(N=\sum_i n_i\)。
- 温度 \(T>0\)，\(\beta=(k_BT)^{-1}\)，fugacity \(z=e^{\beta\mu}\)。能量零点可以平移，故玻色约束应写成 \(\mu\le\varepsilon_0\)，其中 \(\varepsilon_0\) 是最低单粒子能量。
- Boson 单态允许 \(n_i=0,1,2,\ldots\)；fermion 的每个完整单粒子态只允许 \(n_i=0,1\)。若一个空间能级含 spin degeneracy，必须把每个 spin state 分开计数或显式乘简并度。
- 无相互作用、非相对论、均匀三维气体的经典性用 \(n\lambda_{th}^3/g_s\ll1\) 判断；只说“高温”并不充分。
Notation and modeling contract / 符号与建模契约
| Symbol | Meaning | SI unit |
|---|---|---|
| \(\varepsilon_i,\varepsilon_0\) | single-particle energy, ground-state energy | J |
| \(n_i,\bar n_i\) | occupation number, mean occupation | 1 |
| \(N,n=N/V\) | particle number, number density | 1, \(\mathrm{m^{-3}}\) |
| \(\beta\) | inverse thermal energy, \(1/(k_BT)\) | \(\mathrm{J^{-1}}\) |
| \(\mu\) | chemical potential | J |
| \(z\) | fugacity, \(e^{\beta\mu}\) | 1 |
| \(g_s\) | internal-state degeneracy included in density of states | 1 |
| \(\lambda_{th}\) | \(h/\sqrt{2\pi m k_BT}\) | m |
| \(f_{BE},f_{FD},f_{MB}\) | mean occupation functions | 1 |
Knowledge 1 — Indistinguishability and occupation numbers / 不可分辨性与占据数
Knowledge ID: P03-U07-L01-K01
1.1 相同粒子的交换不是一个新实验结果
Classical labeled particles 可由轨迹标签区分；identical quantum particles 的标签不是 observable。交换两粒子后，多体波函数对 bosons 保持对称，对 fermions 改变符号，但所有概率不变。自然的微观态语言因此不是“粒子 1 在哪、粒子 2 在哪”，而是每个单粒子态被占据多少次：
\[\lvert n_1,n_2,\ldots\rangle, \qquad N=\sum_i n_i, \qquad E=\sum_i n_i\varepsilon_i.\]
图 1：occupation-number language 删除不可观测的粒子标签；右图两个圆点只表示同一态有两个 bosons。
1.2 允许占据决定组合计数
把 \(N\) 个相同粒子放入 \(g\) 个等价单粒子态：
\[W_{BE}=\binom{N+g-1}{N}, \qquad W_{FD}=\binom{g}{N}\quad(N\le g).\]
第一式允许任意 pile-up；第二式要求每态至多一个。经典可分辨粒子的 \(g^N\) 是不同问题，不能与量子计数混用。
图 2：Bosons 可以共享完整单粒子态；fermions 的 Pauli restriction 作用于包括 spin quantum number 在内的完整态。
Worked Example 1 — 三种计数为何不同
已知： \(N=2\) 个粒子、\(g=3\) 个等价单粒子态。分别求 bosons、spin-polarized fermions 与可分辨经典粒子的微观态数。
系统与模型： 能量固定且只比较状态计数；fermions 已 spin-polarized，所以每个空间态只有一个允许 internal state。
方程： 使用上面的 combinations；经典粒子保留标签。
\[W_{BE}=\binom{4}{2}=6, \qquad W_{FD}=\binom{3}{2}=3, \qquad W_{labelled}=3^2=9.\]
因此
\[\boxed{W_{BE}:W_{FD}:W_{labelled}=6:3:9.}\]
检查： 三者皆无量纲；若 \(N>g\)，spin-polarized Fermi count 必为零，而 Bose count 仍有限。差异来自可容许多体态，而不是不同的粒子间作用力。
Knowledge 2 — Bose–Einstein distribution / 玻色—爱因斯坦分布
Knowledge ID: P03-U07-L01-K02
2.1 单态 grand partition function 给出平均占据
对一个 bosonic state，\(n_i=0,1,2,\ldots\)。令 \(x_i=e^{-\beta(\varepsilon_i-\mu)}\)，收敛要求 \(x_i<1\)，于是
\[\mathcal Z_i^{BE}=\sum_{n_i=0}^{\infty}x_i^{n_i} =\frac{1}{1-x_i}.\]
平均占据为
\[\bar n_i=x_i\frac{\partial}{\partial x_i}\ln\mathcal Z_i^{BE} =\boxed{\frac{1}{e^{\beta(\varepsilon_i-\mu)}-1}}.\]
分母的 minus sign 允许低能态出现宏观占据；但只有在指定 density of states、维度与 particle-number constraint 后，才能判断是否真的发生 condensation。
图 3：Bose occupation 在 \(x\to0^+\) 时发散；\(x\le0\) 会使单态几何级数不收敛。
2.2 简并态要乘 density of states，不是改分布
若能量壳 \([\varepsilon,\varepsilon+d\varepsilon]\) 含 \(D(\varepsilon)d\varepsilon\) 个单粒子态，粒子数为
\[dN=D(\varepsilon)f_{BE}(\varepsilon)d\varepsilon.\]
\(f_{BE}\) 是每态平均占据，\(D\) 是态数密度；二者单位和物理角色不同。连续积分常需把离散 ground state 单独留下，否则正是可能宏观占据的那一项会被 measure-zero approximation 丢掉。
Worked Example 2 — 一个 bosonic mode 的平均占据
已知： \(T=1000\ \mathrm K\)，\(\varepsilon-\mu=0.150\ \mathrm{eV}\)。求该完整单粒子态的 Bose 平均占据。
系统与能量约定： grand-canonical ideal boson mode；给定的是能差，故能量零点不影响答案。
常数： \(k_B=8.617333262\times10^{-5}\ \mathrm{eV\,K^{-1}}\)。
\[x=\frac{0.150}{(8.617333262\times10^{-5})(1000)}=1.7407,\]
\[\bar n=\frac{1}{e^{1.7407}-1}=0.2127.\]
所以
\[\boxed{\bar n=0.213\ \text{particle per state}.}\]
检查： \(x\) 无量纲且为正，满足收敛条件；平均占据可以小于一且不必为整数。若 \(\varepsilon-\mu\) 增大，结果单调减小。
Knowledge 3 — Fermi–Dirac distribution / 费米—狄拉克分布
Knowledge ID: P03-U07-L01-K03
3.1 Pauli restriction 把几何级数截成两项
对一个 fermionic state，\(n_i=0,1\)，故
\[\mathcal Z_i^{FD}=1+e^{-\beta(\varepsilon_i-\mu)}.\]
同样取 logarithmic derivative，
\[\boxed{ \bar n_i=f_{FD}(\varepsilon_i) =\frac{1}{e^{\beta(\varepsilon_i-\mu)}+1} }.\]
因此 \(0\le f_{FD}\le1\)。在任何 \(T>0\)，\(f_{FD}(\mu)=1/2\)；在 \(T\to0\) 时它成为 step function，低于 Fermi energy 的态填满，高于它的态空出。
图 4：有限温度只在宽度约数个 \(k_BT\) 的窗口内把 Fermi edge smoothed；深能级仍近乎填满。
3.2 “每态至多一个”不等于“每个能量至多一个”
一个能量值可能有 orbital、momentum 与 spin degeneracy。总粒子数应写成
\[N=\sum_i f_{FD}(\varepsilon_i) \quad\text{或}\quad N=\int D(\varepsilon)f_{FD}(\varepsilon)d\varepsilon,\]
其中 \(i\) 遍历完整单粒子态。若 \(D\) 已含 \(g_s=2\)，不能再额外乘二。
图 5：低温 Fermi gas 的热响应来自 Fermi surface 附近的薄能量壳，而不是所有电子同时获得 \(k_BT\)。
Worked Example 3 — Fermi edge 两侧的对称占据
已知： 某 fermionic state 的能量分别为 \(\varepsilon_- =\mu-k_BT\ln3\) 与 \(\varepsilon_+=\mu+k_BT\ln3\)。求两态平均占据。
系统与模型： equilibrium ideal fermions；两个 state 完整包含 internal quantum numbers。
方程： 直接用 Fermi–Dirac distribution。
\[f_- =\frac{1}{e^{-\ln3}+1}=\frac{3}{4}=0.750,\]
\[f_+ =\frac{1}{e^{\ln3}+1}=\frac{1}{4}=0.250.\]
因此
\[\boxed{f_-+f_+=1.000.}\]
检查： 两能量关于 \(\mu\) 对称，Fermi function 满足 \(f(\mu-\Delta)+f(\mu+\Delta)=1\)；结果均在 \([0,1]\) 内。这里的 0.750、0.250 是 ensemble averages，不表示单次测量出现 fractional fermion。
Knowledge 4 — Classical limit of quantum statistics / 量子统计的经典极限
Knowledge ID: P03-U07-L01-K04
4.1 稀占据让 plus/minus one 都变得不重要
当
\[e^{\beta(\varepsilon-\mu)}\gg1,\]
两种分布都退化为
\[f_{BE}\approx f_{FD}\approx e^{-\beta(\varepsilon-\mu)}=z e^{-\beta\varepsilon} \equiv f_{MB}.\]
这要求 relevant states 的平均占据远小于一；不是把 \(\hbar\) 设为零，也不是只看单个能级间隔。
图 6：在 \(x\gg1\) 的 dilute tail，BE 与 FD 的差异相对于 leading Boltzmann factor 消失。
4.2 相空间重叠给出可操作判据
对非相对论三维自由粒子，
\[\lambda_{th}=\frac{h}{\sqrt{2\pi m k_BT}}.\]
若 internal degeneracy 已记作 \(g_s\)，经典极限通常要求
\[\boxed{\frac{n\lambda_{th}^3}{g_s}\ll1.}\]
\(n\lambda_{th}^3\) 衡量一个 thermal wavepacket volume 内平均有多少粒子；升温、降密度或增大质量都会减弱波包重叠。
图 7：Degeneracy parameter 是连续 crossover 指标；\(\sim1\) 不是对所有体系都相同的精确阈值。
图 8：判断经典性必须同时给 mass、temperature、density 与 internal degeneracy，再选择分布。
Worked Example 4 — 氦气在室温是否量子简并
已知： \(^4\mathrm{He}\) 原子质量 \(m=6.646\times10^{-27}\ \mathrm{kg}\)，\(T=300\ \mathrm K\)，\(n=2.00\times10^{25}\ \mathrm{m^{-3}}\)，取 \(g_s=1\)。
系统与模型： 均匀、非相对论、弱相互作用的理想化气体；只判断统计 regime。
方程与计算： 使用 exact SI constants \(h\) 与 \(k_B\)。
\[\lambda_{th} =\frac{6.62607015\times10^{-34}} {\sqrt{2\pi(6.646\times10^{-27})(1.380649\times10^{-23})(300)}} =5.04\times10^{-11}\ \mathrm m.\]
\[n\lambda_{th}^3 =(2.00\times10^{25})(5.04\times10^{-11})^3 =2.56\times10^{-6}.\]
所以
\[\boxed{n\lambda_{th}^3\ll1,\quad \text{Maxwell--Boltzmann regime}.}\]
检查： 参数无量纲；波长约 \(0.0504\ \mathrm{nm}\)，远小于平均粒子间距 \(n^{-1/3}\approx3.68\ \mathrm{nm}\)。真实 dense helium 的相互作用仍须单独检查，但本给定密度下交换效应很小。
Knowledge 5 — Chemical potential for quantum gases / 量子气体的化学势
Knowledge ID: P03-U07-L01-K05
5.1 \(\mu\) 同时控制粒子数与分布可定义域
Grand-canonical weight 为
\[P\propto e^{-\beta(E-\mu N)}.\]
对 particle-number-conserving gas，\(\mu\) 由约束
\[N=\sum_i\frac{1}{e^{\beta(\varepsilon_i-\mu)}\mp1}\]
隐式确定，上方 minus 对 bosons、plus 对 fermions。\(\mu\) 不是任意 fit parameter，也不能在改变能量零点时保持数值不变。
对 bosons，所有单态几何级数收敛要求
\[\boxed{\mu\le\varepsilon_0.}\]
有限体系且 \(T>0\) 时通常是严格小于；thermodynamic-limit condensation below \(T_c\) 常把 \(\mu\) pin 到 \(\varepsilon_0\)。若约定 \(\varepsilon_0=0\)，才可简写 \(\mu\le0\)。Fermi gas 没有同样的上界；在 \(T=0\)，\(\mu=\varepsilon_F\)。
图 9：Bosonic chemical potential 的物理条件是相对于最低单粒子能量的上界，而不是绝对的零。
5.2 不守恒 quasiparticles 在平衡时通常有 \(\mu=0\)
在 cavity radiation 或 harmonic crystal 中，photons/phonons 可由物质吸收、产生，平衡约束不固定其数目，因此
\[\mu_\gamma=0, \qquad \mu_{ph}=0.\]
这是 equilibrium、number-nonconserving 条件的结果。若外部 pumping、近似 conserved quasiparticle number 或 separate reservoirs 使粒子数在观测时间内受约束，则可以出现 nonzero effective chemical potential；不可把 \(\mu=0\) 当作所有光场或晶格振动的运动学恒等式。
图 10：chemical potential 的首要决策是“该粒子数是否在平衡约束中守恒”，之后才选择 Bose 或 Fermi domain。
Lesson synthesis / 全课整合
量子统计的逻辑链是：相同粒子没有可观测标签 → 用 occupation numbers 记录多体态 → exchange symmetry 决定允许占据 → grand-canonical sum 给出 BE 或 FD distribution → \(N\) constraint 决定 \(\mu\) → 稀占据极限回到 MB。分布中的 plus/minus sign 不是记忆装饰，而是“单态可否重复占据”的代数痕迹。
建模流程： 明确粒子与完整单粒子态 → 判定 boson/fermion 与 internal degeneracy → 判断 particle number 是否守恒 → 指定能量零点与 \(\varepsilon_0\) → 由 \(N\) 求 \(\mu\) → 乘 density of states 积分 → 检查 \(n\lambda_{th}^3/g_s\)、分布范围和极限。
Formula-and-decision table / 公式—决策表
| Task | Relation | Use only when | Key caution |
|---|---|---|---|
| Bose occupation | \(f_{BE}=1/(e^{\beta(\varepsilon-\mu)}-1)\) | equilibrium ideal bosonic states | require \(\mu\le\varepsilon_0\) |
| Fermi occupation | \(f_{FD}=1/(e^{\beta(\varepsilon-\mu)}+1)\) | equilibrium ideal fermionic states | one per complete state, not per energy |
| classical occupation | \(f_{MB}=ze^{-\beta\varepsilon}\) | sparse occupation | verify degeneracy parameter |
| thermal wavelength | \(\lambda_{th}=h/\sqrt{2\pi m k_BT}\) | nonrelativistic free particles | convention contains \(2\pi\) |
| classicality test | \(n\lambda_{th}^3/g_s\ll1\) | homogeneous 3D ideal-gas estimate | crossover, not universal sharp boundary |
| boson chemical-potential cap | \(\mu\le\varepsilon_0\) | stable equilibrium Bose gas | \(\mu\le0\) only if \(\varepsilon_0=0\) |
| unconserved quanta | \(\mu=0\) | equilibrium photons or phonons with free creation/annihilation | driven systems may differ |
Model boundaries / 模型边界
- BE/FD distributions assume equilibrium quasifree occupations; interactions can alter spectra, lifetimes and equation of state even when exchange statistics remains fixed.
- Spin–statistics connection is relativistic quantum-field-theory content; this lesson uses its consequence rather than proving it.
- \(n\lambda_{th}^3/g_s\) is a uniform nonrelativistic three-dimensional criterion; traps、lattices、lower dimensions 与 relativistic gases require adapted density of states.
- \(\mu\le\varepsilon_0\) applies to stable ideal bosons; condensed interacting systems use a thermodynamic chemical potential that also contains interaction shifts.
- Grand-canonical number fluctuations of the ideal condensate can be anomalous; canonical treatment may be necessary for finite isolated systems.
- Photon/phonon \(\mu=0\) presumes full thermal and chemical equilibrium with a reservoir able to create and annihilate them.
Misconception clinic / 易错概念诊所
- “不可分辨意味着两个粒子在同一位置。” 错；它表示交换标签不产生新物理态，与空间重叠不是同一句话。
- “Pauli principle 说同一能量只能有一个电子。” 错；同能量可有多个不同完整态，每个完整态最多一个 fermion。
- “Bose occupation 发散说明任意 boson gas 必凝聚。” 错；凝聚还取决于密度、维度、density of states 与热力学极限。
- “高温必定经典。” 不完整；还要比较密度、质量和 \(g_s\)，即检查 \(n\lambda_{th}^3/g_s\)。
- “Boson 的 \(\mu\) 永远小于零。” 错；不变量是 \(\mu\le\varepsilon_0\)，能量零点平移会改变两者数值。
- “平均占据 0.25 表示四分之一个粒子。” 错；它是 ensemble probability/average，单次 fermion occupation 仍为 0 或 1。
Practice and retrieval / 练习与检索
1. Occupation state（K01）
写出三个 bosons 分布在两个单粒子态的所有 occupation-number states，并给出数量。
Answer
状态为 \(\lvert3,0\rangle\)、\(\lvert2,1\rangle\)、\(\lvert1,2\rangle\)、\(\lvert0,3\rangle\)，共 \(\binom{4}{3}=4\) 个；交换粒子标签不增加状态。
2. Pauli count（K01, K03）
四个 spin-polarized fermions 能否放进三个完整单粒子态？若不能，微观态数是多少？
Answer
不能；每态至多一个，\(N=4>g=3\)，故 \(W_{FD}=0\)。
3. Bose occupation（K02）
若 \(\varepsilon-\mu=k_BT\ln2\)，求 \(f_{BE}\)。
Answer
\(f_{BE}=1/(e^{\ln2}-1)=1\)。这是每态平均占据，不要求每次测量均为一个。
4. Fermi midpoint（K03）
证明任何 \(T>0\) 都有 \(f_{FD}(\varepsilon=\mu)=1/2\)。
Answer
指数为零，故 \(f_{FD}=1/(e^0+1)=1/2\)；该结论与 temperature 数值无关，但 transition width 随 \(k_BT\) 改变。
5. Internal degeneracy（K01, K03）
一个 spatial orbital 有两个 spin states。它最多能容纳多少电子？这是否违反 Pauli principle？
Answer
最多两个，分别占据 spin-up 与 spin-down 两个不同完整单粒子态；不违反 Pauli principle。
6. Classical tail（K04）
当 \(x=\beta(\varepsilon-\mu)=6.00\) 时，估算 \(f_{MB}\)；说明 BE 与 FD 为何都接近它。
Answer
\(f_{MB}=e^{-6}=2.48\times10^{-3}\ll1\)。此时分母的 \(e^x\) 远大于 \(1\)，所以 plus/minus one 的相对影响很小。
7. Temperature scaling（K04）
固定 \(n,m,g_s\)，温度升高四倍后 \(n\lambda_{th}^3/g_s\) 变为原来的多少？
Answer
\(\lambda_{th}\propto T^{-1/2}\)，故参数 \(\propto T^{-3/2}\)；升温四倍变为 \(4^{-3/2}=1/8\)。
8. Boson chemical potential（K05）
最低单粒子能量是 \(\varepsilon_0=2.00\ \mathrm{meV}\)。判断 \(\mu=2.10\ \mathrm{meV}\) 是否允许。
Answer
不允许；\(\mu>\varepsilon_0\) 会使最低态 Bose denominator 为负且 grand sum 不收敛。允许条件是 \(\mu\le2.00\ \mathrm{meV}\)。
9. Energy-zero shift（K05）
所有 \(\varepsilon_i\) 与 \(\mu\) 同时加常数 \(C\)，occupation 是否变化？
Answer
不变，因为分布只依赖 \(\varepsilon_i-\mu\)。这也说明“\(\mu\) 的绝对正负”没有脱离能量零点的独立意义。
10. Photon exception（K05）
为什么 blackbody photon gas 在平衡时取 \(\mu_\gamma=0\)？给出一项模型边界。
Answer
腔壁可吸收和发射 photons，平衡约束不固定 photon number，因此其 chemical potential 为零。若光场被 pumped 或 photon number 在观测时间近似守恒，就要重新建模。
Exit ticket / 课末检验
- 从单态允许的 \(n_i\) 值分别写出 \(\mathcal Z_i^{BE}\) 与 \(\mathcal Z_i^{FD}\)，并解释分母符号的物理来源。
- 给定 \(m,n,T,g_s\) 与 energy zero，写出选择 MB/BE/FD、确定 \(\mu\) 并检查结果的完整流程。
References and source notes / 参考资料与来源说明
- MIT OCW 5.62, Lecture 8 — Boltzmann, Fermi–Dirac, and Bose–Einstein Statistics — indistinguishability、quantum-state counting 与三种统计的直接本科讲义。
- MIT OCW 2.57, Complete Lecture Notes, Lecture 10 — grand-canonical weight、BE/FD occupation 与 common Boltzmann limit。
- UC San Diego Physics 140B, Chapter 5 — Quantum Statistics — quantum grand ensemble、classical limit、photon/Bose/Fermi gas 的课程结构与约定。
- University of Maryland, Gould and Tobochnik, Statistical and Thermal Physics — density of states、ideal Fermi/Bose gas、chemical potential 与低温边界。
- NIST, CODATA Recommended Values of the Fundamental Physical Constants: 2022 — \(h\)、\(k_B\) 与 SI numerical constants。
以上链接均于 2026-08-25 实际打开并核对到具体页面或 PDF。正文按统一的 \(\varepsilon_i-\mu\) convention 重新组织；全部图依据本文方程原创，不复制来源图像。
Source-to-content map
| Source | Used for | Knowledge ID / Section | Model boundary or convention |
|---|---|---|---|
| MIT 5.62 Lecture 8 | indistinguishability and counting | K01 | occupation-state rather than particle-label language |
| MIT 2.57 Lecture 10 | BE/FD distributions and MB limit | K02–K04 | equilibrium noninteracting modes |
| UCSD Physics 140B Chapter 5 | fugacity, quantum gases, density of states | K02–K05 | grand-canonical convention |
| UMD statistical physics text | ideal Bose/Fermi gas and chemical potential | K03–K05 | homogeneous thermodynamic limit |
| NIST CODATA 2022 | exact constants in Worked Examples | K04 | current SI definitions |
复习与拓展
- 下一课把 distributions 与 photon、phonon、electron density of states 结合，导出 Planck、Debye、Fermi 与 BEC results。
- Condensed-matter physics 将把 ideal occupations 改写为 quasiparticle spectra，并讨论 interactions 与 collective modes。
- Quantum field theory 将证明 spin–statistics connection，并用 creation/annihilation operators 系统表达 occupation states。
- Nonequilibrium statistics 将讨论 driven photons、effective chemical potential 与不满足平衡分布的情形。