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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Constraints on the Muon Fraction and Density Profile in Neutron Stars

Muons in neutron stars (NSs) play especially important roles in addressing several interesting new physics questions associated with detecting as well as understanding interactions and astrophysical effects of muonphilic dark matter particles. The key model inputs for studying the latter are the total muon mass M μ , the muon mass fraction M μ /M NS over the NS mass M NS , and the muon radial density profile ρ μ (r) in NSs of varying masses. We investigate these quantities within a minimum model for the core of NSs consisting of neutrons, protons, electrons, and muons using an explicitly isospin-dependent parametric equation of state (EOS) constrained by available nuclear laboratory experiments and the latest astrophysical observations of NS masses, radii, and tidal deformabilities. Here, we found that the absolutely maximum muon mass M μ and its mass fraction M μ /M NS in the most massive NSs allowed by causality are about 0.025 M ⊙ and 1.1%, respectively. For the most massive NS of mass 2.14 M ⊙ observed so far, they reduce to about 0.020 M ⊙ and 0.9%, respectively. We also study respective effects of individual parameters describing the EOS of high-density neutron-rich nucleonic matter on the muon contents in NSs with varying masses. We found that the most important but uncertain nuclear physics ingredient for determining the muon contents in NSs is the high-density nuclear symmetry energy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Extending a scaling equation of state to QCD

Whether quantum chromodynamics (QCD) exhibits a phase transition at finite temperature and density is an open question. It is important for hydrodynamic modeling of heavy ion collisions and neutron-star mergers. Lattice QCD simulations have definitively shown that the transition from hadrons to quarks and gluons is a crossover when the baryon chemical potential is zero or small. Here, we combine the parametric scaling equation of state, usually associated with the three-dimensional Ising model, with a background equation of state based on a smooth crossover from hadrons to quarks and gluons. Comparison to experimental data from the Beam Energy Scan II at the Relativistic Heavy Ion Collider or in heavy ion experiments at other accelerators may allow the critical exponents and amplitudes in the scaling equation of state to be determined for QCD if a critical point exists.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Impact of the newly revised gravitational redshift of x-ray burster GS 1826-24 on the equation of state of supradense neutron-rich matter

Thanks to the recent advancement in producing rare isotopes and measuring their masses with unprecedented precision, the updated nuclear masses around the waiting-point nucleus 64 Ge in the rapid-proton capture process have led to a significant revision of the surface gravitational redshift of the neutron star (NS) in GS 1826-24 by refitting its x-ray burst light curve using Modules for Experiments in Stellar Astrophysics (MESA). The resulting NS compactness ξ is between 0.183 and 0.259 at 95% confidence level, and its upper boundary is significantly smaller than the maximum ξ previously known. Incorporating these new data within a comprehensive Bayesian statistical framework, we investigate its impact on the Equation of State (EOS) of supradense neutron-rich matter and the required spin frequency for GW190814’s minor m 2 with mass 2.59 ± 0.05⁢M ⊙ to be a rotationally stable pulsar. We found that the EOS of high-density symmetric nuclear matter (SNM) has to be softened significantly while the symmetry energy at supersaturation densities stiffened compared to our prior knowledge from earlier analyses using data from both astrophysical observations and terrestrial nuclear experiments. In particular, the skewness J 0 characterizing the stiffness of high-density symmetric nuclear matter (SNM) decreases significantly, while the slope L, curvature K sym , and skewness Jsym of nuclear symmetry energy all increase appreciably compared to their fiducial values. Here, we also found that the most probable spin rate for the m 2 to be a stable pulsar is very close to its mass-shedding limit once the revised redshift data from GS 1826-24 is considered, making the m 2 unlikely the most massive NS observed so far.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Properties of first-order hadron-quark phase transition from inverting neutron star observables

