Tunable plasmons in atomically thin gold

Full text

Turn on search term navigation

ARTICLE

Received 27 Nov 2013 | Accepted 4 Mar 2014 | Published 27 Mar 2014

The ability to modulate light at high speeds is of paramount importance for telecommunications, information processing and medical imaging technologies. This has stimulated intense efforts to master optoelectronic switching at visible and near-infrared frequencies, although coping with current computer speeds in integrated architectures still remains a major challenge. As a partial success, mid-infrared light modulation has been recently achieved through gating patterned graphene. Here we show that atomically thin noble metal nanoislands can extend optical modulation to the visible and near-infrared spectral range. We nd plasmons in thin metal nanodisks to produce similar absorption cross-sections as spherical particles of the same diameter. Using realistic levels of electrical doping, plasmons are shifted by about half their width, thus leading to a factor-of-two change in light absorption. These results, which we substantiate on microscopic quantum theory of the optical response, hold great potential for the development of electrical visible and near-infrared light modulation in integrable, nanoscale devices.

DOI: 10.1038/ncomms4548

Tunable plasmons in atomically thin gold nanodisks

A. Manjavacas1,w & F.J. Garca de Abajo2,3

1 IQFRCSIC, Serrano 119, 28006 Madrid, Spain. 2 ICFOInstitut de Ciencies Fotoniques, Mediterranean Technology Park, 08860 Castelldefels (Barcelona), Spain. 3 ICREAInstitucio Catalana de Recerca i Estudis Avancats, Passeig de Lluis Companys 23, 08010 Barcelona, Spain. w Present address: Department of

Physics and Astronomy and Laboratory for Nanophotonics, Rice University, Houston, Texas 77005, USA. Correspondence and requests for materials should be addressed to F.J.G.d.A. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 5:3548 | DOI: 10.1038/ncomms4548 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4548

Surface plasmons, the collective oscillations of conduction electrons in metallic structures, allow us to conne light down to deep subwavelength volumes1. In addition, they

couple strongly to electromagnetic elds2. Because of these properties, plasmons are excellent tools to engineer nanoscale devices for manipulating optical signals, without the limitation imposed by diffraction in far-eld set-ups. This has triggered a number of applications in areas as diverse as ultrasensitive biosensing3, improved photovoltaics4, plasmon-enhanced photodetection5 and photothermal cancer therapy6. The design of plasmonic structures with suitable spectral characteristics involves a careful choice of geometry and composition. In recent years, a vast amount of work has been devoted to producing nanostructures made of noble metals with controlled size and morphology, using in particular colloidal methods7 and lithography8.

Despite these advances in the control over the static characteristics of plasmons, the dynamical modulation of their frequencies and spatial proles remains ellusive, particularly in the visible and near-infrared (vis-NIR) parts of the spectrum. In this context, slow mild changes of the plasmon frequency have been produced by electrochemically injecting electrons in metal nanoparticles9, by electrically driving liquid crystals containing plasmonic particles10, and through controllable metamaterial designs11,12. Magneto-optical modulation has also been explored to control plasmons in noble metal structures13. Hybrids of plasmonic and conductive oxides have been proposed14,15, as well as colloids based on different materials16. However, we still need to devise new methods to produce larger and faster control over plasmons, as required for nanoscale optical commutation and light modulation at high speeds.

Recently, the emergence of graphene17 as a novel plasmonic material1820 has opened up new paths towards the design of dynamically tunable plasmonic devices. Electrically doped graphene supports surface plasmons whose frequency can be efciently varied by changing the level of doping2123. Consequently, the resulting modulation is intrinsically fast because it can be driven by charge-carrier injection using conventional electric gating technology. This promising material has been so far shown to support mid-infrared and lower-frequency plasmons2125, while vis-NIR modes are being pursued by reducing the size of the structures26,27 and increasing the level

of doping23. The search for plasmon modulation in the vis-NIR is thus still ongoing, as these are spectral regions of utmost importance for sensing and optical signal processing technologies.

The origins of the excellent tunability of plasmons in graphene can be found in both the atomic thickness and the peculiar electronic structure of this material. The latter is characterized by a linear dispersion relation, which leads to a vanishing of the density of states at the Fermi level, so that a relatively small density of injected charge carriers produces substantial optical gaps in which collective plasmon modes emerge17,19. Although this unique feature cannot be easily transported to conventional plasmonic materials, such as gold, we can still mimic graphene plasmonics by going to atomically thin noble metals, whose optical response should be more susceptible to doping than traditional thicker layers. In particular, monolayer gold, the synthesis of which has been mastered for a long time in the context of surface science28, presents the advantage of having a plasma frequency compatible with the existence of plasmons in the vis-NIR29.

Here we show that single-monolayer gold disks (SMGDs) with diameters of the order of 10 nm support surface plasmons with large cross-sections comparable to their geometrical areas. The frequencies of these excitations lie in the vis-NIR and can be efciently modulated using attainable concentrations of doping charge carriers, which can be provided via electrical doping using, for example, backgating technology. We also analyse the optical response of periodic arrays of SMGDs, for which we predict an absorbance B25% for metal layer lling fractions B40%.

