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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Detectability of subsolar mass neutron stars through a template bank search

Here, we study the detectability of gravitational-wave signals from subsolar-mass binary neutron star systems by the current generation of ground-based gravitational-wave detectors. We find that finite size effects from large tidal deformabilities of the neutron stars and lower merger frequencies can significantly impact the sensitivity of the detectors to these sources. By simulating a matched-filter based search using injected binary neutron star signals with tidal deformabilities derived from physically motivated equations of state, we calculate the reduction in sensitivity of the detectors. We conclude that the loss in sensitive volume can be as high as 78.4% for an equal mass binary system of chirp mass 0.17M ⊙ , in a search conducted using binary black hole template banks. We use this loss in sensitive volume, in combination with the results from the search for subsolar-mass binaries conducted on data collected by the LIGO-Virgo observatories during their first three observing runs, to obtain a conservative upper limit on the merger rate of subsolar-mass binary neutron stars. Since the discovery of a low-mass neutron star would provide new insight into formation mechanisms of neutron stars and further constrain the equation of state of dense nuclear matter, our result merits a dedicated search for subsolar-mass binary neutron star signals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Nonequilibrium effects on stability of hybrid stars with first-order phase transitions

The stability of hybrid stars with first-order phase transitions as determined by calculating fundamental radial oscillation modes is known to differ from the predictions of the widely used Bardeen-Thorne-Meltzer criterion. We consider the effects of out-of-chemical-equilibrium physics on the radial modes and hence stability of these objects. For a barotropic equation of state, this is done by allowing the adiabatic sound speed to differ from the equilibrium sound speed. We show that doing so extends the stable branches of stellar models, allowing stars with rapid phase transitions to support stable higher-order stellar multiplets similarly to stars with multiple slow phase transitions. We also derive a new junction condition to impose on the oscillation modes at the phase transition. Termed the reactive condition, it is physically motivated, consistent with the generalized junction conditions between two phases, and has the common rapid and slow conditions as limiting cases. Unlike the two common cases, it can only be applied to nonbarotropic stars. Here, we apply this junction condition to hybrid stellar models generated using a two-phase equation of state consisting of nuclear matter with unpaired quark matter at high densities joined by a first-order phase transition and show that, like in the slow limiting case, stars that are classically unstable are stabilized by a finite chemical reaction speed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The importance of few-nucleon forces in chiral effective field theory

We study the importance of few-nucleon forces in chiral effective field theory for describing many-nucleon systems. A combinatorial argument suggests that three-nucleon forces-which are conventionally regarded as next-to-next-to-leading order-should accompany the two-nucleon force already at leading order (LO) starting with mass number $A$ ≃10–20. We find that this promotion enables the first realistic description of the 16 O ground state based on a renormalization-group-invariant LO interaction. We also performed coupled-cluster calculations of the equation of state for symmetric nuclear matter and our results indicate that LO four-nucleon forces could play a crucial role for describing heavy-mass nuclei. In conclusion, the enhancement mechanism we found is very general and could be important also in other many-body problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Equation of State and Progenitor Dependence of Stellar-mass Black Hole Formation

The core collapse of a massive star results in the formation of a proto-neutron star (PNS). If enough material is accreted onto a PNS, it will become gravitationally unstable and further collapse into a black hole (BH). We perform a systematic study of failing core-collapse supernovae in spherical symmetry for a wide range of pre-supernova progenitor stars and equations of state (EOSs) of nuclear matter. We analyze how variations in progenitor structure and the EOS of dense matter above nuclear saturation density affect the PNS evolution and subsequent BH formation. Comparisons of core collapse for a given progenitor star and different EOSs show that the path traced by the PNS in mass-specific entropy phase space ${M}_{\mathrm{grav}}^{\mathrm{PNS}}-\tilde{s}$ is well correlated with the progenitor compactness and is almost EOS independent, apart from the final end point. Furthermore, BH formation occurs, to a very good approximation, soon after the PNS overcomes the maximum gravitational mass supported by a hot NS with constant specific entropy equal to $\tilde{s}$. These results show a path to constraining the temperature dependence of the EOS through the detection of neutrinos from a failed galactic supernova.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2020

Optical Follow-Up of the Neutron Star-Black Hole Mergers S200105ae and S200115j

LIGO and Virgo’s third observing run revealed the first neutron star–black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements creating optical/near-infrared ‘kilonova’ emission. The joint gravitational wave and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter, and independently measure the local expansion rate of the Universe. Here, we present the optical follow-up and analysis of two of the only three high-significance NSBH merger candidates detected to date, S200105ae and S200115j, with the Zwicky Transient Facility. The Zwicky Transient Facility observed ~48% of S200105ae and ~22% of S200115j’s localization probabilities, with observations sensitive to kilonovae brighter than −17.5 mag fading at 0.5 mag d−1 in the g- and r-bands; extensive searches and systematic follow-up of candidates did not yield a viable counterpart. We present state-of-the-art kilonova models tailored to NSBH systems that place constraints on the ejecta properties of these NSBH mergers. We show that with observed depths of apparent magnitude ~22 mag, attainable in metre-class, wide-field-of-view survey instruments, strong constraints on ejecta mass are possible, with the potential to rule out low mass ratios, high black hole spins and large neutron star radii.

Shreya Anand↗
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↗