CONSTRAINING THE SYMMETRY PARAMETERS OF THE NUCLEAR INTERACTION

The nuclear symmetry energy is intimately connected with nuclear astrophysics. This contribution focuses on the estimation of the symmetry energy from experiment and how it is related to the structure of neutron stars. The most important connection is between the radii of neutron stars and the pressure of neutron star matter in the vicinity of the nuclear saturation density ns. This pressure is essentially controlled by the nuclear symmetry energy parameters Sv and L , the first two coefficients of a Taylor expansion of the symmetry energy around ns. We discuss constraints on these parameters that can be found from nuclear experiments. We demonstrate that these constraints are largely model-independent by deriving them qualitatively from a simple nuclear model. We also summarize how recent theoretical studies of pure neutron matter can reinforce these constraints. To date, several different astrophysical measurements of neutron star radii have been attempted. Attention is focused on photospheric radius expansion bursts and on thermal emissions from quiescent low-mass X-ray binaries. While none of these observations can, at the present time, determine individual neutron star radii to better than 20% accuracy, the body of observations can be used with Bayesian techniques to effectively constrain them to higher precision. These techniques invert the structure equations and obtain estimates of the pressure-density relation of neutron star matter, not only near ns, but up to the highest densities found in neutron star interiors. The estimates we derive for neutron star radii are in concordance with predictions from nuclear experiment and theory.

We calculate the equation of state of neutron matter with realistic two- and\nthree-nucleon interactions using quantum Monte Carlo techniques, and illustrate\nthat the short-range three-neutron interaction determines the correlation\nbetween neutron matter energy at nuclear saturation density and higher\ndensities relevant to neutron stars. Our model also makes an experimentally\ntestable prediction for the correlation between the nuclear symmetry energy and\nits density dependence -- determined solely by the strength of the short-range\nterms in the three neutron force. The same force provides a significant\nconstraint on the maximum mass and radius of neutron stars.\n

In this paper, the properties of the symmetry energy in asymmetric nuclear matter, including the quadratic and quartic symmetry energies as well as their density slopes at the saturation density, are studied by the radioactivity of unstable nuclei. The symmetry energy was initially introduced as a quantum correction energy in the Weizsacker-Bethe liquid-drop model, which was proposed to describe the binding energies of finite nuclei. In order to systematically study the properties of the symmetry energy, the concept of the nuclear matter was introduced, which consists of neutrons and protons with finite neutron and proton densities but infinite neutron and proton numbers. By comparing the properties of the isospin symmetric nuclear matter, in which the neutron density is equal to the proton density, and those of the isospin asymmetric nuclear matter, in which the two densities are unequal, the symmetry energy is found to approximately represent the difference of the energy per nucleon in them. That is to say, the symmetry energy is an evaluation of the energy cost to convert protons in symmetric nuclear matter to neutrons in asymmetric nuclear matter. Since the symmetry energy characterizes the isospin asymmetric effect in nuclear matter, it plays critical roles in not only nuclear physics but also astrophysics. For example, the symmetry energy is closely related to the nuclear masses, the structures and properties of nuclei far away from stability line and near drip lines, the mechanism of heavy-ion reactions, the structures and component of neutron stars, the masses and radii of neutron stars, the cooling process of neutron stars, and so on. Besides of the quadratic symmetry energy, the quartic symmetry energy can have obvious effects on the properties of neutron stars as well. The quartic symmetry energy can significantly affect both the proton fraction in β-stable neutron stars and the critical density for the direct Urca process, which leads to faster cooling of neutron stars. In addition, the quartic symmetry energy is found to be very important for the location of the inner edge of crusts and the core-crust transition density and pressure in neutron stars. So the quartic symmetry energy can also have an influence on the structure of neutron stars. Because of the very important roles of the symmetry energy playing in not only nuclear physics but also astrophysics, it is of great significance to study the properties of both the quadratic and quartic symmetry energies. Based on the fundamental Hugenholtz-Van Hove (HVH) theorem, the properties of the symmetry energy are found to be closely related to the isoscalar and the isovector parts of the nucleon single-particle potentials, where the latter two parts can be extracted from the nuclear radioactivity. So by linking both the symmetry energy and the nuclear radioactivity to the nucleon single-particle potentials, the symmetry energy can be directly studied by nuclear radioactivity. In this paper, both the heavy-cluster radioactivity and the proton radioactivity are employed to constrain the properties of not only the quadratic symmetry energy but also the quartic symmetry energy. Besides, the effect of both the spectroscopic factor of proton radioactivity and the deformation in the daughter nucleus on the extracted results are also discussed, which are subsequently found to be very limited.

