Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stochastic Processes”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Fundamental Interactions of Bimetallic Cu x Pd y ( x + y = 4) Clusters Supported on the α-WC(0001) Surface and Their Performance for CO 2 Adsorption and Dissociation

The tungsten carbide α-WC(0001) surface, an active system for the activation of H 2 and important hydrogenation processes involving unsaturated hydrocarbons, can serve as a support of bimetallic clusters to produce materials with unique catalytic properties, opening routes for a wide range of technical applications. In particular, Cu x Pd y clusters are of particular interest because they combine metals with different properties. A stochastic method was applied to obtain the geometry of Cu x Pd y (x + y = 4) bare clusters, evaluating thousands of possibilities to obtain stable structures, yielding one isomer for Cu 4 , Cu 2 Pd 2 , Cu 1 Pd 3 , and Pd 4 and two isomers for Cu 3 Pd 1 . These clusters were supported on C and W terminations of the tungsten carbide (0001) surface, exploring all of the binding possibilities. The adsorption energies on the C and W terminations are in the ranges from −2.51 to −3.02 eV and from −2.26 to −3.30 eV, respectively. The strongest and weakest binding was seen for monometallic Cu 4 and Pd 4 clusters on both C and W terminations, while the Cu-Pd bimetallics have intermediate adsorption energies but lack a clear trend in terms of composition. The location of Cu x Pd y clusters over the (0001) surface induces a decrease in the work function relative to the pristine surface, while the cluster-surface Bader charge transfer and variations in the partial density of states point to changes in the electronic structure of the carbide atoms upon binding of the metallic clusters. The d-band center of the Cu x Pd y deposited on WC(0001) indicates an intermediate reactivity among Cu(111) and Pd(111) surfaces, modulating the reactivity with small numbers of Cu and Pd atoms, i.e., atom economy in catalyst design. The likelihood of existence of the most stable Cu x Pd y (x + y = 4) clusters in the temperature range of 298-400 K is 100%. The composite Cu x Pd y /α-WC(0001) (x + y = 4), is a nontrivial system since 22 isomers are needed to completely describe its structural properties. Among the isomers, seven structures are necessary to represent Cu 3 Pd 1 /α-WC(0001), five for Pd 4 /α-WC(0001), two for Cu 4 /α-WC(0001), and four for Cu 2 Pd 2 /α-WC(0001) and Cu 1 Pd 3 /α-WC(0001). The large number of cluster isomers supported on the tungsten carbide surface opens the door for several applications in the heterogeneous catalysis of the Cu x Pd y /α-WC(0001) composite, with the possibility of modulating the geometric, electronic, and chemical properties according to a desired application. Test studies for the adsorption of CO 2 indicate that the Cu x Pd y /α-WC(0001) composites are highly active for the adsorption and decomposition of the molecule, with bimetallic and admetal-carbide interactions playing a key role in the binding performance. In conclusion, this high activity indicates that these systems should be useful as catalysts for the conversion of CO 2 to oxygenates or light alkanes.

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

Experimental demonstration of accelerating a beam with a large transverse emittance ratio in the relativistic heavy ion collider for the electron-ion collider

The electron-ion collider (EIC), to be constructed at Brookhaven National Laboratory, will collide polarized high-energy electron beams with hadron beams, achieving luminosities of up to 1.0 × 10 34 cm −2 s −1 in the center-of-mass energy range of 20–140 GeV. To reach such high luminosity, the EIC will employ small, flat beams at the interaction point. According to the design of the EIC hadron storage ring (HSR), hadron beams with a large transverse emittance ratio of 11:1 will be generated at the injection energy using an electron cooling technique and then accelerated to high energies for collisions. Accelerating hadron beams with such a large emittance ratio had never been demonstrated elsewhere—until our recent beam experiment at the relativistic heavy ion collider (RHIC). In this experiment, we successfully generated a large transverse emittance ratio of 13:1 with a gold-ion beam at 31 GeV/nucleon using stochastic cooling. We then accelerated this beam, with a transverse emittance ratio of 11:1, from 31 to 100 GeV/nucleon. Thanks to RHIC’s high-performance orbit, tune, and decoupling feedback systems, the large emittance ratio was well maintained throughout the 5-min-long acceleration process. This experiment fully validated the EIC/HSR design assumptions—namely, that large-emittance-ratio hadron beams can be generated at injection energy and then accelerated to high energies for collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Combining High-Throughput Experiments and Active Learning to Characterize Deep Eutectic Solvents

