## Abstract

Dielectric photonics platform provides unique possibilities to control light scattering via utilizing high-index dielectric nanoantennas with peculiar optical signatures. Despite the intensively growing field of all-dielectric nanophotonics, it is still unclear how surrounding media affect scattering properties of a nanoantenna with complex multipole response. Here, we report on light scattering by a silicon cubic nanoparticle embedded in lossless media, supporting optical resonant response. We show that significant changes in the scattering process are governed by the electro-magnetic multipole resonances, which experience spectral red-shift and broadening over the whole visible and near-infrared spectra as the indices of media increase. Most interestingly, the considered nanoantenna exhibits the broadband forward scattering in the visible and near-infrared spectral ranges due to the Kerker-effect in high-index media. The revealed effect of broadband forward scattering is essential for highly demanding applications in which the influence of the media is crucial such as health-care, e.g., sensing, treatment efficiency monitoring, and diagnostics. In addition, the insights from this study are expected to pave the way toward engineering the nanophotonic systems including but not limited to Huygens-metasurfaces in media within a single framework.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

## 1. Introduction

Tuning the scattering resonances of plasmonic or dielectric nanoparticles has enabled scientists and engineers to localize light at nanoscale and enhance peculiar optical phenomena over the last decade [1, 2]. Specifically, changing the shape and dimensions of high-index dielectric nanoparticles leads to the variation of their optical properties opening a door for a plethora of interesting phenomena [2–7]. Nanoparticles scatter light and support the excitation of the geometrical (Mie) electric and magnetic multipolar resonances [2–8] which in turn enhance light-matter interaction. Despite the fact that the influence of nanoparticles size, geometry, and material on their optical properties is being extensively explored [3, 9, 10] for cloaking [11] and spectroscopy [12], the influence of media in which the nanoparticles are embedded is still obscured. Here, therefore, we bridge this gap in knowledge by exploring the optical properties of high-index nanoparticles depending on the optical properties of a surrounding medium. One of the auxiliary tools for understanding the interaction of light with nanoparticles is multipole decomposition [13].

The role of multipole contributions to scattering effect is of high importance for a variety of applications, including but not limited to nanoantennas [14–18], sensors [19–21], solar cells [22], multi-functional metasurfaces [23–25,] and cloaking devices [11,26]. Dielectric nanoparticles could be used for a drug delivery, and as probes for an electromagnetic diagnosis in emerging medical applications [27–29]. The particles’ resonant behavior in dielectric media like, e.g., the human body environment is another very interesting question [30].

Forward scattering phenomenon has been a subject of extensive studies both theoretically and experimentally [31–37]. Most of such studies are based on engineering the multipolar resonances to meet the requirement of the well-known Kerker effect [31]. In a core-shell nanoparticle, for instance, the dielectric shell shifts the electric dipole resonance of the metal core so both core and shell resonances overlap and contribute to the forward scattering due to constructive interference [33]. A high-index substrate on another hand provides a leakage medium and results in a broadening of the multipoles resonances, therefore supporting broadband forward scattering and allowing a design of efficient dielectric Huygens’ metasurfaces [32,34].

Recent analytics [35,37] followed by experimental verification [36] has been carried out for low-index Mie spheres in vacuum for obtaining broadband backscattering suppression. Artificially created low-index materials enable comparable amplitudes for all multipole moments so exhibit the duality effect. However, for realistic applications in visible and near infrared frequencies, it is convenient to use the nanostructures embedded in homogeneous dielectric media.

We show that increasing the refractive index of a medium induces high-order multipole resonances in the scatterer. In addition, induced resonances experience a spectral broadening and incoherent red shift. It also allows to use smaller nanoparticles to design compact optical devices based on high-index particles.

The proposed structure supports overlapping of several multipole resonances leading to the strong asymmetric broadband forward scattering. We found that subwavelength silicon cubic particle embedded in a dielectric medium is capable to support a considerable scattering efficiency while maintaining broadband forward scattering in the visible spectrum. Insights from this study are expected to pave the way towards engineering the nanophotonic systems in media within a single framework.

## 2. Theoretical background

