Differential cross section and analysing power of the pp ->{pp\}_s pi0 reaction at 353 MeV

In order to establish links between p-wave pion production in nucleon-nucleon collisions and low energy three-nucleon scattering, an extensive programme of experiments on pion production is currently underway at COSY-ANKE. The final proton pair is detected at very low excitation energy, leading to an S-wave diproton, denoted here as {pp}_s. We now report on measurements of the differential cross section and analysing power of the pol{p}p->{pp}_s pi^0$ reaction at 353 MeV. Both observables can be described in terms of s- and d-wave pion production and, by using the phase information from elastic pp scattering, unique solutions can be obtained for the corresponding amplitudes. This information is vital for the partial wave decomposition of the corresponding pn->{pp}_s pi^- reaction and hence for the extraction of the p-wave terms.

Within the context of chiral perturbation theory, a significant step forward in our understanding of pion physics at low energies would be to establish that the same shortranged N N → N N π vertex contributes to p-wave pion production, to low energy three-nucleon scattering [1][2][3], γd → πN N [4,5] and πd → γN N [6], as well as in weak reactions like tritium beta decay [7][8][9][10]. The relevant transition amplitude, which connects N N S-waves in the initial and final state with a p-wave pion, contributes to both pp → π + d(π + pn) and pn → ppπ − . However, the extensive data for π + production is of limited use in this context, be- * Corresponding author.
cause the p-wave amplitudes are completely dominated by the 1 D 2 initial state, which hinders a reliable extraction of the 1 S 0 initial state [3].
There is a programme at the COSY-ANKE facility of the Forschungszentrum Jülich to perform a complete set of measurements on N N → {pp} s π at low energy [11]. Here the {pp} s denotes a proton-proton system with very low excitation energy, E pp . At ANKE we select events with E pp < 3 MeV and, under these conditions, the diproton is overwhelmingly in the 1 S 0 state with antiparallel proton spins. This simplifies enormously the spin structure: a partial wave analysis for pp → ppπ 0 without the E pp cut would require twelve additional P -wave final pp spin-triplet states [12][13][14]. The cut also allows one to extract the full information on the production amplitudes without having to make measurements of the final proton polarisations.
Whereas the {pp} s final state is isotriplet, the isosinglet np initial state also contributes to pn → ppπ − . In order to isolate this, which contains the amplitudes of interest, the isotriplet channel needs to be well understood. As the first part of the outlined larger programme, we therefore report here on measurements of the cross section and proton analysing power in the pp → {pp} s π 0 reaction at T p = 353 MeV.
For a spin-singlet diproton, the spin structure of the pp → {pp} s π 0 or np → {pp} s π − reaction is that of Parity and angular momentum conservation require that the initial nucleon-nucleon pair to have spin S = 1. The pion orbital angular momentum ℓ and the initial nucleon-nucleon isospin I are then linked by ℓ + I = odd so that, for the pp → {pp} s π 0 reaction, only even pion partial waves are allowed. As a consequence, the unpolarised cross section for π 0 production, and this times the proton analysing power A y , must be of the form where θ π is the pion c.m. production angle with respect to the direction of the polarised proton beam. Here p is the incident c.m. momentum and k that of the produced pion which, at 353 MeV, have values p = 407 MeV/c and k ≈ 94 MeV/c, where the latter represents an average over the 3 MeV E pp range. The only detailed measurements of the pp → {pp} s π 0 differential cross section over the whole angular range were carried out with the PROMICE-WASA apparatus at CEL-SIUS at a series of energies from 310 to 450 MeV, using the same standard 3 MeV cut on E pp [15]. Throughout this energy range, significant anisotropies were found in the angular distributions which were attributed to interferences between pion s and d waves. On the other hand, there were no corresponding measurements of the proton analysing power, which might also be driven by a strong s-d interference.
We have previously reported measurements of the pp → {pp} s π 0 differential cross section at several energies and small angles [16,17]. Since these were carried out using the ANKE spectrometer [18] under conditions that were similar to the current ones, the description here can be quite brief. ANKE is placed at an internal beam station of the COSY cooler synchrotron. Fast charged particles, resulting from the interaction of the stored transversally polarised proton beam with the hydrogen cluster-jet target [19] and passing through the analysing magnetic field, were recorded in the forward detector (FD) system. The FD, which was the only detector used in this experiment, includes multiwire proportional chambers for tracking and a scintillation counter hodoscope for energy loss and timing measurements.
To start the identification of the pp → {pp} s π 0 reaction, proton pairs were first selected from all the registered twotrack events using the measured momenta of the both particles and the difference in their time-of-flight [20]. The resolution σ(E pp ) in the diproton excitation energy was better than 0.6 MeV, which allowed the E pp < 3 MeV cut to be applied reliably. After selecting the 1 S 0 final state, the kinematics of the pp → {pp} s X process could be reconstructed on an eventby-event basis to obtain a missing-mass M X spectrum. A two-dimensional distribution of M 2 X versus the c.m. polar angle of the diproton θ cm pp is presented in Fig. 1. This demonstrates the large angular acceptance of the apparatus for the pp → {pp} s π 0 reaction at 353 MeV and shows a clean π 0 signal with an almost negligible background. Simulations indicate that the c.m. angular resolution is better than 5 • .
The polarization asymmetry is defined by where N ↑ and N ↓ are the numbers of pp → {pp} s π 0 events with beam proton spin up and down, corrected for dead time, and L ↑ and L ↓ are the corresponding luminosities. The relative luminosity L ↑ /L ↓ ≈ 0.985 ± 0.015 was estimated using events at very small polar angles, where the polarization asymmetry should be negligible. This procedure adds about a 3% systematic error to the values of ε.
The analysing power A y is connected to the asymmetry through: where P is the transverse polarization of the beam and cos φ pp the average over the diproton azimuthal angular distribution. Since the cos φ pp acceptance is concentrated near 1, all the events in the regions analysed contribute usefully to the A y measurement. The polarization of the proton beam was flipped between "spin-up" to "spin-down" (perpendicular to the plane of the the pp → dπ + reaction. The predictions of the SAID program [21] have been scaled to agree with the experimental data and these give average COSY proton beam polarisations of (a) P = 0.687 ± 0.008 and (b) P = 0.668 ± 0.016. In neither case was the uncertainty in the SAID prediction included.
accelerator) every six minutes and no measurements were made with an unpolarized beam. The value of P was estimated from proton-proton elastic scattering and the pp → dπ + reaction that were measured in parallel. The analysing powers for these reactions were taken from the SAID analysis program, solutions SP07 for pp → pp and SP96 for pp → dπ + [21]. The results of the two methods shown in Fig. 2 agreed within measurement errors and gave an average polarization of P = 0.68 ± 0.03, where the error includes the uncertainties arising from the calibration reactions.
A simulation was undertaken of the two-dimensional acceptance in terms of the pp excitation energy E pp and its c.m. polar angle θ pp . This took into account the geometry of the setup and the sensitive areas of the detectors, the efficiency of the multiwire proportional chambers and the track reconstruction algorithm. In order to avoid potential problems arising near the limits of the acceptance, cuts were made around the edges of the exit window of the spectrometer magnet in both the experimental data and simulation. This is only a challenge at the larger angles, 80 • < θ π < 100 • , where a compromise had to be made regarding the acceptance ambiguities and this introduces an extra 4% systematic uncertainty in this angular region.
The numbers of detected π 0 events were then corrected on an event-by-event basis for acceptance, dead time and relative luminosity L ↑ /L ↓ . The latter were important be-cause, in the absence of data with an unpolarised beam, an average has to be evaluated.
The luminosity in the experiment was estimated from measurements of pp elastic scattering carried out in parallel. The numbers of detected events, corrected for the dead time, were compared with a simulation that used a generator which included the differential cross section obtained from the SAID analysis program [21]. Although this program does not furnish error bars, experimental data at nearby energies suggests that the associated uncertainty is about 2%, to which must be added 3% arising from acceptance and similar systematic effects. At this level the statistical error is negligible and the resulting total luminosity was estimated to be 544 ± 22 nb −1 . At this energy the pp → dπ + cross section data are less precise than those of pp elastic scattering but, on the basis of the SAID predictions, one obtains the completely consistent luminosity estimate of 547 nb −1 . The errors shown here are statistical together with a 4% systematic contribution in the 80 • < θπ < 100 • region coming from the acceptance ambiguity discussed in the text. The overall systematic uncertainty is about 4%. Open (red) circles are CELSIUS data obtained at 360 MeV [15]. It should be noted that the latter data represent averages of measurements taken in both hemispheres. The curve is a linear fit in cos 2 θπ to our data.
The differential cross section results are presented in Fig. 3, where they are compared to those obtained at 360 MeV at CELSIUS [15]. Within the 10% luminosity uncertainty in these data, the overall agreement is very good. However, the CELSIUS data at this energy level off a little around 90 • . This seems to be a feature only of the 360 MeV results since, at the other energies, linear fits in cos 2 θ π all have good values of χ 2 /NDF [15].
Fitting our data with a polynomial in cos 2 θ π , as in Eq. (1), gives parameters a 0 = 4.05 ± 0.08 µb/sr, a 2 = −2.31 ± 0.14 µb/sr, Apart from the acceptance uncertainties at the larger angles, the error bars quoted here are purely statisti-cal; the ±4% systematic uncertainty from the luminosity and acceptance largely cancels in the ratio a 2 /a 0 . Since χ 2 /NDF = 23/20, there is clearly no compelling evidence for any cos 4 θ π dependence, i.e., a non-zero a 4 coefficient, and this possibility has been omitted from the curve in Fig. 3.
[nb/sr] The results for the analysing power of the pp → {pp} s π 0 reaction are displayed in Fig. 4, with A y (dσ/dΩ) being shown in panel a and A y in panel b. These observables must be antisymmetric about 90 • and the crossing of the data through zero around this angle is some confirmation of our estimation of L ↑ /L ↓ . These data are subject to the overall uncertainties associated with the luminosity and acceptance evaluation, though these are not relevant for the A y in Fig. 4b. There remains, however, the ±3% arising from the uncertainty in the value of L ↑ /L ↓ .
The resulting fits are shown in Fig. 4. In order to understand the significance of the results reported here, we must attempt a partial wave description of the data. The most general form of the reaction amplitude is where S is the polarisation vector of the initial pp spintriplet state.p andk are unit vectors in the c.m. frame along the directions of the incident proton and final pion, respectively. The observables studied here are expressed in terms of the two scalar amplitudes A and B through [14] dσ The experimental data show no evidence for high partial waves at 353 MeV and so we model these results with only ℓ = 0 and ℓ = 2 contributions. The latter can arise from initial L = 1 or L = 3 waves so that, in total, there are three possible transitions, 3 P 0 → 1 S 0 s, 3 P 2 → 1 S 0 d, and 3 F 2 → 1 S 0 d, see e.g. Ref. [22] for the explicit form of the spin-angular structures. We denote the corresponding amplitudes by M P s , M P d , and M F d , respectively. Expanding the scalar amplitudes in terms of these partial waves gives Equations (8) and (9) then allow one to relate the measured observables of Eqs. (1) and (2) to the partial wave amplitudes. For consistency, since we have neglected any possible effects arising from s-g interference, we shall also drop terms that are bilinear in d-wave production amplitudes. In this approximation and so the data only provide three relations between the three complex amplitudes. The transverse spin correlation parameters contain no extra information since A y,y = 1 and this is also true for A x,x up to d-d interference terms. If the longitudinal-transverse spin correlation parameter A x,z were measured, this would provide one further relation but this would still not be sufficient for an unambiguous partial wave decomposition. For this we need information about the phases of the production amplitudes. The 3 P 0 partial wave is uncoupled and, at the energy where the experiment was performed, its inelasticity is very small. Under these conditions the Watson theorem, which fixes the phase induced by the initial state interaction to that of the elastic proton-proton scattering, applies [23]. Thus we take M P s = |M P s |e iδ3 P 0 , with δ3 P0 = −14.8 • [21]. Note that we do not include any phase associated with the 1 S 0 final pp state because it is common for all partial waves and therefore does not affect the observables.
For coupled channels, such as 3 P 2 − 3 F 2 , the strict conditions of the Watson theorem do not apply. However, at our energy the mixing parameter, as well as the inelasticities, are still negligibly small and thus to a good approximation we may also use the Watson theorem here. Further evidence for the smallness of the channel coupling is to be found in two potential models [24,25]. In both models the T -matrix for the transition from the 3 F 2 to the 3 P 2 wave is almost real; the phase of M P d is driven by δ3 P2 = 17.9 • , whereas the phase of M F d can be neglected. The quality of this approximation was also checked by explicit calculations of the d-wave production amplitudes within chiral effective field theory. These showed that up to order m π /m N the above phase assumptions should be valid to within ±2 • .
Using the phase information in this way, we find that The values quoted here were obtained by considering also our np → {pp} s π − data though the numbers would change but marginally if one included only the pp → {pp} s π 0 results in the fit. The error bars quoted here are statistical and do not include the overall systematic uncertainties. However, changing the normalisations of the differential cross section and analysing powers by 3% and 4%, respectively, leads to changes that are comparable to the quoted errors.
On the other hand, we could not investigate the less tangible ones associated with the neglect of the channel coupling and the truncation in the partial wave expansion. The weakness of pion production from the initial 3 F 2 waves at 353 MeV, in addition to being in agreement with theoretical prejudices, is also consistent with the low inelasticity found for this wave [21].
In summary, we have measured the differential cross section and analysing power of the pp → {pp} s π 0 reaction at 353 MeV. The angular distributions of A y and dσ/dΩ are both well represented by retaining only pion s and d waves in a phenomenological description. The values of dσ/dΩ agree well with the results obtained at CELSIUS [15] over most of the angular range. However, at this energy these data flatten off around the middle of the angular distribution and, if this effect were correct, it would signal a large contribution from 3 F 2 or even higher partial waves.
By making plausible assumptions on the coupling between the nucleon-nucleon channels and invoking the Watson theorem it was possible to estimate the partial wave amplitudes with their phases. These could be checked through a future measurement of the spin-correlation parameter A x,z , though this would require the installation of a Siberian snake in COSY.
In an associated letter [26], the isospin-0 amplitudes are investigated through the measurements of the differential cross section and analysing power of the quasi-free pn → {pp} s π − reaction in this energy domain. The extraction of the p-wave amplitudes from the data of Ref. [26] required a knowledge of the s-and d-wave amplitudes of the type provided here. In addition, data have already been taken on the transverse spin correlation parameter A x,x for this reaction [27]. The full collection of these results will lead to very useful constraints on the parameters of the chiral effective field theory that link pion production to the threenucleon interaction [1,2].