# Anisotropic dynamics of two-photon ionization: An attosecond movie of photoemission

Alice Autuori, Dominique Platzer, Mariusz Lejman, Guillaume Gallician, Lucie Maëder, Antoine Covolo, Lea Bosse, Malay Dalui, David Bresteau, Jean-François Hergott, et al.

► To cite this version:

Alice Autuori, Dominique Platzer, Mariusz Lejman, Guillaume Gallician, Lucie Maëder, et al. Anisotropic dynamics of two-photon ionization: An attosecond movie of photoemission. Science Advances, 2022, 8 (12), pp.eab17594. 10.1126/sciadv.abl7594. hal-03629043

HAL Id: hal-03629043

https://hal.science/hal-03629043

Submitted on 4 Apr 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Anisotropic dynamics of two-photon ionization: An attosecond movie of photoemission

Alice Autuori<fnref n="1" />†, Dominique Platzer<fnref n="1" />†, Mariusz Lejman<fnref n="1" />, Guillaume Gallician<fnref n="1" />, Lucie Maëder<fnref n="1" />, Antoine Covolo<fnref n="1" />, Lea Bosse<fnref n="1" />, Malay Dalui<fnref n="1" />, David Bresteau<fnref n="1" />, Jean-François Hergott<fnref n="1" />, Olivier Tcherbakoff<fnref n="1" />, Hugo J. B. Marroux<fnref n="1" />, Vincent Loriot², Franck Lépine², Lionel Poisson<fnref n="1" />,³, Richard Taëib⁴, Jérémie Caillat⁴*, Pascal Salières<fnref n="1" />*

Imaging in real time the complete dynamics of a process as fundamental as photoemission has long been out of reach because of the difficulty of combining attosecond temporal resolution with fine spectral and angular resolutions. Here, we achieve full decoding of the intricate angle-dependent dynamics of a photoemission process in helium, spectrally and anisotropically structured by two-photon transitions through intermediate bound states. Using spectrally and angularly resolved attosecond electron interferometry, we characterize the complex-valued transition probability amplitude toward the photoelectron quantum state. This allows reconstructing in space, time, and energy the complete formation of the photoionized wave packet.

Copyright © 2022
The Authors, some
rights reserved;
exclusive licensee
American Association
for the Advancement
of Science. No claim to
original U.S. Government
Works. Distributed
under a Creative
Commons Attribution
License 4.0 (CC BY).

## INTRODUCTION

Photoemission has played a key role in the development of quantum physics theory. Besides its fundamental aspects, it has become a widespread technique for the investigation of matter, from chemical analysis in condensed matter to electronic structure determination in the gas phase. The development of bright sources of extreme ultraviolet (XUV) radiation such as synchrotrons has allowed the measurement of angular distributions of photoelectrons with high precision both in the laboratory and molecular frames. Pioneering complete experiments (1, 2) gave access to the relative weights and phases of the partial waves contributing to the wave function of the liberated electron. These measurements, however, do not give access to photoemission in the time domain because they are performed independently for each final electron energy and thus lack the relative phase between the spectral components.

The advent of broadband coherent XUV sources based on high-order harmonic generation (HHG) from intense infrared (IR) pulses (3, 4) has triggered the development of unique interferometric schemes to measure spectral phases, leading to unprecedented insight into the temporal dynamics of elementary processes in a broad range of chemical species with attosecond resolution (1 as = $10^{-18}$ s). Both the streaking (5) and the Reconstruction of Attosecond Beating By Interference of Two-Photon Transitions (RABBIT) technique (6–9) were applied to revisit photoemission in the time domain, through the first measurements of attosecond photoemission delays in atoms (10, 11), molecules (12, 13), nanoparticles (14), liquids (15), and solids (16). Access to higher levels of details was gradually achieved using tunable sources (17–21) or by increasing the spectral resolution at detection with the Rainbow RABBIT technique (22–25). However, all the aforementioned studies relied on measurements averaged

over the direction of photoemission. While attosecond photoelectron interferometry has been investigated in momentum space in earlier experiments, e.g., in (26, 27), including the pioneering study (28), it is only recently that orientation-resolved spectral phase measurements could be performed, using cold target recoil ion momentum spectroscopy (COLTRIMS) (29–34) or velocity-map imaging spectroscopy (VMIS) (35–38). Resonant photoemission dynamics was studied angularly with RABBIT in the vicinity of, e.g., autoionizing states (30) or shape resonances (34). Using a two-color above-threshold interferometric scheme, autoionizing dynamics in chiral molecules could be angularly resolved in the laboratory frame (35).

In the present study, we combine attosecond spectral interferometry with momentum spectroscopy to record the modulus and phase variations of the photoelectron quantum state with high spectral resolution and angular sensitivity. The potential of this complete quantum phase spectroscopy is demonstrated in the test case of two-photon XUV + IR photoionization of helium through the intermediate resonant states 1s3p and 1s4p. The resulting structured photoelectron wave packet is fully characterized by measuring quasi-continuously the spectral and spatial variations of the final quantum state over a 0.8-eV spectral range. This allows reconstructing the attosecond photoemission dynamics strongly affected by the sudden phase jumps of up to $\pi$ rad measured in both dimensions.

## RESULTS

The concept of the technique is illustrated in Fig. 1A, the corresponding experimental setup being detailed in the Supplementary Materials. We aim at fully characterizing the resonant (15 + 1) two-photon transition in He using a reference (17 – 1) transition. Their interference determines the photoelectron angular and spectral distribution in the sideband peak (SB16) that is measured by VMIS (see Fig. 1B). The SB16 intensity $I_{16}$ thus depends on three parameters: (i) the XUV-IR delay $\tau$, (ii) the photoelectron energy $E$, and (iii) the angle $\theta$ between the electron emission and the shared XUV-IR polarization axis $z$, in the generic form

