Estimating the variation of neutron star observables by symmetric dense nuclear matter properties

Recent multi-channel astrophysics observations and the soon-to-be published new measured electromagnetic and gravitation data provide information on the inner structure of the compact stars. These macroscopic observations can significantly increase our knowledge on the neutron star enteriors, providing constraints on the microscopic physical properties. On the other hand, due to the masquarade problem, there are still uncertainties on the various nuclear-matter models and their parameters as well. Calculating the properties of the dense nuclear matter, effective field theories are the most widely-used tools. However the values of the microscopical parameters need to be set consistently to the nuclear and astrophysical measurements. In this work we investigate how uncertainties are induced by the variation of the microscopical parameters. We use a symmetric nuclear matter in an extended $\sigma$-$\omega$ model. We calculate the dense matter equation of state and give the mass-radius diagram. We present that the Landau mass and compressibility modulus of the nuclear matter have definite linear relation to the maximum mass of a Schwarzschild neutron star.


Introduction
The investigation of the structure of compact astrophysical objects like neutron stars, magnetars, quark-or hybrid stars etc. is an active novel research area as a child of astrophysics, gravitational theory and experiment and nuclear physics. So far the extreme dense state of the matter can not be produced in today's Earth-based particle accelerators, thus only celestial objects can be used for tests. Electromagnetic measurements, such as X-ray-and gamma satellites, aim to measure properties of these objects more and more accurately [1][2][3][4]. In parallel, radio array data [5] and the newly discovered gravity waves provide a new way to probe their inner structure [6][7][8]. These observations are particularly important inputs for the theoretical studies of dense nuclear matter [9,10].
From the theoretical point of view, first principle calculations based on lattice field theory are still challenging at high chemical potentials present in compact stars [11][12][13]. Thus, effective theories play an important role in studying the properties of cold dense nuclear matter [14,15]. Recent studies show the importance of the correct handling of the bosonic sector in effective theories of nuclear matter [16,17], moreover, applying the functional renormalization group (FRG) method on the simplest non-trivial nuclear matter, the effect of the microscopical parameters on neutron star observables were shown in Refs. [18,19].
In this paper we study the connection between the parametrizations of effective nuclear models and measurable properties of compact stars in three differently extended versions of the σ-ω model. All these include symmetric nuclear matter with various interaction terms in the bosonic sector. After calculating the equation of state (EoS) corresponding to different parametrizations of these models, the mass-radius (M-R) diagrams are calculated by solving the Tollmann -Oppenheimer -Volkoff (TOV) equations. We show how sensitive the mass-radius relation is to differences in the bosonic sector.
The dependence of particular properties of compact stars (maximum mass and radius) is presented, influenced by different saturation parameters of the symmetric nuclear matter.

The extended σ-ω in the mean field approximation
Here we apply the most common mean field model of the dense nuclear matter, formulating the extended σ-ω model [20,21] with the Lagrange-function taken from Refs. [22,23], where Ψ is the fermionic nucleon field, N f = 2 is the number of nucleons, and m N , m σ , and m ω are the nucleon, sigma, and omega masses, respectively, for the usual scalar and vector fields. We introduced the ω µν = ∂ µ ω ν − ∂ ν ω µ , and the Yukawa coupling corresponding to the σ-nucleon and ω-nucleon interactions is given by g σ and g ω . We denote the general bosonic interaction terms with U i (σ) which can have thee different forms as the considered modified model cases for certain i, In the mean field (MF) approximation the kinetic terms are zero for the mesons and only the fermionic path integral has to be calculated at finite chemical potential and temperature. We consider here the symmetric nuclear matter to be in equilibrium which includes the baryon number conservation. Taking into account this, the standard procedure were applied by minimizing the free energy of the infinite symmetric nuclear matter at the zero temperature limit, where for the proton (n p ) and neutron (n n ) number densities are equal, such as the proper chemical potentials, µ p and µ n respectively.
After applying this for all the three cases in eq. (2) substituting into eq. (1), numerical solution can be obtained after parameter fitting.

