Localized Surface Plasmons in Nanostructured ... - ACS Publications

Desalegn T. Debu , Stephen J. Bauman , David French , Hugh O. H. Churchill .... Sajid M. Choudhury , Di Wang , Krishnakali Chaudhuri , Clayton DeVault...
0 downloads 0 Views 491KB Size
Subscriber access provided by University of Sussex Library

Communication

Localized Surface Plasmons in Nanostructured Monolayer Black Phosphorus Zizhuo Liu, and Koray Aydin Nano Lett., Just Accepted Manuscript • DOI: 10.1021/acs.nanolett.5b05166 • Publication Date (Web): 06 May 2016 Downloaded from http://pubs.acs.org on May 9, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Nano Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Localized Surface Plasmons in Nanostructured Monolayer Black Phosphorus Zizhuo Liu and Koray Aydin∗ Department of Electrical Engineering and Computer Science, Northwestern University, Evanston, IL 60208, United States E-mail: [email protected]

Abstract Plasmonic materials provide electric-field localization and light confinement at subwavelength scales due to strong light-matter interaction around resonance frequencies. Graphene has been recently studied as an atomically thin plasmonic material for infrared and terahertz wavelengths. Here, we theoretically investigate localized surface plasmon resonances (LSPR) in a monolayer, nanostructured black phosphorus (BP). Using finite-difference time-domain simulations, we demonstrate LSPRs at mid-infrared and far-infrared wavelength regime in BP nanoribbon and nanopatch arrays. Due to strong anisotropic in-plane properties of black phosphorus emerging from its puckered crystal structure, black phosphorus nanostructures provide polarization dependent, anisotropic plasmonic response. Electromagnetic simulations reveal that monolayer black phosphorus nanostructures can strongly confine infrared radiation in an atomically thin material. Black phosphorus can find use as a highly anisotropic plasmonic devices. ∗

To whom correspondence should be addressed

1

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Keywords Black Phosphorus, Plasmonics, LSPR, 2D Materials, Anisotropy

Introduction Atomically thin two-dimensional (2D) materials have received burgeoning amount of interest in recent years due to exciting physical properties. They exhibit interesting electronic, optical, mechanical and thermal properties due to their mono to few-layer thicknesses. 1–4 2D materials offer a new class of platform to achieve novel electronic and photonic properties in ultra-compact sizes. Although such materials have inherently monolayer to few-layer thicknesses, it has been showed that 2D materials can interact strongly with the incident light. In particular, graphene emerged as a monolayer plasmonic platform. 5–8 Theoretical and experimental studies confirmed that both propagating 9,10 and localized 11–13 surface plasmon modes can be excited in an atomically thin, nanostructured graphene. Until the discovery of the graphene as a plasmonic material, noble metals such as silver and gold have been the material of choice both for fundamental and applied research. However, 2D plasmonic materials provide a unique opportunity by confining plasmons in an extremely thin optical material and also represent a great challenge for light-matter interactions due to their inherent atomically thin thickness. With such highly localized electric fields in 2D plasmonics, it is possible to build and integrate devices in smaller sizes that exhibit novel electronic and optical properties compared with traditional bulk plasmonic materials. By patterning periodic nanostructures in 2D materials, resonances can be observed optically by measuring the extinction spectra. 14,15 Black phosphorus (BP) is a recently emerging 2D layered material that can be exfoliated to few layers and monolayer. 16–23 In a monolayer BP, the phosphorus atoms form a hexagonal lattice with a puckered structure resulting in in-plane anisotropic properties. Black 2

ACS Paragon Plus Environment

Page 2 of 16

Page 3 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

phosphorus has been recently investigated for potential applications including field effect transistors, 24–27 heterojunction p-n diode, 28 photovoltaic devices, 29 and photodetectors. 30 The Raman scattering spectra and photoluminescence measurements from BP revealed its highly anisotropic properties. 31 A recent theoretical paper by Low et al. reports collective plasmonic excitations in BP. 32 It has been shown that mono to few-layer BP supports anisotropic plasmonic dispersion due to different effective mass along different crystal directions. Here, we propose and numerically demonstrate that localized surface plasmons can be excited in nanostructured monolayer black phosphorus. We present the study of the electromagnetic response of periodically patterned black phosphorus sheet. We consider the parameters varied in a wide range and the different absorption spectrum is compared due to the inherent anisotropy in BP. In particular, we analyze the confinement of the electromagnetic field in both BP nanoribbon and nanopatch arrays. We demonstrate the anisotropic behavior by replacing BP nanoribbons in x direction (armchair direction) and y direction (zigzag direction). The performance of the structures may allow for the realization of new plasmonic devices which will take advantage of the directional dependence of its plasmon properties.

