Controllable multiple plasmonic bending beams via polarization of incident waves

Plasmonic bending beams, which preserve their spatial shapes while propagating along curved trajectories in metal-dielectric interface, offer important applications in the fields of fiber sensor, optical trapping, and micro-nano manipulation. In this work, circular hole array, as a local point-like sources of surface plasmon polaritons, is designed on the metal film to generate multiple plasmonic bending beams. The electric field intensity of multiple plasmonic bending beams is controlled by polarization angle of input light. In addition, the electric filed intensity of multiple plasmonic bending beams relies on circle hole radius. These findings provide guidance in the design and optimization of plasmonic bending beam generators. © 2017 Optical Society of America under the terms of the OSA Open Access Publishing Agreement OCIS codes: (260.2110) Electromagnetic optics; (240.6680) Surface plasmons; (060.5060) Phase modulation; (250.5403) Plasmonics; (240.5440) Polarization-selective devices. References and links 1. C. Genet and T. W. Ebbesen, “Light in tiny holes,” Nature 445(7123), 39–46 (2007). 2. W. L. Barnes, A. Dereux, and T. W. Ebbesen, “Surface plasmon subwavelength optics,” Nature 424(6950), 824– 830 (2003). 3. M. L. Brongersma and V. M. Shalaev, “The case for plasmonics,” Science 328(5977), 440–441 (2010). 4. Y. Liu, H. Xu, F. Stief, N. Zhitenev, and M. Yu, “Far-field superfocusing with an optical fiber based surface plasmonic lens made of nanoscale concentric annular slits,” Opt. Express 19(21), 20233–20243 (2011). 5. M. L. Juan, M. Righini, and R. Quidant, “Plasmon nano-optical tweezers,” Nat. Photonics 5(6), 349 (2011). 6. C. Min, Z. Shen, J. Shen, Y. Zhang, H. Fang, G. Yuan, L. Du, S. Zhu, T. Lei, X. Yuan, and X. C. Yuan, “Focused plasmonic trapping of metallic particles,” Nat. Commun. 4, 2891 (2013). 7. X. Zhu, A. Schülzgen, H. Wei, K. Kieu, and N. Peyghambarian, “White light Bessel-like beams generated by miniature all-fiber device,” Opt. Express 19(12), 11365–11374 (2011). 8. J. Baumgartl, M. Mazilu, and K. Dholakia, “Optically mediated particle clearing using Airy wavepackets,” Nat. Photonics 2(11), 675–678 (2008). 9. J. Arlt, V. Garces-Chavez, W. Sibbett, and K. Dholakia, “Optical micromanipulation using a Bessel light beam,” Opt. Commun. 197(4), 239–245 (2001). 10. 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). 11. A. Salandrino and D. N. Christodoulides, “Airy plasmon: a nondiffracting surface wave,” Opt. Lett. 35(12), 2082–2084 (2010). 12. A. E. Klein, A. Minovich, M. Steinert, N. Janunts, A. Tünnermann, D. N. Neshev, Y. S. Kivshar, and T. Pertsch, “Controlling plasmonic hot spots by interfering Airy beams,” Opt. Lett. 37(16), 3402–3404 (2012). 13. W. Liu, D. N. Neshev, I. V. Shadrivov, A. E. Miroshnichenko, and Y. S. Kivshar, “Plasmonic Airy beam manipulation in linear optical potentials,” Opt. Lett. 36(7), 1164–1166 (2011). 14. L. Li, T. Li, S. M. Wang, C. Zhang, and S. N. Zhu, “Plasmonic Airy beam generated by in-plane diffraction,” Phys. Rev. Lett. 107(12), 126804 (2011). 15. C. Guan, M. Ding, J. Shi, P. Hua, P. Wang, L. Yuan, and G. Brambilla, “Experimental observation and analysis of all-fiber plasmonic double Airy beams,” Opt. Express 22(15), 18365–18371 (2014). 16. L. Li, T. Li, S. M. Wang, and S. N. Zhu, “Collimated plasmon beam: nondiffracting versus linearly focused,” Phys. Rev. Lett. 110(4), 046807 (2013). 17. L. Li, T. Li, S. M. Wang, and S. N. Zhu, “Steering plasmon beam from a point source,” Opt. Lett. 37(24), 5091– 5093 (2012). 18. A. Arie, I. Epstein, and Y. Lilach, “Shaping plasmonic light beams with near-field plasmonic holograms,” J. Opt. Soc. Am. B 31(7), 1642–1647 (2015). 19. I. Epstein and A. Arie, “Arbitrary bending plasmonic light waves,” Phys. Rev. Lett. 112(2), 023903 (2014). Vol. 25, No. 24 | 27 Nov 2017 | OPTICS EXPRESS 29659 #307596 https://doi.org/10.1364/OE.25.029659 Journal © 2017 Received 20 Sep 2017; revised 6 Nov 2017; accepted 6 Nov 2017; published 13 Nov 2017 20. I. Epstein and A. Arie, “Dynamic generation of plasmonic bottle-beams with controlled shape,” Opt. Lett. 39(11), 3165–3168 (2014). 21. A. Libster-Hershko, I. Epstein, and A. Arie, “Rapidly accelerating Mathieu and Weber surface plasmon beams,” Phys. Rev. Lett. 113(12), 123902 (2014). 22. C. Sheng, R. Bekenstein, H. Liu, S. Zhu, and M. Segev, “Wavefront shaping through emulated curved space in waveguide settings,” Nat. Commun. 7, 10747 (2016). 23. C. Sheng, H. Liu, Y. Wang, S. N. Zhu, and D. A. Genov, “Trapping light by mimicking gravitational lensing,” Nat. Photonics 7(11), 902–906 (2013). 24. J. Lin, J. Dellinger, P. Genevet, B. Cluzel, F. D. Fornel, and F. Capasso, “Cosine-Gauss plasmon beam: a localized long-range nondiffracting surface wave,” Phys. Rev. Lett. 109(9), 093904 (2012). 25. S. Wei, J. Lin, Q. Wang, G. Yuan, L. Du, R. Wang, L. Xu, M. Hong, C. Min, and X. Yuan, “Singular diffraction-free surface plasmon beams generated by overlapping phase-shifted sources,” Opt. Lett. 38(7), 1182– 1184 (2013). 26. P. Genevet, J. Dellinger, R. Blanchard, A. She, M. Petit, B. Cluzel, M. A. Kats, F. de Fornel, and F. Capasso, “Generation of two-dimensional plasmonic bottle beams,” Opt. Express 21(8), 10295–10300 (2013). 27. K. Xiao, S. Wei, C. Min, G. Yuan, S. W. Zhu, T. Lei, and X. C. Yuan, “Dynamic cosine-Gauss plasmonic beam through phase control,” Opt. Express 22(11), 13541–13546 (2014). 28. I. Epstein, R. Remez, Y. Tsur, and A. Arie, “Generation of intensity-controlled two-dimensional shapepreserving beams in plasmonic lossy media,” Optica 3(1), 15–19 (2016). 29. J. Lin, J. P. B. Mueller, Q. Wang, G. Yuan, N. Antoniou, X. C. Yuan, and F. Capasso, “Polarization-controlled tunable directional coupling of surface plasmon polaritons,” Science 340(6130), 331–334 (2013). 30. L. Cao, W. Li, Y. K. Wang, X. J. Tian, G. A. Li, H. R. Zheng, and Z. Y. Zhang, “Manipulating surface plasmon polaritons using F-shaped nanoslits array,” IEEE Photonics Technol. Lett. 26(12), 1247–1250 (2014). 31. O. You, B. Bai, X. Wu, Z. Zhu, and Q. Wang, “A simple method for generating unidirectional surface plasmon polariton beams with arbitrary profiles,” Opt. Lett. 40(23), 5486–5489 (2015). 32. E. Y. Song, S. Y. Lee, J. Hong, K. Lee, Y. Lee, G. Y. Lee, H. Kim, and B. Lee, “A double-lined metasurface for plasmonic complex-field generation,” Laser Photonics Rev. 10(2), 299–306 (2016). 33. G. Y. Lee, S. Y. Lee, H. Yun, H. Park, J. Kim, K. Lee, and B. Lee, “Near-field focus steering along arbitrary trajectory via multi-lined distributed nanoslits,” Sci. Rep. 6, 33317 (2016). 34. S. Y. Lee, K. Kim, S. J. Kim, H. Park, K. Lee, and B. Lee, “Plasmonic meta-slit: shaping and controlling nearfield focus,” Optica 2(1), 6 (2015). 35. O. Avayu, I. Epstein, E. Eizner, and T. Ellenbogen, “Polarization controlled coupling and shaping of surface plasmon polaritons by nanoantenna arrays,” Opt. Lett. 40(7), 1520–1523 (2015). 36. S. Xiao, F. Zhong, H. Liu, S. Zhu, and J. Li, “Flexible coherent control of plasmonic spin-Hall effect,” Nat. Commun. 6, 8360 (2015). 37. Q. Xu, X. Zhang, Y. Xu, Ch. Ouyang, Zh. Tian, J. Gu, J. Li, Sh. Zhang, J. Han, and W. Zhang, “Polarization‐controlled surface plasmon holography,” Laser Photonics Rev. 11(1), 1600212 (2017). 38. Y. G. Chen, Y. H. Chen, and Z. Y. Li, “Direct method to control surface plasmon polaritons on metal surfaces,” Opt. Lett. 39(2), 339–342 (2014). 39. H. Li, Y. Z. Xu, G. Wang, T. Fu, L. Wang, and Z. Y. Zhang, “Converting surface plasmon polaritons into spatial bending beams through graded dielectric rectangles over metal film,” Opt. Commun. 383, 423–429 (2017). 40. L. Froehly, F. Courvoisier, A. Mathis, M. Jacquot, L. Furfaro, R. Giust, P. A. Lacourt, and J. M. Dudley, “Arbitrary accelerating micron-scale caustic beams in two and three dimensions,” Opt. Express 19(17), 16455– 16465 (2011). 41. H. Li, Y. Z. Xu, Y. Qu, T. Fu, B. Zhang, and Z. Y. Zhang, “Enhancement of focusing properties by interfering spatial bending beams,” Europhys. Lett. 116(5), 55002 (2017). 42. M. A. Ordal, L. L. Long, R. J. Bell, S. E. Bell, R. R. Bell, R. W. Alexander, Jr., and C. A. Ward, “Optical properties of the metals al, co, cu, au, fe, pb, ni, pd, pt, ag, ti, and w in the infrared and far infrared,” Appl. Opt. 22(7), 1099 (1983). 43. Z. J. Zhang, R. W. Peng, Z. Wang, F. Gao, X. R. Huang, W. H. Sun, Q. J. Wang, and M. Wang, “Plasmonic antenna array at optical frequency made by nanoapertures,” Appl. Phys. Lett. 93(17), 171110 (2008).


