Arxiv:1308.2130V1 [Nucl-Th] 9 Aug 2013 Eea H He-Aynrplin(TBR)

Home , Pomeron

arXiv:1308.2130v1 [nucl-th] 9 Aug 2013 eea h he-aynrplin(TBR). repulsion three-baryon called the baryons, as recovered exis- three here the among be is repulsions can concern universal basic of EoS Our tence the assumption. stiffening this for by clearly TNR the of role ueta h N-yerplin okuieslyfor universally work conjec- repulsions the introduced TNR-type YN [2] the (Y.Y.) that authors ture the Takatsuka of In Nishizaki, one problem, mass. and serious maximum this the avoid for to effect order EoS, TNR the the of cancels softening which remarkable the about brings matter inpoet rcsl,i sncsayt nrdc also phenomenologically. introduce part to this necessary treat We is by satura- TNA. it nuclear account precisely, the into reproduce property to taken tion order In is the this modeling. within TBR means ESC (MPP) ESC potential work, exchange Hereafter, multi-pomeron this the In 4]. [3, ESC08c version. model ESC as of named version latest is one- The usual no models. the and exchange from boson differently explicitly included meson- account is and into boson two-meson effective taken Extended Here, are exchanges the models. pair by (ESC) accomplished Core Soft been has development eto tr.Telreosre au fneutron-star of value observed 1 large mass The of problem stars. maximum-mass neutron the on impact great a brought bevdmxmmms fnurnstars. neutron the of assuring for mass matter, role maximum neutron-star of observed essential EoS an the in of plays contribution stiffening matter values TNR high neutron The a high-density to incompressibility: leads nuclear region of high-density the increasing in contribution rapidly repulsive TNR the the and Especially, (TNA) (TNR). part part attractive the of composed ae nudryn ula neatos nessential ( an three-nucleon the interactions, by played nuclear is role underlying on matter. neutron-star based of (EoS) state of equation the of nmdln for modeling In h e bevto ftesse 11-20[]has [1] J1614-2230 system the of observation new The oee,tehprn( hyperon the However, o elzto ftencerstrto property saturation nuclear the of realization For and . 97 YY M solar scluae ysligteTVeuto.Temxmmms i mass maximum contr The MPP the equation. TOV 1 including one the obtained observed solving is by matter calculated neutron is in sc nucleus-nucleus EoS the univers The analyzing works by which determined matter, is strength neutron-star in indicated sion ASnmes 16.d 57.z 26.60.Kp 25.70.-z, 21.65.Cd, numbers: PACS swl sfor as well as ut-oeo xhneptnil(P)i rpsda mo a as proposed is (MPP) potential exchange multi-pomeron A .INTRODUCTION I. ie eeecniinfrtestiffness the for condition severe a gives YN 1 cioeiNtoa olg fTcnlg,Ihnsk,Iw Ichinoseki, Technology, of College National Ichinoseki ut-oeo euso n h eto-trmass Neutron-star the and repulsion Multi-pomeron .Yamamoto Y. and 1 n hmclRsac RKN,Wk,Siaa 5-18 Ja 351-0198, Saitama, Wako, (RIKEN), Research Chemical and . ihn etrfrAclrtrBsdSine nttt f Institute Science, Accelerator-Based for Center Nishina 97 NN 4 Y M YY MP,Uiest fNjee,Njee,TeNetherlands The Nijmegen, Nijmegen, of University IMAPP, iigi neutron-star in mixing ) solar hysoe htthe that showed They . 3 eateto hsc,CiaIsiueo Technology, of Institute Chiba Physics, of Department neatos important interactions, -- hbzn aahn,Cia2502,Japan 275-0023, Chiba Narashino, Shibazono 2-1-1 N ntebsso h xeietldata. experimental the of basis the on neato (TNI) interaction ) 1 .Furumoto T. , 2 .Yasutake N. , r sdt nlz the The analyze MPP: to used in contributions are MPP constants including potentials coupling folding the G-matrix determine to here aino h incompressibility the on mation ex- the reproduce to as data. so perimental adjusted are MPP of strengths eatos nrf[] twssonceryta h TNR ( the of scattering that tic distributions clearly angular shown in was appeared effect it ref.[6], In teractions. o,eog ogv h bevdnurnsa mass. neutron-star observed of stiff the a EoS give in to result the enough to constructing expected EoS, is in which included matter, neutron-star MPP are G-matrix The model the with theirs: folding analyses from by development determined a contributions can as approach considered Our be 8]. [7, interactions density-dependent n aacnb sda nidpnetmto odeter- to method independent of an value as the used mine be scatter- can that nuclear data out ing high-precision pointed was of it remain analyses hand, still folding-model other EoS the the On for inconclusive. results In the however, collisions. cases, heavy-ion central many high-energy in formed ter evbe r ieyrpoue ihtecmlxG- complex the free-space from with derived potentials reproduced folding matrix nicely ob- are scattering servables nucleus-nucleus and nucleon-nucleus that CON the with theory matter G-matrix choice. nuclear lowest-order high-density the of using properties study can one n ihrhl-iecnrbtosu o3 to up density contributions normal includ- hole-line results higher the well ing simulate to particle shown single were intermediate potentials the for with choice calculations (CON) continuous G-matrix nucleon lowest-order of the case the matter, In interactions. renor- G-matrix are into components malized tensor and free-space short-range correlations by of the induced basis Here, the models. on interaction baryon-baryon systems many-body of studies ayatmt aebe aet xrc oeinfor- some extract to made been have attempts Many h motn etr forMPi hti works it that is MPP our of feature important The n fgetscesso h -arxapoc is approach G-matrix the of successes great of One h -arxter ie odsatn on for point starting good a gives theory G-matrix The teigwt h -arxfligmodel. folding G-matrix the with attering lyaogtre n orbros Its four-baryons. and three- among ally 3 n hA Rijken Th.A. and , bandt elre hnthe than larger be to obtained s bto.Tenurnsa mass neutron-star The ibution. E/A t,0181,Japan 021-8511, ate, ρ e o h he-oyrepul- three-body the for del 0 K 7 MeV), =70 5.O h ai fti recognition, this of basis the On [5]. hc a eosrtdwt s of use with demonstrated was which , rPhysical or 16 O+ pan 16 41 cteig n hnthe then and scattering, O etc K hi nlssi used is analysis Their . fhg-est mat- high-density of 16 ∼ O+ ie of times 4 16 NN elas- O in- 2 universally not only in NNN states but also YNN, YYN elastic scatterings were analyzed with the use of double- and YYY states. It will be shown in our future work that folding potentials. In this work, the TNR effect was as- inclusions of MPP’s in hypernuclear calculations lead to sumed to come from the medium effect on the vector- reasonable results. meson masses differently from the present MPP model. When a scattering energy is high enough, the frozen- density approximation gives a good prescription. In this II. MULTI-POMERON POTENTIAL approximation, G-matrices at about two times of normal density contribute to folding potentials. Then, the We introduce the universal TBR consistently with the TNR effects were shown to appear clearly in angular dis- ESC modeling, assuming that the dominant mechanism tributions in the case of 16O+16O elastic scattering at is triple and quartic pomeron exchange. For the N-tuple E/A = 70 MeV, and also in the cases of 16O elastic scat- pomeron vertex generally, we take the Lagrangian tering by the 12C, 28Si and 40Ca at E/A = 93.9 MeV, and 12C+12C elastic scattering at E/A = 135 MeV. Here, the (N) 4−N N N = gP σP (x)/N! (1) same analyses are performed so that the MPP strengths L M (3) (4) The N-body local potential by pomeron exchange is gP and gP are determined to reproduce the experimental data with the use of the G-matrix folding potential x′ x′ x x x x N x′ x V ( 1, ..., N ; 1, ..., N ) V ( 1, ..., N ) Πi=1δ( i i), derived from ESC including MPP. ≡ − On the other hand, the nuclear saturation property 3 cannot be reproduced only by adding MPP to ESC. N d ki − k ·x x x ( ) N i i i V ( 1, ..., N )= g gP e Then, we introduce also a TNA part phenomenologically P Z (2π)3 as a density-dependent two-body interaction N 3 k N k2 4−3N 0 2 (2π) δ( i)Π exp , (2) VTNA(r; ρ)= V exp( (r/2.0) ) ρ exp( 4.0ρ) . (4) × i=1 − i M TNA − − Xi=1 Here, because the functional form is not determined where the (low-energy) pomeron propagator is the same within our analysis, it is fixed to be similar to the TNA as used in the two-body pomeron potential. Since the 0 part given in [12]. Only the strength VTNA is treated as pomeron is an SU(3)-singlet, MPP’s work universally an adjustable parameter. among baryons. This TNA is assumed to work only in even states: An The effective two-body potential in a baryonic odd-state TNA contributes to the 16O+16O potential at medium is obtained by integrating over the coordinates E/A = 70 MeV more than the nuclear-matter energy, x x 3, ..., N . This gives because G-matrices in the former are with higher mo- menta than those in the latter. Therefore, they cannot (N) x x N−2 3 3 x x x Veff ( 1, 2)= ρNM d x3... d xN V ( 1, 2, ..., N ) be reproduced consistently, for instance, with use of a Z Z Wigner-type TNA. Also the Fujita-Miyazawa TNA used N−2 3 (N) N ρNM 1 mP 1 2 2 in [6] is specified by dominant contributions from even- = g gP − exp mP r ,(3) P 3N 4 π√π √2 −2 12 state components. M On the basis of G-matrix calculations, strengths of ρNM being a nuclear-matter density. We restrict our- (3) (4) 0 the TNR part (gP and gP ) and the TNA part (VTNA) selves here to the triple and quartic pomeron couplings. 16 16 (3) are determined so as to reproduce the O+ O angu- Values of MPP strengths g can be estimated from P lar distribution E/A = 70 MeV and the minimum value the experimental cross sections of the process pp pX → 16 MeV of the energy per nucleon at normal density (diffractive production of showers of particles) at very ∼− in symmetric matter. Because the ratio of g(3) and g(4) high energies [9]: The estimated values are g(3)/g = P P P P cannot be determined in our analysis, it is taken ade- 0.15 0.20. Using the value of gP = 3.67 used in ESC, ∼ quately referring the above estimation based on [9, 10]. we have (3) (4) 0 Then, we have (g , g )= (2.34, 30.0) and VTNA = 36 (3) P P − gP =1.95 2.6 . MeV. This set is called as MP1a. In order to see a role ∼ of g(4) especially in high-density region, we take here an- (4) P In the Reggeon field theory [10], the value of gP can be (3) (4) other set (gP , gP )= (2.94, 0.0) with the same value of estimated as 0 (4) VTNA. This set without gP is called MP2a. (4) (3) 2 gP = (8.8 60)(gP ) = 33 228 . In Fig.1, the calculated results of the differential cross ∼ ∼ sections for the 16O+16O elastic scattering at E/A = 70 MeV are compared with the experimental data [11]. III. MPP IN NUCLEUS-NUCLEUS Here, the dotted curve is obtained with the G-matrix SCATTERING folding potential derived from MP0 meaning ESC without MPP, which deviates substantially from the data. In [6] the complex G-matrix interaction was derived Solid and dashed curves are for MP1a and MP2a, re- from the ESC model including TNR, and nucleus-nucleus spectively, which are found to fit the data nicely. Though 3

100 10

10−2 E/A[MeV] Ruth. -10

d -20

0.1 0.2 0.3 −4 16O + 16O elastic scattering 10 [fm-3] E/A = 70 MeV FIG. 2: Energy per particle of symmetric matter as a function 0 10 20 θ of ρN . Solid and dashed curves are for MP1a and MP2a), c.m. (deg) respectively. Dotted curve is for MP0. The box shows the empirical value. FIG. 1: Diﬀerential cross sections for 16O+16O elastic scattering at E/A = 70 MeV calculated with the G-matrix folding potentials. Solid and dashed curves are for MP1a and MP2a, in nice agreement to the values E = 32.5 0.5 MeV respectively. Dotted curve is for MP0. sym and L = 70 15 MeV determined recently on± the basis of experimental± data [15]. Thus, the nuclear saturation property derived from MP1a or MP2a is quite reason- reduction factors are often multiplied on the imaginary able in comparison with the empirical values, as seen in parts in the folding model analyses [6], such a reduction obtained values of K, E and L. factor is not used in the present analysis. The necessity sym to include the quartic pomeron coupling has to appear in 140 the diﬀerence between results for MP1a and MP2a, but it cannot be found in the present analyses for nucleus- 120 nucleus scattering. 100 80

