SEARCH · Engineering Papers
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.
We present a generalizable database architecture ESAMP that captures the complete provenance associated with a material. We demonstrate this architecture and provenance based machine learning on one of the largest experimental materials databases.
Solutions containingn-alkylammonium polyselenides constitute a new reactive solvent system. This reactive dissolution chemistry is broadly generalizable and can enable solution processed synthesis of a wide range of metal chalcogenide semiconductors, both as thin films and nanoparticles.
A machine learning model for reliable director fields calculation from raw experimental images of active nematics. The model is accurate, robust to noise and generalizable, enhancing analysis such as the detection and tracking of topological defects.
To reliably deploy lithium-ion batteries, a fundamental understanding of cycling aging behavior is critical. Battery aging consists of complex and highly coupled phenomena, making it challenging to develop a holistic interpretation. In this work, we generate a diverse battery cycling dataset with a broad range of degradation trajectories, consisting of 359 high energy density commercial Li(Ni,Co,Al)O 2 /graphite + SiO x cylindrical 21 700 cells cycled across 207 unique cycling protocols. We consolidate aging via 16 mechanistic state-of-health (SOH) metrics, including cell-level performance metrics, electrode-specific capacities/state-of-charges (SOCs), and aging trajectory metrics. We develop a framework using interpretable machine learning and explainable features to generate an aging matrix that visually deconvolutes the complex battery degradation behavior. This generalizable data-driven mechanistic framework simplifies the complex interplay between cycling conditions, degradation modes, and SOH, acting as a hypothesis-generation tool to aid battery users in identifying key degradation regimes for further study and experimentation.
High-throughput density functional theory (DFT) calculations have become a vital element of computational materials science, enabling materials screening, property database generation, and training of “universal” machine learning models. While several software frameworks have emerged to support these computational efforts, new developments such as machine learned force fields have increased demands for more flexible and programmable workflow solutions. This manuscript introduces atomate2, a comprehensive evolution of our original atomate framework, designed to address existing limitations in computational materials research infrastructure. Key features include the support for multiple electronic structure packages and interoperability between them, along with generalizable workflows that can be written in an abstract form irrespective of the DFT package or machine learning force field used within them. Our hope is that atomate2's improved usability and extensibility can reduce technical barriers for high-throughput research workflows and facilitate the rapid adoption of emerging methods in computational material science.
Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.
Interfacial instability between lithium metal and solid-state electrolytes limits the performance of all-solid-state lithium metal batteries (ASSLBs), leading to parasitic reactions, non-uniform Li+ flux, and dendrite growth. Here, we develop a composite interlayer composed of the anti-perovskite Li2OHCl0.75Br0.25 (AP) and carbon nanotubes (CNTs) to enhance both interfacial stability and ionic transport. The AP–CNT interlayer exhibits enhanced Li+ conductivity arising from interfacial electron transfer from AP to CNTs, which generates a built-in electric field that facilitates Li+ migration. Lithium symmetric cells incorporating this interlayer achieve a high critical current density of 2.4 mA cm−2 at 55 °C. This design integrates chemical robustness with coupled ion–electron transport, offering a generalizable strategy for safe, dendrite-free, and high-performance ASSLBs.
Accurately predicting solid electrolyte interphase (SEI) formation requires explicitly resolving the electric double layer (EDL) structure, which deviates significantly from that of the bulk electrolyte. Although an established molecular dynamics (MD) and Density Functional Theory (DFT) framework can model SEI formation by evaluating reduction reactions of local clusters in the EDL, it suffers from a combinatorial computational bottleneck. To overcome this limitation, we introduce a machine-learning-accelerated simulation workflow (MD–DFT–ML), integrating a gradient-boosted regression model trained on EDL composition data to efficiently predict reduction potentials. We apply this framework to seven fluorinated electrolytes comprising fluorinated anions, a fluorinated ester solvent, two types of diluent (ion-solvating ester vs. non-solvating ether), and an FEC additive. The analysis shows that the EDL selectively accumulates cation-binding species; consequently, the non–cation-binding ether diluent rarely enters the EDL and makes minimal contributions to SEI formation. DFT calculations on statistically representative EDL clusters provide reduction potentials and fluorine-release pathways, while the ML model, which substantially reduces the DFT workload, predicts cluster reduction energies with a mean absolute error of 0.1 eV. The combined MD–DFT–ML approach also quantifies contributions from different sources to LiF formation in the SEI. This methodology establishes a generalizable route for multiscale modeling electrolyte and interphase design for next-generation electrochemical energy-storage systems.
Poly(ethylene terephthalate) (PET) oligomers derived from glycolysis depolymerization were converted into multiblock copolymers through diisocyanate-mediated coupling with dihydroxy-terminated oligomers, enabling precise control over copolymer sequence distribution, connectivity, and mechanical performance. Here, we demonstrate that telechelic PET oligomers isolated directly from depolymerized consumer waste can serve as reactive building blocks for the formation of segmented multiblock copolymers, eliminating the need to revert to monomeric feedstocks. Dihydroxy-terminated PET oligomers (Mw ≈ 8 kg mol−1) were coupled with poly(ethylene oxide) (PEO, Mw ≈ 4 kg mol−1) to form PET-PEO multiblock copolymers with high molar mass (Mw ≈ 160 kg mol−1). We further show that the timing of end-capping reactions provides a key control parameter that governs the competition between chain extension and termination, thereby dictating the resulting multiblock architecture and molecular-weight evolution. Evaluation of their mechanical properties reveals that virgin PET exhibits high modulus (∼3 GPa) and strength (43 MPa) but limited ductility (<10% elongation). In contrast, PET–PEO multiblock copolymers retain comparable tensile strength (44 MPa), albeit while exhibiting dramatically enhanced ductility (>90% elongation), forming tougher materials with efficient stress transfer between the rigid PET domains and the flexible PEO segments. When incorporated into PET/PEO blends at low loadings, the multiblock copolymers serve as effective compatibilizers, yielding materials with an intermediate modulus (1.1–1.2 GPa) and improved elongation compared to uncompatibilized blends. Furthermore, the presence of PEO blocks increases water uptake and gas permeability relative to virgin PET, reflecting the tunability of molecular transport through the copolymeric blocks. To our knowledge, this represents the first report of PET–PEO multiblock copolymers derived from post-consumer PET for gas transport applications. These results demonstrate that multiblock copolymer formation from telechelic PET oligomers provides a versatile platform for tailoring the mechanical and transport behavior of polyester-based materials through controlled macromolecular design and establishes a generalizable strategy for transforming consumer plastic waste into functional segmented polymers without requiring complete depolymerization to monomers.
Polymers of intrinsic microporosity (PIMs) offer exceptional gas permeability but remain brittle and susceptible to physical aging, limiting their durability in separation applications. Here, we introduce a reconfigurable microporous polymer network that uniquely integrates permanent PIM microporosity with autonomous, intrinsic self-healing driven by imidazolium-based ionic motifs. Spirobisindane units generate the intrinsic free-volume architecture, while an imidazolium-containing polyamide ionene supplies dynamic ionic and hydrogen-bonding interactions that reorganize under mild activation. Incorporation of imidazolium-based ionic liquids further tunes cohesion, mobility, and densification, enabling the network to relax, re-associate, and retain microporosity without structural collapse. Through a comprehensive multiscale approach combining spectroscopy, scattering, thermal and mechanical characterization with all-atom molecular dynamics and density functional theory calculations, we elucidate how ionic content, as a single control parameter that reshapes free-volume distributions, modulates local coordination environments, and governs relaxation and healing kinetics. At intermediate ionic loadings, the networks achieve rapid, repeatable self-healing while maintaining CO$_2$ selectivity, demonstrating an optimal balance between segmental mobility and structural integrity. By establishing how hierarchical ionic interactions couple structure, dynamics, and transport in microporous ionene networks, this work provides generalizable design rules for adaptive soft-matter systems that require simultaneous mechanical resilience, reconfigurability, and selective gas transport.
Ternary Zintl phosphides are promising light-absorbing semiconductors for thin-film optoelectronic applications, but strategies for controlling their microstructure and optoelectronic quality remain underexplored. Here, we report the synthesis of phase-pure SrZn2P2 thin films using radio-frequency co-sputtering in a PH3 + Ar atmosphere and investigate the impact of post-growth processing on their structural and optical properties. Grazing-incidence X-ray scattering and Raman spectroscopy confirm the formation of crystalline SrZn2P2 films over a finite compositional window. Optical measurements reveal clear absorption near the direct-band-gap energy (~1.8 eV) and near-band-edge photoluminescence. Further, we have studied the effects of chemically compatible halide-assisted annealing. It is found that SrI2 treatments lead to pronounced grain growth and reduced diffraction peak broadening while preserving phase purity, in contrast to rapid thermal or forming-gas annealing. Notably, annealing with SrI2 at 450 degrees C significantly enhances both the intensity and spatial uniformity of the photoluminescence, thus connecting the observed microstructural consolidation with improved radiative recombination. Our study demonstrates that halide-assisted annealing provides an effective pathway for microstructural control in SrZn2P2 thin films and highlights a generalizable processing strategy for advancing Zintl phosphide semiconductors toward optoelectronic applications.
We present a Bayesian calibration algorithm for cosmic microwave background (CMB) observations as implemented within the global end-to-end BEYONDPLANCK framework and applied to the Planck Low Frequency Instrument (LFI) data. Following the most recent Planck analysis, we decomposed the full time-dependent gain into a sum of three nearly orthogonal components: one absolute calibration term, common to all detectors, one time-independent term that can vary between detectors, and one time-dependent component that was allowed to vary between one-hour pointing periods. Each term was then sampled conditionally on all other parameters in the global signal model through Gibbs sampling. The absolute calibration is sampled using only the orbital dipole as a reference source, while the two relative gain components were sampled using the full sky signal, including the orbital and Solar CMB dipoles, CMB fluctuations, and foreground contributions. We discuss various aspects of the data that influence gain estimation, including the dipole-polarization quadrupole degeneracy and processing masks. Comparing our solution to previous pipelines, we find good agreement in general, with relative deviations of -0.67% (-0.84%) for 30 GHz, 0.12% (-0.04%) for 44 GHz and -0.03% (-0.64%) for 70 GHz, compared to Planck PR4 and Planck 2018, respectively. We note that the BEYONDPLANCK calibration was performed globally, which results in better inter-frequency consistency than previous estimates. Additionally, WMAP observations were used actively in the BEYONDPLANCK analysis, which both breaks internal degeneracies in the Planck data set and results in an overall better agreement with WMAP. Finally, we used a Wiener filtering approach to smoothing the gain estimates. We show that this method avoids artifacts in the correlated noise maps as a result of oversmoothing the gain solution, which is difficult to avoid with methods like boxcar smoothing, as Wiener filtering by construction maintains a balance between data fidelity and prior knowledge. Although our presentation and algorithm are currently oriented toward LFI processing, the general procedure is fully generalizable to other experiments, as long as the Solar dipole signal is available to be used for calibration.
Continuous streamflow prediction is crucial in many applications of water resources planning and management. However, streamflow prediction is challenging, particularly in data-scarce regions. Here, we demonstrate an approach to regionalize the flow duration curve for predicting daily streamflow in the data-scare region of the central Himalayas. We developed a regression-based model to estimate streamflow at various segments of a flow duration curve by incorporating basin characteristics and climate variables. This study analyzes the sensitivities of proximity and characteristics between the donor (gauged) and receptor (ungauged) basins for time-series streamflow prediction. Our results show that regionalization techniques perform better in low to medium flows over high flows. Our findings are significant in the central Himalayan regional context to inform operational and management decisions in water sector projects like hydropower plants, which generally rely on low-to-medium streamflow information. Although the quantitative results are region-specific, the approach and insights are generalizable to the Himalayan region.
Aging water system infrastructure, including drinking water, wastewater, and stormwater, poses a growing challenge for utilities and municipalities. These water systems have well documented challenges with respect to their age, condition, and level of service. ASCE annual report cards consistently rate these infrastructure systems in the United States as underfunded, overcapacity, or past service life (ASCE 2025 Report Card). For example, Chini and Stillwell (2017) estimated that the mean water loss in drinking water systems, i.e., non-revenue water, is approximately 16% across the United States. These concerns are not just relegated to the United States, with Courtenay, British Columbia, identifying 17% of their water main pipes as in a ‘poor’ condition state, defined as a category condition 5 out of 5 (City of Courtenay, 2024). These cases illustrate the challenges utilities are facing to manage extensive networks of infrastructure to deliver a consistent and high level of service. For buried infrastructure such as water systems, studies suggest that preventative interventions can lead to lower maintenance costs and fewer service disruptions (Mazumder et al, 2018; Li et al, 2014). The demonstrated need and benefit of appropriately applied asset management is juxtaposed against the relatively sparse literature that evaluates water systems within an asset management construct. Since 2020, just 37 papers specifically reference asset management in the Journal of Water Resources Planning and Management. Of those, only a few specifically look to develop strategies for improved asset management. Therefore, we highlight four key research areas that represent opportunities for advancement of asset management research for water systems. First, advances in condition assessment and forecasting are needed to better estimate asset deterioration using diverse datasets. Second, machine learning (ML) and artificial intelligence (AI) hold promise for predictive maintenance and investment prioritization, though questions of generalizability and model transparency remain. Third, applying a value of information framework can guide utilities in making cost-effective sensor deployment and data collection decisions, to direct monitoring strategies towards data-informed asset management decisions. Finally, integrated infrastructure management is critical, requiring coordinated planning with other infrastructure systems and stakeholder engagement to reduce costs and enhance service delivery.
In this manuscript, we develop multiple machine learning (ML) models to accelerate a scheme for parameterizing site-based models of exciton dynamics from all-atom configurations of condensed phase sexithiophene systems. This scheme encodes the details of a system’s specific molecular morphology in the correlated distributions of model parameters through the analysis of many single-molecule excited-state electronic-structure calculations. These calculations yield excitation energies for each molecule in the system and the network of pair-wise intermolecular electronic couplings. In this work, we demonstrate that the excitation energies can be accurately predicted using a kernel ridge regression (KRR) model with Coulomb matrix featurization. We present two ML models for predicting intermolecular couplings. The first one utilizes a deep neural network and bi-molecular featurization to predict the coupling directly, which we find to perform poorly. The second one utilizes a KRR model to predict unimolecular transition densities, which can subsequently be analyzed to compute the coupling. We find that the latter approach performs excellently, indicating that an effective, generalizable strategy for predicting simple bimolecular properties is through the indirect application of ML to predict higher-order unimolecular properties. Such an approach necessitates a much smaller feature space and can incorporate the insight of well-established molecular physics.
Here we present a novel method for learning reduced-order models of dynamical systems using nonlinear manifolds. First, we learn the manifold by identifying nonlinear structure in the data through a general representation learning problem. The proposed approach is driven by embeddings of low-order polynomial form. A projection onto the nonlinear manifold reveals the algebraic structure of the reduced-space system that governs the problem of interest. The matrix operators of the reduced-order model are then inferred from the data using operator inference. Numerical experiments on a number of nonlinear problems demonstrate the generalizability of the methodology and the increase in accuracy that can be obtained over reduced-order modeling methods that employ a linear subspace approximation.
To evaluate the high temperature sensing properties of metal oxide and perovskite materials suitable for use in combustion environments, it is necessary to understand the temperature dependence of their bandgaps. Although such temperature-driven changes can be calculated via the Allen–Heine–Cardona (AHC) theory, which assesses electron–phonon coupling for the bandgap correction at given temperatures, this approach is computationally demanding. Another approach to predict bandgap temperature-dependence is the O’Donnell model, which uses analytical expressions with multiple fitting parameters that require bandgap information at 0 K. This work employs data-driven Gaussian process regression (GPR) to predict the parameters employed in the O’Donnell model from a set of physical features. We use a sample of 54 metal oxides for which density functional theory has been performed to calculate the bandgap at 0 K, and the AHC calculations have been carried out to determine the shift in the bandgap at non-zero temperatures. As the AHC calculations are impractical for high-throughput screening of materials, the developed GPR model attempts to alleviate this issue by predicting the O'Donnell parameters purely from physical features. To mitigate the reliability issues arising from the very small size of the dataset, we apply a Bayesian technique to improve the generalizability of the data-driven models as well as quantify the uncertainty associated with the predictions. The method captures well the overall trend of the O’Donnell parameters with respect to a reduced feature set obtained by transforming the available physical features. Quantifying the associated uncertainty helps us understand the reliability of the predictions of the O’Donnell parameters and, therefore, the bandgap as a function of temperature for any novel material.
To study the mechanical behavior of mock high explosives, an experimental and simulation program was developed to calibrate, with quantified uncertainty, a material model of the bonded granular material Idoxuridine and nitroplasticized Estane-5703. This paper reports on the efficacy of such a framework as a generalizable methodology for calibrating material models against experimental data with uncertainty quantification. Additionally, this paper studies the effect of two manufacturing temperatures and three initial granular configurations on the unconfined compressive behavior of the resulting bonded granular materials. In each of these cases, the same calibration framework was used; in that, hundreds of high-fidelity direct numerical simulations using a new, graphics processing unit-enabled, high-performance finite element method software, Ratel, were run to calibrate a finite-strain phase-field fracture model against experimental data. It was found that manufacturing temperature influenced the elastic response of the mock high explosives, with higher temperatures yielding a stiffer response. By contrast, it was found that the initial configuration of the grains had a negligible impact on the overall behavior of the mock high explosives though it remains possible that local damage accumulation within the specimens could be altered by the initial configurations. Overall, the calibration framework was successful at creating well-calibrated models, showing its usefulness as an engineering and scientific tool.