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

Frustrated tunneling ionization in the elliptically polarized strong laser fields

Open Access Open Access

Abstract

We theoretically investigated frustrated tunneling ionization (FTI) in the interaction of atoms with elliptically polarized laser pulses by a semiclassical ensemble model. Our results show that the yield of frustrated tunneling ionization events exhibits an anomalous behavior which maximizes at the nonzero ellipticity. By tracing back the initial tunneling coordinates, we show that this anomalous behavior is due to the fact that the initial transverse velocity at tunneling of the FTI events is nonzero in the linear laser pulses and it moves across zero as the ellipticity increases. The FTI yield maximizes at the ellipticity when the initial transverse momentum for being trapped is zero. Moreover, the angular momentum distribution of the FTI events and its ellipticity dependence are also explored. The anomalous behavior revealed in our work is very similar to the previously observed ellipticity dependence of the near- and below-threshold harmonics, and thus our work may uncover the mechanism of the below-threshold harmonics which is still a controversial issue.

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

1. Introduction

Tunneling ionization is a fundamental process in the intense laser-atom/molecule interactions. When the electric field strength of the laser pulse becomes comparable to the binding Coulomb field inside the atom, a potential barrier is formed through which an electron can tunnel out of the atom [1,2]. The tunneled electron is then driven by the laser field and could return back to the parent ion [3], initiating various interesting processes, such as high-order above-threshold ionization [4,5], high-order harmonic generation [611], enhanced multiple ionization [12,13], etc. These processes have drawn considerable attention during the past decades and their applications have been widely explored. For example, the recollision of the tunneled electron could be used for attosecond pulse generation [14,15], molecular orbital tomography [16], laser-induced electron diffraction [17] and time-resolved strong-field photoelectron holography [1823].

When atoms are exposed to the intense laser fields, it has also been observed that a large population of neutral atoms survive in Rydberg states [2426]. The underlying mechanism for creation of these Rydberg states has attracted much interest. Two mechanisms for forming the Rydberg state have been proposed. In the multiphoton ionization region, the Rydberg states are created via AC Stark-shifted multiphoton resonant excitation during the laser pulse [24,25,2729]. Recent experiments with coincident measurement of the correlated ion and electron provided solid evidence of this multiphoton resonance excitation for the Rydberg atoms by the UV light [30]. In the tunneling ionization regime, it was demonstrated that the Rydberg states are formed through the frustrated tunneling ionization (FTI) process [26]. In this process, the tunneled electron does not gain enough drift energy from the laser field thus it is finally recaptured by the Coulomb field of the parents atomic [26,31] or molecular ion [3238]. Initially, it was proposed that rescattering is responsible for the recapture of the tunneled electron [26]. Subsequent studies showed that rescattering is not necessary for FTI and the electron could directly launch into the elliptical orbit and finally stays in the highly excited state after the end of the laser pulses [31,39]. This FTI model has very successfully reproduced the n-distribution (energy) of the Rydberg states [26,40,41] and has initiated a variety of studies in laser-atom interactions, such as strong-field acceleration of neutral atoms [4244] and the survival of Rydberg states [45] and so on [46,47]. With the semiclassical model, the details of the FTI process have been revealed [26,39,48]. This FTI scenario is supported by experiments on the ellipticity dependence of the excited neutral atom yield [26,39].

There is an another comparable phenomena in the interaction of a strong laser field with atoms and molecules, known as below-threshold harmonic generation [49]. Different from the high-order harmonics wherein the underlying dynamics has been well understood, the mechanism for the below-threshold harmonics is still under debate. During the past decade, it has drawn increasing theoretical as well as experimental attention, partly due to its potential application for the coherent extreme-ultraviolet light sources [50,51]. Recently, the correspondence between the below-threshold harmonic generation and FTI has been investigated with a semiclassical model [52]. It is shown that these two processes share the same series of trajectories. By solving the time-dependent Schrödinger equation, it is shown that FTI and below-threshold harmonics have a similar carrier-envelope phase dependence on the driving pulses [52], indicating that FTI should be responsible for the below-threshold harmonic generation. Two decades ago, an intriguing phenomena in below- and near-threshold harmonics was reported [5356]. It was found that the yield of those harmonics exhibited an anomalous dependence on the ellipticity of the driving field, i.e., the yield maximized at a nonzero ellipticity. For FTI, this anomalous behavior should also exist if it shares the same process as the below-threshold harmonic generation. However, previous studies have shown that the Rydberg atom yield decreases monotonously with the ellipticity [26,39].

In this work, we study FTI in the interaction of the atom with the elliptically polarized laser pulses. We observed an anomalous ellipticity dependence of the Rydberg atom yield, which maximizes at a nonzero ellipticity. By tracing back the initial tunneling coordinates, we show that this anomalous behavior is due to the fact that the initial transverse velocity at tunneling of the FTI events is nonzero in the linear laser pulses and it moves across zero as the ellipticity increases. The FTI yield maximizes at the ellipticity when the initial transverse momentum for being trapped is zero. The presence of this anomalous behavior depends on the laser wavelength, intensity and pulse duration. This explains why this anomalous behavior has not been observed in previous studies. Moreover, the angular momentum distribution of the FTI events and its ellipticity dependence are also explored. The anomalous behavior revealed in our work is very similar to the previously observed ellipticity dependence of the near- and below-threshold harmonics, thus our work may uncover the mechanism of the below-threshold harmonics.

2. Theoretical model

In this work, we employ the classical-trajectory Monte Carlo (CTMC) ensemble model [26,48] to study the FTI driven by elliptically polarized laser pulses. This CTMC ensemble model has been well established and has been widely used in studying FTI [26,39,48,57]. In this model, FTI process includes two steps: tunneling ionization and classical evolution of the tunneled electron in a combination of the laser field and Coulomb potential of the ion.

In the first step, i.e., tunneling ionization, the electron is launched at the tunnel exit with a zero initial longitudinal momentum and a Gaussian-like initial transverse momentum distribution. The weight of the trajectory was determined by the tunneling theory [5860] (atomic units are used unless stated otherwise)

$$P(t_0,p_\perp)\propto \exp[-\frac{2(2I_p)^{3/2}}{3|F(t_0)|}] \exp[-p_\perp^2\frac{\sqrt{2I_p}}{|F(t_0)|}].$$
Here, $F(t_0)$ is the instantaneous electric field of the laser field at tunneling time $t_0$, $p_\perp$ is the initial transverse momentum, which is perpendicular to the transient electric field $\textbf {F}(t_0)$ of the laser, and $I_p$ is the ionization energy. After tunneling, the evolution of the electron’s trajectory is determined by Newton’s equation of motion
$$\ddot{\textbf{r}}(t)=-\textbf{F}(t)-{\nabla} V(\textbf{r}),$$
where $V(r)=-1/\sqrt {r^2+a^2}$ is the interaction potential between the ion and electron, and we choose the soft core parameter $a= 0.01$ to avoid the Coulomb singularity. $F(t)$ is the electric field of the laser pulses. In this paper, we study FTI in the elliptically polarized laser pulses where the electric field is written as
$$\textbf{F}(t)=\cos^2(\frac{\omega t}{2N})\frac{E_0}{\sqrt{\xi^2+1}}[\cos(\omega t)\widehat{\textbf{x}}+\xi \sin(\omega t)\widehat{\textbf{y}}].$$
Here $E_0$, $\omega$, N, and $\xi$ are the electric field amplitude, frequency, pulse duration (N denotes the number of optical cycles), and ellipticity of the laser pulses, respectively. Following previous studies [39,48], we use the 2D classical model, i.e., the motion of the electron is confined in the laser polarization plane.

In our calculations, we consider FTI of Ar($I_p=0.59$ a.u.). The tunneling exit is approximately estimated by $\textbf {r}_0=-I_p/\textbf {F}(t_0)$. The ensemble is obtained by sampling four millions of classical trajectories over the tunneling ionization time and initial transverse momentum. In our calculations, we only consider the events where tunnel ionization occurs around one electric field peak of the laser pulse, $-T/4<t_0<T/4$, where T is the optical cycle of the laser field. For the long pulses, the electric fields of the laser are similar from cycle to cycle and thus the electron dynamics in different cycles should be similar. So, the results demonstrated in this paper will not change when more laser cycles for tunneling ionization are considered. We examine the energy of the electron when the laser pulse is turned off. The electrons with the final energy $E_f =v^2/2-1/r>0$ stand for the ionized ones. Those with $E_f<0$ correspond to the tunneled electrons recaptured by the ion, which are the FTI events that we focus on in this study. We estimate the main quantum number (n) of the bound electron with the Rydberg formula $E_f=-0.5/n^2$ [26]. The angular momentum is calculated as $\textbf {L}_f =\textbf {r}_f\times \textbf {P}_f$. In this study, we focus on the ellipticity dependence of yield and the angular momentum distribution of the FTI events.

3. Ellipticity dependence of the Rydberg-atom yields

Figure 1 displays the yields of the excited atoms as a function of the ellipticity of the laser pulses. Here the pulse contains 30 optical cycles (i.e., N=30). We performed the calculations at different laser wavelengths and intensities. For the 1600 nm field at the intensity of $8\times 10^{13}$ W/cm$^2$ (curve c), the yield decreases rapidly when the ellipticity increases. This result is consistent with previous theoretical and experimental results [26,39]. This ellipticity dependence is a solid evidence that FTI is responsible for the excited atoms in strong laser fields. For the wavelength of 800 nm pulse with the intensity of $1.6\times 10^{14}$ W/cm$^2$ (curve b), the yield almost does not change as the ellipticity varies from 0 to 0.2 and then decrease sharply as the ellipticity further increases. At the lower laser intensity, as displayed by the blue circles (curve a), the yield first increases and then decreases with the ellipticity. The yield curve shows a surprising peak at the ellipticity of 0.2.

 figure: Fig. 1.

Fig. 1. The yields of excited atoms as a function of the ellipticity of the laser pulses. The wavelengths and the laser intensities are specified in the legend.

Download Full Size | PDF

In order to understand the underlying dynamics for this anomalous ellipticity dependence, we take the advantage of the CTMC model and trace back to the initial coordinates of the FTI events. In Fig. 2, we show the distributions of the tunneling ionization time and initial transverse momentum of the FTI events at different ellipticities, correspond to the anomalous curve a in Fig. 1. It can be seen that the area of the Rydberg-state atoms is crescent-shaped. The ionization time and the initial transverse momentum for these recaptured electrons are correlated. For the linear laser field ($\xi$=0), the electrons emitted around the crest electric field should possess a proper nonzero initial transverse momentum to be recaptured, and these with nearly zero initial transverse momentum are emitted before the peak of the electric field, as shown in Fig. 2(a). Here we divided the Rydberg area into three regions, AI, AII and B. Regions AI and AII correspond to the nearly vertical distribution in the initial coordinates and Region B corresponds to the nearly horizontal area. As the ellipticity increases, it is shown that the shape of the Rydberg area keeps almost unchanged while the location shifts to the right side gradually as the ellipticity changes from 0 to 0.2. Thus, the probability of the distribution in region AII and region B decrease because the required initial transverse momentum increases, while the probability in region AI increases due to the reduced initial transverse momentum requirement. The competition of the increasing probability in region AII and the decreasing probability in regions AII and B results in the different ellipticity-dependence of the three FTI yield curves in Fig. 1. As the ellipticity further increases, region AI passes through the $p_\perp =0$ plane, and thus the probability decreases in all of the regions, resulting in the rapid decreasing of the excited atom yield for ellipticity larger than 0.2.

 figure: Fig. 2.

Fig. 2. Probability distributions of the FTI events on the coordinates of the tunneling ionization time $t_0$ and initial transverse momentum $P_\perp$ at different ellipticities (from 0 to 0.4) and same laser parameters. The laser intensity (I=$8\times 10^{13}$W/cm$^2$) and wavelength (800nm) are same with the curve a of Fig. 1. The colorbar on the right represents the recapture probability. The increasing pulse ellipticity leads to the crescent-shaped Rydberg area shifting to the right gradually. The black crosses indicate the coordinate of ($t_0$=0, $P_\perp$=0). In Fig. 2(a), the dashed red line divides the Rydberg area into three regions, vertical regions AI, AII and horizontal region B.

