Abstract
Glass-ceramics are promising materials for use as solid electrolytes in all-solid-state batteries because of their high Li+ conductivity. This study investigated an interfacial model for glass-ceramics consisting of 75Li2S-25P2S5 glass and beta-Li3PS4 with seven different beta-Li3PS4 planes using molecular dynamics simulations. The equilibrated structure at 300 K and crystal growth at 500 K were analyzed in terms of interfacial enthalpy and crystallization rates. A novel method for identifying PS43- units as either glassy or crystalline was established based on the rotational motion of PS43- units. The rotational motion of PS43- units in the interfacial model was quantitatively indexed using sulfur trajectories within a 100 ps time window. The crystallization rate calculated using crystal growth simulations was found to depend on the crystalline plane forming the interface. The predicted crystal shape, derived from the interfacial enthalpy and growth rate, was analyzed using Wulff's theorem. The Li+ diffusivity in the crystalline, glassy, and interfacial regions of the interfacial model was evaluated through the time integration of velocity--velocity autocorrelation functions. Notably, the Li+ diffusivity in the crystalline regions varied depending on the interfacial model and differed from that of bulk beta-Li3PS4. This variation was attributed to differences in the degree of PS43- rotational motion in the crystalline regions of the interfacial model. The analytical methodologies developed in this study and the insights into the effects of interfaces on PS43- rotational motion, make valuable contributions to the understanding and design of glass-ceramics with interfaces.
Supplementary materials
Title
Supporting Information
Description
• Initial cell dimensions and number of atoms used for crystal growth simulations.
• ql distributions (l = 1, 2, 3, and 4) calculated from bulk β-Li3PS4 and 75Li2S–25P2S5
glass based on 100 ps MD data at 500 K.
• Snapshot of the equilibrated structure at 300 K showing no crystal growth for the model
with the (010) plane.
• Initial density and lattice constants of the cell equilibrated at 300 K.
• Four unique alignments (i–iv) considered in the interfacial models with seven crystalline planes.
• Snapshots of all models with saturated crystal growth at 500 K.
• Relationship between C(t) and time progression in crystal growth simulations for all
models.
• Table of crystallization kinetics, interfacial energy, and equilibrium crystallinity for the interfacial model with seven crystalline planes.
• Li+ autocorrelation function and integrated results for evaluating the Li+ diffusion coefficient.
Actions
Title
A movie of crystal growth
Description
A movie of crystal growth for all interfacial models at 500 K
Actions