Bifunctional hydroformylation on heterogeneous Rh-WOx pair site catalysts.

Metal-catalysed reactions are often hypothesized to proceed on bifunctional active sites, whereby colocalized reactive species facilitate distinct elementary steps in a catalytic cycle1-8. Bifunctional active sites have been established on homogeneous binuclear organometallic catalysts9-11. Empirical evidence exists for bifunctional active sites on supported metal catalysts, for example, at metal-oxide support interfaces2,6,7,12. However, elucidating bifunctional reaction mechanisms on supported metal catalysts is challenging due to the distribution of potential active-site structures, their dynamic reconstruction and required non-mean-field kinetic descriptions7,12,13. We overcome these limitations by synthesizing supported, atomically dispersed rhodium-tungsten oxide (Rh-WOx) pair site catalysts. The relative simplicity of the pair site structure and sufficient description by mean-field modelling enable correlation of the experimental kinetics with first principles-based microkinetic simulations. The Rh-WOx pair sites catalyse ethylene hydroformylation through a bifunctional mechanism involving Rh-assisted WOx reduction, transfer of ethylene from WOx to Rh and H2 dissociation at the Rh-WOx interface. The pair sites exhibited >95% selectivity at a product formation rate of 0.1 gpropanal cm-3 h-1 in gas-phase ethylene hydroformylation. Our results demonstrate that oxide-supported pair sites can enable bifunctional reaction mechanisms with high activity and selectivity for reactions that are performed in industry using homogeneous catalysts.

High Yield of L-Sorbose via D-Glucose Isomerization in Ethanol over a Bifunctional Titanium-Boron-Beta Zeolite.

D-Glucose-to-L-sorbose isomerization on Lewis acidic zeolite is a highly attractive avenue for sorbose production. But the L-sorbose yield is limited by the competing D-glucose-to-D-fructose isomerization and reaction equilibrium. In this work, it is suggested that ethanol directs the glucose conformation for selective D-glucose-to-L-sorbose isomerization. It also reacts with the produced L-sorbose to form ethyl-sorboside, which allows the D-glucose-to-L-sorbose isomerization to proceed beyond the thermodynamic equilibrium limit.  It is shown that a bifunctional zeolite Beta containing framework titanium (Ti) and boron (B) is a selective catalyst for this tandem reaction: Lewis acidic framework Ti catalyzes the D-glucose-to-L-sorbose isomerization via an intramolecular 5,1-hydride shift process as confirmed by isotopic tracing experiments followed by 13C-NMR, while weak Brønsted acid framework B selectively promotes the sorbose ketalization with ethanol. A remarkably high yield of L-sorbose with a high fraction of sugar (>95%: 27% unreacted glucose, 11.4% fructose, 57% sorbose) was obtained after the mixture produced in ethanol was hydrolyzed.

Tuning Active Site Flexibility by Defect Engineering of Graphene Ribbon Edge-hosted Fe-N3 Sites.

Nitrogen-doped, carbon-supported transition metal catalysts are excellent for several reactions. Structural engineering of M-Nx sites to boost catalytic activity is rarely studied. Here, we demonstrate that the structural flexibility of Fe-N3 site is vital for tuning the electronic structure of Fe atoms and regulating the catalytic transfer hydrogenation (CTH) activity. By introducing carbon defects, we construct Fe-N3 sites with varying Fe-N bond lengths distinguishable by X-ray absorption spectroscopy. We investigate the CTH activity by density-functional theory and microkinetic calculations and reveal that the vertical displacement of the Fe atom out of the plane of the support, induced by the Fe-N3 distortion, raises the Fe 3 d z 2 ${3{d}_{{z}^{2}}{\rm \ }}$ orbital and strengthens binding. We propose that the activity is controlled by the relaxation of the reconstructed site, which is further affected by Fe-N bond length, an excellent activity descriptor. We elucidate the origin of the CTH activity and principles for high-performing Fe-N-C catalysts by defect engineering.

Propane Dehydrogenation on PtxZny Active Sites in Silicalite-1.