$$I_{16}(\tau; E, \theta) = A_{16}(E, \theta) + B_{16}(E, \theta) \cos [2 \omega_0 \tau - C_{16}(E, \theta)] \quad (1)$$

<footnote n="1">¹Université Paris-Saclay, CEA, CNRS, LIDYL,91191 Gif-sur-Yvette, France. ²Université de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, 69622 Villeurbanne, France. ³Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay,91405 Orsay, France. ⁴Sorbonne Université, CNRS, Laboratoire de Chimie Physique-Matière et Rayonnement, 75005 Paris, France.</footnote>

<footnote>*Corresponding author. Email: pascal.salieres@cea.fr (P.S.); jeremie.caillat@sorbonne-universite.fr (J.C.)</footnote>

<footnote>†These authors contributed equally to this work.</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Principle of spectrally and angularly resolved photoelectron interferometry. (A) Helium is ionized with a comb of coherent odd harmonics of a titanium:sapphire laser combined with a weak fraction of the fundamental laser. The two-photon electron wave packet created in the continuum by absorption of H15 and an IR photon is highly structured because of the intermediate 1s3p and 1s4p resonant states. This wave packet, (15 + 1), is a coherent superposition of the s and d partial waves that are populated by the two-photon transition. A reference wave packet, (17 – 1), is created at the same energy by absorption of H17 and stimulated emission of an IR photon. (B) Momentum map of the delay-integrated SB16 from an Abel-inverted VMIS measurement. It exhibits two main spectral components centered around 1s3p + IR and 1s4p + IR.</figcaption></figure>

where $\omega_0$ is the laser frequency. Our study focuses on the amplitude $B_{16}$ and phase $C_{16}$, of which the complex anisotropic photoemission dynamics is decoded. Both quantities are accessed through a Fourier transform of the three-dimensional (3D) spectrogram $I_{16}(\tau; E, \theta)$ with respect to $\tau$ at each sampled energy $E$. They are then calibrated to extract the intrinsic two-photon transition amplitude $M(E, \theta)$ associated with the probed (15 + 1) path independently of the characteristics of the XUV exciting pulses (see the Supplementary Materials).

The modulus and phase of $M(E, \theta)$ obtained from the experimental data are shown in Fig. 2 (A and B). They display many structures related to the buildup of the wave packet through intermediate resonances. At $E \approx E_{1snp} + \hbar\omega_0$, the modulus is enhanced and a $\approx\pi$ rad smooth spectral phase drop occurs, reminiscent of the observations of (17). In contrast, between resonances, the signal substantially drops and a sharp $\approx\pi$ rad spectral phase jump occurs. A most remarkable feature is the strong angular dependence of both modulus and phase (jumps up to $\pi$ rad) over the covered spectral range. This is highlighted by fast changes in shape and phase of the polar representation of $M(E, \theta)$, as shown in Fig. 2C for few selected representative energies.

To support the analysis of our measurements, we performed numerical simulations based on the time-dependent Schrödinger equation (TDSE) with a model potential for helium and pulse

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Probability amplitude associated with the investigated (15 + 1) two-photon transition. Spectral and angular variations of the modulus and phase of $M(E, \theta)$ resulting from a Rainbow RABBIT analysis of the measured (A and B, respectively) and simulated (D and E, respectively) spectrogram $I_{16}(\tau; E, \theta)$. The phases are shown with a negative sign for readability. Polar representation of $M(E, \theta)$ at specific energies labeled from (i) to (viii), where the $\theta$-dependent radius and color correspond respectively to the modulus and phase from experimental measurements (C) and from TDSE simulations (F).</figcaption></figure>
characteristics corresponding to the experimental ones. The results, shown in Fig. 2 (D and E), are in good agreement with the experimental ones, the latter showing slightly attenuated phase jumps that are down-shifted by ~80 meV around E = 0.32 eV and a stronger angular dependence around the 1s4p resonance.

DISCUSSION

The observed rich spectral and angular features are the signature of a structured anisotropic photoemission dynamics, fully encoded in the differential transition probability amplitudes M(E, θ). They can be understood through a partial-wave expansion of the latter

$$M(E, \theta) = M_0(E) Y_{00}(\theta) + M_2(E) Y_{20}(\theta) \quad (2)$$

where $Y_{\ell 0}$ are the spherical harmonics and $M_{\ell}(E)$ are the matrix elements associated with the final angular momenta $\ell$.

For monochromatic fields, in our spectral range where the XUV photon frequency $\Omega$ is nearly resonant with the 1snp states, the partial amplitudes can be approximated by (39)

$$M_{\ell}(E) \approx \lim_{\epsilon \to 0} \sum_{n=3}^{5} \frac{\langle 1sE\ell | z | 1snp \rangle \langle 1snp | z | 1s^2 \rangle}{E_{1s^2} - E_{1snp} + \hbar\Omega + i\epsilon} \quad (3)$$

where $|1s^2\rangle$ is the initial ground state of energy $E_{1s^2}$, $|1sE\ell\rangle$ is the considered final partial wave with photoelectron energy $E = E_{1s^2}$ +

$\hbar\Omega + \hbar\omega_0$. In our experiments and simulations, the finite duration/ bandwidth of the XUV and IR fields results in a smoothing of the outcomes from Eq. 3 (40, 41), which accounts for the spectral variations of $M_{\ell}(E)$ shown in Fig. 3 (C and D).