By inverting the observational data of several neutron star observables in the three-dimensional parameter space of the constant speed of sound (CSS) model while fixing all hadronic equation-of-state parameters at their currently known most probable values, we constrain the three parameters of the CSS model and their correlations. Using two lower radius limits of $R$ $2.01$ = $11.41$ km and $R$ $2.01$ = $12.2$ km for PSR J0740+6620 obtained from two independent analyses using different approaches by the Neutron Star Interior Composition Explorer (NICER) Collaboration, the speed of sound squared $c$ $^2_{QM}$ in quark matter is found to have a lower limit of 0.35 or 0.43 in units of $c$ $2$ , above its conformal limit of $c$ $^2_{QM}$ < $1/3$. Furthermore, an approximately linear correlation between the first-order hadron-quark transition density $ρ$ $t$ and its strength $Δε$ is found. Moreover, the presence of twin stars is deemed improbable by the present work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Symmetry-Energy Dependence of the Bulk Viscosity of Nuclear Matter

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

Neutron stars & pulsars↗
DOE OSTI2026

r -process Nucleosynthesis from Hyperaccreting Neutron Stars in Common Envelopes

We investigate nuclear reactions and feedback in hyperaccreting neutron star environments, considering accretion rates in the range 0.3–3 × 10 4 M ⊙ yr −1 , typical of short-period compact-object binaries in common envelopes. Our models account for weak reactions, neutrino energy loss, nuclear energy release, pair production, degenerate equations of state, and general relativistic hydrodynamics. Depending on the accretion rates, these systems can develop both proton- and neutron-rich atmospheres with strong convective instabilities linking the neutrino sphere to the outgoing accretion shock inside the radiation trapping zone. Convection drives nucleons through multiple heating and cooling cycles, with photodisintegration dominating during the heating phase and heavy element synthesis during the cooling phase, ejecting material with abundances that depend on the accretion rate and depth of the final decompression trajectory. The turbulent nature of convective currents is conducive to creating a wide range of nuclear products through a variety of effects, including nuclear statistical equilibrium freeze-out and the r-, p-, and γ-processes. We also observe a novel multistep process in reheated trajectories, consisting of proton-capture and photodissociation reactions operating on r-process seeds, producing overall neutron-deficient isotopes. A significant amount of infalling gas experiences high entropy and short (millisecond) freeze-out timescales capable of making r-process elements with high overabundances through a disequilibrium effect between neutrons and α-particles that does not require an excess of neutrons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of 208 Pb with Minimal Modeling Assumptions

The symmetry energy and its density dependence are crucial inputs for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of 0.283 ± 0.071 fm for the neutron-skin thickness of 208 Pb, implying a slope parameter L = 106 ± 37 MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy S 0 , L, and $R^{^{208}Pb}_{skin}$ directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER clearly favor smaller values of the neutron skin and L, as well as negative symmetry incompressibilities. Finally, combining astrophysical data with PREX-II and chiral effective field theory constraints yields S 0 = $33.0^{+2.0}_{-1.8}$ MeV , L = $53^{+14}_{-15}$ MeV, and $R^{^{208}Pb}_{skin} = 0.17^{+0.04}_{-0.04}$ fm.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Frequency and Time Domain Nuclear-Electronic Orbital Equation-of-Motion Coupled Cluster Methods: Combination Bands and Electronic-Protonic Double Excitations

The accurate description of excited vibronic states is important for modeling a wide range of photoinduced processes. The nuclear-electronic orbital (NEO) approach, which treats specified protons on the same level as the electrons, can describe excited electronic-protonic states. Herein the multicomponent equation-of-motion coupled cluster with singles and doubles (NEO-EOM-CCSD) method and its time-domain counterpart, TD-NEO-EOM-CCSD, are developed and implemented. The application of these methods to the HCN molecule highlights their capabilities. These methods predict qualitatively reasonable energies and intensities for a combination band corresponding to simultaneous excitation of two vibrational modes, as well as an overtone. These methods also describe states with double excitation character, such as excited electronic-protonic states corresponding to the simultaneous excitation of an electron and a proton. The ability of the NEO-EOM-CCSD method and its time-dependent counterpart to describe combination bands, overtones, and double excitations will enable a wide range of photochemical applications.

