Newer
Older
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% Scattering Theory %
% %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\chapter{Scattering Theory}\label{ch:scattering_theory}
Details of scattering theory have been well understood since the middle of the previous century, when they were summarized in a review article by Lane and Thomas [AL58]\cite{lane_thomas_1958}. A wealth of additional reference material is available to the student of scattering theory; only a few are listed here. The text by Foderaro [AF71]\cite{foderaro_1971} provides a more elementary introduction to the subject. One publication by Fr\"{o}hner [FF80]\cite{frohner_1978} is based on lectures presented at the International Centre for Theoretical Physics (ICTP) Winter Courses on Nuclear Physics and Reactors, 1978; this is a comprehensive and useful guide to applied neutron resonance theory. It includes a variety of topics, including preparation of data, various approximations to scattering theory, Doppler broadening, experimental complications, data-fitting procedures, and statistical tests. Another Fr\"{o}hner paper [FF00]\cite{jeff18_frohner} is somewhat more theoretical, and covers many aspects of data fitting in the resonance region.
The particular aspect of scattering theory with which we are concerned is the R-matrix formalism. A summary of the underlying principles is given here.
R-matrix theory is a mathematically rigorous phenomenological description of what is actually seen in an experiment (i.e., the measured cross section). The theory is not a model of neutron-nucleus interaction, in the sense that it makes no assumptions about the underlying physics of the interaction. Instead it parameterizes the measurement in terms of quantities such as the interaction radii and boundary conditions, resonance energies and widths, and quantum numbers; values for these parameters may be determined by fitting theoretical calculations to observed data. The theory is mathematically correct, in that it is analytic, unitary, and rigorous; nevertheless, in practical applications, the theory is always approximated in some fashion.
R-matrix theory is based on the following assumptions\footnote{In practical applications two of these four assumptions may be violated in one degree or another: (1) The theory may be used for relativistic neutron energies, and corrected for relativistic effects; nevertheless, non-relativistic quantum mechanics is assumed. (2) A fission experiment with more than two final products is treated as a two-step process. That is, the immediate result of the neutron-nuclide interaction is assumed to be limited to two final products, at least one of which decays prior to detection.}:
(1) the applicability of non-relativistic quantum mechanics;
(2) the absence or unimportance of all processes in which more than two product nuclei are formed;
(3) the absence or unimportance of all processes of creation or destruction; and
(4) the existence of a finite radial separation beyond which no nuclear interactions occur, although Coulomb interactions are given special treatment.
R-matrix theory is expressed in terms of channels, where a channel is defined as a pair of (incoming or outgoing) particles, plus specific information relevant to the interaction between the two particles. A schematic depicting entrance and exit channels is shown in Fig. \ref{scattering-theory_rmatrix_channel_diagram}. Note that entrance channels can also occur as exit channels, but some exit channels (e.g., fission channels) do not occur as entrance channels. Two interacting particles are shown in the portion of the figure that is labeled ``Interior Region''; here the particles are separated by less than the interaction radius $a$.
\begin{figure}
\centering
\includegraphics[width=0.4\textwidth]{figures/rmatrix_inc_channel.pdf} \\[\smallskipamount]
\includegraphics[width=0.4\textwidth]{figures/rmatrix_interior_region.pdf}\hfill
\includegraphics[width=0.4\textwidth]{figures/rmatrix_exit_channel.pdf}
\caption{Schematic of entrance and exit channels as used in scattering theory. For the interior region (with separation distance $r < a$), no assumptions are made about the nature of the interaction. In the figure, $m$, $i$, and $z$ refer to the mass, spin, and charge of the incident particle while $M$, $I$ and $Z$ refer to the target particle. Orbital angular momentum is denoted by $l$ and velocity by $v$. Primes are used for post-collision quantities.}
\label{scattering-theory_rmatrix_channel_diagram}
\end{figure}
In Section \ref{sec:equations-for-scattering-theory}, general equations of scattering theory are presented and their derivations discussed. The fundamental R-matrix equations are presented. Section II.A.1\ref{} gives a detailed derivation of the equations for a simple case. Section II.A.2\ref{} shows the relationship between the R-matrix and the A-matrix, which is another common representation of scattering theory.
The approximations to R-matrix theory available in the SAMMY code are detailed in Section II.B\ref{}. The recommended choice for most applications is the Reich-Moore approximation, described in Section II.B.1\ref{}. For some applications, the Reich-Moore approximation is inadequate; for those cases, a method for using SAMMY's Reich-Moore approximation to mimic the full (exact) R-matrix is presented Section II.B.2\ref{}. Two historically useful but now obsolete approximations are single-level and multilevel Breit Wigner (SLBW and MLBW), discussed in Section II.B.3\ref{}. Provisions for including non-compound (direct) effects are discussed in Section II.B.4\ref{}.
In Section II.C\ref{}, details are given for the SAMMY nomenclature and other conventions,
for transformations to the center-of-momentum system, and for the calculation of penetrability,
shift factors, and hard-sphere phase shifts in both Coulomb and non-Coulomb cases.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Equations For Scattering theory} \label{sec:equations-for-scattering-theory}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
In this section, equations for scattering theory are presented but not derived. Specifics for
the R-matrix formulation of scattering theory are presented in Section II.A.1\ref{}, which provides a
discussion of an alternative formulation (the A-matrix). Readers interested in the derivation of
the equations for scattering theory are referred to the Lane and Thomas article \cite{lane_thomas_1958} for a
detailed derivation in the general case, or to Section II.A.2\ref{} of this document for a simplified
version.
In scattering theory, a channel may be defined by $c = (\alpha, l, s, J)$, where the following
definitions apply:
\begin{itemize}
\item $\alpha$ represents the two particles making up the channel; $\alpha$ includes mass ($m$ and $M$), charge ($z$ and $Z$), spin ($i$ and $I$ ) with associated parities, and all other quantum numbers for each of the two particles, plus the Q-value (equivalent to the negative of the threshold energy in the
center of momentum system).
\item $l$ is the orbital angular momentum of the pair, and the associated parity is given by $(-1)^l$.
\item $s$ represents the channel spin (including the associated parity); that is, $s$ is the quantized
vector sum of the spins of the two particles of the pair: $\vec{s} = \vec{i} + \vec{I}$
\item $J$ is the total angular momentum (and associated parity); that is, $J$ is the quantized vector sum
of $l$ and $s$: $\vec{J} = \vec{l} + \vec{s}$.
\end{itemize}
Only $J$ and its associated parity $\pi$ are conserved for any given interaction. The other quantum numbers may differ from channel to channel, as long as the sum rules for spin and parity are obeyed. Within this document and within the SAMMY code, the set of all channels with the same $J$ and $\pi$ are called a ``spin group.''
In all formulae given below, spin quantum numbers (e.g., $J$ ) are implicitly assumed to include the associated parity. Quantized vector sum rules are implicitly assumed to be obeyed. Readers unfamiliar with these sum rules are referred to Section II.C.1.a\ref{} for a mini-tutorial on the subject.
Let the angle-integrated cross sections from entrance channel $c$ to exit channel $c'$ with total angular momentum $J$ be represented by $\sigma_{cc'}$. This cross section is given in terms of the scattering matrix $U_{cc'}$ as
\begin{equation}\label{eq:sigma-ccprime}
\sigma_{cc'} = \frac{\pi}{k_\alpha^2}g_{J\alpha}\left|e^{2iw_c}\delta_{cc'}-U_{cc'}\right|^2\delta_{JJ'} \:,
\end{equation}
\noindent
where $k_\alpha$ is the wave number (and $K_\alpha = \hbar k_\alpha =$ center-of-mass momentum) associated with incident particle pair $\alpha$, $g_{J\alpha}$ is the spin statistical factor, and $w_c$ is the Coulomb phase-shift difference. Note that $w_c$ is zero for non-Coulomb channels. (Details for the charged-particle case are presented in Section II.C.4.\ref{}) The spin statistical factor $g_{J\alpha}$ is given by
\begin{equation}\label{eq:spin-stat-factor}
\end{equation}
\noindent
and center-of-mass momentum $K_\alpha$ by
\begin{equation}\label{eq:center-of-mass-momentum}
K_\alpha^2 = (\hbar k_\alpha)^2 = \frac{2mM^2}{(m+M)^2}E \:.
\end{equation}
\noindent
Here $E$ is the \textbf{laboratory} kinetic energy of the incident (moving) particle. A derivation of this value for $K_\alpha$ is given in Section II.C.2\ref{}.
The scattering matrix $U$ can be written in terms of matrix $W$ as
\begin{equation}\label{eq:scat-matrix}
\end{equation}
\noindent
where $\Omega$ is given by
\begin{equation}\label{eq:omega}
\end{equation}
\noindent
Here again, $w_c$ is zero for non-Coulomb channels, and the potential scattering phase shifts for non-Coulomb interactions $\phi_c$ are defined in many references (e.g., \cite{lane_thomas_1958}) and shown in Table \ref{}. The matrix $W$ in Eq. \ref{eq:scat-matrix} is related to the R-matrix (in matrix notation with indices suppressed) via
\begin{equation}\label{eq:W-matrix}
W = P^{1/2}\left(I-RL\right)^{-1}\left(I-RL^*\right)P^{-1/2} \:.
The quantity $I$ in this equation represents the identity matrix, and superscript $*$ indicates a complex conjugate. The form of the R-matrix is given in Section IIA.1\ref{} in general Section II.B\ref{} for the versions used in SAMMY. The quantity $L$ in Eq. \ref{eq:W-matrix} is given by
\end{equation}
\noindent
with $P$ being the penetration factor (penetrability) $S$ the shift factor, and $B$ the arbitrary boundary constant at the channel radius $a_c$. $P$ and $S$ are functions of energy $E$, and also depend on the orbital angular momentum $l$ and the channel radius $a_c$. Formulae for $P$ and $S$ are found in many references (see, for example Eq. (2.9) in \cite{lynn_1968}).
For non-Coulomb interactions, the penetrability and shift factor have the form
\begin{equation}\label{eq:pen-shift-func-of-rho}
P\rightarrow P_l(\rho) \qquad \text{and} \qquad S \rightarrow S_l(\rho) \:,
\end{equation}
\noindent
where $\rho$ is related to the center-of-mass momentum which in turn is related to the laboratory energy of the incident particle $(E)$. For arbitrary channel $c$ with a particle pair $\alpha$, orbital angular momentum $l$, and channel radius $a_c$, $\rho$ has the form
\begin{equation}\label{eq:de-Broglie-radius}
\rho = k_\alpha a_c = \frac{1}{\hbar} \sqrt{\frac{2m_\alpha M_\alpha}{m_\alpha+M_\alpha} \frac{M}{m+M}} \sqrt{(E-\Xi_\alpha)}\: a_c \:,
\end{equation}
as shown in Section II.C.2\ref{}. Here $\Xi_\alpha$ is the energy threshold for the particle pair $\alpha$, $m_\alpha$ and $M_\alpha$ are the masses of the two particles of particle pair $\alpha$, and $m$ and $M$ are the masses of the incident particle and target nuclide, respectively.
Appropriate formulae for $P$, $S$, and $\phi$ in the non-Coulomb case are shown in Table IIA.1. For two charged particles, formulae for the penetrabilities are given in Section II.C.4\ref{tab:penetrabilities}.
The energy dependence of fission and capture widths is negligible over the energy range of these calculations. Therefore, a penetrability of unity may be used.
% multiline cell: \begin{tabular}{@{}c@{}} line1 \\ line2 \end{tabular}
\begin{threeparttable}
\centering
\caption{Hard-sphere penetrability (penetration factor) $P$, level shift factor $S$, and potential-scattering phase shift $\phi$ for orbital angular momentum $l$, wave number $k$, and channel radius $a_c$, with $\rho=ka_c$} \label{tab:penetrabilities}
\begin{tabular}{ l c c c }
\hline\hline
$\mathbf{l}$ & $\mathbf{P_l}$ & $\mathbf{S_l}$ & $\mathbf{\phi_l}$ \cr
\hline\hline
0 & $\rho$ & 0 & $\rho$ \cr
1 & $\rho^3/(1+\rho^2)$ & $-1/(1+\rho^2)$ & $\rho-tan^{-1}\rho$ \cr
2 & $\rho^5/(9+3\rho^2+\rho^4)$ & $-(18+3\rho^2)/(9+3\rho^2+\rho^4)$ & $\rho-tan^{-1}\left[3\rho/(3-\rho^2)\right]$ \cr
3 & \begin{tabular}{@{}c@{}}$\rho^7/$ \\ $(225+45\rho^2+6\rho^4+\rho^6)$\end{tabular} & \begin{tabular}{@{}c@{}} $-(675+90\rho^2+6\rho^4)/$ \\ $(225+45\rho^2+6\rho^4+\rho^6)$ \end{tabular} & $\rho-tan^{-1}\left[\rho(15-\rho^2)/(15-6\rho^2)\right]$ \cr
\hline
$l$ & $\frac{\rho^2P_{l-1}}{(1-S_{l-1})^2+P_{l-1}^2}$ & $\frac{\rho^2(l-S_{l-1})}{(1-S_{l-1})^2+P_{l-1}^2}-l$ & $\phi_{l-1} - tan^{-1}\left(P_{l-1}/(l-S_{l-1})\right)$\tnote{\textdagger} \cr
\hline\hline
\end{tabular}
\begin{tablenotes}
\item[\textdagger] \footnotesize{The iterative formula for $\phi_l$ could also be defined by $B_l = (B_{l-1}+X_l)/(1-B_{l-1}X_l)$ where $B_l=tan(\rho-\phi_l)$ and $X_l=P_{l-1}/(l-S_{l-1})$}
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
Formulae for a particular cross section type can be derived by summing over the terms in Eq. (II A.1)\ref{eq:sigma-ccprime}. For the total cross section, the sum over all possible exit channels and all spin groups gives
\begin{equation}\label{eq:sigma-tot}
\begin{aligned}
\sigma^{total} & = \sum_{\substack{incident \\ channels \\ c }} \sum_{\substack{all \\ channels \\ c' }} \sum_J \frac{\pi}{k_\alpha^2}g_\alpha|\delta_{cc'}-U_{cc'}|^2 \\
& = \frac{\pi}{k_\alpha^2} \sum_{J}g_J \sum_{\substack{incident \\ channels \\ c }} \sum_{\substack{all \\ channels \\ c' }} \left( \delta_{cc'}-U_{cc'}\delta_{cc'}-U^*_{cc'}\delta_{cc'}+\left|U_{cc'}\right|^2 \right) \\
& = \frac{2\pi}{k_\alpha^2} \sum_J g_J \sum_{\substack{incident \\ channels \\ c }} \left( 1-\text{Re}\left(U_{cc}\right) \right).
\end{aligned}
\end{equation}
For non-charged incident particles, the elastic (or scattering) cross section is given by
\begin{equation}\label{eq:sigma-scat}
\sigma_{\alpha\alpha} = \frac{\pi}{k_\alpha^2} \sum_J g_J \sum_{\substack{ c=incident \\ channel }} \left( 1-2\text{Re}\left(U_{cc}\right) + \sum_{\substack{c'=incident \\ channel }} \left|U_{cc'}\right|^2 \right).
\end{equation}
Similarly, the cross section for any non-elastic reaction can be written
\begin{equation}\label{eq:sigma-reaction}
\sigma_\alpha^{reaction} = \frac{\pi}{k_\alpha^2} \sum_J g_J \sum_{\substack{ c=incident \\ channel }} \sum_{\substack{ c'=reaction \\ channel }} \left| U_{cc'} \right|^2.
\end{equation}
In particular, the capture cross section could be written as the difference between the total and all
other cross sections,
\begin{equation}\label{eq:sigma-capture}
\sigma_\alpha^{reaction} = \frac{\pi}{k_\alpha^2} \sum_J g_J \sum_{\substack{ c=incident \\ channel }} \left( 1 - \sum_{\substack{ c'=all\: channels \\ except\: capture }} \left| U_{cc'} \right|^2 \right).
\end{equation}
(This form will be used later, in Section II.B.1.a, when the capture channels are treated in an
approximate fashion.)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{R-Matrix and A-Matrix Equations}\label{subsec:r-matrix-a-matrix}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The R-matrix was introduced in Eq. \ref{eq:W-matrix} as
\begin{equation} % don't ref, give \ref{eq:W-matrix}
W = P^{1/2}\left(I-RL\right)^{-1}\left(I-RL^*\right)P^{-1/2} \:,
\end{equation}
but the formula for the R-matrix was not given there. If $\lambda$ represents a particular resonance (or energy level), then the general form for the R-matrix is
\begin{equation}\label{eq:r-matrix}
R_{cc'} = \sum_\lambda \frac{ \gamma_{\lambda c}\gamma_{\lambda c'} }{ E_\lambda - E } \delta_{J,J'} \:,
\end{equation}
where $E_\lambda$ represents the energy of the resonance, and the reduced width amplitude $\gamma$ is related to the partial width $\Gamma$ by
\begin{equation}\label{eq:reduced-widths}
\Gamma_{\lambda c} = 2 \bf{P_c}\gamma_{\lambda c}^2.
\end{equation}
The sum in Eq. (II A1.2)\ref{eq:r-matrix} contains an infinite number of levels. All channels, including the ``gamma channel'' for which one of the particles is a photon, are represented by the channel indices.
The R-matrix is not the only possibility for parameterization of the scattering matrix. In the R-matrix formulation, equations are expressed in terms of channel-channel interactions. It is also possible to formulate scattering theory in terms of level-level interactions; this formulation uses what is called the A-matrix, which is defined as
\begin{equation}\label{eq:a-matrix}
A_{\mu\lambda}^{-1} = \left( E_\lambda - E \right)\delta_{\mu\lambda} - \sum_c \gamma_{\mu c}L_c\gamma_{\lambda c} .
\end{equation}
\noindent
To see the relationship of the A-matrix to the R-matrix, we begin by multiplying both
sides of Eq. (II A1.4)\ref{eq:a-matrix} by $A$ and summing over $\lambda$:
\begin{equation}
\begin{aligned}
\sum_\lambda A_{\mu\lambda}^{-1}A_{\lambda\nu} & = \sum_\lambda \left(E_\lambda-E\right)\delta_{\mu\lambda}A_{\lambda\nu} - \sum_c \gamma_{\mu c} L_c \gamma_{\lambda c} A_{\lambda \nu}, \\
& \text{or} \\
\delta_{\mu\nu} & = \left(E_\mu-E\right)A_{\mu\nu} - \sum_c \gamma_{\mu c}L_c \sum_\lambda \gamma_{\lambda c}A_{\lambda\nu} .
\end{aligned}
\end{equation}
\noindent
Dividing by $(E_\mu - E)$, multiplying on the left by $\gamma_{\mu c'}$ and on the right by $\gamma_{\nu c''}$, and summing over $\mu$ puts this equation into the form
\begin{aligned}
\sum_\mu \gamma_{\mu c'}\left(E_{\mu}-E\right)^{-1}\delta_{\mu\nu}\gamma_{\nu c''} = & \sum_\mu \gamma_{\mu c'}\left(E_\mu-E\right)^{-1}\left(E_\mu-E\right)A_{\mu\nu}\gamma_{\nu c''} \\
& - \sum_\mu \gamma_{\mu c'} \left(E_{\mu}-E\right)^{-1} \sum_c \gamma_{\mu c} L_c \sum_\lambda \gamma_{\lambda c} A_{\lambda \nu}\gamma_{\nu c''} ,
\end{aligned}
\end{equation}
\noindent
which can be reduced to
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
\begin{aligned}
\gamma_{\nu c'}\left(E_{\nu}-E\right)^{-1}\gamma_{\nu c''} = & \sum_{\mu} \gamma_{\mu c'} A_{\mu\nu}\gamma_{\nu c''} \\
& - \sum_c \left[ \sum_\mu \gamma_{\mu c'}\left(E_{\mu}-E\right)^{-1}\gamma_{\mu c} \right] L_c \sum_\lambda \gamma_{\lambda c} A_{\lambda\nu} \gamma_{\nu c''}.
\end{aligned}
\end{equation}
\noindent
Summing over $\nu$ puts this into the form
\begin{equation}
\begin{aligned}
\left[ \sum_\nu \gamma_{\nu c'}\left(E_{\nu}-E\right)^{-1}\gamma_{\nu c''} \right] = & \sum_{\mu\nu} \gamma_{\mu c'} A_{\mu\nu}\gamma_{\nu c''} \\
& - \sum_c \left[ \sum_\mu \gamma_{\mu c'}\left(E_{\mu}-E\right)^{-1}\gamma_{\mu c} \right] L_c \sum_{\lambda\nu} \gamma_{\lambda c} A_{\lambda\nu} \gamma_{\nu c''},
\end{aligned}
\end{equation}
\noindent
in which we can replace the quantities in square brackets by the R-matrix, giving
\begin{equation}
\begin{aligned}
R_{c'c''} = & \sum_{\mu\nu} \gamma_{\mu c'} A_{\mu\nu}\gamma_{\nu c''} - \sum_c R_{c'c} L_c \sum_{\lambda\nu} \gamma_{\lambda c} A_{\lambda\nu} \gamma_{\nu c''}, \\
= & \sum_c \left[ \delta_{c'c} - R_{c'c}L_c \right] \sum_{\lambda\nu}\gamma_{\lambda c} A_{\lambda\nu} \gamma_{\nu c''}.
\end{aligned}
\end{equation}
\noindent
Solving for the summation, this equation can be rewritten as
\begin{equation}\label{eq:r-matrix-to-a-matrix}
\left[(I-RL)^{-1}R\right]_{cc''} = \sum_{\lambda\nu} \gamma_{\lambda c} A_{\lambda\nu} \gamma_{\nu c''}.
\end{equation}
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
\noindent
To relate this to the scattering matrix, we note that Eq. (II A.6)\ref{eq:W-matrix} can be rewritten using Eq. (II A.7)\ref{eq:L-matrix} into the form
\begin{equation}\label{eq:w-matrix-transform}
\begin{aligned}
W & = P^{1/2}\left(I-RL\right)^{-1}\left(I-RL^*\right)P^{-1/2} \\
& = P^{1/2}\left(I-RL\right)^{-1}\left(I-RL+2iRP\right)P^{-1/2} \\
& = P^{1/2}\left[\left(I-RL\right)^{-1}\left(I-RL\right) + 2i\left(I-RL\right)^{-1}RP\right]P^{-1/2} \\
& = P^{1/2}P^{-1/2} + 2iP^{1/2}\left(I-RL\right)^{-1}RP P^{-1/2} \\
& = I + 2iP^{1/2}\left(I-RL\right)^{-1}RP^{1/2} .
\end{aligned}
\end{equation}
\noindent
Comparing Eq. (II A1.10)\ref{eq:r-matrix-to-a-matrix} to Eq. (II A1.11)\ref{eq:w-matrix-transform} gives, in matrix form,
\begin{equation}
W = I + 2iP^{1/2}\gamma A\gamma P^{1/2} .
\end{equation}
\noindent
These equations are exact; no approximations have been made.
One common approximation should be discussed here: the ``eliminated channel'' approximation, for which one particular type of channel is treated in aggregate and assumed to not interfere from level to level. This is most easily understood in the A-matrix definition, Eq. (II A1.4)\ref{eq:a-matrix}; assuming no level-level interference for the gamma channels (for example), this equation can be approximated as
\begin{equation}\label{eq:a-matrix-rm-approx}
A_{\mu\lambda}^{-1} \approx (E_\lambda-E)\delta_{\mu\lambda} - \left[\sum_{\substack{\gamma=gamma \\ channels}} \gamma_{\mu\gamma}L_{\gamma}\gamma_{\lambda\gamma}\right]\delta_{\mu\lambda} - \sum_{\substack{c=particle \\channels}} \gamma_{\mu c}L_c\gamma_{\lambda c} .
\end{equation}
\noindent
The quantity in square brackets corresponds to those channels for which the level-level
interference is to be neglected; that is, only the interactions within one level are important. For
gamma channels, $L=S+iP$ reduces to $L=i$, so Eq. (II A1.13)\ref{eq:a-matrix-rm-approx} becomes
\begin{equation}\label{eq:a-matrix-rm-reduced}
A_{\mu\lambda}^{-1} \approx \left(E_\lambda-E-i\overline{\Gamma}_{\lambda\gamma}/2\right)\delta_{\mu\lambda} - \sum_{\substack{c=particle \\ channels}} \gamma_{\mu c}L_c\gamma_{\lambda c} .
The bar over $\overline{\Gamma}_{\lambda\gamma}$ is used to indicate the special treatment for this channel.
In this form, our expression for $A$ is analogous to the exact expression in Eq. (II A1.4)\ref{eq:a-matrix} with two modifications: the additional imaginary term is added to the energy difference, and the sum over the channels includes only the ``particle channels'' (non-eliminated channels). It is therefore possible to immediately write the R-matrix formula for the eliminated-channel approximation as
R_{cc'} = \sum_{\lambda}\frac{\gamma_{\lambda c}\gamma_{\lambda c'}}{E_\lambda-E-i\overline{\Gamma}_{\lambda\gamma}/2} \delta_{JJ'} ,
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
\end{equation}
where the channel indices c and c ' refer only to particle channels, not to the gamma channels.
This formula for the R-matrix is the Reich-Moore approximation and is the form which is used in
the SAMMY code. See Section II.B.1 for more about this formulation of R-matrix theory.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Derivation of Scattering Theory Equations}\label{subsec:derivation-of-scat-theory}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Many authors have given derivations of the equations for the scattering matrix in terms of the R-matrix. Sources for the derivation shown here are unpublished lecture notes of Fr\"{o}hner [FF02]\cite{frohner_lecture_notes_02}, presented at the SAMMY workshop in Paris in 2002, and Foderaro [AF71]\cite{foderaro_1971}. This derivation is valid for only the simple case of spinless projectiles and target nuclei, assuming only elastic scattering and absorption. For the general case, the reader is referred to Lane and Thomas [AL58]\cite{lane_thomas_1958}.
\noindent
\textbf{Schr\"{o}dinger equation}
The Schr\"{o}dinger equation with a complex potential is
\begin{equation}
\left( \frac{-\hbar^2}{2m}\nabla^2 + V + iW \right)\psi = E\psi,
\end{equation}
\noindent
in which one can consider that $V$ causes scattering and $W$ causes absorption. The wave function
can be expanded in the usual fashion,
\begin{equation}
\psi\left(r,\text{cos}(\theta)\right) = \sum_{l=0}^\infty \frac{u_l(r)}{r} P_l(\text{cos}(\theta)) ,
\end{equation}
\noindent
for which the radial portion obeys the equation
\frac{d^2u_l}{dr^2} + \left[ k^2-\frac{2m}{\hbar^2}(V+iW)-\frac{l(l+1)}{r^2} \right]u_l = 0,
\end{equation}
\noindent
subject to the conditions that $|\psi|^2$ is everywhere finite and that
\begin{equation}
u_l(r=0) = 0.
\end{equation}
In the external region, $r>a$, the nuclear forces are zero $(V=W=0)$, so the solution has
the form
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
u_l(r) = I_l(r) - U_l O_l(r) .
\end{equation}
$I_l$ represents an incoming free wave, and $O_l$ represents an outgoing free wave. $U_l$ is the ``collision function'' or ``S function'' that describes the effects of the nuclear interaction, giving both the attenuation and the phase shift of the outgoing wave:
\begin{equation}
\begin{aligned}
|U_l|^2 = 1 & \text{ for } W = 0, \\
|U_l|^2 < 1 & \text{ for } W \neq 0.
\end{aligned}
\end{equation}
\noindent
Our goal is to determine an appropriate analytic form for $U_l$.
\noindent
\textbf{Orthogonal eigenvectors in interior region}
For the interior region $r<a$ , we define eigenfunctions $w_{\lambda l}(r)$ and eigenvalues $E_\lambda$,
\begin{equation}
E_\lambda = \frac{\hbar^2k_\lambda^2}{2m},
\end{equation}
\noindent
for the wave equation without absorption $(W = 0)$,
\begin{equation}\label{eq:wave-func-without-abs}
\frac{d^2w_{\lambda l}}{dr^2} + \left[k_\lambda^2-\frac{2m}{\hbar^2}V-\frac{l(l+1)}{r^2}\right]w_{\lambda l} = 0,
\end{equation}
\noindent
for which the boundary conditions are
\begin{equation}\label{eq:schrod-bound-conditions}
w_{\lambda l}(r=0) = 0 \;\;\;\;\; \text{and} \;\;\;\;\; \frac{a}{w_{\lambda l}(a)}\frac{dw_{\lambda l}}{dr}\Bigr|_{r=a} = B_l.
\end{equation}
\noindent
Note that $w_{\lambda l}(r)$ is real if the boundary parameter $B_l$ is chosen to be real. The eigenfunctions are orthogonal, since
\begin{equation}\label{eq:eigen-func-orthog}
\begin{aligned}
\int_0^a\left( \frac{d^2w_{\lambda l}}{dr^2}w_{\mu l} - w_{\lambda l}\frac{d^2w_{\mu l}}{dr^2} \right) dr & = \int_0^a \frac{d}{dr} \left( \frac{dw_{\lambda l}}{dr}w_{\mu l} - w_{\lambda l}\frac{dw_{\mu l}}{dr} \right) dr \\
& = \left[\frac{dw_{\lambda l}}{dr}w_{\mu l} - w_{\lambda l}\frac{dw_{\mu l}}{dr}\right]_0^a \\
& = \frac{dw_{\lambda l}}{dr}\Bigr|_{r=a} w_{\mu l}(a) - w_{\lambda l}(a)\frac{dw_{\mu l}}{dr}\Bigr|_{r=a} - [0] \\
& = \frac{B_l}{a}\left[ w_{\lambda l}(a)w_{\mu l}(a)-w_{\lambda l}(a)w_{\mu l}(a) \right] = 0,
\end{aligned}
\end{equation}
\noindent
in which both equations of (II A2.9)\ref{eq:schrod-bound-conditions} have been invoked. The integral in Eq. (II A2.10)\ref{eq:eigen-func-orthog} can also be evaluated using Eq. (II A2.8)\ref{eq:wave-func-without-abs}, giving
\begin{equation}\label{eq:reform-orthog}
\begin{aligned}
\int_0^a\left( \frac{d^2w_{\lambda l}}{dr^2}w_{\mu l} - w_{\lambda l}\frac{d^2w_{\mu l}}{dr^2} \right) dr & = \int_0^a\left( \left[-k_\lambda^2-\frac{2mV}{\hbar^2}\right]w_{\lambda l}w_{\mu l} - w_{\lambda l}\left[-k_\lambda^2-\frac{2mV}{\hbar^2}\right]w_{\lambda l} \right)dr \\
& = \int_0^a\left( -k_\lambda^2w_{\lambda l}w_{\mu l} + k_\mu^2w_{\lambda l}w_{\mu l} \right)dr \\
& = -(k_\lambda^2-k_\mu^2)\int_0^a w_{\lambda l}w_{\mu l}dr .
\end{aligned}
\end{equation}
\noindent
Equating Eq. (II A2.10)\ref{eq:eigen-func-orthog} to Eq. (II A2.11)\ref{eq:reform-orthog} gives
\begin{equation}
(k_\lambda^2-k_\mu^2)\int_0^a w_{\lambda l}w_{\mu l}dr = 0.
\end{equation}
\noindent
For $\lambda \neq \mu$, assuming no degenerate states, it therefore follows that
\begin{equation}
\int_0^a w_{\lambda l}w_{\mu l}dr = 0 \;\;\;\;\; \text{if} \;\;\;\;\; \lambda\neq\mu.
\end{equation}
The orthogonality of the eigenvectors is therefore established. We assume that these wave
functions are normalized such that
\begin{equation}\label{eq:normalized-wave-func}
\int_0^a w_{\lambda l}w_{\mu l}dr = \delta_{\lambda\mu}.
\end{equation}
\noindent
\textbf{Matching at the surface}
The internal wave function for the true potential (including the imaginary part $iW$ ) can be
expanded in terms of the eigenfunctions as
\begin{equation}\label{eq:internal-wave-func}
u_l(r) = \sum_{\lambda} c_{\lambda l}w_{\lambda l}(r) \;\;\;\;\; \text{for} \;\; r\leq a,
\end{equation}
\noindent
with
c_{\lambda l} = \int_0^a u_l w_{\lambda l} dr.
\end{equation}
\noindent
This equation for $c_{\lambda l}$ is derived by multiplying Eq. (II A2.15)\ref{eq:internal-wave-func} by $u_{\lambda l}(r)$ , integrating, and applying Eq. (II A2.14)\ref{eq:normalized-wave-func}.
\noindent
Consider now the integral
\begin{equation}
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
\int_0^a\left( \frac{d^2u_{l}}{dr^2}w_{\lambda l} - u_{l}\frac{d^2w_{\lambda l}}{dr^2} \right) dr,
\end{equation}
\noindent
which can be expanded by use of Eqs. (II A2.3)\ref{eq:radial-wave-equation} and (II A2.8)\ref{eq:wave-func-without-abs} to give
\begin{equation}\label{eq:schroding-plugin}
\begin{aligned}
& \int_0^a \left( \frac{d^2u_{l}}{dr^2}w_{\lambda l} - u_{l}\frac{d^2w_{\lambda l}}{dr^2} \right) dr \\
& = \int_0^a\left( \left[k^2-\frac{2m}{\hbar^2}(V+iW)-\frac{l(l+1)}{r^2}\right]u_lw_{\lambda l} + u_l\left[k_\lambda^2-\frac{2m}{\hbar^2}V-\frac{l(l+1)}{r^2}\right]w_{\lambda l} \right)dr \\
& = (k_\lambda^2 - k^2) \int_0^a u_lw_{\lambda l}dr + \frac{2m}{\hbar}\int_0^aWu_lw_{\lambda l}dr.
\end{aligned}
\end{equation}
\noindent
Defining $\overline{W}_{\lambda l}$ as
\begin{equation}
\overline{W}_{\lambda l} = \frac{ \int_0^a Wu_lw_{\lambda l}dr }{ \int_0^a u_lw_{\lambda l}dr }
\end{equation}
\noindent
permits rewriting Eq. (II A2.18)\ref{eq:schroding-plugin} in the form
\begin{equation}\label{eq:schroding-plugin-reduced}
\int_0^a \left( \frac{d^2u_{l}}{dr^2}w_{\lambda l} - u_{l}\frac{d^2w_{\lambda l}}{dr^2} \right) dr = \left( k_\lambda^2 - k^2 + i\frac{2m}{\hbar^2}\overline{W}_{\lambda l} \right)\int_0^a u_lw_{\lambda l}dr.
\end{equation}
\noindent
Integrating the left-hand side of this equation gives
\begin{equation}\label{eq:schroding-left}
\begin{aligned}
\int_0^a \left( \frac{d^2u_{l}}{dr^2}w_{\lambda l} - u_{l}\frac{d^2w_{\lambda l}}{dr^2} \right) dr & = \left[ \frac{du_{l}}{dr}w_{\lambda l} - u_{l}\frac{dw_{\lambda l}}{dr} \right]_0^a = \left[ \frac{du_{l}}{dr}w_{\lambda l} - u_{l}\frac{dw_{\lambda l}}{dr} \right]_{r=a} \\
& = \left[ \frac{du_{l}}{dr}w_{\lambda l} - u_{l}\frac{B_l}{a}w_{\lambda l} \right]_{r=a} = \left[ a\frac{du_{l}}{dr} - u_{l}B_l \right]_{r=a} \frac{w_{\lambda l(a)}}{a},
\end{aligned}
\end{equation}
\noindent
in which we have again made use of the boundary condition of Eq. (II A2.9)\ref{eq:schrod-bound-conditions}. Integrating the right-hand side of Eq. (II A2.20)\ref{eq:schroding-plugin-reduced} by applying Eq. (II A2.16)\ref{eq:c-lambda-l} gives
\begin{equation}\label{eq:schroding-right}
\left( k_\lambda^2 - k^2 + i\frac{2m}{\hbar^2}\overline{W}_{\lambda l} \right)\int_0^a u_lw_{\lambda l}dr = \left( k_\lambda^2 - k^2 + i\frac{2m}{\hbar^2}\overline{W}_{\lambda l} \right)c_{\lambda l}.
\end{equation}
\noindent
Equating Eqs. (II A2.21)\ref{eq:schroding-left} and (II A2.22)\ref{eq:schroding-right} therefore gives
\begin{equation}
\begin{aligned}
\left[ a\frac{du_l}{dr}-u_lB_l \right]_{r=a}\frac{w_{\lambda l}}{a} & = \left( k_\lambda^2 - k^2 + i\frac{2m}{\hbar^2}\overline{W}_{\lambda l} \right)c_{\lambda l}, \\
\left[ a\frac{du_l}{dr}-u_lB_l \right]_{r=a}\frac{w_{\lambda l}}{a} & = \left( E_\lambda - E + i\overline{W}_{\lambda l} \right)\frac{2mc_{\lambda l}}{\hbar^2},
\end{aligned}
\end{equation}
\noindent
or
\begin{equation}
c_{\lambda l} = \frac{\hbar^2w_{\lambda l}(a)}{2ma\left(E_\lambda-E-i\overline{W}_{\lambda l}\right)} \left[a\frac{du_l}{dr}-u_lB_l\right]_{r=a}.
\end{equation}
\noindent
Inserting this into Eq. (II A2.15)\ref{eq:internal-wave-func} gives
\begin{equation}
u_l(r) = \sum_{\lambda} w_{\lambda l}(r)\frac{\hbar^2w_{\lambda l}(a)}{2ma\left(E_\lambda-E-i\overline{W}_{\lambda l}\right)} \left[a\frac{du_l}{dr}-u_lB_l\right]_{r=a},
\end{equation}
\noindent
which when evaluated at $r = a$ , becomes
\begin{equation}
u_l(a) = \sum_{\lambda} \frac{\hbar^2w_{\lambda l}^2(a)}{2ma\left(E_\lambda-E-i\overline{W}_{\lambda l}\right)} \left[a\frac{du_l}{dr}-u_lB_l\right]_{r=a},
\end{equation}
\noindent
Rearranging, this becomes
\begin{equation}\label{eq:u-l-eval-at-bound}
\begin{aligned}
u_l(a) & = \left[a\frac{du_l}{dr}-u_lB_l\right]_{r=a} \sum_{\lambda} \frac{\hbar^2w_{\lambda l}^2(a)/2ma}{\left(E_\lambda-E-i\overline{W}_{\lambda l}\right)}
& = \left[a\frac{du_l}{dr}-u_lB_l\right]_{r=a} \sum_{\lambda} \frac{\gamma_{\lambda l}^2}{\left(E_\lambda-E-i\overline{\Gamma}_{\lambda l}/2\right)},
\end{aligned}
\end{equation}
\noindent
in which the decay amplitude $\gamma_{\lambda l}$ is defined as
\begin{equation}
\gamma_{\lambda l} \equiv \sqrt{\frac{\hbar^2w_{\lambda l}^2(a)}{2ma}}
\end{equation}
\noindent
and the absorption width $\Gamma_{\lambda l}$ as
% NOTE: I'm changing from Nancy's formula here with overline{W}
\begin{equation}
\Gamma_{\lambda l} \equiv 2 \overline{W}_{\lambda l} .
\end{equation}
\noindent
If we then define the R-function as
\begin{equation}
R_l = \sum_\lambda \frac{\gamma_{\lambda l}^2}{\left(E_\lambda-E-i\overline{\Gamma}_{\lambda l}/2\right)},
\end{equation}
\noindent
then Eq. (II A2.27)\ref{eq:u-l-eval-at-bound} can be written in the form
\begin{equation}\label{eq:u-l-r-matrix}
u_l = \left(a\frac{du_l}{dr}-u_lB_l\right)R_l,
\end{equation}
\noindent
in which everything is evaluated at the matching radius a.
\noindent
\textbf{Scattering matrix in terms of R-matrix (neutrons only)}
Equation (II A2.31)\ref{eq:u-l-r-matrix} can be converted into the usual R-matrix formulae by inserting Eq. (II A2.5)\ref{eq:in-out-wave-funcs},
\begin{equation}
u_l = I_l - U_l O_l,
\end{equation}
\noindent
yielding
\begin{equation}\label{eq:u-l-expanded}
I_l - U_l O_l = \left[a\left(\frac{dI_l}{dr}-U_l\frac{dO_l}{dr}\right)-B_l(I_l-U_lO_l)\right]R_l,
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
\noindent
in which everything is again evaluated at the matching radius $a$. Solving Eq. (II A2.33)\ref{eq:u-l-expanded} for $U$ gives
\begin{equation}\label{eq:u-l-i-l-matching}
U_l\left[-O_l+R_l\left(a\frac{dO_l}{dr}-B_lO_l\right)\right] = I_l-R_l\left(a\frac{dI_l}{dr}-B_lI_l\right),
\end{equation}
\noindent
or
\begin{equation}
U_l = \frac{ I_l-R_l\left(a\frac{dI_l}{dr}-B_lI_l\right) }{ \left[-O_l+R_l\left(a\frac{dO_l}{dr}-B_lO_l\right)\right] } = \frac{I_l}{O_l}\frac{ 1-R_l\left(\frac{a}{I_l}\frac{dI_l}{dr}-B_l\right) }{ 1-R_l\left(\frac{a}{O_l}\frac{dO_l}{dr}-B_l\right) }.
\end{equation}
\noindent
We define $L_l$ as
\begin{equation}
L_l \equiv \frac{a}{O_l(a)}\frac{dO_l}{dr}\Bigr|_{r=a} \equiv S_l + iP_l.
\end{equation}
\noindent
For spinless particles, $I_l^* = O_l$, so that
\begin{equation}
\frac{a}{I_l(a)} \frac{dI_l}{dr}\Bigr|_{r=a} = L_l^* = S_l-iP_l
\end{equation}
\noindent
and
\begin{equation}
\frac{I_l}{O_l} = \frac{O_l^*}{O_l} = \frac{|O|e^{-i\phi}}{|O|e^{i\phi}} = e^{-2i\phi}.
\end{equation}
\noindent
Therefore Eq. (II A2.34)\ref{eq:u-l-i-l-matching} becomes
\begin{equation}
U_l = e^{-2i\phi}\frac{1-R_l(L_l^*-B_l)}{1-R_l(L_l-B_l)},
\end{equation}
\noindent
which is the usual form for the scattering matrix in terms of the R-matrix in this simple case.
\noindent
\textbf{Relating the scattering matrix to the cross sections}
The relationship between the scattering matrix $U$ and the cross section $\sigma$ is also described by many authors; see, for example, [AF71]\cite{foderaro_1971}. Here we provide a summary for the simplest case. The wave function for a spinless particle far from the scattering source can be written as
\begin{equation}
\psi(r,\theta) = e^{ikz} + \frac{e^{ikr}}{r}f(\theta),
\end{equation}
\noindent
where $f$ has the form
\begin{equation}\label{f-theta}
f(\theta) = \frac{1}{2ik}\sum_l(2l+1)\left[U_l-1\right]P_l(\cos\theta).
\end{equation}
\noindent
The cross section is then given by
\begin{equation}\label{eq:angle-diff-xs-to-f-theta}
\frac{d\sigma}{d\Omega} = \left|f(\theta)\right|^2.
\end{equation}
For angle-integrated cross sections, the equation found by inserting Eq. (II A2 a.2)\ref{f-theta} into Eq. (II A2 a.3)\ref{eq:angle-diff-xs-to-f-theta} can be integrated to give
\begin{equation}
\begin{aligned}
\sigma =& \int \left[-\frac{1}{2ik}\sum_l(2l+1)\left[U_l^*-1\right]P_l(\cos\theta)\right] \\
& \times\left[\frac{1}{2ik}\sum_{l'}(2l'+1)\left[U_{l'}-1\right]P_{l'}(\cos\theta)\right]d(\cos\theta)d\phi \\
= & \frac{1}{4k^2}\sum_{ll'}(2l+1)(2l'+1)[U_l^*-1][U_{l'}-1]\int_0^{2\pi}d\phi\int_{-1}^1P_l(\cos\theta)P_{l'}(\cos\theta)d(\cos\theta) \\
= & \frac{1}{4k^2}\sum_{ll'}(2l+1)(2l'+1)[U_l^*-1][U_{l'}-1]2\pi\frac{2}{2l+1}\delta_{ll'} \\
= & \frac{\pi}{k^2}\sum_l(2l+1)|U_l-1|^2 .
\end{aligned}
\end{equation}
This is analogous to the ``standard'' scattering theory equation shown in Eq. (II A.1)\ref{eq:sigma-ccprime}.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Versions of Multilevel R-matrix Theory} \label{sec:versions-of-r-matrix}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Many representations of multilevel R-matrix theory have been developed over the years.
For a summary of the more common versions, the reader is referred to the works of Fr\"{o}hner
[FF80, FF00]\cite{frohner1980applied,jeff18_frohner}.
Four versions of R-matrix theory are available in SAMMY: the Reich-Moore
approximation (Section II.B.1)\ref{subsec:r-matrix-RM}, the single-level (SLBW) and multilevel Breit-Wigner (MLBW) approximations (Section II.B.3)\ref{subsec:r-matrix-BW}, and a variant on the Reich Moore which mimics the full R-matrix (Section II.B.2)\ref{subsec:r-matrix-full}. An option to include a direct capture component is also provided (Section II.B.4)\ref{subsec:direct-capture}.
The Reich-Moore approximation is the preferred method for most modern evaluations; it
is the default formalism for SAMMY runs. Fr\"{o}hner, in fact, suggests that the Reich-Moore approximation is universally applicable to all cases: ``Experience has shown that with this approximation [Reich Moore] all resonance cross section data can be described in detail, in the windows as well as in the peaks, even the weirdest multilevel interference patterns...It works equally well for light, medium-mass and heavy nuclei, fissile and nonfissile.'' \cite{jeff18_frohner}(pg. 60). For most purposes, Reich Moore is indeed indistinguishable from the exact formulation. Notable exceptions are situations where interference effects exist between capture and other channels. For those cases, small modifications to the SAMMY input will permit the user to mimic the effect of the non-approximated R-matrix; see Section II.B.2\ref{subsec:r-matrix-full} for details. Occasionally it is not possible to properly describe a cross section within the confines of R-matrix theory, because the reaction includes a direct component. SAMMY has provisions for the user to provide a numerical description of this component; see Section II.B.4\ref{subsec:direct-capture} for details.
Also available within SAMMY are both the SLBW and the MLBW formulations (Section II.B.3)\ref{subsec:r-matrix-BW}; these are included for the sake of completeness, for comparison purposes, and because many of the evaluations in the nuclear data files were performed with Breit-Wigner formulae. However, it is strongly recommended that only Reich Moore be used for new evaluations, for several reasons: MLBW is often inadequate; SLBW is almost always inadequate. When it is correct, MLBW gives identical results to Reich Moore. ``Ease of Programming'' is no longer a valid excuse for using MLBW, since the programming has already been accomplished. Similarly, a slow computer is no longer a legitimate excuse, since modern computers can readily handle the more rigorous formulae.
Finally, it should be noted that SAMMY's implementation of MLBW does not correspond to the usual definition of MLBW. Instead, SAMMY uses the ENDF [ENDF-102]\cite{endf8} convention in which only the elastic cross section is truly multilevel, and all other types of cross section are single level.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Reich-Moore Approximation to Multilevel R-Matrix Theory} \label{subsec:r-matrix-RM}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Simulation of Full R-Matrix} \label{subsec:r-matrix-full}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Breit-Wigner Approximations} \label{subsec:r-matrix-BW}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Direct Capture Component} \label{subsec:direct-capture}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%