Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Eikonal calculation of (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) cross sections for neutron-rich nuclei

M. Gómez-Ramos mgomez40@us.es Departamento de FAMN, Universidad de Sevilla, Apartado 1065, 41080 Sevilla, Spain. Institut für Kernphysik, Technische Universität Darmstadt, D-64289 Darmstadt, Germany
(June 4, 2024)
Abstract

In this work, we present the first, to our knowledge, theoretical description of two-proton removal reactions with proton target (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) for medium-mass nuclei at intermediate energies and present cross sections for the different bound states of the residual nucleus with two fewer protons. The description of the reaction assumes two sequential “quasifree” collisions between the target and removed protons and considers eikonal propagation in between. The formalism is applied to the reactions C12(p,3p)10BesuperscriptC12superscript𝑝3𝑝10Be{}^{12}\mathrm{C}(p,3p)^{10}\mathrm{Be}start_FLOATSUPERSCRIPT 12 end_FLOATSUPERSCRIPT roman_C ( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 10 end_POSTSUPERSCRIPT roman_Be, Mg28(p,3p)26NesuperscriptMg28superscript𝑝3𝑝26Ne{}^{28}\mathrm{Mg}(p,3p)^{26}\mathrm{Ne}start_FLOATSUPERSCRIPT 28 end_FLOATSUPERSCRIPT roman_Mg ( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 26 end_POSTSUPERSCRIPT roman_Ne and Ca54(p,3p)52ArsuperscriptCa54superscript𝑝3𝑝52Ar{}^{54}\mathrm{Ca}(p,3p)^{52}\mathrm{Ar}start_FLOATSUPERSCRIPT 54 end_FLOATSUPERSCRIPT roman_Ca ( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 52 end_POSTSUPERSCRIPT roman_Ar, finding reasonable agreement to experimental data for the 12C target and an overestimation of a factor 3similar-toabsent3\sim 3∼ 3 for the more neutron-rich and 54Ca, which is similar to the results found in two-proton knockout experiments with 9Be and 12C targets.

I Introduction

Two-proton knockout reactions from neutron-rich nuclei using 9Be and 12C at intermediate energies have been shown to proceed as a direct reaction Bazin et al. (2003) and were shown to be able to populate very exotic nuclei via the removal of two protons from already proton-deficient nuclei Bazin et al. (2003); Fridmann et al. (2005). Equivalently, two-neutron knockout reactions from proton-rich nuclei have been used to study very neutron-rich nuclei Yoneda et al. (2006). The analysis of these reactions using an eikonal sudden description Tostevin et al. (2004a, b) has yielded significant results on the structure of these nuclei Fridmann et al. (2006); Gade et al. (2006, 2007); Bastin et al. (2007); Adrich et al. (2008) and on the effect of nuclear correlations on the observables of the reactions, and therefore their value as a probe of these correlations Tostevin et al. (2004b); Simpson et al. (2009); Simpson and Tostevin (2009); Yoneda et al. (2006); Wimmer et al. (2012). The development of new hydrogen-target detectors, such as active targets Beceiro-Novo et al. (2015), or MINOS Obertelli et al. (2014), where a thick liquid-hydrogen target is coupled to a vertex tracker for the recoil protons, has opened the use of proton-induced reactions as reliable probes to explore exotic nuclei, where the reaction mechanism can be explored and understood thanks to the tracking of the paths of the outgoing particles. Therefore, two-proton removal reactions with proton targets, or (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ), appear as an appealing probe to produce exotic nuclei by removing two-protons from already proton-deficient species, and to be able to study their properties thanks to the simpler reaction mechanism and the possibility for proton-tracking. Unfortunately, as was the case for one-neutron removal, the reaction models used with heavier targets Hansen and Tostevin (2003) are not applicable for the case of the proton target, due to the significant recoil of the target proton. Models considering a “quasi-free” interaction between removed and target protons have been more successful in the description of the experimental data for reactions with proton targets (p,2p)𝑝2𝑝(p,2p)( italic_p , 2 italic_p ) Jacob and Maris (1966); Aumann et al. (2013); Moro (2015); Ogata et al. (2015), so a similar approach for (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) reactions seems promising and is required in order to fully exploit (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) reactions as a spectroscopic tool using the experimental possibility that hydrogen active targets provide. This need has been indicated in previous publications Taniuchi et al. (2019) where the lack of such a theory has hindered the analysis of the experimental data.

This work aims to provide a theoretical formalism for (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) reactions, based on the assumption of “quasi-free” collision between the target and removed protons, and is structured as follows: Section II presents the theoretical formalism and briefly shows its derivation, section III presents calculations of (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) reactions for the stable 12C target as validation of the theory and results for the neutron-rich targets 28Mg and 54Ca. Finally, section IV presents the conclusions and summary of this work as well as future extensions.

II Theoretical framework

Refer to caption
Figure 1: Schematic description of the reaction model considered in this work, two sequential collisions of the incoming proton with the removed protons. The labels used in this work are presented in black, indicating as well the wave-numbers of the protons at each step of the process. The elements appearing in the T𝑇Titalic_T-matrix for the reaction can be associated to the different steps of the trajectory of the protons and are shown close to them in purple and white.

We consider a process p+(A+2)3p+A𝑝𝐴23𝑝𝐴p+(A+2)\rightarrow 3p+Aitalic_p + ( italic_A + 2 ) → 3 italic_p + italic_A, where a projectile proton collides with the target nucleus A+2𝐴2A+2italic_A + 2 removing two protons from it, with the remaining nucleus A𝐴Aitalic_A remaining bound. For the derivation, we will assume an infinite mass for A𝐴Aitalic_A. Following the results from Frotscher et al. (2020), we describe the process for this collision as two sequential and independent collisions between the projectile proton and two protons of the target, with the residual nucleus A𝐴Aitalic_A remaining as an inert spectator. We will also assume that the reaction occurs fast enough for the internal degrees of freedom of A𝐴Aitalic_A to remain frozen during the collision, so that the removal of the two protons does not alter the state of A𝐴Aitalic_A. For ease of description, in this derivation, the protons will be treated as distinguishable, since, following Goldberger and Watson Goldberger and Watson (1967), it is sufficient to consider their antisymmetrization in the proton-proton interaction Vppsubscript𝑉𝑝𝑝V_{pp}italic_V start_POSTSUBSCRIPT italic_p italic_p end_POSTSUBSCRIPT. As such, p0subscript𝑝0p_{0}italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT corresponds to the incoming proton, with momentum 𝐤0Planck-constant-over-2-pisubscript𝐤0\hbar\mathbf{k}_{0}roman_ℏ bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. p0subscript𝑝0p_{0}italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is then assumed to collide with the first proton p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, which is expelled with an asymptotic momentum 𝐤1Planck-constant-over-2-pisubscript𝐤1\hbar\mathbf{k}_{1}roman_ℏ bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT. Then p0subscript𝑝0p_{0}italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT collides with a second proton p2subscript𝑝2p_{2}italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, and both escape the nucleus with momenta 𝐤0fPlanck-constant-over-2-pisubscript𝐤0𝑓\hbar\mathbf{k}_{0f}roman_ℏ bold_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT and 𝐤2Planck-constant-over-2-pisubscript𝐤2\hbar\mathbf{k}_{2}roman_ℏ bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT respectively. Through momentum conservation the residual nucleus A𝐴Aitalic_A is left with momentum 𝐤A=(𝐤0𝐤0f𝐤1𝐤2)Planck-constant-over-2-pisubscript𝐤𝐴Planck-constant-over-2-pisubscript𝐤0subscript𝐤0𝑓subscript𝐤1subscript𝐤2\hbar\mathbf{k}_{A}=\hbar(\mathbf{k}_{0}-\mathbf{k}_{0f}-\mathbf{k}_{1}-% \mathbf{k}_{2})roman_ℏ bold_k start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT = roman_ℏ ( bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - bold_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT - bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ). A diagram of the process is shown in Figure 1. Schematically, the transition matrix for this process can be expressed as:

Tdξi(d𝐫i)χp0(𝐫0f,𝐤0f)χp1(𝐫1,𝐤1)χp2(𝐫2,𝐤2)ΦA(ξA)[V02GV01](𝐫0f,𝐫1,𝐫2,ξA;𝐫0,ξ)Φ(A+2)(ξ)χp0(𝐫0,𝐤0),similar-to-or-equals𝑇d𝜉subscriptproduct𝑖dsubscript𝐫𝑖subscriptsuperscript𝜒subscript𝑝0subscript𝐫0𝑓subscript𝐤0𝑓subscriptsuperscript𝜒subscript𝑝1subscript𝐫1subscript𝐤1subscriptsuperscript𝜒subscript𝑝2subscript𝐫2subscript𝐤2subscriptΦ𝐴subscript𝜉𝐴delimited-[]subscript𝑉02𝐺subscript𝑉01subscript𝐫0𝑓subscript𝐫1subscript𝐫2subscript𝜉𝐴subscript𝐫0𝜉subscriptΦ𝐴2𝜉subscript𝜒subscript𝑝0subscript𝐫0subscript𝐤0\begin{split}T\simeq\int&\mathrm{d}\xi\prod_{i}\left(\mathrm{d}\mathbf{r}_{i}% \right)\chi^{*}_{p_{0}}(\mathbf{r}_{0f},\mathbf{k}_{0f})\chi^{*}_{p_{1}}(% \mathbf{r}_{1},\mathbf{k}_{1})\chi^{*}_{p_{2}}(\mathbf{r}_{2},\mathbf{k}_{2})% \Phi_{A}(\xi_{A})\\ &[V_{02}GV_{01}](\mathbf{r}_{0f},\mathbf{r}_{1},\mathbf{r}_{2},\xi_{A};\mathbf% {r}_{0},\xi)\Phi_{(A+2)}(\xi)\chi_{p_{0}}(\mathbf{r}_{0},\mathbf{k}_{0}),\end{split}start_ROW start_CELL italic_T ≃ ∫ end_CELL start_CELL roman_d italic_ξ ∏ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( roman_d bold_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) italic_χ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT ) italic_χ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_χ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) roman_Φ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ( italic_ξ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL [ italic_V start_POSTSUBSCRIPT 02 end_POSTSUBSCRIPT italic_G italic_V start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT ] ( bold_r start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_ξ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_ξ ) roman_Φ start_POSTSUBSCRIPT ( italic_A + 2 ) end_POSTSUBSCRIPT ( italic_ξ ) italic_χ start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , end_CELL end_ROW (1)

where d𝐫idsubscript𝐫𝑖\mathrm{d}\mathbf{r}_{i}roman_d bold_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT denotes all radial variables involved: 𝐫0f,𝐫1,𝐫2,𝐫0subscript𝐫0𝑓subscript𝐫1subscript𝐫2subscript𝐫0\mathbf{r}_{0f},\mathbf{r}_{1},\mathbf{r}_{2},\mathbf{r}_{0}bold_r start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT; χpisubscript𝜒subscript𝑝𝑖\chi_{p_{i}}italic_χ start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT indicate the wavefunctions for the relative motion between the corresponding particles and residual nucleus A𝐴Aitalic_A (with infinite mass) for the corresponding asymptotic momenta, ξ𝜉\xiitalic_ξ and ξAsubscript𝜉𝐴\xi_{A}italic_ξ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT correspond to the internal coordinates of (A+2)𝐴2(A+2)( italic_A + 2 ) and A𝐴Aitalic_A respectively, Vijsubscript𝑉𝑖𝑗V_{ij}italic_V start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT is the interaction potential between protons i𝑖iitalic_i and j𝑗jitalic_j and G𝐺Gitalic_G is the propagator of the system. Due to the spectator approximation for A𝐴Aitalic_A its internal coordinates ξAsubscript𝜉𝐴\xi_{A}italic_ξ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT are not modified during the reaction, so neither V𝑉Vitalic_V nor G𝐺Gitalic_G depend on them. We can expand ξ(A+2)=ξA,𝐫1i,𝐫2isubscript𝜉𝐴2subscript𝜉𝐴subscript𝐫1𝑖subscript𝐫2𝑖\xi_{(A+2)}=\xi_{A},\mathbf{r}_{1i},\mathbf{r}_{2i}italic_ξ start_POSTSUBSCRIPT ( italic_A + 2 ) end_POSTSUBSCRIPT = italic_ξ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 italic_i end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 italic_i end_POSTSUBSCRIPT , where 𝐫jisubscript𝐫𝑗𝑖\mathbf{r}_{ji}bold_r start_POSTSUBSCRIPT italic_j italic_i end_POSTSUBSCRIPT is the position of proton j𝑗jitalic_j “within” (A+2)𝐴2(A+2)( italic_A + 2 ). Thus, we can compute the overlap between (A+2)𝐴2(A+2)( italic_A + 2 ) and A𝐴Aitalic_A:

dξAA(ξA)|(A+2)(ξA,𝐫1i,𝐫2i)=ϕ12(𝐫1i,𝐫2i).differential-dsubscript𝜉𝐴inner-product𝐴subscript𝜉𝐴𝐴2subscript𝜉𝐴subscript𝐫1𝑖subscript𝐫2𝑖subscriptitalic-ϕ12subscript𝐫1𝑖subscript𝐫2𝑖\int\mathrm{d}\xi_{A}\braket{A(\xi_{A})}{(A+2)(\xi_{A},\mathbf{r}_{1i},\mathbf% {r}_{2i})}=\phi_{12}(\mathbf{r}_{1i},\mathbf{r}_{2i}).∫ roman_d italic_ξ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ⟨ start_ARG italic_A ( italic_ξ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ) end_ARG | start_ARG ( italic_A + 2 ) ( italic_ξ start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 italic_i end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 italic_i end_POSTSUBSCRIPT ) end_ARG ⟩ = italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 italic_i end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 italic_i end_POSTSUBSCRIPT ) . (2)

We note that ϕ12(𝐫1i,𝐫2i)subscriptitalic-ϕ12subscript𝐫1𝑖subscript𝐫2𝑖\phi_{12}(\mathbf{r}_{1i},\mathbf{r}_{2i})italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 italic_i end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 italic_i end_POSTSUBSCRIPT ) is independent of the reaction and can be obtained via structure calculations such as nuclear shell model Tostevin et al. (2004a). This results in the T𝑇Titalic_T-matrix:

Td𝐫0fd𝐫1d𝐫2d𝐫0d𝐫1id𝐫2iχp0(𝐫0f,𝐤0f)χp1(𝐫1,𝐤1)×χp2(𝐫2,𝐤2)[V02GV01](𝐫0f,𝐫1,𝐫2;𝐫0,𝐫1i,𝐫2)×χp0(𝐫0,𝐤0)ϕ12(𝐫1i,𝐫2i).similar-to-or-equals𝑇dsubscript𝐫0𝑓dsubscript𝐫1dsubscript𝐫2dsubscript𝐫0dsubscript𝐫1𝑖dsubscript𝐫2𝑖subscriptsuperscript𝜒subscript𝑝0subscript𝐫0𝑓subscript𝐤0𝑓subscriptsuperscript𝜒subscript𝑝1subscript𝐫1subscript𝐤1subscriptsuperscript𝜒subscript𝑝2subscript𝐫2subscript𝐤2delimited-[]subscript𝑉02𝐺subscript𝑉01subscript𝐫0𝑓subscript𝐫1subscript𝐫2subscript𝐫0subscript𝐫1𝑖subscript𝐫2subscript𝜒subscript𝑝0subscript𝐫0subscript𝐤0subscriptitalic-ϕ12subscript𝐫1𝑖subscript𝐫2𝑖\begin{split}T\simeq\int&\mathrm{d}\mathbf{r}_{0f}\mathrm{d}\mathbf{r}_{1}% \mathrm{d}\mathbf{r}_{2}\mathrm{d}\mathbf{r}_{0}\mathrm{d}\mathbf{r}_{1i}% \mathrm{d}\mathbf{r}_{2i}\chi^{*}_{p_{0}}(\mathbf{r}_{0f},\mathbf{k}_{0f})\chi% ^{*}_{p_{1}}(\mathbf{r}_{1},\mathbf{k}_{1})\\ \times&\chi^{*}_{p_{2}}(\mathbf{r}_{2},\mathbf{k}_{2})[V_{02}GV_{01}](\mathbf{% r}_{0f},\mathbf{r}_{1},\mathbf{r}_{2};\mathbf{r}_{0},\mathbf{r}_{1i},\mathbf{r% }_{2})\\ \times&\chi_{p_{0}}(\mathbf{r}_{0},\mathbf{k}_{0})\phi_{12}(\mathbf{r}_{1i},% \mathbf{r}_{2i}).\end{split}start_ROW start_CELL italic_T ≃ ∫ end_CELL start_CELL roman_d bold_r start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 1 italic_i end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 2 italic_i end_POSTSUBSCRIPT italic_χ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT ) italic_χ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL italic_χ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) [ italic_V start_POSTSUBSCRIPT 02 end_POSTSUBSCRIPT italic_G italic_V start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT ] ( bold_r start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 italic_i end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL italic_χ start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 italic_i end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 italic_i end_POSTSUBSCRIPT ) . end_CELL end_ROW (3)

We now approximate propagator G𝐺Gitalic_G as Joachaim (1975):

G=G01+GV01G01G01,𝐺subscript𝐺01𝐺subscript𝑉01subscript𝐺01similar-to-or-equalssubscript𝐺01G=G_{01}+GV_{01}G_{01}\simeq G_{01},italic_G = italic_G start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT + italic_G italic_V start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT italic_G start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT ≃ italic_G start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT , (4)

which corresponds to:

