Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reaction Rates”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

A kinetic model for the thermal nitridation of SiO2/Si

To explain the observed nitrogen distributions in thermally nitridated SiO2 films, a kinetic model is proposed in which the nitridation process is simulated, using the first-order chemical kinetics and Arrhenius dependence of the diffusion and reaction rates on temperature. The calculations show that initially, as the substrate reacts with diffusing nitrogen, a nitrogen-rich oxynitride forms at the SiO2-Si interface, while at nitridation temperatures above 1000 C, an oxygen-rich oxynitride subsequently forms at the interface, due to reaction of the substrate with an increasingly concentrated oxygen displaced by the slower nitridation of the SiO2. This sequence of events results in a nitrogen distribution in which the peak of the interfacial nitrogen concentration occurs away from the interface. The results are compared with the observed nitrogen distribution. The calculated results have correctly predicted the positions of the interfacial nitrogen peaks at the temperatures of 800, 1000, and 1150 C. To account for the observed width of the interfacial nitrogen distribution, it was found necessary to include in the simulations the effect of interfacial strain.

Vasquez, R. P.↗
NASA NTRS1992

Solar control of the upper atmosphere of Triton

If the upper atmosphere and ionosphere of Triton are controlled by precipitation of electrons from Neptune's magnetosphere as previously proposed, Triton could have the only ionosphere in the solar system not controlled by solar radiation. However, a new model of Triton's atmosphere, in which only solar radiation is present, predicts a large column of carbon atoms. With an assumed, but reasonable, rate of charge transfer between N2(+) and C, a peak C(+) abundance results that is close to the peak electron densities measured by Voyager in Triton's ionosphere. These results suggest that Triton's upper atmospheric chemistry may thus be solar-controlled. Measurement of key reaction rate constants, currently unknown or highly uncertain at Triton's low temperatures, would help to clarify the chemical and physical processes occurring in Triton's atmosphere.

Lyons, James R.↗
NASA NTRS1978

Absolute rate of the reaction of bromine atoms with ozone from 200 to 360 K

The rate constant for the reaction Br + O3 yields BrO + O2 has been measured from 200 to 360 K by the technique of flash photolysis coupled to time resolved detection of bromine atoms by resonance fluorescence (FP-RF). Br atoms were produced by the flash photolysis of CH3Br at a wavelength of 165 nm. O3 concentration was monitored continuously under reaction conditions by absorption at 253.7 nm. At each of five temperatures the results were independent of substantial variations in O3 concentration, total pressure (Ar), and limited variations in flash intensity (i.e., initial Br concentration). The measured rate constants obey the Arrhenius expression, k = (7.74 plus or minus 0.50) x 10 to the -12th exp(-603 plus or minus 16/T) cu cm/molecule/sec, where the error quoted is two standard deviations.

Michael, J. V.↗
DOE OSTI2023

An ab initio free energy study of the reaction mechanism and rate-limiting steps of CO 2 capture by aqueous glycine

Aqueous amino acids are promising absorbents for direct air capture (DAC) of CO 2 . Herein, we investigate the possibility of kinetic control of CO 2 absorption with aqueous anionic glycine (GLY - ) by employing extensive ab initio molecular dynamics simulations, free energy analysis, and reaction rate theory. We find that first GLY - binds to CO 2 by overcoming a barrier (7.4 kcal/mol) to form a zwitterion intermediate, which then releases a proton by overcoming a similar barrier. Despite the similarity in the barrier, zwitterion formation appears to be the rate-limiting step because it is two orders of magnitude slower (microseconds) than the proton release step. This is predominantly due to stronger nonequilibrium solvent effects for the former that cause many barrier-recrossing events and effectively slow down the reaction rate. Such a detailed fundamental understanding of the amino acid-based CO 2 -absorption mechanism and rates is key to improving the kinetic efficiency of DAC technology.

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

Striations in moderate pressure dc driven nitrogen glow discharge

