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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Design and synthesis of aryl-functionalized carbazole-based porous coordination cages.

A subset of coordination cages have garnered considerable recent attention for their potential permanent porosity in the solid state. Herein, we report a series of functionalized carbazole-based cages of the structure type M12(R-cdc)12 (M = Cr, Cu, Mo) where the functional groups include a range of aromatic substituents. Single-crystal X-ray structure determinations reveal a variety of intercage interactions in these materials, largely governed by pi-pi stacking. Density functional theory for a subset of these cages was used to confirm that the nature of the increased stability of aryl-functionalized cages is a result of inter-cage ligand interactions.

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.

Understanding solvent effects on adsorption and protonation in porous catalysts.

Solvent selection is a pressing challenge in developing efficient and selective liquid phase catalytic processes, as predictive understanding of the solvent effect remains lacking. In this work, an attenuated total reflection infrared spectroscopy technique is developed to quantitatively measure adsorption isotherms on porous materials in solvent and decouple the thermodynamic contributions of van der Waals interactions within zeolite pore walls from those of pore-phase proton transfer. While both the pore diameter and the solvent identity dramatically impact the confinement (adsorption) step, the solvent identity plays a dominant role in proton-transfer. Combined computational and experimental investigations show increasingly favorable pore-phase proton transfer to pyridine in the order: water < acetonitrile < 1,4 - dioxane. Equilibrium methods unaffected by mass transfer limitations are outlined for quantitatively estimating fundamental thermodynamic values using statistical thermodynamics.

Design and synthesis of capped-paddlewheel-based porous coordination cages.

To leverage the structural diversity of metal-organic frameworks, the ability to controllably terminate them for the isolation of porous coordination cages is advantageous. However, the strategy has largely been limited to ligand termination methods, particularly for paddlewheel-based materials. Here, we show a paddlewheel-capping strategy can be employed to afford previously unattainable coordination cage structures that are mimetic of metal-organic framework pores.

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.

Nickel supported on nitrogen-doped carbon nanotubes as hydrogen oxidation reaction catalyst in alkaline electrolyte.

The development of a low-cost, high-performance platinum-group-metal-free hydroxide exchange membrane fuel cell is hindered by the lack of a hydrogen oxidation reaction catalyst at the anode. Here we report that a composite catalyst, nickel nanoparticles supported on nitrogen-doped carbon nanotubes, has hydrogen oxidation activity similar to platinum-group metals in alkaline electrolyte. Although nitrogen-doped carbon nanotubes are a very poor hydrogen oxidation catalyst, as a support, it increases the catalytic performance of nickel nanoparticles by a factor of 33 (mass activity) or 21 (exchange current density) relative to unsupported nickel nanoparticles. Density functional theory calculations indicate that the nitrogen-doped support stabilizes the nanoparticle against reconstruction, while nitrogen located at the edge of the nanoparticle tunes local adsorption sites by affecting the d-orbitals of nickel. Owing to its high activity and low cost, our catalyst shows significant potential for use in low-cost, high-performance fuel cells.

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.

Mechanism of Brønsted acid-catalyzed glucose dehydration.

We present the first DFT-based microkinetic model for the Brønsted acid-catalyzed conversion of glucose to 5-hydroxylmethylfurfural (HMF), levulinic acid (LA), and formic acid (FA) and perform kinetic and isotopic tracing NMR spectroscopy mainly at low conversions. We reveal that glucose dehydrates through a cyclic path. Our modeling results are in excellent agreement with kinetic data and indicate that the rate-limiting step is the first dehydration of protonated glucose and that the majority of glucose is consumed through the HMF intermediate. We introduce a combination of 1) automatic mechanism generation with isotopic tracing experiments and 2) elementary reaction flux analysis of important paths with NMR spectroscopy and kinetic experiments to assess mechanisms. We find that the excess formic acid, which appears at high temperatures and glucose conversions, originates from retro-aldol chemistry that involves the C6 carbon atom of glucose.

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.

Origin of 5-hydroxymethylfurfural stability in water/dimethyl sulfoxide mixtures.

In the present work, we combined vibrational spectroscopy with electronic structure calculations to understand the solvation of HMF in DMSO, water, and DMSO/water mixtures and to provide insights into the observed hindrance of HMF rehydration and aldol condensation reactions if it is dissolved in DMSO/water mixtures. To achieve this goal, the attenuated total reflection FTIR spectra of a wide composition range of binary and ternary mixtures were measured, analyzed, and compared to the findings of ab initio DFT calculations. The effect of solvent on the HMF C-O and O-H vibrational modes reveals significant differences that are ascribed to different intermolecular interactions between HMF and DMSO or water. We also found that DMSO binds to HMF more strongly than water, and interactions with the HMF hydroxyl group are stronger than those with the HMF carbonyl group. We also showed the preferential solvation of HMF C-O groups by DMSO if HMF is dissolved in DMSO/water mixed solvent. Frontier molecular orbital theory was used to examine the influence of the solvent on side reactions. The results show that HMF solvation by DMSO increases its LUMO energy, which reduces its susceptibility to nucleophilic attack and minimizes undesirable hydration and humin-formation reactions. This result, together with the preferential solvation of HMF by DMSO, provide an explanation for the enhanced HMF stability in DMSO/water mixtures observed experimentally.