Net electron capture in collisions of multiply charged projectiles with biologically relevant molecules

A model for the description of proton collisions from molecules composed of atoms such as hydrogen, carbon, nitrogen, oxygen and phosphorus (H, C, N, O, P) was recently extended to treat collisions with multiply charged ions with a focus on net ionization. Here we complement the work by focusing on net capture. The ion-atom collisions are computed using the two-center basis generator method. The atomic net capture cross sections are then used to assemble two models for ion-molecule collisions: an independent atom model (IAM) based on the Bragg additivity rule (labeled IAM-AR), and also the so-called pixel-counting method (IAM-PCM) which introduces dependence on the orientation of the molecule during impact. The IAM-PCM leads to significantly reduced capture cross sections relative to IAM-AR at low energies, since it takes into account the overlap of effective atomic cross sectional areas. We compare our results with available experimental and other theoretical data focusing on water vapor (H2O), methane (CH4) and uracil (C4H4N2O2). For the water molecule target we also provide results from a classical-trajectory Monte Carlo approach that includes dynamical screening effects on projectile and target. For small molecules dominated by a many-electron atom, such as carbon in methane, or oxygen in water we find a saturation phenomenon for higher projectile charges (Q = 3) and low energies, where the net capture cross section for the molecule is dominated by the net cross section for the many-electron atom, and the net capture cross section is not proportional to the total number of valence electrons.


I. INTRODUCTION
Collisions of multiply charged ions with biologically relevant molecules are recognized as being important for future developments in radiation medicine and related fields. There are numerous experimental efforts which focus largely on differential electron emission, but total ionization, electron capture, and excitation are also relevant since they contribute to energy deposition, and stopping power. A recent example is the study of net ionization in C 6+ − CH 4 collisions [1].
On the theoretical side one finds several methods that attempt to explain the experimental data. Often provided alongside with the experimental work is the continuum distorted wave with eikonal initial state (CDW-EIS) method. Differential electron emission can be obtained also directly from a classical-trajectory Monte Carlo (CTMC) method, which for water molecule targets was enhanced recently to include dynamical screening effects [2]. The CTMC approach takes the (frozen) molecular orientation into account during the collision, and total ionization cross sections, as well as some charge-state correlated cross sections have been compared to experimental data [3].
Another approach is to use collision information obtained for the atomic constituents and to combine them into molecular cross sections, most notably independent atom models (IAM), which either follow the simple Bragg additivity rule, or more sophisticated versions that take the molecular structure of the target into account, and allow for the fact that the effective cross section should be reduced due to overlap effects. This pixel counting method (IAM-PCM) tested originally for proton impact [4][5][6][7] was used recently to investigate scaling behavior of net ionization as a function of projectile charge [8].
The CDW-EIS work also relies on ion-atom collision calculations, and includes some molecular effects on the basis of a Mulliken population analysis. For the ionization problem we found that CDW-EIS and CTMC calculations with frozen potentials generally are close to IAM-AR results, while the CTMC calculations with dynamical screening [3] are closer to IAM-PCM absolute cross sections [8], which are lower in the vicinity of the ionization maximum and merge with IAM-AR results only at very high energies when the projectile charge is high, i.e., Q > 1. The reduction in ionization cross sections (whether in the IAM-PCM results versus IAM-AR or CTMC with versus without dynamical screening) can easily reach a factor of two.

2
In the present paper we focus on an analogous problem at lower collision energies, namely the net capture problem. Here the ion-atom cross sections grow quickly with projectile charge, and the overlap and orientation problem again becomes important for the molecular targets. The experimental situation is more scarce in the case of capture as compared to ionization, and therefore theoretical support is even more important.
The paper is organized as follows. In Sect. II we introduce the theoretical basis for the current work. Sect. II A presents new results for two-center basis generator method (TC-BGM) ion-atom calculations for He 2+ and Li 3+ projectiles; for atomic hydrogen targets the results are compared with theory and experiment. In Sect. II B we provide a summary of the IAM-PCM methodology and briefly describe the CTMC time-dependent mean field approach which is applied to the water molecule target. Sect. III serves to provide a detailed comparison with experimental data and a selection of other theoretical work: in Sect. III A we present results for water, where we also compare with CTMC results with and without dynamical screening; in Sect. III B we present results for methane for which we demonstrate a saturation effect in the case of projectile charges Q = 3, and Sect. III C contains results for uracil. The paper ends with a few concluding remarks in Sect. IV. Atomic units, characterized by = m e = e = 4π 0 = 1, are used unless otherwise stated.

