Next Article in Journal
Quadrupole Effects in the Photoionisation of Sodium 3s in the Vicinity of the Dipole Cooper Minimum
Previous Article in Journal
The Double-Nozzle Technique Equipped with RF-Only Funnel and RF-Buncher for the Ion Beam Extraction into Vacuum
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Single Ionization of Helium by Protons of Various Energies in the Parabolic Quasi-Sturmians Approach

by
Alexander S. Zaytsev
1,
Darya S. Zaytseva
1,
Sergey A. Zaytsev
1,*,
Lorenzo U. Ancarani
2,
Ochbadrakh Chuluunbaatar
3,4,
Konstantin A. Kouzakov
5 and
Yury V. Popov
6,7
1
School of Fundamental and Computer Sciences, Pacific National University, 680035 Khabarovsk, Russia
2
Département de Physique, Université de Lorraine, CNRS, LPCT, 57000 Metz, France
3
Meshcheryakov Laboratory of Information Technologies, Joint Institute for Nuclear Research, 141980 Dubna, Russia
4
Institute of Mathematics and Digital Technology, Mongolian Academy of Sciences, Ulaanbaatar 13330, Mongolia
5
Department of Nuclear Physics and Quantum Theory of Collisions, Faculty of Physics, Lomonosov Moscow State University, 119991 Moscow, Russia
6
Skobeltsyn Institute of Nuclear Physics, Lomonosov Moscow State University, 119991 Moscow, Russia
7
Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 141980 Dubna, Russia
*
Author to whom correspondence should be addressed.
Atoms 2023, 11(10), 124; https://doi.org/10.3390/atoms11100124
Submission received: 30 August 2023 / Revised: 16 September 2023 / Accepted: 25 September 2023 / Published: 28 September 2023
(This article belongs to the Section Atomic, Molecular and Nuclear Spectroscopy and Collisions)

Abstract

:
The parabolic quasi-Sturmian approach, recently introduced for the calculation of ion–atom ionizing collisions, is adapted and applied here to the single ionization of helium induced by an intermediate-energy proton impact. Within the method, the ionization amplitude is represented as the sum of the products of the basis amplitudes associated with the asymptotic behavior of the continuum states of the two noninteracting hydrogenic subsystems ( e , He + ) and ( p + , He + ) . The p e interaction is treated as a perturbation in the Lippmann–Schwinger-type (LS) equation for the three-body system ( e , He + , p + ) . This LS equation is solved numerically using separable expansions for the proton–electron potential. We examine the convergence behavior of the transition amplitude expansion as the number of terms in the representation of the p e interaction is increased and find that, for some kinematic regimes, the convergence is poor. This difficulty, which is absent for a higher proton energy impact, is solved by varying the momentum of the auxiliary proton plane wave introduced into the basis function. Fully differential cross-sections are calculated and compared with experimental data for 75 keV protons and the results obtained with the 3C model.

1. Introduction

