Prediction of Feasibility of Polaronic OER on (110) Surface of Rutile TiO2.

The polaronic effects at the atomic level hold paramount significance for advancing the efficacy of transition metal oxides in applications pertinent to renewable energy. The lattice-distortion mediated localization of photoexcited carriers in the form of polarons plays a pivotal role in the photocatalysis. This investigation focuses on rutile TiO2, an important material extensively explored for solar energy conversion in artificial photosynthesis, specifically targeting the generation of green H2 through photoelectrochemical (PEC) H2O splitting. By employing Hubbard-U corrected and hybrid density functional theory (DFT) methods, we systematically probe the polaronic effects in the catalysis of oxygen evolution reaction (OER) on the (110) surface of rutile TiO2. Theoretical understanding of polarons within the surface, coupled with simulations of OER at distinct titanium (Ti) and oxygen (O) active sites, reveals diverse polaron formation energies within the lattice sites with strong preference for bulk and surface bridge (Ob) oxygen sites. Moreover, we provide the evidence for the facilitative role of polarons in OER. We find that hole polarons situated at the equatorial oxygen sites near the Ti-active site, along with bridge site hole polarons distal from the Ob active site yield a small reduction in OER overpotential by ~0.06 eV and ~0.12 eV, respectively. However, subsurface, equatorial, and bridge site hole polarons significantly reduce the Ti-active site OER overpotential by ~0.4 eV through the peroxo-type oxygen pathway. We also observe that the presence of hole polarons stabilizes the *OH, *O, and *OOH intermediate species compared to the scenario without hole polarons. Overall, this study provides a detailed mechanistic insight into polaron-mediated OER, offering a promising avenue for improving the catalytic activity of transition metal oxide-based photocatalysts catering to renewable energy requisites.

Prediction of Feasibility of Polaronic OER on the (110) Surface of Rutile TiO2.

The front cover artwork is provided by Dr. Hori Pada Sarker from Dr. Frank Abild-Pedersen's research group at the SLAC National Accelerator Laboratory. The image shows the generation of photoexcited carriers (electrons and holes) and the subsequent formation of hole polaron in rutile TiO2 during oxygen evolution reaction (OER). Read the full text of the Research Article at 10.1002/cphc.202400060.

Interpretable Machine Learning Models for Practical Antimonate Electrocatalyst Performance.

Computationally predicting the performance of catalysts under reaction conditions is a challenging task due to the complexity of catalytic surfaces and their evolution in situ, different reaction paths, and the presence of solid-liquid interfaces in the case of electrochemistry. We demonstrate here how relatively simple machine learning models can be found that enable prediction of experimentally observed onset potentials. Inputs to our model are comprised of data from the oxygen reduction reaction on non-precious transition-metal antimony oxide nanoparticulate catalysts with a combination of experimental conditions and computationally affordable bulk atomic and electronic structural descriptors from density functional theory simulations. From human-interpretable genetic programming models, we identify key experimental descriptors and key supplemental bulk electronic and atomic structural descriptors that govern trends in onset potentials for these oxides and deduce how these descriptors should be tuned to increase onset potentials. We finally validate these machine learning predictions by experimentally confirming that scandium as a dopant in nickel antimony oxide leads to a desired onset potential increase. Macroscopic experimental factors are found to be crucially important descriptors to be considered for models of catalytic performance, highlighting the important role machine learning can play here even in the presence of small datasets.

Achieving industrial ammonia synthesis rates at near-ambient conditions through modified scaling relations on a confined dual site.

