Electron-impact total ionization cross sections of hydrocarbon ions - Statistical Data Included

Journal of Research of the National Institute of Standards and Technology, Jan, 2002 by Karl K. Irikura, Yong-Ki Kim, M.A. Ali

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).