Abstract

One of the goals in studies of double ionization (DI) of simple atoms by electron or ion impact is to elucidate and assess the different mechanisms that lead to this atomic process. In this work we present an attempt to model the mechanisms beyond the first order in DI of helium by highly charged projectiles. To this end we employ the continuum distorted wave-eikonal initial state (CDW-EIS) formalism joint with a Monte Carlo event generator (MCEG). The MCEG allows us to generate theoretical event files that represent the counterpart of the data obtained from a kinematically complete experiment. Starting from these event files, a new data analysis tool used to contrast theory and experiment in DI, the four-body Dalitz plots, is easily produced. The higher order mechanisms are simulated by considering DI as a sequential process: a single ionization of a helium atom as a first step and another single ionization of a single-charged helium ion as a second step. Some of the features in the experimental data are very well reproduced by these simulations.

1. Introduction

Among the atomic breakup processes, double ionization (DI) of simple atoms by charged particles represents a particularly interesting case to study the few body dynamics. The main difference between DI and its simplest counterpart, single ionization (SI), is the introduction of an ingredient that plays an important role under certain kinematics conditions: the electron-electron correlation. Additionally, DI proceeds by different pathways, and one of the goals of the theory is to assess the role of each of them for a given set of the parameters. It is common to consider that the total cross section for DI is mainly dominated by three different mechanisms. Two of them are first-order mechanisms and the third one is a higher-order mechanism. The two first-order processes are the two-step-one projectile-electron interaction (TS-1) and the so-called shake-off (SO). In TS-1 only one electron is released through a direct interaction with the incoming projectile. This ejected electron then collides in the second step with the other electron liberating it to the continuum. The SO is very similar to TS-1 in that only one electron is ejected by a direct interaction with the projectile. However, here the second electron gets ejected through rearrangement of the electronic wave function adjusting to a new potential resulting from the ejection of the first electron. Both TS-1 and SO mechanisms are first-order processes in the projectile-target atom interaction and in both some form of electron-electron correlation is required for DI to occur. For simplicity in the following we refer to both of them as TS-1. Finally, in the third mechanism, dubbed two-step-two projectile-electron interaction (TS-2), both electrons are ejected through two independent interactions with the projectile. This is a higher-order process which does not require neither does it rule out electron-electron correlation [1].

Recently, it was demonstrated that a pure form of TS-1 is never realized [2]. Previously, it was assumed that TS-1 is dominant in very fast collisions (except for highly charged ions) and TS-2 for slow collisions and/or highly charged ion impact. However, later it was demonstrated that even for very fast proton impact the most important DI mechanism is a hybrid process between TS-1 and TS-2. This channel, called TS-1-EL, can be viewed as a TS-1 process followed by elastic scattering between the projectile and one of the electrons already ejected to the continuum.

In recent years there has been a steady advance of experimental techniques in atomic physics. Nowadays kinematically complete experiments on diverse fragmentation processes are feasible, for example, single and double ionizations by charged particle impact [3–8] or mutual ionization of both collision partners [9]. Especially Cold Target Recoil Ion Momentum Spectroscopy (COLTRIMS), also known as Reaction Microscopes, has opened a new era of atomic collision experiments [10]. The data obtained from such experiments are very complex with several measured quantities being afflicted with experimental resolution and acceptance limitations. To fully account for these experimental restrictions in theoretical calculations, Monte Carlo Event Generators (MCEGs) are a very powerful tool.

In atomic collisions, MCEG have been developed and applied only recently [11]. In further developments, MCEG have been extended to treat DI of helium and also were used in order to analyze the DI with a new and powerful analysis tool: the four-body Dalitz plots (4D) [12]. In a kinematically complete experiment, the cross sections are usually extracted from a data file known as event file. Such file contains all three momentum components of all collision fragments recorded for a large number of fragmentation events. Theoretical approaches usually do not provide such an event file obtained in experiment. Rather, the probability that an event with a certain combination of momentum components of the collision fragments occurs is computed directly. An important shortcoming of this approach is that this makes a direct comparison between measured data and theory accounting for the experimental conditions difficult. As long as the cross sections depend only on one or two quantities which are afflicted by experimental resolution, this can be accounted for by analytically convoluting the calculation with the resolution. However, fully differential cross sections (FDCSs) can depend on a large number of parameters and a corresponding analytic convolution is simply not feasible. A convolution by numerical integration is usually not a realistic option either because it would make the computation prohibitively time-consuming. These problems can be avoided if an event file with the same structure as the one obtained in the experiment is generated using a certain theoretical approach. The experimental resolution can then be accounted for by adding event-by-event a random number, distributed according to the spectrometer function of the apparatus (typically a Gauss distribution), to each momentum component [11].

Another important advantage of using a theoretical event file is that additional physics, which may be problematic to accurately account for in the calculation, can be added retroactively. For example, the interaction between the nuclei of both collision partners (NN interaction) was often neglected in previous studies. However, convoluting the results of theory with classical elastic scattering between the nuclei using the MCEG technique was recently demonstrated to yield much better agreement between experiment and theory for single ionization in fast C6++He impact [13] and intermediate energy 𝑝+H collisions [14].

Finally, one of the most important benefits of employing theoretical event files is the capability of computing cross sections which are numerically too intensive using conventional methods. Usually, theoretical cross sections are obtained by appropriate integration of the FDCS. For example, double ionization entails 8 independent momentum components (12 components of the 4 final-state particles minus 4 kinematic conservation laws). The computation of a single differential cross section thus requires a sevenfold integration. For sophisticated theoretical models to a large extent this integration often has to be done numerically. However, the symmetry of the problem can considerably reduce the dimensionality of the integration so that the computation of the single differential cross section is not too problematic. On the other hand, if the cross section of interest only entails limited symmetry, or none at all, the numeric efforts required by the integration can quickly make the calculation unfeasible. Using an event file, in contrast, any cross section measured in the experiment can be obtained simply by sorting the histograms precisely the same way as done for the measured data (for computational details see [15]). A prominent example is the so-called four-particle Dalitz (4-D) plots [12], as mentioned above.

The successful usage of these new tools requires reliable and robust theories to account for the underlying ionization dynamics. Over the years continuum distorted wave-based-formalisms have been used extensively for single ionization by electron and ion impact (see, e.g., [16] for a recent review), but, to the best of our knowledge, there is no such approach for double ionization due to the difficulties present to develop a formal theory beyond the First Born approximation. Although the TS-2 simulations do not represent a rigorous continuum distorted wave theory for DI, we demonstrate that nevertheless important qualitative information can be extracted from a comparison with measured data.

The article is organized as follows. The next section describes briefly the theory behind the TS-2 simulations and the computational tools used to produce them. In Section 3 we present a comparison between our theory and experiment. Specifically, we investigate the sensitivity of the theoretical results with respect to the simulation of elastic scattering between the heavy particles. Conclusions and an outlook are offered in Section 4. Atomic units are used throughout the article unless otherwise stated.

2. The Monte Carlo Event Generator (MCEG) and the TS-2 Simulations

As a first attempt to simulate the TS-2 mechanisms we employ the Monte Carlo Event Generator (MCEG) [11] along with the CDW-EIS theory [17]. The MCEG allows us to generate a theoretical event file that will be later used to produce the 4-body Dalitz (4-D) plots [2, 12, 18, 19]. The outcome of the MCEG is a data file, based on a calculation of FDCS for the process of interest (see [15] for computational details). These FDCSs are computed using the CDW-EIS approach which takes into account higher-order contributions to the ejection of each electron.

In double ionization (DI), the number of momentum components per event contained in the event file depends both on the theoretical scheme we apply and the type of processes we want to model, namely, first-order processes or higher-order mechanisms. We model the TS-2 process as a combination of two independent SI steps. Consequently the event file contains the independent momentum components of both SI steps, namely, (π‘˜1,π‘₯,π‘˜1,𝑦,π‘˜1,𝑧,π‘ž1,π‘₯) for the first SI step, that is, the single ionization of a neutral helium atom, and (π‘˜2,π‘₯,π‘˜2,𝑦,π‘˜2,𝑧,π‘ž2,π‘₯) for the second step, that is, the single ionization now of an He+ ion. Here (π‘˜π‘–,π‘₯,π‘˜π‘–,𝑦,π‘˜π‘–,𝑧), with 𝑖=1,2, are the cartesian components of the electron momentum and (π‘žπ‘–,π‘₯,0) the single component of the transverse momentum transfer for each SI step.

