Photoionization and Electron-Ion Recombination of n = 1 to Very High n -Values of Hydrogenic Ions

: Single electron hydrogen or hydrogenic ions have analytical forms to evaluate the atomic parameters for the inverse processes of photoionization and electron-ion recombination (H I + h ν ↔ H II + e) where H is hydrogen. Studies of these processes have continued until the present day (i) as the computations are restricted to lower principle quantum number n and (ii) to improve the accuracy. The analytical expressions have many terms and there are numerical instabilities arising from cancellations of terms. Strategies for fast convergence of contributions were developed but precise computations are still limited to lower n . This report gives a brief review of the earlier precise methodologies for hydrogen, and presents numerical tables of photoionization cross sections ( σ PI ), and electron-ion recombination rate coefﬁcients ( α RC ) obtained from recombination cross sections ( σ RC ) for all n values going to a very high value of 800. σ PI was obtained using the precise formalism of Burgess and Seaton, and Burgess. α RC was obtained through a ﬁnite integration that converge recombination exactly as implemented in the uniﬁed method of recombination of Nahar and Pradhan. Since the total electron-ion recombination includes all levels for n = 1 − ∞ , the total asymptotic contribution of n = 801 − ∞ , called the top-up, is obtained through a n − 3 formula. A FORTRAN program “hpxrrc.f” is provided to compute photoionization cross sections, recombination cross sections and rate coefﬁcients for any nl . The results on hydrogen atom can be used to obtain those for any hydrogenic ion of charge z through z -scaling relations provided in the theory section. The present results are of high precision and complete for astrophysical modelings.


Introduction
Photoionization and electron-ion recombination are inverse processes which can be described as where a hydrogenic ion X ion of charge +z is being ionized by a photon of energy hν.e is the ejected photoelectron or recombining ion.As a single electron system, there is no doubly excited autoionizing state and hence no resonances in these two processes.These are the two most dominant radiative processes along with photoexcitations and electron-impact excitations in astrophysical plasmas.Astrophysical modeling for spectral analysis for abundances, ionization fractions, diagnostics for the physical and chemical conditions, plasma opacity, etc. require electron-ion recombination rate coefficients and photoionization cross sections, particularly if the plasma is around or near a radiative source, such as a star where plasma is being ionized by photons.In the near-empty space, the emission lines of hydrogen from high-lying n-levels are produced by electron-ion recombination as an electron combines with the ion and cascades down to lower levels.The process also increases the intensity of lines from low lying levels by populating the levels.In dense plasmas, while photoexcitations, deexcitations through bound-bound transitions appear as lines, photoionization and electron-ion recombination play important role in the shape and enhancement in the spectral background.
As the most abundant element in the universe, precise values for the atomic processes are needed for all applications.Photoionization and electron-ion recombination of a hydrogen atom have continued to be studied by many for almost a century, such as earlier in the past by Kramers [1] and recently by Kotelnikov and Milstein [2], Rosmej et al. [3].The solution for photoionization cross sections (σ PI ) and recombination rate coefficients (α RC ) for hydrogenic ions can be expressed exactly.Studies by several investigators, such as Sugiura [4], Gordon [5], Gaunt [6], Menzel and Perkis [7], Bethe and Saltpeter [8], Burgess [9], Seaton [10], Burgess and Seaton [11] worked for the exact general expression for the total photoionization cross section σ PI (n, z, hν).However, a major difficulty is faced in the evaluation of the electric dipole matrix elements for the bound-free transitions nl − kl , where k is the momentum of the ejected photoelectron.The transition integral leads to Gaunt factor which has infinite terms and cancellation of contributions of terms which cause slow computation and numerical instabilities.
Recurrence relations were introduced for faster convergence of the infinite number of terms in the bound-free integral (e.g., [12]).However, implementation has remained limited.For higher accuracy with faster approaches for photoionization cross sections and electron-ion recombination rates of hydrogenic ions, the investigation has continued.Omidvar and Guimaraes [13] derived the formulations for optical oscillator strengths for bound-free transitions, f (n, l, n 2 ) by analytic continuation of the bound-bound oscillator strengths.They provided tabulated optical oscillator strengths for n ≤ 50.The achievement of high precision has remained a challenge.Computations have been limited to cover a fairly small range of n or are introduced with some simplifying approximation.
The co-developer of the astrophysical modeling code, MAPPINGS, Ralph Sutherland of the Australian National University wrote: "...I am having trouble with hydrogen (!), and especially the higher Balmer (n > 10) lines as observed in HII regions in deep observations Mike (Dopita) made.They are much too bright in the observations, and even give negative reddening solutions in contradiction to other reddening indicators.... MAPPINGS has the poor situation ...I cannot change collisions, or l-mixing rates, include pumping, or indeed simply run a fully self-consistent large n − l hydrogen model.I'm hoping to make a multi-level atom model with better collisions and opacities to resolve both high n lines as well as nebula continuum in even far from equilibrium conditions" (private communication, 2021).
The present work aims to present state-specific photoionization cross sections, electronion recombination rate coefficient of hydrogen ion for a very large number of n going up to 800, total recombination cross sections, and relevant codes that can be incorporated to astrophysical modelings and extend to hydrogenic ions of any charge z.

