Analogy of electromagnetically induced transparency in plasmonic nanodisk with a square ring resonator

Analogy of electromagnetically induced transparency in plasmonic nanodisk with a square ring resonator

JID:PLA AID:23488 /SCO Doctopic: Plasma and fluid physics [m5G; v1.161; Prn:27/10/2015; 10:41] P.1 (1-6) Physics Letters A ••• (••••) •••–••• 1 C...

2MB Sizes 1 Downloads 70 Views

JID:PLA

AID:23488 /SCO Doctopic: Plasma and fluid physics

[m5G; v1.161; Prn:27/10/2015; 10:41] P.1 (1-6)

Physics Letters A ••• (••••) •••–•••

1

Contents lists available at ScienceDirect

2

67 68

3

Physics Letters A

4 5

69 70 71

6

72

www.elsevier.com/locate/pla

7

73

8

74

9

75

10

76

11 12 13

Analogy of electromagnetically induced transparency in plasmonic nanodisk with a square ring resonator

14 15 16 17 18 19

81

Guangdong Provincial Key Laboratory of Nanophotonic Functional Materials and Devices, School for Information and Optoelectronic Science and Engineering, South China Normal University, Guangzhou 510006, China

84

26 27

30 31 32

83 85 87

a r t i c l e

i n f o

a b s t r a c t

88 89

Article history: Received 14 June 2015 Accepted 20 October 2015 Available online xxxx Communicated by F. Porcelli

28 29

82

86

23 25

79

Xianping Li, Zhongchao Wei, Yuebo Liu, Nianfa Zhong, Xiaopei Tan, Songsong Shi, Hongzhan Liu, Ruisheng Liang

21

24

78 80

20 22

77

Keywords: Surface plasmonic Electromagnetic optics Waveguide Integrated optics devices

We have demonstrated the analogy of electromagnetically induced transparency in plasmonic nanodisk with a square ring resonator. A reasonable analysis of the transmission feature based on the temporal coupled-mode theory is given and shows good agreement with the Finit-Difference Time-Domain simulation. The transparency window can be easily tuned by changing the geometrical parameters and the insulator filled in the resonator. The transmission of the resonator system is close to 80% and the full width at half maximum is less than 46 nm. The sensitivity of the structure is about 812 nm/RIU. These characteristics make the new system with potential to apply for optical storage, ultrafast plasmonic switch and slow-light devices. © 2015 Published by Elsevier B.V.

90 91 92 93 94 95 96 97 98

33

99

34

100

35

101

36 37

102

1. Introduction

38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Electromagnetically induced transparency (EIT) as a thrilling and counterintuitive phenomenon shows a narrow spectral transparency region in a broad absorption regime, which occurs in three energy level atomic systems due to the quantum destructive interference between the excitation path-ways to the atomic upper level [1,2]. This phenomenon has received much attention due to its interesting physics [3,4] and potential applications such as the transfer of quantum correlations [5], slow light propagation [6], and so on. However, the demands of stable gas lasers and low temperature environment severely hamper the implementation of EIT in chip-scale application. Recently, EIT-like has been observed in many systems such as photonic crystal nanocavity [7], hybrid plasmonic waveguides [8,9], symmetric planar metamaterial [10], and resonator [11–13]. Among these systems, the EIT-like effect based on the surface plasmon polaritons (SPPs) has been favored more and more by researchers. Surface plasmon polaritons are a kind of electromagnetic waves propagating along the metal–dielectric interface with an exponential decaying field on both sides [14]. Due to the capabilities of overcoming the classical diffraction limit and manipulating light in the nanoscale domain, SPPs have been considered as one of the most promising energy and information carriers [15,16]. So

61 62 63 64 65 66

E-mail address: [email protected] (Z. Wei). http://dx.doi.org/10.1016/j.physleta.2015.10.035 0375-9601/© 2015 Published by Elsevier B.V.

