Bimetallic catalysts provide opportunities to overcome scaling laws governing selectivity of CO2 reduction (CO2R). Cu/Au nanoparticles show promise for CO2R, but Au surface segregation on particles with sizes ≥7 nm prevent investigation of surface atom ensembles. Here we employ ultrasmall (2 nm) Cu/Au nanoparticles as catalysts for CO2R. The high surface to volume ratio of ultrasmall particles inhibits formation of a Au shell, enabling the study of ensemble effects in Cu/Au nanoparticles with controllable composition and uniform size and shape. Electrokinetics show a nonmonotonic dependence of C1 selectivity between CO and HCOOH, with the 3Au:1Cu composition showing the highest HCOOH selectivity. Density functional theory identifies Cu2/Au(211) ensembles as unique in their ability to synthesize HCOOH by stabilizing CHOO* while preventing H2 evolution, making C1 product selectivity a sensitive function of Cu/Au surface ensemble distribution, consistent with experimental findings. These results yield important insights into C1 branching pathways and demonstrate how ultrasmall nanoparticles can circumvent traditional scaling laws to improve the selectivity of CO2R.
We report scanning tunneling microscopy (STM) studies of individual adatoms deposited on an InSb(110) surface. The adatoms can be reproducibly dropped off from the STM tip by voltage pulses, and impact tunneling into the surface by up to â¼100×. The spatial extent and magnitude of the tunneling effect are widely tunable by imaging conditions such as bias voltage, set current and photoillumination. We attribute the effect to occupation of a (+/0) charge transition level, and switching of the associated adatom-induced band bending. The effect in STM topographic images is well reproduced by transport modeling of filling and emptying rates as a function of the tip position. STM atomic contrast and tunneling spectra are in good agreement with density functional theory calculations for In adatoms. The adatom ionization effect can extend to distances greater than 50 nm away, which we attribute to the low concentration and low binding energy of the residual donors in the undoped InSb crystal. These studies demonstrate how individual atoms can be used to sensitively control current flow in nanoscale devices.
Dopants have the potential to locally modify water-olivine interactions, which can impact geological processes, such as weathering, CO2 sequestration, and abiotic hydrocarbon generation. As a first step in understanding the role of dopants on the water structure and chemistry at water-olivine interfaces, water monomer adsorption on alkaline earth (AE) and transition metal (TM) doped forsterite(010) [Mg2SiO4(010)] surfaces was studied using density functional theory (DFT). Dopants that occur in olivine minerals were considered and consisted of Ca, Sr, and Ba for the AE dopants and Cr, Mn, Fe, Co, and Ni for the TM dopants. The water molecule adsorbs on the olivine surface through a metal-water bond (Me-Ow) and a hydrogen bond with an adjacent surface lattice oxygen (Ox-Hw). A frontier orbital analysis reveals that the 1b2, 3a1, and 1b1 (HOMO) of the water molecule are involved in the bonding. All of the TM dopants show strong net Me-Ow covalent bonding between 3a1 and 1b1 water orbitals and TM d states, while the AE dopants except for Mg2SiO4(010) show negligible Me-Ow covalent bonding. Both the AE and TM dopants show similar hydrogen bonding features involving both the 1b2 and 3a1 orbitals. While the AE cations show an overall lower Me-Ow covalent interaction, the AE dopants have strong electrostatic interactions between the positive metal cation and the negatively charged water dipole. A bonding model incorporating a linear combination of the covalent Me-Ow bond, the Ox-Hw hydrogen bond, the electrostatic interaction between the dopant cation and the H2O molecule, and the surface distortion energy is needed to capture the variation in the DFT adsorption energies on the olivine surfaces. The bonding analysis is able to identify the dominant contributions to water-dopant interactions and can serve as a basis for future studies of more realistic water-olivine interfaces.
We investigated the adsorption and oxidation of n-propane and cyclopropane (C3H8 and c-C3H6) on the IrO2(110) surface using temperature programmed reaction spectroscopy (TPRS) and density functional theory (DFT) calculations. We find that the activation of both C3H8 and c-C3H6 is facile on IrO2(110) at low temperature, and that the dissociated alkanes oxidize during TPRS to produce CO, CO2 and H2O above â¼400 K. Propane conversion to propylene is negligible during TPRS for the conditions studied. Our results show that the maximum yield of alkane that oxidizes during TPRS is higher for c-C3H6 compared with C3H8 (â¼0.30 vs. 0.18 monolayer) and that pre-hydrogenation of the surface suppresses c-C3H6 oxidation to a lesser extent than C3H8. Consistent with the experimental results, DFT predicts that C3H8 and c-C3H6 form σ-complexes on IrO2(110) and that C-H bond activation of the complexes as well as subsequent dehydrogenation are highly facile via H-transfer to Obr atoms (bridging O-atoms). Our calculations predict that propane conversion to gaseous propylene is kinetically disfavored on IrO2(110) because HObr recombination makes Obr atoms available to promote further dehydrogenation at lower temperatures than those needed for the adsorbed C3H6 intermediate to desorb as propylene. We also present evidence that that the ability for c-C3H6 to activate via ring-opening is responsible for cyclopropane attaining higher reaction yields during TPRS and exhibiting a weaker sensitivity to surface pre-hydrogenation compared with n-propane.
Realizing the efficient and selective conversion of ethane to ethylene is important for improving the utilization of hydrocarbon resources, yet remains a major challenge in catalysis. Herein, ethane dehydrogenation on the IrO2(110) surface is investigated using temperature-programmed reaction spectroscopy (TPRS) and density functional theory (DFT) calculations. The results show that ethane forms strongly bound σ-complexes on IrO2(110) and that a large fraction of the complexes undergo C-H bond cleavage during TPRS at temperatures below 200 K. Continued heating causes as much as 40% of the dissociated ethane to dehydrogenate and desorb as ethylene near 350 K, with the remainder oxidizing to COx species. Both TPRS and DFT show that ethylene desorption is the rate-controlling step in the conversion of ethane to ethylene on IrO2(110) during TPRS. Partial hydrogenation of the IrO2(110) surface is found to enhance ethylene production from ethane while suppressing oxidation to COx species. DFT predicts that hydrogenation of reactive oxygen atoms of the IrO2(110) surface effectively deactivates these sites as H atom acceptors, and causes ethylene desorption to become favored over further dehydrogenation and oxidation of ethane-derived species. The study reveals that IrO2(110) exhibits an exceptional ability to promote ethane dehydrogenation to ethylene near room temperature, and provides molecular-level insights for understanding how surface properties influence selectivity toward ethylene production.
Methane undergoes highly facile C-H bond cleavage on the stoichiometric IrO2(110) surface. From temperature-programmed reaction spectroscopy experiments, we found that methane molecularly adsorbed as a strongly bound σ complex on IrO2(110) and that a large fraction of the adsorbed complexes underwent C-H bond cleavage at temperatures as low as 150 kelvin (K). The initial dissociation probability of methane on IrO2(110) decreased from 80 to 20% with increasing surface temperature from 175 to 300 K. We estimate that the activation energy for methane C-H bond cleavage is 9.5 kilojoule per mole (kJ/mol) lower than the binding energy of the adsorbed precursor on IrO2(110), and equal to a value of ~28.5 kJ/mol. Low-temperature activation may avoid unwanted side reactions in the development of catalytic processes to selectively convert methane to value-added products.
We investigated the molecular adsorption of methane, ethane, propane and n-butane on stoichiometric and oxygen-rich RuO2(110) surfaces using temperature-programmed desorption (TPD) and dispersion-corrected density functional theory (DFT-D3) calculations. We find that each alkane adsorbs strongly on the coordinatively-unsaturated Ru (Rucus) atoms of s-RuO2(110), with desorption from this state producing a well-defined TPD peak at low alkane coverage. As the coverage increases, we find that alkanes first form a compressed layer on the Rucus atoms and subsequently adsorb on the bridging O atoms of the surface until the monolayer saturates. DFT-D3 calculations predict that methane preferentially adsorbs on top of a Rucus atom and that the C2 to C4 alkanes preferentially adopt bidentate configurations in which each molecule aligns parallel to the Rucus atom row and datively bonds to neighboring Rucus atoms. DFT-D3 predicts binding energies that agree quantitatively with our experimental estimates for alkane σ-complexes on RuO2(110). We find that oxygen atoms adsorbed on top of Rucus atoms (Oot atoms) stabilize the adsorbed alkane complexes that bind in a given configuration, while also blocking the sites needed for σ-complex formation. This site blocking causes the coverage of the most stable, bidentate alkane complexes to decrease sharply with increasing Oot coverage. Concurrently, we find that a new peak develops in the C2 to C4 alkane TPD spectra with increasing Oot coverage, and that the desorption yield in this TPD feature passes through a maximum at Oot coverages between â¼50% and 60%. We present evidence that the new TPD peak arises from C2 to C4 alkanes that adsorb in upright, monodentate configurations on stranded Rucus sites located within the Oot layer.
Achieving selective C-H bond cleavage is critical for developing catalytic processes that transform small alkanes to value-added products. The present study clarifies the molecular-level origin for an exceptionally strong preference for propane to dissociate on the crystalline PdO(101) surface via primary C-H bond cleavage. Using reflection absorption infrared spectroscopy (RAIRS) and density functional theory (DFT) calculations, we show that adsorbed propane σ-complexes preferentially adopt geometries on PdO(101) in which only primary C-H bonds datively interact with the surface Pd atoms at low propane coverages and are thus activated under typical catalytic reaction conditions. We show that a propane molecule achieves maximum stability on PdO(101) by adopting a bidentate geometry in which a H-Pd dative bond forms at each CH3 group. These results demonstrate that structural registry between the molecule and surface can strongly influence the selectivity of a metal oxide surface in activating alkane C-H bonds.
Proton reduction is one of the most fundamental and important reactions in nature. MoS2 edges have been identified as the active sites for hydrogen evolution reaction (HER) electrocatalysis. Designing molecular mimics of MoS2 edge sites is an attractive strategy to understand the underlying catalytic mechanism of different edge sites and improve their activities. Herein we report a dimeric molecular analogue [Mo2 S12 ](2-) , as the smallest unit possessing both the terminal and bridging disulfide ligands. Our electrochemical tests show that [Mo2 S12 ](2-) is a superior heterogeneous HER catalyst under acidic conditions. Computations suggest that the bridging disulfide ligand of [Mo2 S12 ](2-) exhibits a hydrogen adsorption free energy near zero (-0.05â eV). This work helps shed light on the rational design of HER catalysts and biomimetics of hydrogen-evolving enzymes.
Metal oxide films can form on late transition-metal catalysts under sufficiently oxygen-rich conditions, and typically cause significant changes in the catalytic performance of these materials. Several investigations using sensitive in situ surface characterization techniques reveal that the CO oxidation activity of Pd and other late transition-metal catalysts increases abruptly under conditions at which metal oxide structures begin to develop. Findings such as these provide strong motivation for developing atomic-scale descriptions of oxidation catalysis over oxide films of the late transition-metals. Surface oxygen vacancies can play a central role in mediating oxidation catalysis promoted by metal oxides. In general, adsorbed reactants abstract oxygen atoms from the lattice of the oxide surface, thereby creating oxygen vacancies, while gaseous O2 replenishes the reactive surface oxygen atoms and eliminates oxygen vacancies. Oxygen vacancies also represent a distinct type of surface site on which the binding and reactivity of adsorbed species can differ compared with sites on the pristine oxide surface. Detailed characterization of vacancy-mediated rate processes is thus essential for developing reliable mechanistic descriptions of oxidation catalysis over reducible metal oxide films. Careful measurements performed in ultrahigh vacuum (UHV) using well-defined oxide surfaces in combination with molecular simulations afford the capability to isolate and characterize such reaction steps, and thus provide information that is needed for developing mechanistic models of oxidation catalysis over metal oxides. In this Account, we discuss vacancy-mediated processes that are involved in the oxidation of CO on the PdO(101) surface as determined from UHV surface science experiments and density functional theory (DFT) calculations. These studies show that CO binds strongly on Pd atoms that are located next to surface oxygen vacancies, and diffuses rapidly to these sites during reduction of the oxide surface by CO. The enhanced binding also raises the energy barriers for desorption and oxidation of CO, but the difference in these barriers remains nearly identical to that for CO adsorbed on the pristine PdO(101) surface. These recent studies also show that oxygen from the subsurface efficiently eliminates surface oxygen vacancies during CO oxidation at temperatures as low as 400 K, and thereby reveal a facile pathway by which PdO(101) surface domains can be maintained during oxide reduction.
Density functional theory was used to study CO adsorption on thin Pt metal films supported on SrO- and TiO2-terminated SrTiO3(100) surfaces. Regardless of substrate-termination, significant enhancement in CO binding occurred on the Pt monolayer compared to the bulk Pt(100) surface. We also observed CO-coverage dependent shifting of Pt atoms, influenced by the nature of underlying oxide atoms. These oxide-induced effects become negligible after depositing more than 2 monolayers of Pt. Evaluating the electronic structures of oxide-supported Pt showed that the interaction of filled Pt dxz+yz and empty Pt dz(2) states with CO molecular orbitals can be directly related to CO adsorption on the Pt/SrTiO3(100) surface. A hybrid d-band model is able to capture the CO adsorption trends for systems that do not show large lateral distortion except for the case of Pt adsorbed above the Sr atom on the SrO-termination. For this case, charge transfer from adjacent Pt atoms leads to a large filled dz(2) peak below the Fermi level that weakens the Pt-CO σ bonding due to Pauli repulsion.
Copper cathodes, at sufficiently negative potentials, are selective for hydrocarbon production during the electrochemical reduction of carbon dioxide. Other metals, such as Pt, Fe, Ni and Co, produce low to zero hydrocarbons. We employ density functional theory to examine the coverage of reaction intermediates under CO2 electroreduction conditions. A detailed thermodynamic analysis suggests that a high coverage of adsorbed CO at relevant reduction potentials blocks the metal surface sites for H adsorption, preventing C-H bond formation. The potential-dependent energetics of H adsorption and CO formation are highly sensitive to the surface coverage of the adsorbed species. The formation of surface carbon as a competing adsorption intermediate is also explored at relevant reduction potentials. CO2 electroreduction to hydrocarbons over metals active for the thermal reduction process (Fe, Ni, Co, Pt) would require a H supply for C-H bond formation that is competitive with CO* and C* at the surface.
Advances in the fundamental understanding of alkane activation on oxide surfaces are essential for developing new catalysts that efficiently and selectively promote chemical transformations of alkanes. In this tutorial review, we discuss the current understanding of alkane activation on crystalline metal oxide surfaces, and focus mainly on summarizing our findings on alkane adsorption and C-H bond cleavage on the PdO(101) surface as determined from model ultrahigh vacuum experiments and theoretical calculations. These studies show that alkanes form strongly-bound σ-complexes on PdO(101) by datively bonding with coordinatively-unsaturated Pd atoms and that these molecularly adsorbed species serve as precursors for C-H bond activation on the oxide surface. In addition to discussing the binding and properties of alkane σ-complexes on PdO(101), we also summarize recent advances in kinetic models to predict alkane dissociation rates on solid surfaces. Lastly, we highlight computations which predict that the formation and facile C-H bond activation of alkane σ-complexes also occurs on RuO2 and IrO2 surfaces.
We used conventional density functional theory (DFT) and dispersion-corrected DFT (DFT-D3) calculations to investigate C-H bond activation pathways for methane and ethane σ-complexes adsorbed on the PdO(101) surface. The DFT-D3 calculations predict lower and more physically realistic values of the apparent C-H bond cleavage barriers, which are defined relative to the gas-phase energy level, while giving nearly the same energy differences between stationary states as predicted by conventional DFT for a given reaction pathway. For the stable CH4 η(2) complex on PdO(101), DFT-D3 predicts that the C-H bond cleavage barriers are 55.2 and 16.1 kJ∕mol relative to the initial molecularly adsorbed and gaseous states, respectively. We also predict that dehydrogenation of the resulting CH3 groups and conversion to CH3O species are significantly more energetically demanding than the initial C-H bond activation of CH4 on PdO(101). Using DFT-D3, we find that an η(2) and an η(1) ethane complex can undergo C-H bond cleavage on PdO(101) with intrinsic energy barriers that are similar to that of the methane complex, but with apparent barriers that are close to zero. We also investigated the dissociation kinetics of methane and ethane on PdO(101) using microkinetic models, with parameters derived from the DFT-D3 relaxed structures. We find that a so-called 3N - 2 model, in which two frustrated adsorbate motions are treated as free motions, predicts desorption pre-factors and alkane dissociation probabilities that agree well with estimates obtained from the literature. The microkinetic simulations demonstrate the importance of accurately describing entropic contributions in kinetic simulations of alkane dissociative chemisorption.
We used temperature programmed reaction spectroscopy (TPRS) and molecular beam reflectivity measurements to investigate the initial dissociation of n-butane isotopologues on PdO(101) and determine kinetic parameters governing the selectivity of initial C-H(D) bond cleavage. We observe differences in the reactivity of the n-butane isotopologues on PdO(101) due to kinetic isotope effects, and find that the initial dissociation probability decreases with increasing surface temperature for each isotopologue. We performed an analysis of the dissociation kinetics using a model that is based on a precursor-mediated mechanism for n-butane dissociation and enables quantification of kinetic parameters for selective C-H bond cleavage by considering differences in the reactivity among the n-butane isotopologues. From the analysis, we estimate that 49% of the n-butane molecules which react during TPRS do so through 1° C-H bond cleavage when the initial coverage of n-butane lies between â¼40% and 100% of the saturation coverage of the molecular precursor state. For dissociation in the limit of zero coverage, we estimate that the conditional probability for 1° C-H bond cleavage is equal to â¼87% and varies only weakly with surface temperature from 300 K to 400 K. Analysis of the temperature dependent rate data further predicts that the barrier for 1° C-H bond cleavage is 3.5 kJ mol(-1) lower than that for 2° C-H bond cleavage for n-butane dissociation on PdO(101) in the limit of zero coverage. Our results provide evidence that the selectivity for 1° C-H bond cleavage on PdO(101) increases as the n-butane coverage decreases below â¼40% of the saturation value. We speculate that intermolecular interactions among the n-butane species are responsible for the apparent coverage dependence of the C-H bond selectivity for n-butane dissociation on PdO(101).
On the right path: Based on DFT calculations (incorporating the role of water solvation) of the activation barriers of elementary steps, a new path that leads to methane and ethylene for CO(2) electroreduction on Cu(111) was identified. Methane formation proceeds through reduction of CO to COH (path II, see picture), which leads to CH(x) species that can produce both methane and ethylene, as observed experimentally.
We used dispersion-corrected density functional theory (DFT-D3) calculations to investigate the initial C-H bond cleavage of propane σ-complexes adsorbed on the PdO(101) surface. The calculations predict that propane molecules adsorbed in η(1) configurations can undergo facile C-H bond cleavage on PdO(101), where the energy barrier for C-H bond activation is lower than that for desorption for each molecular complex. The preferred pathway for propane dissociation on PdO(101) corresponds to cleavage of a primary C-H bond of a so-called staggered p-2η(1) complex which initially coordinates with the surface by forming two H-Pd dative bonds, one at each CH(3) group. Among all of the adsorbed propane complexes, the staggered p-2η(1) complex has the highest binding energy and must overcome the lowest energy barrier for C-H bond scission. Analysis of the atomic charges reveals that propane C-H bond cleavage occurs heterolytically on PdO(101), and suggests that primary C-H bond activation is favored because a more stabilizing charge distribution develops within the 1-propyl transition state structures. Lastly, we conducted kinetic simulations using microkinetic models derived from the DFT-D3 structures, and find that the models reproduce the apparent activation energy for propane dissociation on PdO(101) to within 14% of that determined experimentally. We show that the entropic contributions of the adsorbed transition structures greatly exceed those predicted by the harmonic oscillator model, and that quantitative agreement with the apparent dissociation pre-factor may be obtained by approximating two of the frustrated adsorbate motions as free motions while treating the remaining modes as harmonic vibrations.
We investigated the molecular binding of n-alkanes on Pd(111) and PdO(101) using conventional density functional theory (DFT) and the dispersion-corrected DFT-D3 method. In agreement with experimental findings, DFT-D3 predicts that the n-alkane desorption energies scale linearly with the molecule chain length on both surfaces, and that n-alkanes bind more strongly on PdO(101) than on Pd(111). The desorption energies computed using DFT-D3 are slightly higher than the measured values for n-alkanes on Pd(111), though the agreement between computation and experiment is a significant improvement over conventional DFT. The measured desorption energies of n-alkanes on PdO(101) and the energies computed using DFT-D3 agree to within better than 2.5 kJ/mol (< 5%) for chain lengths up to n-butane. The DFT-D3 calculations predict that the molecule-surface dispersion energy for a given n-alkane is similar in magnitude on Pd(111) and PdO(101), and that dative bonding between the alkanes and coordinatively unsaturated Pd atoms is primarily responsible for the enhanced binding of n-alkanes on PdO(101). From analysis of the DFT-D3 results, we estimate that the strength of an alkane η(2)(H, H) interaction on PdO(101) is ~16 kJ/mol, while a single η(1) H-Pd dative bond is worth about 10 kJ/mol.
We investigated regioselectivity in the initial C-H bond activation of propane σ-complexes on the PdO(101) surface using temperature programmed reaction spectroscopy (TPRS) experiments. We observe a significant kinetic isotope effect (KIE) in the initial C-H(D) bond cleavage of propane on PdO(101) such that the dissociation yield of C(3)H(8) is 2.7 times higher than that of C(3)D(8) at temperatures between 150 and 200 K. Measurements of the reactivity of (CH(3))(2)CD(2) and (CD(3))(2)CH(2) show that deuteration of the methyl groups is primarily responsible for the lower reactivity of C(3)D(8) relative to C(3)H(8), and thus that 1° C-H bond cleavage is the preferred pathway for propane activation on PdO(101). By analyzing the rate data within the context of a kinetic model for precursor-mediated dissociation, we estimate that 90% of the propane σ-complexes which dissociate on PdO(101) during TPRS do so by 1° C-H bond cleavage.
Density functional theory at the level of the local density approximation with the projector augmented wave method is used to determine the structure of 180° domain walls in tetragonal ferroelectric PbTiO(3). In agreement with previous studies, it is found that PbO-centered {100} walls have lower energies than TiO(2)-centered {100} walls, leading to a Peierls potential barrier for wall motion along <010> of â¼36 mJ m(-2). In addition to the Ising-like polarization along the tetragonal axis, it is found that near the domain wall, there is a small polarization in the wall-normal direction away from the domain wall. These Néel-like contributions to the domain wall are analyzed in terms of the Landau-Ginzburg-Devonshire phenomenological theory for ferroelectrics. Similar characteristics are found for {110} domain walls, where OO-centered walls are energetically more favorable than the PbTiO-centered walls.
