Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equations 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Addition of tabulated equation of state and neutrino leakage support to illinoisgrmhd

Here we have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to illinoisgrmhd, an open-source GRMHD code for dynamical spacetimes in the einstein toolkit. These new features are provided by two new, nrpy+-based codes: nrpyeos, which performs highly efficient EOS table lookups and interpolations, and nrpyleakage, which implements a new, adaptive mesh refinement (AMR)-capable neutrino leakage scheme in the einstein toolkit. We have performed a series of strenuous validation tests that demonstrate the robustness of these new codes, particularly on the Cartesian AMR grids provided by carpet. Furthermore, we show results from fully dynamical GRMHD simulations of single unmagnetized neutron stars, and magnetized binary neutron star mergers. This new version of illinoisgrmhd, as well as nrpyeos and nrpyleakage, is pedagogically documented in jupyter notebooks and fully open source. The codes will be proposed for inclusion in an upcoming version of the einstein toolkit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Optimization to Generate Equations of State for Hydrogen Production

On a high level, the larger project in question, HydroGEN, aims to develop software used for finding equations of state (EOS) to optimize catalyst configuration for H 2 production through water splitting. In particular, this summer project focused on solving the nonlinear equations used in fitting the equations. This problem involved using Python to solve a linear system with nonlinear constraints. In order for this to be achieved, Pyomo was used to build a model and the solver Ipopt, interior point optimizer, was used. Pyomo is a Python-based language developed at Sandia; it is an optimization modeling language. Rather than solving the entire problem at once, a toy problem was created, simplifying the problem down to the most important focus. This problem had a known solution, comparable to the calculated solution to assess accuracy and as progress was made towards finding solutions, complexity was gradually added to the problem. After building and solving the toy problem, it was found that it gave reasonably accurate solutions, better compared to the two existing solvers previously used with this project in terms of functionality. The solver is now ready for implementation into the project’s main software.

08 HYDROGEN↗
DOE OSTI2021

Normalizing Flows for Microscopic Many-Body Calculations: An Application to the Nuclear Equation of State

We report that normalizing flows are a class of machine learning models used to construct a complex distribution through a bijective mapping of a simple base distribution. We demonstrate that normalizing flows are particularly well suited as a Monte Carlo integration framework for quantum many-body calculations that require the repeated evaluation of high-dimensional integrals across smoothly varying integrands and integration regions. As an example, we consider the finite-temperature nuclear equation of state. An important advantage of normalizing flows is the ability to build highly expressive models of the target integrand, which we demonstrate enables precise evaluations of the nuclear free energy and its derivatives. Furthermore, we show that a normalizing flow model trained on one target integrand can be used to efficiently calculate related integrals when the temperature, density, or nuclear force is varied. This work will support future efforts to build microscopic equations of state for numerical simulations of supernovae and neutron star mergers that employ state-of-the-art nuclear forces and many-body methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Toward an accurate equation of state and B1-B2 phase boundary for magnesium oxide up to terapascal pressures and electron-volt temperatures

By applying auxiliary-field quantum Monte Carlo, we calculate the equation of state (EOS) and B1-B2 phase transition of magnesium oxide (MgO) up to 1 TPa. Here, the results agree with available experimental data at low pressures and are used to benchmark the performance of various exchange-correlation functionals in density functional theory calculations. We determine PBEsol is an optimal choice for the exchange-correlation functional and perform extensive phonon and quantum molecular-dynamics calculations to obtain the thermal EOS. Our results provide a preliminary reference for the EOS and B1-B2 phase boundary of MgO from zero up to 10,500 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Neutron Star Equation of State in Light of GW190814

The observation of gravitational waves from an asymmetric binary opens the possibility for heavy neutron stars, but these pose challenges to models of the neutron star equation of state. Here, we construct heavy neutron stars by introducing nontrivial structure in the speed of sound sourced by deconfined QCD matter, which cannot be well recovered by spectral representations. Their moment of inertia, Love number, and quadrupole moment are very small, so a tenfold increase in sensitivity may be needed to test this possibility with gravitational waves, which is feasible with third generation detectors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

