Nuclear Science and Techniques

Special Section on the Celebration for 55 Years’ Dedicated Research on Heavy-Ion Physics of Natowitz and his 80th Birthday

Constraints on the skewness coefficient of symmetric nuclear matter within the nonlinear relativistic mean field model

Bao-Jun Cai，

Lie-Wen Chen

Vol.28, No.12

Article number 185

Published in print 01 Dec 2017

Available online 05 Dec 2017

DOI：10.1007/s41365-017-0329-1

7902

Within the nonlinear relativistic mean field (NL-RMF) model, we show that both the pressure of symmetric nuclear matter at supra-saturation densities and the maximum mass of neutron stars are sensitive to the skewness coefficient, J₀, of symmetric nuclear matter. Using experimental constraints on the pressure of symmetric nuclear matter at supra-saturation densities from flow data in heavy ion collisions and the astrophysical observation of a large mass neutron star PSR J0348+0432, with the former favoring a smaller J₀ while the latter favors a larger J₀, we extract a constraint of -494 MeV≤J₀≤-10 MeV based on the NL-RMF model. This constraint is compared with the results obtained in other analyses.

Equation of state of nuclear matterHeavy-ion collisionsNeutron stars

1 Introduction

Determination of the equation of state (EOS) of asymmetric nuclear matter (ANM) is one of fundamental questions in contemporary nuclear physics and astrophysics. The exact knowledge on the EOS of ANM provides important information on the in-medium nuclear effective interactions which play a central role in understanding the structure and decay properties of finite nuclei as well as the related dynamical problems in nuclear reactions ([1-16]). The EOS of ANM also plays a decisive role in understanding a number of important issues in astrophysics including the structure and evolution of neutron stars as well as the mechanism of supernova explosion ([17-23]). Conventionally, the EOS of ANM is given by the binding energy per nucleon as functions of nucleon density, ρ, and isospin asymmetry, δ, i.e., E(ρ, δ), and some bulk characteristic parameters defined at the saturation density, ρ₀, of symmetric nuclear matter (SNM) are usually introduced to quantitatively characterize the EOS of ANM. For example, the energy, E₀(ρ₀), and incompressibility, K₀, of SNM, as well as the symmetry energy, E_sym(ρ₀), and its slope parameter, L, are the four famous lower-order bulk characteristic parameters of EOS of ANM. These bulk parameters defined at ρ₀ provide important information on both sub- and supra-saturation density behaviors of the EOS of ANM ([24, 25]).

Based on the empirical liquid-drop-like model analyses of high precision data about nuclear masses, the E₀(ρ₀) is well known to be about -16 MeV. The incompressibility has been determined to be K₀ = 240 ± 40 MeV from analyzing experimental data of nuclear giant monopole resonances (GMR) ([1, 26-30]). For E_sym(ρ₀) and L, the existing constraints extracted from terrestrial laboratory measurements and astrophysical observations are found to be essentially consistent with E_sym(ρ₀) = 32.5 ± 2.5 MeV and L = 55 ± 25 MeV (see, e.g., Refs. [31, 32]). While these lower-order bulk characteristic parameters have been relatively well determined or in significant progress, our knowledge on the higher-order bulk characteristic parameters remains very limited. Following E₀(ρ₀), K₀, E_sym(ρ₀), and L, the next bulk characteristic parameter should be the skewness coefficient, J₀, (also denoted as K’ or Q₀ in some literature) of SNM, which is related to the third-order density derivative of the binding energy per nucleon of SNM at ρ₀. The higher-order bulk characteristic parameter J₀ is expected to be important for the high density behaviors of nuclear matter EOS and thus may play an essential role in heavy ion collisions (HIC), the structure and evolution of neutron stars, supernova explosion, and gravitational wave radiation from merging of compact stars. To our best knowledge, there is very little experimental information on the J₀ parameter, and it is thus of great interest and critical importance to constrain the J₀ parameter, which is the main motivation of the present work.

Within the nonlinear relativistic mean field (RMF) model, we demonstrate in this work that the pressure of SNM at supra-saturation densities and the maximum mass of neutron stars provide good probes of the skewness coefficient, J₀. In particular, combining the experimental constraints on the pressure of SNM at supra-saturation densities from flow data in HIC and the recent astrophysical observation of a large mass neutron star, PSR J0348+0432, one can obtain a strong constraint on the J₀ parameter.

2 The skewness coefficient J₀ in nonlinear RMF model

2.1 Nuclear matter characteristic parameters

The EOS of isospin asymmetric nuclear matter, namely E(ρ,δ), can be expanded as a power series of even-order terms in δ as

E (ρ, δ) ≃ E_{0} (ρ) + E_{sym} (ρ) δ^{2} + O (δ^{4}),

(1)

where E₀(ρ)=E(ρ,δ =0) is the EOS of symmetric nuclear matter, and the symmetry energy is expressed as

