Search NASA⌕ Search

SEARCH · Search NASA

Results for “predictive thermodynamics”

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 433 records · Page 24

Large Destabilization of (TiVNb)-Based Hydrides via (Al, Mo) Addition: Insights from Experiments and Data-Driven Models

High-entropy alloys (HEAs) represent an interesting alloying strategy that can yield exceptional performance properties needed across a variety of technology applications, including hydrogen storage. Examples include ultrahigh volumetric capacity materials (BCC alloys → FCC dihydrides) with improved thermodynamics relative to conventional high-capacity metal hydrides (like MgH 2 ), but still further destabilization is needed to reduce operating temperature and increase system-level capacity. Here, in this work, we demonstrate efficient hydride destabilization strategies by synthesizing two new Al 0.05 (TiVNb) 0.95–x Mo x (x = 0.05, 0.10) compositions. We specifically evaluate the effect of molybdenum (Mo) addition on the phase structure, microstructure, hydrogen absorption, and desorption properties. Both alloys crystallize in a bcc structure with decreasing lattice parameters as the Mo content increases. The alloys can rapidly absorb hydrogen at 25 °C with capacities of 1.78 H/M (2.79 wt %) and 1.79 H/M (2.75 wt %) with increasing Mo content. Pressure-composition isotherms suggest a two-step reaction for hydrogen absorption to a final fcc dihydride phase. The experiments demonstrate that increasing Mo content results in a significant hydride destabilization, which is consistent with predictions from a gradient boosting tree data-driven model for metal hydride thermodynamics. Furthermore, improved desorption properties with increasing Mo content and reversibility were observed by in situ synchrotron X-ray diffraction, in situ neutron diffraction, and thermal desorption spectroscopy.

36 MATERIALS SCIENCE↗

Spontaneous Symmetry–Breaking of Nonequilibrium Steady–States Caused by Nonlinear Electrical Transport

Negative differential resistance (NDR) in certain materials has been attributed to spontaneous emergence of symmetry-breaking electrical current density localization from a previously homogeneous distribution, which is postulated to occur due to the nonequilibrium thermodynamic force of minimization of entropy production. However, this phenomenon has not been quantitatively predicted based on intrinsic material properties and an applied electrical stimulus. Herein an instability criterion is derived for localization of current density and temperature from a thermal fluctuation in a parallel conductor model of a thin film that is subject to Newton's law of cooling. The conditions for steady–state electro-thermal localization is predicted, verifying a decrease in entropy production upon localization. Electro-thermal localization accompanied by a decrease of entropy production is confirmed in a multiphysics simulation of current flow in a thin film. The instability criterion predicts conditions for spontaneous current density localization, relating symmetry breaking fundamentally to dynamical instability via Local Activity theory.

36 MATERIALS SCIENCE↗

Theory-guided design of duplex-phase multi-principal-element alloys

Density-functional theory (DFT) is used to identify phase-equilibria in multi-principal-element and high-entropy alloys (MPEAs/HEAs), including duplex-phase and eutectic microstructures. Here, a combination of composition-dependent formation energy and electronic-structure-based ordering parameters were used to identify a transition from FCC to BCC favoring mixtures, and these predictions experimentally validated in the Al-Co-Cr-Cu-Fe-Ni system. A sharp crossover in lattice structure and dual-phase stability as a function of composition were predicted via DFT and validated experimentally. The impact of solidification kinetics and thermodynamic stability was explored experimentally using a range of techniques, from slow (castings) to rapid (laser remelting), which showed a decoupling of phase fraction from thermal history, i.e., phase fraction was found to be solidification rate-independent, enabling tuning of a multi-modal cell and grain size ranging from nanoscale through macroscale. Strength and ductility tradeoffs for select processing parameters were investigated via uniaxial tension and small-punch testing on specimens manufactured via powder-based additive manufacturing (directed-energy deposition). This work establishes a pathway for design and optimization of next-generation multiphase superalloys via tailoring of structural and chemical ordering in concentrated solid solutions.

36 MATERIALS SCIENCE↗

Predictive Chemical Kinetic Modeling: Where We Succeed, Where We Struggle, and What Comes Next

