Scaling law of single ion–atom impact ionization cross sections of noble gases from He to Xe at strong perturbative energies*
Ren Ping-Yuana), Zou Xian-Rongb)†, Shao Jian-Xiongb), Wang Shi-Yaob), Zhou Manb), Zhou Wangb), Yang Ai-Xiangb), Yan Peng-Xuna), Chen Xi-Mengb)‡
School of Physical Science and Technology, Lanzhou University, Lanzhou 730000, China
School of Nuclear Science and Technology, Lanzhou University, Lanzhou 730000, China

Corresponding author. E-mail: zouxr@lzu.edu.cn

Corresponding author. E-mail: chenxm@lzu.edu.cn

*Project supported by the National Natural Science Foundation of China (Grants Nos. 11174116, and 11175075).

Abstract

We extend our previous work of a classical over-barrier ionization (COBI) model to calculate the single ionization cross sections of noble gases ranging from He to Xe at strong perturbative energies. The calculation results are in good agreement with extensive experimental data. The scaling law of single ion–atom impact ionization cross sections of noble gases on projectile charge q and energy E, also on target ionization energy I is drawn from the model.

Keyword: 34.50.Fa; 34.80.Dp; scaling law; single ionization; noble gases; strong perturbation
1. Introduction

Ion– atom ionizing collisions, which are of importance in atomic physics[1] and many applications in the fields of fusion, [2] accelerator, [3] and plasma physics, [4] have received much attention for several decades.

On the theoretical side, there are several approaches to calculating the cross sections.[57] These approaches include the classical trajectory Monte Carlo (CTMC), [8] the classical deposition models like the DEPOSIT code, [9] and quantum-mechanical calculations based on the first-order perturbation theory.[1015] The ionizations of noble atoms impacted by a projectile at high and low energies have already been understood well both in theory and in experiment.[12, 1620] In low energy range (E < 10  keV/amu, amu is the abbreviation for atomic unit per nucleon), the collision is adiabatic and the electron circulates many times around both nuclei. The electron energy in such a ‘ quasi-molecule’ is determined as a function of the inter-nuclear distance R. Adiabatic approximation such as the well-known over-barrier (OB) model proposed by Bohr et al., [21] can well describe the electron capture process. The OB treatment is applicable to electron capture by the projectile, [22, 23] but is not valid for the direct ionization process.

In a high energy range (E > MeV/amu), the electron interaction with the projectile takes place within a very short time, and the first-order perturbation method is valid. Born approximation, [12, 24] Bethe formula, [25] and the CDW-EIS theory[5] can be used for calculating the ionization cross sections and the comparison with available experimental data is shown in Refs.  [12] and [16]– [20]. All these theories indicate that for E > MeV/amu only collisions with sufficiently small impact parameters result in ionization, the energy is transferred to electrons via violent binary-encounters.

It is evident that the perturbation theories are not valid at strong perturbative energies (several tens to hundreds keV/amu), where the assumptions of free electron and binary encounter did not hold true. Though the calculation results given by the first principles quantum theories[2628] are in good agreement with experimental results, the knowledge of the underlying physics of ionization is still inadequate. On the application side, a universal formula or scaling law is needed to help us to understand the experimental data. Kaganovich et al. have reviewed the process of scaling law for He target in a wide energy range, [16] whereas the scaling law for many kinds of noble gases in strong perturbative energy range has not been reported yet.

The aim of this work is to obtain the projectile charge q and projectile velocity v, also the dependence of single ionization cross section of the target atom, through a classical calculation for the whole set of noble gas targets. The calculated single ionization cross sections from He to Xe impacted by various projectiles at energies ranging from 10  keV/amu to 800  keV/amu are in good agreement with the experimental data. Based on the calculation model, we propose a useful scaling law of single ionization cross sections for noble gases in a strong perturbative regime. The details of the model are presented in Section 2. The calculated results are compared with experimental results and the scaling law is discussed in Section 3. The atomic unit (a.u.) is used throughout the present paper.

2. Model

The model we used to calculate the single ionization cross sections of He, Ne, Ar, Kr, and Xe targets is based on our previous work of COBI[2931] model and also on Bohr’ s OB.[21] Now, we will review them briefly.

