Search NASA⌕ Search

SEARCH · Search NASA

Results for “reaction 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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Unravelling the mechanism of pressure induced polyamorphic transition in an inorganic molecular glass

The atomic structure of a germanium doped phosphorous selenide glass of composition Ge 2.8 P 57.7 Se 39.5 is determined as a function of pressure from ambient to 24 GPa using Monte-Carlo simulations constrained by high energy x-ray scattering data. The ambient pressure structure consists primarily of P 4 Se 3 molecules and planar edge shared phosphorus rings, reminiscent of those found in red phosphorous as well as a small fraction of locally clustered corner-sharing GeSe 4 tetrahedra. This low-density amorphous phase transforms into a high-density amorphous phase at ~6.3 GPa. The high-pressure phase is characterized by an extended network structure. The polyamorphic transformation between these two phases involves opening of the P 3 ring at the base of the P 4 Se 3 molecules and subsequent reaction with red phosphorus type moieties to produce a cross linked structure. The compression mechanism of the low-density phase involves increased molecular packing, whereas that of the high pressure phase involves an increase in the nearest-neighbor coordination number while the bond angle distributions broaden and shift to smaller angles. The entropy and volume changes associated with this polyamorphic transformation are positive and negative, respectively, and consequently the corresponding Clapeyron slope for this transition would be negative. This result has far reaching implications in our current understanding of the thermodynamics of polyamorphic transitions in glasses and glass-forming liquids.

42 ENGINEERING↗
DOE OSTI2023

Learning effective stochastic differential equations from microscopic simulations: Linking stochastic numerics to deep learning

We identify effective stochastic differential equations (SDEs) for coarse observables of fine-grained particle- or agent-based simulations; these SDEs then provide useful coarse surrogate models of the fine scale dynamics. We approximate the drift and diffusivity functions in these effective SDEs through neural networks, which can be thought of as effective stochastic ResNets. The loss function is inspired by, and embodies, the structure of established stochastic numerical integrators (here, Euler–Maruyama and Milstein); our approximations can thus benefit from backward error analysis of these underlying numerical schemes. They also lend themselves naturally to “physics-informed” gray-box identification when approximate coarse models, such as mean field equations, are available. Existing numerical integration schemes for Langevin-type equations and for stochastic partial differential equations can also be used for training; we demonstrate this on a stochastically forced oscillator and the stochastic wave equation. Our approach does not require long trajectories, works on scattered snapshot data, and is designed to naturally handle different time steps per snapshot. We consider both the case where the coarse collective observables are known in advance, as well as the case where they must be found in a data-driven manner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Surrogate models to optimize plasma-assisted atomic layer deposition in high aspect ratio features

In this work, we explore surrogate models to optimize plasma enhanced atomic layer deposition (PEALD) in high aspect ratio features. In plasma-based processes such as PEALD and atomic layer etching (ALE), surface recombination can dominate the reactivity of plasma species with the surface, which can lead to unfeasibly long exposure times to achieve full conformality inside nanostructures like high aspect ratio vias. Using a synthetic dataset based on simulations of PEALD, we train artificial neural networks to predict saturation times based on cross section thickness data obtained for partially coated conditions. The results obtained show that just two experiments in undersaturated conditions contain enough information to predict saturation times within 10% of the ground truth. A surrogate model trained to determine whether surface recombination dominates the plasma–surface interactions in a PEALD process achieves 99% accuracy. This demonstrates that machine learning can provide a new pathway to accelerate the optimization of PEALD processes in areas such as microelectronics. Our approach can be easily extended to ALE and more complex structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Impact of aqueous solution pH on network structure of corrosion–induced surface layers of boroaluminosilicate glass

