Tunable Spectral Ordering of Magnetic Plasmon Resonances in Noble

May 18, 2018 - However, based on analytic modeling and numerical simulation, we show that angle-resolved cathodoluminescence offers an approach to ...
0 downloads 0 Views 3MB Size
Subscriber access provided by READING UNIV

Article

Tunable Spectral Ordering of Magnetic Plasmon Resonances in Noble Metal Nanoclusters Nicholas P. Montoni, Steven C. Quillin, Charles Cherqui, and David J. Masiello ACS Photonics, Just Accepted Manuscript • DOI: 10.1021/acsphotonics.8b00519 • Publication Date (Web): 18 May 2018 Downloaded from http://pubs.acs.org on May 19, 2018

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 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 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.

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 28 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

ACS Photonics

Tunable Spectral Ordering of Magnetic Plasmon Resonances in Noble Metal Nanoclusters Nicholas P. Montoni, Steven C. Quillin, Charles Cherqui, and David J. Masiello

Department of Chemistry, University of Washington, Seattle, WA 98195

E-mail: [email protected]

Abstract Experimental characterization of the optical and magnetic properties of noble metal nanoclusters has exposed a size regime in between single-ring and extended two-dimensional nanocluster networks where the behavior of magnetic plasmon resonances is not well understood. In this intermediate size regime, individual electric dipole plasmons on each nanoparticle within the cluster can hybridize into delocalized magnetic modes that couple to either the electric or magnetic eld of light with seemingly arbitrary energy ordering. Here, using a coupled-dipole model that includes fully-retarded interactions between electric dipole plasmons, we show that magnetic plasmon resonance energies can be controllably tuned and even made to cross as a function of nanocluster size. Experimental conrmation of this prediction is challenging because optical selection rules dictate the simultaneous excitation of many spectrally-overlapping magnetic plasmon modes. However, based on analytic modeling and numerical simulation, we show that angle-resolved cathodoluminescence oers an approach to cleanly isolate these magnetic modes and observe the size-dependence of their spectrum. Not only does our work clarify the rich optical and magnetic behavior of noble metal nanoclusters in 1

ACS Paragon Plus Environment



ACS Photonics 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

the intermediate-size regime, but it also suggests strategies to design negative-index plasmonic metamaterials with multiple tunable resonances in the visible spectrum.

Keywords magnetic plasmons, plasmon metamolecules, plasmon oligomers, angle-resolved cathodoluminescence, retardation eects Plasmon oligomers, or metamolecules, are clusters of noble metal nanoparticles (MNPs) that exhibit anomalously strong magnetic character even though their individual nanoparticle building blocks do not. Known as magnetic plasmons, the collective cyclic oscillation of conduction-band electrons within each cluster has been shown to couple to and enhance the magnetic eld of light and to hybridize with other clusters analogously to electric plasmons in individual MNPs. 17 Since their discovery in 2004, 8 such magnetic plasmon resonances have been the focus of a body of applied and fundamental research from negative index metamaterials 9,10 to quantum tunneling, 11 sensing, 12,13 and Fano interferences. 6,1416 The magnetic properties of single oligomers, i.e., those clusters composed of only one ring of MNPs, are now well understood based upon quasistatic modeling, numerical simulation, and optical and electron-beam characterization. 11,14,17,18 Similar approaches have been employed to study toroidal dipole moments in nanoclusters. 12,16,19,20 Innitely extended two-dimensional networks of oligomers are also well understood theoretically. 21,22 However, recently, a disconnect between these extremes and methods has been discovered, 5,7 revealing a lack of understanding of how oligomer assemblies behave in the so-called intermediate size regime. 37,2325 Central to this disconnect is the question of when and how retardation eects matter. For example, in the case of the two oligomer cluster, the quasistatic approximation deviates from full-wave numerical electrodynamics simulations in that it predicts a dierent energyordering of the two magnetic plasmon resonances. 5 It is well known that retardation eects can aect the electric plasmon energetics in single MNPs, 26 dimers, 27,28 and innite chains

2

ACS Paragon Plus Environment

Page 2 of 28

Page 3 of 28 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

ACS Photonics

and arrays. 21,21,2939 However, to date, no theoretical analysis has been reported on the role of retardation in the description of magnetic plasmon resonances, which are an intrinsically relativistic phenomenon. It is the purpose of this paper to show that the incorporation of retardation corrects the aforementioned disagreement between model and simulation which is pronounced in the intermediate-size regime, and further highlights the ability to tune the magnetic plasmon spectrum by changing the oligomer size. The latter observation is particularly important because it enables manipulation of the system's eective permittivity and permeability, thereby facilitating the design of future tunable optical metamaterials. In the following, we present a coupled-dipole model to demonstrate the size-dependent evolution of the magnetic plasmons in intermediate-size N -ring oligomers, called N -mers, by including the eects of retardation in the mutual interaction of their underlying electric dipole plasmon resonances. Interestingly, as a function of size, we nd that magnetic plasmon resonances switch their ordering in analogy to the way that electric plasmons in hybridized nanoparticles 28 and the giant electric and magnetic dipole plasmons of raspberrylike metamolecules 24 switch order. We further investigate the distinct radiative properties of these magnetic plasmons by computing the angle-resolved cathodoluminescence (CL) spectrum. 4044 CL is ideal over light scattering because the selection rules of the electron probe allow for the preferential excitation of specic magnetic responses, in distinction from light scattering where many broad higher-order electric responses are also excited. Building from the intuition gained, we nally analyze the size-dependent evolution of intermediate-size plasmonic nanoclusters similar to those characterized recently in Ref. 7 We nd that the spectral ordering depends not only on aggregate size but also upon morphology, oering a second handle to modify the spectrum. Such understanding may someday lead to advances in superlensing and electromagnetic cloaking at optical frequencies, as these applications benet directly from the tunability of the system's electric and magnetic response. 4552

3

ACS Paragon Plus Environment

ACS Photonics 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 28

Model Theory

a.

b.

N

N

N

S

|B|

c.

e.

d.

