Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Mapping initial transverse momenta of tunnel-ionized electrons to rescattering double ionization in nondipole regimes

Open Access Open Access

Abstract

We investigate the double ionization of a model Neon atom in strong middle infrared laser pulses by simulating the classical trajectories of the electron ensemble. After one electron tunnels out from the laser-dressed Coulomb barrier, it might undergo different returning trajectories depending on its initial transverse momentum, which in this wavelength may propagate along or deviate from the polarization direction. This initial transverse momentum determines the rescattering time, and thus some trajectories can have returning time longer than one optical cycle. These late-returning trajectories determine the correlated electron-electron momentum distribution for double ionization and allow us to disentangle each double ionization event from the final momentum distribution. The description of these trajectories allow us also to understand how the nondipole effects modify the correlated electron-electron momentum distribution in double ionization.

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Nonsequential double ionization (NSDI) driven by strong laser fields is one of the most fascinating and fundamental processes in ultrafast physics. It has attracted a lot of attentions both theoretically and experimentally in past decades [13]. Experimentally, many studies of NSDI were performed using Ti:Sapphire lasers in the near-infrared (near-IR) (800 nm) [410]. In this wavelength and in the tunneling regime, an electron is ionized at the field maximum where it can gain energy from the field and return to the parent ion [11]. During this revisit, the bound electron may be directly kicked off, if the incoming electron carries enough energy, or jump to an excited state to be ionized later by the field. For each of the two cases, the energy sharing between electrons is different and these two scenarios present distinct final correlated electron-electron momentum distributions (EMDs) [12].

The recent developments of middle infrared (mid-IR) laser technology are now able to produce few-cycle laser pulses with the wavelength longer than 3 micrometers [13]. With these pulses new strong field phenomena have been observed such as high-energy photons in high-harmonic generation [14,15], fork-like structures in photoelectron momentum distributions [16], and the partition of the photon momentum [17,18].

For the case of NSDI the observations show much stronger energy sharing asymmetry in the correlated EMD as compared to the near-IR wavelengths [19,20]. Another important difference with the near-IR case is the importance of the magnetic field induced by the laser field. Since the magnetic field may modify the rescattering and induce the asymmetry in the forward and backward propagation of the laser pulse, the high harmonic generation [21] as well as the high-order above threshold ionization [22] are modified. Dynamics such as the nondipole gating [23] or the unexpectedly high sum momenta in the propagation direction [24] are observed even when the laser intensity is weaker than the previously predicted criteria for the onset of magnetic-field effects, i.e. $\beta _{0}\approx U_{p}/(2\omega c)\approx 1$ [25,26], with $\beta =v/c$, $U_p=\frac {E_0^2}{4\omega ^2}$ the ponderomotive energy, where $\omega$, $E_0$ are the frequency and the electric field amplitude respectively, and $c$ is the light velocity.

Due to the complexity of multiple-electron systems, the fully dimensional time-dependent Schrödinger equation (TDSE), which affords complete and accurate information, is mainly limited to the double ionization induced by EUV and UV fields [2730]. The simulation of the double ionization driven by mid-IR laser fields is still a too heavy burden for the state-of-the-art computer. Therefore, reduced-dimensionality TDSE models and an extension of the time-dependent surface flux method are developed to qualitatively explain electron-electron correlation [3134]. Alternatively, the strong field approximation, which neglects the Coulomb action on final states, plays important roles for explaining the correlated EMDs [3539]. On the other hand, the purely classical method was developed by solving the Newton equation for the classical ensemble in combined laser and Coulomb fields when over-barrier ionization dominates [4044]. In order to include the tunneling process, a semi-classical method has been proposed, which describes the single ionization by the Ammosov-Delone-Krainov (ADK) theory [45] and treats the rest part purely classically, which reproduces and explains a lot of experimental phenomena, such as the knee structure in the plot of double-ionization yields versus laser intensity [4648], the double-hump structure of the longitudinal ion-momentum distribution [49], the finger-like structure in the correlated EMD [50,51], and the yield of anticorrelated electrons with high energy [52].

In this paper we use a semi-classical method to study the double ionization of Neon atoms and to describe the differences in the final EMD for the case of near-IR and mid-IR. We discuss the asymmetric energy sharing that is observed in the longer wavelengths and analyze three kinds of trajectories that contribute to different areas in the correlated EMD. Based on this understanding we continue to explain the nondipole effects over these trajectories. In this case, the modification of electron trajectories caused by the magnetic field or a slight ellipticity of the mid-IR laser pulse changes the relative proportion of different rescattering scenarios, which gives the explanation for the variation in the correlated EMD and the reduction of dication yields. The rest of the paper is organized as following. In Sec. II, we introduce the numerical method. The simulation results are shown in Sec. III. The paper ends with a conclusion in Sec. IV.

2. Simulation method

We simulate the double ionization of a Neon atom, in which two electrons are active and all other electrons freeze with the nucleus using a semi-classical approach. In this model, the exit point of the tunneling electron is determined by the effective potential in parabolic coordinates (atomic units are used throughout unless otherwise specified): $-\frac {1}{4\eta }-\frac {1}{8\eta ^{2}}-\frac {E(t)\eta }{8}=-\frac {I_{p1}}{4}$ with the first ionization energy of the atom $I_{p1}=0.8$ a.u. and $E(t)$ the instantaneous electric field. The $\eta$ is the parabolic coordinate relating to the Cartesian coordinate through $\eta =-2x$, with the positive $x$ direction assumed to be the instantaneous polarization direction, and the exit point is $(-\frac {\eta }{2},0,0)$ in Cartesian coordinates [53]. The corresponding Keldysh parameters $\gamma =\sqrt {I_{p1}/2U_{p}}$ are 0.21 (3000 nm) and 0.78 (800 nm) for single ionization of Ne atoms when the laser intensity is $3\times 10^{14}$ W/cm$^{2}$ [54]. The tunneling ionization rate is described by the ADK [45] formula $W(t)=\frac {4(2I_{p1})^2}{|E(t)|}\exp [-\frac {2(2I_{p1})^{3/2}}{3|E(t)|}]$. The longitudinal momentum of the tunneling electron is assumed to be zero, while the transverse momentum distribution is formulated as

