## New Research In

### Physical Sciences

### Social Sciences

#### Featured Portals

#### Articles by Topic

### Biological Sciences

#### Featured Portals

#### Articles by Topic

- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology

# Spatiotemporal coupling of attosecond pulses

Contributed by Anne L’Huillier, January 18, 2019 (sent for review October 12, 2018; reviewed by David Attwood and Margaret M. Murnane)

### See related content:

- QnAs with Anne L’Huillier- Mar 01, 2019

## Significance

In most optics textbooks, one writes the electric field describing an optical wave as a product of temporal and spatial amplitudes. This approximation often breaks down for short optical pulses. An example of such spatiotemporal coupling is chromatic aberrations, where the focal properties of the radiation vary with frequency over the pulse bandwidth. In this work, we point out significant chromatic aberrations of attosecond pulses, which depend on the geometry of the generation process. These aberrations are intrinsic to the generation process and need to be eliminated in applications requiring attosecond pulses to be focused over a small region.

## Abstract

The shortest light pulses produced to date are of the order of a few tens of attoseconds, with central frequencies in the extreme UV range and bandwidths exceeding tens of electronvolts. They are often produced as a train of pulses separated by half the driving laser period, leading in the frequency domain to a spectrum of high, odd-order harmonics. As light pulses become shorter and more spectrally wide, the widely used approximation consisting of writing the optical waveform as a product of temporal and spatial amplitudes does not apply anymore. Here, we investigate the interplay of temporal and spatial properties of attosecond pulses. We show that the divergence and focus position of the generated harmonics often strongly depend on their frequency, leading to strong chromatic aberrations of the broadband attosecond pulses. Our argument uses a simple analytical model based on Gaussian optics, numerical propagation calculations, and experimental harmonic divergence measurements. This effect needs to be considered for future applications requiring high-quality focusing while retaining the broadband/ultrashort characteristics of the radiation.

- attosecond pulse
- high-order harmonic generation
- Gaussian optics
- spatiotemporal coupling
- focusing of XUV radiation

Electromagnetic waves are usually mathematically described by a product of purely spatial and purely temporal terms. This approximation often fails for broadband femtosecond laser pulses (ref. 1 and references therein), and spatiotemporal couplings need to be considered. Spatiotemporal couplings for visible or infrared (IR) light may be introduced by refractive and dispersive elements, such as lenses, gratings, or prisms. The noncollinear amplification in optical parametric crystals may also potentially lead to spatiotemporal couplings, and it is important to develop characterization methods to measure and reduce their effects (2⇓–4). In some cases, these couplings may be advantageously used, as, for example, demonstrated by Vincenti and Quéré (5) for the so-called “lighthouse” effect (6, 7).

The shortest light pulses, generated by high-order harmonic generation (HHG) in gases, are in the extreme UV (XUV)/soft X-ray region and in the range of 100 as (8⇓⇓–11), with bandwidths of a few tens or even hundreds of electronvolts (12, 13). These pulses are generated in a three-step process (14, 15). When an atom is exposed to a strong laser field, an electron in the ground state can tunnel through the atomic potential bent by the laser field, propagate in the continuum, and recombine back to the ground state when (and if) returning close to the ionic core. In this process, an XUV photon is emitted, with energy equal to the ionization energy plus the electron kinetic energy at return. Two main families of trajectories leading to the same photon energy can be identified. They are characterized by the “short” or “long” time of travel of the electron in the continuum (16, 17). Interferences of attosecond pulses emitted at each laser half-cycle leads to a spectrum of odd-order harmonics.

The investigation of spatiotemporal coupling of attosecond pulses requires measurements of their spatial properties, as a function of time or, equivalently, frequency. Wavefronts of high-order harmonics have been measured by several groups, using different techniques such as Spectral Wavefront Optical Reconstruction by Diffraction (18⇓–20), lateral shearing interferometry (21), point-diffraction interferometry (22), and Hartmann diffraction masks (23, 24). In particular, Frumker et al. (25) pointed out that the variation of wavefront and intensity profile with harmonic order leads to spatiotemporal coupling of the attosecond pulses, with temporal properties depending on where they are measured.