Download Full Size | PDF

To explore the ellipticity dependence of yield for different laser parameters shown in Fig. 1, we display the distributions of the ionization time and the initial transverse momentum at different laser wavelengths and intensities in Fig. 3. The laser parameters in Figs. 3(a)–3(c) are the same as these for the three curves in Fig. 1, respectively. As explained in Fig. 2, the location of the Rydberg area shifts toward the right side when the ellipticity increases. Thus, the yield from region AI increases while the yields in region B and AII decrease. In Fig. 3(a), regions AI and AII have the main contribution to the total Rydberg atom yield. As the ellipticity increases, region AI moves toward $p_\perp$=0 and regions AII and B move away from $p_\perp$=0, and thus the capture probability in region AI increases and those in regions AII and B decrease. The increase of the capture probability in AI is stronger than the decrease in regions AII and B, resulting in the increase of the total FTI yield. As the ellipticity further increases, the yield in region AI decreases and thus the total yield decreases, resulting in an anomalous peak in the yield curve.

 figure: Fig. 3.

Fig. 3. Probability distributions of the FTI events in the initial tunneling coordinates for the linear laser pulses($\xi$=0). The laser wavelengths and intensities are (a) 800 nm and $8\times 10^{13}$ W/cm$^2$, (b) 800 nm and $1.6\times 10^{14}$ W/cm$^2$, (c) 1600 nm and 8$\times 10^{13}$ W/cm$^2$. The colorbar on the right represents the recapture probability. The dashed black lines are the level curves of the probability of FTI events. The dashed red line divides the Rydberg area into three regions, vertical regions AI, AII and horizontal region B.

Download Full Size | PDF

As the laser intensity increases, as shown in Fig. 3(b), the location of the region B is more close to $t_0$=0 and thus the weight of region B increases. So, region AI and region B, AII have the similar contribution to the total yield. As the ellipticity increases, the yield increases in region AI and decreases in regions B, AII. The competition of these three regions leads to total yield almost unchanged with the ellipticity between 0$-$0.2. At the longer wavelength, as shown in Fig. 3(c), the location of region B is very close to the peak of the electric field, and thus it has the dominant contribution to the total Rydberg atom yield. The yield in region B decreases monotonously with the ellipticity and thus the total yield exhibits the sharp decreasing behavior, as shown by the squares (curve c) in Fig. 1.

We mention that the behavior of the Rydberg atom yield on ellipticity in strong field ionization does not only depend on the laser intensity and wavelength [61], but also on the pulses duration [57]. We repeated our calculation at a short pulse duration with N=10 and N=4, as shown in Fig. 4. For the shorter pulse duration, the area of region B increases and thus the contribution to the total yield increase, as shown in Fig. 4(a) and Fig. 4(c). Consequently, the yield decreases with the laser ellipticity, as shown in Fig. 4(b) and Fig. 4(d). This is in consistent with previous theoretical and experimental results [26,39].

 figure: Fig. 4.

Fig. 4. (a, c) The distribution of FTI events in the initial tunneling coordinates and (b, d) the yield of FTI events for laser fields with shorter pulse duration. The laser intensity and wavelength are $8\times 10^{13}$W/cm$^2$ and 800 nm. The pulse duration in (a) (b) is N=10 and in (c) (d) is N=4.

Download Full Size | PDF

The calculations above show that the anomalous behavior is more obvious at the nonadiabatic tunneling region. In nonadiabatic tunneling, it has been shown that the most likely transverse momentum in tunneling by the elliptical laser pulse is nonzero [62]. However, this most probable nonzero initial transverse momentum is much less than the corresponding initial transverse momentum of region AI. So, as ellipticity increases, region AI will move towards the most probable initial transverse momentum and the probability increases. Therefore, the yield of the Rydberg atoms increases and the anomalous behavior will persist when the nonadiabatic effect is considered in our calculations. Of course, the exact ellipticity for the location of the peak will change slightly due to the nonadiabatic effect.

We notice that in the below- and near-threshold harmonic generation, a similar anomalous ellipticity dependence of the yield has been experimentally observed 20 years ago [5356]. The mechanism of the below- and near-threshold harmonic generation is still a hot and controversial topic. Recently, it was suggested that there is a close correspondence between the below-threshold harmonics and FTI [52]. Our results revealed another surprising resemblance on the ellipticity dependence between FTI and near- and below-threshold harmonics. According to the analysis above, the anomalous ellipticity dependence results from the nonzero initial transverse momentum of the responsible trajectories in the linear laser pulses, which moves towards zero momentum as the ellipticity increases. Thus, our results indicate that the emission of the below- and near-threshold harmonics in the linear laser pulse might be also related to the trajectories of the nonzero initial transverse momentum at tunneling. In previous studies, the dynamics for the near- and below-threshold harmonics have been analyzed with the one-dimensional classical trajectories model [63]. Our results indicate that the initial transverse momentum should be considered and thus the one-dimensional analysis should be inadequate. Because of the nonzero initial transverse momentum, the electron trajectory should be an ellipse even in the linear driving pulses (as will be detailedly shown in the subsection below). It implies that the near- and below-threshold harmonics should possess nonzero ellipticity even in the linear driving pulses, which could be easily verified in further experiments. Our results on FTI may help us to recover the mechanism of below-and near-threshold harmonics.

4. Angular momentum distributions

Most of the previous studies on FTI are focused on the yields and the energy distribution of Rydberg atoms. The angular momentum of the excited atoms is also a very important quantity of Rydberg atoms. The angular momentum of Rydberg states is of great importance for both fundamental and applied physics. For example, the low angular-momentum Rydberg states were used for optical quantum information manipulations [64,65] and long-lived circular Rydberg levels with high angular momentum are tools in the exploration of cavity quantum electrodynamic effects [66,67]. Here, we study the angular momentum distribution of the Rydberg atoms in FTI with the semiclassical model.

