P-wave Positron-hydrogen Scattering, Annihilation, and Positronium Formation

In a previous paper (Bhatia A.K. 2016), a hybrid theory for the scattering of positrons from hydrogen atoms was applied to calculate S-wave phase shifts, annihilation, and positronium formation cross sections. This approach is now being applied to calculate P-wave positron-hydrogen scattering. The present results, obtained using short-range correlation functions along with long-range correlations in the Schrödinger equation at the same time, agree very well with the results obtained in an earlier calculation by Bhatia et al. (1974), using the Feshbach projection operator formalism. In these earlier calculations, the correction due to the long-range correlations was applied to the variational results. In spite of the fact that this ad hoc correction destroyed the variational bound, the final results have been considered accurate. Annihilation cross-sections, positronium formation cross-sections, calculated in the distorted-wave approximation, are also presented.


Introduction
Various approximations have been used to study scattering of electrons and positrons by hydrogenic systems.The distortion produced in the target by the incident particle has been taken into account in the ansatz for the wave function for the scattering process by the method of polarized orbitals of Temkin (1959) [1].However, this method is not variational and therefore provides no bounds on the calculated phase shifts.The alternate method is to introduce separate short-range and long-range correlation functions and amalgamate them into a scattering problem via an optical potential in the scattering equation, thereby replacing the many-body Schrödinger equation with a single-particle Schrödinger equation.This approach was carried out in the formalism of the hybrid theory of Bhatia (2007Bhatia ( , 2008Bhatia ( , 2016) ) [2][3][4].The phase shifts obtained are rigorous lower bounds to the exact phase shifts.We now apply this approach to the P-wave scattering of positrons by hydrogen atoms.The details of this formalism have been given in previous publications of Bhatia (2007Bhatia ( , 2008) ) [2,3] and therefore, we only briefly describe the method.The wave function is given by Since positrons and electrons are distinguishable particles, the Pauli principle need not be invoked.However, positronium formation and annihilation of a positron with the bound electron can take place during the scattering.The positron coordinate is given by → r 1 ; the electron coordinate by → r 2 ; and the summation over λ is from 1 to N, the number of terms in the expansion; C λ are the unknown coefficients; and Φ λ L is the correlation function for the angular momentum L.
We use Rydberg units throughout [e 2 = 2, = 1, m e = 0.5, and a 0 = 1].In Equation ( 2), H is the Hamiltonian and ε(λ) is the expectation value of H.We have H in these units: Z is the nuclear charge in the above equation.In order to include polarization of the target, the effective target function can be written as: The function u 1s→p is given by: The target function is given by: In Equation ( 4), χ β (r 1 ) is the smooth cutoff function which guarantees that χ β (r 1 )/r 2 1 → 0 for r 1 → 0 .This function is given by χ β (r 1 ) = (1 − e −β r 1 ) n (7) where n ≥ 3.This function ensures that the polarization takes place throughout the range, rather than only for r 1 > r 2 as in [1].Furthermore, the function gives us another nonlinear parameter β which is a function of k, the incident positron momentum, along with the exponent n.The cutoff function χ ST used in the S-wave calculation [4] could have been used here too.However, the cutoff function given in Equation ( 7) allows flexibility to optimize the phase shifts.The nonlinear parameter β and the exponent n can be varied.Therefore, the cutoff function is different for different k which is not true of the cutoff function χ ST .The angle θ 12 is the angle between → r 1 and → r 2 .In Equation (1), L is the angular momentum; u L is the scattering function; and the function Φ L is the correlation function which can be written in terms of the generalized "radial" functions, which depend upon the radial coordinates and the Euler angles introduced by Bhatia and Temkin (1964) [5]: In the above equation, D κ,ε (ε = +or−) are the modified spherical harmonics which depend on the Euler angles θ, ϕ, and ψ.Here, κ is equal to 1 and the radial functions are given by: In Equations ( 9) and ( 10), C lmn and D lmn are independent eigenvectors because the target is in an s-state when the incident positron is in a P-wave, giving the total angular momentum equal to 1; while the target can be in a p-state when the incident positron is in an s-state, giving the total angular momentum equal to 1.Both states must be considered for the P-wave scattering.The second terms in Equations ( 9) and ( 10) can have different nonlinear parameters than the first term.However, we find that the nonlinear parameters γ, δ, and β are sufficient to obtain converged phase shifts.The radial function in the above equations is given by: R lmn (r 1 , r 2 , r 12 ; γ, δ) = e −γr 1 −δr 2 r l 1 r m 2 r n 12 (11) where l ≥ 1, m ≥ 0, and n ≥ 0. The sum l + m + n is equal to ω = 1, 2, 3, 4, 5, and 6 as given in Pekeris (1958) [6].These values correspond to 4, 10, 20, 35, 56, and 84 correlation terms in each symmetry.The variation with respect to C λ , as indicated by Bhatia (2007) [2], in the functional gives: In the above equation E = k 2 − Z 2 , we get the equation for the scattering function u(r): We give the various quantities: The direct potential is given by: and The quantities in the above equation are The polarizabilty of the target is obtained from the following expression as r → ∞ : where α(r) = 4.5 − e −2Zr ((Zr) 4 + 5(Zr) 3 + 9(Zr) 2 + 9(Zr) + 4.5) The optical potential is given by: Here, Ψ 0 is the wave function given in Equation ( 1) without the correlation term.The phase shifts are inferred from the scattering function u L ≡ u for r tending to infinity: In Table 1, we show the convergence of phase shifts η (radians) for k = 0.4 with the number of terms in the correlation function given in Equation (1).The total number of terms is 2N because there are two independent eigenvectors C lmn and D lmn .The nonlinear parameters γ, δ, and β are also given, while keeping the exponent n fixed at 3. By N = 56, the phase shift has converged to three significant figures.However, calculations have been carried out to N = 84 for all values of k, the incident momentum, in order to have converged results and to use the resulting wave functions for the calculation of Z e f f .In Table 2, results, which have not been extrapolated, are given for k = 0.1 to 0.7, the exponent n = 3 for all values of k.For k = 0.3, phase shift for 2N = 168 is not given because there is a loss of accuracy in the numerical calculation: Instead of increasing with the number of terms, the phase shift decreases to 0.065236.It is seen that the present results are higher than those of Bhatia et al. (1974) [7], obtained without extrapolation.As indicated above, the present results include short-range and long-range correlations at the same time and have variational lower bounds to the exact phase shifts.Over the years, reviews on interactions of positrons with atoms and ions have been published.Among them are reviews by Drachman (1971) [8], Ghosh et al. (1982) [9], Bhatia (2014) [10], and Kadyrov and Bray (2016) [11].

Z eff
In addition to the scattering, there is a possibility of annihilation of an incoming positron and an atomic electron with the emission of two gamma rays as given by Ferrell (1956) [13]: where α = e 2 / c is the fine-structure constant and a 0 is the Bohr radius.Positronium annihilation has been observed in solar flares by Crannell et al. (1976) [14], from the center of our galaxy by Leventhal et al. (1978) [15], and from gamma ray bursts by Cline (1978) [16].The quantity Z eff , which measures the overlap of the target electron with the positron, as indicated by DiRienzi and Drachman (2003) [17], approaches Z, the number of electrons for a free positron.For hydrogen: The normalization of u L (r 1 ) in the above Ψ( r 1 , r 2 ) for r 1 → ∞ is a plane-wave normalization: For angular momentum L, this is equivalent to In Table 3, we present values for Z eff for L = 1 as a function of k below the positronium formation threshold.The present results are compared with the previous ones of Bhatia et al (1977) [18] and also with those obtained by Humberston (1977) [19].The agreement is quite good.

Positronium Formation
Formation of positronium, the bound state of an electron and a positron, takes place when the incident positron captures an electron of the hydrogen atom before the positron and electron annihilate each other: e + + H(1s) → Ps + P (34 where Ps is the positronium atom and P is the proton.Using the distorted-wave approximation, Khan and Ghosh (1983) [20] have carried out calculations using the method of polarized orbitals of Temkin (1959) [1].The present calculations are similar to those of Khan and Ghosh.However, the scattering wave function is calculated variationally, Bhatia (2007) [2].The details of the calculation for S-wave have been given in Bhatia (2016) [4] which can be easily generalized to higher partial waves, and are not repeated here.Similarly, various experiments carried out on the formation of positronium are given in Bhatia (2016) [4] and are not repeated here.It should be pointed out that the continuum function has a plane-wave normalization as in Equation (33).We have used µ i µ f = 1/2 in Rydberg units.Cross-sections obtained without and with polarization are given in Table 4 and compared with those obtained by Khan and Ghosh (1983) [20].The present cross-sections have a maximum at k 2 = 0.75.The cross-sections increase when the polarization of the target is taken into account.
It is seen that the present results are lower than those of Khan and Ghosh (1983) [20].It could be due to the fact that they have incorrect sign [negative sign instead of positive sign] in Equation ( 4), as explained in Bhatia (2016) [4].The present calculations and those of Khan and Ghosh do not include correlation terms.The integrals occurring in the positronium formation calculation have been carried out using the Fourier transforms, as indicated in [4].The correlation terms have high powers of r 2 and r 12 and it becomes very difficult to calculate all the Fourier transforms that are required.

Conclusions
Using the hybrid theory of scattering of Bhatia (2007) [2], we have calculated phase shifts which are lower bounds to the exact phase shifts.This calculation includes the contribution of the long-range interaction −1/r 4 variationally and our results agree with the previous calculations of Bhatia et al. (1974) [7] and those of Armstead (1968) [12].The scattering functions have been used to calculate Z e f f and positronium formation cross-sections in the distorted wave approximation.The hybrid theory has been used previously to calculate accurate phase shifts, Bhatia (2007) [2]; resonance parameters, Bhatia (2008) [3]; and photoabsorption cross sections, Bhatia (2013) [21].It is expected that the present results are accurate to at least four decimal places.

Table 1 .
Convergence of phase shifts η (radians) for e + -H scattering with respect to the number of terms N for k = 0.4.

Table 2 .
Phase shifts η (radians) for e + -H scattering for 2N = 168 for various k and comparison with previous calculations.

Table 3 .
(1)ff(1)as a function of the incident positron momentum k, obtained in this calculation and compared with previous calculations.The results from this publication are for 2N = 168, except for k = 0.3.

Table 4 .
[20]arison of the present results for positronium formation cross-section (π a 2 0 ) in the partial wave L = 1 with those obtained in other calculations.Present results, without polarization; B: present results with polarization; C:Khan and Ghosh (1983)[20]results without polarization; and D: their results with polarization. A: