## Abstract

The physics of high harmonics has led to the generation of attosecond pulses and to trains of attosecond pulses. Measurements that confirm the pulse duration are all performed in the far field. All pulse duration measurements tacitly assume that both the beam’s wavefront and intensity profile are independent of frequency. However, if one or both are frequency dependent, then the retrieved pulse duration depends on the location where the measurement is made. We measure that each harmonic is very close to a Gaussian, but we also find that both the intensity profile and the beam wavefront depend significantly on the harmonic order. Thus, our findings mean that the pulse duration will depend on where the pulse is observed. Measurement of spectrally resolved wavefronts along with temporal characterization at one single point in the beam would enable complete space-time reconstruction of attosecond pulses. Future attosecond science experiments need not be restricted to spatially averaged observables. Our approach paves the way to recovery of the single molecule response to the strong field.

© 2012 OSA

## 1. Introduction

High harmonic generation is a highly coherent process in which many photons are converted into one, leaving the generating atoms in the nonlinear medium unchanged. High harmonic generation is often described with a semiclassical model in which an electron tunnel ionizes, follows a semiclassical trajectory and then recombines by emitting a single photon. In a gas of many atoms a phase matched output results in laser-like coherence of the harmonics, which has been confirmed by interfering spectrally resolved replica beams [1] or different spatial regions of the same beam [2]. Several techniques for high-harmonic wavefront measurements were demonstrated as well [3, 4]. These measurements however, do not address differences in the amplitude or curvature of wavefronts between harmonics, even though we know that the wavefront depends upon the trajectory [5–7]. Such a differences will be important for all attosecond pulse and pulse train applications.

Recently we demonstrated a technique capable of measuring spectrally and spatially resolved amplitude and phases of high harmonics (or SWORD-Spectral Wavefront Optical Reconstruction by Diffraction) [8]. SWORD is a linear measurement that utilizes direct wavefront reconstruction algorithm (i.e. no iterative computations are required for the reconstruction, eliminating potential ambiguities). Spectrally resolved wavefronts can also be measured by producing two sources of harmonics in the generating jet and using lateral shearing interferometry [9].

Our SWORD measurements reveal that the short-trajectory harmonics have an essentially Gaussian spatial profile and that the root-mean-square deviations from a parabolic wavefront are only ≲ 0.2 radians. They also reveal that different harmonics have a significantly different wavefront curvature upon emerging from the generating region. Since we measured both amplitude and phase of each harmonic, we are able to propagate the fields and determine their spatial properties anywhere along the direction of propagation. Specifically we determine their waist position (position of the flat phase) and size - the virtual origin of the coherent beams. We find that this virtual origin of each harmonic is at significantly different position. Particularly important, the distance between the waist positions of some harmonics can be larger than their respective Rayleigh ranges. Since the harmonics are generated in the gas jet, this observation means that they are produced with different radii-of-curvature, depending on the harmonic order. It also means that, if the harmonics are focused to measure their temporal structure, the duration of the pulse will depend on where the measurement is made.

## 2. Experimental details

In our experiment we measure high harmonics produced by focusing an 800nm, 35fsec, 400*μ*J, laser beam onto a Nitrogen gas jet. The laser beam from Ti:Sapphire amplifier was spatially filtered by propagating it through hollow-core fiber in vacuum and the beam was then focused with an *f*/# = 80 lens onto the pulsed gas jet. The gas jet in which the harmonics were generated has 250*μm* aperture and was operated with a backing pressure of 2.7atm.

The Rayleigh range of the fundamental beam (*Z _{R}* ≃ 10mm) was much longer than the interaction region of the jet (

*d*≃ 0.4mm) where harmonics were generated. The intensity of the fundamental beam was kept moderate to prevent excessive ionization and plasma creation. This configuration is particularly useful when we aim to study single atom/molecular response and has been widely used in attosecond science [10].

For SWORD [8] a small 20*μ*m wide slit is scanned across the XUV beam. The XUV spectrometer is built with variable groove spacing flat field Hitachi imaging grating (model 001-0266). The system resolves the spectrum in horizontal direction, and allows essentially free-space field propagation in vertical direction. The harmonic signal was recorded using a 40mm diameter Burle imaging microchannel plate (MCP) in the imaging plane of the spectrometer and the back-side phosphor screen was imaged onto the cooled CCD camera. Two-dimensional images are taken for each position of the scanning slit. Using this diffractogram, we reconstruct spectral wavefront phase and amplitudes of each harmonic.