In this paper, we consider the optical system including the single silicon nanocube suspended in a lossless medium. Figure 1 shows the schematics of the studied system which is a nanocube with edge of *H* = 250 nm, made of polycrystalline silicon [38] and suspended in media with the indexes n of 1 ≤ *n* ≤ 2 having increment of *dn* = 0.2. Note that the wavelengths values throughout the article are considered for the incident light wave in air.

To analyze the considered system, we use the semi-analytical multipole decomposition approach reported in [39]. This approach is based on the multipole expansions frequently used in electrodynamic theory [39–41]. Here we use the expressions of the dynamic multipole moments up to the magnetic quadrupoles as defined in [41]. Note that for scatterers located in a medium the wavenumber, introducing in the expressions for multipoles [41], should be taken in this medium. In addition, we consider the quasistatic electric octupole moment as derived in [39]. We validated that the considered set of multipole moments is sufficient for the proper description of the system we study. In contradiction, the quasistatic approximation is not suitable, because long-wavelength approximation cannot be applied. The comparison between the approximation and the rigorous approach has been performed in [41]. The full electric field in the system is numerically calculated by the finite elements method (FEM) implemented in the COMSOL Multiphysics commercial package [42,43]. The calculated electric field is used for obtaining of the multipole moments and their contributions in the scattering cross sections [39–41].

#### 2.1. Asymmetry parameter

The main goal of our work is studying of the directional scattering effect in lossless
media. Such effects like Kerker-effect are of interest due to
asymmetric light scattering in the forward and backward directions.
For the convenient description of such effects, we introduce the
asymmetry parameter. The asymmetry parameter is the efficiency factor
that quantifies the directivity of light scattering. To derive and
calculate the asymmetry parameter for an arbitrarily shaped particle,
we apply the Cartesian multipole decomposition method for the
scattered field. Then, the scattered electric field
**E*** _{sca}* by a particle in a
homogenous host medium is the superposition of the induced multipole
moments contributions [39]:

**n**is in the direction of scattering vector

**r**, ${\epsilon}_{d}={n}_{d}^{2}$ is the relative dielectric permittivity of the surrounding medium,

*ε*

_{0}is the vacuum electric permittivity, ${v}_{d}=c/\sqrt{{\epsilon}_{d}}$ is the light speed in the surrounding medium and

*c*is the light speed in the vacuum;

*k*

_{0}and

*k*are the wavenumbers in vacuum and in the surrounding medium, correspondingly.

_{d}**m**is the magnetic dipole moment (MD) of a particle;

**D**is the total electric dipole moment (TED);

*Q̂*,

*M̂*and

*Ô*are the electric quadrupole moment tensor (EQ), the magnetic quadrupole moment tensor (MQ) and the tensor of electric octupole moment (OCT), respectively. Note that these tensors are symmetric and traceless and in tensor notation e.g.

*Q̂*is equal to

*Q*, where subscript indices denote components (e.g

_{αβ}*α*=

*x*,

*y*,

*z*) [39]. The scattering cross-section can be presented as (see [39] for details):

**E**

*is the electric field amplitude of the incident light wave. The total scattering cross-section is obtained through the integration of the Pointing vector over a closed surface in the far-field zone and the normalization to the incident field intensity [39]. The asymmetry parameter is the ratio of the cosine-weighted scattering cross section over the total scattering cross section*

_{inc}*C*and can be calculated as

_{sca}*d*Ω =

*sinθdθdϕ*, we find

*I*(1 −

_{inc}*g*)

*C*is the net rate of momentum transferring to the particle in the direction of the propagation where

_{sca}*I*is the irradiance of the incident light beam. Upon inspection of the asymmetry parameter formula presented above, we found that all the multipoles cross-terms scattering light symmetrically along the polar angle

_{inc}*θ*= 90° have vanished. Therefore, the asymmetry parameter is a measure of light directivity; meaning that, its zero value corresponds to the equal scattering in both half spaces. Predominating scattering to the lower half space (forward) and the upper half space (backward) results in positive and negative values, respectively. Furthermore, the aforementioned Kerker effect for dipoles and Kerker-like effects for high-order multipoles can be understood with the asymmetry parameter equations [44,45]. If we consider the case of

