Structure and Chemical Reactivity of Y-Stabilized ZrO2 Surfaces: Importance for the Water-Gas Shift Reaction.

The surface structure and chemical properties of Y-stabilized zirconia (YSZ) have been subjects of intense debate over the past three decades. However, a thorough understanding of chemical processes occurring at YSZ powders faces significant challenges due to the absence of reliable reference data acquired for well-controlled model systems. Here, we present results from polarization-resolved infrared reflection absorption spectroscopy (IRRAS) obtained for differently oriented, Y-doped ZrO2 single-crystal surfaces after exposure to CO and D2O. The IRRAS data reveal that the polar YSZ(100) surface undergoes reconstruction, characterized by an unusual, red-shifted CO band at 2132 cm-1. Density functional theory calculations allowed to relate this unexpected observation to under-coordinated Zr4+ cations in the vicinity of doping-induced O vacancies. This reconstruction leads to a strongly increased chemical reactivity and water spontaneously dissociates on YSZ(100). The latter, which is an important requirement for catalysing the water-gas-shift (WGS) reaction, is absent for YSZ(111), where only associative adsorption was observed. Together with a novel analysis Scheme these reference data allowed for an operando characterisation of YSZ powders using DRIFTS (diffuse reflectance infrared Fourier transform spectroscopy). These findings facilitate rational design and tuning of YSZ-based powder materials for catalytic applications, in particular CO oxidation and the WGS reaction.

Anharmonic Correction to Free Energy Barriers from DFT-Based Molecular Dynamics Using Constrained Thermodynamic Integration.

For the calculation of anharmonic contributions to free energy barriers, constrained thermodynamic λ-path integration (λ-TI) from a harmonic reference force field to density functional theory is presented as an alternative to the established Blue Moon ensemble method (ξ-TI), in which free energy gradients along the reaction coordinate ξ are integrated. With good agreement in all cases, the λ-TI method is benchmarked against the ξ-TI method for several reactions, including the internal CH3 group rotation in ethane, a nucleophilic substitution of CH3Cl, a retro-Diels-Alder reaction, and a proton transfer in zeolite H-SSZ-13. An advantage of λ-TI is that one can use virtually any reference state to compute anharmonic contributions to reaction free energies or free energy barriers. This is particularly relevant for catalysis, where it is now possible to compute anharmonic corrections to the free energy of a transition state relative to any reference, for example, the most stable state of the active site and the reactants in the gas phase. This is in contrast to ξ-TI, where free energy barriers can only be computed relative to an initial state with all reactants coadsorbed. Finally, the Bennett acceptance ratio method combined with λ-TI is demonstrated to reduce the number of required integration grid points with tolerable accuracy, favoring thus λ-TI over ξ-TI in terms of computational efficiency.

Metal-Support Interactions in Heterogeneous Catalysis: DFT Calculations on the Interaction of Copper Nanoparticles with Magnesium Oxide.

Oxide supports play an important role in enhancing the catalytic properties of transition metal nanoparticles in heterogeneous catalysis. How extensively interactions between the oxide support and the nanoparticles impact the electronic structure as well as the surface properties of the nanoparticles is hence of high interest. In this study, the influence of a magnesium oxide support on the properties of copper nanoparticles with different size, shape, and adsorption sites is investigated using density functional theory (DFT) calculations. By proposing simple models to reduce the cost of the calculations while maintaining the accuracy of the results, we show using the nonreducible oxide support MgO as an example that there is no significant influence of the MgO support on the electronic structure of the copper nanoparticles, with the exception of adsorption directly at the Cu-MgO interface. We also propose a simplified methodology that allows us to reduce the cost of the calculations, while the accuracy of the results is maintained. We demonstrate in addition that the Cu nanowire model corresponds well to the nanoparticle model, which reduces the computational cost even further.

Rhodium Single-Atom Catalyst Design through Oxide Support Modulation for Selective Gas-Phase Ethylene Hydroformylation.

A frontier challenge in single-atom (SA) catalysis is the design of fully inorganic sites capable of emulating the high reaction selectivity traditionally exclusive of organometallic counterparts in homogeneous catalysis. Modulating the direct coordination environment in SA sites, via the exploitation of the oxide support's surface chemistry, stands as a powerful albeit underexplored strategy. We report that isolated Rh atoms stabilized on oxygen-defective SnO2 uniquely unite excellent TOF with essentially full selectivity in the gas-phase hydroformylation of ethylene, inhibiting the thermodynamically favored olefin hydrogenation. Density Functional Theory calculations and surface characterization suggest that substantial depletion of the catalyst surface in lattice oxygen, energetically facile on SnO2 , is key to unlock a high coordination pliability at the mononuclear Rh centers, leading to an exceptional performance which is on par with that of molecular catalysts in liquid media.

