## Brazilian Journal of Physics

*versão impressa* ISSN 0103-9733

### Braz. J. Phys. v.35 n.1 São Paulo mar. 2005

#### http://dx.doi.org/10.1590/S0103-97332005000100008

**REGULAR ARTICLES**

**An efficient screening approach to be used in plasma modeling and ion-surface collision experiments **

**J. Pomarico; D. I. Iriarte; H. O. Di Rocco **

Instituto de Física Arroyo Seco Facultad de Ciencias Exactas, UNCPBA Pinto 399, (7000) Tandil, Buenos Aires, Argentina

**ABSTRACT**

In this work we show that the Layzer theory for atomic calculations provides a theoretical framework and also a powerful computational approach if correct rules for the calculation of the screening parameters are given. Using the virial as a model for potential energy and splitting of two-body operators as sum of onebody operators, a neat definition of screening is given, satisfying diverse physically indispensable properties. Many different experimental and theoretical results are reproduced with high accuracy, with no fitting procedure involving energy levels or numerical potentials. A *C*++ code and an executable file are available upon request.

**1 Introduction**

The use of screened hydrogenic models (SHM) for modeling plasmas and surface ion collision experiments has been well documented in recent works [1], [2]. Different methods are available for systematic numerical calculations of atomic parameters, based on the self-consistent Hartree-Fock procedure [3] [4]. However, it is important to have an efficient analytical approach, suitable for the calculation of energies of a large number of excited states. The purpose of this article is to give explicit formulas for the average energies of arbitrary complex configurations, specially of medium to highly ionized atoms. However, also good results for neutral and almost neutral atoms are obtained.

The starting point is the Layzer's formulation of the Z-dependent theory of the many-electron atom truncated to second order [4][5] [6]. Many authors pointed out in the past that the screening approximation is not very useful; however, that statement is not right provided that correct rules are given for calculating the screening. Moreover, relativistic corrections to the energy can be easily made and configuration interaction effects can be properly considered. It is important to point out that several published screening rules produce negative ionization potentials (IP) for neutral or few ionized systems [2][7]. Other authors obtain the screening parameters either in order to reproduce the IP's calculated by Hartree-Fock [9] [10] or by adjusting analytical potentials with numerical results [11].

In this work, no empirical adjustments are made nor negative IP's emerge from our calculations. For the screening parameters, we use the approach due to Kregar [12] generalized by one of us [13], based on the splitting of the two-body potential energy operators into the sum of effective one-body operators. Total binding energies, ionization potentials and energies of complex excited configurations are in good or very good agreement with both, theoretical and experimental values.

In the following, we use the atomic units, where the electron charge and mass and the Dirac constant are equal to one (*e = m _{e}* = = 1) and the energies are measured in Hartrees: 1 Ht º 27.21

*eV*. A numerical code written in

*C*++ and an executable file are available upon request.

**2 A Brief Summary of the Theory**

**2.1 The Z^{-1} expansion**

Layzer's formulation of the Z-dependent theory of the many-electron atom is the starting point of our work [4], [5] [6]. The approximate non-relativistic Hamiltonian in atomic units given by

Introducing the new variable r = *Zr*, it can be shown thatthe total *non-relativistic* energy can be written as the expansion [4]

where, *exactly*

and

*E*_{1} is given by the sum of the average Coulomb energy for electron pairs {*ij*} (see below, Equation 10), and therefore in terms of the Slater's integrals *F ^{k}* and

*G*, evaluated with hydrogenic wavefunctions with

^{k}*Z*= 1 [3][4]. For non-equivalent orbitals:

whereas for the equivalent ones:

the coefficients *g _{k}* and

*f*can be found in the above cited references.

_{k} If we restrict the expansion (2) up to *E*_{2}, the energy can be written as

where *w _{nl}* is the number of electrons in the

*n,l*shell and s

*is the corresponding screening parameter. Comparing Equations (2) and (7), we find that the s*

_{nl}*'s satisfy*

