Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non-equilibrium state”

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.

1,398 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Turbo-Design: Open-Source Radial Equilibrium Turbomachinery Solver: Part I - Turbines

Advances in 3D Geometrical Designs and Cooling have played a significant role in improving the efficiency of turbomachinery. However, these advancements must be effectively translated back to the modeler. Machine learning can facilitate this transition. Specifically, machine learning–based loss models can bridge the gap between 3D and 1D designs, enabling modelers not only to predict velocity triangles but also to extract additional geometric features. Currently, the design tools used at NASA have not been updated to support such integration—until now. TurboDesign is an open-source, Python-based framework that replaces TD2 (LEW-11029-1) and AXOD2 (LEW-16323-1), both of which are radial equilibrium solvers for axial turbines. The goal of this update is to enable the integration of machine learning loss models into radial equilibrium equations. Additionally, TurboDesign is designed to support radial machines. This paper presents the governing equations, the assumptions underlying the code, the integration of legacy loss models, an example of machine learning model integration, and a validation comparison with CFD. All code, tutorials, and documentation are available at: https://www.github.com/nasa/turbo-design

Radial Equilibrium↗
DOE OSTI2026

Polynomial-time preparation of low-temperature Gibbs states for two-dimensional toric code

In this work, we propose a polynomial-time algorithm for preparing the Gibbs state of the two-dimensional toric code Hamiltonian at any temperature, starting from any initial state, significantly improving upon prior estimates that suggested exponential scaling with inverse temperature. We prove that fast mixing at low temperature for the two-dimensional toric code can be achieved by augmenting local jump operators with simple global jump operators, which enable efficient transitions between logical sectors. To establish tight lower bounds on the spectral gap, we introduce a new reduction method that eventually maps the problem to estimating the spectral gap of a perturbed graph Laplacian on a stair graph. Our proof also shows that the Lindblad dynamics with a digitally implemented low-temperature local Davies generator is able to efficiently drive the quantum state toward the ground state manifold.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Realization of two-dimensional discrete time crystals with anisotropic Heisenberg coupling

A discrete time crystal (DTC) is an out-of-equilibrium phase of matter that spontaneously breaks discrete time-translation symmetry. Previous studies have been limited to a set of models with Ising-like couplings - and mostly only in one dimension - thus precluding our understanding of the existence (or not) of DTCs in models with more realistic interactions. In this work, by combining the latest generation of IBM quantum processors with state-of-the-art tensor network methods, we demonstrate the existence of a DTC in a two-dimensional system governed by anisotropic Heisenberg interactions. We uncover a rich phase diagram encompassing spin-glass, ergodic, and time-crystalline phases, and identify the interplay of initialization, interaction anisotropy, and driving protocols in stabilizing the DTC phase. By extending the study of Floquet matter beyond simplified models, we lay the groundwork for exploring how driven systems bridge the gap between quantum coherence and emergent non-equilibrium thermodynamics.

Phase transitions and critical phenomena↗
DOE OSTI2025

Keldysh tuning of photoluminescence in a lead halide perovskite crystal

In 1964, Keldysh laid the groundwork for strong-field physics in atomic, molecular, and solid-state systems by delineating a ubiquitous transition from multiphoton absorption to quantum electron tunneling under intense AC driving forces. While both processes in semiconductors can generate carriers and result in photon emission through electron-hole recombination, the low quantum yields in most materials have hindered direct observation of the Keldysh crossover. Leveraging the large quantum yields of photoluminescence in lead halide perovskites, we show that we can not only induce bright light emission from extreme sub-bandgap light excitation but also distinguish between photon-induced and electric-field-induced processes. Our results are rationalized by the Landau-Dykhne formalism, providing insights into the non-equilibrium dynamics of strong-field light-matter interactions. These findings open new avenues for light upconversion and sub-bandgap photon detection, highlighting the potential of lead halide perovskites in advanced optoelectronic applications.

FOS: Physical sciences↗
DOE OSTI2026

Field-tailoring quantum materials via magneto-synthesis: metastable metallic and magnetically suppressed phases in a trimer iridate

We demonstrate that applying modest magnetic fields (< 0.1 T) during high-temperature crystal growth can profoundly alter the structure and ground state of a spin-orbit-coupled, antiferromagnetic trimer lattice. Using BaIrO₃ as a model system, whose ground state is intricately dictated by the trimer lattice, we show that magneto-synthesis , a field-assisted synthesis approach, stabilizes a structurally compressed, metastable metallic and magnetically suppressed phases inaccessible via conventional methods. These effects include a 0.85% reduction in unit cell, 4-order-of-magnitude decrease in resistivity, a 10-fold enhancement of the Sommerfeld coefficient, and the collapse of long-range magnetic order -- all intrinsic and bulk in origin. First-principles calculations confirm that the field-stabilized structure lies substantially above the ground state in energy, highlighting its metastable character. These large, coherent and correlated changes across multiple bulk properties, unlike those caused by dilute impurities, defects or off-stoichiometry, point to an intrinsic field-induced mechanism. The findings establish magneto-synthesis as a powerful new pathway for accessing non-equilibrium quantum phases in strongly correlated materials.

