Li diffusion in sulfide solid electrolyte

Introducition

Ion-conducting solid electrolytes have been known for a long time, but have undergone remarkable development in the past decade or so. In particular, sulfide-based solid electrolytes have dramatically improved their ionic conductivity and are expected to be used in all-solid-state lithium-ion batteries.

Li10GeP2S12 (LGPS) belongs to the category that shows the highest Li-ion conductivity among the existing sulfide solid electrolytes, and its unique crystal structure contributes to the high performance.

In this section, we compute the diffusion coefficient of Li-ion in LGPS by molecular dynamics calculation using Matlantis.

The sample script is uploaded on github:

image

Calculation models and methods

In this section, we will optimize the crystal structure of LGPS and calculate the diffusion coefficient of Li.

The crystal structure of LGPS can be obtained from the original paper [1] or Materials Project [2]. Based on the obtained structure, we will analyze the behavior of lithium ion by molecular dynamics after confirming the reproducibility by structural optimization.

We adopt NVT ensemble using Langevin dynamics simulation.

image

Results and Discussion

image
image

It is known that in LGPS, diffusion in the c-axis (=z direction) is most prominent and limited in the perpendicular directions. The MSD plot clearly shows such feature. The diffusion coefficients can be computed from the MSD and shown in the Arrhenius plot.

The activation energy is found to be in good agreement with the DFT values as well as the experimental value reported in the literature.[3-4]

This type of MD simulations are often performed at very high temperatures when a conventional DFT is used. However, it can be performed at much lower temperatures using Matlantis because of the low computational cost, and the results can be more readily compared with the experimental data.

Computational Details

image

References

[1] N. Kamaya, et. al., Nature Mater 10, 682–686 (2011). https://www.nature.com/articles/nmat3066 [2] https://materialsproject.org/ [3] Mo et al. Chem. Mater. (2012) 24, 15-17 https://pubs.acs.org/doi/10.1021/cm203303y [4] Y. Kato, et. al. Nat. Energy 1, 16030. https://www.nature.com/articles/nenergy201630

tag

what's new

NEW

SiO₂ Dry Etching Simulation Using LightPFP

semiconductor

Analysis of Surface Reaction Mechanism between ALD Precursor and Substrate

ALDSemiconductor

Catalyst screening for Ammonia Synthesis Catalyst

Catalysts

Analysis of CO2 adsorption dynamics of MOF using NEB method

AdsorbentMOF

Calculation of Hildebrand solubility parameters

Hildebrand solubility parameterSP value