Fingerprints of slingshot non-sequential double ionization on two-electron probability distributions
Georgios Petros Katsoulis, Agapi Emmanouilidou

TL;DR
This paper investigates the mechanisms of non-sequential double ionization in helium under specific laser conditions, identifying distinctive two-electron distribution signatures of the slingshot-NSDI pathway.
Contribution
It reveals the properties of two-electron escape in slingshot-NSDI and identifies observable signatures in electron energy distributions.
Findings
Slingshot-NSDI dominates at low intensities.
Electrons escape with large, roughly equal energies in slingshot-NSDI.
Distinct two-electron probability distribution features characterize the pathway.
Abstract
We study double ionization of He driven by a near-single-cycle laser pulse at low intensities at 400 nm. Using a three-dimensional semiclassical model, we identify the pathways that prevail non-sequential double ionization (NSDI). We focus mostly on the delayed pathway, where one electron ionizes with a time-delay after recollision. We have recently shown that the mechanism that prevails the delayed pathway depends on intensity. For low intensities slingshot-NSDI is the mechanism that prevails. Here, we identify the differences in two-electron probability distributions of the prevailing NSDI pathways. This allows us to identify properties of the two-electron escape and thus gain significant insight into slingshot-NSDI. Interestingly, we find that an observable fingerprint of slingshot-NSDI is the two electrons escaping with large and roughly equal energies.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Fingerprints of slingshot non-sequential double ionization on two-electron probability distributions
G. P. Katsoulis
Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
A. Emmanouilidou
Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
Abstract
We study double ionization of He driven by a near-single-cycle laser pulse at low intensities at 400 nm. Using a three-dimensional semiclassical model, we identify the pathways that prevail non-sequential double ionization (NSDI). We focus mostly on the delayed pathway, where one electron ionizes with a time-delay after recollision. We have recently shown that the mechanism that prevails the delayed pathway depends on intensity. For low intensities slingshot-NSDI is the mechanism that prevails. Here, we identify the differences in two-electron probability distributions of the prevailing NSDI pathways. This allows us to identify properties of the two-electron escape and thus gain significant insight into slingshot-NSDI. Interestingly, we find that an observable fingerprint of slingshot-NSDI is the two electrons escaping with large and roughly equal energies.
keywords:
Keyword1, Keyword2, Keyword3
Introduction
Non-sequential double ionization (NSDI) of atoms driven by intense laser fields is a fundamental process which has attracted considerable theoretical and experimental interest [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 18, 17, 19, 20, 21, 22, 23]. The three-step model underlies NSDI [1]. First, one electron tunnels through the field-lowered Coulomb potential. This tunnel-ionizing electron can return to the core and transfer energy to the other electron. For high intensities, the main pathway of NSDI is the direct one. The recolliding electron transfers sufficient energy to the other electron for both electrons to ionize shortly after recollision. For low intensities, the delayed pathway of NSDI prevails. The recolliding electron transfers enough energy to the bound electron for only one of the two electrons to ionize shortly after recollision. The other electron transitions to an exited state of the remaining ion.
Until recently, the delayed pathway of NSDI was generally accepted as being equivalent to recollision-induced excitation with subsequent field ionization (RESI) [7, 9]. According to RESI, the electron that transitions to an excited state after recollision, ionizes later in time at extrema of the laser field. It does so mainly with the assistance of the field. However, recently, we have shown that RESI does not always prevail the delayed pathway of NSDI. For He driven at 400 nm by near-single-cycle pulses, at intensities below the recollision threshold, we have shown that a new mechanism overtakes RESI in the delayed pathway of NSDI [24]. We labeled this mechanism slingshot-NSDI. We have shown that in slingshot-NSDI, following the transition to an excited state, the electron subsequently undergoes a Coulomb slingshot motion due to the attractive force of the nucleus. Coulomb slingshot is similar to the gravitational slingshot motion that alters the motion of a spacecraft around a planet. During Coulomb slingshot the electron ionizes mostly around the second extremum of the field with the assistance of both the nucleus and the laser field. Due to the electron that ionizes last undergoing slingshot motion, the two electrons escape in opposite directions along the laser field. Owning to this anti-correlated two-electron escape, slingshot-NSDI offers an alternative explanation to multiple recollisions. Multiple recollisions is a mechanism that was previously put forth in the context of RESI to explain anti-correlated two-electron escape [25, 26, 27, 28, 29]. This two-electron escape pattern has been found to prevail NSDI of several atoms driven by intense long duration pulses and has been the object of many theoretical and experimental studies [25, 26, 27, 28, 29, 30, 31, 32].
Here, we show that anti-correlated two-electron escape is not the only feature that distinguishes slingshot-NSDI from the other NSDI pathways. The main NSDI pathways, which are considered in the current study, are the direct pathway, RESI and the double delayed pathway. In the latter pathway both electrons ionize with a delay following recollision. We show that slingshot-NSDI has very distinct fingerprints in both energy and angular two-electron probability distributions. These features can be observed experimentally. Hence, this study paves the way for identifying slingshot-NSDI by kinematically complete experiments that employ carrier envelope phase (CEP)-controlled near-single-cycle pulses [29, 33, 34, 35].
Method
We consider He driven by a near-single-cycle laser pulse at intensities 51014 W/cm2 and 71014 W/cm2 at 400 nm. Both intensities are below the recollision threshold, which corresponds to an intensity of 8.61014 W/cm2. The latter corresponds to the maximum energy of the electron returning to the core being equal to the energy needed to transition to the first excited state of the remaining ion. The maximum energy of the electron returning to the core is 3.17 [1], which is equal to 23.7 eV at 51014 W/cm2 and 33.1 eV at 71014 W/cm2; and are the strength and frequency of the field.
We use a laser field of the form
[TABLE]
where is the CEP, and fs is the full-width-half-maximum of the pulse duration in intensity. We employ atomic units, unless otherwise stated.
We use a three-dimensional (3D) semiclassical model that is formulated in the framework of the dipole approximation [36]. Previous successes of this model include verifying that electron backscattering from the nucleus accounts for the finger-like structure in NSDI of He driven by long laser pulses at higher intensities [36]. This finger-like structure was predicted theoretically [37] and obtained experimentally [17, 18]. Moreover, it was explained in a classical framework [36, 38]. In addition, using this 3D model, we investigated the direct versus the delayed pathway of NSDI for He driven by a long duration laser pulse at 400 nm [39]. For intensities ranging from below- to above-the-recollision threshold, we achieved excellent agreement with fully ab-initio quantum mechanical calculations. In addition, using this model we obtained very good agreement with experimental results for several observables of NSDI for Ar when driven by near-single-cycle laser pulses at 800 nm [40]. These observables were obtained as a function of CEP for intensities ranging from below- to above-the-recollision threshold.
In this model, one electron (recolliding) tunnel-ionizes through the field-lowered Coulomb-barrier with a tunnel-ionization rate that is described by the quantum mechanical Ammosov-Delone-Krainov (ADK) formula [41, 42]. We select the tunnel-ionization time, t0, using importance sampling [43] in the time interval the field is present, that is, [-2,2]. The importance sampling distribution is given by the ADK ionization rate. The exit point of the recolliding electron is along the laser-field direction and is computed using parabolic coordinates [44]. The electron momentum is taken to be equal to zero along the laser field while the transverse momentum is given by a Gaussian distribution [41, 42]. The initially bound electron is described by a microcanonical distribution [45]. The weight of each classical trajectory i that we propagate in time is given by
[TABLE]
where
[TABLE]
is the ADK ionization rate[41, 42] at the time of tunnel-ionization. The effective principal quantum number, , is given by , while and is the first ionization potential. The weight for electron 1 to have a transverse velocity equal to at the time is denoted by and is given by
[TABLE]
Once the initial conditions are specified at time , the position and momentum of each electron are propagated classically in time. We do so using the three-body Hamiltonian of the two electrons with the nucleus kept fixed. All Coulomb forces are accounted for: the interaction of each electron with the nucleus and the laser field and the electron-electron interaction are all included in the time propagation. We also account for the Coulomb singularity by using regularized coordinates [46]. During the time propagation each electron is interacting with the nucleus with charge . A trajectory is labeled as a doubly-ionized event if asymptotically, i.e. , the energies of both electrons are positive. The double ionization probability is given by
[TABLE]
where and N are the numbers of doubly-ionized and all events, respectively.
We identify the main pathways of energy transfer in each double ionization event. To do so we compute the time difference between the recollision time and the ionization time ti of each electron, with . We label the electron that ionizes first as electron 1 and the one that ionizes last as electron 2. To identify these times, for each classical trajectory, we compute the pair potential energy and obtain its maximum which corresponds to the time of minimum approach of the two electrons. We define this time of minimum approach as the recollision time. Moreover, we define the ionization time for each electron, , as the time when the compensated energy becomes positive and remains positive thereafter [47], with and ; is the vector potential and . We compare the time difference between the recollision time and the ionization time of each electron with the time interval tdiff where the electron pair potential energy undergoes a sharp change due to recollision. For the laser field intensities considered in this work, we find tdiff to be roughly equal to 1/8 laser cycle (T). We list in Table 1 the conditions satisfied by the direct, delayed or double delayed double ionization events. We find that in the delayed pathway, the probability for electron 2 to be the recolliding electron increases with decreasing intensity.
For the results presented in this work, we consider the intensities 51014 W/cm2 and 71014 W/cm2. For both intensities, 12 CEPs are considered ranging from to in steps of 30*∘. For each , at 71014* W/cm2 we obtain roughly 104 doubly-ionized events as a result of running 500, 12-hour jobs, while at 51014 W/cm2 we obtain 5103 doubly-ionized events as a result of running 4000 jobs, 12-hour jobs; one job corresponds to 1 CPU. For the results presented regarding total double ionization the average has been taken over all CEPs for each intensity.
Results
Slingshot-NSDI
First, we identify the prevailing pathways of double ionization when He is driven by a 2 fs laser pulse at 400 nm and at intensities of 51014 W/cm2 and 71014 W/cm2. We find that the pathways prevailing NSDI ionization are the delayed and the direct ones at 71014 W/cm2 and the delayed and the double delayed ones at 51014 W/cm2. For the delayed pathway, we note a significant change that takes place with decreasing intensity. The mechanism that prevails switches from RESI to slingshot-NSDI. Since both RESI and slingshot-NSDI are mechanisms of the delayed pathway, in both mechanisms electron 1 ionizes soon after recollision and electron 2 transitions to an excited state of He*+*. Slingshot-NSDI and RESI differ in the mechanism that underlies the subsequent ionization of electron 2. In RESI this mechanism involves electron 2 ionizing with the assistance of the laser field at later times at extrema of the field. In contrast, in slingshot-NSDI we find that electron 2 ionizes with the assistance of both the nucleus and the laser field. Next, we briefly outline some of the properties of slingshot-NSDI [24].
In particular, Fig. 1(a1) shows that in slingshot-NSDI electron 2 undergoes a Coulomb slingshot motion around the nucleus, see motion enclosed by the black arrows in Fig. 1(a1). In Fig. 1, we consider a double ionization event corresponding to CEP equal to zero but similar results hold for other CEPs. As a result of the slingshot motion, electron 2 escapes opposite to electron 1 and its momentum, pz,2, undergoes a large change. The momentum of electron 2 is large and points along the direction of the force from the laser field at the start and at the end of the slingshot motion, see Fig. 1(a2). We show that this large change in the momentum of electron 2 is due to the nucleus, by expressing the total momentum, pz,2, as the sum of the momentum changes due to the interaction with the nucleus and electron 1, , and with the laser field, as follows:
[TABLE]
In the delayed pathway, the repulsive force between the two electrons is roughly zero shortly after recollision, thus, contributing only a constant term to . We plot the momentum change due to the nucleus of electron 2 in Fig. 1(a2). It is clear that, during the Coulomb slingshot motion, the sharp change of the total momentum of electron 2 is mainly due to the term .
Next, we show that, during the slingshot motion, the laser field provides sufficient energy to electron 2 to ionize due to the large change in the electron momentum. Shortly after recollision, at time t, the repulsive force between the two electrons is significantly smaller than during recollision. Hence, after this time, we can safely assume that the energy of electron 2 changes due to the work done mainly by the field as follows:
[TABLE]
where F is the force from the laser field and is the rate of change of the energy of electron 2. During Coulomb slingshot, the close encounter of electron 2 with the nucleus at t takes place past a zero of the laser field. At this time both the momentum of electron 2 and the force from the laser field point along the +-axis. Roughly half a laser cycle later, in the time interval [0.75, 1.25]T, the slingshot motion is concluded with both the momentum of electron 2 and the force of the laser field pointing along the --axis. Thus, during the Coulomb slingshot, is mostly positive in the first half cycle [0.25, 0.75]T and the second one [0.75, 1.25]T after recollision, see red-shaded area in Fig. 1(a3). Fig. 1(a3) clearly shows that is positive due to being positive, see blue-shaded area in Fig. 1(a3). Thus, the large change in the momentum of electron 2 due to Coulomb slingshot is the reason the rate of change of the energy provided to electron 2 by the laser field is positive. This leads to ionization of electron 2 around the second extremum of the field after recollision, i.e. in the time interval [0.75, 1.25]T. In contrast, in RESI electron 2 can ionize at any extremum of the laser field.
As we increase the laser intensity, it is more probable that the force exerted on electron 2 from the laser field is larger than the attractive Coulomb force exerted from the nucleus. As a result, in most delayed double ionization events, electron 2 does not return at time close to the nucleus and therefore does not undergo Coulomb slingshot. That is, following the transition of electron 2 to an excited state, the effect of the nucleus on the motion of this electron decreases with increasing intensity. This results in RESI becoming more important than slingshot-NSDI with increasing intensity.
Two-electron momentum probability distributions
In Fig. 2, we plot the correlated electron momenta. That is, we plot the double ionization probability as a function of each electron momentum along the direction of the laser field at intensities 51014 W/cm2 (top row) and 71014 W/cm2 (bottom row). As the intensity decreases, the two-electron escape changes from correlated, with both electrons escaping in the same direction along the laser field (Fig. 2(b1)), to anti-correlated, see Fig. 2(a1). Two are the main reasons for the transition from anti-correlated two-electron escape at 51014 W/cm2 to correlated one at 71014 W/cm2. At both intensities we find that two pathways of NSDI prevail, with the delayed pathway being one of the two most important pathways in both cases. However, the other prevailing pathway is the direct one at the higher intensity while it is the double delayed one at the lower intensity. In the direct pathway, both electrons ionize soon after recollision around a zero of the laser field. As a result, the two electrons escape with large final electron momenta, roughly given by -, in the same direction along the laser field, see Fig. 2(b2). In the double delayed pathway, the two electrons ionize with a delay after recollision resulting in smaller momenta along the laser field compared to direct events. In the double delayed events, the electrons ionize mostly in opposite directions, see Fig. 2(a2).
In addition, as the intensity decreases from 71014 W/cm2 to 51014 W/cm2, the mechanism underlying the delayed pathway changes from RESI to slingshot-NSDI. At 71014 W/cm2, once electron 2 transitions to an excited state of the remaining ion, the electron subsequently ionizes mostly due to the laser field around extrema of the field. As a result electron 2 escapes with a small final momentum. The electron that ionizes first escapes with large momentum, roughly equal to -. Indeed, this is the pattern exhibited by the correlated electron momenta for RESI both at 51014 W/cm2 and at 71014 W/cm2, see Fig. 2(a3) and (b3). At 51014 W/cm2, through Coulomb slingshot, the nucleus plays a significant role in ionizing electron 2 after it transitions to an excited state. This effect results in both electrons escaping with large momenta in opposite directions along the laser field, see the two-electron probability distribution enclosed by the squares in Fig. 2(a4). The role of the nucleus diminishes with increasing intensity even for slingshot-NSDI events. Indeed, at the higher intensity, for slingshot-NSDI events, electron 2 does not have quite as large momentum as the first to ionize electron, see Fig. 2(b4) and compare with Fig. 2(a4).
In Fig. 3, we plot the double ionization probability as a function of the difference of the electron momenta along the laser field and of the projection of the perpendicular momentum of one electron along the direction of the perpendicular momentum of the other electron, denoted by in Fig. 3(a1). At 71014 W/cm2, Fig. 3(b1) shows that for NSDI events with similar electron momenta along the laser field, i.e. , electron-electron repulsion results in the two electrons escaping with opposite momenta in the direction perpendicular to the laser field. This pattern is due to the direct pathway of NSDI, as Fig. 3(b2) clearly shows. In contrast, electron-electron repulsion does not significantly contribute to RESI and slingshot-NSDI. This is clearly seen in Fig. 3(b3) (RESI) and (b4) (slingshot-NSDI) at 71014 W/cm2 and in Fig. 3(a3) (RESI) and (a4) (slingshot-NSDI) at 51014 W/cm2. At the lower intensity, a comparison of Fig. 3 (a2), (a3) and (a4) shows that electron-electron correlation plays a more important role for the double delayed events rather than for the RESI and the slingshot-NSDI events. This is consistent with both electrons ionizing with a delay in the double delayed events.The authors in ref. [23], use a similar two-electron distribution for all double ionization events for He driven by long laser pulses at small intensities at 400 nm. They find that electron-electron repulsion is smaller for anti-correlated two-electron escape rather than for a correlated one. This is consistent with our finding that electron-electron correlation is larger for direct events compared to RESI and slingshot-NSDI events.
Two-electron angular probability distributions
In Fig. 4, we plot the double ionization probability as a function of the inter-electronic angle between the momenta of the two escaping electrons and the angle formed by the momentum of one of the two electrons with respect to the axis of the laser field. Comparing Fig. 4(a1) and (b1), we find that for most double ionization events the inter-electronic angle is close to 180*∘* at 51014 W/cm2 and close to 0∘ at 71014 W/cm2. This is consistent with anti-correlated two-electron escape prevailing at 51014 W/cm2 and correlated prevailing at 71014 W/cm2, as we have already seen in Fig. 2(a1) and (b1), respectively. At 51014 W/cm2, Fig. 4(a1) shows that for most double ionization events, in addition to the two electrons escaping opposite to each other, one of the two electrons escapes along the polarization direction, i.e. is 0*∘* or 180*∘. This latter pattern of two-electron escape is more pronounced for slingshot-NSDI, see Fig. 4(a4) and compare with Fig. 4(a2) and (a3). We find that the electron escaping along the polarization axis is the one that ionizes first after recollision, giving rise to the probability distribution enclosed by the parallelograms in Fig. 4(a4). The remaining wedge-like shape in Fig. 4(a4) is accounted for by electron 2 forming an angle with the polarization axis while, for , electron 1 forms an angle 180∘* giving rise to and, for , electron 1 forms an angle 0*∘* giving rise to . At the higher intensity, the prevailing correlated two-electron escape is mostly due to the direct pathway of double ionization, see Fig. 4(b2) and compare with Fig. 4(b3) and (b4).
Two-electron energy probability distributions
Besides anti-correlated two-electron escape, we find that another hallmark of slingshot-NSDI is both electrons escaping with large energy. This is clearly seen at the smaller intensity of 51014 W/cm2 in Fig. 5(a4). In contrast, RESI gives rise to an unequal energy sharing between the two escaping electrons. This can be seen at 51014 W/cm2 in Fig. 5(a3) and at 71014 W/cm2 in Fig. 5(b3). Both in RESI and slingshot-NSDI electron 1 escapes with large energy. This is mainly due to the large momentum of electron 1 along the direction of the laser field at the recollision time, which is roughly equal with -. The main difference between RESI and slingshot-NSDI is the influence of the nucleus on electron 2 following its transition to an excited state after recollision. For RESI, the nucleus has a very small effect on electron 2. This electron ionizes with the help of the laser-field around field extrema, resulting in mostly small final energy of electron 2. As a result, unequal energy sharing prevails in RESI, see Fig. 5(a3) and (b3). In contrast, in slingshot-NSDI the nucleus plays a major role on electron 2. As we have previously discussed, this electron undergoes a Coulomb slingshot motion around the nucleus gaining a large amount of energy. As a result, roughly equal energy sharing prevails in slingshot-NSDI, see Fig. 5(a4). This pattern gives rise to a high concentration of NSDI probability around the diagonal at large energies in the two-electron energy probability distribution for all NSDI events, see Fig. 5(a1).
The effect of the nucleus on electron 2, following its transition to an excited state, decreases with increasing intensity. This is manifested in RESI overtaking slingshot-NSDI as the prevailing mechanism of the delayed pathway. This diminishing effect of the nucleus on electron 2 in the delayed pathway is also evident in slingshot-NSDI events. Indeed, the two electrons share the energy more unequally at 71014 W/cm2 in Fig. 5(b4) compared to 51014 W/cm2 in Fig. 5(a4). The reduced contribution of slingshot-NSDI to the delayed pathway as well as the reduced effect of the nucleus on electron 2 in slingshot-NSDI at 71014 W/cm2, accounts for the more unequal energy sharing of the two electrons for all NSDI events in Fig. 5(b1) versus more equal energy sharing at 51014 W/cm2 in Fig. 5(a1).
Conclusions
Using a 3D semiclassical model we investigate how two-electron probability distributions change for different NSDI pathways as well as for different intensities below-the-recollision-threshold. We do so for He driven by near-single-cycle laser pulses at 400 nm. This study allows us to gain significant insight into slingshot-NSDI, which is a new mechanism that prevails the delayed pathway of NSDI at low intensities, as we have recently shown [24]. We find that in slingshot-NSDI the electron that ionizes first does so along the polarization direction of the laser field. The electron that ionizes last escapes opposite to the other electron. This is due to electron 2 undergoing a Coulomb slingshot motion around the nucleus. Indeed, we have shown that, following recollision, unlike direct events of NSDI where electron-electron correlation has a significant effect, electron repulsion has a very small effect on both RESI and slingshot-NSDI. Moreover, we find that in slingshot-NSDI the two electrons escape both with large energy giving rise to a distinct equal energy pattern in the two-electron energy probability distribution for all double ionization events. Thus, in addition to the anti-correlated two-electron escape, which we have previously identified as a hallmark of slingshot-NSDI [24], two-electron escape with roughly equal energy sharing and large energies is yet another trademark of slingshot-NSDI which can be measured by future experiments.
Acknowledgements
A.E. acknowledges the EPSRC grant no. N031326 and the use of the computational resources of Legion at UCL.
Author contributions statement
G. P. Katsoulis performed the analysis of the computations involved and contributed to the ideas involved. A.E. provided the codes used for the computations and the analysis and was responsible for the main ideas involved in the theoretical analysis. All authors reviewed the manuscript.
Additional information
Competing financial interests: The authors declare no competing financial interests.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Corkum, P. B. Plasma perspective on strong field multiphoton ionization. Phys. Rev. Lett. 71 , 1994 (1993).
- 2[2] L’Huillier, A., Lompre L. A., Mainfray G., & Manus C. Multiply charged ions induced by multiphoton absorption in rare gases at 0.53 μ m 𝜇 m \mathrm{\mu m} . Phys. Rev. A 27 , 2503 (1983).
- 3[3] Lambropoulos, P. Mechanisms for Multiple Ionization of Atoms by Strong Pulsed Lasers. Phys. Rev. Lett. 55 , 2141 (1985).
- 4[4] Walker, B. et al. Precision Measurement of Strong Field Double Ionization of Helium. Phys. Rev. Lett. 73 , 1227 (1994).
- 5[5] Larochelle, S. F. J., Talebpour, A. & Chin, S. L. Non-sequential multiple ionization of rare gas atoms in a Ti:Sapphire laser field. J. Phys. B 31 , 1201 (1998).
- 6[6] Moshammer, R. et al. Momentum Distributions of Ne n + superscript Ne limit-from n \mathrm{Ne^{n+}} Ions Created by an Intense Ultrashort Laser Pulse. Phys. Rev. Lett. 84 , 447 (2000).
- 7[7] Kopold, R., Becker, W., Rottke, H. & Sandner, W. Routes to Nonsequential Double Ionization. Phys. Rev. Lett. 85 , 3781 (2000).
- 8[8] Weber, T. et al. Correlated electron emission in multiphoton double ionization. Nature (London) 405 , 658 (2000).