*E*polarized incident light, then Kerker condition of TED and MD fulfills when the ratio

_{x}*v*/

_{d}D_{x}*m*= 1 with

_{y}*g*= 0.5. Kerker-like condition of TED and EQ interplay fulfills at 6

*iD*/

_{x}*k*= 1 while for the MD and MQ interplay it fulfills when 2

_{d}Q_{xz}*im*/

_{y}*k*= 1. Consequently, the EQ and MQ interplay scatter in forward directions when

_{d}M_{yz}*v*/3

_{d}Q_{xz}*M*= 1 and finally for EQ and EO coherence the condition is 5

_{yz}*Q*/4

_{xz}*ik*= 1. Generally, the synchronous multipole couplings enhance the field directivity, and play a key role for the suppression of the backscattering. Here, the total electric fields and the corresponding induced polarization current in the scatterers are calculated numerically using COMSOL Multiphysics. Using the calculated polarization current, the multipole moments and their contributions to scattering cross-sections and the asymmetry parameter are obtained by the numerical integration.

_{d}O_{zzx}In the next section we analyze in detail the results obtained for the different surrounding media and the multipoles behavior in the system.

## 3. Results and discussion

#### 3.1. Multipoles spectra evolution

To analyze the optical properties of the considered nanoparticle in different media, we draw the multipole decomposition of the scattering cross-section spectra for the six different cases of 1 ≤ *n* ≤ 2 (see Fig. 2). To prove that our multipole decomposition approach can be applied, we compare the scattering cross-sections obtained as the sum of multipole contributions (Sum Scat) and the scattering cross-section obtained with the direct calculation (Total scat (COMSOL)) for every case. As discussed above, this comparison shows the good coincidence between the calculation methods when we use the dynamic multipole moments.

We start our analysis with the case of the particle embedded in air. The case of *n* = 1 (see Fig. 2(a)) has been considered in details in [3] for the similar nanoantenna. One can note the pronounced resonant peaks in the scattering cross-section spectrum; as can be noted from Fig. 2, these peaks are associated with the resonant excitations of the total electric dipole (TED) moment, the magnetic dipole (MD) moment, the electric quadrupole (EQ) moment and the magnetic quadrupole (MQ) moment. Figure 2(b) shows that starting from *n* = 1.2 the scattering cross-section resonant peaks start to merge with each other. At *λ* = 640 nm, the scattering cross-section peak is now due to the interaction between MD and EQ moments forming the one wider scattering peak. Similarly, the EQ and MQ moments resonant peaks begin to merge into the single peak around *λ* = 750 nm, but for *n* = 1.2 they still can be distinguished. At last, the TED and MD resonances in the near-infrared spectral range experience some broadening and become a bit smoother. However, the spectra for *n* = 1.2 is still close to the case of air medium.

For *n* = 1.4 (see Fig. 2(c)), the MQ and EQ resonant peaks merge to the single scattering cross-section peak at *λ* = 790 nm. At the same time, the TED and MD resonances in the near-infrared spectral range continue their broadening; the TED resonant excitation no longer leads to a separate scattering cross-section peak. Finally, the EQ and MQ resonances continue to merge and provide the scattering cross-section peak together at *λ* = 650 nm. In the case of *n* = 1.4 only three scattering peaks can be distinguished.

Next, we consider *n* = 1.6 (see Fig. 2(d)). The EQ moment still experiences the bigger red shift in comparison with the magnetic multipole moments MQ and MD. For *n* = 1.6, one can note the resonant excitation of electric octupole moment (OCT) at *λ* = 620 nm. We also expect the excitation of higher-order multipole moments in this wavelength region, but the moments taken into account are enough for a proper qualitative description of the system. It is worth noting that OCT contributes to the scattering peak formed by EQ and MD at *λ* = 650 nm. At the region of longer wavelengths, the spectral positions of MQ and EQ resonances match each other at *λ* = 795 nm and the resulting scattering cross-section peak is well pronounced in Fig. 2(d). The MD-induced scattering cross-section peak can be still noted at *λ* = 1160 nm. Such intermediate case still shows three scattering cross-section peaks however slightly broadened.

