Analyses of Low-Energy π −-12 C Elastic Scattering Data

A new updated simple local optical potential is proposed for analyzing low-energy π−-12C elastic scattering data at 80 MeV and below. This potential is composed of two real terms and an imaginary term. The nature of the real part of the potential is repulsive at smaller radii and attractive at larger ones. In fact, the height of the repulsive term is found to change linearly with the incident pion kinetic energy. On the other hand, the imaginary part of the potential is attractive, shallow and non-monotonic with a dip at about 1.6 fm. Such a nature of the potential makes it feasible to predict π−-12C cross sections at other energies in the energy region considered herein. Coulomb effects are incorporated by following Stricker’s prescription. This study will serve positively in studying both pionic atoms and the role of negative pions in radiotherapy.


Introduction
Pions play a major role in studying the nucleus and, as such, they are crucial in nuclear physics and other disciplines [1] [2].For this reason, many pion facilities have been constructed to provide pions with both charges and different energies [3].In the delta resonance region (T π = 200 ± 100 MeV), the pion has a very small mean free path of less than 1 fm and usually faces a complete absorption at the surface of the nucleus.Such a scattering process of pions from different nuclei is usually described by simple optical potentials [4].These potentials are usually reevaluated, and other theoretical models are also used to allow for achieving better results [5].Even though, these potentials are criticized for several drawbacks mainly, not being able to account for large angle data [6].
Recently Shehadeh et al. [7] have proposed a new simple local optical potential which proves to be successful in explaining the angular distributions of elastically scattered charged pions from different nuclei over the whole angular range in the delta resonance region.This potential relies on extracting potential points from available phase shifts using inverse scattering theory (IST), as a guide, and the full Klein-Gordon (K-G) equation with all necessary relativistic kinematical effects.Such a preferable success forms an inducement to extend the use of this potential in explaining the low-energy pion-nucleus data.
In the low energy region (T π < 100 MeV), the pion has a large mean free path of few Fermi which enables it to penetrate deeply in the nucleus.As such, it can be used as a probe for studying the nuclear structure and discovering subtler aspects [8] of the pion-nucleus interaction.The accumulation of large amount of reliable low-energy pion-nucleus elastic scattering data compels theoreticians to build theories and formulate theoretical models in order to explain these data.Although the suggested theoretical potentials showed a reasonable success, but indeed it is not complete.In fact, these potentials have been criticized as difficult, impossible to interpret in physical terms, and didn't provide a completely adequate description of the elastic scattering process [9] [10].In view of this we'll address here the use of our new potential, and testing the extent of its success, in analyzing low-energy π − - 12 C elastic scattering data.
In Section 2, the theory is presented.Section 3 is mainly concerned with results and discussion.The last section draws the conclusions.

Theory
The analytical form of our potential, usually used, has the following form: ( ) Which consists of Satchler's potential form [6] supplemented by the second term in Equation (1).In this study, and for the scattering of low-energy negative pions from a light nucleus ( 12 C-nucleus), one may disregard the W 2 term, i.e.W 2 = 0, in Equation (1).As such, Equation (1) becomes ( ) The parameters in this potential form are determined by implementing the IST outlined in Ref. [11] and will not be repeated here.Then a minimal number of parameters are adjusted to give a general nice fit to the measured cross sections, without losing the good match between the obtained potential points, real and imaginary, using IST from available phase shifts and the analytical forms, real and imaginary, of the potential, respectively.For charged pions, this can theoretically be tested by inserting this spherical symmetric potential in the radial part of K-G equation: In Equation (3), ( ) nl R r is the r times the radial part of the wave function, also k 2 and

