Charge-Transfer Electronic Structure of NiX$_2$ (X = S, Se)

技術要約：NiX$_2$ (X = S, Se) の電荷移動電子構造

問題提起
ニッケルジカルコゲナイド（NiS$_2$ および NiSe$_2$）の電子的な性質については、電子相関、共有結合、およびカルコゲン二量体（ダイマー）の具体的な役割の相互作用に関して、依然として議論の対象となっている。NiS$_2$ は伝統的に電荷移動（CT）絶縁体（小さなギャップ ~0.2 eV）に分類され、NiSe$_2$ は金属とされるが、固溶体 NiS$_{2-x}$Se$_x$ における金属－絶縁体転移（MIT）を駆動するメカニズムは完全には解明されていない。従来の密度汎関数理論と動的平均場理論を組み合わせた研究（DFT+DMFT）では、硫黄－セレンの反結合性二量体状態がニッケルの 3d 上部ハバードバンド（UHB）の下方に位置するとされていた。しかし、近年の実験的な分光データはこれらのモデルとの不一致を示唆しており、電荷移動エネルギーが著しく小さいか、あるいは負の値である可能性、および電子状態の順序が異なる可能性を示唆している。ニッケル 3d 状態とカルコゲン $p$ 状態の正確な相対的位置は、二重計上補正（$\mu_{dc}$）によって制御されており、実験的な価電子帯、コアレベル、および吸収スペクトルと整合させるためには再評価が必要である。

手法
著者らは、NiS$_2$ および NiSe$_2$ の電子構造を調査するために DFT+DMFT フレームワークを用いた。

• 計算セットアップ: 実験的な結晶構造を用い、WIEN2k パッケージ内の局所密度近似（LDA）を用いて DFT 計算を行った。wien2wannier および wannier90 を用いて、Ni 3d、カルコゲン $p_\pi$、および $p_\sigma$ 軌道を含むタイトバインディング（TB）モデルを構築した。
• 相互作用パラメータ: Ni 3d シェルの相互作用は、ハバード $U = 5.0$ eV およびハウンド結合 $J = 1.0$ eV で扱った。不純物問題は、580 K においてハイブリダイゼーション展開連続時間量子モンテカルロ（CT-QMC）ソルバーを用いて解いた。
• 二重計上補正戦略: 固定されたスキームを用いた従来の研究とは異なり、著者らは二重計上補正（$\mu_{dc}$）を調整可能なパラメータとして扱った。実験的な価電子帯光電子分光（PES）データを再現するように、エネルギー差パラメータ $\Delta_{dp}$（Ni 3d とカルコゲン $p$ 軌道の間のエネルギー差と定義）を系統的に変化させることで $\mu_{dc}$ を最適化した。
• 実験による検証: 理論スペクトルは以下のデータと比較された：
  • SPring-8 および Photon Factory で測定された硬 X 線 PES（HX-PES, $h\nu = 6.5$ keV）および軟 X 線 PES（SX-PES, $h\nu = 1.2$ keV）。
  • Ni 2p コアレベル X 線光電子分光（XPS）および Ni $L_{2,3}$ 端における X 線吸収分光（XAS）。
• モデルのバリエーション: 二量体の反結合状態の役割を孤立させるため、「フルモデル」（カルコゲン $p_\sigma$ 軌道を含む）と、縮小された「Ni-$p_\pi$ モデル」（$p_\sigma$ を除外）を比較した。

主な結果

1. 最適化された電荷移動エネルギー: 計算された価電子帯スペクトルを実験的な PES データに一致させることにより、最適な $\Delta_{dp}$ は NiS$_2$ で 1.82 eV、NiSe$_2$ で 1.58 eV であると決定された。これらの値は、初期の DFT+DMFT 研究で想定されていたよりも小さな電荷移動エネルギーを示しており、近年のクラスターモデル解析と一致している。
2. スペクトル特徴の再現: 洗練されたモデルは、以下を正常に再現した：
   • 軌道依存の断面積の変化を正確に捉えた、両方の HX-PES および SX-PES に対する価電子帯スペクトル。
   • Ni 2p コアレベル PES および NiS$_2$ の XAS スペクトル。これには、最適な $\Delta_{dp}$ 値の近傍でのみ捉えられた Ni 2$p_{3/2}$ 成分における特定のサテライト特徴（St1 および St2）が含まれる。
3. 電子構造と電荷移動:
   • 最適化されたパラメータは、実質的なリガンドから Ni への電荷移動を明らかにしており、その結果、Ni 3d 占有数は約 8.5（NiO の 8.2 よりも大幅に高い）となる。
   • $|t_6e_2\rangle$ と $|t_6e_3p_\pi\rangle$ の構成間の混合が増大することで、Ni の局所モーメントはスクリーニングされた値である 1.38 $\mu_B$ まで減少する。これは低温での実験的観察と一致している。
4. 金属－絶縁体転移（MIT）のメカニズム:
   • NiS$_2$: この系は絶縁体であり、UHB は CT バンドから分離されている。S-S 反結合性二量体状態（$p_\sigma$）はフェルミ準位（$E_F$）より上に位置し、最適な深い Ni $d$ 状態のシナリオにおいては、主に受動的な傍観者として機能する。
   • NiSe$_2$: Se-Se 二量体の分裂はより弱く、反結合性 $p_\sigma$ 状態が NiS$_2$ と比較して低いエネルギーに位置している。特定のパラメータでは Ni-$p_\pi$ 亜系のみではギャップが開いたままとなる可能性があるが、Se $p_\sigma$ 反結合バンドとの弱いハイブリダイゼーションが、系を金属化させるために必要な「押し（nudge）」を提供する。Se $p_\sigma$ バンドが $E_F$ を横切り、ギャップを閉じる。
   • 本研究は、NiS$_{2-x}$Se$_x$ シリーズにおける MIT が、単純なバンド幅制御のみではなく、Ni 上部ハバードバンドと反結合性カルコゲン二量体状態の間の相対的な位置関係によって駆動されることを明らかにしている。

意義および主張
本論文は、二重計上補正の最適化を通じて、実験的な分光データと従来の理論モデルとの間の不一致を解消することを主張している。著者らは、自らのモデルが、低エネルギー電子構造、コアレベル励起、および磁気モーメントの統一的な記述を提供すると断言している。

主な意義は以下の通りである：

• CT エネルギー・スケールの修正: 以前の仮定に対し、UHB と価電子帯上端との間のエネルギー分離を覆し、実験的なコアレベルおよび価電子帯データと一致する、より小さな電荷移動エネルギーを確立したこと。
• MIT メカニズムの解明: NiS$_2$ から金属的な NiSe$_2$ への転移が、カルコゲン二量体の反結合状態の具体的なエネルギー配置と Ni 3d UHB との関係によって支配されていることを示したこと。本研究は、Se $p_\sigma$ 反結合バンドがフェルミ準位と重なることが、二量体分裂の減少と特定の電荷移動エネルギーによって促進され、NiSe$_2$ を金属状態へと導くことを強調している。
• 手法の妥当性確認: $\mu_{dc}$ を実験的な PES に対して調整することが、相関電子材料の電子構造を決定するための堅牢な戦略であることを証明した。これは、価電子帯、コアレベル、および吸収スペクトルの同時再現によって示されている。

著者らは、NiS$_2$ と NiSe$_2$ は非常に類似した Ni-カルコゲン ハイブリダイゼーション強度を持ち、NiS$_2$ は MIT の直近の絶縁体側に位置している一方で、Se 反結合状態の特定の配置が NiSe$_2$ を金属領域へと駆動していると結論付けている。