Dimer-mediated cation diffusion in the stoichiometric ionic conductor Li3N.

Non-equilibrium molecular dynamics has been used to model cation diffusion in stoichiometric Li3N over the temperature range 50 < T/K < 800. The resulting diffusion coefficients are in excellent agreement with the available experimental data. We present a detailed atomistic account of the diffusion process. Contrary to the conclusions drawn from previous studies, our calculations show that it is unnecessary to invoke the presence of a small concentration of intrinsic defects in order to initiate diffusion. The structure can be considered to consist of alternating layers of composition Li2N and Li. As the temperature increases an increasing number of cations leave the Li2N layers and migrate either to the interlayer space or to the Li layer. Those that move into the interlayer space form Li2 dimers with cations in the Li2N layers and those that move into the neighboring layer form dimers with cations therein. The two types of dimer are aligned parallel and perpendicular to [001], respectively and have lifetimes of ∼3 ps. The vacancies so created facilitate rapid diffusion in the Li2N layers and the interlayer cation motion results in slower diffusion perpendicular to the layers.

Diffusion in Li2O studied by non-equilibrium molecular dynamics for 873 < T/K < 1603.

The use of non-equilibrium molecular dynamics facilitates the calculation of the cation diffusion constant of Li2O at temperatures too low to be accessible by other methods. Excellent agreement with experimental diffusion coefficients has been obtained over the temperature range 873 < T/K < 1603. Diffusion below 1200 K was shown to be dominated by a concerted nearest-neighbour hopping process, whereas in the high-temperature superionic region an additional mechanism involving a six-coordinate interstitial cation site in the anti-fluorite structure becomes increasingly dominant. Our model thus accounts for the transition from the superionic regime to the non-superionic regime.