$$\varpi(p_y,p_z,t)=\frac{2\sqrt{2I_{p1}}}{|E(t)|}\exp\left[-\frac{({p_z}^2+{p_y}^2)\sqrt{2I_{p1}}} {|E(t)|} \right].$$
The initial state of the bound electron is depicted by a microcanonical distribution by setting the ionization energy of Ne$^+$ to be $I_{p2}=$1.5 a.u.. More than $10^{9}$ trajectories are calculated in order to obtain converged simulation results. After the tunneling ionization, the subsequent evolution of the two electrons is governed by the Newton equation:
$$\frac{d^{2}\textbf{r}_i}{dt^{2}}={-}\nabla_{\textbf{r}_i}(V_{ne}^i+V_{ee})-\textbf{E}(t)-s\textbf{v}_{i}\times\textbf{B},$$
where $\textbf {v}_i$ and $\textbf {r}_i$ denote the velocity and displacement of the $i$th electron respectively. The Coulomb interactions between the particles are respectively denoted as $V_{ne}^i=-\frac {2}{\sqrt {\textbf {r}_{i}^{2}+a^{2}}}$ and $V_{ee}=\frac {1}{\sqrt {(\textbf {r}_{1}-\textbf {r}_{2})^2+b^{2}}}$. The above equation is solved by using the standard fourth- to fifth-order Runge-Kutta algorithm, where the soft-core parameters $a=0.1$ and $b=0.01$ were employed in the Coulomb interaction to avoid the hard-Coulomb singularities. The electric field is $\textbf {E}(t)=f(t)\frac {E_{0}}{\sqrt {1+\epsilon ^2}}[\hat {\textbf {x}}\cos (\omega t+\phi )-\hat {\textbf {y}}\epsilon \sin (\omega t+\phi )]$, where $f(t)$ denotes the pulse envelope with two cycles plateau followed by one cycle ramp off, the carrier envelope phase $\phi =0$, and $\epsilon$ is the ellipticity. The magnetic field connects to the electric field via the light speed $c$, i.e., $\textbf {B}(t)=\hat {\textbf {z}}\times \textbf {E}(t)$. The parameter $s$ is used to switch on or off the Lorentz force if $s=1$ or $s=0$, respectively. A double-ionization event is counted if the total energy of each electron, i.e., the addition of kinetic energy, electron-nuclei Coulomb potential energy, and half electron-electron Coulomb energy, is positive. Since the electron velocities are small compared to the speed of light in our simulation, the modification caused by the spatially inhomogeneity of the electric field is trivial, and we treat the electric field spatial homogeneity when we calculate the electron trajectory.

3. Simulation results and discussions

When the driving laser wavelength is 800 nm, the correlated EMD of double ionization mainly locates in the first and third quadrants, as shown in Fig. 1(a), and the finger-like structure described in the literature is well reproduced [6,7]. The correlated EMD almost does not change once the magnetic field is included in the simulation, as shown in Fig. 1(b). The correlated EMD driven by the 3000 nm laser pulse, as shown in Fig. 1(c), is quite different from the 800 nm case. The double ionization events locate close to the coordinate axes instead of the diagonal line, and some events invade into the second and fourth quadrants, which is well in agreement with the experimental observations [19,20]. The ponderomotive energy for the driving laser with 3000 nm wavelength ($U_p=9.25$ a.u.) is much higher than 800 nm ($U_p=0.66$ a.u.) when the laser intensities are same ($3\times 10^{14}$ W/cm$^{2}$). The energy sharing between the two electrons is more asymmetric, which gives the explanation that the finger-like structure is still visible but with a very large open angle [44]. Comparing Fig. 1(c) and Fig. 1(d), most electrons are in the first and third quadrants for both panels, however, the near-axis signals are more pronounced in Fig. 1(d). The distributions in each panel have been normalized by their own maxima.

 figure: Fig. 1.

Fig. 1. The correlated EMD when the laser wavelengths are 800 nm (a, b) and 3000 nm (c, d), and magnetic fields of the laser pulses are neglected (a, c) and included (b, d). The laser intensities are $3\times 10^{14}$ W/cm$^{2}$ for 800 nm and 3000 nm laser pulses. Both lasers are linearly polarized.

Download Full Size | PDF

We first investigate the correlated EMD for the case that the magnetic field is not included. Our analysis of the case of 3000 nm starts with Fig. 2. The upper panel shows the NSDI yields where the first electron is tunneling ionized around $t=0.5$T. The bars show the number of particles which rescatter in each of the three events leading to NSDI, the rescattering time windows are marked as shadow regions. The first rescattering is around $t=1.4$T and corresponds to the usual rescattering trajectory, the second is around $t=1.8$T and the third is around at $t=2.25$T, where the recollision time is defined as the instant when two electrons are closest after the tunneling of the first electron. The short trajectories (rescattering times smaller than 1.2T) are not presented for their low proportion (less than 2.8$\%$). One may see that the multipe rescattering has larger probabilities than the first rescattering, which is general when the driving laser pulse has very long wavelengths [44]. The significance of the late returns has been studied by Shvetsov-Shilovski $et$ $al$. using elliptically polarized laser pulses [55], and later observed by Kang $et$ $al$. in the experiment [56]. In Fig. 2(b), the rescattering energy as a function of the rescattering time is presented. The rescattering time and energy (green solid, dashed and dash-dotted lines in Fig. 2(b)) are obtained by solving the Newtonian equation $\frac {d^{2}\textbf {r}_i}{dt^{2}}=\frac {\partial A(t)}{\partial t}$ ,where $\textbf {r}_i$ denotes the electron displacement, $A$(t) is the laser vector potential, and the Coulomb and the magnetic force are neglected. The shapes of these three curves are quite similar , which indicates a precise map between the ionization time and the rescattering energy. However, the late rescattering events have distinct features compared to the first rescattering. Overlay in Fig. 2(b) is the incoming angle, which is defined as the cross angle between the moving direction of the incoming electron and the polarization axis at the moment $\Delta t$ = 3 a.u. just before recollision [57]. The angle is almost zero for the first rescattering, which means that the initial transverse momentum is also zero and the electrons travel along the laser polarization only. For the second and third rescattering events, the incoming angles are different from zero and change as a function of time. As the parallel momentum component is smaller at the edges of the rescattering time window, the incoming angle is larger. The colorscale in these curves represents the numbers of the rescattering electrons at that instant. From this analysis we observe that electrons collide with different times, angles and energies in these three rescattering time windows. This is an interesting extension to the well-known relation between ionization times and rescattering energy [58].

 figure: Fig. 2.