In order to produce the 4-D plots for DI of helium by ion impact, typically a few 100,000 events are required. A summary of the procedure followed to generate a theoretical event file will be presented below. The described cycle is repeatedly executed until the number of events needed is reached. (i)A set of ten random numbers is chosen, eight of which represent the final-state momentum components of the fragments for which the FDCS will be calculated, that is, π‘˜1,π‘₯,π‘˜1,𝑦,π‘˜1,𝑧,π‘ž1,π‘₯,π‘˜2,π‘₯,π‘˜2,𝑦,π‘˜2,𝑧,π‘ž2,π‘₯. The ranges of each electron momentum component and of the momentum transfer components, that is, the intervals of the corresponding random number distributions, are established in advance according to the experimental conditions present in the particular DI process. The ninth and tenth quantities, 𝑒 and 𝑀, are random numbers uniformly distributed in the interval [0,1].(ii)The FDCSs are calculated at the randomly selected point. For the case of TS-2 simulation, we employ in each SI step the CDW-EIS model of Crothers and McCann [17, 20]. This model represents one of the most successful theories for SI of atoms by single and multiply charged ion impact in the range of intermediate to large projectile energies [21]. Since our aim is to model higher-order processes, the CDW-EIS appears to be the most adequate choice.(iii)Once the FDCSs for each SI steps are obtained from the CDW-EIS model, the randomly selected event is considered valid and stored if both of the following conditions are fulfilled𝑒<ξ‚πœŽSI,𝑀<ξ‚πœŽSI+,(1) otherwise the event is discarded. ξ‚πœŽSI (ξ‚πœŽSI+) corresponds to the normalized cross sections for SI of neutral helium (He+ ion) calculated using ξ‚πœŽSI=𝜎SI/𝜎maxSI (ξ‚πœŽSI+=𝜎SI+/𝜎maxSI+) where 𝜎maxSI (𝜎maxSI+) is the maximum value of the cross section within the selected phase space. As the cross section for both SI steps drop rapidly towards larger values of electron momenta and momentum transfers, the simple rejection method described above becomes quite inefficient and consequently many events are rejected so that the cycles (i)–(iii) have to be repeated hundreds of millions or billions of times. Fortunately, using a simple parallelization scheme, we can take advantage of the independence of this kind of Monte Carlo-based calculations, in the sense that the momentum space can be sampled independently and in this way a concurrent computation with high efficiency is easily implementable (see [15] for further details).

The FDCSs of each SI, that is,. 𝜎SI and 𝜎SI+, are calculated using the CDW-EIS approach of Crothers and McCann [17] and can be written as𝜎SI,SI+=d5𝜎d𝐀dπͺβŸ‚βˆ||π‘‡βˆ’CDW-EIS𝑖𝑓||2π›Ώξ€·πΈπ‘“βˆ’πΈπ‘–ξ€Έ,(2) where we have ignored the constants because they cancel out during the normalization procedure explained above. In (2), the electron momentum is 𝐀=(π‘˜π‘–,π‘₯,π‘˜π‘–,𝑦,π‘˜π‘–,𝑧) and the transverse momentum transfer πͺβŸ‚=(π‘žπ‘–,π‘₯,0) with 𝑖=1,2 for each SI step. We employ nonorthogonal Jacobi coordinates (𝐫𝑃,𝐫𝑇) to describe the collision process. These coordinates represent the position of the active electron of the neutral helium (or the single electron of the He+ ion) with respect to the projectile (𝐫𝑃) and the target ion (𝐫𝑇), respectively. 𝐑𝑇 is also needed, representing the position of the incoming ion with respect to the CM of the subsystem e-T. Neglecting terms of orders 1/𝑀𝑇 and 1/𝑀𝑃, where 𝑀𝑇 is the mass of the target ion nucleus and 𝑀𝑃 is the corresponding mass of the incident heavy ion, we can write 𝐑𝑇=π«π‘‡βˆ’π«π‘ƒ and in this way the transition matrix π‘‡βˆ’CDW-EIS𝑖𝑓 is separable in the variables 𝐫𝑃 and 𝐫𝑇 as it will become apparent below.

Within prior CDW-EIS, the transition amplitude can be computed asπ‘‡βˆ’CDW-EIS𝑖𝑓=ξ‚¬πœ’βˆ’CDW𝑓||π‘Šπ‘–||πœ’+EIS𝑖,(3) where the initial (final) state distorted wave πœ’+𝑖 (πœ’βˆ’π‘“) is an approximation to the initial (final) state which satisfies outgoing-wave (+) (incoming-wave (βˆ’)) boundary conditions. For the initial state, the asymptotic form of the Coulomb distortion (eikonal phase) is used in the electron-projectile interaction. To model the initial bound-state wavefunction, we employ a semianalytical Roothaan-Hartree-Fock for the description of the neutral helium atom [22] and a hydrogenic 1s state with 𝑍=2 for the He+ ion. Consequently, the full initial state πœ’+EIS𝑖 of (3) can be written asπœ’+EIS𝑖=(2πœ‹)βˆ’3/2ξ€·expiπŠπ‘–β‹…π‘π‘‡ξ€Έπœ“He,He+𝐫𝑇ℰ+𝑣𝐫𝑃,(4) where β„°+𝑣(𝐫𝑃) isβ„°+𝑣𝐫𝑃𝑍=expβˆ’i𝑃𝑣lnπ‘£π‘Ÿπ‘ƒβˆ’π―β‹…π«π‘ƒξ€Έξ‚Ά,(5) and πœ“He,He+(𝐫𝑇) is the initial bound state of the neutral helium atom or the He+ ion.

