Modeling the noble metal/TiO2 (110) interface with hybrid DFT functionals: A periodic electrostatic embedded cluster model study

Theoretical Investigation of the Electrochemical Oxidation of H2 and CO Fuels on a Ruddlesden-Popper SrLaFeO4-δ Anode

The electrochemical oxidation of H₂ and CO fuels have been investigated on the Ruddlesden-Popper layered perovskite SrLaFeO_4-δ (SLF) under anodic solid oxide fuel cell conditions using periodic density functional theory and microkinetic modeling techniques. Two distinct FeO₂-plane-terminated surface models differing in terms of the underlying rock salt layer (SrO or LaO) are used to identify the active site and limiting factors for the electro-oxidation of H₂, CO, and syngas fuels. Microkinetic modeling predicted an order of magnitude higher turnover frequency for the electro-oxidation of H₂ compared to CO for SLF at short-circuit conditions. The surface model with an underlying SrO layer was found to be more active with respect to H₂ oxidation than the LaO-based surface model. At an operating voltage of less than 0.7 V, surface H₂O/CO₂ formation was found to be the key rate-limiting step, and the surface H₂O/CO₂ desorption was the key charge transfer step. In contrast, the bulk oxygen migration process was found to affect the overall rate at high cell voltage conditions above 0.9 V. In the presence of syngas fuel, the overall electrochemical activity is derived mainly from H₂ electro-oxidation and CO₂ is chemically shifted to CO via the reverse water-gas shift reaction. Substitutional doping of a surface Fe atom with Co, Ni, and Mn revealed that the H₂ electro-oxidation activity of FeO₂-plane terminated anodes with an underlying LaO rock salt layer can be improved with dopant introduction, with Co yielding a three orders of magnitude higher activity relative to the undoped LaO surface model. Constrained ab initio thermodynamic analysis furthermore suggested that the SLF anodes are resistant toward sulfur poisoning both in the presence and absence of dopants. Our findings reflect the role of various elements in controlling the fuel oxidation activity of SLF anodes that could aid the development of new Ruddlesden-Popper phase materials for fuel cell applications.

Aldehyde Hydrogenation by Pt/TiO2 Catalyst in Aqueous Phase: Synergistic Effect of Oxygen Vacancy and Solvent Water

The aldehyde hydrogenation for stabilizing and upgrading biomass is typically performed in aqueous phase with supported metal catalysts. By combining density functional theory calculations and ab initio molecular dynamics simulations, the model reaction of formaldehyde hydrogenation with a Pt/TiO₂ catalyst is investigated with explicit solvent water molecules. In aqueous phase, both the O vacancy (Ov) on support and solvent molecules could donate charges to a Pt cluster, where the Ov could dominantly reduce the Pt cluster from positive to negative. During the formaldehyde hydrogenation, the water molecules could spontaneously protonate the O in the aldehyde group by acid/base exchange, generating the OH* at the metal-support interface by long-range proton transfer. By comparing the stoichiometric and reduced TiO₂ support, it is found that the further hydrogenation of OH* is hard on the positively charged Pt cluster over stoichiometric TiO₂. However, with the presence of Ov on reduced support, the OH* hydrogenation could become not only exergonic but also kinetically more facile, which prohibits the catalyst from poisoning. This mechanism suggests that both the proton transfer from solvent water molecules and the easier OH* hydrogenation from Ov could synergistically promote aldehyde hydrogenation. That means, even for such simple hydrogenation in water, the catalytic mechanism could explicitly relate to all of the metal cluster, oxide support, and solvent waters. Considering the ubiquitous Ov defects in reducible oxide supports and the common aqueous environment, this synergistic effect may not be exclusive to Pt/TiO₂, which can be crucial for supported metal catalysts in biomass conversion.

Dynamic interplay between metal nanoparticles and oxide support under redox conditions