In Fig. 5, we show the angular momentum number distribution of Rydberg-state atoms at different ellipticities. Here,the final angular momentum number is calculated as $L =[\textbf {L}_f] =[\textbf {r}_f\times \textbf {P}_f$], where [x] denotes the nearest integer of x. The laser wavelength is 800 nm and the intensity is $8\times 10^{13}$ W/cm$^2$. It is shown that the angular momentum number distribution exhibits a symmetric triple-hump structure for this linear laser field ($\xi$=0). As the ellipticity increases, the right and the middle peak are suppressed and the distribution exhibits a single peak in the left side.

 figure: Fig. 5.

Fig. 5. The angular momentum number distribution of Rydberg-state atoms at different ellipticities. The laser intensity (I=$8\times 10^{13}$W/cm$^2$) and wavelength (800nm) are same with the curve a of Fig. 1. It exhibits a symmetric triple-hump structure for the linear laser field ($\xi$=0).

Download Full Size | PDF

In order to understand the dependence of the angular momentum distribution on the ellipticity, we trace back to the initial coordinates of the FTI events. Fig. 6(a) shows the angular momentum number distribution in the coordinates of initial transverse momentum and tunneling time at $\xi$=0. Note that here we show the absolute value of L. It is shown that angular quantum number distributed exhibits stripes nearly parallel to the coordinate of tunneling time. It means that the final angular momentum is sensitive to the initial transverse velocity but is almost independent on the tunneling time. Thus, it can be expected that there was a correlation between the final angular momentum and the initial transverse momentum. In Fig. 6(b), we also show the principal quantum number (n) distribution in the initial tunneling coordinates. This distribution exhibits the crescent shape with the smaller principle quantum numbers located at the inner edge, the larger n at outer edge of the crescent. The reason for this crescent shape has been explained in a recent paper [57]. The distributions for the angular momentum and the energy in the initial tunneling coordinates are totally different and this difference results in the different ellipticity dependence. As shown in Fig. 2, the location of the distribution moves towards right with the shape unchanged as the ellipticity increases. Therefore, the weight of Rydberg state atomic in the left side of Fig. 6(a) increases and this in the center and right side decreases. Consequently, the center and right peak in Fig. 5 are suppressed as ellipticity increases. On the contrary, the energy distribution keeps almost unchanged with the laser ellipticity (not shown here) but can be controlled via the intensities, wavelengths and duration of the laser field [28,57,6870].

 figure: Fig. 6.

Fig. 6. (a) The distribution of angular momentum number L in the coordinates of initial transverse momentum and tunneling time at $\xi$=0. (b) The distribution of principal quantum number n in the coordinates of initial transverse momentum and tunneling time at $\xi$=0. Different colors represent different numbers L and n.

Download Full Size | PDF

Moreover, the relationship between the final angular momentum and the initial transverse momentum can be easily established using the ’guiding center’ model [71]. The classical trajectory dynamics is governed by the Hamiltonian

$$H(\textbf{r},\textbf{p},t)=\frac{|\textbf{p}|^2}{2}+V(\textbf{r})+\textbf{r}\cdot\textbf{E}(t),$$
At the lowest order of the perturbative expansion, the electron coordinates are of the form [71]
$$\textbf{r}=\textbf{r}_g+\frac{\textbf{E}(t)}{\omega^2},$$
$$\textbf{p}=\textbf{p}_g+\textbf{A}(t),$$
where ($\textbf {r}_g,\textbf {p}_g$) are the canonically conjugate variables of the guiding center, and $\textbf {A}(t)$ is the vector potential of the laser field. The guiding-center dynamics is governed by the averaged Hamiltonian [71]
$$\bar{H}(\textbf{r}_g,\textbf{p}_g)=\frac{|\textbf{p}_g|^2}{2}+V_{eff}(\textbf{r}).$$
This Hamiltonian no longer depends on time, as a result of averaging. At the lowest order in the perturbative expansion, $V_{eff}(\textbf {r})=V(\textbf {r})$. Consequently, its energy $\varepsilon =\bar {H}(\textbf {r}_g,\textbf {p}_g)$ and the angular momentum of the guiding center is conserved.

In Fig. 7, we show two classical trajectories of the FTI events with the initial coordinated indicated by the black cross in Figs. 2(f) and 2(g). The guiding center trajectory which averages the disturbance of the electric field is in good agreement with the CTMC trajectory after the laser field turned off. Thus we can use the guiding-center trajectory to estimate the angular momentum. Due to the conservation of angular momentum of the guiding center trajectory, we can predict the relationship between the final angular momentum and the initial transverse momentum. We obtain the initial conditions of the corresponding guiding center according Eqs. (5-6) [71], when t0=0

$$\begin{aligned} \textbf{r}_g & =(\frac{I_p}{E_0}+\frac{E_0}{\omega^2})\widehat{\textbf{x}} = \frac{E_0}{\omega^2}(\frac{\gamma^2}{2}+1)\widehat{\textbf{x}}\\ \textbf{p}_g & =P_\perp \widehat{\textbf{y}}. \end{aligned}$$
Then, we get the final angular momentum
$$\begin{aligned} \textbf{L}_f & =\textbf{r}_f\times\textbf{P}_f =\textbf{r}_g\times\textbf{P}_g\\ & =\frac{E_0}{\omega^2}(\frac{\gamma^2}{2}+1)P_\perp \widehat{\textbf{z}}, \end{aligned}$$
where $\gamma =\sqrt {\frac {I_p}{2U_p}}$ is the Keldysh parameter [1] and $U_p=\frac {E_0^2}{4\omega ^2}$ is the ponderomotive energy. As shown in Fig. 8, this expression agrees well with the final angular momentum distribution versus the initial transverse momentum at tunneling.

 figure: Fig. 7.

Fig. 7. (a) A sample classical trajectories of the FTI events (the blue solid curve) and the corresponding guiding center trajectory (the dashed cyan curve) at $\xi =0.25$. (b) The same trajectories at $\xi =0.30$. (a) and (b) have the same initial conditions, $t_0=0, P_\perp =0$. The black crosses show the initial position of classical trajectory, and the black dot marks the zero point of the coordinate.

Download Full Size | PDF

 figure: Fig. 8.