60 IV. EQUATION OF STATE E/A[MeV] 40

20 In Fig.2, energies per nucleon in symmetric nuclear 0 matter are drawn as a function of nucleon density ρN . 0.1 0.2 0.3 0.4 0.5 The box in the figure show the area where nuclear satu- [fm-3] ration is expected to occur empirically. The dotted curve n is obtained for MP0. Then, the saturation density and FIG. 3: Energy per particle of neutron matter as a function the minimum energy in symmetric matter are found to of ρn. Solid and dashed curves are for MP1a and MP2a, be deviated substantially from the box. Solid and dashed respectively. Dotted curve is for MP0. Dot-dashed curve is curves are for MP1a and MP2a, respectively. As seen for UIX [16]. clearly, saturation densities and minimum values in these cases are close nicely to the empirical value shown by the box. Then, obtained values of incompressibility K are 280 MeV and 270 MeV at minimum point in the cases V. NEUTRON STAR MASS of including MP1a and MP2a, respectively. which are comparable to the empirical value 240 20 MeV [13, 14]. ± In Fig.3, we show energy curves of neutron mat- Recently, the important roles of the three-neutron ter, namely binding energy per neutron as a function (3n) repulsion are studied for the EoS of neutron matter and the neutron-star mass [16]. Among their interac- of neutron density ρn. Solid and dashed curves are for MP1a and MP2a, respectively. Dotted curve is for tions, the stiffest EoS is given by the UIX model (the MP0. The difference between the energy curves for neu- Argonne AV8’ model combined with the Urbana IX). tron matter and symmetric matter gives the symmetry This EoS is stiff enough to explain the large neutron-star energy E (ρ), and its slope parameter is defined by mass 1.97Msolar [1]. Their result for the UIX model is sym shown by the dot-dashed curve in Fig.3. Our results for ∂Esym(ρ) L = 3ρ0 ∂ρ . The values of Esym at normal den- MP1a (MP2a) turns out to be stiffer than (comparable sity areh 32.5, 33.4i and 33.5 MeV in the cases of MP0, to) theirs. MP1a and MP2a, respectively, and the values of L are Using the EoS of pure neuron matter, we solve the 71, 72 and 73 MeV correspondingly. These values are Tolmann-Oppenheimer-Volkoff (TOV) equation for the 4 hydrostatic structure of a spherical non-rotating star, and 3.0 obtain the mass and radius of neutron stars. From the 2.5 G-matrix calculations, we obtain the energy per neutron as a function of density ρ. According to [16], the energy 2.0

