arxiv: v1 [nucl-th] 14 Sep 2009

Transcription

1 Resonant formation ofλ(145) by stopped-k absorption in d Jafar Esmaili,a,b, Yoshinori Akaishi a,c, Toshimitsu Yamazaki a,d a RIKEN, Nishina Center, Wako, Saitama , Japan b Department of Physics, Isfahan University of Technology, Isfahan , Iran c College of Science and Technology, Nihon University, Funabashi, Chiba , Japan d Department of Physics, University of Tokyo, Bunkyo-ku, Tokyo , Japan arxiv:99.573v1 [nucl-th] 14 Sep 9 Abstract We show that the current debate on the position of the quasi-bound K p state, namely, Λ(145) orλ (14), can be solved by measuring experimental observables that involve the transition matrix element T 1 = T Σπ KN. As one of such ways, we propose to observe theσπ invariant-mass spectrum in stopped-k absorption in deuterons, since the spectrum, reflecting the deuteron momentum distribution, is expected to have a narrow direct-capture component with an upper edge of M= 143 MeV/c, followed by a significant high-momentum tail toward the lower mass region, where a resonant formation of Λ(145) is revealed as an isolated peak. 1. Introduction Where is the position of the I= L= K p quasi-bound state? This question is directly connected to the strength of the s-wave I = KN interaction. Traditionally, the Λ(145) resonance is identified to this state, and a strongly attractive KN interaction is indicated, which is compatible with theoretical predictions based on meson-exchange [1] and on chiral dynamics []. The attractive KN interaction was evidenced by an anomalously large subthreshold production of K in heavy-ion reactions [3]. Starting from this Λ(145) Ansatz, we constructed a phenomenological KN potential by a coupled-channel procedure, and predicted deeply bound dense kaonic nuclear systems [4, 5, 6, 7, 8]. Recently, on the other hand, another theoretical framework including a double-pole hypothesis has been proposed, claiming that the K p bound state is located around 14 MeV [9, 1, 11, 1, 13, 14]. We call such a hypothetical state Λ (14). The predicted shallow resonance regime (we call this chiral weak, whereas we call the earlier chiral prediction [] chiral strong ) leads to a much less-attractive KN interaction, and thus, no deeply bound states are expected. Since these two alternative regimes are pertinent to the problem of the dense kaonic system caused by a super strong nuclear force [15, 16] (ultimately toward kaon condensation), it is vitally important to distinguish betweenλ(145) andλ (14) experimentally, but there seem to be lots of confusing statements concerning the strategy as to how valid experimental evidence can be obtained. To clarify the situation let us take a simple and sufficient coupled-channel regime with KN (=1) and Σπ (=) channels. The imaginary part of the transition matrix T 11 = T KN KN, which corresponds to the imaginary part of the forward scattering amplitude, has a peak at mass M (either 145 or 14 MeV/c ), but, since it is located below the K+N binding thresh- Corresponding author; jesmaili@riken.jp old, it cannot be detected directly. On the other hand, the transition matrix T 1 = T Σπ KN is responsible for the observation ofσπ pairs following K capture (or equivalently, following associated production together with K + ). The chiral-weak theories predict two poles. The 1st pole, mainly coupled to KN, corresponds to a peak around 14 MeV/c, as shown in Fig. 4 of [14] (The 1st-pole position is even higher, around 143 MeV, but its peak position is pushed down by the effect of the KN threshold). The nd pole, mainly coupled toσπ, is broadly distributed around 14 MeV, which can appear mainly in T = T Σπ Σπ. Contrary to the prevailing belief, the nd pole cannot be detected directly. It will be shown later that any M Σπ spectrum involving T 1 is dominated by the pole of K p; it is not much affected by the presence of a second pole in theσπ channel. Thus, the issue ofλ(145) vsλ (14) is a well-defined problem, and can be solved by experimental observables involving T 1. In our recent work [17] we have shown that theσπ invariantmass spectra (M Σπ ) in stopped-k absorption in 4 He, 3 He and d are governed by the spectator momentum distributions of t, d and n, respectively, but do reflect the resonant formation of a quasi-bound K p state, contrary to the past interpretation in terms of the direct-capture process without resonance. Thus, the issue of the location of the K p resonance state can be examined experimentally by a quantitative comparison of an observed M Σπ spectrum with predicted theoretical distributions including resonant formation. We made aχ analysis of old bubble-chamber data of M Σπ ( 4 He) [18] by varying the mass (M) and width (Γ) of an assumed resonance and found a significant minimum in the M Γcontour presentation ofχ at M= MeV/c andγ= MeV. (1.1) Thus, theλ (14) Ansatz is excluded by more than 99.99% confidence. However, the Λ(145) signal does not appear as a separate peak, but as a small component of the steeply falling Preprint submitted to Elsevier September 14, 9