N S

N S S

N N N

Figure 1: Signed magnetic eld magnitude of the 2-mer's (a,b) and 3-mer's (c,d,e) magnetic plasmon resonances. Each system supports a number of closed-loop magnetic plasmons equal to the number of rings (N ) in the oligomer. Magnetic plasmons that have a net electric dipole moment (a,c,d) show a node in their magnetic eld and have no net magnetic dipole moment. Oppositely, the nodeless magnetic modes (b,e) possess a net magnetic dipole moment that interacts with the magnetic eld of light. Fig. 1 displays two model oligomer systems, the 2-mer and 3-mer, comprising n = 10 and n = 13 silver nanospheres in analogy to the molecules naphthalene and phenalene. The electric plasmons of each MNP can be mapped onto a set of proxy harmonic oscillators as described in Ref. 5,53 This approach forms the basis for modeling the collective magnetic plasmon responses described herein. In the limit where each nanoparticle, here restricted to a spherical geometry of radius a, possesses only an electric dipole plasmon resonance, the hybridized plasmon modes of the oligomer are the solutions of the equations of motion 2 q¨ i + ωsp qi =

e X Ej (ri ), msp j6=i 4

ACS Paragon Plus Environment

(1)

Page 5 of 28 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

ACS Photonics

where ωsp =

q ωp2 /(ε∞ + 2εb ) − (γ/2)2 is the complex Mie resonance frequency and msp =

2 e2 /αsp ωsp is the plasmon's eective mass dened in terms of the static polarizability αsp =

3a3 /(ε∞ + 2εb ). 5 Here, ωp is the metal's bulk plasma frequency, ε∞ and εb are the highfrequency and background dielectric functions, and γ is the sum of the Drude electronion scattering and radiation damping rates; the latter would stem from the addition of a ... radiation reaction force in the equations of motion of the form Frad,i = (2e2 /3msp c3 ) q i . 54 While these frictional forces are not explicitly included as velocity-dependent terms in the equations of motion, friction does contribute to the red-shift and linewidth broadening of the Mie resonance with increasing MNP size. The electric dipole moment di = eqi of each plasmon oscillator is restricted to lie in the plane of the oligomer (taken to be the

x − y plane) at position ri , and Ej (ri ) is the fully retarded electric eld of the j th dipole evaluated at the position of the ith dipole. The electric eld is dened by Ej (ri ) = Λij ·   dj = (3ˆnij nˆ ij − 1ij ) 1/rij3 − ik/rij2 − (ˆnij nˆ ij − 1ij ) k2 /rij eikrij · dj /εb , where nˆ ij is the √ unit vector connecting two dipoles separated by distance rij = |ri − rj | and k = εb ω/c. −3 −2 The eld is decomposed into three partsthe near- (rij ), intermediate- (rij ), and far-eld −1 (rij )which together include all dipolar retardation eects. 55 Interestingly, the latter term

carries the opposite sign of the former two. As will be shown in the following, it is because of both the relative magnitude and sign of these terms as well as the oscillatory behavior of

eikrij that the interactions can change character and switch from energy-lowering to energyraising. Since magnetic plasmons only involve the hybridization of (x, y)-oriented electric dipole plasmons that are orthogonal to the z -direction, only those 2n (of the 3n) electric dipole plasmons in the x − y plane are explicitly accounted for in Eq. (1). As described in the Methods Section, there are as many hybridized modes of each N mer as there are electric dipole plasmons. Of these hybrid modes, some are electric and some are magnetic in character, with each N -mer studied here having N magnetic plasmon resonances. Fig. 1 shows the magnetic plasmon modes of the 2-mer and 3-mer overlaid with their associated signed magnetic eld magnitude. Each mode is named for its particular 5

ACS Paragon Plus Environment

ACS Photonics 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 6 of 28

magnetic eld distribution after the poles of a magnet, either all-North (aN) or North-South (NS). The top right-hand corner of each panel indicates the net electric or magnetic dipole moment magnitude. The black arrows denote the electric dipole plasmons on each MNP within each magnetic mode, while the magenta arrow indicates the mode's net electric dipole moment. In planar geometries, the magnetic dipole moments are always perpendicular to the in-plane electric dipole moments. As a result, they scatter light with dierent directionality as will be discussed.

Spectral Ordering of Magnetic Plasmon Modes in IntermediateSize

N -mers

To determine the impact of incorporating retardation eects, it is rst useful to consider the quasistatic limit. Here individual nanoparticle plasmon resonance frequencies (ωsp ) are size-independent and the coupling strength depends only upon the scale sij = rij /a and not upon the overall oligomer size. That is, if this ratio sij of interparticle distance rij to nanoparticle radius a remains constant, the collective resonance frequencies do not change as the nanoparticles increase in size. This is evident by considering the ka  1 limit of the interaction on the right-hand side of Eq. (1),

e2 e Ej (ri ) = Λij · qj ka→0 msp msp 2 ωsp 3ˆ nij nˆ ij − 1ij = αsp · qj 3 εb rij  2  ωsp 3 3ˆ nij nˆ ij − 1ij = · qj . εb ε∞ + 2εb s3ij lim

(2)

Retardation eects, however, remove this invariance and introduce radius-dependence into both the Mie frequency and the electric eld. Consequently, in the following, we choose the individual MNP radius a as the independent variable. To dene the oligomer geometry, we

6

ACS Paragon Plus Environment

Page 7 of 28 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

ACS Photonics

x the nearest-neighbor distance to 3a, such that the distance between each MNP remains a xed multiple of a for any value of a. Fig. 2 shows the predicted magnetic plasmon resonance energies from both Eq. (1) and full-wave numerical electrodynamics simulation of the 2-mer and 3-mer with increasing