Two important ion– atom interaction distances are introduced in OB. The first distance is the release distance Rr, where the electron can be released from the target nucleus. The release distance Rr satisfies

where I is the ionization energy of the target electron, q is the projectile charge, and Z is the nuclear charge of the target atom.

The second distance is the capture distance RC, where the electron can be captured by the projectile, and it satisfies the following equation:

where v is the velocity of the projectile.

Equations (1) and (2) mean that only if the release occurs within the capture distance, can the released electron be captured, otherwise, the release is temporary. In OB, [21] the electrons which are released but not captured will go back to the target atom after collision. The energy level of the released electron has already been estimated by Bohr[21] and Niehaus, [22] and it is I + q/Rr, i.e., the ionization energy of the target electron plus the Stark-energy by the ion.

In COBI, we consider that the released but not captured electrons will be accelerated by the approaching ion. At some distance RI, the electrons will get enough kinetic energy from the ion to escape from the target atom, meanwhile, these electrons will not be captured, and they will be ionized. The ionization distance RI satisfies

which means that when the Stark-energy transferred to the kinetic energy of the electron is larger than the ionization energy of the quasi-molecular states, the electron will escape from the target atom. If it is not captured, it will be ionized. That is, in our model, the electrons released within the distance RC will be captured; those released electrons between Rr and RC will not be ionized until the ion enters the distance RI. For the given collision parameter b and velocity v, the probability that the electron will be released is

i.e., the ratio of collision duration of release process to the period of the bound electron T = 2π · n3/Z2, for bRr. The capture probability is

i.e., the ratio of collision duration of capture process to the period of the bound electron T, for bRC. Then, the ionization probability is PI (b) = Pr (b) − PC (b), for bRI.

Multiple electron probability will be extracted from these one-electron probabilities by using the independent event model (IEVM). For a multi-electron target, it has been found that IEVM is more adequate[32, 33] than the independent particle model (IPM). IPM ignores the variation of considerable ionization potential with the charge state, which leads to wrong probabilities for the high stages of ionization.

Fig.  1. Single ionization cross sections of He target impacted by various heavy ions. Lines represent calculation results and dots refer to experimental data.[27, 3441]

When we construct the cross section of a multi-electron system, it is necessary to distinguish these electrons with subscript i. They satisfy the following equations:

where qi and Zi are the effective charges of the incident ion and the target, felt by the i-th released electron of the target. for the released electron loses the ability to shield target and projectile nucleus. The values of PCi (b) and PIi (b) are calculated from Rri, RCi, and RIi. With the statistical method of multinomial distribution, the related single ionization probability can be obtained as

The cross section is obtained by integrating the probability over all impact parameters and expressed as

The calculation results are plotted in Fig.  1 together with the experimental data within the same energy range.

3. Results and discussion

We calculate the single ion– atom impact ionization cross sections of He impacted by various kinds of heavy ions in the energy range from tens to hundreds keV/amu. The calculation results are shown in Figs.  1 and 2 and compared with extensive experimental data. It is clear that the model results are in good agreement with the experimental data.

Fig.  2. Single ionization cross section of He target impacted by U90+ , the solid line represents the calculation results and square dots refer to experimental data.[42]

The model well describes the q- and E-dependent single ionization cross sections. The calculation results and the data both show that (i) the single ionization cross sections of helium increase very rapidly at low energies (below about 100  keV/amu) and reach their maxima at energies ranging from 100  keV/amu to 400  keV/amu, and then drop off very slowly at higher energies; (ii) the maximal ionization cross sections increase significantly with the increase of projectile charge q, while the maximum positions increase with q slightly.

The behavior of single ionization cross section on q and E can be understood qualitatively by the model, for the dependence of the single ionization cross section on q and E is similar to that of the single ionization probability PI, and we will discuss the ionization probability PI(0) at b = 0.

Let h = (2Rr)/v · 1/T and g = (RrRC)/Rr, where h is the ratio of collision duration 2Rr/v to the period T of the target electron. Obviously, h is the release probability during the collision, and it is equal to the unity at low energies and decreases according to 1/v for high energies. g is the uncaptured fraction of the released electrons, g = 1 − (2q/Rr) · (1/v2), which increases strongly according to v2 and is close to the unity for high v.