2 tail, and this discrimination seems to be delicate. One would hope to have a more clear-cut case. In the present paper, we point out that the use of a deuteron target in the reaction, stopped K + d X+ n, (1.) X Σ+π, (1.3) can provide a more decisive conclusion. Since the deuteron wavefunction is composed of low- and high-momentum components, the dominant quasi-free shape of M Σπ is narrow, whereas its tail, resulting from the high-momentum component of d, extends to the region where a resonant formation of Λ(145) (but not ofλ (14)) can be revealed as a separate peak. We investigate this problem in detail. arbitrary scale - M inv (Σπ<-KN) AY HW AY HW HW -(1/π) ImT11 HW AY HW AY C T 1 C T - M miss (KN) k k Number of events Width [MeV] CL 68% 9% 99% M Σπ [MeV/c ] Zychor data Gaussian Chi^= Chi^= Chi^= Chi^=16. Events per bin [MeV] M inv (Σπ<-Σπ) Escape Conversion s 1/ [MeV] DD91 (K-Matrix) T Hyodo-Weise T T T Total cm energy [MeV] Hemingway's data Zychor's data Figure 1: (Upper) ( 1/π) Im T 11, T 1 k and T k curves in the chiral-weak model of Hyodo-Weise [14], in comparison with the Λ(145) Ansatz of A-Y [4]. (Lower) The experimental M Σπ spectrum of Hemingway [19], when fitted by T of the phenomenological model of DD91 [1] and the chiral-weak model of HY. Zychor s data [3] are also shown.. Coupled-channel analysis of Λ(145) data.1. Observables of KN-Σπ coupled channels The transition-matrix elements of a coupled KN(= 1) Σπ(= ) system are T 11, T 1, T 1 and T. Among them the experimentally observable quantities below the K + N threshold are M Λ* [MeV/c ] Figure : (Upper) The COSY-ANKE spectrum of M(Σ π ) compared with the best-fit theoretical curve as well as in the two typical cases. (Lower) Contour mapping of likelihood in the plane of M vs Γ. 1/π Im T 11, T 1 k and T k, where k is theπσ relative momentum. The first one, corresponding to a KN missing-mass spectrum, is proportional to the imaginary part of the scattering amplitude, the peak position of which is just our concern. The second one is the Σπ invariant-mass spectrum from conversion process, KN Σπ (we call this the T 1 invariant mass ). The third one is a Σπ invariant-mass spectrum from a scattering process,σπ Σπ (we call this the T invariant mass ). It should be emphasized that the T 1 invariant-mass spectrum coincides with the KN missing-mass spectrum below the K+N threshold, as can be proved from an optical relation of Im T 11 = T 1 Im G in notations of Section 3. Thus, the observation of the T 1 invariant-mass spectrum is nothing but the observation of the imaginary part of the scattering amplitude given in Fig. 15 of Hyodo-Weise [14]. Figure 1 (Upper) shows three observable quantities calculated with Hyodo-Weise s chiral-weak two-channel model. The T 1 and the T invariant masses have peaks at different positions of 14 MeV and 145 MeV, respectively. In comparison, the figure also shows of the T 1 and the T invariant masses in theλ(145) Ansatz of AY [4]. It is to be noted that the issue here is to discriminate whether the peak position of the

