Research Article | Open Access

Maria Hoflund, Jasper Peschel, Marius Plach, Hugo Dacasa, Kévin Veyrinas, Eric Constant, Peter Smorenburg, Hampus Wikmark, Sylvain Maclot, Chen Guo, Cord Arnold, Anne L'Huillier, Per Eng-Johnsson, "Focusing Properties of High-Order Harmonics", *Ultrafast Science*, vol. 2021, Article ID 9797453, 8 pages, 2021. https://doi.org/10.34133/2021/9797453

# Focusing Properties of High-Order Harmonics

#### Abstract

Many applications of the extreme ultraviolet (XUV) radiation obtained by high-order harmonic generation (HHG) in gases require a small focus area in order to enable attosecond pulses to reach a high intensity. Here, high-order harmonics generated in Ar with a multiterawatt laser system in a loose focusing geometry are focused to a few micrometers using two toroidal mirrors in a Wolter configuration with a high demagnification factor. Using a knife-edge measurement technique, we determine the position and size of the XUV foci as a function of harmonic order. We show that the focus properties vary with harmonic order and the generation conditions. Simulations, based on a classical description of the harmonic dipole phase and assuming that the individual harmonics can be described as Gaussian beams, reproduce the experimental behavior. We discuss how the generation geometry affects the intensity and duration of the focused attosecond pulses.

#### 1. Introduction

High-order harmonic generation (HHG) is a highly nonlinear process that converts intense low-frequency light into attosecond pulse trains in the extreme ultraviolet (XUV) spectral region. In spite of a low conversion efficiency [1–3], the choice of appropriate gas targets combined with the use of loose focusing geometries has led to high pulse energies on the order of microjoules [4–8]. The high intensity that can potentially be achieved by focusing this radiation opens the way to multiphoton ionization processes in the XUV range [7, 9–12], and in particular to nonlinear XUV pump-XUV probe studies with attosecond resolution [13–15]. While XUV and X-ray radiation from free-electron lasers provide the required intensities for these types of studies [16, 17], they have only recently started to demonstrate the generation of attosecond pulses [18]. Using these sources for experiments is still very challenging, whereas HHG-based sources already provide a mature technology for nonlinear pump-probe experiments in the attosecond range.

To obtain both high intensity and attosecond pulse duration in the interaction with a low-density gas target, it is important to focus the broadband radiation on a small spot. This requires focusing the different harmonic components at the same position. Focusing of high-order harmonics has been achieved using ellipsoidal mirrors [19, 20], pairs of toroidal mirrors at grazing incidence [21–23], or spherical mirrors at normal incidence for low photon energy [8, 24]. Focal spot sizes of the order of a few *μ*m for the emitted harmonics in a certain spectral range are typically obtained. The methods used for measuring these spot sizes vary from direct techniques like the knife-edge method [20, 25–27] or microscopy [21] to indirect ones based on wavefront determination [23, 28–30]. Recently, focus sizes of the order of *μ*m have been measured for harmonics between 10 and 20 nm (i.e., 60-120 eV) focused with an ellipsoidal mirror [31].

Previous studies [30, 32–36] have shown that high-order harmonics have different wavefronts at the exit of the nonlinear medium, which can be convergent or divergent [35], corresponding to a real or virtual focus, respectively. The wavefront curvature of the harmonics is the sum of that of the driving laser and an additional contribution originating from the process behind the generation of high-order harmonics. The resulting harmonic radius of curvature and spatial profile [37] depends on harmonic order, laser intensity, and generation position. When they are refocused, this results in separated waist positions along the propagation direction, thus affecting the focusing properties of the harmonics and the peak intensity and pulse duration of the focused attosecond pulses. Although theoretically predicted in [34, 35], and experimentally shown for the total harmonic beam in [38], no experiments have yet quantitatively examined the influence of the generation geometry on the focus position, waist, and intensity of individual harmonics and attosecond pulses.

In the present study, we generate high-order harmonics in argon using a 45 mJ, 40 fs, Ti:Sapphire laser in a loose focusing geometry. We refocus the XUV radiation with a pair of toroidal mirrors in a Wolter configuration [23, 39] in order to achieve broadband focusing with a high demagnification factor. A knife-edge technique [20, 25–27] is used to measure the individual harmonic focus positions for different locations of the gas target relative to the laser focus in the generation chamber. In some cases, the harmonics are found to be refocused at positions that evolve strongly with the harmonic order. The observed experimental trends can be explained by a simple model which takes into account the intensity dependence of the spectral phase of the attosecond pulses, also called the dipole phase [34, 35]. We find that the chromatic aberrations induced by the dipole phase affect strongly the temporal profile and intensity of the focused attosecond pulses.