Fig. 2. (a) The laser electric field (black solid line) and the laser vector potential (amaranth dashed line). The electron tunnels out during the time interval [0.495T, 0.54T] will rescatter the parent ion in the time interval [1.25T, 1.45T] (the first rescattering), [1.55T, 1.9T] (the second rescattering) and [2.05T, 2,45T] (the third rescattering), and the corresponding double ionization probability are represented by the bars. (b) The incoming angle distribution for three rescattering events at different time. The variation in rescattering energy of these three events is depicted by green solid, dashed and dash-dotted lines. The laser wavelength is $\lambda =3000$ nm, and the laser intensity is $I=3\times 10^{14}$ W/cm$^2$.

Download Full Size | PDF

To understand the origin of the incoming angle, we analyze the trajectories related to these late rescattering events. The classical trajectory Monte Carlo simulations allow us to trace the electron trajectories and thus view the rescattering details.

Once the electron tunnels out the laser-dressed Coulomb potential, it travels in the combined laser-Coulomb fields. Figure 3 shows the electron trajectories and kinetic momenta along the $z$-axis and $x$-axis for electrons released around $t=0.5$T, and these trajectories will induce double ionization by rescattering. The three columns from left to right are for the events associated with the first, second, and third rescattering, respectively. The first and third rows show the time-dependent displacement $x(t)$ and $z(t)$. The corresponding momenta, $p_x(t)$ and $p_z(t)$, are shown in the second and fourth rows. One may clearly see from the first column that the electron travels along the laser polarization direction with a large excursion ($\sim 800$ a.u.) and rescatters the parent ion at $t=1.44$T. In the transverse direction, the electron trajectories stay close to zero during this time. The electron momentum $p_x(t)$ mirrors the laser vector potential, while $p_z(t)$ keeps constant before rescattering. Therefore, the first rescattering electrons travel along the laser polarization and their incoming angles are close to zero degree.

 figure: Fig. 3.

Fig. 3. The time-evolved electron displacement (the first row) and electron kinetic momentum (the second row) along the laser polarization direction. The third and fourth rows are the time-evolved electron displacement and electron kinetic momentum along the laser propagation direction. The three columns from left to right are for the first, second and third rescattering, respectively.

Download Full Size | PDF

For the second rescattering shown in the second column, the electron propagates in the laser field and backs to $x=0$ at $t=1.26$T. However, the rescattering does not happen at this moment because the electron is away from the ion in the transverse direction. As the electron passes very quickly through $x=0$ around $t=1.26$T, the Coulomb force reverses the trajectory of the electron in the transverse direction. The sudden change in the transverse momentum can be seen in both the displacement $z(t)$ and $p_z(t)$. The transverse momentum changes signs due to the strong attraction of the Coulomb field. At a later time $t=1.9$T both transverse and parallel coordinates are close to zero ($x(t)=z(t)=0$) and the rescattering occurs.

The trajectories associated with the third rescattering are presented in the right column of Fig. 3. The electron misses the parent ion twice at $t$ = 1.44T and $t$ = 1.56T. Different from the scenario of the second rescattering, the electron is slow when it passes through $x=0$, and thus it may spend more time close to the parent ion. Thereby, the Coulomb field may act longer on this electron and reverse its momentum in the transverse direction even though it is up to $\sim 20$ a.u. away in the $z$ transverse direction. At a later time $t=2.44$T, both the $x(t)$ and $z(t)$ are close to zero and the third rescattering occurs almost two cycles later than the tunnel ionization event. Similar to the second rescattering, the synchronous approaching to zero in both $x$- and $z$-axes explains the rescattering time, energy and incoming angle.

Some typical trajectories are depicted in Fig. 4 in the three-dimensional space. The trajectories for tunnel-ionized electrons with the initial transverse momenta equal to zero ($p_z=p_y=0$), the first circle and the second circle in Fig. 5(a) produce rotationally symmetric (in the dipole case) surfaces corresponding to the first, second and third rescattering. On each surface, a typical trajectory determined by the Coulomb focusing and laser driving is plotted by a solid line. These trajectories also identified as glory trajectories [59] are the only trajectories that can rescatter in the combined effect of the laser field and the Coulomb potential. As pointed by Huang $et$ $al$., even higher-order rescattering may happen using a multiple-cycle laser pulse [44]. Here, the fourth rescattering is truncated because a few-cycle laser pulse is used. We may parenthetically point that the high-order returns are overestimated in this semi-classical model since the wave packet expansion is not comprised.

 figure: Fig. 4.

Fig. 4. The typical trajectories for the first (green), second (red) and third (blue) rescattering without considering the magnetic field. The laser wavelength is $3000$ nm, and the intensity is $3\times 10^{14}$ W/cm$^2$. Surfaces are formed by rotating each trajectories with respect to the $x$-axis. The dashed curve is the trajectory starting with the identical displacement and momentum of the first rescattering trajectory but including the magnetic effect.

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. (a) The initial transverse momentum distribution of the tunneling electron associated with the later double ionization. The magnetic field is not included. (b), (c) and (d): The correlated EMDs associated with the initial transverse momenta around 0, in the first and second circles, respectively. The laser wavelength is 3000 nm, and the intensity is $3\times 10^{14}$ W/cm$^{2}$. The laser is linearly polarized.