curve is parameterized as solar 1.5

α β M/M ρ ρ 1.0 E(ρ)/A = a + b . (5) ρ0 ρ0 0.5 The ﬁtting parameters are given in Table I in the cases 0.0 of MP0, MP1a and MP2a. These EoS’s are used for 10 12 14 16 ρ>ρcrust = ρ0/10. Below ρcrust we use the EoS of the R [km] crust obtained in [17, 18]. FIG. 4: Neutron-star masses as a function of the radius R. Solid and dashed curves are for MP1a and MP2a, respectively. TABLE I: Fitting parameters for the neutron matter EoS, a Dotted curve is for MP0. and b being in MeV. a α bβ MP0 10.8 0.442 5.31 1.88

MP1a 11.4 0.481 7.70 2.41 3.0 MP2a 6.79 0.192 11.6 1.93 2.5

The EoS’s for MP1a and MP2a violate causality and 2.0 predict sound speeds over the speed of light above a crit- ical density. Then, we adopt the approximation where solar 1.5 the EoS is replaced by the causal EoS above this density M/M 1.0 in the same way as the treatment in [16]. In Fig.4 and Fig.5, the calculated star masses are 0.5 given as a function of the radius R and the central den- 0.0 sity ρc, respectively. Solid and dashed curves are for 0.0 0.5 1.0 1.5 -3 MP1a and MP2a, respectively. Dotted curves are for c [fm ] MP0. The EoS’s in the cases of including MPP contributions are found to be stiff enough to give maximum FIG. 5: Neutron-star masses as a function of the central den- masses larger than the observed mass Mobs =1.97Msolar sity ρc. Also see the caption of Fig.4. for J1614-2230. The maximum masses for MP1a are substantially larger than that for MP2a. The difference is due to the quartic-pomeron exchange term included in the former. The strengths of the effective two-body interaction derived from triple- and quartic-pomeron exchanges are proportional to matter density ρ and ρ2, re- be determined by fitting the observed angular distribu- spectively. Then, the contribution of the latter become tion of 16O+16O elastic scattering at E/A = 70 MeV sizeable in the high-density region, which makes the max- with use of the G-matrix folding potential derived from imum mass large. However, the difference between MP1a ESC+MPP+TNA and adjusting to reproduce the min- and MP2a contributions cannot be seen in Fig.1, which imum value 16 MeV of the energy per nucleon at means that the quartic-pomeron exchange contribution normal density∼ − in symmetric nuclear matter. Then, demonstrated in neutron-star masses cannot be found by the empirical values of K, Esym and L are reproduced analyses of nucleus-nucleus scattering data. reasonably. The EoS of neutron-star matter obtained from ESC+MPP+TNA is stiff enough to give the large neutron-star mass over 1.97Msolar. It should be noted VI. CONCLUSION that our stiff EoS is determined with use of the experimental data for nucleus-nucleus scattering and nuclear In order to explain the observed maximum mass saturation properties. Our MPP contributions exist uni- of neutron stars, three-body repulsions are considered versally in every baryonic system. It is our future subject to work universally among three baryons. The multi- to show that our stiff EoS with universal MPP repulsions pomeron potential (MPP) is introduced as a model for is free from the softening effect induced by hyperon mix- such a universal three-body repulsion. Furthermore, the ing to neutron-star matter. Such a work will be comple- three-nucleon attraction (TNA) is added phenomenolog- mentary to the analyses based on the relativistic mean ically so as to reproduce the nuclear saturation prop- field models [19–21] for a massive neutron star including erty precisely. The strengths of MPP and TNA can hyprons. 5

