arxiv: v2 [physics.atom-ph] 3 Aug 2014

Transcription

1 Positron scattering and annihiation on nobe gas atoms D. G. Green, J. A. Ludow, and G. F. Gribakin Department of Appied Mathematics and Theoretica Physics, Queen s University Befast, Befast BT7 1NN, Northern Ireand, United Kingdom arxiv: v2 [physics.atom-ph] 3 Aug 214 Positron scattering and annihiation on nobe gas atoms is studied ab initio using many-body theory methods for positron energies beow the positronium formation threshod. We show that in this energy range the many-body theory yieds accurate numerica resuts and provides a nearcompete understanding of the positron-nobe-gas-atom system. It accounts for positron-atom and eectron-positron correations, incuding the poarization of the atom by the positron and the nonperturbative process of virtua positronium formation. These correations have a arge ect on the scattering dynamics and resut in a strong enhancement of the annihiation rates compared to the independent-partice mean-fied description. Computed eastic scattering cross sections are found to be in good agreement with recent experimenta resuts and Kohn variationa and convergent cose-couping cacuations. The cacuated vaues of the annihiation rate parameter Z (ective number of eectrons participating in annihiation) rise steepy aong the sequence of nobe gas atoms due to the increasing strength of the correation ects, and agree we with experimenta data. PACS numbers: 34.8.Uv, 78.7.Bj I. INTRODUCTION The scattering of ow-energy positrons from nobe gas atoms has been the subject of theoretica studies for many decades [1]. For an overview of the fied of owenergy positron scattering, see the review [2]. Athough the exchange interaction is absent, positron scattering from atoms is consideraby more chaenging to treat theoreticay than the reated probem of eectron-atom scattering. For positrons, the static interaction with the atom is repusive. At ow incident positron energies the attractive poarization potentia induced by the positron on the atom overcomes the static repusion, eading to a deicate baance between the opposing potentias. A key roe is thus payed by positron-atom and positron-eectron correations. In addition, phenomena unique to positrons occur, namey, positronium formation (virtua or rea) and positron annihiation. Positronium (Ps) formation is a process in which a positron captures an atomic eectron into a bound state. It occurs when the positron energy exceeds the Psformation threshod ε Ps = I+E 1s (Ps) = I 6.8 ev, where I is the ionization potentia of the atom and E 1s (Ps) is the ground-state energy of Ps. In positron-atom coisions this is usuay the first ineastic channe to open, and it has a pronounced ect on positron scattering [2]. Ps formation aso affects the positron-atom interaction at energies beow ε Ps, where it is virtua. Besides eastic scattering, another channe open at a positron energies is positron annihiation. For atomic and moecuar Address correspondence to dermot.green@baio.oxon.org; Present address: Joint Quantum Centre (JQC) Durham- Newcaste, Department of Chemistry, Durham University, South Road, Durham, DH1 3LE, United Kingdom. Present address: AquaQ Anaytics, Suite 5, Sturgeon Buiding, 9-15 Queen Street, Befast, BT1 6EA, UK. g.gribakin@qub.ac.uk targets the positron annihiation cross section is traditionay parameterized as σ a = πr(c/v)z 2, where r is the cassica eectron radius, c is the speed of ight, v is the incident positron veocity, and Z is the ective number of eectrons participating in the annihiation process. For Z = 1 this formua gives the basic eectronpositron annihiation cross section in the nonreativistic Born approximation [3]. For many-eectron targets Z may naivey be expected to be cose to the number of eectrons in the atom. However, the positron-atom interaction and eectron-positron correations have a strong ect on the annihiation rates [4, 5]. Experimenta studies of positron annihiation in heavier nobe gases yied Z vaues that are orders of magnitude greater than those obtained in a simpe static-fied approximation [6, 7]. In this paper we use diagrammatic many-body theory to describe the interaction of positrons with nobe gas atoms. Many-body theory aows one to understand and quantify the roe and magnitude of various correation ects. Scattering phase shifts, differentia and tota eastic scattering cross sections, and Z are cacuated ab initio with proper incusion of the correations [8]. Exceent agreement with experimenta resuts and the resuts of other sophisticated theoretica approaches is found. This work, taken together with the many-body theory cacuations of γ-spectra and rates for annihiation on core eectrons of nobe gases [9], forms a comprehensive study that provides a near-compete understanding of the positron-nobe-gas-atom system at positron energies beow the Ps-formation threshod. Many-body theory [1] provides a natura framework for the incusion of eectron-eectron and eectronpositron correations. It uses the apparatus of quantum fied theory to deveop a perturbative expansion for the ampitudes of various processes. The abiity to show various contributions pictoriay by means of Feynman diagrams makes the theory particuary transparent and heps one s intuition and understanding of many-body