Fig. 8. The distribution of the final angular momentum ($L_f$) versus the initial transverse momentum at tunneling ($P_\perp$). The dashed magenta line shows the analytical formula (9). The black cross marks the zero point of the coordinate. The colorbar on the right represents the recapture probability.

Download Full Size | PDF

5. Conclusions

Using CTMC simulations we found that the yield of frustrated tunneling ionization events exhibits an anomalous behavior which maximizes at the nonzero ellipticity. This anomalous behavior is very similar to the previously observed ellipticity dependence of the near- and below-threshold harmonics. By tracing back the initial tunneling coordinates, we show that this anomalous behavior is because the initial transverse velocity at tunneling of the FTI events is nonzero in the linear laser pulses and it moves across zero as the ellipticity increases. The FTI yield maximizes at the nonzero ellipticity when the initial transverse momentum for being trapped is zero. We show that the presence of this anomalous behavior depends on the laser intensity, wavelength and pulse duration.

Moreover, the angular momentum distribution of the FTI events and its ellipticity dependence are also explored. The angular momentum distribution of Rydberg state changes from a symmetric triple-hump structure at the linear laser field to a single peak structure as the ellipticity increases. By tracing back the initial tunneling coordinates, we found that the angular quantum number does not depend on the tunneling time but is very sensitive to the initial transverse momentum at tunnel. Using guiding center theory, we obtained the correlation between the final angular momentum and the initial transverse momentum.

Funding

National Natural Science Foundation of China (NSFC) (11604108, 11622431, 11627809, 11874163); Program for HUST Academic Frontier Youth Team.

Acknowledgments

Numerical simulations presented in this paper were carried out using the High Performance Computing Center experimental testbed in SCTS/CGCL (see http://grid.hust.edu.cn/hpcc).

References

1. L. V. Keldysh, “Ionization in the Field of a Strong Electromagnetic Wave,” Sov. Phys. JETP 20, 1307 (1965).

2. S. Luo, M. Li, W. Xie, K. Liu, Y. Feng, B. Du, Y. Zhou, and P. Lu, “Exit momentum and instantaneous ionization rate of nonadiabatic tunneling ionization in elliptically polarized laser fields,” Phys. Rev. A 99(5), 053422 (2019). [CrossRef]  

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

4. G. G. Paulus, W. Nicklich, H. Xu, P. Lambropoulos, and H. Walther, “Plateau in above threshold ionization spectra,” Phys. Rev. Lett. 72(18), 2851–2854 (1994). [CrossRef]  

5. W. Becker, F. Grasbon, R. Kopold, D. B. Miloševic, G. G. Paulus, and H. Walther, “Above-threshold ionization: From classical features to quantum effects,” Adv. At. Mol. Opt. Phys. 48, 35–98 (2002).

6. M. Lewenstein, P. Balcou, M. Y. Ivanov, A. L’Huillier, and P. B. Corkum, “Theory of high-harmonic generation by low-frequency laser fields,” Phys. Rev. A 49(3), 2117–2132 (1994). [CrossRef]  

7. F. Krausz and M. Ivanov, “Attosecond physics,” Rev. Mod. Phys. 81(1), 163–234 (2009). [CrossRef]  

8. J. Li, Q. Zhang, L. Li, X. Zhu, T. Huang, P. Lan, and P. Lu, “Orientation dependence of high-order harmonic generation in nanowire,” Phys. Rev. A 99(3), 033421 (2019). [CrossRef]  

9. Y. He, L. He, P. Lan, B. Wang, L. Li, X. Zhu, W. Cao, and P. Lu, “Direct imaging of molecular rotation with high-order-harmonic generation,” Phys. Rev. A 99(5), 053419 (2019). [CrossRef]  

10. L. Liang, P. Lan, X. Zhu, T. Huang, Q. Zhang, M. Lein, and P. Lu, “Reciprocal-Space-Trajectory Perspective on High-Harmonic Generation in Solids,” Phys. Rev. Lett. 122, 193901 (2019). [CrossRef]  

11. P. Meng, D. Zhao, and W. Liu, “Topological plasmonic modes in graphene-coated nanowire arrays,” Opt. Quantum Electron. 51, 156 (2019). [CrossRef]  

12. W. Becker, X. J. 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]  

13. A. Tong, Q. Li, X. Ma, Y. Zhou, and P. Lu, “Internal collision induced strong-field nonsequential double ionization in molecules,” Opt. Express 27(5), 6415 (2019). [CrossRef]  

14. G. Sansone, E. Benedetti, F. Calegari, C. Vozzi, L. Avaldi, R. Flammini, L. Poletto, P. Villoresi, C. Altucci, R. Velotta, S. Stagira, S. De Silvestri, and M. Nisoli, “Isolated Single-Cycle Attosecond Pulses,” Science 314(5798), 443–446 (2006). [CrossRef]  

15. E. Goulielmakis, M. Schultze, M. Hofstetter, V. S. Yakovlev, J. Gagnon, M. Uiberacker, A. L. Aquila, E. M. Gullikson, D. T. Attwood, R. Kienberger, F. Krausz, and U. Kleineberg, “Single-Cycle Nonlinear Optics,” Science 320(5883), 1614–1617 (2008). [CrossRef]  

16. J. Itatani, J. Levesque, D. Zeidler, H. Niikura, H. Pépin, J. C. Kieffer, P. B. Corkum, and D. M. Villeneuve, “Tomographic imaging of molecular orbitals,” Nature 432(7019), 867–871 (2004). [CrossRef]  

17. C. I. Blaga, J. Xu, A. D. DiChiara, E. Sistrunk, K. Zhang, P. Agostini, T. A. Miller, L. F. DiMauro, and C. D. Lin, “Imaging ultrafast molecular dynamics with laser-induced electron diffraction,” Nature 483(7388), 194–197 (2012). [CrossRef]  

18. 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. Y. 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]  

19. Y. Zhou, O. I. Tolstikhin, and T. Morishita, “Near-Forward Rescattering Photoelectron Holography in Strong-Field Ionization: Extraction of the Phase of the Scattering Amplitude,” Phys. Rev. Lett. 116(17), 173001 (2016). [CrossRef]  

