The Equation of State of Nuclear Matter : from Finite Nuclei to Neutron Stars

{\it Background.} We investigate possible correlations between neutron star observables and properties of atomic nuclei. Particularly, we explore how the tidal deformability of a 1.4 solar mass neutron star, $M_{1.4}$, and the neutron skin thickness of ${^{48}}$Ca and ${^{208}}$Pb are related to the stellar radius and the stiffness of the symmetry energy. {\it Methods.} We examine a large set of nuclear equations of state based on phenomenological models (Skyrme, NLWM, DDM) and {\it ab-initio} theoretical methods (BBG, Dirac-Brueckner, Variational, Quantum Monte Carlo). {\it Results.} We find strong correlations between tidal deformability and NS radius, whereas a weaker correlation does exist with the stiffness of the symmetry energy. Regarding the neutron skin thickness, weak correlations appear both with the stiffness of the symmetry energy, and the radius of a $M_{1.4}$. {\it Conclusion.} The tidal deformability of a $M_{1.4}$ and the neutron-skin thickness of atomic nuclei show some degree of correlation with nuclear and astrophysical observables, which however depends on the ensemble of adopted EoS.


Introduction
The equation of state (EoS) of isospin asymmetric nuclear matter plays a major role in many different realms of modern physics, being the fundamental ingredient for the description of heavy-ion collision dynamics, nuclear structure, static and dynamical properties of neutron stars (NS), core-collapse supernova and binary compact-star mergers [1,2]. In principle, it can be expected that in heavy ion collisions at large enough energy nuclear matter is compressed at density a few times larger than the nuclear saturation density, and that, at the same time, the two collision partners produce flows of matter, which should be connected with the nuclear EoS. In the physics of compact objects, the central density likely reached in the inner core of a NS may reach values up to one order of magnitude larger than the saturation density, and this poses several theoretical problems because a complete theory of nuclear interactions at arbitrarily large values of density, temperature and asymmetry, should in principle be derived from the quantum chromodynamics (QCD), and this is a very difficult task which presently cannot be realised. Therefore, theoretical models and methods of the nuclear many-body theory are required to build the EoS, which has to be applied and tested in terrestrial laboratories for the description of ordinary nuclear structure, and in astrophysical observations for the study of compact objects.
Among possible observables regarding NS, the mass and radius are the most promising since they encode unique information on the EoS at supranuclear densities. Currently the masses of several NSs are known with good precision [3][4][5][6][7], but the information on their radii is less accurate [8,9]. The recent observations of NICER [10] have reached a larger accuracy for the radius, but future planned missions like eXTP [11] will allow us to statistically infer NS mass and radius to within a few percent.
A big step forward is represented by the recent detection by the Advanced LIGO and Virgo collaborations of gravitational waves emitted during the GW170817 NS merger event [12][13][14]. This has provided important new insights on the mass and radii of these objects by means of the measurement of the tidal deformability [15,16], and allowed to deduce upper and lower limits on it [13,17].
In this paper we analyze the constraints on the nuclear EoS obtained from the analysis of the NS merger event GW170817, and try to select the most compatible EoS chosen among those derived from both phenomenological and ab-initio theoretical models. We also examine possible correlations among properties of nuclear matter close to saturation with the observational quantities deduced from GW170817 and nuclear physics experiments. In particular, we concentrate on the tidal deformability of NS, and the neutron skin thickness in finite nuclei, thus connecting astrophysical observables with laboratory nuclear physics.
The paper is organized as follows. In Sect. 2 we give a schematic overview of NS phenomenology, whereas in Sect. 3 we explain the role of the equation of state in determining the main properties of NS, and illustrate the ones we adopt in the present study. The experimental constraints on the nuclear EoS are presented in Sect. 3.1 whereas the astrophysical ones are discussed in Sect. 3.2. A brief overview of different EoS of β-stable matter is given in Sect. 4, along with numerical results. In Sect. 5 we briefly discuss the NS tidal deformability, and its connection to the neutron skin thickness in Sect. 6. Conclusions are drawn in Sect. 7.