ResultsElectrically tunable optical response. The system under study is depicted in Fig. 1a. We consider a gold nanodisk of diameter D, extracted from a single (111) atomic layer of gold. We take the thickness of the gold monolayer to be equal to the separation between (111) atomic planes in bulk gold (that is, a0=

p , where

a0 0.408 nm is the atomic lattice constant). Incidentally, our

results are rather independent of the choice of disk thickness when this is small compared with the diameter, as long as the total valence charge is preserved (see Supplementary Fig. 1). As a rst step in our analysis, we simulate the optical response of a SMGD classically by modelling it as a thin disk described by

a b c

Doping charge density (1013cm2)

Sphere

3.5

6.9

1.7

13.8

3.5

Doping charge density (1013cm2)

Disk

13.8

6.9

Disk diameter, D

4 0 2 4

60 30 0 30 60

Potential (V)

1.0

0.5

0.0

0.80.40.0 0.4 0.8

Potential (V)

Extinction, [afii9846]/area

p/FWHM

D = 20 nm

Gold

111 plane

D = 20 nm

100 50 0 50 100

0.9 1.2 2.5 3.0 16 0 8 16

Energy (eV)

Figure 1 | Optical response and electrical tunability of SMGDs. (a) We consider a SMGD of diameter D carved from a single (111) atomic layer. The disk thickness is a0=

p , where a0 0.408 nm is the atomic lattice constant. (b) Extinction cross-section of a D 20 nm SMGD for different doping

charge-carrier densities (see upper legend). A doping density of 13.8 1013 cm 2 corresponds to a total of 443 additional charge carriers in the disk.

For comparison, we also plot the cross-section for a gold nanosphere of the same diameter and total doping charge, clearly showing an almost negligible degree of tunability. The particles are assumed to be homogenously doped and described classically through the local dielectric function tabulated from measured optical data (solid curves)29. Results obtained from a Drude dielectric function (equation (1)) are shown for comparison (broken curves).(c) Plasmon frequency shift relative to the width of the plasmon resonance for the disk (orange) and the sphere (green) of panel b as a function of doping density. The small scales indicate the potential at the disk/sphere surface for different disk doping densities.

2 NATURE COMMUNICATIONS | 5:3548 | DOI: 10.1038/ncomms4548 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4548 ARTICLE

a frequency-dependent homogenous dielectric function E(o). More precisely, we calculate the extinction cross-section s by solving Maxwells equations using the boundary-element method30. Interestingly, for a diameter D 20 nm, the cross-

section is dominated by a NIR plasmon at an energy B1 eV and it exceeds the geometrical area of the disk (see left part of Fig. 1b).

It is convenient to separate the contribution from s-band electrons in the dielectric function as a Drude term,

Eo Eb

o2p

o o ig

; 1

where Eb accounts for the effect of background screening due to d-band electrons, :op 9.06 eV is the classical plasmon energy

associated with s-valence electrons (see Methods), and :g 71

meV is an inelastic width (we adopt this value of the damping throughout this work). As explained below, we introduce doping in the classical model by changing op in equation (1).

In general, the description of the response of gold in the vis-NIR region including interband transitions requires to use either experimental data29 or a sophisticated multi-Lorentzian model31 for E(o), which yields a o-dependent background Eb Eo o2p= o o ig

. However, we show in Fig. 1b that a

simple Drude model for equation (1) (dashed curves), consisting in xing Eb 9.5 for all frequencies, produces a satisfactory level

of accuracy at the observed relatively low disk-plasmon energies compared with the results obtained from tabulated optical data29 (solid curves). In addition, as the Drude model (that is, constant Eb) provides a natural connection with the quantum mechanical approach described below, we use it for disks in what follows.

We consider next the effect of electrical doping. The spatial distribution of additional charge carriers depends on the doping conguration, as it can be, for example, homogenous for disks connected to a non-absorbing gate (for example, indium tin oxide) or inhomogenous in self-standing charged disks, although the plasmon energies and spatial proles are expected to be similar in both cases based on our experience with graphene disk plasmons32. For simplicity, we assume homogenously doped disks in what follows. The additional doping charge density n adds up to the undoped s-band density n0 meo2p=4pe2, which yields an areal electron density

n0a0=

p 1:4 1015 cm 2 that is rather close to the s-band

areal electron density in neutral monolayer gold, 4=

p a20. The

doping charge is thus introduced by changing the bulk plasma frequency to op [(4pe2/me)(n0 n)]1/2 in equation (1). Now,

the addition of a moderate amount of doping electrons (B510% of n0) results in signicant blue shifts and increase in the strength of the plasmon resonance (cf. purple and green curves of Fig. 1b). Obviously, the injection of similar amounts of holes produces the opposite effects (Fig. 1b, orange curve). The small thickness of the disk is a key factor in producing such dramatic modications in the optical response using realistic doping densities. In fact, repeating this operation with a gold nanosphere of the same diameter, we also observe a prominent plasmon (Fig. 1b, B2.5 eV region), but it remains unchanged when adding similar amounts of doping charges. In the sphere, the doping charges pile up in the outermost atomic layer33, but this produces the same extinction cross-section as if the charges where homogenously distributed over its entire volume and, therefore, the change in bulk charge density is substantially reduced with respect to the disk.

