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

Constraints on the skewness coefficient of symmetric nuclear matter within the nonlinear relativistic mean field model 1 first-author 1 corresp Corresponding author, lwchen@sjtu.edu.cn Corresponding author, lwchen@sjtu.edu.cn lwchen@sjtu.edu.cn School of Physics and Astronomy and Shanghai Key Laboratory for Particle Physics and Cosmology, Shanghai Jiao Tong University, Shanghai 200240, China School of Physics and Astronomy and Shanghai Key Laboratory for Particle Physics and Cosmology, Shanghai Jiao Tong University, Shanghai 200240, China 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 0 , 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 0 while the latter favors a larger J 0 , we extract a constraint of -494 MeV≤ J 0 ≤-10 MeV based on the NL-RMF model. This constraint is compared with the results obtained in other analyses. Equation of state of nuclear matter Heavy-ion collisions Neutron stars 1 Introduction 1 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 1 - 16 16 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 17 - 23 23 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, ρ 0 , of symmetric nuclear matter (SNM) are usually introduced to quantitatively characterize the EOS of ANM. For example, the energy, E 0 ( ρ 0 ), and incompressibility, K 0 , of SNM, as well as the symmetry energy, E sym ( ρ 0 ), and its slope parameter, L , are the four famous lower-order bulk characteristic parameters of EOS of ANM. These bulk parameters defined at ρ 0 provide important information on both sub- and supra-saturation density behaviors of the EOS of ANM ([ 24 24 , 25 25 ]). Based on the empirical liquid-drop-like model analyses of high precision data about nuclear masses, the E 0 ( ρ 0 ) is well known to be about -16 MeV. The incompressibility has been determined to be K 0 = 240 ± 40 MeV from analyzing experimental data of nuclear giant monopole resonances (GMR) ([ 1 1 , 26 26 - 30 30 26-30 ]). For E sym ( ρ 0 ) and L , the existing constraints extracted from terrestrial laboratory measurements and astrophysical observations are found to be essentially consistent with E sym ( ρ 0 ) = 32.5 ± 2.5 MeV and L = 55 ± 25 MeV (see, e.g., Refs. [ 31 31 , 32 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 0 ( ρ 0 ), K 0 , E sym ( ρ 0 ), and L , the next bulk characteristic parameter should be the skewness coefficient, J 0 , (also denoted as K ’ or Q 0 in some literature) of SNM, which is related to the third-order density derivative of the binding energy per nucleon of SNM at ρ 0 . The higher-order bulk characteristic parameter J 0 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 0 parameter, and it is thus of great interest and critical importance to constrain the J 0 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 0 . 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 0 parameter. 2 The skewness coefficient J 0 in nonlinear RMF model 1 2.1 Nuclear matter characteristic parameters 2 The EOS of isospin asymmetric nuclear matter, namely E ( ρ , δ ), can be expanded as a power series of even-order terms in δ as <math display="block" id="M1"><mrow><mi>E</mi><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>δ</mi><mo stretchy="false">)</mo><mo>≃</mo><msub><mi>E</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mi>E</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><msup><mi>δ</mi><mn>2</mn></msup><mo>+</mo><mi mathvariant="script">O</mi><mo stretchy="false">(</mo><msup><mi>δ</mi><mn>4</mn></msup><mo stretchy="false">)</mo><mn>,</mn></mrow></math> where E 0 ( ρ )= E ( ρ , δ =0) is the EOS of symmetric nuclear matter, and the symmetry energy is expressed as <math display="block" id="M2"><mrow><msub><mi>E</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mrow><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><mi>E</mi><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>δ</mi><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><msup><mi>δ</mi><mn>2</mn></msup></mrow></mfrac></mrow><mo>|</mo></mrow></mrow><mrow><mi>δ</mi><mo>=</mo><mn>0</mn></mrow></msub><mn>.</mn></mrow></math> Around the saturation density, ρ 0 , the E 0 ( ρ ) can be expanded, e.g., up to 3rd-order in density, as, <math display="block" id="M3"><mrow><msub><mi>E</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mi>E</mi><mn>0</mn></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mn>0</mn></msub><mo stretchy="false">)</mo><mo>+</mo><mfrac><mrow><msub><mi>K</mi><mn>0</mn></msub></mrow><mrow><mn>2</mn><mo>!</mo></mrow></mfrac><msup><mi>χ</mi><mn>2</mn></msup><mo>+</mo><mfrac><mrow><msub><mi>J</mi><mn>0</mn></msub></mrow><mrow><mn>3</mn><mo>!</mo></mrow></mfrac><msup><mi>χ</mi><mn>3</mn></msup><mo>+</mo><mi mathvariant="script">O</mi><mo stretchy="false">(</mo><msup><mi>χ</mi><mn>4</mn></msup><mo stretchy="false">)</mo><mn>,</mn></mrow></math> where χ =( ρ - ρ 0 )/3 ρ 0 is a dimensionless variable characterizing the deviations of the density from the saturation density, ρ 0 . The first term E 0 ( ρ 0 ) on the right-hand-side of Eq. (3) is the binding energy per nucleon in SNM at ρ 0 and the coefficients of other terms are <math display="block" id="M4"><mrow><msub><mi>K</mi><mn>0</mn></msub><mo>=</mo><msub><mrow><mrow><mrow><mn>9</mn><msubsup><mi>ρ</mi><mn>0</mn><mn>2</mn></msubsup><mfrac><mrow><msup><mtext>d</mtext><mn>2</mn></msup><msub><mi>E</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>d</mtext><msup><mi>ρ</mi><mn>2</mn></msup></mrow></mfrac></mrow><mo>|</mo></mrow></mrow><mrow><mi>ρ</mi><mo>=</mo><msub><mi>ρ</mi><mn>0</mn></msub></mrow></msub><mn>,</mn><mi> </mi><mi> </mi></mrow></math> <math display="block" id="M5"><mrow><msub><mi>J</mi><mn>0</mn></msub><mo>=</mo><msub><mrow><mrow><mrow><mn>27</mn><msubsup><mi>ρ</mi><mn>0</mn><mn>3</mn></msubsup><mfrac><mrow><msup><mtext>d</mtext><mn>3</mn></msup><msub><mi>E</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>d</mtext><msup><mi>ρ</mi><mn>3</mn></msup></mrow></mfrac></mrow><mo>|</mo></mrow></mrow><mrow><mi>ρ</mi><mo>=</mo><msub><mi>ρ</mi><mn>0</mn></msub></mrow></msub><mn>,</mn></mrow></math> where K 0 is the well-known incompressibility coefficient of SNM and J 0 is the skewness coefficient of SNM, i.e. the 3rd-order incompressibility coefficient of SNM ([ 24 24 , 25 25 ]). Similarly, one can expand the E sym ( ρ ) around an arbitrary reference density, ρ r , as <math display="block" id="M6"><mrow><msub><mi>E</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mi>E</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>r</mtext></msub><mo stretchy="false">)</mo><mo>+</mo><mi>L</mi><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>r</mtext></msub><mo stretchy="false">)</mo><msub><mi>χ</mi><mtext>r</mtext></msub><mo>+</mo><mi mathvariant="script">O</mi><mo stretchy="false">(</mo><msubsup><mi>χ</mi><mtext>r</mtext><mn>2</mn></msubsup><mo stretchy="false">)</mo><mn>,</mn></mrow></math> with χ r =( ρ - ρ r )/3 ρ r , and the slope parameter of the symmetry energy at ρ r is expressed as [ 33 33 ] <math display="block" id="M7"><mrow><mi>L</mi><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>r</mtext></msub><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mrow><mrow><mn>3</mn><msub><mi>ρ</mi><mtext>r</mtext></msub><mfrac><mrow><mtext>d</mtext><msub><mi>E</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>d</mtext><mi>ρ</mi></mrow></mfrac></mrow><mo>|</mo></mrow></mrow><mrow><mi>ρ</mi><mo>=</mo><msub><mi>ρ</mi><mtext>r</mtext></msub></mrow></msub><mn>.</mn></mrow></math> For ρ r = ρ 0 , the L ( ρ r ) is reduced to the conventional slope parameter <math id="M8"><mrow><mi>L</mi><mo>≡</mo><mn>3</mn><msub><mi>ρ</mi><mn>0</mn></msub><mtext>d</mtext><msub><mi>E</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>/</mo><mtext>d</mtext><mi>ρ</mi><msub><mo>|</mo><mrow><mi>ρ</mi><mo>=</mo><msub><mi>ρ</mi><mn>0</mn></msub></mrow></msub></mrow></math> . 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 0 are on the same order-3, i.e., δ 2 χ for L and δ 0 χ 3 for J 0 . In this sense, E 0 ( ρ 0 ) is on the order-0, and, K 0 and E sym ( ρ 0 ) are on the order-2. To see the role of J 0 in the EOS of SNM, one can re-write Eq. (3) in a slightly different form as <math display="block" id="M9"><mrow><msub><mi>E</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>≃</mo><msub><mi>E</mi><mn>0</mn></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mn>0</mn></msub><mo stretchy="false">)</mo><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>K</mi><mn>0</mn></msub><msup><mi>χ</mi><mn>2</mn></msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mfrac><mrow><mi>χ</mi><msub><mi>J</mi><mn>0</mn></msub></mrow><mrow><mn>3</mn><msub><mi>K</mi><mn>0</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mn>.</mn></mrow></math> Assuming J 0 has roughly the same magnitude as K 0 , one can see that the contribution from the J 0 term to the EOS of SNM becomes comparable with that from the K 0 term if the baryon density is larger than about 3 ρ 0 , corresponding to the typical densities inside a neutron star. On the other hand, the J 0 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 0 parameter. 2.2 Nuclear matter characteristic parameters in nonlinear RMF model 2 The nonlinear RMF model has made great success during the last decades in describing many nuclear phenomena (see, e.g., [ 34 34 - 46 46 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 0 . The interacting Lagrangian of the nonlinear RMF model supplemented with couplings between the isoscalar and the isovector mesons reads ([ 47 47 - 51 51 47-51 ]) <math display="block" id="M10"><mtable columnalign="left"><mtr><mtd><mi mathvariant="script">L</mi><mo>=</mo><mover accent="true"><mi>ψ</mi><mo stretchy="true">¯</mo></mover><mrow><mo>[</mo><mrow><msub><mi>γ</mi><mi>μ</mi></msub><mo stretchy="false">(</mo><mi>i</mi><msup><mo>∂</mo><mi>μ</mi></msup><mo>−</mo><msub><mi>g</mi><mi>ω</mi></msub><msup><mi>ω</mi><mi>μ</mi></msup><mo>−</mo><msub><mi>g</mi><mi>ρ</mi></msub><msup><mover accent="true"><mi>ρ</mi><mo>⇀</mo></mover><mi>μ</mi></msup><mo>⋅</mo><mover accent="true"><mi>τ</mi><mo>⇀</mo></mover><mo stretchy="false">)</mo><mo>−</mo><mo stretchy="false">(</mo><mi>M</mi><mo>−</mo><msub><mi>g</mi><mi>σ</mi></msub><mi>σ</mi><mo stretchy="false">)</mo></mrow><mo>]</mo></mrow><mi>ψ</mi></mtd></mtr><mtr><mtd><mo>−</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msubsup><mi>m</mi><mi>σ</mi><mn>2</mn></msubsup><msup><mi>σ</mi><mn>2</mn></msup><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mo>∂</mo><mi>μ</mi></msub><mi>σ</mi><msup><mo>∂</mo><mi>μ</mi></msup><mi>σ</mi><mo>−</mo><mi>U</mi><mo stretchy="false">(</mo><mi>σ</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msubsup><mi>m</mi><mi>ω</mi><mn>2</mn></msubsup><msub><mi>ω</mi><mi>μ</mi></msub><msup><mi>ω</mi><mi>μ</mi></msup><mo>−</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><msub><mi>ω</mi><mrow><mi>μ</mi><mi>ν</mi></mrow></msub><msup><mi>ω</mi><mrow><mi>μ</mi><mi>ν</mi></mrow></msup><mo>+</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><msub><mi>c</mi><mi>ω</mi></msub><msup><mrow><mo>(</mo><mrow><msub><mi>g</mi><mi>ω</mi></msub><msub><mi>ω</mi><mi>μ</mi></msub><msup><mi>ω</mi><mi>μ</mi></msup></mrow><mo>)</mo></mrow><mn>2</mn></msup></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msubsup><mi>m</mi><mi>ρ</mi><mn>2</mn></msubsup><msub><mover accent="true"><mi>ρ</mi><mo>⇀</mo></mover><mi>μ</mi></msub><mo>⋅</mo><msup><mover accent="true"><mi>ρ</mi><mo>⇀</mo></mover><mi>μ</mi></msup><mo>−</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><msub><mover accent="true"><mi>ρ</mi><mo>⇀</mo></mover><mrow><mi>μ</mi><mi>ν</mi></mrow></msub><mo>⋅</mo><msup><mover accent="true"><mi>ρ</mi><mo>⇀</mo></mover><mrow><mi>μ</mi><mi>ν</mi></mrow></msup></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msubsup><mi>g</mi><mi>ρ</mi><mn>2</mn></msubsup><msub><mover accent="true"><mi>ρ</mi><mo>⇀</mo></mover><mi>μ</mi></msub><mo>⋅</mo><msup><mover accent="true"><mi>ρ</mi><mo>⇀</mo></mover><mi>μ</mi></msup><msub><mi>Λ</mi><mtext>V</mtext></msub><msubsup><mi>g</mi><mi>ω</mi><mn>2</mn></msubsup><msub><mi>ω</mi><mi>μ</mi></msub><msup><mi>ω</mi><mi>μ</mi></msup><mn>,</mn></mtd></mtr></mtable></math> where <math id="M11"><mrow><msub><mi>ω</mi><mrow><mi>μ</mi><mi>ν</mi></mrow></msub><mo>≡</mo><msub><mo>∂</mo><mi>μ</mi></msub><msub><mi>ω</mi><mi>ν</mi></msub><mo>−</mo><msub><mo>∂</mo><mi>ν</mi></msub><msub><mi>ω</mi><mi>μ</mi></msub></mrow></math> and <math id="M12"><mrow><msub><mi>ρ</mi><mrow><mi>μ</mi><mi>ν</mi></mrow></msub><mo>≡</mo><msub><mo>∂</mo><mi>μ</mi></msub><msub><mover accent="true"><mi>ρ</mi><mo>⇀</mo></mover><mi>ν</mi></msub><mo>−</mo><msub><mo>∂</mo><mi>ν</mi></msub><msub><mover accent="true"><mi>ρ</mi><mo>⇀</mo></mover><mi>μ</mi></msub></mrow></math> are are strength tensors for ω field and ρ field, respectively. ψ , σ , ω μ , <math id="M13"><mrow><msub><mover accent="true"><mi>ρ</mi><mo>⇀</mo></mover><mi>μ</mi></msub></mrow></math> 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 /3+ c σ ( g σ σ ) 4 /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., <math id="M14"><mrow><mover accent="true"><mi>σ</mi><mo>¯</mo></mover><mo>→</mo><mi>σ</mi></mrow></math> , <math id="M15"><mrow><msub><mover accent="true"><mi>ω</mi><mo>¯</mo></mover><mn>0</mn></msub><mo>→</mo><msub><mi>ω</mi><mi>μ</mi></msub></mrow></math> , <math id="M16"><mrow><msubsup><mover accent="true"><mi>ρ</mi><mo>¯</mo></mover><mn>0</mn><mrow><mo stretchy="false">(</mo><mn>3</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>→</mo><msub><mover accent="true"><mi>ρ</mi><mo>⇀</mo></mover><mi>μ</mi></msub></mrow></math> , 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 <math display="block" id="M17"><mrow><msubsup><mi>m</mi><mi>σ</mi><mn>2</mn></msubsup><mover accent="true"><mi>σ</mi><mo>¯</mo></mover><mo>=</mo><mi>g</mi><mi>σ</mi><mrow><mo>[</mo><mrow><msub><mi>ρ</mi><mtext>S</mtext></msub><mo>−</mo><msub><mi>b</mi><mi>σ</mi></msub><mi>M</mi><msup><mrow><mrow><mo>(</mo><mrow><msub><mi>g</mi><mi>σ</mi></msub><mover accent="true"><mi>σ</mi><mo>¯</mo></mover></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><mo>−</mo><msub><mi>c</mi><mi>σ</mi></msub><msup><mrow><mo stretchy="false">(</mo><msub><mi>g</mi><mi>σ</mi></msub><mover accent="true"><mi>σ</mi><mo>¯</mo></mover><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow><mo>]</mo></mrow><mo>,</mo></mrow></math> <math display="block" id="M18"><mtable columnalign="left"><mtr><mtd><msubsup><mi>m</mi><mi>ω</mi><mn>2</mn></msubsup><msub><mover accent="true"><mi>ω</mi><mo>¯</mo></mover><mn>0</mn></msub><mo>=</mo><msub><mi>g</mi><mi>ω</mi></msub><mrow><mo>[</mo><mrow><mi>ρ</mi><mo>−</mo><msub><mi>c</mi><mi>ω</mi></msub><msup><mrow><mrow><mo>(</mo><mrow><msub><mi>g</mi><mi>ω</mi></msub><msub><mrow><mover accent="true"><mi>ω</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn></msub></mrow><mo>)</mo></mrow></mrow><mn>3</mn></msup><mo>−</mo><msub><mi>Λ</mi><mi>V</mi></msub><msub><mi>g</mi><mi>ω</mi></msub><msub><mrow><mover accent="true"><mi>ω</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn></msub><msup><mrow><mrow><mo>(</mo><mrow><msub><mi>g</mi><mi>ρ</mi></msub><msubsup><mrow><mover accent="true"><mi>ρ</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn><mrow><mo stretchy="false">(</mo><mn>3</mn><mo stretchy="false">)</mo></mrow></msubsup></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup></mrow><mo>]</mo></mrow><mn>,</mn></mtd></mtr><mtr><mtd><mrow/></mtd></mtr></mtable></math> <math display="block" id="M19"><mrow><msubsup><mi>m</mi><mi>ρ</mi><mn>2</mn></msubsup><msubsup><mrow><mover accent="true"><mi>ρ</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn><mrow><mo stretchy="false">(</mo><mn>3</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><msub><mi>g</mi><mi>ρ</mi></msub><mrow><mo>[</mo><mrow><msub><mi>ρ</mi><mtext>p</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mtext>n</mtext></msub><mo>−</mo><msub><mi>Λ</mi><mtext>V</mtext></msub><msub><mi>g</mi><mi>ρ</mi></msub><msubsup><mrow><mover accent="true"><mi>ρ</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn><mrow><mo stretchy="false">(</mo><mn>3</mn><mo stretchy="false">)</mo></mrow></msubsup><msup><mrow><mrow><mo>(</mo><mrow><msub><mi>g</mi><mi>ω</mi></msub><msub><mrow><mover accent="true"><mi>ω</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn></msub></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup></mrow><mo>]</mo></mrow><mn>,</mn></mrow></math> where <math display="block" id="M20"><mrow><mi>ρ</mi><mo>=</mo><mo>〈</mo><mover accent="true"><mi>ψ</mi><mo stretchy="true">¯</mo></mover><msup><mi>γ</mi><mn>0</mn></msup><mi>ψ</mi><mo>〉</mo><mo>=</mo><msub><mi>ρ</mi><mtext>n</mtext></msub><mo>+</mo><msub><mi>ρ</mi><mtext>p</mtext></msub><mn>,</mn><mi> </mi><mi> </mi><msub><mi>ρ</mi><mi>S</mi></msub><mo>=</mo><mo>〈</mo><mover accent="true"><mi>ψ</mi><mo stretchy="true">¯</mo></mover><mi>ψ</mi><mo>〉</mo><mo>=</mo><msub><mi>ρ</mi><mrow><mtext>S,n</mtext></mrow></msub><mo>+</mo><msub><mi>ρ</mi><mrow><mtext>S,p</mtext></mrow></msub><mn>,</mn></mrow></math> are the baryon density and scalar density, respectively, with the latter given by <math display="block" id="M21"><mtable columnalign="left"><mtr><mtd><msub><mi>ρ</mi><mrow><mi>S</mi><mn>,</mn><mi>J</mi></mrow></msub><mo>=</mo><mfrac><mn>2</mn><mrow><msup><mrow><mo stretchy="false">(</mo><mn>2</mn><mi>π</mi><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow></mfrac><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><msubsup><mi>k</mi><mi>F</mi><mi>J</mi></msubsup></mrow></msubsup></mrow></mstyle><mtext>d</mtext><mi>k</mi><mfrac><mrow><msup><mi>M</mi><mo>∗</mo></msup></mrow><mrow><msqrt><mrow><mo>|</mo><mi>k</mi><msup><mo>|</mo><mn>2</mn></msup><mo>+</mo><msup><mi>M</mi><mo>∗</mo></msup><msup><mrow/><mn>2</mn></msup></mrow></msqrt></mrow></mfrac></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mrow><msup><mi>M</mi><mo>∗</mo></msup></mrow><mrow><mn>2</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mrow><mo>[</mo><mrow><msubsup><mi>k</mi><mtext>F</mtext><mi>J</mi></msubsup><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mrow><mi>J</mi><mo>∗</mo></mrow></msubsup><mo>−</mo><msup><mi>M</mi><mrow><mo>∗</mo><mn>2</mn></mrow></msup><mi>ln</mi><mrow><mo>(</mo><mrow><mfrac><mrow><msubsup><mi>k</mi><mtext>F</mtext><mi>J</mi></msubsup><mo>+</mo><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mrow><mi>J</mi><mo>∗</mo></mrow></msubsup></mrow><mrow><msup><mi>M</mi><mo>∗</mo></msup></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow><mn>,</mn><mi>J</mi><mo>=</mo><mtext>p,n</mtext><mn>.</mn></mtd></mtr></mtable></math> In the above expression, we have <math id="M22"><mrow><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mrow><mi>J</mi><mo>∗</mo></mrow></msubsup><mo>=</mo><msqrt><mrow><msubsup><mi>k</mi><mtext>F</mtext><mrow><mi>J</mi><mn>2</mn></mrow></msubsup><mo>+</mo><msup><mi>M</mi><mrow><mo>*</mo><mn>2</mn></mrow></msup></mrow></msqrt></mrow></math> and the nucleon Dirac mass is defined as <math display="block" id="M23"><mrow><msup><mi>M</mi><mo>*</mo></msup><mo>≡</mo><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mo>∗</mo></msubsup><mo>=</mo><mi>M</mi><mo>−</mo><msub><mi>g</mi><mi>σ</mi></msub><mover accent="true"><mi>σ</mi><mo stretchy="true">¯</mo></mover><mn>.</mn></mrow></math> <math id="M24"><mrow><msubsup><mi>k</mi><mi>F</mi><mi>J</mi></msubsup><mo>=</mo><msub><mi>k</mi><mtext>F</mtext></msub><msup><mrow><mo stretchy="false">(</mo><mn>1</mn><mo>+</mo><msubsup><mi>τ</mi><mn>3</mn><mi>J</mi></msubsup><mi>δ</mi><mo stretchy="false">)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup></mrow></math> is the Fermi momentum with <math id="M25"><mrow><msubsup><mi>τ</mi><mn>3</mn><mtext>n</mtext></msubsup><mo>=</mo><mo>+</mo><mn>1</mn></mrow></math> for neutrons and <math id="M26"><mrow><msubsup><mi>τ</mi><mn>3</mn><mtext>p</mtext></msubsup><mo>=</mo><mo>−</mo><mn>1</mn></mrow></math> for protons, and k F =(3π 2 ρ /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 <math display="block" id="M27"><mrow><msup><mi mathvariant="script">T</mi><mrow><mi>μ</mi><mi>ν</mi></mrow></msup><mo>=</mo><mover accent="true"><mi>ψ</mi><mo stretchy="true">¯</mo></mover><mi>i</mi><msup><mi>γ</mi><mi>μ</mi></msup><msup><mo>∂</mo><mi>ν</mi></msup><mi>ψ</mi><mo>+</mo><msup><mo>∂</mo><mi>μ</mi></msup><mi>σ</mi><msup><mo>∂</mo><mi>ν</mi></msup><mi>σ</mi><mo>−</mo><msup><mi>ω</mi><mrow><mi>μ</mi><mi>η</mi></mrow></msup><msup><mo>∂</mo><mi>ν</mi></msup><msub><mi>ω</mi><mi>η</mi></msub><mo>−</mo><msup><mi>ρ</mi><mrow><mi>μ</mi><mi>η</mi></mrow></msup><msup><mo>∂</mo><mi>ν</mi></msup><msub><mi>ρ</mi><mi>η</mi></msub><mo>−</mo><mi mathvariant="script">L</mi><msup><mi>g</mi><mrow><mi>μ</mi><mi>ν</mi></mrow></msup><mn>,</mn></mrow></math> 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., <math display="block" id="M28"><mtable columnalign="left"><mtr><mtd><mi>ε</mi><mo>=</mo><mo>〈</mo><msup><mi mathvariant="script">T</mi><mrow><mn>00</mn></mrow></msup><mo>〉</mo></mtd></mtr><mtr><mtd><mo>=</mo><msubsup><mi>ε</mi><mrow><mtext>kin</mtext></mrow><mi>n</mi></msubsup><mo>+</mo><msubsup><mi>ε</mi><mrow><mtext>kin</mtext></mrow><mi>p</mi></msubsup><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>[</mo><mrow><msubsup><mi>m</mi><mi>σ</mi><mn>2</mn></msubsup><msup><mrow><mover accent="true"><mi>σ</mi><mo stretchy="true">¯</mo></mover></mrow><mn>2</mn></msup><mo>+</mo><msubsup><mi>m</mi><mi>ω</mi><mn>2</mn></msubsup><msubsup><mrow><mover accent="true"><mi>ω</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>m</mi><mi>ρ</mi><mn>2</mn></msubsup><msup><mrow><mrow><mo>(</mo><mrow><msubsup><mrow><mover accent="true"><mi>ρ</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn><mrow><mo stretchy="false">(</mo><mn>3</mn><mo stretchy="false">)</mo></mrow></msubsup></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup></mrow><mo>]</mo></mrow></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><msub><mi>b</mi><mi>σ</mi></msub><msup><mrow><mo stretchy="false">(</mo><msub><mi>g</mi><mi>σ</mi></msub><mover accent="true"><mi>σ</mi><mo>¯</mo></mover><mo stretchy="false">)</mo></mrow><mn>3</mn></msup><mo>+</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><msub><mi>c</mi><mi>σ</mi></msub><msup><mrow><mo stretchy="false">(</mo><msub><mi>g</mi><mi>σ</mi></msub><mover accent="true"><mi>σ</mi><mo>¯</mo></mover><mo stretchy="false">)</mo></mrow><mn>4</mn></msup><mo>+</mo><mfrac><mn>3</mn><mn>4</mn></mfrac><msub><mi>c</mi><mi>ω</mi></msub><msup><mrow><mo stretchy="false">(</mo><msub><mi>g</mi><mi>ω</mi></msub><msub><mover accent="true"><mi>ω</mi><mo>¯</mo></mover><mn>0</mn></msub><mo stretchy="false">)</mo></mrow><mn>4</mn></msup></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mn>3</mn><mn>2</mn></mfrac><msup><mrow><mo>(</mo><mrow><msub><mi>g</mi><mi>ρ</mi></msub><msubsup><mrow><mover accent="true"><mi>ρ</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn><mrow><mo stretchy="false">(</mo><mn>3</mn><mo stretchy="false">)</mo></mrow></msubsup></mrow><mo>)</mo></mrow><mn>2</mn></msup><msub><mi>Λ</mi><mtext>V</mtext></msub><msup><mrow><mo stretchy="false">(</mo><msub><mi>g</mi><mi>ω</mi></msub><msub><mover accent="true"><mi>ω</mi><mo>¯</mo></mover><mn>0</mn></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mn>,</mn></mtd></mtr></mtable></math> where <math display="block" id="M29"><mtable columnalign="left"><mtr><mtd><msubsup><mi>ε</mi><mrow><mtext>kin</mtext></mrow><mi>J</mi></msubsup><mo>=</mo><mfrac><mn>2</mn><mrow><msup><mrow><mo stretchy="false">(</mo><mn>2</mn><mi>π</mi><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow></mfrac><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><msubsup><mi>k</mi><mtext>F</mtext><mi>J</mi></msubsup></mrow></msubsup></mrow></mstyle><mtext>d</mtext><mi>k</mi><msqrt><mrow><mo>|</mo><mi>k</mi><msup><mo>|</mo><mn>2</mn></msup><mo>+</mo><msup><mi>M</mi><mrow><mo>∗</mo><mn>2</mn></mrow></msup></mrow></msqrt></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mn>1</mn><mrow><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><msubsup><mi>k</mi><mtext>F</mtext><mi>J</mi></msubsup></mrow></msubsup></mrow></mstyle><msup><mi>k</mi><mn>2</mn></msup><mtext>d</mtext><mi>k</mi><msqrt><mrow><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><msup><mi>M</mi><mrow><mo>∗</mo><mn>2</mn></mrow></msup></mrow></msqrt></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><mrow><mo>[</mo><mrow><mn>3</mn><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mrow><mi>J</mi><mo>∗</mo></mrow></msubsup><msub><mi>ρ</mi><mi>J</mi></msub><mo>+</mo><msup><mi>M</mi><mo>∗</mo></msup><msub><mi>ρ</mi><mrow><mtext>S</mtext><mn>,</mn><mi>J</mi></mrow></msub></mrow><mo>]</mo></mrow><mn>,</mn><mi> </mi><mi> </mi><mi>J</mi><mo>=</mo><mtext>p,n,</mtext></mtd></mtr></mtable></math> 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., <math display="block" id="M30"><mtable columnalign="left"><mtr><mtd><mi>P</mi><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover></mstyle><mo>〈</mo><msup><mi mathvariant="script">T</mi><mrow><mi>j</mi><mi>j</mi></mrow></msup><mo>〉</mo><mo>=</mo><msubsup><mi>P</mi><mrow><mtext>kin</mtext></mrow><mi>n</mi></msubsup><mo>+</mo><msubsup><mi>P</mi><mrow><mtext>kin</mtext></mrow><mi>p</mi></msubsup></mtd></mtr><mtr><mtd><mo>−</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>[</mo><mrow><msubsup><mi>m</mi><mi>σ</mi><mn>2</mn></msubsup><msup><mrow><mover accent="true"><mi>σ</mi><mo stretchy="true">¯</mo></mover></mrow><mn>2</mn></msup><mo>−</mo><msubsup><mi>m</mi><mi>ω</mi><mn>2</mn></msubsup><msubsup><mrow><mover accent="true"><mi>ω</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn><mn>2</mn></msubsup><mo>−</mo><msubsup><mi>m</mi><mi>ρ</mi><mn>2</mn></msubsup><msup><mrow><mrow><mo>(</mo><mrow><msubsup><mrow><mover accent="true"><mi>ρ</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn><mrow><mo stretchy="false">(</mo><mn>3</mn><mo stretchy="false">)</mo></mrow></msubsup></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup></mrow><mo>]</mo></mrow></mtd></mtr><mtr><mtd><mo>−</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><msub><mi>b</mi><mi>σ</mi></msub><msup><mrow><mo stretchy="false">(</mo><msub><mi>g</mi><mi>σ</mi></msub><mover accent="true"><mi>σ</mi><mo>¯</mo></mover><mo stretchy="false">)</mo></mrow><mn>3</mn></msup><mo>−</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><msub><mi>c</mi><mi>σ</mi></msub><msup><mrow><mo stretchy="false">(</mo><msub><mi>g</mi><mi>σ</mi></msub><mover accent="true"><mi>σ</mi><mo>¯</mo></mover><mo stretchy="false">)</mo></mrow><mn>4</mn></msup><mo>+</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><msub><mi>c</mi><mi>ω</mi></msub><msup><mrow><mo stretchy="false">(</mo><msub><mi>g</mi><mi>ω</mi></msub><msub><mover accent="true"><mi>ω</mi><mo>¯</mo></mover><mn>0</mn></msub><mo stretchy="false">)</mo></mrow><mn>4</mn></msup></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>(</mo><mrow><msub><mi>g</mi><mi>ρ</mi></msub><msubsup><mrow><mover accent="true"><mi>ρ</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn><mrow><mo stretchy="false">(</mo><mn>3</mn><mo stretchy="false">)</mo></mrow></msubsup></mrow><mo>)</mo></mrow><mn>2</mn></msup><msub><mi>Λ</mi><mtext>V</mtext></msub><msup><mrow><mo stretchy="false">(</mo><msub><mi>g</mi><mi>ω</mi></msub><msub><mover accent="true"><mi>ω</mi><mo>¯</mo></mover><mn>0</mn></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mn>,</mn></mtd></mtr></mtable></math> where the kinetic part of pressure is given by <math display="block" id="M31"><mrow><msubsup><mi>P</mi><mrow><mtext>kin</mtext></mrow><mi>J</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><mrow><mn>3</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><msubsup><mi>k</mi><mtext>F</mtext><mi>J</mi></msubsup></mrow></msubsup></mrow></mstyle><mtext>d</mtext><mi>k</mi><mfrac><mrow><msup><mi>k</mi><mn>4</mn></msup></mrow><mrow><msqrt><mrow><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><msup><mi>M</mi><mo>∗</mo></msup><msup><mrow/><mn>2</mn></msup></mrow></msqrt></mrow></mfrac><mn>,</mn><mi> </mi><mi> </mi><mi>J</mi><mo>=</mo><mtext>p,n</mtext><mtext>.</mtext></mrow></math> The EOS of ANM can be calculated through the energy density, ε ( ρ , δ ), by <math display="block" id="M32"><mrow><mi>E</mi><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>δ</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mi>ε</mi><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>δ</mi><mo stretchy="false">)</mo></mrow><mi>ρ</mi></mfrac><mo>−</mo><mi>M</mi><mn>.</mn></mrow></math> The EOS of SNM is just <math id="M33"><mrow><msub><mi>E</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>≡</mo><mi>E</mi><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>δ</mi><mo>=</mo><mn>0</mn><mo stretchy="false">)</mo></mrow></math> , and the characteristic parameters K 0 and J 0 can be obtained from the following expressions <math display="block" id="M34"><mrow><msub><mi>K</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>≡</mo><mn>9</mn><msup><mi>ρ</mi><mn>2</mn></msup><mfrac><mrow><msup><mtext>d</mtext><mn>2</mn></msup><msub><mi>E</mi><mn>0</mn></msub></mrow><mrow><mtext>d</mtext><msup><mi>ρ</mi><mn>2</mn></msup></mrow></mfrac><mo>=</mo><mo>−</mo><mfrac><mrow><mn>9</mn><mi>ρ</mi><msubsup><mi>g</mi><mi>σ</mi><mn>2</mn></msubsup><msup><mi>M</mi><mrow><mo>∗</mo><mn>2</mn></mrow></msup></mrow><mrow><msub><mi>Q</mi><mi>σ</mi></msub><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mrow><mo>∗</mo><mn>2</mn></mrow></msubsup></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>9</mn><mi>ρ</mi><msubsup><mi>g</mi><mi>ω</mi><mn>2</mn></msubsup></mrow><mrow><msub><mi>Q</mi><mi>ω</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>3</mn><msubsup><mi>k</mi><mtext>F</mtext><mn>2</mn></msubsup></mrow><mrow><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mo>∗</mo></msubsup></mrow></mfrac><mo>−</mo><mn>6</mn><msub><mi>L</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mn>,</mn></mrow></math> <math display="block" id="M35"><mtable columnalign="left"><mtr><mtd><msub><mi>J</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>≡</mo><mn>27</mn><msup><mi>ρ</mi><mn>3</mn></msup><mfrac><mrow><msup><mtext>d</mtext><mn>3</mn></msup><msub><mi>E</mi><mn>0</mn></msub></mrow><mrow><mtext>d</mtext><msup><mi>ρ</mi><mn>3</mn></msup></mrow></mfrac><mo>=</mo><mo>−</mo><mfrac><mrow><mn>3</mn><msubsup><mi>k</mi><mtext>F</mtext><mn>2</mn></msubsup></mrow><mrow><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mo>∗</mo></msubsup></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>3</mn><msubsup><mi>k</mi><mtext>F</mtext><mn>4</mn></msubsup></mrow><mrow><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mrow><mo>∗</mo><mn>3</mn></mrow></msubsup></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>27</mn><msubsup><mi>g</mi><mi>σ</mi><mn>2</mn></msubsup><msup><mi>M</mi><mrow><mo>∗</mo><mn>2</mn></mrow></msup><msup><mi>ρ</mi><mn>2</mn></msup></mrow><mrow><msub><mi>Q</mi><mi>σ</mi></msub><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mrow><mo>∗</mo><mn>3</mn></mrow></msubsup></mrow></mfrac><mo>×</mo></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mrow><mfrac><mrow><mn>3</mn><msup><mi>π</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><msub><mi>k</mi><mtext>F</mtext></msub><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mo>∗</mo></msubsup></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>2</mn><msubsup><mi>g</mi><mi>σ</mi><mn>2</mn></msubsup></mrow><mrow><msub><mi>Q</mi><mi>σ</mi></msub></mrow></mfrac><mo>−</mo><mfrac><mrow><msub><mi>g</mi><mi>σ</mi></msub><msup><mi>M</mi><mo>∗</mo></msup><mi>η</mi></mrow><mrow><msubsup><mi>Q</mi><mi>σ</mi><mn>2</mn></msubsup></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>2</mn><msubsup><mi>g</mi><mi>σ</mi><mn>2</mn></msubsup><msup><mi>M</mi><mrow><mo>∗</mo><mn>2</mn></mrow></msup></mrow><mrow><msub><mi>Q</mi><mi>σ</mi></msub><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mrow><mo>∗</mo><mn>2</mn></mrow></msubsup></mrow></mfrac><mo>+</mo><mfrac><mrow><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mo>∗</mo></msubsup><mi>ϕ</mi></mrow><mrow><mn>2</mn><msub><mi>Q</mi><mi>σ</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>−</mo><mfrac><mrow><mn>162</mn><msub><mi>c</mi><mi>ω</mi></msub><msubsup><mi>g</mi><mi>ω</mi><mn>7</mn></msubsup><msub><mrow><mover accent="true"><mi>ω</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn></msub><msup><mi>ρ</mi><mn>2</mn></msup></mrow><mrow><msubsup><mi>Q</mi><mi>ω</mi><mn>3</mn></msubsup></mrow></mfrac><mo>−</mo><mn>9</mn><msub><mi>K</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mn>,</mn></mtd></mtr></mtable></math> with <math display="block" id="M36"><mtable columnalign="left"><mtr><mtd><msub><mi>L</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>≡</mo><mn>3</mn><mi>ρ</mi><mfrac><mrow><mtext>d</mtext><msub><mi>E</mi><mn>0</mn></msub></mrow><mrow><mtext>d</mtext><mi>ρ</mi></mrow></mfrac><mo>=</mo><mn>3</mn><mo stretchy="false">[</mo><mfrac><mrow><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mo>∗</mo></msubsup></mrow><mn>4</mn></mfrac><mo>−</mo><mfrac><mrow><msup><mi>M</mi><mo>∗</mo></msup><msub><mi>ρ</mi><mtext>S</mtext></msub></mrow><mrow><mn>4</mn><mi>ρ</mi></mrow></mfrac><mo>+</mo><msub><mi>g</mi><mi>ω</mi></msub><msub><mover accent="true"><mi>ω</mi><mo stretchy="true">¯</mo></mover><mn>0</mn></msub></mtd></mtr><mtr><mtd><mo>−</mo><mfrac><mn>1</mn><mi>ρ</mi></mfrac><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><msubsup><mi>m</mi><mi>σ</mi><mn>2</mn></msubsup><msup><mrow><mover accent="true"><mi>σ</mi><mo stretchy="true">¯</mo></mover></mrow><mn>2</mn></msup><mo>+</mo><mi>U</mi><mo stretchy="false">(</mo><mover accent="true"><mi>σ</mi><mo stretchy="true">¯</mo></mover><mo stretchy="false">)</mo><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msubsup><mi>m</mi><mi>ω</mi><mn>2</mn></msubsup><msubsup><mrow><mover accent="true"><mi>ω</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn><mn>2</mn></msubsup><mo>+</mo><mfrac><mn>3</mn><mn>4</mn></mfrac><msub><mi>c</mi><mi>ω</mi></msub><msubsup><mi>g</mi><mi>ω</mi><mn>4</mn></msubsup><msubsup><mrow><mover accent="true"><mi>ω</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn><mn>4</mn></msubsup></mrow><mo>)</mo></mrow><mo stretchy="false">]</mo><mn>.</mn></mtd></mtr></mtable></math> In the above expressions, we have <math display="block" id="M37"><mrow><msub><mi>Q</mi><mi>σ</mi></msub><mo>=</mo><msubsup><mi>m</mi><mi>σ</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>g</mi><mi>σ</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mrow><mfrac><mrow><mn>3</mn><msub><mi>ρ</mi><mi>S</mi></msub></mrow><mrow><msup><mi>M</mi><mo>∗</mo></msup></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>3</mn><mi>ρ</mi></mrow><mrow><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mo>∗</mo></msubsup></mrow></mfrac></mrow><mo>)</mo></mrow><mo>+</mo><mn>2</mn><msub><mi>b</mi><mi>σ</mi></msub><mi>M</mi><msubsup><mi>g</mi><mi>σ</mi><mn>3</mn></msubsup><mover accent="true"><mi>σ</mi><mo stretchy="true">¯</mo></mover><mo>+</mo><mn>3</mn><msub><mi>c</mi><mi>σ</mi></msub><msubsup><mi>g</mi><mi>σ</mi><mn>4</mn></msubsup><msup><mrow><mover accent="true"><mi>σ</mi><mo stretchy="true">¯</mo></mover></mrow><mn>2</mn></msup><mn>,</mn></mrow></math> <math display="block" id="M38"><mrow><msub><mi>Q</mi><mi>ω</mi></msub><mo>=</mo><msubsup><mi>m</mi><mi>ω</mi><mn>2</mn></msubsup><mo>+</mo><mn>3</mn><msub><mi>c</mi><mi>ω</mi></msub><msubsup><mi>g</mi><mi>ω</mi><mn>4</mn></msubsup><msubsup><mrow><mover accent="true"><mi>ω</mi><mo stretchy="true">¯</mo></mover></mrow><mn>0</mn><mn>2</mn></msubsup><mn>,</mn></mrow></math> and <math display="block" id="M39"><mrow><mi>η</mi><mo>=</mo><mn>3</mn><msubsup><mi>g</mi><mi>σ</mi><mn>3</mn></msubsup><mrow><mo>(</mo><mrow><mfrac><mrow><mn>2</mn><msub><mi>ρ</mi><mi>S</mi></msub></mrow><mrow><msup><mi>M</mi><mrow><mo>∗</mo><mn>2</mn></mrow></msup></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>3</mn><mi>ρ</mi></mrow><mrow><msup><mi>M</mi><mo>∗</mo></msup><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mo>∗</mo></msubsup></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mi>M</mi><mo>∗</mo></msup><mi>ρ</mi></mrow><mrow><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mrow><mo>∗</mo><mn>3</mn></mrow></msubsup></mrow></mfrac></mrow><mo>)</mo></mrow><mo>−</mo><mn>2</mn><msub><mi>b</mi><mi>σ</mi></msub><mi>M</mi><msubsup><mi>g</mi><mi>σ</mi><mn>3</mn></msubsup><mo>−</mo><mn>6</mn><msub><mi>c</mi><mi>σ</mi></msub><msubsup><mi>g</mi><mi>σ</mi><mn>4</mn></msubsup><mover accent="true"><mi>σ</mi><mo stretchy="true">¯</mo></mover><mn>,</mn></mrow></math> <math display="block" id="M40"><mrow><mi>ϕ</mi><mo>=</mo><mfrac><mrow><mn>2</mn><msubsup><mi>g</mi><mi>σ</mi><mn>2</mn></msubsup><msubsup><mi>k</mi><mtext>F</mtext><mn>2</mn></msubsup></mrow><mrow><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mrow><mo>∗</mo><mn>3</mn></mrow></msubsup></mrow></mfrac><mn>,</mn></mrow></math> with <math id="M41"><mrow><msubsup><mi mathvariant="script">E</mi><mtext>F</mtext><mo>∗</mo></msubsup><mo>=</mo><mo stretchy="false">(</mo><msubsup><mi>k</mi><mtext>F</mtext><mn>2</mn></msubsup><mo>+</mo><msup><mi>M</mi><mrow><mo>∗</mo><mn>2</mn></mrow></msup><msup><mo stretchy="false">)</mo><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup></mrow></math> . To our best knowledge, Eq. (23) gives, for the first time [ 52 52 ], the analytical expression of the J 0 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 53 ]. 3 Results and discussions 1 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 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., ρ 0 , E 0 ( ρ 0 ), K 0 , J 0 , <math id="M42"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup><mo>≡</mo><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mo>∗</mo></msubsup><mo stretchy="false">(</mo><msub><mi>ρ</mi><mn>0</mn></msub><mo stretchy="false">)</mo></mrow></math> , 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 33 ]). Then, by varying individually these macroscopic quantities within their known ranges, one can examine transparently the correlatio n 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 54 , 33 33 ]) and giant resonances of finite nuclei ([ 30 30 , 55 55 ]), the higher-order bulk characteristic parameters of ANM ([ 25 25 ]), and the relationship between the nuclear matter symmetry energy and the symmetry energy coefficient in the mass formula ([ 56 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 Fig. 1 the pressure of SNM P ( ρ ) at ρ =3 ρ 0 from the nonlinear RMF model based on the FSUGold interaction ([ 49 49 ]) by varying individually ρ 0 , E 0 ( ρ 0 ), <math id="M43"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> , K 0 , and J 0 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 =0.148 fm -3 , E 0 ( ρ 0 )=-16.3 MeV, <math id="M44"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup><mo>=</mo><mn>0.61</mn><mi>M</mi></mrow></math> , K 0 =230 MeV, J 0 =-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 Fig. 1 that the pressure of SNM P ( ρ ) at ρ =3 ρ 0 increases with ρ 0 , <math id="M45"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> , K 0 , and J 0 while it decreases with E 0 ( ρ 0 ). In particular, the pressure of SNM P ( ρ ) at ρ =3 ρ 0 displays a specially strong correlation with J 0 . We note that the pressure of SNM P ( ρ ) at other supra-saturation densities exhibits similar correlations with ρ 0 , E 0 ( ρ 0 ), <math id="M46"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> , K 0 , and J 0 . These features indicate that the pressure of SNM at supra-saturation densities is sensitive to the J 0 value, and thus the experimental constraints on the pressure of SNM at supra-saturation densities may provide important information on the J 0 value. Since the pressure of SNM at supra-saturation densities is sensitive to the J 0 value, the maximum mass, M max , of static neutron stars is also expected to be sensitive to the J 0 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 <math id="M48"><mrow><mi>P</mi><mo>=</mo><mi>a</mi><mo>+</mo><mi>b</mi><msup><mi mathvariant="script">E</mi><mrow><mn>4</mn><mo>/</mo><mn>3</mn></mrow></msup></mrow></math> because of our poor knowledge about its EOS from first principle ([ 57 57 , 58 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 58 ]. In this work, the ρ t is determined self-consistently within the nonlinear RMF model using the thermodynamical method (see, e.g., [ 51 51 ] for the details). In the outer crust with 6.93× 10 -13 fm -3 < ρ < ρ out , we use the EOS of BPS ([ 59 59 , 60 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 59 ]). Similarly, as in Fig. 1 Fig. 1 , we plot in Fig. 2 Fig. 2 the maximum mass, M max , of static neutron stars from the nonlinear RMF model based on the FSUGold interaction by varying individually ρ 0 , E 0 ( ρ 0 ), <math id="M49"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> , K 0 , J 0 , 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 0 parameter. In addition, the M max exhibits weak positive correlation with the <math id="M50"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> and K 0 , and weak negative correlation with the ρ 0 and E 0 ( ρ 0 ). 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.74 M ⊙ . These features indicate that the observed largest mass of neutron stars may put important constraint on the J 0 value. We would like to point out that the interaction FSUGold is only used in Figs. 1 Figs. 1 and 2 2 as a reference for the correlation analyses and using other RMF interactions will not change our conclusion. Experimentally, the pressure of SNM at supra-saturation densities (from 2 ρ 0 to about 5 ρ 0 ) has been constrained by measurements of collective flows in HIC ([ 3 3 ]), which is shown as a band in the left window of Fig. 3 Fig. 3 . In the nonlinear RMF model, if one only changes the J 0 value while the other 6 macroscopic quantities are kept at their values in the FSUGold interaction, one can find that the J 0 value should be in the range of -985 MeV ≤ J 0 ≤ -327 MeV to be consistent with the flow data in HIC ([ 3 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 0 value from the flow data in HIC will also depends on the values of ρ 0 , E 0 ( ρ 0 ), <math id="M52"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> , and K 0 , 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 =0.153±0.008 fm -3 , E 0 ( ρ 0 )=-16.2±0.3 MeV, <math id="M53"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>*</mo><mn>0</mn></mrow></msubsup><mo>/</mo><mi>M</mi><mo>=</mo><mn>0.61</mn><mo>±</mo><mn>0.04</mn></mrow></math> , and K 0 =230±20 MeV, which represent the typical uncertain ranges known or predicted from different interactions in the nonlinear RMF model ([ 50 50 ]). Based on the pressure of SNM constrained by flow data in HIC ([ 3 3 ]), to extract the upper limit of the J 0 value, one should use the values of ρ 0 , E 0 ( ρ 0 ), <math id="M56"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> , and K 0 that make the resulting pressure of SNM as small as possible when J 0 is fixed. This can be obtained by using ρ 0 =0.145 fm -3 , E 0 ( ρ 0 )=-15.9 MeV, <math id="M57"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mo>*</mo></msubsup><mo>/</mo><mi>M</mi><mo>=</mo><mn>0.57</mn></mrow></math> , and K 0 =210 MeV, denoted as set “S", since the pressure of SNM P ( ρ / ρ 0 ) at supra-saturation densities increases with ρ 0 , <math id="M58"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> , and K 0 while decreases with E 0 ( ρ 0 ), as shown in Fig. 1 Fig. 1 . With the set “S" for ρ 0 , E 0 ( ρ 0 ), <math id="M59"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> , and K 0 , one can find the upper limit of J 0 = -10 MeV for the J 0 value, which is indicated by a solid line in the left window of Fig. 3 Fig. 3 . For J 0 > -10 MeV, the model would over-predict the pressure of SNM constrained by flow data in HIC ([ 3 3 ]). Similarly, one can obtain the lower limit of the J 0 value by using the values of ρ 0 , E 0 ( ρ 0 ), <math id="M60"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> , and K 0 that make the resulting pressure of SNM as large as possible when J 0 is fixed, and this can be obtained with ρ 0 =0.161 fm -3 , E 0 ( ρ 0 )=-16.5 MeV, <math id="M61"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup><mo>/</mo><mi>M</mi><mo>=</mo><mn>0.65</mn></mrow></math> , and K 0 =250 MeV, denoted as set “H". Using the set “H" for ρ 0 , E 0 ( ρ 0 ), <math id="M62"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> , and K 0 , one can extract the lower limit of J 0 = -1280 MeV, which is indicated by dashed line in the left window of Fig. 3 Fig. 3 . The model would under-predict the pressure of SNM constrained by flow data in HIC ([ 3 3 ]) if J 0 < -1280 MeV. Therefore, from the pressure of SNM constrained by flow data in HIC ([ 3 3 ]), one can extract the constraint of -1280 MeV ≤ J 0 ≤ -10 MeV. Recently, a new neutron star, PSR J0348+0432, with a mass of 2.01±0.04 M ⊙ was discovered ([ 61 61 ]), and this neutron star is only the second pulsar with a precisely determined mass around 2 M ⊙ after PSR J1614-2230 ([ 62 62 ]) and sets a new record of the maximum mass of neutron stars. The lower mass limit of 1.97 M ⊙ for PSR J0348+0432 thus may set a lower limit of the J 0 value below which the model cannot predict a neutron star with mass equal or above 1.97 M ⊙ . To extract the lower limit of the J 0 value from the observed heaviest neutron star, PSR J0348+0432, one can use the values of ρ 0 , E 0 ( ρ 0 ), <math id="M63"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> , K 0 , E sym ( ρ c ), and L ( ρ c ) that make the resulting maximum mass of neutron stars as large as possible when J 0 is fixed, and from Fig. 2 Fig. 2 this can be obtained with ρ 0 =0.145 fm -3 , E 0 ( ρ 0 )=-16.5 MeV, K 0 =250 MeV, <math id="M64"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup><mo>/</mo><mi>M</mi><mo>=</mo><mn>0.65</mn></mrow></math> , E sym ( ρ c )=26 MeV, and L ( ρ c )=60 MeV, denoted as set “NS-H". This leads to a lower limit of J 0 = -494 MeV for the J 0 value as shown in the right window of Fig. 3 Fig. 3 where the maximum mass of neutron stars is plotted as a function of J 0 when the other 6 macroscopic quantities are fixed at their values as in set “NS-H". For J 0 < -494 MeV, the maximum mass of static neutron stars predicted in the nonlinear RMF model would be always smaller than 1.97 M ⊙ . 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 0 value would be necessary to obtain a neutron star with mass of 1.97 M ⊙ . Therefore, including the new degrees of freedom in neutron stars will be consistent with the constraint of J 0 ≥ -494 MeV. Combining the constraint of -1280 MeV ≤ J 0 ≤ -10 MeV from the pressure of SNM constrained by flow data in HIC ([ 3 3 ]) which favors a smaller J 0 value and the constraint of J 0 ≥ -494 MeV from the recently discovered heaviest neutron star PSR J0348+0432 ([ 61 61 ]) which favors a larger J 0 value, one can extract the following constraint for the J 0 parameter <math display="block" id="M65"><mrow><mo>−</mo><mn>494</mn><mtext> </mtext><mtext>MeV</mtext><mo>≤</mo><msub><mi>J</mi><mn>0</mn></msub><mo>≤</mo><mo>−</mo><mn>10</mn><mtext> </mtext><mtext>MeV</mtext><mn>.</mn></mrow></math> It should be emphasized that the constraint -494 MeV≤ J 0 ≤-10 MeV represents a conservative extraction based on flow data in HIC ([ 3 3 ]) and the recently discovered heaviest neutron star, PSR J0348+0432 ([ 61 61 ]) in the nonlinear RMF model. This is because we have not considered the possible correlations that existed among the ρ 0 , E 0 ( ρ 0 ), <math id="M66"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> , K 0 , 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 ρ 0 , E 0 ( ρ 0 ), <math id="M67"><mrow><msubsup><mi>M</mi><mrow><mtext>dirac</mtext></mrow><mrow><mo>∗</mo><mn>0</mn></mrow></msubsup></mrow></math> , K 0 , E sym ( ρ c ) and L ( ρ c ) should further narrow the constraint -494 MeV≤ J 0 ≤-10 MeV, and it will be interesting to see the quantitative constraint on the J 0 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 0 ≤-10 MeV obtained in the present work indicates that, if the J 0 value is out of the region -494 MeV≤ J 0 ≤-10 MeV, the nonlinear RMF model either cannot predict the pressure of SNM constrained by flow data in HIC ([ 3 3 ]) or cannot describe the recently discovered heaviest neutron star PSR J0348+0432 ([ 61 61 ]). It is worth mentioning that the constraint on the J 0 depends on our knowledge of the other six quantities. Any improvement on these six macroscopic quantities will make the range for the J 0 constraint narrower. In addition, the extracted constraint on the J 0 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 Fig. 4 , we show the comparison of the J 0 constraint obtained in our analysis with those obtained with other analyses and/or other methods ([ 22 22 , 25 25 , 63 63 , 64 64 ]), including the constraint of J 0 = -700 ± 500 MeV obtained by Ref. [ 63 63 ] from the analysis of nuclear GMR, the constraint of <math id="M68"><mrow><msub><mi>J</mi><mn>0</mn></msub><mo>=</mo><mo>−</mo><msubsup><mrow><mn>280</mn></mrow><mrow><mo>−</mo><mn>410</mn></mrow><mrow><mo>+</mo><mn>72</mn></mrow></msubsup></mrow></math> ( <math id="M69"><mrow><mo>−</mo><msubsup><mrow><mn>500</mn></mrow><mrow><mo>−</mo><mn>290</mn></mrow><mrow><mo>+</mo><mn>170</mn></mrow></msubsup></mrow></math> ) MeV obtained by Ref. [ 64 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 0 = -355 ± 95 MeV deduced by Ref. [ 25 25 ] based on a correlation analysis method within SHF energy density functional, and the constraint of J 0 = -390 ± 90 MeV deduced by Ref. [ 22 22 ] who used a similar method as Ref. [ 25 25 ]. It is seen that the constrained region of J 0 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 0 parameter should be larger than about -500 MeV. 4 Summary 1 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 0 , 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 0 while the latter requires a larger J 0 , we have extracted a constraint of -494 MeV≤ J 0 ≤-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 0 parameter cannot be too small, namely, it should be larger than about -500 MeV. The present constraint on the J 0 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 J.P. Blaizot , Nuclear compressibilities . Phys. Rep . 64 , 171 - 248 ( 1980 ). doi: 10.1016/0370-1573(80)90001-0 http://doi.org/10.1016/0370-1573(80)90001-0 Nuclear compressibilities B.A. Li , C.M. Ko , W. Bauer , Isospin physics in heavy-Ion collisions at intermediate energies . Int. J. Mod. Phys. E 7 , 147 ( 1998 ). doi: 10.1142/S0218301398000087 http://doi.org/10.1142/S0218301398000087 Isospin physics in heavy-Ion collisions at intermediate energies P. Danielewicz , R. Lacey , W.G. Lynch , Determination of the equation of state of dense matter science , 298 , 1592 - 1596 ( 2002 ). doi: 10.1126/science.1078070 http://doi.org/10.1126/science.1078070 Determination of the equation of state of dense matter V. Baran , M. Colonna , V. Greco , et al ., Reaction dynamics with exotic nuclei , Phys. Rep . 410 , 335 - 466 ( 2005 ). doi: 10.1016/j.physrep.2004.12.004 http://doi.org/10.1016/j.physrep.2004.12.004 Reaction dynamics with exotic nuclei A.W. Steiner , M. Prakash , J.M. Lattimer , et al ., Isospin asymmetry in nuclei and neutron stars . Phys. Rep . 411 , 325 - 375 ( 2005 ). doi: 10.1016/j.physrep.2005.02.004 http://doi.org/10.1016/j.physrep.2005.02.004 Isospin asymmetry in nuclei and neutron stars L.W. Chen , C.M. Ko , B.A. Li , et al ., Probing the nucle- ar symmetry energy with heavy-ion reactions induced by neutron-rich nuclei . Front. Phys. China 2 , 327 - 357 ( 2007 ). doi: 10.1007/s11467-007-0037-0 http://doi.org/10.1007/s11467-007-0037-0 Probing the nucle- ar symmetry energy with heavy-ion reactions induced by neutron-rich nuclei B.A. Li , L.W. Chen , C.M. Ko , Recent progress and new challenges in isospin physics with heavy-ion reactions . Phys. Rep . 464 , 113 - 281 ( 2008 ). doi: 10.1016/j.physrep.2008.04.005 http://doi.org/10.1016/j.physrep.2008.04.005 Recent progress and new challenges in isospin physics with heavy-ion reactions J. B. Natowitz , G. Röpke , S. Typel , et al ., Symmetry energy of dilute warm nuclear matter . Phys. Rev. Lett . 104 , 202501 ( 2010 ). doi: 10.1103/PhysRevLett.104.202501 http://doi.org/10.1103/PhysRevLett.104.202501 Symmetry energy of dilute warm nuclear matter B.M. Tsang , J. R. Stone , F. Camera , et al ., Constraints on the symmetry energy and neutron skins from experiments and theory . Phys. Rev. C 86 , 015803 ( 2012 ). doi: 10.1103/PhysRevC.86.015803 http://doi.org/10.1103/PhysRevC.86.015803 Constraints on the symmetry energy and neutron skins from experiments and theory W. Trautmann , H.H. Wolter , Elliptic flow and the symmetry energy at supra-saturation density . Int. J. Mod. Phys. E 21 , 1230003 ( 2012 ). doi: 10.1142/S0218301312300032 http://doi.org/10.1142/S0218301312300032 Elliptic flow and the symmetry energy at supra-saturation density L.W. Chen , C.M. Ko , B.A. Li , et al ., Probing isospin- and momentum-dependent nuclear effective interactions in neutron-rich matter . Eur. Phys. J. A 50 , 29 ( 2014 ). doi: 10.1140/epja/i2014-14029-6 http://doi.org/10.1140/epja/i2014-14029-6 Probing isospin- and momentum-dependent nuclear effective interactions in neutron-rich matter C.J. Horowitz , E.F. Brown , Y. Kim , et al ., A way forward in the study of the symmetry energy: experiment, theory, and observation . J. Phys. G 41 , 093001 ( 2014 ). doi: 10.1088/0954-3899/41/9/093001 http://doi.org/10.1088/0954-3899/41/9/093001 A way forward in the study of the symmetry energy: experiment, theory, and observation B.A. Li , A. Ramos , G. Verde , et al ., Topical issue on nuclear symmetry energy . Eur. Phys. J. A 50 , 9 ( 2014 ). doi: 10.1140/epja/i2014-14009-x http://doi.org/10.1140/epja/i2014-14009-x Topical issue on nuclear symmetry energy X.Q. Liu , M.R. Huang , R. Wada , et al ., Symmetry energy extraction from primary fragments in intermediate heavy-ion collisions . Nucl. Sci. Tech . 26 , S20508 ( 2015 ) doi: 10.13538/j.1001-8042/nst.26.S20508 http://doi.org/10.13538/j.1001-8042/nst.26.S20508 ; Symmetry energy extraction from primary fragments in intermediate heavy-ion collisions F. F. Duan , X.Q. Liu , W.P. Lin , et al ., Investigation on symmetry and characteristic properties of the fragmenting source in heavy-ion reactions through reconstructed primary isotope yields . Nucl. Sci. Tech. 27 , 131 ( 2016 ). doi: 10.1007/s41365-016-0138-y http://doi.org/10.1007/s41365-016-0138-y Investigation on symmetry and characteristic properties of the fragmenting source in heavy-ion reactions through reconstructed primary isotope yields M. Baldo , G.F. Burgio , The nuclear symmetry energy . Prog. Part. Nucl. Phys . 91 , 203 - 258 ( 2016 ). doi: 10.1016/j.ppnp.2016.06.006 http://doi.org/10.1016/j.ppnp.2016.06.006 The nuclear symmetry energy B.A. Li , Nucl. Phys. News, in press , ( 2017 ) [ arXiv:1701.03564 ] Nucl. Phys. News, in press N.K. Glendenning , Compact Stars , 2 nd edition, Spinger-Verlag New York, Inc ., 2000 . Compact Stars J.M. Lattimer and M. Prakash , The physics of neutron stars . Science 304 , 536 - 542 ( 2004 ). doi: 10.1126/science.1090720 http://doi.org/10.1126/science.1090720 ; The physics of neutron stars J.M. Lattimer and M. Prakash , Neutron star observations: Prognosis for equation of state constraints . Phys. Rep. 442 , 109 - 165 ( 2007 ). doi: 10.1016/j.physrep.2007.02.003 http://doi.org/10.1016/j.physrep.2007.02.003 Neutron star observations: Prognosis for equation of state constraints J.M. Lattimer , The nuclear equation of state and neutron star masses . Annu. Rev. Nucl. Part. Sci . 62 , 485 - 515 ( 2012 ). doi: 10.1146/annurev-nucl-102711-095018 http://doi.org/10.1146/annurev-nucl-102711-095018 The nuclear equation of state and neutron star masses K. Kotake , K. Sato , K. Takahashi , Explosion mechanism, neutrino burst and gravitational wave in core-collapse supernovae . Rep. Prog. Phys. 69 , 971 - 1143 ( 2006 ) doi: 10.1088/0034-4885/69/4/R03 http://doi.org/10.1088/0034-4885/69/4/R03 ; Explosion mechanism, neutrino burst and gravitational wave in core-collapse supernovae H.-Th. Janka , K. Langanke , A. Marek , et al ., Theory of core-collapse supernovae . Phys. Rep . 442 , 38 - 74 ( 2007 ). doi: 10.1016/j.physrep.2007.02.002 http://doi.org/10.1016/j.physrep.2007.02.002 Theory of core-collapse supernovae M. Hempel , T. Fischer , J. Schaffner-Bielich , et al ., New equations of state in simulations of core-collapse supernovae . Astrophys. J . 748 , 70 ( 2012 ). doi: 10.1088/0004-637X/748/1/70 http://doi.org/10.1088/0004-637X/748/1/70 New equations of state in simulations of core-collapse supernovae M. Meixner , J.P. Olson , G. Mathews , et al ., The NDL equation of state for supernova simulations . arXiv:1303.0064 , ( 2013 ). The NDL equation of state for supernova simulations M. Oertel , M. Hempel , T. Klähn , et al ., Equations of state for supernovae and compact stars . Rev. Mod. Phys . 89 , 015007 ( 2017 ). doi: 10.1103/RevModPhys.89.015007 http://doi.org/10.1103/RevModPhys.89.015007 Equations of state for supernovae and compact stars L.W. Chen , B.J. Cai , C.M. Ko , et al ., Higher-order effects on the incompressibility of isospin asymmetric nuclear matter . Phys. Rev. C 80 , 014322 ( 2009 ). doi: 10.1103/PhysRevC.80.014322 http://doi.org/10.1103/PhysRevC.80.014322 Higher-order effects on the incompressibility of isospin asymmetric nuclear matter L.W. Chen , Higher order bulk characteristic parameters of asymmetric nuclear matter . Sci. China: Phys. Mech. Astron . 54 , s124 - s129 ( 2011 ). doi: 10.1007/s11433-011-4415-9 http://doi.org/10.1007/s11433-011-4415-9 Higher order bulk characteristic parameters of asymmetric nuclear matter D.H. Youngblood , H.L. Clark , Y.-W. Lui , Incompressibility of nuclear matter from the giant monopole resonance . Phys. Rev. Lett . 82 , 691 - 694 ( 1999 ). doi: 10.1103/PhysRevLett.82.691 http://doi.org/10.1103/PhysRevLett.82.691 Incompressibility of nuclear matter from the giant monopole resonance S. Shlomo , V.M. Kolomietz , and G Colò , Deducing the nuclear-matter incompressibility coefficient from data on isoscalar compression modes . Eur. Phys. J. A 30 , 23 - 30 ( 2006 ). doi: 10.1140/epja/i2006-10100-3 http://doi.org/10.1140/epja/i2006-10100-3 Deducing the nuclear-matter incompressibility coefficient from data on isoscalar compression modes G. Colò , Constraints, limits and extensions for nuclear energy functionals . AIP Conf. Proc . 1128 , 59 ( 2009 ). doi: 10.1063/1.3146221 http://doi.org/10.1063/1.3146221 Constraints, limits and extensions for nuclear energy functionals J. Piekarewicz , Do we understand the incompressibility of neutron-rich matter? J. Phys. G 37 , 064038 ( 2010 ). doi: 10.1088/0954-3899/37/6/064038 http://doi.org/10.1088/0954-3899/37/6/064038 Do we understand the incompressibility of neutron-rich matter? L.W. Chen , J.Z. Gu , Correlations between the nuclear breathing mode energy and properties of asymmetric nuclear matter . J. Phys. G 39 , 035104 ( 2012 ). doi: 10.1088/0954-3899/39/3/035104 http://doi.org/10.1088/0954-3899/39/3/035104 Correlations between the nuclear breathing mode energy and properties of asymmetric nuclear matter L.W. Chen , Recent progress on the determination of the symmetry energy . Nuclear Structure in China 2012 , 43 - 54 ( 2013 ). doi: 10.1142/9789814447485_0007 http://doi.org/10.1142/9789814447485_0007 arXiv:1212.0284. Recent progress on the determination of the symmetry energy B.A. Li , L.W. Chen , F.J. Fattoyev , et al ., Probing nuclear symmetry energy and its imprints on properties of nuclei, nuclear reactions, neutron stars and gravitational waves . J. Phys.: Conf. Ser . 413 , 012021 ( 2013 ) doi: 10.1088/1742-6596/413/1/012021 http://doi.org/10.1088/1742-6596/413/1/012021 Probing nuclear symmetry energy and its imprints on properties of nuclei, nuclear reactions, neutron stars and gravitational waves Z. Zhang , L.W. Chen , Constraining the symmetry energy at subsaturation densities using isotope binding energy difference and neutron skin thickness . Phys. Lett. B 726 , 234 - 238 ( 2013 ). doi: 10.1016/j.physletb.2013.08.002 http://doi.org/10.1016/j.physletb.2013.08.002 Constraining the symmetry energy at subsaturation densities using isotope binding energy difference and neutron skin thickness B.D. Serot and J.D. Walecka , Advances in Nuclear Physics . Vol. 16 , J.W. Negele , E. Vogt , Eds., Plenum , New York ( 1986 ); Advances in Nuclear Physics B.D. Serot , J.D. Walecka , Recent progress in quantum hadrodynamics . Int. J. Mod. Phys. E 6 , 515 ( 1997 ). doi: 10.1142/S0218301397000299 http://doi.org/10.1142/S0218301397000299 Recent progress in quantum hadrodynamics P.-G. Reinhard , The relativistic mean-field description of nuclei and nuclear dynamics . Rep. Prog. Phys . 52 , 439 - 514 ( 1989 ). doi: 10.1088/0034-4885/52/4/002 http://doi.org/10.1088/0034-4885/52/4/002 The relativistic mean-field description of nuclei and nuclear dynamics P. Ring , Relativistic mean field theory in finite nuclei . Prog. Part. Nucl. Phys . 37 , 193 - 263 ( 1996 ). doi: 10.1016/0146-6410(96)00054-3 http://doi.org/10.1016/0146-6410(96)00054-3 Relativistic mean field theory in finite nuclei J. Meng , H. Toki , S.G. Zhou , et al ., Relativistic continuum Hartree Bogoliubov theory for ground-state properties of exotic nuclei . Prog. Part. Nucl. Phys . 57 , 470 - 563 ( 2006 ). doi: 10.1016/j.ppnp.2005.06.001 http://doi.org/10.1016/j.ppnp.2005.06.001 Relativistic continuum Hartree Bogoliubov theory for ground-state properties of exotic nuclei Y. Sugahara , H. Toki , Relativistic mean-field theory for unstable nuclei with non-linear σ and ω terms . Nucl. Phys. A 579 , 557 - 572 ( 1994 ). doi: 10.1016/0375-9474(94)90923-7 http://doi.org/10.1016/0375-9474(94)90923-7 Relativistic mean-field theory for unstable nuclei with non-linear σ and ω terms Z.Z. Ren , Z.Y. Zhu , Y.H. Cai , et al ., Relativistic mean-field study of Mg isotopes . Phys. Lett. B 380 , 241 - 246 ( 1996 ). doi: 10.1016/0370-2693(96)00462-5 http://doi.org/10.1016/0370-2693(96)00462-5 Relativistic mean-field study of Mg isotopes G.A. Lalazissis , J. König , P. Ring , New parametrization for the Lagrangian density of relativistic mean field theory . Phys. Rev. C 55 , 540 ( 1997 ). doi: 10.1103/PhysRevC.55.540 http://doi.org/10.1103/PhysRevC.55.540 New parametrization for the Lagrangian density of relativistic mean field theory W.H. Long , J. Meng , N. Van Giai , et al ., New effective interactions in relativistic mean field theory with nonlinear terms and density-dependent meson-nucleon coupling . Phys. Rev. C 69 , 034319 ( 2004 ). doi: 10.1103/PhysRevC.69.034319 http://doi.org/10.1103/PhysRevC.69.034319 New effective interactions in relativistic mean field theory with nonlinear terms and density-dependent meson-nucleon coupling W.Z. Jiang , Z.Z. Ren , T.T. Wang , et al ., Relativistic mean-field study for Zn isotopes . Eur. Phys. J. A 25 , 29 - 39 ( 2005 ). doi: 10.1140/epja/i2004-10235-1 http://doi.org/10.1140/epja/i2004-10235-1 Relativistic mean-field study for Zn isotopes W.Z. Jiang , Effects of the density dependence of the nuclear symmetry energy on the properties of superheavy nuclei . Phys. Rev. C 81 , 044306 ( 2010 ). doi: 10.1103/PhysRevC.81.044306 http://doi.org/10.1103/PhysRevC.81.044306 Effects of the density dependence of the nuclear symmetry energy on the properties of superheavy nuclei F.J. Fattoyev , C.J. Horowitz , J. Piekarewicz , et al ., Relativistic effective interaction for nuclei, giant resonances, and neutron stars . Phys. Rev. C 82 , 055803 ( 2010 ). doi: 10.1103/PhysRevC.82.055803 http://doi.org/10.1103/PhysRevC.82.055803 Relativistic effective interaction for nuclei, giant resonances, and neutron stars B.K. Agrawal , A. Sulaksono , P.-G. Reinhard , Optimization of relativistic mean field model for finite nuclei to neutron star matter . Nucl. Phys. A 882 , 1 - 20 ( 2012 ). doi: 10.1016/j.nuclphysa.2012.03.004 http://doi.org/10.1016/j.nuclphysa.2012.03.004 Optimization of relativistic mean field model for finite nuclei to neutron star matter F.J. Fattoyev , J. Carvajal , W.G. Newton , et al ., Constraining the high-density behavior of the nuclear symmetry energy with the tidal polarizability of neutron stars . Phys. Rev. C 87 , 015806 ( 2013 ). doi: 10.1103/PhysRevC.87.015806 http://doi.org/10.1103/PhysRevC.87.015806 Constraining the high-density behavior of the nuclear symmetry energy with the tidal polarizability of neutron stars H. Müller , B. D. Serot , Relativistic mean-field theory and the high-density nuclear equation of state . Nucl. Phys. A 606 , 508 - 537 ( 1996 ). doi: 10.1016/0375-9474(96)00187-X http://doi.org/10.1016/0375-9474(96)00187-X Relativistic mean-field theory and the high-density nuclear equation of state C.J. Horowitz , J. Piekarewicz , Neutron star structure and the neutron radius of 208Pb . Phys. Rev. Lett 86 , 5647 ( 2001 ). doi: 10.1103/PhysRevLett.86.5647 http://doi.org/10.1103/PhysRevLett.86.5647 ; Neutron star structure and the neutron radius of 208Pb C.J. Horowitz , J. Piekarewicz , Neutron radii of 208Pb and neutron stars . Phys. Rev. C 64 , 062802 (R) ( 2001 ). doi: 10.1103/PhysRevC.64.062802 http://doi.org/10.1103/PhysRevC.64.062802 ; Neutron radii of 208Pb and neutron stars C. J. Horowitz , J. Piekarewicz , Constraining URCA cooling of neutron stars from the neutron radius of 208Pb . Phys. Rev. C 66 , 055803 ( 2002 ). doi: 10.1103/PhysRevC.66.055803 http://doi.org/10.1103/PhysRevC.66.055803 Constraining URCA cooling of neutron stars from the neutron radius of 208Pb B.G. Todd-Rutel , J. Piekarewicz , Neutron-rich nuclei and neutron stars: A new accurately calibrated interaction for the study of neutron-rich matter . Phys. Rev. Lett . 95 , 122501 ( 2005 ). doi: 10.1103/PhysRevLett.95.122501 http://doi.org/10.1103/PhysRevLett.95.122501 Neutron-rich nuclei and neutron stars: A new accurately calibrated interaction for the study of neutron-rich matter L.W. Chen , C.M. Ko , B.A. Li , Isospin-dependent properties of asymmetric nuclear matter in relativistic mean field models . Phys. Rev. C 76 , 054316 ( 2007 ). doi: 10.1103/PhysRevC.76.054316 http://doi.org/10.1103/PhysRevC.76.054316 Isospin-dependent properties of asymmetric nuclear matter in relativistic mean field models B.J. Cai , L.W. Chen , Nuclear matter fourth-order symmetry energy in the relativistic mean field models . Phys. Rev. C 85 , 024302 ( 2012 ). doi: 10.1103/PhysRevC.85.024302 http://doi.org/10.1103/PhysRevC.85.024302 Nuclear matter fourth-order symmetry energy in the relativistic mean field models Equation (19) was given in the first version of the present paper, i.e. , arXiv:1402.4242v1 [nucl-th] , in February, 2014 . Equation (19) was given in the first version of the present paper, i.e. B.J. Cai , L.W. Chen , Lorentz covariant nucleon self-energy decomposition of the nuclear symmetry energy . Phys. Lett. B 711 , 104 - 108 ( 2012 ). doi: 10.1016/j.physletb.2012.03.058 http://doi.org/10.1016/j.physletb.2012.03.058 Lorentz covariant nucleon self-energy decomposition of the nuclear symmetry energy L.W. Chen , C.M. Ko , B.A. Li , et al ., Density slope of the nuclear symmetry energy from the neutron skin thickness of heavy nuclei . Phys. Rev. C 82 , 024321 ( 2010 ). doi: 10.1103/PhysRevC.82.024321 http://doi.org/10.1103/PhysRevC.82.024321 Density slope of the nuclear symmetry energy from the neutron skin thickness of heavy nuclei Z. Zhang , L.W. Chen , Constraining the density slope of nuclear symmetry energy at subsaturation densities using electric dipole polarizability in 208Pb . Phys. Rev. C 90 , 064317 ( 2014 ). doi: 10.1103/PhysRevC.90.064317 http://doi.org/10.1103/PhysRevC.90.064317 Constraining the density slope of nuclear symmetry energy at subsaturation densities using electric dipole polarizability in 208Pb L.W. Chen , Nuclear matter symmetry energy and the symmetry energy coefficient in the mass formula . Phys. Rev. C 83 , 044308 ( 2011 ). doi: 10.1103/PhysRevC.83.044308 http://doi.org/10.1103/PhysRevC.83.044308 Nuclear matter symmetry energy and the symmetry energy coefficient in the mass formula J. Carriere , C.J. Horowitz , J. Piekarewicz , Low-mass neutron stars and the equation of state of dense matter . Astrophys. J . 593 , 463 - 471 ( 2003 ). doi: 10.1086/376515 http://doi.org/10.1086/376515 Low-mass neutron stars and the equation of state of dense matter J. Xu , L.W. Chen , B.A. Li , et al ., Locating the inner edge of the neutron star crust using terrestrial nuclear laboratory data . Phys. Rev. C 79 , 035802 ( 2009 ). doi: 10.1103/PhysRevC.79.035802 http://doi.org/10.1103/PhysRevC.79.035802 ; Locating the inner edge of the neutron star crust using terrestrial nuclear laboratory data J. Xu , L.W. Chen , B.A. Li , et al . Nuclear constraints on properties of neutron star crusts . Astrophys. J. 697 , 1549 - 1568 ( 2009 ). doi: 10.1088/0004-637X/697/2/1549 http://doi.org/10.1088/0004-637X/697/2/1549 Nuclear constraints on properties of neutron star crusts G. Baym , C. Pethick , P. Sutherland , The ground state of matter at high densities: Equation of state and stellar models . Astrophys. J . 170 , 299 ( 1971 ). doi: 10.1086/151216 http://doi.org/10.1086/151216 The ground state of matter at high densities: Equation of state and stellar models K. Iida , K. Sato , Spin-down of neutron stars and compositional transitions in the cold crustal matter . Astrophys. J . 477 , 294 - 312 ( 1997 ). doi: 10.1017/S0074180900115451 http://doi.org/10.1017/S0074180900115451 Spin-down of neutron stars and compositional transitions in the cold crustal matter J. Antoniadis , P.C.C. Freire , N. Wex , et al ., A massive pulsar in a compact relativistic binary . Science 340 , 1233232 ( 2013 ). doi: 10.1126/science.1233232 http://doi.org/10.1126/science.1233232 A massive pulsar in a compact relativistic binary P. Demorest , T. Pennucci , S. Ransom , et al ., A two-solar-mass neutron star measured using Shapiro delay . Nature 467 , 1081 - 1083 ( 2010 ). doi: 10.1038/nature09466 http://doi.org/10.1038/nature09466 A two-solar-mass neutron star measured using Shapiro delay M. Farine , J.M. Pearson , F. Tondeur , Nuclear-matter incompressibility from fits of generalized Skyrme force to breathing-mode energies . Nucl. Phys. A 615 , 135 - 161 ( 1997 ). doi: 10.1016/S0375-9474(96)00453-8 http://doi.org/10.1016/S0375-9474(96)00453-8 Nuclear-matter incompressibility from fits of generalized Skyrme force to breathing-mode energies A.W. Steiner , J.M. Lattimer , E.F. Brown , The equation of state from observed masses and radii of neutron stars . Astrophys. J . 722 , 33 - 54 ( 2010 ). doi: 10.1088/0004-637X/722/1/33 http://doi.org/10.1088/0004-637X/722/1/33 The equation of state from observed masses and radii of neutron stars 10.1007/s41365-017-0329-1 This work was supported in part by the Major State Basic Research Development Program (973 Program) in China (Nos. 2013CB834405 and 2015CB856904), the National Natural Science Foundation of China (Nos. 11625521, 11275125 and 11135011), the Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning, Key Laboratory for Particle Physics, Astrophysics and Cosmology, Ministry of Education, China, and the Science and Technology Commission of Shanghai Municipality (11DZ2260700). 2016-09-29 2016-10-15 2016-10-19 2017-12-05 2017-12-01 2022-08-18 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2017 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-017-0329-1 http://dx.doi.org/10.1007/s41365-017-0329-1 . 2017 Nuclear Science and Techniques 12 28