In contrast to the structureless (17 – 1) reference path, dominated by the isotropic s wave, the studied transition matrix element for each of the $\ell$ channels undergoes a series of $\pi$ rad phase jumps accompanied by a strong modulus increase, around 0.06- and 0.71-eV photoelectron energies. They are reminiscent of a zero crossing of the denominator in Eq. 3, wherever $\Omega$ is resonant with one of the 1snp ($n$ = 3 and 4) levels. Furthermore, between consecutive resonances, two contributions with opposite signs dominate the sum and cancel each other at a given energy, depending on the relative values of the $\ell$-dependent numerators. This results in a local minimum in the transition amplitude and a very sharp $\pi$ phase jump at 0.47 (0.37) eV for the s (d) channel. These $\ell$-specific cancelations are highly sensitive to the resonance energies and relative strengths. Their precise positions are also affected by the real-valued amplitudes associated with the “IR-first” two-photon paths (see fig. S3). Despite being relatively small, these background contributions become nonnegligible where the effective amplitude associated with all the other paths vanishes. They have no other consequence on our data analysis and interpretation. The cancelations are measured at 0.54 (0.30) eV for the s (d) channel in the experimental decomposition in Fig. 3 (A and B) that otherwise displays a good agreement with the theoretical prediction.

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Partial-wave decomposition of the investigated two-photon transition. Modulus and phase of the partial-wave transition amplitudes $M_{\ell}(E)$ ($\ell$ = 0 and 2) obtained from experiments (A and B, respectively) and from simulations (C and D, respectively). The simulated partial-wave amplitudes for the reference transition are also shown in (C) and (D), with arg $M_0^{\text{ref}}(E)$ taken as phase origin.</figcaption></figure>
The relative imprint of these $\ell$ contributions over the whole energy range is clearly visible in the polar representation of the final quantum state in Fig. 2 (C and F), which again shows a good correspondence between experiments and simulations and slight shifts in energy. In the 1s3p resonance region, the dominant d wave results in a minimum and change of sign at an angle $\theta_0 \approx 60^\circ$ (37, 42) (red to cyan in plot i), close to the node of $Y_{20}(\theta)$. As the energy increases, the d amplitude becomes less dominant, $\theta_0 \to 90^\circ$ and the sign change disappears (ii). The sharp d-wave cancelation results in a (almost) spherical s state (iii), experimentally distorted by a residual $M_4$ contribution. Conversely, the s-wave cancelation leads to a typical d state (vi). Between these energies, the $\pi$ relative phase of the s and d components produces a strong destructive interference at small angles, resulting in a doughnut shape (iv) evolving toward a d shape (v). In the 1s4p resonance region (vii and viii), one recovers the initial shapes (i and ii). A possible source of discrepancies between the measurements and the theoretical results is the absence of electron correlation (43, 44) in our single-active electron model. This may notably affect the highly sensitive positions of the destructive interferences, but it has no practical consequence on our analysis.

Our measurements give direct access to the complete, angularly resolved dynamics of the two-photon transition leading to photoemission. A way of characterizing this dynamics is provided by the transition delay, a quantity specific to multiphoton processes defined as the local spectral derivative of the transition phase [see (41) and references therein], for each emission angle $\theta$

$$ \tau_{\text{tran}}(E, \theta) = \frac{\partial \arg M(E, \theta)}{\partial E} \qquad (4) $$

The results, shown in Fig. 4A and in the Supplementary Materials, reveal that, in the resonance regions, the delay is strongly positive and weakly anisotropic, while in the intermediate region, it is strongly negative and varies a lot with angle. The physical interpretation of this transition delay is straightforward at each resonance: It represents an effective time during which the electron is transiently trapped in the intermediate bound state, before completing the transition (41). Its effective value is bounded by the experimental duration of the IR, which acts as a temporal gate on the process. In the intermediate region, between resonances, where the delay takes both positive and negative values, the interpretation is not so intuitive. The dynamics there are exclusively shaped by destructive quantum interferences between the two-photon paths through the 1s3p and 1s4p intermediate states, respectively.

To get a deeper insight on how the different spectral components interfere in the temporal buildup of the photoelectron wave packet, we perform a Gabor analysis of the experimental $M(E, \theta)$ using a sliding 210-meV Gaussian gate as a balanced spectrotemporal resolution compromise. Moreover, the spectral gate is narrow enough to allow a direct connection between the obtained time-domain data and the atomic ionization rate (45). The results for three illustrative angles are shown in Fig. 4 (B to D), overlaid with $\tau_{\text{tran}}$ now as a function of E. As anticipated, the revealed dynamics appear to be globally “delayed” at resonances (i.e., around 0.06 and 0.66 eV). Apart from this, they are strongly shaped by the presence of the spectral phase jumps distributed over the covered range. They result in destructive interferences in the temporal profiles, appearing as “holes” in the Gabor transform. In particular, the holes in the intermediate region between resonances shift toward higher energies

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Photoemission dynamics retrieved from the experimental amplitudes $M(E, \theta)$. (A) Angular dependence of the photoemission delay $\tau_{\text{tran}}$ determined by linear fit over 40 meV around four representative energies indicated by color ticks in Fig. 2 and corresponding uncertainty given by the fit SD. (B, C and D) Spectrotemporal buildup of the process revealed by the spectrally gated ionization rates, obtained by Gabor analysis of $M(E, \theta)$ at three representative angles: 0, 65 and 85°, respectively. The transition delay $\tau_{\text{tran}}$ as a function of E is overlaid on each Gabor representation (white lines).</figcaption></figure>
when θ increases, as they follow the spectral phase jump (e.g., from 0.33 eV at 0° to 0.59 eV at 65°). The spectrotemporal structuring of the dynamics therefore highly depends on the photoemission angle (see movies in the Supplementary Materials).

At each angle, the group delay τ<sub>tran</sub> follows the apparent ridge line skirting around the holes in the ionization rate maps. This emphasizes the irrelevance of any attempt to intuitively interpret negative group delays as durations for wave packets shaped by purely quantum processes such as destructive interferences, which eventually comes down to questioning the principle of causality itself (41, 46), see also, e.g., (47).

Our study thus establishes a new step toward the complete characterization of the photoemission process by providing the 3D movie of the angularly resolved dynamics of the two-photon transition leading to photoionization.