#### 2. Experimental Method

The experimental setup used in this study is sketched in Figure 1. We use a multiterawatt Ti:Sapphire laser system delivering 40 fs, 806 nm, and 45 mJ pulses [7]. A deformable mirror (DM), used together with an infrared wavefront sensor, allows us to correct for spatial aberrations and to control the IR focusing position with respect to the generation medium by introducing a slight curvature to the beam. After aberration correction, the measured laser wavefront has an RMS value of about /40, corresponding to a Strehl ratio of 0.97. After going through an iris, the infrared beam is loosely focused with an 8.7 m spherical mirror into a 1 cm long, pulsed, Ar gas cell, generating high-order harmonics up to 45 eV. We estimate that the IR waists are approximately 340-360 *μ*m with corresponding Rayleigh lengths of 45 to 50 cm for the different focal lengths used in the experiment.

The XUV radiation is reflected, at a 10° grazing incidence angle, by a fused silica (FS) plate, which is antireflection-coated for the IR, and spectrally filtered by a 200 nm Al filter. The FS plate and Al filter allow us to completely eliminate the IR. Further, the broadband beam is refocused by Wolter optics, consisting of a prealigned assembly of two toroidal mirrors from Thales SESO, designed to image the generation medium, after 6 m propagation, with a high demagnification factor and minimal aberrations in the XUV range. The grazing incidence angle on the two toroidal mirrors is 15°, giving a total deviation angle of 60°. The Wolter optics assembly is aligned in our beamline by using a homemade XUV Hartmann wavefront sensor. The tangential focal lengths are 265.3 mm and 545.2 mm for the first and second mirror while the sagittal focal lengths are 265.0 mm and 544.4 mm. The combined focal length of 164.2 mm provides a demagnification factor of 35. A homemade flat-field spectrometer, based on a Hitachi aberration-corrected concave grating, images the focus in the vertical direction onto a microchannel plate (MCP). In this configuration, no entrance slit is needed and the entire beam is characterized. The grating disperses the harmonics along the vertical axis of the MCP, while the far field spatial profile is recorded along the horizontal axis.

For each position of the cell relative to the IR focus along the propagation direction , called “generation position,” we measure the focus position of the harmonics in the application chamber. Harmonic spectra are recorded for different (transverse to the beam) and (along the beam) positions of a knife-edge located close to the focus (see Figure 1). A rough estimate of the focus position for each harmonic can be obtained by looking at the diffraction image position of the knife-edge on the spatial profile recorded by the spectrometer. If the knife is after the focus, the shadow is on the same side as the knife, while if the focus is after, it is on the opposite side. This is illustrated in Figure 2, showing spectra obtained without the knife-edge (a) and for three positions of the knife-edge (b–d). In (b) (respectively, d), the knife is before (after) the focus for all harmonic orders and the shadow is observed on the right (left) sides of the profile. In (c), the shadow of the knife is on the right side for the lower-order harmonics while it is on the left side for the higher-order harmonics. This means that the knife is before the focus for the lower-order harmonics, and after the focus for the higher orders. The harmonics in the middle are less affected by the knife, which indicates that it is located close to their focus. These results show that the harmonics are indeed focused at different positions along the propagation axis.

Mapping out the individual harmonic focus positions requires scanning both and in small steps. Instead, we use a faster and less tedious method based on measurements at only four positions in , outside the Rayleigh range ( mm), two before and two after the focus, which is sufficient to extract the position of each focus. Figure 3 illustrates this method for the 21^{st} harmonic. The harmonic intensity as a function of knife insertion () is recorded for each knife-edge position along the propagation direction () (see Figure 3(a)). Assuming Gaussian beams, this data is then fitted to a normalized error function, , representing the integral of a Gaussian function from to . For each knife-edge position, the beam waist radius () is obtained from the fit. The position of the harmonic focus () can be determined by performing a linear fit of the beam widths, after mirroring those corresponding to the knife-edge position behind the focus (Figure 3(b)). If we assume that the harmonics are ideal Gaussian beams with , the beam waists () can also be estimated as , where is the wavelength of the harmonic and the (half) divergence angle. Note that represents a lower limit to the beam waist and that a better estimation would require a measurement of the quality factor , which has not been done in the present work.

