Supplementary material for " Coherent and Tunable Terahertz Radiation from Graphene Surface Plasmon Polarirons Excited by Cyclotron Electron Beam "

Similar documents
School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

Collection of Formulas

Electric Potential. and Equipotentials

U>, and is negative. Electric Potential Energy

General Physics II. number of field lines/area. for whole surface: for continuous surface is a whole surface

FI 2201 Electromagnetism

Friedmannien equations

On the Eötvös effect

Prof. Anchordoqui Problems set # 12 Physics 169 May 12, 2015

6. Numbers. The line of numbers: Important subsets of IR:

NEGATIVE REFRACTIVE INDEX SEMICONDUCTOR CIRCULAR WAVEGUIDE FIBER

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 20

CHAPTER 18: ELECTRIC CHARGE AND ELECTRIC FIELD

STUDY OF THE UNIFORM MAGNETIC FIELD DOMAINS (3D) IN THE CASE OF THE HELMHOLTZ COILS

Electric Field F E. q Q R Q. ˆ 4 r r - - Electric field intensity depends on the medium! origin

6. Gravitation. 6.1 Newton's law of Gravitation

Chapter 28 Sources of Magnetic Field

Week 10: DTMC Applications Ranking Web Pages & Slotted ALOHA. Network Performance 10-1

Fourier-Bessel Expansions with Arbitrary Radial Boundaries

Algebra Based Physics. Gravitational Force. PSI Honors universal gravitation presentation Update Fall 2016.notebookNovember 10, 2016

3.1 Magnetic Fields. Oersted and Ampere

The Formulas of Vector Calculus John Cullinan

ELECTRO - MAGNETIC INDUCTION

ELECTROSTATICS. 4πε0. E dr. The electric field is along the direction where the potential decreases at the maximum rate. 5. Electric Potential Energy:

Effect of Heat Generation on Quasi- Static Thermal Stresses in a Solid Sphere

Winter 2004 OSU Sources of Magnetic Fields 1 Chapter 32

Answers to test yourself questions

Picking Coordinate Axes

Lecture 11: Potential Gradient and Capacitor Review:

Electricity & Magnetism Lecture 6: Electric Potential

This immediately suggests an inverse-square law for a "piece" of current along the line.

Qualitative Analysis for Solutions of a Class of. Nonlinear Ordinary Differential Equations

NONLOCAL HOMOGENIZATION MODEL FOR THE ANALYSIS OF ABSORBING PROPERTIES OF MUSHROOM STRUCTURES WITH GRAPHENE PATCHES AT MICROWAVES

4.2 Boussinesq s Theory. Contents

(A) 6.32 (B) 9.49 (C) (D) (E) 18.97

Classical Electrodynamics

DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING FLUID MECHANICS III Solutions to Problem Sheet 3

Physics 11b Lecture #11

Chapter 6 Frequency Response & System Concepts

Production Mechanism of Quark Gluon Plasma in Heavy Ion Collision. Ambar Jain And V.Ravishankar

MAGNETIC EFFECT OF CURRENT & MAGNETISM

Ch 26 - Capacitance! What s Next! Review! Lab this week!

Chapter 2. Review of Newton's Laws, Units and Dimensions, and Basic Physics

ρ θ φ δ δ θ δ φ δ φ π δ φ π δ φ π

Two dimensional polar coordinate system in airy stress functions

Scientific Computing & Modelling NV, Vrije Universiteit, Theoretical Chemistry, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands c

Important design issues and engineering applications of SDOF system Frequency response Functions

Prof. Dr. Yong-Su Na (32-206, Tel )

COMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS

PX3008 Problem Sheet 1

CHAPTER 2 ELECTROSTATIC POTENTIAL

9.4 The response of equilibrium to temperature (continued)

RELATIVE KINEMATICS. q 2 R 12. u 1 O 2 S 2 S 1. r 1 O 1. Figure 1

Chapter 25: Current, Resistance and Electromotive Force. Charge carrier motion in a conductor in two parts

Physics 1502: Lecture 2 Today s Agenda

Chapter 4 Kinematics in Two Dimensions

Physics 604 Problem Set 1 Due Sept 16, 2010

Rydberg-Rydberg Interactions

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

DYNAMICS. Kinetics of Particles: Newton s Second Law VECTOR MECHANICS FOR ENGINEERS: Ninth Edition CHAPTER. Ferdinand P. Beer E. Russell Johnston, Jr.

A COMPARISON OF MEMBRANE SHELL THEORIES OF HYBRID ANISOTROPIC MATERIALS ABSTRACT

A Permanent Magnet Device for Measuring the Coercive Force

( ) ( ) ( ) ( ) ( ) # B x ( ˆ i ) ( ) # B y ( ˆ j ) ( ) # B y ("ˆ ( ) ( ) ( (( ) # ("ˆ ( ) ( ) ( ) # B ˆ z ( k )

(a) Counter-Clockwise (b) Clockwise ()N (c) No rotation (d) Not enough information

RELATIONSHIP BETWEEN DESIGN RESPONSE SPECTRA FOR RARE AND FREQUENT EARTHQUAKE LEVELS

Chapter 31 Faraday s Law

ANTENNAS. Vector and Scalar Potentials. Maxwell's Equations. D = εe. For a linear, homogeneous, isotropic medium µ and ε are contant.

[ ] = jω µ [ + jω ε E r

Supporting Information Wedge Dyakonov Waves and Dyakonov Plasmons in Topological Insulator Bi 2 Se 3 Probed by Electron Beams

Chapter 25: Current, Resistance and Electromotive Force. ~10-4 m/s Typical speeds ~ 10 6 m/s

International Journal of Pure and Applied Sciences and Technology

igid nd non-leky two-comptment building. Yu et l [8] developed non-line govening equtions by consideing the effect of bckgound lekge. Howeve, thee e n

Plane Wave Expansion Method (PWEM)

SPA7010U/SPA7010P: THE GALAXY. Solutions for Coursework 1. Questions distributed on: 25 January 2018.

Introduction to Arrays

Comparative Studies of Law of Gravity and General Relativity. No.1 of Comparative Physics Series Papers

D-STABLE ROBUST RELIABLE CONTROL FOR UNCERTAIN DELTA OPERATOR SYSTEMS

Chapter 6 Thermoelasticity

Available online at ScienceDirect. Procedia Engineering 91 (2014 ) 32 36

APEX CARE INSTITUTE FOR TRB, SLET AND NET IN PHYSICS ELECTROMAGNETIC THEORY- SLET / NET QUESTION BANK

Phase Velocities of Three-Dimensional and Axis-Symmetrical Elastic Waves in Isotropic Cylindrical Shell

Modelling of Low Velocity Impact Damage in Laminated Composites

Optimization. x = 22 corresponds to local maximum by second derivative test

Chapter 8. Ch.8, Potential flow

MAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS

UvA-VU Master Course: Advanced Solid State Physics

Superluminal Group Velocity of Electromagnetic Near-fields *

Previously. Extensions to backstepping controller designs. Tracking using backstepping Suppose we consider the general system

12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31,

Chapter 21: Electric Charge and Electric Field

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN

Dynamically Equivalent Systems. Dynamically Equivalent Systems. Dynamically Equivalent Systems. ME 201 Mechanics of Machines

Study of Electromagnetic Wave Propagation in Periodic Dielectric Structure; MathCAD Analysis

Physics 505 Fall 2005 Midterm Solutions. This midterm is a two hour open book, open notes exam. Do all three problems.

(1) It increases the break down potential of the surrounding medium so that more potential can be applied and hence more charge can be stored.

10 m, so the distance from the Sun to the Moon during a solar eclipse is. The mass of the Sun, Earth, and Moon are = =

The geometric construction of Ewald sphere and Bragg condition:

Radial geodesics in Schwarzschild spacetime

SURFACE TENSION. e-edge Education Classes 1 of 7 website: , ,

Transcription:

Suppleenty teil fo " Coheent nd Tunble Tehet Rdition fo Gphene Sufce Plson Polions Excited by Cycloton Electon Be " To Zho,, Sen Gong,, Min Hu,, Renbin Zhong,,Diwei Liu,,Xioxing Chen,, Ping hng,, Xinn Wng, 3, Cho Zhng, 4, Peiheng Wu, 3, Shenggng Liu *,. Tehet Resech Cente, School of Physicl Electonics, Univesity of Electonic Science nd Technology of Chin, Chengdu, Sichun, 654, Chin,. Coopetive nnovtion Cente of Tehet Science, Chengdu, Sichun, 654, Chin 3. School of Electonic Science nd Engineeing, Nnjing Univesity, Nnjing, Jingsu,, Chin 4. School of Physics nd nstitute fo Supeconducting nd Electonic Mteils, Univesity of Wollongong, New South Wles 5, Austli *liusg@uestc.edu.cn

Appendix : The electognetic fields poduced by the cycloton electon be Fig.. The schetic of cycloton electon be, the electon be oves t velocity u v e v e with cycloton tjectoy. The schetic of n cycloton electon be (CEB) is shown in Fig.., the electon be oves t velocity v ve ve with cycloton tjectoy. The vecto nd scl potentils A nd e used to deive the electognetic fields poduced by the CEB. The vecto nd scl potentils obey the following d'alebet equtions, given below in MKS unit, A t A t J (.) (.) With Loent guge: A t Then, the electognetic fields E nd H cn be obtined: A A E H t (.3) (.4) The position nd velocity vectos of the CEB cn be descibed in Ctesin co odinte syste. (.5) t e x t e y t e t e x y (.6) v t e v t e v t

Then, the chge nd cuent density cn be expessed. Whee q is the chge quntity, e q t q x x t y y t t (.7) J v t e v t e v t q x x t y y t t (.8) v is the coponent of velocity, v is the cycloton velocity of the be ound xis. n the cylindicl coodinte syst e, (.9) x t cos ct (.) y t sin ct v t (.) t Whee is the dius of the tjectoy pojection of cycloton electon be in X-Y plne, v c. Siilly, the obseved position of fields cn be expessed in the cylindicl coodinte syste. x cos (.) y sin (.3) (.4) Now the chge nd cuent density cn be expessed. cos c c cos t, t q e t q sin sin t vt cos c cos t J, t v t ev ev q sin sin ct vt (.5) (.6) The field cn be obtined by using the Fouie tnsfo, ssue wve fcto is in the fo jk e j t., 3 jk jt A t A k e e dk (.7) 3 jk e t J, t J k e dk (.8) 3

Subitting Eq. (.7) nd (.8) into Eq. (.) nd (.), the vecto A cn be obtined by Fouie tnsfotion in the wve vecto spce. A, t 3 jk e t qe v eve dk (.9) k k k Bsed on Loent guge, we cn get the scl potentil. A j (.) Bsed on the Eq. (.4), we obtin the coponents of the electic field. Appendix : The theoy of cicul cylindicl onolye gphene stuctue Fig.. Schetic of the cicul cylindicl onolye gphene stuctue with dielectic loding ( ), the dius of the dielectic ediu is, the dius of the tjectoy pojection of cycloton electon be in X-Y plne is, it oves t velocity v ve ve bove the gphene lye with cycloton tjectoy.. Dispesion eqution As shown in Fig.., the schee cn be divided into two egions. Without chge souces, solving Eq. (.) nd (.) togethe with the boundy conditions, the fields in egion nd cn be obtined. The fcto n the egion ( ): e jkjj t is oitted. 4

E A k J ( k ) H A k J ( k ) c c c c E ja k kc J ) A J ) E A k J ) jakc J ) H A J ) ja kkc J ) H A j kcj ) Ak J ) (.) whee: k c k k, J ) [ J ) J )], J ) J ) n the egion ( ): whee: E A k K ( k ), H A k K ( k ) 3 c c 4 c c E jkkc A3 K ) A4 K) E k A3 K) jkca4 K ) H jkkca4 K ) A3 K) k H A4 K) j kca3 K ) kc k k, K ) [ K ) K )], K ) K) (.) Assuing tht onolye gphene is toiclly thin, it cn be egded s conductive sufce with conductivity. The boundy conditions e shown s below: E E, E E ( H H ) E, ( H H ) E g g g (.3) Afte substituting Eq. (.) nd (.) into the bove boundy conditions, the dispesion eqution of cicul cylindicl gphene stuctue with dielectic loding is obtined, 5

k J ( k ) k j K k k J k k J k K k Q c c [ ( c ) c ( c )] g c ( c ) ( c ) kck ) [ k J ( k ) k K ( k ) Q ] g c c c (.4) k [ K ( k c ) Q J ( k c )] [ k K ( k ) Q k J ( k ) j k J ( k )] whee: Q c c c c g c c k ( )( ) c jk J kc kc J ), kc Q kck kc kck kc ( ) ( ). Powe density.. Excittion of fundentl SPPs ode by linely oving electon be Becuse the fundentl SPPs ode (=) is tnsvese gnetic (TM) ode, it cn be excited by the TM evnescent fields poduced by linely oving electon be (the pticul cse of CEB without otting velocity). The electognetic fields poduced by the linely oving electon cn be expessed s below []. E i H Ae i ci J ( k ) H k () jk ci () Ji) H kci j A e k ci jk q v whee A ( / ), c ci J ( k ) H k J k H k () ci () ( ci ) ci Fo the fundentl ode, the TM fields in the stuctue cn be expessed s below. n the egion ( ):, (.5) E A k J ( k ) c c E ja k k J ( k ) c c H A j k J ( k ) c c (.6) whee: k c k k, J ) J ) n the egion ( ): E A k K ( k ) 3 c c E jk k A K ( k ) c 3 c H j k A K( k ) c 3 c (.7) 6

whee: kc k k, K ) K) The boundy conditions cn be witten s E i i ( E E ),( H H H ) E (.8) g By solving the bove boundy conditions, the fields coefficients cn be deteined... Excittion of hybid odes by CEB The electognetic fields in the stuctue nd poduced by CEB e obtined in the pevious pgphs. The boundy conditions cn be witten s, i i E ( E E ), E ( E E ), ( H H H ) E, ( H H H ) E i i g g (.9) Subitting the electognetic fields into the bove boundy conditions, the field coefficients cn be obtined. The excited SPPs e tnsfoed into dition in the dielectic ediu when the Cheenkov dition condition is stisfied. The dition powe density cn be clculted by the following eqution, Re[ * ] P E H dd (.) Appendix : The theoy of cicul cylindicl double-lye gphene stuctue Fig. 3. Schetic of cicul cylindicl double-lye gphene stuctue with dielectic loding, the dius of the dielectic ediu is, the dielectic fil is in the egion b, nd the dius of the tjectoy pojection of electon be is. 7

3. Dispesion eqution As shown in Fig. 3., the schee cn be divided into thee egions. Without chge souces, solving Eq. (.) nd (.) togethe with the boundy conditions, the fields in egion, nd cn be obtined. The fcto n the egion ( ): e jkjj t c c c c is oitted. E A k J ( k ) H A k J ( k ) E ja k kc J ) A J ) E A k J ) jakc J ) H A J ) ja kkc J ) H A j kcj ) Ak J ) (.) whee: k c k k, J ) [ J ) J )], J ) J ) n the egion ( b ): E A k ( k ) A k K ( k ), H A k ( k ) A k K ( k ) 3 c c 4 c c 5 c c 6 c c E jkkc[ A3 ) A4 K )] [ A5 ) A6 K)] E k [ A3 ) A4 K )] jkc [ A5 ) A6 K )] (.) H jkkc[ A5 ) A6 K )] [ A3 ) A4 K)] k H [ A5 ) A6 K)] j kc[ A3 ) A4 K )] whee: kc k k, ) [ ) )], ) ), K ) [ K ) K )], K ) K) n the egion ( ): b 8

E A k K ( k ), H A k K ( k ) 7 c3 c3 8 c3 c3 E jkkc3a7 K 3) A8 K3) E k A7 K3) jkc3a8 K 3) H jkkc3a8 K 3) A7 K3) k H A8 K3) j kc3a7 K 3) whee: kc3 k k, K 3) [ K 3) K 3)], K 3) K3) The boundy conditions cn be witten s, (.3) E E, E E, ( H H ) E, ( H H ) E g g E E, E E, b b b b ( H H ) E, ( H H ) E g g b b b b (.4) (.5) Subitting the electognetic fields into the bove boundy conditions, the dispesion eqution cn be obtined. 3. Powe density 3.. Excittion of fundentl ode by linely oving electon be Fo the fundentl ode (=), the electognetic fields in the stuctue cn be expessed s, n the egion ( ): E A k J ( k ) c c E ja k k J ( k ) c c H A j k J ( k ) c c (.6) whee: k c k k, J ) J ) n the egion ( b ): 9

