arxiv: v1 [nucl-th] 11 Jun 2018

Transcription

1 YITP Quark matter symmetry energy effects on hadron-quark coexistence in neutron star matter Xuhao Wu, 1,, Akira Ohnishi,, and Hong Shen 1, 1 School of Physics, Nankai University, Tianjin 371, China Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 66-85, Japan arxiv: v1 [nucl-th] 11 Jun 18 We examine the effects of the isovector-vector coupling and hypercharge-vector coupling in quark matter on hadron-quark coexistence in neutron star matter. The relativistic mean field theory with the TM1 parameter set and an extended TM1 parameter set are used to descrie hadronic matter, and the Namu-Jona-Lasinio model with scalar, isoscalar-vector, isovector-vector and hyperchargevector couplings is used to descrie deconfined quark matter. The hadron-quark phase transition is constructed via the Gis conditions for phase equilirium. The isovector-vector and hyperchargevector couplings in quark matter enhance the symmetry energy and hypercharge symmetry energy in neutron star matter, while their effects are found to e suppressed at high densities y the strange quarks. As a result, the hadron-quark mixed phase shrinks with only isovector-vector coupling and moves to higher density with isovector-vector and hypercharge-vector couplings. The maximum mass of neutron star increases slightly with isovector-vector and hypercharge-vector couplings. PACS numers: 1.65.Qr, 6.6.Dd, 6.6.Kp, h Keywords: hadron-quark phase transition, equation of state, isovector-vector coupling I. INTRODUCTION Existence of the first-order phase transition in dense neutron star matter is of primary interest in nuclear physics and compact star astrophysics. In the inner core region of massive neutron stars, the aryon density may reach 5 1n n.15 fm 3, and the chiral and deconfinement hadron-quark phase transition may occur [1 3]. At present, the sign prolem prevents us to predict the properties of cold dense matter in the first principles calculations such as the lattice QCD Monte- Carlo simulations, then phenomenological approaches as well as experimental and oservational data are necessary to explore the inner part of neutron stars. One of the theoretical approaches to study the phase transition is to apply chiral effective models such as the Namu-Jona-Lasinio NJL model. The first-order phase transition can occur at the aryon chemical potential µ = 1 1 MeV using the chiral effective models. Unfortunately, the order of the transition depends on the parameters, and the equation of state EOS of nuclear matter at low densities is not well descried y the NJL model. Thus it is more realistic to consider the coexistence of hadronic and quark matter in order to predict the transition density. The QCD phase transition signal during the core collapse supernovae was investigated y using a relativistic mean field RMF model for hadronic matter and a ag model for quark matter [4, 5], and the early collapse [4] or the sec- Electronic address: wuhaoird@gmail.com Electronic address: ohnishi@yukawa.kyoto-u.ac.jp Electronic address: shennankai@gmail.com ond shock [5] was found to signal the transition to quark matter. In this case, the transition is necessarily of the first order and the transition density strongly depends on the ag constant. The hadron-quark coexistence is investigated also y using RMF for hadronic matter and the NJL model for quark matter [6, 7]. In RMF, the isovector-vector meson ρ plays an important role to control the symmetry energy in hadronic matter, and the coupling has een carefully chosen to explain the properties of finite nuclei. By comparison, the isovector-vector coupling in quark matter has een considered less carefully. The isovector-vector coupling constant was chosen to e = 1.5G in Refs. [7, 8] and = G in Ref. [9], where G is the isoscalar-vector coupling, while the isovector-vector coupling was ignored in Ref. [6]. This difference comes from the two indpendent types of the [ chiral SUN f vector coupling terms, qγ µ q and α qγµ λ α q + qiγ µ γ 5 λ α q ] [1]. The second type includes the isovector-vector coupling terms and gives rise to symmetry energy in quark matter, the energy increase from the u- and d- quark imalance. Since the symmetry energy in nuclear matter is known to affect the neutron star properties such as radii, it is expected to e important also in the hadron-quark coexistence. In this work, we examine the role of the isovectorvector coupling in quark matter on the hadron-quark coexistence in neutron star matter. For this purpose, we apply RMF for hadronic matter and the three-flavor NJL model for quark matter, and we compare the coexistence density region with and without the isovector-vector coupling term in NJL. A finite isovector-vector coupling in quark matter enhances the symmetry energy in quark matter, which characterizes the increase of the energy per aryon from unalanced u- and d- quark densities. In addition to this isospin symmetry energy, the hypercharge