#### 3. Theoretical Model

HHG can be described by a three-step process: tunneling ionization, electron acceleration by the laser field, and recombination [40, 41]. In the second step, the electron wave packet accumulates a phase, called the dipole phase, which is then transferred to the harmonic field. Two families of trajectories, related to the time spent by the electron in the continuum (“short” or “long”), have been identified and observed experimentally [42, 43]. The dipole phase is often expressed as , where depends on the harmonic order and electron trajectory, and is the intensity of the driving laser [44]. This expression is derived by using the strong-field approximation to solve the time-dependent Schrödinger equation (TDSE) of the single-atom response. Instead, we use an analytical expression derived by solving Newton’s equation of motion [34, 45]: where refers to the trajectory (short or long), is the harmonic order, and are the fundamental intensity and angular frequency, respectively, is the angular frequency corresponding to the ionization energy, and is a return time, as explained in [34]. We here consider only the contribution from the short trajectory, which is mainly observed in our experiment. In Equation (1), the quantity is independent of the harmonic order and is equal to zero for the short trajectory, while the quantity is equal to s W cm for a driving laser of 800 nm. The first and third terms depend on the laser intensity, leading to an order-dependent wavefront curvature.

We assume that high-order harmonics are generated at position relative to the IR focus (located at ) in a thin medium [46]. For a Gaussian fundamental beam, with a profile , where is the peak intensity, is the radial coordinate, and is the beam width at position , the -dependence of the dipole phase can be approximated by

The generated harmonic field also includes an -dependent phase contribution from the fundamental, which is equal to , where is the fundamental wavevector and is its radius of curvature. The total radius of curvature of the harmonic field, corresponding to one trajectory , is given by where is the speed of light and is the fundamental beam waist [34].

To fully characterize the harmonic beams, the width at the generation position is estimated by assuming that the harmonic intensity () is proportional to the fundamental intensity to the power of () [46–50], where the degree of nonlinearity, , is assumed to be equal to 4, which is confirmed by TDSE calculations [34]. The harmonic width, at the position of the gas target, is then simply half the width of the driving laser. Assuming Gaussian optics for the harmonic beam, knowing the width and radius of curvature at the generation position allows us to determine the beam properties at any position.

Figure 4 shows the calculated focus positions and waist size, with respect to the IR focus, for harmonic 11 to 25 as a function of generation position . The left column shows the (virtual or real) focus positions and waists at generation while the right column shows the foci and waists of the refocused beams. The first row (a, b) shows results obtained without the dipole phase. In this case, the harmonics have the same wavefront, which is equal to that of the fundamental. When they are generated at the IR waist, they have the same focus, which is identical to the IR focus (Figure 4(a)). However, for this generation position (), the harmonics are not refocused at the same focus position in the application chamber (Figure 4(b)). In this very loose focusing geometry, some harmonics have very long Rayleigh lengths. Consequently, the focusing optics is not fully in the far field, so that the curvatures of the harmonics at this position are slightly harmonic-dependent. Therefore, the harmonics will be focused at slightly different positions. In this case, where the dipole phase is not included, the foci in the application chamber are overlapped for all harmonics, when the generation medium is slightly before the IR focal spot. This is a clear deviation from the predictions of geometrical optics, due to the long Rayleigh lengths mentioned above.

The second row (c, d) is calculated including the dipole phase. This leads to an asymmetric focus position distribution for positive and negative generation positions. Furthermore, there is no generation position where all harmonics are focused at the same position. We also show in (e, f) the harmonic waists both in the generation and application chambers, including the dipole phase in the calculations.

Finally, we also simulate the experimental conditions, where the cell position is fixed and where the focal length is slightly varied in order to change the relative position between the generation and the IR focus. As expected, the results, shown as dashed lines in Figures 4(d) and 4(f), present the same qualitative behavior as when the focal length is not varied (solid lines), and only minor quantitative differences. The simulation with the varying focal length is used in the comparison with the experimental results.

#### 4. Results

To describe our experimental geometry, where the generation cell is fixed and the position of the IR focus varied, we introduce , which is the medium position () relative to the IR focus (). Similarly, in the application chamber, we refer to the refocused focus positions in the coordinate , where is the position of the refocused IR, which was estimated from the simulations. The 19^{th} harmonic is removed from the dataset because its position on the detector coincides with a damaged area on the microchannel plate of the spectrometer, making its detection unreliable.

