Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic network”

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 109 records · Page 6

KETCHUP: Parameterizing of large-scale kinetic models using multiple datasets with different reference states

Large-scale kinetic models provide the computational means to dynamically link metabolic reaction fluxes to metabolite concentrations and enzyme levels while also conforming to substrate level regulation. However, the development of broadly applicable frameworks for efficiently and robustly parameterizing models remains a challenge. Challenges arise due to both the heterogeneity, paucity, and difficulty in obtaining flux and/or concentration data but also due to the computational difficulties of the underlying parameter identification problem. Both the computational demands for parameterization, degeneracy of obtained parameter solutions and interpretability of results has so far limited widespread adoption of large-scale kinetic models despite their potential. Herein, we introduce the Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the (non)steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP is benchmarked against previously parameterized large-scale kinetic models demonstrating an at least an order of magnitude faster convergence than the tool K-FIT while at the same time attaining better data fits. Finally, this versatile toolbox accepts different kinetic descriptions, metabolic fluxes, enzyme levels and metabolite concentrations, under either steady-state or instationary conditions to enable robust kinetic model construction and parameterization. KETCHUP supports the SBML format and can be accessed at https://github.com/maranasgroup/KETCHUP.

59 BASIC BIOLOGICAL SCIENCES↗

Mapping the Temperature-dependent and network site-specific onset of spectral diffusion at the surface of a water cluster cage

We explore the kinetic processes that sustain equilibrium in a microscopic, finite system. This is accomplished by monitoring the spontaneous, time-dependent frequency evolution (the frequency autocorrelation) of a single OH oscillator, embedded in a water cluster held in a temperature-controlled ion trap. The measurements are carried out by applying two-color, IR-IR photodissociation mass spectrometry to the D3O+?(HDO)(D2O)19 isotopologue of the “magic number” protonated water cluster, H+?(H2O)21. The OH group can occupy any one of the five spectroscopically distinct sites in the distorted pentagonal dodecahedron cage structure. The OH frequency is observed to evolve over tens of milliseconds in the temperature range (90-120 K). Starting at 100 K, large “jumps” are observed between two OH frequencies separated by ~300 cm-1 indicating migration of the OH group from the bound OH site at 3350 cm-1 to the free position at 3686 cm-1. Increasing the temperature to 110 K leads to partial interconversion among many sites. All sites are observed to interconvert at 120 K such that the distribution of the unique OH group among them adopts the form one would expect for a canonical ensemble. The spectral dynamics displayed by the clusters thus offer an unprecedented view into the molecular-level processes that drive spectral diffusion in an extended network of water molecules.

Yang, Nan↗

Deep learning surrogate model for kinetic Landau-fluid closure with collision

In this work, the kinetic Landau-fluid (LF) closure with collision and periodic boundary condition is used in the development of the deep learning (DL) surrogate model. A classical neural network, namely, feedforward neural network or sometimes termed multilayer perceptron, is constructed and trained to learn the kinetic LF closure in the static limit and arbitrary mean free path in configuration space. The preliminary relation between best hyperparameters and critical parameters for data generation is found. Compared with the numerical approach (non-Fourier method) of the LF closure, the deep learning surrogate model shows an order of magnitude of improvement in terms of accuracy. Perhaps most importantly, the surrogate model closure has been integrated for the first time with fluid simulations. Our DL-enabled fluid simulations, for the first time, give the correct Landau damping rate for a wide range of wave vectors, while the Hammett–Perkins closure cannot produce the correct damping rate. We correctly connect the collisionless Hammett–Perkins closure and collisional Braginskii closure to reproduce the intrinsic nonlocal feature of the heat flux with DL techniques. We address the most concerning error accumulation problem and find that simulations with the deep learning surrogate model are as good as, if not better than, simulations with the analytic closure in terms of long-term numerical stability in the linear Landau damping test.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations

The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parametrized from MD simulations of complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four order of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example in the case of large molecules or solids. In this work, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parametrize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation timescales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parametrized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Segmental Dynamics and Vitrification in Associating Copolymer Melts: Role of Cluster Formation, Microdomains, and Cross-Linking