Download Full Size | PDF

The identification of the mapping between the initial transverse momentum to the rescattering time, energy and incoming angle allow us to disentangle all the rescattering events in the final correlated EMD. The three types of rescattering trajectories with well defined rescattering times can now help us to understand the formation of the correlated EMDs shown in Figs. 5(b), (c) and (d).

According to the Simpleman model [60], where the Coulomb potential and magnetic fields are neglected, the first rescattering occurs at the moment that the vector potential A is maximum, and the rescattering energy is up to 3.17$U_p$. The corresponding kinetic momentum $p=\sqrt {2\times 3.17U_p }=1.26A_0$ [61], where $A_0$ is the amplitude of the laser vector potential. For the mid-IR laser field, the rescattering energy is much larger than the Ne$^+$ ionization potential $I_{p2}=1.5$ a.u., which means the pitcher (the tunneling electron) only transfers very little energy to the catcher (the bound electron). Just after the rescattering, the pitcher and catcher move in the same direction. In the later propagation, the low energetic catcher is first reversed and then accelerated to high speed by the remaining laser field. However, the pitcher with the kinetic momentum up to 1.26$A_0$ is substantially decelerated but does not reverse its moving direction, and finally contribute to the double ionization events in the second and fourth quadrants. For the less energetic pitcher whose momentum just after the rescattering is less than A will finally reverse its direction, and the corresponding double ionization events locate in the first or third quadrants. According to the classical estimation, the critical laser vector potential for producing rescattering double ionization events in the second and fourth quadrants is $A=\sqrt {3.4I_{p2}}$. Differently, the second rescattering brings the rescattering energy less than $1.54U_p$ (the corresponding maximum kinetic momentum is 0.88$A_0$). Therefore, both the pitcher and catcher will be driven to the same direction by the laser field, as showed in Fig. 5(c). For the third rescattering, the maximum rescattering energy is 2.4$U_p$(the corresponding kinetic momentum is 1.1$A_0$), the pitcher moving direction is very likely to be reversed, thereby the double ionization events locating in the second and fourth quadrants are less than that of the first rescattering.

After having understood the rescattering details within the dipole approximation, we are ready to analyze how the magnetic field changes the rescattering processes and modifies the correlated EMD. When the magnetic field is take into consideration, the Lorentz force would affect the rescattering process of electrons and some interesting phenomena have been presented for double ionization [23,24]. The dashed line in Fig. 4 represents the trajectory of the electron starting with the identical initial position and momentum with those of the green solid line for the first rescattering, and now the electron is pushed to the positive $z$ axis by the Lorentz force and the rescattering is suppressed.

Figure 6(a) shows the initial transverse momentum distribution of the tunneling electron that will ultimately induce rescattering double ionization driven by a 3000 nm, $3\times 10^{14}$ W/cm$^2$ linearly polarized laser pulse by including the action of the magnetic field. For compensating the push of the Lorentz force, the initial transverse momentum distribution has been pushed oppositely to the laser propagation direction in order to meet the nucleus, which has been pointed by Walser $et$ $al$ and Daněk $et$ $al$ [62,63]. Besides the central point shifts to $p_z=-0.08$ a.u., the circles changes into ellipses, which is caused by the break of the cylindrical symmetry of the rescattering trajectories.

 figure: Fig. 6.

Fig. 6. The initial transverse momentum distributions of the tunneling electron that will lead to double ionization when (a) the laser field is linearly polarized and the magnetic field is included, (b) the laser field is elliptically polarized (with the ellipticity 0.005) and the magnetic field is not included, (c) same as (b) but the magnetic filed is included. (d) The correlated EMD for the laser condition of (c). The laser wavelength is 3000 nm, and the intensity is $3\times 10^{14}$ W/cm$^2$.

Download Full Size | PDF

Since the initial transverse momentum distribution of the tunnel-ionized electron has a gaussian shape centered at zero, as formulated by Eq. (1), the magnetic deviation will change the proportion of the first, second, and third rescattering, as well as the correlated EMD. In Fig. 1(c), the relative proportions for the first, second, and third rescattering are $6\%$, $45\%$ and $49\%$. In the nondipole case as showed in Fig. 1(d), however, the relative proportions are $9\%$, $37\%$ and $54\%$ for the three rescattering scenarios respectively. Consistently, one may notice that the events in the first and third quardrants in Fig. 1(c), mainly contributed by the second rescattering, are more than those in the first and third quardrants in Fig. 1(d).

If the driving laser pulse is not purely linearly polarized, which is the usual case in experiment, the electric component $E_y$ of the laser pulse will modify the electron trajectory in the $y$ direction. Accordingly, the initial transverse momentum distribution of the tunnel-ionized electron and the later on trajectory are modified, so does the correlated EMD. Figure 6(b) presents the initial transverse momentum distribution of the tunneling electron associated with the rescattering double ionization when the magnetic field is neglected and the ellipticity is $0.005$. Because of the drag along y-axis by the tiny electric field component $E_y$, the concentric circles presented in Fig. 5(a) are pulled apart and two set of broken ellipses are observed as expected. The relative proportions of different rescattering scenarios are also modified, i.e. $17\%$ for the first rescattering, $30\%$ for the second rescattering, and $53\%$ for the third rescattering, and we find (not shown) the corresponding correlated EMD is reshaped as predicted. When the unavoidable magnetic field is further introduced in the simulation, the initial transverse momentum distribution moves downward, as displayed in Fig. 6(c). The corresponding correlated EMD is shown in Fig. 6(d). The relative proportions for the first, second, and third rescattering are $26\%$, $30\%$ and $44\%$ respectively in Fig. 6(d). Since more double ionization events happen through the first and third rescattering, the main structure in the correlated EMD locates in the near-axis area and the pitcher and catcher have quite different momenta. Such correlated EMD looks similar to the experimental measurements [19,20].