2 2 quantum phenomena. This theory is ideay suited to the probem of a partice interacting with a cosed-she atom, with successfu appications to eectron scattering from nobe-gas atoms (see, e.g., Refs. [11 15]). The study of positron-atom scattering using the many-body theory thus shoud have been straightforward. However, progress in this direction was stymied by the difficuty in accounting for virtua Ps formation, as the Ps bound state cannot be accuratey described by a finite number of pertubation terms. The need for this was reaised eary on [16], but a proper soution incuding summation of an infinite series of adder diagrams was achieved ony much ater [17]. The ect of virtua Ps formation neary doubes the strength of the positron-atom correation potentia, as the terms in the series are of the same sign, eading to a arge tota. In contrast, in eectronatom scattering, such series is sign-aternating, giving a sma, often negigibe, overa contribution. The first appication of the many-body theory to positron scattering was for heium [16]. This study accounted for poarization of the target by the positron and demonstrated the importance of virtua Ps formation by using a rather crude approximation (see aso [18]). This approximation was aso used in subsequent studies for heium and other nobe-gas atoms [19, 2]. A more sophisticated approximation to the virtua Ps contribution was deveoped and appied to positron scattering, binding and annihiation in Refs. [21 24] (see aso [25]). It was ater used to cacuate rea Ps formation in nobe-gas atoms [26] and produced mixed resuts. The compete evauation of the adder-diagram series was impemented in the positron-hydrogen study [17] which used B-spine basis sets to discretize the positron and eectron continua. This approach has since been appied to positron binding to the haogen negative ions [27], to positron scattering and annihiation on hydrogen-ike ions [28], to the cacuation of gamma-ray spectra for positron annihiation on the core and vaence eectrons in atoms [9, 29] and in moecues [3]. Another many-body theory technique that aows one to sum the dominant series of diagrams to a orders is the inearized couped-custer method which was used to cacuate positron-atom bound states for a arge number of atoms [31, 32]. Recenty, a series of high-quaity experimenta measurements and convergent cose-couping (CCC) cacuations have been performed for ow-energy positron scattering aong the nobe gas sequence [33 4]. In the ight of these new data, the many-body theory approach deveoped by the authors is appied here to a thorough study of positron interaction with the nobe gas atoms. The rest of the paper is organized as foows. In Secs. II and III we describe the many-body theory and the numerica impementation of this theory. In Sec. IV we present resuts for the scattering phase shifts, differentia and tota eastic scattering cross sections, and the annihiation parameter Z (both energy-resoved and thermay averaged), and compare with existing experimenta and theoretica data. We concude with a brief summary and outook. Agebraic expressions for the many-body diagrams are provided in Appendix A and tabuated numerica resuts are in Appendix B. We use atomic units (a.u.) throughout, uness stated otherwise. II. MANY-BODY THEORY A. Dyson equation and sef-energy The many-body-theory description of a positron interacting with an atom is based on the Dyson equation (see, e.g., Ref. [41]), (H + Σ ε )ψ ε = εψ ε, (1) where ψ ε is the (quasipartice) wave function of the positron, H is the zeroth-order Hamitonian of the positron in the static fied of the atom (usuay described in the Hartree-Fock (HF) approximation), ε is the positron energy, and Σ ε is a nonoca, energy-dependent correation potentia. This potentia is equa to the sef-energy of the singe-partice Green s function of the positron in the presence of the atom [42]. It incorporates the many-body dynamics of the positron-atom interaction. As the potentia Σ ε is nonoca, Dyson s equation is an integra equation, H ψ ε (r) + Σ ε (r, r )ψ ε (r )dr = εψ ε (r). (2) The correation potentia Σ ε can be evauated as an infinite perturbation series in powers of the residua eectronpositron and eectron-eectron interactions. Because of the spherica symmetry of the atomic potentia, Eq. (2) can be soved separatey for each partia wave of the incident positron. The main contribution to the positron sef-energy Σ ε is given by the two diagrams shown in Fig. 1. The second of these in fact represents an infinite subsequence of diagrams which describe virtua Ps formation. For a positron interacting with a one-eectron target (the hydrogen atom or hydrogen-ike ion), the diagrams shown in Fig. 1 constitute a compete expansion of Σ ε. The agebraic expressions for these diagrams can be found in [17]. Diagram 1 (a), Σ (2) ε, accounts for the poarization of the atom by the positron. At arge positron-atom separations this diagram has the asymptotic behaviour, Σ (2) ε (r, r ) α de 2 2r 4 δ(r r ), (3) where α d is the static dipoe poarizabiity of the atom (here, in the HF approximation). This second-order diagram (with its exchange counterparts) is known to provide a good approximation for Σ ε in eectron-atom scattering (e.g., for argon [13] or xenon [14]). However, the same approximation in positron-atom scattering is seriousy deficient [16, 25].

3 3 ε ν ε µ n (a) n ε ν 1 ν 2 ε Γ µ 1 µ 2 n (b) FIG. 1. Main contributions to the positron sef-energy Σ ε. The owest, second-order diagram (a), Σ (2) ε, describes the ect of poarization; diagram (b), Σ (Γ) ε, accounts for virtua Ps formation represented by the Γ-bock. Top ines in the diagrams describe the incident and intermediate-state positron. Other ines with the arrows to the right are excited eectron states, and to the eft, hoes, i.e., eectron states occupied in the target ground state. Wavy ines represent Couomb interactions. Summation over a intermediate positron, eectron and hoe states is assumed. Diagram 1 (b), which we denote Σ (Γ) ε, describes the short-range attraction between the positron and the atom due to virtua Ps formation. The shaded bock Γ represents the sum of eectron-positron adder diagrams, referred to as the vertex function. It satisfies a inear integra equation represented diagrammaticay in Fig. 2 and written in the operator form as Γ = V + V χγ, (4) where Γ is the vertex function, V is the eectron-positron Couomb interaction and χ is the propagator of the intermediate eectron-positron state. With the eectron and positron continua discretized as described in Sec. III, Γ and V become matrices, with χ being a diagona matrix of energy denominators. In this case Eq. (4) is a inear matrix equation, which is easiy soved numericay [17]. Such discretization of the eectron and positron continua is vaid for energies beow the Ps formation threshod, for which the eectron-positron pair cannot escape to infinity. In order to describe the poarization of muti-eectron atoms more accuratey, a set of third-order diagrams is aso incuded in the cacuation of Σ ε. These diagrams, denoted coectivey Σ (3) ε, are shown in Fig. 3. Agrebraic expressions for these diagrams are given in Appendix A. Diagrams 3 (a), (b), (c) and (d) represent corrections to the second-order poarization diagram Fig. 1 (a) due to eectron correations of the type described by the random-phase approximation with exchange [43]. They account for the eectron-hoe interaction and screening of the positron and eectron Couomb fied. Diagram 3 (e) describes the positron-hoe repusion. Instead of computing the sef-energy Σ ε (r, r ) in the coordinate representation, it is more convenient to dea with its matrix eements ε Σ E ε = ϕ ε (r )Σ E (r, r)ϕ ε (r)drdr, (5) with respect to the zeroth-order static-fied positron wave functions ϕ ε with a given orbita anguar momentum. The atter are eigenstates of the zeroth-order Hamitonian, H ϕ ε = εϕ ε, (6) which satisfy the correct boundary conditions and are appropriatey normaized. For true continuous-spectrum positron states the radia wave functions are normaized to a δ-function of energy in Rydberg, δ(k 2 k 2 ). This corresponds to the asymptotic behaviour P () ε (r) (πk) 1/2 sin(kr π/2 + δ () ), where δ () are the static HF-fied phase shifts, and k is the wavenumber reated to the positron energy by ε = k 2 /2. The intermediate states in the diagrams are square-integrabe eectron and positron basis functions eigenstates of H constructed from B-spines in a finite-size box of radius R (see Sec. III A). B. Scattering phase shifts The sef-energy matrix (5) can be used to obtain the phase shifts directy [43, 44]. First, a reducibe sefenergy matrix ε Σ E ε is found from the integra equation, ε ε Σ E ε = ε Σ E ε ε Σ E ε Σ E ε + P E ε dε, (7) where P denotes the principa vaue of the integra. The scattering phase shift is then given by, where δ (k) = δ () (k) + δ (k), (8) tan [ δ (k)] = 2π ε Σ ε ε, (9) determines the additiona phase shift δ (k) due to positron-atom correations described by the sef-energy. The reducibe sef-energy matrix aso aows one to find the positron quasipartice wave function (i.e., soution to the Dyson equation), as ψ ε (r) = ϕ ε (r) + P ϕ ε (r) ε Σ ε ε ε ε dε. (1) Numericay, the integras in Eqs. (7) and (1) are cacuated using an equispaced positron momentum grid of 2 intervas of k =.2. In order for the quasipartice radia wave function to be correcty normaized and have the asymptotic behaviour P ε (r) 1 πk sin ( kr π 2 + δ() + δ ), (11)