We investigate constraints on neutron star structure arising from the assumptions that neutron stars have crusts, that recent calculations of pure neutron matter limit the equation of state of neutron star matter near the nuclear saturation density, that the high-density equation of state is limited by causality and the largest high-accuracy neutron star mass measurement, and that general relativity is the correct theory of gravity. We explore the role of prior assumptions by considering two classes of equation of state models. In a first, the intermediate- and high-density behavior of the equation of state is parameterized by piecewise polytropes. In the second class, the high-density behavior of the equation of state is parameterized by piecewise continuous line segments. The smallest density at which high-density matter appears is varied in order to allow for strong phase transitions above the nuclear saturation density. We critically examine correlations among the pressure of matter, radii, maximum masses, the binding energy, the moment of inertia, and the tidal deformability, paying special attention to the sensitivity of these correlations to prior assumptions about the equation of state. It is possible to constrain the radii of $1.4~\mathrm{M}_{\odot}$ neutron stars to a be larger than 10 km, even without consideration of additional astrophysical observations, for example, those from photospheric radius expansion bursts or quiescent low-mass X-ray binaries. We are able to improve the accuracy of known correlations between the moment of inertia and compactness as well as the binding energy and compactness. We also demonstrate the existence of a correlation between the neutron star binding energy and the moment of inertia.

Nuclear symmetry energy effects in finite nuclei and neutron star

The radii and tidal deformabilities of neutron stars are investigated in the framework of the relativistic mean-field (RMF) model with different density-dependent behaviors of symmetry energy. To study the effects of symmetry energy on the properties of neutron stars, $\omega$ meson and $\rho$ meson coupling terms are included in a popular RMF Lagrangian, i.e., the TM1 parameter set, which is adopted for the widely used supernova equation of state (EoS) table. The coupling constants relevant to the vector–isovector meson, $\rho$, are refitted by a fixed symmetry energy at subsaturation density and its slope at saturation density, while other coupling constants remain the same as the original ones in TM1 so as to update the supernova EoS table. The radius and mass of maximum neutron stars are not so sensitive to the symmetry energy in these family TM1 parameterizations. However, the radii in the intermediate-mass region are strongly correlated with the slope of symmetry energy. Furthermore, the dimensionless tidal deformabilities of neutron stars are also calculated within the associated Love number, which is related to the quadrupole deformation of the star in a static external tidal field and can be extracted from the observation of a gravitational wave generated by a binary star merger. We find that its value at $1.4 \mathrm{M}_\odot$ has a linear correlation to the slope of symmetry energy, unlike that previously studied. With the latest constraints of tidal deformabilities from the GW170817 event, the slope of symmetry energy at nuclear saturation density should be smaller than $60$ MeV in the family TM1 parameterizations. This fact supports the usage of a lower symmetry energy slope for the updated supernova EoS, which is applicable to simulations of neutron star mergers. Furthermore, an analogous analysis is also done within the family IUFSU parameter sets. It is found that the correlations between the symmetry energy slope with the radius and tidal deformability at $1.4 \mathrm{M}_\odot$ have very similar linear relations in these RMF models.

The properties of neutron and strange matter stars are discussed from global and observational perspectives. The global features, i.e., the mass-radius relation, the moment of inertia, and the binding energy, of these objects can be understood by examination of the few known relevant analytic solutions to Einstein's equations. A close connection exists between neutron star radii and the density dependence of the isospin dependence of strong interactions, i.e., the nuclear symmetry energy. Interestingly, a similar relation has been found to exist between the symmetry energy and the neutron skin thickness of neutron-rich nuclei, the object of a new generation of laboratory experiments, although these are 1018 times smaller. Recent observations of neutron star masses and radiation radii are summarized. The status of these observations as predictors of nuclear force properties is examined. The combination of observations, laboratory experiments, and theory is an extremely powerful tool for both nuclear physics and nuclear astrophysics.