Results and discussion The photonic properties of a monolayer BP can be described by employing a simple semiclassical Drude model. The conductivity 32 is given as: iDj σjj = π(ω + iη¯h )

πe2 n Dj = mj

(1)

where j denotes the direction concerned and Dj is the Drude weight. A similar expression has been utilized for graphene 33 with D =

µe2 h2 ¯

. The electron mass along the x direction and

3

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 16

y direction 32 can be described by:

mcx =

h ¯2 2γ 2 ∆

mcy =

+ ηc

h ¯2 2νc

(2)

The parameters are determined by fitting the known anisotropic mass. For monolayer BP, 32 we have γ =

4a π

eVm , ∆ = 2 eV, ηc =

h2 ¯ , 0.4m0

νc =

h2 ¯ . 1.4m0

We choose the electron doping

n = 1013 cm−2 as well as η = 10 meV to describe the relaxation rate. a is the scale length of the BP and

π a

is the width of the Brillouin Zone.

The dielectric function 34 for a thin film reads(with 2D conductivity):

ǫjj = ǫr +

iσjj ǫ0 ωa

(3)

For a monolayer BP, the relative permittivity 35 is given as ǫr = 5.76. In the finite-difference time-domain (FDTD) simulations, BP film thickness (a) is chosen to be 1 nm. This audacious assumption will speed up the simulation time and resemble the real results as long as the simulation mesh is fine enough in the BP film. In our simulation setup, we choose the mesh to be 0.25 nm. The rightness of using finite thickness to simulate 2D film and more details about the simulation are discussed in the Supporting Information. In our designs, periodic monolayer BP nanoribbons are placed on a transparent insulator and an optically thick gold mirror. The metallic mirror is used to reflect light and suppress the transmission. Compared with free standing BP nanoribbons, the dielectric layer between the BP and gold mirror forms a Fabry-Perot cavity, therefore increases the interaction of light with monolayer BP. Dielectric thickness can be chosen to maximize the absorption using Fabry-Perot interference (see Supporting Information). In our simulations, the dielectric thickness is chosen to be 5 µm. We performed full-field electromagnetic simulations with the wavelengths between 20 - 80 µm. For the optical constants in the simulation, we used a constant, non-dispersive refractive index of n = 1.7 for the dielectric spacer and perfect electric conductor model for the metallic mirror. 4

ACS Paragon Plus Environment

Page 5 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Figure 1: (a) Schematics of periodically patterned BP nanoribbons. (b) Anisotropic absorption spectrum can be observed in BP ribbon arrays along x and y direction. The arrows in the inset indicate the polarization of the electric field. (c) Calculated electric field amplitude, Ez for one BP nanoribbon unit cell(along x direction). (d) Calculated total electric field intensity (|E|2 ) for one BP nanoribbon unit cell(along x direction). Electromagnetic field is localized around the edges of the nanoribbon.

5

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Localized plasmon resonances can be excited with an incident light when electric field is perpendicular to the nanoribbons (Fig. 1a). Since metallic mirror suppresses the transmission, absorptivity (A) can be calculated with a simple formula A = 1 - R, where R is the reflectivity from the simulated structure. Since BP exhibits anisotropic optical and electronic properties along x and y directions, we anticipate polarization-dependent localized plasmon resonance behavior for monolayer BP nanoribbons along the armchair and zigzag directions. In order to investigate the effect of anisotropic optical properties on the plasmonic behavior of the BP, we performed electromagnetic (EM) simulations for two different nanoribbon structures as shown in Fig. 1b. The structural parameters for both x (armchair) and y (zigzag) directions are chosen to be the same. Having the period p = 250 nm and the monolayer BP ribbon width w = 150 nm, the absorption spectrum is plotted in Fig. 1b. For both cases, an absorption peak has been observed for s polarized light. Due to the mass anisotropy, the smaller mass along x indicates higher resonance frequency therefore shorter wavelength. Compared with the zigzag direction (red line), the absorption peak in armchair direction is higher and the resonance linewidth is narrower. This result suggests that the black phosphorus is optically more lossy for the zigzag direction at the corresponding resonance wavelength. As mentioned before, since graphene has also been described with the similar Drude model but with different Drude weight, the plasmons in BP naonribbons also behave differently in scales compared with graphene nanoribbons (see Supporting Information). The side-view electric field profile of the vertical component (Ez ) shown in Fig. 1c is plotted with the case of the nanoribbon arrays along x direction at the resonance wavelength (31.5 µm). Positive and negative dipoles can be easily seen indicating a localized plasmon resonance behavior. The edges of the BP nanoribbon distort the in-plane component of the electric field. And for the edge surface plasmon modes, the electric field enhancement is mostly pronounced around the edges, which resemble the edge modes in graphene ribbon structures. 36 Fig. 1d and the supporting video confirm that the total electric field intensity is localized around the two

