This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Fragmentation of the 4He2 dimer by relativistic highly charged projectiles in collisions with small kinetic energy release

B. Najjari ​​​​​ 1, S. F. Zhang ​​​​​ 1, X. Ma ​​​​​ 1 and A. B. Voitkiv ​​2 [email protected] 1 Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China
2 Institute for Theoretical Physics I, Heinrich-Heine-Universität Düsseldorf, Universitätsstrasse 1, 40225 Düsseldorf, Germany
Abstract

We study theoretically the fragmentation of the helium dimer, 4He2, into singly charged ions in collisions with relativistic highly-charged projectiles. We discuss the main mechanisms driving this process with the focus on the fragmentation caused by direct ionization of both atomic sites of the dimer in a single collision with the projectile. This direct mechanism dominates the He2 \to He+ + He+ breakup events with the kinetic energies of the emerging ionic fragments below 4455 eV. We explore the energy and angular distributions of the He+ ions produced in collisions with 11 and 77 GeV/u U92+ projectiles and show that their shape is significantly affected by relativistic and higher-order effects in the interaction between the projectile and the dimer. We also show that the shape of the energy spectrum is quite sensitive to the binding energy of the He2 dimer which can be exploited for its precise determination. The contribution of the direct mechanism to the total cross section for the He2 fragmentation by 11 and 77 GeV/u U92+ was calculated to be 3.653.65 Mb and 2.42.4 Mb, respectively, representing roughly half of this cross section.

I Introduction

The study of fast ion-atom and ion-molecule collisions has already a long history. However, collisions of fast ions with very extended atomic objects, like humongous van der Waals dimers, are very rarely considered, even though compared to the ’normal’ case the physics of such collisions possesses new and interesting features.

In particular, to our knowledge there are just few papers, He2_alpha-particle -He2-rel-HCI-lett , where collisions of fast ions with the He2 and Li-He dimers were studied. In He2_alpha-particle -thesis_He2_S14+ and He2-rel-HCI-lett the fragmentation process He2 \to 22 He++2e{}^{+}+2\,\,e^{-} was explored for collisions with different projectiles and at different impact velocities while in He-Li-rel-HCI single ionization of the Li-He dimer by ultra-fast projectiles involving energy transfer between the dimer atoms was considered.

Both the He2 and Li-He dimers are spectacular pure quantum systems. For instance, in the He2 dimer the interaction between two ground-state helium atoms is so weak that it supports one bound state only, which has an extremely small binding energy (107\simeq 10^{-7} eV, e-l ). The size of this two-atomic bound system is enormous: its average bond length is 50\approx 50 Å  e-l and the dimer extends to the distances of more than 200200 Å, being the largest known ground-state diatomic molecule. The outer classical turning point in the ground state of the dimer is about 1414 Å  14A that is almost four times smaller than its average bond length showing that the He2 dimer is a quantum halo system which spends most of the time in the classically forbidden region. Because of its enormous size the inter-atomic interaction and the dimer binding energy are noticeably influenced by the Casimir-Polder retardation effect Cas-Pol , retard .

When the He2 interacts with external electromagnetic fields a multitude of various processes becomes possible. They include a fragmentation of the He2 dimer into helium ions which may occur when the dimer is irradiated by an electromagnetic wave or bombarded by charged projectiles. Helium ions repeal each other that results in a Coulomb explosion of the system. The kinetic energy of the ionic fragments, which is released in this explosion, depends on their charges and the initial distance between them. The spectra of the kinetic energy release may contain a valuable information about the structure of the initial He2 dimer as well as transient dimers (which may be formed if the time interval between the interaction with the external field and the start of the Coulomb explosion is sufficiently large). Some of these fragmentation processes have been already studied.

The fragmentation of the He2 dimer into He+ ions, caused by absorption of a single photon, was considered in He2_photon - He2_photon_icd-calc . In dimer-binding-exp the process of He2 fragmentation into He+ ions, induced by absorption of two high-frequency photons, was employed to sample the dimer wave function at inter-atomic distances R10R\gtrsim 10 Å  that enabled the authors to accurately measure the binding energy of the dimer.

The process of He2 fragmentation into singly charged ions in collisions with 150150 keV/u alpha particles and 11.3711.37 MeV/u S14+ projectiles was explored in He2_alpha-particle and He2_S14+ -thesis_He2_S14+ , respectively, where the focus was on the fragmentation events with kinetic energies of the He+ ions greatly exceeding 11 eV.

Very recently, the process of He2 fragmentation in collisions with relativistic highly charged projectiles was considered in He2-rel-HCI-lett . Unlike He2_alpha-particle and He2_S14+ -thesis_He2_S14+ , in He2-rel-HCI-lett the fragmentation events with kinetic energies of He+ ions not significantly exceeding 11 eV were explored. In particular, it was shown in He2-rel-HCI-lett that the lower-energy part of the fragmentation spectrum can be very strongly influenced by relativistic effects caused by the collision velocity approaching the speed of light and that the calculated energy spectrum is quite sensitive to the value of the dimer binding energy which can be exploited for measuring this energy.

In the present paper we continue to study the fragmentation reaction He2 \to 22 He+ ++ 22 e- in relativistic collisions with highly charged ions. Like in He2-rel-HCI-lett , we will focus here on the fragmentation mechanism in which the projectile directly ionizes both helium atoms in a single collision. At relativistic collision velocities this mechanism becomes very long-ranged and, therefore, is especially suited to probe the ground state of the He2 at inter-atomic distances R10R\gtrsim 10 a.u., where it already completely dominates the fragmentation process.

The paper is organized as follows. The next section begins with a discussion of the main fragmentation mechanisms which are present in collisions with relativistic highly charged ions. Then we consider a theoretical description of the direct fragmentation mechanism mentioned above. In section III we present our results for the fragmentation cross sections. Section IV contains main conclusions.

Atomic units (=|e|=me=1\hbar=|e|=m_{e}=1) are used throughout unless otherwise is stated.

II General consideration

II.1 The main fragmentation mechanisms

Let the He2 dimer collide with a projectile, which has a charge ZpZ_{p} and moves with a velocity 𝒗{\bm{v}} whose magnitude approaches the speed of light c137c\approx 137 a.u. At vcv\sim c the parameter η=Zp/v\eta=Z_{p}/v, characterizing the effective strength of the projectile field in the collision, is well below 11, which indicates that the field of the projectile in the collision is on overall weak rather than strong. Consequently, the breakup of the He2 will occur with a non-negligible probability only provided the number of ’steps’ in the interaction of the projectile with the constituents of the dimer is reduced to a necessary minimum. Also, in relativistic collisions with very light atoms the processes of radiative and nonradiative electron capture from the atom by the projectile are already extremely weak having cross sections which are by several orders of magnitude smaller that atomic ionization and excitation cross sections el-cpt ; eic ; we-1998 ; He-ioniz .

The above features inherent to high-energy collisions with light targets restrict the number of the main fragmentation mechanisms, which govern the breakup of the He2 dimer into He+ ions

Zp+He2Zp+He++He++2e\displaystyle Z_{p}+\text{He}_{2}\to Z_{p}+\text{He}^{+}+\text{He}^{+}+2e^{-} (1)

by high-energy charged projectiles, to four He2-rel-HCI-lett .

II.1.1 Mechanisms of the ’delayed’ fragmentation

Two of these mechanisms can be termed as ’delayed’ since the average time, which they require for the formation of the He+–He+ system, is determined by the motion of the helium nuclei and, therefore, exceeds by several orders of magnitude the time, which is needed by the projectile to traverse the region of space occupied by the initial He2 dimer. These mechanisms are also characterized by relatively large values of the kinetic energy of the He+ fragments. Both of them involve the interaction between the projectile and just one site of the He2 dimer.

a) In the first of these mechanisms the collision between the projectile and one of the helium atoms results in its double ionization. Then the He2+–He system evolves until the He2+ radiatively captures one electron from the neutral helium atom and the resulting He+– He+ system undergoes a Coulomb explosion. This mechanism – according to its steps – can be denoted as double-ionization–radiative-electron-transfer (DI-RET).

The fragmentation cross section σfrDI-RET\sigma_{\text{fr}}^{\text{DI-RET}} for this mechanism can be evaluated as σfrDI-RET=2σHe2+PRET\sigma_{\text{fr}}^{\text{DI-RET}}=2\,\sigma_{\text{He}}^{2+}\,P_{\text{RET}}, where σHe2+\sigma_{\text{He}}^{2+} is the cross section for double ionization of the helium atom by the projectile and PRETP_{\text{RET}} the probability for radiative electron transfer. In relativistic collisions with very highly charged projectiles (ZpvZ_{p}\sim v) the cross section σHe2+\sigma_{\text{He}}^{2+} reaches several tens of megabarns which in turn results in the magnitude of the σfrDI-RET\sigma_{\text{fr}}^{\text{DI-RET}} of the order of megabarn (see subsection D of section III).

Within the reflection approximation, the kinetic energy EKE_{K} of the ionic fragments, which is released in a Coulomb explosion, is related to the inter-nuclear distance RR at which the explosion started: in particular, for singly charged ions one has EK=1/RE_{K}=1/R. The experimental results of dimer-binding-exp and the theoretical consideration of we-21-reflect-appr show that in the case of He2 the reflection approximation is quite accurate up to very low energies (EK1E_{K}\sim 1 meV) where it starts to be violated by the recoil effects.

Since radiative electron transfer occurs only at small inter-nuclear distances (not exceeding a few atomic units) the DI-RET mechanism leads to fragmentation with relatively large kinetic energy release: the results of He2_alpha-particle and thesis_He2_S14+ show that in collisions with fast charged projectiles the energy spectrum of the fragments, which are produced via this mechanism, is peaked at EK9E_{K}\approx 9 eV corresponding to R3R\approx 3 a.u., and that its intensity rapidly decreases when the energy EKE_{K} departures from this value.

b) In the second ’delayed’ fragmentation mechanism the projectile also interacts with just one helium atom. Now, however, the interaction results in simultaneous ionization-excitation of this atom. In the presence of a neutral helium atom the excited helium ion can de-excite not only via spontaneous radiative decay but also by transferring energy to the neutral He atom that leads to ionization of the latter.

Such an energy transfer acquires a relatively large effective range becoming quite efficient if the He+ ion is produced in an excited state which can decay to the ground state by an electric dipole transition. Then the energy transfer occurs mainly via dipole-dipole two-center electron transitions.

To our knowledge, a radiationless relaxation of an excited atom via energy transfer to a neighbor atom, which at large inter-atomic distances occurs predominantly via the dipole-dipole interaction, was first considered theoretically in i-a-a . In the process, which was studied in i-a-a and termed there inter-atomic Auger decay (IAAD), the exited atom has initially a vacancy which is filled by its electron, but – unlike in ’normal’ (intra-atomic) Auger decay – the energy release is transferred to the neighbor atom ionizing it. The process of IAAD was shown in i-a-a to be so efficient that it can even outperform intra-atomic Auger decay.

In i-a-a also a simple approximate formula was derived for the IAAD rate Γa\Gamma_{a}. Being adapted to the case under consideration, the formula reads

Γa=3α8π(cω)4ΓrσphHe(ω)R6.\displaystyle\Gamma_{a}=\frac{3\alpha}{8\pi}\,\,\left(\frac{c}{\omega}\right)^{4}\,\,\,\frac{\Gamma_{r}\,\sigma_{\text{ph}}^{\text{He}}(\omega)}{\,R^{6}}. (2)

Here, Γr\Gamma_{r} is the rate for spontaneous radiative decay of the excited state of the He+ ion, ω\omega is the energy difference between the excited and ground states of the He+ ion which is transferred to the He atom, σphHe(ω)\sigma_{\text{ph}}^{\text{He}}(\omega) is the cross section for single ionization of the He atom by absorption of a photon with frequency ω\omega, RR is the distance between the He+ and He, and α\alpha is a numerical parameter 1\sim 1 depending on the magnetic quantum number of the excited state of the He+.

The process of IAAD in a system of two atomic particles, where relaxation via intra-atomic Auger decay is not allowed energetically, was computed in icd and was called there interatomic coulombic decay (ICD). Nowadays the term ICD is often used to denote relaxation mechanisms, where an overlap between electronic shells of the interacting particles is not required since the energy is transferred via the exchange of virtual photons between them. Accordingly, the fragmentation mechanism, in which simultaneous ionization-excitation of one of the atoms by the projectiles is followed by ICD decay, can be denoted as IE-ICD.

The process of ICD in the He+–He system competes with spontaneous radiative decay, which does not result in the production of the second helium ion. Using Eq. (2) and assuming for definiteness that the excited He+ ion was produced in a 2p2p–state we obtain for the ratio Γa/Γrβ(10/R)6\Gamma_{a}/\Gamma_{r}\approx\beta\,(10/R)^{6}, where β1\beta\sim 1. This indicates that when the distance RR between the He+ and He significantly exceeds 1010 a.u. the ICD channel becomes inefficient (in the experiment He2_photon_icd on the He2 photo fragmentation the ICD was clearly ’visible’ up to R12R\approx 12 a.u.). According to the results of He2_alpha-particle and thesis_He2_S14+ in fast collisions with charged projectiles the IE-ICD is very efficient in producing He+ ions with energies in the vicinity of EK8E_{K}\approx 8 eV, but is in essence completely inactivated above 1010 eV and below 4455 eV. As rough estimates suggest (see subsection D of section III), in relativistic collisions with very highly charged projectile the sum of the contributions σfrIE-ICD\sigma_{\text{fr}}^{\text{IE-ICD}} and σfrDI-RET\sigma_{\text{fr}}^{\text{DI-RET}} to the total fragmentation cross section may rich few megabarns.

c) The average size of the He2 dimer greatly exceeds not only the effective range of RET but also that of ICD. Therefore, there must occur a very significant contraction of the intermediate (He+)–He and He2+–He dimers before the distance between the nuclei sufficiently decreases in order for the ICD and RET to come into the play. The contraction is possible because the potential of the He2+–He and (He+)–He systems at large inter-nuclear distances is attractive.

II.1.2 Mechanisms of the ’instantaneous’ fragmentation

In collisions with high-energy projectiles there are two more fragmentation mechanisms in which, unlike the DI-RET and IE-ICD, the He+-He+ system emerges ’instantaneously’, i.e. on the so short time scale that during the production of two helium ions the helium nuclei remain essentially at rest.

a). In one of them the projectile interacts with both atoms of the dimer. As a result, each helium atom emits an electron and becomes a singly charged ion. In this fragmentation mechanism the projectile directly forms the He+ - He+ system Sulf which undergoes a Coulomb explosion. Since the size of the He2 is very large, the interactions between the constituents of the dimer play in this mechanism no noticeable role. In He2-rel-HCI-lett this mechanism was termed the direct fragmentation (DF) and we shall use it also here.

We note that at relativistic impact velocities the time interval TT between the collisions of the projectile with the first and second atoms of the He2 dimer does not exceed few atomic units (T1016T\lesssim 10^{-16} s). This is by orders of magnitude smaller than typical nuclear times (1013\sim 10^{-13}101210^{-12} s) in the He+ – He system, showing that the DF mechanism is indeed ’instantaneous’.

b). In the other ’instantaneous’ fragmentation mechanism the projectile interacts with just one atom of the dimer. As a result of this interaction, the atom emits an electron. There is a certain probability that the emitted electron will move towards the other atom and knock out one of its electrons. Thus, this fragmentation mechanism is a combination of single ionization of the helium atom by a high-energy projectile and the so called e-2e process on helium (single ionization by electron impact). Following He2-rel-HCI-lett we shall refer to this mechanism as single ionization – e-2e (SI–e-2e) f1 .

In the SI–e-2e process the velocity of the first emitted electron is by two orders of magnitude smaller than the velocity of the relativistic projectile. Nevertheless, this electron still moves much faster than the helium nuclei. With its typical velocity ve1v_{e}\sim 1 a.u., we obtain that the time, which the electron needs to propagate between the sites of the dimer, does not exceed 101410^{-14} s. Since this time is much shorter than typical nuclear times in the He+ – He system, the SI–e-2e mechanism can be also viewed as ’instantaneous’.

Since in both the DF and SI–e-2e the helium ions emerge at ’frozen’ positions of the helium nuclei, their energy EKE_{K} is very simply (EK=1/RE_{K}=1/R) related to the size RR of the He2 dimer at the collision instant (as long as the recoil effects can be ignored, see subsection B of the present section).

The cross section for the production of two singly charged helium ions by relativistic projectiles via the DF mechanism depends on the transverse size RR_{\perp} of the He2 dimer He2-rel-HCI-lett , scaling approximately as 1/R21/R_{\perp}^{2}. The corresponding cross section for the production via the SI–e-2e depends on the dimer size RR as 1/R21/R^{2} He2-rel-HCI-lett .

Thus, unlike the DI-RET and IE-ICD, the DF and SI–e-2e mechanisms are not only instantaneous, providing a simple correspondence between the kinetic energy release and the size of the He2 dimer, but also possess much longer effective range. Therefore, they are especially suited for probing the ground state of the He2 in a very large interval of the inter-atomic distances RR.

Our calculations show that in collisions of He2 dimers with very highly charged projectiles (Zp/v1Z_{p}/v\sim 1) the DF mechanism is much more efficient than the SI–e-2e yielding cross sections which are larger by two orders of magnitude. Therefore, in the rest of the paper we shall concentrate on the He2 fragmentation via the DF mechanism.

II.2 The direct fragmentation mechanism

We shall consider collisions between the He2 dimer and the projectile using the semi-classical approximation in which the relative motion of the heavy particles (nuclei) is treated classically whereas the electrons are considered quantum mechanically.

We choose a reference frame in which the dimer is at rest and take the nucleus of one of its atoms as the origin. We shall refer to this atom as atom AA whereas the other will be denoted by BB. In this frame, the coordinates of the nucleus of atom BB are given by the inter-nuclear vector 𝑹{\bm{R}} of the dimer and the projectile moves along a classical straight-line trajectory 𝐑p(t)=𝒃+𝒗t\mathbf{R}_{p}(t)=\bm{b}+\bm{v}t, where 𝒃=(bx,by,0)\bm{b}=(b_{x},b_{y},0) is the impact parameter with respect to the nucleus of atom AA and 𝒗=(0,0,v)\bm{v}=(0,0,v) the collision velocity.

Since the size of the He2 dimer is very large, the ionization of atoms AA and BB occurs independently of each other and we can use the independent electron approximation. According to it the probability PA+B+P_{A^{+}B^{+}} for single ionization of atoms AA and BB in the collision where the projectile moves with an impact parameter 𝒃\bm{b} is given by

PA+B+=PA+(𝒃)PB+(𝒃).\displaystyle P_{A^{+}B^{+}}=P_{A^{+}}(\bm{b})P_{B^{+}}(\bm{b}^{\prime}). (3)

Here, PA+(𝒃)P_{A^{+}}(\bm{b}) and PB+(𝒃)P_{B^{+}}(\bm{b}^{\prime}) are the probabilities for single ionization of atoms AA and BB, respectively, and 𝒃=𝒃𝑹{\bm{b}^{\prime}}={\bm{b}}-{\bm{R}}_{\perp} is the collision impact parameter with respect to atom BB, where 𝑹{\bm{R}}_{\perp} is the part of the inter-nuclear vector 𝑹{\bm{R}} of the dimer which is perpendicular to the projectile velocity 𝒗{\bm{v}}.

Within the independent electron approximation the probabilities for single ionization of helium atoms AA and BB read

PA+(𝒃)\displaystyle P_{A^{+}}(\bm{b}) =\displaystyle= 2w(𝒃)(1w(𝒃))\displaystyle 2\,w(\bm{b})\,(1-w(\bm{b}))
PB+(𝒃)\displaystyle P_{B^{+}}(\bm{b}^{\prime}) =\displaystyle= 2w(𝒃)(1w(𝒃)),\displaystyle 2\,w(\bm{b}^{\prime})\,(1-w(\bm{b}^{\prime})), (4)

where w(𝒃)w(\bm{b}) (w(𝒃)w(\bm{b}^{\prime})) is the probability to remove one electron from helium atom in the collision with a projectile having an impact parameter 𝒃\bm{b} (𝒃\bm{b}^{\prime}) new-approach .

The quantity

σA+B+DF\displaystyle\sigma^{\text{DF}}_{A^{+}B^{+}} =\displaystyle= d2𝒃PA+B+(𝒃)\displaystyle\int d^{2}{\bm{b}}\,P_{A^{+}B^{+}}(\bm{b}) (5)
=\displaystyle= d2𝒃PA+(𝒃)PB+(𝒃)\displaystyle\int d^{2}{\bm{b}}\,P_{A^{+}}(\bm{b})\,P_{B^{+}}(\bm{b}^{\prime})
=\displaystyle= d2𝒃PA+(𝒃)PB+(𝒃𝑹)\displaystyle\int d^{2}{\bm{b}}\,P_{A^{+}}(\bm{b})\,P_{B^{+}}({\bm{b}}-{\bm{R}}_{\perp})

represents the cross section for the production of two singly charged helium atom by the projectile in the collision with the He2 dimer at a given inter-nuclear vector 𝑹{\bm{R}} of the latter. Since the probabilities PA+(𝒃)P_{A^{+}}(\bm{b}) and PB+(𝒃)P_{B^{+}}(\bm{b}^{\prime}) depend just on the absolute value of the respective impact parameters, PA+(𝒃)PA+(b)P_{A^{+}}(\bm{b})\equiv P_{A^{+}}(b) and PB+(𝒃)PB+(b)P_{B^{+}}(\bm{b}^{\prime})\equiv P_{B^{+}}(b^{\prime}), the cross section (5) depends only on the absolute value RR_{\perp} of the two-dimensional vector 𝑹{\bm{R}}_{\perp}: σA+B+DF=σA+B+DF(R)\sigma^{\text{DF}}_{A^{+}B^{+}}=\sigma^{\text{DF}}_{A^{+}B^{+}}(R_{\perp}).

In our calculations of the probabilities w(b)w(b) and w(b)w(b^{\prime}) we regarded the helium atom in its initial and final states as an effectively single-electron system, where the ’active’ electron moves in the effective field, created by the ’frozen’ atomic core consisting of the atomic nucleus and the ’passive’ electron. This field is approximated by the potential

V(𝒓)=1r(1+βr)exp(αr)r,\displaystyle V({\bm{r}})=-\frac{1}{r}-\left(1+\beta r\right)\frac{\exp(-\alpha r)}{r},

where rr is the distance between the active electron and the atomic nucleus, and α=3.36\alpha=3.36 and β=1.665\beta=1.665 param . We note that with this choice of α\alpha and β\beta, V(𝒓)V({\bm{r}}) almost coincides with the exact Hartree-Fock potential.

The cross section (5) can be only computed numerically and the computation is quite time consuming. However, a simple estimate can be obtained for this cross section at R1R_{\perp}\gg 1 a.u. He2-rel-HCI-lett

σDFCZp4v4γ2[K12(RRa)+1γ2K02(RRa)].\displaystyle\sigma^{\text{DF}}\approx C\,\frac{Z_{p}^{4}}{v^{4}\gamma^{2}}\,\,\bigg{[}K^{2}_{1}\left(\frac{R_{\perp}}{R_{a}}\right)+\frac{1}{\gamma^{2}}\,K^{2}_{0}\left(\frac{R_{\perp}}{R_{a}}\right)\bigg{]}. (6)

Here, K0K_{0} and K1K_{1} are the modified Bessel function A-S , Ra=γvω¯R_{a}=\frac{\gamma v}{\overline{\omega}} is the adiabatic collision radius, where γ=1/1v2/c2\gamma=1/\sqrt{1-v^{2}/c^{2}} is the collision Lorentz factor, ω¯1.2\overline{\omega}\approx 1.2 a.u. is the mean transition frequency for single ionization of a helium atom, and CC is a free (fit) parameter which is a very slowly varying function of RR_{\perp} param-C .

Taking into account that K0(x)ln(1.12/x)K_{0}(x)\sim\ln(1.12/x) and K1(x)1/xK_{1}(x)\sim 1/x if x<1x<1 whereas K0(x)K1(x)π2xexp(x)K_{0}(x)\sim K_{1}(x)\sim\sqrt{\frac{\pi}{2x}}\exp(-x) at x>1x>1, it follows from Eq. (6) that at 1RRa1\ll R_{\perp}\lesssim R_{a} the cross section decreases with RR_{\perp} relatively slowly, σDFR2\sigma^{\text{DF}}\sim R^{-2}_{\perp}, whereas at R>RaR_{\perp}>R_{a} the decrease is already exponential. This means that the projectile is able to efficiently irradiate both atoms of the dimer only provided that the dimer transverse size RR_{\perp} does not exceed the adiabatic collision radius RaR_{a}. Since RaγvR_{a}\sim\gamma v, an ultrafast projectile possesses a very large effective interaction range and can probe systems with large dimensions. Simple estimates show that, beginning with impact energies of a few GeV/u, the adiabatic collision radius exceeds the size of the He2 dimer He2-rel-HCI-lett .

