Methods of Calculating Ionization Energies of Multielectron (Five or More) Isoelectronic Atomic Ions

We have previously used simple empirical equations to reproduce the literature values of the ionization energies of isoelectronic sequences of up to four electrons which gave very good agreement. We reproduce here a kinetic energy expression with corrections for relativity and Lamb shift effects which give excellent agreement with the literature values. These equations become more complex as the number of electrons in the system increases. Alternative simple quadratic expressions for calculating ionization energies of multielectron ions are discussed. A set of coefficients when substituted into a simple expression produces very good agreement with the literature values. Our work shows that Slater's rules are not appropriate for predicting trends or screening constants. This work provides very strong evidence that ionization energies are not functions of complete squares, and when calculating ionization energies electron transition/relaxation has to be taken into account. We demonstrate clearly that for particular isoelectronic sequences, the ionizing electrons may occupy different orbitals and in such cases more than one set of constants are needed to calculate the ionization energies.


Introduction
We have previously proposed a simple empirical equation to reproduce the literature values of the ionization energies of one-electron [1] and two-electron [2] atomic ions with very good agreement. This was recently extended to calculate ionization energies and first electron affinities of three and four electron ions [3] which also gave very good agreement with the literature values. However, we used a potential energy approach in our equation which has no theoretical basis.
With the development of quantum theory, the twoparticle problem can be solved exactly, and the kinetic energy of the electron in a hydrogen atom or hydrogen-like ion can be calculated using the Schrödinger equation. In 1930, Dirac [4] produced an equation which included a relativistic correction for the electron energy levels. Lamb and Retherford [5][6][7] were able to show that there is a small shift (now known as the Lamb shift) in the energy levels of the hydrogen atom not included in the Dirac equation. When the Lamb shift is taken into account, the actual ionization energies of one electron atoms are slightly less than that calculated by the Dirac equation. The calculation of the Lamb shift now forms an essential part in the theoretical calculations of energy levels and ionization energies of one and two electron atomic ions [8,9].
In this work, we begin by showing the expressions with relativistic corrections to calculate ionization energies which give excellent agreement with the literature values. However, as the number of electrons in the sequence increases, the expressions become more and more complex and less predictable. Hence, we discuss alternative simple equations to calculate energies of isoelectronic sequences in excess of five electrons and show that it is more practical to adopt a simple expression to reproduce ionization energies which still provide very good agreement with generally accepted values. To maintain our aim of simplicity, the expressions/equations only contain fundamental constants or values derived from fundamental constants and fairly simple numbers.