( )
U r are given by ( ) where E, m, c, and ( ) V r are the effective pion energy, effective pion mass, velocity of electro-magnetic wave in vacuum and the complex pion-nucleus potential, respectively.Following Stricker's prescription [12], the Coulomb potential V c is implicitly considered as a constant in ( ) V r .For π − -12 C, V c is taken +3.4MeV for an effective Coulomb radius R c = 2.54 fm.
To sustain the possibility of using available standard optical codes, Zemlyanaya et al. [13] have recently summarized the basic equations needed in transforming the true pion mass m π and the pion bombarding energy in the laboratory system K  to the effective mass of the incident pion M π and the actual beam energy ( ) where m T is the target mass.The calculated values using Zemlyanaya et al.'s method agree very well with the values obtained from Satchler's treatment.Zemlyanaya et al.'s method can be considered as another version of Satchler's treatment.Here we adopt Zemlyanaya et al.'s method, and computer codes are modified accordingly.As such, one can start by calculating the pion's momentum in the laboratory system ab P  using the equivalent form of the well-known relativistic energy-momentum relation: and obtain ( ) The relativistic momentum of the pion in center-of-mass system P cm can then be calculated, ( ) So the total energy of the pion in the center-of-mass system E is given by Comparing with Satchler's treatment, , one can find γ π , and then the relativistic wave number of the pion k in the center-of-mass system, ( ) where the constant (0.707447967) is not more than the ratio between the true rest mass of the pion m π c 2 = 139.6 MeV and c  = 197.329MeV•fm.One can now easily calculate the pion's kinetic energy in the center-of-mass system (E c.m. ), . .

2
where the constant (20.90104704) is ( ) ( ) is the reduced mass of the two interacting particles in atomic mass units (u).In atomic mass units, the mass of the target nucleus m T = 12u for 12 C and 0.1499 m u π ≅ .

Results and Discussion
The drawbacks of several theoretical models, as Kisslinger first-order pion-nuclear optical potential, to give a full understanding of π ± -nucleus elastic scattering data obliges theoreticians to solve this problem.As an example, Gmitro et al. [14] have supplemented the first-order pion-nuclear potential by a phenomenological ρ 2 -dependent term to give a reasonable explanation for the data.This was actually hinted by Johnson et al. [15] and embedded in the work of others [16].Although Gmitro et al.'s full potential showed an improvement in analyzing the data, but it is still incomplete.In this work, we'll continue using IST, as a guide, to obtain a potential which gives a remarkable consistent description of pion-nucleus interaction from 0 to 80 MeV.The IST, used in extracting potential points from available phase shifts, has proved to be successful in determining the analytical form of the potential and its parameters.Such a potential has shown a remarkable success in describing nicely the elastic scattering of charged pions from calcium, calcium isotopes, 54 Fe [17] [18] and 12 C above 100 MeV [19], and from 40 Ca below 100 MeV [20].Here, the same strategy has been used in determining the potential used for analyzing π − - 12 C elastic scattering data at several energies; namely at 80 MeV and below.Fortunately, and benefiting from the general trend and behavior of the extracted potential points from available phase shifts, the imaginary and real parameters of the potential are kept unchanged in the two energy sub-regions, 50 ≤ T π < 100 and T π < 50, except for the strength of the repulsive real part V 1 which is found to change linearly with the incident pion kinetic energy T π .The potential parameters in Equation ( 2) are: V 0 = 37.0 MeV, R 0 = 3.75 fm, a 0 = 0.324 fm, R 1 = 3.00 fm, a 1 = 0.333 fm, W 3 = 100.0MeV, R 3 = 1.70 fm, a 3 = 0.370 fm for T π = 80, 69.5, 50 MeV.For T π = 40, 30, 20 MeV, all these parameters are kept unchanged except R 0 = 3.95 fm and W 3 = 70.0MeV.The parameter V 1 is found to change with T π as V 1 = 18, 35, 70, 90, 110, 140 MeV for T π = 80, 69.5, 50, 40, 30, 20 MeV, respectively.These values were fitted to a linear equation form, as shown in Figure 13, and the following relation is obtained: The theoretical predictions of the differential cross sections deduced using our new complex potential defined in Equation ( 2) are compared with the experimental ones in Figure 2, Figure 4, Figure 6, Figure 8, Figure 10 and Figure 12 at 80, 69.5, 50, 40, 30 and 20, respectively.The potentials used in obtaining the theoretical differential cross sections are plotted in Figure 1, Figure 3, Figure 5, Figure 7, Figure 9 and Figure 11.The analytical forms of these potentials, with their real and imaginary parts, are compared with the potential points obtained from available shifts using IST at 80, 50, 40, and 30 MeV.Nevertheless, the calculated reaction cross sections are 332, 321, 181, 157, 152, and 138 mb at 80, 69.5, 50, 40, 30 and 20 MeV, respectively.They are in good agreements with the experimental ones [21] [22].
The nature of the imaginary part of the potential is the same at all energies.It is non-monotonic and shallow with a minimum at about 1.6 fm.On the other hand, the real part is non-monotonic, shallow with a small repulsive core at 80 MeV which grows linearly as the pion's incident kinetic energy decreases.In fact, it is repulsive at small radii (r ≤ 2.5 fm) and attractive at large ones (r > 2.5 fm).On the other hand, Figure 2, Figure 4, Figure 6, Figure 8 and Figure 10 show clearly that the calculated cross sections describe the data reasonably well.
A similar nature for the potential, real and imaginary parts, was reported by Friedman [28].For the 80 and 50 MeV cases, a reasonable number of phase shifts are available [23], and deduced potential points from available phase shifts using IST are compared with the analytical forms of the potentials in Figure 1 and Figure 5. On the other hand, a small number of phase shifts available at 40 and 30 MeV [10] [23] were legitimately used to extract potential points using IST for a light nucleus ( 12 C-nucleus) case with a radius of 2.54 fm.The extracted potential points are also compared with the analytical potential forms in Figure 7 and Figure 9.Nevertheless, and as Figure 8 and Figure 10 show, the agreement between calculated and measured angular distributions is excellent.In essence, it is very impressive to see the same potentials in each energy sub-region, except for V 1 parameter approximated from Equation (12), are very successful in accounting nicely for the data at all six energies: 80, 69.5, 50, 40, 30, and 20 MeV.Such a nature of the potential makes it feasible to predict the differential cross sections at other energies in a certain energy sub-region.As an example Figure 14 shows the predicted differential cross sections for the 60 MeV case, and the reaction cross section is 301 mb.It is also worthy to mention              that successful analyses of data at several energies reduces the ambiguity and selects a unique potential.The use of this successful potential can be tested at energies below 20 MeV for explaining pionic atom data, and a reasonable link may be established between pionic atom results and the scattering data.In contrast, Friedman et al. [29] have used optical potentials known from fits to pionic atom data to explain the observed cross sections.
As an important application for this investigation, Raju [30] has confirmed that negative pions, with energies considered herein, are of considerable interest and importance in radiotherapy.Low-energy negative pions can be captured by main tissue elements and provide kinetic energy enough to damage hypoxic cells in these tissue elements with no harm to other nearby cells.