Plasma stratification has been studied for more than a century. Despite the many experimental studies reported on this topic, theoretical analyses and numerical modeling of this phenomenon have been mostly limited to rare gases. In this work, a one-dimensional fluid model with detailed kinetics of electrons and vibrationally excited molecules is employed to simulate moderate-pressure (i.e. a few Torrs) dc discharge in nitrogen in a 15.5 cm long tube of radius 0.55 cm. The model also considers ambipolar diffusion to account for the radial loss of ions and electrons to the wall. The proposed model predicts self-excited standing striations in nitrogen for a range of discharge currents. The impact of electron transport parameters and reaction rates obtained from a solution of local two-term and a multi-term Boltzmann equation on the predictions are assessed. In-depth kinetic analysis indicates that the striations result from the undulations in electron temperature caused due to the interaction between ionization and vibrational reactions. Furthermore, the vibrationally excited molecules associated with the lower energy levels are found to influence nitrogen plasma stratification and the striation pattern strongly. A balance between ionization processes and electron energy transport allows the formation of the observed standing striations. Simulations were conducted for a range of discharge current densities from ~0.018 to 0.080 mA cm –2 , for an operating pressure of 0.7 Torr. Parametric studies show that the striation length decreases with increasing discharge current. Furthermore, the predictions from the model are compared against experimental measurements and are found to agree favorably.

42 ENGINEERING↗
NASA NTRS1990

Tropospheric OH in a three-dimensional chemical tracer model - An assessment based on observations of CH3CCl3

The paper introduces the representation of tropospheric chemistry allowing for calculation of chemical reaction rates within a chemical tracer model (CTM) simulation with the accuracy comparable to that of a detailed chemical model. The global distribution of tropospheric OH is calculated as a function of time, atmospheric composition, temperature, and solar irradiance. The globally integrated OH concentration weighted by frequency of reaction with CH3CCl3 is evaluated, along with the feasibility of using the observed latitudinal distribution of CH3CCl3 as a test for computed concentrations of OH. Seasonal variations of CH3CCl3 are also considered as a potential test for computed concentrations of OH. It is demonstrated that annual cycles of CH3CCl3 are controlled, or strongly influenced, by dynamical processes at latitudes northward of 25 deg South.

Spivakovsky, C. M.↗
DOE OSTI2024

A Cellular Automaton Simulation for Predicting Phase Evolution in Solid-State Reactions

New computational tools for solid-state synthesis recipe design are needed in order to accelerate the experimental realization of novel functional materials proposed by high-throughput materials discovery workflows. This work contributes a cellular automaton simulation framework for predicting the time-dependent evolution of intermediate and product phases during solid-state reactions as a function of precursor choice and amount, reaction atmosphere, and heating profile. The simulation captures the effects of reactant particle spatial distribution, particle melting, and reaction atmosphere. Reaction rates based on rudimentary kinetics are estimated using density functional theory data from the Materials Project and machine learning estimators for the melting point and the vibrational entropy component of the Gibbs free energy. The resulting simulation framework allows for the prediction of the likely outcome of a reaction recipe before any experiments are performed. We analyze five experimental solid-state recipes for BaTiO 3 , CaZrN 2 , and YMnO 3 found in the literature to illustrate the performance of the model in capturing reaction selectivity and reaction pathways as a function of temperature and precursor choice. This simulation framework offers an easier way to optimize existing recipes, aid in the identification of intermediates, and design effective recipes for yet unrealized inorganic solids in silico .

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

Electrochemical characterization of the plasma-water interface

Nonthermal atmospheric pressure plasmas generated using noble gases act as gaseous electrodes capable of promoting electrochemical reactions when in contact with aqueous media. Despite overwhelming evidence in favor of the interfacial region between the plasma and solution being a reducing environment, the plasma-liquid interface has not been evaluated from an electrochemical perspective. Herein, we present an electrochemical method to measure the overpotential at the plasma-liquid interface of solutions containing organic redox couples. Negative overpotentials indicate that the interfacial region is a reducing environment, which is consistent with the results of colorimetric assays using redox indicators. The reaction rates are related to the magnitude of the overpotential, as expected from electrochemical kinetic models. The ability to characterize the overpotential of plasma-liquid interfaces adds to the plasma electrochemistry toolkit and contributes to efforts being made to tune gaseous electrodes for the rational control of redox reactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