G01=1EHV01+iϵ==1ET1T2T0V02U0AU1AV2AV12+iϵ,subscript𝐺011𝐸𝐻subscript𝑉01𝑖italic-ϵ1𝐸subscript𝑇1subscript𝑇2subscript𝑇0subscript𝑉02subscript𝑈0𝐴subscript𝑈1𝐴subscript𝑉2𝐴subscript𝑉12𝑖italic-ϵ\begin{split}&G_{01}=\frac{1}{E-H-V_{01}+i\epsilon}=\\ &=\frac{1}{E-T_{1}-T_{2}-T_{0}-V_{02}-U_{0A}-U_{1A}-V_{2A}-V_{12}+i\epsilon},% \end{split}start_ROW start_CELL end_CELL start_CELL italic_G start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG italic_E - italic_H - italic_V start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT + italic_i italic_ϵ end_ARG = end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL = divide start_ARG 1 end_ARG start_ARG italic_E - italic_T start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_T start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_V start_POSTSUBSCRIPT 02 end_POSTSUBSCRIPT - italic_U start_POSTSUBSCRIPT 0 italic_A end_POSTSUBSCRIPT - italic_U start_POSTSUBSCRIPT 1 italic_A end_POSTSUBSCRIPT - italic_V start_POSTSUBSCRIPT 2 italic_A end_POSTSUBSCRIPT - italic_V start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT + italic_i italic_ϵ end_ARG , end_CELL end_ROW (5)

where the infinite mass approximation for A allows us to remove its kinetic energy. It is now our goal to reduce this four-body propagator to a one-body one. For this, we note that (since A𝐴A\rightarrow\inftyitalic_A → ∞) we can group T2+V2AH2Asimilar-to-or-equalssubscript𝑇2subscript𝑉2𝐴subscript𝐻2𝐴T_{2}+V_{2A}\simeq H_{2A}italic_T start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + italic_V start_POSTSUBSCRIPT 2 italic_A end_POSTSUBSCRIPT ≃ italic_H start_POSTSUBSCRIPT 2 italic_A end_POSTSUBSCRIPT, and approximate for the propagator

H2Aϕ12(𝐫1i,𝐫2i)E2Aϕ12(𝐫1i,𝐫2i)S2p2ϕ12(𝐫1i,𝐫2i).similar-tosubscript𝐻2𝐴subscriptitalic-ϕ12subscript𝐫1𝑖subscript𝐫2𝑖subscript𝐸2𝐴subscriptitalic-ϕ12subscript𝐫1𝑖subscript𝐫2𝑖similar-tosubscript𝑆2𝑝2subscriptitalic-ϕ12subscript𝐫1𝑖subscript𝐫2𝑖H_{2A}\phi_{12}(\mathbf{r}_{1i},\mathbf{r}_{2i})\sim E_{2A}\phi_{12}(\mathbf{r% }_{1i},\mathbf{r}_{2i})\sim\dfrac{S_{2p}}{2}\phi_{12}(\mathbf{r}_{1i},\mathbf{% r}_{2i}).italic_H start_POSTSUBSCRIPT 2 italic_A end_POSTSUBSCRIPT italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 italic_i end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 italic_i end_POSTSUBSCRIPT ) ∼ italic_E start_POSTSUBSCRIPT 2 italic_A end_POSTSUBSCRIPT italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 italic_i end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 italic_i end_POSTSUBSCRIPT ) ∼ divide start_ARG italic_S start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 italic_i end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 italic_i end_POSTSUBSCRIPT ) . (6)

This assumption should be reasonable for an energetic beam, as p2subscript𝑝2p_{2}italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT should remain frozen before its collision and E2A<<Emuch-less-thansubscript𝐸2𝐴𝐸E_{2A}<<Eitalic_E start_POSTSUBSCRIPT 2 italic_A end_POSTSUBSCRIPT < < italic_E. We can also group T1+U1A+V12=H1(A+1)subscript𝑇1subscript𝑈1𝐴subscript𝑉12subscript𝐻1𝐴1T_{1}+U_{1A}+V_{12}=H_{1(A+1)}italic_T start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_U start_POSTSUBSCRIPT 1 italic_A end_POSTSUBSCRIPT + italic_V start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT = italic_H start_POSTSUBSCRIPT 1 ( italic_A + 1 ) end_POSTSUBSCRIPT, which can be approximated as

χp1(𝐫1,𝐤1)H1(A+1)=χp1(𝐫1,𝐤1)E1,subscriptsuperscript𝜒subscript𝑝1subscript𝐫1subscript𝐤1subscript𝐻1𝐴1subscriptsuperscript𝜒subscript𝑝1subscript𝐫1subscript𝐤1subscript𝐸1\chi^{*}_{p_{1}}(\mathbf{r}_{1},\mathbf{k}_{1})H_{1(A+1)}=\chi^{*}_{p_{1}}(% \mathbf{r}_{1},\mathbf{k}_{1})E_{1},italic_χ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_H start_POSTSUBSCRIPT 1 ( italic_A + 1 ) end_POSTSUBSCRIPT = italic_χ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_E start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , (7)

and finally U0A+V02U0(A+1)similar-to-or-equalssubscript𝑈0𝐴subscript𝑉02subscript𝑈0𝐴1U_{0A}+V_{02}\simeq U_{0(A+1)}italic_U start_POSTSUBSCRIPT 0 italic_A end_POSTSUBSCRIPT + italic_V start_POSTSUBSCRIPT 02 end_POSTSUBSCRIPT ≃ italic_U start_POSTSUBSCRIPT 0 ( italic_A + 1 ) end_POSTSUBSCRIPT. With these approximations, the propagator reduces to:

G01=1EE1E2AT0U0(A+1)+iϵ==1E0mT0U0(A+1)+iϵ,subscript𝐺011𝐸subscript𝐸1subscript𝐸2𝐴subscript𝑇0subscript𝑈0𝐴1𝑖italic-ϵ1subscript𝐸0𝑚subscript𝑇0subscript𝑈0𝐴1𝑖italic-ϵ\begin{split}G_{01}&=\frac{1}{E-E_{1}-E_{2A}-T_{0}-U_{0(A+1)}+i\epsilon}=\\ &=\frac{1}{E_{0m}-T_{0}-U_{0(A+1)}+i\epsilon},\end{split}start_ROW start_CELL italic_G start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT end_CELL start_CELL = divide start_ARG 1 end_ARG start_ARG italic_E - italic_E start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_E start_POSTSUBSCRIPT 2 italic_A end_POSTSUBSCRIPT - italic_T start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_U start_POSTSUBSCRIPT 0 ( italic_A + 1 ) end_POSTSUBSCRIPT + italic_i italic_ϵ end_ARG = end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL = divide start_ARG 1 end_ARG start_ARG italic_E start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT - italic_T start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_U start_POSTSUBSCRIPT 0 ( italic_A + 1 ) end_POSTSUBSCRIPT + italic_i italic_ϵ end_ARG , end_CELL end_ROW (8)

where E0m=EE1E2Asubscript𝐸0𝑚𝐸subscript𝐸1subscript𝐸2𝐴E_{0m}=E-E_{1}-E_{2A}italic_E start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT = italic_E - italic_E start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_E start_POSTSUBSCRIPT 2 italic_A end_POSTSUBSCRIPT and which can be interpreted as a one-body propagator for p0subscript𝑝0p_{0}italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT between collisions. We will assume the beam energy is high enough that E1subscript𝐸1E_{1}italic_E start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is high enough that the effect of the potential can be ignored, so that E1=mp2+2k12subscript𝐸1superscriptsubscript𝑚𝑝2superscriptPlanck-constant-over-2-pi2superscriptsubscript𝑘12E_{1}=\sqrt{m_{p}^{2}+\hbar^{2}k_{1}^{2}}italic_E start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = square-root start_ARG italic_m start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_ℏ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG (c=1)𝑐1(c=1)( italic_c = 1 ).

We next consider the proton-proton interaction to be of zero range so that

V01(𝐫0,𝐫1;𝐫0,𝐫1)=V~01δ(𝐫0,𝐫1,𝐫0,𝐫1).subscript𝑉01subscriptsuperscript𝐫0subscriptsuperscript𝐫1subscript𝐫0subscript𝐫1subscript~𝑉01𝛿subscriptsuperscript𝐫0subscriptsuperscript𝐫1subscript𝐫0subscript𝐫1V_{01}(\mathbf{r}^{\prime}_{0},\mathbf{r}^{\prime}_{1};\mathbf{r}_{0},\mathbf{% r}_{1})=\tilde{V}_{01}\delta(\mathbf{r}^{\prime}_{0},\mathbf{r}^{\prime}_{1},% \mathbf{r}_{0},\mathbf{r}_{1}).italic_V start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT ( bold_r start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , bold_r start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = over~ start_ARG italic_V end_ARG start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT italic_δ ( bold_r start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , bold_r start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) . (9)

With these approximations, the integral in Eq. 3 is reduced to only two radial coordinates 𝐫1,𝐫2subscript𝐫1subscript𝐫2\mathbf{r}_{1},\mathbf{r}_{2}bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, which can be interpreted as the location of the collisions between p0subscript𝑝0p_{0}italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and between p0subscript𝑝0p_{0}italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and p2subscript𝑝2p_{2}italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT respectively:

Td𝐫1d𝐫2χp0(𝐫2,𝐤0f)χp1(𝐫1,𝐤1)χp2(𝐫2,𝐤2)V~02G0m(E0m,𝐫2;𝐫1)V~01χp0(𝐫1,𝐤0)ϕ12(𝐫1,𝐫2).similar-to-or-equals𝑇differential-dsubscript𝐫1differential-dsubscript𝐫2subscriptsuperscript𝜒subscript𝑝0subscript𝐫2subscript𝐤0𝑓subscriptsuperscript𝜒subscript𝑝1subscript𝐫1subscript𝐤1subscriptsuperscript𝜒subscript𝑝2subscript𝐫2subscript𝐤2subscript~𝑉02subscript𝐺0𝑚subscript𝐸0𝑚subscript𝐫2subscript𝐫1subscript~𝑉01subscript𝜒subscript𝑝0subscript𝐫1subscript𝐤0subscriptitalic-ϕ12subscript𝐫1subscript𝐫2\begin{split}T\simeq&\int\mathrm{d}\mathbf{r}_{1}\mathrm{d}\mathbf{r}_{2}\chi^% {*}_{p_{0}}(\mathbf{r}_{2},\mathbf{k}_{0f})\chi^{*}_{p_{1}}(\mathbf{r}_{1},% \mathbf{k}_{1})\chi^{*}_{p_{2}}(\mathbf{r}_{2},\mathbf{k}_{2})\\ &\tilde{V}_{02}G_{0m}(E_{0m},\mathbf{r}_{2};\mathbf{r}_{1})\tilde{V}_{01}\chi_% {p_{0}}(\mathbf{r}_{1},\mathbf{k}_{0})\phi_{12}(\mathbf{r}_{1},\mathbf{r}_{2})% .\end{split}start_ROW start_CELL italic_T ≃ end_CELL start_CELL ∫ roman_d bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_χ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT ) italic_χ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_χ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL over~ start_ARG italic_V end_ARG start_POSTSUBSCRIPT 02 end_POSTSUBSCRIPT italic_G start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT ( italic_E start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) over~ start_ARG italic_V end_ARG start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT italic_χ start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) . end_CELL end_ROW (10)

We now introduce eikonal expressions for the wavefunctions

χpn(𝐫j,𝐤i)=subscript𝜒subscript𝑝𝑛subscript𝐫𝑗subscript𝐤𝑖absent\displaystyle\chi_{p_{n}}(\mathbf{r}_{j},\mathbf{k}_{i})=italic_χ start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) = 1(2π)3/2eivizjU(𝐛j,z)dzei𝐤0𝐫11superscript2𝜋32superscript𝑒𝑖Planck-constant-over-2-pisubscript𝑣𝑖subscriptsuperscriptsubscript𝑧𝑗𝑈subscript𝐛𝑗𝑧differential-d𝑧superscript𝑒𝑖subscript𝐤0subscript𝐫1\displaystyle\dfrac{1}{(2\pi)^{3/2}}e^{-\frac{i}{\hbar v_{i}}\int^{z_{j}}_{-% \infty}U(\mathbf{b}_{j},z)\mathrm{d}z}e^{i\mathbf{k}_{0}\mathbf{r}_{1}}divide start_ARG 1 end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT end_ARG italic_e start_POSTSUPERSCRIPT - divide start_ARG italic_i end_ARG start_ARG roman_ℏ italic_v start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG ∫ start_POSTSUPERSCRIPT italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - ∞ end_POSTSUBSCRIPT italic_U ( bold_b start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , italic_z ) roman_d italic_z end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_i bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT
=\displaystyle== 1(2π)3/2Sp0(𝐫1,𝐤0)ei𝐤0𝐫11superscript2𝜋32subscript𝑆𝑝0subscript𝐫1subscript𝐤0superscript𝑒𝑖subscript𝐤0subscript𝐫1\displaystyle\dfrac{1}{(2\pi)^{3/2}}S_{p0}(\mathbf{r}_{1},\mathbf{k}_{0})e^{i% \mathbf{k}_{0}\mathbf{r}_{1}}divide start_ARG 1 end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT end_ARG italic_S start_POSTSUBSCRIPT italic_p 0 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT (11)

where z𝑧zitalic_z follows the direction of momentum 𝐤isubscript𝐤𝑖\mathbf{k}_{i}bold_k start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, 𝐛𝐛\mathbf{b}bold_b is the associated impact parameter and U𝑈Uitalic_U is the optical potential between proton and core A𝐴Aitalic_A (As A𝐴A\rightarrow\inftyitalic_A → ∞ we assume UpAUp(A+1)Up(A+2)similar-tosubscript𝑈𝑝𝐴subscript𝑈𝑝𝐴1similar-tosubscript𝑈𝑝𝐴2U_{pA}\sim U_{p(A+1)}\sim U_{p(A+2)}italic_U start_POSTSUBSCRIPT italic_p italic_A end_POSTSUBSCRIPT ∼ italic_U start_POSTSUBSCRIPT italic_p ( italic_A + 1 ) end_POSTSUBSCRIPT ∼ italic_U start_POSTSUBSCRIPT italic_p ( italic_A + 2 ) end_POSTSUBSCRIPT). With this approximation we obtain for the matrix element:

T=1(2π)6d𝐤0md𝐫1d𝐫2S(𝐫2,𝐤0f)S(𝐫2,𝐤2)V02~G0m(E0m,𝐫2;𝐫1)V01~S(𝐫1,𝐤1)S(𝐫1,𝐤0)ϕ12(𝐫1,𝐫2)ei(𝐤0𝐤1)𝐫1ei(𝐤0f+𝐤2)𝐫2.𝑇1superscript2𝜋6differential-dsubscript𝐤0𝑚differential-dsubscript𝐫1differential-dsubscript𝐫2superscript𝑆subscript𝐫2subscript𝐤0𝑓superscript𝑆subscript𝐫2subscript𝐤2~subscript𝑉02subscript𝐺0𝑚subscript𝐸0𝑚subscript𝐫2subscript𝐫1~subscript𝑉01superscript𝑆subscript𝐫1subscript𝐤1𝑆subscript𝐫1subscript𝐤0subscriptitalic-ϕ12subscript𝐫1subscript𝐫2superscript𝑒𝑖subscript𝐤0subscript𝐤1subscript𝐫1superscript𝑒𝑖subscript𝐤0𝑓subscript𝐤2subscript𝐫2\begin{split}T&=\dfrac{1}{(2\pi)^{6}}\int\mathrm{d}\mathbf{k}_{0m}\mathrm{d}% \mathbf{r}_{1}\mathrm{d}\mathbf{r}_{2}S^{*}(\mathbf{r}_{2},\mathbf{k}_{0f})S^{% *}(\mathbf{r}_{2},\mathbf{k}_{2})\\ &\tilde{V_{02}}G_{0m}(E_{0m},\mathbf{r}_{2};\mathbf{r}_{1})\tilde{V_{01}}S^{*}% (\mathbf{r}_{1},\mathbf{k}_{1})S(\mathbf{r}_{1},\mathbf{k}_{0})\phi_{12}(% \mathbf{r}_{1},\mathbf{r}_{2})\\ &e^{i(\mathbf{k}_{0}-\mathbf{k}_{1})\mathbf{r}_{1}}e^{-i(\mathbf{k}_{0f}+% \mathbf{k}_{2})\mathbf{r}_{2}}.\end{split}start_ROW start_CELL italic_T end_CELL start_CELL = divide start_ARG 1 end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT end_ARG ∫ roman_d bold_k start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_S start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT ) italic_S start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL over~ start_ARG italic_V start_POSTSUBSCRIPT 02 end_POSTSUBSCRIPT end_ARG italic_G start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT ( italic_E start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) over~ start_ARG italic_V start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT end_ARG italic_S start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_S ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL italic_e start_POSTSUPERSCRIPT italic_i ( bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i ( bold_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT + bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT . end_CELL end_ROW (12)

Next let us present the formula for the total cross section, where momentum conservation has already been considered to nullify the integral over 𝐤𝐀subscript𝐤𝐀\mathbf{k_{A}}bold_k start_POSTSUBSCRIPT bold_A end_POSTSUBSCRIPT:

σ=12Ji^2S,s0,s0fd𝐤1d𝐤0fd𝐤2δ(EfEi)(2π)4v0|T|2,𝜎12superscript^subscript𝐽𝑖2subscript𝑆subscript𝑠0subscript𝑠0𝑓differential-dsubscript𝐤1differential-dsubscript𝐤0𝑓differential-dsubscript𝐤2𝛿subscript𝐸𝑓subscript𝐸𝑖superscript2𝜋4Planck-constant-over-2-pisubscript𝑣0superscript𝑇2\sigma=\dfrac{1}{2\hat{J_{i}}^{2}}\sum_{S,s_{0},s_{0f}}\int\mathrm{d}\mathbf{k% }_{1}\mathrm{d}\mathbf{k}_{0f}\mathrm{d}\mathbf{k}_{2}\delta(E_{f}-E_{i})% \dfrac{(2\pi)^{4}}{\hbar v_{0}}\left|T\right|^{2},italic_σ = divide start_ARG 1 end_ARG start_ARG 2 over^ start_ARG italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∑ start_POSTSUBSCRIPT italic_S , italic_s start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_s start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT end_POSTSUBSCRIPT ∫ roman_d bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_d bold_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT roman_d bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_δ ( italic_E start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT - italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) divide start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG roman_ℏ italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG | italic_T | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (13)

where Jisubscript𝐽𝑖J_{i}italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is the angular momentum of (A+2)𝐴2(A+2)( italic_A + 2 ) Ji^=2Ji+1^subscript𝐽𝑖2subscript𝐽𝑖1\hat{J_{i}}=\sqrt{2J_{i}+1}over^ start_ARG italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG = square-root start_ARG 2 italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT + 1 end_ARG, the sum over S𝑆Sitalic_S corresponds to the sum over the spin projections of (A+2),A,p1,p2𝐴2𝐴subscript𝑝1subscript𝑝2(A+2),A,p_{1},p_{2}( italic_A + 2 ) , italic_A , italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and the sums for the final and initial spin of proton p0subscript𝑝0p_{0}italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT are left explicit. We will now express d𝐤𝟎𝐟dsubscript𝐤0𝐟\mathrm{d}\mathbf{k_{0f}}roman_d bold_k start_POSTSUBSCRIPT bold_0 bold_f end_POSTSUBSCRIPT in spherical coordinates and will integrate over k0fsubscript𝑘0𝑓k_{0f}italic_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT using δ(EfEi)𝛿subscript𝐸𝑓subscript𝐸𝑖\delta(E_{f}-E_{i})italic_δ ( italic_E start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT - italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ). Thus we obtain

dk0fδ(EfEi)(k0f)=(k¯0f)1Efk0f|k¯0f==12c2k¯0fϵ¯0f+2c2k¯0fϵ¯22c2𝐤¯0f(𝐤0𝐤1𝐤A)k0fϵ¯2(k¯0f)j(𝐤^0f,𝐤1,𝐤A)(k¯0f)differential-dsubscript𝑘0𝑓𝛿subscript𝐸𝑓subscript𝐸𝑖subscript𝑘0𝑓evaluated-atsubscript¯𝑘0𝑓1subscript𝐸𝑓subscript𝑘0𝑓subscript¯𝑘0𝑓1superscriptPlanck-constant-over-2-pi2superscript𝑐2subscript¯𝑘0𝑓subscript¯italic-ϵ0𝑓superscriptPlanck-constant-over-2-pi2superscript𝑐2subscript¯𝑘0𝑓subscript¯italic-ϵ2superscriptPlanck-constant-over-2-pi2superscript𝑐2subscript¯𝐤0𝑓subscript𝐤0subscript𝐤1subscript𝐤𝐴subscript𝑘0𝑓subscript¯italic-ϵ2subscript¯𝑘0𝑓𝑗subscript^𝐤0𝑓subscript𝐤1subscript𝐤𝐴subscript¯𝑘0𝑓\begin{split}&\int\mathrm{d}k_{0f}\delta(E_{f}-E_{i})\mathcal{F}(k_{0f})=% \mathcal{F}(\bar{k}_{0f})\left.\dfrac{1}{\dfrac{\partial E_{f}}{\partial k_{0f% }}}\right|_{\bar{k}_{0f}}=\\ &=\dfrac{1}{\dfrac{\hbar^{2}c^{2}\bar{k}_{0f}}{\bar{\epsilon}_{0f}}+\dfrac{% \hbar^{2}c^{2}\bar{k}_{0f}}{\bar{\epsilon}_{2}}-\dfrac{\hbar^{2}c^{2}\bar{% \mathbf{k}}_{0f}(\mathbf{k}_{0}-\mathbf{k}_{1}-\mathbf{k}_{A})}{{k_{0f}}\bar{% \epsilon}_{2}}}\mathcal{F}(\bar{k}_{0f})\equiv\\ &\equiv j(\hat{\mathbf{k}}_{0f},\mathbf{k}_{1},\mathbf{k}_{A})\mathcal{F}(\bar% {k}_{0f})\end{split}start_ROW start_CELL end_CELL start_CELL ∫ roman_d italic_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT italic_δ ( italic_E start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT - italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) caligraphic_F ( italic_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT ) = caligraphic_F ( over¯ start_ARG italic_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT ) divide start_ARG 1 end_ARG start_ARG divide start_ARG ∂ italic_E start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT end_ARG end_ARG | start_POSTSUBSCRIPT over¯ start_ARG italic_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT end_POSTSUBSCRIPT = end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL = divide start_ARG 1 end_ARG start_ARG divide start_ARG roman_ℏ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over¯ start_ARG italic_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT end_ARG start_ARG over¯ start_ARG italic_ϵ end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT end_ARG + divide start_ARG roman_ℏ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over¯ start_ARG italic_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT end_ARG start_ARG over¯ start_ARG italic_ϵ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG - divide start_ARG roman_ℏ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over¯ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - bold_k start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ) end_ARG start_ARG italic_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT over¯ start_ARG italic_ϵ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG end_ARG caligraphic_F ( over¯ start_ARG italic_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT ) ≡ end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL ≡ italic_j ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ) caligraphic_F ( over¯ start_ARG italic_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT ) end_CELL end_ROW (14)

