In transfer reaction, the kinematics determines the scattering angles and momenta, but it does not tell the probability of the reaction. Not all angle have same cross section, because when a particle is being transfer to or pickup from a nucleus, the momentum has to be matched. It is like a spaceship need a proper angle and speed in order to orbit around the earth, it the incident angle or speed are not good, the spaceship will just go out or crash into the earth.

Suppose our reaction is A(B,1)2, A is incident particle, B is the target, 1, 2 are the outgoing particles. The momentum transfer from A via 2 to B or 1 is

\vec{q} = \vec{p}_2 - \vec{p}_A = \vec{p}_f - \vec{p}_i

where \vec{p}_i is initial momentum and \vec{p}_f is final momentum.

Assume the reaction take place at the surface of nucleus B. The maximum angular momentum it can create is

\vec{L} = \vec{r} \times \vec{q} \rightarrow L = r q

In QM, the angular momentum must be L = \sqrt{l(l+1)} \hbar . Here we pause for unit conversion. In this blog, we usually using nuclear unit, that momentum in [MeV/c] and \hbar = 197.327 \textrm{MeV fm/c}, thus, to calculate the angular momentum, simple multiple the momentum in MeV/c and nuclear radius in fm. The radius of a nucleus is roughly r = 1.25 A^{1/3} \textrm{fm}.

So, if the momentum \sqrt{(l+1)(l+2)} \hbar > r|q| > \sqrt{l(l+1)} \hbar , the possible angular momentum is \sqrt{l(l+1)} \hbar or the l-orbit.

In the transfer reaction, the equation for \vec{p}_k and \vec{p}_i are known. It is straight forward to calculate r q / \hbar, plot against scattering angle, and we look for the value for \sqrt{l(l+1)} .

Or, we can scale the y-axis by

\displaystyle y' = \sqrt{\left(q\frac{r}{\hbar}\right)^2 + \frac{1}{4}} - \frac{1}{2}

The new y-axis is in unit of \sqrt{l(l+1)} .


We can also plot the contour of momentum matching on the p_k versus scatering angle. The momentum matching is

\displaystyle q^2 = \frac{l(l+1)\hbar^2}{r^2} = p_f^2 + p_i^2 - 2 p_f p_i \cos(\theta)

where \theta is the scattering angle of particle 1. Solve for p_f

\displaystyle p_f = p_i \cos(\theta) \pm \sqrt{ \left( \frac{l(l+1) \hbar^2}{r^2} \right) - k_i^2 \sin^2(\theta) }

We can plot the contour plot for difference l. and then overlap with the momentum and the scattering angle of the outgoing particle 2.


Screen Shot 2018-01-30 at 14.31.12.png

a correction, the y-axis should be p_2, momentum of particle 2.