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.

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.

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.

Preparation and coherent manipulation of pure quantum states of a single molecular ion.

Laser cooling and trapping of atoms and atomic ions has led to advances including the observation of exotic phases of matter, the development of precision sensors and state-of-the-art atomic clocks. The same level of control in molecules could also lead to important developments such as controlled chemical reactions and sensitive probes of fundamental theories, but the vibrational and rotational degrees of freedom in molecules pose a challenge for controlling their quantum mechanical states. Here we use quantum-logic spectroscopy, which maps quantum information between two ion species, to prepare and non-destructively detect quantum mechanical states in molecular ions. We develop a general technique for optical pumping and preparation of the molecule into a pure initial state. This enables us to observe high-resolution spectra in a single ion (CaH+) and coherent phenomena such as Rabi flopping and Ramsey fringes. The protocol requires a single, far-off-resonant laser that is not specific to the molecule, so many other molecular ions, including polyatomic species, could be treated using the same methods in the same apparatus by changing the molecular source. Combined with the long interrogation times afforded by ion traps, a broad range of molecular ions could be studied with unprecedented control and precision. Our technique thus represents a critical step towards applications such as precision molecular spectroscopy, stringent tests of fundamental physics, quantum computing and precision control of molecular dynamics.

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.

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.

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.

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.

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.

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.

Dynamical Observation and Detailed Description of Catalysts under Strong Metal-Support Interaction.

Understanding the structures of catalysts under realistic conditions with atomic precision is crucial to design better materials for challenging transformations. Under reducing conditions, certain reducible supports migrate onto supported metallic particles and create strong metal-support states that drastically change the reactivity of the systems. The details of this process are still unclear and preclude its thorough exploitation. Here, we report an atomic description of a palladium/titania (Pd/TiO2) system by combining state-of-the-art in situ transmission electron microscopy and density functional theory (DFT) calculations with structurally defined materials, in which we visualize the formation of the overlayers at the atomic scale under atmospheric pressure and high temperature. We show that an amorphous reduced titania layer is formed at low temperatures, and that crystallization of the layer into either mono- or bilayer structures is dictated by the reaction environment and predicted by theory. Furthermore, it occurs in combination with a dramatic reshaping of the metallic surface facets.

Rendering Photoreactivity to Ceria: The Role of Defects.

The photoreactivity of ceria, a photochemically inert oxide with a large band gap, can be increased to competitive values by introducing defects. This previously unexplained phenomenon has been investigated by monitoring the UV-induced decomposition of N2 O on well-defined single crystals of ceria by using infrared reflection-absorption spectroscopy (IRRAS). The IRRAS data, in conjunction with theory, provide direct evidence that reducing the ceria(110) surface yields high photoreactivity. No such effects are seen on the (111) surface. The low-temperature photodecomposition of N2 O occurs at surface O vacancies on the (110) surface, where the electron-rich cerium cations with a significantly lowered coordination number cause a local lowering of the huge band gap (ca. 6 eV). The quantum efficiency of strongly reduced ceria(110) surfaces in the photodecomposition of N2 O amounts to 0.03 %, and is thus comparable to that reported for the photooxidation of CO on rutile TiO2 (110).

Alcohol amination with ammonia catalyzed by an acridine-based ruthenium pincer complex: a mechanistic study.

The mechanistic course of the amination of alcohols with ammonia catalyzed by a structurally modified congener of Milstein's well-defined acridine-based PNP-pincer Ru complex has been investigated both experimentally and by DFT calculations. Several key Ru intermediates have been isolated and characterized. The detailed analysis of a series of possible catalytic pathways (e.g., with and without metal-ligand cooperation, inner- and outer-sphere mechanisms) leads us to conclude that the most favorable pathway for this catalyst does not require metal-ligand cooperation.

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.

The first catalytic synthesis of an acrylate from CO2 and an alkene-a rational approach.

For more than three decades the catalytic synthesis of acrylates from the cheap and abundantly available C(1) building block carbon dioxide and alkenes has been an unsolved problem in catalysis research, both in academia and industry. Herein, we describe a homogeneous catalyst based on nickel that permits the catalytic synthesis of the industrially highly relevant acrylate sodium acrylate from CO(2), ethylene, and a base, as demonstrated, at this stage, by a turnover number of greater than 10 with respect to the metal.

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.

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.

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).

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.