==== Front Sci Rep Sci Rep Scientific Reports 2045-2322 Nature Publishing Group UK London 79034 10.1038/s41598-020-79034-2 Article Attosecond streaking using a rescattered electron in an intense laser field Kim Yang Hwan 12 Ivanov Igor A. 1 Hwang Sung In 1 Kim Kyungseung 1 Nam Chang Hee 12 Kim Kyung Taec kyungtaec@gist.ac.kr 12 1 grid.410720.00000 0004 1784 4496Center for Relativistic Laser Science, Institute for Basic Science, Gwangju, 61005 Korea 2 grid.61221.360000 0001 1033 9831Department of Physics and Photon Science, Gwangju Institute of Science and Technology, Gwangju, 61005 Korea 16 12 2020 16 12 2020 2020 10 2207523 9 2020 3 12 2020 © The Author(s) 2020Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.When an atom or molecule is exposed to a strong laser field, an electron can tunnel out from the parent ion and moves along a specific trajectory. This ultrafast electron motion is sensitive to a variation of the laser field. Thus, it can be used as a fast temporal gate for the temporal characterization of the laser field. Here, we demonstrate a new type of attosecond streaking wherein a rescattered electron trajectory is manipulated by an ultrashort laser pulse. The vector potential of the laser pulse is directly recorded in the photoelectron spectra of the rescattered electron. In contrast to high harmonic generation methods, our approach has no directional ambiguity in space, leading to complete in situ temporal characterization. In addition, it provides timing information on ionization and re-scattering events. Therefore, our approach can be a useful tool for the investigation of strong-field processes triggered by rescattering, such as non-sequential double ionization and laser-induced electron diffraction. Subject terms Ultrafast photonicsNonlinear opticsInstitute for Basic ScienceIBS-R012-D1IBS-R012-D1IBS-R012-D1IBS-R012-D1IBS-R012-D1Kim Yang Hwan Hwang Sung In Kim Kyungseung Nam Chang Hee Kim Kyung Taec issue-copyright-statement© The Author(s) 2020 ==== Body Introduction The temporal characterization of a laser field is of great interest in many laser applications. There are many temporal characterization techniques that can be applied for femtosecond laser pulses such as FROG and SPIDER1–3. They utilize a response from a nonlinear crystal for the temporal characterization of a laser pulse. The complete characterization of the laser field including the carrier-envelope phase (CEP) has also been demonstrated using the attosecond streaking technique4,5. Recently, it was also shown that the laser pulse can be characterized using ionization in the air6. These pulse characterization techniques are convenient for use in many applications. However, as they require particular nonlinear interactions, they are in general implemented independently. For the temporal characterization of ultrashort laser pulses used for laser–matter interactions, it is ideal to measure a laser pulse at the place where it is being used. A few approaches support such in situ measurements utilizing the ultrafast electron dynamics. They utilize the process of high-harmonic generation (HHG) as a fast temporal gate7–11, which is particularly useful for the characterization of the laser field used in high harmonic spectroscopy. However, there is ambiguity of the direction of the reconstructed laser field since the relative phase between two laser pulses is measured. The direction of the laser field in space cannot be determined by these measurements. In this work, we demonstrate a new attosecond streaking technique in which a laser field is completely determined by manipulating the trajectory of a backward rescattered electron in an intense laser field. Since the electron dynamics in the process of above-threshold ionization (ATI) is utilized, it can be applied for a broad intensity range of 1013∼1015 W/cm2. The direction of the laser field can be completely determined without ambiguity thanks to the stereo-ATI time-of-flight spectrometer used for the photoelectron measurement. In addition, our new streaking technique provides timing information on ionization and backward rescattering events. Therefore, it will become a useful tool for the investigation of physical processes related to rescattering dynamics such as such as HHG12, ATI13,14, non-sequential double ionization15–17, and laser-induced electron diffraction18–22. Results We used the semi-classical model based on strong-field approximation12 to describe the dynamics of an electron in an intense laser field. In the semi-classical model, an electron moves along different classical trajectories depending on its ionization time23. We focus on a backward rescattered electron that gains the maximum energy (called the cutoff electron hereafter). The cutoff electron tunnels out near the local extrema of the electric field and is driven back to the parent ion as the direction of the laser field is reversed, as shown in Fig. 1a. Then, the electron is rescattered backward near the next zero-crossing time of the laser field.Figure 1 Manipulation of a backward rescattered electron trajectory in above-threshold ionization. (a) Classical electron trajectories that yield maximum energy along one direction in an intense laser field with (orange dashed arrow) or without (blue solid arrow) a perturbing laser field. The main laser field and the perturbing laser field are indicated by the red solid line and red dashed line, respectively. The magnified plot from the ionization time (blue cross) to the backward rescattering time (blue circle) is shown in the inset. (b) ATI spectra along the polarization direction of the laser field obtained at zero (blue solid line) and half-cycle (0.5 T, red solid line) delays between the main laser field and the perturbing laser field. (c) ATI spectra along the polarization of the laser field obtained by solving the time-dependent Schrodinger equation in 3D with a model potential of Xe using the main pulse with a peak intensity of 1×1014 W/cm2 and pulse duration of 4.4 fs together with the perturbing pulse, which has the same temporal waveform as the main field with a peak intensity of 5.0×1011 W/cm2. The CEPs of both pulses were 0.5π. (d) The cutoff momentum shift Δpτ obtained by fitting the Airy function to the spectra shown in (c) (blue solid line) and Δpτ obtained by solving the semi-classical model ignoring Coulomb field (red solid line). (e) Spectral response of the streaking method estimated by solving TDSE calculations. The vertical dashed line indicates the center angular frequency ω0 of a main pulse. A 730 nm, 5 fs laser pulse with the intensity of 2.0×1014 W/cm2 was used for the main laser pulse. A single-cycle laser pulse with the center frequency of 8 ω0 was used as a perturbing field. The trajectory of the cutoff electron can be manipulated by adding a weak perturbing laser field24,25. The intensity of the perturbing laser field is significantly weaker than that of the main pulse (see the Methods section). The addition of the perturbing laser field changes the final momentum of the cutoff electron by Δp as shown in Fig. 1a. The cutoff momentum shift Δp can be found as a function of the time delay τ between the main and perturbing laser pulses using the semi-classical model (see the Methods section for derivation), 1 Δpτ=-2Aptr-τ+1tr-ti∫titrApt-τdt, where Ap is vector potential of the perturbing laser field, and ti and tr are ionization time and backward rescattering time in the absence of the perturbing field, respectively. The first term on the right-hand side of Eq. (1) is proportional to the vector potential of the perturbing field at the rescattering time and the second term is proportional to the time-averaged value of the vector potential of the perturbing field between the ionization time and backward rescattering time. Therefore, the vector potential of the perturbing laser field, ionization time, and backward rescattering time can be found in the momentum shift of the cutoff electron Δp(τ) (see the Methods section for details). The cutoff momentum shift Δp(τ) of a backward rescattered electron can be found by analyzing the above-threshold ionization spectrum. While the cutoff energy of an electron directly liberated or rescattered forward reaches only up to 2Up, the cutoff energy of the backward rescattered electron reaches up to approximately 10Up as shown in Fig. 1b. Here, Up is ponderomotive energy of an electron in the laser field. Therefore, the spectrum of the backward rescattered electron is clearly distinguished from that of the direct or forward rescattered electrons. However, the cutoff momentum shift Δp(τ) should be estimated carefully because the perturbing field alters the momentum and tunneling ionization rate at the same time. For example, photoelectron spectra obtained at two time delays by solving the time-dependent Schrödinger equation (TDSE) in 3D for a Xe potential exhibit differences in the amplitude and cutoff momentum as shown in Fig. 1b,c. Therefore, the cutoff momentum shift cannot be simply determined, and both the cutoff momentum shift and the amplitude modulation should be taken into account. It is known that the momentum distribution of a backward rescattered electron along the polarization of the laser field in the cutoff energy region is well described by the Airy function26. Therefore, we fit the Airy function to the ATI spectra at each delay to find the cutoff momentum shift. To validate our approach, we compared the cutoff momentum shift Δp(τ) obtained by fitting the Airy function to the ATI spectra shown in the upper panel of Fig. 1c, with the cutoff momentum shift Δp(τ) obtained by solving the semi-classical model, as shown in Fig. 1d. The cutoff momentum shift induced by the perturbing field was accurately retrieved by fitting the Airy function to the ATI spectra. The cutoff momentum shift in Eq. (1) can be averaged out if the frequency of the perturbing laser field is too high. This limits the frequency of the perturbing laser field that can be measured. This spectral response was estimated by solving TDSE using a single-cycle perturbing laser field whose beandwidth is broader than 10ω0. The spectral amplitude and phase of the cutoff modulation show a broad distribution up to 10ω0 as shown in Fig. 1e where ω0 is the center angular frequency of the main pulse. It means that our streaking method can be applied for a broad spectral range up to 10ω0. The spectral phase modulates within ±0.2π rad only, indicating that the cutoff modulation is similar to the temporal profile of the perturbing laser field. This spectral modulation was reproduced using Eq. (1). Therefore, the temporal profile of the perturbing laser field can be accurately reconstructed using Eq. (1). For the experimental demonstration, we measured ATI spectra using a homebuilt stereo-ATI time-of-flight spectrometer27 as shown in Fig. 2a (see the Methods section). The pulse was split into two pulses (denoted by M: main pulse, and by P: perturbing pulse, as shown in Fig. 2a) and combined with a time delay τ using a mirror with a round hole at the center in an interferometer without any dispersive optical elements. In this manner, the two pulses have the same waveform while the peak intensity of the perturbing pulse was only 0.5% of the main pulse at the focus. ATI spectra for both left and right directions along the polarization direction were measured using the stereo-ATI spectrometer. Slits were placed in the stereo-ATI spectrometer to obtain an energy spectrum with a small opening angle (2°) along the polarization of the laser field. The target cylinder (not shown in Fig. 2) of the spectrometer was filled with Xe gas with a pressure of 2.0×10-4 torr. We scanned the insertion of the wedge pair before we scan the delay. Corresponding stereo-ATI spectra and asymmetry plot are shown in Fig. 2b,c. We chose specific CEPs of -0.46π and 0.54π for delay scans so that the continuum at the cutoff is observed for both directions. The continuum spectrum indicates that the cutoff spectrum is contributed by electrons tunneled out from a single half optical cycle of the laser field. Therefore, Eq. (1) can be directly applied for the cutoff momentum shift.Figure 2 Schematic drawing of the experimental setup and CEP dependence of the momentum distribution. (a) Schematic drawing of the experimental setup. The input laser pulse is split into a main pulse (denoted by M) and a perturbing pulse (denoted by P) using a mirror with a hole at the center. The two laser pulses are combined after controlling the time delay (see the Methods section). (b) Stereo-ATI spectra recorded using the main laser pulse with a duration of 5 fs and a peak intensity of 1×1014 W/cm2 in Xe for different CEPs. The negative sign in energy is for denoting direction (left) opposite to the positive sign (right). (c) Left–right asymmetry obtained from (b). The ATI spectra measured along the right (p||>0) and left (p||<0) directions are shown in the upper and lower panels, respectively, in Fig. 3a,b. They clearly show the modulation of the cutoff momentum as a function of time delay. The modulations in Fig. 3b are consistent with a numerical result shown in Fig. 3c obtained by solving the time-dependent Schrödinger equation in 1D for CEP equal to 0.54π. The cutoff momentum shifts Δpτ, which are denoted by white solid lines in Fig. 3a–c, were obtained by fitting the Airy function to the experimental data. A prominent feature is the half-cycle shift of the cutoff modulation observed along the two opposite directions and for two opposite CEPs as shown in Fig. 3a–c. The maximum amplitudes of the cutoff modulation are observed near the delay of 1.1 fs (black arrows) in the upper panel of Fig. 3a and the lower panels of Fig. 3b,c, while they are observed around -0.3 fs (black arrows) as shown in the lower panel of Fig. 3a and in the upper panels of Fig. 3b,c. The half-cycle shift of the cutoff modulation confirms that the electron liberated through tunneling at the single half-optical cycle contributes to the cutoff modulation along one direction, while the electron liberated at the next half-optical cycle contributes to the cutoff modulation along the opposite direction. These measurements confirm that the trajectories of the cutoff electrons liberated in a half-optical cycle of the main laser pulse can be successfully steered using the perturbing laser field.Figure 3 Complete characterization of the perturbing field using the ATI streaking method for different CEPs. (a–c) ATI streaking spectra for (a) experimental results with CEP equal to -0.46π, (b) 0.54π, and (c) numerical results obtained by solving TDSE with the 5.2-fs pulse with CEP equal to 0.54π, and a peak intensity of 1 × 10–14 W/cm2. ATI spectra in the focal volume were averaged for the numerical results. (d)–(f) Reconstructed temporal waveforms of vector potentials of perturbing fields using ATI spectra shown in each column. (g) The CEP of the reconstructed vector potential of the perturbing fields for different fused silica thickness. The CEP measured in the experiment (orange circles), and the theoretical estimation for a fused silica material (blue solid line) are shown. The waveform of the perturbing laser field can be reconstructed from the cutoff momentum shift Δpτ using Eq. (1) (see the Methods section). The reconstructed waveforms obtained for both directions are consistent with each other as shown in Fig. 3d,f. The absolute phase of the waveforms is opposite for the two opposite CEPs as shown in Fig. 3d,e. We assigned the positive sign to the direction of the right arm of the stereo-ATI spectrometer and hence the positive (or negative) sign of the field corresponds to the field pointing toward the right (or left). To validate the accuracy of the reconstruction method, we reconstructed waveforms from the ATI spectra obtained numerically. The reconstructed waveforms for both directions (solid lines) are consistent with the original vector potential used in the numerical calculations as shown in Fig. 3f. In addition, we measured the waveform for different CEPs as shown in Fig. 3g. The CEP of the laser pulses was controlled by the wedge insertion thickness. The CEP shift expected from the theoretical values for a fused silica material is also shown in Fig. 3g. The slope of the two values, experimental and theoretical values, agrees well, showing the capability of the CEP measurement. Note that the absolute phase in Fig. 3g does not have the directional ambiguity in space, because the direction of the laser field can be fixed from the asymmetry measurement shown in Fig. 2. It is important to see if our approach is sensitive to the chirp of the laser pulse. For chirped laser pulse measurements, the main pulse should be fixed to the chirp free condition in which electrons liberated in a single half optical cycle contributes to the cutoff. The chirp of the perturbing laser field was controlled by a pair of glass wedges. Firstly, we performed a measurement with a chirp free perturbing laser pulse. Then, the ionization and rescattering times obtained for the chirp free perturbing laser pulse were used for the reconstruction of the chirped perturbing laser field (See the Methods section). The results of our measurements are summarized in Fig. 4. The reconstructed temporal waveforms are shown in Fig. 4d–f with corresponding ATI spectra shown in the upper panels. The spectral phases of the reconstructed pulses nicely exhibit the chirp structure of the positively chirped (Fig. 4a,d), chirp-free (Fig. 4b,e), and negatively chirped pulses (Fig. 4c,f) reconstructed for both the left (solid lines) and right (dashed lines) directions. The waveform measured along the two directions matches to each other even for their CEPs. We changed the fused silica insertion by ± 0.6 mm, and the reconstructed group delay dispersions (GDD) for positively chirped, chirp-free, and negatively chirped pulses were + 23.8 fs2, − 1.77 fs2, and − 28.1 fs2, respectively. The magnitude of the GDD difference between the positively chirped pulse and the chirp-free pulse is 25.6 fs2, and that between the chirp-free and the negatively chirped pulse is 26.3 fs2. The differences in GDD agree well with the expected GDD change 25.7 fs2 which can be calculated by group-velocity dispersion of fused silica. These results confirm the accuracy of our new streaking technique.Figure 4 Spectra and temporal waveforms reconstructed by using the streaking method for different dispersions. (a)–(c) Spectra and absolute spectral phases for (a) positively chirped, (b) chirp free, and (c) negatively chirped pulses reconstructed by using ATI spectra recorded at the left (solid lines) and right (dashed lines). (d)–(f) Cutoff modulation of ATI spectra (upper panels for each subfigures) measured at the left arm of the stereo-ATI spectrometer and reconstructed temporal waveforms of (d) positively chirped, (e) chirp-free, and (f) negatively chirped pulses. The ATI spectra were measured as a function of delay, but it is shown here with the sign of the x-axis flipped and shifted in order to match the modulation with the temporal waveform. The temporal waveforms were reconstructed by using ATI spectra measured at the left (black solid line) and right (red dashed lines). Discussion In summary, we demonstrated a new streaking method in which an electron liberated in an intense laser field is steered by adding a weak perturbing laser field. The addition of the perturbing laser field alters the cutoff energy of an ATI spectrum. The vector potential of the perturbing laser pulse has been successfully retrieved from the momentum shift of the ATI spectrum. The temporal waveform of the laser fields has been completely characterized, including its CEP, without directional ambiguity. The temporal waveforms of chirped pulses were also completely characterized, confirming the validity of our streaking method. Our streaking technique can provide the laser field at the interaction region of above-threshold ionization, thus allowing in situ temporal characterization of a laser pulse. Since our approach can be combined with laser-induced electron diffraction experiments18,21,22,28–31 and it provides critical information on ionization and rescattering timings, it will become a versatile tool for analyzing the ultrafast electron dynamics of atoms and molecules in a strong laser field. The exact timing of the rescattering event and the accurate temporal profile of the laser pulse would provide important information in analyzing multi-electron effects in non-sequential double ionization16,17 and in time-resolved strong-field holography20,32–34. Methods Experimental methods A CEP-stabilized, 1-kHz, 30-fs, 800-nm Ti:sapphire laser was used for our experiments. The laser pulse was compressed to 5 fs using a stretched-hollow-core fiber and a set of chirped mirrors. We placed a pair of wedges made of fused silica in the beam path and controlled the wedge insertion to control the CEP of the few-cycle laser pulse. An f-2f. interferometer setup was placed before the wedges and a feedback loop was made to stabilize long-term slow drift. The long-term standard deviation of the CEP of the compressed pulse was around 300 mrad. The beams were split into two and recombined by using a mirror with a round hole at the center at an angle of 45° with a diameter of 2.5 mm. The two beams are then focused into the stereo-ATI spectrometer by a spherical concave mirror with a focal length of 50 cm. In this way, the two pulses had the same temporal waveform at the focus. We measured ATI spectra for both left and right directions with a small opening angle (2°) along the polarization of the laser field. The target cylinder with two slits in the spectrometer was filled with Xe gas with a pressure of 2.0×10-4 torr. Derivation of cutoff momentum shift Δp (Eq. 1) The semi-classical trajectory model is used to derive Eq. (1). We consider an electron that tunnels out from an atom in an intense laser field. The electron is assumed to be driven by the laser electric field in dipole approximation Et=dAt/dt after tunnelling ionization. Here, At is the vector potential of the laser field. The Coulomb field of the parent ion is neglected. The trajectories of the backward rescattered electron satisfy the rescattering condition, i.e. ∫titrAtdt=tr-tiAti, where tr is the backward rescattering time and ti is the corresponding ionization time. The final momentum pr of the backward rescattered electron is -2Atr+Ati. This final momentum is altered when a weak perturbing laser field is added to the main laser field. Because the perturbing laser field is weak, the final momentum of the electron can be expressed as 2 prτ=-2Amtr+Δtr+Amti+Δti-2Aptr+Δtr-τ+Apti+Δti-τ, where Am and Ap are the vector potentials of the main laser field and the perturbing laser field, respectively, and τ is the time delay between the two laser pulses. As the perturbing field is added to the main field, the rescattering and ionization times change slightly by Δtrτ and Δtiτ, respectively. Each term on the right-hand side of Eq. (2) can be expanded into a Taylor series at tr (or ti) because Δtr (or Δti) is small for a weak perturbing field. For example, we can approximate Amtr+Δtr as Amtr-EmtrΔtr and keep linear terms only. Thus, Eq. (2) can be rewritten as 3 prτ=-2Amtr+Amti-2Aptr-τ+1tr-ti∫titrApti-τdt We did not drop electric field terms, for example, EmtrΔtr, but they are cancelled by using a condition for pr to be extremum. Note that tr and ti are the rescattering and ionization times obtained in the absence of the perturbing field, respectively. Therefore, the shift Δp in the cutoff momentum owing to the perturbing laser field can be written as Eq. (1). Reconstruction algorithm We reconstructed the waveform of the perturbing field and the backward rescattering time using Eq. (1) in the main article. The waveforms of the main field and the perturbing field are assumed to be identical, but their amplitudes are different. Equation (1) in the main article can be viewed as a convolution between the spectral response function Sω of the backward rescattering process and the vector potential of the perturbing field. The Sω is a Fourier transform of the temporal response function St, 4 Stt=-2δt-tr+Θt-ti-Θt-trtr-ti, where Θ is the Heaviside step function, and δ is the Dirac delta function. If we find ti and tr, then we can find Sω that can be used for the deconvolution to reconstruct the perturbing field. Because Δp-τ/2 is a good approximation of the vector potential, we use Δp-τ/2 as an initial estimate for Apt and find ti and tr of the backward rescattered trajectory based on the semi-classical trajectory model. Then, a trial spectral response function Sω using Eq. (4) is found and can be used for the deconvolution of Δpτ. The next estimate of the vector potential can be found from the deconvolution result. We repeat the iterating process until the deconvolution results for Ap, ti, and tr converge. Because Δp-τ/2 is a good initial guess for the vector potential of the perturbing field the results converge rapidly within several iterations. The minimum intensity of the perturbing laser field (2×108 W/cm2) that can be measured can be estimated using the energy resolution of the TOF spectrometer (73 meV at 60 eV). The maximum intensity of the perturbing laser field that can be reliably measured were tested by solving the TDSE calculations for different intensities. These calculations showed that the intensity of the perturbing laser field should be weaker than 20% of the peak intensity of the main pulse. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Acknowledgements The authors appreciate W. T. Ra’s technical support. Author contributions Y.H.K., I.A.I, and K.T.K conceived the idea. Y.H.K performed the experiment and analysed the data. S.I.H and K.K. prepared few cycle laser pulses. Y.H.K, C.H.N, and K.T.K wrote the manuscript; and all authors discussed the results and interpretations. Funding Institute for Basic Science (IBS-R012-D1). Competing interests The authors declare no competing interests. ==== Refs References 1. Kane DJ Trebino R Characterization of arbitrary femtosecond pulses using frequency-resolved optical gating IEEE J. Quantum Electron. 1993 29 571 579 10.1109/3.199311 2. Iaconis C Walmsley IA Spectral phase interferometry for direct electric-field reconstruction of ultrashort optical pulses Opt. Lett. 1998 23 792 794 10.1364/OL.23.000792 18087344 3. Miranda M Characterization of broadband few-cycle laser pulses with the d-scan technique Opt. Express 2012 20 18732 10.1364/OE.20.018732 23038514 4. Itatani, J. et al. Attosecond streak camera. Phys. Rev. Lett.88, (2002). 5. Kienberger R Atomic transient recorder Nature 2004 427 817 821 10.1038/nature02277 14985755 6. Park SB Direct sampling of a light wave in air Optica 2018 5 402 408 10.1364/OPTICA.5.000402 7. Kim KT Photonic streaking of attosecond pulse trains Nat. Photon. 2013 7 651 656 10.1038/nphoton.2013.170 8. Kim, K. T. et al. Petahertz optical oscilloscope. Nat. Photon.7 958–962 (2013) 9. Hammond TJ Brown GG Kim KT Villeneuve DM Corkum PB Attosecond pulses measured from the attosecond lighthouse Nat. Photon. 2016 10 171 175 10.1038/nphoton.2015.271 10. Wyatt AS Attosecond sampling of arbitrary optical waveforms Optica 2016 3 303 10.1364/OPTICA.3.000303 11. Carpeggiani P Vectorial optical field reconstruction by attosecond spatial interferometry Nat. Photon. 2017 11 383 389 10.1038/nphoton.2017.73 12. Corkum PB Plasma perspective on strong field multiphoton ionization Phys. Rev. Lett. 1993 71 1994 1997 10.1103/PhysRevLett.71.1994 10054556 13. Agostini P Fabre F Mainfray G Petite G Rahman NK Free-free transitions following six-photon ionization of xenon atoms Phys. Rev. Lett. 1979 42 1127 1130 10.1103/PhysRevLett.42.1127 14. Freeman RR Above-threshold ionization with subpicosecond laser pulses Phys. Rev. Lett. 1987 59 1092 1095 10.1103/PhysRevLett.59.1092 10035138 15. L’Huillier A Lompre LA Mainfray G Manus C Multiply charged ions formed by multiphoton absorption processes in the continuum Phys. Rev. Lett. 1982 48 1814 1817 10.1103/PhysRevLett.48.1814 16. Bergues B Attosecond tracing of correlated electron-emission in non-sequential double ionization Nat. Commun. 2012 3 813 10.1038/ncomms1807 22569361 17. Liu X Nonsequential double ionization at the single-optical-cycle limit Phys. Rev. Lett. 2004 93 263001 10.1103/PhysRevLett.93.263001 15697974 18. Zuo T Bandrauk AD Corkum PB Laser-induced electron diffraction: a new tool for probing ultrafast molecular dynamics Chem. Phys. Lett. 1996 259 313 320 10.1016/0009-2614(96)00786-5 19. Meckel M Laser-induced electron tunneling and diffraction Science 2008 320 1478 1482 10.1126/science.1157980 18556555 20. Huismans Y Time-resolved holography with photoelectrons Science 2011 331 61 64 10.1126/science.1198450 21163963 21. Pullen MG Imaging an aligned polyatomic molecule with laser-induced electron diffraction Nat. Commun. 2015 6 7262 10.1038/ncomms8262 26105804 22. Wolter B Ultrafast electron diffraction imaging of bond breaking in di-ionized acetylene Science 2016 354 308 312 10.1126/science.aah3429 27846561 23. Milošević, D. B., Paulus, G. G., Bauer, D. & Becker, W. Above-threshold ionization by few-cycle pulses. J. Phys. B Atom. Mol. Opt. Phys.39, R203–R262 (2006). 24. Dudovich N Measuring and controlling the birth of attosecond XUV pulses Nat. Phys. 2006 2 781 786 10.1038/nphys434 25. Kim KT Manipulation of quantum paths for space-time characterization of attosecond pulses Nat. Phys. 2013 9 159 163 10.1038/nphys2525 26. Frolov MV Analytic model for the description of above-threshold ionization by an intense short laser pulse Phys. Rev. A 2014 89 063419 10.1103/PhysRevA.89.063419 27. Wittmann T Single-shot carrier–envelope phase measurement of few-cycle laser pulses Nat. Phys. 2009 5 357 10.1038/nphys1250 28. Meckel M Laser-induced electron tunneling and diffraction Science 2008 320 1478 1482 10.1126/science.1157980 18556555 29. Peters, M. et al. Ultrafast molecular imaging by laser-induced electron diffraction. Phys. Rev. A83, 051403 (2011). 30. Spanner, M., Smirnova, O., Corkum, P. B. & Ivanov, M. Y. Reading diffraction images in strong field ionization of diatomic molecules. J. Phys. B Atom. Mol. Opt. Phys.37, L243–L250 (2004). 31. Yurchenko, S. N., Patchkovskii, S., Litvinyuk, I. V., Corkum, P. B. & Yudin, G. L. Laser-induced interference, focusing, and diffraction of rescattering molecular photoelectrons. Phys. Rev. Lett.93, 223003 (2004). 32. Bian, X.-B. & Bandrauk, A. D. Attosecond time-resolved imaging of molecular structure by photoelectron holography. Phys. Rev. Lett.108, 263003 (2012). 33. Porat, G. et al. Attosecond time-resolved photoelectron holography. Nat. Commun.9, 1–16 (2018). 34. Bian, X.-B. et al. Subcycle interference dynamics of time-resolved photoelectron holography with midinfrared laser pulses. Phys. Rev. A84, 043420 (2011).