\section{Introduction}

Nuclear reactions offer an excellent probe into the properties of nuclei, particularly for
short-lived rare isotopes. Nuclear reactions are extremely proliferous:  they provide access
to structure properties that improve our knowledge of the underlying nucleon-nucleon (NN) force, they are used to populate excited states of interest, and through them we can extract astrophysical rates that cannot otherwise be directly measured.  However, a common concern when using nuclear reactions to extract nuclear properties pertains to the simplifications in describing the dynamics and consequent model dependence of the extracted properties.
For this reason, it is a priority that our field develops new methods for reactions that are not
limited by unnecessary  approximations.
In short, the goal is for a theory that includes  an exact treatment of the dynamics for the relevant degrees of freedom and
incorporates the  relevant reaction channels in a consistent framework and on an equal footing. In
addition, this reaction theory should be applicable across the nuclear chart and for a wide
range of energy regimes, so that it is not plagued by irreconcilable systematic
differences~\cite{Atar:2018dhg}.
This is the context of the present study. We focus first on deuteron induced
reactions on the $\alpha$-particle, but note that  the framework is readily extendable to heavier projectiles.

While nuclear reactions are in a fundamental way many-body problems, in direct reactions, when
only a few degrees of freedom play a role, the problem is often reduced to a few-body problem.
As early as the pioneering work by Shanley~\cite{Shanley:1969zza},  three-body approaches  have
been successfully used to simplify the nuclear reaction problem and allow for an exact
treatment of the few-body dynamics. In Refs.~\cite{Shanley:1969zza,Miyagawa:1985js},
three-body solutions of the
Alt-Grassberger-Sandhas (AGS) equations~\cite{ags1967}  for the scattering of deuterons off an
$\alpha$-particle were obtained using  rank-1 separable two-body forces. Both elastic
scattering observables and the total breakup cross sections, $n+p+\alpha$,  were computed and
compared to experiment. Despite its simplicity, the model provided a fair description of the data.

In the last decade, the AGS equations were  applied in the field of nuclear reactions by
Deltuva and Fonseca~\cite{deltuva:06b,deltuva2009a} without the need to employ separable
interactions, given the advances in computational power. The greatest challenge however,
when moving from applications in the few-nucleon sector to reactions with heavy ions concerns
the Coulomb force. Deltuva and Fonseca relied on screening the Coulomb interaction and subsequent renormalization to ensure
the compactness of the Faddeev equations. Over the last decade, there have been many
applications of this approach on a variety of reactions, such as deuteron stripping and
pickup reactions involving halo nuclei, with latest developments
including also excitations of the nuclear core \cite{deltuva:17b}. Albeit the success of the method,
incorporating the Coulomb force still poses unwanted limitations: when it becomes too strong
(i.e. for heavier nuclei and for lower scattering energies), the Coulomb screening radius
needed for an accurate description of the reaction increases to a point that renders the screening method ineffective.
In practice this means that the methods in~\cite{deltuva:06b,deltuva2009a,deltuva:17b} have  been applied so far
to targets with mass $A \le 58$.

An alternative to introducing Coulomb screening in the AGS equations is to cast the AGS
equations in the Coulomb basis.
In Ref.~\cite{akram2012},
Mukhamedzhanov and collaborators derived the AGS equations in the Coulomb basis,  indicating
also the necessity for employing a separable expansion of the forces in the different
subsystems. That work is focused on the formalism and includes no numerical applications. It thus remains to be proven that this alternative method indeed can provide precise solutions for deuteron induced reactions on heavy nuclei.

Our work represents an important stepping stone for using Faddeev techniques for nuclear
reactions with heavy nuclei. Starting from the non-relativistic AGS equations~\cite{ags1967},
we  make use of the formulation derived by Lovelace~\cite{lovelace1964},  and introduce
pair-wise separable interactions as previously
developed~\cite{hlophe2013,hlophe2014,hlophe2017}. For the treatment of the
singularities in the
free three-body propagator above the breakup threshold, we implement a procedure proposed
in~\cite{Witala:2008my,Elster:2008hn}, which allows one to cast the so-called moving singularities
(see e.g.~\cite{Gloecklebook})
into singularities depending only on one variable, amenable to regular subtraction techniques.
In this work, we do not explicitly include target excitation in our formulation. However,
the formalism readily extends to problems in which the two-body subsystems can be described as a coupled-channel problem involving the excitation of one of the bodies~\cite{hlophe2017est}.