Templated encapsulation of platinum-based catalysts promotes high-temperature stability to 1,100 °C.

Stable catalysts are essential to address energy and environmental challenges, especially for applications in harsh environments (for example, high temperature, oxidizing atmosphere and steam). In such conditions, supported metal catalysts deactivate due to sintering-a process where initially small nanoparticles grow into larger ones with reduced active surface area-but strategies to stabilize them can lead to decreased performance. Here we report stable catalysts prepared through the encapsulation of platinum nanoparticles inside an alumina framework, which was formed by depositing an alumina precursor within a separately prepared porous organic framework impregnated with platinum nanoparticles. These catalysts do not sinter at 800 °C in the presence of oxygen and steam, conditions in which conventional catalysts sinter to a large extent, while showing similar reaction rates. Extending this approach to Pd-Pt bimetallic catalysts led to the small particle size being maintained at temperatures as high as 1,100 °C in air and 10% steam. This strategy can be broadly applied to other metal and metal oxides for applications where sintering is a major cause of material deactivation.

Competition between reverse water gas shift reaction and methanol synthesis from CO2: influence of copper particle size.

Converting CO2 into value-added chemicals and fuels, such as methanol, is a promising approach to limit the environmental impact of human activities. Conventional methanol synthesis catalysts have shown limited efficiency and poor stability in a CO2/H2 mixture. To design improved catalysts, crucial for the effective utilization of CO2, an in-depth understanding of the active sites and reaction mechanism is desired. The catalytic performance of a series of carbon-supported Cu catalysts, with Cu particle sizes in the range of 5 to 20 nm, was evaluated under industrially relevant temperature and pressure, i.e. 260 °C and 40 bar(g). The CO2 hydrogenation reaction exhibited clear particle size effects up to 13 nm particles, with small nanoparticles having the lower activity, but higher methanol selectivity. MeOH and CO formation showed a different size-dependence. The TOFCO increased from 1.9 × 10-3 s-1 to 9.4 × 10-3 s-1 with Cu size increasing from 5 nm to 20 nm, while the TOFMeOH was size-independent (8.4 × 10-4 s-1 on average). The apparent activation energies for MeOH and CO formation were size-independent with values of 63 ± 7 kJ mol-1 and 118 ± 6 kJ mol-1, respectively. Hence the size dependence was ascribed to a decrease in the fraction of active sites suitable for CO formation with decreasing particle size. Theoretical models and DFT calculations showed that the origin of the particle size effect is most likely related to the differences in formate coverage for different Cu facets whose abundancy depends on particle size. Hence, the CO2 hydrogenation reaction is intrinsically sensitive to the Cu particle size.

Mechanistic differences between methanol and dimethyl ether in zeolite-catalyzed hydrocarbon synthesis.

Water influences critically the kinetics of the autocatalytic conversion of methanol to hydrocarbons in acid zeolites. At very low conversions but otherwise typical reaction conditions, the initiation of the reaction is delayed in presence of H2O. In absence of hydrocarbons, the main reactions are the methanol and dimethyl ether (DME) interconversion and the formation of a C1 reactive mixture-which in turn initiates the formation of first hydrocarbons in the zeolite pores. We conclude that the dominant reactions for the formation of a reactive C1 pool at this stage involve hydrogen transfer from both MeOH and DME to surface methoxy groups, leading to methane and formaldehyde in a 1:1 stoichiometry. While formaldehyde reacts further to other C1 intermediates and initiates the formation of first C-C bonds, CH4 is not reacting. The hydride transfer to methoxy groups is the rate-determining step in the initiation of the conversion of methanol and DME to hydrocarbons. Thus, CH4 formation rates at very low conversions, i.e., in the initiation stage before autocatalysis starts, are used to gauge the formation rates of first hydrocarbons. Kinetics, in good agreement with theoretical calculations, show surprisingly that hydrogen transfer from DME to methoxy species is 10 times faster than hydrogen transfer from methanol. This difference in reactivity causes the observed faster formation of hydrocarbons in dry feeds, when the concentration of methanol is lower than in presence of water. Importantly, the kinetic analysis of CH4 formation rates provides a unique quantitative parameter to characterize the activity of catalysts in the methanol-to-hydrocarbon process.

Toward Computing Accurate Free Energies in Heterogeneous Catalysis: a Case Study for Adsorbed Isobutene in H-ZSM-5.

