Titus 2015 378
Titus 2015 378
By
Luke Titus
A DISSERTATION
Submitted to
Michigan State University
in partial fulfillment of the requirements
for the degree of
2016
ABSTRACT
By
Luke Titus
Nuclear reactions play a key role in the study of nuclei away from stability. Single-
nucleon transfer reactions involving deuterons provide an exceptional tool to study the single-
particle structure of nuclei. Theoretically, these reactions are attractive as they can be
cast into a three-body problem composed of a neutron, proton, and the target nucleus.
nucleus optical potentials are made local for convenience. The effects of nonlocal potentials
have historically been included approximately by applying a correction factor to the solution
of the corresponding equation for the local equivalent interaction. This is usually referred
to as the Perey correction factor. In this thesis, we have systematically investigated the
effects of nonlocality on (p, d) and (d, p) transfer reactions, and the validity of the Perey
correction factor. We implemented a method to solve the single channel nonlocal equation
for both bound and scattering states. We also developed an improved formalism for nonlocal
interactions that includes deuteron breakup in transfer reactions. This new formalism, the
nonlocal adiabatic distorted wave approximation, was used to study the effects of including
For the (p, d) transfer reactions, we solved the nonlocal scattering and bound state equa-
tions using the Perey-Buck type interaction, and compared to local equivalent calculations.
Using the distorted wave Born approximation we construct the T-matrix for (p, d) transfer
on 17 O, 41 Ca, 49 Ca, 127 Sn, 133 Sn, and 209 Pb at 20 and 50 MeV. Additionally we studied
(p, d) reactions on 40 Ca using the the nonlocal dispersive optical model. We have also in-
cluded nonlocality consistently into the adiabatic distorted wave approximation and have
investigated the effects of nonlocality on on (d, p) transfer reactions for deuterons impinged
on 16 O, 40 Ca, 48 Ca, 126 Sn, 132 Sn, 208 Pb at 10, 20, and 50 MeV.
We found that for bound states the Perry corrected wave functions resulting from the
local equation agreed well with that from the nonlocal equation in the interior region, but
discrepancies were found in the surface and peripheral regions. Overall, the Perey correc-
tion factor was adequate for scattering states, with the exception for a few partial waves.
Nonlocality in the proton scattering state reduced the amplitude of the wave function in the
nuclear interior. The same was seen for nonlocality in the deuteron scattering state, but the
wave function was also shifted outward. In distorted wave Born approximation studies of
(p, d) reactions using the Perey-Buck potential, we found that transfer distributions at the
first peak differed by 15 − 35% as compared to the distribution resulting from local poten-
tials. When using the dispersive optical model, this discrepancies grew to ≈ 30−50%. When
nonlocality was included consistently within the adiabatic distorted wave approximation, the
If only local optical potentials are used in the analysis of experimental (p, d) or (d, p)
transfer cross sections, the extracted spectroscopic factors may be incorrect by up to 50% in
some cases due to the local approximation. This highlights the necessity to pursue reaction
I would like to thank my advisor Prof. Filomena Nunes whose support and guidance has
helped me both professionally and personally. I am grateful for her continued support
other fields and careers. Her consistent positive attitude is a model for me to follow, and
has made the difficult task of completing my degree a pleasure. Undoubtedly I would not
I would like to express my thanks to my guidance committee, Prof. Mark Voit, Prof.
Phil Duxbury, Prof. Morten Hjorth-Jensen, and Prof. Remco Zegers for their invaluable
advice and suggestions along the way. I’d like to thank again Prof. Mark Voit for taking me
under his wing in astronomy for a semester, and for the countless hours of discussions we
had about galaxy clusters and the universe. I would also like to thank Prof. Ian Thompson,
Prof. Ron Johnson, Prof. Pierre Capel, Prof. Jeff Tostevin, and Dr. Gregory Potel, without
I owe a debt of gratitude to the Department of Physics and Astronomy at Michigan State
University, the National Superconducting Cyclotron Laboratory, and the theory group for
providing financial, academic, and technical support. I would also like thank the National
Science Foundation and the Department of Energy for their financial support.
I thank my colleagues and friends at Michigan State University, and especially my current
and past group members: Bich Nguyen, Neelam Upadhyay, Muslema Pervin, Amy Lovell,
Alaina Ross, Terri Poxon-Pearson, Gregory Potel, Jimmy Rotureau, and Ivan Brida.
Last but not least, I thank my parents, Carrie Gossen and Joe Titus, my brother Blake
Titus, and the rest of my family for love and support throughout my entire life, and their
iv
constant encouragement as I developed my interest in science. Without them, I would never
v
TABLE OF CONTENTS
LIST OF FIGURES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . x
Chapter 1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1 Nuclear Interactions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.2 Nonlocality . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.2.1 Microscopic Optical Potentials . . . . . . . . . . . . . . . . . . . . . . 10
1.2.2 Phenomenological Nonlocal Optical Potentials . . . . . . . . . . . . . 13
1.2.3 Solving Nonlocal Equations . . . . . . . . . . . . . . . . . . . . . . . 14
1.3 Motivation for present work . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.4 Outline . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
Chapter 4 Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
4.1 Numerical Details . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
4.1.1 Effects of Neglecting Remnant . . . . . . . . . . . . . . . . . . . . . . 70
4.2 Distorted Wave Born Approximation with the Perey-Buck Potential . . . . . 70
4.2.1 Proton Scattering State . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.2.2 Neutron Bound State . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
4.2.3 (p, d) Transfer Cross Sections - Distorted Wave Born Approximation . 74
4.2.4 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
vi
4.3 Distorted Wave Born Approximation with the Dispersive Optical Model Po-
tential . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
4.3.1 Proton Scattering State . . . . . . . . . . . . . . . . . . . . . . . . . 81
4.3.2 Neutron Bound State . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
4.3.3 (p, d) Transfer Cross Sections - Distorted Wave Born Approximation . 84
4.3.4 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
4.4 Nonlocal Adiabatic Distorted Wave Approximation with the Perey-Buck Po-
tential . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
4.4.1 The Source Term . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
4.4.2 Deuteron Scattering State . . . . . . . . . . . . . . . . . . . . . . . . 90
4.4.3 (d, p) Transfer Cross Sections . . . . . . . . . . . . . . . . . . . . . . 93
4.4.4 Comparing Distorted Wave Born Approximation and the Adiabatic
Distorted Wave Approximation . . . . . . . . . . . . . . . . . . . . . 99
4.4.5 Energy Shift Method . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
4.4.6 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
APPENDICES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
Appendix A Solving the Nonlocal Equation . . . . . . . . . . . . . . . . . . . . . 111
Appendix B Deriving the Perey Correction Factor . . . . . . . . . . . . . . . . . 116
Appendix C Nonlocal Adiabatic Potential . . . . . . . . . . . . . . . . . . . . . . 123
Appendix D Deriving the T-Matrix . . . . . . . . . . . . . . . . . . . . . . . . . 141
Appendix E Checks of the Code NLAT . . . . . . . . . . . . . . . . . . . . . . . 176
Appendix F Mirror Symmetry of ANCs . . . . . . . . . . . . . . . . . . . . . . . 186
Appendix G List of Acronyms . . . . . . . . . . . . . . . . . . . . . . . . . . . . 196
BIBLIOGRAPHY . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 197
vii
LIST OF TABLES
Table 3.1: Potential parameters for the Perey-Buck [1] and TPM [2] nonlocal
potentials. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
Table 4.1: Percent difference of the (d, p) transfer cross section at the first peak
for a calculation including the remnant term relative to a calculation
without the remnant term. . . . . . . . . . . . . . . . . . . . . . . . 70
Table 4.2: Percent difference of the (p, d) transfer cross sections at the first peak
when using the PCF (2nd column), or a nonlocal potential (3rd col-
umn), relative to the local calculation with the LPE potential, for a
number of reactions occurring at 20 MeV. . . . . . . . . . . . . . . . 79
Table 4.3: Percent difference of the (p, d) transfer cross sections at the first peak
when using the PCF (2nd column), or a nonlocal potential (3rd col-
umn), relative to the local calculation with the LPE potential, for a
number of reactions occurring at 50 MeV. . . . . . . . . . . . . . . . 79
Table 4.4: Percent differences of the (p, d) transfer cross sections at the first
peak at the listed beam energies using the DOM potential relative
to the calculations with the phase-equivalent potential. Results are
listed separately for the effects of nonlocality on the bound state, the
scattering state, and the total. . . . . . . . . . . . . . . . . . . . . . 85
Table 4.5: Same as Table 4.4, but now for the Perey-Buck potential. . . . . . . 85
Table 4.6: Percent difference of the (d, p) transfer cross sections at the first peak
when using nonlocal potentials in entrance and exit channels (1st
column), nonlocal potentials in entrance channel only (2nd column),
and nonlocal potentials in exit channel only (3rd column), relative
to the local calculation with the LPE potentials, for a number of
reactions occurring at 10 MeV. . . . . . . . . . . . . . . . . . . . . . 94
Table 4.7: Percent difference of the (d, p) transfer cross sections at the first peak
when using nonlocal potentials in entrance and exit channels (1st
column), nonlocal potentials in entrance channel only (2nd column),
and nonlocal potentials in exit channel only (3rd column), relative
to the local calculation with the LPE potentials, for a number of
reactions occurring at 20 MeV. . . . . . . . . . . . . . . . . . . . . . 97
viii
Table 4.8: Percent difference of the (d, p) transfer cross sections at the first peak
when using nonlocal potentials in entrance and exit channels (1st col-
umn), nonlocal potentials in entrance channel only (2nd column), and
nonlocal potentials in exit channel only (3rd column), relative to the
local calculation with the LPE potentials, for a number of reactions
occurring at 50 MeV. Figure reprinted from [3] with permission. . . 98
Table E.1: The nonlocal adiabatic integral, rhs of Eq.(2.49), calculated with
Mathematica and NLAT using analytic expressions for the wave func-
tions and potentials. . . . . . . . . . . . . . . . . . . . . . . . . . . . 183
Table F.1: Ratio of proton to neutron ANCs for the dominant component: Com-
parison of this work R with the results of the analytic formula R0
Eq.(F.10) and the results of the microscopic two-cluster calculations
RM CM [4, 5] including the Minnesota interaction. The uncertainty
in R account for the sensitivity to the parameters of VBx . . . . . . 191
ix
LIST OF FIGURES
Figure 1.1: The chart of the nuclides. The proton drip line is indicated by the line
above the stable nuclei, and the neutron drip line is indicated below
the stable nuclei. The proton (neutron) drip line indicates where the
addition of a single proton (neutron) will make the resulting nucleus
unbound. Figure reprinted from [6] with permission. . . . . . . . . . 2
Figure 1.4: Angular distributions for elastic scattering of nucleons off 208 Pb at
25 MeV. (a) n+208 Pb (b) p+208 Pb with differential cross section
normalized to Rutherford. . . . . . . . . . . . . . . . . . . . . . . . 8
Figure 2.1: First, second, and all-order couplings within a set of 0+ , 2+ , and 4+
nuclear levels, starting from the ground state. . . . . . . . . . . . . 28
Figure 2.3: The coordinates used to calculate the T-matrix for (d, p) transfer. . 34
Figure 2.4: The first four Weinberg States when using a central Gaussian which
reproduces the binding energy and radius of the deuteron ground
state. The inset shows the asymptotic properties of each state. . . . 40
Figure 2.5: The coordinates used for constructing the neutron nonlocal potential.
The open dashed circle represents the neutron in a different point in
space to account for nonlocality. . . . . . . . . . . . . . . . . . . . . 42
Figure 2.6: The coordinates used for constructing the system wave function for
the d + A deuteron scattering state. . . . . . . . . . . . . . . . . . . 44
x
Figure 3.1: Differential elastic scattering relative to Rutherford as a function of
scattering angle. (a) 48 Ca(p, p)48 Ca at 15.63 MeV with data from [9]
(b) 208 Pb(p, p)208 Pb at 61.4 MeV with data from [10]. . . . . . . . . 57
Figure 3.3: 49 Ca(p, p)49 Caat 50.0 MeV: The solid line is obtained from using
the Perey-Buck nonlocal potential, the open circles are a fit to the
nonlocal solution, and the dotted line is obtained by transforming the
depths of the volume and surface potentials according to Eq. (B.14).
Figure reprinted from [12] with permission. . . . . . . . . . . . . . . 65
Figure 4.1: Real and imaginary parts of the J π = 1/2− partial wave of the scat-
tering wave function for the reaction 49 Ca(p, p)49 Ca at 50.0 MeV:
ψ N L (solid line), ψ P CF (crosses), and ψ loc (dashed line). Top (bot-
tom) panel: absolute value of the real (imaginary) part of the scat-
tering wave function. Figure reprinted from [12] with permission. . . 72
Figure 4.2: Real and imaginary parts of the J π = 11/2+ partial wave of the
scattering wave function for the reaction 49 Ca(p, p)49 Ca at 50.0 MeV.
See caption of Fig. 4.1. Figure reprinted from [12] with permission. . 72
Figure 4.3: Ground state, 2p3/2 , bound wave function for n+48 Ca. φN L (solid
line), φP CF (crosses), and φloc (dashed line). The inset shows the
difference φN L − φP CF . Figure reprinted from from [12] with per-
mission. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
Figure 4.4: Angular distributions for 49 Ca(p, d)48 Ca at 50 MeV: Inclusion of non-
locality in both the proton scattering state and the neutron bound
state (solid), using LPE potentials, then applying the correction fac-
tor to both the scattering and bound states (crosses), using the LPE
potentials without applying any corrections (dashed line), including
nonlocality only in the proton scattering state (dotted line) and in-
cluding nonlocality only in the neutron bound state (dot-dashed line).
Figure reprinted from [12] with permission. . . . . . . . . . . . . . . 75
Figure 4.5: Same as in Fig. 4.4 but for 49 Ca(p, d)48 Ca at Ep = 20 MeV. Figure
reprinted from [12] with permission. . . . . . . . . . . . . . . . . . . 77
Figure 4.6: Same as in Fig. 4.4 but for 133 Sn(p, d)132 Sn at Ep = 20 MeV. Figure
reprinted from [12] with permission. . . . . . . . . . . . . . . . . . . 77
xi
Figure 4.7: Same as in Fig. 4.4 but for 209 Pb(p, d)208 Pb at Ep = 20 MeV. Figure
reprinted from [12] with permission. . . . . . . . . . . . . . . . . . . 78
Figure 4.9: Same as in Fig. 4.8 but for Ep = 50 MeV. Data from [15]. Figure
reprinted from [14] with permission. . . . . . . . . . . . . . . . . . . 82
Figure 4.10: The real and imaginary parts of the J π = 1/2+ partial wave of the
scattering wave function for the reaction 40 Ca(p, p)40 Ca at Ep = 50
MeV. This shows the wave function resulting from the DOM poten-
tial (solid line) and its LPE potential (dotted line), the Perey-Buck
potential (dashed line) and its LPE potential (dot-dashed line). The
top (bottom) panel shows the absolute value of the real (imaginary)
part of the scattering wave function. Figure reprinted from [14] with
permission. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
Figure 4.11: The neutron ground state 1d3/2 bound wave function for n+39 Ca.
Shown is the wave function obtained using the DOM potential (solid
line), the Perey-Buck potential (crosses) and the local interaction
(dashed line). The inset shows the asymptotic properties of each
wave function. Figure reprinted from [14] with permission. . . . . . 84
Figure 4.12: Angular distributions for the 40 Ca(p, d)39 Ca reaction at (a) Ep = 20
MeV, (b) Ep = 35 MeV, and (c) Ep = 50 MeV. In this figure is the
transfer distribution resulting from using the nonlocal DOM (solid
line) and its LPE potential (dotted line), the Perey-Buck potential
(dashed line) and the Perey-Buck LPE potential (dot-dashed line).
Figure reprinted from [14] with permission. . . . . . . . . . . . . . . 86
Figure 4.13: Absolutes value of the d + A source term when nonlocal and local
potentials are used. (a) d+48 Ca at Ed = 50 MeV. (b) d+208 Pb at
Ed = 50 MeV. Both are for the L = 1 and J = 0 partial wave. Figure
reprinted from [3] with permission. . . . . . . . . . . . . . . . . . . . 89
Figure 4.14: Absolute value of the d + A source term when nonlocal and local
potentials are used. (a) d+48 Ca at Ed = 50 MeV. (b) d+208 Pb at
Ed = 50 MeV. Both are for the L = 6 and J = 5 partial wave. Figure
reprinted from [3] with permission. . . . . . . . . . . . . . . . . . . . 90
xii
Figure 4.15: Absolute value of the d+A scattering wave function using the ADWA
theory when nonlocal and local potentials are used. (a) d+48 Ca and
(b) d+208 Pb. Both for the L = 1 and J = 0 partial wave at Ed = 50
MeV in the laboratory frame. Figure reprinted from [3] with permis-
sion. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
Figure 4.16: Absolute value of the d+A scattering wave function using the ADWA
theory when nonlocal and local potentials are used. (a) d+48 Ca and
(b) d+208 Pb. Both for the L = 6 and J = 5 partial wave at Ed = 50
MeV in the laboratory frame. Figure reprinted from [3] with permis-
sion. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
Figure 4.17: Angular distributions for (d, p) transfer cross sections. The insets
are the theoretical distributions normalized to the peak of the data
distribution. (a) 48 Ca(d, p)49 Ca at Ed = 10 MeV with data [16] at
Ed = 10 MeV in arbitrary units. (b) 132 Sn(d, p)133 Sn at Ed = 10
MeV with data [17] at Ed = 9.4 MeV. (c) 208 Pb(d, p)209 Pb at 20
MeV with data [18] (Circles) and [19] (Squares) at Ed = 22 MeV.
Figure reprinted from [3] with permission. . . . . . . . . . . . . . . 95
Figure 4.18: Angular distributions for (d, p) transfer cross sections. The inset is
the theoretical distributions normalized to the peak of the data distri-
bution. (a) 48 Ca(d, p)49 Ca at Ed = 50 MeV with data [20] at Ed = 56
MeV. (b) 132 Sn(d, p)133 Sn at Ed = 50 MeV. (c) 208 Pb(d, p)209 Pb at
50 MeV. Figure reprinted from [3] with permission. . . . . . . . . . 96
Figure 4.19: Comparison of (d, p) transfer cross sections when using the DWBA as
compared to the ADWA. (a) 16 O(d, p)17 O, (b) 48 Ca(d, p)49 Ca with
data from [20]. (c) 132 Sn(d, p)133 Sn. All distributions at Ed = 50
MeV. Figure reprinted from [3] with permission. . . . . . . . . . . . 100
Figure 4.20: Comparison of (d, p) angular distributions when using the energy
shift method of [21, 22]. (a) 16 O(d, p)17 O at Ed = 10 MeV (b)
40 Ca(d, p)41 Ca at E = 10 MeV (c) 208 Pb(d, p)209 Pb at E = 20
d d
MeV. The solid line is when full nonlocality was included in the en-
trance channel, dashed line is when the LPE potential was used, dot-
dashed line when the CH89 potential [23] was used with the additional
energy shift quantified in [21], and the dotted line when the CH89 po-
tential was used at the standard Ed /2 value. Figure reprinted from
[3] with permission. . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
xiii
Figure E.2: Differential elastic scattering as a function of scattering angle. 208 Pb(n, n)208 Pb
at Ep = 14.5 MeV: The solid line is obtained a nonlocal calculation
using NLAT, and the dashed line is the nonlocal calculation published
by Perey and Buck [1]. . . . . . . . . . . . . . . . . . . . . . . . . . 178
Figure E.3: n+48 Ca bound wave function, and the deuteron bound wave function.
n+48 Ca: The solid line is obtained from NLAT, the dotted line is
obtained from NLAT and setting β = 0.05 fm, and the dashed line is
from FRESCO. Deuteron: Dot-dashed line is deuteron bound wave
function obtained from NLAT, and the open circles are obtained with
FRESCO. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 179
Figure E.4: The local adiabatic potential for d+48 Ca at Ed = 20 MeV calculated
with NLAT and with TWOFNR [24]. (a) Real part, (b) Imaginary
part. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 180
Figure E.5: 48 Ca(d, d)48 Ca at Ed = 20 MeV. The solid line is when using the
local adiabatic potential, the dotted line is when doing a nonlocal
calculation with β = 0.1 fm in the nucleon optical potentials, and the
dashed line is a calculation done in FRESCO. . . . . . . . . . . . . 181
Figure E.7: Angular distributions for 208 Pb(d, p)209 Pb at Ed = 50 MeV obtained
by using different step sizes to calculate the rhs of Eq.(2.49). The
solid line uses a step size of 0.01 fm, the dashed line a step size of
0.03 fm, and the dotted line a step size of 0.05 fm. . . . . . . . . . . 184
Figure E.8: Angular distributions for 208 Pb(d, p)209 Pb at Ed = 50 MeV obtained
by using different step sizes and values of a cut parameter (CutL) to
calculate the rhs of Eq.(2.49). The solid line uses a step size of 0.01
fm with CutL=2, the dashed line a step size of 0.01 fm with CutL=3,
and the dotted line a step size of 0.05 fm with CutL=2. . . . . . . . 185
Figure F.1: Neutron and proton spectroscopic factors for 17 O and 17 F, respec-
tively, considering the 16 O core in its 0+ ground state and 2+ first
excited state: (a) 5/2+ ground state and (b) 1/2+ first excited state.
Figure reprinted from [25] with permission. . . . . . . . . . . . . . 193
Figure F.2: Ratio of proton and neutron ANCs for 17 O and 17 F, respectively, in-
cluding 16 O(0+ , 2+ ): (a) 5/2+ ground state and (b) 1/2+ first excited
state. Figure reprinted from [25] with permission. . . . . . . . . . . 193
xiv
Chapter 1
Introduction
Since the dawn of nuclear physics, reaction studies have been performed to investigate the
properties of the nucleus. One of the many reasons these studies have been carried out
is to address the overarching goal of nuclear physics. This is to understand where all the
matter in the universe came from and how it was formed. To solve this problem, we not only
need to understand the environments in which nuclear reactions occur, but we also need to
understand the nature of the nuclei undergoing the reactions. This is a daunting task with
hundreds of stable nuclei, and thousands of unstable nuclei known to exist [6].
In Fig. 1.1 the chart of the nuclides is shown with the corresponding proton and neutron
drip lines. The drip line is the point that separates bound from unbound nuclei. The neutron
drip line, for example, defines the point where the addition of a single neutron will make the
resulting nucleus unbound. While an extraordinary amount of progress has been made in
experimentally measuring unstable nuclei, it is remarkable how far the neutron drip line is
For many decades, intense experimental and theoretical effort has been put into studying
stable nuclei. While experiments aimed at studying stable isotopes are still performed,
the focus in modern times has shifted towards the study of exotic nuclei. In the context
of understanding the origin of the matter in the universe, exotic nuclei play a crucial role.
While exotic nuclei live for a very short period of time, reactions on exotic nuclei are essential
to creating heavy elements [8]. In certain astrophysical environments, nuclei rapidly capture
1
protons or neutrons, pushing them towards the drip line. These unstable nuclei then β decay
back to the valley of stability. To fully understand the path the nucleosynthesis takes, and
the elements that are produced, we must understand the properties of the exotic nuclei very
far from stability, and the reaction mechanisms of neutrons, protons, or heavier elements on
Figure 1.1: The chart of the nuclides. The proton drip line is indicated by the line above
the stable nuclei, and the neutron drip line is indicated below the stable nuclei. The proton
(neutron) drip line indicates where the addition of a single proton (neutron) will make the
resulting nucleus unbound. Figure reprinted from [6] with permission.
For many nuclear reaction experiments, a good reaction theory is required to extract
reliable information. The same can be said about the potentials we put into our theories. In
fact, the two work hand in hand. An excellent model can be held back by the use of poor
interactions, while the best interaction available will provide little insight when used in a
poor model.
An important part of understanding the properties of nuclei is knowing the spin and parity
assignments of the various energy levels. Single nucleon transfer reactions are an excellent
2
tool for understanding these properties. The protons and neutrons inside a nucleus arrange
themselves in an organized way, roughly following the way levels organize themselves in a
harmonic oscillator potential with a spin-orbit interaction. Filling a shell provides additional
stability. Indicated in Fig. 1.2 (right) are the magic numbers corresponding to the number
of neutrons or protons needed to fill in a shell. The ordering shown in Fig. 1.2 provides a
guide to assigning energy levels. As one moves away from stability there is shell reordering
The use of single nucleon transfer reactions such as (d, p) or (p, d) as a probe to study
nuclear structure began in the early 1950s. Butler realized that the spins and parities of
nuclear energy levels can be obtained from angular distributions, without the need to know
properties of excited states [26]. This fact was reiterated by Huby [27, 28], and later followed
up with theoretical calculations by Bhatia and collaborators [29]. While these early studies
relied on the very simple plane wave Born approximation, it drew considerable attention
to (d, p) reactions as a means to study nuclear structure through the analysis of angular
Since these pioneering studies, the shell structure of nuclei has been studied with the aid
of single nucleon transfer reactions. Of particular interest for this work are the stripping
(d, p) or pickup (p, d) reactions. These types of reactions are an excellent tool for measuring
the energy levels of nuclei, as well as the spin and parity assignments of the corresponding
energy levels. It is transfer reactions such as these which provided much of the structure
information of stable isotopes in the early days of nuclear physics [30, 31, 32, 33, 34, 35].
In Fig. 1.3 we show the dependence of the transfer angular distribution on the transferred
angular momentum for 58 Ni(d, p)59 Ni at 10 MeV. The transferred angular momentum has
an influence on the shape of the transfer distribution, as well as the location of the peak of
3
Figure 1.2: Typical Nuclear Shell Structure.
the transfer distribution. It is seen that for the s1/2 state the peak occurs at 0◦ , and for
increasing angular momentum transfer the first peak gets shifted to increasing angles. The
analogous to that of a single slit diffraction pattern. With increasing energy the diffraction
pattern is found to have more oscillations. Also, as the beam energy increases, the transfer
distribution gets shifted to more forward angles. The magnitude of the cross section is related
to the Q-value, or energy mismatch, of the reaction. The magnitude of the cross section is
Modern reaction theories have progressed greatly since the 1950s, allowing for more reli-
able nuclear structure information to be extracted from experimental data. The theoretical
4
Figure 1.3: Dependence of the transfer angular distribution on the transferred angular mo-
mentum for 58 Ni(d, p)59 Ni at 8 MeV, with data from [7]. Reprinted from [8] with permission.
advances of reaction theory coupled with advances in experimental techniques have made the
use of transfer reactions to study exotic nuclei feasible. In the early days of nuclear physics,
transfer reactions were performed by making a target using stable nuclei, and impinging
protons, deuterons, 3 He, or other nuclei on the target to initiate the transfer process. When
studying unstable nuclei, these reactions are done in inverse kinematics [17, 36, 37, 38, 39, 40].
Since exotic nuclei are too short lived to make a target, a deuterated target, for example,
is sometimes used, and a beam of exotic nuclei is impinged on the target to perform the
experiment.
As (d, p) or (d, n) transfer reactions are a useful tool for studying the overlap function of
the final nucleus, these reactions are also a preferred method to extract the normalization
5
of the tail of the overlap function. This quantity is known as the asymptotic normalization
coefficient (ANC), and is defined in Eq.(F.8). At very low energies, the transfer cross section
is dominated by the amplitude of the overlap function in the asymptotic region. Thus,
a (d, n) transfer reaction can provide information on the proton bound state of the final
nucleus. The ANC can be used to determine astrophysically important (p, γ) reaction rates
Making use of the ANC method, transfer reactions have also become a common tool
[42, 43]. Sometimes, the ANC for the system of interest is not accessible directly, while
the mirror system is. When this is the case, charge symmetry of the nuclear force can be
exploited to derive a model independent quantity relating the ratio of ANCs of the two
systems [44]. This has been shown to be a reliable method to indirectly extract an ANC
Whereas ANCs calculated theoretically can be very different depending on the model
that is used, the idea behind the method proposed in [44] suggests that the ratio of ANCs
of mirror pairs is model independent. This method is very useful to extract the ANC of the
proton state, useful in (p, γ) reactions important for astrophysics, by measuring the mirror
partner. In the early stage of my graduate work, we performed a study to test the model
independence of the ratio of ANCs of mirror pairs, and the validity of the analytic formula
6
1.1 Nuclear Interactions
problem. To solve the problem exactly would require the fully anti-symmetrized many-body
wave function that includes the couplings of the elastic channel to all the other non-elastic
channels available (transfer, inelastic scattering, charge exchange, fusion, fission, etc.). This
is a very difficult problem to solve, and in practice, the scattering process is not solved in this
manner. However, the elastic scattering of a particle from some arbitrary potential, U (R),
is well understood [8, 47]. Assuming that the complicated interaction between some particle
and the nucleus can be represented by a complex mean-field is the basis of the optical model.
In Fig. 1.4 we show the angular distributions for elastic scattering of nucleons off 208 Pb at
25 MeV. In panel (a) is n+208 Pb, and in panel (b) is p+208 Pb. Due to the Coulomb potential,
cross section, and always goes to unity at 0◦ . When this is done, the angular distributions
for proton elastic scattering are unitless. The oscillations result from a diffraction pattern
which can be understood qualitatively in a similar way as single slit diffraction. For a larger
target or lower energy, there will be fewer oscillations between 0◦ and 180◦ , and there will
In the optical model, elastic scattering data are fit by varying potential parameters in an
assumed form for U (R). This complex interaction, referred to as the optical potential, is used
to describe the elastic scattering process of the particle off the nucleus, with the imaginary
part taking into account loss of flux to non-elastic channels. Once the optical potential is
defined, it can then be used as an input to a model that describes some other process with the
goal of obtaining an observable other than elastic scattering, such as transfer cross sections.
7
5
10 1.4
1.2
4
10
1.0
dσ/dΩ [mb/sr]
3
10 0.8
σ/σR
2 0.6
10
0.4
1
10
0.2
a) b)
0
10 0 0.0
30 60 90 120 150 180 0 30 60 90 120 150 180
θc.m. θc.m.
Figure 1.4: Angular distributions for elastic scattering of nucleons off 208 Pb at 25 MeV. (a)
n+208 Pb (b) p+208 Pb with differential cross section normalized to Rutherford.
Elastic scattering data for the desired target and energy are often times not available. To
remedy this problem, optical potentials are constructed through simultaneous fits to large
data sets of elastic scattering. These are referred to as global optical potentials. The energy,
target, and projectile dependent parameters are varied to produce a best fit to the entire
data set. The purpose of using a global potential is that one can easily interpolate in order to
obtain a potential for a nucleus in which there is no experimental data available. Obtaining a
feature of using a global potential and is a credit to their success over the decades. It is for
this reason that considerable effort has been put into creating many different global optical
potentials over the years which have received widespread use [48, 23, 49].
Global potentials are a very useful tool for studying nuclear reactions and predicting
observables. However, the way they are constructed leaves out a considerable amount of
physics. Elastic scattering only constrains the normalization of the scattering wave function
outside the range of the interaction. It is not sensitive to the short-range properties of
the wave function. Therefore, the short-range physics is not constrained at all by elastic
8
scattering. Also, much of the elastic scattering data that exists is for stable nuclei. With
the increasing interest of the study of rare isotopes, the extrapolations to exotic nuclei may
not be reliable. It is for this reason that a more physically motivated form for the optical
All widely used global optical potentials are local. However, when derived from a many-
body theory, the resulting optical potential is nonlocal. The strong energy dependence of
global potentials is assumed to account for the nonlocality that is neglected. With increasing
the validity of the local assumption, and develop methods to incorporate nonlocal potentials
1.2 Nonlocality
It has long been known that the optical potential is nonlocal [50]. In the Hartree-Fock the-
ory, the existence of an exchange term introduces an explicit nonlocal potential [51]. For
scattering, the complicated coupling of the elastic channel to all other non-elastic channels
accounts for another significant source of nonlocality [52, 53]. These two sources of nonlo-
cality, anti-symmetrization and channel couplings, have been known and studied for decades
(e.g. [54]).
locate the center of the deuteron relative to the center of the target. Let’s say that the
deuteron breaks up at R0 as it approaches the target. The deuteron can then propagate
through space in its broken up state, then recombine to form the deuteron again at R.
This process is depicted in Fig. 1.5. Such a process would constitute a channel coupling
9
nonlocality. This will result in a potential of the form V (R, R0 ) since the interaction at a
given point is dependent on the value of the potential and the scattering wave function at
Figure 1.5: An example of a channel coupling nonlocality. In this case, the deuteron is im-
pinged on some target. The channel coupling nonlocality results from the deuteron breaking
up as it approaches the nucleus, propagating through space in its broken up state, and then
recombining to form the deuteron again.
As another example, consider a single nucleon scattering off a nucleus. Since the system
wave function is a fully anti-symmetric many-body wave function, it is not guaranteed that
the projectile in the incident channel is the same particle as the one in the exit. The Pauli
principle also plays a role when the projectile is propagating through the nuclear medium,
and most notably has the effect of reducing the amplitude of the wave function in the nuclear
interior. All of these effects will manifest in a potential of the form V (R, R0 ).
This work is not concerned with constructing a microscopic optical potential, but rather
with using current phenomenological nonlocal optical potentials, and studying the effects of
10
amount of effort that has been put forth in recent decades to construct optical potentials
from microscopic theories. In this thesis we will demonstrate that nonlocality is an important
feature of the nuclear potential that must be considered explicitly. Moving forward, the
development of ab-initio many-body theories offer the promise of realistic microscopic optical
potentials. The methods outlined here will be the tools for future studies.
Several studies have been made to construct a microscopically based optical potential. In
the pioneering work of Watson [55, 56], and later refined by Kerman, McManus, and Thaler
(KMT) [57], the theory of multiple scattering was developed, where the optical potential to
describe elastic scattering is constructed in terms of the amplitudes for the scattering of the
incident particle by the individual neutrons and protons in the target nucleus. This theory
for constructing the optical potential is limited to relatively high energies (> 100 MeV).
Deriving the multiple scattering expansion of the KMT optical potential is a complicated
The optical potential can also be identified with the self-energy, as first indicated by Bell
and Squires [50]. Hüfner and Mahaux studied the optical potential in great detail through
use of a systematic expansion of the self-energy within the Greens function approach to the
for the calculation of the binding energy [60]. This formulation of the optical potential in
terms of the self-energy is attractive as it is suitable for both intermediate and high energy
scattering, and reduces to the expressions of multiple scattering theory at high energies.
It is through the connection to the self-energy that Jeukenne, Lejeune, Mahaux (JLM)
formulated their optical model potential for infinite nuclear matter [61]. In infinite nuclear
matter, the concept of a projectile and target lose their meaning. Instead, a potential energy
and a lifetime for a quasiparticle state obtained by creating a particle or hole with momentum
11
k above the correlated ground state is defined. Later, the JLM approach was extended to
The link between the self-energy and the optical potential was further explored by Ma-
haux and Sartor [63]. This implementation is known as the dispersive optical model (DOM).
The advantage of this method is that it provides a link between nuclear reactions and nu-
clear structure through a dispersion relation. In recent years, a local version of the DOM
was introduced for Calcium isotopes [64], and a nonlocal DOM was subsequently developed
for 40 Ca [65, 11]. Transfer reaction studies have shown that the local DOM is able to de-
scribe transfer cross sections as well as or better than global potentials [66], and that the
nonlocal DOM can significantly modify the shell occupancy, or spectroscopic factor, of the
Various other techniques exist which construct an optical potential through the self-
energy using modern advances in nuclear theory. Making use of the progress that has been
achieved, Holt and collaborators constructed a microscopic optical potential from the self-
energy for nucleons in a medium of infinite isospin-symmetric nuclear matter within the
The two sources of nonlocality, channel coupling and anti-symmetrization, have been
studied over the years by numerous authors [54, 68, 69] to name only a few. Many of
these studies derive the nonlocal potential using some microscopic theory, then compare
the potential obtained to commonly used phenomenological nonlocal potentials. Such was
done in [54] where the multichannel algebraic scattering (MCAS) method [70] was used
to obtain the nonlocal potential resulting from channel coupling. The resulting nonlocal
potential was found to be very different from the simple Gaussian nonlocalities assumed
in phenomenological potentials. However, the MCAS method is only suitable for very low
12
energy projectiles, where just a few excited states are relevant to the coupling, and thus, can
The formalism to develop a microscopic optical potential is complicated, and requires con-
for widespread use. The seminal work of Perey and Buck, [1], was the first attempt to con-
strain the parameters of a nonlocal potential through fits to elastic scattering data. This
work was done in the sixties, but it is still the most commonly referenced nonlocal optical
potential. In the late seventies, Giannini and Ricco constructed a phenomenological nonlo-
cal optical potential, [71, 72]. In that work, the potential parameters were constrained with
fits to a local form, then a transformation formula was used to obtain the nonlocal poten-
tial. Very recently, Tian, Pang, and Ma (TPM) introduced a third nonlocal global optical
potential, [2]. These three works are to our knowledge the only attempts to construct a
A common feature of using a nonlocal potential is that the amplitude of the wave function
is reduced in the nuclear interior as compared to the wave function resulting from using a
local potential. Numerous studies have been performed to investigate this effect, and to
find ways to correct for it [73, 74, 75]. These studies were focused on potentials of the
to the nonlocal potential should formally exist. Attempts have been made to find this local
equivalent potential [76, 77]. In nearly all these cases, the Perey-Buck form for the nonlocal
13
1.2.3 Solving Nonlocal Equations
While the theoretical foundation for constructing nonlocal potentials has been around for
many decades, the broad application of nonlocal potentials in the field of nuclear reactions
has never come to fruition. With a nonlocal potential, the Schrödinger equation transforms
forward way to solve the equation is through iterative methods, which dramatically increases
Since the knowledge of nonlocality dates back to the 1950s when computer power was
much more limited than today, the preferred method was to include nonlocality approxi-
mately through a correction factor [73, 74, 75]. However, several methods now exist that
improve the efficiency of the basic iteration scheme. Kim and Udagawa have presented a
rapid method using the Lanczos technique [78, 79]. A method by Rawitscher uses either
Chebyshev or Sturmian functions as a basis to expand the scattering wave function [80].
Computation time is no longer an issue. In this work, we used an iterative method out-
lined in Appendix A to solve the integro-differential equation. This is, by far, the easiest, but
definitely not the most efficient way to solve the equation. Since the increase in computation
time is minimal, pursuing a faster way was not a priority and will be pursued at a later time.
If one desired to construct their own global nonlocal potential by fitting large amounts of
14
1.3 Motivation for present work
In this work, we would like to describe single nucleon transfer reactions involving deuterons
while using nonlocal optical potentials. Ever since the early days of nuclear physics, right
up to the modern day, the distorted wave Born approximation (DWBA) has been a common
theory used to analyze data from transfer reaction experiments [82, 83]. In the DWBA, the
transfer process is assumed to occur in a one-step process, and an optical potential fitted to
deuteron elastic scattering is used to describe the deuteron scattering state. The shortcoming
of the DWBA is that the deuteron is loosely bound, so it is likely that the deuteron will
breakup as it approaches the nucleus. Not taking deuteron breakup into account explicitly
can have a significant effect on transfer cross sections. In all known implementations of the
DWBA to describe transfer cross sections, local deuteron optical potentials have been used.
These deuteron optical potentials were obtained either by fitting a single elastic scattering
angular distribution, or using a global parameterization such as that from Daehnick [84].
In order to include deuteron break up explicitly, it is necessary to include the n−p degrees
approach was introduced in the zero range approximation by Johnson and Soper [85], and
later extended to include finite range effects by Johnson and Tandy [86]. This is known as
the adiabatic distorted wave approximation (ADWA). A recent systematic study of (d, p)
reactions within the formalism of [86] shows the importance of finite range effects [87]. In
these theories the deuteron scattering state is treated as a three-body problem, composed
of n + p + A. The breakup of the deuteron is included explicitly, and the input potentials
are neutron and proton optical potentials, which are much better constrained than deuteron
optical potentials. In this sense, ADWA is a more advanced theory than the DWBA with
15
the added advantage that nucleon optical potentials exist in a nonlocal form. Therefore, in
this work the explicit inclusion of nonlocality in single nucleon transfer reactions within the
As mentioned before, nonlocality in (d, p) transfer reactions has traditionally been in-
cluded approximately through a correction factor. This is the method exploited in com-
monly used transfer reaction codes such as TWOFNR [24]. The bound and scattering wave
functions are calculated using a suitable local potential, normally a global potential for elas-
tic scattering and a mean field reproducing the experimental binding energy for the bound
state. The correction factor used implies that the nonlocality assumed is of the Perey-Buck
form. From microscopic calculations, it is known that a single Gaussian is not sufficient to
take into account the complex nature of nonlocality [54]. Therefore, not only is this method
of including nonlocality not accurate, but it is limited to a form for the nonlocality that may
Recently, some attempts have been made to include nonlocality within the adiabatic
model by introducing an energy shift to the optical potentials used to calculate the scattering
wave functions [22, 21]. This method is very attractive as all local codes which calculate
(d, p) transfer can still be used without modification. However, the adequacy of this energy
shift to take nonlocality into account must be quantified. Another limitation of this method
is that it relies on energy independent nonlocal nucleon optical potentials assumed to have
While the existence of nonlocality in the optical model has been known for many decades,
not many calculations of transfer reactions with the explicit inclusion of nonlocality have ever
been performed. While the approximate ways to correct for nonlocality are common, it is
not known if these approximate methods are sufficient. The method of constructing local
16
optical potentials through fits to elastic scattering data has been practical and useful, but
since elastic scattering does not constrain the short range nonlocalities present in the nuclear
potential, it is unlikely this approach to constructing the optical potential will be reliable
when moving towards exotic nuclei. Also, it must be understood how other observables are
affected due to the way in which the optical potentials are constructed.
The goal of this thesis is to study the explicit inclusion of nonlocality on single nucleon
transfer reactions involving deuterons. Since nonlocality has either been ignored or included
approximately in nearly all reaction calculations for over half a century, the effect of neglect-
ing nonlocaly on reaction observables must be quantified. Also, the quality of the commonly
used approximate techniques need to be assessed. For this purpose we extend the formalism
of the ADWA to include nonlocality. Finally, with renewed interest in microscopic optical
potentials, the formalism must be kept general so that nonlocal potentials of any form can
be used.
In this thesis, we will first test the concept of the correction factor using the Perey-Buck
potential. This will be done by performing DWBA calculations of (p, d) reactions on a wide
range of nuclei and energies. The correction factor will be applied to the proton scattering
state, and the neutron bound state in the entrance channel. We will then include nonlocality
explicitly in the entrance channel in order to quantify the adequacy of the correction factor
to account for nonlocality. For this part of the study, a local deuteron optical potential will
Since it is well known that deuteron breakup plays an important role in describing the
reaction dynamics, it is crucial to incorporate nonlocality into a reaction theory that ex-
plicitly includes deuteron breakup. Thus, we chose to extend the formalism of the ADWA
to include nonlocal potentials. Also, since the Perey-Buck form for the nonlocality is not
17
consistent with microscopic calculations, the formalism was kept general so that it can be
used with a nonlocal potential of any form that may result from a microscopic calculation.
Finally, through a systematic study, the effect of ignoring nonlocality in the optical
potential on transfer observables can be quantified. We will choose a range of nuclei and
energies, and perform calculations of (d, p) transfer reactions using nonlocal potentials in
both the entrance and exit channels. The resulting cross sections will be compared to cross
sections generated from local phase equivalent potentials in order to quantify the effect of
1.4 Outline
This thesis is organized in the following way. In chapter 2 we will present the necessary
theory. We will begin with a discussion of elastic scattering, and the two-body T-matrix.
We will extend the two-body T-matrix to three-bodies. Then we will introduce the adiabatic
distorted wave approximation, and finally extend this theory to include nonlocal potentials.
In chapter 3 we will discuss optical potentials. First we will introduce the concept of a global
optical potential, then turn our attention to nonlocal potentials. We will introduce Perey-
Buck type potentials, and the corresponding correction factor. We will then describe the
Giannini-Ricco potential and the DOM nonlocal potential. Last there will be a discussion
of local equivalent potentials. In Chapter 4 we will present our results beginning in Sec. 4.2
with a discussion of (p, d) reactions using the Perey-Buck potential in the entrance channel
within the DWBA. In Sec. 4.3 we compare the effects of including the DOM potential and
the Perey-Buck potential in the entrance channel of (p, d) reactions using the DWBA. Lastly,
in Sec. 4.4 we study (d, p) transfer reactions within the ADWA while including nonlocality
18
consistently. Finally, in Chapter 5 we will draw our conclusions and discuss the outlook for
future work.
Some of the work developed during this thesis, while critical, is too technical to present
in the main body. We have thus collected that information in the following appendices.
In Appendix A we discuss the method by which we solve the scattering and bound state
nonlocal equations. In Appendix B we derive the correction factor that is applied to wave
functions resulting from a local potential in order to account for the neglect of nonlocality.
In Appendix C we derive the nonlocal adiabatic potential, and in Appendix D we derive the
partial wave decomposition of the T-matrix used to calculate transfer reaction cross sections.
In Appendix E we go over some checks to ensure the accuracy of the code I developed to
we discuss a method to extract astrophysically relevant ANCs using the concept of mirror
symmetry. While Appendix F is a research project of relevance to the field that stands on
its own [25], it does not fit the theme of the thesis. Therefore, we chose to include it as a
separate appendix.
19
Chapter 2
Nucleons
Elastic scattering is the anchor of many reaction theories since elastic scattering wave func-
tions are often times inputs to these theories, and are used to calculate quantities such as
transfer cross sections. Elastic scattering is also the primary means by which we construct
the nuclear potential. Therefore, for reaction theory to make useful predictions, we must
The theoretical study of transfer reactions commonly uses the distorted-wave Born ap-
proximation (DWBA). In this theory, the transfer process is assumed to be a single step, and
the breakup of the deuteron is included implicitly through the deuteron optical potential.
The deuteron is loosely bound, and is likely to breakup during the course of the reaction.
Therefore, not including the breakup of the deuteron explicitly is known to be inaccurate
[88]. Despite breakup not being included explicitly, the DWBA theory is still commonly
used to describe transfer reactions due to its simplicity and the legacy of codes available.
Modern reaction theories that incorporate breakup begin with the three-body picture
of the process. The three bodies are the neutron and the proton making up the incident
deuteron, and the target nucleus. A practical method for including deuteron breakup was
introduced by Johnson and Tandy [86]. This method is usually referred to as the adiabatic
20
distorted-wave approximation (ADWA). The ADWA has been benchmarked with more ad-
vanced techniques [89, 90], and shown to be competitive. In [89], (d, p) angular distributions
for the ADWA and the exact Faddeev method are compared. It was found that the results
from the ADWA are within 10% of the full solution at forward angles, demonstrating that
the ADWA is a reliable and practical method for calculating angular distributions of transfer
An attractive feature of the ADWA is that it includes breakup explicitly, and also relies
on nucleon optical potentials, which are much better constrained than the deuteron optical
potentials used in the DWBA. In all known uses of the ADWA, local nucleon optical po-
tentials were used. However, recent studies have shown that the nonlocality of the nuclear
potential can have a significant impact on transfer cross sections [91, 12, 14]. Thus, it has
become necessary to extend the ADWA formalism to include nonlocal potentials [3].
To describe elastic scattering distributions, we begin by solving the partial wave decomposed
Schrödinger equation
~2 ∂2
L(L + 1)
− 2
− + UN (R) + VC (R) − E ψα (R) = 0, (2.1)
2µ ∂R R2
with UN (R) being some short-range nuclear potential, VC the Coulomb potential, µ the
reduced mass of the projectile target system, and E the projectile kinetic energy in the center
of mass frame. Here, α = {LIp Jp It } is a set of quantum numbers that define each partial
21
wave, where L is the orbital angular momentum between the projectile and the target, Ip
and It are the spin of the projectile and target respectively, and Jp is the angular momentum
resulting from coupling the orbital angular momentum with the spin of the projectile. In
the asymptotic limit where the nuclear potential goes to zero, the scattering wave function
i −
HL (ηL , kR) − Sα HL+ (ηL , kR) ,
ψα (R) = (2.2)
2
scattering matrix element (S-matrix), and H − and H + are the incoming and outgoing Hankel
functions [92], respectively. For neutrons, η = 0. The theoretical scattering amplitude for
with µpi and µti being the projections of the spin of the projectile and target, respectively,
before the scattering process, while µp and µt are the spin projections after the scattering
22
η h i
2 (θ/2)) + 2iσ (η) ,
fc (θ) = − exp −iη ln(sin 0 (2.4)
2k sin2 (θ/2)
The Sα are determined by matching a numerical solution of Eq.(2.1) to the known asymp-
totic form (2.2). This is done by constructing the R-Matrix, which is simply an inverse
logarithmic derivative.
1 HL− − Sα HL+
Rα = (2.5)
Rmatch H − 0 − Sα H + 0
L L
with the primes indicating derivatives with respect to R. The R-Matrix is evaluated at
some matching point outside the range of the nuclear interaction, denoted by Rmatch . The
0
HL− − Rmatch Rα HL−
Sα = 0. (2.6)
HL+ − Rmatch Rα HL+
Once the S-matrix for each partial wave is calculated, the theoretical differential cross
section, which is the quantity that is compared with experiment, is obtained by summing
the squared magnitude of the scattering amplitude over the final m-states, and averaging
23
dσ 1 X 2
= fµp µt ,µp µt (θ) (2.7)
dΩ Iˆpi Iˆti µp µ µp µ i i
t i ti
We would like to find the transition amplitude (T-matrix) for a (d, p) transfer reaction. Before
we get to transfer reactions, let us first consider the T-matrix for two-body scattering, such
as elastic scattering. The discussion of Sec. 2.1 formulated elastic scattering in terms of an
S-matrix. This is the way most codes solve elastic scattering. Another way of formulating
elastic scattering is in terms of the T-matrix, and leads to a natural generalization to three-
We begin with a partial wave decomposed two-body coupled channel equation [8]
~2
2
d L(L + 1)
hα|V |α0 iψα0 (R0 ).
X
− − + Vc (R) − E ψα (R) = − (2.8)
2µ dR2 R2 0 α
The T-matrix is an important quantity as it gives the amplitude of the outgoing wave after
scattering. In Eq.(2.2) we wrote the asymptotic form of the scattering wave function in
terms of the S-matrix. We can write an equivalent expression for the asymptotic form of the
ψααi (R) → δααi FLi (ηL , kR) + Tααi HL+ (ηL , kR), (2.9)
24
where Fα (R) is the regular Coulomb function, and again, H + is the out going Hankel
function. If U = 0 then T = 0. The goal is thus to find an expression for the T-matrix.
Using Green’s function techniques, the T-matrix for two-body scattering is given by [8]
2µ
Tααi = − 2 hφ(−) |V |Ψi, (2.10)
~ k
where φ is the homogeneous solution when no coupling potentials are present, µ is the
reduced mass of the two-body system, and k is the wave number. The (−) superscript
indicates that φ(−) has incoming spherical waves as the boundary condition. φ(−) is thus
the time reverse of φ. The complex conjugation implied in the bra-ket notation cancels the
Often times, we can decompose V into two parts so that V = U1 + U2 . We would like
to calculate the T-matrix for the transition when two potentials are present, and derive
the two-potential formula. We begin by writing the T-matrix substituting in the separated
expression for V
~2 k 1+2
Z
− T = φ (U1 + U2 ) ψdR. (2.11)
2µ
Using these two potentials, we can define various functions. φ is the free field solution,
χ is the solution distorted by U1 only, and ψ is the full solution. These are related to each
25
[E − T ]φ = 0
χ = φ + Ĝ0 U1 χ
ψ = φ + Ĝ0 (U1 + U2 )ψ
= χ + Ĝ1 U2 ψ, (2.12)
Ĝ0 = [E − T ]−1
~2 k 1+2
Z h i
− Tαα = χ(U1 + U2 )ψ − (Ĝ0 U1 χ)(U1 + U2 )ψ dR
2µ i
Z
= [φU1 χ + χU2 ψ]dR
(2.14)
Consider the elastic scattering of protons as an illustrative example. In this case, U1 could
be the Coulomb potential, and U2 could be the nuclear potential. The first term would be the
Coulomb scattering amplitude, fc (θ), and the second term would be the Coulomb-distorted
26
nuclear amplitude fn (θ). Thus, the nuclear scattering amplitude when Coulomb is present
is not simply the amplitude due to the short-ranged nuclear forces alone, but from the effect
of Coulomb on top of nuclear. From these scattering amplitudes we obtain the differential
elastic cross section by calculating |fc (θ) + fn (θ)|2 . This is used in our studies for computing
Using Eq.(2.14), and the implicit form for ψ in Eq.(2.12), we can, by iteration, form what is
(1+2)
Tααi = T(1) + T2(1)
2µ h i
= T(1) − 2 hχ(−) |U2 |χi + hχ(−) |U2 Ĝ1 U2 |χi + . . . . (2.15)
~ k
Truncating the series after the first term is known as the first-order distorted-wave Born
approximation (DWBA).
The DWBA is particularly useful when we are describing some kind of transition. If U1
is a central optical potential for all non-elastic channels, it cannot cause the transition since
central potentials are not able to change the quantum numbers of the scattered particle,
or change their energy. When this is the case, T(1) = 0, and we get an expression for the
27
4+ 4+ 4+
2+ 2+ 2+
0+ 0+ 0+
First-order Second-order All-orders
Figure 2.1: First, second, and all-order couplings within a set of 0+ , 2+ , and 4+ nuclear
levels, starting from the ground state.
2.1 illustrates first, second, and all order couplings between the 0+ ground state, and the 2+ ,
and 4+ excited states. The first-order DWBA can be thought of as a one step process, where
the ground state couples to either the 2+ state or the 4+ state. Similarly, the second-order
DWBA is a two-step process where, for example, the ground state can couple to the 2+
state, and then the 2+ state can either couple to the 4+ state or the 0+ . For a part of the
transfer reaction studies in this thesis, the first-order DWBA was used. From here on out,
We can generalize the above discussion to a three-body system. Consider the collection of
the three bodies n + p + A, with the coordinates appropriate for A(d, p)B given in Fig. 2.2.
The coordinates rnp and RdA refer to the configuration before the transfer occurs, and the
coordinates rnA and RpB are for immediately after the transfer. The Hamiltonian for the
28
three bodies is given by
where UpA (RpA ) is the core-core optical potential. We can equivalently express the two
kinetic energy terms as Trnp + TR = TrnA + TR . This allows us to write two different
dA pB
internal Hamiltonians for the bound states, Hd = Trnp + Vp (rnp ) and HB = TrnA + Vt (rnA ).
Thus, we can write the Hamiltonian in two ways, called the post and the prior form
rnp
n
RpB RdA rnA
p A RpA
d = n+p
B = n+A
where Ui,f are the entrance and exit channel optical potentials, respectively, and the Vi,f
29
Vi = Vt (rnA ) + UpA (RpA ) − Ui (RdA )
For (d, p) reactions it is advantageous to work in the post form. In such a case, we see
that UpA (RpA ) − UpB (RpB ) ≈ 0. This term is called the remnant term and approximately
cancels for all but light targets. This is because the optical potentials between p + A and
p + (A + 1) are not likely to be significantly different. We demonstrate that the remnant can
Just like in the case of elastic scattering, the differential cross section for an A(d, p)B
reaction is found by summing the squared magnitude of the scattering amplitude over the
final m-states, and averaging over initial states. The T-matrix is related to the scattering
amplitude by
µf vf
r
fµ M µp M (kf , ki ) = − T (k , k ), (2.20)
A d B 2π~2 vi µA Md µp MB f i
where the subscript i(f ) represents the initial (final) state, µf is the reduced mass, and v
is the velocity of the projectile. Here, µA , Md , µp , and MB are the projection of the spin
of the target in the entrance channel, the deuteron, the proton, and the target in the exit
channel, respectively.
We would like to find the T-matrix for A(d, p)B reactions. Using the post representation
for the Hamiltonian, Eq.(2.18), and the two-potential formula, Eq.(2.14), we can identify
30
U1 = Uf (RpB ) + Vt (rnA ) and U2 = Vf . Since Uf (RpB ) + Vt (rnA ) produces the elastic
a result, T2(1) is the only non-zero term. In our notation, the exact T-matrix for a given
post µp M µ M
Tµ M µ M (kf , ki ) = hΨk B |Vnp + ∆|ΨkA d i. (2.21)
A d p B f i
The remnant term ∆ = UpA − UpB is negligible for all but light targets. The ket in Eq.
(2.21) for the T-matrix is the full three-body wave function for n + p + A, while the bra
is the product of a proton distorted wave and the n + A bound state wave function. As
a first approximation, we can approximate the ket as a product of a deuteron bound state
and a deuteron distorted wave. This is the well known distorted wave Born approximation
M µ (+)
|Ψk d A i = ΞI µ (ξA )φji (rnp )χi (ki , rnp , RdA , ξp , ξn ), (2.22)
i A A
where ΞI µ (ξA ) is the spin function for the target, with spin IA and projection µA . φji (rnp )
A A
is the radial wave function for the bound state, which in this case is the deuteron, and ji is
the angular momentum resulting from coupling the spin of the fragment in the bound state
(the neutron) to the orbital angular momentum between the fragment and the core. The
(+)
distorted wave, χi , is given by
31
(+) 4π X L iσL JˆPi χLi Jpi (RdA )
χi (ki , rnp , RdA , ξp , ξn ) = i ie i (2.23)
ki
Li JP
Jˆd RdA
i
( )
n o
× ỸLi (k̂i ) ⊗ ΞIp (ξp ) ⊗ Ỹ`i (r̂np ) ⊗ ΞIn (ξn ) ⊗ ỸLi (R̂dA ) ,
ji J
d JP
i Jd Md
where ΞIp (ξp ) and ΞIn (ξn ) are the spin functions for the proton and neutron respectively,
with spin Ip = In = 21 . The spin of the deuteron is given by Jd = 1, and the spin of the
deuteron coupled to the orbital angular momentum between the deuteron and the target, Li ,
gives the total projectile angular momentum, JPi . The spherical harmonics, ỸL , are defined
with the phase convention that has a built in factor of iL . Therefore, ỸL = iL YL with YL
√
defined on p.133 of the book [93]. The hatted quantities are given by Jˆ = 2J + 1. The
" ! #
~2 ∂2 L (L + 1)
− 2
− i i2 + U dA + V1L
SO
J + VC (RdA ) − Ed χLi J (RdA ) = 0,
2µi ∂RdA RdA i Pi Pi
(2.24)
( )∗
µp MB
n o
(−)∗
hΨk | = ΞI (ξA ) ⊗ Ỹ` (r̂nA ) ⊗ ΞIn (ξn ) φjf (rnA )χf (kf , RpB )
f A f jf
JB MB
(2.25)
32
where φjf (r) is the n + A bound state radial wave function, jf is the angular momentum
of the bound state resulting from coupling the spin of the neutron to the orbital angular
momentum of the bound state, `f , while JB is the total angular momentum of the final
χL J (RpB )
(−)∗ 4π X −Lf iσLf ˆ f Pf
χf (kf , RpB , ξp ) = i e JP
kf Iˆp Lf JP
f RpB
f
n o
× ỸL (k̂f ) ⊗ ΞIp (ξp ) ⊗ ỸL (R̂pB ) (2.26)
f f JP
f Ip µp
" ! #
~2 ∂2 Lf (Lf + 1)
− 2
− 2
+ U pB + VISO
p Lf JP
+ VC (RpB ) − Ep χL J (RpB ) = 0
2µf ∂RpB RpB f f Pf
(2.27)
calculate the scattering amplitude, Eq.(2.20), and hence, the cross section, in a numerically
efficient way. We show in Appendix D that for a general `i and `f relative orbital angular
momentum in the initial and final bound states, the partial wave decomposition of the T-
matrix is given by
33
z
p
d
θ
x
B
A
Figure 2.3: The coordinates used to calculate the T-matrix for (d, p) transfer.
K,Li JP Lf JP K,Li JP Lf JP
X X X i f i f
TQM m = C AQM m (k̂f )I , (2.28)
Q f Q f
K Li JP Lf JP
i f
3I +j +J +2j
32π 3 Iˆn (−) p i d f
C = , (2.29)
Iˆp ki kf ĵi ĵf
34
K,Li JP Lf JP (−)K
i f
AQM m (k̂f ) = h`f `i (K)In In (0)KM |`f In (jf )`i In (ji )KM i
Q f K̂
X X 3L +L +` +` i(σL +σL )
× i i f f ie i f L̂ L̂ Jˆ Jˆ
i f Pi Pf
Li JP Lf JP
i f
×hIp Jd (ji )Lf Li (K)jf mf |Ip Lf (JP )Jd Li (JPi )jf mf i
f
X
× hLf Li (g)JP JPi (jf )QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f
g
X QMQ gmg
× Cgmg j m CL m L 0 YL mg (k̂f ),
f f f g i f
m g
(2.30)
K,Li JP Lf JP K,−M KM
(−)MK CL 0L ,−M
X X
I i f = K C` m̃K` m̃
f i K f f i i
MK m̃f m̃i
Z 2
RpB rnA
× φjf (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA )
f Pf i RdA
× YLi ,−M (R̂dA )Y` m̃ (r̂nA )Y`i m̃i (r̂np ) sin θdRpB drnA dθ.
K f f
(2.31)
The 9j symbol, hj1 j2 (j12 )j3 j4 (j34 )jm|j1 j3 (j13 )j2 j4 (j24 )j 0 m0 i, is given on p.334 of [93], while
j m
the Cj 3m 3j m are the Clebsh-Gordan coefficients. The coordinates used to calculate the
1 1 2 2
With this partial wave decomposition, the differential cross section is given by
35
dσ kf µi µf JˆB2 X
∗
= TQM m TQM . (2.32)
dΩ ki 4π 2 ~4 Jˆ2 Jˆ2 ĵ 2 Q f Q mf
d A f mf QMQ
Introducing Eq.(2.28) into Eq.(2.32) we obtain the form for the transfer cross section used
in this work.
We are interested in describing the reaction A(d, p)B where the final nucleus B = A + n is
a bound state. In principle, the scattering state for the deuteron can be modeled as a d + A
two-body problem. This is often done where the d + A optical potential is taken from fits
to deuteron elastic scattering. However, due to the loosely bound nature of the deuteron,
Here TR and Tr are the kinetic energy operators for the center of mass motion and the
n−p relative motion, respectively. Vnp is the neutron-proton interaction, while UpA and UnA
are the proton-target and neutron-target interactions. The wave function Ψ(r, R) describes
A variety of methods exist to solve the three-body problem. The Faddeev approach offers
an exact method to solve the three-body problem for a particular Hamiltonian [94], such
36
as the Hamiltonian given in Eq.(2.33). Faddeev methods are computationally expensive,
and so far current implementations have difficulties with handling heavy systems due to
the Coulomb potential. The Continuum Discretized Coupled Channel (CDCC) method
offers another means of solving the three-body problem [95]. However, this method too is
computationally expensive. The ADWA can provide a reliable description of transfer cross
sections while requiring minimal computation costs. Studies have benchmarked these three
methods and have shown that the ADWA can reliably reproduce transfer cross sections when
compared to the other two more advanced methods in the energy ranges relevant for this
Consider the three-body wave function describing the deuteron scattering state. A formal
Z
Ψ(r, R) = Φd (r)Xd (R) + dkΦk (r)XK (R), (2.34)
where Φd (r) is the deuteron bound state wave function, and Xd (R) is the elastic deuteron
center of mass scattering wave function. Φk (r) describes the relative motion of an n − p
pair, and the continuum components XK (R) describe the motion of the center of mass of
In the DWBA, Ψ(r, R) = Φd (r)Xd (R), so breakup is not included since the second term
in Eq.(2.34) is neglected, which contains all the breakup components. While it is known that
37
breakup is important to the dynamics of deuteron induced transfer reactions, calculating the
In formulating the ADWA, Johnson and Tandy [86] realized that to calculate transfer
cross sections, we need to know the three-body wave function only in the combination Vnp |Ψi,
as is seen in Eq.(2.21) with the remnant term neglected. Therefore, an alternative expansion
should be sought that accurately represents the three-body wave function within the range
of Vnp . The essence of the ADWA method [86] is to expand the three-body wave function
∞
X
Ψ(r, R) = Φi (r)Xi (R). (2.35)
i=0
The Weinberg states are a complete set of states within the range of the Vnp interaction,
Tr + αi Vnp (r) + d Φi (r) = 0, (2.36)
where d is the deuteron binding energy, and ach state is orthogonal by the relation
The first Weinberg component Φ0 (r) occurs when α0 = 1. Therefore, the first component
is simply the deuteron ground state wave function with a different normalization condition.
38
Each successive Weinberg component will contain an additional node. Since each Weinberg
state has the same binding energy, the asymptotic properties of each Weinberg state will be
identical. In Fig. 2.4 we show the first four Weinberg states when using a central Gaussian
which reproduces the binding energy and radius of the deuteron ground state, as in [96].
For the first four states, αi = {1, 5.2, 12.7, 23.4}. The inset shows the asymptotic properties
of each state. Since each Weinberg state has the same binding energy, they decay with the
Since we are only interested in describing the short-ranged properties of the three-body
wave function, having the wrong asymptotics is not a concern. For an effective expansion,
only a finite number of terms should be necessary for an adequate description of the wave
function with the inclusion of breakup. Keeping all the terms in the expansion of Eq.(2.35)
results in a complicated coupled channel set of equations to describe the scattering process.
To eliminate this complication, the typical procedure is to keep only the first term of the
To derive the adiabatic potential, we insert the expansion of the three-body wave func-
tion, Eq.(2.35), into the Schrödinger equation using our particular three-body Hamiltonian,
Eq.(2.33). Since we are keeping only the first term, we will write the wave function as
TR + Tr + UnA (Rn ) + UpA (Rp ) + Vnp − E Φ0 (r)XAD (R) = 0 (2.38)
Here, E is the total system energy given by E = Ed − d , where Ed is the incident deuteron
kinetic energy in the center of mass frame, and d is the deuteron binding energy. Since
39
0.8
0.6
]
-1/2
0.1
|φi(r)| [fm
0.01
]
-1/2
0.4
0.001
0.2
φi(r) [fm
0 5 10 15 20 25 30
-0.2
1st
2nd
-0.4 3rd
4th
-0.6
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
r [fm]
Figure 2.4: The first four Weinberg States when using a central Gaussian which reproduces
the binding energy and radius of the deuteron ground state. The inset shows the asymptotic
properties of each state.
(Tr + Vnp (r))Φ0 (r) = −d Φ0 (r), we can make this replacement giving us
TR + UnA (Rn ) + UpA (Rp ) − Ed Φ0 (r)XAD (R) = 0 (2.39)
We now multiply by hΦ0 |Vnp and use the orthogonality properties of the Weinberg states to
obtain
h i
Loc (R) − E X
TR + UAD d AD (R) = 0, (2.40)
40
Loc (R) = −hΦ |V (U (R ) + U (R ))|Φ i.
UAD (2.41)
0 np nA n pA p 0
It is important to note that XAD (R) is not the same as the elastic scattering wave
function Xd (R) in the plane wave basis of Eq.(2.34). Xd (R) describes elastic scattering,
and the potential used to generate Xd (R) would be a deuteron optical potential obtained by
fitting elastic scattering data. On the other hand, UAD does not describe deuteron elastic
scattering. In fact, the adiabatic potential is only of use to describe transfer reactions.
However, the input optical potentials, UnA and UpA , do describe elastic scattering, and are
obtain by fits to nucleon data. This is an advantageous feature of the ADWA as nucleon
optical potentials are much better constrained than deuteron optical potentials.
tion
We would like to consider the adiabatic potential in Eq.(2.41) when we are using nonlocal
will give an overview of the derivation. As an example, consider first the neutron nonlocal
Z
ÛnA Ψ(r, R) = UnA (Rn , R0n )Ψ(R0n , R0p )δ(R0p − Rp )dR0n dR0p
Z r r
= 8 UnA R − , 2R0 − R − Ψ(r − 2(R0 − R), R0 )dR0 . (2.42)
2 2
41
The coordinates used for calculating the neutron nonlocal potential are shown in Fig. 2.5,
where the open dashed circle represents the neutron in a different point in space to account
for nonlocality. Since we are calculating the optical potential for the neutron interacting
with the target, the proton remains stationary when integrating the neutron coordinate over
all space. Hence, the reason for the delta function in Eq.(2.42).
Figure 2.5: The coordinates used for constructing the neutron nonlocal potential. The open
dashed circle represents the neutron in a different point in space to account for nonlocality.
In Eq.(2.42), the Jacobian for the coordinate transformation is unity, and we integrated
over dr0 to eliminate the delta function. We used the vector definitions Rp,n = R ± 2r , where
Rp uses the “+” sign and Rn uses the “−” sign. A similar expression is found for the proton
nonlocal operator.
Since we are using only the first Weinberg state, we will drop the “0” subscript on the wave
functions, and write the expansion of the three-body wave function as Ψ(r, R) ≈ Φ(r)X(R).
42
Z r r
ÛN A Φ(r)X(R) = 8 UN A R ± , 2R0 − R ± Φ(r ± 2(R0 − R))X(R0 )dR0 .
2 2
(2.43)
Adding and subtracting R in the second argument of UN A and making the definition s =
Z
ÛN A Φ(r)X(R) = 8 UN A Rp,n , Rp,n + 2s Φ(r ± 2s)X(R + s)ds. (2.44)
In Eq.(2.23) we gave the deuteron distorted wave for each projection of angular momen-
tum of the deuteron and target. Now we need the deuteron wave function for relative motion
between d and A for each value and projection of total angular momentum, JT MT . This is
given by
JT M T
χLJ (R)
X p
Ψ(r, R) ≈ Φ(r)X(R) = φ` (r) (2.45)
R
`LJp
nnn o o o
× Ξ1/2 (ξn ) ⊗ Ξ1/2 (ξp ) ⊗ Ỹ` (r̂) ⊗ ỸL (R̂) ⊗ ΞIt (ξt ) .
1 1 Jp JT M T
The description of each term is given after Eq.(2.23). The coordinates for constructing the
system wave function for the deuteron scattering state are given in Fig. 2.6.
43
n
In A
R It
L
Ip
p
Figure 2.6: The coordinates used for constructing the system wave function for the d + A
deuteron scattering state.
h i
T̂R + VC (R) + Uso (R) − Ed Φ(r)X(R) = − ÛnA + ÛpA Φ(r)X(R),
(2.46)
where Uso (R) is the sum of the neutron and proton spin-orbit potentials. To begin the
( )∗
X nnn o o o
φ`0 (r)Vnp (r) Ξ1/2 (ξn ) ⊗ Ξ1/2 (ξp ) ⊗ Ỹ`0 (r̂) ⊗ ỸL0 (R̂) 0 ⊗ ΞIt (ξt )
1 1 Jp
`0 J T MT
and integrate over dr, dΩR , dξn , dξp and dξt . The lhs of the equation becomes
1 ~2
2
∂ L(L + 1) JT MT
− + E d χ LJp (R). (2.47)
R 2µ ∂R2 R2
44
As we only considered ` = 0 deuterons in our calculations, let us make this assumption right
at the beginning of our partial wave decomposition of the rhs. Therefore, the two Ỹ` (r̂)
terms give 1/4π, and the partial wave decomposition of the rhs of Eq. (2.46) is
J M
Z χ T0 0 T (|R + s|)
8 X L Jp
− φ0 (r)Vnp (r)UN A Rp,n , Rp,n + 2s φ0 (|r ± 2s|) (2.48)
4π 0
|R + s|
L0 Jp
n o ∗ n o
× Ξ1 (ξnp ) ⊗ ỸL (R̂) ⊗ ΞIt (ξt ) Ξ1 (ξnp ) ⊗ ỸL0 (R + s) 0 ⊗ ΞIt (ξt )
\
Jp JT MT Jp JT MT
× dsdrdΩR dξt dξn dξp .
Our goal is to couple the integrand up to zero angular momentum. This will be spherically
symmetric so we can use symmetry to reduce the dimensionality of the integral. After several
2 2
~ ∂ L(L + 1) JT MT
2
− 2
− VC (R) − Uso (R) + Ed χLJ (R) (2.49)
2µ ∂R R p
√ Z JT MT
8R π χLJ (|R + s|)
p
= − φ0 (r)Vnp (r) YL0 (R
\ + s)
L̂ |R + s|
× UnA (Rn , Rn + 2s) φ0 (|r − 2s|) + UpA Rp , Rp + 2s φ0 (|r + 2s|) r2 sin θr drdθr ds.
This is ultimately the nonlocal equation we solve to obtain the adiabatic wave that represents
d + A initial scattering to be introduced into the T-matrix, Eq.(2.21). More details of the
45
2.7 Spectroscopic Factors
Transfer reactions are performed not only to extract spin and parity assignments of energy
levels, but also to extract spectroscopic factors. As an example to understand the concept of
a spectroscopic factor, let us consider the 17 O nucleus, which can be modeled as a 16 O core
plus a valence neutron. Let us assume that 16 O contains only a 0+ ground state and a 2+
excited state. The ground state of 17 O is a 5/2+ state. Due to the possible excited states of
the core, the ground state of 17 O can occur in various configurations. Here we consider only
h i
|17 Og.s. i = α1 16 O(0+ ) ⊗ n
1d5/2 + α2 16 O(2+ ) ⊗ n2s (2.50)
1/2 5/2+
5/2+
These two configurations for 17 O are: the ground state of the 16 O core coupled to the valence
nucleon in a 1d5/2 orbital, and the 16 O core in its excited 2+ state coupled to the valence
neutron in a 2s1/2 orbital. Both configurations must correspond to the ground state energy,
which means that the available energy for the neutron in the 1d5/2 orbital is different than
that for the neutron in the 2s1/2 orbital due to core excitation. The spectroscopic factor
tells us how probable it is to find the valence neutron in 17 Og.s. in a 1d5/2 configuration
46
The spectroscopic factor for the configuration with 16 O in the ground state can often be
cleanly extracted from the 16 O(d, p)17 Og.s. . The reason being the fast radial fall off for the
other configurations due to the additional binding of the neutron caused by core excitation.
In a simple theoretical DWBA analysis, only the first configuration of Eq.(2.50) is included in
the calculation. The peak of the transfer distribution corresponds to impact parameters for
the deuteron grazing the surface. Therefore, one expects the transfer process to adequately
be described as a one-step process. Since we left out all of the other configurations, our
T heory
theory assumed that α1 = 1, so most often it will over-predict the transfer cross
section at the peak. By normalizing the theoretical transfer distribution at the first peak to
the experimental distribution at the first peak, we can extract the physical |α1 |2 value. It
is for this reason that we are interested in the magnitude of the transfer cross section at the
47
Chapter 3
Optical Potentials
Effective potentials describing the scattering process are needed when doing calculations of
reactions, and an accurate theoretical description of these reactions is required for the reliable
primarily from elastic scattering data, but sometimes from absorption cross sections and
polarization observables [98, 23, 49, 99, 48]. In all optical potentials, the nuclear potential is
assumed to be complex, where the imaginary part takes into account loss of flux to non-elastic
channels.
In all commonly used global optical potentials, the interaction is assumed to be local.
As a consequence, these potentials all have a strong energy dependence. Inherent in the
local assumption of the potential is a factoring out of the many-body degrees of freedom.
Therefore, the anti-symmetrization of the many-body wave function, and the coupling to
all the non-elastic channels, is not explicitly taken into account, and must be introduced
effectively into the local potential through an energy dependence of the parameters.
Global optical potentials are often used in the analysis of nuclear reactions. Global potentials
are very convenient as they can easily be extrapolated to regions of the nuclear chart where
data is not available, or they can be used at energies where data has not been taken. However,
48
such extrapolations should always be done carefully. A global optical potential constructed
from fits to stable nuclei may not give sensible results when extrapolated to exotic nuclei.
Nonetheless, using a global potential is sometimes the only option available when making
Global optical potentials attempt to describe the nuclear potential across some range
of mass and energy. To do this, some kind of form for the complex mean field must be
assumed. In most constructions of global optical potentials, the real and imaginary parts
contain combinations of Volume (v), Surface (d), and Spin-Orbit (so) terms given by
where
" !#−1
R − rA1/3
f (R, r, a) = 1 + exp . (3.2)
a
Z Z e 2 2
R
1 22 3− 2 if R < Rc
RC
VC (R) = (3.3)
2
Z1 Z2 e
if R ≥ Rc ,
R
49
where the Coulomb radius is given by Rc = rc A1/3 . Given this definition of the optical
potential there are, in principle, 19 free parameters: 3 parameters per term, 6 terms assuming
the volume, surface, and spin-orbit terms are all complex, and the Coulomb radius.
The real part usually comes from the density distribution of the nucleus, ρ(r), which is
typically of a Woods-Saxon form. This justifies the real volume term, and if ρ(r) has surface
ripples, then one would need a real surface potential as well. The imaginary volume term is
responsible for loss of flux from the elastic channel occurring somewhere inside the nucleus.
This term is sometimes used in global optical potentials, and becomes more important at
higher energies. The imaginary surface term is responsible for removing flux due to non-
elastic events occurring at the nuclear surface. This is a very important term and is included
in all global optical potentials because most reactions occur at the surface. The spin-orbit
term is the interaction between the spin of the projectile and its orbital angular momentum
with the target. A real spin-orbit interaction is always included in global optical potentials.
An imaginary spin-orbit term is sometimes included, but the depth of the imaginary part is
Normally, 19 parameters are too many to constrain from just elastic scattering, so one
needs to further constrain the form of the global optical potential. Phenomenology has
guided us towards the basic form of a real volume, an imaginary surface, and a real spin-
orbit term. Some global optical potentials, such as [49], include an imaginary volume and
imaginary spin-orbit term as well. These terms are normally needed when higher energy
Once a functional form is chosen, the free parameters are varied to obtain a best fit
to a large amount of elastic scattering data. In these fits, the depth, and sometimes the
radius and diffuseness, of the various terms can be energy and mass dependent. When
50
using the potential, the target mass, charge, and the projectile energy must be specified.
Then, the value for the depth, radius, and diffuseness of each term is calculated. There are
several global optical potentials on the market, some of the common ones are discussed here
We already discussed in Sec. 1.2 the sources of nonlocality for the effective NA interaction.
Here we provide additional perspective based on Feshbach’s work [52, 53]. When derived
from the many-body problem, the single particle Schrödinger equation describing the motion
equation for the single particle motion can be derived. The formalism of Feshbach uses the
projection operator P to project the many-body wave function onto the channels that are
considered explicitly, and the projection operator Q projects onto all channels left out from
the model space. Consider the case when P projects the many-body wave function, Ψ, onto
the elastic channel. When this is the case, the operators are defined by
where |Ψgs i gives the elastic scattering channel where the target remains in its ground state,
51
Using this projection operator formalism, Feshbach showed that a formal expression for the
1
E − TR − hΦgs |V |Φgs i − hΦgs |V Q QV |Φgs i Xel (R) = 0, (3.6)
E − QHQ
with V being the bare projectile-target interaction. From this equation, we can identify the
optical potential as
1
U =V +VQ QV. (3.7)
E − QHQ
The first term appears local while the second is inherently nonlocal. If we allow for anti-
symmetriztion between the projectile and all the nucleons of the target, whereby the incident
nucleon may not the be same as the exiting nucleon, even the first term becomes nonlocal.
In the Hartree-Fock theory, used for bound state calculations, the naturally arising exchange
While Eq.(3.7) is a formal equation, it gives some physical insight into the nature of
nonlocality. The nucleon begins in the space of elastic scattering, P-space. The system then
couples to some non-elastic channel and propagates through that space, Q-space, before
returning back to the elastic channel at some later location in space. This also gives a
physical justification for the need to have a potential with two arguments, U (R, R0 ). Flux
leaves P-space and goes into Q-space at R0 . The flux propagates through Q-space, before
52
With a nonlocal potential, the Schrödinger equation, Eq.(2.1), gets transformed into an
integro-differential equation
~2 2
Z
∇ Ψ(R) + EΨ(R) = Uo (R)Ψ(R) + U N L (R, R0 )Ψ(R0 )dR0 . (3.8)
2µ
To describe the physics of flux leaving P-space at R0 , propagating through Q-space, and
dependent on the overlap of the wave function and the potential at all other points in space,
Due to the success of local global optical potentials, it would be natural to assume that
similar global parameterizations have been made to nonlocal potentials. Unfortunately, this
is not the case. To our knowledge, there are only three global nonlocal optical potentials
that are constructed phenomenologically from elastic scattering. The seminal paper of Perey
and Buck in 1962 was the first attempt to make a parameter set for a nonlocal model [1].
In the late 70s Giannini and Ricco constructed their potential by fitting a large amount of
data to a local form. They then used an approximate transformation formula to obtain the
nonlocal parameters [71, 72]. Finally, in 2015, Tian, Pang, and Ma (TPM) constructed their
While the existence of nonlocality in the nuclear potential has long been known, there
has historically been great difficulty in specifying the exact form for the nonlocal nuclear
53
potential. A simple form was first proposed by Frahn and Lemmer [100], and later developed
and implemented by Perey and Buck [1]. The Perey-Buck potential is the most commonly
referred to phenomenological nonlocal optical potential due to its simplicity. The Perey-Buck
potential is given by
R + R0
UPNBL (R, R0 ) R − R0 , β ,
= U H (3.9)
2
R+R0
where the function U 2 is of a Woods-Saxon form, and the function H( R − R0 , β)
2
exp − R−R0
β
R − R0
H ,β = 3
. (3.10)
π 2 β3
0
Making the definition p = R+R
2 , U (p) has a form similar to those in local optical model
calculations. For the Perey-Buck nonlocal potential, U (p) consists of a nonlocal real volume,
The parameter that defines the range of the nonlocality is β. As a physical example
is possible for the incident nucleon to not be the same as the scattered nucleon. For the
incident nucleon to “switch places” with one of the nucleons within the target, it is reasonable
to assume that the two nucleons must be relatively close to each other for this to occur.
Typically, nonlocality ranges are of the order of the size of the nucleon. For the Perey-Buck
potential, β is fixed at 0.85 fm. For other nonlocal potentials, such as the TPM, β is an
54
additional parameter in their fit. The resulting value for β in neutron and proton versions
by fitting elastic scattering. Two data sets were used: n+208 Pb at 7.0 and 14.5 MeV. Perey
and Buck assumed that the parameters were energy and mass independent. Therefore, a
single parameter set completely defines the nonlocal potential of Perey and Buck. The
The work of Tian, Pang, and Ma (TPM) was the first modern attempt to find a parameter
set for a nonlocal potential [2]. In their fit, a multitude of data was considered, spanning
energy and mass. This is a great improvement over the two data sets Perey and Buck used
in their fit. A separate potential for protons and neutrons was found for the TPM potential,
unlike Perey and Buck where no protons were used in the fit. As with Perey and Buck, the
parameters in the TPM potential are assumed energy and mass independent. For higher
energy reactions, the TPM potential was found to provide a better χ2 than the Perey-Buck
potential, while at lower energies the two potentials are comparable. The parameter set for
The TPM potential was published after much of the work for this study was completed.
Hence, for this reason, and due to the popularity and widespread use of the Perey-Buck
potential, the potential we used to assess the effects of nonlocality was that of Perey and
Buck. In this study, we are interested in differences between nonlocal and local equivalent
calculations, and not so much on the quality of the nonlocal calculations themselves. When
it becomes necessary to use nonlocal potentials to extract information from experiments, the
With local optical potentials, hundreds of elastic scattering data sets using both protons
55
Perey-Buck TPM TPM
Neutrons Protons
Vv 71.00 70.00 70.95
rv 1.22 1.25 1.29
av 0.65 0.61 0.58
Wv — 1.39 9.03
rwv — 1.17 1.24
aw v — 0.55 0.50
Wd 15.00 21.11 15.74
rwd 1.22 1.15 1.20
awd 0.47 0.46 0.45
Vso 7.18 9.00 8.13
rso 1.22 1.10 1.02
aso 0.65 0.59 0.59
rc 1.22 — 1.34
β 0.85 0.90 0.88
Table 3.1: Potential parameters for the Perey-Buck [1] and TPM [2] nonlocal potentials.
and neutrons scattering off of a range of nuclei at a range of energies are used to constrain the
parameters of the potential. Therefore, it would be expected that the two data sets Perey
and Buck used to constrain their parameters would not be sufficient to reproduce elastic
scattering over a wide range of nuclei and energies. Nonetheless, reasonable agreement with
data is seen despite the simplistic way in which the potential parameters are constrained, as
In order to solve the nonlocal equation, we first need to do a partial wave expansion of
X gL (R, R0 )
UPNBL (R, R0 ) = ∗ (R̂0 )
YLM (R̂)YLM
RR0
LM
X 2L + 1 gL (R, R0 )
= PL (cos θ), (3.11)
4π RR0
L
56
4.0 2.00
2.5 1.25
σ/σR
σ/σR
2.0 1.00
1.5 0.75
a) b)
1.0 0.50
0.5 0.25
0.0 0.00
0 30 60 90 120 150 180 0 20 40 60 80 100 120
θc.m. [deg] θc.m. [deg]
where we defined θ as the angle between R and R0 . Multiplying both sides by PL (cos θ),
integrating over all angles, using the orthogonality of the Legendre polynomials, and solving
Z 1
gL (R, R0 ) = 2πRR0 UPNBL (R, R0 )PL (cos θ)d(cos θ). (3.12)
−1
Now, inserting the Perey-Buck form for the nonlocal potential, replacing 21 R + R0 with
1 0
2 (R + R ), and doing a few lines of algebra outlined in Appendix A, we arrive at
2iL z R2 + R02
1
gL (R, R0 ) = 1
jL (−iz)exp − U 0
(R + R ) (3.13)
β2 2
π β
2
0
with z= 2RR
2 , and jL being spherical Bessel functions. We now have a partial wave equation
β
~ 2 d2
Z
L(L + 1)
− χL (R) + EχL (R) = gL (R, R0 )χL (R0 )dR0 . (3.14)
2µ dr2 R2
57
3.3.1 Correction Factor
Ever since Perey and Buck introduced their potential in 1962, nearly all analytic work
nonlocal potential, assumed the Perey-Buck form for the nonlocality. That is, there is only
one nonlocality parameter, β, and the nonlocal part of the potential takes the form Eq.(3.9).
This is not true for the DOM or the Giannini-Ricco nonlocal potential, discussed later, where
Accounting for the nonlocality through the energy dependence of a local optical potential
is known to be insufficient. One key feature of a nonlocal potential is that it reduces the
amplitude of the wave function in the nuclear interior compared to the wave function from an
equivalent local potential. This is the so-called Perey effect [101]. Physically, the reduction of
the wave function can be understood to result from the repulsion due to the Pauli principle.
Since it wasn’t practical to solve the integro-differential equation with a nonlocal poten-
tial in the 1960s, there was great interest in finding a way to account for this reduction of
amplitude while still keeping the simplicity of solving a local equation. This was first accom-
plished by Austern, who studied the wave functions of nonlocal potentials and demonstrated
the Perey effect in one dimension [73]. Later, Fiedeldey did a similar study for the three
dimensional case [75]. Using a different method, Austern presented a way to relate wave
functions obtained from a nonlocal and a local potential in the three-dimensional case [74].
Since then, nonlocal calculations have been avoided using the Perey correction factor (PCF).
The Perey correction factor is derived in detail in Appendix B. Here we simply outline
the derivation. To derive the PCF, we begin with the three dimensional Schrödinger equation
58
~2 2 N L
Z
∇ Ψ (R) + EΨN L (R) = Uo (R)ΨN L (R) + U N L (R, R0 )ΨN L (R0 )dR0 , (3.15)
2µ
where Uo (R) is the local part of the potential, and typically contains spin-orbit and Coulomb
terms. Let us define a function, F (R) that connects the local wave function Ψloc (R) resulting
from the potential U LE (R) with the wave function resulting from a nonlocal potential,
ΨN L (R):
The potential U LE (R) is defined such that it reproduces the exact same elastic scattering
as the nonlocal potential. Since the local and nonlocal equations describe the same elastic
scattering, the wave functions should be identical outside the nuclear interior. Thus, F (R) →
~2 2 loc
∇ Ψ (R) + EΨloc (R) = U LE (R)Ψloc (R). (3.17)
2µ
µβ 2 h LE i−1/2
F (r) = 1 − 2 U (R) − Uo (R) . (3.18)
2~
59
It should be noted that the PCF is only valid for nonlocal potentials of the Perey-Buck
form. However, there is no reason to expect that the full nonlocality in the optical potential
will look anything like the Perey-Buck form. On physical grounds, the optical potential
must be energy dependent due to nonlocalities arising from channel couplings. While the
specific form chosen for the Perey-Buck potential is convenient for numerical calculations, a
oversimplification.
The Perey-Buck potential remained the only widely known and used nonlocal potential
available for the following 15 years after its development. An attempt by Giannini and
Ricco was made to construct a nonlocal potential of a similar form but with more data
constraining the parameters [71, 72]. Their first work focused on N = Z spherical nuclei,
while their second work made an extension to N 6= Z nuclei. Unfortunately, in doing the
fits, no nonlocal calculations were performed. Instead, the fits were done using a purely local
optical potential, and a transformation formula was used which related nonlocal and local
To construct their potential, Giannini and Ricco first derived a general expression of the
nonlocal potential in the framework of Watson multiple-scattering theory [55]. The form of
the derived nonlocal potential is a guide for the parametrization of the phenomenological
optical potential, whose parameter values are fitted to both elastic scattering and bound state
properties. To do the fit, a local form for the optical potential was chosen. The parameters
were varied to obtain a best fit of the available data, and then the transformation formulas
60
were used to get the nonlocal Giannini-Ricco potential for N = Z nuclei (GR76) [71], and
later for N 6= Z nuclei (GR80) [72]. We chose not to use this potential since the fits to data
were done using local potentials, and an unreliable transformation formula was used to get
An alternative method for obtaining the optical potential is through the self-energy, which
can be calculated microscopically using modern day structure theory. This is the method
by which the Dispersive Optical Model (DOM) is constructed. The DOM makes use of the
Kramers-Kronig dispersion relation that links the imaginary and real parts of the nucleon
self-energy [63, 65]. The optical potential is constrained by this dispersion relation. This
method was first introduced by Mahaux and Sartor [102]. The nuclear mean field is a function
of energy, where for E < 0 it is the shell-model potential that describes single-particle states,
while for E > 0 it is the optical model potential that describes scattering cross sections.
While the nuclear mean field is a continuous function of energy, its behavior as the energy
changes sign is not simple due to the coupling between elastic and inelastic channels. It is
this coupling that gives rise to the energy dependence and the imaginary component in the
optical potential. Through the dispersive relation, the scattering and bound state parts of
the nuclear mean field can be linked. The scattering parameters can be constrained by use of
fitting elastic scattering, and the bound state parameters can be constrained by comparisons
to single particle energies and (e, e0 p) observables. As seen in Fig. 3.2, the nonlocal DOM
potential reproduces experimental elastic scattering data across a wide range of energies.
Using this dispersive relation, a local version of the DOM has been developed [103]. The
61
Figure 3.2: Calculated and experimental elastic scattering angular distributions using the
nonlocal DOM potential. Data for each energy are offset for clarity with the lowest energy
at the bottom and highest at the top. Data references in [11]. Figure reprinted from [11]
with permission.
local DOM was subsequently used to in the analysis of (d, p) transfer reactions on closed
shell nuclei and shown to describe transfer angular distributions with similar adequacy as
Recently, the dispersive optical model formalism has been extended to explicitly include
nonlocality, specifically for 40 Ca [11]. As compared to the Perey-Buck potential, the nonlocal
DOM has very different ranges for the nonlocality, and a different value for the nonlocal
range in the volume and surface absorption terms. The different ranges of nonlocality for
each term in the potential makes the application of a correction factor difficult. However, as
the adequacy of the correction factor to take into account nonlocal effects has been put into
question [12], a correction factor for the DOM potential should not be sought.
62
3.6 Local Equivalent Potentials
To assess the effects of nonlocal potentials, a local phase equivalent (LPE) potential needs
reproduces the same elastic scattering. This definition is chosen since optical potentials are
constructed through fits to elastic scattering data. Therefore, if two potentials are able to
generate the same elastic scattering distribution, then the two potentials are indistinguishable
The downside of this definition is that the short-ranged nonlocal effects are not con-
strained through elastic scattering. To find a LPE potential, we first assume some form
for the LPE potential. This form is normally chosen to mimic the shape of the nonlocal
potential. As an example, the Perey-Buck nonlocal potential has real volume, real spin-
orbit, and imaginary surface terms. Therefore, the LPE potential was chosen to have the
same terms. We calculate the elastic scattering distribution generated from the nonlocal
potential, then vary the parameters of our LPE potential to obtain a best fit to the elastic
scattering distribution. This was done with the code SFRESCO [104] which performs a χ2
minimization.
Another method to obtain a LPE potential is through S-matrix inversion [105, 106].
This has the advantage over elastic scattering fits since the resulting potential will exactly
reproduce the S-matrix elements you started with [107]. However, it is important to note
that the S-matrix is not an observable, so one cannot extract an S-matrix for each partial
Exactly reproducing the S-matrix elements from the nonlocal calculation when doing the
local fit was one difficulty we encountered in this study. While the fits visually looked very
63
good, there were some very minor differences between the S-matrix elements generated with
the nonlocal potential and the fitted LPE potential. There was always particular difficulty
for surface partial waves. Nonetheless, when there were differences, the differences were
Finding a LPE potential is also an attractive way to make very sophisticated calculations
of the optical potential practical, and to assess their validity. Such has been done with
the nonlocal optical potential generated from multiple scattering [108] using the S-matrix
inversion technique. Here, the PCF was calculated by taking the ratio of the wave function
generated from the nonlocal potential with that from the local potential. A similar procedure
was done using the g-folding model for p+12 C scattering at various energies [109]. In that
study, they investigated the energy dependence of the equivalent local potential, showing
that this energy dependence does not take into account the full nonlocality, and that the
While the S-matrix inversion technique was useful to find a LPE potential in these studies,
it may not always be the most attractive way to obtain a LPE potential in practice. A fit
to elastic scattering is a much more practical and natural way to obtain a local potential,
since it is based on an observable for which one may have data. In practice, when an optical
potential is desired, there may sometimes be elastic scattering data available on the nucleus
of interest at the correct energy. When this is the case, a common procedure may be to
fit the elastic scattering data directly, rather than rely on the extrapolations of some global
potential. It is this philosophy we wanted to follow when obtaining local potentials that are
phase equivalent to a given nonlocal potential. However, rather than fitting theory to data,
An example of one such fit is shown in Fig. 3.3. We show the differential cross section
64
over the Rutherford cross section as a function of the scattering angle. The solid line is the
elastic scattering distribution generated using the Perey-Buck nonlocal potential. The open
circles are a fit to the nonlocal solution, and the dotted line is obtained by transforming the
depths of the volume and surface potentials [12]. Notice that the local fit is essentially exact
all the way out to 180◦ . The transformation formulas relied on by Giannini and Ricco to
construct their potential represents the dotted line. The inadequacy of the transformation
formula to reproduce the solution with the nonlocal potential is why the Giannini-Ricco
10
Perey-Buck
Local Fit
Transformed Depths
σ/σRuth
Figure 3.3: 49 Ca(p, p)49 Ca at 50.0 MeV: The solid line is obtained from using the Perey-Buck
nonlocal potential, the open circles are a fit to the nonlocal solution, and the dotted line is
obtained by transforming the depths of the volume and surface potentials according to Eq.
(B.14). Figure reprinted from [12] with permission.
65
Chapter 4
Results
For many years nonlocality has been effectively included in calculations by use of the Perey
correction factor (PCF) [73, 75], as discussed in Sec. 3.3.1. However, the PCF is only
suitable for use with potentials of the Perey-Buck form, and thus not of use for the DOM
potential or a microscopically derived optical potential. In addition, the quality of the PCF
has never been rigorously tested. Therefore, the first part of this study was to investigate
if the correction factor was adequately able to account for the impact of nonlocality on the
wave functions.
The reduction of the wave function can be understood physically in terms of the repulsion
between fermions due to the Pauli principle. Since one major source of nonlocality is due
to anti-symmetrization, this repulsion will naturally have the effect of pushing some of the
wave function out of the interior as compared to an interaction that doesn’t take anti-
The deuteron scattering state is also affected by nonlocality. When using the DWBA, it
is possible to apply a correction factor to the deuteron scattering wave function. However, a
nonlocal global deuteron optical potential does not exist for the purpose of comparison. Also,
as we have discussed, the DWBA does not take deuteron breakup into account explicitly.
Therefore, we would like to use the more advanced ADWA which does consider breakup, and
relies on better constrained nucleon optical potentials, of which nonlocal global potentials
66
The numerical details of the calculations performed in this thesis are presented in Sec.
4.1. The results of this thesis will be presented in three parts. The first part in Sec.
4.2 investigates (p, d) transfer reactions on 17 O, 41 Ca, 49 Ca, 127 Sn, 133 Sn, and 209 Pb at
proton energies of Ep = 20 and 50 MeV. The transfer cross sections were calculated within
the DWBA, and nonlocality in the deuteron channel is not included. In this study, we
investigated the effect of nonlocality on the proton scattering wave function and the neutron
bound state wave function. We also examined the validity of the commonly used PCF to
Ep = 20, 35, and 50 MeV using the nonlocal DOM potential, as well as the Perey-Buck
potential. Once again, the transfer cross section was calculated within the DWBA, and
nonlocality in the deuteron channel is not included. Here we studied hole states rather
than single particle states, as in the previous study. Hence, the goal of this study was to
understand if the effects of nonlocality seen in the previous study could be generalized to
hole states, and to see if the same conclusions can be drawn when using a different form for
Finally, in Sec. 4.4, we studied (d, p) reactions on 16 O, 40 Ca, 48 Ca, 126 Sn, 132 Sn, and
208 Pb at deuteron energies of Ed = 10, 20 and 50 MeV. For these cases, nonlocality was
included explicitly in the deuteron scattering state within the ADWA, as well as in the
proton channel. In all wave functions, the Perey-Buck nonlocal potential was used. This
study sought to quantify the effect of nonlocality when included consistently in calculations
It is important to note that the purpose of this work is not to describe the data. We do not
expect that the Perey-Buck potential, developed in the sixties for n+208 Pb at intermediate
67
energies using two data sets, will do well for a wide range of targets and energies. The
focus should be on the differences between the nonlocal and the local calculations under
the constraint of the same physical input, namely, that both the nonlocal and local optical
To compare the results of the nonlocal calculations, we must compare our results to calcu-
lations using local potentials with the same constraints. Therefore, to constrain the local
nucleon-target optical potentials, we require that they reproduce the same elastic scattering
obtained when using the Perey-Buck or the DOM potential at the relevant energies. For the
proton scattering states, we calculate (p, p) elastic scattering at the relevant energy using the
Perey-Buck or DOM potential, then fit the resulting distribution to a local form. The fitting
of these local phase equivalent (LPE) potentials was performed using the code SFRESCO
[104].
For the deuteron scattering states, the procedure is somewhat different. In Secs. 4.2
and 4.3 we use the local global deuteron optical potential of Daehnick [84] evaluated at the
relevant energy. In Sec. 4.4, we calculate (n, n) and (p, p) elastic scattering at half the
deuteron energy using the Perey-Buck potential, and again found LPE potentials for the
elastic scattering distributions. The local adiabatic potential is then calculated with the
For the neutron bound states, we calculated the nonlocal equation using a real Woods-
Saxon form with a nonlocality range of β = 0.85 fm. We also used a local spin-orbit
interaction with a depth fixed at 6 MeV. For each term we used a radius of r = 1.25 fm
68
and a = 0.65 fm. The depth of the nonlocal real Woods-Saxon form was then adjusted to
reproduce the physical binding energy. Also, in Sec. 4.3 we used the nonlocal DOM potential
to calculate the neutron bound state. The corresponding bound state resulting from local
potentials was obtained by setting β = 0 and adjusting the local real Woods-Saxon depth to
In Secs. 4.2 and 4.3, the calculated wave functions were read into the code FRESCO to
calculate the (p, d) transfer cross sections. We used the Reid soft core interaction [110] in
the (p, d) T-matrix and to calculate the deuteron bound state. In Sec. 4.4, the bound and
scattering states that are calculated are inserted into the (d, p) T-matrix Eq.(2.21). This
was implemented in the code NLAT (NonLocal Adiabatic Transfer). The NN interaction
in this case was a central Gaussian which reproduces the binding energy and radius of the
In Secs. 4.2 and 4.3 the scattering wave functions were solved by using a 0.05 fm radial
step size with a matching radius of 40 fm. For the bound states solutions, we used a
radial step size of 0.02 fm. The matching radius was half the radius of the nucleus under
consideration, and the maximum radius was 30 fm. The cross sections contain contributions
In Sec. 4.4 the scattering wave functions were calculated in steps of 0.01 fm with a
matching radius of 30 fm. The nonlocal adiabatic potential was obtained on a radial grid of
step 0.05 fm. We used linear interpolation to calculate the nonlocal adiabatic potential in
steps of 0.01 fm in order to calculate the adiabatic deuteron wave function with the same
step size. The bound state wave functions were also calculated in steps of 0.01 fm with a
maximum radius of 30 fm and a matching radius of half the radius of the nucleus under
69
16 O(d, p) 40 Ca(d, p) 48 Ca(d, p) 126 Sn(d, p) 132 Sn(d, p) 208 Pb(d, p)
Table 4.1: Percent difference of the (d, p) transfer cross section at the first peak for a calcu-
lation including the remnant term relative to a calculation without the remnant term.
To get an idea of the significance of the remnant term, we show in Table 4.1 the percent
difference at the first peak of the (d, p) transfer cross section for a calculation with the
remnant term relative to a calculation without the remnant term for a wide range of targets.
These DWBA calculations used the deuteron global optical potential of Daehnick [84] to
describe the deuteron scattering state, and the LPE potentials to the Perey-Buck potential
for the proton scattering state, a central Gaussian for the deuteron bound state, and a real
Woods-Saxon form for the neutron bound state that reproduces the experimental binding
energy.
Buck Potential
The first part of this study was to investigate the effect of the Perey-Buck potential on the
entrance channel of (p, d) transfer reactions. For this study, nonlocality was included explic-
itly in the proton scattering state and the neutron bound state. Using the wave functions
generated with the nonlocal potentials, (p, d) transfer reactions were calculated. These cross
sections were compared to those generated with LPE potentials, discussed in Sec. 3.6. Also,
70
wave functions were modified with the PCF, and the corresponding transfer cross sections
were calculated. The goal in this study was to assess the effect of nonlocality on transfer
cross sections when compared to cross sections generated with LPE potentials, as well as to
determine the quality of the PCF and its ability to reproduce the effects of nonlocality.
When doing calculations of (p, d) or (d, p) reactions using the T-matrix formalism of Eq.(2.21),
it is required to calculate a proton elastic scattering wave function in either the entrance
or exit channel. Because some codes, such as TWOFNR [24], allow for nonlocality to be
included through the PCF, this approach has become common practice. However, until re-
cently, the accuracy of this approach was not understood. Using the Perey-Buck potential,
this methodology to include nonlocality has been tested. To do this check, a LPE potential
needed to be found. The open circles in Fig. 3.3 are one such example of a LPE potential.
The LPE potential found from the fit is the U LE term in Eq.(3.18).
As an example, we will use the LPE potential from Fig. 3.3 and consider the scattering
wave function for the reaction 49 Ca(p, p)49 Ca at 50 MeV. The J π = 1/2− partial wave is
shown in Fig. 4.1. As is seen in the figure, the reduction of the wave function resulting
from the nonlocal potential (solid line) relative to the wave function from the LPE potential
(dashed line) is apparent. Also seen is that the wave functions from the nonlocal potential
and the local potential with the PCF applied (crosses) are in good agreement. This was
a general result for most partial waves. However, in all cases that were studied, problems
arose for partial waves corresponding to impact parameters around the surface region, shown
in Fig. 4.2. Since transfer cross sections tend to be most sensitive to the surface region,
the differences for these angular momenta are particularly relevant. We will see how the
71
0.8
0.6 Nonlocal
Corrected Local
|Re(ψ)|
Uncorrected Local
]
0.4
-1/2
ψ(r) [fm
0.2
-|Im(ψ)|
-0.2
-0.4
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
r [fm]
Figure 4.1: Real and imaginary parts of the J π = 1/2− partial wave of the scattering wave
function for the reaction 49 Ca(p, p)49 Ca at 50.0 MeV: ψ N L (solid line), ψ P CF (crosses), and
ψ loc (dashed line). Top (bottom) panel: absolute value of the real (imaginary) part of the
scattering wave function. Figure reprinted from [12] with permission.
inadequacy of the PCF for surface partial wave affects the resulting transfer cross sections
in Sec. 4.2.3.
0.8
0.6
|Re(ψ)|
Nonlocal
0.4 Corrected Local
]
-1/2
Uncorrected Local
0.2
ψ(r) [fm
0
-|Im(ψ)|
-0.2
-0.4
-0.6
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
r [fm]
Figure 4.2: Real and imaginary parts of the J π = 11/2+ partial wave of the scattering
wave function for the reaction 49 Ca(p, p)49 Ca at 50.0 MeV. See caption of Fig. 4.1. Figure
reprinted from [12] with permission.
72
This inability of the PCF to correct surface partial waves is partly due to the way in which
it was derived. When deriving the PCF, terms related to ∇2 F were neglected, such as the
one in Eq.(B.21). This term only contributes around the nuclear surface. In addition, when
performing the local fit, we occasionally found slight differences in the S-matrix elements for
a particular partial wave. Since the scattering wave functions are normalized according to
Eq.(2.2), these small changes in the S-matrix will result in different amplitudes for the real
and imaginary parts of the scattering wave function in the asymptotic region.
We now turn our attention to the neutron bound state that exists in the entrance channel
of (p, d) reactions. In order to investigate the effects of nonlocality on the bound state wave
functions, and the adequacy of the PCF to correct for nonlocality, the PCF was applied to
the local bound state wave function, and the resulting wave function was renormalized to
unity.
To illustrate, the 2p3/2 ground state wave function for n+48 Ca is shown in Fig. 4.3.
Visually, the correction factor does an excellent job correcting for nonlocality in the bound
state. However, it is important to notice that in the surface region (2 − 5 fm), the PCF
does very little to bring the wave function resulting from the local equivalent potential into
agreement with the wave function resulting from the nonlocal potential. The inset, which
As stated in Sec. 4.2.1, the reason for the inadequacy of the PCF in the surface region
goes back to the way in which the PCF was derived. In this case, the bound wave function
has a large slope around the surface resulting in large differences between the wave function
generated using nonlocal interactions and the one generated from local interactions.
73
0.06
0.4
]
-1/2
0.04
[fm
PCF
0.02
0.2
-φ
NL
0.00
φ
]
-1/2
-0.02
0 0 1 2 3 4 5 6 7 8 9 10
φ(r) [fm r [fm]
−0.2
−0.4
Nonlocal
Corrected Local
Uncorrected Local
−0.6
0 1 2 3 4 5 6 7 8 9 10
r [fm]
Figure 4.3: Ground state, 2p3/2 , bound wave function for n+48 Ca. φN L (solid line), φP CF
(crosses), and φloc (dashed line). The inset shows the difference φN L − φP CF . Figure
reprinted from from [12] with permission.
Another important point to note is that nonlocality has the effect of increasing the
normalization of the asymptotic properties of the wave function (the ANC). Since nonlocality
reduces the amplitude of the wave function in the nuclear interior, and the wave function is
always normalized to unity, the ANC must increase. Therefore, the ANC of the bound wave
function resulting from nonlocal potentials was found to always be larger than the ANC from
local potentials, and the ANC of the corrected wave function was somewhere in between.
proximation
Now that we have studied the effect of nonlocality on the scattering and bound state wave
functions, we can investigate the effect nonlocality has on (p, d) transfer reactions when
74
nonlocality is included explicitly in the entrance channel.
As a first example, consider the transfer reaction corresponding to the wave functions we
have been studying in Secs. 4.2.1 and 4.2.2, 49 Ca(p, d)48 Ca at a proton energy of Ep = 50
MeV in the laboratory frame. The separate and combined effects of nonlocality in the bound
and scattering states are shown in Fig. 4.4. The solid line corresponds to when nonlocality
is included in both the proton scattering state and the neutron bound state, the dashed
line corresponds to the distribution obtained when only local equivalent potentials are used,
the crosses correspond to the cross section obtained when the proton scattering state and
neutron bound state wave functions are both corrected with the PCF. Also shown with the
dotted line is the cross section when nonlocality was only added to the scattering state, and
the dot-dashed line when nonlocality is only added to the neutron bound state.
5.0
Nonlocal
4.5 Corrected Local
4.0 Uncorrected Local
Nonlocal Scattering State
dσ/dΩ [mb/sr]
3.0
2.5
2.0
1.5
1.0
0.5
0.0
0 5 10 15 20 25 30 35 40 45 50
θc.m. (deg)
Figure 4.4: Angular distributions for 49 Ca(p, d)48 Ca at 50 MeV: Inclusion of nonlocality in
both the proton scattering state and the neutron bound state (solid), using LPE potentials,
then applying the correction factor to both the scattering and bound states (crosses), using
the LPE potentials without applying any corrections (dashed line), including nonlocality only
in the proton scattering state (dotted line) and including nonlocality only in the neutron
bound state (dot-dashed line). Figure reprinted from [12] with permission.
75
The results of Fig. 4.4 are unique for the cases we considered in that the shape of the
distribution was significantly changed. The reason for the significant changes around zero
degrees can be seen from an analysis of the scattering and bound wave functions. The bound
wave function has a node which occurs at a radius corresponding to the surface region for
49 Ca. Since the bound wave function has a large slope in this region, the percent difference
between the nonlocal and local wave functions can be quite large. For this case, the nonlocal
bound wave function is smaller than the local wave functions in this region, reducing the
cross section at the peak around 20 degrees. On the other hand, the magnitude of the bound
wave function is large in the asymptotic region, which increases the cross section at zero
degrees.
For the scattering wave function, the most significant differences were for partial waves
corresponding to the surface region. Also, the asymptotics of the scattering wave functions
were different due to small differences in the S-matrix, again mostly for surface partial waves.
There is an interplay between the real and imaginary parts of the scattering wave function
which influences the cross section at forward angles. Then, the complex combination of all
these effects produces the interesting behavior of the transfer cross section at forward angles.
Now consider the same reaction but at a lower energy. In Fig. 4.5 we show 49 Ca(p, d)48 Ca
at a proton energy of Ep = 20 MeV in the laboratory frame. This case is more representative
of the general features we saw in this systematic study. The nonlocality in the scattering state
had the effect of reducing the transfer cross section due to the reduction of the scattering
wave function, while the nonlocality in the bound state had the effect of enhancing the cross
section due to the increase of the wave function in the asymptotic region. At this lower
energy, the overall effect was an enhancement of the transfer cross section at the first peak.
In addition, it is seen that the PCF moves the transfer distribution generated with local
76
55
50 Nonlocal
Corrected Local
45 Uncorrected Local
40 Nonlocal Scattering State
dσ/dΩ [mb/sr]
Nonlocal Bound State
35
30
25
20
15
10
5
0
0 5 10 15 20 25 30 35 40 45 50
θc.m. (deg)
Figure 4.5: Same as in Fig. 4.4 but for 49 Ca(p, d)48 Ca at Ep = 20 MeV. Figure reprinted
from [12] with permission.
14 Nonlocal
Corrected Local
Uncorrected Local
12 Nonlocal Scattering State
dσ/dΩ [mb/sr]
Figure 4.6: Same as in Fig. 4.4 but for 133 Sn(p, d)132 Sn at Ep = 20 MeV. Figure reprinted
from [12] with permission.
potentials in the direction of that generated with nonlocal potentials. However, the PCF
was not able to fully take the effects of nonlocality into account.
Next we consider heavier targets, such as 133 Sn and 209 Pb, at Ep = 20 MeV. In both of
these cases the inclusion of nonlocality in the scattering state decreased the cross section by
77
Nonlocal
2.0 Corrected Local
Uncorrected Local
Nonlocal Scattering State
1.6 Nonlocal Bound State
0.8
0.4
0.0
0 10 20 30 40 50 60 70 80 90
θc.m. (deg)
Figure 4.7: Same as in Fig. 4.4 but for 209 Pb(p, d)208 Pb at Ep = 20 MeV. Figure reprinted
from [12] with permission.
a smaller amount than was seen before. This is due to the low energy of the proton and the
high charge of the target. Nonlocality in the proton scattering state reduced the magnitude
of the wave function in the nuclear interior, but the energy of the proton was not high enough
to penetrate deeply due to the large Coulomb barrier. On the other hand, nonlocality in the
bound state is very significant. Since the projectile energy was low, and the charge of the
target was high, the reaction is dominated by the asymptotitcs of the bound wave function,
For 133 Sn(p, d)132 Sn at Ep = 20 MeV, Fig. 4.6, the PCF does a reasonable job tak-
ing nonlocality into account, but there are still discrepancies between the full nonlocal and
corrected local solutions. For 209 Pb(p, d)208 Pb at Ep = 20 MeV, Fig. 4.7, there are discrep-
ancies at forward angles, but the distributions resulting from nonlocal and local potentials
coincidentally agree at the peak. This agreement is accidental, and comes from the nonlocal
78
Corrected Nonlocal
Elab = 20 MeV Relative to Local Relative to Local
17 O(1d
5/2 )(p, d) 7.1% 18.8%
17 O(2s
1/2 )(p, d) 20.1% 26.5%
41 Ca(p, d) 11.4% 21.9%
49 Ca(p, d) 10.4% 17.3%
127 Sn(p, d) 17.5% 17.3%
133 Sn(p, d) 18.2% 24.4%
209 Pb(p, d) 19.4% 20.8%
Table 4.2: Percent difference of the (p, d) transfer cross sections at the first peak when using
the PCF (2nd column), or a nonlocal potential (3rd column), relative to the local calculation
with the LPE potential, for a number of reactions occurring at 20 MeV.
Corrected Nonlocal
Elab = 50 MeV Relative to Local Relative to Local
17 O(1d
5/2 )(p, d) 17.0% 35.4%
17 O(2s
1/2 )(p, d) 0.2% 12.7%
41 Ca(p, d) 2.9% 5.8%
49 Ca(p, d) −16.0% −17.1%
127 Sn(p, d) 10.1% 4.5%
133 Sn(p, d) −6.7% −16.9%
209 Pb(p, d) 8.6% 8.6%
Table 4.3: Percent difference of the (p, d) transfer cross sections at the first peak when using
the PCF (2nd column), or a nonlocal potential (3rd column), relative to the local calculation
with the LPE potential, for a number of reactions occurring at 50 MeV.
The percent difference at the first peak of the transfer distributions for all the cases that
were studied are summarized in Tables 4.7 and 4.8 for the (p, d) reactions at 20 and 50 MeV.
For the lower energy cases, it is seen that the inclusion of nonlocality provided a general
enhancement to the transfer cross section at the first peak. This is due to the enhancement
of the bound state wave function in the asymptotic region playing a more significant role in
the magnitude of the transfer cross section at low energies. At higher energies, there is a
competition between the enhancement due to the bound state, and the reduction due to the
scattering state. In most cases, there is still an enhancement of the cross section, but the
overall effect is less significant than for the lower energy cases.
79
4.2.4 Summary
In this study, the long established Perey correction factor (PCF) and the effects of nonlocal-
ity on the entrance channel of (p, d) reactions were studied. The integro-differential equation
containing the Perey-Buck nonlocal potential was solved numerically for single channel scat-
tering and bound states. A local phase equivalent (LPE) potential was obtained by fitting
the elastic distribution generated by the Perey-Buck potential. The PCF was applied to the
wave functions generated with the LPE potentials or the local equivalent binding potentials,
and the scattering and bound state wave functions were then used in a finite-range DWBA
We found that the explicit inclusion of nonlocality in the entrance channel increased the
transfer distribution at the first peak by 15 − 35%. In most cases, the transfer distribution
from using a nonlocal potential increased relative to the distribution from the local potential.
In all cases, the PCF moved the transfer distribution in the direction of the distribution which
included nonlocality explicitly. However, nonlocality was never fully taken into account with
the PCF. The PCF can be improved by including the surface terms that were neglected, and
not assuming that the local momentum approximation is valid. Such additional corrections
were not pursued since the full nonlocal solution can be calculated.
The results of the previous work by Titus and Nunes [12], covered in Sec. 4.2, demonstrated
that nonlocality is significant in the study of transfer reactions, and that the PCF is not
80
able to fully reproduce the complex effects of nonlocality. However, the previous study only
considered the Perey-Buck nonlocal potential. Recently, the nonlocal DOM potential for
40 Ca was developed [11]. We wanted to see if nonlocality remains an important ingredient to
transfer reactions when a different nonlocal potential is used. Since the DOM was constructed
only for 40 Ca, we were only able to consider 40 Ca(p, d)39 Ca reactions. For this study, we
10
1
σ/σRuth
DOM
DOM LPE
0.1 PB
PB LPE
Data (19.57 MeV)
0.01
0 30 60 90 120 150 180
θc.m. [deg]
Figure 4.8: Angular distributions for elastic scattering normalized to Rutherford for protons
on 40 Ca at Ep = 20 MeV. The elastic scattering with the DOM potential (solid line), the
DOM LPE potential (open circles), the Perey-Buck interaction (dashed line), and the Perey-
Buck LPE potential (open squares). The data (closed diamonds) from [13]. Figure reprinted
from [14] with permission.
In investigating the effect of nonlocality when using the DOM potential, we no longer con-
sidered the PCF since there is no easy generalization of the PCF to the DOM potential.
This is because each term of the DOM potential has a different value for the nonlocality
parameter, β. While it would be possible to construct a PCF for the DOM potential, the
81
10
1
σ/σRuth
0.1 DOM
DOM LPE
PB
PB LPE
Data (51.93 MeV)
0.01
0 30 60 90 120 150 180
θc.m. [deg]
Figure 4.9: Same as in Fig. 4.8 but for Ep = 50 MeV. Data from [15]. Figure reprinted from
[14] with permission.
results of the previous work by Titus and Nunes [12] made pursuing a PCF for the DOM
potential irrelevant.
In Figs. 4.8 and 4.9 we show the elastic distributions generated from the DOM and the
Perey-Buck nonlocal potentials along with the corresponding LPE potentials. We also show
the corresponding elastic scattering data at the closest relevant energy. It is seen from the
distributions that the DOM potential is superior when it comes to describing the data. This
should not be a surprise as the DOM potential was constructed from fits to nucleon elastic
scattering data on 40 Ca while the Perey-Buck potential was constructed from neutron elastic
scattering on 208 Pb at low energy. Nonetheless, the Perey-Buck potential does a reasonably
good job at describing the elastic scattering data for the energy and angular range that the
data is available.
To investigate the scattering wave functions, we consider the J π = 1/2+ partial wave
for scattering at Ep = 50 MeV in Fig. 4.10. For both the DOM and the Perey-Buck
82
0.8
0.6
]
-1/2 0.4 DOM
ψ(R) [fm DOM LPE
0.2 PB
PB LPE
0
-0.2
-0.4
0 2 4 6 8 10 12 14
R [fm]
Figure 4.10: The real and imaginary parts of the J π = 1/2+ partial wave of the scattering
wave function for the reaction 40 Ca(p, p)40 Ca at Ep = 50 MeV. This shows the wave function
resulting from the DOM potential (solid line) and its LPE potential (dotted line), the Perey-
Buck potential (dashed line) and its LPE potential (dot-dashed line). The top (bottom)
panel shows the absolute value of the real (imaginary) part of the scattering wave function.
Figure reprinted from [14] with permission.
nonlocal potential, we see the reduction of the scattering wave function relative to the wave
function from the LPE potential, which is consistent with earlier studies [12, 73, 75]. Since
the two nonlocal potentials describe different elastic scattering distributions, they will have
different S-matrix elements for each partial wave, and hence, the different normalizations in
The neutron 1d3/2 bound state wave functions using the various potentials are shown in Fig.
4.11. The DOM bound wave function was found using the potential defined in [11]. The
same general features of the bound wave functions in [12] are seen here. Nonlocality reduces
the amplitude of the bound wave function and thus pushes the wave function outward.
83
0.6
]
-1/2
0.1
φ(r) [fm
0.5
0.01
]
-1/2
φ(r) [fm 0.4
0 4 8 12
r [fm]
0.3
DOM
0.2 PB
Local
0.1
0
0 1 2 3 4 5 6 7 8 9 10
r [fm]
Figure 4.11: The neutron ground state 1d3/2 bound wave function for n+39 Ca. Shown is
the wave function obtained using the DOM potential (solid line), the Perey-Buck potential
(crosses) and the local interaction (dashed line). The inset shows the asymptotic properties
of each wave function. Figure reprinted from [14] with permission.
proximation
In Fig. 4.12 we show the transfer distributions for the three energies calculated. Shown is
the transfer distribution resulting from the DOM nonlocal potential and its LPE potential,
as well as the Perey-Buck nonlocal potential and its LPE potential. In general we see that
nonlocality for both potentials provides an enhancement of the cross section at the first
peak. This is consistent with the conclusions of [12]. However, at higher energies, there is
not as much cancellation between the scattering and bound states so that the full nonlocal
calculation still resulted in a fairly significant increase in the cross section. The key difference
is that the neutron was bound by 15.6 MeV in this study whereas the neutron was always
bound less than 10 MeV in all cases in the study of Sec. 4.2.
84
Ep (MeV) bound state scattering state full nonlocal
20 27 % -14 % 15 %
35 31 % 10 % 52 %
50 31 % -3 % 29 %
Table 4.4: Percent differences of the (p, d) transfer cross sections at the first peak at the
listed beam energies using the DOM potential relative to the calculations with the phase-
equivalent potential. Results are listed separately for the effects of nonlocality on the bound
state, the scattering state, and the total.
Table 4.5: Same as Table 4.4, but now for the Perey-Buck potential.
Shown in Table 4.4 and 4.5 are the percent difference of the (p, d) transfer cross sections
at the first peak when using the DOM and Perey-Buck potentials, respectively. In the tables
we show the separate effects of the neutron bound state and the proton scattering state, as
In most cases nonlocality in the scattering state had the effect of reducing the transfer
cross section. One exception was for Ep = 35 MeV when using the DOM potential. The
increase was due to the shape of the scattering wave function near the surface region. In this
particular case, obtaining and exact fit to the nonlocal distribution was much more difficult
than in the other cases. All other cases reduced the cross section by a similar amount. Since
the Coulomb barrier is not large for 40 Ca, there was no suppression of the nonlocal effects
in the scattering state as we seen in the previous study for heavier systems, such as 133 Sn
or 209 Pb.
The effect of nonlocality in the bound state at all proton energies was to increase the
cross section. This is because nonlocality shifts the bound wave function towards the surface
region where these transfer reactions are more sensitive. We also note that the nonlocal
85
0.7
0.6
0.5
0.4
0.3
0.2
0.1 a)
0.0
3 DOM
dσ/dΩ [mb/sr]
DOM PE
PB
PB PE
2
1
b)
1 c)
0
0 20 40 60 80
θc.m [deg]
Figure 4.12: Angular distributions for the 40 Ca(p, d)39 Ca reaction at (a) Ep = 20 MeV, (b)
Ep = 35 MeV, and (c) Ep = 50 MeV. In this figure is the transfer distribution resulting
from using the nonlocal DOM (solid line) and its LPE potential (dotted line), the Perey-
Buck potential (dashed line) and the Perey-Buck LPE potential (dot-dashed line). Figure
reprinted from [14] with permission.
86
effects for the Perey-Buck interaction are generally larger than for the previous study [12].
This is because in the previous study we were studying single particle states in closed shell
4.3.4 Summary
In this work we studied the effects of adding non-locality in the entrance channel of transfer
reactions using a nonlocal potential obtained from the dispersive optical model (DOM) and
comparing it to the results from the older Perey-Buck interaction. Our studies focus on the
40 Ca(p, d)39 Ca reaction at Ep = 20, 35 and 50 MeV. We consider the nonlocality in the
proton channel, and solve the integral-differential equation to obtain the proton scattering
and neutron bound state solutions for both nonlocal potentials. We then computed the
transfer matrix element in the DWBA, ignoring nonlocality in the deuteron channel.
Our results show that, irrespective of the details of the potential, nonlocality reduces
the strength of the wave function in the nuclear interior, an effect most noticeable in the
bound states, but also significant in scattering states. Due to the normalization condition,
nonlocality in the bound state also shifts the wave function to the periphery region, causing
an increase in the transfer cross sections. Typically, nonlocality in the scattering state acts
in the opposite direction, reducing the overall effect. When nonlocality is included in both
the bound and scattering states, the transfer cross sections are increased by ≈ 15 − 50% for
the DOM potential, in contrast with ≈ 30 − 50% obtained with the Perey-Buck interaction.
87
4.4 Nonlocal Adiabatic Distorted Wave Approxima-
The previous two studies (Titus and Nunes [12] discussed in Sec. 4.2, and Ross, Titus,
Nunes, et. al. [14] discussed in Sec. 4.3) have demonstrated that the explicit inclusion of
nonlocality, at least in the entrance channel of (p, d) reactions, is very important to take
into account explicitly. They have also shown that commonly used correction factors are not
In both of these studies, a local deuteron optical potential was used to describe the
deuteron scattering state. We will now turn our attention to studying transfer reactions
within the ADWA, discussed in Sec. 2.5, which includes deuteron breakup explicitly. As
all nucleon-target interactions. For this study, we will focus on (d, p) transfer reactions.
It should be noted that due to time reversal invariance, the cross sections for (d, p) and
(p, d) transfer reactions differ only by a statistical constant, assuming that the initial and
final states are the same. The statistical constant can be determined by detailed balance
[8]. Therefore, even though we are considering a different reaction, we are building on the
In order to compare the effect of nonlocality on the adiabatic potential, we define the rhs
N L (R). It is difficult to compare the nonlocal adiabatic potential
of Eq. (2.49) to be SAD
88
N L (R) has the scattering wave function
directly to the local adiabatic potential since SAD
N L (R) with the local corresponding quantity,
built into it. However, we can compare SAD
Loc (R) = U Loc (R)X Loc (R). After a partial wave decomposition, the source terms become
SAD AD AD
functions dependent only on the scalar R for each LJ combination of the deuteron scattering
state.
40 Nonlocal
Local
|S(R)| [MeV/fm ]
1/2
30
20
10
a) b)
0
0 1 2 3 4 5 6 7 8 9 10 11 12 1 2 3 4 5 6 7 8 9 10 11 12
R [fm]
Figure 4.13: Absolutes value of the d + A source term when nonlocal and local potentials are
used. (a) d+48 Ca at Ed = 50 MeV. (b) d+208 Pb at Ed = 50 MeV. Both are for the L = 1
and J = 0 partial wave. Figure reprinted from [3] with permission.
To get an idea of the effect of nonlocality on the adiabatic potential, we make a comparison
in Fig. 4.13 of this radial source term for the angular momentum values of L = 1 and J = 0
of the d+48 Ca and d+208 Pb wave function, both at a beam energy of Ed = 50 MeV. In
Fig. 4.14 we make the same comparison but for the L = 6 and J = 5 partial wave. In
both figures, the solid line corresponds to the nonlocal source term, while the dashed line
is its local equivalent. The magnitude of the nonlocal source term is reduced compared to
the local source term. It is also seen that the source term in the nonlocal case gets shifted
outward relative to the local case. Both these effects imprint themselves on the adiabatic
89
40
Nonlocal
Local
|S(R)| [MeV/fm ]
30
1/2
20
a) b)
10
0
0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 10
R [fm]
Figure 4.14: Absolute value of the d + A source term when nonlocal and local potentials are
used. (a) d+48 Ca at Ed = 50 MeV. (b) d+208 Pb at Ed = 50 MeV. Both are for the L = 6
and J = 5 partial wave. Figure reprinted from [3] with permission.
The necessary formalism for the local implementation of the ADWA and the nonlocal exten-
sion of the ADWA has been addressed in Chapter 2. The radial equation that must be solved
for each partial wave is given by Eq.(2.49). The rhs of this equation acts as a source term,
and the differences between the nonlocal source term and the corresponding local source
Turning our focus to the deuteron scattering wave function, Figs. 4.15 and 4.16 show the
absolute values of the d + A scattering wave function when using the ADWA with nonlocal
and local potentials. The solid line corresponds to the scattering wave function resulting from
using the nonlocal Perey-Buck potential in Eq.(2.49), while the dashed line is the scattering
wave function that results from using the local adiabatic potential, Eq.(2.41), where the
necessary are used for the nucleon optical potentials. Panel (a) is for d+48 Ca while panel
90
0.7
0.6
0.5
]
-1/2
|χ(R)| [fm
0.4
0.3
0.2
Nonlocal
Local
0.1
a) b)
0
0 2 4 6 8 10 12 14 0 2 4 6 8 10 12 14 16
R [fm]
Figure 4.15: Absolute value of the d + A scattering wave function using the ADWA theory
when nonlocal and local potentials are used. (a) d+48 Ca and (b) d+208 Pb. Both for the
L = 1 and J = 0 partial wave at Ed = 50 MeV in the laboratory frame. Figure reprinted
from [3] with permission.
Here, it is important to note that the individual n+A and p+A local optical potentials are
phase equivalent to the nonlocal Perey-Buck, but the nonlocal and local adiabatic potentials
are not phase equivalent. The adiabatic potential is only useful for calculating the deuteron
scattering wave function within the range of the Vnp interaction, and is not applicable for
calculating deuteron elastic scattering. It is for this reason that we chose for the input optical
When compared to the source term that drives this wave function in Fig. 4.13, we see that
both the wave function and the source term are reduced relative to the local counterpart. This
is the same feature that we saw when studying the proton scattering state in (p, d) reactions.
The reduction of the wave function in the interior is a common feature of using nonlocal
potentials, and can be understood physically in terms of the Pauli exclusion principle.
When studying proton scattering states, nonlocality only had the effect of reducing the
amplitude of the scattering wave function. However, differing from the proton scattering
91
state, the d + A scattering wave function is also shifted outward relative to the wave function
resulting from local potentials (see Fig. 4.15). This is analogous to the bound state case
where the wave function was both reduced and shifted outward due to nonlocality. This
shifting outward of the d + A scattering wave function changes the amplitude of the wave
function at the nuclear surface. Since the surface region is where (d, p) cross sections are
most sensitive, the shifting outward can have a significant effect on the cross section. In
fact, as we will see shortly, nonlocality in the deuteron scattering state increases the transfer
0.7
0.6
0.5
]
-1/2
|χ(R)| [fm
0.4
0.3
0.2
0.1
a) b)
0
0 2 4 6 8 10 12 14 0 2 4 6 8 10 12 14 16
R [fm]
Figure 4.16: Absolute value of the d + A scattering wave function using the ADWA theory
when nonlocal and local potentials are used. (a) d+48 Ca and (b) d+208 Pb. Both for the
L = 6 and J = 5 partial wave at Ed = 50 MeV in the laboratory frame. Figure reprinted
from [3] with permission.
The absolute values of the d + A scattering wave functions for the L = 6 and J = 5
partial wave are shown in Fig. 4.16 for d+48 Ca and d+208 Pb at Ed = 50 MeV. We see
similar features as we did for the L = 1 and J = 0 case: the wave function is both reduced
and pushed outward due to nonlocality. For the d+208 Pb case in particular, there is dramatic
shift in the wave function around the nuclear surface (∼ 7.5 − 8.5 fm), which we will find is
92
very important when calculating transfer cross sections.
For the calculation of (d, p) transfer cross sections, the nonlocal Perey-Buck potential was
used for the neutron and proton optical potentials in the entrance and exit channels. The
separate effects of nonlocality in the proton scattering state and the neutron bound state
have already been studied. The results of such a study has been published in our previous
papers [12, 14], and is not discussed here. In addition, since we have already determined
that the PCF is insufficient to take nonlocality into account, we did not investigate the PCF
in this study and focused instead on the effects of explicitly including nonlocality in the
In our analysis, we computed angular distributions for a wide variety of cases from 16 O
to 208 Pb. Some illustrative examples are shown in Fig. 4.17 and 4.18. Extensive tables for
all cases are shown in Tables 4.6, 4.7, and 4.8. In the tables we show the percent difference
between cross sections produced by nonlocal and local interactions, at the peak of the angular
distribution, relative to a purely local calculation. In the first column we include nonlocality
in all nucleon-target interactions. In the second (third) column we include nonlocality in the
In Fig. 4.17 and 4.18 we include the results when nonlocality is included consistently
(solid line), only in the deuteron channel (dashed line), only in the proton channel (dot-
dashed line), and where only a LPE potential is used (dotted line). In Fig. 4.17 we present
(d, p) calculations for deuterons impinging on: (a) 48 Ca at Ed = 10 MeV, (b) 132 Sn at
Ed = 10 MeV, and (c) 208 Pb at Ed = 20 MeV. The same cases are presented in Fig. 4.18,
but for Ed = 50 MeV. When available, we also present data points. The data in Fig. 4.17a
93
is published in arbitrary units. Therefore, this data set was normalized to the peak of the
Table 4.6: Percent difference of the (d, p) transfer cross sections at the first peak when
using nonlocal potentials in entrance and exit channels (1st column), nonlocal potentials
in entrance channel only (2nd column), and nonlocal potentials in exit channel only (3rd
column), relative to the local calculation with the LPE potentials, for a number of reactions
occurring at 10 MeV.
At low energies, nonlocality in the exit channel provides a significant enhancement of the
cross section for all cases, which is due to the neutron bound state. As mentioned before,
the ANC of the bound state resulting from nonlocal potentials is larger than that from local
potentials. Since low energy transfer reactions are primarily sensitive to the asymptotic
properties of the wave functions, this results in an increase of the cross section.
The nonlocality in the proton scattering state is not felt significantly at low energies,
so the reduction of the cross section due to the reduced amplitude of the proton scattering
state is small. This is consistent with the results published in our previous papers, [12, 14].
At higher energies, the nonlocality of the proton scattering state becomes more significant,
and there is a competition between the effects of nonlocality in the neutron bound state to
enhance the cross section, and the effects of nonlocality in the proton scattering state to
reduce the cross section. Nonlocality in the proton scattering state had a larger effect for
94
120
120
100 90
60
80 30
0
60 0 20 40 60
Nonlocal
Nonlocal Entrance
40 Nonlocal Exit
Local
a) Data
20
Data in arbitrary units
0
0 20 40 60 80 100
24
12
20
8
dσ/dΩ [mb/sr]
16 4
12 0
0 30 60 90 120
b)
8
4
0
0 30 60 90 120 150 180
15
6
12 4
2
9
0
0 30 60 90 120
6
3
c)
0
0 20 40 60 80 100 120
θc.m. [deg]
Figure 4.17: Angular distributions for (d, p) transfer cross sections. The insets are the
theoretical distributions normalized to the peak of the data distribution. (a) 48 Ca(d, p)49 Ca
at Ed = 10 MeV with data [16] at Ed = 10 MeV in arbitrary units. (b) 132 Sn(d, p)133 Sn at
Ed = 10 MeV with data [17] at Ed = 9.4 MeV. (c) 208 Pb(d, p)209 Pb at 20 MeV with data
[18] (Circles) and [19] (Squares) at Ed = 22 MeV. Figure reprinted from [3] with permission.
95
14
12 6
10 4
8 2
6 0
0 10 20 30 40 50
4
2 a)
0
30
dσ/dΩ [mb/sr]
25 Nonlocal
Nonlocal Entrance
20 Nonlocal Exit
Local
15 Data
10
5 b)
0
35
30
25
20
15
10
5 c)
0
0 10 20 30 40 50 60
θc.m. [deg]
Figure 4.18: Angular distributions for (d, p) transfer cross sections. The inset is the theo-
retical distributions normalized to the peak of the data distribution. (a) 48 Ca(d, p)49 Ca at
Ed = 50 MeV with data [20] at Ed = 56 MeV. (b) 132 Sn(d, p)133 Sn at Ed = 50 MeV. (c)
208 Pb(d, p)209 Pb at 50 MeV. Figure reprinted from [3] with permission.
96
Nonlocal Nonlocal
Final Nonlocal Entrance Exit Angle
Elab Bound Relative Relative Relative of
20 MeV State to Local to Local to Local Peak
16 O(d, p) 1d5/2 24.9% 2.6% 25.7% 0◦
16 O(d, p) 2s1/2 7.1% −0.7% 6.0% 0◦
40 Ca(d, p) 1f7/2 43.3% 11.0% 34.1% 26◦
48 Ca(d, p) 2p3/2 14.9% 7.1% 12.2% 8◦
126 Sn(d, p) 1h11/2 33.6% 7.7% 26.4 35◦
132 Sn(d, p) 2f7/2 3.2% 2.5% 4.2% 16◦
208 Pb(d, p) 2g9/2 35.0% 12.6% 20.5% 32◦
Table 4.7: Percent difference of the (d, p) transfer cross sections at the first peak when
using nonlocal potentials in entrance and exit channels (1st column), nonlocal potentials
in entrance channel only (2nd column), and nonlocal potentials in exit channel only (3rd
column), relative to the local calculation with the LPE potentials, for a number of reactions
occurring at 20 MeV.
the heavier nuclei due to a larger surface region being probed. This is seen in Fig. 4.18.
Comparing the dot-dashed lined with the dotted line, we see that there is an enhancement of
the cross section for 48 Ca, but a reduction for 132 Sn and 208 Pb. The net effect of nonlocality
in the exit channel depends on a complex interplay of the properties of the bound state (i.e.
number of nodes, binding energy, and orbital angular momentum), as well as the magnitude
of the real and imaginary parts of the scattering wave function near the nuclear surface.
Depending on the case, the shifting outward of the deuteron wave functions seen in Figs.
4.15 and 4.16 did not always have the same effect on the transfer cross sections. Comparing
the dashed and dotted lines in Figs. 4.17 and 4.18, we see that for 48 Ca, nonlocality in
the deuteron scattering state has a similar effect as for the proton scattering state in that
it reduces the cross section. As the size of the target increases, the outward shift of the
wave function becomes more important. This is seen in the comparison of the d+208 Pb and
the d+48 Ca wave functions in Figs. 4.15 and 4.16. The d+208 Pb wave function is shifted
outward more than the d+48 Ca wave function, which changes the amplitude at the nuclear
97
Nonlocal Nonlocal
Final Nonlocal Entrance Exit Angle
Elab Bound Relative Relative Relative of
50 MeV State to Local to Local to Local Peak
16 O(d, p) 1d5/2 22.3% 3.1% 15.8% 10◦
16 O(d, p) 2s1/2 20.7% 0.4% 21.2% 0◦
40 Ca(d, p) 1f7/2 4.8% 4.4% 0.2% 0◦
48 Ca(d, p) 2p3/2 41.9% −8.1% 39.9% 0◦
126 Sn(d, p) 1h11/2 6.9% 6.7% −2.5 13◦
132 Sn(d, p) 2f7/2 −10.9% 20.4% −26.2% 0◦
208 Pb(d, p) 2g9/2 64.8% 86.5% −1.7% 0◦
Table 4.8: Percent difference of the (d, p) transfer cross sections at the first peak when
using nonlocal potentials in entrance and exit channels (1st column), nonlocal potentials
in entrance channel only (2nd column), and nonlocal potentials in exit channel only (3rd
column), relative to the local calculation with the LPE potentials, for a number of reactions
occurring at 50 MeV. Figure reprinted from [3] with permission.
surface, and has a significant impact on the transfer cross section. As seen in Table 4.8,
nonlocality in the deuteron scattering state for 208 Pb(d, p)209 Pb has the most significant
The insets in Fig. 4.17 show that when the theoretical cross sections are normalized to
the data at the peak of the distribution, the low energy data cannot distinguish between
the various models since the shape of the theoretical distributions are similar. However, for
48 Ca(d, p)49 Ca at Ed = 50 MeV, nonlocality significantly improves the description of the
data. In all cases, if one were to extract a spectroscopic factor from the data, the results
including nonlocality would differ considerably from those when only local interactions are
used.
98
4.4.4 Comparing Distorted Wave Born Approximation and the
The DWBA is still the work-horse used in the analysis of most transfer cross sections. The
DWBA is based on a series expansion described in Sec. 2.2.1. This expansion is usually
truncated to first-order so that deuteron breakup is only included implicitly through the
imaginary part of the deuteron optical potential. This is unlike the ADWA which is based
on a three-body model, includes breakup explicitly, and relies on nucleon optical potentials.
Here we show that the differences in the DWBA and ADWA formalism can lead to very
In Fig. 4.19 the angular distributions for three different (d, p) reactions obtained using
the DWBA are compared to those obtained with the ADWA. There was no obvious way
to compare the effect of nonlocality in the entrance channel since DWBA and ADWA treat
the deuteron channel very differently, and a nonlocal global deuteron optical potential does
not exist. Both local and nonlocal potentials were used in the exit channel. For a useful
comparison, the same nonlocal and local potentials are used in the exit channel. For the
ADWA, the LPE potentials obtained from fits to the distribution generated with the Perey-
Buck potential are used, while for the DWBA we used the deuteron optical potential of
Daehnick [84].
We first focus on the local results, and compare in Fig. 4.19 the DWBA (dotted line) with
the ADWA (dashed line). The shapes are significantly different, as well as the magnitude of
the cross section at the first peak. Including nonlocality in the exit channel does not resolve
this discrepancy. We see that introducing nonlocality in the exit channel has the similar
effect of increasing the cross section for both the DWBA and ADWA calculations. We also
99
20
15
10
5
a)
0
ADWA Nonlocal Exit
12 ADWA Local Exit
dσ/dΩ [mb/sr]
3
b)
0
25
20
15
10
5 c)
0
0 10 20 30 40 50 60
θc.m. [deg]
Figure 4.19: Comparison of (d, p) transfer cross sections when using the DWBA as compared
to the ADWA. (a) 16 O(d, p)17 O, (b) 48 Ca(d, p)49 Ca with data from [20]. (c) 132 Sn(d, p)133 Sn.
All distributions at Ed = 50 MeV. Figure reprinted from [3] with permission.
100
show the in Fig. 4.19b. It is clear that for this case, the DWBA is not able to describe the
angular distribution from experiment. This is one example that demonstrates the need to
Since many reaction problems are solved in coordinate based theories, local interactions have
been preferred due to the simplicity of solving the equations. For this reason, Timofeyuk
and Johnson, [21, 22], developed a method to effectively include nonlocality in the deuteron
scattering state within the formalism of the ADWA. Their method relies on local potentials
so that the nonlocal equation does not need to be solved. Assuming the Perey-Buck form
for the nonlocal potential, and through expansions, they find that by shifting the energy at
which the local potentials are evaluated by ∼ 40 MeV from the standard Ed /2 value, one
can capture the effects of nonlocality. Since we are now able to include nonlocality explicitly,
As we are only concerned with nonlocality in the deuteron channel, we fixed the potentials
in the proton channel so we can make a meaningful comparison. In the exit channel, the LPE
potentials found from fits to Perey-Buck proton elastic scattering distributions were used,
along with the local binding potential used to reproduce the experimental binding energy.
In the entrance channel, we used the nonlocal Perey-Buck potential, and the corresponding
LPE potentials. To use the method of [21, 22], we needed an energy dependent local optical
potential. For this we used the CH89 potential [23] evaluated at the standard Ed /2 value,
and with the additional energy shift that was quantified in [21].
The results of this study are shown in Fig. 4.20. We show angular distributions for (d, p)
101
40
Nonlocal Entrance
Local Entrance
30 CH89 + Energy Shift
CH89
20
10 a)
4 b)
12
3
c)
0
0 20 40 60 80 100
θc.m. [deg]
Figure 4.20: Comparison of (d, p) angular distributions when using the energy shift method
of [21, 22]. (a) 16 O(d, p)17 O at Ed = 10 MeV (b) 40 Ca(d, p)41 Ca at Ed = 10 MeV (c)
208 Pb(d, p)209 Pb at E = 20 MeV. The solid line is when full nonlocality was included in
d
the entrance channel, dashed line is when the LPE potential was used, dot-dashed line when
the CH89 potential [23] was used with the additional energy shift quantified in [21], and the
dotted line when the CH89 potential was used at the standard Ed /2 value. Figure reprinted
from [3] with permission.
102
MeV. The solid line is the distribution form nonlocality explicitly included in the deuteron
channel, the dashed line are the results of the local calculations with the LPE potentials
for the deuteron scattering state, the dot-dashed line used the method of [21, 22], and the
dotted line is when the local potentials are evaluated at the standard Ed /2 value.
Our results show that the energy shift method always increases the cross section. How-
ever, the explicit inclusion of nonlocality in the deuteron scattering state can sometimes
decrease the cross section, as is seen in Fig. 4.20a. Often times, the energy shift moves the
transfer distribution towards larger angles, an effect also seen in the full nonlocal calcula-
tions. In some cases, the energy shift method over shoots the full nonlocal calculation, as
in Fig. 4.20b. In Fig. 4.20c, we see an example where the energy shift does a very good
job at reproducing the nonlocal effects. In general, we found that for most cases, the energy
shift captured the qualitative effects of nonlocality, but was unable to provide an accurate
4.4.6 Summary
In this work we studied the effects of nonlocality on (d, p) transfer reactions. An extension of
the ADWA theory was developed to include nonlocality in the deuteron scattering state using
the Perey-Buck nonlocal nucleon optical potential [1]. In the exit channel the Perey-Buck
potential was used to describe the proton scattering state, and its real part was adjusted
for the neutron bound state. For the scattering, a local phase equivalent (LPE) potential
was obtained by fitting the elastic scattering generated from the corresponding nonlocal
potential. Both the local and nonlocal bound states reproduced the same experimental
binding energies.
For the (d, p) reactions studied, we found that the inclusion of nonlocality in both the
103
entrance and exit channels increased the transfer cross section by ∼ 40%. In most cases,
nonlocality in the deuteron scattering state caused a moderate increase in the transfer cross
section. However, for heavy targets at high energies, this increase was large. Nonlocality
in the exit channel caused, almost exclusively, an increase in the transfer cross section,
except for heavy targets at high energies for which the cross sections were reduced. We also
compared our ADWA result with those from DWBA and found the effects of nonlocality in
the final state to be consistent in both formulations, even if quantitatively different. We also
compared our ADWA results with the energy shift method introduced by Timofeyuk and
Johnson [21, 22] and found that method to be qualitatively consistent with our results.
The conclusion of the present study confirm those of [12, 14]. There are important
differences in the transfer cross sections when including nonlocality explicitly as compared to
when using LPE potentials. This highlights the necessity of explicitly including nonlocality
to describe transfer reactions. Since the inclusion of nonlocality normally increases the
cross section, a re-analysis of transfer reaction data will likely reduce currenlty accepted
104
Chapter 5
5.1 Conclusions
In this thesis we studied the effect of nonlocality of the optical potential in transfer reactions.
For this purpose we developed a method for solving the integro-differential equations and
extended the adiabatic distorted wave approximation for transfer (d, p) reactions to include
For the (p, d) reaction study using the Perey-Buck nonlocal potential of Sec. 4.2, we
considered a range of nuclei, and proton energies of Ep = 20 and 50 MeV [12]. We calculated
the transfer matrix element in the distorted wave Born approximation (DWBA), and used a
local optical potential to describe the deuteron scattering state. We found that the explicit
inclusion of nonlocality increased the transfer cross section at the first peak by 15 − 35%,
relative to when local potentials were use. We found that in all cases, the Perey correction
factor traditionally used does not provide a quantitative description. Our results suggest
that such a correction factor to account for nonlocality should not be used.
In Sec. 4.3 we compared the DOM potential and the Perey-Buck potential to study (p, d)
reactions on 40 Ca at Ep = 20, 35, and 50 MeV [14]. We included nonlocality in the entrance
channel, then computed transfer cross sections in the DWBA, ignoring nonlocality in the
deuteron channel. We generated two local phase equivalent (LPE) potentials, one for the
105
DOM and one for the Perey-Buck potential. Both the DOM and the Perey-Buck potential
produced very large increases in the magnitude of the transfer cross section, ≈ 15 − 50% for
the DOM potential, and ≈ 30 − 50% for the Perey-Buck potential. Like in the first study
of (p, d) reactions, when nonlocality was included only in the bound state, large increases in
the magnitude of the transfer cross section were seen. Typically, nonlocality in the proton
scattering state acts in the opposite direction, reducing the transfer cross section.
In the last study of Sec. 4.4 involving (d, p) reactions, nonlocality was included in the
deuteron channel within the adiabatic distorted wave approximation (ADWA), which unlike
the DWBA, includes deuteron breakup explicitly [3]. The formalism for the local ADWA
theory had to be extended to include nonlocal potentials, and was done in Sec. 2.6. We
found that the inclusion of nonlocality increased the transfer cross section by ∼ 40%. In most
cases, nonlocality in the deuteron scattering state caused a modest increase in the transfer
cross section. However, for heavy targets at high energy, the increase due to nonlocality in
the deuteron channel was very large. This is in contrast to when nonlocality is added to the
proton scattering state, which often times reduced the transfer cross section. The reason
for the difference is that nonlocality in the deuteron scattering state had both the effect of
reducing the magnitude of the scattering wave function within the nuclear interior, but also
All these three studies [12, 14, 3] demonstrate the important differences in the transfer
cross sections when including nonlocality explicitly, as compared to when using local phase
Most often, transfer reactions are performed to extract a spectroscopic factor. The analy-
sis is done using local potentials. As explicit inclusion of nonlocality increased the predicted
106
cross sections, one would expect lower spectroscopic factors to result from the analysis if
nonlocal potentials are used. This may well contribute to solving the discrepancy between
spectroscopic factors extracted from knockout and transfer [112], but that study has yet to
be performed.
5.2 Outlook
Going forward, nonlocality must be carefully taken into account in any advanced reaction
theory. It will become increasingly important to construct a modern nonlocal global optical
potential. Work along these lines has already been done by Tian, Pang, and Ma (TPM)
with the recent publication of their nonlocal potential, [2]. However, this potential is still
On physical grounds, the optical potential must be energy dependent due to channel
coupling effects. Preliminary indications from an unpublished study by Bacq, Lovell, Titus,
and Nunes [113] show that there is indeed an energy dependence in nonlocal potentials of the
Perey-Buck form. Work on specifying the precise energy dependence is ongoing. It would
The Perey-Buck form for the potential comprises a single nonlocality parameter, β. This
simple form has been useful for many decades since it allowed for simple implementation.
It is unlikely that a single nonlocality parameter is able to represent the complex nature
of nonlocality in the realistic many-body problem. The DOM potential, for example, has
a different β for each term of the potential. A more sophisticated form for the nonlocal
107
potential should be considered.
There are several methods available to construct a microscopically based optical poten-
tial. Some of these methods were discussed in Chapter 1, and should be pursued now that it
is known that nonlocality is very important and must be included explicitly. Since previous
microscopic calculations of nonlocal potentials have shown that their form does not resem-
ble the simple Gaussian nonlocality of the Perey-Buck potential, it is important to better
With an improved nonlocal global optical potential, existing transfer data can be re-
analyzed. The large discrepancies between spectroscopic factors extracted from the nonlocal
and local calculations in this study demonstrates that the structure information of most
nuclei are likely to be altered when the data is analyzed with nonlocal potentials.
While this thesis focused on (d, p) and (p, d) reactions, the role of nonlocality should be
investigated in other reactions as well. We are currently investigating the role of nonlocality
in (d, n) reactions [114]. Surprisingly, the effects of nonlocality in (d, n) reactions are even
Along with transfer, there are many other reactions that are performed to extract sin-
gle particle structure of nuclei. Nuclear knockout reactions (A(a, bγ)X) are an alternative
method to extract a spectroscopic factor. Such a reaction also requires an optical potential
between the colliding nuclei. Understanding the effect of nonlocality in this case is also im-
portant. Inelastic scattering provides the transition strength between the ground state and a
bound excited state in a nucleus. Since these are obtained by comparing experimental data
to theoretical distributions resulting from a DWBA or coupled channel analysis, one may
again expect that inclusion of nonlocality in the description of the process will improve the
reliability of the extracted transition strength. Three-body models exist to calculate trans-
108
fer, such as the continuum discretized coupled channel method [95]. In this case, coupled
channel integro-differential equations would need to be solved, and new numerical methods
may be needed to accomplish this task. Finally, charge-exchange reactions, such as (p, n),
probe the spin and isospin properties of nuclei. Once again, information from experiments
are sometimes extracted from a DWBA analysis, and including nonlocality in the optical
potential will improve the reliability of theoretical predictions, consequently the reliability
One aspect that was not addressed in this thesis concerns error quantification. Lovell
and Nunes [115] are currently addressing error quantification in direct reaction theories when
using local potentials. An extension of that work will be necessary as the use of nonlocal
nonlocal potentials becomes more widespread. Much work has highlighted the importance
of nonlocality in the reaction dynamics. How then can we constrain nonlocality from data?
The Perey-Buck potential, for example, has a Gaussian form with a nonlocality range of
in this way may not be the best method. Reactions sensitive to short range correlations may
The results of this thesis, along with the advent of microscopic theories to construct
nonlocal potentials, improved phenomenological nonlocal global optical potentials, and ever
increasing computer power, has the potential to elevate reaction theory to a new level.
While we have focused on (p, d) and (d, p) reactions, with our increased understanding of
nonlocality, it has become necessary to update the other formalisms and codes commonly
109
APPENDICES
110
Appendix A
In order to assess the validity of the local approximation we need to solve Eq.(3.14) exactly.
Several methods exist for solving the scattering state, such as [78, 80]. Our approach follows
Perey and Buck [1] where Eq.(3.14) is solved by iteration. For simplicity, we will drop the
local part of the nonlocal potential, Uo (r), in our discussion, although it is included in our
calculations.
To solve the partial wave equation of Eq.(3.14) numerically, we need to find the kernel
function gL (R, R0 ). In order to do this, we first need to do a partial wave expansion of the
nonlocal potential,
X 2L + 1 gL (R, R0 )
UPNBL (R, R0 ) = PL (cos θ), (A.1)
4π RR0
L
where we defined θ as the angle between R and R0 . Inserting the Perey-Buck form for the
nonlocal potential, multiplying both sides by PL (cos θ), integrating over all angles, using the
orthogonality of the Legendre polynomials, and solving for gL (R, R0 ), we find that
R−R 0 2
Z 1 exp − β
0 0 1 0
gL (R, R ) = 2πRR U (R + R ) 3
PL (cos θ)d(cos θ). (A.2)
2 −1 π 2 β3
111
For a moment, consider only the integral on the right:
R−R 0 2 R 2 +R02
Z 1 exp − exp − Z 1
β β2
3
PL (cos θ)d(cos θ) = 3
ei(−izcosθ) PL (cos θ)d(cos θ).
−1 π 2 β3 π 2 β3 −1
(A.3)
Z 1
1
jL (x) = L eixu PL (u)du, (A.4)
2i −1
2iL z R2 + R02
1
gL (R, R0 ) = 1
jL (−iz)exp − 2
U 0
(R + R )
β 2
π2β
0 1 0
= hL (R, R )U (R + R ) . (A.5)
2
the terms, so we need to approximate this function when doing numerical calculations for
large values of the argument. Making use of the asymptotic expression for the spherical Bessel
0 2
function, and neglecting terms proportional to exp − R+R
β , we find that hL (R, R0 ) for
z 1 can be approximated as
112
2
0
1 − R−R
β
hL (R, R0 ) ≈ 1
e for |z| 1. (A.6)
π2β
Scattering solutions are considered first, where the subscript n denotes the nth order
approximation of the correct solution. The iteration scheme starts with an initialization,
~2 d2
L(L + 1)
− χn=0 (R) + [E − Uinit (R)]χn=0 (R) = 0, (A.7)
2µ dr2 R2
where Uinit (R) is some suitable local potential used to get the iteration process started.
~2 d2
L(L + 1)
− χn (R) + [E − Uinit (R)]χn (R)
2µ dr2 R2
Z
= gL (R, R0 )χn−1 (R0 )dR0 − Uinit (R)χn−1 (R), (A.8)
with as many iterations necessary for convergence. The number of iterations depends mostly
on the partial wave being solved for (lower partial waves require more iterations) and the
quality of Uinit (R). It was rare for any partial wave to require more than 20 iterations to
converge, even with a very poor choice for Uinit (R). If the LPE potential is used as Uinit (R),
For the bound state problem, the method is somewhat different. A variety of methods
exist in the literature, some developed specifically to handle non-analytic forms (e.g. [81]).
Our approach may not be the most efficient, but it is straightforward, general, and easy to
implement. To solve the bound state problem with a nonlocal potential we begin by solving
113
Eq.(A.7). Since we are using the wave function from the previous iteration to calculate the
nonlocal integral, we need to keep track of the different normalizations of the inward and
outward wave functions that results from the choice for the initial conditions for each wave
~ 2 d2
`(` + 1) In
2
− 2
φn (r) + [E − Uinit (r)]φIn
n (r)
2µ dr r
Z R
M ax
= g` (r, r0 )φIn 0 0 In
n−1 (r )dr − Uinit (r)φn−1 (r), (A.9)
0
and
~2 d2
`(` + 1) Out
− φn (r) + [E − Uinit (r)]φOut
n (r)
2µ dr2 r2
Z R
M ax
= g` (r, r0 )φOut 0 0 Out
n−1 (r )dr − Uinit (r)φn−1 (r), (A.10)
0
where RM ax is some maximum radius chosen greater than the range of the nuclear interac-
tion. Note that φIn (r) is the wave function for integrating from the edge of the box inward
and has a normalization set by the Whittaker function as the initial condition, while φOut (r)
is the wave function for integrating from the origin outward and has the normalization set
Even though φOut and φIn differ by only a constant, these two equations (Eq.(A.9) and
Eq.(A.10)) are necessary because the value of the normalization constant is only known after
convergence. For a given iteration, φOut and φIn converge when their logarithmic derivatives
agree at the matching point. To keep the proper normalization throughout the entire range
114
[0, RM ax ], we need to retain two versions of the converged wave function for each iteration:
CnOut φOut
n (r) for 0 ≤ r < RMatch
In
φn (r) = (A.11)
φIn (r)
for RMatch ≤ r ≤ RMax
n
φOut
n (r) for 0 ≤ r < RMatch
φOut
n (r) = (A.12)
C In φIn (r)
for RMatch ≤ r ≤ RMax
n n
where
φOut(In) (RMatch )
C In(Out) = (A.13)
φIn(Out) (RMatch )
The full iteration scheme is converged when the binding energy obtained from the previous
iteration agrees with the binding energy from the current iteration within a desired level of
accuracy. Again, although this may not be the most efficient method, it is general (whatever
the form of nonlocality) and is very stable, providing a good option for future studies beyond
115
Appendix B
Here we provide details on the derivation of the PCF, Eq. (3.18). We also include the
derivation of the transformation formulas, as well as the correct radial version of the trans-
formation formulas which could be used to transform the nonlocal radius and diffuseness to
We start from Eq. (3.15). Let us define a function F (R) that connects the local wave
function ΨLoc (R), resulting from the potential U LE (R), with the wave function resulting
Since the local and nonlocal equations describe the same elastic scattering, the wave
By inserting Eq.(B.1) into the nonlocal equation Eq.(3.15) we can reduce the result to the
~2 2 Loc
− ∇ Ψ (R) + U LE (R)ΨLoc (R) = EΨLoc (R), (B.2)
2µ
116
2 ~2 (∇2 F (R))ΨLoc (R)
− ~µ ∇F (R) · ∇ΨLoc (R) − 2µ
U LE (R) =
F (R)ΨLoc (R)
U N L (R, R0 )F (R0 )ΨLoc (R0 )dR0
R
+ + Uo (R).
F (R)ΨLoc (R)
(B.3)
We next consider the numerator of the second term of Eq.(B.3) and introduce the explicit
Z
NL 1
UW S R − s H(s)F (R − s)ΨLoc (R − s)ds
2
Z Z
1
≈ UW S (R)F (R) H(s)Ψ (R − s)ds − F (R)∇UW S (R) · sH(s)ΨLoc (R − s)ds
N L Loc N L
Z 2
N L (R)∇F (R) ·
− UW SH(s)ΨLoc (R − s)ds, (B.4)
S
where
s2
exp − 2
β
H(s) = . (B.5)
π 3/2 β 3
117
2
1 ~
U LE (R) ≈ Loc
− (∇F · ∇ΨLoc )
F (R)Ψ (R) µ
Z
+ UW N L (R)F (R) H(s)ΨLoc (R − s)ds
S
Z
1
− F (R)∇UW S · sH(s)ΨLoc (R − s)ds
N L
2
Z
N L Loc
− UW S (R)∇F · sH(s)Ψ (R − s)ds
~2 ∇2 F (R)
− + Uo (R). (B.6)
2µ F (R)
Consider the four terms in the brackets. All of these terms are divided by ΨLoc , which
has nodes. The first, third, and fourth terms depend on dot products and gradients of ΨLoc .
These terms are unlikely to individually equal zero when ΨLoc in the denominator equals
zero. Thus, we require that these terms sum to zero so that U LE (R) remains finite. As
pointed out in [74], this is not an approximation, but merely a condition for the method to
"R #
H(s)ΨLoc (R − s)ds ~2 ∇2 F (R)
U LE (R) = UW
N L (R)
S + Uo (R) − (B.7)
ΨLoc (R) 2µ F (R)
~2
Z
1
0 = (∇F · ∇Ψ ) + F (R)∇UW S + UW S (R)∇F · sH(s)ΨLoc (R − s)ds.
Loc N L N L
µ 2
(B.8)
Instead of using the local WKB approximations as Austern did [74], we use the operator
118
form of the Taylor expansion to factorize the wave function:
with k = −i∇. This simplifies the integrals in Eq.(B.7) and Eq.(B.8). Consider first the
integral in Eq.(B.7)
Z Z
H(s)ΨLoc (R − s)ds = e−is·k H(s)ds ΨLoc (R)
−k 2 β 2
= exp ΨLoc (R). (B.10)
4
Therefore, assuming the potentials are scalar functions of R, and replacing Eq.(B.10) into
Eq.(B.7) we obtain;
µβ 2 ~2 ∇2 F (R)
U LE (R) = N L (R) exp
UW S
LE
− 2 E − U (R) + Uo (R) − , (B.11)
2~ 2µ F (R)
where we used k 2 = −∇2 in the exponent to first order, and the Schrödinger’s equation.
formula
~2 ∇2 F (R) µβ 2
N L (R)
UW = Loc (R) +
UW exp Loc
E − UW S (R) − Uo (R) .
S S 2µ F (R) 2~2
The ∇2 F term is significant around the surface, but near the origin this term is negligible.
119
Therefore, if we neglect this term, then we must remove the radial arguments, and consider
µβ 2
N L (0)
UW ≈ Loc (0) exp
UW Loc
E − UW S (0) − Uo (0) . (B.12)
S S 2~2
The UW S (R) functions are of a Woods-Saxon form, and have real and imaginary parts
Near the origin, UILoc ≈ 0 so this term can be neglected in the exponent, and UR ≈ −Vv .
While the spin-orbit term diverges at the origin, it rapidly goes to zero away from the origin,
Coulomb potential at the origin for a uniform sphere of charge. Taking the real part of the
120
µβ 2
VvN L = VvLoc exp E − Vc + VvLoc . (B.15)
2~2
µβ 2
UIN L (R) = UILoc (R) exp E − Vc + VvLoc . (B.16)
2~2
While UI (R) ≈ 0 near the origin, the local and nonlocal terms have the same form factor, so
the form factors exactly cancel as long as the radius and diffuseness are identical. Therefore,
µβ 2
WdN L = WdLoc exp E − Vc + VvLoc , (B.17)
2~2
It is important to note that these equations are only valid for transforming the depths of the
potentials, thus Eq.(B.12) should not be used while retaining the radial dependence. Indeed,
Now consider the integral in Eq.(B.8). Using Eq.(B.9) to expand the wave function, and
~2
Loc 1 N L N L
0 = (∇F · ∇Ψ ) + F (R)∇UW S + UW S (R)∇F
µ 2
Z
× s cos(θ)H(s)e −is·k ds ΨLoc (R). (B.18)
121
~2
1 N L N L
0 = ∇F − F (R)(∇UW S ) + UW S (R)(∇F )
µ 2
β 2
µβ 2
× exp − 2 E − U (R) .LE (B.19)
2 2~
If we assume that the local momentum approximation is valid, this equation can be solved
µβ 2 N L
µβ 2 − 12
LE
F (R) = 1 − 2 UW S (R) exp − 2 E − U (R) . (B.20)
2~ 2~
If the local momentum approximation is not valid, then insertion of Eq.(B.20) into the rhs
of Eq.(B.19) will deviate from zero by a term related to the derivative of U LE (R). This
additional term will be significant at the surface, and thus one can expect discrepancies in
1
µβ 2
~ 2 ∇2 F (R) − 2
F (R) = 1 − 2 U LE (R) − Uo (R) + . (B.21)
2~ 2µ F (R)
Neglecting the term containing ∇2 F gives us Eq.(3.18), which is the correction factor of
Austern [74]. The contribution of ∇2 F/F is only important at the surface, and again it is
precisely for these radii that discrepancies can be expected in applying Eq.(3.18).
122
Appendix C
Here we will derive the nonlocal adiabatic potential. We begin with a three-body Schrödinger
Equation:
h i
T̂R + Tr + Vnp (r) + ÛnA + ÛpA − E Ψ(r, R) = 0. (C.1)
∞
X
Ψ(r, R) = Φi (r)Xi (R), (C.2)
i=0
Noting that
123
h i h i
T̂R − Ed Φ(r)X(R) = − ÛnA + ÛpA Φ(r)X(R). (C.5)
Now we evaluate the nucleon nonlocal operator, ÛN A Φ(r)X(R). For a moment, consider
just the neutron potential (with Rp,n = R ± 2r where the “+” sign is for the proton and the
r 0 r0
Z
= J U R − ,R − Ψ(r0 , R0 )δ(R0p − Rp )dr0 dR0
2 2
r 0 r0 r0
Z
r
= J U R − ,R − Ψ(r , R )δ R + − (R + ) dr0 dR0
0 0 0
2 2 2 2
r0
Z
r
U R − , R0 − Ψ(r0 , R0 )δ r0 − (r − 2(R0 − R)) dr0 dR0
= 8J
2 2
∂R0n ∂R0n
1 − 12
J = ∂R0 ∂r0
= = 1, (C.7)
∂R0p ∂R0p 1
1
∂R0 ∂r0 2
which gives us
r 0 r − 2(R0 − R)
Z
Ψ r − 2(R0 − R), R0 dR0
ÛnA Ψ(r, R) = 8 UnA R − , R −
2 2
Z r r
= 8 UnA R − , 2R0 − R − Ψ(r − 2(R0 − R), R0 )dR0 . (C.8)
2 2
124
For the nucleon nonlocal operator, we have
Z r r
ÛN A Φ(r)X(R) = 8 UN A R ± , 2R0 − R ± Φ(r ± 2(R0 − R))X(R0 )dR0 .
2 2
(C.9)
get
r r r
UN A R ± , 2R0 − R ± = UN A Rp,n , 2R0 − R − R + R ±
2 2 2
= UN A Rp,n , 2(R0 − R) + Rp,n
= UN A Rp,n , Rp,n + 2s , (C.10)
where we made the definition, s = R0 − R. Since dR0 → ds, we see that ÛN A Φ(r)X(R)
becomes
Z
ÛN A Φ(r)X(R) = 8 UN A Rp,n , Rp,n + 2s Φ(r ± 2s)X(R + s)ds. (C.11)
The general expansion of the full wave function for a given partial wave with total angular
125
X n o
J M
Ψ(r, R) ≈ Φ(r)X(R) = Φ(r) ⊗ ỸL0 (R̂) 0 ⊗ ΞIt (ξt ) χ T0 0 T (R)
Jp L Jp
L0 Jp0 J T MT
( )
X n o
= ΞI (ξnp ) ⊗ Ỹ`0 (r̂) ⊗ ỸL0 (R̂) ⊗ ΞIt (ξt )
d Jd Jp0
`0 L0 Jp0 J T MT
J M
χ T0 0 T (R)
L Jp
× φ`0 (r) , (C.12)
R
where
n o
ΞI (ξnp ) = ΞIp (ξp ) ⊗ ΞIn (ξn ) . (C.13)
d I d
In these equations, ΞIp (ξp ), ΞIn (ξn ) and ΞIt (ξt ) are the spin functions for the proton, neu-
tron, and target, respectively. Ip = In = 1/2 are the spin of the proton and neutron
respectively, and It is the spin of the target. ỸL (R̂) is the spherical harmonic for the orbital
motion between the projectile and target, while Ỹ` (r̂) is for the internal orbital angular mo-
mentum of the deuteron. We are using the phase convention where there is a built in factor
of iL so that ỸL (R̂) = iL YL (R̂) where YL (R̂) is defined on [93], p.133, Eq.(1). The spin of
the deuteron is given by Id = 1 and the total angular momentum of the deuteron is Jd = 1.
The total angular momentum of the deuteron is coupled to the orbital angular momentum
between the deuteron and target to give a total angular momentum of the projectile, Jp .
The total angular momentum of the projectile is coupled to the spin of the target to give
126
h i
T̂R − Ed Φ(r)X(R) = −ÛN A Φ(r)X(R)
Z
= −8 UN A Rp,n , Rp,n + 2s Φ(r ± 2s)X(R + s)ds
X Z
= − 8 UN A Rp,n , Rp,n + 2s
`0 L0 Jp0
( )
n o
× ± 2s)
ΞI (ξnp ) ⊗ Ỹ`0 (r\ ⊗ ỸL0 (R
\ + s) ⊗ ΞIt (ξt )
d Jd Jp0 JT MT
J M
χ T0 0 T (|R + s|)
L Jp
× φ`0 (|r ± 2s|) ds. (C.14)
|R + s|
We would like to do a partial wave decomposition to get an equation for each LJ combination
( )∗
X n o
ΞI (ξnp ) ⊗ Ỹ` (r̂) ⊗ ỸL (R̂) ⊗ ΞIt (ξt ) φ` (r)Vnp (r), (C.15)
d Jd
` Jp JT MT
and integrate over dr, dΩR , dξnp , and dξt , where dξnp = dξn dξp . Consider first just the lhs
127
Z X (n )∗
o
ΞI (ξnp ) ⊗ Ỹ` (r̂) ⊗ ỸL (R̂) ⊗ ΞIt (ξt ) φ` (r)Vnp (r)
d Jd
` Jp JT MT
( )
h i X n o
× T̂R − Ed ΞI (ξnp ) ⊗ Ỹ`0 (r̂) ⊗ ỸL0 (R̂) ⊗ ΞIt (ξt )
d Jd Jp0
`0 L0 Jp0 JT MT
J M
χ T0 0 T (R)
L Jp
× φ`0 (r) drdΩR
R
X X 1 ~2 ∂2 L0 (L0 + 1)
Z
J T MT
= − 2
− 2
+ Ed χ 0 0 (R) φ` (r)Vnp (r)φ`0 (r)r2 dr
R 2µ ∂R R L Jp
0
` `0 L0 Jp
Z (n )∗
o
× ΞI (ξnp ) ⊗ Ỹ` (r̂) ⊗ ỸL (R̂) ⊗ ΞIt (ξt ) (C.16)
d Jd Jp JT MT
( )
n o
× ΞI (ξnp ) ⊗ Ỹ`0 (r̂) ⊗ ỸL0 (R̂) ⊗ ΞIt (ξt ) dΩr dΩR dξnp dξt .
d Jd Jp0 JT MT
Z (n )∗
o
ΞI (ξnp ) ⊗ Ỹ` (r̂) ⊗ ỸL (R̂) ⊗ ΞIt (ξt )
d Jd Jp J M
( )T T
n o
× ΞI (ξnp ) ⊗ Ỹ`0 (r̂) ⊗ ỸL0 (R̂) ⊗ ΞIt (ξt ) dΩr dΩR dξnp dξt
d Jd Jp0 JT M T
X Z n o ∗
= ΞI (ξnp ) ⊗ Ỹ` (r̂) ⊗ ỸL (R̂)
d Jd Jp Mp
Mp Mp0
n o
× ΞI (ξnp ) ⊗ Ỹ`0 (r̂) ⊗ ỸL0 (R̂) dΩr dΩR dξnp
d Jd Jp0 Mp0
Z
J M J M
Ξ∗I µ (ξt )ΞI µ0 (ξt )dξt .
X
× CJ TM TI µ C T0 0T 0 (C.17)
p p t t Jp Mp It µt t t t t
µ0t µt
Ξ∗I µ (ξt )ΞI µ0 (ξt )dξt = δµ µ0 , summing over µ0t , and evaluating
R
Evaluating the integral,
t t t t t t
128
the complex conjugate using [93], p.62, Eq.(6), gives
X X J M Z n ∗
JT MT
o
= T T
CJ M I µ C 0 0 ΞI (ξnp ) ⊗ Ỹ` (r̂) ⊗ ỸL (R̂)
p p t t Jp Mp It µt d Jd J p Mp
Mp Mp0 µt
n o
× ΞI (ξnp ) ⊗ Ỹ`0 (r̂) ⊗ ỸL0 (R̂) dΩr dΩR dξnp
d Jd Jp0 Mp0
X X J M J M
= CJ TM TI µ C T0 0T
p p t t Jp Mp It µt
Mp Mp0 µt
Z n o
× (−)Jp −Mp ΞI (ξnp ) ⊗ Ỹ` (r̂) ⊗ ỸL (R̂)
d Jd Jp ,−Mp
n o
× ΞI (ξnp ) ⊗ Ỹ`0 (r̂) ⊗ ỸL0 (R̂) dΩr dΩR dξnp . (C.18)
d Jd Jp0 Mp0
n o n o
ΞI (ξnp ) ⊗ Ỹ` (r̂) ⊗ ỸL (R̂) ΞI (ξnp ) ⊗ Ỹ`0 (r̂) ⊗ ỸL0 (R̂)
d Jd d Jd
Jp ,−Mp Jp0 Mp0
(
X SM n o
= C S ΞI (ξnp ) ⊗ Ỹ` (r̂) ⊗ ỸL (R̂)
Jp ,−Mp Jp0 Mp0 d Jd
SMS Jp
n o )
⊗ ΞI (ξnp ) ⊗ Ỹ`0 (r̂) ⊗ ỸL0 (R̂)
d Jd Jp0 SM
( S
n o
00
→ CJ ,−M J 0 M 0 ΞI (ξnp ) ⊗ Ỹ` (r̂) ⊗ ỸL (R̂)
p p p p d Jd Jp
n )
o
⊗ ΞI (ξnp ) ⊗ Ỹ`0 (r̂) ⊗ ỸL0 (R̂) . (C.19)
d Jd Jp0 00
129
n o
1 00 o n
C ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞI (ξnp ) ⊗ ỸL0 (R̂) 0 . (C.20)
4π Jp ,−Mp Jp0 Mp0 d Jp d Jp 00
X 1 ~2 ∂ 2 L0 (L0 + 1)
Z
JT M T
− 2
− 2
+ Ed χ 0 0 (R) φ` (r)Vnp (r)φ`0 (r)r2 dr
R 2µ ∂R R L Jp
0L0 Jp
X X J M J M 1 00
× CJ TM TI µ C T0 0T C (−)Jp −Mp (C.21)
p p t t Jp Mp It µt 4π Jp ,−Mp Jp0 Mp0
Mp Mp0 µt
Z n o n o
× ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞI (ξnp ) ⊗ ỸL0 (R̂) 0 dΩr dΩR dξnp .
d Jp d Jp 00
The integral over dΩr gives 4π and cancels the 1/4π already there, which leaves us with
X 1 ~2 ∂ 2 L0 (L0 + 1)
Z
JT M T
− 2
− 2
+ Ed χ 0 0 (R) φ` (r)Vnp (r)φ`0 (r)r2 dr
R 2µ ∂R R L Jp
0L0 Jp
X X J M J M
× CJ TM TI µ C T0 0T CJ00,−M J 0 M 0 (−)Jp −Mp
p p t t Jp Mp It µt p p p p
Mp Mp0 µt
Z n o n o
× ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞI (ξnp ) ⊗ ỸL0 (R̂) 0 dΩR dξnp (C.22)
d Jp d Jp 00
δJ J 0 δM ,M 0
p p p p
CJ00p ,−Mp Jp Mp = (−)Jp +Mp . (C.23)
Jˆp
130
Inserting Eq.(C.23) into Eq.(C.22), and summing over Jp0 and Mp0 (using [93], p.236, Eq.(8))
we get,
X 1 ~2 ∂ 2 L0 (L0 + 1)
Z
JT M T
− 2
− 2
+ Ed χ 0 (R) φ` (r)Vnp (r)φ`0 (r)r2 dr
R 2µ ∂R R L Jp
0L
XX J M J M 1
× CJ TM TI µ CJ TM TI µ (−)2Jp (C.24)
Mp µt
p p t t p p t t Jˆp
Z n o n o
× ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞI (ξnp ) ⊗ ỸL0 (R̂) dΩR dξnp
d Jp d Jp 00
X 1 ~2 ∂ 2 L0 (L0 + 1)
Z
JT M T
= − 2
− 2
+ Ed χ 0 (R) φ` (r)Vnp (r)φ`0 (r)r2 dr
R 2µ ∂R R L Jp
L0
Z n o
2J p 1 o n
× (−) ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞI (ξnp ) ⊗ ỸL0 (R̂) dΩR dξnp .
Jˆp d Jp d Jp 00
Coupling the spin functions together and the spherical harmonics together, each up to zero
angular momentum, using [93], p.70, Eq. (11), and p.358 Eq.(4),
n o n o
ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞI (ξnp ) ⊗ ỸL0 (R̂)
d Jp d Jp 00
Id L Jp
nn o n o o
ˆ2
= Jp Id L0 Jp ΞI (ξnp ) ⊗ ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ỸL0 (R̂)
d d 0 0 00
0 0 0
Jˆp δLL0 nn o n o o
= ΞI (ξnp ) ⊗ ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ỸL0 (R̂) . (C.25)
Iˆd L̂ d d 0 0 00
131
1 ~2 L0 (L0 + 1)
2 Z
∂ JT MT
− − + Ed χ 0 (R) φ` (r)Vnp (r)φ`0 (r)r2 dr
R 2µ ∂R2 R2 L Jp
Z n Z n
2Jp 1 o 1 o
× (−) ΞI (ξnp ) ⊗ ΞI (ξnp ) dξnp ỸL (R̂) ⊗ ỸL (R̂) dΩR
Iˆd d d 00 L̂ 00
(C.26)
The integral over the two Weinberg states multiplied by Vnp gives −1 by the normalization
1 ~2 L0 (L0 + 1)
2
∂ J M
2
− 2
+ Ed χ T0 T (R) (C.27)
R 2µ ∂R R L Jp
Z n Z n
1 o 1 o
× ΞI (ξnp ) ⊗ ΞI (ξnp ) dξnp ỸL (R̂) ⊗ ỸL (R̂) dΩR .
Iˆd d d 00 L̂ 00
The integral over the spin functions can be worked out with Eq.(11), p.70, and Eq.(4), p.358
of [93]:
Z n o
ΞI (ξnp ) ⊗ ΞI (ξnp ) dξnp
d d 00
Z nn o n o o
= Ξ1/2 (ξn ) ⊗ Ξ1/2 (ξp ) ⊗ Ξ1/2 (ξn ) ⊗ Ξ1/2 (ξp ) dξn dξp
1 1 00
1 1 1
Z
2 2
nn o n o o
= 1̂2 1 1 Ξ1/2 (ξn ) ⊗ Ξ1/2 (ξn ) ⊗ Ξ1/2 (ξp ) ⊗ Ξ1/2 (ξp ) dξn dξp
2 2 1 0 0 00
0 0 0
1̂2
Z Z n o Z n o
= 2
Ξ (ξ
1/2 n ) ⊗ Ξ (ξ )
1/2 n 00 n dξ Ξ (ξ
1/2 p ) ⊗ Ξ1/2 p 00 dξp .
(ξ ) (C.28)
1̂ 1̂
2
132
Doing the integral over the neutron spin functions by using Eq.(13), p.132, and Eq.(1), p.248,
Z n o
Ξ1/2 (ξn ) ⊗ Ξ1/2 (ξn ) dξn
00
Z
00
X
= C1/2,−µ Ξ (ξ )Ξ (ξ )dξn
n ,1/2,µn 1/2,−µn n 1/2,µn n
µn
Z
00 (−)−1/2−µn Ξ∗1/2µn (ξn )Ξ1/2,µn (ξn )dξn
X
= C1/2,−µn ,1/2,µn
µn
X (−)1/2+µn Z
= −1/2−µn Ξ∗1/2,µn (ξn )Ξ1/2,µn (ξn )dξn
d (−)
µn 1
2
√
= 2, (C.29)
A similar procedure is followed for the integral over the proton spin functions. Thus, with
Id = 1
Z n
1 o
ΞI (ξnp ) ⊗ ΞI (ξnp ) dξnp = 1. (C.30)
Iˆd d d 00
The final integral in Eq.(C.27) can be worked out using Eq.(6), p.62, and Eq.(1), p.248, of
[93]
133
Z n Z
1 o 1 X 00
ỸL (R̂) ⊗ ỸL (R̂) dΩR = CL,−M LM ỸL,−M (R̂)ỸLM (R̂)dΩR
L̂ 00 L̂ M
Z
1 X 00 −L+M ∗ (R̂)Ỹ
= CL,−M LM (−) ỸLM LM (R̂)dΩR
L̂ M
1 X (−)L+M
= (−)−L+M = 1. (C.31)
L̂ M L̂
Therefore, introducing Eq.(C.31) into Eq.(C.27) and joining the rhs of Eq.(C.16) we obtain:
1 ~2
2
∂ L(L + 1) JT MT
2
− 2
+ Ed χLJ (R)
R 2µ ∂R R p
( )∗
X X Z n o
= − 8 ΞI (ξnp ) ⊗ Ỹ` (r̂) ⊗ ỸL (R̂) ⊗ ΞIt (ξt )
d jp Jp
`0 L0 Jp0 ` JT MT
× UN A Rp,n , Rp,n + 2s
( )
n o
× ΞI (ξnp ) ⊗ Ỹ`0 (r\
± 2s) ⊗ ỸL0 (R
\ + s) ⊗ ΞIt (ξt )
d jp Jp0 J T MT
J M
χ T0 0 T (|R + s|)
L Jp
× φ` (r)Vnp (r)φ`0 (|r ± 2s|) dsdrdΩR dξt dξnp (C.32)
|R + s|
We now concentrate on the tensor couplings in the rhs of Eq.(C.32). First, we introduce
134
( )∗
n o
ΞI (ξnp ) ⊗ Ỹ0 (r̂) ⊗ ỸL (R̂) ⊗ ΞIt (ξt )
d jp Jp JT MT
( )
n o
× ΞI (ξnp ) ⊗ Ỹ0 (r\
± 2s) ⊗ ỸL0 (R + s)
\ ⊗ ΞIt (ξt )
d jp Jp0 JT MT
∗
1 n o
= ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞIt (ξt )
4π d Jp JT MT
n o
× ΞI (ξnp ) ⊗ ỸL0 (R + s) 0 ⊗ ΞIt (ξt )
\
d Jp JT MT
1 X X J M J M
n o∗
= CJ TM TI µ C T0 0T 0 ΞI (ξnp ) ⊗ ỸL (R̂)
4π p p t t Jp Mp It µt d Jp Mp
Mp Mp0 µt µt 0
n o Z
× ΞI (ξnp ) ⊗ ỸL0 (R + s) 0 0
\ Ξ∗I µ (ξt )ΞI µ0 (ξt )dξt
d Jp Mp t t t t
1 X X J M J M
n o∗
= CJ TM TI µ C T0 0T ΞI (ξnp ) ⊗ ỸL (R̂)
4π p p t t Jp Mp It µt d Jp Mp
0 µ
M p Mp t
n o
× ΞI (ξnp ) ⊗ ỸL0 (R
\ + s) 0 0 .
d Jp Mp
(C.33)
1 ~2
2
∂ L(L + 1) JT M T
2
− 2
+ Ed χLJ (R)
R 2µ ∂R R p
Z
1 X
= − 8 φ0 (r)Vnp (r)UN A Rp,n , Rp,n + 2s
4π 0 L0 Jp
X X J M J M
n o∗
× CJ TM TI µ C T0 0T ΞI (ξnp ) ⊗ ỸL (R̂) (C.34)
p p t t Jp Mp It µt d Jp Mp
Mp Mp0 µt
J M
χ T0 0 T (|R + s|)
n o L Jp
× ΞI (ξnp ) ⊗ ỸL0 (R
\ + s) 0 0 φ0 (|r ± 2s|) dsdrdΩR dξnp .
d J p Mp |R + s|
135
Next we couple the two tensors together up to zero angular momentum using, Eq.(6), p.62,
n o∗ n o
ΞI (ξnp ) ⊗ ỸL (R̂) ΞI (ξnp ) ⊗ ỸL0 (R + s) 0 0
\
d J p Mp d Jp Mp
n o n o
= (−)Jp −Mp ΞI (ξnp ) ⊗ ỸL (R̂) ΞI (ξnp ) ⊗ ỸL0 (R
\ + s) 0 0
d Jp ,−Mp d J p Mp
X SM
= (−)Jp −Mp C S 0
Jp ,−Mp Jp Mp
SMS
n o n o
× ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞI (ξnp ) ⊗ ỸL0 (R
\ + s) 0
d Jp d Jp 00
→ (−)Jp −Mp CJ00,−M J 0 M
p p p p
n o n o
× ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞI (ξnp ) ⊗ ỸL0 (R + s) 0
\
d Jp d Jp 00
δJ J 0 δM M 0
J p −M p Jp +M p p p p p
= (−) (−)
ˆ
Jp
n o n o
× ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞI (ξnp ) ⊗ ỸL0 (R
\ + s) 0 (C.35)
d Jp d Jp 00
δJ J 0 δM M 0 n o
p p p p
o n
= (−)2J p ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞI (ξnp ) ⊗ ỸL0 (R
\ + s) 0 .
Jˆp d Jp d Jp 00
We replace Eq.(C.35) in Eq.(C.34) and sum over Jp0 and Mp0 to obtain:
136
1 ~2
2
∂ L(L + 1) J T MT
2
− 2
+ Ed χLJ (R)
R 2µ ∂R R p
Z
1 X
= − 8 φ0 (r)Vnp (r)UN A Rp,n , Rp,n + 2s
4π 0
L
XX J M J M 1
× CJ TM TI µ CJ TM TI µ (−)2Jp
Mp µt
p p t t p p t t Jˆp
n o n o
× ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞI (ξnp ) ⊗ ỸL0 (R + s)
\
d Jp d Jp 00
J M
χ T0 T (|R + s|)
L Jp
× φ0 (|r ± 2s|) dsdrdΩR dξnp . (C.36)
|R + s|
1 ~2
2
∂ L(L + 1) J T MT
2
− 2
+ Ed χLJ (R)
R 2µ ∂R R p
Z
1 X
= − 8 φ0 (r)Vnp (r)UN A Rp,n , Rp,n + 2s
4π 0
L
n o
2J 1 o n
×(−) p ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞI (ξnp ) ⊗ ỸL0 (R
\ + s)
Jˆp d Jp d Jp 00
J M
χ T0 T (|R + s|)
L Jp
× φ0 (|r ± 2s|) dsdrdΩR dξnp . (C.37)
|R + s|
We now use Eq.(11), p.70, and Eq.(4), p.358, of [93] to couple the spin functions and the
137
n o n o
ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ΞI (ξnp ) ⊗ ỸL0 (R + s)
\
d Jp d Jp 00
Id L Jp
nn o n o o
ˆ2
= Jp Id L0 Jp ΞI (ξnp ) ⊗ ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ỸL0 (R
\ + s)
d d 0 0 00
0 0 0
Jˆp nn o n o o
= δLL0 ΞI (ξnp ) ⊗ ΞI (ξnp ) ⊗ ỸL (R̂) ⊗ ỸL0 (R + s)
\ . (C.38)
Iˆd L̂ d d 0 0 00
1 ~2
2
∂ L(L + 1) JT M T
2
− 2
+ Ed χLJ (R)
R 2µ ∂R R p
Z
1 X
= − 8 φ0 (r)Vnp (r)UN A Rp,n , Rp,n + 2s
4π 0
L
1 n o n o
×(−)2Jp ΞI (ξnp ) ⊗ ΞI (ξnp ) ỸL (R̂) ⊗ ỸL (R
\ + s)
Iˆd L̂ d d 00 00
JT MT
χLJ (|R + s|)
p
× φ0 (|r ± 2s|) dsdrdΩR dξnp
|R + s|
Z
1 X
= − 8 φ0 (r)Vnp (r)UN A Rp,n , Rp,n + 2s (C.39)
4π 0
L
J M
T T (|R + s|)
1n o χLJ
p
× ỸL (R̂) ⊗ ỸL (R + s)
\ φ0 (|r ± 2s|) dsdrdΩR .
L̂ 00 |R + s|
Bringing the 1/R term from the lhs over to the rhs gives us:
138
2 2
~ ∂ L(L + 1) JT MT
2
− 2
+ Ed χLJ (R)
2µ ∂R R p
Z
R 8 n o
= − φ0 (r)Vnp (r)UN A Rp,n , Rp,n + 2s ỸL (R̂) ⊗ ỸL (R
\ + s)
L̂ 4π 00
J T MT
χLJ (|R + s|)
p
× φ0 (|r ± 2s|) dsdrdΩR
|R + s|
Z
R 8 X 00
= − CL,−M LM φ0 (r)Vnp (r)UN A Rp,n , Rp,n + 2s ỸL,−M (R̂)ỸLM (R
\ + s)
L̂ 4π
M
J T MT
χLJ (|R + s|)
p
× φ0 (|r ± 2s|) dsdrdΩR
|R + s|
R 8 X (−)L−M
Z
= − φ0 (r)Vnp (r)UN A Rp,n , Rp,n + 2s ỸL,−M (R̂)ỸLM (R
\ + s)
L̂ 4π M L̂
J T MT
χLJ (|R + s|)
p
× φ0 (|r ± 2s|) dsdrdΩR . (C.40)
|R + s|
Since the integrand is coupled to zero angular momentum, it is spherically symmetric, which
means that it is invariant under rotations of the three vectors R, r, and s. Thus, we can
Integration over dΩR yields a factor of 4π for all other choices for the direction of R. Since
we are fixing r to be in the xz-plane, we get a factor of 2π from each vector to take care of
rotations around the z-axis. Thus, we need to multiply the integral by (4π) ∗ (2π) = 8π 2 .
integral of Eq.(C.40) and the phase iL of the spherical harmonics, we arrive at the simplified
139
2 2
~ ∂ L(L + 1) J T MT
2
− 2
+ Ed χLJ (R)
2µ ∂R R p
Z
R 8 L̂ L 2L 1 2
= − √ (−) i 8π φ0 (r)Vnp (r)UN A Rp,n , Rp,n + 2s YL0 (R
\ + s)
L̂ 4π 4π L̂
JT MT
χLJ (|R + s|)
p
× φ0 (|r ± 2s|) dsr2 dr sin θr dθr
|R + s|
√ Z
8R π
= − φ0 (r)Vnp (r)UN A Rp,n , Rp,n + 2s YL0 (R
\ + s)
L̂
JT MT
χLJ (|R + s|)
p
× φ0 (|r ± 2s|) dsr2 dr sin θr dθr (C.41)
|R + s|
the nucleon optical potential for either the proton or neutron. Making the replacement
UN A → UnA + UpA gives us the nonlocal adiabatic potential, and the resulting partial wave
equation for the deuteron scattering state when using nonlocal potentials within the ADWA.
140
Appendix D
Here we will derive the explicit form for the T-matrix for (d, p) transfer reactions. This
equation was given in the post form in Eq.(2.21), and is repeated here without the remnant
term
µp M µ M
Tµ M µp M (kf , ki ) = hΨk B |Vnp |ΨkA d i. (D.1)
A d B f i
We need to define the explicit partial wave for all wave functions in Eq.(D.1). We begin by
defining the wave function for relative motion between d + A, which is given by:
X X n o
Ψ`i ji = ΞIp (ξp ) ⊗ Φji (rnp , ξn ) ⊗ ỸLi (R̂dA )
Jd JP M P
Li JP MP i i
i i
χLi J (RdA )
Pi
× ΞI µ (ξA )
A A RdA
( )
X X n o
= ΞIp (ξp ) ⊗ Ỹ`i (r̂np ) ⊗ ΞIn (ξn ) ⊗ ỸLi (R̂dA )
ji Jd
Li JP MP J P MP
i i i i
χLi Jp (RdA )
× ΞI µ (ξA )φji (rnp ) i . (D.2)
A A RdA
As in Appendix C, ΞIp (ξp ), ΞIn (ξn ) and ΞI (ξA ) are the spin functions for the proton,
A
141
neutron, and target, respectively, each with projections µp , µn , and µA . Ỹ`i is the spherical
harmonics for the relative motion between the neutron and proton in the deuteron, and ỸLi
is the spherical harmonic for the relative motion between the deuteron and the target. As
in Appendix C, we are defining our tensors with the built in factor of iL so that Ỹ`i = i`i Y`i
with Y`i defined on p.133, Eq.(1), of [93]. As a reminder, φji (rnp ) is the radial wave function
for the bound state, and ji results from coupling the orbital motion of the deuteron bound
state with the spin of the neutron. χLi Jp (RdA ) is the radial wave function for the deuteron
i
scattering state, and Jpi results from coupling the spin of the deuteron, Jd = 1 to the orbital
The incoming distored wave should depend only on the projections of the projectile and
target, and reduce to a plane wave in the limit of zero potential. Therefore, we multiply
J M
Eq.(D.2) by the incoming coefficient 4π Li iσLi P ∗ P i Pi
ki i e Mi0 ỸLi M 0 (k̂i )CJ M L M 0 giving us:
id d i i
Now, we make the following substitutions using Eq.(13), p.132, and Eq.(10), p.245, of [93]
0
ỸL∗ M 0 (k̂i ) = (−)Li +Mi ỸL ,−M 0 (k̂i )
i i i i
JP M P 0 JˆP J M
C i i 0 = (−)Li +Mi i C d d 0 , (D.4)
Jd Md Li Mi Jˆd Li ,−Mi JPi MPi
142
0
and insert (−)2(Li +Mi ) = 1, Eq.(D.3) becomes:
The partial wave for the incoming distorted wave is written as:
143
n o
Ψ` j = ΞI (ξA ) ⊗ Φjf (rnA , ξn )
f f A J B MB
χL J (RpB )
X X n o f Pf
× ΞIp (ξp ) ⊗ ỸL (R̂pB )
f J RpB
Lf JP MP P f M Pf
f f
( )
n o
= ΞI (ξA ) ⊗ Ỹ` (r̂nA ) ⊗ ΞIn (ξn ) φjf (rnA )
A f jf
JB MB
χL J (RpB )
X X n o f P f
× ΞIp (ξp ) ⊗ ỸL (R̂pB ) .
f JP MP RpB
Lf JP MP f f
f f
(D.7)
In this equation `f is the orbital angular momentum between the target and the bound
neutron, jf is the quantum number resulting from coupling `f to the spin of the target, IA .
The total angular momentum of the target is given by JB and results from coupling jf to
IA . The orbital angular momentum between the proton and the target is given by Lf , and
the total angular momentum of the projectile, Jpf results from coupling Lf to the spin of
the proton, Ip .
As we did for the entrance channel, we need to multiply the exit channel wave function
J P MP
L iσL P f f
by the outgoing coefficient: k4π i f e f M 0 Ỹ ∗ 0 (k̂f )C so that the remaining
f f Lf Mf Ip µp Lf M 0
f
quantum numbers are for the projections of the projectile and target in the exit channel:
144
( )
µp M 4π n o
|Ψ` j B i = ΞI (ξA ) ⊗ Ỹ` (r̂nA ) ⊗ ΞIn (ξn )
f f kf A f jf
J B MB
χL J (RpB )
X L iσL f P f
× φjf (rnA ) i fe f
RpB
Lf JP
f
J P MP n o
f f
YL∗ M 0 (k̂f )C
X
× ΞIp (ξp ) ⊗ ỸL (R̂pB ) .
Ip µp Lf M 0 f J
f f f Pf M P f
M 0 MP
f f
(D.8)
We follow the same steps as before, using Eq.(13), p.132, and Eq.(10), p.245, of [93]
L +Mf 0
YL∗ M 0 (k̂f ) = (−) f ỸL ,−M 0 (k̂f )
f f f f
ˆ
JP MP
f f Lf +Mf0 JPf Ip µp
C = (−) C , (D.9)
Ip µp Lf M 0 0
Iˆp Lf ,−Mf JPf MPf
f
2(Lf +Mf0 )
and use (−) = 1. This results in:
145
( )
µp M 4π L iσL
n o X
|Ψf B i = ΞI (ξA ) ⊗ Ỹ` (r̂nA ) ⊗ ΞIn (ξn ) φjf (rnA ) i fe f
kf A f jf
JB MB Lf JP
f
χL J (RpB ) Jˆ
f P f P f X Ip µp
× C Ỹ 0 (k̂ )
RpB Iˆp Lf ,−M 0 JP MP Lf ,−Mf f
f f
M 0 MP f
f f
n o
× ΞIp (ξp ) ⊗ ỸL (R̂pB )
f J Pf M Pf
( )
4π L iσL
n o X
= ΞI (ξA ) ⊗ Ỹ` (r̂nA ) ⊗ ΞIn (ξn ) φjf (rnA ) i fe f
kf A f jf
JB MB Lf JP
f
χL J (RpB ) Jˆ
f f P f P n o
× ỸL (k̂f ) ⊗ ΞIp (ξp ) ⊗ ỸL (R̂pB )
RpB Iˆp f f JP
f I µ
p p
( )
n o
(+)
= ΞI (ξA ) ⊗ Ỹ` (r̂nA ) ⊗ ΞIn (ξn ) φjf (rnA )χf (kf , RpB , ξp ).
A f jf
JB M B
(D.10)
( )∗
µp M
n o
(−)∗
hΨf B | = ΞI (ξA ) ⊗ Ỹ` (r̂nA ) ⊗ ΞIn (ξn ) φjf (rnA )χf (kf , RpB ),
A f jf
JB MB
(D.11)
where the outgoing distorted wave χ(−) (k, R) is the time reverse of χ(+) , so that χ(−) (k, R) =
χ(+) (−k, R)∗ . Therefore, to make this more explicit we use Eq.(2), p.141, of [93]
146
hχ(−) (k, R)| = χ(−)∗ (k, R)
∗
(+)
= χ (−k, R) ∗
= χ(+) (−k, R)
where k → −k gives a factor of (−)L from the spherical harmonics, as seen in Eq.(2), p.141,
and,
χL J (RpB )
(−)∗ 4π X −Lf iσLf ˆ f Pf
χf (kf , RpB , ξp ) = i e JP
kf Iˆp Lf JP
f RpB
f
n o
× ỸL (k̂f ) ⊗ ΞIp (ξp ) ⊗ ỸL (R̂pB ) . (D.14)
f f JP
f Ip µp
147
" ! #
~2 ∂2 L (L + 1)
− 2
− i i2 + U ad + V1L
SO
J + VC (RdA ) − E χLi J (RdA ) = 0
2µi ∂RdA RdA i Pi Pi
(D.15)
" ! #
~2 ∂2 Lf (Lf + 1)
− 2
− 2
+ U pB + VISO
p Lf JP
+ VC (RpB ) − E χL J (RpB ) = 0
2µf ∂RpB RpB f f Pf
(D.16)
with U pB being a nucleon optical potential. In these equations µi and µf are the reduced
mass in the initial and final states, not to be confused with spin projections µA , µp , and µn .
µf
fµ M µp M (kf , ki ) = − T̃ (k , ki )
A d B 2π~2 rµA Md µp MB f
µf vf
= − T (k , k )
2π~ 2 vi µA Md µp MB f i
v
u ~k
u f
µf u u µf µp M µ M
= − 2 t ~k hΨf B |Vnp |Ψi A d i. (D.17)
2π~ i
µi
The differential cross section is obtained, by averaging the mod of the scattering amplitude
148
squared over initial states, and summing over final m-states:
dσ 1 X 2
= f (k
µA Md µp MB f i , k )
dΩ Jˆd2 JˆA
2
µA Md µp MB
kf µi µf 1 X µp MB µA Md 2
= = hΨ f |V np |Ψi i . (D.18)
ki 4π 2 ~4 Jˆ2 Jˆ2
d A µA Md MB µp
µp MB µ M
We now put Eq.(D.5) and Eq.(D.11) into hΨf |Vnp |Ψi A d i:
Z ( )∗
µp M µ M
n o
hΨf B |Vnp |Ψi A d i = ΞI (ξA ) ⊗ Ỹ` (r̂nA ) ⊗ ΞIn (ξn ) φjf (rnA )
A f jf
JB M B
χL J (RpB )
4π X −Lf iσLf ˆ f Pf
× i e JP
kf Iˆp Lf JP
f RpB
f
n o
× ỸL (k̂f ) ⊗ ΞIp (ξp ) ⊗ ỸL (R̂pB ) (D.19)
f f JP
f Ip µp
Breaking the coupling between the target and the final bound state, and grouping the
149
µp MB µ M XZ
hΨf |Vnp |Ψi A d i = Ξ∗I µ0 (ξA )ΞI µ (ξA )dξA
A A A A
µ0A
Z X
J M
n o∗
× C B0B Ỹ` (r̂nA ) ⊗ ΞIn (ξn ) φjf (rnA )
IA µA jf mf f jf mf
mf
χL J (RpB )
4π X −Lf iσLf ˆ f Pf
× i e JP
kf Iˆp Lf JP
f RpB
f
n o
× ỸL (k̂f ) ⊗ ΞIp (ξp ) ⊗ ỸL (R̂pB ) (D.20)
f f JP
f Ip µp
The integral in the first line gives δµ0 µ . Performing the sum over µ0A provides:
A A
150
µp M µ M
Z X
J M
n o∗
hΨf B |Vnp |Ψi A d i = CI Bµ Bj m Ỹ` (r̂nA ) ⊗ ΞIn (ξn ) φjf (rnA )
A A f f f jf mf
m f
χL J (RpB )
4π X −Lf iσLf ˆ f Pf
× i e JP
kf Iˆp Lf JP
f RpB
f
n o
× ỸL (k̂f ) ⊗ ΞIp (ξp ) ⊗ ỸL (R̂pB ) (D.21)
f f JP
f Ip µp
n o
ỸL (k̂f ) ⊗ ΞIp (ξp ) ⊗ ỸL (R̂pB )
f f JP
f Ip µp
( )
n o
× ỸLi (k̂i ) ⊗ ΞIp (ξp ) ⊗ Ỹ`i (r̂np ) ⊗ ΞIn (ξn ) ⊗ ỸLi (R̂dA )
ji Jd
JP
i J d Md
X QMQ n o
= CI µ J M ỸL (k̂f ) ⊗ ΞIp (ξp ) ⊗ ỸL (R̂pB ) (D.22)
p p d d f f JP
QM Q
f Ip
( )
n o
⊗ ỸLi (k̂i ) ⊗ ΞIp (ξp ) ⊗ Ỹ`i (r̂np ) ⊗ ΞIn (ξn ) ⊗ ỸLi (R̂dA ) ,
ji J
d J
Pi J
d QMQ
151
to obtain:
(4π)2
Z X X
µp M µ M J M QMQ
hΨf B |Vnp |Ψi A d i = CI Bµ Bj m CI µ J M
ki kf Jˆd Iˆp mf QM
A A f f p p d d
Q
n o∗ X X
× Ỹ` (r̂nA ) ⊗ ΞIn (ξn ) φjf (rnA )
f jf mf
Li JP Lf JP
i f
χL J (RpB ) χLi Jp (RdA )
−L iσL f P f iσ
× i f e f JˆP V (rnp )φji (rnp )iLi e Li JˆPi i
f RpB RdA
n o
× ỸL (k̂f ) ⊗ ΞIp (ξp ) ⊗ ỸL (R̂pB )
f
f JP
f Ip
( )
n o
⊗ ỸL (k̂i ) ⊗ ΞIp (ξp ) ⊗ Ỹ`i (r̂np ) ⊗ ΞIn (ξn ) ⊗ ỸLi (R̂dA )
i
ji Jd
JP
i Jd QM
Q
× dRpB drnA dξn dξp .
(D.23)
µp MB µ M X X J M QMQ
hΨf |Vnp |Ψi A d i = CI Bµ Bj m CI µ J M TQM m , (D.24)
A A f f p p d d Q f
mf QM
Q
so that
152
µp MB µ M
|Vnp |Ψi A d i|2
X
|hΨf
µA Md MB µp
X X XJB M B JB M B QMQ Q0 MQ0
= CI µ j m C C C
A A f f IA µA jf m0f Ip µp Jd Md Ip µp Jd Md
µA Md MB µp mf QMQ m0 Q0 M 0
f Q
∗
× TQM m TQ0 M 0 m0 . (D.25)
Q f Q f
Now we consider the first pair of Clebsch-Gordans and use Eq.(10), p.245, of [93],
ˆB 2 jf ,−mf jf ,−m0f
!
J M J M I −µ J
CI Bµ Bj m C B B 0 = (−) A A CI µ J ,−M CI µ J ,−M .
A A f f A µA jf mf
I ĵf A A B B A A B B
(D.26)
µp MB µ M
|Vnp |Ψi A d i|2
X
|hΨf
µA Md MB µp
JˆB
2 X X X jf ,−mf jf ,−m0f
= CI µ J ,−M CI µ J ,−M
ĵf2 m QM A A B B A A B B
f Q m0f Q0 MQ
0 µA MB
X QMQ Q0 MQ
0
× ∗
CI µ J M CI µ J M TQM m TQ
p p d d p p d d Q f 0 M 0 m0
µp Md Q f
JˆB
2 X X
∗
= 2
δm m0 δQQ0 δM M 0 TQM m TQ 0 M 0 m0
ĵf m QM f f Q Q Q f Q f
f 0 0 0
Q m f Q MQ
JˆB
2 X
∗
= TQM m TQM , (D.27)
2
ĵf m QM Q f Q mf
f Q
153
and then the differential cross section is
dσ kf µi µf 1 X µp MB µ M 2
= hΨf |Vnp |Ψi A d i
dΩ ki 4π 2 ~4 Jˆ2 Jˆ2
d A µA Md MB µp
kf µi µf JˆB2 X
∗
= TQM m TQM . (D.28)
ki 4π 2 ~4 Jˆ2 Jˆ2 ĵ 2 Q f Q mf
d A f mf QMQ
(D.29)
Our strategy is to couple the spherical harmonics with the argument k̂ together so we can pull
them out of the integral. We want to couple the spherical harmonics with the arguments
r̂ and R̂ together up to zero angular momentum so we can use symmetry to reduce the
dimensionality of the angular integral. Also, we want to couple the spinors with common
154
arguments up to zero angular momentum so we can integrate them out. This is detailed in
We can group the k̂ spherical harmonics together right away. Let us introduce the
definitions
n o
AJ = ΞIp (ξp ) ⊗ ỸL (R̂pB )
Pf f J Pf
( )
n o
BJ = ΞIp (ξp ) ⊗ Ỹ`i (r̂np ) ⊗ ΞIn (ξn ) ⊗ ỸLi (R̂dA ) (D.30)
Pi ji Jd JP
i
( )
ỸL (k̂f ) ⊗ AJ ⊗ ỸLi (k̂i ) ⊗ BJ = |Lf JP (Ip )Li JPi (Jd )QMQ i
f P f P i f
Ip Jd QMQ
X
= |Lf Li (g)JP JPi (h)QMQ ihLf Li (g)JP JPi (h)QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f f
gh
X
= hLf Li (g)JP JPi (h)QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f
gh
n o
× ỸL (k̂f ) ⊗ ỸLi (k̂i ) ⊗ AJ ⊗ BJ , (D.31)
f g Pf Pi
h QMQ
155
j1 j2 j12
0 0
hj1 j2 (j12 )j3 j4 (j34 )jm|j1 j3 (j13 )j2 j4 (j24 )j m i = δjj 0 δmm0 ĵ12 ĵ13 ĵ24 ĵ34 j3 j4 j34 .
j
13 j24 j
(D.32)
156
AJ ⊗ BJ
Pf Pi
h
n o
= ΞIp (ξp ) ⊗ ỸL (R̂pB )
f JP
f
( )
n o
⊗ ΞIp (ξp ) ⊗ Ỹ`i (r̂np ) ⊗ ΞIn (ξn ) ⊗ ỸLi (R̂dA )
ji Jd
JP
i h
= |Ip Lf (JP )Jd Li (JPi )hmh i
f
|Ip Jd (g 0 )Lf Li (h0 )hmh ihIp Jd (g 0 )Lf Li (h0 )hmh |Ip Lf (JP )Jd Li (JPi )hmh i
X
=
f
g 0 h0
hIp Jd (g 0 )Lf Li (h0 )hmh |Ip Lf (JP )Jd Li (JPi )hmh i
X
=
f
g 0 h0
( )
n o
× Ξ (ξp ) ⊗ ΞIp (ξp ) ⊗ Ỹ`i (r̂np ) ⊗ ΞIn (ξn )
Ip ji J
d g0
n o o
⊗ ỸL (R̂pB ) ⊗ ỸLi (R̂dA ) 0 . (D.34)
f h h
157
(4π)2 X X Li −Lf i(σLi +σLf ) ˆ ˆ
TQM m = i e JPi JP
Q f ki kf Jˆd Iˆp L J L J f
i P f P i f
Z φj (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA )
f f Pf i
×
RpB RdA
X
× hLf Li (g)JP JPi (h)QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f
gh
hIp Jd (g 0 )Lf Li (h0 )hmh |Ip Lf (JP )Jd Li (JPi )hmh i
X
×
f
g 0 h0
( )
n o n o
× Ỹ (k̂ ) ⊗ ỸLi (k̂i ) ⊗ Ξ (ξp ) ⊗ ΞIp (ξp ) ⊗ Ỹ`i (r̂np ) ⊗ ΞIn (ξn )
Lf f g Ip ji J
d g0
n o o n o∗
⊗ ỸL (R̂pB ) ⊗ ỸLi (R̂dA ) 0 Ỹ` (r̂nA ) ⊗ ΞIn µn (ξn )
f h h QM f jf mf
Q
× dRpB drnA dξn dξp . (D.35)
The following tensor in Eq.(D.35) can be simplified using, Eq.(27), p.64, and Eq.(8), p.70,
of [93]
( )
n o
ΞIp (ξp ) ⊗ ΞIp (ξp ) ⊗ Ỹ`i (r̂np ) ⊗ ΞIn (ξn )
ji J
d g0
( )
0 n o
= (−)I p +Jd −g ΞIp (ξp ) ⊗ Ỹ`i (r̂np ) ⊗ ΞIn (ξn ) ⊗ ΞIp (ξp )
ji Jd g0
0 0X
Ip ji Jd
= (−)Ip +Jd −g (−)Jd +Ip +g Jˆd q̂
q g0 I
q
p
n o n o
× Ỹ`i (r̂np ) ⊗ ΞIn (ξn ) ⊗ ΞIp (ξp ) ⊗ ΞIp (ξp ) . (D.36)
ji q g0
Since the spin functions must be coupled to zero angular momentum, this implies that q = 0
158
and g 0 = ji . Therefore, Eq.(D.36) simplifies to, with (−)2Jd = 1, and using Eq.(1), p.299, of
[93]
Ip ji Jd
n o n o
2Ip
(−) Jd ˆ Ỹ`i (r̂np ) ⊗ ΞIn (ξn ) ⊗ ΞIp (ξp ) ⊗ ΞIp (ξp )
j I ji 0 j
p 0 i
i
Jˆd n
o n o
2Ip
= (−) (−) I p +j i +Jd Ỹ`i (r̂np ) ⊗ ΞIn (ξn ) ⊗ ΞIp (ξp ) ⊗ ΞIp (ξp ) .
Iˆp ĵi ji 0 j
i
(D.37)
We now take the last two lines of Eq.(D.38), break all the couplings between the pairs, and
159
( (
n o n o n o
× ỸL (k̂f ) ⊗ ỸLi (k̂i ) ⊗ Ỹ`i (r̂np ) ⊗ ΞIn (ξn ) ⊗ ΞIp (ξp ) ⊗ ΞIp (ξp )
f g ji 0 j
i
n o o
⊗ ỸL (R̂pB ) ⊗ ỸLi (R̂dA ) 0
f h h QM
Q
X QMQ n o X hm
n o
= Cgm hm ỸL (k̂f ) ⊗ ỸLi (k̂i ) C h 0 Ỹ`i (r̂np ) ⊗ ΞIn (ξn )
g h f gmg ji m i h m 0 ji mi
mg mh mi mh0 h
n o n o
× ΞIp (ξp ) ⊗ ΞIp (ξp ) ỸL (R̂pB ) ⊗ ỸLi (R̂dA ) 0 . (D.39)
00 f h mh0
n o n o∗
We now consider ỸL (R̂pB ) ⊗ ỸLi (R̂dA ) 0 Ỹ` (r̂nA ) ⊗ ΞIn µn (ξn ) and use Eq.(23),
f h mh0 f jf mf
p.64, of [93] to obtain:
160
j −mf
n o n o
(−) f Ỹ` (r̂nA ) ⊗ ΞIn (ξn ) Ỹ`i (r̂np ) ⊗ ΞIn (ξn )
f jf ,−mf ji mi
( )
jf −mf X KM
n o n o
= (−) Cj ,−m j m Ỹ` (r̂nA ) ⊗ ΞIn (ξn ) ⊗ Ỹ`i (r̂np ) ⊗ ΞIn (ξn )
f f i i f jf ji
KM KM
j −mf X KM
= (−) f Cj ,−m j m |`f In (jf )`i In (ji )KM i
f f i i
KM
j −mf X KM
h`f `i (g 00 )In In (h00 )KM |`f In (jf )`i In (ji )KM i
X
= (−) f Cj ,−m j m
f f i i
KM g 00 h00
× |`f `i (g 00 )In In (h00 )KM i
j −mf X KM
h`f `i (g 00 )In In (h00 )KM |`f In (jf )`i In (ji )KM i
X
= (−) f Cj ,−m j m
f f i i
KM g 00 h00
n o
× Ỹ` (r̂nA ) ⊗ Ỹ`i (r̂np ) 00 ⊗ ΞIn (ξn ) ⊗ ΞIn (ξn ) h00 . (D.41)
f g KM
Since the spins must be coupled up to zero, we see that h00 = 0 and g 00 = K. Imposing this
condition in Eq.(D.41):
161
3I +j +J +j −mf X
(4π)2 (−) p i d f X L −L i(σL +σL )
f Jˆ Jˆ
TQM m = i i fe i
Pi P f
Q f Iˆp2 ki kf ĵi Li JP Lf JP
i f
Z φj (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA )
f f Pf i
×
RpB RdA
X
× hLf Li (g)JP JPi (h)QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f
gh
hIp Jd (ji )Lf Li (h0 )hmh |Ip Lf (JP )Jd Li (JPi )hmh i
X
×
f
h0
QMQ hm
CjKM
X X X
× Cgm hm C h 0
g h j m
i i h mh0 KM f ,−mf ji mi
mg mh mi m 0
h
n o
× h`f `i (K)In In (0)KM |`f In (jf )`i In (ji )KM i Ỹ` (r̂nA ) ⊗ Ỹ`i (r̂np )
f KM
n o n o
× ΞIn (ξn ) ⊗ ΞIn (ξn ) 00 ỸL (k̂f ) ⊗ ỸLi (k̂i ) ΞIp (ξp ) ⊗ ΞIp (ξp )
f gmg 00
n o
× ỸL (R̂pB ) ⊗ ỸLi (R̂dA ) 0 dRpB drnA dξn dξp . (D.42)
f h mh0
Now we couple the r̂ and R̂ spherical harmonics up to zero using Eq.(1), p.248, of [93]
n o n o
ỸL (R̂pB ) ⊗ ỸLi (R̂dA ) 0 Ỹ` (r̂nA ) ⊗ Ỹ`i (r̂np )
f hm 0 f KM
h
X SM nn o n o o
= C 0 S ỸL (R̂pB ) ⊗ ỸLi (R̂dA ) 0 ⊗ Ỹ` (r̂nA ) ⊗ Ỹ`i (r̂np )
h m 0 KM f h f K SMS
SMS h
nn o n o o
→ 00
Ch0 m KM ỸL (R̂pB ) ⊗ ỸLi (R̂dA ) 0 ⊗ Ỹ` (r̂nA ) ⊗ Ỹ`i (r̂np )
h0 f h f K 00
δ δ
h0 −mh0 h0 K mh0 ,−M
nn o n o o
= (−) ỸL (R̂pB ) ⊗ ỸLi (R̂dA ) 0 ⊗ Ỹ` (r̂nA ) ⊗ Ỹ`i (r̂np )
ĥ0 f h f K 00
δ δ X (−)K+MK n
h0 −mh0 h0 K mh0 ,−M
o
= (−) ỸL (R̂pB ) ⊗ ỸLi (R̂dA )
ĥ0 MK
K̂ f K,−MK
n o
× Ỹ` (r̂nA ) ⊗ Ỹ`i (r̂np ) , (D.43)
f KMK
162
to arrive at:
3I +j +J +j −mf X
(4π)2 (−) p i d f X L −L i(σL +σL )
f Jˆ Jˆ
TQM m = i i fe i
Pi Pf
Q f Iˆp2 ki kf ĵi Li JP Lf JP
i f
Z φj (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA )
f f Pf i
×
RpB RdA
X
× hLf Li (g)JP JPi (h)QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f
gh
hIp Jd (ji )Lf Li (h0 )hmh |Ip Lf (JP )Jd Li (JPi )hmh i ΞIn (ξn ) ⊗ ΞIn (ξn ) 00
X
×
f
h0
QMQ hm
n o
CjKM
X X X
× Cgm hm C h 0 ỸL (k̂f ) ⊗ ỸLi (k̂i )
g h ji mi h mh0 f ,−mf ji mi f gmg
mg mh mi m 0 KM
h
n o
× h`f `i (K)In In (0)KM |`f In (jf )`i In (ji )KM i ΞIp (ξp ) ⊗ ΞIp (ξp )
00
δ 0 δ K+M
h0 −mh0 h K mh0 ,−M (−) K
X n o
× (−) ỸL (R̂pB ) ⊗ ỸLi (R̂dA )
ĥ0 MK
K̂ f K,−MK
n o
× Ỹ` (r̂nA ) ⊗ Ỹ`i (r̂np ) dRpB drnA dξn dξp . (D.44)
f KMK
Next we sum over h0 and mh0 , so that h0 = K, and mh0 = −M . Then Eq.(D.44) becomes:
163
3I +j +J +j −mf X
(4π)2 (−) p i d f X L −L i(σL +σL )
f Jˆ Jˆ
TQM m = i i fe i
Pi Pf
Q f Iˆp2 ki kf ĵi Li JP Lf JP
i f
Z φj (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA )
f f Pf i
×
RpB RdA
X
× hLf Li (g)JP JPi (h)QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f
gh
× hIp Jd (ji )Lf Li (K)hmh |Ip Lf (JP )Jd Li (JPi )hmh i ΞIn (ξn ) ⊗ ΞIn (ξn ) 00
f
X QMQ X hm n o
CjKM
X
× Cgm hm Cj mhK−M ,−m j m Ỹ L ( k̂f ) ⊗ Ỹ (k̂
Li i )
g h i i f f i i f gmg
mg mh mi KM
n o
× h`f `i (K)In In (0)KM |`f In (jf )`i In (ji )KM i ΞIp (ξp ) ⊗ ΞIp (ξp )
00
(−) K+M X (−)K+M K n o
× ỸL (R̂pB ) ⊗ ỸLi (R̂dA )
K̂ K̂ f K,−MK
MK
n o
× Ỹ` (r̂nA ) ⊗ Ỹ`i (r̂np ) dRpB drnA dξn dξp . (D.45)
f KMK
The integrals over dξn and dξp give Iˆp Iˆn , so that:
164
3I +j +J +j −mf
(4π)2 (−) p i d f X X L −L i(σL +σL )
TQM m = Iˆn i i fe i f Jˆ Jˆ
Pi Pf
Q f Iˆp ki kf ĵi Li JP Lf JP
i f
Z φj (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA )
f f Pf i
×
RpB RdA
X 1 X
× hLf Li (g)JP JPi (h)QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
K̂ 2 f f
K gh
× hIp Jd (ji )Lf Li (K)hmh |Ip Lf (JP )Jd Li (JPi )hmh i
f
QMQ hm
(−)M Cj mhK−M CjKM
X X
× Cgm hm
g h i i f ,−mf ji mi
mg mh mi M
(D.46)
We next consider the sum over Clebsch-Gordan coefficients, and use Eq.(11), p.245, of [93]:
165
hm
(−)M Cj mhK,−M CjKM
X
i i f ,−mf ji mi
mi M
hm j +mf K̂ ji ,−mi
(−)M Cj mhK,−M (−) f
X
= Cj ,−m K,−M
i i ĵi f f
mi M
hm j +mf K̂ ĵ jf mf
(−)M Cj mhK,−M (−) f (−)K−M i Cj m K,−M
X
=
i i ĵi ĵf i i
mM i
j +mf +K K̂ X hmh jf mf
= (−) f Cj m K,−M Cj m K,−M
ĵf m M i i i i
i
j +mf +K K̂
= (−) f δhj δm m . (D.47)
ĵf f h f
3I +j +J +2j
(4π)2 (−) p i d f ˆ X X Li −Lf i(σLi +σLf ) ˆ ˆ
TQM m = In i e JPi JP
Q f Iˆp ki kf ĵi ĵf Li JP Lf JP
f
i f
Z φj (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA )
f f Pf i
×
RpB RdA
X (−)K X
× hLf Li (g)JP JPi (jf )QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
K̂ f f
K g
X QMQ
× hIp Jd (ji )Lf Li (K)jf mf |Ip Lf (JP )Jd Li (JPi )jf mf i Cgmg j m
f f f
mg
n o
× h`f `i (K)In In (0)KM |`f In (jf )`i In (ji )KM i ỸL (k̂f ) ⊗ ỸLi (k̂i )
f gmg
n o
(−)MK ỸL (R̂pB ) ⊗ ỸLi (R̂dA )
X
×
f K,−MK
MK
n o
× Ỹ` (r̂nA ) ⊗ Ỹ`i (r̂np ) dRpB drnA . (D.48)
f KMK
166
3I +j +J +2j
(4π)2 (−) p i d f ˆ X (−)K
TQM m = In
Q f Iˆp ki kf ĵi ĵf K
K̂
× h`f `i (K)In In (0)KM |`f In (jf )`i In (ji )KM i
X X L −L i(σL +σL )
× i i fe i f Jˆ Jˆ
P i Pf
Li JP Lf JP
i f
× hIp Jd (ji )Lf Li (K)jf mf |Ip Lf (JP )Jd Li (JPi )jf mf i
f
X
× hLf Li (g)JP JPi (jf )QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f
g
X QMQ n o
× Cgmg j m ỸL (k̂f ) ⊗ ỸLi (k̂i ) (D.49)
f f f gmg
mg
Z φj (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA ) X
f f Pf i
× (−)MK
RpB RdA
MK
n o n o
× ỸL (R̂pB ) ⊗ ỸLi (R̂dA ) Ỹ` (r̂nA ) ⊗ Ỹ`i (r̂np ) dRpB drnA .
f K,−M f KM
K K
gmg
n o
L +L X
ỸL (k̂f ) ⊗ ỸLi (k̂i ) = i f i C Y (k̂ )Y (k̂
f gm g Lf M̃f Li M̃i Lf M̃f f Li M̃i i
M̃i M̃f
L +L X gmg L̂
= i f i C YL M̃ (k̂f ) √ i δM̃ m
Lf M̃f Li 0 f f 4π f g
M̃f
L +L gmg L̂
= i f i CL m L 0 YL mg (k̂f ) √ i . (D.50)
f g i f 4π
167
3I +j +J +2j
(4π)3/2 (−) p i d f ˆ X (−)K
TQM m = In
Q f Iˆp ki kf ĵi ĵf K
K̂
× h`f `i (K)In In (0)KM |`f In (jf )`i In (ji )KM i
i(σL +σL )
i2Li e
X X
× i f L̂ Jˆ Jˆ
i Pi Pf
Li JP Lf JP
i f
× hIp Jd (ji )Lf Li (K)jf mf |Ip Lf (JP )Jd Li (JPi )jf mf i
f
X
× hLf Li (g)JP JPi (jf )QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f
g
X QMQ gmg
× Cgmg j m CL m L 0 YL mg (k̂f ) (D.51)
f f f g i f
mg
Z φj (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA ) X
f f Pf i
× (−)MK
RpB RdA
MK
n o n o
× ỸL (R̂pB ) ⊗ ỸLi (R̂dA ) Ỹ` (r̂nA ) ⊗ Ỹ`i (r̂np ) dRpB drnA .
f K,−M f KM
K K
168
3I +j +J +2j
(4π)3/2 (−) p i d f ˆ X (−)K
TQM m = In
Q f Iˆp ki kf ĵi ĵf K
K̂
× h`f `i (K)In In (0)KM |`f In (jf )`i In (ji )KM i
i(σL +σL )
i2Li e
X X
× i f L̂ Jˆ Jˆ
i Pi Pf
Li JP Lf JP
i f
× hIp Jd (ji )Lf Li (K)jf mf |Ip Lf (JP )Jd Li (JPi )jf mf i
f
X
× hLf Li (g)JP JPi (jf )QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f
g
X QMQ gmg
× Cgmg j m CL m L 0 YL mg (k̂f )
f f f g i f
mg
Z φj (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA )
f f Pf i
×
RpB RdA
L +L K,−MK
(−)MK
X X
× i f i CL M L YL M (R̂pB )YLi Mi (R̂dA )
f f i Mi f f
MK Mf Mi
X ` +` KM
× i f i C` m̃K` m̃ Y` m̃ (r̂nA )Y`i m̃i (r̂np )dRpB drnA .
f f i i f f
m̃f m̃i
(D.52)
RpB is another independent variable. We place R̂pB in the ẑ-direction. In that case, Mf = 0,
√
YL M (R̂pB ) = L̂f / 4π, and Mi = −MK . Eq.(D.52) is then simplified to:
f f
169
3I +j +J +2j
4π (−) p i d f ˆ X (−)K
TQM m = In
Q f Iˆp ki kf ĵi ĵf K
K̂
× h`f `i (K)In In (0)KM |`f In (jf )`i In (ji )KM i
X X 3L +L +` +` i(σL +σL )
× i i f f ie i f L̂ L̂ Jˆ Jˆ
i f P i Pf
Li JP Lf JP
i f
× hIp Jd (ji )Lf Li (K)jf mf |Ip Lf (JP )Jd Li (JPi )jf mf i
f
X
× hLf Li (g)JP JPi (jf )QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f
g
X QMQ gmg
× Cgmg j m CL m L 0 YL mg (k̂f )
f f f g i f
mg
Z φj (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA )
f f Pf i
×
RpB RdA
K,−M
(−)MK CL 0L ,−M
X
× K YLi ,−M (R̂dA ) (D.53)
f i K K
MK
X KM 2 dΩ 2
× C` m̃K` m̃ Y` m̃ (r̂nA )Y`i m̃i (r̂np )RpB RpB dRpB rnA sin θdrnA dθdφ.
f f i i f f
m̃f m̃i
Since we are fixing R̂pB in the ẑ-direction, integrating over dΩR results only in a factor of
pB
4π. We also fix the other vectors to be in the xz-plane, which means that the integral over
170
3I +j +J +2j
32π 3 (−) p i d f ˆ X (−)K
TQM m = In
Q f Iˆp ki kf ĵi ĵf K
K̂
× h`f `i (K)In In (0)KM |`f In (jf )`i In (ji )KM i
X X 3L +L +` +` i(σL +σL )
× i i f f ie i f L̂ ˆ ˆ
i L̂f JPi JPf
Li JP Lf JP
i f
× hIp Jd (ji )Lf Li (K)jf mf |Ip Lf (JP )Jd Li (JPi )jf mf i
f
X
× hLf Li (g)JP JPi (jf )QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f
g
X QMQ gmg
× Cgmg j m CL m L 0 YL mg (k̂f )
f f f g i f
mg
Z 2
RpB rnA
× φjf (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA )
f Pf i RdA
K,−MK
(−)MK CL 0L ,−M
X
× YLi ,−M (R̂dA )
f i K K
MK
X KM
× C` m̃K` m̃ Y` m̃ (r̂nA )Y`i m̃i (r̂np ) sin θdRpB drnA dθ. (D.54)
f f i i f f
m̃f m̃i
Eq.(D.54) is valid for a general `i and `f . Since we are interested in applying the formalism
√
to (d, p) we use `i = 0 deuteron. Therefore Y`i m̃i (r̂np ) = 1/ 4π, m̃i = 0, m̃f = MK , and
171
3I +j +J +2j +`
32π 3 (−) p i d f f ˆ
TQM m = √ In
Q f Iˆp 4πki kf ĵi ĵf `ˆf
× h`f 0(`f )In In (0)`f M |`f In (jf )0In (ji )`f M i
X X 3L +L +` +` i(σL +σL )
× i i f f ie i f L̂ ˆ ˆ
i L̂f JPi JPf
Li JP Lf JP
i f
× hIp Jd (ji )Lf Li (`f )jf mf |Ip Lf (JP )Jd Li (JPi )jf mf i
f
X
× hLf Li (g)JP JPi (jf )QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f
g
X QMQ gmg `f ,−MK
(−)MK CL 0L ,−M
X
× Cgmg j m CL m L 0 YL mg (k̂f )
f f f g i f f i K
mg MK
Z 2
RpB rnA
× φjf (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA )
f Pf i RdA
× YLi ,−M (R̂dA )Y` M (r̂nA ) sin θdRpB drnA dθ. (D.55)
K f K
We can further simplify Eq.(D.55) by using, Eq.(5), p.334, Eq.(1), p.357 Eq. (1), and Eq.(1),
p.299, of [93]:
172
h`f 0(`f )In In (0)`f M |`f In (jf )0In (In )`f M i
`f 0 `f
ˆ ˆ
= `f 0̂ĵf In In In 0
j I `
n
f f
jf In `f
ˆ ˆ
= `f 0̂ĵf In `f 0 `f
I I
n n 0
jf In `f
In +`f +`f +In 1
ˆ ˆ
= `f 0̂ĵf In (−)
`ˆf Iˆn
0 ` I
f n
I +` +` +I 1 j +I +` 1
= `ˆf 0̂ĵf Iˆn (−) n f f n (−) f n f
`ˆf Iˆn Iˆn `ˆf
3I +` +j ĵf
= (−) n f f . (D.56)
`ˆf Iˆn
2`
Since `f is an integer, (−) f = 1, and we obtain:
173
3I +3I +j +J +3j
32π 3 (−) p n i d f
TQM m = √
Q f Iˆp 4πki kf ĵi `ˆ2f
X X 3L +L +` +` i(σL +σL )
× i i f f ie i f L̂ L̂ Jˆ Jˆ
i f Pi P f
Li JP Lf JP
i f
× hIp Jd (ji )Lf Li (`f )jf mf |Ip Lf (JP )Jd Li (JPi )jf mf i
f
X
× hLf Li (g)JP JPi (jf )QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i
f f
g
X QMQ gmg `f ,−MK
(−)MK CL 0L ,−M (D.57)
X
× Cgmg j m CL m L 0 YL mg (k̂f )
f f f g i f f i K
mg MK
Z 2
RpB rnA
× φjf (rnA )χL J (RpB )V (rnp )φji (rnp )χLi Jp (RdA )
f Pf i RdA
× YLi ,−M (R̂dA )Y` M (r̂nA ) sin θdRpB drnA dθ.
K f K
Lf Li g
ˆ ˆ
hLf Li (g)JP JPi (jf )QMQ |Lf JP (Ip )Li JPi (Jd )QMQ i = ĝ ĵf Ip Jd JP JP jf
f f f i
I Jd Q
p
Ip Jd ji
ˆ ˆ ˆ
hIp Jd (ji )Lf Li (`f )jf mf |Ip Lf (JP )Jd Li (JPi )jf mf i = ĵi `f JP JPi Lf Li `f ,
f f
J J j
P f Pi f
(D.58)
Finally, inserting Eq.(D.58) into Eq.(D.57) gives us the form for TQM m which we imple-
Q f
ment in NLAT:
174
32π 3 ˆ
3I +3I +j +J +3j ĵf Jd X
TQM m = √ (−) p n i d f
Q f 4πki kf `ˆf L J
i Pi
3L +L +` +` i(σL +σL )
ˆ2 ˆ2
X
× i i f f ie i f L̂
i L̂f JPi JPf
Lf JP
f
Ip Jd ji Lf Li g
X QM
gmg
X
Q
× Lf Li `f ĝ JP JP jf Cgmg j m CL m L 0
f i f f f g i
g
m
g
J
P JPi jf I
p Jd Q
f
`f ,−MK
(−)MK CL 0L ,−M
X
× YL mg (k̂f )
f f i K
MK
Z
× φjf (rnA )χL J (RpB )V (rnp )φji (rnp )
f P f
2
RpB rnA
× χLi Jp (RdA ) Y Li ,−MK (R̂dA )Y`f MK (r̂nA ) sin θdRpB drnA dθ,
i RdA
(D.59)
The observable is the differential cross section, which as we saw in Eq.(D.28), is given by
dσ kf µi µf JˆB
2 X
∗
= TQM m TQM . (D.60)
dΩ 2 4 ˆ2 ˆ
ki 4π ~ J J ĵ2 2 Q f Q mf
d A f mf QMQ
175
Appendix E
To perform the calculations in this thesis, the code “Nonlocal Adiabatic Transfer” (NLAT)
was written to calculate (d, p) transfer reactions with the inclusion of nonlocality. In order
to ensure that the code works properly, multiple checks were performed, and are discussed
in the following sections. When making comparisons to local calculations, we will compare
First, we look at the local elastic scattering distribution. In Fig. E.1 we show this check for
the reaction 209 Pb(p, p)209 Pb at Ep = 50.0 MeV. The solid line is a local calculation using
NLAT, the dotted line is a nonlocal calculation, but with β = 0.05 fm so that it reduces to
the local calculation, and the dashed line is FRESCO. We used β = 0.05 fm rather than
β = 0 fm since we would have numerical problems with dividing by zero if we set β exactly
equal to zero. For these calculations, we used a step size of 0.01 fm, a maximum radius of
30 fm, and included partial waves up to L = 20. These calculations are converged in that a
smaller step size or more partial waves does not change the results of the calculation.
176
0
10
σ/σR
-1
10
NLAT Local
NLAT β=0.05
Fresco
-2
10 0 30 60 90 120 150 180
θc.m. [deg]
Next, we look at the nonlocal elastic scattering distribution. In Fig. E.2 we present
208 Pb(n, n)208 Pb at En = 14.5 MeV. The solid line is a nonlocal calculation with β = 0.85
fm using NLAT. The dashed line is the digitized results of the same calculation from the
paper of Perey and Buck [1]. The two calculations agree quite well, indicating that NLAT
calculates elastic scattering with a nonlocal potential properly. The calculations of Perey
and Buck were digitized from their paper, so any discrepancies between the results shown
177
5
10
NLAT
4 Perey-Buck
10
dσ/dΩ [mb/sr]
3
10
2
10
1
10
0
10 0 30 60 90 120 150 180
θc.m. [deg]
Figure E.2: Differential elastic scattering as a function of scattering angle. 208 Pb(n, n)208 Pb
at Ep = 14.5 MeV: The solid line is obtained a nonlocal calculation using NLAT, and the
dashed line is the nonlocal calculation published by Perey and Buck [1].
Bound States
Next, we examine the bound wave functions. In Fig. E.3 we show the n+48 Ca bound wave
function as well as the deuteron bound wave function. For the n+48 Ca wave functions,
the solid line is obtained from a local calculation with NLAT, the dotted line is a nonlocal
calculation with β = 0.05 fm, and the dashed line is obtained from FRESCO. For the
deuteron bound wave function, the dot-dashed line results from a local calculation using
NLAT, and the open circles are from FRESCO. For all calculations we used a step size of
0.01 fm, a matching radius of 1.5 fm, and a maximum radius of 30 fm. This model space
178
0.8
48
n+ Ca NLAT Local
0.6 48
n+ Ca NLAT β=0.05
48
n+ Ca Fresco Local
0.4 Deuteron NLAT
]
Deuteron Fresco
-1/2
φ(r) [fm 0.2
-0.2
-0.4
-0.6
0 2 4 6 8 10 12 14
r [fm]
Figure E.3: n+48 Ca bound wave function, and the deuteron bound wave function. n+48 Ca:
The solid line is obtained from NLAT, the dotted line is obtained from NLAT and setting
β = 0.05 fm, and the dashed line is from FRESCO. Deuteron: Dot-dashed line is deuteron
bound wave function obtained from NLAT, and the open circles are obtained with FRESCO.
Adiabatic Potential
Next, we check the adiabatic potential. In Fig. E.4 we show the local adiabatic potential
for d+48 Ca at Ed = 20 MeV calculated with the CH89 global optical potential [23]. The
comparison is with the code TWOFNR [24]. Panel (a) is the real part of the adiabatic
In Fig. E.5 we show elastic scattering normalized to Rutherford for 48 Ca(d, d)48 Ca at
Ed = 20 MeV when using the adiabatic potential. While the adiabatic potential is not
suitable for accurately describing elastic scattering, this comparison is to show that NLAT
calculates the adiabatic potential properly, and correctly does the scattering calculation.
The nonlocal calculation used β = 0.1 since accuracy was lost with a smaller β due to inac-
curacies in calculating the nonlocal integral. The agreement between the nonlocal and local
179
0 0
-20 My Code -3
TWOFNR
Potential [MeV]
-40
-6
-60
-9
-80
-12
-100 a) b)
-120 -15
0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10
R [fm] R [fm]
Figure E.4: The local adiabatic potential for d+48 Ca at Ed = 20 MeV calculated with NLAT
and with TWOFNR [24]. (a) Real part, (b) Imaginary part.
calculations demonstrates that the nonlocal adiabatic integral is being calculated properly
since it reduces to the local calculation in the limit of β → 0, as it should. The solid line
is a calculation done with NLAT using the local adiabatic potential, the dotted line is a
nonlocal calculation with β = 0.1 fm in the nucleon optical potentials, and the dashed line
is a local calculation done with FRESCO. For these calculations, we used a 0.01 fm step
Transfer
Next, we check the T-matrix calculation. In Fig. E.6 we show DWBA transfer cross sections
for 132 Sn(d, p)133 Sn at Ed = 50 MeV. The solid line is a calculation using NLAT, and the
dashed line is a calculation done using FRESCO. The agreement with FRESCO demon-
strates that NLAT calculates the T-matrix for (d, p) transfer reactions properly. Therefore,
as long as the wave functions going into the T-matrix are correct, the correct cross section
will be calculated. The previous checks have demonstrated that the nonlocal wave functions
180
2
10
1
10
0
σ/σR 10
-1
10
NLAT Local
-2
NLAT β=0.1
10 Fresco
-3
10 0 30 60 90 120 150 180
θc.m. [deg]
Figure E.5: 48 Ca(d, d)48 Ca at Ed = 20 MeV. The solid line is when using the local adiabatic
potential, the dotted line is when doing a nonlocal calculation with β = 0.1 fm in the nucleon
optical potentials, and the dashed line is a calculation done in FRESCO.
being calculated are correct, so we can trust that the transfer results when using nonlocal
potentials will be correct as well. For this calculation we used a step size of 0.01 fm, a
2
10
10
1 NLAT
Fresco
dσ/dΩ [mb/sr]
0
10
-1
10
-2
10
-3
10 0 30 60 90 120 150 180
θc.m. [deg]
Figure E.6: 132 Sn(d, p)133 Sn at Ed = 50 MeV. Solid line is a local DWBA calculation with
NLAT, the dashed line is a calculation done with FRESCO.
181
Nonlocal Source
Finally, we examine the nonlocal source. In Sec. E we showed that the nonlocal adiabatic
source is calculated accurately for β ≈ 0 fm. For larger β values, a check with Mathematica
[116] was done, and the results of this comparison are shown in Table E.1. For this compar-
ison, we used analytic expressions for the wave functions that mimicked the behavior of the
2 −0.3r
φ(r) = e , (E.1)
r+3
sin(4R) sin(3R)
χ(R) = −i , (E.2)
6R 5R
2
r
− 1.494
Vnp (r) = −72.15e . (E.3)
There is one additional complication, namely, in order to calculate the T-matrix accu-
rately, we would like our d + A scattering wave function to be calculated in steps of 0.01 fm.
To do this, we need to know our source term S(R) (the rhs of Eq.(2.49)) in steps of 0.01 fm
as well. However, it requires a significant amount of computer time to calculate S(R) with
182
L R Mathematica NLAT
β = 0.45 fm 0 0.05 13.70+i16.53 13.71+i16.54
0 2.00 1.69+i0.78 1.69+i0.78
0 5.00 0.41-i0.29 0.41-i0.29
1 0.05 1.67-i2.30 1.70-i2.33
1 2.00 2.86+i1.31 2.86+i1.30
1 5.00 0.71-i0.50 0.71-i0.51
5 0.05 0.00-i0.0002 -0.015+i0.016
5 2.00 4.07+i1.62 4.08+i1.62
5 5.00 1.30-i0.91 1.29-i0.92
β = 0.85 fm 0 0.05 24.00-i23.61 24.01-i23.62
0 2.00 2.96+i1.10 2.95+i1.10
0 5.00 0.71-i0.33 0.71-i0.33
1 0.05 6.52-i6.61 6.55-i6.64
1 2.00 5.08+i1.90 5.08+i1.89
1 5.00 1.23-i0.56 1.23-i0.57
5 0.05 0.007-i0.0007 -0.02+i0.01
5 2.00 8.91+i3.29 8.92+i3.28
5 5.00 2.33-i1.06 2.32-i1.08
Table E.1: The nonlocal adiabatic integral, rhs of Eq.(2.49), calculated with Mathematica
and NLAT using analytic expressions for the wave functions and potentials.
such a fine grid. Therefore, in practice, S(R) is calculated in steps greater than 0.01 fm, and
then linear interpolation is used to construct S(R) in steps of 0.01 fm. To save computer
time, we would like the step size we calculate S(R) with to be as large as possible while
still maintaining the desired level of accuracy. In Fig. E.7 we show 208 Pb(d, p)209 Pb using
various step sizes for S(R). It is seen that the larger two step sizes agree, while the step size
of 0.01 fm disagrees with the other two calculations. In fact, all calculations with a step size
ranging from 0.02 − 0.05 fm agree, and only when a step size of 0.01 fm was used did we find
When calculating the wave function numerically for high values of the angular momen-
tum, L, there are difficulties near the origin due to the large centrifugal barrier. To remedy
this problem, what is often done is the wave function is set equal to zero near the ori-
183
2
10
Step Size = 0.01
10
1 Step Size = 0.03
Step Size = 0.05
dσ/dΩ [mb/sr]
0
10
-1
10
-2
10
-3
10 0 30 60 90 120 150 180
θc.m. [deg]
Figure E.7: Angular distributions for 208 Pb(d, p)209 Pb at Ed = 50 MeV obtained by using
different step sizes to calculate the rhs of Eq.(2.49). The solid line uses a step size of 0.01
fm, the dashed line a step size of 0.03 fm, and the dotted line a step size of 0.05 fm.
gin. This is done with a parameter we will call “CutL”. From the origin to a distance of
(StepSize)×(CutL)×(L), the wave function is set equal to zero. For all the calculations
done in this study, we used CutL=2. It was suspected that the discrepancy between the
calculation with a step size of 0.01 fm and the other two in Fig. E.7 was because CutL was
not big enough. This can cause problems if we try to calculate the wave function below a
very large centrifugal barrier, because numerical inaccuracies will propagate to the rest of
To figure out which calculation in Fig. E.7 is correct, we increased the CutL parameter for
the 0.01 fm step size calculation. The results are shown in Fig. E.8. When we increased the
CutL parameter of the 0.01 fm calculation, from CutL=2 to CutL=3, the resulting angular
distribution is now in agreement with the other two calculations. Therefore, the discrepancy
was indeed due to an insufficient value for CutL. This problem was investigated for all of the
reactions studied in this thesis. In all cases, a step size of 0.05 fm for the rhs of Eq.(2.49)
184
with CutL=2 was sufficient to obtain converged results.
2
10
Step Size = 0.01, CutL=2
1 Step Size = 0.01, CutL=3
10
dσ/dΩ [mb/sr] Step Size = 0.05, CutL=2
0
10
-1
10
-2
10
-3
10 0 30 60 90 120 150 180
θc.m. [deg]
Figure E.8: Angular distributions for 208 Pb(d, p)209 Pb at Ed = 50 MeV obtained by using
different step sizes and values of a cut parameter (CutL) to calculate the rhs of Eq.(2.49).
The solid line uses a step size of 0.01 fm with CutL=2, the dashed line a step size of 0.01 fm
with CutL=3, and the dotted line a step size of 0.05 fm with CutL=2.
185
Appendix F
Direct proton capture at low relative energies needed for astrophysics are always peripheral
due to the Coulomb barrier. At the limits of E → 0 these reactions are uniquely determined
by the asymptotic normalization coefficient (ANC) of the single proton overlap function of
the final nucleus [41]. It is for this reason that the ANC method [41] has been put forth as an
indirect way of obtaining proton radiative-capture cross sections from the ANCs extracted
important reactions may be difficult or impossible since the experiment would require proton-
rich radioactive beams. However, an indirect technique has been proposed [44] which uses
information about the mirror system in order to extract the necessary ANC. The mirror
nucleus is defined as the nucleus with interchanged numbers of protons and neutrons. While
experiments on the mirror system can sometimes be performed with stable beams, and thus,
In [44, 4, 5], the ratio, R, of the proton to neutron ANC squared is determined for a
wide range of light nuclei within a microscopic cluster model (MCM). In [44] an analytic
derivation of the ratio, Ro , is presented. The ratio obtained from the MCM calculations is
in fair agreement with the predictions of the analytic formula [4, 5]. In this work, we want
186
to explore the effects of couplings induced by deformations of the core and core excitations.
The reason for relying on charge symmetry arguments rather than just calculating the
ANC directly is due to large uncertainties in the theoretical predictions of ANCs. The
individual ANCs in [44, 4, 5] are strongly dependent on the N N interaction used, but the
ratio of ANCs was found to be independent of the choice of the N N interaction, within a
few percent.
Theoretical Considerations
We consider the A = B + x model used in [117], which starts from an effective Hamiltonian
where Tr is the relative kinetic energy operator, and HB is the internal Hamiltonian of the
core. VBx is the effective interaction between the core and the valence nucleon which depends
on the B − x relative coordinate, r, and the internal degrees of freedom of the core, ξ. In
−1
r − R(θ, φ)
VBx (r) = −VW S 1 + exp (F.2)
a
where VW S is the depth, and may depend on the orbital angular momentum, `. The radius,
187
Q
X
R(θ, φ) = RW S 1 + βq Yq0 (θ, φ) (F.3)
q=2
where βq characterized the deformation of the core, and thus, the strength of the coupling
with A the mass number of the B + x system. We also include the typical undeformed
The B + x wave function is expanded in eigenstates of the core, ΦI π B , with spin I, parity
πB , and eigenenergy I π B :
X
ΨJ π = ψn`j (r)Y`j (r̂)ΦI π B (ξ). (F.4)
n`jIπB
In this expansion, we factorize the radial part, ψn`j , and the spin-angular part, Y`j for
convenience. The quantum numbers n and j correspond to the principal quantum number
and the angular momentum obtained from coupling the orbital angular momentum, `, with
the spin, s, respectively. With this expansion, the coupled-channels equation for each ψ is
given by [117]:
h i
Tr` + Vii (r) ψi (r) + Vij ψj (r) = xJ π − i ψi (r),
X
(F.5)
j6=i
where i represents all possible (n`jIπB ) combinations, xJ π is the binding energy in the
188
Vij (r) = hΦi (ξ)Yi (r̂)|VBx (r, ξ)|Yj (r̂)Φj (ξ)i. (F.6)
We take Φi from the rotational model with parameters fixed phenomenologically. The
Z ∞
Six = |ψi |2 r2 dr (F.7)
0
p
with κi = 2µBx |J π − i |/~2 and µBx is the reduced mass. Here, W is the Whittaker
As an example to illustrate the model, consider the mirror pair 17 O and 17 F. The core
for both nuclei is 16 O, which has a 0+ ground state, and two low-lying 2+ and 3− states,
which strongly couple to the ground state through E2 and E3 transitions, respectively. If
we include the 0+ and 2+ states of 16 O in our model space, then the ground, 5/2+ state of
17 O and 17 F would contain not only a 1d5/2 valence nucleon coupled to the ground state,
but also, for example, a 2s1/2 nucleon coupled to the excited 2+ state.
189
p
In this study, we compare the proton ANCs, Ci , with the neutron ANCs, Cin , through
the ratio
p 2
C
R = in . (F.9)
Ci
The ratio of ANCs calculated in our model is then compared with the analytic formula
p 2
F` (iκi RN )
Ro = p , (F.10)
κi RN j` (iκn
i RN )
where F` and j` are regular Coulomb and spherical Bessel functions, respectively, [92], and
RN = 1.25A1/3 is the radius of the nuclear interior, of which Ro is not strongly dependent.
We will compare the ratio of the ANCs from our calculations with the value obtained from
this relation.
Results
Since we are interested in the ANCs for each mirror nuclei, and these depend strongly on the
energy of the system relative to threshold, it is important that we reproduce the experimental
separation energies exactly. We do this by adjusting the depths of VBn and VBp to reproduce
exactly the corresponding binding energies. All calculations are performed with the program
190
Nuclei I πB n`j R Ro RM CM
8 Li/8 B 3/2− , 1/2− 1p3/2 1.04 ± 0.04 1.12 1.08
13 C/13 N 0+ , 2+ 1p1/2 1.19 ± 0.02 1.20 1.14
17 O/17 F (g.s.) 0+ , 3− 1d5/2 1.18 ± 0.01 1.22 1.19
17 O/17 F (e.s.) 0+ , 3− 2s1/2 693 ± 16 799 736
17 O/17 F (g.s.) 0+ , 2+ 1d5/2 1.219 ± 0.004 1.22 1.19
17 O/17 F (e.s.) 0+ , 2+ 2s1/2 756 ± 23 799 736
23 Ne/23 Al 0+ , 2+ , 4+ 1d5/2 (1.852 ± 0.014) × 104 2.06 × 104 2.96 × 104
27 Mg/27 P 0+ , 2+ , 4+ 2s1/2 40.1 ± 1.8 43.7 44.3
Table F.1: Ratio of proton to neutron ANCs for the dominant component: Comparison of
this work R with the results of the analytic formula R0 Eq.(F.10) and the results of the
microscopic two-cluster calculations RM CM [4, 5] including the Minnesota interaction. The
uncertainty in R account for the sensitivity to the parameters of VBx .
FACE [120], and all details of the calculations for each mirror pair can be found in [25].
Our results are summarized in Table F.1. From the proton and neutron wave functions
calculated from Eq.(F.5), we determine the ANCs and the ratio R. For each case, R cor-
responds to rW S = 1.25 fm, a = 0.65 fm, and Vso = 6 MeV. The uncertainty reflects the
range obtained with the geometry rW S = 1.2 fm, a = 0.5 fm, and Vso = 8 MeV. Our results
for R are compared to the values obtained with the analytic formula Ro of Eq.(F.10), and
those obtained within the MCM, where they assumed two clusters and used the Minnesota
For nearly all cases, R, Ro , and RM CM are all in fair agreement. However, there
were a few cases where there were discrepancies. For 23 Ne/23 Al, it is important to note
that in our calculations we impose realistic binding energies, whereas in the MCM results,
binding energies can sometimes differ significantly. Since R depends strongly on the binding
energies, this can cause large differences between our values and those of [5]. The values
for Ro presented in Table F.1 also assume the experimental binding energies, therefore,
differences between R and Ro must be related to the failure of the simple analytic relations.
191
Such is the case for 17 O/17 F (e.s.) when the 3− excited state is included in the model space.
The usefulness of the ratio method is that the ratio R should be model independent. This
was demonstrated in Section F where for most cases studied, the ratio obtained from this
study, R, using our simple model was in fair agreement with the ratio obtained with the
much more sophisticated microscopic cluster model, RM CM . In this subsection, we use the
deformation parameter as a free variable to tune the amount of coupling between the various
configurations. With the configurations of 23 Ne/23 Al and 27 Mg/27 P pairs being very similar
to the 17 O/17 F systems in its ground and excited states, respectively, we concentrate of the
We find no significant difference in the ratio R for both the 8 Li/8 B and 13 C/13 N mirror
pairs. In these cases, the main components of the wave function are p waves, even in the
configurations including core excitation. For |β2 | = 0.0−0.7, the resulting range of values for
R are 1.038−1.044 for 8 Li/8 B and 1.201−1.251 for 13 C/13 N. This constancy is obtained even
x
though the variation in β2 leads to significant changes in the spectroscopic factor: S1p
3/2
goes from 1 to 0.75 for 8 Li/8 B, while S1p
x decreases down to 0.32 for 13 C/13 N.
1/2
The situation for 17 O/17 F differs. We consider the separate effects of including the 3−
state and the 2+ state. Let us first consider the inclusion of 16 O(0+ , 3− ). Like for 8 Li/8 B and
13 C/13 N, the variation in R was small, even though over 30% of the 5/2+ ground-state wave
We next consider the inclusion of 16 O(0+ , 2+ ). In this case, the d5/2 ground state admixes
192
1.0 1.0
neutron neutron
proton 0.9 proton
0.9 +
5/2
+
0.8 1/2
S
S
0.7
0.8
0.6
(a) (b)
0.7 0.5
0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8
β β
Figure F.1: Neutron and proton spectroscopic factors for 17 O and 17 F, respectively, consid-
ering the 16 O core in its 0+ ground state and 2+ first excited state: (a) 5/2+ ground state
and (b) 1/2+ first excited state. Figure reprinted from [25] with permission.
1.4 900
+ +
R (a) 5/2 R (b) 1/2
* *
R R
R0 850 R0
1.3
800
1.2
750
1.1 700
0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8
β β
Figure F.2: Ratio of proton and neutron ANCs for 17 O and 17 F, respectively, including
16 O(0+ , 2+ ): (a) 5/2+ ground state and (b) 1/2+ first excited state. Figure reprinted from
[25] with permission.
with an s1/2 component with the core in its excited state. For the 1/2+ excited state of 17 O,
the s1/2 component coupled to the ground-state of the core admixes with d components with
the core in its 2+ excited state. Like in all cases, the energies of the two lowest states in
17 O and 17 F were refitted by simultaneously adjusting the depths of the potential for each
β2 . For both the 5/2+ and 1/2+ states, the spectroscopic factors, of the component with
the core in the ground state experiences a large reduction at large β2 , as is seen in Fig. F.1.
193
For the ground state (5/2+ ), the proton and neutron spectroscopic factors vary together,
while for the excited state (1/2+ ), there is more admixture in the neutron system than for
the proton system. This is reflected in a different behavior of the ANC ratios.
In Fig. F.2 we present the ratio R as well as a modified ratio compensating for the
shown by the horizontal dashed line. For the 5/2+ ground state, neither R nor R∗ deviate
Fig. F.2a. Both of these ratios are close to the analytic prediction, Ro . On the contrary, for
the 1/2+ excited state, R shows a large variation partly caused by the difference between
neutron and proton spectroscopic factors, as seen in Fig. F.2b. The features seen in Fig. F.2
can be extrapolated to 23 Al and 27 P, since, as mentioned before, the former has a structure
very similar to that of 17 F(g.s.), while the latter exhibits the same components as 17 F(e.s.).
In [4, 5] core excitation is explored within the MCM. Even in these studies there was
growing disagreement between RM CM and R0 as more core states were explicitly included
in the model space. This was understood in terms of the long range Coulomb quadruple
term which was added to the Hamiltonian in the proton case, a term not considered in
the derivation of R0 , nor in our present calculations. Here, however, we not only see a
deviation from R0 , but also a strong dependence on the deformation parameter for some
cases. Therefore we conclude the source for deviations from R0 and the breakdown of the
constant ratio concept is induced by the nuclear quadruple term, which is present in both
The surprising results for the 1/2+ mirror states led to several additional tests which
isolated the cause for the large coupling dependence on R. There are three essential ingredi-
ents: low binding, the existence of an s-wave component coupled to the ground state of the
194
core, and a significant admixture with other configurations. It appears that when all three
conditions are met, the differences between the neutron and proton wave functions increase
around the surface, exactly where the nuclear quadruple interaction peaks. This results in a
stronger effect of coupling on the neutron system compared to the proton system, inducing
tests show that the effect is independent of whether the wave functions have a node.
Conclusions
The proposed indirect method for extracting proton capture rates from neutron mirror part-
ners relies on the ratio between asymptotic normalization coefficients of the mirror state
being model independent. In [25], we tested this idea against core deformation and excita-
tion. We considered a core + N model where the core is deformed and allowed to excite, and
applied it to a variety of mirror pairs (8 Li/8 B, 13 C/13 N, 17 O/17 F, 23 Ne/23 Al, and 27 Mg/27 P)
For most cases, the ratio of the ANC of mirror states was found to be independent of the
deformation, and the calculated ratio of ANCs agreed well with the simple analytic formula.
From our investigations we concluded that there are three conditions that need to be met
for the idea of a model-independent ratio to break down with deformation or core excitation:
(i) the proton system should have very low binding, (ii) the main configuration should be
an s-wave component coupled to the ground state of the core, and (iii), there should be
significant admixture with other configurations. This has implications for the application of
the indirect method based on the ANC ratio of reactions relevant to novae, namely pertaining
195
Appendix G
List of Acronyms
196
BIBLIOGRAPHY
197
BIBLIOGRAPHY
[1] F. Perey and B. Buck, “A non-local potential model for the scattering of neutrons by
nuclei,” Nuclear Physics, vol. 32, no. 0, pp. 353 – 380, 1962.
[2] Y. Tian, D.-Y. Pang, and Z.-Y. Ma, “Systematic nonlocal optical model potential for
nucleons,” International Journal of Modern Physics E, vol. 24, no. 01, p. 1550006,
2015.
[5] N. Timofeyuk, P. Descouvemont, and R. Johnson, “Relation between proton and neu-
tron asymptotic normalization coefficients for light mirror nuclei and its relevance for
nuclear astrophysics,” The European Physical Journal A - Hadrons and Nuclei, vol. 27,
no. 1, pp. 269–276, 2006.
[7] M. Chowdhury and H. S. Gupta, “A study of nuclear structure in 59ni from the (d, p)
reaction on 58 Ni,” Nuclear Physics A, vol. 205, no. 3, pp. 454 – 474, 1973.
[9] J. Lombardi, R. Boyd, R. Arking, and A. Robbins, “Nuclear sizes in 40,44,48 Ca,”
Nuclear Physics A, vol. 188, no. 1, pp. 103 – 114, 1972.
198
[11] H. Mahzoon, M. J. Charity, R. H. Dickhoff, W. H. Dussan, and J. Waldecker, S.
“Forging the link between nuclear reactions and nuclear structure,” Phys. Rev. Lett.,
vol. 112, p. 162503, Apr 2014.
[12] L. J. Titus and F. M. Nunes, “Testing the perey effect,” Phys. Rev. C, vol. 89,
p. 034609, Mar 2014.
[13] J. F. Dicello, G. Igo, W. T. Leland, and F. G. Perey, “Differential elastic cross sections
for protons between 10 and 22 MeV on 40 Ca,” Phys. Rev. C, vol. 4, pp. 1130–1138,
Oct 1971.
[18] K. Hirota, Y. Aoki, N. Okumura, and Y. Tagishi, “Deuteron elastic scattering and (d,p)
reactions on 208 Pb at Ed = 22 MeV and j-dependence of {T20} in (d,p) reaction,”
Nuclear Physics A, vol. 628, no. 4, pp. 547 – 579, 1998.
199
[21] N. K. Timofeyuk and R. C. Johnson, “Nonlocality in the adiabatic model of A(d, p)B
reactions,” Phys. Rev. C, vol. 87, p. 064610, Jun 2013.
[26] S. T. Butler, “On angular distributions from (d, p) and (d, n) nuclear reactions,” Phys.
Rev., vol. 80, pp. 1095–1096, Dec 1950.
[27] H. R., “Physics at Liverpool,” Nature, vol. 166, pp. 552–555, Sept 1950.
[28] H. R. and H. Newns, “Inelastic scattering of deuterons,” The London, Edinburgh, and
Dublin Philosophical Magazine and Journal of Science, vol. 42, no. 335, pp. 1442–1445,
1951.
[29] A. Bhatia, K. Huang, R. Huby, and H. Newns, “Angular distribution in (d, p) and
(d. n) reactions,” The London, Edinburgh, and Dublin Philosophical Magazine and
Journal of Science, vol. 43, no. 340, pp. 485–500, 1952.
[30] L. L. Lee and J. P. Schiffer, “Dependence of the angular distribution of the (d, p)
reaction on the total angular-momentum transfer,” Phys. Rev., vol. 136, pp. B405–
B409, Oct 1964.
[31] C. R. Bingham and G. T. Fabian, “Neutron shell structure in 93 Zr, 95 Zr, and 97 Zr by
(d, p) and (α, 3 He) reactions,” Phys. Rev. C, vol. 7, pp. 1509–1518, Apr 1973.
[33] L. L. Lee and J. P. Schiffer, “Spin assignments in the 2p shell from j dependence in
(d, p) angular distributions,” Phys. Rev., vol. 154, pp. 1097–1100, Feb 1967.
200
[34] C. R. Bingham and D. L. Hillis, “Neutron shell structure in 125 Sn by (d, p) and (α,3 He)
reactions,” Phys. Rev. C, vol. 8, pp. 729–736, Aug 1973.
[35] J. J. Kolata and W. W. Daehnick, “Shell-model structure of 210 Bi: 209 Bi(d, p) at 17
MeV,” Phys. Rev. C, vol. 5, pp. 568–578, Feb 1972.
201
[41] H. M. Xu, C. A. Gagliardi, R. E. Tribble, A. M. Mukhamedzhanov, and N. K. Tim-
ofeyuk, “Overall normalization of the astrophysical S factor and the nuclear vertex
constant for 7 Be(p,γ)8 B reactions,” Phys. Rev. Lett., vol. 73, pp. 2027–2030, Oct 1994.
[43] T. Fukui, K. Ogata, and M. Yahiro, “Breakup and finite-range effects on the 8 B(d, n)9 C
reaction,” Phys. Rev. C, vol. 91, p. 014604, Jan 2015.
[45] B. Guo, Z. Li, W. Liu, X. Bai, G. Lian, S. Yan, B. Wang, S. Zeng, J. Su, and Y. Lu,
“The 8 Li(d,p)9 Li reaction and astrophysical 8 B(p,γ)9 C reaction rate,” Nuclear Physics
A, vol. 761, no. 1-2, pp. 162 – 172, 2005.
[47] M. Goldberger and K. Watson, Collision Theory. John Wiley & Sons, 1964.
[49] A. Koning and J. Delaroche, “Local and global nucleon optical models from 1 keV to
200 MeV,” Nucl.Phys., vol. A713, pp. 231–310, 2003.
[50] J. S. Bell and E. J. Squires, “A formal optical model,” Phys. Rev. Lett., vol. 3, pp. 96–
97, Jul 1959.
[51] W. Dickhoff and D. Van Neck, Many-Body Theory Exposed! 2nd ed. World Scientific,
New Jersey, 2008.
[52] H. Feshbach, “Unified theory of nuclear reactions,” Annals of Physics, vol. 5, no. 4,
pp. 357 – 390, 1958.
202
[53] H. Feshbach, “A unified theory of nuclear reactions II,” Annals of Physics, vol. 19,
no. 2, pp. 287 – 313, 1962.
[56] N. C. Francis and K. M. Watson, “The elastic scattering of particles by atomic nuclei,”
Phys. Rev., vol. 92, pp. 291–303, Oct 1953.
[57] A. Kerman, H. McManus, and R. Thaler, “The scattering of fast nucleons from nuclei,”
Annals of Physics, vol. 8, no. 4, pp. 551 – 635, 1959.
[59] J. Hufner and C. Mahaux, “Low density expansion of the nucleon-nucleus optical-model
potential,” Annals of Physics, vol. 73, no. 2, pp. 525 – 577, 1972.
[60] B. D. DAY, “Elements of the brueckner-goldstone theory of nuclear matter,” Rev. Mod.
Phys., vol. 39, pp. 719–744, Oct 1967.
[62] J.-P. Jeukenne, A. Lejeune, and C. Mahaux, “Optical-model potential in finite nuclei
from Reid’s hard core interaction,” Phys. Rev. C, vol. 16, pp. 80–96, Jul 1977.
203
[65] W. H. Dickhoff, D. Van Neck, S. J. Waldecker, R. J. Charity, and L. G. Sobotka,
“Nonlocal extension of the dispersive optical model to describe data below the fermi
energy,” Phys. Rev. C, vol. 82, p. 054306, Nov 2010.
[67] J. W. Holt, N. Kaiser, G. A. Miller, and W. Weise, “Microscopic optical potential from
chiral nuclear forces,” Phys. Rev. C, vol. 88, p. 024614, Aug 2013.
[69] G. Rawitscher, “The microscopic Feshbach optical potential for a schematic coupled
channel example,” Nuclear Physics A, vol. 475, no. 3, pp. 519 – 547, 1987.
[70] K. Amos, L. Canton, G. Pisent, J. Svenne, and D. van der Knijff, “An algebraic solution
of the multichannel problem applied to low energy nucleon-nucleus scattering,” Nuclear
Physics A, vol. 728, no. 1-2, pp. 65 – 95, 2003.
[71] M. Giannini and G. Ricco, “An energy-independent nonlocal potential model for bound
and scattering states,” Annals of Physics, vol. 102, no. 2, pp. 458 – 492, 1976.
[73] N. Austern, “Wave functions of nonlocal potentials: The Perey effect,” Phys. Rev.,
vol. 137, pp. B752–B757, Feb 1965.
[74] N. Austern, Direct Nuclear Reaction Theories. John Wiley & Sons, 1970.
[75] H. Fiedeldey, “The equivalent local potential and the Perey effect,” Nuclear Physics,
vol. 77, no. 1, pp. 149 – 156, 1966.
204
[77] M. Coz, L. G. Arnold, and A. D. MacKellar, “Nonlocal potentials and their local
equivalents,” Annals of Physics, vol. 59, no. 1, pp. 219 – 247, 1970.
[78] B. T. Kim and T. Udagawa, “Method for nonlocal optical model calculations,” Phys.
Rev. C, vol. 42, pp. 1147–1149, Sep 1990.
[79] B. Kim, M. Kyum, S. Hong, M. Park, and T. Udagawa, “NLOM - a program for non-
local optical model calculations,” Computer Physics Communications, vol. 71, no. 1-2,
pp. 150 – 158, 1992.
[81] N. Michel, “A simple and efficient numerical scheme to integrate non-local potentials,”
The European Physical Journal A, vol. 42, no. 3, pp. 523–527, 2009.
[84] W. W. Daehnick, J. D. Childs, and Z. Vrcelj, “Global optical model potential for elastic
deuteron scattering from 12 to 90 MeV,” Phys. Rev. C, vol. 21, pp. 2253–2274, Jun
1980.
[86] R. Johnson and P. Tandy, “An approximate three-body theory of deuteron stripping,”
Nuclear Physics A, vol. 235, no. 1, pp. 56 – 74, 1974.
205
[87] N. B. Nguyen, F. M. Nunes, and R. C. Johnson, “Finite-range effects in (d, p) reac-
tions,” Phys. Rev. C, vol. 82, p. 014611, Jul 2010.
[89] F. M. Nunes and A. Deltuva, “Adiabatic approximation versus exact Faddeev method
for (d, p) and (p, d) reactions,” Phys. Rev. C, vol. 84, p. 034607, Sep 2011.
[91] A. Deltuva, “Three-body direct nuclear reactions: Nonlocal optical potential,” Phys.
Rev. C, vol. 79, p. 021602, Feb 2009.
[92] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions. Dover, New York,
1970.
[96] A. M. Moro, F. M. Nunes, and R. C. Johnson, “Theory of (d, p) and (p, d) reactions
including breakup: Comparison of methods,” Phys. Rev. C, vol. 80, p. 064606, Dec
2009.
206
[98] F. D. Becchetti and G. W. Greenlees, “Nucleon-nucleus optical-model parameters, A
> 40, E < 50 MeV,” Phys. Rev., vol. 182, pp. 1190–1209, Jun 1969.
[101] F. Perey and D. Saxon, “The local energy approximation to nonlocality and finite
range effects,” Physics Letters, vol. 10, no. 1, pp. 107 – 109, 1964.
[102] C. Mahaux and R. Sartor, “Calculation of the shell-model potential from the optical-
model potential,” Phys. Rev. Lett., vol. 57, pp. 3015–3018, Dec 1986.
[106] A. Ioannides and R. Mackintosh, “A method for s-matrix to potential inversion at fixed
energy: (i). method description and evaluation,” Nuclear Physics A, vol. 438, no. 2,
pp. 354 – 383, 1985.
[109] A. Lovell and K. Amos, “Local and nonlocal equivalent potentials for p−12 C scatter-
ing,” Phys. Rev. C, vol. 62, p. 064614, Nov 2000.
207
[110] R. V. Reid, “Local phenomenological nucleon-nucleon potentials,” Annals of Physics,
vol. 50, no. 3, pp. 411 – 448, 1968.
[111] M. B. Tsang, J. Lee, and W. G. Lynch, “Survey of ground state neutron spectroscopic
factors from Li to Cr isotopes,” Phys. Rev. Lett., vol. 95, p. 222501, Nov 2005.
[112] G. Kramer, H. Blok, and L. Lapiks, “A consistent analysis of (e, e‘p) and (d,3 He)
experiments,” Nuclear Physics A, vol. 679, no. 3-4, pp. 267 – 286, 2001.
[114] A. Ross, L. J. Titus, and F. M. Nunes, “Effects of nonlocal potentials on (p, n) transfer
reactions,” In Preparation, 2015.
[117] F. Nunes, I. Thompson, and R. Johnson, “Core excitation in one neutron halo sys-
tems,” Nuclear Physics A, vol. 596, no. 2, pp. 171 – 186, 1996.
[120] I. Thompson, F. Nunes, and B. Danilin, “FaCE: a tool for three body Faddeev cal-
culations with core excitation,” Computer Physics Communications, vol. 161, no. 1-2,
pp. 87–107, 2004.
208