far, numerous EIT-like devices based on SPPs such as Fabry–Perot resonator [17], Fano resonances [18,19], asymmetric dual sidecoupled cavities [20], stub waveguide with ring resonator [21] have been investigated theoretically and demonstrated experimentally. Among multitudinous plasmonic structures, the plasmonic analogous to EIT in the metal–insulator–metal (MIM) waveguide systems [22–24] draw much attention because of the deep-sub-wavelength confinement of light and relatively simple fabrication [25]. An intriguing potential application for such plasmonic analogue of EIT structures is their use in light manipulation and transmission in nanoscale devices. In this paper, a new kind of analogies of EIT in metal–insulator– metal plasmonic waveguide consisting of a square ring resonator coupled with a nanodisk is proposed and numerically investigated. We firstly briefly reviewed the properties of basic geometry, then theoretically analyzed this structure by the temporal coupled-model theory (CMT). Especially, in order to understand the relationship between the transmission characteristics and geometrical parameters of the structure and refractive index of insulator, we acquired the intuitive images by Finit-Difference Time-Domain (FDTD) simulation. The results show that the transparency window can be easily tuned by changing the geometrical parameters of the structure and the EIT-like resonant peak has a linear relationship with refractive index. These results may inspire interest in nanoscale wavelength-filter, slow light devices and optical switching elements in highly integrated optical circuit.

103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132

JID:PLA

AID:23488 /SCO Doctopic: Plasma and fluid physics

[m5G; v1.161; Prn:27/10/2015; 10:41] P.2 (1-6)

X. Li et al. / Physics Letters A ••• (••••) •••–•••

2

1

67

2

68

3

69

4

70

5

71

6

72

7

73

8

74

9

75

10

76

11

77

12

78

13

79

14

80

15

81

16

82

17

83

18

84

19

85

20 21 22

86

Fig. 1. (a) Schematic of MIM waveguide directly coupled to a nanodisk resonator with gap = 0. (b) Transmission spectra for MIM waveguide directly coupled to nanodisk, with d = 50 nm and r = 57 nm.

23 24 25

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

2. Theoretical analysis and discussion of MIM waveguide with one nanodisk resonator To start, let us firstly briefly review the properties of basic geometry as shown in Fig. 1(a), the nanoscale plasmonic resonator system with one nanodisk directly coupled to a bus waveguide. The dielectric layer here is assumed to be air (n = 1) and the background material in blue is silver. The permittivity of the metallic region can be described by the Drude mode [26]:

ω2p εm (ω) = ε∞ − ω(ω + i γ )

(1)

where ε∞ = 3.7 is the dielectric constant at infinite frequency, ω p = 9.1 eV is the bulk plasma frequency of free conduction electrons, γ = 0.018 eV is the electron collision frequency and ω is the angular frequency of incident light in vacuum. Fig. 1(b) shows the transmission spectra for the MIM waveguide directly with one nanodisk resonator. According to CMT [27,28], the transmission of the system supporting a resonant mode can be indicated as:

(ω − ω0 )2 + (1/τ0 ) T= (ω − ω0 )2 + (1/τ0 + 1/τe )

(2)

where τ0 and τe are the decay rate due to the intrinsic loss and waveguide coupling loss, respectively. It is obvious that the mini2

(1/τ )

mum transmission T min = (1/τ +10/τ )2 , when ω = ω0 . As the size e 0 of the structure is much smaller than incident wavelength, the intrinsic loss can be neglected, then T min = 0, which is accordance well with the numerical simulation.

54 55

58 59 60 61 62 63 64 65 66

ring resonator, the transmission spectrum exhibits a spectral dip at the resonance wavelength of 825 nm due to the destructive interference between the incident wave and escaped power from the resonator [29]. However as the nanodisk and square ring are combined into a composite structure as shown in Fig. 2(a), the transmission spectrum exhibits a narrow transparency peak in the transmission dip near the resonant wavelength of the ring resonator, the EIT-like transparency peak is close to 80% and the full width at half maximum (FWHM) is less than 46 nm. In order to have an insight into the physical mechanism behind the EIT-like effect in the new proposed system, we plot the electric field. Figs. 2(c) and (d) show the field distributions at the EIT-like transparency peak of 825 m without and with the square ring resonator. As shown in Fig. 2(c), the incident lights are reflected in the nanodisk-shaped waveguide, so that the whole modes are locked in the nanodisk, whereas they pass through the waveguide with ring resonator, as can be seen in Fig. 2(d). The results are consistent with the spectral response in Fig. 2(b). It is worth noting that the electromagnetic field in the nanodisk is weakened due to the interference effect when the disk couples with the square ring cavity. In order to analyze the EIT-like phenomenon in detail, we introduce a CMT-based transmission line theory. As shown in Fig. 2(a), the coupling coefficient between the bus waveguide and the nanodisk is denoted by γ , β is the coupling coefficient between the nanodisk and the square ring. α , δ are the decay rates due to the internal loss of nanodisk and square ring, respectively. The amplitudes of the incoming and outing waves of disk are denoted by S +11 , S +12 , S −11 , and S −12 , as seen in Fig. 2(a). The temporal evolution of the amplitude A of the single nanodisk and B of the square ring. The characteristic equation can be given:

3. EIT-like in a square ring resonator coupled with a nanodisk

dA

56 57

88 89

26 27

87

The schematic illustration of EIT-like in metal–insulator–metal plasmonic system we proposed is shown in Fig. 2(a), it is composed of a nanodisk coupled with a square ring resonator. The inner and outer length of the square ring are l1 = 75 nm and l2 = 240 nm, respectively. More, in order to explain the subsequent analysis clearly, we set center length of square ring l = (l1 + l2 )/2. The width of bus waveguide is set to be d = 50 nm. The radii of nanodisk is r = 57 nm and the distance between the nanodisk and square ring is g = 11 nm. Fig. 2(b) shows the transmission spectra without and with the square ring resonator. Without the square

dt dB dt





91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121

= ( j ωt − α − β − γ ) A + j γ ( S +11 + S +12 ) + j β B

(3)

 = ( j ωn − β − δ) B + j β A

(4)

122 123 124 125 126

where ωt , ωn are the resonant frequency of the nanodisk and square ring, respectively. Due to energy conservation and the time reversal symmetry, the relationship of the incoming and outing waves in the bus waveguide can be denoted as:

S −11 = S +12 + j

90



127 128 129 130 131

γA

(5)

132

JID:PLA

AID:23488 /SCO Doctopic: Plasma and fluid physics

[m5G; v1.161; Prn:27/10/2015; 10:41] P.3 (1-6)

X. Li et al. / Physics Letters A ••• (••••) •••–•••

3

1

67

2

68

3

69

4

70

5

71

6

72

7

73

8

74

9

75

10

76

11

77

12

78

13

79

14

80

15

81

16

82

17

83

18

84

19

85

20

86

21

87

22

88

23

89

24

90

25

91

26

92

27

93

28

94

29

95

30

96

31

97

32

98

33

99

34

100

35

101

36

102

37 38 39

Fig. 2. (a) Schematic diagram of a nanodisk coupled with the square ring in MIM waveguide system; (b) Transmission spectra without (black curve) and with (red curve) a square resonator, the geometrical parameters of the structure are set as d = 50 nm, r = 57 nm, g = 11 nm, l1 = 75 nm, l2 = 240 nm; (c–d) The magnetic distribution without and with square ring for the incident wavelength of 825 nm (i.e. EIT-like peak wavelength).

40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57

S −12 = S +11 + j



+11

(6)

β

j (ω − ωt ) + α + β + j (ω−ω )+β+δ n β

j (ω − ωt ) + α + β + γ + j (ω−ω )+δ+β n

2   

(7)

From Eq. (7), we can quantitatively and clearly comprehend the EIT-like phenomenon. The peak transmission can be achieved at the resonance frequency ω = ωn = ωt with the value of T = α +β+β/(β+δ) 2 | α +β+ γ +β/(β+δ) | . It is obvious that the EIT-like transparency window is formed at resonance frequency ω = ωt as long as ωt and ωn are close to each other, and the transmission is derived as j (ω −ωn )+β+1 2 T = | j (ω −t ω )+β+ γ +1 | on condition that the intrinsic loss is igt n nored.

58 59 60

4. Transmission properties of the proposed structure with different parameters

61 62 63 64 65 66

104 105 106

γA

When the S +12 = 0, from Eqs. (3)–(6), the transmission coefficient T into the output port at steady state can be determined as:

    S −12 2     T =  = S

103