A microscopic statistical mechanical theory of the structure, self-assembly, and activated segmental relaxation is employed to study associating copolymer melts with high attractive sticker fractions, local clustering, and disordered microphase ordering. The stickers are dynamically pinned in a manner that does not affect equilibrium structure which mimics the much slower physical bond breaking process or postassembly cross-linking of sticky monomers. Local sticker clustering and microdomain spatial correlations significantly modify the activated relaxation of nonstickers and glass transition temperature, T g . A re-entrant glass-melting feature is predicted as sticker attraction strength is initially increased corresponding to a speed up of segmental relaxation, and hence reduction of T g relative to the cross-linked homopolymer network. A mechanistic analysis reveals three competing effects: a purely kinetic slowing down of nonstickers down due to cross-linking, disordering of the nonsticker local cage and weakening of effective forces they experience due to sticker physical clustering, and a longer range impact of microdomain scale correlations that results in nonmonotonic dynamical effects. At high enough attraction strength, a qualitative change emerges corresponding to a sticker fraction dependent elevation of T g , which eventually surpasses that of the cross-linked homopolymer network. Dynamically, the new physics arises from a complex evolution of the amplitude of the collective elastic field that dresses the large amplitude mobile segment hopping within a coupled local-nonlocal description of the alpha relaxation. Furthermore, the results are qualitatively consistent with recent experiments on associating PDMS and PPG telechelics of fixed sticker fraction but with chemically different end groups of variable attraction strengths. Finally, possible tests using simulation and the influence of material or model specific interaction potentials and other real world complications are discussed.

Copolymers↗

An image-driven machine learning approach to kinetic modeling of a discontinuous precipitation reaction

Micrograph quantification is an essential component of several materials science studies. Machine learning methods, in particular convolutional neural networks, have previously demonstrated performance in image recognition tasks across several disciplines (e.g. materials science, medical imaging, facial recognition). Here, we apply these well-established methods to develop an approach to microstructure quantification for kinetic modeling of a discontinuous precipitation reaction in a case study on the uranium-molybdenum system. Prediction of material processing history based on image data (classification), calculation of area fraction of phases present in the micrographs (segmentation), and kinetic modeling from segmentation results were performed. Results indicate that convolutional neural networks represent microstructure image data well, and segmentation using the k-means clustering algorithm yields results that agree well with manually annotated images. Classification accuracies of original and segmented images are both 94% for a 5-class classification problem. Kinetic modeling results agree well with previously reported data using manual thresholding. The image quantification and kinetic modeling approach developed and presented here aims to reduce researcher bias introduced into the characterization process, and allows for leveraging information in limited image data sets.

36 MATERIALS SCIENCE↗

Ligation of Single-Site Ruthenium within Perovskite Oxides for Efficient Conversion of Thermodynamically Stable Molecules

Ru cation ligation in SrTiO 3 perovskite and their migration to the surface through exsolution are investigated for the dry reforming of methane (DRM), a chemistry that requires activation of two thermodynamically stable molecules, CH 4 and CO 2 . Compared to a supported 1 wt % Ru/SrTiO 3 benchmark, doped and exsolved Ru-SrTiO 3 demonstrate ≥3× higher CH 4 turnover rates (873 K), with isolated, ligated Ru exhibiting higher reactivity. Reactor studies assert that CH 4 and CO 2 activation are both kinetically relevant for CH 4 turnover rates on exsolved Ru-SrTiO 3 , unlike for supported Ru systems, where H-abstraction from CH 4 is the sole kinetically relevant step. As such, exsolved and doped Ru architectures are responsive toward co-reactant activation strategies, with their consequent reaction networks showing marked departures from those established for supported Ru catalysts. In situ spectroscopy and kinetic analyses propose distinct sensitivity toward CO 2 on Ru-SrTiO 3 systems, where higher CH4 turnover rates result from O-assisted C─H bond activation pathways. These pathways occur on paired Ru-oxygen vacancy sites that are inherent to the perovskite structure and are not readily accessible on supported Ru catalysts. Here, CO 2 is activated on oxygen vacancies directly adjacent to Ru active sites, which facilitates CH 4 C─H bond activation through a surface methoxy intermediate. The high reactivity of Ru-SrTiO 3 systems enables stable CH 4 turnover rates at milder reaction temperatures (673 K), conditions under which the supported counterpart, Ru/SrTiO 3 , is inactive. Overall, this work demonstrates that ligation of catalytically active cations in perovskite oxides facilitates site engineering toward atom-efficient activation of thermodynamically stable feedstocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CRNT4SBML: a Python package for the detection of bistability in biochemical reaction networks