Understanding of the long-term stability of vitrified nuclear waste glass is critical to its safe storage, but knowledge gaps remain regarding the corrosion mechanism, especially with respect to surface alteration layers whose presence changes the observed corrosion rates. Here, we report the aqueous corrosion of a model boroaluminosilicate glass at ground water relevant pHs as investigated by vibrational spectroscopy to understand the network structure of the alteration layer that we hypothesize controls the corrosion rates. Both specular reflection infrared spectroscopy and infrared variable angle spectroscopic ellipsometry were used to elucidate details of the structure of the alteration layers formed at pH 7 and pH 9. After 2312 days, the thickness of the alteration layer is found to be ~2000 nm for corrosion at pH 7 and ~300 nm for corrosion at pH 9. Here, the interpretation of the vibrational spectra suggests that alteration layers formed at pH 7 and pH 9 both exhibit a porous network with "silica-like" structures, but with differences in the bond lengths and bond angles in the glass network. Moreover, the relative abundance of the hydrous species within the porous network appears to be dependent on pH. Additionally, the spectroscopic data suggest that the B and Ca trapped in the alteration layer during the secondary corrosion of the pH 9 pre-corroded sample in pH 3 solution are hydrated within the nano-porous layer and not bound to the silicate network. However, the precise chemical form remains to be determined. Based on these findings, we propose that the combination of the network structure of the alteration layer and the relative abundance of hydrous species within the alteration layer act to control the transport of ionic species dissolved from the reaction front to the aqueous solution outside the glass.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Jeankempite, Ca 5 (AsO 4 ) 2 (AsO 3 OH) 2 (H 2 O) 7 , a new arsenate mineral from the Mohawk Mine, Keweenaw County, Michigan, USA

Abstract Jeankempite, Ca 5 (AsO 4 ) 2 (AsO 3 OH) 2 (H 2 O) 7 , is a new mineral species (IMA2018-090) discovered amongst coatings of arsenate minerals on oxidised copper arsenides from the Mohawk No. 2 mine, Mohawk, Keweenaw County, Michigan, USA. The new mineral occurs as lamellar bundles of colourless to white plates up to 1 mm wide and is visually indistinguishable from guérinite, with which it forms intergrowths. Jeankempite is transparent to translucent with a waxy lustre and white streak, is non-fluorescent under longwave and shortwave ultraviolet illumination, has a Mohs hardness of ~1.5 and brittle tenacity with uneven fracture. Crystals are flattened on {01 $\bar{1}$ } and exhibit perfect cleavage on {01 $\bar{1}$ }. Optically, jeankempite is biaxial (+), α = 1.601(2), β = 1.607(2), γ = 1.619(2) (white light); 2V meas. = 72(2)° and 2V calc. = 71.0°. The empirical formula is (Ca 4.97 Na 0.013 Mg 0.017 )(As 3.99 S 0.01 ) 4 O 23 H 16 , based on 23 O and 16 H atoms per formula unit. Thermogravimetric analysis indicates that jeankempite undergoes four weight losses totalling 16.82%, close to the expected loss of 16.30%, corresponding to eight H 2 O. Jeankempite is triclinic, P $\bar{1}$ , a = 6.710(6), b = 14.901(14), c = 15.940(15) Å, α = 73.583(12)°, β = 81.984(12)°, γ = 82.754(12)°, V = 1507(2) Å 3 and Z = 3. The final structure was refined to R 1 = 0.0591 for 2781 reflections with I obs > 3σ I . The crystal structure of jeankempite is built from a network of edge- and vertex-sharing CaO 6 , CaO 7 and AsO 4 polyhedra, and we hypothesise that the new mineral has formed due to a topotactic reaction brought on by dehydration of preexisting guérinite.

Mineralogy↗
NASA NTRS2003

Numerical Modeling of Fluid Transient in Cryogenic Fluid Network of Rocket Propulsion System