The improvement of Pt-based catalysts for propane dehydrogenation (PDH) has progressed by recent investigations that have identified Zn as a promising promoter for Pt subnanometer catalysts. It is desirable to gain insights into the structure, stability, and activity of such active sites and the factors that influence them, such as Zn:Pt ratio, Pt coordination and nuclearity. Here, we employ density functional theory and microkinetic simulations to investigate the stability of PtxZny (x=1-3, y=0-3) active sites grafted on silanols of Silicalite-1 and the PDH activity of Pt. We find that the coordination of a Pt atom to a nest of grafted Zn(II) atoms increases the stability of the Pt1Zny sites, whose activity is similar for y=0-2 and drops dramatically for y>2. We further demonstrate, via linear scaling relations and microkinetic simulations, that the turnover frequency obeys a volcano law as a function of propylene binding strength. The Pt2Zn1 and Pt3Zn1 sites are stable and exhibit activity similar to Pt1Zn2, but only Pt1Zn2 manifests reaction kinetics consistent with experimental data, strongly suggesting the active site composition in the synthesized catalyst samples. The methodology presented here suggests a general strategy for deducing active site information such as composition through simple kinetic experiments.

Simultaneously Enhanced Hydrophilicity and Stability of a Metal-Organic Framework via Post-Synthetic Modification for Water Vapor Sorption/Desorption.

With increasing demands for high-performance water sorption materials, metal-organic frameworks (MOFs) have gained considerable attention due to their high maximum uptake capacities. In many cases, however, high overall capacity is not necessarily accomplishing high working capacity under operating conditions, due to insufficient hydrophilicity and/or water stability. Herein, we present a post-synthetic modification (PSM) of MOF-808, with di-sulfonic acids enhancing simultaneously its hydrophilicity and water stability without sacrificing its uptake capacity of ≈30 mmol g-1 . Di-sulfonic acid PSM enabled a shift of the relative humidity (RH) associated with a sharp step in water vapor sorption from 35-40 % RH in MOF-808 to below 25 % RH. While MOF-808 lost uptake capacity and crystallinity over multiple sorption/desorption cycles, the di-sulfonic acid PSM MOF-808 retained >80 % of the original capacity. PSM MOF-808 exhibited good hydrothermal stability up to 60 °C and high swing capacity.

P-Site Structural Diversity and Evolution in a Zeosil Catalyst.

Phosphorus-modified siliceous zeolites, or P-zeosils, catalyze the selective dehydration of biomass derivatives to platform chemicals such as p-xylene and 1,3-butadiene. Water generated during these reactions is a critical factor in catalytic activity, but the effects of hydrolysis on the structure, acidity, and distribution of the active sites are largely unknown. In this study, the P-sites in an all-silica self-pillared pentasil (P-SPP) with a low P-loading (Si/P = 27) were identified by solid-state 31P NMR using frequency-selective detection. This technique resolves overlapping signals for P-sites that are covalently bound to the solid phase, as well as oligomers confined in the zeolite but not attached to the zeolite. Dynamic Nuclear Polarization provides the sensitivity necessary to conduct 29Si-filtered 31P detection and 31P-31P correlation experiments. The aforementioned techniques allow us to distinguish sites with P-O-Si linkages from those with P-O-P linkages. The spectra reveal a previously unappreciated diversity of P-sites, including evidence for surface-bound oligomers. In the dry P-zeosil, essentially all P-sites are anchored to the solid phase, including mononuclear sites and dinuclear sites containing the [Si-O-P-O-P-O-Si] motif. The fully-condensed sites evolve rapidly when exposed to humidity, even at room temperature. Partially hydrolyzed species have a wide range of acidities, inferred from their calculated LUMO energies. Initial cleavage of some P-O-Si linkages results in an evolving mixture of surface-bound mono- and oligonuclear P-sites with increased acidity. Subsequent P-O-P cleavage leads to a decrease in acidity as the P-sites are eventually converted to H3PO4. The ability to identify acidic sites in P-zeosils and to describe their structure and stability will play an important role in controlling the activity of microporous catalysts by regulating their water content.

Stabilizing Porosity in Organic Cages through Coordination Chemistry.