BioRT‐HBV 1.0: A Biogeochemical Reactive Transport Model at the Watershed Scale

Abstract Reactive Transport Models (RTMs) are essential tools for understanding and predicting intertwined ecohydrological and biogeochemical processes on land and in rivers. While traditional RTMs have focused primarily on subsurface processes, recent watershed‐scale RTMs have integrated ecohydrological and biogeochemical interactions between surface and subsurface. These emergent, watershed‐scale RTMs are often spatially explicit and require extensive data, computational power, and computational expertise. There is however a pressing need to create parsimonious models that require minimal data and are accessible to scientists with limited computational background. To that end, we have developed BioRT‐HBV 1.0, a watershed‐scale, hydro‐biogeochemical RTM that builds upon the widely used, bucket‐type HBV model known for its simplicity and minimal data requirements. BioRT‐HBV uses the conceptual structure and hydrology output of HBV to simulate processes including advective solute transport and biogeochemical reactions that depend on reaction thermodynamics and kinetics. These reactions include, for example, chemical weathering, soil respiration, and nutrient transformation. The model uses time series of weather (air temperature, precipitation, and potential evapotranspiration) and initial biogeochemical conditions of subsurface water, soils, and rocks as input, and output times series of reaction rates and solute concentrations in subsurface waters and rivers. This paper presents the model structure and governing equations and demonstrates its utility with examples simulating carbon and nitrogen processes in a headwater catchment. As shown in the examples, BioRT‐HBV can be used to illuminate the dynamics of biogeochemical reactions in the invisible, arduous‐to‐measure subsurface, and their influence on the observed stream or river chemistry and solute export. With its parsimonious structure and easy‐to‐use graphical user interface, BioRT‐HBV can be a useful research tool for users without in‐depth computational training. It can additionally serve as an educational tool that promotes pollination of ideas across disciplines and foster a diverse, equal, and inclusive user community.

Sadayappan, Kayalvizhi↗
NASA NTRS1990

In situ analysis of CO during chemisorption and oxidation on graphite: Supported Pt by FTIR-microspectrometry

For chemisorption and oxidation on Pt/HOPG (highly-orientated pyrolytic) graphite, reflectance Fourier Transform Infrared (FTIR)-microspectrometry reveals a variable state and reactivity for CO. Even for model surface science systems, where surface heterogeneity is minimal, surface diffusion may be too slow relative to the reaction rate to avoid segregation of reactants into surface islands under steady-state conditions. Thus in CO oxidation on Pt (where the relevant surface diffusion coefficients are such that D sub O less than D sub CO) then reactant CO islands exists at the perimeters of which the surface reaction is thought to occur. Furthermore CO can chemisorb on metals in linear and bridge forms to extents which vary with precise faces predominantly exposed coverage, etc. Infrared has long been used to probe the nature of adsorbed CO on model film and heterogeneous surfaces, but it may now be that FTIR-microspectrometry will allow the state of this adsorbate and reactant to be investigated with a spatial resolution of 4.4 microns on model (and real) catalytic surfaces.

Self, Valerie A.↗
NASA NTRS1996

Combustion Processes in Hybrid Rocket Engines

