Electron-impact total ionization cross sections of hydrocarbon ions.

The Binary-Encounter-Bethe (BEB) model for electron-impact total ionization cross sections has been applied to [CH.sup.+.sub.2], [CH.sup.+.sub.3], [CH.sup.+.sub.4], [C.sub.2][H.sup.+.sub.2], [C.sub.2][H.sup.+.sub.4], [C.sub.2][H.sup.+.sub.6], and [H.sub.3][O.sup.+]. The cross sections for the hydrocarbon ions are needed for modeling cool plasmas in fusion devices. No experimental data are available for direct comparison. Molecular constants to generate total ionization cross sections at arbitrary incident electron energies using the BEB formula are presented. A recent experimental result on the ionization of [H.sub.3][O.sup.+] is found to be almost 1/20 of the present theory at the cross section peak.

Key words: [CH.sup.+.sub.2]; [CH.sup.+.sub.3], [CH.sup.+.sub.4]; [C.sub.2][H.sup.+.sub.2]; [C.sub.2][H.sup.+.sub.4]; [C.sub.2][H.sup.+.sub.6]; and [H.sub.3][O.sup.+]; electron-impact ionization; molecular ions.

1. Introduction

Ionization cross sections for atomic and molecular ions are among the critical data needed in modeling plasmas in fusion devices. Hydrocarbon molecules and their ion fragments are formed inside a tokamak in edge plasmas and near a divertor. The Binary-Encounter-Bethe (BEB) model (1) has successfully generated reliable total ionization cross sections of small as well as large molecules (2-6). The BEB model combines a modified form of the Mott cross section with the asymptotic form of the Bethe theory (i.e., high incident energy T) for electron-impact ionization of a neutral atom or molecule. The original BEB model was slightly modified for applications to atomic and molecular ions (7).

In this article we apply the modified BEB formula for ions to hydrocarbon ions of interest to magnetic fusion: [CH.sup.+.sub.2], [CH.sup.+.sub.3], [CH.sup.+.sub.4], [C.sub.2][H.sup.+.sub.2], [C.sub.2][H.sup.+.sub.4], [C.sub.2][H.sup.+.sub.6], and [H.sub.3][O.sup.+]. We Outline the theory in Sec. 2, and our theoretical results are presented in Sec. 3. A recent experiment on the formation of [H.sub.3][O.sup.++] by electron impact (8) is compared to the present theory in Sec. 3.

2. Outline of Theory

The BEB formula for ionizing an electron from a molecular orbital of a neutral molecule by electron impact is (1):

[[sigma].sub.BED] = S/t + u + 1 [ln t/2 (1 - 1/[t.sup.2])+ 1 - 1/t - ln t/t + 1], (1)

where t = T/B, u = U/B, S = 4[pi][a.sup.2.sub.0] N [R.sup.2]/[B.sup.2], [a.sub.0] is the Bohr radius (= 0.5292 A), R is the Rydberg energy (= 13.6057 eV), T is the incident electron energy, and N, B, and U are the electron occupation number, the binding energy, and the average kinetic energy of the orbital, respectively.

In Eq. (1), the terms in the square brackets are based on the Mott theory and the Bethe theory. However, the denominator t + u + 1 is based on a plausible, but less rigorous argument, i.e., the effective kinetic energy of the incident electron seen by the bound target electron should be the incident electron energy T plus the potential energy U + B of the target electron (9). Hence the T in the denominator of the original Mott and Bethe theories was replaced by T + U + B, or t+ u + 1 in Eq. (1), where B is used as the energy unit.

The net effect of using t + u + 1 instead of t in the denominator of Eq. (1) is to reduce substantially the cross section near the ionization threshold. This modification was found not only to be effective but also absolutely necessary to have the theory agree with reliable experimental ionization cross sections near the threshold for many neutral atoms and molecules.

In a previous article (7) for singly charged molecular ions, we have shown that the denominator t + u + 1 is replaced by t + (u + 1)/2 to generate ionization cross sections in good agreement with available experimental data. The modified BEB equation for singly charged ions is:

[[sigma].sub.ion] = S/t + (u+1)/2[ln t/2(1 - 1/[t.sup.2]) + 1 - 1/t - ln t/t + 1]. (2)

Equation (2) is as simple as the BEB formula for neutral targets, Eq. (1), and does not require any more input data than the original BEB formula.

3. Theoretical Results

We present the BEB cross sections from Eq. (2) for [CH.sup.+.sub.2], [CH.sup.+.sub.3], [CH.sup.+.sub.4], [C.sub.2][H.sup.+.sub.4], [C.sub.2][C.sup.+.sub.6], and [H.sub.3][O.sup.+] in Figs. 1-4. The molecular constants B, U, and N for the molecules are listed in Table 1. For all molecular ions except [CH.sup.+.sub.4], molecular geometries were computed using a hybrid density functional (B3LYP) (10,11) with 6-31G(d) basis sets. For [CH.sup.+.sub.4] B3LYP/6-31G(d) gave an incorrect molecular symmetry ([C.sub.2] point group instead of [C.sub.2v]), so the geometry was computed using frozencore, second-order perturbation (MP2) theory with 631G(d) basis sets. The B3LYP or MP2 geometries were used for all subsequent calculations of B and U. Kinetic energies U for all orbitals, and binding energies B for the inner orbitals, were calculated at the Hartree-Fock (HF) level using 6-311+G(d,p) basis sets. More accurate, correlated values of B were obtained for the outer-valence orbitals by using frozen-core Green's functio n (OVGF) methods (12,13) and 6-311+G(d,p) basis sets. For the important threshold ionization, B values were obtained by using frozen-core coupled cluster theory [CCSD(T)], with the single and double excitation operators included iteratively (14) and the contribution from connected triples estimated perturbatively (15). Dunning's correlation-consistent valence-triple-zeta (cc-pVTZ) basis sets (16) were used for the CCSD(T) calculations. The HF calculations were performed using the GAMESS (17) program package; all other calculations employed Gaussian 98 (1) (18).

In general, when an electron collides with a molecular ion we get

[e.sup.-] + [AB.sup.+] [right arrow] A + [B.sup.+] + [e.sup.-], (3)

or [A.sup.+] + B + [e.sup.-]. (4)

[e.sup.-] + [AB.sup.+] [right arrow] [AB.sup.++] + 2[e.sup.-], (5)

or [A.sup.+] + [B.sup.+] + 2[e.sup.-], (6)

or [A.sup.++] + B + 2[e.sup.-], (7)

or A + [B.sup.++] + 2[e.sup.-]. (8)

Processes (3) and (4) are dissociation without ionization, while processes (5) through (8) are the ionizing events described by the BEB model. The model calculates the sum of all processes (5) through (8) that lead to the ejection of a bound electron. Moreover, the model also assumes-erroneously--that all energy transfers from the incident electron to the target molecule that exceed the ionization energy of a given molecular orbital result in ionization. This is an assumption common to all binary-encounter type theories. Although such an assumption may be valid for atoms, molecules may dissociate without ionizing even if energy transfers exceed the orbital binding energies. If processes (3) and (4) are significant for energy transfers above orbital binding energies, then the BEB model will overestimate ionization cross sections. More discussions on this point can be found in Ref. (5).

Experimentally, the production of doubly charged ions can be detected directly when the doubly charged ions have reasonably long lifetimes. In reality, most doubly charged molecular ions quickly dissociate into two singly charged fragments, making it almost impossible to distinguish processes (3) and (4) from the dissociation of doubly charged ions by Coulomb repulsion, i.e., process (6). For this reason, it is difficult to distinguish processes (3) and (4) from process (6) simply by detecting singly charged ions unless coincidence measurements of all products are performed. The usual experimental procedure is to measure the cross section for producing any ion, i.e., (3) through (8). Then, processes (3) and (4) are measured separately, and subtracted from the total ion production cross section. This subtraction introduces large uncertainties in the resulting experimental ionization cross sections.

The direct measurement of only the doubly charged ions-processes (5), (7), and (8)-tends to produce small cross sections compared to the total ionization cross section because of the high probability for the rapid break-up of the doubly charged ions as shown in Ref. (7) for [CO.sup.+]. As another example, Bahati et al. (8) recently reported experimental cross sections for the process of

[e.sup.-] + [H.sub.3][O.sup.+] [right arrow] [H.sub.3][O.sup.++] + 2[e.sup.-] (9)

