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

Robustness of neutron star merger simulations to changes in neutrino transport and neutrino-matter interactions

Binary neutron star mergers play an important role in nuclear astrophysics: Their gravitational wave and electromagnetic signals carry information about the equation of state of cold matter above nuclear saturation density, and they may be one of the main sources of r-process elements in the Universe. Neutrino-matter interactions during and after merger impact the properties of these electromagnetic signals and the relative abundances of the produced r-process elements. Existing merger simulations are, however, limited in their ability to realistically model neutrino transport and neutrino-matter interactions. Here, we perform a comparison of the impact of the use of state-of-the art two-moment or Monte Carlo transport schemes on the outcome of merger simulations, for a single binary neutron star system with a short-lived neutron star remnant [(5–10) ms]. We also investigate the use of different reaction rates in the simulations. While the best transport schemes generally agree well on the qualitative impact of neutrinos on the system, differences in the behavior of the high-density regions can significantly impact the collapse time and the properties of the hot tidal arms in this metastable merger remnant. The chosen interaction rates and transport algorithm as well as recent improvements by Radice et al. to the two-moment algorithms can all contribute to changes at the (10–30)% level in the global properties of the merger remnant and outflows. Furthermore, the limitations of previous moment schemes fixed by Radice et al. also appear sufficient to explain the large difference that we observed in the production of heavy-lepton neutrinos in a previous comparison of Monte Carlo and moment schemes in the context of a low-mass binary neutron star system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Properties of the neutron star crust: Quantifying and correlating uncertainties with improved nuclear physics

A compressible liquid-drop model (CLDM) is used to correlate uncertainties associated with the properties of the neutron star (NS) crust with theoretical estimates of the uncertainties associated with the equation of state (EOS) of homogeneous neutron and nuclear matter. For the latter, we employ recent calculations based on Hamiltonians constructed using chiral effective-field theory (χEFT). Fits to experimental nuclear masses are employed to constrain the CLDM further, and we find that they disfavor some of the χEFT Hamiltonians. The CLDM allows us to study the complex interplay between bulk, surface, curvature, and Coulomb contributions, and their impact on the NS crust. It also reveals how the curvature energy alters the correlation between the surface energy and the bulk symmetry energy. Here our analysis quantifies how the uncertainties associated with the EOS of homogeneous matter implies significant uncertainties for the composition of the crust, its proton fraction, and the volume fraction occupied by nuclei. We find that the finite-size effects impact the crust composition but have a negligible effect on the net isospin asymmetry of matter. The isospin asymmetry is largely determined by the bulk properties and the isospin dependence of the surface energy. The most significant uncertainties associated with matter properties in the densest regions of the crust, the precise location of the crust-core transition, are found to be strongly correlated with uncertainties associated with the Hamiltonians. By adopting a unified model to describe the crust and the core of NSs, we tighten the correlation between their global properties such as their mass-radius relationship, moment of inertia, crust thickness, and tidal deformability with uncertainties associated with the nuclear Hamiltonians.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Reaction plane correlated triangular flow in Au + Au collisions at $\sqrt{s_{NN}}$ = 3 GeV

Here, we measure triangular flow relative to the reaction plane at 3 GeV center-of-mass energy in Au + Au collisions at the BNL Relativistic Heavy Ion Collider. A significant v 3 signal for protons is observed, which increases for higher rapidity, higher transverse momentum, and more peripheral collisions. The triangular flow is essentially rapidity-odd with a slope at midrapidity, dv 3 /dy| (y=0) , opposite in sign compared to the slope for directed flow. No significant v 3 signal is observed for charged pions and kaons. Comparisons with models suggest that a mean field potential is required to describe these results, and that the triangular shape of the participant nucleons is the result of stopping and nuclear geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Effective field theory for dilute Fermi systems at fourth order

