Probing Anharmonic and Heterogeneous Carrier Dynamics Across Sublattice Melting in a Minimal Model Superionic Conductor

Despite decades of research, the microscopic origin of sublattice melting and fast ion transport in superionic conductors remains elusive. Here, we introduce a chemically neutral minimal binary model consisting of a rigid host lattice stabilized by short-range steric repulsion and a soft carrier sublattice interacting via long-range Wigner-type forces. This contrast naturally produces distinct melting temperatures and an intermediate sublattice-melting phase in which carriers become fluidlike while the host remains crystalline. Molecular-dynamics simulations identify three dynamical regimes-crystalline, sublattice-melt, and fully molten-marked by sharp changes in diffusivity, structural correlations, and dynamic heterogeneity. Near sublattice melting, carrier motion is strongly anharmonic and spatially heterogeneous, beyond mean-field hopping descriptions. By tuning the density, we demonstrate that sublattice melting can be continuously controlled, establishing a direct link between lattice softness, anharmonicity, and collective ion transport. This work provides a unified microscopic foundation for designing mechanically robust, high-performance superionic conductors operable near ambient conditions.

💡 Research Summary

The paper addresses the long‑standing question of why certain solid electrolytes exhibit a sublattice‑melting transition, in which the mobile ion sublattice becomes fluid‑like while the host framework remains crystalline, leading to superionic conductivity. To isolate the essential physics, the authors construct a chemically neutral minimal binary model composed of (i) a rigid host lattice stabilized solely by short‑range steric repulsion and (ii) a soft carrier sublattice that interacts via long‑range Wigner‑type (∝1/r) forces. The non‑additive pair potentials give distinct effective length scales for host‑host, host‑carrier, and carrier‑carrier interactions, creating two separate energy scales and consequently two different melting temperatures.

Molecular dynamics simulations are performed in two‑ and three‑dimensional periodic boxes over a wide temperature range. The authors compute mean‑squared displacements (MSD) for host and carrier particles separately, extract self‑diffusion coefficients D_H and D_C, and analyze structural correlations through radial distribution functions g(r). Three dynamical regimes emerge clearly when plotting D versus inverse temperature (1/T): (I) a low‑temperature crystalline regime where both sublattices are ordered and diffusion follows an Arrhenius law with similar activation energies; (II) an intermediate sublattice‑melting regime where the carrier sublattice loses long‑range order (g_CC becomes featureless) while the host retains sharp peaks, and D_C increases by two orders of magnitude relative to D_H; (III) a high‑temperature fully molten regime where both sublattices are disordered and diffusion is rapid for both species.

Crucially, the transition from regime I to II is continuous for the carriers: D_C rises smoothly without a discontinuity, indicating that transport is not governed by isolated hopping events but by collective, correlated motion. To probe this, the authors evaluate anharmonicity through non‑Gaussian parameters and four‑point dynamic susceptibility χ₄. Near the onset of sublattice melting, both NG and χ₄ exhibit pronounced peaks, revealing strong spatial heterogeneity and string‑like cooperative rearrangements of carriers. These signatures are analogous to those observed in glass‑forming liquids, yet they occur while the host lattice remains crystalline, highlighting a unique hybrid state.

The study further explores density (or carrier‑to‑host ratio) as a tuning knob. Lower densities increase the average carrier spacing, weaken carrier‑carrier interactions, and lower the sublattice‑melting temperature T_m. Higher densities strengthen the soft sublattice, pushing T_m upward and even suppressing melting altogether. This systematic control demonstrates a direct link between lattice softness, anharmonic vibrational fluctuations, and the emergence of fast ion transport.

By comparing the minimal model’s predictions with experimental trends in known superionic conductors such as AgI, PbF₂, and Li₁₀GeP₂S₁₂, the authors argue that the essential ingredients for superionicity are (a) a stiff host that preserves structural integrity, (b) a soft, highly polarizable carrier sublattice that can undergo selective melting, and (c) the resulting anharmonic, dynamically heterogeneous environment that enables collective ion migration. The model therefore provides a unified microscopic framework that bridges mean‑field defect theories and recent observations of collective ion dynamics.

In conclusion, the paper demonstrates that (1) non‑additive interactions naturally generate a sublattice‑melting transition, (2) the melting of the carrier sublattice is accompanied by strong anharmonicity and dynamic heterogeneity, and (3) these phenomena can be continuously tuned via density or interaction parameters. This insight offers a practical design principle for engineering robust, high‑performance solid electrolytes that operate near ambient conditions, by deliberately engineering a “soft” ion sublattice within a mechanically stable host framework.