A. Ion-atom collisions
The independent atom models rely on state-of-the-art ion-atom collision calculations. In Fig. 1 we show the results of the two-center basis generator method (TC-BGM) calculations for atomic targets that form the constituents of biologically relevant molecules (hydrogen, carbon, nitrogen, oxygen, and phosphorus). Experimental verification is only available for atomic hydrogen, as is theoretical confirmation by other state-of-the-art methods, namely the convergent close coupling (CCC) approach [9,10], and a two-center atomic orbital expansion method based on Gaussian-type orbitals [11].
The present TC-BGM ion-atom calculations are as described in Ref. [8]: for proton impact a projectile potential W P hierarchy is obtained on the basis of explicitly including shells up to principal quantum number n = 4 on projectile and target, for Q = 2 projectiles we used For atomic hydrogen targets and proton and He 2+ projectiles the results are compared with the convergent close coupling calculation of Ref. [9,10] respectively. For Li 3+ impact we compare with data from Ref. [11]. For proton impact the experimental data are from Ref. [12][13][14][15], for alpha particle impact from Ref. [15][16][17], and for Li 3+ − H collisions from Ref. [18].
n P = 6 and n T = 5, while for Q = 3 the explicitly included basis was expanded to n P = 7 and n T = 5. These explicitly included states are complemented by states that represent the continuum using the BGM approach [5,19].
For the many-electron atoms the theoretical modelling is at the level of exchange-only density functional theory using the optimized potential method. Therefore, electron correlation effects are not included in these calculations [20,21].
The results for proton and alpha particle impact on atomic hydrogen display excellent agreement with experiments and with CCC theory over several orders of magnitude. For projectile charge Q = 3, i.e., Li 3+ projectiles, our results are slightly higher than the experimental data, but are in very good agreement with the recently reported calculations of Ref. [11]. Not shown in Fig. 1 for proton and He 2+ impact are the calculations of Ref. [11].
They agree very well the TC-BGM and CCC results for protons and are slightly lower (at the 10 % level for energies 25 − 100 keV/amu) for He 2+ projectiles.
The structure of these cross sections is simple. Notable is the increase in net capture as one reduces the collision energy towards 10 keV/amu, with cross section values short of 8 Å 2 for protons, 12 Å 2 for Q = 2, and a value of 20 Å 2 reached already at 30 keV/amu for bare lithium projectiles.
For the atoms containing more than one electron capture from the outer shells becomes large even for intermediate energies with patterns that are not totally straightforward when one compares Q = 1, 2, 3. Capture from phosphorus reaches 100 Å 2 for Q = 3 at low energies, and electron capture from carbon atoms is also very strong. This provides a background for interesting effects when combining these cross sections to make predictions for molecular targets.
For phosphorus atoms there are notable structures in the capture cross sections at higher energies which are related to shell effects. For collision energies below 200 keV/amu capture from the M-shell dominates, but at higher energies capture from the L-shell becomes more important. This was investigated, e.g., in Ref. [19] for p − Na collisions.