Their measurement corresponds to process (5) only, and their peak cross section is 0.049 [A.sup.2] at [approximately equal to] 125 eV. The position of the peak is in agreement with the BEB cross section in Fig. 4, but the magnitude is almost 1/20 of the BEB cross section. This discrepancy and the similar discrepancy in [CO.sup.+] (experiment is lower by a factor of 1/12 at the peak) measured by the same group (19) is a strong indication that either dissociative ionization is the dominant process, or the break-up of the doubly charged ions is faster than the experimental capability to detect.

[FIGURE 1 OMITTED]

[FIGURE 2 OMITTED]

[FIGURE 3 OMITTED]

[FIGURE 4 OMITTED]

Table 1

Molecular point group, molecular orbitals (MO), electron binding energey
B, kinetic energy U, and electron occupation number N for
[CH.sup.+.sub.2], [CH.sup.+.sub.3], [CH.sup.+.sub.4],
[C.sub.2][H.sup.+.sub.2], [C.sub.2][H.sup.+.sub.4],
[C.sub.2][H.sup.+.sub.6], and [H.sub.3][O.sup.+]

Molecule (point group)           MO                 B(eV)      U(eV)  N

[CH.sup.+.sub.2] ([C.sub.2v])    1[a.sub.1]        319.07     436.58  2
                                 2[a.sub.1]         33.16      38.98  2
                                 1[b.sub.2]         27.04      29.90  2
                                 3[a.sub.1]         22.17      37.74  1

[CH.sup.+.sub.3] ([D.sub.3h])    1[a.sub.i]        317.93     436.46  2
                                 2[a.sub.i]         33.55      37.57  2
                                 1e'                25.59      29.56  4

[CH.sup.+.sub.4] ([C.sub.2v])    1[a.sub.1]        315.92     436.20  2
                                 2[a.sub.1]         33.67      33.40  2
                                 3[a.sub.1]         25.87      27.79  2
                                 1[b.sub.2]         24.48      28.73  2
                                 1[b.sub.1]         22.08      29.42  1

[C.sub.2][H.sup.+.sub.2]         1[[sigma].sub.g]  316.56     435.87  2
([D.sub.[infinity]h])            1[[sigma].sub.u]  316.48     436.63  2
                                 2[[sigma].sub.g]   34.59      49.00  2
                                 2[[sigma].sub.u]   27.96      35.66  2
                                 3[[sigma].sub.g]   26.06      34.23  2
                                 1[[pi].sub.u]      21.00      31.42  3

[C.sub.2][H.sup.+.sub.4]         1b                315.15     436.08  2
([D.sub.2])                      1[b.sub.1]        315.12     436.44  2
                                 2a                 33.52      40.73  2
                                 2[b.sub.1]         28.22      36.01  2
                                 1[b.sub.2]         24.41      27.06  2
                                 3a                 23.30      35.07  2
                                 1[b.sub.3]         21.81      29.99  2
                                 2[b.sub.3]         19.21      29.54  1

[C.sub.2][H.sup.+.sub.6]         1[a.sub.1g]       314.28     436.27  2
([D.sub.3d])                     1[a.sub.2u]       314.28     436.30  2
                                 2[a.sub.1g]        31.76      33.19  2
                                 2[a.sub.2u]        29.16      38.10  2
                                 1[e.sub.u]         23.07      26.29  4
                                 1[e.sub.g]         21.38      30.07  4
                                 3[a.sub.1g]        19.04      30.07  1

[H.sub.3][O.sup.+] ([C.sub.3v])  1[a.sub.1]        571.51     794.48  2
                                 2[a.sub.1]         45.59      72.42  2
                                 1e                 29.90      52.11  4
                                 3[a.sub.1]         24.7 (a)   65.41  2

(a)Experimental value from Ref. (8).

Acknowledgments

We gratefully acknowledge partial financial support by the Office of Fusion Energy Sciences of the U.S. Department of Energy and by the Advanced Technology Program of NIST.

Accepted: December 18, 2001.

(1.) Certain commercial equipment, instruments, or materials are identified in this paper to foster understanding. Such identification does not imply recommendation or endorsement by the National Institute of Standards and Technology, nor does it imply that the materials or equipment identified are necessarily the best available for the purpose.

4. References

(1.) Y.-K. Kim and M. E. Rudd, Phys. Rev. A 50, 3954 (1994).

(2.) W. Hwang, Y.-K. Kim, and M. E. Rudd, J. Chem. Phys. 104, 2956 (1996).

