能带是从何而来的？晶体电子结构与能带理论的物理起源及交互式Band Structure Lab 玩具模型程序演示

一、引言

能带结构是固体物理的核心概念之一。能带是由晶体周期性和电子 Hamiltonian 的本征值结构产生。对于从事第一性原理计算（如 VASP）或角分辨光电子能谱（ARPES）实验的研究者而言，理解"电子云如何一步步演变为能带"是建立物理直觉的基础。本文从密度泛函理论（DFT）和紧束缚（Tight-Binding, TB）模型两个层面，系统阐述能带的起源，并介绍一个面向初学者的交互式教学程序——Interactive Band Structure Lab。

二、从孤立原子到周期性晶体

2.1 孤立原子的分立能级

孤立原子中的电子处于分立的能级。以碳原子为例，其价电子占据 2s 和 2p 轨道，能级结构为离散谱。当两个碳原子靠近时，原子轨道发生量子力学耦合，原本简并的能级分裂为成键态与反键态。

2.2 一维原子链的能级展宽

将原子数目扩展至 N 个并排列成周期性一维链，每个孤立原子的分立能级将展宽为包含 N 个准连续能级的能带。设晶格常数为 a，引入周期性边界条件后，波矢 k 的取值为 k = 2πn / Na（n = 0, 1, ..., N−1）。每个 k 对应一个量子态，N 个态构成一条能带。

2.3 Bloch 定理

晶体具有平移周期性，势场满足 V(r + R) = V(r)。根据 Bloch 定理，晶体中的单电子波函数可写为

其中 n 为能带指标，k 为晶体动量，u_nk(r) 具有晶格周期性。因此只需在第一布里渊区内求解。

三、密度泛函理论中的能带起源

3.1 电子密度与 Hohenberg–Kohn 定理

真实晶体包含大量电子，直接求解多电子波函数在计算上不可行。DFT 的核心思想由 Hohenberg–Kohn 定理确立：基态电子密度原则上包含体系的全部基态信息。

3.2 Kohn–Sham 方程

Kohn–Sham 方法将多电子问题映射为一组虚构的非相互作用电子在有效势场中的运动，其单粒子方程为

左边第一项为动能算符，第二项为有效势，右边 ε_i 为 Kohn–Sham 本征值。能带即来源于这些本征值。

3.3 有效势的构成

有效势由三部分组成：

其中 V_ion 为离子实产生的势，Hartree 势描述电子间的经典库仑排斥，V_xc 为交换关联势（如 PBE 泛函），是 DFT 中近似最多的部分。

3.4 自洽场（SCF）循环

有效势依赖于电子密度，而电子密度又由 Kohn–Sham 波函数决定，形成如下自洽循环：

程序迭代更新电子密度直至收敛，最终得到收敛的电子密度以及本征函数和本征值。

3.5 能带的本质

在每个 k 点，Kohn–Sham 方程给出一组本征值。沿高对称路径将这些离散点连接，即得到能带结构图。能带不是"把电子云画出来"，而是周期晶体中 Kohn–Sham Hamiltonian 在不同 k 点对角化后的本征能量。

注意：Kohn–Sham 本征值并不严格等于真实多体体系的电子激发能。对于普通金属、半导体等弱关联体系，DFT 能带是良好的近似；但对于磁性体系或重费米子体系，需引入自能修正（如 DMFT）。

3.6 从 DFT 到 TB：Wannier 函数

DFT 计算得到的是离散的能带。通过 Wannier90 等工具，可将能带反变换为一组局域轨道及其跃迁积分：

由此得到紧束缚 Hamiltonian，再对角化即可复现能带。这一流程为：

四、紧束缚模型：理解能带的最简图像

4.1 单原子链

考虑一维单原子链，每个原子提供一个 s 轨道，最近邻跃迁积分为 t。在周期边界条件下，Hamiltonian 为

其解析解为

当 t > 0 时，k = 0（Γ 点）对应能量最低的成键态，k = π/a（X 点）对应能量最高的反键态。增大 t 直接正比增大带宽（约 4|t|）。

4.2 双原子链与带隙

若原胞包含两个不等价原子 A 和 B，其 Hamiltonian 为 2×2 矩阵：

H(k) = ( ε_A t₁+t₂e^−ika / t₁+t₂e^ika ε_B )

对角化后得到

当 ε_A ≠ ε_B 时，布里渊区边界处打开带隙。若 ε_A = ε_B 且 t₁ = t₂，带隙消失。此时体系获得额外的平移对称性，A/B 双原胞描述可以重新写成单原子原胞，因此原来的两条子带发生折叠关系并在相应位置接触；若恢复最小原胞描述，则对应于单原子链的一条能带。

4.3 二维石墨烯的 Dirac 锥

对于蜂窝状晶格（如石墨烯），每个原胞包含两个碳原子。当 ε_A = ε_B = 0 时，能带在 K 点线性交叉，形成 Dirac 锥：

其中 f(k) 为结构因子。在 Dirac 点附近，色散关系为

若打破子晶格对称性（如 ε_A ≠ ε_B），Dirac 点处打开带隙，系统从半金属转变为半导体。

4.4 三维简立方

三维简立方 s 轨道的色散关系为