The spatial and spectral properties of high-order harmonics strongly depend on the geometry of the interaction and, in particular, on whether the gas medium in which the harmonics are generated is located before or after the focus of the driving laser beam (26). The asymmetry between “before” and “after” can be traced back to the phase of the emitted radiation, which is equal to that of the incident laser field multiplied by the process order, as in any frequency upconversion process, plus the dipole phase which is accumulated during the generation and mostly originates from electron propagation in the continuum. While the former is usually antisymmetric relative to the laser focus, the latter depends on the laser intensity and is therefore symmetric (21, 27). The total phase and thus the divergence properties are different before and after the laser focus, leading to a strong dependence of the spatiotemporal properties of the harmonic radiation on the generation conditions. In some conditions, harmonics can be emitted with a flat wavefront (21) or even as a converging beam (28, 29). Another phenomenon leading to an asymmetry of HHG with respect to the generation conditions is ionization-induced reshaping of the fundamental field, which depends on whether the beam is converging or diverging when entering the gas medium (20, 30⇓–32).

In the present work, we show that the frequency components of attosecond pulses generated by HHG in gases have different divergence properties, which depend on the geometry of the interaction and in particular on where the generating medium is located relative to the laser focus. In some conditions, the position of the focus and divergence strongly vary with frequency, leading to chromatic aberrations, as sketched in Fig. 1, similar to the effect that a chromatic lens has on broadband radiation (33, 34). Any imaging optical component will focus the frequency components of the attosecond pulses at different locations, resulting in spatiotemporal couplings. Depending on the position where the pulses are characterized or used, they will have different central frequencies, pulse durations, and spatial widths. We use an analytical expression for the dipole phase (35) combined with traditional Gaussian optics to predict the radius of curvature, position of focus, and divergence of the two trajectory contributions to HHG (29). This model, which assumes generation in a thin slab (36, 37), is validated by using numerical simulations of HHG (38) for both thin and thick generating media. We also present experimental measurements of the harmonic divergence as a function of position of generation relative to the laser focus. Finally, we discuss the implications of our results for the focusing of broadband attosecond pulses.

## Analytical Expression of the Dipole Phase

The single-atom response of HHG is well described by an approximate solution of the time-dependent Schrödinger equation (TDSE) for an atom in a strong laser field, called the strong-field approximation (SFA) (39). This theory leads to a simple analytical expression of the dipole phase, equal to

In this approximation, the second step of the process is described by solving Newton’s equation of motion for a free particle in the laser field. Fig. 2 shows the frequency (Ω) of the emitted XUV radiation as a function of electron return time for two different fundamental field intensities, indicated by the solid and dashed curves. The frequency varies from *s*, ℓ refers to the electron trajectory (short or long), and **3** do not depend on laser intensity. Using ^{−2}, ^{−2}, ^{−1}⋅cm^{2}.

The dipole phase can be approximated for the two families of trajectories by the expansion:

## Wavefront and Spatial Width of XUV Radiation

We now use this analytical expression for the dipole phase together with traditional Gaussian optics to predict the radius of curvature, position of focus, and divergence of the two trajectory contributions to HHG. A similar derivation has been proposed, independently, by Quintard et al. (29) with, however, a different analytical formulation of the dipole phase. We neglect the influence of propagation, considering an infinitely thin homogeneous gas medium (36, 37, 44). Such an approximation is valid in a loose focusing geometry, where the generating medium length is much smaller than the Rayleigh length. We also assume that the fundamental field is Gaussian, with intensity *q*th harmonic field can be approximated by**6**, for **6**, which does not depend on intensity and therefore on space,