4 4 ν 1 ν 2 ν 1 ν 2 ν 1 ν 2 ν 1 ν 2 ν 1 ν 2 ν 1 µ = + µ 1 µ 2 µ 1 µ µ 2 µ 1 µ 2 µ µ µ 1 µ µ 2 µ 1 Γ Γ ν ν 2 FIG. 2. Eectron-positron adder diagram series and its sum, the vertex function Γ (shaded bock). Comparison between the eft- and right-hand sides of the diagrammatic equation shows that Γ can be found by soving a inear equation, Eq. (4). ε ν ε (a) (b) (c) µ 1 µ 2 ε ν ε ε ν ε µ 1 µ 2 µ 1 µ 2 n 1 n 2 n 1 n 2 ε ν ε ε ν 1 ν 2 ε n 2 (d) n 1 µ 1 (e) µ µ 2 n 2 n 1 n 2 FIG. 3. Third-order correction diagrams, Σ (3) ε. Mirror images of the diagrams (c) and (d) are aso incuded. In a diagrams the top, horizonta ines represent the positron. the wave function obtained from Eq. (1) must be mutipied by the factor [ cos δ = 1 + ( 2π ε Σ ε ε ) ] 2 1/2. (12) C. Positron annihiation The annihiation rate λ for a positron in a gas of atoms or moecues with number density n is usuay parameterized by λ = πr 2 cnz, (13) where r is the cassica radius of the eectron, c is the speed of ight, and Z is the ective number of eectrons that participate in the annihiation process [4, 45]. In genera the parameter Z is different from the number of eectrons in the target atom, Z. In particuar, as we sha see, positron-atom and eectron-positron correations can make Z Z. Theoreticay, Z is equa to the average positron density at the ocations of the target eectrons, i.e., N Z = δ(r r i ) Ψ k (r 1,..., r N ; r) 2 dr 1... dr N dr, i=1 (14) where Ψ k (r 1,..., r N ; r) is the tota wave function which describes the scattering of the positron with momentum k by the N-eectron target. This wave function is normaized at arge positron-atom separations to the positron pane wave incident on the ground-state target with the wave function Φ : Ψ k (r 1,..., r N ; r) Φ (r 1,..., r N )e ik r. (15) Equation (14) has the form of an ampitude, with the eectron-positron deta-function acting as a perturbation. Hence, it is possibe to derive a diagrammatic expansion for Z [17, 25, 26]. Figure 4 shows the set of main annihiation diagrams. In addition to the eements found in the sef-energy diagrams, each of the Z diagrams contains one eectron-positron δ-function vertex. The diagrams in Fig. 4 provide a compete description of Z for one-eectron systems, such as hydrogen and hydrogenike ions [28]. Agebraic expressions for these diagrams can be found in [17]. The simpest, zeroth-order diagram, Fig. 4 (a), corresponds to Z () N = ϕ n (r) 2 ψ ε (r) 2 dr, (16) n=1 i.e., the overap of the eectron and positron densities (ϕ n

5 5 being the nth eectron HF ground-state orbita). It gives Z in the independent-partice approximation. For the many-eectron systems considered here, it is aso important to account for eectron screening in the cacuation of Z. A series of annihiation diagrams with two Couomb interactions, simiar to the sef-energy corrections in Fig. 3, are therefore incuded, see Fig. 5. The corresponding agebraic expressions are given in Appendix A. The externa ines in the Z diagrams represent the wave function of the incident positron. In the owest approximation, one can use the positron wave function in the static fied of the HF ground-state atom, i.e., set ψ ε = ϕ ε, negecting the ect of the correation potentia Σ ε on the positron. This ect is in fact quite strong, so to obtain accurate Z one needs to use the quasipartice positron wave function of Eq. (1). Figures 4 and 5 represent the atter case, with doube ines corresponding to the fuy correated (Dyson) positron quasipartice wave function obtained from Eq. (1), i.e., the HF positron wave functions dressed with the positron sef-energy in the fied of the atom (see Figs. 1 and 3), and normaized by the factor (12). In order to impement the correct normaization of the incident positron wave function to the pane wave e ik r at arge distances, it is necessary to mutipy the Z diagrams computed for the positron wave function with anguar momentum, by III. 4π 2 (2 + 1). (17) k NUMERICAL IMPLEMENTATION Beow we outine the numerica impementation of the many-body theory methods described in Sec. II. A. B-spine basis First, the Hartree-Fock ground state of the atom is cacuated with a standard HF package [46]. Using these numerica wave functions, direct and exchange potentias are constructed and the atomic HF Hamitonian for the positron (i.e., without exchange) or eectron (with exchange) is then diagonaized in a B-spine basis [47, 48]. The corresponding eigenvectors are used to construct the positron and eectron wave functions. This provides ectivey compete sets of positron and eectron basis states covering both bound states and the positive-energy continuum [17, 28]. These states are then used to cacuate the Couomb and δ-function matrix eements (A1) and (A3), and to evauate the many-body diagrams by summing over intermediate eectron and positron states. For the cacuations reported here sets of 4 spines of order 6 in a box of size R = 3 a.u. were used. Two outermost subshes are incuded when cacuating the sefenergy and annihiation diagrams (Figs. 1, 3, 4, and 5). The diagrams are evauated at 8 energy points from zero incident positron energy up to the Ps-formation threshod and then interpoated onto the required energies. The contributions of inner shes to Z are cacuated using the diagrams of Fig. 4 [9]. There is a point concerning boundary conditions satisfied by the B-spine basis states that affects the cacuation of the sef-energy matrix ε Σ E ε. The sefenergy matrix is evauated initiay using the B-spine basis states i as i Σ E j. The number of B-spine basis states used in each partia wave ( 15) is much smaer than the number of continuous spectrum states required for an accurate soution of Eq. (7). The change to ε Σ E ε can be made using the ective competeness of the B-spine basis on the interva [, R], ε Σ E ε = ij ε i i Σ E j j ε, (18) where ε i is the overap of the HF state ε with the B- spine basis state i. However, unike the B-spine states which satisfy the boundary condition P i (R) =, the continuous spectrum radia wave function P ε (r) is finite at the boundary r = R. To improve numerica accuracy, a weighting function f(r) = R r is inserted into Eq. (18): ε Σ E ε = ij ε f i i f 1 Σ E f 1 j j f ε, (19) with the weighted sef-energy matrix i f 1 Σ E f 1 j, being cacuated rather than i Σ E j. B. Finite box size In genera, the finite box size may affect the resuts at ow positron momenta kr < 1. In particuar, it imits the range of the poarization potentia (represented by Σ ε ) to distances not exceeding R. This is countered by adding a correction term to the sef energy. At arge distances the correation potentia is oca, energy independent and of the form α d /2r 4, where α d is the static dipoe poarizabiity. The contribution to the sef-energy matrix ε Σ E ε from distances outside the box can then be approximated by R ( P ε (r) α ) d 2r 4 P ε (r) dr, (2) with the radia wave functions given by their asymptotic form, k P ε (r) = r π [ ] j (kr) cos δ () n (kr) sin δ (), (21) where j and n are the spherica Besse and Neumann functions. The correction (2) is added to the sef-energy matrix cacuated using the many-body theory for r R.

