Spontaneous Charge Separation at the Metal-Water Interface.

Reactions at the metal-water interface are essential in a range of fundamental and technological processes. Using Density Functional Theory calculations, we demonstrate that water substantially affects the adsorption of H and O2 on Cu(111), Ag(111), Au(111), Pd(111) and Pt(111). In water, H is found to undergo a spontaneous charge separation, where a proton desorbs to the water solution while an electron is donated to the surface. The reaction is exothermic over Au and Pt and associated with low barriers. The process is facile also over Pd, albeit slightly endothermic. For O2, water is found to increase the metal-to-adsorbate charge transfer, enhancing the adsorption energy and O-O bond length as compared to the adsorption in the absence of water. The magnitudes of the effects are system dependent, which implies that calculations should treat water explicitly. The results elucidate previous experimental results and highlights the importance of charge-transfer effects at the metal-water interface; both to describe the potential energy landscape, and to account for alternative reaction routes in the presence of water.

Kinetic Monte Carlo Simulations of Low-temperature NH3-SCR over Cu-exchanged Chabazite.

Cu-exchanged chabazite (Cu-CHA) is widely applied for ammonia assisted selective catalytic reduction of nitrogen oxides (NH3-SCR). The Cu+ ions are at low temperatures solvated by NH3 forming mobile [Cu(NH3)2]+ complexes. The dynamic behaviour of the complexes is critical as O2 adsorption requires a pair of complexes to form a [Cu2(NH3)4O2]2+ peroxo-species over which NO couples with NH3. Here we introduce a first principles-based kinetic Monte Carlo approach to explore the effect of the Al-distribution on the reaction kinetics of NH3-SCR over Cu-CHA. The method allows us to scrutinize the interplay between the pairing of [Cu(NH3)2]+ complexes and the reaction landscape for the NH3-SCR reaction over the peroxo-complex. The Al-distribution affects the stability of the [Cu(NH3)2]+ pairs as well as the kinetic parameters of the SCR-reaction. The turn-over frequency is determined by the stability of the [Cu(NH3)2]+ pairs and the relative strength of NO and NH3 adsorption once a pair is present. The results establish the hierarchy of effects that influences the performance of Cu-CHA over NH3-SCR and provide a computational basis for further development of the Cu-CHA material.

Site Communication in Direct Formation of H2O2 over Single-Atom Pd@Au Nanoparticles.

Single atom alloy catalysts offer possibilities to obtain turnover frequencies and selectivities unattainable by their monometallic counterparts. One example is direct formation of H2O2 from O2 and H2 over Pd embedded in Au hosts. Here, a first-principles-based kinetic Monte Carlo approach is developed to investigate the catalytic performance of Pd embedded in Au nanoparticles in an aqueous solution. The simulations reveal an efficient site separation where Pd monomers act as active centers for H2 dissociation, whereas H2O2 is formed over undercoordinated Au sites. After dissociation, atomic H may undergo an exothermic redox reaction, forming a hydronium ion in the solution and a negative charge on the surface. H2O2 is preferably formed from reactions between dissolved H+ and oxygen species on the Au surface. The simulations show that tuning the nanoparticle composition and reaction conditions can enhance the selectivity toward H2O2. The outlined approach is general and applicable for a range of different hydrogenation reactions over single atom alloy nanoparticles.

Surface steps dominate the water formation on Pd(111) surfaces.

Water formation is relevant in many technological processes and is also an important model reaction. Although water formation over Pd surfaces is widely studied, questions regarding the active site and the main reaction path (OH* + OH*) or (OH* + H*) are still open. Combining first-principles density functional theory calculations and kinetic Monte Carlo simulations, we find that the reaction rate is dominated by surface steps and point defects over a wide range of conditions. The main reaction path is found to be temperature dependent where the OH* + OH* path dominates at low temperatures, whereas the OH* + H* path is the main path at high temperatures. Steps facilitate the OH* formation, which is the rate limiting step under all conditions. OH* is formed via O* + H* association or OOH* splitting at low temperatures, whereas OH* is exclusively formed via O* + H* association at high temperatures. The results of the first-principles-based kinetic model are in excellent agreement with experimental observations at high and low temperatures as well as different gas-phase compositions.

Multiple Roles of Alkanethiolate-Ligands in Direct Formation of H2 O2 over Pd Nanoparticles.