The number of studies concerning the permanent porosity of molecular materials, especially porous organic cages (POCs) and porous coordination cages (PCCs), have increased substantially over the past decade. The work presented here outlines novel approaches to the preparation of porous molecular structures upon metalation of nonporous, amine-based organic cages. Reduction of the well-known CC3 and CC1 imine-based POCs affords nonporous, highly flexible amine cages. These materials can be endowed with significant levels of structural rigidity via post-synthetic metalation of their ethylenediamine-type binding pockets. The hybrid metal-organic cages accessed through this approach combine aspects of POC and PCC chemistry, with structures of this type providing a potentially promising new direction for the design and development of porous molecular materials with tunability in overall charge, metal cation, porosity, and solubility.

Carbocation-Mediated Cyclization of Trienes in Acid Zeolites.

The mechanism by which acid zeolites catalyze the formation of aromatic species is not fully understood and is important in an array of industrial processes such as the methanol to gasoline reaction. The so-called "carbon pool" mechanism is generally agreed to be the main channel for the formation of hydrocarbons from methanol. There is, however, no agreed sequence of elementary steps that explains how linear intermediates transform to cyclic intermediates, let alone aromatic rings. Recent work suggests the formation of conjugated trienes during zeolite-catalyzed aromatization, but mechanisms involving triene-derived carbocations have never been investigated using modern computational tools. In this work, we propose a new mechanism for cyclization of hexatriene over the Brønsted acid site of faujasite zeolite. Microkinetic models (MKM) using the results of Density Functional Theory (DFT) calculations predict selectivity for neutral 5-membered-ring intermediates over 6-membered-ring intermediates, as suggested by infrared and UV-vis spectroscopic results reported by others. Given that the products of aromatization are 6-membered rings, this result suggests that triene cyclization can only explain how linear hydrocarbons become cyclic intermediates but not the mechanisms that ultimately lead to the aromatic rings seen in industrial zeolite-catalyzed hydrocarbon processes.

Phosphonate-Modified UiO-66 Brønsted Acid Catalyst and Its Use in Dehydra-Decyclization of 2-Methyltetrahydrofuran to Pentadienes.

Phosphorus-modified all-silica zeolites exhibit activity and selectivity in certain Brønsted acid catalyzed reactions for biomass conversion. In an effort to achieve similar performance with catalysts having well-defined sites, we report the incorporation of Brønsted acidity to metal-organic frameworks with the UiO-66 topology, achieved by attaching phosphonic acid to the 1,4-benzenedicarboxylate ligand and using it to form UiO-66-PO3 H2 by post-synthesis modification. Characterization reveals that UiO-66-PO3 H2 retains stability similar to UiO-66, and exhibits weak Brønsted acidity, as demonstrated by titrations, alcohol dehydration, and dehydra-decyclization of 2-methyltetrahydrofuran (2-MTHF). For the later reaction, the reported catalyst exhibits site-time yields and selectivity approaching that of phosphoric acid on all-silica zeolites. Using solid-state NMR and deprotonation energy calculations, the chemical environments of P and the corresponding acidities are determined.

Adsorption in zeolites using mechanically embedded ONIOM clusters.

We have explored mechanically embedded three-layer QM/QM/MM ONIOM models for computational studies of binding in Al-substituted zeolites. In all the models considered, the high-level-theory layer consists of the adsorbate molecule and of the framework atoms within the first two coordination spheres of the Al atom and is treated at the M06-2X/6-311G(2df,p) level. For simplicity, flexibility and routine applicability, the outer, low-level-theory layer is treated with the UFF. We have modelled the intermediate-level layer quantum mechanically and investigated the performance of HF theory and of three DFT functionals, B3LYP, M06-2X and ωB97x-D, for different layer sizes and various basis sets, with and without BSSE corrections. We have studied the binding of sixteen probe molecules in H-MFI and compared the computed adsorption enthalpies with published experimental data. We have demonstrated that HF and B3LYP are inadequate for the description of the interactions between the probe molecules and the framework surrounding the metal site of the zeolite on account of their inability to capture dispersion forces. Both M06-2X and ωB97x-D on average converge within ca. 10% of the experimental values. We have further demonstrated transferability of the approach by computing the binding enthalpies of n-alkanes (C1-C8) in H-MFI, H-BEA and H-FAU, with very satisfactory agreement with experiment. The computed entropies of adsorption of n-alkanes in H-MFI are also found to be in good agreement with experimental data. Finally, we compare with published adsorption energies calculated by periodic-DFT for n-C3 to n-C6 alkanes, water and methanol in H-ZSM-5 and find very good agreement.

