arXiv:submit/1855971 [physics.optics] 6 Apr 2017 oi metamaterial. monic epsilon-near- Keywords: refraction. as negative eliminate metamateria such and plasmonic completely properties, normal a propagation are construct wave and losses tunable mon order combine vides metallic periodic we the a Secondly, in which propagation. graphene plasmon for lossless value to gain resulting a a properties is polariton th pre there plasmon the applications surface how two affects examine dielectrics we outline Firstly, (active) we for investigation. nanoplasmonics addition, theoretical surfaces, in our In metallic important plete on very nature. plasmons is subwavelength with which its couple polariton, cer to plasmon under a Furthermore, able surface modes world”. is quantum “classical light purely a stances, in are exist non-loca plasmons orig not as graphene gas would they appear hand, electron two-dimensional effects other a quantum the in On whereas, as laws, well classical fi as plasmo Wefrom and metals of graphene. bulk gas investigation namely in electron semi-metal, analytical plasmons two-dimensional two-dimensional for famous a most used in the is metals, appro limit phase three-dimensional random wavelength th in The through long metals. context the sc in this electrons of in in valence interest emerges with plasmonics light the of of triggered tion area has The matter time. with long light of interaction The Abstract 1 rpeeadatv eaaeil:theoretical metamaterials: active and Graphene colo niern n ple cecs avr University, Harvard Sciences, Applied and Engineering of School 2 eateto hsc,Uiest fCee ealo 71003, Heraklion Crete, of University Physics, of Department aisMattheakis Marios 3 eateto hsc,HradUiest,Cambridge, University, Harvard Physics, of Department admpaeapoiain rpee andeetis plas- dielectrics, gain graphene, approximation, phase Random ehd n hsclproperties physical and methods * mi drs:[email protected] address: Email abig,Msahsts018 USA 02138, Massachusetts Cambridge, ascuet 23,USA 02138, Massachusetts 1,2,* iro .Tsironis P. Giorgos , Kaxiras Greece 1 1,3 2 n Efthimios and , dw n that find we nd eobehavior, zero corrections. l ec fgain of sence ancircum- tain etssfor ientists htpro- that l hwthat show lyr of olayers interac- e d thus, nd, ximation al in nally tcom- at iga ming u to due inate ns d 1 Introduction
The interaction of light with matter has triggered the interest of scientists for long time. The area of plasmonics emerges in this context through the inter- action of light with electrons in metals while a plasmon is the quantum of the induced electronic collective oscillation. In three-dimensional (3D) metals as well as in a two-dimensional electron gas (2DEG), the plasmon arises classi- cally through a depolarized electromagnetic field generated through Coulomb long range interaction of valence electrons and crystal ions [1]. Under certain circumstances, light is able to couple with plasmons on metallic surfaces, form- ing a surface plasmon polariton (SPP) [2–4]. The SPPs are very important in nanoplasmonics and nanodevices, due to their subwavelength nature, i.e. because their spatial scale is smaller than that of corresponding free electro- magnetic modes. In addition to classical plasmons, purely quantum plasmon modes exist in graphene, the famous two-dimensional (2D) semi-metal. Since we need the Dirac equation to describe the electronic structure of graphene, the resulting plasmons are purely quantum objects [5–8]. As a consequence, graphene is quite special from this point of view, possessing exceptional optical properties, such as ultra-sub-wavelength plasmons stemming from the specifics of the light-matter interaction [7–10]. In this Chapter, we present basic properties of plasmons, both from a classi- cal standpoint but also quantum mechanically using the random phase approx- imation approach. Plasmons in 3D metals as well as in 2DEG originate from classical laws, whereas, quantum effects appear as non-local corrections [11–13]. In addition, we point out the fundamental differences between volume (bulk), surface and two-dimensional plasmons. We show that graphene plasmons are a purely quantum phenomenon and that they would not exist in a “classical world”. We then outline two applications that complete our theoretical investi- gation. Firstly, we examine how the presence of gain (active) dielectrics affects SPP properties and we find that there is a gain value for which the metallic losses are completely eliminated resulting to lossless SPPs propagation [3]. Sec- ondly, we combine monolayers of graphene in a periodic order and construct a plasmonic metamaterial that provides tunable wave propagation properties, such as epsilon-near-zero behavior, normal and negative refraction [9].
2 Volume and surface plasmons in three dimen- sional metals 2.1 Free collective oscillations: Plasmons Plasma is a medium with equal concentration of positive and negative charges, of which at least one charge type is mobile [1]. In a classical approach, metals are considered to form a plasma made of ions and electrons. The latter are only the valence electrons that do not interact with each other forming an ideal negatively charged free electron gas [1,14]. The positive ions, i.e. atomic nuclei,
2 are uniformly distributed forming a constant background of positive charge. The background positive charge is considered to be fixed in space, and, as a result, it does not response to any electronic fluctuation or any external field while the electron gas is free to move. In equilibrium, the electron density (plasma sea) is also distributed uniformly at any point preserving the overall electrical neutrality of the system. Metals support free and collective longitudinal charge oscillation with well defined natural frequency, called the plasma frequency ωp. The quanta of these charge oscillations are plasmons, i.e. quasi-particles with energy Ep = ~ωp, where ~ is the reduced Plank constant. We assume a plasma model with electron (and ion) density n. A uniform charge imbalance δn is established in the plasma by displacing uniformly a zone of electrons (e.g. a small slab in Cartesian coordinates) by a small distance x (Fig. 1). The uniform displacement implies that all electrons oscillate in phase [2]; this is compatible with a long wavelength approximation (λ /α , p → ∞ where λp is the plasmon wavelength and α is the crystal lattice constant); in this case the associated wavenumber q (Fig. 1(b)) is very small compared with | | Fermi wavenumber kF , viz. q/kF 0 [7]. Longitudinal oscillations including finite wave vector q will be taken into→ account later in the context of quantum mechanics. The immobilized ions form a constant charge density indicated by en, where e is the elementary charge. Let x(t) denote the position of the dis- placed electronic slab at time t with charge density given by eδn(t). Due to the electron displacement, an excess positive charge density, is creat− ed that is equal to eδn(t), which in equilibrium, δn = 0, reduces to zero. Accordingly, an elec- tric field is generated and interacts with the positive background via Coulomb interaction, forcing the electron cloud to move as a whole with respect to the immobilized ions, forming an electron density oscillation, i.e. the plasma oscil- lation. The polarized electric field is determined by the first Maxwell equation as E =4πeδn, (1) ∇ · in CGS units1. The displacement x(t) in the electronic gas produces an electric current density J = e(n + δn) ˙x en ˙x (since δn/n 0), related to the − ≈ − → electron charge density via the continuity equation J = e∂tδn. After integration in time we obtain ∇ · −
δn = n x (2) ∇ · Combining the equations (1) and (2), we find the electric field that is induced by the electron charge displacement, i.e.
E =4πenx. (3)
Newtonian mechanics states that an electron with mass m in an electric field E obeys the equation mx¨ = eE, yielding finally the equation of motion − mx¨ +4πe2nx =0, (4)
1 For SI units we make the substitution 1/ε0 = 4π.
3 Figure 1: (a) A charge displacement is established by displacing uniformly a slab of electrons a small distance x, creating a polarized electric field in the solid. (b) A plasma longitudinal oscillation electric field in the bulk of a solid. The arrow indicates the direction of displacement of electrons and of the wavevector q, while the double faced arrow shows the plasmon wavelength λp. indicating that electrons form a collective oscillation with plasma frequency
4πe2n ωp(0) = . (5) r m where ω (0) ω (q = 0). The energy E = ~ω is the minimum energy nec- p ≡ p p p essary for exciting a plasmon. Typical values of plasmon energy Ep at metallic densities are in the range 2 20 eV. Having shown that an electron− gas supports free and collective oscillation modes, we proceed to investigate the dynamical dielectric function ε(q,ω) of the free electron gas. The dielectric function is the response of the electronic gas to an external electric field and determines the electronic properties of the solid [1, 11, 15]. We consider an electrically neutral homogeneous electronic gas and introduce a weak space-time varying external charge density ρext(x,t) [14]. Our goal is to investigate the longitudinal response of the system as a result of the external perturbation. In free space, the external charge density produces an electric displacement field D(x,t) determined by the divergence relation ∇ · D = 4πρext. Moreover, the system responds and generates additional charges (induced charges) with density ρind(x,t) creating a polarization field P(x,t) defined by the expression P = ρ [1]. Because of the polarization, the total ∇· − ind charge density inside the electron gas will be ρtot = ρext + ρind, leading to the screened electric field E, determined by E =4πρtot. The fundamental relation D = E +4πP is derived after combining∇· the aforementioned field equations. The dielectric function is introduced as the linear optical response of the system. According to the linear response theory and taking into account the non-locality in time and space [2, 14], the total field depends linearly on the external field, if the latter is weak. In the most general case we have
∞ D(x,t)= dx′ dt′ε(x x′,t t′)E(x′,t′), (6) − − Z Z−∞
4 where we have implicitly assumed that all length scales are significantly larger than the crystal lattice, ensuring homogeneity. Thence, the response function depends only on the differences between spatial and temporal coordinates [2,8]. In Fourier space the convolutions turn into multiplications and the fields are decomposed into individual plane-wave components of the wavevector q and angular frequency ω. Thus, in the Fourier domain the equation (6) reads D(q,ω)= ε(q,ω)E(q,ω). (7) For notational convenience, we designate the Fourier transformed quantities with the same symbol as the original while they differ in the dependent vari- ables. The Fourier transform of an arbitrary field F(r,t) is given by F(r,t) = F(q,ω)ei(q·r−ωt)dqdt where ω, q represent the Fourier transform quantities. Hence The Fourier transform of the divergence equations of D and E yields R iq D(q,ω)=4πρ (q,ω) (8) − · ext iq E(q,ω)=4πρ (q,ω). (9) − · tot In longitudinal oscillations the electron displacement field is in the direction of q (Fig. 1(b)), thus, q D = qD and q E = qE, where D(q,ω) and E(q,ω) refer to longitudinal fields.· Combining the· equations (7), (8), (9) yields ρ (q,ω) ρ (q,ω)= ext . (10) tot ε(q,ω)
Interestingly enough, in the absence of external charges, ρext(q,ω) = 0, the equation (10) states that non-zero amplitudes of charge oscillation exist, i.e. ρ (q,ω) = 0, under the condition tot 6 ε(q,ω)=0. (11) In other words, in the absence of any external perturbation, free collective charge oscillations exist with dispersion relation ω(q) that satisfies the condition (11). These are plasmon modes and consequently the equation (11) is referred as plas- mon condition. Furthermore, the condition (11) leads to E = 4πP, revealing that at plasmon frequencies the electric field is a pure depolarization− field [1,2]. We note that due to their longitudinal nature, plasmon waves cannot cou- ple to any transverse wave such as electromagnetic waves; as a result volume plasmons cannot be excited by light. On the other hand, moving charged parti- cles can be used for exciting plasmons. For instance, an electron beam passing through a thin metal excites plasmons by transferring part of its energy to the plasmon excitation. As a result, plasmons do not decay directly via electro- magnetic radiation but only through energy transfer to electron-hole excitation (Landau damping) [2, 8, 14].
2.2 Dynamical Dielectric Function Based on the plasmon condition (11), the problem has been reduced in the calculation of the dynamical dielectric function ε(q,ω). Further investigation of
5 ε(q,ω) reveals the plasmon dispersion relation as well as the Landau damping regime, i.e. where plasmons decay very fast exciting electron-hole pairs [8]. Classically, in the long wavelength limit, the dielectric response ε(0,ω) can be calculated in the context of the plasma model [1,11]. Let us consider the plasma model of equation (4) subjected to a weak and harmonic time-varying external field D(t)= D(ω)e−iωt; the equation (4) is modified to read
mx¨(t)+4πe2nx(t)= eD(t). (12) − Assuming also a harmonic in time electron displacement, i.e. x(t)= x(ω)e−iωt, the Fourier transform of equation (12) yields
mω2 +4πe2n x(q,ω)= eD(q,ω). (13) − − Introducing the equation (3) in (13) and using the relation (7), we derive the spatially local dielectric response
ω (0)2 ε(0,ω)=1 p , (14) − ω2 where the plasma frequency ωp(0) is defined in equation (5). Equation (14) verifies that the plasmon condition (11) is satisfied at the plasma frequency. The dielectric function (14) coincides with the Drude model permittivity. Further investigation of the dynamical dielectric function can be performed using quantum mechanics. An explicit form of ε(q,ω) including screening effect has been evaluated in the context of the random phase approximation (RPA) [8, 12–14] and given by
ε(q,ω)=1 υ (q)χ (q,ω) (15) − c 0 where vc(q) is the Fourier transform of the Coulomb potential and χ0(q,ω) is the polarizability function, known as Lindhard formula [8,12–14]. The Coulomb potential in two and three dimensions, respectively, reads
2 2πe (2D) |q|εb υc(q)= 2 (16) 4πe 2 (3D) |q| εb where εb represents the background lattice dielectric constant of the system. In RPA approach, the dynamical conductivity σ(q,ω) reads [8]
iωe2 σ = χ (q,ω), (17) q2 0 revealing the fundamental relation between ε(q,ω) and σ(q,ω) that also de- pends on system dimensions; we have finally
q2v ε(q,ω) = 1 + i c σ(q,ω). (18) ωe2
6 In the random phase approximation the most important effect of interactions is that they produce electronic screening, while the electron-electron interaction is neglected. The polarizability of a non-interacting electron gas is represented by Lindhard formula as follows
2 f(ǫk+q) f(ǫk) χ0(q,ω)= − (19) −V ~ω (ǫk q ǫk) + i~η k + X − − where the factor 2 is derived by spin degeneracy (summation over the two possible values of spin s = , ) [8, 13, 14]. The summation is over all the wavevectors k, V is the volume,↑ ↓ i~η represents a small imaginary number to be brought to zero after the summation and ǫk is the kinetic energy for the wave vector k. The carrier distribution f is given by Fermi-Dirac distribu- −1 tion f(ǫk) = (exp[β(ǫk µ)] + 1) , where µ is the chemical potential and − β =1/kBT with the Boltzmann’s constant denoted by kB and T is the absolute temperature. Equation (19) describes processes in which a particle in state k, which is occupied with probability f(ǫk), is scattered into state k + q, which is empty with probability 1 f(ǫk+q). The equations (15)-(19) consist the basic equations for a detailed− investigation of charge density fluctuations and the screening effect, electron-hole pair excitation and plasmons. With respect to condition (11), the roots of equation (15) determine the plasmon modes. Moreover, the poles of χ0 accounts for electron-hole pair excitation defining the plasmon damping regime [12–14]. For an analytical investigation we split the summation of equation (19) in two parts. We make an elementary change of variables k + q k, in the term →− that includes f(ǫk+q), and assume that the kinetic energy is symmetric with respect to the wavevector, i.e. ǫk = ǫ−k. Therefore, the formula (19) yields
2 f(ǫk) f(ǫk) χ0(q,ω)= (20) V ~z (ǫk q ǫk) − ~z + (ǫk q ǫk) k + k + ! X − − X − where z = ω + iη. At zero temperature the chemical potential is equal to Fermi energy, i.e. µ = EF [8, 11, 14], and the Fermi-Dirac distribution is reduced to Heaviside step function, thus, f(ǫk) T =0 = Θ(EF ǫk). The kinetic energy of each electron of mass m in state k is| given by −
~2 k 2 ǫk = | | , (21) 2m hence ~2 2 ǫk q ǫk = q +2k q . (22) + − 2m | | · At zero temperature, because of the Heaviside step function, the only terms that survive in summation (20) are those with k < kF , where kF is the Fermi wave- | | 2 1/2 number and related to Fermi energy by equation (21) as kF = (2mEF /~ ) .
7 Subsequently, we obtain for the Lindhard formula
4 ǫk+q ǫk χ0(q,ω)= 2 − 2 (23) V (~z) (ǫk+q ǫk) |kX|