The cross section (5) refers to the situation where two helium atoms in the dimer are separated by the vector 𝑹{\bm{R}} with the absolute value of its transverse projection RR_{\perp}. It yields an important information about the collisions giving an insight into the basic physics of the fragmentation process He2-rel-HCI-lett but cannot be (directly) measured in experiment. Therefore, let us ’convert’ the cross section (5) into quantities which can be measured. We start with the expression

dσfrDFd3𝑹=σA+B+DF(R)|Ψi(𝑹)|2,\displaystyle\frac{d\sigma^{\text{DF}}_{\text{fr}}}{d^{3}{\bm{R}}}=\sigma^{\text{DF}}_{A^{+}B^{+}}(R_{\perp})\,\,\big{|}\Psi_{i}({\bm{R}})\big{|}^{2}, (7)

where Ψi(𝑹)\Psi_{i}({\bm{R}}) is the wave function describing the relative motion of the atoms in the ground state of the He2 dimer.

As was already mentioned, we assume that before the collision with the projectile the He2 dimer was at rest. Let 𝑷Arec{\bm{P}}^{\text{rec}}_{A} and 𝑷Brec{\bm{P}}^{\text{rec}}_{B} be the recoil momenta of the helium ions, which they acquire during the collision because of the interaction with the projectile and electron emission, and let 𝑷A{\bm{P}}_{A} and 𝑷B{\bm{P}}_{B} be the final momenta of these ions after the Coulomb explosion. Then the energy conservation for the relative motion of the ionic fragments reads

EK=[12(𝑷Arec𝑷Brec)]22μ+QAQBR.\displaystyle E_{K}=\frac{\big{[}\frac{1}{2}\left({\bm{P}}^{\text{rec}}_{A}-{\bm{P}}^{\text{rec}}_{B}\right)\big{]}^{2}}{2\mu}+\frac{Q_{A}Q_{B}}{R}. (8)

Here, EK=𝑲2/(2μ)E_{K}={\bm{K}}^{2}/(2\mu) is the final kinetic energy of the relative motion of the ionic fragments, 𝑲=(𝑷A𝑷B)/2{\bm{K}}=\left({\bm{P}}_{A}-{\bm{P}}_{B}\right)/2 and μ\mu are the momentum of their relative motion and the reduced mass, respectively, QAQ_{A} and QBQ_{B} (QA=QB=1Q_{A}=Q_{B}=1) are their charges and RR is the distance between the ions when the Coulomb explosion began and which in the case under consideration coincides with the size of the initial He2 dimer at the collision instant. In the energy balance (8) we have neglected the kinetic energy of the He+ ions, which they had due to the nuclear motion before the collision and which is very small because the depth of the potential well in the He2 dimer is just 1\approx 1 meV.

Since in the process of single ionization of helium atoms by relativistic projectiles the recoil momenta of the helium ions do not noticeably exceed 11 a.u. rhci -r-sea , the first term on the right hand side of Eq. (8) is in the meV range. Such values are comparable to the Coulomb potential energy QAQB/RQ_{A}Q_{B}/R at the inter-nuclear distances RR of the order of 10410^{4} a.u. which is much larger than the size of the He2 (and any known) dimer. Therefore, the neglect of the recoil energy, given by the first term on the right-hand side of Eq. (8), does not have any substantial impact on the energy balance. Thus, the relation

EK=QAQBR\displaystyle E_{K}=\frac{Q_{A}Q_{B}}{R} (9)

between the kinetic energy release and the size of the dimer at the collision instant, which neglects the recoil energies, is expected to be very accurate down to energies EK10E_{K}\simeq 10 meV.

Let us now make a more restrictive assumption that not only the recoil energies are much smaller than the Coulomb energy QAQB/RQ_{A}Q_{B}/R but also that the absolute values of the recoil momenta, |𝑷Arec||{\bm{P}}^{\text{rec}}_{A}| and |𝑷Brec||{\bm{P}}^{\text{rec}}_{B}|, of the fragments are significantly less than K=|𝑲|K=|{\bm{K}}|. Since these values do not exceed 11 a.u. the assumption will be fulfilled beginning with K4K\gtrsim 455 a.u. that corresponds to the energies EK60E_{K}\gtrsim 60 meV. Under such conditions the momentum 𝑲\bm{K} will be directed essentially along the inter-nuclear vector 𝑹\bm{R} of the initial He2 dimer and, taking also into account (9), we get

d3𝑹\displaystyle d^{3}\bm{R} =\displaystyle= (QAQB)3μKEK4d3𝑲.\displaystyle\frac{(Q_{A}Q_{B})^{3}}{\mu\,K\,E^{4}_{K}}\,\,d^{3}{\bm{K}}. (10)

Using (7) and (9)-(10) we obtain the fragmentation cross section differential in the relative momentum 𝑲{\bm{K}} of the ionic fragments

dσfrDFd3𝑲\displaystyle\frac{d\sigma^{\text{DF}}_{\text{fr}}}{d^{3}{\bm{K}}} =\displaystyle= (QAQB)3μKEK4σA+B+DF(R(𝑲))\displaystyle\frac{(Q_{A}Q_{B})^{3}}{\mu\,K\,E^{4}_{K}}\,\,\,\sigma^{\text{DF}}_{A^{+}B^{+}}\big{(}R_{\perp}(\bm{K})\big{)}\, (11)
×|Ψi(QAQBEK𝑲^)|2.\displaystyle\times\bigg{|}\Psi_{i}\bigg{(}\frac{Q_{A}Q_{B}}{E_{K}}\hat{\bm{K}}\bigg{)}\bigg{|}^{2}.

where 𝑲^=𝑲/K\hat{\bm{K}}={\bm{K}}/K, R(𝑲)=(QAQBsinΘ𝑲)/EKR_{\perp}(\bm{K})=(Q_{A}Q_{B}\,\sin\Theta_{\bm{K}})/E_{K} and Θ𝑲\Theta_{\bm{K}} is the polar angle of the momentum 𝑲\bm{K} (the zz-axis is along the projectile velocity 𝒗\bm{v}). Since the bound state of the He2 dimer is spherically symmetric, the wave function Ψi\Psi_{i} does not depend on 𝑲^\hat{\bm{K}}: Ψi(QAQBEK𝑲^)=Ψi(QAQBEK)\Psi_{i}\bigg{(}\frac{Q_{A}Q_{B}}{E_{K}}\hat{\bm{K}}\bigg{)}=\Psi_{i}\bigg{(}\frac{Q_{A}Q_{B}}{E_{K}}\bigg{)}.

Taking into account that d3𝑲=K2dKdΩ𝑲=μKdEKsinΘ𝑲d𝑲dφ𝑲d^{3}{\bm{K}}=K^{2}dK\,d\Omega_{\bm{K}}=\mu KdE_{K}\,\sin\Theta_{\bm{K}}d{\bm{K}}d\varphi_{\bm{K}}, where φ𝑲\varphi_{\bm{K}} is the azimuthal angle of 𝑲\bm{K}, and that the right-hand side of Eq. (11) does not depend on φ𝑲\varphi_{\bm{K}} we obtain the fragmentation cross section differential in the kinetic energy release EKE_{K} and the polar angle Θ𝑲\Theta_{\bm{K}}

dσfrDFdEKdΘ𝑲\displaystyle\frac{d\sigma^{\text{DF}}_{\text{fr}}}{dE_{K}d\Theta_{\bm{K}}} =\displaystyle= 2π(QAQB)3EK4|Ψi(QAQBEK)|2\displaystyle 2\pi\,\frac{(Q_{A}Q_{B})^{3}}{E^{4}_{K}}\,\bigg{|}\Psi_{i}\bigg{(}\frac{Q_{A}Q_{B}}{E_{K}}\bigg{)}\bigg{|}^{2} (12)
×\displaystyle\times sinΘ𝑲σA+B+DF(QAQBsinΘ𝑲/EK).\displaystyle\sin\Theta_{\bm{K}}\,\,\sigma^{\text{DF}}_{A^{+}B^{+}}\big{(}Q_{A}Q_{B}\sin\Theta_{\bm{K}}/E_{K}\big{)}.

In our consideration of the direct fragmentation process, its two steps – ’instantaneous’ production of two He+ ions and the consequent Coulomb explosion of the He+–He+ system – are fully disentangled, even though the energy of the emitted electrons and the kinetic energy release have the same source: the energy transferred by the projectile to the He2 dimer. However, since the electrons emitted from helium atoms in collisions with relativistic projectiles have energies extending up to a few tens of eV whereas, as will be seen below, the spectrum of the kinetic energy release in the DF process peaks at EK<1E_{K}<1 eV and essentially vanishes already at EK5E_{K}\simeq 566 eV, such a consideration is expected to be quite accurate.