DFT Study of Solvent Effects in Acid-Catalyzed Diels-Alder Cycloadditions of 2,5-Dimethylfuran and Maleic Anhydride.

Density functional theory electronic structure calculations were used to explore the mechanism for the Diels-Alder reaction between 2,5-dimethylfuran and maleic anhydride (MA). Reaction paths are reported for uncatalyzed and Lewis and Brønsted acid-catalyzed reactions in vacuum and in a broad range of solvents. The calculations show that, while the uncatalyzed Diels-Alder reaction is thermally feasible in vacuum, a Lewis acid (modeled as Na(+)) lowers the activation barrier by interacting with the dienophile (MA) and decreasing the HOMO-LUMO gap of the reactants. A Brønsted acid (modeled as a proton) can bind to a carbonyl oxygen in MA, changing the reaction mechanism from concerted to stepwise and eliminating the activation barrier. Solvation effects were studied with the SMD model. Electrostatic effects play the largest role in determining the solvation energy of the transition state, which tracks the net dipole moment at the transition state. For the uncatalyzed reaction, the dipole moment is largely determined by charge transfer between the reactants, but in the reactions with ionic catalysts, there is no simple relationship between solvation of the transition state and charge transfer between the reactants. Nonelectrostatic contributions to solvation of the reactants and transition state also make significant contributions to the activation energy.

Solvent-induced frequency shifts of 5-hydroxymethylfurfural deduced via infrared spectroscopy and ab initio calculations.

Solvent-induced frequency shifts (SIFS) of the carbonyl stretching vibration ν(C═O) of 5-hydroxymethylfurfural were measured in protic, polar aprotic, and nonpolar solvents. The Gutmann acceptor number (AN) was found to correlate with the measured frequency shifts. The SIFS in six solvents were investigated using ab initio electronic structure calculations, treating the solvent implicitly and with an explicit solvent ligand interacting with the carbonyl. The conductor-polarizable continuum model (CPCM) of solvation predicted that ν(C═O) shifted according with the dielectric constant as (ε - 1)/(2ε + 1), in agreement with the analytical predictions of the Kirkwood-Bauer-Magat (KBM) theory for a dipole in a dielectric continuum, but in disagreement with the experimental trend. The experimental SIFS were best predicted using gas-phase complexes of HMF and explicit solvent-ligand. Natural bond orbital (NBO) analysis and Bader's atoms in molecules theory were used to investigate the electronic structure of these complexes. Strong SIFS were found to arise from stronger H-bonding interactions, as observed in delocalization of carbonyl lone-pair electrons by H-bonding solvent σ*(X-H) orbitals, and an increase in charge density and a decrease in local potential energy at the H-bond (3, -1) critical point. Consequently, by predicting the experimental SIFS and examining the electronic structure, we find the first theoretical evidence for treating Gutmann's solvent AN as a measure of solvent Lewis acidity.

Conformations of polyolefins on platinum catalysts control product distribution in plastics recycling.

The design of catalysts for the chemical recycling of plastic waste will benefit greatly from an intimate knowledge of the interfacial polymer-catalyst interactions that determine reactant and product distributions. Here, we investigate backbone chain length, side chain length, and concentration effects on the density and conformation of polyethylene surrogates at the interface with Pt(111) and relate them to experimental product distributions resulting from carbon-carbon bond cleavage. Using replica-exchange molecular dynamics simulations, we characterize the polymer conformations at the interface by the distributions of trains, loops, and tails and their first moments. We find that the preponderance of short chains, in the range of 20 carbon atoms, lies entirely on the Pt surface, whereas longer chains exhibit much broader distributions of conformational features. Remarkably, the average length of trains is independent of the chain length but can be tuned via the polymer-surface interaction. Branching profoundly impacts the conformations of long chains at the interface as the distributions of trains become less dispersed and more structured, localized around short trains, with the immediate implication of a wider carbon product distribution upon C-C bond cleavage. The degree of localization increases with the number and size of the side chains. Long chains can adsorb from the melt onto the Pt surface even in melt mixtures containing shorter polymer chains at high concentrations. We confirm experimentally key computational findings and demonstrate that blends may provide a strategy to reduce the selectivity for undesired light gases.

Insights into solvent and surface charge effects on Volmer step kinetics on Pt (111).

The mechanism of pH-dependent hydrogen oxidation and evolution kinetics is still a matter of significant debate. To make progress, we study the Volmer step kinetics on platinum (111) using classical molecular dynamics simulations with an embedded Anderson-Newns Hamiltonian for the redox process and constant potential electrodes. We investigate how negative electrode electrostatic potential affects Volmer step kinetics. We find that the redox solvent reorganization energy is insensitive to changes in interfacial field strength. The negatively charged surface attracts adsorbed H as well as H+, increasing hydrogen binding energy, but also trapping H+ in the double layer. While more negative electrostatic potential in the double layer accelerates the oxidation charge transfer, it becomes difficult for the proton to move to the bulk. Conversely, reduction becomes more difficult because the transition state occurs farther from equilibrium solvation polarization. Our results help to clarify how the charged surface plays a role in hydrogen electrocatalysis kinetics.

Highly active, ultra-low loading single-atom iron catalysts for catalytic transfer hydrogenation.

Highly effective and selective noble metal-free catalysts attract significant attention. Here, a single-atom iron catalyst is fabricated by saturated adsorption of trace iron onto zeolitic imidazolate framework-8 (ZIF-8) followed by pyrolysis. Its performance toward catalytic transfer hydrogenation of furfural is comparable to state-of-the-art catalysts and up to four orders higher than other Fe catalysts. Isotopic labeling experiments demonstrate an intermolecular hydride transfer mechanism. First principles simulations, spectroscopic calculations and experiments, and kinetic correlations reveal that the synthesis creates pyrrolic Fe(II)-plN3 as the active center whose flexibility manifested by being pulled out of the plane, enabled by defects, is crucial for collocating the reagents and allowing the chemistry to proceed. The catalyst catalyzes chemoselectively several substrates and possesses a unique trait whereby the chemistry is hindered for more acidic substrates than the hydrogen donors. This work paves the way toward noble-metal free single-atom catalysts for important chemical reactions.

Tuning the reactivity of carbon surfaces with oxygen-containing functional groups.

Oxygen-containing carbons are promising supports and metal-free catalysts for many reactions. However, distinguishing the role of various oxygen functional groups and quantifying and tuning each functionality is still difficult. Here we investigate the role of Brønsted acidic oxygen-containing functional groups by synthesizing a diverse library of materials. By combining acid-catalyzed elimination probe chemistry, comprehensive surface characterizations, 15N isotopically labeled acetonitrile adsorption coupled with magic-angle spinning nuclear magnetic resonance, machine learning, and density-functional theory calculations, we demonstrate that phenolic is the main acid site in gas-phase chemistries and unexpectedly carboxylic groups are much less acidic than phenolic groups in the graphitized mesoporous carbon due to electron density delocalization induced by the aromatic rings of graphitic carbon. The methodology can identify acidic sites in oxygenated carbon materials in solid acid catalyst-driven chemistry.

Understanding solvent effects in the selective conversion of fructose to 5-hydroxymethyl-furfural: a molecular dynamics investigation.