Neutron stars in a nutshell
Neutron stars are a type of stellar compact remnant that can result from the gravitational collapse of an ordinary star with a mass in the range 8 − 25M ⊙ (with M ⊙ ≈ 2 × 10 33 g the mass of the Sun) during a Type II, Ib or Ic supernova event. A supernova explosion will occur when the star has exhausted its possibilities for energy production by nuclear fusion. Then, the pressure gradient provided by the radiation is not sufficient to balance the gravitational attraction, becoming the star unstable and, eventually, collapsing. The inner regions of the star collapse first and the gravitational energy is released and transferred to the outer layers of the star blowing them away.
NS are supported against gravitational collapse mainly by the neutron degeneracy pressure and may have masses in the range M ∼ 1 − 2M ⊙ (with M ⊙ ≈ 2 × 10 33 g the mass of the Sun) and radii of about 10 − 12 km. A schematic cross section the predicted "onion"-like structure of the NS interior is shown in Fig. 1. At the surface, densities are typically ρ < 10 6 g/cm 3 . The outer crust, with densities ranging from 10 6 g/cm 3 to 10 11 g/cm 3 is a solid region where heavy nuclei, mainly around the iron mass number, is a Coulomb lattice coexist in β-equilibrium (i.e., in equilibrium with respect to weak interaction processes) with an electron gas. Moving towards the center the density increases and the electron chemical potentials increases and the electron capture processes on nuclei opens and the nuclei become more and more neutron-rich. At densities of ∼ 4.3 × 10 11 g/cm 3 the only available levels for the neutrons are in the continuum and they start to "drip out" of the nuclei. We have then reached the inner crust region, where matter consist of a Coulomb lattice of very neutron-rich nuclei together with a superfluid neutron gas and an electron gas. In addition, due to the competition between the nuclear and Coulomb forces, nuclei in this region lose their spherical shapes and presenyt more exotic topologies (droplets, rods, cross-rods, salabs, tubves, bubbles) giving rise to what has been called "nuclear pasta" phase [18]. At densities of ∼ 10 14 g/cm 3 nuclei start to dissolve and one enters the outer core. In this region matter is mainly composed of superfluid neutrons with a smaller concentration of superconducting protons and normal electrons and muons. In the deepest region of the star, the inner  11 10 g cm 6 14 Outer core: Hyperons ?
Inner core: core, the density can reach values of ∼ 10 15 g/cm 3 . The composition of this region, however, is not known, and it is still subject of speculation. Suggestions range from a hyperonic matter, meson condensates, or deconfined quark matter. The observation of NS requires different types of ground-based and on-board telescopes covering all bands of the electromagnetic spectrum. Radio observations are carried out with ground-based antennas located in different places of the Earth. Three examples of these radio telescopes are the Arecibo radio telescope in Puerto Rico, the Green Bank Observatory in West Virginia, and the Nançay decimetric radio telescope in France. Observations in the near infrared and the optical bands can be performed with the use of large ground-based telescopes such as the Very Large Telescope (VLT) in the Atacama desert in Chile. The Hubble-Space Telescope (HST) can be used used to cover the optical and ultraviolet regions. Observations in the extreme ultraviolet, X-ray and γ-ray require the use of space observatories such as the Chandra X-ray Observatory (CXO), the X-ray Multi Mirror (XMM-Newton) and the Rossi X-ray Timing Explorer (RXTE) in the case of X-ray observations; and the High Energy Transient Explorer (HETE-2), the International Gamma-Ray Astrophysics Lboratory (INTEGRAL) and the Fermi Gamma-ray Space Telescope (FGST), in the case of γ-ray ones.
Information on the properties of NS additional to that obtained from the observation of their electromagnetic radiation can be provided from the detection of the neutrinos emitted during the supernova explosion that signals the birth of the star. Examples of neutrino observatories are: the under-ice IceCube observatory placed in the South Pole; the under-water projects ANTARES (Astronomy with a Neutrino Telescope and Abyss environmental REsearch) and the future KM3NET (Cubic Kilometre Neutrino Telescope) in the Mediterranean sea; and the underground observatories SNO (Sudbury Neutrino Observatory) located 2100 meters underground in the Vale's Creighton Mine in Canada, and the Kamioka observatory placed at the Mozumi Mine near the city of Hida in Japan.
The detection of gravitational waves, originated during the coalescence of two NS as in the GW170817 event recently detected by the Advanced LIGO and Advanced VIRGO collaborations [12][13][14] or from the oscillation modes of NS, represents nowadays a new way of observing these objects and constitutes a very valuable new source of information. In particular, observations of NS mergers can potentially provide stringent constraints on the nuclear EoS by comparing model predictions with the precise shape of the detected gravitational wave signal. The interested reader is referred to Ref. [19] for a recent review on this hot and exciting topic.

