Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reaction rate theory”

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.

At least 55 records · Page 3

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗

A theory of phonon-induced friction on molecular adsorbates

In this manuscript, we provide a general theory for how surface phonons couple to molecular adsorbates. Our theory maps the extended dynamics of a surface’s atomic vibrational motions to a generalized Langevin equation, and by doing so captures these dynamics in a single quantity: the non-Markovian friction. The different frequency components of this friction are the phonon modes of the surface slab weighted by their coupling to the adsorbate degrees of freedom. Using this formalism, we demonstrate that physisorbed species couple primarily to acoustic phonons while chemisorbed species couple to dispersionless local vibrations. We subsequently derive equations for phonon-adjusted reaction rates using transition state theory and demonstrate that these corrections improve agreement with experimental results for CO desorption rates from Pt(111).

Science & Technology - Other Topics↗

Univariate Prediction of Hammett Parameters and Select Relative Reaction Rates Using Loewdin Atomic Charges

Loewdin charges from density functional theory calculations were used here to obtain general, univariate linear correlations for the prediction of experimental Hammett parameters and relative reaction rates. While previous studies have established that Hirshfeld and CM5 charges perform strongly as univariate predictors, the near-ubiquitous Loewdin charges have not yet been evaluated. To this end, we assess the predictive capability of Loewdin charges for three chemical systems. First, we show that Loewdin charges outperform Hirshfeld and CM5 charges for Hammett parameter prediction. Second, we see that Loewdin charges generally perform comparably to Hirshfeld charges for predicting the relative rates of olefin cleavage by photoexcited nitroarenes. The single case of poor correlation, between relative rates and the Loewdin charges on nitrogen sites, is ameliorated when considering the net charge on the NO 2 group. Third, we show that Loewdin, Hirshfeld, and CM5 charges all perform very well for generating correlations for relative reaction rates for C–H activation of 9-(4-X-phenyl)-9H-fluorene substrates by a transition metal catalyst. The equations generated throughout the study enable the prediction of Hammett parameters and relative reaction rates. Finally, these tools can accelerate synthetic and experimental studies by enabling the in silico prediction of uncharacterized chemical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessment of Long-Term Degradation of Adsorbents for Direct Air Capture by Ozonolysis

Porous adsorbents are a promising class of materials for the direct air capture of CO 2 (DAC). Practical implementation of adsorption-based DAC requires adsorbents that can be used for thousands of adsorption–desorption cycles without significant degradation. We examined the potential degradation of adsorbents by a mechanism that appears to have not been considered previously, namely, ozonolysis by trace levels of ozone from ambient air. We focused on amine-appended metal–organic frameworks, specifically amine-functionalized Mg 2 (dobpdc), as a representative DAC adsorbent. Estimates based on the number of amine sites in these adsorbents and the ozone concentration in air suggest that degradation by ozone may be relevant over thousands of adsorption–desorption cycles if reactions with adsorbed ozone are fast. We used density functional theory calculations to estimate reaction rates for amine groups and carbon–carbon double bonds in amine-functionalized Mg 2 (dobpdc).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bayesian Estimation of the D(p,γ) 3 He Thermonuclear Reaction Rate

Big bang nucleosynthesis (BBN) is the standard model theory for the production of the light nuclides during the early stages of the universe, taking place for a period of about 20 minutes after the big bang. Deuterium production, in particular, is highly sensitive to the primordial baryon density and the number of neutrino species, and its abundance serves as a sensitive test for the conditions in the early universe. The comparison of observed deuterium abundances with predicted ones requires reliable knowledge of the relevant thermonuclear reaction rates, and their corresponding uncertainties. Recent observations reported the primordial deuterium abundance with percent accuracy, but some theoretical predictions based on BBN are at tension with the measured values because of uncertainties in the cross section of the deuterium-burning reactions. In this work, we analyze the S-factor of the D(p,γ) 3 He reaction using a hierarchical Bayesian model. We take into account the results of eleven experiments, spanning the period of 1955--2021; more than any other study. We also present results for two different fitting functions, a two-parameter function based on microscopic nuclear theory and a four-parameter polynomial. Furthermore, our recommended reaction rates have a 2.2\% uncertainty at 0.8~GK, which is the temperature most important for deuterium BBN. Differences between our rates and previous results are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Finite flame thickness effects on Kolmogorov-Petrovsky-Piskunov turbulent burning velocities