Figure 1c compares the modulation of the nanodisk and the sphere. In particular, we plot the frequency shift normalized to the full width at half maximum (FWHM) for the plasmon resonance as a function of doping charge density. In contrast to the negligible tunability of the sphere, the disk allows shifts comparable to the FWHM to be electrically induced. Incidentally,

the doping densities here considered produce realistic values of the electrostatic potential at the surface of these nanoparticles (Fig. 1c, small scales), indicating that they are compatible with currently available gating technology23.

Quantum mechanical effects. For nanoparticles of only a few nanometres in size, the above classical description fails to account for spatial dispersion and quantum connement effects34,35, which generally require models based on a quantum mechanical treatment of valence electrons and their interactions. Here we use the random-phase approximation (RPA) to calculate the optical response of SMGDs (see Methods) within the electrostatic approximation, which should be rather accurate, given the small sizes of the particles under consideration. This allows us to determine the validity of the classical approach and explore the response of nanodisks with smaller diameters.

We use particle-in-a-box states to describe independent s-band electrons. The cylindrical box has the same dimensions as in the classical calculations (see above) and it is surrounded by an innite potential. The RPA susceptibility is then evaluated using these electron states to obtain the induced charge density, which in turn allows us to compute the optical extinction of the disk. In addition, we model screening due to d-band electrons through an array of point dipoles placed at the atomic positions in the (111) layer and with polarizability adjusted to render an effective background permittivity Eb in the bulk material (see Methods).

Apart from the relative position of these dipoles with respect to the disk centre, our quantum description only depends on the three same parameters as the classical Drude theory (that is, Eb,

op and g).

A major assumption we are making is that g takes the same value as in the bulk metal. We use this as a reasonable estimate because s-band electrons are rather delocalized along directions parallel to the layer (that is, similar to the bulk), while they are narrowly conned to the ground state across the transversal direction, so that plasmons result from in-plane motion. However, the actual value of g might depend on the detailed coupling of valence electrons to impurities and to the atomic lattice. Concerning d-band screening, our effective-dipoles approach should provide a more realistic description than a homogenous polarizable background. Although the discreteness of the dipole lattice can have strong effects in small islands, we nd converged results for large islands, which are independent of the alignment of the dipole lattice relative to the disk centre.

Figure 2 shows the extinction cross-section normalized to the disk area for SMGDs of different diameters ranging from 3 to 15 nm. The main conclusions from this gure are as follows:(1) the extinction cross-sections are of the order of the disk area;(2) the plasmon energy increases with decreasing diameter D, exhibiting an approximate / 1=

p dependence, similar to what one nds in graphene nanodisks23; (3) quantum calculations produce energies above those predicted by classical theory, as well as broader plasmon peaks, but the discrepancy between the two models decreases with increasing diameter; (4) in the absense of d-band screening (Fig. 2a, obtained by setting Eb 1), both levels

of description give rise to smooth plasmon peaks, in contrast to the quantum results obtained when d-band screening is switched on (Fig. 2b, with Eb 9.5); (5) d-band screening also leads to a

redshift of the plasmons, which is more pronounced at small sizes. Incidentally, the induced charge associated with the plasmon exhibits a dipolar prole dressed with radial oscillations mimicking those of Friedel oscillations, which are particularly intense for small diameters (see Supplementary Fig. 2).

Similar blue shifts with respect to classical local theory are also found in small noble metal particles36, the origin of which is a

NATURE COMMUNICATIONS | 5:3548 | DOI: 10.1038/ncomms4548 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4548

Including d-band screening

Without d-band screening

Quantum model Classical model

15 nm

4 nm

15 nm

4 nm

Extinction normalized to disk area, [afii9846]/area

14 nm

3 nm

13 nm

12 nm

11 nm

10 nm

9 nm

8 nm

9 nm

8 nm

7 nm

6 nm

5 nm

3 nm

1 2 3

Energy (eV)

Figure 2 | Optical response of individual SMGDs. We plot the extinction cross-section normalized to the geometrical area for different diameters D, calculated using the quantum model (solid curves) and a classical description (dashed curves). Results obtained with and without inclusion of d-band screening are shown in b and a, respectively.

combination of non-local and quantum effects, particularly owing to the surface spill out of s electrons beyond the polarizable background of d-band electrons. In simple metals such as aluminium, the spill out produces smaller electron densities near the surface and, consequently, also smaller surface plasmon frequencies. In contrast, in noble metals, the spill out results in a weaker interaction with the localized d electrons and, thus, it leads to an increase in the observed frequency, which overcomes the redshift due to the smaller electron density37. We incorporate here the nite extension of s electrons across the normal direction, combined with the localization of the effective d-band dipoles, leading to similar blue shifts. Interestingly, our quantum model predicts splitting of the plasmon into multiple peaks for small disks when d-band screening is included (see for example, the D 3 nm spectrum in Fig. 2b). The presence of these peaks,