The production of ammonia through the Haber-Bosch process is regarded as one of the most important inventions of the 20th century. Despite significant efforts in optimizing the process, it still consumes 1 to 2% of the worldwide annual energy for the high working temperatures and pressures. The design of a catalyst with a high activity at milder conditions represents another challenge for this reaction. Herein, we combine density functional theory and microkinetic modeling to illustrate a strategy to facilitate low-temperature and -pressure ammonia synthesis through modified energy-scaling relationships using a confined dual site. Our results suggest that an ammonia synthesis rate two to three orders of magnitude higher than the commercial Ru catalyst can be achieved under the same reaction conditions with the introduction of confinement. Such strategies will open pathways for the development of catalysts for the Haber-Bosch process that can operate at milder conditions and present more economically viable alternatives to current industrial solutions.

Identification of earth-abundant materials for selective dehydrogenation of light alkanes to olefins.

Selective ethane dehydrogenation (EDH) is an attractive on-purpose strategy for industrial ethylene production. Design of an effective, stable, and earth-abundant catalyst to replace noble metal Pt is the main obstacle for its large-scale application. Herein, we report an experimentally validated theoretical framework to discover promising catalysts for EDH, which combines descriptor-based microkinetic modeling, high-throughput computations, machine-learning concepts, and experiments. Our approach efficiently evaluates 1,998 bimetallic alloys by using accurately calculated C and CH3 adsorption energies and identifies a small number of new promising noble-metal-free catalysts for selective EDH. A Ni3Mo alloy predicted to be promising is successfully synthesized, and experimentally proven to outperform Pt in selective ethylene production from EDH, representing an important contribution to the improvement of light alkane dehydrogenation to olefins. These results will provide essential additions in the discovery and application of earth-abundant materials in catalysis.

Toward Carbon Monoxide Methanation at Mild Conditions on Dual-Site Catalysts.

The catalytic carbon monoxide (CO) methanation is an ideal model reaction for the fundamental understanding of catalysis on the gas-solid interface and is crucial for various industrial processes. However, the harsh operating conditions make the reaction unsustainable, and the limitations set by the scaling relations between the dissociation energy barrier and dissociative binding energy of CO further increase the difficulty in designing high-performance methanation catalysts operating under milder conditions. Herein, we proposed a theoretical strategy to circumvent the limitations elegantly and achieve both facile CO dissociation and C/O hydrogenation on the catalyst containing a confined dual site. The DFT-based microkinetic modeling (MKM) reveals that the designed Co-Cr2/G dual-site catalyst could provide 4-6 orders of magnitude higher turnover frequency for CH4 production than the cobalt step sites. We believe that the proposed strategy in the current work will provide essential guidance for designing state-of-the-art methanation catalysts under mild conditions.

Revealing the structure of a catalytic combustion active-site ensemble combining uniform nanocrystal catalysts and theory insights.

Supported metal catalysts are extensively used in industrial and environmental applications. To improve their performance, it is crucial to identify the most active sites. This identification is, however, made challenging by the presence of a large number of potential surface structures that complicate such an assignment. Often, the active site is formed by an ensemble of atoms, thus introducing further complications in its identification. Being able to produce uniform structures and identify the ones that are responsible for the catalyst performance is a crucial goal. In this work, we utilize a combination of uniform Pd/Pt nanocrystal catalysts and theory to reveal the catalytic active-site ensemble in highly active propene combustion materials. Using colloidal chemistry to exquisitely control nanoparticle size, we find that intrinsic rates for propene combustion in the presence of water increase monotonically with particle size on Pt-rich catalysts, suggesting that the reaction is structure dependent. We also reveal that water has a near-zero or mildly positive reaction rate order over Pd/Pt catalysts. Theory insights allow us to determine that the interaction of water with extended terraces present in large particles leads to the formation of step sites on metallic surfaces. These specific step-edge sites are responsible for the efficient combustion of propene at low temperature. This work reveals an elusive geometric ensemble, thus clearly identifying the active site in alkene combustion catalysts. These insights demonstrate how the combination of uniform catalysts and theory can provide a much deeper understanding of active-site geometry for many applications.

Colloidal Platinum-Copper Nanocrystal Alloy Catalysts Surpass Platinum in Low-Temperature Propene Combustion.