III Results and Discussion

III.1 Preliminary remarks

In this section we report our results for the fragmentation cross sections. In our calculations the ground state of the He2 dimer was described by a wave function Ψi\Psi_{i} corresponding to the dimer binding energy IbI_{b} of 139139 neV (this value is very close to that given in 139.2 ): this wave function was used to get results presented below in figures 1-7. In addition, in order to explore the sensitivity of the cross sections to the value of IbI_{b}, also wave functions Ψi\Psi_{i} corresponding to Ib=130I_{b}=130 neV and 148148 neV were applied for obtaining results shown in figure 8.

Four different approximations were employed for computing the probabilities ww for single electron removal from a helium atom (see Eq. (4)). They include: i) the nonrelativistic first Born approximation (nr-FBA); ii) the nonrelativistic symmetric eikonal approximation (nr-SEA) nr-sea ; iii) the relativistic first Born approximation (r-FBA) and iv) the relativistic symmetric eikonal approximation (r-SEA) r-sea .

The differences between results obtained by using the first Born and symmetric eikonal approximations offer an idea about the importance of the higher-order effects in the interaction between the projectile and the atoms of the dimer. On the other hand, the differences between results of the relativistic and nonrelativistic approximations yield the information about the role of relativistic effects in the He2 fragmentation.

In particular, in the nr-FBA and nr-SEA the speed of light cc is assumed to be infinite which means that all relativistic effects caused by very high collision velocities vanish. In the r-FBA and r-SEA c137c\approx 137 a.u. and relativistic effects arising due to high impact energies are taken into account. Since even in very high-energy collisions with helium targets the overwhelming majority of emitted electrons have energies well below 100100 eV we-jpb-2005 , the main relativistic effects, which are due to high impact energies and which may influence the He2 fragmentation, are caused by the deviation of the electric field of the projectile from the (unretarded) Coulomb form.

One should say that, in addition to relativistic effects due to high impact velocities, there are also relativistic effects in the ground state of the free He2 which – according to 139.2 – reduce the binding energy of the dimer by about 14%14\%. This affects the form of its wave function (especially at very large inter-atomic distances) that in turn may have an impact on the shape of the spectrum of the fragments He2-rel-HCI-lett .

III.2 Fragmentation spectra

Refer to caption
Figure 1: The cross section dσfrDFdEKdΘ𝑲\frac{d\sigma^{\text{DF}}_{\text{fr}}}{dE_{K}d\Theta_{\bm{K}}} for the fragmentation of the He2 dimer by 11 GeV/u U92+ projectiles given as a function of the kinetic energy release EKE_{K} and the angle Θ𝑲\Theta_{\bm{K}}. The upper and lower panels shown results obtained using the r-FBA and r-SEA, respectively.
Refer to caption
Figure 2: Same as in figure 1 but for the fragmentation by 77 GeV/u U92+.

In figures 1 and 2 we show the doubly differential cross section dσfrDFdEKdΘ𝑲\frac{d\sigma^{\text{DF}}_{\text{fr}}}{dE_{K}d\Theta_{\bm{K}}} for the fragmentation of the He2 dimer into He+ ions occurring in collisions with 11 GeV/u and 77 GeV/u U92+ projectiles, respectively. The corresponding collision velocities and the Lorentz factor are v=120v=120 a.u., γ=2.08\gamma=2.08 and v=136v=136 a.u., γ=8.5\gamma=8.5, respectively. The cross section is plotted as a function of the kinetic energy release EKE_{K} and the fragmentation angle Θ𝑲\Theta_{\bm{K}} presenting a general picture of the fragmentation spectrum. Taking into account that the spectrum is symmetric with respect to the transformation Θ𝑲πΘ𝑲\Theta_{\bm{K}}\leftrightarrow\pi-\Theta_{\bm{K}}, only the interval of angles 0Θ𝑲π/20\leq\Theta_{\bm{K}}\leq\pi/2 is considered.

It is seen in these figures that the spectrum is predominantly localized in the energy interval 0.10.1 eV EK2\lesssim E_{K}\lesssim 2 eV and at fragmentation angles Θ𝑲\Theta_{\bm{K}} not significantly exceeding 2020^{\circ}. It is also seen that at very small energies EKE_{K} the spectrum is restricted to smaller fragmentation angles Θ𝑲\Theta_{\bm{K}} whereas with increasing the energy the spectrum shifts to larger Θ𝑲\Theta_{\bm{K}}.

The interval of kinetic energies 0.10.1 eV EK2\lesssim E_{K}\lesssim 2 eV corresponds to the instantaneous size RR of the He2 dimer during the collision in the range 1414 a.u. R270\lesssim R\lesssim 270 a.u. The probability to find the He2 dimer in this range is above 90%90\% that is reflected in the dominance of the above energy interval in the spectrum. At energies 5\gtrsim 566 eV the spectrum practically vanishes since they correspond to the inter-nuclear distances R4R\lesssim 4 a.u., where the probability to find the He2 dimer is negligibly small.

The cross section σA+B+DF\sigma^{\text{DF}}_{A^{+}B^{+}} for the production of two helium ions by the projectile decreases with the transverse size RR_{\perp} of the dimer at the collision instant (σA+B+DF1/R2\sigma^{\text{DF}}_{A^{+}B^{+}}\sim 1/R_{\perp}^{2} at R<RaR_{\perp}<R_{a}). Since at a fixed Θ𝐊\Theta_{\bf K} the value of RR_{\perp} is proportional to RR, the fragmentation with very small EKE_{K} may occur at very small Θ𝐊\Theta_{\bf K} only, whereas at significantly larger values of EKE_{K} it becomes possible at not very small Θ𝐊\Theta_{\bf K} as well. Taking also into account that the cross section dσfrDFdEKdΘ𝑲\frac{d\sigma^{\text{DF}}_{\text{fr}}}{dE_{K}d\Theta_{\bm{K}}} is proportional to the geometrical factor sinΘ𝑲\sin\Theta_{\bm{K}} we can explain the ’shift’ of the spectrum to larger Θ𝑲\Theta_{\bm{K}} with increase in EKE_{K}, which is observed in figures 1 and 2.

It also follows from the results presented in figures 1 and 2 that for collisions with so very highly charged projectiles, as U92 ions, the symmetric eikonal approximation predicts a noticeably smaller number of the fragmentation events. This shows that the higher-order effects in the interaction between the projectile and the target remain visible even at very high impact energies. The smaller cross section values, predicted by the symmetric eikonal approximation, reflects the general tendency that in high-energy collisions the higher-order effects reduce cross section values.

Refer to caption
Figure 3: The cross section dσfrDFdEKdΘ𝑲\frac{d\sigma^{\text{DF}}_{\text{fr}}}{dE_{K}d\Theta_{\bm{K}}} for the fragmentation of the He2 dimer by 11 GeV/u U92+ projectiles given as a function of the angle Θ𝑲\Theta_{\bm{K}} at fixed values of the energy EKE_{K}: 6060 meV (aa), 100100 meV (bb), 11 eV (cc) and 22 eV (dd). Dotted, dashed, dash-dotted and solid curves are results of the calculations in the nr-FBA, the nr-SEA, the r-FBA and the r-SEA, respectively.
Refer to caption
Figure 4: Same as in figure 3 but for collisions with 77 GeV/u U92+ projectiles.

Figures 3 and 4 give a more detailed information about the fragmentation process by displaying the cross section dσfrDFdEKdΘ𝑲\frac{d\sigma^{\text{DF}}_{\text{fr}}}{dE_{K}d\Theta_{\bm{K}}} as a function of the fragmentation angle Θ𝑲\Theta_{\bm{K}} at a given value of the kinetic energy release EKE_{K}: 6060 meV, 100100 meV, 11 eV and 22 eV. Some conclusions can be drawn from the results presented in these figures.

First, by comparing results of the different calculations we see that the shape of the He2 fragmentation pattern is in general influenced by both relativistic effects, which arise due to collision velocities approaching the speed of light, and higher-order effects, which originate in a very high charge of the projectile.

Second, the influence of the relativistic effects increases when the kinetic energy release decreases: these effects especially impact the fragmentation events with very small kinetic energy release where they enhance the cross section up to one or two orders of magnitude. On the other hand, the fragmentation with larger values of EKE_{K} is only very modestly influenced by these effects.

