# Calculation of three-body nuclear reactions with angular-momentum and parity-dependent optical potentials

## Abstract

Angular-momentum or parity-dependent nonlocal optical potentials for nucleon- scattering able to fit differential cross section data over the whole angular regime are developed and applied to the description of deuteron- scattering in the framework of three-body Faddeev-type equations for transition operators. Differential cross sections and deuteron analyzing powers for elastic scattering and transfer reactions are calculated using a number of local and nonlocal optical potentials and compared with experimental data. Angular-momentum or parity-dependence of the optical potential turns out to be quite irrelevant in the considered three-body reactions while nonlocality is essential for a successful description of the differential cross section data, especially in transfer reactions.

###### pacs:

24.10.-i, 21.45.-v, 25.45.Hi, 25.40.Hs## I Introduction

Deuteron scattering from composite nuclei is, in general, a many-body problem, but very often it is treated as three-body problem in a system consisting of proton (), neutron (), and inert nuclear core (1); (2); (3). The interactions between nucleons and core are modeled by effective optical potentials (OPs); there is a large variety of phenomenological parametrizations (4); (5); (6); (7); (8). In some cases such models are successful in describing the experimental data, and in some cases fail, calling for improvements in the OP or the dynamical model, e.g., by including the excitation of the core (9). In particular cases the failures can be seen already in the nucleon-nucleus two-body system. For example, standard optical potentials were found to be unable to account for precise large-angle nucleon-nucleus elastic differential cross section data above 20 MeV lab energy for stable closed-shell nuclei such as or (10). This is not very surprising given significant momentum transfer in that regime and possible excitations. It was argued (11); (12) that these effects give rise to parity- () or angular-momentum () dependent components in the OP. Indeed, with such additional terms the description of large-angle differential cross section data was significantly improved (10); (13). Naturally one may raise the question what consequences the improvements and new OP terms have on three-body reaction observables. Unfortunately, explicitly - or -dependent OPs are not suitable for standard practical calculations within distorted wave Born approximation (DWBA), adiabatic wave approximation (ADWA) and continuum discretized coupled channels (CDCC) frameworks. However, exact Faddeev-type theory, implemented in individual partial-wave representations for all three involved pairs, is capable of using such - and -dependent potentials. Thus, the aim of the present work is to study three-body nuclear reactions with - and -dependent OPs. We choose elastic and transfer reactions as a working example. However, all existing - and -dependent OP parametrizations for nucleon- are local (10); (13), while in Ref. (14) it was found that the Perey-Buck-type exchange nonlocality of the OP is important in three-body reactions, especially for transfer. For this reason we create several nonlocal parametrizations of - and -dependent OPs, fitted to the same experimental data (15); (16) as the local ones (10); (13).

## Ii Nucleon-nucleus potentials

### ii.1 Angular-momentum-dependent optical potential

We start with a nonlocal OP of the form proposed by Giannini and Ricco (17) and augment it with an -dependent part, resulting

(1) |

Here and are potential strengths for various real and imaginary terms, while for each term the shape is given by

(2a) | ||||

(2b) | ||||

(2c) | ||||

(2d) | ||||

(2e) |

with radius , diffuseness , and nonlocality parameter . The dependence on the orbital angular momentum in the form of Woods-Saxon function given in Eq. (2d) is taken over from Ref. (10) where local -dependent potential was constructed. We determine potential strength parameters by fitting theoretical predictions to experimental data for the differential cross section , total inelastic cross section, and proton analyzing power from Refs. (15); (16), i.e., same data that constrain the local OP of Ref. (10). In addition, also geometric parameters , , and have been varied, typically within 10% of original values (17), to improve the fit quality. The parameters defining the -dependence turn out to be comparable to the local case (10), i.e., typically and . Note that some local OP parameters in Ref. (10) depend strongly on the collision energy and that dependence is not smooth, reflecting the fact that backward-angle experimental data exhibit non-monotonic energy dependence, probably due to the presence of resonant proton-nucleus states. We emphasize that in three-body reactions the energy of each interacting two-body subsystem formally runs from the available three-body energy to , but, in order to have a single three-body Hamiltonian and thereby preserve a Hamiltonian theory (18), it is preferable to use two-body potentials with fixed sets of parameters; the results in Sec. III are obtained following this strategy. The ability of the OP to account for the two-body reaction data over a broader energy regime may be important for its success in three-body reactions and deserves investigation. The nonlocal OP, at least to some extent, is able to absorb smooth energy dependence of data into nonlocality, but far less a non-smooth behavior. Thus, we have not achieved a single parameter set describing the data over the full angular regime at all energies. In fact, the analyzing power data are only accounted for center-of-mass (c.m.) scattering angles up to about . Nevertheless, when we fit the experimental data at a given energy, the resulting OP reproduces well the data for other energies at and only fails at backward angles where the differential cross section is very small. This is an important improvement as compared to the local OP of Ref. (10) that yields considerably worse description for the data not included in the fit. An example is shown in Figs. 1 and 2 where the predictions using local and nonlocal -dependent OPs, determined solely by the data at proton lab energy MeV, are compared with data at and 34.1 MeV. Thus, Fig. 1 reflects the quality of the fit at a single given energy, while Fig. 2 reflects the predictive power of fixed-energy OPs for energies not included in the fit. In this latter case nonlocal OPs are more successful, indicating weaker energy dependence of their parameters as compared to the local OP. We present also predictions of parity-dependent OP from the next subsection, , and of nonlocal -independent OP with parameters from Ref. (14), labeled . The latter was not properly fitted to the present data failing at backward angles and, for , also at forward angles , but otherwise provides a rough description of the experimental data. To confirm the conclusion on the superiority of nonlocal -dependent OP in the two-body system, we created a number of parametrizations (19).

