1. Introduction
The ITER project aims to demonstrate the feasibility of fusion as a large-scale energy source. It is currently under construction in Cadarache, France, and relies on the efforts of 35 countries to build the world’s largest tokamak. Since the reactor has carbon components in its diverters, accurate data on the collision of electrons with carbon atoms are needed for plasma diagnostics, such as impurity influx studies [
1].
Despite the potential interest in this topic, experimental reports on the ionization of C by electron impact are quite scarce. Total ionization cross-sections measured by Wang and Crawford are contained within an MIT internal report, which dates back to 1971, and were not published [
2]. An accuracy of ±30% was estimated by the authors but practical details were not made explicit. A few years later, in 1978, Brook et al. published total cross-sections using the fast atom-beam and crossed electron-beam technique [
3]. Their results, though consistent in shape, are lower than those of Wang and Crawford from the threshold up to an impact energy of 1 keV. No new data have been reported since then.
From a theoretical point of view, pioneering calculations by Peach, Omidvar et al. and McGuire made use of the Born–Ochkur approximation [
4,
5], the Born approximation [
6], and a generalized oscillator strength formulation of the Born approximation [
7], respectively. While the results of Peach (Refs. [
4,
5]) were consistent with the data of Wang and Crawford, the calculation of McGuire (Ref. [
7]) supported the data of Brook et al. The results of Omidvar et al. (Ref. [
6]) were right in between both datasets. During the last decade, efforts have been devoted to calculating the total ionization cross-section, either from the ground state or from excited states, by means of highly numerically intensive methods, such as the time-dependent close-coupling (TDCC), the time-independent distorted wave method (TIDW), the R-matrix-with-pseudostates (RMPS) and the B-spline R-matrix-with-pseudostates (BSR) [
8,
9]. These methods have provided results that are in fair agreement with the data of Brook et al. but were only extended up to 60 eV (TDCC, TIDW, RMPS) and 100 eV (BSR).
In the present work, we provide a complementary view of this collision system and calculate the total ionization cross-section for collisions of electrons with C by means of the three-body classical trajectory Monte Carlo (CTMC) method [
10,
11]. The CTMC is a non-perturbative method that, in contrast to highly numerically intensive methods, does not rely on huge basis sets, allowing its implementation in an ample impact energy range. Due to the large difference in the ionization potentials, the K-shell ionization cross-sections are negligible compared to those calculated for the outer shell. Therefore, the present calculations are restricted to the L-shell. Results obtained by means of the simple addition rule are contrasted to those provided by the independent electron model (IEL) and the independent event model (IEV). Limitations inherent to the classical treatment are identified and analyzed.
Atomic units are used throughout this work unless otherwise stated.
2. Theoretical Method
In this work, the three-body CTMC model in its microcanonical formulation is employed. Hamilton’s equations are solved for the mutually interacting three-body system, in which the two electrons interact with each other via a Coulomb potential and with the ion core via the potential model derived from Hartree–Fock calculations by Green et al. [
12] that was later on generalized by Garvey et al. [
13],
This functional form allows the incorporation of the electronic screening effect, providing a more solid description of the collision dynamics compared to the use of fixed effective charges together with Coulomb potentials. Specific values for the parameters
N,
ξ, and
η employed throughout the present study to represent the interaction of an electron with the C
+ and C
2+ ions are shown in
Table 1. The variation in the dynamical charge
with the electron-ion core distance is illustrated in
Figure 1. As a general trend, the electron evolving in such model potential sees an asymptotic charge of +1 which gradually increases towards +
Z as the electron-core distance tends to zero.
Collisional events are tracked until the momenta of all the particles are considered to be converged. The ionization probability as a function of the impact parameter
b is determined by the ratio between the number of ionization events
Nion and the total number of evaluated trajectories
Ntot.
We use three different approaches to evaluate the total ionization probability as a function of the impact parameter:
where we add the ionization probabilities from the 2s and 2p orbitals.
- (ii)
Independent electron model (IEL)
where the multielectronic character of the target is considered by weighting the probability of electron emission from a given orbital by the probability of not removing additional electrons from the same or other orbitals of the neutral carbon atom. This model considers that the electron orbitals do not significantly change during the electron removal process.
- (iii)
Independent event model (IEV)
In this case, the probabilities P′ in the weight factor represent probabilities of electron emission from the sequential C+ ion.
Contributions from C(1s) are not explicitly considered in our work. We point out that explicit calculations for the ionization TCS from C(1s) at the punctual impact energies 500 eV, 750 eV, and 1000 eV were found negligible compared to the contributions of the 2s and 2p levels due to the large ionization potential of the C(1s) orbital.
The corresponding total ionization cross-section (TCS) is finally given by
with
k = (i),
(ii),
(iii). Here,
bmax is the maximum impact parameter that can lead to ionization.
3. Results
In
Figure 2 we show the TCS up to an impact electron energy of 1 keV and provide a quick view of our present state of knowledge. As indicated in the Introduction, it can be seen that early theoretical studies were far from conclusive regarding the discrepancy exhibited by the available experimental data. During the last decade, the implementation of numerically intensive techniques, such as TDCC, RMPS, and BSR lent credence to the dataset of Brook and the early theoretical prediction by McGuire. Regarding the present CTMC, results obtained by using the simple addition rule are in agreement with the experimental data and the theoretical results of McGuire at impact energies greater than about 100 eV. In contrast, the maxima of the TCS exhibit a clear shift towards lower impact energies compared to the data, and a clear overestimation of the threshold region is obtained. Similar trends have been reported in CTMC studies of electron impact ionization of hydrogen and helium [
14].
With the objective of determining whether this shift is provided by the consideration of separate contributions from the 2s and 2p orbitals, in
Figure 3 we analyze the IEM and IEV description. We clearly observe that they do not lead to any improvement in this case, neither in terms of peak position or magnitude. These statements are confirmed in
Figure 3b by a detailed inspection of
bPk(b) as a function of b. We observe that the inclusion of the weight factors in Equations (6) and (7) leads to an energy-dependent global scaling factor but does not modify the shape of the distributions.
In order to gain insight into the limitations of the present classical procedure, we follow the line of reasoning developed by Kim and Rudd [
15], Kim [
16], and Kim and Desclaux [
17] in their Binary-encounter-dipole model for electron-impact ionization and excitation of atomic targets.
It was shown that the total cross-sections produced by the plane-wave Born (PWB) approximation, a first-order theory that considers plane waves for the incoming and scattered projectile and approximate wave functions for non-hydrogenic targets, do not correctly describe the experimental data in the low energy range (impact energies lower than 100 eV) [
15,
16]. Such deficiency has been ascribed to the lack of electron exchange with the target electrons, the distortion of plane waves in the vicinity of the atomic target, or the polarization of the target due to the presence of the incident electron. To correct these deficiencies Kim and Rudd introduced a scaling method, called BEB-scaling [
15]. In this scaling, the TCS for direct ionization of an electron in a bound orbital by an electron of incident energy
T is given by
Here T and B are the impact energy and the ionization potential, respectively, R is the Ry constant, t = T/B, and u = U/B, and U is the orbital kinetic energy.
Similarly, a scaling for electron impact excitation, which was denominated BE and consists in multiplying the PWB total cross-section by an energy-dependent factor, was presented in Ref. [
16]:
In this case,
E is the excitation energy. In Ref. [
17], Kim and Desclaux showed that these semiempirical scaling used together are able to reproduce Brook’s data.
Following the spirit of those studies, we multiply our CTMC total ionization cross-sections by energy-dependent factors,
This model, which we denominate B-CTMC, is not designed to provide the ultimate description of the data, but to help detect the limitations of our classical model instead. In this sense, for our purposes, we have explicitly omitted the
E term in the denominator of the scaling factors. As can be seen in
Figure 4, the agreement with the experimental data remarkably improves. While the B-CTMC results converge to the CTMC results at large impact energies, the peak position and the slope of the total cross-section near the threshold, are now in good agreement with the data. We have also added the recent scaled Distorted Wave Born Approximation (DWBA) results from Jonauskas (Ref. [
18]), in which the receding projectile interaction with the target is either described by a potential that considers the neutral target, just as in the initial state or the already ionized target.
These results suggest that for the electron impact case, the CTMC is not providing an accurate description of the collision process at low impact energies. However, the possible issues indicated by Kim and Rudd for the PWB approximation do not seem to apply in this case. This is due to the fact that the CTMC model is a non-perturbative method, and explicitly considers the projectile interaction with the target electron and the ionic core throughout the whole collision process. In addition, electron exchange events take place and gain relevance for decreasing impact energies. To further illustrate this point, in
Figure 5 we show some classical trajectories that lead to ionization and electron exchange for the 2p and 2s orbitals.
Two points are worth mentioning at this stage. On the one hand, the classical atom has no minimum energy for the ground state. This feature becomes especially relevant at low impact energies when the target electron can accommodate a more deeply bound state as the projectile approaches, leading to a transfer of kinetic energy to the projectile [
14].
On the other hand, as the impact energy decreases the chance of having a transient negative ion increases. The negative ion of carbon has two bound states, the ground state 2s
22p
3 4S
0 and the excited state 2s
22p
3 2D
0 [
19]. From the classical point of view, and in contrast to the prediction of quantum mechanics, atoms or ions with two or more bound electrons autoionize without the need for any external agent. Different alternatives have been proposed to avoid this spurious autoionization in target systems, like the Bohr atom [
20], the neglect of the e-e interaction [
21], the backward-forward propagation scheme [
22], and the introduction of momentum-dependent stabilizing potential terms to the Hamiltonians [
23,
24,
25,
26,
27,
28,
29]. In contrast, in this case, the duration of the transient C
− ion and the fraction of such events is expected to depend on the projectile electron energy. Therefore, it is hard to infer firsthand the role that the e-e interaction might have on the electron emission process. Hence, at this stage we inspect the transient lifetime of C
− ions, that is the lapse of time that the three-body system conforms to a C
− ion during collision events. In
Figure 6 we show the obtained results for the 2p and 2s states at an impact energy of 30 eV, which reveal that these negative ions last for less than 1 a.u. of time.
Having determined the time scale, and with the aim of detecting whether or not this particular issue plays a role in the electron emission process, we decided to turn off the e-e interaction during the lapse that the two electrons were bound to the ionic core. In this drastic approach, the projectile and the target electron are only subjected to the field of the target nucleus during the transient lifetime of the C
− ion. The obtained results are shown in
Figure 7. In
Figure 7a we observe that the obtained total cross-section is now in much better agreement with the data, in terms of magnitude as well as peaking position. In addition,
Figure 7b shows that the fraction of collision events in which transient C
− ions are found, increases with decreasing impact energy as expected. Moreover, the number of such events is found to be more relevant for the 2s orbital.
4. Discussion
Accurate data on electron carbon ionizing collisions are needed for diagnostic purposes in the power fusion reaction program. Although highly reliable numerically intensive methods have been used for this collision system during the past decade, these studies have been limited to impact energies lower than 100 eV. This highlights the need for reliable theoretical data at larger impact energies.
In this work, the classical description of the electron impact ionization of carbon has been analyzed as provided by the three-body CTMC method starting from the threshold region up to an impact energy of 1 keV.
Present results suggest that in its standard microcanonical formulation, and within the standard addition rule, the CTMC method is in perfect agreement with the results provided by the generalized oscillator strength formulation of the Born approximation of McGuire at impact energies greater than 100 eV and clearly support the data of Brook. In contrast, it tends to overestimate the electron emission at low impact energies. Furthermore, it predicts a clear shift of the peak of the total ionization cross-section towards lower impact energies.
A semi-empirical model following the spirit of those developed by Kim and Rudd clearly improved the agreement of the CTMC results with the data, suggesting that some physical aspects of the problem might not be accurately described by a classical description, especially at low impact energies.
Our inspection of the collision dynamics revealed the existence of transient C
− ionic states that the presently used Hamiltonian cannot properly describe. In order to test the potential relevance of these C
− ionic states in the electron emission process, the e-e interaction was turned off during their lapse of existence. This strategy led to a clear improvement in our CTMC predictions in the peak region of the TCS and opens a line of study that needs to be explored. Future studies should concentrate on the threshold region, usually denominated the Wannier region, where classical mechanics for the two-electron continuum is expected to apply and lead to a power law dependence of the TCS as a function of the impact energy [
30].