E_{sym} (ρ) = {\frac{1}{2} \frac{\partial^{2} E (ρ, δ)}{\partial δ^{2}} |}_{δ = 0} .

(2)

Around the saturation density, ρ₀, the E₀(ρ) can be expanded, e.g., up to 3rd-order in density, as,

E_{0} (ρ) = E_{0} (ρ_{0}) + \frac{K_{0}}{2!} χ^{2} + \frac{J_{0}}{3!} χ^{3} + O (χ^{4}),

(3)

where χ =(ρ -ρ₀)/3ρ₀ is a dimensionless variable characterizing the deviations of the density from the saturation density, ρ₀. The first term E₀(ρ₀) on the right-hand-side of Eq. (3) is the binding energy per nucleon in SNM at ρ₀ and the coefficients of other terms are

K_{0} = {9 ρ_{0}^{2} \frac{d^{2} E_{0} (ρ)}{d ρ^{2}} |}_{ρ = ρ_{0}},

(4)

J_{0} = {27 ρ_{0}^{3} \frac{d^{3} E_{0} (ρ)}{d ρ^{3}} |}_{ρ = ρ_{0}},

(5)

where K₀ is the well-known incompressibility coefficient of SNM and J₀ is the skewness coefficient of SNM, i.e. the 3rd-order incompressibility coefficient of SNM ([24, 25]).

Similarly, one can expand the E_sym(ρ) around an arbitrary reference density, ρ_r, as

E_{sym} (ρ) = E_{sym} (ρ_{r}) + L (ρ_{r}) χ_{r} + O (χ_{r}^{2}),

(6)

with χ_r=(ρ -ρ_r)/3ρ_r, and the slope parameter of the symmetry energy at ρ_r is expressed as [33]

L (ρ_{r}) = {3 ρ_{r} \frac{d E_{sym} (ρ)}{d ρ} |}_{ρ = ρ_{r}} .

(7)

For ρ_r = ρ₀, the L(ρ_r) is reduced to the conventional slope parameter $L \equiv 3 ρ_{0} d E_{sym} (ρ) / d ρ |_{ρ = ρ_{0}}$ .

If δ and χ are assumed to be small quantities on the same order, nuclear matter bulk characteristic parameters can then be classified accordingly in different orders. For example, L and J₀ are on the same order-3, i.e., δ²χ for L and δ⁰χ³ for J₀. In this sense, E₀(ρ₀) is on the order-0, and, K₀ and E_sym(ρ₀) are on the order-2. To see the role of J₀ in the EOS of SNM, one can re-write Eq. (3) in a slightly different form as

E_{0} (ρ) ≃ E_{0} (ρ_{0}) + \frac{1}{2} K_{0} χ^{2} (1 + \frac{χ J_{0}}{3 K_{0}}) .

(8)

Assuming J₀ has roughly the same magnitude as K₀, one can see that the contribution from the J₀ term to the EOS of SNM becomes comparable with that from the K₀ term if the baryon density is larger than about 3ρ₀, corresponding to the typical densities inside a neutron star. On the other hand, the J₀ term plays a minor role for the EOS of SNM at subsaturation densities relevant for nuclear structure properties. As we will see later, the pressure of SNM at supra-saturation densities and the maximum mass of neutron stars indeed display strong sensitivity on the J₀ parameter.

2.2 Nuclear matter characteristic parameters in nonlinear RMF model

The nonlinear RMF model has made great success during the last decades in describing many nuclear phenomena (see, e.g., [34-46]). In the following, we briefly describe the nonlinear RMF model that we shall adopt in this work and present some useful expressions of nuclear matter characteristic parameters, especially the skewness coefficient, J₀. The interacting Lagrangian of the nonlinear RMF model supplemented with couplings between the isoscalar and the isovector mesons reads ([47-51])

\begin{array}{l} L = \bar{ψ} [γ_{μ} (i \partial^{μ} - g_{ω} ω^{μ} - g_{ρ} {\overset{⇀}{ρ}}^{μ} \cdot \overset{⇀}{τ}) - (M - g_{σ} σ)] ψ \\ - \frac{1}{2} m_{σ}^{2} σ^{2} + \frac{1}{2} \partial_{μ} σ \partial^{μ} σ - U (σ) \\ + \frac{1}{2} m_{ω}^{2} ω_{μ} ω^{μ} - \frac{1}{4} ω_{μ ν} ω^{μ ν} + \frac{1}{4} c_{ω} {(g_{ω} ω_{μ} ω^{μ})}^{2} \\ + \frac{1}{2} m_{ρ}^{2} {\overset{⇀}{ρ}}_{μ} \cdot {\overset{⇀}{ρ}}^{μ} - \frac{1}{4} {\overset{⇀}{ρ}}_{μ ν} \cdot {\overset{⇀}{ρ}}^{μ ν} \\ + \frac{1}{2} g_{ρ}^{2} {\overset{⇀}{ρ}}_{μ} \cdot {\overset{⇀}{ρ}}^{μ} Λ_{V} g_{ω}^{2} ω_{μ} ω^{μ}, \end{array}

