Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear 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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Astrophysical Equation-of-State Constraints on the Color-Superconducting Gap

We demonstrate that astrophysical constraints on the dense-matter equation of state place an upper bound on the color-superconducting gap in dense matter above the transition from nuclear matter to quark matter. Pairing effects in the color-flavor locked quark matter phase increase the pressure at high density, and if this effect is sufficiently large then the requirements of causality and mechanical stability make it impossible to reach such a pressure in a way that is consistent with what is known at lower densities. The intermediate-density equation of state is inferred by considering extensions of chiral effective field theory to neutron star densities, and conditioning these using current astrophysical observations of neutron star radius, maximum mass, and tidal deformability (PSR J⁢0348+0432, PSR J1624-2230, PSR J⁢0740+6620, GW170817). At baryon number chemical potential μ = 2.6 GeV we find a 95% upper limit on the color-flavor locked pairing gap Δ of 457 MeV using overly conservative assumptions and 216 MeV with more reasonable assumptions. Furthermore, this constraint may be strengthened by future astrophysical measurements as well as by future advances in high-density QCD calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Equation of State and Energy Loss of Hot and Dense Quark-Gluon matter from Holographic Black Holes

By using gravity/gauge correspondence, we construct a holographic model, constrained to mimic the lattice QCD equation of state at zero density, to investigate the temperature and baryon chemical potential dependence of the equation of state. We also obtained the energy loss of light and heavy partons within the hot and dense plasma represented by the heavy quark drag force, Langevin diffusion coefficients and jet quenching parameter at the critical point and across the first-order transition line predicted by the model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Identifying the nature of the QCD transition in heavy-ion collisions with deep learning

In this proceeding, we review our recent work using deep convolutional neural network (CNN) to identify the nature of the QCD transition in a hybrid modeling of heavy-ion collisions. Within this hybrid model, a viscous hydrodynamic model is coupled with a hadronic cascade “after-burner”. As a binary classification setup, we employ two different types of equations of state (EoS) of the hot medium in the hydrodynamic evolution. The resulting final-state pion spectra in the transverse momentum and azimuthal angle plane are fed to the neural network as the input data in order to distinguish different EoS. To probe the effects of the fluctuations in the event-by-event spectra, we explore different scenarios for the input data and make a comparison in a systematic way. We observe a clear hierarchy in the predictive power when the network is fed with the event-by-event, cascade-coarse-grained and event-fine-averaged spectra. The carefully-trained neural network can extract high-level features from pion spectra to identify the nature of the QCD transition in a realistic simulation scenario.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

EOSPAC User's Manual: Version 6.5 Second Edition (Rev. 3)

The EOSPAC utility package is a collection of interface routines, which can be used to access the SESAME data library and perform various data adjustments and interpolations on the SESAME data. The SESAME data library contains both thermodynamic (e.g., equation of state) and transport coefficients (e.g., opacity and conductivity). Note, for simplicity, the term EOS (equation of state) used herein includes both thermodynamic variables and transport coefficients. The EOSPAC utility package is designed to be used by physics codes (henceforth ”host codes”) written in multiple languages and on multiple platforms. The remainder of this manual is organized into several sections. Chapter 2 discusses conventions such as data organization and routine names. Chapter 3 provides a general overview of basic theory and models implemented within EOSPAC. Chapter 4 provides a general overview of how to use the EOSPAC interface library. Chapters 5 to 7 describe the public interfaces of EOSPAC in detail. Chapter 8 provides a brief introduction to some related tools, which may be of use to the user. Chapter 9 provides details related to some selected numerical features of EOSPAC. Chapter 10 gives examples for using the interface routines described in chapters 5 to 7. Chapter 11 provides technical support contact information. Chapter 12 contains a brief set of acknowledgments. Chapter 13 contains a list of referenced documents. Finally, chapter 14 lists the “table types: mnemonic conventions”, “table types: grouped by category, sorted by name”, “table types: eospac version 5 cross reference”, “options: setup phase”, “data information parameters”, “meta-data information parameters”, “options: interpolation phase”, and the “error codes”.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

ARES v1.x - Performance Portable Tool to Simulate Supernovae based on Parthenon Framework

