Initial State Interaction for the 20 Ne + 130 Te and 18 O + 116 Sn Systems at 15.3 AMeV from Elastic and Inelastic Scattering Measurements

: Double charge exchange (DCE) reactions could provide experimentally driven information about nuclear matrix elements of interest in the context of neutrinoless double- β decay. To achieve this goal, a detailed description of the reaction mechanism is mandatory. This requires the full characterization of the initial and ﬁnal-state interactions, which are poorly known for many of the projectile-target systems involved in future DCE studies. Among these, we intend to study the 20 Ne + 130 Te and 18 O + 116 Sn systems at 15.3 AMeV, which are particularly relevant due to their connection with the 130 Te → 130 Xe and 116 Cd → 116 Sn double- β decays. We measure the elastic and inelastic scattering cross-section angular distributions and compare them with theoretical calculations performed in the optical model, one-step distorted wave Born approximation, and coupled-channel approaches using the São Paulo double-folding optical potential. A good description of the experimental data in the whole explored range of transferred momenta is obtained provided that couplings with the 2 + 1 states of the projectile and target are explicitly included within the coupled-channel approach. These results are relevant also in the analysis of other quasi-elastic reaction channels in these systems, in which the same couplings should be included.


Introduction
The intrinsic nature of neutrinos is one of major hot cases in fundamental physics. So far, it is unclear whether neutrinos are Dirac or Majorana particles. The clear-cut is the relationship between particle and its own anti-particle. In Dirac's description, elementary neutral fermionic particles and anti-particles are distinctive entities, whereas Majorana indicates that they are the same entity. The hypothesis of the symmetry between particle and anti-particle was presented by Ettore Majorana [1] in a period when neutrinos were also hypothetical particles. Later, the lepton-number conservation, observed in many decays and reactions, placed Majorana's neutrino theory aside for many decades. The discovery of neutrino oscillations brought the topic back to the debate of the scientific community and has been driving theoretical and experimental efforts. Promising experiments to establish the Majorana or Dirac nature of neutrinos are presently conceived of by searching for the neutrinoless double beta (0νββ) decay in nuclei such as 48 Ca, 76 Ge, 116 Cd and 130 Te, 136 Xe, with half-lives higher than 10 25 y [2][3][4]. The observation of the 0νββ decay would prove the Majorana nature of neutrinos, and the measurement of the half-life would allow for the determination of the effective neutrino mass, once appropriate nuclear matrix elements (NMEs) are accurately known [5][6][7].
The nuclear physics community is considering the possibility of extracting data-driven information on 0νββ NME from the study of double charge exchange (DCE) reaction crosssections induced by heavy projectiles. The main purpose of the NUMENand NUREprojects at the INFN-LNS [8][9][10][11] is to provide relevant information on the NME associated with hot cases for the 0νββ decays, by exploring the DCE reactions induced by heavy ions [12]. Other studies at RIKEN [13] and RCNP [14,15] have recently explored DCE reactions in connection to 0νββ and also to populate exotic structures. The NUMEN experimental activity includes measurements of the DCE cross-sections associated with the 0 + → 0 + transition and detection of ejectiles around 0 • at distinct bombarding energies. Furthermore, it is mandatory to assess the contributions of possible competing mechanisms such as twosequential single charge exchange (SCE) reactions and successive multi-nucleon transfers that lead to the same final state. Direct measurements of all steps of these competing reactions are not feasible because most scenarios involve radioactive nuclei. The approach is to measure cross-sections of some steps [16,17] and to fill the gaps by reliable direct reaction calculations.
The development of a consistent microscopic theoretical description of the DCE is essential both for the reaction mechanism and the nuclear structure sides. In the theoretical description of nuclear reactions, the role of the initial-state (ISI) and final-state (FSI) interactions is fundamental for the study of all the reaction channels. They represent the main component of the interacting nuclear potential and give the strongest contribution to the total reaction cross-section. On the contrary, the DCE processes are described as a second order perturbation of the interacting potential. The state-of-the-art reaction theory for DCE and SCE reactions was recently developed within the NUMEN collaboration [18][19][20][21][22]. The distorted wave Born approximation (DWBA) amplitude for the SCE and DCE reactions can be expressed in terms of M αβ = χ (−) β , bB|V|χ (+) α , aA , whereV represents the proper interaction, defined in [18,21] for the SCE and DCE, respectively. The χ (±) αβ symbols denote the incoming and outgoing distorted waves, respectively. The ISI and FSI are responsible for the distortion of the wave functions of the interacting nuclei. Thus, for the purposes of the NUMEN project, a detailed microscopic description of the DCE and SCE reaction mechanisms requires a careful determination of the ISI and FSI for the reactions of interest. These involve 18 O and 20 Ne beams at energies from 15 to 60 AMeV and the target nuclei candidate for the 0νββ decay.
Initial-and final-state ion-ion interactions for DWBA calculations of nuclear reactions can be described by optical potentials (OPs), i.e., by local potentials depending only on the relative coordinate between the colliding nuclei. Elastic and inelastic scattering are the main tools for probing the OPs [23][24][25]. This is usually performed from the analysis of the cross-section angular distributions by fitting the parameters of a complex Woods-Saxon potential within the optical model (OM) calculation. However, this procedure may lead to ambiguities, with different sets of parameters that reproduce the same experimental data equally well [26]. Moreover, such an approach is not suitable when the elastic channel is not experimentally measurable, e.g., for the core-core potentials involved in the multi-nucleon transfer reactions [16]. Alternatively, double-folding OPs are often considered, which fold the frozen densities of the colliding nuclei with a realistic nucleon-nucleon interaction, adopting adjustable normalization factors for the real and imaginary parts.
The presence of a nuclear rainbow structure in the angular distribution of the elastic scattering allows for accurate determination of the OPs. However, a rainbow-like structure appears in α-like systems ( 12 C + 16 O at 330 MeV, for example [27]) or in deformed target nuclei ( 16 O + 27 Al at 100 MeV and 280 MeV [28][29][30][31]) at large scattering angles and, consequently, very low cross-sections. The performance of the OPs can be better evaluated comparing the predictions for non-elastic channels, like inelastic scattering and transfer reactions with experimental data, if available. This can be accurately done by comparing the results of calculations based on the one-step distorted wave Born approximation (DWBA) and the coupled-channels (CCs) methods. In the latter approach, the relevant internal states of the projectile/target systems are explicitly taken into account, and the role of such states in the distorted wave functions of the incoming and outgoing partitions is clarified. The couplings with some inelastic scattering or reaction channels are the other main ingredients that constitute the ISI and FSI.
The present work is part of the NUMEN experimental campaign carried out with 20 Ne and 18 O beams performed at the INFN-LNS at 15.3 AMeV incident energy. Among these, the 20 Ne + 76 Ge [32], 20 Ne + 116 Cd [16,33,34], 18 O + 76 Se [35], and 18 O + 40 Ca [17,36] systems were recently studied. Optical model and coupled-channel calculations are performed to describe the experimental cross-section angular distributions extracted for the elastic and some inelastic transitions. The results obtained using different optical potentials are compared to the data, and a good agreement is reached when the São Paulo double-folding OP [37] is adopted. A relevant conclusion of these studies is that an accurate description of the experimental data requires the inclusion of couplings with the low-lying excited states. In particular, the relevant couplings have the 2 + 1 states of the projectile and target and sometimes, as for example in the 18 O + 40 Ca case [36], also with the 3 − 1 . Here, we discuss new data for the elastic and inelastic scattering of the 20 Ne + 130 Te and 18 O + 116 Sn systems in the same theoretical framework of our previous works [32,[34][35][36]. Both target nuclei are involved in present 0νββ research. In particular, the 130 Te nucleus is a candidate for the 0νββ decay [38], and 116 Sn is the daughter nucleus for the ββ decay of the 116 Cd one [39,40]. As far as we know, no experimental measurements exist at ≈ 15 AMeV incident energy, and no theoretical analyses are present in the literature regarding the case of the 20 Ne + 130 Te system, although some works are focused on elastic scattering with the 20 Ne beam [41,42] or the 130 Te target [43,44]. Regard-ing the 18 O + 116 Sn case, such analyses were performed using lighter projectiles [45] or by fitting procedures with free-parameter Woods-Saxon potentials [46], not useful for the purposes of NUMEN [8].
This paper is organized as follows: the experimental setup and the data reduction technique are described in Section 2; the theoretical framework is presented in Section 3; the obtained results are discussed in Section 4; and the conclusions are given in Section 5.

Experimental Setup and Results
The experiments were performed in Catania (Italy), at the INFN-LNS, using the 20 Ne and 18 O beams accelerated at 15.3 AMeV by the K800 Superconducting Cyclotron [47]. In both cases, ions were fully stripped by crossing a thin carbon foil placed along the beam line. The 20 Ne beam impinged on a 250 ± 10 µg/cm 2 130 Te (99.75% enriched) evaporated on a 42 ± 2 µg/cm 2 carbon backing. A 950 ± 50 µg/cm 2 CH 2 target was placed after the primary one in order to reduce the amount of partially stripped ejectiles emerging from the target [48], which represents a background for other reaction channels measured in the same experimental campaign [8]. The 18 O beam impinged on a 380 ± 20 µg/cm 2 116 Sn target (98% enriched), evaporated onto a thin polyvinyl formal resin backing (4.0 µg/cm 2 thick). The total beam charge was collected by a copper Faraday cup with a 0.8 cm entrance diameter and a 3 cm depth, mounted 15 cm downstream of the target. The accuracy in the charge collection was better than 10%, guaranteed by an electron suppressor polarized at −200 V and a low noise charge integrator circuit.
The 20 Ne and 18 O ejectile momenta were analyzed by the MAGNEX spectrometer [49,50] in separate runs. For the 20 Ne + 130 Te system, the spectrometer optical axis was placed at θ lab = 8 • , 13 • , and 20 • in the laboratory frame, and in the 18 O + 116 Sn system, MAGNEX was placed at θ lab = 10 • and 15 • . The measurements were performed with the spectrometer working in full acceptance mode (Ω ≈ 50 msr). The measured angular range was 3 • < θ lab < 26 • and 5 • < θ lab < 21 • in the laboratory reference frame for 20 Ne + 130 Te and 18 O + 116 Sn, respectively.
For the 20 Ne + 130 Te system, a supplemental run at θ lab = 8 • with a reduced solid angle acceptance (Ω ≈ 34 msr), which intentionally excluded the most forward angles, was performed in order to get a reasonable measurement from the Faraday cup. Indeed, due to the very large cross-section at forward angles and the maximum tolerable rate by the MAGNEX focal plane detector [51,52], the beam current for the run at θ lab = 8 • in the full acceptance condition was limited to about 100 epA, corresponding to an unacceptable signal-to-noise ratio in the measurement of the total beam charge. The details of the particle identification and the data reduction techniques are explained in [33,[53][54][55]. A fully differential algebraic method [56] was implemented to perform the trajectory reconstruction technique, in which the measured horizontal and vertical positions and angles at the focal plane were required as the input. Examples of the obtained energy spectra for the 130 Te( 20 Ne, 20 Ne) 130 Te and 116 Sn( 18 O, 18 O) 116 Sn scattering are shown in Figure 1 in which E x = Q 0 -Q, where Q 0 is the ground-state-to-ground-state Q-value (Q 0 = 0 for elastic scattering). The multiple-fit procedure plotted in Figure 1 was performed to describe the energy spectra. The width of each peak was fixed according to the experimental energy resolution, including the recoil energy broadening due to the in-flight decay of the ejectile for the transitions in which it was found in a bound excited state.
In the case of 130 Te( 20 Ne, 20 Ne) 130 Te scattering (Figure 1a), the obtained energy resolution of ≈ 1 MeV full width at half-maximum (FWHM) is not enough to isolate the transition to the first excited state of 130 Te at 0.839 MeV. A second peak is visible in the spectrum at ≈  The absolute cross-sections were extracted following the procedures discussed in [55], taking into consideration the overall MAGNEX efficiency [57]. The experimental crosssection angular distributions extracted for the quasi-elastic (elastic + 2 + 1 at 0.839 MeV) in the 20 Ne + 130 Te system, elastic in the 18 O + 116 Sn one, and some inelastic transitions are shown in Figures 2-5. The corresponding scale of transferred linear momentum q is also given. The error bars include uncertainties coming from the statistical contribution, solid angle estimation, and fitting procedure.
The elastic scattering data shown in Figure 2b for the 18 O + 116 Sn system are in good agreement with the Rutherford cross-section at very forward angles, without the need for any scaling factor, whereas in the 130 Te( 20 Ne, 20 Ne) 130 Te scattering (Figure 2a), a scale factor equal to 1.5 was applied to all the experimental points to ensure the same agreement. This value is larger than the systematic errors coming from the target thickness and the total charge collected by the Faraday cup (≈10%). This was associated with efficiency losses in the focal plane detector induced by the high counting rate, not recoverable using the standard procedure [57].
The quasi-elastic and elastic scattering cross-sections expressed in terms of the ratio to the Rutherford one are shown in Figures 3 and 5a, respectively. The scattering is dominated by the Coulomb field up to the grazing angle (θ gr ), located at θ gr ≈ 16 • and θ gr ≈ 14 • in the center-of-mass (c.m.) reference frame for the 20 Ne + 130 Te and 18 O + 116 Sn systems, respectively. Beyond θ gr , the distributions show the characteristic fall-off associated with the near-side scattering amplitudes, since in that region, the data are more responsive to the nuclear part of the nucleus-nucleus potential. Differential cross-section angular distributions were extracted also for some inelastic transitions. In particular, Figure 4 shows the angular distributions for the 20 Ne + 130 Te system corresponding to the peak at ≈1.6 MeV (Panel a) and at ≈2.5 MeV (Panel b). The angular distribution for the peak at ≈1.9 MeV in the 18

Theoretical Analysis
The data were compared with theoretical calculations performed in the OM, one-step DWBA, and CCs approaches using the FRESCO code [58]. The double-folding São Paulo potential (SPP) [37] was adopted as the bare potential (V SPP (R)) of the OP, with real and imaginary strength factors N r and N i , respectively. The double-folding V SPP (R) is calculated using two-parameter Fermi distributions for the matter densities with radius and diffuseness parameters obtained from a systematic analysis of electron scattering data and mean field calculations based on the Hartree-Fock-Bogoliubov theory [37]. In particular, the average matter density diffuseness (a m ) is 0.56 fm, although small variations are allowed to accommodate nuclear structure effects that mostly present at the surface. In particular, in the 18 O projectile, to account for the effect generated by the two valence neutrons bounded to the 16 O core, we adopted a larger a m = 0.61 fm, as typically done [35,36,[59][60][61].
In the OM framework, the one-channel scattering equation is solved, and the OP accounts for the average effect from other reaction channels. In a step further, the onestep DWBA solves sets of two equations, consisting of an elastic and an another channel. These channels are connected by a coupling potential, although back coupling effects are not considered. Consequently, the elastic cross-section within this approach does not deviate from the one obtained from the simple one-channel calculation. This is a reasonable approach whenever the coupling potential is weak. Finally, the CC considers the full coupling scheme, consisting of elastic and inelastic channels. This also includes back coupling effects, not present in the one-step DWBA. Therefore, comparison between the results from one-step DWBA and CCs tells about the importance of coupling potentials on the elastic channel.
OM calculations, using the SPP defined in Equation (1) with N r = 1.0 and N i = 0.78, describe the elastic scattering data of many systems [62]. This set of parameters was used also in the OM and one-step DWBA calculations presented here. In the CCs' calculations we adopted a smaller factor, N i = 0.6, since, in such cases, we are explicitly including the effect of some non-elastic channels in the elastic cross-sections. The imaginary part of the OP effectively takes into account the loss of flux to dissipative processes such as deep inelastic collisions. Within the CC and coupled reaction channel calculations, this value has been shown to describe reasonably well the data for elastic and inelastic scattering in 16  The coupling potentials are key ingredients in such calculations and, in macroscopic models, are usually written in terms of multipole decomposition [24]. For Coulomb excitation, each λ-component depends on the corresponding reduced matrix element for the electric operator, which is connected with the reduced transition probability B(Eλ ↑). Nuclear excitations are usually considered as a deformation of the charge or matter distribution of the nucleus. Deformation of the nuclear surface is decomposed into spherical harmonic terms with respect to a spherical shape. The nuclear coupling potentials are determined as: where δ m λ is the matter deformation length of the λ-pole component. Following Khoa and Satchler, the precedent method is referred to as the deformed optical potential [63]. Usually, it is assumed that the optical potential is independent of the channel state and deformed in the same way as the density distribution of the excited nuclei (equal deformation hypothesis). δ m λ is often assumed to be equal for charge and mass distributions. Therefore, the matter deformation length is associated with the charge distribution and determined from the electric reduced transition probability from state I π to I π , B(Eλ; I → I ): where R = 1.2A 1/3 is the nuclear radius. This expression assumes that Z protons are uniformly distributed in a sphere with radius R (the so-called sharp cut-off model), which can be replaced by a more realistic distribution with a finite diffuseness [64]. The equal deformation hypothesis for matter distribution and potential is, in fact, correct if the projectile is a point-like particle and the potential is obtained by folding the density with a zero-range interaction (see Chap. 14 in [24]). For heavy nuclear systems, we adopted a deformation length corrected for the deformed potential to take into account the relative differences in the density and potential radii, as follows: where R pot. is the radius of the optical potential, listed in Table 1 for the two systems together with the volume integrals. The results agree with the typical values [24], confirming a reasonable description of the optical potential properties. In addition, the deformation of the imaginary part is also considered. This is justified as an effective procedure to account for possible couplings of the considered channels to collective states not explicitly included in the coupling scheme [31]. We assume that the imaginary deformation length is equal to the matter deformation (δ Im. λ ≡ δ corr. λ ).  The couplings considered in this work are sketched in Figure 6, Panel (a) for 20 Ne + 130 Te and (b) for 18 O + 116 Sn. The dotted blue and continuous red arrows indicate couplings within the one-step DWBA and CCs' calculations, respectively. The reduced transition probabilities considered for these couplings are listed in Table 2. Table 2. Reduced transition probabilities B(Eλ; I → I ) adopted in this work. The δ corr. 2 is calculated using Equation (4) with the radius of the matter distribution (R) and the radius of the double-folding optical potential (see Table 1).

Nuclei
Transition Initial State Final State B(Eλ; I → I ) δ corr.

The 20 Ne+ 130 Te System
The results of the OM, one-step DWBA, and CCs' calculations for the quasi-elastic and inelastic scattering in 20 Ne + 130 Te are shown in Figures 3 and 4, respectively. In the quasielastic scattering, we expect contributions to the cross-section from the 2 + 1 excited state of 130 Te at 0.839 MeV, which is not separated from the elastic peak due to the energy resolution (see Figure 1a). The results of the theoretical calculations are compared to the scattering data in Figure 3. The elastic channel is calculated within the OM and CC formalisms. The target excitation at 0.839 MeV (2 + 1 in 130 Te) is considered within the one-step DWBA and CCs approaches. According to the calculations, its contribution starts to be relevant compared to the elastic cross-section at larger angles (θ c.m. ≥ 20 • ). Looking at the sum of the two contributions in Figure 3, it appears that the agreement between the data and calculations is good provided that the couplings are fully taken into account in the CCs approach. The DWBA results, indeed, overestimate the experimental cross-section at larger angles, where the effect of couplings start to be more relevant.
The experimental cross-section angular distribution of the peak at ≈1.6 MeV in Figure 1a is compared to the calculations in Figure 4a. The states included in the CCs' calculations are those indicated by the red arrows in the coupling scheme of Figure 6a, i.e., 20 Ne 1.633 (2 + ) + 130 Te g.s. (0 + ) and 20 Ne g.s. (0 + ) + 130 Te 1.632 (4 + ). In the one-step DWBA calculations, only the 2 + 1 state in 20 Ne is present, since it is not possible to reach 4 + transitions within the adopted approach. According to the calculations, the analyzed inelastic cross-section is mainly coming from the excitation of the 2 + 1 state in 20 Ne. Again, the good agreement between the theories and data is achieved when the full coupling effects are explicitly taken into account in the CCs approach. This result is consistent with the conclusions drawn in the 20 Ne + 76 Ge [32] and 20 Ne + 116 Cd [34], in which the CCs' analysis was necessary to reproduce the data, and the inclusion of the 2 + 1 state in 20 Ne gives the strongest contribution in that energy region.
For the analysis of the third peak visible in the 20 Ne + 130 Te spectrum (see Figure 1a), only the 20 Ne 1.633 (2 + ) + 130 Te 0.839 (2 + ) transition was included in the CCs' calculations. The results, shown in Figure 4b, underestimate the experimental data by a factor of ≈2. The contribution of many other excited states of 130 Te (about 30) is expected in this energy region. These channels were not included in our coupling scheme due to a lack of tabulated values on the corresponding B(Eλ; I → I ). Regardless, the contribution of the simultaneous excitation of the 2 + 1 states of the projectile and target is found to give a significant contribution to the total cross-section.

The 18 O+ 116 Sn System
The results obtained for the 18 O + 116 Sn system are shown in Figure 5. The calculations for the elastic transition describe quite well the experimental data both with the OM and CCs approaches (see Figure 5a). In Figure 5b, the cross-section angular distribution for the peak at ≈ 1.9 MeV is shown. This peak, visible in Figure 1b Figure 6b. The contribution of the 2 + 2 state of 116 Sn is found to be negligible, since it is about two orders of magnitude lower than the data. The experimental cross-section is described quite well including the 18 [36], in which the most important contribution to the cross-section was given by the 2 + 1 states of the projectile and target nuclei. As observed for the elastic data, in the explored range of momentum transfer q, the effects of the couplings to the inelastic channels are not visible (see Figure 5b). Indeed, as seen in the 20

Initial State Interactions
In the previous sections, we compared the theoretical CCs' calculations performed with the SPP double-folding potential, defined in Equation (1), with N r = 1.0 and N i = 0.6, to the experimental data of elastic and inelastic scattering, obtaining a quite good agreement for all the cases. The obtained ISI for the two analyzed systems are shown in Figure 7. The nuclear real and imaginary parts of the SPP are shown together with the Coulomb potential as well. The Coulombic radius of the system, R c = 1.2(A 1/3 target + A 1/3 projectile ) fm, is indicated by the green dotted lines. The ISI determined here is a key ingredient for the description of all the other direct reaction channels involving the same projectile-target system at the same incident energy and angular momentum transfer, e.g., multi-nucleon transfer reactions, provided that it is possible to adopt the coupled-channel Born approximation (CCBA) and coupled-reaction channels (CRCs) approaches.
However, the explicit inclusion of the couplings with the relevant excited states is not feasible when more complex reaction nets need to be calculated, such as SCE or DCE reactions, where at most DWBA calculations are typically performed [18,21,22]. Following the approach of [69], it is possible to deduce an effective polarization potential that incorporates the effect of channel coupling in the elastic optical potential. In the present study, this complex polarization potential affects the nuclear potential in the vicinity of the R c , as shown by the continuous and dot-dashed curves. The polarization potential can be added to the determined optical potentials, used to perform the CCs' calculations, to get the coupled-channel equivalent polarization potential (CCEP) [36]. The ISI corresponding to the CCEP for the two explored systems is reported in Figure 7.
The CCEP was successfully used in one-channel calculations giving a quite reasonable description of the CCs' elastic angular distribution [70]. The same approach was recently adopted for DWBA cross-section calculation of SCE reactions [36] and can be extended also to DCE.

Conclusions
This work presents new data for the quasi-elastic and inelastic scattering in the 20 Ne + 130 Te system and for the elastic and inelastic scattering in the 18 O + 116 Sn one at 15.3 AMeV measured within the NUMEN project. A wide range of transferred momenta is explored in a few angular settings, with a good agreement between independent measurements. Up to eight orders of magnitude in the cross-section are explored for the elastic scattering. The uncertainties in the experimental points are small for the whole angular range thanks to a careful tuning of the experimental setup and the advanced techniques adopted during the data reduction.
The experimental results are compared with theoretical calculations performed within the OM, one-step DWBA, and CCs approaches to assess the ISI for the two systems of interest from the NUMEN research perspective. In both cases, a satisfactory description of the elastic and inelastic data is obtained with CCs' calculations, using the São Paulo double-folding potential as the imaginary nucleus-nucleus OP, with N r = 1.0 and N i = 0.6, and deformation lengths corrected to account for the relative differences in the density and potential radii. In the 20 Ne + 130 Te system, we can conclude that the couplings with the 2 + 1 excited state of the projectile are important, whereas for 18 O + 116 Sn, we cannot draw the same conclusion, due to the small momentum transfer explored (q < 4 fm −1 ) in the narrow range of scattering angles. These couplings are also important to determine the ISI in similar systems at similar incident energy and momentum transfer ( 20 Ne + 76 Ge, 20 Ne + 116 Cd, 18 O + 76 Se, and 18 O + 40 Ca).
The determination of the ISI for the 20 Ne + 130 Te and 18 O + 116 Sn systems achieved in this work is fundamental for the description of all the other quasi-elastic reaction channels induced at the same incident energy and angular momentum transfer. We derived two different optical potentials to be used in the full net of reactions of interest for the NUMEN project. When the CCBA and CRC approaches are feasible, the SPP double-folding potential can be used with the explicit inclusion of the relevant couplings investigated here. For the more involved reaction dynamics, such as SCE and DCE reactions, the use of the DWBA approaches with the CCEP can be exploited to incorporate the couplings in an effective way.