We propose the existence of a lower bound on the energy of pure neutron matter (PNM) on the basis of unitary-gas considerations. We discuss its justification from experimental studies of cold atoms as well as from theoretical studies of neutron matter. We demonstrate that this bound results in limits to the density-dependent symmetry energy, which is the difference between the energies of symmetric nuclear matter and PNM. In particular, this bound leads to a lower limit to the volume symmetry energy parameter S 0. In addition, for assumed values of S 0 above this minimum, this bound implies both upper and lower limits to the symmetry energy slope parameter L ,which describes the lowest-order density dependence of the symmetry energy. A lower bound on neutron-matter incompressibility is also obtained. These bounds are found to be consistent with both recent calculations of the energies of PNM and constraints from nuclear experiments. Our results are significant because several equations of state that are currently used in astrophysical simulations of supernovae and neutron star mergers, as well as in nuclear physics simulations of heavy-ion collisions, have symmetry energy parameters that violate these bounds. Furthermore, below the nuclear saturation density, the bound on neutron-matter energies leads to a lower limit to the density-dependent symmetry energy, which leads to upper limits to the nuclear surface symmetry parameter and the neutron-star crust–core boundary. We also obtain a lower limit to the neutron-skin thicknesses of neutron-rich nuclei. Above the nuclear saturation density, the bound on neutron-matter energies also leads to an upper limit to the symmetry energy, with implications for neutron-star cooling via the direct Urca process.

Nonuniform structure of low-density nuclear matter, known as nuclear pasta, is expected to appear not only in the inner crust of neutron stars but also in core-collapse supernova explosions and neutron-star mergers. We perform fully three-dimensional calculations of inhomogeneous nuclear matter and neutron-star matter in the low-density region using the Thomas-Fermi approximation. The nuclear interaction is described in the relativistic mean-field approach with the point-coupling interaction, where the meson exchange in each channel is replaced by the contact interaction between nucleons. We investigate the influence of nuclear symmetry energy and its density dependence on pasta structures by introducing a coupling term between the isoscalar-vector and isovector-vector interactions. It is found that the properties of pasta phases in the neutron-rich matter are strongly dependent on the symmetry energy and its slope. In addition to typical shapes like droplets, rods, slabs, tubes, and bubbles, some intermediate pasta structures are also observed in cold stellar matter with a relatively large proton fraction. We find that nonspherical shapes are unlikely to be formed in neutron-star crusts, because the proton fraction obtained in $\ensuremath{\beta}$ equilibrium is somewhat small. The inner crust properties may lead to a visible difference in the neutron-star radius.

Pulse profile modeling of X-ray data from the Neutron Star Interior Composition Explorer is now enabling precision inference of neutron star mass and radius. Combined with nuclear physics constraints from chiral effective field theory (χEFT), and masses and tidal deformabilities inferred from gravitational-wave detections of binary neutron star mergers, this has led to a steady improvement in our understanding of the dense matter equation of state (EOS). Here, we consider the impact of several new results: the radius measurement for the 1.42 M ⊙ pulsar PSR J0437−4715 presented by Choudhury et al., updates to the masses and radii of PSR J0740+6620 and PSR J0030+0451, and new χEFT results for neutron star matter up to 1.5 times nuclear saturation density. Using two different high-density EOS extensions—a piecewise-polytropic (PP) model and a model based on the speed of sound in a neutron star (CS)—we find the radius of a 1.4 M ⊙ (2.0 M ⊙) neutron star to be constrained to the 95% credible ranges 12.28−0.76+0.50 km ( 12.33−1.34+0.70 km) for the PP model and 12.01−0.75+0.56 km ( 11.55−1.09+0.94 km) for the CS model. The maximum neutron star mass is predicted to be 2.15−0.16+0.14 M ⊙ and 2.08−0.16+0.28 M ⊙ for the PP and CS models, respectively. We explore the sensitivity of our results to different orders and different densities up to which χEFT is used, and show how the astrophysical observations provide constraints for the pressure at intermediate densities. Moreover, we investigate the difference R 2.0 − R 1.4 of the radius of 2 M ⊙ and 1.4 M ⊙ neutron stars within our EOS inference.