6 6 " " " " " " " " " " µ n n n n n (a) (b) (c) (d) (c) (d) " " " " (e) n n (f) FIG. 4. Many-body-theory diagrammatic expansion for Z. The soid circe in the diagrams is the δ-function annihiation vertex, see Eq. (14). The doube ines represent the fuy correated (Dyson) positron quasipartice wave function of Eq. (1), i.e., the HF positron wave function dressed with the positron sef-energy in the fied of the atom. Diagrams (b), (d) and (e) are mutipied by two to account for their mirror images. ε ν ε ε ν 1 ν 2 ε ε µ 1 µ 2 µ n 1 n 2 n 1 n 2 ν ε ε µ 1 µ 2 n 1 n 2 (a) (b) (c) (d) ε ν ε ε ν ε ε µ 1 n 1 µ 2 µ 1 ν ε n 2 ν µ 2 n 2 ε n 2 (e) n 2 (f) (g) FIG. 5. Annihiation diagrams which account for corrections to Z due to eectron screening of the eectron and positron Couomb fied, and eectron-hoe and positron-hoe interactions. The top, horizonta ines represent the positron. The doube ines represent the fuy correated (Dyson) positron quasipartice wave function of Eq. (1). A the diagrams have equa mirror images. C. Anguar momentum convergence The use of B-spine basis sets provides for a fast convergence of the perturbation-theory sums in the sef-energy and Z diagrams with respect to the number of intermediate eectron and positron states in a given partia wave with anguar momentum. However, the sums in the diagrams are restricted to a finite number of partia waves up to a maximum orbita anguar momentum max, and the question of convergence with respect to max needs to be addressed. One soution successfuy tested in Ref. [17] is to use extrapoation described by the asymptotic formuae [49], δ (k) = δ [max] A (k) (k) + ( max + 1/2) 3, (22) Z (k) = Z [max] (k) + B (k) max + 1/2, (23) where δ [max] (k) and Z [max] (k) are the phase shift and annihiation parameter obtained for a given max, and A (k) and B (k) are constants specific to a particuar coision target, positron partia wave, and momentum k. These constants and the extrapoated vaues of the

7 7 phase shift and Z are determined by fitting δ [max] (k) and Z [max] (k) over a range of max to Eqs. (22) and (23), respectivey. The use of Eqs. (22) and (23) to extrapoate the phase shifts and Z vaues to max is iustrated in Fig. 6 for xenon. It shows that the numerica cacuations adhere cosey to the asymptotic form for max = 7 1, aowing a reiabe extrapoation to be made. This aso indicates that athough the extrapoation formuae are derived using perturbation theory [49], their use is aso vaid for the non-perturbative cacuations presented here. Note that extrapoation is particuary important for Z, where it contributes up to 3% of the tota. The roe of high intermediate orbita anguar momenta in Z is arge because the annihiation probabiity is sensitive to sma eectron-positron separations at the point of coaescence. Phase shift (radian) Xe s (+1/2) 3 2 s p p d p d s d Xe D. Partia-wave convergence of eastic scattering cross sections 15 s The eastic scattering cross section is obtained as a sum over the partia waves [51], σ e = 4π k 2 (2 + 1) sin 2 δ. (24) = At ow positron energies ony a few partia waves contribute to σ e, and the contributions decrease quicky with. On the other hand, the contribution of higher partia waves is more important in the differentia eastic cross section, where f(θ) = 1 2ik dσ e dω = f(θ) 2, (25) (2 + 1)(e 2iδ 1)P (cos θ), (26) = is the scattering ampitude and P (cos θ) are Legendre poynomias. Here arge interfere constructivey at sma scattering anges θ (due to the ong-range poarization potentia), producing a characteristic cusp at θ =. As ony s-, p- and d-wave phase shifts have been cacuated in the present work, some way must be found of accounting for the higher partia waves. This is done by noting that for higher partia waves, the dipoe term in Σ (2) ε + Σ (3) ε dominates the sef-energy at ow energies. At arge distances it corresponds to the oca energyindependent poarization potentia α d /2r 4. It aters the ow-energy ective range expansion of the scatter- Z 1 p p 5 d s s d d (+1/2) 1 FIG. 6. Extrapoation of phase shifts and Z for xenon with s-, p and d-wave incident positrons (back, red and bue symbos, respectivey) of momenta: k =.2 a.u. (circes); k =.4 a.u. (squares); and k =.6 a.u. (diamonds). Symbos show vaues obtained with max = 7 1, dotted ines are shown as a guide ony, and dashed ines show extrapoation to max. ing phase shifts [5], [ tan δ ak 1 πα dk 3a 4α dk 2 ( )] 1 αd k n C, 3 4 (27) δ πα d k 2, ( 1), (28) (2 1)(2 + 1)(2 + 3) where a is the scattering ength and C is a positive constant. We appy Eq. (27) to extract the scattering ength from the numerica s-wave phase shifts (see Sec. IV B), and use Eq. (28) for 3 in the cacuations of the differentia and tota eastic cross sections. p