GW190814: Impact of a 2.6 solar mass neutron star on the nucleonic equations of state

Is the secondary component of GW190814 the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system? This is the central question animating this paper. Covariant density functional theory provides a unique framework to investigate both the properties of finite nuclei and neutron stars, while enforcing causality at all densities. By tuning existing energy density functionals we were able to: (i) account for a 2.6 M ⊙ neutron star, (ii) satisfy the original constraint on the tidal deformability of a 1.4 M ⊙ neutron star, and (iii) reproduce ground-state properties of finite nuclei. Yet, for the class of models explored in this work, we find that the stiffening of the equation of state required to support supermassive neutron stars is inconsistent with either constraints obtained from energetic heavy-ion collisions or from the low deformability of medium-mass stars. Furthermore, we speculate that the maximum neutron star mass can not be significantly higher than the existing observational limit and that the 2.6 M ⊙ compact object is likely to be the lightest black hole ever discovered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Development of a multiphase equation of state for gallium with experiments and ab initio free-energy calculations

Here, we present a five-phase equation of state (EOS) for elemental gallium (Ga) that is developed using both experimental data and new theoretical predictions. Four experimentally observed solid phases (Ga-I, Ga-II, Ga-III, and Ga-IV) and one liquid phase are included. To improve our understanding of the thermal behavior of Ga and its phase boundaries under compression, we have performed ab initio density functional theory (DFT) free-energy calculations for the Ga-III and liquid phases, which enables us to determine the melt temperature to be 2214 ± 100 K at 110 GPa, extending significantly beyond the existing experimental melt data, which are limited to 25 GPa. In order to best describe the electron-thermal contribution, which dominates the liquid free energy at high temperatures, we have carried out averaged-atom-in-jellium DFT calculations to cover the entire temperature and density range of the EOS. The resulting multiphase Ga EOS is able to accurately reproduce a diverse variety of data, including known phase boundaries, the principal Hugoniot, low-pressure liquid isobars, and diamond-anvil-cell isotherm measurements at high pressures. It agrees more closely with key experimentally measured properties than other Ga EOS models targeted for high-pressure applications.