As for scattering, the available experimental data (20); (21); (22) are scarcer and less precise. We tried two options for -dependent potential: i) taking over the parameters of the potential, ii) explicitly fitting to experimental data. An example for scattering at MeV neutron energy is presented in Fig. 3. While for the differential cross section the quality is nearly the same in both cases, explicit fitting leads to a better description of the neutron analyzing power. We show results for two sets of parameters to demonstrate large model dependence for backward angle .

### ii.2 Parity-dependent optical potential

The mechanism giving rise to parity-dependent terms in the OP (12); (13); (23). is different from that of -dependent terms, but nevertheless has led to a comparable quality when fitting the data. To get the nonlocal version we replace the term in Eq. (1) by a parity-dependent term, resulting in a nonlocal -dependent OP

(3) |

In addition, we allowed for an imaginary spin-orbit term with strength , following Ref. (13) where -dependent local OP was developed, but fixed nonlocality parameters to original values from Ref. (17). Otherwise the fitting procedure is the same as in previous subsection, and the achieved quality in describing experimental data is comparable to that of -dependent OP as can be seen in Fig. 1. Worth noting are different radial shapes for -dependent and independent terms, in particular significantly smaller diffuseness as compared to , consistently with findings of Ref. (12).

### ii.3 Discussion

- and -dependence of the OP is not really surprizing, since OP describes not a fundamental interaction but an effective one between the nucleon and composite nucleus. Internal degrees of freedom of the nucleus , if taken into account explicitly, would lead to a highly complicated effective two-body nucleon-nucleus interaction. Solving the -nucleon scattering problem with sufficient accuracy is beyond the present capabilities, but attempts have been made to justify the nonstandard OP terms by the effect of simplest internal degrees of freedom of the nuclear core. E.g., Ref. (12) argues that contributions of the core excitation can be approximated by -dependent terms, while Refs. (11); (10) relate the -dependence of the OP to the deuteron channel coupling. For curiosity we verified this concept in a toy model using theoretical results from Ref. (14) for proton elastic scattering on and , but we expect it to be qualitatively valid for any nuclei and . Starting with an -independent OP for , a real binding potential for , and a realistic potential having central, spin-spin, spin-orbit, and tensor terms, results for were obtained by solving exact three-body equations, thereby including breakup and transfer channels to all orders. The resulting three-body elastic cross section could not be fitted well with the two-body standard OP for , but the inclusion of -dependent terms in the OP for , i.e., in case of Ref. (14), significantly improved the fit. Obviously, such an approach is not reliable for a quantitative determination of the OP as it takes into account only one-neutron internal degrees of freedom in the core and relies on the potentials for the and subsystems (that may be -dependent themselves), but it demonstrates that the -dependence of the OP appears. Of course, -independent phenomenological OPs do not exclude the coupling to the deuteron channel, but include it implicitly in an -averaged way by fiting the data. Applying an -dependent potential in a three-body system with present deuteron channel is justified in exact calculations dealing with three pairwise , , and interactions, but may lead to double counting in DWBA-type approaches that generate the wave not through a rigorous solution of the hree-body problem but from a two-body OP. We also expect that if one would attempt to calculate deuteron-nucleus two-body OP starting from a three-body problem with -independent nucleon-nucleus OPs, the resulting OP in a similar way should acquire -dependence.

## Iii Three-body scattering equations

We describe deuteron-nucleus reactions in the framework of exact Faddeev-type three-body scattering equations. We use Alt-Grassberger-Sandhas (AGS) integral equations (24) for transition operators

(4) |

with , free resolvent , three-particle relative motion kinetic energy operator and available energy , and two-body transition matrix

(5) |

