Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complex systems”

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

Simultaneous Degradation‐Depolymerization of Bioderived Comb Copolymers

Poly(lactic acid) (PLA) is the most widely explored biodegradable alternative for polystyrene; however, its low toughness and glass transition temperature may limit its wider adoption as a sustainable replacement. To improve its material and thermal properties, PLA can be chemically or physically combined with other polymers, like poly(methyl methacrylate) (PMMA), though the incorporation of vinyl‐based polymer components complicates chemical recycling and reduces the sustainability of the material. Here, in this study, we synthesized polymethacrylate‐PLA comb copolymers designed to be thermally deconstructed. Our design strategically extends current polymer deconstruction methodologies to more complex macromolecular systems. Thermally labile units within the polymethacrylate backbone permitted depolymerization that was concurrent with PLA side chain degradation during heating. This dual degradation‐depolymerization process enhances the overall sustainability of lactide/vinyl‐based copolymers and demonstrates the synergistic potential of integrating multiple deconstruction pathways into a single system. This report elaborates on the design of advanced, degradable copolymers, contributing to the further development of sustainable polymer materials.

PLA↗
DOE OSTI2024

Improving the reliability of machine learned potentials for modeling inhomogeneous liquids

The atomic-scale response of inhomogeneous fluids at interfaces and surrounding solute particles plays a critical role in governing chemical, electrochemical, and biological processes. Classical molecular dynamics simulations have been applied extensively to simulate the response of fluids to inhomogeneities directly, but are limited by the accuracy of the underlying interatomic potentials. Here, we use neural network potentials (NNPs) trained to ab initio simulations to accurately predict the inhomogeneous responses of two distinct fluids: liquid water and molten NaCl. Although NNPs can be readily trained to model complex bulk systems across a range of state points, we show that to appropriately model a fluid's response at an interface, relevant inhomogeneous configurations must be included in the training data. In order to sufficiently sample appropriate configurations of such inhomogeneous fluids, we develop protocols based on molecular dynamics simulations in the presence of external potentials. We demonstrate that NNPs trained on inhomogeneous fluid configurations can more accurately predict several key properties of fluids—including the density response, surface tension and size-dependent cavitation free energies—for liquid water and molten NaCl, compared to both empirical interatomic potentials and NNPs that are not trained on such inhomogeneous configurations. This work therefore provides a first demonstration and framework to extract the response of inhomogeneous fluids from first principles for classical density-functional treatment of fluids free from empirical potentials.

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

Characterizing climate pathways using feature importance on echo state networks

The 2022 National Defense Strategy of the United States listed climate change as a serious threat to national security. Climate intervention methods, such as stratospheric aerosol injection, have been proposed as mitigation strategies, but the downstream effects of such actions on a complex climate system are not well understood. The development of algorithmic techniques for quantifying relationships between source and impact variables related to a climate event (i.e., a climate pathway) would help inform policy decisions. Data-driven deep learning models have become powerful tools for modeling highly nonlinear relationships and may provide a route to characterize climate variable relationships. In this paper, we explore the use of an echo state network (ESN) for characterizing climate pathways. ESNs are a computationally efficient neural network variation designed for temporal data, and recent work proposes ESNs as a useful tool for forecasting spatiotemporal climate data. However, ESNs are noninterpretable black-box models along with other neural networks. The lack of model transparency poses a hurdle for understanding variable relationships. We address this issue by developing feature importance methods for ESNs in the context of spatiotemporal data to quantify variable relationships captured by the model. We conduct a simulation study to assess and compare the feature importance techniques, and we demonstrate the approach on reanalysis climate data. In the climate application, we consider a time period that includes the 1991 volcanic eruption of Mount Pinatubo. This event was a significant stratospheric aerosol injection, which acts as a proxy for an anthropogenic stratospheric aerosol injection. Furthermore, we are able to use the proposed approach to characterize relationships between pathway variables associated with this event that agree with relationships previously identified by climate scientists.

black-box models↗
DOE OSTI2025

Robust Spectral Anomaly Detection in EELS Spectral Images via 3D Convolutional Variational Autoencoders

Abstract A 3D Convolutional Variational Autoencoder (3D‐CVAE) is introduced for automated anomaly detection in electron energy‐loss spectroscopy spectrum imaging (EELS‐SI) data. This approach leverages the full 3D structure of EELS‐SI data to detect subtle spectral anomalies while preserving both spatial and spectral correlations across the datacube. By employing cross‐entropy loss and training on bulk spectra, the model learns to reconstruct bulk features characteristic of the defect‐free material. In exploring methods for anomaly detection, both the 3D‐CVAE approach and principal component analysis (PCA) are evaluated, testing their performance using FeL‐edge ΔEpeak shifts designed to simulate material defects. These results show that 3D‐CVAE achieves superior anomaly detection and maintains consistent performance across various shift magnitudes. The method demonstrates clear bimodal separation between bulk and anomalous spectra, enabling reliable classification. Further analysis verifies that lower‐dimensional representations are robust to anomalies in the data. While performance advantages over PCA diminish with decreasing anomaly concentration, our method maintains high reconstruction quality even in challenging, noise‐dominated spectral regions. This approach provides a robust framework for unsupervised automated detection of spectral anomalies in EELS‐SI data, particularly valuable for analyzing complex material systems.

Chemistry↗
DOE OSTI2026

Dynamics of McMillan mappings III. Symmetric map with mixed nonlinearity

This article extends the study of the dynamical properties of the symmetric McMillan map, emphasizing its utility in understanding and modeling complex nonlinear systems. Although the map features six parameters, we demonstrate that only two are irreducible: the linearized rotation number at the fixed point and a nonlinear parameter representing the ratio of terms in the biquadratic invariant. Through a detailed analysis, we classify regimes of stable motion, provide exact solutions to the mapping equations, and derive a canonical set of action-angle variables, offering analytical expressions for the rotation number and nonlinear tune shift. We further establish connections between general standard-form mappings and the symmetric McMillan map, using the area-preserving Hénon map and accelerator lattices with thin sextupole magnet as representative case studies. Our results show that, despite being a second-order approximation, the symmetric McMillan map provides a highly accurate depiction of dynamics across a wide range of system parameters, demonstrating its practical relevance in both theoretical and applied contexts.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Imaging the initial condition of heavy-ion collisions and nuclear structure across the nuclide chart

High-energy nuclear collisions encompass three key stages: the structure of the colliding nuclei, informed by low-energy nuclear physics, the initial condition , leading to the formation of quark–gluon plasma (QGP), and the hydrodynamic expansion and hadronization of the QGP, leading to final-state hadron distributions that are observed experimentally. Recent advances in both experimental and theoretical methods have ushered in a precision era of heavy-ion collisions, enabling an increasingly accurate understanding of these stages. However, most approaches involve simultaneously determining both QGP properties and initial conditions from a single collision system, creating complexity due to the coupled contributions of these stages to the final-state observables. To avoid this, we propose leveraging established knowledge of low-energy nuclear structures and hydrodynamic observables to independently constrain the QGP’s initial condition. By conducting comparative studies of collisions involving isobar-like nuclei—species with similar mass numbers but different ground-state geometries—we can disentangle the initial condition’s impacts from the QGP properties. This approach not only refines our understanding of the initial stages of the collisions but also turns high-energy nuclear experiments into a precision tool for imaging nuclear structures, offering insights that complement traditional low-energy approaches. Opportunities for carrying out such comparative experiments at the Large Hadron Collider and other facilities could significantly advance both high-energy and low-energy nuclear physics. Additionally, this approach has implications for the future electron-ion collider. While the possibilities are extensive, we focus on selected proposals that could benefit both the high-energy and low-energy nuclear physics communities. Originally prepared as input for the long-range plan of U.S. nuclear physics, this white paper reflects the status as of September 2022, with a brief update on developments since then.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Accelerating Hamiltonian Monte Carlo for Bayesian inference in neural networks and neural operators

Hamiltonian Monte Carlo (HMC) is a powerful and accurate method to sample from the posterior distribution in Bayesian inference. However, HMC techniques are computationally demanding for Bayesian neural networks due to the high dimensionality of the network’s parameter space and the non-convexity of their posterior distributions. Therefore, various approximation techniques, such as variational inference (VI) or stochastic gradient MCMC, are often employed to infer the posterior distribution of the network parameters. Such approximations introduce inaccuracies in the inferred distributions, resulting in unreliable uncertainty estimates. In this work, we propose a hybrid approach that combines inexpensive VI and accurate HMC methods to efficiently and accurately quantify uncertainties in neural networks and neural operators. The proposed approach leverages an initial VI training on the full network. We examine the influence of individual parameters on the prediction uncertainty, which shows that a large proportion of the parameters do not contribute substantially to uncertainty in the network predictions. This information is then used to significantly reduce the dimension of the parameter space, and HMC is performed only for the subset of network parameters that strongly influence prediction uncertainties. This yields a framework for accelerating the full batch HMC for posterior inference in neural networks. We demonstrate the efficiency and accuracy of the proposed framework on deep neural networks and operator networks, showing that inference can be performed for large networks with tens to hundreds of thousands of parameters. Finally, we show that this method can effectively learn surrogates for complex physical systems by modeling the operator that maps from upstream conditions to wall-pressure data on a cone in hypersonic flow.