Low-temperature removal of noxious environmental emissions plays a critical role in minimizing the harmful effects of hydrocarbon fuels. Emission-control catalysts typically consist of large quantities of rare, noble metals (e.g., platinum and palladium), which are expensive and environmentally damaging metals to extract. Alloying with cheaper base metals offers the potential to boost catalytic activity while optimizing the use of noble metals. In this work, we show that PtxCu100-x catalysts prepared from colloidal nanocrystals are more active than the corresponding Pt catalysts for complete propene oxidation. By carefully controlling their composition while maintaining nanocrystal size, alloys with dilute Cu concentrations (15-30% atomic fraction) demonstrate promoted activity compared to pure Pt. Complete propene oxidation was observed at temperatures as low as 150 °C in the presence of steam, and five to ten times higher turnover frequencies were found compared to monometallic Pt catalysts. Through DFT studies and structural and catalytic characterization, the remarkable activity of dilute PtxCu100-x alloys was related to the tuning of the electronic structure of Pt to reach optimal binding energies of C* and O* intermediates. This work provides a general approach toward investigation of structure-property relationships of alloyed catalysts with efficient and optimized use of noble metals.

Catalytic Performance and Near-Surface X-ray Characterization of Titanium Hydride Electrodes for the Electrochemical Nitrate Reduction Reaction.

The electrochemical nitrate reduction reaction (NO3RR) on titanium introduces significant surface reconstruction and forms titanium hydride (TiHx, 0 < x ≤ 2). With ex situ grazing-incidence X-ray diffraction (GIXRD) and X-ray absorption spectroscopy (XAS), we demonstrated near-surface TiH2 enrichment with increasing NO3RR applied potential and duration. This quantitative relationship facilitated electrochemical treatment of Ti to form TiH2/Ti electrodes for use in NO3RR, thereby decoupling hydride formation from NO3RR performance. A wide range of NO3RR activity and selectivity on TiH2/Ti electrodes between -0.4 and -1.0 VRHE was observed and analyzed with density functional theory (DFT) calculations on TiH2(111). This work underscores the importance of relating NO3RR performance with near-surface electrode structure to advance catalyst design and operation.

Atom-Specific Probing of Electron Dynamics in an Atomic Adsorbate by Time-Resolved X-Ray Spectroscopy.

The electronic excitation occurring on adsorbates at ultrafast timescales from optical lasers that initiate surface chemical reactions is still an open question. Here, we report the ultrafast temporal evolution of x-ray absorption spectroscopy (XAS) and x-ray emission spectroscopy (XES) of a simple well-known adsorbate prototype system, namely carbon (C) atoms adsorbed on a nickel [Ni(100)] surface, following intense laser optical pumping at 400 nm. We observe ultrafast (∼100 fs) changes in both XAS and XES showing clear signatures of the formation of a hot electron-hole pair distribution on the adsorbate. This is followed by slower changes on a few picoseconds timescale, shown to be consistent with thermalization of the complete C/Ni system. Density functional theory spectrum simulations support this interpretation.

Symmetry-resolved CO desorption and oxidation dynamics on O/Ru(0001) probed at the C K-edge by ultrafast x-ray spectroscopy.

We report on carbon monoxide desorption and oxidation induced by 400 nm femtosecond laser excitation on the O/Ru(0001) surface probed by time-resolved x-ray absorption spectroscopy (TR-XAS) at the carbon K-edge. The experiments were performed under constant background pressures of CO (6 × 10-8 Torr) and O2 (3 × 10-8 Torr). Under these conditions, we detect two transient CO species with narrow 2π* peaks, suggesting little 2π* interaction with the surface. Based on polarization measurements, we find that these two species have opposing orientations: (1) CO favoring a more perpendicular orientation and (2) CO favoring a more parallel orientation with respect to the surface. We also directly detect gas-phase CO2 using a mass spectrometer and observe weak signatures of bent adsorbed CO2 at slightly higher x-ray energies than the 2π* region. These results are compared to previously reported TR-XAS results at the O K-edge, where the CO background pressure was three times lower (2 × 10-8 Torr) while maintaining the same O2 pressure. At the lower CO pressure, in the CO 2π* region, we observed adsorbed CO and a distribution of OC-O bond lengths close to the CO oxidation transition state, with little indication of gas-like CO. The shift toward "gas-like" CO species may be explained by the higher CO exposure, which blocks O adsorption, decreasing O coverage and increasing CO coverage. These effects decrease the CO desorption barrier through dipole-dipole interaction while simultaneously increasing the CO oxidation barrier.

