Modeling polarons in density functional theory: lessons learned from TiO2

Machine-Learning-Driven Exploration of Surface Reconstructions of Reduced Rutile TiO2

Titanium dioxide (TiO2) is widely used as a catalyst support due to its stability, tunable electronic properties, and surface oxygen vacancies, which are crucial for catalytic processes such as the reverse water-gas shift (RWGS) reaction. Reduced TiO2 surfaces undergo complex surface reconstructions that endow unique properties but are computationally challenging to describe. In this study, we utilize machine-learning interatomic potentials (MLIPs) integrated with an active-learning workflow to efficiently explore reduced rutile TiO2 surfaces. This approach enabled the prediction of a phase diagram as a function of oxygen chemical potential, revealing a variety of reconstructed phases, including a previously unreported subsurface shear plane structure. We further investigate the electronic properties of these surfaces and validate our results by comparing experimental and theoretical high-resolution transmission electron microscopy (HRTEM). Our findings provide new insights into how extreme surface reductions influence the structural and electronic properties of TiO2, with potential implications for catalyst design.

Assigning Surface Hole Polaron Configurations of Titanium Oxide Materials to Excited-State Optical Absorptions

For water splitting, a comprehensive understanding of the underlying reaction intermediates and pathways is crucial for optimizing catalyst design. Among the most well-known active photoanodes for the oxygen evolution half-reaction are TiO2-based materials. A hole polaron, which consists of a metal-oxide distortion around trapped holes, has been suggested as a local reactive oxygen configuration. While first-principles calculations identify new electronic states in the middle of the band gap and the influence of trapped hole dynamics on transport, an assignment of hole polaron configurations to a measured spectrum has been challenging due to broad optical transitions in the visible regime. Here, we compare the excited-state absorption (ESA) for two titanium oxide materials with a similar electronic structure but differing crystal structure. The ESA maximum for ultrafast time scales (<1 ps) is isolated by a principal component analysis and shifts from 3.1 eV in rutile TiO2 (100) to 2.2 eV in perovskite SrTiO3. Density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations predict the energies of the midgap states for stable hole polarons and their corresponding spectra. The shift in the ESA is rationalized by the transition optical dipole originating from both edge and deeper states in the valence band being bright for certain configurations of hole polarons in rutile TiO2 (100) (terminal O•-) versus STO (lateral Ti2O•-). The spectral assignment of a shifting ESA between two titanium oxide materials informs the assignment of hole polaron configurations for oxygen evolution catalysis and, more generally, photodriven processes.

Study of rutile TiO2(110) single crystal by transient absorption spectroscopy in the presence of Ce4+cations in aqueous environment. Implication on water splitting

Ce4+cations are commonly used as electron acceptors during the water oxidation to O2reaction over Ir- and Ru-based catalysts. They can also be reduced to Ce3+cations by excited electrons from the conduction band of an oxide semiconductor with a suitable energy level. In this work, we have studied their interaction with a rutile TiO2(110) single crystal upon band gap excitation by femtosecond transient absorption spectroscopy (TAS) in solution in the 350-900 nm range and up to 3.5 ns. Unlike excitation in the presence of water alone the addition of Ce4+resulted in a clear ground-state bleaching (GSB) signal at the band gap energy of TiO2(ca. 400 nm) with a time constantt= 4-5 ps. This indicated that the Ce4+cations presence has quenched the e-h recombination rate when compared to water alone. In addition to GSB, two positive signals are observed and are attributed to trapped holes (in the visible region, 450-550 nm) and trapped electrons in the IR region (>700 nm). Contrary to expectation, the lifetime of the positive signal between 450 and 550 nm decreased with increasing concentrations of Ce4+. We attribute the decrease in the lifetime of this signal to electrostatic repulsion between Ce4+at the surface of TiO2(110) and positively charged trapped holes. It was also found that at the very short time scale (<2-3 ps) the fast decaying TAS signal of excited electrons in the conduction band is suppressed because of the presence of Ce4+cations. Results point out that the presence of Ce4+cations increases the residence time (mobility) of excited electrons and holes at the conduction band and valence band energy levels (instead of being trapped). This might provide further explanations for the enhanced reaction rate of water oxidation to O2in the presence of Ce4+cations.