(3.) Y.-K. Kim, W. Hwang, N. M. Weinberger, M. A. Ali, and M. E. Rudd, J. Chem, Phys. 106, 1026 (1997).

(4.) M. A. Ali, Y.-K. Kim, W. Hwang, N. M. Weinberger, and M. E. Rudd, J. Chem. Phys. 106, 9602 (1997).

(5.) H. Nishimura, W. M. Huo, M. A. Ali, and Y.-K. Kim, J. Chem. Phys. 110, 3811 (1999).

(6.) Y.-K. Kim and M. E. Rudd, Comments At. Mol. Phys. 34, 309 (1999).

(7.) Y.-K. Kim, K. K. Irikura, and M. A. Ali, J. Res. Natl. Inst. Stand. Technol. 105, 285 (2000).

(8.) E. M. Bahati, J. J. Jureta, H. Cherkani-Hassani, and P. Defrance, J. Phys. B 34, L333 (2001).

(9.) A. Burgess, Proc. 3rd Int. Conf. on Electronic and Atomic Collisions, London, 1963, M. R. C. McDowell, ed., North Holland, Amsterdam (1964) p. 237; Proc. Symp. on Atomic Collision Processes in Plasmas, Culham, AERE Rept. 4818 (1964) p. 63.

(10.) A. D. Becke, J. Chem. Phys. 98, 5648 (1993).

(11.) P. J. Stephens, F. J. Devlin, C. F. Chabalowski, and M. J. Frisch, J. Phys. Chem. 98, 11623 (1994).

(12.) W. von Niessen, J. Schirmer, and L. S. Cederbaum, Comput. Phys. Rep. 1, 57 (1984).

(13.) V. G. Zakrzewski and J. V. Ortiz, J. Phys. Chem. 100, 13979 (1996).

(14.) G. D. Purvis and R. J. Bartlett, J. Chem. Phys. 76, 1910 (1982).

(15.) K. Raghavachari, G. W. Trucks, J. A. Pople, and M. Head-Gordon, Chem. Phys. Lett. 157, 479 (1989).

(16.) T. H. Dunning, Jr., J. Chem. Phys. 90, 1007 (1989).

(17.) M. W. Schmidt, K. K. Baldridge, J. A. Boatz, S. T. Elbert, M. S. Gordon, J. H. Jensen, S. Koseki, N. Matsunaga, K. A. Nguyen, S. I. Su, T. L. Windus, M. Dupuis, and J. A. Montgomery, J. Comput. Chem. 14, 1347 (1993).

(18.) M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, V. G. Zakrzewski, J. A. Montgomery, Jr., R. E. Stratmann, J. C. Burant, S. Dapprich, J. M. Millam, A. D. Daniels, K. N. Kudin, M. C. Strain, O. Farkas, J. Tomasi, V. Barone, M. Cossi, R. Cammi, B. Mennucci, C. Pomelli, C. Adamo, S. Clifford, J. Ochterski, G. A. Petersson, P. Y. Ayala, Q. Cui, K. Morokuma, D. K. Malick, A. D. Rabuck, K. Raghavachari, J. B. Foresman, J. Cioslowski, J. V. Ortiz, B. B. Stefanov, G. Liu, A. Liashenko, P. Piskorz, I. Komaromi, R. Gomperts, R. L. Martin, D. J. Fox, T. Keith, M. A. Al-Laham, C. Y. Peng, A. Nanayakkara, C. Gonzalez, M. Challacombe, P. M. W. Gill, B. Johnson, W. Chen, M. W. Wong, J. L. Andres, C. Gonzalez, M. Head-Gordon, E. S. Replogle, and J. A. Pople, Gaussian, Inc., Pittsburgh, PA (1998).

(19.) D. S. Belic, D. J. Yu, A. Siari, and P. Defrance, J. Phys. B 30, 5535 (1997).

About the authors: Karl K. Irikura is a chemist in the Physical and Chemical Properties Division in the NIST Chemical Science and Technology Laboratory. Yong-Ki Kim is a physicist in the Atomic Physics Division in the NIST Physics Laboratory. M. Asgar Ali is a professor in the Department of Chemistry at Howard University, and a guest researcher in the Atomic Physics Division at NIST. The National Institute of Standards and Technology is an agency of the Technology Administration, U.S. Department of Commerce.