Historically, codes for simulating supernovae (such as Arepo, FLASH or LEAFS) have been at the forefront of scientific high-performance computing to the immense computational resources required for full 3D simulations. However, given the shift towards heterogenous HPC architectures, many current-generation codes are at the risk of losing their competitiveness as they are only designed to run on homogeneous CPU-only systems. There exist several efforts to enable these codes for GPU’s, however, these efforts only consider specific architectures or vendors (e.g., implement only CUDA or HIP), limiting themselves to a small range of exascale computing systems. Frameworks such as Kokkos aim to provide a framework which is agnostic of the targeted architecture, enabling the development of performant and portable code. In the Ares code, we develop a performance portable tool to simulate supernovae based on the Parthenon Framework, which in turn uses Kokkos in the background. Here, the Parthenon Framework provides an interface to the underlying mesh-refinement routines, which form the backbone of our code. In addition, we incorporate the already existing Singularity-EOS toolkit to provide us with various equations of state, primarily the Helmholtz equation of state. We also include the JINA Reaclib as a basis for our nuclear network solver. Finally, we implement a gravity solver to complete the required physics. This setup will provide us with a minimal code base to simulate supernova in a similar style to the tried-and-tested Arepo code, but in a futureproof performance portable framework.

Lim, Hyun↗
DOE OSTI2021

Strangeness-neutral equation of state for QCD with a critical point

We present a strangeness-neutral equation of state for QCD that exhibits critical behavior and matches lattice QCD results for the Taylor-expanded thermodynamic variables up to fourth-order in µ B /T. It is compatible with the SMASH hadronic transport approach and has a range of temperatures and baryonic chemical potentials relevant for phase II of the Beam Energy Scan at RHIC. We provide an updated version of the software BES-EoS, which produces an equation of state for QCD that includes a critical point in the 3D Ising model universality class. This new version also includes isentropic trejectories and the critical contribution to the correlation length. Since heavy-ion collisions have zero global net-strangeness density and a fixed ratio of electric charge to baryon number, the BES-EoS is more suitable to describe this system. Comparison with the previous version of the EoS is thoroughly discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

General features of the stellar matter equation of state from microscopic theory, new maximum-mass constraints, and causality

The profile of a neutron star probes a very large range of densities, from the density of iron up to several times the density of saturated nuclear matter, and thus no theory of hadrons can be considered reliable if extended to those regions. We emphasize the importance of taking contemporary ab initio theories of nuclear and neutron matter as the baseline for any extension method, which will unavoidably involve some degree of phenomenology. We discuss how microscopic theory, on the one end, with causality and maximum-mass constraints, on the other, set strong boundaries to the high-density equation of state. We present our latest neutron star predictions where we combine polytropic extensions and parametrizations guided by speed of sound considerations. The predictions we show include our baseline neutron star cooling curves.

chiral effective field theory↗
DOE OSTI2020

GW190814's Secondary Component with Mass 2.50–2.67 M ⊙ as a Superfast Pulsar

Using Stergioulas's RNS code for investigating fast pulsars with Equation of States (EOSs) on the causality surface (where the speed of sound equals that of light) of the high-density EOS parameter space satisfying all known constraints from both nuclear physics and astrophysics, we show that one possible explanation for the GW190814's secondary component of mass (2.50-2.67) M ⊙ is that it is a super-fast pulsar spinning faster than 971 Hz about 42% below its Kepler frequency. If confirmed, it would be the fastest pulsar with the highest mass observed presently. There is a large and physically allowed EOS parameter space below the causality surface where pulsars heavier than 2.50 M ⊙ are supported if they can rotate even faster with critical frequencies depending strongly on the high-density behavior of nuclear symmetry energy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A multi-wavelength streaked optical pyrometer for dynamic shock compression measurements above 2500 K

This work presents a concept and implementation of a calibrated streaked spectral pyrometer (SSP) temperature diagnostic used in dynamically driven shock experiments on a two-stage gas gun. This system relies upon measuring the total system response using a tunable monochromator, a NIST-traceable calibrated power meter, and a SSP. The diagnostic performance is validated against previously measured temperatures of shock driven z-cut quartz at 99 GPa and 93 GPa. The results are found to agree with the literature to within 5% and are discussed in this manuscript. The experimental setup utilizes measurements from multiple SSP systems per sample, providing several independent measurements and substantially increasing confidence in the extrapolated shock driven sample temperature.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

New Covariant Density Functionals of Nuclear Matter for Compact Star Simulations