The latter is calculated for each pair with the corresponding two-body potential , where in the odd-man-out notation denotes the interaction within the pair (23) and so on. On-shell matrix elements of taken between the corresponding channel states are reaction amplitudes needed for the calculation of scattering observables.

We solve AGS equations in the momentum-space partial-wave representation. We employ three complete sets of base functions . Here , , or , is magnitude of relative momentum within pair , is magnitude of relative momentum between spectator and c.m. of pair , and are orbital angular momenta associated with and , respectively, and are spins of the corresponding particles. All discrete angular momentum quantum numbers are coupled to total angular momentum with the projection , while , and are angular momenta of intermediate subsystems. Using all three sets , 2, and 3 of these basis states allows the calculation of each potential and transition matrix in its proper basis. Obviously, this enables easy inclusion of - and -dependent potentials, in contrast to CDCC and other approximations, where only one set of base functions is being used.

The proton-nucleus Coulomb force is included via the screening and renormalization method (25); (26); (27). For elastic and transfer reactions we obtain well-converged results with Coulomb screening radius around 10 or 12 fm, and including states with up to 3, 8, and 14 for , , and pairs, respectively. For the interaction we take the realistic CD Bonn potential (28) and use potentials from previous section for nucleon-nucleus pairs.

## Iv Results

Using nucleon-nucleus OPs from Sec. II and Refs. (10); (17) and the realistic neutron-proton CD Bonn potential (28) we study elastic scattering and transfer reactions . In the former case there exist differential cross section and deuteron analyzing power data at MeV deuteron lab energy (29). Comparison of those experimental data and our predictions, including four -dependent and three -dependent nonlocal models, is presented in Fig. 4. One can notice immediately that the local -dependent OP from Ref. (10), although being successful in scattering, fails heavily at large angles in scattering, strongly overpredicting the differential cross section. In contrast, nonlocal models, both - or -dependent, slightly underpredict the differential cross section at large angles, but quite reasonably follow its shape. Properly fitted - and -dependent models, both local and nonlocal, provide a reasonable description of deuteron analyzing powers up to or (in some cases, with exception of where has a deep minimum), but deviate from data and from each other at larger scattering angles. The -independent nonlocal OP (17) accounts for cross section data with a quality comparable to nonlocal - and -dependent OPs, but fails for the deuteron vector analyzing power at . This is expected given that is was not fitted to data and shows there a similar discrepancy. Quite surprisingly, the description of all measured deuteron tensor analyzing powers , , and using this model (17) turns out to be quite similar to - and -dependent models of Sec. II and Ref. (10). This may indicate that deuteron tensor analyzing powers are not well constrained by nucleon-nucleus data.

Next we study transfer reactions. In this case the potential must support bound state for the final nucleus . We therefore take real binding potentials from Ref. (31) in partial wave () when calculating transfer to ground state (excited state). Differential cross section results for both reactions at MeV are shown in Fig. 5. By comparing with the experimental data (30) one notices the failure of the local -dependent OP from Ref. (10) over the whole angular regime. In contrast, all nonlocal OPs, - or -dependent or not, provide quite good description of the experimental data. Thus, - and -dependent terms in the OP appear to be quite irrelevant while the nonlocality of the OP turns out to be essential. Similar findings regarding the OP nonlocality in transfer reactions emerged in Refs. (14); (31) where a broader range of reactions was investigated. Note that there is a difference between present calculations and those of in Ref. (14) in the choice of potential: it was real in Ref. (14) but complex here (except for bound state partial wave). In Fig. 6 we present one more example, i.e., transfer to ground state at MeV, not considered in Refs. (14); (31). Again, the account of the experimental data (30) by all employed nonlocal potentials is quite good, while the calculations of Ref. (18) with local but explicitly energy-dependent potentials heavily failed in reproducing this observable.

## V Summary and conclusions

We developed a number of angular-momentum or parity-dependent optical potentials for nucleon- system. Those nonstandard additional terms enabled to fit elastic nucleon-nucleus scattering data at large angles. However, the parameters turn out to be energy-dependent. Optical potential with energy-independent parameters is able to fit two-body data very well around the chosen energy in the whole angular regime for the differential cross section and up to about for the analyzing power. At more distant energies the description remains good at not too large scattering angles . The local -dependent OP from Ref. (10) turns out to be much stronger energy-dependent, with a fixed parameter set able to account for the data in narrow energy region only. In this respect the potentials of the present work represent a significant improvement.