Herein, we propose a novel computational protocol that enables calculating free energies with improved accuracy by combining the best available techniques for enthalpy and entropy calculation. While the entropy is described by enhanced sampling molecular dynamics techniques, the energy is calculated using ab initio methods. We apply the method to assess the stability of isobutene adsorption intermediates in the zeolite H-SSZ-13, a prototypical problem that is computationally extremely challenging in terms of calculating enthalpy and entropy. We find that at typical operating conditions for zeolite catalysis (400 °C), the physisorbed π-complex, and not the tertiary carbenium ion as often reported, is the most stable intermediate. This method paves the way for sampling-based techniques to calculate the accurate free energies in a broad range of chemistry-related disciplines, thus presenting a big step forward toward predictive modeling.

Single alloy nanoparticle x-ray imaging during a catalytic reaction.

The imaging of active nanoparticles represents a milestone in decoding heterogeneous catalysts' dynamics. We report the facet-resolved, surface strain state of a single PtRh alloy nanoparticle on SrTiO3 determined by coherent x-ray diffraction imaging under catalytic reaction conditions. Density functional theory calculations allow us to correlate the facet surface strain state to its reaction environment­dependent chemical composition. We find that the initially Pt-terminated nanoparticle surface gets Rh-enriched under CO oxidation reaction conditions. The local composition is facet orientation dependent, and the Rh enrichment is nonreversible under subsequent CO reduction. Tracking facet-resolved strain and composition under operando conditions is crucial for a rational design of more efficient heterogeneous catalysts with tailored activity, selectivity, and lifetime.

A Mild One-Pot Reduction of Phosphine(V) Oxides Affording Phosphines(III) and Their Metal Catalysts.

The metal-free reduction of a range of phosphine(V) oxides employing oxalyl chloride as an activating agent and hexachlorodisilane as reducing reagent has been achieved under mild reaction conditions. The method was successfully applied to the reduction of industrial waste byproduct triphenylphosphine(V) oxide, closing the phosphorus cycle to cleanly regenerate triphenylphosphine(III). Mechanistic studies and quantum chemical calculations support the attack of the dissociated chloride anion of intermediated phosphonium salt at the silicon of the disilane as the rate-limiting step for deprotection. The exquisite purity of the resultant phosphine(III) ligands after the simple removal of volatiles under reduced pressure circumvents laborious purification prior to metalation and has permitted the facile formation of important transition metal catalysts.

Anharmonic Correction to Adsorption Free Energy from DFT-Based MD Using Thermodynamic Integration.

Adsorption processes are often governed by weak interactions for which the estimation of entropy contributions by means of the harmonic approximation is prone to be inaccurate. Thermodynamic integration (TI) from the harmonic to the fully interacting system (λ-path integration) can be used to compute anharmonic corrections. Here, we combine TI with (curvilinear) internal coordinates in periodic systems to make the formalism available in computational studies. Our implementation of ab initio molecular dynamics in VASP is independent of the reaction path and can be thus applied to study adsorption processes relative to the gas phase and does hence provide a useful tool for computational catalysis. We discuss the application of the approach on three model systems for which exact semianalytical solutions exist and illustrate and quantify the importance of anharmonic vibrations, hindered rotations, and hindered translations (dissociation). Eventually, we apply the method to study the adsorption of small adsorbates in a zeolite (H-SSZ-13).

Size-controlled nanocrystals reveal spatial dependence and severity of nanoparticle coalescence and Ostwald ripening in sintering phenomena.

A major aim in the synthesis of nanomaterials is the development of stable materials for high-temperature applications. Although the thermal coarsening of small and active nanocrystals into less active aggregates is universal in material deactivation, the atomic mechanisms governing nanocrystal growth remain elusive. By utilizing colloidally synthesized Pd/SiO2 powder nanocomposites with controlled nanocrystal sizes and spatial arrangements, we unravel the competing contributions of particle coalescence and atomic ripening processes in nanocrystal growth. Through the study of size-controlled nanocrystals, we can uniquely identify the presence of either nanocrystal dimers or smaller nanoclusters, which indicate the relative contributions of these two processes. By controlling and tracking the nanocrystal density, we demonstrate the spatial dependence of nanocrystal coalescence and the spatial independence of Ostwald (atomic) ripening. Overall, we prove that the most significant loss of the nanocrystal surface area is due to high-temperature atomic ripening. This observation is in quantitative agreement with changes in the nanocrystal density produced by simulations of atomic exchange. Using well-defined colloidal materials, we extend our analysis to explain the unusual high-temperature stability of Au/SiO2 materials up to 800 °C.

Metal-Specific Reactivity in Single-Atom Catalysts: CO Oxidation on 4d and 5d Transition Metals Atomically Dispersed on MgO.

