Metasurfaces with Electric Quadrupole and Magnetic Dipole Resonant

Apr 13, 2018 - Collective resonances in plasmonic nanoparticle arrays with electric dipole moment oriented along the lattice wave propagation are theo...
6 downloads 7 Views 2MB Size
Subscriber access provided by UNIV OF DURHAM

Metasurfaces with electric quadrupole and magnetic dipole resonant coupling Viktoriia E. Babicheva, and Andrey B. Evlyukhin ACS Photonics, Just Accepted Manuscript • DOI: 10.1021/acsphotonics.7b01520 • Publication Date (Web): 13 Apr 2018 Downloaded from http://pubs.acs.org on April 14, 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 32 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

Metasurfaces with electric quadrupole and magnetic dipole resonant coupling Viktoriia E. Babicheva∗,† and Andrey B. Evlyukhin‡ †ITMO University, St. Petersburg 197101, Russia ‡Laser Zentrum Hannover e.V., Hollerithallee 8, D-30419 Hannover, Germany E-mail: [email protected]

Abstract Collective resonances in plasmonic nanoparticle arrays with electric dipole moment oriented along the lattice wave propagation are theoretically investigated. The role of electric quadrupole (EQ) and magnetic dipole (MD) moments of gold nanoparticles in the resonant features of the arrays is analyzed. We perform both semi-analytical calculations of coupled multipole equations and rigorous numerical simulations varying contributions of the electric and magnetic multipoles by changing particle size and shape (spheres and disks). The arrays in homogeneous and non-homogeneous environments are considered. We find that even very weak non-resonant EQ and MD moments of a single particle are significantly enhanced in the periodic lattice at the wavelength of collective (lattice) resonance excitation. Importantly, we show that in the infinite arrays, the EQ and MD moments of nanoparticles are coupled and affect each other resonant contributions. We also demonstrate that at the lattice-resonance wavelength, the enhanced EQ and MD moments have contributions to reflection comparable to the electric dipole one resulting in a satisfaction of the generalized Kerker condition, reflection suppression, and strong magnetoelectric coupling in the structure.

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

Keywords. plasmonic nanoparticles, lattice resonances, Kerker effect, Rayleigh anomaly, magnetoelectric coupling, induced transparency Recently, periodic arrays of plasmonic nanoparticles have gained a lot of interest as they can produce narrow collective plasmon resonances in light transmission spectra because of the nanoparticle coupling in the array.1-15 These resonances appear at the wavelength of scattered light close to the wavelength corresponding to the Rayleigh anomaly. The strength and width of the resonances are characterized by the spectral distance between the closest Rayleigh anomaly and the single-particle plasmon resonance. The wavelengths corresponding to the Rayleigh anomaly are determined by the array periods, and the wavelength of the single-particle plasmon resonance is determined by nanoparticle size and shape. The spectral sensitivity to the environmental properties of collective resonances makes them very attractive for sensing applications.16 Resonant properties can also be utilized to enhance the performance of light sources,13 nanolasers,17 modulators,18 solar cells,19,20 and narrowband photodetectors,21 and the enabled functionalities are discussed in detail in the recent review.22 Lattice resonances appear in different shapes and have a broad range of quality factors. In the proximity to the resonance of single particle, Rayleigh anomaly significantly modifies the shape of the initial peak, which may result in its increase, decrease, splitting etc.4,8,10 When Rayleigh anomaly is far away and at the longer wavelength than the single-particle resonance, lattice resonance appears to stand alone and narrows down. As explained in the Ref. 8, the peak width is determined by the imaginary part of inverse particle polarizability, while the intersection of the real part of inverse polarizability with lattice sum defines the position of lattice resonance. Various nano-scatterers have been shown to enhance lattice resonances for orthogonal and parallel plasmonic-photonic coupling,23-26 the dipole coupling of several particles in the cell has been studied,27 and multipolar interactions in the surface-lattice resonances in twodimensional arrays of spheres has been analyzed.28,29 In polarization where an electric field of

2

ACS Paragon Plus Environment

Page 2 of 32

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

the incident wave is perpendicular to direction with a periodicity of interest, resonances are attributed to electric dipole-dipole coupling between nanoparticles, and they are well studied. No resonant feature can appear in this dipole approximation, homogeneous environment, and polarization where the electric field of the incident wave is parallel to the direction with a periodicity of interest.3 Thus, recently observed lattice resonance in parallel polarization30 cannot be attributed to the pure dipole-dipole coupling between nanoparticles and more studies are necessary to understand properties of such structures. In a more general case when excitations of magnetic multipoles are possible, e.g. in high-index dielectric nanoparticles, the lattice resonances appear in parallel polarization even for the dipole approximation, that is when only magnetic dipole (MD) resonances of each nanoparticle are excited.9,10 In this work, we restrict ourselves to the case of metal nanoparticles, where MD resonances are weak and do not provide resonant feature without coupling to other multipoles. According to the first Kerker condition,31 if electric and magnetic polarizabilities of a nanoparticle are equal each other in magnitude and phase, light scattering from this nanoparticle is suppressed in the backward direction (Kerker effect). For silicon spherical nanoparticle array, electric dipole (ED) and MD resonances do not overlap, and only non-resonant Kerker effect is possible,9,32-36 but starting with the pioneering work,37 a number of studies suggest obtaining spectral overlap of ED and MD resonances using all-dielectric nanoparticles with a more complex shape, such as disks, cubes, cones, or pyramids.34,38-41 Furthermore, Kerker conditions have been generalized to include contributions of different multipole modes and achieve directional scattering with plasmonic,42,43 dielectric,44 or hyperbolic-material particles.45,46 Recently, it has been shown that resonant Kerker effect is possible upon overlap of ED and MD lattice resonances excited in the periodic array of either silicon or core-shell nanoparticles,10 and spectral position of the resonant directional scattering can be controlled by the lattice periodicity.10,47 In this work, we study lattice resonances in periodic arrays of plasmonic nanoparticles analyzing different shapes (spheres and disks) and sizes of the nanoparticles, that allows

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

varying contributions of ED- and electric quadrupole (EQ) moments, as well as different refractive indices of the surrounding environment (Fig. 1). Small plasmonic nanoparticles (with the radius being smaller than 60 nm) have only ED resonant contribution to the scattering in the visible spectral range, and the contributions of their EQ and MD moments to the scattering are negligibly small.8 However, we show that even in the case of the homogeneous environment, small nanoparticles with small EQ and MD moments enable excitation of well pronounce lattice resonance due to EQ and MD coupling between nanoparticles in the array. Moreover, overlapping EQ and MD resonances with the broad ED response outside its lattice resonance, one can observe a suppression of total reflectance. We show that this decrease of reflectance is achieved because of Kerker condition of the directional scattering is satisfied. We derive its generalized form which takes into account ED, EQ, MD, and their coupling in the lattice. The array particle multipolar moments show magnetoelectric coupling around Rayleigh anomaly wavelength as the strong magnetic response is induced by the electric resonance. y H

x

Superstrate ns

z

E Au

k py

2R px 1st order

Substrate n 0th order

Figure 1: Schematic view of the nanosphere array with different refractive indices of substrate and superstrate (n and ns , respectively). Nanoparticles with radius R are arranged into the array with the periods px and py . Light incidence k is normal to the substrate, and electric field E is along the x-axis (referred as parallel polarization throughout the work). Both sphere and disk shapes of nanoparticles are of interest. On the one hand, from a theoretical point of view, spheres are the most favorable, as they can be described by Mie theory, and exact solutions can be obtained by analytical approach only. On the other hand, for practical applications, the realization of a periodic array of disks is easier than an array of 4

ACS Paragon Plus Environment

Page 4 of 32

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

