“…Furthermore, calculation of barriers to hole hopping are far less affected by finite size effects since the hole remains localized throughout transition. To model polaron hopping, we employ the MEHAM model, which is able to describe the charge-transfer processes in both the adiabatic and non-adiabatic limits and is well suited to describing small polaron hopping in TiO 2 . ,− There are many previous applications of the approach to model polaron hopping in TiO 2 , Fe 2 O 3 , and organic semiconductors as well as electron transfer between point defects in crystals. ,,− ,, Here, we employ MEHAM with the common approximation that the hole hopping process can be described as a one-dimensional process coupling to a single generalized phonon mode. We note that although we consider a single phonon mode, it includes displacements of all atoms and the wavefunction is not constrained in any way, so it is not the case that adjacent atoms are not involved (in general they all are).…”