On the Structure of Oxygen Deficient Amorphous Oxide Films

Understanding defects in amorphous oxide films and heterostructures is vital to improving performance of microelectronic devices, thin-film transistors, and electrocatalysis. However, to what extent the structure and properties of point defects in amorphous solids are similar to those in the crystalline phase are still debated. The validity of this analogy and the experimental and theoretical evidence of the effects of oxygen deficiency in amorphous oxide films are critically discussed. The authors start with the meaning and significance of defect models, such as "oxygen vacancy" in crystalline oxides, and then introduce experimental and computational methods used to study intrinsic defects in amorphous oxides and discuss their limitations and challenges. To test the validity of existing defect models, ab initio molecular dynamics is used with a non-local density functional to model the structure and electronic properties of oxygen-deficient amorphous alumina. Unlike some previous studies, the formation of deep defect states in the bandgap caused by the oxygen deficiency is found. Apart from atomistic structures analogous to crystal vacancies, the formation of more stable defect states characterized by the bond formation between under-coordinated Al ions is shown. The limitations of such defect models and how they may be overcome in simulations are discussed.

Time-Domain View of Polaron Dynamics in Metal Oxide Photocatalysts

The polaron is a fundamental physical phenomenon in transition metal oxides (TMOs), and it has been studied extensively for decades. However, the implication of a polaron on photochemistry is still ambiguous. As such, understanding the fundamental properties and controlling the dynamics of polarons at the atomistic level is desired. In this Perspective, we seek to highlight the recent advances in studying small polarons in TMOs, with a particular focus on nonadiabatic molecular dynamics at the ab initio level, and discuss the implications for photocatalysis from the aspects of the structure, intrinsic physical properties, formation, migration, and recombination of small polarons. Finally, various methods were proposed to advance our understanding of manipulating the small-polaron dynamics, and strategies to design high-performance TMO-based photoelectrodes were discussed.

Competing electronic states emerging on polar surfaces

Excess charge on polar surfaces of ionic compounds is commonly described by the two-dimensional electron gas (2DEG) model, a homogeneous distribution of charge, spatially-confined in a few atomic layers. Here, by combining scanning probe microscopy with density functional theory calculations, we show that excess charge on the polar TaO2 termination of KTaO3(001) forms more complex electronic states with different degrees of spatial and electronic localization: charge density waves (CDW) coexist with strongly-localized electron polarons and bipolarons. These surface electronic reconstructions, originating from the combined action of electron-lattice interaction and electronic correlation, are energetically more favorable than the 2DEG solution. They exhibit distinct spectroscopy signals and impact on the surface properties, as manifested by a local suppression of ferroelectric distortions.

Accurate Determination of Adsorption-Energy Differences of Metalloporphyrins on Rutile TiO2(110) 1 × 1

Understanding the adsorption of organic molecules on surfaces is of essential importance for many applications. Adsorption energies are typically measured using temperature-programmed desorption. However, for large organic molecules, often only desorption of the multilayers is possible, while the bottom monolayer in direct contact to the surface cannot be desorbed without decomposition. Nevertheless, the adsorption energies of these directly adsorbed molecules are the ones of the most interest. We use a layer-exchange process investigated with X-ray photoelectron spectroscopy to compare the relative adsorption energies of several metalated tetraphenylporphyrins on rutile TiO₂(110) 1 × 1. We deposit a mixture of two different molecules, one on top of the other, and slowly anneal above their multilayer desorption temperature. During the slow heating, the molecules begin to diffuse between the layers and the molecules with the stronger interaction with the surface displace the weaker-interacting molecules from the surface and push them into the multilayer. The multilayers eventually desorb, leaving behind a monolayer of strongly interacting molecules. From the ratio of the two different porphyrin molecules in the residual monolayer and the desorbed multilayer, we can calculate the equilibrium constant of the layer-exchange process and thereby the difference in adsorption energy between the two different porphyrin molecules.