Nuclear effects in protonium formation low-energy three-body reaction

A three-charge-particle system (p̄, μ−, p) with an additional matter-antimatter, i.e. p̄−p+, nuclear interaction is the subject of this work. Specifically, we carry out a few-body computation of the following protonium formation reaction: p̄ + (p+μ−)1s → (p̄p)1s + μ−, where p is a proton, p̄ is an antiproton, μ− is a muon, and a bound state of p and its counterpart p̄ is a protonium atom: Pn = (p̄p). The low-energy cross sections and rates of the Pn formation reaction are computed in the framework of a Faddeevlike equation formalism. The strong p̄−p+ interaction is approximately included in this calculation.


Introduction
Obtaining and storing of low-energy antiprotons (p) is of significant scientific and practical interest and importance in current research in atomic and nuclear physics [1][2][3][4].For example, with the use of slow p's it would be possible to make low temperature ground state antihydrogen atoms H1s -a bound state of p and e + , i.e. a positron.The two-particle atom represents the simplest and stable anti-matter species.By comparing the properties of the hydrogen atom H and H it would be possible to test the fundamentals of physics, such as CPT theorem [5].In this connection it is important to mention the metastable antiprotonic helium atoms too, i.e. atomcules, such as p3 He + and p4 He + [6].In the field of the p physics these Coulomb three-body systems play a very important role.For example, with the use of high-precision laser spectroscopy atomclues allow us to measure p's charge-to-mass ratio and other fundamental constants in the standard model [7].
Together with the atomcules and the H atoms there is now a significant interest in protonium (Pn) atom too: a bound state of p and p + [8][9][10].The two-particle system is also named as antiprotonic hydrogen.In the atomic scale it is a heavy and a very small system with strong Coulomb and nuclear interactions.An interplay between these interactions occurs inside the atom.The last circumstance is responsible for interesting resonance and quasi-bound states in Pn [11].Therefore, Pn represents a very useful tool to study, for example, the anti-nucleon−nucleon ( NN) interaction potential [12,13] and annihilation processes [14,15].Additionally, as we already mentioned, the interplay between Coulomb and nuclear forces plays a significant role in the p and p quantum dynamics [16].The p+p + a e-mail: rasultanov@stcloudstate.edu b e-mail: dcguster@stcloudstate.eduelastic scattering has also been considered in many papers, see for example review [14].Additionally, Pn formation is related to charmonium -a hydrogen-like atom (cc), i.e. a bound state of a c-antiquark (c) and c-quark [17].Because of the fundamental importance of protonium and problems related to its formation, as far as bound and quasi-bound states, resonances and spectroscopy, the two-particle atom have attracted much attention last decades.
There are a number of few-body collisions to make Pn atoms at low temperatures.For example, the following three-charge-particle reaction is one of them: First of all, this collision represents a Coulomb three-body system and has been considered in many theoretical works in which different methods and computational techniques have been applied [18][19][20].We also would like to point out, that because in the process (1) a heavy particle, i.e. proton, is transferred from one negative "center", e − , to another, p, it would be difficult to apply a computational method based on an adiabatic (Born-Oppenheimer) approach [21].Additionally, it would be useful to mention here, that experimentalists use another few-body reaction to produce Pn atoms, i.e. a collision between a slow p and a positively charged molecular hydrogen ion, i.e.H + 2 : In the current work, however, we consider another three-body system of the Pn formation reaction.Specifically, we compute the cross-sections and rates of a very low energy collision between p and a muonic hydrogen atom H μ , i.e. a bound state of p + and a negative muon: Here, α=1s, 2s or 2p is Pn's final quantum atomic state.Because of the μ − participation in the reaction (3), at very low energy collisions Pn can be formed in an extremely small size (compact in the atomic scale) ground and close to ground states α.In these states the hadronic nuclear forces between p and p are much stronger than in the reaction (1) and, probably, would be very effective in order to study them.The size of the Pn atom in its ground state is only a 0 (Pn) = 2 /(e 2 0 m p /2) ∼ 50 fm!At such small distances the Coulomb interaction between p and p + becomes extremely strong.The corresponding binding energy in the Pn atom without the inclusion of the nuclear p−p + interaction is E n (Pn) = −e 4 0 m p /2/(2 n 2 ) ∼ − 10 keV.Here, we took n = 1, is the Planck constant, e 0 is the electron charge, and m p is the proton mass.In connection with this we would like to make a comment.The real p−p + binding energy, i.e. with inclusion of the strong interaction, can have a large value.This value may be comparable or even larger than m p .Therefore, it might be necessary to use a relativistic treatment to the reaction (3) in the output channel [22].The situation with very strong Coulomb interaction inside Pn can also be a reason of vacuum polarization forces.The Casimir forces can contribute (it might be quite significant) to the final cross sections and rates of the Pn formation reaction (3).It would be interesting to take into account all these physical effects in the framework of the reaction (3) and compute their influence on its rate.Hopefully, soon experimentalists will be able to carry out high quality measurements of the reaction (3).Thereafter one could compare the new results with corresponding theoretical data and fit (adjust) the p−p + strong interaction in the framework of the theoretical calculations in order to reproduce the laboratory data.This would allow us to better know and understand the strong p−p + interaction and the annihilation processes.
In the first order approximation reaction (3) can be considered as a three-charged-particle system (123) with arbitrary and comparable masses: m 1 , m 2 , and m 3 .It is shown in Fig. 1.The strong p−p + interaction can be included approximately after a solution of the Coulomb three-body problem.A few-body method based on a Faddeev-like equation formalism is applied.In this approach the three-body wave function is decomposed in two independent Faddeev-type components [23,24].Each component is determined by its own independent Jacobi coordinates.Since, the reaction (3) is considered at low energies, i.e. well below the three-body break-up threshold, the Faddeev-type components are quadratically integrable over the internal target variables r 23 and r 13 .They are shown in Fig. 1.
The next section represents the notation pertinent to the three-charged-particle system (123), the basic few-body equations, boundary conditions, detailed derivation of the set of coupled onedimensional integral-differential equations suitable for numerical calculations and the numerical computational approach used in the current paper.The muonic units, i.e. e = = m μ = 1, are used in this work, where m μ = 206.769m e is the mass of the muon.The proton (anti-proton ) mass is m p = m p=1836.152m e .
2 A three-charge-particle system: p, μ − and p + As we already mentioned in Introduction, a quantum-mechanical few-body approach is applied in this work.A coordinate space representation is used.This approach is based on a reduction of the total three-body wave function Ψ on two or three Faddeev-type components [24].When one has two negative and one positive charges, only two asymptotic configurations are possible below the total energy break-up threshold.The situation is explained in Fig. 1 in the case of the title three-body system.Therefore, one can decompose Ψ only on two components and write down a set of two coupled equations [25,26].A modified close coupling method is applied in order to solve these equations [27][28][29].This means to carry out an expansion of the Faddeev-type components into eigenfunctions of the subsystem Hamiltonians.This technique provides an infinite set of one-dimensional integral-differential equations [30,31].Within this formalism the asymptotic of the full three-body wave function contains two parts corresponding to two open channels [32].

An infinite set of coupled integral-differential few-body equations
One could use the following system of units: e = = m 3 = 1.We denote antiproton p by 1, a negative muon μ − by 2, and a proton p + by 3. Before the three-body breakup threshold two cluster asymptotic configurations are possible in the three-body system, i.e. ( 23)−1 and (13)−2.As we mentioned above, by their own Jacobi coordinates { r j3 , ρ k } as shown in Fig. 1: Here r ξ , m ξ are the coordinates and the masses of the particles ξ = 1, 2, 3 respectively.This suggests a Faddeev formulation which uses only two components.A general procedure to derive such formulations is described in work [25].In this approach the three-body wave function is represented as follows: where each Faddeev-type component is determined by its own Jacobi coordinates.Moreover, Ψ 1 ( r 23 , ρ 1 ) is quadratically integrable over the variable r 23 , and Ψ 2 ( r 13 , ρ 2 ) over the variable r 13 .To define |Ψ l , (l = 1, 2) a set of two coupled Faddeev-Hahn-type equations can be written: Here, Ĥ0 is the kinetic energy operator of the three-particle system, V i j (r i j ) are paired interaction potentials (i j = 1, 2, 3), E is the total energy.Now, let us present the equations ( 6)- (7) in terms of the adopted notation In order to separate angular variables, the wave function components Ψ i are expanded over bipolar harmonics: where ρ and r are angular coordinates of vectors ρ and r; C LM λμlm are Clebsh-Gordon coefficients; Y lm are spherical functions [33].The configuration triangle of the particles (123) is presented in Fig. 2 together with the Jacobi coordinates { r 23 , ρ 1 } and { r 13 , ρ 2 } and angles between them.The centreof-mass of the whole three-body system is designated as O.The centre-of-masses of the two-body subsystems ( 23) and ( 13) are O 1 and O 2 respectively.Substituting the following expansion: into (8), multiplying this by the appropriate biharmonic functions and integrating over the corresponding angular coordinates of the vectors r j3 and ρ k , we obtain a set of equations which for the case of the central potentials has the form: ) − l(l + 1) where the following notation has been introduced: To progress from (11) to one-dimensional equations, we apply a modified close coupling method, which consists of expanding each component of the wave function Ψ i ( r j3 , ρ k ) over the Hamiltonian eigenfunctions of subsystems: Thus, following expansions can be applied: where functions R i nl (r j3 ) are defined by the following equation: ) Substituting Eq. ( 14) into (11), multiplying by the corresponding functions R i nl (r j3 ) and integrating over r 2 j3 dr j3 yields a set of integral-differential equations for the unknown functions f i nlλ (ρ k ): where For brevity one can denote α ≡ nlλ (α ≡ n l λ ), and omit LM because all functions have to be the same.The functions f i α (ρ k ) depend on the scalar argument, but this set is still not one-dimensional, as formulas in different frames of the Jacobi coordinates: with the following mass coefficients: clearly demonstrate that the modulus of ρ j depends on two vectors, over which integration on the right-hand sides is accomplished: Therefore, to obtain one-dimensional integraldifferential equations, corresponding to equations ( 16), we will proceed with the integration over variables { ρ j , ρk }, rather than { r j3 , ρk }.The Jacobian of this transformation is γ −3 .Thus, we arrive at a set of one-dimensional integral-differential equations: where functions S ii αα (ρ j , ρ k ) are defined as follows: In the next section we show that fourfold multiple integration in equations ( 21) leads to a onedimensional integral and the expression (21) could be determined for any orbital momentum value L: where D L mm (0, ω, 0) are Wigner functions, ω is the angle between ρ j and ρ k , ν j is the angle between r k3 and ρ j , ν k is the angle between r j3 and ρ k (please see Fig. 2).Finally, we obtain an infinite set of coupled integral-differential equations for the unknown functions f 1 α (ρ 1 ) and f 2 α (ρ 2 ) [31], i.e. f i(i ) α(α ) (ρ i(i ) ) (i i = 1, 2), α and α belong to two different sets of three-body quantum numbers: The total angular momentum of the three-body system is L. Next in Eq. ( 23): where E i n is the binding energy of the subsystem (i 3).Also: and are the reduced masses.Further: D L mm (0, ω, 0) the Wigner functions, C Lm λ0lm the Clebsh-Gordon coefficients, Y lm are the spherical functions, ω is the angle between the Jacobi coordinates ρ i and ρ i , ν i is the angle between r i 3 and ρ i , ν i is the angle between r i3 and ρ i .One can show that:

Angular integrals
The details of the derivation of the angular integrals S ii αα (ρ j , ρ k ) (22) are explained below in this section.The configuration triangle, (123), is determined by the Jacobi vectors ( r j3 , ρ k ) and should be considered in an arbitrary coordinate system OXYZ.In this initial system the angle variables of the three-body Jacobi vectors { r j3 , ρ k } have the following values: r j3 = (θ j , φ j ), ρk = (Θ k , Φ k ), j k = Euler angles (Φ k , Θ k , ε) [33].Taking into account the transformation rule for the bipolar harmonics between new and old coordinate systems, one can write down the following relationships [33]: and where D L Mm (Φ k , Θ k , ε) are the Wigner functions [33].The fourfold multiple angular integration d ρ j d ρk in Eq. ( 21) can be written in the new variables and be symbolically represented as Next, taking into account the normalizing condition for the Wigner functions [33]: one can obtain the following formula: Now, let us make the next transformation of (123) in which the vertex j = 2 of (123) coincides with the centre O of the O X Y Z and O XYZ, however the axis O Z is directed along ρ j and (123) belongs to the plain O X Z .This transformation, which converts the coordinate frame O X Y Z into O X Y Z is characterized by the following Euler angles (0, ω, 0).Therefore the vectors ( r k3 , ρ j ) have the following new variables: r k3 = (ν j , π), ρ j = (0, 0).As a result of this rotation one can write down the following relationship: and obtain the following result: Now by taking into account that Y lm (0, 0) = δ m,0 √ (2l + 1)/4π [33], the bipolar harmonics in (32) are: with the use of these relationships we finally get the convenient for numerical computations Eq. ( 22).