Fluid transients, also known as water hammer, can have a significant impact on the design and operation of both spacecraft and launch vehicles propulsion systems. These transients often occur at system activation and shut down. For ground safety reasons, many spacecrafts are launched with the propellant lines dry. These lines are often evacuated by the time the spacecraft reaches orbit. When the propellant isolation valve opens during propulsion system activation, propellant rushes into lines creating a pressure surge. During propellant system shutdown, a pressure surge is created due to sudden closure of a valve. During both activation and shutdown, pressure surges must be predicted accurately to ensure structural integrity of the propulsion system fluid network. The method of characteristics is the most widely used method of calculating fluid transients in pipeline [ 1,2]. The method of characteristics, however, has limited applications in calculating flow distribution in complex flow circuits with phase change, heat transfer and rotational effects. A robust cryogenic propulsion system analyzer must have the capability to handle phase change, heat transfer, chemical reaction, rotational effects and fluid transients in conjunction with subsystem flow model for pumps, valves and various pipe fittings. In recent years, such a task has been undertaken at Marshall Space Flight Center with the development of the Generalized Fluid System Simulation Program (GFSSP), which is based on finite volume method in fluid network [3]. GFSSP has been extensively verified and validated by comparing its predictions with test data and other numerical methods for various applications such as internal flow of turbo-pump [4], propellant tank pressurization [5,6], chilldown of cryogenic transfer line [7] and squeeze film damper rotordynamics [8]. The purpose of the present paper is to investigate the applicability of the finite volume method to predict fluid transient in cryogenic flow circuits.

Majumdar, Alok↗
DOE OSTI2021

Carbon Fiber Reinforced Epoxy Vitrimer: Robust Mechanical Performance and Facile Hydrothermal Decomposition in Pure Water

Conventional carbon fiber reinforced thermosetting polymers (CFRPs) are neither recyclable nor repairable due to their crosslinked network. The rapid growing CFRP market raises a serious concern of the waste management. In this work, a viable method to develop a readily recyclable CFRP based on epoxy vitrimer is introduced. First, a self-catalytic epoxy prepolymer with built-in hydroxy and tertiary amine groups is designed, which upon reaction with an anhydride formed a catalyst-free epoxy vitrimer. The epoxy prepolymer is synthesized from a diamine and an excess of bisphenol A epoxy resin. The hydroxyls and tertiary amines of the epoxy prepolymer efficiently catalyze both curing and the dynamic transesterification of the crosslinked polymer without the need of a catalyst. Then, the epoxy vitrimer is used as the matrix resin to prepare CFRP. The resulting CFRP exhibited a tensile strength as high as 356 MPa. More interestingly, the matrix of the CFRP is efficiently degraded in pure water at above 160 °C. This is because the built-in tertiary amines catalyze the hydrolysis of the ester bonds of the crosslinked network. The simple method developed in this work provides a framework for the development of recyclable CFRP.

Liu, Tuan↗
DOE OSTI2024

Quantum Simulations of Radiation Damage in a Molecular Polyethylene Analog

Abstract An atomic‐level understanding of radiation‐induced damage in simple polymers like polyethylene is essential for determining how these chemical changes can alter the physical and mechanical properties of important technological materials such as plastics. Ensembles of quantum simulations of radiation damage in a polyethylene analog are performed using the Density Functional Tight Binding method to help bind its radiolysis and subsequent degradation as a function of radiation dose. Chemical degradation products are categorized with a graph theory approach, and occurrence rates of unsaturated carbon bond formation, crosslinking, cycle formation, chain scission reactions, and out‐gassing products are computed. Statistical correlations between product pairs show significant correlations between chain scission reactions, unsaturated carbon bond formation, and out‐gassing products, though these correlations decrease with increasing atom recoil energy. The results present relatively simple chemical descriptors as possible indications of network rearrangements in the middle range of excitation energies. Ultimately, the work provides a computational framework for determining the coupling between nonequilibrium chemistry in polymers and potential changes to macro‐scale properties that can aid in the interpretation of future radiation damage experiments on plastic materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Maize specialized metabolome networks reveal organ-preferential mixed glycosides