(9)

where $ω_{μ ν} \equiv \partial_{μ} ω_{ν} - \partial_{ν} ω_{μ}$ and $ρ_{μ ν} \equiv \partial_{μ} {\overset{⇀}{ρ}}_{ν} - \partial_{ν} {\overset{⇀}{ρ}}_{μ}$ are are strength tensors for ω field and ρ field, respectively. ψ, σ, ωμ, ${\overset{⇀}{ρ}}_{μ}$ are nucleon field, isoscalar-scalar field, isoscalar-vector field, and isovector-vector field, respectively, and the arrows denote the vector in isospin space; U(σ)=bσM(gσσ)³/3+cσ(gσσ)⁴/4 is the self interaction term for σ field. Λ_V represents the coupling constant between the isovector ρ meson and the isoscalar ω meson, and it is important for the description of the density dependence of the symmetry energy. In addition, M is the nucleon mass and mσ, mω, mρ are masses of mesons.

In the mean field approximation, after neglecting effects of fluctuation and correlation, meson fields are replaced by their expectation values, i.e., $\bar{σ} \to σ$ , ${\bar{ω}}_{0} \to ω_{μ}$ , ${\bar{ρ}}_{0}^{(3)} \to {\overset{⇀}{ρ}}_{μ}$ , where subscript “0" indicates zeroth component of the four-vector, superscript “ (3)" indicates third component of the isospin. Furthermore, we also use in this work the non-sea approximation which neglects the effect due to negative energy states in the Dirac sea. The mean field equations are then expressed as

m_{σ}^{2} \bar{σ} = g σ [ρ_{S} - b_{σ} M {(g_{σ} \bar{σ})}^{2} - c_{σ} {(g_{σ} \bar{σ})}^{3}],

(10)

\begin{array}{l} m_{ω}^{2} {\bar{ω}}_{0} = g_{ω} [ρ - c_{ω} {(g_{ω} {\bar{ω}}_{0})}^{3} - Λ_{V} g_{ω} {\bar{ω}}_{0} {(g_{ρ} {\bar{ρ}}_{0}^{(3)})}^{2}], \end{array}

(11)

m_{ρ}^{2} {\bar{ρ}}_{0}^{(3)} = g_{ρ} [ρ_{p} - ρ_{n} - Λ_{V} g_{ρ} {\bar{ρ}}_{0}^{(3)} {(g_{ω} {\bar{ω}}_{0})}^{2}],

(12)

where

ρ = 〈 \bar{ψ} γ^{0} ψ 〉 = ρ_{n} + ρ_{p}, ρ_{S} = 〈 \bar{ψ} ψ 〉 = ρ_{S,n} + ρ_{S,p},

(13)

are the baryon density and scalar density, respectively, with the latter given by

\begin{array}{l} ρ_{S, J} = \frac{2}{{(2 π)}^{3}} \int_{0}^{k_{F}^{J}} d k \frac{M^{*}}{\sqrt{| k |^{2} + M^{*}^{2}}} \\ = \frac{M^{*}}{2 π^{2}} [k_{F}^{J} E_{F}^{J *} - M^{* 2} \ln (\frac{k_{F}^{J} + E_{F}^{J *}}{M^{*}})], J = p,n . \end{array}

(14)

In the above expression, we have $E_{F}^{J *} = \sqrt{k_{F}^{J 2} + M^{* 2}}$ and the nucleon Dirac mass is defined as

M^{*} \equiv M_{dirac}^{*} = M - g_{σ} \bar{σ} .

(15)

$k_{F}^{J} = k_{F} {(1 + τ_{3}^{J} δ)}^{1 / 3}$ is the Fermi momentum with $τ_{3}^{n} = + 1$ for neutrons and $τ_{3}^{p} = - 1$ for protons, and k_F=(3π²ρ /2)^1/3 is the Fermi momentum for SNM at ρ.

The energy-momentum density tensor for the interacting Lagrangian density in Eq. (9) can be written as

T^{μ ν} = \bar{ψ} i γ^{μ} \partial^{ν} ψ + \partial^{μ} σ \partial^{ν} σ - ω^{μ η} \partial^{ν} ω_{η} - ρ^{μ η} \partial^{ν} ρ_{η} - L g^{μ ν},

(16)

where gμν=(+,-,-,-) is the Minkowski metric. In the mean field approximation, the mean value of the time (zero) component of the energy-momentum density tensor is the energy density of the nuclear matter system, i.e.,