Sources of Data
Moore [10][11][12][13] provided very detailed tables of atomic energy levels and ionization potentials in wave numbers (cm −1 ) and values converted from wave numbers to electron volts 2 The Scientific World Journal (eV) (where 1 cm −1 equals 1.2398418 × 10 −4 eV) for atoms and atomic ions with estimated experimental errors and references to original work. These remain the most extensive survey of ionization energies and are still quoted in recent publications. The CRC Handbook of Chemistry and Physics [14] contains comprehensive data from Moore and later sources but does not supply information on estimated uncertainties. The majority of published ionization energy data are now available on the National Institute of Standards and Technology web site (http://www.nist.gov/srd/). These compilations include values of ionization energies that are accurately measured as well as crude estimates.

Ionization Energies of the Hydrogen to Boron Isoelectronic Sequences
When an electron in a multielectron system is ionized, we assume that the ionization energy contains ( ) the kinetic energy term; ( ) the Lamb shift term, ( ) the relaxation/ transition energy term, and any residual interaction ( ) [2]. Therefore, the ionization energy can be represented as [15] where is the principal quantum number and is the reduced mass and are provided in Table 2. , , and represent , , and multiplied by 1/ 2 , respectively. , , and are multiplied by 1/ 2 because the electron ionized occurs in energy level and it revolves around the centre of mass. A list of reduced masses are given in Table 2. It is common to present ionization energies in eV (electron volt). Calculated results in this work are converted to eV from Joules by using the relationship of 1 eV = 1.60217648 × 10 −19 J, any other figures in cm −1 are converted to eV by multiplying them with the value 0.00012398418.

The Relativistic, Lamb Shift, Electron Transition/Relaxation, and Residual Corrections
The energy of an electron moving in a Bohr orbit can be represented by where is the electron rest mass, V is the velocity of the electron, 1 2 stand for the charges of the electron and nucleus, is the permittivity of a vacuum, and 0 is the Bohr radius. The velocity V of the electron in the hydrogen atom can be calculated from the above relationship and is equal to 2.1876913 × 10 6 m/sec. The velocity V of the electron in successive atoms of the one-electron series increases by times where is the atomic/proton number or V = V . When there is more than one electron in the system, we assume that the velocity of the electron changes by ( − ) where is the screening constant for that electron.
The theory of relativity [16] points out that the mass of a moving particle is given by the expression is the rest mass of particle. Expansion of this expression gives therefore it follows that (1/2) The kinetic energy components including relativistic correction for a one-electron atom at an equilibrium position (when the relativistic correction is half that of the maximum) are then The ionization energy of a one-electron atom is then The Lamb shift is usually computed by highly complex formulas which require lengthy computer routines to compute [17]. We assume (without theoretical justification) that the Lamb shift is a relativistic charge, mass, and size ratio effect. The reduced mass calculation [18] implicitly assumes that the electron and nucleus are point charges. But they have a finite size hence there needs to be an extra component in the reduced mass calculation. The energy reduction to take account of reduced mass is not simply /( + )((1/2) V 2 ) but should include a function of the charge and the ratio of nuclear to atomic size, and this component is calculated by the following expression: where is the electron mass, is the proton mass, and the factor for the reduced mass correction for hydrogen is /( + ). The size of the nucleus increases roughly proportional to 1/3 [19]. is the fine structure constant, and ( /2 .67 ) is a crude approximation of the square root of the ratio of nuclear to atomic size for hydrogen, and is the mass number of the atom. In a one electron system, is zero. After an electron is ionized, one or more of the remaining electron(s) is/are attracted more closely to the nucleus. The attractive energy between the proton and the remaining electron(s) changes because of the change in screening experienced by the remaining electron(s) before and after ionization, [15] and this transition/relaxation energy is a function of 2

4
( where 1 is the screening constant for the remaining electron(s) after ionization and 2 is the screening constant for the remaining electron(s) before ionization.
The Scientific World Journal 3 In the helium system, there are two electrons and both occupy the 1s orbital. Since each electron occupies half of the space and each is repelled by only one other electron (here we have assumed that the two electrons act as in a two-particle problem and are equivalent), the screening constant is a half (0.5). After ionization, there is zero screening.
In the lithium series, the electron that is ionized occupies a higher (2s) orbital and is shielded by two 1s electrons, and the screening increases by 1 to 1.5. The extra screening experienced by the two inner electrons in the 1s orbital increases to 0.625 and which is an increase by 1/8 or 0.125 rather than 0.5 because the third electron occupies a different orbital and in a different electron shell (i.e., 0.5 of 0.5 of 0.5). After ionization, only two electrons are left in the system and the screening reduces to just 0.5.
In the beryllium series, the electron that is ionized occupies the 2s orbital. Since there are four electrons and each moving in an elliptical orbit, each may at any time interact differently with the other three electrons. The screening of the fourth electron increases by a half to 2, and the other electron in the 2s orbital experiences a screening of between 1.5 and 2 and is 1.75. After ionization of the fourth electron, the screening experienced by the third electron drops back to 1.5.
For the boron system, the outermost electron occupies a new orbital and the screening increases by 1 to 3. The screening of the fourth electron increases by 0.25 to 2.25, and after ionization the screening reduces back to 2.
We assume that there are two opposite and competing residual electron-electron interactions, and this is discussed in detail in previous work and is not fully reproduced here [15]. In summary, there are temporary asymmetric distributions of electrons (e.g., they may be nearer to each other or further apart than average). The first type (when they are nearer to each other) is residual electron-electron repulsion which reduces slightly the energy required to ionize the electron, and this reduction diminishes very rapidly because with each successive ionization the size of the electron orbit/shell becomes smaller and the electrons become much more tightly bound to the nucleus. The reduction in energy resulting from this interaction is expressed by The opposite and competing electron-electron interaction occurs when there are more than two electrons in the system. In a three electron system, two electrons occupy the 1s orbital and one occupies the 2s orbital. The instantaneous asymmetric distribution of electrons produces an effect which result in the 2s electron being screened slightly less than 0.5 from each of the 1s electrons. The result is that one of the electrons is temporarily attracted to the nucleus a bit more than expected and hence increases the amount of energy required to remove it from the atom/ion. This is represented by and the total residual electron interaction energy change is Symbols and values of constants shown in all the expressions in this work are given in Table 1.
is the number of electron-electron interactions before ionization and is the number of electron-electron interactions after ionization.

Ionization Energies of One and Two Electron Ions
There is only one electron and is 1, the formula for calculating the ionization energy is 4 The Scientific World Journal The one-electron ionization energies calculated by (11) when compared with the ionization energies published in the CRC Handbook of Chemistry and Physics agree to 99.999% or better in the majority of cases. The biggest absolute difference is 0.086 eV (to 3 decimal places) from a calculated value of 5469.95 eV or 0.00164% [15]. The ionization energy of a two electron system is ( Since, for simplicity, we have not applied a relativistic correction to we have made a crude approximation of reducing the relativistic correction in by another 5% to 0.45 (rather than 0.5) and where V is (V ( − 0.5)), so it is When compared to the CRC Handbook of Chemistry and Physics the majority of values differ by less than 0.01%, the largest absolute difference is 0.216 eV from a calculated value of 2437.846 eV or 0.0089%.

Ionization Energies of Three-, Four-, and Five-Electron Ions
The ionization energy of a three-, four-, and five electron system is The Scientific World Journal As with the two-electron system, we have not applied a relativistic correction to but we have made a crude approximation of reducing the relativistic correction in by 5% to 0.45 and where V is (V ( − )), so that the expression becomes The agreement with the literature values of ionization energies is 99% or better in all cases. It is evident from the above discussion that as the number of electrons in the system increases the equations get more complicated and it is difficult to determine the various corrections.

Alternative Simple Equations to Calculate Ionization Energies
Besides the more complicated equations shown above, ionization energies can be calculated with simpler expressions but not so precise. A simple formula is sometimes used to calculate ionization energies. It often takes the form of [20] = ( 1 where is the Rydberg constant for hydrogen ( ∞ , the Rydberg constant for infinite mass is equivalent to 13.6059 eV), * is an "effective" quantum number, is the atomic number, and is the screening constant based on Slater's rules [21] which enable approximations of analytic wave functions to be constructed for rough estimates [22]. Equation (17) makes a very simplistic assumption that when an electron is removed from an atom or ion, the atom/ion remains unchanged except for the removal of that electron. We believe that it is incorrect to use equations like (17) (where the ionization energy is considered as a function of a complete square) to calculate energies in isoelectronic series. For multielectron systems, we need to consider electron transition/relaxation and other smaller components which may influence the ionization energy of an electron. In a multielectron system, the main components of the energy change during ionization are the electron-proton energy or ( − ) 2 and the electron relaxation energy (or where is a constant dependent on the particular shell/orbital, 1 and 2 represent the screening constants of the remaining electrons after and before the electron is ionized). In most cases, the electron-proton energy component alone accounts for 90% or more of the energy change and together with the relaxation energy can represent more than 95% of the total energy change. If factors which in total contribute only a small percentage of the energy change, such as residual repulsion, pairing or exchange energies are excluded, the expression for calculating the ionization energy can be approximated to This can be expanded to become Since in the second half of expression (19) the 2 term cancels out, only 2 ( 1 − 2 ) and can be reduced to , and ( 2 2 − 1 2 ) becomes a constant . This can be rearranged and simplified to where and are constants for each isoelectronic system, is the principal quantum number, and, unlike (17), expression (20) is not a complete square. We have formulated the following rules for working out screening (shielding) constants: (1) for each additional electron in the system, the screening increases by 0.5 unless (2) the electron to be ionized occupies a new orbital such as from beryllium to boron when it increases by 1 and (3) to account for the pairing effect, such as from nitrogen to oxygen, the screening constant increases by 1. For example, for the carbon system, the screening increases by 0.5 units to 3.5 and increases by a further 0.5 to 4 for nitrogen but increases to 5 for oxygen.
We made different estimates of the screening constants 2 and 1 and obtained various values of and . We then selected the values that give the best results and a list of screening constants, coefficients and and reduced masses as shown in Table 2. They are used with equation (20) to calculated the ionization energies of isoelectronic sequences from five electrons.
The values calculated from expression (20) are presented for the first six appropriate members of each series for the following two reasons. Firstly, as we have already shown [14], ionization energies of the first few members of isoelectronic sequences are the most precise. Sometimes, only the first four or five members of a series are experimentally measured and uncertainties increase further along a series. 6 The Scientific World Journal

Ionization Energies from Five-to Eighteen-Electron Isoelectronic Atomic Ions
Ionization energies reported in the CRC Handbook for five to eighteen electronic series (first six members) are given in Table 3. Values of ionization energies calculated using our coefficients are provided in Table 4. Percentage differences between our values and values published by The Scientific World Journal 7 the CRC Handbook as listed in Table 5 show that all values agree to 98% or better. Just over 76%, the calculated values agree to 99% or better.

Ionization Energies from Nineteen-to Twenty-Seven-Electron Isoelectronic Sequences
Treatment of ionization energies of isoelectronic with nineteen electrons or more are different and more complicated. This is because, from atomic number nineteen, the electron to be ionized is in period 4 of the periodic table and such an isoelectronic system includes ionization of s or d electrons [23]. We have demonstrated that for the transition metals or lanthanide [24] elements, the ionization process is complicated and, in the majority of cases, d or f electrons are not removed in the first or second ionization. For example, consider the potassium, calcium, scandium, and titanium series (isoelectronic series with 19, 20, 21, and 22 electrons resp.). For the potassium series (19-electron isoelectronic), the first ionization of potassium and the second ionization of calcium both involve removal of a 4s electron. But the third ionization of scandium and the fourth ionization of titanium (which are in the same isoelectronic series) involve removal of a 3d electron. For the calcium series (20-electron isoelectronic), the first ionization of calcium and the second ionization of scandium involve removal of a 4s electron whereas the third ionization of titanium involves removal of a 3d electron.
Similarly with the scandium and titanium series, for the first two members of the series, the energy change is the energy required to ionize a 4s electron but from the third member of the series onwards the energy change is the energy for ionizing a 3d electron (a more detailed discussion is provided in a previous work [23]). Equation (20) shows that ionization is a function of 1/ 2 . Therefore, for the 19, 20, 21, and 22 isoelectronic series, the correct fraction or decimal for 1/ 2 to use is 1/16 (or 0.0625) for the first two members of the series but 1/9 (or 0.11111) for the remainder of the series. Hence, for the above reasons, it is incorrect to use a single set of coefficients to calculate the energies of an isoelectronic sequence which may involve removal of electrons in different orbitals. We believe that there is little value calculating the first and second ionization energies of sequences beginning with the potassium. From the potassium series onwards, the third ionization energy requires a different set of coefficients. For example, for the potassium series, a set of coefficients is used to calculate the third ionization energy of potassium, the fourth ionization energy of calcium, and the fifth ionization energy of scandium, and so on because a 3p electron is ionized in all cases. For the scandium series, the third ionization of scandium, the fourth ionization of titanium, and the fifth ionization of vanadium and so on can be calculated using one set of coefficients since in all cases a 3d electron is ionized. However, the third ionization energy of the potassium series and the third ionization energy of the calcium series are identical to the fifth ionization of the chlorine and argon series, respectively. Hence, these two series will not be repeated in Tables 6 and 7.
We have calculated the ionization energies up to the cobalt series (which contains five appropriate published values for comparison) because beyond the cobalt series there are fewer and fewer published values available for use in comparison. Ionization energies of isoelectronic series reported in the CRC Handbook for sequences from scandium to cobalt (for each series beginning with the third ionization energy) are given in Table 6. Values of ionization energies calculated using our coefficients are provided in Table 7. Percentage differences between our values and values in the CRC Handbook as listed in Table 8 show that all calculated values agree to 98%, or better and just under 83% of the values agree to 99% or better.

8
The Scientific World Journal

Discussion
We have used a simple quadratic expression in this work. We have not considered exchange and orbital energies (20) and have ignored any residual interactions or relativistic corrections, which for multielectron systems are difficult to apply. Hence, it is not surprising that the agreement with some of the generally accepted values is less than 99%. However, some of the differences between the calculated values and CRC Handbook values are less than the experimental uncertainties. However, as we have shown above, equations for solving ionization energies can be very complicated and the results may be unpredictable as the number of electrons in an isoelectronic series increases. Therefore, we believe that there is a strong case to use a simple quadratic expression rather than trying to create complex equations to calculate ionization energies. Although Slater's rules are still cited in recent publications [25] as adequate for predicting most periodic trends, it has been pointed out that the rules are unreliable when orbitals with a total quantum number of 4 [26] is reached (e.g., a 3p orbital has a principal quantum number of 3, orbital quantum number of 1, and magnetic quantum number of 1, and spin quantum number of 1/2 already has a total quantum number of 5). Equation (17) and Slater's rules are based on simple assumptions and are unable to account for many different features of ionization energies across the periodic table. We Table 10: Ionization energies (eV) calculated using coefficients in Table 9.  have also shown that ionization energies are not functions of simple complete squares [23], and Slater's rules cannot account for the complex patterns in ionization energies shown in our previous work [24].

Conclusion
Ionization energies calculated by a kinetic energy approach with simple relativistic and Lamb shift corrections give remarkable agreement with generally accepted values for one-to five-electron isoelectronic series. However, for multielectron isoelectronic series with five or more electrons, there is no accepted methodology for calculating relativistic corrections. Therefore, it is practical and more manageable to use a simple quadratic expression to calculate ionization energies which still give very good agreement with generally accepted values. We have not attempted to calculate ionization energies of sequences beyond cobalt because there are few sequences where a long series of data are available. We believe that expression (20) can be used to calculate ionization energies of any multielectron isoelectronic series assuming that sensible coefficients are applied. It is evident that an equation based on (17) (or function of a complete square such as ( − ) 2 ) is not the correct representation of the energy change when an electron is ionized. We have also demonstrated that it is incorrect to use a single set of constants/coefficients to calculate the ionization energies of some isoelectronic series such as the potassium, calcium, or scandium series. Since a great many of the experimental measurements of ionization energies were done over half a century ago, and some of the published values are extrapolated/interpolated, or estimates with fairly large uncertainties there is also a strong case for new measurements to be undertaken and reliability of some of the currently accepted values to be reexamined.

10
The Scientific World Journal

Appendix
We believe it is useful to compare our results with other calculated results (using an equation which is a complete square) to show that our simple model of electron ionization is more realistic and reliable. Agmon [20] published a list of coefficients for the following equation to calculate ionization energies of isoelectronic series: where is the Rydberg constant for hydrogen. is the atomic number, is the screening constant, and * is an "effective" quantum number.
The list is shown in Table 9. Ionization energies (for five-to eighteen-electron sequences) calculated by the above equation with the set of coefficients are shown in Table 10. Percentage differences with those in the CRC Handbook are shown in Table 11.
We have shown that, in general, ionization energies of the first few members (usually the first two or three) of an isoelectronic sequence are most accurately determined, uncertainties normally increase across the higher members of a series. It is clear from Table 10 that the biggest differences between the calculated and accepted values occur where the accuracy and reliability of the accepted values are the greatest. Difference between the calculated and generally accepted values for some of ionization energies are greater than 10%. Secondly, only 13.1% of the calculated figures agree with the accepted ones to 99% or better (as compared to a much higher percentage of ours) and agreement of less (or worse) than 95% occur with over 29.7% of the calculated figures (whereas all our results agree to 95% or better). A major factor is that values of ionization energies are not functions of complete squares, and the equation used by Agmon does not take account of the electron relaxation/transition energy.