Microstructural Insights into Fast Ion Transport in Solid Electrolytes via Multiscale Modeling

技術要約：マルチスケールモデリングによる固体電解質中の高速イオン輸送に関する微細構造的知見

問題提起
高性能な全固体電池の開発は、固体電解質の高速イオン輸送を可能にする微細構造的特徴の理解が不完全であることによって阻害されている。歴史的な取り組みはバルクの超イオン伝導性の最適化に焦点を当ててきたが、実用的な固体電解質は多結晶体であり、緻密な粒界（GB）ネットワークを含んでいる。酸化物電解質における粒界の影響はある程度理解されているものの、硫化物電解質における粒界の役割は、形態の特性評価や粒界の寄与をバルク伝導から分離することの実験的困難さから、依然として不明確である。さらに、多結晶におけるイオン輸送をモデリングすることは計算量的に困難である。これは、原子スケールのホッピング機構と連続体スケールのパーコレーションを橋渡しする必要があり、複雑な粒界相互作用を正確に記述できる原子間ポテンシャルと、従来の第一原理手法では到達できない関連スケールの構造的無秩序を捉えるのに十分な大きさのシミュレーションセルを要求するためである。

手法
著者らは、密度汎関数理論（DFT）、機械学習ポテンシャル（MLP）、分子動力学（MD）、および有限要素（FE）シミュレーションを統合した階層的なマルチスケールモデリングフレームワークを提案している。

1. 機械学習ポテンシャル（MLP）： 精度と計算コストのバランスを取るため、著者らはクローズドループ型の能動学習スキームを用いて、軽量なモーメントテンソルポテンシャル（MTP）を開発した。このプロセスでは、粒界を含む完全な微細構造から局所的な原子環境を抽出し、これらの局所構造に対してDFT計算を行い、近・第一原理精度が得られるまで、複雑性を増しながら反復的にMTPを訓練した。
2. 分子動力学（MD）： 学習されたMTPを用いて、双晶、柱状結晶、およびランダムな粒界ネットワークを持つ多結晶を含むサロゲート粒界モデルの、大規模なMDシミュレーション（最大ナノ秒オーダー）を実施した。これらのシミュレーションにより、300 Kから600 Kの温度範囲におけるLiイオンの無秩序、相関、および拡散キネティクスを捉えた。構造的無秩序を定量化するために、自己拡散係数、拡散障壁、および過剰エントピーを主要な指標とした。
3. 有限要素（FE）シミュレーション： MDから得られた自己拡散係数と粒界幅を、連続体スケールのFEシミュレーションの入力値として使用した。スケール分離を扱うために、2つのアプローチを用いた。
   • 分解能を持つ粒界（Resolved GBs）： ナノグレイン（~25 nm）の場合、フーリエ加速ノードソルバー（FANS）を用いて粒界を体積分解能で解いた。
   • 集約された粒界（Collapsed GBs）： マイクログレイン（~1 µm）の場合、粒界を厚みを集約させた有限要素として表現した。
     この手法により、デバイスに関連する結晶粒径を持つ多結晶系のマクロな拡散係数および導電率の予測が可能となった。

主な貢献と結果
本研究は、アルジロライト型固体電解質 $\text{Li}_6\text{PS}_5\text{X}$（$\text{X} \in \{\text{Cl, Br, I}\}$）に適用され、アニオンの秩序化、ハロゲン化物の置換、および粒界形態の影響を系統的に調査した。

