Topology driven and soft phonon mode enabled Na-ion diffusion in quaternary chalcogenides, Na3ZnGaX4 (X = S, and Se)

The compounds Na3ZnGaX4 (X = S, Se) are potential solid electrolyte materials in sodium-based batteries, which have certain advantages over oxide materials and have shown significant ionic conductivity at ambient temperature. In this paper, we bring out atomic-level features of the diffusion process in these new materials using the microscopic techniques of inelastic neutron scattering (INS), quasielastic neutron scattering (QENS), and ab initio molecular dynamics (AIMD) simulations. The insights obtained from these techniques are unique and not available from other macroscopic experiments. Neutron scattering experiments have been performed at temperatures from 100 to 700 K. The simulations have been carried out up to 900 K. We have calculated the phonon spectra and the space–time correlation functions and found good agreement with the results of the neutron scattering experiments. The simulations enable detailed analysis of the atomic-site dependent dynamical information. We observe low-energy phonon modes of ∼6 meV involving the vibrations of certain Na atoms in the lattice. This reveals that the Na at 32g Wyckoff sites (Na2) has sufficiently shallow potential among the two available crystallographic sites. This shallow potential facilitates diffusion. Furthermore, the specific structural topology of the network of interconnected zig-zag chains of the Na2 atomic sites provides the low-barrier energy pathways for diffusion. A small fraction of vacancy defects appears essential for diffusion. We further observe that the Na2 atoms undergo jump-like diffusion to the vacant next or the 2nd next neighbour sites at ∼4 Å. While the QENS experiments reveal the jump-like diffusion and its time scale, detailed analysis of the AIMD simulations shows that the jumps appear mostly along zig-zag chains of the Na2 sites in the tetragonal ab-plane, as well as between the chains along the c-axis. The net diffusion is essentially 3-dimensional, with little anisotropy despite the anisotropy of the tetragonal crystal structure.