Phase-Resolved Observation of the Gouy Phase Shift of Surface

Mar 20, 2017 - The space-resolved Gouy phase shift for focused surface plasmon polaritons (SPPs) in the optical regime is experimentally demonstrated ...
0 downloads 11 Views 3MB Size
Subscriber access provided by HACETTEPE UNIVERSITESI KUTUPHANESI

Article

Phase Resolved Observation of the Gouy Phase Shift of Surface Plasmon Polaritons Tobias Birr, Tim Fischer, Andrey B. Evlyukhin, Urs Zywietz, Boris N. Chichkov, and Carsten Reinhardt ACS Photonics, Just Accepted Manuscript • DOI: 10.1021/acsphotonics.6b00999 • Publication Date (Web): 20 Mar 2017 Downloaded from http://pubs.acs.org on March 26, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

ACS Photonics is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 7

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

Phase Resolved Observation of the Gouy Phase Shift of Surface Plasmon Polaritons Tobias Birr,



Tim Fischer,



Chichkov,

†Laser

Andrey B. Evlyukhin,





Urs Zywietz,

and Carsten Reinhardt



Boris N.

∗, †, ‡

Zentrum Hannover e.V., Hollerithallee 8, D-30419 Hannover, Germany

‡Hochschule

Bremen, Neustadtswall 30, D-28199 Bremen

E-mail: [email protected]

Abstract

vided a simple derivation of the Gouy phase shift for multimode beams, which for fundamental mode in two dimensional case φGouy,2D can be expressed by 3

The space-resolved Gouy phase shift for focused surface plasmon polaritons (SPPs) in the optical regime is experimentally demonstrated for the rst time. SPPs are excited by CW laser radiation focused onto a dielectric ridge placed on a cover glass surface coated by a silver lm. The Gouy phase measurements are performed using interference between a reference plane wave and leakage radiation. The measured Gouy shift has a value of π/2 corresponding to the twodimensional nature of surface plasmon polaritons.

φGouy,2D =

z 1 arctan , 2 zR

(1)

where z is the propagation direction and zR = n π w02 / λ is the Rayleigh length for the Gaussian beam with the wavelength λ, focus width w0 , and the medium refractive index n. The Gouy shift has already been measured and reported for surface plasmons and phonon polaritons in the THz range. 48 The direct phase resolved observation in the optical range is still missing. This paper provides the rst direct observation of the Gouy shift by comparison of the phase resolved focused and plane wave SPP leakage radiation images, as can be seen in Fig. 1. Since SPPs propagate at the interface of a metal and dielectric medium and represent a pure two dimensional system the total measured phase shift of π/2 is observed in this work. The images were recorded with the leakage radiation microscopy technique which enables phase and time sensitive imaging of the SPPs similar to that introduced by Gorodetski et al. 9 Simple theoretical modelling provides straightforward explanation of the observed results and can be applied for all kinds of beam propagation phenomena. Because the Gouy shift occurs in all kinds of focused waves, it plays a fundamental role in

Keywords Gouy phase shift, Surface plasmon-polaritons, Leakage radiation microscopy In the late 19th century, Louis Georges Gouy has theoretically predicted and experimentally veried the existence of a longitudinal phase shift for light waves passing through the focal point with respect to the phase front of collimated light beams. 1,2 This phase shift became widely known as the Gouy phase shift. The origin can be explained with the uncertainty relation in the focus region, which is ∆k∆x = 1/2 for the two dimensional case. 3 Since the convergent waves have a nite spatial extent, when passing a caustic, the wave vector has to consist of a superposition of transversal and longitudinal components. Feng and Winful pro-

ACS Paragon Plus Environment

1

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Exp

Exp

