Search NASA⌕ Search

SEARCH · Search NASA

Results for “THERMODYNAMIC PROPERTY”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Probabilistic and maximum entropy modeling of chemical reaction systems: Characteristics and comparisons to mass action kinetic models

We demonstrate and characterize a first-principles approach to modeling the mass action dynamics of metabolism. Starting from a basic definition of entropy expressed as a multinomial probability density using Boltzmann probabilities with standard chemical potentials, we derive and compare the free energy dissipation and the entropy production rates. We express the relation between entropy production and the chemical master equation for modeling metabolism, which unifies chemical kinetics and chemical thermodynamics. Because prediction uncertainty with respect to parameter variability is frequently a concern with mass action models utilizing rate constants, we compare and contrast the maximum entropy model, which has its own set of rate parameters, to a population of standard mass action models in which the rate constants are randomly chosen. We show that a maximum entropy model is characterized by a high probability of free energy dissipation rate and likewise entropy production rate, relative to other models. We then characterize the variability of the maximum entropy model predictions with respect to uncertainties in parameters (standard free energies of formation) and with respect to ionic strengths typically found in a cell.

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

Equation of state measurement of detonation carbon condensates using optical microscopy and interferometry

Thermochemical models of detonation that estimate performance (e.g., detonation velocity, energy delivery, etc.), are based on assumptions that carbon condensates (soot) formed during detonation is largely similar to bulk carbon. However, soot constituents can range from amorphous carbon to nanodiamond and include other material phases. Since thermodynamic properties of the soot such as compressibility are imperative for accurate thermochemical modeling of detonation reaction chemistry, experimental measurements of the equation of state (EOS) which determine the compressibility are vital. Due to the mixed-phase nature of detonation soot, typical methods to measure the EOS (e.g., x-ray diffraction) are untenable. In this study, the high-pressure EOS up to 20 GPa was determined for detonation soot collected from PBX 9502, Composition B (Comp B), Hexanitrostilbene (HNS), and LX-21 high explosives by employing a direct volume technique using optical microscopy and interferometry in a diamond anvil cell. Comp B soot was determined to be the least compressible [K 0 = 57.9(17) GPa] with HNS soot [K 0 = 53.7(15) GPa], LX-21 soot [K 0 = 45.8(51) GPa], and PBX 9502 soot [K 0 = 28.2(27) GPa] being more compressible, likely due to differences in nanodiamond content as compared to amorphous carbon and graphite content.

Amorphous materials↗
DOE OSTI2025

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗
DOE OSTI2025

Simultaneous inference of equation of state parameters and unknown data errors with uncertainty quantification via hierarchical Bayesian posterior maximization

Equations of state (EOSs) are a key component in running hydrodynamic simulations as they relate the thermodynamic states for the material. The Davis reactants EOS is commonly used for modeling high explosives (HEs), and the EOS model parameters are calibrated using material specific data. The calibrations are often performed with uncertainty quantification via Bayesian inference to account for uncertainty in the data and generate ensembles of likely parameters. However, there are relatively few HE data sets to use for calibration and many are historical and lack error information. In this work, we simultaneously calibrate the Davis reactants EOS model parameters and unknown data error terms for the high explosive PBX 9501. To quantify the uncertainty in the models and the data, we use a Bayesian framework for the calibration and compute the hierarchical Bayesian posterior distribution with both a posteriori maximization approach and Markov Chain Monte Carlo. In general, we find that, given our assumptions, the two approaches result in similar calibrated parameters, posterior covariance matrices, and insights about the parameters but that the posterior maximization requires far less computational resources.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Modeling the behavior of concentrated aqueous HNO 3 using machine learning interatomic potentials

We develop two multi-defect machine learning interatomic potentials (MLIPs) trained at the BLYP-D2 and PBE-D3 density functional theories using the DeepMD-kit, allowing for the investigation of structural and thermodynamic properties of nitric acid over a wide range of concentrations via molecular dynamics (MD) simulations. We directly compute the degree of dissociation, α, and pK a from MD simulations, revealing that HNO 3 behaves as a weaker acid at higher concentrations, noting that our standard-state pK a value is in excellent agreement with the experimental one. In general, good agreement is observed with experimental results such as α and density outside the training dataset, with only modest deviations at low-to-medium concentrations. We benchmark our custom multi-defect DeepMD MLIPs against foundational models MACE-MP0 and MACE-OFF23. The foundation models capture some aspects of HNO 3 /NO 3 − solvation in concentrated nitric acid but show noticeable density errors and miss subtle structural features relevant to spectroscopy, whereas the bespoke DeepMD MLIPs yield more compact solvation shells, reproduce density-concentration trends, and run ∼12–15× faster than MACE-MP0. Although classical FFs are still more efficient and match experimental densities better, they lack chemical reactivity and thus cannot predict α or pK a , underscoring the need for system-specific reactive MLIPs beyond universal MLIPs.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2025

