Accurate Determination of the Neutron Skin Thickness of ^{208}Pb through Parity-Violation in Electron Scattering

An updated nuclear-physics and multi-messenger astrophysics framework for binary neutron star mergers

The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a 1.4M⊙ neutron star to be [Formula: see text] km.

Determination of the Neutron Skin of ^{208}Pb from Ultrarelativistic Nuclear Collisions

Emergent bulk properties of matter governed by the strong nuclear force give rise to physical phenomena across vastly different scales, ranging from the shape of atomic nuclei to the masses and radii of neutron stars. They can be accessed on Earth by measuring the spatial extent of the outer skin made of neutrons that characterizes the surface of heavy nuclei. The isotope ^{208}Pb, owing to its simple structure and neutron excess, has been in this context the target of many dedicated efforts. Here, we determine the neutron skin from measurements of particle distributions and their collective flow in ^{208}Pb+^{208}Pb collisions at ultrarelativistic energy performed at the Large Hadron Collider, which are mediated by interactions of gluons and thus sensitive to the overall size of the colliding ^{208}Pb ions. By means of state-of-the-art global analysis tools within the hydrodynamic model of heavy-ion collisions, we infer a neutron skin Δr_{np}=0.217±0.058  fm, consistent with nuclear theory predictions, and competitive in accuracy with a recent determination from parity-violating asymmetries in polarized electron scattering. We establish thus a new experimental method to systematically measure neutron distributions in the ground state of atomic nuclei.

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.

Energy-Density Modeling of Strongly Interacting Matter: Atomic Nuclei and Dense Stars

We seek a simple but physically motivated model of strongly interacting matter applicable in atomic nuclei and the dense matter in the core of neutron stars. For densities below and somewhat above normal nuclear density, energy density functional (EDF) theory based on nucleonic degrees of freedom is the ideal candidate. We have explored that direction within the KIDS (Korea-IBS-Daegu-SKKU) framework, which we review in this contribution. The formalism for the KIDS-EoS and microscopic KIDS-EDF and optimization options for the EDF are described in a practical way to facilitate further applications. At densities higher than one nucleon per single-nucleon volume, i.e., roughly 0.4 fm−3, nucleonic degrees of freedom are no longer appropriate. The pseudo-conformal symmetry emergent in dense, topologically altered nuclear matter provides a simple expression for the energy per baryon in terms of the baryonic density. Besides resembling a simple EDF for dense matter, the expression has the appeal that it predicts a converged speed of sound at high densities. It can thus be implemented as a special case of the constant speed of sound (CSS) model. Here we consider a matching between representative nucleonic KIDS-EoSs and the CSS model, including the pseudo-conformal EoS, and apply the unified model to describe the mass–radius relation of neutron stars and examine the compatibility of CSS cores with heavy neutron stars. Although an abrupt transition to the pseudo-conformal regime at low densities does not favor heavy neutron stars, intermediate scenarios including a cusp in the speed of sound are not ruled out, while some appear more favorable to heavy stars than purely nucleonic matter.

Trends of Neutron Skins and Radii of Mirror Nuclei from First Principles

The neutron skin of atomic nuclei impacts the structure of neutron-rich nuclei, the equation of state of nucleonic matter, and the size of neutron stars. Here we predict the neutron skin of selected light- and medium-mass nuclei using coupled-cluster theory and the auxiliary field diffusion Monte Carlo method with two- and three-nucleon forces from chiral effective field theory. We find a linear correlation between the neutron skin and the isospin asymmetry in agreement with the liquid-drop model and compare with data. We also extract the linear relationship that describes the difference between neutron and proton radii of mirror nuclei and quantify the effect of charge symmetry breaking terms in the nuclear Hamiltonian. Our results for the mirror-difference charge radii and binding energies per nucleon agree with existing data.

Tomography of ultrarelativistic nuclei with polarized photon-gluon collisions