spheres, as disks can be fabricated by conventional optical lithography. Furthermore, in the present work, we demonstrate that for the non-homogeneous environment of nanoparticle array (substrate and superstrate are of different indices), electromagnetic coupling between nanoparticles causes a band of transparency in the transmission spectrum. We analyze the properties of the nanoparticle array on the substrate with different refractive indices and show that redistribution of energy in different diffraction orders causes peaks and dips in transmission and reflection spectra. From the practical point of view, we consider the small difference in indices of substrate and superstrate. In experiments, fabricated nanostructures are often covered with another medium to achieve index matching of the substrate39,48 and obtain conditions close to the homogeneous environment in order to observe the strong lattice resonances. The first part of the paper is devoted to the investigation of EQ-MD lattice resonances in homogeneous surrounding medium performed both numerically and analytically for gold nanosphere array. Through comparison of both calculations, we show that EQ-MD coupling is critical for correct calculations of the wavelength where generalized Kerker condition of directional scattering is satisfied. The second part considers an array of disks and besides homogeneous environment, we also study an influence of a small dielectric contrast between substrate and superstrate on transmission and reflection properties of the array.

Methods and results Numerical study of nanosphere array Similar to the Ref. 30, we study infinite gold nanoparticle array being on top of the substrate with index n and covered with another material-superstrate-with refractive index ns = 1.47 (fixed throughout the paper) in parallel polarization (Fig. 1). We perform analysis using several approaches: Mie theory for single spherical particles, decomposed discrete dipole approximation37 for single disk particles, the dipole-quadrupole model for analytical investi5

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

gations of sphere nanoparticle arrays, as well as CST Microwave Studio commercial software for general numerical simulations. In this work, sphere radius is varied, i.e. R = 50, 85, or 100 nm; and disk radius is fixed to Rd = 85 nm, and the height H is either 50 or 100 nm. Gold permittivity is taken from experimental data.49 Let us consider nanoparticles in the homogeneous environment. Total scattering cross sections of single spherical nanoparticles, calculated in numerical simulations, are shown in Fig. 2a. Multipole analysis allows associating the resonant maxima with resonant excitation of ED and EQ moments of the nanoparticles (Fig. 2b). For the single sphere with R = 50 nm, ED resonance is at the wavelength 600 nm, and there are no other multipole resonances as the particle is too small (Fig. 2a, red line). For the single sphere with R = 85 nm, ED resonance is at wavelength 720 nm, and EQ resonance is at wavelength 580 nm (Fig. 2a, blue line). In the absorption cross-section of the single particle, profile features are dominated by EQ resonance as absorption in the dipole is relatively weak (Fig. 2c,b). In Supporting Information Fig. S1e,f, we show the results of calculation of an absorption and forward-tobackward scattering ratio for an array of 5 × 4 particles in x- and y- directions, respectively with distances px = 510 nm and py = 250 nm. An array with larger particles of R = 100 nm shows the increased both absorption and forward-to-backward scattering ratio for the wavelength range 700 - 800 nm indicating stronger collective response already for small array (Fig. S1e,f), and further we study infinity-large arrays in detail. We perform simulations for infinite nanoparticle array with periods px = 510 nm and py = 250 nm in the x- and y-directions, respectively. Arranging nanoparticles in the array with px = 510 nm in the medium with n = 1.47 is expected to cause lattice resonances spectrally close to the Rayleigh anomaly (λRA = npx = 750 nm in free space). Simulations show that either for large enough nanoparticles or in the case of a small mismatch between the indices of substrate and superstrate, a strong resonant feature is observed (see for example Fig. 3). We note that periodicity in the perpendicular direction, that is py , cannot cause dipole lattice resonances at the wavelength of interest because λRA = npy = 368 nm of the

6

ACS Paragon Plus Environment

Page 6 of 32

Page 7 of 32

(a) 2

2 × 10 5 Total scattering cross section, nm (1.47/1.47)

(b)

10 5

2

EQ

1.5

Sphere R = 85 nm, 1.47/1.47 ED MD EQ MQ Total

1.5 Sph 100 nm

EQ

ED

ED

0.5

1

ext

1

, nm 2

ED

0.5

Sph 85 nm Sph 50 nm

0 500

600

700

800

900

0 500

1000

600

Wavelength λ , nm

(c) 6

10 4

700

800

Wavelength

(d)

Spheres, 1.47/1.47

6

10 4

900

Sphere, R = 85 nm, 1.47/1.47 Total ED MD EQ MQ EOC

, nm 2

2

4

4

a

2

0 500

1000

, nm

R = 85 nm R = 100 nm

a , nm

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

2

600

700

800

Wavelength

900

0 500

1000

600

, nm

700

Wavelength

800

900

1000

, nm

Figure 2: Scattering and absorption properties of single gold sphere. (a) Total scattering cross-section of the nanospheres σscat under consideration; numbers denote the radius of the spheres. ED denotes electric dipole, EQ denotes electric quadrupole, and no magnetic dipole MD or quadrupole MQ resonances are excited for the single particle of such shape and dimensions.(b) Extinction cross-section σext of the sphere with R = 85 nm and its multipole decomposition calculated with Mie theory. (c) Absorption cross-section σa calculated with CST Microwave Studio. (d) Absorption cross-section σa calculated with Mie theory of the sphere with R = 85 nm. EOC is electric octupole. In all cases, the sphere is surrounded by a homogeneous environment with n = ns = 1.47 denoted as ’1.47/1.47’. See more results in Supporting Information for multipole decompositions for spheres with R = 50 and 100 nm (Fig. S1).

7

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

corresponding Rayleigh anomaly, and it is smaller than the wavelength of dipole resonance of the single particle.4 Thus, the appearance of such pronounced band of transparency motivates us to study the effect in more detail. Here we consider the only normal incidence of the waves, and the case of oblique incidence is beyond the scope of our current work. For the nanoparticle array, we analyze (i) the zeroth diffraction order of transmission/reflection, that is the propagation of the wave into the substrate/superstrate at normal direction; (ii) the first diffraction order of transmission/reflection; and (iii) full transmission/reflection of the array. We also calculate absorbance of nanoparticle array, i.e. the portion of energy absorbed inside nanoparticles. For the nanosphere array with R = 50 nm (Fig. 3a, red line), the resonant feature at Rayleigh anomaly is not pronounced. In contrast to small particles, in absorbance profile of particles with R = 85 nm, one can see a strong narrow peak close to λRA = npx = 750 nm (Fig. 3a, blue line), which corresponds to the excitation of the collective (lattice) resonance. Comparing the single particle scattering cross-section with the absorbance profiles of the array of such nanoparticles, one can see that in the array, EQ resonance appears at wavelength 540 nm while EQ of the single particle is at the wavelength 580 nm. The same strong narrow peak close to λRA = npx = 750 nm is observed for spheres with R = 100 nm, and the lattice resonance in absorption spectra is more pronounced than for spheres with R = 85 nm. As it follows from the symmetry of the system, under normal incidence with electric field E along the x-axis, the electric dipole moment of nanoparticles is excited in the x-direction. In this case, neither dipole coupling normal to the array nor coupling in the direction of dipoles can contribute to lattice mode. In other words, in the homogeneous environment and electric dipole approximation, no resonant feature related to lattice mode can appear, and no feature is observed in numerical simulations of small particles, e.g. spheres with radius R = 50 nm (Fig. 3a). Yet another process dominates for larger particles, e.g. with R = 85 or 100 nm. Even for the homogeneous environment, one can see resonant features

8

ACS Paragon Plus Environment

Page 8 of 32

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

(a)

Absorbance, 1.47 / 1.47 Sph 50 nm Sph 85 nm Sph 100 nm

0.6

(b) In resonance, l= 750 nm EQ distribution

(c) Out of resonance, l= 753 nm ED distribution

0.5 0.4 0.3 0.2 0.1 0 500

600

700 800 Wavelength λ , nm

900

1000

Figure 3: (a) Absorbance of the nanosphere array in the homogeneous environment with n = ns = 1.47. Periods are px = 510 nm and py = 250 nm. (b) In-resonance (λ = 750 nm) and (c) out-of-resonance (λ = 753 nm) field distributions for the array of nanospheres with R = 85 nm (one unit cell, side view). The calculations are performed with CST Microwave Studio. EQ field distribution is well pronounced for the resonant wavelength, and the field distribution drastically changes to ED with the small change of the wavelength. appear at λ > λRA = npx = 750 nm. As one can see from the field distribution in Fig. 3b,c, these resonances correspond to the lattice modes existed with EQs. Comparing inresonance (λ = 750 nm) and out-of-resonance (λ = 753 nm) field distributions for the array of nanospheres with R = 85 nm, one can see that EQ field distribution is well pronounced for the lattice resonant wavelength (Fig. 3b), and the field drastically transforms to ED distribution upon the small spectral shift (Fig. 3c). At the wavelength of lattice resonance, MD resonances are also excited, but their field distribution does not have a pronounced pattern, and it is difficult to distinguish them from other multipoles. The pronounced lattice resonances for the spheres with R = 85 or 100 nm cause strong resonant profiles in reflectance and transmittance spectra of the arrays (Fig. 4). For R = 85 nm, reflectance minimum close to zero is observed. For R = 100 nm, both reflectance maximum and minimum appear. It means that the excited EQ resonance is strong enough to provide destructive interference with ED and to suppress backscattering: at the wavelength of EQ lattice resonance, the reflection spectrum drops to zero which means that Kerker-like effect occurs. In the next subsection, we perform analytical calculations taking into account ED, EQ, and MD contributions and show that indeed, the observed effect results from the 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

Page 10 of 32

satisfying generalized Kerker condition with all three multipoles. (a)

(b)

Reflectance, 0th order, 1.47 / 1.47 Sph 50 nm Sph 85 nm Sph 100 nm

0.8

0.8

0.6

0.6

0.4

0.4

0.2

0.2

0 500

600

700

800

900

Transmittance (total), 1.47 / 1.47

1

1

1000

0 500

Sph 50 nm Sph 85 nm Sph 100 nm

600

Wavelength λ , nm

700 800 Wavelength λ , nm

900

1000

Figure 4: Lattice resonance in an array of nanospheres: (a) Zero-order reflectance and (b) total transmittance. The calculations are performed with CST Microwave Studio. Periods are px = 510 nm and py = 250 nm. The array is surrounded by a homogeneous environment with n = ns = 1.47 denoted as ’1.47/1.47’. See more simulation results in Supporting Information including the case of the non-homogeneous environment, when the refractive indices of the superstrate and substrate are ns = 1.47 and n = 1.5, respectively (Figs. S1-S3).

Dipole-quadrupole model In the dipole-quadrupole approximation, the spherical particles are considered as point electric and magnetic dipoles and electric quadrupoles with scalar dipolar and quadrupolar polarizabilities αp , αm , and αq , respectively. The vectors of the electric dipole pj and magnetic ˆ j of arbitrary partidipole mj moments and the tensor of the electric quadrupole moment Q cle with number j in the structure of N particles are determined by the coupled equations,

10

ACS Paragon Plus Environment

Page 11 of 32 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

which are obtained by extending the equations of Refs. 8 and 9:

p