3 T 1 spectrum is 14 MeV or 145 MeV. The stopped K + d spectrum of this paper is a T 1 invariant-mass spectrum projected with the deuteron momentum distribution. n-spectator process n-auger process HO Potential QF Argonne potential QF n p - K Λ* nn Σ π p n K - I = 1 Λ* π Σ n arbitrary scale arbitrary scale Neutron spectator process Nuclear Auger process 1 * (145-4) 1 * (14-4) M Σπ [MeV/c ] M Σπ [MeV/c ] Figure 4: M Σπ spectra of stopped K on d, adjusted so that the peak heights to be equal. (Left) For the HO case, and (Right) the case where the Argonne v18 SC potential is used as the interaction between p and n M Σπ [MeV/c ] Figure 3: (Upper) Diagrams for the neutron-spectator and the nuclear Auger processes. (Lower) Calculated M Σπ spectra in the two processes... Hemingway s data and the PDG value The second pole or resonance shape ofσπ could not be obtained, unless experimental observables involving the transition matrix T = T Σπ Σπ could be measured. The present-day PDG value of Λ(145) [] depends entirely on theoretical arguments made by Dalitz and Deloff (hereafter called DD91) [1]. They chose exclusively 1 data points below the KN threshold (namely, discarding the data above the threshold) among Hemingway sσ + π invariant mass spectrum [19], and searched for theχ minimum for T as a function of the resonance energy, E R, under a constraint of I = KN scattering length. DD91 expressed a strong preference for the M-matrix model, and recommended a value of (146.5± 4.) i(5±1) MeV, which is taken up as the PDG value. So far, the T analysis of the same experimental data of [19] has been done to materialize the double-pole hypothesis [9, 1, 11, 1, 13, 14]. However, the application of T to this spectrum is highly questionable, because these Σπ pairs are the decay products ofσ + (166), not those in the free scattering ofσand π. In fact, those who use the T formula (including DD91) are destined to predict a large dip at the KN threshold (M 143 MeV/c ), followed by a recovery above the threshold, which is a characteristic of T. This shape disagrees seriously with the observed spectrum (see Fig. 1). Detailed accounts will be given elsewhere []..3. ANKE data Recently, an invariant-mass spectrum ofσ π, a signature of the genuine I = Λ, was observed from the reaction p+ p K + + p+y at COSY-ANKE [3]. We made aχ fitting of the spectrum by theoretical ones with the conversion T 1 involving M andγas parameters. The best-fit theoretical curve as well as two typical curves (Λ and Λ14-4) were compared with the observed spectrum. All of these curves have skewed shapes, since we take into account theσπ emission threshold and the K + p threshold realistically [4]. The obtained likelihood contours are shown in Fig.. The best-fit values are M= MeV/c, Γ=4±8 MeV. (.4) The best-fit mass does not correspond to the apparent peak position of the theoretical shape, which is found to shift downward because the broad resonance is close to the KN threshold. Although the statistical errors are large, the data are in favour of theλ(145) rather than theλ (14) invoked by Geng-Oset [5]. 3. Coupled-channel formulation of the K p resonance As described in [17] in detail, we treat the K p quasi-bound state as a Feshbach resonance, considering two channels, KN andπσ. We employ a set of separable potentials with a Yukawatype form factor, k v i j k =g( k ) U i j g( k), (3.5) g( k) = Λ Λ + k, U i j= 1 π 1 µ i µ j Λ s i j. (3.6) Here i( j) stands for the KN channel, 1, or theπσ channel,, µ i (µ j ) is the reduced mass of channel i( j), and s i j are nondimensional strength parameters. 3