We discuss high-order calculations in perturbative effective field theory for fermions at low energy scales. The Fermi-momentum or k F a s as expansion for the ground-state energy of the dilute Fermi gas is calculated to fourth order, both in cutoff regularization and in dimensional regularization. For the case of spin one-half fermions we find from a Bayesian analysis that the expansion is well converged at this order for |k F a s |≲ 0.5. Furthermore, we show that Padé-Borel resummations can improve the convergence for |k F a s |≲ 1. Here, our results provide important constraints for nonperturbative calculations of ultracold atoms and dilute neutron matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Uncertainty quantification for neutrino opacities in core-collapse supernovae and neutron star mergers

We perform an extensive study of the correlations between the neutrino-nucleon inverse mean free paths (IMFPs) and the underlying equation of states (EoSs). Strong interaction uncertainties in the neutrino mean free path are investigated in different density regimes. The nucleon effective mass, the nucleon chemical potentials, and the residual interactions in the medium play an important role in determining neutrino-nucleon interactions in a density-dependent manner. Here we study how the above quantities are constrained by an EoS consistent with (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. We then study the uncertainties of both the charged current and the neutral current neutrino-nucleon inverse mean free paths due to the variation of these quantities, using the Hartree-Fock+random phase approximation method. Finally, we calculate the Pearson correlation coefficients between (i) the EoS-based quantities and the EoS-based quantities; (ii) the EoS-based quantities and the IMFPs; (iii) the IMFPs and the IMFPs. We find a strong impact of residual interactions on neutrino opacity in the spin and spin-isospin channels, which are not well constrained by current nuclear modelings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Investigating signatures of phase transitions in neutron-star cores

Neutron stars explore matter at the highest densities in the universe such that their inner cores might undergo a phase transition from hadronic to exotic phases, e.g., quark matter. Such a transition could be associated with nontrivial structures in the density behavior of the speed of sound, such as jumps and sharp peaks. Here, in this work, we employ a physics-agnostic approach to model the density dependence of the speed of sound in neutron stars and study to what extent the existence of nontrivial structures can be inferred from existing astrophysical observations of neutron stars. For this, we exhaustively study different equations of state, including as well those with explicit first-order phase transitions. We obtain a large number of different equation of states which reproduce the same astrophysical observations and obey the same physical constraints, such as mechanical stability and causality. Among them, some have nontrivial structures in the sound speed whereas others do not. We conclude that astrophysical information to date do not require the existence of a phase transition to quark matter in the density range explored in the core of neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Phase transition phenomenology with nonparametric representations of the neutron star equation of state

Astrophysical observations of neutron stars probe the structure of dense nuclear matter and have the potential to reveal phase transitions at high densities. Most recent analyses are based on parametrized models of the equation of state with a finite number of parameters and occasionally include extra parameters intended to capture phase-transition phenomenology. However, such models restrict the types of behavior allowed and may not match the true equation of state. We introduce a complementary approach that extracts phase transitions directly from the equation of state without relying on, and thus being restricted by, an underlying parametrization. We then constrain the presence of phase transitions in neutron stars with astrophysical data. Current pulsar mass, tidal deformability, and mass-radius measurements disfavor only the strongest of possible phase transitions (latent energy per particle ≳ 100 MeV). Weaker phase transitions are consistent with observations. Here we further investigate the prospects for measuring phase transitions with future gravitational-wave observations and find that catalogs of O (100) events will (at best) yield Bayes factors of ~ 10 : 1 in favor of phase transitions even when the true equation of state contains very strong phase transitions. Our results reinforce the idea that neutron star observations will primarily constrain trends in macroscopic properties rather than detailed microscopic behavior. Fine-tuned equation of state models will likely remain unconstrained in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Validity of a finite temperature expansion for dense nuclear matter

In this work we provide a new, well-controlled expansion of the equation of state of dense matter from zero to finite temperatures (𝑇) while covering a wide range of charge fractions (𝑌 𝑄 ), from pure neutron to isospin symmetric nuclear matter. Our expansion can be used to describe neutron star mergers using the equation of state inferred from neutron star observations. We discuss how knowledge from low-energy nuclear experiments and heavy-ion collisions can be directly incorporated into the expansion. We also suggest new thermodynamic quantities of interest that can be calculated from theoretical models or directly inferred by experimental data that can be used to infer the finite temperature equation of state. With our new method, we can quantify the uncertainty in our finite 𝑇 and 𝑌 𝑄 expansions without making assumptions about the underlying degrees of freedom. We can reproduce results from a microscopic equation of state up to 𝑇 = 100 MeV for baryon chemical potential 𝜇 𝐵 ≳ 1100 MeV [≈(1–2)⁢𝑛 sat ] within 5% error, with even better results for larger 𝜇 𝐵 and/or lower 𝑇. We investigate the sources of numerical and theoretical uncertainty and discuss future directions of study.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Radius and equation of state constraints from massive neutron stars and GW190814

Motivated by the unknown nature of the 2.50–2.67M⊙ compact object in the binary merger event GW190814, we study the maximum neutron star mass based on constraints from low-energy nuclear physics, neutron star tidal deformabilities from GW170817, and simultaneous mass-radius measurements of PSR J0030+045 from NICER. Our prior distribution is based on a combination of nuclear modeling valid in the vicinity of normal nuclear densities together with the assumption of a maximally stiff equation of state at high densities, a choice that enables us to probe the connection between observed heavy neutron stars and the transition density at which conventional nuclear physics models must break down. Here we demonstrate that a modification of the highly uncertain suprasaturation density equation of state beyond 2.64 times normal nuclear density is required in order for chiral effective field theory models to be consistent with current neutron star observations and the existence of 2.6M⊙ neutron stars. We also show that the existence of very massive neutron stars strongly impacts the radii of ≈ 2.0M⊙ neutron stars (but not necessarily the radii of 1.4 M⊙ neutron stars), which further motivates future NICER radius measurements of PSR J1614 – 2230 and PSR J0740+6620.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Influence of the symmetry energy on the nuclear binding energies and the neutron drip line position

A clear connection can be established between properties of nuclear matter and finite-nuclei observables, such as the correlation between the slope of the symmetry energy and the dipole polarizability, or between compressibility and the isoscalar monopole giant resonance excitation energy. Establishing a connection between realistic atomic nuclei and an idealized infinite nuclear matter leads to a better understanding of underlying physical mechanisms that govern nuclear dynamics. In this work, we aim to study the dependence of the binding energies and related quantities (e.g., location of drip lines, the total number of bound even-even nuclei) on the symmetry energy S 2 (ρ). The properties of finite nuclei are calculated by employing the relativistic Hartree-Bogoliubov model, assuming even-even axial and reflection symmetric nuclei. Calculations are performed by employing two families of relativistic energy density functionals, based on different effective Lagrangians, constrained to a specific symmetry energy at the saturation density J within the interval of 30–36 MeV. Nuclear binding energies and related quantities of bound nuclei are calculated between 8 ≤ Z ≤ 104 from the two-proton to the two-neutron drip line. As the neutron drip line is approached, the interactions with stiffer J tend to predict more bound nuclei, resulting in a systematic shift of the two-neutron drip line towards more neutron-rich nuclei. Consequentially, a correlation between the number of bound nuclei N nucl and S 2 (ρ) is established for a set of functionals constrained using the similar optimization procedures. Furthermore the direction of the relationship between the number of bound nuclei and the symmetry energy highly depends on the density under consideration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Combining Electromagnetic and Gravitational-Wave Constraints on Neutron-Star Masses and Radii

In this work, we perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass x-ray binaries (QLMXBs), photospheric radius expansion x-ray bursting sources, and x-ray timing observations of J0030+0451. With this dataset, the form of the prior distribution still has an impact on the posterior mass-radius curves and equation of state (EOS), but this impact is smaller than recently obtained when considering QLMXBs alone. We analyze the consistency of the electromagnetic data by including an “intrinsic scattering” contribution to the uncertainties, and find only a slight broadening of the posteriors. This suggests that the gravitational-wave and electromagnetic observations of neutron-star structure are providing a consistent picture of the neutron-star mass-radius curve and the EOS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Bayesian uncertainty quantification for nuclear matter incompressibility

Within a Bayesian statistical framework using the standard Skyrme-Hartree-Fock model, the maximum a posteriori (MAP) values and uncertainties of nuclear matter incompressibility and isovector interaction parameters are inferred from the experimental data of giant resonances and neutron-skin thicknesses of typical heavy nuclei. Furthermore, with the uncertainties of the isovector interaction parameters constrained by the data of the isovector giant dipole resonance and the neutron-skin thickness, we have obtained K 0 = 223$^{+7}_{–8}$ MeV at 68% confidence level using the data of the isoscalar giant monopole resonance in 208 Pb measured at the Research Center for Nuclear Physics (RCNP), Japan, and at the Texas A&M University, USA. Although the corresponding 120 Sn data gives a MAP value for K 0 about 5 MeV smaller than the 208 Pb data, there are significant overlaps in their posterior probability distribution functions.

190 ≤ A ≤ 219↗
DOE OSTI2021

Progress in Constraining Nuclear Symmetry Energy Using Neutron Star Observables Since GW170817

The density dependence of nuclear symmetry energy is among the most uncertain parts of the Equation of State (EOS) of dense neutron-rich nuclear matter. It is currently poorly known especially at suprasaturation densities partially because of our poor knowledge about isovector nuclear interactions at short distances. Because of its broad impacts on many interesting issues, pinning down the density dependence of nuclear symmetry energy has been a longstanding and shared goal of both astrophysics and nuclear physics. New observational data of neutron stars including their masses, radii, and tidal deformations since GW170817 have helped improve our knowledge about nuclear symmetry energy, especially at high densities. Based on various model analyses of these new data by many people in the nuclear astrophysics community, while our brief review might be incomplete and biased unintentionally, we learned in particular the following: (1) The slope parameter L of nuclear symmetry energy at saturation density ρ0 of nuclear matter from 24 new analyses of neutron star observables was about L≈57.7±19 MeV at a 68% confidence level, consistent with its fiducial value from surveys of over 50 earlier analyses of both terrestrial and astrophysical data within error bars. (2) The curvature Ksym of nuclear symmetry energy at ρ0 from 16 new analyses of neutron star observables was about Ksym≈−107±88 MeV at a 68% confidence level, in very good agreement with the systematics of earlier analyses. (3) The magnitude of nuclear symmetry energy at 2ρ0, i.e., Esym(2ρ0)≈51±13 MeV at a 68% confidence level, was extracted from nine new analyses of neutron star observables, consistent with the results from earlier analyses of heavy-ion reactions and the latest predictions of the state-of-the-art nuclear many-body theories. (4) While the available data from canonical neutron stars did not provide tight constraints on nuclear symmetry energy at densities above about 2ρ0, the lower radius boundary R2.01=12.2 km from NICER’s very recent observation of PSR J0740+6620 of mass 2.08±0.07M⊙ and radius R=12.2–16.3 km at a 68% confidence level set a tight lower limit for nuclear symmetry energy at densities above 2ρ0. (5) Bayesian inferences of nuclear symmetry energy using models encapsulating a first-order hadron–quark phase transition from observables of canonical neutron stars indicated that the phase transition shifted appreciably both L and Ksym to higher values, but with larger uncertainties compared to analyses assuming no such phase transition. (6) The high-density behavior of nuclear symmetry energy significantly affected the minimum frequency necessary to rotationally support GW190814’s secondary component of mass (2.50–2.67) M⊙ as the fastest and most massive pulsar discovered so far. Overall, thanks to the hard work of many people in the astrophysics and nuclear physics community, new data of neutron star observations since the discovery of GW170817 have significantly enriched our knowledge about the symmetry energy of dense neutron-rich nuclear matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Beyond $\rho ^{2/3}$ scaling: microscopic origins and multimessengers of high-density nuclear symmetry energy

The nature and Equation of State (EOS) of dense neutron-rich matter are still very poorly known, while they have broad impacts on many interesting issues in both astrophysics and nuclear physics. In particular, the nuclear symmetry energy $E_{\textrm{sym}}(\rho )$ encoding the cost to make nuclear matter more neutron-rich has been the most uncertain component of the EOS of dense neutron-rich nucleonic matter. It significantly affects the radii, tidal deformations, cooling rates, and frequencies of various oscillation modes of isolated neutron stars as well as the strain amplitude and frequencies of gravitational waves from their mergers, besides its many effects on the structures of nuclei as well as the dynamics and observables of their collisions. Siemens (1970s) observed that $E_{\textrm{sym}}(\rho )$ scales as $(\rho /\rho _0)^{2/3}$ near the saturation density $\rho _0$ of nuclear matter, since both the kinetic part and the potential contribution (quadratic in momentum) exhibit this dependence. The scaling holds if: (1) the nucleon isoscalar potential is quadratic in momentum, and (2) the isovector interaction is weakly density-dependent. After examining many empirical evidences and understanding theoretical findings in the literature, we conclude that: (1) Siemens’ $\rho ^{2/3}$ scaling is robust and serves as a valuable benchmark for both nuclear theories and experiments up to $2\rho _0$ but breaks down at higher densities, (2) Experimental and theoretical findings about $E_{\textrm{sym}}(\rho )$ up to $2\rho _0$ are broadly consistent, but uncertainties remain large for its curvature $K_{\textrm{sym}}(\rho )$ and higher-order parameters, (3) Above $2\rho _0$ , uncertainties grow due to poorly constrained spin-isospin-dependent tensor and three-body forces as well as the resulting nucleon short-range correlations. Looking forward, combining signals from both observations of neutron stars and terrestrial heavy-ion reaction experiments is the most promising path to finally constraining the high-density $E_{\textrm{sym}}(\rho )$ and the EOS of supradense neutron-rich matter. Multiple examples of community efforts to further constrain the high-density $E_{\textrm{sym}}(\rho )$ using both real and mocked data of present and future high-precision observations of neutron stars, as well as heavy-ion collisions involving high-energy rare isotopes, are briefly reviewed.

Li, Bao-An [East Texas A & M Univ., Commerce, TX (↗
DOE OSTI2025

Neutron matter from local chiral effective field theory interactions at large cutoffs

Neutron matter is an important many-body system that provides valuable constraints for the equation of state (EOS) of neutron stars. Neutron-matter calculations employing chiral effective field theory (EFT) interactions have been extensively used for this purpose. Among the various many-body methods, quantum Monte Carlo (QMC) methods stand out due to their nonperturbative nature and the achievable precision. However, QMC methods require local interactions as input, which leads to the appearance of stronger regulator artifacts compared to nonlocal interactions. To circumvent this, we employ large-cutoff interactions derived within chiral EFT (400 MeV ≤ Λ 𝑐 ≤ 700MeV) for studies of pure neutron matter. These interactions have been adjusted to nucleon-nucleon scattering phase shifts, the triton binding energy, as well as the triton 𝛽-decay half-life. We find that regulator artifacts significantly decrease with increasing cutoff, leading to a significant reduction of uncertainties in the neutron-matter EOS. We discuss implications for the symmetry energy and demonstrate how our new calculations lead to a reduction in the theoretical uncertainty of predicted neutron-star radii by up to 30% for low-mass stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Hot and dense matter equation of state probability distributions for astrophysical simulations

We add an ensemble of nuclei to the equation of state for homogeneous nucleonic matter to generate a new set of models suitable for astrophysical simulations of core-collapse supernovae and neutron star mergers. We implement empirical constraints from (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. Our model is also guided by microscopic many-body theory calculations based on realistic nuclear forces, including the zero-temperature neutron matter equation of state from quantum Monte Carlo simulations and thermal contributions to the free energy from finite-temperature many-body perturbation theory. We ensure that the parameters of our model can be varied while preserving thermodynamic consistency and the connection to experimental or observational data, thus providing a probability distribution of the astrophysical hot and dense matter equation of state. Furthermore, we compare our results with those obtained from other available equations of state. While our probability distributions indeed represent a large number of possible equations of state, we cannot yet claim to have fully explored all of the uncertainties, especially with regard to the structure of nuclei in the hot and dense medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