a. 40 At small a, the magnetic modes preserve the quasistatic ordering predicted in previous work, 5 which is to be expected when the oligomer system ts within an optical wavelength and the near-eld dominates the electric eld. As the system size increases, the resonance energies decrease and cross at a ≈ 7 nm, and once more at a ≈ 20 nm. The second crossing recovers the same spectral ordering as at small sizes, which might lead to the erroneous conclusion that the quasistatic approximation is correct again. Rather, both crossings are due entirely to retardation eects as can be seen from the total interaction energy U = P −(1/2) ij di · Ej (ri ) between each pair of electric dipole plasmons within each magnetic plasmon mode. It is the inclusion of the far-eld that causes the mode switching, as it carries the opposite sign of the near- and intermediate-eld terms. Fig. 2 shows the competition between these eld components as well as their relative contributions to the total interaction energy. Inspection of Fig. 2c reveals that the crossing points occur where the splitting in the near- and intermediate-elds is equal and opposite to that of the far-eld. Perhaps more interestingly, in any dielectric medium with εb > 1 the crossing points are shifted towards smaller scale due to wavelength contraction. This can be explained through the impact of dielectric medium upon the interaction energy U, which increases with increasing εb due mostly to the far-eld term's dependence upon k 2 . 56 While previous work has explored the oscillatory nature of hybridized electric dipole plasmon resonances and radiation damping as a function of separation in noble metal dimers, 28 our work shows that magnetic plasmon resonances also oscillate with interparticle separation and can be understood by analysis of the dierent contributions to the electric dipole-electric dipole interaction.

7

ACS Paragon Plus Environment

ACS Photonics

a. 3.6

-0.3

c.

0

-0.3

1

All-North North-South

M od el Sim ula tio n

Energy (eV)

Sim ode ula l tio n

30

3.1 1

All-North North-South

30

Energy (eV)

0 Total Interaction

1

All-North North-South

30

-0.3 0

Far Near + Intermediate

All-North North-South

1

Total Interaction

Energy (eV)

Energy (eV)

b. 0

3.6

M

Energy (eV)

3.1

Energy (eV)

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 8 of 28

30

-0.3

1

All-North North-South

30

Far Near + Intermediate

All-North North-South

1

30

Figure 2: (a) Magnetic plasmon resonance energies and (b,c) electric dipole plasmon interaction energies of the 2-mer and 3-mer. For MNP radii a . 7 nm, the magnetic modes are ordered as predicted by quasistatic theory, with the NS mode (orange) lower in energy than the aN mode (blue). This implies that for small oligomers, the quasistatic approximation is accurate. However, at sizes 7 . a . 20 nm, the magnetic plasmon resonances switch spectral order. This is due to the relative strength of the far-eld interaction (squares) in comparison to the near- and intermediate-eld contributions (circles). Finally, for a & 20 nm, the magnetic plasmon resonances switch order again. Also shown in panel a are the corresponding simulated resonance energies of each magnetic plasmon mode (dashed lines), indicating excellent agreement with the presented coupled-dipole model.

8

ACS Paragon Plus Environment

Page 9 of 28 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

ACS Photonics

Radiative Properties of Intermediate-Size

N -mers

Only those magnetic resonances that have net electric or magnetic dipole moments will couple strongly to the radiation eld. The magnetic dipole moment of the aN magnetic mode is perpendicular to the electric dipole moment(s) of the NS mode in both the 2-mer and 3-mer (see Fig. 3). This means that these modes radiate with dierent directionality. It has been predicted and shown experimentally that 1-mers scatter light anisotropically when both their magnetic and electric dipole moments are mutually excited, 6,14 as is evident in dierential scattering spectra and radiation proles. Both observations stem from the time-averaged dierential power, 57,58

# " X X c 2 dP (x) ˆ · Re B∗J (x) = r n EI (x) × dΩ 8π J I

(3)

observed at the point x = rˆ n, where I, J = 1, . . . , N label the 1-mer unit cells within each

N -mer. EI and BI are produced from the sum of eective electric (dI ) and magnetic (µI ) dipole moments, assumed to be co-located at the center of the I th ring and oscillating in time as e−iωt . In this approximation, the dierential power becomes 9

! !# " X X dP (x) ck 4 = nˆ · Re (ˆ n × dI ) × nˆ − nˆ × µI × nˆ × d∗J + (ˆn × µ∗J ) × nˆ dΩ 8π J I " 4 X ck = Re dI · d∗J − (ˆn · dI )(ˆn · d∗J ) + µI · µ∗J − (ˆn · µI )(ˆn · µ∗J ) (4) 8π IJ # ˆ · (dI × µ∗J + d∗J × µI ) . +n Notice that the rst four terms involving d and µ are due to the radiation produced by electric and magnetic dipoles, while the last term involving both is due to their interference. To build intuition, we x the magnetic dipole moment µ to point in the z -direction and the electric dipole d to point in the y -direction for all I, J. The eective electric and magnetic P ˆ and µI = dipoles can be written in terms of their constituent dipoles as dI = i dI,i y 9

ACS Paragon Plus Environment

ACS Photonics 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

(k/2)

P

i

ˆ = (knRdI,i /2)ˆz, where R is the radius of the 1-mer ring consisting of n ri × dI,i φ i

nanoparticles, and dI,i is the electric dipole moment of the ith nanoparticle within the I th unit cell. For each I , all dI,i are equal. Fig. 3 displays the dierential scattering power proles associated with the aN (a) and NS (b) magnetic plasmon resonances as well as their interference (c) for any N -mer with electric and magnetic dipoles oriented as described above. Since the aN mode has only a z -oriented magnetic dipole, its radiation pattern is orthogonal to that of the y -directed electric dipole of the NS mode. This means that asymmetry in the radiation pattern is due to the excitation of both modes. Thus, examination of the radiation directionality allows one to determine the net electric or magnetic dipole character. Such forward/backward scattering asymmetries have been predicted and demonstrated in metal-semiconductor coreshell nanoparticle assemblies, that similarly exploit the interference between electric and magnetic plasmons to direct light. 59 Recent studies have demonstrated the utility of angle-resolved CL spectroscopy to characterize plasmon resonances in noble metal nanoparticles. 4244 Using an electron microscope tted with a parabolic mirror, emitted CL radiation can be collected across a large portion of the backward scattering hemisphere. 4244 Here we show that with the selection rules imposed by the location of the electron beam and the specic directionality of light emission, the size-dependent spectral order of the N -mers' magnetic plasmon modes can be observed in angle-resolved CL spectra. The previous section built up intuition on the directionality of light scattered by the magnetic modes of the 2-mer and 3-mer. This intuition informs the particular angles at which to collect CL spectra. Fig. 4 shows simulated point angle-resolved CL spectra of the 2-mer and 3-mer, with the electron beam position chosen so as to excite the aN (blue) and NS (orange) modes. Note however, due to the symmetry of the oligomers, the blue electron beam position excites not only the aN mode, but also weakly the NS mode. Since the radiation prole of the aN mode is orthogonal to that of the NS mode, it is possible 10

