Protonium Formation in Antiproton Hydrogen Collisions

Transcription

1 WDS'8 Proceedings of Contributed Papers, Part III, , 28. ISBN MATFYZPRESS Protonium Formation in Antiproton Hydrogen Collisions J. Eliášek Charles University in Prague, Faculty of Mathematics and Physics, Czech Republic. Abstract. We present a computation of the cross section of protonium formation in antiproton-hydrogen collisions using the nonlocal resonance model. It is based on restoration of the Born-Oppenheimer approximation, when the adiabatic electronic basis is replaced with the so called diabatic one. This allows us to separate the electron and hadron parts of our problem at the cost of a nonlocal potential in hadron dynamics. Introduction Many groups used different methods to describe collisions of a hydrogen atom and negatively charged heavy particles. The first article about the formation of an atom which has a heavy particle instead of an electron was written by Fermi and Teller. In their work (studying the capture of heavy negative particles in matter) [Fermi and Teller, 1947] they also deduced the critical distance between the proton and the heavy particle now called the Fermi-Teller radius. If the distance is smaller then this radius then the potential is not strong enough to hold the electron bound. This means that we can not use the adiabatic approximation, because all electronic bound states disappear in continuum (see figure 1). More sophisticated methods must be used. Older works are summarized in [Cohen, 24]. Advanced adiabatic methods were used by Ovchinnikov and Macek [25] and by Žďánská et al. [24]. Hesse et al. [24] employed the hyperspherical close coupling method. The wave packet evolution and semiclassical approach were used by Sakimoto [22a, 21, 22b]. Yamanaka and Ichimura [26] solved Faddeev equations for this three body system and calculated probabilities of protonium formation. Some of the approaches scale down the mass of the antiproton to be able to use their method and get some qualitative results, because the antiproton is captured in states with n 4 and l remains approximately constant in the collision. The reaction p + H p p + e is also a prototype of the associative detachment reaction studied in atomic physics [Bieniek and Dalgrano, 1979; Čížek et al., 1998]. As a one electron system, free of electron correlation, it is vital for the understanding of errors that are caused by some approximations used in such calculations. Next, we introduce the theoretical background of our calculation of hadron dynamics during the collision and show some results for a preliminary fit to complex level shift which we took from our previous calculation for electron motion with fixed hadrons. Theory The nonlocal resonance model for collisions of negative particles with an atom was developed elsewhere (see [Domcke, 1991; Čížek, 1999]) and here we summarize it. The Hamiltonian of our problem is (we use atomic units unless stated otherwise) H = T p p + H el, T p p = P2 2µ, (1) H el = p2 2 1 r p e r p e R, where p and P denote the momenta of the electron and hadrons in the center of mass system respectively. We introduce an expansion of the wave function Ψ = i ψ i (R)φ i (r; R), (2) where R is the p- p distance, r is the position of the electron and ψ i and φ i are the nucleonic and electronic parts of the wave function respectively. The electronic part depends parametrically on R. An expansion similar to this is used in the adiabatic approximation: there, φ i are eigenfunctions of the electronic part 168

2 -.2 Potential energy [au] Proton-antiproton distance [au] m= m=1 m=2 Coulombic attraction Polarization Figure 1. The solid line indicates the Coulombic interaction between the proton and the antiproton (continuum threshold). Bound states of the electron in the proton-antiproton potential with the projection of angular momentum on p- p axis m between and 2, i.e. adiabatic potential curves, are shown. Polarization has the same asymptotic behaviour as the ground state of the electron in the field of the proton and the antiproton. We can clearly see that bound states are converging to electronic continuum with the decrease of proton-antiproton distance. of the Hamiltonian H el for each fixed R. In our approach we build this basis differently. We fix one state φ d called the discrete state, in our case the ground state of the hydrogen atom, and the rest is assumed to form a continuum φ k orthogonal to φ d. By doing this we separate our Hilbert space into two parts a part corresponding to the discrete state (we will call it Q) and a part corresponding to the rest (P). Let us define projection operators to these subspaces Q = φ d φ d 1 p p, P = φ k φ k 1 p p = 1 Q. (3) This procedure ensures, with properly chosen discrete state, that φ i will change slowly with R and that the main idea of the Born-Oppenheimer approximation may be used, namely that Q commutes with the kinetic energy operator T p p. Now let us focus on the nucleon dynamics. We write the Hamiltonian in the form H = H + H I, H = PHP + QHQ, H I = PHQ + QHP, (4) so that the Lippmann-Schwinger equation reads Ψ (+) = Ψ i + [E H P H Q + iǫ] 1 [PHQ + QHP] Ψ (+), (5) where H P and H Q stand for PHP and QHQ respectively. If we multiply this equation by the projection operators P, Q and eliminate P Ψ (+), we get an inhomogeneous Schrödinger equation with the effective Hamiltonian H eff = H Q + H QP G P (E)H PQ. (6) In the Q-space we can write the effective Hamiltonian in the form H eff = T p p + V d + F(E), (7) where V d = φ d H el φ d is the energy of the system in the discrete state and R F(E) R = dedωv de (R)(E T p p V e + iǫ) 1 Vde(R ) (8) 169