20. J. Tan, Y. Zhou, M. He, Y. Chen, Q. Ke, J. Liang, X. Zhu, M. Li, and P. Lu, “Determination of the Ionization Time Using Attosecond Photoelectron Interferometry,” Phys. Rev. Lett. 121(25), 253203 (2018). [CrossRef]  

21. M. Li, H. Xie, W. Cao, S. Luo, J. Tan, Y. Feng, B. Du, W. Zhang, Y. Li, Q. Zhang, P. Lan, Y. Zhou, and P. Lu, “Photoelectron Holographic Interferometry to Probe the Longitudinal Momentum Offset at the Tunnel Exit,” Phys. Rev. Lett. 122(18), 183202 (2019). [CrossRef]  

22. J. Tan, Y. Zhou, M. He, Q. Ke, J. Liang, Y. Li, M. Li, and P. Lu, “Time-resolving tunneling ionization via strong-field photoelectron holography,” Phys. Rev. A 99(3), 033402 (2019). [CrossRef]  

23. Y. Liu, J. Tan, M. He, H. Xie, Y. Qin, Y. Zhao, M. Li, Y. Zhou, and P. Lu, “Photoelectron holographic interferences from multiple returning in strong-field tunneling ionization,” Opt. Quantum Electron. 51(5), 145 (2019). [CrossRef]  

24. M. P. de Boer and H. G. Muller, “Observation of large populations in excited states after short-pulse multiphoton ionization,” Phys. Rev. Lett. 68(18), 2747–2750 (1992). [CrossRef]  

25. R. R. Jones, D. W. Schumacher, and P. H. Bucksbaum, “Population trapping in Kr and Xe in intense laser fields,” Phys. Rev. A 47(1), R49–R52 (1993). [CrossRef]  

26. T. Nubbemeyer, K. Gorling, A. Saenz, U. Eichmann, and W. Sandner, “Strong-Field Tunneling without Ionization,” Phys. Rev. Lett. 101(23), 233001 (2008). [CrossRef]  

27. G. N. Gibson, L. Fang, and B. Moser, “Direct femtosecond laser excitation of the $2p$ state of H by a resonant seven-photon transition in ${\textrm{H}}_{2}^{+}$,” Phys. Rev. A 74(4), 041401 (2006). [CrossRef]  

28. Q. Li, X. Tong, T. Morishita, C. Jin, H. Wei, and C. D. Lin, “Rydberg states in the strong field ionization of hydrogen by 800, 1200 and 1600 nm lasers,” J. Phys. B 47(20), 204019 (2014). [CrossRef]  

29. H. Zimmermann, S. Patchkovskii, M. Ivanov, and U. Eichmann, “Unified Time and Frequency Picture of Ultrafast Atomic Excitation in Strong Laser Fields,” Phys. Rev. Lett. 118(1), 013003 (2017). [CrossRef]  

30. W. Zhang, X. Gong, H. Li, P. Lu, F. Sun, Q. Ji, K. Lin, J. Ma, H. Li, J. Qiang, F. He, and J. Wu, “Electron-nuclear correlated multiphoton-route to Rydberg fragments of molecules,” Nature Communications 10, 757 (2019).

31. H. Liu, Y. Liu, L. Fu, G. Xin, D. Ye, J. Liu, X. T. He, Y. Yang, X. Liu, Y. Deng, C. Wu, and Q. Gong, “Low Yield of Near-Zero-Momentum Electrons and Partial Atomic Stabilization in Strong-Field Tunneling Ionization,” Phys. Rev. Lett. 109(9), 093001 (2012). [CrossRef]  

32. B. Manschwetus, T. Nubbemeyer, K. Gorling, G. Steinmeyer, U. Eichmann, H. Rottke, and W. Sandner, “Strong Laser Field Fragmentation of ${\textrm{H}}_{2}$: Coulomb Explosion without Double Ionization,” Phys. Rev. Lett. 102(11), 113002 (2009). [CrossRef]  

33. T. Nubbemeyer, U. Eichmann, and W. Sandner, “Excited neutral atomic fragments in the strong-field dissociation of N2 molecules,” J. Phys. B 42(13), 134010 (2009). [CrossRef]  

34. J. McKenna, S. Zeng, J. J. Hua, A. M. Sayler, M. Zohrabi, N. G. Johnson, B. Gaire, K. D. Carnes, B. D. Esry, and I. Ben-Itzhak, “Frustrated tunneling ionization during laser-induced D2 fragmentation: Detection of excited metastable D* atoms,” Phys. Rev. A 84(4), 043425 (2011). [CrossRef]  

35. J. Wu, A. Vredenborg, B. Ulrich, L. P. H. Schmidt, M. Meckel, S. Voss, H. Sann, H. Kim, T. Jahnke, and R. Dörner, “Multiple Recapture of Electrons in Multiple Ionization of the Argon Dimer by a Strong Laser Field,” Phys. Rev. Lett. 107(4), 043003 (2011). [CrossRef]  

36. J. McKenna, A. M. Sayler, B. Gaire, N. G. Kling, B. D. Esry, K. D. Carnes, and I. Ben-Itzhak, “Frustrated tunnelling ionization during strong-field fragmentation of ${\textrm{D}}_{3}^{+}$,” New J. Phys. 14(10), 103029 (2012). [CrossRef]  

37. A. von Veltheim, B. Manschwetus, W. Quan, B. Borchers, G. Steinmeyer, H. Rottke, and W. Sandner, “Frustrated Tunnel Ionization of Noble Gas Dimers with Rydberg-Electron Shakeoff by Electron Charge Oscillation,” Phys. Rev. Lett. 110(2), 023001 (2013). [CrossRef]  

38. W. Zhang, Z. Yu, X. Gong, J. Wang, P. Lu, H. Li, Q. Song, Q. Ji, K. Lin, J. Ma, H. Li, F. Sun, J. Qiang, H. Zeng, F. He, and J. Wu, “Visualizing and Steering Dissociative Frustrated Double Ionization of Hydrogen Molecules,” Phys. Rev. Lett. 119(25), 253202 (2017). [CrossRef]  

39. A. S. Landsman, A. N. Pfeiffer, C. Hofmann, M. Smolarski, C. Cirelli, and U. Keller, “Rydberg state creation by tunnel ionization,” New J. Phys. 15(1), 013001 (2013). [CrossRef]  