ACS Paragon Plus Environment

Page 10 of 28

Page 11 of 28

a.

z

b.

NS mode z

y y

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

ACS Photonics

c.

aN mode

aN and NS interference Figure 3: Dierential radiative power proles of the (a) NS and (b) aN magnetic plasmon resonances of the N -mer as well as (c) their interference as dictated by Eq. (4). In all cases, the N -mer is oriented to lie in the x − y plane so that the net electric dipole moment of the NS mode points along the y axis and the net magnetic dipole moment of the aN mode points along the z axis.

11

ACS Paragon Plus Environment

ACS Photonics

a.

b.

1

0 3.2

3.4

3.6

Normalized CL Probability

1

Normalized CL Probability

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 12 of 28

3.8

1

0 2.8

3.0

3.2

Photon Energy (eV)

3.4

0 3.2

3.4

3.6

3.8

3.0

3.2

3.4

1

0 2.8

Photon Energy (eV)

Figure 4: Angle-resolved cathodoluminescence spectra of the 2-mer (a) and 3-mer (b). Choosing the electron beam positions (×) together with the light collection angles indicated allows the individual magnetic plasmon resonances and their size-dependent spectral switching to be observed. Simulated CL spectra with a = 15 nm show the aN modes lower in energy than the NS modes, as predicted. With a = 30 nm, the modes switch in accordance with the presented coupled-dipole model. Spectral peaks located at higher energies correspond to higher-order non-magnetic plasmon modes that are not explicitly studied here.

12

ACS Paragon Plus Environment

Page 13 of 28 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

ACS Photonics

to distinguish them by analyzing the angular distribution of emitted CL radiation, even though their splitting is < 0.05 eV. The angle-resolved CL spectra at two such angles are reported in Fig. 4. These simulated spectra conrm the predictions of the model displayed in Fig. 2; specically, between a = 15 − 30 nm, the aN and NS modes switch spectral order. The simulation also shows mode switching at a ≈ 7 nm, however, the splitting energies are too small to be experimentally detectable with conventional CL. Nevertheless, experimental conrmation of the secondary mode splitting is possible, and would verify both our analytical and numerical predictions.

Magnetic Plasmon Resonances in Hexagonally-Packed Nanoclusters Ref. 7 synthesizes, fabricates, and optically characterizes hexagonally-packed nanoclusters composed of gold nanoparticles in the intermediate size regime, between few particle nanoclusters and innite arrays. Here, we apply our analytical modeling and numerical simulation to investigate the behavior of magnetic plasmon resonances in silver nanoclusters of the same geometry as a function of size, with an eye toward engineering magnetic resonances at particular resonance energies. Such tunability may be important in the design of future negative index materials with spectral tunability. Specically, we consider hexagonally packed clusters composed of 13, 19, and 31 MNPs displayed in Fig. 5 as well as the insets of Fig. 6. Fig. 5 displays the signed magnetic eld magnitude of the aN (panel a) and NS (panel b) magnetic plasmon resonances of the 13-, 19-, and 31-particle nanoclusters. The black arrows denote the electric dipole moments of the individual MNPs within each magnetic plasmon mode, the magenta arrows denote the net electric dipole moment of the nanocluster, and the background color denotes the strength of the associated magnetic eld. The net electric or magnetic dipole magnitude is denoted below each panel. The lower rows of panels a and b display the energy-ltered electron energy-loss (EEL) proles of these modes. It is important 13

ACS Paragon Plus Environment

ACS Photonics 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

a.

|B|

b.

EEL

|B|

EEL Figure 5: Signed magnetic eld magnitudes and EEL mode maps of the aN (a) and NS (b) magnetic plasmon modes of 13-, 19-, and 31-particle nanoclusters. The clusters are based on those characterized in Ref., 7 but are composed of silver. The magnetic elds are computed from the coupled-dipole model, while EEL mode maps result from simulation. As with the N -mers, only the NS modes have a net electric dipole moment indicated by the magenta arrow. The EEL mode maps indicate those regions in space where it is most probable to excite each magnetic plasmon.

14

ACS Paragon Plus Environment

Page 14 of 28

Page 15 of 28

to note that each aN mode is weakly excitable from the middle MNP (as indicated by the low EEL probability surrounding it), while the NS mode is strongly excitable from the same location (as indicated by the high EEL probability surrounding it). In all cases, the NS mode is doubly degenerate as would be true of any structure with sixfold symmetry. In analogy to the oligomers, these magnetic modes are chosen because they are the simplest (i.e., lowest order) and their radiation patterns are spatially orthogonal.

3.4 3.2 3.0 2.8

3.5 3.4

b.0.02

20

30

40

50

1

Energy (eV)

Total Interaction

0

-0.02

-0.04 -0.06

3.5 3.4 3.3

10

20

30

40

50

Total Interaction

0.10

0

-0.05

20

30

40

50

30

40

50

40

50

Total Interaction

0.10

0.05

10

1

Energy (eV)

10

All-North North-South

3.6

3.3

1

Energy (eV)

All-North North-South

3.6

Energy (eV)

All-North North-South

Energy (eV)

Energy (eV)

a. 3.6

0.05 0

-0.05

-0.08

c.

0.02

20

30

40

-0.02

-0.04 -0.06 -0.08 1

10

20

30

40