Ionization in the collisions of ions with simple atoms is one of the fundamental processes of atomic dynamics and has been the focus of study for several decades. However, direct measurements of fully differential cross-sections (FDCS), which contain complete collision information, became available only after the development of cold target recoil ion momentum spectroscopy (COLTRIMS) [1,2,3]. For the single ionization of the ground state of helium by high-impact proton energy (1 MeV), FDCS have been recently measured [4] and triggered a number of theoretical calculations [4,5,6,7,8,9,10,11,12]. Despite some unexplained shifts in the angular distributions, overall satisfactory agreement was found. Moreover, those experimental data being on a relative scale, the comparison of absolute theoretical FDCS does not permit us to reach a conclusive statement on the merit of each approach. At a lower impact energy, the absolute FDCS of singly ionizing 75-keV p H e collisions have been measured in various kinematical regimes [13]. For these FDCS, theoretical predictions were obtained based upon the well-known continuum distorted-wave-eikonal initial-state (CDW-EIS) model [14,15], (a) including the projectile and residual target ion interaction (CDW-EIS PI) [16]; (b) the two-Coulomb final wave (2C) + the initial Coulomb projectile wave [17]; (c) the three-Coulomb wave (3C) model [18]; (d) the modified Coulomb–Born approximation (MCB-PI) [19]; (e) the continuum correlated wave model including the interaction between the projectile and the residual target ion (CCW-PT) [20]; (f) including the projectile target core (NN) interaction model convolved with experimental resolution (CDW-EIS NN) [21]. Despite the sophistication of the latest approaches, none of them has managed to reach complete agreement with the experimental data. The most noticeable theory–experiment discrepancies are observed at large values of the transferred momentum. In particular, the proposed theoretical models fail to adequately describe the observed double-peak structure of the FDCS in the perpendicular plane. In the lower-impact energy regime, the electron capture channel may be expected to play an important role and may partially explain such discrepancies. Indeed, the FDCS could be affected, for example, by the well-known Thomas effect [22] that appears when the electron after the second collision with the target nucleus moves in parallel with the same speed as the incident ion. In such a kinematic domain, the particles interact strongly, leading to a resonant increase in the FDCS. Note that the same effect is postulated in a recent work on the Compton disintegration of positronium [23] and is actively discussed in a publication [24] directly related to our study. As far as we know, and although it should be properly taken into account, so far, none of the theoretical approaches modeling the 75 keV measurements [13] have included this channel. Thus, further theoretical developments to describe ion–atom ionizing collisions are required.
Recently [5], some of us proposed the parabolic convoluted quasi-Sturmian (CQS) approach to the FDCS calculation for the ionization process under consideration. The method was successfully tested at high projectile energy, namely with protons at 1 MeV. The results were in good agreement with other theories, and especially with the absolute scale of the WP-CCC approach [25]. In this work, we wish to explore other kinematical situations. We extend the absolute scale comparison with the theoretical predictions of the WP-CCC approach [25] for protons incident at 0.5 and 2 MeV and momentum transfers varying from 0.5 to 1.75 a.u. More importantly, the main aim of this manuscript is to test the applicability and robustness of the CQS approach to a lower incident energy regime (75 keV) for which the absolute measurements [13] offer a serious challenge. In its present version, the CQS approach does not yet include the capture channel, so its application in the lower-energy regime will allow us to compare the cross-sections with those of other theoretical approaches having the same weakness. Clearly, if our CQS approach yields an overall similar picture and if unexplained features (in shape and/or absolute scale) subsist, the inclusion of the capture channel will be identified as the next necessary ingredient to be included in our approach.
Briefly, in the CQS approach [5], one treats the problem in parabolic coordinates with the z ^ axis chosen along the incident proton momentum K 0 . The ionization amplitude is represented as an expansion in the so-called basis amplitudes associated with the asymptotic behavior of the Green’s functions of the two subsystems ( e , He + ) and ( p + , He + ) . In both cases, the object of action of the Green’s function operator is the orthogonal complements to the square integrable Sturmian basic functions of the parabolic coordinates. The transition amplitude expansion coefficients are obtained as a solution to the Lippmann–Schwinger (LS) type equation for the final-state three-body system ( e , He + , p + ) , in which the proton–electron interaction is considered as a perturbation. In order to properly take into account the p e potential in the intermediate-energy regime under consideration, we introduce into the Sturmian basis functions an auxiliary proton plane wave with a momentum Q = ε K 0 , with ε 1 . Contrary to the original implementation proposed in [5], the quantity Q is treated here as a variational parameter that is chosen according to the kinematical values under scrutiny. This modification in the CQS approach is necessary to obtain converged FDCS that can be compared to the 75 keV measurements.
The paper is organized as follows. In Section 2, taking into account the asymptotic behavior of the Green’s function of the Coulomb three-body system ( e , He + , p + ) and making the approximation of the frozen-core model for the residual helium ion, we present the main equation of the method. Section 2.1 proposes a representation of its solution in the form of an expansion in terms of convolutions of quasi-Sturmian functions corresponding to noninteracting subsystems ( e , He + ) and ( p + , He + ) . From this expansion, which presumably converges in the high-energy regime, we deduce the corresponding expansion of the ionization amplitude in terms of the basis amplitudes. In Section 2.2, the matrix representation of the equation for the expansion coefficients of the ionization amplitude is discussed. In particular, the role of the additional pole, which owes its origin to the modification of the exponents of the Sturmian basis functions by introducing the auxiliary proton plane wave Q in the calculation of the matrix element of the Green’s function, is clarified. We show that by varying Q, we can influence the accuracy of the description of the proton–electron interaction within the framework of this approach. In Section 3, the results of our numerical calculations are presented. First, we examine the convergence issues of the differential ionization cross-sections, particularly with respect to the number of terms in the representation of the proton–electron interaction and the stabilizing role played by the variational parameter Q. Then, we make a comparison with the 75 keV experimental data and theoretical cross-sections obtained by other authors. We also present CQS predictions for FDCS for the cases of 0.5 and 2 MeV protons and compare them with the results of the WP-CCC calculations [25]. Finally, Section 4 provides a summary of this work.
Atomic units (a.u.) in which = e = m e = 1 are used throughout, unless otherwise specified.

2. Theory

We wish to study the ionization process
He ( 1 s 2 ) + p + e + p + + He + ( 1 s )
in which a proton with momentum K 0 ionizes the target, ejecting an electron with momentum k e and energy k e 2 / 2 . The helium nucleus will be considered at rest, and R will denote the relative coordinate of the proton, which, in the final state, has momentum K . In our model, we take a good ground state representation Φ ( 0 ) ( r , r ) of the helium target and make the frozen-core approximation for the residual helium ion. In the parabolic CQS approach, the calculation of the amplitude of the ionization process reduces to solving the inhomogeneous equation
E H ^ Φ ( + ) = K 0 , F .
The right-hand side of (2) is the product of the plane wave
R K 0 = e i K 0 · R
for the projectile and the matrix element
F ( R , r ) = ψ 1 s H e + V ^ i Φ ( 0 )
of the incident channel interaction
V i = 2 R 1 | R r | 1 | R r |
between the frozen electron wave function ψ 1 s H e + ( r ) and a helium wave function Φ ( 0 ) ( r , r ) . The final three-body channel Hamiltonian
H ^ = H ^ 0 + U ^
is split into a separable part
H ^ 0 = 1 2 m p R 2 1 2 r 2 + 1 R 1 r
and the proton–electron interaction
U ^ = 1 R r
considered as a perturbation.
The ionization amplitude T K , k e is contained in the leading asymptotic form (for large values of the hyperradius ρ = m p R 2 + r 2 ) of the solution to the driven Equation (2):
Φ ( + ) R , r m p ( 2 π ) 2 ( 2 E ) 3 / 4 e i π 4 ( 2 π ) 1 / 2 exp i 2 E ρ + W 0 ( R , r ) ρ 5 / 2 1 2 T K , k e ,
where W 0 is the Coulomb phase defined by [26]
W 0 ( R , r ) = ρ 2 E 1 R 1 r 1 R r ln 2 2 E ρ .