Despite the scientific and economic importance of maize, little is known about its specialized metabolism. Here, five maize organs were profiled using different reversed-phase liquid chromatography-mass spectrometry methods. The resulting spectral metadata, combined with candidate substrate-product pair (CSPP) networks, allowed the structural characterization of 427 of the 5,420 profiled compounds, including phenylpropanoids, flavonoids, benzoxazinoids, and auxin-related compounds, among others. Only 75 of the 427 compounds were already described in maize. Analysis of the CSPP networks showed that phenylpropanoids are present in all organs, whereas other metabolic classes are rather organ-enriched. Frequently occurring CSPP mass differences often corresponded with glycosyl- and acyltransferase reactions. The interplay of glycosylations and acylations yields a wide variety of mixed glycosides, bearing substructures corresponding to the different biochemical classes. For example, in the tassel, many phenylpropanoid and flavonoid-bearing glycosides also contain auxin-derived moieties. The characterized compounds and mass differences are an important step forward in metabolic pathway discovery and systems biology research.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A Brief Overview of Reactive Transport Codes Used in CO 2 Applications

Reactive transport modelling consists of computational and numerical models that describe the coupled physical, chemical, mechanical, and biological processes interacting with each other over a broad range of spatial and temporal scales. This report reviews some of the reactive transport codes available in the literature for carbon dioxide (CO 2 ) applications and can assist the scientific community with the applicability of each code to solve a specific problem. Approaches implemented in these codes can be categorized in three groups: (i) continuum-scale, (ii) pore-scale, and (iii) hybrid-scale approaches. Theoretical foundations, numerical implementations, and application examples using the models are described. In the continuum approach, flow and transport are formulated in terms of a representative control volume of the medium and coupled with geochemical reactions. The medium is characterized by bulk parameters such as porosity, permeability, or reactive surface area. In the pore-scale approach, each point of space in the pore network is occupied by either a fluid or solid phase. The pore-scale approach requires an exact knowledge of the spatial and temporal phase distribution. In the hybrid-scale approach, media are described by multiple characteristic length-scales, with some regions using a pore-scale approach while others used a continuum approach.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Astromers: Nuclear Isomers in Astrophysics

In this work, we develop a method to compute thermally mediated transition rates between the ground state and long-lived isomers in nuclei. We also establish criteria delimiting a thermalization temperature above which a nucleus may be considered a single species and below which it must be treated as two separate species: a ground-state species and an astrophysical isomer ("astromer") species. Below the thermalization temperature, the destruction rates dominate the internal transition rates between the ground state and the isomer. If the destruction rates also differ greatly from one another, the nuclear levels fall out of or fail to reach thermal equilibrium. Without thermal equilibrium, there may not be a safe assumption about the distribution of occupation probability among the nuclear levels when computing nuclear reaction rates. In these conditions, the isomer has astrophysical consequences and should be treated as a separate astromer species which evolves separately from the ground state in a nucleosynthesis network. We apply our transition-rate methods and perform sensitivity studies on a few well-known astromers. We also study transitions in several other isomers of likely astrophysical interest.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Rule-based simulation models

Procedural modeling systems, rule based modeling systems, and a method for converting a procedural model to a rule based model are described. Simulation models are used to represent real time engineering systems. A real time system can be represented by a set of equations or functions connected so that they perform in the same manner as the actual system. Most modeling system languages are based on FORTRAN or some other procedural language. Therefore, they must be enhanced with a reaction capability. Rule based systems are reactive by definition. Once the engineering system has been decomposed into a set of calculations using only basic algebraic unary operations, a knowledge network of calculations and functions can be constructed. The knowledge network required by a rule based system can be generated by a knowledge acquisition tool or a source level compiler. The compiler would take an existing model source file, a syntax template, and a symbol table and generate the knowledge network. Thus, existing procedural models can be translated and executed by a rule based system. Neural models can be provide the high capacity data manipulation required by the most complex real time models.

Nieten, Joseph L.↗
DOE OSTI2021

Ecosystems and Networks Integrated with Genes and Molecular Assemblies (ENIGMA) (Final Scientific/Technical Report, Subcontract Award No. 6953691)