to the phase front of the exciting beam at the ridge surface. 21 When the laser focus is far away from the excitation ridge, the beam curvature can be neglected and the paraxial part of the SPP beam can be treated as a plane wave. When the laser beam focus is slightly behind the excitation ridge, the excited SPP beam will be converging and focussed. 22,23 The LRM setup, shown in Fig. 2, consists of a 40 x microscope objective with 0.65 NA for focusing the laser beam on the sample, an 40 x oil immersion microscope objective with 1.3 NA for image generation, and a CMOS-Camera. An additional 4-f setup is implemented to provide a possibility for blocking out the excitation laser beam on the recorded images. For enabling phase resolved measurements, about 10 % of the excitation beam intensity is coupled out to a bypass line as a probe beam (PB), before entering the LRM setup. The bypass line oers the same optical path length as the LRM setup. At the output both beams are combined in front of the CMOS. For variation of the optical path length in the bypass line, two mirrors are mounted on a linear micrometer stage. A further mirror behind the linear stage in the bypass line is used for the adjustment of the incident angle of the probe beam onto the CMOS, to align it to the incident angle of the leakage radiation (LR). For maintaining a homogeneous illumination of the whole CMOS, the PB is expanded with an Galilei telescope. In this conguration, the PB and the LR will interfere directly on the CMOS. Thus, the recorded intensity resembles a time integrated superposition of the leakage radiation and the probe beam. This can be expressed by

1.0

5µm

Page 2 of 7

-0.3

Figure 1: Experimental demonstration of the spatially resolved two dimensional Gouy phase shift for SPPs. Leakage radiation images for focused and plane wave SPPs are shown in the upper and lower parts, respectively. The phase resolution was achieved by interference with a probe laser beam. Regions of positive and negative interference are shown by red and blue colors. modern physics. In optical systems, this effect is important for the correct dimensioning of laser cavities, 10 nonlinear eects, 11 attosecond pulse generation, 12 and optical coherence tomography. 13 In SPP systems it should be taken into account for observations of localized nonlinear eects and development of integrated circuits, bio-sensors, and data storage devices. 14

Experimental details For experimental observations of the Gouy phase shift, the SPP system, consisting of a 55 nm thin silver lm thermally evaporated on a standard microscopic cover glass (Menzel), is investigated by leakage radiation microscopy (LRM). 15,16 To excite SPPs, a TM laser beam is focussed on a polymer ridge with a width and height of about 300 nm underneath the silver lm. The polymer ridge structures are produced from Ormosil 17 by means of microscopeprojection-photolithography. 1820 As excitation laser source a titanium-sapphire oscillator (Kapteyn Murnane Labs) is used in CW mode at 780 nm center wavelength. When the excitation laser beam has Gaussian shape, the excited SPP beam will resemble a two dimensional projection of this shape corresponding

ICMOS (x, z, t) ∝ Z

τfov

|ELR (x, z, t) + EPB (t)|2 dt,

(2)

0

where z is the SPP propagation direction and x is perpedicular to z . Because the exposure time of the CMOS is much longer than the SPP propagation time for the eld of view (fov) τSPP = lLfov /cSPP , the integral gives a time average signal of the interference pattern. Since

ACS Paragon Plus Environment

2

Page 3 of 7

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

Figure 2: Scheme of the used setup for phase resolved measurements: about 10 % of the excitation beam intensity is coupled out to a bypass line as a probe beam (PB), before entering the Leakage radiation microscope setup and is recombined with the leakage radiation (LR) in front of the CMOS camera. For maintaining a homogeneous illumination of the whole CMOS, the PB is expanded with an Galilei telescope. An additional 4-f setup is implemented to provide the possibility for blocking out the excitation laser beam on the recorded images. In numerical calculations, the expressions for the respective electric elds are rst superposed as shown in Eq. 2. The square of the absolute values is than directly computed into threedimensional (x, z, t) matrices, which are nally integrated over the time axis. This method allows to use the superior vector computation of numPy (Python), Matlab or Octave. 25,26 As spatial steps of the mentioned matrices, the resolution of the experimental taken LR images are used. To calculate a sucient time discretization, the Courant factor, normally applied for FDTD simulations, 27 is taken into account.

the probe beam is illuminating the CMOS homogeneously, it can be considered as a time dependent plane wave (3)

EPB (t) ∝ exp (−i ω t),

