Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equation of 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.

At least 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Stochastic relativistic viscous hydrodynamics from the Metropolis algorithm

We propose an algorithm for simulating stochastic relativistic fluid dynamics based on Metropolis updates. Each step of the algorithm begins with an update based on ideal hydrodynamics. This is followed by proposing random (spatial) momentum transfers between fluid cells, keeping the total energy fixed. These proposals are then accepted or rejected using the change in entropy as a statistical weight. The algorithm reproduces relativistic viscous hydrodynamics in the “density frame,” which is a formulation of viscous hydrodynamics we review and clarify here. This formulation is first order in time and requires no auxiliary dynamical fields such as Π 𝜇⁢𝜈 . The only parameters are the shear and bulk viscosities and the equation of state. Here, by adopting the 3+1 split of general relativity, we extend the Metropolis algorithm to general space-time coordinates, such as Bjorken coordinates, which are commonly used to simulate heavy-ion collisions.

Hydrodynamic noise↗
DOE OSTI2025

Feasibility of measuring the speed of sound of the quark-gluon plasma from the multiplicity and mean 𝑝 𝑇 of ultracentral heavy-ion collisions

The mean transverse momentum ⟨𝑝 𝑇 ⟩ of hadrons has been observed experimentally and in numerical simulations to have a power-law dependence on the hadronic multiplicity 𝑁 in ultracentral relativistic heavy-ion collisions: ⟨𝑝 𝑇 ⟩∝𝑁 𝑏 UC . It has been put forward that this exponent 𝑏 UC is the speed of sound of quark-gluon plasma measured at a temperature determined from ⟨𝑝 𝑇 ⟩. We study step by step the connection between (i) the energy and entropy of hydrodynamic simulations and (ii) experimentally measurable observables. We show that an argument based on energy and entropy should yield an exponent equal to the pressure over energy density 𝑃/ɛ, rather than the speed of sound 𝑐$_s^2$; however, we also observe that ⟨𝑝 𝑇 ⟩ and 𝑁 are not sufficiently accurate proxies for the energy and entropy to make this possible in practice. From simulations, we find that the exponent 𝑏 UC is significantly different whether the “effective volume” is strictly constant or not, a condition that cannot be enforced experimentally. Additional tests using a modified equation of state find that the exponent 𝑏 UC exhibits a variable degree of correlations with the speed of sound and with 𝑃/ɛ, but is not an accurate measurement of either quantity in general.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Onset of Constituent Quark Number Scaling in Heavy-Ion Collisions at RHIC

Partonic collectivity is one of the necessary signatures for the formation of quark-gluon plasma in high-energy nuclear collisions. Number of constituent quarks (NCQ) scaling has been observed for hadron elliptic flow 𝑣2 in top energy nuclear collisions at the Relativistic Heavy Ion Collider and the LHC, and this has been theoretically suggested as strong evidence for partonic collectivity. In this Letter, a systematic analysis of 𝑣 2 of 𝜋 ± , 𝐾 ± , 𝐾$^{0}_{𝑆}$, 𝑝, and Λ in Au+Au collisions at $\sqrt{s_{NN}}$ = 3.2, 3.5, 3.9, and 4.5 GeV, with the STAR experiment at the Relativistic Heavy Ion Collider, is presented. NCQ scaling is markedly violated at 3.2 GeV, consistent with a hadronic-interaction dominated equation of state. However, as the collision energy increases, a gradual evolution to NCQ scaling is observed. In conclusion, this beam-energy dependence of 𝑣 2 for all hadrons studied provides evidence for the onset of dominant partonic interactions by $\sqrt{s_{NN}}$ = 4.5 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Exploring the role of 𝑑* hexaquarks on quark deconfinement and hybrid stars