Third, the situation with the higher-order effects is just opposite: their role rises with increasing EKE_{K} and they mostly affect the fragments having larger energies. However, even for such fragments these effects remains modest, reducing the cross section values not significantly more than 1010-20%20\%. The higher-order effects also somewhat influence the low-energy part of the spectrum: in this case, however, their effect is limited to the very small angles Θ𝑲\Theta_{\bm{K}} only.

The above qualitative features of the fragmentation process can be understood by noting that the cross section σA+B+DF(R)\sigma^{\text{DF}}_{A^{+}B^{+}}(R_{\perp}) is most profoundly influenced by the relativistic effects when the magnitude of RR_{\perp} is large (Ra/γ<RRaR_{a}/\gamma<R_{\perp}\lesssim R_{a}) whereas the higher-order effects are most significant in collisions where the projectile passes close to both atoms of the dimer (the impact parameters bb and bb^{\prime} do not exceed a few atomic units) that is possible only at relatively small values of RR_{\perp}.

Very large values of RR_{\perp} imply that the instantaneous dimer size RR in the collision has to be very large and, in addition, the dimer orientation angle Θ𝐑\Theta_{\bf R} should not be small. Since EK=1/RE_{K}=1/R and Θ𝐑Θ𝐊\Theta_{\bf R}\approx\Theta_{\bf K}, it follows that the relativistic effects are most significant for fragmentation events, which are characterized by a very small kinetic energy release and not too small angles Θ𝑲\Theta_{\bm{K}}.

On the other hand, collisions occurring at small values of RR_{\perp}, where the higher-order effects can be substantial, are characterized by either small values of the instantaneous dimer size RR or very small angles Θ𝑹Θ𝑲\Theta_{\bm{R}}\approx\Theta_{\bm{K}} (or a combination of both). The kinetic energy release EKE_{K} of 0.060.06 eV, 0.10.1 eV, 11 eV and 22 eV corresponds to R453R\approx 453 a.u., 272272 a.u., 27.227.2 a.u. and 13.613.6 a.u., respectively. In order to have substantial higher-order effects at the first two values of RR the dimer has to be almost parallel to the collision velocity. As a result, in collisions with very small EKE_{K} such effects can become visible only at the fragmentation angles Θ𝑲\Theta_{\bm{K}} close to zero whose contribution, due to the geometric factor sinΘ𝑲\sin\Theta_{\bm{K}}, is very small. For collisions with larger EKE_{K} the restrictions on the angle Θ𝑲\Theta_{\bm{K}} are not so strict and the higher-order effects may be noticeable for the whole range 0Θ𝑲900^{\circ}\leq\Theta_{\bm{K}}\leq 90^{\circ}.

The energy EKE_{K} of 0.060.06 eV, 0.10.1 eV, 11 eV and 22 eV correspond to the relative momentum KK of 4\approx 4 a.u., 5.25.2 a.u., 16.416.4 a.u. and 23.223.2 a.u., respectively. The first two momentum values are not very large. Therefore, in the case of EK=0.06E_{K}=0.06 eV and 0.10.1 eV the very narrow peak at very small angles (see panels (aa) and (bb) of figures 34) can be noticeably smeared out by the recoil effects.

Refer to caption
Figure 5: The angular distribution, dσfrDFdΘ𝑲\frac{d\sigma^{\text{DF}}_{\text{fr}}}{d\Theta_{\bm{K}}}, of the He+ fragments produced in collisions with 11 GeV/u (two upper panels) and 77 GeV/u (two lower panels) U92+ projectiles. The results shown in panels (aa), (cc) and (bb),(dd) were obtained by integrating over the energy intervals 0 eV EK2\leq E_{K}\leq 2 eV and 0 eV EK7\leq E_{K}\leq 7 eV, respectively. Dotted, dash-dotted, dash and solid curves represent results of the calculations performed in the nr-FBA, nr-SEA, r-FBA and r-SEA approximations, respectively.

In figure 5 we display the angular distribution of the He+ fragments represented by the cross section dσfrDFdΘ𝑲\frac{d\sigma^{\text{DF}}_{\text{fr}}}{d\Theta_{\bm{K}}}. The cross section was obtained by integrating either over the energy interval 0 eV EK2\leq E_{K}\leq 2 eV, in which the direct fragmentation mechanism is by far the dominant one, or over the broader range 0 eV EK7\leq E_{K}\leq 7 eV, which covers essentially all energies which are possible in the fragmentation events caused by this mechanism (see figures 1-2).

Figure 5 shows that the angular distribution of the He+ fragments becomes significantly weaker dependent on the angle Θ𝑲\Theta_{\bm{K}} when their energy EKE_{K} increases. This feature (which is already seen in figures 34) can be understood by taking into account that, due to the relation EK=1/RE_{K}=1/R, the He+ fragments with larger energies emerge when the projectiles hit the He2 dimers with smaller size RR. In such collisions the DF mechanism is more efficient causing the He2 breakup with significant probabilities even when the dimer is oriented at large angles with respect to the impact velocity.

Refer to caption
Figure 6: The energy spectrum, dσfrDFdEK\frac{d\sigma^{\text{DF}}_{\text{fr}}}{dE_{K}}, in the fragmentation of the He2 dimer in collisions with 11 GeV/u (the upper panels) and 77 GeV/u (the lower panels) U92+ projectiles. The results shown in panels (aa), (cc) and (bb), (dd) were obtained by integrating over the angles 0Θ𝑲1800^{\circ}\leq\Theta_{\bm{K}}\leq 180^{\circ} and 60Θ𝑲12060^{\circ}\leq\Theta_{\bm{K}}\leq 120^{\circ}, respectively. Dotted, dash-dotted, dash and solid curves represent results of the calculations performed in the nr-FBA, nr-SEA, r-FBA and r-SEA, respectively.

The energy spectrum, dσfrDF/dEKd\sigma^{\text{DF}}_{\text{fr}}/dE_{K}, of He+ fragments produced by bombarding the He2 dimer by 11 and 77 GeV/u U92+ projectiles is displayed in figure 6. The spectrum shown in the left panels of this figure was obtained by integrating over all possible fragmentation angles, 0Θ𝑲1800^{\circ}\leq\Theta_{\bm{K}}\leq 180^{\circ}, whereas that given in the right panels includes only the fragments moving at large angles, 60Θ𝑲12060^{\circ}\leq\Theta_{\bm{K}}\leq 120^{\circ}, with respect to the projectile velocity.

The results shown in figure 6 confirm the correspondences between the relativistic and higher-order effects and the energy EKE_{K} and angle Θ𝑲\Theta_{\bm{K}} of the fragments which were discussed above. In particular, we again observe that at a given collision velocity the relativistic effects rise when the energy EKE_{K} decreases, becoming especially strong in collisions where the He+ fragments possess quite small energies EKE_{K} and move at large angles with respect to the projectile velocity.

Refer to caption
Figure 7: The energy spectrum, dσfrDFdEK\frac{d\sigma^{\text{DF}}_{\text{fr}}}{dE_{K}}, for the He2 breakup into He+ via the DF mechanism calculated within the r-SEA. Results for collisions with 11 GeV/u U92+, 77 GeV/u U92+ and 11.3711.37 MeV/u S14+ are shown by solid, dashed and dotted curves, respectively. The results shown in (aa) and (bb) were obtained by integrating over 0Θ𝑲1800^{\circ}\leq\Theta_{\bm{K}}\leq 180^{\circ} and 45Θ𝑲13545^{\circ}\leq\Theta_{\bm{K}}\leq 135^{\circ}, respectively.

In figure 7 we compare the energy spectra of the He+ ions produced in collisions with 11 and 77 GeV/u U92+ and 11.3711.37 MeV/u S14+ (v=21.2v=21.2 a.u., γ=1.01\gamma=1.01). It is seen that the shape of the spectrum profoundly varies with an increase in the impact energy. In particular, the maximum of the energy distribution is shifted from EK1.7E_{K}\approx 1.7 eV in collisions with 11.3711.37 MeV/u S14+ to EK0.7E_{K}\approx 0.7 eV in collisions with 77 GeV/u U92+ where the low-energy part of the spectrum is strongly enhanced (see figure 7a) and this enhancement becomes even much more pronounced if the fragmentation events with large angles Θ𝑲\Theta_{\bm{K}} are selected (figure 7b).

In collisions with 11.3711.37 MeV/u S14+ and 77 GeV/u U92+ the effective strength of the projectile field, given by the ratio Zp/vZ_{p}/v, is essentially the same (0.680.68 and 0.660.66, respectively). Therefore, the strong enhancement of the lower-energy part of the fragmentation spectrum, which reflects the corresponding strong increase in the breakup of dimers with very large instantaneous size RR, is caused by a much higher impact energy. One should note, however, that since He2 is a very light target, a large increase in the collision velocity has a more profound overall effect on the spectrum shape and the total number of events than an increase in the Lorentz factor γ\gamma.

III.3 The energy spectrum shape versus the He2 binding energy

Refer to caption
Figure 8: The cross section dσfrDFdEK\frac{d\sigma^{\text{DF}}_{\text{fr}}}{dE_{K}}, versus the binding energy IbI_{b} of the He2 in collisions with 77 GeV/u U92+ projectiles. Solid, dashed and dotted curves correspond to Ib=139I_{b}=139 neV, Ib=148I_{b}=148 neV and Ib=130I_{b}=130 neV, respectively. The results for Ib=148I_{b}=148 neV and Ib=130I_{b}=130 neV were normalized to those for Ib=139I_{b}=139 neV at EK=60E_{K}=60 meV.

The values for the binding energy IbI_{b} of the He2 dimer, reported in the literature, vary between Ib=44.8I_{b}=44.8 neV 44.8 and Ib=161.7I_{b}=161.7 neV 161.7 . The distribution of the probability density ρ(R)=|Ψi(R)|2\rho(R)=|\Psi_{i}(R)|^{2} in the ground state of the dimer depends on the value of IbI_{b} and a variation ΔIb\Delta I_{b} in the binding energy would be reflected by the corresponding variation Δρ(R)\Delta\rho(R) of the probability density ρ(R)\rho(R). However, if ΔIb\Delta I_{b} is relatively small, ΔIbIb\Delta I_{b}\ll I_{b}, a large range of RR has to be spanned in order for the changes in the shape of ρ(R)\rho(R) to become noticeable.

In this respect, as it was already emphasized in He2-rel-HCI-lett , the fragmentation by ultra-fast projectiles can be of especial interest. Indeed, since the adiabatic collision radius RaR_{a} increases with the impact energy as vγ\sim v\gamma, they possess a very large effective interaction range that is a great advantage in probing the structure of such enormous objects like the He2 dimer.

In figure 8 we display the energy spectrum dσfrDFdEK\frac{d\sigma^{\text{DF}}_{\text{fr}}}{dE_{K}} for the He2 fragmentation by 77 GeV/u U92+. In the range of EKE_{K}, shown in the figure, the direct fragmentation is by far the dominant breakup channel and the cross section dσfrDFdEK\frac{d\sigma^{\text{DF}}_{\text{fr}}}{dE_{K}} represents the total energy spectrum.

The spectrum was calculated for three values of the dimer binding energy IbI_{b}: 130130 neV, 139139 neV and 148148 neV. For a better visibility of the variation of the spectrum shape with IbI_{b}, the spectra for Ib=130I_{b}=130 neV and Ib=148I_{b}=148 neV were normalized to the spectrum for Ib=139I_{b}=139 neV at an energy EK=60E_{K}=60 meV.

In the energy interval, 6060 meV EK1\leq E_{K}\leq 1 eV, the ratios dσfrDF(Ib=148)dEK/dσfrDF(Ib=139)dEK\frac{d\sigma^{\text{DF}}_{\text{fr}}(I_{b}=148)}{dE_{K}}/\frac{d\sigma^{\text{DF}}_{\text{fr}}(I_{b}=139)}{dE_{K}}, dσfrDF(Ib=139)dEK/dσfrDF(Ib=130)dEK\frac{d\sigma^{\text{DF}}_{\text{fr}}(I_{b}=139)}{dE_{K}}/\frac{d\sigma^{\text{DF}}_{\text{fr}}(I_{b}=130)}{dE_{K}} and dσfrDF(Ib=148)dEK/dσfrDF(Ib=130)dEK\frac{d\sigma^{\text{DF}}_{\text{fr}}(I_{b}=148)}{dE_{K}}/\frac{d\sigma^{\text{DF}}_{\text{fr}}(I_{b}=130)}{dE_{K}} vary by about 18%18\%, 20%20\% and 42%42\%, respectively. The corresponding relative variations, ΔIb/Ib\Delta I_{b}/I_{b}, of the binding energy are approximately 6.5%6.5\%, 7%7\% and 14%14\%, respectively. This means that in the energy interval under consideration the shape of the energy spectrum ’magnifies’ the variation of the binding energy by about three times.

A more detailed analysis of the cross sections shows that the range of quite small EKE_{K} (EK0.2E_{K}\lesssim 0.2 eV) gives the main contribution to the variation of the cross section ratios. This is not very surprising since the small energy range,   6060 meV EK0.2\leq E_{K}\lesssim 0.2 eV,   corresponds to a very large interval, 136136 a.u. R454\lesssim R\lesssim 454 a.u., of the interatomic distances in the He2 dimer, where the shape of its ground state becomes sensitive to even a very modest variation of the binding energy. The projectile is able to span in the collision so large interval of distances since at 77 GeV/u the adiabatic collision radius Ra103R_{a}\simeq 10^{3} a.u. already significantly exceeds the size of the dimer.

III.4 The total fragmentation cross sections

According to our r-SEA calculations, the contribution σfrDF\sigma^{\text{DF}}_{\text{fr}} of the DF mechanism to the total cross section for the breakup of the He2 dimer in collisions with 11 and 77 GeV/u U92+ projectiles amounts to 3.653.65 Mb and 2.392.39 Mb, respectively.

The relativistic effects in the direct fragmentation are most significant when the ’instantaneous’ transverse size RR_{\perp} of the He2 dimer in the collision is very large. Consequently, they most profoundly influence the fragmentation events with small EKE_{K} and large Θ𝑲\Theta_{\bm{K}}. However, such events have relatively small probabilities and, as a result, the cross section σfrDF\sigma^{\text{DF}}_{\text{fr}} is very weakly influenced by the relativistic effects. In particular, for collisions with 11 and 77 GeV/u U92+ the difference between our results for σfrDF\sigma^{\text{DF}}_{\text{fr}}, obtained using relativistic and nonrelativistic approaches, is about merely 1%\simeq 1\%.

The field of the projectile acting on the dimer is strongest when the impact parameters with respect to both atomic sites of the dimer are small, i.e. in collisions where the ’instantaneous’ transverse size RR_{\perp} of the dimer is not large. Since such collisions give a more significant contribution to σfrDF\sigma^{\text{DF}}_{\text{fr}}, than those with very large RR_{\perp}, the influence of the higher-order effects in the interaction between the projectile and the dimer on the magnitude of σfrDF\sigma^{\text{DF}}_{\text{fr}} is substantially larger reaching about 14%14\% for the fragmentation by 11 and 77 GeV/u U92+.

In the present paper the He2 fragmentation via the IE-ICD and DI-RET mechanisms was not calculated. Nevertheless, we can still obtain a rough estimate for the contribution to the total fragmentation cross section due to these two mechanisms. According to the experimental results of thesis_He2_S14+ , in the breakup of the He2 dimer by 11.3711.37 MeV S14+ the ratio of the summed contributions of the IE-ICD and DI-RET channels to the contribution of the DF channel is approximately equal to 1.61.6. Taking into account that our calculation for the DF contribution in collisions with 11.3711.37 MeV/u S14+ yields 1.81.8 Mb, we obtain that the summed contributions of the IE-ICD and DI-RET to the fragmentation by 11.3711.37 MeV S14+ is 1.8×1.6=2.9\approx 1.8\times 1.6=2.9 Mb.

The cross sections for the fragmentation via the IE-ICD and DI-RET are proportional to the cross sections for simultaneous ionization-excitation and double ionization, respectively, of the helium atom. Both these cross section scale similarly with ZpZ_{p} and vv. Therefore, by calculating them for collisions with 11.3711.37 MeV/u S14+ and 11 and 77 GeV/u U92+ projectiles and using our results for the DF cross sections, we have estimated that the summed contribution of the IE-ICD and DI-RET mechanisms to the fragmentation in collisions with 11 and 77 GeV/u U92+ is roughly equal to 4.34.3 Mb and 33 Mb, respectively.

IV Conclusions