By accessing simultaneously the spectral and angular phase variations, the method allows detailed studies in both dimensions and, in particular, gives direct access to the exact photoemission delays as the local derivative of the spectral phase. Our straightforward analysis of the measured data to retrieve the complex transition amplitude is fully experimental (including the calibration procedure) and thus requires no theoretical input. This makes it applicable to a broad variety of atoms, (possibly laser-aligned) molecules, and even solid-state systems using hemispherical analyzers. The detailed information about the quantum photoemission processes will provide a stringent test for theories, in particular, aiming at describing correlated ultrafast multielectronic and vibronic dynamics, a general endeavor common to many fields. Last, the advanced characterization of resonance-enhanced multiphoton ionization provided here opens new prospects for the widely used REMPI technique (48) in a broad range of applications and chemical systems.

## MATERIALS AND METHODS

### Experimental setup

The experimental campaign took place on the 1-kHz attosecond beamline of the ATTOLab facility, pumped by a titanium:sapphire laser providing up to 15-mJ, 25-fs, 804-nm pulses.

The setup is shown on fig. S1. Briefly, it consists of a Mach-Zehnder interferometer in which frequency conversion is performed and two photoelectron spectrometers [a magnetic bottle electron spectrometer (MBES) and a VMIS] in a two-focus geometry (49, 50). For these experiments, 3- to 5-mJ IR pulses are sent on a 90:10 beam splitter. The most intense IR beam is focused with a 2-m focal length lens in a 2-cm-long gas cell to induce HHG. The HHG process yields a comb of mutually coherent odd high harmonics in the XUV range that are first reflected on an IR-antireflective-coated silica plate that transmits most of the intense fundamental radiation. After a 200-nm Al foil that filters out the remaining IR radiation, the XUV beam is refocused into the MBES by a 500-mm focal length gold toroidal mirror at grazing incidence angle (11.5°) in a 2f-2f geometry.

The less intense IR beam, also called dressing beam, can be controlled in energy independently from the generating IR beam using an attenuator composed of a half-wave plate (λ/2) and a polarizer. The energy is set so that an intensity of ≈10<sup>11</sup> W/cm<sup>2</sup> is reached in the MBES and VMIS foci. The dressing delay is controlled by a piezoelectric translation stage (Piezo Jena; 5-nm resolution) (50). The dressing IR is then focused into the MBES using a combination of a 1140-mm lens and a 400-mm lens. This double focus enables to

keep the same oddity of foci number in both arms of the interferometer, ensuring thermal drifts in same direction.

The XUV- and IR-dressing beams are recombined thanks to a drilled mirror with a 3-mm-diameter hole. Their foci are positioned in an argon effusive jet placed in the source volume of the MBES where they produce two-photon XUV-IR ionization. The MBES detection in argon [ionization potential *I*<sub>p</sub>(Ar) = 15.76 eV] allows for a detailed characterization of the XUV radiation, in particular, of harmonic orders 15 and 17 and of their relative phase provided by a RABBIT trace. The detection conditions in the MBES were identical in all the scans performed, allowing direct comparison of the results and spectral calibration of the VMIS data from the XUV ionizing radiation (see the Supplementary Materials). This was verified by comparing the spectral phases of the nonresonant sidebands (SB18 and SB20) measured in the MBES and in the VMIS that showed very similar linear behavior for the three scans.

Both XUV- and IR-dressing beams are then refocused by a second toroidal mirror identical to the first one in a helium gas jet in the source volume of the VMIS (51). This VMIS was specifically designed for attosecond spectroscopy with an original electrostatic lens and a single microchannel plate detector to maximize the spatial and angular resolutions (52). The measurements were taken at low repeller electrode voltage (1 kV) to maximize the resolution at low energy. Electrons with energy up to 10 eV, corresponding to SB22, were detected. When HHG was performed in xenon, all photoelectrons were detected since harmonic 21 (H21) is already in the XUV pulse cutoff. This was not the case for HHG in argon, but we verified that high-energy electrons were not polluting the measurements by recording high repeller traces that compared well to the low repeller ones.

The RABBIT scans were registered simultaneously in the MBES and VMIS over 161 steps of 50 nm, sampling 20 periods of the 400-nm sideband oscillation. Each image was accumulated over 5000 shots. Each scan took around 20 min to be registered, making the interferometer passive stability sufficient to avoid significant temporal drift during the scan. The angular distributions were obtained from the raw images by Abel inversion performed thanks to the direct algorithm for velocity-map imaging system (53).

## Harmonic generation conditions and results

To cover the He response over a broad range (0.85 eV), three scans were measured under three different HHG conditions corresponding to three differently blue-shifted XUV radiations. It is well known that the propagation of the laser pulse in an ionizing medium results

<table><caption>**Table 1. Spectroscopic data of actual and model helium atom.** Relevant eigen-energies (in electron volts) of our model atom and reference values provided by the National Institute of Standards and Technology (NIST) database. The states are labeled with the corresponding configuration of the actual He. The first ionization threshold is taken as the origin of the energy scale, *E* = 0.</caption><thead><tr><th rowspan="2">State</th><th colspan="2">Energy (eV)</th></tr><tr><th>Model</th><th>NIST (57)</th></tr></thead><tbody><tr><td>1s<sup>2</sup></td><td>-24.588</td><td>-24.587</td></tr><tr><td>1s3p</td><td>-1.523</td><td>-1.491</td></tr><tr><td>1s4p</td><td>-0.856</td><td>-0.836</td></tr><tr><td>1s5p</td><td>-0.547</td><td>-0.542</td></tr></tbody></table>
in a blue shift of its central frequency due to the time-varying dispersion of the free electron density. This in turn results in a blue shift (and blue broadening) of the generated harmonics (54). The magnitude of blue shift depends on the generating gas, laser intensity, medium length, and pressure. To vary it, we played on various parameters. First, we used either argon, generating less blue-shifted harmonic combs up to H29, or xenon, enabling stronger blue shift thanks to its low ionization potential, generating harmonics up to H21. Second, we varied the intensity and phase matching conditions by playing both on a first attenuator placed before the beam splitter and on an iris placed just in front of the generation lens to control the focal diameter and energy. Third, the pressure in the 2-cm gas cell was varied as needed.

