Search

The growing concern over the escalating levels of anthropogenic CO2 emissions necessitates effective strategies for its conversion to valuable chemicals and fuels. In this research, we embark on a comprehensive investigation of the nature of zirconia on a copper inverse catalyst under the conditions of CO2 hydrogenation to methanol. We employ density functional theory calculations in combination with the Grand Canonical Basin Hopping method, enabling an exploration of the free energy surface including a variable amount of adsorbates within the relevant reaction conditions. Our focus centers on a model three-atom Zr cluster on a Cu(111) surface decorated with various OH, O, and formate ligands, noted Zr3Ox (OH)y (HCOO)z/Cu(111), revealing major changes in the active site induced by various reaction parameters such as the gas pressure, temperature, conversion levels, and CO2/H2 feed ratios. Through our analysis, we have unveiled insights into the dynamic behavior of the catalyst. Specifically, under reaction conditions, we observe a large number of composition and structures with similar free energy for the catalyst, with respect to changing the type, number, and binding sites of adsorbates, suggesting that the active site should be regarded as a statistical ensemble of diverse structures that interconvert.

A combination of density functional theory (DFT) and experiments with atomically size-selected Ptn clusters deposited on indium-tin oxide (ITO) electrodes was used to examine the effects of applied potential and Ptn size on the electrocatalytic activity of Ptn (n = 1, 4, 7, and 8) for the hydrogen evolution reaction (HER). Activity is found to be negligible for isolated Pt atoms on ITO, increasing rapidly with Ptn size such that Pt7/ITO and Pt8/ITO have roughly double the activity per Pt atom compared to atoms in the surface layer of polycrystalline Pt. Both the DFT and experiment find that hydrogen under-potential deposition (Hupd) results in Ptn/ITO (n = 4, 7, and 8) adsorbing ∼2H atoms/Pt atom at the HER threshold potential, equal to ca. double the Hupd observed for Pt bulk or nanoparticles. The cluster catalysts under electrocatalytic conditions are hence best described as a Pt hydride compound, significantly departing from a metallic Pt cluster. The exception is Pt1/ITO, where H adsorption at the HER threshold potential is energetically unfavorable. The theory combines global optimization with grand canonical approaches for the influence of potential, uncovering the fact that several metastable structures contribute to the HER, changing with the applied potential. It is hence critical to include reactions of the ensemble of energetically accessible PtnHx/ITO structures to correctly predict the activity vs Ptn size and applied potential. For the small clusters, spillover of Hads from the clusters to the ITO support is significant, resulting in a competing channel for loss of Hads, particularly at slow potential scan rates.

The rearrangement of Cu surfaces under electrochemical conditions is known to play a key role in the surface activation for major electrocatalytic reactions. Despite the extensive experimental insights into such rearrangements, from surface-sensitive spectroscopy and microscopy, the spatial and temporal resolution of these methods is insufficient to provide an atomistic picture of the electrochemical interface. Theoretical characterization has also been challenged by the diversity of restructuring configurations, surface stoichiometry, adsorbate configurations, and the effect of the electrode potential. Here, atomistic insight into the restructuring of the electrochemical interface is gained from first principles. Cu(100) restructuring under varying applied potentials and adsorbate coverages is studied by grand canonical density functional theory and global optimization techniques, as well as ab initio molecular dynamics and mechanistic calculations. We show that electroreduction conditions cause the formation of a shifted-row reconstruction on Cu(100), induced by hydrogen adsorption. The reconstruction is initiated at 1/6 ML H coverage, when the Cu-H bonding sufficiently weakens the Cu-Cu bonds between the top- and sublayer, and further stabilized at 1/3 ML when H adsorbates fill all the created 3-fold hollow sites. The simulated scanning tunneling microscopy (STM) images of the calculated reconstructed interfaces agree with experimental in situ STM. However, compared to the thermodynamic prediction, the onsets of reconstruction events in the experiment occur at more negative applied voltages. This is attributed to kinetic effects in restructuring, which we describe via different statistical models, to produce the potential- and pH-dependent surface stability diagram. This manuscript provides rich atomistic insight into surface restructuring in electroreduction conditions, which is required for the understanding and design of Cu-based materials for electrocatalytic processes. It also offers the methodology to study the problem of in situ electrode reconstruction.

Dynamic restructuring of bimetallic catalysts plays a crucial role in their catalytic activity and selectivity. In particular, catalyst pretreatment with species such as carbon monoxide and oxygen has been shown to be an effective strategy for tuning the surface composition and morphology. Mechanistic and kinetic understanding of such restructuring is fundamental to the chemistry and engineering of surface active sites but has remained challenging due to the large structural, chemical, and temporal degrees of freedom. Here, we combine time-resolved temperature-programmed infrared reflection absorption spectroscopy, ab initio thermodynamics, and machine-learning molecular dynamics to uncover previously unidentified timescale and kinetic parameters of in situ restructuring in Pd/Au(111), a highly relevant model system for dilute Pd-in-Au nanoparticle catalysts. The key innovation lies in utilizing CO not only as a chemically sensitive probe of surface Pd but also as an agent that induces restructuring of the surface. Upon annealing in vacuum, as-deposited Pd islands became encapsulated by Au and partially dissolved into the subsurface, leaving behind isolated Pd monomers on the surface. Subsequent exposure to 0.1 mbar CO enabled Pd monomers to repopulate the surface up to 373 K, above which complete Pd dissolution occurred by 473 K, with apparent activation energies of 0.14 and 0.48 eV, respectively. These restructuring processes occurred over the span of ∼1000 s at a given temperature. Such a minute-timescale dynamics not only elucidates the fluxional nature of alloy catalysts but also presents an opportunity to fine-tune the surface under moderate temperature and pressure conditions.

Elucidation of reaction mechanisms and the geometric and electronic structure of the active sites themselves is a challenging, yet essential task in the design of new heterogeneous catalysts. Such investigations are best implemented via a multipronged approach that comprises ambient pressure catalysis, surface science, and theory. Herein, we employ this strategy to understand the workings of NiAu single-atom alloy (SAA) catalysts for the selective nonoxidative dehydrogenation of ethanol to acetaldehyde and hydrogen. The atomic dispersion of Ni is paramount for selective ethanol to acetaldehyde conversion, and we show that even the presence of small Ni ensembles in the Au surface results in the formation of undesirable byproducts via C-C scission. Spectroscopic, kinetic, and theoretical investigations of the reaction mechanism reveal that both C-H and O-H bond cleavage steps are kinetically relevant and single Ni atoms are confirmed as the active sites. X-ray absorption spectroscopy studies allow us to follow the charge of the Ni atoms in the Au host before, under, and after a reaction cycle. Specifically, in the pristine state the Ni atoms carry a partial positive charge that increases upon coordination to the electronegative oxygen in ethanol and decreases upon desorption. This type of oxidation state cycling during reaction is similar to the behavior of single-site homogeneous catalysts. Given the unique electronic structure of many single-site catalysts, such a combined approach in which the atomic-scale catalyst structure and charge state of the single atom dopant can be monitored as a function of its reactive environment is a key step toward developing structure-function relationships that inform the design of new catalysts.

The restructuring of interfaces plays a crucial role in materials science and heterogeneous catalysis. Bimetallic systems, in particular, often adopt very different compositions and morphologies at surfaces compared to the bulk. For the first time, we reveal a detailed atomistic picture of long-time scale restructuring of Pd deposited on Ag using microscopy, spectroscopy, and novel simulation methods. By developing and performing accelerated machine-learning molecular dynamics followed by an automated analysis method, we discover and characterize previously unidentified surface restructuring mechanisms in an unbiased fashion, including Pd-Ag place exchange and Ag pop-out as well as step ascent and descent. Remarkably, layer-by-layer dissolution of Pd into Ag is always preceded by an encapsulation of Pd islands by Ag, resulting in a significant migration of Ag out of the surface and a formation of extensive vacancy pits within a period of microseconds. These metastable structures are of vital catalytic importance, as Ag-encapsulated Pd remains much more accessible to reactants than bulk-dissolved Pd. Our approach is broadly applicable to complex multimetallic systems and enables the previously intractable mechanistic investigation of restructuring dynamics at atomic resolution.

The reaction of CO and O2 with submonolayer and multilayer CoOx films on Pt(111), to produce CO2, was investigated at room temperature in the mTorr pressure regime. Using operando ambient pressure X-ray photoelectron spectroscopy and high pressure scanning tunneling microscopy, as well as density functional theory calculations, we found that the presence of oxygen vacancies in partially oxidized CoOx films significantly enhances the CO oxidation activity to form CO2 upon exposure to mTorr pressures of CO at room temperature. In contrast, CoO films without O-vacancies are much less active for CO2 formation at RT, and CO only adsorbed in the form of carbonate species that are stable up to 260 °C. On submonolayer CoOx islands, the carbonates form preferentially at island edges, deactivating the edge sites for CO2 formation, even while the reaction proceeds inside the islands. These results provide a detailed understanding of CO oxidation pathways on systems where noble metals such as Pt interact with reducible oxides.