Coadsorbed organic species including thiolates can promote direct synthesis of hydrogen peroxide from H2 and O2 over Pd particles. Here, density functional theory based kinetic modeling, augmented with activity measurements and vibrational spectroscopy are used to provide atomistic understanding of direct H2 O2 formation over alkylthiolate(RS) Pd. We find that the RS species are oxidized during reaction conditions yielding RSO2 as the effective ligand. The RSO2 ligand shows superior ability for proton transfer to the intermediate surface species OOH, which accelerates the formation of H2 O2 . The ligands promote the selectivity also by blocking sites for unselective water formation and by modifying the electronic structure of Pd. The work rationalizes observations of enhanced selectivity of direct H2 O2 formation over ligand-funtionalized Pd nanoparticles and shows that engineering of organic surface modifiers can be used to promote desired hydrogen transfer routes.

Structure-Dependent Strain Effects.

Density functional theory calculations of atomic and molecular adsorption on (111) and (100) metal surfaces reveal marked surface and structure dependent effects of strain. Adsorption in three-fold hollow sites is found to be destabilized by compressive strain whereas the reversed trend is commonly valid for adsorption in four-fold sites. The effects, which are qualitatively explained using a simple two-orbital model, provide insights on how to modify chemical properties by strain design.

Hydrogen adsorption on In2O3(111) and In2O3(110).

In2O3-Based catalysts have been measured to have a high activity for CO2 hydrogenation to H3COH. Here, we use density functional theory calculations with and without Hubbard-U corrections in combination with ab initio thermodynamics to investigate the dissociative adsorption of H2 over In2O3(111) and In2O3(110). H2 is found to dissociate heterolytically with a moderate barrier on both facets. Diffusion of hydrogen leads to the preferred homolytic adsorption configuration. Vacancy formation by water formation is thermodynamically preferred at high hydrogen coverages. Both surfaces are found to be hydroxylated at typical reaction conditions with the highest coverage predicted for In2O3(110). O 1s core level shifts are calculated for different coverages. The hydroxylated surfaces show two distinct shifts corresponding to different types of OH-groups. The presence of surface oxygen vacancies is not visible in the O 1s signatures. The results show that hydroxylation of the surfaces results in changes of the oxidation state of In-ions, which suggests that the redox properties on In2O3 are important for catalytic reduction of CO2 to added value chemicals.

Structure of two-dimensional Fe3O4.

We have investigated the structure of an ultrathin iron oxide phase grown on Ag(100) using surface x-ray diffraction in combination with Hubbard-corrected density functional theory (DFT+U) calculations. The film exhibits a novel structure composed of one close-packed layer of octahedrally coordinated Fe2+ sandwiched between two close-packed layers of tetrahedrally coordinated Fe3+ and an overall stoichiometry of Fe3O4. As the structure is distinct from bulk iron oxide phases and the coupling with the silver substrate is weak, we propose that the phase should be classified as a metastable two-dimensional oxide. The chemical and physical properties are potentially interesting, thanks to the predicted charge ordering between atomic layers, and analogy with bulk ferrite spinels suggests the possibility of synthesis of a whole class of two-dimensional ternary oxides with varying electronic, optical, and chemical properties.

Selective Acetylene Hydrogenation over Single-Atom Alloy Nanoparticles by Kinetic Monte Carlo.

Single-atom alloys, which are prepared by embedding isolated metal sites in host metals, are promising systems for improved catalyst selectivity. For technical applications, catalysts based on nanoparticles are preferred thanks to a large surface area. Herein, we investigate hydrogenation of acetylene to ethylene using kinetic Monte Carlo simulations based on density functional theory and compare the performance of Pd/Cu nanoparticles with Pd(111) and Pd/Cu(111). We find that embedding Pd in Cu systems strongly enhances the selectivity and that the reaction mechanism is fundamentally different for nanoparticles and extended surfaces. The reaction mechanism on nanoparticles is complex and involves elementary steps that proceed preferentially over different sites. Edge and corner sites on nanoparticles are predicted to lower the selectivity, and we infer that a rational design strategy in selective acetylene hydrogenation is to maximize the number of (111) sites in relation to edge sites for Pd/Cu nanoparticles.

CO2 adsorption on hydroxylated In2O3(110).

Catalytic synthesis of methanol from CO2 is one route to produce added-value chemicals from a greenhouse gas. Here, density functional theory calculations and ab initio thermodynamics are used to study CO2 adsorption on In2O3(110) in the presence of H2 and H2O. We find that the surface is heavily hydroxylated by either H2 or H2O and that hydroxylation promotes H2-induced vacancy formation. Moreover, CO2 adsorbs rather in a CO2- configuration on hydroxylated In2O3(110) than on oxygen vacancy sites. The results suggest that hydroxylation-induced oxidation-state changes of In-ions play a significant role in CO2 adsorption and activation during methanol synthesis.

A comparative test of different density functionals for calculations of NH3-SCR over Cu-Chabazite.

