参数设置
光子和声子中 μ = 0（粒子数非守恒）。有质量的玻色子气体（^87Rb、^4He 等）中 μ < 0，降温时 μ → 0，BEC 开始发生。
分布函数 n(E)
横轴＝(E−μ)/(kT)／纵轴＝占有数 n／蓝实线＝BE分布，橙虚线＝MB分布／黄竖线＝当前观测点／E→μ 时BE发散（BEC迹象）
温度依赖性 n_BE(T)
固定 E 和 μ，在 1〜1000 K 范围扫描 T。低温时 n_BE 急剧增长，展现 BEC 迹象。黄圆＝当前 T。
理论与主要公式
整数自旋的玻色子在热平衡下占据能量为 $E$ 的状态的平均粒子数遵循玻色-爱因斯坦分布。
占有数（$\mu$ 为化学势，$k$ 为玻尔兹曼常数，$T$ 为绝对温度，$g$ 为简并度）：
$$n_{BE}(E) = \frac{g}{\exp\!\left[\dfrac{E - \mu}{kT}\right] - 1}$$
古典极限（$E - \mu \gg kT$）接近麦克斯韦-玻尔兹曼分布：
$$n_{MB}(E) = g\,\exp\!\left[-\dfrac{E - \mu}{kT}\right]$$
热能与典型值：
$$kT = 8.617\times 10^{-5}\,T\ \text{[eV]},\quad kT_{10\text{K}} \approx 0.862\ \text{meV}$$
当 $E \to \mu$ 时，$n_{BE} \to \infty$，大量玻色子凝聚到单一状态，发生 BEC。在 $E - \mu \ll kT$ 区域，$n_{BE} \approx kT/(E - \mu)$。
玻色-爱因斯坦分布模拟器简介
🙋
「玻色-爱因斯坦分布」与费米-狄拉克分布有什么区别呢？
🎓
两者都是量子统计，但适用对象不同。费米-狄拉克适用于电子这样的半整数自旋「费米子」，受泡利不相容原理限制，一个量子态最多只能容纳 1 个粒子。而玻色-爱因斯坦适用于整数自旋的「玻色子」（光子、声子、^4He、冷却的 ^87Rb 等），同一态可容纳任意多个粒子。公式上，玻色-爱因斯坦是 $n_{BE}(E) = 1/(\exp((E-\mu)/kT)-1)$，费米-狄拉克 是 $n_{FD}(E) = 1/(\exp((E-\mu)/kT)+1)$，分母符号的区别导致完全不同的物理行为。
🙋
默认情况下 E=0.001 eV，μ=−0.001 eV，T=10 K 时，n_BE = 0.109。与 MB 的差异达到 10.9%，这个差值很大呢。
🎓
对，这正说明此时已经进入「量子效应显著」的区域。计算一下，(E−μ)/kT = 0.002/0.000862 ≈ 2.32，exp(2.32) ≈ 10.18，所以 BE 分布给出 1/(10.18−1) = 0.109，而 MB 分布为 exp(−2.32) = 0.0982，相差 10.9%。这说明玻色子有「聚集效应」——倾向于占据相同的量子态。当 (E−μ)/kT = 5 时差异只有 0.7%，趋近古典行为；但当 (E−μ)/kT = 0.5 时，差异超过 30%，量子效应非常强；而当 E→μ 时，BE 分布发散到无穷，这就是玻色-爱因斯坦凝聚的数学信号。
🙋
点击「扫描温度」按钮，低温时 n_BE 会急剧增长呢！
🎓
这正是 BEC（玻色-爱因斯坦凝聚）的「迹象」。在这个模拟器中，E 和 μ 固定，只改变 T，所以当 T 降低时 kT 减小，(E−μ)/kT 变大……等等，反而应该使占有数减小。实际上在真实的物理系统中，当温度降低时，化学势 μ 会上升，趋近基态能量。一旦 μ 接近基态，基态的占有数就会爆炸性地增长，大量粒子「凝聚」到基态形成巨观量子态。1995 年，Cornell、Wieman 等人通过激光冷却和蒸发冷却，把 87Rb 原子冷却到 170 nanoKelvin，首次观察到 BEC，随后获得了 2001 年的诺贝尔物理学奖。
🙋
当 μ 接近 0（比如 −0.0001 eV）时，n_BE 会变得极其巨大…
🎓
没错，这就是 BEC 的数学体现。当 x = (E−μ)/kT 变得很小时，$n_{BE} \approx kT/(E - \mu)$，从指数衰减变成反比例衰减，这是一个质的变化。真实实验中，冷原子在光阱和磁陷阱中被困禁，通过蒸发冷却（类似于咖啡吹凉的原理，逐步除去高能粒子）达到 nanoKelvin 或甚至 picoKelvin 的温度。此时有数千到数百万个原子占据同一个量子态，形成「物质波」，可以干涉、衍射，展现出量子力学最奇妙的一面。
常见问题
化学势 μ 是与粒子数守恒相关的拉格朗日乘数，定义为 $\mu = \partial F/\partial N$（F 为自由能）。光子和声子在黑体辐射和晶格振动中可以自由生成和湮灭，粒子数 N 不守恒。因此在热力学平衡的自由能最小化条件下，必须有 μ = 0。这解释了普朗克黑体辐射公式 $u(\nu) \propto \nu^3/(\exp(h\nu/kT)-1)$ 分母中为什么没有额外的常数项。而对于有质量的粒子（如 ^4He、^87Rb 原子），粒子数受守恒律约束，μ 会随温度和密度变化，低温时 μ → E_0（基态能量），BEC 就此发生。
分母符号的差异导致根本不同的物理。FD 为 $f_{FD}(E) = 1/(e^x+1)$，总是在 0〜1 之间（泡利不相容原理），每态最多 1 个粒子。BE 为 $n_{BE}(E) = 1/(e^x-1)$，可取任意正值（包括无穷大），同一态可容纳多个粒子。当 x ≡ (E−μ)/(kT) ≫ 1 时，两者都趋向麦克斯韦-玻尔兹曼分布 $\propto e^{-x}$。但在 x ≪ 1 区域，FD 趋向 f → 1（线性饱和），BE 趋向 n → 1/x（反比例发散）。这个差异本质上反映了对称性的不同：费米子遵循反对称波函数（排斥），玻色子遵循对称波函数（吸引）。
对于三维自由玻色子气体，给定粒子数密度 n 和粒子质量 m，BEC 临界温度为 $T_c = \dfrac{2\pi\hbar^2}{m k_B}\left(\dfrac{n}{\zeta(3/2)}\right)^{2/3}$，其中 $\zeta(3/2) \approx 2.612$ 是黎曼 zeta 函数。以 ^87Rb 原子为例，n ~ 10^14 cm^-3 时，T_c ~ 100 nanoKelvin；对于超流体 ^4He（n ~ 2.2×10^22 cm^-3），T_c ~ 3.13 K，与实验测得的 λ 点 2.17 K 接近。实际实验中，还需要考虑陷阱势的影响和相互作用修正。本模拟器主要展示分布函数的理想行为。
简并度 g 是指具有相同能量 E 的独立量子态的数目。例如，光子有两个偏振自由度（g=2），^4He 原子（自旋 0）在无外场时 g=1，^87Rb 原子的超微细结构中，F=1 状态三重简并（g=3），F=2 状态五重简并（g=5）。计算单个能级的平均占有粒子数时用 g=1；计算能级簇（考虑简并态）的总占有数时用实际的 g 值。本模拟器中 n_BE 直接乘以 g，这是一个线性缩放，不改变分布函数的形状，只改变数值大小。
实际应用
冷原子和玻色-爱因斯坦凝聚：1995 年 Cornell、Wieman 等通过激光冷却和蒸发冷却技术，将 ^87Rb 原子冷却到 170 nanoKelvin，首次实现 BEC，获得 2001 年诺贝尔物理学奖。BEC 是一种巨观量子现象，由单一波函数描述的几千到几百万个原子形成量子流体，表现出零黏性、量子涡旋等奇异性质。光学陷阱中的 BEC 可用于原子激光、物质波干涉、拓扑缺陷研究，是探索强关联量子系统的最重要平台。
普朗克黑体辐射与宇宙微波背景：光子是 μ = 0 的玻色子，能谱密度为 $u(\nu) \propto h\nu^3/(\exp(h\nu/kT)-1)$，这就是普朗克分布——1900 年量子论的起点。宇宙微波背景辐射（CMB）是大爆炸遗留的光子气体，温度约 2.725 K，其能谱与普朗克分布的吻合度在 $10^{-5}$ 以下，是大爆炸宇宙学最强有力的证据。COBE 卫星（1989）的精密测量获得 1974 年诺贝尔物理学奖。
超流体 ^4He 和超导：液态 ^4He 在 2.17 K（λ 点）以下进入超流相，表现为零黏性、爬壁等奇异流体动力学。这是 ^4He（玻色子）的 BEC 现象，通过 Landau 准粒子理论解释。超导电体（如 Nb、NbTi、铜氧化物等）中，电子通过相互作用形成库珀对（两个电子的束缚态），这些对充当有效玻色子，低温时凝聚形成无阻超导态。MgB_2（T_c=39 K）、铜酸化物高温超导体、最近发现的水合物 LaH_10（高压下 T_c~250 K）都是玻色凝聚思想的应用。
光学和激光技术：激光光束是大量光子凝聚到单个模态的相干态，被称为「相干态」。这是 1917 年爱因斯坦提出的诱导发射原理的直接应用——玻色子的「聚集效应」使得一旦有光子被放出，更多同样的光子更容易被放出，形成级联放大。半导体激光、光通信、LIDAR、光学时钟等现代技术的核心都依赖于对玻色子统计的深刻理解。
常见误解与注意事项
最常见的误解是认为 「BEC 就是液化或凝聚，像水蒸气变成水滴」。实际上 BEC 完全不同——这不是实空间的相变，而是「动量空间」（或能量特征态空间）的凝聚。原子们仍然是分散的气体状态，但它们的量子波函数变成了一个单一的、可以用一个波函数描述的宏观量子态。这是一种纯粹的量子现象，在 nanoKelvin 温度下，de Broglie 波长变得很长（~ μm），许多原子的波函数开始重叠，整个系统变成一个「超原子」。
第二个常见误解是试图输入 μ ≥ E 的参数。玻色-爱因斯坦分布 n_BE(E) = 1/(exp((E−μ)/kT)−1) 的物理有效范围仅限于 E > μ。若 E ≤ μ，分母变为零或负数，占有数变成无穷大或负值，失去物理意义。本模拟器通过约束 μ < 0、E > 0 来保证 E−μ > 0。在真实的粒子数守恒系统中，降低温度时 μ 会逐渐上升，直到某个临界温度时 μ 趋近基态能量 E_0，BEC 就此开始。
第三个误解是 「所有玻色子都能发生 BEC」。BEC 是三维或更高维自由玻色气体的现象，受 Mermin-Wagner 定理的限制：在二维系统中，有限温度不存在真正的凝聚，只有 Berezinski-Kosterlitz-Thouless 转变等拓扑相变；一维更不可能。而且，相互作用也会改变 BEC 的性质，需要用 Bogoliubov 理论描述准粒子。本模拟器展示的是理想无相互作用玻色气体的分布函数，实际应用中需要根据具体系统进行修正。