Ti/Zr/Hf-Bi/Sb蜂窝单层中的量子反常霍尔效应：第一性原理预测

摘要

二维材料研究不断揭示新拓扑相机遇。量子反常霍尔效应（QAHE）虽在磁掺杂拓扑绝缘体中实现，但仍局限于超低温。利用密度泛函理论，我们预测由过渡金属Ti、Zr、Hf与磷族元素Bi、Sb组成的单层蜂窝结构（MPn，M=Ti、Zr、Hf；Pn=Bi、Sb）具有本征QAHE。研究体系中，平面TiBi和HfBi在适度拉伸应变下分别以15 meV和7 meV带隙呈现QAH绝缘体特性。揭示晶格常数和屈曲距离变化驱动拓扑相变，观察自旋-轨道耦合诱导的带反转，并通过Wannier导出紧束缚模型在锯齿纳米带中确认手性边缘模式。研究拓展现实QAH材料目录，提示低能自旋电子器件可行平台。

研究背景

二维拓扑绝缘体（2D TIs）展现自旋极化无能隙边缘态同时保持体绝缘，是量子自旋霍尔（QSH）效应标志[1]。当时间反演对称性（TRS）被打破——由本征磁性或外部扰动——螺旋边缘通道转化为手性通道，产生量子反常霍尔效应（QAHE）[3]。QAHE无耗散、无磁场电荷输运使其成为超低功耗自旋电子学吸引人平台[4]，但实验实现局限于稀磁拓扑绝缘体[8]。理论工作表明将铁磁性和强自旋-轨道耦合（SOC）引入QSH材料可打开拓扑带隙产生非零陈数[9,10]。尽管多种QSH体系——如IV族（Sn）、V族（Bi、Sb）单层——深入研究，过渡金属基蜂窝结构用于QAHE潜力仍 largely unexplored。本工作填补空白，探索MPn（M=Ti、Zr、Hf；Pn=Bi、Sb）蜂窝结构，建立在III-V蜂窝可宿主QSH相和过渡金属掺杂可诱导磁性的近期预测上[17,25-27]。采用第一性原理计算评估Ti、Zr、Hf在Bi/Sb蜂窝上取代能否稳定本征QAH相。

结果与讨论

在原始蜂窝晶格中用过渡金属（Ti、Zr、Hf）替换半数Bi原子生成MPn单层，可采用低屈曲或平面几何。Fig. 1a显示M-Pn晶格俯视图，Fig. 1b,c为屈曲和侧视图，Fig. 1d为第一布里渊区高对称点。

通过变化晶格常数弛豫原子位置评估结构稳定性和应变效应，针对屈曲和平面两相。Fig. 2a-c绘制TiBi、ZrBi、HfBi不同晶格常数下总能量曲线。低屈曲MPn结构能量恒低于平面对应物。然而QAHE主要出现在应变平面区，仅屈曲HfBi在窄晶格常数窗口支持QAH相（Fig. 2c）。

表1列出所有MPn体系平衡晶格常数、磁矩、带隙、计算陈数。非零整数陈数识别QAH绝缘体：平面TiBi（C=+1，Δ=15 meV）和平面HfBi（C=-1，Δ=7 meV）为最有希望候选。TiBi相变可通过减小屈曲距离驱动（Fig. 3），HfBi由拉伸应变诱导（Fig. 4）。

平衡几何下电子带结构见图5。平面几何中TiBi展现清晰带反转和15 meV带隙，平面HfBi显示7 meV带隙且C=-1。屈曲HfBi为半金属，高陈数（C=-3）。其他构型如屈曲ZrBi、TiSb为平庸绝缘体。

聚焦平面TiBi（a=4.76 Å）剖析SOC和磁性作用。Fig. 6显示非磁计算得金属态，铁磁有序（磁矩≈1.05 μ_B/胞，Ti主导）产生半金属谱。引入SOC打开15 meV带隙，确认带反转由SOC驱动并导致QAHE。

边缘态是拓扑相关键指纹。利用Wannier导出紧束缚哈密顿量，构建127 Å宽锯齿HfBi纳米带（Fig. 7）。带结构显示手性边缘模式穿越费米面奇数次，等于|C|=1，确认平面HfBi QAH本质。

声子计算显示自由站立单层有小虚频，提示实验实现可能需要衬底。此外，2×2超胞中铁磁构型在屈曲情况下比反铁磁序能量更优，平面体系FM和AFM态近简并。

结论

第一性原理研究表明将Ti、Zr、Hf掺入Bi或Sb蜂窝晶格可诱导本征QAHE。屈曲几何虽能量优选，但适度拉伸应变或减小屈曲距离驱动向平面构型转变，宿主非零陈数。平面TiBi和HfBi分别以15 meV和7 meV带隙作为稳健QAH绝缘体浮现。这些材料增补二维QAH候选名册，为低功耗自旋电子器件提供有前景平台。

材料与方法

所有计算在密度泛函理论框架内进行，采用广义梯度近似（GGA）和投影增强波（PAW）势，VASP 5.3实现。平面波截断350 eV，结构弛豫至残余力<5×10⁻³ eV/Å。电子自洽循环收敛至10⁻⁶ eV。z轴20 Å真空间距防止层间相互作用，24×24×1 Gamma中心Monkhorst-Pack网格采样二维布里渊区。用WANNIER90生成最大局域Wannier函数，构建紧束缚模型计算边缘态。拓扑不变量通过Z2Pack包获得，追踪混合Wannier电荷中心得陈数。