The dynamic interactions between noble metal particles and reducible metal-oxide supports can depend on redox reactions with ambient gases. Transmission electron microscopy revealed that the strong metal-support interaction (SMSI)-induced encapsulation of platinum particles on titania observed under reducing conditions is lost once the system is exposed to a redox-reactive environment containing oxygen and hydrogen at a total pressure of ~1 bar. Destabilization of the metal-oxide interface and redox-mediated reconstructions of titania lead to particle dynamics and directed particle migration that depend on nanoparticle orientation. A static encapsulated SMSI state was reestablished when switching back to purely oxidizing conditions. This work highlights the difference between reactive and nonreactive states and demonstrates that manifestations of the metal-support interaction strongly depend on the chemical environment.

General embedded cluster protocol for accurate modeling of oxygen vacancies in metal-oxides

The O vacancy (Ov) formation energy, EOv, is an important property of a metal-oxide, governing its performance in applications such as fuel cells or heterogeneous catalysis. These defects are routinely studied with density functional theory (DFT). However, it is well-recognized that standard DFT formulations (e.g., the generalized gradient approximation) are insufficient for modeling the Ov, requiring higher levels of theory. The embedded cluster method offers a promising approach to compute EOv accurately, giving access to all electronic structure methods. Central to this approach is the construction of quantum(-mechanically treated) clusters placed within suitable embedding environments. Unfortunately, current approaches to constructing the quantum clusters either require large system sizes, preventing application of high-level methods, or require significant manual input, preventing investigations of multiple systems simultaneously. In this work, we present a systematic and general quantum cluster design protocol that can determine small converged quantum clusters for studying the Ov in metal-oxides with accurate methods, such as local coupled cluster with single, double, and perturbative triple excitations. We apply this protocol to study the Ov in the bulk and surface planes of rutile TiO2 and rock salt MgO, producing the first accurate and well-converged determinations of EOv with this method. These reference values are used to benchmark exchange–correlation functionals in DFT, and we find that all the studied functionals underestimate EOv, with the average error decreasing along the rungs of Jacob’s ladder. This protocol is automatable for high-throughput calculations and can be generalized to study other point defects or adsorbates.

Understanding and controlling the formation of surface anion vacancies for catalytic applications

Systematic computational efforts aimed at calculating surface anion vacancy formation energies as important descriptors of catalytic performance are summarized.

Activation of Molecular O2 on CoFe2 O4 (001) Surfaces: An Embedded Cluster Study

Dioxygen activation pathways on the (001) surfaces of cobalt ferrite, CoFe2 O4 , were investigated computationally using density functional theory and the hybrid Perdew-Burke-Ernzerhof exchange-correlation functional (PBE0) within the periodic electrostatic embedded cluster model. We considered two terminations: the A-layer exposing Fe2+ and Co2+ metal sites in tetrahedral and octahedral positions, respectively, and the B-layer exposing octahedrally coordinated Co3+ . On the A-layer, molecular oxygen is chemisorbed as a superoxide on the Fe monocenter or bridging a Fe-Co cation pair, whereas on the B-layer it is adsorbed at the most stable anionic vacancy. Activation is promoted by transfer of electrons provided by the d metal centers onto the adsorbed oxygen. The subsequent dissociation of dioxygen into monoatomic species and surface reoxidation have been identified as the most critical steps that may limit the rate of the oxidation processes. Of the reactive metal-O species, [FeIII -O]2+ is thermodynamically most stable, while the oxygen of the Co-O species may easily migrate across the A-layer with barriers smaller than the associative desorption.

Core-Level Binding Energy Reveals Hydrogen Bonding Configurations of Water Adsorbed on TiO_{2}(110) Surface

Using x-ray photoelectron spectroscopy of the oxygen 1s core level, the ratio between intact (D_{2}O) and dissociated (OD) water in the hydrated stoichiometric TiO_{2}(110) surface is determined at varying coverage and temperature. In the submonolayer regime, both the D_{2}O∶OD ratio and the core-level binding energy of D_{2}O (ΔBE) decrease with temperature. The observed variations in ΔBE are shown with density functional theory to be governed crucially and solely by the local hydrogen bonding environment, revealing a generally applicable classification and details about adsorption motifs.

Potential and pH Dependence of the Buried Interface of Membrane-Coated Electrocatalysts