50

20

30

40

0.10 0.05 Far

0

-0.05

Near + intermediate

10

1

Energy (eV)

Far

0

50

50

10

20

30

20

0.10 0.05 Far

0

-0.05

Near + intermediate 1

10

1

Energy (eV)

10

1

Energy (eV)

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

ACS Photonics

Near + intermediate 40

50

1

10

20

30

Figure 6: Magnetic plasmon resonance energies (a) and electric-dipole interaction energies (b,c) of 13-, 19-, and 31-particle nanoclusters similar to those characterized in Ref., 7 but composed of silver. Due to the geometrical dierences between these nanoclusters and the N -mers, the aN modes remain lowest in energy at small MNP radii a. However, as a increases, one or even multiple magnetic plasmon resonance crossings become possible for the same reasons described earlier for the N -mers. Fig. 6 displays the evolution of the aN and NS magnetic resonances with size a as 15

ACS Paragon Plus Environment

ACS Photonics 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

dictated by the coupled-dipole model. Here the size-dependent behavior of the nanocluster resonances is dierent due to the hexagonal packing as opposed to the open-ring structure of the N -mer. Nevertheless, the model is capable of predicting the magnetic spectrum. Most surprising is the oscillatory behavior of the nanoclusters' magnetic resonances with size as opposed to the expected monotonic red-shifting of those of the N -mer (see Fig. 2). Panels b and c display the total interaction energy U , and its near-, intermediate-, and far-eld contributions, showing that this oscillatory behavior is due to retardation eects through the eikrij factor in the electric eld of each electric dipole. Fig. 7 shows simulated angle-resolved CL spectra for the 13-, 19-, and 31-particle nanoclusters acquired at the two indicated beam positions and collection angles. Inspection of the EEL maps of the NS and aN magnetic resonances in Fig. 5 indicates the optimal electron beam positions to excite each magnetic plasmon mode, i.e., those positions with highest EEL probability. In particular, the NS mode (Fig. 5b) is most easily excited near the central MNP, while the aN mode (Fig. 5a) is more easily excited from the outer MNPs surrounding the central MNP. Upon selecting the appropriate excitation location, the radiative prole of each mode further helps identify the two magnetic resonances of interest. Because the aN mode radiates strongly in the x − y plane, and weakly in the z -direction, it is detectable at angles between 0 < θ ≤ π/2 for all φ. Alternatively, by alignment with the x-axis, the NS mode is detectable at angles between 0 ≤ θ < π/2 and 0 < φ < π, and π < φ < 2π. The aN beam position also weakly excites the degenerate NS modes, meaning the corresponding CL spectrum will provide a secondary signature of the spectral location of that mode. This fact is evident in multiple simulated CL spectra displayed in Figs. 4 and 7, demonstrating again the ability of angle-resolved CL to distinguish between resonances split by . 0.05 eV. It should be noted that synthesis methods like those employed in Ref. 7 result in nearly spherical nanoparticles. However, if electron-beam lithography is used, the resulting nanoparticles will be more disklike in shape. This will have the eect of contributing a net redshift to the collective, in-plane plasmonic modes of the nanocluster, 60 while having negligible impact on 16

ACS Paragon Plus Environment

Page 16 of 28

Page 17 of 28

Normalized CL Probability

b.

3.2

3.4 3.6 Energy (eV)

3.8

1

0 3.0

c.

3.2

3.4 3.6 Energy (eV)

3.8

1

0 3.0

3.2

3.4 3.6 Energy (eV)

3.8

Normalized CL Probability

0 3.0

Normalized CL Probability

1

Normalized CL Probability

Normalized CL Probability

a.

Normalized CL Probability

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

ACS Photonics

1

0 2.00

2.50 3.25 Energy (eV)

3.75

2.50 3.25 Energy (eV)

3.75

2.50 3.25 Energy (eV)

3.75

1

0 2.00 1

0 2.00

Figure 7: Simulated angle-resolved cathodoluminescence spectra of the (a) 13-, (b) 19-, and (c) 31-particle nanoclusters from Ref., 7 but composed of silver. Blue (orange) spectra are acquired at the beam positions marked with a blue (orange) × to preferentially excite the aN (x-polarized NS) mode. The dashed lines in each spectrum indicate the resonance locations of the aN (blue) and NS (orange) modes. Each panel also displays the light collection angles associated with each spectrum. Again, all unlabeled resonance peaks correspond to higherorder plasmon modes of either electric or magnetic character, which are not studied herein.

17

ACS Paragon Plus Environment

ACS Photonics 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 18 of 28

the qualitative behavior of their energy ordering.

Conclusion In conclusion, we have presented a coupled-dipole model of MNP oligomers and hexagonallypacked nanoclusters that includes the eects of retardation and use it to show that the spectral ordering of their magnetic plasmon resonances can be controllably tuned as a function of cluster size and morphology. Full-wave numerical simulation is used to conrm the prediction that magnetic plasmon resonances can even be made to switch energy order, offering a relatively facile way to modify the system's magnetic responses in an experiment. The presented analysis is especially important in elucidating the intermediate-size regime in between individual single-ring and extended two-dimensional nanocluster networks, where retardation eects can signicantly modify the spectrum. We further study the angulardependence of the power emitted from the nanoclusters and demonstrate that the aN and NS magnetic plasmon modes, as well as their size- and morphology-dependent spectral order, can be detected in angle-resolved CL. Taken together, our work not only claries the rich optical and magnetic behavior of noble metal nanoclusters in the intermediate-size regime, but also suggests a strategy to design negative-index plasmonic metamaterials with multiple tunable resonances in the visible spectrum.

Methods Coupled-dipole model.

Refs. 5,53 describe how to map the solutions of Maxwell's equations

for electric plasmon resonances onto mechanical oscillators with eective masses. Based on this approach, the resulting Hamiltonian for a system of n frictionless coupled oscillators is

H=

X p2 1 e X i 2 2 + msp ωsp qi − qi · Ej (ri ), 2m 2 2 sp i ij(i6=j)

18