Neutron stars are born when massive stars run out of their nuclear fuel and undergo gravitational collapse. Neutron stars belong to the most compact objects in the observable Universe. Macroscopic properties of neutron stars like their masses and radii are sensitive to the microscopic properties of the nuclear equation of state of dense matter. The equation of state is determined by the strong interaction among the constituents. The underlying theory is quantum chromodynamics that is, however, highly non-perturbative in the physics regime relevant for neutron stars. Moreover, neutron stars provide an interplay between nuclear physics and astrophysics. Astrophysical observations like the detection of 2 solar mass neutron stars have a major impact on the equation of state. Radii are, however, inherently difficult to measure due to systematic uncertainties. Other observables like the moment of inertia or the tidal deformability present promising alternatives. The double neutron star system PSR J0737-3039 constitutes an outstanding system as it provides the prospect of a moment of inertia measurement for the first time. A new era stated with the pioneering observation of gravitational waves from a binary neutron star merger. The analysis of the gravitational wave signal of GW170817 provides a range for the tidal deformability of typical neutron stars. Moreover, the current NICER mission will provide simultaneous mass-radius measurements. In this thesis, we use state-of-the-art chiral effective field theory interactions to describe the equation of state at nuclear densities. In the high-density regime beyond nuclear saturation density, we use different extrapolation approaches. First, we utilize the established ansatz of piecewise polytropic equations of state which provides a direct parametrization. However, piecewise polytropic equations of state possess unphysical behavior such as discontinuities in the speed of sound. Second, we use a physically motivated parametrization of the speed of sound inside the neutron star from which we derive the equation of state. Both methods allow us to probe the equation of state over a large range of densities. We further impose general constraints on the equation of state such as the requirement of causality at all densities and the support of at least 2 solar mass neutron stars. From the equations of state compatible with the constraints, we determine diverse neutron star observables. We begin with non-rotating neutron stars and focus on their masses and radii. We study correlations among properties of the equation of state at nuclear densities and observables of typical neutron stars. Moreover, we explore the impact of hypothetical, simultaneous measurements of masses and radii of neutron stars on the equation of state. Applying both simple compatibility cuts and the framework of Bayesian statistics, we investigate the sensitivity of the inference on the chosen parametrization of the equation of state. We extend then our considerations to slowly rotating neutron stars and study the moment of inertia. Assuming hypothetical moment of inertia measurements, we determine constraints for the radius of neutron stars and thus the equation of state. In addition, we extend our considerations of isolated neutron stars to binary neutron star systems. In particular, we treat the tidal field of the companion as a small perturbation. This allows us to determine the tidal deformability. By applying higher orders in the metric perturbation, we calculate the quadrupole moment of neutron stars. Although the structure of neutron stars is sensitive to the equation of state, relations between the moment of inertia, the tidal deformability, and the quadrupole moment are remarkably insensitive. We investigate the properties of neutron stars in binary systems and ultimately confront the results of our models with the gravitational wave constraints from a binary neutron star merger.

The density dependence of nuclear symmetry energy remains poorly constrained. Starting from precise empirical values of the nuclear volume and surface symmetry energy coefficients and the nuclear saturation density, we show how in the ambit of microscopic calculations with different energy density functionals, the value of the symmetry energy slope parameter L along with that for neutron skin can be put in tighter bounds. The value of L is found to be L=64±5 MeV. For 208Pb, the neutron skin thickness comes out to be 0.188±0.014 fm. Knowing L, the method can be applied to predict neutron skin thicknesses of other nuclei.