KPP (Kolmogorov-Petrovsky-Piskunov) solutions of the reaction-diffusion equation have application in various physical phenomena occurring in biology, ecology, and reacting flows. In particular, these solutions are commonly used in turbulent combustion to scale turbulent burning velocities. Subject to certain conditions on reaction rate profile through the flame brush and turbulent diffusivity, this theory relates the turbulent burning velocity to the derivative of the reaction rate ($\tilde{ω}$) at the leading edge of the flame brush ($d\tilde{ω}$/$d\tilde{c}$ | $\tilde{c} =$0 ). Such waves are often referred to as “pulled fronts.” However, turbulent flames never actually satisfy the KPP conditions for a pulled front, as the turbulent flame brush, parametrized here by the thickness δ t , consists of an ensemble of laminar flamelets of thickness δ, where ε=δ/δ t $\ll$ 1 is very small, but nonzero, and $d\tilde{ω}$ /$d\tilde{c}$ tends to zero at the brush leading edge for high activation energy, combustion-type kinetics. Here, this paper analyzes these effects on KPP wave solutions, parametrized by ε=δ/δ t and Zeldovich number Ze focusing on whether turbulent flames retain their pulled front character and what the correction to the KPP wave speed is. Variational solutions of the reaction-diffusion equation show that the solution can be expanded in powers of 1/|ln⁡ε|. Both numerical and asymptotic results are presented, showing that the wave still exhibits pulled front solutions but with significant corrections to the KPP result. The leading order correction is of the form |ln⁡ε| –2 and independent of Ze. Higher order corrections are function of both ε and Ze. However, the dominant factor influencing the wave speed correction is due to the finite ε, with Ze exhibiting a weaker effect.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computing Reaction Kinetics with MC-PDFT–OPESf: Combining Multireference Electronic Structure Theory and Enhanced Sampling

Accurate rate constants are crucial for understanding and optimizing catalytic reactions mediated by enzymes, metalloproteins, and heterogeneous catalysts. These systems frequently present a dual computational challenge. Multiconfigurational reaction sites require multireference techniques for the accurate treatment of the electronic structure, and high activation barriers prevent efficient sampling of unbiased reactive transitions. In this work, we combine multiconfiguration pair-density functional theory (MC-PDFT) as an accurate and efficient multireference electronic structure method with on-the-fly probability-enhanced sampling flooding (OPESf) as an enhanced sampling method capable of accelerating reactive transitions. We demonstrate the approach on the Diels–Alder [4+2] cycloaddition between cis-butadiene and ethene as a reaction characterized by a large activation barrier and multireference character. MC-PDFT–OPESf provides reaction rates in agreement with experiments at a fraction of the computational cost required by conventional unbiased ab initio calculations. Here, we propose MC-PDFT–OPESf as an efficient approach for computing kinetics in strongly correlated molecular systems.

Chemical calculations↗

The effect of temperature fluctuations of reaction rate constants in turbulent reacting flows

Current models of turbulent reacting flows frequently use Arrhenius reaction rate constants obtained from static or laminar flow theory and/or experiments, or from best fits of static, laminar, and turbulent data. By treating the reaction rate constant as a continuous random variable which is temperature-dependent, the present study assesses the effect of turbulent temperature fluctuations on the reaction rate constant. This model requires that a probability density function (PDF) describing the nature of the fluctuations be specified. Three PDFs are examined: the clipped Gaussian, the beta PDF, and the ramp model. All the models indicate that the reaction rate constant is greater in a turbulent flow field than in an equivalent laminar flow. In addition, an amplification ratio, which is the ratio of the turbulent rate constant to the laminar rate constant, is defined and its behavior as a function of the mean temperature fluctuations is described

Chinitz, W.↗

Estimation of the reaction rate for the formation of CH3O from H + H2CO - Implications for chemistry in the solar system

