六方晶窒化ホウ素における三角形空孔構造：サイズ依存的安定性、磁性、エッジ対欠損効果

概要

密度汎関数理論（DFT）を用いて、単層六方晶窒化ホウ素（h-BN）における三角形空孔の原子構造、電子構造、磁気的振る舞いを調べた。計算の結果、三角形空孔の緩和幾何構造はそのサイズに強く依存し、すべての安定構造がN終端ジグザグエッジを持つことが明らかになった。空孔エッジからホウ素-窒素（BN）対が除去されると、結果として得られる構造の安定性と磁気モーメントは欠損対の位置に応じて変化し；面領域から対が除去された場合に最も安定な構造となり、最小磁気モーメントを生じる。

背景

六方晶窒化ホウ素（h-BN）は、グラフェンに類似した二次元絶縁体だが、BとN原子が交互に構成されている。合成過程で、h-BNは必然的に点欠陥、特に空孔を含み、これらは電子的・機械的特性を変化させ、結果としてh-BNベースナノデバイスの性能に影響を与える。格子が2原子種を含むため、空孔エッジはN終端またはB終端となり得る；前者がエネルギー的に優先され、しばしばジグザグモチーフを採用する。先行理論研究では、三角形空孔の安定性と磁性がエッジ終端と空孔サイズに依存することが示されている。実験では、電子線照射下での三角形空孔形成が確認されており、空孔は三角形を維持しながら成長する。最近の研究では、h-BN原子が空孔エッジで個別ではなく束となって放出されることを観察し、空孔成長における協調的メカニズムを示唆している。本稿では、その観察を拡張し、三角形空孔構造、サイズ依存的緩和、BN対除去の影響について包括的なDFT解析を提供する。

計算手法

すべての計算は、VASP（ウィーンab initioシミュレーションパッケージ）を用い、平面波カットオフ400 eV、プロジェクタ拡張波（PAW）疑似ポテンシャルで実施した。交換相関エネルギーは、グリムDFT-D2ファンデルワールス補正を加えたパードュー・バーク・エルンツェルホフの一般化勾配近似（GGA）で処理し、分散相互作用を捕捉した。原子幾何構造は残留力が0.01 eV/Å以下になるまで緩和した；ブリルアンゾーンサンプリングは、（9×9）と（15×15）スーパーセルでΓ点のみを用いた。計算格子定数2.56 Åは実験値と一致する。

結果と考察

h-BN中の三角形空孔

N終端エッジに焦点を当て、いくつかの三角形空孔サイズを調査した。N終端三角形には2つの異なる緩和幾何構造が現れた：頂点原子がほぼ変化しない対称的なN-symm構造、および各頂点でN-N準同種結合が形成される歪んだNN-bond構造である。最小空孔（V_1B）ではN-symm構造のみが見つかり、N-N距離2.66 Åは、原始格子（2.48 Å）に比べてわずかに拡大している。空孔サイズがV_3B+1NとV_6B+3Nに増加すると、N-symmとNN-bondの両構造が局所的に安定であり、先行研究と一致する。

より大きな空孔（V_10B+6N以上）では、NN-bond構造が唯一の安定幾何構造となる。この領域では、空孔周辺のB-N結合は本質的に変化せず、頂点にN-N結合が持続する。表1（ここでは未掲載）は、空孔サイズ増加に伴いN-symmとNN-bond構造のエネルギー差が縮小し、頂点再構築の全格子への影響が減少することを示す。

これらの空孔の磁気モーメントは、エッジ原子終端と相関する。N-symm構造では、各ダングリングN原子が1 μ_B寄与し、エッジN原子数に等しいモーメントを生む。対照的に、NN-bond構造では、N-N結合がダングリング結合性を消失させるか保持するかにより、磁性が減少または増強される。例えば、NN-bond状態のV_3B+1N、V_6B+3N、V_10B+6Nはそれぞれ0、3、6 μ_Bのモーメントを持つ。

空孔エッジでのBN対除去

実験観察では、空孔成長がエッジからBN対の放出により進行することが多い。これを模倣するため、大きなV_21B+15N空孔（15×15スーパーセル）のエッジの様々なサイトからBN対を除去した。6つの異なる欠損位置を検討した：2つの頂点、2つの頂点付近の面、2つの中央面。緩和後、最適化幾何構造は3クラスに分類される：頂点欠損、面欠損、中央欠損、それぞれ特徴的な再構築を示す。

頂点欠損構造は、わずかに歪んだBN環が形成される以外は最小限の変化しか示さない。対照的に、面欠損構造は空孔に隣接してN二量体と五員環を形成し、このN二量体が安定性と磁性に大きく影響する。相対エネルギーは、三角形中心に近い位置で対が欠損した構成がエネルギー的に優位であることを示す（表2）。磁気モーメントは、頂点欠損の12 μ_Bから、特定の面欠損ケースの9 μ_Bまで変化し、不対電子の再分配を反映する。

スピン密度マップ（図3）は、不対スピンがエッジN原子、特に歪んだBN環とN二量体に局在することを確認する。状態密度解析（図4）は、原始h-BNのバンドギャップ内に欠陥状態が存在し、フェルミ準位付近にスピン非対称特徴を示す。3つの代表的欠損対構成のバンドギャップはそれぞれ0.35、0.24、0.36 eVである。

結論

本第一原理研究は、h-BNにおける三角形空孔の安定性と磁性が、空孔サイズ、エッジ終端、BN対除去の正確な位置によって支配されることを実証した。大きな空孔は頂点でNN-bond構造を優先的に採用し、これが磁気モーメントを決定する。大きな三角形ホールのエッジからBN対が除去されると、最も安定な結果はN-N結合を形成する面欠損構造であり、その磁気特性と局所状態密度は欠損位置に敏感である。これらの知見は、電子デバイスおよびスピントロニクスデバイスにおける欠陥工学的h-BNナノ構造の調整のための微視的基盤を提供する。