ACS Paragon Plus Environment

(5)

Page 19 of 28 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

ACS Photonics

where qi and pi are the coordinate and momentum of the ith plasmon oscillator located at

ri with eective mass msp and resonance frequency ωsp , dened previously. Ej (ri ) is the full electric dipole eld. The resulting equations of motion are, 2 q¨ i = −ωsp qi +

e X Ej (ri , ω) msp j6=i

(6)

which is Eq. (1). Assuming sinusoidal oscillation in time and rewriting the electric eld in terms of the plasmon coordinate results in the equations of motion 2 (ωsp − ω 2 )qi −

e2 X Λij (ω) · qj = 0. msp j6=i

(7)

It is this system of equations for the plasmon oscillator coordinates that we solve to determine the magnetic plasmon resonances. Considering n = 2 identical, coupled, collinear oscillators as an example and projecting out the directional-dependence, Eq. (7) becomes 2 (ωsp − ω 2 )q1 − g12 (ω)q2 = 0 2

(8)

2

(ωsp − ω )q2 − g12 (ω)q1 = 0, or, equivalently,

     2 2 ωsp − ω −g12 (ω)  q1  0    =  . 2 2 −g12 (ω) ωsp − ω q2 0

(9)

−3 −2 ikr12 where g12 = −2e2 [r12 − ikr12 ]e /εb msp . Note that the far-eld contribution to g12 van-

ishes here due to the assumed collinear arrangement of the dipoles. Solution of this system produces the hybridized plasmon resonance frequencies and modes of the dimer. Specically, they are

ω± =

q 2 ± g (ω) ωsp 12

19

ACS Paragon Plus Environment

(10)

ACS Photonics 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 20 of 28

and

1 q± = √ (q1 ± q2 ) . 2

(11)

Note that the coupling term g12 (ω) is frequency dependent, meaning that the eigenvalue problem must be solved iteratively until convergence. The magnetic plasmon resonances described herein are obtained as a straightforward extension of this dimer example. Simulation details.

Optical and electron beam simulations were performed using the MNPBEM

package. 40,41 Each spherical nanoparticle was discretized into 144 surface elements, each assigned the following Drude model parameters for silver: ε∞ = 3.77, ~γ = 0.05 eV, and

~ωp = 9.15 eV. The electron kinetic energy was chosen to be 200 keV and all calculations were performed in vacuum, i.e., εb = 1. The angle-resolved CL spectra were computed using standard functions within MNPBEM together with numerical angular integration through

±25◦ about the angles (θ, φ) indicated in the gures.

Acknowledgement This work was supported by the U.S. Department of Energy Basic Energy Sciences under award number DE-SC0018040 (D.J.M.) and by the State of Washington through the University of Washington Clean Energy Institute and via funding from the Washington Research Foundation (N.P.M.). This work was facilitated through the use of advanced computational, storage, and networking infrastructure provided by the Hyak supercomputer system and funded by the STF at the University of Washington. N.P.M. would like to thank Dr. Niket Thakkar and Harrison Goldwyn for helpful discussions and advice.

References (1) Liu, H.; Genov, D. A.; Wu, D. M.; Liu, Y. M.; Steele, J. M.; Sun, C.; Zhu, S. N.; Zhang, X. Magnetic plasmon propagation along a chain of connected subwavelength

20

ACS Paragon Plus Environment

Page 21 of 28 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

ACS Photonics

resonators at infrared frequencies.

Phys. Rev. Lett.

2006, 97, 243902.

(2) Liu, H.; Genov, D. A.; Wu, D. M.; Liu, Y. M.; Liu, Z. W.; Sun, C.; Zhu, S. N.; Zhang, X. Magnetic plasmon hybridization and optical activity at optical frequencies in metallic nanostructures.

Phys. Rev. B

2007, 76, 073101.

(3) Liu, N.; Mukherjee, S.; Bao, K.; Brown, L. V.; Dorfmüller, J.; Nordlander, P.; Halas, N. J. Magnetic plasmon formation and propagation in articial aromatic molecules. Nano Lett.

2011, 12, 364369.

(4) Liu, N.; Mukherjee, S.; Bao, K.; Li, Y.; Brown, L. V.; Nordlander, P.; Halas, N. J. Manipulating magnetic plasmon propagation in metallic nanocluster networks. Nano

ACS

2012, 6, 54825488.

(5) Cherqui, C.; Bigelow, N. W.; Vaschillo, A.; Goldwyn, H.; Masiello, D. J. Combined tight-binding and numerical electrodynamics understanding of the STEM/EELS magneto-optical responses of aromatic plasmon-supporting metal oligomers. tonics

ACS Pho-

2014, 1, 10131024.

(6) Cherqui, C.; Wu, Y.; Li, G.; Quillin, S. C.; Busche, J. A.; Thakkar, N.; West, C. A.; Montoni, N. P.; Rack, P. D.; Camden, J. P. .; Masiello, D. J. STEM/EELS imaging of magnetic hybridization in symmetric and symmetry-broken plasmon oligomer dimers and all-magnetic fano interference.

Nano Lett.

2016, 16, 66686676.

(7) Greybush, N. J.; Liberal, I. n.; Malassis, L.; Kikkawa, J. M.; Engheta, N.; Murray, C. B.; Kagan, C. R. Plasmon resonances in self-assembled two-dimensional Au nanocrystal metamolecules.

ACS Nano

2017, 11, 29172927.

(8) Sarychev, A. K.; Shalaev, V. M. Magnetic resonance in metal nanoantennas. SPIE

2004, 5508 .

21

ACS Paragon Plus Environment

Proc.

ACS Photonics 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 22 of 28

(9) Alù, A.; Salandrino, A.; Engheta, N. Negative eective permeability and left-handed materials at optical frequencies.

Opt. Express

2006, 14, 15571567.

(10) Alù, A.; Engheta, N. Dynamical theory of articial optical magnetism produced by rings of plasmonic nanoparticles.

Phys. Rev. B

2008, 78, 085112.