We use a Taylor expansion close to the center of the beam to approximate **8**). To determine the harmonic wavefront, we only keep the terms proportional to **8**, to which we add the *i*), intensity at focus (**9** is illustrated in Fig. 3*A*, representing the wavefronts induced by the fundamental (black) and due to the dipole phase for the short trajectory (green) as a function of the generation position. The fundamental wavefront changes from convergent to divergent through the focus, while that induced by the dipole phase is always divergent and independent of the generation position (z).

Using the reduced coordinate **9** can be written as*B* for the 23rd harmonic of 800-nm radiation generated in Ar, with *I*_{0} = 3 × 10^{14} W⋅cm^{−2} . In these conditions, we have *B* presents the radius of curvature in reduced units **11** is

To estimate in a simple way the spatial width of the harmonic field at the generation position, we assume that its amplitude is proportional to the fundamental amplitude to a power p (36, 37, 44, 46, 47). This exponent is quite constant in the plateau region and typically of the order of 4, as confirmed by our TDSE calculations presented below. The harmonic width is then simply equal to

## Focus Position and Beam Waist

Knowing the beam radius of curvature and width at a given position z, it is a simple exercise within Gaussian optics to determine the position of the focus and the corresponding waist (e.g., ref. 45). The position of focus relative to the generation position z is given by**13** can be written as

Fig. 4 shows the position of the harmonic focus (*A*) and the normalized far-field divergence *B*) for the two trajectories, short (blue solid line) and long (red solid line). The color plots represent harmonic intensities obtained from a simulation presented in *Numerical Calculations*. The divergence of the fundamental *A*, zero on the vertical scale means that the focus of the harmonic field coincides with the generation position. The focus position and divergence strongly vary with z and quite differently for the two trajectories. In both cases, the focus position changes sign, and the divergence goes through a minimum when the radius of curvature goes to infinity (Fig. 3).

For the short trajectory and

The same reasoning applies for the long trajectory contribution, except that ^{−2}, the harmonic focus is located just before the generation position, and the divergence is much larger than that of the short trajectory contribution.

At the positions where the radius of curvature diverges (indicated by the dashed line in Fig. 4 for the short trajectory), the harmonics are generated with a flat wavefront and with a large focus (low divergence). In contrast, harmonics generated far away from the divergence minima will inherit a curvature from the fundamental and the dipole phase contribution which corresponds to a significantly smaller beam waist in the real or virtual focus and thus in a significantly larger divergence. The variation of the divergence with generation position is due partly to the dipole phase contribution, but also to the mismatch between the harmonic order q and the amplitude variation here described by a power law with exponent **16**).

## Numerical Calculations

To validate the Gaussian model presented in this work, we performed calculations based on single-atom data obtained by solving the TDSE for a single active electron in Ar exposed to a constant intensity. The time-dependent dipole response was calculated for 5,000 intensity points. This allows us, for each harmonic frequency, to precisely unwrap the amplitude and phase variation as a function of intensity, and thus to accurately describe the interferences of the trajectories. The complex electric-field distribution at a given harmonic frequency is obtained by integrating in time the polarization induced by the fundamental field in an arbitrarily thin sheet of homogeneous Ar gas. The field is then propagated to different positions relative to the generation position by calculating the diffraction integral in Fresnel approximation using Hankel transforms. The influence of ionization is not taken into account. This procedure is repeated for different gas target positions relative to the laser focus. We use a fundamental wavelength of 800 nm, a pulse duration of 45 fs, and a peak intensity of ^{−2}.

Fig. 4*A* presents a color plot of the 23rd harmonic on-axis intensity for different generation positions (horizontal axis). The regions with the warmest colors (i.e., toward red) represent the focal regions. The small regions with high peak intensity (dark red, like that labeled II) correspond to the smallest focus. The agreement between the numerical predictions and those of the Gaussian model is striking. When

Fig. 4*A* presents a series of interference structures, some vertical, others almost horizontal. To identify the physical reason for these structures, we have performed simulations allowing us to separate the contributions of the trajectories, using the thin medium approximation and harmonic fields as in our model. Instead of Gaussian optics, however, we used diffraction integrals for the propagation. These simulations show that the horizontal fringes (e.g., between regions II and III) concern the short trajectory contribution and come from the fact that the harmonic phase front and intensity profile are not those of a Gaussian beam (37). The vertical features (e.g., between I and II), however, are a manifestation of quantum path interferences (41, 42), since they only appear when both contributions are coherently added.

