1. INTRODUCTION
Ultrashort and intense laser pulses are particularly well suited for the generation of electron and other charged particle beams, both in the relativistic and in the non-relativistic regimes of laser-solid interactions (Malka, Reference Malka2002; Irvine et al., Reference Irvine, Dechant and Elezzabi2004; Gupta & Suk, Reference Gupta and Suk2007; Karmakar & Pukhov, Reference Karmakar and Pukkov2007; Nickles et al., Reference Nickles, Ter-Avetisyan, Schnuerer, Sokollik, Sandner, Schreiber, Hilscher, Jahnke, Andreev and Tikhonchuk2007; Flippo et al., Reference Flippo, Hegelich, Albright, Yin, Gautier, Letzring, Schollmeier, Schreiber, Schulze and Fernandez2007; Niu et al., Reference Niu, He, Qiao and Zhou2008; Chen et al., Reference Chen, Unick, Vafaei-Najafabadi, Tsui, Fedosejevs, Naseri, Masson-Laborde and Rozmus2008, and references therein). One of the methods to generate well-behaved electron beams with relatively low-intensity light pulses is surface plasmon enhanced electron acceleration. This phenomenon was discovered recently and it was experimentally demonstrated that it is suitable for the production of relatively high-energy, quasi-monoenergetic electron beams with the usage of simple femtosecond lasers (Zawadzka et al., Reference Zawadzka, Jaroszynski, Carey and Wynne2001; Kupersztych et al., Reference Kupersztych, Monchicourt and Raynaud2001; Irvine et al., Reference Irvine, Dechant and Elezzabi2004). In this scheme, the evanescent electric field of surface plasmon polaritons (SPPs) accelerates photo-emitted electrons away from the surface. This process can be highly efficient so that keV kinetic energy levels can be reached with nJ-class oscillators without having to rely on external direct current fields (Irvine et al., Reference Irvine, Dechant and Elezzabi2004; Irvine & Elezzabi, Reference Irvine and Elezzabi2006). This method seems particularly advantageous for the generation of well-behaved femtosecond electron beams that can be used later on for infrared pump/electron probe methods, such as ultrafast electron diffraction or microscopy (Siwick et al., Reference Siwick, Dwyer, Jordan and Miller2003; Lobastov et al., Reference Lobastov, Srinivasan and Zewail2005). These novel time-resolved methods utilizing electron beams can gain importance in the future by enabling both high spatial and high temporal resolution material characterization at the same time. They will become particularly interesting if the attosecond temporal resolution domain becomes within reach, as suggested recently (Stockman et al., Reference Stockman, Kling, Krausz and Kleineberg2007; Varró & Farkas, Reference Varró and Farkas2008; Fill et al., Reference Fill, Veisz, Apolonski and Krausz2006). Moreover, studying the spectral properties of femtosecond electron bunches has the potential to reveal ultrafast excitation dynamics in solids and provide the basis for a single-shot measurement tool of the carrier-envelope phase of ultrashort laser pulses, as we proposed recently (Irvine et al., Reference Irvine, Dombi, Farkas and Elezzabi2006; Dombi & Rácz, Reference Dombi and Rácz2008). Other, carrier-envelope phase sensitive laser-solid interactions already demonstrated (Apolonski et al., Reference Apolonski, Dombi, Paulus, Kakehata, Holzwarth, Udem, Lemell, Torizuka, Burgdörfer, Hänsch and Krausz2004; Dombi et al., Reference Dombi, Apolonski, Lemell, Paulus, Kakehata, Holzwarth, Udem, Torizuka, Burgdörfer, Hänsch and Krausz2004, Reference Dombi, Krausz and Farkas2006; Fortier et al., Reference Fortier, Roos, Jones, Cundiff, Bhat and Sipe2004; Mücke et al., Reference Mücke, Tritschler, Wegener, Morgner, Kärtner, Khitrova and Gibbs2004) suffer from low experimental contrast, therefore it is necessary to look for higher contrast tools for direct phase measurement.
Motivated by these perspectives, it was shown numerically and partly also experimentally that surface plasmonic electron sources can be ideally controlled with ultrashort laser pulses so that they deliver highly directional, monoenergetic electron beams readily synchronized with the pump pulse (Irvine et al., Reference Irvine, Dechant and Elezzabi2004, Reference Irvine, Dombi, Farkas and Elezzabi2006; Dombi & Rácz, Reference Dombi and Rácz2008). We developed a simple semiclassical approach for the simulation of this process that is analogous to the three-step model of high harmonic generation (Kulander et al., Reference Kulander, Schafer and Krause1993; Corkum, Reference Corkum1993). In this paper, we present the basic elements of this model and show that it provides the same results as a much more complicated treatment of the problem based on the rigorous, but computationally time-consuming solution of Maxwell's equations. We also provide new insight into the spatiotemporal dynamics of SPP enhanced electron acceleration which is also important if one intends to realize adaptive emission control methods (Aeschlimann et al., Reference Aeschlimann, Bauer, Bayer, Brixner, García de Abajo, Pfeiffer, Rohmer, Spindler and Steeb2007).
2. THREE-STEP MODEL OF SPP ENHANCED ELECTRON ACCELERATION
SPP enhanced electron acceleration involves distinct physical processes such as (1) the coupling of the incident light and surface plasmonic electromagnetic fields, (2) the photoinjection of the electrons into vacuum from the metal layer, and (3) the subsequent acceleration of free electrons by the evanescent SPP field on the vacuum side of the surface. These steps represent an analogy with the semiclassical three-step model of high harmonic generation that can be used very efficiently to predict the outcome of atomic recollision processes (Kulander et al., Reference Kulander, Schafer and Krause1993; Corkum, Reference Corkum1993). This motivated us to adapt a more or less analogous model to the SPP environment where instead of a single atom, a solid surface is involved, which determines the conditions for recollision.
As a first step, we gained analytic formulae for SPP fields instead of the computationally intensive complete numerical solution of Maxwell's equations in the Kretschmann-Raether SPP coupling configuration used in previous studies (Irvine et al., Reference Irvine, Dechant and Elezzabi2004). Based on the well-known fact that SPP fields decay exponentially away from the surface (Raether, Reference Raether1988) we took an analytic expression for the SPP field components on the vacuum side of the metal layer in the form of
where E 0 is the field amplitude, Eenv(x, t) is an envelope function determined by the temporal and spatial beam profiles of the incoming Gaussian pulse, η is the field enhancement factor resulting from SSP coupling (Raether, Reference Raether1988), k SPP is the SPP wave vector, ω0 is the carrier frequency, ϕ0 is the carrier-envelope phase of the laser pulse, and α is the inverse decay length of the plasmonic field in vacuum. For the accurate determination of the field, we used the evanescent decay parameter α = 247 nm−1 resulting from previous non-approximate studies carried out with the same input parameters for laser pulses having a central wavelength of 800 nm (Irvine et al., Reference Irvine, Dombi, Farkas and Elezzabi2006). We used the value of a = 0.3 according to the notion that the amplitudes of the x- and y-components of the plasmonic field have this ratio according to the numerical solution of Maxwell's equations (Irvine et al., Reference Irvine, Dombi, Farkas and Elezzabi2006). We concluded that the field given by Eqs. (1) and (2) approximates the exact SPP field with a very good accuracy by comparing our results to those of Irvine et al. (Reference Irvine, Dombi, Farkas and Elezzabi2006). One can examine the distribution of the field amplitude in the vicinity of the surface in Figure 1 (false color representation) showing very good agreement with the above-mentioned calculation and we also succeeded in reproducing the vector representation of the field depicted in Figure 3 in Irvine et al. (Reference Irvine, Dombi, Farkas and Elezzabi2006) with this method. The vector field that can be calculated with our model is depicted in the inset of Figure 1.
(Color online) Illustration of the setup for the generation of electron beams by surface plasmon enhanced electron acceleration with field vectors and electron trajectories. See text for further details.
We then placed a point array along the prism surface and examined the spatial and temporal distribution of the photoemission (induced by the SPP field) along the surface, assuming field emission taking place at higher intensities. To this end, we applied the Fowler-Nordheim equation routinely used in studies involving electron emission from metal nanotips (Hommelhoff et al., Reference Hommelhoff, Kealhofer and Kasevich2006; Ropers et al., Reference Ropers, Solli, Schulz, Lienau and Elsaesser2007). This describes the instantaneous tunneling current based on the fact that plasmonic fields carry substantial field enhancement factors (up to ×100) compared to the generating field (Hommelhoff et al., Reference Hommelhoff, Kealhofer and Kasevich2006; Ropers et al., Reference Ropers, Solli, Schulz, Lienau and Elsaesser2007). This way we gained a spatially and temporally resolved map of tunneling probabilities determined by the SPP field.
Third, we scrutinized the vacuum electron trajectories of photoemitted electrons in the plasmonic field for each point in the above-mentioned array and for several emission instants by solving free-electron equations of motions numerically. Some of these trajectories can be seen in Figure 1 (red curves). In some cases, the electron trajectories involved a recollision with the metal surface and when this happened, no electron emission was assumed. In all other cases, the final kinetic energies and directions of the photoemitted and photoaccelerated electrons were placed in a matrix for each emission point in space and for each emission instant. Figure 2 illustrates the temporal distribution of the final kinetic energies as a function of the electron “birth” instant for a maximum plasmonic field strength of 1.9 × 1011 V/m and for electrons emitted from the central part of the illuminated surface. The figure demonstrates similarities to the corresponding kinetic energy distributions of atomic electrons after being accelerated by the ionizing laser field (Reider, Reference Reider2004). Here, it is important to note that roughly only one-fourth of the possible emission instants can contribute to the acceleration process. This is due to the symmetry breaking of the metal surface and the associated electron recollision and reabsorption processes.
(Color online) Surface plasmon accelerated electron energy as a function of the birth instant of the electrons (scatter plots). The plasmon generating 5 fs laser pulse (solid and dashed lines) has either a “cosine” (dashed curve) or “minus cosine” waveform (solid curve). Electron energies for the “cosine” waveform are depicted as circles, whereas for the “minus cosine” waveform, as squares. For further pulse parameters see text.
3. DETERMINING MACROSCOPIC PROPERTIES OF THE ELECTRON BEAM
By integrating the above described emission maps along the spatial and/or temporal coordinates, macroscopic emission distributions and electron spectra can be calculated. First, we checked whether these reproduce former measurement and simulation results (Irvine et al., Reference Irvine, Dechant and Elezzabi2004, Reference Irvine, Dombi, Farkas and Elezzabi2006) to gain confidence in the simplified three-step model that we used. To this end, we carried out simulations for the same parameters as those published in Irvine and Elezzabi (2006a) and Irvine et al. (Reference Irvine, Dombi, Farkas and Elezzabi2006). For these simulations, we assumed multi-photon-induced electron emission (previously used in other simulations), which does not necessarily hold for higher intensities. However, our purpose in this case was to reproduce former results, therefore the temporal distribution of photoemission was described by j(t, x) ~ I n(t, x), where j is the photocurrent from the surface, I is the temporal intensity distribution of the SPP field along the surface, and n is the order of the multi-photon process. Here, n = 3 is used according to the 4.5–5 eV work function of most metals and the 1.5 eV photon energy at 800 nm. Figure 3a depicts macroscopic electron spectra gained with our model for peak plasmonic fields of 1.9 × 1011 V/m, 2.7 × 1011 V/m, and 3.7 × 1011 V/m, respectively (the full width at half maximum (FWHM) duration of the input Gaussian laser pulse was 30 fs with a central wavelength of 800 nm). Thereby, this figure can be directly compared to Figure 5 in Irvine et al. (Reference Irvine, Dombi, Farkas and Elezzabi2006). It can be seen that the characteristics of the electron spectra are very well reproduced, as well as the linear scaling of the kinetic energies of the most energetic electrons with intensity. Slight differences in the peak and cut-off positions can be attributed to the approximate nature of the SPP field expression used in our case as opposed to the numerical field solution used in Irvine et al. (Reference Irvine, Dombi, Farkas and Elezzabi2006).
(Color online) (a) Macroscopic electron spectra at peak plasmonic fields of 1.9 × 1011 V/m (solid line), 2.7 × 1011 V/m (dashed line), and 3.7 × 1011 V/m (dotted line) for a Gaussian input laser pulse of 30 fs FWHM duration with a central wavelength of 800 nm). (b) Electron spectrum for a 5 fs generating pulse with a peak plasmonic field strength of 1.9 × 1011 V/m, assuming multi-photon photoemission.
In another comparative simulation, we changed the input pulse length to 5 fs FWHM, and assumed that this pulse is focused to a spot with 60 µm diameter on the prism surface. The field peak was 1.9 × 1011 V/m. In Figure 3b, we can see that the spectrum of the electron beam gained this way reproduces the spectrum computed with other methods, such as the one in Figure 1b in Irvine et al. (Reference Irvine, Dombi, Farkas and Elezzabi2006). Slight differences in the cut-off positions can still be observed, however, all spectral features and the position of the main peak are exactly the same. Therefore the applicability of the analytic field expressions (Eqs. (1) and (2)) and the robustness of our approach was confirmed by these examples.
We can now turn attention to modeling electron spectra by assuming field emission from the metal surface, which is a more realistic assumption for higher intensity input beams, approaching the damage threshold of thin metal films. High repetition-rate, ultrafast laser output delivering focused intensity in this range is achievable with simple Ti:sapphire oscillators with an extended cavity, as we demonstrated it recently (Naumov et al., Reference Naumov, Fernandez, Graf, Dombi, Krausz and Apolonski2005; Dombi & Antal, Reference Dombi and Antal2007; Dombi et al., Reference Dombi, Antal, Fekete, Szipöcs and Várallyay2007). We then used the Fowler-Nordheim formula, as described above, and both angularly and spectrally resolved the photoaccelerated electron beam, assuming a maximum input field of 5.8 × 1010 V/m, which is a rather realistic maximum value considering the damage threshold of gold and silver films. We also assumed a tunneling time of 600 attosecond, which in our model, describes the delay between the actual distortion of the potential by the field and the corresponding appearance of the electron in the continuum. The simulation results can be seen in Figure 4a with false color representation of the electron count value as a function of energy and emission angle. It can be seen that the electron beam is highly directional, most of it being emitted to the 95°–105° angular range. The beam has also acquired a quasi-monoenergetic feature as it can be seen in Figure 4b. This spectrum is generated by integrating the angular-energy distribution along the angular axis. The quasi-monoenergeticity of the spectrum is a significant difference as compared to the case in Figure 3b, and can be attributed exclusively to the different emission regimes (multiphoton versus tunneling) involved. The sharp temporal distribution of the tunneling peaks located at the field maxima favor the emission of electrons at those time instants when they can gain significant kinetic energy. The sharp spectral drop corresponds to the highest energy electrons in the multiphoton case. But it is mostly these electrons that are represented in the field emission case; therefore, it appears as a sharp cut-off in the spectrum. The low-energy wing of the spectrum, however, displays a broader feature, making the source less suitable for ultrafast applications. Figures 4c and 4d display calculations assuming a tunneling time of 150 as, otherwise, the parameters are unchanged. It can be seen that high-energy electrons are considerably suppressed and the advantageous quasi-monoenergetic feature of the spectrum also disappears. Hence, the actual value of the tunneling time (depending on the field strength, the work function, and the length of the optical cycle) has a decisive influence on the acceleration process. As it is a significant challenge to measure this delay in the field emission process (Eckle et al., Reference Eckle, Pfeiffer, Cirelli, Staudte, Dörner, Muller, Büttiker and Keller2008) and in spite of the rapid development of attosecond techniques, there is no corresponding measurement for metals that we are aware of, we can only state that this value lies closer to one-fourth of the optical cycle duration for the field strengths concerned. This value can be derived with density functional methods (Lemell et al., Reference Lemell, Tong, Krausz and Burgdörfer2003) as well as from simple quantum mechanical considerations.
(Color online) Angular and kinetic energy distribution (a and c) and macroscopic electron spectrum (b and d) of surface plasmon accelerated electrons for 5 fs FWHM laser pulses and a maximum field strength of 5.8 × 1010 V/m, assuming tunneling emission from the surface with 600 as (a and b) and 150 as delay (c and d) with respect to the SPP electric field.
To generate spectra with higher monoenergeticity, we suggest the application of spatial confinement of the emission area on the metal surface. Experimentally, this can be carried out by various nanofabrication techniques, e.g., by depositing a dielectric layer on top of the metal with a nanoscale opening where the dielectric overlayer is absent and the metal surface is exposed to vacuum. Another possibility is roughening a small rectangular area on top of the metal surface thereby enhancing the emission from that portion of the film. These potential schemes were taken into account in our simulations by selecting only smaller areas of the surface illuminated by the laser beam and considered only those photoelectrons that were emitted from this area. Then, we determined macroscopic spectra as a function of the radius of the nano-opening, and depicted these normalized distributions in false color representation in Figure 5. One sees that limiting the emission area clearly improves the monoenergetic nature of the electron beam, as depicted in Figure 5. According to the figure, we can determine that the ideal size of the opening is around 200–500 nm. The spectral monoenergeticity is highly improved in this range and the escape area is not too small to allow a significant photoelectron escape.
(Color online) Kinetic energy distribution of SPP accelerated electrons emitted from a nano-confined window for a maximum field strength of 5.8 × 1010 V/m as a function of the size of the window (y-axis, in meters). See text for further details.
4. SUMMARY
We modeled few-cycle pulse-induced, SPP enhanced electron acceleration with a semiclassical approach similar to the three-step model of high harmonic generation. Our analytic approximation for the plasmonic field enables the reproduction of former results gained by the direct, but very processor time consuming solution of Maxwell's equations. We have shown that this phenomenon can serve as a basis for an ultrafast, monoenergetic electron source with a spectrum peaking at hundreds of eVs. It is anticipated that even higher degrees of monoenergeticity can be reached by applying a nano-sized window on the emission area. Thereby, the surface plasmonic acceleration scheme can provide a solution for novel time-resolved diffraction methods and serve as a basis for time-resolved studies of ultrafast surface processes.