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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Lifetime measurements of excited states in 15 O

The CNO cycle is the main energy source in stars more massive than our Sun, it defines the energy production and the cycle time that lead to the lifetime of massive stars, and it is an important tool for the determination of the age of globular clusters. One of the largest uncertainties in the CNO chain of reactions comes from the uncertainty in the 14 N(p,γ) 15 O reaction rate. This uncertainty arises predominantly from the uncertainty in the lifetime of the subthreshold state in 15 O at E x = 6792 keV. Previous measurements of this state's lifetime are significantly discrepant. Here, we report on a new lifetime measurement of this state, as well as the excited states in 15 O at E x = 5181 keV and E x = 6172 keV, populated via the 14 N(p,γ) 15 O reaction at proton energies of E p = 1020 keV and E p = 1570 keV. Furthermore, the lifetimes have been determined with the Doppler-Shift Attenuation Method (DSAM) with three separate, nitrogen-implanted targets with Mo, Ta, and W backing. We obtained lifetimes from the weighted average of the three measurements, allowing us to account for systematic differences between the backing materials. For the 6792 keV state, we obtained a τ = 0.6 ± 0.4 fs. To provide cross-validation of our method, we measured the known lifetimes of the states at 5181 and 6172 keV to be τ = 7.5 ± 3.0 and τ = 0.7 ± 0.5 fs, respectively, in good agreement with previous measurements.

6 ≤ A ≤ 19↗
DOE OSTI2020

Material erosion measurements and expected operational lifetime of a deployable photon sieve payload

Spacecraft operating in low-Earth orbit are subjected to a number of hazardous environmental constituents that can lead to decreased system performance and reduced operational lifetimes. Due to their thermal, optical, and mechanical properties, polymers are used extensively in space systems; however they are particularly susceptible to material erosion and degradation as a result of exposure to the LEO environment. The focus of this research is to examine the material erosion and mass loss experienced by a custom Kaptonlike polyimide due to exposure in a simulated low-earth orbit environment. The deployable membrane telescope design discussed in this research is named Peregrine and is the scientific payload for FalconSat-7, a 3U cubesat designed and developed at the United States Air Force Academy (USAFA) for the purpose of demonstrating the capability of deployable membrane telescope technology on a nanosatellite platform. In addition to the polymer samples, chrome, silver and gold specimens will be examined to measure the oxidation rate and act as a control specimen, respectively. A magnetically filtered atomic oxygen plasma source has previously been developed and characterized for the purpose of simulating the low-Earth orbit environment. The plasma source can be operated at a variety of discharge currents and gas flow rates, of which the plasma parameters downstream of the source are dependent. The characteristics of the generated plasma were examined as a function of these operating parameters to optimize the production of O + ions with energy relevant to LEO applications. The erosion yield of the Kapton-like polyimide was experimentally determined to be 2.84 x 10 -24 cm 3 per atom which is in close agreement when compared to the on-orbit measurement for reaction efficiency of Kapton HN. We report the experimentally determined reaction rate for the Kapton-like polyimide was used to estimate the operational lifetime of the photon sieve during the solar conditions expected beginning in May 2019. The effective lifetime is estimated between 126 and 153 days.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Weathering Incongruence in Mountainous Mediterranean Climates Recorded by Stream Lithium Isotope Ratios

Lithium isotope ratios (δ 7 Li) of rivers are increasingly serving as a diagnostic of the balance between chemical and physical weathering contributions to overall landscape denudation rates. Here, we show that intermediate weathering intensities and highly enriched stream δ 7 Li values typically associated with lowland floodplains can also describe small upland watersheds subject to cool, wet climates. This behavior is revealed by stream δ 7 Li between +22.4 and +23.5‰ within a Critical Zone observatory located in the Cévennes region of southern France, where dilute stream solute concentrations and significant atmospheric deposition otherwise mask evidence of incongruence. The water-rock reaction pathways underlying this behavior are quantified through a multicomponent, isotope-enabled reactive transport model. Using geochemical characterization of soil profiles, bedrock, and long-term stream samples as constraints, we evolve the simulation from an initially unweathered granite to a steady state weathering profile which reflects the balance between (a) fluid infiltration and drainage and (b) bedrock uplift and soil erosion. Enriched stream δ 7 Li occurs because Li is strongly incorporated into actively precipitating secondary clay phases beyond what prior laboratory experiments have suggested. Chemical weathering incongruence is maintained despite relatively slow reaction rates and moderate clay accumulation due to a combination of two factors. First, reactive primary mineral phases persist across the weathering profile and effectively “shield” the secondary clays from resolubilization due to their greater solubility. Second, the clays accumulating in the near-surface profile are relatively mature weathering byproducts. These factors promote characteristically low total dissolved solute export from the catchment despite significant input of exogenous dust.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Biocatalytic Membranes for Carbon Capture and Utilization