which are rapidly coalescing into a single plasmon resonance for D48 nm, is a manifestation of the discrete character of the interaction with d-band electrons, which is more pronounced for small Ds. The jumps observed in the FWHM also shares a similar origin. It should be noted that these effects could be sensitive to the exact form of the s-electron transversal wave function in the smallest islands under consideration, which require a more atomistic analysis, based, for example, on density functional theory3840. Similarly, the spectra for the smallest disks depend on the alignment of the d-band dipole lattice relative to the edges. In practice, one-dimensional faceting of the edges becomes an important source of anisotropy, which can contribute to broaden the spectra for Do5 nm.

The convergence of the quantum model to the classical description for increasing diameter is clearly observed in Fig. 3, which summarizes the plasmon energies and widths observed in the spectra of Fig. 2. Here we dene the FWHM as the frequency interval around the peak maximum that contains half of its area; this denition coincides with the standard FWHM for individual Lorentzian resonances, but it can be applied to multiple resonances as well to yield an overall width (in particular to the lower quantum-model spectra of Fig. 2b). Within the electrostatic

limit under consideration, the FWHM predicted by the classical model is independent of diameter and equals the damping energy :g 71 meV (see equation (1)). In contrast, the quantum model

leads to a signicant increase in the FWHM for small diameters, essentially as a consequence of Landau damping, which involves inelastic decay of plasmons to electronhole pairs for momentum transfers Bo/uF, where uF is the Fermi velocity (B106 m s 1, see

Supplementary Fig. 3a). As the momentum transfer provided by the breaking of translational invariance in a disk is B2p/D, the onset of Landau damping is expected to occur at DB2puF/

oB4 nm, in qualitative agreement with the results shown in Fig. 3c,d. An intuitive estimate can be stablished from the electron mean free path uF/gB10 nm, which determines the rate of collisions with the edges (that is, events that provide the noted momentum), and is also in agreement with the trends observed in Fig. 3c,d, although the value of g regarded as a parameter simply produces an additional contribution to the FWHM that is independent of D, and the ultimate origin of broadening for small sizes can be found in Landau damping.

As discussed in Fig. 1, the optical response of SMGDs can be modied through the addition of small amounts of charge carriers to relatively large disks (D 20 nm), for which classical

theory is rather accurate (see Figs 2 and 3). Using smaller SMGDs, we obtain qualitatively similar results as for larger disks (see Fig. 4 for an analysis of a D 8 nm doped disk). Given the

small disk size, we compare classical (Fig. 4df) and quantum (Fig. 4ac) results, showing again a blue shift and plasmon broadening in the latter relative to the former.

In contrast to the nearly linear plasmon shift with doping charge density predicted by classical theory (Fig. 4e), the quantum model leads to initially smaller modulation at low doping (Fig. 4b), which increases to a faster pace than the classical results for higher doping. This nonlinear dependence of the plasmon energy on the doping density could be exploited for improved light modulation by operating around a highly doped SMGD conguration. In particular, the plasmon shift can be as large as the FWHM when the density of s-band electrons is changed by (510)%.

Interestingly, the nonlinear plasmon shifts observed in the quantum model become oscillatory when examining the maximum extinction cross-section and the FWHM (Fig. 4c). The oscillations of these two quantities are out of phase, as required to satisfy the f-sum rule41, and can be traced back to the discreteness of the electronic energies. Importantly, in all cases the maximum cross-section is of the order of the disk area (Fig. 4c,f), thus providing good coupling to light for potential applications to modulation devices.

Periodic arrangement of single-monolayer gold nanodisks. The large optical strength and degree of electrical tunability discussed above for SMGDs can be exploited to modulate light that is either transmitted or reected by a periodic array of such structures. We consider a hexagonal array of D 8 nm disks in Fig. 5 with dif

ferent values of the array spacing d. Given the large mismatch between D and the resonant light wavelength (B830 nm), we approximate the disks as point dipoles of polarizability extracted as explained in Methods (see also Supplementary Fig. 4, where we show that higher-order multipoles play only a small role for the relative distances under consideration). Following previous analytical methods42 to compute the absorbance A, we nd remarkably large values (for example, A 25% for d 1.5D, see

Fig. 5b), given the small amount of gold in the structure (submonoatomic layer lm). The fractional change in absorbance driven by electrical doping (Fig. 5c) is rather independent of lattice spacing and reaches B70% for a 10% variation in the

4 NATURE COMMUNICATIONS | 5:3548 | DOI: 10.1038/ncomms4548 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4548 ARTICLE

Without d-band screening

Including d-band screening

3.0

Plasmon energy (eV)

2.5

1.0

Plasmon energy (eV) FWHM (meV)

2.5

2.0

1.5

Quantum model

Classical model

Disk diameter, D