B. Ion-molecule collisions
The present work reports on calculations of several theoretical models. On the one hand we compare two independent atom models: (i) in the naive case, namely the additivity rule based model (IAM-AR) the atomic capture cross sections described in Section II A are simply added together and completely ignore molecular structure; (ii) in the more sophisticated pixel counting method (IAM-PCM) molecular structure is introduced by a geometric procedure, which will be briefly reviewed below. For the water molecule target we also report results from two classical-trajectory models.
The IAM-PCM has been described and illustrated in detail in previous papers, most notably Refs. [4,7]. Ref. [7] highlights the role played by the contributions from the time evolution of occupied orbitals from different atomic species to the net cross sections. The model is compared there to a methodology employed by the community that uses the CDW-EIS approach, which also uses ion-atom calculations, and incorporates molecular eigenenergies on the basis of quantum molecular structure calculations while employing complete neglect 5 of differential overlap (CNDO), i.e., a Mulliken population analysis [22,23].
In essence the IAM-PCM approach is based on an interpretation of the atomic net cross sections as geometrical areas, e.g., areas that correspond to the net capture (or net ionization) cross section. Rather than summing up all these areas (as assumed by the additivity rule, i.e., the IAM-AR) a pixel counting method is used to measure the effective cross-sectional area that emerges when one eliminates overlaps between cross sections from different atoms encountered by the projectile. Therefore, the effective molecular cross sectional area is defined as a function of molecular orientation. The latter is chosen randomly, and is sampled to obtain converged cross sections. A critical discussion of the merits of this procedure (potential emphasis on atoms encountered first) is found in Ref. [5], where it is argued that for net cross sections the method should certainly be appropriate. An illustration of the method is given below in Section III B where we observe a saturation behavior of the cross section for Q ≥ 3 projectiles colliding with methane (CH 4 ).
For a detailed description of the other method for which we show results for water molecule vapor targets, namely the CTMC model we also refer to previous literature. Calculations with frozen target potential on the basis of a three-center potential are reported in Refs. [24,25].
This effective potential used in the classical statistical ensemble simulation is drawing information from quantum structure calculations. It yields accepted values for the orbital energies and in addition to CTMC calculations was also used for numerical grid solutions of the time-dependent Schrödinger equation [26]. An extension of this model which can be considered a semiclassical ( = 0) approximation to the quantum problem is presented in Ref. [2] where the potential parameters were allowed to vary dynamically as a function of the average (net) ionization state of the water molecule. This resulted in a substantial reduction of the net ionization cross section (by up to a factor of two) compared to the static potential model results. A further extension was carried out recently where the projectile potential was also allowed to be varied as a function of the average (net) charge state of the projectile [3].
Both the static-potential CTMC and the dynamically screened model are sensitive to the orientation of the molecule during the collision. Therefore, the comparison with IAM-PCM and IAM-AR calculations is of great interest. Concerning capture data there is one problem that needs to be addressed in the CTMC approach: for high energies capture from inner shells in target atoms (or molecules) becomes problematic in a classical-trajectory model, because it becomes possible to capture into orbits well below the allowed 1s-level, 6 i.e., orbits with binding energies larger than Q 2 /2 begin to occur. These classically allowed capture contributions need to be removed from the analysis. For proton impact (Q = 1) this correction becomes significant at impact energies of 200 keV/amu and higher, for Q > 1 this point moves to higher energies. The correction procedure is based on a prescription to associate principal quantum numbers n with energy ranges [27].

III. COMPARISON WITH EXPERIMENTS
A. Collisions with water vapor (H 2 O) In Fig. 2 we compare the theoretical model results for proton impact on water vapor with experiment. We note that the IAM-PCM results are corrected compared to those shown in Refs. [4,5] (the wrong bond length was applied in those calculations), and they agree now well within error bars at low energies with the experimental data of Rudd et al. [28]. The reaching a factor-of-two discrepancy at 100 keV/amu. This shortfall of the experimental data goes hand-in-hand with an observed shortfall in double-ionization contributions in this energy range as compared to TC-BGM calculations for pure ionization (cf. Fig. 7 of Ref. [32]) and transfer ionization (cf. Fig. 8 of Ref. [32]). Note that these transfer ionization channels contribute both to net capture and net ionization. This shortfall in ionized electron flux in the experimental data is difficult to understand in the context of modelling the projectile charge state dependence both in the IAM-PCM [8] and the time-dependent mean-field CTMC model [3], and thus we can only ask for additional experimental work in this context. Luna 16 FIG. 4. Net electron capture cross section for Li 3+ -water collisions. The curves follow the same pattern as described for Fig. 2. The experimental data are from Ref. [32].