magneto-synthesis↗
NASA NTRS1995

RTE: A Computer Code for Rocket Thermal Evaluation

The numerical model for a rocket thermal analysis code (RTE) is discussed. RTE is a comprehensive thermal analysis code for thermal analysis of regeneratively cooled rocket engines. The input to the code consists of the composition of fuel/oxidant mixture and flow rates, chamber pressure, coolant temperature and pressure. dimensions of the engine, materials and the number of nodes in different parts of the engine. The code allows for temperature variation in axial, radial and circumferential directions. By implementing an iterative scheme, it provides nodal temperature distribution, rates of heat transfer, hot gas and coolant thermal and transport properties. The fuel/oxidant mixture ratio can be varied along the thrust chamber. This feature allows the user to incorporate a non-equilibrium model or an energy release model for the hot-gas-side. The user has the option of bypassing the hot-gas-side calculations and directly inputting the gas-side fluxes. This feature is used to link RTE to a boundary layer module for the hot-gas-side heat flux calculations.

Mohammad H N Naraghi↗
DOE OSTI2025

Lepto-axiogenesis with light right-handed neutrinos

We study lepto-axiogenesis in theories where the right-handed neutrino is light enough that its dynamics affect the determination of the baryon asymmetry. When compared with theories of high-scale lepto-axiogenesis where the Majorana neutrino mass may be treated as an effective dimension-five operator, we find that the predicted saxion mass is lower. Two distinct scenarios emerge. In the first, processes that generate the baryon asymmetry are in equilibrium down to the mass of the right-handed neutrino. In the second, the relevant processes never reach equilibrium; the baryon number freezes in. We comment on implications for supersymmetric spectra and discuss constraints on late decays of supersymmetric relics and from dark radiation. In contrast to high-scale lepto-axiogenesis, which predicts superpartners with masses of 10–100 TeV or more, we find this scenario is consistent with a wider range of superpartner masses, all the way down to current direct search bounds.

Axions and ALPs↗
DOE OSTI2026

Phonon-mediated ultrafast dynamics in self-assembled monolayers of 4-mercaptobenzoic acid on gold

Non-equilibrium interactions between plasmonic metals and adsorbed molecules lie at the heart of emerging applications such as plasmonic photocatalysis and sensing, though the ultrafast charge and energy transfer mechanisms arising from these interactions are not well understood. Herein, we investigate the ultrafast dynamics of Au nano-islands tethered with a self-assembled monolayer (SAM) of electron-withdrawing 4-mercaptobenzoic acid (4MBA) molecules. Ultrafast UV-visible transient absorption spectroscopy following excitation of the interband transition in Au reveals three well-known, characteristic time constants that quantify electron–electron (el–el), electron–phonon (el–ph) and phonon–phonon (ph–ph) scattering lifetimes. When comparing the dynamics of bare Au and 4MBA-Au, we find that the el–ph and ph–ph scattering lifetimes are notably longer in 4MBA-Au. Density functional perturbation theory calculations ascribe the elongation in el–ph lifetimes in 4MBA-Au to the significant coupling of acoustic phonon modes of Au with certain molecular vibrations of 4MBA, leading to decreased spatial overlap between carrier electronic states and the acoustic modes. We speculate that the elongation of ph–ph scattering lifetimes in 4MBA-Au arises due to poor thermal conductivity of the SAM which disrupts efficient energy dissipation from Au to the environment, thus slowing down the thermalization of phonons. This work provides a glimpse into how molecular adsorbates modify the charge carrier and phonon dynamics of Au and sets the stage for further systematic exploration of plasmonic metal–molecule interactions.

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

The impact of metastability on the high-pressure behavior of cerium