1.0

240

FWHM (meV)

160

71 meV

3 6 12 15

Disk diameter (nm) Disk diameter (nm)

Figure 3 | Comparison of quantum and classical plasmon energies and widths. Energy (a,b) and FWHM (c,d) of the plasmonic resonance of individual SMGDs as a function of disk diameter, calculated from quantum (green circles) and classical (orange triangles) models. Results obtained with and without inclusion of d-band screening are shown in b and d, and a and c, respectively. The dashed curves in c and d indicate the intrinsic broadening :g 71 meV introduced in the Drude formula (equation (1)) and in the RPA susceptibility (equation (4)).

1.36

0.8

Extinction, [afii9846]/area

[afii9846]/area

Doping charge density (1013cm2)

14 7 07 14

Doping charge density (1013cm2)

Quantum model Classical model

D = 8 nm

0.8

Extinction, [afii9846]/area

0.4

14 7 07 14

0.4

0.0

1.0 1.5 2.0

Energy (eV)

1.0 1.5 2.0

Energy (eV)

1.68

1.44

1.49

1.28

Plasmon energy (eV)

1.60

FWHM

1.42

FWHM

1.52

1.35

Plasmon energy

140

0.7

0.6

0.4 80

120 100

FWHM (meV)

0.960.88

[afii9846]/area

FWHM (meV)

0.5

0.80

14 0 14

7 7

14 0 14

7 7

Doping charge density (1013cm2) Doping charge density (1013cm2)

Figure 4 | Quantum versus classical analysis of the electrically tunable optical response. (a,d) Extinction cross-section normalized to the geometrical area for a D 8 nm SMGD calculated with different doping charge densities. (b,e) Plasmon energy as a function of doping charge density.

We indicate the FWHM of the plasmon resonances as shadowed regions. (c,f) Optical extinction cross-section at the plasmon peak energy (green curves and symbols, left scale) and FWHM (orange, right scale) as a function of doping density. Quantum mechanical calculations (ac) are compared with classical results (df), including d-band screening in all cases.

s-band electron density. The potential of patterned gold monolayers for electro-optical modulation in the NIR is thus excellent.

DiscussionIn summary, we have simulated both classically and quantum mechanically the plasmonic response and performance in

NATURE COMMUNICATIONS | 5:3548 | DOI: 10.1038/ncomms4548 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4548

0.26

0.75

d = 1.5 D

d = 2 D d = 3 D

d = 1.5 D d = 2 D

Absorbance,

d = 3 D

0.13

A/A max

0.50

0.25

0.00

1.0 1.5 2.0

Energy (eV)

14 0 14

7 7

Doping charge density (1013cm2)

Figure 5 | Electrical modulation of the absorbance of an hexagonal periodic arrangement of SMGDs. (a) Scheme of the system under study. (b) Absorbance spectrum of undoped nanodisks (diameter D 8 nm) for different values of the array spacing d. (c) Modulation of the absorbance

relative to the undoped state as a function of doping charge density for different array spacings.

electro-optical modulation of gold nanodisks carved from a single (111) atomic layer. Our RPA calculations incorporate the wave functions of free s-valence electrons evolving in a circular box, as well as an adjusted distribution of dipoles to account for d-band screening. Despite the atomic thickness of the disks, this quantum mechanical description converges smoothly to the results of classical dielectric theory, based on the bulk, frequency-dependent dielectric function of gold. This is a remarkable result, which can be intuitively understood from the fact that the electron current associated with the plasmons ows along the gold layer and, thus, it is rather insensitive to electron connement within the small lm thickness. Nontheless, nonlocality plays a crucial role, leading to strong plasmon blue shifts, as well as splitting due to the complex interaction with the d band. We estimate that non-local effects become dominant when the disk diameter is below B10 nm.

Remarkably, the disks interact strongly with light, giving rise to extinction cross-sections exceeding their geometrical areas in the vis-NIR. We have also shown that the optical response of SMGDs can be efciently modulated through the addition or removal of realistic concentrations of doping charge carriers using, for example, gating technology. In particular, periodic patterns of monolayer gold appear to be a suitable solution for combining strong plasmonic response and high doping, for example, using an electrical backgate, because the average charge density of the layer is simply determined by capacitor theory for a xed distance from the gate and, thus, the actual doping charge density in the metal scales with the inverse of the areal lling fraction occupied by the gold. Similar results are expected for lms consisting of only a few atomic layers, although the degree of modulation is then reduced because the doping charge has to be shared across the increased thickness. Other plasmonic metals such as silver and copper should nd similar degree of tunability (see Supplementary Fig. 5). In particular, the small plasmon width of silver compared with gold makes it an attractive candidate to drive plasmon shifts well beyond the FWHM. In addition, the lower d-band screening in this material should result in higher plasmon energies, reaching the visible in small disks or, equivalently, the NIR for larger disk diameters.