A linearly polarized photon can be quantized from the Lorentz-boosted electromagnetic field of a nucleus traveling at ultrarelativistic speed. When two relativistic heavy nuclei pass one another at a distance of a few nuclear radii, the photon from one nucleus may interact through a virtual quark-antiquark pair with gluons from the other nucleus, forming a short-lived vector meson (e.g., ρ⁰). In this experiment, the polarization was used in diffractive photoproduction to observe a unique spin interference pattern in the angular distribution of ρ⁰ → π⁺π^- decays. The observed interference is a result of an overlap of two wave functions at a distance an order of magnitude larger than the ρ⁰ travel distance within its lifetime. The strong-interaction nuclear radii were extracted from these diffractive interactions and found to be 6.53 ± 0.06 fm (¹⁹⁷Au) and 7.29 ± 0.08 fm (²³⁸U), larger than the nuclear charge radii. The observable is demonstrated to be sensitive to the nuclear geometry and quantum interference of nonidentical particles.

Probing neutron skin and symmetry energy with relativistic isobar collisions

In these proceedings, we present the three proposed observables to probe the neutron skin and symmetry energy with relativistic isobar collisions, namely, the isobar ratios of the produced hadron multiplicities (Nch), the mean transverse momenta (〈p⊥〉), and the net charge multiplicities (ΔQ). Our findings suggest potentially significant improvement to neutron skin and symmetry energy determination over traditional low energy methods.

Extended p_{3/2} Neutron Orbital and the N=32 Shell Closure in ^{52}Ca

The one-neutron knockout from ^{52}Ca in inverse kinematics onto a proton target was performed at ∼230  MeV/nucleon combined with prompt γ spectroscopy. Exclusive quasifree scattering cross sections to bound states in ^{51}Ca and the momentum distributions corresponding to the removal of 1f_{7/2} and 2p_{3/2} neutrons were measured. The cross sections, interpreted within the distorted-wave impulse approximation reaction framework, are consistent with a shell closure at the neutron number N=32, found as strong as at N=28 and N=34 in Ca isotopes from the same observables. The analysis of the momentum distributions leads to a difference of the root-mean-square radii of the neutron 1f_{7/2} and 2p_{3/2} orbitals of 0.61(23) fm, in agreement with the modified-shell-model prediction of 0.7 fm suggesting that the large root-mean-square radius of the 2p_{3/2} orbital in neutron-rich Ca isotopes is responsible for the unexpected linear increase of the charge radius with the neutron number.

Combined Theoretical Analysis of the Parity-Violating Asymmetry for ^{48}Ca and ^{208}Pb

The recent experimental determination of the parity violating asymmetry A_{PV} in ^{48}Ca and ^{208}Pb at Jefferson Lab is important for our understanding on how neutrons and protons arrange themselves inside the atomic nucleus. To better understand the impact of these measurements, we present a rigorous theoretical investigation of A_{PV} in ^{48}Ca and ^{208}Pb and assess the associated uncertainties. We complement our study by inspecting the static electric dipole polarizability in these nuclei. The analysis is carried out within nuclear energy density functional theory with quantified input. We conclude that the simultaneous accurate description of A_{PV} in ^{48}Ca and ^{208}Pb cannot be achieved by our models that accommodate a pool of global nuclear properties, such as masses and charge radii, throughout the nuclear chart, and describe-within one standard deviation-the experimental dipole polarizabilities α_{D} in these nuclei.

Ab initio predictions link the neutron skin of 208Pb to nuclear forces

Heavy atomic nuclei have an excess of neutrons over protons, which leads to the formation of a neutron skin whose thickness is sensitive to details of the nuclear force. This links atomic nuclei to properties of neutron stars, thereby relating objects that differ in size by orders of magnitude. The nucleus 208Pb is of particular interest because it exhibits a simple structure and is experimentally accessible. However, computing such a heavy nucleus has been out of reach for ab initio theory. By combining advances in quantum many-body methods, statistical tools and emulator technology, we make quantitative predictions for the properties of 208Pb starting from nuclear forces that are consistent with symmetries of low-energy quantum chromodynamics. We explore 109 different nuclear force parameterizations via history matching, confront them with data in select light nuclei and arrive at an importance-weighted ensemble of interactions. We accurately reproduce bulk properties of 208Pb and determine the neutron skin thickness, which is smaller and more precise than a recent extraction from parity-violating electron scattering but in agreement with other experimental probes. This work demonstrates how realistic two- and three-nucleon forces act in a heavy nucleus and allows us to make quantitative predictions across the nuclear landscape.

Precision Determination of the Neutral Weak Form Factor of ^{48}Ca