E(k) = ε₀ + 2t[cos(k_xa) + cos(k_ya) + cos(k_za)]

Γ 点取最小值，R 点取最大值。

五、Interactive Band Structure Lab 程序介绍

5.1 设计思想

该程序是一个面向初学者的晶体电子结构物理沙盒，核心目标不是追求计算精度，而是强调物理可解释性。用户通过调节参数，可实时观察晶格结构、Hamiltonian 矩阵、能带、态密度（DOS）和费米能级的变化，并阅读对应的物理解释。

5.2 支持的模型

模型	维度	原胞原子数	核心公式
1D 单原子链	1	1	E(k) = ε₀ + 2t cos(ka)
1D 双原子链	1	2	2×2 复厄米矩阵对角化
2D 正方晶格	2	1	E(k) = ε₀ + 2t[cos(k_xa) + cos(k_ya)]
2D 蜂窝（石墨烯）	2	2	Dirac 锥 / 带隙化
3D 简立方	3	1	三维 cos 色散
1D 周期势（NFE）	1	—	平面波展开，BZ 边界开隙

5.3 核心功能

晶格可视化：实时显示实空间原子排布与布里渊区形状。
Hamiltonian 检查器：用户拖动 k 点滑块，实时显示 H(k) 矩阵元、本征值与本征向量，并用柱状图展示轨道权重。
能带与 DOS：沿高对称路径绘制 E_n(k)，同步计算高斯展宽的态密度，自动判定费米能级位置。
电子填充分析：根据用户设定的每原胞电子数 N_e。
有效质量与群速度：点击能带图上任意一点，计算该处的 m* 与 v = (1/ℏ)∂E/∂k。
参数扫描：选择某一参数（如 t、ε、a），在设定范围内自动演化，生成二维色图展示能带变化。
物理故事模式：9 步交互式教学，从孤立原子逐步推导至能带结构。

5.4 高级物理开关

功能	物理机制	程序实现
自旋极化	自旋向上/向下电子独立计算占据数	勾选后 DOS 权重因子调整为 1（原为 2）
Hubbard U	电子关联导致的同轨道双占据能量代价	平均场近似：，SCF 迭代至收敛
磁性	铁磁（FM）/ 反铁磁（AFM）交换劈裂
SOC	自旋–轨道耦合导致能带杂化	，在简并点打开 avoided crossing

5.5 预设实验

核心教学功能层级：

物理结构定义级： 支持一维单原子链、一维双原子链、二维正方晶格以及二维石墨烯蜂窝晶格。用户可自定义晶格常数、原子在位能（Onsite Energy）及轨道跳跃积分（Hopping）。
动态哈密顿量检查器 (Hamiltonian Inspector)： 程序的教学核心。随着用户拖动进度条扫描布里渊区路径，界面同步显示当前点的确切复数矩阵形式，并实时求解输出本征能量及各轨道成分的权重（Orbital Composition）。
高级物理微扰介入： 集成平均场近似下的 Hubbard U 参数（触发铁磁/反铁磁相变与自旋劈裂）及 SOC 微扰项（观察狄拉克点处能隙的开启机制）。
高分辨积分算法： 态密度（DOS）模块摒弃了粗糙的一维线积分，内置高效排序算法，在二维全布里渊区内实施数万个致密 k 点的高斯核密度估计（KDE），精确还原尖锐的范霍夫奇点（Van Hove Singularity）。
故事驱动模式 (Story Mode)： 内置悬浮导览卡片，按步骤控制界面参数，自动演示“孤立原子能级 → 跃迁展宽 → 晶格对称性破缺产生带隙”的连续物理相变过程。

参数	符号	物理意义	对能带的影响
在位能	ε₀, ε_A, ε_B	原子轨道的单电子能量	决定能带整体的能量基准；ε_A ≠ ε_B 时打开带隙
跃迁积分	t, t_ij	相邻原子轨道间的量子力学耦合振幅	\|t\| 越大，能带色散越强，带宽越大
晶格常数	a	相邻原子间距	a 增大 → 轨道重叠减少 → t 减小 → 带宽收缩
Hubbard U	U	同轨道双电子占据的库仑排斥	平均场下产生自旋依赖的在位能，可能导致自旋劈裂与关联带隙
交换劈裂	Δ, J_ex	磁性原子间的交换作用	自旋向上/向下能带整体平移，产生 exchange splitting
SOC 强度	λ	自旋与轨道角动量的耦合	在简并点打开 avoided crossing，可能产生拓扑能隙
周期势幅	V₀	近自由电子模型中的周期势强度	V₀ 越大，BZ 边界处布拉格反射越强，能隙越大
电子填充	N_e	每原胞价电子数	决定费米能级位置；半满为金属，满带为绝缘体
DOS 展宽	σ	高斯展宽宽度	σ 越小，DOS 峰越尖锐；σ 越大，峰越平滑
有效质量	m*	能带极值附近的曲率倒数	曲率越大，m* 越小，载流子迁移率越高
群速度	v	电子波包的传播速度	能带斜率直接对应电子速度；斜率为零处速度为零
直接带隙	E_g^dir	导带底与价带顶位于同一 k 点的能量差	决定光学跃迁效率
间接带隙	E_g^ind	导带底与价带顶位于不同 k 点的能量差	需要声子参与的跃迁