The final step of the study is presented in Figs. 2(e) and 2(f) for *n* = 1.8 and *n* = 2. In the wavelength range 600 ≤ *λ* ≤ 700, the resonant excitation of OCT provides the incremental contribution to the scattering cross-section. This contribution, together with the EQ and MD resonant excitations and the non-resonant MQ contribution, forms the scattering cross-section peak at the considered region. The EQ moment resonance continues its shifting to the red zone with respect to the MQ resonant area; hence, their joint scattering cross-section peak becomes smoother (*λ* = 810 nm). The MD-induced scattering peak at *λ* ≈ 1170 nm is almost smoothed out for *n* = 2. The multipole analysis performed in this way is very important to understand the origins of scattering process. Obtained insights are necessary to study broadband forward scattering enchantment and explain it in detail.

To provide the better visibility of the multipole evolution in Fig. 3, and to improve the consequent analysis of the scattering asymmetry we show the evolution of each multipole contributions as *n* changes. The evolution of the dipole moment contributions is shown in Figs. 3(a) and 3(b) for the TED and MD correspondingly; the evolution of the quadrupole moment contributions is shown in Figs. 3(c) and 3(d) for the EQ and MQ correspondingly. Figure 3(a) clearly shows that the TED resonant area becomes broader, but the peak value decreases as *n* rises. Moreover, the TED resonant region experiences a relatively strong redshift as n rises. Similarly, both MD resonances at *λ* ≈ 660 nm and *λ* ≈ 1150 nm (see Fig. 3(b)) experiences broadening, but the maximum contribution decreases as n rises. In addition, both resonant areas experience the little redshift.

Interestingly, the quadrupole resonances (see Figs. 3(c) and 3(d)) behave slightly differently. As shown in Fig. 3(c), both EQ resonances experience the redshift; it is especially strong for the resonant excitation at longer wavelengths, where the EQ resonant peak shifts from *λ* = 700 nm for *n* = 1 to *λ* = 810 nm for *n* = 2. The maximum value of the contribution to the scattering cross-section for this resonance decreases slower in comparison with TED moment as *n* rises. Both EQ resonant regions experience some broadening. Further, the MQ moment has the only one resonant area for the considered system (see Fig. 3(d)). This resonance experiences the noticeable redshift, but it is weak in comparison with the redshift of the EQ moment. The MQ resonance shows the broadening as *n* rises, but the maximum value of the MQ contribution decreases. However, this decrease is slower in comparison with MD moment. Thus, we can associate the behavior of certain multipole moments with a change in refractive index of the media.

Such comparative analysis allows the study of the evolution of each multipole contribution to the scattering cross-section and to compare their behavior with each other. First, we prove that electric multipole resonances experience a stronger redshift than their magnetic counterparts as *n* increases. Next, we show that the relative contributions of quadrupole moments to the scattering cross-section get stronger with *n*. These insights can be used to design specific combinations of the multipole moments for tuning the direction of the light scattering.

We note that the spectral contributions of the multipole moments significantly change with *n*. In addition, the separated peaks merge, hence, the scattering peaks for higher *n* smooth out and broaden, which, as we will show in the next section, leads to the significant amplification of the forward scattering in a broad wavelength range. This enables us to design novel types of nanoantennas and antireflective coatings to apply in high-index surroundings and provides a room of opportunities for the design of new Huygens-type metasurfaces [32].

#### 3.2. Far-field scattering and electric field distribution inside the particle

