Optimization of ionic conductivity in doped ceria using density functional theory and kinetic lattice Monte Carlo

Fuel cells, particularly solid oxide fuel cells (SOFC), are important for the future of greener and more efficient energy sources. Although SOFCs have been in existence for over fifty years, they have not been deployed extensively because they need to be operated at a high temperature (∼1000 °C), are expensive, and have slow response to changes in energy demands. One important need for commercialization of SOFCs is a lowering of their operating temperature, which requires an electrolyte that can operate at lower temperatures. Doped ceria is one such candidate. For this dissertation work I have studied different types of doped ceria to understand the mechanism of oxygen vacancy diffusion through the bulk. Doped ceria is important because they have high ionic conductivities thus making them attractive candidates for the electrolytes of solid oxide fuel cells. In particular, I have studied how the ionic conductivities are improved in these doped materials by studying the oxygen-vacancy formations and migrations. In this dissertation I describe the application of density functional theory (DFT) and Kinetic Lattice Monte Carlo (KLMC) simulations to calculate the vacancy diffusion and ionic conductivities in doped ceria. The dopants used are praseodymium (Pr), gadolinium (Gd), and neodymium (Nd), all belonging to the lanthanide series. The activation energies for vacancy migration between different nearest neighbor (relative to the dopant) positions were calculated using the commercial DFT code VASP (Vienna Ab-initio Simulation Package). These activation energies were then used as inputs to the KLMC code that I co-developed. The KLMC code was run for different temperatures (673 K to 1073 K) and for different dopant concentrations (0 to 40%). These simulations have resulted in the prediction of dopant concentrations for maximum ionic conductivity at a given temperature.