1. Introduction
One of the most fundamental processes in nature is the scattering of charged particles off an atom or molecule. Atomic scattering occurs throughout the universe, from nuclear fusion, that powers stars, to hadron therapy for cancer, that uses ions to destroy deep-seated tumours [
1,
2]. For this reason, detailed studies of the underlying principles of ion–atom collisions remain an active area of research [
3,
4]. It is well known that, on the energy scale where these collisions occur, their dynamics is governed by quantum mechanics. A wave function satisfying Schrödinger’s equation of motion defines the development of the state of the particles through space and time. While the solutions to the Schrödinger equation of the hydrogen atom are analytically known, adding additional particles significantly complicates the problem so that exact solution is no longer possible. This makes scattering of another particle on a quantum-mechanical few-body system extremely complicated to model.
Significant differences still exist between theory and experiment for some of the most fundamental collision systems. In particular, calculation of the fully differential ionisation cross section remains a challenging problem. For example, the state-of-the-art quantum-mechanical convergent close-coupling approach to
collisions [
5] resulted in differences in the features of the fully differential emission when compared to the three-dimensional pattern experimentally observed by Schulz et al. [
6]. These calculations were performed in the high-energy regime where electron capture is believed to play a negligible role in the collision dynamics, significantly simplifying the theoretical problem. Currently available calculations of the doubly and triply differential cross sections for ionisation in
collisions also show inconsistent agreement with experimental measurements; see, for instance, refs. [
7,
8,
9]. Our ultimate goal was to address these discrepancies using the two-centre formulation of the wave-packet convergent close-coupling (WP-CCC) method that includes strong coupling effects between various reaction channels.
In this work we focus on the energy and angular distributions of electrons produced in proton collisions with molecular hydrogen. This system represents one of the most fundamental ion–molecule collision systems and describing this scattering process accurately is a key step towards developing theories for calculating ion scattering on more complex targets, such as water and other molecules of biological relevance. Accurate theories for calculating stopping power cross sections for such molecules are urgently required for Monte-Carlo simulation for treatment planning in hadron therapy [
1,
10]. For proton collisions with
, the ionisation cross section has a maximum around a projectile energy of 70 keV and becomes the dominant energy loss mechanism at higher incident energies. This makes calculations particularly difficult in this region since there are many competing processes that can occur as a result of scattering, and reaction channels are strongly coupled to each other. Additionally, the ionisation cross section differential, in both the energy and emission angle of the ejected electrons, provides a more thorough description of the outcome of scattering than the total or singly differential cross sections could do. For that reason, calculating the doubly differential cross section represents a more stringent test on theory.
Experimentally,
collisions are easier to investigate than the (theoretically simpler)
system due to the difficulties in preparing atomic hydrogen target atoms [
11]. As a result, the molecular target has been experimentally studied extensively [
12,
13]. The angular and energy distributions of ejected electrons were first measured simultaneously by Kuyattand Jorgensen [
14] across an incident energy range from 50 to 100 keV. Later, it was discovered that their low-energy data near the ionisation threshold was not accurate due to inefficiencies in the apparatus [
11]. Disagreement between their data and the first-order Born approximation (FBA) calculations available at the time prompted further investigations into this process. Ruddand Jorgensen [
15] measured the doubly differential cross section for ionisation at 100 keV over a wider range of ejected-electron energies and angles. They also concluded that although the energy distribution of the ejected electrons was represented reasonably by the FBA calculations, agreement between theory and experiment for the angular distribution was very poor. They extended their measurements to an incident energy of 300 keV [
16] where the measured total ionisation cross section demonstrated improved agreement with the FBA one. However, significant disagreement was still observed between theory and experiment for the doubly differential cross section for ionisation. Additionally, when the differential cross sections measured in these early experiments were integrated the resulting total ionisation cross section (TICS) was 35–
larger than reported by direct measurements. Toburen and Wilson [
17] measured the doubly differential cross section for ionisation in
collisions from 300 keV up to
MeV. They reported an uncertainty of
for their measurements of low-energy electrons and, in general, their results agreed well with those of Rudd et al. [
16] at 300 keV. However, their data deviated from that of Rudd et al. [
16] at ejection energies less than 50 eV for ejection angles greater than
. The most significant difference was seen at the highest measured ejection energies and for emission in the backward direction, where the cross sections were very small and the signal to noise ratio in their measurements was low. They found that the first Born approximation agreed well with their measured doubly differential cross section (DDCS) for ionisation as a function of ejection energy only if the long-range electron–proton interaction was included. At the larger impact energies considered by Toburen and Wilson [
17], the binary-encounter peak becomes visible for ejection of high-energy electrons at small emission angles. This secondary peak in the DDCS is the result of essentially binary collisions between the projectile and active electron. While binary-encounter theories based on this classical collision model have been developed, they cannot be applied to cross sections that are differential in the angle of emission. This is due to the fact that these cross sections are sensitive to the presence of the target nucleus, which is ignored in the binary-encounter theories. Thus, these types of approaches have been confined to calculation of the total cross section and the singly differential cross section as a function of the ejected electron energy. The DDCS for ionisation was also measured by Gibsonand Reid [
18] at 50 keV over a range of ejection energies from 5 to 150 eV and angles from 0 to
. However, the TICS found that integrating their results differed from direct measurements by
. More recent experiments called into question the accuracy of their results due to potential spreading of their proton beam [
19]. The greatest difficulty in obtaining accurate integrated cross sections for ionisation from measured DDCSs is the accuracy of the cross section data for very low energy electrons. This is where the DDCS is largest. Experiments can only measure to a finite minimum ejection energy, typically around 5 eV, and the value of the DDCS from this point to threshold must be estimated. The latest, and possibly, most accurate experimental measurements of this type were performed by Gealy et al. [
11] for incident protons with energies from 20 to 114 keV. They reported the DDCS for ionisation as a function of both ejection energy and angle from electron energies as small as
eV to 400 eV at the highest incident energy. Their angular range extended from 15 to
. This experiment used a new apparatus specifically designed to measure the flux of low-energy ejected electrons accurately. Furthermore, the integrated cross sections in this work agreed well with independently measured total ionisation cross sections for the
collision system.
Calculating the doubly differential cross section for ionisation as a function of ejected electron energy and emission angle in
collisions has proven to be a significant challenge for theorists. There are only a few attempts to calculate this cross section available in the literature. Moreover, the majority of calculations have used the FBA that is only applicable at high projectile energies. Additionally, scattering on the molecular hydrogen target is often treated with either the so-called Bragg additivity rule or an energy-scaling procedure. Both of these approaches begin by calculating cross sections for scattering on atomic hydrogen. Then, results for the
target are obtained by multiplying by a factor of 2, or applying the energy-scaling approach suggested in refs. [
20,
21]. The inability of these approaches to accurately model the underlying physics was demonstrated by the poor agreement with the experimental data of refs. [
14,
15,
16,
17]. A different approach was taken by Macek [
22], that used the first term in the Neuman expansion of Faddeev’s equations for the final state of the projectile-electron-residual ion system to calculate the DDCS for ionisation. The calculations were performed at an incident energy of 300 keV. The results showed improved agreement with the measured cross section compared to FBA calculations but still demonstrated significant deviation from the experimental data. In particular, the results of Macek [
22] significantly overestimate the cross section near small ejection energies. This is possibly due to the fact that the Faddeev equation approach is not valid for breakup in a system of three charged particles. In this case, the kernel of the Faddeev equations becomes non-compact. The doubly differential cross section for ionisation as a function of ejected electron energy and angle were calculated by Galassi et al. [
23], using both the FBA and the continuum distorted wave–eikonal initial state (CDW–EIS) method. As with the FBA, the CDW–EIS is a perturbative approach and only applicable at sufficiently high impact energies. Two methods for determining the cross section for the
target were used. The first one is based on the simple Bragg additivity rule, which multiplies the CDW–EIS results for
scattering by a factor of 2. The second one treats the target as two independent
atoms separated by the equilibrium inter-nuclear distance of the
molecule. The cross section for proton scattering on
is then given by the integral over all orientations of the product of the square of the magnitude of the amplitude for scattering on an effective atomic hydrogen target and an orientation-dependent factor resulting from the two scattering centres. The authors refer to the latter method as the continuum distorted wave–eikonal initial state–orbital (CDW–EIS–MO) method. Both approaches are compared to the experimental results of Gealy et al. [
11] at 114 keV as a function of ejected electron energy at four different ejection angles. Their results show that Bragg’s additivity rule is inadequate at providing accurate doubly differential cross section for ionisation, overestimating the results of ref. [
11], particularly for ejection at angles greater than
. Nevertheless, overall, the CDW–EIS–MO calculations demonstrate improved agreement with experiment.
Recently, we applied the two-centre wave-packet convergent close-coupling approach to calculate the energy and angular distribution of electrons emitted in collisions of 300-keV protons with
[
24], providing the most accurate theoretical description of this process to-date. The results demonstrated excellent agreement with the experimental data of Rudd et al. [
16] and Toburen and Wilson [
17]. Here, we provide details of the method and extend the calculations to more challenging lower incident energies, where the speed of the projectile is comparable to that of the target electrons. Specifically, we present the results at seven different projectile energies from 48 to 200 keV. In addition, we provide full details of the approach.
The layout of this paper is as follows. First, we provide a detailed review of the relevant literature. Then, we describe the presented method and give details of the calculations. Next, we present our results and provide a detailed discussion in the context of the experimental and theoretical data presently available in the literature. Unless specified otherwise, atomic unit (a.u.) is used throughout this manuscript.
3. Results
In this section, we present the results of our calculations for the doubly differential cross section for ionisation in collisions. Each figure below shows the DDCS at a fixed impact energy and the individual panels correspond to different emission angles. At the top of each panel we also show the speed of the ejected electron , in terms of the projectile speed. In the forward direction the ECC component peaks when . The smallest emission angle considered here is . Nevertheless, the presence of the ECC peak is still observed, though to a lesser extent, even for emission. We find that a basis containing bound states for each included orbital quantum number ℓ up to was sufficient to obtain converged results at all impact energies considered here. The continuum was discretised with 25 bins and the electron–momentum cutoff varied from a.u. to a.u., depending on the impact energy. The z-grid was extended from to a.u. The results obtained from Equation (29) are labelled as WP-CCC.
There are many combinations of ejection angle and ejection energy at different incident energies for which experimental data exist. Here, we present the results for a wide range of ejection energies and ejection angles for impact energies from 48 to 200 keV, and compared these with the available experimental data. Where possible, we also compared to previous calculations existing in the literature. In most of the panels we also show the DI and ECC components of the cross section. We omitted these lines in the panels that contain other calculations to prevent crowding of the figures.
Error bars for the experimental data are only given by Kuyattand Jorgensen [
14], and only for a few of their data points. However, this is the oldest data available and more recent works call into question the accuracy of their low-energy points [
11], suggesting a potentially larger uncertainty in their results than estimated in the original work. There are no error bars available for the other experimental data shown in the figures presented below, but a discussion of the estimated uncertainties is given in the original sources. Specifically, the data of Rudd et al. [
16] have an estimated uncertainty of
. The experimental data by Gealy et al. [
11] has an uncertainty of
at
eV, decreasing to
at ejection energies greater than 10 eV. However, at the highest ejection energies recorded the authors estimated at least a
uncertainty, due to the low count rates recorded in the experiment.
In
Figure 1 and
Figure 2 we present our results of the DDCS for ionisation as a function of the energy of the ejected electron at various emission angles for projectile energies of 48 and 67 keV. Here, the projectile speeds were
and
a.u., respectively. The only available data for the DDCS of ionisation at these energies were the measurements by Gealy et al. [
11]. To our best knowledge there are no other calculations at these collision energies. For ejection angles of
and smaller we observed excellent agreement between our results and the experimental data. At 130 and
our DDCS agreed well with experiments for low-energy electrons but for high-energy electrons the WP-CCC approach overestimated the experimentally measured cross section at backward emission angles. There could be two possible explanations for this. First, this is the region in which the DDCS is several orders of magnitude smaller than for low-energy electrons and, therefore, the experimental measurements have the greatest uncertainty due to low counting rates. Second, ionisation resulting in emission in the backward direction is dominated by the DI mechanism. The DI component of the ionisation cross section is given by the product of the amplitude for excitation into a positive-energy pseudo-state and the overlap between the true continuum wave and corresponding wave-packet pseudo-state. The pseudo-states and continuum wave are both found by numerically solving the Schrödinger equation of the model target potential. At asymptotically large distances this potential tends towards a Coulomb attraction with charge
a.u. However, in reality the electron would be on the opposite side of the residual target ion from the projectile proton, resulting in an effective Coulomb interaction at asymptotically large distances with charge
a.u. Hence, in our approach to calculating DI the emitted electron does not account for the additional attraction, due to the projectile proton. Therefore, the current approach may somewhat overestimate the probability of emission in the backward direction.
Figure 3 shows the DDCS for ionisation at 75 keV for ejection angles from 10 to
. For the smallest emission angle we found that the ECC component of the DDCS peaked when the speed of the ejected electron was slightly less than that of the projectile, meaning that the likelihood of the electron staying slightly behind the outgoing projectile was higher. The projectile speed was
a.u. at this impact energy. The ECC peak was strongest for emission at
. Our calculations suggested it had an influence on the cross section for emission, even at
. At
, the DI and ECC components were comparable for most of the considered electron energies. Then, for emission into angles of
and greater, the DI became the dominant mechanism contributing to the total DDCS for ionisation. Only for very high electron energies did the ECC component have a significant effect on the overall results at 120, 140 and
. For emission into the backwards direction the DDCS was entirely dominated by DI, as expected.
In
Figure 4 we present our results for the DDCS for ionisation at 95 keV for ejection angles from 20 to
. The projectile speed at this energy was
a.u. Our results were compared with the experimental data of Gealy et al. [
11]. To our best knowledge there are no other experimental measurements or calculations at this collision energy. We observed very good agreement between the WP-CCC calculations and the experimental data. At the smallest considered ejection angle of
, the interplay between the DI and ECC contributions near the matching speed resulted in an overall smoother cross section that followed the experimental data more closely than observed at lower projectile energies. At this incident energy we also saw improved agreement between our results and the experimental data for high-energy electrons emitted at 130 and
degrees, compared to 48 and 67 keV.
Figure 5 shows the DDCS for ionisation at 100 keV. Here, we compared our results with the experimental data of Ruddand Jorgensen [
15] and Kuyattand Jorgensen [
14]. Ruddand Jorgensen [
15] also performed the FBA calculations. The FBA results significantly underestimated the cross section at an emission angle of
, whereas the WP-CCC calculations showed excellent agreement with the experimental data across the entire measured energy range from threshold to 300 eV. We also saw a clear shoulder, due to the ECC peak for emission at
, in both the experimental data and our calculations. At this impact energy, the matching speed was
a.u. At
the WP-CCC results agreed well with the experimental data of Kuyattand Jorgensen [
14], except for small energies of the emitted electrons. However, these experimental data were unreliable in that region, due to inefficiencies in the apparatus at detecting low-energy electrons [
11]. In this energy region the DDCS was dominated by the DI mechanism. At energies above 30 eV, ECC was the most significant contributor to the ionisation cross section. At
, the FBA calculations continued to underestimate the experimental data, although less significantly, whereas the WP-CCC calculations followed the experiment very closely at all ejection energies. For emission into an angle of
, we found good agreement with the experimental results of Kuyattand Jorgensen [
14], except for the narrow region near the ejection energy of 30 eV, where the DI and ECC components contributed equally to the ionisation cross section. Here, our calculation underestimated the experiment. At an emission angle of
, we also observed a similar difference between our result and the experimental data of Ruddand Jorgensen [
15]. However, the difference was smaller than at
. For ionisation resulting in electron emission at
, the WP-CCC results agreed well with the experimental data, although as the ejection energy increased our results began to underestimate the measured cross section. The FBA calculations by Ruddand Jorgensen [
15] showed similar behaviour, except beyond 100 eV, they fell off more steeply, deviating from the experiment and the WP-CCC calculations. The FBA results showed similar behaviour also at
, whereas the WP-CCC results displayed excellent agreement with the experiment emission energies. Lastly, at
emission the FBA overestimated the experimental cross section near the ionisation threshold and then fell off too steeply with increasing energy. Our calculations agreed well with the data of Ruddand Jorgensen [
15] for low-energy electron emission, slightly underestimating the experimental data above 50 eV.
Figure 6 shows the DDCS at an incident energy of 114 keV, along with the experimental data of Gealy et al. [
11] and the FBA and CDW-EIS-MO calculations by Galassi et al. [
23]. The projectile speed was
a.u. At an emission angle of
, the WP-CCC calculations agreed very well with the experimental data below 250 eV emission energy. Above 250 eV, however, our results slightly overestimated the experiment. The same was observed at 30 and
. The FBA calculations of Galassi et al. [
23] suggested a shoulder at 150 eV for
emission, which was not present in the other calculations or experimental data. At 30 and
, the CDW-EIS-MO results agreed well with the experiment. At
and above, the WP-CCC results agreed well with the experimental data for smaller ejection energies, then tended to overestimate the data as the emission angle and energy increased. At emission angles of 90 and
, the FBA results overestimated the cross section for low-energy electron emission, but agreed well with the experimental data at high ejection energies. The CDW–EIS–MO calculations showed the opposite behaviour, agreeing well with the experiment at small ejection energies, but underestimating it at larger emission energies. We concluded that the WP–CCC method reproduced the experimental cross section more consistently than the FBA and CDW–EIS–MO methods.
In
Figure 7 we present our results for a collision energy of 200 keV in comparison with the experimental data of Rudd et al. [
16]. For emission into
, the experiment showed a clear shoulder near the matching speed of
a.u., which was perfectly reproduced by the WP-CCC results. We found good agreement between the WP-CCC results and experimental data across all ejection angles and energies for which experimental data was available. At the emission angles of 10, 30 and
, our calculations slightly deviated from the experimental data when moderate-energy electrons were emitted. At this impact energy, we saw the binary-encounter peak clearly resolved at small emission angles. In particular, at
, our calculations showed a wide, but pronounced, peak at 350 eV, reproducing the shape of the measured cross section. There was no experiment at
. Here, the WP-CCC results are given for completeness.
At all impact energies considered we also calculated the DDCS via numerical integration over the fully differential cross section; see Equations (
23) and (
24). In all cases, the results were very similar to those found from analytical integration. At all impact energies considered, there was a slight difference at large ejection energies, where the “analytic” WP-CCC result was, generally, marginally larger than the “numerical” WP-CCC calculations, especially for ejection into angles greater than
. However, this happened only at the largest emission energies, where the corresponding DDCS was extremely small. Thus, overall excellent agreement between the “numerical” and “analytical” WP-CCC results provided a reassuring self-consistency check for our calculations. Additionally, the accuracy of the numerical integration was checked by increasing the density of integration points until convergence was obtained in the result. Thus, we could conclude that the small differences observed between the numerical and analytical WP-CCC calculations were the result of the interpolation of the ECC amplitudes when transforming them into the common laboratory frame. To obtain the numerical WP-CCC results, we interpolated the expansion coefficients
, as discussed in
Section 2.3. However, to obtain the ECC component of the cross section for the analytical WP-CCC results, we interpolated the integral appearing in the second term of Equation (29). This quantity, known as the density matrix, is also smooth but has a different structure. Hence, we concluded that the difference between our numerical and analytical WP-CCC results provided an estimate of the uncertainty associated with the interpolation procedure used to bring the ECC amplitudes into a common frame with the DI part. The difference was insignificant.
The results given in
Figure 1,
Figure 2,
Figure 3,
Figure 4,
Figure 5,
Figure 6 and
Figure 7 represent the doubly differential ionisation cross section as a function of two variables. Integrating them over either the emission angle or energy of the ejected electron resulted in the singly differential cross section for ionisation as a function of the remaining variable. We performed this integration numerically and compared the results to our earlier calculations of the SDCS for ionisation presented in ref. [
30]. Both types of the SDCS obtained from the DDCS presented here agreed perfectly with the previously published results. Furthermore, we integrated each of these two SDCS and found that the result was in very good agreement with the TICS calculated directly from the expansion coefficients [
32].
Let us now outline the overall picture. In general, we found excellent agreement between the WP-CCC results and the experimental data of Rudd et al. [
16], but our calculated cross section was slightly smaller than the experimental points from Rudd [
40]. Agreement between our calculations and the data of Gealy et al. [
11] was best at small ejection angles and larger impact energies. We observed the largest difference at high ejection energies for ejection in the backward direction. We saw a small dip in the calculated cross section at low ejection energies, where the DI and ECC components contributed equally, that was not reflected in the experimental data. We compared our results to the FBA and CDW–EIS–MO calculations of Galassi et al. [
23] at ejection angles of 30, 50, 90, and
, available for an impact energy of 114 keV. Agreement between the FBA results and experimental data of Gealy et al. [
11] varied considerably between different ejection angles. At
, the FBA calculations underestimated the experimental data around 60 eV, then overestimated it starting from 150 eV. At larger ejection angles the overall trend agreed better with the experimental points, except at small ejection energies, where the FBA results significantly overestimated the cross section at backward ejection angles. The CDW–EIS–MO calculations are the most sophisticated calculations currently available in the literature. They agreed slightly better with the experimental data at
than the WP-CCC results, but, at larger ejection angles, they began to underestimate the experimental cross section as the ejection energy increased. Thus, we concluded that, overall, the present WP-CCC method could accurately model the underlying physics and provide reliable DDCS for ionisation in the
collisions.
Across all impact energies considered we found that, for small ejection angles, DI and ECC dominated the DDCS for small and large emission energies, respectively. The DI and ECC components were comparable to one another at an emission angle of about for all projectile energies. For ejection at and above, DI contributed the most to the DDCS. ECC was significant only for very large electron energies when the emission angle was greater than . For emission into angles greater than , our results overestimated the experimental data for high-energy electrons. Further investigation, with an improved target description, is required to understand this discrepancy between theory and experiment.
4. Conclusions
We applied the wave-packet convergent close-coupling approach to proton collisions with molecular hydrogen to calculate the ionisation cross section differential in both the angle and energy of the ejected electrons. A wide range of experimental results exist, spanning from the initial measurements of Kuyattand Jorgensen [
14] to the latest data of Gealy et al. [
11]. However, theorists have struggled to accurately reproduce these experiments. All currently available theoretical results are obtained using perturbative approaches that are applicable only at sufficiently high projectile energies. The fundamental assumptions that these methods are based on become questionable at the incident energies considered in this work, where the orbital speed of the target electrons is comparable to, or somewhat greater than, the projectile speed. The present WP-CCC method is the first non-perturbative approach to differential ionisation that is applicable at all considered incident energies. Furthermore, it is the only method applied to this collision problem that accounts for the strong coupling between reaction channels, including capture into the continuum. Our results demonstrate excellent agreement with all available experimental data at projectile energies from 48 to 200 keV and electron emission angles less than
. We observed a consistent overestimation of the experimental data for emission into backward angles greater than
. It is suggested that this may be due to the large uncertainty in the experimental data in this region. In part, this could also be due to the approach to calculating the total breakup amplitude used herein that assumes a weaker effective charge experienced by the emitted electron at asymptotically large distances from the origin than it actually might be in this particular kinematic regime.
Our results agree particularly well with the most recent experiment by Gealy et al. [
11] and the high projectile energy data by Rudd et al. [
16]. At the highest impact energies considered the WP-CCC method reproduces both the main peak in the DDCS near the ionisation threshold and the binary-encounter peak. Furthermore, we explicitly calculated direct ionisation and electron capture into the continuum contributions to the DDCS. This allowed us to analyse the roles these mechanisms play in emission of electrons in a wide range of kinematic regimes.
This work paves the way for applying the WP-CCC approach to calculations of the other types of DDCS for ionisation, such as that where the differential is in both the scattering angle of the projectile and the ejected electron energy [
41,
42]. Next, we plan to extend our calculations to ionisation resulting in emission in the forward direction, where ECC dominates the DDCS. Calculating the DDCS for ionisation in the forward direction has proven a challenging problem. Theoretical investigations are limited to the narrow region where the emitted electron velocity matches the projectile one. So far, there is no theory that properly describes the entire energy range of electrons ejected in the forward direction, including the low-energy and velocity-matching regions. The WP-CCC approach may be able to provide new insight into this problem.
We also plan to apply the two-centre WP-CCC method to calculate the doubly and fully differential cross sections for ionisation of both
and
by proton impact. We recently calculated doubly differential cross sections for ionisation in proton–helium collisions and obtained excellent agreement with the available experimental data [
43,
44]. Currently, there is no theory that can describe the experimental data for the fully differential cross section at 75 keV, recently measured consistently by Schulz et al. [
45] and Dhital et al. [
9] in all kinematic regimes. Existing theoretical methods deviate from one another significantly and cannot reproduce the experiment well. The single-centre version of the WP-CCC has already been applied to the fully differential cross sections for ionisation in proton–helium collisions [
46,
47] at sufficiently high impact energies, where the electron-capture channels are believed to be negligible. The results agreed very well with the recent high-resolution experiment by Gassert et al. [
48]. However, the approach predicted a slight, but unexpected, shift in the position of the binary peak compared to the experimental data. Furthermore, recent measurements of the fully differential cross section for ionisation in
collisions in various geometries also demonstrate significant differences with currently available distorted-wave theories [
9]. We intend to apply the two-centre WP-CCC method to calculate the fully differential cross sections in proton collisions with helium and molecular hydrogen to determine if higher-order effects, resulting from the presence of the second centre, are able to explain the aforementioned discrepancies.