As we know, the transmission characteristics of the plasmonics waveguide system can be affected by the structure parameters. First, we calculated the transmission spectra for different coupling distance g between the nanodisk and square ring. Fig. 3(a) exhibits transmission spectra with different coupling gap g (9–14 nm)

at d = 50 nm, r = 57 nm, l1 = 80 nm, l2 = 250 nm. It is found that the transmission peak and the full width at half-maximum (FHWM) decrease with increasing distance g, so we can select g to get suitable transmission and FHWM. The length of the square ring resonator is also an important factor influencing the transmission properties. As a ring resonator, the standing wave will be excited in the square ring according to the theory of ring resonator [30,31]. The character equation can be given by

J n (kl1 ) J n (kl2 )



N n (kl1 ) N n (kl2 )

107 108 109 110 111 112 113 114 115

=0

where k = ε c (εr )1/2 ,

(8)

εr = n2neff is the effective relative permittiv-

ity, nneff is the effective refractive index of square ring. J n and J n stand for the first kind Bessel function with the order n and its derivative. N n and N n are the second Bessel function with order n and its derivative. Fig. 3(b) exhibits the transmission spectra with different l, it is obvious that the transmission peak of various l stays at a high level and the transparency peak performs a redshift with increasing l. In Fig. 3(c), the wavelength of peak spectra exhibits an approximately linear relation with l, which is in accord with the solution of Eq. (8). The transmission spectra with different refractive indices of nanodisk and square ring are investigated and other parameters are the same as those for Fig. 2(a). As shown in Fig. 4(a), the peak exhibits a redshift with the increase of the refractive index in

116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132

JID:PLA

AID:23488 /SCO Doctopic: Plasma and fluid physics

[m5G; v1.161; Prn:27/10/2015; 10:41] P.4 (1-6)

X. Li et al. / Physics Letters A ••• (••••) •••–•••

4

1

67

2

68

3

69

4

70

5

71

6

72

7

73

8

74

9

75

10

76

11

77

12

78

13

79

14

80

15

81

16

82

17

83

18

84

19

85

20

86

21

87

22

88

23

89

24

90

25

91

26

92

27

93

28

94

29

95

30

96

31

97

32 33 34

Fig. 3. (a) Transmission spectra with different g of the square ring resonator, the other parameters d = 50 nm, r = 57 nm, g = 11 nm, the thickness of the square ring waveguide is fixed as 170 nm. (b) EIT-like peak wavelength with different l. (c) Transmission spectra with different distances l, the other geometrical parameters are the same as Fig. 2(a).

35 36 37 38 39 40 41 42 43

46 47 48 49 50 51

FOM =

S

=

δλ/δn

λ

(9)

FWHM where S = δλ/δn is the refractive index sensitivity (i.e. spectral shifts per refractive index). λ is the resonance line width as the full FWHW centered at the resonance wavelength λ. The refractive index sensitivity and FOM of the new system is as high as 812 nm/RIU and 17.6, respectively.

56 57

60 61 62 63 64 65 66

ng =

c

υg

c

c dϕ (ω)

L

L

= τg = .



102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 118 119

As we know, EIT-like phenomenon is always accompanied by slow-light effect. While the slow-light effect can be depicted as the group index n g , which is denoted by [29]:

58 59

100

117

5. Slow light effect of this new structure

54 55

Fig. 5(a) shows the phase shift from the optical source to the detector. We can see that the phase spectrum has two frequency segments due to anomalous and normal phase dispersion and the phase slope is steepest at the location of the EIT-like peak. Examining Fig. 5(a), we note that the phase spectrum has glitch of phase. Based on the careful analysis, we think the phase dithering results from reduplicative reflection of light between the narrow square ring. It is obvious that nair < nsilver , the reflective light has a dispersion of phase comparing with incident light. As shown in the red area of Fig. 5(b), there is an optical delay as high as 0.14 ps in the new system at the EIT-like peak. Based on the above analysis of phase glitch and Eq. (10), it is easy to understand the dithering of group delay and group index. Moreover, the corresponding group index reach 41.1 at the EIT-like peak, as shown in the red area of Fig. 5(c). 6. Conclusion

52 53

99 101

