Excitations of the nS States of Atomic Hydrogen by Electron Impact, Excitation Rate Coefﬁcients, and Phase Shifts: Comparison with Positron Impact Excitation

: The excitation cross-sections of the nS states of atomic hydrogen, n = 2 to 6, by electron impact on the ground state of atomic hydrogen were calculated using the variational polarized-orbital method at various incident electron energies in the range 10 to 122 eV. Converged excitation cross-sections were obtained using sixteen partial waves ( L = 0 to 15). Excitation cross-sections to 2 S state, calculated earlier, were calculated at higher energies than before. Results obtained using the hybrid theory (variational polarized orbital method) are compared to those obtained using other approaches such as the Born–Oppenheimer, close-coupling, R-matrix, and complex-exterior scaling methods using only the spherical symmetric wave functions. Phase shifts and elastic cross-sections are given at various energies and angular momenta. Excitation rate coefﬁcients were calculated at various electron temperatures, which are required for plasma diagnostics in solar and astrophysics to infer plasma parameters. Excitation cross-sections are compared with those obtained by positron impact excitation.


Introduction
Cecilia Payne, a student of Eddington, suggested that hydrogen is the most abundant element in the universe.The electron impact excitation of hydrogen to the 2p level resulting in Lyman-alpha radiation in astrophysics is well known.There is a possibility of exciting other levels, such as 2s excitation resulting in two-photon decay.Excitation by positron impact is also possible.It is known that positrons are produced inside the Sun when protons collide with protons, contributing to the opacity of the atmosphere of the Sun, in addition to electrons.Their presence is indicated in the Sun as well as in the interstellar space by the detection of the 0.511 MeV line from the center of the galaxy [1].In earlier publications [2,3], cross-sections for excitation of the atomic hydrogen by electrons and positrons were presented.Cross-sections for excitation are a two-channel problem.They must include the initial state as well as the final state.Burke et al. [4], using the closecoupling approximation, carried out such a calculation.Their wave function includes both 1S and 2S states for 2S excitation.Rudge [5] calculated excitation to the 2S state at 3 incident energies using the Born-Oppenheimer approximation and compared his results with those obtained using the close-coupling approximation [4].Bartlett and Stelbovics [6] used the Temkin-Poet models to study electron-hydrogen scattering and calculated excitation cross-sections to various S states.Jones and Stelbovics [7] calculated excitation to S states numerically using the Temkin-Poet model.Callaway [8] carried out calculations using 11-state expansion including 7 pseudostates.Scott et al. [9] carried out calculations using the R-matrix close-coupling approximation.There are other calculations in which excitation has been treated as a single-channel problem and where the initial wave function is considered to be exact.The earlier calculations [2] are now extended to excitation of the higher S states of the hydrogen atom.Since the hydrogenic functions are exact, different approximations can be used to judge the possibility of obtaining accurate results.The present calculation was carried out using the variational polarized orbital method of Bhatia [10].The asymptotic polarizability obtained in this approximation is 4.5 a 3 0 , the exact value.Even though the calculations were carried out as a single-channel approximation, the phase shifts obtained have lower bounds to the exact phase shifts and agree well with those obtained using other approximations.This method has been applied to calculate parameters of Feshbach resonances [11] and photoabsorption crosssections [12].The results obtained are accurate, compare well with results obtained in other calculations, and agree with the experimental results.It seems justified to assume that the scattering functions required in the calculation of excitation cross-sections are accurate.Calhoun et al. [13], using the distorted-wave approximation, calculated excitation to 2S state at 4 incident energies.Excitation to 2S state was reported earlier [2].It is expected that excitation cross-sections obtained using the present approach will be accurate for nS states, for n greater than 2 as well.In electron impact excitation experiments, the sum of excitation cross-sections and elastic scattering cross-sections might be measured.The present investigation would be useful then.Rydberg units are used throughout (electron mass m e = 0.5, e 2 = 2, h = 1).Distances are measured using Bohr radius = 1, energies are in Rydberg = 13.605eV, and a (−b) represents a × 10 (−b) .