Disorder-induced spin excitation continuum and spin-glass ground state in the inverse spinel CuGa 2 O 4

Spinel-structured compounds serve as prototypical examples of highly frustrated systems and are promising candidates for realizing the long-sought quantum spin liquid (QSL) state. However, structural disorder is inevitable in many real QSL candidates and its impact remains a topic of intense debate. In this work, we conduct comprehensive investigations on CuGa 2 ⁢ O 4 , a spinel compound with significant structural disorder, focusing on its thermodynamic properties and spectroscopic behaviors. No long-range magnetic order is observed down to ∼80 mK, as evidenced by magnetic susceptibility, specific-heat, and elastic neutron scattering measurements. More intriguingly, inelastic neutron scattering experiments reveal a broad gapless continuum of magnetic excitations around the Brillouin zone boundary, resembling the magnetic excitation spectra expected for a QSL. Nevertheless, a spin-freezing transition at 𝑇 f ≈ 0.88 K is identified from the cusp in the dc susceptibility curves, where a bifurcation between zero-field-cooling and field-cooling curves occurs. Furthermore, ac susceptibility measurements show a peak close to 𝑇 f at low frequency, which shifts to higher temperature with increasing frequency. These results show that CuGa 2 ⁢ O 4 has a spin-glass ground state, consistent with the establishment of short-range order inferred from the specific-heat measurements. Collectively, these results illustrate the crucial role of disorder in defining the excitation spectrum out of the disordered ground state. Furthermore, our findings shed light onto the broader class of 𝐴⁢𝐵 2 ⁢O 4 spinels and advance our understanding of the spin dynamics in magnetically disordered systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Vortex motion induced losses in tantalum resonators

Tantalum (Ta)-based superconducting circuits have been demonstrated to enable ultrahigh qubit quality factors (𝑄), motivating a careful study of the microscopic origin of the remaining losses that limit their performance. We have recently shown that the losses in Ta-based resonators are dominated by two-level systems at low microwave powers and millikelvin temperatures. We also observe that some devices exhibit loss that is exponentially activated at a lower temperature inconsistent with the superconducting critical temperature (T 𝑐 ) of the constituent film. Specifically, dc resistivity measurements show a T 𝑐 of over 4 K, while microwave measurements of resonators fabricated from these films show losses that increase exponentially with temperature with an activation energy as low as 0.3 K. Here, in this study, we present a comparative study of the structural and thermodynamic properties of Ta-based resonators and identify vortex motion-induced loss as the source of thermally activated microwave loss. Through careful magnetoresistance and x-ray diffraction measurements, we observe that the increased loss occurs for films that are in the clean limit, where the superconducting coherence length (𝜉) is shorter than the mean free path (𝑙). Vortex motion-induced losses are suppressed for films in the dirty limit, which show evidence of structural defects that can pin vortices. We verify this hypothesis by explicitly pinning vortices via patterning, and we find that we can suppress the loss by microfabrication.

36 MATERIALS SCIENCE↗
DOE OSTI2025

String-based model with Hagedorn temperature of 𝑇 𝐻 ∼ 300 MeV describes the spectrum of mesons and glueballs

We consider the thermodynamics of a color-confined phase of quantum chromodynamics (QCD) and pure gauge theory within a string-inspired corresponding to a physical spatial dimension, 𝑑 = 3. We show that the physical mass spectrum of massive mesons—in both the strange and nonstrange sectors separately—is reasonably well described and extended by the exponential mass spectrum of open strings, 𝜌⁡(𝑚), characterized by a unique Hagedorn temperature 𝑇 𝐻 = $\sqrt3⁢𝜎/2⁢𝜋$, expressed by the string tension 𝜎. This 𝑇 𝐻 is the value appropriate for 𝑑 = 3 spatial dimensions and is of the order of 𝑇 𝐻 ∼ 300 MeV for typical values of the string tension. It is much larger than the values of 𝑇 𝐻 which have been phenomenologically extracted so far to describe the meson spectrum. Glueball states in pure gauge theory, modeled by closed strings, exhibit a similarly large Hagedorn temperature, highlighting a universal feature of the exponential spectrum. We further analyze the thermodynamic properties of the equation of state at finite temperature and demonstrate that, in the confined phase, the string models agree with lattice QCD results. This lends further support to the recent interpretation of the QCD phase diagram that incorporates strings as relevant degrees of freedom.

