Self-focusing of femtosecond surface plasmon polaritons

We study the propagation of femtosecond pulses in nonlinear metal-dielectric plasmonic waveguiding structures by employing the finitedifference time-domain numerical method. Self-focusing of plasmon pulses is observed for defocusing Kerr-like nonlinearity of the dielectric medium due to normal dispersion. We compare the nonlinear propagation of plasmon pulses along a single metal-dielectric interface with the propagation within a metal-dielectric-metal slot waveguide and observe that nonlinear effects are more pronounced for the single surface where longer propagation length may compensate for lower field confinement. © 2013 Optical Society of America OCIS codes: (190.4350) Nonlinear optics at surfaces; (190.5530) Pulse propagation and temporal solitons; (240.6680) Surface plasmons. References and links 1. Plasmonic Nanoguides and Circuits, Ed. S.I. Bozhevolnyi (Pan Stanford Publ, Singapore 2009). 2. A.G. Litvak and V.A. Mironov, “Surface waves on the separation boundary between nonlinear media,” Radiophys. Quant. Electron. 11, 1096-1101 (1968). 3. V.M. Agranovich, V.S. Babichenko, and V.Ya. Chernyak, “Nonlinear surface polaritons,” Sov. Phys. JETP Lett. 32, 512-515 (1980) [In Russian: Pisma JETF 32, 532-535 (1980)]. 4. G.I. Stegeman, C.T. Seaton, J. Ariyasu, R.F. Wallis, and A.A. Maradudin, “Nonlinear electromagnetic waves guided by a single interface,” J. Appl. Phys. 58, 2453-2459 (1985). 5. A.D. Boardman, A.A. Maradudin, G.I. Stegeman, T. Twardowski, and E.M. Wright, “Exact theory of nonlinear p-polarized optical waves,” Phys. Rev. A 35, 1159-1164 (1987). 6. J. Ariyasu, C.T. Seaton, G.I. Stegeman, A.A. Maradudin, and R.F. Wallis, “Nonlinear surface polaritons guided by metal films,” J. Appl. Phys. 58, 2460-2466 (1985). 7. A.D. Boardman, G.S. Cooper, A.A. Maradudin, and T.P. Shen, “Surface-polariton solitons,” Phys. Rev. B 34, 8273-8278 (1986). 8. E. Feigenbaum and M. Orenstein, “Plasmon-soliton,” Opt. Lett. 32, 674-676 (2007). 9. A.R. Davoyan, I.V. Shadrivov, and Yu.S. Kivshar, “Self-focusing and spatial plasmon-polariton solitons,” Opt. Express17, 21732-21737 (2009). 10. A. Marini, D. Skryabin, and B. Malomed, “Stable spatial plasmon solitons in a dielectric-metal-dielectric geometry with gain and loss,” Opt. Express 19, 6616-6622 (2011). 11. Yu.S. Kivshar and G.P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (Academic, San Diego, 2003), p. 540. 12. S.O. Demokritov, B. Hillebrands, and A.N. Slavin, “Brillouin light scattering studies of confined spin waves: linear and nonlinear confinment,” Phys. Rep. 348, 441-489 (2001). 13. O. Büttner, M. Bauer, S. O. Demokritov, B. Hillebrands, Yuri S. Kivshar, V. Grimalsky, Yu. Rapoport, and A. N. Slavin, “Linear and nonlinear diffraction of dipolar spin waves in yttrium iron garnet films observed by spaceand time-resolved Brillouin light scattering,” Phys. Rev. B 61, 11576-11587 (2000). 14. P.B. Johnson and R.W. Christy, “Optical constants of the noble metals,” Phys. Rev. B 6, 4370-4379 (1972). #175981 $15.00 USD Received 11 Sep 2012; revised 12 Oct 2012; accepted 13 Oct 2012; published 10 Jan 2013 (C) 2013 OSA 14 January 2013 / Vol. 21, No. 1 / OPTICS EXPRESS 1121 15. A. Taflove and S.C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method (Artech House, Boston 2005). 16. I. S. Maksymov, A. A. Sukhorukov, A. V. Lavrinenko, and Yu. S. Kivshar, “Comparative study of FDTD-adopted numerical algorithms for Kerr nonlinearities,” IEEE Antennas Wireless Propag. Lett. 10, 143-146 (2011). 17. Z.L. Sámson, P. Horak, K.F. MacDonald, and N.I. Zheludev, “Femtosecond surface plasmon pulse propagation,” Opt. Lett.36, 250-252 (2011). 18. M. S. Petrovic, M. R. Belic, C. Denz, and Yu. S. Kivshar, “Counterpropagating optical beams and solitons,” Laser Photonics Rev. 5, 214-233 (2010). 19. A. Pusch, J. M. Hamm, and O. Hess, “Controllable interaction of counterpropagating solitons in three-level media,” Phys. Rev. A82, 023805 (2010). 20. A. R. Davoyan, I. V. Shadrivov, A. A. Zharov, D. K. Gramotnev, and Yu. S. Kivshar, “Nonlinear nanofocusing in tapered plasmonic waveguides,” Phys. Rev. Lett. 105, 116804 (2010).