where j(𝐤^0f,𝐤1,𝐤A)𝑗subscript^𝐤0𝑓subscript𝐤1subscript𝐤𝐴j(\hat{\mathbf{k}}_{0f},\mathbf{k}_{1},\mathbf{k}_{A})italic_j ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ) corresponds to these kinematic factors and (k0f)subscript𝑘0𝑓\mathcal{F}(k_{0f})caligraphic_F ( italic_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT ) is meant to denote the rest of integral in Eq. 13 and where we have applied that k22=k0f22𝐤0f(𝐤0𝐤1𝐤A)+|𝐤0𝐤1𝐤A|2superscriptsubscript𝑘22superscriptsubscript𝑘0𝑓22subscript𝐤0𝑓subscript𝐤0subscript𝐤1subscript𝐤𝐴superscriptsubscript𝐤0subscript𝐤1subscript𝐤𝐴2k_{2}^{2}=k_{0f}^{2}-2\mathbf{k}_{0f}(\mathbf{k}_{0}-\mathbf{k}_{1}-\mathbf{k}% _{A})+|\mathbf{k}_{0}-\mathbf{k}_{1}-\mathbf{k}_{A}|^{2}italic_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 2 bold_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT ( bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - bold_k start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ) + | bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - bold_k start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and 𝐤¯0f,ϵ¯0f,ϵ¯2subscript¯𝐤0𝑓subscript¯italic-ϵ0𝑓subscript¯italic-ϵ2\bar{\mathbf{k}}_{0f},\bar{\epsilon}_{0f},\bar{\epsilon}_{2}over¯ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , over¯ start_ARG italic_ϵ end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , over¯ start_ARG italic_ϵ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT correspond to the momentum for p0fsubscript𝑝0𝑓p_{0f}italic_p start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT and energies of p0fsubscript𝑝0𝑓p_{0f}italic_p start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT and p2subscript𝑝2p_{2}italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT given by energy conservation.

To proceed with the derivation, we will now introduce the strong approximation that all factors in Eqs (12) and (13) are slowly varying functions of the momentum compared to the exponentials in Eq. 12 and can be replaced by a suitable average value which may depend on 𝐫𝟏,𝐫𝟐subscript𝐫1subscript𝐫2\mathbf{r_{1}},\mathbf{r_{2}}bold_r start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT:

S(𝐫i,𝐤j)S(𝐤^0f,𝐫1,𝐫2)j(𝐤^0f,𝐤1,𝐤A)j(𝐤^0f,𝐫1,𝐫2),similar-to𝑆subscript𝐫𝑖subscript𝐤𝑗𝑆subscript^𝐤0𝑓subscript𝐫1subscript𝐫2𝑗subscript^𝐤0𝑓subscript𝐤1subscript𝐤𝐴similar-to𝑗subscript^𝐤0𝑓subscript𝐫1subscript𝐫2\begin{split}S(\mathbf{r}_{i},\mathbf{k}_{j})&\sim S(\hat{\mathbf{k}}_{0f},% \mathbf{r}_{1},\mathbf{r}_{2})\\ j(\hat{\mathbf{k}}_{0f},\mathbf{k}_{1},\mathbf{k}_{A})&\sim j(\hat{\mathbf{k}}% _{0f},\mathbf{r}_{1},\mathbf{r}_{2}),\end{split}start_ROW start_CELL italic_S ( bold_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) end_CELL start_CELL ∼ italic_S ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL italic_j ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_k start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ) end_CELL start_CELL ∼ italic_j ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) , end_CELL end_ROW (15)

which is a general expression, although some of the S𝑆Sitalic_S-matrices may not depend on all variables. With this approximation, Eq. 12 can be expressed as:

T1(2π)6d𝐫1d𝐫2S0f~(𝐤^0f,𝐫2,𝐫1)S2~(𝐤^0f,𝐫2,𝐫1)×G~0m(E0m,𝐫2;𝐫1)S1~(𝐫1,𝐫2)S0~(𝐫1)|ϕ12|2(𝐫1,𝐫2)V01~(𝐫1;𝐫2)V02~(𝐫1;𝐫2)ei(𝐤0𝐤1)𝐫1ei(𝐤0f+𝐤2)𝐫2,similar-to-or-equals𝑇1superscript2𝜋6differential-dsubscript𝐫1differential-dsubscript𝐫2superscript~subscript𝑆0𝑓subscript^𝐤0𝑓subscript𝐫2subscript𝐫1superscript~subscript𝑆2subscript^𝐤0𝑓subscript𝐫2subscript𝐫1subscript~𝐺0𝑚subscript𝐸0𝑚subscript𝐫2subscript𝐫1superscript~subscript𝑆1subscript𝐫1subscript𝐫2~subscript𝑆0subscript𝐫1superscriptsubscriptitalic-ϕ122subscript𝐫1subscript𝐫2~subscript𝑉01subscript𝐫1subscript𝐫2~subscript𝑉02subscript𝐫1subscript𝐫2superscript𝑒𝑖subscript𝐤0subscript𝐤1subscript𝐫1superscript𝑒𝑖subscript𝐤0𝑓subscript𝐤2subscript𝐫2\begin{split}T&\simeq\dfrac{1}{(2\pi)^{6}}\int\mathrm{d}\mathbf{r}_{1}\mathrm{% d}\mathbf{r}_{2}\tilde{S_{0f}}^{*}(\hat{\mathbf{k}}_{0f},\mathbf{r}_{2},% \mathbf{r}_{1})\tilde{S_{2}}^{*}(\hat{\mathbf{k}}_{0f},\mathbf{r}_{2},\mathbf{% r}_{1})\\ \times&\tilde{G}_{0m}(E_{0m},\mathbf{r}_{2};\mathbf{r}_{1})\tilde{S_{1}}^{*}(% \mathbf{r}_{1},\mathbf{r}_{2})\tilde{S_{0}}(\mathbf{r}_{1})|\phi_{12}|^{2}(% \mathbf{r}_{1},\mathbf{r}_{2})\\ &\tilde{V_{01}}(\mathbf{r}_{1};\mathbf{r}_{2})\tilde{V_{02}}(\mathbf{r}_{1};% \mathbf{r}_{2})e^{i(\mathbf{k}_{0}-\mathbf{k}_{1})\mathbf{r}_{1}}e^{-i({% \mathbf{k}_{0f}+\mathbf{k}_{2})\mathbf{r}_{2}}},\end{split}start_ROW start_CELL italic_T end_CELL start_CELL ≃ divide start_ARG 1 end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT end_ARG ∫ roman_d bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT over~ start_ARG italic_S start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT end_ARG start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) over~ start_ARG italic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL over~ start_ARG italic_G end_ARG start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT ( italic_E start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) over~ start_ARG italic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) over~ start_ARG italic_S start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) | italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL over~ start_ARG italic_V start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) over~ start_ARG italic_V start_POSTSUBSCRIPT 02 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i ( bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i ( bold_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT + bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT , end_CELL end_ROW (16)

and now we can perform the integrals over 𝐤1subscript𝐤1\mathbf{k}_{1}bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and 𝐤Asubscript𝐤𝐴\mathbf{k}_{A}bold_k start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT, noting that T𝑇Titalic_T and Tsuperscript𝑇T^{*}italic_T start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT would produce equivalent and complex conjugate exponentials with spatial variables (𝐫1,𝐫2)subscript𝐫1subscript𝐫2(\mathbf{r}_{1},\mathbf{r}_{2})( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) and (𝐫1,𝐫2)subscriptsuperscript𝐫1subscriptsuperscript𝐫2(\mathbf{r}^{\prime}_{1},\mathbf{r}^{\prime}_{2})( bold_r start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ):

d𝐤1d𝐤2[ei𝐤1(𝐫1𝐫1)ei𝐤2(𝐫2𝐫2)]==(2π)6δ(𝐫2𝐫2)δ(𝐫1𝐫1),dsubscript𝐤1dsubscript𝐤2delimited-[]superscript𝑒𝑖subscript𝐤1subscript𝐫1subscriptsuperscript𝐫1superscript𝑒𝑖subscript𝐤2subscript𝐫2subscriptsuperscript𝐫2superscript2𝜋6𝛿subscript𝐫2superscriptsubscript𝐫2𝛿subscript𝐫1subscriptsuperscript𝐫1\begin{split}\int&\mathrm{d}\mathbf{k}_{1}\mathrm{d}\mathbf{k}_{2}\left[e^{-i% \mathbf{k}_{1}(\mathbf{r}_{1}-\mathbf{r}^{\prime}_{1})}e^{i\mathbf{k}_{2}(% \mathbf{r}_{2}-\mathbf{r}^{\prime}_{2})}\right]=\\ &=(2\pi)^{6}\delta(\mathbf{r}_{2}-\mathbf{r}_{2}^{\prime})\delta(\mathbf{r}_{1% }-\mathbf{r}^{\prime}_{1}),\end{split}start_ROW start_CELL ∫ end_CELL start_CELL roman_d bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_d bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT [ italic_e start_POSTSUPERSCRIPT - italic_i bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - bold_r start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_i bold_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - bold_r start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_POSTSUPERSCRIPT ] = end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL = ( 2 italic_π ) start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT italic_δ ( bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_δ ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - bold_r start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) , end_CELL end_ROW (17)

which produces the nice result that the integral reduces to an incoherent sum over the two collision points 𝐫1subscript𝐫1\mathbf{r}_{1}bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and 𝐫2subscript𝐫2\mathbf{r}_{2}bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, which is consistent with the semiclassical assumption.

σ=12J^i2S,s0,s0fd𝐤^0fd𝐫1d𝐫2×1v0(2π)2j(𝐤^0f,𝐫1,𝐫2)|S0f~|2(𝐤^0f,𝐫2,𝐫1)|S2~|2(𝐤^0f,𝐫2,𝐫1)×V02~(𝐫1;𝐫2)|2|G~0m|2(E0m,𝐫2;𝐫1)×|S1~|2(𝐫1,𝐫2)|V01~(𝐫1;𝐫2)|2|S0~(𝐫1)|2|ϕ12|2(𝐫1,𝐫2).\begin{split}\sigma&=\dfrac{1}{2\hat{J}_{i}^{2}}\sum_{S,s_{0},s_{0f}}\int% \mathrm{d}\hat{\mathbf{k}}_{0f}\mathrm{d}\mathbf{r}_{1}\mathrm{d}\mathbf{r}_{2% }\\ \times&\dfrac{1}{\hbar v_{0}(2\pi)^{2}}j(\hat{\mathbf{k}}_{0f},\mathbf{r}_{1},% \mathbf{r}_{2})|\tilde{S_{0f}}|^{2}(\hat{\mathbf{k}}_{0f},\mathbf{r}_{2},% \mathbf{r}_{1})|\tilde{S_{2}}|^{2}(\hat{\mathbf{k}}_{0f},\mathbf{r}_{2},% \mathbf{r}_{1})\\ \times&\tilde{V_{02}}(\mathbf{r}_{1};\mathbf{r}_{2})|^{2}|\tilde{G}_{0m}|^{2}(% E_{0m},\mathbf{r}_{2};\mathbf{r}_{1})\\ \times&|\tilde{S_{1}}|^{2}(\mathbf{r}_{1},\mathbf{r}_{2})|\tilde{V_{01}}(% \mathbf{r}_{1};\mathbf{r}_{2})|^{2}|\tilde{S_{0}}(\mathbf{r}_{1})|^{2}|\phi_{1% 2}|^{2}(\mathbf{r}_{1},\mathbf{r}_{2}).\end{split}start_ROW start_CELL italic_σ end_CELL start_CELL = divide start_ARG 1 end_ARG start_ARG 2 over^ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∑ start_POSTSUBSCRIPT italic_S , italic_s start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_s start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT end_POSTSUBSCRIPT ∫ roman_d over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL divide start_ARG 1 end_ARG start_ARG roman_ℏ italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( 2 italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_j ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | over~ start_ARG italic_S start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT end_ARG | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) | over~ start_ARG italic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL over~ start_ARG italic_V start_POSTSUBSCRIPT 02 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | over~ start_ARG italic_G end_ARG start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_E start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL | over~ start_ARG italic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | over~ start_ARG italic_V start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | over~ start_ARG italic_S start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) . end_CELL end_ROW (18)

We will now use the eikonal expression for the propagator Kawai and Weidenmüller (1992); Goldberger and Watson (1967):

G0m(E0m,𝐫2;𝐫1)=ϵ0m2π2c2eiv0mz1z2U(𝐛,z)dz|𝐫𝟐𝐫𝟏|×eik0m|𝐫𝟐𝐫𝟏|,subscript𝐺0𝑚subscript𝐸0𝑚subscript𝐫2subscript𝐫1subscriptitalic-ϵ0𝑚2𝜋superscriptPlanck-constant-over-2-pi2superscript𝑐2superscript𝑒𝑖Planck-constant-over-2-pisubscript𝑣0𝑚superscriptsubscriptsubscript𝑧1subscript𝑧2𝑈𝐛𝑧differential-d𝑧subscript𝐫2subscript𝐫1superscript𝑒𝑖subscript𝑘0𝑚subscript𝐫2subscript𝐫1\displaystyle\begin{split}G_{0m}(E_{0m},\mathbf{r}_{2};\mathbf{r}_{1})=&-% \dfrac{\epsilon_{0m}}{2\pi\hbar^{2}c^{2}}\dfrac{e^{-\frac{i}{\hbar v_{0m}}\int% _{z_{1}}^{z_{2}}U(\mathbf{b},z)\mathrm{d}z}}{|\mathbf{r_{2}-r_{1}}|}\\ \times&e^{-ik_{0m}|\mathbf{r_{2}-\mathbf{r_{1}}}|},\end{split}start_ROW start_CELL italic_G start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT ( italic_E start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = end_CELL start_CELL - divide start_ARG italic_ϵ start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_π roman_ℏ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG italic_e start_POSTSUPERSCRIPT - divide start_ARG italic_i end_ARG start_ARG roman_ℏ italic_v start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT end_ARG ∫ start_POSTSUBSCRIPT italic_z start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_U ( bold_b , italic_z ) roman_d italic_z end_POSTSUPERSCRIPT end_ARG start_ARG | bold_r start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT - bold_r start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT | end_ARG end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL italic_e start_POSTSUPERSCRIPT - italic_i italic_k start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT | bold_r start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT - bold_r start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT | end_POSTSUPERSCRIPT , end_CELL end_ROW (19)

where we assume the energy E0msubscript𝐸0𝑚E_{0m}italic_E start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT is large enough that the integral can be performed along the straight line between r2subscriptr2\mathrm{r_{2}}roman_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and r1subscriptr1\mathrm{r_{1}}roman_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and k0mPlanck-constant-over-2-pisubscript𝑘0𝑚\hbar k_{0m}roman_ℏ italic_k start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT is the momentum associated to E0msubscript𝐸0𝑚E_{0m}italic_E start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT

We note that a similar method for calculating the (p,2p)𝑝2𝑝(p,2p)( italic_p , 2 italic_p ) total cross section was presented in Aumann et al. (2013). We must now choose the average values to be used for each quantity. For the S-matrices and propagator G0msubscript𝐺0𝑚G_{0m}italic_G start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT Eqs. (11) and (19) provide eikonal expressions. For |V01~(𝐫1;𝐫2)|2superscript~subscript𝑉01subscript𝐫1subscript𝐫22|\tilde{V_{01}}(\mathbf{r}_{1};\mathbf{r}_{2})|^{2}| over~ start_ARG italic_V start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT we will assume on-shell quasifree collisions between the protons and also assume that the removed protons are initially at rest, although their binding energy, taken as S2p/2subscript𝑆2𝑝2S_{2p}/2italic_S start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT / 2, will be taken into account. Therefore in the first collision, proton p0subscript𝑝0p_{0}italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT collides with momentum 𝐤0Planck-constant-over-2-pisubscript𝐤0\hbar\mathbf{k}_{0}roman_ℏ bold_k start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT with proton p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and both are emitted with momenta 𝐤0mPlanck-constant-over-2-pisubscript𝐤0𝑚\hbar\mathbf{k}_{0m}roman_ℏ bold_k start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT and 𝐤1Planck-constant-over-2-pisubscript𝐤1\hbar\mathbf{k}_{1}roman_ℏ bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT respectively. We note that the direction of 𝐤0msubscript𝐤0𝑚\mathbf{k}_{0m}bold_k start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT is parallel to (𝐫2𝐫1)subscript𝐫2subscript𝐫1(\mathbf{r}_{2}-\mathbf{r}_{1})( bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ), which is consistent with eikonal propagation. As such, positions 𝐫1subscript𝐫1\mathbf{r}_{1}bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and 𝐫2subscript𝐫2\mathbf{r}_{2}bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT along with energy and momentum conservation in the quasi-free collision define k0msubscript𝑘0𝑚k_{0m}italic_k start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT and 𝐤1subscript𝐤1\mathbf{k}_{1}bold_k start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT. For the second collision, assuming the proton impacts with momentum 𝐤0mPlanck-constant-over-2-pisubscript𝐤0𝑚\hbar\mathbf{k}_{0m}roman_ℏ bold_k start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT, the direction of 𝐤0fsubscript𝐤0𝑓\mathbf{k}_{0f}bold_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT restricts the value of 𝐤𝟐subscript𝐤2\mathbf{k_{2}}bold_k start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT and the modulus of k0fsubscript𝑘0𝑓k_{0f}italic_k start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT, to be used in S2subscript𝑆2S_{2}italic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and S0fsubscript𝑆0𝑓S_{0f}italic_S start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT, taking 𝐤Asubscript𝐤𝐴\mathbf{k}_{A}bold_k start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT as 0. Finally, |V01~(𝐫1;𝐫2)|2superscript~subscript𝑉01subscript𝐫1subscript𝐫22|\tilde{V_{01}}(\mathbf{r}_{1};\mathbf{r}_{2})|^{2}| over~ start_ARG italic_V start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT can be obtained from the free pp𝑝𝑝p-pitalic_p - italic_p cross section (omitting Coulomb interaction) assuming the Born approximation:

|V01~(𝐫𝟏;𝐫𝟐)|2|T01~(𝐫𝟏;𝐫𝟐)|2.similar-to-or-equalssuperscript~subscript𝑉01subscript𝐫1subscript𝐫22superscript~subscript𝑇01subscript𝐫1subscript𝐫22|\tilde{V_{01}}(\mathbf{r_{1}};\mathbf{r_{2}})|^{2}\simeq|\tilde{T_{01}}(% \mathbf{r_{1}};\mathbf{r_{2}})|^{2}.| over~ start_ARG italic_V start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≃ | over~ start_ARG italic_T start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (20)

Now we move to the NN center of mass system, introducing the adequate Möller factor

|T01~(𝐫1;𝐫2)|2=|T01~(𝐫1;𝐫2)|NN2ϵcmi,0ϵcmi,1ϵcmf,12(mNc2S2p/2)ϵ0ϵ1ϵ0m,superscript~subscript𝑇01subscript𝐫1subscript𝐫22superscriptsubscript~subscript𝑇01subscript𝐫1subscript𝐫2𝑁𝑁2subscriptitalic-ϵ𝑐𝑚𝑖0subscriptitalic-ϵ𝑐𝑚𝑖1subscriptsuperscriptitalic-ϵ2𝑐𝑚𝑓1subscript𝑚𝑁superscript𝑐2subscript𝑆2𝑝2subscriptitalic-ϵ0subscriptitalic-ϵ1subscriptitalic-ϵ0𝑚|\tilde{T_{01}}(\mathbf{r}_{1};\mathbf{r}_{2})|^{2}=|\tilde{T_{01}}(\mathbf{r}% _{1};\mathbf{r}_{2})|_{NN}^{2}\dfrac{\epsilon_{cmi,0}\epsilon_{cmi,1}\epsilon^% {2}_{cmf,1}}{(m_{N}c^{2}-S_{2p}/2)\epsilon_{0}\epsilon_{1}\epsilon_{0m}},| over~ start_ARG italic_T start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = | over~ start_ARG italic_T start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG italic_ϵ start_POSTSUBSCRIPT italic_c italic_m italic_i , 0 end_POSTSUBSCRIPT italic_ϵ start_POSTSUBSCRIPT italic_c italic_m italic_i , 1 end_POSTSUBSCRIPT italic_ϵ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_c italic_m italic_f , 1 end_POSTSUBSCRIPT end_ARG start_ARG ( italic_m start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_S start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT / 2 ) italic_ϵ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ϵ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_ϵ start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT end_ARG , (21)

where, since the collision happens between a proton and a pseudoproton of mass (mNc2S2p/2)subscript𝑚𝑁superscript𝑐2subscript𝑆2𝑝2(m_{N}c^{2}-S_{2p}/2)( italic_m start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_S start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT / 2 ), the center of mass energies of both particles are not the same before the collision (ϵcmi,0subscriptitalic-ϵ𝑐𝑚𝑖0\epsilon_{cmi,0}italic_ϵ start_POSTSUBSCRIPT italic_c italic_m italic_i , 0 end_POSTSUBSCRIPT being the one of the proton and ϵcmi,1subscriptitalic-ϵ𝑐𝑚𝑖1\epsilon_{cmi,1}italic_ϵ start_POSTSUBSCRIPT italic_c italic_m italic_i , 1 end_POSTSUBSCRIPT the one of the pseudoproton) but they are the same after the collision ϵcmfsubscriptitalic-ϵ𝑐𝑚𝑓\epsilon_{cmf}italic_ϵ start_POSTSUBSCRIPT italic_c italic_m italic_f end_POSTSUBSCRIPT, since both outgoing particles are protons. The relation between ϵcmisubscriptitalic-ϵ𝑐𝑚𝑖\epsilon_{cmi}italic_ϵ start_POSTSUBSCRIPT italic_c italic_m italic_i end_POSTSUBSCRIPT and ϵcmfsubscriptitalic-ϵ𝑐𝑚𝑓\epsilon_{cmf}italic_ϵ start_POSTSUBSCRIPT italic_c italic_m italic_f end_POSTSUBSCRIPT can be obtained through the Mandelstan variable s𝑠sitalic_s. We can relate the NN center of mass T matrix to the angular differential free NN cross section Ogata et al. (2015), which we will approximate by the total cross section divided by 4π4𝜋4\pi4 italic_π:

|T01~(𝐫1;𝐫2)|NN2=4(2π)24c4ϵcmf,12dσdΩ=24(2π)24c4ϵcmf,12σNN4πsuperscriptsubscript~subscript𝑇01subscript𝐫1subscript𝐫2𝑁𝑁24superscript2𝜋2superscriptPlanck-constant-over-2-pi4superscript𝑐4superscriptsubscriptitalic-ϵ𝑐𝑚𝑓12d𝜎dΩ24superscript2𝜋2superscriptPlanck-constant-over-2-pi4superscript𝑐4superscriptsubscriptitalic-ϵ𝑐𝑚𝑓12subscript𝜎𝑁𝑁4𝜋|\tilde{T_{01}}(\mathbf{r}_{1};\mathbf{r}_{2})|_{NN}^{2}=\dfrac{4(2\pi)^{2}% \hbar^{4}c^{4}}{\epsilon_{cmf,1}^{2}}\dfrac{\mathrm{d}\sigma}{\mathrm{d}\Omega% }=2\dfrac{4(2\pi)^{2}\hbar^{4}c^{4}}{\epsilon_{cmf,1}^{2}}\dfrac{\sigma_{NN}}{% 4\pi}| over~ start_ARG italic_T start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = divide start_ARG 4 ( 2 italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_ℏ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG italic_ϵ start_POSTSUBSCRIPT italic_c italic_m italic_f , 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG roman_d italic_σ end_ARG start_ARG roman_d roman_Ω end_ARG = 2 divide start_ARG 4 ( 2 italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_ℏ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG italic_ϵ start_POSTSUBSCRIPT italic_c italic_m italic_f , 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG italic_σ start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG start_ARG 4 italic_π end_ARG (22)

(factor 2 due to antisymmetrization) and we get:

|V01~(𝐫1;𝐫2)|24(2π)4c4σNN(𝐫1;𝐫2)×ϵcmi,0ϵcmi,1(mNc2S2p/2)ϵ0ϵ1ϵ0mσNN(𝐫1;𝐫2)f(𝐫1,𝐫2),similar-to-or-equalssuperscript~subscript𝑉01subscript𝐫1subscript𝐫2242𝜋superscriptPlanck-constant-over-2-pi4superscript𝑐4subscript𝜎𝑁𝑁subscript𝐫1subscript𝐫2subscriptitalic-ϵ𝑐𝑚𝑖0subscriptitalic-ϵ𝑐𝑚𝑖1subscript𝑚𝑁superscript𝑐2subscript𝑆2𝑝2subscriptitalic-ϵ0subscriptitalic-ϵ1subscriptitalic-ϵ0𝑚subscript𝜎𝑁𝑁subscript𝐫1subscript𝐫2𝑓subscript𝐫1subscript𝐫2\begin{split}&|\tilde{V_{01}}(\mathbf{r}_{1};\mathbf{r}_{2})|^{2}\simeq 4(2\pi% )\hbar^{4}c^{4}\sigma_{NN}(\mathbf{r}_{1};\mathbf{r}_{2})\\ \times&\dfrac{\epsilon_{cmi,0}\epsilon_{cmi,1}}{(m_{N}c^{2}-S_{2p}/2)\epsilon_% {0}\epsilon_{1}\epsilon_{0m}}\equiv\sigma_{NN}(\mathbf{r}_{1};\mathbf{r}_{2})f% (\mathbf{r}_{1},\mathbf{r}_{2}),\end{split}start_ROW start_CELL end_CELL start_CELL | over~ start_ARG italic_V start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≃ 4 ( 2 italic_π ) roman_ℏ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_σ start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL divide start_ARG italic_ϵ start_POSTSUBSCRIPT italic_c italic_m italic_i , 0 end_POSTSUBSCRIPT italic_ϵ start_POSTSUBSCRIPT italic_c italic_m italic_i , 1 end_POSTSUBSCRIPT end_ARG start_ARG ( italic_m start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_S start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT / 2 ) italic_ϵ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ϵ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_ϵ start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT end_ARG ≡ italic_σ start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_f ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) , end_CELL end_ROW (23)

where the factors multiplying the nucleon-nucleon cross section σNNsubscript𝜎𝑁𝑁\sigma_{NN}italic_σ start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT are included in f(𝐫1,𝐫2)𝑓subscript𝐫1subscript𝐫2f(\mathbf{r}_{1},\mathbf{r}_{2})italic_f ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ). We note that this approximation for the nucleon-nucleon interaction makes it spin-independent. Therefore the projection of the spin of p0subscript𝑝0p_{0}italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT must be the same in the incoming and outgoing channels, so that s0,s0f=2subscriptsubscript𝑠0subscript𝑠0𝑓2\sum_{s_{0},s_{0f}}=2∑ start_POSTSUBSCRIPT italic_s start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_s start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 2, which cancels with the factor 1/2 in Eq. (13). This results in the following expression for the cross section:

σ=1Ji2^Sd𝐤^0fd𝐫1d𝐫2×1v0(2π)2j(𝐤^0f,𝐫1,𝐫2)f02(𝐤^0f,𝐫1,𝐫2)f01(𝐫1,𝐫2)×|S0f~|2(𝐤^0f,𝐫2,𝐫1)|S2~|2(𝐤^0f,𝐫2,𝐫1)σ02(𝐫1;𝐫2)G~0m(𝐫2,𝐫1)|2×|S1~|2(𝐫1,𝐫2)σ01(𝐫1;𝐫2)|S0~(𝐫1)|2|ϕ12|2(𝐫1,𝐫2),\begin{split}\sigma&=\dfrac{1}{\hat{J_{i}^{2}}}\sum_{S}\int\mathrm{d}\hat{% \mathbf{k}}_{0f}\mathrm{d}\mathbf{r}_{1}\mathrm{d}\mathbf{r}_{2}\\ \times&\dfrac{1}{\hbar v_{0}(2\pi)^{2}}j(\hat{\mathbf{k}}_{0f},\mathbf{r}_{1},% \mathbf{r}_{2})f_{02}(\hat{\mathbf{k}}_{0f},\mathbf{r}_{1},\mathbf{r}_{2})f_{0% 1}(\mathbf{r}_{1},\mathbf{r}_{2})\\ \times&|\tilde{S_{0f}}|^{2}(\hat{\mathbf{k}}_{0f},\mathbf{r}_{2},\mathbf{r}_{1% })|\tilde{S_{2}}|^{2}(\hat{\mathbf{k}}_{0f},\mathbf{r}_{2},\mathbf{r}_{1})% \sigma_{02}(\mathbf{r}_{1};\mathbf{r}_{2})\tilde{G}^{m}_{0}(\mathbf{r}_{2},% \mathbf{r}_{1})|^{2}\\ \times&|\tilde{S_{1}}|^{2}(\mathbf{r}_{1},\mathbf{r}_{2})\sigma_{01}(\mathbf{r% }_{1};\mathbf{r}_{2})|\tilde{S_{0}}(\mathbf{r}_{1})|^{2}|\phi_{12}|^{2}(% \mathbf{r}_{1},\mathbf{r}_{2}),\end{split}start_ROW start_CELL italic_σ end_CELL start_CELL = divide start_ARG 1 end_ARG start_ARG over^ start_ARG italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG ∑ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ∫ roman_d over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL divide start_ARG 1 end_ARG start_ARG roman_ℏ italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( 2 italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_j ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_f start_POSTSUBSCRIPT 02 end_POSTSUBSCRIPT ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_f start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL | over~ start_ARG italic_S start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT end_ARG | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) | over~ start_ARG italic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_σ start_POSTSUBSCRIPT 02 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) over~ start_ARG italic_G end_ARG start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL | over~ start_ARG italic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_σ start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | over~ start_ARG italic_S start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) , end_CELL end_ROW (24)

which allows for a very classical interpretation: proton p0subscript𝑝0p_{0}italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT penetrates inside the nucleus with probability |S~0|2superscriptsubscript~𝑆02|\tilde{S}_{0}|^{2}| over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT until it collides with cross section σ01subscript𝜎01\sigma_{01}italic_σ start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT with one of the protons p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT in the nucleus at 𝐫1subscript𝐫1\mathbf{r}_{1}bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and then propagates to 𝐫2subscript𝐫2\mathbf{r}_{2}bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT through |G~0m|2superscriptsubscript~𝐺0𝑚2|\tilde{G}_{0m}|^{2}| over~ start_ARG italic_G end_ARG start_POSTSUBSCRIPT 0 italic_m end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT where it collides with the second proton p2subscript𝑝2p_{2}italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT with cross section σ02subscript𝜎02\sigma_{02}italic_σ start_POSTSUBSCRIPT 02 end_POSTSUBSCRIPT, with the probability of finding both protons at these positions given by |ϕ12|2superscriptsubscriptitalic-ϕ122|\phi_{12}|^{2}| italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT for a certain state of the residual core. Afterwards, all three protons must escape the nucleus with probabilities |S~0f|2,|S~1|2superscriptsubscript~𝑆0𝑓2superscriptsubscript~𝑆12|\tilde{S}_{0f}|^{2},|\tilde{S}_{1}|^{2}| over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , | over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and |S~2|2superscriptsubscript~𝑆22|\tilde{S}_{2}|^{2}| over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, for them to be detected. See Figure 1 for a schematic with the factors associated to the paths of the protons. For a more compact expression, which is also more comparable with the ones in the literature, we will include the spin average and sum in |ϕ12|2superscriptsubscriptitalic-ϕ122|\phi_{12}|^{2}| italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, noting that it is the only term which is sensitive to the spins of the protons:

|ϕ~12|2(𝐫1,𝐫2)=1Ji^2S|ϕ12|2(𝐫1,𝐫2)superscriptsubscript~italic-ϕ122subscript𝐫1subscript𝐫21superscript^subscript𝐽𝑖2subscript𝑆superscriptsubscriptitalic-ϕ122subscript𝐫1subscript𝐫2|\tilde{\phi}_{12}|^{2}(\mathbf{r}_{1},\mathbf{r}_{2})=\dfrac{1}{\hat{J_{i}}^{% 2}}\sum_{S}|\phi_{12}|^{2}(\mathbf{r}_{1},\mathbf{r}_{2})| over~ start_ARG italic_ϕ end_ARG start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) = divide start_ARG 1 end_ARG start_ARG over^ start_ARG italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∑ start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT | italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) (25)

resulting in:

σ=d𝐤^0fd𝐫1d𝐫21v0(2π)2j(𝐤^0f,𝐫1,𝐫2)f02(𝐤^0f,𝐫1,𝐫2)×f01(𝐫1,𝐫2)|S0f~|2(𝐤^0f,𝐫2,𝐫1)|S2~|2(𝐤^0f,𝐫2,𝐫1)σ02(𝐫1;𝐫2)×G~0m(𝐫2,𝐫1)|2|S1~|2(𝐫1,𝐫2)σ01(𝐫1;𝐫2)|S0~(𝐫1)|2|ϕ~12|2(𝐫1,𝐫2).\begin{split}\sigma&=\int\mathrm{d}\hat{\mathbf{k}}_{0f}\mathrm{d}\mathbf{r}_{% 1}\mathrm{d}\mathbf{r}_{2}\dfrac{1}{\hbar v_{0}(2\pi)^{2}}j(\hat{\mathbf{k}}_{% 0f},\mathbf{r}_{1},\mathbf{r}_{2})f_{02}(\hat{\mathbf{k}}_{0f},\mathbf{r}_{1},% \mathbf{r}_{2})\\ \times&f_{01}(\mathbf{r}_{1},\mathbf{r}_{2})|\tilde{S_{0f}}|^{2}(\hat{\mathbf{% k}}_{0f},\mathbf{r}_{2},\mathbf{r}_{1})|\tilde{S_{2}}|^{2}(\hat{\mathbf{k}}_{0% f},\mathbf{r}_{2},\mathbf{r}_{1})\sigma_{02}(\mathbf{r}_{1};\mathbf{r}_{2})\\ \times&\tilde{G}^{m}_{0}(\mathbf{r}_{2},\mathbf{r}_{1})|^{2}|\tilde{S_{1}}|^{2% }(\mathbf{r}_{1},\mathbf{r}_{2})\sigma_{01}(\mathbf{r}_{1};\mathbf{r}_{2})|% \tilde{S_{0}}(\mathbf{r}_{1})|^{2}|\tilde{\phi}_{12}|^{2}(\mathbf{r}_{1},% \mathbf{r}_{2}).\end{split}start_ROW start_CELL italic_σ end_CELL start_CELL = ∫ roman_d over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_d bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT divide start_ARG 1 end_ARG start_ARG roman_ℏ italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( 2 italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_j ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_f start_POSTSUBSCRIPT 02 end_POSTSUBSCRIPT ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL italic_f start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | over~ start_ARG italic_S start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT end_ARG | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) | over~ start_ARG italic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( over^ start_ARG bold_k end_ARG start_POSTSUBSCRIPT 0 italic_f end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_σ start_POSTSUBSCRIPT 02 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL over~ start_ARG italic_G end_ARG start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | over~ start_ARG italic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_σ start_POSTSUBSCRIPT 01 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | over~ start_ARG italic_S start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | over~ start_ARG italic_ϕ end_ARG start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) . end_CELL end_ROW (26)

II.1 Overlap function

We present briefly the calculation of the square of the overlap function |ϕ12|2superscriptsubscriptitalic-ϕ122|\phi_{12}|^{2}| italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT from nuclear structure inputs: two-nucleon amplitudes (TNA) and single-particle wavefunctions. We start from the expression from Simpson et al. (2009):

|ϕ12(𝐫1,𝐫2)|2=1Ji2MiMfΨJiMi(F)|ΨJiMi(F)=ααICαJiJfICαJiJfIDαDαI2×m1m2m1m2<j1m1j2m2|Iμ><j1m1j2m2|Iμ>×[(ϕj1m1|ϕj1m1)(ϕj2m2|ϕj2m2)+(ϕj2m2|ϕj2m2)(ϕj1m1|ϕj1m1)(ϕj1m1|ϕj2m2)(ϕj2m2|ϕj1m1)(ϕj2m2|ϕj1m1)(ϕj1m1|ϕj2m2)],superscriptsubscriptitalic-ϕ12subscript𝐫1subscript𝐫221superscriptsubscript𝐽𝑖2subscriptsubscript𝑀𝑖subscript𝑀𝑓inner-productsubscriptsuperscriptΨ𝐹subscript𝐽𝑖subscript𝑀𝑖subscriptsuperscriptΨ𝐹subscript𝐽𝑖subscript𝑀𝑖subscript𝛼superscript𝛼𝐼superscriptsubscript𝐶𝛼subscript𝐽𝑖subscript𝐽𝑓𝐼superscriptsubscript𝐶superscript𝛼subscript𝐽𝑖subscript𝐽𝑓𝐼subscript𝐷𝛼subscript𝐷superscript𝛼superscript𝐼2subscriptsubscript𝑚1subscript𝑚2superscriptsubscript𝑚1superscriptsubscript𝑚2inner-productsubscript𝑗1subscript𝑚1subscript𝑗2subscript𝑚2𝐼𝜇inner-productsuperscriptsubscript𝑗1superscriptsubscript𝑚1superscriptsubscript𝑗2superscriptsubscript𝑚2𝐼𝜇delimited-[]|subscriptsuperscriptitalic-ϕsuperscriptsubscript𝑚1superscriptsubscript𝑗1subscriptsuperscriptitalic-ϕsubscript𝑚1subscript𝑗1|subscriptsuperscriptitalic-ϕsuperscriptsubscript𝑚2superscriptsubscript𝑗2subscriptsuperscriptitalic-ϕsubscript𝑚2subscript𝑗2|subscriptsuperscriptitalic-ϕsuperscriptsubscript𝑚2superscriptsubscript𝑗2subscriptsuperscriptitalic-ϕsubscript𝑚2subscript𝑗2|subscriptsuperscriptitalic-ϕsuperscriptsubscript𝑚1superscriptsubscript𝑗1subscriptsuperscriptitalic-ϕsubscript𝑚1subscript𝑗1|subscriptsuperscriptitalic-ϕsuperscriptsubscript𝑚1superscriptsubscript𝑗1subscriptsuperscriptitalic-ϕsubscript𝑚2subscript𝑗2|subscriptsuperscriptitalic-ϕsuperscriptsubscript𝑚2superscriptsubscript𝑗2subscriptsuperscriptitalic-ϕsubscript𝑚1subscript𝑗1|subscriptsuperscriptitalic-ϕsuperscriptsubscript𝑚2superscriptsubscript𝑗2subscriptsuperscriptitalic-ϕsubscript𝑚1subscript𝑗1|subscriptsuperscriptitalic-ϕsuperscriptsubscript𝑚1superscriptsubscript𝑗1subscriptsuperscriptitalic-ϕsubscript𝑚2subscript𝑗2\begin{split}&|\phi_{12}(\mathbf{r}_{1},\mathbf{r}_{2})|^{2}=\\ &\dfrac{1}{J_{i}^{2}}\sum_{M_{i}M_{f}}\left\langle\Psi^{(F)}_{J_{i}M_{i}}|\Psi% ^{(F)}_{J_{i}M_{i}}\right\rangle=\sum_{\alpha\alpha^{\prime}I}\dfrac{C_{\alpha% }^{J_{i}J_{f}I}C_{\alpha^{\prime}}^{J_{i}J_{f}I}D_{\alpha}D_{\alpha^{\prime}}}% {I^{2}}\\ \times&\sum_{m_{1}m_{2}m_{1}^{\prime}m_{2}^{\prime}}<j_{1}m_{1}j_{2}m_{2}|I\mu% ><j_{1}^{\prime}m_{1}^{\prime}j_{2}^{\prime}m_{2}^{\prime}|I\mu>\\ \times&[(\phi^{m_{1}^{\prime}}_{j_{1}^{\prime}}|\phi^{m_{1}}_{j_{1}})(\phi^{m_% {2}^{\prime}}_{j_{2}^{\prime}}|\phi^{m_{2}}_{j_{2}})+(\phi^{m_{2}^{\prime}}_{j% _{2}^{\prime}}|\phi^{m_{2}}_{j_{2}})(\phi^{m_{1}^{\prime}}_{j_{1}^{\prime}}|% \phi^{m_{1}}_{j_{1}})\\ &-(\phi^{m_{1}^{\prime}}_{j_{1}^{\prime}}|\phi^{m_{2}}_{j_{2}})(\phi^{m_{2}^{% \prime}}_{j_{2}^{\prime}}|\phi^{m_{1}}_{j_{1}})-(\phi^{m_{2}^{\prime}}_{j_{2}^% {\prime}}|\phi^{m_{1}}_{j_{1}})(\phi^{m_{1}^{\prime}}_{j_{1}^{\prime}}|\phi^{m% _{2}}_{j_{2}})],\end{split}start_ROW start_CELL end_CELL start_CELL | italic_ϕ start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ( bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL divide start_ARG 1 end_ARG start_ARG italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∑ start_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT end_POSTSUBSCRIPT ⟨ roman_Ψ start_POSTSUPERSCRIPT ( italic_F ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT | roman_Ψ start_POSTSUPERSCRIPT ( italic_F ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT ⟩ = ∑ start_POSTSUBSCRIPT italic_α italic_α start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_I end_POSTSUBSCRIPT divide start_ARG italic_C start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT italic_I end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT italic_α start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT italic_I end_POSTSUPERSCRIPT italic_D start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT italic_D start_POSTSUBSCRIPT italic_α start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT end_ARG start_ARG italic_I start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL ∑ start_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT < italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT | italic_I italic_μ > < italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT | italic_I italic_μ > end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL [ ( italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT | italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) ( italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT | italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) + ( italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT | italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) ( italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT | italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - ( italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT | italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) ( italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT | italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) - ( italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT | italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) ( italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT | italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) ] , end_CELL end_ROW (27)

where ΨJiMi(F)subscriptsuperscriptΨ𝐹subscript𝐽𝑖subscript𝑀𝑖\Psi^{(F)}_{J_{i}M_{i}}roman_Ψ start_POSTSUPERSCRIPT ( italic_F ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_M start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT is the overlap function A|A+2inner-product𝐴𝐴2\braket{A}{A+2}⟨ start_ARG italic_A end_ARG | start_ARG italic_A + 2 end_ARG ⟩ with Ji,Misubscript𝐽𝑖subscript𝑀𝑖J_{i},M_{i}italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_M start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT and Jf,Mfsubscript𝐽𝑓subscript𝑀𝑓J_{f},M_{f}italic_J start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT , italic_M start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT being the total and magnetic angular momenta of A+2𝐴2A+2italic_A + 2 and A𝐴Aitalic_A respectively. CαJiJfIsuperscriptsubscript𝐶𝛼subscript𝐽𝑖subscript𝐽𝑓𝐼C_{\alpha}^{J_{i}J_{f}I}italic_C start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT italic_I end_POSTSUPERSCRIPT corresponds to the TNA, α=l1,j1,l2,j2𝛼subscript𝑙1subscript𝑗1subscript𝑙2subscript𝑗2\alpha=l_{1},j_{1},l_{2},j_{2}italic_α = italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, the orbital and angular momenta of the two removed protons, misubscript𝑚𝑖m_{i}italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT the corresponding magnetic angular momentum, Dα=1/2(1+δαα)subscript𝐷𝛼121subscript𝛿𝛼superscript𝛼D_{\alpha}=1/\sqrt{2(1+\delta_{\alpha\alpha^{\prime}})}italic_D start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT = 1 / square-root start_ARG 2 ( 1 + italic_δ start_POSTSUBSCRIPT italic_α italic_α start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ) end_ARG, 𝐈=𝐣𝟏+𝐣𝟐𝐈subscript𝐣1subscript𝐣2\mathbf{I}=\mathbf{j_{1}}+\mathbf{j_{2}}bold_I = bold_j start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT + bold_j start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT and the bracket (ϕj2m2|ϕj2m2)conditionalsubscriptsuperscriptitalic-ϕsuperscriptsubscript𝑚2superscriptsubscript𝑗2subscriptsuperscriptitalic-ϕsubscript𝑚2subscript𝑗2(\phi^{m_{2}^{\prime}}_{j_{2}^{\prime}}|\phi^{m_{2}}_{j_{2}})( italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT | italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) corresponds to:

(ϕj1m1|ϕj1m1)=λ1λ1σ<l1λ1sσ|j1m1><l1λ1sσ|j1m1>ul1j1(r1)Yl1λ1(𝐫𝟏^)ul1j1(r1)Yl1λ1(𝐫𝟏^),conditionalsubscriptsuperscriptitalic-ϕsuperscriptsubscript𝑚1superscriptsubscript𝑗1subscriptsuperscriptitalic-ϕsubscript𝑚1subscript𝑗1subscriptsubscript𝜆1superscriptsubscript𝜆1𝜎inner-productsubscript𝑙1subscript𝜆1𝑠𝜎subscript𝑗1subscript𝑚1inner-productsuperscriptsubscript𝑙1superscriptsubscript𝜆1𝑠𝜎superscriptsubscript𝑗1superscriptsubscript𝑚1subscript𝑢subscript𝑙1subscript𝑗1subscript𝑟1subscript𝑌subscript𝑙1subscript𝜆1^subscript𝐫1subscriptsuperscript𝑢superscriptsubscript𝑙1superscriptsubscript𝑗1subscript𝑟1subscriptsuperscript𝑌superscriptsubscript𝑙1superscriptsubscript𝜆1^subscript𝐫1\begin{split}(\phi^{m_{1}^{\prime}}_{j_{1}^{\prime}}|\phi^{m_{1}}_{j_{1}})=&% \sum_{\lambda_{1}\lambda_{1}^{\prime}\sigma}<l_{1}\lambda_{1}s\sigma|j_{1}m_{1% }><l_{1}^{\prime}\lambda_{1}^{\prime}s\sigma|j_{1}^{\prime}m_{1}^{\prime}>\\ &u_{l_{1}j_{1}}(r_{1})Y_{l_{1}\lambda_{1}}(\hat{\mathbf{r_{1}}})u^{*}_{l_{1}^{% \prime}j_{1}^{\prime}}(r_{1})Y^{*}_{l_{1}^{\prime}\lambda_{1}^{\prime}}(\hat{% \mathbf{r_{1}}}),\end{split}start_ROW start_CELL ( italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT | italic_ϕ start_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) = end_CELL start_CELL ∑ start_POSTSUBSCRIPT italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_σ end_POSTSUBSCRIPT < italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_s italic_σ | italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT > < italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_s italic_σ | italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT > end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL italic_u start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_Y start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( over^ start_ARG bold_r start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT end_ARG ) italic_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_Y start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( over^ start_ARG bold_r start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT end_ARG ) , end_CELL end_ROW (28)

where uljsubscript𝑢𝑙𝑗u_{lj}italic_u start_POSTSUBSCRIPT italic_l italic_j end_POSTSUBSCRIPT is the radial wavefunction, evaluated at S2p/2subscript𝑆2𝑝2S_{2p}/2italic_S start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT / 2 Simpson et al. (2009), Ylλsubscript𝑌𝑙𝜆Y_{l\lambda}italic_Y start_POSTSUBSCRIPT italic_l italic_λ end_POSTSUBSCRIPT the spherical harmonic with magnetic angular momentum λ𝜆\lambdaitalic_λ and s𝑠sitalic_s and σ𝜎\sigmaitalic_σ are the spin of the proton and its projection. Note that in Eq. (27), the first bracket in each term is evaluated for 𝐫1subscript𝐫1\mathbf{r}_{1}bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT while the second bracket is evaluated for 𝐫2subscript𝐫2\mathbf{r}_{2}bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT.

After some angular momentum algebra, we reach the following expression:

|ϕ~12|2(𝐫𝟏,𝐫𝟐)=ααICαJiJfICαJiJfIDαDαl1^l1^l2^l2^j1^j2^j1^j2^16π2×()s1s2+j2+j2LPL(𝐫𝟏𝐫𝟐)×[(()I2j2<l10,l10|L0><l20,l20|L0>{l1s1j1j1Ll1}×{l2s2j2j2Ll2}{j1Ij2j2Lj1})×(ul1j1(r1)ul1j1(r1)ul2j2(r2)ul2j2(r2))+×+ul2j2(r1)ul2j2(r1)ul1j1(r2)ul1j1(r2))(<l10,l20|L0><l20,l10|L0>{l1s1j1j2Ll2}{l2s2j2j1Ll1}×{j1Ij2j1Lj2})(ul2j2(r1)ul1j1(r1)ul1j1(r2)ul2j2(r2))+×+ul1j1(r1)ul2j2(r1)ul2j2(r2)ul1j1(r2))],\begin{split}&|\tilde{\phi}_{12}|^{2}(\mathbf{r_{1}},\mathbf{r_{2}})=\sum_{% \alpha\alpha^{\prime}I}C_{\alpha}^{J_{i}J_{f}I}C_{\alpha^{\prime}}^{J_{i}J_{f}% I}D_{\alpha}D_{\alpha^{\prime}}\dfrac{\hat{l_{1}}\hat{l_{1}^{\prime}}\hat{l_{2% }}\hat{l_{2}^{\prime}}\hat{j_{1}}\hat{j_{2}}\hat{j_{1}^{\prime}}\hat{j_{2}^{% \prime}}}{16\pi^{2}}\\ \times&(-)^{-s_{1}-s_{2}+j_{2}+j_{2}^{\prime}}\sum_{L}P^{L}(\mathbf{r_{1}}% \cdot\mathbf{r_{2}})\\ \times&\left[\left((-)^{I-2j_{2}^{\prime}}<l_{1}0,l_{1}^{\prime}0|L0><l_{2}0,l% _{2}^{\prime}0|L0>\begin{Bmatrix}l_{1}&s_{1}&j_{1}\\ j_{1}^{\prime}&L&l_{1}^{\prime}\end{Bmatrix}\right.\right.\\ \times&\left.\begin{Bmatrix}l_{2}&s_{2}&j_{2}\\ j_{2}^{\prime}&L&l_{2}^{\prime}\end{Bmatrix}\begin{Bmatrix}j_{1}&I&j_{2}\\ j_{2}^{\prime}&L&j_{1}^{\prime}\end{Bmatrix}\right)\\ \times&(u_{l_{1}j_{1}}(r_{1})u^{*}_{l_{1}^{\prime}j_{1}^{\prime}}(r_{1})u_{l_{% 2}j_{2}}(r_{2})u^{*}_{l_{2}^{\prime}j_{2}^{\prime}}(r_{2}))+\\ \times&+u_{l_{2}j_{2}}(r_{1})u^{*}_{l_{2}^{\prime}j_{2}^{\prime}}(r_{1})u_{l_{% 1}j_{1}}(r_{2})u^{*}_{l_{1}^{\prime}j_{1}^{\prime}}(r_{2}))-\\ -&\left(<l_{1}0,l_{2}^{\prime}0|L0><l_{2}0,l_{1}^{\prime}0|L0>\begin{Bmatrix}l% _{1}&s_{1}&j_{1}\\ j_{2}^{\prime}&L&l_{2}^{\prime}\end{Bmatrix}\begin{Bmatrix}l_{2}&s_{2}&j_{2}\\ j_{1}^{\prime}&L&l_{1}^{\prime}\end{Bmatrix}\right.\\ \times&\left.\begin{Bmatrix}j_{1}&I&j_{2}\\ j_{1}^{\prime}&L&j_{2}^{\prime}\end{Bmatrix}\right)(u_{l_{2}j_{2}}(r_{1})u^{*}% _{l_{1}^{\prime}j_{1}^{\prime}}(r_{1})u_{l_{1}j_{1}}(r_{2})u^{*}_{l_{2}^{% \prime}j_{2}^{\prime}}(r_{2}))+\\ \times&+\left.u_{l_{1}j_{1}}(r_{1})u^{*}_{l_{2}^{\prime}j_{2}^{\prime}}(r_{1})% u_{l_{2}j_{2}}(r_{2})u^{*}_{l_{1}^{\prime}j_{1}^{\prime}}(r_{2}))\right],\end{split}start_ROW start_CELL end_CELL start_CELL | over~ start_ARG italic_ϕ end_ARG start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT , bold_r start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT ) = ∑ start_POSTSUBSCRIPT italic_α italic_α start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_I end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT italic_I end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT italic_α start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_J start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT italic_I end_POSTSUPERSCRIPT italic_D start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT italic_D start_POSTSUBSCRIPT italic_α start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT divide start_ARG over^ start_ARG italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG over^ start_ARG italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG over^ start_ARG italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG over^ start_ARG italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG over^ start_ARG italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG over^ start_ARG italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG over^ start_ARG italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG over^ start_ARG italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG 16 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL ( - ) start_POSTSUPERSCRIPT - italic_s start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_s start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT italic_P start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT ⋅ bold_r start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT ) end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL [ ( ( - ) start_POSTSUPERSCRIPT italic_I - 2 italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT < italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT 0 , italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT 0 | italic_L 0 > < italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT 0 , italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT 0 | italic_L 0 > { start_ARG start_ROW start_CELL italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL italic_s start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL start_CELL italic_L end_CELL start_CELL italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL end_ROW end_ARG } end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL { start_ARG start_ROW start_CELL italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL start_CELL italic_s start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL start_CELL italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL start_CELL italic_L end_CELL start_CELL italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL end_ROW end_ARG } { start_ARG start_ROW start_CELL italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL italic_I end_CELL start_CELL italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL start_CELL italic_L end_CELL start_CELL italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL end_ROW end_ARG } ) end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL ( italic_u start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_u start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ) + end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL + italic_u start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_u start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ) - end_CELL end_ROW start_ROW start_CELL - end_CELL start_CELL ( < italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT 0 , italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT 0 | italic_L 0 > < italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT 0 , italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT 0 | italic_L 0 > { start_ARG start_ROW start_CELL italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL italic_s start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL start_CELL italic_L end_CELL start_CELL italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL end_ROW end_ARG } { start_ARG start_ROW start_CELL italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL start_CELL italic_s start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL start_CELL italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL start_CELL italic_L end_CELL start_CELL italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL end_ROW end_ARG } end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL { start_ARG start_ROW start_CELL italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL italic_I end_CELL start_CELL italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL start_CELL italic_L end_CELL start_CELL italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL end_ROW end_ARG } ) ( italic_u start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_u start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ) + end_CELL end_ROW start_ROW start_CELL × end_CELL start_CELL + italic_u start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_u start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_u start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ) ] , end_CELL end_ROW (29)

