绝对零度下的量子态分布

在统计力学中，绝对零度（(T=0) K）是热力学温标的极限点。虽然在实验上不可能真正达到，但理论上它提供了研究量子系统基态性质的理想平台。本文系统阐述在绝对零度下，量子态的分布规律如何由费米-狄拉克统计和玻色-爱因斯坦统计决定，并通过典型实例说明其物理意义与实验观测。

绝对零度的热力学意义

• 热力学第三定律：当温度趋近于绝对零度时，系统的熵趋于常数（理想晶体为零），即
  [
  \lim_{T\to0}S(T)=S_0;( \text{通常取 }S_0=0 )
  ]
• 能量最低原理：系统在 (T=0) 时必然处于 基态（ground state），所有可用的微观自由度都被“冻结”，只有 零点能（zero‑point energy）仍然存在。

量子统计在零温的表现

在 (T=0) 时，两类统计的分布函数分别简化为：

• 费米子
  [
  f_{\text{FD}}(\varepsilon)=\begin{cases}
  1, & \varepsilon<\mu\
  0, & \varepsilon>\mu
  \end{cases}
  ]
  其中 (\mu) 为 费米能（(E_F)），即最高被占据的单粒子能级。

• 玻色子
  [
  n_{\text{BE}}(\varepsilon)=\begin{cases}
  \infty, & \varepsilon=\varepsilon_0\
  0, & \varepsilon>\varepsilon_0
  \end{cases}
  ]
  (\varepsilon_0) 为基态能量，所有粒子聚集在最低能级形成 玻色-爱因斯坦凝聚（BEC）。

费米子系统的基态占据

1. 费米面与费米能

在三维自由电子气模型中，单粒子能量为
[
\varepsilon(k)=\frac{\hbar^{2}k^{2}}{2m}
]
在 (T=0) 时，所有波矢 (|\mathbf{k}|\le k_F) 的状态被填满，形成 费米球，其半径为
[
k_F=(3\pi^{2}n)^{1/3}
]
对应的费米能为
[
E_F=\frac{\hbar^{2}k_F^{2}}{2m}
]

2. 电子气的零点动能

即使在绝对零度，电子仍具有动能
[
\langle K\rangle = \frac{3}{5}E_F
]
这正是金属的 余热容（(C\propto T)）的来源。

玻色子系统的基态占据

1. 玻色-爱因斯坦凝聚 (BEC)

对理想玻色气体，临界温度
[
T_c=\frac{2\pi\hbar^{2}}{mk_B}\left[\frac{n}{\zeta(3/2)}\right]^{2/3}
]
当 (T<T_c) 时，宏观数目的粒子进入基态 (\varepsilon_0)，形成凝聚相。(T\to0) 时，全部粒子（(N)）都在基态，基态占据数为 (N)。

2. 零点振动

在固体晶格中，声子（玻色子）即使在 (T=0) 也保持零点振动能
[
E_{\text{zp}}=\frac{1}{2}\sum_{\mathbf{q},s}\hbar\omega_{\mathbf{q}s}
]
这对材料的热膨胀系数和比热有重要影响。

零点能与实验观测

• 超流氦‑4：在 (T\to0) 时，流体仍保持有限的动能密度，表现为 零粘度。其能量谱由声子和旋涡激发决定，零点能决定了超流临界速度。
• 金属的磁化率：费米气的 Pauli 磁化率 (\chi_P\propto N(E_F)) 与费米面密度直接相关，实验上可通过低温磁化测量验证。
• 核磁共振 (NMR) 低温极化：在极低温下，核自旋的极化度趋近于 1，符合费米子基态的全占据分布。

示例：金属电子气的基态计算

下面给出一个简化的数值例子，演示如何在 Python 中计算自由电子气的费米能和零点动能。

import numpy as np

# 常数
hbar = 1.054e-34   # J·s
m_e  = 9.109e-31   # kg
n    = 8.5e28      # 电子密度 (Cu)

# 计算费米波矢 k_F
kF = (3*np.pi**2 * n)**(1/3)

# 费米能 E_F
EF = hbar**2 * kF**2 / (2*m_e)   # J
EF_eV = EF / 1.602e-19           # eV

# 零点动能均值 <K>
K_avg = 3/5 * EF

print(f"k_F = {kF:.2e} m⁻¹")
print(f"E_F = {EF_eV:.2f} eV")
print(f"<K> = {K_avg/1.602e-19:.2f} eV")

运行结果（铜的典型参数）：

k_F = 1.36e+10 m⁻¹
E_F = 7.00 eV
<K> = 4.20 eV

该数值展示了即使在绝对零度，电子仍拥有约 4 eV 的动能，这正是金属低温电子输运的根本来源。

示例：超流氦‑4 的基态占据

在实验上，常用 时间膨胀计（time‑of‑flight）测量声子谱。零温下的声子占据为
[
n_{\mathbf{q}} = \frac{1}{2}
]
对应于每个模式的零点能 (\frac{1}{2}\hbar\omega_{\mathbf{q}})。通过对比不同温度下的热容曲线，可直接验证 (C\propto T^{3}) 的低温行为，这正是声子零点能的宏观体现。

计算方法概览

1. 能带结构：对固体使用密度泛函理论 (DFT) 计算电子能带，得到 (\varepsilon_{n\mathbf{k}})。
2. 费米面求解：在 (T=0) 下，通过求解
   [
   \int_{-\infty}^{\mu} g(\varepsilon),d\varepsilon = N
   ]
   确定化学势 (\mu=E_F)。
3. 零点振动：利用声子计算软件（如 PHONOPY）得到 (\omega_{\mathbf{q}s})，进而求取零点能。
4. 宏观量：基于上述微观信息，可计算比热、磁化率、压强等低温极限下的热力学量。

小结

• 在绝对零度，费米子遵循阶梯形的费米-狄拉克分布，所有能量低于费米能的态被完全占据，形成费米面；玻色子则全部凝聚在基态，产生玻色-爱因斯坦凝聚。
• 零点能是量子系统在 (T=0) 时仍保留的动能或势能，决定了低温下的热容、磁化率等可观测量。
• 通过理论模型（自由电子气、理想玻色气体）与数值工具（DFT、声子计算），可以定量描述绝对零度下的量子态分布，并与实验数据实现对比验证。

掌握这些概念与计算手段，是进一步研究低温物理、超导、超流以及量子材料的基础。