4 arbitrary scale A-Y single pole H-W two poles , A-Y single pole A-Y single pole M [MeV/c ] Σπ Figure 5: M Σπ spectra of stopped-k + d for differentλ(145) parameters including the Hyodo-Weise double-pole parameters. The spectra depend only on the 1st pole position. As before, we obtain s 11 and s 1 from the M andγvalues of an arbitrarily chosen K p state to be used to calculate theσπ invariant masses. We adopt s =.66, which gives U /U 11 = 4/3 forλ(145) as in a chiral model, andλ= 3.9/fm. The coupled-channel transition matrix, satisfies k t i j k =g( k ) T i j g( k), (3.7) T i j = U i j + U il G l T l j (3.8) with a loop function G l. The solution is given in a matrix form by T= [1 UG] 1 U. (3.9) In our treatment, the loop function is considered to be µ j Λ (UG) l j = s l j µ l (Λ i k j ), (3.1) where k j is the relative momentum in channel j. In the case of the reaction p+ p K + + p+y, the (Σπ) invariant mass distribution from Y decay is given by dγ d(m (Σπ) c ) g(k) T 1(E ) g(κ) k (3.11) 4 with E = µ k = M Σπ c M Σ c m π c, E 1 = E M= µ 1 κ, (3.1) where M = m K + M p m π M Σ is the threshold mass difference. 4.Σπ invariant mass from stopped K on d 4.1. The neutron-spectator and the nuclear Auger processes We showed [17] that the calculated M Σπ spectra in 4 He for the s-orbit capture are favoured. Now, the situation of a liquiddeuteron target is very similar to 4 He. When negative mesons or antiprotons are stopped in liquid-hydrogen targets, the exotic atoms formed in large-n atomic orbits behave like a neutral object and undergo violent collisions, which lead to s-orbit capture after Stark mixing decays [7]). Thus, we calculate M Σπ spectra from s-orbit absorption of K. The theoretical framework for the spectator model was given in detail in [17, 8]. The momentum distribution of the decay particles in the K absorption is given as: d Γ = (π)3 ψ atom dk Σ dk n c nlm () g(k ) T 1 (E ) g(k n /) k Σ k n E π F(k n ), (4.13) E = (E i E n ) c kn M Σc m π c, (4.14) whereψ atom nlm () is a K atomic wave function; k Σ, k π and k n are the momenta of theσ,πand the neutron, respectively, and the T 1 matrix involves theλ resonance effect. The kinematical constraints among the various momenta and invariant-mass spectra are given in [17]. With a similar approach, we have formulated the nuclear Auger process for K absorption at rest on d. In this process (Fig. 3), K is scattered first by one nucleon in the I= 1 channel, and subsequently resonant formation ofλ occurs in the I= channel. Here, the neutron is not a spectator but a participant of the reaction. The neutron is kicked out from the bound orbit to the continuum with a large momentum, like an Auger electron in an atomic system. We calculated the M Σπ spectra of the neutron spectator process from stopped-k absorption in d from the s-orbit for resonant capture to formλ(145) andλ (14) as well as of the nuclear Auger process. Figure 3 shows a comparison of the absolute value of the M Σπ spectra for the neutron-spectator and the nuclear Auger processes. In the neutron spectator process the M Σπ spectrum has a cusplike shape close to the threshold (M=143 MeV/c ), which results from a small momentum distribution of the spectator. The spectrum (4.13) is governed and projected by the spectator momentum distribution, F(k n ), because it is sharper than the resonance shape, as reflected in T 1 (E ). TheΛ (14) has an effect near the threshold. We notice that the Λ(145) resonance,

