title: A comprehensive analysis of differential cross sections and analyzing powers in the proton-deuteron break-up channel at 135 MeV
authors: H. Tavakoli-Zaniani M. Eslami-Kalantari H. R. Amir-Ahmadi M. T. Bayat A. Deltuva J. Golak N. Kalantar-Nayestanaki St. Kistryn A. Kozela H. Mardanpour J. G. Messchendorp M. Mohammadi-Dadkan A. Ramazani-Moghaddam-Arani R. Ramazani-Sharifabadi R. Skibinski E. Stephan H. Witala
text: 
  
  
  
  A comprehensive analysis of differential cross sections and analyzing powers
  
  
  
  
  
                      
  
           
  
  KVI-CART, University of Groningen, Groningen, the Netherlands Department of Physics, School of Science, Yazd University, Yazd, Iran Institute of Theoretical Physics and Astronomy, Vilnius University, Lithuania M.Smoluchowski Institute of Physics, Jagiellonian University, Kraków, Poland Institute of Physics, Jagiellonian University, Kraków, Poland Institute of Nuclear Physics, PAS, Kraków, Poland Department of Physics, University of Sistan and Baluchestan, Zahedan, Iran Department of Physics, Faculty of Science, University of Kashan, Kashan, Iran Department of Physics, University of Tehran, Tehran, Iran Institute of Physics, University of Silesia, Chorzów, Poland
  
  
   
  
   
  
  
  
  
  
  
  A selection of measured cross sections and vector analyzing powers, A_x and A_y, are presented for the pd break-up reaction. The data are taken with a polarized proton beam energy of 135 MeV using the Big Instrument for Nuclear-polarization Analysis (BINA) at KVI, the Netherlands. With this setup, A_x is extracted for the first time for a large range of energies as well as polar and azimuthal angles of the two outgoing protons. For most of the configurations, the results at small and large relative azimuthal angles differ in behavior when comparing experimental data with the theoretical calculations. We also performed a more global comparison of our data with theoretical calculations using a chi-square (χ^2) analysis. The cross-section results show huge values of χ^2/d.o.f.. The absolute values of χ^2/d.o.f. for the components of vector analyzing powers, A_x and A_y, are smaller than the ones for the cross section, partly due to larger uncertainties for these observables. However, also for these observables no satisfactory agreement is found for all angular combinations. This implies that the present models of a three-nucleon force are not able to provide a satisfactory description of experimental data. 
  
  
    
        21.30.-xNuclear forces   21.45.+vFew-body systems  24.70.+sPolarization phenomena in reactions  25.45.DeElastic and inelastic scattering  
        
   
  
  A comprehensive analysis of differential cross sections and analyzing powers in the proton-deuteron break-up channel at 135 MeVSupplementary material in the form of a pdf file available from the journal web page at http://***** 
      H. Tavakoli-Zaniani1,2h.tavakoli.zaniani@rug.nl M. Eslami-Kalantari2meslami@yazd.ac.irH. R. Amir-Ahmadi1M. T. Bayat1A. Deltuva3J. Golak4N. Kalantar-Nayestanaki1St. Kistryn5A. Kozela6H. Mardanpour1J. G. Messchendorp1j.g.messchendorp@rug.nlM. Mohammadi-Dadkan1,7A. Ramazani-Moghaddam-Arani8R. Ramazani-Sharifabadi1,9R. Skibiński4E. Stephan10H. Witała4
  
      Received: date / Revised version: date
  ==================================================================================================================================================================================================================================================================================================================================================================
  
  
  
  
  
  
  
  § INTRODUCTION
  
  
  
  
  Although the nucleon-nucleon (2N) interaction has been studied extensively in the past using proton-proton and proton-neutron scattering data, the role of higher-order forces, such as the three-nucleon force (3NF) remains mysterious. The need for an additional three-nucleon potential became evident when comparing three-body scattering observables <cit.> and binding energies of light nuclei with state-of-the-art calculations <cit.>. The two nucleon force models such as CD Bonn, Argonn V18, Reid93, Nijmegen I and Nij­me­gen II <cit.> are able to describe the two nucleon systems very well below the pion-production threshold. The next step would be to significantly extend the world database in the three-nucleon scattering system as a benchmark to eventually have a better understanding of the structure of the three-nucleon interaction. For almost all observables in nucleon-deuteron elastic scattering, the calculations which only include two nucleon forces (2NFs) fail to a large extent to describe the data, in particular at energies above 60 MeV and at large center-of-mass scattering angle. In addition to the elastic channel, the deuteron break-up reaction offers rich spectrum of kinematical configurations and as such provides a good testing ground for understanding the structure of the nuclear force <cit.>. Several theoretical approaches have been developed, such as a dynamic Δ-isobar <cit.> and the Tucson-Melbourne <cit.> 3NFs and these have been embedded within rigorous calculations using the Faddeev-type equations by, for example, Bochum-Kraków <cit.> and Hannover-Lisbon <cit.> groups. Besides these phenomenological approaches, also two- and three-nucleon for­ces 3NF have been constructed from chiral perturbation theory (ChPT). The leading 3NF in ChPT shows significant contributions to the nuclear force <cit.>, but the most advanced nowadays complete and consist chiral calculations at the third order of chiral expansion <cit.> deliver the 3NF data description of the similar quality to the one from semi-phenomenological models. 
  
  The proton-deuteron break-up is a suitable reaction to study three nucleon systems one can measure various observables in a large part of the available phase space of this reaction. In this paper, the cross sections and vector analyzing powers, A_x and A_y, for d(p, pp)n reaction at 135 MeV are extracted from configurations where the two final-state protons scatter at small polar angles between 14^∘-30^∘ . The data taken at other scattering angles have been reported in Ref. <cit.>.
  
  
  
  § EXPERIMENTAL SETUP
  
  
  The pd break-up reaction was studied using a polarized proton beam of 135 MeV impinging on a liquid deuterium target which was located at the center of BINA (Big Instrument for Nuclear-polarization Analysis).
  The po­la­rized beam is provided with POLIS (POLarized Ion Sou­rce) <cit.>. The beams of (polarized) protons and deuterons are accelerated by AGOR (Accelerateur Groningen ORsay) <cit.> at KVI, the Netherlands. The proton-deuteron break-up reaction was studied with BINA. The BINA detector is particularly suited to study the elastic and break-up reactions at intermediate energies. BINA is composed of two major parts, the forward-wall and the backward-ball. 
  The forward-wall measures the energy and scattering angles of final-state particles in the range 10^∘-37^∘. The forward-wall is composed of three main parts, Energy scintillators (E-scintillators), Δ E-scintillators, and a Multi-Wire Proportional Chamber (MWPC).
  The backward-ball is made of 149 small cut pyramid-shaped scintillator detectors by a ball-shaped detector which covers the rest of the polar angles up to 165^∘. 
  Therefore, the BINA detector covers almost the complete phase space of the break-up and elastic reactions. For a more detailed description of the detector, we refer to <cit.>. In this work, we present the results of break-up configurations in which the final-state protons are registered in coincidence by the forward-wall.
  
  
  
  § DATA ANALYSIS
  
  
  The data analysis of pd break-up reaction, taken with a proton-beam energy of 135 MeV, was performed with the goal of measuring the vector analyzing powers, A_x and A_y, and the differential cross sections.
  
  Events of the break-up reaction are identified by reconstructing the scattering angles and energies of the two final-state protons. During data taking, a hardware trigger was used requiring at least two of the ten E-scintillators to give a signal above the threshold (∼1 MeV). These events were further processed offline by combining the information of the MWPC with the corresponding E-scintillators. In this way, two proton candidate tracks were re­cons­truc­ted for further analysis. The E-scintillators were calibrated by matching their raw charge-to-digital converter (QDC) information with the expected energy correlation of break-up events. Details of the analysis can be found in Ref. <cit.>.
  
  The energy correlation between the two outgoing protons, E_2 versus E_1, after the calibration for a particular configuration (θ_1, θ_2, ϕ_12)=(24^∘±2^∘, 24^∘±2^∘, 180^∘±5^∘) is shown in the left panel of Fig. <ref>, whereby θ_1 and θ_2 are the polar angles of two outgoing protons and ϕ_12 is their relative azimuthal opening angle.
   The solid line shows the kinematical S-curve calculated for the central values of the angular bins. The kinematic variable S corresponds to the arc-length along the kinematic curve with S=0 at the point where E_1 is at its minimum. To measure the break-up observables, at the first step, we make several slices along the kinematical S-curve with a window of ∼9.5 MeV. We note that the energy resolution, ∼4 MeV, is signficantly smaller than this window size. The projection of the indicated region on the line perpendicular to the S-curve (D-axis) is shown in the right panel of Fig. <ref>. 
  The peak around zero corresponds to break-up events. Most of the events on the left-hand side of the peak are also due to break-up events. In these cases, the protons have lost energy due to hadronic interactions inside the detector. The amount of accidental background is small as can been seen from the small amount of events on the right-hand side of the peak. We fit this spectrum by using a third-order polynomial, representing the hadronic interactions and the accidental background, and a Gaussian function, representing the signal. The extracted number of signal events was corrected by the data-acquisition dead-time and the down-scaling factor. This number is subsequently used to measure  the cross sections and vector analyzing powers.   
   
   The cross section of the break-up reaction can be obtained by: 
  
  
      d^5σ/dΩ_1dΩ_2dS=N/QtϵΔΩ_1ΔΩ_2Δ S,
  
  where N is, the number of break-up events in each slice along the S-curve corrected for the down-scaling factor and the dead-time, Q is the total integrated charge, t is the number of the scattering centers, ϵ is the multiplication of all the efficiencies including the MWPC efficiency, hadronic correction and geometrical efficiency, ΔΩs are the solid angles for the two outgoing protons and Δ S is, the width of the selected window in each slice along S-curve <cit.>. 
  We studied various sources that we identified as the main contributors to the systematic uncertainty in the cross section measurements. In the following, we briefly summarize each of them and we give a description on how magnitudes of corresponding errors have been estimated.
  
  The first source we identified as a contributor to the systematic error is related to uncertainties in the determination of the effective target thickness. Taking into account the bulging of the target, we estimated an effective target thickness of 3.85±0.20 mm. The resulting error in this measurement (5%) is assigned as a systematic error in the cross section measurements. This value has been estimated by earlier cross section studies of the elastic proton-proton scattering process using similar targets by comparing data with precision calculations of this reaction <cit.>.
   
  The second systematic uncertainty that we considered is related to the error in estimating the fraction of events that suffered from a hadronic interaction in the scintillators of BINA. Since in the calculation of the number of break-up events, we only account for those events for which the energy of both protons are well reconstructed, one needs to correct for the hadronic interaction effect. To determine this effect, we used Monte Carlo studies that are based on the interaction models provided by the GEANT-3 simulation package <cit.>. Typically, we found that about 12% of all break-up events suffered from hadronic interactions. The uncertainty of this value (6%) is assigned as a source of systematic uncertainty. It has been estimated by taking the difference between the number of hadronic background events derived from simulations with the value estimated from a fit of the measured D spectrum (right panel of Fig. <ref>) <cit.>. 
  
  The third source of systematic uncertainty is associated with the trigger efficiency. This efficiency has been studied using Monte Carlo simulations based on GEANT-3. It was found that for break-up events whereby ϕ_12 is larger than 20^∘, the trigger efficiency is about 98% and that it drops to 88% for selected events associated with ϕ_12=20^∘. To be conservative, we assigned a systematic error due to the trigger efficiency by taking the observed inefficiencies using the Monte Carlo results, therefore 2% for ϕ_12>20^∘ and 12% for ϕ_12=20^∘ <cit.>. 
  
  The fourth source of systematic error is due to uncertainties in the efficiency determination of the MWPC. Proton tracks from the elastic proton-deuteron scattering process were identified using the information of the E and ΔE detectors. The E-ΔE hodoscope provides a grid that is used to map onto the MWPC. This allows us to measure the MWPC efficiency for protons at various locations corresponding to every E-ΔE hodoscopy. Typically, we found an efficiency of about (92±1)% for each proton, whereby the error corresponds to statistical fluctuations of the unbiased data sample that is used in this study. We associated a systematic error due to uncertainties of the MWPC efficiency for the cross section measurements by summing up the efficiency errors of the two final-state protons, i.e. 2% <cit.>.
  
  The total systematic uncertainties for the cross sections at small relative azimuthal angles (≤20^∘) are about 14% and for the larger relative azimuthal angles (> 20^∘) is about 9%. For this, we added up, quadratically, the systematic errors of the various sources assuming them to be independent.
  
  To measure the vector analyzing powers, the number of break-up events were normalized to the collected beam charge for the two polarization states (up and down). The relation between the normalized number of events with the polarized beam, N_ξ,ϕ_12^s(ϕ), and unpolarized beam, N_ξ,ϕ_12^0, is given by <cit.>:
  
  
  
      N_ξ ,ϕ_12^s(ϕ)=      N_ξ ,ϕ_12^0(1+p_z^sA_y(ξ ,ϕ_12)cosϕ
            -p_z^sA_x(ξ ,ϕ_12)sinϕ),
  
  where s indicates the spin of the beam and ξ defines a given kinematical point (θ_1, θ_2, S). The component of the vector polarization of the beam is given by p_z and the vector analyzing powers are indicated by A_x and A_y. Here, ϕ is the angle between
  quantization axis for the polarization and the normal to the scattering plane of the first nucleon in the laboratory frame of reference, with ϕ_1=0. 
  Since the statistics obtained with an unpolarized beam was limited, we extracted the spin observables by solely using N_ξ, ϕ_12^↑(ϕ) and N_ξ, ϕ_12^↓(ϕ), corresponding to the normalized number of events for the spin-up and spin-down polarized beams, respectively. The analyzing powers A_x and A_y are extracted using the following relation:
  
  
      f_ξ ,ϕ_12(ϕ)     = N_ξ ,ϕ_12^↑(ϕ)-N_ξ ,ϕ_12^↓(ϕ)/N_ξ ,ϕ_12^↑(ϕ)p_z^↓-N_ξ ,ϕ_12^↓(ϕ)p_z^↑
          =A_y(ξ ,ϕ_12)cosϕ-A_x(ξ ,ϕ_12)sinϕ,
  
  where p_z^↑ and p_z^↓ are the values of up (0.57±0.03) and down (-0.70±0.04) beam polarizations. The polarization of the proton beam is defined as:
  
  
      p_z = N^+-N^-/N^++N^-,
  
  where N^+,- are the number of particles with a particular spin (up or down). The beam polarization has been determined using the in-Beam Polarimeter (IBP) <cit.> that was installed at the high-energy beam at KVI. The IBP measured regularly the beam polarization by recording the azimuthal asymmetries of the H(p⃗, pp) reaction. The vector analyzing power of the proton-proton scattering process was used as input to the polarization measurements and its uncertainty is the main source of error.
  Parity conservation imposes the following restrictions on the components of the vector analyzing powers <cit.>:
  
  
      A_x(ξ ,-ϕ_12)= -A_x(ξ ,ϕ_12); 
      
       A_y(ξ ,-ϕ_12)= A_y(ξ ,ϕ_12),
  
  where for ϕ_12=180^∘, we expect A_x=0. By taking the sum and difference of f_ξ, ϕ_12(ϕ) and f_ξ,-ϕ_12(ϕ) in combination with the results of Eq. <ref>, the following combination of asymmetries for mirror configurations (ξ,ϕ_12) and (ξ,-ϕ_12) can be obtained <cit.>:
  
  
      g_ξ ,ϕ_12(ϕ)    = f_ξ ,ϕ_12(ϕ)+f_ξ ,-ϕ_12(ϕ)/2,  
         =A_y(ξ ,ϕ_12)cosϕ;
      
       h_ξ ,ϕ_12(ϕ)    = f_ξ ,ϕ_12(ϕ)-f_ξ ,-ϕ_12(ϕ)/2,  
         =-A_x(ξ ,ϕ_12)sinϕ.
  
  The components of vector analyzing-power values, A_x and A_y, are obtained from the fits of Eq. <ref> for various kinematical configurations. 
  
  
  
  
   
  
  
  
  The error of the beam polarization is about 6%. For instance, the beam polarization for the down-mode
  has been measured at a value of ∼0.70±0.04, which gives rise to 6% systematic uncertainty in the beam polarization. We estimated the impact of the polarization uncertainty on the analyzing powers by recalculating both analyzing powers with an input polarization that differs by +6% (-6%) for the spin-up (down) mode. The difference with the results using the nominal values of the beam polarizations is used as an estimate of the corresponding systematic error. We also considered a systematic error due to asymmetries that are induced by rate- or polarization-dependent differences in detection efficiencies that do not cancel in Eq. <ref>. This uncertainty has been estimated by exploiting data at particular kinematical configurations for which the vector analyzing powers are known or constrained. For A_y, we have analyzed various symmetric configurations for which both protons scatter to the same polar angle with a relative azimuthal angle of 180^∘. By taking the average vector analyzing power for the covered S-range, one expects a value of zero. We have performed a fit with a free offset value to the data and we used the corresponding offset as a measure of the systematic uncertainty for A_y. To estimate the systematic error for A_x, we analyzed the data for a relative azimuthal angle of 180^∘ for which A_x should be zero, and extracted the corresponding value for A_x as a function of S. Subsequently, these results are fitted with a zeroth-order polynomial and its value is used as an estimate for the corresponding systematic error. This error and the error in the polarization are added in quadrature assuming them to be independent, to form the total systematic uncertainty.
  
  
  
  
  
  Figures <ref>–<ref> show the cross-sections and vector analyzing powers as a function of S for symmetric and asymmetric configurations at small, intermediate and large relative azimuthal angles. The results of our analysis are indicated as black dots. The error bars indicate statistical uncertainties which are in some cases smaller than the symbol sizes. The cyan bands depict the systematical uncertainty whereby the width corresponds to 2σ. The various lines present the results of Faddeev calculations using 2NF and 2N+3NF models. The results show a different behavior between the data and theoretical calculations at small and large relative azimuthal angles. 
  
  The results of the cross sections and vector analyzing powers as a function of S for about hundred configurations (with 14^∘<θ_1<30^∘, 14^∘<θ_2<30^∘ and 0^∘<ϕ_12<180^∘) for incident proton energy of 135 MeV were extracted. A small subset of vector analyzing power data for selected symmetric configurations was presented in Ref. <cit.>.
  An extensive overview of all the results can be found in the supplementary material associated with this paper <cit.>. In general, we observe similar patterns for other scattering angles comparable to the ones shown in Figs. <ref>–<ref> with respect to small and large relative azimuthal angles. Similar behaviours were also observed in the measurements of the same observables at 190 MeV <cit.>.
  
  To have a more efficient study and to compare globally the theoretical predictions with the complete data set, a global analysis is performed with averages of observables <cit.>. In this paper, we compare our data with the theoretical predictions using a chi-square analysis for several (θ_1, θ_2) and as a function of ϕ_12 as another approach to perform a global analysis. The quantity χ^2 per degree of freedom is defined by 
  
  
      χ_m^2/d.o.f.=1/N-1∑_i=1^N{O_i-T_i^m/σ_i} ^2,
  
  where N is the number of specific configuration in (S, θ_1, θ_2, ϕ_12), O_i is one of the observables (A_x, A_y, or d^5σ/dΩ_1dΩ_2dS), σ_i is the measured statistical error of a data point and T_i^m is the results of the theoretical calculation whereby the sup-index m refers to a specific model, namely, CDB, CDB+TM', CDB+Δ and CDB+Δ+Cou­lo­mb, the sub-index i refers a specific configuration in (S, θ_1, θ_2, ϕ_12) and N is the number of specific configuration in (S, θ_1, θ_2) in the right panels of Figures <ref> and <ref> or the number of specific configuration in (S, ϕ_12) in the left panels of Figures <ref> and <ref>. 
  Figures <ref> and <ref> show the results of χ^2/d.o.f. for sum over (S, ϕ_12) for specific (θ_1, θ_2) versus the angular combination (θ_1 and θ_2) (left panels) and for sum over (S, θ_1, θ_2) for specific ϕ_12 versus ϕ_12 (right panels) for different observables (A_x, A_y and d^5σ/dΩ_1dΩ_2dS) for data taken with a proton-beam energy of 135 MeV and 190 MeV <cit.>, respectively. The asymmetric error bars reflect the systematic uncertainty of the data with respect to one of the theoretical calculations, namely CDB+Δ+Coulomb. These errors were obtained by adding and subtracting the estimated total systematic error to and from the data resulting in two alternative chi-square values corresponding to the edges of the error bars. 
  
  
  
  § DISCUSSION
  
  
  As observed in Fig. <ref>, at small azimuthal opening angles the results are closer to the predictions of the theoretical approach that deploy CDB+Δ+Coulomb potential. This demonstrates that the Coulomb effect is sizeable for this observable at these configurations. Note that in this case the relative energy between the two protons is small. For large relative azimuthal angles, the model based on the CDB+TM' potential appears to be the closest to the experimental data, albeit the differences between the various models are in general small. For intermediate values of ϕ_12, the predicted shape of the cross sections differ significantly with the data. In Figs. <ref> and <ref>, the measurements of A_x for relative azimuthal angles of 180^∘ is found to be consistent with zero as expected from Eq. <ref>. This demonstrates that our procedure to extract the analyzing powers does not suffer from experimental asymmetries. This is also confirmed by our estimate of the systematic uncertainty, which is found to be small. Our polarization observables are reasonably well described by the calculations for kinematical configurations at which the three-nucleon force effect is predicted to be small. However, striking discrepancies are observed at specific configurations, in particular in cases where the relative azimuthal angle between the two outgoing protons becomes small. In this range, the measured values of A_y is close to the results of the 2NF calculation. Although, the disagreement is still significant, the effects of the Coulomb force are very small. The addition of the TM'3NF makes the agreement even worse. Therefore, the origin of this discrepancy must lie in the treatment of 3NFs. The same behavior was observed for the data taken at a beam energy of 190 MeV <cit.>. Possibly, the modeling of short-range 3NF must be significantly improved as e.g. chiral perturbation theory, for which the data presented in this paper and those in Ref. <cit.> can be used as a benchmark.
  
  The results of the analyzing powers for different combinations of (θ_1=θ_2, ϕ_12) for small ϕ_12 which corresponds to d(p,^2He)n for 135 MeV proton beam energy were also compared to the results using a proton beam with an energy of 190 MeV <cit.> to study the spin-isospin sensitivity of the 3NF models <cit.>. 
  
  By inspecting Fig. <ref>, we note that the absolute values of the χ^2/d.o.f. for the analyzing powers A_x and A_y appear to behave better than the ones for the cross section. However, also for analyzing powers, there are clear trends to be observed in which all the model predictions deviate, beyond statistical and systematic uncertainties, from the data.
  For A_y, the trend observed in the plots as a function of polar angle combination looks similar to what is observed for the cross section. The models show a larger discrepancy towards larger angles. The trends as a function of ϕ_12 are vastly different compared to the ones observed in the cross section. Although A_y features a worse agreement towards small ϕ_12 (and partly large ϕ_12), the observable A_x is well predicted at small ϕ_12 except for CDB+TM'. By comparing all the model predictions, the calculation based on the CDB+Δ+Coulomb model is the most compatible with the data.
  
  By comparing the results between the two energies, see Figs. <ref> and <ref>, in general, similar trends as a function of ϕ_12 are observed for the cross section and A_y. For the cross section data taken at 190 MeV, the calculation that is based on CDB+Δ+Coulomb potentials shows the worst agreement, in particular for large values of ϕ_12. The sensitivity to 3NF effects appears to be larger at the higher energy. For both energies, it is clear that the inclusion of the TM' 3NF is by far not sufficient to remedy the observed discrepancies. We also note that the Coulomb effect is very small for both spin observables. By globally reviewing the chi-square data, we note that CDB+TM' gives the worst description of the data for analyzing powers. 
  
  
  
  
  § SUMMARY AND CONCLUSIONS
  
  
  Finding a suitable theory of nuclear forces is one of the main challenges in nuclear physics. 
  To study the three nucleon systems, the reaction d(p, pp)n was studied at KVI using a polarized proton beam.
  In this paper, the results of the vector analyzing powers, A_x and A_y, and the cross section for data taken with a proton-beam energy of 135 MeV are presented. Moreover, we performed a global review of a rich set of cross section and vector analyzing-power data taken with proton-beam energies of 135 MeV and 190 MeV. The results were compared with theoretical Faddeev calculations using 2N and 2N+3NF models such as CD-Bonn, CDB+Δ, CDB+Δ+Coulomb and CDB+TM' <cit.> for the kinematics in which both protons scatter to polar angles  smaller than 30^∘ and with a relative azimuthal opening angle varying between 20^∘ and 180^∘. The results of cross sections and analyzing powers, A_x ans A_y, as a function of S for different configurations (θ_1, θ_2, ϕ_12) are shown in Figs. <ref>–<ref>. At small azimuthal opening angles, the calculation, which is based on the extended CDB+Δ and with Coulomb corrections, CDB+Δ+Coulomb, shows a smaller discrepancy with the data than the other calculations. The results show that there is a general disagreement between the data and the calculations including a 3NF. In particular, predictions for the vector analyzing powers show a systematic deficiency at small relative azimuthal angles, which corresponds to small relative energies. In this range, the data for A_y is closest to the three-body calculation that is based on a 2N potential. The addition of 3NF makes the agreement even worse. 
  
  The results of the global review show very large values of χ^2/d.o.f. for the cross sections at specific scattering and relative azimuthal angles. The deviations, independent of the model and beam energy, appear to increase towards large values of ϕ_12. This implies that the present models are not able to provide a reasonable description of the data. The absolute χ^2/d.o.f. for the analyzing powers A_x and A_y are much closer to unity than the ones observed for the cross sections. However, also for these observables no satisfactory agreement is found for all the angular combinations (θ_1, θ_1) and ϕ_12.
  
  
  
  
  § ACKNOWLEDGEMENT
  
  The authors acknowledge the work by AGOR cyclotron
  and ion-source groups at KVI for delivering the highquality
  polarized beam. This work was partly supported by the Polish National Science Centre under Grant No. 2012/05/B/ST2/02556 and 2016/22/M/ST2/00173. The numerical calculations of Bochum-Krakow group were partially performed on the supercomputer cluster of the JSC, Jülich, Germany.
  
  
  
  epj
  
  
  
  
  
  