Besides reshaping the correlated EMDs, the magnetic field also modifies the dication yields via altering the trajectories of rescattering electrons. Figure 7 shows the dication yield ratio of the nondipole and the dipole approximation cases as a function of laser wavelengths when the intensity is fixed at $3\times 10^{14}$ W/cm$^{2}$. The movement of the initial transverse momentum distribution away from the center of the gaussian shape given by Eq. (1) weakens the rescattering double ionization, which reduces the yield of dications. The magnetic field must be included though it is still far from relativistic regimes.

 figure: Fig. 7.

Fig. 7. The ratio of the dication yields for the nondipole and dipole approximation cases as a function of the laser wavelength. The driving lasers are linearly polarized and have the fixed intensity $3\times 10^{14}$ W/cm$^{2}$.

Download Full Size | PDF

4. Conclusions

To summarize, by following the electron trajectory in the combined Coulomb field and laser field of the driving mid-infrared laser, we found the first released electron may rescatter the parent ion within one optical cycle, or around 1.25 cycles, or 1.75 cycles after its tunneling. These three different rescattering scenarios are mainly governed by the initial transverse momentum at tunneling. Different rescattering events lead to well defined regions on the correlated EMD and help to understand the asymmetric energy sharing leading to double ionization. The inclusion of the unavoidable magnetic field and the slight ellipticity of the mid-IR laser pulse modifies the electron trajectories, weakens the rescattering double ionization, and reshapes the correlated EMD. Physically, due to the limitation of the numerical model, for example, the missing of the wave packet expansion during the electron propagation, the second and third rescattering events are overestimated. Nevertheless, our study provides qualitative understanding of a set of new electron trajectories leading to double ionization driven by mid-IR laser pulses, which are now available in advanced laboratories.

Funding

National Key Research and Development Program of China (2018YFA0404802); National Natural Science Foundation of China (11574205, 11721091, 91850203); Shanghai Municipal Education Commission (2017-01-07-0002-E00034); Ministerio de Economía y Competitividad (FIS2016-75652-P).

Disclosures

The authors declare that there are no conflicts of interest related to this article.

References

1. W. Becker, X. liu, P. J. Ho, and J. H. Eberly, “Theories of photoelectron correlation in laser-driven multiple atomic ionization,” Rev. Mod. Phys. 84(3), 1011–1043 (2012). [CrossRef]  

2. D. N. Fittinghoff, P. R. Bolton, B. Chang, and K. C. Kulander, “Observation of nonsequential double ionization of helium with optical tunneling,” Phys. Rev. Lett. 69(18), 2642–2645 (1992). [CrossRef]  

3. B. Walker, B. Sheehy, L. F. DiMauro, P. Agostini, K. J. Schafer, and K. C. Kulander, “Precision Measurement of Strong Field Double Ionization of Helium,” Phys. Rev. Lett. 73(9), 1227–1230 (1994). [CrossRef]  

4. Th. Weber, H. Giessen, M. Weckenbrock, G. Urbasch, A. Staudte, L. Spielberger, O. Jagutzki, V. Mergel, M. Vollmer, and R. Dörner, “Correlated electron emission in multiphoton double ionization,” Nature (London) 405(6787), 658–661 (2000). [CrossRef]  

5. R. Moshammer, B. Feuerstein, W. Schmitt, A. Dorn, C. D. Schröter, J. Ullrich, H. Rottke, C. Trump, M. Wittmann, G. Korn, K. Hoffmann, and W. Sandner, “Momentum Distributions of Nen+ Ions Created by an Intense Ultrashort Laser Pulse,” Phys. Rev. Lett. 84(3), 447–450 (2000). [CrossRef]  

6. A. Rudenko, V. L. B. de Jesus, Th. Ergler, K. Zrost, B. Feuerstein, C. D. Schröter, R. Moshammer, and J. Ullrich, “Correlated Two-Electron Momentum Spectra for Strong-Field Nonsequential Double Ionization of He at 800 nm,” Phys. Rev. Lett. 99(26), 263003 (2007). [CrossRef]  

7. A. Staudte, C. Ruiz, M. Schöffler, S. Schössler, D. Zeidler, Th. Weber, M. Meckel, D. M. Villeneuve, P. B. Corkum, A. Becker, and R. Dörner, “Binary and Recoil Collisions in Strong Field Double Ionization of Helium,” Phys. Rev. Lett. 99(26), 263002 (2007). [CrossRef]  

8. Y. Liu, S. Tschuch, A. Rudenko, M. Dürr, M. Siegel, U. Morgner, R. Moshammer, and J. Ullrich, “Strong-Field Double Ionization of Ar below the Recollision Threshold,” Phys. Rev. Lett. 101(5), 053001 (2008). [CrossRef]  

9. B. Bergues, M. K uel, N. G. Johnson, B. Fischer, N. Camus, K. J. Betsch, O. Herrwerth, A. Senftleben, A. M. Sayler, T. Rathje, I. Ben-Itzhak, R. R. Jones, G. G. Paulus, F. Krausz, R. Moshammer, J. Ullrich, and M. F. Kling, “Attosecond tracing of correlated electron-emission in non-sequential double ionization,” Nat. Commun. 3(1), 813 (2012). [CrossRef]  

10. H. P. Kang, K. Henrichs, M. Kunitski, Y. Wang, X. Hao, K. Fehre, A. Czasch, S. Eckart, L. Ph. H. Schmidt, M. Schöffler, T. Jahnke, X. Liu, and R. Dörner, “Timing Recollision in Nonsequential Double Ionization by Intense Elliptically Polarized Laser Pulses,” Phys. Rev. Lett. 120(22), 223204 (2018). [CrossRef]  

