Vacancy Diffusion Across FeCrAl Alloy Composition Space for Accident-Tolerant Fuel Cladding

技術要約：事故耐性燃料被覆材としてのFeCrAl合金組成空間における空孔拡散

問題提起
鉄-クロム-アルミニウム（FeCrAl）合金は、ジルカロイと比較して優れた高温酸化耐性を有することから、軽水炉用事故耐性燃料（ATF）被覆材の有力な候補となっている。しかし、これら被覆材の炉内寿命は、放射線誘起微細構造進化、特に点欠陥輸送（具体的には空孔拡散）によって支配されている。保護的な酸化層の形成についてはよく解明されているが、照射下におけるFeCrAlの基礎的な欠陥物理は依然として理解が進んでいない。既存の理論的扱いは、理想化された、あるいは組成的に簡略化された環境に依存することが多く、濃縮固溶体特有の複雑な局所的化学揺らぎを考慮できていない。これらの局所的な環境は、空孔の結合、移動経路、およびクラスター形成やトラッピングといった集団的挙動を著しく変化させる。したがって、熱および照射荷重の複合条件下における長期的な材料進化、相安定性、およびクリープを予測するために、構成元素（Fe、Cr、Al）の全体的な空孔輸送への個別の寄与を分離することが極めて重要である。

手法
著者らは、一定の5 at.% Al濃度を持つ、系統的なFe/Cr組成範囲（FeリッチからCrリッチまで）にわたる種別解像型キネティックモンテカルロ（KMC）フレームワークを用いている。手法は以下のステップで進行する：

1. ポテンシャルおよびデータ生成： 原子間相互作用は、密度汎関数理論（DFT）データベースに適合させた埋め込み原子法（EAM）ポテンシャルを用いて評価される。空孔移動データは「Hop-Decorate」ワークフローを用いて生成される。ここでは、BCC Fe参照構造を、目標とする全組成に従ってFe、Cr、Al原子でランダムにデコレートする。最小エネルギー経路は、Fe$_{80}$Cr$_{15}$Al$_5$、Fe$_{47.5}$Cr$_{47.5}$Al$_5$、およびFe$_{15}$Cr$_{80}$Al$_5$の全組成に対して、Climbing-Image Nudged Elastic Band（CI-NEB）法を用いて計算される。
2. サロゲートモデリング： KMC実行中のオンザフライでのNEB計算の計算コストを回避するため、生成された障壁データベースに基づき、線形サロゲートモデルを学習させる。各ホップの局所的な化学環境は、遷移状態（サドルポイント）および移動原子の周囲の固定された動径殻内のFe、Cr、Alの原子数をカウントする記述子を用いてエンコードされる。各移動種（Fe、Cr、Al）に対して個別の線形回帰モデルを適合させ、種固有の障壁依存性を捉える。
3. KMCシミュレーション： 学習されたサロゲートモデルにより、500 Kから1500 Kの温度範囲で滞留時間KMCシミュレーションを実行する。シミュレーションでは、$10^5$ステップにわたって単一空孔の軌跡を追跡し、平均二乗変位（MSD）を計算して、拡散係数（$D$）、活性化エネルギー（$E_a$）、および前指数因子（$D_0$）を抽出する。

主要な結果

• 移動障壁と局所環境： 移動障壁は局所的な化学環境に対して非常に敏感であることを本研究は確認している。純粋なFeの障壁は約0.63 eVであるが、Crは約0.95 eV、Alは約0.69 eVである。合金中では、特にCrとAlにおいて障壁の分布が著しく広がり、局所的な化学的非対称性が移動経路を調節していることを示している。
• アルミニウムの役割： 低い濃度（5 at.%）であるにもかかわらず、Alは移動障壁に対して不釣り合いなほど大きな影響を及ぼす。3つのサロゲートモデルすべてにおいて、第1最近接（1NN）サドルシェルにAl原子が存在すると、エネルギー障壁が著しく上昇する（Feに対して+0.048 eV、Crに対して+0.073 eV、Alに対して+0.151 eV）。
• 拡散係数の組成依存性：
  • Feリッチから等原子比： Fe$_{80}$Cr$_{15}$Al$_5$からFe$_{32.5}$Cr$_{62.5}$Al_5}の組成範囲において、空孔拡散係数は全組成に対して比較的鈍感である。活性化エネルギーは0.70–0.84 eVの範囲に収束しており、これはこれらの領域ではFeの移動が輸送を支配していることを示唆している。
  • Crリッチ領域： Crリッチ組成（Fe$_{15}$Cr$_{80}$Al$_5$）は、1.135 eVという著しく高い活性化エネルギーを示す。その結果、そのマクロな拡散係数は、低温域（例：570 Kにおいて$10^{-17}$ m$^2$ s$^{-1}$）において、Feリッチ合金（~$10^{-13}$ m$^2$ s$^{-1}$）よりも数桁低い。
• パーコレーションと「スラグ・ディフュージョン（遅い拡散）」： 結果は、低障壁経路（主にFeのホップによって支配される）の連結ネットワークが、Cr含有量の増加に伴って崩壊するというパーコレーション閾点の存在を示唆している。Crリッチ領域では、空孔移動は高障壁のCrホップによって支配され、「スラグ・ディフュージョン」をもたらす。著者らは、化学的に粗いランドスケープにおける空孔の軌跡が、トラッピングポケットからの脱出を伴うことを観察しており、拡散領域に達する前に非線形な平均二乗変位挙動を示すことを確認している。

意義および主張
本論文は、FeCrAl合金のマクロな拡散係数は、純粋な端成分の障壁の単純な組成加重平均によって予測されるのではなく、局所的なエネルギーランドスケープの不均一性と、低障壁移動経路の特定の連結性によって決定されると主張している。

• 設計への示唆： 本研究は、空孔の移動度が抑制される明確な組成領域（Crリッチ）を特定している。この抑制は、Crリッチの$\alpha'$相（放射線下で形成されることが知られている）において空孔の移動度を低下させ、空孔・格子間原子の再結合速度を減少させ、欠陥蓄積のキネティクスを変化させる可能性がある。
• 事故耐性： 温度が1200 K以上に上昇する冷却材喪失事故（LOCA）条件下において、拡散係数の組成依存性は顕著になる。著者らは、高い活性化エネルギーと抑制された拡散性を持つCrリッチ組成が、既知の酸化耐性を補完する形で、事故耐性被覆材のための有利な空孔キネティクスを提供する可能性があると示唆している。
• 手法的な貢献： 本研究は、NEBデータベースによってパラメータ化された種別解像型線形サロゲートモデルを用いた、低次化KMCフレームワークが、複雑な組成空間にわたる欠陥輸送を効果的にマッピングできることを実証している。このアプローチは、経験的な試行錯誤を超え、原子論的プロセスと燃料性能コードで使用されるメゾスケールモデルとの間のメカニズム的な繋がりを提供するものである。

著者が指摘する限界事項
著者らは、本モデルが格子を静的かつ凍結したものとして扱っており、局所的な緩和や時間の経過に伴う放射線誘起組成進化（偏析）を無視していることを明示的に述べている。また、線形サロゲートモデルは方向性や対称性の破れ、あるいは磁気的な寄与によるエネルギーランドスケープの変化を捉えることはできない。さらに、本研究は空孔輸送のみに焦点を当てており、格子間原子クラスターの輸送や、転位ループの放射線誘起変態（これらも欠陥シンク密度に影響を与える可能性がある）については扱っていない。