We report a precise measurement of the parity-violating (PV) asymmetry A_{PV} in the elastic scattering of longitudinally polarized electrons from ^{48}Ca. We measure A_{PV}=2668±106(stat)±40(syst) parts per billion, leading to an extraction of the neutral weak form factor F_{W}(q=0.8733  fm^{-1})=0.1304±0.0052(stat)±0.0020(syst) and the charge minus the weak form factor F_{ch}-F_{W}=0.0277±0.0055. The resulting neutron skin thickness R_{n}-R_{p}=0.121±0.026(exp)±0.024(model)  fm is relatively thin yet consistent with many model calculations. The combined CREX and PREX results will have implications for future energy density functional calculations and on the density dependence of the symmetry energy of nuclear matter.

Sensitivity of Parity-Violating Electron Scattering to a Dark Photon

We explore the sensitivity of the parity-violating electron scattering (PVES) asymmetry in both elastic and deep-inelastic scattering to the properties of a dark photon. Given advances in experimental capabilities in recent years, there are interesting regions of parameter space where PVES offers the chance to discover new physics in the near future. There are also cases where the existence of a dark photon could significantly alter our understanding of the structure of atomic nuclei and neutron stars as well as parton distribution functions.

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.

New Measurements of the Beam-Normal Single Spin Asymmetry in Elastic Electron Scattering over a Range of Spin-0 Nuclei

We report precision determinations of the beam-normal single spin asymmetries (A_{n}) in the elastic scattering of 0.95 and 2.18 GeV electrons off ^{12}C, ^{40}Ca, ^{48}Ca, and ^{208}Pb at very forward angles where the most detailed theoretical calculations have been performed. The first measurements of A_{n} for ^{40}Ca and ^{48}Ca are found to be similar to that of ^{12}C, consistent with expectations and thus demonstrating the validity of theoretical calculations for nuclei with Z≤20. We also report A_{n} for ^{208}Pb at two new momentum transfers (Q^{2}) extending the previous measurement. Our new data confirm the surprising result previously reported, with all three data points showing significant disagreement with the results from the Z≤20 nuclei. These data confirm our basic understanding of the underlying dynamics that govern A_{n} for nuclei containing ≲50 nucleons, but point to the need for further investigation to understand the unusual A_{n} behavior discovered for scattering off ^{208}Pb.

Determination of the ^{27}Al Neutron Distribution Radius from a Parity-Violating Electron Scattering Measurement

We report the first measurement of the parity-violating elastic electron scattering asymmetry on ^{27}Al. The ^{27}Al elastic asymmetry is A_{PV}=2.16±0.11(stat)±0.16(syst)  ppm, and was measured at ⟨Q^{2}⟩=0.02357±0.00010  GeV^{2}, ⟨θ_{lab}⟩=7.61°±0.02°, and ⟨E_{lab}⟩=1.157  GeV with the Q_{weak} apparatus at Jefferson Lab. Predictions using a simple Born approximation as well as more sophisticated distorted-wave calculations are in good agreement with this result. From this asymmetry the ^{27}Al neutron radius R_{n}=2.89±0.12  fm was determined using a many-models correlation technique. The corresponding neutron skin thickness R_{n}-R_{p}=-0.04±0.12  fm is small, as expected for a light nucleus with a neutron excess of only 1. This result thus serves as a successful benchmark for electroweak determinations of neutron radii on heavier nuclei. A tree-level approach was used to extract the ^{27}Al weak radius R_{w}=3.00±0.15  fm, and the weak skin thickness R_{wk}-R_{ch}=-0.04±0.15  fm. The weak form factor at this Q^{2} is F_{wk}=0.39±0.04.

Nuclear Potentials Relevant to the Symmetry Energy in Chiral Models

We employ the extended Nambu-Jona-Lasinio (NJL), linear-σ models, and the density-dependent model with chiral limits to work out the mean fields and relevant properties of nuclear matter. To have the constraint from the data, we re-examine the Dirac optical potentials and symmetry potential based on the relativistic impulse approximation (RIA). Unlike the extended NJL and the density-dependent models with the chiral limit in terms of the vanishing scalar density, the extended linear-σ model with a sluggish changing scalar field loses the chiral limit at the high-density end. The various scalar fields can characterize the different Schrödinger-equivalent potentials and kinetic symmetry energy in the whole density region and the symmetry potential in the intermediate density region. The drop in the scalar field due to the chiral restoration results in a clear rise of the kinetic symmetry energy. The chiral limit in the models gives rise to the softening of the symmetry potential and thereof the symmetry energy at high densities.

