2. 紧束缚模型基础与一维单原子链

紧束缚（Tight-Binding, TB）模型的核心思想是：晶体中的电子在大部分时间里被束缚在原子附近，其行为可以用原子轨道（如 s、p、d 轨道）的线性组合来近似描述。当原子聚集成晶体时，相邻原子的轨道发生交叠，电子得以在格点之间“跳跃”，从而形成能带。这一图像既保留了原子轨道的物理直觉，又能以极低的计算代价给出能带结构，是理解固体电子结构的第一块基石。

一、哈密顿量的构造

考虑一维单原子链，晶格常数为 a，每个原胞含一个原子，原子轨道记为 |n⟩（n 为格点编号）。在二次量子化语言下，紧束缚哈密顿量写为：

H = Σ_n ε₀ |n⟩⟨n| − t Σ_⟨n,m⟩ ( |n⟩⟨m| + |m⟩⟨n| )

其中 ε₀ 为格点轨道能（on-site energy），t 为最近邻跳跃积分（hopping integral），⟨n,m⟩ 表示只对最近邻求和。负号是约定俗成的写法，使 t > 0 时低能态对应成键组合。若考虑次近邻，可再加入 −t₂ 项，但单原子链的最小模型只需 ε₀ 与 t 两个参数。

需要强调符号约定：本文取跳跃项为 −t，t > 0。若读者采用 +t 约定，则色散关系整体变号，但物理结论（带宽、带底带顶位置）不变。全文与图示统一采用 −t 约定。

二、布洛赫定理与能带解析解

一维周期性势场下，布洛赫定理保证本征态可写为 ψ_k(x) = e^{ikx} u_k(x)，其中 u_k 具有晶格周期性。在紧束缚基下，取布洛赫和：

|k⟩ = (1/√N) Σ_n e^{ikna} |n⟩

代入哈密顿量并利用 ⟨n|H|m⟩ 只在 n=m 与最近邻时非零，可得色散关系：

E(k) = ε₀ − 2t cos(ka),   k ∈ [−π/a, π/a]

由此可直接读出几个关键结论：带宽 W = 4t（从 ε₀−2t 到 ε₀+2t）；带底在 k=0（成键态），带顶在 k=±π/a（反键态）；态密度在带边发散（一维范霍夫奇点），在带中心取极小值。这些解析结果将成为后续数值离散化方法的校验基准。

三、数值离散化方法

解析解只对理想无限链成立。实际仿真中，我们通常把链截断为 N 个格点，并施加边界条件。最直接的做法是构造 N×N 的实空间哈密顿矩阵：对角元全为 ε₀，次对角元全为 −t，其余为零。用 Python 的 NumPy 可几行写出：

import numpy as np

N, eps0, t, a = 200, 0.0, 1.0, 1.0
H = np.zeros((N, N))
for i in range(N):
    H[i, i] = eps0
    if i < N-1:
        H[i, i+1] = -t
        H[i+1, i] = -t

# 开边界：直接对角化
evals = np.linalg.eigvalsh(H)
print(evals.min(), evals.max())   # 应接近 -2t 与 +2t

开边界（有限链）下，本征值不再连续，而是 N 个离散能级，且带边附近出现明显的边界态修正。若改用周期性边界条件，只需补上 H[0, N−1] = H[N−1, 0] = −t，此时本征值应严格落在 [ε₀−2t, ε₀+2t] 内，且与解析色散在 k = 2πm/(Na) 处采样一致。这是检验代码正确性的第一道关卡。

另一种更贴近能带图像的离散化是k 空间采样：直接在布里渊区内取 M 个 k 点，对每个 k 计算 E(k) = ε₀ − 2t cos(ka)，得到色散曲线与态密度。两种方法互为验证：实空间对角化的本征值排序后，应与 k 空间采样值逐一对应（在 PBC 下）。

四、从解析到数值的对照检查

下表给出一个可直接用于自检的对照表（示例参数：ε₀=0，t=1，a=1，N=200）：

需要说明，上表数值为按模型参数推算的示例结果，用于说明检查逻辑，并非某次实际运行的实测记录。读者在自己的环境中运行后，数值应在同一量级附近浮动。

五、常见陷阱与约定

• 符号约定不一致：−t 与 +t 两种写法都常见，混用会导致带底带顶颠倒。建议在代码注释中显式写明约定。
• 边界条件遗漏：忘记补 PBC 的角元，会把周期链算成开链，带边出现虚假能级。
• 单位与量纲：t 与 ε₀ 应使用同一能量单位（如 eV），a 与 k 的单位需匹配（k 用 1/a 为单位时可直接取 [−π, π]）。
• N 太小：N < 50 时边界修正显著，带宽与带边位置偏差可达百分之几，不宜用于定量对照。

掌握了一维单原子链的哈密顿量构造、布洛赫定理与两种离散化路径后，后续的双原子链、缺陷链以及更复杂的晶格模型都可以在此基础上扩展：只需改变基组数目、跳跃矩阵元与边界条件，而“解析给出基准、数值给出实现、两者互相对照”的工作方式将贯穿整个仿真流程。