It should be stressed that, while the synthesis of single-layer gold is a mature eld28, the fabrication of laterally conned thin gold nanostructures represents a technical challenge, which could benet from advances in lithography and self-assembly. Alternatively, one could use a continuous gold layer, which also exhibits large electrical tunability of its propagating plasmons (see Supplementary Fig. 6), coupled to external light by decorating it with dielectric colloids to bridge the light-plasmon momentum mismatch (that is, this is essentially what nanostructuration does in the SMGDs that we study above). The resulting planar structures hold great potential for light modulation at vis-NIR frequencies, which could be the basis of a new generation of

electrically tunable optical devices with applications ranging from sensing to nanoscale spectral photometry.

Methods

Quantum mechanical RPA simulations. We consider small disk sizes compared with the light wavelength, so that we work in the electrostatic limit. Within this approximation, assuming an overall monochromatic time dependence e iot with frequency o, the induced charge density rind can be expressed in terms of the self-consistent potential f as

rindr; o Z d3r0w0r; r0; ofr0; o w0 f; 2

where w0 is the non-interacting susceptibility associated with the s-valence electrons of gold, and the last identity denes a matrix notation in which matrix multiplication involves integration over space coordinates. We obtain w0 within the RPA41 in which a one-electron picture is assumed and only individual electronhole pair excitations are explicitly considered. We further approximate the wave functions of valence electrons by the solutions of a cylindrical box with the same diameter as the nanodisk and a thickness corresponding to the separation between (111) atomic planes in bulk gold (that is, a0=

p 0:236 nm, see Fig. 1a).

More precisely,

clmr NlmJm QlmR

eimjg1z; 3

where Nlm is a normalization constant, Qlm 2zlm/D, zlm is the lth zero of the

Bessel function Jm and g1z sin p 3

p z=a0

yields the dependence on the coordinate z normal to the disk. For simplicity, we are assuming that the z dependence is separable in the wave function, so that electron diffraction effects at the disk edges are not important. Moreover, we assume that the electrons remain in the ground state of the vertical cavity of thickness a0=

p , which is a reasonable approximation if we consider that the rst excited state (that is, g2z sin 2p 3

p z=a0 )

lies B20 eV above g1, well beyond the Fermi and vacuum levels.

With the wave functions of equation (3), we can write the susceptibility as

w0 r; r0; o

2e2

l;l ;m;m

fl m flm

clmrc lmr0c lm rcl m r0

o elm el

ig=2 ; 4

where spin degeneracy is simply included through an overall factor of 2, elm 2Q2lm=2me is the energy of state clm (notice that the energy associated

with z motion cancels out in equation (4), so we disregard it), g is an intrinsic relaxation time, which we take from a Drude t (equation (1)) to measured optical data29 (:g 71 meV), and flm {exp[(elm EF)/kBT] 1} 1 is the FermiDirac

distribution function, evaluated here at T 0. The method used to ll the energy

levels in the disk has been discussed in Supplementary Fig. 3.

The total potential f is the sum of the external potential fext and the potential produced by the induced charges

f fext u rind; 5

where u(rr0) 1/|rr0| is the Coulomb interaction. From here and equation (2),

we solve the induced charge density as

rind w0 1 u w0

1 fext:

As the polarization along the direction normal to the disk is expected tobe negligible, we focus on parallel components and write fext Reij (that is,

we focus on solutions with m 1 azimuthal symmetry), where j is the azimuthal

angle of (x,y) and R

x2 y

2 . This allows us to obtain the in-plane polarizability by calculating

Z d3r Re ij rind r; o

6 NATURE COMMUNICATIONS | 5:3548 | DOI: 10.1038/ncomms4548 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4548 ARTICLE

Finally, the extinction cross-section is obtained from

so 4po=c

Im ao

f g:

Inclusion of d-band screening. Deeper electrons in the d band are relatively localized in the gold atoms and, therefore, we model them by assuming a background of polarizable point particles at the atomic positions in the (111) layer (see lower part of Fig. 1a). The polarizability ab of these particles is adjusted to t the experimentally measured bulk dielectric function of gold Eexp. That is, if we

subtract the Drude s-band contribution from Eexp (see equation (1)), we obtain

the background permittivity Eb Eexp o2p=oo ig, where o2p 4pe2n0=me is

determined by the s-band electron density n0 4=a30 5:9 1028m 3. This

yields :opE9.01 eV, which is slightly different from the best t of equation (1) to measured data29 (9.06 eV), from which we also nd Eb 9.5. Now, the

ClausiusMossotti relation43 leads to

4pn0

2 :

Using dyadic notation, the susceptibility tensor of the background dipoles

reduces to w0br; r0 P

r abdr rjdr0 rj r

!0, where j runs over the

positions of the metal atoms, whereas r (r

!0) acts on r-dependent (r0-dependent) functions to the left (right) of operator w0b. As the charge induced through both s and d bands contribute together to the full potential, we can rewrite equation (2) as

rind w0 w0b

1 u w0 w0b

1 fext;from which we calculate the disk polarizability and the extinction cross-section as explained above.

to take into account the effect of d-band screening. Using this expression together with equation (5), the total induced charge density becomes

