1 Introduction
The ionization of atoms and molecules in strong laser pulses can trigger a variety of strong field phenomena, such as high-order above-threshold ionization [
1,
2], high-order harmonic and attosecond pulse generation [
3,
4], and nonsequential double ionization [
5,
6]. Compared with a linearly polarized laser pulse, there is an additional transverse electric field component for an elliptically polarized laser field, which suppresses the rescattering of the ionized electron with its parent ion. The resulting strong field phenomena are dramatically different from those in linearly polarized laser fields [
7–
11]. Therefore, the ionization in elliptically polarized laser fields has attracted much interest in recent years, and it has been shown to be important for the attoclock principle [
12,
13], which is a powerful scheme to resolve attosecond electron dynamics in atoms and molecules. In the attoclock method [
14,
15], precisely calibrating the ellipticity is a difficult task in the experiment [
16]. Moreover, most of previous attoclock studies focused on the ionization rate of the electron wave packet (WP) in elliptically polarized laser fields [
13,
17–
19] and its phase distribution was ignored. The phase distribution of the electron WP [classical action of the strong-field approximation (SFA), which is a function of the electron momentum] is the key to understand the wavelike property of the electron [
20] and thus it is the basis for matter wave interferometry. As far as we know, the phase distribution of the electron WP was only studied in circularly polarized laser fields [
21,
22], and was not studied in elliptically polarized laser fields before.
In this work, we use an interferometic scheme to extract the phase distribution of the electron WP from above-threshold ionization (ATI) in an elliptically polarized laser field. We use another electron WP released from a circularly polarized laser field as a reference wave to interfere with the electron WP ionized by the elliptically polarized laser field, generating a photoelectron vortex structure in the photoelectron momentum distribution (PMD). Initially, this photoelectron vortex structure, an unusual kind of Ramsey interference [
23] of laser induced electron WPs, was discovered and investigated in the single photon ionization by a pair of time-delayed counter-rotating circularly polarized laser pulses [
24]. Subsequently, the vortex structures were studied in various processes, including in multiphoton single ionization of He atoms [
25], in single-photon double ionization of He [
26] and H
2 [
27], in two-color single ionization of molecular ions [
28,
29], and in ionization of negative hydrogen ions [
30]. Experimentally, spiral vortices were confirmed in the ionization of potassium atoms [
21,
31] and sodium atoms [
32]. Using the vortex-shaped momentum distribution, the electronic displacement in a ultrashort laser pulse was accurately extracted [
22].
We find that, different from the vortex structure in a pair of circularly polarized laser pulses, the interference fringe of the vortex structure produced by a circularly polarized laser pulse and a time-delayed counter-rotating elliptically polarized laser pulse is unevenly spaced in the photoelectron angular distribution. From the modulation of the fringe spacing with the electron emission angle, we can extract the angle-dependent phase distribution of the electron WP in the elliptically polarized laser pulse with high precision. The results from the SFA are in a good agreement with the numerical solutions of the time-dependent Schördinger equation (TDSE) using a short-range potential. By artificially screening the ionic potential at different ranges when solving the TDSE, we find that the angle-dependent phase distribution in the elliptically polarized laser fields reveals a clear angular shift with respect to the SFA, whose value is the same as the attoclock shift. We further find that the amplitude of the angle-dependent phase distribution is sensitive to the laser ellipticity, which could be used to calibrate the pulse ellipticity with high precision.
2 Theoretical methods
2.1 SFA theory
In the SFA simulation, we use two time-delayed counter-rotating laser pulses to ionize a hydrogen atom, which consists of a left circularly polarized pulse (defined from the point of view of the receiver)
and a time-delayed right elliptically polarized pulse
where E0 denotes the peak amplitude, is the ellipticity of the time-delayed elliptically polarized pulse, is the time delay between the two pulses, and ω is the angular frequency of the laser pulse ( a.u. corresponds to the laser wavelength of 400 nm). The laser pulse envelope is employed with a duration of , where T0 is the period of the 400-nm laser field. We use atomic units (a.u.) throughout unless specified otherwise.
In the SFA, the transition amplitude from the ground state to the Volkov state with final momentum
p is written as [
33]
where
with being the ionization potential of the hydrogen atom. In the simulation, Eq. (3) is directly integrated without relying on the saddle-point approximation.
The final doubly differential momentum distribution of the photoelectron is given by
where px and py represent the final momenta along the x and y directions, respectively. The intensity of the laser field is set to be 3×1013 W/cm2 in the simulation.
2.2 Numerical solutions of the TDSE
We also numerically solve the two-dimensional TDSE in the single-active-electron approximation. Using the length gauge, the TDSE can be written as
where is the electron wave function, denotes the electron position, and is the electric field of the laser pulse.
In this paper, we use a soft-core Coulomb potential and a short-range potential in the TDSE calculation, with a being the soft parameter and b being the screening length of the short-range potential. As the screening length tends to infinity, the short-range potential becomes the Coulomb potential. For the case of the Coulomb potential, a is chosen to be 0.64 to match the ionization potential of the hydrogen atom. For the short-range potential, a is varied for different b to ensure the same binding energy for the Coulomb potential and short-range potential.
The split-operator spectral method is used on a Cartesian grid [
34] ranges from −500 to 500 a.u. to numerically solve the two-dimensional TDSE. The time step is fixed at
a.u. and the spatial discretization is
a.u.. The wave function in the whole two-dimensional space is divided into the inner wave function
and the outer wave function
smoothly by a splitting technique [
35]. The wave function in the momentum space
is calculated by the Fourier transformation of the outer wave function. Finally, we obtain the final momentum distribution related to the sum of the wave function in the momentum space
for each time step
τ,
Here,
F is electron energy associated with
p,
θ is the angle of the emitted electron with respect to the
x axis, and
is the vector potential of the laser field. The initial wave function is prepared by imaginary-time propagation [
36].
3 Results and discussion
3.1 Extracting the phase distribution of the electron wave packet
We first show in Fig. 1(a) the PMD of the hydrogen atom in two time-delayed counter-rotating circularly polarized pulses calculated by the SFA. One can see many ring-like ATI structures centered at around zero momentum [marked by the white Arabic numerals in Fig. 1(a)], which come from the intercycle interference of the electron WPs [
37,
38]. Notably, there is a left-helicity momentum vortex structure appearing within each ATI ring. From the view of multiphoton ionization, those electron vortex structures can be interpreted as a result of the interference of the free-electron WPs released by the two time-delayed laser pulses with different magnetic quantum numbers
m [
31]. In the left circularly polarized laser pulse, the phase distribution of the ionized electron WP as a function of the emission angle
can be written as
(
for right circularly polarized laser pulse) [
22], where
is the electron emission angle with respect to the
x axis [shown by the white arrow in Fig. 1(a)], and
N is the number of photons absorbed by the atom for the corresponding ATI ring [i.e.,
, where
Up is the ponderomotive potential,
]. Because the carrier frequency is 0.114 a.u. in our work, five photons are absorbed for the first order ATI. As shown in Fig. 1(a), there are ten spiral arms for the first order ATI ring in the PMD, which doubles the number of the photons absorbed in the ionization process [
22,
24]. For the second-order ATI, the number of the spiral arms for the vortex structure is 12, which also doubles the number of the absorbed photons.
To study the phase distribution of the ionized electron WP in the elliptically polarized laser field, we take the electron WP released from a circularly polarized laser pulse as a reference, which interferes with the electron WP ionized by the time-delayed counter-rotating elliptically polarized laser pulse. The resulting PMD calculated by the SFA is shown in Fig. 1(b). Compared with Fig. 1(a), the ATI spectra show maximal ionization rate along the y axis (minor axis of the laser ellipse of the elliptically polarized laser pulse). There are also many vortex structures within each order ATI. The number of the spiral arms for the vortex structure within each order ATI is the same as that in Fig. 1(a). We next focus on the vortex structure within the first-order ATI ring. Figure 1(c) presents the photoelectron angular distribution within the first-order ATI ring at a.u., as marked by a red-dashed circle in Fig. 1(b). Clearly, there is an oscillation with ten minima in the photoelectron angular distribution, which corresponds to five-photon absorption for the first order ATI ring. Thus, the angular difference [shown by the red dash-dotted line in Fig. 1(c)] between two neighboring minima (adjacent angular fringe spacing) should be closed to . However, a closer inspection reveals that the angular fringe spacing is larger around than that around . The nonuniform spaced distribution means that the angle-dependent phase distribution in the elliptically polarized laser field does not linearly depend on .
Following the scheme described in Ref. [
22], we can extract the angle-dependent phase distribution
of the ionized electron WP from the nonuniform spaced vortex patterns. We denote the phase difference between the elliptically and circularly polarized laser fields as
, where
is the angle-dependent phase distribution of the electron WP in a circularly polarized laser field. Because there is only a difference of
between
and
,
is also referred to as the angle-dependent phase distribution in the following text. For circularly polarized laser field,
. The phase difference for the vortex interference structure between
and
is given by
, where
and
represent the phase distribution of the electron WP ionized by the first and second laser pulses, respectively. In Fig. 1(c), the spacing of adjacent interference minima of the photoelectron angular distribution corresponding to a
value change. So the phase difference satisfies
, where the subscript
k labels the number of the interference minimum in the photoelectron angular distribution. Therefore, the adjacent angular fringe spacing
can be expressed as
where is the phase difference of the electron WPs for adjacent interference minima. From the adjacent angular fringe spacing in Fig. 1(c), we can obtain the phase distribution with respect to according to Eq. (8). Thus the angle-dependent phase distribution in the elliptically polarized laser field is directly obtained from the PMD. The extracted angle-dependent phase distribution from Figs. 1(a) and (b) are shown by the orange dotted line and the blue dashed line in Fig. 1(d), respectively. In the SFA method, we can also obtain the angle-dependent phase distribution from the transition amplitude MP [Eq. (3)] for a single elliptically polarized laser field. Note that in a real experiment, the transition amplitude in the SFA is not a direct observable. The angle-dependent phase distribution obtained from the transition amplitude of the SFA is shown by the red solid line in Fig. 1(d). One can see that the angle-dependent phase distribution extracted from the PMD agrees well with that obtained from the transition amplitude of a single elliptically polarized laser field. For comparison, in Fig. 1(d), we also show the electron yield of the first-order ATI (integrated over a momentum range of 0.26 a.u.<p<0.44 a.u.) in the elliptically polarized laser field by the green solid line with arbitrary units. It is interesting that the oscillation of the angle-dependent phase distribution show a clear phase shift of with respect to the oscillation of the angle-dependent electron yield.
To validate our interferometic scheme, we compare the angle-dependent phase distributions extracted from the PMD and obtained from the transition amplitude for different orders of ATI and different laser ellipticities using the SFA method. In Fig. 2, the angle-dependent phase distributions in the left column are extracted from the PMD while in the right column are obtained from the transition amplitude of the SFA. Figures 2(a) and (b) present the angle-dependent phase distributions for the first-order and second-order ATI ring, which are denoted by the blue dash-dotted and red dotted curves, respectively. Clearly, the amplitude of the angle-dependent phase distribution decreases with the increase of the ATI order. In Figs. 2(c) and (d), the angle-dependent phase distributions in the elliptically polarized laser field with the ellipticities of 0.93, 0.95 and 0.97 are represented by the red solid, blue dash-dotted and green dotted curves, respectively. We find that, as the laser ellipticity increases, the amplitude of the angle-dependent phase distribution decreases, and tends to zero (the case of circular polarization). For different orders of ATI and different laser ellipticities, the angle-dependent phase distributions extracted from the vortex-shaped PMDs are in a good agreement with those obtained from the transition amplitude.
To interpret the angular phase shift between the oscillation of the photoelectron yield and the oscillation of angle-dependent phase distribution in Fig. 1(d), we express the phase distribution of the ionized electron WP in an elliptically polarized laser field by the quasiclassical action [
39],
where t0 is the ionization time of the electron. We take a half cycle of the laser field as an example. When (major axis of the laser ellipse), the phase distribution is the same as that of the circularly polarized laser field, i.e., the angle-dependent phase distribution is equal to zero, while the angle-dependent yield distribution reveals a maximum. When the phase distribution S is maximized, the ionization time obtained by solving Eq. (9) is around . This means that there is an angular phase shift of between the oscillation of the photoelectron yield and the oscillation of angle-dependent phase distribution, which is consistent with Fig. 1(d).
3.2 The effect of the ionic potential on the angle-dependent phase distribution
In shed light on the effect of the ionic potential on the phase distribution of the ionized electron wave packet, we show in Figs. 3(a) and (b) the vortex-shaped momentum distribution by solving the TDSE with the short-range potential (
b=3) and the Coulomb potential. As shown in Fig. 3(a), the result of the short range potential agrees well with the SFA. Compared with Fig. 3(a), the whole momentum spectra in Fig. 3(b) rotate counterclockwise under the influence of the Coulomb potential [
40,
41]. Figure 3(c) shows the normalized electron yield of the first-order ATI (the data are integrated over a momentum range of 0.26 a.u.<
p<0.44 a.u.) with respect to the angle
for the Coulomb potential and the short-range potential. Obviously, with considering the effect of the Coulomb potential, the maximum of the electron yield exhibit a significant angular shift, i.e., attoclock shift [indicated by two dashed vertical lines in Fig. 3(c)], which is the main observable in the attoclock experiment [
12,
13]. Figure 3(d) shows the photoelectron angular distribution within the first-order ATI ring at
a.u. for the Coulomb potential and the short-range potential with
, which extracted from the vortex-shaped PMD of Figs. 2(a) and (b). From Fig. 3(d), one can see that the adjacent angular fringe spacing is changed under the influence of the Coulomb potential, which means that the ionic potential has a significant effect on the angle-dependent phase distribution in the elliptically polarized laser field. This is different from the case of circular polarization, in which the angle-dependent phase distribution is not influenced by the ionic potential due to the rotational symmetry of the laser field [
22].
The extracted angle-dependent phase distributions for the first-order ATI ring at a.u. from the vortex-shaped PMD are shown by the blue dotted lines in Fig. 4. The PMD were obtained by solving the TDSE at different screening lengths. It is clear from Fig. 4 that there is an increasing angular shift for the angle-dependent phase distribution with the increase of the screening length. To quantitatively study this angular shift, we define as the angular shift [indicated by two dashed vertical lines in Fig. 4(c)] between the maxima of the angle-dependent phase distributions calculated by the TDSE (blue dashed lines) and the SFA (the red solid lines). The angular shift for the angle-dependent phase distribution with respect to the screening length b is shown by the blue solid circles in Fig. 5(a). One can see that the angular shift increases with the increase of the screening length. For comparison, we also show the attoclock shift for the first-order ATI as a function of the screening length b by the red solid line with triangles. One notices that the angular shift of the angle-dependent phase distribution is the same as the attoclock shift. Moreover, in Fig. 4, we can see that the amplitude of the angle-dependent phase distribution shows a dependence on the screening length. For [Fig. 4(a)], The angle-dependent phase distribution calculated by the TDSE agrees with the SFA calculation. As can be seen from Figs. 4(b)–(d), the amplitude of the angle-dependent phase distribution [Figs. 4(b) and (c)] become larger for the screening lengths of 5 a.u. and 9 a.u., while it becomes smaller for the Coulomb potential [Fig. 4(d)] as compared to the SFA result. In Fig. 5(b), we show the amplitude of the angle-dependent phase distribution as a function of the screening length b for the ellipticity values of 0.91, 0.93, 0.95, and 0.97. One finds that the amplitude of the angle-dependent phase distribution first increases and then decreases with increasing the screening length of the ionic potential for all ellipticity values.
To qualitatively understand how the ionic potential influences the amplitude of angle-dependent phase distribution, we study the quasiclassical action along the electron trajectory. Without considering the effect of the ionic potential, the action of the electron WP can be expressed as
, while with considering the effect of the ionic potential, the action is replaced by
, where
and
. The ionic potential affects both terms of
V and
T, and the contribution of
V to the action
S is opposite to that of
T [
42].
For the short-range potential with small value of b (), the effect of the ionic potential is ignored, so the amplitude of angle-dependent phase distribution is the same as the SFA calculation [Fig. 4(a)]. As the screening length increases, the predominant contribution (T or V) to the action S is changed, which might lead to the nonmonotonic dependence of the amplitude on the screening length. Furthermore, in Fig. 5(b), though the ellipticity is slightly changed, the amplitude of the angle-dependent phase distribution has a considerable change for both short-range potential and Coulomb potential. This means the amplitude of the angle-dependent phase distribution depends sensitively on the laser ellipticity, which provides an alternative way to precisely calibrate the laser ellipticity in the attoclock measurement.
4 Conclusion
In conclusion, we use an interferometic scheme to extract the phase distribution of the electron WP from the ATI in the elliptically polarized laser fields by modulating adjacent angular fringe spacing from the vortex-shaped photoelectron interference in the PMD. The results from the SFA are in a good agreement with the numerical solutions of the TDSE with a short range potential. Furthermore, we reveal the significance of the ionic potential on the angle-dependent phase distribution by artificially screening the ionic potential at different ranges. We find that the angle-dependent phase distribution shows an angular shift with respect to the SFA, and the value agrees well with the attoclock shift. Furthermore, the amplitude of the angle-dependent phase distribution is found to be sensitive to the ellipticity of the laser pulse. Our work provides a feasible tool for calibrating the ellipticity of the laser pulse and for understanding the wavelike property of the particle. Moreover, our results suggest that accurate descriptions of the strong-field ionization require complete characterization of both amplitude and phase of the continuum WP. This will impact the analysis of the time delay in strong-field ionization, which is believed to be closely related to the ionization phase. Looking ahead, because the ionization phase distribution is usually sensitive to the electron motion in atoms and molecules, accurate knowledge of the phase distribution of the continuum WP will pave the way to controlling the electron motion on the natural scale.