Bayesian inference↗
DOE OSTI2024

Multi-level Monte Carlo methods in chemical applications with Lennard-Jones potentials and other landscapes with isolated singularities

We describe and compare outcomes of various Multi-Level Monte Carlo (MLMC) method variants, motivated by the potential of improved computational efficiency over rejection based Monte Carlo, which scales poorly with problem dimension. With an eye toward its application to computational chemical physics, we test MLMC's ability to sample trajectories on two problems — a familiar double-well potential, with known stationary distributions, and a Lennard-Jones solid potential (a Galton Board). By sampling Brownian motion trajectories, we are able to compute expectations of observable averages. These multi-basin potential energy problems capture the essence of the challenges with using MLMC, namely, maintaining correspondence of sample paths as time-resolution is varied. Addressing this challenge properly can lead to MLMC significantly outperforming standard Monte Carlo path sampling. We describe the essence of this problem and suggest strategies that circumvent diverging multilevel sample paths for an important class of problems. In the tests we also compare the computational cost of several, “adaptive,” variants of MLMC. Our results demonstrate that MLMC overcomes the collision, time scale limitation of the more familiar Brownian path MC samplers, and our implementation provides tunable error thresholds, making MLMC a promising candidate for application to larger and more complex molecular systems.

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

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory↗
DOE OSTI2025

Automatic building energy model development and debugging using large language models agentic workflow

Building energy modeling (BEM) is a complex process that demands significant time and expertise, limiting its broader application in building design and operations. While Large Language Models (LLMs) agentic workflow have facilitated complex engineering processes, their application in BEM has not been specifically explored. This paper investigates the feasibility of automating BEM using LLM agentic workflow. Here, we developed a generic LLM-planning-based workflow that takes a building description as input and generates an error-free EnergyPlus building energy model. Our robust workflow includes four core agents: 1) Building Description Pre-Processing, 2) IDF Object Information Extraction, 3) Single IDF Object Generator Suite, and 4) IDF Debugging Agent. These agents divide the complex tasks into manageable sub-steps, enabling LLMs to generate accurate and reliable results at each stage. The case study demonstrates the successful translation of a building description into an error-free EnergyPlus model for the iUnit modular building at the National Renewable Energy Laboratory. The effectiveness of our workflow surpasses: 1) naive prompt engineering, 2) other LLM-based workflows, and 3) manual modeling, in terms of accuracy, reliability, and time efficiency. The paper concludes with a discussion on the interplay between foundational models and LLM agent planning design, advocating for the use of fine-tuned, specialized models to advance this field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A score-based diffusion model approach for adaptive learning of stochastic partial differential equation solutions

In this paper, we propose a novel framework for adaptively learning the time-evolving solutions of stochastic partial differential equations (SPDEs) using score-based diffusion models within a recursive Bayesian inference setting. SPDEs play a central role in modeling complex physical systems under uncertainty, but their numerical solutions often suffer from model errors and reduced accuracy due to incomplete physical knowledge and environmental variability. To address these challenges, we encode the governing physics into the score function of a diffusion model using simulation data and incorporate observational information via a likelihood-based correction in a reverse-time stochastic differential equation. This enables adaptive learning through iterative refinement of the solution as new data becomes available. To improve computational efficiency in high-dimensional settings, we introduce the ensemble score filter, a training-free approximation of the score function designed for real-time inference. Numerical experiments on benchmark SPDEs demonstrate the accuracy and robustness of the proposed method under sparse and noisy observations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A stress-sensitive precipitate nucleation model beyond classical nucleation theory