The high tunability of deep eutectic solvents (DESs) stems from the ease of changing their precursors and relative compositions. However, measuring the physicochemical properties across large composition and temperature ranges, necessary to properly design target-specific DESs, is tedious and error-prone and represents a bottleneck in the advancement and scalability of DES-based applications. As such, active learning (AL) methodologies based on Gaussian processes (GPs) were developed in this work to minimize the experimental effort necessary to characterize DESs. Owing to its importance for large-scale applications, the reduction of DES viscosity through the addition of a low-molecular-weight solvent was explored as a case study. A high-throughput experimental screening was initially performed on nine different ternary DESs. Then, GPs were successfully trained to predict DES viscosity from its composition and temperature, showcasing the ability of these stochastic, nonparametric models to accurately describe the physicochemical properties of complex mixtures. Finally, the ability of GPs to provide estimates of their own uncertainty was leveraged through an AL framework to minimize the number of data points necessary to obtain accurate viscosity modes. This led to a significant reduction in data requirements, with many systems requiring only five independent viscosity data points to be properly described.

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

Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas

Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.

Plasma kinetic theory↗
DOE OSTI2025

The role of the droplet interface in controlling the multiphase oxidation of thiosulfate by ozone

Predicting reaction kinetics in aqueous microdroplets, including aerosols and cloud droplets, is challenging due to the probability that the underlying reaction mechanism can occur both at the surface and in the interior of the droplet. Additionally, few studies directly measure the surface activities of doubly charged anions, despite their prevalence in the atmosphere. Here, deep-UV second harmonic generation spectroscopy is used to probe surface affinities of the doubly charged anions thiosulfate, sulfate, and sulfite, key species in the thiosulfate ozonation reaction mechanism. Thiosulfate has an appreciable surface affinity with a measured Gibbs free energy of adsorption of -7.3 ± 2.5 kJ mol -1 in neutral solution, while sulfate and sulfite exhibit negligible surface propensity. The Gibbs free energy is combined with data from liquid flat jet ambient pressure X-ray photoelectron spectroscopy to constrain the concentration of thiosulfate at the surface in our model. Stochastic kinetic simulations leveraging these novel measurements show that the primary reaction between thiosulfate and ozone occurs at the interface and in the bulk, with the contribution of the interface decreasing from ~65% at pH 5 to ~45% at pH 13. Additionally, sulfate, the major product of thiosulfate ozonation and an important species in atmospheric processes, can be produced by two different pathways at pH 5, one with a contribution from the interface of >70% and the other occurring predominantly in the bulk (>98%). The observations in this work have implications for mining wastewater remediation, atmospheric chemistry, and understanding other complex reaction mechanisms in multiphase environments. Future interfacial or microdroplet/aerosol chemistry studies should carefully consider the role of both surface and bulk chemistry.

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

The Importance of an Ensemble Approach for Modeling Aerosol‐Convection Interactions

Aerosol-convection interactions modulate cloud microphysics, thermodynamics, and updraft intensity, contributing to climate-scale aerosol-radiative forcing. However, quantifying aerosol indirect effects in mixed-phase deep convection remains challenging due to uncertainties in parameterized physics and initial conditions driving nonlinear evolution of convective processes. This study investigates the convective updraft sensitivity to random initial temperature perturbations using an idealized ensemble modeling framework informed by in situ thermodynamic and aerosol observations from the DOE TRACER field campaign. We analyze the impact of small-scale initial perturbations on updraft velocity and supersaturation, and determine the ensemble size required to minimize stochastic internal variability. Results show that minor thermodynamic perturbations can produce updraft variability comparable to aerosol-induced changes reported in prior work. An ensemble of 10 members sufficiently reduces variability, enabling robust investigation of aerosol-related updraft invigoration and informing the design of future ensemble-based aerosol-convection interaction studies in environments with significant spatiotemporal mesoscale thermodynamic and aerosol heterogeneity.

aerosol-convection interactions↗
DOE OSTI2025

AI-powered topic modeling: comparing LDA and BERTopic in analyzing opioid-related cardiovascular risks in women