The final state wavefunction is cast into the form [23]πœ’βˆ’CDW𝑓=(2πœ‹)βˆ’3/2ξ€·expiπŠπ‘“β‹…π‘π‘‡ξ€Έπœ’βˆ’π‘‡ξ€·π«π‘‡ξ€ΈπΆβˆ’π‘ƒξ€·π«π‘ƒξ€Έ,(6) where πΆβˆ’π‘ƒ represents the Coulomb distortion of the ejected electron wave function due to the projectile.πΆβˆ’π‘ƒξ€·π«π‘ƒξ€Έξ€·πœˆ=𝑁𝑃1F1ξ€·βˆ’iπœˆπ‘ƒ,1,βˆ’iπ‘˜π‘ƒπ‘Ÿπ‘ƒβˆ’i𝐀𝑃⋅𝐫𝑃,(7) where πœˆπ‘ƒ=𝑍𝑃/π‘˜π‘ƒ is the Sommerfeld parameter, 𝐀𝑃 is the relative momentum of the e-P subsystem and 𝑁(πœˆπ‘ƒ) is the usual Coulomb factorπ‘ξ€·πœˆπ‘ƒξ€Έξ€·=Ξ“1βˆ’iπœˆπ‘ƒξ€Έξ‚€expπœ‹πœˆπ‘ƒ2.(8) On the other hand, πœ’βˆ’π‘‡(𝐫𝑇) is the wave function for the ejected electron in the field of the target residual ion.πœ’βˆ’π‘‡ξ€·π«π‘‡ξ€Έ=(2πœ‹)βˆ’3/2ξ€·expiπ€π‘‡β‹…π«π‘‡ξ€Έπ‘ξ€·πœˆπ‘‡ξ€ΈΓ—1F1ξ€·βˆ’iπœˆπ‘‡,1,βˆ’iπ‘˜π‘‡π‘Ÿπ‘‡βˆ’i𝐀𝑇⋅𝐫𝑇(9) with πœˆπ‘‡=𝑍𝑇/π‘˜π‘‡ and where 𝐀𝑇 is the relative momentum of the e-T subsystem. For the neutral helium atom, we use an effective charge 𝑍𝑇=1.69 in order to consider the partial screening of the passive electron. On the other hand, for the SI of the He+ ion, 𝑍𝑇=2 is employed.

The perturbation potential π‘Šπ‘– in (3) is defined by [24]ξ€·π»π‘–βˆ’πΈπ‘–ξ€Έπœ’+𝑖=π‘Šπ‘–πœ’+𝑖,(10) where 𝐻𝑖 is the full electronic initial Hamiltonian (neglecting the total CM motion) and 𝐸𝑖 is the total initial energy of the system in the CM frame, respectively. π‘Šπ‘– is composed by two differential operators that can be writtenπ‘Šπ‘–=12βˆ‡2π«π‘ƒβˆ’βˆ‡π«π‘‡β‹…βˆ‡π«π‘ƒ.(11) In the computation of the FDCS (2), we have neglected the interaction between the heavy nuclei, the so-called NN interaction. This interaction, however, is incorporated later in the processing of the event file (see, e.g., [13] for details) and is taken into account in our simulations.

3. Results and Discussions

In Figures 1–3, 4-D plots are shown for DI in 158 MeV Au33++He collisions. Similar plots were discussed in [19], but here we place emphasis on an analysis of the sensitivity of the theoretical results with respect to the simulation of elastic scattering between the heavy particles. Details on how to read 4-D plots were presented previously [12, 15, 18]. In short, a tetrahedral coordinate system is employed where each tetrahedron plane represents one of the four final-state fragments of the collision (front and bottom planes: ejected electrons; backplane: target nucleus; right plane: momentum transfer from projectile to target atom). The distance of a given data point to the four planes represents a set of the relative squared momenta of the final-state particles πœ‹π‘–=𝑝2𝑖/βˆ‘π‘2𝑗. In Figure 1, the comparison between the experimental (panel a) and theoretical 4-D plots, calculated within the TS-2 CDW-EIS model (panel b), shows only marginal, qualitative agreement as far as the boundary of the area containing any data points at all is concerned. However, the intensity distribution within these boundaries is poorly reproduced. This lack of agreement is not entirely surprising considering the very large perturbation parameter πœ‚ (projectile charge to speed ratio) of 5.8. For such extreme πœ‚, the cross sections are dominated by higher-order contributions which are extremely difficult to accurately account for theoretically. It is not even clear whether a perturbative expansion still converges. As a result, the calculations become very sensitive to the details of the description of the various interactions. A simple example is illustrated in panels (b) through (d) in Figure 1. There, different assumptions were made for 𝑏𝑃(𝑏) for double ionization (where 𝑏 is the impact parameter and 𝑃(𝑏) is the double ionization probability), which enter in the convolution of the FDCS with elastic scattering between the heavy particles [13]. In panel (b), 𝑏𝑃(𝑏) was assumed to maximize at 𝑏=10 a.u., in panel (c) at 𝑏=6 a.u., and in panel (d) the FDCSs were not convoluted at all with elastic scattering. The significant qualitative changes in the corresponding 4-D plots are quite apparent.