Scan (1) was taken with HHG in Ar under standard pressure and intensity conditions (the iris diameter was 17 mm and the power at the entrance of the generation chamber was 1.7 W, hence a $1.6 \times 10^{14}$ W/cm$^2$ intensity at focus); the effective generation central photon energy $E_{eff}$ was 1.550 eV, which is slightly higher than the expected photon energy at 804 nm, $E_0 = 1.546$ eV; the length of the generation medium favored blue shift even at relatively low-generation intensity. H15 had a 580-meV full width at half maximum and a central photon energy of 25.25 eV, which virtually brings the atom 150 meV above its 1s3p state. Scan (2) was taken with slightly more blue-shifted conditions, $E_{eff} = 1.560$ eV by generating in high-pressure xenon at a relatively high intensity (22-mm iris diameter, 3-W power, $4.9 \times 10^{14}$ W/cm$^2$ intensity at focus); H15 width was 640 meV. Last, scan (3) was the more blue-shifted case, $E_{eff} = 1.578$ eV with generation in high-pressure xenon and high intensity (30-mm iris diameter, 2-W power, $6.0 \times 10^{14}$ W/cm$^2$ intensity); H15 width was 730 meV.

## Simulation toolbox

All the theoretical results were obtained using a single active electron model of He. The effective electron-core potential reads, in atomic units (55)

$$V(r) = -\frac{1}{r}(1.23 + 2r) e^{-4r} \quad (5)$$

where r is the electron-nucleus distance. The first few eigen-energies of interest of this model atom, computed numerically, are reported in Table 1 and compared to the experimental ones. We see that the agreement between the two sets is very good (<35 meV).

The simulated photoemission data were obtained by solving the TDSE numerically with the Peaceman-Rachford algorithm (56) using an integration time step $\Delta t = 27 \times 10^{-3}$ atomic units (a.u.). We chose the velocity gauge to represent the interaction between the electron and the electric fields. The temporal envelope for the XUV field is a Gaussian, while a $\sin^2$ shape is used for the IR field. The field characteristics (intensity and duration) were adjusted to mimic the experimental conditions, notably their intensities lying safely within the two-photon perturbative regime. A partial-wave expansion of the wave function over few (eight) angular momenta was sufficient to achieve convergence, with a radial grid extending up to $r_{max} = 4000$ a.u. by $\Delta r = 0.125$ a.u. steps. We used a $\cos^{1/8}$ absorber at the edge of the grid (56) to prevent spurious boundary reflections of faster photoelectrons possibly released by any higher-order process. Electron spectra were computed out of the wave function upon propagation using a sliding spectral filter (window operator) with a $\delta E = 12.5$ meV bandwidth (56).

To mimic the experimental results, three sets of TDSE simulations were run under conditions as close as possible to the experimental ones, in particular, using the same effective IR central energies (1.553, 1.560, and 1.578 eV) and XUV bandwidths. The simulated complex amplitudes were then calibrated spectrally in the same way as the experimental ones, as described in the Supplementary Materials.

## SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at https://science.org/doi/10.1126/sciadv.abl7594

## REFERENCES AND NOTES

1. N. A. Cherepkov, G. Raseev, J. Adachi, Y. Hikosaka, K. Ito, S. Motoki, M. Sano, K. Soejima, A. Yagishita, K-shell photoionization of CO: II. Determination of dipole matrix elements and phase differences. *J. Phys. B* **33**, 4213–4236 (2000).

2. M. Lebech, J. C. Houver, A. Lafosse, D. Dowek, C. Alcaraz, L. Nahon, R. R. Lucchese, Complete description of linear molecule photoionization achieved by vector correlations using the light of a single circular polarization. *J. Chem. Phys.* **118**, 9653–9663 (2003).

3. A. McPherson, G. Gibson, H. Jara, U. Johann, T. S. Luk, I. A. McIntyre, K. Boyer, C. K. Rhodes, Studies of multiphoton production of vacuum-ultraviolet radiation in the rare gases. *J. Opt. Soc. Am. B* **4**, 595–601 (1987).

4. M. Ferray, A. L'Huillier, X. F. Li, L. A. Lompre, G. Mainfray, C. Manus, Multiple-harmonic conversion of 1064 nm radiation in rare gases. *J. Phys. B.* **21**, L31–L35 (1988).

5. M. Hentschel, R. Kienberger, C. Spielmann, G. A. Reider, N. Milosevic, T. Brabec, P. Corkum, U. Heinzmann, M. Drescher, F. Krausz, Attosecond metrology. *Nature* **414**, 509–513 (2001).

6. V. Véniard, R. Taieb, A. Maquet, Phase dependence of (N+1)-color (N > 1) IR-UV photoionization of atoms with higher harmonics. *Phys. Rev. A* **54**, 721–728 (1996).

7. P. M. Paul, E. S. Toma, P. Breger, G. Mullot, F. Augé, P. Balcou, H. G. Muller, P. Agostini, Observation of a train of attosecond pulses from high harmonic generation. *Science* **292**, 1689–1692 (2001).

8. H. G. Muller, Reconstruction of attosecond harmonic beating by interference of two-photon transitions. *Appl. Phys. B* **74**, s17–s21 (2002).

