Attosecond field emission

News

HomeHome / News / Attosecond field emission

Jun 27, 2023

Attosecond field emission

Nature volume 613, pages 662–666 (2023)Cite this article 14k Accesses 8 Citations 195 Altmetric Metrics details Field emission of electrons underlies great advances in science and technology, ranging

Nature volume 613, pages 662–666 (2023)Cite this article

14k Accesses

8 Citations

195 Altmetric

Metrics details

Field emission of electrons underlies great advances in science and technology, ranging from signal processing at ever higher frequencies1 to imaging of the atomic-scale structure of matter2 with picometre resolution. The advancing of electron microscopy techniques to enable the complete visualization of matter on the native spatial (picometre) and temporal (attosecond) scales of electron dynamics calls for techniques that can confine and examine the field emission on sub-femtosecond time intervals. Intense laser pulses have paved the way to this end3,4 by demonstrating femtosecond confinement5,6 and sub-optical cycle control7,8 of the optical field emission9 from nanostructured metals. Yet the measurement of attosecond electron pulses has remained elusive. We used intense, sub-cycle light transients to induce optical field emission of electron pulses from tungsten nanotips and a weak replica of the same transient to directly investigate the emission dynamics in real time. Access to the temporal properties of the electron pulses rescattering off the tip surface, including the duration τ = (53 as ± 5 as) and chirp, and the direct exploration of nanoscale near fields open new prospects for research and applications at the interface of attosecond physics and nano-optics.

The interaction of atoms and molecules with intense laser fields gives rise to attosecond electron pulses10 that can study the structure and dynamics of these systems on recollision with their parent ion11. Attosecond techniques12 can now gain access to the temporal profile of the recolliding electron pulses and concomitant structural dynamics13,14 in their parent ions by measuring the transient properties of high harmonics15 emitted during the interaction. Studies of the interaction of intense laser fields with nanostructured metals over the past two decades have suggested that the semiclassical concepts3,4,5,6,7,8,9,16,17 developed earlier to describe electron dynamics in atoms can afford a central role in the understanding of the optical field electron emission. In analogy to atoms, electrons set free from the apex of a nanotip at the field crest of an intense laser pulse should also form ultrashort electron pulses (Fig. 1a, inset), which—on recollision with the tip surface about three-quarters of the laser period (T ≈ 2 fs) later—could examine both dynamics and structure. Owing to the ultrashort time interval between emission and recollision events, and in contrast to other emerging electron-pulse technologies18,19,20,21, the electron-pulse wave packet shall undergo a negligible temporal spread, allowing its confinement to sub-cycle timescales.

a, Simplified schematic of the experimental setup. A sub-cycle pulse (orange curve) is spatially separated and focused by a dual concave Ni mirror module. A time delay between the pulses reflected by the inner and outer mirrors is introduced by a piezo stage. Tungsten nanotips (apex radius around 35 nm) or a gas jet of neon atoms can be positioned in the laser focus. Emitted electron spectra are recorded by a time-of-flight spectrometer (acceptance angle about 6°) placed roughly 3 mm downstream of the electron source and aligned along the laser polarization axis. The inset shows that electrons marked by the green shaded curve are (i) set free and accelerated by the intense laser field to form an electron pulse, which, on recollision with the nanotip surface, (ii) can investigate both dynamics as well as structure. On backscattering off the tip surface, (iii) the electron pulse is further accelerated by the laser to escape the interaction area. b, Total electron yield per pulse as a function of increasing peak intensity of the driving laser pulses (yellow points) and its linear fitting (purple line) on a logarithmic scale. c, Electron spectra from the tungsten nanotip (false colour in logarithmic scale) versus peak intensity. Stars and dots denote the cutoff energies. Black and grey dashed lines show the linear fitting of the cutoff energy versus incident peak intensity of the laser pulse. d, Optical emission electron spectra from the tungsten nanotip (red curve) and neon atoms (blue curve) for nearly identical peak intensity (about 40 TW cm−2). e, Same as in c but for neon atoms.

Source data

However the real-time tracking of electron pulses generated in the optical field emission has remained challenging. Although ordinary attosecond streaking techniques can be used to map the temporal structure of the extreme ultraviolet (EUV) electron emission from metal surfaces22 and nanotips23, they cannot directly investigate electron pulses emerging in the optical field emission. Furthermore, the absence of high harmonic emission from laser-driven nanotips constrains the applicability of in situ attosecond techniques24,25,26,27 for examining the structure of the electron pulses in these systems.

Harnessing electron pulses emerging in the optical field emission for realizing new spectroscopies that combine attosecond temporal and nanometre spatial resolution calls for essential advancements on both their generation and their measurement methodologies. On the generation side, the driving laser pulses should be both short and intense so as to confine the tunnelling of electrons into a sub-femtosecond window (<1 fs), as well as to impart the recolliding electron pulses with a de Broglie wavelength (<2.74 Å, >20 eV) that allows atomic-scale investigation of the parent surface28. On the measurement side, in situ attosecond metrologies24,25,26,27 shall be extended to incorporate temporal gating of the optical field emission without relying on the concomitant high harmonic radiation. Measurements of this kind have so far permitted access to the driving-field waveform of light waves by tracking the spectrally integrated currents induced in the bulk of solids29,30 or the cutoff energy variation of rescattered electrons in atoms31 and nanotips32, but a direct time-resolved measurement of attosecond electron pulses in the optical field emission has remained beyond reach.

Guided by the above provisions, we studied the field emission in tungsten nanotips (work function ϕ ≈ 4.5–5.0 eV) using intense (around 1013 W cm−2), sub-cycle (about 1.9 fs) optical transients (centroid energy around 1.8 eV) generated in a light-field synthesizer33,34. The experiments were performed in a multifunctional experimental setup (see Methods and Extended Data Fig. 1) that combines photoemission spectroscopy of atoms and solids, optical pump–probe measurement methodologies and EUV attosecond streaking35 for the sampling of the driving-field waveforms.

We began our experiments by analysing the nonlinearity of the electron yield to the driving-field intensity. A plate detector (see Extended Data Fig. 2) recorded the total number of emitted electrons from the electrically grounded nanotip (yellow points, Fig. 1b) for a range of peak intensities of the impinging laser transient. For the highest peak intensity of around 42 TW cm−2, we recorded up to about 1,000 electrons per pulse. Evaluation of the slope of the total electron yield versus peak intensity (purple line in Fig. 1b) showed an emission nonlinearity of about 1.18 ± 0.09, which is well below the multiphoton threshold ϕ/ħωL ≈ 3 and supports the notion that field-driven electron tunnelling dominates the ionization of tungsten over the entire range of the studied intensities.

Next we conducted a spectral-domain study of the emission under precisely characterized driving fields and interrogated compatibility of our findings with the predictions of semiclassical, single-electron models. Given the multielectron nature of the emission (Fig. 1b), this step is essential for applying semiclassical single-electron concepts later in this work to scrutinize the temporal structure of the electron emission. As space–charge interactions can substantially distort the emission spectra, in comparison with semiclassical predictions, a detailed interrogation of the cutoff energies of the emission can shed light onto potential multielectron contributions.

Electron spectra recorded as a function of the peak intensity of the driving pulse (Fig. 1c) showed the formation of two well-discernible cutoffs (black points, grey points and false-colour plot in Fig. 1c, respectively), whose energies scale linearly with peak intensity. Evaluation of the corresponding slopes for high and low cutoff energies (Ec) by linear fitting (black and grey dashed lines in Fig. 1c) of the experimental data yielded \({s}_{\text{W},\text{high}}^{(\exp )}={\rm{d}}{E}_{{\rm{c}}}/{\rm{d}}{U}_{{\rm{p}}}=118\pm 5.1\) and \({s}_{\text{W},\text{low}}^{(\exp )}={\rm{d}}{E}_{{\rm{c}}}/{\rm{d}}{U}_{{\rm{p}}}=24.1\pm 1.35\), respectively. Here Up is the ponderomotive energy. By taking the ratio between high and low cutoff energy slopes (about 4.91 ± 0.35), we find that it is compatible with that anticipated for backscattered and direct electrons (10Up/2Up ≈ 5) in the single-electron approximation36.