N k02 X ˆ p l i ˆ Q (Q ˆ l nlj )] = αp [gjl × ml ] + G [Gjl p + jl ε0 l6=j ck0

j

+αp E0 (rj ) , N X ˆ p ml + c [gjl × pl ] = αm k02 [εS G jl ik0 l6=j c ˆ l nlj )]] + αm H0 (rj ) , + [qjl × (Q ik0 αq [∇E0 (rj ) + E0 (rj )∇] = 2 N αq k02 X ˆ p pl ) + ( G ˆ p pl )∇j ] + [∇j (G jl jl 2ε0 l6=j

mj

ˆj Q

(1)

N αq k02 X ˆ Q (Q ˆ l nlj )) + (G ˆ Q (Q ˆ l nlj ))∇j ] + [∇j (G jl jl 2ε0 l6=j N αq k02 i X [∇j [gjl × ml ] + [gjl × ml ]∇j ] , + 2ε0 ck0 l6=j

where ∇j is the nabla operator with respect to rj , E0 (rj ) and H0 (rj ) are the incident fields at the point of the particles with number j (under condition of normal incidence, which is considered in the paper, these fields are the same for all particles in the structures), k0 is the wave number in vacuum, c = (ε0 µ0 )−1/2 is the vacuum speed of light (ε0 and µ0 are the vacuum permittivity and permeability, respectively), εS is the dielectric constant of the host surrounding medium. The dipole and quadrupole Green’s tensors of the medium without particles are (

! i 1 ≡ 1+ − Uˆ 2 kS Rjl kS2 Rjl ! ) i3 3 eikS Rjl −1 − + 2 2 nlj nlj , kS Rjl kS Rjl 4πRjl ˆp G jl

+

(2)

11

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

and ! 6 i6 i3 ˆQ ≡ + + 3 3 Uˆ G −1 − jl 2 kS Rjl kS2 Rjl kS Rjl ! ) ikS eikS Rjl i6 15 i15 1+ − 2 2 − 3 3 nlj nlj . kS Rjl kS Rjl kS Rjl 24πRjl (

+

(3)

Here Rjl = |rj − rl |, nlj = (rj − rl )/|rj − rl | is the unit vector directed from rl to rj , nlj nlj √ ˆ j is is the dyadic product, kS = k0 εS , Uˆ is the 3 × 3 unit tensor, i is the imaginary unit, Q traceless and symmetric tensor. The vectors eikS Rjl gjl = 4πRjl

ikS 1 − 2 Rjl Rjl

! Rjl

(4)

and k 2 eikS Rjl qjl = S 2 24πRjl

3i 3 1+ − 2 2 kS Rjl kS Rjl

! Rjl .

(5)

Dipole-quadrupole coupling in the infinite array In the previous works,8,9 it has been shown that, in infinite periodic 2D array with single nanoparticle in the elementary cell, ED does not couple to either EQ (Ref. 8) or MD (Ref. 9), and all nanoparticles of the array have the same induced ED, MD, and EQ moments at the normal incidence of external light waves. Here, however, we show that the coupling of EQ and MD is realized in such infinite arrays. In all derivations below, we consider monochromatic x-polarized incident light with fields (Ex , Hy , 0) and the time dependence exp(−iωt). We also work in the approximations of ˆ = Q0 (ˆ spherical particles for which p = (p0 xˆ +0ˆ y +0ˆ z ), m = (0ˆ x +m0 yˆ +0ˆ z ), and Q xzˆ + zˆxˆ), where xˆ, yˆ and zˆ are the unit vectors of the Cartesian coordinate system. At these conditions the general system of Eqs. (1) is written (after straightforward transformations) in the

12

ACS Paragon Plus Environment

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

following form

p0 =

αp [ε0 Ex (r0 ) + S1 p0 ] , ε0

k0 c S3 Q0 ] , i αq ik0 = [ε0 ikS Ex (r0 ) + S4 Q0 + S5 m0 ] , 2ε0 c

m0 = αm [Hy (r0 ) + S2 m0 + Q0

(6)

where r0 is the position of arbitrary particle with number of zero in the array (below we will consider that this particle is located in the origin of Cartesian coordinate system), the polarizabilities αp , αm , and αq are expressed through the scattering coefficients a1 , b1 , and a2 of Mie theory8,9

αp = i

S1

S2

6πε0 εS a1 , kS3

αm = i

6π b1 , kS3

αq = i

120πε0 εS a2 , kS5

 k02 X eikS rl i ≡ = 1+ 4π l6=0 rl kS rl l6=0  1 x2l i3x2l 3x2l − 2 2− 2 − + , kS rl rl kS rl3 kS2 rl4 k02

X

Gpxx (0, rl )

 i kS2 X eikS rl ≡ 1+ = 4π l6=0 rl kS rl l6=0  1 i3yl2 3yl2 y2 − 2 2 − l2 − + , kS rl rl kS rl3 kS2 rl4 kS2

X

(7)

Gpyy (0, rl )

1 X ˆ l nl0 )]y [q0l × (Q Q0 l6=0   2 ikS rl X 3i 3 2 kS e = (−xl ) 1+ − , 24πrl3 kS rl kS2 rl2 l6=0

(8)

S3 ≡

13

ACS Paragon Plus Environment

(9)

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 14 of 32

and8

S4

 6 + kS2 x2l ik02 kS X eikS rl −2 − i = 24π l6=0 rl2 kS rl 12 + 7kS2 x2l 12 + 27kS2 x2l + i kS2 rl2 kS3 rl3  60x2 i60x2 − 2 4l − 3 5l , kS rl kS rl

+

(10)

and

S5 = 6S3 =

X x2 k 2 eikS rl  l

l6=0

S

4πrl3

3i 3 + 2 2 −1 − kS rl kS rl

 .

(11)

Most importantly, the coefficients S3 and S5 in Eqs. (6) are not equal to zero for infinite arrays providing coupling between MD (m0 ) and EQ (Q0 ) moments. Solving of the system Eqs. (6) one can find effective polarizabilities of the particles in the array. These polarizabilities eff/coup

are determined by the expressions: αpeff = p0 /Ex ; αm

eff/coup

= m0 /Hy ; αq

= 2Q0 /(ikS Ex ).

Thus we obtain 1 eff/coup

αm

1 eff/coup

αq

eff S3 αqeff k02 /(2ε0 ) 1 − S5 αm , = eff αm [1 + S3 αqeff k02 /(2ε0 )]

(12)

eff 1 − S5 αm S3 αqeff k02 /(2ε0 ) = , eff ] αqeff [1 + S5 αm

(13)

1 1 S1 = − , eff αp αp ε0 where 1 1 = − S2 , eff αm αm

1 1 S4 = − . eff αq αq 2ε0

14

ACS Paragon Plus Environment

(14)

Page 15 of 32 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

Reflection and transmission in the dipole-quadrupole model Estimation of the reflection and transmission coefficients can be obtained if we consider total electric field in the far-field region for z < 0 and z > 0.

E = E0 + Ep + Em + Eq

(15)

The electric fields generated by ED, MD, and EQ of nanoparticles are k02 Ex αpeff (Grxx , 0, 0) ε0

Ep =

Em = − Eq =

ik0 eff/coup r Hy αm (gz , 0, 0), cε0

k02 ikS Ex eff/coup Q,r αq (Gx , 0, 0), ε0 2

(16)

(17) (18)

respectively. For the case of wavelength larger than the array lattice periods, non-zero total electric field component Exf in the far field approximation is Exf

k02 eff r eff/coup r α G − ikS αm gz ε0 p xx  k02 eff/coup ikS Q,r Gx α + ε0 q 2 

= Ex eikS z +

(19)

where9 Grxx =

i 2SL kS

e∓ikS z ,

gzr = ±

1 ∓ikS z e 2SL

and (the calculation method as in Ref. 9)

GQ,r =∓ x

1 ∓ikS z e , 12SL

SL is the area of the lattice unit cell and the upper sing corresponds to z < 0 and the lower sing for the case when z > 0. In this approach the reflection (for electric field) and

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

Page 16 of 32

transmission coefficients are (compare with Ref. 9)   ikS 1 eff k02 eff/coup eff/coup r= α − αm − α , 2SL ε0 εS p 12ε0 q

(20)

  k02 eff/coup ikS 1 eff eff/coup + . t=1+ α + αm α 2SL ε0 εS p 12ε0 q

(21)

The intensity reflection R0 and transmission T0 coefficients are R0 = |r|2 ,

T0 = |t|2 .

The expressions above are derived for an infinitely large array (which makes particle identical), but an implementation of numerical calculations with Eqs. (7)-(11) include a summation over the finite number of particles, and this may result in artifact oscillations in the spectral profiles. To decrease such numerical artifacts, we perform summation for different array sizes and find result as an average of that; in particular, we use array dimensions a × a, where a is the randomly assigned number between 500 and 2000, and we average over 70 arrays. eff/coup

Calculations of the contributions αm

eff/coup

and k02 αq

/(12ε0 ) show that their magni-

tudes at resonance are comparable to αpeff /(ε0 εS ) (Fig. 5a,b). One can also note that without eff coupling, αm does not possess any resonant features. In the case of coupling, both EQ and

MD lattice resonances are excited spectrally close to Rayleigh anomaly, while without coupling (or in the case of zero MD response) EQ lattice resonance is shifted further away from RA similar to well-known ED lattice resonance shifts. Magnetic dipole induced by EQ indicates the magnetoelectric coupling in the array. In Fig. 5d,e, we demonstrate an enlarged view of spectra for spheres with R = 100 nm, and we show that reflectance minimum and transmittance maximum almost coincide with the maximum of absorbance, where EQ-MD lattice resonance is excited. Comparison of analytical calculations with Eqs. (20), (21) with numerical modeling performed with CST Microwave Studio demonstrates a very good quan-

16

ACS Paragon Plus Environment

Page 17 of 32

titative agreement (Fig. 5d,e). It also shows that taking into account EQ-MD coupling is critical for obtaining the correct profile and position of Kerker effect (Fig. 5d). The phase change of transmittance coefficient experiences a rise in the proximity of lattice resonance (Fig. 5c). Thus, we observe resonant Kerker effect enabled by destructive interference between the ED resonance of a single particle in the array and the EQ-MD lattice resonance of the particle with the non-resonant contribution of EQ and MD moments.

m3 eff 2 k /12], q 0

0.015

eff p,m

or

0.010

0.005

760

770

Wavelength [nm]

(d)

|t|2 ,|r| 2

1

0 740

750

100 80 60 40 740

770

|t|2 no coup. |r| 2 no coup.

0.5

|t|2 with coup. |r| 2 with coup.

750

760

770

750

760

770

Wavelength [nm]

Wavelength [nm]

Analytics, nanospheres R = 100 nm

0 740

760

(e) Numerical modeling, sph. R =100 nm Absorb., Reflect., Transmit.

750

Analytics, nanospheres R = 100 nm

120

Im[

-0.005 740

p-eff m-eff m-eff/coup q-eff q-eff/coup

Phase of t, deg.

m3

0.005

eff p,m or

0.010

eff 2 q k 0 /12],

Analytics, nanospheres R = 100 nm

0

(c)

Analytics, nanospheres R = 100 nm

(b)

(a)

Re[

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

Refl. Transm. Abs. RA

0.5

0 740

Wavelength [nm]

750

Wavelength

760

770

, nm

eff/coup

Figure 5: Nanosphere array. (a) Real and (b) imaginary parts of αm and 2 eff/coup k0 αq /(12ε0 ) with and without coupling of EQ and MD (see Eqs. (12) and (13)). Legend is the same for both (a) and (b) panels. Polarizability of the ED is not affected by other multipoles (see Eq. (14)). RA line denotes the wavelength of Rayleigh anomaly. (c) Change of phase of transmittance coefficient t calculated by Eq. (21). (d) Analytical calculations of reflectance and transmittance for the cases without EQ and MD coupling ’no coup.’ and with EQ and MD coupling ’with coup.’ (e) Reflectance, transmittance, and absorbance calculated in numerical modeling with CST Microwave Studio. One can see that the reflectance minimum and transmittance maximum almost coincide with the maximum of absorbance, where EQ lattice resonance is excited. The nanoparticle array with spheres of R = 100 nm is in the homogeneous environment with n = ns = 1.47. Periods are px = 510 nm and py = 250 nm. Analytical calculations are performed with the dipole-quadrupole model.

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

Nanodisk array Considering nanodisks with dimensions comparable to spheres (Rd = 85 nm and H = 50 or 100 nm), first, we analyze their scattering and absorption cross-sections under normal light incidence (Fig. 6a-c). Comparison of total scattering, extinction, and absorption cross-sections shows good agreement between calculations with decomposed discrete dipole approximation37,50 and numerical simulations with CST Microwave Studio (Fig. 6a,c). Multipole decomposition of scattering cross-sections for single disks with decomposed discrete dipole approximation37,50 confirms that for the thick disk with H = 100 nm, the total scattering cross-sections have a broad ED-resonant peak at wavelengths 800-900 nm and a narrow EQ-resonant peak at 630 nm wavelength as well as minuscule MD peak at 650 nm (Fig. 6b). For the disk with the smallest height H = 50 nm, only ED resonance is present, and EQ and MD contributions are very weak (Fig. S5 in Supporting Information). In contrast to spheres, excitation of different modes of nanodisk depends on light illumination conditions, and disks’ polarizability is anisotropic. Analytical calculations similar to sphere case require an extension of the model to include polarizability tensor. Therefore, for acquiring reflectance, transmittance, and absorbance spectra of the nanodisk array, here we use only numerical simulations with CST Microwave Studio. For the nanodisks with H = 50 nm arranged in the periodic array (Fig. 6d), reflectance peak and transmittance dip are at wavelength 750 nm (Fig. 6e,f), and the absorbance profile has both broad and narrow peaks (Fig. 6g). For the disks with H = 100 nm, there is also a broad dip in transmittance. In contrast, in absorbance, a narrow feature is the strongest in comparison to the spheres of all sizes and disks with H = 50 nm. Thus, even though the magnitudes of EQ and MD moments are not pronounced in comparison to ED resonance, their excitation is responsible for the lattice resonance (Fig. 6h,i). However, its spectral width is very narrow, which corresponds to extremely long propagation length of the excited lattice mode.6 Similar to nanospheres array, lattices of nanodisks demonstrate a suppressed reflection at the EQ-MD lattice resonance (Figs. 6e,j), which corresponds to destructive interference of EQ, MD, and 18

ACS Paragon Plus Environment

Page 18 of 32

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

ED scattering in a backward direction mimicking the resonant Kerker effect. The reflectance and transmittance profiles are similar to those of nanospheres (compare Figs. 6j and 5e).

Substrate-induced transparency Further, we study nanoparticle array in a non-homogeneous environment that is a structure where particles (disks with Rd = 85 nm and H = 100 nm) are placed on the substrate with varied n, while for the superstrate, ns = 1.47 is fixed (Fig. 7). In this case, we also observe the collective resonances, and due to the differences between the diffraction conditions in the substrate and superstrate, the resonant signature in the transmission and absorption spectra is split into two features. One can see that the transmittance spectra include narrow transmission bands between two Rayleigh anomalies that correspond to the substrate and superstrate. Strong absorption in the structure exists in homogeneous surrounding and disappears with increasing of difference between n and ns (Fig. 7a). To clarify the reason for the narrow transmittance band in the case of non-homogeneous systems, we analyze the role of EQ resonance excitation at each Rayleigh anomaly in the spectral profiles (Fig. 7e,f). Weak EQ resonances are excited in the cases of n = 1.4 and n = 1.6 at the exact position of Rayleigh anomalies, but well pronounced EQ is observed for the wavelength where the reflectance is almost zero (Fig. 7f). There, reflections from two EQ and two MD resonances compensate dipole reflection, and the Kerker effect is observed in the region between resonances. We have performed numerical simulations of the experimental structure in Ref. 30: disks with Rd = 85 nm, H = 50 nm, periods are px = 510 nm and py = 250 nm, indices ns = 1.47, n = 1.5, and 20-nm ITO layer with nITO = 1.9 appears below the disks. At λ ≈ ns px ≈ 750 nm, both reflection and absorption have local maximum, and field distribution corresponds to the ED resonance (Fig. 7g). At the same time, at λ ≈ npx ≈ 765 nm, both reflection and absorption have minima, reflection is nearly zero, and field distribution corresponds to EQ resonance. This confirms that directional scattering and Kerker effect appear at the 19

ACS Paragon Plus Environment

ACS Photonics

( m2 )

( m2 )

0.1

0 500

600

700

800

Wavelength

900

0.1

0.05

0 500

1000

600

700

800

Wavelength

, nm

(e)

900

Absorption cross-section, H = 100 nm 0.15

Total ED MD EQ MQ EOC

scat

Scat, Discr. Ext, Discr. Scat, CST Ext, CST

0.05

(c)

Scattering scross-section: Total + Decomp. 0.15 H = 100 nm

( m2 )

(b)

Scattering and extinction cross-section 0.15 H = 100 nm

abs

(a)

scat, ext

0.1

0.05

Discr. CST

0 500

1000

600

700

800

Wavelength

, nm

(f)

Reflectance, 0th order, 1.47 / 1.47

1

1

0.8

0.8

0.6

0.6

0.4

0.4

900

1000

, nm

Transmittance (total), 1.47 / 1.47

(d) x

y H k

Superstrate ns

z

E py

2Rd

Au H

px 1st order

Disk H = 100 nm

Substrate n

0.2

0th order

0.2

Disk H = 100 nm

0 500

600

700

800

900

0 500

1000

600

Wavelength λ , nm

(g)

Absorbance, 1.47 / 1.47

(h) In resonance, l= 751 nm EQ distribution

0.6 0.5 0.4 Disk H = 100 nm

(i) Out of resonance, l= 758 nm ED distribution

0.2 0.1

900

1000

Numerical modeling, disks H = 100 nm

1

0 600

700 800 Wavelength λ , nm

Refl. Transm. Abs. RA

0.5

0.3

0 500

(j) Absorb., Reflect., Transmit.

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 32

700

800

Wavelength λ , nm

900

1000

745

750

755

760

Wavelength , nm

Figure 6: (a) Scattering σscat and extinction σext cross-sections for a single: comparison of numerical simulations and decomposed discrete dipole approximation calculations. (b) Total scattering cross-sections σscat and its multipole decomposition of the single nanodisks calculated using a decomposed discrete dipole approximation.37 EOC is electric octupole, and MQ is magnetic quadrupole, and the multipole moment values are near-zero and almost coincide with abscissa axis. (c) Absorption cross-sections σa for a single nanodisk: comparison of numerical simulations and decomposed discrete dipole approximation calculations. Inset in (c): Discretized nanodisk: calculations are performed with the step 7 nm. In panels (a)-(c), nanodisk has Rd = 85 nm and H = 100 nm. (d) Periodic array of the nanodisks with radius Rd and height H. (e) Zero order reflection, (f) total transmission, and (g) absorbance of the nanodisk array with Rd = 85 nm obtained by numerical calculations with CST Microwave Studio. (h) In-resonance (λ = 751 nm) and (i) out-of-resonance (λ = 758 nm) field distributions for the array of nanodisks with H = 50 nm (one unit cell, side view). EQ field distribution is well pronounced for the resonant wavelength, and the field distribution drastically changes to ED with the small change of wavelength. (j) Reflectance, transmittance, and absorbance of the nanodisk array, H = 100 nm: One can see that the reflectance and transmittance profiles are similar to those of nanospheres (Fig. 5(e)). RA line denotes the wavelength of Rayleigh anomaly. In panels (d)-(j), periods are px = 510 nm and py = 250 nm. In all cases, the nanodisk array is surrounded by a homogeneous environment with n = ns = 1.47 denoted as ’1.47/1.47’. See more simulation results in Supporting Informa20 tion including the case of the non-homogeneous environment, when the refractive indices of ACS Paragon Plus Environment the superstrate and substrate are ns = 1.47 and n = 1.5, respectively (Figs. S6 and S7).

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

(a)

(b)

Absorbance

2

n = 1.60 1.5

4

600

700

800

Wavelength λ , nm

900

1000

2

1

n = 1.40

600

700

800

900

n = 1.47

n = 1.47

1

n = 1.40

n = 1.50

n = 1.50

n = 1.47

0.5

n = 1.60

3

2

n = 1.47

Reflectance (total)

n = 1.60

n = 1.50

1

0 500

4

3

n = 1.50

0.5

(d)

Transmittance (total)

n = 1.60

1.5

1

0 500

(c)

Transmittance, 1st order

2

0 1000 500

Wavelength λ , nm

n = 1.40

n = 1.40

600

700

800

900

Wavelength λ , nm

0 1000 500

600

700

800

900

1000

Wavelength λ , nm

Figure 7: Lattice resonance in an array of nanodisks when substrate and superstrate have different refractive indices. (a) The absorbance, (b) first order transmission, (c) total transmission, and (d) total reflection of the array of nanodisks (Rd = 85 nm and H = 100 nm) for different n denoted in the plots. Periods are px = 510 nm and py = 250 nm, ns = 1.47. The dotted color lines mark a wavelength where the resonances in the substrate are expected, i.e. λRA = npx , and the dotted black line corresponds to the resonance in superstrate, i.e. λRAs = ns px . Each plot is shifted by either 0.5 or 1 with respect to the previous one. See more simulation results for this array in Supporting Information (Figs. S8). Similar behavior has been obtained for disks with H = 50 nm (Fig. S9) and other periods of the array (not shown here). Analysis of absorption spectra for the other set of parameters (px = 420 nm and 1.4 < n < 2.15) is presented in the Supporting Information (Fig. S10). (e)-(g) Field distributions at different parts of the reflectance profiles (top parts): (e) for n = 1.4; (f) for n = 1.6; and (g) corresponds to the experimental structure in Ref. 30. Quadrupole field distribution is the most pronounced for near-zero reflectance for n = 1.6 and experimental structure in Ref. 30.

21

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 22 of 32

Transmittance (total), Disk, Rd = 85 nm, H = 50 nm

(d)

1.47 / 1.5 1.47 / 1.48 1.47 / 1.47

0.6 0.5 0.4 0.3 0.2 0.1 700

750

800

850

Wavelength λ , nm

Figure 8: (a)-(c) Schematic view of transmission and reflection in the zeroth (T0 and R0 ) and the first diffraction orders (T±1 and R±1 ) for different wavelengths (shown in panels) in the case ns < n. While Rayleigh anomaly is a complex process and results in non-monotonic spectra, for simplicity, we use notations λRA = npx and λRAs = ns px . For the case ns < n and λ < λRAs , both reflection and transmission into the first diffraction order are allowed and non-zero (panel (a)). Next, for λRAs < λ < λRA , the first diffraction order of reflection is suppressed (panel (b)). Finally, for λ > λRA , both the first diffraction orders of reflection and transmission vanish (panel (c)). When the first order disappears in the superstrate, a lattice mode of the structure is excited and it couples to the substrate, which results in an increase of transmission and suppression of reflection. (d) The transmittance of the array of nanodisks with Rd = 85 nm and H = 50 nm with a small difference between the refractive indices of the substrate and the superstrate (denoted in the legend). Periods are px = 510 nm and py = 250 nm.

22

ACS Paragon Plus Environment

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

wavelength of strong EQ and MD excitations. Figure 8a-c schematically shows transmission and reflection in the zeroth and first orders of diffraction and their changes that depend on the incident wavelength with respect to Rayleigh anomaly in the substrate and the superstrate. We calculated the cases with a small difference between the refractive indices of the substrate and the superstrate, in particular, ns = 1.47 and n = 1.48. There, the substrate-induced transparency is narrower and stronger in comparison to the case n = 1.5 (Fig. 8d). The smaller the difference between indices the stronger and the narrower the transmittance peak, and it appears between the two Rayleigh anomalies associated with the substrate and superstrate. For equal indices, the transmittance peak vanishes. The substrate-induced transparency effect is essentially different from the effects of the substrate observed before.4,51-54 In contrast to the previous works, where the large difference between medium indices causes smearing out of the lattice resonances, in our work, the difference is small, resonance conditions in each medium are similar, and because of the effective "broadening", the lattice resonances are more pronounced. In the experimental work Ref. 30, substrate and index-matching superstrate refractive indices have a difference about 0.1, which is associated with fabrication process, and particularly this difference has given a rise to well-pronounce resonances features.

Conclusions It is well known that dipole coupling in one- and two-dimensional plasmonic nanoparticle arrays can produce narrow collective plasmon resonances in light transmission spectra, and the wavelengths are determined by the array periods. In the electric dipole approximation, the collective lattice resonances involve only dipole moments of the nanoparticles oriented perpendicular to the lattice wave propagation. We have studied transmittance, reflectance, and absorbance of periodic nanoparticle arrays in parallel polarization, where only EQ and MD moments are involved in lattice resonances, but the ED moment is broad and out of its

23

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

lattice resonance. It has been shown that these arrays can support excitation of lattice mode resonances in both homogeneous and non-homogeneous environments. The effect cannot be described using only electric dipole approximations, and higher multipoles coupling between nanoparticles are necessary to explain pronounced resonant features. Even for very small, but non-zero, EQ and MD moments, their lattice resonances are strong enough to provide destructive interference with the electric dipole. This results in zero reflectance associated with a resonant Kerker effect where the generalized condition of directional scattering is satisfied. The most importantly, we have found that in the infinite periodic nanoparticle lattice even in the homogeneous environment, EQ and MD moments of nanoparticles are coupled and affect each other resonant contributions. In this work we have studied the case of infinite arrays with identical nanoparticles under normal incidence of light, and one of the main results of the work is a demonstration of coupling between electric quadrupole and magnetic dipole moments resulting in lattice resonances. Frequently, the lattice resonances are very narrow, and as one can expect, the finite size of the realistic array, experimental imperfections, uncertainties, defects, and so on cause resonance broadening, decreasing, and smearing out. Previous works on plasmonic nanoparticle arrays had shown the influence of the array size on the lattice resonances profile.8,55 For instance, lattice resonances appear already in limited-number nanoparticle array, such as ∼50 particles55 and variations in angle of incidence related to the experimental measurements. Furthermore, in our previous work,10 we have shown that arrays with ∼50 particles not only support lattice resonances but also allow they overlap with resonances of other multipoles and lattice Kerker effect. Lattice resonances have also been shown to persist even upon shuffling and removing particles in the array56 and lasing enabled by lattice resonances can be observed even when 99% of nanoparticle are absent. Finally, we note that the present our work has been motivated by experimental observation of lattice resonances in parallel polarization;30 and well-pronounced lattice resonances in the experiment confirm the feasibility of the effect observation and its practical utilizing. In the finite-size arrays because of

24

ACS Paragon Plus Environment

Page 24 of 32

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

the boundary effects, the coupling between particle multipoles can be more pronounced in a broad spectral range.8 Although in the experimental realization with finite size particles arrays, the excitation of the lattice mode is difficult, the resonance broadening is possible in the case of the nonhomogeneous environment. Owing to an excitation of collective multipole resonances in the nanoparticle arrays in the non-homogeneous environment, e.g. with different substrate and superstrate indices, a narrow transparency band can appear between two Rayleigh anomalies that correspond to disappearing of the reflection and transmission first diffraction orders in the substrate and superstrate, respectively. For the small difference between indices, strong substrate-induced transparency can be achieved along with a narrow peak in absorbance profile. We have concluded that experimentally observed resonances in parallel polarization in Ref. 30 are the result of the presence of substrate and superstrate of different refractive indices and excitation of resonant lattice modes supported by the multipole coupling. In terms of practical applications, because of the high spectral sensitivity of collective resonances to the environment, these lattice resonances can be used for sensing.

Acknowledgement The authors acknowledge financial support from the Deutsche Forschungsgemeinschaft (Germany), the project EV 220/2-1. The numerical study has been supported by the Russian Science Foundation (Russian Federation), the project 16-12-10287.

Supporting Information Available Additional information is provided on transmittance, reflectance, and absorbance profiles for different parameters of the nanoparticle array. This material is available free of charge via the Internet at http://pubs.acs.org/. 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

References (1) Zou, S.; Janel N.; Schatz, G. C. Silver nanoparticle array structures that produce remarkably narrow plasmon lineshapes. J. Chem. Phys. 2004, 120, 10871. (2) Zou, S.; Schatz, G. C. Theoretical studies of plasmon resonances in one-dimensional nanoparticle chains: narrow lineshapes with tunable widths. Nanotechnology 2006, 17, 2813-2820. (3) Markel, V.A. Divergence Of Dipole Sums And The Nature Of Non-Lorentzian Exponentially Narrow Resonances In One-Dimensional Periodic Arrays Of Nanospheres. Journal of Physics B: Atomic, Molecular and Optical Physics 2005, 38, L115-L121. (4) Auguie, B.; Barnes, W.L. Collective Resonances In Gold Nanoparticle Arrays. Phys. Rev. Lett. 2008, 101, 143902. (5) Kravets, V.G.; Schedin, F.; Grigorenko, A.N. Extremely narrow plasmon resonances based on diffraction coupling of localized plasmons in arrays of metallic nanoparticles. Phys. Rev. Lett. 2008, 101, 087403. (6) Vecchi, G.; Giannini, V.; Gomez Rivas, J. Surface Modes In Plasmonic Crystals Induced By Diffractive Coupling Of Nanoantennas. Phys. Rev. B 2009, 80, 201401(R). (7) Zhou, W.; Odom, T. Tunable subradiant lattice plasmons by out-of-plane dipolar interactions. Nature Nanotechnol. 2011, 6, 423-427. (8) Evlyukhin, A.B.; Reinhardt, C.; Zywietz, U.; Chichkov B.N. Collective resonances in metal nanoparticle arrays with dipole-quadrupole interactions. Phys. Rev. B 2012, 85(24), 245411. (9) Evlyukhin, A.B.; Reinhardt, C.; Seidel, A.; Lukyanchuk, B.S.; Chichkov, B.N. Optical Response Features Of Si-Nanoparticle Arrays. Phys. Rev. B 2010, 82, 045404.

26

ACS Paragon Plus Environment

Page 26 of 32

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

(10) Babicheva, V.E.; Evlyukhin, A.B. Resonant lattice Kerker effect in metasurfaces with electric and magnetic optical responses. Laser and Photonics Reviews 2017, 11, 1700132. (11) Nikitin, A.G.; Kabashin, A.V.; Dallaporta, H. Plasmonic resonances in diffractive arrays of gold nanoantennas: near and far field effects. Opt. Express 2012, 20(25), 2794127952. (12) Rodriguez, S.R.K.; Abass, A.; Maes, B.; Janssen, O.T.A.; Vecchi, G.; Gomez Rivas, J. Coupling Bright and Dark Plasmonic Lattice Resonances. Phys. Rev. X 2011, 1, 021019. (13) Rodriguez, S.R.K.; Lozano, G.; Verschuuren, M.A.; Gomes, R.; Lambert, K.; De Geyter, B.; Hassinen, A.; Van Thourhout, D.; Hens, Z.; Gomez Rivas, J. Quantum rod emission coupled to plasmonic lattice resonances: A collective directional source of polarized light. Appl. Phys. Lett. 2012, 100(11), 111103. (14) Mousavi, S.H.; Khanikaev, A.B.; Neuner, B.; Fozdar, D.Y.; Corrigan, T.D.; Kolb, P.W.; Drew, H.D.; Phaneuf, R.J.; Alu, A.; Shvets, G. Suppression Of Long-Range Collective Effects In Meta-Surfaces Formed By Plasmonic Antenna Pairs. Optics Express 2011, 19, 22142. (15) Lozano, G.; Louwers, D.J.; Rodriguez, S.R.; Murai, S.; Jansen, O.T.; Verschuuren, M.A.; Gomez Rivas, J. Plasmonics For Solid-State Lighting: Enhanced Excitation And Directional Emission Of Highly Efficient Light Sources. Light: Science and Applications 2013, 2, e66. (16) Offermans, P.; Schaafsma, M.C.; Rodriguez, S.R.K.; Zhang, Y.; Crego-Calama, M.; Brongersma, S.H.; Gomez Rivas, J. Universal scaling of the figure of merit of plasmonic sensors. ACS Nano 2011, 5(6), 5151-5157. (17) Zhou, W.; Dridi, M.; Suh, J.Y.; Kim, C.H.; Co, D.T.; Wasielewski, M.R.; Schatz, G.C.;

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

Odom, T.W. Lasing Action In Strongly Coupled Plasmonic Nanocavity Arrays. Nature Nanotechnology 2013, 8, 506-511. (18) Thackray, B.D.; Thomas, P.A.; Auton, G.H.; Rodriguez, F.J.; Marshall, O.P.; Kravets, V.G.; Grigorenko, A.N. Super-Narrow, Extremely High Quality Collective Plasmon Resonances At Telecom Wavelengths And Their Application In A Hybrid GraphenePlasmonic Modulator. Nano Lett. 2015, 15, 3519-3523. (19) Zhukovsky, S.V.; Babicheva, V.E.; Uskov, A.V.; Protsenko, I.E.; Lavrinenko, A.V. Enhanced Electron Photoemission By Collective Lattice Resonances In Plasmonic Nanoparticle-Array Photodetectors And Solar Cells. Plasmonics 2013, 9, 283-289. (20) Zhukovsky, S.V.; Babicheva, V.E.; Uskov, A.V.; Protsenko, I.E.; Lavrinenko, A.V. Electron Photoemission In Plasmonic Nanoparticle Arrays: Analysis Of Collective Resonances And Embedding Effects. Appl. Phys. A 2014, 116, 929-940. (21) Sobhani, A.; Knight, M.W.; Wang, Y.; Zheng, B.; King, N.S.; Brown, L.V.; Fang, Z.; Nordlander, P.; Halas, N.J. Narrowband photodetection in the near-infrared with a plasmon-induced hot electron device. Nat. Commun. 2013, 4(3), 1643. (22) Wang, W.; Ramezani, M.; Vakevainen, A.I.; Torma, P.; Rivas, J.G.; Odom, T.W. The Rich Photonic World Of Plasmonic Nanoparticle Arrays. Materials Today 2017, https://doi.org/10.1016/j.mattod.2017.09.002. (23) Lin, L.; Yi, Y. Lattice plasmon resonance in core-shell SiO2/Au nanocylinder arrays. Opt. Lett. 2014, 39, 4823-4826. (24) Lin, L.; Yi, Y. Orthogonal and parallel lattice plasmon resonance in core-shell SiO2/Au nanocylinder arrays. Opt. Express 2015, 23, 130-142. (25) Lin, L.; Zheng, Y. Multiple plasmonic-photonic couplings in the Au nanobeaker arrays: enhanced robustness and wavelength tunability. Opt. Lett. 2015, 40, 2060-2063. 28

ACS Paragon Plus Environment

Page 28 of 32

Page 29 of 32 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

(26) Lin, L.; Zheng, Y. Engineering of parallel plasmonic-photonic interactions for onchip refractive index sensors. Nanoscale 2015, 7, 12205-12214. (27) Baur, S.; Sanders, S.; Manjavacas, A. Hybridization Of Lattice Resonances. ACS Nano 2018, 12, 1618-1629. (28) Swiecicki, S.D.; Sipe, J.E. Surface-lattice resonances in two-dimensional arrays of spheres: Multipolar interactions and a mode analysis. Phys. Rev. B 2017, 95, 195406. (29) Swiecicki, S.D.; Sipe, J.E. Periodic Green functions for 2D magneto-electric quadrupolar arrays: explicitly satisfying the optical theorem. J. Opt. 2017, 19, 095006. (30) Nikitin, A.G. Diffraction-induced subradiant transverse-magnetic lattice plasmon modes in metal nanoparticle arrays. Applied Physics Letters 2014, 104, 061107. (31) Kerker, M.; Wang, D.; Giles, C. Electromagnetic Scattering by Magnetic Spheres. J. Opt. Soc. Am. 1983, 73, 765. (32) Fu, Y. H.; Kuznetsov, A. I.; Miroshnichenko, A. E.; Yu, Y. F.; Lukyanchuk, B. Directional Visible Light Scattering by Silicon Nanoparticles. Nat. Commun. 2013, 4, 1527. (33) Person, S.; Jain, M.; Lapin, Z.; Saenz, J. J.; Wicks, G.; Novotny, L. Demonstration of Zero Optical Backscattering from Single Nanoparticles. Nano Lett. 2013, 13(4), 18061809. (34) van de Groep, J.; Polman, A. Designing Dielectric Resonators on Substrates: Combining Magnetic and Electric Resonances. Opt. Express 2013, 21, 26285-26302. (35) Baryshnikova, K.V.; Petrov, M.I.; Babicheva, V.E.; Belov, P.A. Plasmonic and silicon spherical nanoparticle antireflective coatings. Scientific Reports 2016, 6, 22136. (36) Babicheva, V.; Petrov, M.; Baryshnikova, K.; Belov, P. Reflection compensation mediated by electric and magnetic resonances of all-dielectric metasurfaces [Invited]. Journal of the Optical Society of America B 2017, 34, D18-D28. 29

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

(37) Evlyukhin, A. B.; Reinhardt, C.; Chichkov, B. N. Multipole Light Scattering by Nonspherical Nanoparticles in the Discrete Dipole Approximation. Phys. Rev. B 2011, 84, 235429. (38) Staude, I.; Miroshnichenko, A. E.; Decker, M.; Fofang, N. T.; Liu, S.; Gonzales, E.; Dominguez, J.; Luk, T. S.; Neshev, D. N.; Brener, I.; Kivshar, Y. Tailoring Directional Scattering through Magnetic and Electric Resonances in Subwavelength Silicon Nanodisks. ACS Nano 2013, 7, 7824-7832. (39) Decker, M.; Staude, I.; Falkner, M.; Dominguez, J.; Neshev, D. N.; Brener, I.; Pertsch, T.; Kivshar, Y. S. High-Efficiency Dielectric Huygens Surfaces. Adv. Opt. Mater. 2015, 3, 813-820. (40) Sikdar, D.; Cheng, W.; Premaratne, M. Optically Resonant Magneto-Electric Cubic Nanoantennas for Ultra-Directional Light Scattering. J. Appl. Phys. 2015, 117, 083101. (41) Terekhov, P.D.; Baryshnikova, K.V.; Artemyev, Y.A.; Karabchevsky, A.; Shalin, A.S.; Evlyukhin, A.B. Multipolar Response Of Nonspherical Silicon Nanoparticles In The Visible And Near-Infrared Spectral Ranges. Physical Review B 2017, 96, 035443. (42) Alaee, R.; Filter, R.; Lehr, D.; Lederer, F.; Rockstuhl, C. A Generalized Kerker Condition For Highly Directive Nanoantennas. Optics Letters 2015, 40, 2645. (43) Pors, A.; Andersen, S.K.H.; Bozhevolnyi, S.I. Unidirectional scattering by nanoparticles near substrates: generalized Kerker conditions. Opt. Express 2015, 23, 28808-28828. (44) Yang, Y.; Miroshnichenko, A.E.; Kostinski, S.V.; Odit, M.; Kapitanova, P.; Qiu, M.; Kivshar, Y. Multimode Directionality In All-Dielectric Metasurfaces. Physical Review B 2017, 95, 165426. (45) Babicheva, V.E. Multipole resonances and directional scattering by hyperbolic-media antennas. arXiv preprint arXiv:1706.07259, 2017. 30

ACS Paragon Plus Environment

Page 30 of 32

Page 31 of 32 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

(46) Babicheva, V.E. Directional scattering by the hyperbolic-medium antennas and silicon particles, MRS Advances, 10.1557/adv.2018.112 (2018). (47) Yang, Ch.-Y.; Yang, J.-H.; Yang, Z.-Y.; Zhou, Z.-X.; Sun, M.-G.; Babicheva, V.E.; Chen, K.-P. Nonradiating Silicon Nanoantenna Metasurfaces As Narrowband Absorbers. ACS Photonics 2018, 10.1021/acsphotonics.7b01186. (48) Yu, Y.F.; Zhu, A.Y.; Paniagua-Dominguez, R.; Fu, Y.H.; Lukyanchuk, B.; Kuznetsov, A.I. High-Transmission Dielectric Metasurface With 2π Phase Control At Visible Wavelengths. Laser and Photonics Reviews 2015, 9, 412-418. (49) Johnson, P.B.; Christy, R.W. Optical constants of the noble metals. Phys. Rev. B 1972, 6, 4370-4379. (50) Evlyukhin, A.B.; Reinhardt, C.; Evlyukhin, E.; Chichkov, B.N. Multipole Analysis of Light Scattering by Arbitrary-Shaped Nanoparticles on a Plane Surface. J. Opt. Soc. Am. B 2013, 30, 2589. (51) Vitrey, A.; Aigouy, L.; Prieto, P.; Garcia-Martin, J.M.; Gonzalez, M. U. Parallel Collective Resonances in Arrays of Gold Nanorods. Nano Lett. 2014, 14, 2079-2085. (52) Auguie, B.; Bendana, X.M.; Barnes, W.L.; Garcia de Abajo, F.J. Diffractive arrays of gold nanoparticles near an interface: Critical role of the substrate. Phys. Rev. B 2010, 82, 155447. (53) Nikitin, A.G.; Nguyen, T.; Dallaporta, H. Narrow plasmon resonances in diffractive arrays of gold nanoparticles in asymmetric environment: Experimental studies. Applied Physics Letters 2013, 102, 221116. (54) Lamprecht, B.; Schider, G.; Lechner, R.T.; Ditlbacher, H.; Krenn, J.R.; Leitner, A.; Aussenegg, F.R. Metal Nanoparticle Gratings: Influence of Dipolar Particle Interaction on the Plasmon Resonance. Phys. Rev. Lett. 2000, 84, 4721. 31

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 32 of 32

(55) Rodriguez, S.R.K.; Schaafsma, M. C.; Berrier, A.; Gomez Rivas, J. Collective resonances in plasmonic crystals: Size matters. Physica B: Condensed Matter 2012, 407(20), 4081-4085. (56) Schokker, A.H.; Koenderink, A.F. Statistics of Randomized Plasmonic Lattice Lasers. ACS Photonics 2015, 2, 1289.

For Table of Contents Use Only Journal: ACS Photonics Manuscript ID: ph-2017-01520y.R1 Title: "Metasurfaces with electric quadrupole and magnetic dipole resonant coupling" Author(s): Babicheva, Viktoriia; Evlyukhin, Andrey

Au

R = 85 nm

EQ ED

Reflectance o

EQ & MD coupling Generalized Kerker effect

600

700

800

900

Wavelength λ , nm Figure 9: The decrease of reflection from the periodic array at the lattice resonance enabled by the coupling of EQ and MD moments of the nanoparticles

32

ACS Paragon Plus Environment