6

ACS Paragon Plus Environment

Page 6 of 16

Page 7 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

edges of the BP nanoribbon.

Figure 2: (a,b) Schematics and perspective side view of the BP ribbon arrays along x and y direction separately. (c,d) Absorption map of two directions with fixed period 250 nm and various widths. (e,f) Absorption spectra for various widths of two directions. We have also performed additional simulations for different nanoribbon widths as well as different periodicities in order to understand the effect of the nanoribbon array parameters on the observed plasmon resonance behavior. In Fig. 2, we summarize our results for monolayer BP nanoribbon arrays with various widths. In these simulations, the periodicity of the arrays is kept at 250 nm, as shown from the side-view in Fig. 2a and Fig. 2b for ribbons along both x and y directions. The width is changed from 0 nm (no BP) to 250 nm (continuous BP). 7

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Absorption spectra as a function of wavelength and varying widths for nanoribbon arrays along x and y directions are plotted in Fig. 2c and Fig. 2d, respectively. It is clear that the absorption peak position, therefore the resonance wavelength increases to longer wavelengths with increasing nanoribbon width, which follow the scaling behavior of localized plasmons. The wider BP nanoribbons support the electromagnetic oscillations at longer wavelengths. The shift trends are similar for both armchair and zigzag directions, but the absorption peak is much larger for armchair direction (x). In Fig. 2e and Fig. 2f, we plot the absorption spectra for different nanoribbon widths along x and y directions separately. It can be seen that with increasing width, the resonance peaks become broader due to the increasing optical losses at longer wavelengths. One interesting property we note is that there is an absorption maximum of 0.47 for 190 nm width nanoribbon with a periodicity of 250 nm. Additional simulations are performed for different periodicities changing from p = 150 nm to 500 nm for fixed nanoribbon width of w = 150 nm. Results are summarized in Fig. 3. For small periodicities (periodicity is closer to the nanoribbon width), there is a strong coupling between neighboring nanoribbons resulting in a change for the resonance wavelength positions. The absorption increases for reduced periodicities mainly due to the smaller gap between the nanoribbons and stronger inter-ribbon coupling. However, for larger periodicities, where the interaction between neighboring BP nanoribbons are weak, the position of the absorption peak shifts to shorter wavelengths. And the total absorption decreases due to the smaller volume of the BP in certain areas. For the large period case, the resonance almost stays at the same wavelength as expected for a localized plasmon resonance, and the plasmon wave vector follows the trend q ∼

π 37 . w

Previous results in graphene ribbons have shown the

existence of phase shift of the plasmons upon reflection, 13,38 which might also occur in BP nanoribbons. More discussions could be found in the the Supporting Information. Inspired by the anisotropic plasmonic response of the BP nanoribbon arrays along x and y direction, we can construct BP nanopatch arrays (square shaped isolated BP nanostructures) as shown in Fig. 4a. For consistency and ease of comparison with earlier discussions, we

8

ACS Paragon Plus Environment

Page 8 of 16

Page 9 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Figure 3: (a,b) Absorption map of two directions with fixed width 150 nm and various periods. (c,d) Absorption spectra for various periods of two directions.

9

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

