Search NASA⌕ Search

SEARCH · Search NASA

Results for “equilibrium constant”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Measurement of Palladium Hydride and Palladium Deuteride Isotherms Between 130 K and 393 K

The data collected in the present work extend the measured phase diagram for palladium hydride and palladium deuteride to a region that has been sparsely reported in open literature. Absorption isotherms were measured using a 2.5-g bed of palladium powder at temperatures between 130 K and 393 K and pressures less than 1.3 x 10 5 Pa. Such low-pressure and low-temperature measurements are useful for characterizing palladium beds used for tritium pumping and storage. For tritium storage, pressures are kept below a few millibar for safety reasons. Low temperatures increase the tritium storage capacity of palladium. The measured absorption isotherms show the well-documented, two-phase behavior for this system: two solubility regions and a mixed, hydride-forming region. The isotherms show that an increased quantity of hydride is formed at lower temperatures, as marked by an increase in the hydride-forming region. This region exceeds hydrogen to metal ratios of 0.75 for T ≤ 273 K. Equilibrium pressures in the mixed region decrease with decreasing temperatures until a critical temperature is reached for each isotope. Below these critical temperatures, the rate of pressure decrease with decreasing temperature is significantly reduced. This change in trend suggests hydrogen isotopes are adsorbed onto the palladium surface, rather than forming a hydride. Using the equilibrium pressures recorded at temperatures between 236 K and 393 K for protium and between 211 K and 354 K for deuterium, the van’t Hoff constants were calculated to be A = –36 ± 1 kJ/mol and B = 88 ± 3 J/K for protium and A = -32 ± 2 kJ/mol and B = 88 ± 9 J/K for deuterium. As a result, these constants agree favorably with literature in the range where the temperatures of the measured isotherms overlap.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Revealing the Drag Instability in One-fluid Nonideal Magnetohydrodynamic Simulations of a 1D Isothermal C-shock

C-type shocks are believed to be ubiquitous in turbulent molecular clouds thanks to ambipolar diffusion. We investigate whether the drag instability in 1D isothermal C-shocks, inferred from the local linear theory of Gu & Chen, can appear in nonideal magnetohydrodynamic simulations. Two C-shock models (with narrow and broad steady-state shock widths) are considered to represent the typical environment of star-forming clouds. The ionization-recombination equilibrium is adopted for the one-fluid approach. In the 1D simulation, the inflow gas is continuously perturbed by a sinusoidal density fluctuation with a constant frequency. The perturbations clearly grow after entering the C-shock region until they start being damped at the transition to the post-shock region. We show that the profiles of a predominant Fourier mode extracted locally from the simulated growing perturbation match those of the growing mode derived from the linear analysis. Moreover, the local growth rate and wave frequency derived from the predominant mode generally agree with those from the linear theory. Therefore, we confirm the presence of the drag instability in simulated 1D isothermal C-shocks. We also explore the nonlinear behavior of the instability by imposing larger-amplitude perturbations to the simulation. We find that the drag instability is subject to wave steepening, leading to saturated perturbation growth. Issues concerning local analysis, nonlinear effects, one-fluid approach, and astrophysical applications are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

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↗
DOE OSTI2020

Ion Conductivity and Correlations in Model Salt-Doped Polymers: Effects of Interaction Strength and Concentration

Correlated anion and cation motion can greatly reduce the overall ion conductivity in electrolytes versus the ideal conductivity calculated based on the diffusion constants alone. Using coarse-grained molecular dynamics simulations, we calculate conductivity and the degree of uncorrelated ion motion in salt-doped homopolymers and block copolymers as a function of concentration and interaction strengths. Calculating conductivity from ion mobility under an applied electric field increases accuracy versus the typical use of fluctuation dissipation relationships in equilibrium simulations. In typical electrolytes, correlation in cation-anion motion is often expected to be reduced at low ion concentrations. Yet, for these polymer electrolytes with strong ion-polymer and ion-ion interactions, we find cation-anion motion is more correlated at lower concentrations when other variables are held constant. We show this phenomenon is related to the slower ion cluster relaxation rate at low concentrations rather than the static spatial state of ion aggregation or the fraction of free ions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Toward Rapid Actinium-225 Purification via Membrane Adsorbers with Covalently Tethered Diglycolamide Ligands

Extractive diglycolamide (DGA) resins are used in several state-of-the-art techniques for purifying 225 Ac, a promising radiometal for targeted alpha therapy. Unfortunately, separation processes that rely on resins are often limited to slow flow rates, high elution volumes, and long processing times. Membrane adsorbers functionalized with DGA ligands are an alternative separation material that may overcome these challenges. This work presents (1) the synthesis of an aminated tetrahexyldiglycolamide ligand, (2) the covalent tethering of the ligand to electrospun poly(vinylbenzyl chloride) fiber mats, and (3) the adsorption and desorption of La(III) and 225 Ac. Chemical and physical characterization supports the covalent tethering of the ligand to the fiber mat, as well as the preservation of the fiber surface area and porosity after functionalization. Equilibrium adsorption experiments were performed with stable La(III) and radioactive 225 Ac. Trends in affinity are consistent between commercial resins and the synthesized membrane adsorbers; however, the Langmuir constants and the maximum binding capacity of the membrane adsorbers were generally lower than the resins. Despite these differences, the modeled selectivity for an equimolar solution of La(III)/ 225 Ac in 10 M nitric acid is 57. Furthermore, 225 Ac is rapidly desorbed from the fibers in 10 M nitric acid (<20 min). The La(III)/ 225 Ac selectivity and rapid 225 Ac desorption indicate this class of materials is promising for rapid radioanalytical separations.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

A Numerical Modeling Framework for Flocculation and Cohesive Sediment Transport in the Wave Bottom Boundary Layer

Flocculation, a critical process in coastal and estuarine systems, plays a significant role in sediment transport, nutrient cycling, and ecological health. This study develops a cohesive sediment transport modeling framework tailored to the wave bottom boundary layer under dilute and equilibrium conditions, explicitly incorporating flocculation effects via a Population Balance Equation (PBE). Using Direct Numerical Simulation, six baseline cases, each with a distinct sediment concentration profile resulting from a constant settling velocity and critical erosion shear stress, are generated to drive the PBE flocculation model for given floc yield strength and stickiness. Results reveal that flocculation significantly influences sediment concentration profiles promoting three distinct stages, well‐mixed, transition to lutocline, and well‐developed lutocline. At low concentrations with well‐mixed profiles, cohesive floc properties are less significant, and turbulence is a main flocculation driver. In contrast, as concentration increases, cohesive floc properties become crucial, facilitating lutocline formation. Here, the analysis also highlights limitations of depth‐averaged settling velocity as a parameterization. It is suitable for well‐mixed and transitional profiles but fails in well‐developed lutoclines, where empirical formulations that explicitly incorporate turbulent shear rate and sediment concentration better capture variability. This study underscores the necessity of incorporating flocculation effects into sediment transport models to enhance predictions of sediment dynamics in wave bottom boundary layers.