Innovative carbon capture technologies that capture CO 2 from large point sources and directly from air are urgently needed to combat the climate crisis. Likewise, corresponding technologies are needed to convert this captured CO 2 into valuable chemical feedstocks and products that replace current fossil-based materials to close the loop in creating viable pathways for a renewable economy. Biocatalytic membranes that combine high reaction rates and enzyme selectivity with modularity, scalability, and membrane compactness show promise for both CO 2 capture and utilization. This review presents a systematic examination of technologies under development for CO 2 capture and utilization that employ both enzymes and membranes. CO 2 capture membranes are categorized by their mode of action as CO 2 separation membranes, including mixed matrix membranes (MMM) and liquid membranes (LM), or as CO 2 gas–liquid membrane contactors (GLMC). Because they selectively catalyze molecular reactions involving CO 2 , the two main classes of enzymes used for enhancing membrane function are carbonic anhydrase (CA) and formate dehydrogenase (FDH). Small organic molecules designed to mimic CA enzyme active sites are also being developed. CO 2 conversion membranes are described according to membrane functionality, the location of enzymes relative to the membrane, which includes different immobilization strategies, and regeneration methods for cofactors. Parameters crucial for the performance of these hybrid systems are discussed with tabulated examples. Progress and challenges are discussed, and perspectives on future research directions are provided.

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

Carbide-Modified Pd on ZrO2 as Active Phase for CO2-Reforming of Methane—A Model Phase Boundary Approach

Starting from subsurface Zr0-doped “inverse” Pd and bulk-intermetallic Pd0Zr0 model catalyst precursors, we investigated the dry reforming reaction of methane (DRM) using synchrotron-based near ambient pressure in-situ X-ray photoelectron spectroscopy (NAP-XPS), in-situ X-ray diffraction and catalytic testing in an ultrahigh-vacuum-compatible recirculating batch reactor cell. Both intermetallic precursors develop a Pd0–ZrO2 phase boundary under realistic DRM conditions, whereby the oxidative segregation of ZrO2 from bulk intermetallic PdxZry leads to a highly active composite layer of carbide-modified Pd0 metal nanoparticles in contact with tetragonal ZrO2. This active state exhibits reaction rates exceeding those of a conventional supported Pd–ZrO2 reference catalyst and its high activity is unambiguously linked to the fast conversion of the highly reactive carbidic/dissolved C-species inside Pd0 toward CO at the Pd/ZrO2 phase boundary, which serves the role of providing efficient CO2 activation sites. In contrast, the near-surface intermetallic precursor decomposes toward ZrO2 islands at the surface of a quasi-infinite Pd0 metal bulk. Strongly delayed Pd carbide accumulation and thus carbon resegregation under reaction conditions leads to a much less active interfacial ZrO2–Pd0 state.

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

Kinetics of 2-Methylfuran Acylation with Fatty Acid Anhydrides for Biorenewable Surfactants

Friedel–Crafts acylation using long-chain fatty acid derivatives and biomass-derived furans is the key reaction to produce “alkyl furan ketones”, an important precursor for biorenewable oleo-furan surfactants. In this work, the steady-state acylation kinetics and reaction mechanism were investigated using a model system of 2-methylfuran and n-octanoic anhydride in a fixed-bed tubular reactor using Al-MCM-41, a mesoporous aluminosilicate. An apparent activation energy (15.5 ± 1.4 kcal mol –1 ) was obtained for the formation of the acylated product, 2-octanoyl-5-methylfuran (2O5MF), for a temperature range of 348–408 K. Furthermore, the apparent reaction rate orders were ~0.6 and ~0.5 in the 2-methylfuran and anhydride concentrations, respectively, while near-zero apparent rate orders were measured in the product concentrations, indicating negligible product inhibition. An Eley–Rideal catalytic acylation mechanism was proposed to explain the experimentally observed apparent rate orders.

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

FeS colloids – formation and mobilization pathways in natural waters

