Positron Scattering from Atoms and Molecules

A review on the positron scattering from atoms and molecules is presented in this article. The focus on positron scattering studies is on the rise due to their presence in various fields and application of cross section data in such environments. Positron scattering is usually investigated using theoretical approaches that are similar to those for electron scattering, being its anti-particle. However, most experimental or theoretical studies are limited to the investigation of electron and positron scattering from inert gases, single electron systems and simple or symmetric molecules. Optical potential and polarized orbital approaches are the widely used methods for investigating positron scattering from atoms. Close coupling approach has also been used for scattering from atoms, but for lighter targets with low energy projectiles. The theoretical approaches have been quite successful in predicting cross sections and agree reasonably well with experimental measurements. The comparison is generally good for electrons for both elastic and inelastic scatterings cross sections, while spin polarization has been critical due to its sensitive perturbing interaction. Positron scattering cross sections show relatively less features than that of electron scattering. The features of positron impact elastic scattering have been consistent with experiment, while total cross section requires significant improvement. For scattering from molecules, utilization of both spherical complex optical potential and R-matrix methods have proved to be efficient in predicting cross sections in their respective energy ranges. The results obtained shows reasonable comparison with most of the existing data, wherever available. In the present article we illustrate these findings with a list of comprehensive references to data sources, albeit not exhaustive.


Introduction
A positron, the antiparticle of the electron, has the same mass, electric charge (but positive) and spin (1/2) as that of an electron. Like other antiparticles, positrons were produced during the period of baryogenesis when the universe was extremely hot and dense, but now they exist in much lower numbers than its counter part, the electrons. Although not found in normal conditions, they are produced at the galatic center or supernovae events and are found in copious amount in cosmic ray showers and in the ionosphere. Positrons are created naturally in β + radioactive decays such as from K-40, particle reactions or by pair production from a sufficiently energetic photon interacting with the atomic nuclei in a material. Nevertheless, a small percent of potassium (0.0117%) K-40 is the single most abundant radioisotope in the human body and produces about 4000 natural positrons per day. However, soon after its creation, they annihilate with electrons or form the exotic atom, positronium (Ps), with a very short lifetime, finally decaying to 2 or 3 gamma rays each with energy 511 keV. Ps has to calculate the polarizability of atoms [53,54] and the scattered wave was obtained from a potential scattering problem. The method showed good agreement with measured cross sections for for He and Ne atoms [17,21]. McEachran et al. [55] implemented their method successfully for other atoms as well.
One of the most rigorous approach for positron-ion scattering has been the Feshbach projection operator method [56], where the usual Hartree-Fock and exchange potentials are augmented by an optical potential [47]. However, the method employs correlation functions that are of Hylleras type and hence do not include long range correlation.
Kohn Variational Principle (KVP) (e.g., Reference [57]) is usually applied to low-energy positron scattering to obtain elastic and Ps formation cross sections. In this method a two-component trial wave function is chosen, having the correct asymptotic form with enough flexibility to describe all the short-range distortions and correlations of the positron-atom system. This wave function is then used in the Kohn functional, which can be written in terms of the K-matrix elements. From the K-matrix, cross sections can be calculated.
The many-body-theory for a positron-atom interacting system (e.g., References [58,59]) is based on the Dyson equation. This is solved by the representation of eigenfunctions of the Hartree-Fock Hamiltonian. The formulated self-energy matrix gives the phase shifts. This approach is used to study low energy positron scattering from atoms.
Schwinger Multichannel Calculations (SMC) is a well-known method to study low energy scattering (e.g., References [60,61]). The backbone of the method is the computation of variational expression for the scattering amplitude. SMC describes target polarization through single virtual excitations of the target wave function, explicitly considered in the expansion of the scattering wave function. The Lippmann-Schwinger scattering equations are then solved to obtain the cross sections [62].
The other method of interest is the close-coupling (CC) approach for the scattered wave and the R-matrix method [63,64] to solve the coupled set of integro-differential equations. Jones et al. [26] used the convergent close-coupling (CCC) approximation, where they solve the equations with a different set of codes than standard R-matrix codes. They use multi-configuration Dirac-Fock program of Grant et al. [65] to obtain the target wave functions. Their results for positrons scattering from Ne and Ar indicate that, while both polarized orbital method and CCC approximation showed good agreement with experiment in general, the polarized orbital method yielded slightly better cross sections. In the standard close-coupling formalism for molecules (e.g., Reference [66]), the convergence in the expansion of the three-body wave function is obtained using the exact discrete eigenstates of the atomic target. The technique relies on the expansion of the total wave function in the set of target states of the atom and Ps. The CCC method allows for the examination of the effect of virtual excitation to the continuum as well [67].
For molecules, the optical potential method (e.g., References [49,68]), Born approximations (e.g., References [69,70]), and so forth are the most common quantum mechanical perturbative theories used for electron scattering presently. The positron collision studies is an extension of the optical potential method [71]. In case of non-perturbative theories, the close-coupling or the grid-based method for solving Schrödinger equation are employed. Irrespective of whether the method is perturbative or non-perturbative, the positron-molecule scattering is an extension to the positron-atom interaction technique. One has to consider multi-centre approach to deal with the projectile-target interaction due to the absence of spherical symmetry and due to the complexity of molecules. Here we will discuss few of the most commonly used theoretical methods to deal with positron-atom/molecule interaction.
Theoretical methods to investigate positron scattering and various target molecules studied by these approaches along with references are listed in Table 1. The list given below is not exhaustive, but gives an overall picture of various studies done so far. Further, this review will elaborate the most common approach, the optical potential method, to study electron and positron scattering from atoms and molecules with reasonable success. Table 1.