Semipermeable silica membranes are attractive as protective coatings for metal electrocatalysts such as platinum, but their impact on the catalytic properties has not been fully understood. Here, we develop a first-principles formalism to investigate how silica membranes interact with the surface of platinum metal electrocatalysts to develop a better understanding of the membrane-metal interplay. By generalizing the concept of Pourbaix diagrams to electrochemical solid-solid interfaces, we establish which bonds are formed between the SiO₂ membrane and Pt(111) surface in aqueous electrolytes for different pH values and potential biases. We find that the membrane termination changes as a function of the pH and potential, which affects the adhesion strength and the energy requirements for partial membrane detachment, controlling the Pt surface area that is accessible for reactant species. The charge transfer between the Pt surface and SiO₂ membrane is also pH- and potential-dependent and results in changes of the Pt surface d-band states, which are known to correlate with catalytic activity. Our analysis reveals the complex response of a buried interface to the electrochemical environment and identifies trends that are expected to apply also to other membrane-coated electrocatalysts.

Final-State Simulations of Core-Level Binding Energies at Metal-Organic Hybrid Interfaces: Artifacts Caused by Spurious Collective Electrostatic Effects

Core-level energies are frequently calculated to explain the X-ray photoelectron spectra of metal-organic hybrid interfaces. The current paper describes how such simulations can be flawed when modeling interfaces between physisorbed organic molecules and metals. The problem occurs when applying periodic boundary conditions to correctly describe extended interfaces and simultaneously considering core hole excitations in the framework of a final-state approach to account for screening effects. Since the core hole is generated in every unit cell, an artificial dipole layer is formed. In this work, we study methane on an Al(100) surface as a deliberately chosen model system for hybrid interfaces to evaluate the impact of this computational artifact. We show that changing the supercell size leads to artificial shifts in the calculated core-level energies that can be well beyond 1 eV for small cells. The same applies to atoms at comparably large distances from the substrate, encountered, for example, in extended, upright-standing adsorbate molecules. We also argue that the calculated work function change due to a core-level excitation can serve as an indication for the occurrence of such an artifact and discuss possible remedies for the problem.

Oxidative dehydrogenation of propane on the oxygen adsorbed edges of boron nitride nanoribbons

A N₂O/NOx-type site on the boron nitride nanoribbon edge is proposed as highly selective, active site for the oxidative dehydrogenation of propane.

σ-Holes and σ-lumps direct the Lewis basic and acidic interactions of noble metal nanoparticles: introducing regium bonds

Using local DFT-based probes for electrostatic as well as charge transfer/polarization interactions, we are able to characterize Lewis basic and acidic sites on copper, silver and gold nanoparticles.

The local electronic properties of individual Pt atoms adsorbed on TiO2(110) studied by Kelvin probe force microscopy and first-principles simulations

Noble metal nanostructures dispersed on metal oxide surfaces have applications in diverse areas such as catalysis, chemical sensing, and energy harvesting. Their reactivity, chemical selectivity, stability, and light absorption properties are controlled by the interactions at the metal/oxide interface. Single-atom metal adsorbates on the rutile TiO₂(110)-(1 × 1) surface have become a paradigmatic model to characterize those interactions and to understand the unique electronic properties of these supported nanostructures. We combine Kelvin probe force microscopy (KPFM) experiments and density functional theory (DFT) calculations to investigate the atomic-scale variations in the contact potential difference of individual Pt atoms adsorbed on a hydroxylated (h) TiO₂(110)-(1 × 1) surface. Our experiments show a significant drop in the local contact potential difference (LCPD) over Pt atoms with respect to the TiO₂ surface, supporting the presence of an electron transfer from the Pt adsorbates to the substrate. We have identified two characteristic regimes by LCPD spectroscopy. At far tip-sample distances, LCPD values show a weak distance dependence and can be attributed to the intrinsic charge transfer from Pt to the oxide support. Beyond the onset of short-range chemical interactions, LCPD values exhibit a strong distance dependence that we ascribe to the local structural and charge rearrangements induced by the tip-sample interaction. These findings also apply to other electropositive adsorbates such as potassium and the hydrogen atoms forming the OH groups that are present on the h-TiO₂(110) surface, promoting KPFM as a suitable tool for the understanding of electron transfer in catalytically active materials.