where we note that the angular dependence on 𝐫^1subscript^𝐫1\hat{\mathbf{r}}_{1}over^ start_ARG bold_r end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and 𝐫^2subscript^𝐫2\hat{\mathbf{r}}_{2}over^ start_ARG bold_r end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is reduced to the Legendre polynomial PL(𝐫𝟏𝐫𝟐)superscript𝑃𝐿subscript𝐫1subscript𝐫2P^{L}(\mathbf{r_{1}}\cdot\mathbf{r_{2}})italic_P start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT ( bold_r start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT ⋅ bold_r start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT ), which allows for more efficient treatment of the overlap wavefunction, which can be stored as a function of r1subscript𝑟1r_{1}italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, r2subscript𝑟2r_{2}italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (only the modulus) and L𝐿Litalic_L, instead of 𝐫1subscript𝐫1\mathbf{r}_{1}bold_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and 𝐫2subscript𝐫2\mathbf{r}_{2}bold_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT.

III Numerical results

For all the following calculations the optical potential between proton and nucleus has been taken from the global Dirac parametrization Hama et al. (1990); Cooper et al. (1993), while the proton-proton elastic cross section has been taken from the parametrization by Bertulani and De Conti Bertulani and De Conti (2010). The single-particle wavefunctions have been computed using Woods-Saxon potentials with diffusivity a=0.7𝑎0.7a=0.7italic_a = 0.7 fm and a radius adjusted to reproduce the rms radius from Hartree-Fock calculations using the SkX interaction Alex Brown (1998). The integration over r1subscript𝑟1r_{1}italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and r2subscript𝑟2r_{2}italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is extended up to the radius where the single-particle wavefunctions is reduced to 1/1000110001/10001 / 1000 of their maximum value (the maximum radius among the wavefunctions) and the step of integration is taken as 0.2 fm, which was found to give converged results to similar-to\sim5%.

