We have used density-functional theory [generalized gradient approximation (GGA)] to study lithium intercalation at low concentration into anatase TiO2. To describe the defect states produced by Li doping a Hubbard “+U” correction is applied to the Ti d states (GGA+U). Uncorrected GGA calculations predict LixTiO2 to be metallic with the excess charge distributed over all Ti sites, whereas GGA+U predicts a defect state 0.96 eV below the conduction band, in agreement with experimental photoelectron spectra. This occupied defect state corresponds to charge strongly localized at a single Ti 3d site neighboring the intercalated lithium with a magnetization of 1 μB. This polaronic state produces a redshifted optical absorption spectrum, which is compared to those for the native O-vacancy and Ti-interstitial defects. The strong localization of charge at a single Ti center lowers the symmetry of the interstitial geometry relative to that predicted by GGA. The intercalated lithium sits close to the center of the octahedral site, occupying a single potential energy minimum with respect to displacement along the  direction. This challenges the previous interpretation of neutron diffraction data that there exist two potential energy minima separated by 1.6 Å along the  direction within each octahedron. Nudged elastic band calculations give barriers to interoctahedral diffusion of ∼0.6 eV, in good agreement with experimental data. These barrier heights are found to depend only weakly on the position of the donated electron. The intercalation energy is 2.14 eV with GGA and 1.88 eV with GGA+U, compared to the experimental value of ∼1.9 eV. Li-electron binding energies have also been calculated. The [Li•i-Ti′Ti] complex has a binding energy of 56 meV, and a second electron is predicted to be bound to give [Li•i-2Ti′Ti] with a stabilization energy of 30 meV, indicating that intercalated lithium will weakly trap excess electrons produced during photoillumination or introduced by additional n-type doping.