First-Principles Free-Energy Barriers for Photoelectrochemical Surface Reactions: Proton Abstraction at TiO_{2}(110)

We explicitly calculate the free-energy barrier for the initial proton abstraction in the water splitting reaction at rutile TiO_{2}(110) through ab initio molecular dynamics. Combining solid-state embedding, an energy based reaction coordinate and state-of-the-art free-energy reconstruction techniques renders the calculation tractable at the hybrid density-functional theory level. The obtained free-energy barrier of approximately 0.2 eV, depending slightly on the orientation of the first acceptor water molecule, suggests a hindered reaction on the pristine rutile surface.

Embedded-cluster calculations in a numeric atomic orbital density-functional theory framework

We integrate the all-electron electronic structure code FHI-aims into the general ChemShell package for solid-state embedding quantum and molecular mechanical (QM/MM) calculations. A major undertaking in this integration is the implementation of pseudopotential functionality into FHI-aims to describe cations at the QM/MM boundary through effective core potentials and therewith prevent spurious overpolarization of the electronic density. Based on numeric atomic orbital basis sets, FHI-aims offers particularly efficient access to exact exchange and second order perturbation theory, rendering the established QM/MM setup an ideal tool for hybrid and double-hybrid level density functional theory calculations of solid systems. We illustrate this capability by calculating the reduction potential of Fe in the Fe-substituted ZSM-5 zeolitic framework and the reaction energy profile for (photo-)catalytic water oxidation at TiO2(110).

The effect of the morphology of supported subnanometer Pt clusters on the first and key step of CO2 photoreduction

Other than the chemistry of the Pt cluster, the cluster morphology also determines CO₂ binding, attributed to structural fluxionality and bonding competitions among Pt atoms and CO₂.

First principles calculations of the structure and V L-edge X-ray absorption spectra of V2O5 using local pair natural orbital coupled cluster theory and spin-orbit coupled configuration interaction approaches

A detailed study of the electronic and geometric structure of V2O5 and its X-ray spectroscopic properties is presented. Cluster models of increasing size were constructed in order to represent the surface and the bulk environment of V2O5. The models were terminated with hydrogen atoms at the edges or embedded in a Madelung field. The structure and interlayer binding energies were studied with dispersion-corrected local, hybrid and double hybrid density functional theory as well as the local pair natural orbital coupled cluster method (LPNO-CCSD). Convergence of the results with respect to cluster size was achieved by extending the model to up to 20 vanadium centers. The O K-edge and the V L2,3-edge NEXAFS spectra of V2O5 were calculated on the basis of the newly developed Restricted Open shell Configuration Interaction with Singles (DFT-ROCIS) method. In this study the applicability of the method is extended to the field of solid-state catalysis. For the first time excellent agreement between theoretically predicted and experimentally measured vanadium L-edge NEXAFS spectra of V2O5 was achieved. At the same time the agreement between experimental and theoretical oxygen K-edge spectra is also excellent. Importantly, the intensity distribution between the oxygen K-edge and vanadium L-edge spectra is correctly reproduced, thus indicating that the covalency of the metal-ligand bonds is correctly described by the calculations. The origin of the spectral features is discussed in terms of the electronic structure using both quasi-atomic jj coupling and molecular LS coupling schemes. The effects of the bulk environment driven by weak interlayer interactions were also studied, demonstrating that large clusters are important in order to correctly calculate core level absorption spectra in solids.

Reactivity of metal oxide clusters with hydrogen peroxide and water--a DFT study evaluating the performance of different exchange-correlation functionals