As the fundamental beam size in the interaction region is much smaller than the gas jet size, we can assume a uniform gas density across the beam in the generation region and therefore rotational symmetry of the generated harmonic beam. Figure 1(a) shows examples of the reconstruction for harmonics 15 and 21.

The relative vertical position of the diffraction pattern’s centroid is proportional to the wavefront slope of the sampled wavefront slice. We reconstruct the wavefront phase in Fig. 1(a) using the following relation:
${\delta}_{j}=\delta z\cdot \frac{{y}_{j}-{z}_{j}}{d}$, where *y _{j}* - is the relative position of the centroid of the diffracted pattern,

*z*- relative slit position,

_{j}*d*- is the distance between the scanning slit and the imaging plane,

*δz*- is the scanning slit iteration step, and

*δ*- is the optical path difference across the slit. Then we determine the sampled wavefront phase profile - ${\varphi}_{N\mathit{exp}}^{i}$ for the harmonic wavelength -

_{j}*λ*across the scanning direction using ${\varphi}_{N\mathit{exp}}^{i}=\frac{2\pi}{{\lambda}_{N}}{\sum}_{j=0}^{i}{\delta}_{j}$. We determine the wavefront amplitude ${I}_{\mathit{exp}}^{i}$ at each sampling point by integrating intensity of the corresponding diffraction pattern. Further details of SWORD are discussed in Ref. [8].

_{N}The experimental wavefront phase {
${\varphi}_{N\mathit{exp}}^{i}$} sampled at positions {
${x}_{\mathit{exp}}^{i}$} was fit with parabolic function *ϕ _{N}* (

*x*) =

*a*(

_{N}*x*–

*x*

_{0})

^{2}for each harmonic. The Nonlinear Least Square method was used for the fit. To minimize the noise contribution, the phase is fit within the range where measured intensity is larger than 5% relative to the peak intensity. The wavefront intensity distribution { ${I}_{\mathit{exp}}^{i}$} measured at the scanning slit was fit to the Gaussian profile $I(x)={I}_{0}\text{exp}\left[-2{\left(x-{x}_{0}\right)}^{2}/{w}_{sN}^{2}\right]+{I}_{b}$, where

*I*

_{0}is maximum intensity,

*w*- Gaussian width at the slit position, and

_{sN}*I*- background signal. The measured parabola curvature -

_{b}*a*, and the Gaussian width -

_{N}*w*for each harmonic (H13-H25) were used to retrieve the virtual waist position

_{sN}*Z*and waist size

_{N}*w*

_{0}

*. A sensitivity analysis showed that the*

_{N}*Z*is highly sensitive to

_{N}*a*and depends weakly on

_{N}*w*(ex.: 10% change in the Gaussian width causes only ∼ 0.1% change in

_{sN}*Z*). We use the confidence bound of 95% in determination of the error bars for the relative waist position

_{N}*δZ*with respect to the scanning slit as shown in Fig. 2.

We have measured short trajectory harmonics. The divergence of the short trajectories mainly depends on the position of the jet relative to focal position of the fundamental beam. For all jet positions we find that the virtual origin of the harmonics differs by more than their confocal parameter across the spectrum.

The intensity distribution that we measure for different harmonics in Fig. 1(a) is Gaussian, within the ∼ 15% error bars caused by fluctuations in the laser power and gas jet density. We extract the Gaussian beam parameters by fitting the measured intensity profile for each harmonic.

The wavefront phase *ϕ _{N}* (

*x*,

*z*), where

_{slit}*z*is the slit position along the beam propagation and

_{slit}*x*is the coordinate across the slit, is parabolic for all harmonic orders

*N*(Fig. 1(a)). Because the noise in the phase measurements is small, the error bars are not visible on the scale of Fig. 1(a). To better estimate the phase error we subtract the fitted parabolic phase from the measured phase as shown in Fig. 1(b). Across the central part of the beam, containing 95% of the beam energy, the standard deviation of measured phase from ideal parabolic phase is 0.12rad. (better than

*λ*/50) for the 15th harmonic and 0.22rad (better than

*λ*/28) for the 21st harmonic.