Motivation: Signaling pathways capable of switching between two states are ubiquitous within living organisms. They provide the cells with the means to produce reversible or irreversible decisions. Switchlike behavior of biological systems is realized through biochemical reaction networks capable of having two or more distinct steady states which are dependent on initial conditions. Investigation of whether a certain signaling pathway can confer bistability involves a substantial amount of hypothesis testing. The cost of direct experimental testing can be prohibitive. Therefore, constraining the hypothesis space is highly bene?cial. One such methodology is based on Chemical Reaction Network Theory which uses computational techniques to rule out pathways that are not capable of bistability regardless of kinetic constant values and molecule concentrations. Although useful, these methods are complicated from both pureandcomputationalmathematicsperspectives.Thus,theiradoptionisverylimitedamongstbiologists. Results: We brought Chemical Reaction Network Theory approaches closer to experimental biologists by automating all the necessary steps in CRNT4SMBL. The input is based on SBML format which is the community standard for biological pathway communication. The tool parses SBML and derives C-graph representations of the biological pathway with mass action kinetics. Next steps involve an ef?cient search for potential saddle-node bifurcation points using an optimization technique. This type of bifurcation is important as it has the potential of acting as a switching point between two steady states. Finally, if any bifurcation points are present, numerical continuation analysis extends the equilibria branches for generating the diagram. Presence of an S-shaped bifurcation diagram indicates that the pathway acts as a bistable switch for the given optimization parameters. Availability: CRNT4SBML is available via the Python Package Index. The documentation can be found at https://crnt4sbml.readthedocs.io. CRNT4SBML is licensed under the Apache Software License 2.0. Contact: vladislav.petyuk@pnnl.gov

Reyes, Brandon C.↗

Pyridine Catalysis of Anhydride Hydrolysis within Carbodiimide‐Driven Reaction Networks

Carbodiimide-fueled reaction networks offer a versatile platform for nonequilibrium chemical systems. Typically, the carbodiimide converts a carboxylic acid to its anhydride, called “activation,” which subsequently undergoes hydrolysis, called “deactivation.” Here, we investigate pyridines with variable nucleophilicity as catalysts to control deactivation (pyridine, 4-methylpyridine, 4-methoxypyridine, and 4-dimethylaminopyridine). Reactions have been monitored by NMR spectroscopy. Although this reaction network is simple, determination of well-defined rate constants from kinetic modeling is challenging because of correlation between the parameters. This issue can be addressed by analyzing the anhydride hydrolysis independently. The rate of attack of the pyridines on the anhydride follows expected nucleophilicity trends, although this is offset by increased protonation of more-nucleophilic pyridines at typical pH's. The optimized parameters can be used to model the full carbodiimide-driven process, although the presence of the common carbodiimide EDC has unanticipated effects on the anhydride hydrolysis rate. The results offer context for controlling carbodiimide-fueled reaction networks through the choice of suitable catalysts and pH.

Anhydrides↗

Deep learning with plasma plume image sequences for anomaly detection and prediction of growth kinetics during pulsed laser deposition

Abstract Materials synthesis platforms that are designed for autonomous experimentation are capable of collecting multimodal diagnostic data that can be utilized for feedback to optimize material properties. Pulsed laser deposition (PLD) is emerging as a viable autonomous synthesis tool, and so the need arises to develop machine learning (ML) techniques that are capable of extracting information from in situ diagnostics. Here, we demonstrate that intensified-CCD image sequences of the plasma plume generated during PLD can be used for anomaly detection and the prediction of thin film growth kinetics. We develop multi-output (2 + 1)D convolutional neural network regression models that extract deep features from plume dynamics that not only correlate with the measured chamber pressure and incident laser energy, but more importantly, predict parameters of an auto-catalytic film growth model derived from in situ laser reflectivity experiments. Our results demonstrate how ML with in situ plume diagnostics data in PLD can be utilized to maintain deposition conditions in an optimal regime. Further, the predictive capabilities of plume dynamics on the kinetics of film growth or other film properties prior to deposition provides a means for rapid pre-screening of growth conditions for the non-expert, which promises to accelerate materials optimization with PLD.

36 MATERIALS SCIENCE↗

A mathematical analysis of carbon fixing materials that grow, reinforce, and self-heal from atmospheric carbon dioxide