(11) Scholl, J. A.; García-Etxarri, A.; Aguirregabiria, G.; Esteban, R.; Narayan, T. C.; Koh, A. L.; Aizpurua, J.; Dionne, J. A. Evolution of plasmonic metamolecule modes in the quantum tunneling regime.

ACS Nano

2016, 10, 13461354.

(12) Ye, Q. W.; Guo, L. Y.; Li, M. H.; Liu, Y.; Xiao, B. X.; Yang, H. L. The magnetic toroidal dipole in steric metamaterial for permittivity sensor application. Phys.

Scripta

2013, 88, 055002. (13) Chen, J.; Fan, W.; Zhang, T.; Tang, C.; Chen, X.; Wu, J.; Li, D.; Yu, Y. Engineering the magnetic plasmon resonances of metamaterials for high-quality sensing.

Opt. Express

2017, 25, 36753681. (14) Sheikholeslami, S. N.; García-Etxarri, A.; Dionne, J. A. Controlling the interplay of electric and magnetic modes via Fano-like plasmon resonances.

Nano Lett.

2011,

,

11

39273934. (15) Hentschel, M.; Dregely, D.; Vogelgesang, R.; Giessen, H.; Liu, N. Plasmonic oligomers: The role of individual particles in collective behavior.

ACS Nano

2011, 5, 20422050.

(16) Wang, J.; Wang, N.; Hu, J.; Jiang, R. Toroidal dipole-induced absorption and scattering dip in (dielectric core)@(plasmonic shell) nanostructures.

Opt. Express

2017,

,

25

2893528945. (17) Brandl, D. W.; Mirin, N. A.; Nordlander, P. Plasmon modes of nanosphere trimers and quadrumers.

J. Phys. Chem. B

2006, 110, 1230212310.

22

ACS Paragon Plus Environment

Page 23 of 28 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

ACS Photonics

(18) Darvishzadeh-Varcheie, M.; Guclu, C.; Capolino, F. Magnetic nanoantennas made of plasmonic nanoclusters for photoinduced magnetic eld enhancement. Phys.

Rev. Appl.

2017, 8, 024033. (19) Ögüt, B.; Talebi, N.; Vogelgesang, R.; Sigle, W.; van Aken, P. A. Toroidal plasmonic eigenmodes in oligomer nanocavities for the visible. Nano

Letters

2012, 12, 52395244.

(20) Watson, D. W.; Jenkins, S. D.; Ruostekoski, J.; Fedotov, V. A.; Zheludev, N. I. Toroidal dipole excitations in metamolecules formed by interacting plasmonic nanorods. Rev. B

Phys.

2016, 93, 125420125432.

(21) Haynes, C. L.; McFarland, A. D.; Zhao, L.; Van Duyne, R. P.; Schatz, G. C.; Gunnarsson, L.; Prikulis, J.; Kasemo, B.; Käll, M. Nanoparticle optics: the importance of radiative dipole coupling in two-dimensional nanoparticle arrays.

J. Phys. Chem. B

2003, 107, 73377342. (22) Weick, G.; Woollacott, C.; Barnes, W. L.; Hess, O.; Mariani, E. Dirac-like plasmons in honeycomb lattices of metallic nanoparticles.

Phys. Rev. Lett.

2013, 110, 106801.

(23) Qian, Z.; Hastings, S. P.; Li, C.; Edward, B.; McGinn, C. K.; Engheta, N.; Fakhraai, Z.; Park, S.-J. Raspberry-like metamolecules exhibiting strong magnetic resonances. Nano

ACS

2015, 9, 12631270.

(24) Li, C.; Lee, S.; Qian, Z.; Woods, C.; Park, S.-J.; Fakhraai, Z. Controlling magnetic dipole resonance in raspberry-like metamolecules.

J. Phys. Chem. C

2018, 122, 6808

6817. (25) Parker, J.; Gray, S.; Scherer, N. F. Optical magnetism in core-satellite nanostructures excited by vector beams.

Proc. SPIE

2018, 10541 .

(26) Turner, M. D.; Hossain, M. M.; Gu, M. The eects of retardation on plasmon hybridization within metallic nanostructures.

New J. Phys.

23

2010, 12, 083062.

ACS Paragon Plus Environment

ACS Photonics 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 24 of 28

(27) Nordlander, P.; Oubre, C.; Prodan, E.; Li, K.; Stockman, M. I. Plasmon hybridization in nanoparticle dimers.

Nano Lett.

2004, 4, 899903.

(28) Dahmen, C.; Schmidt, B.; von Plessen, G. Radiation damping in metal nanoparticle pairs.

Nano Lett.

2007, 7, 318322.

(29) Clippe, P.; Evrard, R.; Lucas, A. A. Aggregation eect on the infrared absorption spectrum of small ionic crystals.

Phys. Rev. B

1976, 14, 17151721.

(30) Aravind, P.; Nitzan, A.; Metiu, H. The interaction between electromagnetic resonances and its role in spectroscopic studies of molecules adsorbed on colloidal particles or metal spheres.

Surf. Sci.

1981, 110, 189  204.

(31) Kottmann, J. P.; Martin, O. J. F. Retardation-induced plasmon resonances in coupled nanoparticles.

Opt. Lett.

2001, 10961098.

(32) Prodan, E.; Radlo, C.; Halas, N. J.; Nordlander, P. A hybridization model for the plasmon response of complex nanostructures.

Science

2003, 302, 419422.

(33) Prodan, E.; Nordlander, P. Plasmon hybridization in spherical nanoparticles. Chem. Phys.

J. of

2004, 120, 54445454.

(34) Rechberger, W.; Hohenau, A.; Leitner, A.; Krenn, J.; Lamprecht, B.; Aussenegg, F. Optical properties of two interacting gold nanoparticles.

Opt. Commun.

2003,

,

220

137141. (35) Bouhelier, A.; Bachelot, R.; Im, J. S.; Wiederrecht, G. P.; Lerondel, G.; Kostcheev, S.; Royer, P. Electromagnetic interactions in plasmonic nanoparticle arrays. J. Phys. Chem. B