We generate three families of extended covariant density functionals of nuclear matter that have varying slope of symmetry energy and skewness at nuclear saturation density, but otherwise share the same basic parameters (symmetry energy, compressibility, saturation parameters, etc.) with the standard DDME2, DD2, and MPE functionals. Tables of the parameters of these new density functionals are given, which can be straightforwardly used in DDME2, DD2, and MPE parameterization-based codes. Furthermore, we provide tables of a large number of equations of state (81 for each family) that can be used in astrophysical simulations to assess the impact of variations of the not-well-known slope of symmetry energy and skewness of nuclear systems on the astrophysics of compact objects. We also provide tables of computed integral parameters (mass, radius, and tidal deformability) that can be used, e.g., for modeling gravitational waveforms. Finally, for the extended DDME2-based parameterization, we implement a first-order phase transition to quark matter to obtain a family of equations of state that accommodates a phase transition to quark matter. Analogous tables of the equations of state and integral parameters are provided for this case as well.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

QCD Equilibrium and Dynamical Properties from Holographic Black Holes

By using gravity/gauge correspondence, we employ an Einstein-Maxwell-Dilaton model to compute the equilibrium and out-of-equilibrium properties of a hot and baryon rich strongly coupled quark-gluon plasma. The family of 5-dimensional holographic black holes, which are constrained to mimic the lattice QCD equation of state at zero density, is used to investigate the temperature and baryon chemical potential dependence of the equation of state. We also obtained the baryon charge conductivity, and the bulk and shear viscosities with a particular focus on the behavior of these observables on top of the critical end point and the line of first order phase transition predicted by the model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Viscosities of the Baryon-Rich Quark-Gluon Plasma from Beam Energy Scan Data

Here, this work presents the first Bayesian inference study of the (3 + 1)D dynamics of relativistic heavy-ion collisions and quark-gluon plasma viscosities using an event-by-event (3 + 1)D hydrodynamics + hadronic transport theoretical framework and data from the Relativistic Heavy Ion Collider Beam energy scan program. Robust constraints on initial state nuclear stopping and the baryon chemical potential-dependent shear viscosity of the produced quantum chromodynamic (QCD) matter are obtained. The specific bulk viscosity of the QCD matter is found to exhibit a preferred maximum around $\sqrt{s_{NN}}$ = 19.6 GeV . This result allows for the alternative interpretation of a reduction (and/or increase) of the speed of sound relative to that of the employed lattice-QCD based equation of state for net baryon chemical potential μ B ~0.2(0.4) GeV .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Connecting GRBs from binary neutron star mergers to nuclear properties of neutron stars

The fate of the binary neutron star (NS) merger remnants hinges sensitively upon the NS equation of state and the threshold mass, M ls , that separates a long-lived from a short-lived NS remnant. The nature of the electromagnetic counterparts is also influenced by the remnant type, particularly in determining whether a gamma-ray burst from a compact binary merger (cbGRB) is of short or long duration. We propose a novel approach to probe Mls by linking it to the estimated observed ratio of long to short cbGRBs. We find that current observations broadly favor a relatively high value for this transition, M ls ≃ 1.3M TOV , for which M TOV ≲ 2.6M ⊙ , consistent with numerical simulations, as also shown here. Furthermore, our results disfavor nuclear physics scenarios that would lead to catastrophic pressure loss at a few times nuclear density and temperatures of tens of MeV, leading to a rapid gravitational collapse of binaries with total mass M ≲ 1.3M TOV . Future individual gravitational wave events with on-axis cbGRBs can further bound Mls.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Bayesian Exploration of Phenomenological EoS of Neutron/Hybrid Stars with Recent Observations

The description of the stellar interior of compact stars remains as a big challenge for the nuclear astrophysics community. The consolidated knowledge is restricted to density regions around the saturation of hadronic matter ρ 0 =2.8 × 10 14 g cm -3 , regimes where our nuclear models are successfully applied. As one moves towards higher densities and extreme conditions up to the quark/gluons deconfinement, little can be said about the microphysics of the equation of state (EoS). Here, we employ a Markov Chain Monte Carlo (MCMC) strategy to access the variability at high density regions of polytropic piecewise models for neutron star (NS) EoS or possible hybrid stars, i.e., a NS with a small quark-matter core. With a fixed description of the hadronic matter for low density, below the nuclear saturation density, we explore a variety of models for the high density regimes leading to stellar masses near to 2.5 M ⊙ , in accordance with the observations of massive pulsars. The models are constrained, including the observation of the merger of neutrons stars from VIRGO-LIGO and with the pulsar observed by NICER. In addition, we also discuss the possibility of the use of a Bayesian power regression model with heteroscedastic error. The set of EoS from the Laser Interferometer Gravitational-Wave Observatory (LIGO) was used as input and treated as the data set for the testing case.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Modeling Neutron Star Matter in the Age of Multimessenger Astrophysics