Carbon fixing materials are a new class of self-healing, self-reinforcing materials we have introduced that utilize ambient CO 2 to chemically add to an ever extending carbon backbone. This class of materials can utilize biological or non-biological photocatalysts and support a wide range of potential backbone chemistries. However, there is no analysis to date that describes their fundamental limits in terms of chemical kinetics and mass transfer. In this computational study, we employ a reaction engineering, and materials science analysis to answer basic questions about the maximum growth rate, photocatalytic requirements and limits of applicable materials. Our proposed mathematical framework envelops three main functions required for carbon fixing materials: (1) adsorption of CO 2 from air, (2) photocatalytic reduction of CO 2 into selected monomers, and (3) polymerization of CO 2 -derived monomers. First, by performing a Damköhler number analysis, we derive criteria for the cross over from kinetic control to mass transfer limited growth, setting upper limits on performance of a potential photocatalyst. Next, we analyze photocatalytic reduction of CO 2 to single carbon products, using known catalytic pathways and kinetic data. We identify formaldehyde as a C1 intermediate having unique potential for incorporation into the material backbone of carbon fixing materials. As an example, we find that a diamond-shaped reaction network graph for CO 2 reduction, passing through CO and HCOOH intermediates, accurately describes kinetic data for cobalt-promoted TiO 2 nanoparticles at room temperature and 1 atm CO 2 . Finally, as an applied case study, we consider and analyze the example of a carbon fixing poly(oxymethylene) system with embedded catalyst promoting the photocatalytic reduction of CO 2 to formaldehyde. The latter then trimerizes to a trioxane monomer which subsequently polymerizes to polyoxymethylene. This reaction engineering analysis introduces benchmarks for carbon fixing materials with respect to achievable rates of photocatalysis, adsorption, and polymerization. These results should prove valuable for the design, evaluation, and benchmarking of this emergent and new class of environmentally sustainable materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine Learning Approach to Modeling of Neutral Particles Transport in Plasma

A propagator‐based approach is investigated for Monte‐Carlo (MC) modeling of neutral particle transport in fusion boundary plasmas. The propagator is based on a Green's function for the neutral kinetic equation, which depends on the plasma profiles. A neural network (NN)‐based model for the propagator provides a fast and accurate solution for the neutral distribution function in plasma. Preliminary results from a small 1D test problem look encouraging. The proposed approach, a propagator‐based NN model for neutral transport in plasma, has potential for generalization to higher dimensions and efficient coupling with plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Synergistic Coupling Effect of Electronic Conductivity and Interphase Compatibility on High-Voltage Na 3 V 2 (PO 4 ) 2 F 3 Cathodes

Na 3 V 2 (PO 4 ) 2 F 3 (NVPF) has been considered an up-and-coming cathode material candidate for sodium (Na) ion batteries in light of its high specific capacity and working voltage. However, an erratic cathode/electrolyte interface layer is inevitably formed, accompanied by continuous electrolyte decomposition on the NVPF surface, when the voltage exceeds 4.2 V vs Na + /Na. Herein, the interphase features of NVPF are obviously enhanced owing to the ameliorated electronic conductivity obtained by combining it with carbon nanotubes (CNT). The NVPF with 3 wt % CNT (NVPF@3% CNT) reduces the Na + diffusion kinetic energy barrier and electron transport resistance. Furthermore, the conducting network formed by CNT with sturdy structure strength can promptly accommodate the volumetric changes during sequential Na + extraction/insertion and thus effectively improve the long-term cyclic performance of NVPF/hard carbon full cells. The initial discharge capacity approaches 105 mA h g –1 at 0.5C, and it retains 94% capacity retention after 200 cycles at the temperature of –10 °C. The cathode/electrolyte interphase characterization results further demonstrate that the interphase layer on the NVPF@3% CNT cathode is thinner and more compact compared with pristine samples. Here, this research provides a competitive strategy to facilitate the interfacial compatibility between the NVPF and electrolytes and accelerate the commercialization of high-performance Na-ion batteries.

25 ENERGY STORAGE↗

Tandem synthesis of tetrahydroquinolines and identification of the reaction network by operando NMR

The study of the reaction mechanism and complex network for heterogeneously catalyzed tandem reactions is challenging but can guide reaction design and optimization. Here, we describe a case study using bifunctional metal–organic framework supported Pd nanoparticles (Pd/UiO-66(HCl)) for the one-pot tandem synthesis of substituted tetrahydroquinolines via the Claisen–Schmidt condensation and reductive intramolecular cyclization. The directly observed evolution of intermediates and products, including reactive species containing hydroxylamine group and unstable intermediate 2-phenyl-3,4-dihydroquinoline, was enabled by operando magic angle spinning nuclear magnetic resonance studies under 50 bar H 2 . The reaction network of the tandem reaction is deduced based on reaction kinetic information obtained from the operando study. The optimized procedure was applied to various acetophenone and nitrobenzaldehyde derivatives carrying different functional groups, and eight valuable substituted tetrahydroquinolines were obtained in moderate to good yields. This work provides a molecular-level understanding of the catalytic system and brings up new opportunities for efficient and sustainable synthesis of medicinally relevant building blocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Symbolic regression development of empirical equations for diffusion in Lennard-Jones fluids

