An official website of the United States government
The chemical pathway for synthesizing covalent organic frameworks (COFs) involves a complex medley of reaction sequences over a rippling energy landscape that cannot be adequately described using existing theories. Even with the development of state-of-the-art experimental and computational tools, identifying primary mechanisms of nucleation and growth of COFs remains elusive. Other than empirically, little is known about how the catalyst composition and water activity affect the kinetics of the reaction pathway. Here, for the first time, we employ time-resolved in situ Fourier transform infrared spectroscopy (FT-IR) coupled with a six-parameter microkinetic model consisting of ∼10 million reactions and over 20 000 species. The integrated approach elucidates previously unrecognized roles of catalyst p K a on COF yield and water on growth rate and size distribution. COF crystalline yield increases with decreasing p K a of the catalysts, whereas the effect of water is to reduce the growth rate of COF and broaden the size distribution. The microkinetic model reproduces the experimental data and quantitatively predicts the role of synthesis conditions such as temperature, catalyst, and precursor concentration on the nucleation and growth rates. Furthermore, the model also validates the second-order reaction mechanism of COF-5 and predicts the activation barriers for classical and non-classical growth of COF-5 crystals. The microkinetic model developed here is generalizable to different COFs and other multicomponent systems.
more »
« less
Implementing the Blowers–Masel Approximation to Scale Activation Energy Based on Reaction Enthalpy in Mean-Field Microkinetic Modeling for Catalytic Methane Partial Oxidation
Mean-field microkinetic modeling is a powerful tool for catalyst design and the simulation of catalytic processes. The reaction enthalpies in a microkinetic model often need to be adjusted when changing species’ binding energies to model different catalysts, when performing thermodynamic sensitivity analyses, and when fitting experimental data. When altering reaction enthalpies, the activation energies should also be reasonably altered to ensure realistic reaction rates. The Blowers–Masel approximation (BMA) relates the reaction barrier to the reaction enthalpy. Unlike the Brønsted–Evans–Polani relationship, the BMA requires less data because only one parameter, the intrinsic activation energy, needs to be determined. We validate this application of BMA relations to model surface reactions by comparing against density functional theory data taken from the literature. By incorporating the BMA rate description into the open-source Cantera software, we enable a new workflow, demonstrated herein, allowing rapid screening of catalysts using linear scaling relationships and BMA kinetics within the process simulation software. For demonstration purposes, a catalyst screening for catalytic methane partial oxidation on 81 hypothetical metals is conducted. We compared the results with and without BMA-corrected rates. The heat maps of various descriptors (e.g., CH4 conversion, syngas yield) show that using BMA rates instead of Arrhenius rates (with constant activation energies) changes which metals are most active. Heat maps of sensitivity analyses can help identify which reactions or species are the most influential in shaping the descriptor map patterns. Our findings indicate that while using BMA-adjusted rates did not markedly affect the most sensitive reactions, it did change the most influential species.
more »
« less
- Award ID(s):
- 1931389
- PAR ID:
- 10513625
- Publisher / Repository:
- American Chemical Society
- Date Published:
- Journal Name:
- ACS Catalysis
- Volume:
- 14
- Issue:
- 10
- ISSN:
- 2155-5435
- Page Range / eLocation ID:
- 8013 to 8029
- Subject(s) / Keyword(s):
- mean-field microkinetic modeling Blowers–Masel approximation sensitivity analysis linear scaling catalytic methane partial oxidation
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Hydrogen gas is a promising renewable energy storage medium when produced via water electrolysis, but this process is limited by the sluggish kinetics of the anodic oxygen evolution reaction (OER). Herein, we used a microkinetic model to investigate promoting the OER using programmable oxide catalysts (i.e., forced catalyst dynamics). We found that programmable catalysts could increase current density at a fixed overpotential (100-600× over static rates) or reduce the overpotential required to reach a fixed current density of 10 mA cm-2 (45-140% reduction vs static). In our kinetic parametrization, the key parameters controlling the quality of the catalytic ratchet were the O*-to-OOH* and O*-to-OH* activation barriers. Our findings indicate that programmable catalysts may be a viable strategy for accelerating the OER or enabling lower-overpotential operation, but a more accurate kinetic parametrization is required for precise predictions of performance, ratchet quality, and resulting energy efficiency.more » « less
-
Per- and polyfluoroalkyl substances (PFAS) are a group of synthetic organofluorine surfactants that are resistant to typical methods of degradation. Thermal techniques along with other novel, less energy-intensive techniques are currently being investigated for the treatment of PFAS-contaminated matrices. Non-equilibrium plasma is one technique that has shown promise for the treatment of PFAS-contaminated water. To better tailor non-equilibrium plasma systems for this application, knowledge of the energy required for mineralization, and in turn the roles that plasma reactive species and heat can play in this process, would be useful. In this study, fundamental thermodynamic equations were used to estimate the enthalpies of reaction (480 kJ/mol) and formation (−4640 kJ/mol) of perfluorooctanoic acid (PFOA, a long-chain legacy PFAS) in water. This enthalpy of reaction estimate indicates that plasma reactive species alone cannot catalyze the reaction; because the reaction is endothermic, energy input (e.g., heat) is required. The estimated enthalpies were used with HSC Chemistry software to produce a model of PFOA defluorination in a 100 mg/L aqueous solution as a function of enthalpy. The model indicated that as enthalpy of the reaction system increased, higher PFOA defluorination, and thus a higher extent of mineralization, was achieved. The model results were validated using experimental results from the gliding arc plasmatron (GAP) treatment of PFOA or PFOS-contaminated water using argon and air, separately, as the plasma gas. It was demonstrated that PFOA and PFOS mineralization in both types of plasma required more energy than predicted by thermodynamics, which was anticipated as the model did not take kinetics into account. However, the observed trends were similar to that of the model, especially when argon was used as the plasma gas. Overall, it was demonstrated that while energy input (e.g., heat) was required for the non-equilibrium plasma degradation of PFOA in water, a lower energy barrier was present with plasma treatment compared to conventional thermal treatments, and therefore mineralization was improved. Plasma reactive species, such as hydroxyl radicals (⋅OH) and/or hydrated electrons (e−(aq)), though unable to accelerate an endothermic reaction alone, likely served as catalysts for PFOA mineralization, helping to lower the energy barrier. In this study, the activation energies (Ea) for these species to react with the alpha C–F bond in PFOA were estimated to be roughly 1 eV for hydroxyl radicals and 2 eV for hydrated electrons.more » « less
-
null (Ed.)Palladium catalyzed cross-coupling reactions represent a significant advancement in contemporary organic synthesis as these reactions are of strategic importance in the area of pharmaceutical drug discovery and development. Supported palladium-based catalysts are highly sought-after in carbon–carbon bond forming catalytic processes to ensure catalyst recovery and reuse while preventing product contamination. This paper reports the development of heterogeneous Pd-based bimetallic catalysts supported on fumed silica that have high activity and selectivity matching those of homogeneous catalysts, eliminating the catalyst's leaching and sintering and allowing efficient recycling of the catalysts. Palladium and base metal (Cu, Ni or Co) contents of less than 1.0 wt% loading are deposited on a mesoporous fumed silica support (surface area SA BET = 350 m 2 g −1 ) using strong electrostatic adsorption (SEA) yielding homogeneously alloyed nanoparticles with an average size of 1.3 nm. All bimetallic catalysts were found to be highly active toward Suzuki cross-coupling (SCC) reactions with superior activity and stability for the CuPd/SiO 2 catalyst. A low CuPd/SiO 2 loading (Pd: 0.3 mol%) completes the conversion of bromobenzene and phenylboronic acid to biphenyl in 30 minutes under ambient conditions in water/ethanol solvent. In contrast, monometallic Pd/SiO 2 (Pd: 0.3 mol%) completes the same reaction in three hours under the same conditions. The combination of Pd with the base metals helps in retaining the Pd 0 status by charge donation from the base metals to Pd, thus lowering the activation energy of the aryl halide oxidative addition step. Along with its exceptional activity, CuPd/SiO 2 exhibits excellent recycling performance with a turnover frequency (TOF) of 280 000 h −1 under microwave reaction conditions at 60 °C. Our study demonstrates that SEA is an excellent synthetic strategy for depositing ultra-small Pd-based bimetallic nanoparticles on porous silica for SCC. This avenue not only provides highly active and sintering-resistant catalysts but also significantly lowers Pd contents in the catalysts without compromising catalytic activity, making the catalysts very practical for large-scale applications.more » « less
-
The catalytic effects of iridium pincer complexes for the hydrogenation of carbon dioxide were investigated with the Unified Reaction Valley Approach (URVA), exploring the reaction mechanism along the reaction path traced out by the reacting species on the potential energy surface. Further details were obtained with the Local Mode Analysis performed at all stationary points, complemented by the Natural Bond Orbital and Bader’s Quantum Atoms in Molecules analyses. Each of the five reaction paths forming the catalytic cycle were calculated at the DFT level complemented with DLPNO-CCSD(T) single point calculations at the stationary points. For comparison, the non-catalytic reaction was also investigated. URVA curvature profiles identified all important chemical events taking place in the non-catalyzed reaction and in the five reactions forming the catalytic cycle, and their contribution to the activation energy was disclosed. The non-catalytic reaction has a large unfavorable activation energy of 76.3 kcal/mol, predominately caused by HH bond cleave in the H2 reactant. As shown by our study, the main function of the iridium pincer catalyst is to split up the one–step non-catalytic reaction into an energy efficient multistep cycle, where HH bond cleavage is replaced by the cleavage of a weaker IrH bond with a small contribution to the activation energy. The dissociation of the final product from the catalyst requires the cleavage of an IrO bond, which is also weak, and contributes only to a minor extent to the activation energy. This, in summary, leads to the substantial lowering of the overall activation barrier by about 50 kcal/mol for the catalyzed reaction. We hope that this study inspires the community to add URVA to their repertoire for the investigation of catalysis reactions.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1021/acscatal.3c05436
