Skip to main content
U.S. flag

An official website of the United States government

Official websites use .gov
A .gov website belongs to an official government organization in the United States.

Secure .gov websites use HTTPS
A lock ( ) or https:// means you’ve safely connected to the .gov website. Share sensitive information only on official, secure websites.

Kinetic Lattice Monte Carlo Model for Oxygen Vacancy Diffusion in Praseodymium-Doped Ceria: Applications to Materials Design



Pratik Dholabhai, Shahriar Anwar, James B. Adams, Peter A. Crozier, Renu Sharma


Abstract We present a kinetic lattice Monte Carlo model for investigating the oxygen vacancy diffusion in praseodymium-doped ceria. The current approach uses a database of activation energies for oxygen vacancy migration, calculated using first-principles, for various migration pathways in praseodymium-doped ceria. Also, since our first-principles calculations revealed significant vacancy-vacancy repulsion, we investigate the importance of that effect by conducting simulations with and without a repulsive interaction. The two models show distinct maximum for ionic conductivity as a function of dopant concentration. Initially, as dopant concentrations increase, vacancies increase, so conductivity increases. However, at higher concentrations, vacancies interfere and repel one another, and dopants trap vacancies, creating a "traffic jam" that decreases conductivity, which is consistent with the experimental findings. The activation energy for vacancy migration increased with increasing dopant concentration in agreement with the experimental data. This effect is primarily due to the increasing likelihood of finding two Pr ions near an oxygen vacancy, requiring higher energy to overcome these barriers. There is also a smaller effect of vacancies impeding the movement of other vacancies due to repulsive forces. The current methodology comprising a blend of first-principle calculations and kinetic lattice Monte Carlo techniques provides a very powerful fundamental tool for predicting the optimal dopant concentration in ceria related materials.
Journal of Solid State Chemistry


Doped ceria, vacancy migration, kinetic lattice Monte Carlo mdoeling


Dholabhai, P. , Anwar, S. , Adams, J. , Crozier, P. and Sharma, R. (2011), Kinetic Lattice Monte Carlo Model for Oxygen Vacancy Diffusion in Praseodymium-Doped Ceria: Applications to Materials Design, Journal of Solid State Chemistry, [online], (Accessed June 20, 2024)


If you have any questions about this publication or are having problems accessing it, please contact

Created February 11, 2011, Updated October 12, 2021