Calculation of protonHe elastic scattering between 7 and 35 MeV
Abstract
 Background

Theoretical calculations of the fourparticle scattering above the fourcluster breakup threshold are technically very difficult due to nontrivial singularities or boundary conditions. Further complications arise when the longrange Coulomb force is present.
 Purpose

We aim at calculating protonHe elastic scattering observables above three and fourcluster breakup threshold.
 Methods

We employ Alt, Grassberger, and Sandhas (AGS) equations for the fournucleon transition operators and solve them in the momentumspace framework using the complexenergy method whose accuracy and practical applicability is improved by a special integration method.
 Results

Using realistic nuclear interaction models we obtain fully converged results for the protonHe elastic scattering. The differential cross section, proton and He analyzing powers, spin correlation and spin transfer coefficients are calculated at proton energies ranging from 7 to 35 MeV. Effective three and fournucleon forces are included via the explicit excitation of a nucleon to a isobar.
 Conclusions

Realistic protonHe scattering calculations above the fournucleon breakup threshold are feasible. There is quite good agreement between the theoretical predictions and experimental data for the protonHe scattering in the considered energy regime. The most remarkable disagreements are the peak of the proton analyzing power at lower energies and the minimum of the differential cross section at higher energies. Inclusion of the isobar reduces the latter discrepancy.
pacs:
21.45.v, 25.10.+s, 21.30.x, 24.70.+sI Introduction
ProtonHe () scattering is one of the most commonly used experiments to study the fournucleon system tilley:92a : It involves two charged particles that are stable and easy to detect with acceptable precision, there are no competing channels until MeV proton lab energy, and beyond that only three and fourcluster breakup takes place up to the pionproduction threshold. Much like neutronH () scattering, is dominated by isospin , but has the deciding experimental advantage of having a proton beam and a nonradioactive He target. On the contrary, from the theoretical point of view, is more difficult to calculate than due to the longrange Coulomb force between protons () that gives rise to complicated boundary conditions in the coordinatespace and noncompact kernel in the momentumspace representation. Nevertheless, these difficulties have been solved below the threecluster breakup threshold using three different theoretical frameworks, namely, the hyperspherical harmonics (HH) expansion method viviani:01a ; kievsky:08a , the FaddeevYakubovsky (FY) equations yakubovsky:67 for the wave function components in coordinate space lazauskas:04a ; lazauskas:09a , and the Alt, Grassberger and Sandhas (AGS) equations grassberger:67 for transition operators that were solved in the momentum space deltuva:07a ; deltuva:07b . A good agreement between these methods has been demonstrated in a benchmark for and elastic scattering observables viviani:11a using realistic nucleonnucleon (NN) potentials.
Recently we extended the AGS calculations to energies above three and fourcluster breakup thresholds deltuva:12c ; deltuva:13a . The complex energy method kamada:03a ; uzu:03a was used to deal with the complicated singularities in the fourparticle scattering equations; its accuracy and practical applicability was greatly improved by a special integration method deltuva:12c . This allowed us to achieve fully converged results for elastic scattering and neutronneutrondeuteron recombination into using realistic NN interactions. We note that the FY calculations of elastic scattering have been recently extended as well to energies above the fournucleon breakup threshold lazauskas:12a , however, using a semirealistic NN potential limited to waves.
In the present work we extend the method of Ref. deltuva:12c to calculate the elastic scattering above breakup threshold and compare with existing data for cross sections and spin observables over a wide range of proton beam energies up to MeV. The Coulomb interaction is treated as in Refs. deltuva:05a ; deltuva:07b using the method of screening of the Coulomb potential followed by the phase renormalization of transition amplitudes taylor:74a ; alt:80a . Thus, standard AGS scattering equations with shortrange potentials are applicable. At energies bellow threecluster threshold our results agree with those obtained by other methods as mentioned in Ref. viviani:11a . Compared to our previous scattering calculations above the breakup threshold deltuva:12c , the most serious complication for is the convergence of the partialwave expansion that requires a larger number of states due to the longer range of the screened Coulomb potential.
Ii 4N scattering equations
We use the symmetrized AGS equations deltuva:07a as appropriate for the fournucleon system in the isospin formalism. They are integral equations for the fourparticle transition operators , i.e.,
(1a)  
(1b) 
Here, corresponds to the partition (12,3)4 whereas corresponds to the partition (12)(34); there are no other distinct twocluster partitions in the system of four identical particles. The free resolvent at the complex energy is given by
(2) 
with being the free Hamiltonian. The pair (12) transition matrix
(3) 
is derived from the potential ; for the pair includes both the nuclear and the screened Coulomb potential . Our calculations are done in momentum space; however, we start with the configurationspace representation
(4) 
where is the true Coulomb potential, is the fine structure constant, is the screening radius, and controls the smoothness of the screening. All transition operators acquire parametric dependence on but it is suppressed in our notation, except for the scattering amplitudes. The symmetrized 3+1 or 2+2 subsystem transition operators are obtained from the respective integral equations
(5) 
The basis states are antisymmetric under exchange of two particles in the subsystem (12) for the partition and in (12) and (34) for the partition. The full antisymmetry of the fournucleon system is ensured by the permutation operators of particles and with and .
The scattering amplitude with nuclear plus screened Coulomb interactions at available energy is obtained from the onshell matrix element in the limit . Here is the Faddeev component of the asymptotic state in the channel , characterized by the bound state energy MeV and the relative momentum , being the average nucleon mass. Due to energy conservation .
The amplitude is decomposed into its longrange part , being the twobody onshell transition matrix derived from the screened Coulomb potential of the form (4) between the proton and the center of mass (c.m.) of , and the remaining Coulombdistorted shortrange part. Renormalizing by the phase factor taylor:74a ; alt:80a ; deltuva:07b , in the limit, yields the full transition amplitude in the presence of Coulomb
(6) 
where the first term is obtained from that converges, in general as a distribution, to the exact Coulomb amplitude between the proton and the c.m. of the nucleus taylor:74a . The renormalization factor is defined in Refs. deltuva:05a ; deltuva:07b . The second term in Eq. (6), after renormalization by , represents the Coulombmodified nuclear shortrange amplitude. It has to be calculated numerically, but, due to its shortrange nature, the limit is reached with sufficient accuracy at finite screening radii . Since the convergence with is faster at higher energies, the required screening radii are smaller than in our lowenergy calculations deltuva:07b . We found that ranging from 8 to 10 fm leads to wellconverged results in the energy regime considered in the present paper. Furthermore, we take a sharper screening with such that at short distances the screened Coulomb approximates the full Coulomb better than with used in Ref. deltuva:07b and at the same time vanishes more rapidly at thereby accelerating the partialwave convergence.
We solve the AGS equations (1) in the momentumspace partialwave framework. The states of the total angular momentum with the projection are defined as for the configuration and for the . Here and are the fourparticle Jacobi momenta in the convention of Ref. deltuva:12a , , , and are the associated orbital angular momenta, and are the total angular momenta of pairs (12) and (34), is the total angular momentum of the (123) subsystem, and are the spins of nucleons 3 and 4, and , , and are channel spins of two, three, and fourparticle system. A similar coupling scheme is used for the isospin. We include a large number of fournucleon partial waves, up to , , and , such that the results are well converged. In fact, lower cutoffs are sufficient for lower , e.g., and are sufficient for . Furthermore, for most observables or even are enough for the convergence; yields small but still visible effect only above MeV.
The numerical calculations are performed for complex energies, i.e., with finite . The limit needed for the calculation of the amplitude is obtained by the extrapolation of finite results as proposed in Ref. kamada:03a . A special integration method developed in Ref. deltuva:12c is used to treat the quasisingularities of the AGS equations (1). We obtain accurate results by using ranging from 1 to 2 MeV at lowest considered energies and from 2 to 4 MeV at highest considered energies. Grid points for the discretization of each momentum variable range from 30 (at lower energies) to 35 (at higher energies). Further details on the other numerical techniques for solving the fournucleon AGS equations can be found in Ref. deltuva:07a .
Iii Results
We study the scattering using realistic highprecision NN potentials, namely, the Argonne (AV18) potential wiringa:95a , the insidenonlocal outsideYukawa (INOY04) potential by Doleschall doleschall:04a ; lazauskas:04a , the chargedependent Bonn potential (CD Bonn) machleidt:01a , and its extension CD Bonn + deltuva:03c allowing for an excitation of a nucleon to a isobar and thereby yielding effective three and fournucleon forces (3NF and 4NF). The binding energy calculated with AV18, CD Bonn, CD Bonn + , and INOY04 potentials is 6.92, 7.26, 7.54, and 7.73 MeV, respectively; the experimental value is 7.72 MeV. Therefore most of our predictions correspond to INOY04 as it is the only potential that nearly reproduces the experimental binding energy of . The calculations with other potentials are done at fewer selected energies.
In Figs. 1  2 we show the differential cross section for elastic scattering as a function of the c.m. scattering angle at a number of proton energies ranging from to 35.0 MeV. This observable decreases rapidly with the increasing energy and also changes the shape; the calculations describe the energy and angular dependence of the experimental data fairly well. Below MeV the experimental data are slightly underpredicted at forward angles as happens also at energies below the threecluster breakup threshold viviani:11a ; deltuva:07b . At the minimum the predictions scale with the binding energy: the weaker the binding the lower the dip of that is located between and . The scaling is more pronounced at higher . For the INOY04 potential that fits the binding energy, one gets an good agreement in the whole angular region up to MeV but, as the energy increases, the calculated cross section underpredicts the data at the minimum much like what happens in nucleondeuteron elastic scattering witala:98a ; nemoto:98c ; deltuva:05a but for nucleon energies above 60 MeV. In line with the conjectures that were made 15 years ago for the threenucleon system witala:98a ; nemoto:98c , this underprediction of the data at the minimum of may be a sign for the need to include the 3NF.
In Fig. 3 we show the proton analyzing power for elastic scattering at proton energies ranging from 7.0 to 35.0 MeV. We observe that the sensitivity to the nuclear force model and energy is considerably weaker as compared to the regime below threecluster threshold deltuva:07a ; deltuva:07b . Most remarkably, in contrast to low energies where the famous  puzzle exists viviani:01a ; fisher:06 ; deltuva:07b , the peak of around 120 degrees is described fairly well but there is a discrepancy in the minimum region. This is similar to the energy evolution of the puzzle in the threenucleon system gloeckle:96a .
In Fig. 4 we show the analyzing power for elastic scattering at proton energies ranging from 7.0 to 35.0 MeV. varies slowly with energy but is slightly more sensitive to the NN potential. Contrary to , calculated is in better agreement with data over the whole energy range, particularly when the INOY04 potential is used.
The experimental data are scarcer for double polarization observables. In Fig. 5 we show the spin correlation coefficient for elastic scattering at 7.03, 8.52, 10.03, and 19.4 MeV, and in Fig. 6 we show for elastic scattering at 8.52 and 19.4 MeV proton energy. Calculated exhibits some sensitivity to the NN potential model and describes the data reasonably well; the agreement with data is the best when the INOY04 interaction is used. The same happens for but for the single data set we know of.
Finally in Fig. 7 we show the proton spin transfer coefficients , , and for elastic scattering at 8.52, 10.77, and 16.23 MeV proton energy. Note that 8.52MeV predictions are compared to experimental data taken at 8.82 MeV but, given the weak energy dependence of these observables, the comparison is appropriate. The calculated spin transfer coefficients show a rich angular structure and follow the data reasonably well but cannot be fully tested by the available data confined to the angular region below . In contrast to other shown spin observables, the spin transfer coefficient around , i.e., in the region of the differential cross section minimum, varies quite rapidly with the energy. There is little sensitivity to NN interaction model, except for around at MeV. Theoretical results and experimental data for are close to 1 up to , but data sets at different energies seem to be inconsistent as they show different angular dependence. In contrast, theoretical predictions at the three considered energies show nearly the same angular dependence for .
As already mentioned, above MeV the minimum of the elastic differential cross section is underpredicted. In order to establish the importance of the 3NF as a means to cure this discrepancy we study the effect of the isobar excitation on both the differential cross section and proton analyzing power at 30 MeV proton energy. This has been done before at energies below threecluster breakup threshold deltuva:08a and we follow here the same procedure. The results are shown in Fig. 8 in a way that one can single out the isobar effects of 2N nature, the socalled 2N dispersion, and of 3N and 4N nature, the 3NF and 4NF. The competition between 2N dispersion and 3NF, often found in the 3N system, is well seen also here for the differential cross section. As shown in Fig. 8 dispersive effect increases the discrepancy with data (dasheddouble dotted curves) while 3NF and 4NF effects reverse that trend for the differential cross section (dotted curves). Nevertheless, when the two effects are put together the net result is an improvement towards the data (solid curves) but not quite enough to bridge the original gap. For only the dispersive effect around the minimum is visible; it moves the predictions away from data.
Iv Summary
In summary, we performed fully converged proton elastic scattering calculations with realistic potentials above the three and fourcluster breakup threshold. The symmetrized Alt, Grassberger, and Sandhas fourparticle equations were solved in the momentumspace framework. We used the complex energy method whose accuracy and the efficiency is greatly improved by the numerical integration technique with special weights. The Coulomb interaction was included rigorously using the method of screening and renormalization.
The differential cross section exhibits rapid energy dependence and, in the minimum region around , also sensitivity to the NN interaction model. The calculations using the INOY04 potential describe the experimental data well up to MeV but underpredict the differential cross section in the minimum at higher energies; other potential models fail even more. In contrast, most of the calculated spin observables show little sensitivity to the interaction model, and also the dependence on the beam energy is weaker than below the threecluster breakup threshold. The overall agreement with the experimental data for the spin observables is quite good, considerably better than in the lowenergy  scattering which is affected by wave resonances. In particular, the peak of the proton analyzing power that is strongly underpredicted at low energies, is reproduced fairly well above MeV but there is discrepancy in the minimum. The observed sensitivity to the NN interaction model seems to be mostly due to different predictions of the binding energy; the calculations using the INOY04 potential with correct binding provide the best description of the experimental data.
We also studied the effect of three and fournucleon forces through the explicit inclusion of the isobar excitation. We found that generated manynucleon forces significantly improve the description of the differential cross section but have almost no effect for the proton analyzing power. However, there are also quite strong dispersive isobar effects that often reduce or even reverse the effect of 3NF. Therefore the total isobar effect, although beneficial, is not large enough to bridge the gap between the differential cross section data and calculations. It even increases the discrepancy in the minimum of . It might be possible that using the standard approach of including static 3NF one might be able to explain the data at higher energies, particularly using Effective Field Theory generated interactions epelbaum:00a ; machleidt:11a . Extension of the method to other reactions in the fournucleon system is in progress.
References
 (1) D. R. Tilley, H. Weller, and G. M. Hale, Nucl. Phys. A541, 1 (1992).
 (2) M. Viviani, A. Kievsky, S. Rosati, E. A. George, and L. D. Knutson, Phys. Rev. Lett. 86, 3739 (2001).
 (3) A. Kievsky, S. Rosati, M. Viviani, L. E. Marcucci, and L. Girlanda, J. Phys. G 35, 063101 (2008).
 (4) O. A. Yakubovsky, Yad. Fiz. 5, 1312 (1967), [Sov. J. Nucl. Phys. 5, 937 (1967)].
 (5) R. Lazauskas and J. Carbonell, Phys. Rev. C 70, 044002 (2004).
 (6) R. Lazauskas, Phys. Rev. C 79, 054007 (2009).
 (7) P. Grassberger and W. Sandhas, Nucl. Phys. B2, 181 (1967), ; E. O. Alt, P. Grassberger, and W. Sandhas, JINR report No. E46688 (1972).
 (8) A. Deltuva and A. C. Fonseca, Phys. Rev. C 75, 014005 (2007).
 (9) A. Deltuva and A. C. Fonseca, Phys. Rev. Lett. 98, 162502 (2007).
 (10) M. Viviani, A. Deltuva, R. Lazauskas, J. Carbonell, A. C. Fonseca, A. Kievsky, L. E. Marcucci, and S. Rosati, Phys. Rev. C 84, 054010 (2011).
 (11) A. Deltuva and A. C. Fonseca, Phys. Rev. C 86, 011001(R) (2012).
 (12) A. Deltuva and A. C. Fonseca, Phys. Rev. C 87, 014002 (2013).
 (13) H. Kamada, Y. Koike, and W. Glöckle, Prog. Theor. Phys. 109, 869L (2003).
 (14) E. Uzu, H. Kamada, and Y. Koike, Phys. Rev. C 68, 061001(R) (2003).
 (15) R. Lazauskas, Phys. Rev. C 86, 044002 (2012).
 (16) A. Deltuva, A. C. Fonseca, and P. U. Sauer, Phys. Rev. C 71, 054005 (2005).
 (17) E. O. Alt and W. Sandhas, Phys. Rev. C 21, 1733 (1980).
 (18) J. R. Taylor, Nuovo Cimento B 23, 313 (1974), ; M. D. Semon and J. R. Taylor, Nuovo Cimento A 26, 48 (1975).
 (19) A. Deltuva, Phys. Rev. A 85, 012708 (2012).
 (20) R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. Rev. C 51, 38 (1995).
 (21) P. Doleschall, Phys. Rev. C 69, 054001 (2004).
 (22) R. Machleidt, Phys. Rev. C 63, 024001 (2001).
 (23) A. Deltuva, R. Machleidt, and P. U. Sauer, Phys. Rev. C 68, 024005 (2003).
 (24) T. B. Clegg, A. C. L. Barnard, J. B. Swint, and J. L. Weil, Nucl. Phys. 50, 621 (1964).
 (25) R. L. Hutson, N. Jarmie, J. L. Detch, and J. H. Jett, Phys. Rev. C 4, 17 (1971).
 (26) B. T. Murdoch, D. K. Hasell, A. M. Sourkes, W. T. H. van Oers, P. J. T. Verheijen, and R. E. Brown, Phys. Rev. C 29, 2001 (1984).
 (27) D. G. McDonald, W. Haeberli, and L. W. Morrow, Phys. Rev. 133, B1178 (1964).
 (28) H. Witała, W. Glöckle, D. Hüber, J. Golak, and H. Kamada, Phys. Rev. Lett. 81, 1183 (1998).
 (29) S. Nemoto, K. Chmielewski, S. Oryu, and P. U. Sauer, Phys. Rev. C 58, 2599 (1998).
 (30) M. T. Alley and L. D. Knutson, Phys. Rev. C 48, 1890 (1993).
 (31) N. Jarmie and J. H. Jett, Phys. Rev. C 10, 57 (1974).
 (32) S. D. Baker, T. Cahill, P. Catillon, J. Durand, and D. Garreta, Nucl. Phys. A160, 428 (1971).
 (33) J. Birchall, W. T. H. van Oers, J. W. Watson, H. E. Conzett, R. M. Larimer, B. Leemann, E. J. Stephenson, P. von Russen, and R. E. Brown, Phys. Rev. C 29, 2009 (1984).
 (34) B. M. Fisher, C. R. Brune, H. J. Karwowski, D. S. Leonard, E. J. Ludwig, T. C. Black, M. Viviani, A. Kievsky, and S. Rosati, Phys. Rev. C 74, 034001 (2006).
 (35) W. Glöckle, H. Witała, D. Hüber, H. Kamada, and J. Golak, Phys. Rep. 274, 107 (1996).
 (36) D. Müller, R. Beckmann, and U. Holm, Nucl. Phys. A311, 1 (1978).
 (37) R. H. McCamis, P. J. T. Verheijen, W. T. H. van Oers, P. Drakopoulos, C. Lapointe, G. R. Maughan, N. T. Okumusoglu, and R. E. Brown, Phys. Rev. C 31, 1651 (1985).
 (38) W. Weitkamp, W. Grüebler, V. König, P. Schmelzbach, R. Risler, and B. Jenny, Nucl. Phys. A311, 29 (1978).
 (39) R. A. Hardekopf and D. D. Armstrong, Phys. Rev. C 13, 900 (1976).
 (40) A. Deltuva, A. C. Fonseca, and P. U. Sauer, Phys. Lett. B 660, 471 (2008).
 (41) E. Epelbaum, W. Glöckle, and U.G. Meissner, Nucl. Phys. A671, 295 (2000).
 (42) R. Machleidt and D. R. Entem, Phys. Rep. 503, 1 (2011).