The overarching objective the Baliga lab is to develop a mechanistic understanding of the field relevant metabolic processes underlying community partitioning of respiration pathways and physiological state shifts. In order to dissect how environmental perturbations influence microbial regulatory networks, we have constructed data-driven models that capture the dynamic interplay between abiotic and biotic factors during laboratory simulations. Specifically, we have developed multiple computational models (i.e., cMonkey2, EGRIN 2.0, and multiple genome-scale metabolic models) connecting environmental influence to genome-encoded regulatory programs. In addition, we have generated new tools such as Live Anaerobic Cell Sorting (LAnCS) and Fluidized Bed Reactors (FBRs) to investigate the dynamics of community partitioning of metabolic processes supporting sulfate or nitrate respiration, which are two critical activities observed from the field site at Oak Ridge National Labs (ORNL). This framework of connecting field phenomena with laboratory simulations and computational modeling has created the Environmental Simulations and Modelling (EnvSim) campaign within ENIGMA. The campaign, which spans multiple laboratories across ENIGMA and is led by the Baliga group, has established synthetic communities (SynComs) to discover, characterize, and dissect key microbial processes relevant to field observations, such as emissions of the greenhouse gas nitrous oxide (N 2 O). For example, the EnvSim campaign so far has deduced four potential mechanisms that may account for the N 2 O emissions at the FRC and are currently being investigated by labs across ENIGMA. Denitrification may be driven by complete denitrifiers, however, their NosZ enzymes, which catalyze the final step in denitrification by converting N 2 O to N 2 , may be sensitive at a lower pH. Additionally, the denitrification process could be partitioned among organisms and some may have pH-sensitive NosZ genes. Another factor contributing to variable N 2 O emissions relates to the metal co-factors involved in the denitrification pathway. For instance, excess Cu, Al, Mn, U, Ni, Co, Cu, and/or Cd may have inhibitory effects on multiple enzymatic steps during denitrification, while the enzymatic production of nitrite, the precursor of N 2 O, can be limited by the essential metal Mo. Lastly, abiotic production of N 2 O may occur as a result of chemodenitrification, in which metals like Fe, Mn, and some organic compounds can drive redox reactions that convert nitrogen cycle intermediates to N 2 O under the right conditions. Current investigations and analyses have generated multiple transcriptomic profiles aiding the refinement of new metabolic and gene regulatory network model for our denitrifying SynCom. In sum, the previous funding cycle has generated 33 peer-reviewed publications spanning predictive network biology to co-operativity of mutualistic interspecies interactions of evolved communities. Our integration of multiple omics datasets and the development of new network modelling tools has revealed important insights that can explain field observed phenomenon and generated new hypotheses to be investigated in the field.

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

Effective many-body interactions in reduced-dimensionality spaces through neural network models

Accurately describing properties of challenging problems in physical sciences often requires complex mathematical models that are unmanageable to tackle head on. Therefore, developing reduced-dimensionality representations that encapsulate complex correlation effects in many-body systems is crucial to advance the understanding of these complicated problems. However, a numerical evaluation of these predictive models can still be associated with a significant computational overhead. To address this challenge, in this paper we discuss a combined framework that integrates recent advances in the development of active-space representations of coupled cluster (CC) downfolded Hamiltonians with neural network approaches. The primary objective of this effort is to train neural networks to eliminate the computationally expensive steps required for evaluating hundreds or thousands of Hugenholtz diagrams, which correspond to multidimensional tensor contractions necessary for evaluating a many-body form of downfolded effective Hamiltonians. Using small molecular systems (the H 2 O and HF molecules) as examples, we demonstrate that training neural networks employing effective Hamiltonians for a few nuclear geometries of molecules can accurately interpolate or extrapolate their forms to other geometrical configurations characterized by different intensities of correlation effects. We also discuss differences between effective interactions that define CC downfolded Hamiltonians with those of bare Hamiltonians defined by Coulomb interactions in the active spaces. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Deep learning-enabled probing of irradiation-induced defects in time-series micrographs

Modeling time-series data with convolutional neural networks (CNNs) requires building a model to learn in batches as opposed to training sequentially. Coupling CNNs with in situ or operando techniques opens the possibility of accurately segmenting dynamic reactions and mass transport phenomena to understand how materials behave under the conditions in which they are used. In this article, in situ ion irradiation transmission electron microscopy (TEM) images are used as inputs into the CNN to assess the defect generation rate, defect cluster density, and saturation of defects. We then use the output segmentation maps to correlate with conventional TEM micrographs to assess the model’s ability to detail nanoscale interactions. Next, we discuss the implications of preprocessing and hyperparameters on model variability, accuracy when expanded to other datasets, and the role of regularization when controlling model variance. Ultimately, we eliminate human bias when extrapolating physical metrics, speed up analysis time, decouple reactions that happen at 100 ms intervals, and deploy models that are both accurate and transferable to similar experiments.