Boundary conditions, numerics, cross sections and the reaction rates
To find a unique solution to Eqs. ( 23) appropriate boundary conditions depending on the specific physical situation need to be considered.First we impose: Next, for the three-body charge-transfer problems we apply the well known K−matrix formalism.This method has already been applied for solution of three-body problems in the framework of the Schrődinger equation [34,35] and coordinate space Faddeev equation [36].For the present rearrangement scattering problem with i+( j3) as the initial state, in the asymptotic region, it takes two solutions to Eq.( 23) to satisfy the following boundary conditions: where K ij are the appropriate coefficients, and v i (i = 1, 2) is a velocity in channel i.With the following change of variables in Eq. ( 23): (i=1, 2) we get two sets of inhomogeneous equations which are solved numerically.The coefficients K i j can be obtained from a numerical solution of the FH-type equations.The cross sections are given by the following expression: where (i, j = 1, 2) refer to the two channels and D = K 11 K 22 − K 12 K 21 .Also, from the quantummechanical unitarity principle one can derive that the scattering matrix K = K 11 K 12 K 21 K 22 has the following important feature: In this work, the relationship ( 39) is checked for all considered collision energies in the framework of the 1s, 1s+2s and 1s+2s+2p modified close coupling approximation Eqs. ( 10) and ( 14).The solution of the Eqs.( 6)-( 7) involving both components Ψ 1(2) required that we apply the expansions ( 10) and ( 14) over the angle and the distance variables respectively.However, to obtain a numerical solution for the set of coupled Eqs.(23) we only include the -s and -p waves in the expansion (10) and limit n up to 2 in the Eq. ( 14).As a result we arrive at a truncated set of six coupled integral-differential equations, since in Ψ 1(2) only 1s, 2s and 2p target two-body atomic wave-functions are included.This method represents a modified version of the close coupling approximation with six expansion functions.The set of truncated integral-differential Eqs.(23)  CNR * 15 9004 procedure, i.e. on the right side of the equations the integrals over ρ 1 and ρ 2 are replaced by sums using the trapezoidal rule [37] and the second order partial derivatives on the left side are discretized using a three-point rule [37].By this means we obtain a set of linear equations for the unknown coefficients w j S (12)  αα (ρ 1i , ρ 2 j ) f (2) w j S (21)  α α (ρ 2i , ρ 1 j ) f (1) Here, coefficients w j are weights of the integration points ρ 1i and ρ 2i (i = 1, N p ), N s is the number of quantum states which are taken into account in the expansion (14).Next, D 2 i j is the three-point numerical approximation for the second order differential operator: where Δ is a step of the grid Δ = ρ i+1 − ρ i .The vector B 21 α (i) is: B (21)  α (i) = M 2 /γ 3 N p j=1 w j S (21)  α 1s0 (i, j) sin(k 1 ρ j ), and in symbolic-operator notations the set of linear Eqs. ( 40)-(41) has the following form: The discretized equations are subsequently solved by the Gauss elimination method [38].As can be seen from Eqs. ( 40)-(41) the matrix A should have a so-called block-structure: there are four main blocks in the matrix: two of them related to the differential operators and other two to the integral operators.Each of these blocks should have sub-blocks depending on the quantum numbers α = nlλ and α = n l λ .The second order differential operators produce three-diagonal sub-matrixes [31].However, there is no need to keep the whole matrix A in computer's operating (fast) memory.The following optimization procedure shows that it would be possible to reduce the memory usage by at least four times.Indeed, the numerical equations ( 40)-(41) can be written in the following way: 12 and S 21 are sub-matrixes of A. Now one can determine that: Thereby one can obtain a reduced set of linear equations which are used to perform the calculations: To solve the coupled integral-differential equations ( 23) one needs to first compute the angular integrals Eqs.(22).They are independent of energy E. Therefore, one needs to compute them only once and then store them on a computer's hard drive (or solid state drive) to support future computation of other observables, i.e. the charge-transfer cross-sections at different collision energies.
The sub-integral expressions in ( 22) have a very strong and complicated dependence on the Jacobi coordinates ρ i and ρ i .To calculate S (ii )  αα (ρ i , ρ i ) at different values of ρ i and ρ i an adaptable algorithm has been applied together with the following mathematical substitution: The angle dependent part of the equation can be written as the following one-dimensional integral: 2L + 1 . This figure shows our numerical result for the low-energy proton transfer reaction integral cross section σ tr in the three-charge-particle collision p + H μ → (pp) α + μ − , where H μ is a muonic hydrogen atom: a bound state of a proton and a negative muon.In this work only the reaction final channel with α=1s in considered in the framework of the 1s, 1s+2s and 1s+2s+2p modified close-coupling approximation (MCCA) approach.
An adaptive algorithm which is incorporated in a FORTRAN subroutine from [39] is used in this work in order to carry out the angle integration in (43).This recursive computer program, QUADREC, is a better, modified version of the well known program QUANC8 [38].QUADREC provides a much higher quality, stable and more precise integration than does QUANC8 [39].The expression (43) differs from zero only in a narrow strip, i.e. when ρ i ≈ ρ i .This is because in the considered threebody system the coefficient β i is approximately equal to one.Therefore, in order to obtain numerically reliable converged results it is necessary to adequately distribute a very large number of discretization points (up to 6000) between 0 and ∼80 muonic units.