The reaction of CO and O2 with submonolayer and multilayer CoOx films on Pt(111), to produce CO2, was investigated at room temperature in the mTorr pressure regime. Using operando ambient pressure X-ray photoelectron spectroscopy and high pressure scanning tunneling microscopy, as well as density functional theory calculations, we found that the presence of oxygen vacancies in partially oxidized CoOx films significantly enhances the CO oxidation activity to form CO2 upon exposure to mTorr pressures of CO at room temperature. In contrast, CoO films without O-vacancies are much less active for CO2 formation at RT, and CO only adsorbed in the form of carbonate species that are stable up to 260 °C. On submonolayer CoOx islands, the carbonates form preferentially at island edges, deactivating the edge sites for CO2 formation, even while the reaction proceeds inside the islands. These results provide a detailed understanding of CO oxidation pathways on systems where noble metals such as Pt interact with reducible oxides.

Unlike nanostructured metal catalysts, supported single-atom catalysts (SACs) contain only atomically dispersed metal atoms, hinting at much more pronounced metal-support effects. Herein, we take a series of polyoxometalate-supported Pt catalysts as examples to quantitatively investigate the stability of Pt atoms on oxide supports and how the Pt-support interaction influences the catalytic performance. For this entire series, we show that the Pt atoms prefer to stay at a 4-fold hollow site of one polyoxometalate molecule and that the least adsorption energy to obtain sintering-resistant Pt SACs is 5.50 eV, which exactly matches the cohesive energy of bulk Pt metal. Further, we compared their catalytic performance in several hydrogenation reactions and simulated the reaction pathways of propene hydrogenation by density functional theory (DFT) calculations. Both experimental and theoretical approaches suggest that despite the Pt1-support interactions being different, the reaction pathways of various Pt1-polyoxometalate catalysts are very similar and their effective reaction barriers are close to each other and as low as 24 kJ/mol, indicating the possibility of obtaining SACs with improved stability without compromising activity. DFT calculations show that all reaction elementary steps take place only on the Pt atom without involving neighboring O atoms and that hydrogenation proceeds from the molecularly adsorbed H2 species. Pt SACs give a weaker H2 adsorption energy than Pt clusters or surfaces, resulting in small adsorption equilibrium constants and small apparent activation barriers, which agree between experiment and theory.

The interaction of water with metal oxides controls their activity and stability in heterogeneous catalysis and electrocatalysis. In this work, we combine density functional theory calculations and infrared reflection absorption spectroscopy (IRAS) to identify the structural motifs formed upon interaction of water with an atomically defined Co3O4(111) surface. Three principal structures are observed: (i) strongly bound isolated OD, (ii) extended hydrogen-bonded OD/D2O structures, and (iii) a third structure which has not been reported to our knowledge. In this structure, surface Co2+ ions bind to three D2O molecules to form an octahedrally coordinated Co2+ with a "half hydration shell". We propose that this hydration structure represents an important intermediate in reorganization and dissolution on oxide surfaces which expose highly unsaturated surface cations.

In alkene metathesis, while group 6 (Mo or W) high-oxidation state alkylidenes are accepted to be key reaction intermediates for both homogeneous and heterogeneous catalysts, it has been proposed that low valent species in their +4 oxidation state can serve as precatalysts. However, the activation mechanism for these latter species-generating alkylidenes-is still an open question. Here, we report the syntheses of tungsten(IV)-oxo bisalkoxide molecular complexes stabilized by pyridine ligands, WO(OR)2py3 (R = CMe(CF3)2 (2a), R = Si(O tBu)3 (2b), and R = C(CF3)3 (2c); py = pyridine), and show that upon activation with B(C6F5)3 they display alkene metathesis activities comparable to W(VI)-oxo alkylidenes. The initiation mechanism is examined by kinetic, isotope labeling and computational studies. Experimental evidence reveals that the presence of an allylic CH group in the alkene reactant is crucial for initiating alkene metathesis. Deuterium labeling of the allylic C-H group shows a primary kinetic isotope effect on the rate of initiation. DFT calculations support the formation of an allyl hydride intermediate via activation of the allylic C-H bond and show that formation of the metallacyclobutane from the allyl "hydride" involves a proton transfer facilitated by the coordination of a Lewis acid (B(C6F5)3) and assisted by a Lewis base (pyridine). This proton transfer step is rate determining and yields the metathesis active species.

