We propose to use two lasers of the same frequency to achieve interference modulation of photoelectron emission. Using a quantum mechanical model, we study the time-resolved photoelectron energy spectra and emission current modulation under different laser and dc fields. We find that strong interference modulation of photoemission can be easily achieved with two lasers of the same frequency, due to, on one hand, the straightforward access to the single-frequency laser pairs in experiments and, on the other hand, the low threshold value of the ratio of the laser fields for large modulation depth even with a strong dc field. Our study demonstrates the capability of using interference modulation by single-frequency laser pairs for practical measurements of time-resolved photoelectron energy spectra.
Ultrafast electron emission from nanostructures driven by lasers is of large interest, due to its wide applications, such as time-resolved photoelectron microscopy,1–5 electron diffraction imaging,6–9 ultrafast electron sources,10–13 free-electron lasers,14,15 attosecond electronic devices,16–18 and emerging vacuum nanodevices.19–21 Recently, two-color laser-induced electron emission from nanoemitters22–28 has drawn strong interest as it provides an attractive platform for modulating photoelectron emission by the relative phase difference between the two lasers. It also shows promise for the potential application of time-resolved photoelectron spectroscopy.28,29 However, the two-color laser system typically relies on the generation of higher order harmonics of a fundamental laser,22–26 which, because of its stringent requirements on the experimental setup and its relatively low efficiency, greatly limits the accessibility of the two-color laser system. For higher intensity lasers, harmonic generation becomes increasingly complex and difficult to control.30,31
In this work, we propose to utilize two lasers of the same frequency to modulate the photoelectron emission by their relative phase delay. This is motivated by the simple experimental implementation of single-frequency laser pairs, e.g., via a beam splitter with various coating materials to control the reflection and transmission of incident light.32–35 The two same-frequency lasers may be tuned to have a virtually arbitrary ratio of intensities (in contrast to a small harmonic-to-fundamental intensity ratio in the two-color laser system22–26), thus providing a much larger parameter space to assess the interference effect of the two lasers and the induced photoelectron emission. Using a quantum mechanical model,36 we study the photoemission modulation properties for a dc-biased metal cathode illuminated by two laser fields with the same frequency. We investigate the modulation of photoemission current and the dynamics of multiphoton excited states for different laser fields, wavelengths, cathode materials, and dc bias. Our study demonstrates the capability of measuring the time-resolved photoelectron energy spectra using single-frequency laser pairs.
Under the action of two laser fields F1cos(ωt) and F2cos(ωt + θ) and a dc electric field F0, the time-dependent potential barrier near the surface of the cathode reads27,28,36,37
where EF is the Fermi energy of the metal cathode, is the effective work function with the Schottky effect,36 with W being the nominal work function, e is the elementary charge, ε0 is the free space permittivity, x is the distance away from the cathode surface (x = 0), and F is the magnitude of the total laser field due to the two laser fields F1cos(ωt) and F2cos(ωt + θ),
From Eq. (2), it is clear that the magnitude of the total laser field depends strongly on the phase delay of the two lasers θ, which is expected to provide a similar current modulation to that in the two-color laser setup.27,28 The resultant phase φ = arcsin(F2 sin θ/F), the effect of which becomes important for photoemission only in very short laser pulses when the carrier–envelope phase matters. For laser pulses longer than about 10 cycles, it can be well approximated by continuous-wave excitation for photoemission.36 Thus, in this study, we ignore the effects of the absolute phase and set φ = 0 without loss of generality. Based on the exact quantum theory of photoemission,36,37 the time-averaged normalized emission current density, defined as the time-averaged ratio of the transmitted probability current density over the incident probability current density, ⟨w(ε, x, t)⟩ = ⟨Jt/Ji⟩, can be obtained as
where ⟨wn⟩ denotes the normalized emission current density through the nth channel with emitted electron energy ε + nℏω due to the n-photon contribution, ℏ is the reduced Planck constant, m is the electron mass, and Tn represents the transmission coefficient of electron wave functions, which is calculated from
where δ(l) is the Dirac delta function, l and n are integers, and are the Fourier coefficients, with and , where r(αn) = Ai(αn) − iBi(αn), , and with En = ε + nℏω − EF − Weff − Up. Here, Ai and Bi are the Airy functions of the first kind and second kind, respectively, Up = e2F2/4mω2 is the ponderomotive energy, and a prime denotes the derivative with respect to the argument. For the special case of zero dc field F0 = 0, the time-averaged normalized emission current density becomes36,37
where Tn is still calculated from Eq. (4) with Pnl and Qnl unchanged, but with and . These analytical results in Eqs. (3)–(5) are obtained by solving the time-dependent Schrödinger equation exactly with the potential energy given in Eq. (1) (see Ref. 36 for the detailed derivation). We would like to point out that although the calculation here is equivalent to photoemission using a single laser with varying intensity, our proposed method of using two lasers provides more degrees of freedom for experimental control and also offers the capability of measuring the emitted electron spectra vs time after emission, i.e., time-resolved photoelectron energy spectra.
In Fig. 1, we plot the calculated photoelectron energy spectra as a function of the phase difference between the two lasers θ for different dc fields F0. The wavelength of both lasers is 800 nm (ℏω = 1.55 eV). The metal is assumed to be tungsten,10,17,22 with a Fermi energy EF = 7 eV and a work function W = 4.31 eV. Since most of the electrons emitted from sources are located near the Fermi level,36,38–41 we choose the electron initial energy ε = EF for simplicity. Unless mentioned otherwise, these are the default values for the calculations in this paper. Note with laser fields F1 = 1.8 V/nm and F2 = 0.3 V/nm for the special case of θ = 0, the total normalized emission current density in Fig. 1 is ⟨w⟩ = 6.67 × 10−7 and 8.71 × 10−5, for the DC field F0 = 0 V/nm and 0.8 V/nm, respectively. Using free-electron theory of metal,37 we find that the corresponding emission current density is 5.74 × 102 A/cm2 and 6.75 × 104 A/cm2, respectively.
Time-resolved photoelectron energy spectra. (a) and (b) Energy spectra as a function of the phase difference between the two lasers θ, for dc field (a) F0 = 0 and (b) F0 = 0.8 V/nm. (c) and (d) Projections of the spectra in (a) and (b) on the θ–⟨w⟩n plane, respectively. (e) and (f) Projections of the spectra in (a) and (b) on the n–⟨w⟩n plane, respectively. Here, the laser fields F1 = 1.8 V/nm and F2 = 0.3 V/nm (experimental parameters in Ref. 22).
Time-resolved photoelectron energy spectra. (a) and (b) Energy spectra as a function of the phase difference between the two lasers θ, for dc field (a) F0 = 0 and (b) F0 = 0.8 V/nm. (c) and (d) Projections of the spectra in (a) and (b) on the θ–⟨w⟩n plane, respectively. (e) and (f) Projections of the spectra in (a) and (b) on the n–⟨w⟩n plane, respectively. Here, the laser fields F1 = 1.8 V/nm and F2 = 0.3 V/nm (experimental parameters in Ref. 22).
When the dc field F0 is turned off, the dominant emission process is three-photon absorption (n = 3) [cf. Figs. 1(a) and 1(e)]. This is consistent with the ratio of the work function of tungsten over the photon energy, W/ℏω ≈ 2.8. By changing the phase difference θ between the two lasers, the electron emission varies sinusoidally [cf. Figs. 1(a) and 1(c)]. When applying a large dc field F0 to the cathode, the tunneling emission channels (n ≤ 2) are opened up, as shown in Figs. 1(b) and 1(f). This is because the dc field adequately narrows the surface potential barrier, in addition to the Schottky effect induced barrier lowering, enabling the tunneling emission process. In the case of F0 = 0.8 V/nm, the dominant emission process is shifted to two-photon absorption. From Figs. 1(c) and 1(d), it is found that the multiphoton excited states (n ≥ 3) vary with respect to the phase delay θ sinusoidally in the same way, with the maximum at θ = 0 and the minimum at θ = ±π, for both values of dc bias F0. This is in contrast to the two-color laser induced photoemission, where the dynamics of multiphoton excited states changes under different dc bias F0 [cf. Figs. 10(g) and 10(f) in Ref. 28]. The one-photon (n = 1) absorption and direct tunneling (n = 0) process are almost independent of θ for the case of F0 = 0.8 V/nm, as shown in Figs. 1(d) and 1(f). The sinusoidal modulation in the total emission current density ⟨w⟩ is shown in Figs. 2(a) and 2(b), for the case of F0 = 0 V/nm. When the laser field ratio F2/F1 increases, the maximum emission current ⟨w⟩max at θ = 0 increases, while the minimum emission current ⟨w⟩min at θ = π decreases, due to the more profound interference of the two lasers. Figure 2(c) shows the modulation depth, Γ=(⟨w⟩max − ⟨w⟩min)/(⟨w⟩max + ⟨w⟩min), as a function of laser field ratio F2/F1 under different dc fields F0. For a given F0, Γ increases as F2/F1 increases, and it reaches the maximum value of 100% when F1 = F2. It is important to note that, in order to reach a large modulation depth (Γ ≥ 90%), only a small laser field ratio F2/F1 is needed even with a strong dc field, e.g., F2/F1 ≤ 0.4 when F0 = 1 V/nm. The dependence of Γ on the dc field F0 [Fig. 2(c)] is not monotonic and will be examined further in Fig. 3.
Photoemission current modulation. (a) Normalized total time-averaged emission current density ⟨w⟩ as a function of phase difference θ for different F2/F1, when the dc field F0 = 0. (b) Semilog plot of ⟨w⟩ in (a). (c) Current modulation depth Γ (solid lines) as a function of the laser field ratio F2/F1 for different dc fields F0. Dotted (dashed) lines in (c) are for the maximum (minimum) emission current density ⟨w⟩ at θ = 0 (θ = π). Here, the laser field F1 is fixed as 1.8 V/nm.
Photoemission current modulation. (a) Normalized total time-averaged emission current density ⟨w⟩ as a function of phase difference θ for different F2/F1, when the dc field F0 = 0. (b) Semilog plot of ⟨w⟩ in (a). (c) Current modulation depth Γ (solid lines) as a function of the laser field ratio F2/F1 for different dc fields F0. Dotted (dashed) lines in (c) are for the maximum (minimum) emission current density ⟨w⟩ at θ = 0 (θ = π). Here, the laser field F1 is fixed as 1.8 V/nm.
(a) Emission current modulation depth Γ (solid lines) as a function of the dc field F0 with and without the image-charge-induced potential barrier lowering (or the Schottky effect), for laser fields F1 = 1.8 V/nm and F2 = 0.3 V/nm. The case without the Schottky effect is calculated by replacing Weff with the nominal work function of metal W in Eq. (1). (b) Modulation depth Γ (solid lines) as a function of F0 for different laser field ratios F2/F1, with the effective work function Weff. F1 is fixed at 1.8 V/nm in (b). In (a) and (b), the dotted (dashed) lines are for the maximum (minimum) emission current density ⟨w⟩ at θ = 0 (θ = π).
(a) Emission current modulation depth Γ (solid lines) as a function of the dc field F0 with and without the image-charge-induced potential barrier lowering (or the Schottky effect), for laser fields F1 = 1.8 V/nm and F2 = 0.3 V/nm. The case without the Schottky effect is calculated by replacing Weff with the nominal work function of metal W in Eq. (1). (b) Modulation depth Γ (solid lines) as a function of F0 for different laser field ratios F2/F1, with the effective work function Weff. F1 is fixed at 1.8 V/nm in (b). In (a) and (b), the dotted (dashed) lines are for the maximum (minimum) emission current density ⟨w⟩ at θ = 0 (θ = π).
Besides making the surface potential barrier narrower, the dc bias induces a reduction of the barrier height via the image charge effect (or the Schottky effect), which strongly influences the photoemission processes.28,36 In Fig. 3(a), we compare the emission current modulation depth Γ as a function of the dc field F0 with and without the Schottky effect. When the Schottky effect is not considered, Γ gradually decreases with F0. It is clear that the Schottky effect greatly alters the dependence of modulation depth Γ on the dc field F0 [cf. solid lines in Fig. 3(a)]. The change in Γ originates from the change in the maximum (minimum) values of emission current ⟨w⟩ with the Schottky effect, as shown by dotted (dashed) lines in Fig. 3(a). As F0 varies, the effective potential barrier Weff changes, which induces an increase (decrease) in the emission current when the ratio Weff/ℏω becomes closer to (further away from) an integer, where resonant n-photon absorption occurs (cf. Figs. 6 and 7 of Ref. 36). This resonant emission process causes the nonlinear behavior of Γ as a function of dc field F0.
Figure 3(b) shows the modulation depth Γ as a function of the dc field F0 for different laser field ratios F2/F1 with fixed F1 = 1.8 V/nm. As F2/F1 approaches 1, the modulation depth Γ gradually approaches the maximum value of 1 for the full range of dc field F0 from 0 V/nm to 1 V/nm. This is consistent with the observation in Fig. 2(c). Note that when F0 is increased from 0 V/nm to 1 V/nm, the total emission current density can be increased by orders of magnitude.
We also examine the photoemission current modulation depth Γ for cathode materials with different work functions in Fig. 4(a) and for various incident laser wavelengths in Fig. 4(b). We fix the dc field F0 = 0.8 V/nm and laser fields F1 = 1.8 V/nm and F2 = 0.3 V/nm. Under the same illumination condition, the electron emission current depends strongly on the work function; however, the modulation depth Γ varies only slightly. This is because Γ is predominantly determined by the ratio of the laser field strengths. Figure 4(b) shows the effect of the laser wavelength on both emission current and modulation depth for a tungsten cathode. The nonlinear dependence may also be attributed to the change in the ratio Weff/ℏω near resonant n-photon processes.36
Normalized total time-averaged emission current density ⟨w⟩ as a function of the phase difference θ, for various (a) cathode materials and (b) incident wavelengths. In (a), the laser wavelength λ = 800 nm (ℏω = 1.55 eV). The nominal work function of different materials is WAg = 4.26 eV,42 Ww = 4.31 eV,22,25 WMo = 4.6 eV,42 WCu = 4.65 eV,42 and WAu = 5.1 eV.36,42 In (b), the metal is tungsten. Here, the dc field F0 is 0.8 V/nm, and the laser fields F1 and F2 are fixed at 1.8 V/nm and 0.3 V/nm, respectively.
Normalized total time-averaged emission current density ⟨w⟩ as a function of the phase difference θ, for various (a) cathode materials and (b) incident wavelengths. In (a), the laser wavelength λ = 800 nm (ℏω = 1.55 eV). The nominal work function of different materials is WAg = 4.26 eV,42 Ww = 4.31 eV,22,25 WMo = 4.6 eV,42 WCu = 4.65 eV,42 and WAu = 5.1 eV.36,42 In (b), the metal is tungsten. Here, the dc field F0 is 0.8 V/nm, and the laser fields F1 and F2 are fixed at 1.8 V/nm and 0.3 V/nm, respectively.
In summary, we propose to utilize two lasers of the same frequency to modulate the photoelectron emission by their phase delay. Compared to the two-color laser configuration, single-frequency laser pairs can be more easily implemented in experiments since they relax the requirement of higher order harmonic generation, which becomes increasingly difficult in the high laser intensity regimes. The intensity ratio of the single-frequency laser pairs can be tuned over a much wider range than the two-color laser system. We find that a strong current modulation (>90%) can be achieved with a moderate ratio of the laser fields (<0.4) even under a strong dc bias. The nonlinear effects of the dc field, cathode materials, and laser wavelength on both the emission current level and modulation depth are examined. The strong dependence of photoelectron energy spectra on the phase delay of the two lasers demonstrates promising potential for the application of time-resolved photoelectron spectroscopy using single-frequency laser pairs. Our study may inspire new routes toward many applications requiring both high photoemission current and strong current modulation, such as tabletop particle accelerators, photoelectron microscopy, and x-ray sources. Future research may consider the effects of cathode geometry, surface states and materials (e.g., semiconductor and two-dimensional materials), and different polarization and shape of laser pulses on photoelectron emission modulation. Further optical engineering techniques, such as those which may be provided by spatial light modulators (SLMs), might allow synthesis of even more complex current waveforms in a compact experimental package. As SLMs have been demonstrated to operate at both high power and high energy, implementation in a variety of laser systems is possible.43
DATA AVAILABLITY
The data that support the findings of this study are available within the article.
ACKNOWLEDGMENTS
This work was supported by Air Force Office of Scientific Research (AFOSR) YIP Grant No. FA9550-18-1-0061 and Office of Naval Research (ONR) YIP Grant No. N00014-20-1-2681.