\begin{array}{l} ε = 〈 T^{00} 〉 \\ = ε_{kin}^{n} + ε_{kin}^{p} + \frac{1}{2} [m_{σ}^{2} {\bar{σ}}^{2} + m_{ω}^{2} {\bar{ω}}_{0}^{2} + m_{ρ}^{2} {({\bar{ρ}}_{0}^{(3)})}^{2}] \\ + \frac{1}{3} b_{σ} {(g_{σ} \bar{σ})}^{3} + \frac{1}{4} c_{σ} {(g_{σ} \bar{σ})}^{4} + \frac{3}{4} c_{ω} {(g_{ω} {\bar{ω}}_{0})}^{4} \\ + \frac{3}{2} {(g_{ρ} {\bar{ρ}}_{0}^{(3)})}^{2} Λ_{V} {(g_{ω} {\bar{ω}}_{0})}^{2}, \end{array}

(17)

where

\begin{array}{l} ε_{kin}^{J} = \frac{2}{{(2 π)}^{3}} \int_{0}^{k_{F}^{J}} d k \sqrt{| k |^{2} + M^{* 2}} \\ = \frac{1}{π^{2}} \int_{0}^{k_{F}^{J}} k^{2} d k \sqrt{k^{2} + M^{* 2}} \\ = \frac{1}{4} [3 E_{F}^{J *} ρ_{J} + M^{*} ρ_{S, J}], J = p,n, \end{array}

(18)

is the kinetic part of the energy density. Similarly, the mean value of space components of the energy-momentum density tensor corresponds to the pressure of the system, i.e.,

\begin{array}{l} P = \frac{1}{3} \sum_{j = 1}^{3} 〈 T^{j j} 〉 = P_{kin}^{n} + P_{kin}^{p} \\ - \frac{1}{2} [m_{σ}^{2} {\bar{σ}}^{2} - m_{ω}^{2} {\bar{ω}}_{0}^{2} - m_{ρ}^{2} {({\bar{ρ}}_{0}^{(3)})}^{2}] \\ - \frac{1}{3} b_{σ} {(g_{σ} \bar{σ})}^{3} - \frac{1}{4} c_{σ} {(g_{σ} \bar{σ})}^{4} + \frac{1}{4} c_{ω} {(g_{ω} {\bar{ω}}_{0})}^{4} \\ + \frac{1}{2} {(g_{ρ} {\bar{ρ}}_{0}^{(3)})}^{2} Λ_{V} {(g_{ω} {\bar{ω}}_{0})}^{2}, \end{array}

(19)

where the kinetic part of pressure is given by

P_{kin}^{J} = \frac{1}{3 π^{2}} \int_{0}^{k_{F}^{J}} d k \frac{k^{4}}{\sqrt{k^{2} + M^{*}^{2}}}, J = p,n .

(20)

The EOS of ANM can be calculated through the energy density, ε(ρ,δ), by

E (ρ, δ) = \frac{ε (ρ, δ)}{ρ} - M .

(21)

The EOS of SNM is just $E_{0} (ρ) \equiv E (ρ, δ = 0)$ , and the characteristic parameters K₀ and J₀ can be obtained from the following expressions

K_{0} (ρ) \equiv 9 ρ^{2} \frac{d^{2} E_{0}}{d ρ^{2}} = - \frac{9 ρ g_{σ}^{2} M^{* 2}}{Q_{σ} E_{F}^{* 2}} + \frac{9 ρ g_{ω}^{2}}{Q_{ω}} + \frac{3 k_{F}^{2}}{E_{F}^{*}} - 6 L_{0} (ρ),

(22)

\begin{array}{l} J_{0} (ρ) \equiv 27 ρ^{3} \frac{d^{3} E_{0}}{d ρ^{3}} = - \frac{3 k_{F}^{2}}{E_{F}^{*}} - \frac{3 k_{F}^{4}}{E_{F}^{* 3}} + \frac{27 g_{σ}^{2} M^{* 2} ρ^{2}}{Q_{σ} E_{F}^{* 3}} \times \\ (\frac{3 π^{2}}{2 k_{F} E_{F}^{*}} + \frac{2 g_{σ}^{2}}{Q_{σ}} - \frac{g_{σ} M^{*} η}{Q_{σ}^{2}} - \frac{2 g_{σ}^{2} M^{* 2}}{Q_{σ} E_{F}^{* 2}} + \frac{E_{F}^{*} ϕ}{2 Q_{σ}}) \\ - \frac{162 c_{ω} g_{ω}^{7} {\bar{ω}}_{0} ρ^{2}}{Q_{ω}^{3}} - 9 K_{0} (ρ), \end{array}

(23)

with

