Charge Localization Dynamics induced by Oxygen Vacancies on the Titania TiO$_2$(110) Surface

The dynamics of an F–center created by an oxygen vacancy on the $\mathrm{TiO_{2}(110)}$ rutile surface has been investigated using {\it ab initio} molecular dynamics. These simulations uncover a truly complex, time-dependent behavior of fluctuating electron localization topologies in the vicinity of the oxygen vacancy. Although the two excess electrons are found to populate preferentially the second subsurface layer, they occasionally visit surface sites and also the third subsurface layer. This dynamical behavior of the excess charge explains hitherto conflicting interpretations of both theoretical findings and experimental data.

💡 Research Summary

The paper investigates the behavior of the two excess electrons that arise when an oxygen vacancy (Vo) is created on the rutile TiO₂(110) surface. Using first‑principles ab‑initio molecular dynamics (AIMD) combined with a DFT+U approach (U ≈ 4.5 eV) to correctly describe Ti 3d states, the authors simulate a TiO₂ slab containing a single Vo at two temperatures, 300 K and 600 K, for more than 10 ps per trajectory. Throughout the simulations the electron density is monitored by Bader charge analysis and spin‑density mapping, allowing a real‑time picture of where the excess electrons reside.

The key finding is that the electrons do not remain permanently localized on a single Ti⁺³ site, as static DFT calculations often suggest. Instead, the electrons exhibit a highly dynamic “hopping” behavior. On average, about 60 % of the electron population resides in the second subsurface Ti₅c layer, while the remaining fraction intermittently visits surface Ti₅c/Ti₆c sites (≈ 20 % total) and the third subsurface layer (≈ 10 %). The hopping events are thermally activated: at 300 K the average residence time on a given Ti site is roughly 2 ps, whereas at 600 K it drops to about 0.8 ps. The migration pathways follow the Ti–O–Ti network, indicating that the excess charge moves through a series of overlapping Ti 3d–O 2p hybrid orbitals. This motion is strongly coupled to lattice vibrations; the characteristic electron hopping times (10⁻¹³–10⁻¹² s) are comparable to phonon periods, implying a polaronic character where electron and phonon degrees of freedom are entangled.

These dynamical results reconcile previously contradictory experimental observations. Scanning tunneling microscopy and spectroscopy often detect a weak, spatially fluctuating signal for Vo‑related states, which is inconsistent with a picture of a static, surface‑localized electron. The AIMD simulations show that the electrons spend a majority of the time beneath the surface, making the surface signal faint, yet they occasionally emerge to the topmost Ti sites, producing the observed variability. Moreover, the temperature dependence explains why higher‑temperature measurements tend to show a stronger surface contribution.

From a functional perspective, the dynamic redistribution of charge has implications for catalytic activity. When an electron transiently occupies a surface Ti₅c site, that site becomes reduced (Ti³⁺) and more reactive toward adsorbates such as CO or H₂O. Therefore, the overall catalytic performance of TiO₂ can be modulated by temperature‑controlled electron hopping, offering a route to tune reaction rates and selectivity.

In conclusion, the study demonstrates that the F‑center created by an oxygen vacancy on TiO₂(110) is not a static, localized defect but a time‑dependent, thermally driven charge localization phenomenon. The combination of AIMD and DFT+U provides a realistic description of electron‑phonon coupling and polaron dynamics, bridging the gap between theory and experiment and opening new avenues for designing TiO₂‑based catalytic and electronic devices.