Roman Shindin Experiment: Quasi-Elastic Reaction at Energy Tn.=0.55-2.0 GeV

1. The nd→p(nn) quasi elastic reactionat energy Tn = 0.55 – 2.0 GeV Roman Shindin Experiment Delta-Sigma Our research program intends to obtain a complete np data set at the zero angle: the measurements of total cross section differencesΔσL(np) and ΔσT(np) for the longitudinal (L) or transverse (T) beam and target polarizations and spin-correlation parameters A00kk(np) and A00nn(np)as well as unpolarized measurements of values Δσ0,tot(np), dσ/dt(np→pn) and Rdpratio of yields of quasi elastic and elastic charge exchange reactions using the D2 and H2 targets Rdp= dσ/dt (nd→p(nn)) / dσ/dt (np→pn) . Main task of these studies is to determine the Real and Imaginary parts of npamplitudes over the energy region 1.2 – 3.7GeV. The knowledge of Rdp could provide an additional constraint and will allow one of some sign uncertainties to be eliminated for the direct reconstruction of the Re parts of the scattering amplitudes. It has appeared that the nd→p(nn) reaction are related with a think effect of Fermi momentum distribution in deuteron. About this subject and also about method of determination such effect will be talking in my presentation

2. Magnetic spectrometer for detection of protons from np elastic or nd quasi elastic charge exchange reactions under 0°Lab • SP-94 - analyzing magnetic dipole • Gx, Gy, 1x, 2x, 3x,y, 4x,y – multiwire propotional chambers • MPT-polarized proton target or liquid D2/H2 target, surrounded by DTS system for detecting of ∆ recoils which decay to π, p and γ • Trigger counters – A, S1, ST1,2,3 and time-of-flight system – S1, TOF1,2

3. Resolution

4. Rejection the np→dπ° andnd→d+π+X background Using the 2-dimensional plot in coordinates of momentum and time-of-flight and with the help of hyperbola cut the background np→dπ° ornd→d+π+X yields rejected from the proton momenta spectra.

5. Rejection the np→dπ° andnd→d+π+X background The capture reactionnp→dπ°andnd→d+π+Xbegin appear after the threshold of meson production 290 MeV. Their yields achieve a maximum aboutenergy 500 - 600 MeV. After that this background decreases.

6. Rejection the np→dπ° andnd→d+π+X background After the threshold Tn > 580 MeV the ∆-resonance excitation appears and become more considerable for large energy. Inelastic events from ∆ production begin to close the elastic peak on the left side, and the calculation of number of elastic ones receive an additional problem.

7. Breit-Wigner + Gaussfit The separation of inelastic and elastic peaks can be performed by the Breit-Wigner and Gauss functions. In this way the mass scale is translated to the momentum.

8. Breit-Wigner + Gaussfit An additional this method was tested by the Monte-Carlo simulation and we obtain a good χ2 ~ 1. But we fined also its error about 15% of elastic events calculation that is related with two cases of ∆-resonance excitation with the target or beam nucleons. For the case of target ∆ the momentum spectrum (see blue histogram) well approximated using the Breit-Wigner function. In the second case the ∆ have any spatial angleθ∆ and the proton from ∆ have itself θp angle too. As a result we obtain more complex distribution (see green) which did not describe by this fit.

9. Suppression inelastic events by the DTS system This equipment is intended for the registration of recoil particle from∆-resonances such as a protons, mesons or gamma quanta. According to our calculation the DTS efficiency for charged and neutral modes of ∆-decay should be equal 92% и67% respectively. It means that the inelastic yields will be suppressed with a factor of 5.

10. Suppression inelastic events by the DTS system The blueshaded histograms present the protons spectra without the DTS signal. The transparent histograms with errors bars show the same as ones but the DTS works in anticoincidence with the spectrometer trigger. The problem of inelastic background is solved for the np→pnreaction at hydrogen target. For the nd→p(nn) charge exchange process obtained at deuterium target the number of inelastic events under left side of elastic peak remains significant.