8 8 IV. RESULTS: SCATTERING A. Phase shifts Eastic scattering phase shifts for the nobe gas atoms are tabuated in appendix B. The genera features of the phase shifts as functions of the positron momentum k are iustrated for krypton in Fig. 7. In the HF (staticfied) approximation, the phase shifts are negative, indicating a repusive eectrostatic fied, as is expected for positrons. Incusion of the second-order correation potentia Σ (2) ε eads to an attractive positron-atom potentia at ong range, making the phase shifts positive for ow k (dashed curves in Fig. 7). The asymptotic form of Σ (2) ε (i.e., α d /2r 4 ) eads to terms quadratic in k in the owenergy expansions (27) and (28). As a resut, the s-wave phase shift reaches a maximum and then fa off with increasing k, passing through zero (Ramsaur-Townsend ect) to negative phase shifts at higher k. The higherorder contributions to the correation potentia (Σ (3) ε and Σ (Γ) ε ) have opposing ects. Incusion of third-order screening diagrams Σ (3) ε decreases the strength of the positron-atom potentia and reduces the phase shifts. The contribution of virtua positronium formation Σ (Γ) ε is the greater of the two. It increases the strength of the positron-atom potentia, giving a particuary arge contribution at higher k for p and d waves. In Fig. 7 we aso show the phase shifts from the poarized-orbita cacuations of McEachran et a. [52] (squares). The poarized-orbita approximation makes a number of drastic assumptions. It considers a inear response of the target to the fied of a stationary (i.e., infinitey massive) partice, and drops the monopoe poarization term in the potentia. Unike Σ ε, the poarizedorbita potentia is oca and energy independent, and it does not account for the (nonperturbative) contribution of virtua Ps formation. It is thus remarkabe, and ikey fortuitous, that the poarized orbita cacuations give the s-wave phase shift in such cose agreement with the many-body cacuation. For p and d waves, however, the distinct ect of virtua Ps formation for k >.4 a.u. produces phase shifts that are 1 2% greater than those from the poarized-orbita cacuation. This behaviour of the phase shifts, incuding the Ramsauer-Townsend minimum in the s-wave scattering, is observed for a nobe gas atoms. Quantitativey, the correation ects [i.e., the contribution of δ (k) to the phase shift (8)] become progressivey arger from He to Xe. It was this increase in positron-atom correationa attraction that ed to predictions of positron binding to neutra atoms [22]. Phase shift (rad.) Phase shift (rad.) Phase shift (rad.) Kr s-wave Positron momentum k (a.u.) Kr p-wave Positron momentum k (a.u.) Kr d-wave Positron momentum k (a.u.) FIG. 7. Scattering phase shifts for s-, p- and d-wave positrons on Kr in various approximations: HF (static-fied) approximation (dotted curve); HF pus second-order correation potentia Σ (2) ε (dashed curve); HF + Σ (2) ε + Σ (3) ε (dot-dotdash curve); HF + Σ (2) ε + Σ (Γ) ε (dot-dashed curve); and tota: HF + Σ (2) ε + Σ (3) ε + Σ (Γ) ε (soid curve). Squares are theoretica resuts cacuated using the poarized orbita method [52].

9 9 B. Scattering engths A singe quantity that characterizes the strength of positron-atom attraction at ow energies is the scattering ength a. It can be extracted from the ective-range expansion (27) of the s-wave phase shift, written as δ (k) ak πα dk 2 3 4aα d αd k 3 k3 n + Dk 3, (29) 4 where D is a constant. Equation (29) is convenient for systems in which the scattering ength is not too arge. It works we for heium, neon, argon and krypton, and we use it as a two-parameter fit over the range of momenta k =.2.6 a.u. [53]. The corresponding vaues of the scattering ength are shown in Tabe I. TABLE I. Positron scattering engths a in a.u. for the nobe gas atoms. He Ne Ar Kr Xe a a a b a c a d.48 a e a f 4.9 ± ± ± 18.4 a Present many-body cacuations. b Previous many-body cacuations [23]. c Poarized orbita cacuations [52]. d Kohn variationa cacuations [54]. e Convergent cose-couping cacuations [4]. f Experiment: Ar [35], Kr [34], Xe [36]. Compared with other nobe-gas atoms, the s-wave phase shift for xenon is arge (see Tabe VII), indicating a much greater scattering ength. In this case it is more convenient to anayse the behaviour of the phase shift using Eq. (27) in the form k cot δ 1 a + πα dk 3a 2 + 4α dk 2 ( n C 3a αd k 4 ). (3) Figure 8 shows the dependence of k cot δ on the positron momentum k, together with a two-parameter fits using Eq. (3), and a three-parameter fit in which a cubic term Dk 3 is added on the right-hand-side of Eq. (3). The two-parameter fit gives a = 86.6 a.u., whie the three parameter fit gives a = 82.4 a.u. Given the uncertainty of the fitting procedure, our predicted scattering ength for Xe is a = 84.5 ± 2 a.u. In fact, the exact vaue of the scattering ength for xenon is very sensitive to the positron-atom correation potentia. When the scattering ength is arge, its reciproca κ = 1/a is known to vary ineary with the strength of the potentia [51]. We can thus compare our theoretica prediction κ =.118 a.u. with a typica vaue κ.5 a.u., compatibe with the radius of the Xe atom. (For exampe, the positron scattering k cot δ Positron momentum (a.u.) FIG. 8. s-wave scattering phase shift for Xe. Circes show vaues from the present many-body theory cacuation; soid ine is a two-parameter fit by Eq. (3) at k =.2 and.4 a.u.; dashed ine is a three-parameter fit [Eq. (3) with a cubic term Dk 3 added] using k =.2.6 a.u. ength in the static fied of the Xe atom is a = 1.93 a.u., which corresponds to κ =.518 a.u.) This shows that predicting κ (and hence, the scattering ength) with 1% accuracy requires a better than.1% accuracy in the cacuation of the correation potentia, which is probaby beyond any method for such a compex many-eectron target as Xe. If the scattering ength a is arge and negative then a virtua state exists at the energy ɛ = 2 /2ma 2, where m is the positron mass [51]. It gives rise to enhanced eastic scattering and annihiation ( a 2 ) at ow positron energies. Tabe I compares the scattering engths extracted from the s-wave phase shifts with other theoretica and experimenta resuts. Heium and neon dispay cose scattering engths in the present cacuation, with the resut for neon in good agreement with previous many-body theory cacuations [23], athough somewhat ower than other theoretica predictions [4, 52, 54]. For argon and krypton, there is cose accord between the present manybody theory, poarized orbita [52], and CCC [4] cacuations, and the experimenta resuts of Zecca et a. [34, 35]. The scattering ength increases across the nobe-gas atom sequence, giving rise to a virtua s-wave eve for xenon with the energy of approximatey 2 mev. As mentioned above, this indicates that scattering cacuations in the ow-energy region for xenon wi dispay a high sensitivity to the representation of the correation potentia. The many-body theory scattering ength cacuated here for xenon is a factor of two greater than that from the poarized-orbita cacuations [52], but somewhat ower than those obtained in the earier many-body cacuations [23] and the CCC cacuations [4]. However, it is compatibe with the experimenta vaue determined by extrapoating the measured ow-energy cross section with the aid of CCC cacuations [36].

