Abstract
Ab initio molecular dynamics (MD) simulations and NMR GIPAW (gauge including projector augmented wave) calculations have been used to analyze the coordination and mobility of Li ions in LiTi(2) (PO(4) )(3) (rhombohedral), LiZr(2) (PO(4) )(3) (triclinic), and LiZr(2) (PO(4) )(3) (rhombohedral) phases. Significant discrepancies are observed between static calculations of (7) Li quadrupolar parameters and experimental values. The dynamical origin of this disagreement is demonstrated by incorporating in the calculations thermal vibrations and local motion of atoms with MD simulations. For LiTi(2) (PO(4) )(3) , the quadrupolar constant associated with Li ions grows with temperature because the local symmetry of the system decreases, whereas for the Zr phases, the quadrupolar constant decreases because thermal vibrations reduce the anisotropy of the interaction. Finally, for both Zr phases, MD yields Li distributions that compare well with disorder reported from diffraction studies.