arxiv: v1 [nucl-ex] 9 Oct 2007

Transcription

1 Kaon Pair Production in Proton Proton Collisions arxiv: v [nuclex] 9 Oct 007 Y.Maeda,, M.Hartmann,, I.Keshelashvili,,3 S.Barsov, 4 M.Büscher, M.Drochner, 5 A.Dzyuba,,4 V.Hejny, A.Kacharava, 3,6 V.Kleber, 7 H.R.Koch, V.Kotev, 4 P.Kulessa, 8 B.Lorentz, T.Mersmann, 9 S.Mikirtytchiants,,4 A.Mussgiller, 6 M.Nekielov,,4 H.Ohm, D.Prasuhn, R.Schleichert, H.J.Stein, H.Ströher, Yu.Valdau,,4 C.Wilkin, 0 and P.Wüstner 5 Research Center for Nuclear Physics, Osaka University, Ibaraki, Osaka , Jaan Institut für Kernhysik, Forschungszentrum Jülich, 545 Jülich, Germany 3 High Energy Physics Institute, Tbilisi State University, 086 Tbilisi, Georgia 4 High Energy Physics Deartment, Petersburg Nuclear Physics Institute, Gatchina, Russia 5 Zentralinstitut für Elektronik, Forschungszentrum Jülich, 545 Jülich, Germany 6 Physikalisches Institut II, Universität Erlangen Nürnberg, 9058 Erlangen, Germany 7 Physikalisches Institut, Universität Bonn, 535 Bonn, Germany 8 H. Niewodniczański Institute of Nuclear Physics PAN, 334 Kraków, Poland 9 Institut für Kernhysik, Universität Münster, 4849 Münster, Germany 0 Physics and Astronomy Deartment, UCL, Gower Street, London WCE 6BT, UK (Dated: November 3, 08) The differential and total cross sections for kaon air roduction in the K K reaction have been measured at three beam energies of.65,.70, and.83gev using the ANKE magnetic sectrometer at the COSY Jülich accelerator. These near threshold data are searated into airs arising from the decay of the φ meson and the remainder. For the non φ selection, the ratio of the differential cross sections in terms of the K and K invariant masses is strongly eaked towards low masses. This effect can be described quantitatively by using a simle ansatz for the K final state interaction, where it is seen that the data are sensitive to the magnitude of an effective K scattering length. When allowance is made for a small number of φ events where the K rescatters from the roton, the φ region is equally well described at all three energies. A very similar henomenon is discovered in the ratio of the cross sections as functions of the K and K invariant masses and the identical final state interaction model is also very successful here. The world data on the energy deendence of the non φ total cross section is also reroduced, excet ossibly for the results closest to threshold. PACS numbers: 3.60.Le, 4.40.Aq, 5.40.E I. INTRODUCTION There are severalmechanisms that can lead to the roduction of kaon antikaon airs in nucleon nucleon collisions near threshold. These can be divided mainly into two general classes involving the roduction of a non strange meson that subsequently decays into K K, and the associated roduction of KY, where the K is formed through the decay of the hyeron Y. The most rominent meson that is found almost exclusively in the K K /K 0 K0 channels is, of course the φ, and there have been measurements of the φ K K reaction from DISTO [] and COSY ANKE [], as well as of the n dφ dk K reaction [3]. These exeriments show strength also away from the region of the φ eak that could arise from other mesons decaying into K K. The lightest of the non strange scalar mesons are the a 0 (I = ) and f 0 (I = 0), which decay mainly into ηπ and ππ, resectively. Both have masses around 980MeV/c and widths of the order of 50 00MeV/c [4]. Since ymaeda@rcn.osakau.ac.j M.Hartmann@fzjuelich.de the central mass values fall almost exactly at the K K threshold, the strong couling to this channel distorts significantly the uer arts of the mass sectra [5]. It has even been suggested that both resonances might be mainly molecular in nature [6]. Since the K K system is a mixture of isosin I = 0 and I =, one can consider whether the non φ events observed in K K might be linked to the roduction and decay of the a 0 and/or the f 0 mesons. The I = channel is isolated cleanly by looking at the dk K0 reaction [7]. On the other hand, we can also exect kaon airs to be created in the K Y reaction, where the hyeron decays through Y K. There are, of course, several excited hyerons that could contribute to such a rocess. Of articular interest for low energy roduction are the Σ(385) and the Λ(405). Though nominally lying below the sum of the K and masses, their large widths ( 50MeV/c ) ensure that they overla the K threshold [4]. Results have recently been resented on the roduction of these states in the K Σ 0 (385) and K Λ(405) reactions, where the hyerons were detected in the Λπ 0 and Σ 0 π 0 channels, resectively [8]. In addition to being involved in a direct roduction mechanism, the Λ(405) and the Σ(385) might also lead to a significant K final state interaction (fsi) in the