The color plot in Fig. 4*B* is the 23rd harmonic radial intensity at a distance of *B* can be attributed to the long trajectory, while the red–yellow–bright green regions are due to the short trajectory.

An important question is whether these results are still valid after propagation in a finite medium. We used the single-atom data described as input in a propagation code based on the slowly varying envelope and paraxial approximations (38). We present in Fig. 5 results obtained for a 5.4-mm-long (*A*), 30-mm-long (*B*), and 60-mm-long (*C*) homogeneous medium, using a 2-mbar gas pressure and a fundamental waist size of *A* compares very well with the results shown in Fig. 4*A*, as expected, Fig. 5 *B* and *C* shows clear effects of propagation, related to ionization-induced defocusing of the fundamental laser beam. In fact, two different phase-matching regimes appear: one similar to what is present in absence of propagation and which agrees well with the predictions of the Gaussian model (compare regions I and III in Fig. 5 *A* and *B*), and a second one, which also follows a similar model but for a fundamental focus moved to the left (see regions I′ and III′ in Fig. 5*B*), as expected for a fundamental beam that is defocused due to partial ionization of the medium (20, 30⇓–32). To examine in more details the effect of propagation goes beyond the scope of this paper.

## Experimental Divergence Measurements

Experiments were performed at the intense XUV beamline of the Lund Laser Centre (48, 49), by using a multiterawatt 45-fs titanium–sapphire laser operating at a 10-Hz repetition rate. The beam was (slightly) apertured to 27 mm and focused by using a spherical mirror with focal length

The widths of the 13th–19th harmonics are shown in Fig. 6*A* and compared with theoretical predictions in *B*–*D*, obtained by using a laser waist of 220 μm and a maximum intensity of ^{−2} (the Rayleigh length is estimated to 0.2 m). The harmonic widths were calculated as *B* presents results of numerical calculations based on solving the TDSE and the propagation equations, using parameters mimicking the experimental conditions as well as possible. In Fig. 6*C*, the results of the Gaussian model for the short trajectory are shown, while in Fig. 6*D*, a “truncated” Gaussian model is used, where the expressions for the beam waist, radius of curvature, and intensity variation of the fundamental beam now include the effect of a circular aperture (50, 51), taken to be equal to the experimental one. Going from the left to the right in all of the plots in Fig. 6, the harmonic widths first decrease (or stay approximately constant for the highest orders) and then increase. The harmonic widths vary more strongly in the Gaussian model than in the other calculations and in the experiment. We investigated the reason for this difference by varying the parameters used in the propagation simulations, such as medium length, gas pressure, aperture diameter, and pulse energy. Unlike the conditions used in Fig. 5 *B* and *C*, effects due to propagation, e.g., induced by ionization-defocusing, are negligible, and the main reason for the difference between Fig. 6 *B* and *C* is the beam truncation due to the aperture, as confirmed by Fig. 6*D*. Effects due to propagation in the nonlinear medium, which become nonnegligible at higher laser intensity, actually lead to faster variation of the beam divergence on both sides of the laser focus.

## Chromatic Aberrations of Attosecond Pulses

Finally, we study the variation of the focus position and beam waist over a large spectral bandwidth. To obtain a broad spectral region, we consider generation of high-order harmonics in neon atoms. HHG spectra obtained in Ne (52) are broader and flatter than those in Ar, which exhibit a strong modulation due to a Cooper minimum at ∼45 eV. Fig. 7 shows the predictions of the Gaussian model, for the 31st to the 71st harmonics of 800-nm radiation, at an intensity of 5^{−2}. We only consider here the contribution from the short trajectory. The Gaussian model is used here for simplicity. It should be reasonably accurate in conditions of thin medium, not too high degree of ionization, which is the case for Ne atoms at the intensity used, and high spatial quality, nontruncated, fundamental beam.