Parameter fitting in the extended σ-ω model
As the general procedure, all the models considered cases in eq. (2) need to fit to the nucleon saturation data found in e.g. Refs. [22,24]. In parallel to the effective mass, we introduced the definition of the Landau mass Where k = k F the Fermi-surface and E k is the dispersion relation of the nucleons. The Landau mass (m L ) and the effective mass (m * ) are not independent in relativistic mean field theories, This is the reason why the Landau mass and the effective mass of the nucleons can not be fitted simultaneously in the models we consider [24]. In this paper we deal with this problem in the following way. We fit all of the models two times: using the effective mass value from Table 1 and one calculated from eq. (5) to reproduce the Landau mass value from Table 1. If the models with U 3 and U 4 type interaction terms are used, than there is not enough free parameters to fit the data in Table 1. In these cases the nucleon effective mass, saturation density, and binding energy are fitted and the compression modulus is a prediction, given by which has a simple connection to the thermodynamical compressibility at the saturation density n 0 .
In the case of U 34 , all four parameters can be fitted simultaneously, and there is another way to incorporate data regarding both Landau and effective mass. For this model we consider a third fit, where the value of the effective mass is chosen in a way that minimizes the error coming from not fitting the two types of masses correctly. Technically this value of the effective mass minimized the χ 2 of the fit, with value Since the incompressibility is different for the model cases with different interactions terms, we compared them in Table 2. For model cases with U 3 and U 4 there are two fits, for Landau and effective mass which produce different incompressibility values, because they do not have enough free parameters to fit the correct value. However, for U 34 there are enough parameters to fit the incompressibility, so it has the same value for all three fits: for the Landau mass, for the effective mass and for the optimal mass. As Table 2 presents incompresibility values for U 3 with Landau fit is quite close and U 34 result provide the best fit with the saturation nuclear matter parameters in Table 1. These models differ in their predictions for higher densities of nuclear matter which complicates the description of the compact star interior.

Properties of nuclear matter in the extended σ-ω model
The nuclear properties of the different model cases were compared with all the possible parameter fits at the equation of state (EOS) level. We used all three types of interaction terms: U 3 , U 4 , and U 34 , and each was considered with two parametrizations corresponding to Landau and effective mass fits. All results were cross-checked with the original σ-ω model parameters. In case of the model characterized by U 34 interactions, we used a fit which reproduces (7).
The energy density, pressure and density were calculated in all of these models. The equation of states corresponding to these model and fit cases are shown on Similar grouping can be seen on Figure 2 presenting the binding energy B as the function of the nuclear density n. However, near the minimal B, the value of incompressibility governs the curves, because it determines the curvature of the curves around minimum. Effective mass fits with open symbols has steeper rise with the increasing density, n, while Landau-mass fit curves with full symbols are wider. We note that below the saturation density curves are getting more independent on the model choice and parameter fits.

Compact star observables from the effective models
The Tolman -Oppenheimer -Volkoff equations provide the general relativistic description of the compact stars assuming spherically symmetric and time independt space-time structure [22], Here P(r) and ε(r) are the pressure and energy density as functions of the radius of the star, G is the gravitational constant while m(r) is the mass of star which is included in a shell with radius r. To integrate the equations, one needs a connection between P(r) and ε(r) at given r which is provided by the nuclear matter equation of state in the form of the relation P(r) = P(ε(r)). To start the integration one has to choose a central energy density value ε c for the star as an initial condition. After solving the equations (8) using the EoS from the model cases with various fits, the mass (M) and radius (R) of a compact star with a given energy density can be determined. Results corresponding to different energy densities in a given model are summarized on a mass-radius M -R diagram on Figure 3.  The model variants inherited the behaviour as on Figs. 1 and 2: the curves are grouped based on whether or not they parametrized to reproduce the effective mass or Landau mass. Models with smaller effective mass (open symbols) systematically produce higher maximum mass stars compared to their parametrization with larger effective mass (Landau mass, with full symbols). Moreover, all models fitted for the effective mass value in Table 1 produce higher maximum star mass than the ones fitted for the Landau mass. Since the Landau mass and effective mass are not independent as in eq. (4), the above statement is equivalent to saying that higher effective mass produces smaller maximum star mass. This picture is supported by the curve corresponding to the model case parametrized by the optimal mass in eq. (7) which is the best fit of the model. The maximum star mass in this case is between the values produced by parametrizations described by smaller and larger effective nucleon mass. It is interesting to note that parametrizations corresponding to the effective mass and to the optimal mass value may be ruled out by observations, as the produce much larger maximum star mass and radius than the most recent measurements and theoretical predictions suggest [28].
The value of incompressibility in these models seems to influence very little the maximum star mass, however, at the same time it has an effect on the compactness by influencing the radius of the star. Thus can be seen on Fig. 3 by looking at the curves belonging to the fits with the same effective mass (Landau mass). In these model cases all of the parameters are the same which listed in Table 1 apart form the incompressibility.