Topic modeling is a crucial technique in natural language processing (NLP), enabling the extraction of latent themes from large text corpora. Traditional topic modeling, such as Latent Dirichlet Allocation (LDA), faces limitations in capturing the semantic relationships in the text document although it has been widely applied in text mining. BERTopic, created in 2022, leveraged advances in deep learning and can capture the contextual relationships between words. In this work, we integrated Artificial Intelligence (AI) modules to LDA and BERTopic and provided a comprehensive comparison on the analysis of prescription opioid-related cardiovascular risks in women. Opioid use can increase the risk of cardiovascular problems in women such as arrhythmia, hypotension etc. 1,837 abstracts were retrieved and downloaded from PubMed as of April 2024 using three Medical Subject Headings (MeSH) words: “opioid,” “cardiovascular,” and “women.” Machine Learning of Language Toolkit (MALLET) was employed for the implementation of LDA. BioBERT was used for document embedding in BERTopic. Eighteen was selected as the optimal topic number for MALLET and 23 for BERTopic. ChatGPT-4-Turbo was integrated to interpret and compare the results. The short descriptions created by ChatGPT for each topic from LDA and BERTopic were highly correlated, and the performance accuracies of LDA and BERTopic were similar as determined by expert manual reviews of the abstracts grouped by their predominant topics. The results of the t-SNE (t-distributed Stochastic Neighbor Embedding) plots showed that the clusters created from BERTopic were more compact and well-separated, representing improved coherence and distinctiveness between the topics. Our findings indicated that AI algorithms could augment both traditional and contemporary topic modeling techniques. In addition, BERTopic has the connection port for ChatGPT-4-Turbo or other large language models in its algorithm for automatic interpretation, while with LDA interpretation must be manually, and needs special procedures for data pre-processing and stop words exclusion. Therefore, while LDA remains valuable for large-scale text analysis with resource constraints, AI-assisted BERTopic offers significant advantages in providing the enhanced interpretability and the improved semantic coherence for extracting valuable insights from textual data.

Research & Experimental Medicine↗
DOE OSTI2024

Nanometer-Scale Acoustic Wave Packets Generated by Stochastic Core-Level Photoionization Events

We demonstrate that the absorption of femtosecond hard x-ray pulses excites quasispherical, high-amplitude, and high-wave-vector coherent acoustic phonon wave packets using an all hard-x-ray pump-probe scattering experiment. The time- and momentum-resolved diffuse scattering signal is consistent with an ensemble of 3D strain wave packets induced by the rapid electron cascade dynamics following photoionization at uncorrelated excitation centers. We quantify key parameters of this process, including the localization size of the stress field and the photon energy conversion efficiency into elastic energy. The parameters are determined by the photoelectron and Auger electron cascade dynamics, as well as the electron-phonon interaction. In particular, we obtain the localization size of the observed strain wave packet to be 1.5 and 2.5 nm for bulk SrTiO 3 and KTaO 3 single crystals, respectively. The results provide crucial information on the mechanism of x-ray energy deposition into matter and shed light on the shortest collective length scales accessible to coherent acoustic phonon generation using x-ray excitation. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Comprehensive Analysis of Uncertainties in Warm-Rain Parameterizations in Climate Models Based on In Situ Measurements

Abstract Because of the coarse grid size of Earth system models (ESMs), representing warm-rain processes in ESMs is a challenging task involving multiple sources of uncertainty. Previous studies evaluated warm-rain parameterizations mainly according to their performance in emulating collision–coalescence rates for local droplet populations over a short period of a few seconds. The representativeness of these local process rates comes into question when applied in ESMs for grid sizes on the order of 100 km and time steps on the order of 20–30 min. We evaluate several widely used warm-rain parameterizations in ESM application scenarios. In the comparison of local and instantaneous autoconversion rates, the two parameterization schemes based on numerical fitting to stochastic collection equation (SCE) results perform best. However, because of Jessen’s inequality, their performance deteriorates when grid-mean, instead of locally resolved, cloud properties are used in their simulations. In contrast, the effect of Jessen’s inequality partly cancels the overestimation problem of two semianalytical schemes, leading to an improvement in the ESM-like comparison. In the assessment of uncertainty due to the large time step of ESMs, it is found that the rainwater tendency simulated by the SCE is roughly linear for time steps smaller than 10 min, but the nonlinearity effect becomes significant for larger time steps, leading to errors up to a factor of 4 for a time step of 20 min. After considering all uncertainties, the grid-mean and time-averaged rainwater tendency based on the parameterization schemes is mostly within a factor of 4 of the local benchmark results simulated by SCE.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