Chemical kinetic modeling plays a foundational role in fields ranging from energy to environmental science, pharmaceuticals, and advanced materials. The past two decades have seen remarkable progress, particularly in modeling gas-phase reactions for thermochemical processes, leading to impactful industrial applications such as steam cracking and air quality management. However, new challenges are emerging. The successful development of systematic methodologies for the description of gas-phase kinetics opens the possibility to apply the same approach to the study of more challenging systems. Here, we review recent advances, including ab initio transition state theory-based master equation estimation of elementary rates, automated mechanism generation, machine-learning-assisted kinetics, and uncertainty quantification, and discuss the advances needed to apply the same methodological approach in areas such as heterogeneous catalysis, electrochemistry, liquid-phase and solid-state reactivity, and multiscale model integration. We advocate for the development of targeted tools, especially methods that go beyond empirical tuning toward first-principles-based predictions. We highlight the need for accessible software and AIaugmented workflows to democratize modeling for industry and academia alike. In this perspective, we call attention to not only what has worked but also what remains unsolved, advocating to avoid overemphasizing successes in scientific works at the expense of realism. The next decade should focus on predictive capability, physical accuracy, and community infrastructure (e.g., databases and services) to enable innovation across diverse fields. We argue that kinetic modeling, properly equipped, can accelerate discovery far beyond its traditional domains.

ab initio calculations↗

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid–base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Robustness and Eventual Slow Decay of Bound States of Interacting Microwave Photons in the Google Quantum AI Experiment

Integrable models are characterized by the existence of stable excitations that can propagate indefinitely without decaying. This includes multimagnon bound states in the celebrated 𝑋⁢𝑋⁢𝑍 spin-chain model and its integrable Floquet counterpart. A recent Google Quantum AI experiment [A. Morvan et al., Nature 612, 240 (2022)] realizing the Floquet model has demonstrated the persistence of such collective excitations even when the integrability is broken: this observation is at odds with the expectation of ergodic dynamics in generic nonintegrable systems. Here, we study the spectrum of the model realized in the experiment using exact diagonalization and physical arguments. We find that isolated bands corresponding to the descendants of the exact bound states of the integrable model are clearly observable in the spectrum for a large range of system sizes. However, our numerical analysis of the localization properties of the eigenstates suggests that the bound states become unstable in the thermodynamic limit. A perturbative estimate of the decay rate agrees with the prediction of an eventual instability for large system sizes.

Exact diagonalization↗

FROMP-from-Monomer

This repository provides the 1D reaction-diffusion simulation engine for the multi-scale framework. It takes monomer-level kinetic and thermodynamic descriptors as inputs, approximated from quantum chemical calculations, to simulate and predict the macroscopic propagation behavior of Frontal Ring-Opening Metathesis Polymerization (FROMP). By solving a coupled 1D reaction-diffusion PDE system, this tool predicts front propagation speed, front temperature, and degree of conversion. The solver explicitly captures the three fundamental ROMP steps, initiator activation/inhibition, chain initiation, and propagation, while accounting for high-temperature cycloreversion as a competing pathway. The numerical implementation uses FEniCS to handle the coupled heat diffusion and reaction kinetics.

