Breakdown of Soft Anharmonic Phonons Heralds Fast Ionic Diffusion in Lithium Argyrodite

POSTER

Abstract

A fundamental understanding of the atomic structure and dynamics enabling fast ionic transport in solids is essential for the development of next-generation solid-state electrolytes (SSE). Focusing on the promising SSE candidate Li6PS5Cl with argyrodite structure, we resolve the coupling between fast diffusion of Li+ and vibrational dynamics of the host framework through extensive inelastic and quasielastic neutron scattering measurements, combined with machine-learned molecular dynamics (MLMD) simulations based on first-principles data. Our results establish that host lattice vibrations enable an order-of-magnitude increase in Li+ diffusivity at ambient temperature. Our experiments and simulations both show a clear overlap and interplay of hopping dynamics and vibrational frequencies in the terahertz regime, with a continuous spectral evolution from harmonic phonons to strongly anharmonic overdamped vibrations, and fast Li+ diffusion. We identify the key degrees-of-freedom enabling fast Li diffusion as low-frequency dynamics of PS43- polyanions, which are distinct from the commonly assumed "paddle-wheel" scenario. Bringing together neutron measurements and large-scale MLMD simulations, our results build a "beyond phonons" picture of complex atomic dynamics in SSEs in terms of overdamped spectral functions. These results offer microscopic insights into the mechanism of fast Li+ diffusion in lithium argyrodites and provide guidance for the design of future SSE materials.

*Neutron scattering data collection, MD simulations and analysis by JD was supported by a U.S. National Science Foundation project NSF DMREF under award DMR-2119273. OD acknowledges Duke startup funds for the support of MKG. Sample synthesis by CR and WGZ was supported by the Deutsche Forschungsgemeinschaft under grant number ZE 1010/4-1. The use of Oak Ridge National Laboratory's Spallation Neutron Source was sponsored by the Scientific User Facilities Division, Office of Basic Energy Sciences, U.S. DOE. Theoretical calculations were performed using the National Energy Research Scientific Computing Center, a U.S. DOE Office of Science User Facility supported by the Office of Science of the U.S. DOE under Contract No. DE-AC02-05CH11231.

Presenters

  • Jingxuan Ding

    • Harvard University
    • Duke University

Authors

  • Jingxuan Ding

    • Harvard University
    • Duke University