Modeling Potential-Dependent Electrochemical Activation Barriers: Revisiting the Alkaline Hydrogen Evolution Reaction.

Accurate theoretical simulation of electrochemical activation barriers is key to understanding electrocatalysis and guides the design of more efficient catalysts. Providing a detailed picture of proton transfer processes encounters several challenges: the constant potential requirement during charge transfer, the different time scales involved in the processes, and the thermal fluctuation of the solvent. Hence, it is prohibitively expensive computationally to apply density functional theory (DFT) calculations in modeling the potential-dependent activation barrier at the electrode-solvent interface, and the results are dubious. To address these challenges, we have developed an analytical approach based on charge conservation and decoupled potential energy surfaces to compute charge transfer barriers. The method makes it possible to simulate an electrochemical process at different potentials and explicitly include thermal fluctuations of the solvent at the electrode-solvent interface. We use the Pt-catalyzed alkaline hydrogen evolution reaction (HER) as our benchmark reaction, and we model the microkinetics of HER with consideration of the spatial fluctuations between the metal surface and the first solvent layer at room temperature. The distribution of water-metal distances has a large effect on the barriers of the charge transfer processes, and an accurate account of the statistical fluctuation in the reaction network leads to a several orders of magnitude increase in HER current as compared to transfer from a static solvent. The trends of the different reaction mechanisms in HER were successfully simulated with our model, and the theoretical I-V curves obtained are in good qualitative agreement with experimental results.

Ultrafast Adsorbate Excitation Probed with Subpicosecond-Resolution X-Ray Absorption Spectroscopy.

We use a pump-probe scheme to measure the time evolution of the C K-edge x-ray absorption spectrum from CO/Ru(0001) after excitation by an ultrashort high-intensity optical laser pulse. Because of the short duration of the x-ray probe pulse and precise control of the pulse delay, the excitation-induced dynamics during the first picosecond after the pump can be resolved with unprecedented time resolution. By comparing with density functional theory spectrum calculations, we find high excitation of the internal stretch and frustrated rotation modes occurring within 200 fs of laser excitation, as well as thermalization of the system in the picosecond regime. The ∼100 fs initial excitation of these CO vibrational modes is not readily rationalized by traditional theories of nonadiabatic coupling of adsorbates to metal surfaces, e.g., electronic frictions based on first order electron-phonon coupling or transient population of adsorbate resonances. We suggest that coupling of the adsorbate to nonthermalized electron-hole pairs is responsible for the ultrafast initial excitation of the modes.

Combining artificial intelligence and physics-based modeling to directly assess atomic site stabilities: from sub-nanometer clusters to extended surfaces.

