7. Wannier 轨道与降维投影
本课位置: 第 7/10 课。先完成前一课的推导;使用下方练习检查理解。
1. 目标与假设
Wannier 函数把 Bloch 能带与局域轨道联系起来。先在每个波矢的选定能带子空间内作酉变换,再进行傅里叶变换。不同规范可产生不同轨道形状和跳跃范围,却保持被表示的能带。因此局域化是精心选择的表示,不是唯一原子轨道的额外实验测量。
孤立能带组在各波矢处维数固定;纠缠能带则需从较大窗口选择平滑子空间。外窗口与冻结窗口会影响所得模型。应在独立波矢点比较插值能带和轨道性质,不能只在训练网格上判断。平滑规范帮助局域化,但在强制某些对称条件时,拓扑障碍可能使完全局域表示复杂化。
降维投影可以通过分块代数消去不关心的自由度。将系数分为保留的 A 区和消去的 B 区,先由 B 方程求出 cB,再代回 A 方程,得到依赖能量的有效哈密顿量。若目标远离被消去能级,可以围绕参考能量展开并得到近似能量无关模型。如果靠近消去态的共振,矩阵逆变大,低能近似会失败。
选择子空间后再截断实空间跳跃,是另一层独立近似。应检查跳跃范围收敛,并保证正反向矩阵的厄米关系。复现本征值对许多用途是必要条件,却未必足以保证全部观测量:跃迁矩阵元、Berry 相位和位置算符也需要一致变换。窄窗口内准确的模型,不能未经验证就外推到高激发或显著变化的结构。
2. 逐步推导
逐行阅读等式及其假设。DFT 方程默认采用原子单位,另有说明除外;紧束缚参数保持明确能量和长度单位。匕首表示共轭转置,坐标上的撇号表示独立积分变量。
2.1. 傅里叶插值有明确约定
按选定实空间约定对跳跃乘 Bloch 相位求和重建矩阵。胞内轨道位置可能带额外相位,取决于基定义。跳跃文件、速度和光学矩阵元必须使用同一约定。截断若不均匀删除对称相关项,会破坏对称。应在独立密网格检查费米面附近误差;局域轨道图不证明响应所需矩阵元都变换一致。
3. 计算示例
保留能级零、消去能级二电子伏特,耦合为 0.2 时,低能修正约为 −0.02。完整双能级最低本征值约 −0.0198。两者接近是因为耦合远小于能级差,而不是降维结果普遍不依赖能量。
4. 习题与解释性答案
习题. 解释接近被消去能级时修正为何变大,以及如何处理。
解释性答案. 分母趋近零,微扰消去失去尺度分离。应把该能级纳入保留子空间,或求解依赖能量的问题,不能用拟合常数隐藏共振。
进一步检查. 写出算例各参数的单位和适用范围。改变一个假设,指出推导哪一步必须重做。正确解释应说明受到影响的约束、算符或边界条件,而不只是说结果会变化。
5. 误区与局限
Wannier 规范不唯一。发布模型时应报告窗口、初始投影、规范约束和跳跃截断。
图为原创教学示意,不是数值材料模拟输出。
6. 联系与来源
- Slater–Koster, 轨道几何原始构造 (1954)
- Su–Schrieffer–Heeger, 聚乙炔中的孤子 (1979)
- Marzari et al., 最大局域 Wannier 函数综述 (2012)
解释、推导和假设算例为原创教学综合。所引文献建立相关理论,不是这些教学数值的来源。
7. 关联理论与实践
量子力学 · 密度泛函理论 · ABACUS · DFTB+