Constraints on the density dependence of the symmetry energy

A Statistical Approach to Neutron Stars' Crust-Core Transition Density and Pressure

In this paper, a regression model between neutron star crust-core pressure and the symmetry energy characteristics was estimated using the Akaike information criterion and the adjusted coefficient of determination Radj2. The most probable value of the transition density, which should characterize the crust-core environment of the sought physical neutron star model, was determined based on the obtained regression function. An anti-correlation was found between this transition density and the main characteristic of the symmetry energy, i.e., its slope L.

Constraining the Nuclear Symmetry Energy with Multimessenger Resonant Shattering Flares

Much effort is devoted to measuring the nuclear symmetry energy through neutron star (NS) and nuclear observables. Since matter in the NS core may be nonhadronic, observables like radii and tidal deformability may not provide reliable constraints on properties of nucleonic matter. By performing the first consistent inference using ensembles of core and crust equations of state from astrophysical and nuclear data, we demonstrate that coincident timing of a resonant shattering flare (RSF) and gravitational wave signal during binary NS inspiral probes the crust-core transition region and provides constraints on the symmetry energy comparable to terrestrial nuclear experiments. We show that nuclear masses, RSFs, and measurements of NS radii and tidal deformabilities constrain different density ranges of the equation of state, providing complementary probes.

The Symmetry Energy: Current Status of Ab Initio Predictions vs. Empirical Constraints

Infinite nuclear matter is a suitable laboratory to learn about nuclear forces in many-body systems. In particular, modern theoretical predictions of neutron-rich matter are timely because of recent and planned experiments aimed at constraining the equation of state of isospin-asymmetric matter. For these reasons, we have taken a broad look at the equation of state of neutron-rich matter and the closely related symmetry energy, which is the focal point of this article. Its density dependence is of paramount importance for a number of nuclear and astrophysical systems, ranging from neutron skins to the structure of neutron stars. We review and discuss ab initio predictions in relation to recent empirical constraints. We emphasize and demonstrate that free-space nucleon–nucleon data pose stringent constraints on the density dependence of the neutron matter equation of state, which essentially determines the slope of the symmetry energy at saturation.

Constraining neutron-star matter with microscopic and macroscopic collisions

Interpreting high-energy, astrophysical phenomena, such as supernova explosions or neutron-star collisions, requires a robust understanding of matter at supranuclear densities. However, our knowledge about dense matter explored in the cores of neutron stars remains limited. Fortunately, dense matter is not probed only in astrophysical observations, but also in terrestrial heavy-ion collision experiments. Here we use Bayesian inference to combine data from astrophysical multi-messenger observations of neutron stars1-9 and from heavy-ion collisions of gold nuclei at relativistic energies10,11 with microscopic nuclear theory calculations12-17 to improve our understanding of dense matter. We find that the inclusion of heavy-ion collision data indicates an increase in the pressure in dense matter relative to previous analyses, shifting neutron-star radii towards larger values, consistent with recent observations by the Neutron Star Interior Composition Explorer mission5-8,18. Our findings show that constraints from heavy-ion collision experiments show a remarkable consistency with multi-messenger observations and provide complementary information on nuclear matter at intermediate densities. This work combines nuclear theory, nuclear experiment and astrophysical observations, and shows how joint analyses can shed light on the properties of neutron-rich supranuclear matter over the density range probed in neutron stars.

Ranking Love Numbers for the Neutron Star Equation of State: The Need for Third-Generation Detectors

Gravitational-wave measurements of the tidal deformability in neutron-star binary coalescences can be used to infer the still unknown equation of state (EOS) of dense matter above the nuclear saturation density. By employing a Bayesian-ranking test, we quantify the ability of current and future gravitational-wave observations to discriminate among families of nuclear-physics based EOS which differ in particle content and ab initio microscopic calculations. While the constraining power of GW170817 is limited, we show that even twenty coalescences detected by LIGO-Virgo at design sensitivity are not enough to discriminate between EOS with similar softness but distinct microphysics. However, just a single detection with a third-generation detector such as the Einstein Telescope or Cosmic Explorer will rule out several families of EOS with very strong statistical significance and can discriminate among models which feature similar softness, hence, constraining the properties of nuclear matter to unprecedented levels.

