玻色—爱因斯坦统计

玻色—爱因斯坦统计是统计力学中用于对全同玻色子的状态进行计数并确定其平衡占据数的理论框架。与经典粒子计数不同，它将粒子视为不可分辨的，并允许单粒子态的占据数为任意非负整数。其特有的分布描述了处于热力学平衡的理想玻色子系统，也是热辐射和玻色—爱因斯坦凝聚的理论基础。这些统计规则普遍适用于玻色子，但常见的占据数公式以粒子间无相互作用或可采用适当的独立模式描述为前提。(damtp.cam.ac.uk)

历史发展

1924年，萨特延德拉·纳特·玻色提出了一种新的光量子计数方法，重现了黑体辐射的光谱。阿尔伯特·爱因斯坦在1924年和1925年发表的论文中，将这一方法推广到实物粒子。这一推广预言，当温度降至足够低时，宏观数量的粒子可以占据最低能态。这成为玻色—爱因斯坦凝聚的理论基础。(nobelprize.org)

1995年，实验首次在稀薄原子气体中实现了凝聚。埃里克·康奈尔、卡尔·威曼和沃尔夫冈·凯特勒因在稀薄碱金属原子气体中实现凝聚并研究凝聚态的性质，获得了2001年诺贝尔物理学奖。这些实验以尤为直接的方式，实现了数十年前理论所预言的量子统计行为。(nist.gov)

不可分辨性与状态计数

在量子力学中，交换两个全同玻色子不会改变其多粒子波函数。因此，玻色子态具有交换对称性。玻色子的自旋角动量量子数为整数，其统计行为与费米子不同：费米子态在交换下具有反对称性，并遵循泡利不相容原理。(ocw.mit.edu)

玻色子的多粒子组态由占据数 $n_i=0,1,2,\ldots$ 指定，而不是通过给粒子赋予可区分的标号来指定。若某一能级包含 $g_i$ 个不同的单粒子态和 $N_i$ 个玻色子，则可能的占据组态数为

W_i=\frac{(N_i+g_i-1)!}{N_i!\,(g_i-1)!}.

在各能级粒子数给定时，总组态数就是这些因子的乘积。这种计数与经典计数不同：交换粒子标号不会产生新的物理组态。多个粒子可以占据同一状态，但并非必须如此；多重占据本身也不意味着粒子之间存在吸引力。(ocw.mit.edu)

平衡占据数分布

对于彼此独立的玻色子态，能量为 $\epsilon_i$ 的状态的平均占据数为

\bar n_i= \frac{1}{\exp[(\epsilon_i-\mu)/(k_{\mathrm B}T)]-1}.

其中， $T$ 为绝对温度， $k_{\mathrm B}$ 为玻尔兹曼常数， $\mu$ 为化学势。该公式给出的是平均粒子数，而不是对各能级归一化的概率。若某能级的简并度为 $g_i$ ，则该能级的平均总占据数为 $g_i\bar n_i$ 。对各态的占据数求和可得平均总粒子数；对 $\epsilon_i\bar n_i$ 求和可得平均总能量。(damtp.cam.ac.uk)

这一分布可以从巨正则系综自然地导出。令 $\beta=1/(k_{\mathrm B}T)$ ，则单个状态的巨配分函数为几何级数

\Xi_i=\sum_{n=0}^{\infty}e^{-\beta n(\epsilon_i-\mu)} =\frac{1}{1-e^{-\beta(\epsilon_i-\mu)}}.

为保证级数收敛，在有限的理想系统中， $\mu$ 必须低于最低单粒子能量。当 $\mu$ 从下方趋近该能量时，基态的占据数会变得很大。(homepages.ucl.ac.uk)

占据数涨落与经典极限

对于一个独立模式，令 $q=e^{-\beta(\epsilon-\mu)}$ ，则占据数的概率分布为

P(n)=(1-q)q^n,\qquad n=0,1,2,\ldots.

其方差满足

\operatorname{Var}(n)=\bar n(1+\bar n).

与均值相同的泊松分布相比，二次项使涨落更强。这一结果适用于处于热平衡的独立模式；额外的约束或关联可能改变整个系统的涨落。(ocw.mit.edu)

当 $\exp[(\epsilon-\mu)/(k_{\mathrm B}T)]\gg1$ 时，分母中的减一项可以忽略，从而得到经典的指数形式占据规律。因此，在稀薄、低占据数极限下，玻色—爱因斯坦统计趋于麦克斯韦—玻尔兹曼统计。相比之下，费米—狄拉克统计的分母中为加一项，并将单个状态的平均占据数限制为不超过一。(damtp.cam.ac.uk)

玻色—爱因斯坦凝聚

对于仅有一种内部状态的均匀、三维、非相对论性玻色子理想气体，在给定温度下，当化学势趋近基态能量时，激发态所能容纳的粒子数密度存在有限上限。此时，额外的粒子会以宏观数量积聚在基态，形成玻色—爱因斯坦凝聚态。(ocw.mit.edu)

对于粒子数密度 $\rho=N/V$ 和粒子质量 $m$ ，有

T_c=\frac{2\pi\hbar^2}{mk_{\mathrm B}} \left[\frac{\rho}{\zeta(3/2)}\right]^{2/3}, \qquad \frac{N_0}{N}=1-\left(\frac{T}{T_c}\right)^{3/2} \quad(T<T_c).

其中， $\hbar=h/(2\pi)$ ， $h$ 为普朗克常数， $\zeta$ 为黎曼ζ函数。这些关系描述的是热力学极限下的理想均匀气体，并非适用于所有受俘获或存在相互作用的系统。(web.mit.edu)

辐射、振动与适用范围

对于通常处于平衡态的热辐射，光子可以被产生和吸收，因此其数目不守恒，化学势为零。其占据数变为

\bar n(\omega)=\frac{1}{e^{\hbar\omega/(k_{\mathrm B}T)}-1},

即普朗克分布。它确定的是热平衡下各模式的占据数，而不要求总光子数固定。(farside.ph.utexas.edu)

同样的占据因子也描述固体中的简谐声子模式。(oden.utexas.edu) 相互作用需要进一步处理：玻色子的交换对称性仍然成立，但存在相互作用的系统未必具有理想气体的占据数分布。同样，凝聚与超流性相互关联，但并不等同；仅靠无相互作用的凝聚态模型，无法完整描述超流行为。(damtp.cam.ac.uk)