Conclusion
The new updated potential proves to be successful in explaining low energy π − -12 C elastic scattering data.In fact, the calculated differential and reaction cross sections are in nice agreements with the measured ones.It is interesting to see that the imaginary part of the potential, in each energy sub-region, is energy independent; and only the strength of the repulsive part V 1 of the real part decreases linearly as the incident pion kinetic energy T π increases.This, and for the first time, establishes a relation between V 1 and T π given by Equation (12).It is worth noting that when T π ≈ 0, V 1 ≈ 172 MeV; and when T π > 87 MeV, V 1 becomes negative.The nature of the potential makes it possible to predict the differential and reaction cross sections at other energy in the energy range of 20 -80 MeV.Moreover, this investigation emphasizes the strength of IST in predicting our new updated potential which proves to be very successful in explaining low-energy π − -nucleus elastic scattering data.In addition, it will serve positively in both pionic atom studies and radiotherapy applications.

Figure 1 .
Figure 1.Real and imaginary parts of the potential are compared with the extracted potential points from available phase shifts[23] using inverse scattering theory.

Figure 2 .
Figure 2.The calculated angular distributions (dashed line) using the potential given by Equation (2) are compared with the experimental data[9] (empty triangles) as a function of center of mass angle (θ c.m. ) at T π = 80 MeV for negative pions.

Figure 3 .
Figure 3. Real and imaginary parts of the potential are presented.Unfortunately, no phase shifts are available.

Figure 4 .
Figure 4. Same as Figure 2 but for 69.5 MeV.The experimental data are from Ref. [24].

Figure 6 .
Figure 6.Same as Figure 2 but for 50 MeV.The experimental data are from Ref. [25].

Figure 7 .
Figure 7. Same as Figure 1 but for 40 MeV.The phase shifts are taken from Ref. [10].

Figure 8 .
Figure 8. Same as Figure 2 but for 40 MeV.The experimental points are taken from Ref. [26].

Figure 12 .
Figure 12.Same as Figure 2 but for 20 MeV.The solid triangles[26] and the empty triangles[27] are the experimental differential cross sections.

Figure 13 .
Figure 13.Height of repulsive core versus incident pion kinetic energy.

Figure 14 .
Figure 14.The predicted differential cross sections for 60 MeV incident negative pions on 12 C-nucleus.