Elastic Differential Cross Sections for Space Radiation Applications
Abstract
The eikonal, partial wave (PW) Lippmann-Schwinger, and three-dimensional Lippmann-Schwinger (LS3D) methods are compared for nuclear reactions that are relevant for space radiation applications. Numerical convergence of the eikonal method is readily achieved when exact formulas of the optical potential are used for light nuclei (), and the momentum-space representation of the optical potential is used for heavier nuclei. The PW solution method is known to be numerically unstable for systems that require a large number of partial waves, and, as a result, the LS3D method is employed. The effect of relativistic kinematics is studied with the PW and LS3D methods and is compared to eikonal results. It is recommended that the LS3D method be used for high energy nucleon-nucleus reactions and nucleus-nucleus reactions at all energies because of its rapid numerical convergence and stability.
pacs:
24.10.Ht, 24.10.Cn, 24.10.JvI Introduction
The space radiation environment is composed of solar particle emissions and ions produced from supernovae distributed throughout the galaxy Ackermann et al. (2013); Benton and Benton (2001). Solar particle events, including both coronal mass ejections and solar flares, are composed of mostly protons with energies that can exceed several hundred MeV. Galactic cosmic rays originate from the shock waves of supernovae and consist of protons and heavier ions with energies that reach hundreds of GeV per nucleon. Radiation transport codes are used to describe the transport of ions, and secondary particles produced from nuclear collisions, from the space radiation environment through shielding materials. Space radiation transport codes require cross sections for the numerous nuclear reactions that occur as a result of collisions of nuclei in the space radiation environment with nuclei in the shield. NASA’s deterministic transport code, HZETRN Slaba et al. (2010a, b); Wilson et al. (1991), currently transports all ions up to nickel—where, thereafter, incident particle fluxes are negligible Simpson (1983)—with energies that extend from MeV to hundreds of GeV per nucleon through shielding materials. Efficient, accurate codes are needed for the computation of nuclear cross sections due to the large number of nuclear reactions that occur at these energies.
The Lippmann-Schwinger (LS) equation is an expression for the scattering transition amplitude Joachain (1983). Scattering amplitudes can be obtained by either solving the LS equation or by employing some approximation, such as the eikonal method. The elastic differential cross section is computed from the absolute square of the scattering amplitude, and the total cross section is related to the imaginary part of the forward scattering amplitude. The elastic cross section is obtained by performing the angular integration of the elastic differential cross section, and the reaction cross section is found from the difference between the total and elastic cross sections.
The input into the elastic scattering equation is the optical potential, which can be expressed in an infinite series of nucleon-nucleon (NN) transition amplitudes, . If the transition matrix is written for ground states of the projectile and target, then, in the factorization approximation, the optical potential is proportional to and the nuclear densities of the projectile and target Picklesimer et al. (1984a, b); Wolfe (1983). The model of used in the present work is parameterized to NN total cross sections, slope parameters, and the real to imaginary ratios of the transition amplitude. Nuclear charge density distributions are obtained from electron scattering experiments De Jager et al. (1974); De Vries et al. (1987). Matter densities of nuclei are found from nuclear charge densities by factoring out the charge distribution of the proton. The internal charge structure of the proton is not taken into account in this analysis; instead, nucleons are treated as point particles. Harmonic well densities are typically used for lighter nuclei because of the Gaussian-like decay of the nuclear charge density as a function of radial distance. Wood-Saxon densities, also known as two-parameter and three-parameter Fermi densities, are better suited for heavier nuclei, where the nuclear charge density is relatively constant before decreasing to zero at larger radial distances.
The two most common ways of solving the LS equation are to use the eikonal approximation or the method of partial wave (PW) decomposition Joachain (1983). The eikonal approximation was first introduced by Moliere and systematically developed by Glauber in the treatment of many-body nuclear reactions with a quantum collision theory of composite objects Joachain (1983); Miller et al. (2007). The eikonal approximation can be derived by assuming high energy and small angle scattering, which leads to a linearized propagator in the LS equation from which the eikonal scattering wave function may be obtained Joachain (1983). The scattering amplitude is determined from eikonal phase factor, which is a function of the the optical potential Wilson et al. (1991); Townsend et al. (1983); Wilson (1974); Wilson and Townsend (1981); Townsend et al. (1982); Townsend (1983).
Besides being an approximation, a drawback of the eikonal approximation is that it may be numerically inefficient for the evaluation of the cross sections for a given optical potential. In the position-space representation, the optical potential, , is given by a 6-dimensional integration for heavy ion collisions Wilson et al. (1991). Therefore, the eikonal phase factor depends on a 6-dimensional integral in the position representation of the optical potential and an additional integration variable over a coordinate in the scattering plane. The numerical integration over 7-dimensions in the position space representation is inefficient when an analytic expression of the optical potential is not known. It is desirable to use exact formulas for the optical potential when analytic expressions of the optical potential can be found. The current work implements expressions of the optical potential for nucleon-nucleus (NA) and nucleus-nucleus (AA) scattering utilizing harmonic well nuclear matter densities for light nuclei (), and the optical potential is expressed in momentum space for cases where no analytic expression can be found () Werneth et al. ().
The LS equation may also be solved via the method of partial wave decomposition Joachain (1983); Werneth et al. (2013), where the transition amplitude is expanded in an infinite series of functions of relative momenta and angular dependent spherical harmonics or Legendre polynomials. After integrating over the angular dependence, the transition amplitude is solved for a given partial wave. Once the partial wave solutions are found, the full solution for the transition amplitude is found by re-summing the series, which is terminated when some pre-defined tolerance of precision is reached.
The PW method is known to become numerically unstable for reactions that require many partial waves Joachain (1983), which is not only limited to high energy NA reactions (GeV/n) but also includes AA reactions at relatively low energy per nucleon (hundreds of MeV/n). The numerical instability can be traced back to highly oscillating Legendre polynomials in the PW expansion and large on-shell momenta for elastic reactions, where contributions to the transition amplitude tend to be localized.
Although there are numerical limitations associated with the PW method, the full three-dimensional Lippmann-Schwinger (LS3D) solution method circumvents the necessity of using highly oscillating Legendre polynomials Elster et al. (1997a); Fachruddin et al. (2000); Rodriguez-Gallardo et al. (2008); Veerasamy et al. (2013); Liu et al. (2005). Most of the LS3D studies have consisted of NN interactions Elster et al. (1997a); Fachruddin et al. (2000); Veerasamy et al. (2013) with the exception of Rodriguez–Gallardo et al. Rodriguez-Gallardo et al. (2008) who studied NA and AA reactions at relatively low energies and Liu et al. Liu et al. (2005) who studied three-body reactions. This demonstrates the validity of the method and can be compared to results generated with the PW method since few partial waves are needed for such reactions. In the present work, the LS3D method is compared to the PW and eikonal methods for NA and AA reactions with energies extending from 150 MeV/n to 20 GeV/n.
The eikonal method is a non-relativistic approximation; however, when energies become sufficiently high, relativistic effects will be manifested in the elastic differential cross section. Relativistic kinematics are needed for high energy reactions and are easily incorporated into the momentum space-representation of the PW and LS3D equations, where the momentum is simply a number instead of a spatial derivative operator, as in the position space-representation. At relativistic energies, the PW and LS3D models will agree if convergence of the partial wave solution is reached, but both methods should differ from the eikonal results, which are non-relativistic. In the low energy limit, the eikonal method should break down and begin to diverge from the PW and LS3D results, since small angle scattering is not appropriate for such reactions. To examine the effect of kinematics, model results are compared for various nuclear reactions at relativistic and non-relativistic energies.
At relativistic energies, the inner structure of the nucleons may be probed. The multiple scattering theory (MST) upon which the model of interaction is based and the NN transition amplitude do not account for the inner structure of the nucleons. The complications associated with the inner structure of the nucleons are assumed to be included in the parameterizations to experimental NN transition amplitudes.
In this paper, exact formulas of the optical potential in the position space representation are used for light nuclei (), and the momentum-space representation of the eikonal phase factor is used for heavier nuclei. The PW and LS3D methods are solved with non-relativistic and relativistic kinematics, and comparisons of the models are made for reactions that are relevant to space radiation. Based on the results presented herein, it is recommended that the LS3D method be used for high energy NA reactions and AA reactions at all energies because of its rapid numerical convergence and stability. The effect of the kinematics for projectiles and targets with equal masses and extensive comparisons to experimental data will be communicated in subsequent manuscripts.
The present work is organized as follows. In section II, a theoretical overview of the LS equation, MST, the elastic scattering equation, and the optical potential are reviewed. This is followed by a discussion of the eikonal, PW, and LS3D solution methods in section III. Comparisons of model results and experimental data are given in section IV. The conclusions are stated in section V.
Ii Theoretical Framework
The LS equation is an expression for the scattering transition operator —the fundamental quantity that is used to evaluate the elastic differential, elastic, reaction, and total cross sections for nuclear reactions—and is given as
(1) |
where V is the sum of residual two-body interactions for the projectile-target system, and is the unperturbed two-body propagator Maung et al. (2007). Using projection operators, Eq. 1 can be expressed as a coupled system of equations Feshbach (1958, 1962)
(2) | ||||
(3) |
where Eq. 2 is the elastic scattering equation, and is the optical potential. The ground state projector is defined as , and the excited state projectors are defined , where is the projectile wave-vector, and is the target wave-vector.
In the non-relativistic multiple scattering theory (MST), the free Hamiltonian can be separated from the residual interaction, . If the interaction is expressed as the sum of two-body projectile and target nucleon interactions, , then the Watson series for the optical potential is given by Watson (1953); Werneth and Maung (2013),
(4) |
with
(5) |
where is the number of nucleons in the projectile () or target (), and are the Watson- operators that are expressed as . The Watson- operators are often approximated by the free two-body transition amplitudes (impulse approximation) given by
(6) |
where is the free NN Green’s function. The current work uses the first order (single scattering) approximation for the optical potential, which is given by
(7) |
Note that even in the first order approximation, represents an infinite series in terms of —which can be seen by iteration of equation (6)—but, in practice, is parameterized to experimental data.
The elastic scattering equation is written
(8) |
where is the initial (final) momentum in the center of momentum (CM) frame, is the energy, is the relative on-shell momentum, and is imposed to ensure outward scattering boundary conditions.
The optical potential is found by taking the matrix element of equation (7),
(9) |
where using the Watson Watson (1953) convention. Following the work in references Picklesimer et al. (1984a); Wolfe (1983); Elster et al. (1997b), the optical potential can be expressed as
(10) |
where , , is the nuclear matter density, is the NN transition amplitude, is the NN CM energy, and is the Möller frame transformation factor Möller (1945); Joachain (1983) used to transform from the AA to NN CM frame. Nuclear charge densities and the NN transition amplitude are parameterized to experimental data Werneth et al. (); Werneth et al. (2013); Wilson et al. (1991); De Jager et al. (1974); De Vries et al. (1987); Townsend et al. (1982); Townsend (1983).
The elastic scattering amplitude is related to the transition matrix by Joachain (1983)
(11) |
where , is the CM scattering angle, the density of states, , is given by
(12) |
and is the energy. For non-relativistic (NR) kinematics, , where is the reduced mass, is the mass of the projectile, and is the mass of the target. When using relativistic (REL) kinematics, . Elastic differential cross sections are determined from the scattering amplitude by using
(13) |
Iii Solution Methods
The Lippmann Schwinger equation was solved with two approximate methods and a full three-dimensional approach. Approximate solutions include the eikonal method, which employs a forward scattering approximation, and the PW method in the which the transition amplitude is expanded in an infinite series of Legendre polynomials. This section outlines the solution methods and numerical techniques used to solve for the transition matrix and scattering amplitude.
iii.1 Lippmann-Schwinger Partial Wave Solution Method
The LS equation is often solved with partial wave decomposition, a well-known method that is described in standard texts Joachain (1983); Landau (1996); Werneth et al. (2013). In this method, the transition matrix is decomposed into a complete orthonormal set of momenta dependent functions and angular dependent Legendre polynomials. For elastic scattering
(14) |
where , , are the Legendre Polynomials, , is the angle between and , and . The angular dependence is integrated, and the solution to the Lippmann-Schwinger equation is found for each partial wave,
(15) |
where
(16) |
Equation (15) is expressed in terms of its principal value integral, and Gaussian quadrature is used for the momentum integration variable. Sloan’s method Sloan (1968) is employed for the principal value integral, and the transition amplitude is expressed as a matrix equation for each partial wave, which is solved. The number of partial waves needed for an acceptable tolerance of convergence is not known a priori. Partial waves must be generated until such a tolerance is reached.
In the current work, the authors use a finite summation formula for the transition amplitude, which is given by Werneth et al. (2013)
(17) |
where represents a finite angular momentum that is reached when according to a pre-defined tolerance of %.
iii.2 Lippmann-Schwinger 3D Solution Method
The three-dimensional Lippmann-Schwinger (LS3D) solution method avoids the numerical difficulties associated with the PW method and has been used for relatively low energy reactions Elster et al. (1997a); Fachruddin et al. (2000); Rodriguez-Gallardo et al. (2008); Veerasamy et al. (2013); Liu et al. (2005). This section outlines the LS3D equation and the solution methods.
If one considers only central potentials in equation (8), then both and are scalar functions; that is, for some function , where represents the unit vector associated with . The possible scalar products of the LS equation are as follows Elster et al. (1997a); Rodriguez-Gallardo et al. (2008):
(18) | |||||
The incoming momentum, , is taken to be in the direction of the -axis, and the azimuthal angle between and is set to zero: ; therefore, may be expressed as a function of , , and Elster et al. (1997a); Rodriguez-Gallardo et al. (2008),
(19) |
and the LS3D equation is given by Elster et al. (1997a)
(20) | ||||
where Elster et al. (1997a)
(21) |
The numerical implementation of the LS3D method proceeds in the same manner as the PW method, but there are now two additional integration variables over azimuthal and polar angles. The azimuthal dependence only occurs in the potential and is integrated with 40 Gaussian quadrature points. As was seen with the PW method, the principal value integral over momenta is handled with Sloan’s method Sloan (1968), and the transition amplitude is expressed as a matrix equation, which is solved. The solution, , corresponds to the transformed Gaussian quadrature points associated with the integral. These results are substituted back into equation (20) to obtain the transition amplitude at the specified final momentum () and angle ().
It has been observed that the transition amplitude for reactions with large on-shell momenta—including high energy NA reactions and AA reactions at every energy—do not converge efficiently if the integration ranges of both momenta and polar angles are not restricted to regions that give significant contributions to the LS equation. The momenta which give non-zero contributions are estimated from the range of the optical potential and tend to be localized near the on-shell momentum, . For numerical efficiency and to ensure convergence, the integrations are truncated accordingly. The number of Gaussian quadrature points for the LS3D solution method was increased to a maximum of 44 points such that the total elastic cross sections changed less than 1% for all reactions with energies up to 100 GeV/n.
iii.3 Eikonal Solution Method
The eikonal approximation is used for high energy, small angle scattering to calculate elastic, reaction, total, and elastic differential cross sections Joachain (1983); Glauber (1959). To compute cross sections with the eikonal method, one solves for the eikonal scattering amplitude, , which is given as Joachain (1983)
(22) |
where is the relative momentum of the projectile-target system in the CM frame, is the ordinary cylindrical Bessel function, is the scattering angle in the CM frame, is the impact parameter, and is the eikonal phase shift function, the latter of which is obtained by integrating over the optical potential, Joachain (1983):
(23) |
The -integration is taken to be in the same direction as the initial wave vector of the incident projectile. The optical potential in equation (23) is the Fourier transform of equation (10) Wilson et al. (1991); Townsend et al. (1983)
The numerical evaluation of equation (23) is inefficient when the six-dimensional position-space integral of the optical potential is solved. In the present work, formulas of the optical potential are used for light nuclei () Werneth et al. (), and the eikonal phase function is written in the momentum-space representation for heavier nuclei Werneth et al. (),
(24) |
The advantage of equation (24) is that the optical potential is in the momentum-space representation, and the -integration need not be performed. Instead, the 7-dimensional integral for has been reduced to 2-dimensions over the magnitude of the momentum transfer, , and the angle, , between the momentum transfer and the impact parameter. This result significantly increases the efficiency for the numerical evaluation of .
Although the momentum space method for the eikonal phase function is much more efficient than the position space calculation, additional interpolation over the impact parameter and momentum transfer was performed for additional numerical efficiency. Convergence of the total elastic cross sections was used to establish the number of Gaussian quadrature points used for integration. The number of Gaussian quadrature points for the eikonal solution method was increased to a maximum of 100 points such that the total elastic cross sections changed less than 1% for all reactions with energies up to 100 GeV/n.
Iv Results
In the results that follow, each model uses the same set of fundamental parameterizations for the nuclear matter densities and the NN transition amplitude. Harmonic well and two-parameter Fermi (Wood-Saxon) nuclear charge data are taken from references De Jager et al. (1974); De Vries et al. (1987) and are normalized to matter densities as described in reference Wilson et al. (1991). When data are not available for the two-parameter Fermi densities, a nuclear droplet model Meyers and Swiatecki (1974) is used for parameter estimates. Nuclei are assumed to be near the beta stability curve. The NN transition amplitude used in the current work is described in reference Werneth et al. () and depends on parameterizations of the NN cross sections, slope parameter, and real to imaginary ratio of the transition amplitude. The NN cross sections are taken from reference Norbury (), and the slope parameter is from reference Werneth et al. ().
In Figs. 1-4, NA and AA elastic differential cross sections are shown at energies that are relevant to space radiation applications, including p + O, p + Fe, He + O, and C + Fe reactions at lab projectile kinetic energies of 150, 500, 1000, and 20,000 MeV/n. Results are indicated non-relativistic by (NR) and relativistic by (REL). LS3D (REL) results are given as a solid red line; a dashed, black line is used for eikonal results, denoted (Eik); a green square represents the PW (NR) results; a solid blue circle indicates LS3D (NR) results; and a violet asterisk is for PW (REL) results. Note that the Coulomb interaction has not been included in this analysis.
Excellent agreement between PW and LS3D results are seen in Figs. 1-4 for each kinematic selection for energies greater than 150 MeV. The p + O and p + Fe reactions at 150 MeV/n in Figs. 1 and 2 show slight disagreements between NR PW and LS3D codes and eikonal results. This is likely the result of the forward scattering approximation used in the eikonal method, since very light projectiles may deviate from forward scattering at low energy. The slight disagreement between the eikonal and NR PW and LS3D codes is not observed for the heavier nuclei in Figs. 3 and 4, where the small angle scattering approximation is more appropriate.
The next obvious feature is that of the relativistic shift observed in Figs. 1-4. The magnitude of the differential cross section is larger at smaller angles as compared to the NR cases. The effect is more pronounced at higher energies, as expected, but is also driven by projectile and target mass differences. A comparison of Figs. 1 and 3 shows that the relativistic effect is more pronounced for the p + O reaction, which has larger mass difference than the He + O system. Ultimately, the relativistic effects can be tracked back to kinematic differences in the relative on-shell momentum.
As an example of the LS3D method and illustration of the relativistic shift, comparisons to experimental data Alkhazov et al. (1972); Hoffmann et al. (1981); Lombard et al. (1981); Alkhazov et al. (1977) are performed. Fig. 5 shows the elastic differential cross sections of the following reactions: (a) p + S at GeV Alkhazov et al. (1972) (b) p + Ca at MeV Hoffmann et al. (1981) (c) p + Ni at GeV Lombard et al. (1981) and (d) He + Ca at MeV/n Alkhazov et al. (1977). NR results are indicated with a solid red line, and REL results are indicated with a solid blue line. In each case, there is better agreement with experiment when relativistic kinematics are used. Since the fundamental parameterizations are based on small-angle scattering data, the results are in better agreement with the measured differential cross section data at forward scattering angles. Also note that spin-dependence and medium effects have not been included, which may account for differences between the model and the experimental data in Fig. 5 (d).
V Conclusions
The eikonal, PW, and LS3D methods have been compared for NA and AA reactions for reactions relevant for space radiation applications. Numerical convergence of the eikonal method is readily achieved when formulas of the optical potential are used for light nuclei () and the momentum-space representation of the optical potential is used for heavier nuclei Werneth et al. (). The LS formalism has an advantage over the eikonal method in that relativistic kinematics are easily included.
The PW solution method is numerically unstable for reactions that have large on-shell momenta, including both high energy reactions and relatively low energy AA reactions, due to the highly oscillatory Legendre polynomials needed for convergence of these systems. To circumvent this difficulty, the LS3D solution method was implemented. Convergence of the LS3D equation can be achieved quickly after identifying the integration range for momenta and polar angles that give non-zero contributions to the LS equation. This numerical method is also useful for obtaining convergence for the partial wave analysis; however, numerical instability still exists because of the Legendre polynomial oscillations.
It was shown that the NR PW and NR LS3D methods agree with the eikonal method, except at very low energies for projectile nucleons, where the eikonal method is not well-suited. As the lab energy is increased, relativistic effects are seen as a shift in differential cross section resonances toward higher magnitudes and lower angles. Although some comparisons to experimental data were performed, the aim of this was manuscript was to demonstrate that (1) all three models agree in the appropriate energy regimes (2) there is a noticeable shift in the elastic differential cross section when relativistic kinematics are used (3) the LS3D method can be used for high energy reactions, where PW methods are numerically unstable.
Based on the results presented herein, it is recommended that the LS3D method be used for high energy NA and AA reactions at all energies because of its rapid numerical convergence and stability. The effect of equal mass kinematics on differential cross sections and extensive comparisons to experimental data will be elucidated in subsequent manuscripts.
Acknowledgments
The authors would like to thank Drs. Steve Blattnig, Ryan Norman, Jonathan Ransom, and Francis Badavi for reviewing this manuscript. This work was supported by the Human Research Program under the Human Exploration and Operations Mission Directorate of NASA and NASA grant number NNX13AH31A. Khin Maung Maung would like to thank Alexander Maung for his helpful conversations.
References
- M. Ackermann et al., Science 339, 807 (2013).
- E. R. Benton and E. V. Benton, Nucl. Instr. Meth. B 184, 255 (2001).
- T. C. Slaba, S. R. Blattnig, and F. F. Badavi, J. Comput. Phys. 229, 9397 (2010a).
- T. C. Slaba, S. R. Blattnig, M. S. Clowdsley, S. A. Walker, and F. F. Badavi, Adv. Space Res. 46, 800 (2010b).
- J. W. Wilson, L. W. Townsend, W. Schimmerling, G. S. Khandelwal, F. Khan, J. E. Nealy, F. A. Cucinotta, L. C. Simonsen, J. L. Shinn, and J. W. Norbury, Transport methods and interactions for space radiations, NASA Reference Publication 1257 (1991).
- J. A. Simpson, Ann. Rev. Nucl. Part. Sci. 33, 323 (1983).
- C. J. Joachain, Quantum Collision Theory (American Elsevier, New York, 1983).
- A. Picklesimer, P. C. Tandy, R. Thaler, and D. Wolfe, Phys. Rev. C 30, 1861 (1984a).
- A. Picklesimer, P. C. Tandy, R. Thaler, and D. Wolfe, Phys. Rev. C 29, 1582 (1984b).
- D. H. Wolfe, Ph.D. thesis, Kent State University (1983).
- C. W. De Jager, H. De Vries, and C. De Vries, Atom. Data Nucl. Data 14, 479 (1974).
- H. De Vries, C. W. De Jager, and C. De Vries, Atom. Data Nucl. Data 36, 495 (1987).
- M. L. Miller, K. Reygers, S. J. Sanders, and P. Steinberg, Annu. Rev. Nucl. Part. Sci. 57, 205 (2007).
- L. W. Townsend, H. B. Bidasaria, and J. W. Wilson, Can. J. Phys. 61, 867 (1983).
- J. W. Wilson, Phys. Lett. B 52, 419 (1974).
- J. W. Wilson and L. W. Townsend, Can. J. Phys. 59, 1569 (1981).
- L. W. Townsend, J. W. Wilson, and H. B. Bidasaria, Can. J. Phys. 60, 1514 (1982).
- L. W. Townsend, Can. J. Phys. 61, 93 (1983).
- C. M. Werneth, K. M. Maung, W. P. Ford, J. W. Norbury, and M. D. Vera, Elastic Differential Cross Sections, NASA Technical Publication 2014 (submitted).
- C. M. Werneth, K. M. Maung, L. R. Mead, and S. R. Blattnig, Nucl. Instr. Meth. B 308, 40 (2013).
- C. Elster, J. H. Thomas, and W. Glockle, arXiv:9708017v1 [nucl-th] (1997a).
- I. Fachruddin, C. Elster, and W. Glockle, Phys. Rev. C 62, 044002 (2000).
- M. Rodriguez-Gallardo, A. Deltuva, E. Cravo, R. Crespo, and A. C. Fonseca, Phys. Rev. C 78, 034602 (2008).
- S. Veerasamy, C. Elster, and W. N. Polyzou, Few Body Systems 54, 2207 (2013).
- H. Liu, C. Elster, and W. Glockle, Phys. Rev. C 72, 054003 (2005).
- K. M. Maung, J. W. Norbury, and T. Coleman, J. Phys. G: Nucl. Part. Phys. 34, 1861 (2007).
- H. Feshbach, Ann. Phys. 5, 357 (1958).
- H. Feshbach, Ann. Phys. 19, 287 (1962).
- K. M. Watson, Phys. Rev. 89, 575 (1953).
- C. M. Werneth and K. M. Maung, Can. J. Phys. 91, 424 (2013).
- C. Elster, S. P. Weppner, and C. R. Chinn, Phys. Rev. C 56, 2080 (1997b).
- C. Möller, K. Dan. Vidensk. Selsk. Mat.-Fys. Medd. 23, 1 (1945).
- R. H. Landau, Quantum Mechanics II: A Second Course in Quantum Theory (John Wiley & Sons Inc., New York, 1996).
- I. H. Sloan, J. Comput. Phys. 3, 332 (1968).
- R. J. Glauber, Lectures in Theoretical Physics (Interscience Publishers Inc., New York, 1959).
- W. D. Meyers and W. J. Swiatecki, Ann. Phys. 84, 186 (1974).
- J. W. Norbury, Total Nucleon-Nucleon Cross Section, NASA Technical Publication 2008-215116.
- G. D. Alkhazov, G. M. Amalsky, S. L. Belostotsky, A. A. Vorobyov, O. A. Domchenkov, and Y. V. Dotsenko, Phys. Lett. B 42, 121 (1972).
- G. W. Hoffmann et al., Phys. Rev. Lett. 47, 1436 (1981).
- R. M. Lombard, G. D. Alkhazov, and O. A. Domchenkov, Nucl. Phys. A 360, 233 (1981).
- G. D. Alkhazov et al., Nucl. Phys. A 280, 365 (1977).