rind w0 w0b

References

1. Novotny, L. & Hecht, B. Principles of Nano-Optics (Cambridge Univ. Press,

2006).

2. Halas, N. J., Lal, S., Chang, W., Link, S. & Nordlander, P. Plasmons in strongly coupled metallic nanostructures. Chem. Rev. 111, 39133961 (2011).3.lvarez-Puebla, R. A., Liz-Marzn, L. M. & Garca de Abajo, F. J. Light concentration at the nanometer scale. J. Phys. Chem. Lett. 1, 24282434 (2010).

4. Atwater, H. A. & Polman, A. Plasmonics for improved photovoltaic devices. Nat. Mater. 9, 205213 (2010).

5. Knight, M. W., Sobhani, H., Nordlander, P. & Halas, N. J. Photodetection with active optical antennas. Science 332, 702704 (2011).

6. ONeal, D. P., Hirsch, L. R., Halas, N. J., Payne, J. D. & West, J. L. Photo-thermal tumor ablation in mice using near infrared-absorbing nanoparticles. Cancer Lett. 209, 171176 (2004).

7. Grzelczak, M., Prez-Juste, J., Mulvaney, P. & Liz-Marzn, L. M. Shape control in gold nanoparticle synthesis. Chem. Soc. Rev. 37, 17831791 (2008).

8. Stokes, N., McDonagh, A. M. & Cortie, M. B. Preparation of nanoscale gold structures by nanolithography. Gold Bull. 40/4, 310320 (2007).

9. Mulvaney, P., Prez-Juste, J., Giersig, M., Liz-Marzn, L. M. & Pecharromn, C. Drastic surface plasmon mode shifts in gold nanorods due to electron charging. Plasmonics 1, 6166 (2006).

10. Chu, K. C., Chao, C. Y., Chen, Y. F., Wu, Y. C. & Chen, C. C. Electrically controlled surface plasmon resonance frequency of gold nanorods. Appl. Phys. Lett. 89, 103107 (2006).

11. Boardman, A. D. et al. Active and tunable metamaterials. Laser Photon. Rev. 10, 00012 (2010).

12. Liu, A. Q., Zhu, W. M., Tsai, D. P. & Zheludev, N. I. Micromachined tunable metamaterials: a review. J. Opt. 14, 114009 (2012).

13. Armelles, G., Cebollada, A., Garca-Martn, A. & Gonzlez, M. U. Magnetoplasmonics: combining magnetic and plasmonic functionalities. Adv. Opt. Mater. 1, 1035 (2013).

14. Feigenbaum, E., Diest, K. & Atwater, H. A. Unity-order index change in transparent conducting oxides at visible frequencies. Nano Lett. 10, 21112116 (2010).

15. Abb, M., Albella, P., Aizpurua, J. & Muskens, O. L. All-optical control of a single plasmonic nanoantenna-ITO hybrid. Nano Lett. 11, 24572463 (2011).

16. Comin, A. & Manna, L. New materials for tunable plasmonic colloidal nanocrystals. Chem. Soc. Rev. doi:http://dx.doi.org/10.1039/C3CS60265F

Web End =10.1039/C3CS60265F (2014).

17. Castro Neto, A. H., Guinea, F., Peres, N. M. R., Novoselov, K. S. & Geim, A. K. The electronic properties of graphene. Rev. Mod. Phys. 81, 109162 (2009).

18. Jablan, M., Buljan, H. & Soljai, M. Plasmonics in graphene at infrared frequencies. Phys. Rev. B 80, 245435 (2009).

19. Koppens, F. H. L., Chang, D. E. & Garca de Abajo, F. J. Graphene plasmonics: a platform for strong light-matter interactions. Nano Lett. 11, 33703377 (2011).

20. Grigorenko, A. N., Polini, M. & Novoselov, K. S. Graphene plasmonics. Nat. Photon. 6, 749758 (2012).

21. Chen, J. et al. Optical nano-imaging of gate-tunable graphene plasmons. Nature 487, 7781 (2012).

22. Fei, Z. et al. Gate-tuning of graphene plasmons revealed by infrared nano-imaging. Nature 487, 8285 (2012).

23. Fang, Z. et al. Gated tunability and hybridization of localized plasmons in nanostructured graphene. ACS Nano 7, 23882395 (2013).

24. Ju, L. et al. Graphene plasmonics for tunable terahertz metamaterials. Nat. Nanotechnol. 6, 630634 (2011).

25. Brar, V. W., Jang, M. S., Sherrott, M., Lopez, J. J. & Atwater, H. A. Highly conned tunable mid-infrared plasmonics in graphene nanoresonators. Nano Lett. 13, 25412547 (2013).

26. Manjavacas, A., Thongrattanasiri, S. & Garca de Abajo, F. J. Plasmons driven by single electrons in graphene nanoislands. Nanophotonics 2, 139151 (2013).

27. Manjavacas, A. et al. Tunable molecular plasmons in polycyclic aromatic hydrocarbons. ACS Nano 7, 36353643 (2013).

