-
PDF
- Split View
-
Views
-
Cite
Cite
Shoya Ogawa, Ryo Horinouchi, Masakazu Toyokawa, Takuma Matsumoto, Microscopic optical potentials including breakup effects for elastic scattering, Progress of Theoretical and Experimental Physics, Volume 2019, Issue 12, December 2019, 123D04, https://doi.org/10.1093/ptep/ptz128
- Share Icon Share
Abstract
We construct a microscopic optical potential including breakup effects for the elastic scattering of weakly binding projectiles within the Glauber model, in which a nucleon–nucleus potential is derived using the |$g$|-matrix folding model. The derived microscopic optical potential is referred to as the eikonal potential. For |$d$| scattering, calculation with the eikonal potential, i.e., the EP model, reasonably reproduces the result of the exact calculation with the continuum-discretized coupled-channels method at intermediate energies. From the properties of the EP model, the inaccuracy of the eikonal approximation used in the Glauber model is partially excluded. Also analyzed with the eikonal potential is the scattering of |$^6$|He from |$^{12}$|C, and its applicability to scattering with many-body projectiles is shown.
1. Introduction
A microscopic understanding of the nucleon–nucleus (|$NA$|) and nucleus–nucleus (|$AA$|) optical potentials is one of the most important issues in the theory of nuclear reactions. Optical potentials are necessary for describing not only elastic scattering but also reactions including higher-order processes. For example, the distorted-wave Born approximation and the continuum-discretized coupled-channels method (CDCC) [1–3] require the optical potentials between constructs of a projectile (P) and a target (T) to describe inelastic scattering, breakup, and transfer reactions.
The |$g$|-matrix folding model has been widely used as a reliable method for describing optical potentials [4–10]. The |$g$| matrix is an effective nucleon–nucleon interaction [11–28] in nuclear matter and depends on the density |$\rho_{\rm m}$| of the nuclear matter. In the |$g$|-matrix folding model, optical potentials are derived by folding the |$g$| matrix with the target density |$\rho_{\rm T}$| for |$NA$| scattering and with |$\rho_{\rm T}$| and the projectile one |$\rho_{\rm P}$| for |$AA$| scattering. The procedures are referred to as the single-folding (SF) model for |$NA$| scattering and the double-folding (DF) model for |$AA$| scattering.
In the SF model, |$\rho_{\rm T}$| is regarded as being |$\rho_{\rm m}$| in the |$g$| matrix via the local density approximation. The |$NA$| optical potential thus obtained becomes nonlocal by taking into account knock-on exchange processes. The nonlocal potential is impractical in many applications but can be localized with good accuracy by using the Brieva–Rook approximation [14–16]. The localized potential is quite successful in reproducing experimental data systematically when reliable |$g$| matrices such as the Melbourne [23], complex effective Gaussian [19,20,25–27], and chiral effective field theory [28] ones are used.
On the other hand, it is more difficult to calculate the optical potential with the DF model for |$AA$| scattering than it is with the SF model for |$NA$| scattering. The main problem with the DF model is how to treat |$\rho_{\rm m}$| in the |$g$| matrix for |$AA$| scattering. In general, the sum of |$\rho_{\rm P}$| and |$\rho_{\rm T}$| is taken as |$\rho_{\rm m}$| in the |$g$| matrix, a procedure known as the frozen-density approximation (FDA). To reproduce experimental data, the optical potential derived using the DF model with the FDA often requires a normalization factor for its real and imaginary parts. Thus, the choice of |$\rho_{\rm m}$| in the |$g$| matrix for the DF model is a longstanding open problem.
In previous work on the scattering of |$^{3,4}$|He from various targets [29,30], we proposed the target-density approximation (TDA) in the DF model, where |$\rho_{\rm m}$| is estimated using only |$\rho_{\rm T}$|. The DF model with the TDA (the DF-TDA model) does well at reproducing experimental data for |$^{3,4}$|He scattering at intermediate energies with no adjustable parameter, and it has been validated theoretically by using multiple-scattering theory [31–33]. Furthermore, as a practical version of the DF-TDA model, we also proposed the double–single folding (DSF) model for |$^{3,4}$|He scattering, in which the optical potential between |$^{3,4}$|He and the target is derived by folding the localized |$NA$| optical potential of the SF model with the |$^{3,4}$|He densities.
However, the optical potential calculated by the DSF model does not account for breakup effects in elastic scattering. In fact, as shown in Fig. 3(a) in Ref. [30], results with the DSF model differ slightly from experimental data for |$^3$|He scattering at low energies. This discrepancy is due to the breakup effects of |$^3$|He and can be corrected by applying a CDCC calculation, which can describe breakup effects accurately. The CDCC has been used successfully to analyze the scattering of weakly binding two- and three-body projectiles. However, although the CDCC is reliable for analyzing breakup reactions, it is difficult to apply it to the scattering of projectiles with more than five bodies (e.g., |$^{8}$|He) because of the high computational cost.
The Glauber model [34] has been proposed as a means of treating P and T excitation effects for many-body systems. In the Glauber model, the excitation effects are described under the adiabatic and eikonal approximations, which work well for scattering at high incident energies. An advantage of the Glauber model is that it allows an optical potential to be constructed that includes excitation effects. The obtained optical potential is referred to as the eikonal potential (EP) [34–36]. Because the eikonal potential is based on the adiabatic and eikonal approximations, it is expected to be highly valid for scattering at high incident energies. However, the validity of the eikonal potential for scattering at low incident energies is yet to be investigated, and the corresponding validity of the adiabatic and eikonal approximations has hardly been discussed to date.
In this paper, we investigate the validity of the eikonal potential by comparing with CDCC results for |$d$| scattering at 20–200 MeV/nucleon, which is described by a |$p$| + |$n$| + T three-body model. How the adiabatic and eikonal approximations affect the elastic scattering is also discussed. The optical potentials for |$p$|–T and |$n$|–T systems are derived using the SF model with the Melbourne |$g$| matrix. In the analysis, we also discuss the validity of the adiabatic and eikonal approximations for low-energy scattering. Furthermore, to show the applicability of the eikonal potential to many-body systems, we use it to analyze |$^6$|He scattering that is described by a |$^4$|He + |$n$| + |$n$| + T four-body model. Herein, we compare the EP model and the four-body CDCC for the first time.
The paper is organized as follows. In Sect. 2, we describe the theoretical framework for deriving the eikonal potential for |$d$| scattering. In Sect. 3, we present numerical results and discuss the properties of the eikonal potential. Finally, we summarize the paper in Sect. 4.
2. Theoretical framework
2.1. Glauber model
Note here that although |$f_{\rm el}$| of Eq. (14) depends on |$a_{\rm C}$|, the differential cross section defined as |$|f_{\rm el}^{\rm G}(\boldsymbol{q})|^2$| does not. In this paper, we refer to the calculation with |$f_{\rm el}^{\rm G}(\boldsymbol {q})$| as the Glauber model.
2.2. Eikonal potential model
The details of |$U_{\rm E}$| are given in Refs. [34–36], and |$U_{\rm E}$| thus obtained is referred to as the “eikonal potential”.
In the eikonal potential model (EP model), we solve Eq. (16) without the eikonal approximation, namely, with a fully quantum-mechanical approach. Meanwhile, solving Eq. (16) with the eikonal approximation leads to the same result as that of the Glauber model by definition. Therefore, the difference between the Glauber model and the EP model is whether the eikonal approximation is used, and we expect that the inaccuracy of the eikonal approximation in the EP model may be absent from the Glauber model even though |$U_{\rm E}$| is based on the adiabatic and eikonal approximations. The validity of the conjecture for the EP model will be investigated by analyzing |$d$|–T scattering.
2.3. |$NA$| optical potential
In the present study, we use the Melbourne |$g$| matrix [23] constructed from Bonn-B nucleon–nucleon interaction [41]. The Melbourne |$g$| matrix reproduces well the nucleon elastic scattering either with or without localization of the exchange term [23,42]. For |$^{208}$|Pb, the matter density is calculated using the spherical Hartree–Fock method with the Gogny-D1S interaction [43], and the spurious center-of-mass motion is removed using the standard procedure [10]. For |$^{12}$|C, we take the phenomenological proton density determined from electron scattering [44]; here, the finite-size effect of the proton charge is unfolded using the standard procedure [45], and the neutron density is assumed to have the same geometry as the proton one.
3. Results and discussion
3.1. Differential cross section for |$d$| scattering
We consider the scattering of |$d$| from |$^{12}$|C and |$^{208}$|Pb at |$E_{\rm in}/A_{\rm P}=$|20–200 MeV. In the present analysis, the CDCC results are regarded as the exact solution. The CDCC has been used to analyze the scattering of |$d$|, and it reproduced the experimental data well. Figure 1 shows the differential cross section as the ratio to Rutherford cross section for |$d$| scattering calculated using the CDCC and the DSF model, where the latter neglects the breakup effects of |$d$|. The solid and dashed lines are the CDCC and DSF results, respectively. Clearly, the present CDCC calculations reproduce the experimental data well. The minor differences between the data and the calculation around the dips are due to neglecting the spin–orbit interaction between |$d$| and the target.
![Differential cross section as the ratio to Rutherford cross section depending on transfer momentum $q$ for $d$ scattering from $^{12}$C targets at $E_{\rm in}/A_{\rm P} = 20$–$80$ MeV and $^{208}$Pb targets at $E_{\rm in}/A_{\rm P} = 40$ MeV. The solid line is the result of the CDCC and the dashed line is that of the DSF model. The experimental data are from Refs. [46–48]. Each cross section is multiplied by the factor shown in the figure, as are the corresponding experimental data.](https://oup.silverchair-cdn.com/oup/backfile/Content_public/Journal/ptep/2019/12/10.1093_ptep_ptz128/4/m_ptz128f1.jpeg?Expires=1750235732&Signature=i-Dtct1bxhMSqtkIwLCPEvVv2B9veDmJ2a1ZKxI920tu5KTzJvX8jIxYtvov5zraBmCOUTgPmkK-vIl5cf6x6FkhihpgkvhPKW7tdluMTPE~s8YIptfw8MkjAVjvogUpbkFIxDIoslVeSMOjFKtcV9qmqIIy9-LVtOqehr-veD2jDrNm1RDM1J5STKa32WeBjs1cdRp0bd5AwlTAQo9spel6nhQ1ZOg03YBCsllRElmeU7nnho6tLmOwLGkx4sEIBLSDKePQjxnd0Mu6UdoUC~R-pWzUgEbvs4QSPuiS1FrrYzdC41rb4AQmmNHJVEeuZfywH-nrMzY34NCHFBOAGA__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
Differential cross section as the ratio to Rutherford cross section depending on transfer momentum |$q$| for |$d$| scattering from |$^{12}$|C targets at |$E_{\rm in}/A_{\rm P} = 20$|–|$80$| MeV and |$^{208}$|Pb targets at |$E_{\rm in}/A_{\rm P} = 40$| MeV. The solid line is the result of the CDCC and the dashed line is that of the DSF model. The experimental data are from Refs. [46–48]. Each cross section is multiplied by the factor shown in the figure, as are the corresponding experimental data.
In Fig. 2, we compare the results of the CDCC with those of the EP model. The solid, dashed, and dot-dashed lines are the results of the EP model, DSF, and CDCC, respectively. The EP model reproduces the CDCC results well at high incident energies, thereby showing that the EP model can treat the breakup effects of |$d$| reasonably. At low incident energies, there are some deviations between the results of the CDCC and EP model due to the inaccuracies of the adiabatic and eikonal approximations used in the EP model. We will discuss the validity of the adiabatic and eikonal approximations for |$d$| scattering.

Differential cross section as the ratio to Rutherford cross section depending on transfer momentum |$q$| for |$d$| scattering from (left) |$^{12}$|C and (right) |$^{208}$|Pb targets at |$E_{\rm in}/A_{\rm P} = 20$|–|$200$| MeV. The solid line is the result of the EP model, the dashed line is that of the DSF model, and the dot-dashed line is that of the CDCC. Each cross section is multiplied by the factor shown in the figure.
Figure 3 shows the results of the Glauber model, adiabatic-CDCC, and CDCC for |$d$| scattering. For the adiabatic-CDCC, we solve the CDCC equation with the adiabatic approximation, where the excited energy of each state is replaced by the ground-state energy. The solid, dashed, and dot-dashed lines are the results of the CDCC, adiabatic-CDCC, and Glauber model, respectively. The difference between the CDCC and adiabatic-CDCC results is due to the inaccuracy of the adiabatic approximation, and that between the Glauber model and adiabatic-CDCC results is due to the inaccuracy of the eikonal approximation. Clearly, the three calculations are consistent at high incident energies, where both the adiabatic and eikonal approximations work well. For the scattering at low incident energies, the results of the Glauber model are distinct from those of the adiabatic-CDCC and CDCC while the latter two are almost the same. This shows that the eikonal approximation worsens quickly as the incident energy decreases, whereas the adiabatic approximation still holds even at |$E_{\rm in}/A_{\rm P}=$|20 MeV. Furthermore, it is difficult to apply the eikonal approximation to the scattering of a heavy target. In fact, the results of the Glauber model for the scattering of |$d$| from |$^{208}$|Pb are inconsistent with the CDCC results even at forward angles. This is because the eikonal approximation cannot describe exactly the scattering with the Coulomb interaction, as shown in Sect. 2.1.

Differential cross section as the ratio to Rutherford cross section depending on transfer momentum |$q$| for |$d$| scattering from (left) |$^{12}$|C and (right) |$^{208}$|Pb targets at |$E_{\rm in}/A_{\rm P} = 20$|–|$200$| MeV. The solid line is the result of the CDCC, the dashed line is that of the adiabatic-CDCC, and the dot-dashed line is that of the Glauber model. Each cross section is multiplied by the factor shown in the figure.
In Fig. 4, we compare the results of the EP model with those of the Glauber model and adiabatic-CDCC for |$d$| scattering at 20–40 MeV. Clearly, the results of the EP model reproduce those of the adiabatic-CDCC well. In particular, the EP model is more consistent with the adiabatic-CDCC around |$q=2$| fm|$^{-1}$| than with the CDCC. This shows that the deviations between the results of the CDCC and EP model at low energies as shown in Fig. 2 come mainly from the inaccuracy of the adiabatic approximation. On the other hand, the inaccuracy of the eikonal approximation in the EP model is partially excluded. Thus we conclude that the eikonal potential can be regarded as an optical potential that includes breakup effects within the adiabatic approximation only, even though the eikonal potential is based on both the eikonal and adiabatic approximations.

Differential cross section as the ratio to Rutherford cross section depending on transfer momentum |$q$| for |$d$| scattering from (left) |$^{12}$|C and (right) |$^{208}$|Pb targets at |$E_{\rm in}/A_{\rm P} = 20$|–|$40$| MeV. The solid line is the result of the EP model, the dashed line is that of the Glauber model, and the dot-dashed line is that of the adiabatic-CDCC. Each cross section is multiplied by the factor shown in the figure.
3.2. Differential cross section of |$^6$|He scattering
Here, |$\boldsymbol {\xi}$| is the internal coordinate of |$^6$|He and |$U_n$| (|$U_\alpha$|) is the optical potential between |$n$| (|$^4$|He) and T. For |$U_n$|, we use the same potential as that used in the analysis of |$d$| scattering, and |$U_\alpha$| is obtained by folding |$U_n$| and |$U_p$| with the |$^4$|He density. The ground-state wave function |$\Phi_0$| of |$^6$|He is calculated using the Gaussian expansion method (GEM) with a |$^4$|He + |$n$| + |$n$| three-body model. In the GEM, the model Hamiltonian of |$^6$|He and the parameter set of the Gaussian basis functions are the same ones used in Ref. [52].
Figure 5 shows the results of the EP model (solid line), DSF model (dotted line), and CDCC (dot-dashed line) for |$^6$|He scattering from |$^{12}$|C at |$E_{\rm in}/A_{\rm P}=20$|–200 MeV. The dots are the experimental data [53]. Clearly, the breakup effects of |$^6$|He as represented by the differences between the results of the DSF model and CDCC are stronger than those for |$d$| scattering. Furthermore, the EP model reproduces the CDCC results well, even at |$E_{\rm in}/A_{\rm P}\leq40$| MeV. One reason for this is that the condition represented by Eq. (9) is satisfied well because |$K$| of the |$^6$|He scattering is larger than that of the |$d$| scattering for the same |$E_{\rm in}/A_{\rm P}$|. Thus we expect the EP model to work well for scattering with heavier projectiles.
3.3. Dynamical polarization potential
Finally, we discuss the breakup effects of the optical potential for |$^6$|He scattering. Figure 6 shows the optical potentials derived using the EP and DSF models for |$^6$|He scattering from |$^{12}$|C at |$E_{\rm in}/A_{\rm P} = 40$| MeV. The solid line corresponds to the eikonal potential and the dashed line corresponds to the DSF potential, the latter excluding as it does the breakup effects of |$^6$|He. The dot-dashed line corresponds to the DP potential |$U_{\rm DP}$| defined as Eq. (20), which represents the breakup effects of |$^6$|He on the optical potential between |$^6$|He and |$^{12}$|C. In Fig. 6, the strength of the DP potential is multiplied by four. In the EP model, the DP potential is independent of the orbital angular momentum |$L$|, although the DP potential calculated using the CDCC with the equivalent approach [1,49] depends on |$L$| and has singularities at the nodes of the total wave function. As another approach to obtain the |$L$|-independent DP potential from the coupled-channels calculation, the weight mean approach has been proposed [54]. It will be interesting to compare the DP potentials in future work.

|$R$| dependence of potentials for |$^6$|He scattering from |$^{12}$|C at |$E_{\rm in}/A_{\rm P}=40$| MeV. The solid and dashed lines represent the eikonal potential and the DSF potential, respectively. The dot-dashed line represents the DP potential and is multiplied by four.
In Fig. 6, the real part of the DP potential is repulsive in the peripheral region and the imaginary part is absorptive in the whole region. The behavior of the DP potential is almost the same as the previous result of the |$L$|-dependent DP potential in Ref. [49]. To understand the behavior, we calculate the |$S$| matrices for the scattering. In the eikonal approximation as shown in Eqs. (29) and (30), the real part of the optical potential is related to the argument of the |$S$| matrix, and the imaginary part of the optical potential is the absolute value of the |$S$| matrix. Figure 7 shows the |$S$| matrices for the scattering. The crosses, circles, and squares are the results of the EP, DSF models, and the CDCC, respectively. The |$S$| matrices are plotted from the grazing angular momentum |$L=25$|. Here, this momentum is defined as the position of the peak in the partial breakup cross section. The differences between the results of the DSF model and the CDCC represent the breakup effects on the |$S$| matrix, and the results of the EP model agree well with those of the CDCC. Clearly, the absolute value of the |$S$| matrix for the CDCC is smaller than that for the DSF model, and the argument of the |$S$| matrix for the CDCC is also smaller than that for the DSF model. This shows that the breakup effects make the real part of the optical potential shallow and the imaginary one deep.

|$S$| matrices for |$^6$|He scattering from |$^{12}$|C at |$E_{\rm in}/A_{\rm P}=40$| MeV plotted from the grazing angular momentum. The crosses, circles, and squares are the results of the EP, DSF models, and the CDCC, respectively. |$L$| is the orbital angular momentum and increases toward |$S=(1,0)$|.
4. Summary
We constructed a microscopic optical potential including the excitation effects of P and T based on the Glauber model, namely, the eikonal potential. To investigate the validity of the eikonal potential, we analyzed deuteron scattering from |$^{12}$|C and |$^{208}$|Pb at 20–200 MeV and compared the results of the EP model with those of the CDCC. The EP model reproduced the CDCC results well with increasing incident energy, from which we found the eikonal potential to be a reliable potential including the breakup effects of |$d$| at intermediate energies. Meanwhile, for low-energy scattering, there were some deviations between the results of the EP model and those of the CDCC. Those deviations came from the inaccuracy of the adiabatic approximation in the EP model. On the other hand, the eikonal approximation could be partially excluded from the EP model. Thus, the EP model provides results similar to those of the adiabatic approximation, and the eikonal potential is regarded as an optical potential that includes breakup effects within the adiabatic approximation.
We also used the EP model to analyze |$^6$|He+|$^{12}$|C scattering at 20–200 MeV, which was described by a |$^4$|He + |$n$| + |$n$| + |$^{12}$|C four-body model. In the analyses, the EP model also reproduced well the results of the four-body CDCC, which has been successful in analyzing the breakup reactions of |$^6$|He. This shows the applicability of the EP model to many-body systems. In future work, we will use the EP model to analyze |$^8$|He scattering, which is described by a |$^4$|He + |$n$| + |$n$| + |$n$| + |$n$| + T six-body model.
On the other hand, the EP model has some problems, such as the inaccuracy of the adiabatic approximation and the treatment of the Coulomb interaction, which are important for various breakup reactions. Some approaches have been proposed to solve these problems [55,56]. For example, an eikonal-CDCC could reduce the eikonal potential without the inaccuracy of the adiabatic approximation.
Acknowledgements
This work was supported in part by a Grant-in-Aid for Scientific Research (No. JP18K03650) from the Japan Society for the Promotion of Science (JSPS).