The dynamic evolution of precipitates and second phases dictates the strength and stability of most engineering alloys. By design, or as a consequence of thermo-mechanical aging, engineering metals and alloys often form precipitates of second phases when subjecting to diverse thermal and mechanical loads. Precipitation is governed by several factors, including the alloy’s composition, processing/operating temperature, and stresses — either as a result of external loads or from residual stresses. However, state-of-the-art models for precipitate nucleation (i.e., classical nucleation theory) typically lacks consistent method to capture the effects of externally applied and/or internal stresses on nucleation; thereby severely limiting the applicability of these models to complex materials systems and to representative loading scenarios. Here, in this work, we extend upon classical nucleation theory to account for the effect of stresses on precipitation kinetics and thermodynamics. This is achieved via the use of an Eshelbian micromechanics framework keeping track of (i) the stress build up resulting from second phase formation as a function of mechanical load and, (ii) the effects of dislocations on precipitate formation. This new model is applied to σ precipitate in Fe–Cr binary alloys and M 23 C 6 precipitate in 316H stainless steel (SS). Simulations demonstrate the important role of both the remotely applied loads and dislocation pile ups on precipitate nucleation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Capturing the Page curve and entanglement dynamics of black holes in quantum computers

Quantum computers are emerging technologies expected to become important tools for exploring various aspects of fundamental physics in the future. Therefore, we pose the question of whether quantum computers can help us to study the Page curve and the black hole information dynamics, which has been a key focus in fundamental physics. In this regard, we rigorously examine the qubit transport model, a toy qubit model of black hole evaporation on IBM’s superconducting quantum computers, to shed light on this question. Specifically, we implement the quantum simulation of the scrambling dynamics in black holes using an efficient random unitary circuit. Furthermore, we employ the swap-based many-body interference protocol and the randomized measurement protocol to measure the entanglement entropy of Hawking radiation qubits in this model. Finally, by incorporating quantum error mitigation techniques into our challenging implementation of entanglement entropy measurement protocols on the IBM quantum hardware, we accurately determine the Rényi entropy in the qubit transport model, thus showcasing the utility of quantum computers for future investigations of complex quantum systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Coupling Microdroplet-Based Sample Preparation, Multiplexed Isobaric Labeling, and Nanoflow Peptide Fractionation for Deep Proteome Profiling of the Tissue Microenvironment

There is increasing interest in developing in-depth proteomic approaches for mapping tissue heterogeneity in a cell-type-specific manner to better understand and predict the function of complex biological systems such as human organs. Existing spatially resolved proteomics technologies cannot provide deep proteome coverage due to limited sensitivity and poor sample recovery. Herein, we seamlessly combined laser capture microdissection with a low-volume sample processing technology that includes a microfluidic device named microPOTS (microdroplet processing in one pot for trace samples), multiplexed isobaric labeling, and a nanoflow peptide fractionation approach. The integrated workflow allowed us to maximize proteome coverage of laser-isolated tissue samples containing nanogram levels of proteins. We demonstrated that the deep spatial proteomics platform can quantify more than 5000 unique proteins from a small-sized human pancreatic tissue pixel (∼60,000 μm2) and differentiate unique protein abundance patterns in pancreas. Furthermore, the use of the microPOTS chip eliminated the requirement for advanced microfabrication capabilities and specialized nanoliter liquid handling equipment, making it more accessible to proteomic laboratories.

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

Real-Time Detection of Hydrogen and Ammonia Isotopologues for Impurity Removal and Recovery of Tritium

To accommodate gas measurements for impurity removal and recovery of tritium, a silver-coated optical or waveguide is employed for collecting Raman scattered signals to determine relative hydrogen and ammonia isotopologue populations in real time. The data and results presented here demonstrate an analytical methodology for the analysis of four ammonia and three hydrogen isotopologues in a hydrogen–deuterium exchange reaction by gas phase Raman spectroscopy. Standard chemometric modeling techniques effectively unravel the signatures of the isotopologues involved observed here; however, a sophisticated quantum chemical approach supports the spectral assignments. An interpretation of the data presented here can emphasize the practicality and reliability of the gaseous monitoring system in complex chemical environments for the hydrogen fuel economy as well as the more distant energy source from a facility that handles tritium. There are still considerable concerns about the measurement of tritium in isotope separation and radiological impurities from gas processing. A common impurity in gas processing is ammonia, which can form readily in the presence of nitrogen and tritium. Substituted ammonia (NQ 3 ), where Q = H, D, or T, is traditionally removed through getters or diffusers along with other non-hydrogen contaminants. A preferable analytical approach is noninvasive and can be deployed for real-time process evaluation in radiological environments.

Ammonia↗
DOE OSTI2025

A Synchrotron-Based Vacuum Ultraviolet Photoionization Mass Spectrometer-Coupled Microreactor To Probe Thermocatalysis