A meshless stochastic method for Poisson–Nernst–Planck equations

A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.

Chemistry↗
DOE OSTI2024

From Coherence to Function: Exploring the Connection in Chemical Systems

The role of quantum mechanical coherences or coherent superposition states in excited state processes has received considerable attention in the last two decades largely due to advancements in ultrafast laser spectroscopy. These coherence effects hold promise for enhancing the efficiency and robustness of functionally relevant processes, even when confronted with energy disorder and environmental fluctuations. Understanding coherence deeply drives us to unravel mechanisms and dynamics controlled by order and synchronization at a quantum mechanical level, envisioning optical control of coherence to enhance functions or create new ones in molecular and material systems. In this frontier, the interplay between electronic and vibrational dynamics, specifically the influence of vibrations in directing electronic dynamics, has emerged as the leading principle. Here, two energetically disparate quantum degrees of freedom work in-sync to dictate the trajectory of an excited state reaction. Moreover, with the vibrational degree being directly related to the structural composition of molecular or material systems, new molecular designs could be inspired by tailoring certain structural elements. In the realm of chemical kinetics, our understanding of the dynamics of chemical transformations is underpinned by fundamental theories, such as transition state theory, activated rate theory, and Marcus theory. These theories elucidate reaction rates by considering the energy barriers that must be overcome for reactants to transform into products. Those barriers are surmounted by the stochastic nature of energy gap fluctuations within reacting systems, emphasizing that the reaction coordinate, the pathway from reactants to products, is not rigidly defined by a specific vibrational motion but encompasses a diverse array of molecular motions. While less is known about the involvement of specific intramolecular vibrational modes, their significance in certain cases cannot be overlooked. In this Account, we summarize key experimental findings that offer deeper insights into the complex electronic–vibrational trajectories encompassing excited states afforded from state-of-the-art ultrafast laser spectroscopy in three exemplary processes: photoinduced electron transfer, singlet–triplet intersystem crossing, and intramolecular vibrational energy flow in molecular systems. We delve into the rapid decoherence, or loss of phase and amplitude correlations, of vibrational coherences along promoter vibrations during subpicosecond intersystem crossing dynamics in a series of binuclear platinum complexes. This rapid decoherence illustrates the vibration-driven reactive pathways from the Franck–Condon state to the curve crossing region. We also explore the generation of new vibrational coherences induced by impulsive reaction dynamics rather than by the laser pulse in these systems, which sheds light on specific energy dissipation pathways and thereby on the progression of the reaction trajectory in the vicinity of the curve crossing on the product side. Another property of vibrational coherences, amplitude, reveals how energy can flow from one vibration to another in the electronic excited state of a terpyridine–molybdenum complex hosting a nonreactive dinitrogen substrate. In conclusion, a slight change in vibrational energy triggers a quasi-resonant interaction, leading to constructive wavepacket interference and ultimately intramolecular vibrational redistribution from a Franck–Condon active terpyridine vibration to a dinitrogen stretching vibration, energizing the dinitrogen bond.

Electrical energy↗
DOE OSTI2024

Task-oriented machine learning surrogates for tipping points of agent-based models

We present a machine learning framework bridging manifold learning, neural networks, Gaussian processes, and Equation-Free multiscale approach, for the construction of different types of effective reduced order models from detailed agent-based simulators and the systematic multiscale numerical analysis of their emergent dynamics. The specific tasks of interest here include the detection of tipping points, and the uncertainty quantification of rare events near them. Our illustrative examples are an event-driven, stochastic financial market model describing the mimetic behavior of traders, and a compartmental stochastic epidemic model on an Erdös-Rényi network. We contrast the pros and cons of the different types of surrogate models and the effort involved in learning them. Importantly, the proposed framework reveals that, around the tipping points, the emergent dynamics of both benchmark examples can be effectively described by a one-dimensional stochastic differential equation, thus revealing the intrinsic dimensionality of the normal form of the specific type of the tipping point. This allows a significant reduction in the computational cost of the tasks of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Exploring single-cell biosynthetic noise and dynamics for enhanced betaxanthin production in Escherichia coli