Method
Target (

Scattering Parameters: Cross Sections and Spin Polarizations for Atoms
The characteristic features of the scattering can be observed in the cross section and in spin polarization caused by the projectile. While cross section can be obtained by solving non-relativistic Schrödinger or relativistic Dirac or Dirac-Fock equations, the latter provides accurate treatment for spin polarization parameters. In the present review, we will present relativistic single particle Dirac approach, which has been successful in reproducing the observed scattering phenomena.
The relativistic Dirac equation for a projection of rest mass m o and velocity v traveling in a central field V(r) is given by (e.g., Reference [180,181]), where α and β are the usual 4 × 4 Dirac matrices and ψ is a four-component (spinor) function, ψ = (ψ 1 , ψ 2 , ψ 3 , ψ 4 ). (ψ 1 , ψ 2 ) are the large components and (ψ 3 , ψ 4 ) are the small components of ψ. Defining γ = (1 − v 2 /c 2 ) −1/2 , the total energy is E = mc 2 = m o γc 2 = E + m o c 2 where E is the kinetic energy, and writing the radial function of the large large component as G l = √ ηg l (r)/r, the equation for the large component can be rewritten as the Dirac equation reduced to the form similar to Schrödinger equation (e.g., Reference [9]) as where the effective Dirac potentials due to spin up and spin down respectively are, and The prime and double primes represent the first-and the second-order derivatives with respect to r, In atomic unit, m o = e =h = 1, 1/c = α, where α is the fine structure constant and hence γ = (1 + α 2 K 2 ) 1/2 and E = γc 2 = γ/α 2 . The proper solution of the Schrödinger like Dirac equation behaves asymptotically as, where j l and n l are spherical Bessel functions of the first and second kind respectively, and δ ± i are the phase shifts due to collisional interactions. The plus sign corresponds to the incident particles with spin up and the minus sign to those with spin downs. δ ± i indicates the shifts in the phase of the radial wave function due to the effect of interaction potentials in the scattering. The radial wave function will be "pushed out" if the potential is repulsive and vice-versa with respect to the incoming free radial wave. So from this quantity, we can determine various microscopic quantities like cross section. The values of δ ± i may be obtained from the values of g ± l at two adjacent points r and r + h (h << r) as .
The wave functions g ± l can be obtained by numerical integration of g ±" l using Numerov method and the spherical Bessel functions as described in Nahar and Wadehra [1]. Schrödinger/Dirac equation, can be solved by various other approaches mentioned above, such as, Kohn-Feshbach variational method [43], polarized orbitals method [2][3][4]46,47,51], close-coupling approximation [26,64] for the scattered wave function from which the phase shift is determined.
The generalized scattering amplitude for the collision process is given by [9], where andn is the unit vector perpendicular to the scattering plane. The differential cross section (DCS) for the scattering of the spin-1/2 particles by the spin zero neutral atom is given by where χ ν represents a spin state and P i =< χ ν |σ|χ ν > is the incident-beam polarization, which is assumed to be zero. The integrated elastic cross section for the unpolarized incident beam can be obtained as and the momentum transfer cross section by The integrated total cross section given by where S ± l = exp(2iδ ± l ). The integrated absorption cross section can be obtained from σ abs = σ tot − σ el Since the spin-orbit interaction is a short-range interaction, the phase shifts of the spin-up and the spin-down particles are equal (δ + l = δ − l ) for the large angular momenta lh. Hence for large l, g(θ) = 0 and the contribution to the scattering amplitude comes only from f (θ). If Born approximation is used for higher partial wave with l > M, f (θ) can be written as [9], where f B (K, θ) is the Born amplitude, S Bl = exp(2iδ Bl ) and δ Bl is the Born phase shift. The number of exact phase shifts to be evaluated depends on the impact energy before use of Born approximation. The contribution due to Born approximation should be small. At large distance the interaction potential V(r) is dominated by the long range part V LR (r) = −α d /2r 4 of the polarization potential and Born phase shift δ Bl . The interaction potential between the spin of the electron or positron and the orbital angular momentum L, which depends on the velocity and position vector with respect to the target atom, can cause the spin to orient. Hence, even with an unpolarized incident beam the orientation in a preferred direction can give a net spin polarization in the scattered beam. The amount of polarization produced due to the collision in the scattered beam is given by [182], The other two spin polarization parameters, T and U giving the angle of the component of the polarization vector in the scattering plane are given by [182], The three polarization parameters are interrelated through the condition P 2 + T 2 + U 2 = 1.

e ± and Target Atom Interaction Potential
To calculate the scattering parameters, we define the projectile-target interaction and a method to determine the respective wave functions. The scattering can be described in two general categories, elastic (where the total kinetic energy is conserved and the interaction potential is real) and inelastic (where part of the energy is lost due to absorption). For the inelastic processes such as excitation, ionization, positronium formation through electron capture and so forth the absorption potential is developed, which forms the imaginary part of the total complex potential. The total interaction potential between a neutral target (or a single atomic electron) and a projectile electron or positron is assumed to be symmetric or central, V(r), which depends on r only. The general form of V(r) is, where the real part V R (r) represents the elastic scattering and the imaginary part V A (r) represents the absorption of energy through the inelastic channels. When the total kinetic energy is conserved the imaginary part, V A (r), is zero. The absorption potential is negative and typically depends on the local density function.
V R (r) has several components: the averaged static potential V S (attractive for positrons and repulsive for electrons), V P polarization potential (attractive for both electrons and positrons) and an electron-electron exchange potential V ex (only for electrons). For a positron there is no exchange probability. The total real potential is represented as, The static potential, V S , is obtained by averaging the projectile-target interaction over the target wave function as, where ψ T is the asymmetric Hartree-Fock target wave function, Φ nlm (r) = φ nl (r)Y lm (r) are the partial atomic orbitals, e p is the projectile charge and N nlm is the occupancy number of the orbital (n, l, m). The radial part φ(r) of an orbital can be an analytic expansion, for example, tables of Clementi and Roetti [183] or in numerical form obtained from configuration interaction atomic structure calculation (e.g., Reference [65] or a Hylleras type wave function expansion (e.g., Reference [26]). Use of configuration interaction form is common in close coupling approximation. The polarization potential usually has a short and a long range part, where r c is the point where the two forms cross each other for the first time. The long range behavior is known to be of the form α o /r 4 where α o is the polarizability of the target. The short range form can vary. For the electrons scattering from a neutral atom, it could be the parameter free potential, for example, that given by O'Connel and Lane [184]. They developed the potential on the basis of energy dependent free-electron gas exchange potential and the energy-independent electron-gas correlation potential smoothly joining to the long-range polarization interaction and is given by, is the undistorted electronic charge density of the target. ρ(r) for the spherically symmetric atom is given by, where N nl is the occupancy number of the orbital (nl). For the positrons scattering the polarization potential can also be parameter free, such as that by Jain [185]. It is based on correlation energy of a single positron in a homogeneous electron gas with an asymptotic behavior of the long range polarization potential, and is given by, The long range form of V p ± (r) is given by V LR (r) = −α d /(2r 4 ), where α d is the static electric dipole polarizability. In polarized orbital method, the distortion due to polarization is incorporated in the wave function.
The exchange potential, V ex (r) is due to exchange between the projectile electron and the target electrons. One of the common form is given by Riley and Truhlar [186], where V D = V S + V P − is the direct interaction potential and ρ(r) is the radial density of the target. Chen et al. [187] introduced another type of potential which was used for elastic scattering from heavy inert gas, Kr, with reasonable success. When the impact energy becomes accessible for the inelastic processes (such as excitations of the target, positronium formation, etc.) absorption potential is introduced. The total absorption of energy has been represented by various model potentials with poor to good success for certain atoms (e.g., for Ar [9]). One major issue was the inclusion of various threshold energies for excitations and electron capture in case of positrons to form positronium. One successful absorption potential model, especially for electron scattering, has been the semi-empirical potential of Staszewska et al. [188,189], which are based on qualitative features of an absorption potential at short and long ranges in order to predict accurate differential cross sections. Their later model [189] has been in use considerably, and is given by, v is the local velocity of the projectile for (E − V R ) ≥ 0, ρ is the target electron density per unit volume andσ b is the average quasifree binary cross section for Pauli allowed electron-electron collisions and is obtained non-empirically by using the free-electron gas model for the target as, where p is the incident momentum of the projectile and k F is the target Fermi momentum. In their third version of V A , V.3, they define the parameters α and β as, where ∆ is the threshold energy for inelastic scattering and I is the ionization potential. The factor 1/2 in the equation is introduced to account for the exchange between the incident electrons and the atomic electrons of the target. The same absorption potential can be used for the positron scattering with the factor 1/2 removed, since there is no exchange effect during the positron scattering. The earlier version of Staszewska et al. [188] has also shown fair to good representation of absorption potential in reproducing the total cross sections. Various other absorption potential models are also available in literature, but have been only partly successful and hence need improvement. Figure 1 demonstrate the general features of various components of the real part V R (r) of the total projectile-target interaction potential V(r). The components are static potential (repulsive for electron and attractive for positrons), polarization potential (attractive for both e ± ), exchange potential (only for electrons and attractive) for e ± scattered by the cadmium atoms [8]. The static potential was obtained using Slater-type orbitals of Roothan-Hartree-Fock wave functions given by Clementi and Roetti [183].
The same orbital functions were used to obtain the electron density in the absorption potentials. As expected, the static potential dominates near to the nucleus and exchange potential starts away from it, but moves toward it with increasing energy of the projectile.

e ± Scattering from Molecules
Two of the most commonly used approaches to study electron/positron interaction from molecules are described below:

Optical Potential Approach
As mentioned above, the positron-target interaction is represented by an optical potential, which has a real part representing elastic processes and the imaginary part taking care of the loss of flux due to different inelastic channels. The calculations are generally carried out under the fixed nuclei (FN) approximation. This method has been successful in investigating not only atoms, but also larger molecules under multicentre approach. Various additivity rules are applied to find the cross section of a molecule from its atomic constituents. The independent atom model (IAM) is the simplest among them, where the cross section of the constituent atoms of the molecule is simply added to find the cross section for the molecule [72]. This approach was further modified and screening correction (IAM-SCAR) was applied to the cross section including the interference term and rotational excitation cross sections [190].
The methodology to obtain radial charge density and interaction potentials for atomic systems was already described earlier. The parametric form for atomic ρ and V st are then used to compute electron/positron impact cross section for atoms in the spherical complex optical potential (SCOP) approach [49]. In this method, the radial part of Schrödinger equation is solved through partial wave analysis to obtain the asymptotic solution, where δ l (k) is the phase shift, which carries the signature of interaction on the solution. To compute phase shift numerically, the radial Schrödinger equation is solved in the region r < a, where a is the finite range of the interaction potential. By applying the boundary condition, r = a we have, and its logarithmic derivative, where j l and η l are the spherical Bessel and Neumann functions respectively. The 'prime' means derivative with respect to 'r'. By inverting the equation as above we can obtain, The phase shift, δ l (k) obtained from above is then used to construct elastic and inelastic cross sections as [182], However, for molecules there is no direct method to obtain ρ and V st to be used to compute cross sections. Hence, various indirect techniques were adopted to approximate the charge density and static potential of molecules. For closed molecules like CH 4 , a single center expansion (SCE) is used, where the charge density of lighter atom is expanded from the center of mass of the heavier atom [191,192]. However, for larger non-spherical molecules, this is not possible. Antony and co-workers [50,193,194] have extensively used group additive rule to efficiently compute molecular charge density (and static potential). In this method the molecule is assumed to be made of different centers/groups having a group of atoms, each of them scattering positron independently. The grouping of atoms are performed based on the assumption that the charge density of one group do not overlap with the other. Even though such an approximation is a simplification of the actual scenario, it is sufficient to produce reasonable results in the present case (intermediate to high energy projectiles) without the loss of generality. Thus, for each scattering center the ρ and V st are developed and the desired cross section is computed for that center. Finally, cross sections calculated from each scattering center is added to get the molecular cross section. The total ionization cross section Q ion is obtained from the Q inel using complex scattering potential-ionization contribution (CSP-ic) technique [195,196]. Further details on the method is presented in the references.
In 1996 Reid and Wadehra [197] introduced an absorption model for positron collisions. The form of this potential is given as, In the above equation, a 0 is the Bohr radius and R is the Rydberg constant, δ = ∆ E F , = E i E F , f (x) = 2 δ x 3 + 6x + 3 ln −x +x and E F is the Fermi energy. ∆ defines a threshold below which all inelastic processes such as excitation and ionization are energetically forbidden. Determination of ∆ is important while dealing with positrons. Reid and Wadehra [197] in their original model used ∆ = ∆ p = I − 6.8 eV (where ∆ p is the Ps formation threshold energy). However, this cannot be used as a general rule, since excitations may occur below this energy. The parameter free model developed to find ∆ is given as [198], where E i is the incident energy and ∆ e is the lowest electronic excitation energy.

ab initio R-Matrix Method
The R-matrix method was developed to study scattering problems in nuclear physics by Wigner and others [199,200] in the late 1940s. Burke et al. [201] modified it to include electron-atom scattering processes and subsequently adapted to carry out accurate calculations for diatomic molecules [202]. This ab initio theory is primarily based on the concept of constructing two divisions in coordinate space (inner and outer regions) to adequately and efficiently represent the target and scattering regions. These regions are separated by a spherical boundary of radius 'a' containing the target molecule with N + 1 indistinguishable electrons (where +1 is due to the scattering electron). This N + 1 collision system acts as a bound state, where the short-range correlation effect and electron exchange among the electrons are dominant. This makes the inner region scattering problem quite complex. However, calculation is needed to be done only once as the problem is independent of energy. For an accurate solution of the inner region, Quantum chemistry codes are employed generating an energy independent wave function for the N + 1 electrons system. The results obtained in the inner region are fed to the outer region, which contains only the scattering electron. In the outer region, the scattering electron is far away from the center and thus exchange and correlation between electrons can be ignored. Therefore, only the long range multipolar interaction between the scattering electron and the target electron is considered here. For the simplicity of the numerical problem, single center approximation is employed, which converges quite rapidly without losing the generality of the problem. Since the calculations in the outer region are relatively simple, they are repeated for each electron energy on a finite grid to find desired scattering properties such as cross sections. In general, the inner region radius is varied from around 10 to 15 a.u. depending on the size of the target, while the outer region is kept infinite (100 a.u.). The value of R-matrix radius is chosen to accommodate total wave within the sphere.
The molecular orbitals are developed from atomic ones, which are expressed as basis functions centered on respective nuclei. The orbitals are a representation of the molecular charge density distribution, which must be negligibly small at the R-matrix boundary. The electrons are assumed to be in certain combinations of target orbitals to produce configuration state functions (CSF) in various total symmetries. A configuration interaction (CI) molecular wave function is expressed as a linear combination of CSFs in CI expansion. The molecular orbital representation in occupied and virtual orbitals are constructed using Hartree-Fock Self-Consistent Field (HF-SCF) method. For this, Gaussian-type orbitals (GTO) are used and the continuum orbitals of Faure et al. [203] are included up to g (l = 4). For dipole-forbidden excitations, ∆J = 1, J is rotational constant without spin coupling, the convergence is rapid for partial waves. However, for dipole-allowed excitations, ∆J = ±1, convergence is quite slow due to the long range dipole interactions.
The closed coupling expansion (CC) of inner region wave function under fixed nuclei approximation is expressed, where 'A' is the anti-symmetrization operator. u ij (x) are target continuum orbitals and Φ i (x) are target wave functions (together called 'target+continuum' configuration). Φ i are constructed from basic HF method and CI expansion is used for the other. The χ i are two-center quadratically integrable (L 2 ) functions. These are formulated with target occupied and virtual orbitals. Care has to be taken in choosing the correct basis function as final results are strongly depended on it. Since the continuum orbitals do not vanish at the R-matrix boundary (r = a), the Hamiltonian is modified by adding a surface or Bloch operator [204] to make it Hermitian inside the R-matrix sphere. The form of Bloch operator is, Then by diagonalizing the modified Hamiltonian (H N+1 + L N+1 ) the inner region electron wave function is obtained. The inner region solutions are used to set up the R-matrix at the boundary of the sphere. Then it is propagated to far region (a out ), where it is matched with the asymptotic functions obtained from the Gailitis expansion [205]. The numerical value of a out depends on the nature of the target. The quality of any scattering model depends on a satisfactory representation of the target. Therefore getting the right values of target parameters viz. ground state energy, vertical excitation energies, ground state dipole moment and rotational constant is a necessary precursor to scattering calculations.
The exchange and correlation effects are negligible in the outer region and the physical interactions dominate the scattering processes. The scattering electron moves under the influence of long range multipole potentials of the target. Hence, a single center expansion of the scattering wave function sufficient here. This is given by, where x j = (r j , σ j ) is the position and spin of the jth target electron, the functionsφ i are formed by coupling the scattering electron spin σ N+1 with the target state φ i and F i are reduced radial wave functions. By appropriate substitutions, a set of coupled, homogeneous, differential equations for the reduced radial wave function are obtained. The solution of this equation is obtained by propagating the R-matrix from boundary to sufficiently large distances, where the interaction between scattering electron and target molecule tends to zero [206]. Thereafter, asymptotic expansion techniques are used to solve for the outer region functions [207]. In the limit r → ∞ above equation have different linearly independent standing wave asymptotic solution j for each energetically open channels i given as, The coefficients K ij define the real, symmetric K-matrix, which contains all the scattering information. The eigenphase sum δ is used for the detection and parametrization of resonances, obtained directly from the diagonalized K-matrix, K D ij as δ = ∑ i arctan K D ij , where the summation is over the open channels. The scattering matrix S is a transformation of the K-matrix given by, The integral cross section for the excitation from states i to i [208] is given as, where S i is the spin angular momentum of the ith target state, S is the total spin angular momentum, Γ runs over the symmetry and l and l are orbital angular momentum quantum numbers corresponding to i to i states respectively. K-matrices form the input of POLYDCS program [209], from which the scattering observables are evaluated. For further details on R-matrix methodology to study positron collisions, please refer to References [115,118] and references therein.

Results and Discussions
Positron and electron scattering from neutral atoms/molecules and one or two electron ions have been investigated extensively as seen in the partial list of references in this article. For consistent features and feasible experimental set-ups, the targets have been mainly hydrogen atom, inert gases, He + and Li +,2+ ions, symmetric closed shell molecules, and few other simple molecules. Theoretical methods have been developed with significant accuracy for benchmarking the experiments. However, issues remain unresolved for better representation of the interaction and the wave function for sensitive cases. Examples of various features of the scattering cross sections and spin polarization for e ± from atoms and molecules are illustrated below.

Cross Sections and Spin Polarization for e ± Scattering from Atoms
Since the theoretical basis for e ± scattering from atoms is similar, reproduction of electron scattering cross sections gives a measure of accuracy of the method to study positron scattering [210,211]. However, its worth noting that the presence of Ps formation channel strongly affects the shape of the cross section in the energy range from its threshold to around 200 eV. Figure 2 left panel gives an example of accuracy in the features of differential cross sections (DCS) for scattering of electrons from argon, where measured DCS from 7 different experiments [18][19][20][22][23][24]212] are compared with theoretical predictions in nonrelativistic [1] and relativistic [9] approximations. Nahar and Wadehra obtained the static potential using analytic orbitals obtained by Clementi and Roetti [183] and absorption potential from Staszewska et al. [189]. Very good agreement in features, shape and magnitude benchmarks the theory. It may be noted that use of absorption potential has little effect on the DCS values as expected for elastic scattering and that relativistic corrections have removed some discrepancy at the minimum for DCS at 300 eV.
The potentials used for the cross sections, were used to study the angular dependence of spin polarization of the elastically scattered electrons from argon as demonstrated on the right panel of Figure 2 [9] at several electron impact energies. The spin-orbit polarization is highly sensitive to interaction potential. We see very good agreement between theory and experiment at 40 eV projectile energy. However, comparison between the theoretical and experiment values show only similar qualitative features at 50 and 100 eV. Spin polarization can be a sensitive accuracy indicator for the benchmark of both theory and experiment. Figure 3 [9] presents the integrated relativistic cross section, where the left panel shows the elastic and the right panel shows total scattering cross section for electrons from argon. The trend shows cross section rising from very low energy to form a peak around 15 eV and then decaying. The predicted cross sections from relativistic Dirac equation [9] and from non-relativistic Schrödinger equation [1], both using the optical potential, compare very well with the measured values in the low energy region for both elastic and total scattering. The predicted integrated elastic scattering cross sections using V R do not show any difference between relativistic and non-relativistic approach. They fall almost on the same curve (solid and dot-dashed) and agree well with the measured values from 6 sets of data from different sources-open triangle [22], cross [20], diamond [23], plus [24], open circle [34], asterisk [35]. Cross sections with inclusion of absorption potential [189] (dotted curve) improves the agreement slightly. However, cross sections with inclusion of an older form absorption potential [188] (dashed curve) show better agreement with couple of experiments, for example, open squares [18]. The earlier version of the absorption potential [188] seems to produce total cross sections closer to the experiments than those by version 3 where the measured values are from [25] (open square), [37] (open circle), [38] (cross), [36] (asterisk), [39] (open triangle). , and 300 eV illustrating features from interference effects and comparing theoretical predictions from optical potential approach with experimental measurements. Theory: solid curve represents use of only real part of the the potential and dot that including both real and absorption potentials [9], dash represents nonrelativistic approach with only real potential [1]). Experimental data corresponds to 7 different measurements: at 20 eV, filled circle [22], open circle [18], at 100 eV, open circle [18], solid circle [20], at 300 eV, filled circle [23], open circle [24], asterisks [19]. Right: Angular dependence of spin polarization P for elastic scattering of electrons from Ar at various electron impact energies. Theoretical prediction from optical potential approach [9] is compared with measured values from 3 set-ups-open circles at 40 eV [40], solid circles at 40 and 50 eV [41], cross at 40 eV [42]. Number of measurements for elastic and total scattering of positrons from atoms are relatively lower than that of electrons. Differential cross sections for positron scattering from argon was first studied experimentally by Hyder et al. [27]. The theoretical study for DCS by Nahar and Wadehra [9] using the the optical potential approach described in the theory section was compared (Figure 4) with this single measurement available at that time of reporting. To compare the impact of various absorption potentials, the same set of experimental values was normalized at 90 • with different symbols to differentiate the type of curves produced by different absorption potentials. The comparison shows that the relative features of the measured DCS values agree very well with the predicted values when only the real potential (solid curve) or the complex potential that included version 3 of absorption potential given by Staszewska et al. [189] (dotted curve) were used. It may be noted that in contrast to electron scattering, positron scattering shows less interference structures and decays smoothly with impact energy. Use of other absorption potentials introduce dips and humps in the DCS curves not seen in the experiment.   [189], measured integrated cross sections are seen to be higher than the predicted values (dotted curve). Use of their earlier version of the absorption potential [188] yielded cross sections (solid curve) that show better agreement in the lower energy, but remains higher overall at higher energies. These indicate that absorption potentials for positron scattering have partial success in representing the positron-atom interaction. The right panel presents a later study of Jones et al. [26] who carried out both the experiment and calculated the total cross sections for positrons scattering from argon and neon. They claim significant improvement in agreement between the two. They measured the grant total cross sections (σ GT ) which includes excitations and formation of positronium from zero to 60 eV in contrast to earlier existing experiments, where the energy goes up to 300 eV. They do not show the broad rising peak at low energy as seen Figure 5. Their numbers for the cross sections appear to be different from those in the left panel. The reason for this difference is not clear. For the theoretical predictions, they consider two approaches (i) relativistic optical potential (ROP) method, which reduces to their relativistic polarized orbital method below the first excitation threshold and (ii) convergent close coupling (CCC) method. Their ab initio absorption potential correspond to inelastic scattering due to excitations and positronium formation channels. There is an overall improvement of their approaches in agreeing with the measured values. It is difficult to determine which approach worked better, since their CCC shows better agreement at lower energies with their own experiment, while ROP is better with other measurements.  [21], open circles [25], asterisks [28], and cross [29]. The solid curve corresponds to use of absorption potential of Reference [188] and dotted curve to version 3 potential of Reference [189]. Although sold curve shows agreement in the low energy cross sections the rest of the measured integrated cross sections remain in-between the two predictions. Right: Total integrated cross sections for scattering of positrons from Ar at various energies by Jones et al. [26]: (a) Measured grand total cross section (σ GT ) (filled circles) below the ionization threshold for positronium formation (E Ps are compared with various theoretical approximations, (b) Measured σ GT is compared with other experiments, (c) Measured σ GT above (E Ps are compared with other measured values and theories, (d) Measured cross section σ GT−Ps where cross section for positronium has been subtracted above (E Ps are compared with their ROP prediction.) In contrast to Ar, optical potential approach has shown better agreement with experiments in reproducing the integrated elastic scattering cross sections for Mg (e.g., Reference [10]) and total scattering cross sections for Na atoms (e.g., Reference [12]) by positrons. The left panel of Figure 6 shows integrated cross sections for elastic scattering of positrons from magnesium by Hossain et al. [10] (solid curve) using an optical potential similar to that used for Ar by Nahar and Wadehra [9] and obtained very good agreement with the measured values of Stein et al. [32,33]). The right panel in the figure presents total cross sections for positron scattering from Na atom [12] where optical potential approach [12] was found to be higher in magnitude than other results, but with similar shape and within the experimental errors of Kaupilla et al. [30] and Kwan et al. [31]. The cross-sections of Reid and Wadehra [13], Hewitt et al. [14] and Lugovskoy et al. [15] falls below the SCOP data. However, the data of Cheng et al. [16] show a different nature, even though falls very close to other values.

e ± Scattering from Molecules
We present some illustrative examples of electron and positron scatterings from molecules. Figure 7 presents positron impact total cross section from HCl molecule. The result obtained employing the SCOP method is compared with the theoretical values of Sun et al. [213] and measurements of Hamada et al. [214]. The theories agrees reasonably well in the high energy range. However, below 40 eV the data of Sun et al. [213] underestimate SCOP. This is probably because they performed the calculation using additivity rule to compute the total cross sections for the molecule. Further, they have taken ∆ p as the inelastic threshold, which results in larger magnitude of cross sections at low energies. On the other hand, SCOP values show excellent agreement the measurements of Hamada et al. [214], both in shape and magnitude. Their values are corrected for forward angle scattering, which strengthen our case in terms of reliability of data obtained using SCOP method. Figure 8 shows Q tot for electron/positron scattering by C 4 H 4 N 2 molecule. For both cases SCOP method is in very good agreement with the theoretical values reported by other groups [93,94,215,216]. However, for positron scattering data of Sanz et al. [93] falls faster that other cross sections. The measurements [217][218][219] also show reasonable agreements with theories, except with Reference [220] for electron scattering and Reference [94] for positron scattering. Both these measurements underestimates other values significantly [221]. Figure 9 shows the comparison of electron/positron scattering total and ionization cross sections respectively for propene. In Figure 10 the ratios of cross section is plotted to clearly show the variation in cross section. As expected, the total cross section for electrons is larger than positrons in the intermediate energies. This is very clear in the ratio in Figure 10. However, for the ionization curve, the trend is reversed at intermediate energies. At high energies, typically around 1 keV, both curves tends to converge. This is because at such energies the scattering potential is weak compared to the kinetic energy of the projectile and the scattered wave is approximated as plane wave, in accordance with the first Born approximation.   [222]. Right: Total cross section for e + -pyrimidine molecule interaction [222]. Figure 9 shows the comparison of electron/positron scattering total and ionization cross sections respectively for propene. In Figure 10 the ratios of cross section is plotted to clearly show the variation in cross section. As expected, the total cross section for electrons is larger than positrons in the intermediate energies. This is very clear in the ratio in Figure 10. However, for the ionization curve, the trend is reversed at intermediate energies. At high energies, typically around 1 keV, both curves tends to converge. This is because at such energies the scattering potential is weak compared to the kinetic energy of the projectile and the scattered wave is approximated as plane wave, in accordance with the first Born approximation.  [223]. Right: Total ionization cross section for e − /e + -propene interaction [223]. Positron impact total cross sections are computed for many atomic and molecular targets by Antony and his group using SCOP formalism. The list of the targets studied, with references to articles published, are presented in the Table 2. Table 2. Positron scattering cross section predicted through optical potential method with references (Q tot : Total CS, Q ps : positronium formation CS, Q el : elastic CS, Q ion : direct ionization CS, Q mtcs : momentum transfer CS, Q inel : ine;elastic CS, Q tion : total ionization CS, Q di f f : differential CS).

Target Cross Sections Reference
Inert gases Q tot , Q ion , Q ps , Q tion [11]

Conclusions
We present a review of the elastic and inelastic scattering of electrons and positrons from atoms and molecules. The present conclusion can be given by the following points: • Current theoretical methods for electron scattering from inert atoms are well established. However, their predictions for spin polarization may provide a measure of their effective representation. • Current theoretical method for elastic scattering of positrons from inert gases, single and double electrons systems and molecules have also been developed well, but further studies are needed.

•
There is a critical need for improved absorption potential for positron scattering from atoms and molecules.

•
Further study is needed to expand the scope for scattering from ions that are abundant in astrophysical plasmas. Funding: This research received no external funding.