• 粒界のエネルギー論と構造： アニオン秩序化された $\text{Li}_6\text{PS}_5\text{Cl}$ において、粒界形成エネルギーは低く（一般的な粒界で $12\text{--}22\text{ meV \AA}^{-2}$、$\Sigma3$ で $\sim 5\text{ meV \AA}^{-2}$）、これは合成試料において粒界が豊富かつランダムに存在することを示唆している。粒界は構造的無秩序を誘起し、それは過剰エントピーの増加とサブナノメートル規模の空隙によって示されるが、これらの空隙が拡散率に与える影響は最小限であった。
• アニオンの秩序化と置換の影響：
  • アニオンの無秩序化： アニオン無秩序型の $\text{Li}_6\text{PS}_5\text{Cl}$ では、粒界は障壁として作用し、Liイオンの自己拡散率をバルクに対して約60%低下させる。しかし、特徴的なホッピング周波数は依然として高く（$\sim 10^8\text{ Hz}$）、これは実験的に観察される低い粒界抵抗と一致している。
  • ハロゲン化物の置換： 秩序相において、Liイオン拡散障壁とアニオン半径の間に線形スケーリング関係が確立された。アニオン半径を大きくする（$\text{Cl} \to \text{Br} \to \text{I}$）と、拡散障壁は $0.1\text{ \AA}$ の増加につき約 $80\text{ meV}$ 上昇する。
  • 相反する粒界効果： 粒界の影響はバルクの状態に依存する。高速拡散相（アニオン無秩序型の $\text{Li}_6\text{PS}_5\text{Cl}$）では、粒界は拡散を抑制する。一方、低速拡散相（アニオン秩序型の $\text{Li}_6\text{PS}_5\text{I}$）では、粒界は拡散を最大60倍まで促進し、接続チャネルとして機能する。
• 非アレニウス型輸送： シミュレーションは、アニオン秩序型の $\text{Li}_6\text{PS}_5\text{I}$ で観察される非アレニウス型輸送挙動を再現することに成功し、輸送メカニズムがバルク支配（高温域）から粒界支配（低温域）へと移行する転移温度（~250 K）を特定した。この挙動は、この特定の系における粒界の活性化エネルギーがバルクよりも低いことに起因する。
• 結晶粒径の効果： 連続体シミュレーションにより、結晶粒径をナノスケールに微細化することで、非アレニウス挙動が緩和されることが明らかになった。$\text{Li}_6\text{PS}_5\text{Cl}$ の場合、結晶粒の微細化は超イオン伝導性を損なうことなく粒間接触を改善する。一方、$\text{Li}_6\text{PS}_5\text{I}$ の場合、ナノサイズ化は本来の超イオン伝導性を欠く系において輸送を活性化させる戦略となり得る。

意義と主張
本論文は、そのマルチスケールモデリングフレームワークが、原子スケールのメカニズムとマクロなデバイス性能との間の厳密かつ定量的なリンクを提供し、多結晶固体電解質におけるイオン輸送の長年の課題を解決することを主張している。

• メカニズム的知見： 本研究は、粒界が普遍的にイオン輸送を妨げるわけではないことを明確にした。むしろ、粒界の影響（促進か抑制か）は、バルク相の固有の拡散率によって決定される。
• 設計指針： 結果は、微細構造設計のための a priori（先行的）なガイダンスを提供する。$\text{Li}_6\text{PS}_5\text{Cl}$ のような高速拡散硫化物の場合、結晶粒の微細化は、導電性を犠牲にすることなく接触を改善するために有益である。$\text{Li}_6\text{PS}_5\text{I}$ のような低速拡散系については、ナノサイズ化が輸送を活性化させるための戦略となり得る。
• 組成パラメータ： 本研究は、アニオン半径が秩序化されたアルジロライトにおける拡散障壁を支配する主要な組成パラメータであることを特定し、格子極性化のみが支配的な要因であるという概念に異を唱えた。
• 実験的検証： 予測された非アレニウス挙動および転移温度は、$\text{Li}_6\text{PS}_5\text{I}$ における実験的観察と一致しており、粒界がこれらの材料におけるこの現象の起源であることを示唆している。

著者らは、このアプローチが、理想化されたバルクモデルを超えて、現実的な多結晶表現を明示的に考慮することで、高度な固体電解質の設計を可能にすると結論付けている。