Introduction
The study of pulse and beam propagation in plasmonic waveguides is of great interest for telecommunication, computing, and information processing [1].It is expected that nonlinear effects may even further enhance the capabilities of plasmonic devices also opening the door to the observation of novel effects such as plasmon self-focusing, self-trapping, and frequency conversion.
Early pioneering studies of nonlinear propagation of surface plasmon polaritons employed simple approximations for the analysis of Kerr nonlinearities in a single metal-dielectric interface geometry [2][3][4][5], and in a metal film embedded into a nonlinear Kerr medium [6].However, only very few works have considered so far nonlinear guided waves in metal-dielectric structures, reporting on the studies of temporal effects and the formation of temporal surfacepolariton solitons [7].More recently, several groups have analyzed spatial localization of plasmon beams and the formation of spatial plasmon solitons [8][9][10].
One of the main challenges associated with the study of nonlinear effects in plasmonic waveguides, such as self-focusing and the formation of optical spatial and temporal solitons [11], is the presence of strong characteristic losses in metals at optical frequencies.Indeed, such losses can even be much stronger than effects of nonlinearity and dispersion, and they may consequently overshadow any manifestation of nonlinear self-action effects.Notwithstanding, the study of strong self-focusing of nonlinear waves in lossy media is well documented in other fields such as magnetic spin waves [12].In particular, it is well established that nonlinear evolution of magnetostatic spin waves in thin YIG magnetic films may demonstrate strong spatiotemporal self-focusing and the formation of localized waves in the form of spin-wave solitons and bullets [13]; the effects that can be readily observed even in the presence of relatively strong losses.Moreover, self-action manifests itself as a transient effect at the initial period of the pulse evolution when the attenuation is weaker than the nonlinearity.
In this paper, we study the propagation of femtosecond pulses in nonlinear metal-dielectric plasmonic waveguiding structures by employing the finite-difference time-domain numerical method.We demonstrate that in the presence of a defocusing Kerr-like nonlinear response of the dielectric medium, femtosecond plasmonic pulses undergo self-focusing similar to the pulses in optical fibres [11] but with a transient formation of temporal plasmon-solitons, and then experience broadening due to strong losses in metal.We consider two major types of plasmonic waveguides with a nonlinear dielectric, namely a single metal-dielectric interface and a metaldielectric-metal (MDM) slot waveguide.We find that nonlinear effects are more pronounced for plasmon pulses propagating along a single metal-dielectric interface where longer propagation length may compensate for lower field confinement.Finally, we study the interaction of the plasmonic pulses with the counter-propagating continuous wave (CW) plasmons and show that, by changing the amplitude of incoming CW, one can control the phase of the pulses.