In recent years, there has been a resurgence of interest in the development of hybrid rocket engines for advanced launch vehicle applications. Hybrid propulsion systems use a solid fuel such as hydroxyl-terminated polybutadiene (HTPB) along with a gaseous/liquid oxidizer. The performance of hybrid combustors depends on the convective and radiative heat fluxes to the fuel surface, the rate of pyrolysis in the solid phase, and the turbulent combustion processes in the gaseous phases. These processes in combination specify the regression rates of the fuel surface and thereby the utilization efficiency of the fuel. In this paper, we employ computational fluid dynamics (CFD) techniques in order to gain a quantitative understanding of the physical trends in hybrid rocket combustors. The computational modeling is tailored to ongoing experiments at Penn State that employ a two dimensional slab burner configuration. The coordinated computational/experimental effort enables model validation while providing an understanding of the experimental observations. Computations to date have included the full length geometry with and with the aft nozzle section as well as shorter length domains for extensive parametric characterization. HTPB is sed as the fuel with 1,3 butadiene being taken as the gaseous product of the pyrolysis. Pure gaseous oxygen is taken as the oxidizer. The fuel regression rate is specified using an Arrhenius rate reaction, which the fuel surface temperature is given by an energy balance involving gas-phase convection and radiation as well as thermal conduction in the solid-phase. For the gas-phase combustion, a two step global reaction is used. The standard kappa - epsilon model is used for turbulence closure. Radiation is presently treated using a simple diffusion approximation which is valid for large optical path lengths, representative of radiation from soot particles. Computational results are obtained to determine the trends in the fuel burning or regression rates as a function of the head-end oxidizer mass flux, G=rho(e)U(e), and the chamber pressure. Furthermore, computation of the full slab burner configuration has also been obtained for various stages of the burn. Comparisons with available experimental data from small scale tests conducted by General Dynamics-Thiokol-Rocketdyne suggest reasonable agreement in the predicted regression rates. Future work will include: (1) a model for soot generation in the flame for more quantitative radiative transfer modelling, (2) a parametric study of combustion efficiency, and (3) transient calculations to help determine the possible mechanisms responsible for combustion instability in hybrid rocket motors.

Venkateswaran,S.↗
DOE OSTI2024

Verification and validation of detonation-shock-dynamics relations for explosives described by general equation of state and chemical reaction models

Detonation shock dynamics is a powerful method to model the behaviour of High Explosives (HE). However in order to use this method, the underlying relationship between the local radius of curvature and the detonation speed must be known. Previous work has developed methods to calculate this effect using simple, single-step Arrhenius and polytropic gas, models for the chemical reaction and the equation of state, respectively. In recent years, more complex models for both reaction rates and equations of state have been developed which show better agreement with experimental data than these simple models, especially when considering condensed phase explosives.. This work presents the governing equations for solving these problems in a way that is generalised to use arbitrary equations of state as well as reaction models which may have more than a single step and multiple product species. This implementation is verified against exact solutions, demonstrating that the equations were implemented properly. The verified algorithm is then validated against experimental data and high fidelity simulations, showing that it is able to make accurate predictions in a regime where the underlying assumptions of the governing equations are valid. Importantly, this approach has many applications: from creating equivalent detonation shock dynamics models for existing reactive burn calibrations for HE; to developing new functional forms and calibrations of reactive burn models for condensed phase high explosives.

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

Degradation of Poly- and Perfluoroalkyl Substances (PFAS) in Water via High Power, Energy-Efficient Electron Beam Accelerator

The goal of the 2-year workplan was to see if electron beam (EB) could be used to break down a sub-set of the larger chemical family of per and polyfluoroalkylated substances (PFAS) in an energy efficient and economical manner when compared to conventional water treatment technologies. Year one (Y1) work focused on sample EB treatment work in the Fermi National Accelerator Laboratory’s (FNALs) Accelerator Applications Demonstration and Development (A2D2) EB accelerator. While there are reportedly thousands of types of PFAS, for the point of most of the work herein, a small subset was examined, typically perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA). PFOA and PFOS are two of the most well studied PFAS and are studied for baseline evaluations and are considered most useful. The work from Y1 provided information about the optimal operating parameters and additives to use when treating PFOS and PFOA via EB. The data were then used to see where in a water treatment system an EB accelerator would be best suited to treat PFAS. A conventional water treatment technology, GAC, was then compared to e-beam treatment technology with respect to energy and costs for treatment. In year two (Y2), several conventional e-beam accelerator designs, and FNAL’s developmental compact SRF accelerator design, were evaluated for their suitability in PFAS treatment, from an energy efficiency and cost standpoint. Several EB parameters were evaluated and optimized for the removal of PFOA and PFOS from water at normal pressure and temperature, measured as total PFAS removal. Under the optimized test conditions both PFOA showed complete destruction to inorganic fluoride, and PFOS to inorganic fluoride and sulfate, with mass balance. The effect on PFAS removal relative to solution pH, total EB dose, EB dose rate, dissolved oxygen concentration (DO), temperature, and initial PFAS concentration were evaluated. In general, PFOA was easier to destroy than PFOS. Degradation products, typically observed under less-than-optimal EB conditions, provided insight to degradation mechanisms. Products were identified to rule out possible deleterious biproduct formation. The water radiolysis radical reaction kinetics with PFOS and PFOA were not dependent on the initial concentration over 5-orders of magnitude from 2 μg/L to 20 mg/L. This is thought to be because there was an overabundance of the reactive water radiolysis radicals relative to PFAS molecules and largely attributed to aqueous electrons. The reaction rates appeared to be diffusion limited. Testing at higher concentrations (100-200 mg/L) showed a decrease in removal efficiency, suggesting alternative kinetics, possibly second order rates, at higher concentrations. In all, we successfully defined a set of optimal EB parameters to treat PFOA and PFOS at concentrations of 20 mg/L in water with destruction efficiencies near 100%. We further tested the optimized EB parameters with other types of PFAS, including shorter and longer fluorocarbon chain homologs of PFOA and PFOS, and PFAS with alternative functional groups such as sulfonamides. Based on our results EB can be optimized as an effective destructive technology for removing PFAS from water. The conditions optimized for PFOA and PFOS were less effective with ultra-short fluorocarbon compounds like TFMS, PFES, PFPS and PFBS, and likely require re-optimization of parameters to them. In all, it was determined that from a cost and energy efficiency standpoint, EB would be best applied to waste streams with relatively high concentrations of PFOS and PFOA and is not as cost effective as GAC treatment for removing low concentrations of PFAS from water. Higher concentrations of PFAS can be found in the wastewater of conventional treatment processes such as RO and IE and therefore EB may be used to supplement such treatment technologies. Some real-world IE regeneration wash water and RO reject water containing higher concentrations of PFAS and obtained from pilot scale industrial wastewater treatment system at a fluorochemical manufacturing facility, showed that EB could remove PFAS from such types of wastewaters. The IE regenerant wash water appeared to be the most efficient of the two types of wastewaters tested. However, some further optimization of the EB parameters for the specific PFAS types present in those wastewaters may be required. Also, the effects of co-present TOC and mineral salts should be considered during such optimization efforts. From the experimental Y1 results it was seen that the aqueous electron drives degradation of the PFAS. In a hypothetical water treatment skid using EB for PFAS destruction the parameters of the system should be optimized to promote aqueous electron production. Before EB treatment, the PFAS should be preconcentrated when possible, the pH should be raised to pH 10 or higher to enhance aqueous electron production, and the water should be nitrogen purged to remove dissolved oxygen to minimize aqueous electron scavenging. An excel spreadsheet was created that calculates optimal conditions based on inlet PFAS concentration and desired outlet concentration, by optimizing the accelerator power, dose rate, water treatment rate, pH and dissolved oxygen levels to reach the desired endpoint. Given this information on accelerator operating conditions five different EB accelerator systems were compared. One EB system was a continuous-wave, linear superconducting accelerator being designed at Fermilab. Three other EB systems (IMPELA at 5% and 25% duty factor and the ILU-14) were normal conducting pulsed linear accelerators. The fifth system was an IBA Rhodotron which is a normal conducting, circular, continuous-wave accelerator. The accelerator efficiency (% of the incoming power that is used in water treatment) was the dominating factor in accelerator choice. The radio frequency (RF) power supply and the accelerator design (superconducting versus warm technology) drive the accelerator efficiency. The IBA Rhodotron was seen to be the most energy efficient commercially available technology with a wall-plug (total) power efficiency of 43% at 400 kW. The Fermilab design, with a prototype for a different application currently being fabricated, was the most energy efficient at 55% when driven by a Klystron RF power supply and as high as 77% when powered by a magnetron. As the Fermilab design was the most energy efficient by approximately 10-30%, further design work was done on the accelerator and beam delivery system specific to the destruction of PFAS in water. The Fermilab design is unique from industrial accelerators in that is superconducting. Superconducting technology allows for the acceleration of electrons without losses. The accelerator must be cooled to below the point where it is superconducting and is operated around 4 degrees Kelvin. The bulk of the design work for the accelerator is on making the accelerator as energy efficient as possible so that it does not require liquid helium and can be cooled with conduction cooling via cryocoolers. Final design work resulted in an EB accelerator that would operate at minimally 200 kW and 10 MeV. Prototype construction would cost $\$ $7.8 million dollars when driven by a Klystron power supply. A second version of the same accelerator would cost $\$ $5.5 million dollars when driven by a magnetron that is still under development. The commercially available 300 kW IBA Rhodotron cost was estimated at approximately $\$ $9 million. While it is hard to directly compare, an operational GAC system used by 3M for groundwater treatment capital cost (2022 dollars) was estimated to cost $\$ $3.3 million. While the capital expense of the EB accelerator systems was higher than GAC, the accelerator EB treatment would result in destruction of the PFAS and not just sequestration of PFAS to form a new waste stream that requires further treatment or disposal. The operating cost to destroy the PFAS via 400 kw EB system was less than $\$ $1000/kg of PFAS destroyed when treating at a 20 mg/L PFAS concentration, compared to GAC with operating costs that calculated at $\$ $27,530 per kg of PFAS sequestered when treating 100 μg/L PFOA and PFOS combined concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1989

