Open access peer-reviewed chapter

Nanoplasmonics in Metallic Nanostructures and Dirac Systems

By Hari P. Paudel, Alireza Safaei and Michael N. Leuenberger

Submitted: October 17th 2016Reviewed: January 31st 2017Published: June 21st 2017

DOI: 10.5772/67689

Downloaded: 1085

Abstract

In this book chapter, we review some of the progress made in nanoplasmonics and related optoelectronics phenomena in the field of two-dimensional (2D) materials and the recent 3D Weyl semimetals. We give a brief overview of plasmonics for three-dimensional (3DEG) and two-dimensional electron gases and draw comparisons with graphene, 3D topological insulators, 3D Weyl semimetals, and nanoplasmonics in nanogeometries. We discuss the decay of plasmons into electron-hole pairs and the subsequent thermalization and cooling of the hot carriers. We present our recent results in the fields of plasmonics in different nanostructures made of noble metals, such as Silver, and plasmonics in Dirac systems such as graphene and 3D topological insulators. We show a possibility of dynamically shifting the plasmon resonances in hybrid metal-semiconductor nanostructures. Plasmonics in 3D topological insulator and 3D Weyl semimetals have been least explored in nanoplasmonics although it can provide a variety of interesting physical phenomena involving spin plasmonics and chirality. Due to the inherent large spin-orbit coupling, locked spin-momentum oscillations can exist under special conditions and in the presence of an external laser field. We explore symmetric and antisymmetric modes in a slab of 3D TIs and present their dependences on the thickness of the slab.

Keywords

  • surface plasmon polaritons
  • metal
  • Dirac fermions

1. Introduction

Electromagnetic properties of metal-dielectric interfaces have attracted a vast amount of research efforts. Ever since, the work of Mie [1] for small particles and Ritchie [2] for flat interfaces, a wide variety of scientists ranging from physicist, chemists, material scientists to biologists have explored plasmonics-based phenomena and their potential applications in practical life. In nanostructures, under the right circumstances, light waves propagating at metal-dielectric interface excite collective modes of electrons at the metal surface, resulting in the generation of charge density waves called surface plasmons (SPs), which can be divided into localized modes called localized surface plasmons (LSPs) and propagating modes called surface plasmon polaritons (SPPs), propagating along the interface like ripples across the surface of water with an effective wavelength much less than that of incident electromagnetic wave. Free electrons respond collectively by oscillating in resonance with the light waves. In optics, scientists have investigated methods to use plasmonics for concentrating, channeling, and changing the phase of light using subwavelength metallic structures. This would lead to miniaturized plasmonic circuits with length scales much smaller than those in current use [313]. A creatively engineered metal-dielectric interface can generate surface plasmons with the same frequency as the outside electromagnetic waves but with several times shorter wavelength. This interesting phenomenon can be utilized in a way in which surface plasmons carry information in microprocessors faster than current electronic transistors [14]. Plasmonics holds promise for a higher information density than conventional electronics [15]. While this proposed application needs still to be proven feasible, metallic nanostructures much smaller than the wavelength of light have already been successfully built for amplifying signals in surface enhanced Raman spectroscopy (SERS), providing a powerful method to detect a single molecule [1627]. Plasmonically enhanced electric fields are already being used for sensing biomolecules [2836]. The efficient heating property of plasmonic fields can be used for photothermal cancer treatment [37, 38] and also for thermally assisted magnetic recording [3941]. Plasmonic lasers are able to achieve ultrafast dynamics with sub-wavelength mode confinement [4250].

Metallic nanoparticles can be made in different sizes and shapes, and the distance among them can be controlled as well. These parameters can be used to tune the plasmon resonance frequency [5162]. The use of optically excited plasmons as a tunable frequency source that can be mixed with a laser through Raman scattering enables dynamical shifting of the wavelength of light in a controlled manner [63]. Fluegel et al. [64] used a continuous laser beam of a few microwatts power to excite carriers in a carefully engineered narrow GaAs quantum well. These photogenerated carriers are selectively passed through a thick barrier of AlAs into a wide GaAs quantum well in which the two-dimensional electron gas (2DEG) supports high charge density waves with collective electron motion normal to the layer. A signal laser beam operating at a different wavelength undergoes inelastic Raman scattering from the plasmon-phonon (longitudinal optical phonon) modes in the two-dimensional (2D) quantum well. The result is that a signal with 13–15 nm redshifted frequency is generated. Plasmonics can also be exploited in optical tweezers to confine nanoparticles to small dimensions [63]. Grigorenko et al. [65, 66] have made electromagnetically coupled gold pillars. A trapping beam simultaneously excites gap plasmon (GP) modes in the gap between the pillars when they are used as a substrate, resulting in an enhanced field at the trapping site.

In a metallic nanoshell containing a core of dielectric material, due to different dielectric environments in the core and outside of the nanoshell, plasmons are excited in the inner and outer surface of the metallic nanoshell, which can be tuned by varying the ratio of inner to the outer radii of the nanoshell and by altering the dielectric materials [67, 68]. Figure 1 shows that the plasmon resonance energy can be tuned over a wide range by varying the material and/or shaping parameters. Figure 1 compares the plasmon resonance wavelengths of gold, silver, graphene, and topological insulator nanostructures. For a particle diameter much smaller than the wavelength of light, the light-matter interaction leads to an oscillating homogeneous polarization of the particle volume, resulting in an oscillating dipole field. For spherical Au and Ag nanoparticles, the dipole plasmon resonance occurs in the visible part of the spectrum. In particle ensembles, additional shifts are noticed due to electromagnetic coupling between LSP modes. For example, in case of a dimer, plasmons can be viewed as bonding and antibonding combinations, in analogy to molecular orbitals, i.e., hybridization of the individual nanoparticle LSPs occurs, giving rise to resonance shifts following the 1/d3 interaction between two classical dipoles [69]. The observed shifts in LSP resonance and the linewidth narrowing can be understood by means of the Mie scattering theory [1]. Due to the coupling of LSPs, the arrays support one longitudinal and two transverse modes of propagating polarization waves. The transport of energy along such a chain is the main idea for transmitting information using plasmonics [70, 71]. Large heat dissipation has so far slowed down the progress in this field.

Figure 1.

Comparison of the plasmon resonance wavelength spectra of gold, silver, graphene, and topological insulator (TI) nanostructures. The lifetime of plasmon for gold and silver ( τ=1−100 fs) is shorter than that of graphene and TI ( τ=10−100 fs).

Nevertheless, we proposed a method of guiding electromagnetic waves along a chain below the diffraction limit in a controlled manner using a chain of nanoshells [72, 73]. The dynamic control over the plasmon resonances and their coupling gave rise to the idea of an optically controlled plasmonic switch. In this chapter, we present a concept of dynamical control over the plasmon resonances that can be obtained by controlling the dielectric environment of the LSPs using a pump probe technique. In Ref. [74], we show that a shift up to 125 nm can be achieved in an Ag core-TiO2 coated nanostructure.

In addition to that, we also review plasmonics in Dirac systems. Like in the case of a two-dimensional electron gas (2DEG), the surface plasmon mode in 2D materials such as graphene is tunable by changing the gate voltage through shifting the Fermi energy. Exciting surface plasmons give rise to light absorption enhancement in graphene, which can be utilized for photodetectors based on surface plasmon polaritons [75], optical switching of infrared plasmon polariton [76], and THz plasmonic lasing [77]. Plasmonics in Dirac systems show interesting features due to massless electrons around the Dirac nodes. In particular, we focus on surface plasmons in graphene, in Bismuth-based 3D topological insulators (3DTIs), and in 3D Weyl semimetals. Graphene, 3DTIs, and 3D Weyl semimetals are interesting due to their special electronic and optical properties arising from the linear dispersion relation around the Dirac cones in the Brillouin zone [78]. Around these points, energy dispersion of electrons can be described by a low energy Dirac Hamiltonian: HG(k) = h(k) ⋅ σ, where h(k) = [hx(k), hy(k), 0] and σ = [σx, σy, σz]. For a small wave vector, q = k - K, we have h(k) = ħυFq, showing that the electronic states are helical in nature, therefore allowing for “one-way traffic.” This helical nature of electronic states significantly enhances the charge mobility, resulting in reduced plasmonic loss, and therefore high-quality factor in doped graphene.

Due to the excitation of surface plasmons, it is now possible to engineer the behavior of light on nanometer length scales and to increase the light-matter interaction [7982]. This interaction is an outcome of the near-field enhancement close to the metal surface, which also leads to plasmon damping through radiative decay and through nonradiative decay inside the material, due to Landau damping, i.e., creation of electron-hole pairs via interband or intraband transitions, electron-phonon interaction, and boundary effects. The intraband transition happens in the conduction band and the interband transition occurs between other bands (such as the d-band) and the conduction band, as shown schematically in Figure 2a [8385].

Figure 2.

