Rate Equation for the Transfer of Interstitials across Interfaces between Equilibrated Crystals

技術要約：平衡結晶間の介在物移動に関する速度式

問題提起
本論文は、金属水素化物（Pd、PdPtなど）に代表される、介在型固溶体の充放電に関する実験的観察結果と既存の動力学モデルとの間の乖離に対処している。実験によれば、組成範囲が相変態または混合ギャップの臨界点を含む場合、組成変化の特性時間（$T$）は劇的に（数桁単位で）増大する。電気化学的に制御された水素化動力学を解析するために広く用いられている標準的なバトラー・フォルマー（BV）式は、この観察結果を説明できない。BVの枠組みは、サイトの占有率と基底状態エネルギーが既知の固定値であると仮定し、駆動力のみを化学ポテンシャルの差（$\Delta \mu_S$）として扱う。しかし、ナノスケールの金属水素化物においては、内部拡散が界面移動よりも数桁速いため、バルク相は内部平衡に準即時的に達する。その結果、界面移動は、占有率や基底状態エネルギーが変動する拡張されたエルゴード的な溶質粒子のアンサンブル間で行われることになるが、これは標準的なジャンプ確率モデルでは捉えきれないシナリオである。

手法
著者は、現象論的なジャンプ確率ではなく、統計力学と詳細釣合いの原理に基づいた新しい動力学速度式を導出している。系は、エネルギー障壁（遷移状態 $T$）によって隔てられた2つの介在型固溶体相（AおよびB）の複合体としてモデル化されている。

1. 統計的枠組み: 各相は、溶質粒子の数（$N_S$）とマトリックスサイトの数（$N_M$）が保存された、一定温度（$T$）における正準系として扱われる。バルク拡散は、各相内での内部エルゴード性を維持するのに十分な速さであると仮定される。
2. 遷移状態解析: 「転移構成（$C_p$）」、すなわち粒子がソースサイトにおいて活性化状態にあり、かつターゲットサイトが空である確率を、制約された系の全微視的状態の総和として計算する。
3. 自由エネルギーの分解: 転移構成に達するために必要な自由エネルギー変化（$\delta F$）は、以下の4ステップの仮想プロセスを通じて計算される：
   • ソース相における溶質アンサンブルを圧縮してサイトを解放する（エントロピー的コスト）。
   • ターゲット相における溶質アンサンブルを圧縮してサイトを解放する（エントロピー的コスト）。
   • ソース相から粒子を基準状態へ取り除く（化学ポテンシャルによる仕事）。
   • 粒子を特定の遷移状態へと挿入する（エネルギーコスト）。
4. フラックスの導出: 計算された転移構成の確率を用いて、順方向および逆方向のフラックスの差から、正味のフラックス密度（$J_B$）を導出する。

主要な貢献と結果
主要な貢献は、電荷中性粒子の界面フラックスに関する新しい速度式（式15）の導出である：
$$J_B = \nu_0 \rho_T (1 - x_A)(1 - x_B) \exp\left(-\frac{\epsilon_T}{k_B T}\right) \left( \exp\left(\frac{\mu_S^A}{k_B T}\right) - \exp\left(\frac{\mu_S^B}{k_B T}\right) \right)$$

この結果の主な特徴は以下の通りである：

• 化学ポテンシャルへの明示的な依存性: 化学ポテンシャルの差（$\Delta \mu_S$）に依存するBV式とは異なり、この速度式は2つの相の個別の化学ポテンシャル（$\mu_S^A$ および $\mu_S^B$）に明示的に依存する。
• 空孔分率への依存性: フラックスは空孔分率 $(1-x_A)$ および $(1-x_B)$ によって明示的に変調される。これは、エルゴード的なアンサンブルにおいて空のターゲットサイトを見出す統計的確率を反映している。
• 詳細釣合い: この式は、$\mu_S^A = \mu_S^B$ のときに消失するという詳細釣合いの原理を自然に満たす。
• 動力学的遅延の説明: 混合ギャップの臨界点またはスピノーダル付近では、溶質感受率（$\chi_S = dx/d\mu_S$）が発散する。特性時間 $T$ は、交換フラックス密度と感受率の積に反比例するため（$T^{-1} \propto J_{eq}^\mu / \chi_S$）、$\chi_S$ の発散は $T$ の発散をもたらす。これは、臨界点付近で観察される充電時間の劇的な遅延に対する理論的な説明を与える。
• 非対称性: 導出された法則は、臨界点の両側において特性時間 $T$ の非対称性を予測する。これは実験的に観察されている特徴であるが、対称的なBV式では再現されない。

意義
本論文は、導出された速度式が、バルク相が界面よりも急速に平衡化する系における界面移動の、より物理的に一貫した記述を提供すると主張している。転移サイトを固定されたサイトではなく、エルゴード的な正準系に結合したものとして扱うことで、相互作用する固溶体に固有の、変動する基底状態エネルギーと占有率を考慮に入れている。本研究は、金属水素化物における動力学的異常を説明できないバトラー・フォルマー式の失敗が、平衡組成への明示的な依存性と、空孔の利用可能性の統計的性質の無視に起因することを示唆している。これらの結果は、転送係数や障壁エネルギーへのアドホックな調整を行うことなく、相変態付近におけるエネルギー貯蔵材料の充放電動力学を理解するための理論的基礎を提供するものである。