_{nl}

and

Then, to second-order approximation in the non-relativistic context, the average energy of a configuration is given by

where *w _{i}* is a short notation for , the number of electrons in the

*n*shell.

_{i},l_{i} Comparing Equation 8 with Equation 10, we see that *E*_{1} can be calculated by two ways: in terms of the Slater integrals *F ^{k}* and

*G*or in terms of the ss. Indeed, if we are capable of select a good method for s's calculations, the two values must be nearly equal. One of the merits of this work is the algorithm used for this calculus.

^{k} **2.2 The calculation and desirable properties of the screening parameters**

The concept of screening (and screened orbitals) is of old data and it is impossible to give a short account in this paper. In the past, screening parameters were obtained using spectroscopic data, numerical calculations and theoretical approaches. A short review can be found in the paper from one of the authors [6].

Screening parameters must have several indispensable properties: i) it must be derived theoretically in a neat form, ii) it must be dependent on *Z* and *N*, accounting for the contraction of the orbitals, iii) the sum given by Equation 8 must be *very* nearly to the value å* _{i}w_{i}*(

*w*– 1){

_{i }*ii*} + å

*{*

_{i,j}w_{i}w_{j}*ij*} ; this property is verified and it is exact for

*Z*® ¥, when s ® s

_{0}(see below, Eq.11), iv) the corresponding values of

*E*

_{2}must follow a law of the type

*E*

_{2}~

*a+b/Z*in order to reproduce the non-relativistic calculations using the Cowan's code (see below, section Results), v) it must satisfy the

*Z*expansion

We use the point of view introduced by Kregar generalized by one of us, based on the virial as potential energy and the splitting of the two-body potential energy operators into the sum of effective one-body operators (see below Equation12) [12][13]. Briefly, in the zero approximation, disregarding for a moment the exchange and sub-shell corrections, as can be viewed in Equations 5 and 6, the average Coulomb energy of electron pairs is (the *i-th* electron is equally or more strongly bound that the *j-th* one)