5 Number of events Number of events Γ [MeV] Theory : Simulation M Σπ [MeV/c ] Theory : 14-3 Simulation M Σπ [MeV/c ] a) b) c) 68% CL 9% CL 99% CL 99.9% CL Λ145-4 Λ M Λ* [MeV/c ] Λ14-4 Λ14-3 Figure 6: a) and b) Examples of simulation spectra assuming Λ(145) and Λ (14) withγ=3 MeV. The spectra above 14 MeV/c are unobserved because of the small efficiency in an actual experiment. c) Confidence contours obtained after aχ fitting of the four simulation spectra in the M-Γ plane. 5 which is far from the threshold, can be populated as a nearly isolated peak. On the other hand, in the nuclear Auger process, no such quasi-elastic peak exists, and only a resonantly formedλ peak shows up. However, the nuclear Auger process has a much smaller intensity than the neutron spectator process by orders of magnitude. It is negligible here, but becomes important in in-flight capture reactions. Thus, hereafter we adopt the neutron-spectator process as the main process to calculate the M Σπ spectra. 4.. M inv distributions Since d is a loosely bound system, which has a long tail in coordinate space, the low-momentum component of d becomes very important to the spectrum shape of M Σπ. Figure 4 shows a comparison between the harmonic oscillator (HO) and the Av18 SC potential [9] cases. In the HO case, the M Σπ spectra are broader and the contributions ofλ(145) andλ (14) are hardly distinguishable. On the other hand, in the realistic Argonne v18 SC case, the direct-capture shape arising from the low-momentum component is much narrower, whereas the high-momentum component produces a long tail in the spectrum on which resonant formation of Λ(145) appears. Thus, various models onλ can be distinguishable. Especially, in the case of our prediction forλ(m 146 MeV/c,Γ 6 MeV), the peak stands out as a large isolated bump. We calculated the M Σπ spectra for stopped K on d with s- orbit absorption for various models and parameters of Λ(145). Figure 5 shows the single-pole cases with parameters of and In addition, the calculated spectrum shape based on the double-pole model of Hyodo-Weise ( and ) is shown together with the single-pole case of The double-pole case is in excellent agreement with the corresponding single-pole case, clearly indicating that the second pole does not have any effect on M Σπ. We also examined this point using a model of Révai-Shevchenko [3], which provides both single-pole and double-pole cases. There is virtually no change of spectra, indicating that the second pole is irrelevant to experimental observables Effects of theσ(1385) resonance We considered the effect of the population of theσ (1385) resonance, which is known to decay toσ π + with a branching of about 5% []. It is shown that up to % mixing of theσ (1385) population forλ(145) andλ (14) withγ=4 MeV, the shapes of the spectrum are nearly unchanged. As expected, this smallσ (1385) mixing can influence on the shape of the spectrum just for lower mass of around 139 MeV Data simulation and fitting An experiment is planned at J-PARC [31]. The invariantmass spectrum can be obtained as a missing-mass spectrum of M(Kd, n) from observations of n together with Σπ events. The mass resolution depends on the momentum resolution of n. It is expected to be better than MeV, good enough for