For this study, we have used synchrotron-based X-ray absorption spectroscopy (structure of Fe–S clusters), transmission electron microscopy (solid-phase crystallinity), Fourier-transform ion-cyclotron-resonance mass spectrometry (identity and composition of natural organic carbon compounds), inductively coupled plasma optical emission spectrometry (total aqueous Fe), and the revised ferrozine method (aqueous Fe( II ) and Fe( III ) concentrations) to determine the stability and nature of colloids generated by sulfidation of ferrihydrite nanoparticles in the absence and presence of organic compounds. We observed that reductive dissolution of ferrihydrite by aqueous sulfide generates nm-scale FeS clusters. Their subsequent aggregation, which promotes settling of FeS aggregates into the solid fraction, was directly correlated with sulfide/Fe ratio. At sulfide/Fe ratios ≤0.5, FeS clusters and larger colloids remained in suspension for at least 14 days (and up to several months). At sulfide/Fe ratios >0.5, sulfidation reaction rates were rapid and FeS cluster aggregation was accelerated. Moreover, the presence of organic compounds increased the time of suspension of FeS colloids, whereas increased ionic strength inhibited the generation of FeS colloids. We present a general conceptual model to predict when and where FeS colloids can form and enhance or inhibit the mobility of contaminants and nutrients associated with them. Our study indicates that in lowsalinity fresh groundwater systems poor in sulfate (i.e. low sulfidation; lakes, floodplains, peatlands etc.), the ferrihydrite sulfidation reaction generates aqueous FeS clusters and larger colloids that remain suspended over long time periods, thus mobilizing a substantial fraction of the total aqueous Fe and S.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A Mesopore-Confined and Graphene Oxide-Localized Ruthenium Catalyst Increases Rates of Mid-Chain Polyolefin Hydrogenolysis

A catalyst architecture with mesoporous silica coating Ru nanoparticles on reduced graphene oxide (mSiO 2 /Ru/rGO) localizes the 2 nm Ru solely at the closed bottoms of 2.9 nm diameter mesopores. mSiO 2 /Ru/rGO catalyzes the rapid, selective hydrogenolysis of polyolefins at wax formation rates (ν wax ) up to 1700 g wax ·g Ru –1 ·h –1 and a turnover frequency (TOF) for C–C bond cleavage of 130 ± 8 min –1 . The ν wax and TOF for mSiO 2 /Ru/rGO are 23× and 16× those of nonporous Ru/rGO, respectively, indicating that faster chain cleavage is also more selective inside mesopores. The methane yield from mSiO 2 /Ru/rGO is ca. 30% of the value obtained from Ru/rGO. Methane decreases, and wax selectivity improves with narrow-pore (2.3 nm) mSiO 2 /Ru/rGO. The reaction rate decreases at lower and higher H 2 pressures from its maximum at 37 bar. Log(TOF)-log(P H2 ) plots reveal rates ∝ [P H2 ] 2.8 or [P H2 ] −2.0 in the lower or higher pressure ranges, respectively. Corresponding plots for Ru/C give rates ∝ [P H2 ] 1 or [P H2 ] −2.4 . Ru is hydrocarbon-covered at low pressure; thus, the higher H 2 order for mSiO 2 /Ru/rGO indicates that mesopores increase the density of cleavable unsaturated moieties on the Ru surface. The lower H 2 order on Ru/C implies a lower density of cleavable groups because saturated segments occupy a larger fraction of the Ru surface. Higher surface occupancy by the hydrocarbon reactant correlates with more methane formation. Therefore, pore confinement, narrower pore diameter in mSiO 2 /Ru/rGO, or higher H 2 pressure leads to lower methane yields. Ru nanoparticles in mSiO 2 /Ru/rGO maintain equivalent activity and selectivity over five recycling tests.

Catalysts↗
DOE OSTI2022

Inferring small-sample-reactivity worth from oscillation experiment with auto-corrected local flux

To estimate reactivity changes as a function of time from the signal of detectors count rates, one usually uses point kinetic equations. However, when detectors are located inside the core, it is often necessary to compute a correction factor to account for the flux-shape deformation (i.e. the change of detector efficiency that may have been induced by neutron absorbers movements). Various simulation methods have been developed depending on transients observed. They range from the resolution of static equations up to full-time-dependent neutron transport in 3D. However, those methods require high computational-cost calculations of local reaction rates that rely on nuclear data evaluations. Therefore, a 'measured' reactivity is never determined through pure measurements. Taking advantage of the known oscillation movement of a small sample at the center of the MINERVE reactor, the innovative algorithm presented in this paper derives the correction factor needed to get no time drift on reactivity steps. Derived correction factors are in good agreement with those obtained independently with TRIPOLI-4 through sensitivity analysis. Low sensitivity to the sample composition and moderate sensitivity to the detectors position reinforce the reliability of those results. This application paves the way to almost pure nuclear data reactivity measurements. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Revealing the interplay between “intelligent behavior” and surface reconstruction of non-precious metal doped SrTiO 3 catalysts during methane combustion