\begin{array}{l} L_{0} (ρ) \equiv 3 ρ \frac{d E_{0}}{d ρ} = 3 [\frac{E_{F}^{*}}{4} - \frac{M^{*} ρ_{S}}{4 ρ} + g_{ω} {\bar{ω}}_{0} \\ - \frac{1}{ρ} (\frac{1}{2} m_{σ}^{2} {\bar{σ}}^{2} + U (\bar{σ}) + \frac{1}{2} m_{ω}^{2} {\bar{ω}}_{0}^{2} + \frac{3}{4} c_{ω} g_{ω}^{4} {\bar{ω}}_{0}^{4})] . \end{array}

(24)

In the above expressions, we have

Q_{σ} = m_{σ}^{2} + g_{σ}^{2} (\frac{3 ρ_{S}}{M^{*}} - \frac{3 ρ}{E_{F}^{*}}) + 2 b_{σ} M g_{σ}^{3} \bar{σ} + 3 c_{σ} g_{σ}^{4} {\bar{σ}}^{2},

(25)

Q_{ω} = m_{ω}^{2} + 3 c_{ω} g_{ω}^{4} {\bar{ω}}_{0}^{2},

(26)

and

η = 3 g_{σ}^{3} (\frac{2 ρ_{S}}{M^{* 2}} - \frac{3 ρ}{M^{*} E_{F}^{*}} + \frac{M^{*} ρ}{E_{F}^{* 3}}) - 2 b_{σ} M g_{σ}^{3} - 6 c_{σ} g_{σ}^{4} \bar{σ},

(27)

ϕ = \frac{2 g_{σ}^{2} k_{F}^{2}}{E_{F}^{* 3}},

(28)

with $E_{F}^{*} = (k_{F}^{2} + M^{* 2})^{1 / 2}$ . To our best knowledge, Eq. (23) gives, for the first time [52], the analytical expression of the J₀ parameter in the nonlinear RMF model. In addition, we would like to point out that the general expression for E_sym(ρ) and L(ρ) in the nonlinear RMF model has been derived in Ref. [53].

3 Results and discussions

For the Lagrangian in Eq. (9), the properties of infinite nuclear matter is uniquely determined by fσ=gσ, bσ, cσ, fω=gω/mω, cω, fρ=gρ/mρ, Λ_V, and M. If the nucleon mass in vacuum is set to be M = 939 MeV, one then has seven total parameters to determine the properties of infinite nuclear matter in the nonlinear RMF model. Following the correlation analysis method proposed in Ref. [54] within the Skyrme-Hartree-Fock (SHF) approach, instead of directly using the seven microscopic parameters, i.e., fσ, bσ, cσ, fω, cω, fρ, and Λ_V, one can determine their values explicitly in terms of seven macroscopic quantities, i.e., ρ₀, E₀(ρ₀), K₀, J₀, $M_{dirac}^{* 0} \equiv M_{dirac}^{*} (ρ_{0})$ , E_sym(ρ_c), and L(ρ_c) where ρ_c is the cross density whose value is fixed in this work at 0.11 fm^-3 ([33]). Then, by varying individually these macroscopic quantities within their known ranges, one can examine transparently the correlation of nuclear matter properties with each individual macroscopic quantity. Recently, this simple correlation analysis method has been successfully applied to study the neutron skin ([54, 33]) and giant resonances of finite nuclei ([30, 55]), the higher-order bulk characteristic parameters of ANM ([25]), and the relationship between the nuclear matter symmetry energy and the symmetry energy coefficient in the mass formula ([56]). We would like to point out although the seven macroscopic quantities defined above coherently act on the maximum mass of neutron stars and the pressure of the SNM, they are independent with each other in our analysis since we vary one quantity by keeping other six quantities fixed. This is one of the main advantages of our approach since the physics of these macroscopic quantities is different.

To examine the correlation of pressure of SNM at supra-saturation densities with each macroscopic quantity, we show in Fig. 1 the pressure of SNM P(ρ) at ρ=3ρ₀ from the nonlinear RMF model based on the FSUGold interaction ([49]) by varying individually ρ₀, E₀(ρ₀), $M_{dirac}^{* 0}$ , K₀, and J₀ within their empirical uncertain ranges, namely, varying one quantity at a time while keeping all others at their default values in FSUGold for which we have ρ₀=0.148 fm^-3, E₀(ρ₀)=-16.3 MeV, $M_{dirac}^{* 0} = 0.61 M$ , K₀=230 MeV, J₀=-522.6 MeV, E_sym(ρ_c)=27.11 MeV, and L(ρ_c)=49.97 MeV. It should be mentioned that the pressure of SNM is independent of the values of E_sym(ρ_c) and L(ρ_c). It is seen from Fig. 1 that the pressure of SNM P(ρ) at ρ=3ρ₀ increases with ρ₀, $M_{dirac}^{* 0}$ , K₀, and J₀ while it decreases with E₀(ρ₀). In particular, the pressure of SNM P(ρ) at ρ=3ρ₀ displays a specially strong correlation with J₀. We note that the pressure of SNM P(ρ) at other supra-saturation densities exhibits similar correlations with ρ₀, E₀(ρ₀), $M_{dirac}^{* 0}$ , K₀, and J₀. These features indicate that the pressure of SNM at supra-saturation densities is sensitive to the J₀ value, and thus the experimental constraints on the pressure of SNM at supra-saturation densities may provide important information on the J₀ value.

