Calculations of positron scattering from atomic oxygen

The single-center convergent close coupling (CCC) method has been applied to positron scattering from the oxygen atom. Cross sections have been calculated for total, elastic, momentum transfer, total ionization, inelastic, and excitation scattering processes from threshold to 5000 eV. Due to their relevance in transport studies, we also present stopping power and mean excitation energy results. Low-energy studies have been conducted to calculate scattering length, the hidden Ramsauer–Townsend minimum, and the energy of the positron virtual state. CCC-scaled complex optical potential calculations are utilized to calculate positronium-formation, direct ionization, and values between the positronium-formation and ionization thresholds. Good agreement is generally observed at high energies with past theory, electron-atomic oxygen experiment, and halved electron/positron O2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document} experiment. Large discrepancies, however, have been found between current and previous calculations for positronium-formation cross section and for low-energy elastic scattering cross section.


Introduction
The oxygen atom (O) is a fundamental part of life as it composes molecular oxygen (O 2 ), water (H 2 O), and many important biomolecules, such as the bases of DNA.Consequently, atomic oxygen is the main component of the human body, comprising approximately 65% of its mass.There is great interest in positron interactions within the human body due to the application of positrons in the biomedical industry, most notably in positron emission tomography (PET) scans [1] and positherapy [2].Quantifiable data of the scattering processes occurring in these procedures are currently an active area of research [3].To model positron collisions with the complex molecules present in the body, theoretical methods rely upon accurate cross sections of constituent atoms [4][5][6][7].Therefore, accurate calculations of positron scattering from the oxygen atom are a necessary step required to facilitate such modeling.
Atomic and molecular oxygen are a fundamental part of the Earth's atmosphere.Due to the high chemical reactivity of atomic oxygen, it does not exist in substantial quantities close to the Earth's surface.Instead, it is predominantly found in the upper parts of the atmosphere, where abundant ultraviolet radiation continuously dissociates O 2 .Accordingly, studies have found atomic oxygen to be the main component of the atmosphere between 200 and 650 km above the Earth's surface [8,9], with atmospheric models predicting that, by a e-mail: nicolasamori@protonmail.com (corresponding author) volume, 96% of the atmosphere in low earth orbit consists of atomic oxygen [10,11].As a result, atomic oxygen is essential for cooling the atmosphere through collisions with other atmospheric molecules or direct radiative cooling to the ground state [9].
In the same region where atomic oxygen is dominant, positrons are readily produced through cosmic ray interactions with the atmosphere, terrestrial gamma ray flashes, or via nuclear reactions [12][13][14][15].These positrons are either destroyed within the atmosphere, escape to space, or can become trapped and accumulated by the inner magnetosphere into a positron radiation belt [13][14][15][16].Therefore, positron collisions with atomic oxygen are likely most common in the upper atmosphere.Accurate cross sections will, therefore, be essential for modeling positron transport through this region, or in the atmospheres of oxygen-rich bodies such as Venus [17] and Europa [18].
For positron scattering on atomic oxygen, several theoretical methods have been used to obtain a variety of important cross sections for this system.Total cross sections have been calculated by Singh et al. [19] using a spherical complex optical potential (SCOP), Reid and Wadehra [20] using a parameter-free model potential, and by Pindariya et al. [21] using a quantum mechanical formulation based upon a complex atomic spherical potential.Singh and Antony [22] and Pindariya et al. [21] have also obtained results for the total ionization cross section (TICS) through utilization of the complex scattering potential-ionization contribution (CSPic) technique [23].Pindariya et al. [21] obtained results for this cross section only above the ionization threshold, but their approach could not distinguish between positronium-formation and direct ionization.Singh et al. [22], on the other hand, were able to calculate results for positronium-formation and direct ionization separately.
For the elastic cross section, calculations have been completed by Reid and Wadhera [20] and Pindariya et al. [21], using the same methods as before, and also by Dapor and Miotello [24] for incident positron energies between 500 and 4000 eV.Dapor and Miotello [24] have also calculated elastic differential cross sections (DCS) for each of these energies.Total inelastic and the total bound excitation cross sections have been presented solely by Pindariya et al. [21].Stopping power results have been presented for this system by Gumus et al. [25] from calculations by a generalized oscillator strength (GOS) model and the penelope program [26].No experimental studies have been conducted for positron scattering from atomic oxygen.
As cross sections for incident positron and electrons are expected to become equivalent at higher energies, direct comparisons can be made for energies above approximately 500 eV.Electron scattering from atomic oxygen has a number of experimental and theoretical results; however, there has been no recent experimental research for the e − −O system, even after numerous requests from researchers [27][28][29][30].Measurements have been taken for total cross sections below the ionization threshold [31][32][33][34], ionization from threshold to 2000 eV [35][36][37], and multiple excitation and optical emission cross sections which are described in several review articles [29,30,38].This system is well documented theoretically, with calculations existing for excitations, optical emission, stopping power, total ionization, and total cross sections from multiple different approaches [27,28,[39][40][41][42][43][44][45].
We have utilized the single-center CCC method, previously applied to positron scattering from carbon [46], to calculate e + -O scattering.To facilitate calculations of positronium-formation, direct ionization, and cross sections between the positroniumformation and ionization thresholds we utilize a complex model potential calculation scaled to our singlecenter CCC results.We have addressed discrepancies found for low-energy positronium-formation results for positron scattering on atomic hydrogen [46] within this approach by including a scaled component that is expected to increase the accuracy at energies below the positronium-formation maximum.Therefore, we have obtained total, direct and total ionization, positroniumformation, elastic, total bound-state excitation, stopping power, and elastic differential cross sections for incident energies between threshold and 5000 eV.Quantities of mean excitation energy, scattering length, and the energy of the positron-oxygen virtual state are also obtained.

Atomic structure calculation
We utilize a configuration interaction (CI) representation of the atomic wave function where φ i are antisymmetrized N e -electron functions, C i are CI coefficients, N is the number of configurations, and x i = (r i , σ i ) represents the spatial (r i ) and spin (σ) coordinates of electron i.
To obtain these CI coefficients the target Hamiltonian (H) is diagonalized with a finite N number of square integrable basis functions.We obtain N target pseudostates (Φ N n ) which, for large N , form a nearcomplete basis and satisfy where N n is the pseudostate energy and δ n,m is the Kronecker delta symbol.

Scattering calculation
The total Hamiltonian of the positron scattering from the target atom, where E is the energy of the collision system, V is the positron-atom interaction potential, and Ψ (+) is the total scattering wave function.
To solve equation (4) the single-center CCC method relies upon expanding the total scattering wave functions with the previously obtained target pseudostates where n (r 0 ) is a positron channel function.This expansion is substituted into the Schrödinger equation which can then be transformed through the Green's function approach into the coupled Lippmann-Schwinger equations for the T matrix Here |k n is a plane wave with energy k 2 n /2 and V is the interaction potential between the positron and the target.
We solve coupled Lippmann-Schwinger Eq. ( 6) by expanding in J partial waves of total orbital angular momentum and solving each partial wave separately.Analytical Born completion is also used to increase the efficiency of calculations, as the Born approximation is applicable for high partial waves [55].When this equation is solved, we obtain the T -matrix elements for each of the pseudostates from which relevant cross sections can be extracted.Mass stopping power and mean excitation energy can be calculated following the procedure described in [53,54].

Calculation details
Orbitals 1s to 2p were obtained from a HF calculation of O + , optimized on its ground state.All other orbitals were obtained from a Laguerre basis, where L 2 +1 k−1 are the associated Laguerre polynomials, α are exponential fall-off parameters, and N is the number of functions for each .For this calculation, N = 18 − with α = 0.4 for 0 ≤ ≤ 4 and α = 1.0 for ≤ 8.The present close-coupling expansion includes 2s 2 np 4 , 2snp 5 , and np 6 .The 2s 2 2p 3 n continuum is included for ≤ 8.We also include the 2s2p 4 n continuum, but with ≤ 4, and 2s2p 3 n 2 configurations for orbitals n between 3s and 5s.To limit the size of the calculation the only 2s2p 3 n n configurations included were 2s2p 3 4 s4p, 2s2p 3 4s4d, 2s2p 3 4d4f , 2s2p 3 4p4d, and 2s2p 3 4d5s.The restriction of these configurations is expected to slightly decrease the accuracy of the current structure model.States were included with energies up to 75 eV above ionization, resulting in a total of 1543 states.
For energies above 500 eV, calculations were completed with max = 4, which is sufficient to obtain convergent results at these energies.Below the positroniumformation threshold, a 1409 state model with α = 5.0 and N =25 − containing all states with energies up to 300 eV was used.For energies above 100 eV, TICS and stopping power were extrapolated using a Born model with the same configuration structure as the Fig. 1 Convergence study for the total ionization cross section of O where max represents the maximum in which calculations are completed CCC model but with N = 25 − , the 2s2p 4 n continuum extended to ≤ 8, and all 3863 generated states included.Calculations are completed to J = 10 for all energies, except for presented DCS results which were obtained with J ≤ 50 and for low energies where convergence can be achieved for J ≤ 5.
A convergence study for the TICS is presented in Fig. 1.Calculations are fully converged for max = 7 for energies ≥ 20 eV, max = 6 for energies ≥ 100 eV, and max = 4 for energies ≥ 300 eV.

CCC-scaled complex model potential
A full description of the complex potential we utilize for these calculations is provided in [46].Therefore, only a summary is provided here.We calculate the optical potential (V opt ) through, where V st is the static potential, V pol the polarization potential, and V abs the absorption potential.For V abs , we utilize the absorption potential of Staszeweska et al. [56], Here, Δ is the absorption threshold, k F (r) is the magnitude of the Fermi wave vector, T loc is the local kinetic energy of the positron, and θ(x) is the Heaviside step function.The dynamic functions A 1 , A 2 , and A 3 are provided by [56] and are unmodified in the current formulation.
182 Page 4 of 14 Eur. Phys.J. D (2023) 77 :182 As positronium-formation cannot be explicitly included in the absorption potential, we rely upon the delta variational technique to calculate positronium-formation cross sections.This entails modifying the absorption threshold with the form given by Chiari et al. [57], here Δ e and Δ p are the electronic-excitation and positronium-formation threshold energies, respectively.For the adjustable parameter E m , we have used the energy at which the single-center CCC TCS has a maximum.
To solve the scattering equations for the complex model potential, we solve the Lippmann-Schwinger equations for a potential scattering system to obtain the T -matrix elements.These on-shell partial-wave Tmatrix elements are then utilized to obtain the TCS, elastic, and inelastic cross sections.The positroniumformation cross section (σ Ps ) is given by the difference between the inelastic cross sections obtained from calculations with Δ(E) = Δ e and Δ(E) given by Eq. ( 10).For energies below the maximum of the positroniumformation, we have modified this approach as described in Sect.2.5.
We then scale the cross section for inelastic scattering minus positronium-formation to reproduce the total cross section σ tot obtained from the single-center CCC calculations at high energies.Following this, we scale the positronium-formation cross section σ Ps to reproduce the single-center σ tot at the maximum cross section between the ionization threshold and 10 eV above this threshold.The direct ionization is obtained by subtracting the calculated positronium-formation values from our CCC TICS results.To calculate direct ionization for energies below 10 eV above the ionization threshold, we have followed the CSP-ic method [23,46].Excitation cross sections are then found by subtracting the direct ionization cross section from the total direct inelastic cross section.Calculations that have utilized this approach to obtain a range, or all, of their values are referred to as CCC-pot.

Low-energy positronium-formation scaling
Through analysis of calculations for atomic hydrogen, available in appendix of [46], we found that the CCC-pot approach significantly underestimates twocenter calculations for positronium-formation at energies below its maximum.As atomic oxygen has a positronium-formation threshold almost equal to that of atomic hydrogen, we expect this same deficit to occur for this calculation.To accommodate for this, we have extended the CCC-pot calculation to include a scaling factor which increases the positronium-formation component of V abs for energies between the positroniumformation threshold and the maximum positroniumformation cross section calculated in our CCC-pot method.The values for this factor were determined to produce positronium-formation equal to that of the  two-center calculation for atomic hydrogen for energies below the positronium-formation maximum.The unscaled and scaled results for atomic hydrogen are shown in Fig. 2 alongside the two-center CCC calculation [58].
To allow for calculation at any energy below the maximum of positronium-formation, this factor was then fitted with a function of the form where Here, E is the incident energy and E H Ps is the positroniumformation threshold of atomic hydrogen, with both units in eV.The parameters for Eq.11 are shown in Table 1.
To apply this to atoms beyond atomic hydrogen, we have formulated an approach based upon the difference in energy between the positronium-formation threshold and maximum for each atom.This replaces Eq. 12 with, where Here, E A Ps is the positronium-formation threshold for atom A and E A max refers to the energy that the maximum positronium-formation is found for this atom, with E H max representing this for the hydrogen atom.For Eq. 13 units are in eV, whereas Eq. 14 is dimensionless.It is expected that this approach will decrease errors at lower incident energies and better approximate twocenter CCC calculations.

O structure
Excitation energies for triplet target states of the atomic oxygen atom are presented in Table 2. Agreement between the current calculation and past theory and experiment is satisfactory, with differences between NIST [59] and CCC ranging from 0.002 to 0.278 eV.Oscillator strengths for several transitions between triplet bound states are shown in Table 3. Good agreement is found with NIST and past theory for all presented transitions except for the 2p 3 3s 3 D 0 transition, though this is a weak transition with a small oscillator strength.This discrepancy is expected to mainly result from the restrictions imposed on the included configurations.An important quantity for scattering, mostly at low incident energies, is the dipole polarizability (α D ).For O, this is predicted by experiment to be 5.2±0.4 a 3 0 [60].Our largest 3863 state Born model had α D = 5.50 a 3 0 , which is within experimental error and slightly higher than other theory [40].Our 1543 state model had a dipole polarizability of 6.00 a 3 0 and our small-energy 1490 state model of 5.31 a 3 0 , which is within 3% of the largest model and within experimental uncertainty.

Total cross section
The total cross section arising in the CCC calculation of atomic oxygen is presented in Fig. 3 alongside the theoretical calculations of Singh et al. [19], Reid and Wadehra [20], and Pindariya et al. [21].In this case, the CCC-pot is equivalent to the 1490 state CCC calcu-182 Page 6 of 14 Eur.Phys.J. D (2023) 77 :182 lation below the positronium-formation threshold and the 1543 state model above 23.5 eV.The calculation is equivalent to the CCC-scaled complex model potential between these two energies.Due to the lack of existing experimental data for electron and positron scattering from atomic oxygen, we compare against experimental results for positron and electron results for O 2 , which are divided by 2, in accordance with the additivity rule.This rule is expected only to be relevant for higher energies.Therefore, only qualitative comparisons can be made for lower energies.
CCC results are higher than the halved positron O 2 experimental results of Charlton et al. [62] and Dababneh et al. [63] for energies above 40 eV.However, comparison with the more recent measurements of Chiari et al. [57], shows near-perfect agreement.For high incident energies, positron and electron results are expected to become equal to the Born approximation, which is independent of projectile charge.We observe this for energies above 1000 eV, with excellent agreement found with the halved electron O 2 measurements of García et al. [64].
The peak cross section of the CCC TCS is of similar magnitude to the calculation of Pindariya et al. [21] and the halved positron O 2 experiment of [57].However, the results of Singh et al. [19] are much larger than our CCC result at the peak TCS and also depicts significantly different behavior for lower energies.At low energies, our results behave similarly to O 2 experiment, with results flat from 1 eV to the positronium-formation threshold.On the other hand, Singh et al. [19] predict the cross section to decrease significantly from 1 eV to the positronium threshold and the experiment of Dababneh et al. [63] predicts the cross section to increase over this energy range.Consequently, this calculation predicts a value at 1 eV over five times higher than the CCC result.We expect that CCC results converge to exact results at energies below the positroniumformation threshold, given in Fig. 3. Therefore, the discrepancy is likely due to the polarization potential being too high in the model potential or from the exclusion of virtual positronium-formation in the SCOP approach.A similar discrepancy was noted between the CCC and results of Singh et al. [19] for the carbon atom [46], but not as pronounced.

Total ionization cross section
The total ionization (electron-loss) cross section (TICS) for positron scattering on oxygen, which is equivalent to the sum of direct ionization and positronium-formation, is shown in Fig. 4. As with the total cross section, for energies below 23.5 eV the CCC-pot results are from the CCC-scaled complex model potential, whereas, for energies above this limit they are from the single-center CCC calculation.The results of [21] are in close perfect agreement with those of CCC from near the direct ionization threshold through to high energies.The results of [22] are considerably higher for energies above 10 eV.  4 Total ionization (electron-loss) cross section for positron scattering on oxygen.Current CCC results are presented alongside the calculations of [19] and [21].The CCCpot calculations incorporate an estimate of the positroniumformation cross section contribution from its threshold, see text Fig. 5 Direct ionization cross section for positron scattering on oxygen.Theoretical results include CCC-pot and Singh and Antony [22] calculations.Experimental measurements for electron scattering from Zipf [36], Brook et al. [35], and Thompson et al. [37] are expected to converge to positron scattering ones at high energies

Direct ionization and positronium-formation cross section
The e + −O direct ionization cross section is presented in Fig. 5.This includes the current CCC-pot results alongside the theoretical calculation of Singh and Antony [22] and e − − O experiment.At high energies, electron and positron results are expected to become equal.We find this to occur by 1000 eV, with excellent agreement between the CCC and electron experiment above this energy.Similar to the TICS, the calculation of Singh and Antony [22]  are presented in Fig. 6.Current CCC-pot results are presented together with the only other theoretical calculation, by Singh and Antony [22].The experimental results are for the O 2 , and have therefore been divided by two.The positronium-formation threshold energy of O 2 is lower than that of O, which is why near-threshold energies, there is a discrepancy in the behavior of the O 2 experiment and theoretical values for O. Compared to the halved O 2 experiment, we find agreement with the measurements of Marler and Surko [65] for energies above 40 eV, and Griffith [67] above 100 eV.Apart from the lower energies the results of Singh and Antony [22] are generally considerably above those of CCCpot.A similar disagreement was found comparing these approaches for atomic carbon [46].The absence of existing experiment for atomic oxygen makes it difficult to judge the true accuracy of the CCC-pot calculations, particularly at the lower energies.However, as our calculation is scaled directly to the ab-initio CCC calculation and is in agreement with halved O 2 experiment at the higher energies, we believe lower energy CCC-pot cross section uncertainty is no more than 20%.

Inelastic cross section
Figure 7 presents the total inelastic cross section for positron scattering on oxygen from its threshold to 5000 eV.Above 500 eV, close agreement is found between the CCC and the calculation of Pindariya et al. [21].The CCC is slightly above the results of Pindariya et al. [21] between 100 and 500 eV.Below this, Pindariya et al. [21] are notably larger than the CCC results, except for below 20 eV, where they are smaller.The calculation of Reid and Wadehra [20] overestimates other theory over their calculated energy range of 100-5000 eV.

Elastic and momentum transfer cross section
The integrated elastic cross section is shown in Fig. 8 alongside the theoretical calculations of Reid and Wadehra [20], Pindariya et al. [21], Dapor and Miotello [24], and, for the electron case, NIST.Above 500 eV, we find near-perfect agreement with the calculations of Dapor and Miotello [24] and above 3000 eV with the NIST electron results.Comparisons are also made with the halved electron experimental results of Dapor and Miotello [68] for O 2 .At 1000 eV, the electron and positron theory lies within the experimental uncertainty.
In Fig. 9 we present the momentum transfer cross section for this system.Current results are compared with the calculations of Dapor and Miotello [24], the electron calculations of NIST, and halved electron-O 2 experimental data of Iga et al. [68].As with the elastic cross section, for energies above 500 eV excellent agreement is observed between the CCC and the calculations of Dapor and Miotello [24].The NIST results for the electron case are larger than the positron theoretical results, with no agreement viewed by 5000 eV.The experimental data are found to lie between the positron and electron calculations.
The DCS for various energies between 500 and 3500 eV are shown in Fig. 10.As with the integrated cross section, excellent agreement is viewed with the calculations of Dapor and Miotello [24] for each of these energies.In Fig. 11 we present a selection of elastic DCS for incident energies between 1 and 100 eV.A technique similar to Green et al. [69] is used to ensure convergence in the partial wave expansion, resulting in a cusp at low scattering angles for incident energies above 1 eV.As with carbon [46], this feature is not present at lower energies due to the triplet ground state resulting in pwave scattering being dominant at low energies.[20], Pindariya et al. [21], and Dapor and Miotello [24].Experimental results for positrons for electrons are from Iga et al. [68].The vertical dotted lines represent the positronium-formation and ionization thresholds  9 Momentum transfer cross section for positron scattering on oxygen.CCC results are presented alongside the positron calculation of [24] and the NIST electron results [59].Halved O2 experimental results for electrons are from Iga et al. [68].The vertical dotted lines represent the positronium-formation and ionization thresholds

Low-energy study
Elastic and momentum transfer cross sections for energies below 1 eV are shown in Fig. 12.Although the IAM approximation does not hold for low energies, we find good agreement with the halved O 2 experiment from Chiari et al. [57]; this experiment does not have forward-scattering corrections so likely underestimates the actual result for O 2 .However, similar behavior exhibited provides some evidence of the accuracy of the current results.At lower energies, as scattering becomes isotropic, the elastic and momentum transfer cross sections are expected to become equal.We observe this for energies below 10 −3 eV.
where A is the scattering length, and the asymptotic value of our elastic cross section at 10 −5 eV, we determine A = −1.862for this system, where the sign is derived from our phase shift at low energies [46].As |A| is larger than the mean radius of O, the scattering cross section is greater than the geometric size of the atom at low energies.This results from the elastic cross section being enhanced due to a virtual level of the positron projectile [70].The energy of this can be calculated through which for this system is = 3.92 eV.
As the scattering length is negative this implies a Ramsauer-Townsend minimum that occurs at incident Using this equation, we obtain an energy of E min = 1.53 eV, which (due to neglected higher order terms in Eq. 17) is 0.5 eV lower than the value of 2.05 eV observed in the s-wave scattering component of the elastic cross section in Fig. 13.In this figure, the s-, p-, d-, and higher-wave cross sections refer to the angular momentum of the incident positron L for which their contributions to the elastic cross section are extracted from the relevant partial-waves J.The Ramsauer-Townsend minimum is not observed in the total elastic cross section as the magnitude of the s-wave component becomes small relative to other partial waves as the incident energy increases.This process can be observed in Fig. 13, with the peak of the p-wave and increasing values of the d-and higher-wave cross sections supplementing the decreased s-wave contribution at its minimum.As we observed in carbon [46], although the effect of the Ramsauer-Townsend minimum is hidden for the elastic cross section, its impact can be observed for the momentum transfer cross section, which has a minimum at 1 eV due to this effect.

Bound-state excitation cross section
There is no other existing positron theory or experiment for specific bound-state excitations.Therefore, we compare against existing electron theory and experiment, which are expected to be equivalent to positron results for high incident energies.To obtain stable results for energies below 20 eV, values were obtained for calculations with max = 2 and only bound states included.Figure 14 presents the excitation of the ground state to 2p 3 3s 3 S o .Excellent agreement is found for energies above 500 eV with the recommended electron theoreti- Fig. 14 The 2p 4 3 P → 2p 3 3s 3 S o excitation cross section for positron scattering on oxygen.Theoretical results for incident electrons are from Tayal and Zatsarinny [39] and the recommended values from Johnson et al. [72].Experimental electron measurements are from Johnson et al. [61], Kanik et al. [73], and Vaughan and Doering [74] cal results of Johnson et al. [72].At lower energies, our positron calculation lies within the experimental error of existing electron measurements and is close to the BSR calculation of [39]. Figure 15 presents the 2p 4 3 P → 2p 3 3 s 3 S o excitation cross section.Comparing with the BSR results, we find that differences decrease with increasing energy and that they become close by 100 eV.The CCC results lie within the uncertainty of the experiment of Gulcicek et al. [75] below 20 eV.
The 2p 3 3d 3 D o excitation cross section is shown in Fig. 16.Due to the differences observed with the oscillator strength for this transition and those of NIST, we have applied OOS scaling, which is the scaling of the CCC cross section by the ratio between these two oscillator strengths.Agreement is observed after this scaling with the recommended electron results of Johnson Fig. 16 The 2p 4 3 P → 2p 3 3d 3 D o excitation cross section for positron scattering on oxygen.Theoretical results for incident electrons are from Tayal and Zatsarinny [39] and recommended from Johnson et al. [72].Experimental electron measurements are from Kanik et al. [73] et al. [72] for energies above 500 eV.Closer agreement is observed with the BSR calculation for the unscaled result of this excitation.
Figure 17 shows the 2p 3 3s 3 D o excitation cross section.After OOS scaling we find excellent agreement with the recommended electron results of Johnson et al. [72] at high energies.At lower energies our positron calculation is found to be lower than the electron theory, but within the uncertainty of electron experiment.
The total bound-state excitation cross section is presented in Fig. 18.The CCC-pot calculation is equivalent to the CCC-scaled complex model potential calculation for energies below 23.5 eV and the CCC for energies above this.Little agreement is observed with the only existing positron calculation of Pindariya et al. [21], which is significantly larger than the CCC result for  18 Total bound-state excitation cross section for positron scattering on oxygen.Theoretical results are from the CCC, Pindariya et al. [21], and, for the electron case, Joshipura and Patel [41] and Johnson et al. [72] energies below 300 eV.However, both calculations find a maximum at ≈ 20 eV.For energies above 500 eV, the electron results of Joshipura and Patel [41] agree with the positron results of Pindariya et al. [21] above 500 eV, but both are lower than the CCC results.The presented calculation of Johnson et al. [72] is a sum of only the three previously presented excitations.Therefore, it is expected to underestimate the total as it does not contain enough states to approximate the full excitation spectrum.Consequently, the summed recommended values follow the shape of the CCC calculation, which contains 28 bound states, for energies above 100 eV but is much lower even at 1000 eV.As the summed Johnson et al. [72] results are higher than the other presented electron and positron calculation at Theoretical calculations for this system are from the current CCC-pot calculation, the penelope code [26], and Gumus et al. [25].Calculations for e − -O by Gupta et al. [45] and semi-empirical results for the e − -O2 system by Williart et al. [77] are also presented.The positronium-formation component of the stopping power (PosF) is also shown for the CCC-pot calculation this energy, this suggests they underestimate the total bound-state excitation cross section at higher energies.

Mean excitation energy and stopping power
In Fig. 19, we present the stopping power for a positron incident on atomic oxygen.The calculation of stopping power requires the cross section and excitation energy for each reaction channel.The stopping power for the positronium-formation channel was obtained from the total positronium-formation cross sections, as discussed in Ref. [53].Reaction channels corresponding to boundstate excitation and ionization are treated in the following way.Below 23.5 eV, for bound-state excitation, we have used the present estimates for the cross sections.
For ionization, we have used cross sections for positiveenergy pseudostates from the single-center CCC model scaled to obtain the present direct ionization cross section.For energies above 23.5 eV, the stopping power is obtained directly from the single-center CCC calculation.
For energies above 200 eV, CCC, the GOS model of Gumus et al. [25] and the penelope code [26] predicts similar behavior, with the CCC results slightly below these other models.Below 200 eV, the CCC-pot calculation is in significant disagreement with these other approaches and predicts a substantially lower stopping power.As the models used for these positron calculations apply simple approximations, do not account for processes such as positronium-formation, and do not directly model the atom's target structure, the current CCC calculation is expected to be more accurate for the presented energies.
According to Bragg's additivity rule [76], which is accurate for high energies, the stopping power of O 2 is expected to be equal to that of O. CCC calculations of stopping power for H and H 2 [53] targets found this rule accurate for energies above 100 eV.As electron and positron results are equivalent for high energies, we also present the e − -O calculation of Gupta et al. [45] and the e − -O 2 semi-empirical calculations of Williart et al. [77].For energies above 1000 eV, we find excellent agreement between the current CCC results and those of Williart et al. [77].The results of Gupta et al. [45] are in close agreement with those of the CCC and Williart et al. [77] for energies above 3000 eV.The mean excitation energy, directly obtained from the stopping power and the inelastic cross section, is shown in Fig. 20.There are no existing calculations of mean excitation energies for positron scattering on atomic oxygen or O 2 .Williart et al. [77] only provide experimental measurements and their uncertainty for incident energies of 300 and 600 eV for the e − -O 2 system.As their stopping power result was obtained with their experimental measurements of mean excitation energy, we have extracted the experimental values for other energies using their presented stopping power and inelastic cross section.There is no discussion provided of the experimental uncertainty for these other incident energies.
The mean excitation energy between atomic oxygen and O 2 is expected to be similar at higher energies.For CCC calculations of e − and e + scattering on H 2 [53], the e + mean excitation energy was found to follow similar behavior as the e − result but at a magnitude ≈ 20% higher.A similar trend is found here but with a smaller difference between the e + and e − result.At 600 eV, the CCC results lies just within the electron experimental uncertainty.Above this energy, as our results decrease, the difference between them and this experiment decreases to be within 10% at 5000 eV.This experiment and the CCC results exhibit similar behavior, with the mean excitation energy increasing to 600 eV.Above this energy, the mean excitation energy becomes nearly constant in the experiment and slowly decreases in the CCC-pot calculation.The sharp rise observed in the CCC-pot at threshold is due to positronium-formation.

Conclusion
A comprehensive set of cross sections and other quantities have been calculated for positron scattering from atomic oxygen within the single-center CCC approach.The target structure model for atomic oxygen adopted in the present calculations is in good agreement with the energy levels and oscillator strengths of NIST and previous theory, reflecting its accuracy.The CCC-scaled complex model potential results was used to obtain cross sections between the positronium-formation and ionization thresholds, positronium-formation, and direct ionization cross sections.A modification is introduced to this approach that is expected to result in more accurate positronium-formation cross sections between its threshold and maximum.
The total, elastic, momentum transfer, direct ionization, positronium-formation, total ionization, inelastic, excitation, and stopping power cross sections were calculated between threshold and 5000 eV.Quantities such as the scattering length, mean excitation energy, hidden Ramsauer-Townsend minimum, and the energy of the virtual positron-oxygen state were also calculated.Good agreement is observed between current results and high-energy electron atomic oxygen experiment and theory for total ionization, total, elastic, and bound-state excitations.Halved O 2 experiments for both positron and electrons are also in good agreement with the present calculations.For past positron theoretical calculations, agreement is observed at high energies with some available approaches.However, large differences are observed at lower energies and for the positronium-formation cross section.Future work, particularly experimental, is recommended on this system to address the discrepancies found.

Fig. 2
Fig.2Impact of scaling factor on CCC-pot positroniumformation cross section compared against two-center CCC results[58]

Fig. 3
Fig.3 Total cross section for positron scattering on oxygen.Theoretical CCC results are shown alongside optical potential calculations by Singh et al.[19], Reid and Wadehra[20], and Pindariya et al.[21].Experimental O2 results for positron scattering are from Charlton et al.[62], Chiari et al.[57], and Dababneh et al.[63].Electron O2 results are by García et al.[64].The vertical dotted lines represent the positronium-formation and ionization thresholds Fig.4 Total ionization (electron-loss) cross section for positron scattering on oxygen.Current CCC results are presented alongside the calculations of[19] and[21].The CCCpot calculations incorporate an estimate of the positroniumformation cross section contribution from its threshold, see text

Fig. 7
Fig.7 Total inelastic cross section for positron scattering on oxygen.CCC results are shown alongside those of Pindariya et al.[20,21]

Fig. 8
Fig.8 Total elastic cross section for positron scattering on oxygen.Theoretical results include the CCC, Reid and Wadehra[20], Pindariya et al.[21], and Dapor and Miotello[24].Experimental results for positrons for electrons are from Iga et al.[68].The vertical dotted lines represent the positronium-formation and ionization thresholds Fig.9 Momentum transfer cross section for positron scattering on oxygen.CCC results are presented alongside the positron calculation of[24] and the NIST electron results[59].Halved O2 experimental results for electrons are from Iga et al.[68].The vertical dotted lines represent the positronium-formation and ionization thresholds

Fig. 10 Fig. 11
Fig.10 High-energy elastic differential cross sections for positron scattering on oxygen.CCC results are presented alongside the theoretical calculations of Dapor and Miotello[24]

Fig. 12
Fig.12 Elastic and momentum transfer cross sections for positron scattering on atomic oxygen at energies below 1 eV.Halved O2 experimental results are from Chiari et al.[57]

Fig. 13
Fig.13 Elastic cross sections for positron scattering on atomic oxygen and its s-, p-, d-, and higher -wave contributions.Here, L refers to the angular momentum of the incident positron

Fig. 17
Fig.17 The 2p4 3 P → 2p 3 3s 3 D o excitation cross section for positron scattering on oxygen.Theoretical results for incident electrons are from Tayal and Zatsarinny[39] and recommended results from Johnson et al.[72].Experimental electron measurements are from Vaughan and Doering[74] and Kanik et al.[73]

Fig. 19
Fig.19 Stopping power for positron scattering on oxygen.Theoretical calculations for this system are from the current CCC-pot calculation, the penelope code[26], and Gumus et al.[25].Calculations for e − -O by Gupta et al.[45] and semi-empirical results for the e − -O2 system by Williart et al.[77] are also presented.The positronium-formation component of the stopping power (PosF) is also shown for the CCC-pot calculation

Fig. 20
Fig.20 Mean excitation energy for positron scattering on oxygen.CCC-pot theoretical results are presented alongside experimental results for the e − -O2 system from Williart et al.[77]

Table 1
Fit parameters

Table 2
Excitation energies (eV) for atomic oxygen triplet target bound states

Table 3
Oscillator strengths for bound triplet states of O