Decoupling Many-Body Interactions in CeO2 (111) Oxygen Vacancy Structure: Insights from Machine-Learning and Cluster Expansion

技術要約：CeO₂ (111) 酸素空孔構造における多体相互作用のデカップリング

問題提起
酸素空孔（$V_O$）は、セリア（$CeO_2$）の酸化還元特性および触媒特性において極めて重要である。しかし、$CeO_2(111)$ 表面におけるこれらの空孔の分布と安定性を正確に予測することは、電子配置の膨大な数に起因して計算量的に極めて困難である。空孔濃度が増加するにつれ、関連する $Ce^{3+}$ イオン（ポラロン）の数も増加し、構成可能な状態数は指数関数的に増大する。これらの構造の安定性は、$Ce^{3+}$-$Ce^{3+}$、$V_O$-$Ce^{3+}$、$V_O$-$V_O$ の反発および引力を含む複雑な多体相互作用に依存している。従来の計算手法では、これらの構成を効果的にサンプリングすることが難しく、これまでの研究の多くは孤立した空孔や低濃度に焦点を当てており、サブサーフェス（下層）を超えた深い酸素層における空孔の潜在的な安定性を見落としてきた。

手法
著者らは、多体相互作用をデカップリング（分離）するために、第一原理計算、機械学習（ML）、およびクラスター展開（CE）を組み合わせたハイブリッド理論フレームワークを開発した。

1. 第一原理計算: スピン分極を考慮した DFT+U 計算（PBE 関数、$U=5.0$ eV）を $CeO_2(111)$ スラブに対して行った。研究対象は、上位3つの酸素層と上位2つのカチオン層とした。$Ce^{3+}$ の同定は、局所磁気モーメント（$\sim 1 \mu_B$）および $Ce$-$O$ 結合長に基づいて行った。
2. クラスター展開モデル: いかなる構成の形成エネルギーも、有効クラスター相互作用（$J_\alpha$）の総和として表す CE モデルを構築した。初期のプールには、8 Å のカットオフ内にある最大4点（$V_O$ または $Ce^{3+}$ を表す）を含む 3,144 個のクラスターが含まれた。
3. 機械学習 (LASSO): LASSO（Least Absolute Shrinkage and Selection Operator）回帰を用いて、最適なクラスターのサブセットを選択し、過学習を防ぎつつ、最も重要な相互作用項を特定した。訓練セットは、1000 K でのメトロポリス・モンテカルロ（MC）サンプリングを用いた能動学習ループによって反復的に精緻化された。
4. モンテカルロ・サンプリング: CE モデルが収束した後（予測誤差の閾値を $V_O$ あたり <0.05 eV と設定）、このモデルを用いて $8 \times 8$ スーパーセルにおける広範な MC シミュレーション（$10^8$ ステップ以上）を実行した。これにより、様々な濃度（$\Theta = 1/16, 1/8, 3/16, 1/4$）および温度（300 K から 1100 K）における空孔分布を探索した。

主な貢献と結果

• 相互作用のデカップリング: 本研究は、初期の 3,144 個のクラスターから 15 個の主要なクラスターを特定することに成功し、多体相互作用を効果的にデカップリングした。モデルは、訓練セットで最大絶対誤差（MAE）0.07 eV/$V_O$、交差検証で 0.08 eV/$V_O$ を達成した。
• モチーフの特定: 安定性を支配する特定の局所モチーフが明らかになった。
  • 孤立空孔: サブサーフェスの空孔は、表面の空孔よりも $\sim 0.2$ eV 低いエネルギーで安定する。
  • 相互作用: 異なる層に存在する隣接空孔間には強い反発が存在する。逆に、特定のパターン、例えばサブサーフェスの空孔と第3層の空孔が対になったもの（モチーフ $VO_2VO_3\_5N$）は、強い引力（$-0.16$ eV）を示す。
  • ポラロンの局在化: $Ce^{3+}$ イオンは一般に表面サイトを好むが、次近接（NNN）空孔に対して強い引力を示す。
• 濃度依存的な安定性:
  • 低濃度（$\Theta \le 1/8$）では、空孔は主にサブサーフェス層を占有しており、これは先行文献と一致している。
  • 高濃度（$\Theta \ge 3/16$）では、空孔は第3酸素層にも入り込み始める。$\Theta = 1/4$ における最も安定な構成（A2233 と命名）は、サブサーフェスに2つ、第3層に2つの空孔を持つ。この構造は、以前に提案されていた全サブサーフェス型の $(2 \times 2)$ 秩序構造よりも 0.09 eV 安定である。
• 安定化メカニズム: 高濃度における第3層空孔の安定性は、広範な幾何学的緩和によって駆動されている。結合切断エネルギー（$E_b$）は第3層の空孔にペナルティを与えるが、A2233 構成における緩和エネルギー（$E_r$）（表面の $Ce^{4+}$ イオンが反対側の結合の不在により、より自由に上方へ緩和できる）がそのペナルティを上回る。対照的に、すべての空孔が第3層にある構成（A3333）は、過大な結合切断エネルギーのため不安定である。
• 温度効果: MC シミュレーションによれば、温度と濃度が上昇するにつれて、表面層に対する第3層における空孔の存在確率が高まることが示された。ただし、高温時には表面空孔も時折出現する。

意義
本論文は、多体相互作用をデカップリングすることで、金属酸化物における複雑な空孔構造を解明するための新しい理論的枠組みを導入したと主張している。主な意義は、還元された $CeO_2(111)$ における酸素空孔は、高濃度において表面やサブサーフェスに留まるのではなく、幾何学的緩和を最大化するために第3酸素層へと移動するという発見にある。この知見は、第3層を考慮していなかった従来のセリア表面化学の解釈に再考を迫る可能性がある。機械学習とクラスター展開を組み合わせた本アプローチは、明確な格子構造を持つ他の金属酸化物を研究するための堅牢な手法として提示されており、酸化還元および触媒メカニズムをより正確に理解するための道筋を示すものである。