11. P. B. Corkum, “Plasma perspective on strong field multiphoton ionization,” Phys. Rev. Lett. 71(13), 1994–1997 (1993). [CrossRef]  

12. B. Feuerstein, R. Moshammer, D. Fischer, A. Dorn, C. D. Schröter, J. Deipenwisch, J. R. Crespo Lopez-Urrutia, C. Höhr, P. Neumayer, J. Ullrich, H. Rottke, C. Trump, M. Wittmann, G. Korn, and W. Sandner, “Separation of Recollision Mechanisms in Nonsequential Strong Field Double Ionization of Ar: The Role of Excitation Tunneling,” Phys. Rev. Lett. 87(4), 043003 (2001). [CrossRef]  

13. T. Fuji, N. Ishii, C. Y. Teisset, X. Gu, Th. Metzger, A. Baltuška, N. Forget, D. Kaplan, A. Galvanauskas, and F. Krausz, “Parametric amplification of few-cycle carrier-envelope phase-stable pulses at 2.1μm,” Opt. Lett. 31(8), 1103 (2006). [CrossRef]  

14. J. Tate, T. Auguste, H. G. Muller, P. Salières, P. Agostini, and L. F. DiMauro, “Scaling of Wave-Packet Dynamics in an Intense Midinfrared Field,” Phys. Rev. Lett. 98(1), 013901 (2007). [CrossRef]  

15. T. Popmintchev, M. Chen, D. Popmintchev, P. Arpin, S. Brown, S. Ališauskas, G. Andriukaitis, T. Balčiunas, O. D. Mücke, A. Pugzlys, A. Baltuška, B. Shim, S. E. Schrauth, A. Gaeta, C. Hernández-García, L. Plaja, A. Becker, A. Jaron-Becker, M. M. Murnane, and H. C. Kapteyn, “Bright Coherent Ultrahigh Harmonics in the keV X-ray Regime from Mid-Infrared Femtosecond Lasers,” Science 336(6086), 1287–1291 (2012). [CrossRef]  

16. Y. Huismans, A. Rouzée, A. Gijsbertsen, J. H. Jungmann, A. S. Smolkowska, P. S. W. M. Logman, F. Lépine, C. Cauchy, S. Zamith, T. Marchenko, J. M. Bakker, G. Berden, B. Redlich, A. F. G. van der Meer, H. G. Muller, W. Vermin, K. J. Schafer, M. Spanner, M. Yu. Ivanov, O. Smirnova, D. Bauer, S. V. Popruzhenko, and M. J. J. Vrakking, “Time-Resolved Holography with Photoelectrons,” Science 331(6013), 61–64 (2011). [CrossRef]  

17. C. T. L. Smeenk, L. Arissian, B. Zhou, A. Mysyrowicz, D. M. Villeneuve, A. Staudte, and P. B. Corkum, “Partitioning of the Linear Photon Momentum in Multiphoton Ionization,” Phys. Rev. Lett. 106(19), 193002 (2011). [CrossRef]  

18. A. Ludwig, J. Maurer, B. W. Mayer, C. R. Phillips, L. Gallmann, and U. Keller, “Breakdown of the Dipole Approximation in Strong-Field Ionization,” Phys. Rev. Lett. 113(24), 243001 (2014). [CrossRef]  

19. B. Wolter, M. G. Pullen, M. Baudisch, M. Sclafani, M. Hemmer, A. Senftleben, C. D. Schröter, J. Ullrich, R. Moshammer, and J. Biegert, “Strong-Field Physics with Mid-IR Fields,” Phys. Rev. X 5(2), 021034 (2015). [CrossRef]  

20. T. Shaaran, N. Camus, J. Dura, L. Fechner, A. Thai, A. Britz, M. Baudisch, T. Steinle, A. Senftleben, C. D. Schröter, J. Ullrich, T. Pfeifer, C. H. Keitel, J. Biegert, K. Z. Hatsagortsyan, and R. Moshammer, “Role of high ponderomotive energy in laser-induced nonsequential double ionization,” Phys. Rev. A 99(2), 023421 (2019). [CrossRef]  

21. E. Pisanty, D. D. Hickstein, B. R. Galloway, C. G. Durfee, H. C. Kapteyn, M. M. Murnane, and M. Ivanov, “High harmonic interferometry of the Lorentz force in strong mid-infrared laser fields,” New J. Phys. 20(5), 053036 (2018). [CrossRef]  

22. S. Brennecke and M. Lein, “High-order above-threshold ionization beyond the electric dipole approximation,” J. Phys. B 51(9), 094005 (2018). [CrossRef]  

23. A. Emmanouilidou, T. Meltzer, and P. B. Corkum, “Non-dipole recollision-gated double ionization and observable effects,” J. Phys. B 50(22), 225602 (2017). [CrossRef]  

24. A. Emmanouilidou and T. Meltzer, “Recollision as a probe of magnetic-field effects in nonsequential double ionization,” Phys. Rev. A 95(3), 033405 (2017). [CrossRef]  

25. H. R. Reiss, “Limits on Tunneling Theories of Strong-Field Ionization,” Phys. Rev. Lett. 101(4), 043002 (2008). [CrossRef]  

26. H. R. Reiss, “The tunnelling model of laser-induced ionization and its failure at low frequencies,” J. Phys. B 47(20), 204006 (2014). [CrossRef]  

27. J. Colgan and M. S. Pindzola, “Core-Excited Resonance Enhancement in the Two-Photon Complete Fragmentation of Helium,” Phys. Rev. Lett. 88(17), 173002 (2002). [CrossRef]  

28. Kitzler Markus, Milosevic Nenad, Scrinzi Armin, Krausz Ferenc, and Brabec Thomas, “Quantum Theory of Attosecond XUV Pulse Measurement by Laser Dressed Photoionization,” Phys. Rev. Lett. 88(17), 173904 (2002). [CrossRef]  