The impact of surface reconstruction of a model perovskite, SrTiO 3 (STO), on CH 4 activation for combustion and oxidative coupling was previously revealed that the reaction rate was proportional to the creation of Srterminated step sites. Doped perovskites (SrTi 1-x M x O 3 , M=metal dopant) present yet another form of reconstruction throughout the surface and the bulk, where the metal dopant can migrate in and out of the perovskite lattice, also known as "intelligent behavior". In this work, understanding the interplay between perovskite surface reconstruction (surface termination) and the "intelligent behavior" is tackled for the first time, and the catalytic consequences are probed with CH 4 combustion as a model reaction. A set of experimental techniques, including XRD, Raman spectroscopy, X-ray adsorption spectroscopy, kinetic measurements, as well as DFT calculations were used to understand the catalytic behavior of the reconstructed surfaces of Ni and Cu-doped STO for methane combustion. Here, we found that during methane oxidation, the diffusion of Ni and Cu into the lattice due to the "intelligent behavior" is accompanied by Sr enrichment on the surface of the perovskite. This Srenrichment process is reversible when Cu or Ni species exsolute as clusters/nanoparticles upon H 2 treatment. Such a surface reconstruction is found to greatly impact the catalytic activity of doped perovskites towards methane combustion.

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

Role of Silica Support in Phosphoric Acid Catalyzed Production of p -Xylene from 2,5-Dimethylfuran and Ethylene

p-Xylene is a commodity chemical of industrial importance for terephthalic acid production, for which renewable sourcing from naturally abundant lignocellulosic biomass is highly desired. Previous work demonstrated that phosphoric acid stabilized on siliceous zeolite supports (e.g., P-BEA, P-SPP) exhibits high selectivity toward p-xylene (>97%) from 2,5-dimethylfuran (DMF) and ethylene. However, the effect of the support and the contribution of heterogeneous versus homogeneous phosphoric acid on the observed catalytic behavior in the solvated reaction system have not been addressed. Here, we determine the phosphoric acid catalytic activity for DMF conversion and its selectivity to p-xylene when it is supported on a silica support as well as in the absence of a support. Specifically, phosphoric acid catalysis was studied in three different scenarios: (1) phosphoric acid was added in the liquid reaction mixture in the absence of any solid support, (2) phosphoric acid was added in the liquid reaction mixture along with inert silica support including siliceous zeolite (i.e., allowing for phosphoric acid–support assembly to proceed in the reaction mixture), and (3) phosphoric acid was first impregnated on the siliceous zeolite support and then the preassembled supported phosphoric acid catalyst was added in the liquid reaction mixture. We found that the reaction rate and selectivity to p-xylene are different in the above scenarios reflecting the effect of the solid support on the catalytic performance of phosphoric acid. In scenario 1, a low concentration of phosphoric acid (1.7 mM) in the absence of any solid support exhibited high selectivity to p-xylene (80% selectivity to p-xylene at 60% conversion of DMF), which decreased with increasing acid concentration. The selectivity to p-xylene and activity of phosphoric acid significantly increased by adding a silica support into the reaction system (scenario 2). Furthermore, this improvement was attributed to phosphoric acid partial association with the surface of the silica support under the reaction conditions (in situ catalyst assembly). Phosphoric acid predeposited on siliceous zeolite supports (e.g., P-BEA, P-SPP) synthesized via impregnation prior to the reaction (scenario 3) catalyzed the reaction heterogeneously without noticeable leaching and exhibited the highest activity and selectivity to p-xylene, suggesting an important role of the silica support and the need to ensure that phosphoric acid acts as a heterogeneous catalyst in order to accomplish selective conversion of DMF to p-xylene.

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

Water‐mediated Decomposition Pathways of ε‐Hexanitrohexaazaisowurtzitane