Abstract The interpretation of available and forthcoming data obtained from multimessenger astrophysical observations—potentially providing unprecedented access to neutron star properties—will require the development of novel, accurate theoretical models of dense matter. Of great importance, in this context, will be the capability for devising a description of thermal effects applicable to the study of quantities other than the equation of state, such as the transport coefficients and the neutrino mean free path in the nuclear medium. A formalism based on correlated basis states and the cluster expansion technique has been previously employed to derive a well-behaved effective interaction—suitable for use in standard perturbation theory—from a state-of-the-art nuclear Hamiltonian, including phenomenological two- and three-nucleon potentials. Here, we provide a comprehensive and self-contained account of the extension of this approach to the treatment of finite-temperature effects, and report the results of numerical calculations of a number of properties of nuclear matter with arbitrary neutron excess and temperature up to 50 MeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Ab initio calculations of the rate of carrier trapping and release at dopant sites in NaI: Tl beyond the harmonic approximation

We present ab initio calculations of electron capture and release coefficients as a function of temperature for Tl dopants in the widely used scintillator NaI, which is a soft ionic crystal. The modelled capture and release events occur by transitions mediated by multiple phonon absorption and emission between states with significantly different local geometries around the trapping site. We demonstrate that such transitions are not well-described by the normal harmonic approximation to the nuclear dynamics. We avoid the harmonic approximation by numerically solving the vibrational Schrödinger equation for the motion along/of a single effective phonon coordinate. The localized trapped state is not correctly described by semilocal density functionals, so we employ a hybrid functional tuned to reproduce the band gap and lattice constant of the host material. Our calculations, which combine plane-wave and cluster electronic structure methods in a novel way, are in reasonable agreement with available experimental data for detrapping and predict an unusual temperature dependence for trapping rates.

Prange, Micah P.↗
DOE OSTI2025

Impact of Magnetic-field-driven Anisotropies on the Equation of State Probed in Neutron Star Mergers

Binary neutron star mergers can produce extreme magnetic fields, some of which can lead to strong magnetar-like remnants. While strong magnetic fields have been shown to affect the dynamics of outflows and angular momentum transport in the remnant, they can also crucially alter the properties of nuclear matter probed in the merger. In this work, we provide a first assessment of the latter, determining the strength of the pressure anisotropy caused by Landau-level quantization and the anomalous magnetic moment. To this end, we perform the first numerical relativity simulation with a magnetic polarization tensor and a magnetic-field-dependent equation of state using a new algorithm we present here, which also incorporates a mean-field dynamo model to control the magnetic field strength present in the merger remnant. Our results show that—in the most optimistic case—corrections to the anisotropy can be in excess of 10% and are potentially largest in the outer layers of the remnant. This work paves the way for a systematic investigation of these effects.

General relativity↗
DOE OSTI2023

Principal Component Analysis of azimuthal flow in intermediate-energy heavy-ion reactions

Principal Component Analysis (PCA) via Singular Value Decomposition (SVD) of large datasets is an adaptive exploratory method to uncover natural patterns underlying the data. Several recent applications of the PCA-SVD to event-by-event single-particle azimuthal angle distribution matrices in ultra-relativistic heavy-ion collisions at RHIC-LHC energies indicate that the sine and cosine functions chosen a priori in the traditional Fourier analysis are naturally the most optimal basis for azimuthal flow studies according to the data itself. We perform PCA-SVD analyses of mid-central Au+Au collisions at $E$ $beam$ / $A$ =1.23 GeV simulated using an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model to address the following two questions: (1) if the principal components of the covariance matrix of nucleon azimuthal angle distributions in heavy-ion reactions around 1 GeV/nucleon are naturally sine and/or cosine functions and (2) what if any advantages the PCA-SVD may have over the traditional flow analysis using the Fourier expansion for studying the EOS of dense nuclear matter. In conclusion, we find that (1) in none of our analyses the principal components come out naturally as sine and/or cosine functions, (2) while both the eigenvectors and eigenvalues of the covariance matrix are appreciably EOS dependent, the PCA-SVD has no apparent advantage over the traditional Fourier analysis for studying the EOS of dense nuclear matter using the azimuthal collective flow in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