2.1. Parabolic Sturmians

We consider the representation of Equation (2) in the square-integrable Sturmians defined by [27]
ξ , η , ϕ | n , m , ϰ = e i ϰ ϕ 2 π φ n | ϰ | ( ξ ) φ m | ϰ | ( η ) ,
where
φ n λ ( ρ ) = 2 b n ! ( n + λ ) ! ( 2 b ρ ) λ / 2 e b ρ L n λ ( 2 b ρ ) , λ = | ϰ | ,
with the basis scale parameter b. The basis vectors are represented by products
Q , N Q , n m
of an auxiliary proton plane wave Q and two Sturmians
n n 1 , n 2 , ϰ , m m 1 , m 2 , ϰ
of the parabolic coordinates ξ 1 , η 1 , ϕ 1 and ξ 2 , η 2 , ϕ 2 corresponding to the proton R and electron r position vectors, respectively. Note that the introduced projectile plane wave (it may be compared with the impact parameter model [28], where the proton part of the wave function is approximated in a similar way) is a key ingredient of the approach, as it allows us—albeit partially—to compensate for the oscillating term on the right-hand side of Equation (2). In addition, as will be shown below, we exploit the freedom to choose the value of Q in such a way as to optimize our numerical treatment of the p e interaction.
Further, we propose to look for the solution Φ ( + ) in the form of an expansion
Φ ( + ) = N C N S Q , N ( + )
in terms of CQS functions whose number is determined by some limits M and N to the ranges | ϰ | M and n j , m j < N , j = 1 , 2 . The CQS functions are constructed in such a way as to convey the asymptotic behavior (9) to the solution as completely as possible, namely
S Q , N ( + ) = G ^ 0 ( + ) ( E ) Q , N ˜ ,
where G ^ 0 ( + ) E = E H ^ 0 1 and Q , N ˜ = Q , n ˜ m ˜ is the orthogonal complement to the basis vectors Q , N .
Here, the orthogonal complement l ˜ to each Sturmian
l = n , m , ϰ ,
is defined by
ξ , η , ϕ l ˜ = 4 ξ + η ξ , η , ϕ n , m , ϰ .
By comparing the CQS (16) asymptotic behavior with (9), it has been found that the amplitude T K , k e is expressed (up to a phase factor) in terms of so-called basis amplitudes A n ( + ) and A m ( + ) [5,29]:
T K , k e = 2 2 π N C N e i ϰ ( ϕ p ϕ e ) A n ( + ) ( β p , K , Q ; θ p ) A m ( + ) ( β e , k e , 0 ; θ e ) ,
where β p = m p K and β e = 1 k e are the Sommerfeld parameters for the two subsystems ( p + , He + ) and ( e , He + ) . Explicit analytical expressions for these basis amplitudes are provided in [5,29]. The main numerical task consists then of calculating the coefficients C N .

2.2. Matrix Equation for the Coefficients C N