The nuclear equation of state
The theoretical description of nuclear matter under extreme density conditions is a very challenging task. Theoretical predictions in this regime are diverse, ranging from purely nucleonic matter with high neutron-proton asymmetry, to baryonic strange matter or a quark deconfined phase of matter. In this work we adopt a conventional description by assuming that the most relevant degrees of freedom are nucleons. Theoretical approaches to determine the nuclear EoS can be classified in two categories: phenomenological and microscopic (ab-initio).
Phenomenological approaches, either non-relativistic or relativistic, are based on effective interactions that are frequently built to reproduce the properties of nuclei [20]. Skyrme interactions [21,22] and relativistic mean-field (RMF) models [23,24] are among the most used ones. Many of such interactions are built to describe finite nuclei in their ground state, i.e. close to the isospin symmetric case and, therefore, predictions at high isospin asymmetries should be taken with care. For instance, most Skyrme forces are, by construction, well behaved close to nuclear saturation density ρ 0 ≈ 0.15 − 0.16 fm −3 and moderate values of the isospin asymmetry, but predict very different EoS for pure neutron matter, and therefore give different predictions for NS observables. In this work we use the 27 Skyrme forces that passed the restrictive tests imposed by Stone et al. in Ref. [20] over almost 90 existing Skyrme parametrizations. These forces are: GS and Rs [25], SGI [26], SLy0-SLy10 [27] and SLy230a [28,29] of the Lyon group, the old SV [30], SkI1-Sk5 [31] and SkI6 [32] of the SkI family, SkMP [33], SkO and SkO' [34], and SkT4 and SkT5 [35].
Similarly, relativistic mean-field models are based on effective Lagrangian densities where the interaction between baryons is described in terms of meson exchanges. The couplings of nucleons with mesons are usually fixed by fitting masses and radii of nuclei and the bulk properties of nuclear matter, whereas those of other baryons, like hyperons, are fixed by symmetry relations and hypernuclear observables. In this work we consider two types of RMF models: models with constant meson-baryon couplings described by the Lagrangian density of the nonlinear Walecka model (NLWM), and models with density-dependent couplings [hereafter referred to as density-dependent models (DDM)]. In particular, within the first type, we consider the models GM1 and GM3 [36], TM1 [37], NL3 and NL3-II [38], and NL-SH [39]. For the DDM, we consider the models DDME1 and DDME2 [40], TW99 [41], and the models PK1, PK1R and PKDD of the Pekin group [42].
Microscopic approaches, on other hand, are based on realistic two-and three-body forces that describe nucleon scattering data in free space and the properties of the deuteron. These interactions are based on meson-exchange theory [43,44] or, very recently, on chiral perturbation theory [45][46][47][48]. Then one has to solve the complicated many-body problem [49] in order to obtain the nuclear EoS. The main difficulty is the treatment of the short-range repulsive core of the nucleon-nucleon interaction. Different many-body approaches have been devised for the construction of the nuclear matter EoS, e.g., the Brueckner-Hartree-Fock (BHF) [50] and the Dirac-Brueckner-Hartree-Fock (DBHF) [51][52][53] theories, the variational method [54], the correlated basis function formalism [55], the self-consistent Green's function technique [56,57], the V low k approach [58] or Quantum Monte Carlo techniques [59,60].
As far as the microscopic approaches are concerned, in this paper we adopt several BHF EoS based on different nucleon-nucleon potentials, namely the Bonn B (BOB) [43,61], the Nijmegen 93 (N93) [62,63], and the Argonne V 18 (V18) [64]. In all those cases, the two-body forces are supplemented by nucleonic three-body forces (TBF), which are needed in order to reproduce correctly the saturation properties of nuclear matter. Currently a complete ab-initio theory of TBF is not available yet, and therefore we adopt either phenomenological or microscopic models [65][66][67][68]. The microscopic TBF employed in this paper are described in detail in Refs. [68,69], whereas a phenomenological approach based on the Urbana model [66,70,71], is also adopted. In this case the corresponding EoS is labelled UIX in Table 1. Within the same theoretical framework, we also studied an EoS based on a potential model which includes explicitly the quark-gluon degrees of freedom, named FSS2 [72,73]. This reproduces correctly the saturation point of symmetric matter and the binding energy of few-nucleon systems without the need of introducing TBF. In the following we use two different EoS versions labelled respectively as FSS2CC and FSS2GC. Moreover, we compare these BHF EoSs with the often-used results of the Dirac-BHF method (DBHF) [53], which employs the Bonn A potential, the APR EoS [54] based on the variational method and the V 18 potential, and a parametrization of a recent Auxiliary Field diffusion Monte Carlo (AFDMC) calculation of Gandolfi et al. given in Ref. [74].