The structures adopted by solids under pressure are often assumed to reflect thermodynamic equilibrium, yet in many materials phase selection is strongly influenced by kinetic pathways and microstructural inheritance. Elemental cerium (Ce) exemplifies this challenge, with decades of conflicting reports describing different high-pressure crystal structures emerging under nominally identical conditions. Here we use neutron diffraction from large ( ~ 60 mm 3 ) sample volumes to follow the structural evolution in ultra-high-purity Ce during controlled pressure-temperature cycling between 85 and 295 K and up to 8 GPa. We find that the crystal structure formed at high pressure depends on the compression pathway: slow compression ( ~ 0.25 GPa hr −1 ) at room temperature favors an orthorhombic phase (α'), whereas slow ( ~ 0.25 GPa hr −1 ) and also moderately faster ( ~ 0.5 GPa hr −1 ) compression at low temperature stabilizes a pure monoclinic phase (α"). The low-temperature phase persists metastably over a wide temperature range but transforms irreversibly upon heating above ~ 280 K, or modest pressure cycling. Remarkably, the lower-pressure γ phase remains trapped far beyond the expected stability field, persisting to the highest pressures studied. These observations show that phase selection in Ce is governed by kinetics and microstructural memory rather than equilibrium thermodynamics or sample morphology alone, establishing path dependence as a defining feature of its high-pressure behavior.

Ridley, Christopher J. [Oak Ridge National Laborat↗
DOE OSTI2026

Ultrafast Raman probe of the photoinduced superconducting to normal state transition in the cuprate Bi 2 ⁢Sr 2 ⁢CaCu 2 ⁢O 8+𝛿

Here, we report an ultrafast Time-Resolved Raman scattering study of the out-of-equilibrium photoinduced dynamics across the superconducting to normal state phase transition of the cuprate Bi 2 ⁢Sr 2⁢ CaCu 2 ⁢O 8+𝛿 . Using the polarization-resolved momentum space selectivity of Raman scattering, we track the superconducting condensate destruction and recovery dynamics with subpicoseconds time resolution in the antinodal region of the Fermi surface where the superconducting gap is maximum. Leveraging ultrafast Raman thermometry, we find a significant dichotomy between the superconducting condensate and the quasiparticle temperature dynamics near the antinodes, which cannot be framed in terms of a single effective electron temperature. The present work demonstrates the ability of Time-Resolved Raman scattering to selectively probe out-of-equilibrium pathways of different electronic subdegrees of freedom during a photoinduced phase transition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRSDate not supplied

A Machine Learning Framework for Error Compensation in Radiative Transfer Calculations

Radiative heat transfer influences the amount of heat flux transferred to the surface of the hypersonic vehicle, which is essential to evaluate the performance of thermal protection systems. The radiative heat flux is found to be computationally prohibitive while accounting for the variation in spatial, angular, and spectral domains. A new methodology has been recently developed to alleviate the cost of computation in the spectral domain by constructing flow-agnostic reduced-order models (ROMs). The developed spectral ROM databases provide grouping strategies that account for non-equilibrium absorption and emission as well as interaction between disparate species due to spectral overlap in associated radiative processes. However, the developed ROMs need to be optimized for a specific combination of interacting gas species and would need to re-calibrated in case individual species are added/omitted. In this work, we use various machine learning (ML) techniques to approximate the radiative intensities determined by a ROM optimized for a specific gas mixture. The ML model relies on the ROM databases developed for a single species which ignores any spectral overlap. Thus, radiation evaluation starts with a simple summation of radiative intensities predicted using these non-calibrated ROMs for the contributing species. The ML framework then provides a correction to account for the interplay in the frequency, i.e., emission of photons by one species and absorption by another, and yields mixture-specific radiation fields. Once trained on the individual ROM databases, the ML framework offers instantaneous corrections that serves as a time/cost effective alternative to the optimization of ROMs for a specific gas mixture. The ML framework is trained on both the high fidelity and ROM evaluated line of sight (LOS) data from Orion, Stardust, and FIRE II cases to obtain a general purpose correction model for earth re-entry scenarios when radiation contributions from both atomic nitrogen and atomic oxygen are considered. A geometric length scale parameter is used in the training process to account for errors introduced in the ROM databases as a consequence of high optical thickness. The efficacy of the ML framework is underscored through extensive analysis of train and test errors with respect to all the re-entry scenarios. The applicability of such an ML framework was further corroborated by embedding it in a state-of-the-art US3D - NERO system for determining the radiative heat flux transferred to the hypersonic vehicle surface.

Radiation↗
NASA NTRS2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λc, that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λc-based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

Amanda Stricklan↗
NASA NTRS2000

Chalcogenide-Based Non-Volatile Memory Technology

Chalcogenide is a proven phase change material used in re-writeable CDs and DVDs. This material changes phases, reversibly and quickly, between an amorphous state that is dull in appearance and electrically high in resistance, and a polycrystalline state that is highly reflective and low in resistance. The application of this commercially proven technology to create dense, high-speed, non-volatile semiconductor memories is discussed.

J Maimon↗
DOE OSTI2026

Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling

Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .

Close-clearance impellers↗
DOE OSTI2025

Entanglement Structure of Non-Gaussian States and How to Measure It

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.

Fermi gases↗
DOE OSTI2026

Sparse non-Markovian Noise Modeling of Transmon-Based Multi-Qubit Operations

The influence of noise on quantum dynamics is one of the main factors preventing current quantum processors from performing accurate quantum computations. Sufficient noise characterization and modeling can provide key insights into the effect of noise on quantum algorithms and inform the design of targeted error protection protocols. However, constructing effective noise models that are sparse in model parameters, yet predictive can be challenging. In this work, we present an approach for effective noise modeling of multi-qubit operations on transmon-based devices. Through a comprehensive characterization of seven devices offered by the IBM Quantum Platform, we show that the model can capture and predict a wide range of single- and two-qubit behaviors, including non-Markovian effects resulting from spatiotemporally correlated noise sources. The model’s predictive power is further highlighted through multi-qubit dynamical decoupling demonstrations and an implementation of the variational quantum eigensolver. As a training proxy for the hardware, we show that the model can predict expectation values within a relative error of 0.5%; this is a sevenfold improvement over default hardware noise models. Through these demonstrations, we highlight key error sources in superconducting qubits and illustrate the utility of reduced noise models for predicting hardware dynamics.

open quantum systems & decoherence↗