In conclusion, we have studied the fragmentation of the helium dimer into singly charged helium ions by relativistic highly charged projectiles in collisions with relatively low kinetic energy release EK3E_{K}\lesssim 344 eV. Such breakup events occur solely due to the direct fragmentation mechanism in which the projectile ionized both dimer’s atoms in a single collision. In this mechanism the two helium ions are produced within the so short time interval (1016\lesssim 10^{-16} s) that during the production the helium nuclei remain essentially at rest. Consequently, the energy EKE_{K} of the ionic fragments, produced via the direct mechanism, is very simply (EK=1/RE_{K}=1/R) related to the size RR of the He2 dimer at the collision instant.

We have investigated in detail the energy and angular spectra of the He+ ions produced in collisions with 11 and 77 GeV/u U92+ projectiles. Our main findings can be summarized as follows.

In collisions with γ1\gamma\gg 1 the fragmentation events, in which the He+ ions move with low kinetic energies (EK0.1E_{K}\lesssim 0.1 eV) under large angles (20Θ𝑲16020^{\circ}\lesssim\Theta_{\bm{K}}\lesssim 160^{\circ}) with respect to the projectile velocity, are strongly affected by the relativistic effects.

The shape of the energy spectrum of the He+ ions is quite sensitive to the binding energy of the He2 dimer which can be exploited for its precise determination. Here the relativistic effects also play the important role since, by significantly enhancing the lower-energy part of the spectrum, they enable one to span a substantially broader range of the inter-atomic distances RR in the dimer, effectively increasing the sensitivity of the spectrum shape to the variation of the dimer binding energy.

The relativistic effects, having a strong impact on the fragments with low EKE_{K} and large Θ𝑲\Theta_{\bm{K}} and making the energy spectrum shape more sensitive to the variation in the dimer binding energy, influence nevertheless very weakly the total amount of the fragmentation events caused by the direct mechanism because they affect only their minor part. In this respect, a large increase in the collision velocity has a much stronger overall effect on the spectrum shape and the total number of events.

In contrast to relativistic effects, the role of the higher-order effects in the projectile-dimer interaction rises with increasing energy of the He+ fragments becoming substantial at EK1E_{K}\gtrsim 1 eV. Since such events give a more significant contribution to the fragmentation than those with very small EKE_{K}, the influence of the higher-order effects is substantially larger reaching about 14%14\% for the total cross section in the direct fragmentation by 11 and 77 GeV/u U92+.

According to our calculations the contribution of the direct mechanism to the total cross section for the He2 fragmentation by 11 and 77 GeV/u U92+ amounts to 3.653.65 Mb and 2.392.39 Mb, respectively. A rough estimate for the total fragmentation cross section in these collisions, which takes into account all the fragmentation mechanisms, suggests that it is about two times larger than the above values.

Acknowledgements.
BN, SFZ and XM gratefully acknowledge the support from the ‘National Key Research and Development Program of China’ (Grant No. 2017YFA0402300) and the CAS President’s Fellowship Initiative. Our numerical results were obtained using the facilities of the Supercomputer Center HIRFL at the Institute of Modern Physics (Chinese Academy of Sciences).

References

  • (1) J. Titze, M. S. Schöffler, H.-K. Kim, et al, Phys. Rev. Lett. 106, 033201 (2011).
  • (2) H.-K. Kim, H. Gassert, J. N. Titze, et al, Phys. Rev. A 89 022704 (2014);
  • (3) H.-K. Kim, PhD Thesis, (Frankfurt University, 2014) (unpublished).
  • (4) A. Jacob, C. Müller and A: B. Voitkiv, Phys. Rev. A 103, 042804 (2021).
  • (5) B. Najjari, Z. Wang and A. B. Voitkiv, Phys. Rev. Lett. 127, 203401 (2021).
  • (6) R. E. Grisenti, W. Schöllkopf, J. P. Toennies, G. C. Hegerfeldt, T. Köhler, and M. Stoll, Phys Rev Lett 85, 2284 (2000).
  • (7) B. A. Friedrich, Physics 6, 42 (2013).
  • (8) H. B. G. Casimir and D. Polder, Phys. Rev. 73, 360 (1948).
  • (9) F. Luo, G. Kim, G. C. McBane, C. F. Giese, and W. R. Gentry, J. Chem. Phys. 98, 9687 (1993); M. Jeziorska, W. Cencek, K. Patkowski, B. Jeziorski, and K. Szalewicz, J. Chem. Phys. 127, 124303 (2007).
  • (10) T. Havermeier, T. Jahnke, K. Kreidi, et al, Phys. Rev. Lett. 104, 153401 (2010).
  • (11) T. Havermeier, T. Jahnke, K. Kreidi, et al, Phys. Rev. Lett. 104, 133401 (2010).
  • (12) N. Sisourat , N. V. Kryzhevoi1 , P. Kolorenc, S. Scheit, T. Jahnke and L. S. Cederbaum, Nature Physics 6, 508 (2010).
  • (13) S. Zeller, M. Kunitski, J. Voigtsberger, et al, www.pnas.org/cgi/doi/10.1073/pnas.1610688113 (2016).
  • (14) For example, in collisions of 11 GeV/u U92+ with He the radiative electron capture cross section σrec1023\sigma_{rec}\simeq 10^{-23} cm2 (see e.g. results for radiative recombination on p. 272 of eic ) is about 88 orders of magnitude smaller than the cross section for single ionization of He σi1015\sigma_{i}\approx 10^{-15} cm2 (see e.g. He-ioniz ). The process of nonradiative electron capture is characterized by even much weaker cross sections.
  • (15) J. Eichler and W. E. Meyerhof, Relativistic Atomic Collisions (Academic Press, San Diego, 1995).
  • (16) A. B. Voitkiv and A. V. Koval, J. Phys. B 31, 499 (1998).
  • (17) A. B. Voitkiv, B. Najjari, R. Moshammer, and J. Ullrich, Phys. Rev. A 65, 032707 (2002).
  • (18) B. Najjari and A. B. Voitkiv, Phys. Rev. A 104 033104 (2021).
  • (19) J. Matthew and Y. Komninos, Surf. Sci. 53, 716 (1975);
  • (20) L. S. Cederbaum, J. Zobeley, and F. Tarantelli, Phys. Rev. Lett. 79, 4778 (1997); R. Santra, J. Zobeley, L. S. Cederbaum, and N. Moiseyev, Phys. Rev. Lett. 85, 4490 (2000).
  • (21) This mechanism was initially considered in He2_S14+ for collisions of 11.3711.37 MeV/u S14+ with He2 dimers having the ’instantaneous’ internuclear distance RR in the range 55 a.u. R10\leq R\leq 10 a.u.
  • (22) For the He2 fragmentation by photo absorption a similar mechanism was explored experimentally and theoretically in He2_photon and we-21-reflect-appr , respectively.
  • (23) In He2-rel-HCI-lett the probabilities PA+(𝒃)P_{A^{+}}(\bm{b}) and PB+(𝒃)P_{B^{+}}(\bm{b}^{\prime}) were defined as PA+(𝒃)=2w(𝒃)P_{A^{+}}(\bm{b})=2\,w(\bm{b}) and PB+(𝒃)=2w(𝒃)P_{B^{+}}(\bm{b}^{\prime})=2\,w(\bm{b}^{\prime}). The inclusion of the factors (1w(𝒃))(1-w(\bm{b})) and (1w(𝒃))(1-w(\bm{b}^{\prime})) somewhat reduces the cross section values compared to those reported in He2-rel-HCI-lett .
  • (24) F. Martin and A. Salin, Phys. Rev. A 55, 2004 (1997).
  • (25) M. Abramowitz and I. Stegun, Handbook of Mathematical Functions (New York: Dover, 1964).
  • (26) For instance, in the interval 5050 a.u. R500\lesssim R_{\perp}\lesssim 500 a.u. the value of CC varies by just about 16%16\% (between 0.560.56 and 0.650.65 if RR_{\perp}, RaR_{a}, ZpZ_{p} and vv are given in a.u. while the cross section is in Mb).
  • (27) R. Moshammer et al, Phys. Rev. Lett. 79, 3621 (1997).
  • (28) A. B. Voitkiv and B. Najjari, J. Phys. B 37, 4831 (2004).
  • (29) M. Przybytek et al, Phys. Rev. Lett. 104, 183003 (2010).
  • (30) P. D. Fainstein and R. D. Rivarola, J. Phys. B 20, 1285 (1987).
  • (31) A. B. Voitkiv, B. Najjari and J. Ullrich, J. Phys. B 38, L107 (2005).
  • (32) R. Feltgen, H. Kirst, K. A. Köhler et al, J. Chem. Phys. 76(5), 2360 (1982).
  • (33) A. R. Janzen, R.A. Aziz, J. Chem. Phys. 107(3), 914 (1997).