Engineering topics
Fushimi, Rebecca
Publications and source records attributed to Fushimi, Rebecca.
Understanding Catalyst ‘Volcano’ Dependence Through Fermi-Level Controlled Kinetics Using Electronic Theory
The ubiquitous two-step Michaelis–Menten and Temkin–Boudart reaction mechanisms are extended to include the influence of the catalyst electronic subsystem in a 5-step mechanism. The resulting kinetic equation provides an alternative explanation for the well-known volcano-shaped dependence found in catalysis. The equilibrium constants of fast electronic steps are highlighted for their influence on adsorption and desorption through the relative concentration of charged versus neutral intermediates. This generalized concept can be widely applied to determine the optimal catalyst, based on the Fermi level of the material, for reactions proceeding via this universal reaction.
Model-Based Design of Experiments for Temporal Analysis of Products (TAP): A Simulated Case Study in Oxidative Propane Dehydrogenation
Not Available
Simulating Catalysis with Realistic Pellet Geometries Using Mesoflow: A Case Study of Catalytic Propane Dehydrogenation
We present a case study of catalytic propane dehydrogenation with our open-source multiphysics solver, Mesoflow. The solver was developed to simulate reactive flow coupled to heterogeneous catalytic reactions and deactivation in the context of complex, mesoscale geometry. The method leverages cartesian block-structured adaptive mesh refinement to capture realistic catalyst microstructural features acquired directly from X-ray computed tomography data. A kinetic model for propane dehydrogenation and catalyst deactivation was developed based on temporal analysis of products (TAP) reactor experiments. The TAP reactor experiments allow for precise characterization of intrinsic kinetic reaction steps which are implemented into Mesoflow simulations to model the spatial and temporal evolution of reactants, products, and catalyst active sites. The short-term and long-term deactivation behavior is studied by using XCT data collected from fresh and aged catalyst pellets, which exhibit different microstructural features. This study employs time-splitting algorithms to connect disparate reaction and flow timescales, enabling the simulations to achieve realistic deactivation timescales on the order of minutes while the flow time-scales for small particles (100 microns) are several milliseconds. We also introduce a flexible automated python script that writes the necessary files to construct a Mesoflow simulation from user-created chemical mechanisms. We will also introduce a few new features that are added to Mesoflow such as higher order schemes, implicit chemistry integrators and the ability to run on AMD and NVIDIA graphics-processing-units.
Carbon-Assisted, Continuous Syngas Production in a Chemical Looping Scheme
In the current energy and environment scenario, it is imperative to develop energy efficient routes for chemical manufacturing that also pave the way for mitigation of greenhouse gas emissions. This work presents an efficient pathway for continuous syngas production via a chemical looping conversion of the two most potent greenhouse gases—CH 4 , and CO 2 . The well-known dry-reforming process of converting CH 4 , and CO 2 to syngas is energy-intensive and suffers from catalyst deactivation. The chemical looping approach, on the other hand, provides avenues for mitigating catalyst deactivation and enabling improved energy efficiency. The key to such process enhancements lies in the intricate structure–function relationships of the catalyst and its correlation to the process variables. We present the reduction and oxidation characteristics of 5 wt.% Ni/Ce 1-x Zr x O 2 -based catalysts (x = 0, 0.4, and 0.625). We demonstrate low temperature CH4 activation over Ni-promoted samples as opposed to pure Ce 1-x Zr x O 2 . Moreover, our results depict an optimum regeneration of these catalysts when oxidized by CO 2 , and H 2 O, which allows for chemical looping operation of steam reforming of methane as well. Process variables were tuned to optimize the CH 4 conversion (over 80%), and H 2 /CO ratio at 650 °C. The critical surface reactions—carbon accumulation and gasification, and thermocatalytic CO 2 splitting were investigated to elucidate the dynamic nature of the catalyst surface. In conclusion, the impact of this work lies in showcasing the opportunities to design chemical looping reactors for energy efficient syngas production from waste greenhouse gases.
Deciphering the Mechanistic Role of Individual Oxide Phases and Their Combinations in Supported Mn–Na 2 WO 4 Catalysts for Oxidative Coupling of Methane
Oxidative coupling of methane (OCM) is an attractive direct route for upgrading methane to valuable chemicals. In this study, Temporal Analysis of Products (TAP) and steady state experiments are conducted to understand the role of individual oxide phases and their combinations in supported Mn-Na 2 WO 4 /SiO 2 catalysts for OCM. The results from TAP transient kinetic studies indicate that Mn plays an important role in promoting gas phase oxygen activation, while NaO x /SiO 2 and WO x /SiO 2 are relatively inert towards gas phase oxygen and methane activation. However, the supported catalyst combining Na and W in the form Na 2 WO 4 show enhanced gas phase oxygen activation exhibiting a much lower oxygen activation energy (148 kJ/mol) and enhanced activity toward methane activation as compared to the individual supported oxide catalysts. Addition of Mn to Na 2 WO 4 /SiO 2 further decreases the oxygen activation energy by 40 kJ/mol. Moreover, methane activation is also enhanced with CH 3 as the main intermediate but, with increasing Mn content, more CH 2 intermediates are observed. Different forms of oxygen (both dioxygen and atomic) are detected on the catalyst surface using isotopic pump/probe pulsing and their distribution is found to depend on the catalyst composition. An optimal Mn content in the Na 2 WO 4 /SiO 2 catalyst system is needed to enhance the amount of dioxide surface species (e.g., superoxide 16 O 2 - or peroxide 16 O 2 2- ) associated with the Na 2 WO 4 leading to high C 2 selectivity for OCM. When the Mn content is too high, the larger MnOx domains are shown to contribute to the formation of higher concertation of monoxide surface species that lead to nonselective OCM pathways. This insight from transient kinetic characterization using TAP combined with conventional steady state studies, provides a deeper understanding of the role of individual oxide phases and their combination on supported catalysts toward the formation of intermediate surface species and their impact on the OCM reaction mechanism. This knowledge is critical toward designing superior catalyst formulations for OCM.
TAPsolver: A Python package for the simulation and analysis of TAP reactor experiments
Not Available
Formation of Surface Impurities on Lithium–Nickel–Manganese–Cobalt Oxides in the Presence of CO 2 and H 2 O
Surface impurities involving parasitic reactions and gas evolution contribute to the degradation of high Ni content LiNi x Mn y Co z O 2 (NMC) cathode materials. The transient kinetic technique of Temporal Analysis of Products (TAP), density functional theory, and infrared spectroscopy have been used to study the formation of surface impurities on varying nickel content NMC materials (NMC811, NMC622, NMC532, NMC433, NMC111) in the presence of CO 2 and H 2 O. CO 2 reactivity on a clean surface as characterized by CO 2 conversion rate in the TAP reactor follows the order: NMC811 > NMC622 > NMC532 > NMC433 > NMC111. The capacity of CO 2 uptake follows a different order: NMC532 > NMC433 > NMC622 > NMC811 > NMC111. Moisture pretreatment slows down the direct CO 2 adsorption process and creates additional active sites for CO 2 adsorption. Electronic structure calculations predict that the (012) surface is more reactive than the ($10\bar{1}4$) surface for CO 2 and H 2 O adsorption. CO 2 adsorption leading to carbonate formation is exothermic with formation of ion pairs. The average CO 2 binding energies on the different materials follow the CO 2 reactivity order. Water hydroxylates the (012) surface and surface OH groups favor bicarbonate formation. Water creates more active sites for CO 2 adsorption on the ($10\bar{1}4$) surface due to hydrogen bonding. The composition of surface impurities formed in ambient air exposure is dependent on water concentration and the percentage of different crystal planes. Different surface reactivities suggest that battery performance degradation due to surface impurities can be mitigated by precise control of the dominant surfaces in NMC materials.
Pt-Assisted Carbon Remediation of Mo2C Materials for CO Disproportionation
Using the CO disproportionation (Boudouard) reaction as a probe reaction, an in-depth analysis of temperature-programmed pulse response data shows that the addition of Pt to Mo2C mitigates deactivation of Mo active sites by acting as a carbon collector. CO2 production on Mo2C and Pt/Mo2C materials is dependent on both the activation energy and the CO surface concentration. Detailed plane-wave density functional theory calculations of the CO adsorption and disproportion reactions on Mo2C-supported Pt nanoparticles (NPs) are reported. The Mo2C was modeled by the ß-Mo2C (100) surface, and the Pt/Mo2C interface was modeled by the addition of 12 Pt atoms to the Mo2C (100) surface (12Pt@Mo2C). The potential energy surfaces of the Boudouard reaction were calculated on pure Mo2C, 12Pt@Mo2C, and Pt (111) surfaces. CO dissociation readily occurs on the Mo2C (100) surface, but not on the Pt (111) surface, with the former being exothermic and the latter being endothermic. At the Pt/Mo2C interface, CO dissociation is still exothermic, but with a larger energy barrier. The Boudouard reaction takes place on the Mo2C region, where CO2 is formed from a surface O atom dissociated from one CO molecule in reaction with another CO molecule, leaving one C atom on the surface. C adsorption is preferential on the Pt site in comparison to the Mo site. The supported Pt domains can collect the remaining C atoms, facilitating further CO2 formation on the active Mo sites. A Bader charge analysis shows that the surface metal-carbon bond is a mixture of covalent and ionic bonds, whereas the surface metal-oxygen bond is ionic. Electron localization function (ELF) and partial charge density calculations agree well with the Bader charge analysis. These computational results are consistent with experimental observations of the interaction of CO with Mo2C nanotube supported Pt domains in the transient regime under far from equilibrium conditions. The Boudouard reaction is an important side reaction, and the unexpected role found for Pt as a carbon collector, with Mo serving as a disproportionation site, provides a unique vantage point for understanding carbon and coke formation on catalytic materials.
Rate/Concentration Kinetic Petals: A Transient Method to Examine the Interplay of Surface Reaction Processes
Transient pulse response experiments are used to construct rate/concentration (RC) kinetic dependencies—RC petals—and provide a new method to distinguish the timing and interplay of adsorption, surface reaction, and product formation on complex (industrial) materials. A petal shape arises as the dynamic “reaction-diffusion” experiment forces the concentration and reaction rate to return to zero. In contrast to the typical steady-state “Langmuir-type” RC dependence, RC petals have two branches, which arise as a result of decoupled gas and surface concentrations in the non-steady-state regime. To demonstrate this approach, the characteristics of petal shapes using ammonia decomposition as a probe reaction are presented. Ammonia, hydrogen, and nitrogen transformation rates are compared on three simple materials: Fe, Co, and a bimetallic CoFe preparation when ammonia is pulsed at 550°C in a low-pressure diffusion reactor. All materials demonstrate a two-branch kinetic RC dependence for ammonia adsorption, and rate constants are quantified in the low-coverage regime. We found that H 2 and N 2 product formation was dependent on the concentration of surface intermediates for all materials with one exception: for Co, an additional fast hydrogen generation process was observed; the rate of which coincided with ammonia adsorption. Nitrogen generation was only significant for CoFe and Co and on the CoFe catalyst, a self-inhibition property was observed. A method for estimating the number of active sites based on the RC petals is presented and was applied to the Fe and CoFe samples. The surface coverage and rate of formation/conversion of surface intermediates are interpreted from the examination of shape characteristics of the RC petals for each material.
Recent advances in dynamic chemical characterization using Temporal Analysis of Products
Not Available
Pulse response analysis using the Y-procedure: A data science approach
We report that pulse response experiments can provide a wealth of data for transient kinetic studies of gas phase heterogeneous reactions, with the ability to paint a detailed picture of intrinsic reaction kinetics pertaining to complex catalytic mechanisms. A novel framework called the Y-Procedure was developed by Yablonsky et al. to extract information from non-steady-state kinetic data (Yablonsky et al., 2007). The Y-Procedure reconstructs the chemical transformation rate and gas composition in the catalytic zone with no a priori assumptions on the kinetic model and is used to relate changes in catalytic reaction rates and kinetic parameters to changes in a catalyst’s surface composition. This paper provides a detailed description of the data science work flow implemented in the Y-procedure through an experimental example of oxygen adsorption on platinum. Lastly, this powerful method, along with a robust list of pre-processing methodology, will enable deeper understanding of kinetic processes and mechanism to be derived from the TAP pulse response technique.
Extracting Knowledge from Data through Catalysis Informatics
Catalysis informatics is a distinct sub-field that lies at the intersection of cheminformatics and materials informatics, but with unique challenges arising from the dynamic, surface-sensitive, multi-scale nature of heterogeneous catalysis. The ideas behind catalysis informatics can be traced back decades, but the field is only recently emerging due to advances in data infrastructure, statistics, machine learning, and computational methods. In this work we review the field from early works on expert systems and knowledge engines to more recent approaches utilizing machine-learning and uncertainty quantification. The data-information-knowledge hierarchy is introduced and used to classify various developments. The chemical master equation and micro-kinetic models are proposed as a quantitative representation of catalysis knowledge, which can be used to generate explanative and predictive hypotheses for the understanding and discovery of catalytic materials. We discuss future prospects for the field, including improved quantitative coupling of experiment/theory, advanced micro-kinetic models, and the development of open-source software tools. As a result, integration of existing chemical and physical models with emerging statistical and computational tools presents a promising route toward the automated design, discovery, and optimization of heterogeneous catalytic processes.
Precise non-steady-state characterization of solid active materials with no preliminary mechanistic assumptions
This paper presents a straightforward and user-friendly procedure for extracting a reactivity characterization of catalytic reactions on solid materials under non-steady-state conditions, particularly in temporal analysis of products (TAP) experiments. The kinetic parameters derived by this procedure can help with the development of detailed mechanistic understanding. The procedure consists of the following two major steps: 1) Three “Laplace reactivities” are first determined based on the moments of the exit flow pulse response data; 2) Depending on a select kinetic model, kinetic constants of elementary reaction steps can then be expressed as a function of reactivities and determined accordingly. In particular, we distinguish two calculation methods based on the availability and reliability of reactant and product data. The theoretical results are illustrated using a reverse example with given parameters as well as an experimental example of CO oxidation over a supported Au/SiO 2 catalyst. The procedure presented here provides an efficient tool for kinetic characterization of many complex chemical reactions.
Active site densities, oxygen activation and adsorbed reactive oxygen in alcohol activation on npAu catalysts
The activation of molecular O 2 as well as the reactivity of adsorbed oxygen species is of central importance in aerobic selective oxidation chemistry on Au-based catalysts. Herein, we address the issue of O 2 activation on unsupported nanoporous gold (npAu) catalysts by applying a transient pressure technique, a temporal analysis of products (TAP) reactor, to measure the saturation coverage of atomic oxygen, its collisional dissociation probability, the activation barrier for O 2 dissociation, and the facility with which adsorbed O species activate methanol, the initial step in the catalytic cycle of esterification. The results from these experiments indicate that molecular O 2 dissociation is associated with surface silver, that the density of reactive sites is quite low, that adsorbed oxygen atoms do not spill over from the sites of activation onto the surrounding surface, and that methanol reacts quite facilely with the adsorbed oxygen atoms. In addition, the O species from O 2 dissociation exhibits reactivity for the selective oxidation of methanol but not for CO. The TAP experiments also revealed that the surface of the npAu catalyst is saturated with adsorbed O under steady state reaction conditions, at least for the pulse reaction.