Ti/Zr/Hf-Bi/Sbハニカム単層における量子異常ホール効果：第一原理予測

要旨

二次元材料研究の進展は新規トポロジカル相の機会を継続的に明らかにしています。量子異常ホール効果（QAHE）は最近磁気ドープトポロジカル絶縁体で実証されましたが、極低温に限定されています。密度汎関数理論を用いて、遷移金属Ti、Zr、Hfとnictogen元素Bi、Sb（MPn、M=Ti, Zr, Hf; Pn=Bi, Sb）からなる単層ハニカム構造が本質的QAHEを宿すことを予測しました。検討系の中で、平面TiBiとHfBiが控えめな引張歪み下でQAH絶縁体として現れ、それぞれ15 meVと7 meVのバンドギャップを持ちます。格子定数とバックリング距離の変化がトポロジカル相転移を駆動し、スピン軌道結合によるバンド反転を観察し、ジグザグナノリボンでのカイラルエッジモードをWannier導出タイトバインディングモデルで確認しました。これらの知見は現実的なQAH材料のカタログを拡充し、低エネルギースピントロニクスデバイスの実現可能なプラットフォームを示唆します。

背景

二次元トポロジカル絶縁体（2D TIs）はスピン偏極した無ギャップエッジ状態を示しながら絶縁バルクを維持し、量子スピンホール（QSH）効果の特徴です[1]。時間反転対称性（TRS）が固有磁性や外部摂動で破れると、らせんエッジチャネルがカイラル化し、量子異常ホール効果（QAHE）が生じます[3]。QAHEの散逸なし、磁場フリー電荷輸送は超低消費電力スピントロニクスに魅力的ですが[4]、実験的実現は希薄磁性トポロジカル絶縁体に限定されています[8]。理論研究では、QSH材料に強磁性と強いスピン軌道結合（SOC）を導入するとトポロジカルバンドギャップが開き、非ゼロチェルン数が生成されると示されています[9, 10]。IV族（Sn）やV族（Bi, Sb）単層など多くのQSH系が広く研究されてきた一方、遷移金属ベースハニカムのQAHEポテンシャルは未探索領域でした。ここでは、Bi/SbハニカムへのTi、Zr、Hf置換が本質的QAH相を安定化できるかを第一原理計算で探索し、このギャップを埋めます。

結果と考察

プリスタインハニカム格子のBi原子半分を遷移金属（Ti, Zr, Hf）で置換すると、MPn単層が得られ、低バックルまたは平面幾何学を採り得ます。図1aにM-Pn格子の上面図、図1b, cにバックルと平面構造の側面図、図1dに第1ブリルアンゾーン（BZ）と高対称点を示します。

格子定数を変化させ原子位置を緩和することで、バックルと平面両相の構造安定性と歪み効果を評価しました。TiBi, ZrBi, HfBiの全エネルギー曲線を図2a-cにプロットしました。低バックルMPn構造が平面対応物より一貫して低エネルギーです。しかし、QAHEは主に歪み平面領域で出現し、バックルHfBiでのみ狭い格子定数窓でQAH相を支持します（図2c参照）。

表1に全MPn系の平衡格子定数、磁気モーメント、バンドギャップ、計算チェルン数をまとめます。非ゼロ整数チェルン数がQAH絶縁体を特定します：平面TiBi（C=+1, Δ=15 meV）と平面HfBi（C=-1, Δ=7 meV）が最も有望です。TiBiではバックリング距離減少（図3）、HfBiでは引張歪み（図4）で相転移が駆動されます。

平衡幾何学での電子バンド構造を図5に示します。平面幾何学でTiBiは明確なバンド反転と15 meVギャップを示し、平面HfBiは7 meVギャップでC=-1を示します。一方バックルHfBiは高チェルン数（C=-3）の半金属です。バックルZrBiやTiSbなど他の構成は自明な絶縁体のままです。

SOCと磁性の役割を解明するため、平面TiBi（a=4.76 Å）に焦点を当てます。図6に示すように、非磁性計算では金属状態ですが、強磁性秩序（セルあたり磁気モーメント約1.05 μB、Ti主導）で半金属スペクトルになります。SOC導入で15 meVギャップが開き、バンド反転がSOC駆動でQAHEを担うことを確認しました。

エッジ状態はトポロジカル相の本質的指紋です。Wannier導出タイトバインディングハミルトニアンを用いて、127 Å幅ジグザグHfBiナノリボンを構築しました（図7）。バンド構造はフェルミレベルを奇数回（|C|=1に等しい）横切るカイラルエッジモードを示し、平面HfBiのQAH性を確認しました。

フォノン計算でフリースタンディング単層に小さな虚数周波数が見られ、実験実現には基板が必要かもしれません。また、2×2超胞内で強磁性配置がバックルケースで反強磁性秩序よりエネルギー的に有利で、平面系ではFMとAFM状態が準縮退です。

結論

本第一原理研究は、BiまたはSbハニカム格子へのTi, Zr, Hf置換が本質的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ガンマ中心Monkhorst-Packグリッドで2Dブリルアンゾーンサンプリング。最大局在化Wannier関数をWANNIER90で生成し、タイトバインディングモデル構築とエッジ状態計算を可能にしました。トポロジカル不変量はZ2Packパッケージで、ハイブリッドWannier電荷中心追跡によりチェルン数を取得しました。