Laboratory constraints on the nuclear EoS
Around saturation density ρ 0 and isospin asymmetry δ ≡ (N − Z)/(N + Z) = 0 [being N(Z) the number of neutrons (protons)], the nuclear EoS can be characterized by a set of few isoscalar (E 0 , K 0 ) and isovector (S 0 , L, K sym ) parameters. These parameters can be constrained by nuclear experiments and are related to the coefficients of a Taylor expansion of the energy per particle of asymmetric nuclear matter as a function of density and isospin asymmetry where x ≡ (ρ − ρ 0 )/3ρ 0 , E 0 is the energy per particle of symmetric nuclear matter at ρ 0 , K 0 the incompressibility and S 0 ≡ E sym (ρ 0 ) is the symmetry energy coefficient at saturation. The parameters L and K sym characterize the density dependence of the symmetry energy around saturation. These parameters are defined as The incompressibility K 0 gives the curvature of E(ρ) at ρ = ρ 0 , whereas S 0 determines the increase of the energy per nucleon due to a small asymmetry δ.
Properties of the various considered EoSs are listed in Table 1, namely, the value of the saturation density ρ 0 , the binding energy per particle E 0 , the incompressibility K 0 , the symmetry energy S 0 , and its derivative L at ρ 0 . Measurements of nuclear masses [78] and density distributions [79] yield E 0 = −16 ± 1 MeV and ρ 0 = 0.14 − 0.17 fm −3 , respectively. The value of K 0 can be extracted from the analysis of isoscalar giant monopole resonances in heavy nuclei, and results of Ref. [80] suggest K 0 = 240 ± 10 MeV, whereas in Ref.
[81] a value of K = 248 ± 8 MeV is reported. Even heavy ion collision experiments point to a "soft" EoS, i.e., a low value of K 0 , though the constraints inferred from heavy ion collisions are model dependent because the analysis of the measured data requires the use of transport models [82]. Experimental information on the symmetry energy at saturation S 0 and its derivative L can be obtained from several sources such as the analysis of giant [83] and pygmy [84,85] dipole resonances, isospin diffusion measurements [86], isobaric analog states [87], measurements of the neutron skin thickness in heavy nuclei [88][89][90][91][92] and meson production in heavy ion collisions [93]. However, whereas S 0 is more or less well established (≈ 30 MeV), the values of L (30 MeV < L < 87 MeV), and especially those of K sym (−400 MeV < K sym < 100 MeV) are still quite uncertain and poorly constrained [94,95]. The reason why the isospin dependent part of the nuclear EoS is so uncertain is still an open question, very likely related to our limited knowledge of the nuclear forces and, in particular, to its spin and isospin dependence.
From Table 1, we notice that all the adopted EoSs in this work agree fairly well with the empirical values. Marginal cases are the slightly too low E 0 and K 0 for V18, too large S 0 for N93, and too low K 0 for UIX and FSS2GC. We notice that the L parameter does not exclude any of the microscopic EoSs, whereas several phenomenological models predict too large L values.