The results of our knife-edge measurements are shown in Figures 5(a)–5(e), for different generation positions, as shown at the top of the figure. The experimental data are compared with simulations (dashed line) using the model presented above. The different harmonics are found to be focused at different positions, with a separation of up to 1.5 mm, which is significant since the harmonic Rayleigh lengths are typically around 1 mm. This effect decreases as the generation position moves towards and beyond the laser focus. When the medium is before the focus (see Figures 5(a)–5(c)), the lower-order harmonics are focused approximately one mm after the higher-order harmonics, on different sides of the IR focus. When the generation medium is centered at the infrared focus or beyond (see Figures 5(d) and 5(e)), all of the harmonics are focused much closer to each other, at a position before the IR focus. Simulations reproduce the main features of the experimental data, which indicates that, despite its simplicity, our model describes the physics behind the focusing properties of high-order harmonics well [34, 35].

In Figures 6(a)–6(e), we present the measured (solid line) and simulated (dashed line) beam waists (), for a few harmonics (13, 15, 17, 21, and 23), as a function of generation position. As explained in the experimental section, the measured beam waists should be considered lower limits. We find that the harmonic waist sizes decrease with order. For harmonics 13, 15, and 17 (a–c), the generation position leading to the maximum waist moves progressively from around to negative values. The measured values for the harmonic beam waists are consistent with (but slightly below) previous spectrally integrated measurements, which were performed by wavefront analysis of the total beam combined with numerical back-propagation and gave *μ*m FWHM ( *μ*m half-width at ) [23]. The simulations in Figure 6 (a–e, dashed lines) follow qualitatively the experimental results, with slightly larger waists. This difference might be due to a deviation of the experimental laser (and consequently harmonic) beams from a Gaussian behavior (i.e., ). Another possible cause of the discrepancy between the experimental results and our simulations could be that the dipole phase is more complex than the model dipole or that phase-matching impacts also the wavefront of the XUV beam. Despite the approximation of our model, the agreement is very satisfactory.

Figures 6(f)–6(j) show the focus positions (blue) and yield (orange) for the same harmonics, as a function of generation position. Simulations, indicated by dashed lines, agree well with the experimental measurements. As shown in Figure 5, the harmonic focus position varies with the generation geometry, especially for the lower-order harmonics. The harmonic yield, on the other hand, is relatively insensitive to the generation position over the studied range.

#### 5. Discussion

Our results can be used to estimate not only the intensity of the refocused individual harmonics but also that of the attosecond pulses in the train, in which case the relative phase of the harmonics has to be included [8]. The electric field describing the attosecond pulse train can be expressed as where are the relative experimental harmonic intensities. The phases, , are the sum of the Gouy phase obtained by using the experimentally determined position of the focus and Rayleigh length for each harmonic and the inherent phase due to the generation process (see Equation (1)), also called attochirp [51].

In Figures 7(a)–7(c), we show the estimated attosecond pulse intensity on axis in the application chamber, as a function of time and position along the propagation axis (), for three generation positions, mm (a), (b), and mm (c). This is obtained using Equation (4) for each and taking the absolute square to get the intensity. The peak of the attosecond pulse intensity is extracted from these results and shown in Figure 7(d) in black as a function of , for different generation positions (). The intensity variation for each harmonic is marked in color going from the 11^{th} (red) to 25^{th} (blue). Figure 7(d) also shows the maximum XUV peak intensity (green) and the total yield (red) for each generation position. Furthermore, the intensity distribution as a function of the transverse coordinate and is shown for mm (e), (f), and mm (g).

As seen in Figure 7, the distribution of harmonic focus positions is clearly broader when the cell is before the IR focus (negative ), as shown in panels (a) and (e). In contrast, the foci are closer to each other along the propagation axis when the cell is located at the IR focus or beyond, as shown in panels (b, c) and (f, g). The maximum intensity is obtained for and mm, as seen in panels (b) and (f). When , the harmonics are focused at different positions, while when , the generation efficiency decreases. The total yield and the maximum intensity vary differently with respect to the generation position, as shown in Figure 7(d).

