Anisotropic Li diffusion in pristine and defective ZnO

We study the Li interstitial diffusion in pristine and defective ZnO bulk by means of first-principles density functional theory (DFT) coupled with Nudged Elastic Band (NEB) calculations. We consider three types of point defects, i.e., oxygen vacancy (Ovac), Zn vacancy (Znvac) and ZnO vacancy pair (ZnOvac-pair) and investigate their individual effect on the energy barrier of Li interstitial diffusion. Our results predict that Ovac and Znvac lower the Li diffusion energy barrier as compared to the pristine ZnO case. However, we further find that Li interstitial on the other hand may possibly be trapped inside the Znvac subsequently forming the LiZn substitutional type of defect. The similar behavior also observed for Li interstitial in the vicinity of Zn-Ovac_pair though with less change of Li diffusion barriers as compared to the other two cases. Our results indicate that among the three considered defects only Ovac shows possible enhancement of the kinetics of Li diffusion inside ZnO bulk.