Kinetic Monte Carlo Simulations for Solid State Ionics: Case Studies with the MOCASSIN Program

Abstract

Kinetic Monte Carlo simulations are a useful tool to predict and analyze the ionic conductivity in crystalline materials. We present here the basic functionalities and capabilities of our recently published Monte Carlo software for solid state ionics called MOCASSIN, exemplified by simulations of several model systems and real materials. We address the simulation of tracer correlation factors for various structures, the correlation in systems with complex migration mechanisms like interstitialcy or vehicle transport, and the impact of defect interactions on ionic conductivity. Simulations of real materials include a review of oxygen vacancy migration in doped ceria, oxygen interstitial migration in La-rich melilites, and proton conduction in acceptor doped fully hydrated barium zirconate. The results reveal the impact of defect interactions on the ionic conductivity and the importance of the defect distribution. Combinations of these effects can lead to unexpected transport behavior in solid state ionic materials, especially for multiple mobile species. Kinetic Monte Carlo simulations are therefore useful to interpret experimental data which shows unexpected behavior regarding the dependence on temperature and composition.