Troe's (1977) approximate theory is presently used in conjunction with transition state theory to estimate the rate coefficient of the reaction by which CO is reduced to CH4; attention is given to the role that may be played in the reduction process by the formation of the CH3O radical from H + H2CO. Attention is given to the implications of such a reaction (1) for the CO chemistry on Jupiter and within the solar nebula, (2) for the interpretation of such experimental results as those of Bar-Nun and Shaviv (1975) and Bar-Nun and Chang (1983), and (3) for organic synthesis in the prebiotic terrestrial atmosphere.

Yung, Yuk L.↗

Degrees of Rate Control in Interconnected Reaction Networks

Overall reactions in interconnected networks exhibit net, forward, and reverse rates that are governed by both constitutive elementary steps in the pathway of interest and branching elementary steps that lead to alternative products. Accordingly, steps in branching pathways exhibit negative net, forward, and reverse degrees of rate control, as they reduce reaction flux to the desired product. We here contextualize the forward and reverse degrees of rate control in terms of kinetic resistances (inverse of rates) and leverage the additive nature of kinetic resistance to decouple kinetic driving forces contributed by constitutive elementary steps and branching points (nodal species) in interconnected networks. Regardless of the network connectivity, forward and reverse degrees of rate control are shown to converge at equilibrium. Away from equilibrium, we identify two critical features of interconnected networks: stoichiometric regularity─condition where all stoichiometric numbers are unity─and pathway symmetry around nodal species─condition where branching pathways share the same rate constants, stoichiometry, and species concentrations/activities─that result in (i) equal forward, reverse, and consequently net degrees of rate control and (ii) forward and reverse degrees of rate control that exhibit constant offsets, respectively, across all extents of reaction. Furthermore, our discourse further provides a mathematical description for the influence of stoichiometric irregularity and pathway asymmetry on forward and reverse degrees of rate control. Altogether, the presented work details the effects of network (inter)connectivity and stoichiometry on reaction kinetics and, in doing so, establishes general protocols for capturing these effects as additive terms in the formulation of forward and reverse degrees of rate control.

10 SYNTHETIC FUELS↗

Revealing the structure of a catalytic combustion active-site ensemble combining uniform nanocrystal catalysts and theory insights

Supported metal catalysts are extensively used in industrial and environmental applications. To improve their performance, it is crucial to identify the most active sites. This identification is, however, made challenging by the presence of a large number of potential surface structures that complicate such an assignment. Often, the active site is formed by an ensemble of atoms, thus introducing further complications in its identification. Being able to produce uniform structures and identify the ones that are responsible for the catalyst performance is a crucial goal. Here, we utilize a combination of uniform Pd/Pt nanocrystal catalysts and theory to reveal the catalytic active-site ensemble in highly active propene combustion materials. Using colloidal chemistry to exquisitely control nanoparticle size, we find that intrinsic rates for propene combustion in the presence of water increase monotonically with particle size on Pt-rich catalysts, suggesting that the reaction is structure dependent. We also reveal that water has a near-zero or mildly positive reaction rate order over Pd/Pt catalysts. Theory insights allow us to determine that the interaction of water with extended terraces present in large particles leads to the formation of step sites on metallic surfaces. These specific step-edge sites are responsible for the efficient combustion of propene at low temperature. This work reveals an elusive geometric ensemble, thus clearly identifying the active site in alkene combustion catalysts. These insights demonstrate how the combination of uniform catalysts and theory can provide a much deeper understanding of active-site geometry for many applications

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of magnetic field on reduction of nickel oxide

Results of observations recorded during constant temperature reduction of NiO, a paramagnetic substance, to Ni, a ferromagnetic element, are presented. The application of a large magnetic field (4,200 oersted) does not result in an acceleration of the reduction rate. To explain earlier observations that 500- and 1,400-oersted magnetic fields increase the reduction rate of iron ore which contains hematite (Fe2O3), Peters (1973) had suggested that thermodynamics theory could predict the acceleration in reaction rate when product and reactant species differ widely in their magnetic properties.

Rowe, M. W.↗

Neutral gas pressure dependence of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations

In this computational study, we describe a self-consistent trajectory simulation approach to capture the effect of neutral gas pressure on ion–ion mutual neutralization (MN) reactions. The electron transfer probability estimated using Landau–Zener (LZ) transition state theory is incorporated into classical trajectory simulations to elicit predictions of MN cross sections in vacuum and rate constants at finite neutral gas pressures. Electronic structure calculations with multireference configuration interaction and large correlation consistent basis sets are used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of background gas on ion transport. For H + – H – and Li + – H(D) – , our approach quantitatively agrees with measured speed-dependent cross sections for up to ~10 5 m/s. For the ion pair Ne + – Cl – , our predictions of the MN rate constant at ~1 Torr are a factor of ~2 to 3 higher than the experimentally measured value. Similarly, for Xe + – F – in the pressure range of ~20 000–80 000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. Here, the paradigm of using trajectory simulations to self-consistently capture the effect of gas pressure on MN reactions advanced here provides avenues for the inclusion of additional nonclassical effects in future work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CN radical reactions with hydrogen cyanide and cyanogen - Comparison of theory and experiment

The method of laser photolysis/laser-induced fluorescence is used to obtain absolute rate constants for CN radical reactions with HCN and C2N2. The rate constants were found to be temperature-dependent in the range 300-740 K and pressure independent in the range 100-600 Torr. Rice-Remsperger-Kassel-Marcus theory for both reactions employing the transition state parameters obtained by the BAC-MP4 method are made. These calculations yielded reasonable results for the CN + HCN reaction, predicting both the temperature dependence and pressure independence. No pressure effect was observed in the pressure range 100-1000 Torr at temperatures below 900 K, confirming the experimental results.

Yang, D. L.↗

Semiclassical transition state theory through the lens of the restricted partition function

The wide adoption of transition state theory in resolving the rates of molecular processes relies on the simplification from reducing the formal and numerical expense of dynamics by a geometric constraint. Such a reduction is at odds with the uncertainty in localization that the uncertainty principle requires. While many forms of semiclassical transition state theory (SCTST) have been aimed at addressing this challenge, a popular approach has relied on resolving the underlying phase space structure of the exact rate formula to leverage Bohr–Sommerfeld quantization. Here, the Hernandez–Miller SCTST reframed the thermal rate formula into an integral of the so-called restricted partition function (RPF) over the action associated with the reaction. The density-of-state SCTST has reframed the rate formula in terms of the instanton’s density of states (DoS). Here, we show the relationship between the RPF-SCTST and the DoS-SCTST and derive the latter from the former. In this way, we help unify these branches of SCTST and provide a clearer formalism for future advances.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Model associated with: "Thermodynamic control on the decomposition of organic matter across different electron acceptors"

This model data package is associated with the publication “Thermodynamic control on the decomposition of organic matter across different electron acceptors” submitted to Soil Biology and Biochemistry (Zheng et al., 2023; https://doi.org/10.1016/j.soilbio.2024.109364).In this research, a thermodynamic modeling framework is built to flexibly incorporate both organic matter (OM) molecules and electron acceptors for estimating potential free energy release from various redox reactions and to further predict reaction rates based on Microbial Transition State Theory. The model package includes scripts for thermodynamic modeling and postprocessing. Input Fourier-transform ion cyclotron resonance (FTICR) data are from a previous experimental study (Boye et al., 2018), and model outputs are free energy predictions and stoichiometric coefficients associated with all possible redox reactions.This data package is associated with the project GitHub repository found at MM_bioenergetic_modeling.This data package contains four folders (Input_FTICR, Model, Output, and Output_processing), a file-level metadata (FLMD) csv, and a data dictionary (dd) csv. Please see Zheng_bioenergetic_modeling_flmd.csv for a list of all files contained in this data package and descriptions for each. The Zheng_bioenergetic_modeling_dd.csv file describes the csv column headers. The “Model” folder contains scripts to run energy balance calculations for each electron acceptor. The “Output” folder contains csv files with stoichiometric information from model simulations. And the "Output_processing" folder contains scripts for reaction rate calculations and to generate plots.

54 ENVIRONMENTAL SCIENCES↗