2 K K reaction. Evidence that the K fsi is indeed imortant is to be found in the measurement of the K K cross section at energies below the φ threshold [9]. It was shown that the invariant masses of the K final system were on average lower than those of the K airs and it was suggested that this is due to the strong attraction in the K system. It is therefore an interesting but oen question whether the K K data from the non φ region should be viewed, in first aroximation, as scalar resonance roduction via a 0 /f 0 or hyeron roduction through K Λ(405)/Σ 0 (385). It is one aim of the resent aer to study this cometition by using the data obtained when investigating φ meson roduction in roton roton collisions at ANKE []. A rogram to investigate the roduction of kaon airs in nucleon nucleon collisions has been initiated at the COSY storage ring of the Forschungszentrum Jülich. Although the methods of the measurement and analysis of the K K reaction have already been given in Ref. [], for the sake of clarity the main oints are reeated in Sect II. Due to the finite width of the φ, distinguishing between events corresonding to the roduction of this meson and the rest requires a detailed modelling of the K K sectrum, as described in Sect. III. When more than one air of articles interact strongly in the final state, there is no reliable rescrition to evaluate a corresonding enhancement factor. Our aroach to this roblem in terms of an effective K scattering length a is described in Sect. IV, with a heuristic justification being given in Aendix A. The relative distributions of the K and K invariant masses are shown and discussed in Sect. V. The ratio of the differential cross sections as functions of the K and K invariant masses indicates that the K attraction is very strong. This confirms the earlier COSY findings [9] but with higher statistics over much wider ranges of energy and mass. The events that come from the non φ region in the.65gev data allow us to fix the arameter a necessary to describe the K /K ratio within our final state interaction model. The data are mainly sensitive to the magnitude a and the value obtained is not dissimilar to that required to describe free K scattering. The aroach reroduces also our results at the other two energies as well as the COSY data. If allowance is made for the small ossible contribution to this ratio from φ events, of the order of 0%, the two article ratio is also well described for all K K mass intervals. In addition, we find a comletely analogous effect in the threebody K system, with the K /K ratio being strongly biased towards the lowest values of the invariant mass. We show in Sect. VI that these distributions are also redicted quantitatively by the same K final state interaction model with the same value of the arameter a. The total cross sections for φ and non φ roduction are resented in Sect. VII. The energy deendence of the non φ total cross section is strongly influenced by both the and K fsi but, even after including these effects, the lowest energy data are not well reroduced. This leads to the suggestion that the effect might be connected with the two or three oints in the K K sectrum at very low invariant masses, which lie consistently above the simulations. Our conclusions and views on the outlook are resented in Sect. VIII. II. EXPERIMENT AND DATA ANALYSIS The study of kaon air roduction has been erformed using data taken with the ANKE magnetic sectrometer [0]. This is laced at an internal target station of COSY, which is the COoler SYnchrotron of the Research Center Jülich []. We have already reorted an initial analysis of the K K data, from which values of the φ cross sections were extracted []. The measurements were erformed at three roton kinetic energies of T =.65GeV,.70GeV and.83gev, corresonding to excess energies of ε = 5MeV, 67MeV and 08MeV with resect to the K K threshold. A dense hydrogen cluster jet target [] rovided areal densities of cm and the average luminosity during the exeriment was tyically cm s. The detection system of the three diole magnetic sectrometer ANKE registers simultaneously slow ositively and negatively charged articles, with fast ositive articles being measured in the forward system [0, 3, 4]. Kaons roduced in the K K reaction at.65 and.70 GeV have laboratory momenta between 0. and.0gev/c and this full range was covered by the detection system for both K and K. For the.83gev run, the K were detected only u to 0.6GeV/c, which reduced the geometrical accetance. The momentum range of the rotons emerging from the reaction is between 0.8 and.6gev/c, for which the forward detector has full accetance. The momenta of both the kaons and rotons were determined from the track information furnished by sets of multiwire roortional chambers. The tracking efficiency for kaons fell smoothly from 98% to 93% as the momentum increased, whereas that for rotons was about 80% over the whole momentum range. The efficiency ma, including its deendence on the article momenta and osition, was generated for the track reconstruction and used in the analysis. The data taking efficiency was around 90% and, taken in combination with the tracking efficiency, resulted in a total detector efficiency of about 60%. In order not to reduce the geometric accetance unnecessarily, only one forward going roton was measured in coincidence with the charged kaon air. The K K reaction was then identified by requiring that the missing mass in the reaction was consistent with that of the non observed roton. As a first ste, ositive kaons were selected through a rocedure, described in

