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.
This project, funded by the U.S. Department of Energy under Award DE-EE0009248, advances the domestic development of fuel-cell stack technologies for heavy-duty transportation applications. The work directly supports the Million-Mile Fuel Cell Truck (M2FCT) initiative by improving performance, durability, and manufacturability—key requirements for next-generation hydrogen fuel-cell systems.
Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.
Through an Expert Match technical assistance project, Ford County requested the National Laboratory of the Rockies evaluate solar PV, wind, and battery storage options at the Western State Bank Expo Center, assessing cost savings, on-site generation, and outage resilience. This slide deck presents the findings of this analysis.
Alpha (α) eigenvalues, which describe the logarithmic time derivative of the neutron population in a multiplying system, are integral to time-dependent behavior and diagnostic applications. However, uncertainties in the evaluated nuclear data can significantly impact the accuracy of transport simulations for such quantities. This work explores the use of machine learning models to predict two key outputs, α-eigenvalues and keff bias, using input features derived from α-eigenvalue sensitivities to nuclear data. The criticality safety benchmark models used in this study come from the International Handbook of Evaluated Criticality Safety Benchmark Experiments. Three models, random forest, XGBoost, and NGBoost, are trained on both energy-resolved and energy-summed α sensitivities. For the α-eigenvalue bias prediction, NGBoost achieved the highest R 2 (0.9476) using energy-resolved features, while XGBoost performed best using summed sensitivities. In contrast, when predicting the keff bias, all the models showed moderate predictive capability (best R 2 ≈ 0.72), as the mapping from the static α-sensitivities to the static keff bias was less direct. SHAP (SHapley Additive exPlanations) analysis was used to interpret the model predictions. Across both prediction tasks, the features associated with neutron capture [H-1 (n, γ)], uranium scattering reactions (such as 235 U elastic/inelastic), and actinide capture/fission reactions (such as 239 Pu and 234 U) were consistently identified as the most impactful. This highlights the key role of specific nuclear reactions and energy ranges in shaping both time-dependent and steady-state criticality behavior. These results demonstrated that α-sensitivities, despite being computed for time-dependent metrics, can provide valuable insights for predicting both α-eigenvalues and the keff bias. Moreover, machine learning models offer a promising pathway for uncovering important nuclear data dependencies and guiding future data evaluation efforts.
Crosslinking polymers is a common approach to create mechanically stable solid materials such as polymer electrolytes for lithium batteries. In conventional liquid electrolytes, the solvent molecules move freely to accommodate the field-induced motion of ions. However, in crosslinked polymer electrolytes, the rearrangement of polymer segments is constrained by the deformation limits of the network. Herein, we develop a new transport model that accounts for both the formation of concentration gradients and the elasticity of the electrolyte. The elasticity is incorporated by adding an additional term related to the entropy of crosslinked strands to the electrochemical potential of the salt. The resulting Crosslink Model contains two adjustable parameters: $\mathcal{N}$, the average number of monomers in a strand, and λ crit , the maximum strain the network can sustain. These solid-like constraints produce singularities in the governing transport equations, fundamentally altering the concentration profiles. Plateaus in salt concentrations emerge near the electrodes, and network elasticity introduces a strain overpotential. When compared to a Baseline Model ($\mathcal{N}$ → ∞, equivalent to concentrated solution theory), which predicts steepest gradients near the electrodes, both models yield similar current–voltage relationships. Model predictions are compared to electrochemical data for a poly(ethylene oxide)-based crosslinked polymer electrolyte.
Rapidly growing capabilities of quantum simulators to probe quantum many-body phenomena require new methods to characterize increasingly complex states. Here, we present a protocol that constrains quantum states using experimentally measured correlation functions. This method enables measurement of a quantum state’s entanglement structure, opening a new route to study entanglement-related phenomena. Our approach extends Gaussian state parameterizations by systematically incorporating higher-order correlations. We show the protocol’s usefulness in conjunction with current and forthcoming experimental capabilities, focusing on weakly interacting fermions as a proof of concept. Here, the lowest nontrivial expansion quantitatively predicts early time thermalization dynamics, including signaling the onset of quantum chaos indicated by the entanglement Hamiltonian.
Superdiffusive spin transport in the one-dimensional (1D) Heisenberg model is a key theoretical discovery in nonequilibrium quantum many-body physics. Although extensively studied in 1D systems, the breakdown and sustenance of superdiffusion in two-dimensional (2D) lattices with integrability-breaking terms, as found in real materials, remains an open question. To address this, we develop a toy model that extends the 1D Heisenberg model with a representative set of 2D interaction types and tunable strengths. Our model exhibits varying degrees of superdiffusion breakdown depending on the interaction type, spanning ballistic to diffusive regimes. We establish and justify a hierarchy of 2D interactions based on their resilience against superdiffusion breakdown: Heisenberg >𝑋𝑋 > Ising. This precise control over the superdiffusive behavior also enables rigorous benchmarking of quantum hardware, and our simulations on IBM's Heron devices confirm the hardware's ability to accurately capture these many-body nonequilibrium phenomena. Overall, our results are relevant not only to simulating superdiffusion in real materials, such as the 1D Heisenberg compound KCuF3, which contains modest nonintegrable 2D terms, but also to extending superdiffusive behavior to larger 2D qubit lattices and other 2D materials.
Vegetation processes are fundamentally limited by nutrient and water availability, the uptake of which is mediated by plant roots in terrestrial ecosystems. While tropical forests play a central role in global water, carbon, and nutrient cycling, we know very little about tradeoffs and synergies in root traits that respond to resource scarcity. Tropical trees face a unique set of resource limitations, with rock-derived nutrients and moisture seasonality governing many ecosystem functions, and nutrient versus water availability often separated spatially and temporally. Root traits that characterize biomass, depth distributions, production and phenology, morphology, physiology, chemistry, and symbiotic relationships can be predictive of plants’ capacities to access and acquire nutrients and water, with links to aboveground processes like transpiration, wood productivity, and leaf phenology. In this review, we identify an emerging trend in the literature that tropical fine root biomass and production in surface soils are greatest in infertile or sufficiently moist soils. We also identify interesting paradoxes in tropical forest root responses to changing resources that merit further exploration. For example, specific root length, which typically increases under resource scarcity to expand the volume of soil explored, instead can increase with greater base cation availability, both across natural tropical forest gradients and in fertilization experiments. Also, nutrient additions, rather than reducing mycorrhizal colonization of fine roots as might be expected, increased colonization rates under scenarios of water scarcity in some forests. Efforts to include fine root traits and functions in vegetation models have grown more sophisticated over time, yet there is a disconnect between the emphasis in models characterizing nutrient and water uptake rates and carbon costs versus the emphasis in field experiments on measuring root biomass, production, and morphology in response to changes in resource availability. Closer integration of field and modeling efforts could connect mechanistic investigation of fine-root dynamics to ecosystem-scale understanding of nutrient and water cycling, allowing us to better predict tropical forest-climate feedbacks.
Cerium hydride has a variety of interesting properties, including a known lattice contraction and densification with increasing hydrogen content. However, precise stoichiometric control is not experimentally straightforward and ab initio approaches are not computationally feasible for many properties such as melting and low temperature diffusion. Therefore, we develop a machine-learned interatomic potential for cerium hydride that is valid for H to Ce ratios from 2.0 to 3.0. A query-by-committee active learning approach is used to develop the training set. Leveraging classical molecular dynamics simulations, we assess a range of properties and provide fundamental mechanisms for the trends with stoichiometry. Finally, a majority of the properties follow the trend of lattice contraction, being governed by the stronger lattice binding induced by adding octahedral atoms.
High-energy hadronic collisions generate environments characterized by temperatures above 100 MeV (refs.), about 100,000 times hotter than the centre of the Sun. At present, it is therefore unclear how light (anti)nuclei with mass number A of a few units, such as the deuteron, 3 He or 4 He, each bound by only a few MeV, can emerge from these collisions. Here, the ALICE Collaboration reports that deuteron–pion momentum correlations in proton–proton (pp) collisions provide model-independent evidence that about 90% of the observed (anti)deuterons are produced in nuclear reactions following the decay of short-lived resonances, such as the Δ(1232). These findings, obtained by the ALICE Collaboration at the Large Hadron Collider, resolve a gap in our understanding of nucleosynthesis in ultrarelativistic hadronic collisions. Apart from offering insights on how (anti)nuclei are formed in hadronic collisions, the results can be used in the modelling of the production of light and heavy nuclei in cosmic rays and dark-matter decays.
Low-temperature electrochemical cells reducing iron oxides to metal in alkaline electrolytes can support fully electrified steelmaking processes. Previous studies on these cells have primarily focused on high-surface-area hematite, Fe 2 O 3 , reactants whereas attempts to reduce suspensions of magnetite, Fe 3 O 4 —one of the two feedstocks for existing ironmaking reactors—have generally been limited to low rates of reaction (<30 mA cm −2 ). Here, in this study, we control the crystalline domain size of Fe 2 O 3 and Fe 3 O 4 particles in 10 M NaOH electrolytes to study how the nanoscale morphology of oxides controls the rate of electrochemical ironmaking. Rotating-ring disk electrode measurements of Fe 2+ , in situ Raman spectroscopy of the electrode surface, and ex situ electron microscopy were consistent with a hypothesized passivation process at Fe 3 O 4 surfaces that may prevent the continuous formation of soluble intermediates. Sufficiently small (<100 nm diameter) oxide particles yielded Fe partial current densities >160 mA cm −2 , a fivefold increase relative to previously reported rates for Fe 3 O 4 suspensions and comparable to active Fe 2 O 3 . Electron microscopy revealed that electrodeposited films were composed of micron-scale crystalline Fe domains with a porous film of Fe 3 O 4 nanoparticles and supports a model where Fe is grown primarily from soluble Fe 2+ intermediates. Based on these insights, inactive blast-furnace-grade iron-oxides were transformed into high surface area nanoparticles (1.6 to 229.2 m 2 g −1 ) via reprecipitation, leading to a ninefold enhancement in faradaic efficiency and an Fe partial current density of 120 mA cm −2 . When integrated with chlor-iron cells producing reagents for reprecipitation, this approach could lead to a cost-competitive process for electrochemical ironmaking from industrially relevant feedstocks.
Solving power system capacity expansion planning (CEP) problems at realistic spatial resolutions is computationally challenging. Thus, a common practice is to solve CEP over zonal models with low spatial resolution rather than over full-scale nodal power networks. Due to improvements in solving large-scale stochastic mixed integer programs, these computational limitations are becoming less relevant, and the assumption that zonal models are realistic and useful approximations of nodal CEP is worth revisiting. Here, this work is the first to conduct a systematic computational study on the assumption that spatial aggregation can reasonably be used for ISO-scale CEP. By considering a realistic, large-scale test network based on the state of California with over 8000 buses, we find that well-designed small spatial aggregations can yield good approximations but that coarser zonal models may result in large distortions of investment decisions, e.g., capacity under-investment of up to 41% for the lowest resolution model considered.
Abstract This paper reports global nonlinear gyrokinetic simulations that couple meso-scale reversed shear Alfvén eigenmodes (RSAEs) driven by energetic particles (EPs) and ion temperature gradient (ITG) microturbulence driven by thermal plasma, using equilibrium and profiles from DIII-D discharge #159243. In simulations focusing only on the ITG, electrostatic ITG drives a huge thermal ion heat transport, which is reduced by a factor of 10 to a level close to the experimental value in electromagnetic simulation due to finite β effect. In the simulations coupling the RSAE and ITG, ITG can scatter the resonant EP nonlinearly trapped by the RSAE and damp the zonal flows generated by the RSAE. The regulation of the RSAE by the ITG greatly reduces the initial saturation amplitude of the RSAE but increases the RSAE amplitude and associated EP transport to experimental levels in the quasi-steady state. The RSAE effects on the ITG, specifically the stronger zonal flows generated by the RSAE and the RSAE frequency modulation of the ITG-induced thermal ion heat transport, in turn, leads to a reduction of the thermal ion heat transport by more than a factor of 2 . For a stronger background ITG, the regulation of the RSAE by the ITG is stronger, while the RSAE effects on the ITG are weaker. This work highlights the importance of cross-scale coupling in the dynamics of the AE turbulence and EP transport.
Magnetic kagome materials provide a platform for exploring magneto-transport phenomena, symmetry breaking and charge ordering driven by the intricate interplay among electronic structure, topology and magnetism. Yet geometric frustration in conventional kagome magnets limits their tunability. Here we propose a design strategy for interweaving quasi-one-dimensional magnetic Tb zigzag chains with non-magnetic Ti-based kagome bilayers in TbTi 3 Bi 4 . Comprehensive spectroscopic analyses reveal coexisting elliptical-spiral magnetic and spin-density-wave orders accompanied by a large ~90 meV band-folding gap. The combined magnetic and electronic state leads to a giant anomalous Hall conductivity of 10 5 Ω −1 cm −1 , which exceeds that observed in frustrated kagome analogues. These results establish TbTi 3 Bi 4 as a model system of magnetic kagome metals with strong electron–magnetism interactions and underscore the necessity of interweaving designed magnetic and charge layers separately to achieve tunable transport properties. This design strategy will enable the discovery of emergent quantum states and next-generation electronic materials.
A controls system for a range-extended electric vehicle comprising an overall system control unit, an engine control module configured to control a range extender of the range-extended electric vehicle, power electronics configured to control a generator of the range-extended electric vehicle, and a supervisory control module coupled between the overall system control unit and the engine control module and the power electronics, the supervisory control module configured to receive information from the overall system control unit and provide commands to the engine control module and the power electronics.
Glancing Angle Deposition (GLAD) has emerged as a versatile and powerful nanofabrication technique for developing next-generation gas sensors by enabling precise control over nanostructure geometry, porosity, and material composition. Through dynamic substrate tilting and rotation, GLAD facilitates the fabrication of highly porous, anisotropic nanostructures, such as aligned, tilted, zigzag, helical, and multilayered nanorods, with tunable surface area and diffusion pathways optimized for gas detection. This review provides a comprehensive synthesis of recent advances in GLAD-based gas sensor design, focusing on how structural engineering and material integration converge to enhance sensor performance. Key materials strategies include the construction of heterojunctions and core–shell architectures, controlled doping, and nanoparticle decoration using noble metals or metal oxides to amplify charge transfer, catalytic activity, and redox responsiveness. GLAD-fabricated nanostructures have been effectively deployed across multiple gas sensing modalities, including resistive, capacitive, piezoelectric, and optical platforms, where their high aspect ratios, tailored porosity, and defect-rich surfaces facilitate enhanced gas adsorption kinetics and efficient signal transduction. These devices exhibit high sensitivity and selectivity toward a range of analytes, including NO2, CO, H2S, and volatile organic compounds (VOCs), with detection limits often reaching the parts-per-billion level. Emerging innovations, such as photo-assisted sensing and integration with artificial intelligence for data analysis and pattern recognition, further extend the capabilities of GLAD-based systems for multifunctional, real-time, and adaptive sensing. Finally, current challenges and future research directions are discussed, emphasizing the promise of GLAD as a scalable platform for next-generation gas sensing technologies.
Abstract A detailed investigation addressing the effects of functionalizing conjugated polymers with oligo(ethylene glycol) (EG n ) sidechains on the performance and polymer‐electrolyte compatibility of electrochromic devices (ECDs) is reported. The electrochemistry for a series of donor‐acceptor copolymers having near‐infrared (NIR)‐optical absorption, where the donor fragment is 3,4‐ethylenedioxythiophene (EDOT) or an EG n functionalized bithiophene (g2T) and the acceptor fragment is diketopyrrolopyrrole (DPP) functionalized with branched alkyl or EG n sidechains, is extensively probed. ECDs are next fabricated and it is found that EG n sidechain incorporation must be finely balanced to promote polymer‐electrolyte compatibility and provide efficient ion exchange. Proper electrolyte‐cation pairing and polymer structural tuning affords a 2x increase in optical contrast (from 12% to 24%) and >60x reduction in switching time (from 20 to 0.3 s). Atomic force microscopy (AFM)/grazing incidence wide‐angle X‐ray scattering (GIWAXS) characterization of the polymer film morphology/microstructure reveals that an over‐abundance of EG n sidechains generates large polymer crystallites, which can suppress ion exchange. Lastly, time‐of‐flight secondary ion mass spectrometry (ToF‐SIMS) indicates sidechain/electrolyte identity does not influence the electrolyte penetration depth into the films, and EG n sidechain inclusion increases electrolyte cation uptake. The material structural design insight and guidelines regarding the polymer‐electrolyte ion insertion/expulsion dynamics reported here should be of significant utility for developing next‐generation mixed ionic‐electronic conducting materials.
Lithium- and manganese-rich layered oxides (LMR) stand out as next-generation lithium-ion cathode chemistries, which harness both transition-metal and lattice-oxygen redox processes to deliver exceptional capacity and energy density. However, their full potential is hindered by intrinsic oxygen instability and structural degradation, resulting in pronounced voltage fade and capacity decay. Here, we present a complex-concentrated anion-doping paradigm in which multiple anions, F, Br, and S, are incorporated into the oxygen sublattice to enhance oxygen-redox and structural stability. X-ray absorption spectroscopy and aberration-corrected scanning transmission electron microscopy confirm ultra-stable local oxygen coordination environments during long-term cycling, with detrimental phase transformations and oxygen-loss-induced cavitation dramatically inhibited. Notably, we show that the characteristic LiTM6 transition metal (TM) honeycomb ordering is preserved even after electrochemical cycling. Concurrently, this strategy yields an unprecedented volume change of only 0.63% upon charging to 4.8 V vs. Li+/Li, achieving the first zero-strain LMR cathode. The resulting LMR cathode delivers ultralow voltage fade (1 mV per cycle during the first 100 cycles and becomes negligible in subsequent cycles) and outstanding energy retention (93% after 200 cycles) in a pouch cell configuration. Our complex-concentrated anion-doping concept establishes a broadly applicable strategy for resolving chemo-mechanical failure mechanisms in ceramic intercalation electrodes for next-generation energy storage.