3 , -Γ/2 [au] R=.5 R=1. R= R=.5 R=1. R= E[au] Figure 2. A fit of the data from [Eliášek, 27] used in tuning calculations. In the upper part there is Γ(E,R) 2 and in lower part there is (E, R) for a few values of the parameter R. The input data are displayed by dots, the fit by lines. We used three terms in the series (11). is the R-representation of the nonlocal potential with V being the interaction potential between the hadrons, i.e. V = 1/R and V de is the coupling between the discrete state and the continuum state φ k with the electron energy e = k 2 /2. The nonlocal potential is often split into Hermitian and antihermitian parts F(E, R, R ) = (E T p p V, R, R ) i 2 Γ(E T p p V, R, R ). (9) The electronic part of the Hamiltonian is fully parametrized by the functions V (R), V d (R), Γ and. To calculate the nonlocal potential (9) we used the energy shift (E, R) and the decay width Γ(E, R) from previous work [Eliášek, 27]. There we calculated F(E, R) complex level shift in fixed nuclei approach, using the formula F(E, R) = φ d H el PḠ(E)PH el φ d, (1) where Ḡ(E) is Green s function of the electron in a fixed nuclei potential calculated in P part of the full Hilbert space. This problem is separable in spheroidal coordinates and was solved in [Eliášek, 27]. For the numerical solution of nucleon dynamics we assume that and Γ can be written as series of separable terms Γ(E, R) = f i (E)gi 2 (R), (E, R) = f i (E)g i 2 (R). (11) i i The equation (9) is then rewritten as F(E, R, R ) = ν R ν g i (R)[f i (E E ν) i 2 f i(e E ν )]g i (R ) ν R (12) where E ν are the eigenvalues and ν are the eigenvectors of T p p + V (R). In our problem there is an infinite number of bound states. Our approach to this problem was to enclose the protonium in a box. This limits the number of bound states to a finite value. We can easily add more states by changing the size of the box and check the convergence. Practically, we are searching for zero points in solutions of a differential equation for Couloubic problem at the border of the box. Once we set up the nonlocal potential F, we can write the final Lippmann-Schwinger equation ψ (+) = ψ + G (+) Q (E)F(E) ψ(+), (13) 17

4 Potential energy[au] R FT R[au] V V d V loc Polarization Figure 3. The potential curves of our problem. V is the potential of the proton-antiproton interaction, V d is the energy of the discrete state and V loc is the adiabatic potential curve of the ground state within our model. where G (+) Q Model is the Green s function for Hamiltonian H Q, and calculate cross section σ = 2(2π) 5 µ 2 ν V de f ψ (+) 2 (14) We fit the data from [Eliášek, 27] for the level shift function F(E, R) by the series (11) shown in figure 2. For the testing of the numerical parameters of the calculation we used three terms in this series (lines in figure 2). The accuracy is about 1% for small R and decreases for larger R. The most important region for nucleon dynamics is around the Fermi-Teller radius R FT =.639 au, where the lowest bound state vanishes. As stated previously, V d, V and F contain full information about the electronic problem. Specially, it is possible to find the energy V loc (R) of the bound state for R > R FT, i.e. the adiabatic potential. As shown by Domcke [1991], V loc = E res + V, where E res is the solution of the equation E res (R) = V d (R) V + (E res, R, R). (15) This potential curve is shown in figure 3. The function should be the same as the lowest potential curve in figure 1. We see that it does not have the same asymptotic behaviour (as polarization) and it joins the continuum at R greater than the Fermi-Teller radius, which is due to inaccuracy of our fit for F(E, R). Results and future prospects We calculated the integral cross section of the formation of protonium (fig.4) from this model using equations (13), (14). Values of the total cross section are of the same order as results of calculations of Tong et al. [27]. We can clearly see the threshold at E D =.5 au that corresponds to the opening of an ionization channel. The structures at lower energies correspond to the possibility of creation of protonium with different principal quantum numbers (of this hydrogen-like atom), i.e. they are excitation thresholds. Below the lowest excitation threshold (roughly.5 ev) the cross section exhibits a minor deviation from monotonic behaviour. We believe that this has no physical significance and is probably caused by an inaccurate fit for F(E, R). Now we are working on a better fit including more terms in (11) to get better agreement for Γ and and to repair correct polarization asymptote and the point of disappearance of the bound state. Next, we would like to include final electronic states with non-zero angular momentum (now we include only an outgoing electronic s-wave). This leads to a system of coupled Lippmann-Schwinger equations. And finally we would like to use more discrete states in this problem. This will introduce a multicomponent wave function and enable the possibility of studying antiproton collisions with an excited hydrogen atom. Acknowledgments. This work is supported by GAUK through Grant No 116/ I would like to thank M. Čížek for substantial support during the preparation of this paper. 171