10 1 C. Eastic scattering cross sections The eastic scattering cross sections aong the nobe gas sequence are shown in Figs. 9 13, where we compare them with existing experimenta and theoretica data. The numerica cross sections are tabuated in appendix B. For heium (Fig. 9) our many-body theory cacuations agree cosey with the variationa cacuations [55], the convergent cose-couping resuts [56] and the earier experimenta measurements of Refs. [57, 58], as we as the most recent ones [33]. A consensus on the eastic scattering cross section appears to have been reached. The poarized orbita resuts of [52] are not in agreement with other theoretica and experimenta data, and the magnetic-fied-free experimenta measurements of Nagumo et a. [61] are arger than a theory and experimenta resuts. For neon (Fig. 1), examining theoretica data first, the present many-body theory resuts agree most cosey with the reativistic poarized orbita resuts [38]. At energies beow 2 ev, the present many-body theory resuts agree we with the previous many-body theory cacuations [23], but trend ower above this energy. The convergent cose-couping resuts of [4] are higher than the other cacuations, whie the poarized orbita resuts of [52] are consideraby ower than both experiment and theory at energies above 2 ev. Comparing with experimenta data, the present many-body theory resuts agree most cosey with the recent measurements of Suivan et a. [37] above 2 ev, but are ower than Suivan s resuts beow this energy, where they agree better with the measurements of Stein et a. [57]. For argon (Fig. 11) the many-body theory resuts agree we with the convergent cose-couping resuts [4], but are higher than the nonreativistic and reativistic poarized orbita resuts [39, 52] above 2 ev. Comparing to experiment, the many-body theory resuts are in good agreement with the measurements of Sinapius et a. [62], Zecca et a. [35] and Suivan et a. [37], athough the atter two give sighty higher vaues at most energies shown. Note that the Ramsauer-Townsend minimum, which is very prominent in He and Ne, is not visibe in argon. This is a resut of the shift of the minimum in the s-wave scattering cross section towards higher energies, where it gets fied with higher partia wave contributions. This interpay of the contributions of different partia waves produces a characteristic pateau in the cross section, which stretches from 2 ev to the Ps formation threshod at 8.96 ev. For krypton, the convergent cose-couping resuts of [4] are in good agreement with the many-body theory resuts, whie the poarized orbita resuts (both reativistic [52] and nonreativistic [38]) are ower than other theoretica predictions above 1 ev. The experimenta data of Sinapius et a. [62], and the more recent measurements by Zecca et a. [34] and Suivan et a. [38] (above 2 ev) are in good agreement with each-other and with the present many-body theory resuts. The measurements of Dababneh et a. [66] and Coeman et a. [6] ie beow the manybody theory resuts. Compared with argon, the arger s-wave cross section at ow energies and the increased contributions of higher partia waves produces a steadiy decreasing tota cross section, with no trace of the underying Ramsauer-Townsend minimum. For xenon, the present many-body theory resuts are in exceent agreement with the convergent cose-couping resuts of [4], ower than the previous many-body theory cacuations of [23] (which empoyed an approximate treatment of virtua Ps formation), and somewhat higher than the nonreativistic and reativistic variants of the poarized orbita method [39, 52]. Compared to experiment the present resuts agree most cosey with the measurements of Sinapius et a. [62] and Suivan et a, especiay towards higher energies. The experimenta resuts of Dababneh et a. [66] and Coeman et a. [6] ie beow the other experimenta resuts and the present manybody theory cacuations. A possibe reason for this may ie in the fact that the differentia scattering cross section for positron scattering on heavier nobe gas atoms is strongy forward peaked. Poor detection of forwardscattered positrons [63] wi cause an underestimate in the cross section. Some genera trends can be seen in the eastic scattering cross sections across the nobe gas sequence. The many-body theory resuts are in good agreement with non-perturbative convergent cose-couping resuts [4], apart from neon where there is a discrepancy that wi need further investigation. The poarized orbita resuts [52] are seen to underestimate the cross sections at higher energies, ikey due to the negect of virtua positronium formation and the use of energy-independent correation potentia. Agreement with recent experimenta measurements [33 39] is generay cose. D. Differentia cross sections A quantity more sensitive to the accuracy of a scattering cacuation is the differentia cross section, Eq. (25). By making use of Eq. (28) for higher partia waves, the scattering ampitude (26) can be written as, f(θ) = = [ e 2iδ 1 πα d k 2 ] (2 + 1) P (cos θ) πα dk sin θ 2ik (2 1)(2 + 1)(2 + 3) 2 2, (31)

11 11 Cross Section (1-16 cm 2 ) He Energy (ev) FIG. 9. Eastic scattering cross section for He. Theory: soid, present many-body theory; dot-dashed, poarized orbita, [52]; dashed, MBPT [23]; dotted, CCC [56]; dot-dot-dashed, Kohn variationa [55]. Experiment: squares, Stein et a. [57]; circes, Mizogawa et a. [58]; stars, Karwasz et a. [59]; diamonds, Coeman et a. [6]; crosses, Suivan et a. [33]; trianges, Nagumo et a. [61]. 1.2 Cross Section (1-16 cm 2 ) Ne Energy (ev) FIG. 1. Eastic scattering cross section for Ne. Theory: soid, present many-body theory; dot-dashed, poarized orbita, [52]; dashed, MBPT [23]; dotted, CCC [4]; dot-dot dashed, reativistic poarized orbita [37]. Experiment: squares, Stein et a. [57]; circes, Sinapius et a. [62]; diamonds, Coeman et a. [6]; crosses, Suivan et a. [37]; trianges, Nagumo et a. [85].