Inserting the expansion (15) into (2), and projecting by Q , N , gives the following matrix equation:
N δ N , N M N U ^ M Q , M ˜ G ^ 0 ( + ) ( E ) Q , N ˜ C N = N K 0 Q , F .
In our approach, the p e potential (8) is treated as a perturbation (in the high-energy regime) and is approximated by a truncated Sturmian basis set (13) expansion
U ^ N 0 = N , M N ˜ N U ^ M M ˜ ,
whose size N 0 = ( 2 M 0 + 1 ) N 0 4 is determined by the limits M 0 and N 0 to the ranges | ϰ | M 0 and n j , m j < N 0 ( j = 1 , 2 ) .
It is intuitively clear that, in the general case, the LS-type equation is inapplicable when taking into account the interaction U ^ , which is negligible compared to the energy E. On the other hand, the presence of the proton plane wave Q in the basis function (13) leads to the appearance of the term Q 2 2 m p in the matrix element of the operator E H ^ 0 . Roughly speaking, the use of these modified basis functions transforms the energy E E e f f = E Q 2 2 m p ; as a consequence, the values of the matrix elements of the Green’s function operator are increased. Thus, one might expect that there is an optimal value of Q for which the p e potential can be properly described when using the LS equation.
The matrix elements of the Green’s function are evaluated (numerically) employing the convolution integral representation [30,31]
Q , N ˜ G ^ 0 ( + ) ( E ) Q , N ˜ = 1 2 π i C d E m ˜ G ^ C ( + ) ( 1 , 1 ; E ) m ˜ × Q , n ˜ G ^ C ( + ) ( 1 , m p ; E E ) Q , n ˜ ,
where the path of integration C is obtained by a rotation of the contour C 0 through a negative angle φ about some point E 0 on the positive real axis (see Figure 1). In (22), G ^ C ( + ) denotes the Green’s function G ^ C ( + ) Z , M ; E = E H ^ C 1 for a Coulomb system with the Hamiltonian
H ^ C = 1 2 M X 2 + Z X ,
where X is the position vector.
The integration contours C and C in the electron and proton complex energy planes E e and E p , respectively, are presented in Figure 1. Note that the use of the modified exponents b i Q 2 and b + i Q 2 in the proton basis functions φ n 1 λ ( ξ 1 ) and φ n 2 λ ( η 1 ) , respectively, leads to the appearance of pole singularities in the matrix elements of the proton Green’s function at momenta P ± = ± Q 2 i b . The question arises as to how one should deal with the contribution to the contour integral (22) from the superfluous ‘resonance’ pole at
E = P 2 2 m p = Q 2 4 b 2 2 m p i 2 Q b m p ,
which lies on the unphysical sheet: 2 π arg ( E p ) < 0 . The criterion here can be the condition that the matrix of the Green’s function operator G ^ 0 ( + ) must be inverse to the matrix representation of the operator E H ^ 0 in the basis set (13). Our numerical calculations show that this condition is satisfied if the residue of this pole is included in the integral. The simplest way to ensure this is to choose a sufficiently small absolute value of the contour rotation angle:
| φ | < tan 1 4 Q b Q 2 4 b 2 2 m p ( E E 0 )
(see Figure 1). We found that in the energy regime considered in the numerical applications, for Q K 0 , the residue of the pole constitutes the lion’s share of the integral (22).
From the above, it follows that the optimal choice for Q, which provides the most accurate account of the proton–electron interaction achievable within our approach, is the admissible value closest to K 0 . Above such upper bound Q 0 , the p e potential can no longer be considered as a perturbation (compared with E e f f ), thus resulting in the very poor convergence of the FDCS as more terms are included in the Sturmian representation (21).

2.3. Numerical Scheme

In summary, the CQS approach to calculating FDCS follows schematically the following steps. First, we set the limits M 0 and N 0 for the expansion (21) used to represent the perturbation U ^ . We also calculate, for a given energy E, the Green’s function matrix elements (22) by adequately choosing Q. With these elements, we solve the matrix Equation (20) and obtain the coefficients C N . Next, we set the limits M and N for the expansion (15) in CQS and, for given kinematical and geometrical configurations, we sum the basis amplitudes according to (19) and obtain the ionization amplitude T K , k e . Finally, we calculate the FDCS, which, in the laboratory frame, reads
d 5 σ d E e d Ω e d Ω p = k e m p 2 ( 2 π ) 5 K K 0 T K , k e 2 .

3. Results

All FDCS presented hereafter are on absolute scale.
We start by discussing our CQS results for the FDCS for a proton incident energy of 75 keV and ejected electrons of 5.4 eV, corresponding to the experimental situation [13]. The experimental FDCS values were measured in absolute scale for electron ejection both into the scattering and into the perpendicular plane at different values of the transverse momentum transfer [13], noted here η as in the paper [21].
The CQS calculations have been carried out using a quite accurate ground-state wave function Φ ( 0 ) ( r , r ) obtained by diagonalizing the helium Hamiltonian matrix representation in the complete square-integrable Laguerre basis set [32]. In this case, the ground-state energy is E 0 = 2.9033 a.u. For the value of the basis (13) scale parameter b = 1 , the limits M = 3 and N = 21 to the ranges | ϰ | M and n j , m j < N , j = 1 , 2 are sufficient to reach convergence for the 2C-like amplitude, i.e., the ionization amplitude (19) calculated with the p e potential switched off. In doing this, we put Q = 0 . Then, we switch on the p e interaction and set the value Q = K 0 : it was found that for the transverse transferred momenta η = 0.13 , 0.41 , 0.73 , and 1.38 , the FDCS diverges with the increasing size of the matrix representation (21) of the proton–electron interaction. We thus searched for a more adequate Q value, for each transferred momentum. Figure 2 presents the FDCS for M 0 = 3 and increasing values of N 0 from 5 to 9. It can be seen that the FDCS convergence for η = 0.13 and η = 0.41 is achieved at K 0 Q = 4 , while, for η = 0.73 , it is required to increase K 0 Q to 5. Finally, for the largest η = 1.38 , a further increase in K 0 Q is necessary. In particular, Figure 2 shows the convergence behavior for η = 1.38 at K 0 Q = 7 .
Figure 3 and Figure 4 show our converged CQS results for the ionization FDCS for the collision and perpendicular planes in comparison with the experimental data [13] and the 3C calculations [33]. These 3C results were obtained with a strongly correlated function (CF) [33]; however, we recall that it was found that the 3C predictions for the FDCS depend weakly on the wave function of the helium ground state. The comparison between our CQS and the 3C model calculations is made here in order to focus on the effect of the phase factor (corresponding to the Coulomb proton–electron interaction), which is explicitly present in the 3C wave function.
One can observe that, for all the cases, both the experimental FDCS and the theoretical ones (see also [16,18,19,20,21]) in the scattering plane exhibit a strong peak near the direction of the momentum transfer (indicated in each panel by the arrow). At the same time, the theoretical prediction for the binary peak position is shifted towards smaller angles relative to the experimentally observed value for all momentum transfers, except for η = 1.38 .
In the perpendicular plane, the double-peak structure of the measured FDCS at η = 1.38 is reproduced in the 3C model calculations only qualitatively (see also [18]). Similar results were obtained in the MDC-PI [19] and CCW-PT [20] calculations. On the other hand, our CQS calculations (see also Figure 4) do not predict this FDCS feature.
It should be noted that for a proper comparison between the theory and experiments, and between theories, in Figure 2, Figure 3 and Figure 4, all theoretical calculations and the measurements are presented on an absolute scale. This is not the case in some other papers in which scaling factors were purposely used to focus on shapes’ comparison. One should also keep in mind that the calculations do not account for the uncertainties due to the finite energy and angular resolutions of the experiment. Overall, the magnitude of the binary peak in the eight panels is not unreasonable. Another observation is the fact that the recoil to binary peak ratios for the CQS approach are always larger than with the 3C model, except for η = 1.38 and in the scattering plane (see also [18,21]). The absence of experimental points in this angular region, however, does not permit us to favor one approach with respect to the other.
To conclude, we briefly return to the choice of the proton Q value and the convergence problem discussed above. We suspect that the difficulty arises at relatively large proton scattering angles θ p and not only at large transferred momenta. For example, the value η = 1.38 in the case considered in the experiments corresponds to θ p = 0.43 mrad. If we consider the kinematic ionization regimes with impact energies of 0.5 and 2 MeV (as in [25]), then θ p does not exceed 0.21 and 0.11 mrad, respectively. In this case, we found convergence of the FDCS already at Q = K 0 . Thus, within the framework of our approach, it seems that it is the value of the proton scattering angle that is critically important for an adequate description of the proton–electron interaction. Figure 5 and Figure 6 show the CQS and WP-CCC [25] results at, respectively, 2 and 0.5 MeV protons, for a selection of ejected energies and transferred momenta, noted here q = K 0 K . Since this is not always the case when comparing calculations of different theoretical models, it is worth noting that, here, the absolute scale is very similar; the same was observed at 1 MeV [5]. The largest differences are observed at large values of the transferred momentum. It is especially noteworthy that the FDCS obtained by the WP-CCC method does not give the usual sharp two-peak structure, e.g., for q = 1.75 at E = 2.5 or 5 or 10 eV. In the CQS results, the main peak is always surrounded by two clear minima (which are practically zero FDCS). Note that from the discussion presented above, in our calculations, the p e potential is taken into account with the best accuracy achievable within the framework of our approach.

4. Summary

In this paper, the CQS approach is applied to the calculation of the FDCS for proton–helium ionization at an impact energy of 75 keV. In contrast to the case of high projectile energy, in this intermediate regime, we are faced with the problem of adequately describing the p e interaction when using the LS-type equation for the three-body system ( e , He + , p + ) . Specifically, in the general case, the transition amplitudes show a lack of convergence as the number of terms in the Sturmian representation of the proton–electron potential is increased. In order to improve the convergence rate for small η and achieve convergence in the case of large η values, we use the auxiliary proton plane wave momentum Q as a variable parameter, the optimal choice of which allows us to solve this problem.
From the performed calculations, we deduced that the critically important parameter affecting the rate of convergence of the cross-section with an increase in the size of the p e interaction matrix representation is not the transferred momentum η itself, but the value of the proton scattering angle θ p . To test this hypothesis, we also calculated the FDCS for the incident proton energies of 0.5 and 2 MeV explored in [25]. For these values, in spite of rather significant transferred momenta, there was no need to resort to varying the parameter Q in the CQS approach, since the FDCS convergence is already achieved at Q = K 0 . As for the 1 MeV case analyzed in [5], we compared our CQS cross-sections for 0.5 and 2 MeV with the WP-CCC predictions; globally, the shapes and the magnitudes are quite similar, but, for specific kinematic regimes, some discrepancies are observed, particularly in the recoil region.
For the 75 keV kinematics, considerable discrepancies in the shape of the angular distributions are observed between the experimental data and our calculated FDCS, both in the scattering and in the perpendicular plane. A similar overall picture is found when comparing the measured FDCS with other theoretical predictions. It follows from our calculations that the proper treatment of the Coulomb proton–electron potential plays a key role at intermediate energies of the incident proton. In particular, at such energies, the discussed ionization channel competes with the electron capture channel. This indicates that some loosely accounted features of the dynamics in the ( e , p + ) subsystem, such as contributions to the FDCS from virtual proton–electron bound states, might be responsible for the observed discrepancies between our CQS calculations and experiment. Thus, the goal of our further research will be to develop our approach by improving the description of proton–electron dynamics. The first step in this direction will be the explicit inclusion of the proton–electron states, similar to what the authors of Ref. [34] intend to do. Specifically, we are going to redesign the present formulation so that the role of ‘perturbation’ is played by the interaction between the projectile and the residual ion, rather than the projectile–electron potential (8), while the spectra of subsystems ( e , He + ) and ( e , p + ) are taken into account completely in terms of the corresponding CQS functions.