The uncaptured fraction of the released electrons g increases according to v2; this is the main reason for the increase of ionization cross sections at low energies. The interaction time h decreases according to causing the ionization cross sections to decrease very slowly for high energies. At intermediate energies, the uncaptured fraction g increases while the interaction time h decreases. The competition between them leads to the occurrance of maximum ionization cross sections.

Moreover, capture and ionization are both dependent on the projectile charge. However, increasing collision velocity reduces the capture cross section more rapidly than the ionization cross section. In order to reduce the capture cross section, the effect of collision velocity must overcome the capture process, which means that higher and higher velocities are needed for a significant ionization contribution when the projectile charge state is increased. This is the reason why the ionization cross section maximum shifts toward higher velocities for higher projectile charges. Approximately, the maximum position, where the maximum ionization cross section is reached, can be obtained by setting the derivative of Eq.  (7) to be zero. The estimated maximum position Emax is found to be proportional to 6q/Rr. In our calculation, we find that the values of Rr can be fitted well by the curve of q0.4 (see Fig.  3). Thus, the maximum position will show the q0.6 dependence. This is why the maximum position of Emax does not vary significantly with q.

Fig.  3. Calculation results of Rr and RI and their fitting curves.

On the other hand, the maximum ionization cross section is shown to strongly increase with q. According to Eq.  (3), the maximum ionization cross section can be estimated roughly at The calculation result of RI shows the q0.65 dependence (see Fig.  3). So, we can expect that the maximum ionization cross section will be dependent on q in the form of about q1.3.

In order to obtain the target dependence of the single ionization cross section, we calculate the single ionization cross section of noble gases ranging from He to Xe impacted by He+ ion. The calculated results are shown in Fig.  4 and compared with the available experimental data. We can see that the model gives good estimation to the experimental data. Based on the present model, we can estimate the target dependence of the maximum ionization cross section qualitatively. According to Eq.  (3), the dependence of RI on target ionization energy I is 1/I, and then the maximum ionization cross sections of different targets will show the (1/I)2 dependence.

Fig.  4. Single ionization cross sections of noble gases from He to Xe impacted by He+ ion. Lines represent the calculation results and dots refer to the experimental data.[43]

Now, we suggest that the q, v, and I dependences of the maximum ionization cross section are σ maxq1.3/I2 and Emaxq0.6. Based on such dependences, we can obtain a scaling law of the single ionization cross sections of noble gases. As a result, the single ionization cross section of a target X having the first ionization energy I and being impacted by a projectile with charge q and energy E can be estimated by

and E* ≈ E · q0.6. For clarifying the validity of such a scaling law drawn from the model, we scale the single ionization cross sections of the whole set of noble gases impacted by various kinds of heavy ions in Fig.  5. In this figure, we can find that the extensive experimental data can be fitted well by the scaling law. This scaling law builds a convenient connection between the usual proton data and other data for different collision systems at energies ranging from tens to hundreds keV/amu.

Fig.  5. Scaled single ionization cross sections of noble gases ranging from He to Xe impacted by various kinds of heavy ions.[27, 3443] The dashed line is used only to guide the eyes.

4. Conclusions

The single ionization cross sections of noble gases ranging from He to Xe, are calculated at strong perturbative energies. The calculation results are in good agreement with the extensive experimental data. Based on the model, the scaling law is drawn to describe the q, E, and I dependences of single ionization cross sections at such energies. It is found that the extensive experimental data can be fitted well by the scaling law of

and E* ≈ E · q0.6. This work provides a convenient way to estimate the single ionization cross sections of noble gas targets impacted by various kinds of heavy ions from the easily obtained proton data.

