BELFASTBelfast Data Base on Atomic and Molecular Physics Department of Computer Science The Queen's University of Belfast Belfast BT7 1NN, U.K.ORNL-CFADCControlled Fusion Atomic Data Center Oak Ridge National Laboratory M.S. 6373, P.O. Box 2008 Oak Ridge, TN 38731-6372 USAPPPLPrinceton Plasma Physics Laboratory P.O. Box 451 Princeton, NJ 08544 USARIC-NAGOYAResearch Information Center Institute of Plasma Physics Nagoya University Chikusa-ku, Nagoya 464, JapanIAEA-AMD-UNITAtomic and Molecular Data Unit Nuclear Data Section International Atomic Energy Agency P.O. Box 100 Wagramer Strasse 5 A-1400 Vienna, AustriaJAERI-NDCNuclear Data Center Japanesse Atomic Energy Research Institute (JAERI) Tokai-Mura Naka-Gun Ibaraki-Ken 319-11, JapanNIST-AEL-DCAtomic Energy Levels Data Center National Institute of Standards and Technology United States Department of Commerce Physics Building, Room A167 Gaithersburg Maryland 20899, USANIST-TP-DCTransition Probabilities Data Center National Institute of Standards and Technology United States Department of Commerce Gaithersburg Maryland 20899, USACRAAMDChinese Research Association for Atomic and Molecular Data P.O. Box 8009 Telex: 222499 IAPCM CN Beijing The People's Republic of ChinaDOC=H-HE-PLASMA"Elementary Processes in Hydrogen-Helium Plasmas," R.K. Janev, W.D. Langer, K. Evans and D.E. Post, Jr., Springer Series on Atoms and Plasmas, Springer-Verlag (Berlin, Heidelberg, New York) 1987.DOC=IPPJ-AM-52"Energy Dependence of Ion-induced Sputtering Yields of Monatomic solids in the Low Energy Region", N. Matsunami, Y. Yamamura, N. Itoh, H. Tawara and T. K. Kawamura, Report IPPJ-AM-52, Institute of Plasma Physics, Nagoya University, July, 1987.DOC=IPPJ-AM-56"Total and Partial Cross Sections for Electron Capture for C[q+] (q=6-2) and O[q+] (q=8-2) Ions in Collisions with H, H[2] and He Atoms," H. Tawara, Report IPPJ-AM-56, Institute of Plasma Physics, Nagoya University, October, 1987.DOC=JPCRD-12/83"Recommended Cross Sections and Rates for Electron Ionisation of Light Atoms and Ions," K.L. Bell, H.B. Gilbody, J.G. Hughes, A.E. Kingston and F.J. Smith, J. Phys. Chem. Ref. Data 12, 891-916 (1983).DOC=JPCRD-17/88"Recommended Data on the Electron Impact Ionization of Atoms and Ions: Fluorine to Nickel," M.A. Lennon, K.L. Bell, H.B. Gilbody, J.G. Hughes, A.E. Kingston, M.J. Murray and F.J. Smith, J. Phys. Chem. Ref. Data 17, 1285-1363 (1988).DOC=NUC-FUS-SUPP/87"Recommended Data on Atomic Collision Processes Involving Iron and its Ions," C. Bottcher, D.C. Griffin, H.T. Hunter, R.K. Janev, A.E. Kingston, M.A. Lennon, R.A. Phaneuf, M.S. Pindzola and S.M. Younger Special Supplement to Nuclear Fusion, 1987.DOC=ORNL-6090Atomic Data for Fusion, Volume 5 "Collisions of Carbon and Oxygen Ions with Electrons, H, H[2] and He" R.A. Phaneuf, R.K. Janev and M.S. Pindzola Oak Ridge National Laboratory Report ORNL-6090, February, 1987.DOC=CLM-R294/1989"Atomic and Molecular Data for Fusion, Part 3 "Recommended Cross Sections and Rates for Electron Impact Ionization of Atoms and Ions: Copper to Uranium ", M.J. Higgins, M.A. Lennon, J.G. Hughes, K.L. Bell, H.B. Gilbody, A.E. Kingston and F.J. Smith, Culham Laboratory Report, CLM-R294, UKAEA, Culham Laboratory, Abingdon, Oxon. OX14 3DB, England.DOC=ADNDT,33,149,(1985)"Electron-Impact Cross Sections and Rate Coefficients for Excitations of Carbon and Oxygen Ions," Y. Itikawa, S. Hara, T. Kato, S. Nakazaki, M. S. Pindzola, D. H. Crandall, Atomic Data and Nuclear Data Tables, 33,149,(1985).DOC=ADNDT,42,313,(1989)"Recommended Data for Excitation Rate Coefficients of Helium Atoms and Helium-like Ions by Electron Impact," T. Kato and S. Nakazaki, Atomic Data and Nuclear Data Tables, 42,313,(1989).DOC=ORNL-6086Atomic Data for Fusion, Volume 1, "Collisions of H, H{2}, He and Li Atoms and Ions with Atoms and Molecules," C.F. Barnett, Oak Ridge National Laboratory Report ORNL-6086, July, 1990.DOC=IAEA-AM-REP,249,(1991)"Particle Reflection from Surfaces - A Recommended Data Base", E.W. Thomas, R. K. Janev , J. J. Smith, International Atomic Energy Agency, Atomic and Molecular Data Unit Report IAEA-AM-REP,249 June 1991.#MEWEGENERALIZED GAUNT FACTOR LINE EXCITATION CALCULATION REFERENCE: R. MEWE, ASTRON. AND ASTROPHYS. 20, 215 (1972) INPUT: ELECTRON TEMPERATURE (keV) RETURNS: SPECTRAL LINE EXCITATION RATE COEFFICIENT (CM3/S) #MEWE DATA TYPE COEFFICIENTS ARE (IN ORDER): 1) EXCITATION WAVELENGTH (ANGSTROMS) 2) EXCITATION ENERGY (EV) (IF PEIJ=0, ASSUME A RESONANCE LINE AND CALCULATE EIJ FROM PWVL) 3) OSCILLATOR STRENGTH 4) BRANCHING RATIO 5-8) A, B, C, D COEFFICIENTS FOR MEWE GAUNT FACTOR FORMULA#TAB1DGeneral one-dimensional tabular data type Arbitary number of (X,Y) data pairs are supplied as coefficients X1, Y1, X2, Y2, ... These (x,y) data pairs must be entered in order with X increasing. Evaluation of #TAB1D data uses linear interpolation between the (X,Y) data pairs to return a Y value for a specified X. If the given X is below the first X data value or larger than the last X data value, an error return occurs.#CHEBORNL Chebyshev polynomial fitting subroutine. The parameter list consists of nine polynomial coefficients followed by two additional parameters which define the lower and upper limits for the independent variable over which the fit is valid. For a reference, see page 40 of the document DOC=NUC-FUS-SUPP/87.#KINGIAEA subroutine to evaluate the polynomial fit used to represent electron-impact excitation rate coefficients from Kingston et al. (See DOC=NUC-FUS-SUPP/87, page 47 ). The first three coefficient data values are 1. The excitation energy for the transition in eV. 2. The lower limit of validity for the fit. 3. The upper limit of validity of the fit. The remaining data values, of maximum number 9, are the fitting coefficients for the polynomial.#PHACXIAEA subroutine to evaluate the analytic function used to represent the charge transfer cross sections by R. A. Phaneuf et al. (See DOC=NUC-FUS-SUPP/87, page 14 ). The coefficient data values are (in order); 1. The lower limit of the range of validity for the scaled projectile energy (in eV/(u*sqrt(q), where q is the charge state of the projectile ion). 2. The upper limit of the range of validity for the scaled projectile energy (in eV/(u*sqrt(q)). 3-6. The four parameters for the analytic fit. As the seventh parameter in the coefficient array passed to #PHACX the charge state of the projectile ion must be passed. For the data referenced in DOC=NUC-FUS-SUPP/87 the limits for q are +5 <= q <= +26.#BELIIAEA subroutine to evaluate the analytic fit for the electron-impact ionization cross sections from Bell et al. (See DOC=JPCRD-12/83). The number of fitting parameters specified in the coefficient data varies depending on the number of terms taken in the numerical fit and on whether allowance for excitation-autoionization (E-A) has been made. For cross sections which include E-A contributions, two seperate fits are defined. One from the ionisation threshold to the threshold for E-A and a second fit for energies above the E-A threshold. (The number of parameters (in any entry) is given in the variable KNCF). The coefficient data values are For the first region 1. The ionisation potential (eV) of the target. 2-7. The fitting parameters for the fit. (A, B1, B2 etc of Bell et al) For the second region (for cross sections with E-A only) 8. The E-A threshold energy (eV). 9. An ionisation potential (eV) for the region. 10-15. The fitting parameters for the fit. (A, B1, B2 etc of Bell et al)#JAN1IAEA subroutine to evaluate the 9-term polynomial fit used in the book of Janev et al. (see DOC=H-HE-PLASMA). The first and second parameters of the coefficient data are the lower and upper range of validity of the fit, for cross sections Emin to Emax and for rate coefficients Tmin to Tmax. This range data is followed by the nine coefficients for the polynomial. The Boolean labels included with this evaluation function are a. For cross sections, values at Emin in label XS(EMIN) while for reaction rate coefficients the value at Tmin in label RC(Tmin). b. For cross sections, the maximum valure of the cross section in label XSMAX, while for rates the maximunm rate in label RCMAX. c. The fitting error, ERROR. (see DOC=H-HE-PLASMA for meaning).#ETABSQIAEA subroutine to calaculate the yields (atoms/ions) of ion-induced sputtering of monatomic solids of N. Matsunami et al. (see DOC=IPPJ-AM-52). The data coefficients that must be passed to #ETABSQ are: 1. The atomic number of the incident ion 2. The atomic mass of the incident ion 3. The atomic number of the target 4. The atomic mass of the target 5. Sublimation energy of the target (in eV) 6. The Q value, defined by in IPPJ-AM- 52, for the target The coefficient data for this evaluation function are read by the evaluation function #RTABSQ.#RTABSQIAEA subroutine to read and interpret the entry data for the yields (atoms/ions) of ion-induced sputtering of monatomic solids of N. Matsunami et al. (see DOC=IPPJ-AM-52). This routine requires as input (in order) the specification of three atomic species of the incident ion and of the target. Ther subroutine returns the coefficient data required by #ETABSQ to determine the sputtering yield.#JANRDIAEA subroutine to evaluate the double 9-term polynomial fit used in the book of Janev et al. (see DOC=H-HE-PLASMA) for reaction rate coefficients (cm[3]/s) as a function of the particle impact energy (eV) and the plasma temperature (eV) for heavy particle reactions. The first and second parameters of the coefficient data are the lower and upper range of validity for the plasma temperature (in eV) of the fit, while the third and fourth give the range of validity for the energy of the incident particle (in eV). This range data is followed by the 81 coefficients for the polynomial. The order of apperance of the coefficiennts in the coefficient array is all 9 coefficients for each incident particle energy, ie A(E=0,T=0-9), A(E=1,T=0-9),... , A(E=9,T=0-9). The coefficient array is temporarily updated whren calling JANRD. PCF(86) = E, the energy of the incident particle. PCF (87-95) stores temperature independent data (ln (E)) ** i , where is 0 to 9. See DOC=H-HE-PLASMA and evaluation function for details. The Boolean labels which may be included with this evaluation function are 1. The fitting error, ERROR. (see DOC=H-HE-PLASMA for meaning). 2. The error class for the fit. (see DOC=H-HE-PLASMA for meaning).#JBORN1IAEA subroutine to evaluate the empirical Born-Bethe type fit Born,1 used for electron impact cross sections (cm[2]) as a function of electron energy in the book of Janev et al., (see DOC=H-HE-PLASMA). This fit is derived using uniform weighting in the least-squares procedure. The coefficient data values are 1. The threshold energy for the reaction in eV. 2. The minimum energy , Emin, for application of the fitting formula. 3. The coefficient aBorn,1 used in DOC=H-HE-PLASMA. 4. The coefficient nBorn,1 used in DOC=H-HE-PLASMA. The Boolean labels include a. The cross section at Emin in label XS(EMIN) b. The fitting error, ERROR. (see DOC=H-HE-PLASMA for meaning).#JBORN2IAEA subroutine to evaluate the empirical Born-Bethe type fit Born,2 used for electron impact cross sections (cm[2]) as a function of electron energy in the book of Janev et al., (see DOC=H-HE-PLASMA). This fit is derived using a weighting which does not take into account the points at energies beloew the peak in the cross section in the least-squares procedure. The coefficient data values are 1. The threshold energy for the reaction in eV. 2. The minimum energy , Emin, for application of the fitting formula. 3. The coefficient aBorn,2 used in DOC=H-HE-PLASMA. 4. The coefficient nBorn,2 used in DOC=H-HE-PLASMA. The Boolean labels include a. The cross section at Emin in label XS(EMIN) b. The fitting error, ERROR. (see DOC=H-HE-PLASMA for meaning).#EEXCH2IAEA subroutine to evaluate the analytic function used to represent the electron impact excitation cross section (cm[2]) as a function of electron energy (eV) of L. Vriens and H. M. Smeets. (See DOC=VRIENS,PRA,22,1980). The coefficient data values are 1. The principal quantum number (n) of the initial state. 2. The principal quantum number (m) of the final state. On the first call to the evaluation function the coefficient array is updated with the following data, where the integer refers to the index for the coefficient array , 3. The threshold energy (eV) for the transition. 4. Parameter A(n,m) as quoted in the formulae of DOC=H-HE-PLASMA. 5. Parameter B(n,m) as quoted in the formulae of DOC=H-HE-PLASMA. 6. Parameter delta(n,m) as quoted in the formulae of DOC=H-HE-PLASMA. 7. Parameter gamma(n,m) as quoted in the formulae of DOC=H-HE-PLASMA.#EIONH2IAEA subroutine to evaluate the analytic function used to represent the electron impact ionization cross section (cm[2]) as a function of electron energy (eV) of L. Vriens and H. M. Smeets. (See DOC=VRIENS,PRA,22,1980). The coefficient data values are 1. The principal quantum number (n) of the initial state. On the first call to the evaluation function the coefficient array is updated with the following data, where the integer refers to the index for the coefficient array , 2. The threshold energy (eV) for the transition.#EIONHRIAEA subroutine to evaluate the analytic function used to represent the electron impact ionization rate coefficient (cm[3]/s) as a function of plasma temperature (eV) of L. Vriens and H. M. Smeets. (See DOC=VRIENS,PRA,22,1980). The coefficient data values are 1. The principal quantum number (n) of the initial state.#JRREC1IAEA subroutine to evaluate the analytic function used to represent the radiative recombination rate coefficient (cm[3]/s), as a function of temperature (eV), into the n=1 and n=2 levels as used in reaction 2.1.8 of Janev et al. (See DOC=H-HE-PLASMA). The coefficient data values are 1. The threshold energy (eV) for the transition. 2. Parameter A(n,l) as quoted in the formulae of DOC=H-HE-PLASMA. 3. Parameter chi (n,l) as quoted in the formlulae of DOC=H-HE-PLASMA.#JRREC3IAEA subroutine to evaluate the analytic function used to represent the radiative recombination rate coefficient (cm[3]), as a function of temperature (eV), into the levels with principal quantum number, n>=3 as used in reaction 2.1.8 of Janev et al. (See DOC=H-HE-PLASMA). The evaluation function is only defined for the reactions covered in DOC=H-HE-PLAMSA. It may need to be modified to cover reactions not contained in the above database. The coefficient data values are 1. Z, the atomic number of the target. 2. Q, the charge state of the target ion. (for Hydrogen Q=1, for Helium Q= 1 or 2). 3. The principal quantum number (n) of the final state. (n>=3)#JEEXC1IAEA subroutine to evaluate the analytic function used to represent electron impact excitation cross sections (cm[2]), as a function of the electron energy (eV), for dipole allowed transitions The analytic expression is defined in reaction 2.3.7 in DOC=H-HE-PLASMA. The coefficient data values are 1. Parameter beta(n) as quoted in the formulae of DOC=H-HE-PLASMA. 2. Parameter gamma(n) as quoted in the formlulae of DOC=H-HE-PLASMA. 3. Parameter delta as quoted in the formlulae of DOC=H-HE-PLASMA. 4. Principal quantum number (n) of the initial state. 5. Orbital angular momentum (l) of the initial state. 6. Orbital angular momentum (l') of the final state. 7. Spin multiplicity (2S+1) of the initial state. 8. The total angular momentum of the initial state, only required for fine structure transitions. 9. Principal quantum number (n') of the final state. On the first call to the evaluation function the coefficient array is updated with the following data, where the integer refers to the index for the coefficient array, 10. The threshold energy (eV) for the transition. 11. The oscillator strength for the transition. 12. Coefficient xi of the formulae in DOC=H-HE-PLASMA.#JEEXC2IAEA subroutine to evaluate the analytic function used to represent electron impact excitation cross sections (cm[2]) as a function of the electron energy (eV), for dipole forbidden transitions without a spin change. The analytic expression is defined in reaction 2.3.2 in DOC=H-HE-PLASMA. The coefficient data values are 1. Parameter beta(n,l) as quoted in the formulae of DOC=H-HE-PLASMA. 2. Parameter gamma(n,l) quoted in the formlulae of DOC=H-HE-PLASMA. 3. Parameter delta(n,l) quoted in the formlulae of DOC=H-HE-PLASMA. 4. Principal quantum number (n) of the final state. 5. Orbital angular momentum (l) of the final state. On the first call to the evaluation function the coefficient array is updated with the following data, where the integer refers to the index for the coefficient array, 6. The threshold energy (eV) for the transition. 7. The energy of ionization of the (n,l) state. 8. Coefficient t of the formulae in DOC=H-HE-PLASMA.#JEEXC3IAEA subroutine to evaluate the analytic function used to represent electron impact excitation cross sections (cm[2]) as a function of electron energy (eV), for dipole forbidden transitions with a spin change. The analytic expression is defined in reaction 2.3.3 in DOC=H-HE-PLASMA. The coefficient data values are 1. Parameter a(n,l) as quoted in the formulae of DOC=H-HE-PLASMA. 2. Parameter b quoted in the formlulae of DOC=H-HE-PLASMA. 3. Principal quantum number (n) of the final state. 4. Orbital angular momentum (l) of the final state. On the first call to the evaluation function the coefficient array is updated with the following data, where the integer refers to the index for the coefficient array, 5. The threshold energy (eV) for the transition. 6. The energy of ionization of the (n,l) state.#JEEXC4IAEA subroutine to evaluate the analytic function used to represent electron impact excitation cross sections (cm[2]) as a function of electron energy (eV), for dipole forbidden transitions with a spin change. The analytic expression is defined in reaction 2.3.4 in DOC=H-HE-PLASMA. The coefficient data values are 1. Parameter Q as quoted in the formulae of DOC=H-HE-PLASMA. 2. Orbital angular momentum (l) of the initial state. 3. Spin multiplicity (2S+1) of the initial state. 4. Orbital angular momentum (l') of the final state. 5. Spin multiplicity (2S'+1) of the final state. On the first call to the evaluation function the coefficient array is updated with the following data, where the integer refers to the index for the coefficient array, 6. The threshold energy (eV) for the transition.#JEEXC6IAEA subroutine to evaluate the analytic function used to represent electron impact excitation cross sections as a function of electron energy (eV), for transitions between states defined in terms of the initial and final principal quantum numbers alone. The analytic expression is defined in reaction 2.3.8 in DOC=H-HE-PLASMA. The coefficient data values are 1. Principal quantum number (n) of the initial state. 2. Principal quantum number (m) of the final state. On the first call to the evaluation function the coefficient array is updated with the following data, where the integer refers to the index for the coefficient array, 3. The threshold energy (eV) for the transition. 4. Coefficient xi of the formulae in DOC=H-HE-PLASMA. 5. The oscillator strength for the transition, averaged over the orbital angular momentum and spin quantum numbers.#JIONHEIAEA subroutine to evaluate the analytic function used to represent the electron impact ionization cross sections (cm[2]), as a function of electron energy (eV), as defined in reaction 2.3.10 in DOC=H-HE-PLASMA. The coefficient data values are 1. An integer, ITRANS, that defines the type of initial excited target state. For (n,n') transitions ITRANS = 1 while for (nl,n'l') transitions ITRANS = 2. 2. Coefficient beta of the formulae in DOC=H-HE-PLASMA. 3. Principal quantum number (n) of the initial state. For (nl,n'l') transitions the following coefficients are also required, 4. Orbital angular momentum (l) of the initial state. 5. Spin multiplicity (2S+1) of the initial state. On the first call to the evaluation function the coefficient array is updated with the following data, where the integer refers to the index for the coefficient array, 6. The threshold energy (eV) for the transition.#JEION5IAEA subroutine to evaluate the analytic function used to represent the electron impact ionization cross sections (cm[2]), as a function of electron energy (eV), as defined in reaction 2.3.10 in DOC=H-HE-PLASMA. The coefficient data values are 1. Coefficient A of the formulae in DOC=H-HE-PLASMA. 2. Principal quantum number (n) of the initial state.#JRREC2IAEA subroutine to evaluate the analytic function used to represent the radiative recombination rate coefficient (cm[3]) as a function of temperature (eV), into n=2 levels as used in reaction 2.3.13 of Janev et al. (See DOC=H-HE-PLASMA). The coefficient data values are 1. The threshold energy (eV) for the transition. 2. Parameter A(n,l) as quoted in the formulae of DOC=H-HE-PLASMA.#JEEXC5IAEA subroutine to evaluate the analytic function used to represent electron impact excitation cross sections (cm[2]) as a function of electron energy (eV), for a target of singly charged Helium. The basic form of the analytic expression is given in reactions 2.3.17 in DOC=H-HE-PLASMA. The coefficient data values are 1. An integer, TTYPE. If TTYPE = 1 the expression contains a natural logarithm involving the reduced energy u, if TTYPE=2 this term is replaced by unity. 2. Normalisation coefficient (leading number in formula of DOC=H-HE-PLASMA) 3. Parameter C(n) as quoted in the formulae of DOC=H-HE-PLASMA. 4. Parameter phi(n) as quoted in the formulae of DOC=H-HE-PLASMA. 5. Principal quantum number (n) of the initial state. 6. Principal quantum number (m) of the final state. On the first call to the evaluation function the coefficient array is updated with the following data, where the integer refers to the index for the coefficient array, 7. The threshold energy (eV) for the transition.#JEEXC7IAEA subroutine to evaluate the analytic function used to represent electron impact excitation cross sections (cm[2]), as a function of electron energy (eV), for a target of singly charged Helium for (n,n') transitions with n>=3. The analytic expression is given in reactions 2.3.18 in DOC=H-HE-PLASMA. The coefficient data values are 1. Principal quantum number (n) of the initial state. 2. Principal quantum number (m) of the final state. On the first call to the evaluation function the coefficient array is updated with the following data, where the integer refers to the index for the coefficient array, 3. The threshold energy (eV) for the transition. 4. The oscillator strength for the transition, averaged over the orbital angular momentum and spin quantum numbers.#JDREC1IAEA subroutine to evaluate the analytic function used to represent the reaction rate coefficient for dielectronic recombination. (cm[3]) as a function of temperature (eV), The analytic expression is given in reactions 2.3.22 in DOC=H-HE-PLASMA. The coefficient data values are 1. Principal quantum number (n) of the initial state. 2. Principal quantum number (m) of the final state.#JAN3IAEA subroutine which is a general form for a cross section with simple logarithmic and inverse energy dependence, which is of the form cross section = A * LN ( B*E + C) / E where A, B and C are coefficients and E is the energy. For instance, as used in reaction 2.1.3 in DOC=H-HE-PLASMA. The coefficient data values are 1. Threshold energy for the transition (eV). 2. Coefficient A. 3. Coefficient B. 4. Coefficient C.#PEXCH4IAEA subroutine to evaluate the analytic function used to represent the heavy particle impact excitation cross sections (cm[2]), as a function of impact energy (eV), from the dipole close-coupling approximation. The analytic expression is defined in reaction 3.1.5 in DOC=H-HE-PLASMA. The coefficient data values are 1. The atomic number of the target. 2. The atomic number of the projectile. 3. An integer, ITRANS, that defines the type of initial excited target state. For (n,n') transitions ITRANS = 1 while for (nl,n'l') transitions ITRANS = 2. 4. Principal quantum number (n) of the initial state. 5. Principal quantum number (n') of the final state. For (nl,n'l') transitions the following coefficients are also required, 6. An integer, SUMEN, that defines the choice of ionization energy to be used. If SUMEN = 0 the ionization energy for the specific (n,l,s) state is taken, while if SUMEN = 1 the ionization energy averaged over possible (l,s) values is taken. 7. Orbital angular momentum (l) of the initial state. 8. Orbital angular momentum (l') of the final state. 9. Spin multiplicity (2S+1) of the initial state. 10. The total angular momentum of the initial state, only required for fine structure transitions. On the first call to the evaluation function the coefficient array is updated with the following data, where the integer refers to the index for the coefficient array, 11. The threshold energy (eV) for the transition. 12. The energy difference between the initial and final states. 13. Coefficient lambda(n,n').#IONBEAIAEA subroutine to evaluate the analytic function used to represent the heavy particle impact ionization cross sections (cm[2]), as a function of impact energy (eV), from the binary encounter approximation. The analytic expression is defined in reaction 3.1.7 in DOC=H-HE-PLASMA. The coefficient data values are 1. The atomic number of the target. 2. The charge of the core of the target. 3. The binding energy of the target atom. 4. The charge of the projectile ion. 5. An integer, ITRANS, that defines the type of initial excited target state. For (n,n') transitions ITRANS = 1 while for (nl,n'l') transitions ITRANS = 2. 6. An integer, SUMEN, that defines the choice of ionization energy to be used. If SUMEN = 0 the ionization energy for the specific (n,l,s) state is taken, while if SUMEN = 1 the ionization energy averaged over possible (l,s) values is taken. 7. Principal quantum number (n) of the initial state. 8. Orbital angular momentum (l) of the initial state. 9. Spin multiplicity (2S+1) of the initial state. On the first call to the evaluation function the coefficient array is updated with the following data, where the integer refers to the index for the coefficient array, 11. The threshold energy (eV) for the transition. 12. The "effective" value of the principal quantum number used in evaluating the cross section.#JCX1IAEA subroutine to evaluate the analytic function used to represent the charge transfer cross sections (cm[2]), as a function of impact energy (eV), from the analytic expression defined in reaction 3.1.12 in DOC=H-HE-PLASMA. The coefficient data values are 1. Principal quantum number (n) of the initial state.#JCX2IAEA subroutine to evaluate the analytic function used to represent the charge transfer cross sections (cm[2]), as a function of impact energy (eV), from the analytic expression defined in reaction 6.1.7 in DOC=H-HE-PLASMA. The coefficient data values are 1. Principal quantum number (n) of the initial state. 2. Principal quantum number (n') of the final state.#DBETAIAEA function to evaluate the D(beta) function defined in Janev R.K. and Presnyakov L.P., J. Phys. B., 13, 4233, (1980) and also as defined on page 319 in DOC=H-HE-PLASMA. The argument passed to DBETA is the value of beta.#OSCSTH1IAEA subroutine to return the oscillator stength, F(nl,ml') for hydrogen. The transitions included are those contained in table A.1 on page 314 in DOC=H-HE-PLASMA. The argument call of OSCSTH1 is, (in order): 1. Principal quantum number (n) of the initial state. (Integer input) 2. Orbital angular momentum (l) of the initial state. (Integer input) 3. Principal quantum number (m) of the initial state. (Integer input) 4. Orbital angular momentum (l') of the final state. (Integer input) 5. Oscillator strength, F(nl,ml'), (Real output)#OSCSTH2IAEA subroutine to evaluate the oscillator stength, F(n,m) for hydrogen from the formula of Johnson L. C., Astrophys. J. 174, 227, (1972). The argument call of OSCSTH2 is: 1. Principal quantum number (n) of the initial state. (Integer input) 2. Principal quantum number (m) of the final state. (Integer input) 3. Oscillator strength, F(n,m). (Real output)#OSCSTHEIAEA subroutine to return the oscillator stength, F(nl,ml') for helium. The transitions included are those contained in table A.6 on page 317 in DOC=H-HE-PLASMA. The argument call of OSCSTHE is, (in order): 1. Principal quantum number (n) of the initial state. (Integer input) 2. Orbital angular momentum (l) of the initial state. (Integer input) 3. Principal quantum number (m) of the initial state. (Integer input) 4. Orbital angular momentum (l') of the final state. (Integer input) 5. Spin multiplicity (2S+1) of the state. (Integer input) 6. The total angular momentum of the initial state, only required for fine structure transitions (Integer input). 7. Oscillator strength, F(nl,ml'), (Real output) 8. Error Indicator. (Character output field)#EXINTIAEA function to evaluate the exponential integral function defined in "Handbook of Mathematical Functions", Abramowitz M. and Stegun I. A, National Bureau of standard, Washington DC, page 231, 1964.#HEIONENIAEA subroutine to return the ionization energy for excited states of helium taken from tables A.4 and A.5 on page 316 in DOC=H-HE-PLASMA. The argument call of HEIONEN is 1. Principal quantum number (n) of the excited electron. (Integerinput) 2. Orbital angular momentum (l) of the excited electron. (Integer input) 3. Spin multiplicity (2S+1) of the state. (Integer input) 4. An input integer, SUMEN, that defines the choice of ionization energy to be returned. If SUMEN = 0 the ionization energy for the specific (n,l,s) state is taken, while if SUMEN = 1 the ionization energy averaged over possible (l,s) values is returned. 5. The ionization energy (eV). (Real output field). 6. Error Indicator. (Character output field)#HEEXCENIAEA subroutine to return the excitation energy for excited states of helium taken from tables A.4 and A.5 on page 316 in DOC=H-HE-PLASMA. The argument call of HEEXCEN is 1. Principal quantum number (n) of the excited electron. (Integer input) 2. Orbital angular momentum (l) of the excited electron. (Integer input) 3. Spin multiplicity (2S+1) of the state. (Integer input) 4. An input integer, SUMEN, that defines the choice of excitation energy to be returned. If SUMEN = 0 the excitation energy for the specific (n,l,s) state is taken, while if SUMEN = 1 the excitation energy averaged over possible (l,s) values is returned. 5. The ionization energy (eV). (Real output field). 6. Error Indicator. (Character output field)#FNAGEXIAEA subroutine to return the electron impact excitation reaction rate coefficients for a specified electron temperature in eV, using the analytic expressions of Y. Itikawa et al, See DOC=ADNDT,33,149,(1985). In this reference two type of analytic fits are employed to represent the reaction rate coefficients. The basic forms are a power-log type expression and an exponential type expression. See DOC=ADNDT,33,149,(1985) for details. The argument call of FNAGEX is 1. Index for the type of fit, either 1 or 2, Variable ITYPE 2. The excitation energy (eV) 3. The lower limit of the fit to the reduced energy, X 4. The upper limit of the fit to the reduced energy, X 5. statistical weight of the initial state (2s+1)*(2l+1) 6. Parameter A, 7. Parameter B, 8. Parameter C, 9 Parameter D 10. Parameter E If ITYPE = 1 this can be followed by three more parameters for the region which contains resonaces represented by a linear term 11. Parameter P, 12. Parameter Q, 13. Parameter X1 If ITYPE = 2 this can be followed by 11. Parameter F and possibly three more parameters for the region which contains resonances represented by a linear term. 12. Parameter P, 13. Parameter Q, 14. Parameter X1#REFL1IAEA subroutine to return the general form for the particle or energy reflection coefficient as a function of projectile energy in keV. DOC=IAEA-AM-REP-249,(1991) The coefficient array passed to REFL1 is the six fitting coefficients needed to determine the reflection coefficient.#REFL2IAEA subroutine to return the particle or energy reflection coefficient as a function of projectile energy in keV for the extended energy fit. DOC=IAEA-AM-REP-249,(1991) The coefficient array passed to REFL2 is the eight fitting coefficients needed to determine the reflection coefficient.#REFL3IAEA subroutine to return the self particle or energy reflection coefficient as a function of projectile energy in keV. DOC=IAEA-AM-REP-249,(1991) The coefficient array passed to REFL3 is the nine fitting coefficients needed to determine the reflection coefficient.