Author Contributions

Conceptualization, L.U.A. and Y.V.P.; methodology, K.A.K. and Y.V.P.; software, O.C., A.S.Z. and D.S.Z.; formal analysis, S.A.Z. and K.A.K.; writing—original draft preparation, S.A.Z.; writing—review and editing, L.U.A. and Y.V.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Higher Education of the Russian Federation (project no. 0818-2020-0005) and the Ministry of Education and Science of Mongolia (grant no. ShuG 2021/137).

Data Availability Statement

All data is published in open sources.

Acknowledgments

This work The research was carried out using the equipment of the Shared Facility Center “Data Center of FEB RAS” (Khabarovsk, Russia).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
LSLippmann–Schwinger equation
FDCSFully differential cross-section
3CThree-Coulomb waves
CQSConvolution of quasi-Sturmians

References

  1. Dörner, R.; Mergel, V.; Jagutzki, O.; Spielberger, L.; Ullrich, J.; Moshammer, R.; Schmidt-Böcking, H. Cold Target Recoil Ion Momentum Spectroscopy: A ‘momentum microscope’ to view atomic collision dynamics. Phys. Rep. 2000, 330, 95. [Google Scholar]
  2. Ullrich, J.; Moshammer, R.; Dorn, A.; Dörner, R.; Schmidt, L.P.H.; Schmidt-Böcking, H. Recoil-ion and electron momentum spectroscopy: Reaction-microscope. Rep. Prog. Phys. 2003, 66, 1463. [Google Scholar] [CrossRef]
  3. Ullrich, J.; Moshammer, R.; Dörner, R.; Jagutzki, O.; Mergel, V.; Schmidt-Böcking, H.; Spielberger, L. Recoil-ion momentum spectroscopy. J. Phys. B 1997, 30, 2917. [Google Scholar] [CrossRef]
  4. Gassert, H.; Chuluunbaatar, O.; Waitz, M.; Trinter, F.; Kim, H.-K.; Bauer, T.; Laucke, A.; Müller, C.; Voigtsberger, J.; Weller, M.; et al. Agreement of Experiment and Theory on the Single Ionization of Helium by Fast Proton Impact. Phys. Rev. Lett. 2016, 116, 73201. [Google Scholar] [CrossRef]
  5. Zaytsev, A.S.; Zaytseva, D.S.; Zaytsev, S.A.; Ancarani, L.U.; Kouzakov, K.A. Parabolic quasi-Sturmian approach to proton-impact ionization of helium. Phys. Rev. A 2022, 105, 62818. [Google Scholar] [CrossRef]
  6. Voitkiv, A.B. Single ionization of helium by 1-MeV protons. Phys. Rev. A 2017, 95, 32708. [Google Scholar] [CrossRef]
  7. Abdurakhmanov, I.B.; Kadyrov, A.S.; Bray, I.; Bartschat, K. Wave-packet continuum-discretization approach to single ionization of helium by antiprotons and energetic protons. Phys. Rev. A 2017, 96, 22702. [Google Scholar] [CrossRef]
  8. Chuluunbaatar, O.; Zaytsev, S.A.; Kouzakov, K.A.; Galstyan, A.; Shablov, V.L.; Popov, Y.V. Fully differential cross sections for singly ionizing 1-MeV p+He collisions at small momentum transfer: Beyond the first Born approximation. Phys. Rev. A 2017, 96, 42716. [Google Scholar] [CrossRef]
  9. Sarkadi, L. Fully differential cross sections for the single ionization of helium by fast ions: Classical model calculations. Phys. Rev. A 2018, 97, 42703. [Google Scholar] [CrossRef]
  10. Igarashi, A.; Gulyás, L. Effective fully differential cross-section in single ionization of helium by fast ions. J. Phys. B At. Mol. Opt. Phys. 2019, 52, 245203. [Google Scholar] [CrossRef]
  11. Chuluunbaatar, O.; Kouzakov, K.A.; Zaytsev, S.A.; Zaytsev, A.S.; Shablov, V.L.; Popov, Y.V.; Gassert, H.; Waitz, M.; Kim, H.-K.; Bauer, T.; et al. Single ionization of helium by fast proton impact in different kinematical regimes. Phys. Rev. A 2019, 99, 62711. [Google Scholar] [CrossRef]
  12. Amiri Bidvari, S.; Fathi, R. Applying a four-body continuum-distorted-wave formalism in calculating the fully differential cross section for single ionization of helium in collision with high-energy protons. Eur. Phys. J. Plus 2021, 136, 190. [Google Scholar] [CrossRef]
  13. Schulz, M.; Hasan, A.; Maydanyuk, N.V.; Foster, M.; Tooke, B.; Madison, D.H. Kinematically complete experiment on single ionization in 75-keV p+He collisions. Phys. Rev. A 2006, 73, 62704. [Google Scholar] [CrossRef]
  14. Belkic, D.; Gayet, R.; Salin, A. Electron capture in high-energy ion-atom collisions. Phys. Rep. 1979, 56, 279. [Google Scholar] [CrossRef]
  15. Crothers, D.S.F.; Dube, L.J. Continuum Distorted Wave Methods in Ion–Atom Collisions. Adv. At. Mol. Opt. Phys. 1992, 30, 287–337. [Google Scholar]
  16. Ciappina, M.F.; Cravero, W.R.; Schulz, M. Post-collisional effects on single ionization in 75 keV p+He collisions. J. Phys. B At. Mol. Opt. Phys. 2007, 40, 2577. [Google Scholar] [CrossRef]
  17. Jana, D.; Samaddar, S.; Purkait, K.; Purkait, M. Fully differential cross sections for single ionization of helium by proton impact. Eur. Phys. D 2021, 75, 164. [Google Scholar] [CrossRef]
  18. Ma, X.Y.; Li, X.; Sun, S.Y.; Jia, X.F. Structures in fully differential cross-sections for proton impact ionization of helium. Europhys. Lett. 2012, 98, 53001. [Google Scholar] [CrossRef]
  19. Duan, Y.H.; Sun, S.Y.; Jia, X.F. Internuclear interaction on fully differential cross-sections for single ionization of helium by proton impact. Europhys. Lett. 2015, 110, 13001. [Google Scholar] [CrossRef]
  20. Niu, X.; Sun, S.; Wang, F.; Jia, X. Dynamic correlation effects in fully differential cross sections for 75-keV proton-impact ionization of helium. Phys. Rev. A 2017, 96, 22703. [Google Scholar] [CrossRef]
  21. Gulyás, L.; Egri, S.; Igarashi, A. Theoretical investigation of the fully differential cross sections for single ionization of He in collisions with 75-keV protons. Phys. Rev. A 2019, 99, 032704. [Google Scholar] [CrossRef]
  22. Thomas, L.H. On the Capture of Electrons by Swiftly Moving Electrified Particles. Proc. R. Soc. Lond. Ser. B 1927, 114, 561. [Google Scholar]
  23. Stepantsov, I.S.; Volobuev, I.P.; Popov, Y.V. Comparative analysis of the Compton ionization of hydrogen and positronium. Eur. Phys. D 2022, 76, 30. [Google Scholar] [CrossRef]
  24. Schulz, M.; Vajnai, T.; Gaus, A.D.; Htwe, W.; Madison, D.H.; Olson, R.E. Absolute doubly differential single-ionization cross sections in p + He collisions. Phys. Rev. A 1996, 54, 2951. [Google Scholar] [CrossRef] [PubMed]
  25. Abdurakhmanov, I.B.; Kadyrov, A.S.; Alladustov, S.U.; Bray, I.; Bartschat, K. Fully differential cross sections for single ionization of helium by energetic protons. Phys. Rev. A 2019, 100, 62708. [Google Scholar] [CrossRef]
  26. Merkuriev, S.P.; Faddeev, L.D. Quantum Scattering Theory for Several Particle Systems, 1st ed.; Springer: Dordrecht, The Netherlands, 1993. [Google Scholar]
  27. Zaytsev, S.A.; Ancarani, L.U.; Zaytsev, A.S.; Kouzakov, K.A. A parabolic quasi-Sturmian approach to quantum scattering by a Coulomb-like potential. Eur. Phys. J. Plus 2020, 135, 655. [Google Scholar] [CrossRef]
  28. Bransden, B.H.; McDowell, M.R.C. Charge Exchange and the Theory of Ion-Atom Collisions; Oxford University Press: New York, NY, USA, 1992. [Google Scholar]
  29. Zaytsev, S.A.; Zaytsev, A.S.; Nasyrov, V.V.; Zaytseva, D.S.; Ancarani, L.U.; Kouzakov, K.A. Single ionization of He by energetic protons in a parabolic quasi-Sturmians approach. Atoms 2022, 10, 13. [Google Scholar] [CrossRef]
  30. Baz, A.; Zeldovich, Y.; Perelomov, A. Scattering, Reactions, and Decays, in Nonrelativistic Quantum Mechanics; Nauka: Moscow, Russia, 1966. (In Russian) [Google Scholar]
  31. Shakeshaft, R. Integral representation of the Coulomb Green function derived from the Sturmian expansion. Phys. Rev. A 2004, 70, 42704. [Google Scholar] [CrossRef]
  32. Zaytsev, A.S.; Zaytseva, D.S.; Ancarani, L.U.; Zaytsev, S.A. Double ionization of helium with a convoluted quasi Sturmian approach. Eur. Phys. J. D 2019, 73, 111. [Google Scholar] [CrossRef]
  33. Chuluunbaatar, O.; Puzynin, I.V.; Vinitsky, P.S.; Popov, Y.V.; Kouzakov, K.A.; Cappello, C.D. Role of the cusp conditions in electron-helium double ionization. Phys. Rev. A 2006, 74, 14703. [Google Scholar] [CrossRef]
  34. Bastola, S.; Dhital, M.; Majumdar, S.; Hasan, A.; Lomsadze, R.; Davis, J.; Lamichhane, B.; Borbély, S.; Nagy, L.; Schulz, M. Interference effects in fully differential ionization cross sections near the velocity vatching in p+He collisions. Atoms 2022, 10, 119. [Google Scholar] [CrossRef]