Burns, Kory (ORCID:0000000158019909)↗
DOE OSTI2023

Grid-Integrated Production of Fischer-Tropsch Synfuels from Nuclear Power

Idaho National Laboratory (INL) investigates the relative economic profitability of an integrated energy system (IES) coupling an NPP with a synfuel production process at selected case study locations across the United States. In the synfuel IES, a high-temperature steam electrolysis (HTSE) plant is thermally and electrically coupled with an NPP to produce zero-carbon hydrogen. The synthetic fuel is produced from this H 2 combined with a CO 2 supply using the reverse water gas shift process followed by the Fischer-Tropsch (FT) reaction. This analysis considers a system in which the CO 2 is sourced from regional CO 2 emitters via the construction and operation of pipeline-based CO 2 supply networks. Locating the FT plant at the same site as the NPP and HTSE plants enables the NPP to provide zero-carbon heat and power to the HTSE plant and zero-carbon power to the FT plant as well as avoid the requirement for long-distance H 2 product transport from the HTSE plant to the FT plant. Hydrogen storage is used to enable the NPP to dispatch power to the electrical grid (instead of the HTSE plant) when grid demand increases, thus enabling the FT plant to continue to operate in a steady-state production mode. The ability to cease hydrogen production for several hours within each day enables the NPP to provide power to the grid to balance the electricity market during peak periods and maximize revenues for the nuclear synfuel IES. The FT process design considered has a 99% carbon conversion efficiency. The use of nuclear energy and nuclear energy-derived hydrogen enables synfuel production to achieve this high level of carbon utilization. Additionally, the life-cycle carbon emissions of the nuclear-based synfuel production process are very low, with WTW emissions of approximately 25 gCO 2 e/MJ, including steam credits (generated from FT process excess heat), and approximately 7 gCO 2 e/MJ, if steam credits are excluded. This compares favorably with the WTW emissions of 90.5 gCO 2 e/MJ for a compression-ignition, direct injection (CIDI) vehicle with a fuel economy of 31.6 miles per gallon gasoline equivalent (MPGGE), using low-sulfur diesel produced using conventional petroleum production and refining processes. Several NPPs in various regions of the U.S. are considered as case study analyses. Supply locations and transportation via pipeline of the CO 2 feedstock to the NPP site are analyzed through the National Energy Technology Laboratory (NETL) CO 2 Transport Cost model. The team finds that the amount of CO 2 generated by different sectors is sufficient for the synfuel production process at all locations considered. The CO 2 transportation costs are functions of the distance of the source to the NPP location, the CO 2 capture cost at the source, and the quantity of CO 2 transported. Historical electricity prices for the NPP case study locations are collected and analyzed. Monthly average prices, price range, and duration of negative-price periods vary among these locations. For each location, an auto-regressive moving average (ARMA) model is trained on historical electricity price data. ARMA validation is done to ensure the synthetic price distributions represent one of historical prices with high fidelity. Synthetic time series from these ARMA models are used in a coupled dispatch and system optimization in the Holistic Energy Resource Optimization Network (HERON) to compute the differential net present value (NPV) of the IES. The team finds that this econometric is positive, ranging from $14M–1.3bn (2020) depending on the location. The optimal synfuel IES configuration to obtain this increase in NPV often maximizes the size of the synfuel production process with regards to the size of the NPP. However, the team shows that the NPP still plays a stabilizing role for the grid: In periods of high prices and high loads, more electricity from the NPP is sent to the grid. A high variability of electricity prices and extreme maximum prices tend to drive up electricity production. While it requires significant investment, the synfuel IES could increase the economic profitability for the existing fleet of LWRs across the country while still maintaining the grid stabilizer role of NPPs. During its lifetime, the main costs for the nuclear synfuel IES are the carbon feedstock transportation costs, followed by the capital expenses (CAPEX) and operation and maintenance (O&M) costs while the revenue comes first from the IRA H 2 production tax credit (PTC) and then from the sales of synfuel products. The profitability of the synfuel IES is most sensitive to the value of the hydrogen PTC and the synfuel products as well as the cost of the carbon feedstock, highlighting the importance of governmental incentives regarding hydrogen, carbon emissions, and synfuel in driving the deployment of future nuclear synfuel IESs.