Pavosevic, Fabijan↗
DOE OSTI2021

Detailed examination of astrophysical constraints on the symmetry energy and the neutron skin of 208 Pb with minimal modeling assumptions

We report the symmetry energy and its density dependence are pivotal for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of 0.283 ± 0.071 fm for the neutron-skin thickness of 208 Pb, $R^{^{208}Pb}_{skin}$, implying a symmetry-energy slope parameter L of 106 ± 37 MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy S 0 , L, and $R^{^{208}Pb}_{skin}$ directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER favor smaller values of the neutron skin and L, as well as negative symmetry incompressibilities. Combining astrophysical data with chiral effective field theory (χEFT) and PREX-II constraints yields S 0 = 33.0$^{+2.0}_{-1.8}$ MeV, L = 53$^{+14}_{-15}$ MeV, and $R^{^{208}Pb}_{skin}$ = 0.17 $^{+0.04}_{-0.04}$ fm. We also examine the consistency of several individual χEFT calculations with astrophysical observations and terrestrial experiments. We find that there is only mild tension between χ EFT, astrophysical data, and PREX-II's $R^{^{208}Pb}_{skin}$ measurement (p value = 12.3%) and that there is excellent agreement between χEFT, astrophysical data, and other nuclear experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Large and massive neutron stars: Implications for the sound speed within QCD of dense matter

Here, the NASA telescope NICER has recently measured x-ray emissions from the heaviest of the precisely known two-solar mass neutron stars, PSR J0740 + 6620. Analysis of the data suggests that PSR J0740 + 6620 has a radius in the range of R 2.0 ≈ (11.4 – 16.1) km at the 68% credibility level. In this article, we study the implications of this analysis for the sound speed in the high-density inner cores by using recent chiral effective field theory (χEFT) calculations of the equation of state at next-to-next-to-next-to-leading order to describe outer regions of the star at modest density. We find that the lower bound on the maximum speed of sound in the inner core, min {c$^{2}_{s,max}$}, increases rapidly with the radius of massive neutron stars. If χ EFT remains an efficient expansion for nuclear interactions up to about twice the nuclear saturation density, R 2.0 ≥ 13 km requires min {c$^{2}_{s,max}$} ≥ 0.562 and 0.442 at the 68% and 95% credibility level, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Spin-dependent interactions and heavy-quark transport in the quark-gluon plasma

We extend a previously constructed T-matrix approach to the quark-gluon plasma (QGP) to include the effects of spin-dependent interactions between partons. Following earlier work within the relativistic quark model, the spin-dependent interactions figure as relativistic corrections to the Cornell potential. When applied to the vacuum spectroscopy of quarkonia, in particular their mass splittings in S- and P-wave states, the issue of the Lorentz structure of the confining potential arises. We confirm that a significant admixture of a vector interaction (to the previously assumed scalar interaction) improves the description of the experimental mass splittings. The temperature corrections to the in-medium potential are constrained by results from thermal lattice-QCD for the equation of state (EoS) and heavy-quark (HQ) free energy in a self-consistent set-up for heavy- and light-parton spectral functions in the QGP. We then deploy the refined in-medium heavy-light T-matrix to compute the charm-quark transport coefficients in the QGP. The vector component of the confining potential, through its relativistic corrections, enhances the friction coefficient for charm quarks in the QGP over previous calculations by tens of percent at low momenta and temperatures, and more at higher momenta. Our results are promising for improving the current phenomenology of open heavy-flavor observables at Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Symmetry reduction of tensor networks in many-body theory: I. Automated symbolic evaluation of SU(2) algebra