B. Collisions with methane (CH 4 )
IAM-PCM and IAM-AR results for net capture of electrons from methane by proton impact were described previously in Ref. [7], where they were compared with other theoretical works, namely CDW-EIS and CNDO calculations. The CDW-EIS method was applied to higher projectile charges in Ref. [23].
In Fig. 5 we compare the net capture cross sections from IAM-AR and IAM-PCM for projectile charges Q = 1, 2, 3. For the methane target the difference between the two models at the lowest energy reported exceeds a factor of two. For alpha particle impact the experimental data fall below the present IAM-PCM results, and merge with the proton The curves merge at lower energies for proton impact (shown in blue), followed by alpha-particle impact (shown in green), and even higher energies for Li 3+ projectiles (shown in red). The experimental data are from Refs. [29,31,33,34].
impact data at 10 keV/amu collision energy.
For Li 3+ projectiles we find that the net capture cross section at low energies becomes 60 Å 2 , which is the same value that is reached by Li 3+ − C collisions (cf. Fig. 1). This can be called a saturation behavior in the sense that it is not the total number of valence electrons that determines the size of net cross sections for high projectile charges. We first investigate this phenomenon in detail for methane, and then make some comments further below as to which other target molecules will be affected in a similar way. Looking at the equivalent The saturation behavior in the net capture cross section as one goes to low impact energies and higher projectile charges can be illustrated as an overlap effect in the IAM-PCM.
A geometric condition for the saturation behavior in CH 4 is obtained from the following consideration: if r C and r H are the radii of the carbon and hydrogen cross-sectional disks with r X = σ X /π and b the bond length of C-H, then saturation happens if r C ≥ r H + b.
For molecules which involve large bond lengths saturation is much less likely to occur.
As shown in Fig. 1 net electron capture from carbon atoms for low energies exceeds net capture from atomic hydrogen by more than a factor of two for Q = 3 projectiles. The saturation condition r C ≥ r H + b is fulfilled up to E ≈ 80keV/amu. As a result of saturation the IAM-PCM cross section becomes independent of the orientation of the methane molecule.
This overlap effect is demonstrated for methane in Fig. 6 and analogous figures can be drawn for water vapor. To further illustrate the saturation phenomenon we show in Fig. 7  Clearly, saturation in the sense that the cross section of just one atom becomes an effective bound for the molecular cross section is not achievable for large molecules. If, however, 13 saturation does occur (as was observed in methane for carbon atoms vs hydrogen) then one can generalize, namely one finds σ CH 4 = σ C = σ CH 3 = σ CH 2 , i.e., the cross sections of the functional groups CH 3 and CH 2 are also limited by the cross section of carbon. For pure hydrocarbons (alkanes, alkenes, alkynes, aromatics, etc.) this means that they can be described as carbon clusters with molecular geometry. We begin with a comparison of the net capture cross sections for our three projectiles with charges Q = 1, 2, 3 in Fig. 8. The structure of the cross sections as a function of collision energy is similar to that obtained for methane (cf. Fig. 5 larger compared to that for methane by a factor of 1.5. This is in contrast with the IAM-AR results, where this factor is about four. We are not aware of experimental data for these cross sections with Q > 1, and therefore we cannot verify this aspect of the IAM-PCM prediction.
In Fig. 9 we compare our present results with the experimental data and with other theoretical data for proton impact. The experimental data of Tabet et al. The experimental data are from Ref. [35], the other theoretical data are from Refs. [36][37][38] Champion et al. [36]. The distorted-wave CNDO calculations of Purkaitet al. [37,38] cross the CDW results, while being lower below 100 keV and higher above the transition point in energy where the IAM-AR and IAM-PCM results merge. At low energies the IAM-PCM results are clearly the lowest by about a factor of three.

IV. CONCLUSIONS
As a follow-up to our previous work on net ionization of biologically relevant molecules by highly charged ions [8] we have presented results for net capture for three molecules. The choice of water and methane was motivated by the existence of some experimental data, while uracil is an example of a substantially larger molecule. Comparison with other theories shows that a number of them yield cross sections that are in the vicinity of Bragg additivity rule results (IAM-AR). For all molecules considered the IAM-PCM yields significantly reduced net capture cross sections at low energies with a particularly strong effect for higher projectile charges.
For the water molecule two classical trajectory based simulations were used, the standard CTMC model was generally found to be closer to the IAM-AR results, while the recently