Two-dimensional analysis of two-phase reacting flow in a firing direct-injection diesel engine

The flow field, spray penetration, and combustion in two-stroke diesel engines are described. Fuel injection begins at 345 degrees after top dead center (ATDC) and n-dodecane is used as the liquid fuel. Arrhenius kinetics is used to calculate the reaction rate term in the quasi-global combustion model. When the temperature, fuel, and oxygen mass fraction are within suitable flammability limits, combustion begins spontaneously. No spark is necessary to ignite a localized high temperature region. Compression is sufficient to increase the gaseous phase temperature to a point where spontaneous chemical reactions occur. Results are described for a swirl angle of 22.5 degrees.

Nguyen, H. Lee↗
DOE OSTI2023

Probing reaction channels via reinforcement learning

Abstract Chemical reactions are dynamical processes involving the correlated reorganization of atomic configurations, driving the conversion of an initial reactant into a result product. By virtue of the metastability of both the reactants and products, chemical reactions are rare events, proceeding fleetingly. Reaction pathways can be modelled probabilistically by using the notion of reactive density in the phase space of the molecular system. Such density is related to a function known as the committor function, which describes the likelihood of a configuration evolving to one of the nearby metastable regions. In theory, the committor function can be obtained by solving the backward Kolmogorov equation (BKE), which is a partial differential equation (PDE) defined in the full dimensional phase space. However, using traditional methods to solve this problem is not practical for high dimensional systems. In this work, we propose a reinforcement learning based method to identify important configurations that connect reactant and product states along chemical reaction paths. By shooting multiple trajectories from these configurations, we can generate an ensemble of states that concentrate on the transition path ensemble. This configuration ensemble can be effectively employed in a neural network-based PDE solver to obtain an approximation solution of a restricted BKE, even when the dimension of the problem is very high. The resulting solution provides an approximation for the committor function that encodes mechanistic information for the reaction, paving a new way for understanding of complex chemical reactions and evaluation of reaction rates.

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

A Pd III Sulfate Dimer Initiates Rapid Methane Monofunctionalization by H Atom Abstraction

An electrogenerated Pd III 2 species in fuming sulfuric acid is competent for rapid and concurrent methane monohydroxylation to methyl bisulfate (CH 3 OSO 3 H) and methane sulfonation to methanesulfonic acid (CH 3 SO 3 H). In situ NMR at 50 °C is used to track methane transformation exclusively to CH 3 OSO 3 H and CH 3 SO 3 H at high conversions. Integrating a set of kinetic and computational studies, the mechanism of methane monofunctionalization by Pd III 2 is examined. Here, experimental rate laws and common kinetic isotope effects for CH 3 OSO 3 H and CH 3 SO 3 H formation suggest that both transformations proceed via a common rate-limiting C-H activation step. Introduction of O 2 or Pd II,III 2 suppresses CH 3 SO 3 H generation, indicating a radical chain sequence. Although the metal-metal bonded Pd III 2 complex is a net two-electron oxidant, our aggregate kinetic data point to a mechanistic model that features rate-limiting H atom abstraction by the Pd III 2 complex to generate a methyl radical intermediate. The CH 3 • intermediate then recombines with Pd II,III 2 to furnish a CH 3 Pd III 2 intermediate that reductively eliminates CH 3 OSO 3 H. Alternatively, the CH 3 intermediate can enter a chain reaction with SO 3 to generate CH 3 SO 3 H. DFT computations support the radical-based C-H activation by Pd III 2 and delineate H atom abstraction pathways with computed reaction barriers and kinetic isotope effects (KIEs) that are consistent with experimental data. These mechanistic investigations challenge the paradigm of electrophilic C-H activation and highlight H atom abstraction as a potent pathway for selective methane C-H oxidative functionalization at high reaction rates.

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