Figure 1:

Pressure of SNM at ρ=3ρ₀ from the nonlinear RMF model based on the FSUGold interaction by varying individually ρ₀ (a), E₀(ρ₀) (b),

M_{dirac}^{* 0}

(c), K₀ (d), and J₀ (e).

Since the pressure of SNM at supra-saturation densities is sensitive to the J₀ value, the maximum mass, M_max, of static neutron stars is also expected to be sensitive to the J₀ value. The mass and radius of static neutron stars can be obtained from solving the Tolman-Oppenheimer-Volkoff (TOV) equations with a given neutron star matter EOS. A neutron star generally contains core, inner crust, and outer crust from the center to surface. In this work, for the core where the baryon density is larger than the core-curst transition density, ρ_t, we use the EOS of β-stable and charge neutral, npeμ matter obtained from the nonlinear RMF model. In the inner crust with densities between ρ_out and ρ_t where the nuclear pastas may exist, we construct its EOS (pressure, P, as a function of energy density, 𝒠) according to $P = a + b E^{4 / 3}$ because of our poor knowledge about its EOS from first principle ([57, 58]). The ρ_out=2.46× 10^-4 fm^-3 is the density separating the inner from the outer crust. The constants a and b are then easily determined by the pressure and energy density at ρ_t and ρ_out [58]. In this work, the ρ_t is determined self-consistently within the nonlinear RMF model using the thermodynamical method (see, e.g., [51] for the details). In the outer crust with 6.93× 10^-13 fm^-3<ρ <ρ_out, we use the EOS of BPS ([59, 60]), and in the region of 4.73×10^-15 fm^-3<ρ <6.93× 10^-13 fm^-3, we use the EOS of FMT ([59]).

Similarly, as in Fig. 1, we plot in Fig. 2 the maximum mass, M_max, of static neutron stars from the nonlinear RMF model based on the FSUGold interaction by varying individually ρ₀, E₀(ρ₀), $M_{dirac}^{* 0}$ , K₀, J₀, E_sym(ρ_c), and L(ρ_c) within their empirical uncertain ranges. Indeed, one can see that the M_max displays a very strong positive correlation with the J₀ parameter. In addition, the M_max exhibits weak positive correlation with the $M_{dirac}^{* 0}$ and K₀, and weak negative correlation with the ρ₀ and E₀(ρ₀). It is interesting to see that the M_max is essentially independent of the values of E_sym(ρ_c) and L(ρ_c), implying that, in the nonlinear RMF model, the M_max is basically determined by the isoscalar part of the nuclear matter EOS. Since the seven microscopic parameters change with the macroscopic quantities, it is thus not surprising to see that the maximum mass of a neutron star based on the FSUGold interaction by varying macroscopic quantities may be totally different from the default one from FSUGold, which is about 1.74M_⊙. These features indicate that the observed largest mass of neutron stars may put important constraint on the J₀ value. We would like to point out that the interaction FSUGold is only used in Figs. 1 and 2 as a reference for the correlation analyses and using other RMF interactions will not change our conclusion.

Figure 2:

(Color online) Maximum mass of static neutron stars from the nonlinear RMF model based on the FSUGold interaction by varying individually ρ₀ (a), E₀(ρ₀) (b),

M_{dirac}^{* 0}

(c), K₀ (d), J₀ (e), E_sym(ρ_c) (f), and L(ρ_c) (g).

Experimentally, the pressure of SNM at supra-saturation densities (from 2ρ₀ to about 5ρ₀) has been constrained by measurements of collective flows in HIC ([3]), which is shown as a band in the left window of Fig. 3. In the nonlinear RMF model, if one only changes the J₀ value while the other 6 macroscopic quantities are kept at their values in the FSUGold interaction, one can find that the J₀ value should be in the range of -985 MeV ≤ J₀≤ -327 MeV to be consistent with the flow data in HIC ([3]). However, keeping the other 6 macroscopic quantities at their values in the FSUGold interaction is obviously a strong assumption because the extraction of the J₀ value from the flow data in HIC will also depends on the values of ρ₀, E₀(ρ₀), $M_{dirac}^{* 0}$ , and K₀, which can be varied within their empirical uncertain ranges. For the nonlinear RMF model, we use, in this work, the following empirical uncertain ranges for these macroscopic quantities, i.e., ρ₀=0.153±0.008 fm^-3, E₀(ρ₀)=-16.2±0.3 MeV, $M_{dirac}^{* 0} / M = 0.61 \pm 0.04$ , and K₀=230±20 MeV, which represent the typical uncertain ranges known or predicted from different interactions in the nonlinear RMF model ([50]).