Translating Neutron Star Observations to Nuclear Symmetry Energy via Deep Neural Networks

One of the most significant challenges involved in efforts to understand the equation of state of dense neutron-rich matter is the uncertain density dependence of the nuclear symmetry energy. In particular, the nuclear symmetry energy is still rather poorly constrained, especially at high densities. On the other hand, detailed knowledge of the equation of state is critical for our understanding of many important phenomena in the nuclear terrestrial laboratories and the cosmos. Because of its broad impact, pinning down the density dependence of the nuclear symmetry energy has been a long-standing goal of both nuclear physics and astrophysics. Recent observations of neutron stars, in both electromagnetic and gravitational-wave spectra, have already constrained significantly the nuclear symmetry energy at high densities. The next generation of telescopes and gravitational-wave observatories will provide an unprecedented wealth of detailed observations of neutron stars, which will improve further our knowledge of the density dependence of nuclear symmetry energy, and the underlying equation of state of dense neutron-rich matter. Training deep neural networks to learn a computationally efficient representation of the mapping between astrophysical observables of neutron stars, such as masses, radii, and tidal deformabilities, and the nuclear symmetry energy allows its density dependence to be determined reliably and accurately. In this work, we use a deep learning approach to determine the nuclear symmetry energy as a function of density directly from observational neutron star data. We show, for the first time, that artificial neural networks can precisely reconstruct the nuclear symmetry energy from a set of available neutron star observables, such as masses and radii as measured by, e.g., the NICER mission, or masses and tidal deformabilities as measured by the LIGO/VIRGO/KAGRA gravitational-wave detectors. These results demonstrate the potential of artificial neural networks to reconstruct the symmetry energy and the equation of state directly from neutron star observational data, and emphasize the importance of the deep learning approach in the era of multi-messenger astrophysics.

Information Content of the Parity-Violating Asymmetry in ^{208}Pb

The parity-violating asymmetry A_{PV} in ^{208}Pb, recently measured by the PREX-2 Collaboration, is studied using modern relativistic (covariant) and nonrelativistic energy density functionals. We first assess the theoretical uncertainty on A_{PV} which is intrinsic to the adopted approach. To this end, we use quantified functionals that are able to accommodate our previous knowledge on nuclear observables such as binding energies, charge radii, and the dipole polarizability α_{D} of ^{208}Pb. We then add the quantified value of A_{PV} together with α_{D} to our calibration dataset to optimize new functionals. Based on these results, we predict a neutron skin thickness in ^{208}Pb r_{skin}=0.19±0.02  fm and the symmetry-energy slope L=54±8  MeV. These values are consistent with other estimates based on astrophysical data and are significantly lower than those recently reported using a particular set of relativistic energy density functionals. We also make a prediction for the A_{PV} value in ^{48}Ca that will be soon available from the CREX measurement.

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_{skin}^{^{208}Pb} 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. Combining astrophysical data with PREX-II and chiral effective field theory constraints yields S_{0}=33.0_{-1.8}^{+2.0}  MeV, L=53_{-15}^{+14}  MeV, and R_{skin}^{^{208}Pb}=0.17_{-0.04}^{+0.04}  fm.

Charge Radius of Neutron-Deficient ^{54}Ni and Symmetry Energy Constraints Using the Difference in Mirror Pair Charge Radii

The nuclear root-mean-square charge radius of ^{54}Ni was determined with collinear laser spectroscopy to be R(^{54}Ni)=3.737(3)  fm. In conjunction with the known radius of the mirror nucleus ^{54}Fe, the difference of the charge radii was extracted as ΔR_{ch}=0.049(4)  fm. Based on the correlation between ΔR_{ch} and the slope of the symmetry energy at nuclear saturation density (L), we deduced 21≤L≤88  MeV. The present result is consistent with the L from the binary neutron star merger GW170817, favoring a soft neutron matter EOS, and barely consistent with the PREX-2 result within 1σ error bands. Our result indicates the neutron-skin thickness of ^{48}Ca as 0.15-0.21 fm.