The explicit angular-momentum or parity dependence of the OP can be handled in the Faddeev/AGS three-body scattering equations solved in the momentum-space partial-wave representation where each two-body potential and the corresponding transition matrix is calculated in its proper basis. In an energy-independent form the nonlocal - or -dependent potentials were used to calculate differential cross section and deuteron analyzing powers for elastic scattering and transfer reactions. To isolate the effect of - or -dependence and nonlocality, same observables were calculated using local -dependent (10) and nonlocal -independent (17) potentials. In all considered reactions nonlocal OPs provide quite similar and reasonable description of differential cross section data. In contrast, the predictions using the local -dependent OP (10) strongly deviate from the data and all nonlocal OPs for in elastic scattering and in the whole angular regime for transfer reactions. Based on this fact we conclude that - and -dependent terms in the OP may be quite irrelevant for three-body scattering but the nonlocality plays a major role, especially in transfer reactions; the latter finding is in accordance with Refs. (14); (31); (32); (33). The comparison of predictions and data for deuteron analyzing powers in elastic scattering is less conclusive. The agreement is reasonable for all properly fitted models at not too large scattering angles , but beyond the predictions may deviate from data and from each other. Furthermore, a proper fit to two-body analyzing power data appears to be relevant for deuteron vector analyzing power , but not for tensor analyzing powers , , and .

This work was supported by Lietuvos Mokslo Taryba (Research Council of Lithuania) under Contract No. MIP-094/2015.

### References

- R. C. Johnson and P. J. R. Soper, Phys. Rev. C 1, 976 (1970).
- N. Austern, Y. Iseri, M. Kamimura, M. Kawai, G. Rawitscher, and M. Yahiro, Phys. Rep. 154, 125 (1987).
- A. Deltuva, A. M. Moro, E. Cravo, F. M. Nunes, and A. C. Fonseca, Phys. Rev. C 76, 064602 (2007).
- B. A. Watson, P. P. Singh, and R. E. Segel, Phys. Rev. 182, 978 (1969).
- J. J. H. Menet, E. E. Gross, J. J. Malanify, and A. Zucker, Phys. Rev. C 4 1114, (1971).
- F. D. Becchetti Jr. and G. W. Greenlees, Phys. Rev. 182, 1190 (1969).
- R. L. Varner, W. J. Thompson, T. L. McAbee, E. J. Ludwig, and T. B. Clegg, Phys. Rep. 201, 57 (1991).
- A. J. Koning and J. P. Delaroche, Nucl. Phys. A713, 231 (2003).
- A. Deltuva, Phys. Rev. C 88, 011601(R) (2013).
- A. M. Kobos and R. S. Mackintosh, J. Phys. G 5, 97 (1979).
- R. S. Mackintosh and A. M. Kobos, J. Phys. G 5, 359 (1979).
- G. H. Rawitscher and D. Lukaszek, Phys. Rev. C 69, 044608 (2004).
- S. G. Cooper, Nucl. Phys. A 618, 87 (1997).
- A. Deltuva, Phys. Rev. C 79, 021602(R) (2009).
- J. M. Cameron, J. R. Richardson, W. T. H. van Oers, and J. W. Verba, Phys. Rev. 167, 908 (1968).
- H. B. Eldridge, S. N. Bunker, J. M. Cameron, J. R. Richardson, and W. T. H. van Oers, Phys. Rev. 167, 915 (1968).
- M. M. Giannini and G. Ricco, Ann. Phys. (NY) 102, 458 (1976).
- A. Deltuva and A. C. Fonseca, Phys. Rev. C 79, 014606 (2009).
- Parameter tables and computer codes can be obtained from authors upon request.
- P. Grabmayr, J. Rapaport, and R. Finlay, Nucl. Phys. A 350, 167 (1980).
- S. T. Lam et al., Phys. Rev. C 32, 76 (1985).
- J. P. Delaroche, M. S. Islam, and R. W. Finlay, Phys. Rev. C 33, 1826 (1986).
- S. G. Cooper and R. S. Mackintosh, Nucl. Phys. A 592, 338 (1995).
- E. O. Alt, P. Grassberger, and W. Sandhas, Nucl. Phys. B2, 167 (1967).
- J. R. Taylor, Nuovo Cimento B 23, 313 (1974); M. D. Semon and J. R. Taylor, Nuovo Cimento A 26, 48 (1975).
- E. O. Alt and W. Sandhas, Phys. Rev. C 21, 1733 (1980).
- A. Deltuva, A. C. Fonseca, and P. U. Sauer, Phys. Rev. C 71, 054005 (2005).
- R. Machleidt, Phys. Rev. C 63, 024001 (2001).
- N. Matsuoka et al., Nucl. Phys. A455, 413 (1986).
- M. D. Cooper, W. F. Hornyak, and P. G. Roos, Nucl. Phys. A218, 249 (1974).
- A. Deltuva, Phys. Rev. C 92, 064613 (2015).
- N. K. Timofeyuk and R. C. Johnson, Phys. Rev. C 87, 064610 (2013).
- L. J. Titus, F. M. Nunes, and G. Potel, Phys. Rev. C 93, 014604 (2016).