From the estimated spatiotemporal intensities in Figure 7, we find for the conditions where the maximum intensity is obtained ( and mm), an average pulse duration of the individual attosecond pulses of 200 as and a beam waist of 1.5 *μ*m. Combining this with the XUV pulse energy, measured using an X-ray CCD camera (Andor) to be 1.6 nJ, and a harmonic pulse duration of 20 fs (half of the IR pulse duration, according to the same power law used for the harmonic beam widths above), gives an estimate of the maximum XUV peak intensity of W/cm. This corresponds to an intensity of 1 on the vertical axis in Figure 7(d).

#### 6. Conclusion

By inserting a knife before and after the focus and scanning it through the XUV beam, we have mapped out the focus position and estimated the waist size of each individual harmonic, as a function of the position of the gas cell relative to the IR focus. The key conclusion of this work is that the harmonics’ focusing properties are found to depend on the order and the generation position. We interpret this effect as a consequence of the intensity-dependent dipole phase. The experimental data agree with relatively simple simulations based on an analytical expression for the dipole phase [34, 35], the assumption that the harmonics can be represented as Gaussian beams [33], and without taking into account propagation effects in the nonlinear medium. These measurements allow us to estimate the focused intensity of the individual harmonics and the attosecond pulses in the train and to define the optimal conditions to achieve high-intensity attosecond illumination. We find that these optimal conditions do not necessarily coincide with those leading to the highest yield. Even when the yield is approximately constant, the harmonic focusing properties can change significantly and thus affect the pulse duration and intensity of the attosecond pulses.

#### Data Availability

All data discussed in the article will be made available upon request.

#### Disclosure

Hugo Dacasa’s present address is Department of Physical Chemistry, Complutense University of Madrid, Ciudad Universitaria, 28040 Madrid, Spain. Kévin Veyrinas’ present address is Amplitude Laser Group, 11, avenue de Canteranne, Cité de la Photonique, 33600 Pessac, France. Hampus Wikmark’s present address is Department of Physics and Astronomy, Uppsala University, Box 516, SE-751 20, Uppsala, Sweden. Sylvain Maclot’s present address is Department of Physics, Gothenburg University, Box 100, 405 30 Gothenburg, Sweden.

#### Conflicts of Interest

The authors declare no conflict of interest.

#### Authors’ Contributions

M.H, J.P., M.P., and H.D. performed the experiments. P.E.-J, J.P., and C.A. designed the experiment. K.V., E.C., H.W., and S.M. contributed to building the experimental setup and preparing the experiment. M.H. did the data analysis and the simulations based on theory that P.S., H.W., C.G., C.A., and A.L. developed. M.H., A.L., P.E.-J., J.P., and M.P. wrote the article, with feedback from all the authors.

#### Acknowledgments

The authors acknowledge support from the Swedish Research Council, the European Research Council (advanced grant QPAP), the Knut and Alice Wallenberg Foundation, and the Crafoord Foundation. The research leading to these results has received funding from LASERLAB-EUROPE (grant agreement no. 654148, European Union’s Horizon 2020 research and innovation programme). S.M. acknowledges financial support from the COST Action CA18212-Molecular Dynamics in the GAS phase (MD-GAS), supported by COST (European Cooperation in Science and Technology).

#### References