Theory and Calculations
The method for calculating excitation cross-sections is discussed in detail in Ref.
[2] and it is only briefly described here.
The excitation cross section from state 'I' to 'f' state is given by In the above equation, k i and k f are the initial and final momenta and the transition matrix T f i is given by The potential involved in the transition is given by In the above expression, r 1 is the distance of the incident electron from the nucleus, assumed to be infinitely heavy, r 12 is the inter-electron distance, and Z is the nuclear charge.In principle, the initial state wave is supposed to be an exact solution of the Schrödinger equation.It is given by where u( → r 1 ) is the scattering function and 2 ) is the target wave function polarized due to the incident electron which produces a change in the energy equal to − 1 2 αE 2 , where α is the polarizability and E is the electric field produced by the incident electron.Temkin and Lamkin [14] showed that in the presence of the incident electron the target wave function can be written as χ St is a smooth cutoff function introduced by Shertzer and Temkin [15], it is given by The function u 1s→p (r 2 ) in Equation ( 5) is given by The unperturbed target function in the initial state is given by The scattering function in the plane-wave normalization is given by The spherical harmonic, Y L0 , is a function of the angular momentum L depending on spherical polar angles θ 1 and φ 1 .The function u(r 1 ) satisfies the Schrödinger equation: In the above equation, E = E target + k 2 is the total energy, where k 2 ≡ k 2 i is the incident energy.H is the Hamiltonian given by After integrating over d → r 2 and angles θ 1 and φ 1 , the resulting integro-differential equation, which is quite complicated, is solved by a non-iterative method [16].The method being used determine the initial state wave function in the present calculation is the variational polarized orbitals method [10], in which the phase shifts for scattering have lower variational bounds.The final wave function is given by The excited state target functions are φ ns ( → r 2 ), n = 2 to 6.These functions φ ns ( → r ) are given in Appendix A, and k f is the final electron momentum after excitation.

Results
The excitation cross-sections were calculated at the incident momenta k = 0.87 to 3.0 using 16 partial waves (L = 0 to 15).Convergence with the partial waves is shown in Table 2 of reference [2].The 2S excitation cross sections converge by L = 6 or 7, excitations to the higher S states require more partial waves to get converged cross-sections.
Cross-sections for 1S excitation to 2S, having 2 peaks, are shown again in Figure 1.It was indicated in reference [2] that the absolute value of 2S excitation measured by Kauppila [17] is 0.163 ± 0.020 πa 2 0 at 11 ± 0.2 eV, while the present calculation gives 0.137 πa 2 0 , in spite of including contributions of the higher partial waves (up to L max = 15).The excitation cross sections for 2S, 3S, 4S, 5S, and 6S are shown in Figures 1-5.These figures are for L max = 12 by which cross-sections have essentially converged.The excitation cross-sections 1S to 3S state have a peak at k = 1.1, excitation to 4S state at k = 1.02, excitation to 5S state at k = 1.05, and excitation to 6S state at k = 1.05.The excitation cross-sections to 6S state have another peak at k = 1.70.This peak arises because for L = 0, the cross-section is maximum at k = 1.70.This fact has been verified by calculating cross-sections below and above k = 1.70, but close to this k.
The excitation cross sections for 2S, 3S, 4S, 5S, and 6S are shown in Figures 1-5.These figures are for Lmax = 12 by which cross-sections have essentially converged.The excitation cross-sections 1S to 3S state have a peak at k = 1.1, excitation to 4S state at k = 1.02, excitation to 5S state at k = 1.05, and excitation to 6S state at k = 1.05.The excitation cross-sections to 6S state have another peak at k = 1.70.This peak arises because for L = 0, the cross-section is maximum at k = 1.70.This fact has been verified by calculating cross-sections below and above k = 1.70, but close to this k.