Transient cycling of nitrogen, organic carbon and oxygen within the free-flowing Columbia River corridor: Linking exposure time dependent biogeochemical reactions to river stage fluctuations (Final project report)

The objective of this project was to develop next-generation techniques for representing the transformations of complex reaction networks like those impacting transient river-corridors, then apply them to develop upscaling tools. The project contributed significantly to the development of Lagrangian “mass transfer particle tracking” (MTPT) tools that are the first numerical framework that explicitly separates mixing and spreading processes; a crucial distinction for accurate representations of reaction rates across scales. The MTPT approach uses an operator splitting scheme where physical transport processes are simulated using classical random walk methods and each particle is treated as a “container” that carries with it any number of chemical species. Mixing between particles is simulated using a colocation-probability based mass transfer kernel and reactions are evaluated on each particle after mixing. The method is stable, accurate, and also allows the explicit modeling of residence time distributions of the individual containers of mass. The MTPT scheme allowed us to explore several research questions related to reactive transport in river corridors and complex reaction networks. We found that it is possible to quantify the exposure time distributions (ETDs) of reactants and that this does suggest some pathways to upscaling. However, comparable approximations can be obtained using simplified (reduced dimensionality) MTPT simulations in less time with similar uncertainty, so development of extensive ETD-based methods was not productive. The major findings are that i) the MTPT schemes are robust and highly accurate across a range of arbitrarily complex reactions, ii) the methods can be efficiently parallelized and the parallel performance characteristics of MTPT are predictable, and iii) the MTPT tools allow simultaneous tracking of residence time. The development of these research tools into efficient software packages continues. The methods are already available to the community because open-source, working examples have been included with all publications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Microscopic Investigation of H 2 Reduced CuO x /Cu(111) and ZnO/CuO x /Cu(111) Inverse Catalysts: STM, AP-XPS, and DFT Studies

Understanding the reduction mechanism of ZnO/CuO x interfaces by hydrogen is of great importance for advancing the performance of industrial catalysts for CO 2 hydrogenation to methanol. Here, the reduction of pristine and ZnO-modified CuO x /Cu(111) by H 2 was investigated using ambient pressure scanning tunnelling microscopy (AP-STM), ambient pressure X-ray photoelectron spectroscopy (AP-XPS) and density functional theory (DFT). The morphological changes and reaction rates seen for the reduction of CuO x /Cu(111) and ZnO/CuO x /Cu(111) are very different. On CuO x /Cu(111), perfect "44" and "29" structures displayed a very low reactivity towards H 2 at room temperature. A long induction period associated with an autocatalytic process was observed to enable the reduction by the removal of chemisorbed non-lattice oxygen initially and lattice oxygen sequentially at the CuO x -Cu interface, which led to formation of oxygen deficient "5-7" hex and honeycomb structures. In the final stages of the reduction process, regions of residual oxygen species and metallic Cu were seen. The addition of ZnO particles to CuO x /Cu(111) opened new reaction channels. On the ZnO sites, the dissociation of H 2 was fast and H adatoms easily migrated to adjacent regions of copper oxide. This hydrogen spillover substantially enhanced the rate of oxygen removal, resulting in the rapid reduction of the copper oxide located in the periphery of the zinc oxide islands with no signs for the reduction of ZnO. The deposited ZnO completely modified the dynamics for H 2 dissociation and hydrogen migration, providing an excellent source for CO 2 hydrogenation processes on the inverse oxide/metal system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