(a) Different types for decay of plasmon. The plasmon is damped by radiative (left) or nonradiative decay through Landau damping (right). (b) Schematic representation of the surface plasmon excitation followed by damping through hot-carrier generation. Schematics for the surface-assisted excitation of hot electron-hole (e-h) pair with energy ħωp. (c) An external electric field incident on a metal nanoparticle with characteristic size of L excites surface plasmon. (d) It leads to geometry-assisted hot e-h excitation (e) followed by momentum relaxation through surface scattering in a region with size of the order of vF/ω.

The radiative decay part of plasmon damping is due to the direct photon emission by coherent electron oscillation. As the size of the nanoparticle increases, the radiative decay of the plasmon is more significant. For larger nanoparticle elements, the radiative decay component is the main reason of plasmon resonance broadening and weakening of the dipole strength. In contrast, decreasing the size of the nanoparticle lets the nonradiative component dominate the plasmon decay. For applications in information technology, a slow dephasing of optical polarization by electron oscillation is essential, which is characterized by the dephasing time (T2). The dephasing time can be found experimentally by measuring the plasmon linewidth Г = 2ħ/T2, where T2 is related to the time constant of the inelastic decay of the plasmons T1 and elastic scattering time constant T* by the authors of Refs. [8284].

1T2=12T1+1T.E1

T1 is a combination of time constants related to the radiative and nonradiative decay processes T1 = T1,r + T1,nr. Nonradiative decay of surface plasmon through Landau damping gives rise to electron-hole pair generation. When the energy of these excited carriers is much larger than the thermally excited electron-hole pairs at ambient temperature, they are called hot carriers. This phenomenon of hot-carrier generation finds lots of applications in energy harvesting, photo detection, and photocatalysis. The hot carriers can be used directly to drive chemical reactions [86] in metals and to drive photocurrents in a semiconductor for the photovoltaic effect or for photoelectrochemical processes [87, 88]. The energy distribution of the generated hot carriers is sensitive to the detailed band structure of the metal. For example, the energy of the excited holes in copper and gold is higher than for electrons; however, in aluminum and silver, the electron and hole energies are almost equal. The advantage of the asymmetry between the produced hot electron and holes by the surface plasmon helps to collect more hot carriers before inelastic scattering occurs. The energy distribution of the carriers generated by decay of plasmons, the subsequent elastic and inelastic scatterings and their transport are important for the design of the hot-carrier-based devices. The inelastic scatterings due to the Coulomb interaction lead to the thermalization of the hot carriers among each other, thereby reaching their carrier equilibrium temperature. The inelastic scatterings due to the electron-phonon interaction result in the cooling of the carriers, thereby reaching the equilibrium temperature of the lattice of the material. For example, in graphene, thermalization takes place on a time scale of 100 fs, while cooling is much slower and happens on a time scale of 10 ps.

This book chapter is organized as follows: in Section 2, we discuss the SP resonances in hybrid metal-semiconductor nanostructures. By altering the dielectric environment of nanostructures dynamically using pump-probe techniques, we show that it is possible to shift the SP resonance wavelength. In Section 3, we discuss in detail the plasmon excitations and their damping pathways in a three-dimensional electron gas (3DEG). We discuss both the radiative and nonradiative damping mechanisms of SPs in 3DEGs. In Section 4, we give a brief overview of plasmons in a two-dimensional electron gas (2DEG). In Section 5, we present the size dependent properties of the SPs in nanostructures. Graphene plasmonics and losses are discussed in Section 6. Section 7 is dedicated to the description of the SPs in Dirac systems. We focus on the SPs in 3DTI materials and 3D Weyl semimetals, and we discuss graphene plasmons as a limiting case of the 3DTI plasmons in the limit where the thickness of the 3DTI slab d → 0. Finally, we wrap up the book chapter with the conclusion in Section 8.

2. Surface plasmon resonances in metal nanostructures

A nanoparticle shows tunable optical properties under controlled variation of its geometry. In a pure Ag spherical nanoparticle in vacuum, for example, the plasmon resonance occurs at 320 nm. These plasmon modes are shifted if the nanoparticle is coated with dielectric materials. It has been shown that with increasing shell thickness, the local electric field enhancement factor peak increases and redshifts for ε2 > ε3, whereas the local field enhancement factor peak decreases and blueshifts nonlinearly for ε2 < ε3, where ε2 and ε3 are dielectric functions of the shell and the surrounding materials, respectively [89]. The electric field enhancement factor is defined by EF = |E/E0|2. It can be easily calculated by means of the finite-difference time domain (FDTD) technique. Figure 3 shows the example of an Au nanoparticle with a radius of 50 nm in air with index of refraction n = 1.

Figure 3.

The resonant excitation of the Au nanoparticle causes large local electric fields close to the particle surface, which is obtained by means of finite-difference time domain (FDTD) calculations. Here, EF(x,y) is shown in the xz plane at 494 nm.

The quasi-static approximation provides a good estimate for a nanoparticle size of around 1/10 or smaller of the incident light wavelength. For larger nanoparticles, due to the finite speed of light, retardation effects lead to a redshift of the plasmon resonance [90]. In Ref. [91], authors have found an analytical expression for a spheroid that takes into account the depolarization factors and that gives a good approximation for nanogeometries of size up to 150 nm. Figure 4 shows our results for the local field enhancement in the presence of an Ag nanocube. As expected from electrostatics, the largest enhancement occurs at the vertices of the Ag nanocube.

Figure 4.

Local field enhancement EF for an Ag nanocube of 20 nm size. The surrounding medium has a refractive index of n = 1.65, which is the case for boron nitride (BN). The green arrow points in the propagation direction of the electromagnetic wave, and the red double arrow shows the polarization of the electromagnetic wave.

The optical resonances of a nanoshell exhibit enhanced sensitivity to its local dielectric environment relative to the solid nanoparticle, as shown in Figure 5. For a particle diameter less than the wavelength of light, the light-matter interaction leads to an oscillating homogeneous polarization of the particle volume, resulting in a dipole field. Figure 5a and b shows the dependence of plasmon resonances on the shell thickness and the size of a hybrid metal-semiconductor nanostructure of 15 nm diameter with an Ag core coated by TiO2 shell obtained in the quasi-static approximation. The expressions are given in Ref. [67].

Figure 5.