nanodisk and square ring. In Fig. 4(b), it is obvious that the EIT-like resonant peak nearly has a linear relation with refractive index, the transmission peaks of transmittances are about 80% at 835, 851, 867, 884, and 900 nm when n = 1.00, 1.02, 1.04, 1.06, and 1.08, respectively, so we can use the new system as a sensor when the square ring and nanodisk is filled with changed dielectric. For sensing application, figure of merit (FOM) is usually applied to further evaluate the sensing performance as the following [32],

44 45

98

(10)

where c is the velocity of light in the free space, υ g is the group velocity in the plasmonic waveguide systems, τ g and ϕ (ω) represent the optical delay time and transmission phase shift, respectively, L stands for length of the new structure. The slow-light effect is investigated numerically in this system with d = 50 nm, r = 57 nm, g = 11 nm, l1 = 250 nm, l2 = 80 nm.

A new kind of analogies of EIT in metal–insulator–metal plasmonic waveguide consisting of a square ring resonator coupled with a nanodisk is proposed and numerically investigated. The transparency window can be easily tuned by changing the geometrical parameters of the structure and the insulator filled in resonator. The CMT-based transmission line theory is well in accordance with the FDTD results, what’s more, we find the EIT-like resonant peak still has a linear relationship with refractive index and the refractive index sensitivity and FOM of the new system is as high as 812 nm/RIU and 17.6, respectively. While group index reach 41.1 at the EIT-like peak with the transmission is close to 80% at d = 50 nm, r = 57 nm, g = 11 nm, l1 = 250 nm, l2 = 80 nm. These advantages of our analogy of electromagnetically induced

120 121 122 123 124 125 126 127 128 129 130 131 132

JID:PLA

AID:23488 /SCO Doctopic: Plasma and fluid physics

[m5G; v1.161; Prn:27/10/2015; 10:41] P.5 (1-6)

X. Li et al. / Physics Letters A ••• (••••) •••–•••

5

1

67

2

68

3

69

4

70

5

71

6

72

7

73

8

74

9

75

10

76

11

77

12

78

13

79

14

80

15

81

16

82

17 18

Fig. 4. (a) Transmission spectra with different refractive index of a nanodisk and a square ring. The other parameters are d = 50 nm, r = 57 nm, g = 11 nm, l1 = 250 nm, l2 = 80 nm. (b) EIT-like peak wavelength with different refractive index.

83 84

19

85

20

86

21

87

22

88

23

89

24

90

25

91

26

92

27

93

28

94

29

95

30

96

31

97

32

98

33

99

34

100

35

101

36

102

37

103

38

104

39

105

40

106

41

107

42

108

43

109

44

110

45

111

46

112

47

113

48

114

49

115

50

116

51

117

52

118

53 54 55 56 57 58 59 60 61 62 63 64

Fig. 5. Simulated transmission phase shift (a), optical delay time (b), group index (c) with d = 50 nm, r = 57 nm, g = 11 nm, l1 = 250 nm, l2 = 80 nm.

transparency based on plasmonic waveguide resonator system provide a potential application in highly integrated optical circuits and signal processing. Acknowledgements This work was supported by the National Natural Science Foundation of China (NSFC) under Grants Nos. 61275059 and 11374107. References

65 66

[1] C.L. Garrido Alzar, M.A.G. Martinez, P. Nussenzveig, Am. J. Phys. 70 (2002) 37.

J. Eberly, Opt. Express 2 (1998) 173. L. Maleki, A.B. Matsko, A.A. Savchenkov, V.S. Ilchenko, Opt. Lett. 29 (2004) 626. S.E. Harris, Phys. Today 50 (1997) 36. M.D. Lukin, S.F. Yelin, M. Fleischhauer, Phys. Rev. Lett. 84 (2000) 4232. V. Yannopapas, E. Paspalakis, N.V. Vitanov, Phys. Rev. B 80 (2009) 035104. X. Yang, M. Yu, D.-L. Kwong, C.W. Wong, Phys. Rev. Lett. 102 (2009) 173902. D.A. Ketzaki, O. Tsilipakos, T.V. Yioultsis, E.E. Kriezis, J. Appl. Phys. 114 (2013) 113107. [9] H. Lu, X. Liu, Y. Gong, D. Mao, L. Wang, Appl. Opt. 50 (2011) 1307. [10] A. Ourir, B. Gallas, L. Becerra, J. de Rosny, P. Dahoo, J. Biophotonics 2 (2015) 308. [11] H.-m. Li, S.-b. Liu, S.-y. Liu, S.-y. Wang, H.-f. Zhang, B.-r. Bian, X.-k. Kong, Appl. Phys. Lett. 106 (2015) 114101. [2] [3] [4] [5] [6] [7] [8]