9. Y. Mairesse, A. de Bohan, L. J. Frasinski, H. Merdji, L. C. Dinu, P. Monchicourt, P. Breger, M. Kovačev, R. Taieb, B. Carré, H. G. Muller, P. Agostini, P. Salières, Attosecond synchronization of high-harmonic soft x-rays. *Science* **302**, 1540–1543 (2003).

10. M. Schultze, M. Fieß, N. Karpowicz, J. Gagnon, M. Korbman, M. Hofstetter, S. Neppl, A. L. Cavalieri, Y. Komninos, T. Mercouris, C. A. Nicolaides, R. Pazourek, S. Nagele, J. Feist, J. Burgdörfer, A. M. Azzeer, R. Ernstorfer, R. Kienberger, U. Kleineberg, E. Goulielmakis, F. Krausz, V. S. Yakovlev, Delay in photoemission. *Science* **328**, 1658–1662 (2010).

11. K. Klünder, J. M. Dahlström, M. Gisselbrecht, T. Fordell, M. Swoboda, D. Guénot, P. Johnsson, J. Caillat, J. Mauritsson, A. Maquet, R. Taieb, A. L'Huillier, Probing single-photon ionization on the attosecond time scale. *Phys. Rev. Lett.* **106**, 143002 (2011).

12. S. Haessler, B. Fabre, J. Higuet, J. Caillat, T. Ruchon, P. Breger, B. Carré, E. Constant, A. Maquet, E. Mével, P. Salières, R. Taieb, Y. Mairesse, Phase-resolved attosecond near-threshold photoionization of molecular nitrogen. *Phys. Rev. A* **80**, 011404 (2009).

13. M. Huppert, I. Jordan, D. Baykusheva, A. von Conta, H. J. Wörner, Attosecond delays in molecular photoionization. *Phys. Rev. Lett.* **117**, 093001 (2016).

14. L. Seiffert, Q. Liu, S. Zherebtsov, A. Trabattoni, P. Rupp, M. C. Castrovilli, M. Galli, F. Süßmann, K. Wintersperger, J. Stierle, G. Sansone, L. Poletto, F. Frassetto, I. Halfpap, V. Mondes, C. Graf, E. Rühl, F. Krausz, M. Nisoli, T. Fennel, F. Calegari, M. F. Kling, Attosecond chronoscopy of electron scattering in dielectric nanoparticles. *Nat. Phys.* **13**, 766–770 (2017).

15. I. Jordan, M. Huppert, D. Rattenbacher, M. Peper, D. Jelovina, C. Perry, A. von Conta, A. Schild, H. J. Wörner, Attosecond spectroscopy of liquid water. *Science* **369**, 974–979 (2020).

16. A. L. Cavalieri, N. Müller, T. Uphues, V. S. Yakovlev, A. Baltuška, B. Horvath, B. Schmidt, L. Blümel, R. Holzwarth, S. Hendel, M. Drescher, U. Kleineberg, P. M. Echenique, R. Kienberger, F. Krausz, U. Heinzmann, Attosecond spectroscopy in condensed matter. *Nature* **449**, 1029–1032 (2007).

17. M. Swoboda, T. Fordell, K. Klünder, J. M. Dahlström, M. Miranda, C. Buth, K. J. Schafer, J. Mauritsson, A. L'Huillier, M. Gisselbrecht, Phase measurement of resonant two-photon ionization in helium. *Phys. Rev. Lett.* **104**, 103003 (2010).

18. S. B. Schoun, R. Chirla, J. Wheeler, C. Roedig, P. Agostini, L. F. DiMauro, K. J. Schafer, M. B. Gaarde, Attosecond pulse shaping around a cooper minimum. *Phys. Rev. Lett.* **112**, 153001 (2014).

19. M. Kotur, D. Guénot, A. Jiménez-Galán, D. Kroon, E. W. Larsen, M. Louisy, S. Bengtsson, M. Miranda, J. Mauritsson, C. L. Arnold, S. E. Canton, M. Gisselbrecht, T. Carette, J. M. Dahlström, E. Lindroth, A. Maquet, L. Argenti, F. Martín, A. L'Huillier, Spectral phase
measurement of a fano resonance using tunable attosecond pulses. *Nat. Commun.* **7**, 10566 (2016).

20. L. Barreau, C. L. M. Petersson, M. Klinker, A. Camper, C. Marante, T. Gorman, D. Kiesewetter, L. Argenti, P. Agostini, J. González-Vázquez, P. Salières, L. F. DiMauro, F. Martín, Disentangling spectral phases of interfering autoionizing states from attosecond interferometric measurements. *Phys. Rev. Lett.* **122**, 253203 (2019).

21. L. Drescher, T. Witting, O. Kornilov, M. J. J. Vrakking, Phase dependence of resonant and antiresonant two-photon excitations. *Phys. Rev. A* **105**, L011101 (2022).

22. V. Gruson, L. Barreau, Á. Jiménez-Galán, F. Risoud, J. Caillat, A. Maquet, B. Carré, F. Lepetit, J.-F. Hergott, T. Ruchon, L. Argenti, R. Taieb, F. Martin, P. Salières, Attosecond dynamics through a fano resonance: Monitoring the birth of a photoelectron. *Science* **354**, 734–738 (2016).

23. M. Isinger, R. J. Squibb, D. Busto, S. Zhong, A. Harth, D. Kroon, S. Nandi, C. L. Arnold, M. Miranda, J. M. Dahlström, E. Lindroth, R. Feifel, M. Gisselbrecht, A. L'Huillier, Photoionization in the time and frequency domain. *Science* **358**, 893–896 (2017).

24. D. Busto, L. Barreau, M. Isinger, M. Turconi, C. Alexandridi, A. Harth, S. Zhong, R. J. Squibb, D. Kroon, S. Plogmaker, M. Miranda, Á. Jiménez-Galán, L. Argenti, C. L. Arnold, R. Feifel, F. Martin, M. Gisselbrecht, A. L'Huillier, P. Salières, Time–frequency representation of autoionization dynamics in helium. *J. Phys. B* **51**, 044002 (2018).