A general challenge in density functional theory calculations is to simultaneously account for different types of bonds. One such example is reactions in zeolites where both van der Waals and chemical bonds should be described accurately. Here, we use different exchange-correlation functionals to explore O2 dissociation over pairs of Cu(NH3)2+ complexes in Cu-Chabazite. This is an important part of selective catalytic reduction of NOx using NH3 as a reducing agent. The investigated functionals are PBE, PBE+U, PBE+D, PBE+U+D, PBE-cx, BEEF and HSE06+D. We find that the potential energy landscape for O2 activation and dissociation depends critically on the choice of functional. However, the van der Waals contributions are similarly described by the functionals accounting for this interaction. The discrepancies in the potential energy surface are instead related to different descriptions of the Cu-O chemical bond. By investigating the electronic, structural and energetic properties of reference systems including bulk copper oxides and (Cu2O2)2+ enzymatic crystals, we find that the PBE+U approach together with van der Waals corrections provides a reasonable simultaneous accuracy of the different bonds in the systems.

Steps Control the Dissociation of CO2 on Cu(100).

CO2 reduction reactions, which provide one route to limit the emission of this greenhouse gas, are commonly performed over Cu-based catalysts. Here, we use ambient pressure X-ray photoelectron spectroscopy together with density functional theory to obtain an atomistic understanding of the dissociative adsorption of CO2 on Cu(100). We find that the process is dominated by the presence of steps, which promote both a lowering of the dissociation barrier and an efficient separation between adsorbed O and CO, reducing the probability for recombination. The identification of steps as sites for efficient CO2 dissociation provides an understanding that can be used in the design of future CO2 reduction catalysts.

A Chemical View on X-ray Photoelectron Spectroscopy: the ESCA Molecule and Surface-to-Bulk XPS Shifts.

In this paper we remind the reader of a simple, intuitive picture of chemical shifts in X-ray photoelectron spectroscopy (XPS) as the difference in chemical bonding between the probed atom and its neighbor to the right in the periodic table, the so called Z+1 approximation. We use the classical ESCA molecule, ethyl trifluoroacetate, and 4d-transition metals to explicitly demonstrate agreement between core-level shifts computed as differences between final core-hole states and the approach where each core-ionized atom is replaced by a Z+1 atom. In this final state, or total energy picture, the XPS shift arises due to the more or less unfavorable chemical bonding of the effective nitrogen in the carbon geometry for the ESCA molecule. Surface core level shifts in metals are determined by whether the Z+1 atom as an alloy segregates to the surface or is more soluble in the bulk. As further illustration of this more chemical picture, we compare the geometry of C 1s and O 1s core-ionized CO with that of, respectively, NO+ and CF+ . The scope is not to propose a new method to compute XPS shifts but rather to stress the validity of this simple interpretation.

Se-C Cleavage of Hexane Selenol at Steps on Au(111).

Selenols are considered as an alternative to thiols in self-assembled monolayers, but the Se-C bond is one limiting factor for their usefulness. In this study, we address the stability of the Se-C bond by a combined experimental and theoretical investigation of gas-phase-deposited hexane selenol (CH3(CH2)5SeH) on Au(111) using photoelectron spectroscopy, scanning tunneling microscopy, and density functional theory (DFT). Experimentally, we find that initial adsorption leaves atomic Se on the surface without any carbon left on the surface, whereas further adsorption generates a saturated selenolate layer. The Se 3d component from atomic Se appears at 0.85 eV lower binding energy than the selenolate-related component. DFT calculations show that the most stable structure of selenols on Au(111) is in the form of RSe-Au-SeR complexes adsorbed on the unreconstructed Au(111) surface. This is similar to thiols on Au(111). Calculated Se 3d core-level shifts between elemental Se and selenolate in this structure nicely reproduce the experimentally recorded shifts. Dissociation of RSeH and subsequent formation of RH are found to proceed with high barriers on defect-free Au(111) terraces, with the highest barrier for scissoring R-Se. However, at steps, these barriers are considerably lower, allowing for Se-C bond breaking and hexane desorption, leaving elemental Se at the surface. Hexane is formed by replacing the Se-C bond with a H-C bond by using the hydrogen liberated from the selenol to selenolate transformation.

MonteCoffee: A programmable kinetic Monte Carlo framework.

Kinetic Monte Carlo (kMC) is an essential tool in heterogeneous catalysis enabling the understanding of dominant reaction mechanisms and kinetic bottlenecks. Here we present MonteCoffee, which is a general-purpose object-oriented and programmable kMC application written in python. We outline the implementation and provide examples on how to perform simulations of reactions on surfaces and nanoparticles and how to simulate sorption isotherms in zeolites. By permitting flexible and fast code development, MonteCoffee is a valuable alternative to previous kMC implementations.