An experimental assessment of the near-field enhancement in the vicinity of the nanotip and its comparison with the theoretical predictions could further attest to the validity of single-electron semiclassical pictures to describe emission under our experimental conditions. To this end, we compared the emission spectra of electrons from the nanotip with those in a low-density neon gas. Figure 1d contrasts electron spectra recorded from W nanotip (red curve) and neon atoms (blue curve) under identical field waveforms. Figure 1e shows electron spectra emanating from Ne over a broad range of driving-pulse intensities. A linear fitting (Fig. 1e, dashed line) of the cutoff energy (Fig. 1e, points) of the neon spectra versus intensity yielded a slope \({s}_{{\rm{Ne}}}^{(\exp )}=9.88\pm 0.34\) that agrees well with the semiclassical predictions under our driving waveforms (sNe = 10.8) (see Methods). The field-enhancement factor in the vicinity of the tungsten nanotip, which is now evaluated as \(f=\sqrt{\frac{{s}_{{\rm{W}},{\rm{high}}}^{(\exp )}}{{s}_{{\rm{Ne}}}^{(\exp )}}}=3.46\pm 0.10\), shows a fair agreement with the theoretical prediction, fth ≈ 3.8 (Methods). This comparison further attests to the compatibility of the emission processes from the nanotip with single-electron, semiclassical concepts and suggests that multielectron charge interactions remain negligible.

Detailed semiclassical simulations based on the experimentally derived quantities (see Methods) further support the above perspective. The simulations accurately reproduce the experimental spectra (Extended Data Fig. 4) and the association of high and low cutoff energies of the emission to backscattered and direct electrons, respectively. The inclusion of multielectron interactions (Extended Data Fig. 6) showed the presence of weak-only effects of space–charge distortions on the emitted spectra and the spectral phase of the electron emission further corroborating the appropriateness of single-electron models to describe the emission process.

The absence of discernible cutoffs at intermediate energies (60–120 eV) in both experiments (Fig. 1c) and simulations (Extended Data Figs. 4 and 6) supports the notion that, under sub-cycle driving, the high-energy part of the emitted spectrum is associated with the recollision of a single electron pulse at the tip surface. A weak, low-energy backscattered emission indicated by the simulations (Extended Data Fig. 4a) is not directly resolved in the experiments (Fig. 1c). Yet the presence of such emission channels will become apparent later in this work by time resolving the optical field emission.

As vacuum is dispersive to electron pulses, a temporal characterization has a concrete meaning at a specific point in space. Because the generated electron pulses scrutinize the ‘sample’ during the recollision with the parent surface as to their temporal structure at this point, it is mostly relevant for harnessing the power of these pulses in time-resolved applications.

To understand how, we revisit the process of strong field recollision of an electron wave packet under an intense optical waveform (Fig. 2). Set free around the peak of a laser field crest, an electron pulse will recollide with the tip surface at an instance tr (Fig. 2a) with an energy of about 3Up (refs. 11,37). An attosecond measurement of the electron pulse entails access into its waveform ψr(t) at the surface of the nanotip or equivalently, into its associated complex spectral amplitude \({\widetilde{\psi }}_{{\rm{r}}}(p)\), in which p is the recollision momentum of the electron. However, as \({\widetilde{\psi }}_{{\rm{r}}}(p)\) is not directly accessible in measurements, it is important to link it to other measurable quantities. Following backscattering off the tip surface, the wave packet acquires further phase from both its interaction with the driving field (Volkov phase) as well owing to its free-space propagation. If we define an auxiliary terminal wave packet \({\widetilde{\psi }}_{{\rm{t}}}(p)\), the spectral intensity \(I(p)={| {\widetilde{\psi }}_{{\rm{t}}}(p)| }^{2}\) (refs. 38,39,40) (as marked in Fig. 2a) can be expressed as (in atomic units):

in which \(S(p,{t}_{{\rm{r}}}\,;{A}_{{\rm{p}}}(t))={\int }_{{t}_{{\rm{r}}}}^{{\rm{\infty }}}\left[\frac{1}{2}\right.{[p+{A}_{{\rm{p}}}(t)]}^{2}-\left.\frac{1}{2}{p}^{2}\right]\,{\rm{d}}t\) denotes the Volkov phase imparted to the electron wave packet by the vector potential Ap(t) of the intense driving pulse (hereafter referred to as pump, see red curve in Fig. 2a) after rescattering at time tr. Hence reconstruction of the recolliding wave packet ψr(t) should be possible if, other than I(p), which is a directly measurable quantity (that is, the spectrum of the electron emission), the phase of \({\widetilde{\psi }}_{{\rm{t}}}(p)\) as well as the waveform Ap(t) are accessed.

a, An electron wave packet ψr(t), sampled on recollision at the nanotip surface, is represented by its momentum distribution \({\widetilde{\psi }}_{{\rm{r}}}(p)\)and corresponding spectral intensity \({\left|{\widetilde{\psi }}_{{\rm{r}}}(p)\right|}^{2}\). Following backscattering, further acceleration of the wave packet by the vector potential of the pump field Ap(t) (red curve) gives rise to a terminal wave packet \({\widetilde{\psi }}_{{\rm{t}}}(p)\) with corresponding spectral intensity \(I(p)={\left|{\widetilde{\psi }}_{{\rm{t}}}(p)\right|}^{2}\). The coherent superposition (dashed yellow curve) of the pump field with a weak, time-delayed (τ) replica (the gate) results in a momentum shift, −eAHAS(t + τ) of the final momentum distribution \(I(p,\tau )={\left|{\widetilde{\psi }}_{{\rm{t}}}(p,\tau )\right|}^{2}\). b, Simulated delay-dependent homochromatic attosecond streaking (HAS) spectrograms for electron pulses without chirp (left panel), positive chirp of 3.5 × 103 as2 (middle panel) and third-order chirp of 1 × 105 as3 (right panel). The dashed white curves highlight the resulting energy-dependent shifts of the intensity modulation.

Source data

Access to the phase of \({\widetilde{\psi }}_{{\rm{t}}}(p)\) is possible by temporally gating the electron emission with a weak replica of the pump pulse (hereafter referred to as gate with vector potential Ag(t)) when \(\eta \equiv {\left|{A}_{{\rm{g}}}(t)\right|}^{2}/{\left|{A}_{{\rm{p}}}(t)\right|}^{2}\ll 1\) (refs. 27,31). In this limit, the pump pulse is solely responsible for releasing the electron wave packet, whereas the gate pulse primarily alters its phase. This is manifested by the shift and reshaping of the terminal photoelectron spectra at the end of the driving pulse (Fig. 2a). If the delay τ between pump and gate pulses is varied (see Methods), the terminal spectral distribution of the released electron, \(I(p,\tau )={\left|{\widetilde{\psi }}_{{\rm{t}}}(p,\tau )\right|}^{2}\) can be approximated as:

in which ψt(t) is the inverse Fourier transform of \({\widetilde{\psi }}_{{\rm{t}}}(p)\) and AHAS(t + τ) represents an effective vector potential, which is explicitly related to the incident vector potential Ag(t) of the gate pulse as shown in Methods and which accounts for the momentum an electron accumulates from the instance of its birth to the detection. Equation (2) implies: (i) a variation of τ permits the composition of a streaking-like spectrogram whose reconstruction can allow retrieving the phase of \({\widetilde{\psi }}_{{\rm{t}}}(p)\) and (ii) the momentum variation of the electron distribution follows AHAS(t). Whereas implication of (ii) (see Methods) allows sampling of the waveform of a light pulse31, (i) is essential for mapping the dynamics of the field emission. To distinguish from conventional attosecond streaking, we refer to this approach as homochromatic attosecond streaking (HAS), highlighting that the carrier frequency of the pump and gate fields is identical.

Shown in Fig. 2b are simulated HAS spectrograms under conditions pertinent to the experiments presented here. Notably, and in close analogy to ordinary attosecond streaking41,42, different types of chirp of the recolliding electron pulse ψr(t) yield distinct visual signatures in the spectrogram, manifested as shifts of the spectral intensity modulation versus delay and energy, as highlighted by the dashed white curves in Fig. 2b.

In our experiments, we derived pump and gate pulses by the spatiotemporal division of the sub-cycle optical transients using a dual-mirror module as shown in Fig. 1a. Figure 3a plots a HAS spectrogram recorded by our apparatus. A ratio η ≈ 6.3 × 10−3 warrants that the gate pulse is sufficiently weak to serve as a phase gate (see Methods), whereas a remaining weak amplitude modulation (5–10%) of the spectrogram versus delay was useful for evaluating the timing between pump and gate pulses, and therewith to clock the recollision instance with respect to the waveform of the driving pump pulse.

a, Homochromatic attosecond streaking (HAS) spectrogram of the optical field emission from a tungsten nanotip consisting of 60 individual electron spectra recorded as a function of the delay (step size of 200 as) between an intense (Ip ≈ 25.4 TW cm−2) single-cycle pulse and a 6.3 × 10−3 times weaker gate pulse. b, Vector potential Ag(t) evaluated by the photoelectron energy cutoff variation in the HAS spectrogram in panel a as a function of the delay (blue) and from an ordinary streaking spectrogram (red) (see Extended Data Fig. 9) in absolute units. c, Simulated HAS spectrogram using the experimentally recorded field waveforms of pump and gate pulses. The dashed vertical lines in a and c denote different rescattering events. d, Terminal energies of direct (green) and backscattered electron emissions (purple) associated (dashed horizontal lines) with corresponding cutoff energies in the spectrogram (c), as well as release times (dashed vertical lines) near-field pump waveform (red curve). EUV-AS, extreme ultraviolet attosecond streaking.