Figure 3:

(Color online) Left window: Pressure of SNM as a function of baryon density. The solid (dashed) line is the prediction from the nonlinear RMF model with J₀=-10 (-1280) MeV and the set “S (H)" for ρ₀, E₀(ρ₀),

M_{dirac}^{* 0}

, and K₀. The band represents the constraints from flow data in HIC [3]. Right window: The maximum mass of static neutron stars as a function of J₀ in the nonlinear RMF model with the set “NS-H" for ρ₀, E₀(ρ₀),

M_{dirac}^{* 0}

, K₀, E_sym(ρ_c), and L(ρ_c). The band represents mass 2.01± 0.04M_⊙ for PSR J0348+0432 [61].

Based on the pressure of SNM constrained by flow data in HIC ([3]), to extract the upper limit of the J₀ value, one should use the values of ρ₀, E₀(ρ₀), $M_{dirac}^{* 0}$ , and K₀ that make the resulting pressure of SNM as small as possible when J₀ is fixed. This can be obtained by using ρ₀=0.145 fm^-3, E₀(ρ₀)=-15.9 MeV, $M_{dirac}^{*} / M = 0.57$ , and K₀=210 MeV, denoted as set “S", since the pressure of SNM P(ρ / ρ₀) at supra-saturation densities increases with ρ₀, $M_{dirac}^{* 0}$ , and K₀ while decreases with E₀(ρ₀), as shown in Fig. 1. With the set “S" for ρ₀, E₀(ρ₀), $M_{dirac}^{* 0}$ , and K₀, one can find the upper limit of J₀= -10 MeV for the J₀ value, which is indicated by a solid line in the left window of Fig. 3. For J₀ > -10 MeV, the model would over-predict the pressure of SNM constrained by flow data in HIC ([3]). Similarly, one can obtain the lower limit of the J₀ value by using the values of ρ₀, E₀(ρ₀), $M_{dirac}^{* 0}$ , and K₀ that make the resulting pressure of SNM as large as possible when J₀ is fixed, and this can be obtained with ρ₀=0.161 fm^-3, E₀(ρ₀)=-16.5 MeV, $M_{dirac}^{* 0} / M = 0.65$ , and K₀=250 MeV, denoted as set “H". Using the set “H" for ρ₀, E₀(ρ₀), $M_{dirac}^{* 0}$ , and K₀, one can extract the lower limit of J₀= -1280 MeV, which is indicated by dashed line in the left window of Fig. 3. The model would under-predict the pressure of SNM constrained by flow data in HIC ([3]) if J₀ < -1280 MeV. Therefore, from the pressure of SNM constrained by flow data in HIC ([3]), one can extract the constraint of -1280 MeV ≤ J₀ ≤ -10 MeV.

Recently, a new neutron star, PSR J0348+0432, with a mass of 2.01±0.04M_⊙ was discovered ([61]), and this neutron star is only the second pulsar with a precisely determined mass around 2M_⊙ after PSR J1614-2230 ([62]) and sets a new record of the maximum mass of neutron stars. The lower mass limit of 1.97M_⊙ for PSR J0348+0432 thus may set a lower limit of the J₀ value below which the model cannot predict a neutron star with mass equal or above 1.97M_⊙. To extract the lower limit of the J₀ value from the observed heaviest neutron star, PSR J0348+0432, one can use the values of ρ₀, E₀(ρ₀), $M_{dirac}^{* 0}$ , K₀, E_sym(ρ_c), and L(ρ_c) that make the resulting maximum mass of neutron stars as large as possible when J₀ is fixed, and from Fig. 2 this can be obtained with ρ₀=0.145 fm^-3, E₀(ρ₀)=-16.5 MeV, K₀=250 MeV, $M_{dirac}^{* 0} / M = 0.65$ , E_sym(ρ_c)=26 MeV, and L(ρ_c)=60 MeV, denoted as set “NS-H". This leads to a lower limit of J₀= -494 MeV for the J₀ value as shown in the right window of Fig. 3 where the maximum mass of neutron stars is plotted as a function of J₀ when the other 6 macroscopic quantities are fixed at their values as in set “NS-H". For J₀ < -494 MeV, the maximum mass of static neutron stars predicted in the nonlinear RMF model would be always smaller than 1.97M_⊙. It should be pointed out that here the interior of neutron stars has been assumed to be npeμ matter. New degrees of freedom, such as hyperons or/and quark matter that could be present in the interior of neutron stars, usually soften the EOS of neutron star matter and thus a larger J₀ value would be necessary to obtain a neutron star with mass of 1.97M_⊙. Therefore, including the new degrees of freedom in neutron stars will be consistent with the constraint of J₀≥ -494 MeV.

