Quantifying Symmetry Breaking as a Design Variable for Giant Altermagnetic Spin Splitting

問題提起

アルターマグネット（交代磁性体）は、反転対称性や並進対称性ではなく、回転対称性によって駆動される、正味の磁化がゼロで運動量依存のスピン分裂（SSE）が1 eVを超える、第三のコロニアル磁性体のクラスである。磁気点群解析は、ある材料がアルターマグネティズムを示す「可能性」があるかどうかを厳密に決定する（二値分類）が、スピン分裂の「大きさ」を予測するメカニズムは提供しない。デバイス性能を左右するSSEは、化学的に多様な材料間で約10 meVから約1 eVまで連続的に変化する。現在の探索ワークフローは、候補となるすべての材料に対してフルスピン偏極密度汎関数理論（DFT）に依存しており、評価コストが発見コストと同等になるという循環依存を生んでいる。結晶構造と組成をSSEの大きさに直接マッピングするための、連続的かつ解釈可能なフレームワークが欠如しているため、この分野は既知の構造タイプの列挙に留まっており、新しい設計仮説を生成するには至っていない。

手法

著者らは、二値的な対称性分類を、連続的で物理的に解釈可能な記述子に置き換える、データ駆動型の逆設計フレームワークを提案している。

1. データセットの生成:

   • ソース: MatterGen拡散モデルを用いて生成された3,851個の二元構造（323個の親化合物を含む）。
   • ラベル付け: 占有ウィンドウ $[E_F - 2 \text{ eV}, E_F]$ 内における最大スピンチャネルエネルギー分離として定義される、DFT計算によるSSEで構造にラベルを付与した。
   • スクリーニング: 強固なアルターマグネティック挙動を持つ候補をフィルタリングした：正味の磁化がほぼゼロ（$<0.005 \mu_B$）、$\Gamma$点でのスピン分裂が無視できる程度、および有意な非相対論的SSE（$\ge 1 \text{ meV}$）。
   • 拡張: 化学組成を変えずに記述子の分布を広げるため、各親化合物に対して等方性および一軸歪みサンプリング（$\pm 2.5\%, \pm 5\%$）による拡張を行った。

2. 記述子の開発:
   著者らは、結晶座標と元素の同一性のみから導出される、52個のDFTフリーなモチーフレベルの記述子を導入した。以下の3つの主要な記述子が支配的であると特定された。

   • モチーフ対称性の破れ指数 (MSBI): 並進（$tT$）および反転（$PT$）の関係によってスピン縮退が強制される、二つの反平行な磁気モチーフからの偏差を定量化する連続的なスカラー値。これは、以下の2つの無次元の偏差の積として計算される：
     • $D_d$: 二つのサブ格子間のリガンドベクトル集合間の、結合長正規化RMSD（ハンガリアン・アルゴリズムによりマッチング）。
     • $D_I$: 一方のリガンド集合を点反転させた後の偏差。
     • $MSBI = D_d \times D_I$。この乗法形式により、両方の縮退保護対称性が破れた場合にのみMSBIが大きくなることが保証される。
   • モチーフ充填率 (MPF): リガンドケージが単位格子の体積（または平面モチーフの場合は面積）に占める割合であり、超交換相互作用のスケール（$J \sim t^2/U$）を制御する。
   • p/d 電子比: 金属のd電子に対するリガンドのp電子の比であり、Zaanen–Sawatzky–Allenの枠組みにおける共有結合性とp–d混成の強さのプロキシ（代理指標）として機能する。

3. 機械学習と最適化:

   • サロゲートモデル: 52個の記述子を用い、3,851個のDFTラベル付き構造を用いてXGBoostモデルを学習させた。このモデルは、グループ交差検証において $R^2 = 0.70$、平均絶対誤差（MAE）123.1 meVを達成した。
   • 解釈性: 特徴量の重要度をランク付けし、SSEを駆動する物理的メカニズムを特定するために、SHAP（SHapley Additive exPlanations）分析を用いた。
   • 逆設計: ベイジアン最適化（BO）を、SSEを最大化する候補を記述子空間内で探索するためにXGBoostサロゲートと結合させた。
   • 検証: 上位の候補は、記述子空間における余弦類似性を用いて実在の構造プロトタイプにマッピングされ、スピン偏極DFTを用いて独立して検証された。