The dominant initial decomposition mechanism of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (HNIW) is understood from the literature to be the N─NO 2 bond scission unimolecular mechanism. Here we have found a water-mediated HONO release initial decomposition mechanism energetically competitive with the N─NO 2 bond scission mechanism modeled using ab initio nudged elastic band simulations. The activation energy of the water-mediated HONO mechanism was calculated to be 138.9 kJ/mol with the Perdew, Burke, and Ernzerhof (PBE) functional (177.8 with PBE0), while the double NO 2 release mechanism is 171.9 kJ/mol (203.3 with PBE0). Branching secondary, tertiary, and quaternary decomposition steps were also discovered including oxidation of HONO or H 2 O causing ring-opening that leads to C─N bond breaking and N 2 O or NO release. The reaction rate of HONO oxidation is much faster than HONO release, making the HONO concentration low. The Helmholtz free energy barrier of the double NO 2 release mechanism is lower than the water-mediated HONO release above 438 K (165°C) due to the vibrational contribution to the free energy. This helps explain why the literature reports more NO 2 release at higher temperatures. The present study reveals a physical mechanism for how water can catalyze the decomposition of HNIW at low enough temperatures, and reveals secondary, tertiary, and quaternary reaction mechanisms leading to NO x release. In conclusion, as water is realistically always present, the reported barriers are important for building kinetic models needed to ensure the safety and reliability of HNIW.