We have performed a density functional theory (DFT) investigation of the interactions of H2O2, H2O and HO radicals with clusters of ZrO2, TiO2 and Y2O3. Different modes of H2O adsorption onto the clusters were studied. In almost all the cases the dissociative adsorption is more exothermic than molecular adsorption. At the surfaces where H2O has undergone dissociative adsorption, the adsorption of H2O2 and the transition state for its decomposition are mediated by hydrogen bonding with the surface HO groups. Using the functionals B3LYP, B3LYP-D and M06 with clusters of 26 and 8 units of ZrO2, the M06 functional performed better than B3LYP in describing the reaction of decomposition of H2O2 and the adsorption of H2O. Additionally, we investigated clusters of the type (ZrO2)2, (TiO2)2 and (Y2O3) and the performance of the functionals B3LYP, B3LYP-D, B3LYP*, M06, M06-L, PBE0, PBE and PWPW91 in describing H2O2, H2O and HO˙ adsorption and the energy barrier for decomposition of H2O2. The trends obtained for HO˙ adsorption onto the clusters are discussed in terms of the ionization energy of the metal cation present in the oxide. In order to correctly account for the existence of an energy barrier for the decomposition of H2O2, the functional used must include Hartree-Fock exchange. Using minimal cluster models, the best performance in describing the energy barrier for H2O2 decomposition was obtained with the M06 and PBE0 functionals - the average absolute deviations from experiments are 6 kJ mol(-1) and 5 kJ mol(-1) respectively. With the M06 functional and a larger monoclinic (ZrO2)8 cluster model, the performance is in excellent agreement with experimental data. For the different oxides, PBE0 was found to be the most effective functional in terms of performance and computational time cost.

The role of reducible oxide-metal cluster charge transfer in catalytic processes: new insights on the catalytic mechanism of CO oxidation on Au/TiO2 from ab initio molecular dynamics

To probe metal particle/reducible oxide interactions density functional theory based ab initio molecular dynamics studies were performed on a prototypical metal cluster (Au20) supported on reducible oxides (rutile TiO2(110)) to implicitly account for finite temperature effects and the role of excess surface charge in the metal oxide. It is found that the charge state of the Au particle is negative in a reducing chemical environment whereas in the presence of oxidizing species coadsorbed to the oxide surface the cluster obtained a net positive charge. In the context of the well-known CO oxidation reaction, charge transfer facilitates the plasticization of Au20, which allows for a strong adsorbate induced surface reconstruction upon addition of CO leading to the formation of mobile Au-CO species on the surface. The charging/discharging of the cluster during the catalytic cycle of CO oxidation enhances and controls the amount of O2 adsorbed at oxide/cluster interface and strongly influences the energetics of all redox steps in catalytic conversions. A detailed comparison of the current findings with previous studies is presented, and generalities about the role of surface-adsorbate charge transfer for metal cluster/reducible oxide interactions are discussed.

Displacement of hexanol by the hexanoic acid overoxidation product in alcohol oxidation on a model supported palladium nanoparticle catalyst

This work characterizes the adsorption, structure, and binding mechanism of oxygenated organic species from cyclohexane solution at the liquid/solid interface of optically flat alumina-supported palladium nanoparticle surfaces prepared by atomic layer deposition (ALD). The surface-specific nonlinear optical vibrational spectroscopy, sum-frequency generation (SFG), was used as a probe for adsorption and interfacial molecular structure. 1-Hexanoic acid is an overoxidation product and possible catalyst poison for the aerobic heterogeneous oxidation of 1-hexanol at the liquid/solid interface of Pd/Al(2)O(3) catalysts. Single component and competitive adsorption experiments show that 1-hexanoic acid adsorbs to both ALD-prepared alumina surfaces and alumina surfaces with palladium nanoparticles, that were also prepared by ALD, more strongly than does 1-hexanol. Furthermore, 1-hexanoic acid adsorbs with conformational order on ALD-prepared alumina surfaces, but on surfaces with palladium particles the adsorbates exhibit relative disorder at low surface coverage and become more ordered, on average, at higher surface coverage. Although significant differences in binding constant were not observed between surfaces with and without palladium nanoparticles, the palladium particles play an apparent role in controlling adsorbate structures. The disordered adsorption of 1-hexanoic acid most likely occurs on the alumina support, and probably results from modification of binding sites on the alumina, adjacent to the particles. In addition to providing insight on the possibility of catalyst poisoning by the overoxidation product and characterizing changes in its structure that result in only small adsorption energy changes, this work represents a step toward using surface science techniques that bridge the complexity gap between fundamental studies and realistic catalyst models.