Quantifying the Effect of Initial Fluctuations on Isospin-Sensitive Observables from Heavy-Ion Collisions at Intermediate Energies

Initial fluctuation is one of the ingredients that washes fingerprints of the nuclear symmetry energy on observables in heavy-ion collisions. By artificially using the same initial nuclei in all collision events, the effect of the initial fluctuation on isospin-sensitive observables, e.g., the yield ratio of free neutrons with respect to protons Nn/Np, 3H/3He yield ratio, the yield ratio between charged pions π−/π+, and the elliptic flow ratio or difference between free neutrons and protons v2n/v2p (v2n-v2p), are studied within the ultrarelativistic quantum molecular dynamics (UrQMD) model. In practice, Au + Au collisions with impact parameter b = 5 fm and beam energy Elab = 400 MeV/nucleon are calculated. It is found that the effect of the initialization on the yields of free protons and neutrons is small, while for the yield of pions, the directed and elliptic flows are found to be apparently influenced by the choice of initialization because of the strong memory effects. Regarding the isospin-sensitive observables, the effect of the initialization on Nn/Np and 3H/3He is negligible, while π−/π+ and v2n/v2p (v2n-v2p) display a distinct difference among different initializations. The fingerprints of symmetry energy on π−/π+ and v2n/v2p can be either enhanced or reduced when different initializations are utilized.

Symmetry Energy and the Pauli Exclusion Principle

In this article we present a classical potential that respects the Pauli exclusion principle and can be used to describe nucleon-nucleon interactions at intermediate energies. The potential depends on the relative momentum of the colliding nucleons and reduces interactions at low momentum transfer mimicking the Pauli exclusion principle. We use the potential with Metropolis Monte Carlo methods and study the formation of finite nuclei and infinite systems. We find good agreement in terms of the binding energies, radii, and internal nucleon distribution of finite nuclei, and the binding energy in nuclear matter and neutron star matter, as well as the formation of nuclear pastas, and the symmetry energy of neutron star matter.

Unified Equation of State for Neutron Stars Based on the Gogny Interaction

The effective Gogny interactions of the D1 family were established by D. Gogny more than forty years ago with the aim to describe simultaneously the mean field and the pairing field corresponding to the nuclear interaction. The most popular Gogny parametrizations, namely D1S, D1N and D1M, describe accurately the ground-state properties of spherical and deformed finite nuclei all across the mass table obtained with Hartree–Fock–Bogoliubov (HFB) calculations. However, these forces produce a rather soft equation of state (EoS) in neutron matter, which leads to predict maximum masses of neutron stars well below the observed value of two solar masses. To remove this limitation, we built new Gogny parametrizations by modifying the density dependence of the symmetry energy predicted by the force in such a way that they can be applied to the neutron star domain and can also reproduce the properties of finite nuclei as good as their predecessors. These new parametrizations allow us to obtain stiffer EoS’s based on the Gogny interactions, which predict maximum masses of neutron stars around two solar masses. Moreover, other global properties of the star, such as the moment of inertia and the tidal deformability, are in harmony with those obtained with other well tested EoSs based on the SLy4 Skyrme force or the Barcelona–Catania–Paris–Madrid (BCPM) energy density functional. Properties of the core-crust transition predicted by these Gogny EoSs are also analyzed. Using these new Gogny forces, the EoS in the inner crust is obtained with the Wigner–Seitz approximation in the Variational Wigner–Kirkwood approach along with the Strutinsky integral method, which allows one to estimate in a perturbative way the proton shell and pairing corrections. For the outer crust, the EoS is determined basically by the nuclear masses, which are taken from the experiments, wherever they are available, or by HFB calculations performed with these new forces if the experimental masses are not known.