40. H. Zimmermann, J. Buller, S. Eilzer, and U. Eichmann, “Strong-Field Excitation of Helium: Bound State Distribution and Spin Effects,” Phys. Rev. Lett. 114(12), 123003 (2015). [CrossRef]  

41. S. Larimian, S. Erattupuzha, C. Lemell, S. Yoshida, S. Nagele, R. Maurer, A. Baltuška, J. Burgdörfer, M. Kitzler, and X. Xie, “Coincidence spectroscopy of high-lying Rydberg states produced in strong laser fields,” Phys. Rev. A 94(3), 033401 (2016). [CrossRef]  

42. U. Eichmann, T. Nubbemeyer, H. Rottke, and W. Sandner, “Acceleration of neutral atoms in strong short-pulse laser fields,” Nature 461(7268), 1261–1264 (2009). [CrossRef]  

43. S. Eilzer and U. Eichmann, “Steering neutral atoms in strong laser fields,” J. Phys. B 47(20), 204014 (2014). [CrossRef]  

44. H. Zimmermann and U. Eichmann, “Atomic excitation and acceleration in strong laser fields,” Phys. Scr. 91(10), 104002 (2016). [CrossRef]  

45. U. Eichmann, A. Saenz, S. Eilzer, T. Nubbemeyer, and W. Sandner, “Observing Rydberg Atoms to Survive Intense Laser Fields,” Phys. Rev. Lett. 110(20), 203002 (2013). [CrossRef]  

46. D. G. Arbó, C. Lemell, and J. Burgdörfer, “Signatures of tunneling and multiphoton ionization by short-laser pulses: The partial-wave distribution,” J. Phys.: Conf. Ser. 635, 012003 (2015).

47. S. Larimian, C. Lemell, V. Stummer, J.-W. Geng, S. Roither, D. Kartashov, L. Zhang, M.-X. Wang, Q. Gong, L.-Y. Peng, S. Yoshida, J. Burgdörfer, A. Baltuška, M. Kitzler, and X. Xie, “Localizing high-lying Rydberg wave packets with two-color laser fields,” Phys. Rev. A 96(2), 021403 (2017). [CrossRef]  

48. N. I. Shvetsov-Shilovski, S. P. Goreslavski, S. V. Popruzhenko, and W. Becker, “Capture into rydberg states and momentum distributions of ionized electrons,” Laser Phys. 19(8), 1550–1558 (2009). [CrossRef]  

49. W. Xiong, L. Peng, and Q. Gong, “Recent progress of below-threshold harmonic generation,” J. Phys. B 50(3), 032001 (2017). [CrossRef]  

50. D. C. Yost, T. R. Schibli, J. Ye, J. L. Tate, J. Hostetter, M. B. Gaarde, and K. J. Schafer, “Vacuum-ultraviolet frequency combs from below-threshold harmonics,” Nat. Phys. 5(11), 815–820 (2009). [CrossRef]  

51. H. Yun, J. H. Mun, S. I. Hwang, S. B. Park, I. A. Ivanov, C. H. Nam, and K. T. Kim, “Coherent extreme-ultraviolet emission generated through frustrated tunnelling ionization,” Nat. Photonics 12(10), 620–624 (2018). [CrossRef]  

52. W.-H. Xiong, X.-R. Xiao, L.-Y. Peng, and Q. Gong, “Correspondence of below-threshold high-order-harmonic generation and frustrated tunneling ionization,” Phys. Rev. A 94(1), 013417 (2016). [CrossRef]  

53. N. H. Burnett, C. Kan, and P. B. Corkum, “Ellipticity and polarization effects in harmonic generation in ionizing neon,” Phys. Rev. A 51(5), R3418–R3421 (1995). [CrossRef]  

54. K. Miyazaki and H. Takada, “High-order harmonic generation in the tunneling regime,” Phys. Rev. A 52(4), 3007–3021 (1995). [CrossRef]  

55. M. Kakehata, H. Takada, H. Yumoto, and K. Miyazaki, “Anomalous ellipticity dependence of high-order harmonic generation,” Phys. Rev. A 55(2), R861–R864 (1997). [CrossRef]  

56. H. Soifer, P. Botheron, D. Shafir, A. Diner, O. Raz, B. D. Bruner, Y. Mairesse, B. Pons, and N. Dudovich, “Near-Threshold High-Order Harmonic Spectroscopy with Aligned Molecules,” Phys. Rev. Lett. 105(14), 143904 (2010). [CrossRef]  

57. L. Ortmann, C. Hofmann, and A. S. Landsman, “Controlling the quantum number distribution and yield of Rydberg states via the duration of the laser pulse,” arXiv:1810.07164 (2018).

58. A. M. Perelomov, V. S. Popov, and M. V. Terent’ev, “Ionization of Atoms in an Alternating Electric Field: II,” Sov. Phys. JETP 24, 207 (1967).

59. 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).

60. N. B. Delone and V. P. Krainov, “Energy and angular electron spectra for the tunnel ionization of atoms by strong low-frequency radiation,” J. Opt. Soc. Am. B 8(6), 1207 (1991). [CrossRef]  

61. L. Ortmann, C. Hofmann, and A. S. Landsman, “Dependence of Rydberg-state creation by strong-field ionization on laser intensity,” Phys. Rev. A 98(3), 033415 (2018). [CrossRef]  

62. V. D. Mur, S. V. Popruzhenko, and V. S. Popov, “Energy and momentum spectra of photoelectrons under conditions of ionization by strong laser radiation (The case of elliptic polarization),” J. Exp. Theor. Phys. 92(5), 777–788 (2001). [CrossRef]  

63. J. A. Hostetter, J. L. Tate, K. J. Schafer, and M. B. Gaarde, “Semiclassical approaches to below-threshold harmonics,” Phys. Rev. A 82(2), 023401 (2010). [CrossRef]  

64. T. Peyronel, O. Firstenberg, Q. Liang, S. Hofferberth, A. V. Gorshkov, T. Pohl, M. D. Lukin, and V. Vuletic, “Quantum nonlinear optics with single photons enabled by strongly interacting atoms,” Nature 488(7409), 57–60 (2012). [CrossRef]  