Given the amplitude and the phase of the electric field of each harmonic at the scanning slit plane, we propagate the complex electromagnetic field backwards in time. Gaussian beams with parabolic wavefront phase remain Gaussian as they propagate. They can be characterized by their beam waist position *Z*_{N} (relatively to the scanning slit) and their beam waist size *w*_{0N}. Figure 2(a) shows the reconstructed waist positions as a function of the harmonic order. The differences in the position *Z*_{N} for different harmonics quantitatively show the different curvatures for each harmonic. The relative waist position *δZ* (Fig. 2(a)) varies by more than their Rayleigh range (*Z*_{0N}) which monotonically decrease from 11.8mm for H13 to 2.7mm for H25 (
${Z}_{0N}=\frac{\pi {w}_{0N}^{2}}{{\lambda}_{N}}$, where *λ _{N}* is the wavelength of

*N*harmonic). The corresponding beam (Fig. 2(b)) waists vary from 15.2

^{th}*μm*for H13 to 5.3

*μm*for H25.

## 3. Implications

The dependence of wavefront structure on harmonic order means that it is impossible to focus the attosecond pulse into one diffraction limited spot by achromatic optics (by using mirrors for example). Figure 3 shows how the harmonics would focus if an f=100mm mirror was placed at the scanning slit position (about 250mm from the jet).

With different harmonics focusing at different longitudinal positions it is impossible to maintain a Fourier-limited attosecond pulses as the beam propagates. As soon as a particular harmonic component passes through the focus - it acquires a *π* phase shift (the Gouy phase). Figure 3(f) shows the Gouy phase across the focal region for different harmonics.

Spectrally resolved wavefronts open the door for complete spatial-temporal reconstruction of attosecond pulses. Given the measured spectral wavefront (spectral wavefront phase - *ϕ _{N}* (

*r⃗,z*), and spectral wavefront amplitude - |

*E*(

*r⃗*,

*z*,

*ω*)|, where

_{N}*r⃗*- transverse coordinate across the beam,

*z*- is the coordinate along the beam propagation) at any plane

*z*=

*z*

_{0}along the beam propagation, it is enough to have temporal information at

*any single point r⃗*=

*r⃗*

_{0}across the beam to determine the temporal structure everywhere.

Since it is sufficient to measure the spectral phase *φ _{N}* (

*r⃗*

_{0},

*z*

_{0}) for any (

*r⃗*

_{0},

*z*

_{0}), we could select the point, for example, by centering a small aperture on high-harmonic beam and measuring the spectral phase of the transmitted radiation using RABBIT, FROG-CRAB [11, 12] or any other suitable technique. Then, the complete spatial (3D)+ temporal (1D) field profile of the beam can be reconstructed -

*ω*

_{0}is the angular frequency of the fundamental beam and the summation is over all emitted harmonic orders. Note, while the spectral wavefronts |

*E*(

*r⃗,z*,

*ω*)

_{N}*|*exp

*i*[

*ϕ*(

_{N}*r⃗*,

*z*)] is measured across one plane (

*z*=

*z*

_{0}), it can be found across any plane along the beam propagation using Fresnel-Kirchhoff diffraction formula [13] for each harmonic.

We apply this principle to our case, when the wavefront structure of each harmonic is taken from our measurements, but we arbitrary assume that the attosecond pulses within the train are transform limited at the center of the focusing mirror (Fig. 3(d)) i.e. *φ _{N}* =

*Nω*

_{0}

*t*, where

_{e}*t*is harmonic emission time. As we move along or across the beam, the relative spectral amplitude and phase of the harmonics changes. Consequently the temporal structure of each attosecond pulse in the pulse train depends on where the pulse is measured. Using the measured spectral wavefronts in our experiment, the temporal profiles is shown for selected points along the beam in Fig. 3(a)–3(c). Even in cases where the atto-chirp [14, 15] has been corrected, it can only be corrected strictly speaking at one point.