Phase transitions↗
DOE OSTI2026

Direct estimation of the density of states for fermionic systems

Simulating time evolution is one of the most natural applications of quantum computers and is thus one of the most promising prospects for achieving practical quantum advantage. Here, we develop quantum algorithms to extract thermodynamic properties by estimating the density of states (DOS), which is a central object in quantum statistical mechanics. We introduce several key innovations that significantly improve the practicality and extend the generality of previous techniques. First, our approach allows one to estimate the DOS only for a specific subspace of the full Hilbert space. This is crucial for fermionic systems, since both canonical and grand canonical ensemble thermal equilibrium properties depend on subspaces of fixed number. Second, in our approach, by time evolving very simple, random initial states, such as randomly chosen computational basis states, we can exactly recover the DOS on average. Third, due to circuit-depth limitations, we only reconstruct the DOS up to a convolution with a Gaussian window—thus all imperfections that shift the energy levels by less than the width of the convolution window will not significantly affect the estimated DOS. For these reasons, we find the approach is a promising candidate for early quantum advantage as even short-time, noisy dynamics can yield a semiquantitative reconstruction of the DOS (convolution with a broad Gaussian window), while early fault-tolerant devices will likely enable higher-resolution DOS reconstruction through longer time evolutions. We demonstrate the practicality of our approach in representative Fermi-Hubbard and spin models and indeed find that our approach is highly robust against algorithmic errors in the time evolution and against gate noise. We further demonstrate that our approach is compatible with noisy intermediate-scale quantum (NISQ) computing NISQ-friendly variational techniques, introducing and leveraging a technique for variational time evolution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Manipulating Cu vacancy in Yb triangular lattice Materials: Investigating the quantum disordered compound CuYbSe 2 and spin glass Cu 2.73 ⁢YbP 2

Yb-based triangular lattice materials have attracted significant attention due to their potential for hosting exotic magnetic quantum states. In this paper, we report the structural, magnetic, and thermodynamic properties of two Yb triangular lattice compounds with vacancy disorder at nonmagnetic Cu sites: CuYbSe 2 and Cu 2.73 ⁢YbP 2 . Single-crystal x-ray diffraction determines notable Cu-site vacancies in both materials. Unlike a single Cu site with 0.5 occupancy in CuYbSe 2 , Cu 2.73 ⁢YbP 2 only presents 27% vacancies on one of three Cu sites. Magnetic susceptibility measurements indicate dominant antiferromagnetic interactions between Yb 3+ ions. The effective magnetic moment is consistent with an effective spin-1/2 state of Yb 3+ . CuYbSe 2 exhibits strong magnetic anisotropy between in-plane and out-of-plane directions with the estimated exchange interactions 𝐽 𝑎 /𝑘 B = 6.48 K and 𝐽 𝑐 /𝑘 B = 2.55 K. No magnetic ordering is observed down to 0.4 K, supporting its quantum disordered ground state and candidacy as a quantum spin liquid. In contrast, Cu 2.73 ⁢YbP 2 shows clear evidence of spin freezing at 4.5 K in both magnetization and heat capacity data, indicating a spin glass state. The field-induced magnetic orders are found in CuYbSe 2 . These results highlight the critical role of structural vacancies in perturbing Yb-based triangular lattices. In conclusion, we demonstrate that while the intrinsic Yb triangular lattices are comparable, the presence of Cu vacancies can fundamentally alter the magnetic ground state—tuning the system between quantum disordered and spin glass behavior.

Frustrated magnetism↗
DOE OSTI2024

Mechanisms and models of the turbulent boundary layers at transcritical conditions