Introduction
Surface plasmon polaritons (SPPs) are electromagnetic waves in two-dimensional system, that consist of collective electron oscillations in metal-dielectric interface decay, thus, their diffraction limits can be surpassed [1,2], thereby offering important applications in the fields of photonics and electronics [3][4][5].Plasmonic bending beams, as a specific SPP waves, which preserve their shapes while propagating along desired curved trajectories in metal-dielectric interface, are attracting increased attention.They offer potential applications in the fields of optical trapping [4][5][6] and micro-nano manipulation [7][8][9].
Controlling plasmonic complex fields by polarization-sensitive nanostructure arrays are important for enriching the components of nanophotonic devices to analyze the polarization [29,30].Various polarization-sensitive nanostructure arrays, such as aperture arrays [29], Fshaped nanoslits [30], Δ-shaped nanoantennas [31], and double-lined distribution nanoslits arrays [32,33], were designed to manipulate propagation directional of SPP waves.Multilined distributed nanoslits arrays were designed to achieve focused SPPs [32][33][34].Crossshaped nanoantennas were designed to generate double focused SPPs [35].Multiple rings of nanoslits with special orientations was designed to control SPP orbitals [36].Ring-shaped hologram consisting of slit-pair array was designed to achieve switchable SPPs [37].All of these reports demonstrate that phase and amplitude of the generation SPP beams can be independently controlled by polarization-sensitive nanostructures.These reports also show that these tunable and compact beams have extensive potential applications [32].However, these beams were mostly controlled by adjusting the structural orientation angles.If users want to control structural orientation angles, they need rebuild the structures, which is not convenient.
In this study, a circular hole array is used as local point-like sources of SPPs to generate SPP wavefront.The desired multiple plasmonic bending beams are generated with propagation along the desired trajectory.The intensity of the electric fields of multiple plasmonic bending beams is controlled by the polarization angle of the input light.In addition, the dependence of the intensity of the electric fields of multiple plasmonic bending beams on the radius of circular hole is studied.All these findings can provide deep applications in designing polarization-controlled plasmonic bending beam generators or polarization sensing devices.

