SEARCH · Search NASA
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.
Symmetry energy of infinite nuclear matter calculated for cases of constant density and minimum saturation energy by density variation, considering finite compressibility effects
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.
A review is made of constraints on the nuclear symmetry energy parameters arising from nuclear binding energy measurements, theoretical chiral effective field predictions of neutron matter properties, the unitary gas conjecture, and measurements of neutron skin thicknesses and dipole polarizabilities. While most studies have been confined to the parameters 𝑆 𝑉 and L, the important roles played by, and constraints on 𝐾 sym , or, equivalently, the neutron matter incompressibility 𝐾 𝑁 , are discussed. Strong correlations among 𝑆 𝑉 ,𝐿, and 𝐾 𝑁 are found from both nuclear binding energies and neutron matter theory. However, these correlations somewhat differ in the two cases, and those from neutron matter theory have smaller uncertainties. To 68% confidence, it is found from neutron matter theory that 𝑆 𝑉 = 32.0 ± 1.1 MeV, 𝐿 = 51.9 ± 7.9 MeV and 𝐾 𝑁 = 152.2 ± 38.1 MeV. Theoretical predictions for neutron skin thickness and dipole polarizability measurements of the neutron-rich nuclei 48 Ca, 120 Sn, and 208 Pb are compared to recent experimental measurements, most notably the CREX and PREX neutron skin experiments from Jefferson Laboratory. By themselves, PREX I+II measurements of 208 Pb and CREX measurement of 48 Ca suggest 𝐿 = 121 ± 47 MeV and 𝐿 = −5 ± 40 MeV, respectively, to 68% confidence. However, we show that nuclear interactions optimally satisfying both measurements imply 𝐿 = 53 ± 13 MeV, nearly the range suggested by either nuclear mass measurements or neutron matter theory, and is also consistent with nuclear dipole polarizability measurements. This small parameter range implies 𝑅 1.4 = 11.6 ± 1.0 km and Λ 1.4 = 228$^{+148}_{−90}$, which are consistent with NICER X-ray and LIGO/Virgo gravitational wave observations of neutron stars.
The nuclear symmetry energy is a key quantity in nuclear (astro)physics. It describes the isospin dependence of the nuclear equation of state, which is commonly assumed to be almost quadratic. Here, in this work, we confront this standard quadratic expansion of the equation of state with explicit asymmetric nuclear-matter calculations based on a set of commonly used Hamiltonians including two- and three-nucleon forces derived from chiral effective-field theory. We study, in particular, the importance of nonquadratic contributions to the symmetry energy, including the nonanalytic logarithmic term introduced by Kaiser [Phys. Rev. C 91, 065201 (2015)]. Our results suggest that the nonquadratic contribution to the symmetry energy can be systematically determined from the various Hamiltonians employed, and we obtain 0.74$^{+0.11}_{−0.08}$ MeV (or −1.02$^{+0.11}_{−0.08}$ MeV for the potential term with the effective-mass contribution) at nuclear saturation density, while the logarithmic contribution to the symmetry energy is relatively small and model-dependent. We also employ the meta-model approach to study the impact of the higher-order contributions on the neutron-star crust-core transition density, and find a 5% correction.
The density dependence of nuclear symmetry energy E sym (ρ) remains the most uncertain aspect of the equation of state (EOS) of supradense neutron-rich nucleonic matter. Utilizing an isospin-dependent parameterization of the nuclear EOS, we investigate the implications of the observational crustal fraction of the neutron star (NS) moment of inertia ΔI/I for the E sym (ρ). We find that symmetry energy parameters significantly influence the ΔI/I, while the EOS of symmetric nuclear matter has a negligible effect. In particular, an increase in the slope L and skewness J sym of symmetry energy results in a larger ΔI/I, whereas an increase in the curvature K sym leads to a reduction in ΔI/I. Moreover, the ΔI/I is shown to have the potential for setting a lower limit of symmetry energy at densities exceeding 3 ρ 0 , particularly when L is constrained to values less than 60 MeV, thereby enhancing our understanding of supradense NS matter.
Background: The nuclear symmetry energy E sym (ρ) encodes information about the energy necessary to make nuclear systems more neutron-rich. While its slope parameter L at the saturation density ρ 0 of nuclear matter has been relatively well constrained by recent astrophysical observations and terrestrial nuclear experiments, its curvature K sym characterizing the E sym (ρ) around 2ρ 0 remains largely unconstrained. Over 520 calculations for E sym (ρ) using various nuclear theories and interactions in the literature have predicted several significantly different K sym –L correlations. Purpose: If a unique K sym –L correlation of E sym (ρ) can be firmly established, it will enable us to progressively better constrain the high-density behavior of E sym (ρ) using the available constraints on its slope parameter L. Here, we investigate if and by how much the different K sym –L correlations may affect neutron star observables. Method: A meta-model of nuclear Equation of States (EOSs) with three representative K sym –L correlation functions is used to generate multiple EOSs for neutron stars. We then examine effects of the K sym –L correlation on the crust-core transition density and pressure as well as the radius and tidal deformation of canonical neutron stars. Results: The K sym –L correlation affects significantly both the crust-core transition density and pressure. It also has strong imprints on the radius and tidal deformability of canonical neutron stars especially at small L values. The available data from LIGO/VIRGO and NICER set some useful limits for the slope L but can not distinguish the three representative K sym –L correlations considered.
The neutron skin thickness Δr np in heavy nuclei has been known as one of the most sensitive terrestrial probes of the nuclear symmetry energy E sym (ρ) around $\frac{2}{3}$ of the saturation density ρ 0 of nuclear matter. Existing neutron skin data mostly from hadronic observables suffer from large uncertainties and their extraction from experiments are often strongly model dependent. While waiting eagerly for the promised model-independent and high-precision neutron skin data for 208 Pb and 48 Ca from the parity-violating electron scattering experiments (PREX-II and CREX at JLab as well as MREX at MESA), within the Bayesian statistical framework using the Skyrme-Hartree-Fock model we infer the posterior probability distribution functions (PDFs) of the slope parameter L of the nuclear symmetry energy at ρ 0 from imagined Δr np ( 208 Pb)=0.15, 0.20, and 0.30 fm with a 1σ error bar of 0.02, 0.04, and 0.06 fm, respectively, as well as Δr np ( 48 Ca)=0.12, 0.15, and 0.25 fm, with different 1σ error bar of 0.01 and 0.02 fm, respectively. The results are compared with the PDFs of L inferred using the same approach from the available Δr np data for 116, 118, 120, 122, 124, 130, 132 Sn from hadronic probes. They are also compared with results from a recent Bayesian analysis of the radius and tidal deformability data of canonical neutron stars from GW170817 and NICER. The neutron skin data for Sn isotopes gives L = 45.5 $^{+ 26.5}_{-21.6}$ MeV surrounding its mean value or L = 53 . 4 $^{+ 18.6}_{ -29.5}$ MeV surrounding its maximum a posteriori value, respectively, with the latter smaller than but consistent with the L = 66 $^{+ 12}_{-20}$ MeV from the neutron star data within their 68% confidence intervals. We found that Δr np = 0.17 –0.18 fm in 208 Pb with an error bar of about 0.02 fm leads to a PDF of L compatible with that from analyzing the Sn data. To provide additionally useful information on L extracted from the Δr np of Sn isotopes, the experimental error bar of Δr np in 208 Pb should be at least smaller than 0.06 fm aimed by some current experiments. In addition, the Δr np ( 48 Ca) needs to be larger than 0.15 fm but smaller than 0.25 fm to be compatible with the Sn and/or neutron star results. To further improve our current knowledge about L and distinguish its PDFs in the examples considered, even higher precisions leading to significantly less than ±20 MeV error bars for L at 68% confidence level are necessary.
By directly inverting several neutron star observables in the three-dimensional parameter space for the Equation of State of super-dense neutron-rich nuclear matter, we show that the lower radius limit for PSR J0740 + 6620 of mass 2.08 ± 0.07 M ⊙ from Neutron Star Interior Composition Explorer (NICER)'s very recent observation sets a much tighter lower boundary than previously known for nuclear symmetry energy in the density range of (1.0 ~ 3.0) times the saturation density ρ 0 of nuclear matter. Furthermore, the super-soft symmetry energy leading to the formation of proton polarons in this density region of neutron stars is clearly disfavoured by the first radius measurement for the most massive neutron star observed reliably so far.
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.
Recent beam energy scan experiments at the BNL Relativistic Heavy Ion Collider by the STAR Collaboration found that hadronic interactions dominate the collective flow and the proton cumulant ratios are driven by baryon number conservation in a region of high baryon density in √ s NN = 3 GeV Au + Au reactions, indicating that the dense medium formed in such collisions is likely hadronic matter. Within an updated a relativistic transport model with momentum dependent isoscalar and isovector single-nucleon mean-field potentials corresponding to different symmetry energies at suprasaturation densities, the n/p, π – /π + , $K$ $^{0}_{s}$/K + , Σ – /Σ + , and Ξ – /Ξ 0 ratios are studied for central Au + Au collisions at √ s NN = 3 GeV, where the maximum central density reaches about (3.6–4.0) ρ 0 . The doubly strange Ξ – /Ξ 0 ratio is found to have the strongest sensitivity to the variation of high-density nuclear symmetry energy. Furthermore, the Ξ – /Ξ 0 ratio in relativistic heavy-ion reactions at √ s NN~3 GeV may help probe sensitively the poorly known symmetry energy of dense neutron-rich matter critically important for understanding various properties of neutron stars.
While ample evidence for the so-called empirical parabolic law of the Equation of State (EOS) of isospin asymmetric nuclear matter (ANM) has been obtained in many studies within both non-relativistic and relativistic nuclear many-body theories using various interactions, it has been unclear if there is any fundamental physics reason for the small quartic symmetry energy compared to the quadratic one even as the ANM approaches pure neutron matter. Within both relativistic and non-relativistic Free Fermi Gas (FFG) models in coordinate spaces of arbitrary dimension d with and without considering Short-Range Correlations (SRC) as well as non-linear Relativistic Mean Field (RMF) models, we study effects of relativistic kinematics, dimensionality, interactions and SRC on the ratio Ψ(ρ) of quartic over quadratic symmetry energies in ANM EOSs. We found that the ratio Ψ(ρ) in the FFG model depends strongly on the dimension d. While it is very small already in the normal 3D space, it could be even smaller in spaces with reduced dimensions for sub-systems of particles in heavy-ion reactions and/or whole neutron stars due to constraints, collectivities and/or symmetries. Here, we also found that the ratio Ψ(ρ) could theoretically become very large only at the ultra-relativistic limit far above the density reachable in neutron stars. On the other hand, nuclear interaction directly and/or indirectly through SRC-induced high-momentum nucleons affect significantly the density dependence of Ψ(ρ) compared to the relativistic FFG model prediction. The SRC affects significantly not only the kinetic energy of symmetric nuclear matter but also the ratio Ψ(ρ) while the relativistic corrections are found negligible. The results may help better understand the EOS of dense neutron-rich matter.
The rapidity distributions of hydrogen isotopes emitted from central collisions of neutron-rich 132 Sn+ 124 Sn and neutron-deficient 108 Sn+ 112 Sn systems at 270 MeV/nucleon were investigated at RIKEN-RIBF. The data are compared with antisymmetrized molecular dynamics (AMD) calculations and the rapidity distributions can be reproduced after adjusting the in-medium nucleon-nucleon cross sections. The double ratios between the two reaction systems taken for the relative yields of deuteron to proton (d/p) and triton to proton (t/p) are further examined in the midrapidity domain, where the adjustments in the AMD calculations do not affect much on them. The d/p and t/p double ratios at midrapidity agree well with the ratio of the system neutron numbers and its squared value, respectively, and the rapidity dependence of these double ratios is consistent with a picture of partial mixing of colliding nuclei. By comparing with the AMD model which shows a strong symmetry energy dependence of the t/p double ratio, the experimental result in the midrapidity domain favors the calculation with a symmetry-energy slope parameter around L = 46 MeV rather than L = 108 MeV.
Not Available
Nuclear symmetry energy E sym (ρ) at density ρ is normally expanded or simply parameterized as a function of χ = (ρ – ρ 0 )/3ρ 0 in the form of Esym(ρ) ≈ S + Lχ + 2 –1 K sym χ 2 + 6 –1 J sym χ 3 + … using its magnitude S, slope L, curvature K sym and skewness J sym at the saturation density ρ0 of nuclear matter. Much progress has been made in recent years in constraining especially the S and L parameters using various terrestrial experiments and astrophysical observations. Furthermore, such kind of expansions/parameterizations do not converge at supra-saturation densities where χ is not small enough, hindering an accurate determination of high-density E sym (ρ) even if its characteristic parameters at ρ0 are all well determined by experiments/observations. By expanding the E sym (ρ) in terms of a properly chosen auxiliary function Π sym (χ,Θ sym ) with a parameter Θsym fixed accurately by an experimental E sym (ρ r ) value at a reference density ρr, we show that the shortcomings of the χ-expansion can be completely removed or significantly reduced in determining the high-density behavior of E sym (ρ). In particular, using two significantly different auxiliary functions, we show that the new approach effectively incorporates higher χ-order contributions and converges to the same E sym (ρ) much faster than the conventional χ-expansion at densities ≲3ρ 0 . Several quantitative demonstrations using Monte Carlo simulations are given.
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.
Here, we introduce a new framework for quantifying correlated uncertainties of the infinite-matter equation of state derived from chiral effective field theory (𝜒EFT ). Bayesian machine learning via Gaussian processes with physics-based hyperparameters allows us to efficiently quantify and propagate theoretical uncertainties of the equation of state, such as 𝜒EFT truncation errors, to derived quantities. We apply this framework to state-of-the-art many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the 𝜒EFT expansion. This produces the first statistically robust uncertainty estimates for key quantities of neutron stars. We give results up to twice nuclear saturation density for the energy per particle, pressure, and speed of sound of neutron matter, as well as for the nuclear symmetry energy and its derivative. At nuclear saturation density, the predicted symmetry energy and its slope are consistent with experimental constraints.
Nonequilibrium and equilibrium schemes studies of nuclear symmetry energy of finite nuclei, discussing binding energy of nuclei with neutron excess
Mass formulas differing in treatment of nuclear symmetry energy with very neutron-rich nuclei being bound in one formula and unbound in other