Reference
1 Datz S, Drake G W F, Galagher T F, Kleinpoppen H and Zu Putlitz G 1999 Rev. Mod. Phys. 71 S223 DOI:10.1103/RevModPhys.71.S223 [Cited within:1]
2 Logan B G, Celata C M, Kwan J W, Lee E P, Leitner M, Seidl P A, Yu S S, Barnard J J, Friedman A, Meier W R and Davidson R C 2002 Laser Part. Beams 20 369 DOI:10.1017/S0263034602203018 [Cited within:1] [JCR: 2.016]
3 Beyer H F and Shevelko V P 1999 Atomic Physics with Heavy Ions Berlin Springer [Cited within:1]
4 Govil R, Leemans W P, Backhaus E Yu and Wurtele J S 1999 Phys. Rev. Lett. 83 3202 DOI:10.1103/PhysRevLett.83.3202 [Cited within:1] [JCR: 7.943]
5 Janev R K, Presnyakov L P and Shevelko V P 1999 Physics of Highly Charged Ions Berlin Springer [Cited within:2]
6 McDowell M R C and Coleman J P 1970 Introduction to the Theory of Ion–Atom Collisions Amsterdam North-Holland [Cited within:1]
7 Bransden B H and McDowell M R C 1992 Charge Exchange and the Theory of Ion–Atom Collisions Oxford Clarendon [Cited within:1]
8 Olson R E, Ullrich J and Schmidt-Böcking H 1989 Phys. Rev. A 39 5572 DOI:10.1103/PhysRevA.39.5572 [Cited within:1] [JCR: 3.042]
9 Shevelko V P, Litsarev M S and Tawara H 2008 J. Phys. B: At. Mol. Opt. Phys. 41 115204 DOI:10.1088/0953-4075/41/11/115204 [Cited within:1]
10 Land au L D and Lifshitz E M 1985 Quantum Mechanics Reading, MA Addison-Wesley [Cited within:1]
11 Bohr N 2010 K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 18 8 [Cited within:1]
12 Haugen H K, Andersen L H, Hvelplund P and Knudsen H 1982 Phys. Rev. A 26 1950 [Cited within:3] [JCR: 3.042]
13 Theodosiou C E 1988 Phys. Rev. A 38 4923 DOI:10.1103/PhysRevA.38.4923 [Cited within:1] [JCR: 3.042]
14 Sand ner W and Theodosiou C E 1990 Phys. Rev. A 42 5208 DOI:10.1103/PhysRevA.42.5208 [Cited within:1] [JCR: 3.042]
15 McGuire J H 1982 Phys. Rev. A 26 143 DOI:10.1103/PhysRevA.26.143 [Cited within:1] [JCR: 3.042]
16 Kaganovich I D, Startsev E and Davidson R C 2006 New J. Phys. 8 278 DOI:10.1088/1367-2630/8/11/278 [Cited within:3] [JCR: 4.063]
17 Schlachter A S, Berkner K H, Graham W G, Pyle R V, Schneider P J, Stalder K R, Stearns J W, Tanis J A and Olson R E 1981 Phys. Rev. A 23 2331 DOI:10.1103/PhysRevA.23.2331 [Cited within:1] [JCR: 3.042]
18 Berg H, Ullrich J, Bernstein E, Unverzagt M, Spielbergert L, Euler J, Schardt D, Jagutzki O, Schmidt-Bocking H, Mann R, Mokler P H, Hagmann S and Fainstein P D 1992 J. Phys. B: At. Mol. Opt. Phys. 25 3655 DOI:10.1088/0953-4075/25/17/017 [Cited within:1] [JCR: 2.031]
19 Gilespie G H 1981 Phys. Rev. A 24 608 DOI:10.1103/PhysRevA.24.608 [Cited within:1] [JCR: 3.042]
20 Knudsen H and Reading J F 1992 Phys. Rep. 212 107 DOI:10.1016/0370-1573(92)90013-P [Cited within:2] [JCR: 22.929]
21 Bohr N and Lindhard J 1954 K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 28 7 [Cited within:4]
22 Knudsen H, Haugen H K and Hvelplund P 1981 Phys. Rev. A 23 597 DOI:10.1103/PhysRevA.23.597 [Cited within:2] [JCR: 3.042]
23 Niehaus A 1986 J. Phys. B: At. Mol. Opt. Phys. 19 2925 DOI:10.1088/0022-3700/19/18/021 [Cited within:1] [JCR: 2.031]
24 Garcia J D 1970 Phys. Rev. A 1 280 DOI:10.1103/PhysRevA.1.280 [Cited within:1] [JCR: 3.042]
25 Bethe H 1930 Ann. Phys. 397 325 DOI:10.1002/(ISSN)1521-3889 [Cited within:1] [JCR: 0.25]
26 Fainstein P D, PonceV H and Rivarola R D 1991 J. Phys. B: At. Mol. Opt. Phys. 24 3091 DOI:10.1088/0953-4075/24/14/005 [Cited within:1] [JCR: 2.031]
27 Wu W, Deveney E F, Datz S, Desai D D, Krause H F, Sand ers J M, Vane C R, Cocke C L and Giese J P 1996 Phys. Rev. A 53 2367 DOI:10.1103/PhysRevA.53.2367 [Cited within:1] [JCR: 3.042]
28 Wu W, Cocke C L, Giese J P, Melchert F, Raphaelian M L A and Stöckli M 1995 Phys. Rev. Lett. 75 1054 DOI:10.1103/PhysRevLett.75.1054 [Cited within:1] [JCR: 7.943]
29 Shao J X, Chen X M, Liu Z Y, Qi R and Zou X R 2008 Phys. Rev. A 77 042711 DOI:10.1103/PhysRevA.77.042711 [Cited within:1] [JCR: 3.042]
30 Shao J X, Chen X M, Zou X R, Wang X A and Lou F J 2008 Phys. Rev. A 78 042701 DOI:10.1103/PhysRevA.78.042701 [Cited within:1] [JCR: 3.042]
31 Shao J X, Zou X R, Chen X M, Zhou C L and Qiu X Y 2011 Phys. Rev. A 83 022710 DOI:10.1103/PhysRevA.83.022710 [Cited within:1] [JCR: 3.042]
32 Mendez L, Errea L F, Illescas C, Rabadan I, Pons B and Riera A 2008Radiation Damage in Biomolecular Systems: Proceedings of 5th International Conference (RADAM 2008), AIP Conf. Proc. No.  1080 51 [Cited within:1]
33 Kabachnik N M, Kondratyev V N, Roller-Lutz Z and Lutz H O 1997 Phys. Rev. A 56 2848 DOI:10.1103/PhysRevA.56.2848 [Cited within:1] [JCR: 3.042]
34 Shah M B and Gilbody H B 1985 J. Phys. B: At. Mol. Opt. Phys. 18 899 DOI:10.1088/0022-3700/18/5/010 [Cited within:1] [JCR: 2.031]
35 Shah M B, McCallion P and Gilbody H B 1989 J. Phys. B: At. Mol. Opt. Phys. 22 3037 DOI:10.1088/0953-4075/22/19/018 [Cited within:1] [JCR: 2.031]
36 DuBois R D 1987 Phys. Rev. A 36 2585 DOI:10.1103/PhysRevA.36.2585 [Cited within:1] [JCR: 3.042]
37 Woitke O, Zavodszky P A, Ferguson S M, Houck J H and Tanis J A 1998 Phys. Rev. A 57 2692 DOI:10.1103/PhysRevA.57.2692 [Cited within:1] [JCR: 3.042]
38 Knudsent H, Andersent L H, Hvelplundt P, Astner G, Cederquistl H, Danared H, Liljeby L and Rensfelt K G 1984 J. Phys. B: At. Mol. Opt. Phys. 17 3545 [Cited within:1] [JCR: 2.031]
39 Datz S, Hippler R, Anderson L H, Dittner P F, Knudsen H, Krause H F and Miller P D 1990 Phys. Rev. A 41 3559 DOI:10.1103/PhysRevA.41.3559 [Cited within:1] [JCR: 3.042]
40 Damsgaard H, Haugen H K, Hvelplund P and Knudsen H 1983 Phys. Rev. A 27 112 DOI:10.1103/PhysRevA.27.112 [Cited within:1] [JCR: 3.042]
41 Shinpaugh J L, Sand ers J M, Hall J M, Lee D H, Schmidt-Bocking H, Tipping T N, Zouros T J M and Richard P 1992 Phys. Rev. A 45 2922 DOI:10.1103/PhysRevA.45.2922 [Cited within:1] [JCR: 3.042]
42 Berg H, Jagutzki O, Dorner R, DuBois R D, Kelbch C and Schmidt-Bocking H 1992 Phys. Rev. A 46 5539 DOI:10.1103/PhysRevA.46.5539 [Cited within:1] [JCR: 3.042]
43 Santos A C F, Melo W E, Sant’Anna M M, Siguad G M and Montenegro E C 2001 Phys. Rev. A 63 062717 DOI:10.1103/PhysRevA.63.062717 [Cited within:1] [JCR: 3.042]