Understanding and tuning the catalytic properties of metals atomically dispersed on oxides are major stepping-stones toward a rational development of single-atom catalysts (SACs). Beyond individual showcase studies, the design and synthesis of structurally regular series of SACs opens the door to systematic experimental investigations of performance as a function of metal identity. Herein, a series of single-atom catalysts based on various 4d (Ru, Rh, Pd) and 5d (Ir, Pt) transition metals has been synthesized on a common MgO carrier. Complementary experimental (X-ray absorption spectroscopy) and theoretical (Density Functional Theory) studies reveal that, regardless of the metal identity, metal cations occupy preferably octahedral coordination MgO lattice positions under step-edges, hence highly confined by the oxide support. Upon exposure to O2-lean CO oxidation conditions, FTIR spectroscopy indicates the partial deconfinement of the monatomic metal centers driven by CO at precatalysis temperatures, followed by the development of surface carbonate species under steady-state conditions. These findings are supported by DFT calculations, which show the driving force and final structure for the surface metal protrusion to be metal-dependent, but point to an equivalent octahedral-coordinated M4+ carbonate species as the resting state in all cases. Experimentally, apparent reaction activation energies in the range of 96 ± 19 kJ/mol are determined, with Pt leading to the lowest energy barrier. The results indicate that, for monatomic sites in SACs, differences in CO oxidation reactivity enforceable via metal selection are of lower magnitude than those evidenced previously through the mechanistic involvement of adjacent redox centers on the oxide carrier, suggesting that tuning of the oxide surface chemistry is as relevant as the selection of the supported metal.

The transformation of cuboctahedral to icosahedral nanoparticles: atomic structure and dynamics.

The rearrangement of transition metal nanoparticles from cuboctahedral to icosahedral structures is studied for up to 923 atoms. The atomic structure and temperature dependence of the transition are investigated with a well-defined collective variable. This collective variable describes the folding of the square fcc(100) facets into two triangular facets through a linear combination of the diagonals of all fcc(100) facets of all shells of the particle. Activation barriers are determined through harmonic transition state theory and constrained molecular dynamics simulations based on force field potentials. These calculations predict an activation entropy larger than 1 meV K-1, leading to strongly temperature dependent activation barriers. Density functional theory calculations were additionally performed both as single point calculations and as full optimizations. Cu, Ag, Au and Ni clusters show low barriers for concerted, symmetric transition up to the 309-atomic clusters. In contrast, for Pd, Pt, Rh and Ir higher barriers are required, already for the 147-atomic clusters. With increasing barriers, an asymmetric but still concerted rearrangement becomes energetically more favorable than the fully symmetric transformation. The material-dependence of the transition can be correlated with the melting point of the bulk metals.

How Accurately Do Approximate Density Functionals Predict Trends in Acidic Zeolite Catalysis?

Density functional theory (DFT) is increasingly used for computational screening procedures with the aim of finding new catalysts. To achieve this, it is critical that relative differences between materials are predicted with high accuracy. How DFT at the generalized gradient approximation (GGA) level performs in this respect is investigated here for catalytic reactions employing acidic zeotypes using highly accurate DLPNO-CCSD(T) calculations as the reference. This is studied for 65 reaction energies and 130 reaction barriers related to zeolite catalysis. Our results obtained for the PBE-D3 and BEEF-vdW functionals show that while these functionals are prone to large errors, they predict trends occurring from one catalyst to another with an accuracy of about 5 kJ/mol, strongly supporting the widespread use of DFT calculations for the computational screening and design of new catalytic materials.

Prospects of Heterogeneous Hydroformylation with Supported Single Atom Catalysts.

The potential of oxide-supported rhodium single atom catalysts (SACs) for heterogeneous hydroformylation was investigated both theoretically and experimentally. Using high-level domain-based local-pair natural orbital coupled cluster singles doubles with perturbative triples contribution (DLPNO-CCSD(T)) calculations, both stability and catalytic activity were investigated for Rh single atoms on different oxide surfaces. Atomically dispersed, supported Rh catalysts were synthesized on MgO and CeO2. While the CeO2-supported rhodium catalyst is found to be highly active, this is not the case for MgO, most likely due to increased confinement, as determined by extended X-ray absorption fine structure spectroscopy (EXAFS), that diminishes the reactivity of Rh complexes on MgO. This agrees well with our computational investigation, where we find that rhodium carbonyl hydride complexes on flat oxide surfaces such as CeO2(111) have catalytic activities comparable to those of molecular complexes. For a step edge on a MgO(301) surface, however, calculations show a significantly reduced catalytic activity. At the same time, calculations predict that stronger adsorption at the higher coordinated adsorption site leads to a more stable catalyst. Keeping the balance between stability and activity appears to be the main challenge for oxide supported Rh hydroformylation catalysts. In addition to the chemical bonding between rhodium complex and support, the confinement experienced by the active site plays an important role for the catalytic activity.