5 8 7 6 cross section [Å 2 ] E D E[eV] Figure 4. Dependence of the total cross section of the formation of protonium in antiproton-hydrogen collision on H- p kinetic energy. Dissociative threshold is indicated by the arrow at 13.6 ev. Thresholds corresponding to the possibility of protonium creation in different internal states (different principal quantum number and therefore energy - protonium is hydrogen-like atom) can be recognized at lower energies. For reason of lucidity, we don t depict cross section value greater than 8 Å 2. The energy dependence in this region exhibits almost monotonic behaviour with limiting a value of 16 Å 2, although small oscillations can be observed below the lowest excitation threshold (roughly.5 ev). We believe that these structures are of no physical importance being probably caused by an inaccurate fit for complex level shift. References Bieniek, R. J. and A. Dalgrano, Associative detachment in collicons of H and H, Astrophys. J., 228, Cohen, J. S., Capture of negative exotic particles by atom, ions and molecules, Rep. Prog. Phys. 67, , 24. Čížek, M., Resonant Processes in atomic collisions - theoretical considerations and calculations, PhD Thesis, Charles University, Prague Čížek, M., J. Horáček and W. Domcke, Nuclear dynamics of H 2 collicion complex beyond the local approximation: associative detachment and dissociative attachment to rotationally and vibrationally excited molecules, J. Phys. B: At. Mol. Opt. Phys. 31, , Domcke, W., Theory of resonance and threshold effects in electron-molecule collisions: The projection-operator approach, Phys. Rep. 28, 9788, Eliášek, J., Tvorba protonia při srážce antiprotonu s atomem vodíku, Master Thesis, 27. Fermi, E., E. Teller, The Capture of Negative Mesotrons in Matter, Phys. Rev. 72, , Hesse, M., A. T. Le, and C. D. Lin, Protonium formation in the p-h collision at low energies by a diabatic approach, Phys. Rev. A 69, 52712, 24. Ovchinnikov, S. Yu., J. H. Macek, Protonium formation in antiprotonhydrogen-atom collisions, Phys. Rev. A 71, 52717, 25. Sakimoto, K., Full quantum-mechanical study of protonium formation in slow collisions of antiprotons with hydrogen atoms, Phys. Rev. A 65, 1276, 21. Sakimoto, K., Quantum mechanical and semiclassical calculations of ionization processes in slow collisions between a negative muon and a hydrogen atom, J. Phys. B: At. Mol. Opt. Phys. 35, , 22a. Sakimoto, K., Antiproton, kaon, and muon capture by atomic hydrogen, Phys. Rev. A 66, 3256, 22b. Tong, X. M., K. Hino, and N. Toshima, State-Specified Protonium Formation in Low-Energy AntiprotonHydrogen-Atom Collisions, Phys. Rev. Lett. 97, 24322, 26. Tong, X. M., T. Shirahama, K. Hino, and N. Toshima, Time-dependent approach to three-body rearrangement collisions: Application to the capture of heavy negatively charged particles by hydrogen atoms, Phys. Rev. A 75, 52711, 27. Yamanaka, N., A. Ichimura, Channel-coupling array theory for Coulomb three-body dynamics: Antiproton capture by hydrogen atoms, Phys. Rev. A 74, 1253, 26. Žďánská, P. R., H. R. Sadeghpour and N. Moiseyev, Non-Hermitian diabatic formulation of antiproton collision: protonium formation, J. Phys. B: At. Mol. Opt. Phys. 37, L35-L41,