investigation on the special smith purcell radiation from a nano scale rectangular metallic grating

Investigation on the special Smith-Purcell radiationfrom a nano-scale rectangular metallic grating Weiwei Li, Weihao Liu,aand Qika Jia National Synchrotron Radiation Laboratory, Universi

Investigation on the special Smith-Purcell radiation from a nano-scale rectangular metallic grating

Weiwei Li, Weihao Liu, and Qika Jia

Citation: AIP Advances 6, 035202 (2016); doi: 10.1063/1.4943502

View online: http://dx.doi.org/10.1063/1.4943502

View Table of Contents: http://aip.scitation.org/toc/adv/6/3

Published by the American Institute of Physics

Articles you may be interested in

Theory of the special Smith-Purcell radiation from a rectangular grating

AIP Advances 5, 127135127135 (2015); 10.1063/1.4939538

Investigation on the special Smith-Purcell radiation

from a nano-scale rectangular metallic grating

Weiwei Li, Weihao Liu,aand Qika Jia

National Synchrotron Radiation Laboratory, University of Science and Technology of China,

Hefei, Anhui, 230029, China

(Received 23 January 2016; accepted 24 February 2016; published online 4 March 2016)

The special Smith-Purcell radiation (S-SPR), which is from the radiating eigen modes of a grating, has remarkable higher intensity than the ordinary Smith-Purcell radiation Yet in previous studies, the gratings were treated as perfect conductor without considering the surface plasmon polaritons (SPPs) which are of signif-icance for the nano-scale gratings especially in the optical region In present paper, the rigorous theoretical investigations on the S-SPR from a nano-grating with SPPs taken into consideration are carried out The dispersion relations and radiation characteristics are obtained, and the results are verified by simulations According to the analyses, the tunable light radiation can be achieved by the S-SPR from a nano-grating, which offers a new prospect for developing the nano-scale light sources C 2016 Author(s) All article content, except where other-wise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).[http://dx.doi.org/10.1063/1.4943502]

I INTRODUCTION

Since the first day of its observation, Smith-Purcell radiation (SPR), occurring as electron beam passes over a periodic surface, has been one of the most attractive research topics because

of its tremendous applications.17 The most promising application in the present time may be its potential in generating terahertz (THz) wave radiation8,9 which can not be easily obtained by other means.10–16 The metallic rectangular gratings are commonly used to get the SPR in these applications.17,18The representative characteristic of the SPR is its famous dispersion relation:

λ = −L

where λ is radiation wavelength, θ indicates radiation direction, L is the structural period, β is the ratio of beam velocity to light speed in vacuum, and n is a negative integer that indicates the harmonic order

By reducing the groove width of the rectangular grating, Ref.20recently uncovered an inter-esting special kind of Smith-Purcell radiation (S-SPR) which is a monochromaticity radiation with much higher intensity in specified direction The theoretical analyses of the S-SPR were carried out

in Ref 21, which demonstrated that the S-SPR is exactly from the radiating eigen modes of the grating The remarkable advantages in intensity and tunability make the S-SPR a promising way for tunable wave generation and beam diagnostic

In Ref.21, the metal was treated as perfect conductor, namely, the oscillating fields in metal had been ignored Based on the rigorous theory of metal, this simplification is not so accurate, especially for a nano-scale structure operating in the optical frequency region The more general properties

of metal are governed by the Drude-Sommerfeld’s free-electron-gas model, according to which the metal should be described by the following dielectric function22,23:

a Electronic mail: liuwhao@ustc.edu.cn

035202-2 Li, Liu, and Jia AIP Advances 6, 035202 (2016)

εm(ω) = ε∞− ω2

p

where ε∞is a constant representing the influence of the metal atoms, ωpis the effective bulk plasma frequency of the free electrons in the metal, and γ indicates the electrons’ collision frequency From

Eq (2) one can see that the metal is actually a dispersive (frequency dependent) dielectric medium Based on Drude’s model, the evanescent electromagnetic mode, terminologically the surface plas-mon polaritons (SPPs), will present on the metal surface in the optical region The physical base of SPPs are the collective oscillations of the conduction electrons in metal In past years, SPPs have attracted extensive interests for their great promising applications, especially with the advancement

of the nano-scale fabrication techniques.24 – 26When the SPPs are considered in the metallic grating, the electromagnetic properties as well as the radiation properties of the grating will greatly be changed

The SPR from the metallic surface with SPPs taken into consideration was preliminarily analyzed in Ref.27, which showed that the radiation intensity can be enhanced in specific direc-tions Years later, Ref 28 analytically predicted a large emission rate from a chain of nano aluminum spheres via the SPR effect Ref.29simulated the SPR from a rectangular metallic grat-ing, and it showed the radiations were composed by three primary components: the ordinary SPR, the original SPPs, and the so-called mimic-SPPs Most recently, Refs.30and31demonstrated that the SPPs on a metallic nanowire-array can be transformed into radiation via an exotic Smith-Purcell

effect In mechanism, it was the collective radiation from the SPPs modes in each unit of the array The radiation power density can be significantly higher than the ordinary SPR All those results illustrate the copious interesting phenomena when the SPPs are considered

In present paper, we will carry out the rigorous theoretical investigation on the S-SPR from the rectangular metallic grating with the SPPs taken into consideration More interesting results will be obtained, which indicate new potential applications The paper is organized as follows: Section II

presents the derivations of the theoretical formulas, Section IIIshows the calculation results and simulation verifications, and SectionIVgives the conclusion

II FORMULATION OF THE PROBLEM

The schematic diagram of the problem that we are going to study is shown in Fig.1 An elect-ron beam passes over a metallic grating with rectangular groves When the SPPs are considered, the fields will penetrate into the metallic grating, which means the whole space will be filled with the oscillation fields This will greatly increase the complexity of the problem as will be shown

To handle the problem, we apply the mode matching method and divide the whole space into four regions according to the boundaries as shown in the figure The fields in different regions can be obtained by solving the homogeneous or non-homogeneous Maxwell equations

In each region, all electromagnetic field components can be expressed as functions of the z-directional electric field component Ez, which should be obtained first In the upper half space

FIG 1 Schematic diagram of the problem.

(region I), Ezcomponent satisfies the following non-homogeneous equation in the Cartesian coordi-nate system32:

∂2Ez

∂x2 +∂2Ez

∂z2 − 1

c2

∂2Ez

∂t2 =ε1

0

∂ ρ

∂z +µ0

∂ Jz

where c is the light speed in vacuum, ε0and µ0are respectively the permittivity and permeability

of the vacuum, ρ and Jzrespectively denote the charge density and current density For the line electron beam with quantity of q and velocity of v, they can be expressed as:

where x0denotes the x-directional position of the electron beam as shown in the figure

With the help of Fourier transformation method, Eq (3) can be solved The solutions are composed by two parts: the special solution which represents the incident wave from the electron beam and the general solution signifying the diffraction wave from the grating To satisfy the peri-odic boundary conditions of the grating, the diffraction wave should be expanded into infinite space harmonics Then the total field in region I can be expressed as:





















Eω,zI = EI, i

ω,z+ EI,r ω,z

4πωε0

e− j kc |x−x 0 |− jk z z+

∞



n =−∞

An

j kI

x n

ωε0

e−kIx n xe− j kz n z

Hω, yI = HI, i

ω, y+ HI,r

ω, y

= s q 4πe

− j kc|x−x 0 |− jk z z+

∞



n =−∞

Ane−kIx n xe− j kz n z

(6)

where k2= k2

0− k2= (ω

c)2−(ωv)2, kI

x n 2

= k2

z n− k2

0, kz n= kz+2nπ

L , s= −1 for x > x0 and s= 1 for x < x0 In Eq (6), the Excomponent was omitted for simplification as it will not be used in the mode matching method Note that if kz n> k0, kIx nis real, namely, the diffraction waves are evanes-cent in x direction and can not radiate into space; while if kz n< k0, kIx nis imaginary, which means the diffraction waves can propagate in x direction It’s easy to find that only the negative harmonic waves (n < 0) can radiate into upper space The radiation wavelength and direction should satisfy the SPR relation of Eq (1)

For the fields in regions II and IV, the homogeneous equation of Eq (3) should be solved according to their own boundary conditions Here we only consider the fundamental modes, which

is reasonable since the gap width is much less than the wavelength Then the fields in regions II and

IV can respectively be expressed as:





















Eω,zII = kIIx

ωεII

rε0

(

B1ej kIIx x− C1e− j kIIx x) cos kIIzz

cos kII

zd/2

Eω, xII = − j kzII

ωεII

rε0

(

B1ej kIIx x+ C1e− j kIIx x) sin kII

zz
cos kII

zd/2

Hω, yII =(B1ej kIIx x+ C1e− j kIIx x) cos kII

zz
cos kII

zd/2

and





















Eω,zIV = kIVx

ωεIV

r ε0

(

B2ej kIVx x− C2e− j kIVx x)

*

,

ekIVz(d/2−|z|)+ ekIVz(|z|−L+d/2)

1+ ekIVz(d−L) +

-Eω, xIV = j k

IV

z sign(z)

−ωεIV

r ε0

(

B2ej kIVx x+ C2e− j kIVx x)

*

,

ekIVz(d/2−|z|)− ekIVz(|z|−L+d/2)

1+ ekIVz(d−L) +

-Hω, yIV =(B2ej kIVx x+ C2e− j kIVx x)

*

ekIVz (d/2−|z|)+ ek IV

z (|z|−L+d/2)

1+ ekIVz(d−L)

+

035202-4 Li, Liu, and Jia AIP Advances 6, 035202 (2016)

In above equations, εII

rk2

0= kII2

x +kII2

z , εIV

r k2

0= kIV2

x +kIV2

z , εII

ris the dielectric constant of the medium filled in region II (in following calculations its value is unit), and εIV

r = εmis the dielectric function

of metal expressed by Eq (2)

Following the case of region I, fields in region III can be written as:













Hω, yIII =

n

DnekIIIx n x− j kz nz

Eω,zIII = −

n

Dn

j kIII

x n

εIII

rε0ωe

k III

x n x− j kz nz (9)

in which εIII

r k2

0= kIII2

x n+kIII2

zn In all previous equations, the Ans, B, C, and Dnsare the coefficients to

be determined by the boundary conditions

Now we consider the boundary conditions At the boundary of adjacent gaps z= ±d/2, the continuity conditions for the tangential fields in regions II and IV can be expressed as:









HyII

z =±d/2= HIV

y

z =±d/2

EIIxz=±d/2= EIV

x



z =±d/2

Substituting Eqs (7), (8) into Eqs (10) we can get

kzII

εII r

tan kzIId/2
= kIVz

εIV r

1 − ekIVz(d−L)

1+ ekIVz(d−L), (11) which defines the relation between frequency and kI I

x So it is actually the dispersion equation concerning the waves propagating in the x direction along the interfaces of the metal and the gaps.33

At the boundaries of x=0 and x=-h, the following mode matching conditions can be used34:













Eω,zI, i 

x =0+ EI,r ω,zx=0=







Eω,zII 

x =0, |z| < d/2

Eω,zIVx=0, d/2 < |z| < L/2

Hω, yI, i

x =0+ HI,r

ω, y

x =0= HII

ω, y

x =0, |z| < d/2,

(12)

and













Eω,zIII 

x =−h=







Eω,zII x=−h, |z| < d/2

Eω,zIV x=−h, d/2 < |z| < L/2

Hω, yIII

x =−h= HII

ω, y

x =−h, |z| < d/2

(13)

Substituting Eqs (6)-(9) into Eqs (12) and (13) and doing complicate but straightforward arithmet-ical operation, coefficients (Ans) of the diffraction waves can be solved The radiation power can then be obtained Following Ref.35, the radiation power per solid angle is

dIn

2L2ε0

n2sin2θ (1/ β − cos2θ)3|Rn|2exp

(

−2x0

λe

)

where Rnis the radiation factor of the grating given by:

|Rn|2=

Anej kc x0 q

4π

2

Letting the incident waves vanish in Eq (12), we can get the dispersion equation of the structure:

*

,

n

M1, nN1, n+

-*

,

jtan kIIxs


n

M2, nN1, n+ N2+

n

M2, nN1, n− tan kIIxs
N2+N2

where M1, n= Mn/ jkI

x nL
, N1, n= d sin c(k z n d

2

) ,

cos kII

zd/2

 d/2

−d/2

and

FIG 2 Dispersion curves (a) and radiation spectra (b) versus the grove depth h, here L =300 nm, d=20 nm Dispersion curves (c) and radiation spectra (d) versus the grove width d, here L =300 nm, h=30 nm The shaded regions of the figure indicate the surface modes, and curves above the light lines signify the radiating modes.

035202-6 Li, Liu, and Jia AIP Advances 6, 035202 (2016)

εII

r cos kII

zd/2

d/2



−d/2

ej kz n zcos kIIzd/2
dz

εIV r

(

1+ ek IV

z (d−L))

L/2



d/2

(

ekIVz (d/2−|z|)+ ek IV

z (|z|−L+d/2))cos(kz nz) dz

It is a complex transcendental equation which should be solved by the numerical method The detailed calculations will be given in the next section For the case that the SPPs are not considered,

we only need to let εIII

r, εIV

III DISPERSION RELATIONS AND RADIATION PROPERTIES

The electromagnetic dispersion properties of the grating essentially determines the radiation characteristics According to Ref.21, there are two kinds of eigen modes that can be supported by the grating: the surface eigen modes and radiating eigen modes To get the S-SPR which has much higher intensity than the ordinary SPR, the radiating eigen modes should be effectively excited by the electron beam whereas the surface eigen modes should not be excited It is essentially different from the ordinary SPR, in which the electron beam principally interacts with the surface wave.18 , 19

Here we will show that when the SPPs are taken into consideration for a nano-grating, both the dispersion and radiation properties of the grating will be changed

Based on the equations derived in Section II, the dispersion curves and radiation spectra of the grating can be calculated Figure2shows the calculation results for the gratings with different structure parameters Here the dielectric parameters of the metal are set as: ωp= 1.39 × 1016r ad/s,

γ = 3.2 × 1013H z, and ε∞= 4.5, indicating that the silver grating is considered The exciting beam energy is 100 keV One can observe that the radiating eigen modes are effectively excited, and the S-SPRs with obvious peak frequencies are generated The radiation frequencies change gradually with the width ‘d’ or depth ‘h’ of the groves (increase with d while decrease with h) The radiation from the first order radiating modes can be tuned from 500 THz to 700 THz And frequency of the second order radiating modes can extend to be higher than 900 THz In other words, the frequency can be tuned from infrared to ultraviolet region by adjusting the structure parameters For compar-isons, the cases without considering the SPPs are also calculated The results show that neither surface modes nor radiating modes can be excited by the electron beam, and the radiation intensity

is significantly lower than that with SPPs as shown in the figure This is because that there is no wave can be supported by the grating in this frequency region when the SPPs are not considered In other words, under the influences of the SPPs, both the dispersion curves and the radiation spectra are significantly changed

Then we exam the electrical tunability of the radiation According to the dispersion curves obtained above, we choose the structure parameters to be L=300 nm, h=20 nm, d=20 nm, and the beam energy changes from 60 keV to 100 keV The calculated radiation spectra are shown in Fig.3

FIG 3 Radiation spectra versus the beam voltage, here L =300 nm, h=20 nm, and d=20 nm.

FIG 4 Comparisons for the theoretical and simulation results of the radiation frequency depending on (a) h, (b) d, and (c) beam energy Other structure parameters are the same as that given in Fig 2

One can observe that the radiation frequency increases with the beam energy, and the tuning range is from 600 THz to 700 THz

To verify the theoretical results, we perform the simulations by applying the fully electromag-netic particle-in-cell code.36 , 37The simulation model and conditions follow that given in Refs 20

and 30, in which the simulation obtained radiation snapshots and spectra are presented The comparisons for the theoretical and simulation results of the radiation frequencies are shown in Fig.4 Here only the first order radiating modes are presented One can observe that their results are reasonably agreed with each other

IV CONCLUSION

The investigation on the special Smith-Purcell radiation from the rectangular metallic nano grating with surface plasmon polaritons taken into consideration was carried out by applying the fundamental electromagnetic theory and numerical simulations By adjusting the structure and beam parameters, the tunable light radiations are obtained, which offers a new prospect for develop-ing the nano-scale light sources

This work is supported by Natural Science Foundation of China (Grant No 61471332) and Anhui Provincial Natural Science Foundation (Grant No 1508085QF113)

1 S J Smith and E M Purcell, Phys Rev 92, 1069-1070 (1953).

2 S E Korbly, A S Kesar, J R Sirigiri, and R J Temkin, Phys Rev Lett 94, 054803 (2005).

3 D Li, M Hangyo, Z Yang, M R Asakawa, S Miyamoto, Y Tsunawaki, K Takano, and K Imasaki, Nucl Instrum Methods Phys Res A 637(1), 135-137 (2011).

4 J Urata, M Goldstein, M F Kimmitt, A Naumov, C Platt, and J E Walsh, Phys Rev Lett 80(3), 516-519 (1998).

5 T Simpei, I Koji, and S Eiichi, Opt Express 15(24), 16222-16229 (2007).

6 H L Andrews, C A Brau, and J D Jarvis, Phys Rev Special Topics 8, 110702 (2005).

7 J Gardelle, P Modin, and J T Donohue, in Plasma Sciences (ICOPS) held with 2014 IEEE International Conference on High-Power Particle Beams (BEAMS), 2014 IEEE 41st International Conference on (IEEE, 2014), p 1.

8 O Akiko and M Tatsunosuke, Opt Express 22(14), 17490-17496 (2014).

9 K Tantiwanichapan, X Wang, A K Swan et al., Appl Phys Let 105, 24 (2014).

035202-8 Li, Liu, and Jia AIP Advances 6, 035202 (2016)

10 S Son, S J Moon, and J Y Park, Opt Let 38(22), 4578-4580 (2013).

11 C Seco-Martorell, V Lopez-Dominguez, G Arauz-Garofalo, A Redo-Sanchez, J Palacios, and J Tejada, Opt Exp 21(15), 17800-17805 (2013).

12 R Eichholz, H Richter, M Wienold, L Schrottke, R Hey, H T Grahn, and H.-W Hubers, Opt Exp 21(26), 32199-32206 (2013).

13 H He and X.-C Zhang, Appl Phys Lett 100, 061105 (2012).

14 Q Qin, B S Williams, S Kumar, J L Reno, and Q Hu, Nat Photonics 3, 732 (2009).

15 W Liu and Z Xu, J Appl Phys 115, 014503 (2014).

16 W Liu, S Gong, Y Zhang, J Zhou, P Zhang, and S liu, J Appl Phys 111, 063107 (2012).

17 M Cao, W Liu, Y Wang, and K Li, Phys Plasmas 21, 023116 (2014).

18 M Cao, W Liu, Y Wang, and K Li, J Appl Phys 116, 103104 (2014).

19 P V D Berg, JOSA 64(3), 325-328 (1974).

20 W Liu and Z Xu, New J Phys 16, 073006 (2014).

21 W Liu, W Li, Z He, and Q Jia, AIP adv 5, 127135 (2015).

22 P B Johnson and R W Christy, Phys Rev B 6, 4370 (1972).

23 S A Maier and H A Atwater, J Appl Phys 98, 011101 (2005).

24 J M Romo-Herrera, R A Alvarez-Puebla, and L M Liz-Marzán, Nanoscale 3(4), 1304-1315 (2011).

25 S I Bozhevolnyi1, V S Volkov, E Devaux, J Y Laluet, and T W Ebbesen, Nature 440, 508C511 (2006).

26 S Thongrattanasiri and F J García de Abajo, Phys Rev Lett 110, 187401 (2013).

27 S L Chuang and J A Kong, J Opt Soc Am A 1, 672 (1984).

28 F J García de Abajo, Phys Rev E 61, 5743 (2000).

29 S Taga, K Inafune, and E Sano, Opt Express 15(24), 16222-16229 (2007).

30 W Liu, Opt Lett 40(20), 4579-4582 (2015).

31 W Liu, Photonics J IEEE 7(6), 1-8 (2015).

32 P V D Berg, JOSA 63(12), 1588-1597 (1973).

33 W Liu, R Zhong, J Zhou, Y Zhang, M Hu, and S Liu, J Phys D: Appl Phys 45, 205104 (2012).

34 K Zhang and D Li, Electromagnetic Theory in Microwave and Optoelectronics (2008), pp 426-430.

35 A Gover, P Dvorkis, and U Elisha, JOSA 1(5), 723-728 (1984).

36 B Goplen, L Ludeking, D Smith, and G Warren, Computer Phys Commun 87, 54-86 (1995).

37 W Liu, Opt Lett 40(17), 3974 (2015).

IV CONCLUSION

The investigation on the special Smith- Purcell radiation from the rectangular. .. rectangular metallic nano grating with surface plasmon polaritons taken into consideration was carried out by applying the fundamental electromagnetic theory and numerical simulations By adjusting... both the dispersion and radiation properties of the grating will be changed

Based on the equations derived in Section II, the dispersion curves and radiation spectra of the grating can be