Atomistic modeling of bulk and grain boundary diffusion in solid electrolyte Li$_6$PS$_5$Cl using machine-learning interatomic potentials
Yongliang Ou,Yuji Ikeda,Lena Scholz,Sergiy Divinski,Felix Fritzen,Blazej Grabowski
2024-07-31
Abstract:Li$_6$PS$_5$Cl is a promising candidate for the solid electrolyte in all-solid-state Li-ion batteries. In applications, this material is in a polycrystalline state with grain boundaries (GBs) that can affect ionic conductivity. While atomistic modeling provides valuable information on the impact of GBs on Li diffusion, such studies face either high computational cost (\textit{ab initio} methods) or accuracy limitations (classical potentials) as challenges. Here, we develop a quality-level-based active learning scheme for efficient and systematic development of \textit{ab initio}-based machine-learning interatomic potentials, specifically moment tensor potentials (MTPs), for large-scale, long-time, and high-accuracy simulations of complex atomic structures and diffusion mechanisms as encountered in solid electrolytes. Based on this scheme, we obtain MTPs for Li$_6$PS$_5$Cl and investigate two tilt GBs, $\Sigma3(1\bar{1}2)[110]$, $\Sigma3(\bar{1}11)[110]$, and one twist GB, $\Sigma5(001)[001]$. All three GBs exhibit low formation energies of less than \SI{20}{meV/\angstrom\textsuperscript{2}}, indicating their high stability in polycrystalline Li$_6$PS$_5$Cl. Using the MTPs, diffusion coefficients of the anion-ordered and anion-disordered bulk, as well as the three GBs, are obtained from molecular dynamics simulations of atomistic models. At \SI{300}{\kelvin}, the GB diffusion coefficients fall between the ones of the anion-ordered bulk structure (\SI{0.012e-7}{cm^2/s}, corresponding ionic conductivity about \SI{0.2}{mS/cm}) and the anion-disordered bulk structure (\SI{50}{\percent} Cl/S-anion disorder; \SI{2.203e-7}{cm^2/s}, about \SI{29.8}{mS/cm}) of Li$_6$PS$_5$Cl. Experimental data fall between the Arrhenius-extrapolated diffusion coefficients of the investigated atomic structures.
Materials Science
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the influence of grain boundaries (GBs) on lithium - ion diffusion in the solid - state electrolyte Li₆PS₅Cl. Specifically, the research aims to develop high - quality interatomic potentials efficiently and systematically through machine - learning interatomic potentials (MLIPs), especially moment tensor potentials (MTPs), in order to perform high - precision simulations of complex atomic structures and diffusion mechanisms on a large - scale and long - time scale. This helps to overcome the limitations of traditional methods (such as ab - initio methods and classical force fields) in terms of computational cost and accuracy.
### Research Background
- **All - solid - state lithium - ion batteries**: Compared with traditional lithium - ion batteries, all - solid - state lithium - ion batteries have higher safety due to the use of solid rather than flammable liquid electrolytes.
- **Li₆PS₅Cl material**: This material is considered a potential solid - state electrolyte candidate for all - solid - state lithium - ion batteries because of its high ionic conductivity.
- **Influence of grain boundaries**: In practical applications, Li₆PS₅Cl usually exists in a polycrystalline form, and grain boundaries may affect its ionic conductivity. Therefore, it is crucial to understand the influence of grain boundaries on lithium - ion diffusion.
### Research Objectives
- **Develop high - quality MTPs**: Systematically develop MTPs suitable for large - scale, long - time, high - precision simulations by proposing an active learning scheme based on quality levels.
- **Study the influence of grain boundaries on diffusion**: Use the developed MTPs to study the influence of three different types of grain boundaries (Σ3(1̄12)[110], Σ3(̄111)[110] and Σ5(001)[001]) on lithium - ion diffusion.
- **Provide accurate diffusion coefficients**: Obtain the diffusion coefficients of these grain boundaries as well as the ordered and disordered bulk phases through molecular dynamics simulations, thereby enhancing the understanding of experimental data.
### Method Overview
- **Grain boundary construction**: Generate periodic grain boundary structures based on the cubic symmetry of Li₆PS₅Cl and the coincidence - site lattice theory.
- **Self - diffusion coefficient calculation**: Calculate the self - diffusion coefficient through the mean square displacement (MSD) in molecular dynamics simulations.
- **Training and validation of machine - learning potentials**: Gradually improve the quality level of MTP through an active learning scheme, and validate its fitting error in energy and force.
### Main Results
- **Formation energy of grain boundaries**: The formation energies of the three grain boundaries are all lower than 20 meV/Ų, indicating that they have high stability in polycrystalline Li₆PS₅Cl.
- **Diffusion coefficient**: At 300 K, the diffusion coefficient of the grain boundaries is between that of the ordered bulk phase (1.2×10⁻⁹ cm²/s) and that of the disordered bulk phase (2.2×10⁻⁷ cm²/s).
- **Consistency with experimental data**: The experimental data fall within the Arrhenius extrapolated diffusion coefficient range predicted theoretically, supporting the quantitative computer prediction.
Through these studies, the paper provides an important theoretical basis for understanding the influence of grain boundaries on lithium - ion diffusion in Li₆PS₅Cl, which is helpful for optimizing the design and performance of all - solid - state lithium - ion batteries.