Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Structured light for focusing surface plasmon polaritons

Open Access Open Access

Abstract

We propose a structureless method for focusing surface plasmon polaritons (SPPs) on a flat metal film under illumination of radially polarized cogwheel-like structured light beams. Without metal structures, the locally induced SPPs can further be propagated following the predefined patterns to form symmetric focal spots with dimensions beyond diffraction limit. Benefiting from the radial polarization, this method can be employed to pattern various center-symmetric evanescent distributions for generating SPPs reconfigurably. The SPPs will be propagating and focusing in radial directions.

©2010 Optical Society of America

1. Introduction

Surface plasmon polaritions (SPPs) [1] were highly patterned and confined near the structured metal surface and intrinsically localized in a small volume in many previous works [2]. Recently, study of focusing SPPs has attracted great attention due to the important focusing properties. Tightly focused laser beams and appropriate plasmonic lens such as nanometer-sized protrusions or holes have been applied for focusing SPPs [37]. As pioneering work, SPPs excited by a radially polarized beam created an evanescent Bessel-like standing wave on the axially symmetric plasmonic structure. The degrees of freedom for dynamic SPPs patterns were also an important point in recent research on SPPs focusing. The focused beam technique of localizing and controlling SPPs are potentially utilized for various applications such as optical probe for high-resolution imaging, sensing, and optical tweezers [8,9].

Complementary to the patterned metal surface technique, in this paper, we propose a structureless method of focusing dynamic SPPs patterns on a homogeneous gold (Au) metal film by a highly focused radially polarized structured light, “cogwheel” beam. The structured light with modified wave front (amplitude or phase) offers new opportunities for subwavelength resolution nanooptics [10]. The proposed model was described in detail based on the focused beam technique. We theoretically discuss the physics involved in the beyond-diffraction-limit focusing phenomenon and the unique characteristics of dynamic focused SPPs patterns using the vectorial diffraction theory. By performing the 3D-FDTD simulations, the proposed focusing model was verified by comparing the difference of field distributions with and without the metal film.

2. Generation and confinement of focused SPPs by structured light

Solving the Maxwell’s equations under the appropriate boundary conditions yields the SPPs dispersion relation, i.e. the frequency-dependent SPPs wave-vector,

ksp(ω)=ω/c*ε0εm(ω)/(ε0+εm(ω))
where εm(ω)and ε0are dielectric constants of metal and air respectively, and ω is the frequency of incident light. Surface plasmon resonance (SPR) excitation occurs when the wave vector of the evanescent wave klight matches that of the SPPs, which can be expressed as:
klight=k0εdsinθsp=ksp
where θsp is the resonance angle, εdis the dielectric constant of the targeted medium. In this work, a 50 nm thickness Au film (εm=5.588+2.215i at 532 nm) was deposited on the surface of glass substrate (n2 = 1.5). The medium above Au layer was chosen to be air (n1 = 1). The gap between the objective lens and the glass cover slip was filled with an index matching immersion oil (n2 = 1.5). The SPR angle θsp is about 47°.

In this paper, the focused SPPs patterns are generated by a radially polarized “cogwheel” beam. The structured light, “cogwheel” beam is the collinear superposition mode of Laguerre-Gaussian (LG) modes described by [11,12]