25. M. Turconi, L. Barreau, D. Busto, M. Isinger, C. Alexandridi, A. Harth, R. J. Squibb, D. Kroon, C. L. Arnold, R. Feifel, M. Gisselbrecht, L. Argenti, F. Martin, A. L'Huillier, P. Salières, Spin–orbit-resolved spectral phase measurements around a fano resonance. *J. Phys. B* **53**, 184003 (2020).

26. T. Remetter, P. Johnsson, J. Mauritsson, K. Varju, Y. Ni, F. Lepine, E. Mansten, M. Kling, J. Khan, R. Lopez-Martens, K. J. Schafer, M. J. J. Vrakking, A. L'Huillier, Attosecond electron wave packet interferometry. *Nat. Phys.* **2**, 323–326 (2006).

27. J. Mauritsson, T. Remetter, M. Swoboda, K. Klünder, A. L'Huillier, K. J. Schafer, O. Ghafur, F. Kelkensberg, W. Siu, P. Johnsson, M. J. J. Vrakking, I. Znakovskaya, T. Uphues, S. Zherebtsov, M. F. Kling, F. Lépine, E. Benedetti, F. Ferrari, G. Sansone, M. Nisoli, Attosecond electron spectroscopy using a novel interferometric pump-probe technique. *Phys. Rev. Lett.* **105**, 053001 (2010).

28. S. A. Aseyev, Y. Ni, L. J. Frasinski, H. G. Muller, M. J. J. Vrakking, Attosecond angle-resolved photoelectron spectroscopy. *Phys. Rev. Lett.* **91**, 223902 (2003).

29. S. Heuser, Á. J. Galán, C. Cirelli, C. Marante, M. Sabbar, R. Boge, M. Lucchini, L. Gallmann, I. Ivanov, A. S. Kheifets, J. Marcus Dahlström, E. Lindroth, L. Argenti, F. Martín, U. Keller, Angular dependence of photoemission time delay in helium. *Phys. Rev. A* **94**, 063409 (2016).

30. C. Cirelli, C. Marante, S. Heuser, C. L. M. Petersson, Á. Jiménez-Galán, L. Argenti, S. Zhong, D. Busto, M. Isinger, S. Nandi, S. Maclot, L. Rading, P. Johnsson, M. Gisselbrecht, M. Lucchini, L. Gallmann, J. M. Dahlström, E. Lindroth, A. L'Huillier, F. Martin, U. Keller, Anisotropic photoemission time delays close to a fano resonance. *Nat. Commun.* **9**, 955 (2018).

31. J. Vos, L. Cattaneo, S. Patchkovskii, T. Zimmermann, C. Cirelli, M. Lucchini, A. Kheifets, A. S. Landsman, U. Keller, Orientation-dependent stereo Wigner time delay and electron localization in a small molecule. *Science* **360**, 1326–1330 (2018).

32. J. Joseph, F. Holzmeier, D. Bresteau, C. Spezzani, T. Ruchon, J. F. Hergott, O. Tcherbakoff, P. D'Oliveira, J. C. Houver, D. Dowek, Angle-resolved studies of XUV-IR two-photon ionization in the RABBITT scheme. *J. Phys. B.* **53**, 184007 (2020).

33. J. Fuchs, N. Douguet, S. Donsa, F. Martin, J. Burgdörfer, L. Argenti, L. Cattaneo, U. Keller, Time delays from one-photon transitions in the continuum. *Optica* **7**, 154–161 (2020).

34. S. Heck, D. Baykusheva, M. Han, J.-B. Ji, C. Perry, X. Gong, H. J. Wörner, Attosecond interferometry of shape resonances in the recoil frame of CF₄. *Sci. Adv.* **7**, eabj8121 (2021).

35. S. Beaulieu, A. Comby, A. Clergerie, J. Caillat, D. Descamps, N. Dudovich, B. Fabre, R. Géneaux, F. Légaré, S. Petit, B. Pons, G. Porat, T. Ruchon, R. Taieb, V. Blanchet, Y. Mairesse, Attosecond-resolved photoionization of chiral molecules. *Science* **358**, 1288–1294 (2017).

36. D. M. Villeneuve, P. Hockett, M. J. J. Vrakking, H. Niikura, Coherent imaging of an attosecond electron wave packet. *Science* **356**, 1150–1153 (2017).

37. D. Busto, J. Vinbladh, S. Zhong, M. Isinger, S. Nandi, S. Maclot, P. Johnsson, M. Gisselbrecht, A. L'Huillier, E. Lindroth, J. M. Dahlström, Fano's propensity rule in angle-resolved attosecond pump-probe photoionization. *Phys. Rev. Lett.* **123**, 133201 (2019).

38. V. Loriot, A. Marciniak, S. Nandi, G. Karras, M. Hervé, E. Constant, E. Plésiat, A. Palacios, F. Martín, F. Lépine, High harmonic generation-2ω attosecond stereo-photoionization interferometry in N₂. *J. Phys. Photonics* **2**, 024003 (2020).

39. C. Cohen-Tannoudji, J. Dupont-Roc, G. Grynberg, *Atom-Photon Interactions: Basic Processes and Applications* (Atomic and Molecular Physics, Wiley, 2008).

40. Á. Jiménez-Galán, F. Martín, L. Argenti, Two-photon finite-pulse model for resonant transitions in attosecond experiments. *Phys. Rev.* **A** **93**, 023429 (2016).

41. M. Vacher, R. Gaillac, A. Maquet, R. Taieb, J. Caillat, Transition dynamics in two-photon ionisation. *J. Opt.* **19**, 114011 (2017).

42. A. S. Kheifets, A. W. Bray, Rabbitt phase transition across the ionization threshold. *Phys. Rev. A* **103**, L011101 (2021).