In few-nucleon physics, the use of separable interactions is widespread and has been proven to
be accurate in describing neutron-deuteron processes below the pion
threshold~\cite{nemoto1998,Cornelius:1990zz}.
It is important that these benchmarks be performed for nuclear reactions, where the interaction
has larger complexity. Recently, we performed a benchmark for the bound state of $^{6}$Li~\cite{hlophe2017} and showed that one can obtain 4 digit accuracy in the binding energy if
including rank-8 potentials.  In the present work, we perform a similar benchmark but now for
deuteron-alpha elastic scattering and breakup reactions.
We compare the results obtained in the new framework using separable interactions with the solutions obtained when no separable expansion is introduced~\cite{deltuva:06b}.
The current work  includes only  nuclear interactions to allow for a careful inspection of the treatment of the short range aspects of the problem. We will focus on the inclusion of Coulomb in a subsequent study.

The paper is organized in the following way.
In Sec.~II, a brief summary of the theory is provided, highlighting essential ingredients for
our calculations. Section III introduces the two-body input to our calculations and presents
deuteron-alpha scattering observables for elastic as well as breakup scattering together with
essential convergence tests. We conclude in Sec. IV. A detailed description of our treatment
of the three-body breakup singularities and the calculation of the kinematical S-curve for
three particles with different masses is given in the appendices.

\section{Formal Considerations }

For the description of the scattering of a deuteron off an alpha particle, we employ the Alt,
Grassberger, and Sandhas (AGS) equations~\cite{AGS}, which are Faddeev-type integral equations in
momentum space for three-particle transition operators
\begin{equation}
 U^{ij} = \bar\delta_{ij}\;G_0^{-1} +\sum\limits_{k} \bar\delta_{ik} t_k G_0 \; U^{kj}.