Surface plasmons (SPs) field enhancement in a hybrid nanoshell with a Ag core and TiO2 shell of size of 15 nm at different shell thickness using the quasi-static field approximation (a), and different positions inside the shell (b) [see Ref. [74]. The inset in Figure 5a shows a scheme of a core-shell nanoparticle.

Plasmon resonances in a nanoshell can be tuned dynamically by letting a pump laser pulse of energy equal to the band gap or above generate electron-hole pairs in a semiconducting material surrounding the nanoshell. A probe laser pulse at a plasmon resonance frequency is used to excite plasmons on the metal surface. The generation of free electron-hole pairs alters the dielectric function of the surrounding semiconducting material. Due to the reduced dielectric function caused by the excitation of the electron-hole pairs, the excitation of surface plasmons by a probe pulse requires a higher energy. The frequency of the probe pulse is smaller than that of the pump pulse ensuring that no excitons are excited in the semiconductor during the probing. The change in the dielectric function of the surrounding medium due to the pump pulse can be calculated using Fermi’s golden rule:

Im(ε(ω))=4π2e2m2ω2k,c,v|e·Pc,v|2f(Ev(k))[1f(Ec(k))]δ(Ec(k)Ev(k)ħω),E2

where Pc,v is a transition dipole matrix element between conduction and the valence band, e is the direction of the polarization and f(E) is the Fermi-Dirac distribution function. At zero temperature, the factor f(Ev(k))[1f(Ec(k))]=1before the optical excitation and becomes less than 1 after optical excitation.

Figure 6 shows shifts in the resonance peak of the surface plasmons occurring at around 620 nm before the generation of excitons in a nanoshell structure (as shown in inset of Figure 5a) with diameter of 15 nm. After the pump pulse, depending on the density of the excitations, the plasmon resonance peaks are excited by the probe pulse shift. The larger the density of excited free electron-hole pairs in the semiconductors, the larger is the blueshift of the plasmon resonance peak. For a density of excitation of 5 × 1021 cm−3, a resonance shift of up to 125 nm can be achieved.

Figure 6.

Local field enhancement factor EF on the surface of a nanoshell of 15 nm diameter. The core is made of Ag with a diameter of 7 nm. At higher density of excitation, the resonance peak is gradually blueshifted.

3. Plasmon excitation and damping for a three-dimensional electron gas (3DEG)

In this section, we discuss the plasmon excitation and their damping pathways in 3D materials made of metal. The plasmonic damping pathways in 3D materials include radiative decay, Landau damping, and resistive loss, as depicted in Figure 7. During Landau damping, plasmon quasi-particles lose their energy by exciting hot electron-hole pairs via direct interband or phonon/geometry-assisted intraband transitions. In the case of geometry-assisted intraband transition, the translational symmetry is broken due to electric field confinement or boundaries of the material [92]. In the case of resistive loss, single carriers, electrons or holes that are the building blocks of the plasmon quasi-particle, are kicked out of the phase-coherent collective plasma oscillation through electron-electron or electron-phonon scattering, giving rise to plasmon damping. In 2D materials, the plasmons follow similar damping pathways. Figure 7 shows the stages of the plasmon decay, the initial nonequilibrium configuration after the excitation, the thermalization, and the cooling of the hot carriers. Some of these damping pathways can be used to inject hot carriers into other materials. For example, at a metal-semiconductor interface, hot electron-hole pairs can be separated by means of the Schottky barrier for the purpose of energy harvesting. In metal-graphene or metal-MoS2 junctions, the surface plasmons can generate hot electron-hole pairs, thereby injecting electrons/holes into n/p-doped 2D materials, giving rise to hot carrier-induced doping [93] or even insulator-to-metal phase transitions [94].

Figure 7.

The stages of energy relaxation of excited plasmon. (a) The incident electric field on a graphene nanostructure excites the surface plasmon. (b) The plasmon quasi-particles lose their energies by different scattering processes to excite hot carriers, which have nonthermal distribution on the order of tens to hundreds of femtoseconds. (c) Electron-electron and other scattering processes redistribute the hot carriers and thermalized them in the order of one to hundreds of picoseconds to generate ling-lived hot carrier distribution. (d) Cooling stages and energy transferring to the substrate happen by phonon and geometry-assisted scattering processes over a longer time scale (hundreds of picoseconds to tens of nanoseconds) for relaxation of the hot carrier distribution to the equilibrium.

The initial distribution of the hot carriers can be estimated using the jellium model for metal nanoparticles and nanoshells [95, 96], but this approach cannot explain the material dependence of this process because the specific band structure of the metal is completely neglected. In order to capture the material properties, it is necessary to combine FDTD calculations for obtaining the plasmon modes with ab initio density functional theory (DFT) calculations to determine the initial energy distribution of generated hot carriers by decay of surface plasmon. Figure 7b shows the schematic of plasmonic hot-carrier generation [97].

The plasmon resonance frequency and dispersion can be obtained by evaluating the dynamic polarizability in the presence of the carrier-carrier Coulomb interaction. The dynamic polarizability in the random phase approximation (RPA) is given by

χRPA(q,ω)=χ0(q,ω)εRPA(q,ω),E3

where χ0(q,ω) is the noninteracting (zeroth-order) polarizability (single pair bubble) and εRPA(q,ω)=εmvc(q)χ0(q,ω), with εm being the permittivity of the environment, and vc(q)=e2/2ε0qthe Coulomb interaction between the carriers. The RPA method corresponds to the expansion of 1/εRPA(q,ω), leading to an infinite power series over the bubble diagrams.

If optical phonons are also considered, the effective dielectric function in the RPA expansion takes the form [98, 99]

εRPA(q,ω)=εmvc(q)χ0(q,ω)εmlvsph,l(q,ω)χ0(q,ω)εmvoph(q,ω)χj,j0(q,ω),E4

where εm is the average of dielectric constants of the environment. The second term represents the effective Coulomb interaction of electrons in graphene, and vc(q)=e2/2qε0is the direct Coulomb interaction. The third term is the effective dielectric function for different phonon modes (l) coming from electron-electron interaction mediated by substrate optical phonons, which couple to the electrons by means of the Fröhlich interaction, i.e.,

vsph,l(q,ω)=|Msph|2Gl0(ω),E5

where |Msph|2 is the scattering and Gl0is the free phonon Green’s function. The last term of Eq. (4) corresponds to the optical phonon mediated electron-electron interaction.

voph(q,ω)=|Moph|2G0(ω).E6

Here |Moph|2 defines the scattering matrix element and Go(ω) is the free phonon Green’s function. In Eq. (4), χj,j0(q,ω)is the current-current correlation function. This description is very general and can be applied to any metallic system.

For the simplest case, when only the carrier-carrier Coulomb interaction is present, we can derive the dynamical plasmon dispersion relation following standard textbooks [100]. The first step is to calculate noninteracting dynamical polarizability.

χ0(q,ω)=1Ωk,sn0k+q,sn0k,sεk+qεkħωiħη,E7

where Ω is the volume of the sample, n0k,sis the equilibrium electron density, and εk is the single-particle energy dispersion of the electrons. Assuming |q| ≪ kF, we can expand χ0 (q, ω) in q by taking advantage of the expansions εk+q=εk+qkεk+and n0k+q=n0k+n0εqkεk+. At zero temperature, we have n0εk=δ(εkεF)and the velocity can be approximated by the Fermi velocity, i.e., kεk=ħvkħvF, leading to the approximation

χ0(q,ω)2(2π)3d3kqvFδ(εkεF)qvFωiη=2(2π)211d(cosθ)0k2dkqvFcosθqvFcosθωiηδ(εkεF)=kF22π2ħvF[23(qvFω+iη)2+25(qvFω+iη)4+]kF3q23π2m(ω+iη)2[1+35q2vF2(ω+iη)2]=n0q2m(ω+iη)2[1+35q2vF2(ω+iη)2].E8

Using the equation kF=(3π2n0)1/3n0=kF33π2for 3D metals, one obtains

χ0(q,ω)ρ0q2m(ω+iη)2[1+35q2vF2(ω+iη)2]=ρ0q2m(ω+iη)2R(q,ω)2E9

with

R(q,ω)2=[1+35q2vF2(ω+iη)2].E10

Consequently, the dynamical polarizability in the RPA and long wavelength regime is

χRPA(q,ω)=χ0(q,ω)ε(q,ω)=ρ0q2m(ω+iη)2R(q,ω)214πe2ρ0m(ω+iη)2R(q,ω)2=ρ0q2R(q,ω)2m(ω+iη)24πe2ρ0R(q,ω)2=ρ0q2R(q,ω)2mωp{1ω+iηωpR(q,ω)1ω+iη+ωpR(q,ω)},E11

where the quantum plasma frequency of a 3D metal in the free electron approximation is given by the Pines and Bohm equation [101, 102]

ωp=4πn0e2m,E12

with n0 being the equilibrium electron density and m the electron mass. Using limη0+1ziη=P(1z)+iπδ(z), it follows that

Im[χRPA(q,ω)]πρ0q2R(q,ω)2mωp{δ(ωωpR(q,ωp))δ(ω+ωpR(q,ωp))},E13

which yields the bulk plasmon dispersion relation for 3D metals, i.e.,

ωq=ωpR(q,ωp)=ωp[1+310q2vF2(ωp+iη)2+].E14

The same result can be obtained by solving εRPA(q,ω)=0. Note that the slope of the parabolic dependence at q=kFis

3vF2kF5ωp=3vF2kF54πe2n0m=9πmvF2kF20e2kF3=9πmvF220e2(mvFħ)2=9πħ220e2m=0.67mscεd.E15

This difference in slopes between the 3D bulk plasmon and the photon-like surface plasmon polariton is clearly visible in Figure 8 (see below).

Figure 8.

Dispersion curves of the RPA bulk plasmon resonance (red), the semi-classical surface plasmon polaritons (blue), and the RPA 2DEG plasmons. The shaded area marks the electron-hole continuum. When the bulk plasmon resonance enters the shaded area, the plasmon decays into electron-hole pairs which is called Landau damping.

It is well known that in 3D metals Landau damping occurs when the plasmon resonance energy enters the electron-hole continuum, which is determined by the condition

ħω=εk+qεk=ħ2(k+q)2ħ2k22m=ħ2q2+2ħ2kq2m=ħ2q2+2ħ2kqcosθ2m=ħ2q(q+2kcosθ)2m0.E16

Since 1cosθ1and 0kkF, the electron-hole continuum is given by the gray shaded area in Figure 8. The plasmon can decay into an electron-hole pair when the plasmon dispersion curve enters the electron-hole continuum limit. This decay corresponds to intraband Landau damping.

The presence of a planar boundary for a 3D metal adds a new mode known as surface plasmon, which propagates at the metal-dielectric interface. Since the electron charge density of a metal leaks outside the interface into the dielectric in the order of 1/kF1Å, a macroscopic description based on Maxwell’s equations is sufficient to understand the surface plasmon qualitatively. Taking the boundary conditions into account, the dispersion relation of the surface plasmon is determined by Ritchie and Eldridge [103]

qx=ωcεmεdεm+εd=ωcω2ωp2(1+1εd)ω2ωp2εdE17

in the Drude model, where εm(ω)=1ωp2/ω2. This biquadratic equation can be solved analytically, yielding

ωsp(qx)=±12εdc2qx2(1+εd)+εdωp2±c4qx4(1+εd)2+2c2qx2(1+εd)εdωp2+εd2ωp4.E18

Since ωsp(qx)0, there are only two physical solutions, which are drawn in Figure 8. While the upper branch corresponds to the photon-like plasmon polariton, the lower branch represents the plasmon-like surface plasmon polariton (SPP) with the following asymptotics:

ωsp(qx)={cεdqx[112(cωp)2qx2+],qxωpcωp1+εd[112(ωpc)2(εd1+εd)2qx2+],qxωpcE19

which in the nonretarded regime reduces to Ritchie’s equation [2],

ωsp=ωp1+εd=ωp2=2πn0e2m=(4e22m3πħ3)1/2εF3/4E20

for wave vectors q in the range of ωs/cqqF, where qF is the Fermi wave vector [104]. Since kF=(3π2n0)1/3in 3D metals, the bulk plasmon frequency exhibits a εF3/4n01/2dependence.

To estimate the plasmon decay rate, the band structure of the materials should be calculated by means of DFT to find out the exact quasi-particle orbitals and energies. Different electronic structures can be used to calculate the band structure of noble metals. To estimate the decay rate, the electronic states and energies of the metals resulting from PBEsol + U DFT method are used in Eq. (2). To find the plasmon decay distribution, energy conservation is used to determine the regions of the first Brillouin zone, which are responsible for the electron-hole generation [97]. Since aluminum, which is a nearly free electron metal, has a band crossing close to the Fermi level at the W point, interband transitions from the valence band with energies ħω from zero to the plasmon energy are allowed for creating a symmetric amount of both hot electrons and holes. There is an extra possibility for interband transitions at the K point, which mostly contribute to the generation of hot holes with energies >2 eV, giving rise to a modest asymmetry between hot electron and hole distributions. For copper, the interband transitions happen at the X, L and K, points, and for gold they occur at the X and L points [97]. In these two noble metals, since all the transitions originate from the “d” band, the hot hole energies are larger than 2 eV and consequently they have larger kinetic energies than the electrons. In metal-semiconductor interface, this asymmetry between hot hole and electron energy distribution plays a dominant role in the hot-carrier collection efficiency. For gold and copper, the holes have more kinetic energy than the electrons. This means that in a metal-semiconductor Schottky junction such as gold-n type gallium arsenide junctions, the hot electrons do not have enough energy to pass the Schottky barrier and therefore need to tunnel or get more thermal energy to overcome Schottky barrier. However, in gold-p type gallium arsenide junctions, the hot holes have sufficient kinetic energy to overcome the Schottky barrier. This means for copper and gold, the metal-p type semiconductor junction is more efficient for the collection of hot holes. But for aluminum and silver, the efficiency is almost the same for hot holes and electrons. According to the Fowler theory [96], the angular distribution of the hot-carriers should be isotropic. However, experiments show that none of the noble metals exhibit an isotropic distribution [97].

In nanoconfined structures, because of lack of translational symmetry, the crystal momentum q is not a good quantum number and there is no momentum conservation anymore. Therefore the transition creating an electron-hole pair is allowed to be nonvertical. This is the case for geometry-assisted intraband transition, where no phonon is required. A vector potential for one quantum of each symmetric and antisymmetric modes can be used in Fermi’s golden rule to calculate the decay rate of an LSP to an electron-hole pair. By decreasing the thickness, the effect of confinement is to be more pronounced and more hot electrons can be generated by intraband transitions. This effect is more noticeable for the antisymmetric mode, because a smaller part of the light is located inside the gold film, which leads to reduction of interband transitions [97].

In metals, plasmons decay is not only in the ultraviolet and visible spectrum but also in the infrared and microwave regimes [105, 106]. Due to the conservation of momentum in infinite crystal lattices, the direct interband transition induced by plasmon decay is only possible for energies larger than the band gap energy. However, for energies below the visible spectrum, typically phonon-assisted and surface-assisted intraband electron-hole pair generations are able to bypass this selection rule [92].

The plasmon decay rate is related to the imaginary part of the dielectric tensor Imε¯(ω)[92], i.e.,

Γ(r)=12πħE(r)Imε¯(ω)E(r).E21

Let us consider a 3D semi-infinite metal slab extending in the negative z-direction with a single surface in the xy-plane at position z = 0. Using the electric field profile of a single quantum of the surface plasmon polariton with wave vector k and frequency ω and integration over the space gives the total nonradiative (nr) decay rate.

Γnr(r)=12L(ω)|γ(z<0)|λImε¯(ω)λE22

The decay rate of plasmon as a function of frequency can be calculated by substituting the experimental data for the complex dielectric function measured by ellipsometry. Within the random phase approximation (RPA), the nonradiative decay rate induced by direct interband transition is [97]

Γnr(direct)=2π2e2L(ω)|γ(z<0)|me2ω2BZdq(2π)3nn(fqnfqn)δ(εqnεqnħω)×|λpnnq|2E23

where pnnq=Ψqnσ(r)iħΨqnσ(r)are momentum matrix elements between the quasi-particle orbitals Ψqnσ(r)and Ψqnσ(r). fqn are the Fermi occupation functions of the quasi-particles, λk^z^k/γ(z<0)is the polarization vector, and γ(z) is the imaginary part of the wave vector in z-direction. Eq. (23) can also be obtained in first quantization and is based on the electric dipole approximation, which means the plasmon mode function slowly changes on the atomic scale. To include the effect of electron-electron and electron-phonon scattering on the decrease in the lifetime of carriers, the δ function is replaced by a Lorentzian with half width ImΣqn+ImΣqn, where ImΣqnis the total carrier lifetime due to the scatterings.

The contribution of the electron-phonon interaction to the plasmon decay through intraband transitions can be calculated in second-order perturbation theory [107, 108], giving

Γnr(phonon)=2π2e2L(ω)|γ(z<0)|me2ω2BZdqdq(2π)6nnα±(fqnfqn)(Nqq,α+1212)×δ(εqnεqnħωħωqq,α)|λn1(gqn,qn1qq,αpn1nqεqn1εqnħω+iη+pnn1qgqn1,qnqq,αεqn1εqnħωqq,α+iη)|2,E24

where ħωk,αis the energy of a phonon with wave vector k and polarization index α, Nk,α is the Bose occupation factor for the phonons, gqn,qnk,αis the electron-phonon matrix element with q and n being the wave vector and band index of the electronic states, and the ± signs distinguish between absorption and emission of phonons.

Excitation of surface plasmon generates a strong field confinement on the surface with the exponential decay length (|γ(z<0)|1)inside the metal, which creates a Lorentzian distribution for the momentum of the plasmon in the z-direction. This momentum distribution leads to diagonal intraband transition and introduces a “surface-assisted” plasmon decay [92]

Imεsurface(ω)=ωp2ω3×(34|γ(z<0)|vF)2k2k2+|γ(z<0)|2,E25

where ωp=4πne2/meis the bulk plasma frequency of the metal, vF=(ħ/me)(3π2n)3/2is the Fermi velocity, and n is the bulk carrier density of the metal.

Numerical studies based on the free electron jellium model show that in nanostructures, due to the localization of electronic states and the nonconservation of the crystal momentum, intraband transitions are enhanced [95, 96, 109]. Using Fermi’s golden rule together with the free electron eigen states and the dipole field profile, the nanoconfinement contribution is

Imεsphere(ω)=ωp2ω3×6vFπ2R,E26

where R is the radius of the spherical nanoparticle.

In the Landau damping theory, the lowest-order processes consisting of direct, phonon-assisted, and surface/geometry-assisted electron-hole pair excitation contribute to the decay of plasmons [102, 110]. Higher-order processes leading to the excitation of many electron-hole pairs or many phonons are suppressed due to the phase-space factors at small energies [111]. Only at large energies, the higher-order processes become significant. Being completely different from Landau damping, another source of plasmon damping is the resistive loss in the metal, which can be calculated by means of the linearized Boltzmann equation in the relaxation time approximation [92], giving

Imεresist(ω)=4πσ0ω(1+ω2τ2),E27

where σ0 is the zero frequency conductivity and τ is the momentum relaxation time, which can be derived from the carrier-carrier Coulomb and electron-phonon matrix elements.

4. Plasmons in a two-dimensional electron gas (2DEG)

According to Ritchie [2] and Stern [112], the dispersion relation for 2DEG plasmons is

ω2DEG=nse22m*ε0ε¯(q,ω)qx,E28

where ns is the surface carrier density. A typical dispersion curve is shown in Figure 8. In the case of a typical metal-oxide-semiconductor (MOS) system, the effective dielectric function is given by

ε¯(q,ω)=12[εSC(ω)+εOX(ω)coth(dqx)],E29

where εSC(ω)and εOX(ω)are the dielectric constants of the semiconductor and oxide layers, respectively. The oxide layer has a thickness of d. For example, Heitmann [113] and Wilkinson [114] were able to observe this type of 2DEG plasmons in AlGaAs-GaAs heterostructures. This result shows clearly that the dispersion curve can be tuned by changing the thickness of the oxide layer and by altering the gate voltage. Similarly, the plasmon dispersion can be tuned in atomically thin 2D materials. Since the dispersion relation in Eq. (28) exhibits a qxdependence for a 2DEG, it is plausible to assume that atomically thin metallic 2D materials follow a similar dependence, which is indeed the case, as we discuss below in Sections 6 and 7. Since in a 2DEG

ns=1Ω|k|kF,s1=2d2k(2π)21=2πkF2(2π)2=kF22π=m*εFπħ2,E30

the 2DEG plasmon frequency exhibits a nsεFdependence. Note that the electron-hole continuum in a 2DEG is determined by a similar formula to Eq. (16). Therefore the gray shaded area in Figure 8 has similar forms for 3DEGs and 2DEGs.

5. Static geometry of metallic objects and environment

When the size of the nanoparticle is much smaller than the wavelength of the incident light, the particle exhibits a dipolar oscillation mode (Fröhlich mode). As the diameter of the nanoparticle is increased, the electrostatic limit is not a good approximation anymore and the multipolar oscillation modes start to appear. Excitation of these modes gives rise to the broadening of the resonance [115], as seen from Figure 9b. The LSP resonance redshifts with increasing diameter of the sphere, which is due to retardation effects [90]. Decreasing the size of nanoparticle less than mean free path of electrons moves the material band structure and dielectric function away from the bulk properties and increases the surface scattering that gives rise to broadening of the absorption spectrum [115], as shown in Figure 9a. The internal field enhancement of an illuminated spherical nanoparticle is

Figure 9.

Calculated extinction coefficient per unit volume for a spherical Aluminum nanoparticle of various diameters. (a) For the case when the size of the nanoparticle becomes smaller than the mean free path of the electron, the linewidth becomes wider due to increased surface scattering. (b) For the case when the size of the nanoparticle becomes much larger than the mean free path of the electron, excitation of multipole modes leads to broadening of the linewidth and the retardation effect leads to a redshift of the LSP resonance peak.

Ei(ω)E0=3εmε(ω)+εm,E31

where ε and εm are the dielectric functions of the nanoparticle and the medium, respectively. An ellipsoidal nanoparticle with a shape defined by (x/a)2+(y/b)2+(z/c)2=1supports an internal field enhancement given by

Ei(ω)E0=εmεm+L(ε(ω)εm),E32

where L is the shape factor, whose value lies between 0 and 1. For the resonance frequency, the real part of the denominator is zero, which means the resonance occurs at frequencies when

ε(ω)=εm(11L)E33

is satisfied. The shape factor for the axis parallel to the polarization of the incident light (e.g., the axis a) is given by

L=abc20dq(a2+q)f(q),E34

where

f(q)=[(a2+q)(b2+q)(c2+q)]12.E35

For the core-shell spherical nanoparticles the resonance condition is given by

εc=2εs[εs(1f)+εh(2+f)εs(1+2f)+2εh(1f)],E36

where f=(rc/rs)3is the core fill fraction and εc, εs, and εh are the core, shell, and host dielectric functions, respectively. The resonance frequency is tunable by modifying the shell thickness and the radius of the core. As shown in Figure 10, due to the hybridization of the plasmon modes, even in the electrostatic limit there are multiple resonance frequencies [116].

Figure 10.

Energy-level diagram of plasmon hybridization in metallic core-shell nanoparticles. The hybridization results from the coupling between the cavity and sphere surface plasmons.

6. Plasmon theory for graphene

Graphene is a two-dimensional (2D) material comprised of a single layer of carbon atoms in a honeycomb lattice. It has unique electrical, optical, and mechanical properties due to its tunable band dispersion relation and atomic thickness. Because of its unique band structure graphene possesses a very high mobility and a fast carrier relaxation time [117121], making it an attractive candidate for ultrafast electronics and optoelectronics. Exciting surface plasmons on graphene is a distinct technique to increase absorption with low damping rate. The surface plasmon couples the electromagnetic (EM) wave to the conductive medium, giving rise to direct absorption of light by monolayer graphene and providing the opportunity of electrical tunability of the plasmon resonance frequency, high degree of electric field confinement, and low plasmon damping rate [122125]. The increased light-matter interaction results in an enhanced spontaneous emission rate close to the nanostructure edges [126, 127]. Recent experiments have achieved an absorption of 90% in the mid-IR range by connecting graphene with high carrier mobility to a silicon diffractive grating [128] and designing graphene nanoribbons [128, 129], nanodisks [130], and antidot array [125] theoretically. These high carrier mobilities can be achieved only for mechanically exfoliated graphene. Exciting plasmon on CVD-grown monolayer graphene with lower mobility than the mechanically exfoliated one reduces the absorbance to 19% and 28% for graphene nanoribbons [99, 127, 131] and nanodisks [132, 133], respectively. We show in experiments that the coupling of a patterned CVD-grown graphene sheet to an optical cavity amplifies the excited LSPs and enhances the light absorption to a current world record of 45% [134]. We also show that the theoretically achievable enhancement is 60% for a square lattice of holes [134].

The electric current of graphene in the interaction picture is given by

jxx(ω)=σ(ω)E(ω),E37

where E(ω) is the in-plane electric field and σ(ω) is the interband optical conductivity

σinter(ω)=ie2ωAkvx2ρvvρccεc,kεv,kħ(ω+iη)E38

with A being the cross-section area, vx the velocity operator along x-direction, ρvv and ρcc the digonal elements of the density matrix, εc,k /εv,k the conduction/valence band energy, and η an infinitesimal value. Using the energy dispersion, εn,k=nvFħ|k|with n=+1/1for conduction/valence band and vF=106m/s, and replacing the summation by the integration in Eq. (38) leads to

σinter(ω)=e2vf22ωg4π02πkdk(ρvvρcc)δ(εc,kεv,kħω)=e2vf2ω0kdk(f(εv)f(εc))δ(εc,kεv,kħω)=e24ħ[f(ħω2)f(ħω2)],E39

Where g = 4 is the graphene degeneracy and ρ is replaced by the Fermi distribution function f. By inserting the explicit expression of the Fermi distribution, Eq. (39) is reduced to

σinter(ω)=e24ħsinh(ħω2kBT)cosh(εF2kBT)+cosh(ħω2kBT),E40

where εF is the Fermi energy of graphene, kB is the Boltzmann constant and T is the temperature. Eq. (40) is the graphene optical conductivity due to interband transition, which is valid in the visible light spectrum. This result can be derived from the Kubo formula too. When the energy of the incident light (ħω) and kBT are of the same order and εF is much larger than kBT, the interband optical conductivity is constant, i.e.,

σinter(ω)=e24ħ.E41

For infrared and THz radiation, the Fermi energy can be tuned to become much larger than the incident photon energy, and therefore due to Pauli blocking there are only intraband transitions. According ot the Boltzmann equation and under the relaxation time approximation, the carrier distribution in the presence of a constant electric field with x-polarization is given by

f(k)=f0(k)qExvx(k)τ(k)f0(k)ε,E42

where f 0(k) is the carrier distribution in the absence of the electric field, q is the carrier charge, and τ(k) is the relaxation time in the presence of electron-phonon interaction, electron-impurity interaction, and other scattering processes. By using Eq. (42), the electric current is given by

jx=qAkf(k)νx(k)=q2ExAkνx2(k)τ(k)f0(k)ε.E43

In the absence of the external electric field, the net electric current is zero, so the summation of the first part of Eq. (43) in the electric current is zero. Since f0(k)εis so narrow around the Fermi surface, only the wave vectors near the Fermi energy contribute to the integration, and τ(k) is approximately constant. By replacing νx2(k)with νF2/2, the intraband optical conductivity becomes

σintra=q2τAνF22kf0(k)ε=q2τvF22π20dk2πkf0(k)ε.E44

For each k there are two energy values in the conduction and the valence bands. Using the linear dispersion relation, the intraband optical conductivity is reduced to

σintra=q2τπħ2(0εdεf0(ε)ε+0εdεf0(ε)ε)=q2τπħ2|ε|dεf0(ε)ε=e2τπħ22kBTln(2cosh(εF2kBT)),E45

where q is replaced by the electron charge (−e). Assuming εFkBT, which is usually the case, the intraband optical conductivty simplifies to

σintra=e2πħ2εF.E46

In the presense of an oscillating electric field, the relaxation time is a complex function (τ1τ1iω)and the intraband conductivity is given by

σintra(ω)=e2πħ22kBTτ1iωln(2cosh(εF2kBT)),E47

which in the case of εFkBTis reduced to [118, 120]

σintra(ω)=e2πħ2εFτ1iω,E48

where τ is determined by impurity scattering and electron-phonon interaction τ1=τimp1+τe-ph1.

The dielectric function of graphene can be obtained via its AC conductivity by means of [135]

ε(ω)=2.5+iσ(ω)ε0ωdE49

where εg = 2.5 is the dielectric constant of graphite. Here d is the thickness of graphene.

The bandstructure of graphene is linear in the tight-binding approximation, as shown in Figure 11b. The dispersion relation for the TM mode in the geometry depicted in Figure 11a, which consists of graphene surrounded by dielectrics with constants εr1 and εr2, is given by [122]

εr1q2εr1ω2c2+εr2q2εr2ω2c2=σ(q,ω)iωε0,E50

Figure 11.

(a) Schematic of the patterned graphene layer surrounded by two media of dielectrics εr1 and εr1. (b) The sketch of intraband (green arrow) and interband (red arrow) transitions after absorbing a photon on the Dirac cone. (c) The graphene plasmon loss function for the structure shown in part (a) without plasmon-optical phonon coupling. This function shows the graphene plasmon dispersion relation. (d) The graphene dispersion relation in presence of substrate/graphene optical phonons.

where σ(ω,q) is the conductance of graphene. The semiclassical conductivity σ(ω) does not depend on q. In the nonretarded regime (q = ω/c), the dispersion relation of the TM mode is reduced to [122]

qε0εr1+εr222iωσ(q,ω).E51

The larger the q (smaller plasmon wavelength), the higher is the transversal confinement of the TM mode due to excitation of the surface plasmon (SP) on graphene, which is similar to the case of regular metals. The main difference of SPs on metal and 2D materials is due to their band structures. While normal metals have parabolic bands, graphene has Dirac cones, which leads to a novel dispersion relation for the SPs on graphene, as shown in Figure 11c. The interactions between the substrate and graphene optical phonons and the graphene plasmons open a band gap in the graphene plasmon dispersion relation, as shown in Figure 11d, showing the plasmon loss function in the presence of plasmon-optical phonons interaction given by Eqs. (4)–(6) [134]. The electrons in 2D materials are bound in the normal direction, which leads to more confinement of the coupled electromagnetic wave due to plasmon excitation. In graphene, thermalization takes place on a time scale of 100 fs, while cooling is much slower and happens on a time scale of 10 ps, as seen from Figure 11 [136, 137].

There are two different approaches to obtain the dispersion relation of the surface plasmon. In the semiclassical approximation, the Drude-like conductivity is used to obtain the plasmon dispersion relation [122, 138]

q(ω)=πħ2ε0(εr1+εr2)e2εF(1+iτω)ω2E52

and the plasmon loss

ReqImq=ωτ=2πcτλair.E53

Another approach that works for both intraband and interband regimes is based on the selfconsistent linear response theory, known as random phase approximation (RPA) along with the relaxation time approximation defined by Mermin [139]. The dispersion relation of the plasmon can be obtained by solving

εRPA(q,ω)=0,E54

with the complex wave vector q=q1+iq2. Considering the Coulomb interaction of electrons in graphene and the medium dielectric, first and second terms of Eq. (4), one obtains the plasmon dispersion relation

εr1+εr12+e22ε0q1Reχ(q1,ω)=0E55

and the plasmon loss relation [122]

q2=Imχ(q1,ω)+1τωReχ(q1,ω)+1ωτReχ(q1,ω)[1χ(q1,ω)χ(q1,0)]1q1Reχ(q1,ω)q1Reχ(q1,ω).E56

The plasmon losses (Req/Imq), wave localization (λair/λp), and the group velocity of the graphene plasmon (νg=ω/q)are calculated by means of the semiclassical approximation and RPA, as shown in Figure 12 for the Fermi level of εF=0.135eVand relaxation time of τ=1.35×1013s[122]. Below the interband regime, the plasmon loss, the wave localization, and the group velocity calculated by RPA are in very good agreement with the semiclassical approach.

Figure 12.

(a) Plasmon loss, (b) filed confinement, and (c) the group velocity of the graphene surface plasmon for εF=0.135eV and the carrier mobility of μ=10000cm2/Vs. The solid and dashed line represent the RPA and semiclassical results, respectively, and the rose shaded area is the high-interband loss region.

For 2D materials such as graphene, the amount of plasmon loss can be calculated by using the effective dielectric function of the material. The dynamical polarizability

χ(q,iωn)=1A0βdτeiωnτTρq(τ)ρq(0)E57

determines several important quantities, such as the effective electron-electron interaction, plasmon spectra, phonon spectra, and Friedel oscillations. ωn=2πnβare the Matsubara frequencies, ρq is the density operator in q-space and A denotes the area. This quantity is calculated in the canonical ensemble for both of the sub-lattice density operators (ρ = ρa + ρb). The dynamical polarizability in zeroth-order in the electron-electron interaction in the long wavelength limit is [140]

χ0(q,iωn)=gsgv4π2d2ks,s=±fss(k,q)f(εs(k))f(εs(|k+q|))εs(k)εs(|k+q|)+iħωn,E58

where gs=gv=2are the spin and valley degeneracy, f(εs(k))is the Fermi distribution function, and εs(k)=sħvFkεFis the graphene energy. The band overlap of the wavefunctions

fss(k,q)=12(1+ssk+qcosφ|k+q|)E59

is a specific property of graphene, where φ signifies the angle between k and q. At zero temperature, the Fermi-Dirac distribution functions are step functions. In that case, using ωnω+iδ, Eq. (58) is simplified to

χT=00±(q,ω)=gsgv4π2ħd2kα=±αf±(k,q)ω+αvF(k|k+q|)+iδ,E60

where the + and – signs denote the intra and interband transitions, respectively. Integration over φ and k gives the retarded polarizability or charge-charge correlation function

χ0(q,ω)=χεF=00(q,ω)+Δχ0(q,ω),E61

where

χεF=00(q,ω)=iπI(q,ω)ħ2νF2,E62

and

Δχ0(q,ω)=gεF2πħ2vF2+I(q,ω)ħ2vF2{(ħω+2εFħvFq)Θ(2εFħωħvFq1)×((2εFħωħvFq)iπ)Θ(ħω2εFħvFq+1)×(ħω2εFħvFq)}E63

The two functions I(q,ω)and (x)are defined by

I(q,ω)=g16πħvF2q2ω2νF2q2E64

and

(x)=xx21ln(x+x21),E65

Where g=gsgv=4. For ω>qvFand in the long wavelength limit q → 0, x=|(ħω±2εF)/ħvFq|1, thus x21x2and (x)x22ln(x). We derive here the dynamical polarizability in Eq. (63) in these regimes. The expansion of I(q,ω)gives

I(q,ω)=g16πħvF2q2ω(1vF2q2ω2)1/2g16πħvF2q2ω(1+vF2q22ω2).E66

In this condition and for intraband transitions ħω<2εF( ħω<2μ)

(ħω+2εFħvFq)(2εFħωħvFq)=8ħωεFħ2vF2q2+2ln|2εFħω2εF+ħω|.E67

As a result, Δχ0(q,ω)reduces to

Δχ0(q,ω)=gεF2πħ2vF2+I(q,ω)ħ2vF2{8ħωεFħ2vF2q2+2ln|2εFħω2εF+ħω|+iπ}=gq28πħω{2εFħω+12ln|2εFħω2εF+ħω|+iπ2}.E68

If 2εFħω, then

Δχ0(q,ω)=gq28πħω{2εFħω+iπ2}.E69

By taking the decay rate ω ω+iτ1into account and inserting Eq. (62) into Eq. (61), the dynamical polarizability reduces to

χ0(q,ω)εFq2πħ2(ω+iτ1)2.E70

In the presence of optical phonons, the effective dielectric function in the RPA regime is given by Eq. (4), which we restate here for convenience [98, 99]:

εRPA(q,ω)=εmvc(q)χ0(q,ω)εmlvsph,l(q,ω)χ0(q,ω)εmvoph(q,ω)χj,j0(q,ω),E71

where εm=(ε1+ε2)/2is the average of dielectric constants of graphene’s environment. The collective oscillation modes of the electrons can be obtained by solving εRPA(q,ω)=0. The extinction function is identified by Z=δTT0, or for the plasmonic structure coupled to an optical cavity Z=δRR0, where δR=RR0and R/R0 is the reflectance with/without plasmon excitation, which corresponds to the enhanced absorbance at resonance frequencies [129, 131]

ZIm(1εRPA).E72

The loss function represents the amount of energy dissipated by exciting the plasmon coupled to the substrate and graphene optical phonons. The collective oscillation modes of the electrons can be obtained by solving εRPA(q,ω)=0. Considering the first two terms in Eq. (71) gives the plasmon dispersion relation of graphene

ħωp(q)=gαGεFħvFq2ε0,E73

where the graphene fine structure constant is given by

αGe2εħvF=cεvFα,E74

with α=e2/ħcbeing the fine structure constant. In the case of graphene on a SiO2 substrate, εr1=1for air and εr2=3.9for SiO2, which yield αG=0.9. For suspended graphene, εr1=εr2=1and therefore αG=2.2. Thus, it is possible to tune αG by altering the dielectric materials surrounding the graphene sheet [141].

It is also interesting to note that in contrast to the εF3/4n01/2dependence of the bulk plasmon resonance frequency in Eq. (20) and the εF1/2ns1/2dependence of the 2DEG plasmon, the graphene plasmon exhibits a εF1/2ns1/4dependence.

In order to show the results for εF>0, we define two regimes for (x), i.e.,

>(x)=xx21cosh1(x),for x>1<(x)=x1x2cosh1(x),for x<1E75

Using Eq. (60), the imaginary part of dynamical polarization is given by [140]

ImΔχ0(q,ω)=I(q,ω)ħ2vF2×{>(2εFħωħvFq)>(2εF+ħωħvFq),1Aπ1B>(2εF+ħωħvFq),2A<(ħω2εFħvFq),2B0,3A0.3BE76

and the real part is obtained by the Kramers-Kronig relation, yielding

ReΔχ0(q,ω)=gεF2πħ2vF2+I(q,ω)ħ2vF2×{π,1A>(2εFħωħvFq)+>(2εF+ħωħvFq),1B<(ħω2εFħvFq),2A>(2εF+ħωħvFq),2B<(ħω2εFħvFq)+<(2εF+ħωħvFq),3A>(2εF+ħωħvFq)>(ħω2εFħvFq).3BE77

Different regions are shown in Figure 13. As we mentioned, the plasmon dispersion relation is determined by ε(q,ωpiγ)=0, where γ is the decay rate of the plasmon. For weak damping, the plasmon dispersion relation and damping are given by

Figure 13.

Different regions related to the dynamical polarization. The regions are separated by the straight line ω=vFq (solid), ħω=ħvFq−2εF (dashed), and ħω=2εF−ħvFq (dotted). The red-shaded area depicts the region of interband Landau damping while the green-shaded area marks the region of intraband Landau damping. The radiative interband Landau damping originating from direct (vertical) transitions occurs at q = 0. The red-shaded area with q≠0 corresponds to the region of nonradiative interband Landau damping.

ε0=vcReΔχ0(q,ωp),γ=ImΔχ0(q,ωp)(/ω)ReΔχ0(q,ω)|ωp.E78

The solution of the first part of Eq. (78) exists only for ReΔχ0(q,ωp)>0, which is valid only for finite graphene doping and ω>vFq. Interestingly, a plasmon does not decay if ImΔχ0(q,ωp)=0, which is the case in region 1B shown in Figure 13. For the finite doping (εF>0), the acoustic phonon at long wavelength is inside the region 1A in Figure 13.

7. Surface plasmon resonances in 3D topological insulators

So far, we have discussed nanoplasmonics in 3D metals and graphene and their damping mechanisms. Now we focus on the RPA theory of nanoplasmonics in 3D topological insulator (TI) materials. In particular, we are going to identify the graphene plasmons as a limiting case of 3DTI plasmons in the case when the thickness of the 3DTI slab becomes atomically thin. Bismuth selenide (Bi2Se3) is the prime example of a 3DTI material that has a rhombohedral crystal structure and consists of five atomic layers arranged along the z-direction, known as quintuple layers (QLs). Two such QLs are coupled by a weak interaction, predominantly of the van der Waals type. At the vicinity of the Γ point in the Brillouin zone of Bi2Se3, the low energy of Dirac particles in the xy-plane can be described by the effective Hamiltonian [142, 143], HTI=ħvf(σ×kz^)EF, where z is the unit vector in a z-direction. In the case of graphene, it is the pseudospin, representing the two sub-lattice indices, which is locked to the momentum of the electron. The chiral nature of the electronic states with the spin locked in perpendicular direction to the momentum results in a unique type of collective excitations which are, indeed, accompanied by transverse spin fluctuations [142, 143].

Let us consider a slab of a 3DTI material of thickness d > 5 nm, which is large enough to suppress any overlap of the single-electron states between the top and the bottom layers. In contrast, long-range Coulomb interaction exists and couples the opposite surfaces as in 2D electron plasmas.

We consider a Hamiltonian that describes the properties of collective oscillations on the top and the bottom surfaces of a 3DTI material as [142]

Hl/l'=ħvfα,β,k,l,l'ak,l,α+(τ×k)·z^ak,l,β+12Sq,l,l'υl,l'(q)nq,lnq,l',E79

where l(l')indicates the top (bottom) surface, α, βare spin indices, a is the annihilation operator, and S is the surface area. nq,l=k,αakq,l,α+ak,l,αis the density operator, and τ={τx,τy}are the Pauli matrices. With a long-range Coulomb interaction within the RPA, the response function can be written as

χRPA2x2=χ02x2(q,ω)1υ(q)χ02x2(q,ω),E80

where χ02x2(q,ω)=(χ1100χ22)is a 2 × 2 matrix. υ(q)is a 2 × 2 matrix whose diagonal and off diagonal elements are the two-dimensional Fourier transforms of Coulomb potentials and can be obtained by solving Poisson’s equation [144]. Similar to the RPA dielectric function for graphene shown in Eq. (71), the equivalent equation for the RPA dielectric function in 3DTIs is

εRPA2x2(q,ω)=1v(q)χ02x2(q,ω).E81

Compared to graphene, the main difference is that υ(q)is a 2 × 2 matrix accounting for the intrasurface and intersurface interactions. Here we write these expressions as

υ11/22(q)=4πe2ε2(cosh(qd)+(ε3/1/ε2)sinh(qd))q(ε2(ε1+ε3)cosh(qd)+(ε1ε3+ε2)sinh(qd))E82

for the intrasurface interaction, and

υ12(q)=4πe2ε2q(ε2(ε1+ε3)cosh(qd)+(ε1ε3+ε2)sinh(qd))E83

for the intersurface interaction. εi (i = 1, 2, 3) are the dielectric functions for the top, middle, and bottom layers, respectively. In the limit of equal dielectric constant, ε1=ε2=ε3=ε, V11=V222π2e2/εq, and V12=V21(2π2e2/εq)eqd. This shows that for a given wave vector the intrasurface Coulomb interaction depends on the dielectric constant of the materials, the intersurface Coulomb interaction not only depends on the dielectric constant, but is also modulated by the negative distance-dependent exponential factor. For a sufficiently thick slab, i.e., qd, the top and the bottom surfaces are completely decoupled and therefore we obtain the limit of V12=V210, and thus we only have a nonzero intrasurface potential. This is the case that is equivalent to the Coulomb interaction in a sheet of graphene.

The response function χ provides important information about the collective states that are excited at small transferred momentum. The collective mode frequencies of the system can be obtained by solving Det[1υ(q)χ02x2(q)]=0. In the region vfq<ω<2Ef/ħvfq, collective modes of oscillations are undamped. Beyond that limit, such modes are not observed because the energy of the modes is transferred to the particle-hole excitations. In Section 6, we derived the dynamical polarizability for the case of a graphene sheet. For vfq<ωwith Im[χ]=0, the linear response function can be written as [144]

χnnl/l'(q,ω)=EFl/l'2πħ2υF2+116πħq2ω2(υFq)2×[G(2EFl/l'+ħωħυFq)G(2EFl/l'ħωħυFq)],E84

where G(x)=xx21cosh1(x). Using the relation for the relative permittivity εr=1+χ, the linear response function given by Eq. (84) is obtained from Eq. (61) in the long wavelength limit. This charge-charge response function describes explicitly the dispersion of the collective modes exist on the two sides of a slab geometry of a 3DTI slab. Using Eq. (80), the pole of Det[1υ(q)χ02x2(q)]=0can be solved with the help of Eq. (35). The potential υ(q) is a 2×2 block diagonal (corresponding to spin and charge) for two surfaces. In the limit of an infinitely thick slab, we recover the intrasurface potential for υ(q)and the corresponding mode for graphene. The two modes in the long wavelength limit are: the symmetric (photon-like) mode and the antisymmetric (plasmon-like) mode. At a thickness below the Thomas-Fermi screening length, the symmetric mode is highly Landau damped. In Figure 14, we show the antisymmetric mode for a slab of thickness 120 nm at different Fermi energies of 0.3, 0.2, and 0.1 eV, taking into account the inter and intrasurface potentials given in Eqs. (82) and (83) in the q → 0 limit of the RPA response function. This result can be compared with the experimental observation of the Dirac plasmon dispersion in Bi2Se3 [145].

Figure 14.

Plasmon dispersion in Bi2Se3 slab of thickness 120 nm at different Fermi energy with top and bottom dielectric constant of 1 and 9.3, respectively, taking full potential and q→0 limit. Only the antisymmetric mode is shown here.

In the limit of small thickness of the slab, when the condition qd → 0 is satisfied, the antisymmetric photon-like mode behaves as

ωant2=vf2kfαDε1+ε3[(1+ξ2)1/2+(1ξ2)1/2]qE85

and the symmetric plasmon-like mode as

ωsym2=vf2kfαD2ε2[(1ξ)1/2(1+ξ)1/2(1ξ)1/2+(1+ξ)1/2]dq2,E86

where ξ=(n1n2)/(n1+n2)is the density polarization and kf=4π(n1+n2)is the total Fermi wave vector. αD=e2/ħvfis the fine structure constant for the Dirac system. It is to be noted that indices 1 and 2 indicate two different surfaces. For equal Fermi energy at the top and the bottom surfaces, the density polarization reduces to zero, resulting in an equilibrium situation. When the thickness of the slab tends to zero, we recover only the antisymmetric mode, which is the mode obtained for a sheet of graphene. As the slab thickness goes to infinity, the two surfaces interact weakly and the intersurface potential falls rapidly as the thickness increases. In this case, we obtain that the two antisymmetric modes each correspond to a single sheet of graphene. Using a series expansion of the frequency, it is possible to obtain a more accurate solution for the symmetric mode, which is derived elsewhere [144]. The solution is given by

ωsym=(2ε2+αDkfd)vfq/(2ε2(ε2+2αDkfd))E87

For a size of the slab smaller than the mode wavelength, i.e., d<λp, the antisymmetric mode is independent of thickness, and the symmetric mode has a d2 dependence. The plasmon dispersions for the limit of qd → 0 for symmetric and antisymmetric modes with 10 and 20 nm of slab width at different dielectric environments are shown in Figure 15a. The electric field profiles along (Ex) and perpendicular (Ez) to the plasmon polarization direction are shown in Figure 15b for a thin slab of thickness 10 nm with ε1=ε3=3.8for an antisymmetric mode. We also show the Landau damping regions with a blue shade given by the line segments ωυFqand ω2EFυFq, in complete analogy to Figure 13. The distribution of excited “hot” carriers is limited to a narrow region around the Fermi level. In quasi-equilibrium, the hot carrier distribution can affect the plasmon lifetime.

Figure 15.

(a) Symmetric (dotted) and antisymmetric (solid) modes of plasmon dispersion in 3DTI in the limit of qd→0 for the slab thickness of 10 and 20 nm (Eqs. (7) and (9)). The Landau damping region is shown with the blue shade, determined by ω≤υFq and ω≥2EF−υFq. (b) The electric field profile of surface plasmons along and perpendicular direction of polarization for antisymmetric mode (plasmon-like mode) in a slab of thickness 10 nm.

The topological surface states are extended inside the bulk with a localization length given by lo=ħvf/Egwhere Eg is the bulk bandgap. lo is equal to 1.5 nm (~2 quintuple layer (QL)) in Bi2Se3. Due to the quantum confinement of the carriers, two different conductivities arise: one due to the topological surface states that extend to 2 QL and the other one due to the 2DEG trapped below the surface, which extends up to 8 QL [146]. Theoretically, it should be possible to separate out the two different effects and write the response functions as χ=χTI+χ2DEGfor two independent channels. The contribution of the 2DEG to the plasmon dispersion can be calculated separately [147].

Surface plasmons in 3DTI surfaces do not consist only of charge density waves but are also accompanied by spin density waves due to inherent nature of spin-momentum locking. This can be qualitatively understood by calculating the surface current in terms of the spin and charge quantity by J=vfσ×z^. The continuity equation nq/t=Jshows that the charge density nq oscillates in phase with the transverse spin sT=z^(q^×s), and the amplitude of their ratio is given by sT/nq=1/ħvf[EFαDvf/2q]1/2, which can be much larger than unity in the long wavelength limit [142]. To understand the electromagnetic response of spin and charge oscillations in a selfsustained mode, we write the induced potential Vind=υ(q)χ04x4(q,ω)Vindunder vanishing external perturbation [147]. One can solve again for the pole (1υ(q)χ04x4(q,ω))Vind=0in order to obtain the explicit expression for the induced potential. This time, both the spin and charge part of the response functions are included, i.e., the solutions include the explicit information of the spin and charge oscillations. The solutions to the induced potential proportional to

V±ind(10y±0),E88

where y+=1 and y=EFlEFl'. The induced charge density is calculated from the response function ρ=χ04x4(q,ω)Vind. Spins are inherently attached to the momentum as accounted by the complete response function given by Eq. (80). The symmetric and antisymmetric modes are purely spin-like and purely charge-like for a slab of 3DTI. This can be seen in the following expressions in the limit of qd0as

(φ1φ2)=ρ(1αD(EFl+EFl')ħυFq(εl+εl')EFlEFl'αD(EFl+EFl')ħυFq(εl+εl')EFlEFl')E89

for the antisymmetric mode, and

(φ1φ2)=ρ(1αDd(EFl×EFl')ħυFqεTI(EFl+EFl')1αDd(EFl×EFl')ħυFqεTI(EFl+EFl'))E90

for the symmetric mode, where φ1=(ρlsl)and φ2=(ρl'sl'), and the indices l and l′ can be chosen to represent the top and the bottom surfaces. Clearly, charges oscillate in phase and spins oscillate out of phase for antisymmetric mode, whereas spins oscillate in phase and charges oscillate out of phase for the symmetric mode on the top and the bottom surfaces.

It has been observed that the phonon modes in Bi2Se3 have frequencies of around 2 THz. Since the phonon modes overlap with the plasmon modes in energy space, there is a possibility of mode hybridization, similar to the case of graphene. In a micro-ring structure of Bi2Se3, there are bonding (lower) and antibonding (upper) plasmon modes [148]. The antibonding plasmon mode frequency is much larger than the phonon mode frequency, and the bonding plasmon mode overlaps with the phonon mode, resulting in hybridization in two branches with an interaction frequency of around 0.35 THz.

Recently, in Ref. [149], authors investigated the possibility of obtaining SPPs in Weyl semimetals. In Weyl semimetals such as TaAs, NbAs, YbMnBi2, and Eu2Ir2O7 [150152] the valence and the conduction band touch in isolated points of the Brillouin zone close to the chemical potential, and their dispersion is described by an equation similar to the one for Dirac metals. A pair of Weyl nodes appears with opposite chiralities with a distance of b in the reciprocal space. Their topological properties are described by Sθ=e24πħcdtd3r θ E.B, where θ=2(brbot). Sθis an important parameter that alters the electrodynamics of the bulk response to an applied electric field. In such a system, the authors found that the SPPs with b0have properties similar to those in ordinary metal SPPs in the presence of an external magnetic field, despite their unique origin. Therefore, the θ term that contributes to the dielectric tensor originating from the anomalous Hall displacement current induces an anomalous surface magnetoplasmon. Interestingly, for certain orientations of the surface of the Weyl semimetal, the authors predict the nonreciprocal dispersion of SPPs, i.e., the propagation of the SPP is unidirectional.

8. Conclusion

In this book chapter, we present descriptions of the plasmonic properties of metal nanostructures of different geometry, their size dependence, and applications in modern nanotechnology. We show dynamic control over the plasmon resonances where a shift up to 125 nm at a density of 1022 cm−3 can be achieved using a pump-probe technique. This provides the opportunity to utilize plasmonics in modern information processing devices. In addition to plasmonics in 3D metal nanostructures, we present a description of graphene and 3DTI plasmonics using classical and quantum perspectives. Using the RPA theory, we obtain symmetric and antisymmetric modes in a slab of 3DTI, which reduces to a graphene plasmonics in the limit of zero thickness of the slab. Surface plasmon (SP) damping mechanisms are interesting due to their potential applications for enhanced current density that comes from SPs nonradiative damping in nanostructures smaller than the skin depth. We present a quantum theory of SPs damping in metals and layered materials like graphene.

There are several potential applications of the nanoplasmonics in graphene and 3DTI. Graphene and 3DTI are potential candidates for nanospasers that utilize Dirac fermions, unlike the massive electrons or holes in the originally proposed spasing scheme by Bergman and Stockman in 2003 [153]. The spaser is a nanoplasmonic counterpart of a laser, where photons are not emitted. In Ref. [154], the authors have proposed a scheme of nanospasing using a sheet of graphene with an electrically pumped cascaded quantum well structure working as a gain medium. For the range of Fermi energy, it was shown that spasing could be potentially obtained for a plasmon relaxation time of several femtoseconds in the mid-IR range. A similar scheme with an optically pumped nanospaser for a slab 10 nm of 3DTI was proposed in Ref. [155]. It is advantageous to use a 3DTI, such as Bi2S3, as a nanospaser due to the possibility of using its bulk as a gain medium and the surface as a medium that supports SPPs. This configuration avoids the use of a separate gain medium to provide the feedback for the SPPs. Therefore, a 3DTI nanospaser can be truly nanoscopic and can be used for various applications in physics, chemistry, and biology.

Acknowledgments

We acknowledge the support provided by NSF grant CCF-1514089 and DARPA grant HR0011-16-1-0003.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Hari P. Paudel, Alireza Safaei and Michael N. Leuenberger (June 21st 2017). Nanoplasmonics in Metallic Nanostructures and Dirac Systems, Nanoplasmonics - Fundamentals and Applications, Gregory Barbillon, IntechOpen, DOI: 10.5772/67689. Available from:

chapter statistics

1085total chapter downloads

1Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Ultrafast Nonlinear Optical Effects of Metal Nanoparticles Composites

By Kwang-Hyon Kim, Anton Husakou and Joachim Herrmann

Related Book

First chapter

Thermodynamic Properties of Nano-Silver and Alloy Particles

By Wangyu Hu, Shifang Xiao, Huiqiu Deng, Wenhua Luo and Lei Deng

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us