Extraction of the K K isovector scattering length from pp dk + K0 data near threshold arxiv:hep-ph/0508118v1 10 Aug 2005 M. Büscher 1, A. Dzyuba 2, V. Kleber 1, S. Krewald 1, R. H. Lemmer 1, 3, and F. P. Sassen 1 1 Institut für Kernphysik, Forschungszentrum Jülich, D 52425 Jülich, Germany 2 High Energy Physics Department, Petersburg Nuclear Physics Institute, 188350, Gatchina, Russia 3 Max Planck Institut für Kernphysik, Postfach 10 39 80, D 69029, Heidelberg, Germany. (Dated: February 19, 2008) Abstract The results of a recent experiment measuring the reaction pp dk + K0 near threshold are interpreted in terms of a spectator model that encapsulates the main features of the observed K + K0 invariant mass distribution. A χ 2 fit to this data determines the imaginary part of the isovector scattering length in the K K channel as Im(a 1 ) = (0.61 ± 0.05) fm. We then use the Flatté representation of the scattering amplitude to infer a value Re(a 1 ) = (0.02 ± 0.03) fm for the real part under the assumption that scaling (Baru et al.[17]) is approximately satisfied. PACS numbers: 13.75.-n, 14.40.Cs, 25.40.Ve Corresponding author: R H Lemmer Electronic address: rh lemmer@mweb.co.za 1
1. Introduction The reaction pp dk + K0 has been measured at an excess energy of Q = 46 MeV with the spectrometer ANKE at the cooler synchrotron COSY-Jülich[1]. This experiment has achieved a resolution of 1 8 MeV for the invariant mass spectrum of the K K system which considerably improves the data base. Therefore it is timely to investigate what information regarding the K K scattering length can be extracted from this data. A model-independent analysis of the angular distributions shows a dominance of an s wave between the two kaons accompanied by a p wave deuteron with respect to the mesons. This has been interpreted as evidence for a dominant production of the a 0 (980) near threshold [1, 2]. Detailed calculations of the reaction pp dk + K0 performed by Grishina et al.[3] confirm this conclusion. In particular the non resonant K + K0 production is shown to be small below excess energies of Q = 100 MeV. The final state interaction in the Kd system reduces the cross section, but leaves the invariant mass distributions virtually unaffected [3]. A. Sibirtsev et al. [4] have extracted limits for the Kd scattering length from the data of Ref. [1] and found the magnitudes of both the real and imaginary parts of this quantity to be less than 1.3 fm. A chiral approach to the reaction pp dk + K0 has been formulated by Oset, Oller and Meissner [5] that suggests a potentially large sensitivity of the data to the Kd scattering length. A recent study [6] based on this formalism finds an optimal fit to the COSY data for a purely imaginary Kd scattering length of a Kd = (0.0 + i1.0) fm. The relatively small influence of the antikaon deuteron final state interaction on the shape of the measured invariant mass distributions as indicated in Refs. [3, 4, 6] encourages us to explore the working assumption that the deuteron s role is simply that of a spectator, and that the K + K0 production in pp dk + K0 is dominated by final state interactions in the meson meson channels only, and so can proceed via the resonant intermediate mode π + η a + 0 K + K0 with a cross section proportional to the I = 1, I 3 = 1 cross section σ a (π + η K + K0 ). This is our basic assumption. The approximation should be particularly appropriate in the vicinity of the K K threshold. Note further that this model cannot make any predictions about the three additional measured angular distributions or the d K 0 invariant mass distribution reported in Ref.[1]. However it does allow one to extract the imaginary part of the K K isovector scattering length a 1 as described below. 2. Spectator model calculation The low energy isospin one amplitude for K K scattering at low relative momentum k in the CM system has the form f(θ) = a 1 (1 ika 1 + ) 3c 1 k 2 cosθ + (1) where a 1 and c 1 are the (complex) s and p wave scattering length and volume respectively. Using the optical theorem, the s and p wave annihilation cross sections are obtained from and σ (s) a (K + K0 π + η) = σ (s) tot σ (s) elas 4π k Im(a 1)[1 + 2Im(a 1 )k] (2) σ (p) a (K+ K0 π + η) 12πkIm(c 1 ) (3) 2
in this approximation after dropping the O(k 4 ) elastic correction from the latter expression for consistency. It is clear from these expressions that the spectator model is limited to fixing only Im(a 1 ) and Im(c 1 ) in the effective range approximation to this order in k. Due to phase space restrictions the kaon pair invariant mass window that is sampled in the experiment is quite restrictive. Hence the k dependence featured in Eqs. (2) and (3) turns out to be quite adequate to describe the important features in the experimental data. Using two channel unitarity, one now has p 2 σ (l) a (π + η K + K0 ) = k 2 σ (l) a (K + K0 π + η), l = 0, 1 (4) where p is the CM momenta in the πη channel. The relative momenta p and k are related to the invariant K K mass P 0 according to k = 1 2 P 0 p = 1 2 P 0 [ (M K + M 1 K 0)2 ] 1/2 [ (M K + + M 1 K 0)2 ] 1/2 (5) P 2 0 [ (M π M η ) 2 ] 1/2 [ (M π + M η ) 2 ] 1/2 1 1 (6) P 2 0 We now assume that the l th partial wave single differential cross section dσ (l) /dp 0 for the pp dk + K three body decay is given by the proportionality dσ (l) /dp 0 tmax t min P 2 0 dtσ (l) a (π+ η K + K0 ) = C l k 2 p 2 σ(l) a (K+ K0 π + η) t (7) with t equal to the difference in limits of the phase space integration, t = t max t min P 2 0 (s (m D + P 0 ) 2 )(s (m D P 0 ) 2 ); (8) s = (p A + p B ) 2 is the invariant mass of the colliding protons, and C l a constant of proportionality that will depend on specific assumptions about the transfer reaction mechanism. In principle C l could also depend on P 0 : we take it to be a constant for simplicity. For the three body collision the K K invariant mass is thus restricted to the range from threshold (P 0 ) min = (M K + + M K0) to (P 0 ) max = ( s m D ) so that t vanishes at the upper limit; m D is the deuteron mass. At the COSY proton beam energy of T p = 2.65 GeV (p p = 3.46 GeV/c) at which the experiment was performed, these limits are (P 0 ) min = 0.9914 GeV and (P 0 ) max = 1.038 GeV resp. This restricts the K + K0 and π + η CM momenta to lie in the very limited intervals 0 < k < 0.15 GeV and 0.33 < p < 0.36 GeV respectively. With the aid of Eq. (2) the s wave contribution to the cross section can be reformulated as dσ (s) /dp 0 = C s k t p [1 2β sk]; β 2 s = Im(a 1 ) > 0 (9) after redefining C s C s = 4πIm(a 1)C s. The P 0 dependence of the cross section is thus completely determined by the variation of k, p and t with P 0 as fixed by Eqs. (5), (6) and (8). The corresponding expression for p waves is k 3 dσ (p) /dp 0 C p t (10) p2 3
with C p = 12πIm(c 1)C p. The scale constants C s[µb GeV 2 ] and C p[µb GeV 4 ] control the absolute value of the cross section; the parameter β s = Im(a 1 ) determines the skewness of the s wave contribution. The fits to the various data shown in Ref.[1] employ a six parameter transition matrix[7] fitted simultaneously to the two invariant mass distributions K + K0 and d K 0, plus three angular distributions. In order to compare the cross section given by the spectator model with their result for the K + K0 invariant mass distribution, we have determined the pair of s wave constants [C s, β s] by performing a χ 2 fit to the data points in Fig. 1 that lie above threshold, using the sum of Eqs. (9) and (10) with the scale constant C p adjusted to reproduce the p wave partial cross section given in Ref.[1]. One then finds C s = (2.39 ± 0.32)µb GeV 2, together with an imaginary part of the scattering length of β s = Im( a 1 ) = (0.61 ± 0.05) fm at χ 2 min = 0.81, while C p = 16.03µb GeV 4. Note that the uncertainties in the β s parameter also include the systematic uncertainties of the data of Ref.[1] not shown in Fig. 1. The fit to experiment for this parameter set is shown in Fig. 1. By integrating out over the allowed range of P 0, one obtains a total cross section of σ(pp dk + K0 ) = [36 ± 9] nb, in agreement with the experimental cross section [38 ± 2(stat) ± 14(syst)] nb reported in [1]. 3. Flatté parametrization for scattering length As pointed out below Eq. (3) the spectator model only fixes the imaginary part of the scattering length: it does not allow one to extract the real part directly. In order to make further progress one therefore has to introduce specific model dependent assumptions regarding the structure of the K K scattering amplitude in the isovector channel that describe the formation and decay of the a 0 (980). There is little experimental information available on the K K scattering length. However, within the last decade the production of πη pairs have been studied experimentally in proton antiproton annihilation, pion proton reactions and in Φ decay, and analysed using either Flatté distributions[8, 9, 10, 11], or other approaches[12, 13, 14]. The Flatté scattering amplitude[15] near the K K threshold takes on the form S K K 1 iḡ K Kk ǫ a0 + i 2 (ḡ K Kk + ḡ πη p) (11) where k and p are the CM momenta defined in Eqs. (5) and (6), while the a 0 (980) partial decay widths Γ πη = ḡ πη p and Γ K K = ḡ K Kk define the dimensionless coupling constants (ḡ πη, ḡ K K); ǫ a0 = M a0 (M K + + M K0) is the (real or virtual) binding energy of the a 0 (980) relative to the K K threshold. Eq. (11) leads to the following result for the K K isovector scattering length[16] a 1 = 1 2ḡK K ǫ a0 + i 2 Γ πη = R α + i p 1 th (12) when rewritten in terms of the dimensionless ratios[17] R = ḡ K K/ḡ πη and α = 2ǫ a0 /Γ πη instead of the individual coupling constants and energy parameters. Here p th = p[(p 0 ) min ] = 0.33 GeV = 1.65 fm 1 is the threshold momentum for K K production in the πη channel. The R and α are seen to be invariant under a scaling ḡ K K λḡ K K etc. of the original Flatté parameters. Table I shows a compendium of mass, width and coupling constant 4
values for a 0 (980) taken from the recent literature[8, 9, 10, 11, 12, 13, 14] and as partially summarised in [16, 17]. We have also used Eq. (12) to extend the calculations given in [16] for the scattering lengths associated with these paramaters. The results are displayed in order of increasing α in Table I. From Eq. (12) one sees that there is a linear relation between the real and imaginary parts of the (Flatté) scattering lengths, Re( a 1 ) = Im( a 1 )α (13) If we now take Im( a 1 ) = 0.61 fm from the spectator model fit, this fixes the slope of the straight line in Eq. (13). The result can then be compared with the values of Re( a 1 ) calculated from Eq. (12) at the values of α and R given in Table I. This comparison is shown in Fig. 2. The spectator model fit is seen to be consistent with all of the data (except perhaps that of Ref. [13]) once the uncertainties in α and hence Re( a 1 ) due to the experimental error in the a 0 (980) mass are taken into account. It has also been pointed out[17] that direct measurements of the K K elastic scattering cross section near threshold would serve to determine both R and α independently, and thus the scattering length as given by Eq. (12). Such data is not available yet. However, by assuming that scale invariance is valid, one can infer a value for α to be used in Eq. (13) by arguing as follows: the combinations α in Table I would have appeared as a common constant α c had the data satisfied scaling exactly. It is therefore reasonable to estimate this unknown constant by the mean of the α s given in Table I, weighted by the inverse of the average error introduced into each α by the uncertainty in the a 0 (980) mass value used to calculate it. The resultant mean and its standard deviation is then α c = 0.031 ± 0.054. Inserting this estimate into Eq. (13) one obtains Re( a 1 ) = (0.02 ± 0.03) fm. The final result for our extracted isovector scattering length thus reads a 1 = [(0.02 ± 0.03) + i(0.61 ± 0.05)]fm (14) Note that this result relies on two very different sources of information: (i) a direct determination of Im(a 1 ) based on the spectator model fit of the pp dk + K0 data, and (ii), an inferred estimate of Re(a 1 ) based on the assumption that scaling is approximately valid for the experimental parameters extracted assuming a Flatté distribution. Acknowledgments One of us (RHL) would like to thank the Ernest Oppenheimer Memorial Trust for research support in the form of a Harry Oppenheimer Fellowship. Thanks are also due to the Institut für Kernphysik, Forschungszentrum, Jülich, and the Max Planck Institut für Kernphysik, Heidelberg for their kind hospitality. The authors would also like to acknowledge C. Hanhart for critical discussions. 5
References ( ) Present address: Physikalisches Institut, Universität Bonn, Nussallee 12,D 53115, Bonn, Germany. ( ) Permanent address: School of Physics, University of the Witwatersrand, Johannesburg, Private Bag 3, WITS 2050, South Africa. [1] V. Kleber et al., Phys. Rev. Lett. 91, (2003) 172304 1. [2] C. Hanhart, Phys. Rep. 397, (2004) 155. [3] V. Yu. Grishina et al., Eur. Phys. J. A21, (2004) 507. [4] A. Sibirtsev et al., Phys. Lett. B601, (2004) 132. [5] E. Oset, J. A. Oller and U. G. Meissner, Eur. Phys. J. A12, (2001) 435. [6] F. P. Sassen, PhD dissertation, University of Bonn (2004); to be published. [7] V. Yu. Grishina et al., Phys. Lett. B 521, (2001) 217; V. Chernyshev et al., nucl th/0110069; A. E. Kudryavtsev et al., Phys. Rev. C 66, (2002) 015207. [8] A. Aloisio et al. (The KLOE Collaboration), Phys. Lett. B536, (2002) 209. [9] N. N. Achasov and A. N. Kiselev, Phys. Rev. D 68, (2003) 014006. [10] S. Teige et al., Phys. Rev. D 59, (1999) 021001. [11] M. N. Achasov et al., Phys. Lett. B479, (2000) 53. [12] D. V. Bugg, V. V. Anisovich, A. Sarantsev, and B. S. Zou, Phys. Rev. D 50, (1994) 4412. [13] M. Bargiotti et al. (The OBELIX Collaboration), Eur. Phys. J C 26, (2003) 371. [14] A. Abele, et al., Phys. Rev. D 57, (1998) 3860. [15] S. Flatté, Phys. Lett. 63B, (1976) 224. [16] V. Baru et al., Phys. Lett. B 586, (2004) 53. [17] V. Baru, J. Haidenbauer, C. Hanhart, A. Kudryavtsev, and U. G. Meissner, Eur. Phys. J. A23, (2005) 523. 6
FIG. 1: The invariant mass cross section for the production of K + K0 in pp dk K. The solid line denotes a χ 2 fit to the data[1] using the spectator model. The broken curve shows the corresponding fit calculated in [1] using a six parameter transition amplitude squared. Shown separately are the p wave contributions to the cross section obtained from the spectator model, Eq. (10) (solid curve), and the fit (broken curve) given in [1]. The error bars show the statistical uncertainties of the data. 7
FIG. 2: Plot of Re( a 1 ) (filled circles) versus the dimensionless parameter α = 2ε 0 /Γ πη. The plotted points have been tagged by the corresponding references in Table I. The horizontal and vertical error bars show the errors introduced into α, and hence Re( a 1 ), by the uncertainties in the associated a 0 (980) mass values. The heavy straight line shows the calculated value of Re( a 1 ) using Eq. (13) for the particular value Im( a 1 ) = 0.61 fm extracted from the spectator model fit. 8
TABLE I: Summary of calculated values of the isovector scattering length a 1 taken from Eq. (12), arranged in order of increasing α = 2ε 0 /Γ πη, for the various mass, width and coupling constants for the a 0 (980) as extracted[17] from the literature. Ref. M a0 (MeV) ε a0 (MeV) Γ πη (MeV) ḡ πη ḡ K K R α a 1 (fm) [8] 984.8 +1.2 1.2 6.6 121.5 0.373 0.412 1.10 0.108 +0.071 0.66i [9] 992 +14 7 0.6 145.3 0.446 0.560 1.252 0.008 0.006 0.76i [11] 995 +52 10 3.6 125 0.384 1.414 3.685 0.058 0.13 2.23i [12] 999 +5 5 7.6 143 0.439 0.516 1.176 0.106 0.075 0.70i [9] 1003 +32 13 11.6 153 0.470 0.834 1.776 0.152 0.16 1.05i [14] 999 +3 3 7.6 69 0.212 0.222 1.048 0.220 0.13 0.61i [10] 1001 +3 3 9.6 70 0.215 0.224 1.043 0.274 0.16 0.59i [13] 998 +10 10 6.6 46 0.141 0.476 3.376 0.287 0.54 1.89i 9