The performance of functional materials is dictated by chemical and structural properties of individual atomic sites. In catalysts, for example, the thermodynamic stability of constituting atomic sites is a key descriptor from which more complex properties, such as molecular adsorption energies and reaction rates, can be derived. In this study, we present a widely applicable machine learning (ML) approach to instantaneously compute the stability of individual atomic sites in structurally and electronically complex nano-materials. Conventionally, we determine such site stabilities using computationally intensive first-principles calculations. With our approach, we predict the stability of atomic sites in sub-nanometer metal clusters of 3-55 atoms with mean absolute errors in the range of 0.11-0.14 eV. To extract physical insights from the ML model, we introduce a genetic algorithm (GA) for feature selection. This algorithm distills the key structural and chemical properties governing the stability of atomic sites in size-selected nanoparticles, allowing for physical interpretability of the models and revealing structure-property relationships. The results of the GA are generally model and materials specific. In the limit of large nanoparticles, the GA identifies features consistent with physics-based models for metal-metal interactions. By combining the ML model with the physics-based model, we predict atomic site stabilities in real time for structures ranging from sub-nanometer metal clusters (3-55 atom) to larger nanoparticles (147 to 309 atoms) to extended surfaces using a physically interpretable framework. Finally, we present a proof of principle showcasing how our approach can determine stable and active nanocatalysts across a generic materials space of structure and composition.

Predicting metal-metal interactions. I. The influence of strain on nanoparticle and metal adlayer stabilities.

Strain-engineering of bimetallic nanomaterials is an important design strategy for developing new catalysts. Herein, we introduce an approach for including strain effects into a recently introduced, density functional theory (DFT)-based alloy stability model. The model predicts adsorption site stabilities in nanoparticles and connects these site stabilities with catalytic reactivity and selectivity. Strain-based dependencies will increase the model's accuracy for nanoparticles affected by finite-size effects. In addition to the stability of small nanoparticles, strain also influences the heat of adsorption of epitaxially grown metal-on-metal adlayers. In this respect, we successfully benchmark the strain-including alloy stability model with previous experimentally determined trends in the heats of adsorption of Au and Cu adlayers on Pt (111). For these systems, our model predicts stronger bimetallic interactions in the first monolayer than monometallic interactions in the second monolayer. We explicitly quantify the interplay between destabilizing strain effects and the energy gained by forming new metal-metal bonds. While tensile strain in the first Cu monolayer significantly destabilizes the adsorption strength, compressive strain in the first Au monolayer has a minimal impact on the heat of adsorption. Hence, this study introduces and, by comparison with previous experiments, validates an efficient DFT-based approach for strain-engineering the stability, and, in turn, the catalytic performance, of active sites in bimetallic alloys with atomic level resolution.

Predicting metal-metal interactions. II. Accelerating generalized schemes through physical insights.

Operando-computational frameworks that integrate descriptors for catalyst stability within catalyst screening paradigms enable predictions of rates and selectivity on chemically faithful representations of nanoparticles under reaction conditions. These catalyst stability descriptors can be efficiently predicted by density functional theory (DFT)-based models. The alloy stability model, for example, predicts the stability of metal atoms in nanoparticles with site-by-site resolution. Herein, we use physical insights to present accelerated approaches of parameterizing this recently introduced alloy-stability model. These accelerated approaches meld quadratic functions for the energy of metal atoms in terms of the coordination number with linear correlations between model parameters and the cohesive energies of bulk metals. By interpolating across both the coordination number and chemical space, these accelerated approaches shrink the training set size for 12 fcc p- and d-block metals from 204 to as few as 24 DFT calculated total energies without sacrificing the accuracy of our model. We validate the accelerated approaches by predicting adsorption energies of metal atoms on extended surfaces and 147 atom cuboctahedral nanoparticles with mean absolute errors of 0.10 eV and 0.24 eV, respectively. This efficiency boost will enable a rapid and exhaustive exploration of the vast material space of transition metal alloys for catalytic applications.

Understanding Structure-Property Relationships of MoO3-Promoted Rh Catalysts for Syngas Conversion to Alcohols.