Cell-to-cell variability often limits the efficiency of microbial bioproduction, yet how individual cells fluctuate over time and how these fluctuations shape population-level output remain unclear. To address this issue, we tracked a heterologous betaxanthin pathway in Escherichia coli using microfluidics-assisted time-lapse microscopy, allowing simultaneous measurement of fluctuations in betaxanthin, its biosynthetic enzyme DOD and growth across generations. Here we show that over 50% of high betaxanthin producers become medium or low producers after two divisions. Betaxanthin variation primarily originates from DOD noise, with a smaller contribution from growth rate fluctuations. We further develop a stochastic model to explore various control circuits and find that pathway enzyme or metabolite-based growth selection strategies are most effective in enhancing production. We experimentally validate the model by coupling enzyme expression to nutrient availability, which enriches high producers and boosts titer by 4.4-fold. Our results highlight key sources of metabolic heterogeneity and provide a framework for designing robust microbial processes.

Bacterial systems biology↗
DOE OSTI2025

Elucidation of Local Ordering and Atomic-Scale Structure in Polymer-Derived SiOC

Silicon oxycarbide (SiOC) is a versatile ceramic material with tunable microstructure and compositions that can be modulated through precursor chemistry and processing conditions. Though there are several noteworthy uses of SiOC across a range of application spaces, the difficulties in elucidating the short- to medium-range order within these materials have limited the maturation of strategies to precisely control SiC x O 4–x compositions for user-tailored applications. In this contribution, we implement a range of synchrotron scattering and spectroscopy methods coupled with stochastic modeling techniques to elucidate changes in local chemistry and structure associated with the pyrolysis of a commercially available SiOC polymer precursor. Stochastic modeling approaches provide valuable insights into decoupling local Si–O and Si–C environments while confirming predominate heterogeneous phases in materials. Using pyrolysis temperatures between 250 to 800 °C results in a heterogeneous material predominately composed of SiOC and amorphous SiO 2 domains. At 1100 °C, redistribution of Si–C pairs in the SiOC network and Si–O from the SiO 2 domains create a more ordered SiOC phase with local cubic SiC-like ordering. In addition, residual carbon leads to a detectable carbon phases at 1100 °C that persist at higher temperatures. These efforts address the difficulties of obtaining atomic-scale insights into the local structure and nanoscale heterogeneities in SiOC, providing pathways toward establishing structure–property relationships for future materials development.

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

Elastic Stochastic Full Waveform Inversion (eSFWI)

This collaboration between Lawrence Livermore National Security, LLC (LLNS) as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Chevron USA Inc., acting through its Chevron Technical Center division, aimed at developing next-generation computational methods for the Elastic Stochastic Full Waveform Inversion (eSFWI). Seismic imaging is heavily used in the oil and gas industry for identifying and operating subsurface reservoirs. Improved seismic imaging methods can improve productivity, lower costs, and improve operational and environmental safety. This CRADA demonstrated that new high-performance computing (HPC) architectures being rolled out over the next five years can enable unprecedented seismic imaging resolution when using eSFWI techniques to process active seismic data. An open-source computational mini-application was developed, capable of demonstrating near-peak performance for eSFWI algorithms on CPU and GPU enabled HPC platforms. Performance was demonstrated on LLNL HPC systems such as Lassen, as well as on Chevron systems. This project benefited Chevron USA Inc. by demonstrating the potential computational efficiency of their full waveform inversion capabilities used to characterize oil/gas reservoirs, which in turn benefits the public through potential increases in capabilities to perform analysis of leasing sites.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2025

Introduction: Neuromorphic Materials

The explosive growth in data collection and the need to process it efficiently, as well as the desire to automate increasingly complex tasks in transportation, medical care, manufacturing, security and many other fields have motivated a growing interest in neuromorphic computing. Unlike the binary, transistorbased ON/OFF logic gates and separate logic and memory functionalities employed in digital computing, neuromorphic computing is inspired by animal brains that use interconnected synapses and neurons to perform processing, storage and transmission of information at the same location, while only consuming ~20 W or less of power. Motivated by the brain’s efficiency, adaptability, self-learning and resiliency qualities, neuromorphic computing can be broadly defined as an approach to processing and storing information using hardware and algorithms inspired by models of biological neural systems. Present research in neuromorphic computing encompasses approaches that vary significantly in their degree of neuro-inspiration, from systems that only incorporate features such as asynchronous, event-driven operation or use crossbar arrays of non-volatile memory (NVM) elements to accelerate deep neural networks (DNNs), to designs that embrace the extreme parallelism, sparsity, reconfigurability, adaptability, complexity and stochasticity observed in nervous systems. The term ‘neuromorphic’ computing is often credited to Carver Mead, who in the 1980s investigated Si-based analog electronics to replicate functions of the animal retina. Earlier important advances in this field include the work of Frank Rosenblatt, who proposed the concept of the perceptron, Bernard Widrow, who used this concept to build one of the first analog neural networks, the Adaline and many other researchers (see ref. 6 for an historical perspective on neuromorphic computing). With the recent increase in the use of artificial intelligence and large language models, and rising concerns over the associated energy costs, interest in neuromorphic hardware has expanded rapidly. According to some estimates, driven largely by the drastic growth in the training use of artificial intelligence (AI) models using the current computing architectures, the energy cost of computing is projected to reach the energy supply worldwide by 2045. Furthermore, while this is not a realistic outcome, it means that, if more efficient computing technologies are not developed -- soon -- the world will soon become one where demand for energy and market constraints limit the continued increase of societal access to AI and cloud services from data centers. Data centers used for training and use of these models consume hundreds of terawatt hours of electricity, already past 4% of the US electricity demand.