Existing Values for Photoionization Cross Sections
Green et al. [14] calculated photoionization cross sections σ PI of hydrogen and provided earlier values of σ PI in tabular form.Karzas and Latter [15] provided the Gaunt factors relevant to the ground state photoionization of hydrogen.Using them, Ditchburn and Öpik [16] reported the first accurate photoionization cross section for the ground state of hydrogen.They obtained σ PI at ionization threshold to be 6.3 Mb.This was verified experimentally by Palenius et al. [17] and Norwood and Ng [18].
Various codes use central field approximation to compute background feature of σ PI for complex atomic systems (e.g., [19,20]).They use distorted wave approximation to describe the continuum electron.Such programs can provide σ PI for hydrogen to a good approximation.However, use of Gaunt factors are needed for precise values.Using Gaunt factor and recurrence relations of Burgess [12], Seaton computed the first set of high accuracy photoionization cross sections of nl levels, where n = 1 -10 and l = n − 1, of hydrogen under the Opacity Project (OP) (The Opacity Project Team [21]).OP carried out the first systematic detailed study of photoionization of the ground and many excited states of the astrophysically important multi-electron atomic systems, from hydrogen to iron.Seaton's cross sections were not published, but are available at the OP database, TOPbase [22].Atomic data for the same levels with n = 1 -10 are also available with a finer energy mesh at the database NORAD-Atomic-Data [23,24].

Existing Values for Electron-Ion Recombination of Hydrogenic Ions
The electron-ion recombination coefficients for n-shells, which are summed over llevels, for hydrogen were obtained by Seaton [10], and later on for nl levels by Burgess [12].Using the expressions on the Kramers-Gaunt factor g I I by Seaton [10], Flower and Seaton [25] developed the code which produced accurate radiative recombination rates for hydrogen.Brocklehurst [26] extended the code further to compute radiative recombination rates of hydrogenic ions for the study of level populations in gaseous nebulae.Brocklehurst's work was used later by Hummer and Storey [27] for computation of radiative recombination rates of hydrogenic ions in their study of recombination lines.Hummer [28] reported total recombination and energy loss coefficients for hydrogenic ions which are being used widely till today.In Ferland et al. [29] for the study of broad line regions of active galactic nuclei, Nahar computed the level specific recombination rate coefficients of hydrogen up to n = 10, but introduced a finite limit integral approach to obtain the precise values for the rates.The approach was implemented for the top-up contributions of nl levels with n = 4 -10 and l = 3 -9 of the unified method of electron-ion recombination of Nahar and Pradhan [30,31] and Nahar [32] for the initial applications of the method.All these work do not provide the level-specific recombination rate coefficients and high-n data needed for various astrophysical applications.
Following the study of electron-ion recombination of various ions by Nahar [33], importance of the contributions of the high-n recombination at low temperature was realized.Nahar computed those contributions from hydrogen atom and used the data to compute recombination of various hydrogenic ions through z-scaled relation.The contributions were incorporated to revise values for recombination rate coefficients in [33].The present work extends those works and develops the code "hpxrrc.f" to compute photoionization and electron-ion recombination for any nl level of hydrogen and study the characteristic features.

Theory of Photoionization and Electron-Ion Recombination of Hydrogenic Ions
Brief descriptions of the theoretical background of the two inverse processes are given below for the general guidance of the present work.