Source data

The evaluated Ag(t) (blue curve in Fig. 3b) by HAS exhibits an excellent waveform matching that derived by EUV attosecond streaking (red curve in Fig. 3b), as attested by the degree of similarity43 (around 0.95). In this case, the HAS-evaluated vector potential represents the near field of the gate pulse in the vicinity of the tip. Moreover, the ratio between the absolute amplitudes of the vector potentials evaluated by the two methods provides a direct measurement of the field-enhancement factor f = 3.73 ± 0.25, which is in close agreement with the result of the methodology highlighted in Fig. 1c (conducted with a different nanotip) and the theoretical estimations.

How time-domain measurement of electron emission in nanostructured materials benefits the intuitive understanding of the process and its real-time clocking can be best appreciated by a close inspection of the oscillating phases of high and low cutoff energies in Fig. 3a. For instance, an apparent delay of the maxima of the corresponding oscillations (highlighted by red and blue dashed curves, respectively, in Fig. 3a) versus delay indicates a retardation in the electron emission at lower energies by approximately a laser cycle (TL ≈ 2.3 fs). This feature, also well reproduced in our simulation (Fig. 3c), shows that the low-energy emission consists of a mixture of direct electrons emerging within the main half-cycle of the driving field (green dots in Fig. 3d) and low-energy backscattered electrons generated approximately a cycle of the driving field later (purple dots in Fig. 3d). Moreover, the uniform amplitude and energy modulation of the photoelectron spectrum (Fig. 3a) over a broad range of energies (50–150 eV) offers compelling evidence that emission is confined to a single electron pulse generated within a field half-cycle.

For the reconstruction of the electron wave packet, we are primarily interested in the properties of the recolliding electron pulses at energies typically higher than 20 eV, that is, where this pulse could serve high-resolution/atomic-scale, spatial analysis of surfaces28. Considering the ratio of the terminal energy of electrons (about 10Up) and that at the recollision instance (about 3.17Up), we conclude that the relevant information is encoded at the high-energy end (>80 eV) of the spectrogram of Fig. 3a, as isolated in Fig. 4a. Figure 4b shows the numerical reconstruction of data in Fig. 4a based on equations (1) and (2), the retrieved field parameters Ap(t) and Ag(t) (Fig. 3b), the absolute time delay τ (Extended Data Fig. 8) and the numerical algorithm detailed in Methods.

a,b, Experimentally recorded (a) and reconstructed (b) HAS spectrograms from a W nanotip. c, Retrieved backscattered electron-pulse spectrum (magenta fill) and its spectral phase (red curve). d, Intensity profile of the electron pulse (magenta fill) and its temporal phase (red line). e, Near field (red curve) and its timing with respect to the attosecond electron pulse (magenta fill). f, Time–frequency analysis of the attosecond electron pulse (false colour is intensity in arbitrary units) and retrieved release times (black curve). The dashed red curve denotes release times calculated semiclassically for the near-field light waveform in e. Error bars represent standard deviations of the corresponding values indicated.

Source data

The retrieved electron-pulse profiles in spectral and time domains are shown in Fig. 4c,d, respectively. The spectrum extends over the range 20–80 eV (Fig. 4c, magenta-shaded curve) and is temporally confined to about 53 ± 5 as, as measured at the full width at half maximum of its intensity profile (Fig. 4d, magenta-shaded curve). Figure 4e compares the retrieved electron pulse and the near-field profiles, and shows that the electron wave packet revisits the surface of the nanotip (magenta-shaded curve) at times close to the zero transition of the driving field (red curve). This observation is compatible with the semiclassical understanding of strong-field rescattering in atoms11,36.

A closer inspection of the retrieved spectral and temporal phases (red curves in Fig. 4c,d) also shows a negligible temporal spread of the electron pulse compared with its Fourier-limited duration (about 50 as). The compatibility of this finding with the semiclassical model can be best appreciated by comparing the group delay of the electron release (solid black line in Fig. 4f) as evaluated by the time–frequency analysis (Fig. 4f, false colour) of the attosecond electron-pulse waveform in Fig. 4d with the semiclassical release times (dashed red curve in Fig. 4f) calculated using the near-field waveform of the driving pulse (red curve in Fig. 4e). The comparison highlights that, as anticipated by the semiclassical recollision model, electrons of energy close to the cutoff, that is, where short and long trajectories merge, carry a minor temporal chirp. This feature, common in near-cutoff high harmonic emission in atoms24,41, further strengthens the link between strong-field optics in atoms and nano-optics systems.

Direct measurement of attosecond electron pulses in the optical field emission extends the repertoire of ultrafast-science tools. The technique presented here could offer a broad range of capabilities for tracking ultrafast collective or correlated electron dynamics in various materials deposited or attached to tungsten nanotips44,45, as well as in gas-phase systems. The recolliding electron pulses shall also enable attosecond-resolved, laser-induced electron diffraction46 experiments at nanotips, as well as nanodiffraction ex situ experiments in specimens placed at nanometric distance from the tip apex47. Such possibilities could enable new routes for exploring the structure and dynamics of condensed matter in four dimensions.

For attosecond EUV streaking measurements (Extended Data Fig. 1a), the sub-cycle transients are focused onto the neon gas jet to generate EUV pulses by high-harmonic generation. The collinearly propagating EUV and optical pulses are spatially separated by a Zr filter, which also acts as an EUV high-pass spectral filter (>70 eV), enabling the isolation of a single attosecond pulse. The EUV and optical pulses are reflected off a dual-mirror assembly, which consists of a Mo/Si inner mirror (centred at around 85 eV) and a nickel outer mirror, respectively. Inner and outer mirrors can be delayed with nanometric resolution (Extended Data Fig. 1a). EUV and optical pulses are focused onto a second Ne gas jet. Photoelectron spectra recorded as a function of the delay between the inner and outer mirrors allow the composition of attosecond streaking spectrograms, which allow the detailed characterization of the attosecond EUV pulse and, notably for these experiments, the field waveform of the optical pulse. Details about the relevant techniques can be found in refs. 35,41,48.

HAS measurements are performed on the same setup by means of (1) automated removal of the Ne gas used to generate high harmonics and the Zr filter, (2) automated replacement of the inner mirror in the dual-mirror module by a Ni-coated one of the same focal lengths (Extended Data Fig. 1b) and (3) the streaking gas nozzle is replaced by an electrically grounded tungsten nanotip. The above setup modifications are executed in a fraction of a minute and warrant identical conditions for all relevant measurements. In the HAS configuration of the setup in Extended Data Fig. 1b, Ni-coated inner and outer mirrors spatially divide the optical pulse into pump (inner mirror beam) and gate (outer mirror) pulses. A delay between the pulses reflected off the inner and outer mirrors, respectively, is introduced by a precision transitional stage (see inset in Extended Data Fig. 1b). The HAS measurements as well as the intensity-scaling measurements of the electron spectrum and yields were performed with driving-field intensities below a critical intensity, at which the irreversible optical damage of the nanotips occurs. This upper limit of intensity was experimentally identified by observing sudden and irreversible drops of electron cutoff energy in intensity-scaling measurements.

For the measurement of the total electron counts per pulse generated in our setup, a thin electrode (size roughly 5 mm × 5 mm) is introduced about 2 mm above the nanotip (Extended Data Fig. 2). This configuration allows detection of released electrons over a solid angle Ω > π steradians. The induced voltage on the thin plate is measured by a lock-in amplifier at the reference frequency of the repetition rate of the driving laser (about 3 KHz). The electronic current is evaluated by dividing the induced voltage by the system impedance (10 MΩ). The total electron count per pulse is in turn obtained by dividing the current by the repetition rate of the laser and the electron charge.

Tungsten nanotips were exposed to intensities of up to I ≈ 45 TW cm−2 without noticing any damage. This was verified by performing the intensity–cutoff energy study shown in Fig. 1b, of both increasing as well as decreasing intensity, and recorded identical curves. When the critical intensity reached I ≈ 45 TW cm−2, the tip is damaged and the cutoff energy irreversibly drops to much lower values without the possibility for recovery unless a new tip is installed.

To explore short-term and long-term stability of our system, we recorded electron spectra as a function of time under identical conditions and for intensities typically higher than those used in the HAS measurements. Data shown in Extended Data Fig. 3 suggest an excellent stability of cutoff energy and electron yield, implying the structural robustness of the nanotip on the timescale of typical measurements (a few minutes) as well as over several hours.

