Apparent softening of the symmetry energy with the inclusion of non-nucleonic components in nuclear matter

Apparent softening of the symmetry energy with the inclusion of non-nucleonic components in nuclear matter 1 corresp first-author Corresponding author. E-mail address:wzjiang@seu.edu.cn Corresponding author. E-mail address:wzjiang@seu.edu.cn wzjiang@seu.edu.cn Department of Physics, Southeast University, Nanjing 211189, China Department of Physics, Southeast University, Nanjing 211189, China Apparent softening of the symmetry energy with the inclusion of hyperon and quark degrees of freedom is demonstrated by the fact that the phase transition causes the change of the interaction and the suppression of nucleon fractions. The demonstration is fulfilled in the relativistic mean-field model. Symmetry energy Nuclear matter Relativistic mean-field model 1 Introduction 1 The nuclear symmetry energy of isospin asymmetric nuclear matter is important to understand the structure of neutron- or proton-rich nuclei, the reaction dynamics of heavy-ion collisions [ 1 1 , 2 2 , 3 3 ] , and many astrophysical issues [ 4 4 , 5 5 , 6 6 ] as well. In the past four decades, the properties of symmetric matter have been constrained rather satisfactorily, while just in the recent decade, appreciable progress has been achieved on constraining the symmetry energy at saturation and subsaturation densities either through the extraction based on astrophysical observations or in terms of terrestrial data [ 7 7 , 8 8 , 9 9 , 10 10 , 11 11 ] . However, the density dependence of the symmetry energy is still poorly known especially at supra-normal densities [ 3 3 , 12 12 , 13 13 , 14 14 ] . The density dependence of the symmetry energy is rather diverse according to theoretical predictions. It may increase nonlinearly or linearly with the density. Strikingly, some non-relativistic models even predict that the symmetry energy soon goes to negative values at densities several times normal density. This can probably cause severe problems in stabilize the structure of neutron stars. On the other hand, similarly diverse symmetry energy was extracted from analyzing the FOPI/GSI data on the π − / π + ratio in relativistic heavy-ion collisions with various transport models [ 12 12 , 13 13 , 14 14 ] . Recently, the required coincidence with data from the ALADIN-2000 collaboration, analyzed by Kumar et al. [ 15 15 ] , suggested a soft symmetry energy that differs from the super-soft one obtained from the analysis of FOPI/GSI data. More discussions on the status quo of the experimental extraction can be found in a recent work [ 16 16 ] . These extractions really reflect the fact that the density dependence of the symmetry energy is still very illusory, provided the equation of state of symmetric matter used in the transport models is well constrained. Regardless of the inconsistency in the experimental extractions, the theoretical uncertainty of high-density symmetry energy is regarded to be associated with the tensor force that originates from the exchange terms [ 17 17 , 18 18 ] . In the ladder approximation, the exchange terms can be well treated in the Brueckner theory either in the relativistic or non-relativistic frameworks [ 19 19 , 20 20 ] . However, the mean-field approximation without exchange terms works more sophisticatedly in high-density matter. It is thus incomprehensible that the tensor force can completely explain the marvellous softening of the high-density symmetry energy. In this work, we decline the strategy using the tensor force to solve the symmetry energy divergence but explore the effect of the non-nucleonic degrees of freedom on the symmetry energy. The results reported here are performed in the relativistic Hartree approximation. 2 Modelling 1 In this work, we consider the non-nucleonic degrees of freedom which include hyperons and quarks. On the hadron level, we make use of the relativistic mean-field (RMF) models. The original Lagrangian of the RMF model was first proposed by Walecka 40 years ago [ 21 21 ] . The Walecka model and its improved versions were characteristic of the cancellation between the big attractive scalar field and the big repulsive vector field. The success of the RMF models is partially attributed to its dynamical description for the spin-orbit interaction. On the quark level, we adopt the MIT bag model [ 22 22 ] to portray the quark phase. The mixed phase is built on the mechanical and chemical equilibriums according to Gibbs conditions. In the parabolic approximation, the energy per nucleon in isospin asymmetric nuclear matter can be written as <math display="block" id="M1"><mrow><mi>ε</mi><mo>/</mo><mi>ρ</mi><mo>=</mo><mi>E</mi><mo>/</mo><mi>A</mi><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><mi mathvariant="normal">sym</mi></mrow></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><msup><mi>δ</mi><mn>2</mn></msup><mo>,</mo></mrow></math> where <math id="M2"><mrow><msub><mi>e</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow></math> is the energy per nucleon in symmetric nuclear matter, <math id="M3"><mrow><msub><mi>E</mi><mrow><mi mathvariant="normal">sym</mi></mrow></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow></math> is the density dependence of symmetry energy, and <math id="M4"><mrow><mi>δ</mi><mo>=</mo><mo stretchy="false">(</mo><msub><mi>ρ</mi><mi>n</mi></msub><mo>−</mo><msub><mi>ρ</mi><mi>p</mi></msub><mo stretchy="false">)</mo><mo>/</mo><mi>ρ</mi></mrow></math> is the isospin asymmetry. The symmetry energy in the RMF models reads, <math display="block" id="M5"><mrow><msub><mi>E</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msubsup><mi>C</mi><mi>ρ</mi><mn>2</mn></msubsup><mi>ρ</mi><mo>+</mo><mfrac><mrow><msubsup><mi>k</mi><mi>F</mi><mn>2</mn></msubsup></mrow><mrow><mn>6</mn><msub><mi>E</mi><mi>F</mi></msub></mrow></mfrac><mo>,</mo></mrow></math> where <math id="M6"><mrow><msub><mi>C</mi><mi>ρ</mi></msub></mrow></math> is the ratio of the ρ coupling constant with the nucleon to the ρ meson effective mass, k F and E F are the nucleon Fermi momentum and energy, respectively. In the presence of hyperons, it reads <math display="block" id="M7"><mrow><msub><mi>E</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mi>B</mi></msub><mo stretchy="false">)</mo><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msubsup><mi>C</mi><mi>ρ</mi><mn>2</mn></msubsup><mfrac><mrow><msubsup><mi>ρ</mi><mi>N</mi><mn>2</mn></msubsup></mrow><mrow><msub><mi>ρ</mi><mi>B</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><msubsup><mi>k</mi><mi>F</mi><mn>2</mn></msubsup></mrow><mrow><mn>6</mn><msub><mi>E</mi><mi>F</mi></msub></mrow></mfrac><mfrac><mrow><mi>ρ</mi><mrow/><mi>N</mi></mrow><mrow><msub><mi>ρ</mi><mi>B</mi></msub></mrow></mfrac><mo>,</mo></mrow></math> where <math id="M8"><mrow><msub><mi>ρ</mi><mi>B</mi></msub></mrow></math> and <math id="M9"><mrow><msub><mi>ρ</mi><mi>N</mi></msub></mrow></math> are the baryon and nucleon density, respectively. As the quarks appear, the system is first in the mixed phase and then in pure quark phase at very high densities. The symmetry energy can be defined according to the similar parabolic approximation of the equation of state, <math id="M10"><mrow><mi>ε</mi><mo>/</mo><msub><mi>ρ</mi><mi>B</mi></msub><mo>=</mo><msub><mi>e</mi><mn>0</mn></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mi>B</mi></msub><mo>,</mo><mi>Y</mi><mo stretchy="false">)</mo><mo>+</mo><msubsup><mi>E</mi><mrow><mtext>sym</mtext></mrow><mi>H</mi></msubsup><mo stretchy="false">(</mo><msub><mi>ρ</mi><mi>B</mi></msub><mo>,</mo><mi>Y</mi><mo stretchy="false">)</mo><msubsup><mi>δ</mi><mi>H</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>E</mi><mrow><mtext>sym</mtext></mrow><mi>Q</mi></msubsup><mo stretchy="false">(</mo><msub><mi>ρ</mi><mi>B</mi></msub><mo>,</mo><mi>Y</mi><mo stretchy="false">)</mo><msubsup><mi>δ</mi><mi>Q</mi><mn>2</mn></msubsup><mo>,</mo></mrow></math> where H and Q denote hadrons and quarks, respectively, and Y is the quark phase proportion. 3 Results 1 For simplicity, we just consider the Λ hyperon which usually occupies a majority of hyperons in hyeronized matter. The symmetry energy for various Λ fractions is calculated in symmetric matter at δ =0. The effect of Λ hyperons on the nuclear symmetry energy with the RMF model NL3 [ 23 23 ] is illustrated in Fig.1 Fig.1 . It is shown in Fig.1 Fig.1 that the symmetry energy is softened clearly with the increase of the Λ fraction. Compared with Eqs.(2) and (3), we see that the softening is dominated by the suppression factor <math id="M11"><mrow><msub><mi>ρ</mi><mi>N</mi></msub><mo>/</mo><msub><mi>ρ</mi><mi>B</mi></msub></mrow></math> . However, even without this suppression factor, the symmetry energy in hyperonized matter is still modified by the isoscalar Λ hyperons provided there exists the isoscalar-isovector coupling Λ v (for model details, see Ref.[ 2 2 ]). This is clearly seen in the inset of the lower panel of Fig.1 Fig.1 where the potential part of the symmetry energy is displayed. Similarly, if the charged hyperons are included, the potential part of the symmetry energy can be modified even without the isoscalar-isovector coupling. This can be verified numerically, while it is beyond the scope of the present work. Nevertheless, we may infer that the symmetry energy may significantly be modified by taking into account the hyperons. With the increase of density, the hadron-quark phase transition may occur. In this work, quark matter, regarded as the free fermion gas without interactions, is described by the MIT bag model [ 22 22 ] . The mixed phase consists of high-density quark matter and low-density nuclear matter with the quark phase proportion Y being obtained according to Gibbs conditions. The quark phase proportion Y depends on the isospin asymmetry. In this way, the symmetry energy in the mixed phase obtained in symmetric matter cannot simply be used to predict the properties of asymmetric matter because the quark phase proportion changes with the isospin asymmetry in asymmetric matter. Nevertheless, the symmetry energy obtained in symmetric matter is instructive to exhibit its variation tendency in the mixed phase. Shown in Fig.2 Fig.2 is the nuclear symmetry energy as a function of baryon density using the RMF models SLC and SLCd [ 24 24 , 25 25 , 26 26 ] combined with the MIT bag model with the bag constant B =(160MeV) 4 (upper panel) and B =(180MeV) 4 (lower panel). Apparent decrease of the symmetry energy can be observed after the hadron-quark phase transition occurs. With the increase of density, the nucleon proportion decreases, this causes an apparent reduction of the nuclear symmetry energy. As the nucleon proportion reduces to zero, the nuclear symmetry energy vanishes. Appreciably, it is seen that the change of the symmetry energy is very sensitive to the bag constant. Moreover, in Fig.2 Fig.2 we compare two cases with and without hyperons. It is found that the inclusion of hyperons further softens the symmetry energy as the hadron-quark phase transition occurs. It is worth noting that the softening of the symmetry energy in the mixed phase is mostly apparent because once the quark phase proportion can be identified at given densities the nuclear symmetry energy would be extracted appropriately by singling out the effect of suppression factor (1− Y ). However, the determination of the Y is strongly model-dependent and far from experimental feasibility. Thus, the extraction of the high-density symmetry energy for pure nucleonic matter is not well grounded once the hadron-quark phase transition takes place. Most likely, the high-density symmetry energy extracted from the heavy-ion collisions would be as soft as that presented in this work because of the absence of a reliable discrimination or calibration for dynamically evolutional matter. 4 Conclusion 1 We have reported the effect of Λ hyperons and quarks on the nuclear symmetry energy at high densities with the RMF models combined with the MIT bag model. The softening of the nuclear symmetry energy is observed in nuclear matter at given Λ fractions. In the presence of the hadron-quark phase transition, the nuclear symmetry energy obtained in the mixed phase reduces quickly with the rise of quark phase proportion. We argue that this softening is apparent but is most likely realistic according to the capacity of the experimental detection. References Li B A . Phys Rev Lett , 2000 , 85 : 4221 - 4224 ; Phys Rev Lett Li B A . Phys Rev Lett , 2002 , 88 : 192701 . Phys Rev Lett Jiang W Z and Zhao Y L . Phys Lett B , 2005 , 617 : 33 - 39 . Phys Lett B Li B A , Chen L W , Ko C M . Phys Rep , 2008 , 464 : 113 - 281 . Phys Rep Lattimer J M and Prakash M . Phys Rep , 2000 , 333 : 121 - 146 ; Phys Rep Lattimer J M and Prakash M . Phys Rep , 2007 , 442 : 109 - 165 . Phys Rep Horowitz C J and Piekarewicz J . Phys Rev Lett , 2001 , 86 : 5647 - 5610 . Phys Rev Lett Steiner A W , Lattimer J M , Brown E F . Astrophys J , 2010 , 722 : 33 - 54 . Astrophys J Chen L W , Ko C M , Li B A , et al . Phys Rev C , 2010 , 82 : 024321 . Phys Rev C Newton W , Gearheart M , Li B A . 2011 , arXiv:1110.4043 . arXiv:1110.4043 Tsang M B , Stone J R , Camera F , et al . Phys Rev C , 2012 , 86 : 015803 . Phys Rev C Wen D H , Newton W G , Li B A . Phys Rev C , 2012 , 85 : 025801 . Phys Rev C Lattimer J M and Lim Y . 2012 , ArXiv:1203.4286v1 . ArXiv:1203.4286v1 Xiao Z G , Li B A , Chen L W , et al . Phys Rev Lett , 2009 , 102 : 062502 . Phys Rev Lett Feng Z Q and Jin G M . Phys Lett B , 2010 , 683 : 140 - 144 . Phys Lett B Russotto P , Trauntmann W , Li Q F , et al . Phys Lett B , 2011 , 697 : 471 - 476 . Phys Lett B Kumar S and Ma Y G . Phys Rev C , 2012 , 86 : 051601(R) . Phys Rev C Kumar S , Ma Y G , Zhang G Q , et al . Phys Rev C , 2012 , 85 : 024620 . Phys Rev C Xu C and Li B A . Phys Rev C , 2010 , 81 : 064612 . Phys Rev C Vidana I , Polls A , Providencia C . Phys Rev C , 2011 , 84 : 062801(R) . Phys Rev C Brockmann R and Machleidt R . Phys Rev C , 2000 , 42 : 1965 - 1980 . Phys Rev C Song H Q , Baldo M , Giansiracusa G , et al . Phys Rev Lett , 1998 , 81 : 1584 - 1587 . Phys Rev Lett Walecka J D . Ann Phys (NY) , 1974 , 83 : 491 - 501 . Ann Phys (NY) Chodos A , Jaffe R L , Johnson K , et al . Phys Rev D , 1974 , 9 : 3471 - 3495 . Phys Rev D Lalazissis G A , Konig J , Ring P . Phys Rev C , 1997 , 55 : 540 - 543 . Phys Rev C Jiang W Z , Li B A , Chen L W . Phys Rev C , 2007 , 76 : 054314 . Phys Rev C Jiang W Z , Li B A , Chen L W . Phys Lett B , 2007 , 653 : 184 - 189 . Phys Lett B Jiang W Z , Li B A , Chen L W . Astrophys J , 2012 , 756 : 56 . Astrophys J 10.13538/j.1001-8042/nst.2013.05.007 Supported by National Natural Science Foundation of China (NSFC) projects (Nos.10975033 and 11275048) . 2013-06-27 2013-10-01 2021-04-11 2013 Nuclear Science and Techniques 5 24