The Nuclear Matter Density Functional under the Nucleonic Hypothesis

A Bayesian analysis of the possible behaviors of the dense matter equation of state informed by recent LIGO-Virgo as well as NICER measurements reveals that all the present observations are compatible with a fully nucleonic hypothesis for the composition of dense matter, even in the core of the most massive pulsar PSR J0740+6620. Under the hypothesis of a nucleonic composition, we extract the most general behavior of the energy per particle of symmetric matter and density dependence of the symmetry energy, compatible with the astrophysical observations as well as our present knowledge of low-energy nuclear physics from effective field theory predictions and experimental nuclear mass data. These results can be used as a null hypothesis to be confronted with future constraints on dense matter to search for possible exotic degrees of freedom.

A Neutron Star Is Born

A neutron star was first detected as a pulsar in 1967. It is one of the most mysterious compact objects in the universe, with a radius of the order of 10 km and masses that can reach two solar masses. In fact, neutron stars are star remnants, a kind of stellar zombie (they die, but do not disappear). In the last decades, astronomical observations yielded various contraints for neutron star masses, and finally, in 2017, a gravitational wave was detected (GW170817). Its source was identified as the merger of two neutron stars coming from NGC 4993, a galaxy 140 million light years away from us. The very same event was detected in γ-ray, X-ray, UV, IR, radio frequency and even in the optical region of the electromagnetic spectrum, starting the new era of multi-messenger astronomy. To understand and describe neutron stars, an appropriate equation of state that satisfies bulk nuclear matter properties is necessary. GW170817 detection contributed with extra constraints to determine it. On the other hand, magnetars are the same sort of compact object, but bearing much stronger magnetic fields that can reach up to 1015 G on the surface as compared with the usual 1012 G present in ordinary pulsars. While the description of ordinary pulsars is not completely established, describing magnetars poses extra challenges. In this paper, I give an overview on the history of neutron stars and on the development of nuclear models and show how the description of the tiny world of the nuclear physics can help the understanding of the cosmos, especially of the neutron stars.

Nuclear Physics and Astrophysics Constraints on the High Density Matter Equation of State

(1) This review has been written in memory of Steven Moszkowski who unexpectedly passed away in December 2020. It has been inspired by our many years of discussions. Steven’s enthusiasm, drive and determination to understand atomic nuclei in simple terms of basic laws of physics was infectious. He sought the fundamental origin of nuclear forces in free space, and their saturation and modification in nuclear medium. His untimely departure left our job unfinished but his legacy lives on. (2) Focusing on the nuclear force acting in nuclear matter of astrophysical interest and its equation of state (EoS), we take several typical snapshots of evolution of the theory of nuclear forces. We start from original ideas in the 1930s moving through to its overwhelming diversity today. The development is supported by modern observational and terrestrial data and their inference in the multimessenger era, as well as by novel mathematical techniques and computer power. (3) We find that, despite the admirable effort both in theory and measurement, we are facing multiple models dependent on a large number of variable correlated parameters which cannot be constrained by data, which are not yet accurate, nor sensitive enough, to identify the theory closest to reality. The role of microphysics in the theories is severely limited or neglected, mostly deemed to be too difficult to tackle. (4) Taking the EoS of high-density matter as an example, we propose to develop models, based, as much as currently possible, on the microphysics of the nuclear force, with a minimal set of parameters, chosen under clear physical guidance. Still somewhat phenomenological, such models could pave the way to realistic predictions, not tracing the measurement, but leading it.

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

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

Implications of PREX-2 on the Equation of State of Neutron-Rich Matter

Laboratory experiments sensitive to the equation of state of neutron rich matter in the vicinity of nuclear saturation density provide the first rung in a "density ladder" that connects terrestrial experiments to astronomical observations. In this context, the neutron skin thickness of ^{208}Pb (R_{skin}^{208}) provides a stringent laboratory constraint on the density dependence of the symmetry energy. In turn, an improved value of R_{skin}^{208} has been reported recently by the PREX collaboration. Exploiting the strong correlation between R_{skin}^{208} and the slope of the symmetry energy L within a specific class of relativistic energy density functionals, we report a value of L=(106±37)  MeV-which systematically overestimates current limits based on both theoretical approaches and experimental measurements. The impact of such a stiff symmetry energy on some critical neutron-star observables is also examined.