Methods
To calculate the dynamics of the propagating TM polarized plasmon pulses we use a finitedifference time-domain (FDTD) numerical technique.For TM polarization, Maxwell's equations reduce to the following form, To model the metallic response of silver in the telecommunication regime we employ the Drude model with ω p = 1.39 • 10 16 rad/s and Γ = 3.2 • 10 13 rad/s, which fits the experimental data of Johnson and Christy [14].The Kerr-type response of the nonlinear dielectric is calculated [15,16] via a recursive application of the nonlinear material relation x + E 2 y , and we also take ε d = 4.9.In our analysis we consider two structures.The first structure is a simple metal-dielectric interface [see Fig. 1(a)] supporting the propagation of a surface plasmon with the magnetic field distribution of the form, The complex propagation constant β can be found from the boundary conditions as where k m,d are transverse wavenumbers: , and k 0 = ω/c.The propagation length of the surface plasmon ([2Im(β )] −1 ) at telecommunication wavelength (λ 0 = 1550 nm) is evaluated to 105 µm, and the corresponding imaginary part of the wavenumber is 4764 m −1 .
The second structure investigated here [see Fig. 1(b)] is a metal-dielectric-metal (MDM) slot waveguide.Here the permittivities are given as where we use the thickness of the dielectric layer a = 100 nm.The dispersion equation for the propagation constant is then given as We calculate the propagation length for the slot waveguide as 20 µm, substantially shorter than the propagation length of the surface plasmon at a single interface.Both structures show normal dispersion at telecommunication wavelengths, so that a negative ε NL is necessary to achieve the effective pulse compression [11].
To inject the correct pulse form into the FDTD computation domain, the profile of the linear mode is first calculated and then used as input for a total-field/scattered-field plane [15], in which both electric and magnetic fields are staggered with a phase difference corresponding to the propagation constant of the mode.For the study of the propagation of a single surface plasmon, we use a grid cell size of ∆x = 20 nm and a total length of the computation domain along x-axis of L = 400 µm.The dielectric extends for 3 µm in the y-direction to minimize the effects of the boundaries, whereas the metal extends for 1 µm.
In the calculation of the MDM slot waveguide the resolution is increased to ∆x = 10 nm in order to sufficiently resolve the large field gradients across the structure.As a consequence of that, the computational domain has to be smaller with L = 220 µm, which corresponds to more than 10 times the propagation length, and thus it is long enough to capture all nonlinear effects.In both the cases perfectly matched layer boundary conditions are used at the boundaries of the simulation domains.Using the above resolutions in a two-dimensional FDTD simulation results in some numerical dispersion that is, however, still much smaller than the dispersion introduced by the Drude model of free electrons in the metal and the wave-guiding geometry.Note that throughout this work, the intensities are normalized to ε NL , and we monitored the maximum nonlinear correction to the dielectric permittivity δ ε max that is reached in the nonlinear medium.We use this parametrization since the nonlinear effects that can be achieved in plasmonic systems are ultimately limited by the saturation of the nonlinearities.

Results
Figures 2(a) and 2(b) show schematically the propagation dynamics of nonlinear and linear pulses in lossy metal-dielectric plasmonic structures studied numerically.In the nonlinear case, the pulses are first compressed (or self-focused) until dissipation has reduced the pulse amplitude into the linear regime and dispersive pulse broadening takes over [see Fig. 2 (a)].In the linear regime pulses experience continuous broadening [see Fig. 2(b)], as also observed in experiment [17].
In Figs.
We can see that the pulses are compressed more quickly in the slot waveguide because of the higher field confinement.The high losses (short propagation length) in the slot waveguide, however, lead to the plasmon pulse quickly reaching the linear regime, where it starts to broaden again.Thus the surface plasmon can, when excited with the same maximal field amplitude, experience substantially stronger compression, as can be seen in Fig. 3(a).Here, the compression only slows down considerably after about 200 µm (twice the propagation length).
One of the fascinating properties of optical solitons is their ability to interact during collisions [18].The inherent weakness of the Kerr nonlinearity, however, limits the effective interaction strength for short pulses, unlike, e.g. , in self-induced transparency solitons where even femtosecond pulses can interact very strongly and induce soliton birth or breakup [19].Therefore, we now analyse the effect that a counter-propagating quasi-CW excitation has on a short femtosecond pulse.
Figure 4 shows the maximum phase shift that can be accumulated by a pulse subjected to the interaction with a counter-propagating CW excitation on the surface (black triangles) or the slot waveguide (red diamonds).The achievable phase shift depends mainly on the maximum change in the permittivity before nonlinearity saturation sets in.A CW plasmon entering from the right [see Fig. 4(a) for a schematic view] attenuates upon propagation and ultimately it is too weak to induce a phase shift.Again, it is the shorter propagation lengths in the slot waveguide that limit the nonlinear effect.The phase shift is thus larger for the surface plasmon on a single surface.Here, the same (but sign-reversed) effect could be achieved with a positive values of the third order susceptibility.Such a nonlinear interaction setup can be employed for controlling the phase of the plasmonic pulses.We note that for the plasmon propagating distances and characteristic times that it takes for the plasmons to propagate these distances, we didn't observe any onsets of modulational instability of the CW.