29. J. S. Parker, B. J. S. Doherty, K. T. Taylor, K. D. Schultz, C. I. Blaga, and L. F. DiMauro, “High-Energy Cutoff in the Spectrum of Strong-Field Nonsequential Double Ionization,” Phys. Rev. Lett. 96(13), 133001 (2006). [CrossRef]  

30. J. Feist, S. Nagele, R. Pazourek, E. Persson, B. I. Schneider, L. A. Collins, and J. Burgdörfer, “Nonsequential two-photon double ionization of helium,” Phys. Rev. A 77(4), 043420 (2008). [CrossRef]  

31. M. Lein, E. K. U. Gross, and V. Engel, “Intense-Field Double Ionization of Helium: Identifying the Mechanism,” Phys. Rev. Lett. 85(22), 4707–4710 (2000). [CrossRef]  

32. C. Ruiz, L. Plaja, L. Roso, and A. Becker, “Ab initio Calculation of the Double Ionization of Helium in a Few-Cycle Laser Pulse Beyond the One-Dimensional Approximation,” Phys. Rev. Lett. 96(5), 053001 (2006). [CrossRef]  

33. J. S. Prauzner-Bechcicki, K. Sacha, B. Eckhardt, and J. Zakrzewski, “Time-Resolved Quantum Dynamics of Double Ionization in Strong Laser Fields,” Phys. Rev. Lett. 98(20), 203002 (2007). [CrossRef]  

34. A. Zielinski, V. P. Majety, and A. Scrinzi, “Double photoelectron momentum spectra of helium at infrared wavelength,” Phys. Rev. A 93(2), 023406 (2016). [CrossRef]  

35. A. Becker and F. H. M. Faisal, “Correlated Keldysh-Faisal-Reiss theory of above-threshold double ionization of He in intense laser fields,” Phys. Rev. A 50(4), 3256–3264 (1994). [CrossRef]  

36. R. Kopold, W. Becker, H. Rottke, and W. Sandner, “Routes to Nonsequential Double Ionization,” Phys. Rev. Lett. 85(18), 3781–3784 (2000). [CrossRef]  

37. S. V. Popruzhenko and S. P. Goreslavskii, “Photoelectron momentum distribution for double ionization in strong laser fields,” J. Phys. B 34(8), L239–L246 (2001). [CrossRef]  

38. C. Figueira de Morisson Faria, X. Liu, W. Becker, and H. Schomerus, “Coulomb repulsion and quantum-classical correspondence in laser-induced nonsequential double ionization,” Phys. Rev. A 69(2), 021402 (2004). [CrossRef]  

39. T. Shaaran, K. Z. Hatsagortsyan, and C. H. Keitel, “Coulomb effect in laser-induced recollision excitation,” Phys. Rev. A 98(2), 023410 (2018). [CrossRef]  

40. R. Panfili, J. H. Eberly, and S. L. Haan, “Comparing classical and quantum dynamics of strong-field double ionization,” Opt. Express 8(7), 431 (2001). [CrossRef]  

41. S. L. Haan, J. S. Van Dyke, and Z. S. Smith, “Recollision Excitation, Electron Correlation, and the Production of High-Momentum Electrons in Double Ionization,” Phys. Rev. Lett. 101(11), 113001 (2008). [CrossRef]  

42. F. Mauger, C. Chandre, and T. Uzer, “From Recollisions to the Knee: A Road Map for Double Ionization in Intense Laser Fields,” Phys. Rev. Lett. 104(4), 043005 (2010). [CrossRef]  

43. Y. Zhou, C. Huang, Q. Liao, and P. Lu, “Classical Simulations Including Electron Correlations for Sequential Double Ionization,” Phys. Rev. Lett. 109(5), 053004 (2012). [CrossRef]  

44. C. Huang, Y Zhou, and P. Lu, “Recollision dynamics in nonsequential double ionization of atoms by long-wavelength pulses,” Opt. Express 24(25), 28361 (2016). [CrossRef]  

45. M. V. Ammosov, N. B. Delone, and V. P. Krainov, “Tunnel ionization of complex atoms and of atomic ions in an alternating electromagnetic field,” Sov. Phys. JETP 64, 1191 (1986).

46. T. Brabec, M. Y. Ivanov, and P. B. Corkum, “Coulomb focusing in intense field atomic processes,” Phys. Rev. A 54(4), R2551–R2554 (1996). [CrossRef]  

47. V. R. Bhardwaj, S. A. Aseyev, M. Mehendale, G. L. Yudin, D. M. Villeneuve, D. M. Raymer, M. Y. Ivanov, and P. B. Corkum, “Few Cycle Dynamics of Multiphoton Double Ionization,” Phys. Rev. Lett. 86(16), 3522–3525 (2001). [CrossRef]  

48. G. L. Yudin and M. Y. Ivanov, “Physics of correlated double ionization of atoms in intense laser fields: Quasistatic tunneling limit,” Phys. Rev. A 63(3), 033404 (2001). [CrossRef]  

49. J. Chen, J. Liu, L. B. Fu, and W. M. Zheng, “Interpretation of momentum distribution of recoil ions from laser-induced nonsequential double ionization by semiclassical rescattering model,” Phys. Rev. A 63(1), 011404 (2000). [CrossRef]  

50. A. Emmanouilidou, “Recoil collisions as a portal to field-assisted ionization at near-uv frequencies in the strong-field double ionization of helium,” Phys. Rev. A 78(2), 023411 (2008). [CrossRef]  

51. D. F. Ye, X. Liu, and J. Liu, “Classical Trajectory Diagnosis of a Fingerlike Pattern in the Correlated Electron Momentum Distribution in Strong Field Double Ionization of Helium,” Phys. Rev. Lett. 101(23), 233003 (2008). [CrossRef]  

52. G. P. Katsoulis, A. Hadjipittas, B. Bergues, M. F. Kling, and A. Emmanouilidou, “Slingshot Nonsequential Double Ionization as a Gate to Anticorrelated Two-Electron Escape,” Phys. Rev. Lett. 121(26), 263203 (2018). [CrossRef]  