Acknowledgments

We wish to thank Prof. T. Uesaka for giving us the motivation of this work.

[1] P.B. Cemorest, T. Pennucci, S.M. Ransom, M.S.E. [11] F. Nuoffer, et al., Nuovo Cimento A111, 971 (1998). Roberts, and J.W. Hessels, Nature (London) 467, 1081 [12] I.E. Lagaris and V.R. Pandharipande, Nucl. Phys. A359, (2010). 349 (1981). [2] S. Nishizaki, Y. Yamamoto, and T. Takatsuka, Prog. [13] G. Colò, N. Van Giai, J. Meyer, K. Bennaceur, and P. Theor. Phys. 105, 607 (2001); 108, 703 (2002). Bonche, Phys. Rev. C70, 024307 (2004). [3] Th.A. Rijken, M.M. Nagels, and Y. Yamamoto, Prog. [14] Li-Gang Cao, H. Sagawa, and G. Colo, Phys. ReV. C86, Theor. Phys. Suppl. 185, 14 (2010). 054313 (2012). [4] Th.A. Rijken, M.M. Nagels, and Y. Yamamoto, in Pro- [15] P. Möller, W.D. Myers, H. Sagawa and S. Yoshida, Phys. ceedings of the International Workshop on Strangeness Rev. Lett. 108, 052501 (2012). Nuclear Physics, Neyagawa 2012, Genshikaku Kenkyu [16] S. Gandolfi, J. Carlson, and Sanjjay Reddy, Phys. Rev. 57, Suppl.3, 6 (2013). C85, 032801(R) (2012). [5] M. Baldo, A. Fiasconaro, H.Q. Song, G. Giansiracusa, [17] G. Baym, A. Bethe, and C. Pethick, Nucl. Phys. A175, and U. Lombardo, Phys. Rev. C65, 017303 (2001). 225 (1971). [6] T. Furumoto, Y. Sakuragi, and Y. Yamamoto, Phys. Rev. [18] G. Baym, C.J. Pethick, and P. Sutherland, Astrophys. C79, 011601(R) (2009); C80, 044614 (2009). J.170, 299 (1971). [7] D.T. Khoa, et al., Phys. Rev. Lett. 74, 34 (1995). [19] S. Weissenborn, D. Chatterjee, and J. Schaffner-Bielich, [8] D.T. Khoa, G.R. Satchler, and W.von Oertzen, Phys. Nucl. Phys. A881, 62 (2012). Rev. C56, 954 (1997). [20] I. Bednarek, P. Haensel, J.L. Zdunik, M. Bejger, and R. [9] A.B. Kaidalov and K.A. Ter-Materosyan, Nucl. Phys. 74, Manka, Astronomy & Astrophysics A157, 543 (2012). 471 (1974). [21] Wei-Zhou Jiang, Bao-An Li, and Lie-Wen Chen, Astro- [10] J.B. Bronzan and R.L. Sugar, Phys. Revs. D16, 466 phys. J.756, 56 (2012). (1977).