Steele, Brad A. [Lawrence Livermore National Labor↗
DOE OSTI2024

Probing time-resolved plasma-driven solution electrochemistry in a falling liquid film plasma reactor: Identification of HO$^−_2$ as a plasma-derived reducing agent

Many applications involving plasma–liquid interactions depend on the reactive processes occurring at the plasma–liquid interface. In this article, we report on a falling liquid film plasma reactor allowing for in situ optical absorption measurements of the time-dependence of the ferricyanide/ferrocyanide redox reactivity, complemented with ex situ measurement of the decomposition of formate. We found excellent agreement between the measured decomposition percentages and the diffusion-limited decomposition of formate by interfacial plasma-enabled reactions, except at high pH in thin liquid films, indicating the involvement of previously unexplored plasma-induced liquid phase chemistry enabled by long-lived reactive species. We also determined that high pH facilitates a reduction-favoring environment in ferricyanide/ferrocyanide redox solutions. In situ conversion measurements of a 1:1 ferricyanide/ferrocyanide redox mixture exceed the measured ex situ conversion and show that conversion of a 1:1 ferricyanide/ferrocyanide mixture is strongly dependent on film thickness. We identified three dominant processes: reduction faster than ms time scales for film thicknesses >100 µm, •OH-driven oxidation on time scales of <10 ms, and reduction on 15 ms time scales for film thickness <100 µm. We attribute the slow reduction and larger formate decomposition at high pH to HO$^−_2$ formed from plasma-produced H 2 O 2 enabled by the high pH at the plasma–liquid interface as confirmed experimentally and by computed reaction rates of HO$^−_2$ with ferricyanide. Overall, this work demonstrates the utility of liquid film reactors in enabling the discovery of new plasma-interfacial chemistry and the utility of atmospheric plasmas for electrodeless electrochemistry.

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

Understanding the Origin of Negative Temperature Dependence and Activity of N-Coordinated Cobalt Sites During Ethylene Dimerization

The on-demand production of short-chain linear alpha olefins (LAOs; C4-C8) via C2H4 dimerization and oligomerization is industrially attractive, prompting extensive research on designing active, selective, and stable catalysts for industrial use. Cobalt supported on ammoniated carbon (Co(NH3)x/C) catalysts have shown remarkable activity and selectivity in this process. However, critical aspects such as the active phase, active site structure, the role of the catalyst support, cobalt loading effects, and the inverse correlation of the reaction rate with temperature remain inadequately understood. This study systematically explores these factors using a combination of steady-state differential catalytic tests, in situ molecular characterization including diffuse reflectance UV-Vis (DR-UV-Vis), Infrared, and Raman spectroscopies, and ex situ X-ray diffraction (XRD) and high annular aberration-corrected dark field transmission electron microscopy (HAADF-STEM). Various supports (SiO2, Al2O3, NH4-ZSM-5, g-C3N4, and C) and cobalt loadings (1.0-3.0 Co nm-2) were studied to determine the optimal catalyst composition and identify the active phase and sites. Carbon-supported catalysts uniquely produce C4-8 LAOs during C2H4 dimerization, with site-time-yield remaining constant (~10-3 s-1) for 1.0-4.0 Co nm-2 at prolonged reaction times (24-48?h time-on-stream). At higher loadings of 6.0 Co nm-2, the formation of crystalline CoO and Co3O4 phases reduces catalytic activity and LAO selectivity. Our findings show that active catalysts lack crystalline cobalt oxides and instead feature dispersed Co2+ sites, tetra-coordinated to a mix of N/NH3 and O/H2O ligands, which catalyze C2H4 dimerization via the Cossee-Arlman mechanism, exhibiting 1st order dependence on C2H4 concentration. The observed inverse rate-temperature correlation is attributed to compensation effects (i.e., presence of Cremer-Constable relationship) linked to changes in adsorption enthalpic and entropic factors.

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

Rapid Construction of Validated Chemistry Models for Advanced Biofuels (Co-Optima Final Report)

Computer methods were developed and demonstrated that make it possible to predict the combustion and pyrolysis chemistry of proposed biofuels, even before the biofuels are ever manufactured or experimentally tested. The methods are based on first-principles calculations of reaction rates and equilibria. A large number of complicated calculations are required, so we developed a computer workflow that automates many of the steps. With the chemical structure of the biofuel and the system conditions (temperature, pressure, compositions, etc.) as the input, the computer methods can construct a detailed kinetic model with refined parameters without intense human intervention. The automation of model construction makes it practical for a single engineer to make the predictions in a reasonable time period. The accuracy of the computer predictions was tested by measuring the time-profiles of various reaction intermediates (carbon monoxide, carbon dioxide, water, formaldehyde, methane, and ethylene) formed as the biofuels pyrolyze or ignite using sophisticated laser shock tube methods. Several biofuels, including methyl propyl ether, cyclopentanone, butyl acetate isomers, etc., were carefully studied to demonstrate the validity of the proposed computational and experimental methodology. This research project opens up the possibility for assessing some proposed biofuels on the computer without expensive experimental campaigns. The new computer methods are open source and well-documented, allowing a knowledgeable modeler to easily access the developed methods and implement them in their work routine. Moreover, they are also convenient and powerful, promising for developing the fuel chemistry models needed for use in engine simulators, facilitating co-optimization of new engine designs with new fuels.

09 BIOMASS FUELS↗
DOE OSTI2022

Multitask Machine Learning of Collective Variables for Enhanced Sampling of Rare Events

Computing accurate reaction rates is a central challenge in computational chemistry and biology because of the high cost of free energy estimation with unbiased molecular dynamics. In this work, a data-driven machine learning algorithm is devised to learn collective variables with a multitask neural network, where a common upstream part reduces the high dimensionality of atomic configurations to a low dimensional latent space and separate downstream parts map the latent space to predictions of basin class labels and potential energies. Here, the resulting latent space is shown to be an effective low-dimensional representation, capturing the reaction progress and guiding effective umbrella sampling to obtain accurate free energy landscapes. This approach is successfully applied to model systems including a 5D Müller Brown model, a 5D three-well model, the alanine dipeptide in vacuum, and an Au(110) surface reconstruction unit reaction. It enables automated dimensionality reduction for energy controlled reactions in complex systems, offers a unified and data-efficient framework that can be trained with limited data, and outperforms single-task learning approaches, including autoencoders.

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

River Dynamics Control Transit Time Distributions and Biogeochemical Reactions in a Dam-Regulated River Corridor

Residence Time Distributions (RTDs) exert an important control on biogeochemical translation in the hydrological systems. Previous tracer studies have revealed that RTDs often followed time-invariant exponential, lognormal, or heavy-tailed RTDs that have power-law behavior for long tails in headwater or low-order streams. However, there is a recognition that RTDs can be more complicated and time-variable in response to dynamic hydrological forcing. Here, we use particle tracking to estimate RTDs along the Hanford Reach of the Columbia River and to quantify the influences of river stage fluctuations. Particle tracking is conducted using the velocity field from high-resolution 3D groundwater flow simulations. The effects of hydrological forcing on the residence time distribution were evaluated by varying river flow boundary conditions and releasing particles in different time windows. Our results revealed that dynamic stage fluctuations created rapidly changing losing-gaining conditions in the river and led to highly transient RTDs, which contributes to multiple modes of RTDs. Dam-induced high-frequency (sub-daily) flow variation contributes to the short-time (sub-daily) component of the RTDs. Deviation of the reactant consumption under the single-mode assumption compared to the multimodal RTDs is relatively small (~5%) and appears when the Damköhler number was close to one.More specifically, high-frequency (daily) dam-induced stage variations have higher impacts on biogeochemical reactions with faster reaction rate, and can potentially increase the reactant consumption by 26.68%. Based on these findings, we suggest that current river basin models could be improved by including bank storage and more complex RTDs influenced by both short and long term river stage fluctuations.

54 ENVIRONMENTAL SCIENCES↗