Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical kinetics and dynamics”

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 91 records · Page 5

Kinetic Deep Learning v0.1

Here, we present a method that uses protein levels to predict times series of metabolite concentrations. Understanding this type of pathway dynamics is important in order to predict the behavior of the pathway and, more pragmatically, to be able to design biological systems (such as strains bioengineered to produce chemical products) reliably. Typically, for this purpose, kinetic models consisting of differential equations based on the Michaelis-Menten dynamics have been used in the past. However, these methods can rarely produce good fits to measured data time series. Possibly, this happens because the kinetic constants are unknown or are different from the ones measured in vivo, or perhaps because Michaelis-Menten dynamics is not a satisfactory description. In order to improve the predictive nature of these kinetic models we have eliminated the Michaelis-Menten description of pathway dynamics and we have substituted it by algorithms that automatically learn these dynamics from previously obtained metabolomics and proteomics data using machine learning approaches. Specifically, kinetic deep learning uses deep learning to map proteomics time series to metabolite concentration time series, instead of learning the first metabolite derivative and integrating in (as in the first version of kinetic learning). This approach is shown to provide good to excellent results with a data set specifically collected for this purpose.

Garcia Martin, Hector [Joint BioEnergy Institute (↗

Surface chemistry models for GaAs epitaxial growth and hydride cracking using reacting flow simulations

Hydride vapor phase epitaxy (HVPE) is a promising technology that can aid in the cost reduction of III-V materials and devices manufacturing, particularly high-efficiency solar cells for space and terrestrial applications. However, recent demonstrations of ultrafast growth rates (~500 µm/h) via uncracked hydrides are not well described by present models for the growth. Therefore, it is necessary to understand the kinetics of the growth process and its coupling with transport phenomena, so as to enable fast and uniform epitaxial growth. In this work, we derive a kinetic model using experimental data and integrate it into a computational fluid dynamics simulation of an HVPE growth reactor. We also modify an existing hydride cracking model that we validate against numerical simulations and experimental data. Here, we show that the developed growth model and the improved cracking model are able to reproduce experimental growth measurements of GaAs in an existing HVPE system.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗

Catalytic site ensembles: A context to reexamine the Langmuir-Hinshelwood kinetic description

The Langmuir-Hinshelwood formalism describes catalytic reactions of Langmuirian surface species under the assumption that all adsorbates are randomly-distributed – enabling adjacency of surface-bound intermediates to be determined solely by coverages of single-site occupants. We demonstrate herein that this approximation is inappropriate even for simple catalytic reactions (e.g. A + A → A 2 ) and manifestly neglects islanding of slowly-consumed species and partitioning of highly-reactive species inherently engendered by ≥ two-site elementary steps (e.g. A*–A* → A 2(g) + *–*). Rigorous description of kinetically-consequential islanding/partitioning phenomena requires explicit description of the coverage and chemical dynamics of all multi-site ensembles. Higher-order, ensemble-specific rate terms identify the particular microscopic events relevant to each ensemble, and, in doing so, reveal that each elementary step (e.g. A (g) adsorption) describes not one event (e.g. A (g) + * → A*) , but a sum over all ensemble-specific paths (e.g. A (g) + *–* → A*–* and A (g) + A*–* → A*–A*). De-convoluting each elementary step into its constituent multi-site paths proffers kinetic detail otherwise inaccessible – enabling (i) identification of rate- and selectivity-determining site ensembles, (ii) calculation of rates and degrees of rate control of ensemble-specific elementary steps, (iii) incorporation of adsorbate surface diffusion, (iv) incorporation of lateral adsorbate interactions, and (v) quantitative description of catalysis of multi-site-occupying intermediates (e.g. *C n H m * species in hydrocarbon (de-)hydrogenation and C-C bond coupling/cleavage reactions) which we demonstrate here is inaccessible to the Langmuir-Hinshelwood formalism even if adsorbate surface diffusion is infinitely-fast.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomic dynamics of gas-dependent oxide reducibility

Understanding oxide reduction is critical for advancing metal production, catalysis and energy technologies. Although carbon monoxide (CO) and hydrogen (H 2 ) are widely used reductants, the mechanisms by which they work are often presumed to be similar, both involving lattice oxygen removal. However, because of growing interest in replacing CO with H 2 to lower CO 2 emissions, distinguishing gas-specific reduction pathways is critical. Yet, capturing these atomic-scale processes under reactive gas and high-temperature conditions remains challenging. Here we use environmental transmission electron microscopy, which is capable of real-time, atomic-resolution imaging of gas–solid redox reactions to directly visualize the gas-dependent oxide reduction dynamics in NiO. We show that CO drives surface nucleation and the growth of metallic Ni islands, leading to self-limiting surface metallization. Conversely, H 2 activates a coupled surface-to-bulk transformation, where protons from dissociated H 2 infiltrate the oxide lattice to promote the inward migration of surface-generated oxygen vacancies and enabling bulk metallization. By contrast, oxygen vacancies formed by CO remain confined near the surface, where they rapidly form a metallic Ni layer that inhibits further reduction. Furthermore, these results reveal distinct atomistic pathways for CO and H 2 and provide insights that may guide metallurgical processes and catalyst design.

36 MATERIALS SCIENCE↗

Solvated Electrons Have Multiple Personalities in Molten Salts

P-Crosscut: Solvated Electrons Have Multiple Personalities in Molten Salts [EFRC – MSEE] Alejandro Ramos-Ballesteros;1 Hung H. Nguyen;2 Kazuhiro Iwamatsu;3 Santanu Roy;4 Vyacheslav Bryantsev;4 Michael E. Woods;1 Ruchi Gakhar;1 Phillip Halstenberg;4 Bobby Layne;5 Jay A. LaVerne;6 Claudio J. Margulis:2* and James F. Wishart5* 1Idaho National Laboratory; 2The University of Iowa; 3Hunter College; 4Oak Ridge National Laboratory; 5Brookhaven National Laboratory; 6University of Notre Dame Abstract: The solvated (eS–) is a powerful reducing agent and one of the primary products of molten salt radiolysis. In these extreme high-temperature environments, the eS– can initiate cascades of redox processes that significantly alter the physical and chemical properties of a molten salt, posing challenges for molten salt reactor (MSR) performance and longevity. Consequently, mastering the fundamental behavior of the eS– could enable the design of specific molten salt mixtures with tailored Lewis acidities for controlling the speciation and chemical reactivity of the eS–, thereby mitigating its overall impact on MSR technologies. Here, we use time-resolved electron pulse radiolysis techniques for determining chemical kinetics and transient absorption spectra, combined with ab initio molecular dynamics simulations to explore the influence of multivalent metal cations on the fundamental speciation and distribution of eS– coordination environments.

38 - RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCL↗

An out-of-distribution-aware autoencoder model for reduced chemical kinetics

While detailed chemical kinetic models have been successful in representing rates of chemical reactions in continuum scale computational fluid dynamics (CFD) simulations, applying the models in simulations for engineering device conditions is computationally prohibitive. To reduce the cost, data-driven methods, e.g., autoencoders, have been used to construct reduced chemical kinetic models for CFD simulations. Despite their success, data-driven methods rely heavily on training data sets and can be unreliable when used in out-of-distribution (OOD) regions (i.e., when extrapolating outside of the training set). In this paper, we present an enhanced autoencoder model for combustion chemical kinetics with uncertainty quantification to enable the detection of model usage in OOD regions, and thereby creating an OOD-aware autoencoder model that contributes to more robust CFD simulations of reacting flows.Here, we first demonstrate the effectiveness of the method in OOD detection in two well-known datasets, MNIST and Fashion-MNIST, in comparison with the deep ensemble method, and then present the OOD-aware autoencoder for reduced chemistry model in syngas combustion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Intercalation of Lithium into Graphite: Insights from First-Principles Simulations

Understanding ion intercalation at electrode–electrolyte interfaces is key to the development of energy storage and water desalination. Here, we investigate Li + kinetics at a prototypical interface between graphite anodes and an organic electrolyte, and we elucidate key factors that determine ion transport, using first-principles methodology coupling ab initio molecular dynamics simulations with a solvation model. We show that surface chemical composition significantly influences the kinetics of ion intercalation from the liquid into graphite. We find that this is partly related to the ion desolvation process, which varies notably for different graphite surface chemical terminations. In addition, interfacial polarization is found to play an important role in determining energy barriers for ion transfer. We also discuss the impact of electrode potentials, which is often neglected in conventional first-principles calculations despite being a key factor in device configurations. Our study provides insights into the coupling of electronic and ionic effects of interfacial chemistry on ion transport at complex electrode–electrolyte interfaces.

25 ENERGY STORAGE↗

Coherence mapping to identify the intermediates of multi-channel dissociative ionization

Identifying the short-lived intermediates and reaction mechanisms of multi-channel radical cation fragmentation processes remains a current and important challenge to understanding and predicting mass spectra. We find that coherent oscillations in the femtosecond time-dependent yields of several product ions following ultrafast strong-field ionization represent spectroscopic signatures that elucidate their mechanism of formation and identify the intermediate(s) they originate from. Experiments on endo-dicyclopentadiene show that vibrational frequencies from various intermediates are mapped onto their resulting products. Aided by ab initio methods, we identify the vibrational modes of both the cleaved and intact molecular ion intermediates. These results confirm stepwise and concerted fragmentation pathways of the dicyclopentadiene ion. This study highlights the power of tracking the femtosecond dynamics of all product ions simultaneously and sheds further light onto one of the fundamental reaction mechanisms in mass spectrometry, the retro-Diels Alder reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time-Gated Raman Spectroscopy (Intern Final Report)

Time-gated Raman spectroscopy is a powerful tool for rapid non-destructive analysis of chemical compounds, particularly when combined with other investigative methods. Time-gating the signal obtained from a sample during Raman spectroscopy allows for the removal of background emission and black body radiation, improving clarity and accuracy of scans. The main objectives of this project were to prove the lab’s time-gated Raman spectroscopy system works as proof of concept and to make select improvements to the in-house LabVIEW software used to run the time-gated system. This time-gated system will be added to the CTK’s TAP reactors in the future, but this is out of the scope of this internship project. This research is important, as understanding the mechanisms behind catalysts leads to designing better systems/materials and the perfection of catalysts, and energy-efficient chemical manufacturing is key when designing a greener future. The research for this project was done with samples at both ambient conditions in sample holders and at operando conditions (i.e., high temperature, unsteady-state) inside a small-volume chemical reactor. Numerous spectro-kinetic scans of the catalysts were obtained, but the advanced dynamics of these systems are still being interpreted. There is data that demonstrates that the time-gated system is correctly rejecting emission from samples, and one can observe real-time changes in the signal from the sample, due to coking or changes induced by redox chemistry. Chemical manufacturers and the planet are the main beneficiaries of catalysis research, as better catalysts will reduce carbon emissions.

36 MATERIALS SCIENCE↗

Complex and Unusual Excited-State Relaxation Dynamics of 9,9′-Bifluorenylidene Revealed by Comprehensive Time-Resolved Spectroscopy and MRSF-TDDFT Calculations

Design of novel photochemical molecular motors often requires molecular building blocks that exhibit rather unusual photoactivity, for which conventional analyses of spectroscopic data can lead to conflicting interpretations. We here systematically investigated the excited-state relaxation dynamics of one such molecule, 9,9′-bifluorenylidene (BF), through comprehensive and complementary integration of ultrafast transient absorption (TA) and femtosecond stimulated Raman (FSRS) spectroscopies and first-principles mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT). TA and FSRS identified two sequentially formed transients following photoexcitation. The decay kinetics of the two intermediates differ in response to excitation wavelengths and the viscosity/polarity of solvents. MRSF-TDDFT calculations reveal a direct, barrierless internal-conversion pathway from the bright Franck−Condon state to a dark S1 minimum, where the excited-state population is transiently trapped, accounting for the first transient species observed in spectroscopic experiments. Further tracking down along the PES with MRSF-TDDFT mapped out two nonradiative relaxation pathways via conical intersections that connect the dark S 1 state to three configurations in the ground-state manifolds, within which a ring structure with a C8−C8′ bond and the vibrationally excited ground-state BF were identified from spectroscopic and kinetic data. The complexity of relaxation kinetics was attributed to the flexible torsional and twisting motions about the C9−C9′ bridge bond enabled by the diradical character of the S 1 state. These findings clarify unusual photoactive relaxation dynamics stemming from a novel correlation between structural flexibility and shifting electronic characteristics, and they demonstrate the importance of integrating spectroscopic and advanced electronic structure calculation studies for judicious clarification of complex, competing relaxation pathways of excited states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Universal progression of structure and dynamics in colloidal nanocrystal gels during salt-accelerated aging

Controlling the structure and function of colloidal gels requires a detailed understanding of how the various components govern network formation and aging. In particular, molecular additives like salts are widely used to tune interparticle interactions, yet their influence on gelation pathways in complex systems such as colloidal nanocrystal gels remains inadequately understood. Here, we investigate how noncoordinating salts modulate the evolution of gels formed using chemically linked tin-doped indium oxide nanocrystals. Through combined structural, dynamic, and kinetic analyses, we demonstrate that increasing salt concentration accelerates gelation. When rescaled by salt-dependent characteristic times, the evolution collapses onto universal trajectories, revealing a time-salt superposition principle. The universality extends across length scales, suggesting a consistent salt-dependent mechanism that controls both local structuring and macroscopic network formation. This observed salt modulation of structure and dynamics provides a predictive basis for controlling the kinetics of nonequilibrium nanocrystal gel assembly, enhancing the rational design of functional nanomaterials with tunable properties.

36 MATERIALS SCIENCE↗

Influence of Nd(II) and Nd(III) Ions on the Speciation and Kinetic Dynamics of Radiolytic Transients in Molten LiCl-KCl Eutectic Salt Mixtures

In the pyrochemical reprocessing of used nuclear fuel, understanding the effects of ionizing radiation on the chemical behavior of fission products in molten chloride eutectic salt is crucial. Here, this study investigates the impact of radiation on the chemical and kinetic speciation of Nd ions in molten LiCl-KCl eutectic, which are challenging to separate due to their multivalent oxidation states and comproportionation reactions. Electron pulse radiolysis techniques were used to determine transient absorption spectra, chemical kinetics, and activation parameters for the reaction of Nd ions with the primary products of molten LiCl-KCl eutectic radiolysis: the solvated electron (e S - ) and dichlorine radical anion (Cl 2 •- ). Initially, Nd(III) reacted rapidly with e S - , forming Nd(II) with a second-order rate coefficient of k = (4.54 ± 0.07) × 10 10 M -1 s -1 at 673 K. The resulting Nd(II) ions then reacted slowly with Cl 2 •- , regenerating Nd(III) with an estimated rate of k = (1.72 ± 0.04) × 10 10 M -1 s -1 at 673 K. Additionally, the comproportionation reaction of Nd(III) and Nd(s) to form Nd(II) was monitored chemically and electrochemically. This work suggests that the comproportionation reaction of Nd(III) with Nd(s) to form Nd(II) is unstable in molten LiCl-KCl eutectic.

36 - MATERIALS SCIENCE↗

Evaluating Chemical Kinetics Predictions for Propane Using 3-D and 0-D Models in a Boosted Spark-Ignited Engine

Propane has been shown to be a promising alternative fuel to reduce emissions while simultaneously achieving high efficiencies in medium- and heavy-duty engines. These high-power density applications require boosted engines which, combined with high compression ratio, can lead to auto-ignition and knock. While three-dimensional (3-D) computational fluid dynamics (CFD) models are often used for resolving the complex fluid flow in engines, these models can become computationally expensive when simulating detailed chemical kinetics. Likewise, zero-dimensional (0-D) models are computationally concise enough for kinetics development, but lack any flow-field information which governs the flame propagation processes in spark ignition (SI) engines. This work presents a comprehensive comparison between 3-D and 0-D closed cycle simulations at knocking conditions in a high compression ratio high stroke-to-bore ratio propane engine. In order to initialize the flow-field for the 3-D closed cycle (intake valve closing, (IVC) to exhaust valve opening, (EVO)) simulation, a motored multi-cycle 3-D model was run using Converge to create a map at IVC, reducing the computational time. The map allowed a non-homogeneous 3-D closed cycle simulation to be satisfactorily validated against experiments, while a homogeneous case using only the turbulence field mapping was also simulated, mimicking 0-D modeling. The 3-D simulations were used to prescribe the initial conditions (e.g., IVC thermodynamics, speciation, burn-rate profile) for a 2-zone 0-D SI engine model in Chemkin Pro for both cases. It was found that 2-zone 0-D modeling underpredicted the knock onset timing, likely due to the lack of thermal stratification in the unburned gas region. Future work will carry multi-zone 0-D modeling to capture the fuel auto-ignition in the unburned region.

Douvry-Rabjeau, Julien [Oakland University, Roches↗

Nonequilibrium Solvent Effects during Photodissociation in Liquids: Dynamical Energy Surfaces, Caging and Chemical Identity

In the gas phase, potential energy surfaces can be used to provide insight into the details of photochemical reaction dynamics. In solution, however, it is unclear what potential energy surfaces, if any, can be used to describe even simple chemical reactions such as the photodissociation of a diatomic solute. In this paper, we use mixed quantum/classical (MQC) molecular dynamics (MD) to study the photodissociation of $Na$ $^+_2$ in both liquid Ar and liquid tetrahydrofuran (THF). We examine both the gas-phase potential surfaces and potentials of mean force (PMF), which assume that the solvent remains at equilibrium with the solute throughout the photodissociation process and show that neither resemble a nonequilibrium dynamical energy surface that is generated by taking the time integral of work. For the photodissociation of $Na$ $^+_2$ in liquid Ar, the dynamical energy surface shows clear signatures of solvent caging, and the degree of caging is directly related to the mass of the solvent atoms. For $Na$ $^+_2$ in liquid THF, local specific interactions between the solute and solvent lead to changes in chemical identity that create a kinetic trap that effectively prevents the molecule from dissociating. Here, the results show that nonequilibrium effects play an important role even in simple solution-phase reactions, requiring the use of dynamical energy surface to understand such chemical events.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Flame acceleration and transition to detonation in a pre-/main-chamber combustion system

Numerical simulations are performed to study the mechanism of deflagration to detonation transition (DDT) in a pre-/main-chamber combustion system with a stoichiometric ethylene–oxygen mixture. A Godunov algorithm, fifth-order in space, and third-order in time, is used to solve the fully compressible Navier–Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The two-dimensional simulation shows that a laminar flame grows in the pre-chamber and then develops into a jet flame as it passes through the orifice. A strong shock forms immediately ahead of the flame, reflecting off the walls and interacting with the flame front. The shock–flame interactions are crucial for the development of flame instabilities, which trigger the subsequent flame development. The DDT arises due to a shock-focusing mechanism, where multiple shocks collide at the flame front. A chemical explosive mode analysis (CEMA) criterion is developed to study the DDT ignition mode. Preliminary one-dimensional computations for a laminar propagating flame, a fast flame deflagration, and a Chapman–Jouguet detonation are conducted to demonstrate the validity of CEMA on the chemical-diffusive model, as well as to determine the proper conditioning value for CEMA diagnostic. The two-dimensional analysis with CEMA indicates that the DDT initiated by the shock-focusing mechanism can form a strong thermal expansion region at the flame front that features large positive eigenvalues for the chemical explosive mode and dominance of the local autoignition mode. Thus, the CEMA criterion proposed in this study provides a robust diagnostic for identifying autoignition-supported DDT, of which the emergence of excessive local autoignition mode is found to be a precursor. The effect of grid size, initial temperature, and orifice size are then evaluated, and results show that although the close-chamber DDT is highly stochastic, the detonation initiation mechanism remains robust.

42 ENGINEERING↗

Comprehensive uncertainty quantification (UQ) for full engineering models by solving probability density function (PDF) equation

This report details a new method for propagating parameter uncertainty (forward uncertainty quantification) in partial differential equations (PDE) based computational mechanics applications. The method provides full-field quantities of interest by solving for the joint probability density function (PDF) equations which are implied by the PDEs with uncertain parameters. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known apriori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the driving (dynamics) term of the PDEs based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the dynamics term of the governing PDEs. The efficacy of the method is demonstrated by applying it to two target problems: solid mechanics explicit dynamics with uncertain material model parameters, and reacting hypersonic fluid mechanics with uncertain chemical kinetic rate parameters. A minimally invasive implementation of the method for representative codes SPARC (reacting hypersonics) and NimbleSM (finite- element solid mechanics) and associated software details are described. For solid mechanics demonstration problems the method shows order of magnitudes improvement in accuracy over traditional stochastic collocation. For the reacting hypersonics problem, the method is implemented as a streamline integration and results show very good accuracy for the approximate sample solutions of re-entry flow past the Apollo capsule geometry at Mach 30.

42 ENGINEERING↗

Thermo-mechanical deterioration and molecular degradation of 3D-printed methacrylate-based polymer in various chemical environments

The advancement of additive manufacturing (AM) has accelerated the development of stereolithographic (SLA) photo-curable resins, particularly methacrylate-based polymers, due to the ability to their high-resolution, robust mechanical properties, and suitability. However, their long-term performance in chemical environments remains poorly understood. This study investigates the extent and mechanisms of degradation on an SLA-printed methacrylate-based polymer subjected to various chemicals, including polar and non-polar solvents, as well as strongly acidic aqueous solutions over a 12-week accelerated aging period. A comprehensive analytical approach incorporating swelling kinetics, surface morphology, tensile and dynamic mechanical analysis (DMA), Fourier-transform infrared (FTIR) spectroscopy, and mass spectrometry (MS) was employed to characterize chemical absorption, structural integrity, and leached products. Results reveal that degradation severity is governed by both the polarity and reactivity of the chemical environment. Notably, exposure to 6 mol L -1 HNO₃ induced the most severe deterioration, with over threefold higher swelling compared to other media, and significant reductions in tensile strength, tensile modulus, and glass transition temperature (T g ). In contrast, specimens aged in non-polar solvents (xylene and dodecane) exhibited negligible chemical interaction and retained mechanical performance. FTIR and MS analyses identified acid-catalyzed hydrolysis of ester groups as prominent degradation pathways in acidic media, while diffusion-controlled plasticization prevailed in polar solvents. Furthermore, this study provides valuable insights into the chemical stability of SLA-printed polymers and develops predictive degradation profiles that are crucial for designing durable polymer systems for advanced industrial use.

36 MATERIALS SCIENCE↗