Fig. 1 .
Fig. 1.Schematic of the pulse propagation along (a) metal-dielectric interface and in (b) metal-dielectric-metal slot waveguide.(c,d) Pulse profiles shown as squared magnetic field carried by the pulses.

Fig. 2 .
Fig.2.Schematic of pulse propagation in (a) nonlinear and (b) linear case.The pulse energy always decreases during propagation.In the nonlinear case (a), the pulse width w first decreases before a linear propagation regime is reached, and then it increases again.In the linear regime (b), a steady increase of the pulse duration is observed.

Fig. 3 .
Fig. 3. Evolution of the pulse widths for Gaussian pulses with propagation distance for (a) single interface for the nonlinear (red triangles) and linear (black circles) dielectrics, and (b) in the slot waveguide with nonlinear (red triangles) and linear (black circles) dielectrics in the slot.The maximum nonlinear change of the dielectric permittivity at x = 0 is found as δ ε max = −ε NL E 2 0 = 0.03.
Figures 2(a) and 2(b) show schematically the propagation dynamics of nonlinear and linear pulses in lossy metal-dielectric plasmonic structures studied numerically.In the nonlinear case, the pulses are first compressed (or self-focused) until dissipation has reduced the pulse amplitude into the linear regime and dispersive pulse broadening takes over [see Fig.2 (a)].In the linear regime pulses experience continuous broadening [see Fig.2(b)], as also observed in experiment[17].In Figs.3(a) and 3(b) the qualitative behaviour indicated in Figs.2(a) and 2(b) is quantified for both the single interface and slot waveguide.The initial Gaussian pulses are given byE(t) = E 0 cos [ω(t − t 0 )] exp[−(t − t 0 ) 2 /τ 2 ] with τ = 36 fs.We can see that the pulses are compressed more quickly in the slot waveguide because of the higher field confinement.The high losses (short propagation length) in the slot waveguide, however, lead to the plasmon pulse quickly reaching the linear regime, where it starts to broaden again.Thus the surface plasmon can, when excited with the same maximal field amplitude, experience substantially stronger compression, as can be seen in Fig.3(a).Here, the compression only slows down considerably after about 200 µm (twice the propagation length).One of the fascinating properties of optical solitons is their ability to interact during collisions[18].The inherent weakness of the Kerr nonlinearity, however, limits the effective interaction strength for short pulses, unlike, e.g. , in self-induced transparency solitons where even femtosecond pulses can interact very strongly and induce soliton birth or breakup[19].Therefore, we now analyse the effect that a counter-propagating quasi-CW excitation has on a short femtosecond pulse.Figure4shows the maximum phase shift that can be accumulated by a pulse subjected to the interaction with a counter-propagating CW excitation on the surface (black triangles) or the slot waveguide (red diamonds).The achievable phase shift depends mainly on the maximum

Fig. 4 .
Fig. 4. (a) Schematic of the nonlinear soliton collision.A pulse coming from the left (black line) collides with a quasi-CW excitation coming from the right (light blue line).Note that the CW excitation loses energy during propagation.(b) Maximum phase-shifts attainable for a pulse colliding with a CW excitation evaluated for various initial nonlinear permittivity changes δ ε max due to the CW field.Black triangles denote the data for a single interface, and red diamonds -the data for the MDM slot waveguide.