2 symmetry energy appears with N f = 3 and controls the s- quark contriution. We also examine the effects of nuclear matter symmetry energy slope L. We use the TM1 parameter set L = 11 MeV [11] and an extended TM1 parameter set TM1e [1], where the symmetry energy at the density of n =.11 fm 3 is fixed and the symmetry energy slope is tuned to e L = 5 MeV. We find that the quark matter symmetry energy increases the starting density of the hadron-quark coexistence and enhances the maximum mass of neutron stars. The quark matter symmetry energy effects on the hadron-quark coexistence are suppressed y the s- quarks and are smaller than those of the nuclear matter symmetry energy. This article is organised as follows. In Sec. II, we introduce the RMF model and the NJL model with the isovector-vector coupling for hadronic and quark matter, and we shortly descrie the Gis conditions used as the equilirium conditions in the hadron-quark mixed phase. In Sec. III, we show the numerical results of the hadronquark coexistence in neutron star matter and discuss the impact of the isovector-vector coupling. Section IV is devoted to a summary. II. HADRONIC MATTER, QUARK MATTER AND HADRON-QUARK COEXISTENCE A. Hadronic matter We adopt the RMF theory to descrie the hadronic phase, in which aryons interact y exchanging the isoscalar scalar σ, isoscalar vector ω and isovector vector ρ mesons. These mesons are treated as classical field under the mean-field approximation. For neutron star matter, we use the Lagrangian given as L RMF = ψ i {iγ µ µ M +g σ σ i=p,n γ µ [ g ω ω µ + g ρ τ aρ aµ]} ψ i + 1 µσ µ σ 1 m σ σ 1 3 g σ g 3σ W µνw µν + 1 m ω ω µω µ c 3ω µ ω µ 1 4 Ra µνr aµν + 1 m ρρ a µρ aµ +Λ v g ω ω µ ω µ gρρ a µρ aµ + ψ l iγ µ µ m l ψ l, 1 l=e,µ which contains the contriutions of aryons n and p and leptons e and µ. W µν and R aµν are the antisymmetric field tensors for ω µ and ρ aµ, respectively. The parameters in the Lagrangian are usually determined y fitting nuclear matter saturation properties and groundstate properties of finite nuclei. We use the TM1 parameter set [11] and an extended TM1 parameter set [1], referred to as the TM1e parameter set in later discussions. In TM1e, the symmetry energy slope parameter is tuned to e L = 5 MeV at saturation density, as listed in Tale I. For the homogeneous matter system, the meson field equations have the form as following: m σσ +g σ +g 3 σ 3 = g σ n s p +n s n, m ω ω +c 3ω 3 +Λ v g ω g ρ ρ ω = g ω n p +n n, 3 m ρ ρ+λ vgω g ρ ω ρ = g ρ n p n n, 4 where n s i and n i represent the scalar and vector densities of the i-th aryon i = n, p, respectively. The equations of motion for nucleons give the standard relations etween the densities and chemical potentials, µ p = k p F +M +g ω ω + g ρ ρ, 5 µ n = kf n +M +g ω ω g ρ ρ, 6 where M = M +g σ σ is the effective nucleon mass, and kf i is the Fermi momentum of species i, which is related to the numer density y n i = 3/3π kf i. For neutron star matter in β equilirium, the chemical potentials satisfy the relations µ p = µ n µ e and µ µ = µ e, wherethe chemical potentials of leptons are given y µ l = kf l +m l. In neutron star matter, the total energy density and pressure are given y ε = i=p,n 1 π k i F k +M k dk + 1 m σσ g σ g 3σ m ωω c 3ω m ρ ρ +3Λ v g ω ω g ρ ρ +ε l, 7 P = i=p,n 1 k i F 3π 1 k +M k4 dk 1 m σσ 1 3 g σ g 3σ m ωω c 3ω m ρρ +Λ v g ω ω g ρ ρ +P l, 8 where ε l and P l l = e,µ are the energy density and the pressure from leptons, respectively. B. Quark matter We adopt the three flavor NJL model to descrie the deconfined quark phase. The Lagrangian is given y L NJL = q iγ µ µ m 8 q +G S [ qλ a q + qiγ 5 λ a q ] a= K{det[ q1+γ 5 q]+det[ q1 γ 5 q]}+l V, 9