Reactivity studies on catalytic transition metal clusters are usually performed on a single global minimum structure. With the example of a Pt13 cluster under a pressure of hydrogen, we show from first-principle calculations that low energy metastable structures of the cluster can play a major role for catalytic reactivity and that hence consideration of the global minimum structure alone can severely underestimate the activity. The catalyst is fluxional with an ensemble of metastable structures energetically accessible at reaction conditions. A modified genetic algorithm is proposed to comprehensively search for the low energy metastable ensemble (LEME) structures instead of merely the global minimum structure. In order to reduce the computational cost of density functional calculations, a high dimensional neural network potential is employed to accelerate the exploration. The presence and influence of LEME structures during catalysis is discussed by the example of H covered Pt13 clusters for two reactions of major importance: hydrogen evolution reaction and methane activation. The results demonstrate that although the number of accessible metastable structures is reduced under reaction condition for Pt13 clusters, these metastable structures can exhibit high activity and dominate the observed activity due to their unique electronic or structural properties. This underlines the necessity of thoroughly exploring the LEME structures in catalysis simulations. The approach enables one to systematically address the impact of isomers in catalysis studies, taking into account the high adsorbate coverage induced by reaction conditions.

A theoretical study of cohesion, structural deformation, inclusion and dynamics in porous hydrogen-bonded molecular networks was carried out. Density functional theory (DFT) calculations show that ability of these networks to undergo east deformations arises from the flexibility of the individual molecules and distortions of the directional interactions that hold them together.

The oxygen-promoted reconstruction of the Ag{111} surface is of direct significance to an important industrial catalytic process, the partial oxidation of ethene to ethene epoxide, C(sub 2)H(sub 4)O. A thorough electronic and structural analysis of the possible configurations has enabled us to understand the site selectivity on this surface.

A detailed study of hetergeneous catalysis through subsurface sites using a theoretical approach is presented.

The characteristics of the Rh-diamine precursor in the catalytic cycle of the asymmetric hydride transfer reduction of carbonyl compounds was theoretically and empirically investigated. Density functional equations were used on various (RhH(NH3)n(C2H4)4-n) complexes in the trigonal bipyramidal (TBP) and square pyramidal forms. Findings showed that the only stable complex corresponds to n = 2 in a TBP form, with the ethylenes in the equatorial plane. Meanwhile, mass analyses of the synthesized complexes indicated that their composition was (Rh(COD)(diamine))+ X- (X- = Cl- or PF6-).

An extended Huckel molecular orbital analysis is used as the framework for the comparison between the eta1 and eta2 coordinations of aldehydes and ketones on some types of organometallic fragments. The interaction between the hydrocarbon fragment and the organometallic center has two components, which include stabilizing two electron and destabilizing four electron interactions. The four electron interactions between low-lying orbitals determine the preferred structure. These interactions are stronger for the eta2 mode, making the structure less energetically favored.

The stability and reactivity of [epsilon], x, and [theta] iron carbide phases in Fischer-Tropsch synthesis (FTS) catalysts is examined as a function of relevant reaction conditions. The studies have disclosed the thermodynamic and kinetic factors inducing the [epsilon]-x-[theta] carbide transformation as a function of the carbon chemical potential [[mu].sub.C].

The periodic and cluster simulations are used to calculate the infrared OH stretching frequencies of the various types of hydroxyl groups on MgO surfaces. The IR frequency of OH groups is affected by the hydrogen bonding among the OH groups and location of OH groups.

A combination of experimental data and density functional theory calculations are used to study the nature and structure of alumina single site catalysts. The grafting mechanisms for these complexes on alumina are compared and are focused on the existence of cationic metallic species and on the unique character of the formed surface entity.

A study is conducted which invalidates the usual empirical rule stating that the C=C bond is more reactive than the C=O group for catalytic hydrogenations for unsaturated aldehydes. Density functional theory calculations of all the competitive hydrogenation routes of acrolein on Pt(111) reveal conversely that the attack at the C=O bond is systematically favored.