2005, 109, 31953198.

(36) Myroshnychenko, V.; Rodriguez-Fernandez, J.; Pastoriza-Santos, I.; Funston, A. M.; Novo, C.; Mulvaney, P.; Liz-Marzan, L. M.; Garcia de Abajo, F. J. Modelling the optical response of gold nanoparticles.

Chem. Soc. Rev.

24

ACS Paragon Plus Environment

2008, 37, 17921805.

Page 25 of 28 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

ACS Photonics

(37) Davis, T. J.; Vernon, K. C.; Gómez, D. E. Designing plasmonic systems using optical coupling between nanoparticles.

Phys. Rev. B

2009, 79, 155423.

(38) Kinnan, M. K.; Chumanov, G. Plasmon coupling in two-dimensional arrays of silver nanoparticles: II. Eect of the particle size and interparticle distance.

J. Phys. Chem.

2010, 114, 74967501.

C

(39) Kravets, V. V.; Pinchuk, A. O. Near-, middle-, and far-eld dipolar interactions in gold nanoparticle arrays.

Proc. SPIE

2016, 9724, 97240B97240B8.

(40) Hohenester, U.; Trügler, A. MNPBEM - A matlab toolbox for the simulation of plasmonic nanoparticles.

Comp. Phys. Comm.

2012, 183, 370381.

(41) Hohenester, U. Simulating electron energy loss spectroscopy with the MNPBEM toolbox.

Comp. Phys. Comm.

2014, 185, 11771187.

(42) Coenen, T.; Vesseur, E. J. R.; Polman, A.; Koenderink, A. F. Directional emission from plasmonic Yagi-Uda antennas probed by angle-resolved cathodoluminescence spectroscopy.

Nano Lett.

2011, 11, 37793784.

(43) Coenen, T.; Vesseur, E. J. R.; Polman, A. Angle-resolved cathodoluminescence spectroscopy.

Appl. Phys. Lett.

2011, 99, 143103.

(44) Coenen, T.; Bernal Arango, F.; Femius Koenderink, A.; Polman, A. Directional emission from a single plasmonic scatterer.

Nat. Comm.

2014, 5, 3250.

(45) Pendry, J. B. Focus issue: negative refraction and metamaterials.

Opt. Express

2003,

, 639639.

11

(46) Fang, N.; Lee, H.; Sun, C.; Zhang, X. Sub-diraction-limited optical imaging with a silver superlens.

Science

2005, 308, 534537.

(47) Cai, W.; Chettiar, U. K.; Kildishev, A. V.; Shalaev, V. M. Optical cloaking with metamaterials.

Nature Photonics

2007, 1, 224. 25

ACS Paragon Plus Environment

ACS Photonics 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 26 of 28

(48) Pinchuk, A. O.; Schatz, G. C. Metamaterials with gradient negative index of refraction. J. Opt. Soc. Am. A

2007, 24, A39A44.

(49) Shalaev, V. M. Transforming light.

Science

2008, 322, 376379.

(50) Valentine, J.; Zhang, S.; Zentgraf, T.; Ulin-Avila, E.; Genov, D. A.; Bartal, G.; Zhang, X. Three-dimensional optical metamaterial with a negative refractive index. Nature

2008, 455, 376379.

(51) Ferrari, J. A.; Perciante, C. D. Superlenses, metamaterials, and negative refraction. Opt. Soc. Am. A

J.

2009, 26, 7884.

(52) Chen, S.; Zhang, Y.; Shih, T.-M.; Yang, W.; Hu, S.; Hu, X.; Li, J.; Ren, B.; Mao, B.; Yang, Z.; Tian, Z. Plasmon-induced magnetic resonance enhanced Raman spectroscopy. Nano Lett.

Article ASAP,

(53) Cherqui, C.; Thakkar, N.; Li, G.; Camden, J. P.; Masiello, D. J. Characterizing localized surface plasmons using electron energy-loss spectroscopy. Ann. Rev. Phys. Chem. 2016, , 331357.

67

(54) Draine, B. T.; Goodman, J. Beyond Clausius-Mossotti - Wave propagation on a polarizable point lattice and the discrete dipole approximation.

Astrophys. J.

1993, 405,

685697. (55) Purcell, E. M.; Pennypacker, C. R. Scattering and absorption of light by nonspherical dielectric grains.

Astrophys. J.

1973, 186, 705714.

(56) Jain, P. K.; El-Sayed, M. A. Noble metal nanoparticle pairs: eect of medium for enhanced nanosensing. (57) Jackson, J. D.

Nano Lett.

2008, 8, 43474352. , 3rd ed.; Wiley: New York, NY, 1999.

Classical Electrodynamics

(58) Schwinger, J.; Deraad, L.; Milton, K.; Tsai, W.; Norton, J. Advanced book program; Avalon Publishing, 1998. 26

ACS Paragon Plus Environment

Classical Electrodynamics

;

Page 27 of 28 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

ACS Photonics

(59) Liu, W.; Miroshnichenko, A. E.; Neshev, D. N.; Kivshar, Y. S. Broadband unidirectional scattering by magneto-electric core-shell nanoparticles. ACS

Nano

2012, 6, 54895497.

(60) Noguez, C. Surface plasmons on metal nanoparticles: The inuence of shape and physical environment.

J. Phys. Chem. C

2007, 111, 38063819.

27

ACS Paragon Plus Environment

ACS Photonics 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 28 of 28

Graphical TOC Entry

For Table of Contents Use Only Tunable Spectral Ordering of Magnetic Plasmon Resonances in Noble Metal Nanoclusters Nicholas P. Montoni, Steven C. Quillin, Charles Cherqui, and David J. Masiello Magnetic plasmon nanoclusters oer the ability to harness both the electric

and

magnetic

components

of

light.

If

incorporated

into

plasmon-enhanced solar-light harvesting devices, such spectrally tunable nanostructures would more eciently capture photons across the optical spectrum than traditional electric-plasmon-based architectures.

28

ACS Paragon Plus Environment