Helicity reversion in high harmonic generation driven by bichromatic counterrotating circularly polarized laser fields
We investigate the polarization properties of high harmonics generated with the bichromatic counterrotating circularly polarized (BCCP) laser fields by numerically solving time-dependent Schrödinger equation (TDSE). It is found that, the helicity of the elliptically polarized harmonic emission is reversed at particular harmonic orders. Based on the time-frequency analysis and the classical three-step model, the correspondence between the positions of helicity reversions and the classical trajectories of continuum electrons is established. It is shown that, the electrons ionized at one lobe of laser field can be divided into different groups based on the different lobes they recombine at, and the harmonics generated by adjacent groups have opposite helicities. Our study performs a detailed analysis of high harmonics in terms of electron trajectories and depicts a clear and intuitive physical picture of the HHG process in BCCP laser field.
pacs:42.50.Hz, 32.80.Qk, 33.80.Wz, 32.80.Wr
High-harmonic generation (HHG) occurs in the interaction of intense femtosecond laser field with atoms or molecules PBCorkum (); FKrausz (). HHG has been intensively studied in experimental and theoretical investigations Zhu2 (); Xin (); Zhu3 (), because it promises both applications for generating attosecond pulses Sansone (); LM (); Cavalieri () and for probing the molecular structure Itatani (); Haessler (); Kraus (); Zhai () and ultrafast dynamics Smirnova (); Baker (); HJW (); Travnikova (); Boeglin (). The HHG process can be nicely explained by the three-step model Corkum (); Schafer (). First, under the influence of the laser field, an electron escapes from the atom or molecule and reaches the continuum. Then it is accelerated in the laser field. Finally the accelerated electron is driven back to the parent ion and a high energy photon is emitted. According to this picture, only laser fields, which are able to drive the electrons back to the ion, can be used to efficiently generate high harmonics.
Very recently, the bicircualr counter-rotating circularly polarized (BCCP) driving laser field has been a hot spot due to its unique features and applications. A BCCP laser field consists of two coplanar counter-rotating circularly polarized fields with the fundamental field frequency and its second harmonic Long (); Eichmann (); DB (). The time-dependent electric field performs a Lissajous figure on its polarization plane with a three-fold symmetry, i.e. it consists of three segments or lobes per cycle. Therefore, the laser field peaks three times one cycle of the fundamental, producing three ionization bursts. The electrons ionized to the continuum near the laser peaks can successfully recombine with the parent ion, radiating an attosecond pulse train with three pulses one cycle of the fundamental LM (); DB2 (); DB3 (); DB4 ().
This attractive laser field has been demonstrated to be an effective tool for fully controlling the high harmonic polarization from linear through elliptical to circular polarization Fleischer (). Importantly, the conversion efficiency for harmonics with arbitrary polarization is comparable to that of the standard HHG process, which yields linearly polarized high harmonics driven by a linearly polarized laser pulse LM (); DB (); Fleischer (). When the BCCP laser pulse is acting on an isotropic medium, in the frequency domain, a comb of high efficiency harmonic with frequency (3n1) is radiated LM (); Kfir (); Reich (); AF (). According to the selection rules Fleischer (); Alon (); Pisanty (); DB5 (); Mauger (), the harmonics with frequency (3n+1) have the same helicity as the fundamental while the (3n-1) harmonics have the same helicity as the second harmonic pulse (3n is parity forbidden). Selecting harmonics with one kind of helicity would generate attosecond pulses with circular polarization. Some efforts have been made for the selection, e.g., optimizing the phase-matching conditions Kfir () or using crossed beams of counter-rotating circularly polarized driving lasers Hickstein (). This circularly polarized light source is very useful for ultrafast spectroscopy, especially for the ultrafast circular dichroism studies AF (); Hergenhahn (); NB (); Powis (). It was also shown that the bicircular HHG driven by this BCCP laser pulse can be a new method to probe dynamical symmetries of atoms and molecules and their evolution in time Baykusheva (); Reich2 (). These attractive features of HHG in BCCP laser fields are related to the symmetry of the laser field and the corresponding three-step dynamics of the continuum electrons. To better understand the features of HHG in the BCCP laser field, the correspondence between the harmonic polarization and the electron trajectories needs to be established.
In our study, the polarization properties of high harmonics generated in 2D model Ne atoms by the BCCP laser fields are investigated. We focus on the helicity reversions in the harmonic spectrum, where the rotation direction of the harmonic electric field reverses at specific orders. To explain this phenomenon, we depict the three-step picture of HHG in BCCP fields and reveal the correspondence between the time-frequency distribution of harmonic degree of circular polarization and the classical trajectories. Our study helps people better understand the HHG process and depicts a clear picture for the HHG process in the BCCP laser field.
Ii Theoretical model
The high harmonic spectrum generated from the interaction between the BCCP laser field and Ne atoms is calculated by numerically solve the two-dimensional time-dependent Schrödinger equation (TDSE) in length-gauge M. (); Rui (). The 2D model potential of Ne atom is given by Barth (); Zhu1 () (atomic units(a.u.) are used throughout this paper):
where and a=2.88172 to obtain the ionization potential of Ne atom =0.793 a.u. for 2p orbitals. denotes the electron position in the two-dimensional x-y plane.
The bicircular electric field is obtained by combining left-circularly polarized and right-circularly polarized laser fields. The laser vector in the polarized plane is defined by LM ()
where denotes the integer multiple of the fundamental frequency. The laser field has laser lobes and reflects fold symmetry. The field with helicity +1 rotates counterclockwise while the second field with helicity -1 rotates clockwise. is the trapezoidal envelop with 2 cycle rising and 2 cycle falling edges and 5 cycle plateau (in units of fundamental). We set the same laser intensity between the field and the field. is the amplitude of the laser field.
We use the split-operator method to solve Eq.(1). To avoid spurious reflections from the spatial boundaries, the electron wave function is multiplied by a “mask function” with the following form Yuan ()
at each time step, in which
For all results reported here, we set the width of the “absorbing” area as =60 a.u. and =60 a.u.. In our calculation, we set the range of the “mask function” from -240 a.u. to 240 a.u. in the Cartesian coordinates, i.e., ==240 a.u.. The space step is 0.1 a.u.. The time step is a.u..
The time-dependent dipole acceleration can be obtained by means of the Ehrenfest theorem
Then, the harmonic spectrum is obtained by Fourier transforming the dipole acceleration
The intensity of left- and right- polarized harmonic components can be obtained by
which can be obtained from the intensity measurements. The ellipticity has a one-to-one correspondence with the DCP as DB (). The two parameters both vary in the interval from -1 to 1 and have the same sign.
The rotation direction of the electric field can be quantitatively described by the helicity, which is defined as:
The helicity takes the values -1 and +1, indicating the two opposite rotation directions.
Iii Results and discussions
Firstly, we present the harmonic intensity and the DCP as a function of the harmonic order with different laser parameters in Fig.1. We choose the example of Ne atoms with the 2p ground state. In this case, the two components of the bicircular field with =2 have equal intensities. According to the selection rules, the harmonics 3n are absent, and the harmonics (3n+1) and (3n-1) have opposite helicities. The harmonic spectra for left-circularly polarized harmonics (3n+1) and right-circularly polarized harmonics (3n-1) are presented with light red lines and blue lines respectively. The envelopes of these harmonic spectra are presented with the dark red solid lines and blue dashed lines correspondingly. For equal intensities of the two laser field components, the position of the harmonic cutoff is obtained by DB ()
The ponderomotive energy of the entire field is defined by . The cutoff obtained by Eq.(11) is almost coincident with that from the simulated results.
The harmonic spectrum driven by the bicircular laser field with the fundamental laser wavelength nm and laser intensity W/cm is shown in Fig. 1(a). From this figure, one can see that, the intensities of the harmonics (3n-1), having the helicity , are different from that of the harmonics (3n+1), for which . There are two helicity reversions where the red-solid and blue-dashed lines intersect. One is at about the ionization threshold (22nd-order harmonic), the other is at about 107th-order harmonic (labelled by the dash-dot line). Below the ionization threshold, the intensities of (3n+1)th-order harmonics are higher than those of (3n-1)th-order. Between the ionization threshold 22nd-order and the 107th-order harmonic, the intensities of (3n-1)th-order harmonics turn to be higher than those of the (3n+1)th-order. While between 107th-order and the cutoff 120th-order obtained by the Eq.(11), an opposite result to that between the ionization threshold 22nd-order and the 107th-order harmonic order appears. That is, the intensities of the harmonics (3n+1)th-order co-rotating with the field are higher than those of the (3n-1)th-order counter-rotating with the field. In order to clearly show the helicity reversion phenomenon in the harmonic spectra, we display the DCP in Fig.1(e). According to the helicity definition, the helicity reversion appears at the position where the DCP changes its sign. From Fig.1(e), The DCP crosses zero at ionization threshold 22nd-order and 107th-order harmonic, which is coincident with the discussion above.
To check whether the phenomenon is general for the bicircular field interaction with Ne atoms and uncover the reason for the reversion, we apply different laser parameters for fundamental laser wavelength nm , 1600 nm and laser intensity W/cm, W/cm. The harmonic spectra are presented in Figs.1(b-d) and the DCP are plotted in Figs.1(f-h). It is shown that the intensities of (3n+1)th-order harmonics are higher than those of the (3n-1)th-order harmonics below the ionization threshold, and correspondingly, the values of the DCP are negative. There is a helicity reversion at the ionization threshold for different laser parameters. We have also calculated the harmonic spectra for other model atoms with different two-dimensional soft-core potentials and different ionization potentials. It is shown that the reversion at the threshold is not a general phenomenon and is dependent on the target chosen. Except the helicity reversion at the ionization threshold, the helicity reversion above the threshold is found in all these spectra. The reversion position is labelled by the dash-dot line. In our study, we mainly investigate the helicity reversions in the above-threshold region.
In order to explain the phenomenon of the helicity reversion and elucidate the physics behind it, we present the time-frequency distribution of the harmonic DCP and the classical ionization and recombination paths in Fig.2. The laser parameters are the same as in Fig.1(a). Here we analyze the time-frequency properties of the harmonic DCP with the Gabor transform, which has been proved to be a very powerful tool to analyze the emission times of HHG and to clearly discriminate the electron trajectories. Considering the dipole acceleration of Eq.(6), the Gabor transform is defined as CC ()
where is the frequency of the high harmonics. In our study, we use . The DCP distribution is obtained by , in which . In Fig.2(a), we present the DCP distribution of harmonics between 1.5 optical cycle and 8 optical cycle. Due to the symmetry of the BCCP laser pulse, three radiation bursts per fundamental laser pulse cycle can be found in this figure, as predicted in previous works LM (); DB2 (); DB3 (); DB4 (). There are two helicity reversions (labelled by the sold line and the dashed line respectively) and three distinct regions can be identified: (I) Below threshold 22nd-order (labelled by the solid line), the DCP is mostly negative. The harmonic emissions counter-rotating with the laser field are radiated. (II) The middle region between the threshold (22nd-order) and the second reversion (about 107th-order, labelled by the dashed line), the DCP is high and positive. The harmonic emissions are emitted with the left-circular polarization. (III) The upper region between the second reversion (about 107th-order) and the cutoff (120th-order), the DCP is almost negative opposed to the second region. In the second region, the highest DCP is at about the 30th-order harmonic, corresponding to the position of the maximum intensity difference between the left-circularly and right-circularly polarized harmonic components in Fig.1(a). Then the pulse DCP decreases gradually along with the increase of the harmonic order. At about 107th-order denoted by the dashed line, the DCP transits from positive value to negative value. These are directly reflected in Fig.1(e).
To give a clear interpretation of the quantum results above, we explore the classical behavior of the electron path by solving Newton’s equation under the BCCP laser field driven Meng (); Zhuo (); Yin (); Tong (); Min (). In Fig.2(b), we present the electrons classical ionization times labelled by the blue diamonds, as well as, the recombination times labelled by the red pentagrams. The solid markers are used for the short trajectories and the hollow markers are used for the long trajectories. Note that HHG in bicircular fields is dominated by short trajectories DB (). We focus on the short trajectories of the electrons ionized at one laser field lobe. As shown in this figure, three ionization bursts and three recombination bursts in one fundamental laser cycle reflect the threefold symmetry of the counter-rotating circularly polarized laser field Long (); Eichmann (); DB (). The electrons ionized in one laser lobe can recombine with the parent ion in the following two lobes and two peaks marked as and are generated respectively at laser lobe and (see Fig.3(b)). The maxima of the peak and peak are 111th-order and 122nd-order harmonics respectively, which approximately correspond to the reversion position and the cutoff position. The slight deviation between the quantum and classical model originates from the quantum effect and the absence of harmonic 3n in quantum model. But it does not influence greatly on our conclusion. By comparing the time-frequency analysis and the classical calculations, one can see that, the DCP between the ionization threshold and about 40th-order harmonic is positive, which is only contributed by . However, between the maximum of and the harmonic cutoff, the DCP contributed by only is negative. Besides, at the overlapping part of the two peaks and (40th-order harmonic to 111th-order harmonic), the DCP decreases along with the increase of harmonic order. The correspondence relationship between the time-frequency analysis and the classical calculations indicates that the harmonics generated by are left-circularly polarized while the harmonics generated by are right-circularly polarized. Therefore, the harmonics, which are contributed by the electrons ionized at one field lobe and recombining at the following adjacent lobes, have opposite helicities.
To further illustrate our idea, we present the electron classical trajectories and the plot of the laser field vector of the field in Fig.3. The dash-dot green line in Fig.3(a) is the amplitude of the laser field with three peaks per optical cycle. Every peak corresponds to one laser lobe. In order to increase the readability of the figure, the second lobe marked as and the third lobe marked as are superimposed by the light red background and light blue background respectively. We calculate the short trajectories of the electrons ionized at first laser lobe. It is found that, these trajectories can be classified by the laser lobe they recombine at. That is, the trajectories of the electrons recombining at are classified as group (the red lines), and the trajectories of the electrons recombining at are classified as group (the blue dot lines). To give a clear description of the trajectory characteristics, the polar plot of the electric field vector of the - field is shown in Fig.3(b). A clover-leaf structure that has a threefold symmetry is presented. The laser field symmetry is reflected by an attosecond pulse train with three pulses generated in one optical cycle in Fig.2. The electronic field vector at the ionization time at the first laser lobe is marked as bold purple line . The bold red line at and bold blue line at denote the laser field vector at the recombination times for electrons of groups and respectively. The electrons ionized at can recombine with the nucleus at and through the path and . Note that the trajectories and correspondingly contribute to short trajectories of and in Fig.2(b). Then, based on the discussion about Fig.2, the harmonics radiated from and have opposite helicities.
Next, we also calculate the relative phase of the , harmonic components using the strong-field approximation (SFA) model Zhou (); Lewenstein (). Fig.4(a) and (b) respectively exhibit the calculation results of the harmonic emissions generated by electrons recombining at and . Based on the classical calculations in Fig.2(b), we choose the recombination time period and integral over a ionized time period for each recombination time chosen. From this figure, one can see that the relative phase of the two components radiated at the laser lobe varies between 0.4 to , which indicates that the radiated harmonics have the positive DCP. However, the relative phase of the two components radiated at the lobe changes from 0 to -0.5 which indicates that the harmonics have the negative DCP. The result of Fig.4 coincides with that of Fig.2. Therefore, our conclusion that the harmonics contributed by electrons ionized at the first lobe and recombining at the following two lobes have opposite helicities is confirmed by using the classical model and the SFA model.
Finally, we have also calculated the HHG driven by BCCP laser field with . The fundamental laser wavelength and laser intensity are nm and W/cm, respectively. In Fig.5(a), the harmonic spectra for left-circularly polarized harmonics (4n+1) and right-circularly polarized harmonics (4n-1) are presented with light red and blue lines respectively. The envelopes of these harmonic spectra are presented with the dark red solid lines and blue dashed lines correspondingly. In this spectrum, four helicity reversions are obviously observed. The first one is also at the ionization threshold, which is determined by the target. The remaining three reversions are at about 69th-order (the dashed line), 80th-order and 95th-order (the dash-dot line). In order to clearly display the phenomenon of helicity reversion of the harmonics, we plot the time-frequency distribution of the harmonic DCP. Four attosecond pulses are generated per optical cycle, which reflects the fourfold symmetry of the laser field. The helicity reversions at about 69th-order and 95th-order are labelled by the dashed line and the dash-dot line. Similarly, the classical electron paths in one optical cycle are exhibited in Fig.5(c). We only study the trajectories of the electrons ionized at one laser lobe and recombining at the following three lobes. The recombination times of these electrons are labelled by the dark red pentagrams. Others are labelled by the light red pentagrams. Three peaks and are presented in Fig.5(c). Comparing Fig.5(b) with Fig.5(c), it is shown that, the DCP in the region between the ionization threshold 22nd-order and the reversion position 69th-order is positive and decreases with harmonic order. This region is dominated by . The region between the reversion position 95th-order (the maximum of peak ) and the harmonic cutoff 109th-order (the maximum of peak ) is only contributed by . The helicity of is negative, which is opposite to that of . The helicity reversions at the 69th-order and the 95th-order in Figs.5(a) and 5(b) appear at the maxima of and in Fig.5(c). In the overlapping region of and between 69th-order and 95th-order, the harmonic emissions have both positive and negative DCP and an additional helicity reversion appears. This will be discussed later. For a clear analysis of the classical electron trajectories in the laser field, the laser field vector is plotted in Fig.5(d). The electronic field vector at the ionization time at the first laser lobe is marked as bold purple line . The recombination times of the short trajectories of , and are marked as the red bold line in laser lobe , cyan bold line in lobe and blue bold line in lobe respectively. It is found that the peaks , and are contributed by the electrons recombining at different laser lobes. Based on the above comparison between Fig.5(b) and Fig.5(c), the harmonic emissions generated in have opposite helicities with those in , and the helicities of harmonics generated in are opposite to those in , which is consistent with the conclusion discussed above for laser field.
To understand the appearance of the additional helicity reversion in the overlapping region of and , the weights of the contributing electron trajectories should be obtained. Fig. 6 shows the weights of the short trajectories in both the and laser fields. For each trajectory, the weight is calculated by . is the ionization rate at time obtained from the TDSE simulation by calculating the derivative of the time-dependent ionization probability. is the weight with a Gaussian distribution for the tunnel electron with transverse initial velocity Delone (). For the laser field as shown in Fig. 6(a), the spectrum range can be divided into four regions. In the first region below the 69th-order, the weight of peak is much greater than those of peaks and . In the second region from the 69th-order to the 79th-order, peak disappears and the weight of is greater than . Then from about 79th-order to 95th-order, the weight of becomes greater than in the third region. In the last region above the 95th-order, only peak remains contributing to the harmonic emissions. In each region, the helicity of the generated harmonics is determined by helicity of the strongest contributing peak, which can well explain the results in Figs. 5(a-b). Owing to the intersection between the weights of peaks and , the additional helicity reversion in the overlapping region of and appears. For the laser field as shown in Fig. 6(b), the spectrum range can be divided by the 111th-order into two regions. The weight of peak is much greater than that of peak in the first region. In the second region, only peak contributes to the harmonic emissions. Since the weights of the two peaks have no intersection, there is no extra helicity reversion in the BCCP laser field.
Based on all the discussions above, a picture of the electron classical trajectories in HHG driven by BCCP laser field is revealed. Besides, the correspondence between the harmonic polarization and the classical trajectories has also been established. Specifically, the electron ionized at one laser lobe can recombine with the nucleus at the following laser lobes. These recombining trajectories can be classified into different groups based on the laser lobe they recombining at. Moreover, the harmonics emitted from these groups at adjacent laser lobes have opposite helicities. For instance, in laser field, groups and are both ionized at one laser lobe . And they recombine with the nucleus at the following laser lobes and and form the peaks and , respectively. The harmonics from and are radiated with opposite helicities and the helicity reversion occurs at the maximum of . In laser field, electrons ionized at one laser lobe recombine at the following lobes , , and form the peaks , , , respectively. The harmonics from peaks and have positive helicity while those from has negative helicity. The helicity reversion occurs at the position where the dominantly contributing peak is changed to its neighbouring peak, i.e. at the maxima of and and at the intersection between the weights of and .
In conclusion, under the interaction of bicircular laser field with the 2D model Ne atoms, a harmonic spectrum with different helicities in different spectrum ranges is generated. It is shown that, the helicity of the elliptically polarized harmonic emission is reversed at particular harmonic orders. Based on the time-frequency analysis and the classical three-step model, the relationship between the reversion position and the electron trajectories is obtained. The electrons ionized at one lobe of laser field can recombine with the nucleus at the following laser lobes radiating harmonics. These harmonic emissions generated from adjacent laser lobes have opposite helicities. Besides, we confirm our conclusion through the calculation of relative phase by the SFA model. Our study performs a detailed analysis for the HHG driven by the BCCP laser field in terms of the classical trajectories and established the correspondence between the harmonic polarization and the classical trajectories.
This work was supported by the National Natural Science Foundation of China under Grants No. 11404123, 11234004, 61275126, 11422435 and 11574101. 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).
- (1) P. B. Corkum and F. Krausz, Nat. Phys. 3, 381–387 (2007).
- (2) F. Krausz and M. Ivanov, Rev. Mod. Phys. 81, 163–234 (2009).
- (3) X. Zhu, X. Liu, Y. Li, M. Qin, Q. Zhang, P. Lan, and P. Lu, Phys. Rev. A 91, 043418 (2015).
- (4) L. He, P. Lan, Q. Zhang, C. Zhai, F. Wang, W. Shi, and P. Lu, Phys. Rev. A 92, 043403 (2015).
- (5) X. Zhu, M. Qin, Q. Zhang, W. Hong, Z. Xu, and P. Lu, Opt. Express 20, 16275 (2012).
- (6) 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, Science 314, 443¨C446 (2006).
- (7) L. Medišauskas, J. Wragg, H. van der Hart, and M. Y. Ivanov, Phys. Rev. Lett. 115, 153001 (2015).
- (8) A. L. Cavalieri, N. Müller, T. Uphues, V. S. Yakovlev, A. Baltuška, B. Horvath, B. Schmidt, L. Blümel, R. Holzwarth, S. Hendel, M. Drescher, U. Kleineberg, P. M. Echenique, R. Kienberger, F. Krausz, and U. Heinzmann, Nature (London) 449, 1029 (2007).
- (9) J. Itatani, J. Levesque, D. Zeidler, H. Niikura, H. P¨¦pin, J. C. Kieffer, P. B. Corkum, and D. M. Villeneuve, Nature 432, 867¨C871 (2004).
- (10) S. Haessler, J. Caillat, W. Boutu, C. Giovanetti-Teixeira, T. Ruchon, T. Auguste, Z. Diveki, P. Breger, A. Maquet, B. Carré, R. Taïeb, and P. Salières, Nat. Phys. 6, 200 (2010).
- (11) P. M. Kraus et al., Science 350, 790 (2015).
- (12) C. Zhai, L. He, P. Lan, X. Zhu, Y. Li, F. Wang, W. Shi, Q. Zhang, and P. Lu, Scientific Reports 6, 23236 (2016); M. Qin, and X. Zhu, Optics and Laser Technology 87, 79-86 (2017).
- (13) O. Smirnova, Y. Mairesse, S. Patchkovskii, N. Dudovich, D. Villeneuve, P. Corkum, and M. Y. Ivanov, Nature (London) 460, 972 (2009).
- (14) S. Baker, J. Robinson, C. Haworth, H. Teng, R. A. Smith, C. C. Chiril, M. Lein, J.W. G. Tisch, and J. P. Marangos, Science 312, 424 (2006).
- (15) H. J. Wörner, J. B. Bertrand, D. V. Kartashov, P. B. Corkum, and D. M. Villeneuve, Nature (London) 466, 604 (2010).
- (16) O. Travnikova, J.-C. Liu, A. Lindblad, C. Nicolas, J. Söderström, V. Kimberg, F. Gelmukhanov, and C. Miron, Phys. Rev. Lett. 105, 233001 (2010).
- (17) C. Boeglin, E. Beaurepaire, V. Halt¨¦, V. L¨®pez-Flores, C. Stamm, N. Pontius, H. A. D¨¹rr, and J.-Y. Bigot, Nature (London) 465, 458 (2010).
- (18) P. B. Corkum, Phys. Rev. Lett. 71, 1994¨C1997 (1993).
- (19) K. J. Schafer, B.Yang, L. F. DiMauro, and K. C.Kulander, Phys. Rev. Lett. 70, 1599 (1993).
- (20) S. Long, W. Becker, and J. K. McIver, Phys. Rev. A 52, 2262 (1995).
- (21) H. Eichmann, A. Egbert, S. Nolte, C. Momma, B. Wellegehausen, W. Becker, S. Long, and J. K. McIver, Phys. Rev. A 51, R3414 (1995).
- (22) D. B. Milošević, W. Becker, and R. Kopold, Phys. Rev. A 61, 063403 (2000);
- (23) D. B. Milošević, Opt. Lett. 40, 2381 (2015).
- (24) D. B. Milošević, Phys. Rev. A 92, 043827 (2015).
- (25) D. B. Milošević and W. Becker, Phys. Rev. A 62, 011403 (2000).
- (26) A. Fleischer, O. Kfir, T. Diskin, P. Sidorenko, and O. Cohen, Nat. Photonics 8, 543 (2014).
- (27) O. Kfir, P. Grychtol, E. Turgut, R. Knut, D. Zusin, D. Popmintchev, T. Popmintchev, H. Nembach, J. M. Shaw, A. Fleischer, H. Kapteyn, M. Murnane, and O. Cohen, Nat. Photonics 9, 99 (2015).
- (28) D. M. Reich and L. B. Madsen, Phys. Rev. A 93, 043411 (2016).
- (29) A. Ferré, C. Handschin, M. Dumergue, F. Burgy, A. Comby, D. Descamps, B. Fabre, G. A. Garcia, R. Géneaux, L. Merceron, E. Mével, L. Nahon, S. Petit, B. Pons, D. Staedter, S.Weber, T. Ruchon, V. Blanchet, and Y. Mairesse, Nat. Photonics 9, 93 (2014).
- (30) O. E. Alon, V. Averbukh, and N. Moiseyev, Phys. Rev. Lett. 80, 3743¨C3746 (1998).
- (31) E. Pisanty, S. Sukiasyan, and M. Ivanov, Phys. Rev. A 90, 043829 (2014).
- (32) D. B. Milošević, J. Phys. B 48, 171001 (2015).
- (33) F. Mauger, A. D. Bandrauk and T. Uzer, J. Phys. B: At. Mol. Opt. Phys. 49, 10LT01 (5pp) (2016).
- (34) D. D. Hickstein, F. J. Dollar, P. Grychtol, J. L. Ellis, R. Knut, C. Hernández-García, D. Zusin, C. Gentry, J. M. Shaw, T. Fan, K. M. Dorney, A. Becker, A. Jaroń-Becker, H. C. Kapteyn, M. M. Murnane, and C. G. Durfee, Nat. Photonics 9, 743-750 (2015).
- (35) U. Hergenhahn, E. E. Rennie, O. Kugeler, S. Marburger, T. Lischke, I. Powis, and G. Garcia, J. Chem. Phys. 120, 4553 (2004).
- (36) N. Böwering, T. Lischke, B. Schmidtke, N. Müller, T. Khalil, and U. Heinzmann, Phys. Rev. Lett. 86, 1187 (2001).
- (37) I. Powis, J. Chem. Phys. 112, 301 (2000).
- (38) D. Baykusheva, M. S. Ahsan, N. Lin, and H. J. Wörner, Phys. Rev. Lett. 116, 123001 (2016).
- (39) D. M. Reich and L. B. Madsen, arXiv:1606.05524v1.
- (40) M. D. Feit, J. A. Fleck, JR and A. Steiger, J. Comput. Phys. 47, 412-433 (1982).
- (41) M. He, Y. Li, Y. Zhou, M. Li, and P. Lu, Phys. Rev. A 93, 033406 (2016).
- (42) I. Barth and M. Lein, J. Phys. B 47, 204016 (2014).
- (43) X. Zhu, P. Lan, K. Liu, Y. Li, X. Liu, Q. Zhang, I. Barth, and P. Lu, Opt. Express 24, 4196 (2016).
- (44) K. J. Yuan and A. D. Bandrauk, Phys. Rev. A 83, 063422 (2011).
- (45) P. Antoine, B. Carre, A. L’Huillier, and M. Lewenstein, Phys. Rev. A 55, 1314 (1997).
- (46) E. Hecht, Optics Fourth Edition.
- (47) M. Born and E. Wolf, Principles of the Optics.
- (48) S. Odžak, and D. B. Milošević, Phys. Rev. A 92, 053416 (2015).
- (49) C. C. Chirilǎ, I. Dreissigacker, E. V. van der Zwan, and M. Lein, Phys. Rev. A 81, 033412 (2010).
- (50) X. Ma, Y. Zhou, and P. Lu, Phys. Rev. A 93, 013425 (2016).
- (51) Z. Wang, M. Li, Y. Zhou, Y. Li, P. Lan, and P. Lu, Phys. Rev. A 93, 013418 (2016).
- (52) Y. Chen, Y. Zhou, Y. Li, M. Li, P. Lan, and P. Lu, The Journal of Chemical Physics 144 024304 (2016).
- (53) A. Tong, Y. Zhou, and P. Lu, Opt. Express 12, 15774 (2015).
- (54) M. Li, P. Zhang, S. Luo, Y. Zhou, Q. Zhang, P. Lan, and P. Lu, Phys. Rev. A 92, 063404 (2015).
- (55) Y. Zhou, O.I. Tolstikhin, and T. Morishita, Phys. Rev. Lett. 116, 173001 (2016).
- (56) M. Lewenstein, Ph. Balcou, M. Yu. Ivanov, A. L’Huillier, and P. B. Corkum, Phys. Rev. A 49, 2117 (1994).
- (57) N. B. Delone, V. P. Krainov, J. Opt. Soc. Am. B 8, 1207 (1991).