安全性などの観点から次世代電池として全固体電池が期待され,既に実用化されてきているが,未だに液系に匹敵するイオン伝導度を持つ固体電解質は実用化されていない.焼結された固体電解質における粒界でのイオン伝導抵抗が大きいことが一因だが,粒界抵抗と組成・構造の関係は十分に解明されていない.本研究では,組成・構造が与えられたときに粒界においてどのような原子構造が実現し,イオン伝導抵抗にどう影響するかを分子動力学(MD)シミュレーションを用いて原子スケールから解明することを目的としている. さまざまな組成・構造の候補材料のMDシミュレーションを行うため,量子論に基づく高精度ではあるが計算負荷の大きな MDシミュレーションから計算負荷の低いが十分な近似精度の原子間力場を非常に効率的に構築する手法を開発した.本手法は,従来法とは異なり,高精度MDから動的な物理量の分布関数を参照データとして原子間力場のパラメータを最適化することで,MDシミュレーションを行った際に安定な力場を構築することができる.また,パラメータ探索範囲を自動的に更新する新手法を用いてメタヒューリスティクスを行うことで,従来は試行錯誤を必要としていた作業を効率化することに成功し た.本研究で開発した手法を用いることで,非常に効率的に高精度な原子間力場を構築することができるようになり,新規固体電解質材料であってもイオン伝導挙動のMDシミュレーションを即座に行うことが可能となった. また,多結晶材料のMDシミュレーションから粒界におけるイオン伝導経路を抽出する新手法を開発し,次世代電解質材料候補であるNASICON型LZPの粒界におけるイオン伝導抵抗の原子スケールの機構を明らかにした.
|