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 [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text] and H3O(+). 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 H3O(+) is found to be almost 1/20 of the present theory at the cross section peak.

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: , , , , , , and H3O+. We outline the theory in Sec. 2, and our theoretical results are presented in Sec. 3. A recent experiment on the formation of H3O++ 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]: where t = T/B, u = U/B, , a0 is the Bohr radius (= 0.5292 Å), 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: 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 , , , , , , and H3O+ in Figs. 1–4. The molecular constants B, U, and N for the molecules are listed in Table 1. For all molecular ions except , molecular geometries were computed using a hybrid density functional (B3LYP) [10,11] with 6-31G(d) basis sets. For , B3LYP/6-31G(d) gave an incorrect molecular symmetry (C2 point group instead of C2), so the geometry was computed using frozen-core, second-order perturbation (MP2) theory with 6-31G(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-311G(d,p) basis sets. More accurate, correlated values of B were obtained for the outer-valence orbitals by using frozen-core Green’s function (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 981 [18].

Fig. 1

Electron-impact ionization cross sections for and .

Fig. 2

Electron-impact ionization cross sections for and .

Fig. 3

Electron-impact ionization cross section for and .

Fig. 4

Electron-impact ionization cross section for H3O+.

Table 1

Molecular point group, molecular orbitals (MO), electron binding energy B, kinetic energy U, and electron occupation number N for , , , , , , and H3O+

Molecule (point group)	MO	B(eV)	U(eV)	N
CH2+ (C2v)	1a1	319.07	436.58	2
	2a1	33.16	38.98	2
	1b2	27.04	29.90	2
	3a1	22.17	37.74	1
CH3+ (D3h)	1a1′	317.93	436.46	2
	2a1′	33.55	37.57	2
	1e′	25.59	29.56	4
CH4+ (C2v)	1a1	315.92	436.20	2
	2a1	33.67	33.40	2
	3a1	25.87	27.79	2
	1b2	24.48	28.73	2
	1b1	22.08	29.42	1
C_2H2+ (D∞h)	1σg	316.56	435.87	2
	1σu	316.48	436.63	2
	2σg	34.59	49.00	2
	2σu	27.96	35.66	2
	3σg	26.06	34.23	2
	1σu	21.00	31.42	3
C_2H4+ (D2)	1a	315.15	436.08	2
	1b1	315.12	436.44	2
	2a	33.52	40.73	2
	2b1	28.22	36.01	2
	1b2	24.41	27.06	2
	3a	23.30	35.07	2
	1b3	21.81	29.99	2
	2b3	19.21	29.54	1
C_2H6+ (D3d)	1a1g	314.28	436.27	2
	1a2u	314.28	436.30	2
	2a1g	31.76	33.19	2
	2a2u	29.16	38.10	2
	1eu	23.07	26.29	4
	1eg	21.38	30.07	4
	3a1g	19.04	30.07	1
H3O+ (C3v)	1a1	571.51	794.48	2
	2a1	45.59	72.42	2
	1e	29.90	52.11	4
	3a1	24.7a	65.41	2

Experimental value from Ref. [8].

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

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+. As another example, Bahati et al. [8] recently reported experimental cross sections for the process of Their measurement corresponds to process (5) only, and their peak cross section is 0.049 Å2 at T ≈ 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+ (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.