The diffusion of Li in bulk Si and crystalline LiSi is investigated over a wide range of temperatures employing first-principles calculations based on density functional theory, transition state theory, and the kinetic Monte Carlo method. Nuclear quantum effects are incorporated by computing the vibrational spectrum and its effect on the effective energy barrier. The Li diffusion coefficient in bulk Si calculated with such quantum effects is ∼33% lower than the classical limit near room temperature due to higher effective energy barrier and tends to the classical limit at a high temperature (>1000 K). The presence of anharmonicity, estimated by the quasiharmonic approximation and the cBω model, increases the diffusion coefficient by ∼60%. For Li diffusion in LiSi with multiple vacancy jumps, we obtain an effective diffusion barrier of 0.27 eV ± 0.01 eV. In the Li-Si system, the quantum mechanical effects are only marginally significant at room temperature.