Thin water films and particle morphology evolution in nanocrystalline MgO.

A key question in the field of ceramics and catalysis is how and to what extent residual water in the reactive environment of a metal oxide particle powder affects particle coarsening and morphology. With X-ray Diffraction (XRD) and Transmission Electron Microscopy (TEM), we investigated annealing-induced morphology changes on powders of MgO nanocubes in different gaseous H2O environments. The use of such a model system for particle powders enabled us to describe how adsorbed water that originates from short exposure to air determines the evolution of MgO grain size, morphology, and microstructure. While cubic nanoparticles with a predominant abundance of (100) surface planes retain their shape after annealing to T = 1173 K under continuous pumping with a base pressure of water p(H2O) = 10-5 mbar, higher water partial pressures promote mass transport on the surfaces and across interfaces of such particle systems. This leads to substantial growth and intergrowth of particles and simultaneously favors the formation of step edges and shallow protrusions on terraces. The mass transfer is promoted by thin films of water providing a two-dimensional solvent for Mg2+ ion hydration. In addition, we obtained direct evidence for hydroxylation-induced stabilization of (110) faces and step edges of the grain surfaces.

The Site-Assembly Determines Catalytic Activity of Nanoparticles.

Heterogeneous catalysts are often designed as metal nanoparticles supported on oxide surfaces. Here, the relation between particle morphology and reaction kinetics is investigated by scaling relation kinetic Monte Carlo simulations using CO oxidation over Pt nanoparticles as a model reaction. We find that different particle morphologies result in vastly different catalytic activities. The activity is strongly affected by kinetic couplings between sites, and a wide site distribution generally enhances the activity. The present study highlights the role of site-assemblies as a concept that, in addition to isolated active sites, can be used to understand catalytic reactions over nanoparticles.

2D-3D structural transition in sub-nanometer PtN clusters supported on CeO2(111).

Transition metal particles dispersed on oxide supports are used as heterogeneous catalysts in numerous applications. One example is platinum clusters supported on ceria which is used in automotive catalysis. Although control at the nm-scale is desirable to open new technological possibilities, there is limited knowledge both experimentally and theoretically regarding the geometrical structure and stability of sub-nanometer platinum clusters supported on ceria. Here we report a systematic, Density Functional Theory (DFT) study on the growth trends of CeO2(111) supported PtN clusters (N = 1-10). Using a global optimization methodology as a guidance tool to locate putative global minima, our results show a clear preference for 2D planar structures up to size Pt8. It is followed by a structural transition to 3D configurations at larger sizes. This remarkable trend is explained by the subtle competition between the formation of strong Pt-O bonds and the cluster internal Pt-Pt bonds. Our calculations show that the reducibility of CeO2(111) provides a mechanism to anchor PtN clusters where they become oxidized in a two-way charge transfer mechanism: (a) an oxidation process, where Osurface atoms withdraw charge from Pt atoms forming Pt-O bonds, (b) surface Ce4+ atoms are reduced, leading to Ce3+. The active role of the CeO2(111) support in modifying the structural and eventually the chemical properties of sub-nanometer PtN clusters is computationally demonstrated.

Tuning the Reactivity of Ultrathin Oxides: NO Adsorption on Monolayer FeO(111).

Ultrathin metal oxides exhibit unique chemical properties and show promise for applications in heterogeneous catalysis. Monolayer FeO films supported on metal surfaces show large differences in reactivity depending on the metal substrate, potentially enabling tuning of the catalytic properties of these materials. Nitric oxide (NO) adsorption is facile on silver-supported FeO, whereas a similar film grown on platinum is inert to NO under similar conditions. Ab initio calculations link this substrate-dependent behavior to steric hindrance caused by substrate-induced rumpling of the FeO surface, which is stronger for the platinum-supported film. Calculations show that the size of the activation barrier to adsorption caused by the rumpling is dictated by the strength of the metal-oxide interaction, offering a straightforward method for tailoring the adsorption properties of ultrathin films.

Methane Oxidation over PdO(101) Revealed by First-Principles Kinetic Modeling.

The catalytic oxidation of methane to carbon dioxide and water over PdO(101) is investigated with first-principles based microkinetic modeling. Extensive exploration of the reaction landscape allows for determination of preferred pathways at different reaction conditions. The predicted kinetic behavior is in good agreement with a range of experimental findings including reaction orders in methane, water, and oxygen as well as apparent activation energies. The results consolidate the role of the PdO(101) surface in the activity of PdO catalysts and offer starting points for computational design of materials with improved catalytic activity. Moreover, the study demonstrates the predictive power of first-principles based kinetic modeling for oxide surfaces when hybrid functionals are applied in conjugation with kinetic models that go beyond the mean-field approximation.