We investigate the impact of the d*( 2380) hexaquark on the equation of state (EOS) of dense matter within hybrid stars (HSs) using the chiral mean-field model (CMF). The hexaquark is included as a new degree of freedom in the hadronic phase, and its influence on the deconfinement transition to quark matter is explored. We reparametrize the CMF model to ensure compatibility with recent astrophysical constraints, including the observation of massive pulsars and gravitational wave events. Our results show that the presence of d* significantly modifies the EOS, leading to a softening at high densities and a consequent reduction in the predicted maximum stellar masses. Furthermore, we examine the possibility of a first-order deconfinement phase transition within the context of the extended stability branch of slow stable HSs (SSHSs). We find that the presence of hexaquarks can delay the deconfinement phase transition and reduce the associated energy density gap, affecting the structure and stability of HSs. Our results suggest that, as the hexaquark appearance tends to destabilize stellar configurations, fine-tuning of model parameters is required to obtain both the presence of hexaquarks and quark deconfinement in these systems. In this scenario, the SSHS branch plays a crucial role in obtaining HSs with hexaquarks that satisfy current astrophysical constraints. Our work provides new insights into the role of exotic particles like d* in dense matter and the complex interplay between hadronic and quark degrees of freedom inside compact stellar objects.

Nuclear astrophysics↗
DOE OSTI2025

Neutron stars with exceptionally light QCD axions

We present a comprehensive study of axion condensed neutron stars that arise in models of an exceptionally light axion that couples to quantum chromodynamics (QCD). These axions solve the strong-charge-parity (𝐶⁢𝑃) problem, but have a mass-squared lighter than that due to QCD by a factor of 𝜖 < 1. Inside dense matter, the axion potential is altered, and much of the matter in neutron stars resides in the axion condensed phase where the strong-𝐶⁢𝑃 parameter 𝜃 = 𝜋 and 𝐶⁢𝑃 remains a good symmetry. In these regions, masses and interactions of nucleons are modified, in turn changing the equation of state (EOS), structure, and phenomenology of the neutron stars. We take the first steps toward the study of the EOS of neutron star matter at 𝜃 = 𝜋 within chiral effective field theory and use relativistic mean field theory to deduce the resulting changes to nuclear matter and the neutron star low-density EOS. We derive constraints on the exceptionally light axion parameter space based on observations of the thermal relaxation of accreting neutron stars, isolated neutron star cooling, and pulsar glitches, excluding the region up to 5 ×10 −7 ≲ 𝜖 ≲ 0.2 for 𝑚 𝑎 ≳ 2 × 10 −9 eV. We comment on potential changes to the neutron star mass-radius relationship, and discuss the possibility of novel, nuclear-density compact objects with 𝜃 = 𝜋 that are stabilized not by gravity but by the axion potential.

axions↗
DOE OSTI2026

Pressure-induced 𝐵⁢1 to 𝐵⁢2 phase transition in CeN studied by 𝑎⁢𝑏 𝑖⁢𝑛⁢𝑖⁢𝑡⁢𝑖⁢𝑜 correlation matrix renormalization theory calculations

We apply correlation matrix renormalization theory (CMRT) to cerium nitride (CeN) under pressure. For 𝐵⁢1 (NaCl-type) phase, CMRT gives an equation of state consistent with ambient pressure experiments. It produces electronic density-of-state (DOS) characterized by a sharp 4⁢𝑓 quasiparticle resonance peak pinned at the Fermi level and two subbands formed by strong hybridization between the localized Ce-4⁢𝑓 electrons and the itinerant Ce-5⁢𝑑 and N-2⁢𝑝 electrons below the Fermi level, consistent with x-ray photoemission spectroscopy experiments. Upon compression, CMRT predicts a first-order 𝐵⁢1 → 𝐵⁢2 (CsCl-type) transition with ∼11% volume collapse in agreement with experiments. Across the transition, the 4⁢𝑓 spectrum broadens, the 4⁢𝑓 orbital occupancy increases, and the hybridization with conduction states enhances, signaling a crossover from partially localized to more itinerant 4⁢𝑓 behavior. Furthermore, these features are in excellent agreement with experimental observations, demonstrating that CMRT provides a parameter-free description and prediction of correlation-driven structural and electronic transitions in rare-earth compounds.

Ab initio calculations↗
DOE OSTI2026

SPT-3G D1: CMB temperature and polarization power spectra and cosmology from 2019 and 2020 observations of the SPT-3G main field

