density functional theory (DFT) simulations in conjunction with experimental nuclear magnetic resonance (NMR) in order to investigate the Li ion migration mechanisms. Disorder in the anion sublattice can impact the material properties and the diffusion of point defects [14,15]. This is a common future in energy related materials where intrinsic disorder and doping can influence the formation (i.e. concentration of point defects mediating diffusion) and the migration (i.e. the energy barriers for diffusion) [16][17][18][19]. Controlling point defects such as vacancies is important for semiconductors, superconductors, and oxides [20][21][22][23][24][25][26][27][28]. A way to defect engineer the concentration and clustering of point defects in the anion sublattice (for example oxygen vacancies) is via the introduction of dopants. This is effectively to maintain charge balance in the lattice [29]. For example the introduction of two trivalent dopants in the tetravalent cerium site in CeO 2 can be charge balanced by the formation of an oxygen vacancy. Therefore the introduction of trivalent dopants in CeO 2 is an efficient way to form oxygen vacancies at concentrations higher than the equilibrium concentration [29].Atomistic simulation is an efficient and powerful way to understand the energetics of point defects in energy materials [30][31][32]. The main aim of the present study is to systematically investigate the intrinsic defect processes and impact of doping in Li 2 ZrO 3 using DFT. In particular, we consider here divalent (Mg, Zn, Ca, Cd, Sr, Ba), trivalent (Al, Ga, Sc, In, Y) and tetravalent (Si, Ge, Ti, Sn, Pb, Ce) substitutionals and their association with oxygen vacancies.Abstract Lithium zirconate (Li 2 ZrO 3 ) is an important material that is considered as an anode in lithium-ion batteries and as a nuclear reactor breeder material. The intrinsic defect processes and doping can impact its material properties. In the present study we employ density functional theory calculations to calculate the defect processes and doping in Li 2 ZrO 3 . Here we show that the lithium Frenkel is the dominant intrinsic defect process and that dopants substituting in the zirconium site strongly associate with oxygen vacancies. In particular, it is calculated that divalent dopants more strongly bind with oxygen vacancies, with trivalent dopants following in binding energies and even tetravalent dopands having significant binding energies. The results are discussed in view of the application of Li 2 ZrO 3 in energy applications.