Anisotropic Li diffusion in pristine and defective ZnO bulk and (10<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" altimg="si1.svg"><mml:mover accent="true"><mml:mn>1</mml:mn><mml:mo stretchy="true">¯</mml:mo></mml:mover></mml:math>0) surface

We study the Li interstitial diffusion in pristine and defective ZnO bulk and (10 1 ¯ 0) surface by means of first-principles density functional theory (DFT) coupled with the Nudged Elastic Band (NEB) calculations. We consider three types of point defects, i.e. , oxygen vacancy (O vac ), Zn vacancy (Zn vac ) and ZnO vacancy pair (ZnO vac-pair ) and investigate their individual effect on the energy barrier of Li interstitial diffusion. Our results predict that O vac and Zn vac 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 Zn vac subsequently forming the Li Zn substitutional type of defect. The similar behavior also observed for Li interstitial in the vicinity of Zn-O vac_pair though with less change of Li diffusion barriers as compared to the other two cases. For Li diffusion on the ZnO surface, our results indicate that O vac imposes similar effect as in the bulk case; it lowers the energy barriers for Li to diffuse outward (inward) from (to) ZnO sub-surface lattice. Our results indicate that among the three considered defects only O vac shows possible enhancement of the kinetics of Li diffusion inside the bulk and on the surface of ZnO. • O, Zn and Zn O pair vacancies effect on Li diffusion in ZnO was theoretically studied. • O vacancy possibly enhances the kinetic of Li mobility in ZnO. • Zn and Zn-O pair vacancies possibly hinder the Li mobility.