Connection between maximum star mass and nuclear parameters
Our results present strong connection between the nuclear effective mass, m L and maximum star mass, M max M . To further study this phenomena, the model variant with interaction term U 34 is used and presented below. The couplings in the model are calculated to reproduce all values listed in Table 1 except for the nuclear effective mass and Landau mass which are kept as external parameters. To study the maximum star mass dependence on effective nucleon mass, many different fits of this model are considered all with different values of effective nucleon mass. The M -R diagrams corresponding to these parametrizations are calculated, and the maximum star mass is determined in each case. This procedure makes it possible to determine the dependence of maximum star mass on nucleon effective mass in this model, if everything else is kept constant. The results are summarized in Fig. 4. The connection between maximum star mass and nucleon effective mass is well described by a linear connection, which gives the best fit of the numerical data, where M max M is given in units of M Sun and m L is given in MeV.
Although the above equation is derived considering only the model with interaction term U 34 and taking into account parametrizations which differ only in the value of nucleon effective mass, it generalizes well and it approximates the maximum star mass corresponding to the other model variants we consider with very high accuracy. This seems to indicate that the linear connection is a good approximation regardless of the bosonic interaction term used, it holds across different model variants.
This observations is further supported by the connection between the radius of the maximum mass star (R max M ) and nuclear Landau mass. This connection is also present strong linear dependence and holds across all model variants we considered in this paper. The results are show on Fig. 4, where the plotted linear relation is described by applying R is given in km units and m L is given in MeV.
Since besides the nucleon Landau (effective) mass the only parameter which is different in our model variants is the incompressibility, it is worth to study the relation between incompressibility and maximum star mass while keeping every other parameter constant. For this end in the model with U 34 interaction term, the Landau mass fixed for the value which is listed in Table 1     where M max M is measured in units of M Sun , K in MeV units, and R max M is in km. The slope of both linear functions is tiny and even double the K produces only a few percent of variation in the maximum mass star mass and radius. These results provide a heuristic understanding of the previous results on the sense matter properties: the maximum mass star radius and mass depends strongly on nucleon Landau mass, since it's dependence on incompressibility is negligible compared to how strongly it depends on the value of effective mass and all other parameters are kept constant.

Summary
We investigate the macroscopical observables of neutron stars with a symmetric dense nuclear and matter in the exterior. We compared three cases with different general bosonic interactions, we used variations of the nuclear parameter fit method. We investigated how these bosonic interactions matter in the nuclear equation of state, and especially how the effective mass and Landau mass parameters play role in the nuclear potentials.
With variation of these microscopical parameters we explored the mass-radius groups in case of a Schwarzschild neutron star model. We presented that the Landau mass of the nuclear matter have definite linear relation to the maximum mass and maximum radius of the neutron star, considered within the best U 34 model case. We also obtained that the effect of the incompressibility is negligible on these macroscopical parameters, with M max M ≈ 1.78 and R max M = 7.87 within the regime K =200-400, meaning 5% variation for both cases.
These results support the idea, that evolution equation of the strongly interacting matter saves the magnitude of the uncertainty originated from the microscopical physical parameters, which appears as deterministic variation of the macroscopical observables as suggested in Refs. [18,19].