so that the operator 1/*r _{ij }*can be written as

where *g _{ij}* and

*f*are the external and internal screening parameters, respectively. When

_{ji}*i = j*,

*f*º

_{ii}= g_{ii}*k*. For a given value of

_{ii}*r*, inner screening is associated with electrons at

*r*

^{'}__<__

*r*, outer screening is due to a shell of charge located at

*r*>

^{'}*r*,adding a constant to the potential

*V*(

*r*). As it was shown in the References [12] and [13], the two parameters can be expressed in simple form in terms of

*y = Z*;

_{i}n_{j}/Z_{j}n_{i}*Z*and

_{i}*Z*are the screened charge seen by each electron. A neat demonstration of this result can be found in reference [13].

_{j}To add the exchange and sub-shell corrections implies the changes

and can be accounted for as the ratio between the second term of the r.h.s. of Equations 5 and 6 and *F*^{0}.

Then, taking into account the exchange interaction and subshell effects, effective charges are given by

where *q _{i}, q_{j}* are the occupation number of the respective subshells. As (see details in [13]), in terms of

*y = Z*

_{i}n_{j}/Z_{j}n_{i}

a short and simple iterative procedure must be used for their determination. The iteration can be reduced to one or two steps if, as zero approximation, Slater-like rules are used as starting point for *Z _{i}* (see below, section 3). Explicit expressions for the coefficients

*a*and

_{k}*b*can be found in reference [13].

_{k} Both the Layzer and the Kregar formulations ensures exact energy levels for H-like ions in the non-relativistic regime; for He-like ions, s_{0} = 5/16, as deduced for the variational principle.

**2.3 Relativistic corrections**

From the Breit-Pauli Hamiltonian the relativistic shift operator *H _{RS}* has five terms [4]:

that signifies respectively: mass correction, one- and two-body Darwin, spin-spin contact and orbit-orbit terms. Far, the more important are the two first ones, that scale as *Z*^{4}. These energy shifts can be estimated very simple and correctly as the sum of contributions

The remaining terms scale as *Z*^{3} and will be ignored in this work. An analysis is presented in the work by Safronova et al [2].

**3 Results and comparison with experiments and other calculations**

The comparison with experiment, in order to see the *Z*-dependence of experimental binding energies it is difficult, because there are relativistic effects dependent of *Z*^{4} and *Z*^{3} [8] and higher powers of *Z*. Therefore, we compare in first place our expansion with non-relativistic calculations. We verified that the non-relativistic values from the Cowan's code follow a dependence of the type

moreover, the Cowan values with relativistic corrections are very similar to the results from the Dirac-Hartree-Fock values from Huang et al. [15]

In Table 1, we show the center-of-gravity binding energies, including relativistic corrections for ground configurations, for atoms with closed (sub)shells. As we can see, our values are in very good agreement both with Hartree-Fock ([3], [15]) and experimental ones (up to *Z* = 20, because there are not experimental results for all elements and ionization degrees).

In Fig.1 we show the first ionization potential for the above cited elements, jointly with experimental values. The relation is of the order of two and better values can not be expected. In fact, the Layzer theory is valid, in principle, for large *Z*. Therefore this, the most rigorous test for our choice of the screening evaluation, is very promissory.

In Table 2 are shown the ionization energies for the *Ar* isonuclear sequence; the relative values are near one.

In Table 3, are shown the ionization energies of the *Ar* isoelectronic sequence; as it is expected, the relation between our values and the Cowan values are in excellent agreement.

In Fig.2 shows the binding energies of the internal sub-shells for *Xe*; our values are compared with those from the theory and the experiment [15]; as can be seen, the agreement is very good.

As other example of complex transitions, we show in Table 4 transitions of the type F_{0}F_{1} – 2*p*^{–1}3*d*F_{0}F_{1} as defined in Reference [2]: F_{0} = 1*s*^{2}2*s*^{2}2*p*^{6}, F_{1} = 3*d ^{n}* F

_{0}. Both values are within 1%, but it must be mentioned that the approach of Safronova et al. produces very bad values (negatives!) of the ionization energies for neutral and almost-neutral elements.

**3.1 Calculation of terms and energy levels**

The results shown up to this point allow the calculation of the configurations average energies [3]. This may be enough for calculating, for example, the array oscillator strength, which is frequently used for statistical treatments of dielectronic recombination and collisional excitation, in plasma physics and astrophysical applications. For the calculation of the multiplet oscillator strength, the *terms* are needed and, consequently, the Slater integrals, *F ^{k}*(

*ij*) and

*G*(

^{k}*ij*), which scale with

*Z*. Thus,

_{eff}*F*(

^{k}*ij;Z*) = (

*Z*– s

*) ×*

_{nl}*F*(

^{k}*ij;Z*= 1) ,

*G*(

^{k}*ij;Z*) = (

*Z*– s

*) ×*

_{nl}*G*(

^{k}*ij;Z*= 1). Values of these integrals for

*Z*= 1 can be found in the book by Condon-Odabasi [14]. In particular, and as an example,

*G*

^{1}(4

*s*4

*p;Z*= 1) = 0.02444565

*Ry*.

Table 5 shows the values of *G*^{1}(4*s*4*p;Z*) = 0.02444565(*Z* – s* _{nl}*) for the configuration 4

*s*4

*p*of

*Zn – like*with

*Z*= 30 – 74. These results are compared with the relativistic Hartree-Fock approximation (HXR). It is very important take into account that in practice, in order to give a better accordance with experimental spectroscopy, the theoretical Slater integrals must be multiplied by ''fudge'' factors lesser than 1 [3].

**3.2 Slater-like rules and constant values for the g_{ij}'s and the f_{ji}'s for Z ® ¥.**

As stated above (see equation 11), a desirable property for s is that for *Z *® ¥, s ® s_{0}. This is verified by our approach and implies that for *Z* ® ¥, the values of *f _{ij}* and

*g*trend themselves to constant values, which we can denominate and . An example is given in Table 6 for the particular case of the orbitals of

_{ij}*Kr*-like ions. These values, once known, can be stored and reliably reused either in calculations for high ionization stagesor, as input parameter to reduce the number of iterations needed to achieve convergency in Eq. 14. For example, using Eq. 13, s

_{1s}= 0.4973, s

_{2s}= 5.0253, etc.

**4 Conclusions**

In this work we showed that the Layzer expansion is not only a theoretical framework but also a very convenient quantitative method when a correct approach for screening parameter calculations is given. We select the point of view introduced by Kregar, generalized by one of us, taking into account subshell and exchange corrections. With our approach, we always obtain positive ionization potentials, even for neutral atoms, the most difficult test. This fact is easily understood because the Layzer theory is an expansion valid for large *Z*. In practice, for slightly ionized atoms (as low as three times), results are very accurate accordingly with the main scope of this work: plasma modelling and collision experiments. It it important to point out that in the present work no fitting procedure is used although, if we want to extend the method for neutral and almost neutral atoms, a correction factor (very close to 1) for the s^{'}s could be used. Moreover, it should be mentioned that, although a set of *f _{ji}* and

*g*

_{ij}*constants*could be used (which we have called and ), better results are obtained if the

*f*and

_{ji}*g*are handled as

_{ij}*parameters*. Our approach assure this property taking into account the variation with

*Z*and

*N*. However, it is true that for sufficiently ionized atoms, the use of constants provides a very good approximation, useful for plasma modeling and ion-surface collision experiments [2].

**References**

[1] J.G. Rubiano, R. Rodríguez, J.M. Gil, F. H. Ruano, P. Maretel, and E. Mínguez, J. Quant. Spectrosc. Radiat. Transfer, **72**, 575 (2002). [ Links ]

[2] U. I. Safronova et al., Phys. Scripta **47**, 364 (1993); [ Links ]M.W. Clark et al, Phys. Rev. A **47**, 3983 (1993). [ Links ]

[3] R. D. Cowan, *The Theory of Atomic Structure and Spectra* (University of California Press, Berkeley and Los Angeles, 1981). [ Links ]

[4] C. Froese Fischer, T. Brage, and P. Jonsson, *Computational Atomic Structure* (IoP Publishing, Bristol, 1997). [ Links ]

[5] D. Layzer, Ann. Phys. **8**, 271 (1959); [ Links ]D. Layzer et al., Ann. Phys. **29**, 101 (1964). [ Links ]

[6] H.O. Di Rocco, Il Nuovo Cimento D **20**, 131 (1998). [ Links ]

[7] G. Faussurier, C. Blancard, and A. Decoster, J. Quant. Spectrosc. Radiat. Transfer **58**, 233 (1997). [ Links ]

[8] D. Layzer, Int. J. Quantum Chemistry **1 S**, 45 (1967). [ Links ]

[9] R. M. More, J. Quant. Spectrosc. Radiat. Transfer **27**, 345 (1982). [ Links ]

[10] F. Perrot, Phys. Scripta **39**, 332 (1989). [ Links ]

[11] P. Martel, J.G. Rubiano, J.M. Gil, L. Doreste, and E. Mínguez, J. Quant. Spectrosc. Radiat. Transfer, **60**, 623 (1998). [ Links ]

[12] M. Kregar, Phys. Scripta **29**, 438 (1984), Phys. Scripta **31**, 246 (1985). [ Links ]

[13] H. O. Di Rocco, Braz. J. Phys. **22**, 227 (1992). [ Links ]

[14] E. U. Condon and H. Odabasi, *Atomic Structure* (University Press, Cambridge, 1980). [ Links ]

[15] K.-N. Huang et al., At. Data Nucl. Data Tables **18**, 249 (1976). [ Links ]

[16] H. O. Di Rocco, Revista Mexicana de Física **48**, 76 (2002). [ Links ]

Received on 22 November, 2003