Selective conversion of fructose to 5-hydroxymethyl-furfural (HMF) involves the participation of high-boiling solvents like dimethyl sulfoxide (DMSO). In order to replace DMSO with low-boiling solvents, it is imperative that we understand the effect of DMSO solvation in protecting (i) HMF from rehydration and humins formation reactions and (ii) fructose from side reactions, other than its dehydration to HMF. In the present work, molecular dynamics simulations of HMF and fructose in water and in water-DMSO mixtures are carried out using the OPLS-AA force field. Radial pair distribution functions, coordination numbers and the hydrogen-bond network between the HMF/fructose molecule and the solvent molecules are analysed. The local 3-dimensional picture of the arrangement of solvent molecules around the solute, which cannot be accessed from pair distribution functions, is also computed. We show preferential coordination of DMSO around HMF and explain how this could provide a shielding effect to the HMF molecule, thus protecting it from further rehydration to levulinic acid and formic acid and from humins formation. In the case of fructose, the presence of DMSO also reduces the number of water molecules in the immediate vicinity of fructose. Though fewer water molecules coordinate around fructose, they are bound strongly to it. Analysis of the local 3-dimensional arrangement of DMSO molecules suggests that it protects the fructose molecule from side reactions that would lead to condensation or reversion products. However, the presence of DMSO molecules does not hamper the water molecules coming into contact with the oxygen atom of the hydroxyl groups of fructose, which is required for a proton transfer from water to fructose, to initiate the dehydration reaction to HMF.

Hybrid quantum mechanics/molecular mechanics-based molecular dynamics simulation of acid-catalyzed dehydration of polyols in liquid water.

We use the conversion of protonated glycerol to acrolein for a case study of the mechanism of acid-catalyzed dehydration of polyols in aqueous environments. We employ hybrid Quamtum Mechanics/Molecular Mechanics Molecular Dynamics (QM/MM MD) simulations with biased sampling and perform free energy calculations for the elementary steps of the reaction. We investigate the effects of solvent dynamics and in particular the role of quantum mechanical water in the dehydration mechanism. We present results supporting a mechanism that proceeds via water-mediated proton transfers and thus through an enol intermediate. We find that the first dehydration may take place by two, low-energy pathways requiring, respectively, 20.9 and 18.8 kcal/mol of activation free energy. The second dehydration requires 19.9 kcal/mol of activation free energy while for the overall reaction we compute a free energy change of -8 kcal/mol.

What controls au nanoparticle dispersity during growth?

Nanoparticles are critical in diverse nanotechnological areas. Their physical and chemical properties are essentially determined by their size, shape, and dispersity, but the underlying mechanisms that determine these are generally poorly understood. Herein, we focus on the electronic structure and solvation properties of the Au-thiolate capping agents (Au(n)MBA(n)) used in the synthesis of gold nanoparticles and whose accepted function is to stump growth. By DFT calculations we show that they exhibit "magic" number, electronic configuration, and stability for n = 4 and 8. By Molecular Dynamics simulations, we learn that in water-methanol solution the dominant species is the Au(4)MBA(4) complex, in a conformation stabilized by hydrogen bonds. We argue that the solvent has a multifunctional role: the water acts via hydrogen bonding as a "molecular locker" stabilizing particular conformers, while the methanol prevents polymeric growth of the Au(n)MBA(n) complexes. A plausible growth mechanism conducive to the formation of monodisperse gold nanoparticles is proposed.

Understanding mixing of Ni and Pt in the Ni/Pt(111) bimetallic catalyst via molecular simulation and experiments.

Molecular dynamics (MD) simulations employing embedded atom method potentials and ultrahigh vacuum (UHV) experiments were carried out to study the mixing process between the Ni and Pt atoms in the Ni/Pt(111) bimetallic system. The barrier for a Ni atom to diffuse from the top surface to the subsurface layer is rather high (around 1.7 eV) as calculated using the nudged elastic band (NEB) method. Analysis of the relaxation dynamics of the Ni atoms showed that they undergo diffusive motion through a mechanism of correlated hops. At 600 K, all Ni atoms remain trapped on the top surface due to large diffusion barriers. At 900 K, the majority of Ni atoms diffuse to the second layer and at 1200 K diffusion to the bulk is observed. We also find that smaller Ni coverages and the presence of Pt steps facilitate the Ni-Pt mixing. By simulated annealing simulations, we found that in the mixed state, the Ni fraction oscillates between layers, with the second layer being Ni-richer at equilibrium. The simulation results at multiple time scales are consistent with the experimental data.