In Figure 2 experimental and theoretical 4-D plots are shown for the same collision system and in the same order as in Figure 1. However, this time only the transverse momentum components of each particle in the scattering plane (spanned by the initial and final projectile momenta) were used in the computation of the πœ‹π‘–. Remarkably, in this case the experimental data are reasonably well reproduced by theory. This improved agreement goes along with a much reduced sensitivity of the calculation on the assumptions regarding 𝑏𝑃(𝑏). Furthermore, the experimental 4-D plots for the transverse components are very similar to those for DI in 6 MeV 𝑝+ He collisions, corresponding to a very small πœ‚ of about 0.06. The insensitivity to both 𝑏𝑃(𝑏) and πœ‚ suggests that the 4-D plots for the transverse components are not very well suited to study the DI dynamics. On the other hand, the observation that theory is capable of reproducing the experimental plots, regardless of πœ‚, shows that in the transverse direction the 4-D plots offer a good test of the theoretical codes for errors.

In contrast to the 4-D plots for the transverse components, those for the longitudinal components are very sensitive to the collision dynamics. Here, the plots look qualitatively very different for small and large πœ‚ [19]. After testing the codes using the transverse components, a comparison between experiment and theory using the longitudinal components should therefore provide detailed information about the DI dynamics. This comparison is shown in Figure 3 (experiment: panel (a); TS-2 CDW-EIS: panel (b)). It should be noted that in the longitudinal direction elastic scattering between the heavy particles has essentially no effect because the longitudinal momentum transfer is to a very good approximation directly proportional to the inelasticity. Qualitatively, the data are very well reproduced by the TS-2 CDW-EIS calculation. In contrast, a TS-1 calculation based on the First Born Approximation (FBA) [18], shown in panel (c) of Figure 3, is in rather poor agreement with experiment. The TS-2 CDW-EIS and TS-1 calculations shown in Figure 3 differ fundamentally from each other in two aspects: first, both calculations obviously treat different DI mechanisms. Second, the TS-2 CDW-EIS model treats already the ejection of a single electron to higher order, while the TS-1 model is a first-order calculation even for the complete DI process. A comparison between the TS-1 results and a TS-2 calculation in which the ejection of each electron is treated within the FBA [2] (TS-2 FBA, panel d of Figure 3) shows that different DI mechanisms lead to nonnegligible differences in the longitudinal 4-D plots. However, earlier we demonstrated in detail that most of the characteristic features in the experimental plot can be explained by higher-order effects in the ejection of a single electron. More specifically, the postcollision interaction (PCI) between the outgoing projectile and one or both of the ejected electrons plays a very important role [2].

4. Conclusions and Outlook

We have presented a first attempt to model higher-order processes in double ionization of helium by highly charged particles. The approach treats double ionization as a sequential process where the projectile first ionizes a neutral helium atom and subsequently removes the remaining electron of the He+ ion. Each step on this sequence is a single ionization process and can be modeled using the CDW-EIS theory. We have compared our theoretical results on 4-body Dalitz plots with the experimental data. These plots are generated using a theoretical event file, that was computed using the Monte Carlo Event Generator modified accordingly to consider double ionization as a sort of convolution of two independent single ionization steps. Good agreement with the experiments for the case of double ionization of helium by impact of 158 MeV of Au33+ projectiles was observed in the 4-D plots generated for the longitudinal and transverse components. We have also proven that the 4-D plots produced using the transverse components are quite insensitive to the charge to velocity ratio and in this way allow us to validate our different theoretical approaches. Among the next steps in order to overcome the remaining discrepancies, we plan to develop and implement a more formal theory for double ionization where both first and higher-order processes are treated coherently.