Astrophysical constraints on the nuclear EoS
The main astrophysical constraints on the nuclear EoS are those arising from the observation of NS. An enormous amount of data on different NS observables have been collected after fifty years of NS observations. These observables include: masses, radii, rotational periods, surface temperatures, gravitational redshifts, quasi-periodic oscillations, magnetic fields, glitches, timing noise and, very recently, gravitational waves. In the next lines we briefly review how masses and radii are measured. Observational constraints derived from the recent observation of the gravitational wave signal from the merger of two NS detected by the Advanced LIGO and Advanced VIRGO collaborations [12][13][14] will be discussed in detail in Sect. 5.
NS masses can be directly measured from observations of binary systems. There are five orbital parameters, also known as Keplerian parameters, which can be precisely measured. They are the projection of the pulsar's semi-major axis (a 1 ) on the line of sight (x ≡ a 1 sin i/c, where i is inclination of the orbit), the eccentricity of the orbit (e), the orbital period (P b ), and the time (T 0 ) and longitude (ω 0 ) of the periastron. With the use of Kepler's Third Law, these parameters can be related to the masses of the NS (M p ) and its companion (M c ) though the so-called mass function where v 1 = 2πa 1 sin i/P b is the projection of the orbital velocity of the NS along the line of sight. The individual masses of the two components of the system cannot be obtained if only the mass function is determined. Additional information is required. Fortunately, deviations from the Keplerian orbit due to general relativity effects can be detected. The relativistic corrections to the orbit are parametrized in terms of one or more parameters called post-Keplerian. The most significant ones are: the combined effect of variations in the transverse Doppler shift and gravitational redshift around an elliptical orbit (γ), the range (r) and shape (s) parameters that characterize the Shapiro time delay of the pulsar signal as it propagates through the gravitational field of its companion, the advance of the periastron of the orbit (ω) and the orbital decay due to the emission of quadrupole gravitational radiation (Ṗ b ). These post-Keplerian parameters can be written in terms of measured quantities and the masses of the star and its companion (see e.g., Ref. [100] for specific expressions). The measurement of any two of these post-Keplerian parameters together with the mass function f is sufficient to determine uniquely the masses of the two components of the system. As the reader can imagine NS radii are very difficult to measure, the reason being that NS are very small objects and are very far away from us (e.g., the closest NS is the object RX J1856.5-3754 in the constellation Corona Australis which is about 400 light-years from the Earth). That is why there not exist direct measurements of NS radii yet. Nevertheless, it is possible to determine them by using the thermal emission of low-mass X-ray binaries (systems where one of the components is a NS and the companion a less massive object (M c < M ⊙ ) which can be a main sequence star, a red giant or a white dwarf). The observed X-ray flux (F) and estimated surface temperature (T) together with a determination of the distance (D) of the star, can be used to obtain the radius of the NS through the implicit relation Here σ is the Stefan-Boltzmann constant and M the mass of the NS. The major uncertainties in the measurement of the radius through Eq. (10) come from the determination of the temperature, which requires the assumption of an atmospheric model, and the estimation of the distance of the star. However, the analysis of present observations from quiescent low-mass X-ray binaries is still controversial (see e.g., Refs. [101,102]). We notice that the simultaneous measurement of both mass and radius of the same NS would provide the most definite observational constraint on the nuclear EoS. Very recently the NICER (Neutron Star Interior Composition Explorer) mission has reported a Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030+0451 [10]. The values inferred from the analysis of the collected data are 1.34 +0. 15 −0.16 M ⊙ and 12.71 +1.14 −1.19 km, respectively.