Here, we present measurements of the temperature and E-mode polarization angular power spectra of the cosmic microwave background (CMB) from observations of 4% of the sky with SPT-3G, the current camera on the South Pole Telescope (SPT). The maps used in this analysis are the deepest used in a CMB TT / TE / EE analysis to date. The maps and resulting power spectra have been validated through blind and unblind tests. The measurements of the lensed EE and TE spectra are the most precise to date at ℓ = 1800 – 4000 and ℓ = 2200 – 4000 , respectively. Combining our TT / TE / EE spectra with previously published SPT-3G CMB lensing results, we find parameters for the standard Λ cold dark matter ( Λ CDM ) model consistent with Planck and ACT DR6 with comparable constraining power. We report a Hubble constant of H 0 = 66.66 ± 0.60 km s - 1 Mpc - 1 from SPT-3G alone, 6.2 σ away from local measurements from SH0ES. For the first time, combined ground-based ( SPT + ACT ) CMB primary and lensing data have reached Planck’s constraining power on some parameters, a milestone for CMB cosmology. The combination of these three CMB experiments yields the tightest CMB constraints to date, with H 0 = 67.19 ± 0.38 km s - 1 Mpc - 1 , and the amplitude of clustering σ 8 = 0.8137 ± 0.0037 . CMB data alone show no evidence for physics beyond Λ CDM , however, we observe a 2.8 σ difference in Λ CDM between CMB and baryon acoustic oscillation (BAO) results from DESI-DR2, which is relaxed in extended models. The combination of CMB and BAO yields 2 − 3 σ shifts from Λ CDM in the curvature of the Universe, the amplitude of CMB lensing, or the dark energy equation of state. It also drives mild preferences for models that address the Hubble tension through modified recombination or variations in the electron mass in a nonflat Universe. This work highlights the growing power of ground-based CMB experiments and lays a foundation for further cosmological analyses with SPT-3G.

Camphuis, Etienne [Sorbonne Université; CNRS; Inst↗
DOE OSTI2025

Interface Modes in Inspiralling Neutron Stars: A Gravitational-Wave Probe of First-Order Phase Transitions

At the extreme densities in neutron stars, a phase transition to deconfined quark matter is anticipated. Yet masses, radii, and tidal deformabilities offer only indirect measures of a first-order phase transition, requiring many detections to resolve or being ineffective observables if the discontinuity exists at lower densities. Here, we report on a smoking-gun gravitational-wave signature of a first-order transition: the resonant tidal excitation of an interface mode. Using relativistic perturbation theory with an equation-of-state family informed by chiral effective field theory, we show that such a resonance may be detectable with next-generation interferometers and possibly already with LIGO A+ for sufficiently loud events.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

SPT-3G D1: Axion early dark energy with CMB experiments and DESI observations

We present the most up-to-date constraints on axion early dark energy (AEDE) from cosmic microwave background (CMB) and baryon acoustic oscillation (BAO) measurements. In particular, we assess the impact of data from ground-based CMB experiments, the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT)—both with and without Planck—on constraints on AEDE. We also highlight the impact that BAO information from the Dark Energy Spectroscopic Instrument (DESI) has on these constraints. From CMB data alone, we do not find statistically significant evidence for the presence of AEDE, and we find only moderate reduction in the Hubble tension. From the latest SPT data alone, we find the maximal fractional contribution of AEDE to the cosmic energy budget is f EDE < 0.12 at 95% confidence level (CL), and the Hubble tension between the SPT and Supernovae, H0, for the Equation of State of Dark Energy (SH0ES) results is reduced to the 2.3 σ level. When combining the latest SPT, ACT, and Planck datasets, we find f EDE < 0.070 at 95% CL and the Hubble tension at the 3.6 σ level. In contrast, adding DESI data to the CMB datasets results in mild preference for AEDE and, in some cases, non-negligible reduction in the Hubble tension. From SPT + DESI , we find f EDE = 0.08 1 - 0.052 + 0.037 at 68% CL, and the Hubble tension reduces to 1.5 σ . From the combination of DESI with all three CMB experiments, we get f EDE = 0.05 5 - 0.047 + 0.024 at 68% CL and a weak preference for AEDE over Λ CDM . This data combination, in turn, reduces the Hubble tension to 2.6 σ . We highlight that this shift in parameters when adding the DESI dataset is a manifestation of the discrepancy currently present between DESI and CMB experiments in the concordance model Λ CDM .