The time-dependent ionization probability from a tungsten nanotip was approximated by the Fowler–Nordheim formula as5,17,49,50:

in which E(t) is the electric-field waveform of the driving pulses, ϕ is the work function of the metal and m, ħ and e are electron mass, reduced Planck’s constant and electron charge, respectively. We calculated electron trajectories using the classical equations of motion in the single-electron limit as7,8,11,17,36:

Here i is the index of each individual trajectory and f0 is the field-enhancement factor. At the end of the pulse, an electron spectrum is evaluated by a spectral binning of the energies of all trajectories weighted by the ionization rate at the instances of their births.

For the experiments described here, we simulated electron spectra from tungsten (ϕ = 4.53 eV) versus peak field intensity of the driving pulse (Extended Data Fig. 4). The driving field (red curve in Extended Data Fig. 4a) used in our simulations was measured by an EUV attosecond streaking setup35,48 integrated in our apparatus. A field-enhancement factor of f0 = 3.46 used in these simulations was derived experimentally as described in the main text. The quiver length of electrons in our intensity range (the longest length of about 1.7 nm at the highest intensity of about 41 TW cm−2) is considerably shorter than the decay length of the near field (about 30 nm). Hence the released electrons experience a nearly homogeneous near field and therefore quenching effects17 owing to the near-field decay can be neglected.

In agreement with the data of Fig. 1c, the simulated electron spectra exhibit two well-discernible energy cutoffs (red and blue dashed lines in Extended Data Fig. 4b) associated with the backscattered (purple line in Extended Data Fig. 4a) and the direct (green line in Extended Data Fig. 4a) electrons. The slopes of high (\({s}_{{\rm{W}},{\rm{high}}}^{({\rm{th}})}={{\rm{d}}E}_{{\rm{c}}}/{{\rm{d}}U}_{{\rm{p}}}=130.3\)) and low (\({s}_{{\rm{W}},{\rm{low}}}^{({\rm{th}})}=26.1\)) cutoff energies agree well with those in our measurements (Fig. 1c). The theory shows further emission cutoffs at energies lower than that of the direct electrons. Because these are relatively weak, they do not leave any direct signatures in the photoelectron spectra. Yet such contributions become visible in HAS spectrograms, as discussed in Fig. 3.

To theoretically estimate the near-field enhancement in the vicinity of the tungsten nanotip, we numerically solved Maxwell’s equations by means of three-dimensional finite-difference time-domain (FDTD) simulations. The nanotip was modelled as shown in Extended Data Fig. 5a as a sphere with radius of 35 nm that smoothly transitions to a cone with an opening angle (single side) of 15° and considering optical properties for tungsten51. The simulations predict a peak field-enhancement factor of about 3.8 close to the surface at the tip apex. For comparison, the spatial distribution of the enhancement at a respective tungsten nanosphere (that is, excluding the cone) is shown in Extended Data Fig. 5b, with a slightly lower peak enhancement factor of approximately 2.7.

To inspect whether charge interaction substantially affects the electron emission dynamics for the considered parameters, we performed semiclassical trajectory simulations using the Mie–mean field–Monte Carlo (M3C) model52. The latter has been used extensively for the study of strong-field ionization in dielectric nanospheres53,54,55 and has recently also been adopted for the description of metallic nanotips56. In brief, we mimic the apex of the nanotip as one-half of a sphere with corresponding radius. The near field is evaluated as the combined linear polarization field owing to the incident pulse (evaluated by means of the Mie solution of Maxwell’s equations) and an extra nonlinear contribution owing to charge interaction treated as a mean field in electrostatic approximation (evaluated by high-order multipole expansion). The latter includes Coulomb interactions among the emitted electrons as well as their image charges (that is, an extra sphere polarization caused by the free electrons). Photoelectron trajectories are generated by Monte Carlo sampling of ionization events at the sphere surface, at which we evaluate tunnelling probabilities within WKB approximation by integration through the barrier provided by the local near field. Trajectories are propagated in the near field by integration of classical equations of motion and accounting for electron–atom collisions through respective scattering cross-sections for electrons moving within the material. To mimic the slightly higher peak enhancement of the linear response near field at a tungsten tip (≈3.8) as compared with a sphere (≈2.7), see Extended Data Fig. 5, we rescaled the incident laser intensity by a factor of 1.4. The performed M3C simulations predict about three times fewer emitted electrons than the experiments, which we attribute to contributions of slow electrons originating from the shank of the nanotip. This is substantiated by comparison of the total electron yields predicted for the nanotip and the half-sphere, obtained through integration of the local ionization rates over the respective surface areas and the pulse duration. However the charge densities at the pole of the sphere and the tip apex are comparable, enabling to inspect the impact of charge interactions within the simplified simulation model. Extended Data Fig. 6a compares multielectron spectra simulated with (solid curves) and without (dashed curves) accounting for charge interactions among the individual electrons and for four representative settings of the peak intensities (and corresponding electron yields) of the driving field, as indicated in the legend. The presence of charge interaction is primarily manifested by a noticeable decline of the yield of direct electrons (<50 eV), which—in accordance with earlier works52,57—can be attributed to a partial trapping of these electrons in the vicinity of the tip. Quasi-static electric fields generated by the trapping of low-energy electrons in turn affect the dynamics of the recolliding electrons and give rise to an increase of the terminal electron-energy cutoff52,57. For the highest intensity and corresponding electron yield, this shift is approximately 8% in energy.

Although transparent manifestations of such effects cannot be discerned in our experiments, it is useful to understand potential implications of charge interactions on the characterization of attosecond electron pulses using HAS. To this end, we extended our analysis to the time domain for pulses of intensity (about 31 TW cm−2) and electron yield (about 600 electrons per pulse), mimicking the experimental conditions in our HAS measurements.

Extended Data Fig. 6b shows the recollision energy distribution of the electron ensemble versus release times without (top) and with (bottom) inclusion of charge interactions. The instantaneous release energy is evaluated by taking the first momentum of the time-resolved spectra (black dashed and solid curves) and the temporal phases (blue dashed and solid curves) by the temporal integration of the instantaneous energy. Their comparison, shown in Extended Data Fig. 6b, suggests marginal differences and, thus, corresponding negligible effects on the temporal characteristics of the electron pulse at the time of recollision.

When the time-domain analysis is extended to the terminal energies of the recolliding electrons versus release time, the space–charge interactions are manifested by a uniform upshift of terminal energies by about 8% but leave the temporal phase of the terminal electron wave packet unaffected (Extended Data Fig. 6c).

To investigate how this energy shift could potentially affect the retrieval of the electron pulse at the instance of recollision, we applied the simulated phase effects on our experimental data (see the section ‘HAS reconstruction methodology’) and evaluated once again the spectral and temporal properties of the recolliding attosecond electron pulse. Key observations associated with space–charge interactions include a weak, uniform shift of the central energy of the recolliding electron by a few electronvolts (Extended Data Fig. 6d) and a subtle change in the temporal profile of the electron pulse (Extended Data Fig. 6e), resulting in an approximately 4-as elongation of its duration at the full width at half maximum, which is within the error of the experimental reconstruction (about 5 as).

The key objective of HAS is to retrieve the temporal structure of an attosecond electron pulse wave packet ψr(t) at the moment of its recollision on its parent surface. As this wave function is not directly accessible, it is necessary to link it to other quantities that are either directly measured in the experiments (such as the terminal spectral intensity \(I(p)={\left|{\widetilde{\psi }}_{{\rm{t}}}(p)\right|}^{2}\) at a detector) or can be reconstructed from the experimental data.

Considering an electron released from and driven back to a surface by a strong pump field Ep(t), its recolliding wave packet ψr(t) can be linked to its terminal spectral amplitude \({\widetilde{\psi }}_{{\rm{t}}}(p)\) at the end of interaction with the driving pulse. The time-dependent recollision wave packet ψr(t) is expressed by means of its Fourier representation \({\widetilde{\psi }}_{{\rm{r}}}(p)={\int }_{-\infty }^{\infty }{\psi }_{{\rm{r}}}({t}_{{\rm{r}}})\exp [{\rm{i}}\frac{{p}^{2}}{2}{t}_{{\rm{r}}}]{\rm{d}}{t}_{{\rm{r}}}\) and, following the recollision, the spectral amplitude is transformed to the terminal form38,39,40,58 (in atomic units):

Here S is the Volkov phase imparted to the electron wave packet only by the vector potential Ap(t) of the pump pulse after recollision at an instance tr, at which the general form of the Volkov phase accumulated from a time instance t1 to a later instance t2 by an electric field with vector potential A(t) is expressed as59:

Note that equation (6) excludes free-space propagation, that is, it vanishes in the absence of the field, and equation (5) reflects the momentum-dependent wave function at the surface, including phases accumulated only by the pump field.

Earlier semiclassical theories of strong-field emission39,40,58 in atoms have suggested that the recolliding wave packet ψr(t) can be expressed by integration over ionization instance t′ before recollision at time tr and over canonical momenta p′ in terms of the ionization amplitude, dictated by the dipole transition Ep(t′)d(p′ + Ap(t)), the scattering amplitude typically described as g(p′ + Ap(tr)) and the Volkov phase that the electron accumulates from t′ to tr as:

Here d and g are the dipole and scattering matrix element, respectively, as defined in refs. 39,40,58, and eiϕt reflects the extra phases acquired during the time evolution of the bound state before ionization. In our experimental setting of HAS in which the Keldysh parameter of γ ≈ 0.38 suggests the tunnelling regime60, those three processes (ionization, propagation and backscattering) can be independently treated without non-adiabatic corrections on the Volkov phase, S (refs. 61,62).

Equation (5) implies that access to ψr(t) is possible if \({\widetilde{\psi }}_{{\rm{t}}}(p)\) and Ap(t) are known. Therefore our goal is to describe how these quantities can be accessed using a phase-gating process of the optical field emission by a weak replica of the driving pulse (gate pulse). Now we inspect the effects of adding the gate pulse on the dynamics of the electron described in equations (5)–(7). We define the gate pulse by a field Eg(t + τ) and its vector potential Ag(t + τ), in which τ is the delay between the pump and gate pulses as described above. By replacing the pump fields and its vector potentials by the superposition of pump and gate pulses in equations (5)–(7) as Ep(t) → Ep(t) + Eg(t + τ) and Ap(t) → Ap(t) + Ag(t + τ), the terminal spectral amplitude perturbed by the gate can be rewritten as:

in which \({\psi }_{{\rm{r}}}^{({\rm{g}})}({t}_{{\rm{r}}},\tau )\) denotes the recolliding electron wave packet perturbed by the extra gate pulse as marked by the superscript (g) to be distinguished from the gate-free counterpart ψr(t) (see equation (7)). Because the gate-free quantity ψr(t) is of interest in this discussion, the subject in this section is how to express \({\widetilde{\psi }}_{{\rm{t}}}(p,\tau )\) in terms of ψr(t) with phase terms introduced by the gate.

First we investigate the perturbed recolliding electron wave packet \({\psi }_{{\rm{r}}}^{({\rm{g}})}({t}_{{\rm{r}}},\tau )\) in equation (8) and how to link it with the gate-free electron wave packet ψr(t). If the gate field is sufficiently weak, that is, \(\eta ={\left|{A}_{{\rm{g}}}(t)\right|}^{2}/{\left|{A}_{{\rm{p}}}(t)\right|}^{2}\ll 1\), the dipole transition and scattering amplitudes can be considered invariant, that is, \(\left[{E}_{{\rm{p}}}(t)+{E}_{{\rm{g}}}(t+\tau )\right]d(p+{A}_{{\rm{p}}}(t)\,+\)\({A}_{{\rm{g}}}(t+\tau ))\approx {E}_{{\rm{p}}}(t)g(p+{A}_{{\rm{p}}}(t))\) and \(g(p+{A}_{{\rm{p}}}(t)+{A}_{{\rm{g}}}(t+\tau ))\approx g(p+{A}_{{\rm{p}}}(t))\) in the expression of the recolliding electron wave packet (equation (7)). In such case, the gate only modifies the phase imparted on the wave packet between ionization and recollision:

In view of these considerations, the variation of the phase (see equation (6)) can be expressed as:

Here the square term of Ag(t + τ) is ignored, as its contribution is negligible compared with the other terms. Equation (10) implies that the gate field introduces an extra phase of \(\Delta S(p+{A}_{{\rm{p}}}(t),{t}_{2},{t}_{1}\,;{A}_{{\rm{g}}}(t+\tau ))\) to the gate-free case. As a result, equation (9) can be rewritten as:

Note that if the last phase term e−iΔS was missing, the expression would be identical to equation (7) and, thus, the gate-free recollision wave function ψr(tr). Hence it would obviously be convenient to remove the e−iΔS term from the integrals, as this would enable to express the recolliding electron wave packet \({\psi }_{{\rm{r}}}^{({\rm{g}})}({t}_{{\rm{r}}},\tau )\) in equation (8) through the unperturbed counterpart ψr(tr) and an extra phase. To proceed further with this idea, we consider two approximations.

First, following the famous saddle-point approximation, the dominant contribution in the integration over the canonical momenta p′ is provided by the kinetic momentum p′ + Ap(t) that equals the kinetic momentum pr of the recolliding electron at the surface. Hence the kinetic momentum term in ΔS may be approximated as \(\Delta S({p}^{{\prime} }{+A}_{{\rm{p}}}(t),{t}_{{\rm{r}}},{t}^{{\prime} }{\rm{;}}{A}_{{\rm{g}}}(t+\tau ))\approx \Delta S({p}_{{\rm{r}}},{t}_{{\rm{r}}},{t}^{{\prime} }{\rm{;}}{A}_{{\rm{g}}}(t+\tau ))\).

Second, because the exponent term of the extra phase e−iΔS is oscillating slowly compared with e−iS in the time integration over t′ in equation (11), the extra phase ΔS can be approximated by a time average \(\overline{\Delta S}\) within a time window Δt:

As our main goal is to reconstruct attosecond electron wave packets that contribute to the spectral cutoff, we choose Δt as the time interval between ionization and recollision of the classical backscattering trajectory that results in the highest final kinetic energy. To evaluate the averaged phase \(\overline{\Delta S}\) and thereby simplify the analytical form of equation (12), ΔS can be expressed as:

By inserting equation (13) into equation (12), the effective (averaged) phase variation \(\overline{\Delta S}\,\) can now be evaluated:

in which \({\bar{A}}_{{\rm{g}}}^{({\rm{b}})}(t)\) is defined as:

Using the above-described approximations now enables to pull the extra phase term out of the integrations in equation (11) and, considering the sign flip \(\overline{\Delta S}\to -\overline{\Delta S}\) at the backscattering instance, the perturbed recolliding electron wave packet can be expressed through the gate-free wave packet and the extra averaged phase term \({\psi }_{{\rm{r}}}^{({\rm{g}})}({t}_{{\rm{r}}},\tau )\approx {\psi }_{{\rm{r}}}({t}_{{\rm{r}}}){{\rm{e}}}^{{\rm{i}}\overline{\Delta S}}\) in equation (8).

We now move on to discuss how the electron wave packet can be described at the end of the interaction in the presence of the pump and gate fields (equation (8)). Taking the results of equations (11)–(14) and restoring pr with the kinetic momentum p + Ap(t), equation (8) can be rewritten as:

Because the integration ranges for S and ΔS are identical (from tr to ∞), the two phases can be merged into a single equation (S′ = S + ΔS),

in which AHAS(t) is hereafter referred to as the effective HAS vector potential and reads:

This expression of the effective HAS vector potential is compatible with the classical momentum accumulation during the excursion from the ionization to the detection, \(\Delta p=-\,{\rm{e}}\left[2A({t}_{{\rm{r}}})-A({t}_{{\rm{r}}}-\Delta t)\right]\), under the rescattering condition, \({\int }_{{t}_{{\rm{r}}}-\Delta t}^{{t}_{{\rm{r}}}}{A}_{{\rm{g}}}(t){\rm{d}}t=\Delta t{A}_{{\rm{g}}}({t}_{{\rm{r}}}-\Delta t)\)(refs. 11,36,63). Using equation (17), the terminal electron amplitude \({\widetilde{\psi }}_{{\rm{t}}}(p,\tau )\) can be expressed as:

Equation (19) implies that the gate also contributes to the terminal momentum of the electron wave packet by AHAS(tr + τ), which depends on the time delay τ. As described in equation (5), the momentum contribution Ap(t) of the pump field is already incorporated in the gate-free terminal electron spectral amplitude \({\widetilde{\psi }}_{{\rm{t}}}(p)\), whose intensity is directly accessible in experiments. Therefore it is convenient for the analysis of HAS data to express equation (19) with the terminal form of \({\widetilde{\psi }}_{{\rm{t}}}(p).\) We decompose the phase in equation (19) by means of \(S(p,{t}_{2},{t}_{1}\,;{A}_{{\rm{p}}}(t)+{A}_{{\rm{g}}}(t+\tau ))\approx S(p,{t}_{2},{t}_{1}\,;{A}_{{\rm{p}}}(t))+\Delta S(p+{A}_{{\rm{p}}}(t),{t}_{2},{t}_{1}\,;{A}_{{\rm{HAS}}}(t+\tau ))\) and rewrite equation (19) as:

In analogy to equation (11), if the e−iΔS term vanishes, the above equation is identical to equation (5), which links the recolliding electron wave packet ψr(t) to the terminal spectral amplitude \({\widetilde{\psi }}_{{\rm{t}}}(p)\). Here, similar to the Fourier representation of the recollision wave packet, we define the Fourier pair of the terminal electron wave packet as:

Note that the terminal electron wave packet ψt(t) is an auxiliary electron wave packet that contains time–structure information of the recolliding electron wave packet ψr(t) at the recollision surface with the momenta translated by the Volkov propagation with the exponent \(\exp \left[-{\rm{i}}S(p,\infty ,{t}_{{\rm{r}}}{\rm{;}}{A}_{{\rm{p}}}(t))\right]\) (see equations (1) and (5)), but without the phase from space propagation to the detection. Using the terminal electron wave packet ψt(t) (equation (21)), the terminal electron spectral amplitude (equation (20)) can be further simplified as:

under the condition that the variation of the vector potential is weak during the time window of the recollision. The HAS spectrogram equation then reads:

Equation (23) describes a spectrogram whose reconstruction allows access to the final electron wave packet ψt(t) and, correspondingly, \({\widetilde{\psi }}_{{\rm{t}}}(p)\) as well as AHAS(t).

An explicit relationship between the incident gate vector potential Ag(t) and the effective HAS vector potential AHAS(t) can be best understood in the Fourier domain. Using the Fourier expansion, \({A}_{{\rm{g}}}(t)={\int }_{-{\rm{\infty }}}^{{\rm{\infty }}}{\rm{d}}\omega {\widetilde{A}}_{{\rm{g}}}(\omega ){{\rm{e}}}^{{\rm{i}}\omega t}\), the effective HAS vector potential can be expressed as,

in which the newly introduced multiplier \(\widetilde{g}(\omega )\) is defined as:

As shown in equation (24), the Fourier components of the effective HAS vector potential \({\widetilde{A}}_{{\rm{HAS}}}(\omega )\) is related to those of the incident gate vector potential \({\widetilde{A}}_{{\rm{g}}}(\omega )\) by multiplication of \(\widetilde{g}(\omega )\)

The gate multiplier \(\widetilde{g}(\omega )\) is independent from \({\widetilde{A}}_{{\rm{g}}}(\omega )\). This allows the possibility of the complete characterization of Ag(t) from AHAS(t) imprinted in a HAS spectrogram.

To better visualize the concept of AHAS(t) and to verify the validity of the assumptions used in the above derivation, a semiclassical simulation of a HAS spectrogram was performed using single-cycle pulses. The photoelectron spectrum cutoff energy variation evaluated by the HAS spectrogram is compared with the effective HAS vector potential AHAS(t) calculated using equation (24) (Extended Data Fig. 7). The multiplier \(\widetilde{g}(\omega )\) depends on the excursion time Δt between ionization and the backscattering event of the highest-energy electron. On the basis of the well-established recollision model, 0.685 times the central excursion period64,65 (central period of \(E(\omega )/\left.{\omega }^{2}\right)\)), which corresponds to about 0.85TL, was used for Δt to evaluate \(\widetilde{g}(\omega )\). Here TL is the centroid period of the laser pulse. Extended Data Fig. 7b–d shows that the cutoff energy variation in a HAS spectrogram closely follows AHAS(t) (black curve), as calculated by the unmodified vector potential of the incident pulse (dashed red curve), regardless of the carrier-envelope phase.

The above discussion suggests that, by tracing the variation of the cutoff energy in a HAS spectrogram, we can obtain AHAS(t) (red curves in Extended Data Fig. 8a,b). Therefore access to the Fourier components of the effective HAS vector potential allows the characterization of the vector potential of the incident gate: \({\widetilde{A}}_{{\rm{g}}}(\omega )={\widetilde{g}}^{-1}(\omega ){\widetilde{A}}_{{\rm{HAS}}}(\omega )\) (Extended Data Fig. 8c,d). The retrieved incident vector potential Ag(t) is shown in blue in Extended Data Fig. 8b.

Identification of the zero delay between pump and gate pulses in a HAS spectrogram can be obtained with various methods. Here we opted for a method that allows the absolute delay to be derived directly from the HAS spectrogram. Even though the difference between the intensities of pump and gate pulses is more than two orders of magnitude, discernible modulations (approximately 5–10%) of the spectral amplitude of the spectrogram remain. In a HAS spectrogram, the total photoelectron yield variation can be evaluated by spectral integration at each delay point (Extended Data Fig. 8e). The absolute zero-delay point can be found as the delay point at which the yield is maximally varied (vertical dashed line in Extended Data Fig. 8e).

EUV attosecond streaking provides access to the detailed field waveform of a pulse33,34,35. Because this technique of field characterization is integrated in our experimental setup, it allows us to benchmark HAS as a field-characterization method.

Extended Data Fig. 9a,b shows the HAS and EUV attosecond streaking measurements, respectively. The vector potential waveform of the incident gate pulse retrieved from the cutoff analysis in HAS (red curve in Extended Data Fig. 9c) and that from EUV attosecond streaking (blue curve in Extended Data Fig. 9c) show excellent agreement, as verified by the degree of similarity of about 0.95 (ref. 43), and support the notion that the gate pulse indeed acts as a phase gate.

The compact description of HAS as a spectrogram implied by equations (2) and (23) assumes that the weak replica of the pump field acts as a nearly pure phase gate on the electrons released by the pump. In other words, it can modify the momentum of electrons released by the pump field but does not greatly influence the process of electron ionization. Yet, unless the ionization nonlinearities are well understood (for instance, in atoms), a theoretical estimate of the required ratio between pump and gate pulses for attaining a sufficiently pure phase gate requires experimental validation.

To identify safe limits within which the above condition is met, we performed HAS measurements under different gate strengths and compared the vector potential waveforms extracted from HAS to those characterized by EUV attosecond streaking. As shown in Extended Data Fig. 10a,b, the vector potential waveforms from two techniques achieve best agreement at low gate/pump intensity ratio (η < 10−2). At higher intensity ratios, we observe a gradually increasing disagreement between the reconstructed waveforms with the two methods (Extended Data Fig. 10c,d), implying that the gate pulse no longer serves as a weak perturbation. These measurements suggest that, for the studied system, HAS measurements require a gate pulse whose intensity is about 10−2 lower than the pump intensity.

At the first stage of the reconstruction of the HAS spectrogram, the terminal electron wave packet ψt(t) is retrieved, as its spectral intensity \({\left|{\widetilde{\psi }}_{{\rm{t}}}(p)\right|}^{2}\) can be directly obtained by a gate-free photoelectron spectrum. Therefore the reconstruction problem is reduced to retrieval of the spectral phase.

As motivated in the main text, we isolated a spectral area of interest (AOI) from 80 to 230 eV (Extended Data Fig. 11a,b). With this region of interest, the terminal wave packet can be expressed as:

Here φ(ω) is the spectral phase of the electron wave packet modelled as a polynomial series up to the sixth order,

in which Dn and ωc are the nth order dispersion and central frequency, respectively. The reconstruction is based on a least-squares algorithm written in MATLAB, which aims at the total minimization of the difference among the experimental (Fig. 4a and Extended Data Fig. 11a) and reconstructed spectrogram (Fig. 4b and Extended Data Fig. 11b). To further increase the fidelity of the reconstruction, we also simultaneously fit the differential map D(E, τ) of a HAS spectrogram I(E, τ), which is defined as:

The differential map is useful because it can eliminate the unmodulated intensity along the delay axis and allows the retrieval algorithm to reconstruct fine details of the experimental trace (Extended Data Fig. 11c–e). As an initial guess for the phase, zero phase was used. The retrieval of the terminal electron wave packet is shown in Extended Data Fig. 11f,g.

In a next stage of the reconstruction, the recolliding electron pulse, which is the key quantity in this work, is evaluated by the inverse Volkov propagation of the retrieved terminal electron wave packet as:

The Volkov basis is reconstructed by precisely measuring the pump-field waveform and its timing with respect to the emission. The retrieved recolliding electron pulse is shown in Fig. 4.

The data supporting the conclusions in the paper are available from the corresponding authors on reasonable request. Source data are provided with this paper.

Han, J.-W., Moon, D.-I. & Meyyappan, M. Nanoscale vacuum channel transistor. Nano Lett. 17, 2146–2151 (2017).

Article ADS CAS Google Scholar

Brodusch, N., Demers, H. & Gauvin, R. Field Emission Scanning Electron Microscopy: New Perspectives for Materials Characterization (Springer, 2018).