The variation of the focus position as a function of generation position strongly depends on the process order. This is due to the frequency dependence of Eq. **9** and, in particular, depends on whether the radius of curvature diverges. Since **9** move closer to each other, as is clear in Fig. 7*A* by comparing, e.g., harmonics 41 and 43 (*A*, *B*.

To estimate the consequence of these spatial properties on the spectral characteristics of the attosecond pulses (25), we examine the variation of the on-axis spectrum at different positions (labelled ① to ④), for the generation position indicated by the dashed line. This is equivalent to examining the properties of the generated radiation after refocusing as illustrated in Fig. 1, for different “detection positions” in the focal region. We here assume equal strength of the generated harmonics, but account for the frequency variation in beam waist size (Fig. 7*B*) and position (Fig. 7*A*). The harmonic spectra shown in Fig. 7*C* are found to be strongly dependent on the detection position, with, in some cases, strong bandwidth reduction and displacement of the central frequency. In contrast, the on-axis harmonic phase is dominated by the attosecond chirp, given by the last term in Eq. **6** and indicated by the black line in Fig. 7*C*, and does not vary much with the detection position.

## Spatiotemporal Coupling of Attosecond Pulses

Finally, we estimate the influence of the chromatic aberrations on the temporal properties of the attosecond pulses. We consider a flat spectrum between harmonics 31 and 71 at the generation position indicated by the dashed line in Fig. 7, in the same conditions as in the previous paragraph. We propagate the harmonic fields using diffraction integrals and coherently add them to obtain the resulting attosecond pulse train in space and time at different detection positions. We take into account the different focus positions and divergences of the frequency components of the attosecond pulses, as well as the attosecond positive chirp.

Fig. 8 shows the spectral (*A*) and temporal (*B*) intensity (in color) of the generated attosecond pulse on-axis as function of the detection position relative to the generation position, here equal to *A*). The highest intensity is obtained before the generation position, while the shortest pulse is obtained afterward, as follows from Fig. 8*B*. The attosecond pulse is not the shortest at the generation position, where the spectral bandwidth is the largest, because the attosecond chirp stretches the pulse in time. Fig. 8 *A* and *B* strikingly show that the shortest pulse and the highest intensity of the attosecond pulse are obtained in different positions, illustrating the difficulty of refocusing high-order harmonics, particularly for applications requiring high intensity.

Finally, Fig. 8 *C* and *D* shows the spatiotemporal intensity profiles of the attosecond pulse at the positions where is it spectrally broadest (*C*, ②) and where it is most intense (*D*, ③). The difference between the two cases is a signature of the strong spatiotemporal couplings of the generated attosecond pulses. These couplings, here studied at

## Conclusion

In this work, we examine the focusing properties of high-order harmonics generated in gases. We use a simple Gaussian optics model, valid for a thin generating medium, assuming a fundamental Gaussian beam, and based on an analytical expression of the frequency- and intensity-dependent dipole phase. This model allows us to predict the focus and divergence of the two trajectory contributions to HHG. We validate the predictions of the model by numerical calculations based on solving the TDSE for the single-atom response and propagation equations for the response from the macroscopic medium. Experimental divergence measurements performed at the intense XUV beamline of the Lund Laser Centre show similar trends as those predicted by the numerical calculations, as well as by an extension of the Gaussian model, which includes the effect of a circular aperture. We also discuss the consequences of the fact that the harmonics have different focus positions and beam waists on the resulting spectra and pulse durations. The relative harmonic amplitudes are found to vary with the detection position, thus strongly affecting the spatiotemporal properties of the corresponding attosecond pulses.

The effects investigated in the present work have a strong impact on applications of attosecond pulses, requiring a small focal spot (e.g., to reach a high XUV intensity) over a broad bandwidth or during a short (attosecond) duration. These spatiotemporal couplings may be reduced by locating the generation medium after the laser focus and/or by minimizing the influence of the dipole phase, using a shaped fundamental beam (53, 54) or generating in waveguides (capillaries) (47, 55).

## Acknowledgments

This research was supported by the Swedish Research Council; the Swedish Foundation for Strategic Research; European Research Council Grant 339253 PALP; the Knut and Alice Wallenberg Foundation; and National Science Foundation Grant PHY-1713761. This project received funding from the European Union’s Horizon 2020 research and innovation program under Marie Skłodowska-Curie Grant Agreements 641789 MEDEA and 793604 ATTOPIE.