A Modern View of the Equation of State in Nuclear and Neutron Star Matter

Background: We analyze several constraints on the nuclear equation of state (EOS) currently available from neutron star (NS) observations and laboratory experiments and study the existence of possible correlations among properties of nuclear matter at saturation density with NS observables. Methods: We use a set of different models that include several phenomenological EOSs based on Skyrme and relativistic mean field models as well as microscopic calculations based on different many-body approaches, i.e., the (Dirac–)Brueckner–Hartree–Fock theories, Quantum Monte Carlo techniques, and the variational method. Results: We find that almost all the models considered are compatible with the laboratory constraints of the nuclear matter properties as well as with the largest NS mass observed up to now, 2.14−0.09+0.10M⊙ for the object PSR J0740+6620, and with the upper limit of the maximum mass of about 2.3–2.5M⊙ deduced from the analysis of the GW170817 NS merger event. Conclusion: Our study shows that whereas no correlation exists between the tidal deformability and the value of the nuclear symmetry energy at saturation for any value of the NS mass, very weak correlations seem to exist with the derivative of the nuclear symmetry energy and with the nuclear incompressibility.

Probing the Neutron Skin with Ultrarelativistic Isobaric Collisions

Particle production in ultrarelativistic heavy ion collisions depends on the details of the nucleon density distributions in the colliding nuclei. We demonstrate that the charged hadron multiplicity distributions in isobaric collisions at ultrarelativistic energies provide a novel approach to determine the poorly known neutron density distributions and thus the neutron skin thickness in finite nuclei, which can in turn put stringent constraints on the nuclear symmetry energy.

How Well Do We Know the Neutron-Matter Equation of State at the Densities Inside Neutron Stars? A Bayesian Approach with Correlated Uncertainties

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.

Constraining the Neutron Star Equation of State Using Multiband Independent Measurements of Radii and Tidal Deformabilities

Using a Bayesian approach, we combine measurements of neutron star macroscopic observables obtained by astrophysical and gravitational observations to derive joint constraints on the equation of state (EOS) of matter at supranuclear density. In our analysis, we use two sets of data: (i) the masses and tidal deformabilities measured in the binary neutron star event GW170817, detected by LIGO and Virgo, and (ii) the masses and stellar radii measured from observations of nuclear bursts in accreting low-mass x-ray binaries. Using a phenomenological parametrization of the equation of state, we compute the posterior probability distributions of the EOS parameters, using which we infer the posterior distribution for the radius and the mass of the two neutron stars of GW170817. The constraints we set on the radii are tighter than previous bounds.

Novel halos in light kaonic nuclei as an indicator of nuclear equation of state at supra-normal densities

The sensitive correlations between the low-density halo structure and the high-density properties of the nuclear equation of state (EOS) are constructed in light kaonic nuclei with the relativistic mean-field theory. More specifically, the 1p _1/2 halo spreads out linearly with increasing the pressure and sound velocity square at supra-normal densities and decreasing the incompressibility at saturation density. These results suggest that the novel halo in light kaonic nuclei can serve as a sensitive indicator of the nuclear EOS of symmetric matter at supra-normal densities. The experimental production and detection of the light kaonic nuclei, yet to be available, is discussed in some details at last.

Dynamic Isovector Reorientation of Deuteron as a Probe to Nuclear Symmetry Energy

We present the calculations on a novel reorientation effect of deuteron attributed to isovector interaction in the nuclear field of heavy target nuclei. The correlation angle determined by the relative momentum vector of the proton and the neutron originating from the breakup deuteron, which is experimentally detectable, exhibits significant dependence on the isovector nuclear potential but is robust against the variation of the isoscaler sector. In terms of sensitivity and cleanness, the breakup reactions induced by the polarized deuteron beam at about 100 MeV/u provide a more stringent constraint to the symmetry energy at subsaturation densities.

Neutron skin of (208)pb from coherent pion photoproduction

