Energy dependence of forward 1S0 diproton production in the pp ->pp pi0 reaction

The pp ->{pp}_s pi0 differential cross section has been measured with the ANKE spectrometer at COSY-Juelich for seven proton beam energies T_p between 0.5 and 1.97 GeV. By selecting proton pairs with an excitation energy of less than 3 MeV it is ensured that the final {pp}_s system is in the 1S0 state. In the measured region of theta_{pp}^{cm}<18 deg, the data reveal a forward dip for T_p =<1.4 GeV whereas a forward peaking is seen at 1.97 GeV. The energy dependence of the forward cross section shows a broad peak in the 0.6-0.8 GeV region, probably associated with Delta(1232) excitation, and a minimum at 1.4 GeV. Some of these features are similar to those observed for the spin-isospin partner reaction, pp ->d pi+. However, the ratio of the forward differential cross sections of the two reactions shows a significant suppression of single pion production associated with a spin--singlet final nucleon pair.

Single pion production in nucleon-nucleon collisions, NN → NNπ, is one of the principal tools used in the investigation of NN dynamics at intermediate energies [1][2][3]. Because of the large momentum transfers involved, even close to threshold, such a meson production process is sensitive to the short-distance part of the NN interaction.
The pp → dπ + reaction has been the subject of extensive experimental study with the measurement of many spin observables, as well as of the differential cross section from threshold up to several GeV. However, the information that this provides is restricted to final NN states with spin S = 1 and isospin T = 0. On the other hand, the pp → {pp} s π 0 process is kinematically very similar to this provided that the excitation energy in the final proton pair is very small. In this case, due the the Pauli principle, the protons must be in the singlet 1 S 0 state, i.e. have quantum numbers (J P , T, L) = (0 + , 1, 0) compared to the (1 + , 0, 0 and 2) of the deuteron. Pion production in the two cases therefore involves different transition matrix elements so that a combined study of the two processes should yield greater insight into the reaction dynamics.
If we consider the pp → dπ + reaction as the limit of triplet {pn} t production in the pp → {pn} t π + channel, where the strong final state interaction produces the deuteron [4], the ratio of π 0 to π + cross sections will provide information on the relative strength of spin-singlet to spin-triplet production, i.e. give information on the relative probability of pion production at short distances in channels with different spin orientation of the final nucleons. Because of the smallness of the signal, attempts to identify spin-singlet production directly from data on the pp → pnπ + cross section have only yielded upper limits [5][6][7].
A small value of the singlet-triplet ratio is expected near threshold since swave isovector pion rescattering is absent for pp → ppπ 0 and heavy (ω) meson exchange provides the largest driving term [8]. A small value of the ratio is also predicted for energies around 0.4-0.6 GeV since the S-wave ∆(1232)N intermediate state that dominates the pp → dπ + cross section is forbidden by conservation laws in the pp → {pp} s π 0 case. The theoretical situation at higher energies is largely open. The position is rather similar on the experimental side since, away from the low energy domain, the only other published data in the 1 S 0 conditions were limited to energies T p ≤ 0.425 GeV [9].
We have previously reported a measurement of the pp → {pp} s π 0 differential cross section obtained using the ANKE spectrometer [10] for the single beam energy of 0.8 GeV [11] . Here the two protons were detected at small angles with respect to the incident beam and cuts were made such that the excitation energy E pp in the final pp system was less than 3 MeV. Under these conditions we expect the final {pp} s pair to be almost purely in the singlet 1 S 0 state. It was found for this beam energy that the angular variation was rather strong but that at all the measured angles the cross section was orders of magnitude smaller than that for pp → dπ + . To study the energy dependence of these effects, we present here further measurements taken over a wide range of energies, 0.508, 0.625, 0.700, 1.100, 1.400, and 1.970 GeV.
The magnetic spectrometer ANKE is placed at an internal beam station of the COSY cooler synchrotron of the Forschungszentrum Jülich. Fast charged particles, resulting from the interaction of the proton beam with the hydrogen cluster-jet target [12] and passing through the analysing magnetic field, were recorded in the forward detector (FD) system [13]. The FD system includes multiwire proportional chambers for tracking and a scintillation counter hodoscope for energy loss and timing measurements. The triggers employed required the crossing of the two planes of the hodoscope by at least one charged particle (SP-trigger) or by two particles (DP-trigger) [14]. When the DPtrigger was used for data taking, the SP-trigger ran (prescaled) in parallel for luminosity measurement and calibration purposes. The tracking system provided momentum resolution σ p /p ≈ 1% in the range of interest and resolution in the excitation energy of σ(E pp ) ≈ 0.2 − 0.8 MeV for events with E pp < 3 MeV.
Additional details of the experimental setup and the measurement procedure are to be found in refs. [11,13,15,16]. The data at 0.8 GeV and above were taken during a single beam-time run, whereas the lower energy results were obtained from other ANKE calibration runs.
From measurements of the momenta of two charged particles in ANKE, the pp → ppπ 0 channel was isolated by determining the missing mass M x in the reaction. In more than 80% of cases where two fast particles were detected, the tracks passed through different counters of the forward hodoscope. For these events the particles could be clearly identified as protons on the basis of the timing information. The difference of the arrival times of the two particles measured by the counters was compared with the time-of-flight difference deduced from the measured momenta assuming that the particles both had the mass of the proton. For the remaining ≈ 20% of events, the hypothesis was made that the two particles were indeed protons. As already shown for the 0.8 GeV data [11], the missing-mass distributions for both classes of events are very similar with only a slightly enhanced background when no timing information was available. The two sets were therefore combined in the subsequent analysis. Figure 1 presents examples of the measured missing-mass distributions in the π 0 region for events with E pp < 3 MeV. In addition to the π 0 peak, a rise of counts is seen on the right hand side due to two-pion production which seems to be largest at 1.1 and 1.4 GeV. This background gives little contribution in the π 0 region and the data were fitted as a sum of a Gaussian and straight line. The region close to M 2 x = 0 was excluded from the fit since, as discussed for the 0.8 GeV data [11], there is the possibility here of some contribution from the pp → {pp} s γ reaction. In all cases the peak position was consistent with m π 0 to within the experimental uncertainties of about ±10 MeV/c 2 . Events within ±2σ of the central value were retained for the determination of the pp → ppπ 0 differential cross section. The numbers of π 0 events deduced in this way are given for the different energies in Table 1.
The luminosities recorded in Table 1 were obtained by measuring in parallel elastic proton-proton scattering using the SP trigger. The ANKE setup detects the fast proton produced by this reaction for cms angles between about 10 • and 30 • . However, to avoid regions where the acceptance changes rapidly with angle, only data from the range 15 • < θ cm p < 24 • were retained for 1.4 GeV and below while at 1.97 GeV the interval 17 • < θ cm p < 27 • was selected. After corrections for acceptance, the numbers of events in several bins of θ cm p were compared with the values of the elastic differential cross section taken from the SP07 solution provided by the SAID phase shift analysis [17]. The precision of these predictions was checked by looking at experimental data at small angles from which it was assessed to be typically about ±4%, though a little larger  Table 1 Summary of experimental conditions: L int is the integrated luminosity at the beam energy T p and N bg /N π 0 is the fractional background under the peak, where N π 0 is the number of pp → {pp} s π 0 events registered in the angular region θ cm pp < 15 • .
A full Monte Carlo simulation of the ANKE spectrometer has been developed within the framework of the GEANT program [18]. This allowed us to estimate the acceptance factors for different bins of the kinematic variables and hence, on the basis of the luminosities given in Table 1, to evaluate differential cross sections.
The experimental distributions in the final proton-proton variables for our data at all energies are very similar to those published at 0.8 GeV [11]. The excitation energy E pp spectrum is well described in terms of an S-wave pp final state interaction [4], provided that the Coulomb force is included. For E pp < 3 MeV the angular distribution of the pp relative momentum vector in the final pp rest frame is consistent with isotropy, as expected for the production of a 1 S 0 pair.
Although the coverage in the diproton angle θ cm pp with respect to the beam direction is rather limited in the ANKE spectrometer, we already noted that the data at 0.8 GeV showed a strong dependence on this angle [11]. The same is true for the other energies shown in Fig. 2 where, since the two initial protons are identical, the data are plotted as functions of cos 2 θ cm pp . For energies of 1.4 GeV and below, the results show a forward dip whereas at the highest energy the cross section is maximal in the forward direction. This is perhaps an indication that the reaction mechanism changes with energy as different intermediate nucleon isobars are excited.
and the values obtained for the forward differential cross section σ 0 and the slope parameter σ 1 are recorded in Table 2. Figure 3 summarises the results from this and a previous experiment [9] by showing the energy dependence of the forward differential cross section and the slope parameter σ 1 . Although the statistical error on the 0.508 GeV point is very large, the pp → {pp} s π 0 cross section data suggest a significant rise from this energy to a maximum in the 0.6-0.7 GeV region, where ∆(1232) production is expected to be particularly strong, followed by a monotonic decline. However, the 2 GeV point lies much above this trend and this, together with the slope information shown in the lower panel, indicates that the data have entered here a different domain. A similar tendency is seen in the pp → dπ + forward cross sections, also shown on the figure, where the cross section displays a broad minimum around 1.4 GeV after leaving the region influenced by the ∆(1232) isobar. However, π + production rises at lower energies due to contributions from the excitation of the ∆(1232 Energy dependence of the forward differential cross section for the pp → {pp} s π 0 reaction with E pp < 3 MeV. The closed circles represent the results from the present experiment while the triangles show the low energy CELSIUS data [9]. For comparison we show also the corresponding cross section for the pp → dπ + reaction. For energies up to 1.4 GeV this is represented by the dashed line taken from the SAID parameterisation [20] whereas at higher energies the Akemoto et al. data are shown as open circles [21]. Lower panel: Slope parameter σ 1 of the pp → {pp} s π 0 data, as defined by Eq. (1). For both panels the solid curve corresponds to the predictions from the N N/∆(1232)N model of Niskanen [19].
wave, which is excluded on spin-parity grounds in the π 0 case.
The behaviour of the slope parameter σ 1 shown in the lower panel of Fig. 3 is rather different, with a much smoother variation and a change of sign between 1.40 and 1.97 GeV. It is interesting to note that data on the pp → dπ + differential cross section also show a small forward dip in the 0.6-1.4 GeV energy range [20] and, furthermore, that the sign of the forward slope changes in the  Table 2 The forward differential cross section σ 0 and slope parameter σ 1 of Eq. (1) for the pp → {pp} s π 0 reaction with E pp < 3 MeV. The results at 0.425 GeV and below are from Ref. [9]. Also shown are the forward cross sections (σ π + ) for pp → dπ + [20,21] and the ratio R(π 0 /π + ) of the forward pp → {pp} s π 0 to pp → {pn} t π + cross sections integrated up to 3 MeV excitation energy.
The only theoretical estimates of the cross section and slope parameter for the pp → {pp} s π 0 reaction for the small E pp kinematical domain have been made by Niskanen [19]. The cross section predictions shown in Fig. 3, which reflect the delicate interferences between contributions from NN and ∆(1232)N intermediate states, do illustrate a displacement of the ∆(1232)N peak upwards by about 100 MeV compared to pp → dπ + . Although this agrees with our observations, Niskanen's predicted minimum around 0.7 GeV is in complete contradiction to our findings of a maximum in this region. Furthermore we see no sign of any big peak in the slope σ 1 around 0.6 GeV so that any agreement with our 0.8 GeV data seems fortuitous.
To facilitate the comparison of π + and π 0 production, the values of both the pp → {pp} s π 0 and pp → dπ + forward cross sections are given in Table 2. However, to compare the strengths of pion production leading to the S-wave spin-triplet and singlet NN states, we require rather data on the pp → {pn} t π + channel integrated up to 3 MeV excitation energy. Using finalstate-interaction theory, this can be approximated for low pn excitation energies by the pp → dπ + data through [23] dσ dΩ where k is the relative momentum in the final NN system, k 2 max /m N = 3 MeV, and α 2 t = m N B, with B being the deuteron binding energy. Though there are some deviations from this, due in part to the pn tensor force, the approximation gives a plausible description of experimental data [7]. A common feature of the R(π 0 /π + ) ratio seen in Fig. 4 is its relatively low value. Below 1.4 GeV it is typically a few per cent, which is significantly smaller than the trivial spin-statistics factor of 1/3. These conclusions are not changed significantly if the difference between the spin-singlet and triplet final state interactions are taken into account [24].
The dip in the ratio in the 0.5 GeV region is in part a reflection of the fact that the S-wave ∆(1232)N intermediate state is forbidden for π 0 production whereas it plays a vital role in the case of the π + . However it is not at all clear why R(π 0 /π + ) is small in the 1 GeV region where the contribution from the ∆(1232) maximum is much reduced. On the other hand, the much larger value at 1.97 GeV is yet another indication that the reaction mechanisms might be different at higher energies.
In summary, we have measured the differential cross section for the pp → {pp} s π 0 reaction at seven beam energies from 0.5 to 1.97 GeV under the specific kinematic conditions where the proton-proton excitation energy is below 3 MeV and the cms angle between the diproton momentum and the beam axis is less than 18 0 . The observed form of the E pp spectra and isotropy of the angular distribution in the diproton cms are consistent with the assumption that the two final protons are in the 1 S 0 state. Except for the highest energy, the data all show a sharp minimum in the forward direction and the ratio there of the differential cross sections for the pp → {pp} s π 0 and pp → dπ + reactions is below 1%. The situation changes radically at 2 GeV when a forward maximum is observed with a much enhanced value of R(π 0 /π + ). This must be a reflection of a different reaction dynamics at high energies. A broad peak around 0.6-0.7 GeV is observed in the energy dependence of the zero-degree differential cross section. The data are consistent with the predicted displacement of the ∆(1232) maximum to higher energies for π 0 production but the other features of our results are not explained by a microscopic model [19]. This could, for example, be due to phase differences between intermediate NN and ∆(1232)N contributions and much theoretical work still remains to be done.
We have attempted to study the relative amplitudes for the production of spin-singlet and spin-triplet final NN states by using a simplistic model that links the pp → dπ + and pp → {pn} t π + cross sections. The resulting R(π 0 /π + ) ratio shows that singlet production in the forward direction remains small at energies even above the ∆(1232) excitation region.
Further data are clearly required in the 0.5-0.6 GeV region where the Niskanen model should be more reliable than at higher energies to see if there is indeed any trace of his predicted structure there. The transition region above 1.4 GeV is also of interest. Does the spin-singlet suppression continue when higher isobars enter?
For a 1 S 0 final state there are only two spin-amplitudes for the pp → {pp} s π 0 reaction. Since we observe a very strong angular variation over our acceptance region, it is then likely that the analysing power will also be significant there. It is therefore planned that in the near future measurements of A y will be carried out at COSY [25] and eventually spin-correlation studies might be achievable [26]. This ensemble of data should provide a significant challenge for theory in this, one of the simplest pion-production reactions.
for their help with this experiment and to the COSY crew for providing such good working conditions. J.A. Niskanen kindly provided the numerical values of the calculations of Ref. [19].