Density Functional Theory and Reaction Kinetics Studies of the Water-Gas Shift Reaction on Pt-Re Catalysts
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Several recent observational and laboratory studies of processes involved in polar stratospheric ozone loss have prompted a reexamination of aspect of out understanding for this key indicator of global change. To a large extent, our confidence in understanding and projecting changes in polar and global ozone is based on our ability to to simulate these process in numerical models of chemistry and transport. These models depend on laboratory-measured kinetic reaction rates and photlysis cross section to simulate molecular interactions. In this study we use a simple box-model scenario for Antarctic ozone to estimate the uncertainty in loss attributable to known reaction kinetic uncertainties. Following the method of earlier work, rates and uncertainties from the latest laboratory evaluation are applied in random combinations. We determine the key reaction and rates contributing the largest potential errors and compare the results to observations to evaluate which combinations are consistent with atmospheric data. Implications for our theoretical and practical understanding of polar ozone loss will be assessed.
Several recent observational and laboratory studies of processes involved in polar stratospheric ozone loss have prompted a reexamination of aspects of our understanding for this key indicator of global change. To a large extent, our confidence in understanding and projecting changes in polar and global ozone is based on our ability to simulate these processes in numerical models of chemistry and transport. The fidelity of the models is assessed in comparison with a wide range of observations. These models depend on laboratory-measured kinetic reaction rates and photolysis cross sections to simulate molecular interactions. A typical stratospheric chemistry mechanism has on the order of 50- 100 species undergoing over a hundred intermolecular reactions and several tens of photolysis reactions. The rates of all of these reactions are subject to uncertainty, some substantial. Given the complexity of the models, however, it is difficult to quantify uncertainties in many aspects of system. In this study we use a simple box-model scenario for Antarctic ozone to estimate the uncertainty in loss attributable to known reaction kinetic uncertainties. Following the method of earlier work, rates and uncertainties from the latest laboratory evaluations are applied in random combinations. We determine the key reactions and rates contributing the largest potential errors and compare the results to observations to evaluate which combinations are consistent with atmospheric data. Implications for our theoretical and practical understanding of polar ozone loss will be assessed.
A combined theoretical and experimental approach is presented that uses a comprehensive mean-field microkinetic model, reaction kinetics experiments, and scanning transmission electron microscopy imaging to unravel the reaction mechanism and provide insights into the nature of active sites for formic acid (HCOOH) decomposition on Au/SiC catalysts. All input parameters for the microkinetic model are derived from periodic, self-consistent, generalized gradient approximation (GGA-PW91) density functional theory calculations on the Au(111), Au(100), and Au(211) surfaces and are subsequently adjusted to describe the experimental HCOOH decomposition rate and selectivity data. It is shown that the HCOOH decomposition follows the formate (HCOO) mediated path, with 100% selectivity toward the dehydrogenation products (CO21H2) under all reaction conditions. An analysis of the kinetic parameters suggests that an Au surface in which the coordination number of surface Au atoms is 4 may provide a better model for the active site of HCOOH decomposition on these specific supported Au catalysts.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The dissociation and isomerization kinetics of the methyl ester combustion intermediates methoxycarbonyl radical (CH3O$\dot{C}$(=O)) and (formyloxy)methyl radical ($\dot{C}$H2OC(=O)H) are investigated theoretically using high-level ab initio methods and Rice–Ramsperger–Kassel–Marcus (RRKM)/master equation (ME) theory. Geometries obtained at the hybrid density functional theory (DFT) and coupled cluster singles and doubles with perturbative triples correction (CCSD(T)) levels of theory are found to be similar. We employ high-level ab initio wave function methods to refine the potential energy surface: CCSD(T), multireference singles and doubles configuration interaction (MRSDCI) with the Davidson–Silver (DS) correction, and multireference averaged coupled-pair functional (MRACPF2) theory. MRSDCI+DS and MRACPF2 capture the multiconfigurational character of transition states (TSs) and predict lower barrier heights than CCSD(T). The temperature- and pressure-dependent rate coefficients are computed using RRKM/ME theory in the temperature range 300–2500 K and a pressure range of 0.01 atm to the high-pressure limit, which are then fitted to modified Arrhenius expressions. Dissociation of CH3O$\dot{C}$(=O) to $\dot{C}$H3 and CO2 is predicted to be much faster than dissociating to CH3$\dot{O}$ and CO, consistent with its greater exothermicity. Isomerization between CH3O$\dot{C}$(=O) and $\dot{C}$H2OC(=O)H is predicted to be the slowest among the studied reactions and rarely happens even at high temperature and high pressure, suggesting the decomposition pathways of the two radicals are not strongly coupled. The predicted rate coefficients and branching fractions at finite pressures differ significantly from the corresponding high-pressure-limit results, especially at relatively high temperatures. Finally, because it is one of the most important CH3$\dot{O}$ removal mechanisms under atmospheric conditions, the reaction kinetics of CH3$\dot{O}$ + CO was also studied along the PES of CH3O$\dot{C}$(=O); the resulting kinetics predictions are in remarkable agreement with experiments.
Supported metal single atom catalysts have shown unique activity and selectivity for several reactions. Investigating the reaction mechanism on single atom catalysts and unambiguous assignment of the activity and selectivity requires catalysts with exclusively single atoms. However, it is challenging to prepare pure single atom catalysts with conventional impregnation methods, which typically result in a mixture of single atoms and nanoparticles. Here, we show that owing to different reaction mechanisms on single atoms and nanoparticles, reaction kinetics can serve as a surface sensitive characterization technique for quantifying their surface site fractions. As a case study, we use CO oxidation kinetics on Ir/MgAl2O4 to quantify the surface site fractions of single atoms and nanoparticles and the results are consistent with aberration-corrected scanning transmission electron microscopy, X-ray absorption fine structure and infrared spectroscopies. Additionally, by carefully choosing the reaction conditions, the activity of single atoms (or nanoparticles) can be made dominant, enabling detection of a small fraction, or allowing a study of their reaction mechanism on a catalyst containing a mixture of single atoms and nanoparticles. These results are general and could be applied to other systems where two types of sites have different reaction mechanisms.
The unimolecular dissociation and isomerization kinetics of the three methyl propanoate (MP) radicals, CH 3 CH 2 C(=O)O$\dot{C}$H 2 (MP-m), CH 3 $\dot{C}$HC(=O)OCH 3 (MP-α), and $\dot{C}$H 2 CH 2 C(=O)OCH 3 (MP-β), are theoretically investigated using high-level ab initio methods and the Rice–Ramsperger–Kassel–Marcus (RRKM)/master equation (ME) theory. Stationary-point energies are obtained using the coupled cluster singles and doubles with perturbative triples correction (CCSD(T)), multi-reference singles and doubles configuration interaction (MRSDCI) with the Davidson-Silver (DS) correction, and multi-reference averaged coupled pair functional (MRACPF2) theories. The isomerization barriers between the three radicals are predicted to be generally lower than the corresponding bond dissociation channels, leading to a strongly coupled reaction system in subsequent kinetics studies. The phenomenological temperature- and pressure-dependent rate coefficients are computed using the RRKM/ME theory over a temperature range of 500 to 2000 K and at a pressure range of 0.01 atm to the high-pressure limit, which are then fitted to modified Arrhenius expressions. The β-scission rate coefficients of MP-α to CH 3 CHC(=O) and CH 3 $\dot{O}$ are predicted to be the smallest because of its highest activation energy among all studied unimolecular reactions channels. Analysis of branching fractions shows that both MP-m and MP-α radicals mainly decompose to the bimolecular products CH 3 CH 2 $\dot{C}$(=O) and H 2 CO, whereas the MP-β radical primarily decomposes via cleavage of a C–C bond to form C 2 H 4 and CH 3 O$\dot{C}$(=O). The isomerization channels dominate at low temperatures, the branching fractions of which decrease with increasing temperature and become very minor at about 2000 K. Our accurate rate coefficients and branching fractions help to illuminate the unique combustion properties of MP.
Not Available
Development of mathematical models of heat and mass transfer in glass-melting furnaces began in the 1970s and progressed rapidly with advances in sophisticated experimental/numerical techniques and increasing computational power. Today, practically all newly built or rebuilt furnaces are optimized with these models to meet stringent quality requirements, reduce the unit costs of manufacturing, or control emissions. One remaining hurdle is to model the batch-to-glass conversion accurately enough to reliably assess the glass production rate. This article summarizes two key aspects of the batch-conversion modeling—the heat transfer and the kinetics of conversion—and reviews the current state-of-the-art approaches to simulating them. We critically examine the advantages of the commonly used heat transfer approach, but also explain that its predictive capabilities are significantly restricted by the dependence of batch thermal properties on the time-temperature history. We argue that kinetic approaches to the batch-conversion modeling would offer a significant improvement when coupled with the heat transfer approach. Finally, we summarize key areas requiring further research on the way toward a realistic model of the batch blanket.
A collaborative effort involving experiments, kinetic modeling, and computational fluid dynamics (CFD) was used to understand co-gasification of coal-biomass mixtures. The overall goal of the work was to determine the key reactive properties for coal-biomass mixed fuels. Sub-bituminous coal was mixed with biomass feedstocks to determine the fluidization and gasification characteristics of hybrid poplar wood, switchgrass and corn stover. It was found that corn stover and poplar wood were the best feedstocks to use with coal. The novel approach of this project was the use of a red mud catalyst to improve gasification and lower gasification temperatures. An important results was the reduction of agglomeration of the biomass using the catalyst. An outcome of this work was the characterization of the chemical kinetics and reaction mechanisms of the co-gasification fuels, and the development of a set of models that can be integrated into other modeling environments. The multiphase flow code, MFIX, was used to simulate and predict the hydrodynamics and co-gasification, and results were validated with the experiments. The reaction kinetics modeling was used to develop a smaller set of reactions for tractable CFD calculations that represented the experiments. Finally, an efficient tool was developed, MCHARS, and coupled with MFIX to efficiently simulate the complex reaction kinetics.
Precipitation of silver ions simultaneously with chloride and bromide ions in ethylene glycol at a mild temperature (e.g., 60 degrees C) has been successfully demonstrated for the synthesis of silver chlorobromide (AgClxBr1-x, 0 < x < 1) nanoparticles, which is realized by injecting a AgNO3 solution into a solution containing both halogen ions. The injection rate of the AgNO3 solution has been determined to be critical for controlling the uniformity of AgClxBr1-x nanoparticles. Time-resolved in situ high-energy synchrotron X-ray diffraction has been applied, for the first time, to quantitatively monitor the reaction kinetics of nanocrystal formation. The real-time results shed light on the fact that the injection rate of AgNO3 solution significantly influences the nucleation and growth processes, and thus the quality of resulting AgClxBr1-x nanoparticles. Specifically, fast injection enables the complete addition of AgNO3 solution to the reaction solution before the nucleation process starts, leading to a good separation of nucleation and growth and thus the formation of uniform AgClxBr1-x nanocubes with well-defined composition and narrow size distribution. By contrast, slow injection results in a continuous addition of AgNO3 solution to the reaction solution even after nucleation starts, leading to continuous multiple nucleation/growth processes and thus the formation of AgClxBr1-x nanoparticles with broad dimensional and morphological distributions.
The thermoreversible Diels–Alder (DA) reaction involving furan and maleimide precursors has been extensively exploited to develop reversible thermosets for circular economy and additive manufacturing. Recently, applications of photothermal nanoparticles, which can absorb light and generate nanoscale heat, have gained significant attention in driving the Diels–Alder (DA) reaction more efficiently. This approach can be utilized by various applications, for example, solar light-enhanced recycling, self-healing at targeted area, and photothermal 3-D printing. Here, in this work, we address an important fundamental question on the photothermal approach: how does the reaction kinetics of the photothermally driven DA reaction compared to that of the conventional heat driven one? We found that the forward DA reaction kinetics varies dramatically depending on a type of heat sources, i.e., light-induced photothermal heat or conventional heat, despite of a similar bulk temperature. We observed that the light-induced photothermal heat significantly enhances the DA reaction rate in comparison to the bulk heating in an inert environment. On the contrary, the photothermally driven forward DA reaction is significantly slower than that of the heat driven reaction in the presence of oxygen. This is attributed to the photoexcited electrons in photothermal nanoparticles that can generate singlet oxygen ( 1 O 2 ) in the presence of visible light. This singlet oxygen may facilitate an endoperoxide side reaction on the furan group, ultimately, retard the DA reaction from the furan group. These results reveal the effectiveness of photothermal nanoparticles in accelerating DA reaction while underscoring the need to minimize oxygen exposure to prevent unfavorable side reactions.
X-ray reflectivity and transmission electron microscopy (TEM) were used to characterize the morphological changes in thin film electrodes with alternating Ni and NiO layers during lithiation as a function of the Ni buffer layer thickness. Complete lithiation of the active NiO layers occurs only when the thickness of the Ni/NiO bilayers are less than 75 Å – a threshold value that is determined by the sum of the Ni quantity in the Ni/NiO bilayer of the multilayer stack. Thicker Ni/NiO bilayers present a kinetic barrier for lithium ion diffusion inside the stack resulting in partial lithiation of the multilayer electrodes in which only the top NiO layer lithiates. Lithiation of NiO layers in a multilayer stack also leads to an interface-specific reaction that is observed to increase the thicknesses of adjacent Ni layers by 3–4 Å and is associated with the formation of a low-density Li 2 O layer, corresponding to an interfacially-driven phase separation of the NiO. Rate dependent cyclic voltammetry studies reveal a linear relation between the peak current and scan rate suggesting that the lithiation kinetics are controlled by charge transfer resistance at the liquid–solid interface.
Not Available
In the past decade the notable effort placed on improving intrinsic electrochemical kinetics of platinum group metal (PGM)-free electrocatalysts for the oxygen reduction reaction (ORR) has led to a significant improvement in both performance and understanding of this class of electrocatalysts. However, a limited amount of this development and understanding has been undertaken using operando electrochemical diagnostics at the membrane electrode assembly (MEA) level. In this work, the operando ORR kinetics on an atomically dispersed iron-nitrogen-carbon ((AD)Fe-N-C) PGM-free electrocatalyst have been examined to extract the reaction order and the activation energy of the ORR. The experiments were carefully designed to ensure the stability/predictability of the electrocatalyst during the data collection process and thus validate the relevance of the values obtained for the aforementioned parameters. A kinetic model that considers a potential-dependent availability of active sites (θ) is proposed. Active site availability is shown to be a function of both the change in the oxidation state (nR) and the redox potential at which the metal center transitions from a higher oxidation state to a lower one (). The resulting model fitting parameters for nR and (0.71 and 0.788 V, respectively) obtained from the analysis of operando data correlate well with those from in situ X-ray absorption near edge structure measurements (nR=0.57) and in situ cyclic voltammetry measurements (0.75 V < <0.8 V) in the MEA environment. The resulting model provides an excellent fit of MEA performance across the range of pressures, temperatures, and potentials under which the data were collected.