Information on the size and shape of the neutron skin on (208)Pb is extracted from coherent pion photoproduction cross sections measured using the Crystal Ball detector together with the Glasgow tagger at the MAMI electron beam facility. On exploitation of an interpolated fit of a theoretical model to the measured cross sections, the half-height radius and diffuseness of the neutron distribution are found to be c(n)=6.70±0.03(stat.)  fm and a(n)=0.55±0.01(stat.)(-0.03)(+0.02)(sys.)  fm, respectively, corresponding to a neutron skin thickness Δr(np)=0.15±0.03(stat.)(-0.03)(+0.01)(sys.)  fm. The results give the first successful extraction of a neutron skin thickness with an electromagnetic probe and indicate that the skin of (208)Pb has a halo character. The measurement provides valuable new constraints on both the structure of nuclei and the equation of state for neutron-rich matter.

Constraints on the Skyrme equations of state from properties of doubly magic nuclei

I use properties of doubly magic nuclei to constrain symmetric nuclear matter and neutron matter equations of state. I conclude that these data determine the value of the neutron equation of state at a density of ρ(on)=0.10  nucleons/fm3 to be 11.4(10) MeV. The slope at that point is constrained by the value of the neutron skin. Analytical equations are given that show the dependence of the Skyrme equations of state on the neutron skin.

Has a thick neutron skin in 208Pb been ruled out?

The Lead Radius Experiment has provided the first model-independent evidence in favor of a neutron-rich skin in 208Pb. Although the error bars are large, the reported large central value of 0.33 fm is particularly intriguing. To test whether such a thick neutron skin in 208Pb is already incompatible with laboratory experiments or astrophysical observations, we employ relativistic models with neutron-skin thickness in 208Pb ranging from 0.16 to 0.33 fm to compute ground-state properties of finite nuclei, their collective monopole and dipole response, and mass-versus-radius relations for neutron stars. No compelling reason was found to rule out models with large neutron skins in 208Pb from the set of observables considered in this Letter.

Probing the equation of state of nuclear matter via neutron star asteroseismology

We general-relativistically calculate the frequency of fundamental torsional oscillations of neutron star crusts, where we focus on the crystalline properties obtained from macroscopic nuclear models in a way that is dependent on the equation of state of nuclear matter. We find that the calculated frequency is sensitive to the density dependence of the symmetry energy, but almost independent of the incompressibility of symmetric nuclear matter. By identifying the lowest-frequency quasiperiodic oscillation in giant flares observed from soft gamma-ray repeaters as the fundamental torsional mode and allowing for the dependence of the calculated frequency on stellar models, we provide a lower limit of the density derivative of the symmetry energy as L≃50  MeV.

Measurement of the neutron radius of 208Pb through parity violation in electron scattering

We report the first measurement of the parity-violating asymmetry A(PV) in the elastic scattering of polarized electrons from 208Pb. A(PV) is sensitive to the radius of the neutron distribution (R(n)). The result A(PV)=0.656±0.060(stat)±0.014(syst) ppm corresponds to a difference between the radii of the neutron and proton distributions R(n)-R(p)=0.33(-0.18)(+0.16) fm and provides the first electroweak observation of the neutron skin which is expected in a heavy, neutron-rich nucleus.

Connecting neutron star observations to three-body forces in neutron matter and to the nuclear symmetry energy

Using a phenomenological form of the equation of state of neutron matter near the saturation density which has been previously demonstrated to be a good characterization of quantum Monte Carlo simulations, we show that currently available neutron star mass and radius measurements provide a significant constraint on the equation of state of neutron matter. At higher densities we model the equation of state by using polytropes and a quark matter model. We show that observations offer an important constraint on the strength of the three-body force in neutron matter, and thus some theoretical models of the three-body force may be ruled out by currently available astrophysical data. In addition, we obtain an estimate of the symmetry energy of nuclear matter and its slope that can be directly compared to the experiment and other theoretical calculations.

Properties of the nuclear medium

