Reaction theory: status and perspectives

Reaction theory: Status and perspectives A.M. Moro1and J. G´ omez-Camacho1,2 1Departamento de FAMN, Universidad de Sevilla, Apartado 1065, 41080 Sevilla, Spain 2CNA (U. Sevilla, J. Andaluc´ıa,CSIC),Av.TomasA.Edison7, 41092 Sevilla, Spain Abstract The current status of the reaction theory of nuclear collisions involving weakly-bound exotic nuclei is presented. The problem is addressed within the Continuum Discretized Coupled Channel (CDCC) framework, recalling its foundations and applications, as well as its connection with the Faddeev formalism. Recent developments and improvements of the method, such as core and target excitations and the extension to three-body projectiles, are presented. The use of the CDCC wave function in the calculation of inclusive breakup reactions is also introduced. 1 Introduction Ongoing and planned radioactive beam facilities will produce in near future accurate reaction data covering unexplored areas of the nuclear chart. The extraction of meaningful information from these experiments will demand new advances in nuclear structure and reaction theories. The description of the collision between composite nuclei is a very complicated quantum mechanical many-body problem which, in general, cannot be solved exactly. Although, ideally, the problem could be solved from a many-body viewpoint, starting with some effective nucleon-nucleon interaction, in practice one needs to resort to additional approximations, aimed at making the problem numerically tractable. These approximate models try 201 ,0600 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0600 epj 5 1177 NN 22  © The Authors, published by EDP Sciences SI . This is an open access article distributed under the terms of the Creative Commons Attribution License 4.0 (http://creativecommons.org/licenses/by/4.0/). F - to emphasize specific degrees of freedom, those which are expected to dominate the reaction dynamics, or those in which one is interested, and recast the problem in terms of an effective Hamiltonian in which these selected degrees of freedom appear explicitly. Within the Coupled Channels (CC) method, the Hamiltonian that describes a scattering problem, in which some internal degrees of freedom (from projectile, target or both) can be excited, is: H=Hint(ξ)+ ˆ TR+U(R,ξ),(1) with Hint as the internal Hamiltonian, depending on the internal coordinates ξ,ˆ TRthe kinetic energy operator, and Uthe projectile-target interaction, depending both on their relative coordinate Rand ξ. In CC methods [1], the scattering wave function is written as an expansion of the form ΨCC(ξ,R)= β Φβ(ξ)χβ(R),(2) where Φβ(ξ) are eigenstates of Hint. The index βrepresents as set of quantum numbers characterizing a given projectile or target state. The functions χβ(R) describe the projectile-target relative motion and are obtained by inserting the expansion (2) into the Schr¨odinger equation, [H−E]ΨCC =0. This gives rise to a set of coupled differential equations which, solved with the proper boundary conditions, produce the scattering observables. In practice, the CC method is applied using a partial wave decomposition of the total wave function as explained elsewhere [1,2]. In the following section, we will adapt the general CC method to relevant cases in exotic nuclei scattering. We will identify the relevant internal degrees of freedom ξ, and specify the particular expression of the internal wave function Φβ(ξ). 2 The Continuum Discretized Coupled Channels (CDCC) method When dealing with the scattering of a two-body weakly-bound nucleus on a well bound target, an important degree of freedom that must be taken into account is the dissociation (i.e., breakup) of the projectile into its constituents. Within the CC scheme, this problem can be accurately treated using the CDCC method. This method was originally introduced by G. Rawitscher [3] and later refined by the Pittsburgh-Kyushu collaboration [2, 4] 201 ,0600 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0600 epj 5 1177 NN 22 2 p n r A R 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 111111 111111 111111 111111 111111 111111 111111 111111 111111 111111 111111 111111 111111 111111 111111 111111 111111 111111 111111 111111 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 11111 max H H= 2.22 Me V Il,1 Il,2 Il,2 swaves pwaves dwaves Breakup threshold l=0 l=1 l=2 n=1 n=2 n=3 H  ground state Figure 1: Left: Relevant coordinates for a deuteron-target scattering problem. Right: Illustration of continuum discretization for the same problem. to describe the effect of the breakup channels on the elastic scattering of deuterons. Denoting the reaction by a+A,witha=c+v(referred hereafter as the core and valence particles, respectively), the method assumes the effective three-body Hamiltonian H=Hproj +ˆ TR+UcA(rcA)+UvA(rvA),(3) with Hproj =Tr+Vvc the projectile internal Hamiltonian, ˆ Trand ˆ TRare kinetic energy operators, Vcv the intercluster interaction and UcA and UvA are the core-target and valence-target optical potentials evaluated at the corresponding incident energy per nucleon. In the CDCC method the three-body wave function of the system is expanded in terms of the eigenstates of the Hamiltonian Hproj including both bound and unbound states. Since the latter form a continuum, a procedure of discretization is applied, consisting in representing this continuum by a finite and discrete set of square-integrable functions. So, in this case, ξ≡r, and the internal wave functions in Eq. (2) take the form Φβ(ξ)≡φp n,,jπ(kn,r),(4) where φp n,,jπ(kn,r) is a discrete continuum wave function representing a region of the continuum with orbital (total) angular momentum (j)and nominal wave number kn. The standard CDCC method relies on a strict three-body model of the reaction (c+v+A), and has proven to be rather successful to describe elastic and exclusive breakup cross sections of deuterons and other weakly bound two-body nuclei (like 6,7Li) [2]. However, it has limitations. The assumption of inert bodies is not always justified, since excitations of the projectile constituents (cand v) and of the target (A) may take place in addition to the projectile dissociation. Furthermore, the two-body picture 201 ,0600 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0600 epj 5 1177 NN 22 3 may be inadequate for some nuclei as, for example, in the case of the Borromean systems, i.e., three-body systems in which the binary sub-systems are unbound, such as 6He and 11Li. Some extensions of the CDCC method to deal with the these situations are outlined in the following. 2.1 Inclusion of target excitations The explicit inclusion of target excitation was done by the Kyushu group in the early days of the CDCC method [4]. Target excitation arises from the non-central part of the valence-target and core-target interactions and, consequently, the effective Hamiltonian, Eq. (3), must be generalized to: H=Hproj(r)+Htar(ξt)+ ˆ TR+UcA(rcA,ξ t)+UvA(rvA,ξ t),(5) where Htar(ξt) is the target Hamiltonian, depending on the target degrees of freedom (denoted as ξt). Note that c−Aand v−Ainteractions depend now, in addition to the corresponding relative coordinate, on ξt.So,the relevant degrees of freedom are ξ≡{r,ξ t}and the internal states become Φβ(ξ)≡φp n,,jπ(kn,r)⊗φt m(ξt)JM ,(6) which are just products of the target eigenstates times the continuum discretized projectile wave function. They can be coupled to a total internal angular momentum J, M. In the calculations performed in Ref. [4], the target excitation was treated within the vibrational model and the c−Aand v−Ainteractions represented by deformed optical potentials. They performed calculations for the d+58Ni reaction at Ed= 22 and 80 MeV, including s-wave breakup of the deuteron and excitation of the target to its first excited state (2+). In this case, deuteron breakup was found to dominate over target excitation. Recently, the problem has been also addressed by P. Chau [5] and applied also to deuteron elastic and inelastic scattering, but the method does not seem to have been implemented in publicly available codes. 2.2 Inclusion of core excitations Similarly to the case of target excitation, excitations of the projectile constituents (cand vin our case) are also possible. For example, in the important case of the scattering of two-body halo nuclei, collective excitations of the core can be important. These core excitations will affect not only the reaction dynamics, but also the structure of the projectile. In the inert-core 201 ,0600 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0600 epj 5 1177 NN 22 4 picture, the projectile states will correspond to pure single-particle or cluster states. On the other hand, if the core is allowed to excite, the projectile states will contain in general admixtures of these core excited components. These two effects (dynamic and structure effects) have been recently investigated within extended versions of the DWBA and CDCC methods [6–9]. To this end, the following effective three-body Hamiltonian has been used: H=Hproj(r,ξ c)+ ˆ TR+UcA(rcA,ξ c)+UvA(rvA).(7) Note that the core degrees of freedom (ξc) appear in the projectile Hamiltonian (structure effect) as well as in the core-target interaction (dynamical effect). In the weak coupling limit, the basis functions can be expressed as Φβ(ξ)≡ α [ϕp α(r)⊗ΦI(ξc)]JM ,(8) where ξ≡{ξc,r}.Thelabelαdenotes the set of quantum numbers {, s, j, I},withIthecoreintrinsicspin, j= +s and  J= j+ I.The functions ΦI(ξc)andϕp α(r) describe, respectively, the core states and the valence–core relative motion. For continuum states, a procedure of continuum discretization is used, similar to the case of the standard CDCC. Once the projectile states (8) have been calculated, the three-body wave function is written as an expansion in terms of these states, according to Eq. (2). Calculations using this extended CDCC method (XCDCC) were first performed by Summers et al. [8,10] for 11Be and 17Con9Be and 11Be+p. Later on, De Diego et al. [9] performed calculations for 11Be on 1H, 64Zn and 208Pb. These calculations suggest that, for the proton target case, dynamic core excitations may increase significantly the breakup cross sections whereas, for the heavier targets, the dynamic core excitation mechanism is small but its effect on the structure of the projectile is still important. As an example, we show in Fig. 2 the XCDCC calculations for the reaction 11Be+pat 63.7 MeV/nucleon (adapted from Ref. [9]). A no-recoil DWBA version of this approach has been also proposed in Refs. [6,7]. Calculations performed using this formalism have confirmed the conclusions of the more elaborated XCDCC calculations and have also served to evidence that the interplay between the single-particle and dynamic core excitation mechanisms produce a distinctive effect on the interference pattern of the resonant breakup angular distributions [11]. 201 ,0600 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0600 epj 5 1177 NN 22 5 10 20 30 40 Tc.m. (deg) 0 10 20 30 dV/d:c.m. (mb/sr) (a) Erel=0.0-2.5 MeV XCDCC: with core excitations XCDCC: no core excitation 10 20 30 40 Tc.m. (deg) 0 10 20 30 Tc.m. (deg) (b) Erel=2.5-5 MeV Figure 2: (Color online) Differential breakup cross sections, with respect to the outgoing 11Be∗c.m. scattering angle, for the breakup of 11Be on protons at 63.7 MeV/nucleon. Left and right panels correspond to the neutron-core relative energy intervals Erel=0–2.5 MeV and Erel=2.5–5 MeV, respectively. Adapted from [9]. 2.3 Extension to three-body projectiles To study the scattering of three-body projectiles, such as Borromean nuclei, the Hamiltonian (3) must be generalized in order to take into account the three-body structure of the projectile. For example, for a Borromean system with a structure of the form c+v1+v2one may use the Hamiltonian H=Hproj(r1,r2)+ ˆ TR+UcA +Uv1A+Uv2A,(9) where Hproj(r1,r2) is the projectile (three-body) Hamiltonian (depending now on two relative coordinates), and Uvi(rviA) the valence-target effective interactions. Clearly, the calculation of the projectile states will be much more involved than in the two-body case. The basis functions can be expressed in several ways. For example, using independent coordinates, Φβ(ξ)= α [ϕα(r1,r2)]JM ,(10) where ξ≡{r1,r2}, which are the internal coordinates of the three-body system. The label αdenotes the set of quantum numbers {1, 2,K,...}, required to characterize the three-body state, which may vary depending on the three-body approach used. Note that Φβ(ξ) is an (approximate) eigenstate of the projectile three-body Hamiltonian, with a given energy and angular momentum. 201 ,0600 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0600 epj 5 1177 NN 22 6 0 30 60 90 120 150 Tc.m. (deg) 0 0.5 1 (dV/d:)/(dVR/d:) 9Li+208Pb 11Li+208Pb No continuum 4b-CDCC 11Li+208Pb @ E=29 MeV 9Li 11Li 0 30 60 90 120 150 Tlab(9Li) 1000 2000 3000 4000 dVBU/d: (mb/sr) TRIUMF data 4b-CDCC 11Li+208Pb -> 9Li +X Figure 3: Left: Differential elastic cross section, relative to Rutherford, for the 11Li+208Pb reaction at an incident energy of 29 MeV. The solid (dashed) line is the four-body CDCC calculation including (excluding) 11Li continuum states. The 9Li elastic scattering data are also included for comparison. Right: Angular distribution of 9Li fragments resulting from the same reaction, compared with the four-body CDCC calculation. The data and calculations are from Refs. [14,15]. The formulation and application of this four-body CDCC method was first done by the groups at Kyushu [12] and Seville [13]. As an application of the method, we show in Fig. 3 the elastic (left) and breakup (right) of 11Li on 208Pb at 29 MeV [14,15]. In both panels, the full four-body CDCC calculations are represented by the solid line. 3 Comparison with the Faddeev formalism The CDCC method was originally devised as a physically sound and numerically appealing ansatz for the three-body wave function, rather than as a formally rigorous solution of the three-body problem. The latter was provided by Faddeev in the 60s [16], who showed that the exact solution of this problem can be obtained by solving a system of coupled-differential equations, named the Faddeev equations after him. The numerical solution of the Faddeev equations is very involved but, recently, it has been possible to solve these equations for a number of situations [17,18], thus providing a valuable benchmark to more approximate models. These comparative studies have shown that the elastic and breakup observables calculated with the CDCC method agree in general very well with the Faddeev solution (see Fig. 4). However, there are also kinematical situations [18] in which the agreement is not good, for example, when the incident energy is too low. This result calls for additional studies and, possibly, for extensions or improvements of the CDCC formalism. 201 ,0600 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0600 epj 5 1177 NN 22 7 -60 -40 -20 0 20 40 60 Tp (deg) 100 101 102 103 d4/d:nd:p (mb/sr2) Matsuoka (1982) Faddeev CDCC 12C(d,pn)12C @ Ed=56 MeV Tn=15o Figure 4: Breakup differential cross section for 12C(p,pn)12CatEd=56MeV,as a function of the proton scattering angle and for a fixed neutron detection angle of 15◦. The Faddeev and CDCC calculations are given by the solid and dashed lines, respectively. The circles are the data from Ref. [19]. 4 The problem of inclusive breakup In previous sections, we have considered only exclusive breakup processes in which the three fragments (c,vand A) survive after the collision and are observed in a definite internal state. In particular, when all fragments end up in their ground state, the process is called elastic breakup (EBU). A qualitatively different scenario occurs when the final state of one or more fragments is not specified in the final state. For example, this is the case of reactions of the form A(a, c)X, in which only one of the projectile constituents (the core in this case) is observed. The angular/energy distribution of the cfragments will contain contribution from all possible final states of the v+Asystem. This includes the EBU channel, in which vand Aremain in their ground state, but also vtransfer, breakup accompanied by excitations of A,andv+Afusion [named incomplete fusion (ICF)]. These non-elastic breakup components (NEB) must be added to the EBU component to give the total inclusive breakup. Whereas the calculation of EBU can be accurately done within the CDCC method and other approaches, the calculation of NEB is more involved due to the large number of accessible states. The problem was addressed in the 1980s by several groups, which developed formal techniques to reduce the sum over final states to a closed form. We cite here the theory developed by Ichimura, Austern and Vincent (IAV) [2,20,21], in which the double differential cross section for NEB with 201 ,0600 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0600 epj 5 1177 NN 22 8 respect to the angle and energy of the cfragments is given by d2σ dEcdΩc    NEB =−2 ¯hvi ρc(Ec)ϕv( kc)|Im[UvA]|ϕv( kc)),(11) where ρc(Ec)=kcμc/((2π)3¯h2), UvA is the optical potential describing v+ Aelastic scattering, and ϕv( kc,rvA) is the wave function describing the evolution of the vparticle after dissociating from the projectile, when the core is scattered with momentum  kcand the target remains in its ground state. This function is obtained as a solution of the inhomogeneous equation (E+ v−Kv−UvA)ϕv( kc,rvA)=(χ(−) c( kc)|Vpost|Ψ3b,(12) where Ev=E−Ec,χ(−) c( kc,rcB) is the distorted-wave describing the scattering of the outgoing cfragment with respect to the B≡v+Asystem, obtained with some optical potential UcB,andVpost ≡Vvc +UcA −UcB is the post-form transition operator. This equation is to be solved with outgoing boundary conditions. Although IAV suggest approximating the three-body wave function appearing in the source term of Eq. (11), Ψ3b, by the CDCC one, a simpler choice is to use the DWBA approximation ψ3b v≈χ(+) a(R)φcv(r), where χ(+) a is a distorted wave describing a+Aelastic scattering, obtained with some optical potential, and φcv is the projectile ground state wave function. This DWBA version of the IAV model has been recently revisited by some groups, and applied to several deuteron [22,23] and 6Li induced reactions [24], showing in most cases a very promising agreement with the existing inclusive data. Further applications and developments are in progress. 5 Conclusions The nuclear reaction theory community has developed in recent years quantum mechanical methods, based on the coupled channels approach, which allow to interpret the exclusive scattering experiments of exotic nuclei. These methods have been compared, in some cases, with exact three-body Faddeev calculations. Preliminary extensions of these methods to the calculation of inclusive breakup cross sections is in progress show very encouraging results. References [1] I.J. Thompson, F.M. Nunes, Nuclear Reactions for Astrophysics, Cambridge, UK: Cambridge University Press, 2009 1(2009) 201 ,0600 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0600 epj 5 1177 NN 22 9