Photoionization of Hydrogenic Ions
A brief review of the theory is needed to understand the present approach and the results presented.The first treatment for photoionization process was based on classical approximation.For low electron energy, the motion of the electron in hydrogen can be described by a classical trajectory to a good degree of accuracy.Using semi-classical approach, Kramers [1] obtained the following simple expression for the ground state photoionization cross sections, σ PI , of hydrogen where ω is the photon frequency and n is the principal quantum number.The above equation has been derived later by others, e.g., by Burke [34] and Kogan et al. [35] but differs in format.The above equation is similar to what Burke obtained except multiplying by a 2 o .Kramers cross section at the ionization threshold of hydrogen is known to be about 7.8 ×10 −18 cm 2 , e.g., Rosmej et al. [3].This is not the exact value of the cross section, 6.3 ×10 −18 cm 2 mentioned above, but gives the correct behavior at higher energies.The formula is often used to obtain approximated σ PI at very high energies and for highly excited states of hydrogenic and multi-electron systems e.g., [36].
The energy of the photon, hν, ionizing a hydrogenic ion of charge z can be expressed in terms of ionization threshold energy of hydrogen, I H as hω = hν = hRcz 2 1 where energy of the ejected electron of a hydrogenic ion of charge z is k 2 = z 2 in terms of hydrogen electron energy and Rydberg units (I H ≈ 13.6 eV).Using the above photon energy, the total photoionization cross section for shell n of a hydrogenic ion is given by [10] where g I I is the Kramers-Gaunt factor and is of order unity for hydrogenic ions.This expression is similar to what Kramers [1] obtained for the ground state photoionization of hydrogen except for the gaunt factor g I I .If g I I is set exactly equal to 1, then the expression reduces to that of Kramers [1].For photoionization of hydrogen at the ionization threshold of the ground state, z = 1, n = 1, k = 0.If we choose g I I = 1 and use α = 7.297×10 −03 and a o = 0.529 −08 cm, the above expression gives σ PI (1, 1, I H ) = 7.8×10 −18 cm 2 which is same as that from Kramers formula [1].g I I can be expressed using hypergeometric functions [7], and is difficult to evaluate exactly.The asymptotic expansion for g I I derived by Menzel and Pekeris [7] and corrected by Burgess [9] has the form (given in [10]) where u = n 2 .The expression indicates the complexity in computing the infinite number of terms of the series and the instability that may arise from cancellation of values.Let the ion being ionized be in the quantum state specified by the principal quantum number n and orbital quantum number l. Photoionization cross section of a hydrogenic atom for subshell nl is given by ( [11], details can also be found in [36]) where and the bound-free integral g is where α = 2πe 2 /hc is the fine structure constant, K = k/z, l is the angular momentum quantum number of the ejected electron.4παa 2 o /3 = 8.5594 × 10 −19 cm 2 .P nl (r) and F kl (r) are the initial and final radial wave functions of the ejected electron, which satisfy and are normalize such that ∞ 0 P nl (r)P n l dr = δ nn , and F kl (r) being normalized to asymptotic amplitude k −1/2 .As mentioned earlier, the main numerical challenge arises in calculating the exact bound-free dipole allowed transition matrix nl → kl in contrast to using approximations, such as distorted wave approximation.
To evaluate the general bound-free integral, Burgess and Seaton [11] introduced ρ → zr and K = k/z which can reduce the equations from direct z-dependence, and substituted the function, in the Schrödinger equation and obtained The functions are p nl and f Kl normalize such that Then the bound-free integral [11,12] g(n, l : does not explicitly depend on z.The solutions p nl (ρ) and f Kl (ρ) can be written as ( [5,8]) and where 1 F 1 are confluent hypergeometric functions.Substitution of the functions in the found-free transition integral leads to [5,8], given in [12]) where and In the above equations, η = 1/K and 2 F 1 are hypergeometric functions.The expressions are not suitable to compute in general because of the large number of terms occurring in the series when n − l is large, and also for some ranges of the parameters, there is gross cancellation between the terms of the series [12].
The following subsection gives the recurrence relations for fast computation of the g-integral.They have been implemented in the program for the present photoionization cross sections and can be used for any energy.However, very high energy region can still have some numerical instability.Since the treatment provides σ PI which agrees with those from Kramers formula, σ PI at high energies can easily be obtained using Kramers formula or Seaton's formula without the Gaunt factor from a known value at a lower energy as where σ 1 is the known cross section at photon energy hν 1 .
For very high n, beyond n = 800 as considered in the present work, the value of σ PI may not have any practical importance as they give single point or a few non-zero values.If there is a need, program "hpxrrc.f" can be used for σ PI for n beyond 800.Kramers formula or Seaton's formula without g I I can also be used.
The cross section for each nl level belonging to n, σ pi (nl, k 2 ), is related to the total cross section, σ PI (n, z, hν) through the summed contributions as 2.1.1.Recurrence Relations for the g Bound-Free Transition Integral The instability of the bound-free g integral can be removed through use of its recurrence relations that satisfy the exact matrix elements.Biedenharn et al. [37] and Alder and Winther [38] introduced the recurrence relations of g(n, l; K, l') connecting matrix elements having different values of l (see [12]), Using these relations, Burgess [12] adopted an alternate procedure to calculate the transition integrals.He simplified the relations with approximation of the hypergeometric functions and produced a revised set of recurrence relations satisfied by the exact transition matrix elements and enable them to be evaluated rapidly and to a high accuracy.He defined where the quantities G(n, l; K, l ) satisfy the recurrence relations Since repeated use of the two G relations does not involve division or the evaluation of square roots, this scheme is very suitable for the fast computing.These relations have been used by others, such as [25].

Ground State Photoionization Cross Sections of Hydrogen
The exact analytical form for the ground state photoionization of a hydrogenic ion is available.The electron at ground state 1s can photoionize only to continuum p.The exact expression can be obtained using the bound state wavefunction of 1s and the radial wavefunction of the continuum electron with l = 1 in a Coulomb field.Using the expression for the confluent hypergeometric function F(α, γ, x), photoionization cross section of the ground state of hydrogen can be obtained as (e.g., [16,39,40]) σ o is the threshold cross section for hydrogen and hν 1 = z 2 hν o ≈ 13.6z 2 eV is the ionization threshold energy.The value of σ o was obtained first by Ditchburn and Öpik (1962).Photoionization cross section of hydrogen atom in the ground state, can also given by a simpler formula [41] σ where λ is the photon wavelength in Angstrom .The relevant Gaunt factors, which depends also on λ, for the ground state photoionization, are available from Karzas and Latter (1961).

Electron-Ion Recombination of Hydrogenic Ions
The theoretical model for recombination of hydrogenic ions used in the present work is based on the unified method developed by Nahar and Pradhan [30][31][32].In the past significantly more emphasis was given to calculating the radiative rate coefficient of hydrogen than photoionization.As the inverse processes, recombination cross section σ RC can be obtained from the photoionization cross section σ PI through the principle of detailed balance, where g i and g j are the statistical weight factors of the atom being ionized and of the residual ion.The photon momentum is p hν = hν/c and the photoelectron momentum is p e = mv = √ 2m .Hence, It may be noted that σ RC diverges at photoelectron energy approaches zero.Most astrophysical applications need the temperature dependent quantity, the recombination rate coefficient α RC (T).This is obtained from averaging σ RC over Maxwellian distribution function of electrons at temperature T, Using f (v, T), = 1 2 mv 2 , and d = mvdv, The expression, Equation (33) introduces some uncertainty at and near zero photoelectron energy where σ RC diverges.The uncertainty is removed by replacing σ RC with σ PI through the principle of detailed balance as [31] α where hν = I + and I is the ionization energy.
Although the integration for α RC (T) can be carried over a large energy range as done in other computations, ideally the integration extends to infinity where the background photoionization cross section continues to contribute, although by a very small amount.To include the full contributions, the integral is transformed for finite limits.Introducing x = exp(− /kT), which has limits between 0 and 1, and a slow variation at low temperature and a relatively fast variation at high T, the following fast convergent form of α R (T) can be obtained (author's contribution in [29], also see [31,32]) The evaluation of the integration is carried out over a number of energy regions where the x-mesh changes depending on the variation of the integrand.The mesh is smaller for a fast variation and larger for a slow variation.This ensures higher precision and inclusion of full contribution from photoionization cross sections to the recombination rates.
The state specific recombination rate coefficient, α RC (nl, T) for each nl level can be obtained using the above expression.The program of Nahar [33] has been extended in the present work to compute state specific nl rates.They are needed for determination of level populations and various diagnostics.
The total recombination rate coefficients α RC (T), which is the sum of all state-specific recombination rate coefficients, is needed to obtain various plasma quantities, such as ionization fractions.Since there are an infinite number of bound levels where recombination can take place, α RC (T) requires contribution from infinite number of levels, that is,

Top-Up Contribution from Very High-n Recombination
It is computationally prohibitive to compute α RC (nl, T) for infinite number of terms, particularly when contributions become significantly negligible, for the total rate α RC (T).We can obtain α RC (nl, T) up to a very high n value until the individual contributions become negligible for the total recombination rate.Then, the summed contribution from the rest of the levels, from a high n to infinity, which can be called the "top-up" contribution, can be obtained following the approximation made by Hummer [28].
The top-up contribution, from a high n to ∞, is generally negligible at high temperature, but could be significant in the very low temperature where low energy electrons are slow enough to recombine to the ions.At sufficiently large n, the recombination rate coefficient to an n-shell varies as α RC (n) ≈ 1/n 3 (e.g., [28]).Following Hummer [28], the asymptotic top-up contribution starting at a high n, α 0 (n), is expressed as the top-up contributions to the total α RC is given by is a useful quantity, particularly to see the total spectrum with resonances and is a measurable quantity (e.g., [42]).For hydrogen, which does not resonance, the quantity can add to the total background of a spectrum and hence be used in astrophysical modeling.The total recombination rate coefficients α RC ( ) with respect to photoelectron energy E is also a measurable quantity in a laboratory set-up (e.g.[43]) and hence can provide a test for the accuracy of the recombination method as was done by [43].α RC (E) can be obtained using the following relation where v is the photoelectron velocity.Both σ RC (E) and α RC (E) can be evaluated from recombination collision strength Ω RC which is related to σ RC (E) as (e.g., [36]) Note the advantage of collision strength is that it does not diverge with very low energy, but σ RC (E) does.

Relation between Hydrogen and Hydrogenic Ions
With a single electron orbiting, the electronic properties of hydrogen and hydrogenic ions are similar.Only the nuclear potential increases with higher z.The atomic parameters for hydrogen can be transformed in to those for a hydrogenic ion through z-scaled relations as given below.
(i) Binding and photon energies: (ii) Photoionization cross section at energy (iii) Recombination cross sections at photoelectron energy E PE (z), the same as that of hydrogen, (iv) Recombination collision strengths at photoelectron energy, E PE (z), the same as that of hydrogen, (v) Recombination rate coefficient at photoelectron energy, E PE (z), the same as that of hydrogen, (vi) Temperature T(z) is equivalent to hydrogen temperature z 2 T(H) (vii) Recombination rate coefficient at temperature T(z), It means that the temperature range for α(H, T) should be is wide enough to scale a highly charged hydrogenic ion.

Relativistic Fine Structure Splitting of σ PI and σ RC from LS Coupling
The present results have been obtained in the LS coupling approximation.Hence no relativistic fine structure effects have been taken in to account.Relativistic effects are expected to be important at very n and at very high temperature where slight perturbation, like that from relativistic effects, can make a change in the transition values.However, it is possible to obtain the fine structure components of photoionizatoin as well as recombinatoin of its LS value through an algebraic transformation procedure described by Nahar [44] for oscillator strengths of Fe II.It is also described in [36].These values will not have the contributions from relativistic corrections.
To carry out the fine structure splitting of photoionization or recombination cross sections, we can split the line strength of photoionization transition matrix in LS coupling S(LS) in to its fine structure components S(LSJ) [33,36].The corresponding cross sections σ PI (LSJ) can be obtained from the fine structure line strength S(LSJ) by multiplying the energy and statistical weight factors.Then from σ PI (LSJ), the corresponding recombination cross section σ RC (LSJ) and α RC (LSJ) can be obtained.
Although the algebraic transformation above will provide more accurate values, except for the contribution of the relativistic effects, a simpler approximation of fine structure splitting by statistical weight factors can be adopted.Typical features of σ PI and σ RC for the fine structure components of an LS value show similar pattern with some differences in resonances.For hydrogen there is no resonance.Hence to a good approximation, statistical weight factor approach can be used.The fine structure components are estimated from the LS value by multiplying the latter with the statistical fractions of the level, e.g., for fine structure level k, the fraction is (2J k + 1)/ ∑ i (2J i + 1).For very high n and high temperature, inclusion of relativistic effects will increase the accuracy.

Programs for Photoionization and Electron-Ion Recombination, and the Data Files
The author has written the FORTRAN program code, "hpxrrc.f",which can compute quantities for both photoionization and electron-ion recombination.It computes the atomic parameters for each individual l-level belonging to a n-shell as well as for each individual n.It will be available online.Followings describe the quantities that "hpxrrc.f"computes and the corresponding data files which are available.
1. Photoionization cross sections (σ PI ) for each individual l-level belonging to a n-shell and for each individual n.
The author received the program "hypho.f"from M.J. Seaton (private communication, 1991) which computed the total photoionization cross sections σ PI (n, E) for each shell n.As mentioned earlier, the contents of the program indicate the references for the program are Flower and Seaton [25] and Brocklehust [26].The author revised the program to compute σ PI (nl, E) for individual l-levels of each n along with those for each n.Photoelectron energy can be chosen up to a given highest energy specified by the user or up to a default energy which is five times higher than the ionizing threshold energy.
Data file: Photoionization data file "p0100.1-800.5ry"contains σ PI (nl, E) of l-levels in sets of each n-complex, i.e., σ PI (nl, E) of all l-levels of n at photoelectron energies going up to 5 Ry.

Total recombination cross sections (σ RC (E))
and total electron-ion rate coefficients α RC (E) with respect to photoelectron from already computed photoionization cross sections, such as those given in "p0100.1-800.5ry".These can also be obtained at a desired photoelectron energy E for which σ PI will need to be computed beforehand.
3. Electron-ion recombination rate coefficient, α RC (nl, T), of all individual l-levels of any shell n.The original program written by the author for [31,32] computed α RC (n, T) for each n-shell.It has been extended to compute α RC (nl, T) for all l-levels separately.For this computation, σ PI (nl, E) are computed during the run of the program.Hence no file containing photoionization cross section is needed to evaluate α RC (nl, T).
Results are printed as sets of n-complex, i.e., α RC (nl, T) is printed for various for l-levels belonging to the same n and at specified temperature range.Each n-set is ended with α RC (n, T) of shell n, i.e., the summed values from all l-levels as n-total.
Data file: "rc.0100" contains rate coefficients α RC (n, T) for each n-shell, from 1 to 800.Below them, the file contains the total α RC (T), i.e., the added values from n = 1 to 800, top-up contribution from n = 801 to infinity obtained using n −3 behavior, then the summed total α RC (T) of all contributions (n = 1-∞).
α RC (T) are given for a wide range of temperatures, T = 10 −3 − 10 9 K, at 137 tempera- tures in order to compute rates of hydrogenic ions of low to high charge z.
4. α RC (z, T) for any hydrogenic ion of charge z.For this, the program needs the temperature and the z-value.

Results and Discussions
Results of photoionization and electron-ion recombination for hydrogenic ions are presented for a very large n going up to n = 800.Contributions to total electron-ion recombination from beyond n = 800 to infinity are obtained using n −3 behavior.Characteristic features of the two processes are presented in two separate sections below.Although the features are demonstrated for hydrogen atom, they are similar for hydrogenic ions.The results for hydrogen can be scaled to any hydrogenic ion using the Z-scaled relations given in the Subsection 3 of the theory section.

Photoionization Cross Sections
Photoionization cross sections of hydrogenic ions for lower n are available from several codes which use various approximations.The present σ PI have been obtained from the exact expressions as described in the theory section and using the recurrence relations revised by Burgess [12].Hence the present numbers agree with those obtained for n = 1 going up to 10 by Seaton under the Opacity Project [21].As mentioned above that Seaton's data are unpublished but available at data base TOPbase at CDS [22].
Figure 1a,b from the present work demonstrate typical smooth feature of σ PI of hydrogen for the series of levels ns with n = 1-5 (left) and nd with n = 3-7 (right).σ PI of ns = 1s corresponds to the ground state photoionization.The figure illustrates that σ PI at the ionization threshold increases with n, but the background decreases faster in the higher energy with increasing n.The trend is expected since at ionization threshold, σ PI is inversely proportional to the ionization energy.The ratio of cross sections at zero photoelectron energy (k = 0) could be estimated as where n1, l is the lowest possible n for angular momenta l.Since the background cross section at higher energy deceases as (ν o /ν) 3 where ν o is the photon energy at ionization threshold, cross section decreases faster for lower ionization threshold energy I o .
Figure 2 presents σ PI of all 5 l-levels (l = 0 -4 of shell n = 5.Although all cross sections have the same threshold energy, the peaks at threshold for various levels are different, but the background decays faster with larger l.The bound-free integral in photoionization cross sections depends on the l values as g(n, l; K, l = l ± 1) ∼ ( 4n 1+n 2 K 2 ) l < +1 indicating σ PI (nl, ) variation with on the angular momentum l.However, it is not linear and hence the general trend can not be predicted.This is demonstrated in Figure 2.  The program "hpxrrc.f" can be used for σ PI of any nl levels.However, machine accuracy for very low numbers should be discussed.All computations were carried out in double precision using the clusters at Ohio Supercomputer Center (OSC).Numerical instabilities for very small values of σ PI at very high n can be summarized as: (i) With increasing energy, cross sections reach to a negligible quantity beyond which computation gives zero value.These cross sections may not add to any significant values.(ii) σ PI at threshold continues to dominate while it decreases faster with very high n.(iii) For a very high n, only a few points, sometimes only the threshold value for various l-levels survive before computer prints out zeros even with double precision accuracy.(iv) With a ver high n, such as near n = 700 or 800, even threshold energies become sensitive to the precision of the computer which could yield a lower value.There could be numerical instability in the recurrence relations for these very high n-values in addition to computer numerical tolerence issue.Seaton's relation without the Gaunt factor or Kramers relation could be useful for such situation of very high-n cross sections.

Electron-Ion Recombination
The present work reports recombination rates for all individual shells from n = 1 to 800.Some examples are given for demonstration of features.
An example of α RC (nl) with n = 10 and l = 0 -9 are illustrated in Figure 3.All curves corresponding to l = 0 -9are drawn with different colors and/or with different type (e.g., solid, dashed, or dotted).The sum of all α RC (nl) provides the "n-total" curve (red, solid) which is distinctly separated from all curves for individual l-values below it.The figure shows that at low temperature the height of α RC (nl) curves do not follow the l-order as expected from l-dependence of σ PI at threshold but they follow the order at higher temperature region.α RC (10s, T) of the lowest 10s level does not have the highest peak at lower temperature, but dominates at higher temperature region.The high temperature dominance comes from the slower decrease in σ PI of 10s at higher energies compared to the higher l-levels.Couple of small kinks, visible at logarithmic scale, can be noticed for two l-curves (8 and 9) at very high temperature.They arise from some numerical instability of very small σ PI values.Although α RC (nl, T) for various l-levels belonging to the same n-complex show variations in dominance at lower temperature (Figure 3), the summed total α RC (n, T) exhibit systematic decrease in values with increasing n. Figure 4 illustrates this behavior with α RC (n, T) of shells n = 1 -10.Various colors/type of curves represent different n as specified in the figure.The topmost curve corresponds to the ground state n = 1 and the lower ones are sequentially in order for n = 2 -10.The right panel shows logarithmic behavior of the α RC (n, T) values.To demonstrate the general linear trend of the α RC (n, T) values, the left panel presents the linear plot of the same α RC (n, T) in the same temperature range of T = 10 −3 to a million degree 10 6 .The curves show sharp rise in α RC (n, T) at lower temperature for all n = 1 -10with more clarity than in the log-scale.
As mentioned above, the present work provides a file which contains α RC (n, T) of shells from n = 1 going up n = 800 of hydrogen atom.It does not include components α RC (nl, T) for each n.However, the individual α RC (nl, T) for l-levels which can be computed using the program "hpxrrc.f" as needed.Figure 5a presents the curve for the total recombination rate coefficients of hydrogen, α RC (H, T), which include contributions of individual shells from n = 1 to 800 and the top-up contributions from n = 801 to infinity.Table 1 presents the numerical values of α RC (H, T).The range of temperature for α RC (H, T) is chosen wide, from 10 −3 to 10 9 K.This allows computation of α RC (z, T) for hydrogenic ions over a temperature range of practical need.The temperature for an ion of nuclear charge z is equivalent to z 2 T(H).Hence, obtaining α RC (z, T) of a hydrogenic ion of charge z, such as hydrogen-like neon with z = 10, at temperature 1000 K will need α RC (H) at 10 K.As in illustration of calculating α RC (z, T) of a hydrogenic ion, Figure 5b presents α RC (FeXXV I, T) of Fe XXVI obtained from α RC (H, T) using z-scaling.
The total recombination cross sections, σ RC (E), and recombination rate coefficients α RC (E) of hydrogen with respect to the photoelectron energy E are presented in Figure 6.The energy unit is chosen to be in eV because of use of eV unit for the photoelectron in laboratory measurements.Both σ RC (E) and α RC (E) include contributions of n = 1-∞.
Except for the sharp rise in the very low photoelectron energy, the features show smooth decrease with energy.
Except for the sharp rise in the very low photoelectron energy, the features show smooth decrease with energy.The present work employs a finite limit integral approach that limits between 0 and 1 replacing the direct dependency on the photoelectron energy going out to infinity.This reduces the uncertainty arising from the limit of energy set for computation of the recombination rate coefficients Hence present results are expected to be most accurate compared to the existing values.The accuracy of the recombination atomic data also depends on the values of photoionization cross sections.As explained in the theory section that the analytical form of σ PI (H) has numerical instability which is compensated by the recurrence relations.For very high n. the recurrence relations may have some uncertainties since the cross sections decreases sharply.This may have contributed some uncertaintis in the rates.As explained above that the accuracy of the values also depends on the precision and numerical tolerance of the computer being used.With the best estimation of the factors, that lower-n levels, where values are more stable and accurate, contribute more than the highern levels, use of precise theory and numerically stable recurrence relations, and finite limit integration for recombination, the current results are expected to have an accuracy within 5% for most temperature range.

Conclusions
The present report can be summarized as 1.Study of the two inverse processes of photoionization and electron-ion recombination are presented with a brief review of theory that treat them precisely and more accurately compared to all existing atomic structure codes that use other approximation, mainly distorted wave approximation.2. Detailed features of both the processes are illustrated.Although hydrogen and hydrogenic ions do not have any resonant features, accurate characteristic variation with energy and temperature are crucial for precise astrophysical spectroscopy and modeling.3. The present work provides atomic data files containing σ PI (E) for all l-levels of n from 1 to 800, and α RC (T) of all values of n from 1 to 800.It also provides total α RC (H, T) with temperature, and total α RC (H, E), and total σ RC (E) with photoelectron energy.This is the first time that all these data with very high n have been made available for applications.4. Use of precise theory, numerical methods, high precision computers, as explained at the end of the sections of photoionization, and electron-ion recombination, can predict the accuracy of the present results is within 5% for most energy and temperature ranges.5.The work provides the FORTRAN program "hpxrrc.f" which can generate all these values.It also computes l-level specific α RC (nl, T).The program "hpxrrc.f" can also compute α RC (H, T) for any hydrogenic ion of charge Z using the data of hydrogen.6. Importance of relativistic effects and how to obtain fine structure components for photoionization and electron-ion from their values in the present LS coupling approximation have been discussed.7.All atomic data and the program will be available online at database, NORAD-Atomic-Data ( [24], http://norad.astronomy.ohio-state.edu(accessed on 1 September 2007)).

Figure 2 .
Figure 2. Photoionization cross sections (σ PI ) of the five l-values (0-4) of shell n = 5.They demonstrate rising trend of σ PI for various l-values compared to l = 0 (black curve) at the same ionization threshold.

Figure 3 .
Figure 3.Level specific electron-ion recombination rate coefficients (α RC (nl)) of the 10 individual l-levels (0 -9)of shell n = 10, and the n-total (red curve with highest values) which is the summed total of all 10 individual levels.They demonstrate α RC (nl) peak does not follow values of l-order at low temperature bu follows at higher temperature.

Figure 4 .
Figure 4. Total electron-ion recombination rate coefficients (α RC (n, T)) of shells n = 1-10, which are the summed total of all individual l-levels, Left) in linear scale to see the general trend at low temperature and Right) in log scale for logarithmic features.The curves demonstrate α RC (n, T) decreasing with increasing n.

Figure 5 .Figure 6 .
Figure 5. (a) Total electron-ion recombination rate coefficients α RC (H, T), summed contributions from n = 1 to infinity over a wide temperature range.These values are used to obtain (b) α RC (FeXXV I, T) of Fe XXVI using the Z-scale formula.

Figure 5 .Figure 6 .
Figure 5. (a) Total electron-ion recombination rate coefficients α RC (H, T), summed contributions from n = 1 to infinity over a wide temperature range.These values are used to obtain (b) α RC (FeXXV I, T) of Fe XXVI using the Z-scale formula.
39) 2.2.2.Recombination with Respect to Photoelectron Energy Recombination cross sections for a level σ RC ((n, E) with respect o photoelectron energy E shows the dielectronic satellite lines, if there are resonances.The total recombination cross section σ RC (E) (summed over σ RC ((n, E) of infinite number of recombining n-levels) with respect to photoelectron energy E σ RC

Table 1 .
Total recombination rate coefficients of hydrogen, α RC (H, T), at a wide range of temperatures.