Applying the screened hybrid functional Heyd–Scuseria–Ernzerhof (HSE) method, we studied the polaronic degree of freedom of different charged oxygen vacancies V o in rutile TiO 2. The HSE method not only corrects the band gap, but also allows for correct polaron localization. Due to the important role of phonon in oxygen vacancy associated levels in the gap, we calculated configuration coordinate (CC) potential energy surfaces for all charged V o 's. Our calculated CC diagrams with effective impression on host states, show significant improvement of electron–lattice interaction compared to semi(local) DFT methods. The obtained values of stokes shifts for sequential transitions of charged vacancies agree well with experimental evidences which confirm Ti 3+ centers are responsible for photoluminescence. In addition, we explored the effect of polaron localization on diffusive mechanism of V o along most open [001] direction. Calculated values of migration barriers for [Formula: see text] are found to be in quantitative agreement with experimental migration energy [E. Iguchi and K. Yajima, J. Phys. Soc. Jpn.32 (1971) 1415] of 2.4 eV. These results highlight the small polaronic behavior of V o 's and is consistent with studies suggest the polaronic hopping model for electron transport of n-type conductivity in reduced TiO 2 [J.-F. Baumard and F. Gervais, Phys. Rev. B15 (1977) 2316–2323].