Abstract
We experimentally investigate the propagation of Airy surface plasmon polaritons (Airy SPPs) on a gold film by multiphoton Photoemission Electron Microscopy (PEEM) at different excitation wavelengths and compare the result with rigorous numerical simulations. The typical bent trajectories of the excited two-dimensional beams are observed and analyzed over a wide range of wavelengths. We furthermore investigate the generation bandwidth of the diffraction grating from modal overlap calculations and evaluate the possibility of creating ultrashort Airy plasmon pulses. This provides a viable route to engineer two-dimensional ultrashort non-diffracting pulsed beams in the field of ultrafast nanophotonics.
© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
The concept of Airy waves was first introduced in the seminal work of Berry and Balazs [1]. Thereafter, it became popular in various fields of physics including optics [2–7]. Airy beams are a solution of the 1D paraxial wave equation that describes the propagation of light in a homogeneous dielectric medium, where diffraction is limited to one transverse spatial dimension. Such Airy beams are self accelerating, non-diffracting beams that travel along a parabolic trajectory in a medium without any refractive index gradient. They have the interesting property to regain their Airy wave profile even after colliding with a finite sized obstacle. In contrast to other non-diffracting beams such as transverse Bessel, Mathieu, and Weber beams [8], Airy beams can exist only in a planar system, which renders them particularly attractive for plasmonics or other systems supporting surface waves.
Surface plasmon polaritons (SPPs) are coupled electromagnetic oscillations of photons and electrons, which are bound to a metal-dielectric interface [9]. SPPs with the properties of Airy beams have been proposed and experimentally verified [10–13]. These so called Airy surface plasmon polaritons (Airy SPPs or just Airy plasmons) propagate diffraction free along a parabolic trajectory on a metal-dielectric interface.
2. Experimental characterization of Airy surface plasmons using PEEM
In our work, we have investigated the spatial dynamics of Airy plasmons using multiphoton Photoemission-Electron Microscopy (PEEM) over a wide range of wavelengths. PEEM is an experimental tool to directly visualize and measure the electromagnetic fields at a metal surface through photoemission [14,15]. The detected photoemitted electrons provide a map of the electromagnetic field at the surface. PEEM can be further combined with ultrafast laser excitation to study the spatio-temporal evolution of pulsed Airy plasmons. In this context, we have further investigated the propagation dynamics of pulsed Airy plasmons and their temporal evolution in space based on a numerical model. This concept has potential applications in ultrafast nanophotonics [16] and can be employed in combination with 2D materials [17,18]. One such application can be pulsed Airy plasmons on graphene surface in which active control over the propagation dynamics can be achieved by tuning the Fermi energy of graphene [19].
A diffraction grating (Fig. 1) based on the idea of Minovich et al [12] has been used to excite the Airy plasmons. The grating couples the exciting plane wave, which is incidenting from the free space on the grating, to SPPs and simultaneously imprints the Airy wave profile. In z direction, the grating structure consists of 11 periods of slit arrays whose period $p$ in the z direction has been adjusted to 745 nm. The slit size in the z direction is fixed to 200 nm. In x direction the slit size and therefore also the local period varies and is determined by the distance between the roots of the Airy function. For the slits of the first column the size in x direction is set to 2${x_0}$, where ${x_0}$ is the full width of the main lobe of the Airy profile and is set to $x_0= 0.7\mu$m. The grating pattern was milled by a Focused Ion Beam into a 200 nm thick polycrystalline gold film, which was deposited on a fused silica substrate by evaporation (see inset in Fig. 2(a) for an SEM micrograph).
We used a PEEM from Focus GmbH, Germany to study the spatial dynamics of plasmon polariton propagation from the grating by imaging the photoelectrons emitted from the multiphoton interference of the plasmon polaritons with the exciting electromagnetic field. For the excitation light source we used a home-built tunable optical parametric chirped pulse amplifier, which provides an average output power of 100 mW at the sample position with a repetition rate of 1 MHz, corresponding to a pulse energy of $0.1\,\mu$J. A TM polarized wide Gaussian beam of diameter $80\,\mu$m is used to illuminate the grating with it’s wavevector placed in yz-plane and having an angle $4^{\circ }$ with the y axis. The central wavelength of the excitation was varied between 670 nm and 840 nm. The pulses had a spectral bandwidth of 15 nm and a Gaussian envelope with a pulse length of 500 fs [FWHM]. The total exposure time for one image was 20 minutes. The rate of emitted electrons per laser pulse is small enough in order to suppress the detrimental interaction between the photoemitted electrons and was optimized to give the best signal-to-noise ratio in the experiments.
Large scale 3D simulations have been performed to optimize and theoretically verify the design. Computation is done using a Finite Difference Time Domain (FDTD) method [20]. We have used the commercially available Lumerical FDTD solver [21] for our simulation purpose. We have used a TM polarized plane wave inclined by $4^{\circ }$ from the normal in the yz-plane and the fields are extracted in air 5 nm above the gold-air interface. The dielectric permittivity of gold is taken from Johnson and Christy [22]. The plots in Fig. 2 are shown on a logarithmic scale. Qualitatively, there is a good agreement between the electron yield distribution in the experiment and the simulated total light intensity distribution. Both plots capture similarly the details of the generated Airy plasmons and the interference pattern with inherent phase information.
Figure 2 displays the qualitative comparison of experimental findings and numerically simulated results for an excitation wavelength of $\lambda =700$ nm. An overlaid SEM micrograph of the fabricated grating structure with the experimentally obtained image by the PEEM is shown in Fig. 2(a). Airy plasmons are launched and propagate along the air-gold interface. The generated plasmons simultaneously interfere with the incident inclined plane wave at $4^{\circ }$ with respect to surface normal in the yz-plane. The resulting field is mapped by the photoemitted electrons from the gold surface. The measured electron yield by the multiphoton PEEM identifies the regions of constructive and destructive interference and provides a clear contrast between the higher and lower photoemission yield. The map clearly shows the non-diffracting, self accelerating property of the Airy plasmons. The main lobe of the Airy plasmon propagates along a parabolically curved trajectory over a distance of at least $20\,\mu$m.
Furthermore, it is important to comprehend the effect of a change in the excitation wavelength on the propagation dynamics. The wavelength dependence of the Airy plasmon’s excitation and the spatial propagation dynamics is shown in Fig. 3 over the wavelength range from 670 nm to 770 nm at 20 nm wavelength intervals for the same excitation grating as in Fig. 2. The qualitative comparison of the measured PEEM yield and numerical simulations shows a good agreement. The wavelength tuning changes the trajectory of the generated plasmon polaritons.
Experiments and simulations even agree in minor details. As an example a careful observation near the first column of the grating shows the generation of a weak lobe at a wavelength 670 nm (indicated by an arrow (Fig. 3(a), (d))). It is a signature of a localized second order resonance in the individual slits of the first column of the excitation grating. This second order resonance diminishes for longer excitation wavelengths. We also observe a spatial shift of the Airy plasmon generation towards the longer slits as the wavelength increases. Furthermore, at longer wavelengths the effect of non-paraxiality on beam propagation becomes more evident. This is clearly manifested by the breaking of the main lobe along the propagation direction at wavelength 750 and 770 nm (Fig. 3(k), l) [23].
In the experiments, it can be seen that the electron yield starts decreasing significantly beyond 750 nm (Fig. 3(h), i). This can be explained by the fact that the photoemission in multiphoton PEEM is a nonlinear process and is very sensitive to the surface properties of the material. The work function of the evaporated gold used in our facility is about 5 eV. The incident photon energy is 1.65 eV at 750 nm. This clearly requires a 3 photon absorption process for the electrons to reach the vacuum state. The photoemission intensity at 3 photon process is proportional to the 6th power of the instantaneous total electric field integrated over time [24]. This nonlinear photoemission process results in low electron yields in the experiment. For even longer wavelengths, the photon energy becomes smaller and 4 photon absorption is required for photoemission. As a result, the electron counts in the photoemission becomes so small that we do only faintly observe the propagating plasmons for a wavelength of $\lambda =770$nm (Fig. 3(i)).
3. Generation bandwidth of Airy plasmons
Figure 3 clearly shows that the experimentally observable generation bandwidth is limited at the long wavelength side, which, as explained above, is due to the lower photon energy being insufficient to extract photoelectrons in the PEEM. This sets a limit to our experimental studies of ultra-short pulse excitation of Airy plasmons, which would require a broad bandwidth of the investigated Airy plasmon pulses. However the good correspondence of our experimental finding to the rigorous numerical simulations in Figure 3 justify to extrapolate the studies to a broader excitation bandwidth based on numerical simulation. Thus, we are now going to estimate numerically the bandwidth of the Airy plasmon’s efficient excitation. In addition, we also study the effect of the angle of incidence on the generation efficiency, which can be further utilized in pulsed Airy plasmons investigation in a later section. For this purpose, we calculated quantitatively the generation of Airy plasmons over the wavelength range of 600 nm to 1100 nm under normal and oblique ($4^{\circ }$) incidence. An overlap integral analysis is performed to quantify the Airy plasmon’s generation efficiency ($\eta$). It measures the overlap of our numerical solution with an analytically derived Airy wave profile. This is done using the following overlap integral formula
Analytical electric and magnetic fields are calculated by developing a simple formulation. We assume a TM polarized plane surface wave bound to the gold-air interface in the xz-plane at y=0 and decaying exponentially away from the interface along the the y-direction. The material above the interface (air) has the permittivity $\varepsilon _{I}$ and the material below the interface (gold) has the permittivity $\varepsilon _{II}$. The entire distribution of the TM-polarized surface plasmons in real space is then described as
Figure 4 depicts the resulting dependence of Airy plasmon’s generation efficiency on the wavelength of excitation under normal incidence (blue curve) and $4^{\circ }$ oblique incidence (red curve) calculated by the above model. For the given angles of incidence we investigated the effect of the grating scattering on the Airy plasmon formation. The scattered field from the grating can be assumed to be made of Airy plasmon modes and non-Airy plasmon modes. The idea here is to calculate the Airy plasmons mode’s contribution to the scattered fields. We did the analysis using forward propagated analytical fields [25–27]. For wavelengths from 600 to 1100 nm the overlap is calculated between the numerically simulated fields at a propagation distance of $z=10\,\mu$m and the analytically evolved fields. We find that for normal incidence the maximum generation efficiency is $\eta \approx 58\%$ at a wavelength of 790 nm. For oblique incidence the efficiency increases to $\eta \approx 64\%$ for a wavelength shifted to 820 nm.
4. Spatio-temporal evolution of Airy plasmon pulses
From the discussion above it is clear that the generation bandwidth of the Airy plasmon grating is large. This opens up the possibility to proceed with the investigation of short pulse excitation and to study the spatio-temporal dynamics of Airy plasmons and their evolution over time [2,4,28,29].
Now we will study the spatio-temporal propagation characteristics of Airy plasmon pulses generated at metal-dielectric interfaces by ultrashort pulse excitation with 20 fs (FWHM) Gaussian pulses having a bandwidth of $\approx 40$ nm (FWHM), which easily fits into the excitation bandwidth demonstrated in Fig. 4. These simulations are carried out using again the Lumerical FDTD solver. Since the efficient FDTD simulation of ultrashort pulses with broad spectral width and spatial inclination imposed numerical problems, we investigated the normal incidence case. Nevertheless, from the close correspondence of the results derived from the bandwidth estimation above for inclined and normal incidence we expect that the influence of the inclined excitation on the observed principle physics will be negligible also for short pulses. Thus we set the central wavelength of the 20 fs excitation pulses to 745 nm, which is the center of the plateau of high excitation efficiency for normal incidence in Fig. 4.
The time-integrated intensity of the simulated pulsed Airy SPP excitation and propagation is plotted in Fig. 5(a). It clearly shows that the Airy plasmon pulse preserves the properties of an Airy wave e.g. self accelerating bending of the beam and a non-diffracting nature of the main lobe. The full temporal evolution of such a spatio-temporal Airy wave packet is visualized in the supporting information.
We have taken several spatial location along the main lobe to study the temporal evolution of the Airy plasmon pulse. We also observe a high intensity spot along the main lobe at $z=13.41\,\mu$m along the propagation direction (shown as point 5 in Fig. 5(a)). Beyond this distance the side lobes vanish and the main lobe starts to diffract significantly.
Figure 5(b) shows the time traces of the intensity over time for the positions indicated in Fig. 5(a). The dotted curve is the source signal with a 50.24 fs offset. At position 1 ($0.5\,\mu$m from the grating edge), we see an enhancement in intensity over the source signal due to local resonances in the grating holes. This near-field effect vanishes at position 2 and the intensity peak diminishes. After this location, we see an increase in the peak intensity up to the maximum intensity position of the main lobe (position 5). The appearance of a maximum intensity along the propagation direction can be attributed to an initial phase mismatch of the grating period. After this maximum position the main lobe starts diffracting spatially, which is evident from the decrease in peak intensity at positions 6 and 7 in Fig. 5(a). We can also observe this intensity decreases in the transverse direction from point 7 to point 8, characterizing the divergence of the main lobe.
Furthermore, the spectral intensity distribution shown in the Fig. 5(c) does not change significantly. We have also calculated the time bandwidth product ($\Delta t\cdot \Delta {\nu }$), which becomes minimum at the maximum intensity position 5 (see Table 1).
5. Conclusion
In summary, we investigated the Airy plasmon propagation dynamics for a broad range of wavelengths using multiphoton PEEM and rigorous numerical simulations. The experimental and numerical results show good agreement. The photoelectron yield by PEEM clearly maps the constructive and destructive interference of the total electromagnetic field. We have also developed a semianalytical formulation to quantify the Airy like plasmon generation efficiency from the diffraction grating and it is estimated to be $\eta \,\approx 58\%$ under normal incidence and $\eta \,\approx 62\%$ for $4^{\circ }$ inclined incidence. We, furthermore demonstrated the creation of a spatio-temporal Airy plasmon pulse and estimated the time bandwidth product. The Airy plasmon pulse preserves its non-diffracting property along the propagation. These findings can be further utilized to engineer complexly shaped pulsed surface plasmon polaritons for applications in e.g. nanoparticle sensing or manipulation as well as in high speed on-chip optical signal processing.
Funding
German Academic Exchange Service New Delhi; Deutsche Forschungsgemeinschaft (278747906, SFB 1375).
Acknowledgments
The authors thank Matthias Zilk for his contribution to analytical modeling and numerical analysis in this manuscript.
Disclosures
The authors declare no conflicts of interest.
References
1. M. V. Berry and N. L. Balazs, “Nonspreading wave packets,” Am. J. Phys. 47(3), 264–267 (1979). [CrossRef]
2. G. A. Siviloglou and D. N. Christodoulides, “Accelerating finite energy airy beams,” Opt. Lett. 32(8), 979–981 (2007). [CrossRef]
3. G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Observation of accelerating airy beams,” Phys. Rev. Lett. 99(21), 213901 (2007). [CrossRef]
4. A. Chong, W. H. Renninger, D. N. Christodoulides, and F. W. Wise, “Airy-bessel wave packets as versatile linear light bullets,” Nat. Photonics 4(2), 103–106 (2010). [CrossRef]
5. K. Y. Kim, C. Y. Hwang, and B. Lee, “Slow non-dispersing wavepackets,” Opt. Express 19(3), 2286–2293 (2011). [CrossRef]
6. O. Vallée, Airy Functions and Applications to Physics (Imperial College Press, Chichester, United Kingdom, 2004).
7. M. S. N. K. Efremidis, Z. Chen, and D. N. Christodoulides, “Airy beams and accelerating waves: an overview of recent advances,” Optica 6(5), 686–701 (2019). [CrossRef]
8. M. Mazilu, D. J. Stevenson, F. Gunn-Moore, and K. Dholakia, “Light beats the spread: "non-diffracting" beams,” Laser Photonics Rev. 4(4), 529–547 (2010). [CrossRef]
9. S. A. Maier, Plasmonics: Fundamentals and Applications (Springer, Berlin, 2007).
10. A. Salandrino and D. N. Christodoulides, “Airy plasmon: a nondiffracting surface wave,” Opt. Lett. 35(12), 2082–2084 (2010). [CrossRef]
11. P. Zhang, S. Wang, Y. Liu, X. Yin, C. Lu, Z. Chen, and X. Zhang, “Plasmonic airy beams with dynamically controlled trajectories,” Opt. Lett. 36(16), 3191–3193 (2011). [CrossRef]
12. A. Minovich, A. E. Klein, N. Janunts, T. Pertsch, D. N. Neshev, and Y. S. Kivshar, “Generation and near-field imaging of airy surface plasmons,” Phys. Rev. Lett. 107(11), 116802 (2011). [CrossRef]
13. A. Minovich, A. E. Klein, D. N. Neshev, T. Pertsch, Y. S. Kivshar, and D. N. Christodoulides, “Airy plasmons: non-diffracting optical surface waves,” Laser Photonics Rev. 8(2), 221–232 (2014). [CrossRef]
14. O. Schmidt, M. Bauer, C. Wiemann, R. Porath, M. Scharte, O. Andreyev, G. Schönhense, and M. Aeschlimann, “Time-resolved two photon photoemission electron microscopy,” Appl. Phys. B 74(3), 223–227 (2002). [CrossRef]
15. P. Vasa, C. Ropers, R. Pomraenke, and C. Lienau, “Ultra-fast nano-optics,” Laser Photonics Rev. 3(6), 483–507 (2009). [CrossRef]
16. K. F. MacDonald, Z. L. Sámson, M. I. Stockman, and N. I. Zheludev, “Ultra-fast nano-optics,” Nat. Photonics 3(1), 55–58 (2009). [CrossRef]
17. A. N. Grigorenko, M. Polini, and K. S. Novoselov, “Graphene plasmonics,” Nat. Photonics 6(11), 749–758 (2012). [CrossRef]
18. F. Xia, H. Wang, D. X. M. Dubey, and A. Ramasubramaniam, “Two-dimensional material nanophotonics,” Nat. Photonics 8(12), 899–907 (2014). [CrossRef]
19. C. Guan, T. Yuan, R. Chu, Y. Shen, Z. Zhu, J. Shi, P. Li, L. Yuan, and G. Brambilla, “Generation of ultra-wideband achromatic airy plasmons on a graphene surface,” Opt. Lett. 42(3), 563–566 (2017). [CrossRef]
20. A. Taflove, S. C. Hagness, and A. House, Computational Electrodynamics: The Finite-Difference Time-Domain Method (Artech House, Boston, 2005).
21. https://www.lumerical.com/products/fdtd-solutions/.
22. P. B. Johnson and R. W. Christy, “Optical constants of noble metals,” Phys. Rev. B 6(12), 4370–4379 (1972). [CrossRef]
23. A. V. Novitsky and D. V. Novitsky, “Nonparaxial airy beams: role of evanescent waves,” Opt. Lett. 34(21), 3430–3432 (2009). [CrossRef]
24. T. Kaiser, M. Falkner, J. Qi, A. Klein, M. Steinert, C. Menzel, C. Rockstuhl, and T. Pertsch, “Characterization of a circular optical nanoantenna by nonlinear photoemission electron microscopy,” Appl. Phys. B 122(3), 53 (2016). [CrossRef]
25. H. J. W. M. Hoekstra, “On beam propagation methods for modelling in integrated optics,” Opt. Quantum Electron. 29(2), 157–171 (1997). [CrossRef]
26. R. Scarmozzino and R. M. Osgood, “Comparison of finite difference and fourier transform solutions of the parabolic wave equation with emphasis on integrated optics applications,” J. Opt. Soc. Am. A 8(5), 724 (1991). [CrossRef]
27. E. O. Brigham, The Fast Fourier Transformation (Prentice Hall, Englewood Cliffs, NJ, 1974).
28. C. Lemke, C. Schneider, T. Leißner, D. Bayer, J. W. Radke, P. M. A. Fischer, A. B. Evlyukhin, B. N. Chichkov, C. Reinhardt, M. Bauer, and M. Aeschlimann, “Spatiotemporal characterization of spp pulse propagation in two-dimensional plasmonic focusing devices,” Nano Lett. 13(3), 1053–1058 (2013). [CrossRef]
29. B. E. A. Saleh and M. C. Teich, Fundamentals of Photonics (John Wiley & Sons, New Jersey, 2007).