The equation of state (EOS) of infinite nuclear matter with a small proton/neutron fraction is a crucial input to determine the properties of neutron stars and compare model predictions to astronomical observations. The so-called `symmetry energy' is the part of the EOS accounting for the difference in the number of neutrons and protons. Numerous experiments have been devised to assess the symmetry energy and constrain its functional dependence with the nucleon density. Further constraints follow from a stellar modeling using the EOS to reproduce astronomical observations such as neutron star masses and radii. The recent detection of gravitational waves emitted from neutron star mergers and the nucleosynthesis ensuing from these events caused a surge of interest for such studies. Several types of nuclear reactions have been proposed to study the symmetry energy part of the EOS. Some of them consist in determining the neutron skin in nuclei and exploit its correlation with the slope parameter of the symmetry energy. In this article we explore a particular set of reactions using high energy ($E_{lab} \sim 1$ GeV/nucleon) neutron-rich projectiles. We explore measurements of all reaction fragments (a) in the same isotopic chain, i.e., only by removal of neutrons, (b) in all charge-changing channels, and (c) total interaction cross sections. Using Hartree-Fock-Bogoliubov (HBF) predictions for neutron and proton densities with Skyrme interactions, we explore the sensitivity of these cross sections with the neutron skin in nuclei.

We examine the correlations of neutron star radii with the nuclear matter incompressibility, symmetry energy, and their slopes, which are the key parameters of the equation of state (EoS) of asymmetric nuclear matter. The neutron star radii and the EoS parameters are evaluated using a representative set of 24 Skyrme-type effective forces and 18 relativistic mean field models, and two microscopic calculations, all describing 2$M_\odot$ neutron stars. Unified EoSs for the inner-crust-core region have been built for all the phenomenological models, both relativistic and non-relativistic. Our investigation shows the existence of a strong correlation of the neutron star radii with the linear combination of the slopes of the nuclear matter incompressibility and the symmetry energy coefficients at the saturation density. Such correlations are found to be almost independent of the neutron star mass in the range $0.6\text{-}1.8M_{\odot}$. This correlation can be linked to the empirical relation existing between the star radius and the pressure at a nucleonic density between one and two times saturation density, and the dependence of the pressure on the nuclear matter incompressibility, its slope and the symmetry energy slope. The slopes of the nuclear matter incompressibility and the symmetry energy coefficients as estimated from the finite nuclei data yield the radius of a $1.4M_{\odot}$ neutron star in the range $11.09\text{-}12.86$ km.

Context. Neutron stars are currently studied with an rising number of electromagnetic and gravitational-wave observations, which will ultimately allow us to constrain the dense matter equation of state and understand the physical processes at work within these compact objects. Neutron star global parameters, such as the mass and radius, can be used to obtain the equation of state by directly inverting the Tolman-Oppenheimer-Volkoff equations. Here, we investigate an alternative approach to this procedure. Aims. The aim of this work is to study the application of the artificial neural networks guided by the autoencoder architecture as a method for precisely reconstructing the neutron star equation of state, using their observable parameters: masses, radii, and tidal deformabilities. In addition, we study how well the neutron star radius can be reconstructed using only the gravitational-wave observations of tidal deformability, that is, using quantities that are not related in any straightforward way. Methods. The application of an artificial neural network in the equation-of-state reconstruction exploits the non-linear potential of this machine learning model. Since each neuron in the network is basically a non-linear function, it is possible to create a complex mapping between the input sets of observations and the output equation-of-state table. Within the supervised training paradigm, we construct a few hidden-layer deep neural networks on a generated data set, consisting of a realistic equation of state for the neutron star crust connected with a piecewise relativistic polytropes dense core, with its parameters representative of state-of-the art realistic equations of state. Results. We demonstrate the performance of our machine-learning implementation with respect to the simulated cases with a varying number of observations and measurement uncertainties. Furthermore, we study the impact of the neutron star mass distributions on the results. Finally, we test the reconstruction of the equation of state trained on parametric polytropic training set using the simulated mass–radius and mass–tidal-deformability sequences based on realistic equations of state. Neural networks trained with a limited data set are capable of generalising the mapping between global parameters and equation-of-state input tables for realistic models.