EoS for β-stable matter
In order to study the structure of the NS core, we have to calculate the composition and the EoS of cold, neutrino-free, catalyzed matter. As stated before, we consider a NS with a core of nucleonic matter without hyperons or other exotic particles. We require that it contains charge neutral matter consisting of neutrons, protons, and leptons (e − , µ − ) in β-equilibrium, and compute the EoS for charge neutral and β-stable matter in the following standard way [103]. The output of the many-body calculation is the energy density of lepton/baryon matter as a function of the different densities ρ i of the species i = n, p, e, µ , ǫ(ρ n , ρ p , ρ e , ρ µ ) = (ρ n m n + ρ p m p ) + (ρ n + ρ p )E(ρ n , ρ p ) + ǫ(ρ e ) + ǫ(ρ µ ) , where m i are the corresponding masses, and E(ρ n , ρ p ) is the enegy per particle of asymmetric nuclear matter. We have used ultrarelativistic and relativistic expressions for the energy densities of electrons ǫ(ρ e ) and muons ǫ(ρ µ ), respectively [103]. Since microscopic calculations are very time consuming in the case of these models we have used the parabolic approximation [104-108] of the energy per particle of asymmetric nuclear matter given in Eq. (2) with the symmetry energy calculated simply as the difference between the energy per particle of pure neutron matter E(ρ n = ρ, ρ p = 0) and symmetric nuclear matter Once the energy density (Eq. (11)) is known the various chemical potentials can be computed straightforwardly, and solving the equations for β-equilibrium, (b i and q i denoting baryon number and charge of species i) along with the charge neutrality, allows one to find the equilibrium composition ρ i at fixed baryon density ρ, and finally the EoS, Once the EoS of β-stable matter is known, one can determine the hydrostatical equilibrium configurations just solving the Tolman-Oppenheimer-Volkoff (TOV) [103] equations which describe the structure of a non-rotating spherically symmetric star in general relativity: where G is the gravitational constant, P the pressure, ǫ the energy density, m the mass enclosed within a sphere of radius r. The TOV equations have an easy interpretation. Consider a spherical shell of matter of radius r and thickness dr. The second equation gives the mass in this shell whereas the left hand side of the first one is the net force acting on the surface of the shell by the pressure difference dP(r). The first factor of the right hand side of this equation is the attractive Newtonian force of gravity acting on the shell by the mass interior to it. The remaining three factors result from the correction of general relativity. So the TOV equations express the balance at each r between the internal pressure as it supports the overlying material against the gravitational attraction of the mass interior to r. The integration of the TOV equations gives the mass and radius of the star for a given central density. It turns out that the mass of the NS has a maximum value as a function of radius (or central density), above which the star is unstable against collapse to a black hole. The value of the maximum mass depends on the nuclear EoS, so that the observation of a mass higher than the maximum mass allowed by a given EoS simply rules out that EoS. We now turn to the discussion of some results. We display in Fig.2 the β-stable matter EoS obtained for some of the models illustrated in Table 1, a limited sample of each class being plotted in one single panel. We see that the pressure is a monotonically increasing function of the energy density for all EoS. Each EoS is characterized by a given stiffness, which determines the maximum mass value of a NS: the stiffer the EoS the larger the maximum mass predicted.
The corresponding mass-radius relation is displayed in Fig.3. The observed trend is consistent with the EoS displayed in Fig.2. As expected, when the EoS stiffness increases the NS maximum mass increases as well. The considered EoS are compatible with the largest masses observed up to now, M max > 2.14 +0.10 −0.09 [7] for the object PSR J0740+6620 (cyan hatched area), and PSR J0348+0432 [6], M G = 2.01 ± 0.04 M ⊙ (red hatched area). We notice that recent analysis of the GW170817 event indicate also an upper limit of the maximum mass of about 2.2-2.3 M ⊙ [109][110][111][112], with which most of the models shown in the figure are compatible. The box shows the estimation of the mass ( Figure 3. Mass-radius relation predicted by the different EoS displayed in Fig.2. The observed masses of the millisecond pulsars PSR J0740+6620 [4] and PSR J0348+0432 [6] are also shown. The box shows the constraints inferred from the observations reported by the NICER mission [10]. See text for details.

The neutron star tidal deformability
Recently the tidal deformability λ, or equivalently the tidal Love number k 2 of a NS [113][114][115], has been recognized to provide valuable information and constraints on the related EoS, because it strongly depends on the compactness of the object, i.e. β ≡ M/R. More specifically, the Love number with Λ ≡ λ/M 5 , can be obtained by solving the TOV equations (17), along with the following first-order differential equation [116], with the EoS P(ε) as input, c 2 s = dP/dε the speed of sound, and boundary conditions given by being y R ≡ y(R), and the mass-radius relation M(R) provided by the condition P(R) = 0 for varying central pressure P c . For an asymmetric binary NS system, (M, R) 1 + (M, R) 2 , with mass asymmetry q = M 2 /M 1 , and known chirp mass M c , which characterizes the GW signal waveform, the average tidal deformability is defined bỹ Λ = 16 13 (1 + 12q)Λ 1 + (q + 12)Λ 2 (1 + q) 5 From the analysis of the GW170817 event [12][13][14], a value of M c = 1.186 +0.001 −0.001 M ⊙ was obtained, corresponding to M 1 = M 2 = 1.36 M ⊙ for a symmetric binary system, q = 0.73 − 1 andΛ < 730 from the phase-shift analysis of the observed signal. It turns out that, requiring both NSs to have the same EoS, leads to constraints 70 < Λ 1.4 < 580 and 10.5 < R 1.4 < 13.3 km [13] for a 1.4 solar mass NS.
However the high luminosity of the kilonova AT2017gfo following the NS merger event, imposes a lower limit on the average tidal deformability, Eq. (22),Λ > 400, which was deduced in order to justify the amount of ejected material heavier than 0.05 M ⊙ . This constraint could indicate that R 1. 4 12 km, which was used in Refs. [117][118][119][120] in order to constrain the EoS. This lower limit has to be taken with great care and, in fact, it has been recently revised toΛ 300 [121], but considered of limited significance in Ref. [122].   Table 1. The orange box indicate the experimental and observational constraints on L (see Table 1), Λ 1.4 and R 1.4 [13]. The violet line indicates a linear (quadratic) fit of the EoS data. The values of the corresponding correlation factors are also given. See text for details.
One of the main theoretical issues, following the detection of gravitational waves from NS mergers, regards the possibility of finding correlations between properties of nuclear matter and NS observables [123]. Along this same line, we further explore this issue, and using the set of microscopic EoS and the several Skyrme forces and relativistic models listed in Table 1, in the left panel of Fig. 4 we show the tidal deformability of a 1.4 solar mass NS as a function the symmetry energy parameter L at saturation density. The orange box shows the constraint on Λ 1.4 inferred from the observational data of the GW170817 event [13] together with the experimental limits of L reported in Table 1. We observe some degree of correlation between the tidal deformability and L , for which we can estimate the so-called correlation factor r, defined as being n the number of data pairs,L andΛ 1.4 the mean values of L and Λ 1.4 ; and σ L over the data set and σ Λ 1.4 their standard deviations. We get a value r = 0.817, which indicates a weak correlation. We note that several EoS lie outside the orange observational band. In particular, we notice that all DDM EoS (blue diamonds), except TW99, are not compatible with the data, as well as all the NLWM EoS (red squares). On the other hand, most of the Skyrme interactions lie within the orange band, with a few cases incompatible with observations because the predicted L values lie outside the experimental range, and some other are marginally compatible. As far as microscopic calculations are concerned, they are all in agreement with GW observations, except the DBHF EoS. In the right panel of Fig. 4 we report the tidal deformability as a function of the radius for a 1.4 solar mass NS, R 1.4 , for the same set of EoS. The observational constraints on Λ 1.4 and R 1.4 from GW170817 [13] are shown by the orange box. Contrary to the weak Λ 1.4 − L correlation found, we observe a strong quadratic correlation between Λ 1.4 and R 1.4 the correlation factor being in this case r = 0.986. This strong correlation was already noticed in Ref. [124] using a different set of EoS based again on Skyrme and relativistic mean field models. The behaviour of the microscopic and phenomenological EoS look very similar.

The neutron skin thickness
As stated in the previous Section, correlations between astrophysical observations and microscopic constraints from nuclear measurements, could help to better understand the properties of nuclear matter. For this purpose, the limits derived for the tidal deformability in GW170817 could be very valuable and exploited for studying the neutron skin thickness, defined as the difference between the neutron (R n ) and proton (R p ) root-mean-square radii: δR = r 2 n − r 2 p . It has been shown that this is strongly correlated to both L and to the radius of low-mass NS, since the size of a NS and the neutron skin thickness originate both from the pressure of neutron-rich matter, and are sensitive to the same EoS. As shown by Brown and Typel [88,89], and confirmed later by other authors [90,[125][126][127][128], the neutron skin thickness calculated in mean field models with either non-relativistic or relativistic effective interactions, is very sensitive to the density dependence of the nuclear symmetry energy, and, in particular, to the slope parameter L at normal nuclear saturation density. Using the Brueckner approach and the several Skyrme forces and relativistic models considered here, the authors of Ref. [129] made an estimation of the neutron skin thickness of 208 Pb and 132 Sn, adopting the suggestion of Steiner et al. in Ref. [125], where δR is calculated to lowest order in the diffuseness corrections as δR ∼ 3 5 t, being t the thickness of semi-infinite asymmetric nuclear matter In the above expression E s is the surface energy taken from the semi-empirical mass formula equal to 17.23 MeV, r 0 is obtained from the normalization condition (4πr 3 0 /3)(0.16) = 1, and δ c is the isospin asymmetry in the center of the nucleus taken as δ c = δ/2 according to Thomas-Fermi calculations. In this paper, we use the same prescription for the calculation of the neutron skin thickness of 208 Pb and 48 Ca, and we show the results in Fig. 5. The orange bands represent the predicted data for 48 Ca (left panels) for which the Calcium Radius Experiment (CREX) has not been run yet [131], whereas in the right panels experimental data obtained in the Lead Radius Experiment (PREX) [130] for 208 Pb, δR = 0.33 +0. 16 −0.18 fm, are plotted. In the upper panels, results are shown for the neutron skin thickness as a function of the derivative of the symmetry energy L. We notice that all the theoretical predictions from phenomenological models and some of the microscopic ones show some correlation between the neutron skin thickness and L, as indicated by the linear fits (violet curve) and by the value of the correlation coefficient, r = 0.803 for 48 Ca (r = 0.800 for 208 Pb ). Almost all the microscopic EoS turn out to be compatible with the PREX experimental data, whereas some phenomenological models, e.g. those of the NLWM class, give predictions out of the experimental range. The linear increase of δR with L is not surprising since the thickness of the neutron skin in heavy nuclei is determined by the pressure difference between neutrons and protons, which is proportional to the parameter L, that is, P(ρ 0 , δ) ≈ Lρ 0 δ 2 /3. On

Conclusions
In this work we have analyzed the existence of possible correlations between NS observables and properties of atomic nuclei. In particular, we have examined correlations of the tidal deformability Λ 1.4 of a 1.4M ⊙ NS and the neutron skin thickness δR of 48 Ca and 208 Pb with the stellar radius R 1.4 and the symmetry energy derivative L. To such end we have used a large set of different models for the nuclear equation of state, that include microscopic calculations based on the Brueckner-Hartree-Fock and Dirac-Brueckner-Hartree-Fock theories, the variational method and Quantum Monte Carlo techniques, and several phenomenological Skyrme and relativistic mean field models. We have found a strong quadratic correlation between Λ 1.4 and R 1.4 in agreement with the results of the recent work by Tsang et al. [124]. On the contrary, we have observed a weaker linear correlation between Λ 1.4 and L. Our results have confirmed the existence of a quite linear correlation between the neutron skin thickness of 48 Ca and 208 Pb with L, already pointed out by several authors using nonrelativistic and relativistic phenomenological models. A much weaker correlation has been found between δR and R 1.4 . The existence of these correlations, predicted by models based on approaches of different nature, suggest that their origin goes beyond the mean field character of the models employed.
To select the most compatible EoS among the ones predicted by the different models considered in this work, we have employed the experimental constraints on L and δR together with the observational ones on the mass, radius and tidal deformability imposed by the mass measurement of the millisecond pulsars PSR J1614-2230 [4] and PSR J0348+0432 [6], the GW170817 NS merger event [12][13][14] and the data of the NICER mission [10]. Our results have shown that only five microscopic models (BOB, V18, N93, UIX and DBHF) and four Skyrme forces (SGI, SkMP, SkO and SkO') are simultaneoulsy compatible with the present constraints on L (30 MeV < L < 87 MeV) and the PREX experimental data on the 208 Pb neutron skin thickness. All the NLWM and DDM models and the majority of the Skyrme forces are excluded by these two experimental constraints. We have also found that almost all the models considered are compatible with the largest masses observed up to now, M max > 2.14 +0. 10 −0.09 [7] for the object PSR J0740+6620, and PSR J0348+0432 [6], M G = 2.01 ± 0.04 M ⊙ , and with the upper limit of the maximum mass of about 2.2-2.3 M ⊙ [109][110][111][112] deduced from the analysis of the GW170817 event. Finally, we have seen that the estimation of the mass (1.34 +0. 15 −0.16 M ⊙ ) and equatorial radius (12.71 +1.14 −1.19 km) of the millisecond pulsar PSR J0030+0451 inferred from the Bayesian analysis of the data collected by the NICER mission [10] excludes most of the NLWM EoS considered in this work.
The major experimental, observational and theoretical advances on understanding the nuclear EoS done in the last decades has lead to constrain rather well its isoscalar part. Nevertheless, the isovector part of the nuclear EoS is less well constraint due mainly to our still limited knowledge of the nuclear force and, in particular, of its in-medium modifications and its spin and isospin dependence. Future laboratory experiments being planned in existing or next-generation radioactive ion beam facilities together with further NS observations, particularly a precise simultaneous measurement of the mass and radius of a single object, are fundamental to provide more stringent constraints on the nuclear EoS, and are very much awaited.