12 12 Cross Section (1-16 cm 2 ) Ar Energy (ev) FIG. 11. Eastic scattering cross section for Ar. Theory: soid, present many-body theory; dot-dashed, poarized orbita, [52]; dashed, MBPT [23]; dotted, CCC [4]; dot-dot-dashed, reativistic poarized orbita [37]. Experiment: squares, Stein et a. [64]; circes, Sinapius et a. [62]; stars, Karwasz et a. [65]; diamonds, Coeman et a. [6]; trianges, Zecca et a. [35]; crosses, Suivan et a. [37]. 4 Cross Section (1-16 cm 2 ) Kr Energy (ev) FIG. 12. Eastic scattering cross section for Kr. Theory: soid, present many-body theory; dot-dashed, poarized orbita, [52]; dashed, MBPT [23]; dotted, CCC [4]; dot-dot-dashed, reativistic poarized orbita [38]. Experiment: squares, Dababneh et a. [66]; circes, Sinapius et a. [62]; diamonds, Coeman et a. [6]; trianges, Zecca et a. [34]; crosses, Suivan et a. [38].

13 13 6 Cross Section (1-16 cm 2 ) Xe Energy (ev) FIG. 13. Eastic scattering cross section for Xe. Theory: soid, many-body theory; dot-dashed, poarized orbita [52]; dashed, Gribakin et a. [23]; dotted, CCC [4]; dot-dot-dashed, reativistic poarized orbita [39]. Experiment: squares, Dababneh et a. [66]; circes, Sinapius et a. [62]; crosses, Suivan et a. [39]; diamonds, Coeman et a. [6]. where is the maximum partia wave for which the phase shift has been cacuated expicity. This procedure, or some other way of ectivey summing over a partia waves up to infinity, is necessary to describe the cusp of the differentia cross section at θ =, which is due to the ong-range ( α d /2r 4 ) poarization potentia. Figures 14 and 15 compare differentia cross sections cacuated using our many-body phase shifts for krypton and xenon at an incident positron energy of 2 ev with other theoretica cacuations and experimenta data from the San Diego [67, 68] and ANU groups [38, 39]. It is important to note that these groups use positron traps to accumuate positrons, which are then extracted to form a pused, energy tunabe positron beam. This beam is magneticay guided to the target gas ce, and the differentia cross section is measured by observing the change in the ongitudina positron energy. However, the apparatus is not abe to distinguish between forward-scattered and back-scattered partices, as any back-scattered partices are refected and passed back through the gas ce. The measured differentia cross section for a scattering ange θ is therefore the sum of the cross sections at the ange θ and 18 θ. Theoretica resuts have therefore been foded about 9 where necessary. For krypton (Fig. 14) the present many-body theory cacuations are in good agreement with the CCC resuts [34], nonreativistic and reativistic variants of the poarized orbita method [38, 52], and the experimenta data of [67] and [38] across the anguar range, with the CCC vaues trending sighty higher at sma anges. (The atter is compatibe with the arger absoute vaue of the scattering ength in the CCC cacuation, Tabe I.) However, a arge discrepancy is observed with the poarization cacuations of [69]. In the case of xenon (Fig. 15), there is exceent agreement between the present many-body cacuations and the CCC cacuations [39], with both cacuations being somewhat higher than the reativistic poarized orbita cacuations [39]. As in krypton, the present resuts show a arge difference with the poarization potentia cacuations of [69]. The many-body theory and CCC cacuations are in better agreement with the ANU experiment [39] than the data of Ref. [68], with the two experiments in better accord above 4. The stronger peaking of both many-body theory and CCC resuts at θ = can be reated to the arger vaues of the scattering engths in these cacuations compared to the poarizedorbita resut (Tabe I). The arge scattering ength in Xe is behind the strongy enhanced positron annihiation rates at ow (e.g., room-temperature, therma) energies in

14 14 DCS (1-16 cm 2 /sr) Kr 2 ev Ange (deg) FIG. 14. Differentia eastic scattering cross section for Kr. Theory: soid, present many-body theory; dotted, poarized orbita, [52]; dashed, CCC [34]; dot-dashed, reativistic poarized orbita [38]; dashed-dashed-dotted, poarization potentia [69]. Experiment: circes, Gibert et a. [67]; squares, Suivan et a. [38]. DCS (1-16 cm 2 /sr) Xe 2 ev Ange (deg) FIG. 15. Differentia eastic scattering cross section for Xe. Theory: soid, present many-body theory; dashed, CCC [39]; dot-dashed, reativistic poarized orbita [39]; dashed-dashed-dotted, poarization potentia [69]. Experiment: circes, Marer et a. [68]; squares, Suivan et a. [39].