65. L. Li and A. Kuzmich, “Quantum memory with strong and controllable Rydberg-level interactions,” Nat. Commun. 7(1), 13618 (2016). [CrossRef]  

66. C. Guerlin, J. Bernu, S. Deléglise, C. Sayrin, S. Gleyzes, S. Kuhr, M. Brune, J.-M. Raimond, and S. Haroche, “Progressive field-state collapse and quantum non-demolition photon counting,” Nature 448(7156), 889–893 (2007). [CrossRef]  

67. C. Sayrin, I. Dotsenko, X. Zhou, B. Peaudecerf, T. Rybarczyk, S. Gleyzes, P. Rouchon, M. Mirrahimi, H. Amini, M. Brune, J.-M. Raimond, and S. Haroche, “Real-time quantum feedback prepares and stabilizes photon number states,” Nature 477(7362), 73–77 (2011). [CrossRef]  

68. I. Burenkov, A. Popov, O. Tikhonova, and E. Volkova, “New features of interaction of atomic and molecular systems with intense ultrashort laser pulses,” Laser Phys. Lett. 7(6), 409–434 (2010). [CrossRef]  

69. E. A. Volkova, A. M. Popov, and O. V. Tikhonova, “Ionization and stabilization of atoms in a high-intensity, low-frequency laser field,” J. Exp. Theor. Phys. 113(3), 394–406 (2011). [CrossRef]  

70. M. V. Fedorov, N. P. Poluektov, A. M. Popov, O. V. Tikhonova, V. Y. Kharin, and E. A. Volkova, “Interference Stabilization Revisited,” IEEE J. Sel. Top. Quantum Electron. 18(1), 42–53 (2012). [CrossRef]  

71. J. Dubois, S. A. Berman, C. Chandre, and T. Uzer, “Capturing Photoelectron Motion with Guiding Centers,” Phys. Rev. Lett. 121(11), 113202 (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 (8)

Fig. 1.
Fig. 1. The yields of excited atoms as a function of the ellipticity of the laser pulses. The wavelengths and the laser intensities are specified in the legend.
Fig. 2.
Fig. 2. Probability distributions of the FTI events on the coordinates of the tunneling ionization time $t_0$ and initial transverse momentum $P_\perp$ at different ellipticities (from 0 to 0.4) and same laser parameters. The laser intensity (I=$8\times 10^{13}$W/cm$^2$) and wavelength (800nm) are same with the curve a of Fig. 1. The colorbar on the right represents the recapture probability. The increasing pulse ellipticity leads to the crescent-shaped Rydberg area shifting to the right gradually. The black crosses indicate the coordinate of ($t_0$=0, $P_\perp$=0). In Fig. 2(a), the dashed red line divides the Rydberg area into three regions, vertical regions AI, AII and horizontal region B.
Fig. 3.
Fig. 3. Probability distributions of the FTI events in the initial tunneling coordinates for the linear laser pulses($\xi$=0). The laser wavelengths and intensities are (a) 800 nm and $8\times 10^{13}$ W/cm$^2$, (b) 800 nm and $1.6\times 10^{14}$ W/cm$^2$, (c) 1600 nm and 8$\times 10^{13}$ W/cm$^2$. The colorbar on the right represents the recapture probability. The dashed black lines are the level curves of the probability of FTI events. The dashed red line divides the Rydberg area into three regions, vertical regions AI, AII and horizontal region B.
Fig. 4.
Fig. 4. (a, c) The distribution of FTI events in the initial tunneling coordinates and (b, d) the yield of FTI events for laser fields with shorter pulse duration. The laser intensity and wavelength are $8\times 10^{13}$W/cm$^2$ and 800 nm. The pulse duration in (a) (b) is N=10 and in (c) (d) is N=4.
Fig. 5.
Fig. 5. The angular momentum number distribution of Rydberg-state atoms at different ellipticities. The laser intensity (I=$8\times 10^{13}$W/cm$^2$) and wavelength (800nm) are same with the curve a of Fig. 1. It exhibits a symmetric triple-hump structure for the linear laser field ($\xi$=0).
Fig. 6.
Fig. 6. (a) The distribution of angular momentum number L in the coordinates of initial transverse momentum and tunneling time at $\xi$=0. (b) The distribution of principal quantum number n in the coordinates of initial transverse momentum and tunneling time at $\xi$=0. Different colors represent different numbers L and n.
Fig. 7.
Fig. 7. (a) A sample classical trajectories of the FTI events (the blue solid curve) and the corresponding guiding center trajectory (the dashed cyan curve) at $\xi =0.25$. (b) The same trajectories at $\xi =0.30$. (a) and (b) have the same initial conditions, $t_0=0, P_\perp =0$. The black crosses show the initial position of classical trajectory, and the black dot marks the zero point of the coordinate.
Fig. 8.
Fig. 8. The distribution of the final angular momentum ($L_f$) versus the initial transverse momentum at tunneling ($P_\perp$). The dashed magenta line shows the analytical formula (9). The black cross marks the zero point of the coordinate. The colorbar on the right represents the recapture probability.

Equations (9)

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

P ( t 0 , p ) exp [ 2 ( 2 I p ) 3 / 2 3 | F ( t 0 ) | ] exp [ p 2 2 I p | F ( t 0 ) | ] .
r ¨ ( t ) = F ( t ) V ( r ) ,
F ( t ) = cos 2 ( ω t 2 N ) E 0 ξ 2 + 1 [ cos ( ω t ) x ^ + ξ sin ( ω t ) y ^ ] .
H ( r , p , t ) = | p | 2 2 + V ( r ) + r E ( t ) ,
r = r g + E ( t ) ω 2 ,
p = p g + A ( t ) ,
H ¯ ( r g , p g ) = | p g | 2 2 + V e f f ( r ) .
r g = ( I p E 0 + E 0 ω 2 ) x ^ = E 0 ω 2 ( γ 2 2 + 1 ) x ^ p g = P y ^ .
L f = r f × P f = r g × P g = E 0 ω 2 ( γ 2 2 + 1 ) P z ^ ,
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.