Wu, Christine J. [Lawrence Livermore National Labo↗
DOE OSTI2025

Two lectures on Yang-Lee edge singularity and analytic structure of QCD equation of state

These lecture notes, prepared for the 2024 XQCD PhD, provide an introduction to the analytic structure of an equation of state near a second-order phase transition and its most prominent landmark: the Yang-Lee edge singularity. In addition to discussing general properties, the notes review recent theoretical progress in locating the QCD critical point by tracking the trajectory of the Yang-Lee edge singularity.

Skokov, Vladimir (ORCID:0000000176191796)↗
DOE OSTI2020

Equation of State Constraints from Nuclear Physics, Neutron Star Masses, and Future Moment of Inertia Measurements

We explore constraints on the equation of state (EOS) of neutron-rich matter based on microscopic calculations up to nuclear densities and observations of neutron stars. In a previous work we showed that predictions based on modern nuclear interactions derived within chiral effective field theory and the observation of two-solar-mass neutron stars result in a robust uncertainty range for neutron star radii and the EOS over a wide range of densities. In this work we extend this study, employing both the piecewise polytrope extension from Hebeler et al. as well as the speed of sound model of Greif et al., and show that moment of inertia measurements of neutron stars can significantly improve the constraints on the EOS and neutron star radii.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

Equation of state of ammonia-water liquid - Derivation and planetological applications

The present least-squares fit calculation of the equation of state for ammonia-water liquid has yielded results for the zero-100 wt pct NH3, 170-300 K temperature, and zero-10 kb pressure parameter ranges. In conjunction with solid density and thermodynamic measurements, the present calculated and measured liquid densities are used to yield estimates of density and thermal expansion at 1 bar for the solid phases of ammonia's monohydrate, dihydrate, and hemihydrate between absolute zero and their respective melting points. Attention is given to the implications for icy satellite morphologic and tectonic forms of peritectic ammonia-water liquid that is neutrally buoyant relative to the corresponding solid phases.

Croft, S. K.↗
NASA NTRS1994

Rapidly rotating neutron stars in general relativity: Realistic equations of state

We construct equilibrium sequences of rotating neutron stars in general relativity. We compare results for 14 nuclear matter equations of state. We determine a number of important physical parameters for such stars, including the maximum mass and maximum spin rate. The stability of the configurations to quasi-radial perturbations is assessed. We employ a numerical scheme particularly well suited to handle rapid rotation and large departures from spherical symmetry. We provide an extensive tabulation of models for future reference. Two classes of evolutionary sequences of fixed baryon rest mass and entropy are explored: normal sequences, which behave very much like Newtonian sequences, and supramassive sequences, which exist for neutron stars solely because of general relativistic effects. Adiabatic dissipation of energy and angular momentum causes a star to evolve in quasi-stationary fashion along an evolutionary sequence. Supramassive sequences have masses exceeding the maximum mass of a nonrotating neutron star. A supramassive star evolves toward eventual catastrophic collapse to a black hole. Prior to collapse, the star actually spins up as it loses angular momentum, an effect that may provide an observable precursor to gravitational collapse to a black hole.

Cook, Gregory B.↗
DOE OSTI2025

From Existing and New Nuclear and Astrophysical Constraints to Stringent Limits on the Equation of State of Neutron-Rich Dense Matter

Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star (NS) observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. To constrain the EOS across its entire density range, this information needs to be combined consistently. However, the impact and model dependency of individual observations vary. Given their growing number, assessing the various methods is crucial to compare the respective effects on the EOS and discover potential biases. For this purpose, we present a broad compendium of different constraints and apply them individually to a large set of EOS candidates within a Bayesian framework. Specifically, we explore different ways of how chiral effective field theory and perturbative quantum chromodynamics can be used to place a likelihood on EOS candidates. We also investigate the impact of nuclear experimental constraints, as well as different radio and x-ray observations of NS masses and radii. This is augmented by reanalyses of the existing data from binary neutron star coalescences, in particular of GW170817, with improved models for the tidal waveform and kilonova light curves, which we also utilize to construct a tight upper limit of 2.39 M ⊙ on the TOV mass based on GW170817’s remnant. Our diverse set of constraints is eventually combined to obtain stringent limits on NS properties. We organize the combination in a way to distinguish between constraints where the systematic uncertainties are deemed small and those that rely on less conservative assumptions. For the former, we find the radius of the canonical 1.4 M ⊙ neutron star to be R 1.4 = 12.2 6 − 0.91 + 0.80 km and the TOV mass at M TOV = 2.2 5 − 0.22 + 0.42 M ⊙ (95% credibility). Including all the presented constraints yields R 1.4 = 12.2 0 − 0.48 + 0.50 km and M TOV = 2.3 0 − 0.20 + 0.07 M ⊙ . When comparing these limits to individual data points, we find that the quoted radius of HESS J1731-347 displays noticeable tension with other constraints. Constraining microphysical properties of the EOS proves more challenging. For instance, the symmetry energy slope is restricted to L sym = 48 − 25 + 21 MeV , where this constraint is mainly dominated by our reanalysis of the PREX-II and CREX experiment. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Equations of state, phase relations, and oxygen fugacity of the Ru-RuO 2 buffer at high pressures and temperatures

Experimental studies and measurements of inclusions in diamonds show that ferric iron components are increasingly stabilized with depth in the mantle. To determine the thermodynamic stability of such components, their concentration needs to be measured at known oxygen fugacities. The metal-oxide pair Ru and RuO 2 are ideal as an internal oxygen fugacity buffer in high-pressure experiments. Both phases remain solid to high temperatures and react minimally with silicates, only exchanging oxygen. To calculate oxygen fugacities at high pressure and temperature, however, requires information on the phase relations and equation of state properties of the solid phases. We have made in situ synchrotron X-ray diffraction measurements in a multi-anvil press on mixtures of Ru and RuO 2 to 19.4 GPa and 1473 K with which we have determined phase relations of the RuO 2 phases and derived thermal equations of state (EoS) parameters for both Ru and RuO 2 . Rutile-structured RuO 2 was found to undergo two phase transformations, first at ~7 GPa to an orthorhombic structure and then above 12 GPa to a cubic structure. The phase boundary of the cubic phase was constrained for the first time at high pressure and temperature. We have derived a continuous Gibbs free energy expression for the tetragonal and orthorhombic phases of RuO 2 by fitting the second-order phase transition boundary and P-V-T data for both phases, using a model based on Landau theory. The transition between the orthorhombic and cubic phases was then used along with EoS terms derived for both phases to determine a Gibbs free energy expression for the cubic phase. We have used these data to calculate the oxygen fugacity of the Ru + O 2 = RuO 2 equilibrium, which we have parameterized as a single polynomial across the stability fields of all three phases of RuO 2 . The expression is log 10 f O 2 (Ru – RuO 2 ) = (7.782 – 0.00996P + 0.001932P 2 – 3.76 × 10 –5 P 3 ) + (–13 763 + 592P – 3.955P 2 )/T + (–1.05 × 10 6 – 4622P)/T 2 , which should be valid from room pressure up to 25 GPa and 773–2500 K, with an estimated uncertainty of 0.2 log units. Our calculated fO 2 is shown to be up to 1 log unit lower than estimates that use previous expressions or ignore EoS terms.

Armstrong, Katherine↗
DOE OSTI2024

Enhanced contribution of the pairing gap to the QCD equation of state at large isospin chemical potential

In this paper QCD at large isospin density is studied, which is known to be in the superfluid state with Cooper pairs carrying the same quantum number as pions. The gap equation derived from the perturbation theory up to the next-to-leading-order corrections is solved. The pairing gap at large isospin chemical potential is found to be enhanced compared to the color-superconducting gap at large baryon chemical potential due to the 2 difference in the exponent arising from the stronger attraction in one-gluon exchange in the singlet channel. Then, using the gap function, the contribution of the condensation energy of the superfluid state to the QCD equation of state is evaluated. At isospin chemical potential of a few GeV, where the lattice QCD and the perturbative QCD can be both applied, the effect of the condensation energy becomes dominant even compared to the next-to-leading order corrections to the pressure in the perturbation theory. It resolves the discrepancy between the recent lattice QCD results and the perturbative QCD result. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2020

Equation of state of hot, dense magnesium derived with first-principles computer simulations

Using two first-principles computer simulation techniques, path integral Monte Carlo and density functional theory molecular dynamics, we derive the equation of state of magnesium in the regime of warm dense matter, with densities ranging from 0.43 to 86.11 g cm –3 and temperatures from 20,000 K to 5×10 8 K. These conditions are relevant for the interiors of giant planets and stars as well as for shock compression measurements and inertial confinement fusion experiments. Here, we study ionization mechanisms and the electronic structure of magnesium as a function of density and temperature. We show that the L shell electrons, 2s and 2p energy bands, merge at high densities. This results in gradual ionization of the L-shell with increasing density and temperature. In this regard, Mg differs from MgO, which is also reflected in the shape of its principal shock Hugoniot curve. For Mg, we predict a single broad pressure-temperature region, where the shock compression ratio is approximately 4.9. Mg thus differs from Si and Al plasmas that exhibit two well-separated compression maxima on the Hugoniot curve for L and K shell ionizations. Finally, we study multiple shocks and effects of preheat and precompression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Thermal equation of state for zoisite: Implications for the transportation of water into the upper mantle and the high-velocity anomaly in the Farallon plate

Seismic tomography studies have revealed a high-velocity anomaly at depths between 100 km and 300 km in the Farallon plate. However, the reasons for the high-velocity anomaly continue to be debated. An analysis of the mineral proportions of eclogites exhumed from the Farallon plate shows that the average amount of zoisite in eclogite is ∼16.0 vol%. Therefore, the presence of zoisite eclogite needs to be considered to explain the high-velocity anomaly of the Farallon plate. However, the thermal equation of state and stability of zoisite have not been fully investigated under high pressure−temperature (P-T) conditions. We investigated the high-pressure and high-temperature behavior of natural zoisite utilizing synchrotron single crystal-X-ray diffraction (XRD). The results indicate that zoisite is metastable up to 24.8 GPa and 700 K. We obtained the ambient unit-cell volume V0 = 901.26(3) Å3 by synchrotron single crystal-XRD measurement. We also fitted the pressure-volume-temperature data to a high-temperature Birch-Murnaghan equation of state and obtained the zero-pressure bulk modulus K0 = 134.7(8) GPa, the temperature derivative of the bulk modulus (∂K/∂T)P = −0.011(4) GPa/K, and the thermal expansion coefficients α0V = 1.9(7) × 10−5 K−1 and α1V = 3(2) × 10−8 K−2 when the pressure derivative of bulk modulus (K′0) is fixed at 4. By incorporating the results from previous studies, we calculated the density and bulk sound velocity profiles of zoisite eclogite along the Farallon plate geotherm. Finally, we infer that zoisite could carry water to depths of ∼300 km within cold subducting slabs and that the bulk sound velocity of typical zoisite eclogite with ∼61.0 vol% omphacite, ∼23.0 vol% garnet, and ∼16.0 vol% zoisite could cause the high-velocity anomalies at depths of 100−300 km in the Farallon plate.

Geology↗
DOE OSTI2024

Chiral condensate and the equation of state at nonzero baryon density from the hadron resonance gas model with a repulsive mean field

We study the QCD equation of state and the chiral condensate using the hadron resonance gas model with repulsive mean-field interactions. We find that the repulsive interactions improve the agreement with the lattice results on the derivatives of the pressure with respect to the baryon chemical potential up to eighth order. From the temperature dependence of the chiral condensate we estimate the crossover temperature as a function of baryon chemical potential, T p c ( μ B ) . We find that the chiral crossover line starts to deviate significantly from the chemical freeze-out line already for μ B > 400 MeV . Furthermore, we find that the chiral pseudocritical line can be parametrized as T p c ( μ B ) / T p c ( 0 ) = 1 − κ 2 [ μ B / T p c ( 0 ) ] 2 − κ 4 [ μ B / T p c ( 0 ) ] 4 with κ 2 = 0.0150 ( 2 ) and κ 4 = 3.1 ( 6 ) × 10 − 5 , which are in agreement with lattice QCD results for small values of μ B . For the first time we find a tiny but nonzero value of κ 4 in our study. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Impact of the equation of state on f - and p - mode oscillations of neutron stars

We investigate the impact of the neutron-star matter equation of state on the f- and p1-mode oscillations of neutron stars obtained within the Cowling approximation and linearized general relativity. The f- and p 1 -mode oscillation frequencies, and their damping times are calculated using representative sets of Skyrme Hartree-Fock and relativistic mean-field models, all of which reproduce nuclear systematics and support 2M ⊙ neutron stars. Our study shows strong correlations between the frequencies of f- and p 1 -modes and their damping times with the pressure of β-equilibrated matter at densities equal to or slightly higher than the nuclear saturation density ρ 0 . Such correlations are found to be almost independent of the composition of the stars. The frequency of the p 1 -mode of 1.4M ⊙ star is strongly correlated with the slope of the symmetry energy L 0 and β-equilibrated pressure at density ρ 0 . Compared to GR calculations, the error in the Cowling approximation for the f-mode is about 30% for neutron stars of low mass, whereas it decreases with increasing mass. Here, the accuracy of the p 1 -mode is better than 15% for neutron stars of maximum mass, and improves for lower masses and higher number of radial nodes. I.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