53. L. B. Fu, J. Liu, and S. G. Chen, “Correlated electron emission in laser-induced nonsequence double ionization of helium,” Phys. Rev. A 65(2), 021406 (2002). [CrossRef]  

54. L. V. Keldysh, “Ionization in the field of a strong electromagnetic wave,” Sov. Phys. JETP 20, 1307 (1965).

55. N. I. Shvetsov-Shilovski, S. P. Goreslavski, S. V. Popruzhenko, and W. Becker, “Ellipticity effects and the contributions of long orbits in nonsequential double ionization of atoms,” Phys. Rev. A 77(6), 063405 (2008). [CrossRef]  

56. H. Kang, K. Henrichs, Y. Wang, X. Hao, S. Eckart, M. Kunitski, M. Schöffler, T. Jahnke, X. Liu, and R. Dörner, “Double ionization of neon in elliptically polarized femtosecond laser fields,” Phys. Rev. A 97(6), 063403 (2018). [CrossRef]  

57. S. Haan, Z. Smith, K. Shomsky, and P. Plantinga, “Anticorrelated electrons from weak recollisions in nonsequential double ionization,” J. Phys. B 41(21), 211002 (2008). [CrossRef]  

58. W. Becker, F. Grasbon, R. Kopold, D. Milosevic, G. Paulus, and H. Walther, “Above-threshold ionization: From classical features to quantum effects,” Adv. At., Mol., Opt. Phys. 48, 35–98 (2002). [CrossRef]  

59. Q. Z. Xia, J. F. Tao, J. Cai, L. B. Fu, and J. Liu, “Quantum Interference of Glory Rescattering in Strong-Field Atomic Ionization,” Phys. Rev. Lett. 121(14), 143201 (2018). [CrossRef]  

60. H. B. van Linden van den Heuvell and H. G. Muller, “Multiphoton Processes,” (Cambridge University Press, Cambridge, England, 1988).

61. Z. J. Chen, A. T. Le, T. Morishita, and C. D. Lin, “Quantitative rescattering theory for laser-induced high-energy plateau photoelectron spectra,” Phys. Rev. A 79(3), 033409 (2009). [CrossRef]  

62. M. W. Walser, C. H. Keitel, A. Scrinzi, and T. Brabec, “High Harmonic Generation Beyond the Electric Dipole Approximation,” Phys. Rev. Lett. 85(24), 5082–5085 (2000). [CrossRef]  

63. J. Daněk, K. Z. Hatsagortsyan, and C. H. Keitel, “Analytical approach to Coulomb focusing in strong-field ionization. I. Nondipole effects,” Phys. Rev. A 97(6), 063409 (2018). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. The correlated EMD when the laser wavelengths are 800 nm (a, b) and 3000 nm (c, d), and magnetic fields of the laser pulses are neglected (a, c) and included (b, d). The laser intensities are $3\times 10^{14}$ W/cm$^{2}$ for 800 nm and 3000 nm laser pulses. Both lasers are linearly polarized.
Fig. 2.
Fig. 2. (a) The laser electric field (black solid line) and the laser vector potential (amaranth dashed line). The electron tunnels out during the time interval [0.495T, 0.54T] will rescatter the parent ion in the time interval [1.25T, 1.45T] (the first rescattering), [1.55T, 1.9T] (the second rescattering) and [2.05T, 2,45T] (the third rescattering), and the corresponding double ionization probability are represented by the bars. (b) The incoming angle distribution for three rescattering events at different time. The variation in rescattering energy of these three events is depicted by green solid, dashed and dash-dotted lines. The laser wavelength is $\lambda =3000$ nm, and the laser intensity is $I=3\times 10^{14}$ W/cm$^2$.
Fig. 3.
Fig. 3. The time-evolved electron displacement (the first row) and electron kinetic momentum (the second row) along the laser polarization direction. The third and fourth rows are the time-evolved electron displacement and electron kinetic momentum along the laser propagation direction. The three columns from left to right are for the first, second and third rescattering, respectively.
Fig. 4.
Fig. 4. The typical trajectories for the first (green), second (red) and third (blue) rescattering without considering the magnetic field. The laser wavelength is $3000$ nm, and the intensity is $3\times 10^{14}$ W/cm$^2$. Surfaces are formed by rotating each trajectories with respect to the $x$-axis. The dashed curve is the trajectory starting with the identical displacement and momentum of the first rescattering trajectory but including the magnetic effect.
Fig. 5.
Fig. 5. (a) The initial transverse momentum distribution of the tunneling electron associated with the later double ionization. The magnetic field is not included. (b), (c) and (d): The correlated EMDs associated with the initial transverse momenta around 0, in the first and second circles, respectively. The laser wavelength is 3000 nm, and the intensity is $3\times 10^{14}$ W/cm$^{2}$. The laser is linearly polarized.
Fig. 6.
Fig. 6. The initial transverse momentum distributions of the tunneling electron that will lead to double ionization when (a) the laser field is linearly polarized and the magnetic field is included, (b) the laser field is elliptically polarized (with the ellipticity 0.005) and the magnetic field is not included, (c) same as (b) but the magnetic filed is included. (d) The correlated EMD for the laser condition of (c). The laser wavelength is 3000 nm, and the intensity is $3\times 10^{14}$ W/cm$^2$.
Fig. 7.
Fig. 7. The ratio of the dication yields for the nondipole and dipole approximation cases as a function of the laser wavelength. The driving lasers are linearly polarized and have the fixed intensity $3\times 10^{14}$ W/cm$^{2}$.

Equations (2)

Equations on this page are rendered with MathJax. Learn more.

ϖ ( p y , p z , t ) = 2 2 I p 1 | E ( t ) | exp [ ( p z 2 + p y 2 ) 2 I p 1 | E ( t ) | ] .
d 2 r i d t 2 = r i ( V n e i + V e e ) E ( t ) s v i × B ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.