6 the present purpose (reconstruction of M Σπ fromσandπmomenta gives a much worse resolution). However, the neutron detection efficiency drops down for the mass region above 14 MeV/c. This means that the quasi-free peak region is missing from observations. To examine whether we can obtain useful information on the mass and width ofλ from an experimental missing-mass spectrum, we produce simulation spectra assuming the cases ofλ(145) andλ (14), both withγ=3 MeV, as shown in Fig. 6 a, b), respectively. We attempt to fit the simulation spectra by a theoretical spectrum with assumed parameters (M,Γ) and drawχ confidence contours in Fig. 6 c). Although the simulation spectra lack the region above 14 MeV, the fitting results are remarkable. The data a) (Λ145-3) produce fitting results: M = ±.3 andγ=9.6±.8, which are in excellent agreement with the input data. On the other hand, the data b) (Λ14-3) give M= ± 1. MeV/c andγ=7.9±1. MeV, which deviate from the input data, particularly for Γ. This result can be understood in view of the missing region (M> 14 MeV/c ) of the simulation spectrum. Nevertheless, the present fitting procedure turns out to be capable of extracting the value of M to high precision. 5. Conclusion We have shown that the issue ofλ(145) vsλ (14) can be solved from experimental observables of conversion M Σπ, for which T 1 is responsible. Contrary to the prevailing belief, the effect of the presence of a second pole, if any, is negligible in such spectra. Since a deuteron has both low-momentum and high-momentum components rather separately, the M Σπ spectral shape from stopped-k absorption in d is characterized by a narrow direct-capture component, followed by a long tail, on which the resonant formation ofλ will be revealed, possibly as an isolated peak in the case of Λ(145). A proposed experiment at J-PARC will yield interesting results. [6] A. Doté, H. Horiuchi, Y. Akaishi and T. Yamazaki, Phys. Lett. B 59 (4) 51. [7] A. Doté, H. Horiuchi, Y. Akaishi and T. Yamazaki, Phys. Rev. C 7 (4) [8] T. Yamazaki, A. Doté and Y. Akaishi, Phys. Lett. B 587 (4) 167. [9] E. Oset, A. Ramos and C. Bennhold, Phys. Lett. B 57 () 99. [1] T. Hyodo, S.I. Nam, D. Jido and A. Hosaka, Phys. Rev. C 68 (3) 181. [11] B. Borasoy, R. Nißler and W. Weise, Eur. Phys. J. A 5 (5) 79. [1] B. Borasoy, U.-G. Meißner and R. Nißler, Phys. Rev. C 74 (6) 551. [13] D. Jido, J. A. Oller, E. Oset, A. Ramos and U.-G. Meißner, Nucl. Phys. A 75 (3) 181. [14] T. Hyodo and W. Weise, Phys. Rev. C 77 (8) 354. [15] T. Yamazaki and Y. Akaishi, Proc. Jpn. Acad. Ser. B 83 (7) 144. [16] T. Yamazaki and Y. Akaishi, Phys. Rev. C 76 (7) 451. [17] J. Esmaili,Y. Akaishi and T. Yamazaki, arxiv:96.55v1. [18] B. Riley, I-T. Wang, J.G. Fetkovich and J.M. McKenzie, Phys. Rev. D 11 (1975) 365. [19] R.J. Hemingway, Nucl. Phys. B 53 (1985) 74. [] Particle Data Group, W. M. Yao et al., J. Phys. G 33 (6) 1. [1] R.H. Dalitz and A. Deloff, J. Phys. G: Nucl. Part. Phys. 17 (1991) 89. [] Y. Akaishi et al., to be published. [3] I. Zychor et al., Phys. Lett. B 66 (8) 167. [4] Y. Akaishi, Khin Swe Myint and T. Yamazaki, Proc. Jpn. Acad. Ser. B 84 (8) 64. [5] L.S. Geng and E. Oset, Eur. Phys. J. A 34 (7) 45. [6] T. Yamazaki and Y. Akaishi, Phys. Lett. B 453 (1999) 1. [7] M. Leon and H.A. Bethe, Phys. Rev. 17 (196) 636. [8] T. Yamazaki and Y. Akaishi, Nucl. Phys. A 79 (7) 9. [9] R. B. Wiringa, V. G. J. Stoks and R. Schiavilla, Phys. Rev. C 51 (1995) 38. [3] J. Révai and N.V. Shevchenko, Phys. Rev. C 79 (9) 35. [31] T. Suzuki et al., Proposal for J-PARC 5 GeV Proton Synchrotron P3 (July, 9). 6. Acknowledgements The authors are grateful to Dr. T. Suzuki for his stimulating discussion. One of us (J.E.) would like to thank Professor E. Hiyama for inviting him to RIKEN and hospitality during the time that this work was performed, and thank Professor S.Z. Kalantari for helping him before coming to Japan. This work is supported by the Grant-in-Aid for Scientific Research of Monbu-Kagakusho of Japan. References [1] A. Müller-Groeling, K. Holinde and J. Speth, Nucl. Phys. A 513 (199) 557. [] T. Waas, N. Kaiser and W. Weise, Phys. Lett. B 365 (1996) 1; B 379 (1996) 34; N. Kaiser, P.B. Siegel and W. Weise, Nucl. Phys. A 594 (1996) 35. [3] A. Schröter et al., Z. Phys. A 35 (1994) 11; F. Laue et al., Phys. Rev. Lett. 8 (1999) 164. [4] Y. Akaishi, T. Yamazaki, Phys. Rev. C 65 () 445. [5] T. Yamazaki, Y. Akaishi, Phys. Lett. B 535 () 7. 6