跳转至
课程笔记 · 07

7. Wannier 轨道与降维投影

本课位置: 第 7/10 课。先完成前一课的推导;使用下方练习检查理解。

1. 目标与假设

Wannier 函数把 Bloch 能带与局域轨道联系起来。先在每个波矢的选定能带子空间内作酉变换,再进行傅里叶变换。不同规范可产生不同轨道形状和跳跃范围,却保持被表示的能带。因此局域化是精心选择的表示,不是唯一原子轨道的额外实验测量。

孤立能带组在各波矢处维数固定;纠缠能带则需从较大窗口选择平滑子空间。外窗口与冻结窗口会影响所得模型。应在独立波矢点比较插值能带和轨道性质,不能只在训练网格上判断。平滑规范帮助局域化,但在强制某些对称条件时,拓扑障碍可能使完全局域表示复杂化。

降维投影可以通过分块代数消去不关心的自由度。将系数分为保留的 A 区和消去的 B 区,先由 B 方程求出 cB,再代回 A 方程,得到依赖能量的有效哈密顿量。若目标远离被消去能级,可以围绕参考能量展开并得到近似能量无关模型。如果靠近消去态的共振,矩阵逆变大,低能近似会失败。

选择子空间后再截断实空间跳跃,是另一层独立近似。应检查跳跃范围收敛,并保证正反向矩阵的厄米关系。复现本征值对许多用途是必要条件,却未必足以保证全部观测量:跃迁矩阵元、Berry 相位和位置算符也需要一致变换。窄窗口内准确的模型,不能未经验证就外推到高激发或显著变化的结构。

2. 逐步推导

逐行阅读等式及其假设。DFT 方程默认采用原子单位,另有说明除外;紧束缚参数保持明确能量和长度单位。匕首表示共轭转置,坐标上的撇号表示独立积分变量。

\[ \begin{aligned} |w_{n\mathbf R}\rangle&=N_k^{-1/2}\sum_{\mathbf k}e^{-i\mathbf k\cdot\mathbf R}\sum_m|\psi_{m\mathbf k}\rangle U_{mn}(\mathbf k),\\ H_{mn}(\mathbf R)&=\langle w_{m0}|H|w_{n\mathbf R}\rangle,\\ H_{\mathrm{eff}}(E)&=H_{AA}+H_{AB}(E-H_{BB})^{-1}H_{BA},\\ c_B&=(E-H_{BB})^{-1}H_{BA}c_A. \end{aligned} \]

2.1. 傅里叶插值有明确约定

按选定实空间约定对跳跃乘 Bloch 相位求和重建矩阵。胞内轨道位置可能带额外相位,取决于基定义。跳跃文件、速度和光学矩阵元必须使用同一约定。截断若不均匀删除对称相关项,会破坏对称。应在独立密网格检查费米面附近误差;局域轨道图不证明响应所需矩阵元都变换一致。

3. 计算示例

保留能级零、消去能级二电子伏特,耦合为 0.2 时,低能修正约为 −0.02。完整双能级最低本征值约 −0.0198。两者接近是因为耦合远小于能级差,而不是降维结果普遍不依赖能量。

4. 习题与解释性答案

习题. 解释接近被消去能级时修正为何变大,以及如何处理。

解释性答案. 分母趋近零,微扰消去失去尺度分离。应把该能级纳入保留子空间,或求解依赖能量的问题,不能用拟合常数隐藏共振。

进一步检查. 写出算例各参数的单位和适用范围。改变一个假设,指出推导哪一步必须重做。正确解释应说明受到影响的约束、算符或边界条件,而不只是说结果会变化。

5. 误区与局限

Wannier 规范不唯一。发布模型时应报告窗口、初始投影、规范约束和跳跃截断。

最近邻链及余弦色散

图为原创教学示意,不是数值材料模拟输出。

6. 联系与来源

相关:密度泛函基础 · 多体蒙特卡罗方法

DFTB+

解释、推导和假设算例为原创教学综合。所引文献建立相关理论,不是这些教学数值的来源。

7. 关联理论与实践

量子力学 · 密度泛函理论 · ABACUS · DFTB+


上一课 · 课程目录 · 下一课