Computational models for high-pressure transcritical turbulence in wall-modeled large-eddy simulation typically rely on wall-function models to overcome the need for resolving the boundary layer structures. However, the mechanisms and models of turbulent boundary layer at transcritical conditions remain poorly understood since the near-wall flow and heat transfer are significantly affected by the enhanced fluctuations and steep gradients in thermodynamic properties. Here, to address this issue, we study the mechanisms of transcritical turbulent boundary layers and wall-attached models at transcritical conditions. It is shown that the real-fluid variable-property effects are associated with the coupling between the near-wall cycle and the wall-normal coherent motions in the outer layer, resulting in the amplification of turbulent energy in the outer layer and noticeable energy transfer between the log-layer and the outer layer; hence, turbulence in the log-layer is modulated by the outer layer. Based on this underlying physical principle, we propose the characteristic velocity and length scales for the attached eddy at transcritical conditions, and extend the attached eddy model to transcritical turbulent boundary layers by introducing a mixed scaling that incorporates both inner and outer scalings. We show that the new characteristic scales and extended attached-eddy model perform well in characterizing the structures in transcritical turbulent boundary layers.

Li, Fangbo↗
DOE OSTI2024

Thermal activation signatures of the Anderson insulator and the Wigner solid forming near ν = 1

When interactions overcome disorder, integer quantum Hall plateaus support topological phases with different bulk insulators. In the center of the ν = 1 plateau the bulk is an Anderson-type insulator, while in the flanks of the plateau the bulk is the integer quantum Hall Wigner solid. We find that the activation energy along the ν = 1 plateau exhibits a very dramatic nonmonotonic dependence on the magnetic field, a dependence that is strongly correlated with the stability regions of the two phases. Furthermore, the activation energy has an unexpected minimum at the boundary between the Anderson insulator and the Wigner solid. Our findings constrain the theory of the integer quantum Hall Wigner solid, determine its thermodynamic properties, and reveal unusual behavior at the boundary between the Anderson insulator and the Wigner solid. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamical Correlations and Order in Magic-Angle Twisted Bilayer Graphene

The interplay of dynamical correlations and electronic ordering is pivotal in shaping phase diagrams of correlated quantum materials. In magic-angle twisted bilayer graphene, transport, thermodynamic, and spectroscopic experiments pinpoint a competition between distinct low-energy states with and without electronic order, as well as between localized and delocalized charge carriers. In this study, we utilize dynamical mean-field theory on the topological heavy fermion model of twisted bilayer graphene to investigate the emergence of electronic correlations and long-range order in the absence of strain. We contrast moment formation, Kondo screening, and ordering on a temperature basis and explain the nature of emergent correlated states based on three central phenomena: (i) the formation of local spin and valley isospin moments around 100 K, (ii) the ordering of the local isospin moments around 10 K preempting Kondo screening, and (iii) a cascadic redistribution of charge between localized and delocalized electronic states upon doping. At integer fillings, we find that low-energy spectral weight is depleted in the symmetric phase, while we find insulating states with gaps enhanced by exchange coupling in the zero-strain ordered phases. Doping away from integer filling results in distinct metallic states: a “bad metal” above the ordering temperature, where scattering off the disordered local moments suppresses electronic coherence, and a “good metal” in the ordered states with coherence of quasiparticles facilitated by isospin order. This finding reveals coherence from order as the microscopic mechanism behind the Pomeranchuk effect observed experimentally by Rozen [] and by Saito []. Upon doping, there is a periodic charge reshuffling between localized and delocalized electronic orbitals leading to cascades of doping-induced Lifshitz transitions, local spectral weight redistributions, and periodic variations of the electronic compressibility ranging from nearly incompressible to negative. Our findings highlight the essential role of charge transfer, hybridization, and ordering in shaping the electronic excitations and thermodynamic properties in twisted bilayer graphene and provide a unified understanding of the most puzzling aspects of scanning tunneling spectroscopy, transport, and compressibility experiments. Published by the American Physical Society 2024

Physics↗
DOE OSTI2025

Kerr effective black hole geometries in supergravity

We derive the explicit embedding of the effective Kerr spacetimes, which are pertinent to the vanishing of static Love numbers, soft hair descriptions of Kerr black holes, and low-frequency scalar-Kerr scattering amplitudes, as solutions within 𝑁 = 2 supergravity. These spacetimes exhibit a hidden 𝑆⁢𝐿⁡(2,𝑅) × 𝑈⁡(1) or 𝑆⁢𝑂⁡(4,2) symmetry resembling the so called subtracted geometries with 𝑆⁢𝐿⁡(2,𝑅) × 𝑆⁢𝐿⁡(2,𝑅) symmetry, which accurately represent the near-horizon geometry of Kerr black holes and, as we will argue most accurately represents the internal structure of the Kerr black hole. To quantify the differences among the effective Kerr spacetimes, we compare their physical quantities, internal structures, and geodesic equations. Although their thermodynamic properties, including entropy, match those of Kerr, our study uncovers significant differences in the interiors of these effective Kerr solutions. A careful examination of the internal structure of the spacetimes highlights the distinctions between various effective Kerr geometries and their quasinormal spectra.

quantum aspects of black holes↗
DOE OSTI2024

A Perspective on Scalable AI on High-Performance Computing and Leadership Class Supercomputing Facilities [Industrial and Governmental Activities]

Many scientific applications that support the mission of the US Department of Energy (US-DoE) require modeling complex engineering and/or physical systems. Here, examples of such complex systems arise from: (a) materials science to develop new compounds with exceptional mechanical and thermodynamical properties (e.g., resistance to mechanical stresses and high temperatures), (b) structural and nuclear engineering to model the temporal evolution of the structural damage of concrete shields exposed to continuous neutron and gamma radiations emitted by the nuclear reactor core, (c) urban sciences (e.g., transportation and smart buildings), and (d) power grid systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Effects of Heat Treatment on Microstructure, Thermal Transport and Stability in SRF Niobium

Efficient heat transfer and thermal stability are critical for preventing thermal quench and ensuring the high performance and reliability of superconducting radiofrequency (SRF) cavities, especially under high RF fields. In this article, we investigate thermal conductivity and heat capacity of high-purity cold-worked niobium used in SRF applications, with a particular focus on the effects of high-temperature heat treatments. Our measurements reveal a pronounced sensitivity of thermal transport and thermodynamic properties to the underlying microstructure such as grain size and dislocation density. The results show the critical role of controlled heat treatment processes in optimizing the thermal performance of SRF niobium, providing valuable insights for improving cavity fabrication and processing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thermodynamic Modeling of Point Defects in Triple Conducting Perovskite Ba 0.95 La 0.05 FeO 3- $_δ$ with Incorporation of the Hydride Defect Formation Reaction for Solid Oxide Cells

Distinct from the proton defect, the hydride defect species may be present in certain perovskite materials in reducing environments such as in fuel electrodes of solid oxide cells (SOCs) or on the reducing side of ceramic membranes. A generalized defect thermodynamic model was developed for the triple-conducting perovskites (La,Ba)Fe 1-x M x O 3-δ (M = Y and Zr) to allow inclusion of the hydride defect formation reaction in addition to the other three main defect reactions, namely, the oxygen vacancy formation, hydration, and charge disproportionation reactions. This comprehensive defect model also allows the incorporation of polynomial functional forms of oxygen nonstoichiometry δ to describe the defect reaction energies and entropies and to enable refinements of the defect reaction equilibrium constants in the defect thermodynamic analysis. As a first step, the developed model is applied to the Ba 0.95 La 0.05 FeO 3-δ material as an illustrative system to obtain its Brouwer diagrams with both the proton and hydride defects in relevant SOC conditions, particularly for more reducing environments. In conclusion, the results provide direct guidance on the influence of electronic and ionic defect concentrations upon thermodynamic properties and ultimately on the performance of Ba 0.95 La 0.05 FeO 3-δ and potentially other (La,Ba)Fe 1-x M x O 3-δ perovskite materials involved in SOC applications.

25 ENERGY STORAGE↗
DOE OSTI2026

Direct Measurement of Cl − Activity in Metal Chloride Molten Salts using a Cl 2 /Cl − Electrode

Molten metal chloride salts are promising candidates for advanced heat transfer fluids in generation IV nuclear reactors and beyond. The activity of the chloride ion in the salt has a large influence on the redox characteristics of the corresponding molten salt. There exists a knowledge gap for the direct measurement of chlorobasicity (Lewis basicity) in molten-salt mixtures. Here, the focus of this work was to develop a Cl 2 /Cl − electrode for the direct measurement of chloride activity in molten metal-chloride salt mixtures. The LiCl–KCl eutectic system was utilized as the reference melt to determine the change in chloride ion activity as increasing concentrations of MgCl 2 are added to the LiCl–KCl eutectic working electrode. The results of the measurements indicated a reduction in chloride ion activity as the concentration of magnesium chloride increased, consistent with the complexation of free chlorides by Lewis acidic magnesium cationic species. The Temkin model was used to estimate the thermodynamic properties of MgCl 4 2− complex.

Chloride↗