choose the nanopatch array parameters as p = 250 nm and w = 150 nm in both x and y directions. When light is polarized along the x direction (0◦ ), we observe a high absorption resonance peak at shorter wavelength. For y-polarization, i.e., the polarization angle is 90◦ , the resonance shifts to longer wavelengths as expected due to higher effective mass along the zigzag (y) direction. Compared with Fig. 1b, the resonance wavelengths are kept at the same position as expected. But the absorption intensity is reduced which can be explained due to the discontinuity of the ribbons in the case of nanopatch arrays resulting in lower BP volume. Interestingly, for polarization angle of 45◦ , localized plasmon resonances for armchair and zigzag directions can be excited simultaneously. Absorption intensity reduces at respective plasmon resonance wavelengths due to reduced light intensity at corresponding x and y directions. In order to understand the resonance behavior for monolayer BP nanopatch arrays, we calculated the amplitude and electric field intensity using EM simulations. Fig. 4c plots the amplitude of the Ez for x-polarization (0◦ ) right above the BP film and clearly illustrates the electric dipole behavior expected in structures exhibiting localized plasmon resonances. Total electric field intensity profile shown in Fig. 4d implies that the electric field is mostly localized at the left and right edges corresponding to the polarization of the light. The electric field intensity is highest around four corners of the square due to the sharpness in shape. In conclusion, we propose and numerically demonstrate localized surface plasmons modes in monolayer BP using nanoribbon and nanopatch arrays. Anisotropic behavior in the absorption spectra are predicted when exciting BP plasmons along two different directions (armchair and zigzag). Various structure parameters are examined to analyze the material properties and the localized nature. Furthermore, we find that the electric field is mostly localized around the edges of the structure. It is possible to design a symmetric metasurface with BP and fulfill the polarization conversion due to the anisotropic absorption. And BP could also be used to build a polarization sensor at IR wavelengths even with symmetrical

10

ACS Paragon Plus Environment

Page 10 of 16

Page 11 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Figure 4: (a) Schematics of periodically patterned BP nanopatches. (b) Anisotropic absorption spectrum can be observed when applying light with different angles. (c) Calculated electric field amplitude, Ez for one BP nanopatch unit cell(along x direction). (d) Calculated total electric field intensity (|E|2 ) for one BP nanopatch unit cell(along x direction).

11

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

structures. Although it is still a big challenge for current research to obtain large monolayer BP film with either exfoliation method or chemical vapor deposition method, our proposed structures provide the guidelines for future experimental research on new 2D anisotropic plasmonics with BP and directional dependent plasmonic devices.

Acknowledgement This material is based upon work supported by the Materials Research Science and Engineering Center (NSF-MRSEC) (DMR-1121262) of Northwestern University. We also acknowledge partial support from the Institute for Sustainability and Energy at Northwestern (ISEN) through ISEN Booster Award.

Notes The authors declare no competing financial interests.

Supporting Information Available Details of the simulation method and data, discussions about dielectric thickness as well as additional figures and video for the electrical field can be found in supporting material. This material is available free of charge via the Internet at http://pubs.acs.org/.

References (1) Xia, F.; Wang, H.; Xiao, D.; Dubey, M.; Ramasubramaniam, A. Nat. Photonics 2014, 8, 899–907. (2) Koppens, F.; Mueller, T.; Avouris, P.; Ferrari, A.; Vitiello, M.; Polini, M. Nat. Nanotechnol. 2014, 9, 780–793. (3) Butler, S. Z.; Hollen, S. M.; Cao, L.; Cui, Y.; Gupta, J. A.; Gutiérrez, H. R.; Heinz, T. F.; Hong, S. S.; Huang, J.; Ismach, A. F. et al. ACS Nano 2013, 7, 2898–2926. 12

ACS Paragon Plus Environment

Page 12 of 16

Page 13 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