- P. Agostini and L. F. DiMauro, “The physics of attosecond light pulses,”
*Reports on Progress in Physics*, vol. 67, no. 6, pp. 813–855, 2004. View at: Publisher Site | Google Scholar - J. Tate, T. Auguste, H. G. Muller, P. Salières, P. Agostini, and L. F. DiMauro, “Scaling of wave-packet dynamics in an intense midinfrared field,”
*Physical Review Letters*, vol. 98, no. 1, article 013901, 2007. View at: Publisher Site | Google Scholar - A. D. Shiner, C. Trallero-Herrero, N. Kajumba et al., “Wavelength scaling of high harmonic generation efficiency,”
*Physical Review Letters*, vol. 103, no. 7, article 073902, 2009. View at: Publisher Site | Google Scholar - Y. Kobayashi, T. Ohno, T. Sekikawa, Y. Nabekawa, and S. Watanabe, “Pulse width measurement of high-order harmonics by autocorrelation,”
*Applied Physics B: Lasers and Optics*, vol. 70, no. 3, pp. 389–394, 2000. View at: Publisher Site | Google Scholar - J.-F. Hergott, M. Kovacev, H. Merdji et al., “Extreme-ultraviolet high-order harmonic pulses in the microjoule range,”
*Physical Review A*, vol. 66, no. 2, article 021801, 2002. View at: Publisher Site | Google Scholar - P. Rudawski, C. M. Heyl, F. Brizuela et al., “A high-flux high-order harmonic source,”
*Review of Scientific Instruments*, vol. 84, no. 7, article 073103, 2013. View at: Publisher Site | Google Scholar - B. Manschwetus, L. Rading, F. Campi et al., “Two-photon double ionization of neon using an intense attosecond pulse train,”
*Physical Review A*, vol. 93, no. 6, article 061402, 2016. View at: Publisher Site | Google Scholar - I. Makos, I. Orfanos, A. Nayak et al., “
*Α*10-gigawatt attosecond source for non-linear XUV optics and XUV-pump-XUV- probe studies,”*Scientific Reports*, vol. 10, no. 1, p. 3759, 2020. View at: Publisher Site | Google Scholar - K. Midorikawa, Y. Nabekawa, and A. Suda, “XUV multiphoton processes with intense high-order harmonics,”
*Progress in Quantum Electronics*, vol. 32, no. 2, pp. 43–88, 2008. View at: Publisher Site | Google Scholar - E. J. Takahashi, P. Lan, O. D. Mücke, Y. Nabekawa, and K. Midorikawa, “Attosecond nonlinear optics using gigawatt-scale isolated attosecond pulses,”
*Nature Communications*, vol. 4, no. 1, 2013. View at: Publisher Site | Google Scholar - A. Nayak, I. Orfanos, I. Makos et al., “Multiple ionization of argon via multi-xuv-photon absorption induced by 20-gw high-order harmonic laser pulses,”
*Physical Review A*, vol. 98, no. 2, article 023426, 2018. View at: Publisher Site | Google Scholar - B. Senfftleben, M. Kretschmar, A. Hoffmann et al., “Highly non-linear ionization of atoms induced by intense high-harmonic pulses,”
*Journal of Physics: Photonics*, vol. 2, no. 3, 2020. View at: Publisher Site | Google Scholar - P. A. Carpeggiani, P. Tzallas, A. Palacios, D. Gray, F. Martín, and D. Charalambidis, “Disclosing intrinsic molecular dynamics on the 1-fs scale through extreme-ultraviolet pump-probe measurements,”
*Physical Review A*, vol. 89, no. 2, article 023420, 2014. View at: Publisher Site | Google Scholar - P. Tzallas, D. Charalambidis, N. A. Papadogiannis, K. Witte, and G. D. Tsakiris, “Direct observation of attosecond light bunching,”
*Nature*, vol. 426, no. 6964, pp. 267–271, 2003. View at: Publisher Site | Google Scholar - Y. Nabekawa, H. Hasegawa, E. J. Takahashi, and K. Midorikawa, “Production of doubly charged helium ions by two-photon absorption of an intense sub-10-fs soft X-ray pulse at 42 eV photon energy,”
*Physical Review Letters*, vol. 94, no. 4, article 043001, 2005. View at: Publisher Site | Google Scholar - C. Pellegrini, A. Marinelli, and S. Reiche, “The physics of x-ray free-electron lasers,”
*Reviews of Modern Physics*, vol. 88, no. 1, article 015006, 2016. View at: Publisher Site | Google Scholar - L. Young, K. Ueda, M. Gühr et al., “Roadmap of ultrafast x-ray atomic and molecular physics,”
*Journal of Physics B: Atomic, Molecular and Optical Physics*, vol. 51, no. 3, p. 032003, 2018. View at: Publisher Site | Google Scholar - J. Duris, S. Li, T. Driver et al., “Tunable isolated attosecond x-ray pulses with gigawatt peak power from a free-electron laser,”
*Nature Photonics*, vol. 14, no. 1, pp. 30–36, 2020. View at: Publisher Site | Google Scholar - A. Suda, H. Mashiko, and K. Midorikawa, “Focusing intense high-order harmonics to a micron spot size,” in
*Progress in Ultrafast Intense Laser Science II*, pp. 183–198, Springer Berlin Heidelberg, Berlin, Heidelberg, 2007. View at: Google Scholar - H. Motoyama, T. Sato, A. Iwasaki et al., “Development of high-order harmonic focusing system based on ellipsoidal mirror,”
*Review of Scientific Instruments*, vol. 87, no. 5, article 051803, 2016. View at: Publisher Site | Google Scholar - C. Valentin, D. Douillet, S. Kazamias et al., “Imaging and quality assessment of high-harmonic focal spots,”
*Optics Letters*, vol. 28, no. 12, pp. 1049–1051, 2003. View at: Publisher Site | Google Scholar - F. Frassetto, S. Anumula, F. Calegari, A. Trabattoni, M. Nisoli, and L. Poletto, “Microfocusing beamline for xuv-xuv pump-probe experiments using hh generation,” in
*Adaptive X-Ray Optics III*, vol. 9208, pp. 159–167, September 2014. View at: Publisher Site | Google Scholar - H. Coudert-Alteirac, H. Dacasa, F. Campi et al., “Micro-focusing of broadband highorder harmonic radiation by a double toroidal mirror,”
*Applied Sciences*, vol. 7, no. 11, p. 1159, 2017. View at: Publisher Site | Google Scholar - I. Orfanos, I. Makos, I. Liontos et al., “Non-linear processes in the extreme ultraviolet,”
*Journal of Physics: Photonics*, vol. 2, no. 4, 2020. View at: Publisher Site | Google Scholar - J. A. Arnaud, W. M. Hubbard, G. D. Mandeville, B. de la Clavière, E. A. Franke, and J. M. Franke, “Technique for fast measurement of Gaussian laser beam parameters,”
*Applied Optics*, vol. 10, no. 12, pp. 2775-2776, 1971. View at: Publisher Site | Google Scholar - D. R. Skinner and R. E. Whitcher, “Measurement of the radius of a high-power laser beam near the focus of a lens,”
*Journal of Physics E: Scientific Instruments*, vol. 5, no. 3, pp. 237-238, 1972. View at: Publisher Site | Google Scholar - L. le Déroff, P. Salières, and B. Carré, “Beam-quality measurement of a focused high-order harmonic beam,”
*Optics Letters*, vol. 23, no. 19, pp. 1544–1546, 1998. View at: Publisher Site | Google Scholar - D. G. Lee, J. J. Park, J. H. Sung, and C. H. Nam, “Wave-front phase measurements of high-order harmonic beams by use of point-diffraction interferometry,”
*Optics Letters*, vol. 28, no. 6, pp. 480–482, 2003. View at: Publisher Site | Google Scholar - E. Frumker, G. G. Paulus, H. Niikura, D. M. Villeneuve, and P. B. Corkum, “Frequency-resolved high-harmonic wavefront characterization,”
*Optics Letters*, vol. 34, no. 19, pp. 3026–3028, 2009. View at: Publisher Site | Google Scholar - D. R. Austin, T. Witting, C. A. Arrell et al., “Lateral shearing interferometry of high-harmonic wavefronts,”
*Optics Letters*, vol. 36, no. 10, pp. 1746–1748, 2011. View at: Publisher Site | Google Scholar - H. Motoyama, A. Iwasaki, Y. Takei et al., “Broadband nano-focusing of high-order harmonics in soft x-ray region with ellipsoidal mirror,”
*Applied Physics Letters*, vol. 114, no. 24, article 241102, 2019. View at: Publisher Site | Google Scholar - D. T. Lloyd, K. O’Keeffe, and S. M. Hooker, “Complete spatial characterization of an optical wavefront using a variable-separation pinhole pair,”
*Optics Letters*, vol. 38, no. 7, pp. 1173–1175, 2013. View at: Publisher Site | Google Scholar - E. Frumker, G. G. Paulus, H. Niikura, A. Naumov, D. M. Villeneuve, and P. B. Corkum, “Order-dependent structure of high harmonic wavefronts,”
*Optics Express*, vol. 20, no. 13, article 13870, p. 13877, 2012. View at: Publisher Site | Google Scholar - H. Wikmark, C. Guo, J. Vogelsang et al., “Spatiotemporal coupling of attosecond pulses,”
*Proceedings of the National Academy of Sciences*, vol. 116, no. 11, pp. 4779–4787, 2019. View at: Publisher Site | Google Scholar - L. Quintard, V. Strelkov, J. Vabek et al., “Optics-less focusing of XUV high-order harmonics,”
*Science Advances*, vol. 5, no. 4, article eaau7175, 2019. View at: Publisher Site | Google Scholar - H. Dacasa, H. Coudert-Alteirac, C. Guo et al., “Single-shot extreme-ultraviolet wavefront measurements of high-order harmonics,”
*Optics Express*, vol. 27, no. 3, pp. 2656–2670, 2019. View at: Publisher Site | Google Scholar - C. Hernández-García, A. Picón, J. San Román, and L. Plaja, “Attosecond extreme ultraviolet vortices from high-order harmonic generation,”
*Physical Review Letters*, vol. 111, no. 8, article 083602, 2013. View at: Publisher Site | Google Scholar - S. Chatziathanasiou, S. Kahaly, D. Charalambidis, P. Tzallas, and E. Skantzakis, “Imaging the source of high-harmonics generated in atomic gas media,”
*Optics Express*, vol. 27, no. 7, pp. 9733–9739, 2019. View at: Publisher Site | Google Scholar - H. Wolter, “Spiegelsysteme streifenden Einfalls als abbildende Optiken für Röntgenstrahlen,”
*Annalen der Physik*, vol. 445, no. 1-2, pp. 94–114, 1952. View at: Publisher Site | Google Scholar - P. B. Corkum, “Plasma perspective on strong field multiphoton ionization,”
*Physical Review Letters*, vol. 71, no. 13, pp. 1994–1997, 1993. View at: Publisher Site | Google Scholar - K. J. Schafer, B. Yang, L. F. DiMauro, and K. C. Kulander, “Above threshold ionization beyond the high harmonic cutoff,”
*Physical Review Letters*, vol. 70, no. 11, pp. 1599–1602, 1993. View at: Publisher Site | Google Scholar - M. Bellini, C. Lyngå, A. Tozzi et al., “Temporal coherence of ultrashort high-order harmonic pulses,”
*Physical Review Letters*, vol. 81, no. 2, pp. 297–300, 1998. View at: Publisher Site | Google Scholar - P. Salières, B. Carré, L. le Déroff et al., “Feynman’s path-integral approach for intense-laser-atom interactions,”
*Science*, vol. 292, no. 5518, pp. 902–905, 2001. View at: Publisher Site | Google Scholar - M. Lewenstein, K. C. Kulander, K. J. Schafer, and P. H. Bucksbaum, “Rings in above-threshold ionization: a quasiclassical analysis,”
*Physical Review A*, vol. 51, no. 2, pp. 1495–1507, 1995. View at: Publisher Site | Google Scholar - C. Guo, A. Harth, S. Carlström et al., “Phase control of attosecond pulses in a train,”
*Journal of Physics B: Atomic Molecular and Optical Physics*, vol. 51, no. 3, article 034006, 2018. View at: Publisher Site | Google Scholar - C. Hernández-García, J. S. Román, L. Plaja, and A. Picón, “Quantum-path signatures in attosecond helical beams driven by optical vortices,”
*New Journal of Physics*, vol. 17, no. 9, article 093029, 2015. View at: Publisher Site | Google Scholar - F. Catoire, A. Ferré, O. Hort et al., “Complex structure of spatially resolved high-order-harmonic spectra,”
*Physical Review A*, vol. 94, no. 6, 2016. View at: Publisher Site | Google Scholar - L. Rego, J. S. Román, A. Picón, L. Plaja, and C. Hernández-García, “Nonperturbative twist in the generation of extreme-ultraviolet vortex beams,”
*Physical Review Letters*, vol. 117, no. 16, 2016. View at: Publisher Site | Google Scholar - J. A. Pérez-Hernández, L. Roso, and L. Plaja, “Harmonic generation beyond the Strong-Field approximation: the physics behind the short-wave-infrared scaling laws,”
*Optics Express*, vol. 17, no. 12, pp. 9891–9903, 2009. View at: Publisher Site | Google Scholar - C. G. Durfee, A. R. Rundquist, S. Backus, C. Herne, M. M. Murnane, and H. C. Kapteyn, “Phase matching of high-order harmonics in hollow waveguides,”
*Physical Review Letters*, vol. 83, no. 11, pp. 2187–2190, 1999. View at: Publisher Site | Google Scholar - Y. Mairesse, A. de Bohan, L. J. Frasinski et al., “Attosecond synchronization of high-harmonic soft X-rays,”
*Science*, vol. 302, no. 5650, pp. 1540–1543, 2003. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2021 Maria Hoflund et al. Exclusive Licensee Xi’an Institute of Optics and Precision Mechanics. Distributed under a Creative Commons Attribution License (CC BY 4.0).