Circuits↗
DOE OSTI2025

Ensemble cure kinetics network (ECK-Net): A method to derive cure kinetics of thermosetting resin

This paper introduces an Ensemble Cure Kinetics Network (ECK-Net), a neural network (NN)–based framework for modeling the cure kinetics of thermosetting resins within a phenomenological context. ECK-Net replaces traditional analytic models, which require extensive chemical insight and multiple isothermal/non-isothermal experiments, with a data-driven surrogate that maps nonlinear relationships between temperature, degree of cure, and reaction rate from differential scanning calorimetry data. The proposed approach predicts input-dependent kinetic coefficients of a generalized nth-order reaction equation rather than reaction rates directly, enabling a single unified model to represent various epoxy systems without relying on iso-conversional analysis or predefined functional forms. To ensure robustness, multiple independently trained networks under different random initializations are blended through an ensemble strategy, effectively mitigating the stochastic variability inherent to neural networks. The framework is validated using experimental datasets from multiple resin systems, including aerospace-grade materials (Toray 3900-2, Cycom 5320-1, and Hexcel 8552) and a windmill-grade resin (RIMR 035c). The model accurately reproduces the temporal evolution of the degree of cure under manufacturers’ recommended cure cycles across all tested resins systems, yielding Pearson’s correlation coefficients of 0.992, 0.994, 0.993, 0.997, respectively. To demonstrate process-level applicability, the trained network was implemented within the Abaqus environment to simulate out-of-autoclave (OOA) curing process of the CFRP panel composed of Toray T830H-6K/3900-2D prepreg. The simulation results showed excellent agreement with experimental temperature response (maximum peak temperature, simulation: 189.6 °C, experiment: 188.5 °C) and the final degree of cure (simulation: 0.948, experiment: 0.960 ± 0.013), confirming ECK-Net’s capability as a reliable alternative to conventional cure kinetics modeling methods.

Composite curing↗
DOE OSTI2026

Spatial Correlations of the Poisson Model for Radiation Transport

Characterizing the relationship between bulk physical properties and mixing in randomly heterogeneous media is a central challenge across many areas of science and engineering. A benchmark model for such studies is the Poisson model, a random tessellation of space by a Poisson process of hyperplanes. In radiation transport studies, the lack of exact expressions for the Poisson model’s spatial multipoint functions has led to approximate methods being used, introducing unquantified sources of error. Here, we recently introduced an exact solution for the Poisson model’s multipoint functions and closely related conditional probability functions (CPFs), providing a new opportunity to understand and reduce these sources of error. In this paper, we enable a more rigorous investigation of radiation transport in stochastic media by applying the recently introduced exact solution for the Poisson model’s CPFs. This paper consists of three main contributions. First, we introduce a unified framework for CPFs of the Poisson model, encompassing the recently introduced exact CPFs as well as the previously introduced atomic mix, nearest-neighbor, and combination CPFs. This framework also includes existing pruning techniques for the approximate CPFs, such as angular exclusion, as well as a novel form of angular exclusion suitable for the exact CPFs. Second, we use the exact CPFs to characterize the spatial regions where each approximate three-point CPF is most accurate, thereby explaining the observed hierarchy of accuracy among the approximate models. Finally, we evaluate material transmittance, reflectance, and flux in a three-dimensional test problem using conditional point sampling, demonstrating the relationship between CPF accuracy and transport simulation accuracy.

Poisson model↗