Figure 1. The contour C 0 goes from to just above the real axis in the electron complex energy plane E e . The bound-state poles of G ^ C ( + ) ( 1 , 1 ; E e ) are shown by filled circles. The integration path C in Equation (22) is obtained by a negative angle φ rotation of C 0 , about some point E 0 on the positive real axis. C is the image of the rotated contour C on the proton energy plane E p . The unitarity branch cuts of G ^ C ( + ) ( 1 , 1 ; E e ) and G ^ C ( + ) ( 1 , m p ; E p ) are depicted by the zigzag lines. The symbols × show the superfluous poles E ± in the matrix elements of the proton Green’s function G ^ C ( + ) ( 1 , m p ; E p ) calculated in the basis (13).
Figure 1. The contour C 0 goes from to just above the real axis in the electron complex energy plane E e . The bound-state poles of G ^ C ( + ) ( 1 , 1 ; E e ) are shown by filled circles. The integration path C in Equation (22) is obtained by a negative angle φ rotation of C 0 , about some point E 0 on the positive real axis. C is the image of the rotated contour C on the proton energy plane E p . The unitarity branch cuts of G ^ C ( + ) ( 1 , 1 ; E e ) and G ^ C ( + ) ( 1 , m p ; E p ) are depicted by the zigzag lines. The symbols × show the superfluous poles E ± in the matrix elements of the proton Green’s function G ^ C ( + ) ( 1 , m p ; E p ) calculated in the basis (13).
Atoms 11 00124 g001
Figure 2. Convergence behavior of FDCS as more terms are included in the Sturmian representation (21) of the proton–electron potential. The kinematic conditions are E p = 75 keV, and electrons ejected in the scattering plane with E e = 5.4 eV. The arrow indicates the direction of the momentum transfer q .
Figure 2. Convergence behavior of FDCS as more terms are included in the Sturmian representation (21) of the proton–electron potential. The kinematic conditions are E p = 75 keV, and electrons ejected in the scattering plane with E e = 5.4 eV. The arrow indicates the direction of the momentum transfer q .
Atoms 11 00124 g002
Figure 3. FDCS for electrons with an energy of 5.4 eV ejected into the scattering plane in 75-keV p + H e . Solid lines present the CQS results. Dashed lines present 3C calculations. Experiment: • from Ref. [13]. The arrow indicates the direction of the momentum transfer q .
Figure 3. FDCS for electrons with an energy of 5.4 eV ejected into the scattering plane in 75-keV p + H e . Solid lines present the CQS results. Dashed lines present 3C calculations. Experiment: • from Ref. [13]. The arrow indicates the direction of the momentum transfer q .
Atoms 11 00124 g003
Figure 4. Same as Figure 3 but for electrons ejected into the perpendicular plane.
Figure 4. Same as Figure 3 but for electrons ejected into the perpendicular plane.
Atoms 11 00124 g004
Figure 5. FDCS for single ionization of helium by 2-MeV protons in different kinematical regimes in the scattering plane. The ejected electron energies and momentum transfer q values are indicated in the legend. Solid curves show the results obtained with our CQS method. The dotted lines represent the WP-CCC results [25].
Figure 5. FDCS for single ionization of helium by 2-MeV protons in different kinematical regimes in the scattering plane. The ejected electron energies and momentum transfer q values are indicated in the legend. Solid curves show the results obtained with our CQS method. The dotted lines represent the WP-CCC results [25].
Atoms 11 00124 g005
Figure 6. Same as Figure 5 but for the incident proton energy E p = 0.5 MeV.
Figure 6. Same as Figure 5 but for the incident proton energy E p = 0.5 MeV.
Atoms 11 00124 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zaytsev, A.S.; Zaytseva, D.S.; Zaytsev, S.A.; Ancarani, L.U.; Chuluunbaatar, O.; Kouzakov, K.A.; Popov, Y.V. Single Ionization of Helium by Protons of Various Energies in the Parabolic Quasi-Sturmians Approach. Atoms 2023, 11, 124. https://doi.org/10.3390/atoms11100124

AMA Style

Zaytsev AS, Zaytseva DS, Zaytsev SA, Ancarani LU, Chuluunbaatar O, Kouzakov KA, Popov YV. Single Ionization of Helium by Protons of Various Energies in the Parabolic Quasi-Sturmians Approach. Atoms. 2023; 11(10):124. https://doi.org/10.3390/atoms11100124

Chicago/Turabian Style

Zaytsev, Alexander S., Darya S. Zaytseva, Sergey A. Zaytsev, Lorenzo U. Ancarani, Ochbadrakh Chuluunbaatar, Konstantin A. Kouzakov, and Yury V. Popov. 2023. "Single Ionization of Helium by Protons of Various Energies in the Parabolic Quasi-Sturmians Approach" Atoms 11, no. 10: 124. https://doi.org/10.3390/atoms11100124

APA Style

Zaytsev, A. S., Zaytseva, D. S., Zaytsev, S. A., Ancarani, L. U., Chuluunbaatar, O., Kouzakov, K. A., & Popov, Y. V. (2023). Single Ionization of Helium by Protons of Various Energies in the Parabolic Quasi-Sturmians Approach. Atoms, 11(10), 124. https://doi.org/10.3390/atoms11100124

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop