Reduced-Order Modelling of Defect Transport using Surrogate Kinetics: Application to U$_x$Pu$_ {1-x}$N

技術要約：サロゲート・キネティクスを用いた $U_xPu_{1-x}N$ における欠陥輸送の低次モデル化

問題提起
$U_xPu_{1-x}N$ のような化学的に無秩序な材料における欠陥輸送のモデリングは、大きな課題を伴う。化学的に均一な格子では、欠陥の移動経路は有限の遷移速度の集合として列挙できる。しかし、無秩序系では局所的な化学環境が連続的に変化するため、移動障壁は閉じた集合ではなく、分布として生成される。この化学的複雑さにより、密度汎関数理論（DFT）や分子動力学（MD）による全環境の列挙は計算量的に極めて困難となる。さらに、標準的な連続時間マルコフ連鎖（CTMC）モデルでは、局所的な無秩序によって生成される実質的に無限の種類の遷移速度を扱うことができない。この問題は、核分裂によって $U_xPu_{1-x}N$ 固溶体へと進化する可能性がある、あるいは混合燃料として製造されるウランモノナイト（UN）のような高度な核燃料にとって極めて重要である。燃料の性能（膨潤、核分裂生成物ガスの放出、熱伝導率の低下）を予測するには、組成の変化が欠陥のエネルギー論および輸送にどのように影響するかを理解する必要があるが、膨大な構成空間にわたる直接的な原子論的シミュレーションは不可能である。

手法
著者らは、原子スケールのメカニズムを計算効率の高い関数形式に埋め込むための、低次モデル化ワークフローを提案している。本手法は、主に4つの段階で構成される：

1. データ生成 (Hop-Decorate): LAMMPSおよびASE上に構築されたHop-Decorate (HopDec) ワークフローを用い、約3,000個の移動障壁（$E_b$）およびエネルギー差（$\Delta E$）のデータベースを生成した。これは、参照用UN構造における典型的な最近接（1NN）移動経路を特定し、それを $U_{0.5}Pu_{0.5}N$、$U_{0.1}Pu_{0.9}N$、および $U_{0.9}Pu_{0.1}N$ の全組成において1,000回「再装飾（re-decorating）」することで達成された。原子間エネルギーは、DFTデータによって検証されたFinnis–Sinokai EAMポテンシャルとBuckinghamポテンシャルの組み合わせを用いて計算された。各構成のエネルギー評価には、Climbing-Image Nudged Elastic Band (CI-NEB) 法が用いられた。
2. サロゲートモデルの構築: 局所的な化学環境は、欠陥の初期、終状態、および鞍点（saddle-point）配置における特定の最近接（NN）シェル内のウラン原子の数によって定義された。これらの単純な記述子に基づき、$\Delta E$ と $E_b$ を予測するための線形サロゲート関数をフィッティングした：
   • $\Delta E$ は、極小値シェルの間のUカウントの変化の線形関数としてモデル化された。
   • $E_b$ は、極小値および鞍点シェルのUカウント、および種特異的なベースライン補正（$\beta$）の線形関数としてモデル化された。
   • アクチニド副格子上のU媒介ホップとPu媒介ホップでは、強い種依存性があるため、個別のフィッティングが必要であった。
3. キネティックモンテカルロ（KMC）シミュレーション: サロゲートモデルをカスタム格子ベースKMC（laKMC）シミュレータに統合した。$10 \times 10 \times 10$ nm$^3$ のスーパーセル（64,000原子）を用い、100,000ステップのシミュレーションを実行した。様々な温度および組成における長時間の拡散係数（$D$）を算出するために、平均二乗変位（MSD）を抽出した。統計的平均化は、各条件につき50個の軌跡に対して行われた。
4. トラッピング解析: マクロな輸送の原子スケールにおける起源を解明するため、各組成について、無秩序平均滞留時間（$\langle \tau \rangle_c$）および速度加重平均エネルギー変化（$\langle \Delta E \rangle_c$）を計算することで、速度論的・熱力学的空間を解析した。

主な結果
本研究では、3種類の欠陥種（アクチニド空孔（$V_A$）、窒化物空孔（$V_N$）、およびアクチニド・窒化物二空孔（$V_{A-N}$））を解析した。

• 移動エネルギー論: サロゲートモデルは低い平方根平均二乗誤差（RMSE）を達成した。これは、$U_xPu_{1-x}N$ における移動エネルギー論が、複雑な長距離の集団効果ではなく、主に短距離の化学的同一性によって支配されていることを示している。

  • $V_A$ (アクチニド空孔): 移動は種によって強く制御される。障壁分布は、移動する原子がUかPuかによって分岐する。局所環境による影響はわずかな広がりを与えるのみであり、主要な組成効果は、UホップとPuホップの統計的な重み付けである。
  • $V_N$ (窒化物空孔): 移動は真の環境依存的な変調を示す。障壁分布は組成とともに連続的に進化し、中間組成において、異なる局所的なアクチニド構成に対応する多峰性の特徴を示す。特に、Puに富むマトリックス中の希薄なUは、局所的に有利な低障壁経路を作り出す。
  • $V_{A-N}$ (二空孔): 結合した欠陥は、いずれかの副格子を介して進行する孤立した空孔の挙動に類似した挙動を示す。このペアの存在は、障壁の高さに系統的なシフトをもたらすが、局所的な配位に対する感度を根本的に変えることはない。

• 拡散の傾向:

  • 非線形性: 拡散は、UNとPuNの端成分間の単純な線形補間（ベガード則）に従わない。
  • $V_A$ の挙挙動: 少量のPuの添加により、拡散は急激かつ非線形に増加する。これは、高速なPu媒介ホップのパーコレーションによって駆動される。Pu濃度が増加するにつれ、システムは遅いUホップによって制限される状態から、高速なPuホップのネットワークへと移行する。
  • $V_N$ の挙動: 低濃度Puにおいて、拡散は「停滞拡散（sluggish diffusion）」を示す。ここでは、Puが速度論的なトラップとして機能し、いくつかの低障壁イベントが存在するにもかかわらず、正味の輸送を減少させる。高濃度Puでは、低障壁経路がより普及し接続されるため、拡散が増加する。
  • $V_{A-N}$ の挙動: 二空孔の拡散は、両方の副格子の化学的感度の間の競争を反映しており、その結果、孤立した空孔とは異なる複雑で非線形な組成依存性を示す。

• トラッピング解析: 速度論的・熱力学的マップ（$\langle \tau \rangle_c$ 対 $\langle \Delta E \rangle_c$）により、以下が明らかになった：

  • $V_A$ については、Pu含有量の増加は、エネルギー的なバイアスを大きく変えることなく滞留時間を減少させる。これは、高速経路のパーコレーションを裏付けている。
  • $V_N$ については、少量のPuは、滞留時間は安定している一方でエネルギー的なバイアス（$\langle \Delta E \rangle_c$）を増大させる。これは、相対的なエネルギーを単に増大させるのではなく、エネルギーを変調する化学的トラップの形成を示している。

意義および主張
本論文は、物理的な解釈性を欠く「ブラックボックス」的な機械学習モデルに頼ることなく、化学的に無秩序な材料における欠陥輸送をモデリングするための、一般的かつ計算効率の高い枠組みを提供すると主張している。

• 低次表現: 著者らは、無秩序なアクチニド窒化物の複雑な移動エネルギー論が、局所的な配位数の線形関数を用いた予想外に単純な低次表現によって記述できることを示した。
• メカニズムの洞察: 本研究は、輸送を支配する2つの主要なメカニズムを特定した：種制御の速度論（アクチニド空孔に支配的）および環境依存のトラッピング（窒化物空孔に支配的）。また、マクロな拡散係数は、単なる平均障壁値ではなく、障壁分布、経路の接続性、および滞留時間の間の競争から創発することを強調している。
• スケーラビリティ: 原子論的データをコンパクトなサロゲート関数に凝縮することで、このアプローチは、原子スケールのメカニズムをより高次のシミュレーション（例：クラスター動力学や連続体燃料性能コード）へ統合することを可能にする。これにより、燃焼に伴う組成の変化を効率的にモデリングすることが可能となる。
• 適用可能性: 本フレームワークは $U_xPu_{1-x}N$ に対して実証されているが、ブルートフォースによる局所環境の列挙が困難な、高エントロピー合金やドープされた半導体などの他の組成無秩序材料にも適用可能である。

著者らは、二空孔を最近接ペアに限定していること（解離・再結合を無視している）、および電子効果（例：Puの5f軌道の局在化）を明示的に扱わない経験的ポテンシャルを使用していることなどの限界についても述べている。しかし、現在の手法は、原子論的なメカニズムと、予測的な燃料性能モデリングのための扱いやすい輸送方程式との間のギャップを埋めることに成功していると主張している。