(4) Wang, Q. H.; Kalantar-Zadeh, K.; Kis, A.; Coleman, J. N.; Strano, M. S. Nat. Nanotechnol. 2012, 7, 699–712. (5) Grigorenko, A.; Polini, M.; Novoselov, K. Nat. Photonics 2012, 6, 749–758. (6) Low, T.; Avouris, P. ACS Nano 2014, 8, 1086–1101. (7) Koppens, F. H.; Chang, D. E.; Garcia de Abajo, F. J. Nano Lett. 2011, 11, 3370–3377. (8) Bludov, Y. V.; Ferreira, A.; Peres, N.; Vasilevskiy, M. Int. J. Mod. Phys. B 2013, 27, 1341001. (9) Gao, W.; Shu, J.; Qiu, C.; Xu, Q. ACS Nano 2012, 6, 7806–7813. (10) Nikitin, A. Y.; Alonso-González, P.; Hillenbrand, R. Nano Lett. 2014, 14, 2896–2901. (11) Ju, L.; Geng, B.; Horng, J.; Girit, C.; Martin, M.; Hao, Z.; Bechtel, H. A.; Liang, X.; Zettl, A.; Shen, Y. R. et al. Nat. Nanotechnol. 2011, 6, 630–634. (12) Chen, J.; Badioli, M.; Alonso-González, P.; Thongrattanasiri, S.; Huth, F.; Osmond, J.; Spasenović, M.; Centeno, A.; Pesquera, A.; Godignon, P. et al. Nature 2012, 487, 77– 81. (13) Brar, V. W.; Jang, M. S.; Sherrott, M.; Lopez, J. J.; Atwater, H. A. Nano Lett. 2013, 13, 2541–2547. (14) Yan, H.; Low, T.; Zhu, W.; Wu, Y.; Freitag, M.; Li, X.; Guinea, F.; Avouris, P.; Xia, F. Nat. Photonics 2013, 7, 394–399. (15) Thongrattanasiri, S.; Koppens, F. H.; de Abajo, F. J. G. Phys. Rev. Lett. 2012, 108, 047401. (16) Balendhran, S.; Walia, S.; Nili, H.; Sriram, S.; Bhaskaran, M. Small 2015, 11, 640–652. (17) Qiao, J.; Kong, X.; Hu, Z.-X.; Yang, F.; Ji, W. Nat. Commun. 2014, 5 . 13

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(18) Xia, F.; Wang, H.; Jia, Y. Nat. Commun. 2014, 5 . (19) Low, T.; Rodin, A. S.; Carvalho, A.; Jiang, Y.; Wang, H.; Xia, F.; Castro Neto, A. H. Phys. Rev. B 2014, 90, 075434. (20) Wood, J. D.; Wells, S. A.; Jariwala, D.; Chen, K.-S.; Cho, E.; Sangwan, V. K.; Liu, X.; Lauhon, L. J.; Marks, T. J.; Hersam, M. C. Nano Lett. 2014, 14, 6964–6970. (21) Rodin, A. S.; Carvalho, A.; Castro Neto, A. H. Phys. Rev. Lett. 2014, 112, 176801. (22) Fei, R.; Yang, L. Nano Lett. 2014, 14, 2884–2889. (23) Tran, V.; Soklaski, R.; Liang, Y.; Yang, L. Phys. Rev. B 2014, 89, 235319. (24) Li, L.; Yu, Y.; Ye, G. J.; Ge, Q.; Ou, X.; Wu, H.; Feng, D.; Chen, X. H.; Zhang, Y. Nat. Nanotechnol. 2014, 9, 372–377. (25) Wang, H.; Wang, X.; Xia, F.; Wang, L.; Jiang, H.; Xia, Q.; Chin, M. L.; Dubey, M.; Han, S.-j. Nano Lett. 2014, 14, 6424–6429. (26) Du, Y.; Liu, H.; Deng, Y.; Ye, P. D. ACS Nano 2014, 8, 10035–10042. (27) Buscema, M.; Groenendijk, D. J.; Blanter, S. I.; Steele, G. A.; van der Zant, H. S.; Castellanos-Gomez, A. Nano Lett. 2014, 14, 3347–3352. (28) Deng, Y.; Luo, Z.; Conrad, N. J.; Liu, H.; Gong, Y.; Najmaei, S.; Ajayan, P. M.; Lou, J.; Xu, X.; Ye, P. D. ACS Nano 2014, 8, 8292–8299. (29) Buscema, M.; Groenendijk, D. J.; Steele, G. A.; van der Zant, H. S.; CastellanosGomez, A. Nat. Commun. 2014, 5 . (30) Engel, M.; Steiner, M.; Avouris, P. Nano Lett. 2014, 14, 6414–6417. (31) Zhang, S.; Yang, J.; Xu, R.; Wang, F.; Li, W.; Ghufran, M.; Zhang, Y.-W.; Yu, Z.; Zhang, G.; Qin, Q. et al. ACS Nano 2014, 8, 9590–9596. 14

ACS Paragon Plus Environment

Page 14 of 16

Page 15 of 16

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

(32) Low, T.; Roldán, R.; Wang, H.; Xia, F.; Avouris, P.; Moreno, L. M.; Guinea, F. Phys. Rev. Lett. 2014, 113, 106802. (33) Jablan, M.; Buljan, H.; Soljačić, M. Phys. Rev. B 2009, 80, 245435. (34) Maier, S. A. Plasmonics: fundamentals and applications; Springer Science & Business Media, 2007. (35) Berger, L. I. Semiconductor Materials; CRC press, 1996. (36) Nikitin, A. Y.; Guinea, F.; García-Vidal, F.; Martín-Moreno, L. Phys. Rev. B 2011, 84, 161407. (37) Mikhailov, S. A.; Savostianova, N. Phys. Rev. B 2005, 71, 035320. (38) Nikitin, A. Y.; Low, T.; Martin-Moreno, L. Phys. Rev. B 2014, 90, 041407.

15

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Graphical TOC Entry

16

ACS Paragon Plus Environment

Page 16 of 16