with the the angular frequency ω . The leakage radiation images the SPP electrical eld and can therefore be represented as two dimensional distribution of a Gaussian beam. The electrical eld of the Gaussian SPP beam at the metal surface can be presented by 24

ELR,foc (x, z, t) ∝ 

  exp  −

x2

Results and discussion

 + ikSPP z − iωt 

iz w02 (1 + ) z rR iz 1+ zR

,

The SPP resembles a two dimensional projection of the incident laser beam. The SPP beam can be focused or collimated to a plane wave by altering the distance of the laser beam focus with respect to the polymer ridge at the sample. 22,23 The LR images corresponding to both cases are compared in Fig. 1. Here, a perfect overlap of the wave fronts for the plane and focused SPP waves is clearly visible in front of the focus position. After the focus, one can observe the phase mismatch between the wave fronts

(4)

where w0 is the beam waist in the focus and zR is the Rayleigh length. The SPP wave vector kSPP is assumed to be complex, to include propagation losses into the model. The plane SPP wave can be expressed by

ELR,pl (x, z, t) ∝ exp (i kSPP z − i ω t),

(5)

ACS Paragon Plus Environment

3

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(a)

Exp

5µm

Fit 1.0

Page 4 of 7

Exp

(b) (c)

3D

2D

-0.3

Fit

5µm

Figure 3: Numerical calculations for SPP waves were tted to the observed experimental interference pattern for (a) plane wave and (b) focused SPP beams, where regions of positive and negative interference are shown by red and blue colors. (c) The Gouy phase shift was determined by comparing the distances between the wave front maxima for the plane and focused SPP waves on the propagation axis. The conversion of these distances to the phase shifts in radians are shown as red dots. The theoretically expected arc tangent function (with zR = 1.5 µm) is added as dashed black curves for the two dimensional ( π/2 phase shift) and the three dimensional ( π phase shift) cases. corresponding to the Gouy phase shift. Numerical calculations for plane SPP waves were tted to the observed experimental interference pattern. As a tting parameter, the complex SPP wave vector was used with an initial value of 763.3 nm, given by the semi innite approach 28 and the Lorentz-Drude model 29 for the dielectric function of silver. The result of this t is shown in Fig.3(a) and demonstrates a perfect coincidence between the measured and tted SPP plane waves. This procedure gave for the SPP wavelength the value of 752 nm. The dierence to the calculated value is about 1.5 %, proving a very good accuracy of the applied scaling factor for the microscope images. The determined SPP wave vector was further used for interpretation of experiments with focused SPP Gaussian beams. In this case to t numerical data to the experimental results, the beam waist w0 as the only tting parameter

was used. It was determined by the t to be w0 = 589.7 nm. This is in very good correlation with the expected diraction limited spot size of d = λ/2NA = 600 nm. Comparison of the experimental and numerical results is shown in Fig. 3(b), where again a perfect agreement is visible. The Gouy phase shift was determined by comparing the distances between the wave front maxima for the plane and focused SPP waves on the propagation axis. These distances were converted to the phase shifts in radians shown in Fig. 3(c) as red dots. The theoretically expected arc tangent function, see Eq. 1, with the Rayleigh length zR =1.5 µm is shown as a dashed black curves for the two dimensional (π/2 phase shift) and the three dimensional (π phase shift) cases, respectively. One can observe a perfect agreement between the measured and computed results for the two dimensional case, emphasizing the two dimen-

ACS Paragon Plus Environment

4

Page 5 of 7

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

sional nature of the shown phenomenon.

(5) McGowan, R. W.; Cheville, R. A.; Grischkowsky, D. Direct Observation of the Gouy Phase Shift in THz Impulse Ranging Direct Observation of the Gouy Phase Shift in THz Impulse Ranging. Appl. Phys. Lett. 2000, 76, 670672.

Conclusion Using the leakage radiation microscopy, the space-resolved Gouy phase shift for focused surface plasmon polaritons (SPPs) has been experimentally investigated for the rst time. These measurements have been performed using interference between a reference plane wave and leakage radiation. A perfect agreement between the measured and theoretically predicted arc tangent function has been demonstrated. The measured total Gouy shift of π/2 corresponds to the two-dimensional nature of surface plasmon polaritons. In the development of SPP systems, the Gouy shift has to be taken into account in design of integrated circuits, bio-sensors, data storage and nonlinear devices.