_{e}There are at least two approaches for bringing the harmonics to a common focus and generating Fourier limited attosecond pulses and pulse trains. (1)The phase of any harmonic depends on both the intensity and the phase of generating beam. Optimizing the spatial profile (both phases and amplitudes) of the fundamental beam will allow us to minimize the chromatic wavefront variations between different harmonics. (2)The chromatic aberrations that we measure are reproducible. They can, in principle, be corrected with specially designed linear optics.

## 4. Conclusion

In conclusion, within the signal-to-noise limits of our experiment, we have obtained complete (phase and amplitude) information on wavefront profile of each harmonic. We have shown that spectrally resolved wavefronts allow the field distribution to be found anywhere. If we know the position of the generating medium with sufficient accuracy, this includes the position at the generating medium itself. The wave front amplitude and phase in the generating medium isolate an imprint of the underlying single atom or molecule response. This will solve the long standing problem of intensity and phase focal volume averaging in attosecond science experiments.

Within a 3-step model [16], for a thin, low density gas, the measured wavefront phase *ϕ _{N}* (

*r⃗,z*) is given by the sum of four contributions:

_{jet}*ϕ*

_{fund}(

*r⃗*,

*z*) is the wavefront phase of the fundamental beam at the jet position and can be measured.

_{jet}*ϕ*(

_{r}*ω*,

_{N}*z*) is the phase contribution of the transition moment [17].

_{jet}*ϕ*(

_{r}*ω*,

_{N}*z*) is intensity independent for a single orbital, but intensity dependent for multiple orbitals. The phase that the electron acquires in the continuum is −|

_{jet}*α*|

_{N}*I*(

*r⃗*,

*z*). This phase can be calculated within the strong field approximation and its corrections [5, 6].

_{jet}*ϕ*(

_{t}*I*(

*r⃗*,

*z*),

_{jet}*ω*) is the yet not fully understood tunneling phase [18,19]. It might be intensity dependent. Therefore, wavefront structure will serve as a sensitive measurement of the intensity dependence of these parameters.

_{N}Once the relative phase of each harmonic is determined and spectrally resolved wavefronts are measured, then we have shown that we can resolve the temporal profile of an attosecond pulse anywhere in space. No matter what the pulse temporal structure should prove to be, we have already determined that there will be significantly different temporal profile at the center and at the edges of a beam and the temporal profile must change as the beam passes through the respective foci of the harmonics. Therefore, any spatially extended measurement with attosecond pulses will be affected and any measurement of the pulse itself is also affected. The result of any experiment will be effectively averaged over the different temporal profiles at different positions in the interaction region of the focal volume.

Finally, attosecond science has been restricted to spatially averaged measurements to image orbitals [10], tracing molecular dynamics [20], identify and time resolve tunneling wave packets [18] and to follow Auger decay [21]. Spectrally resolved wavefronts and the complete spatio-temporal characterization of attosecond pulses that they facilitate will allow us much greater experimental precision in all of these experiments, facilitating access to a single molecule response in the strong field.

We thank Mike Spanner, Misha Ivanov, Kyung Teac Kim and Eric Constant for stimulating discussions. We acknowledge financial support of Canada’s National Research Council, National Sciences, Engineering Research Council, the Canada Research Chair Program and MURI grant W911NF-07-1-0475. E.F acknowledges the support of Marie Curie International Outgoing Fellowship and H.N acknowledges the support of JST’s Presto program.

## References and links

**1. **M. Bellini, C. Lynga, A. Tozzi, M. B. Gaarde, T. W. Hänsch, A. L’Huillier, and C. G. Wahlstrom, “Temporal coherence of ultrashort high-order harmonic pulses,” Phys. Rev. Lett. **81**, 297–300 (1998). [CrossRef]

**2. **R. A. Bartels, A. Paul, H. Green, H. C. Kapteyn, M. M. Murnane, S. Backus, I. P. Christov, Y. W. Liu, D. Attwood, and C. Jacobsen, “Generation of spatially coherent light at extreme ultraviolet wavelengths,” Science **297**, 376–378 (2002). [PubMed]

**3. **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,” Opt. Lett. **28**, 480–482 (2003). [CrossRef]

**4. **C. Valentin, J. Gautier, J Goddet, C. Hauri, T. Marchenko, E. Papalazarou, G. Rey, S. Sebban, O. Scrick, P. Zeitoun, G. Dovillaire, X. Levecq, S. Bucourt, and M. Fajardo, “High-order harmonic wave fronts generated with controlled astigmatic infrared laser,” J. Opt. Soc. Am. B **25**, B161–B166 (2008). [CrossRef]