Structure and calculation method
All of motion of electro-magnetic can be described by the Helmholtz equation [19].Plasmonic bending beams, the solutions of the two-dimensional (propagation in the x-y plane) Helmholtz equation, can be expressed [19,38]  ( , ) ( , ) 0, where k sp represents the wave-vector of SPPs on the metal surface.E SP (x, y) represents the distribution of SPP waves in x-y plane.Figure 1 shows the geometric diagram of curve arrangement of circular hole array.Every circular hole is regarded as an independent electric field radiation source.Both the amplitude and phase of SPP are controlled by arrangement of circular hole array.In order to achieve the distribution of plasmonic bending beams, the phase modulation method is introduced.Phase φ(x) of the required plasmonic bending beams was calculated utilizing Legendre transform method [18][19][20].The geometric relationship between the phase, the angle θ, and the desired arbitrary bending trajectory f(y) is defined in Fig. 1.
The phase of the required plasmonic bending beams can be obtained for the case of nonparaxial ( where tan is the first-order derivative of f (y).In this study, in order to deduce simply, the case of paraxial sin θ ≈tan θ was applied.The phase of the required plasmonic bending beams was obtained for paraxial regimes [20,36,[38][39][40] The position of every radiation source was directly located by the phase ( ) x φ .In this study, to avoid near-field interaction between adjacent neighboring holes which can reduce the intensity of bending beams, the distance W is set as 0.6 μm.The quadratic curve (f(y) = -ay 2 ) is chosen, where constant a is 1.13 × 10 −2 .According to Eq. ( 3), the desired phase ( ) is obtained for paraxial regimes [41].The distances

Results and discussion
The electric field intensity distribution in the x-y plane is depicted in Fig. 3

4(c) -4(e)], the electric field intensities of beams I and II increase, and the electric field intensity of beam III decreases.
To understand the relationship between the electric field intensity of SPP waves and the polarization angle β, the radiation properties of a circular hole array are analyzed theoretically.When every circular hole of an array is illuminated by incident light, the surface plasmons are excited and can be regarded as an electric dipole.The electric dipoles can be written as [43]: The first term represents the interaction between in-plane incident electric field and the circular holes.n δ depends on the radius r of the circular hole and the permittivity of a metal ε m .The second term represents the interaction among the circular holes in the array.The third term represents the field from the propagating SPPs [43].In this study, the distance between neighboring circular holes is larger than 0.2 μm, and the interaction among other dipoles (the second and third terms) is relatively weak.Every electric dipole in the array is approximated as P n = δ n E 0 .Thus, the electric distribution of the n-electric dipoles can be approximatively written as .
The far-field distribution of SPP waves E sp (x, y) from each independent radiation source (electric dipoles) can be obtained as a superposition of the electric fields that are radiated from every source, where n is an integer, ( ) ...
The solid curves in Fig. 4(f) result in the fitting of Eq. ( 6).The goodness for beams I, II, and III is 0.969, 0.963, and 0.977, respectively.These results demonstrate that the electric field intensity of multiple plasmonic bending beams (I, II, and III) are controlled by the polarization angle of incident light, which provides a method of dynamic controlling and arbitrarily selecting multiple plasmonic bending beams.The effects of the radius of a circular hole on the multiple plasmonic bending beams are also investigated.6).These results demonstrate that the value of electric field intensity can be controlled by simply selecting suitable radius circular whole arrays.The third plasmonic bending beams (beam III) is smaller than the other two plasmonic bending beams because the electric field intensity contribution of the top part of beam III is smaller than the beams located far from the circular hole.In addition, the paraxial limit (i.e., θ < 20°) is considered.The results are obtained in the approximation case in this study, where sin tan θ θ ≈ .Thus, the approximation conditions (i.e., θ < 20°) cannot be matched at the top part of the target curve.