Rh-based catalysts have shown promise for the direct conversion of syngas to higher oxygenates. Although improvements in higher oxygenate yield have been achieved by combining Rh with metal oxide promoters, details of the structure of the promoted catalyst and the role of the promoter in enhancing catalytic performance are not well understood. In this work, we show that MoO3-promoted Rh nanoparticles form a novel catalyst structure in which Mo substitutes into the Rh surface, leading to both a 66-fold increase in turnover frequency and an enhancement in oxygenate yield. By applying a combination of atomically controlled synthesis, in situ characterization, and theoretical calculations, we gain an understanding of the promoter-Rh interactions that govern catalytic performance for MoO3-promoted Rh. We use atomic layer deposition to modify Rh nanoparticles with monolayer-precise amounts of MoO3, with a high degree of control over the structure of the catalyst. Through in situ X-ray absorption spectroscopy, we find that the atomic structure of the catalytic surface under reaction conditions consists of Mo-OH species substituted into the surface of the Rh nanoparticles. Using density functional theory calculations, we identify two roles of MoO3: first, the presence of Mo-OH in the catalyst surface enhances CO dissociation and also stabilizes a methanol synthesis pathway not present in the unpromoted catalyst; and second, hydrogen spillover from Mo-OH sites to adsorbed species on the Rh surface enhances hydrogenation rates of reaction intermediates.

Accessing the C-C transition state energy on transition metals.

The search for catalysts that can efficiently convert large hydrocarbons has been an active area of research for decades. To gain insight into those reactions, electronic structure calculations are playing an increasing role but the screening efforts are impeded by the complexity of the reaction networks that can contain hundreds of elementary steps, presenting a large number of computationally expensive transition state barrier calculations. A large number of the sub reactions in the network involve C-C bond dissociation, a step that has been identified as rate determining in many studies. The purpose of this article is to present a methodology that allows for accurate and rapid assessment of transition state energies for C-C bond breaking in any hydrocarbon based on a small number of simple calculations. Our model significantly enhances the capability of expanding the search space for new and efficient catalysts.

Understanding trends in C-H bond activation in heterogeneous catalysis.

While the search for catalysts capable of directly converting methane to higher value commodity chemicals and liquid fuels has been active for over a century, a viable industrial process for selective methane activation has yet to be developed. Electronic structure calculations are playing an increasingly relevant role in this search, but large-scale materials screening efforts are hindered by computationally expensive transition state barrier calculations. The purpose of the present letter is twofold. First, we show that, for the wide range of catalysts that proceed via a radical intermediate, a unifying framework for predicting C-H activation barriers using a single universal descriptor can be established. Second, we combine this scaling approach with a thermodynamic analysis of active site formation to provide a map of methane activation rates. Our model successfully rationalizes the available empirical data and lays the foundation for future catalyst design strategies that transcend different catalyst classes.

Nature of Lone-Pair-Surface Bonds and Their Scaling Relations.

We investigate the (surface) bonding of a class of industrially and biologically important molecules in which the chemically active orbital is a 2 p electron lone pair located on an N or O atom bound via single bonds to H or alkyl groups. This class includes water, ammonia, alcohols, ethers, and amines. Using extensive density functional theory (DFT) calculations, we discover scaling relations (correlations) among molecular binding energies of different members of this class: the bonding energetics of a single member can be used as a descriptor for other members. We investigate the bonding mechanism for a representative (H2O) and find the most important physical surface properties that dictate the strength and nature of the bonding through a combination of covalent and noncovalent electrostatic effects. We describe the importance of surface intrinsic electrostatic, geometric, and mechanical properties in determining the extent of the lone-pair-surface interactions. We study systems including ionic materials in which the surface positive and negative centers create strong local surface electric fields, which polarize the dangling lone pair and lead to a strong "electrostatically driven bond". We emphasize the importance of noncovalent electrostatic effects and discuss why a fully covalent picture, common in the current first-principles literature on surface bonding of these molecules, is not adequate to correctly describe the bonding mechanism and energy trends. By pointing out a completely different mechanism (charge transfer) as the major factor for binding N- and O-containing unsaturated (radical) adsorbates, we explain why their binding energies can be tuned independently from those of the aforementioned species, having potential implications in scaling-driven catalyst discovery.