whee: E A k ( k ) A k K ( k ) 3 c c 4 c c E jk k [ A ( k ) A K ( k )] c 3 c 4 c H j k [ A ( k ) A K( k )] c 3 c 4 c (.7) kc k k, ) ), K ( k ) ( ) c K k c n the egion ( ): b E A k K ( k ) 7 c3 c3 E jk k A K ( k ) c3 7 c3 H j k A K( k ) c3 7 c3 (.8) with k k k c3 The boundy conditions cn be witten s, E E, ( H H) E g E ( E E), ( H H H) E (.9) i i b g b b b By subitting the electognetic fields into the bove boundy conditions, the fields coefficients cn be obtined... Excittion of hybid odes by CEB The hybid odes cn be excited by the CEB, the boundy conditions cn be witten s, E E, E E ( H H ) E, ( H H ) E g g (.) i i E E E, E ( E E ) b b b b ( H H H ) E, ( H H H ) E (.) i i g g b b b b Subitting the electognetic fields into the bove boundy conditions, the fields coefficients cn be obtined. Appendix V: nfluence of elxtion tie of gphene on dition pefonce The pefonce of dition fo electon be excited SPPs is inly dependent on the qulity of gphene, especilly its elxtion tie. n this section, we

discuss the influence of on the dition pefonce. n the nuscipt, the vlue of. ps is used bsed on ecent high-qulity gphene []. While, it is known tht the ost CVD gphene is with low-qulity, nd its elxtion tie would futhe decese when plced on dielectic substte due to the extinsic sctteing [3], which leds to the sll elxtion tie of the ode of. ps [4]. Howeve, ecent dvnces in CVD fbiction of high-qulity gphene ke long elxtion tie vilble [5]. Fig. 4. shows the esults of nolied ttenution constnts of gphene SPPs nd Fouie spect of dition intensity s function of fequency fo elxtion tie. ps nd.ps, espectively. And othe petes e the se s those used in Fig. in the nuscipt. The ttenution constnt of elxtion tie.ps is uch lge thn tht of.ps, especilly in the long-wvelength egie whee the dition intensity is vey wek. Fig. 4. () Nolied ttenution constnts of gphene SPPs nd (b) Fouie spect of dition intensity s function of fequency fo elxtion tie.ps nd.ps. Howeve, the sitution is uch bette when the elxtion tie inceses to. ps, s shown in Fig. 4.. The tio of ttenution constnts fo elxtion tie. ps nd. ps is less thn, so the dition intensity is stonge fo.ps

Rdition intensity thn tht fo thitieth of tht fo.ps. But the dition intensity t the dition pek is only. ps. Fig. 4. 3 shows Fouie spect of dition intensity fo diffeent elxtion ties. Highe elxtion tie leds to stonge nd shpe dition pek, which hs been nlyed by [6]. This is cused by stonge esonnt stength nd lowe esonnt dping. While the dition intensity fo sevel ties highe thn those fo obtined in wide elxtion tie nge.. ps is.4 ps. So high pefonce dition cn be Fig. 4. () Nolied ttenution constnts of gphene SPPs nd (b) Fouie spect of dition intensity s function of fequency fo elxtion tie.ps nd.ps. 6 5 4 3. ps. ps.8 ps.6 ps.4 ps 3 4 5 Fequency (TH) Fig. 4. 3 Fouie spect of dition intensity s function of fequency fo diffeent elxtion ties. Refeence:. Liu, S. et. el. Phys. Rev. Lett. 9, 539 ().. Den, C. et l. Boon nitide substtes fo high-qulity gphene electonics.

Ntue Nnotechnol. 5, 7-76 (). 3. Chen, J., Jng, C., Xio, S., shigi, M. & Fuhe, M. Nt. Nnotechnol. ntinsic nd extinsic pefonce liits of gphene devices on SiO. 3, 6 (8). 4. Tssin, P., Koschny, T. & Soukoulis, C. Gphene fo Tehet Applictions. Science. 34, 6 (3). 5. Ho,Y. et l. The ole of sufce oxygen in the gowth of lge single-cystl gphene on coppe. Science. 34, 7 (3). 6. Zhn, T. et l. Tunble tehet dition fo gphene induced by oving electons. Phys. Rev. B. 89, 45434 (4). 3