We review our knowledge on the properties of the nuclear medium that have been studied, over many years, on the basis of many-body theory, laboratory experiments and astrophysical observations. Throughout the presentation particular emphasis is placed on the possible relationship and links between the nuclear medium and the structure of nuclei, including the limitations of such an approach. First we consider the realm of phenomenological laboratory data and astrophysical observations and the hints they can give on the characteristics that the nuclear medium should possess. The analysis is based on phenomenological models, that however have a strong basis on physical intuition and an impressive success. More microscopic models are also considered, and it is shown that they are able to give invaluable information on the nuclear medium, in particular on its equation of state. The interplay between laboratory experiments and astrophysical observations is particularly stressed, and it is shown how their complementarity enormously enriches our insights into the structure of the nuclear medium. We then introduce the nucleon-nucleon interaction and the microscopic many-body theory of nuclear matter, with a critical discussion about the different approaches and their results. The Landau-Fermi liquid theory is introduced and briefly discussed, and it is shown how fruitful it can be in discussing the macroscopic and low-energy properties of the nuclear medium. As an illustrative example, we discuss neutron matter at very low density, and it is shown how it can be treated within the many-body theory. The general bulk properties of the nuclear medium are reviewed to indicate at which stage of our knowledge we stand, taking into account the most recent developments both in theory and experiments. A section is dedicated to the pairing problem. The connection with nuclear structure is then discussed, on the basis of the energy density functional method. The possibility of linking the physics of exotic nuclei and the astrophysics of neutron stars is particularly stressed. Finally, we discuss the thermal properties of the nuclear medium, in particular the liquid-gas phase transition and its connection with the phenomenology on heavy ion reactions and the cooling evolution of neutron stars. The presentation has been taken for non-specialists and possibly for non-nuclear physicists.

Neutron skin of (208)Pb, nuclear symmetry energy, and the parity radius experiment

A precise determination of the neutron skin Δr(np) of a heavy nucleus sets a basic constraint on the nuclear symmetry energy (Δr(np) is the difference of the neutron and proton rms radii of the nucleus). The parity radius experiment (PREX) may achieve it by electroweak parity-violating electron scattering (PVES) on (208)Pb. We investigate PVES in nuclear mean field approach to allow the accurate extraction of Δr(np) of (208)Pb from the parity-violating asymmetry A(PV) probed in the experiment. We demonstrate a high linear correlation between A(PV) and Δr(np) in successful mean field forces as the best means to constrain the neutron skin of (208)Pb from PREX, without assumptions on the neutron density shape. Continuation of the experiment with higher precision in A(PV) is motivated since the present method can support it to constrain the density slope of the nuclear symmetry energy to new accuracy.

Suppression of charge equilibration leading to the synthesis of exotic nuclei

Charge equilibration between two colliding nuclei can take place in the early stage of heavy-ion collisions. A basic mechanism of charge equilibration is presented in terms of the extension of single-particle motion from one nucleus to the other, from which the upper energy limit of the bombarding energy is introduced for significant charge equilibration. The formula for this limit is presented, and is compared to various experimental data. It is examined also by comparison to three-dimensional time-dependent density functional calculations. The suppression of charge equilibration, which appears in collisions at the energies beyond the upper energy limit, gives rise to remarkable effects on the synthesis of exotic nuclei with extreme proton-neutron asymmetry.

Symmetry energy of dilute warm nuclear matter

The symmetry energy of nuclear matter is a fundamental ingredient in the investigation of exotic nuclei, heavy-ion collisions, and astrophysical phenomena. New data from heavy-ion collisions can be used to extract the free symmetry energy and the internal symmetry energy at subsaturation densities and temperatures below 10 MeV. Conventional theoretical calculations of the symmetry energy based on mean-field approaches fail to give the correct low-temperature, low-density limit that is governed by correlations, in particular, by the appearance of bound states. A recently developed quantum-statistical approach that takes the formation of clusters into account predicts symmetry energies that are in very good agreement with the experimental data. A consistent description of the symmetry energy is given that joins the correct low-density limit with quasiparticle approaches valid near the saturation density.