Rate Coefficients
In solar and astrophysics, parameters, namely the electron density and temper of the plasma are required.These parameters can be determined by solving the stat equilibrium equations involving electron collision rate coefficients and radiative tion rates.The excitation rate coefficient is defined by  The excitation cross-sections obtained in the hybrid theory [10] are given in Table 1.The empty spaces indicate that the incident energies are not sufficient to excite these states.There was a typo in a routine in the code for calculating excitation cross-sections from 1S state to 2S state, reported in ref. [2].That correction now has been made.The presently calculated excitations to the 2S state differ in the 4th or 5th decimal place compared to those given in reference [2].The excitation cross-sections to the 2S state were compared to the results of the previous calculations (close-coupling, R-matrix, etc.) and those comparisons are not repeated here again.
The excitation cross-sections to the 3S state have been calculated using the closecoupling [18] and the Born approximations by Whelan et al. [19].Their results are higher than the present results.There are measurements of 3S excitations by Mahan et al. [20].Their results include cascades from higher states.Whelan et al. [19] indicate 0.013πa 2 0 for the cross section at k = 2.0, which is lower than the result of Mahan et al. [20], which includes cascades from higher states.By L max = 15, cross-sections have converged to at least to the 4th decimal place.
At k = 2.1, the cross-section of the 6S is 6.76 (−5)πa 2 0 , and k = 2.6, it is 8.33 (−4) for L = 15.Cross-sections start increasing at k = 2.5.The sum of all these cross sections at each energy is also given in Table 1, including elastic cross-sections.The present results should be accurate because the calculations were carried out using the quadruple precision.Callaway et al. [18] used six-state close-coupling expansion including an optical potential to calculate excitation to 2S and 3S states.The 2S excitation results were already compared to those in [18] and 3S excitations are now compared to those in [18] in Table 2. Bartlett et al. [21] carried out measurements at 12.24 eV and compared their results for the excitation of the 3S, 4S, and 5S states with the results of exterior complex scaling calculations, the experimental uncertainty being 35%.Their interpolated results for 3S state are compared with the present results in Table 2. Their results are lower compared to the present results.Similarly, their 4S results are found to be lower than the present results.The converged close-coupling calculations of Bray et al. [22] gave for 2S cross-sections equal to 0.252a 2 o , 0.212a 2 0 , and 0.139a 2 0 at k equal to 1.60, 2.0, and 2.71, respectively.Their results are lower than those given in Table 1 in units of πa 2 0 .Calculations by Poet [23] using spherically symmetric wave functions and R-matrix calculations of Scott et al. [9] do not have excitation to 3S state.We find that the results of the exterior complex scaling calculation [6] are lower while those obtained using the converged close-coupling calculations [22] are close to the present results.These approaches [6,[18][19][20][21][22] are very complex compared to the present variational polarized orbital method, which is very easy to use and provides fairly accurate results.
Phase shifts at various incident energies are given in Table 3a-e for partial waves L = 0 to 15.They are variationally correct, that is, they are below the exact results.The triplet phase shifts are higher than the singlet phase shifts.In Table 3b, for L = 4, at k = 2.6 and k = 3.0, phase shifts appear to be discontinuous.However, if π radians are added, phase shifts are found to be continuous.These phase shifts have been used to calculate the elastic scattering cross-sections, which are also given in Table 1 along with the excitation cross-sections.The sum of the excitation cross-sections and elastic cross-sections would be useful when the measurements of total cross-sections are carried out.

Rate Coefficients
In solar and astrophysics, parameters, namely the electron density and temperature, of the plasma are required.These parameters can be determined by solving the statistical equilibrium equations involving electron collision rate coefficients and radiative transition rates.The excitation rate coefficient is defined by where Ω ij = g i E i σ ij is the collision strength at the incident energy E i , g i = 2 is the weight factor for the initial state, k is Boltzmann's constant, σ ij are the excitation cross-sections given in Table 1, and the lower limit in the integral is the threshold energy for excitation.The excitation rate coefficients for excitation to the nS states, n = 2-6, are given in Table 4.
Defining E ij = E j − E i , it can be shown from the principle of detailed balance that the de-excitation rate coefficients are given by

Comparison with Positron Impact Excitation
As it was mentioned in the introduction that positrons are as important as electrons to describe the opacity of the Sun and interstellar space, it is worth having a comparison of the excitation cross-sections by electron and positron impact on hydrogen atoms.A comparison with the positron impact excitation [3] is given in Table 5. Positron impact cross-sections are higher than the electron impact cross sections.

Conclusions
Calculations were carried for excitation from 1S state to nS states of the hydrogen atom, n = 2-6, by electron impact, in the distorted-wave approximation using the varia-

Conclusions
Calculations were carried for excitation from 1S state to nS states of the hydrogen atom, n = 2-6, by electron impact, in the distorted-wave approximation using the variational polarized-orbital method [10].An approximation is called a distorted-wave approximation when the channel coupling between the initial and final states is neglected, as in the present calculation.Cross-sections have been calculated from k = 0.87 to 3.0.Present results have been compared to those obtained using other approaches.Rate coefficients for excitation, required for diagnostic purposes, are calculated using the excitation crosssections.A comparison of the electron impact excitation cross sections with the positron impact excitation cross sections was made because positrons are found to be important in solar physics [26].The cross-sections in the threshold region are shown to be proportional to (ln k f ) −2 , as shown in Figure 6.The calculations were carried out in the quadruple precision.Therefore, results obtained are expected to be accurate.
Funding: The research received no external funding.

Figure 3 .
Figure 3. (Color online) Total 1S to 4S excitation cross sections ( 2 0 a π ) vs. the incident momentum using spherically symmetric wave functions and R-matrix calculations of Scott et al

16 Figure 6 .
Figure 6.(Color online) 2S excitation cross sections ( 2 0 a π ) as a function of the incident k (Ryd) in the threshold region, the solid line represents the fit to the squares (calculated cross sections).

Figure 6 .
Figure 6.(Color online) 2S excitation cross sections (πa 2 0 ) as a function of the incident k (Ryd) in the threshold region, the solid line represents the fit to the squares (calculated cross sections).

Table 1 ,
and the lower limit in the integral is the threshold ener excitation.The excitation rate coefficients for excitation to the nS states, n = 2-6, are

Table 1 .
Cross sections (πa 2 0 ) of nS states of atomic hydrogen using the variational polarized orbital method with maximum L max = 15.
a Total: represents the sum of all the excitation cross sections.b Total1: represents the sum of excitation cross sections and the elastic cross sections.

Table 4 .
Excitation rate coefficients for nS states, n = 2 to 6, at various temperatures.

Table 5 .
Comparison of the positron and electron impact excitation of the sum of the n = 2 to 6 states of the atomic hydrogen.