SEARCH · Search NASA
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.
Prediction of bond dissociation energies for charged molecules with a graph neural network enabled by global molecular features and reaction difference features between products and reactants.
Concrete gains its strength from the precipitation of a calcium–alumino–silicate–hydrate (C–A–S–H) colloidal gel, which acts as its binding phase. However, despite concrete’s ubiquity in the building environment, the atomic-scale mechanism of C–A–S–H precipitation is still unclear. Here, we use reactive molecular dynamics simulations to model the early-age precipitation of a C–A–S–H gel. We find that, upon gelation, silicate and aluminate precursors condensate and polymerize to form an aluminosilicate gel network. Notably, we demonstrate that the gelation reaction is driven by the existence of a mismatch of atomic-level internal stress between Si and Al polytopes, which are initially experiencing some local tension and compression, respectively. The polymerization of Si and Al polytopes enables the release of these competitive stresses.
Here we present the analysis of approximation rates of operator learning in Chen and Chen (1995) and Lu et al. (2021), where continuous operators are approximated by a sum of products of branch and trunk networks. In this work, we consider the rates of learning solution operators from both linear and nonlinear advection–diffusion equations with or without reaction. We find that the approximation rates depend on the architecture of branch networks as well as the smoothness of inputs and outputs of solution operators.
Thermoplastic polymers have many desirable properties for consumer applications and are complemented by efficient thermal processing techniques, reducing the cost of manufacturing. Lignin exists as an immense biobased carbon source but has largely been researched for its use in thermoset materials due to its own cross-linked, polyfunctional nature. In this study, a new reaction design is employed to create a thermoplastic polyamide network incorporating lignin that is tested to be 99% biobased carbon by radiocarbon analysis. Chemical analysis reveals the nature of lignin incorporation based on chain extension and cross-linking models. The thermal and rheological properties of the new polymers are thoroughly investigated to demonstrate the higher melt-strength capability of the lignin-based polymers facilitating their use in modern processing equipment. This analysis results in finding an optimal lignin loading ratio in the polymer composition reflected by improved tensile strength and stiffness. The results point to a promising polymer design for applying industrial kraft lignin in high-value thermoplastic polymer applications.
Graph layout algorithms used in network visualization represent the first and the most widely used tool to unveil the inner structure and the behavior of complex networks. Current network visualization software relies on the force-directed layout (FDL) algorithm, whose high computational complexity makes the visualization of large real networks computationally prohibitive and traps large graphs into high energy configurations, resulting in hard-to-interpret “hairball” layouts. Here we use Graph Neural Networks (GNN) to accelerate FDL, showing that deep learning can address both limitations of FDL: it offers a 10 to 100 fold improvement in speed while also yielding layouts which are more informative. We analytically derive the speedup offered by GNN, relating it to the number of outliers in the eigenspectrum of the adjacency matrix, predicting that GNNs are particularly effective for networks with communities and local regularities. Finally, we use GNN to generate a three-dimensional layout of the Internet, and introduce additional measures to assess the layout quality and its interpretability, exploring the algorithm’s ability to separate communities and the link-length distribution. The novel use of deep neural networks can help accelerate other network-based optimization problems as well, with applications from reaction-diffusion systems to epidemics.
Microporous polymers feature shape-persistent free volume elements (FVEs), which are permeated by small molecules and ions when used as membranes for chemical separations, water purification, fuel cells, and batteries. It remains a significant challenge to identify FVEs with analyte specificity, due to difficulties in generating microporous polymer libraries with sufficient diversity for screening their properties. Here, we describe a diversity-oriented synthetic (DOS) strategy for microporous polymer membranes from which we identified those whose FVEs serve as solid solvation cages for lithium ions (Li + ). Furthermore, key elements of our strategy included diversification of bis(catechol)-type monomers via multi-component Mannich reactions to introduce Li+-coordinating functionality within individual FVEs, topology-enforcing polymerizations for generating macromolecular skeletal diversity for networking FVEs into different pore architectures, and several classes of on-polymer reactions for diversifying pore geometries and dielectric properties. Lead candidate polymer membranes featuring explicit ion solvation cages exhibited both higher ionic conductivity and higher cation transference number than control membranes where FVEs were aspecific, which indicates conventional bounds for membrane permeability and selectivity for ion transport can be overcome.4 These advantages are tied to enhanced Li + partitioning from the electrolyte when the cages are present, higher diffusion barriers for anions within the pores, and network-enforced restrictions on the number of solvent molecules bound to Li+ by comparison to the bulk electrolyte, which reduces the effective mass of the working ion. Such membranes show promise as anode-stabilizing interlayers in high-voltage lithium-metal batteries for electric mobility.
The scalability of all-iron redox flow batteries is limited by the plating reaction at the negative electrode on charge. Slurry electrodes have been proposed to allow the plated metal to remain dispersed in the electrolyte and alleviate this limitation. Greater understanding is needed in how faradaic current is distributed between the slurry and the stationary electrode/current collector, especially at relatively low slurry particle loadings. In this work, experiments were conducted using three carbon powders (Kuraray YP-50F, SGL Carbon GFG20HP, and Fisher Scientific GP14736) dispersed in aqueous acidic electrolyte to investigate the effects of slurry particle size, shape, composition, and loading on the behavior of the iron (II/III) reaction. The slurries of YP-50F carbon particles were shown to significantly enhance the faradaic reaction rate due to their high specific surface area but failed to pull the majority of the faradaic current away from the stationary electrode due to the low conductivity of the percolated network. The two graphitic carbons’ slurries did not accelerate the reaction in the same way, but their higher electrical conductivity did succeed in allowing the faradaic reaction to occur primarily on the mobile slurry instead of the stationary current collector.
Thermoset polyurethanes (PUs) have been successfully reprocessed as covalent adaptable networks (CANs) by catalyzing carbamate exchange. Here we extend bond exchange beyond the internal network cross-links to create a dynamic urethane adhesive. Interfacing PU CANs to substrates with nucleophilic functional groups creates adhesives capable of reversible transcarbamoylation with the substrate, which has not been demonstrated previously by CAN adhesives. Two types of thermoset PU films were synthesized, one containing the green carbamate exchange catalyst Zr(tmhd)4 and the other containing no catalyst. Although otherwise identical in chemical and network properties, as indicated by FT-IR spectroscopy and dynamic mechanical thermal analysis (DMTA), the film containing catalyst showed dynamic bond exchange behavior through stress relaxation analysis. When evaluated as an adhesive, the CAN film exhibited self-healing properties and retained its adhesive strength for five cycles, which is attributed to reversible covalent bonding to the glass substrate. Furthermore, this work expands industrially relevant CANs to structural adhesives and demonstrates their potential value in an application that presently employs PUs as single-use materials.
Misfolded tau proteins play a critical role in the progression and pathology of Alzheimer's disease. Recent studies suggest that the spatio-temporal pattern of misfolded tau follows a reaction-diffusion type equation. However, the precise mathematical model and parameters that characterize the progression of misfolded protein across the brain remain incompletely understood. Here, we use deep learning and artificial intelligence to discover a mathematical model for the progression of Alzheimer's disease using longitudinal tau positron emission tomography from the Alzheimer's Disease Neuroimaging Initiative database. Specifically, we integrate physics informed neural networks (PINNs) and symbolic regression to discover a reaction-diffusion type partial differential equation for tau protein misfolding and spreading. First, we demonstrate the potential of our model and parameter discovery on synthetic data. Then, we apply our method to discover the best model and parameters to explain tau imaging data from 46 individuals who are likely to develop Alzheimer's disease and 30 healthy controls. Our symbolic regression discovers different misfolding models f(c) for two groups, with a faster misfolding for the Alzheimer's group, f(c) = 0.23c 3 – 1.34c 2 + 1.11c, than for the healthy control group, f(c) = –c 3 + 0.62c 2 + 0.39c. Our results suggest that PINNs, supplemented by symbolic regression, can discover a reaction-diffusion type model to explain misfolded tau protein concentrations in Alzheimer's disease. Furthermore, we expect our study to be the starting point for a more holistic analysis to provide image-based technologies for early diagnosis, and ideally early treatment of neurodegeneration in Alzheimer's disease and possibly other misfolding-protein based neurodegenerative disorders.
In heterogeneous catalysis, free energy profiles of reactions govern the mechanisms, rates, and equilibria. Energetics are conventionally computed using the harmonic approximation (HA), which requires determination of critical states a priori. Here, we use neural networks to efficiently sample and directly calculate the free energy surface (FES) of a prototypical heterogeneous catalysis reaction—the dissociation of molecular nitrogen on ruthenium—at density-functional-theory-level accuracy. We find that the vibrational entropy of surface atoms, often neglected in HA for transition metal catalysts, contributes significantly to the reaction barrier. The minimum free energy path for dissociation reveals an “on-top” adsorbed molecular state prior to the transition state. While a previously reported flat-lying molecular metastable state can be identified in the potential energy surface, it is absent in the FES at relevant reaction temperatures. Finally, these findings demonstrate the importance of identifying critical points self-consistently on the FES for reactions that involve considerable entropic effects.
A longstanding mystery surrounding electron bifurcation is the significance of inverted (or “crossed”) reduction potentials of the two-electron bifurcating cofactor. Using a many-electron open-system kinetic model, we show that reversible and efficient electron bifurcation is possible without inverted reduction potentials at the bifurcating site if the absolute value of the difference between first and second reduction potentials of the bifurcating species is sufficiently large (on the scale of the redox-potential span of the high- and low-potential branches). Surprisingly, the case with strong, normally ordered potentials at the bifurcating cofactor can produce electron bifurcation that is just as effective as the case with strongly inverted potentials. Lastly, this finding amplifies the puzzle surrounding the recruitment of inverted potentials in the few well-characterized bifurcating systems of nature and suggests that electron bifurcating cofactors without strongly inverted potentials may yet be discovered.
We show that unsupervised machine learning can be used to learn chemical transformation pathways from observational Scanning Transmission Electron Microscopy (STEM) data. To enable this analysis, we assumed the existence of atoms, a discreteness of atomic classes, and the presence of an explicit relationship between the observed STEM contrast and the presence of atomic units. With only these postulates, we developed a machine learning method leveraging a rotationally invariant variational autoencoder (VAE) that can identify the existing molecular fragments observed within a material. The approach encodes the information contained in STEM image sequences using a small number of latent variables, allowing the exploration of chemical transformation pathways by tracing the evolution of atoms in the latent space of the system. The results suggest that atomically resolved STEM data can be used to derive fundamental physical and chemical mechanisms involved, by providing encodings of the observed structures that act as bottom-up equivalents of structural order parameters. The approach also demonstrates the potential of variational (i.e., Bayesian) methods in the physical sciences and will stimulate the development of more sophisticated ways to encode physical constraints in the encoder–decoder architectures and generative physical laws and causal relationships in the latent space of VAEs.
For more than half a century, networks of radiative association, dissociative recombination, and bimolecular reactions have been postulated to drive the low-temperature chemistry of cold molecular clouds. Third-body stabilizations of collision complexes have been assumed to be 'irrelevant' due to short lifetime of such complexes. Here, we conduct crossed molecular beam studies of ground state atomic silicon with diacetylene in combination with electronic structure calculations and microcanonical kinetics models operating under cold molecular cloud conditions. Our combined experimental, electronic structure, and microcanonical kinetics modelling investigations provide compelling evidence that three-body collisions of molecular hydrogen with long-lived reaction intermediates accessed through intersystem crossing are prevalent deep inside molecular clouds. This concept might be exportable to reactions involving polycyclic aromatic hydrocarbons thus affording a versatile machinery to complex organics via third-body stabilizations of bimolecular collision complexes deep inside cold molecular clouds.
TRIMER, Transcription Regulation Integrated with MEtabolic Regulation, is a genome-scale modeling pipeline targeting at metabolic engineering applications. Using TRIMER, regulated metabolic reactions can be effectively predicted by integrative modeling of metabolic reactions with a Transcription Factor (TF)-gene regulatory network (TRN), which is modeled via a Bayesian network (BN). In this paper, we focus on sensitivity analysis of metabolic flux prediction for uncertainty quantification of BN structures for TRN modeling in TRIMER. We propose a computational strategy to construct the uncertainty class of TRN models based on the inferred regulatory order uncertainty given transcriptomic expression data. With that, we analyze the prediction sensitivity of the TRIMER pipeline for the metabolite yields of interest. The obtained sensitivity analyses can guide Optimal Experimental Design (OED) to help acquire new data that can enhance TRN modeling and achieve specific metabolic engineering objectives, including metabolite yield alterations. Here we have performed small- and large-scale simulated experiments, demonstrating the effectiveness of our developed sensitivity analysis strategy for BN structure learning to quantify the edge importance in terms of metabolic flux prediction uncertainty reduction and its potential to effectively guide OED.
Efforts to include more detailed representations of biogeochemical processes in basin-scale water quality simulation tools face the challenge of how to tractably represent mass exchange between the flowing channels of streams and rivers and biogeochemical hotspots in the hyporheic zones. Multiscale models that use relatively coarse representations of the channel network with subgrid models for mass exchange and reactions in the hyporheic zone have started to emerge to address that challenge. Two such multiscale models are considered here, one based on a stochastic Lagrangian travel time representation of advective pumping and one on multirate diffusive exchange. The two models are formally equivalent to well-established integrodifferential representations for transport of non-reacting tracers in steady stream flow, which have been very successful in reproducing stream tracer tests. Despite that equivalence, the two models are based on very different model structures and produce significantly different results in reactive transport. In a simple denitrification example, denitrification is two to three times greater for the advection-based model because the multirate diffusive model has direct connections between the stream channel and transient storage zones and an assumption of mixing in the transient storage zones that prevent oxygen levels from dropping to the point where denitrification can progress uninhibited. By contrast, the advection-based model produces distinct redox zonation, allowing for denitrification to proceed uninhibited on part of the hyporheic flowpaths. These results demonstrate that conservative tracer tests alone are inadequate for constraining representation of mass transfer in models for reactive transport in streams and rivers.
In this work, we present studies of track parameter reconstruction from raw information in CLAS12 detector's Drift Chambers, using Machine Learning (ML). We study the resolution of tracks reconstructed with different types of ML models/algorithms, including Multi-Layer Perceptron (MLP), Extremely Randomized Trees (ERT) and Gradient Boosting Trees (GBT) using simulated data. We find that the resulting ML model is capable of reconstructing track parameters (particle momentum, and polar and azimuthal angles) with accuracy similar to Hit Based (HB) tracking code, but $150$ times faster. Moreover, physics reactions can be identified using the particles reconstructed by the neural network in real-time (with a rate of about $34~kHz$) during experimental data collection. The developed model can be used in numerous applications, such as triggering specific physics reactions in real-time, detector performance monitoring, and real-time detector calibration.
Permeability in reactive porous media evolves in complex ways that are hard to predict. Macroscopic empirical equations are often used to estimate permeability evolution but fail to reflect the impact of pore scale reactions on permeability. This work aims to understand the evolution of permeability in systems where mineral dissolution and precipitation occur simultaneously. Pore network models are used for permeability simulation of a Paluxy sandstone using pore and pore-throat size distributions from the analysis of X-ray CT images. Pore and pore-throat radii are increased and decreased to reflect the effects of dissolution and precipitation reactions. Varying spatial distributions of reactions are simulated and the resulting porosity and permeability values compared with commonly used macroscopic equations. It is observed that when dissolution occurs at the inlet and precipitation at the outlet, porosity increases while permeability decreases. Similar results are also observed in the respective opposite scenario where precipitation occurs at the inlet and dissolution at the outlet. When dissolution and precipitation reactions are randomly distributed throughout the network, porosity increases with little change in permeability that is not captured by the empirical equations. When reactions are controlled by pore and pore-throat sizes, such as dissolution occurring in pores and pore-throats of small size and precipitation in pores and pore-throats of larger size, porosity and permeability decrease. When precipitation occurs in pores and pore-throats of smaller size and dissolution in pores and pore-throats of larger size, porosity and permeability increase and is described relatively well by the Verma-Pruess porosity-permeability equation.