(6) Feurer, T.; Stoyanov, N. S.; Ward, D. W.; Nelson, K. A. Direct Visualization of the Gouy Phase by Focusing Phonon Polaritons. Phys. Rev. Lett. 2002, 88, 257402. (7) Zhu, W.; Agrawal, A.; Nahata, A. Measurements of the Gouy Phase Shift for Surface Plasmons. America (NY). 2007, 1, 56245624. (8) Wang, X.; Sun, W.; Cui, Y.; Ye, J.; Feng, S.; Zhang, Y. Complete Presentation of the Gouy Phase Shift with the THz Digital Holography. Opt. Express 2013, 21, 2387 2393.

ACKNOWLEDGEMENT

The authors acknowledge nancial support of this work from the Deutsche Forschungsgemeinschaft (DFG, SPP1391: "Ultrafast Nanooptics" and SFB/TRR123: "Planar Optronic Systems") and support of the Laboratorium of Nano- and Quantenengineering (LNQE). The authors further acknowledge support from Hannover School of Nanotechnology (HSN).

(9) Gorodetski, Y.; Chervy, T.; Wang, S.; Hutchison, J. A.; Drezet, A.; Genet, C.; Ebbesen, T. W. Tracking Surface Plasmon Pulses Using Ultrafast Leakage Imaging. Optica 2016, 3, 4853. (10) Feng, S.; Winful, H. G. Spatiotemporal Structure of Isodiracting Ultrashort Electromagnetic Pulses. Phys. Rev. E 2000, 61, 862873.

References

(11) Constant, E.; Garzella, D.; Breger, P.; Mével, E.; Dorrer, C.; Le Blanc, C.; Salin, F.; Agostini, P. Optimizing High Harmonic Generation in Absorbing Gases: Model and Experiment. Phys. Rev. Lett. 1999, 82, 16681671.

(1) Gouy, L. G. Sur Une Propriété Nouvelle Des Ondes Lumineuses. Comptes Rendus des Seances L'Académie des Sci. 1890, 1, 12511253. (2) Gouy, L. G. Sur La Propagation Anomale Des Ondes. Ann. Chim. Phys 1891, 24, 145213. (3) Feng, S.; Winful, H. G. Physical Origin of the Gouy Phase Shift. Opt. Lett. 2001, 26, 485487.

(12) Lindner, F.; Paulus, G. G.; Walther, H.; Baltu²ka, A.; Goulielmakis, E.; Lezius, M.; Krausz, F. Gouy Phase Shift for FewCycle Laser Pulses. Phys. Rev. Lett. 2004, 92, 113001.

(4) Run, A. B.; Rudd, J. V; Whitaker, J. F.; Feng, S.; Winful, H. G. Direct Observation of the Gouy Phase Shift with Single-Cycle Terahertz Pulses. Phys. Rev. Lett. 1999, 83, 34103413.

(13) Johnson, J. L.; Dorney, T. D.; Mittleman, D. M. Enhanced Depth Resolution in Terahertz Imaging Using Phase-Shift Interferometry. Appl. Phys. Lett. 2001, 78, 835 837.

ACS Paragon Plus Environment

5

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(14) Radko, I. P.; Bozhevolnyi, S. I.; Evlyukhin, A. B.; Boltasseva, A. Surface Plasmon Polariton Beam Focusing with Parabolic Nanoparticle Chains. Opt. Express 2007, 15, 65766582.

Page 6 of 7

F. J.; Boltasseva, A. Eciency of Local Surface Plasmon Polariton Excitation on Ridges. Phys. Rev. B - Condens. Matter Mater. Phys. 2008, 78, 17. (22) Kiyan, R.; Reinhardt, C.; Passinger, S.; Stepanov, A. L.; Hohenau, A.; Krenn, J. R.; Chichkov, B. N. Rapid Prototyping of Optical Components for Surface Plasmon Polaritons. Opt. Express 2007, 15, 4205 4215.