Abstract The ongoing progress in (nuclear) many-body theory is accompanied by an ever-rising increase in complexity of the underlying formalisms used to solve the stationary Schrödinger equation. The associated working equations at play in state-of-the-art ab initio nuclear many-body methods can be analytically reduced with respect to angular-momentum, i.e. SU (2), quantum numbers whenever they are effectively employed in a symmetry-restricted context. The corresponding procedure constitutes a tedious and error-prone but yet an integral part of the implementation of those many-body frameworks. Indeed, this symmetry reduction is a key step to advance modern simulations to higher accuracy since the use of symmetry-adapted tensors can decrease the computational complexity by orders of magnitude. While attempts have been made in the past to automate the (anti-) commutation rules linked to Fermionic and Bosonic algebras at play in the derivation of the working equations, there is no systematic account to achieve the same goal for their symmetry reduction. In this work, the first version of an automated tool performing graph-theory-based angular-momentum reduction is presented. Taking the symmetry-unrestricted expressions of a generic tensor network as an input, the code provides their angular-momentum-reduced form in an error-safe way in a matter of seconds. Several state-of-the-art many-body methods serve as examples to demonstrate the generality of the approach and to highlight the potential impact on the many-body community.

Tichai, A.↗
DOE OSTI2023

Decoding the flow evolution in Au+Au reactions at 1.23 A GeV using hadron flow correlations and dileptons

We investigate the development of the directed, v 1 , and elliptic flow, v 2 , in heavy ion collisions in mid-central Au+Au reactions at E lab = 1.23A GeV. We demonstrate that the elliptic flow of hot and dense matter is initially positive (v 2 >0) due to the early pressure gradient. This positive v 2 transfers its momentum to the spectators, which leads to the creation of the directed flow v 1 . In turn, the spectator shadowing of the in-plane expansion leads to a preferred decoupling of hadrons in the out-of-plane direction and results in a negative v 2 for the observable final state hadrons. We propose a measurement of v 1 −v 2 flow correlations and of the elliptic flow of dileptons as methods to pin down this evolution pattern. The elliptic flow of the dileptons allows then to determine the early-state EoS more precisely, because it avoids the strong modifications of the momentum distribution due to shadowing seen in the protons. This opens the unique opportunity for the HADES and CBM collaborations to measure the Equation-of-State directly at 2-3 times nuclear saturation density.

Reichert, Tom↗
DOE OSTI2023

Fluid-phase helium: Shock-compression experiments, quantum molecular dynamics simulations, and development of an equation of state

Helium (He) plays a critical role in numerous areas ranging from the study of celestial objects like brown dwarfs and gas giants to modern-day technologies like nuclear energy and rocket propulsion. For many of these applications, it is essential to have a reliable equation of state (EOS) for He that yields an accurate representation of its thermodynamic behavior. To help constrain and develop such EOS models, we have performed a series of shock-compression experiments on cryogenic liquid He to pressures exceeding 100 GPa using a magnetically accelerated flyer plate on Sandia National Laboratories' Z-machine. We have also performed quantum molecular dynamics simulations that are consistent with our shock measurements. None of the previously available EOSs agree with our experimental and simulation results, motivating the development of a fluid-phase He EOS that we present in this study. Here, we show that our EOS yields good agreement with published data that span temperatures and pressures encountered across a diverse array of applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Four-dimensional QCD equation of state with multiple chemical potentials

Here, we construct a four-dimensional version of the equation of state (EoS) model neos, neos-4d, as a function of the temperature and chemical potentials of baryon, electric charge, and strangeness for the hot and dense quantum chromodynamics (QCD) matter created in relativistic nuclear collisions. This EoS enables multiple conserved charge current evolution in a relativistic fluid. Input from lattice QCD simulations and a hadron resonance gas model is considered for constructing the equation of state. We investigate its applicability to the relativistic hydrodynamic description of nuclear collisions and present a method for efficient numerical implementation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