28. Green, A. & Bauer, E. Gold monolayers on silicon single crystal surfaces. Surf. Sci. 103, L127L133 (1981).

29. Johnson, P. B. & Christy, R. W. Optical constants of the noble metals. Phys. Rev. B 6, 43704379 (1972).

30. Garca de Abajo, F. J. & Howie, A. Retarded eld calculation of electron energy loss in inhomogeneous dielectrics. Phys. Rev. B 65, 115418 (2002).

31. Raki, A. D., Djurii, A. B., Elazar, J. M. & Majewski, M. L. Optical properties of metallic lms for vertical-cavity optoelectronic devices. Appl. Opt. 37, 52715283 (1998).

32. Thongrattanasiri, S., Silveiro, I. & Garca de Abajo, F. J. Plasmons in electrostatically doped graphene. Appl. Phys. Lett. 100, 201105 (2012).

33. Lang, N. D. & Kohn, W. Theory of metal surfaces: charge density and surface energy. Phys. Rev. B 1, 45554568 (1970).

34. Zuloaga, J., Prodan, E. & Nordlander, P. Quantum description of the plasmon resonances of a nanoparticle dimer. Nano Lett. 9, 887891 (2009).

35. Savage, K. J. et al. Revealing the quantum regime in tunnelling plasmonics. Nature 491, 574577 (2012).

36. Kreibig, U. & Vollmer, M. Optical Properties of Metal Clusters (Springer-Verlag, 1995).

37. Liebsch, A. Surface-plasmon dispersion and size dependence of Mie resonance: silver versus simple metals. Phys. Rev. B 48, 1131711328 (1993).

38. Onida, G., Reining, L. & Rubio, A. Electronic excitations: density-functional versus many-body Greens-function approaches. Rev. Mod. Phys. 74, 601659 (2002).

39. Sakko, A., Rossi, T. P., Enkovaara, J. & Nieminen, R. M. Atomistic approach for simulating plasmons in nanostructures. Appl. Phys. A doi:http://dx.doi.org/10.1007/s00339-013-8034-3

Web End =10.1007/s00339-013- http://dx.doi.org/10.1007/s00339-013-8034-3

Web End =8034-3 (2013).

40. Lisinetskaya, P. G. & Mitri, R. Ab initio simulations of light propagation in silver cluster nanostructures. Phys. Rev. B 89, 035433 (2014).

41. Pines, D. & Nozires, P. The Theory of Quantum Liquids (W. A. Benjamin Inc., 1966).

42. Thongrattanasiri, S., Koppens, F. H. L. & Garca de Abajo, F. J. Complete optical absorption in periodically patterned graphene. Phys. Rev. Lett. 108, 047401 (2012).

43. Ashcroft, N. W. & Mermin, N. D. Solid State Physics (Harcourt College Publishers, 1976).

Acknowledgements

This work has been supported in part by the European Commission (Graphene Flagship

CNECT-ICT-604391 and FP7-ICT-2013-613024-GRASP). A.M. acknowledges nancial

support from the Spanish MEC through the FPU programme and from the Evans

Attwell-Welch Postdoctoral Fellowship for Nanoscale Research, administered by the

Richard E. Smalley Institute for Nanoscale Science and Technology.

Author contributions

F.J.G.d.A. proposed the study. A.M. and F.J.G.d.A. worked out the theory, discussed the

results and wrote the paper. A.M. performed the numerical calculations.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/

http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/

http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Manjavacas, A. & Garca de Abajo, F.J. Tunable plasmons in

atomically thin gold nanodisks. Nat. Commun. 5:3548 doi: 10.1038/ncomms4548 (2014).

NATURE COMMUNICATIONS | 5:3548 | DOI: 10.1038/ncomms4548 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

Word count: 6702

Show less

Abstract

Translate

The ability to modulate light at high speeds is of paramount importance for telecommunications, information processing and medical imaging technologies. This has stimulated intense efforts to master optoelectronic switching at visible and near-infrared frequencies, although coping with current computer speeds in integrated architectures still remains a major challenge. As a partial success, mid-infrared light modulation has been recently achieved through gating patterned graphene. Here we show that atomically thin noble metal nanoislands can extend optical modulation to the visible and near-infrared spectral range. We find plasmons in thin metal nanodisks to produce similar absorption cross-sections as spherical particles of the same diameter. Using realistic levels of electrical doping, plasmons are shifted by about half their width, thus leading to a factor-of-two change in light absorption. These results, which we substantiate on microscopic quantum theory of the optical response, hold great potential for the development of electrical visible and near-infrared light modulation in integrable, nanoscale devices.

Details

Title

Tunable plasmons in atomically thin gold nanodisks

Author

Manjavacas, A; García De Abajo, Fj

Pages

3548

Publication year

2014

Publication date

Mar 2014

Publisher

Nature Publishing Group

e-ISSN

20411723

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/ncomms4548

ProQuest document ID

1510578268

Tunable plasmons in atomically thin gold nanodisks

Jump to:

Full text

Abstract

Details

Suggested sources