|umn=(|um+exp(iϕ')|un)/2
where |um> and |un> are the fields of LG0m and LG0n modes and φ’ is the relative phase of the interfering modes. The optical vortices mode with a phase singularity exp (-ilφ) has a helical wave front, where l is named as the topological charge. The “cogwheel” mode is a superposition of two vortices with equal but opposite topological charges (l = m = -n), resulting in a periodic intensity modulation around the ring circumference with a large number of intensity maxima corresponding to 2l and can propagate without changing the form in space.

In our proposed model as shown in Fig. 1 , the SPPs with complex structures can be engineered on the Au surface with a focused, radially polarized “cogwheel” beam, instead of protrusions or holes patterned on the metal film. The radial-polarized “cogwheel” beam can be produced by passing a radially polarized beam through an azimuthally modulated phase mask or a computer-generated hologram [11,12]. To a given topological charge, the radius of “cogwheel” ring can be modulated to satisfy the SPR angle requirement [13]. When the normal incident “cogwheel” beam converges toward the geometric focus, it gives rise to many diffraction-limited spots in “cogwheel” form, which contain large spectrum of wavevectors limited by the numerical aperture (NA = 1.49) of the lens. At resonant angles ± θsp, the wavevector-matched SPP waves can be generated and propagated towards the geometric center at all azimuthal directions, and the counterpropagating SPP waves in opposite directions interfere with each other to form the SPP standing waves. The focused incident beam is entirely p-polarized, thus providing cylindrical symmetric focusing. The SPPs interference patterns in the shape of “cogwheel” beam are excited at the focal plane of Au surface. Existence of these counterpropagating waves was recently shown by recording the fluorescence from single fluorophores by epiillumination scanning confocal microscope [14] and using direct measurement method aperture-type scanning near-field optical microscopy [7].

 figure: Fig. 1

Fig. 1 Schematic of dynamic SPPs pattern generated by highly focused radially polarized “cogwheel” beam at SPR condition.

Download Full Size | PDF

3. Intensity distribution of dynamic focused SPPs patterns

In a high NA system, a full vectorial formulation is derived by using vectorial Debye integral of Richards and Wolf [15,16], which inspired follow-up work in the study of high NA focusing property of doughnut beams, Bessel beam, and cylindrical vector beams in the past [17,18]. We investigate the focusing property of a monochromatic radially polarized “cogwheel” beam. The diffracted field distribution near the focal plane above the metal film is expressed by [4,17,18]:

E(r,ψ,z)=iλθ0θmax02πcosθE0(θ,ϕ)tp(θ)exp{i(k12(k2sinθ)2)1/2z}exp{ikrsinθcos(ϕψ)}[cosθcos(ϕψ)r^+sinθz^]sinθdθdϕ
where E0(θ,ϕ)is the incident “cogwheel” beam, when l = m = -n, from Eq. (3) we can derive E0(θ,ϕ)(exp(ilϕ)+exp(ilϕ))RD(l,θ)=cos(lϕ)RD(l,θ), RD(l,θ) represents the radial distribution of “cogwheel” mode. r, ψ, and z are the cylindrical coordinates of an observation points. k1 and k2 are the wave vectors in media above and below the Au layer. φ is the angle with respect to the polarization direction. The upper and lower limits of the integral angle θ denote NA of the objective lensθmax=sin1(NA)and the critical angle(sin1(ε1/ε2)), respectively. tpis the transmission coefficient of the p-polarized light for three-layer configurations, and can be define as [1]

tijkp=tijptjkpexp(ikzjd)1+rijprjkpexp(2ikzjd)

Here, d is the thickness of the Au thin film,tijp, tjkp, rijp, rjkp are the transmission and reflection Fresnel coefficients for the respective two layer interfaces, i = glass/index-matching material, j = Au and k = air. kzj is the longitudinal wave vector.

According to the Eq. (4) and Eq. (5), the electric field distribution above the Au layer, which takes into account of transmission factors for glass/Au/air configuration, can further be expressed as

Er=πiλ(I1+I2)cos(lψ)Ez=πiλI3cos(lψ)
where

I1=(i)l+1θ0θmaxcos1/2(θ)RD(θ)exp{iz(k12(k2sinθ)2)1/2}tp(θ)Jl+1(k2rsinθ)dθI2=(i)l1θ0θmaxcos1/2(θ)RD(θ)exp{iz(k12(k2sinθ)2)1/2}tp(θ)Jl1(k2rsinθ)dθI3=(i)lθ0θmaxcos1/2(θ)RD(θ)exp{iz(k12(k2sinθ)2)1/2}tp(θ)Jl(k2rsinθ)dθ

Integrations over φ have been accomplished using the identity

02πcos(lϕ)exp{ikrsinθcos(ϕψ)}dϕ=(i)n2πcos(lψ)Jl(krsinθ)
Jl(x) is the Bessel function of the first kind with order l. Under the SPR condition,k12(k2sinθsp)2<0 thus it indicates the radial and longitudinal components become evanescent along the z direction.

Intensity distribution of the focused SPPs patterns is proportional to the modulus squared of Eq. (6). The symmetric SPPs field can be focused and structured from all azimuthal directions above the flat metal layer and its shape is perfectly maintained with incident “cogwheel” beam due to cylindrical symmetry polarization as shown in Fig. 2(a) . The number of SPPs spots is different depending on the topological charges. For l = 1 and l = 2 as examples, there are 2l SPPs spots periodically arranging on the primary intensity ring circumference. Due to the dynamic character of optical vortex, it is only necessary to modulate phase values of the incident beam for desired SPPs landscapes. The fringes period is measured as 240 nm ± 5 nm originating from the interference of SPPs with opposite wavevectors, in good agreement with half wavelength of SPPs (λsp/2≈241 nm) by solving the Maxwell’s equations [Eq. (1)] under appropriate boundary conditions. FWHM of the SPPs spot on the primary intensity ring is calculated as being 0.3λ0. The peak intensity of |Ez|2 component is about six times stronger than |Er|2 component, and plays a major role in shaping the overall intensity distribution. The reason is that an angular selectivity near SPR angle is made after passing through the metal film via SPR excitation, and |Ez|2 component is predominantly enhanced due to the field enhancement effect compared to |Er|2 component.

 figure: Fig. 2

Fig. 2 (a) Normalized focused SPPs patterns intensity distributions in the vicinity of focus for radially polarized “cogwheel” beams with l = 1 and l = 2. The transverse coordinates are normalized to wavelength. (b) Transverse profiles of SPPs intensity distributions for radially polarized “cogwheel” beams of different topological charges l.

Download Full Size | PDF

The field strength transverse profiles are plotted in Fig. 2(b). In the previous works on evanescent zeroth-order Bessel beam excited by radial polarized beam, SPP wave is focused into an intense spot on the focal plane [4,5]. Whereas, in our case the focused SPPs field can be modified into the special patterns and corresponds to non-zeroth order Bessel function Jl(k2rsinθsp) derived by Eq. (7). The longitudinal component Ez with topological charge l mainly contributing to the SPPs field distribution is approximately proportional to Jl(k2rsinθsp). It also can be seen from Fig. 2(b) that the focused SPPs patterns give zero field intensity in the central region except for l = 1, this is because in Eq. (7) Jl+1(0) = 0 and Jl-1(0) = 0 for all integer l except for J0(0) = 1 when l = 1. From Eq. (6), it is seen that the focused SPPs components with the phase factor cos() retain the phase characteristic of “cogwheel” mode formed by two vortex beams with the same topological charge but opposite helicity, resulting in cosine-dependent intensity distribution (2l hotspots) in the azimuthal direction.

In the same focusing condition, the field distributions under linear polarization illumination are inhomogeneous and much weaker as shown in Fig. 3 [19]. Since SPPs can only be excited by p-polarized incident beam, the SPPs patterns appear in only one direction and the homogeneous SPPs spots cannot be obtained. So complementary to the linear-polarization, the constructive interference launched by a certain radial-polarized interfered vortex beam allows an SPPs to be focused into the evanescent interfered vortex beam and remain the phase characteristic of the incident beam. The concept presented in this paper would be extended to generate various plasmonic structured light patterns due to symmetric polarization. We can directly modify the phase structure to induce the change and rotation of SPPs patterns in a controlled manner continuously, which has more freedom than the method of patterning metal film.

 figure: Fig. 3

Fig. 3 SPPs intensity distributions for the linearly polarized “cogwheel” beam with l = 3. There are three components contributing to the total field (a): radial component (b), azimuthal component (c) and longitude component (d). The transverse coordinates are normalized to wavelength.

Download Full Size | PDF

4. Comparison of the field distributions with and without the metal

To confirm the formation of an effective focal region for SPPs focusing, we carried out FDTD simulations of SPPs propagation through flat metal film by using the FullWave module of commercial RSOFT software [20] and the results are shown in Fig. 4 . The field intensity upon the metal surface (z0 = 0) is strongly modulated [Fig. 4(a)] and the field pattern decays exponentially in the longitudinal direction [Fig. 4(b)], indicating that counterpropagating SPPs are excited. The counterpropagating SPP waves propagate towards the geometric centre and interference fringes are visible. As mentioned above, owing to the focused SPPs maintaining the singular phase character of optical vortex, the centre region is dark. The period of SPPs interference is measured as 240nm in good agreement with the vectorial diffraction theory results (Fig. 2). As depicted in Fig. 4(b), the total field strength (red curve) along the z axis above the Au film obtained from the Fig. 4(a) shows the field amplitude decreases exponentially as exp(−2|kz||z|). The decay length is estimated at about 182 nm by finding the 1/e2 point.

 figure: Fig. 4

Fig. 4 (a) 3D-FDTD results showing SPPs intensity distributions in the xz plane for radially polarized “cogwheel” beam with l = 6. The grid sizes used are Δx = Δy = 20 nm. The longitudinal nonuniform grid sizes Δz are 5 nm and 20 nm in the metal region and other regions to save the computation memory and time. (b) Normalized total field strength along the z axis (red curve) and the curve of function, exp (−2|kz||z|) (green curve).

Download Full Size | PDF

We also evaluated the electric field intensity distributions along xy plane at z = 50 nm with or without the Au film. It is seen from Fig. 5(a) that the SPPs spots are shaped like many concentric “cogwheel” rings of diminishing intensity. Figure 5(b) shows the focal intensity distribution at the glass/air interface without metal film. FWHM of the spot on the primary intensity ring without the metal film is about 0.6λ while with the metal film it is only 0.3λ which is beyond the conventional diffraction limit. The reason is that the electric field contains both propagating and evanescent components in the case without the metal film, the incident energy below the critical angle can entirely transmit, but with the metal film, the propagating part of the incident energy below the SPR angle is reflected back into the objective lens and is efficiently blocked. The excitation scheme can easily pattern each of the SPPs spots in an area approximately 0.09λ2 from a single beam and has potential applications in biological molecule manipulation and subwavelength optical sorting, since evanescent field including SPPs contains high spatial frequencies and reduced lateral size in comparison with propagating fields.

 figure: Fig. 5

Fig. 5 Electric field intensity distributions in the xy plane at z = 50 nm with the Au film (a) and without the Au film (b) for radially polarized “cogwheel” beams with topological charge 6 by using 3D-FDTD method. The grid sizes used in the xy plane are Δx = Δy = 10 nm.

Download Full Size | PDF

5. Conclusions

In conclusion, the structured light beams are versatile sources for generation of desired focused SPPs patterns confined beyond the diffraction limit, which no protrusion or hole is necessary. The focusing SPPs by “cogwheel” beam in radial polarization can be structured into the shape of concentric “cogwheel” rings at a planar dielectric-metal interface. We demonstrate that cogwheel-like structured light for dynamic SPPs array is a reconfigurable solution among the contenders. Benefiting from the symmetric polarization, the concept presented in this paper could be extended to generate various shapes of plasmonic structured beams through designing different phase masks, which may offer the new opportunity for the high resolution imaging microscopy and optical tweezers. The focused SPPs spots field gradient force for the dielectric particle will be our interest for investigation in future.

Acknowledgements

This work was partially supported by the National Natural Science Foundation of China under Grant No. (10974101 and 60778045), Ministry of Science and Technology of China under Grant no. 2009DFA52300 for China-Singapore collaborations and National Research Foundation of Singapore under Grant No. NRF-G-CRP 2007-01.

References and links

1. H. Raether, Surface-Plasmons on Smooth and Rough Surfaces and on Gratings, Springer Tracts in Modern Physics (Springer, Berlin, 1988).

2. M. Righini, A. S. Zelenina, C. Girard, and R. Quidant, “Parallel and selective trapping in a patterned plasmonic landscape,” Nat. Phys. 3(7), 477 (2007). [CrossRef]  

3. G. M. Lerman, A. Yanai, and U. Levy, “Demonstration of nanofocusing by the use of plasmonic lens illuminated with radially polarized light,” Nano Lett. 9(5), 2139–2143 (2009). [CrossRef]   [PubMed]  

4. Q. Zhan, “Evanescent Bessel beam generation via surface plasmon resonance excitation by a radially polarized beam,” Opt. Lett. 31(11), 1726–1728 (2006). [CrossRef]   [PubMed]  

5. W. B. Chen and Q. Zhan, “Realization of an evanescent Bessel beam via surface plasmon interference excited by a radially polarized beam,” Opt. Lett. 34(6), 722–724 (2009). [CrossRef]   [PubMed]  

6. H. Kano, S. Mizuguchi, and S. Kawata, “Excitation of surface-plasmon polaritons by a focused laser beam,” J. Opt. Soc. Am. B 15(4), 1381–1386 (1998). [CrossRef]  

7. A. Bouhelier, F. Ignatovich, A. Bruyant, C. Huang, G. Colas des Francs, J.-C. Weeber, A. Dereux, G. P. Wiederrecht, and L. Novotny, “Surface plasmon interference excited by tightly focused laser beams,” Opt. Lett. 32(17), 2535–2537 (2007). [CrossRef]   [PubMed]  

8. R. Vander and S. G. Lipson, “High-resolution surface-plasmon resonance real-time imaging,” Opt. Lett. 34(1), 37–39 (2009). [CrossRef]  

9. K. J. Moh, X.-C. Yuan, J. Bu, S. W. Zhu, and B. Z. Gao, “Radial polarization induced surface plasmon virtual probe for two-photon fluorescence microscopy,” Opt. Lett. 34(7), 971–973 (2009). [CrossRef]   [PubMed]  

10. L. David, Andrews, Structured Light and Its Applications, Academic Press is an imprint of Elsevier (Elsevier, USA, 2008).

11. D. W. Zhang and X.-C. Yuan, “Entangled double-helix phase,” Opt. Lett. 28(20), 1864–1866 (2003). [CrossRef]   [PubMed]  

12. A. Jesacher, S. Fürhapter, S. Bernet, and M. Ritsch-Marte, “Size selective trapping with optical “cogwheel” tweezers,” Opt. Express 12(17), 4129–4135 (2004), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-12-17-4129. [CrossRef]   [PubMed]  

13. J. Lin, X.-C. Yuan, S. H. Tao, and R. E. Burge, “Variable-radius focused optical vortex with suppressed sidelobes,” Opt. Lett. 31(11), 1600–1602 (2006). [CrossRef]   [PubMed]  

14. F. D. Stefani, K. Vasilev, N. Bocchio, N. Stoyanova, and M. Kreiter, “Surface-plasmon-mediated single-molecule fluorescence through a thin metallic film,” Phys. Rev. Lett. 94(2), 023005 (2005). [CrossRef]   [PubMed]  

15. B. Richards and E. Wolf, “Electromagnetic diffraction in optical systems II. Structure of the image field in an aplanatic system,” Proc. R. Soc. Lond. A Math. Phys. Sci. 253(1274), 358–379 (1959). [CrossRef]  

16. E. Wolf, “Electromagnetic Diffraction in Optical Systems. I. An Integral Representation of the Image Field,” Proc. R. Soc. Lond. A Math. Phys. Sci. 253(1274), 349–357 (1959). [CrossRef]  

17. K. S. Youngworth and T. G. Brown, “Focusing of high numerical aperture cylindrical-vector beams,” Opt. Express 7(2), 77–87 (2000), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-7-2-77. [CrossRef]   [PubMed]  

18. Q. Zhan and J. R. Leger, “Focus shaping using cylindrical vector beams,” Opt. Express 10(7), 324–331 (2002), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-10-7-324. [PubMed]  

19. Z. J. Hu, X.-C. Yuan, S. W. Zhu, G. H. Yuan, P. S. Tan, J. Lin, and Q. Wang, ““Dynamic surface Plasmon patterns generated by reconfigurable “cogwheel-shaped” beams,” Appl. Phys. Lett. 93(18), 181102 (2008). [CrossRef]  

20. http://www.rsoftdesign.com/

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1
Fig. 1 Schematic of dynamic SPPs pattern generated by highly focused radially polarized “cogwheel” beam at SPR condition.
Fig. 2
Fig. 2 (a) Normalized focused SPPs patterns intensity distributions in the vicinity of focus for radially polarized “cogwheel” beams with l = 1 and l = 2. The transverse coordinates are normalized to wavelength. (b) Transverse profiles of SPPs intensity distributions for radially polarized “cogwheel” beams of different topological charges l.
Fig. 3
Fig. 3 SPPs intensity distributions for the linearly polarized “cogwheel” beam with l = 3. There are three components contributing to the total field (a): radial component (b), azimuthal component (c) and longitude component (d). The transverse coordinates are normalized to wavelength.
Fig. 4
Fig. 4 (a) 3D-FDTD results showing SPPs intensity distributions in the xz plane for radially polarized “cogwheel” beam with l = 6. The grid sizes used are Δx = Δy = 20 nm. The longitudinal nonuniform grid sizes Δz are 5 nm and 20 nm in the metal region and other regions to save the computation memory and time. (b) Normalized total field strength along the z axis (red curve) and the curve of function, exp (−2|kz||z|) (green curve).
Fig. 5
Fig. 5 Electric field intensity distributions in the xy plane at z = 50 nm with the Au film (a) and without the Au film (b) for radially polarized “cogwheel” beams with topological charge 6 by using 3D-FDTD method. The grid sizes used in the xy plane are Δx = Δy = 10 nm.

Equations (8)

Equations on this page are rendered with MathJax. Learn more.

ksp(ω)=ω/c*ε0εm(ω)/(ε0+εm(ω))
klight=k0εdsinθsp=ksp
|umn=(|um+exp(iϕ')|un)/2
E(r,ψ,z)=iλθ0θmax02πcosθE0(θ,ϕ)tp(θ)exp{i(k12(k2sinθ)2)1/2z}exp{ikrsinθcos(ϕψ)}[cosθcos(ϕψ)r^+sinθz^]sinθdθdϕ
tijkp=tijptjkpexp(ikzjd)1+rijprjkpexp(2ikzjd)
Er=πiλ(I1+I2)cos(lψ)Ez=πiλI3cos(lψ)
I1=(i)l+1θ0θmaxcos1/2(θ)RD(θ)exp{iz(k12(k2sinθ)2)1/2}tp(θ)Jl+1(k2rsinθ)dθI2=(i)l1θ0θmaxcos1/2(θ)RD(θ)exp{iz(k12(k2sinθ)2)1/2}tp(θ)Jl1(k2rsinθ)dθI3=(i)lθ0θmaxcos1/2(θ)RD(θ)exp{iz(k12(k2sinθ)2)1/2}tp(θ)Jl(k2rsinθ)dθ
02πcos(lϕ)exp{ikrsinθcos(ϕψ)}dϕ=(i)n2πcos(lψ)Jl(krsinθ)
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.