Conclusion
In conclusion, a special arrangement of circular hole array is designed on the metal film to generate multiple plasmonic bending beams.Every circular hole of the array is applied as a local point-like source of SPPs and regarded as an independent radiation source.The multiple plasmonic bending beams propagate along the desired trajectory.The electric field intensity of multiple plasmonic bending beams is controlled by the polarization angle.In addition, the radius increases with the electric field intensity.These simulated results are consistent with the theoretical calculation.All these findings offer potential applications in polarization sensors, optical tweezers, and polarization-controlled plasmonic bending beam generators. ,

Fig. 1 .
Fig. 1.The geometric diagram of several radiation sources; illustration of an arbitrary bending trajectory f(y) and the tangent line, the dull-red circle points represents circle holes array (independent electric field radiation sources E n ).The configuration of the Au/SiO 2 structure is shown in Fig. 2(a).The thickness of Au film and SiO 2 substrate are fixed at 0.1 μm and 0.3 μm, respectively.Circular hole arrays are designed on the Au film.The blue block diagram illustrates a local amplification of the structure the definitions of structural parameters, where r represents the radius of circular hole, D n represents the distance between adjacent radiation sources, β represents the polarization angle between the in-plane incident electric E 0 and the x axis.W and H n represent the distances in x and y direction between adjacent radiation sources, respectively.In this study, to avoid near-field interaction between adjacent neighboring holes which can reduce the intensity of bending beams, the distance W is set as 0.6 μm.The quadratic curve (f(y) = -ay 2 ) is chosen, where constant a is 1.13 × 10 −2 .According to Eq. (3), the desired phase ( calculated when the number of circular hole is equal to 20 [Fig.2(b)].The propagation properties of the plasmonic bending beams are simulated using finite-difference time-domain method.A plane wave with 0.633 μm wavelength is incident ( + z) into the circular holes.The relative permittivity ε m of the gold film is −9.8 + 1.96i at 0.633 μm[42].

Fig. 2 .
Fig. 2. (a) The configuration of the proposed structure.The blue block diagram illustrates a local amplification of the structure and the definitions of structural parameters.(b) Calculated the distances H n of radiation source.
(a), with the polarization angle β is 0°, and the radius r is 0.2 μm.It shows that the target plasmonic bending beams (corresponding to symbol I) is generated with the propagation along the desired trajectory.The black solid curve depicts the parabolic curve f(y).Furthermore, the second (corresponding to symbol II) and third (corresponding to symbol III) plasmonic bending beams that propagate along the parabolic curve f(y) are observed.The second and third plasmonic bending beams are also satisfied with the proposed design principle.The transverse electric field intensity distributions at y = −5 μm and x = 17 μm are shown in Fig.3(b) and 3(c), respectively.These results demonstrate that multiple plasmonic bending beams are generated by especially designed circular hole arrays.

Fig. 3 .
Fig. 3. (a) Simulation electric field intensity distribution of the proposed structure with angle β = 0° and radius r = 0.2 μm.Transverse electric field intensity distributions at (b) x = 17 μm and (c) y = −5 μm..The electric field intensity distributions with different polarization angle β in x-y plane are depicted in Fig. 4, with the radius r is 0.2 μm.It shows that, as the angle β increases from 0° to 90° [Fig.4(a) -4(c)], the electric field intensities of beams I and II decrease, and the electric field intensity of beam III increases.As the angle β increases from 90° to 180° [Fig.
study, n, W, H n , and angle θ have invariable values.Thus, the electric field intensity distribution of SPP waves have a relationship with angle β.The maximum electric field intensity occurs at β = 0°.Electric fields at selected point A (x = 17 μm, y = −4 μm) of beam I, point B (x = 17 μm, y = 12 μm) of beam II, and point C (x = 2 μm, y = −5 μm) of beam III are studied to unveil this polarization-dependent property.The polar plots of the electric field intensity values of points A, B, and C with varied values of polarization angle β are shown in Fig. 4(f).The intensity of electric field at point B is about 0.14 times of the intensity of incident light wave (magnitude 1V/m).The equation to fit these data is given by 0

Figure 5 (
a) -5(c) show the electric field intensities of different circular hole radii (r = 0.1 μm, 0.15 μm, and 0.2 μm) with a polarization angle β = 0°.The transverse electric field intensity of different circular hole radii r at x = 17 μm and y = −5 μm is shown in Fig. 5(d) and 5(e), respectively.As the radius increases, the electric field intensities of the beams increase.Because as the radius of the circular hole increases, the intensity of the radiation electric field of SPPs increases, thereby resulting in the increasing intensities of incident light in the circular hole.In addition, Fig. 5(d) and 5(e) show that as the radius increases to 0.25 μm, the electric field intensities of the beams decrease.As the radius increases (the distance D n is fixed at 0.6 μm), the interaction among dipoles (i.e., the effect of the second and third terms in Eq. (4)) is enhanced and the radiation electric field intensity of SPPs decreased.The electric dipole mode converts to quartic dipole mode, which results in the change of the orientation and intensity of radiation electric field with the increasing radii.Meanwhile, the polar plots of the electric field intensity values of point A, with varied polarization angle β values and varied circular radius r, are shown in Fig. 5(f).The solid curves represent the fitting values of Eq. (

Fig. 5 .
Fig. 5. Simulation electric field intensity distribution of the proposed structure with different radius r of circle hole array: (a) r = 0.1 μm, (b) r = 0.15 μm, and (c) r = 0.2 μm.Transverse electric field intensity distributions with different radius r of circle hole array at (d) x = 17 μm and (e) y = −5 μm, (f) The polar diagram for electric field intensities at Points A with varied polarization angle β values.