Combining the constraint of -1280 MeV ≤ J₀ ≤ -10 MeV from the pressure of SNM constrained by flow data in HIC ([3]) which favors a smaller J₀ value and the constraint of J₀≥ -494 MeV from the recently discovered heaviest neutron star PSR J0348+0432 ([61]) which favors a larger J₀ value, one can extract the following constraint for the J₀ parameter

- 494 MeV \leq J_{0} \leq - 10 MeV .

(29)

It should be emphasized that the constraint -494 MeV≤J₀≤-10 MeV represents a conservative extraction based on flow data in HIC ([3]) and the recently discovered heaviest neutron star, PSR J0348+0432 ([61]) in the nonlinear RMF model. This is because we have not considered the possible correlations that existed among the ρ₀, E₀(ρ₀), $M_{dirac}^{* 0}$ , K₀, E_sym(ρ_c), and L(ρ_c), and have simply set simultaneously their values in the boundary of their empirical uncertain ranges. Considering the correlations possibly existed among the ρ₀, E₀(ρ₀), $M_{dirac}^{* 0}$ , K₀, E_sym(ρ_c) and L(ρ_c) should further narrow the constraint -494 MeV≤J₀≤-10 MeV, and it will be interesting to see the quantitative constraint on the J₀ parameter based on the data of finite nuclei, neutron stars, and heavy-ion collisions using the exhaustive statistical analysis method although this is beyond the scope of the present work. The conservative constraint -494 MeV≤J₀≤-10 MeV obtained in the present work indicates that, if the J₀ value is out of the region -494 MeV≤J₀≤-10 MeV, the nonlinear RMF model either cannot predict the pressure of SNM constrained by flow data in HIC ([3]) or cannot describe the recently discovered heaviest neutron star PSR J0348+0432 ([61]). It is worth mentioning that the constraint on the J₀ depends on our knowledge of the other six quantities. Any improvement on these six macroscopic quantities will make the range for the J₀ constraint narrower. In addition, the extracted constraint on the J₀ could depend on the form of the energy density functional, and it will be interesting to see how the constraint changes if other energy density functionals are used.

In Fig. 4, we show the comparison of the J₀ constraint obtained in our analysis with those obtained with other analyses and/or other methods ([22, 25, 63, 64]), including the constraint of J₀ = -700 ± 500 MeV obtained by Ref. [63] from the analysis of nuclear GMR, the constraint of $J_{0} = - 280_{- 410}^{+ 72}$ ( $- 500_{- 290}^{+ 170}$ ) MeV obtained by Ref. [64] from analyzing a heterogeneous data set of six neutron stars using a Markov chain Monte Carlo algorithm within a Bayesian framework by assuming r_ph>> R (r_ph= R) where r_ph is the photospheric radius at the time the flux is evaluated and R is the stellar radius, the constraint of J₀ = -355 ± 95 MeV deduced by Ref. [25] based on a correlation analysis method within SHF energy density functional, and the constraint of J₀ = -390 ± 90 MeV deduced by Ref. [22] who used a similar method as Ref. [25]. It is seen that the constrained region of J₀ obtained in the present work has a remarkable overlap with those existing in the literature. In particular, our present constraint from the relativistic model is nicely consistent with the constraints deduced from the non-relativistic SHF approach, and they all indicate that the J₀ parameter should be larger than about -500 MeV.

Figure 4:

(Color online) Comparison between the constraint of J₀ extracted in the present work and those from the existing literature [22, 25, 63, 64].

4 Summary

Within the nonlinear relativistic mean field model, using macroscopic nuclear matter characteristic parameters instead of the microscopic coupling constants as direct input quantities, we have demonstrated that the pressure of symmetric nuclear matter at supra-saturation densities and the maximum mass of neutron stars provide useful probes for the skewness coefficient, J₀, of symmetric nuclear matter. In particular, using the existing experimental constraints on the pressure of symmetric nuclear matter at supra-saturation densities from flow data in heavy ion collisions and the astrophysical observation of recently discovered heaviest neutron star PSR J0348+0432, with the former requiring a smaller J₀ while the latter requires a larger J₀, we have extracted a constraint of -494 MeV≤J₀≤-10 MeV.

We have compared the present constraint with the results obtained in other analyses, and found they are nicely in agreement. In particular, our present constraint from the relativistic model is nicely consistent with the constraints deduced from the non-relativistic Skyrme-Hartree-Fock approach, and they all indicate that the J₀ parameter cannot be too small, namely, it should be larger than about -500 MeV. The present constraint on the J₀ parameter provides important information on the high density behaviors of the EOS of symmetric nuclear matter and also may be potentially useful for the determination of the high density behaviors of the EOS of asymmetric nuclear matter, especially the high density symmetry energy.

References

[1]

J.P. Blaizot,

Nuclear compressibilities

. Phys. Rep. 64, 171-248 (1980). doi: 10.1016/0370-1573(80)90001-0