In this section, we discuss the physical limitations of the forward scattering
enhancement as refractive index contrast between the scatterer and the
host media is reduced. The optical theorem states that the total
energy extinct by a single nanoparticle is in direct relation to the
scattering amplitude in the forward direction [46]. Based on this definition we
underline the energy conservation and physical limitations with the
following explanations. The wave impedance of a medium is
(*Z*_{0}/*n*) where
*Z*_{0} is the wave impedance in free space and
hence high refractive index surrounding medium provide better coupling
for the optical modes leaking out of the particle. Consequently, as
the index contrast between the scatterer and the host media reduces,
the inner reflection from the scatterer-medium interface decreases
according to the Fresnel coefficients [13]. As a result, the induced multipole
moments broadened and become less confined inside the nanoparticle.
Figure 4 illustrates this
effect by plotting the evolution of the field distribution both inside
the nanocube and in the far-field region. The upper row (i.e. Figs. 4(a)–4(c)) in Fig. 4 shows the evolution of the
far-field scattering pattern at the wavelength corresponding to the
first MQ peaks in spectrum for the cases *n* =
1, 1.6 and 2; in turn, the middle row (i.e. Figs. 4(d)–4(f)) presents the field
distribution inside the particle for the same cases. By selecting
these points, we show the smooth transformation of the far-field
radiation pattern. The lower row (i.e. Figs. 4(g)–4(i)) in Fig. 4 shows the near-field distribution in
(*xz*) plane for a better visualization of the leakage
and spreading of the mode to the surrounding dielectric media. In the
case of air medium MQ resonant contribution leads to the simultaneous
scattering in four directions (see Fig. 2); however, the radiation pattern transforms as n rises
and exhibits the well-pronounced forward scattering effect with
*n* = 2. Here we analyse the evolution of the
radiation pattern to emphasize the huge difference between two extreme
cases.

In Figs. 5(a) and 5(b), we present the scattering efficiency (Eq. (2)) normalized to geometrical cross-section of the cube) for the three cases of *n* = 1, *n* = 1.6 and *n* = 2 along with the asymmetry parameter (Eq. (4)). As expected, the multipole broadening afore-explained in Fig. 3 causes the synchronous overlapping of the multipoles of different orders and leads the asymmetry parameter to be dramatically broadened and enhanced as the surrounding medium refractive index increases. Balanced dipole moments fulfilling the Kerker’s condition provide the asymmetry parameter *g* = 0.5; however, when the quadrupole contributions are significant, the asymmetry parameter becomes *g* > 0.6, indicating strong enhancement of the forward directivity. To clarify this point, in Fig. 5(c), we show the 2D radiation patterns for 6 different wavelengths over the spectrum for *n* = 2 and a strong asymmetry remains for all spectral points. Thereby, the forward scattering amplification in the broad range is the direct consequence of the high refractive index of the surrounding medium; the similar nanoantenna does not show such properties in case of *n* = 1 [3]. Importantly, we note that the scattering efficiency for the cubical particle is also enhanced over the short wavelength optical range where the dipole, quadrupole and higher order multipoles contribute. However, for the long wavelength region, the scattering decreases since the higher order multipoles contribution is not enough to compensate the leaked dipoles. This result can be widely implemented to develop nanoantennas and other optical devices for applications in non-air media like liquids and dielectric materials. In addition, the broadening of the multipole moments especially the TED is of an importance for sensing applications. The weakening of the field inside the particle leads the multipole moments to be highly sensitive and exhibit redshift when the surrounding medium refractive index changes.

## 4. Conclusion

To conclude, we reported on the broadband forward scattering effect due to the evolution of multipole moments for silicon cubical nanoparticles embedded in lossless media. We found that the electric multipole moments (TED, EQ) experience stronger redshift compared to the magnetic multipole moments (MD, MQ) as the index of medium increases. In addition, the separate scattering cross-section peaks transform to smoother merged peaks as the index of medium increases. We discovered that with increasing the index of surrounding media, the broadband Kerker effect appears. We confirmed this effect through deriving the asymmetry factor for Cartesian multipoles system. We showed that the scattering intensity of the silicon nanoparticles is highly sensitive to the surrounding medium refractive index. Therefore, the studied system can work as an efficient sensor and is expected to be superior to the wavelength-shift depending plasmonic sensors. Also, the silicon nanocube as a building block in high-index media is expected to be applicable for the design of Huygens metasurfaces and transparent devices operating in the broadband spectral range.

## Funding

Israel Innovation Authority-Kamin Program (62045); Russian Foundation for Basic Research (RFBR) (18-02-00414, 18-52-00005); Ministry of Education and Science of the Russian Federation (GOSZADANIE 3.4982.2017/6.7, 16.7162.2017/8.9); Russian Science Foundation (16-12-10287); Government of the Russian Federation (08-08); Deutsche Forschungsgemeinschaft under Germany’s Excellence Strategy within the Cluster of Excellence PhoenixD (EXC 2122, 390833453).

## Acknowledgments

The research described was performed by Pavel Terekhov as a part of the joint Ph.D. program between the BGU and ITMO. The development of the analytical approach and the calculations of multipole moments have been partially supported by the Russian Science Foundation (16-12-10287).

## References

**1. **A. Karabchevsky, A. Mosayyebi, and A. V. Kavokin, “Tuning the
chemiluminescence of a luminol flow using plasmonic
nanoparticles,” Light: Sci.
Appl. **5**, e16164
(2016). [CrossRef]

**2. **A. I. Kuznetsov, A. E. Miroshnichenko, M. L. Brongersma, Y. S. Kivshar, and B. Lukyanchuk, “Optically resonant dielectric nanostructures,” Science **354**, 2472 (2016). [CrossRef] [PubMed]

**3. **P. D. Terekhov, K. V. Baryshnikova, Y. A. Artemyev, A. Karabchevsky, A. S. Shalin, and A. B. Evlyukhin, “Multipolar response of nonspherical silicon nanoparticles in the visible and near-infrared spectral ranges,” Phys. Rev. B **96**, 35443 (2017). [CrossRef]

**4. **A. B. Evlyukhin, C. Reinhardt, A. Seidel, B. S. Lukyanchuk, and B. N. Chichkov, “Optical response features of Si-nanoparticle arrays,” Phys. Rev. B **82**, 45404 (2010). [CrossRef]

**5. **A. I. Kuznetsov, A. E. Miroshnichenko, Y. H. Fu, J. Zhang, and B. Luk’yanchuk, “Magnetic
light,” Sci. Rep. **2**, 492 (2012). [CrossRef]

**6. **A. B. Evlyukhin, S. M. Novikov, U. Zywietz, R. L. Eriksen, C. Reinhardt, S. I. Bozhevolnyi, and B. N. Chichkov, “Demonstration of magnetic dipole resonances of dielectric nanospheres in the visible region,” Nano Lett. **12**, 3749–3755 (2012). [CrossRef] [PubMed]

**7. **S. Jahani and Z. Jacob, “All-dielectric metamaterials,” Nat. Nanotechnol. **11**, 23 (2016). [CrossRef] [PubMed]

**8. **A. K. Ospanova, A. Karabchevsky, and A. A. Basharin, “Metamaterial
engineered transparency due to the nullifying of multipole
moments,” Opt. Lett. **43**, 503–506
(2018). [CrossRef]

**9. **A. A. Basharin, M. Kafesaki, E. N. Economou, C. M. Soukoulis, V. A. Fedotov, V. Savinov, and N. I. Zheludev, “Dielectric metamaterials with toroidal dipolar response,” Phys. Rev. X **5**, 11036 (2015).

**10. **A. B. Evlyukhin, C. Reinhardt, and B. N. Chichkov, “Multipole light scattering by nonspherical nanoparticles in the discrete dipole approximation,” Phys. Rev. B **84**, 235429 (2011). [CrossRef]

**11. **Y. Galutin, E. Falek, and A. Karabchevsky, “Invisibility cloaking
scheme by evanescent fields distortion on composite plasmonic
waveguides with Si nano-spacer,” Sci.
Rep. **7**, 12076
(2017). [CrossRef]

**12. **A. Katiyi and A. Karabchevsky, “Si nanostrip pptical
waveguide for on-chip broadband molecular overtone spectroscopy in
near-infrared,” ACS Sens. **3**, 618–623
(2018). [CrossRef]

**13. **J. D. Jackson and R. F. Fox, “Classical Electrodynamics,” Am. J. Phys. **67**, 841 (1999). [CrossRef]

**14. **A. E. Krasnok, A. E. Miroshnichenko, P. A. Belov, and Y. S. Kivshar, “All-dielectric optical nanoantennas,” Opt. Express **20**, 20599–20604 (2012). [CrossRef] [PubMed]

**15. **B. Rolly, B. Stout, and N. Bonod, “Boosting the directivity of optical antennas with magnetic and electric dipolar resonant particles,” Opt. Express **20**, 20376–20386 (2012). [CrossRef] [PubMed]

**16. **D. Markovich, K. Baryshnikova, A. Shalin, A. Samusev, A. Krasnok, P. Belov, and P. Ginzburg, “Enhancement of
artificial magnetism via resonant
bianisotropy,” Sci. Rep. **6**, 22546 (2016). [CrossRef]

**17. **K. V. Baryshnikova, A. Novitsky, A. B. Evlyukhin, and A. S. Shalin, “Magnetic field concentration with coaxial silicon nanocylinders in the optical spectral range,” J. Opt. Soc. Am. B **34**, D36–D41 (2017). [CrossRef]

**18. **V. Kozlov, D. Filonov, A. S. Shalin, B. Z. Steinberg, and P. Ginzburg, “Asymmetric backscattering from the hybrid magneto-electric meta particle,” Appl. Phys. Lett. **109**, 203503 (2016). [CrossRef]

**19. **A. Katiyi and A. Karabchevsky, “Figure of merit of all-dielectric waveguide structures for absorption overtone spectroscopy,” J. Light. Technol. **35**, 2902–2908 (2017). [CrossRef]

**20. **N. Bontempi, K. E. Chong, H. W. Orton, I. Staude, D.-Y. Choi, I. Alessandri, Y. S. Kivshar, and D. N. Neshev, “Highly sensitive biosensors based on all-dielectric nanoresonators,” Nanoscale **9**, 4972–4980 (2017). [CrossRef] [PubMed]

**21. **O. Yavas, M. Svedendahl, P. Dobosz, V. Sanz, and R. Quidant, “On-a-chip biosensing based on all-dielectric nanoresonators,” Nano Lett. **17**, 4421–4426 (2017). [CrossRef] [PubMed]

**22. **P. Spinelli, M. A. Verschuuren, and A. Polman, “Broadband omnidirectional antireflection coating based on subwavelength surface Mie resonators,” Nat. Commun. **3**, 692 (2012). [CrossRef] [PubMed]

**23. **P. D. Terekhov, V. E. Babicheva, K. V. Baryshnikova, A. S. Shalin, A. Karabchevsky, and A. B. Evlyukhin, “Multipole analysis of dielectric metasurfaces composed of non-spherical nanoparticles and lattice invisibility effect,” Phys. Rev. B **99**, 045424 (2019). [CrossRef]

**24. **A. Sayanskiy, M. Danaeifar, P. Kapitanova, and A. E. Miroshnichenko, “All-dielectric
metalattice with enhanced toroidal dipole
response,” Adv. Opt. Mater. **6**, 1800302 (2018). [CrossRef]

**25. **P. D. Terekhov, K. V. Baryshnikova, Y. Greenberg, Y. H. Fu, A. B. Evlyukhin, A. S. Shalin, and A. Karabchevsky, “Enhanced absorption in all-dielectric metasurfaces due to magnetic dipole excitation,” Sci. Rep. **9**, 3438 (2019).

**26. **R. Fleury, F. Monticone, and A. Alù, “Invisibility and cloaking: origins, present, and future perspectives,” Phys. Rev. Appl. **4**, 37001 (2015). [CrossRef]

**27. **Q. Zhao, H. Geng, Y. Wang, Y. Gao, J. Huang, Y. Wang, J. Zhang, and S. Wang, “Hyaluronic acid
oligosaccharide modified redox-responsive mesoporous silica
nanoparticles for targeted drug delivery,”
ACS Appl. Mater. Interfaces **6**,
20290–20299
(2014). [CrossRef]

**28. **G. A. Hughes, “Nanostructure-mediated drug delivery,” in *Nanomedicine in Cancer*, (Pan Stanford, 2017), pp. 47–72.

**29. **E. Blanco, H. Shen, and M. Ferrari, “Principles of
nanoparticle design for overcoming biological barriers to drug
delivery,” Nat. Biotechnol. **33**, 941 (2015). [CrossRef]

**30. **E. V. Salomatina, B. Jiang, J. Novak, and A. N. Yaroslavsky, “Optical properties of
normal and cancerous human skin in the visible and near-infrared
spectral range,” J. Biomed.
Opt. **11**, 64026
(2006). [CrossRef]

**31. **M. Kerker, D.-S. Wang, and C. L. Giles, “Electromagnetic scattering by magnetic spheres,” J. Opt. Soc. Am. **73**, 765–767 (1983). [CrossRef]

**32. **M. Decker, I. Staude, M. Falkner, J. Dominguez, D. N. Neshev, I. Brener, T. Pertsch, and Y. S. Kivshar, “High-efficiency dielectric Huygens’ surfaces,” Adv. Opt. Mater. **3**, 813–820 (2015). [CrossRef]

**33. **W. Liu, A. E. Miroshnichenko, D. N. Neshev, and Y. S. Kivshar, “Broadband unidirectional scattering by magneto-electric core-shell nanoparticles,” ACS Nano **6**, 5489–5497 (2012). [CrossRef] [PubMed]

**34. **J. H. Yan, P. Liu, Z. Y. Lin, H. Wang, H. J. Chen, C. X. Wang, and G. W. Yang, “Magnetically induced forward scattering at visible wavelengths in silicon nanosphere oligomers,” Nat. Commun. **6**, 1–9 (2015). [CrossRef]

**35. **M. I. Abdelrahman, C. Rockstuhl, and I. Fernandez-Corbaton, “Broadband suppression
of backscattering at optical frequencies using low permittivity
dielectric spheres,” Sci.
Rep. **7**, 14762
(2017). [CrossRef]

**36. **M. I. Abdelrahman, H. Saleh, I. Fernandez-Corbaton, B. Gralak, J.-M. Geffrin, and C. Rockstuhl, “Experimental demonstration of spectrally broadband huygens sources using low-index spheres,” APL Photonics **4**, 020802 (2019). [CrossRef]

**37. **R. Dezert, P. Richetti, and A. Baron, “Isotropic huygens dipoles and multipoles with colloidal particles,” Phys. Rev. B **96**, 180201 (2017). [CrossRef]

**38. **E. D. Palik, *Handbook of Optical Constants of
Solids*, vol. 3
(Academic,
1998).

**39. **A. B. Evlyukhin, T. Fischer, C. Reinhardt, and B. N. Chichkov, “Optical theorem and multipole scattering of light by arbitrary shaped nanoparticles,” Phys. Rev. B **94**, 205434 (2016). [CrossRef]

**40. **R. E. Raab and O. L. De Lange, *Multipole Theory in Electromagnetism:
Classical, Quantum, and Symmetry Aspects, with Applications*,
vol. 128 (Oxford
University, 2005).

**41. **R. Alaee, C. Rockstuhl, and I. Fernandez-Corbaton, “An electromagnetic multipole expansion beyond the long-wavelength approximation,” Opt. Commun. **407**, 17–21 (2018). [CrossRef]

**42. **D. W. Pepper and J. C. Heinrich, *The finite element method: basic concepts and applications with MATLAB, MAPLE, and COMSOL* (CRC, 2017). [CrossRef]

**43. ** COMSOL
Multiphysics, www.comsol.com.

**44. **P. D. Terekhov, K. V. Baryshnikova, A. S. Shalin, A. Karabchevsky, and A. B. Evlyukhin, “Resonant forward scattering of light by high-refractive-index dielectric nanoparticles with toroidal dipole contribution,” Opt. Lett. **42**, 835–838 (2017). [CrossRef] [PubMed]

**45. **H. K. Shamkhi, K. V. Baryshnikova, A. Sayanskiy, P. Kapitanova, P. D. Terekhov, A. Karabchevsky, A. B. Evlyukhin, P. Belov, Y. Kivshar, and A. S. Shalin, “Transverse scattering
with the generalised Kerker effect in high-index
nanoparticles,” arxiv 1808.10708
(2018).

**46. **M. J. Berg, C. M. Sorensen, and A. Chakrabarti, “Extinction and the optical theorem Part I Single particles,” J. Opt. Soc. Am. A **25**, 1504 (2008). [CrossRef]