Symbolic regression (SR) with a multi-gene genetic program has been used to elucidate new empirical equations describing diffusion in Lennard-Jones (LJ) fluids. Some examples include equations to predict self-diffusion in pure LJ fluids and equations describing the finite-size correction for self-diffusion in binary LJ fluids. The performance of the SR-obtained equations was compared to that of both the existing empirical equations in the literature and to the results from artificial neural net (ANN) models recently reported. It is found that the SR equations have improved predictive performance in comparison to the existing empirical equations, even though employing a smaller number of adjustable parameters, but show an overall reduced performance in comparison to more extensive ANNs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A physics-informed deep learning description of Knudsen layer reactivity reduction

A physics-informed neural network (PINN) is used to evaluate the fast ion distribution in the hot spot of an inertial confinement fusion target. The use of tailored input and output layers to the neural network is shown to enable a PINN to learn the parametric solution to the Vlasov–Fokker–Planck equation in the absence of any synthetic or experimental data. As an explicit demonstration of the approach, the specific problem of Knudsen layer fusion yield reduction is treated. Here, the predictions from the Vlasov–Fokker–Planck PINN are used to provide a non-perturbative solution of the fast ion tail in the vicinity of the hot spot, thus allowing the spatial profile of the fusion reactivity to be evaluated for a range of collisionalities and hot spot conditions. Excellent agreement is found between the predictions of the Vlasov–Fokker–Planck PINN and the results from traditional numerical solvers with respect to both the energy and spatial distribution of fast ions and the fusion reactivity profile, demonstrating that the Vlasov–Fokker–Planck PINN provides an accurate and efficient means of determining the impact of Knudsen layer yield reduction across a broad range of plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling Combustion Reaction ODEs with Neural Networks

The chemistry of combustion reactions is complex as it involves many species and reactions. In practice, such a system is often modeled computationally using an empirically derived chemical mechanism. Given the large range of reaction rates, the system is then time evolved using a stiff ODE solver. However, even when reduced chemical mechanisms are employed, the system can become computationally expensive for two-dimensional or three-dimensional systems. Such systems can also face problems with instability. As such, it is desirable to find a cheaper, stable alternative to solving the reaction system. Neural networks offer the potential to learn these reaction ODEs and time evolve a combustion reaction in a more cost-efficient manner than stiff ODE solvers. In this study, a variety of neural networks are trained on zero-dimensional Cantera simulations of methane combustion with varying initial conditions. Several predictive approaches as well as several neural network architectures (artificial neural network with dropout, ResNet, and Neural ODE) are compared in their ability to predict combustion trajectories. Promising models are then identified.

combustion kinetics↗

Molecular Dynamics Simulation and Theoretical Analysis of Structural Relaxation, Bond Exchange Dynamics, and Glass Transition in Vitrimers

Vitrimers are a class of polymer networks featuring dynamic covalent cross-links that can undergo associative bond exchange. There has been recent interest in these materials due to their promise as recyclable thermosets or self-healing polymers because of the ability of vitrimer networks to rearrange at the molecular level and undergo macroscopic flow. However, the practical use of these materials often occurs in the supercooled regime or glassy state, where the implications of dynamic bonds are complicated by the interplay between slow activated segmental dynamics, cross-link (i.e., bond-exchange) kinetics, and ultimately material properties. Here, in this paper, we combine coarse-grained molecular dynamics simulation and microscopic statistical mechanical theory to understand how cross-linking kinetics affect material dynamics and how this couples to segmental relaxation of the polymeric network strands across a spectrum of length and time scales, especially in the supercooled regime. We characterize the Kuhn segmental alpha relaxation time and bond exchange time for vitrimer systems across various cross-link densities, temperatures, and bond exchange rates. Simulation and theory both exhibit a bending-up behavior for bond exchange time upon cooling, suggesting a coupling between bond exchange dynamics and segmental relaxation that intensifies with faster bond exchange kinetics. We also found bond exchange dynamics have an impact on Kuhn segment alpha relaxation time, which is most significant at higher cross-link densities. Both these effects are most prominent when the bond exchange time is similar to the Kuhn segment alpha relaxation time, and the resulting coupling of these two relaxation processes is tied to both the probability of a free end to find a bonded pair and the time scale of the constraints imposed by the dynamic cross-links. This relationship is reflected by a cross-link dependence of a theoretical parameter which represents the quantitative degree of coupling between bond exchange and segmental dynamics. Overall, the combination of simulation and theory clarifies the intricate interaction between bond kinetics and segmental relaxation and demonstrates the ability to provide molecular-level insights into vitrimer dynamics over a wide temperature range.

dynamic relaxation↗