III.1 C12(p,3p)10BesuperscriptC12superscript𝑝3𝑝10Be{}^{12}\mathrm{C}(p,3p)^{10}\mathrm{Be}start_FLOATSUPERSCRIPT 12 end_FLOATSUPERSCRIPT roman_C ( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 10 end_POSTSUPERSCRIPT roman_Be

We first present a study of the 12C(p,3p)10superscript𝑝3𝑝10(p,3p)^{10}( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 10 end_POSTSUPERSCRIPTBe reaction. Admittedly, the 12C nucleus is rather light while the formalism developed in this work is meant for heavier systems. An estimation of the effect of the finite mass of the nucleus can be obtained from Kawai and Weidenmüller (1992), which presents a very similar result to ours in the context of the study of the total nucleon-nucleus cross section. In Kawai and Weidenmüller (1992), the cross section is scaled by a factor of (AA+1)4superscript𝐴𝐴14\left(\frac{A}{A+1}\right)^{4}( divide start_ARG italic_A end_ARG start_ARG italic_A + 1 end_ARG ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT, which in our case would introduce a factor of 0.7 in the cross sections. We believe the ambiguities in the optical potentials may introduce larger uncertainties in the total cross section, so we will not consider this factor in the following. The nuclear structure of 12C is well known and there exist experimental data for fragmentation of 12C on proton targets to produce 10Be Olson et al. (1983), so we consider this reaction a reasonable one to benchmark our calculations.

The TNA for the calculation have been obtained using the WBT interaction Warburton and Brown (1992) in the shell model code oxbash Etchegoyen et al. (1985) (calculations using the interaction by Cohen and Kurath Cohen and Kurath (1965) show consistent TNA). Those TNA deemed too small have been excluded from the calculation. In Table 1 the used TNA are presented. In order to compare with the total experimental cross section in Olson et al. (1983), the calculation has been performed at 1.05 GeV/A. Only states under the neutron emission threshold for 10Be have been considered and their energies have been taken as the experimental ones.

Jfπsubscriptsuperscript𝐽𝜋𝑓J^{\pi}_{f}italic_J start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT Exsubscript𝐸𝑥E_{x}italic_E start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT (MeV) [1s1/2]2superscriptdelimited-[]1subscript𝑠122[1s_{1/2}]^{2}[ 1 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [1p3/2]2superscriptdelimited-[]1subscript𝑝322[1p_{3/2}]^{2}[ 1 italic_p start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [1p1/21p3/2]delimited-[]1subscript𝑝121subscript𝑝32[1p_{1/2}1p_{3/2}][ 1 italic_p start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT 1 italic_p start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT ] [1p1/2]2superscriptdelimited-[]1subscript𝑝122[1p_{1/2}]^{2}[ 1 italic_p start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [1d5/2]2superscriptdelimited-[]1subscript𝑑522[1d_{5/2}]^{2}[ 1 italic_d start_POSTSUBSCRIPT 5 / 2 end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [1d3/2]2superscriptdelimited-[]1subscript𝑑322[1d_{3/2}]^{2}[ 1 italic_d start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [2s1/2]2superscriptdelimited-[]2subscript𝑠122[2s_{1/2}]^{2}[ 2 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [2s1/21d5/2]delimited-[]2subscript𝑠121subscript𝑑52[2s_{1/2}1d_{5/2}][ 2 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT 1 italic_d start_POSTSUBSCRIPT 5 / 2 end_POSTSUBSCRIPT ] [2s1/21d3/2]delimited-[]2subscript𝑠121subscript𝑑32[2s_{1/2}1d_{3/2}][ 2 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT 1 italic_d start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT ]
0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 0.00 1.461 0.706 -0.064 -0.056 -0.056
2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 3.37 2.060 -0.854 0.038 0.020
2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 5.96 -0.419 1.204 0.028 0.016 -0.020 -0.010
0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 6.18 -0.018 -0.165 0.145 0.016
Jfπsubscriptsuperscript𝐽𝜋𝑓J^{\pi}_{f}italic_J start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT Exsubscript𝐸𝑥E_{x}italic_E start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT (MeV) [1p3/21s1/2]delimited-[]1subscript𝑝321subscript𝑠12[1p_{3/2}1s_{1/2}][ 1 italic_p start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT 1 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT ] [1p1/21s1/2]delimited-[]1subscript𝑝121subscript𝑠12[1p_{1/2}1s_{1/2}][ 1 italic_p start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT 1 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT ] [1d5/21p3/2]delimited-[]1subscript𝑑521subscript𝑝32[1d_{5/2}1p_{3/2}][ 1 italic_d start_POSTSUBSCRIPT 5 / 2 end_POSTSUBSCRIPT 1 italic_p start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT ] [1d5/21p1/2]delimited-[]1subscript𝑑521subscript𝑝12[1d_{5/2}1p_{1/2}][ 1 italic_d start_POSTSUBSCRIPT 5 / 2 end_POSTSUBSCRIPT 1 italic_p start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT ] [1d3/21p3/2]delimited-[]1subscript𝑑321subscript𝑝32[1d_{3/2}1p_{3/2}][ 1 italic_d start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT 1 italic_p start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT ] [1d3/21p1/2]delimited-[]1subscript𝑑321subscript𝑝12[1d_{3/2}1p_{1/2}][ 1 italic_d start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT 1 italic_p start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT ] [2s1/21p3/2]delimited-[]2subscript𝑠121subscript𝑝32[2s_{1/2}1p_{3/2}][ 2 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT 1 italic_p start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT ] [2s1/21p1/2]delimited-[]2subscript𝑠121subscript𝑝12[2s_{1/2}1p_{1/2}][ 2 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT 1 italic_p start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT ]
1superscript11^{-}1 start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT 5.96 -0.247 0.135 -0.117 0.069 -0.077 0.065 0.051
2superscript22^{-}2 start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT 6.26 0.045 -0.048 -0.055 0.105 -0.021
Table 1: TNA used for the 12C(p,3p)10superscript𝑝3𝑝10(p,3p)^{10}( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 10 end_POSTSUPERSCRIPTBe reaction, not including the isospin Clebsch-Gordan factor, obtained using the WBT interaction.

The results are presented in Table. 2. We must note that since the kinematics of the removed protons was not measured in the experiment, another process that could result in 10Be is a (p,2p)𝑝2𝑝(p,2p)( italic_p , 2 italic_p ) reaction producing an excited B11superscriptsuperscriptB11{}^{11}\mathrm{B}^{*}start_FLOATSUPERSCRIPT 11 end_FLOATSUPERSCRIPT roman_B start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT which then decays by emission of a proton and thus competes with the process studied here Olson et al. (1983). In a recent measurement of the 12C(p,2p)𝑝2𝑝(p,2p)( italic_p , 2 italic_p ) reaction at 398 MeV/A, the cross section for the 12C(p,2p)11superscript𝑝2𝑝11(p,2p)^{11}( italic_p , 2 italic_p ) start_POSTSUPERSCRIPT 11 end_POSTSUPERSCRIPTB process was found to be 18.1(2.0) mb while the cross section for the decay of 11B into 10Be (distinguished by imposing a coincidence of a forward proton and 10Be) was 0.8(0.3) Panin (2012). Assuming the same proportion for the reaction at 1.05 GeV/A, given that the 12C(p,2p)11superscript𝑝2𝑝11(p,2p)^{11}( italic_p , 2 italic_p ) start_POSTSUPERSCRIPT 11 end_POSTSUPERSCRIPTB cross section was of 30.9(3.4) mb, this would yield a cross section for excitation-evaporation of 1.4(0.6) mb, leaving 2.0(0.7) mb for the direct removal of the two protons studied in this work, which gives reasonable agreement with the result of our calculations, and thus we consider this result a validation of this work. It should be remarked that in Tostevin et al. (2004b); Simpson and Tostevin (2011), for two-proton removal from 12C using a carbon target at the same beam energy, the excitation-evaporation process was not removed and good agreement between theory and experiment was found still, which could indicate differences between proton and carbon targets when populating proton-unbound excited states of 11B.

Jfπsubscriptsuperscript𝐽𝜋𝑓J^{\pi}_{f}italic_J start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT Exsubscript𝐸𝑥E_{x}italic_E start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT (MeV) σ𝜎\sigmaitalic_σ (mb)
0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 0.00 0.86
2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 3.37 0.68
2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 5.96 0.27
1superscript11^{-}1 start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT 5.96 9.21039.2superscript1039.2\cdot 10^{-3}9.2 ⋅ 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT
0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 6.18 4.11034.1superscript1034.1\cdot 10^{-3}4.1 ⋅ 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT
2superscript22^{-}2 start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT 6.26 5.51045.5superscript1045.5\cdot 10^{-4}5.5 ⋅ 10 start_POSTSUPERSCRIPT - 4 end_POSTSUPERSCRIPT
Total 1.82
Exp.Olson et al. (1983) 3.41(0.54)
Exp. (direct (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p )) 2.0(0.7)*
Table 2: Cross sections for the 12C(p,3p)10superscript𝑝3𝑝10(p,3p)^{10}( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 10 end_POSTSUPERSCRIPTBe(Jfπ)subscriptsuperscript𝐽𝜋𝑓(J^{\pi}_{f})( italic_J start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ) reaction at 1.05 GeV/A
See text

III.2 Mg28(p,3p)26NesuperscriptMg28superscript𝑝3𝑝26Ne{}^{28}\mathrm{Mg}(p,3p)^{26}\mathrm{Ne}start_FLOATSUPERSCRIPT 28 end_FLOATSUPERSCRIPT roman_Mg ( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 26 end_POSTSUPERSCRIPT roman_Ne

In this section, we explore the Mg28(p,3p)26NesuperscriptMg28superscript𝑝3𝑝26Ne{}^{28}\mathrm{Mg}(p,3p)^{26}\mathrm{Ne}start_FLOATSUPERSCRIPT 28 end_FLOATSUPERSCRIPT roman_Mg ( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 26 end_POSTSUPERSCRIPT roman_Ne reaction at 250 MeV/A. Although no experimental data exists for this reaction, the equivalent two-nucleon knockout reaction with beryllium target has been measured Bazin et al. (2003) and thoroughly studied Tostevin et al. (2004a, b); Simpson and Tostevin (2010), so we find this reaction is a good benchmark to compare the (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) and two-proton knockout reactions when populating different excited states. The structure inputs are the same as those from Tostevin et al. (2004a), the TNA are presented in Table 3 for completeness.

Jfπsubscriptsuperscript𝐽𝜋𝑓J^{\pi}_{f}italic_J start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT Exsubscript𝐸𝑥E_{x}italic_E start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT (MeV) [1d3/2]2superscriptdelimited-[]1subscript𝑑322[1d_{3/2}]^{2}[ 1 italic_d start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [1d3/21d5/2]delimited-[]1subscript𝑑321subscript𝑑52[1d_{3/2}1d_{5/2}][ 1 italic_d start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT 1 italic_d start_POSTSUBSCRIPT 5 / 2 end_POSTSUBSCRIPT ] [1d5/2]2superscriptdelimited-[]1subscript𝑑522[1d_{5/2}]^{2}[ 1 italic_d start_POSTSUBSCRIPT 5 / 2 end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [2s1/21d3/2]delimited-[]2subscript𝑠121subscript𝑑32[2s_{1/2}1d_{3/2}][ 2 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT 1 italic_d start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT ] [2s1/21d5/2]delimited-[]2subscript𝑠121subscript𝑑52[2s_{1/2}1d_{5/2}][ 2 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT 1 italic_d start_POSTSUBSCRIPT 5 / 2 end_POSTSUBSCRIPT ] [2s1/2]2superscriptdelimited-[]2subscript𝑠122[2s_{1/2}]^{2}[ 2 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT
0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 0.00 -0.301 -1.047 -0.305
2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 2.02 -0.050 0.374 -0.637 -0.061 -0.139
4+superscript44^{+}4 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 3.50 0.331 1.596
2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 3.70 0.047 -0.072 0.853 0.161 0.176
Table 3: TNA used for the 28Mg(p,3p)26superscript𝑝3𝑝26(p,3p)^{26}( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 26 end_POSTSUPERSCRIPTNe reaction, taken from Tostevin et al. (2004a).

We are interested in whether (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) and two-proton knockout reactions populate equally the excited states of the residual nucleus, as such we present in Table 4 the cross sections for the four considered states for the (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) reaction at 250 MeV/A, a typical energy for radioactive beam facilities, and the two-proton knockout reaction at 82.3 MeV/A Bazin et al. (2003), as well as their ratio to the total cross section. Experimental results for the two-proton knockout reaction are presented as well. The theoretical and experimental results for the two-proton knockout reaction are taken from Simpson and Tostevin (2010) and Tostevin et al. (2004a) respectively.

In general, we find the distribution of the cross section is similar for both reaction, with the ground and 4+superscript44^{+}4 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT excited states taking most of the cross section, although it should be noted that the contribution of the ground state is significantly larger for the (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) reaction than for the two-proton knockout one. To further explore this discrepancy we have performed a calculation in which the interference between configurations with different TNA has been turned off in order to evaluate the effect of this interference in the cross section. The results are shown on the third column of Table. 4, while the fourth column shows the ratio between these calculations and the original ones. It is noteworthy that it is the ground state the one that presents most sensitivity to this interference, suggesting that it lies behind the difference between (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) and knockout, possibly because (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) is more sensitive to the nuclear interior, which presents a different behaviour than the nuclear surface in respect to the interference of different configurations. This suggests that interference between configurations and correlations between the removed protons plays a significant role in two-proton-removal reactions Simpson and Tostevin (2010, 2011), and that (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) and two-proton knockout reactions may be differently sensitive to them. The last two columns of Table 4 show the sum of the square of the TNA and their ratio to the total for each Ne26superscriptNe26{}^{26}\mathrm{Ne}start_FLOATSUPERSCRIPT 26 end_FLOATSUPERSCRIPT roman_Ne final state, to check whether a direct relation between the reaction cross section and the structure observables, akin to spectroscopic factors for one-nucleon removal reactions, can be established for two-nucleon removal ones. As expected, this is not the case, as the contribution of each final state to the total cross section difers significantly from its “weight”, measured as the sum of the square of the TNAs. Therefore the reaction theory is essential to be able to predict the relative population of each final state, although a rough relation can be established between TNA and cross sections, in that states with a small sum of TNA2𝑇𝑁superscript𝐴2TNA^{2}italic_T italic_N italic_A start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT can be expected to be weakly-populated and a state with a large sum is likely to be significantly populated.

Jfπsubscriptsuperscript𝐽𝜋𝑓J^{\pi}_{f}italic_J start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT Exsubscript𝐸𝑥E_{x}italic_E start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT (MeV) σ(p,3p)𝜎𝑝3𝑝\sigma(p,3p)italic_σ ( italic_p , 3 italic_p ) (mb) Ratio σ(p,3p)𝜎𝑝3𝑝\sigma(p,3p)italic_σ ( italic_p , 3 italic_p ) no interf. (mb) (p,3p)nointerf.(p,3p)superscript𝑝3𝑝nointerf𝑝3𝑝\dfrac{(p,3p)^{\mathrm{nointerf.}}}{(p,3p)}divide start_ARG ( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT roman_nointerf . end_POSTSUPERSCRIPT end_ARG start_ARG ( italic_p , 3 italic_p ) end_ARG σ2p𝜎2𝑝\sigma 2pitalic_σ 2 italic_pKO (mb) Ratio Exp 2p2𝑝2p2 italic_pKO (mb) TNA2superscriptTNA2\sum\mathrm{TNA}^{2}∑ roman_TNA start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT Ratio
0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 0.00 0.274 0.48 0.196 0.71 1.19 0.40 0.70(15) 1.28 0.24
2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 2.02 0.048 0.08 0.048 1.00 0.32 0.11 0.09(15) 0.69 0.13
4+superscript44^{+}4 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 3.50 0.181 0.32 0.160 0.88 1.02 0.34 0.58(9) 2.66 0.49
2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 3.70 0.063 0.11 0.061 0.97 0.45 0.15 0.15(9) 0.79 0.15
Total 0.57 0.46 2.98 1.50(10) 5.42
Table 4: Cross sections for the 28Mg(p,3p)26superscript𝑝3𝑝26(p,3p)^{26}( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 26 end_POSTSUPERSCRIPTNe(Jfπ)subscriptsuperscript𝐽𝜋𝑓(J^{\pi}_{f})( italic_J start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ) reaction at 250 MeV/A (3rd column) and the9Be(28(^{28}( start_POSTSUPERSCRIPT 28 end_POSTSUPERSCRIPTMg,26Ne(Jfπ))X(J^{\pi}_{f}))X( italic_J start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ) ) italic_X reaction at 82.3 MeV/A (7th column Simpson and Tostevin (2010), 9th column Tostevin et al. (2004a)), and the ratio of each 26Ne final state to the total cross section (2nd and 6th column respectively). The 3rd column corresponds to calculations where configurations with different TNA are not allowed to interfere and the 4th column is the ratio between this calculation and the full one. The 10th column corresponds to the sum of the square of the TNA considered for each state, and the 11th column to its ratio to the total adding all final states.

III.3 Ca54(p,3p)52ArsuperscriptCa54superscript𝑝3𝑝52Ar{}^{54}\mathrm{Ca}(p,3p)^{52}\mathrm{Ar}start_FLOATSUPERSCRIPT 54 end_FLOATSUPERSCRIPT roman_Ca ( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 52 end_POSTSUPERSCRIPT roman_Ar

Although there is a number of medium-mass nuclei for which the (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) reaction has been measured Frotscher et al. (2020), for these nuclei the nuclear structure calculations are computationally heavy and the experimental data scarce, so there are significant uncertainties in the TNA (admittedly, also in the optical potentials). As such, it is challenging to find neutron-rich nuclei for which to compare experimental cross sections to the results of this work. For the Ca54(p,3p)52ArsuperscriptCa54superscript𝑝3𝑝52Ar{}^{54}\mathrm{Ca}(p,3p)^{52}\mathrm{Ar}start_FLOATSUPERSCRIPT 54 end_FLOATSUPERSCRIPT roman_Ca ( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 52 end_POSTSUPERSCRIPT roman_Ar reaction, we were facilitated unpublished experimental data Liu at 250 MeV/A, which yielded a total cross section of σ=0.047(6)𝜎0.0476\sigma=0.047(6)italic_σ = 0.047 ( 6 ) mb. For <54Ca|52Ar><^{54}\mathrm{Ca}|^{52}\mathrm{Ar}>< start_POSTSUPERSCRIPT 54 end_POSTSUPERSCRIPT roman_Ca | start_POSTSUPERSCRIPT 52 end_POSTSUPERSCRIPT roman_Ar >, TNA have been graciously provided by Prof. Y. Utsuno Utsuno , and have been obtained using an extension of the GXPF1Br Steppenbeck et al. (2013) interaction extended to the sd-pf-sdg space. The ones considered in these calculations are presented in Table 5. When considering which states to include, a significant strength can be found close to the separation energy of the neutron Snsubscript𝑆𝑛S_{n}italic_S start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT in 52Ar. Taking the values from NuDat nud , we find a value for Snsubscript𝑆𝑛S_{n}italic_S start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT of 3.0(7) MeV, obtained by a systematic fit to the nuclear mass, so there is an ambiguity on whether the population of states with excitation energy Ex=3.03.7subscript𝐸𝑥3.03.7E_{x}=3.0-3.7italic_E start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = 3.0 - 3.7 MeV would lead to a bound residual 52Ar (thus contributing to the cross section), or not.

Jfπsubscriptsuperscript𝐽𝜋𝑓J^{\pi}_{f}italic_J start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT Exsubscript𝐸𝑥E_{x}italic_E start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT (MeV) [1d5/2]2superscriptdelimited-[]1subscript𝑑522[1d_{5/2}]^{2}[ 1 italic_d start_POSTSUBSCRIPT 5 / 2 end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [1d5/21d3/2]delimited-[]1subscript𝑑521subscript𝑑32[1d_{5/2}1d_{3/2}][ 1 italic_d start_POSTSUBSCRIPT 5 / 2 end_POSTSUBSCRIPT 1 italic_d start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT ] [1d5/22s1/2]delimited-[]1subscript𝑑522subscript𝑠12[1d_{5/2}2s_{1/2}][ 1 italic_d start_POSTSUBSCRIPT 5 / 2 end_POSTSUBSCRIPT 2 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT ] [1d3/2]2superscriptdelimited-[]1subscript𝑑322[1d_{3/2}]^{2}[ 1 italic_d start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [1d3/22s1/2]delimited-[]1subscript𝑑322subscript𝑠12[1d_{3/2}2s_{1/2}][ 1 italic_d start_POSTSUBSCRIPT 3 / 2 end_POSTSUBSCRIPT 2 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT ] [2s1/2]2superscriptdelimited-[]2subscript𝑠122[2s_{1/2}]^{2}[ 2 italic_s start_POSTSUBSCRIPT 1 / 2 end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT σp3psubscript𝜎𝑝3𝑝\sigma_{p3p}italic_σ start_POSTSUBSCRIPT italic_p 3 italic_p end_POSTSUBSCRIPT (mb) Rssubscript𝑅𝑠R_{s}italic_R start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT
0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 0.00 -0.221 -0.839 -0.222 3.51102absentsuperscript102\cdot 10^{-2}⋅ 10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT
2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 1.719 -0.205 -0.238 -0.277 -1.187 0.650 4.82102absentsuperscript102\cdot 10^{-2}⋅ 10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT
0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 2.356 0.053 0.277 -0.030 2.75103absentsuperscript103\cdot 10^{-3}⋅ 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT
2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 2.364 0.118 0.074 0.192 0.575 -0.492 1.54102absentsuperscript102\cdot 10^{-2}⋅ 10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT
3+superscript33^{+}3 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 2.809 0.029 -0.013 1.35105absentsuperscript105\cdot 10^{-5}⋅ 10 start_POSTSUPERSCRIPT - 5 end_POSTSUPERSCRIPT
2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 3.012 0.167 0.077 0.052 1.378 -0.147 3.32102absentsuperscript102\cdot 10^{-2}⋅ 10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT
3+superscript33^{+}3 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 3.363 -0.077 -0.013 7.62105absentsuperscript105\cdot 10^{-5}⋅ 10 start_POSTSUPERSCRIPT - 5 end_POSTSUPERSCRIPT
4+superscript44^{+}4 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 3.464 0.073 0.318 2.57103absentsuperscript103\cdot 10^{-3}⋅ 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT
4+superscript44^{+}4 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 3.592 0.0511 0.355 3.00103absentsuperscript103\cdot 10^{-3}⋅ 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT
2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 3.639 -0.017 0.108 0.263 -0.678 -0.984 1.89102absentsuperscript102\cdot 10^{-2}⋅ 10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT
0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT 3.670 0.001 0.089 -0.096 2.97104absentsuperscript104\cdot 10^{-4}⋅ 10 start_POSTSUPERSCRIPT - 4 end_POSTSUPERSCRIPT
Total (Exmax=Sn=3.0superscriptsubscript𝐸𝑥maxsubscript𝑆𝑛3.0E_{x}^{\mathrm{max}}=S_{n}=3.0italic_E start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_max end_POSTSUPERSCRIPT = italic_S start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = 3.0 MeV) 1.021011.02superscript1011.02\cdot 10^{-1}1.02 ⋅ 10 start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT 0.46
Total (Exmax=Sn+σSn=3.7superscriptsubscript𝐸𝑥maxsubscript𝑆𝑛subscript𝜎subscript𝑆𝑛3.7E_{x}^{\mathrm{max}}=S_{n}+\sigma_{S_{n}}=3.7italic_E start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_max end_POSTSUPERSCRIPT = italic_S start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT + italic_σ start_POSTSUBSCRIPT italic_S start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 3.7 MeV) 1.591011.59superscript1011.59\cdot 10^{-1}1.59 ⋅ 10 start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT 0.29
Table 5: States, TNA and cross sections considered for the 54Ca(p,3p)52superscript𝑝3𝑝52(p,3p)^{52}( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 52 end_POSTSUPERSCRIPTAr reaction. The last column corresponds to the reduction factor Rs=σexp/σthsubscript𝑅𝑠subscript𝜎expsubscript𝜎𝑡R_{s}=\sigma_{\mathrm{exp}}/\sigma_{th}italic_R start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = italic_σ start_POSTSUBSCRIPT roman_exp end_POSTSUBSCRIPT / italic_σ start_POSTSUBSCRIPT italic_t italic_h end_POSTSUBSCRIPT. States above the double horizontal line lie below the nominal value of Snsubscript𝑆𝑛S_{n}italic_S start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT for 52Ar, while those below lie below its value plus one σ𝜎\sigmaitalic_σ nud .

To check the effect of this ambiguity we have computed the cross section to all states for energies up to 3.73.73.73.7 MeV, and present the cross sections adding up the states up to Ex=3.0subscript𝐸𝑥3.0E_{x}=3.0italic_E start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = 3.0 MeV (which corresponds to the nominal value of Snsubscript𝑆𝑛S_{n}italic_S start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT) and up to Ex=3.7subscript𝐸𝑥3.7E_{x}=3.7italic_E start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = 3.7 MeV (which would also include the states within the 1σ𝜎\sigmaitalic_σ error range). Both values yield results of 0.1020.1020.1020.102 mb and 0.1590.1590.1590.159 mb, which shows that even the uncertainty in Snsubscript𝑆𝑛S_{n}italic_S start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT produces differences of 50%percent5050\%50 % in the results. Better measurements of the masses of medium-mass nuclei or measurements of (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) reactions with gamma coincidence would help alleviate this ambiguity. As for the reduction factors, we find values of Rs=0.46subscript𝑅𝑠0.46R_{s}=0.46italic_R start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0.46 and 0.290.290.290.29, both of which are reasonably compatible to those found in Tostevin and Brown (2006); Tostevin et al. (2013), even when considering the uncertainties in the separation energy. The reduction factors for two-nucleon knockout Simpson and Tostevin (2010); Tostevin and Brown (2006); Tostevin et al. (2013) (also including two-neutron knockout from proton-rich nuclei) and (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) reactions are presented in Figure 2 as a function of the difference between the separation energy of two protons and two neutrons ΔS=S2pS2nΔ𝑆subscript𝑆2𝑝subscript𝑆2𝑛\Delta S=S_{2p}-S_{2n}roman_Δ italic_S = italic_S start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT - italic_S start_POSTSUBSCRIPT 2 italic_n end_POSTSUBSCRIPT for two-proton knockout and (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) and S2nS2psubscript𝑆2𝑛subscript𝑆2𝑝S_{2n}-S_{2p}italic_S start_POSTSUBSCRIPT 2 italic_n end_POSTSUBSCRIPT - italic_S start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT for two-neutron knockout, as a parallel to the renowned figure for single-nucleon removal Gade et al. (2008); Tostevin and Gade (2021a, b). Rssubscript𝑅𝑠R_{s}italic_R start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT for two-proton(neutron) knockout correspond to the blue circles(squares) while the red diamonds correspond to the (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) reactions studied in this work. As can be seen in the figure the tendency is similar for both reactions, showing factors 1similar-toabsent1\sim 1∼ 1 for ΔS0similar-toΔ𝑆0\Delta S\sim 0roman_Δ italic_S ∼ 0 and a significant, roughly constant reduction for more asymmetric nuclei. It is quite remarkable that both reactions show similar trends in the description of the cross sections, while for single-nucleon removal reactions a different trend was found for the reduction factors on ΔS=Sn(p)Sp(n)Δ𝑆subscript𝑆𝑛𝑝subscript𝑆𝑝𝑛\Delta S=S_{n(p)}-S_{p(n)}roman_Δ italic_S = italic_S start_POSTSUBSCRIPT italic_n ( italic_p ) end_POSTSUBSCRIPT - italic_S start_POSTSUBSCRIPT italic_p ( italic_n ) end_POSTSUBSCRIPT for one-nucleon knockout and (p,2p)𝑝2𝑝(p,2p)( italic_p , 2 italic_p ) reactions Aumann et al. (2021), which has been alleged to originate from deficiencies in the reaction mechanism Gómez-Ramos et al. (2023) or in the description of the wavefunctions of the removed nucleons Bertulani et al. (2021); Li et al. (2022); Bertulani (2023). The fact that the trends reconcile for two-nucleon removal reactions could be related to the stronger peripherality of these reactions Simpson et al. (2009), when compared to one-nucleon removal ones, so that the nuclear interior, where the reaction mechanism and the wavefunctions are worse understood, plays a smaller role in the reaction. Further study on these trends is required to clarify these issues.

Refer to caption
Figure 2: Reduction factors for 2-proton (blue circles) and 2-neutron (blue squares) knockout reactions Simpson and Tostevin (2010); Tostevin and Brown (2006); Tostevin et al. (2013) and (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) reactions, as a function of S2pS2nsubscript𝑆2𝑝subscript𝑆2𝑛S_{2p}-S_{2n}italic_S start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT - italic_S start_POSTSUBSCRIPT 2 italic_n end_POSTSUBSCRIPT for two-proton removal and S2sS2psubscript𝑆2𝑠subscript𝑆2𝑝S_{2s}-S_{2p}italic_S start_POSTSUBSCRIPT 2 italic_s end_POSTSUBSCRIPT - italic_S start_POSTSUBSCRIPT 2 italic_p end_POSTSUBSCRIPT for two-neutron removal. It can be seen that the two types of reactions follow similar trends.

IV Summary and conclusions

In this work, we present the first, to our knowledge, theoretical description of the sequential (p,3p)𝑝3𝑝(p,3p)( italic_p , 3 italic_p ) 2-proton removal direct reaction able to produce cross sections for specific states of the final nucleus based on the two-nucleon amplitudes from nuclear structure calculations. This theory relies on the eikonal description of the interaction between protons and nucleus and the assumption of quasi-free collisions between the incoming proton and the removed protons. We have validated this theory through comparison with experimental data for the C12(p,3p)10BesuperscriptC12superscript𝑝3𝑝10Be{}^{12}\mathrm{C}(p,3p)^{10}\mathrm{Be}start_FLOATSUPERSCRIPT 12 end_FLOATSUPERSCRIPT roman_C ( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 10 end_POSTSUPERSCRIPT roman_Be reaction, finding reasonable agreement. We have also studied the Mg28(p,3p)26NesuperscriptMg28superscript𝑝3𝑝26Ne{}^{28}\mathrm{Mg}(p,3p)^{26}\mathrm{Ne}start_FLOATSUPERSCRIPT 28 end_FLOATSUPERSCRIPT roman_Mg ( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 26 end_POSTSUPERSCRIPT roman_Ne reaction finding parallels with the results found in 2-proton knockout reactions with heavier targets and applied our theory to the Ca54(p,3p)52ArsuperscriptCa54superscript𝑝3𝑝52Ar{}^{54}\mathrm{Ca}(p,3p)^{52}\mathrm{Ar}start_FLOATSUPERSCRIPT 54 end_FLOATSUPERSCRIPT roman_Ca ( italic_p , 3 italic_p ) start_POSTSUPERSCRIPT 52 end_POSTSUPERSCRIPT roman_Ar reaction, finding an overestimation of the cross section of a factor similar-to\sim2-3, which is consistent with results found in 2-proton knockout reactions. This contrasts with the discrepancy found in single-nucleon reactions. Further study on the two-nucleon-removal trends may prove illuminating for the still-unknown causes of the one-nucleon-removal reactions discrepancy. Application of this formalism to other targets, such as those measured in Frotscher et al. (2020) or Taniuchi et al. (2019), requires the obtainment of two-nucleon amplitudes, which can prove challenging for these heavy nuclei. Further developments of this work may include the use of microscopic optical potentials for the exotic cases where global parametrizations may prove less reliable, as well as the development of a theory to obtain the momentum distribution of the residual target and the angular correlation between the outgoing protons.

Acknowledgements.
The author thanks Y. Utsuno for providing the two-nucleon amplitudes used for the 54Ca case and H. Liu for the access to the preliminary experimental data for 54Ca. The author would like to thank as well J. Gómez-Camacho, A.M. Moro and A. Obertelli for illuminating discussions and critical reading of the manuscript. The author acknowledges financial support by MCIN/AEI/10.13039/501100011033 under I+D+i project No. PID2020-114687GB-I00 and under grant IJC2020-043878-I (also funded by “European Union NextGenerationEU/PRTR”), by the Consejería de Economía, Conocimiento, Empresas y Universidad, Junta de Andalucía (Spain) and “ERDF-A Way of Making Europe” under PAIDI 2020 project No. P20_01247, by the European Social Fund and Junta de Andalucía (PAIDI 2020) under grant number DOC-01006 and by the Alexander von Humboldt foundation.

References

  • Bazin et al. (2003) D. Bazin, B. A. Brown, C. M. Campbell, J. A. Church, D. C. Dinca, J. Enders, A. Gade, T. Glasmacher, P. G. Hansen, W. F. Mueller, H. Olliver, B. C. Perry, B. M. Sherrill, J. R. Terry, and J. A. Tostevin, Phys. Rev. Lett. 91, 012501 (2003).
  • Fridmann et al. (2005) J. Fridmann, I. Wiedenhöver, A. Gade, L. T. Baby, D. Bazin, B. A. Brown, C. M. Campbell, J. M. Cook, P. D. Cottle, E. Diffenderfer, D.-C. Dinca, T. Glasmacher, P. G. Hansen, K. W. Kemper, J. L. Lecouey, W. F. Mueller, H. Olliver, E. Rodriguez-Vieitez, J. R. Terry, J. A. Tostevin, and K. Yoneda, Nature 435, 922 (2005).
  • Yoneda et al. (2006) K. Yoneda, A. Obertelli, A. Gade, D. Bazin, B. A. Brown, C. M. Campbell, J. M. Cook, P. D. Cottle, A. D. Davies, D.-C. Dinca, T. Glasmacher, P. G. Hansen, T. Hoagland, K. W. Kemper, J.-L. Lecouey, W. F. Mueller, R. R. Reynolds, B. T. Roeder, J. R. Terry, J. A. Tostevin, and H. Zwahlen, Phys. Rev. C 74, 021303 (2006).
  • Tostevin et al. (2004a) J. A. Tostevin, G. Podolyák, B. A. Brown, and P. G. Hansen, Phys. Rev. C 70, 064602 (2004a).
  • Tostevin et al. (2004b) J. Tostevin, P. Batham, G. Podolyák, and I. Thompson, Nuclear Physics A 746, 166 (2004b), proceedings of the Sixth International Conference on Radioactive Nuclear Beams (RNB6).
  • Fridmann et al. (2006) J. Fridmann, I. Wiedenhöver, A. Gade, L. T. Baby, D. Bazin, B. A. Brown, C. M. Campbell, J. M. Cook, P. D. Cottle, E. Diffenderfer, D.-C. Dinca, T. Glasmacher, P. G. Hansen, K. W. Kemper, J. L. Lecouey, W. F. Mueller, E. Rodriguez-Vieitez, J. R. Terry, J. A. Tostevin, K. Yoneda, and H. Zwahlen, Phys. Rev. C 74, 034313 (2006).
  • Gade et al. (2006) A. Gade, R. V. F. Janssens, D. Bazin, R. Broda, B. A. Brown, C. M. Campbell, M. P. Carpenter, J. M. Cook, A. N. Deacon, D.-C. Dinca, B. Fornal, S. J. Freeman, T. Glasmacher, P. G. Hansen, B. P. Kay, P. F. Mantica, W. F. Mueller, J. R. Terry, J. A. Tostevin, and S. Zhu, Phys. Rev. C 74, 021302 (2006).
  • Gade et al. (2007) A. Gade, P. Adrich, D. Bazin, M. D. Bowen, B. A. Brown, C. M. Campbell, J. M. Cook, S. Ettenauer, T. Glasmacher, K. W. Kemper, S. McDaniel, A. Obertelli, T. Otsuka, A. Ratkiewicz, K. Siwek, J. R. Terry, J. A. Tostevin, Y. Utsuno, and D. Weisshaar, Phys. Rev. Lett. 99, 072502 (2007).
  • Bastin et al. (2007) B. Bastin, S. Grévy, D. Sohler, O. Sorlin, Z. Dombrádi, N. L. Achouri, J. C. Angélique, F. Azaiez, D. Baiborodin, R. Borcea, C. Bourgeois, A. Buta, A. Bürger, R. Chapman, J. C. Dalouzy, Z. Dlouhy, A. Drouard, Z. Elekes, S. Franchoo, S. Iacob, B. Laurent, M. Lazar, X. Liang, E. Liénard, J. Mrazek, L. Nalpas, F. Negoita, N. A. Orr, Y. Penionzhkevich, Z. Podolyák, F. Pougheon, P. Roussel-Chomaz, M. G. Saint-Laurent, M. Stanoiu, I. Stefan, F. Nowacki, and A. Poves, Phys. Rev. Lett. 99, 022503 (2007).
  • Adrich et al. (2008) P. Adrich, A. M. Amthor, D. Bazin, M. D. Bowen, B. A. Brown, C. M. Campbell, J. M. Cook, A. Gade, D. Galaviz, T. Glasmacher, S. McDaniel, D. Miller, A. Obertelli, Y. Shimbara, K. P. Siwek, J. A. Tostevin, and D. Weisshaar, Phys. Rev. C 77, 054306 (2008).
  • Simpson et al. (2009) E. C. Simpson, J. A. Tostevin, D. Bazin, and A. Gade, Phys. Rev. C 79, 064621 (2009).
  • Simpson and Tostevin (2009) E. C. Simpson and J. A. Tostevin, Phys. Rev. C 79, 024616 (2009).
  • Wimmer et al. (2012) K. Wimmer, D. Bazin, A. Gade, J. A. Tostevin, T. Baugher, Z. Chajecki, D. Coupland, M. A. Famiano, T. K. Ghosh, G. F. Grinyer, R. Hodges, M. E. Howard, M. Kilburn, W. G. Lynch, B. Manning, K. Meierbachtol, P. Quarterman, A. Ratkiewicz, A. Sanetullaev, E. C. Simpson, S. R. Stroberg, M. B. Tsang, D. Weisshaar, J. Winkelbauer, R. Winkler, and M. Youngs, Phys. Rev. Lett. 109, 202505 (2012).
  • Beceiro-Novo et al. (2015) S. Beceiro-Novo, T. Ahn, D. Bazin, and W. Mittig, Progress in Particle and Nuclear Physics 84, 124 (2015).
  • Obertelli et al. (2014) A. Obertelli, A. Delbart, S. Anvar, L. Audirac, G. Authelet, H. Baba, B. Bruyneel, D. Calvet, F. Château, A. Corsi, P. Doornenbal, J.-M. Gheller, A. Giganon, C. Lahonde-Hamdoun, D. Leboeuf, D. Loiseau, A. Mohamed, J. P. Mols, H. Otsu, C. Péron, A. Peyaud, E. C. Pollacco, G. Prono, J.-Y. Rousse, C. Santamaria, and T. Uesaka, The European Physical Journal A 50, 8 (2014).
  • Hansen and Tostevin (2003) P. Hansen and J. Tostevin, Annu. Rev. of Nucl. and Part. Sci 53, 219 (2003).
  • Jacob and Maris (1966) G. Jacob and T. A. J. Maris, Rev. Mod. Phys. 38, 121 (1966).
  • Aumann et al. (2013) T. Aumann, C. A. Bertulani, and J. Ryckebusch, Phys. Rev. C 88, 064610 (2013).
  • Moro (2015) A. M. Moro, Phys. Rev. C 92, 044605 (2015).
  • Ogata et al. (2015) K. Ogata, K. Yoshida, and K. Minomo, Phys. Rev. C 92, 034616 (2015).
  • Taniuchi et al. (2019) R. Taniuchi, C. Santamaria, P. Doornenbal, A. Obertelli, K. Yoneda, G. Authelet, H. Baba, D. Calvet, F. Château, A. Corsi, A. Delbart, J.-M. Gheller, A. Gillibert, J. D. Holt, T. Isobe, V. Lapoux, M. Matsushita, J. Menéndez, S. Momiyama, T. Motobayashi, M. Niikura, F. Nowacki, K. Ogata, H. Otsu, T. Otsuka, C. Péron, S. Péru, A. Peyaud, E. C. Pollacco, A. Poves, J.-Y. Roussé, H. Sakurai, A. Schwenk, Y. Shiga, J. Simonis, S. R. Stroberg, S. Takeuchi, Y. Tsunoda, T. Uesaka, H. Wang, F. Browne, L. X. Chung, Z. Dombradi, S. Franchoo, F. Giacoppo, A. Gottardo, K. Hadyńska-Klek, Z. Korkulu, S. Koyama, Y. Kubota, J. Lee, M. Lettmann, C. Louchart, R. Lozeva, K. Matsui, T. Miyazaki, S. Nishimura, L. Olivier, S. Ota, Z. Patel, E. Şahin, C. Shand, P.-A. Söderström, I. Stefan, D. Steppenbeck, T. Sumikama, D. Suzuki, Z. Vajta, V. Werner, J. Wu, and Z. Y. Xu, Nature 569, 53 (2019).
  • Frotscher et al. (2020) A. Frotscher, M. Gómez-Ramos, A. Obertelli, P. Doornenbal, G. Authelet, H. Baba, D. Calvet, F. Château, S. Chen, A. Corsi, A. Delbart, J.-M. Gheller, A. Giganon, A. Gillibert, T. Isobe, V. Lapoux, M. Matsushita, S. Momiyama, T. Motobayashi, M. Niikura, H. Otsu, N. Paul, C. Péron, A. Peyaud, E. C. Pollacco, J.-Y. Roussé, H. Sakurai, C. Santamaria, M. Sasano, Y. Shiga, N. Shimizu, D. Steppenbeck, S. Takeuchi, R. Taniuchi, T. Uesaka, H. Wang, K. Yoneda, T. Ando, T. Arici, A. Blazhev, F. Browne, A. M. Bruce, R. Carroll, L. X. Chung, M. L. Cortés, M. Dewald, B. Ding, Z. Dombradi, F. Flavigny, S. Franchoo, F. Giacoppo, M. Górska, A. Gottardo, K. Hadyńska-Klek, Z. Korkulu, S. Koyama, Y. Kubota, A. Jungclaus, J. Lee, M. Lettmann, B. D. Linh, J. Liu, Z. Liu, C. Lizarazo, C. Louchart, R. Lozeva, K. Matsui, T. Miyazaki, K. Moschner, S. Nagamine, N. Nakatsuka, C. Nita, S. Nishimura, C. R. Nobs, L. Olivier, S. Ota, Z. Patel, Z. Podolyák, M. Rudigier, E. Sahin, T. Y. Saito, C. Shand, P.-A. Söderström, I. G. Stefan, T. Sumikama, D. Suzuki, R. Orlandi, V. Vaquero, Z. Vajta, V. Werner, K. Wimmer, J. Wu, and Z. Xu, Phys. Rev. Lett. 125, 012501 (2020).
  • Goldberger and Watson (1967) M. L. Goldberger and K. M. Watson, Collision theory (John Wiley & sons. Inc., 1967).
  • Joachaim (1975) C. J. Joachaim, Quantum collision theory (North Holland Publishing Company, Amsterdam, 1975).
  • Kawai and Weidenmüller (1992) M. Kawai and H. A. Weidenmüller, Phys. Rev. C 45, 1856 (1992).
  • Hama et al. (1990) S. Hama, B. C. Clark, E. D. Cooper, H. S. Sherif, and R. L. Mercer, Phys. Rev. C 41, 2737 (1990).
  • Cooper et al. (1993) E. D. Cooper, S. Hama, B. C. Clark, and R. L. Mercer, Phys. Rev. C 47, 297 (1993).
  • Bertulani and De Conti (2010) C. A. Bertulani and C. De Conti, Phys. Rev. C 81, 064603 (2010).
  • Alex Brown (1998) B. Alex Brown, Phys. Rev. C 58, 220 (1998).
  • Olson et al. (1983) D. L. Olson, B. L. Berman, D. E. Greiner, H. H. Heckman, P. J. Lindstrom, and H. J. Crawford, Phys. Rev. C 28, 1602 (1983).
  • Warburton and Brown (1992) E. K. Warburton and B. A. Brown, Phys. Rev. C 46, 923 (1992).
  • Etchegoyen et al. (1985) A. Etchegoyen, W. D. Rae, N. S. Godwin, W. A. Richter, C. H. Zimmerman, B. A. Brown, W. E. Ormand, and J. S. Winfield, MSU-NSCL Report  (1985).
  • Cohen and Kurath (1965) S. Cohen and D. Kurath, Nuclear Physics 73, 1 (1965).
  • Panin (2012) V. Panin, PhD. Thesis, Technische Universität Darmstadt  (2012).
  • Simpson and Tostevin (2011) E. C. Simpson and J. A. Tostevin, Phys. Rev. C 83, 014605 (2011).
  • Simpson and Tostevin (2010) E. C. Simpson and J. A. Tostevin, Phys. Rev. C 82, 044616 (2010).
  • (37) H. Liu, (private communication).
  • (38) Y. Utsuno, (private communication).
  • Steppenbeck et al. (2013) D. Steppenbeck, S. Takeuchi, N. Aoi, P. Doornenbal, M. Matsushita, H. Wang, H. Baba, N. Fukuda, S. Go, M. Honma, J. Lee, K. Matsui, S. Michimasa, T. Motobayashi, D. Nishimura, T. Otsuka, H. Sakurai, Y. Shiga, P.-A. Söderström, T. Sumikama, H. Suzuki, R. Taniuchi, Y. Utsuno, J. J. Valiente-Dobón, and K. Yoneda, Nature 502, 207 (2013).
  • (40) National nuclear data center, information extracted from the nudat database.
  • Tostevin and Brown (2006) J. A. Tostevin and B. A. Brown, Phys. Rev. C 74, 064604 (2006).
  • Tostevin et al. (2013) J. A. Tostevin, B. A. Brown, and E. C. Simpson, Phys. Rev. C 87, 027601 (2013).
  • Gade et al. (2008) A. Gade, P. Adrich, D. Bazin, M. D. Bowen, B. A. Brown, C. M. Campbell, J. M. Cook, T. Glasmacher, P. G. Hansen, K. Hosier, S. McDaniel, D. McGlinchery, A. Obertelli, K. Siwek, L. A. Riley, J. A. Tostevin, and D. Weisshaar, Phys. Rev. C 77, 044306 (2008).
  • Tostevin and Gade (2021a) J. A. Tostevin and A. Gade, Phys. Rev. C 103, 054610 (2021a).
  • Tostevin and Gade (2021b) J. A. Tostevin and A. Gade, Phys. Rev. C 103, 054610 (2021b).
  • Aumann et al. (2021) T. Aumann, C. Barbieri, D. Bazin, C. Bertulani, A. Bonaccorso, W. Dickhoff, A. Gade, M. Gómez-Ramos, B. Kay, A. Moro, T. Nakamura, A. Obertelli, K. Ogata, S. Paschalis, and T. Uesaka, Progress in Particle and Nuclear Physics 118, 103847 (2021).
  • Gómez-Ramos et al. (2023) M. Gómez-Ramos, J. Gómez-Camacho, and A. Moro, Physics Letters B 847, 138284 (2023).
  • Bertulani et al. (2021) C. A. Bertulani, A. Idini, and C. Barbieri, Phys. Rev. C 104, L061602 (2021).
  • Li et al. (2022) J. Li, C. A. Bertulani, and F. Xu, Phys. Rev. C 105, 024613 (2022).
  • Bertulani (2023) C. Bertulani, Physics Letters B 846, 138250 (2023).