(15) Drezet, A.; Hohenau, A.; Koller, D.; Stepanov, A.; Ditlbacher, H.; Steinberger, B.; Aussenegg, F. R.; Leitner, A.; Krenn, J. R. Leakage Radiation Microscopy of Surface Plasmon Polaritons. Mater. Sci. Eng. B Solid-State Mater. Adv. Technol. 2008, 149, 220229.

(23) Evlyukhin, A. B.; Bozhevolnyi, S. I.; Stepanov, A. L.; Kiyan, R.; Reinhardt, C.; Passinger, S.; Chichkov, B. N. Focusing and Directing of Surface Plasmon Polaritons by Curved Chains of Nanoparticles. Opt. Express 2007, 15, 1666716680.

(16) Sardana, N.; Birr, T.; Schlenker, S.; Reinhardt, C.; Schilling, J. Surface Plasmons on Ordered and Bi-Continuous Spongy Nanoporous Gold. New J. Phys. 2014, 16, 63053.

(24) Kraus, H. G. HuygensFresnelKirchho WaveFront Diraction Formulation: Paraxial and Exact Gaussian Laser Beams. J. Opt. Soc. Am. A 1990, 7, 47.

(17) Ovsianikov, A.; Viertl, J.; Chichkov, B.; Oubaha, M.; MacCraith, B.; Sakellari, I.; Giakoumaki, A.; Gray, D.; Vamvakaki, M.; Farsari, M.; Fotakis, C. Ultra-Low Shrinkage Hybrid Photosensitive Material for Two-Photon Polymerization Microfabrication. ACS Nano 2008, 2, 22572262.

(25) Cai, X.; Langtangen, H. P.; Moe, H. On the Performance of the Python Programming Language for Serial and Parallel Scientic Computations. Sci. Program. 2005, 13, 3156.

(18) Birr, T.; Zywietz, U.; Chhantyal, P.; Chichkov, B. N.; Reinhardt, C. Ultrafast Surface Plasmon-Polariton Logic Gates and Half-Adder. Opt. Express 2015, 23, 3175531765.

(26) Chaves, J. C.; Nehrbass, J.; Guilfoos, B.; Gardiner, J.; Ahalt, S.; Krishnamurthy, A.; Unpingco, J.; Chalker, A.; Warnock, A.; Samsi, S. Octave and Python: HighLevel Scripting Languages Productivity and Performance Evaluation. In HPCMP Users Group Conference, 2006; IEEE, 2006; 429434.

(19) Birr, T.; Zywietz, U.; Fischer, T.; Chhantyal, P.; Evlyukhin, A. B.; Chichkov, B. N.; Reinhardt, C. Ultrafast Surface Plasmon-Polariton Interference and Switching in Multiple Crossing Dielectric Waveguides. Appl. Phys. B 2016, 122, 164.

(27) Inan, U. S.; Marshall, R. A. Numerical Electromagnetics: The FDTD Method; Cambirdge Iniversity Press: Cambridge, 2011.

(20) Love, J. C.; Wolfe, D. B.; Jacobs, H. O.; Whitesides, G. M. Microscope Projection Photolithography for Rapid Prototyping of Masters with Micron-Scale Features for Use in Soft Lithography. Langmuir 2001, 17, 60056012.

(28) Raether, H. Surface Plasmons on Smooth Surfaces; Springer, 1988. (29) Rakic, A. D.; Djurisic, A. B.; Elazar, J. M.; Majewski, M. L. Optical Properties of Metallic Films for Vertical-Cavity Optoelectronic Devices. Appl. Opt. 1998, 37, 52715283.

(21) Radko, I. P.; Bozhevolnyi, S. I.; Brucoli, G.; Martín-Moreno, L.; García-Vidal,

ACS Paragon Plus Environment

6

Experiment

2D Gaussian focus Experiment

Page 7 of 7

1 Calculated fit Experiment 2

ACS Photonics

Experiment Plane wave

3 Gouy ACS Paragon Plus Environment phase shift 4 5 5µm Calculated fit 6