08 HYDROGEN↗
NASA NTRSDate not supplied

Building and Breaking Carbon Composites with REACTER

Carbon composites have become indispensable for aerospace and other high-performance applications, and a detailed picture of their morphology and failure mechanisms remains difficult to obtain through experiment. REACTER is a versatile computational modeling tool for atomistic molecular dynamics designed to model chemical reactions at the speed and length scales of classical force fields.1 In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites, carbon-fiber reinforced polymers (CFRP) and carbon nanotube (CNT) composites. Carbon fiber core morphologies were created by the method of Desai et al.,2 but using the advanced reaction constraints framework of REACTER, their proposed multistep procedure was reduced to a single uninterrupted molecular dynamics simulation. The carbon fiber filler was embedded into a polymer matrix by simulated in situ polymerization of several thermosetting resins, including bismaleimide and polyarylacetylene, to obtain the final CFRP model. To generate the second class of carbon composite, CNT networks were grown dynamically using the new ‘create atoms’ feature of REACTER, and similarly infiltrated with resin to obtain CNT composites. The resulting models were compared directly to experiment using simulated high resolution transmission electron microscopy and x-ray diffraction. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a reaction constraint on the total potential energy of the involved atoms.

Molecular Dynamics↗
DOE OSTI2023

MetaboDirect: an analytical pipeline for the processing of FT-ICR MS-based metabolomic data

Background: Microbiomes are now recognized as the main drivers of ecosystem function ranging from the oceans and soils to humans and bioreactors. However, a grand challenge in microbiome science is to characterize and quantify the chemical currencies of organic matter (i.e., metabolites) that microbes respond to and alter. Critical to this has been the development of Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS), which has drastically increased molecular characterization of complex organic matter samples, but challenges users with hundreds of millions of data points where readily available, user-friendly, and customizable software tools are lacking. Results: Here, we build on years of analytical experience with diverse sample types to develop MetaboDirect, an open-source, command-line-based pipeline for the analysis (e.g., chemodiversity analysis, multivariate statistics), visualization (e.g., Van Krevelen diagrams, elemental and molecular class composition plots), and presentation of direct injection high-resolution FT-ICR MS data sets after molecular formula assignment has been performed. When compared to other available FT-ICR MS software, MetaboDirect is superior in that it requires a single line of code to launch a fully automated framework for the generation and visualization of a wide range of plots, with minimal coding experience required. Among the tools evaluated, MetaboDirect is also uniquely able to automatically generate biochemical transformation networks (ab initio) based on mass differences (mass difference network-based approach) that provide an experimental assessment of metabolite connections within a given sample or a complex metabolic system, thereby providing important information about the nature of the samples and the set of microbial reactions or pathways that gave rise to them. Finally, for more experienced users, MetaboDirect allows users to customize plots, outputs, and analyses. Conclusion: Application of MetaboDirect to FT-ICR MS-based metabolomic data sets from a marine phage-bacterial infection experiment and a Sphagnum leachate microbiome incubation experiment showcase the exploration capabilities of the pipeline that will enable the research community to evaluate and interpret their data in greater depth and in less time. It will further advance our knowledge of how microbial communities influence and are influenced by the chemical makeup of the surrounding system. The source code and User’s guide of MetaboDirect are freely available through (https://github.com/Coayala/MetaboDirect) and (https://metabodirect.readthedocs.io/en/latest/), respectively.

54 ENVIRONMENTAL SCIENCES↗