主な結果

1. 3つの記述子の支配性:
   SHAP分析により、MSBI、MPF、およびp/d電子比がモデルの全説明力の37.4%を占めていることが明らかになった。

   • MSBI は支配的な駆動因子（寄与度23.2%）であり、「準必要条件のゲート」として機能する。大きなSSE（$>0.4$ eV）は、本質的にMSBI $> 0.4$ の領域に限定されている。
   • MPF はSSEと正の相関を示し、密な充填がホッピング積分と超交換相互作用を強化する。
   • p/d 比 は逆の関係を示す。低い比率（$p/d < 1$）は、強いスピン非対称なp–d混成と高いSSEを促進し、$p/d \ge 1$ は分裂を抑制する。

2. 設計ルールとメカニズム:

   • 対称性の破れ: MSBIは、対称性の破れの連続的な度合いを定量化する。CuO（MSBI = 1.018, SSE = 1.634 eV）とFeSi（MSBI = 0.0003, SSE = 0.003 eV）のケーススタディは、指数が分裂の大きさと相関していることを示している。
   • 共有結合性の制御: 同じ $P6_3/mmc$ ホスト格子を共有するVOとCrSbの制御比較では、組成を変更してp/d比を1.0から0.6に下げたことで、SSEが7倍（0.165 eVから1.194 eVへ）増加し、共有結合性のメカニズムを孤立させた。

3. 新しい候補の発見:
   BOワークフローは、独立したDFTによって検証されたいくつかの高SSE候補を特定した：

   • $\alpha$-NiS (NiAs型): 予測SSE 0.893 eV、DFT検証済みSSE 0.823 eV。この材料は、アルターマグネットとして独立して特定されている。
   • CoS (八面体型): 予測SSE 0.874 eV、DFT検証済みSSE 1.103 eV。g波アルターマグネットとして独立して報告されている。
   • FeAs (八面体型): 予測SSE 0.867 eV、DFT検証済みSSE 1.089 eV。アルターマグネティックな同定はこれまでになされていない。
   • FeS (平面四角形型): 予測SSE 1.061 eV、DFT検証済みSSE 1.297 eV。アルターマグネティックな同定はこれまでになされていない。
   • 本フレームワークは、Cu-Sのような候補が非磁性状態へと緩和することをも正しく特定しており、これはサロゲートが対称性や共有結合性に焦点を当てていることを強調している。

意義と主張

本論文は、アルターマグネットの発見を、離散的な対称性分類問題から、連続的で解釈可能な構造・物性設計問題へと再構成することを主張している。

• 対称性のための連続変数: 主要な貢献はMSBIであり、これは「アルターマグネティズムが許容されるか？」という二値の問いを、「それはどの程度強く実現されるのか？」という連続的な変数へと変換する。これにより、勾配ベースおよびサロゲートベースの最適化手法の使用が可能になる。
• 実行可能な設計空間: MSBI、MPF、およびp/d比を支配的な軸として特定することで、著者らは実験的なハンドル（格子歪み、圧力、化学置換）を用いてSSEを調整する手段を提供している。
• 効率性: 本フレームワークは、DFTフリーの記述子を用いて候補をスクリーニングすることで発見コストを削減し、最も有望な設計に対してのみ高価なDFT計算を予約することを可能にする。
• 汎用性: 著者らは、この記述子のロジック（連続的な対称性の破れの指標と、幾何学的および共有結合性の軸の組み合わせ）が、アルターマグネティズムを超えた他の対称性駆動の電子特性にも拡張可能であることを示唆している。

本研究は、熱力学的安定性や合成経路を解決すると主張しているのではない。むしろ、特定の配位モチーフとコロニアルな秩序が安定化できるのであれば、本フレームワークは結果として生じる巨大なスピン分裂を予測できると主張している。$\alpha$-NiSおよびCoSの検証、ならびにFeAsおよびFeSの予測は、本フレームワークが既知の物理を回収し、新しいモチーフレベルの設計仮説を提案する能力を持っていることを示している。