119 120 121 122 123 124 125 126 127 128 129 130 131 132

JID:PLA

6

1 2 3 4 5 6 7 8 9 10 11

AID:23488 /SCO Doctopic: Plasma and fluid physics

[m5G; v1.161; Prn:27/10/2015; 10:41] P.6 (1-6)

X. Li et al. / Physics Letters A ••• (••••) •••–•••

[12] Q. Xu, S. Sandhu, M.L. Povinelli, J. Shakya, S. Fan, M. Lipson, Phys. Rev. Lett. 96 (2006) 123901. [13] L. Zhou, T. Ye, J. Chen, Opt. Lett. 36 (2011) 13. [14] W.L. Barnes, A. Dereux, T.W. Ebbesen, Nature 424 (2003) 824. [15] C. Genet, T.W. Ebbesen, Nature 445 (2007) 39. [16] H.J. Lezec, A. Degiron, E. Devaux, R.A. Linke, L. Martin-Moreno, F.J. Garcia-Vidal, T.W. Ebbesen, Science 297 (2002) 820. [17] Z. Han, S.I. Bozhevolnyi, Opt. Express 19 (2011) 3251. [18] J. Chen, Z. Li, Y. Zou, Z. Deng, J. Xiao, Q. Gong, Plasmonics 8 (2013) 1627. [19] C. Wu, N. Arju, G. Kelp, J.A. Fan, J. Dominguez, E. Gonzales, E. Tutuc, I. Brener, G. Shvets, Nat. Commun. 5 (2014). [20] M. Hu, F. Wang, R. Liang, S. Zhou, L. Xiao, Phys. Lett. A 379 (2015) 581. [21] Y. Cui, C. Zeng, Opt. Commun. 297 (2013) 190.

[22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32]

G. Wang, H. Lu, X. Liu, Opt. Express 20 (2012) 20902. H. Lu, X. Liu, G. Wang, D. Mao, Nanotechnology 23 (2012) 444003. H. Lu, X. Liu, D. Mao, Y. Gong, G. Wang, Opt. Lett. 36 (2011) 3233. J.A. Dionne, L.A. Sweatlock, H.A. Atwater, A. Polman, Phys. Rev. B 73 (2006) 035407. ´ A.B. Djurišic, ´ J.M. Elazar, M.L. Majewski, Appl. Opt. 37 (1998) A.D. Rakic, 5271. Q. Li, F. Liu, Z. Zhang, M. Qiu, Y. Su, J. Lightwave Technol. 26 (2008) 3744. Y. Akahane, T. Asano, B.-S. Song, S. Noda, Opt. Express 13 (2005) 1202. H. Lu, X. Liu, D. Mao, Phys. Rev. A 85 (2012) 053803. J. Wang, W.B. Lu, X.B. Li, Z.H. Ni, T. Qiu, J. Phys. D, Appl. Phys. 47 (2014) 135106. I. Wolff, N. Knoppik, Electron. Lett. 7 (1971) 779. J. Wang, C. Fan, J. He, P. Ding, E. Liang, Q. Xue, Opt. Express 21 (2013) 2236.

67 68 69 70 71 72 73 74 75 76 77

12

78

13

79

14

80

15

81

16

82

17

83

18

84

19

85

20

86

21

87

22

88

23

89

24

90

25

91

26

92

27

93

28

94

29

95

30

96

31

97

32

98

33

99

34

100

35

101

36

102

37

103

38

104

39

105

40

106

41

107

42

108

43

109

44

110

45

111

46

112

47

113

48

114

49

115

50

116

51

117

52

118

53

119

54

120

55

121

56

122

57

123

58

124

59

125

60

126

61

127

62

128

63

129

64

130

65

131

66

132