43. K. L. Ishikawa, K. Ueda, Competition of resonant and nonresonant paths in resonance-enhanced two-photon single ionization of He by an ultrashort extreme-ultraviolet pulse. *Phys. Rev. Lett.* **108**, 033003 (2012).

44. D. I. R. Boll, O. A. Fojón, C. W. McCurdy, A. Palacios, Angularly resolved two-photon above-threshold ionization of helium. *Phys. Rev. A* **99**, 023416 (2019).

45. A. Desrier, A. Maquet, R. Taieb, J. Caillat, Ionization dynamics through a Fano resonance: Time-domain interpretation of spectral amplitudes. *Phys. Rev. A* **98**, 053406 (2018).

46. L. Argenti, Á. Jiménez-Galán, J. Caillat, R. Taieb, A. Maquet, F. Martín, Control of photoemission delay in resonant two-photon transitions. *Phys. Rev. A* **95**, 043426 (2017).

47. M. W. Mitchell, R. Y. Chiao, Causality and negative group delays in a simple bandpass amplifier. *Am. J. Phys.* **66**, 14–19 (1998).

48. C. E. H. Dessent, K. Müller-Dethlefs, Hydrogen-bonding and van der waals complexes studied by KEZE and REMPI spectroscopy. *Chem. Rev.* **100**, 3999–4022 (2000).

49. R. Locher, M. Lucchini, J. Herrmann, M. Sabbar, M. Weger, A. Ludwig, L. Castiglioni, M. Greif, M. Hengsberger, L. Gallmann, U. Keller, Versatile attosecond beamline in a two-foci configuration for simultaneous time-resolved measurements. *Rev. Sci. Instrum.* **85**, 013113 (2014).

50. S. J. Weber, B. Manschwetus, M. Billon, M. Böttcher, M. Bougeard, P. Breger, M. Géléoc, V. Gruson, A. Huetz, N. Lin, Y. J. Picard, T. Ruchon, P. Salières, B. Carré, Flexible attosecond beamline for high harmonic spectroscopy and XUV/near-IR pump probe experiments requiring long acquisition times. *Rev. Sci. Instrum.* **86**, 033108 (2015).

51. A. T. J. B. Eppink, D. H. Parker, Velocity map imaging of ions and electrons using electrostatic lenses: Application in photoelectron and photofragment ion imaging of molecular oxygen. *Rev. Sci. Instrum.* **68**, 3477–3484 (1997).

52. D. Platzer, "Spectroscopie de photoionisation d'atomes et molécules en phase gazeuse aux échelles de temps femtoseconde et attoseconde," thesis, Université Paris-Saclay (2020).

53. G. R. Harrison, J. C. Vaughan, B. Hidle, G. M. Laurent, Davis: A direct algorithm for velocity-map imaging system. *J. Chem. Phys.* **148**, 194010 (2017).

54. C.-G. Wahlström, J. Larsson, A. Persson, T. Starczewski, S. Svanberg, P. Salières, P. Balcou, A. L'Huillier, High-order harmonic generation in rare gases with an intense short-pulse laser. *Phys. Rev. A* **48**, 4709–4720 (1993).

55. R. M. Potvliege, E. Meşe, S. Vučić, Multiphoton ionization and multiphoton resonances in the tunneling regime. *Phys. Rev. A* **81**, 053402 (2010).

56. K. C. Kulander, K. J. Schafer, J. L. Krause, "Atoms in Strong Fields" in *Adv. Atom. Molec. Opt. Phys. Suppl. 1* (Acad. Press, San Diego, 1992), pp. 247–300.

57. A. Kramida, "Atomic energy levels and spectra bibliographic database" (Technical Report, National Institute of Standards and Technology, 2010).

Acknowledgments

**Funding:** This research was supported by Agence Nationale de la Recherche, grant nos. ANR-20-CE30-0007-02-DECAP, ANR-15-CE30-0001-CIMBAAD, ANR-11-EQPX0005-ATTOLAB, and ANR-10-LABX-0039-PALM; COST, grant no. CA18222-AttoChem; and Laserlab-Europe, grant no. EU-H2020-871124. **Author contributions:** Conceptualization: P.S., L.P., R.T., J.C., and F.L. Investigation: A.A., D.P., M.L., G.G., L.M., A.C., L.B., M.D., D.B., and H.J.B.M. Resources: D.B., J.-F.H., and O.T. Software: D.P., D.B., L.P., and R.T. Formal analysis: A.A., D.P., M.L., and A.C. Funding acquisition: P.S., L.P., F.L., R.T., and J.C. Methodology: J.C., R.T., P.S., A.C., and A.A. Supervision: P.S., L.P., and R.T. Validation: V.L. and H.J.B.M. Writing—original draft: A.A., P.S., J.C., and R.T. Writing—review and editing: A.A., H.J.B.M., V.L., F.L., L.P., R.T., J.C., and P.S. **Competing interests:** The authors declare that they have no competing interests. **Data and materials availability:** All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials.

Submitted 3 August 2021

Accepted 31 January 2022

Published 23 March 2022

10.1126/sciadv.abl7594
# Anisotropic dynamics of two-photon ionization: An attosecond movie of photoemission

Alice AutuoriDominique PlatzerMariusz LejmanGuillaume GallicianLucie MaëderAntoine CovoloLea BosseMalay DaluiDavid BresteauJean-François HergottOlivier TcherbakoffHugo J. B. MarrouxVincent LoriotFranck LépineLionel PoissonRichard TaïebJérémie CaillatPascal Salières

Sci. Adv., 8 (12), eabl7594. • DOI: 10.1126/sciadv.abl7594

### View the article online

https://www.science.org/doi/10.1126/sciadv.abl7594

### Permissions

https://www.science.org/help/reprints-and-permissions