3 3 detail in Ref. [5], that used the time of flight (TOF) between START and STOP scintillation counters of a dedicatedk detectionsystem. Inthesecondstage,boththe K and forward going roton were identified from the time of flight differences between the STOP counters in the negative and forward detector systems with resect to the STOP counter in the ositive system that was hit by the K. By comaring this with the time difference exected for a kaon air based uon their measured momenta, a clear searation is achieved, with good events lying in an island where the two time differences are consistent [4]. These two TOF selections, as well as that for the K, were carried out within ±3σ bands. derive cross sections. ) Counts/(MeV/c ) Counts/( MeV/c K K invariant mass [GeV/c FIG. : The K K invariant mass distribution for the K K reaction at.65gev. The rominent φ eak sits on a steadily falling background corresonding to the ANKE accetance for the four article final state. The simulation of these two contributions and their sum are shown, resectively, by the dotted, dashed, and solid histograms. ] Missing mass (K K ) [GeV/c ] FIG. : Missing mass distribution of K K X events at T =.65 GeV. The hatching indicates the range used for the selection of rotons. The side band events scaled by the solid line are indicated by the dotted shading. The sectrum of missing masses with resect to the K K system, shownin Fig.forthe.65GeVdata, is dominated by the roton eak. The corrected numbers of K K eventsobtainedinthiswayarearound300(at T =.65GeV), 300 (.70GeV), and 650 (.83GeV). The estimated background inside the roton cut window is 5%, %, and 8% at these three energies, resectively. This has been subtracted using the side band events scaled by the solid line shown in Fig.. The uncertainty resulting from this rocedure is estimated to be less than 3% for the total cross sections and negligible comared to the statistical errors for the differential distributions. After the identification of the K K events, the data were binned in intervals of the K K invariant mass, with the.65gevresultsbeing resentedin Fig.. This sectrum shows a very rominent φ eak ositioned closeto the nominal massofthe meson. In addition there is a non φ contribution that decreases steadily from low to high invariant masses. This behavior, which is mainly a reflection of the ANKE accetance for non resonant K K roduction, has to be simulated in order to The accetance has been determined through a Monte Carlo simulation based uon GEANT4 [6], where the geometrical accetance, resolution, detector efficiency, and kaon decay robability were taken into account. In the initial ste, the φ roduction distributions reviously ublished [] were taken as the basis of the simulation, whereas four body hase sace was used for the non φ comonent of Fig.. There are seven degrees of freedom for a fourbody final state roduced in an unolarized reaction. Thus for each beam energy, seven indeendent c.m. distributions were generated (five angular distributions, the relative momentum of final rotons, and the K K invariant mass). These were then divided into two grous deending on the value of the K K invariant mass, i.e., φ rich and φ oor regions. All distributions were jointly fitted to the exerimental data and the relative contribution of φ and non φ roduction evaluated for the determination of the accetances. In the φ oor region, most of the c.m. distributions are fairly uniform, though the olar angle of the K K system relative to the beam axis in the overall c.m. system shows deviations. Fitting this with the function αcos θ gives anisotroies of α = 0.5 ± 0. and α =.5 ± 0.5 at.65gev and the two higher energies, resectively. The inclusion of this factor increases the total accetance for non φ roduction by about 0%. In addition, as will be seen in Sects. V and VI, the K ± and K ± invariant masses deviate strongly from hase sace. These deviations have been taken into account in the simulations, using iteratively the theoretical aroach discussed in Sect. IV, in order to converge on the accetance corrected distributions. The simulations of

4 4 the φ and non φ comonents of the uncorrected K K sectrum at.65gev are shown in Fig. together with their sum. The estimated total accetance for kaonair roduction is.% (.65GeV), 0.7% (.70GeV) and 0.% (.83 GeV), and the associated systematic uncertainties in the total cross section are 5%, 7% and %, resectively. The shaes of the differential cross sections, esecially those for the K and K invariant masses, are very stable to this iterative rocedure. As a consequence, the uncertainty that this introduces in the differential cross sections is negligible comared to the statistical errors. The luminosity required to evaluate absolute cross sections was determined by measuring in arallel roton roton elastic scattering for laboratory angles between 5.0 and 8.5. For this urose the momentum of a forward going roton was determined using the ANKE forward detector. After a missing mass analysis, the (small) background underneath the roton eak was subtracted and the remaining events were corrected for efficiency and accetance. ] dσ/dm [µb/gev c III. THE K K INVARIANT MASS DISTRIBUTION 5 (c) Figure 3 shows the K K differential cross sections in terms of the K K invariant masses as obtained at our three energies. These have been corrected for accetance, deadtime, etc., and normalized on the basis of the roton roton elastic scattering data, as discussed in Sect. VII. The resulting non φ contribution is much flatter than that illustrated in Fig.. In order to ensure a clean searation of the φ from the rest, it is necessary to model the shaes as reliably as ossible. In our earlier work [], the non φ contribution was assumedtobedescribedbyafour bodyhasesace. Theφ was taken to have a Breit Wigner form with a constant width of Γ = 4.6MeV/c [4], which was then convoluted with an exerimental resolution of σ = MeV/c. Only for the φ comonent was the sectrum modified to include effects arising from the final state interaction between the emerging rotons. These three assumtions have now to be reconsidered. Since the φ is a wave resonance that decays strongly into K K, its width deends uon the K K and K 0 K0 relative momenta. Due to the narrowness of the state, the eak itself is little affected by this refinement but it does suress slightly the low mass tail of the resonance. Although for low K K effective masses there is no clear evidence for any fsi, only minor effects are to be exected there because of the large amount of energy that is then available to excite the system. Nevertheless, the fsi cannotbe neglectedathigherk K masses, where its inclusion distorts slightly the four body hase sace sectrum. Finally, as will be seen in Sect. V, it is clear that there is a large final state interaction between K K invariant mass [GeV/c FIG. 3: Differential cross section for the K K reaction with resect to the K K invariant mass at.65 GeV,.70 GeV, and (c).83 GeV. The dotted histogram and the dashed curve reresent the simulations for the φ and non φ contributions to the sectrum, with their sum being shown by the solid histogram. For all three energies the simulations fail to reroduce the lowest mass oints. the K and one or both rotons and this also modifies the shae of the non φ comonent, enriching the central art of the K K sectrum. By chance, all three effects go in the same direction. In the fits to Fig. 3, the fractions to be associated with φ roduction are diminished somewhat comared to those reorted in our revious work []. It must be noted that the lowest mass oints lie above the fitted curves for all three data sets of Fig. 3. These deviations from the simulations are significant because, asseen from the comarisonoffigs. and 3a, the ANKE accetance is largest in this region. Although ignored in ]

5 5 the fitting rocess, they were taken into account in the evaluation of the non φ contribution to the total cross section, contributing 4% (.65 GeV), 6% (.70 GeV), and 0% (.83GeV). IV. FINAL STATE INTERACTION ASSUMPTIONS dσ/d [µb/(gev/c)] 0. We wish to describe the K K data at all energies in terms of final state interactions between the two rotons and between the K and either one or both of the rotons. There is no reliable wayofcarryingout such a rogram without solving multibody equations, which is well outside the scoe of the resent aer. As exlained in Aendix A, we make the ad hoc assumtion that the overall enhancement factor is the roduct of the enhancements in the and two K systems, all evaluated at the aroriate relative momenta q: F = F (q ) F K (q K ) F K (q K ). () Similar aroaches are to be found elsewhere in the literature as, for examle, in the descrition of the η final state at low energies [7]. The K enhancement factor is taken in the scattering length aroximation F K (q) = iqa, () where some caution should be exercised in the interretation of the comlex arameter a as an effective scattering length. In rincile a could have an energy deendence, which would corresond to the introduction of an effective range term. However, it is hard from our data to determine the values of any extra free arameters. A roton roton enhancement factor of the form F (q) = q β q α (3) isassumed[8], wheretheositionofthe S 0 virtualstate is well fixed at α = 0.fm. The final state interaction should not influence higher artial waves and, in order to reduce its effects above about 0MeV, we take β = 0.5fm. The relative momentum of the two rotons in the rest frame, which is very sensitive to the form of the fsi, is then reroduced reasonably, as seen from Fig. 4. It is imortant to note that the form of neither the northe K enhancementfactorisvalidforanexcitation energy below about MeV due to the exlicit Coulomb force. However, these regions reresent only very small fractions of the allowed hase sace. V. THE KAON PROTON INVARIANT MASS DISTRIBUTION In data taken with the COSY sectrometer at ε = 8MeV, i.e. at an energy where there is very lit Momentum [GeV/c] FIG. 4: Relative momentum of the final rotons from the K K reaction in the rest frame. The exerimental oints at.65gev are comared to the simulation, using the standard arameters, which is shown as a histogram. tle φ roduction, it was found that the K and K invariant mass M(K) distributions were markedly different [9]. The ratio of the accetance corrected distributions R K = dσ/dm(k ) dσ/dm(k ) showedaverystrongreferencefor lowvalues ofm(k). Figure 5 shows the K sectra obtained at our lowest energy. In order to consider a situation similar to that of COSY, the data have been selected such that only events with K K invariant masses less than.0gev/c are retained. The simulation that takes into accountonlythefsi failstodescribethedata. Furthermore, the ratio of the K to K mass distributions is far from constant, as it would be if there were only the fsi. The eaking of R K to the lowest invariant masses confirms the COSY observation [9] but with higher statistics and over a wider range of K masses. When the K final state interaction factors are also introduced into the four body hase sace simulation of the non φ contribution to the K K reaction, the individual dσ/dm(k ) and dσ/dm(k ) distributions of Fig. 5, as well as their ratio R K, are all well described. The best fit to these data is achieved with an effective K scattering length of a = (0.5i)fm and this value of a will be retained for all the subsequent simulations. Nevertheless, the uncertainties in the real and imaginary arts are quite large and strongly correlated such that a is determined much better than its hase, as illustrated by the χ contour lot of Fig. 6. The free K I = 0 and I = values from Ref. [9] are also shown there along with their one σ limits. From this it is seen that there is an overla with either isosin scattering length. One has also to bear in mind that our

6 6 ] dσ/dm [µb/gev c Im(a K _ ) a 0 a Re(a K ) _ FIG. 6: Contour lot corresonding to the one σ level uncertainty of the fits of the comlex arameter a of the fsi ansatz to the data of Fig. 5. Our best fit is obtained with a =.5ifm (marked with a cross). Also shown are the best fit values (oints) and corresonding contours (dashed lines) for the isoscalar (a 0) and isovector (a ) K scattering lengths derived from the study of free K nucleon scattering [9]. / (c) K invariant mass [GeV/c ] FIG. 5: Differential cross sections for the K K reaction with resect to the invariant mass of the K and the K systems. In anel (c) their ratio is shown on a logarithmic scale. These.65 GeV data have been selected to come from the region M(K K ) <.0GeV/c, where there is only a very small φ contribution. The dashed histograms reresent the results of four body hase sace simulations, where the only distortion is that coming from the fsi. The solid histograms include also the K fsi. data reresent averages over some energy range and that a constant scattering length is necessarily an oversimlification. The situation for higher K K masses is comlicated by the contribution from the φ meson. At our beam energies, the average searation of the K (roduced in the φ decay) from one of the final rotons is about 7fm. This large distance reduces considerably the effects of the K fsi. We would therefore exect that deviations of R K from unity should arise rimarily from the nonφ contributions. The modelling of the data thus deends in the first instance uon the determination of the fraction of φ roduction in the overall K K sectra resented in Fig. 3. However, classical simulations suggest that more than 5% of the φ decays lead to at least one K air with a searation of less than fm, for which the fsi is not negligible. To investigate this oint further, the.65gev data have been divided into four regions of K K invariant mass, viz. M(K K ) < 0.995GeV/c, < M(K K ) <.003GeV/c,.003 < M(K K ) <.00GeV/c, and.00gev/c < M(K K ). Thefirst two intervals contain almost no φ events, whereas the last corresonds mainly to φ roduction. The corresonding R K distributions are shown in Fig. 7, where they are comared to the same simulation as that shown in Fig. 5. The nonφ region is very well described but, for high K K invariant masses, the asymmetry seems to be diluted too much by kaons from the φ decay. The agreement is imroved by assuming that 0% of the φ events are influenced by the K fsi. The same effect could be achieved by artificially reducing the fraction of φ mesons in the eak of Fig. 3a by 0%, but this would be hard to justify. If the strong K final state interaction is resonsible for the distortion in the K /K ratio seen in Fig. 7, similar effects should be seen for the other excess energies. For our ε = 66.6 and 08MeV results shown in Fig. 8, the statistics are lower and we merely divide each data set into two, comrising φ rich and φ oor regions. Without any change in the value of the K effective scattering length, all the distributions are well described, esecially if it is assumed that 0% of the φ events lead to K fsi effects. Data on the K /K ratio were not ublished by

7 7 / /.45.5 (c) / (c) (d) (d) K invariant mass [GeV/c ] K invariant mass [GeV/c ] FIG. 7: Values of the R K ratio for the K K reaction at.65 GeV extracted for four different selections of the K K invariant mass, viz. M(K K ) < 0.995GeV/c, < M(K K ) <.003GeV/c, (c).003 < M(K K ) <.00GeV/c, and (d).00gev/c < M(K K ). The histograms are redictions of the model described in the text, where the K mesons from non φ events are subject to a K final state interaction. In the high mass interval, the effects of assuming a 0% robability of the K from the φ decay being influenced by a fsi is considered and this leads to the dashed histograms. the DISTO collaboration []. However, the same good descrition of R K is found for the COSY results of Fig. 9 where, in view of the low excess energies, there is no ambiguity regarding the φ contribution. Given that the systematic uncertainties in this exeriment are indeendent of those at ANKE, the quantitative agreement FIG. 8: Values of the R K ratio for the φ oor region and the φ rich at.70gev with (c) and (d) showing resectively the analogous results at.83 GeV. The redictions obtained using the standard arameters are shown with no fsi effect for the K from the decay of the φ (solid curve) and 0% (dashed). shown by the 8 MeV data is convincing. Unfortunately, in the ublication [9] the values of the ratio at 0MeV were scaled uwards to make them lie closer to those taken at 8MeV. Since the numbers of K and K events are the same, the weighted average of the ratio must be unity and this allows one to make an informed guess regarding the value of the scaling that has to be reintroduced. No account has been taken of the COSY invariant mass resolution, which will be of greater significance for the 0MeV data.

8 8 Ratio 0 mass deendence which is aarent in Fig. 0 does not necessarily imly that the K forms a bound state with the two rotons [0,, ]. Though the error bars are considerably larger at the two higher energies, as seen in Fig., the simulation is consistent with the results. In all these three article mass distributions the introduction of a small amount of fsi for the roducts of the φ decay seems to be of less imortance than for the K sectra. Ratio K invariant mass [GeV/c ] FIG. 9: The R K ratio for the COSY K K data at ε = 0MeV (uer anel) and ε = 8MeV. The small contribution from the low mass tail of the φ has been neglected in the 8MeV simulation. Since the 0MeV results were arbitrarily scaled in the ublication [9], they have here been divided by a factor of three before being comared to the redictions of the fsi model, using the same values of the arameters as for the ANKE data. VI. THE KAON PROTON PROTON INVARIANT MASS DISTRIBUTION We have argued in Sect. IV and Aendix A that it is natural to take the final state enhancement to be the roduct of the individual factors corresonding to the and two K interactions, which allows the kaon to be attracted simultaneously to both rotons. This in turn suggests that there might also be an enhancement in the three body system at low K invariant masses. This is indeed the case. The size of the effect at.65gev is illustrated in Fig. 0 for the regions with M(K K ).0GeV/c. There is a striking reference for low K masses in both K K intervals and this behavior is reroduced quantitatively by our simulation, which uses the same K effective scattering length a that was chosen to describe the two body R K ratio. It rovides some emirical justification for our choice of the final state interaction enhancement factor as the roduct of the individual two body enhancements. This favors configurations where the K is attracted simultaneously to both rotons. It must, however, be stressed that the strong VII. TOTAL CROSS SECTIONS The total cross section results for both φ and non φ roduction are given in Table I for the three energies. The systematic errors given here arise from the quadratic sum of the uncertainties in the background subtraction, accetance correction, tracking efficiency correction for kaons (3%) and the analysis of the elastic data used for the normalization. Unfortunately, the database is very limited in our region of angle and energy. This is recognized in the recent SP07 udate to the SAID rogramwhereitisstatedthatour solution should be considered at best qualitative between.5 and 3 GeV [3]. The raw SP07 solution generally overestimates the values of those differential cross sections that are available in this range [4]. A normalization otion has to be chosen to ensure agreement with the exerimental data and this tyically leads to reductions ofabout 5% at our two lower energies but almost 0% at.83 GeV. To eliminate any ambiguities, the values that we have assumed for the elastic cross section, integrated over the solid angle defined by 5.0 < θ lab < 8.5, are resented in Table I for our three energies. TABLE I: Total cross sections for the {K K } nonφ and φ reactions at three beam energies T and corresonding excess energies ε KK with resect to the K K threshold. The excess energy with resect to the φ threshold is given by ε φ = ε KK 3.MeV. The φ cross section has been corrected for the 49.% branching ratio [4]. For both channels the first error is statistical and the second systematic. The values thatwe have takenfor σ, the integral of the roton roton elastic cross section over the solid angle range 5.0 < θ lab < 8.5, are also given. The associated ±6% uncertainty has not, however, been comounded with the other errors. T ε K K σ σ nonφ(tot) σ φ (tot) [GeV] [MeV] [mb] [nb] [nb] ±± 33±± ±±3 64±4± ±8±5 33±±7 The value that we obtain for the luminosity at.65gev on the basis of the cross section given in the table could be checked using a novel technique [5]. This deends uon the energy loss of the coasting roton beam caused by its multile traversals through the

9 9 )] dσ/dm [µb/(gev/c / 0. )] dσ/dm [µb/(gev/c / invariant mass [GeV/c K ] K invariant mass [GeV/c ] K invariant mass [GeV/c ] FIG. 0: Differential cross sections for the K K reaction at.65gev with resect to the invariant masses of the K and the K systems and their ratio. The data are divided into a φ oor region with M(K K ) <.0GeV/c (uer anels) and φ rich with M(K K ) >.0GeV/c (lower anels). The redictions obtained using the standard arameters are shown with no fsi effect for the K from the decay of the φ (solid curve) and 0% (dashed). The simulations that describe the two article ratios also reresent well these three body ratios. very thin target. The resultant frequency change can be measured with high accuracy using the sectrum of the Schottky noise. Combined with measurements of the beam current, this yielded luminosities with an exected recision of about ±6%. The two methods give consistent results within this uncertainty and, although the frequency change technique could not be alied as reliably at the higher energies, it is exected that the overall luminosity uncertainty is on about this level. Our total cross section results for non φ roduction are lotted in Fig. along with results taken from DISTO [] and COSY [9, 6, 7]. It is seen that four body hase sace cannot describe simultaneously the K K data at high and low excess energies. The situation is imroved only a little if account is taken of the strong attraction between the two rotons in the S 0 state but a greater imrovement is achieved through the introduction of the fsi in the K system. Nevertheless, the COSY data at excess energies of 0MeV or less still seem to be underestimated. This might be connected to the incomlete descrition of the K K sectra of Fig. 3 at very low invariant masses. However, itmustbestressedthatintheevaluationofthecosy accetance no account was taken of any K fsi. VIII. CONCLUSIONS AND OUTLOOK We have investigated the K K reaction at three energies. In those arts of the K K mass sectra where φ roduction is small, we find a marked difference between the K and K invariant mass distributions. There is a strongeaking in the K /K ratio towards the K threshold for both the COSY and our data. This clearly indicates that the K K reaction cannot be dominated by the undistorted roduction of a single scalar resonance a 0 or f 0. Either a mixture of the two is roduced or one of the kaons from the a 0 /f 0 decay interacts with a final roton. As we have shown, these data can be exlained quantitatively through the introduction of a simle K final state interaction. In kinematic regions where φ roduction is significant, the K /K asymmetry is reduced because the meson may travel some distance before decaying. Our data in such conditions are well fit by assuming that about 0% of the K from the φ decay interact with one or both final rotons. The uncertainty in this is, however, quite large because of the limited recision with which the φ contribution to the K K mass sectrum can be isolated. There is no universal technique for modelling the fsi when more than one air of articles is involved. The

10 0 / [nb] σ tot [MeV] 0 /.4.4 (c) FIG. : Energy deendence of the nonφ contribution to the K K total cross section. In addition to our three oints (closed circles), there is a high energy oint from DISTO (oen circle) [] and four oints close to threshold from COSY (closed squares) [9], (oen triangle) [6], and (oen square) [7]. The dashed curve reresents the energy deendence from four body hase sace. The dotdashed includes the effects of the final state interaction whereas the solid curve contains in addition distortion from the K fsi, as described in the text. (d) K invariant mass [GeV/c ] FIG. : The ratio of the differential cross sections of the K K reaction with resect to the K to K invariant masses for the φ oor region and the φ rich at.70gev with (c) and (d) showing resectively the analogous results at.83 GeV. The redictions obtained using the standard arameters are shown with no fsi effect for the K from the decay of the φ (solid curve) and 0% (dashed). The simulations that describe the two article ratios also reresent well these three body ratios. roduct ansatz has the benefit that the same simulation with the same arameters reroduces also the strong eaking in the ratio of the K /K mass distributions. We find an effective K scattering length that is of the same order of magnitude as that needed to describe the free scattering data [9]. However, the error bars are large and our data are rimarily sensitive to a rather than its hase. It is imortant to stress that we do not know the relative strengths for the roduction of isosin I = 0 and I = K airs in the K K reaction and hence the relative weights of a 0 and a. Some information on this might come from looking at the data on the K Σ 0 (385) and K Λ(405) reactions taken at.83 GeV [8]. However, even this would not give access to the hase between the I = 0 and I = amlitudes. It must also be stressed that we have neglected any energy deendence of the arameter a. Another imortant caveat is that we have not considered any ossible interaction of the K with the rotons. Although there is no evidence for any strong effect of this nature in the K Λ reaction [8], since the K interaction is reulsive, its neglect might be interreted as extra attraction in the K system. We have clearly shown the imortance of the K interaction, which is ossibly related to the Λ(405) or the Σ 0 (385), but we must now ask whether there is any evidence for effects that might be connected with the roduction of the a 0 /f 0 resonances. The simulation of the energy deendence of the total cross section for the roduction of the non φ comonent in the K K reaction underestimates the COSY data very near threshold. Even more intriguing, and ossibly connected, are the lowest mass oints in the K K sectra. At all three energies shown in Fig. 3, these oints lie much higher than the simulations and a similarbehavioris visible in the DISTO results [] aswell as in our data on the n dk K reaction [3]. The

11 mass scale ofthe variationis not that of the widths ofthe scalar resonances, which are quite large. It is temting to suggest that this structure might be due to the oening of the K 0 K0 channel at a mass of 8MeV/c, which induces some cus structure that changes the energy deendence of the total cross section near threshold [9, 9]. This would require a very strong K K K 0 K0 channel couling, which might be driven by the a 0 /f 0 resonances. Thus, although the K K reaction may not be ideal for investigating the roerties of scalar states, their indirect effects might still be crucial. Given that the K is strongly attracted to rotons, it is natural to exect that similar fsi effects should exist for other reactions. The K 0 d/k d ratio in the K K0 d reaction does indeed show a very strong enhancement near the Kd threshold [7]. The interretation is, however, slightly more comlex because sin arity constraints means that there must be at least one wave in the final state, even at low energies. Unfortunately, the MOMO data on d K K 3 He [30] cannot be used to isolate any K 3 He interaction because the charges of the kaons were not identified and only averaged sectra could be studied [3]. There is currently an overall uncertainty in our values of the cross sections because of the limitation in the roton roton elastic scattering database in our region of energy and angle, though we hoe that this will be alleviated by future measurements [5]. With our new estimates for the luminosity, the φ roduction cross sections are reduced comared to those of our first analysis [], which relied uon the non renormalized 004 SAID solution [3]. Furthermore, the better searation of the φ and non φ events goes by chance in the same direction. These two effects lower the ratio of φ to ω roduction in roton roton collisions near threshold so that it is now only about a factor of six above the Okubo Zweig Iizuka limit [3]. Acknowledgments We would like to thank the COSY machine crew for their suort as well as that of other members of the ANKE Collaboration. Discussions with J. Haidenbauer, C. Hanhart and A. Sibirtsev were very helful and those with I. Strakovsksy rovided imortant clarifications on the uncertainties in the SAID rocedure above.5gev. R. Nißler and U. G. Meißner kindly rovided us with the redictions of Ref. [9] in numerical form. Corresondence with R. Bertini and P. Winter is gratefully acknowledged. This work has been suorted in art by the BMBF, DFG, Russian Academy of Sciences, and COSY FFE. APPENDIX A: FINAL STATE INTERACTIONS When two articles interact strongly in the final state, the resulting matrix element involves an average of a roduction oerator with the relative wave function ψ q ( r) of the strongly interacting air. If the interaction is of very short range, the wave function may be evaluated at the origin to leave an enhancement factor F (q) ψ q (0) D(q), (A) where q is the relative momentum in the air and D(q) is the S wave Jost function [8]. In the commonly used scattering length aroximation, one retains only the linear term in q, in which case F(q) = /( iqa), where a is the scattering length of the interacting air. On the other hand, if two or more airs of articles interact in the final state, there is no reliable rescrition to evaluate an analogous enhancement factor since a three body equation then needs to be solved. If we denote the interacting airs as and 3, the corresonding wavefunction will be Ψ q, q 3 ( r, r 3 ). We nowmakethe ad hoc assumtion that this wave function factorizes in the form Ψ q, q 3 ( r, r 3 ) = ψ q ( r ) ψ q3 ( r 3 ). (A) In this case, the three body enhancement factor is simly the roduct of the two body factors, evaluated at the aroriate relative momenta: F 3 (q,q 3 ) = F (q ) F (q 3 ) ( iq a )( iq 3 a 3 ). (A3) It must be stressed that Eq. (A3) is merely an ansatz to try to understand our ensemble of data and, even if it rovides a satisfactory descrition of these, this does not mean that it is alicable more generally. However, it has also been used to treat the η reaction near threshold, where all the final airs of articles interact strongly [7]. The aroach does retain the necessary ole structure when q = i/a and similarly for q 3. Thus article can interact simultaneously with both and 3. Furthermore, to lowest order in the momenta, the ansatz corresonds to the scattering length of article from a comosite 3 with the desired combined scattering length of a = a a 3. Although the S wave interaction is well known, the K interaction is far more comlex because of the channel coulings to Σπ and also because there are two isosins I = 0 and I =. In a recent study, within a chiral SU(3) unitary aroach [9], K data leading to a variety of channels have been fit and values of the (comlex) K scattering lengths deduced for I = 0 and I =. Their full model resulted in a 0 = (.64i0.75)fm, a = ( 0.06i0.57)fm. (A4)

12 These values describe the raid dro in the K elastic amlitudes from threshold. If isosin breaking induced by the significant mass differences is neglected, the K scattering length a K = (a 0 a )/. However, this is not necessarily the arameter that is relevant for the analysis of the K {K } because it is not clear whether the reaction mechanism referentially excites I = 0 or I = {K } statesor some combinationthereof. The values in Eq. (A4) should therefore be considered merely as order of magnitude estimates of what might be exected when making fits to our data. [] F. Balestra et al., Phys. Rev. C 63, (00). [] M. Hartmann et al., Phys. Rev. Lett. 96, 430 (006). [3] Y. Maeda et al., Phys. Rev. Lett. 97, 430 (006). [4] W. M. Yao et al., J. Phys. G 33, (006). [5] See the short note on scalar mesons by S. Sannier and N.A Törnquist in S. Eidelman et al., Phys. Lett. B 59, 506 (004). [6] J. Weinstein and N. Isgur, Phys. Rev. D 4, 36 (990). [7] V. Kleber et al., Phys. Rev. Lett. 9, 7304 (003); A. Dzyuba et al., Eur. Phys. J. A 9, 45 (006). [8] I. Zychor et al., arxiv: [nuclex]. [9] P. Winter et al., Phys. Lett. B 635, 3 (006). [0] S. Barsov et al., Nucl. Instrum. Meth. Phys. Res., Sect. A 46, 364 (00). [] R. Maier et al., Nucl. Instrum. Meth. Phys. Res., Sect. A 390, (997). [] A. Khoukaz et al., Eur. Phys. J. D 5, 75 (999). [3] S. Dymov et al., Part. Nucl. Lett., 40 (004). [4] M. Hartmann et al., Int. J. Mod. Phys. A, 37 (007). [5] M. Büscher et al., Nucl. Instrum. Meth. Phys. Res., Sect. A 48, 378 (00). [6] S. Agostinelli et al., Nucl. Instrum. Meth. A 506, 50 (003); htt://geant4.web.cern.ch/geant4. [7] V. Bernard, N. Kaiser, and U. G. Meißner, Eur. Phys. J. A 4, 59 (999). [8] M.L. Goldberger and K.M. Watson, Collision Theory, (Wiley, N.Y., 964) [9] B. Borasoy, U. G. Meißner, and R. Nißler, Phys. Rev. C 74, 0550 (006). [0] M. Agnello et al., Phys. Rev. Lett. 94, 303 (005). [] T. Yamazaki and Y. Akaishi, arxiv: [nuclth]. [] N.V. Shevchenko, A. Gal, J. Mareš, J. Révai, arxiv: [nuclth]. [3] R.A. Arndt, I.I. Strakovsky, and R.L. Workman, Phys. Rev. C 6, (000); R.A. Arndt, W.J. Briscoe, I.I. Strakovsky, and R.L. Workman, Phys. Rev. C 76, 0509 (007); SAID solution SP07, available from htt://gwdac.hys.gwu.edu. [4] D.Alberset al., Eur.Phys.J.A, 5(004); R.C.Kammerud et al., Phys. Rev. D 4, 309 (97); S.A. Azimov et al., Sov. Phys. JETP 5, 99 (96); I. Ambats et al., Phys. Rev. D 9, 79 (974); T. Fujii et al., Phys. Rev. 8, 836 (96); W.M. Preston, R. Wilson, and J.C. Street, Phys. Rev. 8, 579 (960). [5] H.J. Stein et al., in rearation. [6] M. Wolke, PhD thesis, University of Münster (997). [7] C. Quentmeier et al., Phys. Lett. B 55, 76 (00). [8] S. Abd ElSamad et al., Phys. Lett. B 63, 7 (006). [9] W. Oelert et al., Int. J. Mod. Phys. A, 50 (007). [30] F. Bellemann et al., Phys. Rev. C 75, 0504 (007). [3] V.Yu. Grishina, M. Büscher, and L.A. Kondratyuk, Phys. Rev. C 75, 0508 (007). [3] S.Okubo, Phys. Lett. 5, 65 (963); G.Zweig, CERN reort TH40 (964); J. Iizuka, Prog. Theor. Phys. Sul. 38, (966).