15 15 xenon (see Sec. V). We, therefore, beieve that the manybody theory, CCC and ANU data are more accurate for θ < 3. a fraction of one per cent. In addition, for neon about.3% is due to positron annihiation with the core (1s) eectrons (Tabe IV). V. POSITRON ANNIHILATION A. Energy resoved Z Figure 16 shows partia-wave contributions to Z for positron annihiation on the vaence she eectrons of the nobe gas atoms as functions of the positron momentum. Note that at ow positron momenta k the Wigner threshod aw predicts Z k 2 for the positron with the orbita anguar momentum [51]. As a resut, the s-wave contribution dominates at ow energies. As one moves aong the nobe gas sequence, the s-wave Z becomes increasingy arge and strongy peaked at ow energies. This strong energy dependence is due to the existence of positron-atom virtua eves [5, 23, 25], which is signified by arge scattering engths a for Ar, Kr and Xe. In this case the momentum dependence of Z at ow energies can be described anayticay [23, 25, 7 72], as Z (k) = K κ 2 + A, (32) + k2 where K and A are constants. The first term in Eq. (32) is due to the s-wave contribution enhanced by the virtua abe (i.e., sma κ = 1/a, κ 1/R a 1 a.u., where R a is the atomic radius). The constant term A accounts for the nonresonant background and contributions of higher partia waves to Z. Equation (32) makes it cear that if the therma positron momentum k is smaer than κ then the annihiation rate wi be proportiona to the scattering ength squared. The p-wave Z aso appears to show the formation of a broad shape resonance especiay prominent for Xe. This is supported by the behaviour of the p-wave phase shift, e.g., that for Kr shown in Fig. 7. For more poarizabe targets which generate stronger positron attraction, such as Mg and Zn, p-wave resonances in Z become a prominent feature of the energy dependence of the annihiation parameter [73]. B. Thermay averaged Z Most of the avaiabe experimenta data for Z in nobe gases have been obtained for thermaized positrons at room temperature [74]. As can be seen from Fig. 16, the annihiation rates in heium and neon have a weak energy dependence in the range of therma positron momenta (k.45 a.u. at room temperature). Hence, for these atoms we take the Z vaues at k =.4 a.u. to compare with experiment (see Tabe II). The Z vaues at such ow momenta are primariy due to the s-wave positron annihiation, with the p-wave contributing ony TABLE II. Thermay averaged Z for the nobe gases. a Z b Z c Z d Z He Ne Ar Kr Xe , a Present many-body theory cacuations. b He: Kohn variationa cacuations [75]; Ne, Ar, Kr, Xe: poarized orbita cacuations [52]. c Experiment: He, Ne, Ar [6]; Kr [76], Xe [76, 77]. d Experiment: Ar and Kr [78]; Xe [7]. For argon, krypton and xenon the energy dependence of Z becomes progressivey stronger and we compute therma Z vaues from the Maxweian average Z = Z (k) exp( k2 /2k B T ) (2πk B T ) 3/2 4πk 2 dk, (33) where k B is the Botzmann constant and k B T = a.u. at room temperature T = 293 K. The integration is performed using a fit of the cacuated tota Z vaues (s-, p-, and d-wave, vaence and core), of the form Z (k) = K κ 2 + k 2 + A, (34) + βk4 over the momentum range k =.2.3 a.u. This fit is based on Eq. (32), with an extra k 4 term incuded to improve the accuracy. This is especiay important for xenon where Z has the most vigorous momentum dependence. Figure 17 shows the cacuated Z for Xe (open circes) [79], together with two fits of the form (34), experimenta data from Ref. [74], and Z from the poarized orbita cacuations [52]. In the first fit shown in Fig. 17 by the dashed ine, the parameter κ is fixed by the vaue of the scattering ength, κ = 1/a =.118 a.u., and the vaues of the other parameters are K =.683, A = 3.95 and β = a.u. This fit gives Z = 448 at T = 293 K. In the second fit shown by the soid ine, κ is regarded as a free parameter. The corresponding set of vaues K =.647, κ =.914, A = 26.74, and β = 5.49 a.u. produces an exceent fit of the numerica Z over the whoe momentum range. It gives Z = 458. Given the 2% difference between the two vaues, we report Z = 45 as our best prediction of the room-temperature vaue for Xe. Using Eq. (34) as a four-parameter fit for Ar (K =.22, κ =.148, A = 11.26, and β = a.u.) yieds Z = 26. at T = 293 K. The fit for Kr (K =.3659, κ =.696, A = 16.4, and β = 3.74 a.u.) produces Z = Athough the momentum dependence of Z in Ar and Kr is not neary as steep as in Xe,

16 s-wave 1 2 Z Positron momentum (a.u.) p-wave Z Positron momentum (a.u.) 15 d-wave 1 Z Positron momentum (a.u.) FIG. 16. Contributions of the s-, p- and d-waves to the annihiation rate Z on the vaence subshes of the nobe gases: heium (1s, soid ine); neon (2s + 2p, dotted ine); argon (3s + 3p, dashed ine); krypton (4s + 4p, dash-dotted ine) and xenon (5s + 5p, dash-dash-dotted ine). therma averaging is important for them. For exampe, the cacuated vaues at k =.4 a.u. are Z = 27.1 for Ar and 72.1 for Kr, whie the thermay averaged vaues obtained above are cose to Z (k) at k.48 a.u.

17 Xe 1 3 Z Positron momentum (a.u.) FIG. 17. Cacuated Z for Xe (sum of s-, p- and d-wave vaues incuding both vaence and core contributions, open circes) and their fits (34), in which κ = 1/a =.118 a.u. is fixed (dashed ine), or used as a fitting parameter (soid ine). Soid circes are experimenta vaues from Ref. [74]. Trianges with dot-dashed ine are the poarized orbita resuts [52]. The present many-body theory and the Z vaues obtained above are a significant improvement on the previous many-body theory study [23]. In Tabe II the cacuated Z vaues for the nobe gas atoms are compared with avaiabe experimenta data and some theory vaues. For heium the many-body theory resut is in good agreement with precise variationa cacuations [75] and with the measurements [6], the discrepancy being ess than 5%. For neon, argon and krypton the many-body resuts are in cose agreement with the earier measurements [6], but differ significanty from the positron-trap resuts for Ar and Kr [78], and are ower than the poarized orbita cacuations [52]. Even assuming an error of 5 1% in the many-body cacuations, the resuts of Iwata et a. [78] for Ar and Kr appear to be anomaousy high. For xenon, the many-body resut is much higher than that from the poarized orbita cacuations [52] and is in reasonabe accord with the experiment of [7] and one set of earier measurements [76]. Given the approximations made in the poarized orbita method, the discrepancy between its resuts and our theory are not unexpected. (In fact it is remarkabe that this reativey crude method produces such reasonabe resuts.) To gain a better understanding of the discrepancies between theory and experiment, it is necessary to examine the experimenta techniques used to measure Z. The resuts of Refs. [6, 76, 77] were obtained using positron ifetime spectroscopy. In this technique, positrons from a radioactive source (e.g., 22 Na) are injected into a gas ce, where they thermaize and annihiate. The ifetime is measured as the deay between the nucear gamma ray emitted in the β + decay, and the annihiation gamma ray. By measuring the ifetime of positrons it is possibe to obtain a ifetime spectrum. The annihiation rate can then be found from the exponentia fit of this spectrum. For ighter nobes gases (He to Kr) thermaization occurs much faster than annihiation, which provides for a reiabe measurement of Z with fuy thermaized positrons. In xenon the positron energy oss in momentum-transfer coisions is smaer due to the arge mass of the atom, whie the Z vaues in Xe are higher. As a resut, the annihiation rate is measured for epitherma positrons, resuting in Z vaues ower than expected (e.g., Z = 32 [77]). Adding sma amounts of a ighter, ow-z gas, e.g., He or H 2, to Xe aows for fast thermaization and produces truy thermaized annihiation rates with Z = 4 45 [76]. As seen in Tabe II, these vaues are in good agreement with the present cacuation. The resuts of Refs. [7, 78] were obtained using a Penning-Mamberg positron trap. In this type of experiment positrons from a radioactive source are sowed