Arxiv:2001.11731V2 [Physics.Optics] 15 Jul 2020 Can Be Efficiently Explored Using the Discrete Dipole Proved Accuracy of Such Approach Compared to the Con- Model14–17
Total Page:16
File Type:pdf, Size:1020Kb
Nonlocal response of Mie-resonant dielectric particles Daniel A. Bobylev,1 Daria A. Smirnova,2, 3 and Maxim A. Gorlach1 1Physics and Engineering Department, ITMO University, Saint Petersburg 197101, Russia 2Nonlinear Physics Centre, Australian National University, Canberra ACT 2601, Australia 3Institute of Applied Physics, Russian Academy of Science, Nizhny Novgorod 603950, Russia Mie-resonant high-index dielectric particles are at the core of modern all-dielectric photonics. In many situations, their response to the external fields is well-captured by the dipole model which neglects the excitation of higher-order multipoles. In that case, it is commonly assumed that the dipole moments induced by the external fields are given by the product of particle polarizability tensor and the field in the particle center. Here, we demonstrate that the dipole response of non- spherical subwavelength dielectric particles is significantly more complex since the dipole moments are defined not only by the field in the particle center but also by the second-order spatial derivatives of the field. As we prove, such nonlocal response is especially pronounced in the vicinity of anapole minimum in the scattering cross-section. We examine the excitation of high-index dielectric disk in microwave domain and silicon nanodisk in near infrared applying group-theoretical analysis and retrieving the nonlocal corrections to the dipole moments. Extending the discrete dipole model to include nonlocality of the dipole response, we demonstrate an improved agreement with full-wave numerical simulations. These results provide important insights into meta-optics of Mie-resonant non-spherical particles as well as metamaterials and metadevices based on them. I. INTRODUCTION order of magnitude as the wavelength of light at magnetic or electric dipole resonance2. Therefore, it is not obvious Over the recent years, all-dielectric nanophotonics and a priori that electric and magnetic dipole moments of an meta-optics1,2 have demonstrated a variety of exciting arbitrarily shaped particle are related only to the field in functionalities including strong directional scattering of the particle center. light3,4, flexible phase manipulation of the transmit- In this Article, we assess this important assumption ted signal with transparent metasurfaces5, high-quality and reveal that the dipole moments of the typical Mie- modes of dielectric particles6, enhanced nonlinear phe- resonant disk are governed not only by electric and mag- nomena in the arrays of resonant nanoparticles7,8 and netic fields in its center but also by their second-order precise molecular fingerprinting with all-dielectric meta- spatial derivatives which crucially determine the electro- surfaces9. The physics underlying this plethora of effects magnetic response of the disk in the vicinity of anapole is based on the combination of electric and magnetic re- minimum18 in the scattering cross-section. To high- sponses of dielectric particles10{12, where optical mag- light that the discussed physics is universally valid across netic response is associated with circular displacement the entire electromagnetic spectrum, we examine two currents excited in the scatterer. representative examples: dielectric disk made of high- A key theoretical tool to capture electromagnetic be- permittivity ceramics with resonances in the microwave havior of subwavelength objects is multipole expansion13 domain and silicon nanodisk supporting Mie resonances which presents the field scattered by the particle as a sum in the near infrared. In both cases, we focus on the of different multipoles, each of them being characterized frequency range where electric and magnetic dipole re- by the unique set of polarization and angular dependence sponses provide the dominant contribution to the scatter- of the radiation pattern. In many situations, the domi- ing cross-section separating them from the contributions nant contribution to the scattering cross-section of sub- of higher-order multipoles via multipole decomposition wavelength particles is provided by the lowest-order mul- technique. tipoles, namely, magnetic and electric dipoles, while the Furthermore, we extend the discrete dipole model by contribution of higher-order multipoles can be neglected. incorporating the dependence of the dipole moments on In that case, the physics of complex particle arrays spatial derivatives of the fields and demonstrate an im- arXiv:2001.11731v2 [physics.optics] 15 Jul 2020 can be efficiently explored using the discrete dipole proved accuracy of such approach compared to the con- model14{17. In this approach, the scatterer is viewed as ventional discrete dipole approximation. electric and magnetic dipoles placed in its center, while It should be stressed that our results are conceptu- all information on particle properties is embedded into ally different from the nonlocal effects in small metal- its electric and magnetic polarizability tensors. These lic nanoparticles manifested via size-dependent resonance tensors link the external fields acting on the particle to shifts and linewidth broadening19,20. The latter effects its dipole moments and can be retrieved from full-wave occurring due to electron-electron interactions are en- numerical simulations. hanced as the size of plasmonic nanoparticle decreases. However, dielectric particles used in experimentally In stark contrast, the mechanism we discuss here becomes relevant situations are not that deeply subwavelength, increasingly important as the size of the particle grows and the diameter of the particle has typically the same becoming of the same order of magnitude as the wave- 2 length at the dipole resonance. Kronecker symbols δij and even powers of n vector di- The rest of the article is organized as follows. In rected along the disk axis. Besides that, αij and γijlm Section II we construct a general relation between the are symmetric with respect to the first pair (i; j) of the 25 dipole moment of the particle from one side and inci- indices due to symmetry of kinetic coefficients ; γijlm is dent field with its spatial derivatives from the other us- also symmetric with respect to the last pair (l; m) of the ing the symmetry arguments. To support our analysis indices. The above requirements yield: further, in Sec. III we consider a dielectric disk made of high-permittivity ceramics in the microwave domain αij = α1 δij + α2 ninj ; (3) and a nanodisk made of crystalline silicon in the near γijlm = γ1 δij nlnm + γ2 (ninl δjm + njnm δil infrared. The parameters of both particles are chosen in +n n δ + n n δ ) + γ n n n n (4) such a way that the dipole response dominates the rest j l im i m jl 3 i j l m of multipole contributions in a sufficiently wide frequency +γ4 δijδlm + γ5 (δilδjm + δimδjl) + γ6 ninjδlm ; range. Performing full-wave numerical simulations of the disks' dipole response, we reveal the conditions for the where α1;2 and γ1−6 are some unknown scalar coefficients strongest nonlocal effects. Focusing on the case of ce- which depend on the material and shape of the particle ramic disk, we extract the components of polarizability and the frequency of excitation. tensors as well as the nonlocal corrections to the induced Combining Eqs. (3), (4) with Eqs. (1), (2), we derive dipole moments. Section IV continues with the discus- the expression for the dipole moment of the disk: sion of the generalized discrete-dipole model which in- 2 corporates nonlocal corrections to the dipole moments di = α?"0 Ei + (αk − α?)ni "0 Ez + γ1 "0 Eikz and enables an increased accuracy in the description of 2 +2γ2kikz "0 Ez + γ3ni "0 Ezk ; (5) meta-crystals composed of such Mie-resonant disks. Fi- z nally, Sec. V concludes with a summary and an outlook where we take into account that jkj = k0 = !=c and for future studies. k · E = 0 since the incident field satisfies the condi- tion div E = 0. α? and αk are the standard frequency- dependent components of polarizability tensor of an II. SYMMETRY ANALYSIS OF THE DISK 2 anisotropic particle defined as α? = α1 + γ4 k0 and RESPONSE 2 αk = α? + α2 + γ6 k0. Similar equation holds for the magnetic dipole moment: In the most general scenario, electric and magnetic dipole moments d and m of the particle induced by the (m) (m) (m) (m) 2 mi = α? Hi + (αk − α? )niHz + γ1 Hikz impinging plane wave can be presented as Taylor series (m) (m) 2 with respect to wave vector k: +2γ2 kikz Hz + γ3 niHzkz : (6) di = "0 fαij Ej + βijl Ejkl + γijlm Ejklkm + ::: g ; (1) Quite importantly, if the particle is spherical and has (m) (m) (m) full rotational symmetry, the only nonzero components of mi = αij Hj + βijl Hjkl + γijlmHjklkm + :::; (2) the tensors Eqs. (3), (4) are α1, γ4 and γ5, which ensures that the link between the dipole moment and the field −i! t where SI system of units and e time convention are remains local. th n used. n terms of both expansions scale as (R/λ) rel- In the case of a disk, second-order nonlocal corrections ative to the respective leading-order terms, where R is to the dipole moment of the particle are captured by the particle characteristic size and λ is the wavelength. the three additional terms proportional to γ1, γ2 and γ3. Therefore, assuming that the particle is subwavelength, Note that all of them exhibit characteristic dependence we keep only the first three terms in Eqs. (1), (2). on the direction of the incident wave propagation since Note that in the analysis below we define the dipole they depend on kz. Based on this observation, we con- moment based on the angular dependence of the fields, sider the geometry illustrated in Fig. 1. consistently with Refs.21,22. In the alternative formula- 23,24 TE-polarized wave [Fig. 1(a)] excites dy, mx and mz tions of multipole expansion , however, thus defined components of the dipole moments given by the equations di and mi correspond to the sum of dipole moment, toroidal moment and higher-frequency contributions.