11. Inelastic events under the right side of elastic peak Table 1. Inelastic background in units %. Calculated without DTS system. Table 2. Inelastic background in units %. Calculated using DTS system. The inelastic yields from the Δ resonances does not go further the centroid of elastic or quasi elastic peaks. At the right side of both peaks this background is negligible. For calculation the numbers of events which are belonged to these peaks we should use only their right parts. It is more simple, but in this way the centroid position must be defined very clearly.

12. Relative shift between thenp→pnandnd→p(nn)peaks At low energy the inelastic backgrounds are almost excluded from protons spectra. It is convenient for observation of elastic and quasi elastic peaks. They have a similar shapes and parametersσH2 and σD2are almost identical. But the quasi elastic peaks is shifted to the smaller valuesrelative the elastic one. This shiftδPequals6± 2 MeV/c.

13. Relation between two Gausses Supposed that the both Gauss functions have equivalent shapes but their centroid are differed among themselves by the shift δP. Then we can define the R function If the value δP is more less then the σ parameter the derivative function dR/dP can be written as very simple formula Thus the relation between G1 and G2 function can be approximated by the direct line and the tangent of line inclination will be proportional to the value of shiftδP.

14. Modeling of relation between two Gausses Let’s call this method as the “Line-Shift” δP tg α

15. Determination ofδP shifting Using the δP value the D2 momentum spectrum is shifted and divided by the H2spectrum. The histogram of D2/H2 yield is fitted by a constant at the region of “well plateau”, and it provides the value of Rdp ratio over the angular range 0 ≤ θ≤ 20 mrad.

16. Determination ofδP shifting The “well plateau” is limited by 300-400 MeV/c after the centroid position of elastic peak at H2 target. The Rdp ratio begin to rise before the left side of “well plateau”, that is related with the inelastic background form the ∆-resonances excitation. Itshows also, that the inelastic events in the right half of elastic or quasi-elastic peaks are negligible.

17. Directly fit by R function The “Line-Shift” method supposed above doesn’t take into account the possible difference among the shapes of elastic and quasi elastic peaks. However we can redefine the R function without this simplification It gives five parameters from which the C and δP should be absolutely free but the parameters MH2, σH2 and σD2 are limited in the frames of errors of their values which are defined preliminary. We shall accept also a sound conditionrelative to the shapes of elastic and quasi elastic peaks

18. Directly fit by R function This method have some advantages. For the first we need not to make the many iterations and can define the δPvalue at once (p3 parameter). For the second the Rdp ratio is defined simultaneously Rdp =C·σD2/σH2≈ C(it is p1·p5/p4).

19. Experimental data of relative shift These points have agreement among themselves. All values are near 7 MeV/c. It show that the relative shift is really effect of quasi elastic reaction. But haw can we get this value from theory?

20. Simple explanation Neutron-spectator Is lost during the nd reaction and recoil neutron takes transfer momentum q Using low of energy conservation it gives For the relative shift we define

21. Simple explanation It gives only qualitative agreement and big differences with the values of δP. The transfer momentumq can be shared to the both neutrons but our calculation will give again this discrepancy. We need more fundamental approach.

22. Correct approach In the frame of impulse approximationthedifferential cross section of nd→p(nn)reaction expressesby the Dean formula which used the spin-flipand spin-non-flipcontributionof np→pncharge exchange process Spin d = 1 Spin (nn) = 0 Spin (nn) = 1 If scattering angle of proton equals zero the t closes to zero too and F(t) ≈ 1. First term on the right side of Dean formula is vanished. Therefore the charge exchange reaction nd→p(nn) under 0º is performed by the spin-flip only. At once after the impact the wave function of relative motion of two neutrons will have the next form There the Ψd(r)is the Hulthen expression for normalized S-wave deuteron wave function having the correct asymptotic behavior. The ρ = 4.31 Fmis the deuteron radius, d =1.7Fm is the effective radius of the low energy neutron-proton interaction and α=6.25.

23. Correct approach

24. Correct approach Two neutrons interact among themselves and (nn)-system as a whole receives the transfer momentumq

25. Correct approach

26. Experimental data of Rdp and rnfl/fl ratios Energy dependences of the ratios Rdp(0) and rnfl/fl(0). The PSA solutions VZ40, FA91 and SP07 were taken from the SAID data base as amplitudes for the np backward reaction and transformed to the representation ofcharge exchange forward process by the unitary transition.