## Footnotes

↵

^{1}H.W., C.G., and J.V. contributed equally to this work.↵

^{2}Present address: Max Born Institute, 12489 Berlin, Germany.- ↵
^{3}To whom correspondence should be addressed. Email: anne.lhuillier{at}fysik.lth.se.

This contribution is part of the special series of Inaugural Articles by members of the National Academy of Sciences elected in 2018.

Author contributions: H.W., C.G., J.V., P.W.S., H.C.-A., P.R., P.J., C.L.A., and A.L. designed research; H.W., C.G., J.V., P.W.S., H.C.-A., J.L., J.P., P.R., H.D., S.C., S.M., M.B.G., P.J., C.L.A., and A.L. performed research; C.G., J.V., P.W.S., H.C.-A., S.C., S.M., and M.B.G. contributed with analytic tools; H.W., C.G., J.V., P.W.S., H.C.-A., J.P., P.R., S.M., and C.L.A. analyzed data; and H.W. and A.L. wrote the paper.

Reviewers: D.A., University of California, Berkeley; and M.M.M., University of Colorado and NIST.

The authors declare no conflict of interest.

See QnAs on page 4767

- Copyright © 2019 the Author(s). Published by PNAS.

This open access article is distributed under Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND).

## References

- ↵
- ↵
- Miranda M, et al.

- ↵
- Pariente G,
- Gallet V,
- Borot A,
- Gobert O,
- Quéré F

- ↵
- Harth A, et al.

- ↵
- ↵
- Kim KT, et al.

- ↵
- Louisy M, et al.

- ↵
- Sansone G, et al.

- ↵
- Goulielmakis E, et al.

- ↵
- Li J, et al.

- ↵
- Gaumnitz T, et al.

- ↵
- Popmintchev T, et al.

- ↵
- Cousin SL, et al.

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Lloyd DT,
- O’Keeffe K,
- Anderson PN,
- Hooker SM

- ↵
- Johnson AS, et al.

- ↵
- ↵
- ↵
- Valentin C, et al.

- ↵
- Freisem L,
- Jansen GSM,
- Rudolf D,
- Eikema KSE,
- Witte S

- ↵
- ↵
- ↵
- ↵
- Quintard L

- ↵
- Quintard L, et al.

- ↵
- ↵
- Tamaki Y,
- Itatini J,
- Nagata Y,
- Obara M,
- Midorikawa K

- ↵
- ↵
- ↵
- Bor Z,
- Horváth ZL

- ↵
- Guo C, et al.

- ↵
- Hernández-García C,
- San Román J,
- Plaja L,
- Picón A

- ↵
- Catoire F, et al.

- ↵
- ↵
- ↵
- Varjú K, et al.

- ↵
- Carlström S, et al.

- ↵
- ↵
- Guo C

- ↵
- Rego L,
- San Román J,
- Picón A,
- Plaja L,
- Hernández-García C

- ↵
- Saleh BEA,
- Teich MC

- ↵
- ↵
- ↵
- Manschwetus B, et al.

- ↵
- Coudert-Alteirac H, et al.

- ↵
- Horváth ZL,
- Bor Z

- ↵
- ↵
- ↵
- Boutu W, et al.

- ↵
- ↵
- Popmintchev T,
- Chen MC,
- Arpin P,
- Murnane MM,
- Kapteyn HC

## Citation Manager Formats

## Sign up for Article Alerts

## Article Classifications

- Physical Sciences
- Physics

## Jump to section

- Article
- Abstract
- Analytical Expression of the Dipole Phase
- Wavefront and Spatial Width of XUV Radiation
- Focus Position and Beam Waist
- Numerical Calculations
- Experimental Divergence Measurements
- Chromatic Aberrations of Attosecond Pulses
- Spatiotemporal Coupling of Attosecond Pulses
- Conclusion
- Acknowledgments
- Footnotes
- References

- Figures & SI
- Info & Metrics