Results
In this section we report our computational results.The Pn formation three-body reaction is computed.A Faddeev-like equation formalism ( 6)-( 7) has been applied.The few-body approach is presented in previous sections.In order to solve the coupled equations ( 6)-( 7) two different independent sets of target expansion functions have been used (14).This method allows us to avoid the over-completeness problem and the two targets are treated equivalently.The main goal of this work is to carry out a reliable quantum-mechanical computation of the cross sections and corresponding rates of the title Pn formation reaction at low and very low collision energies.It would also be interesting to include and estimate a contribution of the strong p−p interaction.The coupled integral-differential Eqs.(23) have been solved numerically for the case of the total angular momentum L = 0 in the framework of the two-level 2×(1s), four-level 2×(1s+2s), and six-level 2×(1s+2s+2p) close coupling approximations in Eq. ( 14).The sign "2×" indicates that two different sets of expansion functions are applied.The following boundary conditions (35), (36), and (37) 38) has been applied.The constructed equations satisfy the Schrődinger equation exactly.For the energies below the three-body break-up threshold they exhibit the same advantages as the Faddeev equations [24], because they are formulated for the wave function components with correct physical asymptotes.To solve the equations, a close-coupling method is applied, which leads to an expansion of the system's wave function components into eigenfunctions of the subsystem (target) Hamiltonians providing with a set of one-dimensional integral-differential equations after the partial-wave projection.A further advantage of the Faddeev-type method is the fact that the Faddeev-components are smoother functions of the coordinates than the total wave function.
Below we report our new computational results.We compare some of our findings with the corresponding data from older work [40].The Pn formation cross section in the reaction (3) are shown in Fig. 3.Here we used only 1s, 1s + 2s and 1s + 2s + 2p states within the modified closecoupling approximation (MCCA) approach.One can see that the contribution of the 2s-and 2p-states from each target is becoming even more significant while the collision energy becomes smaller.It would be useful to make a comment about the behaviour of σ tr (ε coll ) at very low collision energies: ε coll ∼ 0. From our calculation we found that the p + transfer cross sections σ tr → ∞ as ε coll → 0. However, the p + transfer rates, λ tr , are proportional to the product σ tr × v c.m. and this trends to a finite value as v c.m. → 0.Here v c.m. = √ 2ε coll /M k is a relative center-of-mass velocity between the particles in the input channel of the three-body reactions, and M k is the reduced mass.To compute the proton transfer rate the following formula is used: Therefore, additionally, for the process (1) we compute the numerical value of the following important quantity: Λ(Pn) = λ tr = σ tr (ε coll → 0)v c.m. ≈ const, which is proportional to the actual Pn formation rate at low collision energies.functions from each target are applied.Therefore, for example, in the case of the 1s MCCA approach two expansion functions are used in our calculations.However, in the case of the 1s + 2s MCCA approach four expansion functions are applied.Additionally, as we mentioned above, the unitarity relationship, i.e.Eq. ( 39), is checked for all considered collision energies E. It is seen, that χ(E) posses a fairly constant value close to one.Figs. 4 and 5 show our convergence results.In these cases we used different number of integration points, namely 75 and 85 per one muonic radius length and different values of the up level of integration, namely 62, 69 and 76 m. a. u.Thus, the maximum amount of the integration number used in this work is N max = 76 × 85 = 6460.It is seen that the results are in a good agreement with each other for the transfer and elastic cross sections.
In the framework of the 2×(1s + 2s + 2p) MCCA approach, i.e. when six coupled Faddeev-Hahntype integral-differential equations are solved, our result for the Pn formation rate has the following value: Λ 1s2s2p (Pn) = σ tr ×v c.m. ≈ 0.32 m.a.u.The corresponding rate from the work [40] is: Λ (Pn) ≈ 0.2 m.a.u.Both results are in a fairly good agreement with each other and were multiplied by 5 as in [40].In conclusion, because of the complexity of the few-body system and the method, in this work only the total orbital momentum L = 0 has been taken into account.It was adequate in the case of the slow and ultraslow collisions discussed above.However, inclusion of the strong proton-antiproton interaction in this work by just shifting the ground state energy level of the Pn atom [13] increased our results by ∼ 50%.In a future work it would be useful to take into account higher atomic target states like 3s+3p+3d+4s+... as well as the continuum spectrum.This calculation would be an interesting investigation.However, at very low energy collisions, which are considered in this paper, all these channels are closed and located far away from the actual collision energies.At the same time the main contribution from s-and p-states (polarization) is taken into account.A direct inclusion of the strong p−p interaction [14,15] would be a very interesting future investigation in the framework of the three-body reaction (3).