\label{eq:2.1}
\end{equation}
Here $\bar\delta_{ij}=1-\delta_{ij}$ is the anti-Kronecker delta, $G_0(E)=(E+i0 -H_0)^{-1}$ is the free resolvent at the available three-particle center-of-mass energy $E$, and $H_0$ is the free three particle Hamiltonian. The three particles with masses $m_i$, $m_j$, and $m_k$ and spins $s_i$, $s_j$, and $s_k$ interact via pairwise forces $v^i \equiv v_{jk}$ ($i,j,k=1,2,3$ and cyclic permutations thereof) so that the operator $t_i(E) = v^i +v^i G_0(E) t_i(E)$ describes the two-body $t$~matrix in the subsystem $i$. The AGS transition amplitudes $U^{ij}$ are represented in their natural Jacobi coordinates $(p_i,q_i)$, where $p_i$ is the relative momentum of the $i^{th}$ pair and $q_i$ the momentum of the $i^{th}$ particle (also called the `spectator') relative to the pair. For solving the AGS equations, we choose a momentum space basis which depends on the magnitude of the momenta and angular momentum eigenstates. To proceed, we define the total spin $S_i=s_j+s_k$, relative orbital angular momentum $l_i$, and the total angular momentum by $J_i=l_i+S_i$ for the $i^{th}$ pair.
The orbital angular momenta and spin of the spectator are represented by $\lambda_i$ and $s_i$ so
that the corresponding total angular momentum is ${\cal J}_i=\lambda_i+s_i$. The third component
of the spectator spin is denoted by $m_{s_i}$. The total angular momentum of the three-particle system is represented by $J$ and the third component by $M_J$. The states of conserved total angular momentum are thus given
as
\begin{equation}
|{p}_i{q}_i \alpha_{i} \rangle = |p_iq_i\;\big(l_i(s_is_j)S_i)J_{i} (\lambda_i s_i){\cal J}_i\big)\;JM_J
\rangle,
\label{eq:2.2}
\end{equation}
and are normalized as
\begin{equation}
\langle {p}'_i{q}'_i \alpha_{i} ' | {p}_i{q}_i \alpha_{i} \rangle=
\frac{\delta(p'_i-p_i)}{p'_ip_i} \frac{\delta(q'_i-q_i)}{q'_iq_i} \delta_{\alpha_{i} '\alpha_{i}}.
\label{eq:2.3}
\end{equation}
The matrix elements of an AGS operator for the transition from an arrangement channel $j$ to $i$ is given by the expectation value $\langle \Phi_i|U^{ij}|\Phi_j \big\rangle$, where
\begin{eqnarray}
 |\Phi_i\rangle\equiv\sum\limits_{l_iS_i}\big|\phi_{{l_iS_i}}^{J_i}\;(l_iS_i)
J_{i}M_{J_i}\big\rangle|{\bf q}_i\;s_im_{s_i}\rangle,
\label{eq:2.4}
\end{eqnarray}
is the asymptotic state in arrangement channel $i$ with $|\phi_{l_iS_i}^{J_i}\big\rangle$ being the two-body bound state  wavefunction.

Two approaches for solving Eqs.~(\ref{eq:2.1}) are adopted.
We shall refer to them as
the `separable' and the `non-separable' approach.  The former involves a separable
expansion of the two-body $t$~matrix while the latter does not. The non-separable method
is taken over from Ref.~\cite{deltuva:06b} and calculates fully off-shell $t$ matrices on the momentum grid, with subsequent interpolation using global spline functions whenever needed.

For the separable expansion method, the subsystem $t$~matrices take the well known form (see e.g.~\cite{Ernst:1973zzb})
\begin{equation}
t_{i}^{\alpha_{i} \alpha_{i}'} ( E_{q_i}) = \sum_{mn} |h^{\alpha_{i}}_m\rangle
\tau_{mn}^{\alpha_{i} \alpha_{i}'}(E_{q_i})
\langle h^{\alpha_{i}'}_n| \ ,
\label{eq:2.5}
\end{equation}
where $|h^{\alpha_{i}}_m\rangle$ are the so-called form factors and $E_{q_i}\equiv E-q_i^2/2M_i$ represents the available two-body energy. Here $M_k$ is the reduced mass of the pair and the spectator. The indices ${m,n}$ represent the rank
of the separable potential, and $i$ stands for the arrangement channel.
If the potential $v^{i}$ supports a bound state, the corresponding wave-function has the form
\begin{equation}
|\phi^i_{\alpha_i}\rangle=\sum\limits_{m}G_0(\epsilon_{\alpha_i})\;|h_{m}^{\alpha_i}\rangle\;c_{m\alpha_i},
\label{eq:2.6}
\end{equation}
where the two-body bound state energy is related to the on-shell
spectator momentum $q_{i0}$ by
\begin{equation}
\epsilon_{\alpha_i} =E_{q_{i0}}= E -\frac{q_{i0}^2}{2M_i}.
\label{eq:2.7}
\end{equation}
The constants $c_{m\alpha_i}$ are determined by substituting Eq.~(\ref{eq:2.6}) into a bound
state Lippmann-Schwinger equation. The partial wave matrix element for a transition from the bound state $|\phi^i_{\alpha_i}\rangle$ to $|\phi^j_{\alpha_j}\rangle$ becomes
\begin{equation}
\langle\phi^i_{\alpha_i}q_i\alpha_i|U^{ij}|\phi^j_{\alpha_j} q_j \alpha_j\rangle
=\sum\limits_{mn} c_{m\alpha_i}^*c_{n\alpha_j}
\langle h_m^{\alpha_i}\;q_i \alpha_i|G_0(\epsilon_{\alpha_i})U_{ij}
G_0(\epsilon_{\alpha_j})|h_n^{\alpha_j}q_j \alpha_j\rangle.
\label{eq:2.8}
\end{equation}
To proceed, one defines the effective two-body AGS transition amplitudes~\cite{AGS}
\begin{equation}
 X_{m\alpha_i,n\alpha_j}^{ij}(q_i,q_j;E)\equiv
\langle h_m^{\alpha_i}\;q_i \alpha_i|G_0(E)U^{ij}
G_0(E)|h_n^{\alpha_j}q_j \alpha_j\rangle ,
\label{eq:2.9}
\end{equation}
whose on-shell values coincide with those of Lovelace~\cite{Lovelace:1964} and
appear on the right-hand side of Eq.~(\ref{eq:2.8}). From Eqs.~(\ref{eq:2.1}), one obtains

\begin{eqnarray}
\langle h_m^{\alpha_i}q_i \alpha_i|G_0(E)U^{ij}G_0(E)|h_n^{\alpha_j}q_j\alpha_j\rangle &=&
\bar\delta_{ij} \;\langle h_m^{\alpha_i}q_i \alpha_i|G_0(E)|h_n^{\alpha_j}q_j\alpha_j\rangle\cr\cr
&+& \sum_k \sum_{\alpha_k\alpha_k'} \bar\delta_{ik} \int dq_k q_k^2\;
 \langle h_m^{\alpha_i}q_i \alpha_i |G_0(E) | q_k \alpha_k\rangle\langle q_k \alpha_k|
t_{k}(E_{q_k})|q_k\alpha_k'\rangle \cr\cr
&& \times \langle q_k\alpha_k' | G_0(E) U^{kj}G_0(E)| h_n^{\alpha_j}q_i\alpha_j\rangle
\label{eq:2.10},
\end{eqnarray}
and substituting the separable two-body $t$~matrix yields

\begin{eqnarray}
\langle h_m^{\alpha_i}q_i \alpha_i|G_0(E)U^{ij}G_0(E)|h_n^{\alpha_j}q_j\alpha_j\rangle &=&
\bar\delta_{ij} \;\langle h_m^{\alpha_i}q_i \alpha_i|G_0(E)|h_n^{\alpha_j}q_j\alpha_j\rangle\cr\cr
&+& \sum_k \sum_{\alpha_k\alpha_k'}\sum_{m'n'} \bar\delta_{ik} \int dq_k q_k^2\;
 \langle h_m^{\alpha_i}q_i \alpha_i |G_0(E) |h_{m'}^{\alpha_k} q_k \alpha_k\rangle\;
\tau^{\alpha_k\alpha_k'}_{m'n'}(E_{q_k})\; \cr\cr
&& \times \langle q_k\alpha_k' h_{n'}^{\alpha_k'} | G_0(E) U^{kj}G_0(E)| h_n^{\alpha_j}q_j\alpha_j\rangle.
\label{eq:2.11}
\end{eqnarray}
Defining the effective two-body `transition potentials'
\begin{align}
Z_{m\alpha_i,n\alpha_j}^{ij}(q_i,q_j,E) = \bar\delta_{ij} \;\langle h_m^{\alpha_i}q_i \alpha_i|G_0(E)|h_n^{\alpha_j}q_j\alpha_j\rangle,
\label{eq:2.12}
\end{align}
and using Eq.~(\ref{eq:2.9}), one can express Eqs.~(\ref{eq:2.11}) in a condensed form
\begin{align}
X_{m\alpha_i,n\alpha_j}^{ij}(q_i,q_j;E) & = Z_{m\alpha_i,n\alpha_j}^{ij}(q_i,q_j,E)
+ \sum_{k}\sum_{\alpha_k\alpha_{k}'}\sum_{m'n'} \int dq_k q_k^2\; Z_{m\alpha_i,m'\alpha_k}^{ik}(q_i,q_k;E) \nonumber \\
&\times \tau_{m'n'}^{\alpha_k\alpha_k'}(E_{q_k}) X_{n'\alpha_k',n\alpha_j}^{kj}(q_k,q_j;E).
\label{eq:2.13}
\end{align}
These equations are solved using iterative Lanzcos-type techniques~\cite{Saad:2003}.
The kernel contains two types of singularities, namely, the bound state pole and the
three-body breakup singularity. One one hand, the former constitutes a simple pole and is removed using standard  subtraction techniques. On the other hand, the three-body breakup pole has a complex structure
and its treatment is consequently more involved. A detailed discussion is provided in Appendix~\ref{sing}.

 To evaluate the transition amplitude
for a breakup process, we first define the on-shell momentum of the pair $p_{i0}=\sqrt{2\mu_i(E-q_i^2/2M_i)}$, where $\mu_i$ is the reduced mass of the pair. The on-shell breakup amplitude is then expressed in terms of the amplitudes for elastic scattering and rearrangement processes using Eq.~(\ref{eq:2.1}) leading to
\begin{eqnarray}
\langle   q_i p_{i0} \alpha_i|U^{0j}|\Phi_j \big\rangle=\sum\limits_{k=1}^3 \langle  p_{0i}\; q_i\alpha_i|t_k\;G_0(E)\;U^{kj}|\Phi_j \big\rangle.
\label{eq:2.14}
\end{eqnarray}
For separable two-body $t$ matrices, one obtains
\begin{eqnarray}
\langle   q_i p_{i0} \alpha_i|U^{0j}|\Phi_j \big\rangle=\sum\limits_{k=1}^3 \sum_{\alpha_k\alpha_k'\alpha_j}
\sum_{mn'n}\;\int dp_k p_k^2\;dq_k q_k^2\;\langle q_i p_{i0} \alpha_i|q_k p_k \alpha_k\rangle\; h_{m}^{\alpha_k}(p_k)\;\tau^{\alpha_k\alpha_k'}_{mn'}(E_{q_k})\;X_{n'\alpha_k',n\alpha_j}^{kj}(q_k,q_j;E)\;c_{n\alpha_j}.~~~~~~
\label{eq:2.15}
\end{eqnarray}
The matrix elements $\langle q_i p_{i0} \alpha_i|q_k p_k \alpha_k\rangle$ describe a transformation between two Jacobi coordinates for $i\ne k$ and are evaluated as described in
Ref.~\cite{Hlophe:2017bkd}. If $i=k$, these matrix elements reduce to the $\delta$-functions of Eq.~(\ref{eq:2.3}). It is thus seen that, once the effective two-body amplitudes $X^{ij}_{mn}$ have been determined, the transition amplitudes for elastic scattering, as well as rearrangement and breakup processes can be readily computed.

\section{Results and Discussion}
\label{discuss}

To demonstrate the accuracy of solving the Faddeev-AGS equations for deuteron-alpha scattering based on a
separable expansion of the forces in the two-body subsystems, the convergence of the
expansion must be tested and finally the converged calculation must be compared to
a numerically converged calculation based on the non-separable version of the same forces.
For this, we first define the three-body Hamiltonian with the forces in the different
two-body subsystems, namely the neutron-proton ($np$) force and the effective
interaction between a neutron or
proton and the alpha-particle, i.e. the $n/p$-$\alpha$ forces.
Then, elastic as well as breakup observables for $d$-$\alpha$ scattering are calculated, and
their convergence with respect to the basis (rank) of the separable expansion is explored.
Finally, the well-converged separable calculations are compared to the corresponding ones
obtained with the non-separable forces in the subsystems.

\subsection{Two-body input}
\label{input}

For the $np$ force, the high precision CD-Bonn~\cite{Machleidt:2000ge} potential is adopted.
As an effective interaction in the  $n/p$-$\alpha$ system, we employ the
Bang potential as given in Ref.~\cite{Thompson2000}
 consisting of an attractive  central Woods-Saxon  and spin-orbit terms.
The two-body model space is restricted to $l_{i}\le 2$ for both the $np$ and $n/p$-$\alpha$ systems. The Coulomb repulsion in the $p$-$\alpha$ system is omitted.
The Bang potential supports a bound state in the $n/p$-$\alpha$ S-wave channel,
which is removed with the projection technique described in~\cite{Hlophe:2017bkd}.

For constructing  separable representations for the two-body $t$ matrices which enter the AGS
equations, we employ the method suggested by Ernst, Shakin, and Thaler~\cite{Ernst:1973zzb} (EST).
The advantage of this scheme is
that the  $t$ matrices calculated at specific energies with a given potential are taken as form factors of the
separable expansion. While in~\cite{Ernst:1973zzb} half-shell $t$ matrices are used, an extension to
using off-shell $t$ matrices is straightforward~\cite{Ernst:1974zza}.
The details on the explicit construction of the separable representations employed here are laid out
in~\cite{Hlophe:2013xca}. The application to the three-body calculation of the $^6$Li bound state is presented in~\cite{Hlophe:2017bkd}. For the convenience of the reader, as well as to establish notation, essential ingredients are briefly repeated.

The EST approach applied in momentum space requires solving a two-body Lippmann-Schwinger (LS) equation
at a specific scattering energy $E_n$ with a given
real (or complex) potential, leading to a $t$ matrix $t(p,p_n;E_n)$. The on-shell $t$ matrix corresponds to momentum $p_n=\sqrt{2\mu E_n}$. For all other values of $p_n$, the $t$ matrix is fully off-shell. We use  $E_n$ (support energies) and $p_n$  (support momenta) independently to construct the separable interactions.
Thus any specific solution of the LS equation is characterized by the pair ($E_n$, $p_n$), which we call
an EST support point. Naturally, when $E_n$ corresponds to a bound state energy in the two-body system, the
$t$ matrix is always fully off-shell.
In Ref.~\cite{Hlophe:2017bkd}, we calculated the three-body binding energy of $^6$Li  and found that by
choosing $p_n$ independently from the bound state energy, we could achieve better
accuracy.

To ensure that the separable expansion converges, it is imperative to calculate observables using
successively increasing ranks. For this reason, we define separable representations of the CD-Bonn and
Bang potentials with ranks ranging from 3 to 7. Table~\ref{table1} lists the EST support points used to
construct the separable representations of the CD-Bonn potential, while  Table~\ref{table2} shows those
for the Bang potential. The separable representations provide a good description of the two-body $t$ matrix over a range of relative two-body
energies $E_{2}$ which, in the context of solving the Faddeev-AGS equations, corresponds to $-\infty < E_{2} \le E$  ($E$ being the three-body energy in the c.m. frame).

\subsection{Deuteron-alpha scattering observables}

The validity of the separable expansion depends on the beam energy for the reactions, thus
we choose three deuteron beam energies, $E_d=10,\;20,$ and$\;50$~MeV, for the benchmark
calculations of the separable vs. non-separable solutions of the AGS equations. These energies cover the typical range for experiments of (d,p) reactions.

\subsubsection{Integrated cross sections}
\label{obsv:intg}

Integrated cross sections are an important test of the calibration of our methods.  We aim to achieve a precision of $\approx 1$\% on this observables which is well below the typical experimental errors and the uncertainties associated with the two-body interactions. Tables~\ref{table3} and \ref{table4} show the integrated cross sections computed using the separable potentials given in Tables~\ref{table1} and~\ref{table2} at 10 and 20~MeV incident
deuteron energy. We show results for separable potentials constructed only with support momenta obeying the constraint $p_n=\sqrt{2\mu|E_n|}$ (NN-EST3-1, NN-EST4-1, NN-EST5-1, NN-EST6-1, NN-EST7-1, NA-EST3-1, NA-EST4-1, NA-EST5-1, NA-EST6-1, and NA-EST7-1) as in~\cite{Cornelius:1990zz}. We also show results for separable potentials for which the support momenta $p_n$ are independent from $E_n$ (NN-EST6-2, NN-EST6-3, NN-EST7-2, NA-EST6-2, NA-EST6-3,
and NA-EST7-2). These results are benchmarked against  the integrated cross sections obtained using the non-separable approach, given in
the bottom row. The first point to make is that a similar convergence rate is observed for both elastic and breakup cross sections. Moreover, the convergence pattern is similar for both the
10 and 20~MeV deuteron beam energies.
The second point worth making is that, in this case, the inclusion of off-shell momenta does not represent an improvement of the restricted $p_n=\sqrt{2\mu|E_n|}$ basis, contrary to our observation for bound states. Thirdly, a rank-6 interaction already provides the desired level of precision.
Finally, and most importantly, the results agree with the cross sections obtained from the non-separable calculation  to $\approx 1$\%.

 It is interesting to contrast our findings with those
 of Ref.~\cite{Cornelius:1990zz} for $nd$ scattering at 10~MeV incident neutron energy. In that work the authors demonstrated
 that an EST rank-3 potential could reproduce the integrated cross sections calculated
 using the non-separable method to about 2~$\%$. From Tables~\ref{table3} and~\ref{table4},
 we see that the discrepancy between results calculated with the separable and non-separable
approach is already better than $\sim 1.5\%$ for the rank-4 potentials, indicating that our findings are
consistent with those of Ref.~\cite{Cornelius:1990zz}.

\subsubsection{Elastic scattering: deuteron angular distributions}
\label{obsv:el}

Next, we consider the differential cross sections for elastic deuteron-alpha
scattering for three deuteron beam energies $E_d$.
In Fig.~\ref{fig1}, the differential cross sections for elastic $d+\alpha$ scattering as a function of
the center-of-mass (c.m.) angle $\theta_{c.m.}$ are shown for
three different incident deuteron energies, $E_d=$~10, 20~MeV, and 50~MeV. The solid lines indicate the
angular distributions evaluated using the non-separable approach. The calculations
 obtained using rank-3, rank-4, and rank-5 separable potentials are indicated by the dash-dot-dotted,
dash-dotted, and dashed lines. All four lines are nearly indistinguishable which demonstrates that the
separable expansion is not only well-converged, but also  yields deuteron angular distributions that are in excellent agreement with those obtained using the non-separable approach. The relative differences between the two results remain at $\approx 1\;$\% throughout the angular
range. This is well below the usual experimental uncertainties for this observable.
It is worth noting that the separable expansion
method is not limited to the low beam energies mentioned above. Increasing the beam energy by tens of
MeV does not introduce any principal technical difficulties except that, obviously, the dependence on the rank of the
separable expansion needs to be reevaluated.

\subsubsection{Deuteron breakup: fivefold differential cross section}
\label{obsv:bu}

 For an even more stringent test of the separable expansion method, we explore the
 convergence of the five-fold breakup differential cross section with respect to the rank
 of the separable expansion. Typically, one proceeds by specifying configurations
 defined by the outgoing proton and alpha particle angles $(\theta_p,\;\phi_p)$ and  $(\theta_\alpha,\;\phi_\alpha)$. The corresponding S-curves are then evaluated
 using energy and momentum conservation as described in Appendix~\ref{scurve}.
 Each configuration is given in the format  $(\theta_\alpha,\;\phi_\alpha;\;\theta_p,\;\phi_p)$, with the angles given in degrees. We consider two configurations, one of which corresponds to the final state interaction (FSI). The FSI configurations are defined such that
$E_{np}\approx 0$ for a single value of the arclength $S$.
  Here $E_{np}$ is the relative energy between the outgoing neutron and proton.

Figure~\ref{fig2} displays the five-fold breakup differential cross section on the S-curve for
$E_d=10$~MeV for two different configurations of the $\alpha+n+p$ system. Panel (a) shows results for the configuration $(25.6^{\circ},\;0^{\circ};\;63.6^{\circ},\;180^{\circ})$ while panel (b) depicts results for the FSI configuration $(31.4^{\circ},\;0^{\circ};\;5.1^{\circ},\;180^{\circ})$.
The cross sections computed with the non-separable method are indicated by the
solid line. Results evaluated using the separable expansion method are depicted
by the dash-dot-dotted, dash-dotted, and dashed lines for the rank-3, rank-4, and rank-5
potentials of Tables~\ref{table1} and \ref{table2}. The FSI point is indicated by the filled square.
The deviation of the dash-dot-dotted line from the other curves in both panels
demonstrates that the rank-3 potential is inadequate for this observable. However, the separable expansion is still rapidly converging so that the rank-4 and rank-5 curves are virtually indistinguishable. We also observe that the converged results from the separable approach are in excellent agreement with those calculated using the non-separable method.

Next, we consider the five-fold breakup differential cross section at 20~MeV incident deuteron energy. Figure~\ref{fig3} is the same as Figure~\ref{fig2} but for
$E_d=20$~MeV. The configurations depicted in panels (a) and (b) are $(29.0^{\circ},\;0^{\circ};\;22.5^{\circ},\;180^{\circ})$ and $(25.6^{\circ},\;0^{\circ};\;1.7^{\circ},\;180^{\circ})$, where
the latter corresponds to FSI. We note that the convergence pattern
is similar to that of $E_d=10$~MeV. Also, the rank-4 and rank-5 results are in very good agreement with those obtained via the non-separable approach.

In Section~\ref{obsv:el} we showed that the agreement between deuteron angular distributions calculated using the separable and non-separable methods remains
excellent when the beam energy is increased to $E_d=$~50~MeV. For completeness,
it is imperative that we also compare the five-fold breakup differential cross sections at this energy.

The comparison is shown in Fig.~\ref{fig4}.
 Panel (a) depicts the configuration $(14.0^{\circ},\;0^{\circ};\;7.2^{\circ},\;180^{\circ})$
while the off-plane configurations $(14.0^{\circ},\;0^{\circ};\;7.2^{\circ},\;120^{\circ})$ and
$(22.2^{\circ},\;0^{\circ};\;104.4^{\circ},\;100^{\circ})$
are shown in panels (b) and (c). The FSI configuration $(22.2^{\circ},\;0^{\circ};\;104.4^{\circ},\;180^{\circ})$  is illustrated in panel (d). The filled square indicates the FSI point. The solid line corresponds to results calculated using the non-separable approach while the dashed lines depicts those computed via the separable expansion method. We observe that the two methods are in excellent agreement and that the level of agreement is consistent with the one obtained for
$E_d=$~10 and 20~MeV.

\section{Summary and Outlook}

Deuteron induced nuclear reactions are a widely used probe in nuclear physics. From a theoretical standpoint, these reactions are often mapped on a three-body problem $n+p+A$. The exact solution of the three-body problem can be obtained using the momentum-space Faddeev AGS framework.
Given the long-range Coulomb force, the standard implementation of the AGS equations
for deuteron induced nuclear reactions relies on the screening and renormalization
method. However this limits the application to lighter targets and/or higher beam
energies. To circumvent this limitation, one can instead use separable interactions
and the Coulomb distorted basis as proposed in \cite{akram2012}.

In this study, we implemented an AGS framework based on EST-like  separable interactions, solved those separable AGS equations for scattering and constructed elastic and breakup observables.  For the purpose of this benchmark, no Coulomb interactions were included. We applied our method to $d-\alpha$ scattering for $E_d=10,\; 20,\; 50$ MeV, taking the same interactions as those used in a previous work where we benchmarked the use of separable interactions for three-body bound states \cite{Hlophe:2017bkd}. We find that the new method converges well. Depending on the desired observable, rank-4 or rank-6 were sufficient for obtaining 1\% precision.

For the benchmark, we also performed the AGS calculations  using the non-separable interactions.  The results obtained with the separable interactions agree very well with those using the standard non-separable AGS method. Specifically, the total cross sections and the elastic differential angular distributions calculated using these two methods agree to $1$\%.  Breakup is typically harder to calculate precisely.
For this reason,  we selected a wide variety of configurations for proton and neutron angles in order to test the method in the most extreme limits. For all cases considered, we found that the five-fold differential breakup cross sections obtained using the separable framework for rank-4 was already in very good agreement with the results using the
non-separable approach.  Moreover, no special adjustment of the separable force was needed when going from the three-body bound state calculation to the scattering application.

Consistent with the conclusion from the study of $^6$Li bound states \cite{Hlophe:2017bkd}, we here demonstrate that the separable formulation  provides a reliable method to solve the three-body scattering problem. It is now possible to follow onto the final step of the process, namely the inclusion of the Coulomb potential by solving the Faddeev-AGS equations in the Coulomb basis.