Chua, Lauren [Massachusetts Inst. of Technology (M↗

The Impact of Nonzenith Elevation Angles on Ground-Based Infrared Thermodynamic Retrievals

Observing thermodynamic profiles within the planetary boundary layer is essential to understanding and predicting atmospheric phenomena because of the significant exchange of sensible and latent heat between the land and atmosphere within that layer. The Atmospheric Emitted Radiance Interferometer (AERI) is a ground-based infrared spectrometer used to obtain the vertical profiles of temperature and water vapor mixing ratio. Most AERIs are only capable of zenith views, although the Marine AERI (M-AERI) has a design that allows it to view various elevation angles. In this study, we quantify the improvement in the information content and accuracy of the retrieved profiles when nonzenith angles are included, as is common with microwave radiometer profilers. The impacts of the additional scan angles are quantified through both a synthetic study and with M-AERI observations from the ARM Cloud Aerosol Precipitation Experiment (ACAPEX) campaign. Furthermore, the simulation study shows that low elevation angles contain more information content for temperature whereas high elevation angles have more information content for water vapor. Outside of very humid environments, the addition of low elevation angles also results in lower root-mean-square errors when compared with high angles for both temperature and water vapor mixing ratio, although this is primarily a result of averaging multiple observations together to reduce instrument noise. Real-world results from the ACAPEX dataset indicate similar results as were found for the simulation study, although not all predicted benefits are realized because of the small sample size and observational uncertainties.

54 ENVIRONMENTAL SCIENCES↗

STAT7 v2.0 User Guide

The STAT7 software was developed to perform steady-state, single-phase thermal hydraulics analysis of plate-fueled reactors based on statistical propagation of uncertainties. Application of the software includes non-power research and test reactor analysis, and it has been used for the conversion to low- enriched uranium fuel of U.S. High Performance Research Reactors such as the Massachusetts Institute of Technology Research Reactor. Since it can be necessary to repeat reactor safety analysis, such as during fuel reloading, STAT7 accommodates flexibility in analyzing many practical aspects of reactor fuel management. STAT7 uses a Monte Carlo approach to model uncertainty in common fuel fabrication parameters and other key reactor operating parameters required for reactor thermal hydraulics analysis. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux, or onset of flow instability. STAT7 supports water properties based on the IAPWS-IF97 functions (The International Association for the Properties of Water and Steam Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam) in addition to fitted functions. STAT7 predicts axial profiles of fuel, cladding, and coolant temperature along a lateral stripe that runs the full length of the fuel plate from the bottom to the top. STAT7 can simultaneously analyze every axial node in each lateral stripe of all fuel plates and coolant channels in every fuel element of an entire reactor core. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. In a single execution, STAT7 can be used to perform full core analysis by analyzing the margin to onset of nucleate boiling and onset of flow instability for each axial node of each stripe of each plate of each fuel element in the core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High entropy oxides prediction and discovery by the Mixed Enthalpy-Entropy Descriptor

The vast, high-dimensional composition space of high-entropy oxides (HEOs) offers exceptional opportunities for functional materials discovery, yet it also poses a fundamental challenge: the rational and efficient prediction of stable, synthesizable compositions and the corresponding structure–property relationships. Despite growing interest, the field still lacks broadly applicable, physically grounded descriptors capable of navigating various large chemical spaces. Here, we introduce a Mixed Enthalpy–Entropy Descriptor (MEED) that enables rapid, first-principles–based prediction of HEOs synthesizability across diverse chemistries. Using MEED, we perform high-throughput screening of two distinct HEO families: rocksalt oxides and perovskite oxides. The predicted top candidates in each family were experimentally validated. MEED reveals unifying thermodynamic and structural principles governing stability across both chemical compositions and polymorphs, providing mechanistic insight into the formation of high-entropy phases. This work significantly broadens the accessible chemical design space for HEOs and establishes a data-efficient framework for accelerating the discovery of next-generation functional materials.

Yu, Liping [University of Central Florida]↗

Interaction effects in graphene in a weak magnetic field

A weak perpendicular magnetic field, B, breaks the chiral symmetry of each valley in the electron spectrum of graphene, preserving the overall chiral symmetry in the Brillouin zone. Here we explore the consequences of this symmetry breaking for the interaction effects in graphene. In particular, we demonstrate that the electron-electron interaction lifetime acquires an anomalous B dependence. Also, the ballistic zero-bias anomaly, δv(ω), where ω is the energy measured from the Fermi level, emerges at a weak B and has the form δv(B) ~ B 2 /ω 2 . Temperature dependence of the magnetic-field corrections to the thermodynamic characteristics of graphene is also anomalous. We discuss experimental manifestations of the effects predicted. The microscopic origin of the B-field sensitivity is an extra phase acquired by the electron wave function resulting from the chirality-induced pseudospin precession.

36 MATERIALS SCIENCE↗

Thermodynamic origin of nonvolatility in resistive memory

Electronic switches based on the migration of high-density point defects, or memristors, are poised to revolutionize post-digital electronics. Despite significant research, key mechanisms for filament formation and oxygen transport remain unresolved, hindering our ability to predict and design device properties. For example, experiments have achieved 10 orders of magnitude longer retention times than predicted by current models. Here, using electrical measurements, scanning probe microscopy, and first-principles calculations on tantalum oxide memristors, we reveal that the formation and stability of conductive filaments crucially depend on the thermodynamic stability of the amorphous oxygen-rich and oxygen-poor compounds, which undergo composition phase separation. Including the previously neglected effects of this amorphous phase separation reconciles unexplained discrepancies in retention and enables predictive design of key performance indicators such as retention stability. Furthermore, this result emphasizes non-ideal thermodynamic interactions as key design criteria in post-digital devices with defect densities substantially exceeding those of today’s covalent semiconductors.

36 MATERIALS SCIENCE↗

Air separation via a two-step solar thermochemical cycle based on (Ba,La) x Sr 1-x FeO 3-δ : Thermodynamic analysis

A two-step solar thermochemical cycle was considered for air separation to produce N2 based on (Ba,La) x Sr 1-x FeO 3-δ perovskite reduction/oxidation (redox) reactions for A-site fractions of 0 ≤ x ≤ 0.2. The cycle steps encompassed (1) thermal reduction and O 2 release via concentrated solar input and (2) re-oxidation with air to uptake O 2 and produce high-purity N 2 . Thermogravimetry at temperatures between 400 and 1100 °C in atmospheres of 0.005 to 90% O 2 /Ar at 1 bar was performed to measure equilibrium nonstoichiometries. The compound energy formalism was applied to model redox thermodynamics for both Ba 2+ and La 3+ substitution. Non-linear regression was used to determine the empirical parameters based on the thermogravimetric measurements. The model was used to define partial molar reaction enthalpies and entropies and predicted equilibrium oxygen nonstoichiometry as functions of oxide stoichiometry, site fraction, temperature, and O 2 partial pressure. The thermodynamic analysis showed the materials are appealing for air separation at temperatures below 800 °C.

42 ENGINEERING↗

Grain growth kinetics of the gamma phase metallic uranium

We report metallic uranium is a leading fuel form for sodium cooled fast reactors as an enabling technology of future nuclear energy systems. Mechanistic understanding of fuel behaviors and kinetics under thermodynamic equilibrium and highly non-equilibrium conditions are essential for evaluating fuel performance. It is important to understand and predict the grain and pore evolutions of metallic fuels under thermal and irradiation conditions. However, very limited data are available on the grain growth kinetics and mechanisms of pure gamma phase uranium. In this paper, the pure gamma uranium pellets with different grain structures were fabricated by combining high-energy ball milling and spark plasma sintering. Isothermal annealing tests were performed to investigate the grain growth behavior of the pure gamma phase uranium with different initial grain sizes. A parabolic relationship in grain growth with time was identified for the submicron-sized (374 nm) sample. In contrast, for the nano-sized (137 nm) sample, the grain growth shows a linear relationship with time. The activation energies of grain growth were determined as 199.5 KJ/mol and 80.6 KJ/mol for nano-sized and submicron-sized grain structures, respectively. For the nano-sized sample, the rate-control step of grain growth is dominated by the triple-junction migration, in which the grain boundary triple junction drags the grain growth, leading to a higher activation energy than the bulk diffusion. The dominating mechanism for the submicron-sized sample is grain boundary diffusion. The mechanistic understanding and critical data obtained on the kinetics of pure uranium phases will be useful to evaluate fuel behavior under thermodynamic equilibrium conditions and develop a high fidelity model to predict fuel performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of a kinetic-thermodynamic model for lime-stabilization of Na-bentonite

This study presents the first kinetic model to predict the solid and pore solution composition of Na-bentonite clay reacting with slaked lime over a period of 720 days. The model successfully accounts for most experimental data using a single kinetic rate constant. The following sequence of reactions was predicted by the model: initial rapid dissolution of portlandite within the first 7 days, leading to a decrease in pH and dissolved calcium, and concurrent formation of calcium silicate hydrates (C-S-H: jennite), calcium aluminate hydrate (C-A-H: C₄AH₁₃), calcium aluminosilicate hydrates (stratlingite) and hydrotalcite. After 7 days, jennite and stratlingite are predicted to transform into tobermorite-II, contributing to strength development up to 28 days. From 28 to 90 days, continued montmorillonite dissolution is predicted, along with minor formation of ettringite, partial tobermorite-II dissolution, and precipitation of secondary phases such as albite and talc. Experimentally, portlandite dissolution was confirmed by TGA and XRD and found to be complete within 7 days, in agreement with model predictions. However, other predicted solid-phase transformations (e.g., tobermorite-II formation and dissolution, ettringite, albite, and talc formation) could not be conclusively verified through experimental techniques. Aqueous phase measurements confirmed that the pH and Ca trends in solution, and that equilibrium was reached by 90 days.

Chemical kinetics↗

Temperature and Solvent Effects on H 2 Splitting and Hydricity: Ramifications on CO 2 Hydrogenation by a Rhenium Pincer Catalyst

The catalytic hydrogenation of carbon dioxide holds immense promise for applications in sustainable fuel synthesis and hydrogen storage. Mechanistic studies that connect thermodynamic parameters with the kinetics of catalysis can provide new understanding and guide predictive design of improved catalysts. Reported here are thermochemical and kinetic analyses of a new pincer-ligated rhenium complex ( tBu POCOP)Re(CO) 2 ( tBu POCOP = 2,6-bis(di-tert-butylphosphinito)phenyl) that catalyzes CO 2 hydrogenation to formate with faster rates at lower temperature. Because the catalyst follows the prototypical “outer sphere” hydrogenation mechanism, comprehensive studies of temperature and solvent effects on the H 2 splitting and hydride transfer steps are expected to be relevant to many other catalysts. Strikingly large entropy associated with cleavage of H2 results in a strong temperature dependence on the concentration of [( tBu POCOP)Re(CO) 2 H] – present during catalysis, which is further impacted by changing the solvent from toluene to tetrahydrofuran to acetonitrile. New methods for determining the hydricity of metal hydrides and formate at temperatures other than 298 K were developed, providing insight into how temperature can influence the favorability of hydride transfer during catalysis. These thermochemical insights guided the selection of conditions for CO 2 hydrogenation to formate with high activity (up to 364 h –1 at 1 atm or 3330 h –1 at 20 atm of 1:1 H 2 CO 2 ). In cases where hydride transfer is the highest individual kinetic barrier, entropic contributions to outer sphere H 2 splitting lead to a unique temperature dependence: catalytic activity increases as temperature decreases in tetrahydrofuran (200-fold increase upon cooling from 50 to 0 °C) and toluene (4-fold increase upon cooling from 100 to 50 °C). Furthermore, the ramifications on catalyst structure-function relationships are discussed, including comparisons between “outer sphere” mechanisms and metal–ligand cooperation mechanisms.

14 SOLAR ENERGY↗

Tasting nuclear pasta made with classical molecular dynamics simulations

Nuclear clusters or voids in the inner crust of neutron stars were predicted to have various shapes collectively nicknamed nuclear pasta. Furthermore, the recent review by López, Dorso and Frank summarized their systematic investigations into properties especially the morphological and thermodynamical phase transitions of the nuclear pasta within a Classical Molecular Dynamics model, providing further stimuli to find more observational evidences of the predicted nuclear pasta in neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Predicting receptor-ligand pairing preferences in plant-microbe interfaces via molecular dynamics and machine learning

Microbiome assembly, structure, and dynamics significantly influence plant health. Secreted microbial signaling molecules initiate and mediate symbiosis by binding to structurally compatible plant receptors. For example, lipo-chitooligosaccharides (LCOs), produced by nitrogen-fixing rhizobial bacteria and various fungi, are recognized by plant lysin motif receptor-like kinases (LysM-RLKs), which activate the common symbiotic pathway. Accurately predicting these molecular interactions could reveal complementary signatures underlying the initial stages of endosymbiosis. Despite the breakthrough in protein-ligand structure prediction with deep learning-based tools, such as AlphaFold3, the large size and highly flexible nature of signaling compounds like LCOs present major challenges for detailed structural characterization and binding-affinity prediction. Typical structure-/physics-based methods of ligand virtual screening are designed for small, drug-like molecules, often rely on high-resolution, experimentally determined structures of the protein receptors, and rarely achieve sufficient sampling to obtain converged thermodynamic quantities with large ligands. In this study, we developed a hybrid molecular dynamics/machine learning (MD/ML) approach capable of predicting binding affinity rankings with high accuracy in systems involving large, flexible ligands, despite limited experimental structural information. Using coarse initial structural models, the predictions using the MD/ML workflow achieved strong alignment with experimental trends, particularly in the top-affinity tier for four legume LysM-RLKs (LYR3) binding to LCOs and a chitooligosaccharide. Furthermore, the MD-based conformation selection protocol provided critical structural insights into substrate specificity and binding mechanisms. This study demonstrates a powerful method to screen for challenging cognate ligand-receptors and advance our understanding of the molecular basis of microbial colonization in plants.

Lipo-chitooligosaccharides↗