3 3 TABLE I: Parameters in the TM1 and TM1e parameter sets. The masses are given in the unit of MeV. Model LMeV M m σ m ω m ρ g σ g ω g ρ g fm 1 g 3 c 3 Λ v TM TM1e with L V = G qγ µ q G V 8 α=1 [ qγ µ λ α q + qiγ µ γ 5 λ α q ], 1 in which q denotes the quark field with three flavors N f =3 and three colors N c =3. The determinant interaction is included in order to take account of the U1 A anomaly. G S, G and G V are the scalar, flavor-singletvector and flavor-octet-vector coupling constants, respectively, and have dimensions of energy. In the mean field approximation, only those terms with diagonal matrix elements in λ α remain, then L V is reduced to [ L V = G qγ µ q qγ µ λ 3 q + qiγ µ γ 5 λ 3 q ] G 8 [ qγ µ λ 8 q + qiγ µ γ 5 λ 8 q ]. 11 In the flavor SU3 limit, the isovector-vector coupling and hypercharge-vector coupling G 8 constants should e the same, = G 8 = G V. In order to discuss the isospin symmetry energy and the hypercharge symmetry energy effects separately, we consider the cases with G 8 as well. In the mean-field approximation, quarks get constituent quark masses y spontaneous chiral symmetry reaking, m i = m i 4G S q i q i +K q j q j q k q k, 1 where q i q i C i denotes the quark scalar density, and i,j,k is a permutation of u,d,s. For charge neutral quark matter containing quarks u, d, and s and leptons e and µ in β equilirium, the chemical potentials satisfy the relations µ s = µ d = µ u +µ e and µ µ = µ e, where the chemical potentials of u, d, and s quarks are given y µ u = ε n u = k u F +m u +G n u +n d +n s + n u n d + 3 G 8n u +n d n s, 13 µ d = ε = k d F +m d n +G n u +n d +n s d n u n d + 3 G 8n u +n d n s, 14 µ s = ε = k s F +m n s +G n u +n d +n s s 4 3 G 8n u +n d n s. 15 The energy density of quark matter is given y ε NJL = [ 3 Λ ] π k +m i k dk i=u,d,s k i F +G S C u +C d +C s 4KCu C d C s +G n u +n d +n s + n u n d G 8n u +n d n s ε, 16 where ε is sutracted to ensure ε NJL = in the vacuum. The total energy density and pressure for quark matter are given y ε = ε NJL +ε l, 17 P = n i µ i ε. 18 i=u,d,s,e,µ We employ the parameter set given in Ref. [8], m u = m d = 5.5 MeV, m s = 14.7 MeV, Λ = 6.3 MeV, G S Λ = 1.835, and KΛ 5 = The vector couplings G,,G 8 are considered as free parameters, and we use G = 5G S,,G 8 =, 1.5, 1 G. For larger G values G > 7G S, the energy in quark matter is found to e always larger than that in hadronic matter at high densities, and there is no phase transition. The parameterchoiceof = G 8 = 1.5G correspondstothe Lagrangian adopted in Refs. [7, 8]. The larger isovectorvector coupling, = G 8 = 1G, roughly gives the symmetry energy slope of L Q 5 MeV. The vector couplings increase the energy per aryon as ε V n =9G n + n δ +3G 8 n δ h, 19 δ = n d n u n, δ h = n n s n = B +S B. There are two types of asymmetry parameter, δ and δ h, which are the isospin asymmetry and the hypercharge Y = B+S fraction. The symmetry energy is defined as thecoefficientofδ, thenthevectorcouplingcontriution to the symmetry energy is given as S V n = n = n +3 n n n 3n. 1 The vector coupling contriution to the slope parameter is L V = 3 n = 6.6 and 44 MeV for = 1.5G and 1G, respectively.

4 4 C. Hadron-quark phase transition For neutron star matter, the β equilirium and the charge neutrality are satisfied. We adopt simple Gis conditions for the mixed phase connecting the pure hadronic phase and the pure quark phase. The mixed phase may appear in the inner core of the neutron stars. In the Gis conditions, the gloal charge neutrality condition is given y un c +1 un HP c =, in which u = V /V + V HP represents the volume fraction of quark matter in the mixed phase. The mechanical equilirium requires P HP µ n,µ e = P µ n,µ e. 3 There are two independent chemical potentials, µ n and µ e, the hadronic and quark phases satisfy the chemical equilirium condition, µ u +µ e = µ d = µ s = µ n 3 + µ e 3. 4 With these equilirium constrains, we can solve the mixed phase self-consistently and otain the properties of the hadron-quark mixed phase. III. ISOVECTOR-VECTOR COUPING DEPENDENCE OF HADRON-QUARK COEXISTENCE IN NEUTRON STAR MATTER We shall now investigate the effects of isovector-vector coupling on the equation of state, the density range of hadron-quark coexistence and the properties of neutron stars. For this intent, we use the RMF and NJL models to descrie hadronic and quark matter, respectively. For RMF, the TM1 and TM1e parameter sets are used. These parameter sets show different symmetry energy slopes L, which show significant effects on the neutron star radius. We include isovector-vector coupling in NJL, which modifies the quark matter symmetry energy. In the following discussions, we fix the isoscalar-vector coupling as G = 5 G S, and we compare the results of /G,G 8 /G =,,1.5,,1.5,1.5,1, and 1,1, referred to as NJL-V, NJL-VR1, NJL-VRY1, NJL-VR and NJL-VRY, respectively. NJL-V and NJL-VRY1 corresponds to models in Ref.[6] and Refs.[7, 8], respectively. We define n 1 and n as the starting and the ending aryon densities of the mixed phase. At the density n 1, the energy per aryon in pure hadronic phase ecomes lower than that of the mixed phase. The volume fraction of quark matter u increases with the aryon numer density n, and it transforms into pure quark phase at the density n under the condition that the mixed phase has larger energy density than the pure quark phase. The model dependence of the phase transition densities n 1 and n are summarized in Tale II. A. Equation of state We first discuss the equation of state. In Fig. 1, we show the pressure as a function of the aryon numer density for hadronic matter, quark matter, and the mixed phase otained under Gis conditions. The top and ottom panels show the results from the TM1 L = 11 MeV and TM1e L = 5 MeV parameter sets, respectively. The left, middle and right panels show the results of NJL-V = G 8 =, NJL-VRY1 = G 8 = 1.5G and NJL-VRY = G 8 = 1 G, respectively. In Fig., we compare the energy per aryon, E/A M, as a function of aryon numer density n for the hadronic, mixed, and quark phases. Open circles squares, triangles show the transition densities of NJL- V NJL-VRY1, NJL-VRY. Although the energy difference etween NJL-V, NJL-VR1 and NJL-VR is small in quark matter, the n 1 difference is very visile, whereas the n difference is rather small. The reason comes from the difference of the isospin and hypercharge asymmetry of quarks, δ Q = 3n d n u /n d + n u + n s and δ Q h = n u +n d n s /n d +n u +n s, in the pure quark phase and in the mixed phase. In the mixed phase, the quark matter part is negatively charged and the numer density difference etween u- quark and d- quark so as hypercharge difference is igger than that in the charge neutral pure quark phase as can e found from the electron chemical potential discussed elow. We show the electron chemical potential as a function of the neutron chemical potential the aryon density in the left right panel of Fig. 3. The electron chemical potential µ e in the mixed phase is significantly larger than that in the pure quark matter at the same aryon numer density n. The electron chemical potential µ e reflects the chemical potential difference of the u- and d- quarks µ d µ u = µ e. It can e seen as a signal of the imalance etween u- and d- quarks. The system ecomes more symmetric with decreasing µ e. The ehavior of µ e can explain why the effect of isovector-vector coupling is more significant at lower densities. For pure quark matter, there exists a maximum value of µ e, which corresponds to the appearance of s quark. This change expresses the imalance etween u- and d- quarks is getting smaller. These trends are the same as those found in Fig.. It would e interesting to discuss the reason why the electron chemical potential is small in pure quark matter. In Fig. 3, we find that µ e in pure quark matter increases at low densities, reaches µ e 1 MeV, and turns to decrease at around µ n = 13 MeV. At this density, the quarkchemicalpotentialsareevaluatedasµ u 37MeV and µ d = µ s 47 MeV. Since the chemical potential of s- quark is close to its threshold value for appearance, then we expect that the appearance of s- quarks would e the mechanism to suppress the electron chemical potential.

5 5 TABLE II: Model dependence of the coexistence densities and neutron star properties, G = 5G S. Model L /G G 8/G n 1 n M max/m n Mmax c R1.4 M MeV fm 3 fm 3 fm 3 km TM1/NJL-V TM1/NJL-VR TM1/NJL-VRY TM1/NJL-VR TM1/NJL-VRY TM1e/NJL-V TM1e/NJL-VR TM1e/NJL-VRY TM1e/NJL-VR TM1e/NJL-VRY ± 17 [14, 15] 1 ± 1 [17] Constraints 4 6 [13] ± 4 [16] ± 1. [18] > 14 [19] P MeV/fm 3 P MeV/fm = 4 HP n L=5 MeV = n P MeV/fm 3 P MeV/fm n L=5 MeV n P MeV/fm 3 P GeV/fm n L=5 MeV n FIG. 1: Color online Pressures as a function of the aryon numer density for the pure hadronic phase dash-dotted lines, the mixed phase dashed lines and the pure quark phase solid lines. We find that the mixed phase shrinks when the isovector-vector coupling and hypercharge-vector coupling G 8 are switched on. Two couplings, and G 8, modify oth of the transition densities, and the shift of n 1 is larger than that of n. This difference comes from the density dependence of the isospin asymmetry, δ 1 un n n p +un d n u /n. In the hadronquark mixed phase, the matter tends to e more symmetric with increasing aryon numer density n. At n n 1 <.9fm 3, the isospinasymmetryδ isstill significant, while δ is almost zero at n n >.3 fm 3. Since the energy from the isovector-vector coupling is proportional to δ, it ecomes small at high aryon densities. This mechanism also applies to the effect of the symmetry energy slope L. B. Particle fraction and asymmetry We now discuss the isovector-vectorcoupling and hypercharge-vector G 8 coupling effects in terms of the isospinasymmetryδandhyperchargefractionδ h. We show δ δ h as a function of the aryon numer density in the top ottom two rows in Fig. 4. We first find that and G 8 oviously suppress δ and δ h in the pure quark phaseandthemixedphase. Wealsonotethatδ decreases

6 6 E/A-M MeV 4 3 =G V G V G V G V G HP 3 = 1 HP L=5 MeV 1. n E/A-M MeV HP n FIG. : Color online Energy density per aryon as a function of the aryon numer density for the hadronic, mixed, and quark phases. The left panel is a local enlargement of the right panel. Open circles squares, triangles lael the transition points for = G 8 = = G 8 = 1.5 G, = G 8 = 1 G. e MeV =G VGV = G V G V HP HP L=5 MeV e MeV HP L=5 MeV HP 1 15 MeV 5 1. n n FIG. 3: Color online The left panel shows the relation etween the chemical potential of neutron and electron. The right panel is the electron chemical potential as a function of aryon numer density for different phases. rapidly when s- quark appears and δ h ecomes smaller than unity. With increasing n, the effect from different symmetryenergyslopelecomessmaller. Finite and G 8 increase the energy and delay the appearance of the mixedphase,i.e. twocouplingspushupn 1. By comparison, n is determined y oth the energy increase and the asymmetry decrease in quark matter. The increase of the energy in quark matter pushes up n. The decrease of δ and thus µ e tends to make coexisting hadronic matter isospin symmetric and positively charged, disfavors the mixed phase, and pushes down n. When oth and G 8 are switched on, effects from the energy increase and the asymmetry decrease seem to cancel, and n is almost the same as that without and G 8. When only is switched on, effects from the asymmetry decrease are larger than those from the energy increase, then n decreases slightly. Let us further discuss the effects of quark matter symmetryenergyon particlefractions. In Fig. 5, we showthe particle numer fractions in hadron matter, quark matter and the mixed phase, as functions of the aryonnumer density. As already mentioned, the electron chemi-

7 7 = HP n = HP L=5 MeV n = = HP n = HP L=5 MeV n = h h n L=5 MeV = n L=5 MeV = h h n n FIG. 4: Color online The isospin asymmetry and hypercharge fraction for /G,G 8/G =,,1.5,,1.5,1.5,1, and 1,1.

8 8 i=n,p,u,d,s = Y u Y d Y s i=n,p,u,d,s i=n,p,u,d,s n n n i=n,p,u,d,s L=5 MeV = n i=n,p,u,d,s L=5 MeV n i=n,p,u,d,s L=5 MeV n FIG. 5: Color online Particle numer fractions for the pure hadronic phase dash-dotted lines, the mixed phase dashed lines and the pure quark phase solid lines. cal potential in the mixed phase is larger than that in pure quark matter, and therefore the fractions of d- and s- quarks in the mixed phase are larger than those in quark matter. When the isovector-vector coupling and hypercharge-vector coupling G 8 are taken into account, the differences etween quark fractions ecome smaller and the quark matter tends to e more SU3 symmetric. When s- quark appears, the differences etween u- quark and d- quark is suppressed furthermore. C. Neutron stars Using the EOS of pure hadronic matter, the TM1 parameter set predicts a maximum neutron-star mass of.18 M, with M eing the solar mass. The neutron star oservations of PSR J [14, 15] and PSR J [16] constrain that the neutron star maximum massneeds to e largerthan M. To examinethe effect of isovector-vector coupling and hyperchargevector coupling G 8 on the properties of neutron stars, we solve the Tolman-Oppenheimer-Volkoff TOV equation y using models listed in Tale II. The mass-radius relation is presented in the left panel of Fig. 6. The right panel shows the neutron star mass as a function of the neutron star central density n c. We find that the isovector-vector coupling enhances the maximum mass of neutron stars slightly, ut the effect is inconspicuous. For the results of TM1e L = 5 MeV, it shows smaller radius than that of TM1 L = 11 MeV. The influence of on neutron star maximum mass ecomes even smaller with TM1e. This is ecause the onset of the mixed phase in TM1e is later than that in TM1. In the right panel of Fig. 6, we can see that when the transition density n 1 is close to the maximum neutron star central density n Mmax c for TM1e case, the hadron-quark coexistence has little effect for the maximum mass neutron stars. IV. SUMMARY Effects of the isovector-vector and hypercharge-vector couplings in quark matter on hadron-quark phase transition and neutron star properties are investigated. In this work, we have used the RMF theory to descrie hadronic matter, and the three flavor NJL model including the isovector-vector coupling has een used for the quark matter. The Gis conditions are applied to descrie the hadron-quark mixed phase. We have found that the mixed phase shrinks with the isovectorvector coupling in quark matter, while the mixed phase moves to higher density with oth isovector-vector and hypercharge-vector couplings included. If only isovectorvector coupling is included in quark matter, the transition density to the mixed phase n 1 increases y 6 5fm 3, and the transition density to the pure quark matter n decreases y 4.14fm 3. The hypercharge-vector coupling delays oth n 1 and n. The inclusion of the isovector-vector and hyperchargevector couplings in quark matter has similar effects as decreasing the symmetry energy slope L in hadronic matter. Both of them can affect the asymmetry of the system. We have found that the isovector-vector and hypercharge-vector couplings suppress the asymmetry of the nuclear-quark matter system. Meanwhile, the softening of the EOS due to the hadron-quark phase transition ecomes weaker and the neutron star maximum mass increases without modifying the neutron star radius around

9 9.5.5 M/M. 1.5 L=5 MeV.5 G 3 = R km M/M L=5 MeV n c FIG. 6: Color online The left panel plots mass-radius relations of neutron stars for different EOS. The solid lines and the dash-dotted lines dashed lines show the results for without and with isovector-vector coupling and hypercharge-vector coupling respectively. The right panel shows maximum mass as the function of neutron star central density. M 1.4M, wherethecentraldensityisn 3n. Acknowledgment This work is supported in part y the Grants-in- Aid for Scientific Research from JSPS Nos. 15K579, 15H3663, 16K535, the Grants-in-Aid for Scientific Research on Innovative Areas from MEXT Nos. 4151, 4158, the National Natural Science Foundation of China No , and y the Yukawa International Program for Quark-hadron Sciences YIPQS. Xuhao Wu acknowledges the financial support provided y the China Scholarship Council CSC. [1] N. K. Glendenning, Phys. Rep. 34, [] H. Heiselerg and M. Hjorth-Jensen, Phys. Rep. 38, 37. [3] F. Weer, Prog. Part. Nucl. Phys. 54, [4] K. Nakazato, K. Sumiyoshi, and S. Yamada, Phys. Rev. D 77, [5] I. Sagert, T. Fischer, M. Hempel, G. Pagliara, J. Schaffner-Bielich, A. Mezzacappa, F.-K. Thielemann and M. Lieendorfer, Phys. Rev. Lett. 1, [6] S. Benic, D. Blaschke, D. E. Alvarez-Castillo, T. Fischer and S. Typel, Astron. Astrophys. 577, A4 15 [7] X. H. Wu and H. Shen, Phys. Rev. C 96, [8] P. Reherg, S. P. Klevansky, and J. Hüfner, Phys. Rev. C 53, [9] H. Ueda, T. Z. Nakano, A. Ohnishi, M. Ruggieri, and K. Sumiyoshi, Phys. Rev. D 88, [1] Peng-Cheng Chu, Bin Wang, Hong-Yang Ma, Yu-Min Dong, Su-Ling Chang, Chun-Hong Zheng, Jun-Ting Liu, and Xiao-Min Zhang, Phys. Rev. D 93, [11] Y. Sugahara and H. Toki, Nucl. Phys. A 579, [1] S. S. Bao, J. N. Hu, Z. W. Zhang, and H. Shen, Phys. Rev. C 9, [13] I. Tews, J. M. Lattimer, A. Ohnishi, and E. E. Kolomeitsev, Astrophys. J. 848, [14] P. B. Demorest, T. Pennucci, S. M. Ranson, M. S. E. Roerts, and J. W. T. Hessels, Nature 467, [15] E. Fonseca et al., Astrophys. J. 83, [16] J. Antoniadis et al., Science 34, [17] J. Nättilä, A. W. Steiner, J. J. E. Kajava, V. F. Suleimanov, and J. Poutanen, Astronomy & Astrophysics 591, A5 16 [18] S.Guillot, and R. E. Rutledge, Astrophys. J. 796, L3 14. [19] V. Suleimanov, J. Poutanen, M. Revnivtsev and K. Werner, Astrophys. J. 74, 1 11.