Vacuum ultraviolet photoionization (VUV-PI) mass spectrometry offers an isomer-selective and universal ionization with minimal fragmentation detection of organics in complex chemical systems such as pyrolysis and combustion. Here, in this study, we report a state-of-the-art experimental setup of a universal catalytic microreactor combined with a molecular beam to investigate the thermocatalytic oxidation of a heterogeneous system relevant for probing reactions at gas–solid interfaces. In strong contrast to traditional off-line analytical methods, this technique is capable of identifying and quantifying short-lived species (radicals) as well as stable products to decipher initial reaction steps via the detection of nascent products. The thermocatalytic oxidative degradation of exo-tetrahydrodicyclopentadiene (JP-10), a high energy-density hydrocarbon fuel, over solid titanium–aluminum–boron reactive mixed metal nanopowder (Ti-Al-B RMNP) is exploited to showcase potential applications. Overall, some 59 nascent gas-phase products are identified via photoionization efficiency (PIE) curves, including oxygenated species and hydrocarbons comprising closed-shell molecules and radicals. The critical temperature for complete oxidative decomposition of JP-10 was lowered by 450 K from 1400 K to 950 K, indicating an efficient thermocatalytic action of Ti-Al-B nanoparticles on JP-10. The enabling of a universal chemical microreactor along with VUV-PI mass spectrometry broadens the applicability of this technique to hydrocarbon fuel oxidation and pyrolysis characterization. This isomer-selective sensitive probing along with the detection of radical transients makes the aforementioned technique superior to other conventional analytical techniques such as microflow tube and pyrolysis-gas chromatography coupled with mass spectrometry for investigating similar pyrolysis reactions and comprehensive quantification.

Paul, Dababrata [Univ. of Hawaii at Manoa, Honolul↗
DOE OSTI2026

Data-Driven Discovery and Experimental Validation of Solvent Polarity Effects on Conjugated Polymer Solution-to-Film Assembly Pathways

Understanding how solvent properties influence the solution-to-film assembly of conjugated polymers remains a critical challenge due to the complex and intertwined nature of polymer–solvent interactions. In this study, we integrate a data-driven framework with experimental validation to identify key parameters influencing the assembly and performance of poly[2,5-(2-octyldodecyl)-3,6-diketopyrrolopyrrole-alt-5,5-(2,5-di(thien-2-yl)thieno[3,2-b]thiophene)] (DPP-DTT) in organic field-effect transistors (OFETs). A machine learning (ML) approach identified the normalized Reichardt polarity parameter (E T N ) as a significant descriptor correlated with DPP-DTT hole mobility (μ). Systematic DPP-DTT devices fabricated using solvents across a wide E T N range revealed that higher E T N solvents yield enhanced μ. To elucidate the structural origins of high μ, we conducted comprehensive analyses using UV–vis–NIR spectroscopy and grazing incidence wide angle X-ray scattering (GIWAXS) measurements. The results revealed that films processed from high E T N solvents exhibit reduced paracrystallinity. By analyzing the solution-state behavior using optical microscopy and solution WAXS, we revealed polymer solubility differences in the various solvents and associated distinct polymer assembly pathways, elucidating why the high E T N solvent produces long-range ordered films. Notably, the high E T N solvent shows a pronounced preference for liquid-crystal (LC)-mediated assembly, providing a mechanistic explanation for the enhanced structural order. Therefore, these results demonstrate that solvent polarity, as evaluated by E T N , serves as an important parameter that plays a significant role in the DPP-DTT assembly pathway and resultant solid-state morphology. This work provides a strategy for integrating data science with experiments to identify critical parameters associated with complex polymer systems and helps guide rational process design for high-performance organic electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Accelerating Instanton Theory with the Line Integral Nudged Elastic Band Method and Gaussian Process Regression

Quantum tunneling plays a fundamental role in many chemical reactions, particularly proton transfer processes. Ring polymer instanton theory offers a practical framework for computing tunneling rates in complex molecular systems. However, applying the ring polymer instanton method with a potential energy surface generated on-the-fly using electronic structure calculations can be computationally demanding. Here, in this work, we present a new efficient implementation of the ring polymer instanton method by combining the Line Integral Nudged Elastic Band (LI-NEB) approach with Gaussian Process Regression (GPR). We benchmarked this method on prototypical ground-state proton transfer systems, including the benchmark gas-phase hydrogen abstraction reaction H + CH 4 → H 2 + CH 3 , malonaldehyde, and Z-3-amino-propenal (aminopropenal). Our results show that this approach is an order of magnitude faster than traditional instanton algorithms while maintaining excellent agreement with their tunneling rates. This development opens the door to studying proton transfer in larger systems with improved efficiency.

chemical physics↗