Figure 1 .
Figure 1.Two asymptotic spacial configurations of the 3-body system (123), or more specifically (p, μ − , p + ) which is considered in this work.The few-body Jacobi coordinates ( ρ i , r jk ), where i j k = 1, 2, 3 are also shown together with the 3-body wave function components Ψ 1 and Ψ 2 : Ψ = Ψ 1 + Ψ 2 is the total wave function of the 3-body system.

Figure 2 .
Figure 2. The title three-charge-particle system p, μ − and p (or p + -proton) and system's configurational triangle (123) are shown together with the few-body Jacobi coordinates (vectors): { ρ 1 , r 23 } and { ρ 2 , r 13 }.Additionally, r 12 is the vector between two negative particles in the system.The needed for detailed few-body treatment geometrical angles between the vectors such as η 1(2) , ν 1(2) , ζ and ω are also presented in this figure.

2 .
Let us adopt a new coordinate system O X Y Z in which the axis O Z is directed over the vector ρ k , (123) belongs to the plain O X Z and the vertex k = 1 of (123) coincides with the origin O of the new OX Y Z .The new angle variables of the Jacobi vectors in the O X Y Z system have now the following values: r j3 = (ν k , π), ρ k = (0, 0), r k3 = (η k , π), ρ j = (ω, π), here k = 1 and j = 2.The spatial rotational transformation from OXYZ to O X Y Z has been done with the use of the following

Figure 5 .
Figure 5. Numerical convergence results for the low-energy proton transfer reaction integral cross section σ tr and elastic scattering cross section σ el multiplied by the collision velocity in p + H μ → (pp) α + μ − , where H μ is a muonic hydrogen atom and α =1s .In these calculations only 1s+2s+2p modified close-coupling approximation (MCCA) approach is used. ) is solved by a discretization have been used.To compute the charge transfer cross Numerical convergence results for the low-energy proton transfer reaction integral cross section σ tr in p + H μ → (pp) α + μ − , where H μ is a muonic hydrogen atom and α =1s .In these calculations only 1s + 2s + 2p modified close-coupling approximation (MCCA) approach is used.

Table 1 .
where H μ is a muonic hydrogen atom and α =1s .In these calculations only 1s+2s+2p modified close-coupling approximation (MCCA) approach is used.The total Pn formation cross section σ tr in the three-body reaction (3), when α = 1, and a product of this cross section and the corresponding center-of-mass velocities v c.m. between p and H μ .The results are prsented together with the corresponding unitarity ratio between the K 12 and K 21 coefficients of the scattering K-matrix.
Table1includes our data for this important parameter together with the Pn formation total cross section.All these results are obtained in the framework of the 2×(1s) and 2×(1s + 2s) approximations.The sign "2 × " means, as we mentioned, that two sets of the expansion