本研究では,全固体電池の性能を左右する電極・電解質界面におけるミクロな反応機構を明らかにするため,分子動力学シミュレーション(MD)の中で電荷移動を伴う原子構造変化やイオンダイナミクスの再現が可能な原子間力場を効率的に構築する手法およびプログラムを開発する.開発した手法を用いてLi金属とNASICON型固体電解質の界面のシミュレーションを可能とする力場を構築し,電極・電解質 界面のMDシミュレーションを行い,ナノスケールの観点から界面構造変化やイオン伝導機構をより深く理解できる方法論を確立することを目的とする. 初年度から前年度にかけて,電荷可変型クーロン力場と環境依存型のMEAM力場のプログラムに加え,効率的パラメータ最適化プログラムを汎用化して,MDプログラムで計算可能なあらゆる物理量を教師データとして最適化することを可能とするプログラムを開発し,公開した. Li金属負極からのデンドライト形成・成長を抑制する可能性のあるコーティング材として期待されているLiFを対象に,Li/LiF界面における反応機構解析を行うための反応力場の構築に取り組んだ.Li単体,LiF単体,F2分子における生成エネルギー差,各構造における動径・角度分布関数をよく再現する力場を構築し,Li/LiFにおける原子構造を第一原理計算と比較して精度検証を行った.電荷可変型とMEAN型力場を用いて,Li/LiFおよび空孔形成エネルギーなど所望の物性全てを十分な精度で再現するのは困難であることが示された. 精度に関して妥協したLi/LiF力場を用い,LiF粒界を解析した結果,粒界近傍ではF空孔が形成されやすいこと,Liイオンは自由体積の大きい粒界に沿って速く拡散することが示され,LiFの粒界を介して高いLiイオン伝導を確保したLi金属コーティング材となる可能性があることを示した.
|