# Full text of "Eta nuclear bound states revisited"

## See other formats

T] nuclear bound states revisited E. Friedman^, A. GaP'*, J. Mares^ ""Racah Institute of Physics, The Hebrew University, 91904 Jerusalem, Israel fT^ I ^Nuclear Physics Institute, 25068 Rez, Czech Republic o (N ^! Abstract < _i. ' The strong energy dependence of the s-wave r]N scattering amphtude at CN ■ and below threshold, as evident in coupled-channels i^-matrix fits and chiral models that incorporate the Su A^*(1535) resonance, is included self consis- tently in 77-nuclear bound state calculations. This approach, applied recently in calculations of kaonic atoms and ^-nuclear bound states, is found to im- pose stronger constraints than ever on the onset of r^-nuclear binding, with a pj ■ minimum value of Re Qr^N ~ 0.9 fm required to accommodate an r^-'^He bound state. Binding energies and widths of ?7-nuclear states are calculated within several underlying r]N models for nuclei across the periodic table, including ^^Mg for which some evidence was proposed in a recent COSY experiment. > 00 l/^ \ Keywords: ^ ' pion-baryon interactions, mesons in nuclear matter, mesic nuclei ^' O cn g '■ PACS: 13.75. Gx, 21.65.Jk, 21.85.+d 1. Introduction /\ . Searches for meson-nuclear bound states have focused on K~ and r] c^ \ mesons, motivated by a general theoretical consensus that the near-threshold ^A^ and r]N attraction generated by the s-wave resonances A(1405) and A^*(1535), respectively, translates into sufficiently attractive A^'-nucleus and ?7-nucleus interactions. A corollary of this resonance dominance is a strong energy dependence of the underlying ^A^ and r]N interactions. Here we apply the lessons gained by handling the strong energy dependence of the near-threshold KN interaction in A'~-nuclear calculations [l| to ry-nuclear * Corresponding author: Avraham Gal, avragal@vnis.huji.ac.il Preprint submitted to Physics Letters B April 25, 2013 bound-state calculations. Early calculations by Haider and Liu |2|, y] pre- dicted r^-nuclear bound states beginning with nuclear mass number A ~ 12. In these, as well as in a follow-up calculation J4|, a fairly weak rjN attrac- tion input was used, with Re a^Tv ^ 0.3 fm, where ar^N is the rjN scattering length. Several versions of coupled-channels chiral models give similar val- ues as well [a, la, LZI] whereas other models, particularly those using i^-matrix methods to fit ttN and 7A^ reaction data in the A^*(1535) resonance re- gion, e.g. [sl, y, |lO|, yield considerably stronger rjN attraction with values of Re ttriN ~ 1 fmo This might suggest that the onset of ?7-nuclear binding occurs already in the He isotopes for which strong final-state interaction pre- cursors have been noted in deuteron-initiated rj production [ll|, ll2|. A robust pattern of ry-nuclear bound states could yield useful information on the size of SU(3) flavor 77 — 77' mixing and about axial U(l) dynamics 13 1300 1400 1500 1600 1700 1800 E (MeV) cm ^ Figure 1: Energy dependence of the s-wave rjN cm scattering amplitude friN{\/s), with scattering length o^at = frjN{^/sth) = 0.91(6) + i0.27(2) fm, as calculated by Green and Wycech (GW) fitting pion reaction and photoproduction data over a wide range of cm energies Ecm (=v^)- The spike of Re / is located at the t]N threshold, whereas Im / peaks close to the A^*(1535) mass. Figure adapted from Ref. [8|. Regardless of the strong model dependence of Re a^^r, all studies of the TjN system near threshold, y^Sth = ^^n + ^^77 ~ 1487 MeV, agree that both real and imaginary parts of the s-wave center-of-mass (cm) scattering am- ^As for Im arjN, it is constrained by ttN — > rjN cross-section measurements, with values of Im a^jv ^ 0.2-0.3 fm in most theoretical analyses. Table 1: rjN scattering length a^^v (in fm) in three coupled-channels models used in the present work. Ml and M2 correspond to versions I and II, respectively, of rjN amplitudes from the recent chiral-model work by Mai et al. [7[, and GW denotes the if-matrix rjN amplitude due to Green and Wycech ^ shown in Fig. [T] Model Re CL-ql^l Im ariN Ml 0.22 0.24 M2 0.38 0.20 GW 0.96 0.26 plitude fj^N decrease steeply in going subthreshold, as illustrated in Fig. [H Since the in-medium rjN interaction relevant to the evaluation of ?7-nuclear bound states involves subthreshold rjN configurations, a procedure for go- ing subthreshold is mandatory. Previous calculations focused on shifting the energy variable of /,,Af(^/i) or its in-medium version by a fixed amount be- low threshold: 5y/s = —30 MeV was found in Ref. |j| to provide a good approximation to a variety of off-shell effects, whereas 6^/s = —B^i, with B^i the r^-nuclear binding energy, was used in Ref. IJ]. The latter proce- dure requires a self-consistent calculation to ensure that the S^ argument of the input /,,Ar coincides with the -B^ output of the binding energy calcula- tion. However, it was shown in our recent studies of ii"~-nuclear dynamics 15L lid . llTl . llSL Il9l . |20| that a more involved self-consistent calculation is re- quired to correctly implement the subthreshold energy dependence, and it is this self-consistent procedure that is applied here to calculate ?7-nuclear bound states. This procedure results in imposing stronger constraints than ever on the onset of 77-nuclear binding. Below we proceed to describe briefly the self-consistent procedure used to handle the subthreshold energy dependence of the rjN amplitude for bound nucleons, and its embeding into a dynamical Relativistic Mean Field (RMF) scheme which allows for the flrst time to consider the polarization of the core nucleus by the bound t] meson. To span a broad range of bound-state scenarios we apply our methodology to three distinct rjN amplitude models, with threshold values listed in Table [1] These amplitudes differ primarily in the value of the real part, while their shape below threshold exhibits a substantial decrease particularly for Im a^Ar, as illustrated for the GW amplitude model |8[ in Fig. [H We have calculated 77-nuclear bound states across the periodic table for these three amphtude models, as reported and discussed here for Is^ states. Finally, we also confront our results with a recent experimental suggestion of a ^Mg bound state [21 . 2. Methodology In close analogy to the latest calculation of K~ nuclear bound states [18 we calculate t] nuclear bound states by solving self consistently the Klein- Gordon (KG) equation [ v^ + ul -ml- n,K, p) ] v^ = 0, (1) where u^^ = Un — ir^/2 and w^ = m^ — B^, with B^ and F^ the binding energy and the width of the r] nuclear bound state. The self-energy operator n^ is related to a density- and energy-dependent optical potential K7 which is given by the following "tp" form: 11^(0;^, p) = 2ujnVr, = -AnFriNiVs, p)p, (2) where s = {E^ + E]^Y — (j)^ + p^Y is the Lorentz invariant Mandelstam variable s which reduces to the square of the total r]N energy in the two- body cm frame and F^at is the in-medium rjN s-wave scattering amplitude in the lab system. Note that for A^ 1 the lab system approximates well the ?7-nucleus cm system. Our in-medium F^jv accounts for Pauli correlations in the Ericson-Ericson multiple-scattering approach, as reformulated in Ref. [22 and used recently in Ref. [20 FM^s,p) = ^-^^^^ , e(p) = ^, (3) i + e(p)V(v^)p ^pV ^ where fr,N{,y/s) = {y/s/mN)frjN{y/s), with the kinematical factor y/s/mj^ transforming / from the two-body cm frame to the lab /, and where pp is the local Fermi momentum corresponding to density p = 2p|./(37r^). Note that Fj^N{y/s,p) — i- fnN{y/s) upon p — )• 0, as required by the low-density limit. Extensions of Eq. ((3)) to coupled channels and inclusion of self energies are discussed elsewhere [23.] . The in-medium y/s = -^/sth + S^/s energy argument of F^at in Eq. ([3]) is density dependent, with a form used recently in K -atom studies [20 6,/^ ^-Bn-- ^nB^^ - ^NTN{-f'^ + ^r^Re K,(5^, p), (4) P Po Po where ^Ar(r;) = f^N {■q) I ij^N + f^ri) , B]\r ^ 8.5 MeV is the average nucleon bind- ing energy, T^r ^ 23 MeV is the average bound-nucleon kinetic energy in the Fermi Gas Model, and po (p) is the maximal (average) nuclear density. The Vrj term in expression (jl]) originates from the nonzero in-medium momentum (squared) of the r] meson in finite nuclei, and is related in the local density approximation to the r^-meson kinetic energy [15]. This contrasts with the common assumption Pr^ = made in nuclear matter calculations. The depen- dence of Vrj on energy through 5^ and on density p is explicitly marked in this expression. Eq. (jl]) satisfies the low-density limit requirement that 6^/s — )■ with p — )■ 0. Note that for an attractive K; and as long as p 7^ 0, the shift of the two-body energy away from threshold is negative definite, 6y/s < 0, even as i?^ — )■ 0. For a given 5^, neither a/s nor Vr^ can be evaluated sepa- rately, implying that V^ is to be constructed self consistently together with y/s, which takes typically about 5 cycles of iteration. Once V^(i?j,,p) has been determined, it is used in the KG Eq. ([1]) to solve for the binding en- ergy eigenvalue Bi^ in the 77-nuclear n(. single-particle state. While varying the value of B^ in this process, the self-consistent requirement Eq. ^ is imposed at each step of the calculation of the eigenvalue. > u i 1 1 1 1 1 1 1 1 1 \ s . 10 — \^^-^ - 20 — \ ^^ '■. _ \ ""-■;••. Ml ^\ ^^' ' • . 30 — \. ^^^ '■-. — \ M2 ""-^■•••: 40 GW ^v 50 - \. - 1,1,1,1, 0.2 0.4 0.6 P/Po 0.8 Figure 2: Subthreshold rjN energies probed by the rj nuclear potential as a function of the relative nuclear RMF density in Ca. Each of the three curves was calculated self- consistently according to Eq. ^ for a specific version of rjN subthreshold amplitude model, see text. In Fig. [2] we show the downward subthreshold energy shift S^/s = E — Et^ as a function of the nuclear RMF density p in Ca, calculated self consistently according to Eq. (jl]) for riN amplitude models Ml, M2 and GW (see Ta- ble [1]). The hierarchy of the three curves reflects the strength of the input Re fr)N{.\fs) in the subthreshold region, with threshold values listed in Ta- ble dJ It is clear that downward energy shifts of up to ~55 MeV are involved in the present self-consistent calculations. 3. Results and discussion The methodology described in the last section was used to solve the KG equation ([1]) for r^-nuclear bound states across the periodic table. In this Letter we highlight the systematics of the Is^ bound state and compare our treatment of subthreshold energy dependence with previous studies. A more detailed account plus extensions are given elsewhere 23j. Three represen- tative rjN amplitude models Ml, M2 and GW (see Table [1]) are employed here in order to span a wide range of tjN interaction strengths. Our main results are shown in Fig. [3] for binding energies i?^ and widths F^^ calculated for Is^ nuclear states in core nuclei from ^^C to ^°'^Pb. RMF equations of motion, along with the KG equation ([T]), are solved self consistently [24 thereby allowing for core polarization by the rj meson (see Ref. [16| for the latest application to K~ mesons). The core polarization effect on 5^ and F^ was found in all cases displayed here to be less than 1 MeV. Therefore, the use of static nuclear densities is acceptable for not-too-light nuclear cores. Inspection of the l.h.s. of Fig. |3] reveals that for each of the three input rjN amplitude models the binding energy increases with A and appears to saturate for large values of A. As in Fig. [2l here too the hierarchy of the three curves reflects the strength of the Re frjNiVs) input in the subthreshold region, with threshold values listed in Table [H The Ml and M2 amplitudes are too weak to produce a Is,, bound state in ^^C, with the onset of binding for the weaker Ml amplitude deferred to around ^°Ca. Of our three represen- tative amplitudes. Ml is the closest one on shell to the Haider-Liu standard amplitude [2|, y, \^ which was used by these authors to argue for ^^C as the approximate onset of r] nuclear binding. In contrast. Re fr,N{,y/s) of the GW model is sufficiently strong to bind the Is^ state in ^^C and in lighter core nuclei, in spite of the suppression it undergoes here by forming its in-medium version and dealing with its energy dependence. The GW amplitude model 75 - Ca Zr 1 Pb ..'^ GW 20 > S 15 m C Mg.~ - 10 .7— ■■ M2 - ,..^- 5 -^.- Ml 1 - 10 20 30 40 Figure 3: Binding energies (left) and widths (right) of Is,, nuclear states across the periodic table calculated self consistently using the Ml, M2 and GW subthreshold ijN scattering amplitudes within a dynamical RMF scheme, see text. even admits a Isrj bound state in ^He with as low a binding energy as 1.2 MeV and a width of 2.3 MeV, both calculated using a static ^He density. Inspection of the r.h.s. of Fig. [3] reveals a trend for the three curves of calculated widths which is opposite to that observed on the l.h.s. for the calculated binding energies. Here, the GW model produces relatively small widths of order 4 MeV uniformly across the periodic table, whereas Ml and M2 give larger widths, particularly Ml with widths of order 20 MeV. This reflects partly the energy dependence of Im friN{y/s) in the subthreshold region, which is quite distinct in each one of the three amplitude models, and partly the difference in the in-medium renormalization arising from the Re friN{\/s) input. The latter point is readily understood by noting in Fig. [2] that the largest values of subthreshold downward energy shift are due to the GW subthreshold amplitude. This causes a particularly large reduction in the strength of the Im friN{y/s) input for the GW amplitude model. Focusing on a given core nucleus, we show in Table |2] results of static- density calculations in models Ml, M2 and GW of the Is,, state in ^^Mg with (Yes) and without (No) employing the in-medium modification of Eq. (|3]). The self-consistency requirement imposed by Eq. (jlj) is used and comparison is made with two other procedures applied in previous studies to incorporate energy dependence. These are (i) a fixed 30 MeV downward shift applied to the free-space r]N amplitude by Haider and Liu (HL) [4]; and (ii) shifting Table 2: Static calculations of Is,, binding energy (B) and width (F) in ^^Mg, using three rjN amplitude models (Ml, M2, GW) with (Yes) and without (No) medium corrections from Eq. ([3]), for several procedures of treating the energy dependence of fr^N- Energies and widths are given in MeV. Eq. §i> subthreshold Ml M2 GW B, r,, B^ r. B, r. No 5^ = -30 - - 3.0 11.2 31.2 10.0 No (5v^ Eq. (iD - - 3.2 10.6 23.8 7.4 Yes 6^/s = -B^ 3.8 21.7 8.3 13.2 19.4 5.8 Yes 5^s Eq. (H - - 2.5 7.4 14.8 3.9 down self-consistently the energy argument of the in-medium f^^^ arnplitude by the resultant i?^, as implemented by Garcia- Redo et al. (GR) [ij]. In the first two lines of the table we compare the HL procedure with ours, using free-space amplitudes. Both do not produce binding for the weakest amplitude Ml and practically agree for M2, while disagreeing significantly for the strongest GW amplitude. In the last two lines we compare the GR procedure with ours, using in-medium amplitudes. The GR procedure is found to give significantly higher binding energies and widths than ours for all amplitude models tested here, particularly for the weaker M amplitudes where it is the only one that produces a Is,, bound state for Ml. Of the three models used by us with in-medium amplitudes in Table [2] (last line), only GW provides B^, which is comparable with S^"P(2jMg) = 13.1 ± 1.6 MeV, r"P(2jMg) = 10.2 ± 3.0 MeV, (5) deduced from the following ^^Mg interpretation of a peak observed by the COSY-GEM Collaboration ^■. p + 27 Al ^ 25^g ^ 3gg ^ (tt- + p) + X, (6) with a decay induced hj r] + n ^ n^ + p. Hence, if this peak assignment to a Isri state is correctp then the underlying threshold value Re ajjN must be ^This has been contested recently by Haider and Liu who offered a different interpre- tation of the observed peak [25| . 8 rather large, close to 1 fm. Other procedures listed in Table [2] for treating the subthreshold r]N energy dependence require considerably smaller values of Re a^Tv- Finally, the relatively small value of width F produced in the GW model should not be viewed as too restrictive since the total width must be larger than given in these models, owing to true rjNN absorption and two-pion production rjN — )■ -r-kN processes that are not accounted for by the models considered in the present work. 4. Conclusions In this work we have demonstrated the importance of, as well as the sub- tleties involved in constructing self-consistent 77-nucleus optical potentials that incorporate the strong subthreshold energy dependence of the under- lying rjN scattering amplitude. Of the three rjN amplitude models studied here self-consistently, even the relatively weak attraction in model Ml with a threshold value Re a^A? ~ 0.2 fm requires going down to about 35 MeV below threshold, as shown in Fig. [21 in order to calculate reliably the ry-nuclear op- tical potential Ky(p) at central nuclear densities. This downward energy shift exceeds by far the downward shifts —B^, with i?^ < 20 MeV encountered in the self- consistent calculations of Garcia- Redo et al. [ij]. The relatively large downward energy shifts in the present approach together with the rapid decrease of the free-space and in-medium rjN amplitudes lead to smaller than ever binding energies calculated in comparable models ^ [ij]. Thus, ^^C bound states are unlikely in models with threshold values Re a^^v ^ 0.5 fm, and as large a value as Re ar^N ~ 0.9 fm is required to reproduce the ^^Mg bound-state hint from the recent COSY-GEM experiment [21 A value of Re a^^v ~ 0.9 fm is likely to yield a near-threshold ^He bound state, as found by us using the GW amplitude model, but it is short of binding ^He. Stretching the limits of optical potential usage down to these light systems is of course questionable, and corresponding few-body calculations are highly needed to resolve such issues. Nevertheless, if one applies our subthreshold self-consistency scheme to r^-^He low-energy scattering, then a very large imaginary part that might indicate a nearby virtual state is found for the ?7-^He cm scattering amplitude when using the GW amplitude model. This large imaginary part might be associated with the strong final- state interaction effects observed for the ry-^He system ll|. In contrast to previous estimates that assigned a value of Re ar/N ~ 0.5 fm to describe such occurrence [26|, in our self-consistent calculations it requires substantially larger values, more likely around 0.9 fm. Acknowledgements We thank Maxim Mai for providing us with numerical tables of r]N scat- tering amplitudes from Ref. |7|. A.G. thanks Colin Wilkin for useful and stimulating correspondence. J.M. acknowledges financial support within the agreement on scientific collaboration between the Academy of Sciences of the Czech Republic and the Israel Academy of Sciences and Humanities. This work was supported by the GACR Grant No. 203/12/2126, as well as by the EU initiative FP7, HadronPhysicsS, under the SPHERE and LEANNIS cooperation programs. References [1] For a recent review at HYP2012, see A. Gal, Nucl. Phys. A (2013) http://dx.doi.Org/10.1016/j.nuclphysa.2013.0L01l] flarXiv:1301.21 45). [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] Q. Haider, L.C. Liu, Phys. Lett. B 172 (1986) 257. L.C. Liu, Q. Haider, Phys. Rev. C 34 (1986) 1845. Q. Haider, L.C. Liu, Phys. Rev. C 66 (2002) 045208. T. Waas, W. Weise, Nucl. Phys. A 625 (1997) 287. T. Inoue, E. Oset, M.J. Vicente Vacas, Phys. Rev. C 65 (2002) 035204. M. Mai, P.C. Bruns, U.-G. MeiBner, Phys. Rev. D 86 (2012) 094033. A.M. Green, S. Wycech, Phys. Rev. C 71 (2005) 014001, and references listed therein to earlier work by these authors. R.A. Arndt, et al., Phys. Rev. C 72 (2005) 045202; see Table I for r]N scattering-length overview. V. Shklyar, H. Lenske, U. Mosel, Phys. Rev. C 87 (2013) 015201. For r/^He: C. Wilkin, et al., Phys. Lett. B 654 (2007) 92; T. Mersmann, et al. (COSY-ANKE Collab.), Phys. Rev. Lett. 98 (2007) 242301; see also F. Pheron, et al. (MAMI), Phys. Lett. B 709 (2012) 21. 10 [12] For r^^He: A. Budzanowski, et al. (COSY-GEM Collab.), Nucl. Phys. A 821 (2009) 193; P. Adlarson, et al. (WASA@COSY Collab.), Phys. Rev. C 87 (2013) 035204. [13' [14 [15 [16 [ir; [18 [19 [20 [21 [22 [23 [24 [25 [26 S.D. Bass, A.W. Thomas, Phys. Lett. B 634 (2006) 368. C. Garcia-Recio, T. Inoue, J. Nieves, E. Oset, Phys. Lett. B 550 (2002) 47. A. Cieply, E. Friedman, A. Gal, D. Gazda, J. Mares, Phys. Lett. B 702 (2011) 402. A. Cieply, E. Friedman, A. Gal, D. Gazda, J. Mares, Phys. Rev. C 84 (2011) 045206. E. Friedman, A. Gal, Nucl. Phys. A 881 (2012) 150. D. Gazda, J. Mares, Nucl. Phys. A 881 (2012) 159. N. Barnea, A. Gal, E.Z. Liverts, Phys. Lett. B 712 (2012) 132. E. Friedman, A. Gal, Nucl. Phys. A 899 (2013) 60. A. Budzanowski, et al. (COSY-GEM Collab.), Phys. Rev. C 79 (2009) 012201(R). T. Waas, M. Rho, W. Weise, Nucl. Phys. A 617 (1997) 449. A. Cieply, E. Friedman, A. Gal, J. Mares, in preparation. J. Mares, E. Friedman, A. Gal, Nucl. Phys. A 770 (2006) 84. Q. Haider, L.C. Liu, J. Phys. G: Nucl. Part. Phys. 37 (2010) 125104. N. Wilhs, et al., Phys. Lett. B 406 (1997) 14. 11