Vibrational Frequencies of Cerium-Oxide-Bound CO: A Challenge for Conventional DFT Methods.

In ceria-based catalysis, the shape of the catalyst particle, which determines the exposed crystal facets, profoundly affects its reactivity. The vibrational frequency of adsorbed carbon monoxide (CO) can be used as a sensitive probe to identify the exposed surface facets, provided reference data on well-defined single crystal surfaces together with a definitive theoretical assignment exist. We investigate the adsorption of CO on the CeO_{2}(110) and (111) surfaces and show that the commonly applied DFT(PBE)+U method does not provide reliable CO vibrational frequencies by comparing with state-of-the-art infrared spectroscopy experiments for monocrystalline CeO_{2} surfaces. Good agreement requires the hybrid DFT approach with the HSE06 functional. The failure of conventional density-functional theory (DFT) is explained in terms of its inability to accurately describe the facet- and configuration-specific donation and backdonation effects that control the changes in the C─O bond length upon CO adsorption and the CO force constant. Our findings thus provide a theoretical basis for the detailed interpretation of experiments and open up the path to characterize more complex scenarios, including oxygen vacancies and metal adatoms.

Identification of the Reaction Sequence of the MTO Initiation Mechanism Using Ab Initio-Based Kinetics.

The initiation of the methanol-to-olefins (MTO) process is investigated using a multiscale modeling approach where more than 100 ab initio computed (MP2:DFT) rate constants for H-SSZ-13 are used in a batch reactor model. The investigated reaction network includes the mechanism for initiation (42 steps) and a representative part of the autocatalytic olefin cycle (63 steps). The simulations unravel the dominant initiation pathway for H-SSZ-13: dehydrogenation of methanol to CO is followed by CO-methylation leading to the formation of the first C-C bond in methyl acetate despite high barriers of >200 kJ/mol. Our multiscale approach is able to shed light on the reaction sequence that ultimately leads to olefin formation and strikingly demonstrates that only with a full reactor model that includes autocatalysis with olefins as cocatalysts is one able to understand the initiation mechanism on the atomic scale. Importantly, the model also shows that autocatalysis takes over long before significant amounts of olefins are formed, thus guiding the interpretation of experimental results.

POCN Ni(ii) pincer complexes: synthesis, characterization and evaluation of catalytic hydrosilylation and hydroboration activities.

A series of iminophosphinite POCN pincer Ni(ii) complexes, (POCN)NiMe and (POCN)NiLn(BX4) (L = CH3CN, n = 0, 1; X = F, Ph, C6F5), have been developed and subjected to catalytic hydrosilylation of alkenes, aldehydes and ketones and hydroboration of carbonyl compounds. The stoichiometric reactivity of (POCN)NiMe and (POCN)Ni(BF4) with PhSiH3 and HBPin suggests that catalytic reactions proceed via the hydride intermediate (POCN)NiH. With regard to reactions with HBPin, efficient and mild hydroboration of a variety of carbonyl compounds, including highly chemoselective hydroboration of benzaldehyde in the presence of other common potent reductive functional groups, such as alkenes, alkynes, esters, amides, nitriles, nitro compounds and even ketones, and the first example of base metal catalyzed hydroboration of amides, including mild direct hydroborative reduction of primary and secondary amides to borylated amines were demonstrated for (POCN)NiMe.

Supported Catalyst Deactivation by Decomposition into Single Atoms Is Suppressed by Increasing Metal Loading.

In the high-temperature environments needed to perform catalytic processes, supported precious metal catalysts severely lose their activity over time. Even brief exposure to high temperatures can lead to significant losses in activity, which forces manufacturers to use large amounts of noble metals to ensure effective catalyst function for a required lifetime. Generally, loss of catalytic activity is attributed to nanoparticle sintering, or processes by which larger particles grow at the expense of smaller ones. Here, by independently controlling particle size and particle loading using colloidal nanocrystals, we reveal the opposite process as a novel deactivation mechanism: nanoparticles rapidly lose activity by high-temperature nanoparticle decomposition into inactive single atoms. This deactivation route is remarkably fast, leading to severe loss of activity in as little as ten minutes. Importantly, this deactivation pathway is strongly dependent on particle density and concentration of support defect sites. A quantitative statistical model explains how for certain reactions, higher particle densities can lead to more stable catalysts.