**5. **M. Lewenstein, P. Salières, and A. L’huillier, “Phase of the atomic polarization in high-order harmonic-generation,” Phys. Rev. A **52**, 4747–4754 (1995). [CrossRef]

**6. **M. B. Gaarde, F. Salin, E. Constant, P. Balcou, K. J. Schafer, K. C. Kulander, and A. L’Huillier, “Spatiotemporal separation of high harmonic radiation into two quantum path components,” Phys. Rev. A **59**, 1367–1373 (1999). [CrossRef]

**7. **P. Salières, A. L’Huillier, and M. Lewenstein, “Coherence control of high-order harmonics,” Phys. Rev. Lett. **74**, 3776–3779 (1995). [CrossRef]

**8. **E. Frumker, G. G. Paulus, H. Niikura, D. M. Villeneuve, and P. B. Corkum, “Frequency-resolved high-harmonic wavefront characterization,” Opt. Lett. **34**, 3026–3028 (2009). [CrossRef]

**9. **D. R. Austin, T. Witting, C. A. Arrell, F. Frank, A. S. Wyatt, J. P. Marangos, J. W. Tisch, and I. A. Walmsley, “Lateral shearing interferometry of high-harmonic wavefronts,” Opt. Lett. **36**, 1746–1748 (2011). [CrossRef]

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

**11. **P. M. Paul, E. S. Toma, P. Breger, G. Mullot, F. Auge, P. Balcou, H. G. Muller, and P. Agostini, “Observation of a train of attosecond pulses from high harmonic generation,” Science **292**, 1689–1692 (2001). [CrossRef]

**12. **Y. Mairesse and F. Quéré, “Frequency-resolved optical gating for complete reconstruction of attosecond bursts,” Phys. Rev. A **71**, 011401 (2005). [CrossRef]

**13. **J. W. Goodman, *Introduction to Fourier Optics* (McGraw-Hill, 2005).

**14. **Y. Mairesse, A. de Bohan, L. J. Frasinski, H. Merdji, L. C. Dinu, P. Monchicourt, P. Breger, M. Kovacev, R. Taieb, B. Carre, H. G. Muller, P. Agostini, and P. Salières, “Attosecond synchronization of high-harmonic soft x-rays,” Science **302**, 1540–1543 (2003). [CrossRef]

**15. **K. T. Kim, C. M. Kim, M. Baik, G. Umesh, and C. H. Nam, “Single sub-50-attosecond pulse generation from chirp-compensated harmonic radiation using material dispersion,” Phys. Rev. A **69**, 051805 (2004). [CrossRef]

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

**17. **A. T. Le, R. R. Lucchese, S. Tonzani, T. Morishita, and C. D. Lin, “Quantitative rescattering theory for high-order harmonic generation from molecules,” Phys. Rev. A **80**, 013401 (2009). [CrossRef]

**18. **O. Smirnova, Y. Mairesse, S. Patchkovskii, N. Dudovich, D. M. Villeneuve, P. B. Corkum, and M. Y. Ivanov, “High harmonic interferometry of multi-electron dynamics in molecules,” Nature **460**, 972–977 (2009). [CrossRef]

**19. **B. K. McFarland, J. P. Farrell, P. H. Bucksbaum, and M. Guhr, “High harmonic generation from multiple orbitals in N-2,” Science **322**, 1232–1235 (2008). [CrossRef]

**20. **W. Li, X. Zhou, R. Lock, S. Patchkovskii, A. Stolow, H. C. Kapteyn, and M. M. Murnane, “Time-resolved dynamics in N_{2}O_{4} probed using high harmonic generation,” Science **322**, 1207–1211 (2008). [CrossRef]

**21. **M. Drescher, M. Hentschel, R. Kienberger, M. Uiberacker, V. Yakovlev, A. Scrinzi, T. Westerwalbesloh, U. Kleineberg, U. Heinzmann, and F. Krausz, “Time-resolved atomic inner-shell spectroscopy,” Nature **419**, 803–807 (2002). [CrossRef]