Hommelhoff, P. & Kling, M. Attosecond Nanophysics: From Basic Science to Applications (Wiley, 2015).

Dombi, P. et al. Strong-field nano-optics. Rev. Mod. Phys. 92, 25003 (2020).

Article CAS Google Scholar

Hommelhoff, P., Sortais, Y., Aghajani-Talesh, A. & Kasevich, M. A. Field emission tip as a nanometer source of free electron femtosecond pulses. Phys. Rev. Lett. 96, 77401 (2006).

Article ADS Google Scholar

Ropers, C., Solli, D. R., Schulz, C. P., Lienau, C. & Elsaesser, T. Localized multiphoton emission of femtosecond electron pulses from metal nanotips. Phys. Rev. Lett. 98, 43907 (2007).

Article ADS CAS Google Scholar

Krüger, M., Schenk, M. & Hommelhoff, P. Attosecond control of electrons emitted from a nanoscale metal tip. Nature 475, 78–81 (2011).

Article Google Scholar

Piglosiewicz, B. et al. Carrier-envelope phase effects on the strong-field photoemission of electrons from metallic nanostructures. Nat. Photonics 8, 37–42 (2014).

Article ADS CAS Google Scholar

Bormann, R., Gulde, M., Weismann, A., Yalunin, S. V. & Ropers, C. Tip-enhanced strong-field photoemission. Phys. Rev. Lett. 105, 147601 (2010).

Article ADS CAS Google Scholar

Niikura, H. et al. Sub-laser-cycle electron pulses for probing molecular dynamics. Nature 417, 917–922 (2002).

Article ADS CAS Google Scholar

Corkum, P. B. Plasma perspective on strong field multiphoton ionization. Phys. Rev. Lett. 71, 1994–1997 (1993).

Article ADS CAS Google Scholar

Krausz, F. & Ivanov, M. Attosecond physics. Rev. Mod. Phys. 81, 163–234 (2009).

Article ADS Google Scholar

Baker, S. et al. Probing proton dynamics in molecules on an attosecond time scale. Science 312, 424–427 (2006).

Article ADS CAS Google Scholar

Wolter, B. et al. Ultrafast electron diffraction imaging of bond breaking in di-ionized acetylene. Science 354, 308–312 (2016).

Article ADS CAS Google Scholar

Ferray, M. et al. Multiple-harmonic conversion of 1064 nm radiation in rare gases. J. Phys. B At. Mol. Opt. Phys. 21, L31–L35 (1988).

Article CAS Google Scholar

Zherebtsov, S. et al. Controlled near-field enhanced electron acceleration from dielectric nanospheres with intense few-cycle laser fields. Nat. Phys. 7, 656–662 (2011).

Article CAS Google Scholar

Herink, G., Solli, D. R., Gulde, M. & Ropers, C. Field-driven photoemission from nanostructures quenches the quiver motion. Nature 483, 190–193 (2012).

Article ADS CAS Google Scholar

Baum, P. & Zewail, A. H. Attosecond electron pulses for 4D diffraction and microscopy. Proc. Natl Acad. Sci. 104, 18409–18414 (2007).

Article ADS CAS Google Scholar

Hassan, M. T., Baskin, J. S., Liao, B. & Zewail, A. H. High-temporal-resolution electron microscopy for imaging ultrafast electron dynamics. Nat. Photon. 11, 425–430 (2017).

Article ADS CAS Google Scholar

Morimoto, Y. & Baum, P. Diffraction and microscopy with attosecond electron pulse trains. Nat. Phys. 14, 252–256 (2018).

Article CAS Google Scholar

Priebe, K. E. et al. Attosecond electron pulse trains and quantum state reconstruction in ultrafast transmission electron microscopy. Nat. Photon. 11, 793–797 (2017).

Article ADS CAS Google Scholar

Cavalieri, A. L. et al. Attosecond spectroscopy in condensed matter. Nature 449, 1029–1032 (2007).

Article ADS CAS Google Scholar

Förg, B. et al. Attosecond nanoscale near-field sampling. Nat. Commun. 7, 11717 (2016).

Article ADS Google Scholar

Dudovich, N. et al. Measuring and controlling the birth of attosecond XUV pulses. Nat. Phys. 2, 781–786 (2006).

Article CAS Google Scholar

Shafir, D. et al. Resolving the time when an electron exits a tunnelling barrier. Nature 485, 343–346 (2012).

Article ADS CAS Google Scholar

Kim, K. T. et al. Manipulation of quantum paths for space–time characterization of attosecond pulses. Nat. Phys. 9, 159–163 (2013).

Article CAS Google Scholar

Kim, K. T. et al. Petahertz optical oscilloscope. Nat. Photon. 7, 958–962 (2013).

Article ADS CAS Google Scholar

Van Hove, M. A., Weinberg, W. H. & Chan, C.-M. Low-Energy Electron Diffraction: Experiment, Theory and Surface Structure Determination Vol. 6 (Springer, 2012).

Schiffrin, A. et al. Optical-field-induced current in dielectrics. Nature 493, 70–74 (2013).

Article ADS Google Scholar

Rybka, T. et al. Sub-cycle optical phase control of nanotunnelling in the single-electron regime. Nat. Photon. 10, 667–670 (2016).

Article ADS CAS Google Scholar

Kim, Y. H. et al. Attosecond streaking using a rescattered electron in an intense laser field. Sci. Rep. 10, 22075 (2020).

Article ADS CAS Google Scholar

Hoff, D. et al. Tracing the phase of focused broadband laser pulses. Nat. Phys. 13, 947–951 (2017).

Article CAS Google Scholar

Wirth, A. et al. Synthesized light transients. Science 334, 195 LP–200 (2011).

Article ADS Google Scholar

Hassan, M. T. et al. Optical attosecond pulses and tracking the nonlinear response of bound electrons. Nature 530, 66–70 (2016).

Article ADS CAS Google Scholar

Goulielmakis, E. et al. Direct measurement of light waves. Science 305, 1267–1269 (2004).

Article ADS CAS Google Scholar

Paulus, G. G., Becker, W., Nicklich, W. & Walther, H. Rescattering effects in above-threshold ionization: a classical model. J. Phys. B At. Mol. Opt. Phys. 27, L703–L708 (1994).

Article ADS CAS Google Scholar

L’Huillier, A. et al. High-order harmonic-generation cutoff. Phys. Rev. A 48, R3433–R3436 (1993).

Article ADS Google Scholar

Lewenstein, M., Balcou, P., Ivanov, M. Y., L’Huillier, A. & Corkum, P. B. Theory of high-harmonic generation by low-frequency laser fields. Phys. Rev. A 49, 2117–2132 (1994).

Article ADS CAS Google Scholar

Lewenstein, M., Kulander, K. C., Schafer, K. J. & Bucksbaum, P. H. Rings in above-threshold ionization: a quasiclassical analysis. Phys. Rev. A 51, 1495–1507 (1995).

Article ADS CAS Google Scholar

Milošević, D. B., Paulus, G. G., Bauer, D. & Becker, W. Above-threshold ionization by few-cycle pulses. J. Phys. B At. Mol. Opt. Phys. 39, R203–R262 (2006).

Article ADS Google Scholar

Goulielmakis, E. et al. Single-cycle nonlinear optics. Science 320, 1614–1617 (2008).

Article ADS CAS Google Scholar

Hofstetter, M. et al. Attosecond dispersion control by extreme ultraviolet multilayer mirrors. Opt. Express 19, 1767–1776 (2011).

Article ADS Google Scholar

Hecht, E. Optics (Pearson, 2017).

Tafel, A., Meier, S., Ristein, J. & Hommelhoff, P. Femtosecond laser-induced electron emission from nanodiamond-coated tungsten needle tips. Phys. Rev. Lett. 123, 146802 (2019).

Article ADS CAS Google Scholar

Li, C. et al. Extreme nonlinear strong-field photoemission from carbon nanotubes. Nat. Commun. 10, 4891 (2019).

Article ADS Google Scholar

Meckel, M. et al. Laser-induced electron tunneling and diffraction. Science 320, 1478–1482 (2008).

Article ADS CAS Google Scholar

Garg, M. & Kern, K. Attosecond coherent manipulation of electrons in tunneling microscopy. Science 367, 411–415 (2020).

Article ADS CAS Google Scholar

Schultze, M. et al. State-of-the-art attosecond metrology. J. Electron. Spectros. Relat. Phenom. 184, 68–77 (2011).

Article CAS Google Scholar

Fowler, R. H. & Nordheim, L. Electron emission in intense electric fields. Proc. R. Soc. Lond. A Math. Phys. Eng. Sci. 119, 173–181 (1928).

ADS CAS MATH Google Scholar

Yalunin, S. V., Gulde, M. & Ropers, C. Strong-field photoemission from surfaces: theoretical approaches. Phys. Rev. B 84, 195426 (2011).