Khalife, A. R. [Paris, Inst. Astrophys.] (ORCID:00↗
DOE OSTI2025

Emulators for Scarce and Noisy Data: Application to Auxiliary-Field Diffusion Monte Carlo for Neutron Matter

Understanding the equation of state (EOS) of pure neutron matter is necessary for interpreting multimessenger observations of neutron stars. Reliable data analyses of these observations require well-quantified uncertainties for the EOS input, ideally propagating uncertainties from nuclear interactions directly to the EOS. This, however, requires calculations of the EOS for a prohibitively larger number of nuclear Hamiltonians, solving the nuclear many-body problem for each one. Quantum Monte Carlo methods, such as auxiliary-field diffusion Monte Carlo (AFDMC), provide precise and accurate results for the neutron matter EOS, but they are very computationally expensive, making them unsuitable for the fast evaluations necessary for uncertainty propagation. Here, we employ parametric matrix models to develop fast emulators for AFDMC calculations of neutron matter and use them to directly propagate uncertainties of coupling constants in the Hamiltonian to the EOS. As these uncertainties include estimates of the effective field theory truncation uncertainty, this approach provides robust uncertainty estimates for use in astrophysical data analyses. In conclusion, this Letter will enable novel applications such as using astrophysical observations to put constraints on coupling constants for nuclear interactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Symmetry-Energy Dependence of the Bulk Viscosity of Nuclear Matter

We clarify how the weak-interaction-driven bulk viscosity ζ and the bulk relaxation time τ Π of neutrino-transparent n p e matter depend on the nuclear symmetry energy. We show that, at saturation density, the equation-of-state dependence of these transport quantities is fully determined by the experimentally constrained nuclear symmetry energy S and its slope L . Variations of L can change the bulk viscosity by orders of magnitude, which can affect both the dissipative and the conservative tidal response of neutron stars. Furthermore, this suggests that both conservative and dissipative effects encoded in the gravitational-wave signatures of binary neutron star inspirals may help constrain nuclear symmetry-energy properties.

Neutron stars & pulsars↗
DOE OSTI2025

Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂⁡(10 2 ) ⁢GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂⁡(10 0 ) ⁢GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Monte Carlo methods↗
DOE OSTI2026

Machine learning-based interatomic potential development and phase transition analysis of ferroelectric hafnium dioxide

The ferroelectric phase (𝑃⁢𝑐⁢𝑎⁢2 1 , which is in orthorhombic symmetry) of hafnium dioxide (HfO 2 ) has gained much attention due to its potential applications in nanoelectronics and advanced memory devices. However, its complex phase behavior under external stimuli, such as pressure and temperature, remains a subject of intense investigation. This study focuses on developing a machine learning-based interatomic potential (MLIP) that is trained with data from density-functional theory (DFT) calculations to simulate phase transitions and mechanical properties of HfO 2 . The developed MLIP predicts lattice parameters, equations of state, bulk and shear moduli, and elastic constants that closely align with DFT predictions for several phases and at various pressures. Once validated, the MLIP is used to investigate the phase transitions of ferroelectric HfO 2 (𝑃⁢𝑐⁢𝑎⁢2 1 ) under both isobaric and constant stress conditions at elevated temperatures ranging from 200 to 2500 K. We used several complementary methods, including local symmetry identification, radial distribution function, and x-ray diffraction characterization, to identify interesting phase transitions among several competitive hafnia phases predicted from our simulations. The suggested methods uniformly reveal that under pure deviatoric condition, the system favors a transition from the orthorhombic 𝑃⁢𝑐⁢𝑎⁢2 1 phase to a tetragonal (𝑃⁢4 2 /𝑛⁢𝑚⁢𝑐) phase, whereas a zero stress condition drives the system from the 𝑃⁢𝑐⁢𝑎⁢2 1 phase to another orthorhombic (𝑃⁢𝑏⁢𝑐⁢𝑛) phase. These findings provide crucial insights into stress and temperature-induced phase behavior of hafnia, guiding future experimental and theoretical studies for optimizing hafnia-based ferroelectric devices.

Ferroelectric HfO2↗
DOE OSTI2026

First-principles simulation of a shocked H-He mixture along the principal Hugoniot

Recent laser-shock experiments on an H–He mixture containing 11% helium (atomic fraction) have suggested the presence of an immiscibility region inside Jupiter. Reflectivity measurements were used as the primary diagnostic of H–He demixing, with discontinuities in the optical reflectivity proposed as a signature of phase separation under conditions relevant to Jupiter's interior. Here, we investigate shock-compressed H–He using ab initio molecular dynamics simulations with optical properties evaluated within the Kubo–Greenwood formalism. The equation of state and ionic configurations were obtained using the thermal Tr 2 SCANL meta-GGA exchange–correlation (XC) functional, while optical properties were computed using the recently developed RS-KDT0 range-separated thermal hybrid XC, which provides state-of-the-art accuracy for band-gap predictions in the warm dense matter regime. The calculated reflectivity shows overall good agreement with experimental measurements; however, no discontinuity is observed at elevated temperatures. Moreover, the reflectivity predictions for the mixed system are consistent with the experimental measurements in the temperature range where the mixture is inferred to be demixed. Furthermore, these results suggest that reflectivity alone may not provide a unique or sensitive diagnostic of H–He demixing at low helium concentrations under these conditions.

Density functional calculations↗
DOE OSTI2025

New Class of Three-Nucleon Forces and Their Implications

Here, we identify a new class of three-nucleon forces that arises in the low-energy effective theory of nuclear interactions including pions. We estimate their contribution to the energy of neutron and nuclear matter and find that it can be as important as the leading-order three-nucleon forces previously considered in the literature. The magnitude of this force is set by the strength of the coupling of pions to two nucleons and is presently not well constrained by experiments. The implications for nuclei, nuclear matter, and the equation of state of neutron matter are briefly discussed.

chiral perturbation theory↗
DOE OSTI2026

Transport-based initial conditions for heavy-ion collisions at finite densities

Here, we employ the SMASH transport model to provide event-by-event initial conditions for the energy-momentum tensor and conserved charge currents in hydrodynamic simulations of relativistic heavy-ion collisions. We study the fluctuations and dynamical evolution of three conserved charge currents (net baryon, net electric charges, and net strangeness) with a four-dimensional lattice-QCD-based equation of state, NEOS-4D, in the hydrodynamic phase. Out-of-equilibrium corrections at the particlization are generalized to finite densities to ensure the conservation of energy, momentum, and the three types of charges. These theoretical developments are integrated within the X-SCAPE code as a unified framework for studying the nuclear matter properties in the Beam Energy Scan program.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Thermodynamics of spin-orbit-coupled bosons in two dimensions from the complex Langevin method

We investigate the thermal properties of interacting spin-orbit-coupled bosons with contact interactions in two spatial dimensions. To that end, we implement the complex Langevin method, motivated by the appearance of a sign problem, on a square lattice with periodic boundary conditions. We calculate the density equation of state nonperturbatively in a range of spin-orbit couplings and chemical potentials. In this work, our results show that mean-field solutions tend to underestimate the average density, especially for stronger values of the spin-orbit coupling. Additionally, the finite nature of the simulation volume induces the formation of pseudocondensates. These have been observed to be destroyed by the spin-orbit interactions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Dissipative superfluid hydrodynamics for the unitary Fermi gas

In this work we establish constraints on the temperature dependence of the shear viscosity η in the superfluid phase of a dilute Fermi gas in the unitary limit. Our results are based on analyzing experiments that measure the aspect ratio of a deformed cloud after release from an optical trap. Here, we discuss how to apply the two-fluid formalism to the unitary gas and provide a suitable parametrization of the equation of state. We show that in expansion experiments the difference between the normal and superfluid velocities remains small and can be treated as a perturbation. We find that expansion experiments favor a shear viscosity that decreases significantly in the superfluid regime. Using an exponential parametrization, we find η(T c /2T F ) ≲ 0.37[η(T c /T F )], where T c is the critical temperature and T F is the local Fermi temperature of the gas.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