Penaloza‐Giraldo, Jorge A. [Oak Ridge National Lab↗
DOE OSTI2025

Dynamic response of a freely rotating butterfly valve in the advanced test reactor − dynamic coefficients modeling

Here, in evaluating the water hammer issue pertaining to the primary-coolant-regulating butterfly valve in the Advanced Test Reactor, the dynamic fluid body interaction (DFBI) approach was implemented in the analysis covered in Part I. Although DFBI modeling accurately and simultaneously solved the dynamic motion of the valve’s disk along with the flow field of the surrounding fluid, it shed little light on the reason behind such motion. For Part II, the reacting torque of the fluid on the disk was decomposed into representations of the dynamic coefficients in terms of stiffness, damping, and added mass. These were evaluated via simulations with steady-state static (stiffness), constant angular speed (damping), and variable angular speed (added mass) disks. Substituting the dynamic coefficients into Newton’s second law enabled the response trajectories to be obtained. Stable (by average) and unstable equilibrium positions and thrust tendencies of the valve were determined based on the stiffness coefficient (or static torque), the response amplitude was dampened or enlarged by the damping coefficient (minorly affected by added mass), and the response frequency was altered by the damping and added mass coefficients. Although the dynamic coefficient approach renders slightly different trajectories, due to the averaging effect of the torque in comparison to the DFBI method, the overall trend of the response aligns with the DFBI simulation, thus confirming the conclusion in Part I that a fix to the current butterfly valve is necessary.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Defluorination Mechanisms and Real-Time Dynamics of Per- and Polyfluoroalkyl Substances on Electrified Surfaces

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants found in groundwater sources and a wide variety of consumer products. In recent years, electrochemical approaches for the degradation of these harmful contaminants have garnered a significant amount of attention due to their efficiency and chemical-free modular nature. However, these electrochemical processes occur in open, highly non-equilibrium systems, and a detailed understanding of PFAS degradation mechanisms in these promising technologies is still in its infancy. To shed mechanistic insight into these complex processes, we present the first constant-electrode potential (CEP) quantum calculations of PFAS degradation on electrified surfaces. These advanced CEP calculations provide new mechanistic details about the intricate electronic processes that occur during PFAS degradation in the presence of an electrochemical bias, which cannot be gleaned from conventional density functional theory calculations. We complement our CEP calculations with large-scale ab initio molecular dynamics simulations in the presence of an electrochemical bias to provide time scales for PFAS degradation on electrified surfaces. Taken together, our CEP-based quantum calculations provide critical reaction mechanisms for PFAS degradation in open electrochemical systems, which can be used to prescreen candidate material surfaces and optimal electrochemical conditions for remediating PFAS and other environmental contaminants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Theoretical Understanding of Effects of Operating Modes on the Performance Durability of Solid Oxide Cells: A Comparison between Potentiostatic and Galvanostatic Operations

In solid oxide cells (SOCs), the choice between galvanostatic (constant current) and potentiostatic (constant voltage) modes does not significantly affect the performance of SOCs as long as the cell components remain unchanged and intact. However, the degradation of cell components, which leads to changes in electrochemical and physical properties of the cell, elevates the importance of the selected operating mode. This paper aims to investigate the effects of galvanostatic and potentiostatic operating modes on the evolving properties of SOCs and their subsequent influence on performance durability. Employing non-equilibrium thermodynamic analysis, a crucial approach for understanding the degradation phenomena within an active electrochemical system, this study aims to provide in-depth insights into how these operating modes affect the longevity and efficacy of SOCs. Key findings include: In cases where oxygen electrode (OE) degradation is accelerated by higher partial pressure of oxygen ( p O 2 ), operating under constant voltage electrolysis can mitigate the high p O 2 at the OE|electrolyte (OE|EL) interface. Conversely, if OE degradation occurs more rapidly under a lower p O 2 , constant current electrolysis is more effective in suppressing degradation by achieving a high p O 2 at the OE|EL interface. For degradation of the fuel electrode (FE) due to higher p O 2 , constant current electrolysis is beneficial for more stable performance, which helps maintain low p O 2 at the FE|EL interface. When FE degradation is accelerated by lower p O 2 , constant voltage electrolysis can avert low p O 2 at the FE|EL interface. In practical scenarios, more complex degradation mechanisms come into play, especially when p O 2 significantly deviates from initial conditions. Degradation in one electrode can influence p O 2 in the other electrode, a phenomenon more pronounced in potentiostatic than in galvanostatic electrolysis.

Electrochemistry↗
DOE OSTI2022

Numerical Investigation of the Spark Discharge Process in a Crossflow

The present study numerically investigates the spark discharge process under crossflow conditions using a thermal equilibrium plasma solver that fully couples the electromagnetic physics and fluid dynamics in a computational framework. Numerical results are validated by the comparison with experimental data. Here, the spark discharge experiment is performed in a constant volume vessel using an inductive coil ignition system for automotive applications, and the evolution of the spark channel is measured using high-speed imaging. The crossflow in the gap between the spark-plug electrodes is generated by a rotating fan with two different fan speeds, and the flow velocity across the gap is characterized by particle image velocimetry (PIV) measurement. A computational fluid dynamics (CFD) solver is employed to simulate the crossflow and provide the flow field variables (velocity, pressure, temperature) to the plasma solver. The crossflow velocity predicted in the flow simulation agrees well with the PIV data in that the non-uniform velocity profiles at monitoring points are reproduced by the CFD code. With the crossflow initialization in the plasma solver, the simulated spark discharge process from the breakdown to spark discharge matches the experimental data, including the voltage and circuit waveforms and the high-speed images of the spark channel evolution. The stretch of spark channel captured by plasma simulations agrees with the measured data. The plasma simulation reveals that the mean temperature of the spark channel is maintained at 5000 K during the discharge phase, and the temperature varies along the spark channel so that the highest value is obtained at the spark root on the center electrode. Overall, the results presented in this paper are meant to provide valuable information about the properties of the plasma generated by the spark discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Visualizing the Active Site Oxyanion Loop Transition Upon Ensitrelvir Binding and Transient Dimerization of SARS-CoV-2 Main Protease

N-terminal autoprocessing from its polyprotein precursor enables creating the mature-like stable dimer interface of SARS-CoV-2 main protease (MPro), concomitant with the active site oxyanion loop equilibrium transitioning to the active conformation (E*) and onset of catalytic activity. Here, through mutagenesis of critical interface residues and evaluating noncovalent inhibitor (ensitrelvir, ESV) facilitated dimerization through its binding to MPro, we demonstrate that residues extending from Ser1 through Glu14 are critical for dimerization. Combined mutations G11A, E290A and R298A (MPro™) restrict dimerization even upon binding of ESV to monomeric MPro™ with an inhibitor dissociation constant of 7.4 ± 1.6 µM. Contrasting the covalent inhibitor NMV or GC373 binding to monomeric MPro, ESV binding enabled capturing the transition of the oxyanion loop conformations in the absence of a reactive warhead and independent of dimerization. Characterization of complexes by room-temperature X-ray crystallography reveals ESV bound to the E* state of monomeric MPro as well as an intermediate approaching the inactive state (E). It appears that the E* to E equilibrium shift occurs initially from G138-F140 residues, leading to the unwinding of the loop and formation of the 3 10 -helix. Finally, we describe a transient dimer structure of the MPro precursor held together through interactions of residues A5-G11 with distinct states of the active sites, E and E*, likely representing an intermediate in the autoprocessing pathway.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Characteristics of the sub-micron ash aerosol generated during oxy-coal combustion at atmospheric and elevated pressures

This paper is concerned with the effect of pressure on the particle size distribution and the size-segregated composition of the sub-micron ash aerosol created during oxy-coal combustion under near practical self-sustaining combustion conditions. The problem is important because pressurized oxy-coal combustion has been proposed as one promising technology to minimize CO 2 emissions. Sub-micron ash plays a major role in ash deposition mechanisms, which, in turn, can control boiler performance. In this work, the same bituminous coal was burned at pressures of 1, 8 and 15 bar in O 2 /CO 2 environments. Here, tests employed a 100 kW (rated) oxy-fuel combustor (OFC) operated at atmospheric pressure (1 bar) and a 300 kW (rated) entrained-flow pressurized reactor (EFPR) at elevated pressures (8 and 15 bar). Although these tests were conducted under near practical combustion conditions, confounding effects of peak flame temperature variations were minimized for the 1 bar and 15 bar tests, allowing the role of elevated pressure to be isolated. For the EFPR tests, a specially designed sampling system was used to sample sub-micron ash aerosols from the pressurized combustor and is described in detail. Results showed that at the same peak temperature but higher pressure, the fractions of ash aerosol partitioned into the PM 0.6 and PM 1 size fractions were greatly diminished. Moreover, elevated pressures caused significant changes in the composition of the (size-segregated) sub-micron aerosol, especially in its alkali content, which increased significantly. Examination of fractions of each element that ended up in the sub-micron fume suggested that, at constant temperature, the effect of pressure on vaporization of semi-volatile metals was very different from that on the release into the gas phase of non-volatile metals and could not be explained by equilibrium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Rate Equations for Reversible Disproportionation Reactions and Fitting to Time-Course Data

Integrated rate equations are straightforward to fit to experimental data to verify a proposed mechanism and to extract kinetic parameters. Such equations are derived for reversible disproportionation/comproportionation reactions with any set of initial concentrations. Extraction of forward and reverse rate constants from experimental data by fitting the rate law to the data is demonstrated for the disproportionation of 2,2,6,6-tetramethyl-1-piperidinyl-N-oxyl (TEMPO) under acidic conditions where the approach to equilibrium is observed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Guest-dependent bond flexibility in UiO-66, a “stable” MOF

We report “flexibility constants”—a conceptual analog to metal–ligand stability constants—of UiO-66, the prototypical “stable” MOF, across a wide temperature range in both vacuum and in the presence of typical guest solvents. With these data, we extract key thermodynamic parameters governing the reversible bond equilibrium and demonstrate that guest molecules strongly favor the reversible dissociation of MOF metal–linker bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Selective electron-phonon coupling strength from nonequilibrium optical spectroscopy: The case of MgB 2

The coupling between quasiparticles and bosonic excitations rules the energy transfer pathways in condensed matter systems. The possibility of inferring the strength of specific coupling channels from their characteristic timescales measured in nonequilibrium experiments is still an open question. Here, in this study, we investigate MgB 2 , in which conventional superconductivity at temperatures as high as 39 K is mediated by the strong coupling between the conduction electrons and the 𝐸 2⁢𝑔 phonon mode. By means of broadband time-resolved optical spectroscopy, we show that this selective electron-phonon coupling dictates the nonequilibrium optical response of MgB 2 at early times (<100 fs) after photoexcitation. Furthermore, based on an effective temperature model analysis, we estimate its contribution to the total electron-boson coupling function extracted from complementary equilibrium spectroscopy approaches, namely, optical reflectivity and angle-resolved photoemission spectroscopy. The coupling strength with the 𝐸 2⁢𝑔 phonon modes is thus estimated to be 𝜆 ≃ 0.56, which is approximately half of the total coupling constant, in agreement with 𝑎𝑏 initio calculations from the literature. As a benchmark, broadband time-resolved optical spectroscopy is performed also on the isostructural and nonsuperconducting compound AlB 2 , showing that the nonequilibrium optical response relaxes on a slower timescale due to the lack of strongly coupled phonon modes. Our findings demonstrate the possibility to resolve and quantify selective electron-phonon coupling from nonequilibrium optical spectroscopy.

Mor, Selene [Università Cattolica del Sacro Cuore,↗
DOE OSTI2021

Iron, magnesium, and titanium isotopic fractionations between garnet, ilmenite, f ayalite, biotite, and tourmaline: Results from NRIXS, ab initio , and study of mineral separates from the Moosilauke metapelite

Interpreting isotopic signatures documented in natural rocks requires knowledge of equilibrium isotopic fractionation factors. Here, we determine equilibrium Fe isotope fractionation factors between several common rock-forming minerals using a comparative approach involving three independent methods: (i) isotopic analyses of natural minerals from a metapelite from Mt. Moosilauke, New Hampshire, for which equilibration temperature and pressure are well constrained to be near the aluminosilicate triple point (T ≃ 500 °C, P ≃ 4 kbar), (ii) Nuclear Resonant Inelastic X-ray Scattering (NRIXS) measurements of Fe force constants of minerals, and (iii) Density Functional Theory (DFT) ab initio calculations of Fe force constants of minerals. The minerals studied for Fe isotopes include, in increasing order of their β-factors: garnet < ilmenite ≈ fayalite < biotite < tourmaline < muscovite ≈ plagioclase. Some of this ordering is affected by the presence of Fe 3+ in the minerals, which tends to form stiffer bonds and be associated with heavy Fe isotope enrichments relative to Fe 2+ . We are, however, able to assess the magnitude of the effect of the ratio Fe 3+ /ΣFe on equilibrium fractionation factors, notably on the ilmenite-hematite solid solution. Equilibrium Fe isotopic fractionation factors between garnet, ilmenite, biotite, tourmaline and fayalite are determined. We also report Mg and Ti isotopic compositions of selected Moosilauke minerals that allow us to better constrain the equilibrium fractionation factors for garnet-biotite-tourmaline (Mg isotopes) and biotite-ilmenite (Ti isotopes). We show how the newly determined equilibrium fractionation factors can be used to address diverse problems in Earth and planetary sciences, notably (i) Fe and Mg isotopic fractionation during anatexis, (ii) Fe isotopic fractionation in lunar ilmenite, and (iii) Ti isotopic fractionation during fluvial transport of minerals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Master equation study of hydrogen abstraction from HCHO by OH via a chemically activated intermediate

The abstraction reaction of hydrogen from formaldehyde by OH radical plays an important role in formaldehyde oxidation. The reaction involves a bimolecular association to form a chemically activated hydrogen-bonded reaction complex followed by a unimolecular reaction of the complex to generate the products. The reaction rate is usually considered to be pressure-independent by assuming equilibrium between the reactants and the complex. However, our nonequilibrium calculations based on the chemically significant eigenmode of the master equation, carried out with our recently developed TUMME program, indicate that the reaction complex makes the rate constant dependent on pressure at low temperatures (T < 200 K). The calculations include anharmonicity, variational effects, and multi-dimensional tunneling. We find that the reaction rate constant reaches a low-pressure limit at pressures below 10 Torr over the whole investigated temperature range (20–1800 K), which explains why the available low-temperature experiments, which are for pressures below 2 Torr, did not observe the pressure dependence. A new extension of the TUMME master-equation program is used to explore the time evolutions of the concentrations of the OH radical and the complex under pseudo-first-order conditions. The time-dependent evolution of the concentrations of the complex at a low temperature provide direct evidence for the stabilization of the reaction complex at high pressures, and it shows the negligible role of the stabilized reaction complex at low pressures. Finally, the picture that emerges is qualitatively consistent with our previous study of the reaction of methanol with OH in that the tunneling in the unimolecular step from the complex to the products affects the phenomenological reaction rate constants differently at high and low pressures and leads to a significant pressure effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