Article ADS Google Scholar

Rakić, A. D., Djurišić, A. B., Elazar, J. M. & Majewski, M. L. Optical properties of metallic films for vertical-cavity optoelectronic devices. Appl. Opt. 37, 5271–5283 (1998).

Article ADS Google Scholar

Seiffert, L., Zherebtsov, S., Kling, M. F. & Fennel, T. Strong-field physics with nanospheres. Adv. Phys. X 7, 2010595 (2022).

Google Scholar

Süßmann, F. et al. Field propagation-induced directionality of carrier-envelope phase-controlled photoemission from nanospheres. Nat. Commun. 6, 7944 (2015).

Article ADS Google Scholar

Seiffert, L. et al. Attosecond chronoscopy of electron scattering in dielectric nanoparticles. Nat. Phys. 13, 766–770 (2017).

Article CAS Google Scholar

Liu, Q. et al. Ionization-induced subcycle metallization of nanoparticles in few-cycle pulses. ACS Photonics 7, 3207–3215 (2020).

Article CAS Google Scholar

Schötz, J. et al. Onset of charge interaction in strong-field photoemission from nanometric needle tips. Nanophotonics 10, 3769–3775 (2021).

Article Google Scholar

Seiffert, L. et al. Trapping field assisted backscattering in strong-field photoemission from dielectric nanospheres. J. Mod. Opt. 64, 1096–1103 (2017).

Article ADS Google Scholar

Suárez, N., Chacón, A., Ciappina, M. F., Biegert, J. & Lewenstein, M. Above-threshold ionization and photoelectron spectra in atomic systems driven by strong laser fields. Phys. Rev. A 92, 63421 (2015).

Article ADS Google Scholar

Wolkow, D. M. Über eine Klasse von Lösungen der Diracschen Gleichung. Z. Phys. 94, 250–260 (1935).

Article ADS MATH Google Scholar

Keldysh, L. V. Ionization in the field of a strong electromagnetic wave. Sov. Phys. JETP 20, 1307–1314 (1965).

Google Scholar

Pedatzur, O. et al. Attosecond tunnelling interferometry. Nat. Phys. 11, 815–819 (2015).

Article CAS Google Scholar

Kneller, O. et al. A look under the tunnelling barrier via attosecond-gated interferometry. Nat. Photon. 16, 304–310 (2022).

Article ADS CAS Google Scholar

Scrinzi, A., Ivanov, M. Y., Kienberger, R. & Villeneuve, D. M. Attosecond physics. J. Phys. B At. Mol. Opt. Phys. 39, R1–R37 (2005).

Google Scholar

Paulus, G. G., Becker, W. & Walther, H. Classical rescattering effects in two-color above-threshold ionization. Phys. Rev. A 52, 4043–4053 (1995).

Article ADS CAS Google Scholar

Becker, W., Goreslavski, S. P., Milošević, D. B. & Paulus, G. G. The plateau in above-threshold ionization: the keystone of rescattering physics. J. Phys. B At. Mol. Opt. Phys. 51, 162002 (2018).

Article ADS Google Scholar

Download references

We thank J. Apportin and B. Kruse for support with the FDTD simulations. This work was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) SFB 1477 ‘Light-Matter Interactions at Interfaces’, project number 441234705.

Open access funding provided by Universität Rostock.

Institut für Physik, Universität Rostock, Rostock, Germany

H. Y. Kim, S. Mandal, L. Seiffert, T. Fennel & E. Goulielmakis

Max Planck Institute for Solid State Research, Stuttgart, Germany

M. Garg

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

H.Y.K., M.G. and S.M. performed experiments and analysed the experimental data. H.Y.K., M.G., L.S. and T.F. performed theoretical modelling and simulation. E.G. conceived and supervised the project. All authors contributed to the preparation of the manuscript.

Correspondence to E. Goulielmakis.

The authors declare no competing interests.

Nature thanks Qing Dai, Michael Krüger and Peng Zhang for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

EUV (a) and HAS (b) configurations. The schematics of the adjustable iris and the dual inner mirror assembly are reproduced with permission by SmarAct GmbH.

Electrons emitted from a tungsten nanotip exposed to single-cycle pulses were captured by a thin electrode placed about 2 mm above the tip. The induced voltage on the electrode was recorded by a lock-in amplifier.

a, Consecutive measurements of photoelectron spectra recorded under identical conditions and over a time range of about 100 s (left). The grey line indicates the cutoff energy. Electron yield evaluated by integration of the individual spectra in the left panel (middle). Individual spectra superimposed for comparison (right). b, The same measurement as in a but over a time period of about 80 min.

Source data

a, Terminal electron energy as a function of the ionization time under the electric field of the incident single-cycle pulses (red curve). A peak field intensity of 41 TW cm−2 was used in accordance with the experiments. Purple and green lines indicate terminal energies of backscattered and direct electrons as functions of birth time, for which the ionization rate is encoded on the line shading in logarithmic scale. b, Electron spectra plotted as a function of the peak field intensity. Red and blue dashed lines indicate the high and low cutoff energies, respectively.

Source data

Evaluated spatial distribution of the field-enhancement factor at the apex of a tungsten nanotip with an apex radius of 35 nm and opening angle of 15° (a) and a similarly sized tungsten sphere (b).

Source data

a, Electron spectra obtained from three-dimensional semiclassical M3C trajectory simulations for the experimental parameters without (dashed curves) and with (solid curves) the inclusion of charge interactions for representative pulse peak intensities. b, Recollision energy distribution as a function of the release time without (top) and with (bottom) charge interactions (false colour denotes spectral intensity). Black and blue curves are the instantaneous energy and the temporal phase, respectively, whereas dashed and solid curves mark the cases without and with the charge interactions, respectively. c, The same analysis as in b but for the terminal energy of the electron ensemble. The charge interaction free instantaneous energy (dashed black curve) is virtually identical to that modified by charge interactions (solid black curve), except for a uniform shift (dashed grey curve) along the energy axis. d, Retrieved spectral intensity (blue) and phase (red) of the recolliding electron pulse with (dashed curves) and without (solid curves) inclusion of space–charge effects on the terminal electron pulse. e, The same analysis as in d but in the temporal domain.

Source data

a, Gate multiplier function \(\widetilde{g}(\omega )\) as defined in the text. b–d, Semiclassically simulated HAS spectrograms at different carrier-envelope phase (CEP) settings of the driving pulse. The red dashed and black solid curves indicate waveforms of the incident gate vector potential Ag(t) and effective HAS gate vector potential AHAS(t) calculated using equation (24), respectively.

Source data

a, Measured HAS spectrogram (Fig. 3a). b, AHAS(t) (red) directly extracted by tracing the cutoff energy variation in the spectrogram and the retrieved incident vector potential Ag(t) (blue). c, Amplitude (blue) and phase (red) of the inverse of the multiplier \(\widetilde{g}{\left(\omega \right)}^{-1}\). d, Spectral intensity of AHAS(t) (red) and incident Ag(t) (blue) vector potentials. e, Identification of timing between gate and pump pulses in a HAS spectrogram. The variation of electron yields is evaluated by the spectral integration of the HAS spectrogram (blue curve) for each delay. The absolute zero delay between gate and pump pulses is identified as the point at which the yield variation is maximized (dashed red line).

Source data

HAS (a) and EUV attosecond streaking (b) spectrograms under identical driving-pulse field. c, Comparison of the vector potential of the gate pulse evaluated from HAS (red) and EUV attosecond streaking (blue) techniques.

Source data

Top panels, HAS spectrograms recorded at different gate/pump intensity ratios, η ≈ 4.5 × 10−3 (a), 8.3 × 10−3 (b), 2.4 × 10−3 (c) and 1.4 × 10−3 (d). Bottom panels, comparison of the vector potential waveform characterized by EUV attosecond streaking (red curve) and HAS (blue curve) for each gate/pump intensity ratio.

Source data

Measured (a) and reconstructed (b) HAS spectrograms in the energy range from 80 to 230 eV. Differential maps of measured (c) and reconstructed (d) spectrograms. e, The phase of the intensity shift in measured (blue) and reconstructed (red) spectrograms. Retrieved terminal electron wave packet in spectral (f) and temporal (g) domains. The purple and red curves denote the intensity and phase, respectively.

Source data

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Kim, H.Y., Garg, M., Mandal, S. et al. Attosecond field emission. Nature 613, 662–666 (2023). https://doi.org/10.1038/s41586-022-05577-1

Download citation

Received: 25 June 2022

Accepted: 18 November 2022

Published: 25 January 2023

Issue Date: 26 January 2023

DOI: https://doi.org/10.1038/s41586-022-05577-1

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Nature (2023)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.